CFD中文网

请问想在计算可压流动叶片的sonicLiquidFoam中添加Maxwell滑移边界与温度跳跃，应该如何操作呢

如下图，我想求这个通道中心线上的温度沿程分布，需要先求每个截面的面积平均温度，再取这条线的温度沿程平均值。这该怎么做到呢，查了一下现有的工具确实做不了呀，求大佬解答，谢谢。
ebab4892-66fb-4b7b-87ce-842c0381a798-image.png
下面是求面平均值，目前还只能求patch的平均，不能求沿程任意截面的

surfaceAve { type surfaceFieldValue; libs ("libfieldFunctionObjects.so"); writeControl writeTime; writeFields true; surfaceFormat raw; regionType patch;//(faceZone patch functionObjectSurface sampledSurface) name inlet; operation areaAverage; fields ( T ); }

大家好，变形网格求解cellDisplacement时，需要用到cellMotion边界条件。我在看其updateCoeffs()函数时（见下），不理解这句：this->operator[](i) = pp[i].average(points, pointMotion); average函数咋有两个输入参数了呢？

const fvPatch& p = this->patch(); const polyPatch& pp = p.patch(); const fvMesh& mesh = this->internalField().mesh(); const pointField& points = mesh.points(); word pfName = this->internalField().name(); pfName.replace("cell", "point"); const GeometricField<Type, pointPatchField, pointMesh>& pointMotion = this->db().objectRegistry::template lookupObject<GeometricField<Type, pointPatchField, pointMesh>> (pfName); forAll(p, i) { this->operator[](i) = pp[i].average(points, pointMotion); }

各位大佬有在foam-extend-4.0 植入 overset Mesh代码的经历吗？

ps：我知道fe41有自带的重叠网格，但是我个人用fe41支持的solids4Foam效果不是很好，因此打算用fe40支持的fsiFoam，但我又想要重叠网格功能。我在fe40的发布页看到fe40在库方面支持重叠网格。

Library-level support for Immersed Boundary Method and Overset Mesh
Release notes for foam-extend-4.0

因此，想咨询一下论坛大佬是否有相关经历，如果可以推荐GitHub上的库就更好了！ ：）

ps：我拿fe41的overset拷到fe40中，发现需要好多库以至于连src/foam/fields/FieldFields都没有。。。

OpenFOAM数字版里的sprayFoam的喷雾子模型里的LPT是分散相的同种性质的液滴在一个parcel，但是如果我想在工质雾化过程中里植入闪沸的过程，雾化过程属于多相流过程，那么气泡在液体的过程可以用高相分数模拟，气泡是离散型，那么离散型模型可以用朗格朗日粒子进行追踪。即里面有气泡生成，那么气泡成核、长大和破碎这个过程里的气泡也是可以包裹在朗格朗日框架下的parcel里面吗？
我看文献的气泡成核的气泡核的数量密度N与过热程度和喷射器配置有关如下是：
QQ20241003-222941.png
那么这个气泡parcel里气泡数量是怎么确定呢？可以随意给定还是用上述公式。
我想简化这个闪沸过程，每个液滴里面一个气泡，那么我就有一个问题，那么是哪个液滴里可以成核一个气泡，是一个液滴parcel里面生成一个气泡吗？那么这个气泡是用朗格朗日框架下的parcel得到的吗？那么可以在拉格朗日框架下计算液滴parcel和气泡parcel吗？通过我的理解是一个parcel里面包含了一个parcel，这样怎么个计算流程呢？可能我的描述不太清楚，也是因为我本人也是有点糊涂，有了解的大佬方便讲解一下吗？感谢感谢！

求大佬看看！
702a48d4-6fc4-4ebb-9b44-3ab220ebebbb-image.png
Added patches in = 0.02 s

Edge intersection testing:
Number of edges : 5008207
Number of edges to retest : 490832
Number of intersected edges : 475794
Decomposition method none [1]

Refinement phase

Found point (10 20 5) in cell 587820 on processor 0

Feature refinement iteration 0

Marked for refinement due to explicit features : 6 cells.
Determined cells to refine in = 9.79 s
Selected for feature refinement : 13 cells (out of 1506875)
Edge intersection testing:
Number of edges : 5008546
Number of edges to retest : 781
Number of intersected edges : 475794
Refined mesh in = 2.28 s
After refinement feature refinement iteration 0 : cells:1506966 faces:5008546 points:2028784
Cells per refinement level:
0 438931
1 42098
2 131721
3 456157
4 438059

Feature refinement iteration 1

Marked for refinement due to explicit features : 0 cells.
Determined cells to refine in = 9.22 s
Selected for feature refinement : 0 cells (out of 1506966)
Stopping refining since too few cells selected.

Surface refinement iteration 0

Marked for refinement due to surface intersection : 27528 cells.
Determined cells to refine in = 0.62 s
Selected for refinement : 29733 cells (out of 1506966)
hexRef8 : Dumping cell as obj to "/mnt/g/City4CFD/City4CFD/examples/TUDCampus/c4c_tst_20240924/cell_813430.obj"

--> FOAM FATAL ERROR: (openfoam-2312)
cell 813430 of level 2 uses more than 8 points of equal or lower level
Points so far:8(467622 506063 542474 660192 812094 865708 953356 1395965)

From Foam::labelListList Foam::hexRef8::setRefinement(const labelList&, Foam::polyTopoChange&) in file polyTopoChange/polyTopoChange/hexRef8/hexRef8.C at line 3604.

FOAM aborting

[stack trace] #1 Foam::error::simpleExit(int, bool) in ~/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libOpenFOAM.so
#2 Foam::hexRef8::setRefinement(Foam::List<int> const&, Foam::polyTopoChange&) in ~/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libdynamicMesh.so
#3 Foam::meshRefinement::refine(Foam::List<int> const&) in ~/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libsnappyHexMesh.so
#4 Foam::meshRefinement::balanceAndRefine(Foam::string const&, Foam::decompositionMethod&, Foam::fvMeshDistribute&, Foam::List<int> const&, double, int) in ~/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libsnappyHexMesh.so
#5 Foam::snappyRefineDriver::surfaceOnlyRefine(Foam::refinementParameters const&, int, int) in ~/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libsnappyHexMesh.so
#6 Foam::snappyRefineDriver::doRefine(Foam::dictionary const&, Foam::refinementParameters const&, Foam::snapParameters const&, bool, Foam::meshRefinement::FaceMergeType, Foam::dictionary const&) in ~/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/lib/libsnappyHexMesh.so
#7 ? in ~/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/bin/snappyHexMesh
#8 ? in /lib/x86_64-linux-gnu/libc.so.6
#9 __libc_start_main in /lib/x86_64-linux-gnu/libc.so.6
#10 ? in ~/OpenFOAM-v2312/platforms/linux64GccDPInt32Opt/bin/snappyHexMesh

Aborted
通过报错信息生成了个obj文件，求问这个是个啥，为什么会生成

OpenFOAM教程 初学者入门 里面提到的自带教程是在哪里呀？f54d69d6-02b5-4ec6-ba96-8400c60d6cd9-image.png

请问老师关于数据驱动的例子训练结束后生成的pth文件如何进一步调用？

大家好呀，李老师之前分享过【OpenFOAM-8 + OpenFOAM-v2106 + OpenFOAM-v2006 + OpenFOAM-6 + OpenFOAM-2.1.1 + foam-extend-4.1】虚拟机呀，但是现在东岳流体官网分享链接已经没了，哪位大佬下载过的，是否方便再分享一下呀，感谢！

bea87e1d-5977-463b-a33d-8a00f20c75e9-image.png
请教下各位老师，我想在流场内某一处加一个气流，来模仿旋翼产生的下洗流场，应该怎么实现呢？加一个速度就行，但是要持续存在。模型是二维的。

气液两相不可压缩，氢气在液相中饱和溶解，不考虑温度影响，通过压力的变化得出液体组分中溶解氢的析出。

我正在模拟三维下的波浪-海流-结构物相互作用，主要用olaFlow，偶尔也用基金会版本的interFoam。然后发现olaFlow的造波边界也好，OF6自带的造波边界也罢，在多面体网格中表现的都不太行，olaFlow的边界尤甚。而且我在网上搜索到的算例，要么是结构网格，要么是笛卡尔网格（切体网格），没怎么看到多面体或者四面体网格。

所以想问一下，是不是OF框架下的边界造波，对多面体网格的支持都不是太好？或者说，相对于横平竖直的结构网格和迪卡尔网格，在多面体网格上应用造波边界需要什么特殊的设置？

超算中的参数有节点数、CPU数、任务数，OpenFOAM自身有个划分网格的数量。怎么设置可以最大化计算效率呢？

各位老师好，我在学习重叠网格案例时遇到了下图中的情况，挖洞区域非常小，离核心区域很近，请问如何避免种情况？
Snipaste_2024-09-26_16-54-08.png Snipaste_2024-09-26_16-55-48.png

各位老师好，

我想实现弹簧固定的模型上再耦合给定的运动，例如固定刚度的刚体在特定的时间范围上指定一个位移。位移是离散的数据，也可以拟合出函数。
弹簧模型可以在sixDof里实现，通过定义pointDisplacement文件中的边界类型为如下

type sixDoFRigidBodyDisplacement;

可以实现弹簧的定义。
但是问题在于特定的时间范围上指定一个位移，我想到的方法也是在pointDisplacement文件中自定义边界条件实现

type codedFixedValue;

这样分开定义都可以实现，叠加在一起就矛盾了。各位老师有相关的经验吗，是否有其他的实现方案。

例如我想提取这个区域中间轴线上的温度值，这个怎么提取，我只知道利用singleGraph提取一条直线上的但是两个椭圆曲线怎么提取呢ca1761ae-450b-4fb2-9ee0-a749ea131f04-image.png

东岳老师，求一个of3.0.0的虚拟机，实验室用的这个版本，想在虚拟机上做点操作试验下:143:

各位前辈，晚上好！请问在reactingFoam中如何实现一维层流燃烧？chemkin有自带的一维火焰传播模型，那reactingFoam是否可以实现，如何修改边界条件呢？

微信截图_20240923172218.png 请教下，算到第二步就报浮点异常，是边界调节设置错了吗？

在使用sprayFoam的时候，原模型设置中燃油只喷射0.00125秒，但是我对比实验设置燃油喷射时间为0.00154秒。想要将flowRateProfile中的右列数据不变，左列数据扩展到0.00125，但是计算出来的结果相差很大。将实验数据中的VIRTUAL INJECTION RATE写入flowRateProfile中，计算结果也有很大差别 想要请教一下flowRateProfile这个列表中数值的意义和量纲，或者sprayfoam有什么原出处的论文或者实验吗？
2024-09-22 15-52-45屏幕截图.png