[stack]

Hello all,

I using rhoCentralFoam to run a blunt body, at Mach 5.

The simulation outputs the error below after five timesteps, with a max Courant number of 0.046. This exact setup does not crash using sonicFoam, but I don’t think that solver is appropriate for this situation.

The boundary conditions on the inlet/outlets are waveTransmissive, for both p and U, and zeroGradient for T.
The BCs for the blunt are zeroGradient for p and T, and slip for U.

I have also attached pictures of the five timesteps, at the upper corner of the body, where I believe the discontinuity is causing it to crash.

If anyone could help me troubleshoot this problem, it would be greatly appreciated, thanks!

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/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::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
at ??:?
Floating point exception (core dumped)

[stack]

I got it working.

I changed fvSchemes, in the "interpolationSchemes" section from:

reconstruct(rho) vanLeer;
reconstruct(U) vanLeerV;
reconstruct(T) vanLeer;

to:

reconstruct(rho) upwind;
reconstruct(U) upwind;
reconstruct(T) upwind;

[yrganiri]

Hi,

I am having the same problem with rhoCentralFoam. Can I see your files (maybe that will help me)?

The thing is: I am running a simulation on a NACA airfoil with sonicFoam (no errors there), but when I change to rhoCentralFoam I get exactly the same error as you:

Mean and max Courant Numbers = 0.00232163 0.295612
deltaT = 3.89625e-05
Time = 0.000532113

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUx, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoUy, Initial residual = 0, Final residual = 0, No Iterations 0
diagonal: Solving for rhoE, Initial residual = 0, Final residual = 0, No Iterations 0
ExecutionTime = 0.78 s ClockTime = 1 s

#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/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::tmp<Foam::GeometricField<doub le, Foam::fvPatchField, Foam::volMesh> > const&) at ??:?
#5
at ??:?
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
at ??:?
Floating point exception

Actually, did you start from an OpenFoam tutorial?