Changes between Version 136 and Version 137 of FluxLimitedDiffusion


Ignore:
Timestamp:
03/31/13 10:52:15 (12 years ago)
Author:
Jonathan
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • FluxLimitedDiffusion

    v136 v137  
    467467\begin{eqnarray}
    468468E^{n+1}_i-E^{n}_i = & \left [ \alpha^n_{i+1/2} \left ( \psi E^{n+1}_{i+1} + \bar{\psi} E^{n}_{i+1}- \psi E^{n+1}_{i} - \bar{\psi} E^n_{i} \right ) - \alpha^n_{i-1/2} \left ( \psi  E^{n+1}_{i} + \bar{\psi} E^{n}_i - \psi E^{n+1}_{i-1} - \bar{\psi}E^{n}_{i-1} \right ) \right ] \\
    469  + & \left [ \zeta^n_{i+1/2} \left ( \psi E^{n+1}_{i+1} + \bar{\psi} E^{n}_{i+1} + \psi E^{n+1}_{i} + \bar{\psi} E^n_{i} \right ) - \zeta^n_{i-1/2} \left ( \psi E^{n+1}_{i} + \bar{\psi} E^n_i + \psi E^{n+1}_{i-1} + \bar{\psi}E^{n}_{i-1} \right ) \right ] \\
    470 - & \frac{1}{1-\psi \phi} \left [ \theta + \epsilon \left ( \psi E^{n+1}_i + \bar{\psi} E^n_i \right ) - \omega v^n_x \left ( \psi E^{n+1}_{i+1} - \psi E^{n+1}_{i-1} + \bar{\psi} E^n_{i+1}- \bar{\psi} E^n_{i-1} \right ) + \xi \left ( \psi E^{n+1}_i + \bar{\psi} E^{n}_i \right ) \right ] \\
     469 + & \left [ \zeta^n_{i+1/2} v^n_{x,i+1/2} \left ( \psi E^{n+1}_{i+1} + \bar{\psi} E^{n}_{i+1} + \psi E^{n+1}_{i} + \bar{\psi} E^n_{i} \right ) - \zeta^n_{i-1/2} v^n_{x,i-1/2}\left ( \psi E^{n+1}_{i} + \bar{\psi} E^n_i + \psi E^{n+1}_{i-1} + \bar{\psi}E^{n}_{i-1} \right ) \right ] \\
     470- & \frac{1}{1-\psi \phi} \left [ \theta + \epsilon \left ( \psi E^{n+1}_i + \bar{\psi} E^n_i \right ) + \omega v^n_x \left ( \psi E^{n+1}_{i+1} + \bar{\psi} E^n_{i+1} - \psi E^{n+1}_{i-1} - \bar{\psi} E^n_{i-1} \right ) - \xi \left ( \psi E^{n+1}_i + \bar{\psi} E^{n}_i \right ) \right ] \\
    471471\end{eqnarray}
    472472}}}
     
    478478and where
    479479
    480 [[latex(\zeta_{i+1/2}= \frac{\Delta t}{\Delta x}\frac{3-R_{2,i+1/2}}{2} v^n_x,{i+1/2} )]]
     480[[latex(\zeta_{i+1/2}= \frac{\Delta t}{\Delta x}\frac{3-R_{2,i+1/2}}{2} )]]
    481481
    482482and
     
    490490and
    491491
    492 [[latex(\theta = f(e^n_i) = \epsilon^n_i \frac{4 \pi}{c} B \left ( T^n_i \right ) )]]
     492[[latex(\theta = \Delta t f(e^n_i) = \epsilon^n_i \frac{4 \pi}{c} B \left ( T^n_i \right ) )]]
    493493
    494494and