Changes between Version 7 and Version 8 of FluxLimitedDiffusion


Ignore:
Timestamp:
03/18/13 13:54:45 (13 years ago)
Author:
Jonathan
Comment:

—

Legend:

Unmodified
Added
Removed
Modified

  • FluxLimitedDiffusion

    v7 v8  
    100100
    101101
    102 This gives matrix coefficiencts
     102This gives matrix coefficients
    103103
    104104[[latex(\left ( 1 + \alpha^n_{i+1/2} + \alpha^n_{i-1/2} + \epsilon^n_i \right ) E^{n+1}_i - \left ( \alpha^n_{i+1/2} \right ) E^{n+1}_{i+1} - \left ( \alpha^n_{i-1/2} \right ) E^{n+1}_{i-1}=E^n_i+\frac{4\pi \epsilon^n_i}{c}B \left (T^n_i \right ) )]]
     105
     106and for 2D the matrix coefficients would be
     107
     108|| [[latex(1+\alpha^n_{i+1/2,j}+\alpha^n_{i-1/2,j}+\alpha^n_{i,j+1/2}+\alpha^n_{i,j-1/2}+\epsilon^n_i )]] || [[latex(E^{n+1}_{i,j})]] ||
     109|| [[latex(-\alpha^n_{i-1/2,j})]] || [[latex(E^{n+1}_{i-1,j})]] ||
     110|| [[latex(-\alpha^n_{i+1/2,j})]] || [[latex(E^{n+1}_{i+1,j})]] ||
     111|| [[latex(-\alpha^n_{i,j-1/2})]] || [[latex(E^{n+1}_{i,j-1})]] ||
     112|| [[latex(-\alpha^n_{i,j+1/2})]] || [[latex(E^{n+1}_{i,j+1})]] ||
     113
     114and the source term would be unchanged from 1D
     115
     116||  [[latex(E^n_i+\frac{4\pi \epsilon^n_i}{c}B \left (T^n_i \right ) )]]