In computational physics, the term advection scheme refers to a class of numerical discretization methods for solving hyperbolic partial differential equations. In the so-called upwind schemes typically, the so-called upstream variables are used to calculate the derivatives in a flow field. That is, derivatives are estimated using a set of data points biased to be more "upwind" of the query point, with respect to the direction of the flow. Historically, the origin of upwind methods can be traced back to the work of Courant, Isaacson, and Rees who proposed the CIR method.
Model equation To illustrate the method, consider the following one-dimensional linear advection equation
∂ u ∂ t + a ∂ u ∂ x = 0 {\displaystyle {\frac {\partial u}{\partial t}}+a{\frac {\partial u}{\partial x}}=0}
which describes a wave propagating along the x {\displaystyle x} -axis with a velocity a {\displaystyle a} . This equation is also a mathematical model for one-dimensional linear advection. Consider a typical grid point i {\displaystyle i} in the domain. In a one-dimensional domain, there are only two directions associated with point i {\displaystyle i} – left (towards negative infinity) and right (towards positive infinity). If a {\displaystyle a} is positive, the traveling wave solution of the equation above propagates towards the right, the left side is called the upwind side and the right side is the downwind side. Similarly, if a {\displaystyle a} is negative the traveling wave solution propagates towards the left, the left side is called downwind side and right side is the upwind side. If the finite difference scheme for the spatial derivative, ∂ u / ∂ x {\displaystyle \partial u/\partial x} contains more points in the upwind side, the scheme is called an upwind-biased or simply an upwind scheme.
First-order upwind scheme
The simplest upwind scheme possible is the first-order upwind scheme. It is given by
where n {\displaystyle n} refers to the t {\displaystyle t} dimension and i {\displaystyle i} refers to the x {\displaystyle x} dimension. (By comparison, a central difference scheme in this scenario would look like
u i n + 1 − u i n Δ t + a u i + 1 n − u i − 1 n 2 Δ x = 0 , {\displaystyle {\frac {u_{i}^{n+1}-u_{i}^{n}}{\Delta t}}+a{\frac {u_{i+1}^{n}-u_{i-1}^{n}}{2\Delta x}}=0,}
regardless of the sign of a {\displaystyle a} .)
Compact form Defining
… excerpt ends here. Continue reading the full article.


