Changes between Version 28 and Version 29 of u/bliu/pnfldiff


Ignore:
Timestamp:
08/16/18 16:50:23 (6 years ago)
Author:
Baowei Liu
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • u/bliu/pnfldiff

    v28 v29  
    6969Discretize the equation of $n_{e}$ in the same way as [https://astrobear.pas.rochester.edu/trac/wiki/ThermalConduction Thermal Conduction]
    7070
    71 $\partial_{t}n_{e} +\partial_{i}(n_{e}u_{i})=I-r$
     71$\partial_{t}n_{e} +\partial_{i}(n_{e}u_{i})=\partial_{i}J_{i}-r$
    7272
    7373Replace $n_{e}$ with $n_{e}^{*}$ where[[BR]]
     
    7676$\Phi=1$ for backward Euler, $\Phi=1/2$ for Crank Nicholson
    7777
    78 $\partial_{t}n_{e} +u_{i}\partial_{i}n_{e}^{*}+n_{e}^{*}\partial_{i}u_{i}=I-r$
     78$\partial_{t}n_{e} +u_{i}\partial_{i}n_{e}^{*}+n_{e}^{*}\partial_{i}u_{i}=\partial_{i}J_{i}-r$
    7979
    80 
    81 
     80$\partial_{t}n_{e}=\frac{n_{e}^{\prime\hat{i}}-n_{e}^{\hat{i}}}{\Delta t} $
     81$-(1-\Phi)u_{i}\partial_{i}n_{e}=(1-\Phi)u_{\hat{i}}\frac{n^{\hat{i}+1}_{e}-n^{\hat{i}-1}_{e}}{2\Delta x}$
     82$-\Phi u_{i}\partial_{i}n^{'}_{e}=-\Phi u_{\hat{i}}\frac{n^{'\hat{i}+1}_{e}-n^{'\hat{i}-1}_{e}}{2\Delta x}$
     83$-(1-\Phi)n_{e}\partial_{i}u_{i}=(1-\Phi)n_{e}^{\hat{i}}\frac{u_{\hat{i}+1}-u_{\hat{i}-1}}{2\Delta x}$
     84$-\Phi n_{e}^{'}\partial_{i}u_{i}=-\Phi\frac{u_{\hat{i}+1}-u_{\hat{i}-1}}{2\Delta x}$
    8285
    8386== Links and References ==