[chenyr10]

I was trying to run a LES case in parallel, the boundary for U at inlet is 'decayingTurbulenceInflowGenerator' from libLEMOS-2.3.x (my openfoam version is 2.3.0).

When I run decomposePar, the error occured like:

Time = 0
#0 Foam::error:rintStack(Foam::Ostream&) in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 ? in "/lib64/libc.so.6"
#3 Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#4 Foam::sqrt(Foam::tmp<Foam::Field<double> > const&) in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#5 Foam::decayingTurbulenceFvPatchVectorField::decayi ngTurbulenceFvPatchVectorField(Foam::fvPatch const&, Foam:imensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libLEMOS-2.3.x.so"
#6 Foam::fvPatchField<Foam::Vector<double> >::adddictionaryConstructorToTable<Foam::decayingT urbulenceFvPatchVectorField>::New(Foam::fvPatch const&, Foam:imensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/lib/libLEMOS-2.3.x.so"
#7 Foam::fvPatchField<Foam::Vector<double> >::New(Foam::fvPatch const&, Foam:imensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/bin/decomposePar"
#8 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::readField( Foam:imensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/bin/decomposePar"
#9 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields(Foam::dictionary const&) in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/bin/decomposePar"
#10 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readFields() in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/bin/decomposePar"
#11 Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/bin/decomposePar"
#12 ? in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/bin/decomposePar"
#13 ? in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/bin/decomposePar"
#14 __libc_start_main in "/lib64/libc.so.6"
#15 ? in "/vol6/home/ccegroup/OpenFOAM/OpenFOAM-2.4.0/platforms/linux64GccDPOpt/bin/decomposePar"

I found that the problem is relevant to the RField in the boundary, if the RField is set to be a uniform value like 'uniform (10 0 0 10 0 10)' then decomposePar is ok. But if I set it to be a list of nonuniform symmTensor, then the error happens.

Could someone please tell me why the error happens and how to solve it?

[xuhuang]

Hi Yiran,

This is probably due to Lund transformation in the boundary condition. You have a negative number while using function sqrt(*). Make sure your nonuniform RField gives you proper values, or you can use limiters.

Cheers,
Xu