CFD中文网

@是嘿嘿牛啊 密度是被动地，是根据流体种类、温度和压力变化的，你可以设置上下两层温差或不同的流体种类

二变量高斯分布：
\begin{equation}
f(u,v)=\frac{1}{2\pi\sigma_1\sigma_2\sqrt{1-\rho^2}}\exp\left(-\frac{1}{2(1-\rho^2)}\left[\frac{(u-\mu_1)^2}{\sigma_1^2}-\frac{2\rho(u-\mu_1)(v-\mu_2)}{\sigma_1\sigma_2}+\frac{(v-\mu_2)^2}{\sigma_2^2}\right]\right)
\end{equation}
MGF为：
\begin{equation}
m_{i,j}=\exp\left(i\mu_1+j\mu_2+0.5(\sigma_1^2i^2+2\rho\sigma_1\sigma_2ij+\sigma_2^2j^2)\right)
\end{equation}
纯矩计算方法为
\begin{equation}
\begin{split}
m_{0,0}&=1\\
m_{1,0}&=\mu\\
m_{2,0}&=\mu^2+\sigma^2\\
m_{3,0}&=\mu^3+3\mu\sigma^2\\
\end{split}
\end{equation}
假设$\mu_1=10,\mu_2=20,\sigma_1=\sigma_2=2,\rho=0.5$，有纯矩：
\begin{split}
m_{0,0}&=1\\
m_{1,0}&=10\\
m_{2,0}&=104\\
m_{3,0}&=1120\\
m_{0,1}&=20\\
m_{0,2}&=404\\
m_{0,3}&=8240\\
\end{split}
同时有混合矩
\begin{equation}
m_{i,j}=\exp\left(36\right)
\end{equation}
不行，混合矩计算方法不对

@LuoS 在 snappyHexMesh对复杂地形添加边界层报“浮点数例外 (核心已转储)”错误 中说：

教程中似乎有个建筑风的案例：
/Tutorials/inco

好的谢谢 我试一下

@tens 好的，谢谢老师

有大佬帮助查看一下case，调一下收敛吗。有偿感谢

@yingqing 不知道为啥，有一次成功了，在弄就不能成功了。节点的位置进行了设置，但网格的节点与我设置的节点不重合。就像这样
177924ad-d438-468f-8cf5-a612fe049465-image.png

866f0481-1f41-4269-96e1-c71fc6d8d2e4-image.png
如题，设定好流入和流出边界后，如果只改变入流的攻角，但不想改变三维网格，在90°与-90°等大攻角情况下是否依然可行？因为此时的入流风已经与流入、流出壁面垂直。
2.在不同攻角情况下，翼段升力系数、阻力系数的参考面积是否为该方向为法向的投影面积。例如升力方法的投影面积为翼段弦长和展向长度形成的矩形面积
3.参数化后能否在每次循环中批量输出xy升力系数曲线等

感谢各位大佬指教！！！

您好，请问我用fluent计算这个算例，为什么整体gas hold up相差一个数量级。入口气体速度设置的0.0213m/s。设置基本和论文里相同。
Snipaste_2025-03-24_20-32-07.png eb5ceff2-0be5-4eb1-acbb-826f3f732768-image.png

@西湖冷月 这就是天秀

各位老师好！
报错信息：
Screenshot from 2021-11-20 14-31-59.png

solidParticleCloud.C的代码：

solidParticle::trackingData td(*this, rhoInterp, UInterp, nuInterp, g.value()); Cloud<solidParticle>::move(*this, td, mesh_.time().deltaTValue()); if(mesh_.time().value()> td.cloud().tInjStart_ && // mesh_.time().value()< td.cloud().tInjEnd_) {this->inject(td);} void Foam::solidParticleCloud::inject(solidParticle::trackingData &td) //. { //. label cellI = -1; //. label tetFaceI = -1; label tetPtI = -1; mesh_.findCellFacePt(td.cloud().posB_, cellI, tetFaceI, tetPtI); if(cellI > 0){ solidParticle* ptr1 = new solidParticle(mesh_, td.cloud().posB_, cellI, tetFaceI, tetPtI,td.cloud().dB_, td.cloud().UB_); Cloud<solidParticle>::addParticle(ptr1); } }

3d408cf3-836a-4aef-bb0b-a370d462f7e5-image.png 如图，左边的上部是大气边界，右边是封闭壁面。如果不在左边画和右边一样的细网格，得到的网格质量很差。如果画了细网格，则计算速度很慢。有更好的处理方式吗

http://www.cfd-china.com/topic/3005 看看这个 我不太熟悉超临界流体

@Albert-Lee 我跟你遇到了相同的问题，你的测试和我的结果也相差不大。目前我也在寻找一个解决方法，我看到的在UDF中获取Level set function 是使用 C_LSF(c,t)宏，其中c = cell index, t = mixture phase thread，原链接：https://www.cfd-online.com/Forums/fluent/144991-access-level-set-function-udf.html

注意：需要额外include "sg_ls.h" 因为他没有包含在udf.h里面

@冠竹 有可能是入口速度太小了

请教各位，我使用fluent计算二维LNG储罐发生密度引起的液层翻滚以及气象空间超压过程。气液之间使用VOF模型，液相内部使用组分输运模型，也就是多相多组分模型。其中气液相变之间使用Lee模型udf进行。

QQ20250311-113118.jpg

在计算的过程中我这么设置：假设流体是甲烷和乙烷的混合物，并且液相的密度是随着温度和组分发生变化，随着乙烷组分的增加密度增加，随着流体温度的增加密度减小。当流体温度小于甲烷饱和温度就二者都发生冷凝，当流体温度大于甲烷饱和温度小于乙烷饱和温度就甲烷发生蒸发乙烷发生冷凝，当流体温度大于乙烷饱和温度就甲烷和乙烷都发生蒸发。

在我模拟的过程中，开始未发生多组分液层混合过程，气相空间的压力是逐渐变大的，但是当混合完成后，虽然气相和液相空间的温度都在上升，液相空间的密度下降，但是气相空间开始压力下降了，我不明白这是为什么，所以想请教下各位大佬。

0DDFF043D7976A246C818A596DF9F6CD.png

同时，在我进行Lee模型使用过程中，我发现有人使用Hashemi-Wesson这个表面蒸发公式进行蒸发率的计算，Lee模型的量纲是kg/(m^{3} \cdot s)，而这个公式的量纲是kg/(m^{3} \cdot s)。那这个纯蒸发的模型是怎么应用在fluent中，我也想请教一下大佬们。

QQ20250311-094104.jpg

@sibo 并行限制2核，网格限制是52万

方便的话加下我微信（18985935252），咱们交流一下啊

已解决

@Gordonaero 感谢，这两篇论文之前没有看到过。
这个动态填充的研究是想预测填充以及降温时间吗？关于这方面管道预冷的研究，之前看过国内的西交大有做这方面的研究，包括实验与数值计算。
按照我的理解，这个均相模型将混合物作为一种赝流体，封闭的关键其实是混合物的状态方程以及如何将混合物拆解成气液两相的比例。

:143: