Changes between Version 29 and Version 30 of u/erica/ExactRiemannSolver

Timestamp:

01/05/13 13:26:25 (12 years ago)

Author:

Erica Kaminski

Comment:

—

u/erica/ExactRiemannSolver

@@ -3 +3 @@
 The Riemann problem is an IVP for the Euler Equations, which consists of 2 constant initial data states separated by a discontinuity between them, say at x=0. For x<x=0, we say the data state is XL, for which there are initial variables WL=(rho_L, p_L, u_L). Similarly for the initial right data state, XR=x>x=0, the initial variables are WR=(rho_R, p_R, u_R). The solution of this IVP consists of 3 nonlinear waves, a left wave that is either a shock or a rarefaction, a center contact discontinuity, and a right wave that is either a shock or a rarefaction wave. Depending on which type of L- or R-wave is present, different expressions exist that describe the change in variables across them. The contact discontinuity is special in that the pressure (p) and velocity (u) are constant across it. 
 
-[http://www.pas.rochester.edu/~erica/solver.f90 here]
 
 == Program Outline for Exact Riemann Solver == 
@@ -10 +9 @@
 
 === Exact Riemann Solver ===
-Attached here is the Exact RS program I wrote. It reproduces the results for the Toro tests, as can be seen in the next section. The program is broken up into the following routines:
+Attached [http://www.pas.rochester.edu/~erica/solver.f90 here] is the Exact RS program I wrote. It reproduces the results for the Toro tests, as can be seen in the next section. The program is broken up into the following routines:
 
 1. Estimate an initial starting P*. This subroutine was taken from Toro almost exactly, as I have not yet covered the details for approximation methods in my studies. Once this P* is found, it is fed into the !FindPstar routine, whereby a Newton Rhapson iteration method corrects this P* until it converges.