CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Main CFD Forum

SIMPLE method with explicit solver?

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

Reply
 
LinkBack Thread Tools Display Modes
Old   July 24, 2010, 18:07
Default SIMPLE method with explicit solver?
  #1
New Member
 
Scott
Join Date: Jul 2010
Posts: 4
Rep Power: 7
Scott2 is on a distinguished road
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
Scott2 is offline   Reply With Quote

Old   July 25, 2010, 00:09
Default
  #2
New Member
 
Join Date: Jul 2010
Posts: 8
Rep Power: 7
startingcfd is on a distinguished road
Related question:

Does anyone have a paper or lecture link that has a very basic and detailed explanation of the SIMPLE algorithm?
startingcfd is offline   Reply With Quote

Old   July 26, 2010, 07:45
Default
  #3
New Member
 
Dhileep
Join Date: Apr 2010
Posts: 10
Rep Power: 7
Dhileep T K is on a distinguished road
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.
Dhileep T K is offline   Reply With Quote

Old   July 27, 2010, 08:45
Default
  #4
New Member
 
Join Date: Jan 2010
Posts: 3
Rep Power: 7
thehexman is on a distinguished road
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.
thehexman is offline   Reply With Quote

Old   July 30, 2010, 12:13
Default Projection Method
  #5
Member
 
Hans Bihs
Join Date: Jun 2009
Location: Trondheim, Norway
Posts: 34
Rep Power: 8
valgrinda is on a distinguished road
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
valgrinda is offline   Reply With Quote

Old   August 7, 2010, 13:44
Default
  #6
New Member
 
Scott
Join Date: Jul 2010
Posts: 4
Rep Power: 7
Scott2 is on a distinguished road
Hans,
Thank you for your reply. It was extreamly helpful.

Scott
Scott2 is offline   Reply With Quote

Old   August 9, 2010, 06:08
Question SHIP FLOW help please
  #7
New Member
 
harsha
Join Date: Aug 2010
Posts: 4
Rep Power: 6
harsha00711 is on a distinguished road
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
harsha00711 is offline   Reply With Quote

Reply

Tags
explicit, implicit, monotone, simple

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
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 & Dual-Time 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


All times are GMT -4. The time now is 16:04.