buoyantPimpleFoam相关问题
-
@李东岳 东子老师，我切换到of8果然可以正常计算了，但算完之后发现一个很奇怪的问题。模型简介：计算域是长方体，floor面是底部进水，ceiling面是顶部出水，fixedwalls是四周的壁面，层流，状态方程我用的是Boussinesq假设，thermophysicalProperties文件的设置如下：
thermoType { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState Boussinesq; specie specie; energy sensibleEnthalpy; } pRef 100000; mixture { specie { molWeight 18; } equationOfState { rho0 988; T0 333; beta 1e-05; } thermodynamics { Cp 4182; Hf 0; } transport { mu 1e-03; Pr 7; } }其中beta值我原先设置的是1e-06，温度云图看起来也比较合理，只是中间的斜温层比较薄，顶部水流出，温度交界面也是水平的。接着，我设置成1e-05，顶部温度交界面出现了非水平的现象，这是不合理的。水的膨胀系数是2.6e-04，我也算了一遍，发现也出现了不合理的现象。
这三个case我只改变了beta值，就出现了很大的区别，我百思不得其解，想请教下您。 -
@李东岳 东子老师，我试了下v2206的buoyantBoussinesqPimpleFoam，发现没有of8出现的非物理现象，我专门去您的网站看了buoyantBoussinesqPimpleFoam的算法，很有收获，还有几个疑问想请教您。
问题1
alphat文件应该是能量方程中的$\alpha_{eff}$吧？这个$\alpha_{eff}$应该就是热扩散系数吧？那为啥internalField是0呢？代码中type，Prt和value是代表啥意思？
internalField uniform 0; boundaryField { floor { type alphatJayatillekeWallFunction; Prt 0.85; value uniform 0; }
问题2
这是p_rgh文件里面的，这里出现一个rhok，我知道它代表的是动量方程的$\rho_{k}$，为啥这里给的值是0？
internalField uniform 0; floor { type fixedFluxPressure; rho rhok; value uniform 0; }
问题3
transportProperties文件里面有5个物性参数，最后一个湍流Pr数我没找到是啥意思。因为我选的是层流模型，湍流Pr应该不起作用吗？
transportModel Newtonian; // Laminar viscosity nu 1e-06; // Thermal expansion coefficient beta 2.6e-04; // Reference temperature TRef 300; // Laminar Prandtl number Pr 7; // Turbulent Prandtl number Prt 0.85;
