[Dorian1504]

Hello,

I would like to modify the nacaAirfoil tutorial case for sonicFoam in a steady state case. Is it possible to directly calculate the steady state without computing the transient regime with sonicFoam?
I naively tried to change the ddtSchemes from Euler to steadyState, as I can see on rhoSimpleFoam. However, sonicFoam return this error :

Courant Number mean: 0.000174184 max: 0.948172
#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::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#4 Foam:perator/(Foam::UList<double> const&, Foam::UList<double> const&) at ??:?
#5 Foam::diagonalSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8 Foam::fvMatrix<double>::solve() at ??:?
#9 ? at ??:?
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11 ? at ??:?
Floating point exception (core dumped)

What should I do to transform a transient case into steady state case ?
Moreover, I use sonicFoam because I heard that rhoSimpleFoam doesn't provide good result for high velocity. Is it true, why ? I would like to simulate a mach 3 simulation.


Thank you for your help.