Skip to topic | Skip to bottom
Home

Start of topic | Skip to actions
Riemann (1860) presented a solution to the one-dimensional Euler equations in which the velocity and pressure are functions of density only. Landau in Section 101 gives the explict relations for a simple wave in a polytropic gas assuming that there is a point in the wave for which v =0
c=c_o \pm \frac{1}{2}(\gamma -1)v
\rho=\rho_o (1 \pm \frac{1}{2}(\gamma -1)v/c_o)^{2/(\gamma-1)}
P=P_o (1 \pm \frac{1}{2}(\gamma -1)v/c_o)^{2\gamma/(\gamma-1)}
.. with the velocity defined implicity by
x=t(\pm c_o + \frac{1}{2}(\gamma +1)v) + f(v)
or alternately
v=F[x-(\pm c_o   + \frac{1}{2}(\gamma +1)v)t]

A Fortran90 code that computes this exact solution using a Newton-Raphson method has been written and maybe found in the repository. The velocity spectra of the exact solution is also computed (using an FFT).

While the code may be easily modified, it currently assumes that the initial velocity profile is a single mode of the form

F[x,t=0] = u_o \sin[\pi x].
For such an initial condition it can be shown that the shock develops at the time
t = \frac{2}{u_o \pi (\gamma+1)}

An example simple wave an .avi animation shows the density up to breaking as computed by the exact solution, while a simulation using the TCD-WENO scheme computes past breaking. The parameters for this example were

 \rho_o = 1, P_o = 1, u_o = 1/4.
This exact solution with the is used to study the convergance of the skew-symetric form of the solver in the Simple Wave Simulations section.

-- DavidHill - 20 Dec 2004


You are here: Amroc > WenoHome > SimpleWave

to top

Copyright © 1997-2024 California Institute of Technology.