汽蝕UDF編寫中遇到的問(wèn)題:interpretd可以通過(guò),但是計(jì)算幾步之后就報(bào)錯(cuò)了。我把cas和udf都發(fā)上來(lái)了,請(qǐng)大家看看
雖然fluent里邊有自帶的模型,但是為了驗(yàn)證新的汽蝕模型,我需要采用UDF進(jìn)行編寫函數(shù)。其中存在三種方法
第一種是采用下圖1中的cavitation mass rate,但是這種方法,在Fluent 14.5里是需要先打開(kāi)singhal 模型,這種模型依賴于湍流模型,所以不太適合我的UDF。因此采用如下的兩種方法
但為了先檢驗(yàn)我編寫udf的方法對(duì)不對(duì),如下的兩種方法采用的都是成熟的singhal 模型,對(duì)應(yīng)的湍流模型是k-ε模型。另外相關(guān)的物質(zhì)是液氮和氮蒸汽,目前和水的汽蝕是相同的。
第二種方法:
使用下圖中的mass transfer mechanisms,圖2中的mechanism選擇user defined,對(duì)應(yīng)的函數(shù)為
DEFINE_MASS_TRANSFER(water_to_vapor,cell,thread,from_index,from_species_index,to_index,to_species_index)
具體程序見(jiàn)附錄或者文件interaction1.c。
interpreted是沒(méi)有問(wèn)題的,但是計(jì)算十多個(gè)迭代后就報(bào)錯(cuò)了。
第三種方法:
是對(duì)cell-zone-fluid中的兩相(分布為liquid和vapor添加源項(xiàng)),如下圖3(以液相為例)所示
對(duì)應(yīng)的函數(shù)為
DEFINE_SOURCE(liq_src,cell,pri_th,dS,eqn) /*液相*/
DEFINE_SOURCE(vap_src, cell, sec_th, dS, eqn) /*氣相*/
具體的函數(shù)請(qǐng)見(jiàn)附錄或者文件fenkaiFCMabs0.c,interpreted也是沒(méi)有問(wèn)題的,但是同樣計(jì)算十多個(gè)迭代后就報(bào)錯(cuò)了。另外如果采用compiled,則在編譯的時(shí)候就報(bào)錯(cuò)了。為了檢驗(yàn)我的編程方法對(duì)不對(duì),編寫的udf的汽蝕模型是成熟的Singhal模型![汽蝕(空化)UDF編寫中遇到的問(wèn)題]()
圖片1.jpg
![汽蝕(空化)UDF編寫中遇到的問(wèn)題-1]()
圖片2.png
![汽蝕(空化)UDF編寫中遇到的問(wèn)題-2]()
圖片3.jpg
[ Last edited by youy624 on 2013-11-20 at 10:41 ] |