你这个fan那面的外圈 留一点空隙 我估计会好一些 要不然质量不守恒

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

@XieXiaoyang 您好，我想请教一下，像这种 DNS / LES 槽道流的结果，我想画速度剖面，进行正则化时候的摩擦速度${u_{\tau}}$和粘性长度${\delta_{\nu}=\nu/u_{\tau}}$是怎样计算的？
我当前的做法是：使用postChannel（就像channel395那个算例一样）将空间时间都做平均之后，可以得到${u(y)}$这样一个分布，然后计算壁面上的速度曲线的斜率，计算壁面摩擦应力，计算摩擦速度：
\begin{equation}
{u_\tau } = \sqrt {\frac{{{\tau _w}}}{\rho }} {\rm{ = }}\sqrt {\nu \frac{{\partial u}}{{\partial y}}}
\end{equation}
但算出来的摩擦速度和预设的相差很大。所以最后的曲线是这样的
8544ca30-1b5b-4552-9e71-567bb4fe6134-channel1000DirectCalUtao.jpg
算了很久，对数率那里也无法和理论解对上。计算结果的整个曲线也与理论解的形状不符合。
感觉我使用的方法不太对，想请问一下您的做法。

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:

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

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

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

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

@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 你好，你这个功能实现了吗，可以请教一下如何实现的吗？

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

谢谢大家！

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

你先备份一下。

然后你wclean一下，看一下这些东西会不会消失。消失的话再编译一下

@李东岳 好滴，感谢李老师，我编一下

@李东岳 在 VOF方法中如何准确计算自由表面曲率？ 中说：

这个方法在函数calculateK()里面，可直接调用。计算结果在3232网格上能得到光滑的结果，但不对。在6464网格上出现了数值振荡。

VOF结果在网格比较细的时候结果不对？这个结论很有挑战性

李老师，我用interIsoFoam算了一下，还是不对。在32$\times$32的网格上可以得到光滑的结果，但是和解析解不一样；在64$\times$64的网格上结果出现了振荡。interfoam也跑了一下，结果和interisofoam不一样，也不对。。。（这个算例的设置参照和解析解在主贴的参考文献[1]中）

算例在这里，大家有兴趣可以跑一下。要是算例设置有问题也可以指出来。另外，我是在v2006环境下跑的(这个链接可能需要注册坚果云，我没找到其他的上传方法:xinlei: )
capillaryWave

@李东岳 并行文件夹里面的2e-08文件夹是有pointProcAddressing文件的，但是后面算出来的2e-02，4e-02就没有这个文件了。所以reconstructPar说找不到这个文件，报错了。

@exthan 您好，我刚刚也下载了OF7版本的makeAxialmesh,但好像没找到编译文件，没有Allmake文件或者wmake文件，请问您是怎么编译的呢d27bf862-45db-4120-8ba1-050e18d46128-image.png

@bestucan 谢谢老师，我试了下，还是不行~

他是由于none量 他把它删了，看来我的问题就是那个omega，只要用

const volScalarField& omega_ = ( this->mesh_.objectRegistry::template lookupObject<volScalarField>("omega") );

就会找不到，估计是结构问题，难道只能重写了。

@李东岳 李老师您好，感谢您提出的意见，现在只是算的有些慢但是似乎问题已经解决了，谢谢。33.png

@jacob-z

mpirun --oversubscribe -np 20 xxFoam -parallel

应该是进程超限锁死，需要添加申请允许超限，就可以了

还有种方法，通过一个工具将OF中颗粒坐标转换成旧版本的格式，然后就能被paraview和tecplot等读取了
https://github.com/blueCFD/lagrangianExtraFunctionObjects

@李东岳 多谢李老师，我的物理问题中速度是一定会求解的。所以放在泊松方程中能够更加稳定，避免震荡。
明白了，多谢您