= =fvOptions(U)是用于处理源项的，暂时可以忽略。

SimpleFoam动量方程参考 http://dyfluid.com/icoFoam.html 中的方程1，只不过没有时间项。我建议你看看这个文档，simpleFoam和icoFoam区别非常小。

@李东岳 好的谢谢您，李老师

在模拟co2、水蒸气的对流扩散的时候，准备用[Fluent theory guide]{http://www.afs.enea.it/project/neptunius/docs/fluent/html/th/node128.htm}的species transport方程，其中对流项和扩散项都没有什么问题，主要是重力作用的加入，因为co2水蒸气都比空气重，在实际中重力作用很明显的
对于重力项之前别人推荐参考buoyantFoam，我也根据[东岳老师]{http://www.dyfluid.com/buoyantPimpleFoam.html}的这篇讲解，在OpenFOAM里面加入speices transport引起的重力变化。主要还是加入了密度场，通过不同密度对速度方程中的浮力项的作用，来实现重力的影响，并没有直接对co2对流扩散有作用。所以主要还是靠的速度变化，引起co2向下流动的现象。这种变化感觉非常不明显，比如人呼出的co2才占整体重量的4%，实际模拟中呼出co2基本不太向下聚集，Fluent在species transport中是不是也加入了源项？0_1519956843102_scalarTransportBuoyant.zip

@zxyf1 any reference?
thank you

@李东岳
man top

Line 2 shows CPU state percentages based on the interval since the last refresh. As a default, percentages for these individual categories are displayed. Where two labels are shown below, those for more recent kernel versions are shown first. us, user : time running un-niced user processes sy, system : time running kernel processes ni, nice : time running niced user processes id, idle : time spent in the kernel idle handler wa, IO-wait : time waiting for I/O completion hi : time spent servicing hardware interrupts si : time spent servicing software interrupts st : time stolen from this vm by the hypervisor

niced user process 表示优先级较高的用户进程。

用google搜索吧。

你现在可以按照提示操作一遍。如果还不行，把你的错误输入到google中去，一般都会找到满意的答案。起码我目前为止是的。

没有google就花钱买v屁恩

solution { active true; coupled true; transient false;//transonic cellValueSourceCorrection off; calcFrequency 10; maxTrackTime 100;

硬生生的debug 2个小时....:zoule:

请问为什么速度散度是0，相分数就是0到1之间呢?

http://www.cfd-china.com/topic/1315/关于interphasechangefoam的疑问

公式30的等号右侧在每一个相分数是0的时候，并不等于0啊。

等于0啊。

请问有详细的讲这方面的书或者文献吗？

这篇文章是我写的最详细的文章了。写这个文章的本意就是把twoPhaseEulerFoam解析写成一个review，但是被拒的可能很大，于是加了个算例变成了regular paper。目前我还没看到更详细的文献。

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

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

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

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

@mohui 如果使用湍流模型，那么就是模化（或部分模化）了流动的脉动。真实流动的脉动越强，湍动能也更大；但是你计算的时候是不会算出来流动的脉动的，只会得到较大的k值。
你在k方程加扰动，对结果影响应该不会很大。

@haining-luo 参考压力应该是只有在没有进出口边界条件的情况下才发挥作用，例如封闭罐体。所以你的算例里设置pRef实际上并没有发挥作用。

自己破案了……
vbox的3D加速关掉就行了
尴尬:joy:
附上链接
ParaFoam OpenGL Warning

但是建议做个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，在我入门的时候给了我很大的帮助

测试一下