Changes between Version 54 and Version 55 of u/erica/RoeSolver


Ignore:
Timestamp:
05/17/13 12:33:16 (12 years ago)
Author:
Erica Kaminski
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • u/erica/RoeSolver

    v54 v55  
    145145= Short-comings of the Roe Solver =
    146146
     147As with all methods of linearizing the Euler equations, discontinuities are resolved, but continuous fan-like structures are not. Thus the Roe solver has trouble solving the RP inside of a sonic (aka transonic) rarefaction wave (see below figure), but less of a problem with contacts and shocks. This problem is refered to as the 'entropy problem', because the method treats the rarefaction wave as a 'rarefaction shock', which is entropy violating. Recall that the entropy condition that must be satisfied for a discontinuous wave jump is that the characteristics ''run into'' each other. That is, Sb>S>Sa, where Sb is the speed BEHIND the wave, S is the speed of the wave, and Sa is the speed ahead of the wave. There was a later 'entropy fix' developed for the Roe solver that corrected this issue.
     148
     149Here is a diagram of the different types of waves:
     150
     151[[Image(waveStructures.png, 50%)]]
     152
     153Here is an example of the Roe solver poorly computing the flow inside of a left rarefaction. As you can see, it computes the solution as if there was a discontinuity inside of the wave, not a continuous fan.
     154
     155[[Image(Test1Roe.png)]]
     156
     157
     158
     159
    147160= The code =
    148161