[zhajingjing]

Hi , I met floating exception error when using interDyMFoam solver to run my cases.
Surprisingly, this error didn't happen until recently. Previously, I just found the gamma values are not calculated well. Please see http://www.cfd-online.com/Forums/openfoam/73494-problem-about-max-gamma-1-a.html. After I improve the
tolerance , that problem seems to be fixed. I run several cases well.
I remember that I make no change recently. So why this error pops up?

Please give me some hint about this problem. Thanks!

Foam::error:rintStack(Foam::Ostream&) in "/home/jingjing/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/home/jingjing/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 __restore_rt at sigaction.c:0
#3 void Foam::MULES::limiter<Foam:neField, Foam::zeroField, Foam::zeroField>(Foam::Field<double>&, Foam:neField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::zeroField const&, Foam::zeroField const&, double, double, int) in "/home/jingjing/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#4 void Foam::MULES::explicitSolve<Foam:neField, Foam::zeroField, Foam::zeroField>(Foam:neField const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, Foam::zeroField const&, Foam::zeroField const&, double, double) in "/home/jingjing/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#5 Foam::MULES::explicitSolve(Foam::GeometricField<do uble, Foam::fvPatchField, Foam::volMesh>&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh>&, double, double) in "/home/jingjing/OpenFOAM/OpenFOAM-1.5/lib/linux64GccDPOpt/libfiniteVolume.so"
#6 main in "/home/jingjing/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/interDyMFoam"
#7 __libc_start_main in "/lib64/libc.so.6"
#8 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/home/jingjing/OpenFOAM/OpenFOAM-1.5/applications/bin/linux64GccDPOpt/interDyMFoam"

[zhajingjing]

If I decrease the tolerance, my case continue to run, but the previous problem that gamma value exceed 1 still exist.
So how to get accuracy on base of stability? As I am making numerical converge test for the wave generated in my wave tank, a proper accurate result should be needed in order to analyse wave dissipation due to viscous effect, on the base of eliminating numerical dissipation and dispersion effect. But I have no experience on numerical analysis. Any advise for improvement will be appreciated.

my fvSchemes file:
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
default Euler;
}

gradSchemes
{
default Gauss linear;
grad(U) Gauss linear;
grad(gamma) Gauss linear;
}

divSchemes
{
div(rho*phi,U) Gauss limitedLinearV 1;
div(phi,gamma) Gauss vanLeer;
div(phirb,gamma) Gauss interfaceCompression;
}

laplacianSchemes
{
default Gauss linear corrected;
}

interpolationSchemes
{
default linear;
}

snGradSchemes
{
default corrected;
}

fluxRequired
{
default no;
pd;
pcorr;
gamma;
}


// ************************************************** *********************** //


my fvSolution file:

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
pcorr PCG
{
preconditioner DIC;
tolerance 1e-7;
relTol 0;
};
pd PCG
{
preconditioner DIC;
tolerance 1e-15;
relTol 0;
};
pdFinal PCG
{
preconditioner DIC;
tolerance 1e-15;
relTol 0;
};
U PBiCG
{
preconditioner DILU;
tolerance 1e-7;
relTol 0;
};
cellDisplacement PCG
{
preconditioner DIC;
tolerance 1e-7;
relTol 0;
};
}

PISO
{
momentumPredictor no;
nCorrectors 2;
nNonOrthogonalCorrectors 0;
nGammaCorr 1;
nGammaSubCycles 2;
cGamma 1;
pRefCell 0;
pRefValue 0;
}


// ************************************************** *********************** //