CFD中文网

@youhaoyu 如果计算过程中发现粒径没有变化，可能是添加的RP方程计算程序没有被调用。参考原生的lagrange库，粒子速度的更新函数calcVelocity（），是在函数calc()里面被调用的。你可以看看你的RP方程函数是否在calc()里面调用了。

粒子性质的定义不仅在KinematicParcel.C和KinematicParcel.H，在其他文件中也涉及相关代码。你可以正则搜索一下，看看那些文件涉及粒子性质的声明和定义。

这个问题我也遇到过，比如对于同一个流场，同样的空化数，来流速度10m/s就不会产生负压，来流速度20m/s就有负压产生。增加SS空化模型的系数n的值或许能解决问题。

@杨英狄 感谢感谢！

利用compressibleInterFoam进行汽液两相流计算(空化)，流场温度是非均的。
在定义thermoPhysicalProperties时，需要选用：

1、transport模型：计算粘度和导热率
2、thermo模型：计算比热容
但是这二者的子模型只能给定物理性质随温度的变化关系，事实上这些物理性质(如粘度、导热率、比热容)也随压力发生变化。
所以如何定义thermoPhysicalProperties，才能让他们在不同的温度和压力下均能计算出正确的物性？
谢谢各位！

@李东岳 李老师您好，我在做空化，需要考虑湍动能对空化的影响，其实是在饱和蒸气压上加一个修正项：pturb=0.39*（密度*湍动能）。所以就需要在求解器里面声明一下湍动能。
感谢李老师。

大家好，想请教一下各位，如何在求解器中，或者说是在createFields.H中声明湍流相关的场，如湍动能、湍流耗散率等。

因为目前需要考虑湍流对汽液传质的影响，需要在压力方程中用到湍动能，所以想知道怎样在求解器中声明。

感谢各位！

@李东岳 请教一下李老师，怎么在creatFields或者求解器中声明湍流的相关信息呢，如湍动能k，谢谢李老师

@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吗？

各位好
若使用DPMFoam对气相粒子在液相中的行为进行计算
气相粒子的碰撞参数取多少较为合适呢

collisionModel pairCollision;

    pairCollisionCoeffs
    {
        maxInteractionDistance  0.0025;

        writeReferredParticleCloud no;

        pairModel pairSpringSliderDashpot;

        pairSpringSliderDashpotCoeffs
        {
            useEquivalentSize   no;
            alpha               0.02;
            b                   1.5;
            mu                  0.10;
            cohesionEnergyDensity 0;
            collisionResolutionSteps 12;
        };

        wallModel wallSpringSliderDashpot;

        wallSpringSliderDashpotCoeffs
        {
            useEquivalentSize no;
            collisionResolutionSteps 12;
            youngsModulus   1e8;
            poissonsRatio   0.23;
            alpha           0.01;
            b               1.5;
            mu              0.09;
            cohesionEnergyDensity 0;
        };

        U     U.air;
    }

谢谢各位！

@卡洛 是直接断掉了，算例结束运行了