Picard-Newton Iterative Method with Time Step Control for Multimaterial Non-Equilibrium Radiation Diffusion Problem. Commun. Comput. Phys. 2011;10(4):844-866. doi:10.4208/cicp.310110.161010a