https://github.com/OpenFOAM/OpenFOAM-2.1.x/blob/master/src/finiteVolume/finiteVolume/fvc/fvcMeshPhi.C

看这个代码

注入质量

如果你指的是往计算域里面添加液体的话，设置进口边界条件不就好了？

这图应该指的是在网格分辨率足够的情况吧，如果分辨率不够，在0cell和1cell之间不会只夹着一个0，1之间的cell

计算完毕后，规定任意一条直线及点数，如何插值求出直线上任意一点的值？

特别是计算区域分为几个小的区域，每个区域内有不同的物性参数的情况如何计算？

@李东岳
的确奇怪，我现在都不知道为什么纯粹进行值之间的交换会出现问题，又没有显示命名的重复，
不过我把发生错误的变量名字改了以后，错误就不在出现，虽然解决了，但还是对于错误的原因存疑。

你好，请问你这个问题解决了吗？

你在MRFProperties里面或者fvoptions里面需要定义旋转的区域，paraview是显示的那个效果那个变不了。你可以把结果当做稳态时候的结果。或者搅拌器叶片转到某一个时刻，流场就是这个样子的。

@wwzhao 谢谢BZ！我试试~

好像以前见过一个类似的帖子，下面那个value 是没用的，主要用于paraview读取，你随便写100也可以。

这个问题通常是由于内存不够引起的崩溃，你机器的内存多大？算例的网格多大？

例如下面的两相界面处的jump condition，在表面张力系数σ不为常值函数时，沿界面方向的表面张力作用不能被忽略
[μ∇u · n] = ∇s σ
其中∇s =(I - n⊗n)∇，表示表面梯度算子
在interFoam中表面张力系数为常值函数，用CSF模型计算法向jump condition，想知道在多数CFD计算中，这种jump condition都会怎么处理?

好的 谢谢 我去看一下

@mohui 我是基于那个方程进行修改添加一些，只不过它本身存在一些问题，当我设置1的时候温度正常

可以，so？

@李东岳 前辈，我一直用pyFoam来观察残差，我猜的这个bounded variable plot出现是因为计算出现了不稳定，如果这个bounded量计算没过界就不会出现。:confused:

另外观测压力等物理量的残差如下，这个就比较好理解了:cheeky:

0_1508035857600_WechatIMG15956.jpeg

@我也不知道吧 好的，我加你。

我认为这个force应该是用来维持湍流的并不是产生湍流，否则没有这一项的话湍流会被粘性耗散掉。如果初始的速度场没有用boxTurb产生随机场的话，而是均匀的速度为零的速度场的话，根本算不出湍流。

这个我还不确定，你可以不用boxturb，然后用dnsFoam算一下他那个tutorial看看有没有湍流出来。

想问一下这一段代码是如何定义force的？

目前只知道基于快速傅里叶变换，更深的代码还没看 :cheeky:

gv只能打开一个eps

@李东岳 除了添加力，还需要求解界面上表活剂输运方程，具体的介绍（James, Ashley J, Lowengrub, et al. A surfactant-conserving volume-of-fluid method for interfacial flows with insoluble surfactant[J]. Journal of Computational Physics, 2004, 201(2):685-722. ）。怎么求解输运方程这一块我不知道怎么处理。

http://www.cfd-china.com/topic/905/几个有意思的研究方向