OpenFOAM

在ofv2312上编译ofv2406代码遇到问题

bestucan — Fri, 17 Oct 2025 03:44:09 GMT

测试回复

老版本OpenFOAM E-L多cloud的解决方式

星星星星晴 — Wed, 15 Oct 2025 14:33:55 GMT

参考coalChemistryFoam（OpenFOAM4）即可。
注意在修改cloud的名称，引用，以及在各个方程中添加对应的源项。

需要注意的是引用顺序，一定要在求解flow之前，不然两个cloud见到的Euler场是变化的。
尤其在2-way coupling中非常重要。不然第二个cloud见到的flow是被改变的。

同时，这个方法也有一个问题，也就是两个cloud相互之间是不可见的。
如果考虑4-way的话，或者涉及到考虑整个cell中的theta()的话，这个方法是不正确的。
因为两个cloud是相互独立的，求解的受力也都是独立的。

关于风场模拟入流ABL边界条件和其适配的网格大小

HITSC30 — Tue, 14 Oct 2025 14:03:13 GMT

会吧，即使体心的速度不改变，太粗糙的网格也会让你的速度近似为多线段。如果你的网格足够细的话，平均风剖面应该是拟合对数率的，只是会过滤掉部分脉动成分的。

自定义边界条件遇到的问题

Joker — Mon, 13 Oct 2025 10:59:43 GMT

@xpqiu 哇，真的是太感谢您了，您一说我就明白了，我确实是没有写这行代码！

求助大佬们标准k-epsilon初始条件

xpqiu — Mon, 13 Oct 2025 05:55:26 GMT

@Noob
也不可行，因为标准 kEpsilon 某些是所谓的 High Re 湍流模型，不能用于求解边界层内的流动。假设按照你说的，y+很小，k 和 nut 都用 lowReWallFunction，并且epsilonWallFunction 开启 lowRe correction，也是错的。你可以画一个平板边界层的算例，通过调整网格的首层厚度，得到 y+不同的网格，然后你用标准 kEpsilon 模型来计算不同网格下的壁面剪切力，就会发现 y+ 小于 30 以后，结果就不对了。
但是kOmegaSST就可以在不同的 y+ 下都得到较好的壁面剪切力结果。

codedfixedvalue实现运动边界速度的指定

Joker — Sat, 11 Oct 2025 11:14:05 GMT

type codedFixedValue; value uniform (0 0 0); // 初始速度 name UM2; // 自定义边界条件名称 codeInclude #{ #include "fvMesh.H" #include "polyPatch.H" #include "fvcMeshPhi.H" #include "addToRunTimeSelectionTable.H" #}; code #{ // 避免重复计算 if (this->updated()) { return; } // 获取网格和边界数据 const fvPatch& p = patch(); const polyPatch& pp = p.patch(); const fvMesh& mesh = internalField().mesh(); // 仅在网格运动时计算（静态网格跳过） if (mesh.moving()) { // 1. 计算上一时刻的面中心坐标 const pointField& oldPoints = mesh.oldPoints(); vectorField oldFc(pp.size()); forAll(oldFc, i) { oldFc[i] = pp[i].centre(oldPoints); // 基于上一时刻节点坐标 } // 2. 计算几何速度（位置变化/时间步长） scalar deltaT = mesh.time().deltaTValue(); vectorField Up = (pp.faceCentres() - oldFc) / deltaT; // 3. 通量修正确保连续性 const volVectorField& U = static_cast(internalField()); scalarField phip = fvc::meshPhi(U, p.index()); // 网格运动产生的通量 vectorField n = p.nf(); // 边界法向量 scalarField magSf = p.magSf(); // 面面积 tmp Un = phip / (magSf + vSmall); // 通量对应的法向速度 // 4. 合成最终速度：几何速度 + 法向修正 // 公式：修正速度 = 原始速度 + 法向*(通量法向速度 - 原始法向分量) vectorField::operator=(Up + n*(Un() - (n & Up))); } // 调用父类方法完成更新 fixedValueFvPatchVectorField::updateCoeffs(); #}; codeOptions #{ -I$(LIB_SRC)/finiteVolume/lnInclude \ -I$(LIB_SRC)/meshTools/lnInclude #}; codeLibs #{ -lfiniteVolume \ -lmeshTools #};

请教一下各位大佬，我想要通过codedfixedvalue实现movingwallvelocity这个边界条件的效果，这个代码是我仿照movingwallvelocity边界条件弄得，我用自己写的这个边界条件算了涡激振动，算出来的结果和movingwallvelocity结果几乎没有差别。计算完成后，我将输出结果导入tecplot，可以正常导入，但是切片时会出现以下报错，目前没有找到问题出在哪儿，应该是我写的这个边界条件tecplot无法识别导致的，但我不知道问题出在哪儿。
641bf902-dad9-4f14-91ab-7f5bc7eda851-QQ20251011-140535.png

OverInterDyFoam计算放大尺度FloatingObject发散的问题

荣光 — Fri, 10 Oct 2025 14:27:37 GMT

在V2412中，我的原始case是计算自由液面上漂浮的物体，可以正常运行。但是一改变浮体的尺寸（对应的质量和惯性矩都改好了，水的密度和粘性不变）就会发散。

之后尝试了修改dynamicMesh里面的accelerationRelaxtion和accelerationDamp，发现修改后者到0.3以下才能正常运行算例
然而当阻尼这么大时已经不符合真实情况了，浮体的运动受影响了。

Github上有一个模组认为FloatingObject算例的added mass instability会造成这种发散：https://github.com/FloatStepper/FloatStepper但其不支持重叠网格。

目前打算从OpenFOAM本身解决该问题，但无从下手...

paraview后处理.vtp文件

XieXiaoyang — Fri, 10 Oct 2025 13:46:03 GMT

我有时候会用matlab处理批量vtp文件，用vtkRead。TransientU_x.png

如何实时获得每个y对应的某个量最大值所在的位置

Rachel0096 — Wed, 08 Oct 2025 13:17:12 GMT

大家晚上好
1136009f-4795-4f3b-8f6f-9fdd3b5900fa-image.png
图中红线对应的是thermicity peak，类似的我想实时获得每个y对应的热释率最大值所在的位置，想问一下有没有大佬有经验呀！

compressibleInterDyMFoam计算近壁空化泡出现奇怪的结果

火山口玩泥巴 — Sat, 04 Oct 2025 08:52:28 GMT

近期尝试使用compressibleInterDyMFoam计算近壁面空化泡的溃灭，但在计算一定时间后出现了一些与实际不符的现象。
算例设置如下：
初始气泡半径R0=0.1mm，距离壁面距离λ=2，原计算域大小为5mm5mm1mm，实际计算域取其1/4。草图如下：
1265c60d-3e45-457e-84df-081f40874abd-image.png
初始网格密度为Δx=Δy=Δz=0.025mm。后续用自适应网格进行加密，dynamicMeshDict：

FoamFile { version 2.0; format ascii; class dictionary; location "constant"; object dynamicMeshDict; } // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // dynamicFvMesh dynamicRefineFvMesh; //dynamicFvMesh staticFvMesh; dynamicRefineFvMeshCoeffs { // How often to refine refineInterval 1; // Field to be refinement on field alpha.water; // Refine field inbetween lower..upper lowerRefineLevel 0; upperRefineLevel 0.999; // If value < unrefineLevel unrefine unrefineLevel 10; // Have slower than 2:1 refinement nBufferLayers 3; // Refine cells only up to maxRefinement levels maxRefinement 5; // Stop refinement if maxCells reached maxCells 6000000; // Flux field and corresponding velocity field. Fluxes on changed // faces get recalculated by interpolating the velocity. Use 'none' // on surfaceScalarFields that do not need to be reinterpolated. correctFluxes ( (phi none) (nHatf none) (rhoPhi none) (alphaPhi none) (alphaPhiUn none) (alphaPhi10 none) (alphaPhi1Corr0 none) (ghf none) ); // Write the refinement level as a volScalarField dumpLevel true; }

计算到10.2μs后，气泡突然从溃灭状重新开始增大，而并没有完全溃灭，如图所示。
10.2μs：
a833b3e2-ebc9-43b4-8363-d153ec8db0e0-10.2um.png
10.4μs：
e248f2ea-78d2-4100-a3a9-551f88c499e3-10.4um.png

想问问大家这是什么原因所导致的，可能会是初始场设置不合理或是边界条件设置不合理？
初始场由于初始网格较疏以及后续加密的缘故，几个时间步后会变成这样：
170ce1d8-861d-4cfb-857b-2e2da4ed06eb-初始场.png

边界条件如下所示：
alpha.water:

boundaryField { wall { type zeroGradient; } outlet { type zeroGradient; } symmetry { type symmetry; } }

boundaryField { wall { type calculated; value uniform 100000; } outlet { type calculated; value uniform 100000; } symmetry { type symmetry; } }

p_rgh:

boundaryField { wall { type fixedFluxPressure; gradient uniform 0; value uniform 100000; } outlet { type fixedFluxPressure; gradient uniform 0; value uniform 100000; } symmetry { type symmetry; } }

boundaryField { wall { type zeroGradient; } outlet { type zeroGradient; } symmetry { type symmetry; } }

T.air:

boundaryField { wall { type zeroGradient; } outlet { type zeroGradient; } symmetry { type symmetry; } }

T.water:

boundaryField { wall { type zeroGradient; } outlet { type zeroGradient; } symmetry { type symmetry; } }

boundaryField { wall { type noSlip; } outlet { type zeroGradient; value uniform (0 0 0); } symmetry { type symmetry; } }

请教大佬们关于后处理中的阻力计算

Noob — Fri, 03 Oct 2025 10:19:00 GMT

@xpqiu 感谢感谢

求助centos安装openfoam

Amadeus — Thu, 02 Oct 2025 06:10:42 GMT

@旺财R CentOS 装2206是可以的, 但是wave2foam之类的就没试过了

网格y+求助

Miraitowa — Thu, 02 Oct 2025 05:20:10 GMT

@coolhhh 是空的计算域，几何部分已经挖掉了

带薄层的非均匀温度边界条件如何植入？

刘文凯 — Mon, 29 Sep 2025 09:58:36 GMT

非均匀温度边界条件可以通过codedFixedValue轻松实现，
带壁厚的温度边界条件可以通过externalWallHeatFluxTemperture实现。
如代码所示：

\\非均匀温度边界条件 wall { type codedFixedValue; name dummy; code #{ const vectorField& Cf = patch().Cf(); const scalar a0 = 0.1; forAll(Cf, faceI) { const scalar y = Cf[faceI].y(); (*this)[faceI] = a0*y; } #}; } \\薄壁温度条件 wall { type externalWallHeatFluxTemperature; mode coefficient; h uniform 1e10; Ta 1620; thicknessLayers (0.015); kappaLayers (4); kappaMethod fluidThermo; value $internalField; }

但是，当我想像下面这样编程使Ta为非均匀值时，会报错，即使我使用了codeInclude codeOptions codeLibs之后，仍会报错 ‘patch’ was not declared in this scope，似乎在externalWallHeatFluxTemperature中，无法访问到网格信息。请问各位大佬有解决办法吗？或者有没有别的办法实现带薄壁的非均匀温度边界呢？

‘vectorField’ does not name a type ‘Cf’ was not declared in this scope Ta coded; code #{ const vectorField& Cf = patch().Cf(); const scalar a0 = 0.1; forAll(Cf, faceI) { const scalar y = Cf[faceI].y(); (*this)[faceI] = a0*y; } #};

fluidsolver和isothermalfluid的继承关系问题

Miraitowa — Sun, 28 Sep 2025 02:39:00 GMT

各位老师好，之前我都在of10中修改求解器，但最近换到了of13，我希望去更改buossinesq假设时的求解器代码，但是对于fluidsolver和isothermalfluid的继承感到困惑，在官网看到fluidsolver继承isothermalfluid，但在fluidsolver中我并未找到相关代码，因此想请教各位在模块化代码下要怎么去修改求解器，谢谢大家！

两相流自由液面处速度场异常

liujm — Fri, 26 Sep 2025 14:53:12 GMT

各位老师好，

最近在使用K-omega SST湍流模型处理两相流的时候发现自由液面处的速度异常，但是不影响波面，如图所示。请问，这种情况是什么引起的，会影响计算结果吗？要如何解决？

万分感谢！

662a54c8-2f33-4bff-ac92-c8c5059e396d-image.png

如何验证网格的二阶收敛精度

lwjetmann — Fri, 26 Sep 2025 03:52:16 GMT

simpleFoam求解的定常问题，用的空间离散格式都是二阶的，审稿人让验证网格的二阶收敛精度。不知有没有人用过paraview将稀网格数据插值到密网格上并算过残差E和收敛阶次p？另外对于很复杂的三维网格，网格的非均匀性对收敛阶次公式的计算又该如何考虑呢？
35fb0dd0-05d4-4fcf-b850-3b3e8d317042-image.png

be5eddf7-3dd1-41c7-b320-e9fa9531fffc-image.png

PIMPLE发散， PISO就发散

星星星星晴 — Thu, 25 Sep 2025 14:29:01 GMT

Hi，

我这现在在做欧拉拉格朗日模拟，2-way。然后在模拟中途就爆掉了。。检查了一下flow，发现整个flow都乱掉了，因为parcel的聚集，毕竟两者是相互影响的。

但是神奇的是，同事把PIMPLE改为了PISO就过了。
理解不了为什么。

flow.0068.png

flow.0069.png

flow.0070.png

flow.0071.png

flow.0072.png

flow.0073.png

下面两个图是根据log文件画的，pimple 的拟合造成了温度场和压力场的变化？理解不了？求解惑

WhatsApp Image 2025-09-25 at 12.08.19.jpeg

WhatsApp Image 2025-09-25 at 12.08.07.jpeg

triSurfaceMesh对象切割cell单元

Amadeus — Thu, 25 Sep 2025 13:24:40 GMT

是想用基于cut CELL的浸没边界法? 可以参考FOAM Extend 4.1 之后的做法, 找ibPatch

如何使用moveDynamicMesh初始化变形网格

chenboyao — Thu, 25 Sep 2025 08:38:38 GMT

@liujm 。。。。。。。。。，除了你画图的时候浮体的重心位置变化了，其他设置都没啥变化的。

rhoCentralFoam求解器有没有更新MRF

yingqing — Thu, 25 Sep 2025 05:30:14 GMT

rhoCetralFoam求解器可不可以增加MRF功能。目前，有没有哪个版本做出更新，希望知道的大佬能给点提示

Rheointerfoam求解不稳定

WeakForm — Thu, 25 Sep 2025 01:31:15 GMT

目前在用open foam9版本的rheotool求解粘弹性单气泡上升问题，
使用PPTLog模型或者FENE-CRLog模型稳定性都很差，Co数缩小到0.05依然发散，GiesekusLog模型能够收敛但是气液界面看起来很怪，不是光滑的水滴型。
a3619796-3abb-4e40-9a05-b1c5fb886877-image.png
仿真主要设置参考的官方case，即theta的散度离散采用高阶形式，
PPT模型和fene模型主要发散都发生在theta计算特征值或者计算应力时候发散，使用clipping办法限制相关量的范围也没有取得很好的效果
各位前辈有没有好的建议

Openfoam全系统在vscode下的版本切换问题

夏伯阳 — Wed, 24 Sep 2025 07:47:00 GMT

@星星星星晴在 Openfoam全系统在vscode下的版本切换问题中说：

把你的默认of改一下因为每开一个terminal 就要读~/.bashrc文件

感谢前辈指点！

OF13 源代码编译找不到scotch.h

星星星星晴 — Tue, 23 Sep 2025 11:54:11 GMT

ubuntu 24.04 OF13 编译找不到scotch.h
解决办法：https://bugs.openfoam.org/view.php?id=4270
安装bison

拉格朗日injection model的自编译

youhaoyu — Tue, 23 Sep 2025 02:52:05 GMT

@星星星星晴谢谢大佬回复，我后面也看看这个头文件如何使用

OpenFOAM

在ofv2312上编译ofv2406代码遇到问题

老版本OpenFOAM E-L多cloud的解决方式

关于风场模拟入流ABL边界条件和其适配的网格大小

自定义边界条件遇到的问题

求助大佬们标准k-epsilon初始条件

codedfixedvalue实现运动边界速度的指定

OverInterDyFoam计算放大尺度FloatingObject发散的问题

paraview后处理.vtp文件

如何实时获得每个y对应的某个量最大值所在的位置

compressibleInterDyMFoam计算近壁空化泡出现奇怪的结果

请教大佬们关于后处理中的阻力计算

求助centos安装openfoam

网格y+求助

带薄层的非均匀温度边界条件如何植入？

fluidsolver和isothermalfluid的继承关系问题

两相流自由液面处速度场异常

如何验证网格的二阶收敛精度

PIMPLE发散， PISO就发散

triSurfaceMesh对象切割cell单元

如何使用moveDynamicMesh初始化变形网格

rhoCentralFoam求解器有没有更新MRF

Rheointerfoam求解不稳定

Openfoam全系统在vscode下的版本切换问题

OF13 源代码编译 找不到scotch.h

拉格朗日injection model的自编译

OF13 源代码编译找不到scotch.h