Stabilized finite element methods for time dependent convection-diffusion equations

STABILIZED FINITE ELEMENT METHODS FOR TIME DEPENDENT CONVECTION-DIFFUSION EQUATIONS In this thesis, enriched finite element methods are presented for both steady and unsteady convection diffusion equations. For the unsteady case, we follow the method of lines approach that consists of first discretizing in space and then use some time integrator to solve the resulting system of ordinary differential equation. Discretization in time is performed by the generalized Euler finite difference scheme, while for the space discretization the streamline upwind Petrov-Galerkin (SUPG), the Residual free bubble (RFB), the more recent multiscale (MS) and specific combination of RFB with MS (MIX) methods are considered. To apply the RFB and the MS methods, the steady local problem, which is as complicated as the original steady equation, should be solved in each element. That requirement makes these methods quite expensive especially for two dimensional problems. In order to overcome that drawback the pseudo approximation techniques, which employ only a few nodes in each element, are used. Next, for the unsteady problem a proper adaptation recipe, including these approximations combined with the generalized Euler time discretization, is described. For piecewise linear finite element discretization on triangular grid, the SUPG method is used. Then we derive an efficient stability parameter by examining the relation of the RFB and the SUPG methods. Stability and convergence analysis of the SUPG method applied to the unsteady problem is obtained by extending the Burman’s analysis techniques for the pure convection problem. We also suggest a novel operator splitting strategy for the transport equations with nonlinear reaction term. As a result two subproblems are obtained. One of which we may apply using the SUPG stabilization while the other equation can be solved analytically. Lastly, numerical experiments are presented to illustrate the good performance of the method.

[1]  J. M. Thomas,et al.  Introduction à l'analyse numérique des équations aux dérivées partielles , 1983 .

[2]  J. Douglas,et al.  Stabilized mixed methods for the Stokes problem , 1988 .

[3]  A. Quarteroni,et al.  Numerical Approximation of Partial Differential Equations , 2008 .

[4]  A. I. Nesliturk,et al.  On the choice of stabilizing sub-grid for convection-diffusion problem on rectangular grids , 2010, Comput. Math. Appl..

[5]  Isaac Harari,et al.  Stability of semidiscrete formulations for parabolic problems at small time steps , 2004 .

[6]  P. Frolkovic,et al.  Numerical Simulations Of Radionuclides Transport In Double Porosity Media With Sorption , 2000 .

[7]  Pavel B. Bochev,et al.  Stability of the SUPG finite element method for transient advection-diffusion problems , 2004 .

[8]  Alessandro Russo,et al.  Further considerations on residual-free bubbles for advective - diffusive equations , 1998 .

[9]  A. Nesliturk A STABILIZING SUBGRID FOR CONVECTION¿DIFFUSION PROBLEM , 2006 .

[10]  Javier de Frutos,et al.  Stabilization of Galerkin Finite Element Approximations to Transient Convection-Diffusion Problems , 2010, SIAM J. Numer. Anal..

[11]  Giancarlo Sangalli,et al.  LINK-CUTTING BUBBLES FOR THE STABILIZATION OF CONVECTION-DIFFUSION-REACTION PROBLEMS , 2003 .

[12]  Gert Lube,et al.  Stabilized finite element methods for singularly perturbed parabolic problems , 1995 .

[13]  Josef Stoer,et al.  Numerische Mathematik 1 , 1989 .

[14]  Javier de Frutos,et al.  Bubble stabilization of linear finite element methods for nonlinear evolutionary convection–diffusion equations , 2008 .

[15]  A. Nesliturk,et al.  Pseudo-multi-scale functions for the stabilization of convection-diffusion equations on rectangular grids , 2013 .

[16]  Volker John,et al.  Error Analysis of the SUPG Finite Element Discretization of Evolutionary Convection-Diffusion-Reaction Equations , 2011, SIAM J. Numer. Anal..

[17]  V. Thomée Galerkin Finite Element Methods for Parabolic Problems (Springer Series in Computational Mathematics) , 2010 .

[18]  T. Hughes,et al.  Stabilized finite element methods. I: Application to the advective-diffusive model , 1992 .

[19]  G. Marchuk Splitting and alternating direction methods , 1990 .

[20]  Alessandro Russo,et al.  Deriving upwinding, mass lumping and selective reduced integration by residual-free bubbles , 1996 .

[21]  B. Sleeman,et al.  Mathematical modeling of capillary formation and development in tumor angiogenesis: Penetration into the stroma , 2001, Bulletin of mathematical biology.

[22]  G. Tanoglu,et al.  An Operator Splitting Approximation Combined With The SUPG Method For Transport Equations With Nonlinear Reaction Term , 2012 .

[23]  J. Geiser Decomposition Methods for Differential Equations: Theory and Applications , 2009 .

[24]  李幼升,et al.  Ph , 1989 .

[25]  Claes Johnson,et al.  Finite element methods for linear hyperbolic problems , 1984 .

[26]  L. Franca,et al.  Stabilized finite element methods. II: The incompressible Navier-Stokes equations , 1992 .

[27]  T. Hughes,et al.  A new finite element formulation for computational fluid dynamics: II. Beyond SUPG , 1986 .

[28]  L. Franca,et al.  On an Improved Unusual Stabilized Finite Element Method for theAdvective-Reactive-Diffusive Equation , 1999 .

[29]  L. Franca,et al.  Stabilized Finite Element Methods , 1993 .

[30]  Alessandro Russo,et al.  CHOOSING BUBBLES FOR ADVECTION-DIFFUSION PROBLEMS , 1994 .

[31]  Erik Burman,et al.  Consistent SUPG-method for transient transport problems: Stability and convergence , 2010 .

[32]  Frédéric Valentin,et al.  Towards multiscale functions : enriching finite element spaces with local but not bubble-like functions , 2005 .

[33]  T. Hughes,et al.  Streamline upwind/Petrov-Galerkin formulations for convection dominated flows with particular emphasis on the incompressible Navier-Stokes equations , 1990 .

[34]  F. Brezzi,et al.  A relationship between stabilized finite element methods and the Galerkin method with bubble functions , 1992 .

[35]  Alessandro Russo,et al.  On the choice of a stabilizing subgrid for convection?diffusion problems , 2005 .

[36]  L. D. Marini,et al.  A Priori Error Analysis of Residual-Free Bubbles for Advection-Diffusion Problems , 1999 .

[37]  Ismael Herrera,et al.  Contaminant transport and biodegradation: 1. A numerical model for reactive transport in porous media , 1989 .

[38]  R. Rannacher,et al.  Finite-element approximations of the nonstationary Navier-Stokes problem. Part IV: error estimates for second-order time discretization , 1990 .

[39]  Frédéric Valentin,et al.  Convergence Analysis of a Multiscale Finite Element Method for Singularly Perturbed Problems , 2005, Multiscale Model. Simul..

[40]  Giancarlo Sangalli,et al.  Global and Local Error Analysis for the Residual-Free Bubbles Method Applied to Advection-Dominated Problems , 2000, SIAM J. Numer. Anal..

[41]  Wulf G. Dettmer,et al.  An analysis of the time integration algorithms for the finite element solutions of incompressible Navier-Stokes equations based on a stabilised formulation , 2003 .

[42]  Thomas J. R. Hughes,et al.  Stabilized finite element methods for steady advection-diffusion with production , 1994 .

[43]  Alessandro Russo,et al.  Applications of the pseudo residual-free bubbles to the stabilization of convection-diffusion problems , 1998 .

[44]  Erik Burman,et al.  ANALYSIS OF THE SPACE SEMI-DISCRETIZED SUPG METHOD FOR TRANSIENT CONVECTION-DIFFUSION EQUATIONS , 2011 .

[45]  Frédéric Valentin,et al.  Enriched finite element methods for unsteady reaction-diffusion problems , 2005 .

[46]  Frédéric Valentin,et al.  Multiscale and Residual-Free Bubble Functions for Reaction-Advection-Diffusion Problems , 2005 .

[47]  B. Sleeman,et al.  Mathematical modeling of capillary formation and development in tumor angiogenesis: Penetration into the stroma , 2001 .

[48]  Giancarlo Sangalli,et al.  Coupling stabilized finite element methods with finite difference time integration for advection–diffusion–reaction problems , 2007 .

[49]  Martin Stynes,et al.  On the Stability of Residual-Free Bubbles for Convection-Diffusion Problemsand their Approximation by a Two-Level Finite Element Method. , 1997 .

[50]  Zahari Zlatev,et al.  Computer Treatment of Large Air Pollution Models , 1995 .