Changes between Version 8 and Version 9 of u/erica/2D_Godunov


Ignore:
Timestamp:
08/07/13 14:50:14 (11 years ago)
Author:
Erica Kaminski
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • u/erica/2D_Godunov

    v8 v9  
    33= Introduction =
    44
    5 The method I chose for practice writing a higher dimension code is the 1st order accurate, split scheme, Godunov + Exact Riemann Solver to solve the 2D Euler equations for a cylindrical explosion. Given the higher dimension of the problem, new types of waves are present in the solution, namely shears. Shears are passively advected with the flow, as can be shown by combining the continuity equation and the corresponding momenta equations.
     5The method I chose for practice writing a higher dimension code is the 1st order accurate, split scheme, Godunov + Exact Riemann Solver to solve the 2D Euler equations for a cylindrical explosion. Given the higher dimension of the problem, new types of waves are present in the solution, namely shears. Shears are passively advected with the flow, as can be shown by combining the continuity equation and the corresponding momenta equations. This means that the velocity of the shear (i.e. velocity in the tangential direction relative to wave propagation) does not change across either shock or rarefaction. This means that the value of the tangential velocity only changes discontinuously across the middle wave of the Riemann problem, i.e. the contact discontinuity. The value of the shear velocity on the left side of the contact is just the tangential velocity of the left state, the value of the shear velocity on the right side of the contact is of the right state. This makes for an extremely straightforward/easy addition to the exact Riemann solver.
    66
    77= Initialization =
    88
    99This method begins by initializing a 2D Cartesian grid, with a circle in the center. The primitive fluid variables (rho, x-velocity, y-velocity, pressure) are set inside and outside of this circle.
     10
     11Here are some plots of the initial condition. You will see a circle of high density/pressure in the center of the mesh, with zero velocity in both x and y.
     12
     13Density:
     14
     15[[Image(rhoInit3d.png, 35%)]]
     16
     17Velocity:
     18
     19[[Image(uInit.png, 35%)]]
    1020
    1121= The essence of split schemes =