Changes between Version 67 and Version 68 of u/erica/RoeSolver
- Timestamp:
- 05/17/13 16:50:12 (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/erica/RoeSolver
v67 v68 16 16 17 17 {{{#!Latex 18 \vec{U}_t + \frac{\partial \vec{F}} {\partial \vec{U}} \vec{U}_x = \vec{0}18 \vec{U}_t + \frac{\partial \vec{F}} {\partial \vec{U}} ~\vec{U}_x = \vec{0} 19 19 }}} 20 20 … … 22 22 23 23 {{{#!Latex 24 \hat{A} = \frac{\partial \vec{F}}{\partial\vec{U}}24 \hat{A} = \frac{\partial ~\vec{F}}~{\partial ~\vec{U}} 25 25 }}} 26 26 where, A-hat is a Jacobian matrix of averaged/constant values, we can derive an expression for the numerical flux in terms of 1) wave strengths (alpha), and the 2) eigenvalues (lambda) and 3) right eigenvectors (K) of the 'averaged' Jacobian: … … 124 124 125 125 {{{#!Latex 126 \vec{F}_{i+1/2} = <\rho_{i+1/2}, u_{i+1/2},p_{i+1/2}>127 }}} 128 129 {{{#!Latex 130 \vec{F}_L = <\rho_L, \rho_L * u_L,0.5\rho_L*u_L^2 + \frac{p_L}{\gamma - 1}>131 }}} 132 133 {{{#!Latex 134 \vec{K}^{(1)} = <1, u-a,H-ua>135 }}} 136 137 {{{#!Latex 138 \vec{K}^{(2)} = <1, u,0.5 u^2>126 \vec{F}_{i+1/2} = <\rho_{i+1/2}, ~u_{i+1/2}, ~p_{i+1/2}> 127 }}} 128 129 {{{#!Latex 130 \vec{F}_L = <\rho_L, ~\rho_L * u_L, ~0.5\rho_L*u_L^2 + \frac{p_L}{\gamma - 1}> 131 }}} 132 133 {{{#!Latex 134 \vec{K}^{(1)} = <1, ~u-a, ~H-ua> 135 }}} 136 137 {{{#!Latex 138 \vec{K}^{(2)} = <1, ~u, ~0.5 u^2> 139 139 }}} 140 140 141 141 {{{#!Latex 142 \vec{K}^{(3)} = <1, u+a,H+ua>142 \vec{K}^{(3)} = <1, ~u+a, ~H+ua> 143 143 }}} 144 144