“…In this article, we add the Soret effect (thermal diffusion), discuss other strategies to solve the stiff system of species equations, and compare the quality of the solutions obtained by the different strategies. As in [3], it is assumed that the velocities, temperatures, pressures, and densities are known in the computational domain. Our work in [3] and this article, differ from work in literature on (laminar) reacting flows, in which, as a result of computational costs and stiffness, the simulations are limited to steady state and mostly 2D (see [1]).…”