Backward similarity solution of the Boussinesq groundwater equation
Abstract
Groundwater flow in an unconfined aquifer resting on a horizontal impermeable layer with a boundary condition of a rapid increase in the source water level is considered in this work. The newly introduced condition, referred to as the backward power-law head condition, represents a situation where the water level in the adjacent water body increases more rapidly than do conventional problems, which can only represent a situation akin to a traveling wave under rising water level conditions, given its consideration of infinite time. This problem admits the similarity transformation which allows the nonlinear partial differential equation to be converted into a nonlinear ordinary differential equation via the Boltzmann transformation. The reduced boundary value problem is interpreted as the initial value problem for a system of ordinary differential equations (ODE), which can be numerically solved via Shampine's method. The numerical solutions are in good agreement with Kalashinikov's special solution, which is also introduced into the Boussinesq equation. We find that the solution is consistent with the limit of the forward power-law head condition. The new approximate analytical solution and the associated wetting front position are derived by assuming that the solution has the form of quadratic polynomials, which enables the description of the time progression of a real front position. The obtained approximation is compared to Shampine's solution to check the accuracy. Furthermore, the finite element method is applied to the original partial differential equation (PDE), which validates Shampine's solution for use as a benchmark.