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打开就可以了

低雷诺数模型(例如v2f，SST k-w等)在精细网格的情况下可以完全求解壁面的流动，那为什么还要壁面函数(比如v2WallFunction、fWallFunction等)？

@CFDngu 那你需要在参会的人里面物色个男朋友 :expressionless: