Changes between Version 111 and Version 112 of FluxLimitedDiffusion
- Timestamp:
- 03/28/13 13:36:51 (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
FluxLimitedDiffusion
v111 v112 474 474 and 475 475 476 where [[latex(\omega_{i} = \frac{\lambda_i \Delta t}{\Delta x} \left ( \frac{\kappa_{0P,i}}{\kappa_{0R,i}}-\frac{1}{2} \right ) = \frac{1}{2}\mathbf{F} \cdot \frac{\mathbf{v}}{c} \left ( \kappa_{0P,i}-\kappa_{0R,i} \right ))]] 477 478 476 where [[latex(\omega_{i} = \frac{\lambda_i \Delta t}{\Delta x} \left ( \frac{\kappa_{0P,i}}{\kappa_{0R,i}}-\frac{1}{2} \right ))]] 477 478 where 479 480 [[latex(\frac{\omega}{\alpha} = \frac{v}{c} \kappa_{0R,i} \Delta x)]] 479 481 === Time Discretization === 480 482 … … 536 538 So we would just modify [[latex(\alpha)]] and zero out the matrix coefficient to the ghost zone 537 539 540 541 538 542 === Constant Slope Boundary === 539 543 Here we want the flux to be constant so energy does not pile up near the boundary. If we cancel all derivative terms on both sides of the cell, this will effectively match the incoming flux with the outgoing flux. This can also be done by setting [[latex(\alpha_g = \alpha_i = 0)]]