A Symmetrical Streamline Stabilization scheme for high advective transport

✍ Scribed by Wendland, E.; Schmid, G.

• Publisher: John Wiley and Sons
• Year: 2000
• Tongue: English
• Weight: 197 KB
• Volume: 24
• Category: Article
• ISSN: 0363-9061
• DOI: 10.1002/(sici)1096-9853(200001)24:1<29::aid-nag51>3.0.co;2-f

✦ Synopsis

An overview of numerical techniques and previous investigations related to the solution of advectiondominated transport processes is presented. In addition a new Symmetrical Streamline Stabilization (S) scheme is introduced. The basis of the technique is to treat the transport equation in two steps. In the "rst step the dispersion part is approximated by a standard Galerkin approach, while in the second step the advection is approximated by a least-squares method. The two parts are reassembled, resulting in one system of equations. The resulting coe$cients' matrix is symmetric. Only half of a sparse matrix needs to be stored. Robust iterative algorithms for symmetrical systems of equations such as the preconditioned conjugate gradient method (PCG) can be successfully used. The new method leads to an implicit introduction of an &arti"cial di!usion' term. Solute transport with high Peclet and Courant numbers does not lead to oscillations due to an inherent upwind damping. The upwind e!ect acts only in #ow direction. The e$ciency of the new formulation in terms of accuracy and computation time is shown in comparison with the Galerkin approach for mesh parallel and mesh oblique high advective solute transport.