【仿真技巧】fluent TUI使用技巧及燃烧室自动化仿真案例

 #include "udf.h"  DEFINE_PROFILE(pressure_profile,t,i) {    real x[ND_ND];    /* this will hold the position vector */    real y;    face_t f;    begin_f_loop(f,t)      {       F_CENTROID(x,f,t);       y = x[1];       F_PROFILE(f,t,i) = 1.1e5 - y*y/(.0745*.0745)*0.1e5;      }    end_f_loop(f,t) }

(define my-iterate    (lambda (mark-temp)        (if (< (pick-a-real "/re/vi/max solid () temperature no" 1) mark-temp)            (begin                  (ti-menu-load-string "/so/dti 1 20")                 (my-iterate mark-temp)))))

Vmax*(5.0*exp(t-0.3 [s]/2.8 [s]))

/file/read-case-data test.case /solve/set/under-relaxation/mom 0.2 /solve/set/under-relaxation/pressure 0.4 /solve/iterate 200  /file/write-case-dat result.cas

;; This journal is writen by wenming jing of Ansys China, plz contact wenming.jing@ansys.com for help;; read case/mesh/file/read-case combustor.msh.h5;; confirm overwrite/file/confirm-overwrite? yes;;choose steady define/models/steady?;;activate energy equationdefine/models/energy? yes yes yes ;;choose rke turbulence model define/models/viscous/ke-realizable? yes
;;define operating conditionsdefine/operating-conditions/operating-pressure 101325 gravity yes 0 -9.8 0 q
;; combustion model settingdefine/models/species/partially-premixed-combustion? yes yes no yes 101325 no no 0.1 315 500 yes "jet-a" 1 0 "n2" 0 0.767 "o2" 0 0.233 , , , yes 15 200 0.25 0.25 20 298 no "combustor.pdf" yes
;;DPM model setting define/models/dpm/interaction/coupled-calculations? yesdefine/models/dpm/interaction/dpm-iteration-interval 30define/models/dpm/interaction/update-dpm-sources-every-flow-iteration? nodefine/models/dpm/unsteady-tracking no/define/models/dpm/numerics/tracking-parameters 50000 no 5

;;DPM injection settingdefine/models/dpm/injections/create-injection , yes droplet yes surface diesel-liquid fuelinlet () no no no no no no jet-a no 0 0 10 0.0001 315 0.01 ;;read injection file;;/file/read-injections injection-file
;; inletair1 settingdefine/boundary-conditions/velocity-inlet inletair1 , , , , , 40 , , , 500 , , , , , , , , , , , , ;;another way inletair1 setting ;;define/boundary-conditions/set/velocity-inlet inletair1 () vmag no 31 temperature no 319 fmean no 0 premixc no 0 fvar no 0 q ;; inletair2 settingdefine/boundary-conditions/velocity-inlet inletair2 , , , , , 30 , , , 500 , , , , , , , , , , , , ;; fuelinlet settingdefine/boundary-conditions/velocity-inlet fuelinlet , , , , , 15 , , , 300 , , , , , , , , 1 , , , ;; outlet settingdefine/boundary-conditions/pressure-outlet outlet , , 0 , 800 , , , , , 5 10 , 0 , , , , , , , , , 
;; wall settingdefine/boundary-conditions/wall mainbody:1 , , , , , , 0 , , , , , , , , , , , , , , , , 
;;hybrid initializationsolve/initialize/hyb-initialization
;;report definitionsolve/report-definitions/add max-temp volume-max zone-names mainbody , field temperature report-type volume-max q
;;add report to plot and filesolve/report-plot/add max-temp plot-instantaneous-values? yes report-defs max-temp () q
solve/report-files/add max-temp report-defs max-temp () write-instantaneous-values? yes file-name max-temp.out q
;;define a surface planedisplay/surface/plane-surface midplane zx-plane 0
;;setting up midplane velocity contourdisplay/objects/create contour midplane-velocity field velocity-magnitude surfaces-list midplane , q
;;setting up midplane velocity vector display/objects/create vector midplane-v-vector field velocity-magnitude surfaces-list midplane , q
;;setting up auto save /file/auto-save/ case-frequency if-case-is-modified data-frequency 1000 retain-most-recent-files yes root-name combustor
;;freeze specify equation/solve/set/equations/premixed no
;;activate permixed equation/solve/set/equations/premixed no
;; solve for 200 iterationsolve/iterate 100
;;activate specify equation/solve/set/equations/premixed yes
;;patch temperaturesolve/patch mainbody () temperature 3000;;patch progress variablesolve/patch mainbody () premixc 1 ;; solve for 1000 iterationsolve/iterate 1000