Changes between Version 191 and Version 192 of FluxLimitedDiffusion


Ignore:
Timestamp:
04/05/13 13:38:05 (12 years ago)
Author:
Jonathan
Comment:

—

Legend:

Unmodified
Added
Removed
Modified

  • FluxLimitedDiffusion

    v191 v192  
    534534or in discretized form
    535535
    536 [[latex(Delta t < \xi \frac{T^n_i \Delta t}{4 \Gamma \left | \left ( -\theta_i + \epsilon_i E^*_i + \omega_i v^n_{x,i} \left ( E^*_{i+1}-E^*_{i-1} \right ) - \xi_i E^* \right ) \right |} )]]
     536[[latex(\Delta t < \xi \frac{T^n_i \Delta t}{4 \Gamma \left | \left ( -\theta_i + \epsilon_i E^*_i + \omega_i v^n_{x,i} \left ( E^*_{i+1}-E^*_{i-1} \right ) - \xi_i E^* \right ) \right |} )]]
    537537
    538538
    … …  
    543543as
    544544
    545 [[latex(g = \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 ) )]]
     545[[latex(g = \epsilon_i \left ( \psi E^{n+1}_i + \bar{\psi} E^n_i \right ) + \omega_i v^n_{x,i} \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_i \left ( \psi E^{n+1}_i + \bar{\psi} E^{n}_i \right ) )]]
    546546
    547547which along with the other terms gives
    … …  
    549549#!latex
    550550\begin{eqnarray}
    551 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 ] \\
    552  & - & \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 ] \\
    553  & - & \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 ] \\
     551E^{n+1}_i-E^{n}_i & = & \left [ \alpha_{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_{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 ] \\
     552 & - & \left [ \zeta_{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_{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 ] \\
     553 & - & \frac{1}{1+\psi \phi_i}  \left [ - \theta_i  + \epsilon_i \left ( \psi E^{n+1}_i + \bar{\psi} E^n_i \right ) + \omega_i 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_i \left ( \psi E^{n+1}_i + \bar{\psi} E^{n}_i \right ) \right ] \\
    554554\end{eqnarray}
    555555}}}
    … …  
    616616#!latex
    617617\begin{eqnarray}
    618  & \left ( 1 + \psi \left ( \alpha^n_{i+1/2} +  \alpha^n_{i-1/2} + \zeta^n_{i+1/2} v^n_{x,i+1/2} - \zeta^n_{i-1/2} v^n_{x,i-1/2} + \frac{\epsilon - \xi}{1+\psi \phi} \right ) \right ) E^{n+1}_i   \\
    619 - & \left ( \psi \left (  \alpha^n_{i+1/2} - \zeta^n_{i+1/2} v^n_{x,i+1/2} - \frac{\omega v^n_x}{1+\psi \phi} \right ) \right ) E^{n+1}_{i+1}  \\
    620 - & \left ( \psi \left (  \alpha^n_{i-1/2} + \zeta^n_{i-1/2} v^n_{x,i-1/2} + \frac{\omega v^n_x}{1+\psi \phi} \right ) \right ) E^{n+1}_{i-1}  \\
    621 = & \left ( 1 - \bar{\psi} \left ( \alpha^n_{i+1/2} + \alpha^n_{i-1/2} + \zeta^n_{i+1/2} v^n_{x,i+1/2} -  \zeta^n_{i-1/2} v^n_{x,i-1/2} + \frac{\epsilon - \xi}{1+\psi \phi} \right ) \right ) E^{n}_i  \\
    622 + & \left ( \bar{\psi} \left (  \alpha^n_{i+1/2} - \zeta^n_{i+1/2} v^n_{x,i+1/2} - \frac{\omega v^n_x}{1+\psi \phi} \right ) \right ) E^{n}_{i+1}  \\
    623 + & \left ( \bar{\psi} \left (  \alpha^n_{i-1/2} + \zeta^n_{i-1/2} v^n_{x,i-1/2} + \frac{\omega v^n_x}{1+\psi \phi} \right ) \right ) E^{n}_{i-1}  \\
    624 + & \frac{\theta}{1+\psi \phi}   \\
     618 & \left ( 1 + \psi \left ( \alpha_{i+1/2} +  \alpha_{i-1/2} + \zeta_{i+1/2} v^n_{x,i+1/2} - \zeta_{i-1/2} v^n_{x,i-1/2} + \frac{\epsilon_i - \xi_i}{1+\psi \phi_i} \right ) \right ) E^{n+1}_i   \\
     619- & \left ( \psi \left (  \alpha_{i+1/2} - \zeta_{i+1/2} v^n_{x,i+1/2} - \frac{\omega_i v^n_x}{1+\psi \phi_i} \right ) \right ) E^{n+1}_{i+1}  \\
     620- & \left ( \psi \left (  \alpha_{i-1/2} + \zeta_{i-1/2} v^n_{x,i-1/2} + \frac{\omega_i v^n_x}{1+\psi \phi_i} \right ) \right ) E^{n+1}_{i-1}  \\
     621= & \left ( 1 - \bar{\psi} \left ( \alpha_{i+1/2} + \alpha_{i-1/2} + \zeta_{i+1/2} v^n_{x,i+1/2} -  \zeta_{i-1/2} v^n_{x,i-1/2} + \frac{\epsilon_i - \xi_i}{1+\psi \phi_i} \right ) \right ) E^{n}_i  \\
     622+ & \left ( \bar{\psi} \left (  \alpha_{i+1/2} - \zeta_{i+1/2} v^n_{x,i+1/2} - \frac{\omega_i v^n_x}{1+\psi \phi_i} \right ) \right ) E^{n}_{i+1}  \\
     623+ & \left ( \bar{\psi} \left (  \alpha_{i-1/2} + \zeta_{i-1/2} v^n_{x,i-1/2} + \frac{\omega_i v^n_x}{1+\psi \phi_i} \right ) \right ) E^{n}_{i-1}  \\
     624+ & \frac{\theta_i}{1+\psi \phi_i}   \\
    625625\end{eqnarray}
    626626}}}