z-logo
Premium
The Laplace Transform Finite Difference Method for Simulation of Flow Through Porous Media
Author(s) -
Moridis George J.,
Reddell Donald L.
Publication year - 1991
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/91wr01190
Subject(s) - discretization , laplace transform , finite difference , finite difference method , matrix (chemical analysis) , partial differential equation , mathematics , transient (computer programming) , time derivative , stability (learning theory) , mathematical analysis , laplace transform applied to differential equations , numerical stability , flow (mathematics) , porous medium , mathematical optimization , numerical analysis , computer science , porosity , geometry , materials science , machine learning , composite material , operating system
A new numerical method, the Laplace transform finite difference (LTFD) method, was developed to solve the partial differential equation (PDE) of transient flow through porous media. LTFD provides a solution which is semianalytical in time and numerical in space by solving the discretized PDE in the Laplace space and numerically inverting the transformed solution vector. The effects of the traditional treatment of the time derivative on accuracy and stability are rendered irrelevant because time is no longer a consideration. For a single time step, LTFD requires no more than eight matrix solutions and an execution time eight times longer than the analogous finite difference (FD) requirement without an increase in storage. This disadvantage is outweighed by an unlimited time step size without any loss of accuracy, a superior accuracy, and a stable, nonincreasing round off error. Thus, a problem in standard FD format may require several hundred time steps and matrix inversions between the initial condition and the desired solution time, but LTFD requires only one time step and no more than eight matrix inversions to achieve a more accurate result.

This content is not available in your region!

Continue researching here.

Having issues? You can contact us here