[jnanabrao]

I am running a Mach 2 jet LES calculation with dynamicKEqn using sonicFoam. It runs on a coarse mesh (200k cells). But when I map it to a fine mesh (25mil cells) with mapFields and same boundary conditions and run it, I get a sigFpe error in calculation of Eddy diffusivity in the correctNut function.

Selecting finite volume options model type limitTemperature
    Source: limitT
    - selecting all cells
    - selected 25291152 cell(s) with volume 0.0341883194
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib64/libc.so.6"
#3  Foam::sqrt(Foam::Field<double>&, Foam::UList<double> const&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::sqrt<Foam::fvPatchField, Foam::volMesh>(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#5  Foam::LESModels::dynamicKEqn<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::Ck(Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const at ??:?
#6  Foam::LESModels::dynamicKEqn<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correctNut(Foam::GeometricField<Foam::SymmTensor<double>, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) at ??:?
#7  Foam::LESModels::dynamicKEqn<Foam::EddyDiffusivity<Foam::ThermalDiffusivity<Foam::CompressibleTurbulenceModel<Foam::fluidThermo> > > >::correctNut() at ??:?
#8  ? at ??:?
#9  __libc_start_main in "/lib64/libc.so.6"
#10  ? at ??:?

The error seems to originate in evaluation of Ck, particularly in evaluating the square root of field KK which is surprising as KK represents the subgrid kinetic energy and is obtained as :

volScalarField KK
    (
        0.5*(filter_(magSqr(this->U_)) - magSqr(filter_(this->U_)))
    );

My question is how do I resolve this? Ideally I would like to know if there is a reason why this works for the coarse mesh but not the fine one.
Any help is deeply appreciated.
JBR