# SIMPLE Algorithm Finite Difference Equations: how to discretize and solve?

 Register Blogs Members List Search Today's Posts Mark Forums Read August 3, 2015, 12:12 SIMPLE Algorithm Finite Difference Equations: how to discretize and solve? #1 New Member   Justin Join Date: Jan 2015 Posts: 28 Rep Power: 10 I'm attempting to plot velocity and pressure profiles for air flow with around a NACA airfoil using the SIMPLE procedure to solve the incompressible Navier-Stokes equations. So far I've generated an elliptic grid and transformed the equations pertinent to the SIMPLE procedure to the computational plane. Where I'm completely confused is how to discretize and solve velocities and pressure. From what I've read, which is a lot of CFD books, I'll be solving the momentum equations for the velocities by the alternating direction implicit (ADI) method. Does this apply to all steps of the SIMPLE procedure (i.e. solving for velocities, solving the pressure correction, and solving velocity corrections)? Or solely the solution to the momentum equations? The part that really has me confused is how to treat the nonlinear terms. For example, the 2-D x-momentum Navier-Stokes equation is Applying the ADI method, I'll first determine the solution for a half time step implicitly in the direction and explicitly in the direction and then a next half time step explicitly in and implicitly in . When I disretize the momentum equations, solving via ADI will leave me with a multidiagonal system. But I'm lost on exactly how to solve because I'll have terms such as and as part of the solution vector. I'm actually lost on how to organize the equation to set up as a matrix equation in the first place. I've written out the equations, which has spanned many pages, yet still do not see how to proceed. My last question is how to discretize the following two terms: and . I've seen in several books where they let but I'm unsure how that gets incorporated into the solution. As for both the nonlinear terms I mentioned, will I basically have and values for each point in the domain? Like I said, I'm really confused. I'm also a newbie to CFD so please take it easy on me. I appreciate any advice. Thanks!  Thread Tools Search this Thread Show Printable Version Email this Page Search this Thread: Advanced Search Display Modes Linear Mode Switch to Hybrid Mode Switch to Threaded Mode 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 OffTrackbacks are Off Pingbacks are On Refbacks are On Forum Rules