但是建议做个sourceMask，能量方程可以写成
\begin{equation}
LHS=-\frac{\partial p}{\partial t}\times (sourceMask-1)-S_c\times sourceMask
\end{equation}

of5.0里面已经有fvOption了，

tmp<fvScalarMatrix> epsEqn ( fvm::ddt(alpha, rho, epsilon_) + fvm::div(alphaRhoPhi, epsilon_) - fvm::laplacian(alpha*rho*DepsilonEff(), epsilon_) == C1_*alpha()*rho()*G*epsilon_()/k_() - fvm::SuSp(((2.0/3.0)*C1_ - C3_)*alpha()*rho()*divU, epsilon_) - fvm::Sp(C2_*alpha()*rho()*epsilon_()/k_(), epsilon_) + epsilonSource() + fvOptions(alpha, rho, epsilon_) );

或者复制模仿realizableKE的格式，赋值于epsilonSource() 或者 kSource()．

@李东岳 在 turbulence里k设置的问题 中说：

你这个可以通过codedfixedvalue来实现，在论坛里有代码你找找这个关键词。我现在手机写不了代码，找不到再回帖。

还有一个方法上是用timeVaryingMappedFixedValue这个boundary condition。

参见tutorials/incompressible/simpleFoam/pitzDailyExptInlet。

需要在constant里建立一个boundaryData的文件，一般只要给定0 时刻分布，就相当于给定constant分布。

boundaryData中点的坐标和具体数值可用python或matlab自行生成。

必须注意的是，boundaryData里面文件的具体格式，不同版本OpenFOAM可能不同。

如果是vof，目前还没有植入能量方程。
如果是TFM，扩散传热系数需要给定。

docker 好像没有这么多问题。

不是，这里很少有用OpenFOAM实现的

@李东岳 确实跟fluent相比有很多地方很麻烦，特别对于非专门搞cfd的人，效率有点低。个人感想，用了openfoam之后感觉一下子从理论层次到了实践层次。收获很大。
另外感谢一下李博士的dyfluid，在我入门的时候给了我很大的帮助

测试一下

@enthusiasm qq168167704 加下好友交流一下

应该是正在写“4.9”，心里打着草稿后面要写安装5，然后5就溜出来。。。一边聊天一边写字很容易把说的听的话写出来:confused:

我刚换virtualBox，这几天Vmware总是有问题，我这个是4K显示器，感觉水土不服。我试试virtualBox怎么样。

sudo apt-get install compizconfig-settings-manager compiz-plugins-extra

然后运行。

@wwzhao 谢谢。这个问题我现在已经搞清楚了，是为了增加矩阵的“对角占优”特性，提高收敛性。经过实验，假如修改成不加nu_t的公式，对于一些Co数较大的case就有可能计算发散。

大家好，我在用动态铺层的网格做模拟，模拟的物理过程为：
在一个用来冷却的容器中，每单位时间内有固定质量且固定温度(1640)的液体进入该容器中（认为不带入动量）并进行冷却。简单的模型如下：
0_1517960921550_图片1.png

为了实现这一过程，我使用了动态铺层的网格模拟液体的增加，并且当有新增网格出现时，将新增网格的internalFieldValue进行人为的定义：
0_1517961307920_1517960362(1).png
从而避免新增网格的值是周围网格值差值计算得来，实现每单位时间内有固定质量且固定温度液体进入的物理过程。
但是，我计算出来的结果显示：流场内的最高温度超过了1640度，并且随着网格数目的增加，最高温度越来越高，我做的是冷却过程啊！
哪位大神知道我这是出什么问题了吗？

这个改动主要是为了区别 const 和 non-const 访问。

@李东岳 谢谢 看来得一步步挖到底层了

邱老师现在玩的大项目，都是上亿网格的工业级别 :cheeky:

OpenFOAM对于这种切应力有没有什么计算的方法。

剪切应力的定义是固定的，用什么软件的方程都是一样的。你要从流场中获得$\tau$？印象中OpenFOAM里面没有计算tau的现成的工具。如果你不确定，你可以和Fluent比较其他量看看。

如果你要计算tau，你可以写个后处理小工具。比如这个

bool Foam::functionObjects::Lambda2::calc() { if (foundObject<volVectorField>(fieldName_)) { const volVectorField& U = lookupObject<volVectorField>(fieldName_); const tmp<volTensorField> tgradU(fvc::grad(U)); const volTensorField& gradU = tgradU(); const volTensorField SSplusWW ( (symm(gradU) & symm(gradU)) + (skew(gradU) & skew(gradU)) ); return store ( resultName_, -eigenValues(SSplusWW)().component(vector::Y) ); } else { return false; } }

https://coding.net/u/dyfluid/p/OpenFOAM-5.x/git/blob/master/src/functionObjects/field/Lambda2/Lambda2.C

- 这个操作符重载后你进行的操作不匹配-操作符的参数，你可以把你的代码贴在这里，我看看。

https://coding.net/u/dyfluid/project

没啥特殊的，就是把github上面的OpenFOAM拽过来了。国内访问比较快。有时候紧急看一下代码，打开网站随便看一下就行了。不用开linux

注册用户可以创建2个私有项目，我用了半个多月，还挺满意。速度非常不错。主要是github私有项目要收费。感兴趣的可以试试这个平台。不了解的可以百度一下版本控制。