@jasper-0 更新一下目前进度，我尝试在KinematicLookupTableInjection.C里将injectors_设置为MUST_READ_IF_MODIFIED,

template<class CloudType> Foam::KinematicLookupTableInjection<CloudType>::KinematicLookupTableInjection ( const dictionary& dict, CloudType& owner, const word& modelName ) : InjectionModel<CloudType>(dict, owner, modelName, typeName), inputFileName_(this->coeffDict().lookup("inputFile")), duration_(this->coeffDict().getScalar("duration")), parcelsPerSecond_(this->coeffDict().getScalar("parcelsPerSecond")), randomise_(this->coeffDict().getBool("randomise")), injectors_ ( IOobject ( inputFileName_, owner.db().time().constant(), owner.db(), IOobject::MUST_READ_IF_MODIFIED, IOobject::NO_WRITE ) ),

但是在运行solver的时候提示不支持自动重新读取，有没有大佬知道这种情况咋解决呢？

GlobalIOList<kinematicParcelInjectionData> kinematicLookupTableInjection constructed with IOobject::MUST_READ_IF_MODIFIED but GlobalIOList<kinematicParcelInjectionData> does not support automatic rereading.

应该是drag太大了，颗粒没弹起来

只能祝收敛..

@一二 在 LES介绍的文章的一个公式 中说：

嗯，下面是我推的
\begin{equation}
\overline{\mathbf{S}}:\mathbf{\tau}+C_{e}\frac{k_{sgs}^{1.5}}{\Delta}=0 \
left = \overline{\mathbf{S}}:[\frac{2}{3}k_{sgs}\mathbf{I}-2\mu_{sgs}\mathbf{dev}(\overline{\mathbf{S}})]+k_{sgs}^{1.5}\frac{C_{\epsilon}}{\Delta}
\end{equation}
因为
\begin{equation}
\mu_{sgs}=C_{k}\Delta\sqrt{k_{sgs}} \
= \overline{\mathbf{S}}:[\frac{2}{3}k_{sgs}\mathbf{I}-2C_{k}\Delta\sqrt{k_{sgs}}]+k_{sgs}^{1.5}\frac{C_{\epsilon}}{\Delta} \
= \sqrt{k_{sgs}}(\frac{C_{\epsilon}}{\Delta}k_{sgs}+\frac{2}{3}\mathbf{tr}(\overline{\mathbf{S}})\sqrt{k_{sgs}}-2C_{k}(\mathbf{dev}(\overline{\mathbf{S}}):\overline{\mathbf{S}})) \
= ak_{sgs}+b\sqrt{k_{sgs}}-c\
= right = 0
\end{equation}
其中
\begin{equation}
a=\frac{C_{\epsilon}}{\Delta} \
b=\frac{2}{3}tr(\overline{\mathbf{S}}) \
c=2 \Delta C_{k}(\mathbf{dev}(\overline{\mathbf{S}}):\overline{\mathbf{S}} \
\sqrt{k_{sgs}}=\frac{-b+\sqrt{b^2+4ac}}{2a}
\end{equation}
当为不可压缩流体时$tr{\overline{\mathbf{S}}}=0$,那么$b=0$、$c=2 \Delta C_{k} (\mathbf{dev}(\overline{\mathbf{S}}):\overline{\mathbf{S}})$,就可以得到$k_{sgs}=\frac{c}{a}=\frac{2C_{k}||\overline{\mathbf{S}}||^{2}{}\Delta^2}{C_{\epsilon}}$

@yang山青 您好，请问您这个问题最后解决了嘛？

你这个看起来就是压力波的传递导致的。近音速可能会这样。openfoam有无反射边界条件，不过用起来效果一般。你可以试试调整出口的边界条件，也可以试试把计算域延长一些

@soulx7 可以参考这篇文章实现：OpenFOAM｜20 自定义边界条件

@dingcy 关掉耦合应该是oneway，也可以这么用，要是DPMFoam满足你使用要求的话就可以。

@李东岳 :
太感谢李老师了，明白了，谢谢。也谢谢曾帮过我的朋友。
一个玩笑：
把CFD大佬比为大城市美女，那咱乡下穷傻小子的心情就如那句歌词：
。。。would it be a sin, when i cannot help falling in love with u ........（有可能记忆有出入）

谢谢李老师

@李东岳 谢谢李老师回复，麻烦您了。我知道我哪里错了是我太粗心了。

@李东岳 可以可以，我找时间试试看。

还真没这么玩过。看看其他大佬怎么说

blastFoam这个求解器里面的算法还挺有意思的。求解器的作者也很厉害。我一直在关注。

@李东岳 效果非常好，非常感谢李老师:140:

自问自答，我粗略的整出来了，如图1，要注意画红圈地方的选择:chitang:
746340cd-5bc6-47e8-a92f-b507dd49d9a8-image.png
图1

@wangfei9088 你好，哈哈昨天晚上搞得太晚神智有点不清了，刚改了编译通过了，多谢

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

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

大家好，目前我可以输出温度场的三维数据，比如
QQ截图20221222151315.png
我想垂直于z轴做100个等间距垂直切面，在每个切面上对温度求积分并取平均值。
\begin{equation}
\int_{A} T dx dy /A
\end{equation}
然后沿着z轴画出一条温度平均值曲线。

谢谢大家！

@李东岳 李老师您好，D()是一个volVectorField。主要是发现建立两个fvMatrix时，前面一个fvMatrix会对后面的fvMatrix产生影响，主要体现在源项上。按理说这是两个独立的fvMatrix对象，不应该互相影响啊