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

cavity case in rhoPimpleFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 6, 2012, 11:14
Default cavity case in rhoPimpleFoam
  #1
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,173
Rep Power: 16
immortality is on a distinguished road
at the cavity case in rhoPimpleFoam folder occurs something odd!
and that is when it is running residual of rho is zero in command shell from start to end and in paraView also the pressure and velocity fields are not like icoFoam test case at all.
why it is so?
can anyone think about and tell me?
thanks.
immortality is offline   Reply With Quote

Old   December 6, 2012, 13:42
Default
  #2
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,173
Rep Power: 16
immortality is on a distinguished road
this is a time step(it seems rho is not solved):

Courant Number mean: 0.1716988131 max: 0.5466230657
Time = 2

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.002584181588, Final residual = 2.30733021e-07, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 0.002904801618, Final residual = 1.334497056e-07, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 0.00263319834, Final residual = 6.165570932e-08, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.004894302097, Final residual = 3.77008551e-05, No Iterations 10
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.064078536e-07, global = -5.170722179e-11, cumulative = -1.229301626e-08
rho max/min : 1.158623339 1.158622586
DICPCG: Solving for p, Initial residual = 0.0003231577867, Final residual = 7.031648618e-07, No Iterations 14
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.856414427e-09, global = -5.098053007e-11, cumulative = -1.234399679e-08
rho max/min : 1.158623339 1.158622586
DILUPBiCG: Solving for omega, Initial residual = 0.0006058053485, Final residual = 2.019874798e-06, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.002108742738, Final residual = 5.150541624e-06, No Iterations 1
ExecutionTime = 1.46 s ClockTime = 2 s

End
immortality is offline   Reply With Quote

Old   December 11, 2012, 09:13
Default
  #3
Senior Member
 
Join Date: Nov 2012
Posts: 168
Rep Power: 3
hz283 is on a distinguished road
Hi,

I found that I have the same problem when I used the solve rhoPimpleFoam . the Initial residual, final residual and iterations are all zero. Do you have got the cause for this issue? Thank you very much.

best regards,
hz
hz283 is offline   Reply With Quote

Old   December 11, 2012, 14:34
Default
  #4
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,066
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
Atleast in two posts in this forum, i saw the answer of this question,
it solves equation in explicit manner, so there is no residual , new rho will be calculated based on previous time step value
(diagonal solver)
__________________
Training Course on OpenFOAM at (http://www.isme.ir/)
My Weblog (http://openfoam.blogfa.com/)
nimasam is offline   Reply With Quote

Old   December 14, 2012, 04:07
Default
  #5
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,173
Rep Power: 16
immortality is on a distinguished road
yes some one that has a lot of kind has answered before.totally here all are kindful together and dont reproach eachother for replying.however then why these residuals are shown at the results when they are all zero.as i know in rhoCentralFoam density should be calculated implicitly from continuity.correct if needs.
immortality is offline   Reply With Quote

Old   December 14, 2012, 04:15
Default
  #6
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,066
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
Quote:
.as i know in rhoCentralFoam density should be calculated implicitly from continuity
why ? from where?
Quote:
solve(fvm::ddt(rho) + fvc::div(phi));
from code, it is clear that it calculates rho explicitly
__________________
Training Course on OpenFOAM at (http://www.isme.ir/)
My Weblog (http://openfoam.blogfa.com/)
nimasam is offline   Reply With Quote

Old   December 14, 2012, 04:49
Default
  #7
Senior Member
 
Join Date: Nov 2012
Posts: 168
Rep Power: 3
hz283 is on a distinguished road
Quote:
Originally Posted by nimasam View Post
why ? from where?

from code, it is clear that it calculates rho explicitly
Hi,

thank you for your help. When the simulation starts, in the first step, where does the code get the initial values of rho and phi(=rho*U*A)? because I found that we did not need to initialize the density in the folder 0/. Thank you very much.
hz283 is offline   Reply With Quote

Old   December 14, 2012, 05:53
Default
  #8
Senior Member
 
Nima Sam
Join Date: Sep 2009
Location: Tehran, Iran
Posts: 1,066
Blog Entries: 1
Rep Power: 14
nimasam is on a distinguished road
Send a message via Yahoo to nimasam
look at creatFields.H,

rho calculate from another class thermo
Quote:
volScalarField rho
(
IOobject
(
"rho",
runTime.timeName(),
mesh,
IOobject::NO_READ,
IOobject::AUTO_WRITE
),
thermo.rho(),
rhoBoundaryTypes
);
and phi is :
Quote:
surfaceScalarField phi("phi", mesh.Sf() & fvc::interpolate(rhoU));
you need to read the code
__________________
Training Course on OpenFOAM at (http://www.isme.ir/)
My Weblog (http://openfoam.blogfa.com/)
nimasam is offline   Reply With Quote

Old   December 14, 2012, 09:55
Default
  #9
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,173
Rep Power: 16
immortality is on a distinguished road
i accept the reality of approach but my argument is on the truth means what should be or is more rational to be in existance.therefore in your opinion what will be difference between a density based scheme and pressure one when both calculate rho explicitly and at the other hand what is the necessity for displaying and printing these zero and zero numbers from starting to end.i know we can disable writing of those by changing solver but what intentions do the providers have got?does exist a rational explanation for this behaviour?
immortality is offline   Reply With Quote

Old   January 9, 2013, 06:22
Default
  #10
Senior Member
 
immortality's Avatar
 
Ehsan
Join Date: Oct 2012
Location: Iran
Posts: 2,173
Rep Power: 16
immortality is on a distinguished road
why sum of continuity cumulative and global are negative?
immortality is offline   Reply With Quote

Reply

Thread Tools
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 On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
bc's of a komegaSST case Zymon OpenFOAM 11 July 25, 2010 09:36
cavity in flat plate and drag prediction Far FLUENT 0 May 19, 2010 14:47
drag of flat plate with cavity Far FLUENT 0 May 18, 2010 04:57
boundary conditions for cavity case sudhar OpenFOAM Pre-Processing 0 March 11, 2010 02:14
ParaFOAM issue with cavity case mschoenberg OpenFOAM Paraview & paraFoam 3 November 21, 2008 03:57


All times are GMT -4. The time now is 19:57.