@LYT 我有个问题要请教你，在更新相分数场的时候就要用到dgdt，此时dgdt的表达式就是setfields中的形式，但是这样的话跟公式是对不上的，这个是什么原因呢，谢谢！

东岳老师在 东岳流体 中的compressibleinterfoam解析现在没有了，还有能看到这个解析的途径吗？

想请教各位一个问题：如果我想研究高温水蒸气在液态水中的冷凝及流动情况，应该使用哪个求解器呢？或者应该在那个求解器的基础上进行修改呢？？谢谢

在OpenFOAM7虚拟机上对求解器进行编译，发现只要将求解器从solvers文件夹换到其他文件夹就好编译出错，报错内容如下：

如果有哪位知道原因可以请教一下吗，谢谢!

@杨英狄 就是说想以这样的方式注入拉格朗日粒子：每个时间步注入流场的粒子的数量、位置都与当前时间步的流场有关，这样的注入方式怎么实现呢？

@yfclark 你好，你这个功能实现了吗，可以请教一下如何实现的吗？

如果想用OpenFOAM实现一个多尺度模型，即比如利用interfoam计算气液两相流，当流场中某个气泡的尺度太小时，则转化为拉格朗日粒子。
那么该如何实现在流场中某一位置注入拉格朗日粒子呢？谢谢！

如果已知流场内某一点的坐标（a，b，c），想通过插值获得这一点的流场物理量，可以使用interpolate（）函数，这一函数的定义如下：

template<class Type>
Type Foam::AveragingMethods::Basic<Type>::interpolate
(
    const barycentric& coordinates,
    const tetIndices& tetIs
) const
{
    return data_[tetIs.cell()];
}

代码中coordinates、tetIs分别表示什么呢？
如果已知点（a，b，c）如何获取coordinates、tetIs呢？
谢谢！

@疏影横斜水清浅 这种注入方式引起数值不稳定我也遇到过，我是把入口边界条件改成卷吸边界条件了，可能在物理上不太合理，但是能够保证数值稳定。

你好，你是通过哪种方式注入的粒子呢？是patchinjection吗？