# SIMPLE method with explicit solver?

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

 July 24, 2010, 18:07 SIMPLE method with explicit solver? #1 New Member   Scott Join Date: Jul 2010 Posts: 4 Rep Power: 15 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: 15 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: 16 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: 16 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. granzer likes this.

 July 30, 2010, 12:13 Projection Method #5 Super Moderator   Hans Bihs Join Date: Jun 2009 Location: Trondheim, Norway Posts: 390 Rep Power: 18 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 divergence-free 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 granzer likes this.

 August 7, 2010, 13:44 #6 New Member   Scott Join Date: Jul 2010 Posts: 4 Rep Power: 15 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: 15 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