@lrl3512 compressibleInterFoam可以参考这个
https://zhuanlan.zhihu.com/p/66444729
我也是看这个学的

你的phi需要乘以密度才行。目前来看你的phi就是网格乘以速度

phi *= fvc::interpolate(rho);

@fangyuanaza :146: :146: :146:

@miying :146: :146: :146:

@李东岳 是的，编译通过了。谢谢李老师！

@jasper-0 我觉得不需要

@李东岳 收到 多谢李老师

@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 ........（有可能记忆有出入）

谢谢李老师

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

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

@tens 不好意思，隔这么久回复你，求解器已经改好了，问题在于矩阵的形式不同，接下来准备验证求解器的可靠性:chigua:

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