CFD中文网

大家好。

我最近在学习cfdem，模拟了小球入水。流体网格尺寸0.003m，小球直径0.0254m，钢制小球。界面捕捉采用VOF方法，并采用浸没边界法IBM刻画小球形状。液体（水）表面张力0.07。

入水后小球位移模拟值与试验值对比良好。使用paraview-counter设置alpha.water=0.5提取液面，发现气液界面附近的空腔结构与试验照片对比较好，但是小球周围并没有出现空腔结构。

请问是算例设置相关问题还是paraview设置问题？

文献（https://doi.org/10.1016/j.oceaneng.2021.110198）
a72bdbf6-06ad-4f91-8cb0-05d6ada9f5cc-image.png

我的模拟结果小球周围没有出现空腔结构
2e351048-86c2-47d4-9f10-b76146ea1e8a-image.png

该书11.8节（中文版305页，英文版411页）中的例2，推导源项的表达式时，高阶格式减去一阶迎风格式，一阶格式的值取上上游的值φu，为什么这里不取φc呢？先谢谢了。
7a735b6e-01c9-4088-971a-359144da8189-image.png

东岳老师，请问您的SCI文章《Simulation of bubbly flows with special numerical treatments of the semi-conservative and fully conservative two-fluid model》里模拟鼓泡塔时的入口速度是如何设置的呢？直接速度入口，表观气体速度和论文里的相同吗？文中的这句话“At the inlet, the velocity of the gas phase was calculated from the superficial velocity in different test cases;”如何理解呢？谢谢！

6c2fddb3-ec46-43b7-9757-5adb6b854726-image.png

各位大佬好，我这边采用OpenFOAM5.x版本进行计算，在场2中存在一个计算变量A需要传递到场1的重叠区域中，目前两个计算域的网格如图片中所示存在重叠，且设置的边界保持一致。请问该如何保证数据传递的可靠性，我目前使用的指令如下：
mapFields $sourcedir$ -case $targetdir$ -noFunctionObjects -sourceTime '30'
指定map方法的三种（mapNearest, interpolate and cellPointInterpolate）我也试过，效果也并不是很好。各位大佬给点建议，有没有相类似的数据处理经验。感谢感谢！

14295afe-5391-425a-88dd-94e8f07b96ab-image.png

各位大佬好，我在用vof模型计算电磁力驱动铝液-电解质两相流动，界面模型选的是尖锐/分散，初始两相界面的高度是0，计算完毕后，流速场的结果比较符合，但是两相界面大部分区域都在0之下，感觉是底部的流体变少了，请问这是什么情况？
下图是模型设置
a0a496e9-dbd7-43b5-ab21-4798a920eb1a-image.png

b200479f-f2af-4ef8-8edc-593b74a77d62-image.png 目前针对幂律流体唯一一个RANS湍流模型是16年的一个论文根据v2f模型开发的一个模型，openfoam里边带v2f模型，可以参考着改，但是phasecompressible文件夹下没有。对比了一下多相和单相的kepsilon模型，多了不少代码，而且针对多相的参考的文章也是针对多相开发的kepsilon而不是原始的kepsilon方程，想问一下如果想把这个方程应用到幂律流体内的气泡流是不是没啥戏了

大佬们在openfoam中想定义一个化学反应，例如CO+0.5O2=CO2，dcaad515-926a-4226-a7c1-d60958b9b5bd-738e45f4117e03a697f70d7b89645cf.png
怎么实现定义反应物浓度的指数呢？
dcd343de-d00e-4700-8c9e-90a9b8bd9640-f2712f0616b60be50bfbaef1929543b.png

有大佬帮助查看一下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是没有报错的，编入的时候出错了