@东岳 谢谢。这个好像是针对2.2版本的。改一下还是可以用在2.3版本。

感觉是个非常简单的小问题，目前你应该已经解决了。。

@散漫守望2016 openfoam3.0.1下，自行设置了一个算例，入口速度2m/s， 利用simpleFoam求解器求解流场的话，得到的流场从数值大小角度来看是合理的；而同样的初始及边界条件（网格也采用同一套网格），换个求解器DPMFoam单纯求解流场（关闭颗粒发射），得到的区域内最大流场速度却达到了惊人的百米量级 ，也用其它版本的of算过，DPMFoam求解出的流场数值都在百米量级左右。结果确实怪异！

DPMFoam求解器求解流场也有什么限制吗？

你可以看下这篇文章
Direct numerical simulation of near nozzle diesel jet evolution with full temporal-spatial turbulence inlet profile

我先做一个从水中剥离一氧化二氮的模拟，准备修改bubbleColumnEvaporatingDissolving。 但是我的收敛有问题，然后我有些不太理解的地方，在初始值的文件夹里water.liquid, water.gas, air.gas, air.liquid到底代表什么？water.liquid + air.liquid = 1?

还有我的T一直不收敛。谢谢大家帮忙

罗宾边界条件是对流扩散方程的一般形式的绝缘边界条件。这里，边界处的对流和扩散通量总和为零：

0_1522589151720_a1f0a086-969a-48dd-ac44-d23c64e49ea4-image.png

其中D是扩散常数，u是边界处的对流速度，c是浓度。第二项是菲克扩散定律的结果。
求大神们帮忙看看！

那我等着毕业加入你们团队@东岳

谢谢哈!

esi最新发布的isoAdvector方法碉堡了

如果你说的是这个，那么在OpenFOAM-1712里面有。

有人碰到过类似的情况吗？有没有什么方法可以缓解一下呢。

@zhuyuxuan https://www.cfd-online.com/Forums/openfoam/76331-understanding-fvm-sp.html

@李东岳 正在重新学C++，OF的程序看起来比普通的C++高级点

大家好，我在CFD online post了一个问题，希望大家支持一下，谢谢！
https://www.cfd-online.com/Forums/openfoam-solving/200287-size-distribution-models-openfoam-4-x.html

http://openfoamwiki.net/index.php/Installation/Linux/OpenFOAM-5.x/CentOS_SL_RHEL#CentOS_7.4_.281708.29
最近也正在安装OpenFOAM，找了很久找到一个比较靠谱的教程，相信你需要的资料上面都会有...
呃...这帖子是16年的....

可以做个clip，切除球内，然后在球内做vector滤镜，就只有液滴内部的矢量了。

@浪迹天大 functions好强大:surprised: ，赞一个，不过题主貌似还想要复杂的表达式改其它文件，这个就麻烦到我宁可手动调了。:mihu:

@txwwbz
OF5指南我从去年10 1就开始弄，一直很零散的弄，导致进度很慢。

2014年出OpenFOAM指南的时候，当时在读博士，无独立科研能力，时间很多，导师不给我活，我就闲着，索性翻译了。

现在有了独立科研能力，一直都在搞新的算法，指南一直没更新。

我看看下周五之前，发布OpenFOAM指南5.0

嗯，之前属于设置上的疏忽，现在已经更改，参见 http://www.cfd-china.com/topic/223/

“CFD中文网”为了知识的可追溯性，在用户发帖10分钟后不可更改以及删除帖子，因此在发帖、回帖之前，请对内容仔细考量，有删帖需求请联系我；

类似制度参考小木虫和cfd-online

@izumi 感谢。

@cyw 可以在场变量更新的地方添加correctBoundaryConditions();如U.correctBoundaryConditions();目的是可以让进程间边界面上的值更新；具体可以参考 https://www.cfd-online.com/Forums/openfoam-programming-development/88050-how-handel-fields-parallel-computation.html 。希望有用

明明是和zeroGradient等边界条件放在一起的，应该是个基础边界条件吧？为什么在使用的时候显示unknown呢

小白试了你更新的代码，运行时出现下面的错误，请问是我编译位置错了吗0_1521630956325_error.png 图片地址)

@sungda 高手，默认为8，调节为4就好了 :quwan:

如题，需要做一个全反射边界，希望有库可以参考

我的意思是在VOF里面无法加入破碎，VOF是直接模拟不需要添加模型。如果要添加破碎模型，不能用传统的VOF。

我在外面，回去看看。

您好！请问您还有当时无反射边界条件那段的代码吗？我想借看一下，理解后改一改。谢谢

@散漫守望2016 好贴，支持

@李东岳 谢谢。我将原代码放入了另外一个求解器中就没有了这种问题，可能是路径的问题

自己定义了一个类，如图，通过该类求解出来的数据在进程间不连续，我应该怎么写进程间的通信，有没有大神碰见过和我一样的问题，或者openfoam内有没有可以参考的MPI案例，求解答

0_1521016902124_TIM截图20180314163901.png

再更贴一下，openfoam中矢量数组可以新建为矢量列表形式，初始化方法为List<vector> A(n,vector::zero),希望对和我用同样疑问的人有帮助

@李东岳 感谢东岳大神的翻牌

参考你的帮助，刚刚自我捅咕解决了。。。3q

@李东岳 和icoFoam本身差不多，我可以传给您看一下，你帮我看看怎么改合适。真的困扰我很久了。 0_1520988896448_demicoFoam2.zip

@李东岳 ok,我来尝试下

哪位老师能解答一下吗

0_1520841411082_error.png

@李东岳
谢谢东岳老师！明白是什么意思了！

但这样你就无法附加除了zeroGradient的边界条件了？

动态接触角边界条件应该可以

@李东岳 谢谢李老师！

@tidedrinker 多谢回复，不过没怎么看明白，我再研究下

我在做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的选择会影响很大

谢谢大家！

= =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方程加扰动，对结果影响应该不会很大。