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

How to construct the variable on faces of control volume

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 5, 2014, 05:19
Default How to construct the variable on faces of control volume
  #1
New Member
 
ZhengJJ
Join Date: Mar 2011
Posts: 7
Rep Power: 15
fengyunmanbu is on a distinguished road
Hello CFD-Onliners~
I am currently studying to program a finite volume solver in cell centered scheme. When calculate the gradient in cell by the method metioned as below.
http://www.cfd-online.com/Wiki/Gradient_computation
It's need to know the variables on the faces of control cell.

May I ask how?
ps: Is there any books or someting else can be recommended
fengyunmanbu is offline   Reply With Quote

Old   November 5, 2014, 08:26
Default
  #2
New Member
 
Parth Thaker
Join Date: Mar 2014
Location: Surat , Gujarat , INDIA
Posts: 26
Rep Power: 12
Parth04 is on a distinguished road
advection schemes are used for finding values at face from values at cell centers . . .

referances - book of anderson, book of patnakar
Parth04 is offline   Reply With Quote

Old   November 5, 2014, 21:53
Default
  #3
New Member
 
ZhengJJ
Join Date: Mar 2011
Posts: 7
Rep Power: 15
fengyunmanbu is on a distinguished road
Quote:
Originally Posted by Parth04 View Post
advection schemes are used for finding values at face from values at cell centers . . .

referances - book of anderson, book of patnakar
thank you for your reply
I find two method
First: set the average of values at two adjacent cell centers as the values at their interface.

Second: calculate the values at grid point by weighted mean of central values, then average the values at points of face as the value of the face.

The first one seems simply, but which one is better?
fengyunmanbu is offline   Reply With Quote

Old   November 6, 2014, 05:11
Default
  #4
New Member
 
Parth Thaker
Join Date: Mar 2014
Location: Surat , Gujarat , INDIA
Posts: 26
Rep Power: 12
Parth04 is on a distinguished road
there are many mathods, cental differencing (the first u said that) , first order upwind , second order upwind, QUICK, it depends on your problem and solver to select the method
Parth04 is offline   Reply With Quote

Old   November 6, 2014, 15:04
Default
  #5
Senior Member
 
sbaffini's Avatar
 
Paolo Lampitella
Join Date: Mar 2009
Location: Italy
Posts: 2,150
Blog Entries: 29
Rep Power: 39
sbaffini will become famous soon enoughsbaffini will become famous soon enough
Send a message via Skype™ to sbaffini
There is a clear difference between computing the face values for the convective fluxes and for the gradients. Simply because the former usually require the latter.

Thus, we speak here for the gradients (and not the fluxes).

According to the Fluent references, the second method (node based) is superior to using a simple average of adjacent cell values. Actually, i tested this one on my own and found it also inconsistent on some tetrahedral grids.

If you don't want to go node based (which indeed is quite involved), you can iterate on the cell based one. You add to the simple mean the contribution of the gradients from both cells. At the first iteration such contributions are null. Then they build up with the iterations. This does not necessarily converge (because the resulting computational stencil is growing) but 2-3 iterations will improve the basic method.

My personal experience is to use:

- least-squares with two cell layers (or maybe one by carefully playng with a weighted least squares)
- Radial Basis Functions. Interpolate one or two cell layers around your cell with a RBF and take its analytical derivative. Its simpler than the least squares (once you have everything set up) and can give you much higher accuracy with, possibly, only a minor overhead (they both use matrices to compute the final derivative).
sbaffini is offline   Reply With Quote

Reply

Tags
gradient fvm

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
how to set periodic boundary conditions Ganesh FLUENT 15 November 18, 2020 07:09
snappyhexmesh remove blockmesh geometry philipp1 OpenFOAM Running, Solving & CFD 2 December 12, 2014 11:58
foam-extend_3.1 decompose and pyfoam warning shipman OpenFOAM 3 July 24, 2014 09:14
[blockMesh] BlockMesh FOAM warning gaottino OpenFOAM Meshing & Mesh Conversion 7 July 19, 2010 15:11
[snappyHexMesh] external flow with snappyHexMesh chelvistero OpenFOAM Meshing & Mesh Conversion 11 January 15, 2010 20:43


All times are GMT -4. The time now is 14:17.