running Fluent on Apocrita Cluster

- 日理万妓 2021-12-24 15:51 171阅读 0赞

two files: code.sh, code.jou

code.sh

#!/bin/bash
    #$ -cwd
    #$ -j y
    #$ -m bea
    #$ -M k.ai@qmul.ac.uk
    #$ -pe smp 4
    #$ -l h_rt=24:0:0
    #$ -l h_vmem=4G
    
    module load ansys
    fluent 3ddp $FLUENT_OPTS -g -t4 -rsh -i code.jou

code.jou

rc e387_wave.cas.gz
    rd e387_2550.dat.gz
    solve/dual-time-iterate 600 20
    solve/set/time-step 0.01
    wd e387_3000.dat.gz
    exit
    yes

More advanced,

rc st_e387_w_hy-14M_open_invisid.cas.gz
    rd st_e387_w_hy-14M_open_invisid.dat.gz
    ;set drag moment and monitor in y coordinate on "blades" surface
    solve/monitors/force/unscaled? yes
    solve/monitors/force/set-drag-monitor cd yes blades () no yes thrust-history no no 0 1 0
    solve/monitors/force/set-moment-monitor moment yes blades () no yes moment-history no no 0 0 0 0 1 0    
    ;Coupled solver, Coupled with Volume Fractions [no] 
    ;/solve/set p-v-coupling 24  no 
    ;/solve/set/p-v-controls [courant-number] [explicit momentum] [explicit pressure] 
    ;/solve/set/p-v-controls  1 0.5 0.5 
    ;/solve/set/under-relaxation/k 0.5           
    ;/solve/set/under-relaxation/epsilon  0.5   
    ;/solve/set/under-relaxation/turb-viscosity 0.5
    (display "Save the residual in a file") (newline)
        (let ((writefile (lambda (p)
        (define np (length (residual-history "iteration")))
        (let loop ((i 0))
        (if (not (= i np))
        (begin (define j (+ i 1))
        (display (list-ref (residual-history "iteration") (- np j)) p) (display " " p)
        (display (list-ref (residual-history "continuity") (- np j)) p) (display " " p)
        (display (list-ref (residual-history "x-velocity") (- np j)) p) (display " " p)
        (display (list-ref (residual-history "y-velocity") (- np j)) p) (display " " p)
        (display (list-ref (residual-history "z-velocity") (- np j)) p) (display " " p)
        (display (list-ref (residual-history "k") (- np j)) p) (display " " p)
        (display (list-ref (residual-history "omega") (- np j)) p)
        (newline p)
        (loop (+ i 1))
        )
        )
        )
        ) )
        (output-port (open-output-file "residual_1000_e387_udf.dat")))
        (writefile output-port)
        (close-output-port output-port))
    solve/iterate 3000
    wd e387_st_hy_14M_open_channel_3k.dat.gz
    exit
    yes

转载于:https://www.cnblogs.com/code-saturne/p/11019408.html