Difference scheme

A system of difference equations (cf. Difference equation) approximating a differential equation and subsidiary (initial, boundary and other) conditions. Approximating an original differential problem by a difference scheme is one of the methods of approximating the original problem by a discrete problem. For this the given region $ G $ of independent variables is replaced by a discrete set of points $ G _ {h} $— the grid, and the derivatives entering into the differential equation are replaced by difference relations on the grid $ G _ {h} $. As a result of such a change one obtains a closed system of a large number of algebraic equations (linear or non-linear, depending on the original differential equation) which makes up the difference scheme. In essence a difference scheme is a family of difference equations depending on the steps of the grid. The solution of a difference scheme also depends parametrically on the steps of the grid. A difference scheme is a multi-parameter and complicated object. Besides the coefficients of the original differential equation it contains also its own characteristic parameters, such as the steps with respect to time and space, weight factors, and others. The influence of these parameters can substantially distort the representation of the behaviour of the original differential problem.

The following questions are being studied in connection with the approximation by difference schemes of differential equations: means of constructing a difference scheme; the convergence under grid refinement of the solution of the difference problem to the solution of the original differential equation; and methods for solving systems of difference equations. All the questions listed are considered in the theory of difference schemes (cf. Difference schemes, theory of). Effective numerical methods for solving typical difference schemes for ordinary and partial differential equations, assuming the use of a high-speed computer, have been developed.

A simple example of a difference scheme is given below. Suppose one is given the differential equation

$$ \tag{1 } \left . \begin{array}{c} u ^ {\prime\prime} ( x) - q ( x) u ( x) = - f ( x),\ \ q ( x) \geq 0,\ \ 0 < x < 1, \\ u ^ \prime ( 0) = \ \sigma u ( 0) - \mu _ {1} ,\ \ u ( 1) = \ \mu _ {2} ,\ \ \sigma > 0. \\ \end{array} \right \} $$

The domain $ G = \{ 0 < x < 1 \} $ is replaced by the grid

$$ G _ {h} = \ \{ {x _ {i} = ih } : { i = 0 \dots N; hN = 1 } \} . $$

A difference scheme for problem (1) has the form

$$ \tag{2 } \left . \begin{array}{c} \frac{y _ {i + 1 } - 2y _ {i} + y _ {i - 1 } }{h ^ {2} } - q _ {i} y _ {i} = - f _ {i} ,\ \ i = 1 \dots N - 1, \\ { \frac{y _ {1} - y _ {0} }{h} } = \ ( \sigma + 0.5hq _ {0} ) y _ {0} - ( \mu _ {1} + 0.5hf _ {0} ) y _ {N} = \mu _ {2} , \\ \end{array} \right \} $$

where $ y _ {i} = y ( x _ {i} ) $, $ q _ {i} = q ( x _ {i} ) $, $ x _ {i} \in G _ {h} $. It can be shown that, as $ h \rightarrow 0 $, the solution to the difference problem (2) converges to the solution of the original problem (1) for sufficiently smooth functions $ q , f $.

The difference scheme (2) has second-order accuracy, that is,

$$ \max _ {0 \leq i \leq N } \ | y _ {i} - u ( x _ {i} ) | \leq \ Mh ^ {2} , $$

where $ M $ is a constant that does not depend on $ h $. The solution to the difference scheme (2) can be found by the shooting method.

References

[1]	A.A. Samarskii, "Theorie der Differenzverfahren" , Akad. Verlagsgesell. Geest u. Portig K.-D. (1984) (Translated from Russian)
[2]	A.A. Samarskii, E.S. Nikolaev, "Numerical methods for grid equations" , 1–2 , Birkhäuser (1989) (Translated from Russian)

References

[a1]	G.E. Forsythe, W.R. Wasow, "Finite difference methods for partial differential equations" , Wiley (1960)
[a2]	P.R. Garabedian, "Partial differential equations" , Wiley (1964)
[a3]	I. Gladwell (ed.) R. Wait (ed.) , A survey of numerical methods for partial differential equations , Clarendon Press (1979)
[a4]	A.R. Mitchell, D.F. Griffiths, "The finite difference method in partial differential equations" , Wiley (1980)
[a5]	R.D. Richtmeyer, K.W. Morton, "Difference methods for initial value problems" , Wiley (1967)
[a6]	G.D. Smith, "Numerical solution of partial differential equations" , Oxford Univ. Press (1977)
[a7]	N.N. Yanenko, "The method of fractional steps: solution of problems of mathematical physics in several variables" , Springer (1971) (Translated from Russian)