80
Views
2
CrossRef citations to date
0
Altmetric
Original Articles

An Alternative Direction Iterative Method With Second-Order Upwind Scheme For Convection-Diffusion Equations

Pages 527-533 | Published online: 15 Sep 2010
 

Consider the following convection diffusion equation, $$\left\{\matrix{{\displaystyle{\partial u \over \partial t}} + b_{1}(x, y){\displaystyle{\partial u \over \partial x}} + b_{2}(y){\displaystyle{\partial u \over \partial y}} - \left(a_{1} {\displaystyle{\partial^{2}u \over \partial x^{2}}} + a_{2} {\displaystyle{\partial^{2}u \over \partial y^{2}}}\right) = f\hfill \hbox{in } \Omega \times J,\hfill\cr u(x, y, t) = \varphi (x, t)\hfill \hbox{on } \partial \Omega \times J,\hfill\cr u(x, y, 0) = u_{0}(x, y)\hfill \hbox{in } \Omega,\hfill\cr}\right. \eqno(1)$$ where z =(0, 1) 2 (0, 1), J =(0, T ), b 1 ( x , y ), b 2 ( y ) are smooth functions and a 1 , a 2 are positive constants. When convection dominates diffusion, i.e. 0< a 1 , a 2 $\ll$ | b |, the general finite difference or finite element methods often result in numerical oscillation [1]. The upwind method is an efficient method but is only first order accurate. For one dimensional stable problem with constant coefficient, [3] presented a high order upwind scheme, but it is difficult to extend the method to variable coefficient problem and two dimensional problem. In this paper we give an alternative direction iterative method combining with one dimensional second order upwind scheme for two dimensional problem. It can be as a high speed algorithm on parallel computer. The maximum principle and the second order uniform norm error estimate are obtained. Finally we give some numerical examples.

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.