Error solving Energy Eqs - Compressible Solver
 

Dear FOAMers,

I am working on a compressible solver (density based) that I created using the pre-existing chtMultiRegionFoam solver. As stated in the title, I get (see hereunder) an error when solving for the following energy equation:
fvScalarMatrix EEqn
(
fvm::ddt(rho, he) + fvm::div(phi, he)
+ fvc::ddt(rho, K) + fvc::div(phi, K)
==
fvm::laplacian(kEff, T)
+ HS
);
HS is just a volumetric source term, and kEff is the effective conductivity.
he is the internal energy (I compute it using an external Fortran code)
while K is set as being: const volScalarField& K = 0.5 * (U & U);

From the log I can tell that the error most probably arises from the laplacian term, and a division by 0 occurs somewhere, but after investigating I cannot find out the source of the error. Thanks a lot for your help !

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2
at sigaction.c:?
#3 Foam::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam::operator/(Foam::UList<double> const&, Foam::tmp<Foam::Field<double> > const&) at ??:?
#5 Foam::compressible::customturbulentTemperatureCoup ledBaffleMixedFvPatchScalarField::updateCoeffs() at ??:?
#6 Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::updateCoef fs() at ??:?
#7 Foam::fvMatrix<double>::fvMatrix(Foam::GeometricFi eld<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::dimensionSet const&) at ??:?
#8 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacianUncorrected(Foam::GeometricFi eld<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#9 Foam::fv::gaussLaplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#10 Foam::fv::laplacianScheme<double, double>::fvmLaplacian(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#11 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) at ~/OpenFOAM/OpenFOAM-2.3.0/src/finiteVolume/lnInclude/fvmLaplacian.C:219 (discriminator 1)
#12 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::laplacian<double, double>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ~/OpenFOAM/OpenFOAM-2.3.0/src/finiteVolume/lnInclude/fvmLaplacian.C:251 (discriminator 5)
#13
at ~/OpenFOAM/OpenFOAM-2.3.0/applications/solvers/heatTransfer/CompressibleHeSolver/./fluid/EEqn.H:24
#14 __libc_start_main at ??:?
#15
at ??:?
Floating point exception (core dumped)

Hi,

After I initialized "kEff" to a default value I was able to get past the first error I described in the previous post. However, now I get the following error:

#0 Foam::error::printStack(Foam::Ostream&) at ??:?
#1 Foam::sigSegv::sigHandler(int) at ??:?
#2
at sigaction.c:?
#3 memcpy at interp.c:?
#4 std::string::append(std::string const&) in "/home/faabid/OpenFOAM/ThirdParty-2.3.0/platforms/linux64/gcc-4.8.1/lib64/libstdc++.so.6"
#5 std::basic_string<char, std::char_traits<char>, std::allocator<char> > std::operator+<char, std::char_traits<char>, std::allocator<char> >(std::basic_string<char, std::char_traits<char>, std::allocator<char> > const&, std::basic_string<char, std::char_traits<char>, std::allocator<char> > const&) at ~/OpenFOAM/ThirdParty-2.3.0/platforms/linux64/gcc-4.8.1/include/c++/4.8.1/bits/basic_string.h:2370
#6
at ~/OpenFOAM/OpenFOAM-2.3.0/src/finiteVolume/lnInclude/fvcDiv.C:240 (discriminator 2)
#7
at ~/OpenFOAM/OpenFOAM-2.3.0/applications/solvers/heatTransfer/CompressibleHeSolver/./fluid/EEqn.H:24 (discriminator 1)
#8 __libc_start_main at ??:?
#9
at ??:?
Segmentation fault (core dumped)

Where fluid/EEqn.H:24 corresponds to:
+ fvc::ddt(rho, K) + fvc::div(phi, K)
I also remind that K is set to: const volScalarField& K = 0.5*magSqr(U);

Could anyone please help me ??