拉格朗日求解器合并块问题
-
拉格朗日DPMFoam求解器使用过程中,我并行结束后合并不了块,出现以下代码
Create time Reconstructing fields for mesh region0 Time = 0.5 Reconstructing FV fields Reconstructing volScalarField::DimensionedInternalFields kinematicCloud:UCoeff Reconstructing volVectorField::DimensionedInternalFields kinematicCloud:UTrans Reconstructing volScalarFields k.air mu.air alpha.air epsilon.air nut.air p rho.air Reconstructing volVectorFields U.air Reconstructing surfaceScalarFields phi.air Reconstructing point fields No point fields Reconstructing lagrangian fields for cloud kinematicCloud --> FOAM FATAL ERROR: cell, tetFace and tetPt search failure at position (0.855908 3.37988 -0.0348682) for requested cell 78838 If this is a restart or reconstruction/decomposition etc. it is likely that the write precision is not sufficient. Either increase 'writePrecision' or set 'writeFormat' to 'binary' From function void Foam::particle::initCellFacePt() in file /home/zly/OpenFOAM/OpenFOAM-3.0.0/src/lagrangian/basic/lnInclude/particleI.H at line 735. FOAM aborting #0 Foam::error::printStack(Foam::Ostream&) at ??:? #1 Foam::error::abort() at ??:? #2 Foam::reconstructLagrangianPositions(Foam::polyMesh const&, Foam::word const&, Foam::PtrList<Foam::fvMesh> const&, Foam::PtrList<Foam::IOList<int> > const&, Foam::PtrList<Foam::IOList<int> > const&) at ??:? #3 ? at ??:? #4 __libc_start_main in "/lib64/libc.so.6" #5 ? at ??:? Aborted (core dumped
然后当我想mpirun -np继续运行时候,他就运行不了了,出现以下报错
Create time Create mesh for time = 7 PIMPLE: Operating solver in PISO mode Reading g Reading field U Reading field p Reading/calculating continuous-phase face flux field phic Creating turbulence model Selecting incompressible transport model Newtonian Creating field alphac Constructing kinematicCloud kinematicCloud Constructing particle forces Selecting particle force ErgunWenYuDrag Selecting particle force gravity Constructing cloud functions none Constructing particle injection models Creating injector: model1 Selecting injection model patchInjection Constructing 3-D injection Choosing nParticle to be a fixed value, massTotal variable now does not determine anything. [6] [6] [6] --> FOAM FATAL ERROR: [6] cell, tetFace and tetPt search failure at position (0.856587 3.37816 -0.0300663) for requested cell 78838 If this is a restart or reconstruction/decomposition etc. it is likely that the write precision is not sufficient. Either increase 'writePrecision' or set 'writeFormat' to 'binary' [6] [6] From function void Foam::particle::initCellFacePt() [6] in file /home/zly/OpenFOAM/OpenFOAM-3.0.0/src/lagrangian/basic/lnInclude/particleI.H at line 735. [6] FOAM parallel run aborting [6] [6] #0 Foam::error::printStack(Foam::Ostream&) at ??:? [6] #1 Foam::error::abort() at ??:? [6] #2 ? at ??:? [6] #3 ? at ??:? [6] #4 ? at ??:? [6] #5 ? at ??:? [6] #6 ? at ??:? [6] #7 __libc_start_main in "/lib64/libc.so.6" [6] #8 ? at ??:? -------------------------------------------------------------------------- MPI_ABORT was invoked on rank 6 in communicator MPI_COMM_WORLD 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. -------------------------------------------------------------------------- ^Cmpirun: abort is already in progress...hit ctrl-c again to forcibly terminate
我认为主要报错原因是这一段话
cell, tetFace and tetPt search failure at position (0.856587 3.37816 -0.0300663) for requested cell 78838 If this is a restart or reconstruction/decomposition etc. it is likely that the write precision is not sufficient. Either increase 'writePrecision' or set 'writeFormat' to 'binary
但是我不理解他意思,是指的这个时间段颗粒在分块的交界线上,所以产生报错现象吗???我感觉of不会出现这种问题把。那出现这种问题有什么补办法么,难道只有修改decomposePar个数,或者把网格画粗一点重新计算么?有没有什么方法能够让他合并成功?
主要为我问题一个是我合并不了块,无法看拉格朗日项,第二个问题是我无法继续进行计算(因为我中途暂停了只是想看一下流场是什么样子的。我能合并流场的块,不能合并拉格朗日项的块很难受)。