
[Sponsors] 
January 13, 2009, 12:50 
To reproduce the behavior:
op

#1 
Senior Member
Maka Mohu
Join Date: Mar 2009
Posts: 305
Rep Power: 9 
To reproduce the behavior:
openFOAM version: 1.3. channel395 tutorial. If we change fvSchemes dictionary to the description at the end of this message to achieve 4th order accuracy and rerun channel395 tutorial, will finish the simulation with max Courant number of: 0.99 while using 2nd order schemes it finishes with 0.35. Q1) Is there something wrong in the selection of schemes made below? Q2) Is the three fold increase in Courant number a bug or such increase is the normal price in terms of stability that we pay for higher order schemes? Q3) Henry before in one thread mentioned that 2nd order schemes should be fine "particularly when the mesh is nonorthogonal". Can any body highlight the relation between the choice of scheme and the nonorthogonal correction in OpenFOAM.  fvSchemes: divSchemes { default none; div(phi,U) Gauss cubic; //linear; div(phi,k) Gauss limitedCubic 1; //limitedLinear 1; div(phi,B) Gauss limitedCubic 1; limitedLinear 1; div(B) Gauss cubic; //linear; div(phi,nuTilda) Gauss limitedCubic 1; //limitedLinear 1; div((nuEff*dev(grad(U).T()))) Gauss cubic; //linear; } laplacianSchemes { default none; laplacian(nuEff,U) Gauss cubic corrected; //linear corrected; laplacian(1A(U),p) Gauss cubic corrected; //linear corrected; laplacian(DkEff,k) Gauss cubic corrected; //linear corrected; laplacian(DBEff,B) Gauss cubic corrected; //linear corrected; laplacian(DnuTildaEff,nuTilda) Gauss cubic corrected; //linear corrected; } interpolationSchemes { default cubic; //linear; interpolate(U) cubic; //linear; } snGradSchemes { default corrected; } fluxRequired { default no; p; } }} Best regards, Maka. 

January 13, 2009, 13:12 
I would guess that the increas

#2 
Senior Member
Join Date: Mar 2009
Posts: 854
Rep Power: 13 
I would guess that the increase in Courant number is a consequence of instabilities in the velocity field. Have you looked at the results in detail? Are there unphysical spikes in the velocity when running cubic differencing?
Note that running such scheme does not provide overall 4thorder accuracy because the Gauss integrals are 2ndorder accurate. H 

January 13, 2009, 15:18 
(1) Yes, there are localized i

#3 
Senior Member
Maka Mohu
Join Date: Mar 2009
Posts: 305
Rep Power: 9 
(1) Yes, there are localized increase of velocity near the core of the channel but due to diffusion it is hard to say if they are physical or not because of the turbulent nature of the flow. I have seen the same behavior in other problems using LES and tried to reproduce it with the the channel tutorial. I agree that in some case the spikes are very clear and are unphysical. If cubic scheme is unbounded then that explains such behavior.
(2) regarding the overall 4thorder accuracy (actually that was my objective): I used to think that when we replace volume integral with surface integral in FVM using Gauss theorem, the resulting equation is exact till we approximate the surface terms. I used to think that Gauss that we write in front of the scheme (which is the only option for divergence and Laplacian schemes) refers to such process of replacing the volume integrals but surface integral. I would be grateful if you explain where such limitation of overall 2nd order accuracy in Gauss integrals come from? Many thanks Henry for your help. Best regards, Maka. 

January 13, 2009, 16:42 
(1) All schemes other than upw

#4 
Senior Member
Join Date: Mar 2009
Posts: 854
Rep Power: 13 
(1) All schemes other than upwind are potentially unbounded and the higher the order the more likely the unboundedness and the larger number of computational modes supported.
(2) Gauss theorem allows us to convert volume integrals into surface integrals which if evaluated exactly would allow us to evaluate the volume integrals exactly. However, currently we use interpolation to obtain values at the facecentres (on regular meshes) which are the Gausspoints for evaluating the surface integrals to secondorder simply by multiplying the faceareas. Hence irrespective of the order of the interpolation the order of the volumeintegrals is second on a regular mesh. Nevertheless it is still often beneficial to use higherorder schemes for interpolation, particularly for convection terms. Note also that irregular meshes the above procedure is not usually even secondorder because the simple interpolation does not provide values at the facecentres but this can be corrected for and we are developing a new set of polynomialfit interpolation schemes which are formally second and higher order on arbitrary meshes. H 

January 14, 2009, 04:08 
Many thanks for your explanati

#5 
Senior Member
Maka Mohu
Join Date: Mar 2009
Posts: 305
Rep Power: 9 
Many thanks for your explanations Henry.
Best regards, Maka. 

January 17, 2009, 04:55 
Prof. Weller,
one of my objec

#6 
Senior Member

Prof. Weller,
one of my objective is to implement a linearized euler solver for acoustic propagation in OpenFOAM. Now, as you said, it's not possible to achieve an integration order higher than 2 because of the Gauss integration. For my purposes, I need 4th order non dispersive schemes, so, with your polynomial fitting it's possible to achieve it in a FVM framework? Thanks, Ivan 

January 19, 2009, 09:13 
It is possible to implement th

#7 
Senior Member
Join Date: Mar 2009
Posts: 854
Rep Power: 13 
It is possible to implement the FVM higher than secondorder in general and in OpenFOAM by interpolating to more than one Gausspoint on the cell faces in order to obtain a higherorder surface integral. However, for our current purposes this has not proved necessary and we can achieve the accuracy we need with highorder interpolation to the cell face centres only. Do you think that this might be suitable for acoustic propagation or do you think that formally higherorder integration is also necessary?
Another option would be to use the polynomial fitting within a finitedifference framework to provide a formally higherorder scheme if local conservation is less important than formal accuracy. H 

January 19, 2009, 09:31 
Dear Prof Weller,
I have allr

#8 
Senior Member

Dear Prof Weller,
I have allready tried to implement an aeroacoustic solver (linearized euler) with the allready existing schemes (for istance, I tried linear and cubic), but I can't get good results, because linear is too dissipative and cubic is unstable and not so good with the dispersion error, so I gave up and wait for new evolutions. The state of art in aeroacoustic is DRP (Dispersion Relation Preserving) 4th order finite difference schemes with filters (huge stencils!) or Discontinuous Galerkin schemes, there are people (Davidson at Chalmer, as I know) that work on finite difference  like FVM schemes, but they need structured cells. So, conservation is not as fundamental as for istance in LES. Can you give me some reference on polynomial fitting schemes? Do you manage to put it in some new release of OF? Thank you 

January 19, 2009, 12:44 
Do you need highorder centred

#9 
Senior Member
Join Date: Mar 2009
Posts: 854
Rep Power: 13 
Do you need highorder centred or upwinded schemes? My guess is that cubic is unstable because it is a centred scheme and some kind of high frequency filtering would have to be applied.
We are developing both highorder centred and upwinded schemes and these will be included in the next release of OpenFOAM. Some information on these schemes is contained in: "Voronoi, Delaunay and Block Structured Mesh Refinement for Solution of the Shallow Water Equations on the Sphere" Hilary Weller, Henry Weller and Aime' Fournier Submitted to Monthly Weather Review, 2008 http://www.met.reading.ac.uk/~sws02h...rnier_2008.pdf and more details can be found in the references. H 

January 19, 2009, 12:58 
Someone uses centered schemes

#10 
Senior Member

Someone uses centered schemes with selective filters, someone uses specific low dissipation upwind schemes. Now, me and my research group have a Fortran code based on Discontinuous Galerkin that is 4th order with an upwind treatment of the interfaces fluxes (LaxFriedrichs fluxes), and with the acoustics test cases it works quite well.
It would be great to have Highorder schemes in the next release, I hope it will be out as soon as possible! 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
2nd order conservative schemes  taw  Main CFD Forum  1  September 16, 2008 07:05 
CFL condition for higher order schemes  Shyam  Main CFD Forum  2  February 14, 2008 15:24 
Highorder discretisation schemes  Hu  Phoenics  2  April 25, 2002 20:41 
Higher Order FV Schemes for unstructured meshes  Apurva Shukla  Main CFD Forum  4  December 15, 2000 10:17 
Gradient Estimation for Higher Order Schemes  jianxia  Main CFD Forum  0  June 6, 2000 19:40 