dynamic adaptive mesh with load balance crashes

[wengzf20]

Hi,

I am using blastReactingFoam for gaseous detonation simulation. The AMR is used together with loadbalance. The setup follows the tutorial case building3DWorkshop. The simulation crashes when mapping vector fields in the load balance step. The log information is attached.

I do notice that the issue has been mentioned in the user guide of blastFoam,

"There are also known problems in the transferring of data for calculated vector and movingWall boundary conditions. The fields used within the standard blastFoam solver have addressed this, but if adaptiveFvMesh is used with a different solver, its stability with balancing is not guaranteed."

Could you give more details how the issue can be solved? Thank you for your attention.

Best,
Z.

Log information attached

Selected 8765 cells for refinement out of 1992252.
...
max(p): 223162.182871, min(p): 6670
max(T): 2932.29448685, min(T): 298
ExecutionTime = 476.15 s ClockTime = 481 s

Refined from 1992252 to 2018547 cells.
Selected 7318 split edges out of a possible 130851.
Unrefined from 2018547 to 1996593 cells.
Maximum imbalance = 208.403765815 %
Distributing the mesh ...
Mapping the fields ...
[23] #0 Foam::error::printStack(Foam::Ostream&)[145] #0 Foam::error::printStack(Foam::Ostream&) at ??:?
[23] #1 Foam::sigFpe::sigHandler(int) at ??:?
[23] #2 ? at ??:?
[145] #1 Foam::sigFpe::sigHandler(int) in "/lib64/libc.so.6"
[23] #3 void Foam::writeEntry<Foam::Vector >(Foam::Ostream&, Foam::Field<Foam::Vector > const&) at ??:?
[23] #4 at ??:?
[145] #2 ?Foam::calculatedFvPatchField<Foam::Vector >::write(Foam::Ostream&) const at ??:?
[23] #5 Foam::GeometricField<Foam::Vector, Foam::fvPatchField, Foam::volMesh>::Boundary::writeEntry(Foam::word const&, Foam::Ostream&) const in "/lib64/libc.so.6"
[145] #3 at ??:?
[23] #6 void Foam::writeEntry<Foam::Vector >(Foam::Ostream&, Foam::Field<Foam::Vector > const&) at ??:?
[145] #4 Foam::Ostream& Foam::operator<< <Foam::Vector, Foam::fvPatchField, Foam::volMesh>(Foam::Ostream&, Foam::GeometricField<Foam::Vector, Foam::fvPatchField, Foam::volMesh> const&)Foam::calculatedFvPatchField<Foam::Vector >::write(Foam::Ostream&) const at ??:?
[23] #7 void Foam::fvMeshDistribute::sendFields<Foam::GeometricField<Foam::Vector, Foam::fvPatchField, Foam::volMesh> >(int, Foam::ListFoam::word const&, Foam::fvMeshSubset const&, Foam::Ostream&) at ??:?
[145] #5 at ??:?
[23] #8 Foam::fvMeshDistribute::distribute(Foam::List const&)Foam::GeometricField<Foam::Vector, Foam::fvPatchField, Foam::volMesh>::Boundary::writeEntry(Foam::word const&, Foam::Ostream&) const at ??:?
[23] #9 Foam::adaptiveBlastFvMesh::balance() at ??:?
[145] #6 at ??:?
[23] #10 Foam::adaptiveBlastFvMesh::refine(bool)Foam::Ostream& Foam::operator<< <Foam::Vector, Foam::fvPatchField, Foam::volMesh>(Foam::Ostream&, Foam::GeometricField<Foam::Vector, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
[23] #11 at ??:?
[145] #7 void Foam::fvMeshDistribute::sendFields<Foam::GeometricField<Foam::Vector, Foam::fvPatchField, Foam::volMesh> >(int, Foam::ListFoam::word const&, Foam::fvMeshSubset const&, Foam::Ostream&)? at ??:?
[23] #12 __libc_start_main at ??:?
[145] #8 Foam::fvMeshDistribute::distribute(Foam::List const&) in "/lib64/libc.so.6"
[23] #13 at ??:?
[145] #9 Foam::adaptiveBlastFvMesh::balance()? at ??:?
[145] #10 Foam::adaptiveBlastFvMesh::refine(bool) at ??:?