| 140 | | My code essentially does the following: 1) sets up left and right local Riemann problems as usual, 2) for each intercell boundary determines if the contact is left or right traveling, 3) Chooses bases on the position of the contact, which wave would cross the intercell boundary (the left wave associated with the u-a Euler eigenvalue for a right traveling contact, and vice versa), 4) decides if this wave is a sonic rarefaction, and if so, applies the Harten- Hyman entropy fix (described below), 5) if the wave is not a sonic rarefaction it determines if it crossed the intercell boundary, and if so, calculates the Roe flux, 6) if the wave has not crossed the intercell boundary, then it can be thought of as not having disturbed the adjacent cell's state, and so we just take the numerical flux to be the physical flux of that adjacent cell. |
| | 140 | My Roe solver + entropy fix essentially does the following: 1) sets up left and right local Riemann problems as usual, 2) for each intercell boundary determines if the contact is left or right traveling, 3) based on the position of the contact, chooses which wave would cross the intercell boundary (the left wave associated with the (u-a) Euler eigenvalue for a right traveling contact, and vice versa), 4) decides if this wave is a sonic rarefaction, and if so, applies the Harten- Hyman entropy fix (described below), 5) if the wave is not a sonic rarefaction it determines if it crossed the intercell boundary, and if so, calculates the Roe flux, 6) if the wave has not crossed the intercell boundary, then it can be thought of as not having disturbed the adjacent cell's state, and so we just take the numerical flux to be the physical flux of that adjacent cell. |
