CFD中文网

第一次在发帖，请问一下各位老师，就是想在openfoam中计算一个场值的laplac的值，用openfoam中自带的laplac显式计算的值的结果容易较大，希望能够调用临近单元的再一个临近单元来进行更高阶的计算来防止数值过大的问题，请问这种方法怎么实现呢，谢谢

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

cylinder.tar.xz

sdf asdf

在计算流固热耦合的两相流时，出现warning：no bounding box for zero-sized point field
虽然依旧可以计算，但是计算结果phase1-volume fraction不演化，因此应该是哪里出现了问题，想请教一下老师们，该问题如何解决，网格使用的是blockMesh画的轴对称网格，两面使用wedge边界条件
644707cd-9109-4ade-802a-773a7c3ec8ce-image.png

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

感谢各位大佬指教！！！

openFoam VOF 如何监测进出口的流量和压差

b768ef65-777f-4eff-827d-64e7647225e3-image.png ca8f10bf-0a79-46a0-868e-3508b1d64b78-image.png两相流模型为VOF，流场按照浸入边界法对流速场进行了插值，颗粒内部流速和颗粒速度一样为0，但是还是有初始颗粒内部虚拟流体流出，各位大佬有什么解决办法么:135:

各位大佬，本人想用rigidbodymotion计算双浮体铰接，铰的位置打算设置在两个浮体中间，但是目前的设置只能铰在两个浮体重心位置，我的浮体1重心位置在全局坐标系中坐标为（0，0，-0.0126），浮体2重心位置在全局坐标系中坐标为（0.5，0，-0.0126），打算把铰设置在在全局坐标系中（0.25，0，-0.0126）的位置，我该如何实现，以下是我的设置1.jpg 2.png 3.png

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

在OpenFOAM11之后，跨音速求解器rhoCentralFoam（或者说求解模块shockFluid）的求解变量变成了原始变量而非守恒变量。我还没有详细的看算法，粗略的看是OpenFOAM基金会自己做的数值改动。

这知识啊，无穷无尽。

新版：https://github.com/OpenFOAM/OpenFOAM-11/tree/master/applications/modules/shockFluid

HLL类：http://dyfluid.com/hll.html

CFD: 可压 + 密度基：http://dyfluid.com/density.html

请问拉格朗日的TomiyamaLift如何进行编写

TomiyamaLift { sigma 0.072； }

我这个写法openfoam报错了

试用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中没有警告信息出现）

33095a1e-6724-4493-b8a0-376e8d4e19a6-图片.png
请问各位老师，这是我做的颗粒对壁面的磨损，但是从paraview中看磨损在碰到壁面时为贯穿式，其中壁面通过SW建模，fluent meshing画网格。有没有解决方法，十分感谢。下面是liggghts代码：
221c09f4-6c30-4abb-9817-98cea3a76378-image.png
e4ce9c6d-8813-406b-a6ba-36e600ff8f26-image.png

66fe34e5-2b9e-4952-ab9b-9c798f0e4be5-1740230382354.png

#include "kinematicMomentumTransportModels.H" #include "transportModel.H" #include "addToRunTimeSelectionTable.H" #include "RASModel.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // makeBaseMomentumTransportModel ( geometricOneField, geometricOneField, incompressibleMomentumTransportModel, IncompressibleMomentumTransportModel, transportModel ); #define makeRASModel(Type) makeTemplatedTurbulenceModel (transportModelIncompressibleTurbulenceModel, RAS, Type) #include "myRAS" makeRASModel(myRAS); 以上是models.C 单独编译myRAS是没有报错的，编入的时候出错了

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

目前想要用多孔介质代替小孔计算浮动owc的重叠网格模型

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

有大佬可以分享一下旧版本中使用view factor模型计算的案例文件么，我装的v11找不到相关的案例了，感谢感谢

当我使用wave2foam对一个浮式平台进行六自由度运动计算时，使用变形网格的情况下，网格的变形是否会影响波浪的作用？下图为初始状态
d6b44b5037b32f6841b94ea87c9ce51.png
下图为为运行50s后液面附近的网格
2fa135c091d48299c532ad50beb1165.png
这种网格的变形是否会造成波浪的衰减呢？

各位大佬好，由于前段时间使用codedFixedValue+AMR运行报错，怀疑是codedFixedValue造成的自适应加密出问题，现在我从网上学习使用其他自定义边界条件设置。
我创建了新的边界类型，编译已经成功，但是运行后出错，附件是我的计算配置文件，真诚请教各位大佬！！
感谢大家！
e8c1f25c-8331-4a81-b430-da6c55e6888c-image.png
边界条件文件：
parabolicVelocityFvPatchVectorField.C parabolicVelocityFvPatchVectorField.H
U
计算设置文件：
controlDict