CFD Online Discussion Forums - rhoCentralFoam - kOmegaSST

Hello everyone,

I don't really have a single question but I am faced with a simulation problem I haven't been able to solve, so I hope someone here can help me.

I am no CFD expert but for my master thesis in mechanical engineering I have to simulate then build a linear aerospike nozzle to accelerate high pressure air to super sonic speed.

For the simulation I created a mesh using gmsh, imported it to foam and I am using the rhoCentralFoam solver with the kOmega turbulence model.

The simulation starts "fine" but crashes after 6 microseconds... The error log I get is this:

Code:

#0  Foam::error::printStack(Foam::Ostream&) at ??:?

#1  Foam::sigFpe::sigHandler(int) at ??:?

#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"

#3  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?

#4  Foam::hePsiThermo<Foam::psiThermo, Foam::pureMixture<Foam::sutherlandTransport<Foam::species::thermo<Foam::hConstThermo<Foam::perfectGas<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?

#5  ? at ??:?

#6  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"

#7  ? at ??:?

Floating point exception (core dumped)

which looks like something is numerically "exploding" in the computation of the temperature (correct me if I am wrong).

Here is a link to the whole simulation case. I copy/pasted the error log from the terminal at the end of the "Log" file. The "spike26deg_mm.geo" and "spike26deg_mm.msh" are respectively the geometry and mesh file from GMSH. I set a write frequency of 1e-06 second because I wanted to see what was "happening" before the crash of the simulation.

https://www.dropbox.com/sh/qd2iz7ifd...GP8kYkcxa?dl=0

You should be able to run the sim just like I did and toy around with the initial and boundary conditions. The inlet pressure is pretty high (6 atmospheres) but that's not the problem since the same problem occurs with 2 atmospheres.

If you want any more details, feel free to ask. I will post any breakthrought I get.

Thanks!

DGW