Recently Active Topics

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

bestucan — Fri, 17 Oct 2025 06:27:17 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是相互独立的，求解的受力也都是独立的。

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

xpqiu — Wed, 15 Oct 2025 12:31:18 GMT

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

paraview后处理.vtp文件

XieXiaoyang — Wed, 15 Oct 2025 12:27:05 GMT

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

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

HITSC30 — Wed, 15 Oct 2025 00:50:58 GMT

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

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

Joker — Mon, 13 Oct 2025 15:02:38 GMT

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

求助centos安装openfoam

Amadeus — Mon, 13 Oct 2025 08:35:15 GMT

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

网格y+求助

Miraitowa — Mon, 13 Oct 2025 03:42:19 GMT

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

爆炸计算压力震荡

xuqiming — Mon, 13 Oct 2025 01:46:59 GMT

@好名字感谢回复，我试一下看看

fluent中的抑制回流

齐彦鹏 — Sun, 12 Oct 2025 07:52:31 GMT

各位大佬，我想请教一个问题，我现在在研究压气机出口存在反传压力波的情况，这道波强度比较大，需要进行瞬态模拟，那我在fluent中数值模拟的时候，出口的抑制回流应该打开还是关闭

几何模型

wu06040000 — Sun, 12 Oct 2025 01:39:52 GMT

怎么在斜面上画出竖直的几何微信图片_20251012093941.jpg

计算作用于cellZone上的压力和剪切力

liujm — Sat, 11 Oct 2025 15:40:37 GMT

@Linzd 多孔介质区域在toposet的时候本身只是选中了cell集合，只有在壁面上才能通过原始代码计算法向力和切向力。

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本身解决该问题，但无从下手...

Cd计算不准的问题-LES-Re3900-三维圆柱绕流-pisoFoam

chon01 — Fri, 10 Oct 2025 13:46:49 GMT

@chon01 NUMERICAL SIMULATIONS ON FLOW AROUND AN INCLINED CIRCULAR CYLINDER AT HIGH REYNOLDS NUMBER

triSurfaceMesh对象切割cell单元

Amadeus — Thu, 09 Oct 2025 12:04:55 GMT

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

关于wave2foam造波问题

chenboyao — Thu, 09 Oct 2025 07:57:06 GMT

@liujm 对的，通常我们计算整个风机都是采用komegaSST湍流模型，毕竟这个湍流模型基本上是做风机的都会用的，很少用其他湍流模型。

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

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

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

离散伴随和连续伴随的优劣

Noob — Wed, 08 Oct 2025 05:19:09 GMT

@木木枭我在做离散伴随，大佬有机会的话交流哇

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

Noob — Tue, 07 Oct 2025 03:23:42 GMT

@xpqiu 感谢感谢