BUG: turbulentDFSEMInlet; Option = planarInterpolation; decomposePar causes Floating Point Exception for Fields R/L/U
Summary
When turbulentDFSEMInlet
with planarInterpolation
option is ON, decomposePar
causes FPE with the following error message:
Time = 0
Turbulent DFSEM patch inlet: interpolating field R from "/home/snoopy2/kuta/OpenFOAM/kuta-v1806/run/channel395DFSEM/constant/boundaryData/inlet/0/R"
#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in /lib64/libc.so.6
#3 Foam::pointToPointPlanarInterpolation::calcCoordinateSystem(Foam::Field<Foam::Vector<double> > const&) const at ??:?
#4 Foam::pointToPointPlanarInterpolation::pointToPointPlanarInterpolation(Foam::Field<Foam::Vector<double> > const&, Foam::Field<Foam::Vector<double> > const&, double, bool) at ??:?
#5 Foam::turbulentDFSEMInletFvPatchVectorField::patchMapper() const at ??:?
#6 Foam::tmp<Foam::Field<Foam::SymmTensor<double> > > Foam::turbulentDFSEMInletFvPatchVectorField::interpolateBoundaryData<Foam::SymmTensor<double> >(Foam::word const&) const at ??:?
#7 Foam::turbulentDFSEMInletFvPatchVectorField::turbulentDFSEMInletFvPatchVectorField(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#8 Foam::fvPatchField<Foam::Vector<double> >::adddictionaryConstructorToTable<Foam::turbulentDFSEMInletFvPatchVectorField>::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#9 Foam::fvPatchField<Foam::Vector<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#10 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::Boundary::readField(Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) at ??:?
#11 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) at ??:?
#12 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields() at ??:?
#13 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&, bool) at ??:?
#14 ? at ??:?
#15 ? at ??:?
#16 __libc_start_main in /lib64/libc.so.6
#17 ? at /home/abuild/rpmbuild/BUILD/glibc-2.18/csu/../sysdeps/x86_64/start.S:125
Floating point exception
Steps to reproduce
Please consider tutorials/incompressible/pimpleFoam/LES/channel395DFSEM
case.
Change the following under 0/U
:
inlet
{
type turbulentDFSEMInlet;
...
mapMethod nearestCell;
...
}
to
inlet
{
type turbulentDFSEMInlet;
...
mapMethod planarInterpolation;
...
}
Then execute the following respectively:
blockMesh
decomposePar
PS: scotch
and hierarchial
were tested with various numbers of numberOfSubdomains
.