OpenFOAM11中使用surfaceFilme 模型进行燃烧仿真,温度失控导致浮点溢出
-
各位老师,师兄,师姐,大佬们好,
在OpenFOAM11中,我使用surfaceFilm模型进行一个3D仿真。仿真流程如下:- 先运行一个关闭surfaceFilm和spray injection 模型,让T_inlet=900k 运行0.05s 初始化流场。
- 激活燃烧,喷雾和薄膜模型,向模型中注入正庚烷C7H16(LPT),然后进行燃烧仿真。
在只使用燃烧和喷雾模型的时候,仿真可以正确的进行并获得有物理意义的结果。
在激活surfaceFilm模型后,模型设置参照openFOAM自带算例的“heatBox”,仿真疑似因为在薄膜位置计算温度时浮点溢出而不能运行。
在"CFDonline"中,我查到在之前的版本中 可以使用“fvOption” 来对温度的上下限进行限制以防温度失控,但是在OpenFOAM11版本中似乎没有该功能。我在中外各个网站查了一圈,有关于该模型的使用案例似乎特别少。
请问这个问题出现的原因是什么?我应该如何解决呢
下附报错图片,我认为比较有用的设置文件截图,以及我的仿真文件
报错如下
可能有用的文件截图:
film的温度
film的fvModel
关闭薄膜模型时的结果
求求大佬们救救孩子吧~~~~~~~
附上我的仿真文件下载链接:
https://www.jianguoyun.com/p/DTKpvfgQlPqQDBi01tUFIAA -
自回一下:
问题已解决,经过多次尝试需要在film的physical properteis 中增加温度限制,粒子数据之类的信息,还需要查询具体的粘度系数和热熔系数公式,我把文件内容放在下面,希望对后来人有用thermoType { type heRhoThermo; mixture pureMixture; transport polynomial; thermo janaf; equationOfState perfectGas; specie specie; energy sensibleInternalEnergy; } mixture { specie { molWeight 100.20557; } thermodynamics { Tlow 200; Thigh 3500; Tcommon 1000; highCpCoeffs ( 20.4565203 0.0348575357 -1.09226846e-05 1.67201776e-09 -9.8102485e-14 -32555.6365 -80.4405017 ); lowCpCoeffs ( 11.1532994 -0.00949419773 0.000195572075 -2.49753662e-07 9.84877715e-11 -26768.8904 -15.9096837 ); } transport { muCoeffs<8> ( -10.378 0.244 -0.000054 0 0 0 0 0 ); //( 1000 -0.05 0.003 0 0 0 0 0 ); kappaCoeffs<8> ( 80.10 3.45542 0.00128822 -0.00000166835 0 0 0 0 ); } equationOfState { rho 684.0; // 庚烷的密度,单位 kg/m³ } specie { Hf 0; } } sigma { type constant; sigma 0.07; } deltaWet 1e-4;