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

Writing conservaion equations with a control volume approach with staggered grids

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

Like Tree2Likes
  • 1 Post By FMDenaro
  • 1 Post By leflix

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 8, 2012, 19:35
Default Writing conservaion equations with a control volume approach with staggered grids
  #1
New Member
 
Fabien
Join Date: Aug 2012
Posts: 3
Rep Power: 14
Fabien is on a distinguished road
I have a very fundamental question. Lets say that I have a medium (1D) and I want to write the momentum conservation equation and continuity (mass balance) equation for this medium. Based on CFD books, first I partition this medium with a grid (lets say uniform) which the boundaries of the girds coincide the left and right boundaries of the medium (the grid with solid lines). This is the gird that I solve pressure at the center point (grid for scalar properties) of it. Then I draw another grid but staggered half a gridblock and write the x-momentum on these grids (the grid with dashed lines). For writing the momentum and mass balance equations, I follow control volume approach which I evaluate the integral over surface boundary of v.n for each surface of gridblcoks. Everything looks fine and I solve for velocity at the boundary of the gridblock (as I need it to compute the convective flux) and I solve for pressure at the center of my original grid.

Now my question is that, DO I CONSERVE MOMENTUM ON THE whole MEDIUM? Because if you look at the grids, we see that there are half gridblcoks on the left and right sides close to the boundary that I don’t write momentum balance for these volumes (the volumes that are shown with a shaded area in the figure)!!! I understand that the effect of boundary is felt (lets say the boundary is wall and u_1=-u_ghost and u_n=-u_ghost ) but we don’t really write momentum blanace equation for these half-grid volumes. Does anyone have an explanation for this?
Please see the below figure for a better illustration.

Thank you very much


http://imageshack.us/photo/my-images/15/67923582.jpg/

Fabien is offline   Reply With Quote

Old   August 9, 2012, 05:05
Default
  #2
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
Quote:
Originally Posted by Fabien View Post
I have a very fundamental question. Lets say that I have a medium (1D) and I want to write the momentum conservation equation and continuity (mass balance) equation for this medium. Based on CFD books, first I partition this medium with a grid (lets say uniform) which the boundaries of the girds coincide the left and right boundaries of the medium (the grid with solid lines). This is the gird that I solve pressure at the center point (grid for scalar properties) of it. Then I draw another grid but staggered half a gridblock and write the x-momentum on these grids (the grid with dashed lines). For writing the momentum and mass balance equations, I follow control volume approach which I evaluate the integral over surface boundary of v.n for each surface of gridblcoks. Everything looks fine and I solve for velocity at the boundary of the gridblock (as I need it to compute the convective flux) and I solve for pressure at the center of my original grid.

Now my question is that, DO I CONSERVE MOMENTUM ON THE whole MEDIUM? Because if you look at the grids, we see that there are half gridblcoks on the left and right sides close to the boundary that I don’t write momentum balance for these volumes (the volumes that are shown with a shaded area in the figure)!!! I understand that the effect of boundary is felt (lets say the boundary is wall and u_1=-u_ghost and u_n=-u_ghost ) but we don’t really write momentum blanace equation for these half-grid volumes. Does anyone have an explanation for this?
Please see the below figure for a better illustration.

Thank you very much


http://imageshack.us/photo/my-images/15/67923582.jpg/



Consider a medium going for x= 0 to L, uniformly discretized as xi=(i-1)dx, i=1,N+1

The, x1 and xN+1 correspond to the physical boundaries of the medium.

Velocity is colocated on such nodes, while pressure is staggered of half-size such that p2 is colocated at x2-dx/2 and so on. You solve pressure for i=2, N+1.

You will solve the momentum equation going from x2 to xN and the "telescopic" property of the FVM allows the conservation. Then, all depends on the physical boundary conditions imposed at x1 and xN+1.
You don't need ghost cells
FMDenaro is offline   Reply With Quote

Old   August 9, 2012, 11:14
Default
  #3
New Member
 
Fabien
Join Date: Aug 2012
Posts: 3
Rep Power: 14
Fabien is on a distinguished road
Quote:
Originally Posted by FMDenaro View Post
Consider a medium going for x= 0 to L, uniformly discretized as xi=(i-1)dx, i=1,N+1

The, x1 and xN+1 correspond to the physical boundaries of the medium.

Velocity is colocated on such nodes, while pressure is staggered of half-size such that p2 is colocated at x2-dx/2 and so on. You solve pressure for i=2, N+1.

You will solve the momentum equation going from x2 to xN and the "telescopic" property of the FVM allows the conservation. Then, all depends on the physical boundary conditions imposed at x1 and xN+1.
You don't need ghost cells

Thank you very much Flippo for kindly responding to my question. I understand the gridding, but my confusion is about the control volumes that I consider for grids x2 to xN. These control volumes do not cover my whole physical region. So, I guess the heart of the answer is that ""telescopic" property of the FVM allows the conservation". Do you know about a good reference which I use to read about this property (with more kind of physical and intuitive discussion).
My question arose, as I try to write the control volume equations for each defined control volumes for density and velocity. As I explained I come up with half control volumes in both ends close to the boundary that I dont know how to write the control volume equations for them. Note that its not about using already developed equations but I want to develop equations from the very beginning (integral forms) by myself.

Thank you very much again for your help
Fabien is offline   Reply With Quote

Old   August 9, 2012, 11:51
Default
  #4
Senior Member
 
Filippo Maria Denaro
Join Date: Jul 2010
Posts: 6,893
Rep Power: 73
FMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura aboutFMDenaro has a spectacular aura about
I understand your doubts. My opinion is that:
- FVM are automatically conservative (in discrete form) since they are derived from the integral form of the governing equation.
- The staggered grid I considered before ensure the (local and global) conservation of mass, which in incompressible flow is mandatory. The momentum equation is locally conserved, also globally conserved in the domain that cover all the FVs (excluding the strips near the boundary).
To ensure full conservation there are many techniques you can use, for example writing an equation on each volume adjacent to the boundary or using deformed volumes ore prescribing suitable BC.

You can read Cap.9 of Versteeg H.K., Malalasekera W. Introduction to computational fluid dynamics, also the book of Peric and Ferziger can help.
See also this paper of Morinishi

http://www.google.it/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&ved=0CEcQ FjAA&url=http%3A%2F%2Fciteseerx.ist.psu.edu%2Fview doc%2Fdownload%3Fdoi%3D10.1.1.32.4898%26rep%3Drep1 %26type%3Dpdf&ei=KsojUI3ENab54QS4mYGAAw&usg=AFQjCN GbVO4FjeXrrVG9zyNJE3pKAf5Bew&sig2=MjhiX25ATWBPSbu-QeryZg
Fabien likes this.
FMDenaro is offline   Reply With Quote

Old   August 9, 2012, 14:51
Default
  #5
New Member
 
Fabien
Join Date: Aug 2012
Posts: 3
Rep Power: 14
Fabien is on a distinguished road
Thank you very much !
Fabien is offline   Reply With Quote

Old   August 9, 2012, 22:11
Default
  #6
Senior Member
 
Join Date: Aug 2011
Posts: 272
Rep Power: 16
leflix is on a distinguished road
[QUOTE=Fabien;376236]
I understand the gridding, but my confusion is about the control volumes that I consider for grids x2 to xN. These control volumes do not cover my whole physical region. So, I guess the heart of the answer is that ""telescopic" property of the FVM allows the conservation".
As I explained I come up with half control volumes in both ends close to the boundary that I dont know how to write the control volume equations for them.


Taking the notation of Filipo and as he said it
U(1) and U(NX+1) are boundary conditions. You discretize the momentum and solve it for nodes i=2 to i=NX.
The drawback you mentioned (the fact that your control volumes do not cover the whole physical domain) comes from the cell centred configuration that you have implicitly chosen.

If you take a node centred configuration your control volumes for momentum will cover your physical domain, but this time the control volumes for scalar will not. Check the sketch below.

But whatever the configuration you chose it does not really matter.
The integration remains the same.
Attached Files
File Type: zip sketch.zip (66.7 KB, 10 views)
Fabien likes this.

Last edited by leflix; August 9, 2012 at 22:26.
leflix is offline   Reply With Quote

Reply

Tags
finite volume method, staggered grid

Thread Tools Search this Thread
Search this Thread:

Advanced Search
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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Guide: Writing Equations in LaTeX on the CFD Online Forums pete Site Help, Feedback & Discussions 27 May 19, 2022 04:19
Calculation of the Governing Equations Mihail CFX 7 September 7, 2014 07:27
air bubble is disappear increasing time using vof xujjun CFX 9 June 9, 2009 08:59
control volume based finite element araz Main CFD Forum 1 November 15, 2003 18:50
Equations for finite difference, element & volume W.L.C. Main CFD Forum 9 April 3, 2001 13:08


All times are GMT -4. The time now is 00:23.