|
[Sponsors] |
How to construct the variable on faces of control volume |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
November 5, 2014, 05:19 |
How to construct the variable on faces of control volume
|
#1 |
New Member
ZhengJJ
Join Date: Mar 2011
Posts: 7
Rep Power: 15 |
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 |
|
November 5, 2014, 08:26 |
|
#2 |
New Member
Parth Thaker
Join Date: Mar 2014
Location: Surat , Gujarat , INDIA
Posts: 26
Rep Power: 12 |
advection schemes are used for finding values at face from values at cell centers . . .
referances - book of anderson, book of patnakar |
|
November 5, 2014, 21:53 |
|
#3 | |
New Member
ZhengJJ
Join Date: Mar 2011
Posts: 7
Rep Power: 15 |
Quote:
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? |
||
November 6, 2014, 05:11 |
|
#4 |
New Member
Parth Thaker
Join Date: Mar 2014
Location: Surat , Gujarat , INDIA
Posts: 26
Rep Power: 12 |
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
|
|
November 6, 2014, 15:04 |
|
#5 |
Senior Member
|
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). |
|
Tags |
gradient fvm |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |