Changes between Version 67 and Version 68 of u/erica/RoeSolver


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

Legend:

Unmodified
Added
Removed
Modified
  • u/erica/RoeSolver

    v67 v68  
    1616
    1717{{{#!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}
    1919}}}
    2020
     
    2222
    2323{{{#!Latex
    24 \hat{A} = \frac{\partial \vec{F}}{\partial \vec{U}}
     24\hat{A} = \frac{\partial ~\vec{F}}~{\partial ~\vec{U}}
    2525}}}
    2626where, 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:
     
    124124
    125125{{{#!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>
    139139}}} 
    140140
    141141{{{#!Latex
    142 \vec{K}^{(3)} = <1, u+a, H+ua>
     142\vec{K}^{(3)} = <1, ~u+a, ~H+ua>
    143143}}}
    144144