Results

• Morel Asymmetric Method
  • Decreasing intensity (like an intensity sink) starts when wave reaches skewed part of the mesh
  • Fix-up: when radiation temperature dips below the temperature floor, use low-order scheme in that cell
  • Fix-up eliminates positive off-diagonals in matrix, which would guarantee a positive solution if done over entire mesh
  • Fix-up was successful: problem runs until steady-state
  • All plots are from this method
• Support Operator Method
  • Instabilities:
    • grow without bound from roundoff values
    • located at the skewed parts of the mesh
    • begin at t $\approx$ 6 x 10^-7  sh, before the wave has reached the area
    • could be a coding error?
  • Fix-up has no effect