Most previous investigations for consolidation-induced solute transport models have been limited to one-dimensional studies in unsaturated porous media and lack systematic parameter sensitivity analysis. This study addresses these gaps by analyzing the effects of hydraulic conductivity (K), shear modulus (G), saturation (Sr), Poisson's ratio (nu), partitioning coefficient (Kd), and anisotropy ratio (KxKz and KyKz) on pore water pressure, soil deformation, and solute transport. The findings reveal that higher Kd values significantly hinder solute migration through enhanced adsorption and reduced vertical transport to deeper layers, while increasing anisotropy ratios primarily enhance horizontal migration, with their effects diminishing beyond a threshold. Additionally, a higher K accelerates pressure dissipation and solute movement, while a lower G increases soil deformation and speeds up solute migration. Saturation has a minor effect on solute concentration, with slight increases under higher Sr. The Poisson ratio significantly impacts the transport of the solute, with smaller nu accelerating and larger nu slowing migration. These insights offer valuable theoretical support for optimizing models in unsaturated porous media.

土体冻结和融化时的水分迁移、相变与传热是一个相互影响的耦合过程。采用基于有限体积法的开源软件Open FOAM,编制描述土体冻融过程的水热耦合计算程序。首先,基于土体水分和热量迁移基本方程、水分相变与温度的平衡方程,同时考虑相变对水分特征参数和热特性参数的影响以及相变潜热对传热过程的影响,建立冻土水热耦合数学模型。然后,采用基于多面体网格的有限体积方法对水热耦合控制方程进行空间离散,采用全隐式向后差分方法对方程进行时间离散,由此编制冻土水热耦合计算程序。该程序具有良好的几何适应性、质量和能量守恒性,具备面向复杂问题的并行计算功能。最后,采用该程序对两组不同温度边界条件的室内土体冻结试验进行数值模拟,并与试验结果进行对比,结果表明该程序可以较为准确地模拟土体冻结过程中温度场和水分场的演化特征。