请教各位，我使用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

在parallel setting里把MPI Type从默认改成msmpi，我的问题解决了

试用fluent 2024r1版本，运行计算时，发现一个时间步内节点0上，DEFINE_DELTA宏要被执行两次，而在2022r1版本中，在节点0上则是正常执行一次。
在DEFINE_DELTA宏中，通过I_AM_NODE_ZERO_P来只在节点0上输出变量值到txt中，理论上是应该一个时间步只执行一次的吧；在2024r1中却被执行两次，且每个时间步都出现了警告信息：Warning: physical time step size after mesh modified = 4.5000e-5s，不知道这个警告信息所说明的问题是不是和执行两次DELTA宏有联系？（同一个算例，2022中没有警告信息出现）

两相流传质有专门的宏，不需要申请UDM，找到Euler两相流的传质宏，挂载到mass transfer那个选项里就行

各位大佬，想请教下模拟搅拌喷气，气液用欧拉双流体模型，固体颗粒用dpm模型，最终颗粒受力是怎么计算的？是气固的力加上液固的力吗？欧拉欧拉双流体模型流体是相互渗透的，那么在计算颗粒受力是到底是哪一相的力还是说两相的力

已解决

如题，我现在申请了服务器节点，想进行fluent的计算，但是服务器不能安装fluent，我能不能在本地fluent软件中提交任务到服务器节点上呀？有没有大佬教教我，或者有相关的视频学一下:135:

fluent LES模拟螺旋列板抑制涡激振动时高度0.1D和0.2D的螺旋列板升力系数幅值在0.1和0.15左右，但0.15D的螺旋列板的升力系数幅值却到了0.5。三个都是同样的网格划分与时间步，为什么差异会这么大？

微小是多大尺度，另外你的物理场里边合并破碎主要由什么效应主导的需要明确一下，luo模型只考虑了小于气泡尺度的湍流和气泡碰撞引起的气泡破碎，这个气泡需要大于柯尔莫戈夫尺度

@gemini 谢谢您，服务fluent版本跟调试用的是一致的，我去试试修改成您说的这样

考古来了 想看看大佬的论文

计算的是气液混输泵的气液两相流计算，air-water(25℃)的混合介质，air为离散相给定体积分数和固定的气泡直径，water为连续项，湍流模型为SST k-w，欧拉双流体模型。计算出来之后总感觉残差是不是波动有点大了。
有没有这方面的大佬教教我，是不是设置有问题还是其他的问题。下面是计算得出的残差
47a7ac58-69fe-4db7-a71d-e8ba6fd07ebc-image.png
800b8c63-4b22-47a5-9131-cc03980b1291-image.png
05ca6207-15dd-45ab-bb5e-9628346b8a9e-image.png
e252f689-2fd6-4d82-b46a-77c7453a624a-image.png

udf有问题，你试试这个

#include "udf.h" DEFINE_PROFILE(inlet_velocity, t, i) { real x[ND_ND],y,vel; face_t f; begin_f_loop(f,t) { F_CENTROID(x, f, t); y = x[1]; if(y>0.0393) vel = 29.5; else vel = 29.5 * pow((y / 0.0393), 1/7); F_PROFILE(f, t, i) = vel; } end_f_loop(f,t) }

@深蓝 请问您找到输入粒径分组的方式了吗？我最近也要做这方面，如果您看到了请您教我一下

@bestucan 厉害，第一次知道这种用法！:146:

请问各位，我在做Re=3900的三维圆柱绕流时发现升力系数和文献差不多，但阻力系数却一直在0.6附近，和文献中1.1差距很大，这是什么原因啊？

想解决IV型激波，入射激波如何通过边界条件设置产生

在使用Fluent 6dof+overset+隐式更新时，发现写出的motion history中的重心坐标与实际明显不符，差的很大，而在udf中通过DT_CG得到的却是正确的重心坐标，不知道有没有人遇到过这样的问题，是什么原因导致的呢？

跟出口速度有关系，延长出口计算域试试？

如题，在FLUENT的用户手册里写到可以在单元区域条件里打开，但是我的2024R2似乎没有或者没找到，我用的是多相流模型

417c3160-424b-42e4-bbf0-0f81e52c3b97-image.png 说