[Six Chen]

Hello everyone


I am a beginner of OpenFoam, when I try to run the simulation of a Axisymmetric Near-Sonic Jet from NASA (LINK: https://turbmodels.larc.nasa.gov/jetnearsonic_val.html). The simulation always stops right after it starts:


--> FOAM Warning :
From function virtual void Foam::wedgePolyPatch::calcGeometry(Foam::PstreamBu ffers&)
in file meshes/polyMesh/polyPatches/constraint/wedge/wedgePolyPatch.C at line 70
Wedge patch 'wedgeBack' is not planar.
At local face at (-0.150876 0.108539 -0.000947207) the normal (0 -0.00872654 -0.999962) differs from the average normal (-5.51924e-09 -0.00872661 -0.999962) by 5.84762e-15
Either correct the patch or split it into planar parts

SIMPLE: Convergence criteria found
p: tolerance 0.0001
U: tolerance 0.0001
"(k|epsilon|omega)": tolerance 0.001

Reading thermophysical properties

Selecting thermodynamics package
{
type hePsiThermo;
mixture pureMixture;
transport sutherland;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}

Reading field U

Reading/calculating face flux field phi

Creating turbulence model

Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
RAS
{
model kEpsilon;
turbulence on;
printCoeffs on;
Cmu 0.09;
C1 1.44;
C2 1.92;
C3 0;
sigmak 1;
sigmaEps 1.3;
}

Creating thermophysical transport model

Selecting thermophysical transport type RAS
Selecting default RAS thermophysical transport model unityLewisEddyDiffusivity
No MRF models present

No fvModels present
Creating fvConstraints from "system/fvConstraints"

Selecting finite volume constraint type limitPressure
Name: limitp
min 10132.5
max 358275


Starting time loop

Time = 0.01s

DILUPBiCGStab: Solving for Ux, Initial residual = 1, Final residual = 0.0877972, No Iterations 26
DILUPBiCGStab: Solving for Uy, Initial residual = 1, Final residual = 0.0445231, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 1, Final residual = 0.090669, No Iterations 3
DILUPBiCGStab: Solving for e, Initial residual = 0.562846, Final residual = 0.0444755, No Iterations 3
GAMG: Solving for p, Initial residual = 1, Final residual = 0.0827434, No Iterations 4
GAMG: Solving for p, Initial residual = 1.13893e-05, Final residual = 5.80668e-07, No Iterations 3
time step continuity errors : sum local = 0.0580252, global = -0.0253458, cumulative = -0.0253458
DILUPBiCGStab: Solving for epsilon, Initial residual = 0.854138, Final residual = 0.00951507, No Iterations 17
bounding epsilon, min: -7940.89 max: 2.52748e+08 average: 3.93132e+07
DILUPBiCGStab: Solving for k, Initial residual = 1, Final residual = 0.0640378, No Iterations 2
bounding k, min: -10.2461 max: 2.16176e+07 average: 411444
ExecutionTime = 1.04329 s ClockTime = 4 s

Time = 0.02s

DILUPBiCGStab: Solving for Ux, Initial residual = 0.995904, Final residual = 0.00997181, No Iterations 1
DILUPBiCGStab: Solving for Uy, Initial residual = 0.999948, Final residual = 0.0612613, No Iterations 1
DILUPBiCGStab: Solving for Uz, Initial residual = 0.894813, Final residual = 0.0781122, No Iterations 1
DILUPBiCGStab: Solving for e, Initial residual = 0.820998, Final residual = 1.8519, No Iterations 1000
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3 ? in "/lib/x86_64-linux-gnu/libm.so.6"
#4 Foam::hePsiThermo<Foam:siThermo::composite, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::calculate() at ??:?
#5 virtual thunk to Foam::hePsiThermo<Foam:siThermo::composite, Foam:ureMixture<Foam::sutherlandTransport<Foam:: species::thermo<Foam::hConstThermo<Foam:erfectGa s<Foam::specie> >, Foam::sensibleInternalEnergy> > > >::correct() at ??:?
#6 ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"
#7 ? in "/lib/x86_64-linux-gnu/libc.so.6"
#8 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#9 ? in "/opt/openfoam10/platforms/linux64GccDPInt32Opt/bin/rhoSimpleFoam"


It would be great if someone could tell me what exactly the problem is and how should I make the change. Thank you so much.

[dlahaye]

Seems that temperature reaches unrealistic high values (judging from error message).

Suggest to:

1/ limit temperature. Search this forum on how to;

2/ replace rhoSimpleFoam (stationary) with rhoPimpleFoam (transient). Check whether using smaller time steps works out;

3/ use finer mesh;

Good luck.