Changes between Version 100 and Version 101 of u/erica/RoeSolver


Ignore:
Timestamp:
05/31/13 14:50:53 (11 years ago)
Author:
Erica Kaminski
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • u/erica/RoeSolver

    v100 v101  
    152152for a left sonic rarefaction. For a right sonic rarefaction, the coefficients are slightly adjusted. Note that the entropy fix requires estimates for the velocity and sound speed of the star region, ustar and astarL. We get estimates of these using the Roe approach for linearizing the Euler equations. They are listed in Toro, section 11.4.3.
    153153
     154Here is a screen shot of my code:
    154155
    155 [[Image(Roe1.png, 50%)]]
    156 [[Image(Roe2.png, 50%)]]
    157 [[Image(Roe3.png, 50%)]]
    158156
     157[[Image(Roe1.png, 100%)]]
     158[[Image(Roe2.png, 100%)]]
     159[[Image(Roe3.png, 100%)]]
     160
     161and corresponding files, modules, etc., are here: [http://www.pas.rochester.edu/~erica/solver.f90 Exact Riemann Solver -module], [http://www.pas.rochester.edu/~erica/God_Roe.f90 Roe Solver - main program], http://www.pas.rochester.edu/~erica/GodRoeproblem.data problem.data file]. To compile, simply type '''''gfortran solver.f90 God_Roe.f90 -o program'''''.
    159162
    160163
    161164= Results, Entropy Fix Code =
    162165
    163 
     166The code works correctly. Here are the 4 tests in Toro verifying this.
    164167
    165168== Test 1 - Sod Test ==
     
    171174== Test 2 ==
    172175
    173 Fails with the Roe Solver.
     176The Roe solver fails with low density flows such as that generated in this test. This is the major short coming of the Roe Solver, that is not present for the other approximation methods we have studied such as PVRS models.
    174177
    175178== Test 3 ==
     
    192195[[Image(RoeUTest5.png, 35%)]]
    193196[[Image(RoePTest5.png, 35%)]]
    194 
    195 = Discussion =
    196 
    197 
    198