Nonlinear convergence, accuracy, and time step control in nonequilibrium radiation diffusion

Abstract We study the interaction between converging the nonlinearities within a time step and time step control, on the accuracy of nonequilibrium radiation diffusion calculations. Typically, this type of calculation is performed using operator-splitting where the nonlinearities are lagged one time step. This method of integrating the nonlinear system results in an “effective” time-step constraint to obtain accuracy. A time-step control that limits the change in dependent variables (usually energy) per time step is used. We investigate the possibility that converging the nonlinearities within a time step may allow significantly larger time-step sizes and improved accuracy as well. The previously described Jacobian-free Newton–Krylov method (JQSRT 63 (1999) 15) is used to converge all nonlinearities within a time step. In addition, a new time-step control method, based on the hyperbolic model of a thermal wave (J. Comput. Phys. 152 (1999) 790), is employed. The benefits and cost of a second-order accurate time step are considered. It is demonstrated that for a chosen accuracy, significant increases in solution efficiency can be obtained by converging nonlinearities within a time step.