Method of lines

Last updated
Method of lines - the example, which shows the origin of the name of method. Method of lines.gif
Method of lines - the example, which shows the origin of the name of method.

The method of lines (MOL, NMOL, NUMOL [1] [2] [3] ) is a technique for solving partial differential equations (PDEs) in which all but one dimension is discretized. By reducing a PDE to a single continuous dimension, the method of lines allows solutions to be computed via methods and software developed for the numerical integration of ordinary differential equations (ODEs) and differential-algebraic systems of equations (DAEs). Many integration routines have been developed over the years in many different programming languages, and some have been published as open source resources. [4]

Contents

The method of lines most often refers to the construction or analysis of numerical methods for partial differential equations that proceeds by first discretizing the spatial derivatives only and leaving the time variable continuous. This leads to a system of ordinary differential equations to which a numerical method for initial value ordinary equations can be applied. The method of lines in this context dates back to at least the early 1960s. [5] Many papers discussing the accuracy and stability of the method of lines for various types of partial differential equations have appeared since. [6] [7]

Application to elliptical equations

MOL requires that the PDE problem is well-posed as an initial value (Cauchy) problem in at least one dimension, because ODE and DAE integrators are initial value problem (IVP) solvers. Thus it cannot be used directly on purely elliptic partial differential equations, such as Laplace's equation. However, MOL has been used to solve Laplace's equation by using the method of false transients. [1] [8] In this method, a time derivative of the dependent variable is added to Laplace’s equation. Finite differences are then used to approximate the spatial derivatives, and the resulting system of equations is solved by MOL. It is also possible to solve elliptical problems by a semi-analytical method of lines. [9] In this method, the discretization process results in a set of ODE's that are solved by exploiting properties of the associated exponential matrix.

Recently, to overcome the stability issues associated with the method of false transients, a perturbation approach was proposed which was found to be more robust than standard method of false transients for a wide range of elliptic PDEs. [10]

Related Research Articles

<span class="mw-page-title-main">Numerical analysis</span> Study of algorithms using numerical approximation

Numerical analysis is the study of algorithms that use numerical approximation for the problems of mathematical analysis. It is the study of numerical methods that attempt to find approximate solutions of problems rather than the exact ones. Numerical analysis finds application in all fields of engineering and the physical sciences, and in the 21st century also the life and social sciences, medicine, business and even the arts. Current growth in computing power has enabled the use of more complex numerical analysis, providing detailed and realistic mathematical models in science and engineering. Examples of numerical analysis include: ordinary differential equations as found in celestial mechanics, numerical linear algebra in data analysis, and stochastic differential equations and Markov chains for simulating living cells in medicine and biology.

<span class="mw-page-title-main">Partial differential equation</span> Type of differential equation

In mathematics, a partial differential equation (PDE) is an equation which computes a function between various partial derivatives of a multivariable function.

<span class="mw-page-title-main">Numerical methods for ordinary differential equations</span> Methods used to find numerical solutions of ordinary differential equations

Numerical methods for ordinary differential equations are methods used to find numerical approximations to the solutions of ordinary differential equations (ODEs). Their use is also known as "numerical integration", although this term can also refer to the computation of integrals.

<span class="mw-page-title-main">Separation of variables</span> Technique for solving differential equations

In mathematics, separation of variables is any of several methods for solving ordinary and partial differential equations, in which algebra allows one to rewrite an equation so that each of two variables occurs on a different side of the equation.

Second-order linear partial differential equations (PDEs) are classified as either elliptic, hyperbolic, or parabolic. Any second-order linear PDE in two variables can be written in the form

<span class="mw-page-title-main">Differential equation</span> Type of functional equation (mathematics)

In mathematics, a differential equation is an equation that relates one or more unknown functions and their derivatives. In applications, the functions generally represent physical quantities, the derivatives represent their rates of change, and the differential equation defines a relationship between the two. Such relations are common; therefore, differential equations play a prominent role in many disciplines including engineering, physics, economics, and biology.

Numerical methods for partial differential equations is the branch of numerical analysis that studies the numerical solution of partial differential equations (PDEs).

In mathematics, a differential-algebraic system of equations (DAE) is a system of equations that either contains differential equations and algebraic equations, or is equivalent to such a system.

<span class="mw-page-title-main">Computational electromagnetics</span> Branch of physics

Computational electromagnetics (CEM), computational electrodynamics or electromagnetic modeling is the process of modeling the interaction of electromagnetic fields with physical objects and the environment using computers.

In numerical analysis, finite-difference methods (FDM) are a class of numerical techniques for solving differential equations by approximating derivatives with finite differences. Both the spatial domain and time domain are discretized, or broken into a finite number of intervals, and the values of the solution at the end points of the intervals are approximated by solving algebraic equations containing finite differences and values from nearby points.

A parabolic partial differential equation is a type of partial differential equation (PDE). Parabolic PDEs are used to describe a wide variety of time-dependent phenomena, including heat conduction, particle diffusion, and pricing of derivative investment instruments.

Dynamic simulation is the use of a computer program to model the time-varying behavior of a dynamical system. The systems are typically described by ordinary differential equations or partial differential equations. A simulation run solves the state-equation system to find the behavior of the state variables over a specified period of time. The equation is solved through numerical integration methods to produce the transient behavior of the state variables. Simulation of dynamic systems predicts the values of model-system state variables, as they are determined by the past state values. This relationship is found by creating a model of the system.

<span class="mw-page-title-main">Finite element method</span> Numerical method for solving physical or engineering problems

The finite element method (FEM) is a popular method for numerically solving differential equations arising in engineering and mathematical modeling. Typical problem areas of interest include the traditional fields of structural analysis, heat transfer, fluid flow, mass transport, and electromagnetic potential.

Finite difference methods for option pricing are numerical methods used in mathematical finance for the valuation of options. Finite difference methods were first applied to option pricing by Eduardo Schwartz in 1977.

<span class="mw-page-title-main">Ordinary differential equation</span> Differential equation containing derivatives with respect to only one variable

In mathematics, an ordinary differential equation (ODE) is a differential equation (DE) dependent on only a single independent variable. As with other DE, its unknown(s) consists of one function(s) and involves the derivatives of those functions. The term "ordinary" is used in contrast with partial differential equations (PDEs) which may be with respect to more than one independent variable, and, less commonly, in contrast with stochastic differential equations (SDEs) where the progression is random.

In mathematics a partial differential algebraic equation (PDAE) set is an incomplete system of partial differential equations that is closed with a set of algebraic equations.

<span class="mw-page-title-main">Finite water-content vadose zone flow method</span>

The finite water-content vadose zone flux method represents a one-dimensional alternative to the numerical solution of Richards' equation for simulating the movement of water in unsaturated soils. The finite water-content method solves the advection-like term of the Soil Moisture Velocity Equation, which is an ordinary differential equation alternative to the Richards partial differential equation. The Richards equation is difficult to approximate in general because it does not have a closed-form analytical solution except in a few cases. The finite water-content method, is perhaps the first generic replacement for the numerical solution of the Richards' equation. The finite water-content solution has several advantages over the Richards equation solution. First, as an ordinary differential equation it is explicit, guaranteed to converge and computationally inexpensive to solve. Second, using a finite volume solution methodology it is guaranteed to conserve mass. The finite water content method readily simulates sharp wetting fronts, something that the Richards solution struggles with. The main limiting assumption required to use the finite water-content method is that the soil be homogeneous in layers.

The soil moisture velocity equation describes the speed that water moves vertically through unsaturated soil under the combined actions of gravity and capillarity, a process known as infiltration. The equation is alternative form of the Richardson/Richards' equation. The key difference being that the dependent variable is the position of the wetting front , which is a function of time, the water content and media properties. The soil moisture velocity equation consists of two terms. The first "advection-like" term was developed to simulate surface infiltration and was extended to the water table, which was verified using data collected in a column experimental that was patterned after the famous experiment by Childs & Poulovassilis (1962) and against exact solutions.

Validated numerics, or rigorous computation, verified computation, reliable computation, numerical verification is numerics including mathematically strict error evaluation, and it is one field of numerical analysis. For computation, interval arithmetic is used, and all results are represented by intervals. Validated numerics were used by Warwick Tucker in order to solve the 14th of Smale's problems, and today it is recognized as a powerful tool for the study of dynamical systems.

References

  1. 1 2 Schiesser, W. E. (1991). The Numerical Method of Lines . Academic Press. ISBN   0-12-624130-9.
  2. ↑ Hamdi, S.; W. E. Schiesser; G. W. Griffiths (2007), "Method of lines", Scholarpedia, 2 (7): 2859, Bibcode:2007SchpJ...2.2859H, doi: 10.4249/scholarpedia.2859
  3. ↑ Schiesser, W. E.; G. W. Griffiths (2009). A Compendium of Partial Differential Equation Models: Method of Lines Analysis with Matlab. Cambridge University Press. ISBN   978-0-521-51986-1.
  4. ↑ Lee, H. J.; W. E. Schiesser (2004). Ordinary and Partial Differential Equation Routines in C, C++, Fortran, Java, Maple and Matlab. CRC Press. ISBN   1-58488-423-1.
  5. ↑ E. N. Sarmin; L. A. Chudov (1963), "On the stability of the numerical integration of systems of ordinary differential equations arising in the use of the straight line method", USSR Computational Mathematics and Mathematical Physics, 3 (6): 1537–1543, doi:10.1016/0041-5553(63)90256-8
  6. ↑ A. Zafarullah (1970), "Application of the Method of Lines to Parabolic Partial Differential Equations With Error Estimates", Journal of the Association for Computing Machinery, vol. 17, no. 2, pp. 294–302, doi: 10.1145/321574.321583 , S2CID   15114435
  7. ↑ J. G. Verwer; J. M. Sanz-Serna (1984), "Convergence of method of lines approximations to partial differential equations", Computing, 33 (3–4): 297–313, doi:10.1007/bf02242274, S2CID   30171258
  8. ↑ Schiesser, W. E. (1994). Computational mathematics in Engineering and Applied Science: ODEs, DAEs and PDEs. CRC Press. ISBN   0-8493-7373-5.
  9. ↑ Subramanian, V.R.; R.E. White (2004), "Semianalytical method of lines for solving elliptic partial differential equations", Chemical Engineering Science, 59 (4): 781–788, Bibcode:2004ChEnS..59..781S, doi:10.1016/j.ces.2003.10.019
  10. ↑ P. W. C. Northrop; P. A. Ramachandran; W. E. Schiesser; V. R. Subramanian (2013), "A Robust False Transient Method of Lines for Elliptic Partial Differential Equations", Chem. Eng. Sci., vol. 90, pp. 32–39, Bibcode:2013ChEnS..90...32N, doi:10.1016/j.ces.2012.11.033


This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.