... the non-constant transport properties are treated by lagging in time leading to decoupled diagonally dominant systems. The consistency, stability and convergence results are derived. The method of manufactured solutions is adopted to numerically verify the theoretical results. The scheme is then applied...

... of the dependent variable in the parabolic PDE for the evaluation of the time derivative. As will be demonstrated later, this key step is responsible for improving the accuracy of ITMOL over its predecessor TMOL. Throughout the paper, the theoretical properties of ITMOL, such as consistency, stability, convergence...