ABSTRACT
ABSTRACT: This work deals with the efficient numerical solution of nonlinear transient flow problems posed on a two-dimensional soil profile of general geometry. The spatial semidiscretization of such problems is achieved by using a cell-centered finite difference scheme on a logically rectangular grid. The resulting nonlinear stiff initial-value problem is then integrated in time by means of a fractional step method, combined with a decomposition of the flow domain into a set of overlapping subdomains. The totally discrete scheme involves a set of nonlinear systems which can be linearized by means of suitable Taylor expansions. Thus, the proposed algorithm reduces the original problem to the solution of several linear systems per time step. Moreover, each one of such systems is a set of uncoupled linear subsystems which can be solved in parallel. A numerical example illustrates the behaviour of the method in the last section of the paper.
