burgers1.f90

Taken from original preamble - wiki to be updated soon

solve the 1d inviscid burgers equation, in conservative form:

u_t + [0.5 * u**2 ]_x=0 u=u(x,t) (user defined at t=0)

using a second order finite volume predictor-corrector scheme and optional gradient limiters. Boundaries are (minimally-tested) zero gradient.

user controled parameters in subroutines "control" and "initial_conditions"

Visually seems to capture correct shock formation (e.g. use IC sinusoid1) and rarefaction (e.g. IC jump with uright>uleft). Not done a more quantatative comparison to analytics yet

I was a bit surprised in jump (shock) problems with no gradient limiter how well it seems to do...

The code is broken into more subroutines and modules than neccesary for such a simple problem. The logic of the subroutines mimics a more sophistic hydrodynamics code whihc would warrent such an approac