Compact Stencil

In mathematics, especially in the areas of numerical analysis called numerical partial differential equations, a compact stencil is a type of stencil that uses only nine nodes for its discretization method in two dimensions. It uses only the center node and the adjacent nodes. For any structured grid utilizing a compact stencil in 1, 2, or 3 dimensions the maximum number of nodes is 3, 9, or 27 respectively. Compact stencils may be compared to non-compact stencils. Compact stencils are currently implemented in many partial differential equation solvers, including several in the topics of CFD, FEA, and other mathematical solvers relating to PDE's.Communications in Numerical Methods in Engineering, Copyright © 2008 John Wiley & Sons, Ltd.

Two Point Stencil Example

The two point stencil for the ''first derivative'' of a function is given by:

$f'(x_0)=\frac + O\left(h^2\right)$.

This is obtained from the Taylor series expansion of the first derivative of the function given by:

$\begin f'(x_0)=\frac -\frach - \frach^2 - \frach^3 + \cdots \end$.

Replacing $h$ with $-h$, we have:

$\begin f'(x_0)=-\frac + \frach - \frach^2 + \frach^3 + \cdots \end$.

Addition of the above two equations together results in the cancellation of the terms in odd powers of $h$:

$\begin 2f'(x_0)= \frac -\frac -2\frach^2 + \cdots \end$.

$\begin f'(x_0)= \frac - \frach^2 + \cdots \end$.

$\begin f'(x_0)= \frac + O\left(h^2\right) \end$.

Three Point Stencil Example

For example, the three point stencil for the ''second derivative'' of a function is given by:

$\begin f^(x_0)= \frac + O\left(h^2\right) \end$.

This is obtained from the Taylor series expansion of the first derivative of the function given by:

$\begin f'(x_0)=\frac -\frach - \frach^2 - \frach^3 + \cdots \end$.

Replacing $h$ with $-h$, we have:

$\begin f'(x_0)=-\frac + \frach - \frach^2 + \frach^3 + \cdots \end$.

Subtraction of the above two equations results in the cancellation of the terms in even powers of $h$:
$\begin 0= \frac +\frac - 2\frach - 2\frach^3 + \cdots \end$.

$\begin f^(x_0)= \frac - 2\frach^2 + \cdots \end$.

$\begin f^(x_0)= \frac + O\left(h^2\right) \end$.

See also

*Stencil (numerical analysis)
*Non-compact stencil
*Five-point stencil

References

{{reflist

Numerical differential equations