Provable stable arbitrary order symmetric interior penalty discontinuous Galerkin (SIP) discretisations of variable viscosity, incompressible Stokes flow utilising @math -- @math elements and hierarchical Legendre basis polynomials are developed and investigated.For solving the resulting linear system, a block preconditioned iterative method is proposed. The nested viscous problem is solved by a @math -multilevel preconditioned Krylov subspace method. For the @math -coarsening, a twolevel method utilising element-block Jacobi preconditioned iterations as a smoother is employed. Piecewise bilinear ( @math ) and piecewise constant ( @math ) @math -coarse spaces are considered. Finally, Galerkin @math -coarsening is proposed and investigated for the two @math -coarse spaces considered. Through a number of numerical experiments, we demonstrate that utilising the @math coarse space results in the most robust @math -multigrid method for variable viscosity Stokes flow. Using this @math coarse space we observe that the convergence of the overall Stokes solver appears to be robust with respect to the jump in the viscosity and only mildly depending on the polynomial order @math . It is demonstrated and supported by theoretical results that the convergence of the SIP discretisations and the iterative methods rely on a sharp choice of the penalty parameter based on local values of the viscosity.
Inf-sup stable discontinuous Galerkin (DG) methods of the interior penalty type for the Stokes equations can be constructed by using the tensor product element pairs @math -- @math , and @math -- @math @cite_23 @cite_39 , as well as the @math -conforming Raviart-Thomas, Brezzi-Douglas-Marini, and Brezzi-Douglas-Fortin-Marini kind element pairs; see @cite_12 @cite_3 @cite_48 @cite_27 .
@cite_48 @cite_3 @cite_39 @cite_27 @cite_23 @cite_12
140909251 1982359839 2225768769 2100903519
Discontinuous Galerkin (DG) methods for solving partial differential equations, developed in the late 1990s, have become popular among computational scientists. This book covers both theory and computation as it focuses on three primal DG methods--the symmetric interior penalty Galerkin, incomplete interior penalty Galerkin, and nonsymmetric interior penalty Galerkin which are variations of interior penalty methods. The author provides the basic tools for analysis and discusses coding issues, including data structure, construction of local matrices, and assembling of the global matrix. Computational examples and applications to important engineering problems are also included. Discontinuous Galerkin Methods for Solving Elliptic and Parabolic Equations: Theory and Implementation is divided into three parts: Part I focuses on the application of DG methods to second order elliptic problems in one dimension and in higher dimensions. Part II presents the time-dependent parabolic problems without and with convection. Part III contains applications of DG methods to solid mechanics (linear elasticity), fluid dynamics (Stokes and Navier Stokes), and porous media flow (two-phase and miscible displacement). Appendices contain proofs and MATLAB code for one-dimensional problems for elliptic equations and routines written in C that correspond to algorithms for the implementation of DG methods in two or three dimensions. Audience: This book is intended for numerical analysts, computational and applied mathematicians interested in numerical methods for partial differential equations or who study the applications discussed in the book, and engineers who work in fluid dynamics and solid mechanics and want to use DG methods for their numerical results. The book is appropriate for graduate courses in finite element methods, numerical methods for partial differential equations, numerical analysis, and scientific computing. Chapter 1 is suitable for a senior undergraduate class in scientific computing. Contents: List of Figures; List of Tables; List of Algorithms; Preface; Part I: Elliptic Problems; Chapter 1: One-dimensional problem; Chapter 2: Higher dimensional problem; Part II: Parabolic Problems; Chaper 3: Purely parabolic problems; Chapter 4: Parabolic problems with convection; Part III: Applications; Chapter 5: Linear elasticity; Chapter 6: Stokes flow; Chapter 7: Navier-Stokes flow; Chapter 8: Flow in porous media; Appendix A: Quadrature rules; Appendix B: DG codes; Appendix C: An approximation result; Bibliography; Index. We propose and analyze a discontinuous Galerkin approximation for the Stokes problem. The finite element triangulation employed is not required to be conforming and we use discontinuous pressures and velocities. No additional unknown fields need to be introduced, but only suitable bilinear forms defined on the interfaces between the elements, involving the jumps of the velocity and the average of the pressure. We consider hp approximations using ℚk′–ℚk velocity-pressure pairs with k′ = k + 2, k + 1, k. Our methods show better stability properties than the corresponding conforming ones. We prove that our first two choices of velocity spaces ensure uniform divergence stability with respect to the mesh size h. Numerical results show that they are uniformly stable with respect to the local polynomial degree k, a property that has no analog in the conforming case. An explicit bound in k which is not sharp is also proven. Numerical results show that if equal order approximation is chosen for the velocity and pressure, no spurious pressure modes are present but the method is not uniformly stable either with respect to h or k. We derive a priori error estimates generalizing the abstract theory of mixed methods. Optimal error estimates in h are proven. As for discontinuous Galerkin methods for scalar diffusive problems, half of the power of k is lost for p and hp pproximations independently of the divergence stability. A multigrid method for the Stokes system discretized with an Hdiv-conforming discontinuous Galerkin method is presented. It acts on the combined velocity and pressure spaces and thus does not need a Schur complement approximation. The smoothers used are of overlapping Schwarz type and employ a local Helmholtz decomposition. Additionally, we use the fact that the discretization provides nested divergence free subspaces. We present the convergence analysis and numerical evidence that convergence rates are not only independent of mesh size, but also reasonably small. We present a compact discontinuous Galerkin (CDG) method for an elliptic model problem. The problem is first cast as a system of first order equations by introducing the gradient of the primal unknown, or flux, as an additional variable. A standard discontinuous Galerkin (DG) method is then applied to the resulting system of equations. The numerical interelement fluxes are such that the equations for the additional variable can be eliminated at the element level, thus resulting in a global system that involves only the original unknown variable. The proposed method is closely related to the local discontinuous Galerkin (LDG) method [B. Cockburn and C.-W. Shu, SIAM J. Numer. Anal., 35 (1998), pp. 2440-2463], but, unlike the LDG method, the sparsity pattern of the CDG method involves only nearest neighbors. Also, unlike the LDG method, the CDG method works without stabilization for an arbitrary orientation of the element interfaces. The computation of the numerical interface fluxes for the CDG method is slightly more involved than for the LDG method, but this additional complication is clearly offset by increased compactness and flexibility. Compared to the BR2 [F. Bassi and S. Rebay, J. Comput. Phys., 131 (1997), pp. 267-279] and IP [J. Douglas, Jr., and T. Dupont, in Computing Methods in Applied Sciences (Second Internat. Sympos., Versailles, 1975), Lecture Notes in Phys. 58, Springer, Berlin, 1976, pp. 207-216] methods, which are known to be compact, the present method produces fewer nonzero elements in the matrix and is computationally more efficient.