
[Sponsors] 
July 24, 2010, 18:07 
SIMPLE method with explicit solver?

#1 
New Member
Scott
Join Date: Jul 2010
Posts: 4
Rep Power: 7 
I have a question about implementing the SIMPLE method using and explicit method for the momentum equations. All of the literature that I have found implies that the momentum equations for each direction be advanced implicitly, before solving the pressure correction equations and applying the velocity corrections. However, I would like to solve momentum equations using a direction split, explicit, monotone method. Is this a reasonable approach? Any help is greatly appreciated.
Thanks, Scott 

July 25, 2010, 00:09 

#2 
New Member
Join Date: Jul 2010
Posts: 8
Rep Power: 7 
Related question:
Does anyone have a paper or lecture link that has a very basic and detailed explanation of the SIMPLE algorithm? 

July 26, 2010, 07:45 

#3 
New Member
Dhileep
Join Date: Apr 2010
Posts: 10
Rep Power: 8 
Refer the following book for SIMPLE algorith.
AN INTRODUCTION TO COMPUTATIONAL FLUID DYNAMICS, BY VERSTEEG AND MALALASEKARA which will be very useful for the beginners. 

July 27, 2010, 08:45 

#4 
New Member
Join Date: Jan 2010
Posts: 3
Rep Power: 8 
Hi,
I have also been trying to implement a very small general flow solver using the SIMPLE/pressure correction method. As simplicity is more important than performance at this stage of my project, I used MATLAB to program it. In the future, I plan to implement an immersed boundary test case to gain some experience with this method. The code uses upwind and central differencing and is explicit in time. This seems to work quite well if the usual restrictions (CFL, Re etc) are upheld. One of the simplest ways to ensure convergence and stability of the simulation is of course to increase the amount of cells and decrease the time step. This has worked quite well, although the actual computation is of course slow. One of the major problems I am experiencing are the boundary conditions of the simulation. Away from the boundaries, 4 neighboring cells are used to compute the pressure in the center cell enclosed. However, computing the pressure at the boundary cells is difficult because there are only 3 neighboring cells and the pressure is unknown, therefore rendering an implementation using ghost cells useless. Does anyone have a good solution to calculating the pressures at the boundaries?  BTW, a good book explaining the pressure correction method is also: Numerical computation of internal and external flows by Charles Hirsch. Unfortunately, he advocates using implicit methods, but the equations can be changed on paper without much ado. 

July 30, 2010, 12:13 
Projection Method

#5 
Member
Hans Bihs
Join Date: Jun 2009
Location: Trondheim, Norway
Posts: 34
Rep Power: 8 
Hei,
When using explicit time discretization, the straightforward and natural treatment of the pressure is Chorin's Projection Method. Here the pressure is calculated after each timestep and then used to correct the velocities in order to achieve divergencefree velocities. This method does not work for implicit time stepping, since all velocities are from the new time step. When the velocities are corrected with the pressure, they are divergence free, but the momentum equation is not fulfilled anymore. The solution to this problem is an iterative procedure, where the velocities and the pressure are solved several times in each timestep, so that the velocities fulfill both the continuity and the momentum equation. The SIMPLE algorithm is build upon this idea, except it does not solve for the pressure, but for the pressure correction. Details can be found in the Ferziger & Peric book, chapter 7. Cheers Hans 

August 7, 2010, 13:44 

#6 
New Member
Scott
Join Date: Jul 2010
Posts: 4
Rep Power: 7 
Hans,
Thank you for your reply. It was extreamly helpful. Scott 

August 9, 2010, 06:08 
SHIP FLOW help please

#7 
New Member
harsha
Join Date: Aug 2010
Posts: 4
Rep Power: 7 
hi
i am using shipflow to model an amphibious tank. i have made the model in SHIPFLOW design, i made a reference plane and then i tried to generate a section group. the section group is not being created. it is giving the following error. 'INFO section calculation failure [boundary coincidence]' someone please help me in this matter. thank u very much 

Tags 
explicit, implicit, monotone, simple 
Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Converting implicit solver to explicit one  dohnie  OpenFOAM Programming & Development  1  October 1, 2009 17:46 
Explicit vs Implicit method in free surface flow  Farhad  Main CFD Forum  4  November 16, 2006 11:42 
simple & DualTime Stepping Method  simple  Main CFD Forum  0  July 29, 2006 08:17 
Low rate of Convergency in SIMPLE method  Ahmad Falahat Pisheh  Main CFD Forum  2  January 10, 2004 04:25 
SIMPLE Algorithm & Conjugate Gradient Method  Abhijit Tilak  Main CFD Forum  7  March 20, 2001 01:01 