|
Cococubed.com
|
| Su-Olson Problem |
Home
Teaching materials
Astronomy research
Astronomy codes
... Stellar equation of states
... EOS with ionization
... EOS for supernovae
... Chemical potentials
... Stellar atmospheres
... Voigt Function
... Polytropic stars
... Cold white dwarfs
... Hotter white dwarfs
... Cold neutron stars
... Stellar opacities
... Neutrino energy loss rates
... Ephemeris routines
... Fermi-Dirac functions
... Galactic chemical evolution
... Nuclear reaction networks
... Nuclear statistical equilibrium
... Laminar deflagrations
... CJ detonations
... ZND detonations
... Fitting to conic sections
... Unusual linear algebra
... Derivatives on uneven grids
... Pentadiagonal solver
... Quadratics, Cubics, Quartics
... Supernova light curves
... Exact Riemann solutions
... 1D PPM Hydrodynamics
... Verification problems
... EZ stellar evolution
... FLASH code
... Mesa code
Astronomy talks
Astronomy images
Outreach
Family album
Bicycle adventures
Artwork
Contact us:
J.D. Maldonado
F.X.Timmes, my vitae
The Su-Olson problem is a one-dimensional, half-space, non-equilibrium Marshak wave problem. The radiative transfer model is a one-group diffusion approximation with a finite radiation source boundary condition, where the radiative and material fields are out of equilibrium. As the energy density of the radiation field increases, energy is transfered to the material. Su & Olson (1996) found a solution, to quadrature, for the distribution of radiative energy and material temperature as a function of spacetime. This solution is useful for verifying time-dependent radiation diffusion codes. The codes below are released under LA-CC-05-101.
 the three gammas |
 the three thetas |
 various integrands |
|
suo.f reproduces the
paper's tables suo02.f gives physical solutions problem  |
 photon and material temperatures |
 evolution movie |
Radiation temperature uniform mesh:
 Trad 200x200 cells |
 Trad 400x400 cells |
 Trad 800x800 cells |
 Trad error 200x200 cells |
 Trad error 400x400 cells |
 Trad error 800x800 cells |
 Trad asymmetry 200x200 |
 Trad asymmetry 400x400 |
 Trad asymmetry 800x800 |
Material temperature adaptive mesh:
 Tmat 200x200 cells |
 Tmat 400x400 cells |
 Tmat 800x800 cells |
 Tmat error 200x200 cells |
 Tmat error 400x400 cells |
 Tmat error 800x800 cells |
 Tmat asymmetry 200x200 |
 Tmat asymmetry 400x400 |
 Tmat asymmetry 800x800 |