Upwind scheme

In computational fluid dynamics, the upwind schemes are any of a class of discretization methods to solve hyperbolic partial differential equations numerically. The wave equation, the advection equation, the Euler equations in fluid dynamics, etc. belongs to hyperbolic PDEs. Upwind schemes use an adaptive or solution-sensitive finite difference stencil to numerically simulate more properly the direction of propagation of information in a flow field. More specifically, upwind schemes attempt to discretize hyperbolic partial differential equations by using differencing biased in the direction determined by the sign of the characteristic speeds. Historically, the origin of upwind methods can be traced back to the work of Courant, Isaacson, and Reeves who proposed the CIR method. [Courant, R., Isaacson, E., and Rees, M. (1952). "On the Solution of Non-Linear Hyperbolic Differential Equations", Comm. Pure Appl. Math., 5, 243.]

Model equation

To illustrate the method, consider the following one-dimensional linear wave equation:$qquad\; frac\{partial\; u\}\{partial\; t\}\; +\; a\; frac\{partial\; u\}\{partial\; x\}\; =\; 0$It describes a wave propagating in the $x$-direction with a velocity $a$. The preceding equationis also a mathematical model for one-dimensional linear advection. Consider a typical grid point $i$ in thedomain. In a one-dimensional domain, there are only two direction associated with point $i$ - left andright. If $a$ is positive the left side is called "upwind" side and right side is the "downwind" side. Similarly, if $a$ is negative the left side is called "downwind" side and right side is the "upwind" side. If the finite difference scheme for the spatial derivative, $partial\; u\; /\; partial\; x$contains morepoints 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 [cite book |last=Patankar |first=S. V. |title=Numerical Heat Transfer and Fluid Flow |publisher=Taylor & Francis |year=1980 |isbn=978-0891165224] :$quad\; (1)\; qquad\; frac\{u\_i^\{n+1\}\; -\; u\_i^n\}\{Delta\; t\}\; +\; a\; frac\{u\_i^n\; -\; u\_\{i-1\}^n\}\{Delta\; x\}\; =\; 0\; quad\; ext\{for\}\; quad\; a\; >\; 0$:Defining:$qquad\; qquad\; a^+\; =\; ext\{max\}(a,0),,\; qquad\; a^-\; =\; ext\{min\}(a,0)$and:$qquad\; qquad\; u\_x^-\; =\; frac\{u\_i^\{n\}\; -\; u\_\{i-1\}^\{n\{Delta\; x\},,\; qquad\; u\_x^+\; =\; frac\{u\_\{i+1\}^\{n\}\; -\; u\_\{i\}^\{n\{Delta\; x\}$the two conditional equations (1) and (2) can be combined and written in a compact form as:$quad\; (3)\; qquad\; u\_i^\{n+1\}\; =\; u\_i^n\; -\; Delta\; t\; left\; [\; a^+\; u\_x^-\; +\; a^-\; u\_x^+\; ight]$Equation (3) is a general way of writing any upwind-type schemes.The upwind scheme is stable if the following Courant–Friedrichs–Lewy condition (CFL) condition is satisfied. [cite book |last=Hirsch |first=C.|title=Numerical Computation of Internal and External Flows|year=1990|publisher=John Wiley & Sons |isbn=978-0471924524] :$qquad\; qquad\; c\; =\; left|\; frac\{aDelta\; t\}\{Delta\; x\}\; ight|\; le\; 1\; .$

A Taylor series analysis of the upwind scheme discussed above will show that it is first-order accurate in space and time. The first-order upwind scheme introduces severe numerical diffusion in the solution where large gradients exists.

econd-order upwind scheme

The spatial accuracy of the first-order upwind scheme can be improved by choosing a more accurate finite difference stencil for the approximation of spatial derivative. For the second-order upwind scheme, $u\_x^-$ in equation (3) is defined as :$qquad\; qquad\; u\_x^-\; =\; frac\{3u\_i^n\; -\; 4u\_\{i-1\}^n\; +\; u\_\{i-2\}^n\}\{2Delta\; x\}$and $u\_x^+$ is defined as :$qquad\; qquad\; u\_x^+\; =\; frac\{-u\_\{i+2\}^n\; +\; 4u\_\{i+1\}^n\; -\; 3u\_i^n\}\{2Delta\; x\}$This scheme is less diffusive compared to the first-order accurate scheme.

Third-order upwind scheme

For the third-order upwind scheme, $u\_x^-$ in equation (3) is defined as:$qquad\; qquad\; u\_x^-\; =\; frac\{2u\_\{i+1\}\; +\; 3u\_i\; -\; 6u\_\{i-1\}\; +\; u\_\{i-2\{6Delta\; x\}$and $u\_x^+$ is defined as :$qquad\; qquad\; u\_x^+\; =\; frac\{-u\_\{i+2\}\; +\; 6u\_\{i+1\}\; -\; 3u\_i\; -\; 2u\_\{i-1\{6Delta\; x\}$This scheme is less diffusive compared to the second-order accurate scheme. However, it is known to introduce slight dispersive errors in the region where the gradient is high.

ee also

* Finite volume method

References