Changes between Version 5 and Version 6 of FluxLimitedDiffusion


Ignore:
Timestamp:
03/18/13 12:54:07 (12 years ago)
Author:
Jonathan
Comment:

Legend:

Unmodified
Added
Removed
Modified
  • FluxLimitedDiffusion

    v5 v6  
    7777For static diffusion, the terms in blue with v^2^/c can be dropped and the system can be split into the usual hydro update (black), radiative source terms (green) using time centered radiation energy, and a coupled implicit solve (red) for the radiation energy density and thermal energy density (ie temperature).  If the opacity is not temperature dependent, then the implicit solve only involves the radiation energy density.  If the opacity is not a linear function of the energy, then some form of sub-cycling would be required.
    7878
     79For now we will assume that [[latex(\kappa_{0P})]] and [[latex(\kappa_{0R})]] are constant over the implicit update.  In this case we can solve the radiation energy equation:
     80
     81[[latex(\frac{\partial E}{\partial t}  = \nabla \cdot \frac{c\lambda}{\kappa_{0R}} \nabla E + \kappa_{0P} (4 \pi B-cE))]]
     82
     83
     84