Context Navigation

Changes between Version 24 and Version 25 of FluxLimitedDiffusion

-              v24
+              v25
+If we ignore the change in the Planck function due to heating during the implicit solve, it is equivalent to replacing the term with [[latex(\psi \phi^n_i e^{n+1}_i)]] with [[latex(\psi \phi^n_i e^n_i)]] in which case the equations simplify to
+||   [[latex(\left [ 1 + \psi \left( \alpha^n_{i+1/2} + \alpha^n_{i-1/2} + \epsilon^n_i \right ) \right ] E^{n+1}_i - \left ( \psi \alpha^n_{i+1/2} \right ) E^{n+1}_{i+1} - \left ( \psi \alpha^n_{i-1/2} \right ) E^{n+1}_{i-1} =\left [ 1 - \bar{\psi} \left( \alpha^n_{i+1/2} + \alpha^n_{i-1/2} + \epsilon^n_i \right ) \right ] E^n_i+\phi^n_i e^n_i + \theta^n_i)]]   ||
+||   [[latex(\left (1 \right ) e^{n+1}_i = \left ( \psi \epsilon^n_i \right )E^{n+1}_i + \left ( 1 - \phi^n_i \right ) e^n_i + \left ( \bar{\psi} \epsilon^n_i \right ) E^n_i-\theta^i_n )]]   ||
 Clearly the second equation is trivial to solve after the first system of equations has been solved.  So if we treat the temperature as being constant we can calculate the local heating/cooling due to radiative emission/absorption