这个求解器应该是调用了原来位置的相对路径下面的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;
}

@WeakForm onegroup

two-group相比计算精度提升大么

这个太不好说了。你要是问Ishii，那肯定提升大，哈哈。感觉跟具体模型有关。

你得壁面函数用的哪个壁面函数？标准壁面函数还是Spalding那个

@李东岳 谢谢李老师。

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

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

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

@李东岳 原文作者是用CFX做的:136:

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

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

@WangQi-0 openfoam里面存在一些多计算域模拟，比如chtMultiRegionFoam。

就是想把界面的以及界面以上的液相电势全部设置为0

你具体这个是算什么？求解液相电势方程么。上面是空气还是其他的？

Info改成Pout试试？

@Hope 多谢多谢！我也研究研究。:duang:

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

phi *= fvc::interpolate(rho);

@fangyuanaza :146: :146: :146:

@miying :146: :146: :146:

各位专家学者，各位老师同学好，我在运行Fluent时，准备大涡模拟将一个计算域的出口参数，作为另一个计算域（也是大涡模拟）的入口条件，请问有教程推荐，或愿意不吝赐教大佬，感激不尽

@李东岳 感谢李老师分享，看了这篇文章主要适用于RANS模型，目前对这块还没有研究。我主要研究的是LES大涡模拟入口湍流生成方法，由于入口湍流生成方法有很多种，生成的风场性质也都不同，湍流发展过程可能也有不一样的特点。针对不同入口风场，源项、或壁面函数、或粗糙元布置方式应该也要对应修改，使得每个瞬时的近地面速度发生预期的调整，最后要求计算得到的平均风和湍流强度符合预期值，目前感觉实现一种通用的设置方式还是比较困难。

计算亚临界水在圆管中的流动和传热，压力为20MPa，流体入口温度300k，管壁623k，想计算二维稳态下的管内的温度分布，用水的常物性参数很快就收敛了，但是分段插值改成变物性后就收敛不了，请教下是什么原因？
1.jpg

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

1.4.4 follow path

https://www.paraview.org/Wiki/Advanced_Animations

@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}}$

@沉剑 Fluent这搞得挺狠啊

fluent等使用VOF守恒形式的通病。

类似分离涡这种大尺度的涡结构，WMLES即可，涉及转拟湍流，湍流结构等，三方向的尺度需满足计算要求。

@firejet 有一篇AMSE的Li 2012年的写的Fluent 求解器的文章，你可以看下，类似2512f2b3-f387-4506-8930-f40312a95f1e-image.png

@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:

@知行 可以。之前就有，但是挂出去一年一个人都没有分享的。大家都对自己的算例非常珍惜。如果你有类似的分享。我可以直接放到 dyfluid.com

epsi=0.09* k*omega
可以

@冠竹 您好，老师。我最近弄的和你是类似的，我是准备弄这个湿蒸汽的凝结模型，在动量方程和能量方程中添加质量传输的原项，以及添加2个用户自定义的标量函数UDS。我比较疑惑这个UDS的书写格式，因为我看书上写的是直接在user-define面板上定义两个UDS就可以，并没有写啥非稳态项和对流项的udf。但运算的时候显示我定义的UDS有问题，所以请教您有啥UDS指导的资料吗或网址吗？

@yan 我想问的是壁面边界上的粘性应力，您说的这个应该不是吧？

@gemini 可能能行，但会不会治标不治本。可能UDF里写了一套逻辑，真正实现的时候给你来另外一套。