CFD中文网

@evensun 能在高校待时间长的，必然都有外快。各有各的路子，有些大家能看见，有些大家看不见。否则待个5 6年必然就走了。

或许可以将代码改为用不可压缩求解器求解的边界上的温度和导热率来获得，边界上的传热量。例如：

if(db().foundObject<volScalarField>("T")) { const volScalarField& T = db().lookupObject<volScalarField>("T"); const fvPatchScalarField& Tp = T.boundaryField()[patchi]; qDot = kappa*Tp.snGrad(); }

其中kappa需要根据你的物性和边界上的温度场来计算获得

@sibo 大佬，请教下，如何在这个buoyantPimpleFoam中添加两种流体，进而形成双扩散自然对流的情形哈？

@东岳 是啊我也看到了，改成压力基说是在low mach时更好。感觉各种拼凑。。

grad(U)和wallGrad(U)有关系么？？

从你的描述来看：

wallGradU.boundaryField()[patchi] =-U.boundaryField()[patchi].snGrad();

wallGrad(U)是面法向梯度，grad(U)是梯度。这个区别比较大。如果新版本没有了wallGrad(U)，会不会grad(U)在处理壁面的时候已经当成面法向梯度了，你研究研究这个角度对不对？

受OpenFOAM中wall function实现方式的影响，贴近壁面处的du/dy并不是真实的速度梯度，因此用du/dy求出来的并不是真实的wall shear stress。

也就是说du/dy的算法只适用于laminar flow，而不适用于turbulent flow。

重新运行了算例, 提高了内存请求(200GB->251GB)。同样的结果。请求内存似乎不是问题所在，峰值 RSS 和 200G内存一样。

JobID Elapsed NCPUS AveRSS MaxRSS AveVMSize ReqMem ------------ ---------- ---------- ---------- ---------- ---------- ---------- 14389 06:00:15 192 251Gn 14389.batch 06:00:15 48 302K 194452190K 113468K 251Gn 14389.extern 06:00:15 192 115K 144K 107948K 251Gn

我在做gas-liquid compressible flow时，需要water的equation of state 去解液态场的压力。
我知道现在有这样几种方法：
the stiffened-gas equation
the Tait equation
Noble-Abel Stiffened-Gas equation of state
想问一下大家有没有这方面的经验，因为有很多系数，不太清楚怎么选择系数。比如：
0_1520637932983_Screen Shot 2018-03-09 at 5.25.10 PM.png
关于gamma和P_infinite的选择会影响很大

谢谢大家！

如何close这里的压力？我看到有些人在直接用空气场的压力，还有的是用stiffened-gas equation of state来解压力

这篇文章求解的也是双流体模型，但是求解策略非传统方法，即采用Operator Split的方式求解，也就是3. Numerical Methods中所说的fractional-step。不同领域的叫法不同，但思想相同。这种方法（暂且称之为新方法）求解的也是动量方程。 和传统方法的区别在于：

传统方法直接求解动量方程 新方法第一步只考虑对流扩散（方程17），然后考虑源项ODE（方程18），类似拉格朗日

新方法中压力通过方程2计算而来。

这和工作习惯有关，我这面接触的课题组里面大家都没怎么用doxygen，导致我也没怎么用，不过我这面有人用Qt，但是我个人用的是VIM。但我知道好像大家说doxygen很好用，Qt在国内也很受欢迎。对于找代码，我非常低级，就是从类的继承关系找..需要哪个打开哪个

是啊，相信CFD在国内也会有不错的发展！

是的你说的很正确。

东岳大神说了开发的问题。我来说说我的看法？清楚物理模型的物理意义貌似是最基本的也是最容易出错的部分。openfoam方面可以参考The Finite Volume Method in Computational Fluid Dynamics。还有就是最笨的办法：读代码