Calculating fvm::laplacian(tensor,scalar) on cylinder
 

I am using OpenFoam 2.2.2. I'm getting the following error when trying to solve laplacianFoam for a cylinder using a slight modification in which laplacianFoam takes a tensor as the coefficient D.

#0 Foam::error::printStack(Foam::Ostream&) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 void Foam::fvc::surfaceIntegrate<double>(Foam::Field<do uble>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#4 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::surfaceIntegrate<double>(Foam::Geometri cField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#5 Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::fvc::div<double>(Foam::GeometricField<double , Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#6 Foam::fv::gaussLaplacianScheme<double, Foam::Tensor<double> >::fvmLaplacian(Foam::GeometricField<Foam::Tensor< double>, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#7 Foam::fv::laplacianScheme<double, Foam::Tensor<double> >::fvmLaplacian(Foam::GeometricField<Foam::Tensor< double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/lib/libfiniteVolume.so"
#8
in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/bin/laplacianFoam"
#9 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#10
in "/home/jim/OpenFOAM/OpenFOAM-2.2.2/platforms/linux64GccDPOpt/bin/laplacianFoam"
Floating point exception (core dumped)

The problem appears to be with the mesh that I am using because the same solver will run just fine for the provided flange example. The solver also works on the cylinder mesh if a scalar is used for D.

I have attached both the model files in cylinder.zip and the slightly modified solver so that it uses a tensor coefficient in laplacianFoamT.zip.

Any insight into how I can get this to work would be much appreciated.

Thanks,
Jim

It appears that the original mesh that I used with blockMesh was not suitable. I have remeshed the cylinder using snappyHexMesh and the swiftSnap tool developed for blender (it works really well!)

The revised case is attached below.

-James