如果单纯的mules求解每一相是无法保证alpha1+alpha2=1的

现在我就想确定，alpha2通过对应方程推导，用MULES求解不可行的吗？

MULES的出发点是保证变量的有界。所以如果你用MULES求解alpha1，那么理论上alpha1是有界的。alpha2的求解没有调用MULES，他是通过1-alpha1算出来的，如果alpha1有界，alpha2也有界。

如果你用算法同时求解alpha1和alpha2，如何处理耦合？这俩个变量是耦合在一起的，就像速度和压力。你不能单独的去分离求解，如果分离求解就就需要迭代。迭代就导致计算速度变慢。因此现存大厂据我所知都是只求解alpha1，然后alpha2=1-alpha1。

我在去你年底验证了MULES，理论上是可以保证有界，但是真实计算的时候，还是会越界。尤其是物理模型比较复杂的时候。个人觉得这方面内容搞出来，绝对是个好文。

邱老师现在玩的大项目，都是上亿网格的工业级别 :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私有项目要收费。感兴趣的可以试试这个平台。不了解的可以百度一下版本控制。

@aeronastro 您好，使用wallgradu求出的结果是个张量，要提取法向的乘以miu吗？还有一个问题请教您，wallgradu的源代码在哪里？

@李东岳 对，我做喷雾燃烧，喷雾用拉格朗日，所以计算量不大，主要的计算量就在化学反应矩阵的求解上。很多都在在化学加速的算法，比如最近一年OF更新的TDAC，就是一个加速的方法，用于几千个反应的化学动力学机理。

@bestucan 多谢指点，已经按这个方法弄好了！
@李东岳 也感谢东岳

找到了浮力项了

哦，这个，OpenFOAM编程，我计划加在OKS课程中。那个pdf可能不会写了，当时读博士时间很多，现在时间很紧 :big_mouth:

paraview中contour中能画虚线么，求大神指导

谢谢啦。

0_1517028217409_捕获.PNG

计算了几个时间步，我感觉你应该看看提高网格质量。如果你不关心这个圆圈部位的数据，应该也没问题。

@李东岳
您好，李老师，我现在的做法就是不用Glyph来显示particles，直接在paraview中load 原始particles，再用clip的box来划分出一个区域，这样在information中，number of cells=number of pints的值就是这个box中 particles的真实数量。
但是现在出现了一个新问题，我知道单个 particle体积，这个box中particles的数量，nparticles的值我也知道，这样通过计算 ：particle的体积particle个数nparticle就计算出这个box中particles的体积数，用这个数除以Box的体积就应该能算出这个box中particle的体积分数，即 alpha（volume fraction of particle). 但是算出来的这个数总是比我初始设定的particle体积分数要小很多，请问您知道是什么原因吗？

计算不准这个问题，你或许可以看看改变下网格。我之前看过一些文献，这种圆柱绕流用的都是那种园网格。比如这个帖子里面的讨论：http://www.cfd-china.com/topic/664/cd计算不准的问题-les-re3900-三维圆柱绕流-pisofoam/6

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

你用的什么版本，试试把U.boundaryField()改为U.boundaryFieldRef()

据我了解OpenFOAM里面不能同时实现俩个，不过实现应该不难，你需要编程把两个类库结合起来。纯编程问题，因为动网格算法已经有了。

@yuan 非常感谢:happy:

@txwwbz 没有了，直接用paraview打开就可以了