Changes between Version 30 and Version 31 of u/erica/MusclHancock
 Timestamp:
 06/18/13 11:47:17 (11 years ago)
Legend:
 Unmodified
 Added
 Removed
 Modified

u/erica/MusclHancock
v30 v31 26 26 [[latex($ \frac{\partial \vec{u}}{\partial t} + \frac{\partial f(\vec{u})}{\partial x} = 0 $)]] 27 27 28 where u is the vector of conserved variables and f is the flux function. Adhering to the IVBP known as the Riemann Problem, we solve this system of equations using the fully discrete, explicit, conservative formula:28 where u is the vector of conserved variables (u=u(x,t)), and f is the flux function. Adhering to the IVBP known as the Riemann Problem (i.e. constant left and right data states separated by an initial discontinuity), we solve this system of equations using the fully discrete, explicit, conservative formula: 29 29 30 30 [[latex($ u_i^ {n+1} = u_i ^n + \frac{\triangle t}{\triangle x}[f_{i1/2}  f_{i+1/2}]$)]] … … 40 40 [[Image(MusclReconstruction.png, 35%)]] 41 41 42 Using these end values, one has the new Riemann problem: 43 44 [[latex($ \frac{\partial \vec{u}}{\partial t} + \frac{\partial f(\vec{u})}{\partial x} = 0 $)]] 45 46 [[latex($ u(x,0) = { \frac{u_i^r, x\<0}{u_{i+1}^l, x\>0} $)]] 47 48 where 0 is the local origin, i.e. any given intercell boundary. 49 50 51 42 52 At this point, we are left with a higher order accurate code, but not one that is free of spurious oscillations near large gradients. To circumvent this, we need to add a TVD measure. 43 53