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

High number iteration of pressure equation

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 11, 2014, 09:12
Default High number iteration of pressure equation
  #1
mgg
New Member
 
Join Date: Nov 2012
Posts: 27
Rep Power: 6
mgg is on a distinguished road
Hallo everyone,

I am doing DNS of heated pipe flow (one phase) with OF 2.2.2. The solver is buoyantpimplefoam. The pipe diameter is 2mm and length of 30D. I use Inflow Generator as inlet U BC (zeroGradient for p_rgh), convective BC (fixedValue for p_rgh) as outlet U BC. On the wall is fixedValue for U and fixedFluxPressure for p_rgh.

Mesh for pipe is O type. Total mesh number is 16 Mio. (y1 plus 0.15 (y1 under 1e6m)), time step is about 2.5e-5. The simulations is running but I observed that solving implicit pressure equation p_rgh withg PCG solver costs extremly high iterations (about 3000). Even with GAMG, the number of iteration reduce, but the simulation time is not shorter. Can anyone tell me the reason? Thanks.

Code:
Courant Number mean: 0.139046 max: 0.382011
Time = 0.18235

diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG:  Solving for Ux, Initial residual = 0.00172453, Final residual = 1.09164e-05, No Iterations 1
DILUPBiCG:  Solving for Uy, Initial residual = 0.0120354, Final residual = 7.28735e-05, No Iterations 1
DILUPBiCG:  Solving for Uz, Initial residual = 0.0120533, Final residual = 7.37222e-05, No Iterations 1
DILUPBiCG:  Solving for h, Initial residual = 0.002042, Final residual = 1.96612e-05, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.277668, Final residual = 0.0276332, No Iterations 15
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.98488e-06, global = 5.99754e-08, cumulative = 0.000204476
DICPCG:  Solving for p_rgh, Initial residual = 0.159447, Final residual = 0.0159405, No Iterations 2306
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.15001e-06, global = 3.65048e-08, cumulative = 0.000204512
PIMPLE: iteration 2
DILUPBiCG:  Solving for Ux, Initial residual = 0.000840124, Final residual = 5.48588e-07, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.00545121, Final residual = 1.47463e-07, No Iterations 3
DILUPBiCG:  Solving for Uz, Initial residual = 0.00546338, Final residual = 1.82867e-07, No Iterations 3
DILUPBiCG:  Solving for h, Initial residual = 2.04831e-05, Final residual = 6.68786e-07, No Iterations 1
DICPCG:  Solving for p_rgh, Initial residual = 0.187843, Final residual = 0.0180666, No Iterations 28
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.87581e-06, global = -3.96611e-08, cumulative = 0.000204472
DICPCG:  Solving for p_rgh, Initial residual = 0.165657, Final residual = 9.80264e-07, No Iterations 3181
diagonal:  Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 3.85402e-10, global = -9.35153e-14, cumulative = 0.000204472
ExecutionTime = 48955.7 s  ClockTime = 49428 s
mgg is offline   Reply With Quote

Old   July 11, 2014, 16:24
Default
  #2
mgg
New Member
 
Join Date: Nov 2012
Posts: 27
Rep Power: 6
mgg is on a distinguished road
I also wanna to mention that, in my DNS pipe flow without heat transfer (Re5400, d=2mm,time step=2.5e-5s,L=5D,cyclic bc, pimpleFoam), the number of pressure iteration is still high,about 300 to 500 iterations.
mgg is offline   Reply With Quote

Old   October 8, 2014, 03:47
Default
  #3
Member
 
Sebastian W.
Join Date: Nov 2012
Location: Saxony, Germany
Posts: 42
Rep Power: 6
nero235 is on a distinguished road
Send a message via ICQ to nero235
What does checkMesh say about your mesh? Also, an O type mesh for a pipe is a really bad idea since you get cells with really small angles in the center of the pipe. For any circular geometry it is a standard to use a OH grid or butterfly topology.

Hope this helps! Sebastian
nero235 is offline   Reply With Quote

Old   June 8, 2016, 12:12
Default
  #4
New Member
 
Lennart
Join Date: Feb 2016
Posts: 21
Rep Power: 2
elmo555 is on a distinguished road
Although this is an old thread, I'd still be interested in an answer.
elmo555 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
decomposePar no field transfert Jeanp OpenFOAM Pre-Processing 0 June 20, 2014 05:59
decomposePar pointfield flying OpenFOAM Running, Solving & CFD 28 December 30, 2013 16:05
Mesh Refinement Luiz Eduardo Bittencourt Sampaio (Sampaio) OpenFOAM Mesh Utilities 41 January 17, 2013 03:43
The correction on pressure equation of SIMPLE algorithm in MRFSimpleFOAM solver renyun0511 OpenFOAM Running, Solving & CFD 0 November 10, 2010 02:47
Constant velocity of the material Sas CFX 15 July 13, 2010 08:56


All times are GMT -4. The time now is 03:40.