FLUENT Scheme 腳本調(diào)用及其案例

2016-08-19  by:CAE仿真在線  來源:互聯(lián)網(wǎng)

FLUENT 中使用 Scheme,如果好好利用,能夠大大提高工作效率。這是自己編寫的第一個(gè) journal 文件,用來設(shè)置一些選項(xiàng)。由于經(jīng)常使用 User-defined Materials Database,在 journal 中添加曾經(jīng)遭遇錯(cuò)誤,后來才發(fā)現(xiàn)是 Material 的名稱出了問題。

現(xiàn)在的計(jì)算平臺(tái)有很多是 Linux Cluster,常用的作業(yè)管理系統(tǒng)采用 PBS。如果要采用 FLUENT 進(jìn)行計(jì)算就需要解決在 Text Mode 下調(diào)用 FLUENT 進(jìn)行求解的問題,而且多數(shù)時(shí)候

需要執(zhí)行并行版本的 FLUENT。

Linux 調(diào)用 FLUENT 的文本界面其實(shí)很簡單,比如

出現(xiàn)了 FLUENT 的 console,這樣就能輸入命令并執(zhí)行了。

不想手工輸入命令的話,寫一個(gè) journal 或 scheme 腳本文件,讓 FLUENT 啟動(dòng)后自動(dòng)執(zhí)行。

最好加上 I/O 重定向。

 

fluent 3d -g -i script.scm > out.log

遠(yuǎn)程連接時(shí)后臺(tái)執(zhí)行就比較重要了。

 

fluent 3d -g -i script.scm > out.log &

以上是之前在 PC 上常用的模式,還會(huì)利用 screen 之類的進(jìn)行監(jiān)視。

那么并行版本又如何呢?

本地并行的話,基本上就是

利用多個(gè)節(jié)點(diǎn)的話,需要做些處理

多了幾個(gè)選項(xiàng),其中

-pethernet 表示各節(jié)點(diǎn)利用局域網(wǎng)連接進(jìn)行通信,發(fā)現(xiàn)用 default 選項(xiàng)的話會(huì)報(bào)錯(cuò)的,至少 fluent12 是這個(gè)樣子。

-cnf=myhosts 表示運(yùn)行的節(jié)點(diǎn)信息

-ssh 表示各節(jié)點(diǎn)采用 ssh 進(jìn)行登錄,默認(rèn)的話采用 rsh,如果沒有啟用 rsh 服務(wù)的話就崩潰了。

到這里基本可以了。

如果要在 PBS 中進(jìn)行調(diào)用,啟動(dòng)命令可以不加 i/O 重定向和后臺(tái)運(yùn)行的符號(hào)

 

1、添加菜單項(xiàng) (Text Mode)

(define ti-flab-iterate
(lambda ()
(display "Hello World!")
)
)
(define flowlab-menu
(make-menu
"flowlab"
("iterate" #t ti-flab-iterate "Iterate.")
)
)
(define add-menu
(lambda (name menu test value help)
(ti-menu-insert-item!
menu
(make-menu-item name test value help)
)
#t
)
)
(add-menu "flowlab/" main-menu #t flowlab-menu "")

---

---

2 輸出當(dāng)前的迭代次數(shù)和參數(shù)

(define (export-niters fname parameter)

(let

(

(niter (%iterate 0))

(out-port (open-file fname "w"))

)

(if (not out-port)

(cx-error-dialog "unable to open output file for exporting niters")

)

(format out-port "$~a = ~a~%" parameter niter)

(flush-output-port out-port)

)

)

---

3 字符串處理

;takes any string as argument 'str', breaks it on the basis of blank-space or
;newline character into multiple strings. These multiple strings are returned
;as a list. 
(define (tokenizer str)
(let
(
(l (string->list str))
(result '())
(delimitters (list #/space #/newline #/) #/())
(started? #f)
(temp '())
)
(for-each
(lambda (c)
(if started?
(if (not (memv c delimitters))
(set! temp (append temp (list c)))
(begin
(set! result (append result (list (list->string temp))))
(set! temp '())
(set! started? #f)
)
)
(if (not (memv c delimitters))
(begin
(set! temp (append temp (list c)))
(set! started? #t)
)
)
)
)
l
)
(if started?
(set! result
(append result (list (list->string temp)))
)
)
result
)
)

---

4 自動(dòng)保存,記錄輸出信息的迭代腳本

;;
;; iterate.scm
;;
(define casename "coal_inj_00")
(define iternum 500)
;
(define hy-start-transcript
(lambda (filename)
(if (file-exists? filename)
(remove-file filename)
#t
)
(ti-start-transcript filename #t)
#t
)
)
;;
(define (hy-stop-transcript)
(if transcript-open?
(stop-transcript #t)
#t
)
#t
)
;;
(define (log-exec-proc exec-proc trn-file)
(begin
(hy-start-transcript trn-file)
(exec-proc)
(newline)
(hy-stop-transcript)
)
)
;;
(define my-iterate
(lambda (num)
(if (number? num)
(begin
(read-case-data casename)
(iterate num)
(write-case-data (format #f "~a-i~a" casename num))
#t
)
#f
)
)
)
;;
(log-exec-proc (lambda () (my-iterate iternum)) (string-append casename ".trn"))

5 讀取 report 文件中的數(shù)值

;;

;;

;;
(define (hy-read-rp-file filename)
(let ((p (open-input-file (format #f "~a" filename))))
(do
((x (read p) (read p)))
(
(or (number? x) (eof-object? x))
(close-input-port p)
(if (number? x) x #f)
)
)
)
)

---

6 讀取數(shù)據(jù)到 list 中

(define hy-read-data-to-list
(lambda (fname)
(let ((p (open-input-file (format #f "~a" fname))))
(let f ((x (read p)) (l '()))
(if (or (eof? p) (eof-object? x))
(begin
(close-input-port p)
l
)
(begin
(display x)
(newline)
(if (number? x)
(set! l (append l (list x)))
)
(f (read p) l)
)
)
)
)
)
)

---

7 導(dǎo)出 bc 信息

(define (hy-export-bc-names)
(for-each
(lambda (name)
(display
(format #f "{~a, /"~a/", /"~a/"},/n"
(zone-name->id name)
name
(zone-type (get-zone name))
)
)
)
(inquire-zone-names)
)
)

---

8 顯示 list 內(nèi)容

;;;
(define hy-display-list
(lambda (l)
(do ((i 0 (+ i 1)))
((>= i (length l)))
(display (format #f "~a --- ~a/n" i (list-ref l i)))
)
)
)

---

9 非穩(wěn)態(tài)計(jì)算時(shí)自動(dòng)保存腳本

;; Scheme file to autosave files at a particular time interval,

 

;;original: FLUENT Online Technical Support "Solution 656 :Autosave at specified time intervals "

;; Usage:

;; To use the scheme file to save files at a given time interval,

;; 0. Read case and data files, then read autosave.scm through File > Read > Scheme

;; 1. Open the Execute Commands panel in the GUI (Solve>Execute Commands)

;; 2. Increase the number of Defined Commands by

;; 3. Click ON next to the most recent command and select Time Step

;; next to the drop down arrow under When

;; 4. Enter the following text in the field under Command:

;; (autowrite "wcd" "autoname-" 0.2)

;; 5. Enter the compression status below. Change only the text between the quotation marks.

;; For compressed files use ".gz"

;; For uncompressed files use ""

(define g_zip "")

;;

;; The value of 0.2 should be replaced with the desired interval for the current case

;; The value of "wcd" should be replaced with

;; "wd " for data files only

;; "wcd " for case and data files

;; The value of "autoname-" should be replaced with the actual pathname+filename :

;; if it is written in current directory,pathname can be omitted.

;;default setting

(define g_filename "auto-")

(define g_wcom "wd")

(define (writefiles current-time)

(let* ((time-string) (writecommand))

(rpsetvar 'autowrite/last-file-time current-time)

(set! writecommand (string-append g_wcom " " g_filename g_zip))

(ti-menu-load-string writecommand)

)

)

(define (autowrite wcom filename interval)

(let* ( (current-time) (overlap) (last-write) (file-write-time) )

(if (not (rp-var-object 'autowrite/last-file-time))

(rp-var-define 'autowrite/last-file-time (rpgetvar 'flow-time) 'real #f)

)

(set! current-time (rpgetvar 'flow-time))

(set! last-write (rpgetvar 'autowrite/last-file-time))

(set! overlap (- last-write (* interval (truncate (/ last-write interval)))))

(set! file-write-time (+ interval (- last-write overlap)))

(if (> current-time file-write-time)

(begin

(cond

((truncate (/ last-write interval)<? 10)

(set! g_filename (string-append filename "000" (number->string (inexact->exact (truncate(/ last-write interval))))))

)

((truncate (/ last-write interval)<? 100)

(set! g_filename (string-append filename "00" (number->string (inexact->exact (truncate(/ last-write interval))))))

)

((truncate (/ last-write interval)<? 1000)

(set! g_filename (string-append filename "0" (number->string (inexact->exact (truncate(/ last-write interval))))))

)

(else

(set! g_filename (string-append filename (number->string (inexact->exact (truncate(/ last-write interval))))))

)

)

(set! g_wcom wcom)

(writefiles current-time)

)

)

)

)

相關(guān)標(biāo)簽搜索：FLUENT Scheme 腳本調(diào)用及其案例 FLUENT Scheme 腳本 fluent Fluent培訓(xùn) Fluent流體培訓(xùn) Fluent軟件培訓(xùn) fluent技術(shù)教程 fluent在線視頻教程 fluent資料下載 fluent分析理論 fluent化學(xué)反應(yīng) fluent軟件下載 UDF編程代做