AbstractA numerical model, called HT2, is presented for one-dimensional heat transfer in saturated incompressible layered soil with steady fluid flow, effective porosity, and possible local thermal nonequilibrium (LTNE) between solid and fluid phases. The model uses a series-parallel approach and accounts for advection, conduction, linear and nonlinear thermal dispersion, and interstitial heat transfer. The key to HT2 is the definition of separate columns for the solid matrix (SM) and mobile pore fluid (MPF). The SM column includes the solid phase and immobile pore fluid and consists of fixed elements. The MPF column consists of moving elements and uses Lagrangian element-tracking to follow the fluid motion, which reduces numerical dispersion and simplifies heat transfer within the column to that of dispersive flux between contiguous elements. The model uses an interstitial heat transfer coefficient to calculate kinetic heat transfer between SM and MPF elements. The development of HT2 is first presented, followed by verification checks and a limited parametric study. Numerical simulations indicate that larger particle size and higher fluid discharge velocity yields greater local temperature disequilibrium between solid and fluid phases during heat transfer through saturated soil.