Changes between Version 184 and Version 185 of FluxLimitedDiffusion

Timestamp:

04/03/13 17:23:32 (12 years ago)

Author:

Jonathan

Comment:

—

FluxLimitedDiffusion

@@ -545 +545 @@
 \begin{eqnarray}
 E^{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 ] \\
- & + & \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 ] \\
- & - & \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 ] \\
+ & - & \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 ] \\
+ & - & \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 ] \\
 \end{eqnarray}
 }}}
@@ -568 +568 @@
 and
 
-[[latex(\theta = \Delta t f(e^n_i) =  \Delta t 4 \pi \kappa^n_{0P,i} B \left ( T^n_i \right ) )]]
+[[latex(\theta = - \Delta t f(e^n_i) =  \Delta t 4 \pi \kappa^n_{0P,i} B \left ( T^n_i \right ) )]]
 
 and