36 | | So the goal is to compute the wave speeds and associated eigenvalues and eigenvectors of the Jacobian matrix. There are 2 methods by which we can do this: 1) The 'Roe' approach, which (non-trivially) constructs an averaged Jacobian directly that must satisfy conditions such as hyperbolicity and conservation, and 2), the newer 'Roe-Pike' approach, which avoids solving for the Jacobian and insteads develops algebraic expressions for the sought quantities based on averages of the initial data. It is the 2nd, more widely used, approach that we will explore here. |
| 36 | So the goal is to compute the wave speeds and associated eigenvalues and eigenvectors of the Jacobian matrix. There are 2 methods by which we can do this: 1) The 'Roe' approach, which (non-trivially) constructs an averaged Jacobian directly that must satisfy conditions such as hyperbolicity and conservation, and 2), the newer 'Roe-Pike' approach, which avoids solving for the Jacobian and instead develops algebraic expressions for the sought quantities based on averages of the initial data. It is the 2nd, more widely used approach, that we will explore here. |
| 37 | |
| 38 | = The Roe-Pike Approach = |
| 39 | |
| 40 | = A Sample Algorithm = |
| 41 | |
| 42 | = Short-comings of the Roe Solver = |
| 43 | |
| 44 | = The code = |
| 45 | |
| 46 | = Results = |
| 47 | |
| 48 | = Discussion = |
| 49 | |
| 50 | |