CFD中文网

郑

@youhaoyu 谢谢，我找到这个文件了，但老哥你是怎么编译makeParcelForces.H的呢？直接编译intermediate里的make嘛？

郑

我对拉格朗日粒子的受力模型中（/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces）添加了一个自己写的受力项（unidirectional），基本内容都是根据已有的模型修改的公式得来的。编译已经成功了，并且在$FOAM_USER_LIBBIN中成功存在libInterfaceForceUnidirectional.so，但是实际计算的时候ParticleForce类中并没有这一项unidirectional。chatgpt各种代码排查都通过了但没啥用，还请指点

郑

@coolhhh 目前是没有出现过类似的报错，感觉应该是起效了。具体方法网上很多你搜一下

郑

@bestucan 十分感谢！

郑

@bestucan 谢谢大佬，所以您的建议是，我们把虚拟内存swap拉高到256，然后buff/cache没有必要清除或者长周期清除一次，是这样嘛？

郑

我们组的硬件是个9节点的机架式服务器，平常算fluent、OpenFOAM都有。一般在一天的计算后，每个节点都会积攒不定大小的buff/cache缓存，如图

本着能跑就别动它的原则，一开始没有管，继续用，后来很多用fluent的人反映算到一半算例炸了，就开始怀疑是缓存太高，留给计算用的不够。于是我自己在每个节点上，写了个定时任务，每天凌晨自动清理缓存，如图

buffccache_clean.sh文件：

不知道这种处理办法是否合理，还请各位交流发表下意见。目前隔三差五还是有同学反应算例内存不够，我怀疑还是一天内就把缓存积累满了。

郑

@李东岳因为我做的是含颗粒液滴的模拟，只是粒子目前还比较正常，就没显示。这个算例倒是没有发散，但是这个速度场最高能到20m/s左右，所以一直有疑惑，并且和文献中的实验相对比，确实会发现液体铺展太快了，您看是不是和接触角有关？

郑

我在做一个下落液滴撞击立柱的case，但是液滴4m/s接触立柱时，会出现一片极大的速度场，最终导致流体铺展速度极快。

跑的是层流，用的MPPICInterFoam，时间步长1e-6

U
type            noSlip;

p_rgh 
type            fixedFluxPressure;
value           uniform 0;

alpha.water
type           constantAlphaContactAngle;
theta0         70;
limit          gradient;
value          uniform 0;

希望大腿们指点！

郑

@Shihang-Chen 在多相流液体铺展行为中说：

您好，我想请教一下，您在原本的矩形网格上是否有尝试过削弱这种铺展的行为呢？

我后来也尝试把固体的压力条件从0gradient改成了fixedFluxPressure，这个更适合密度发生变化的多相流模拟。不过你说的原本矩形网格削弱铺展，往往伴随着液滴接触固体时，形成一个较高的速度场，进而造成了算例发散，还没有得到很好解决。如有兴趣咱们可以深入交流

郑

自己回答一下哈，这个问题应该是出在网格加密方式上，按照文献上扁平的瓦状网格进行加密，最后得到了很好的效果
682413cb-e876-42a1-9639-6bff19bb514e-微信截图_20240611212200.png

郑

我在做单根不锈钢丝切割液滴的模拟，发现液滴触碰到固体后，会快速铺展覆盖固体，这明显与实际情况不符：

求解器用的v2212的MPPICInterFOAM,固体壁面的边界条件为：U noSlip；P fixedFluxPressure；Alpha constantAlphaContactAngle，接触角为155疏水界面

文献和实际可视化拍摄中，液体不会有如此快速的铺展，还请各位指点一二

郑

eb812a45-a24b-435b-a01c-246db6db8afc-250c0b40b51035fb0e43108730efee4.jpg
请教各位：
我想如图设置一个扇形的计算域，中间一段圆环设置为旋转动网格，其中涉及a-a' 与b-b'两个对应面之间设置的问题，因为整个过程中，两组对应面并不是时刻交互的，某些面区域会作为边界，这种情况应该如何设置呢？谢谢

郑

@zhe 谢谢，我按照方法排查尝试了一下，发现并不是粒子的加入影响了液体速度场，二而是MPPICInterFoam对于setFields对于输入速度场就有问题。下图一个是MPPICInterFoam计算的，一个是interFoam计算的，同样设置的5m/s液体下落速度，但结果截然不同。@李东岳也请麻烦李老师能否帮忙看一下，谢谢！

郑

@zhe 您好我想请教个问题，就是我MPPICInteFoam设置液体和粒子初速度都不管用。我是setfields设置的液体初速度，kinematicCloudProperties-injectionModels-kinematicCloudPosition设置的粒子位置，kinematicCloudProperties-injectionModels-U0设置的粒子初速度。计算时发现液体和粒子都不运动，请教！

郑

@李东岳在 interFoam里面的VOF，液膜，CSF模型源于哪些资料？中说：

@dyj19901127
跟Henry直接要的。here：http://pan.baidu.com/s/1cLrAkm

@李东岳在 interFoam里面的VOF，液膜，CSF模型源于哪些资料？中说：

@dyj19901127
跟Henry直接要的。here：http://pan.baidu.com/s/1cLrAkm

@李东岳您好李老师，这个分享链接过期了，可否方便再发一个？谢谢

郑

版本：v2206
soler：MPPICInterFoam
问题：非0时间无法执行decomposePar，提示时间文件夹下没有polymesh/point文件，无法进行网格划分
场景：之前我已经从0并行计算到0.02，并且合并了网格，但是由于存储问题删除了之前的processor们，现在想由之前合并好的0.02重新划分出processor，好像有些问题，请教各位大神！

郑

本人OpenFOAM入门小白，通过自学达到使用一般求解器的水平，对简单的网格、相场生成、求解能够独立完成。
提出这个问题是这样一个契机，我个人课题通过MPPICInterFoam这个求解器能满足，但具体在拉格朗日粒子设置的kinematicCloudProperties文件设置中，关于粒子摄入方式所需的cellZoneInjection方式不甚了解，虽然已经定位了OF源码中的CellZoneInjection.C文件，但阅读起来较为吃力，难以理解。
请教各位大佬，当初在学习阶段，都是需要哪些只是储备与能力，才能实现“查源码就能看懂”这项技能的呢？谢谢

郑

@李东岳试了一下，文件名和“uniform-time”文件里面的时间信息也都更改了，分网格出来仍然没有“point”文件

郑

@李东岳李老师您好
这是中间时间文件decomposePar的结果

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 10
Exec   : decomposePar
Date   : Mar 21 2023
Time   : 16:04:48
Host   : "dyfluid-virtual-machine"
PID    : 3798
I/O    : uncollated
Case   : /home/dyfluid/OpenFOAM/dyfluid-10/run/2Dtry
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time



Decomposing mesh region0

Create mesh

Calculating distribution of cells
Selecting decomposer scotch

Finished decomposition in 0.038637 s

Calculating original mesh data

Distributing cells to processors

Distributing faces to processors

Distributing points to processors

Constructing processor meshes

Processor 0
    Number of cells = 42000
    Number of faces shared with processor 1 = 80
    Number of processor patches = 1
    Number of processor faces = 80
    Number of boundary faces = 85660

Processor 1
    Number of cells = 42000
    Number of faces shared with processor 0 = 80
    Number of processor patches = 1
    Number of processor faces = 80
    Number of boundary faces = 85180

Number of processor faces = 80
Max number of cells = 42000 (0% above average 42000)
Max number of processor patches = 1 (0% above average 1)
Max number of faces between processors = 80 (0% above average 80)

Time = 0.0005s

Processor 0: field transfer
Detected additional non-decomposed files in "/home/dyfluid/OpenFOAM/dyfluid-10/run/2Dtry/0.0005/uniform"
Processor 1: field transfer
Detected additional non-decomposed files in "/home/dyfluid/OpenFOAM/dyfluid-10/run/2Dtry/0.0005/uniform"

End

这个0.0005时间文件，是之前并行计算后reconstructPar出来的
我自己对比了一下，用0文件分网格出来，和0.0005分网格出来，确实后者会少一个“point”文件
个人推测和动网格有关，因为我用dambreak的并行计算对比过，dambreak中途进行并行计算是可行的，但是具体哪里出问题了还不清楚

郑

求解器interFoam，涉及旋转动网格
具体问题是：
从中间时间步（非0时间文件）开始计算的话，decomposePar生成的processor文件的constan-polyMesh文件中不包含“point”文件，导致后续并行计算报错

/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  10
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 10
Exec   : interFoam -parallel
Date   : Mar 20 2023
Time   : 14:31:00
Host   : "dyfluid-virtual-machine"
PID    : 5576
I/O    : uncollated
Case   : /home/dyfluid/OpenFOAM/dyfluid-10/run/2Dtry
nProcs : 2
Slaves : 1("dyfluid-virtual-machine.5577")
Pstream initialised with:
    floatTransfer      : 0
    nProcsSimpleSum    : 0
    commsType          : nonBlocking
    polling iterations : 0
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0.2

Selecting fvMeshMover motionSolver
Selecting motion solver: solidBody
[0] 
[0] 
[0] --> FOAM FATAL ERROR: [1] 
[1] 
[1] --> FOAM FATAL ERROR: 
[1] cannot find file "/home/dyfluid/OpenFOAM/dyfluid-10/run/2Dtry/processor1/constant/polyMesh/points"
[1] 
[1]     From function virtual Foam::autoPtr<Foam::ISstream> Foam::fileOperations::uncollatedFileOperation::readStream(Foam::regIOobject&, const Foam::fileName&, const Foam::word&, bool) const

[0] cannot find file "/home/dyfluid/OpenFOAM/dyfluid-10/run/2Dtry/processor0/constant/polyMesh/points"
[0] 
[0]     From function virtual Foam::autoPtr<Foam::ISstream> Foam::fileOperations::uncollatedFileOperation::readStream(Foam::regIOobject&, const Foam::fileName&, const Foam::word&, bool) const
[0]     in file global/fileOperations/uncollatedFileOperation/uncollatedFileOperation.C at line 539.
[0] 
FOAM parallel run exiting
[0] 
[1]     in file global/fileOperations/uncollatedFileOperation/uncollatedFileOperation.C at line 539.
[1] 
FOAM parallel run exiting
[1] 
--------------------------------------------------------------------------
MPI_ABORT was invoked on rank 1 in communicator MPI COMMUNICATOR 3 SPLIT FROM 0
with errorcode 1.

NOTE: invoking MPI_ABORT causes Open MPI to kill all MPI processes.
You may or may not see output from other processes, depending on
exactly when Open MPI kills them.
--------------------------------------------------------------------------
[dyfluid-virtual-machine:05572] 1 more process has sent help message help-mpi-api.txt / mpi-abort
[dyfluid-virtual-machine:05572] Set MCA parameter "orte_base_help_aggregate" to 0 to see all help / error messages

目前尚不清楚该如何解决，求助

CFD中文网

郑学习

帖子