Changes between Version 28 and Version 29 of u/bliu/pnfldiff
- Timestamp:
- 08/16/18 16:50:23 (6 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
u/bliu/pnfldiff
v28 v29 69 69 Discretize the equation of $n_{e}$ in the same way as [https://astrobear.pas.rochester.edu/trac/wiki/ThermalConduction Thermal Conduction] 70 70 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$ 72 72 73 73 Replace $n_{e}$ with $n_{e}^{*}$ where[[BR]] … … 76 76 $\Phi=1$ for backward Euler, $\Phi=1/2$ for Crank Nicholson 77 77 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$ 79 79 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}$ 82 85 83 86 == Links and References ==