Navier-Stokes spherical discretization

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

 April 14, 2010, 16:52 Navier-Stokes spherical discretization #1 Senior Member   Join Date: Apr 2009 Posts: 118 Rep Power: 10 Hi, I'm trying to discretize the Navier-Stokes equation written in the spherical form for a Finite Volume formulation. I'm having problems due to the 1/r^2 terms appearing in the equation. I'm following the method of Patankar "Numerical heat transfer and fluid flow" where he integrates the equation for the control volume. It is spherically symmetric so I'm only retaining the radial terms in my equations. Does anyone know how to do this? Thanks.

 April 14, 2010, 22:56 #2 Senior Member     p ding Join Date: Mar 2009 Posts: 337 Rep Power: 11 use the control volume center coordiante to compute the 1/r^2.

 April 15, 2010, 05:01 #3 Senior Member   Join Date: Apr 2009 Posts: 118 Rep Power: 10 My problem is actually with integrating the equation. I'm using Favre averaged version of the equation so I have the following convection term (1/r^2)*d/dr(r^2*rho*u*u) because u (the velocity) are functions of r I need to integrate this by parts. I tried multiplying the whole equation by r^2, but then this causes the same problem to other terms that doesn't have the (1/r^2) term. I've tried searching through the literature to find any papers that does finite volume method on spherical Navier-Stokes equation and haven't been able to find any that have done the same thing as what I'm trying to do.

 April 15, 2010, 06:16 #4 Senior Member     p ding Join Date: Mar 2009 Posts: 337 Rep Power: 11 multiply r2 at both sides

 April 15, 2010, 07:48 #5 Senior Member   Join Date: Apr 2009 Posts: 118 Rep Power: 10 Thanks for the reply. As I've already mentioned in my email, I've already tried that method. But the problem is when I multiply by r^2 on both sides I get terms like r^2*dp/dr, which again has to be integrated by parts. But since p is generally a function of r this causes problems in the integration. In addition I have more complex terms that arises due to Favre averaging such as r^2 d/dr(), where <> refers to Favre averages. I use the Boussinesq approximation to model this but still I get terms such as r*u, which again needs to be integrated by parts. Surely, there's another way to do this?

 April 15, 2010, 18:59 #6 Senior Member     p ding Join Date: Mar 2009 Posts: 337 Rep Power: 11 nope do not integrate by part. take r^2 is a constant

 April 16, 2010, 05:19 #7 Senior Member   Join Date: Apr 2009 Posts: 118 Rep Power: 10 r^2 is not a constant, it arises by writing the Navier-Stokes equation in spherical form. r is one of the variable dimensions. So when I'm integrating with respect to r, r^2 can't remain a constant.

 April 16, 2010, 08:57 #8 Senior Member   Join Date: Jul 2009 Posts: 247 Rep Power: 12 I'm not sure I'm following what you are actually trying to do at this point. It sounds as if you are trying to integrate the full Navier-Stokes equations analytically. If you want to do a finite volume discretization, then the 1/r**2 term will be a part of what you discretize - in that case you will use the value at the cell center, as ztdep has suggested. Treat it just like any other term in the discretized equation.

 April 16, 2010, 10:41 #9 Senior Member   Join Date: Apr 2009 Posts: 118 Rep Power: 10 Sorry I should've explained it a bit more clearly So I have the following 1-D Navier-Stokes equation written in spherical coordinates where Now what I want to do is integrate this momentum equation over the control volume and over time. Since it's 1-D I only consider a control volume around a single grid point P, the grid point East of point P is referred to as E and to the West is referred to as W. Let's say I want to integrate the term If the term wasn't there this would have been straightforward, which is the case in the Cartesian coordinates. So one thing I tried to do was to multiply the whole equation by . But then I get terms like which again gives me the same problem. This may be different to what's normally done in Control Volume analysis. I'm trying to use the method of Patankar "Numerical Heat Transfer and Fluid Flow", 1980. Thanks. Last edited by lost.identity; April 16, 2010 at 11:19.

 April 16, 2010, 12:07 #10 Senior Member     p ding Join Date: Mar 2009 Posts: 337 Rep Power: 11 r^2 tem represent the element volume in the sperical coordinate system Dv. as in cartesian coordinate, it is 1 dx dy dz.

 April 16, 2010, 13:41 #11 Senior Member   Join Date: Apr 2009 Posts: 118 Rep Power: 10 Thanks, I think I sort of know what you're saying. Are you saying that when I multiply by , and integrate then terms such as is this true for 1-D case such as mine? and would ? or is where is the radial value at the mid-point of the control volume? Last edited by lost.identity; April 16, 2010 at 15:42.

 April 16, 2010, 21:36 #12 Senior Member     p ding Join Date: Mar 2009 Posts: 337 Rep Power: 11 yes, that is what i mean/ by the way , how do you input the formular?

 April 17, 2010, 17:06 #13 Senior Member   Join Date: Apr 2009 Posts: 118 Rep Power: 10 THanks Just want to confirm, you mean The site uses Latex formatting to type equations. If you click any of the equations you'd see the code.

 April 18, 2010, 12:39 #14 Senior Member   Join Date: Apr 2009 Posts: 118 Rep Power: 10 Actually shouldn't it be where , P and E refers to the left and right hand side faces control volume respectively, and e is the midpoint of the control volume. which would finally give

 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 Palani Velladurai Main CFD Forum 7 September 6, 2013 02:51 Sree Main CFD Forum 0 August 1, 2007 04:55 Jose Choy Main CFD Forum 3 October 24, 2003 02:28 Rajil Saraswat Main CFD Forum 2 June 9, 2003 07:21 Jose Choy Main CFD Forum 2 May 18, 2000 05:45

All times are GMT -4. The time now is 18:56.