;;
;; Initial setup for Fluent
;;
;; run with:
;;   fluent 2d -i <this_file>
;; or rather use another scm file to define some variables first,
;; then (load "initial-setup.scm").

;; Solution model: steady/unsteady-1st-order/unsteady-2nd-order (/unsteady-global-time)
(if (not (symbol-bound? 'solutionmodel (the-environment)))
	(define solutionmodel "unsteady-1st-order") )
(if (not (symbol-bound? 'axisymmetric (the-environment)))
	(define axisymmetric #f))             ;; Use axisymmetric solver? (#t/#f)
(if (not (symbol-bound? 'timestep (the-environment)))
	(define timestep 0.05))               ;; time step for iterations

(if (not (symbol-bound? 'meshfile (the-environment)))
	(define meshfile "test.msh"))         ;; Input mesh

;; Load file, do some checks. Axisymmetric must be set first.
(ti-menu-load-string (format #f "/define/models/axisymmetric ~a" (cond ((eq? axisymmetric #t) "y") (else "n"))))
(ti-menu-load-string (format #f "/file/read-case ~a" meshfile))
(if (ti-menu-load-string "/grid/check") () (begin
	(display "Error: grid check failed") (exit) ))

;;;; Setup simulation parameters
(ti-menu-load-string (format #f "/define/models/~a y" solutionmodel))
;(ti-menu-load-string "/define/models/noniterative-time-advance y")
(ti-menu-load-string "/define/operating-conditions/gravity y 0 -9.81")
(ti-menu-load-string "/define/operating-conditions/operating-density y")
(ti-menu-load-string "/define/materials/copy fluid water-liquid")
;; multiphase 
;;(ti-menu-load-string "/define/models/multiphase vof 2 geo-reconstruct 0.25 n n n")
(ti-menu-load-string "/define/models/multiphase vof 2,,,,,")
(ti-menu-load-string "/define/phases/phase-domain phase-1 phase-1 y air")
(ti-menu-load-string "/define/phases/phase-domain phase-2 phase-2 y water-liquid")
(ti-menu-load-string "/define/phases/interaction-domain n y constant 0.07275 ,") ;; surface tension water-air at 293K
;(ti-menu-load-string "/define/phases/interaction-domain n y none ,")

(ti-menu-load-string "/solve/set/flow-warnings n")
;; monitors, plot windows
(ti-menu-load-string "/solve/monitors/residual/print n")
(ti-menu-load-string "/solve/monitors/residual/plot n")
(ti-menu-load-string "/solve/monitors/residual/window 5")

(ti-menu-load-string "/display/set/rendering-options/double-buffering y")

(ti-menu-load-string "/display/set/contours/clip-to-range y")
(ti-menu-load-string "/display/set/contours/filled-contours y")
(ti-menu-load-string "/display/set/contours/global-range y")
(ti-menu-load-string "/display/set/contours/line-contours n")
(ti-menu-load-string "/display/set/contours/log-scale n")
(ti-menu-load-string "/display/set/contours/n-contour 2")
(ti-menu-load-string "/display/set/contours/node-values y")
(ti-menu-load-string "/display/set/contours/render-grid y")
(ti-menu-load-string "/display/set/contours/surfaces ")
(ti-menu-load-string "/display/set/line-weight 1")


(ti-menu-load-string "/display/open 0")

;; Initialize solution
(init-customfield)
(ti-menu-load-string "/solve/initialize/compute-defaults/all-zones")
(ti-menu-load-string "/solve/initialize/initialize-flow y")
(patch-phase)
(rpsetvar 'physical-time-step timestep)

;; XXX first view
;; Somehow GUI has to be used before it's ok (empty grid points!)
(cx-gui-do cx-activate-item "MenuBar*DisplayMenu*Contours...")
(cx-gui-do cx-activate-item "Contours*PanelButtons*PushButton1(OK)")
(cx-gui-do cx-activate-item "Contours*PanelButtons*PushButton1(Cancel)")
(ti-menu-load-string "/display/set/grid-surfaces (0 1 2 3 4 5 6)") ;; XXX static surface number !!!
(ti-menu-load-string "/display/set/render-grid y")

;; vim:ts=2:sw=2:ft=scheme: