@东岳 感觉还是自适应网格的问题，我去掉了dynamicMeshDict，就能算了 :136:

@Samuel-Tu upwind只有一阶精度，更对得上其他论文的数据，得查看其他原因了

现在很多研究采用的是waves2foam造波工具箱集成到OF中。我用的是OpenFoam1906版本，这套版本中有自带的造波工具箱即原来的ihFoam，所以我就不想再用waves2foam了。但是这个版本中造波后如何提取任意一点处的波高随时间的变化呢，有没有用过的大神。

目前我在controldic文件中添加了一段代码能追踪自由波面，但是提取的文件也是vtk的格式无法打开，恳请会的大佬不吝赐教！感激不尽！

@东岳 在 微分公式OpenFOAM编程 中说：

volVectorField t = fvc::grad(wbeta);
volScalarField t2 = t.component(2);

我有试过求梯度，再求分量，通过源项的方式加入微分方程中：

volVectorField t = fvc::grad(wbeta); volScalarField t2 = t.component(2); fvm::div(phi, beta)==-fvc::su(t2, beta);

这样可以吗:134: 。。

你的连接可以打开 但是我这网速比较慢 只有文章的第一页 明天我再看看吧

文章挺有意思，目前不清楚是不是CFD指导实验设计，还是CFD用来辅助实验验证。看起来确实可以用interFoam来算。

@a602657563 非常感谢，我好好研究一下

谢谢东岳老师回答，每个网格应该都是有粒子存在的，您说的第一种方法扩大出口区域，是像下面这样吗
X)O8BF$DQD@}2~IBNZ5_F2W.png
第二种把P改成P_rgh以及可压缩处理啥意思没太听明白，论坛上也没翻到相关的东西，老师能再讲清楚点吗

@炎冰灬灬 对于耦合面的主动控制我想用RBF方法，所以导致耦合变形的动网格求解器和主动控制的动网格求解器不一样

@Ivo blockMesh太费劲了 商软应该可以

(a23fd0ef-f371-4ba5-82e4-427507498568-图片.png 图片地址)

(8e592cf8-88f0-4dc1-9b93-90b03aa69a0c-图片.png 图片地址)

好的，我去试一试，谢谢前辈指点:chigua:

贴纸们，有个边界条件设置不知道怎么弄，请大家帮忙看看：
给定了一个速度场，我要模拟的是速度场中污染物浓度（即单位体积内污染物的质量）的扩散。此处污染物是用连续体来表征的（不是颗粒）。
现在我的边界条件是在入口设置单位时间内进入的质量。。但是求解的微分方程实际是污染物的浓度建立的。
怎么把单位时间内进入的质量加到入口边界上呢。。是否需要某种转化呢。。

@东岳 在 可压缩求解器为什么密度残差一直为0，迭代步也一直是0 中说：

degenerate

感谢李老师指导！:ok3:

@bestucan :146: 感谢感谢

@东岳 好的东岳老师！ :ok3:
最近也是刚开始试着改写求解器，顿感困难重重，确实需要多锻炼多积累。

谢谢老师

在creatFields.H文件中看到了如下代码，没有查到相关解读，创建类fieldTable的对象fields，然后调用了add(形参)函数，我在multivariateSurfaceInterpolationScheme.H找到了该函数，感觉是往hashTable中插入该变量的相关信息？麻烦有了解的大佬指导一番:xinxin:

Info<< "Creating multi-variate interpolation scheme\n" << endl; multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields; fields.add(h); fields.add(rho);

另外，在multivariateSurfaceInterpolationScheme.H文件中发现如下的类派生关系，两组类的派生关系这么声明是什么含义？第一次见:136: ，还望大佬解读一下:

namespace Foam { /*---------------------------------------------------------------------------*\ Class multivariateSurfaceInterpolationScheme Declaration \*---------------------------------------------------------------------------*/ template<class Type> class multivariateSurfaceInterpolationScheme : public refCount { public: //- fieldTable class fieldTable : public HashTable<const GeometricField<Type, fvPatchField, volMesh>*> { public: fieldTable() {} void add(const GeometricField<Type, fvPatchField, volMesh>& f) { this->insert(f.name(), &f); } }; ...... }；

OpenFOAM 当然更配 paraview。paraview 有个数据筛叫 celltopoint ，相反的也有，pointtocell。转了之后再输出。
tecplot用的不多，不清楚。

@bestucan 学到了，感谢

@东岳 请问李老师您能帮忙看一下这个问题吗？之前也了解到您做了许多这方面的很好的工作。另外，我在openfoam里面试着对比了一下无滑移和自由滑移的速度边界条件，发现二者的差别非常小。
另外不知道openfoam里面是不是也按照式（2）的方式，通过接触角算出界面的法向量，以此来植入接触角边界的？鄙人目前还看不懂openfoam源码，只是简单的应用，惭愧。所以向您请教。