
[Sponsors] 
April 14, 2010, 16:52 
NavierStokes spherical discretization

#1 
Senior Member
Join Date: Apr 2009
Posts: 118
Rep Power: 9 
Hi,
I'm trying to discretize the NavierStokes 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

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: 9 
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 NavierStokes 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

multiply r2 at both sides


April 15, 2010, 07:48 

#5 
Senior Member
Join Date: Apr 2009
Posts: 118
Rep Power: 9 
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(<rho*u''*u''>), 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

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: 9 
r^2 is not a constant, it arises by writing the NavierStokes 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: 225
Rep Power: 11 
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 NavierStokes 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: 9 
Sorry I should've explained it a bit more clearly
So I have the following 1D NavierStokes 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 1D 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

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: 9 
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 1D case such as mine? and would ? or is where is the radial value at the midpoint of the control volume? Last edited by lost.identity; April 16, 2010 at 15:42. 

April 16, 2010, 21:36 

#12 
Senior Member

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: 9 

April 18, 2010, 12:39 

#14 
Senior Member
Join Date: Apr 2009
Posts: 118
Rep Power: 9 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
Filtered navier stokes equation..LES::  Palani Velladurai  Main CFD Forum  7  September 6, 2013 02:51 
Navier Stokes, Pressure BC, Discretization  Sree  Main CFD Forum  0  August 1, 2007 04:55 
Navier stokes compresible viscid flow fea, somebody can help?  Jose Choy  Main CFD Forum  3  October 24, 2003 02:28 
Newbie:Viscoelasticity and Navier stokes equation  Rajil Saraswat  Main CFD Forum  2  June 9, 2003 07:21 
help: I am trying to solve Navier Stokes compressible and viscid flow  Jose Choy  Main CFD Forum  2  May 18, 2000 05:45 