A Novel Multivariate Spectral Local Quasilinearization Method (MV‐SLQLM) for Modelling Flow, Moisture, Heat, and Solute Transport in Soil

Conventionally, the problem of studying the transport of water, heat, and solute in soil or groundwater systems has been numerically solved using finite difference (FD) or finite element (FE) methods. FE methods are attractive over FD methods because they are geometrically flexible. However, recent studies demonstrate that spectral collocation (SC) methods converge exponentially faster than FD or FE methods using a few grid points or on coarse grids. This work proposes and applies a multivariate spectral local quasilinearization method (MV‐SLQLM) to model the transportation and interaction of soil moisture, heat, and solute concentration in a nonbare soil ridge. The MV‐SLQLM uses a quasilinearization method (QLM) to linearize any nonlinear equations and then employs a local linearization method (LLM) to decouple the linearized system of PDEs to form a sequence of equations that are solved in a computationally efficient manner. The MV‐SLQLM is thus an extension of the bivariate spectral local linearization method (BI‐SLLM) that fails to deal with a 2D problem and is a modification of the MV‐SQLM whose efficiency is compromised when operating on high dense solution matrices. We use the residual error norms of the difference between successive iterations to affirm convergence to the expected solution. To illustrate the application and check the solution accuracy, we conduct systematic analyses of the effect of model parameters on distribution profiles. Findings are in good agreement with theory and literature, thereby revealing suitability of the MV‐SLQLM to solve coupled nonlinear PDEs with environmental fluid dynamics applications.