
[Sponsors] 
June 29, 2016, 07:20 
SIMPLE method small time step

#1 
New Member
DSR
Join Date: Sep 2014
Posts: 11
Rep Power: 8 
Hi, I have created a CFD code following the SIMPLE method. It is a colocated grid and uses RhieChow interpolation. The problem is, I need to take rediculously small time steps for it to be stable.
Has anyone ever had this issue? Fixes? Any response will be much appreciated! 

June 29, 2016, 09:18 

#2  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 5,688
Rep Power: 60 
Quote:


June 29, 2016, 09:45 

#3 
New Member
DSR
Join Date: Sep 2014
Posts: 11
Rep Power: 8 
First order implicit. I've tried with multiple Re, but it seems the time step is always limited by a CFL of 1.0. Which seems very low for an implicit solver.


June 29, 2016, 21:33 

#4  
Senior Member
Arjun
Join Date: Mar 2009
Location: Nurenberg, Germany
Posts: 980
Rep Power: 25 
Quote:
Your program has a bug and it could be anywhere. Instablity in cfd solver could creep in from anything from small to big. It is not possible to point out bug from your information. 

June 30, 2016, 07:30 

#5 
New Member
DSR
Join Date: Sep 2014
Posts: 11
Rep Power: 8 
But is it safe to say that there is definitely a bug? That I should be able to run with much higher time steps?


June 30, 2016, 08:05 

#6  
Senior Member
Alex
Join Date: Jun 2012
Location: Germany
Posts: 2,731
Rep Power: 38 
Quote:


June 30, 2016, 08:43 

#7 
New Member
DSR
Join Date: Sep 2014
Posts: 11
Rep Power: 8 
Could it be related to the rhie chow or pressure correction methods which both use the coefficient value from momentum (which depends on the time step)? Or any other suggestions? My temporal discretization method comes straight from the Versteeg textbook.


July 1, 2016, 21:30 

#8 
Senior Member
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 22 
What is the nature of the instability? flow variables going to infinity? If so, which one starts it? Or do linear systems fail to converge? If so, which field? Or does the nonlinear iteration fail to converge?


July 4, 2016, 06:44 

#9 
New Member
DSR
Join Date: Sep 2014
Posts: 11
Rep Power: 8 
The pressure correction equation goes unstable first. If I just do heat conduction, I can use a large time step. The instability seems to occur during velocitypressure coupling.


July 4, 2016, 06:58 

#10  
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 5,688
Rep Power: 60 
Quote:
just check after one time step if your pressure solution appears regular and if the gradients added to the partial velocity field satisfy the continyuity equation. 

July 4, 2016, 07:08 

#11 
New Member
DSR
Join Date: Sep 2014
Posts: 11
Rep Power: 8 
If my time step is too big (not really even that big), my pressure correction equation diverges (first time step, first loop).


July 4, 2016, 09:11 

#12 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 5,688
Rep Power: 60 

July 4, 2016, 09:53 

#13 
New Member
DSR
Join Date: Sep 2014
Posts: 11
Rep Power: 8 
Ok thanks. What should I expect for time step (based on cell size, velocity, etc.)?


July 4, 2016, 11:13 

#14 
Senior Member
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 5,688
Rep Power: 60 

July 4, 2016, 23:54 

#15 
Senior Member
Michael Prinkey
Join Date: Mar 2009
Location: Pittsburgh PA
Posts: 363
Rep Power: 22 
If you are using the SIMPLE method, then the algorithm is completely implicit and (in principle) there is no timestep limit. In fact, taking the limit of dt>inf should lead you to the normal steady SIMPLE formulation. This, of course, assumes that the SIMPLE algorithm (outer iterations) converge at each timestep.
I am aware of no theory...or even heuristic...that gives timestep guidance with regard to the convergence of SIMPLE. There are requirements on pressure and velocity underrelaxation factors for the steady state case under assumptions of simple models...no reactions, no buoyancy, etc. I agree with the others advice. Verify that the face fluxes after pressure correction application results in conservation of mass. If that is not true both globally and in every cell, then your first problem to solve is there. There may be other bugs, but you will not be able to check because the fluxes driving the convection of other properties will be garbage. And remember that because this is colocated, the cell values of velocity are (almost) irrelevant. The face fluxes are the primary metric with regard to pressure correction and mass conservation. 

Tags 
cfd, rhiechow, simple, time step 
Thread Tools  Search this Thread 
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
AMI speed performance  danny123  OpenFOAM  21  October 24, 2020 04:13 
Stuck in a Rut interDyMFoam!  xoitx  OpenFOAM Running, Solving & CFD  14  March 25, 2016 07:09 
[Other] Contribution a new utility: refine wall layer mesh based on yPlus field  lakeat  OpenFOAM Community Contributions  57  February 1, 2015 08:25 
How to write k and epsilon before the abnormal end  xiuying  OpenFOAM Running, Solving & CFD  8  August 27, 2013 15:33 
pisoFoam with kepsilon turb blows up  Some questions  Heroic  OpenFOAM Running, Solving & CFD  26  December 17, 2012 03:34 