@杨英狄 在 可凝气体的冷凝数值模拟 中说：

高温水蒸气在液态水中的冷凝及流动

不太好理解。水蒸气在水里冷凝然后流动？

后来这个解决了吗 :-)

主要是因为OpenFOAM那面多相流的工作量比较大，改进也比较多，新老版本都存在，所以就同时有了这么多东西。

twoPhaseEulerFoam: 最老版本的双流体求解器，稳健性一般 multiphaseEulerFoam: 最老版本的多流体求解器，稳健性不好 reactingTwoPhaseEulerFoam: 末期版本的双流体求解器，稳健性比较好 reactingMultiPhaseEulerFoam: 末期版本的多流体求解器，稳健性比较好 twoLiquidMixingFoam 不属于上面几类

想用IATE和类方法计算不含相变和传热的问题

用reacting类求解器

偶尔误入这个帖子，还挺有意思。

@李东岳 十分感谢李老师，后面课题组有培训经费一定要来上您的编程课:huahua:

你好。想问一下在openfoam中离散过程中对于面插值的中心差分格式的源文件代码openfoam/src/finiteVolume/interpolation/surfaceInterpolation/surfaceInterpolationScheme/surfaceInterpolationScheme.C
其中对于点乘面差值格式的代码形式主要由Foam::surfaceInterpolationScheme<Type>::dotInterpolate 函数决定的，我现在的主要想法是想通过修正插值格式在代码sfi[fi] = Sfi[fi] & (lambda[fi]*(vfi[P[fi]] - vfi[N[fi]]) + vfi[N[fi]] ) 后面加体心梯度项点乘偏移量向量，如图所示：
b2b2b5cc-a725-48e9-ad35-ba8bc36d7208-1675996388057.png

我该如何修改代码，或者说大致的思路，请各位指点，拜托了。

@李东岳 是的李老师，这个功能挺有意思的。我前不久写过Pope和Celik Index的类似后处理代码，看到之后新功能就分享给大家了。当然我实现的比较简单，适用性没有那么好，但是足够用了。您说的写到controlDict的思路也可以，这个其实也就是个简单验证，也不需要多优雅～

nut 在壁面上使用壁面函数 nutUSpaldingWallFunction，会出现上面那样的速度分布图
但另外一个${Re_{\tau}}$也等于1000的算例，得到的速度分布也差不多
我等等重开一个帖子说明下问题

@是嘿嘿牛啊 haha 理解理解

你在机械硬盘跑，数据少的话，不影响跑的速度

000.jpg

如图，是一个立方体空腔，腔体内是空气，立方体顶层有四个入口，每个入口注入的都是同种液体（比如水）。

请问怎样能够比较方便的标识出入口信息呢？就是把inlet1注入的液体标识为1，inlet2注入的液体标识为2，以此类推。（不需要考虑液体混合的问题，只模拟一小段时间，到混合前就停止了）

烦请李老师和各位大佬解答 @李东岳

@李东岳 有两个问题，第一，计算稳态流场的情况下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

捕获.PNG

看起来属于这部分内容，我之前手算过类似的东西（附录A.3 扩散系数的影响）。是，那是要多计算域求解。但是你这个还涉及到一个动边界的问题。其他sci里面怎么处理的？

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

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

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

请教大家：
openfoam2106版本中如何监测某个面上的涡量值
我的思路是在controdict中加入surfaceSampling和vorticity两个函数，发现确实时间文件里面有了涡量的结果，而且也出现了postprocessing文件夹，但是该文件夹的数据里面只有x y z坐标，而没有涡量的值，请问大家有什么办法吗，具体的Function如下

surfaceSampling
{
type surfaces;
libs ( "libsampling.so" );
enabled true;
timeStart 0;
writeControl runTime;
writeInterval 0.06;
interpolationScheme cellPointFace;
surfaceFormat raw;
fields (vorticityField);
surfaces
(
plane
{
type plane;
planeType pointAndNormal;
pointAndNormalDict
{
normal (0 0 1);
point (0 0 0.1) ;
}
bounds (0.1 0.1 0.08) (0.6 0.4 0.12);
}
);
}

vorticity1
{
type vorticity;
libs (fieldFunctionObjects);
field U;
result vorticityField;
region region0;
enabled true;
log true;
timeStart 0.25;
timeEnd 1000;
executeControl timeStep;
executeInterval 5;
writeControl runTime;
writeInterval 0.06;
}

@李东岳 谢谢李老师。

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

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

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

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

pFianl控制连续方程的最终残差来保证速度一进一出是守恒的。在瞬态的情况下要保证时间步准确，pFianl不能添加任何低松弛。PIMPLE通过松弛确实可以实现大时间步的跨越，如果外迭代有足够多的次数，可以保证大时间步后的收敛解。这样可以解释的通。但保守起见，我更倾向把fianl的松弛去掉。

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

@yuhxfoamer

首先pathlib这个模块在最终的代码中并没有被使用到，所以你可以将这个模块和对应的代码移除，这并不会影响代码运行结果。 然后谈谈你这个报错。readData.py这个文件是用pvpython执行的，这个python解释器是paraview自带的，对应的第三方库也是在安装paraview时配置好的，所以错误来源于你的paraview附带的pvpython解释器（可能是安装问题，可能是系统环境问题，也可能是软件版本问题，请自行排查），和系统python解释器无关。你使用pip和conda只是在安装你系统python解释器的第三方库。 尽管最终版本这个模块被弃用，不过这里再补充一下：pathlib在本代码中的用途是创建空文件夹（pathlib.Path），在上面的讨论中也提到了，因此可以使用os.mkdir替代。

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

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

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

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

谢谢李老师

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

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

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

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

刚接触sprayfoam，请教三个个问题：

初始场指定区域怎么布置液滴？setFields？ 液滴的喷射位置和喷射速度怎么设置为跟随inlet边界？ 我模拟的是预气化的液滴，比如C7H16既有液滴也有气体，我在算例里面该怎么区分他们？

谢谢前辈，希望解答或者提供一些解决问题的方向:140:

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

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

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

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

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

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