CFD Online Logo CFD Online URL
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

Problem with pressure equation in rhoSimpleFoam

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

LinkBack Thread Tools Search this Thread Display Modes
Old   July 10, 2017, 09:05
Default Problem with pressure equation in rhoSimpleFoam
New Member
Join Date: Dec 2016
Posts: 24
Rep Power: 6
dariodario32 is on a distinguished road
Hi all,
I'm trying to solve an external field around a 3D wing. The flow is at Mach = 0.78.
The checkMesh application seems to be ok:

Checking geometry...
Overall domain bounding box (-39.54701387 -0.02841032651 -40.63372186) (46.72535601 50.17495664 40.1780381)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (4.795358352e-17 -3.74681294e-17 1.258816993e-16) OK.
Max cell openness = 1.73712416e-15 OK.
Max aspect ratio = 115.5201724 OK.
Minimum face area = 1.338572813e-05. Maximum face area = 59.99552135. Face area magnitudes OK.
Min volume = 6.131690308e-07. Max volume = 49.41030061. Total volume = 174080.4162. Cell volumes OK.
Mesh non-orthogonality Max: 82.83259138 average: 21.70561315
*Number of severely non-orthogonal (> 70 degrees) faces: 101.
Non-orthogonality check OK.
<<Writing 101 non-orthogonal faces to set nonOrthoFaces
Face pyramids OK.
Max skewness = 1.967783537 OK.
Coupled point location match (average 0) OK.

Mesh OK.
by the way, when I try to solve the system, something strange happen with the p-equation: the residuals seems to be already converged at the first iteration (see image attached).

I tried to change the fvscheme and fvsolution, but nothing change. Here is what I'm using:

default steadyState;

default Gauss linear;

default none;

div(phi,U) bounded Gauss upwind;
div(((rho*nuEff)*dev2(T(grad(U))))) Gauss linear;
div(phi,e) bounded Gauss upwind;
div(phi,epsilon) bounded Gauss upwind;
div(phi,k) bounded Gauss upwind;
div(phid,p) bounded Gauss upwind;

div(phi,Ekp) bounded Gauss upwind;
div((phi|interpolate(rho)),p) bounded Gauss upwind;

default Gauss linear corrected;

default linear;

default corrected;
solver GAMG;
tolerance 1e-12;
relTol 0.1;
smoother GaussSeidel;
nCellsInCoarsestLevel 20;

solver GAMG;
tolerance 1e-08;
relTol 0.01;
smoother GaussSeidel;
nCellsInCoarsestLevel 20;

nNonOrthogonalCorrectors 0;
rhoMin 0.1;
rhoMax 2.0;
transonic yes;
consistent yes;

p 1e-5;
U 1e-5;
e 1e-5;

// possibly check turbulence fields
"(k|epsilon|omega)" 1e-4;

p 0.3;
rho 0.7;
T 0.7;
p 0.7;
U 0.7;
e 0.7;
k 0.7;
epsilon 0.7;
T 0.7;
Running the simulation for 100 steps all the others residuals go under 10e-5, while pressure residual is static. Obviously the solution of the pressure is wrong, in paraFoam I see a uniform field for all time steps, and a non uniform field for velocity.
Do you know what could be wrong?

Attached Images
File Type: png Screenshot from 2017-07-10 13-59-33.png (25.2 KB, 21 views)
Attached Files
File Type: txt log.rhoSimpleFoam.txt (9.7 KB, 4 views)
dariodario32 is offline   Reply With Quote


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
Low pressure de Laval simulation convergence problem heksel8i FLUENT 3 July 22, 2013 10:28
About theory:how to get the pressure equation when solving u-p simultaneously? lzw2003 Main CFD Forum 8 April 8, 2012 22:13
error message cuteapathy CFX 14 March 20, 2012 06:45
I have a problem about possion pressure equation rainbow66 Main CFD Forum 4 November 26, 2011 10:25
A fundamental problem about Pressure equation of the potentialFoam solver dbxmcf OpenFOAM Running, Solving & CFD 0 October 6, 2006 11:32

All times are GMT -4. The time now is 22:05.