Changes between Version 111 and Version 112 of FluxLimitedDiffusion


Ignore:
Timestamp:
03/28/13 13:36:51 (12 years ago)
Author:
Jonathan
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • FluxLimitedDiffusion

    v111 v112  
    474474and
    475475
    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 
     476where [[latex(\omega_{i} = \frac{\lambda_i \Delta t}{\Delta x} \left ( \frac{\kappa_{0P,i}}{\kappa_{0R,i}}-\frac{1}{2} \right ))]]
     477
     478where
     479
     480[[latex(\frac{\omega}{\alpha} = \frac{v}{c} \kappa_{0R,i} \Delta x)]]
    479481=== Time Discretization ===
    480482
     
    536538So we would just modify [[latex(\alpha)]] and zero out the matrix coefficient to the ghost zone
    537539
     540
     541
    538542=== Constant Slope Boundary ===
    539543Here 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)]]