@Afr1yne 谢谢您
Y
yjc95
@yjc95
帖子
-
buoyantPimpleFoam相关问题 -
buoyantPimpleFoam相关问题@Afr1yne 大佬，追问一个问题，问题3中的湍流Pr该咋取呀？
-
buoyantPimpleFoam相关问题@李东岳 东子老师说的对，我基础太差了，我该好好去读您网站上给出的书籍。
-
buoyantPimpleFoam相关问题@yjc95 东子老师，您有空看看哈
-
buoyantPimpleFoam相关问题@李东岳 突然看到网站已经被装饰上了彩灯和雪人，才发现马上要过圣诞节了，很有喜庆的氛围，提前祝东子老师圣诞节快乐啊！
-
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; -
buoyantPimpleFoam相关问题@Banbor 谢谢您，我才意识到原来状态方程用boussinesq近似有专门的求解器，我还以为都在buoyantPimpleFoam里面
-
buoyantPimpleFoam计算问题@Afr1yne 东子老师的虚拟机里没有v2112呀，是自己装的嘛
-
buoyantPimpleFoam相关问题@李东岳 老师，我切换到of7，输入buoyantBoussinesqPimpleFoam，提示已经被buoyantPimpleFoam取代了，然后我在of7下输入buoyantPimpleFoam，算出来的结果和上面的图一样。您在 https://cfd-china.com/topic/7624 这里面提到您是把bouyantBoussinesqPimpleFoam移植到openfoam10，没有别的招了，我想试试您这个方法，能教教我嘛？
-
buoyantPimpleFoam计算问题@Afr1yne 这俩我都试过了，of8下没有buoyantFoam，buoyantPimpleFoam在改变膨胀系数beta值的时候会出现非物理现象（https://www.cfd-china.com/topic/7611/buoyantpimplefoam相关问题/63
这个帖子的最后几个回复有图），东子老师说of8的buoyantPimpleFoam有bug，建议我试试buoyantBoussinesqPimpleFoam，我看您在of8下可以用buoyantBoussinesqPimpleFoam，但我试了之后发现用不了。 -
buoyantPimpleFoam计算问题@Afr1yne 大佬，您是用of哪个版本下的bouyantBoussinesqPimpleFoam呀？我用of8提示没有这个求解器
-
buoyantPimpleFoam相关问题@李东岳 东子老师，我在案例文件夹下切换成of8，输入bouyantBoussinesqPimpleFoam
提示找不到命令 “bouyantBoussinesqPimpleFoam”，打开方式不对嘛？
-
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值，就出现了很大的区别，我百思不得其解，想请教下您。 -
buoyantPimpleFoam相关问题@李东岳 您的意思是，我得先输入of8，切换到版本8，然后blockmesh， 最后再buoyantPimpleFoam嘛？
-
buoyantPimpleFoam相关问题@李东岳 东子老师，我把可压+瞬态算法全部看完了，代码也看完了，打算用它来模拟储热过程。我在of8里面找到了buoyantPimpleFoa下面的hotRoom案例，感觉和我想要的储热过程很相似，所以我就把floor设置成了低温水入口，ceiling设置成了高温水出口，具体的设置参数是这样的：
0文件夹下有5个文件，其中：
nutdimensions [0 2 -1 0 0 0 0]; internalField uniform 1e-6; boundaryField { floor { type nutkWallFunction; value uniform 0; } ceiling { type nutkWallFunction; value uniform 0; } fixedWalls { type nutkWallFunction; value uniform 0; } }p
dimensions [1 -1 -2 0 0 0 0]; internalField uniform 1e5; boundaryField { floor { type zeroGradient; } ceiling { type fixedValue; value uniform 0; } fixedWalls { type zeroGradient; } }p_rgh
dimensions [1 -1 -2 0 0 0 0]; internalField uniform 0; boundaryField { floor { type zeroGradient; } ceiling { type zeroGradient; } fixedWalls { type zeroGradient; } }T
dimensions [0 0 0 1 0 0 0]; internalField uniform 363; boundaryField { floor { type fixedValue; value nonuniform 283; } ceiling { type zeroGradient; } fixedWalls { type zeroGradient; } }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 zeroGradient; } fixedWalls { type noSlip; } }constant文件夹下面修改了momentum Transport，pRef和thermophysicalProperties：
momentum TransportsimulationType laminar;pRef
dimensions [1 -1 -2 0 0 0 0]; value 1e5;thermophysicalProperties
thermoType { type heRhoThermo; mixture pureMixture; transport const; thermo hConst; equationOfState Boussinesq; specie specie; energy sensibleEnthalpy; } pRef 100000; mixture { specie { molWeight 18; } equationOfState { rho0 1000; T0 273; beta 3e-05; } thermodynamics { Cp 4182; Hf 0; } transport { mu 1e-03; Pr 7; } }system文件夹下就修改了fvSolution文件，增加了pRefCell，pRefPoint，pRefValue：
PIMPLE { momentumPredictor yes; nOuterCorrectors 1; nCorrectors 2; nNonOrthogonalCorrectors 0; pRefCell 0; pRefPoint 0; pRefValue 0; }做完上述修改之后，输入buoyantPimpleFoa，报错了：
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // Create time Create mesh for time = 0 --> FOAM FATAL IO ERROR: keyword version is undefined in dictionary "/home/dyfluid/OpenFOAM/dyfluid-11/run/hotRoom/constant/polyMesh/points" file: /home/dyfluid/OpenFOAM/dyfluid-11/run/hotRoom/constant/polyMesh/points from line 10 to line 13. From function const Foam::entry& Foam::dictionary::lookupEntry(const Foam::word&, bool, bool) const in file db/dictionary/dictionary.C at line 799. FOAM exiting我觉得可能是边界条件没设置好，检查再三感觉也没啥问题，麻烦您抽空看看，谢谢您！
-
虚拟机相关问题@李东岳 东子老师，重装就可以了
-
虚拟机相关问题@李东岳 东子老师，您网站上的虚拟机是不是没办法使用共享文件夹这个功能呀？我设置了共享文件夹，但是在hgfs文件夹里面没有看到windows共享文件夹里面的文件。
-
buoyantPimpleFoam相关问题@李东岳 老师，这俩没加上体积
-
buoyantPimpleFoam相关问题@李东岳 东子老师，根据公式(36)，红色箭头部分应该缺了$1/V_{p}$吧。
-
buoyantPimpleFoam相关问题@李东岳 那
fvc::div(phi)和of8中的fvc::div(phiHbyA)是代表一个意思，只不过是of版本不同吗？