@李东岳 有两个问题，第一，计算稳态流场的情况下buoyantBoussinesqSimpleFoam可以采用这种方式吗？第二，这种方法是不是没计算离散相方程。

用多重网格算法求解线性方程组Ax=b，算法求解过程中按理只需要给定系数粗细网格的系数矩阵就可以了，但在看算法理论时，都涉及到粗细网格的位置关系的判定来确定插值矩阵、限定矩阵。这个就需要在算法上导入网格具体的物理位置了吧。这样好像不对（按理是只要给定系数矩阵A，b,就能求解Ax=b）对于这点一直想不明白:mihu:

@hy1112006 是的李老师，在cfd-online上找到试了一个方法，虽然有系统警告，但暂且管用：
cafe1fe5-612f-45bc-adac-914ea277d477-image.png

警告如下，目前不影响运行：
7997c7d3-2fcc-47c1-a44b-c92d5a4c78f1-image.png

@李东岳 李老师真是无私啊，正好公司想把两台工作站并联，你这直接有现成的教程！

这个求解器应该是调用了原来位置的相对路径下面的createAlphaFluxes.H文件。你看一下Make/options，换成绝对路径之类的

@Try-harder 请问这个问题解决了吗

@李东岳 谢谢李老师。

这应该是TJunctionFan这个tutorial吧，这个确实能得到局部的速度增加，但是如果管道足够长的话，管道的入口和出口处仍然是相同的速度和流量。

回复您之前的问题，是inlet，outlet这两个地方。我想实现的结果就是在outlet的速度比inlet速度高，这确实是一个不符合物理的问题。。流量不守恒。没有外力的情况下是不可能实现的。硬要做的话，好像是可以通过添加momentum source实现，但是确实是一个非物理的问题，打扰了。。

@李东岳 谢谢老师的回复，我再排查一下错误！

@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山青 您好，请问您这个问题最后解决了嘛？