我在做单根不锈钢丝切割液滴的模拟,发现液滴触碰到固体后,会快速铺展覆盖固体,这明显与实际情况不符:
求解器用的v2212的MPPICInterFOAM,固体壁面的边界条件为:U noSlip;P fixedFluxPressure;Alpha constantAlphaContactAngle,接触角为155疏水界面
文献和实际可视化拍摄中,液体不会有如此快速的铺展,还请各位指点一二
我在做单根不锈钢丝切割液滴的模拟,发现液滴触碰到固体后,会快速铺展覆盖固体,这明显与实际情况不符:
求解器用的v2212的MPPICInterFOAM,固体壁面的边界条件为:U noSlip;P fixedFluxPressure;Alpha constantAlphaContactAngle,接触角为155疏水界面
文献和实际可视化拍摄中,液体不会有如此快速的铺展,还请各位指点一二
请教各位:
我想如图设置一个扇形的计算域,中间一段圆环设置为旋转动网格,其中涉及a-a' 与b-b'两个对应面之间设置的问题,因为整个过程中,两组对应面并不是时刻交互的,某些面区域会作为边界,这种情况应该如何设置呢?谢谢
@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
目前尚不清楚该如何解决,求助
@zhe 您好请问和李老师合作的文章发表了么?不知道文章名称是什么,希望能够参考,谢谢!
具体研究体系是气液固催化反应,催化剂粉末与液体作为拟均相混合,气液不互溶,固体只存在于液相,流动过程中夹杂催化剂的液体会被切碎成小液滴。本人小白看了论坛上很多帖子,目前感觉这个MPPICInterFoam最符合我的体系,不知道有没有用过的前辈可以指点一二
已解决,个人发现在controldict中的 “ writeformat ” 为binary,更改为ascii即可
平时启动paraview采用算例内新建casename.foam文件的方法,正常的cavity pitzdaliy dambreak运行均无问题,但自适应网格算例dambreakWithObstacle出现了不能打开的情况。
求助各位大神
@tens 在 OpenFOAM如何设置球形相域? 中说:
- 球形用sphereToCell
- 在openfoam下\etc\caseDicts\preProcessing\setFieldsDict可以查到
感谢!跑通了!
图为我师兄指导所发
我师兄使用的of6,我在of9上如上输入显示
“
keyword set is undefind in dictionary "boxToCell"
"
所以请问球形相场的关键词该是什么?
并且对于其他形态的相场,我应该去哪个文件目录能找到对应的输入方法呢?谢谢各位大神