# finite difference method for navier-stokes problem

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

 February 14, 2003, 01:14 finite difference method for navier-stokes problem #1 dallybird Guest   Posts: n/a hi all, I have a problem in solving nonlinear PDE with finite difference method(unstaggered scheme).It is an 2D flow through a contraction in a pipe. Finally I ended up with 2 governing equation(x and y momentum equations) and 1 continuity euqation:du/dx+dv/dy=0, where u and v are velocity and x and y directions. The two momentum equations contain velocity terms(linear and non-linear) and pressure terms. To my understanding, they needs iterative methods with an initial guess of velocity(u,v) and pressure(p). here comes the problem: with an initial guess of (u,v) field and pressure(p),I got a new set of (u,v). then how should I continue? what is the function of the continuity equation? I guess I need it to update the pressure(p),and continue using new sets of (u,v,p) to proceed with first two momentum equations , and so on. but I cant see how to get pressure updated with continuity equation. I read the book by H.K Versteeg and W Malalasekera, which used finite volumn methods,but it wont help me with finite difference methods, right? anyone had any suggestion on the problem or any helpful resource? thanks in advance.

 February 14, 2003, 10:26 Re: finite difference method for navier-stokes pro #2 Jim Park Guest   Posts: n/a For 2-d flow through a pipe (variation along the axis and radially from center line to wall edge), I would expect a cylindrical coordinate system to be used. Then the continuity equation would be du/dr + u/r + dv/dz = 0 with r being the radius and z being the axial coordinate. To do an axisymmetric pipe flow in Cartesian (x,y) coordinates, you'll need to do a 3-d calculation with curved boundaries! The extra term, u/r, arises from the cylindrical coordinate system. Similar changes occur in the momentum equations. Patankar's book describes one way to extract pressure from the coupling of momentum and continuity. Earlier work at Los Alamos (MAC, SMAC, the SOLA schemes) show the original details. There are detail differences between finite difference and finite volume implementation of these schemes (and formal accuracy estimates may be different), but either should be workable. All of the schemes mentioned above were developed for a staggered mesh. I can't advise on the details of moving them to an unstaggered mesh, although apparently the work of Rhie & Chou (spelling?) addresses this problem. Hope I haven't misunderstood your problem.

 February 14, 2003, 13:47 Re: finite difference method for navier-stokes pro #3 ryoga Guest   Posts: n/a Here are two references that you may want to consider: 1. Computational Fluid Mechanics and Heat Transfer by Tannehill, Anderson, and Pletcher 2. Finite Analytic Method in Flows and Heat Transfer by Chen, Bernatz, Carlson and Lin The second reference gives pretty specific details for implementing SIMPLE methods on both staggered and non-staggered grids. You might also want to considering "reformulating" the continuity equation into a Poisson pressure equation. This will give you a "definitive" equation for the pressure (just like the u momentum equation is a "definitive" equation for the u velocity, etc.) This formulation has nice convergence and stability properties however, continuity is not always exactly satisfied and large dilation (dilitation) values may appear (so it's a tradeoff!) Hope this helps! Ryoga

 February 14, 2003, 16:15 Re: finite difference method for navier-stokes pro #4 Tony Guest   Posts: n/a You may also search for the Uzawa iterations. A little bit slower than the approximate projection method as Ryoga suggested, but the divergence can be made free to machine accuracy.

 February 17, 2003, 09:26 Re: finite difference method for navier-stokes pro #5 Manosh Guest   Posts: n/a Hi There! Another alternative I can suggest you. It is possible to modify Continuity & Momentum equations from (u,v,p) to stream function and vorticity forms. In this case you only have two unknowns which are stream function and vorticity, no headache about pressure, P. For 2D case this process works nicely even in DNS methods. Manosh

 February 17, 2003, 23:00 Re: finite difference method for navier-stokes pro #6 dallybird Guest   Posts: n/a thanks a bunch guys! I am now learning to use vorticity-stream function approach.

 Thread Tools Display Modes Linear 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 On Pingbacks are On Refbacks are On Forum Rules

 Similar Threads Thread Thread Starter Forum Replies Last Post m-fry Main CFD Forum 1 April 20, 2010 14:40 misabel Main CFD Forum 0 February 10, 2010 07:06 gregmason Main CFD Forum 2 July 8, 2009 11:53 Luciane Grossi Main CFD Forum 1 July 27, 2000 23:28 C-H Kuo Main CFD Forum 3 November 5, 1998 10:07

All times are GMT -4. The time now is 20:06.