buoyantPimpleFoam相关问题
-
@李东岳 东子老师，我试了下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; -
@李东岳 东子老师，新年好呀。我带着新问题来向您请教了。我用v2206版本下的heatTransfer-buoyantBoussinesqPimpleFoam-hotRoom作为模板来模拟储热池的放热过程，我把floor设置成低温水进口，把ceiling设置成高温水出口，采用层流，介质为水，算出来的结果温度场很奇怪，中间的斜温层很厚。
我查阅文献，并且用fluent算出来的斜温层应该非常薄才对。
此外，我把所有物性参数扩大和缩小1000倍，温度场没变化。
具体设置如下：- alphat
FoamFile { version 2.0; format ascii; class volScalarField; object alphat; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { floor { type alphatJayatillekeWallFunction; Prt 0.1; value uniform 0; } ceiling { type alphatJayatillekeWallFunction; Prt 0.1; value uniform 0; } fixedWalls { type alphatJayatillekeWallFunction; Prt 0.1; value uniform 0; } }- nut
FoamFile { version 2.0; format ascii; class volScalarField; object nut; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -1 0 0 0 0]; internalField uniform 0; boundaryField { floor { type nutkWallFunction; value uniform 0; } ceiling { type nutkWallFunction; value uniform 0; } fixedWalls { type nutkWallFunction; value uniform 0; } }- p
FoamFile { version 2.0; format ascii; class volScalarField; object p; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 1e5; boundaryField { floor { type zeroGradient; } ceiling { type zeroGradient; } fixedWalls { type zeroGradient; } }- p_rgh
FoamFile { version 2.0; format ascii; class volScalarField; object p_rgh; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 2 -2 0 0 0 0]; internalField uniform 1e5; boundaryField { floor { type zeroGradient; rho rhok; } ceiling { type zeroGradient; rho rhok; } fixedWalls { type zeroGradient; rho rhok; } }- T
FoamFile { version 2.0; format ascii; class volScalarField; object T; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 0 0 1 0 0 0]; internalField uniform 363; boundaryField { floor { type fixedValue; value uniform 283; } ceiling { type zeroGradient; } fixedWalls { type zeroGradient; } }- U
FoamFile { version 2.0; format ascii; class volVectorField; object U; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dimensions [0 1 -1 0 0 0 0]; internalField uniform (0 0 0); boundaryField { floor { type fixedValue; value uniform (0 0.1 0); } ceiling { type fixedValue; value uniform (0 0.1 0); } fixedWalls { type noSlip; } }- transportProperties
FoamFile { version 2.0; format ascii; class dictionary; object transportProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // transportModel Newtonian; // Laminar viscosity nu 1e-06; // Thermal expansion coefficient beta 3e-04; // Reference temperature TRef 273; // Laminar Prandtl number Pr 7; // Turbulent Prandtl number Prt 0.1;- turbulenceProperties
FoamFile { version 2.0; format ascii; class dictionary; object turbulenceProperties; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // simulationType laminar; /* RAS { RASModel kEpsilon; turbulence on; printCoeffs on; } */麻烦您抽空看看，谢谢您！
"其中的
指的是相邻网格的速度，在这里是用当前时刻已知的速度代入么？同样的，方程(41)中的
应该是要求的控制体体心处的压力吧，
是否是用当前时刻各个面的压力代入，这个压力好像之前并不知道吧，是咋求出
