
[Sponsors] 
April 27, 2003, 16:12 
Godunov scheme

#1 
Guest
Posts: n/a

Hi
I was woundering what is the Godunov scheme ? Also compared to upwind, QUICK, central etc. Thanks in advance. Regards Peter 

April 28, 2003, 14:55 
Re: Godunov scheme

#2 
Guest
Posts: n/a

The Godunov scheme is the exact, but very expensive, solution to the Riemann problem at an interface. It solves the fully nonlinear implicit Jacobian, which is where the expense is entailed. An approximate Riemann solver freezes the Jocobian at the interface, and essentially linearises it, so that the problem may be solved faster, an example is the Roe approximate Riemann solver or the Osher scheme. Nobody bothers to solve the exact Riemann problem at every interface, since it is too expensive to evaluate these nonlinear flux functions.
You now have for instance chosen your scheme; lets take Roe. The terms you have mentioned, are purely the methods of doing the convective differencing, they solve the convective part of the Navierstokes equations. You should be careful that you do not refer to these as a scheme, they are not, they are just an algebraic expression for solving the convective terms of the equations, however, you will see all the time things like an upwind scheme or a central difference scheme, just beware. QUICK is from the MUSCL family, and is formulated for structured grids, it is meant to be second order accurate, and upwind. Upwind means that you weight the value at the interface based on the wind direction of the flow, in a proper TVD edge based solver this can be seen easily, but again as with the term scheme, the term upwind is used less strictly than it should be. Central refers to central difference, second order by construction; you can read up on differencing. I have tried to put the above in some short context, for instance if you want to find out about QUICK, then read about MUSCL, Van Leer was the inventor. I make the cautious warnings because when you write a solver, you will soon become aware that you are not writing an upwind scheme, or a quick scheme, or linear reconstruction, but are using these formulations to calculate the fluxes at an interface, in a higher order manner, and that these in tern will be used in the scheme, Roe, Osher, and so on, to calculate the convective part of the equations as well as being used to work out the correct amount of artificial viscosity that should be added to maintain a stable scheme, this is done via the Jacobian, which you have linearised and is essential to the scheme. In the context of the Roe scheme for example this Jacobian will manifests itself in your code as the Roe averages, at an interface. Hope this helps, rather brief but should provide a starting point. Cheers Andy 

April 29, 2003, 05:28 
Re: Godunov scheme

#3 
Guest
Posts: n/a

Hi,
You say "The Godunov scheme is the exact, but very expensive, solution to the Riemann problem at an interface. It solves the fully nonlinear implicit Jacobian, which is where the expense is entailed. An approximate Riemann solver freezes the Jocobian at the interface, and essentially linearises it, so that the problem may be solved faster, an example is the Roe approximate Riemann solver or the Osher scheme" Where can I read about this....for a beginner? 

April 29, 2003, 05:43 
Re: Godunov scheme

#4 
Guest
Posts: n/a

Try Randall LeVeque's book
Finite Volume Methods for Hyperbolic Problems (Cambridge Texts in Applied Mathematics) Actually the Riemann problem is only exactly solvable in one spacedimension (i.e. the classic dam burst problem in shallow water theory). In higher dimensions you have to resort to approximate solvers. Tom. 

April 29, 2003, 06:24 
Re: Godunov scheme

#5 
Guest
Posts: n/a

An excellent book that will put it in context as well as the approximate solvers like Roe, and completly different approches like ENO and WENO, and HLLC, can be found in Computational Gas Dynamics by Laney (might have spelt that wrongly). He also has a very nice web page where you can down load fortran files for different schemes and try them out.
Tom is right about the oned problem, however, even when you solve the approximate riemann problem you are still solving an approximate oned problem, be clear that you understand why we bother to use approximate riemann solvers, it is only due to the nonlinear implicit nature of the Jacobian. An active area of research is multidimensional upwind where you resolve for the absoulute direction of the flow at an interface, instead of taking it normal. This introduces the multidimensional nature into the upwind. Hope this helps, again the book is excellent. Andy 

May 6, 2003, 05:37 
Re: Godunov scheme

#6 
Guest
Posts: n/a

also E.Toro's book.


November 27, 2009, 11:51 

#7 
New Member
Nicholas F Camus
Join Date: Sep 2009
Location: London, England.
Posts: 21
Rep Power: 9 
"The Godunov scheme is the exact, but very expensive, solution to the Riemann problem at an interface. It solves the fully nonlinear implicit Jacobian, which is where the expense is entailed..."
This is not true. A Godunov scheme is an upwind numerical method that is positively conservative. IT does not solve the Riemann problem. The Riemann problem stems from the use of a Godunov scheme, it is the solution of the Riemann problem that provides you with the intercell Godunov fluxes... The books you suggest however, are correct. Cheers. 

Thread Tools  
Display Modes  


Similar Threads  
Thread  Thread Starter  Forum  Replies  Last Post 
2nd order upwind scheme (Fluent and CFX)  Far  FLUENT  0  May 22, 2011 01:50 
how to understand high resolution scheme and high order scheme  iilw1314  Main CFD Forum  3  September 30, 2010 10:43 
Implementation of QUICK scheme  Romuald Skoda  Main CFD Forum  10  August 3, 2010 12:23 
Definition of limiter function for central dirrerencing scheme  sebastian_vogl  OpenFOAM Running, Solving & CFD  0  January 5, 2009 12:08 
extrapolation in MUSCL scheme  Chandra  Main CFD Forum  6  February 14, 2007 12:21 