CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Problem when starting rhoSimpleFoam

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 19, 2022, 21:51
Default Problem when starting rhoSimpleFoam
  #1
New Member
 
Six Chen
Join Date: Sep 2022
Posts: 1
Rep Power: 0
Six Chen is on a distinguished road
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.
Six Chen is offline   Reply With Quote

Old   September 20, 2022, 06:53
Default
  #2
Senior Member
 
Domenico Lahaye
Join Date: Dec 2013
Posts: 517
Blog Entries: 1
Rep Power: 14
dlahaye is on a distinguished road
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.
dlahaye is offline   Reply With Quote

Reply

Tags
error, rhosimplefoam, rhosimplefoam error

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem with rhoSimpleFoam - injectorPipe Rebeca OpenFOAM Pre-Processing 0 April 22, 2021 05:14
rhoSimpleFoam convergence problem - bounding omega inf.vish OpenFOAM Running, Solving & CFD 1 October 20, 2020 09:20
[GAMBIT] HELP!I have a problem when starting the unsteady iteration usst139 ANSYS Meshing & Geometry 3 January 21, 2010 05:04
Problem starting Gambit ViHAR MALViYA FLUENT 0 November 10, 2006 11:26
problem with starting ansyslmd Mike CFX 0 March 2, 2005 08:59


All times are GMT -4. The time now is 16:58.