Discontinuous Galerkin Methods for Isogeometric Analysis for Elliptic Equations on Surfaces

Abstract

We propose a method that combines isogeometric analysis (IGA) with the discontinuous Galerkin (DG) method for solving elliptic equations on 3-dimensional (3D) surfaces consisting of multiple patches. DG ideology is adopted across the patch interfaces to glue the multiple patches, while the traditional IGA, which is very suitable for solving partial differential equations (PDEs) on (3D) surfaces, is employed within each patch. Our method takes advantage of both IGA and the DG method. Firstly, the time-consuming steps in mesh generation process in traditional finite element analysis (FEA) are no longer necessary and refinements, including $h$-refinement and $p$-refinement which both maintain the original geometry, can be easily performed by knot insertion and order-elevation (Farin, in Curves and surfaces for CAGD, 2002). Secondly, our method can easily handle the cases with non-conforming patches and different degrees across the patches. Moreover, due to the geometric flexibility of IGA basis functions, especially the use of multiple patches, we can get more accurate modeling of more complex surfaces. Thus, the geometrical error is significantly reduced and it is, in particular, eliminated for all conic sections. Finally, this method can be easily formulated and implemented. We generally describe the problem and then present our primal formulation. A new ideology, which directly makes use of the approximation property of the NURBS basis functions on the parametric domain rather than that of the IGA functions on the physical domain (the former is easier to get), is adopted when we perform the theoretical analysis including the boundedness and stability of the primal form, and the error analysis under both the DG norm and the ${\mathcal {L}}^{2}$ norm. The result of the error analysis shows that our scheme achieves the optimal convergence rate with respect to both the DG norm and the ${\mathcal {L}}^{2}$ norm. Numerical examples are presented to verify the theoretical result and gauge the good performance of our method.

Keywords

Elliptic equation on (3D) surface / Isogeometric analysis / Discontinuous Galerkin method / Multiple patches / Non-conforming patches / The optimal convergence rate

References

[1] Farin G. Curves and Surfaces for CAGD. 2002 5 San Francisco: Morgan Kaufmann Publishers

[2] Hughes TJR, Cottrell JA, Bazilevs Y. Isogeometric analysis: CAD, finite elements, NURBS, exact geometry, and mesh refinement. Comput. Methods Appl. Mech. Eng.. 2005, 194 4135-4195

[3] Bazilevs Y, de Veiga B, Cottrell JA, Hughes TJR, Sangalli G. Isogeometric analysis: approximation, stability and error estimates for h-refined meshes. Math. Models Methods Appl. Sci.. 2006, 16 07 1031-1090

[4] Bazilevs Y, Calo VM, Cottrell JA, Hughes TJR, Reali A, Scovazzi G. Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows. Comp. Methods Appl. Mech. Eng.. 2007, 197 173-201

[5] Bazilevs Y, Calo VM, Hughes TJR, Zhang Y. Isogeometric fluid-structure interaction: theory, algorithms, and computations. Comput. Mech.. 2008, 43 3-37

[6] Cottrell JA, Reali A, Bazilevs Y, Hughes TJR. Isogeometric analysis of structural vibrations. Comput. Methods Appl. Mech. Eng.. 2006, 195 5257-5296

[7] Cottrell JA, Hughes TJR, Reali A. Studies of refinement and continuity in isogeometric analysis. Comput. Methods Appl. Mech. Eng.. 2007, 196 4160-4183

[8] Cottrell JA, Hughes TJR, Bazilevs Y. Isogeometric Analysis: Toward Integration of CAD and FEA. 2009 Chichester: Wiley

[9] Reed, W.H., Hill, T.R.: Triangular mesh methods for the neutron transport equation, Los Alamos Report LA-UR-73-479 (1973)

[10] Cockburn B, Shu C-W. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws. II. General framework. Math. Comput.. 1989, 52 186 411-435

[11] Cockburn B, Lin SY, Shu C-W. TVB Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws III: One-dimensional systems. J. Comput. Phys.. 1989, 84 1 90-113

[12] Cockburn B, Hou S, Shu C-W. The Runge-Kutta local projection discontinuous Galerkin finite element method for conservation laws IV: The multidimensional case. Math. Comput.. 1990, 54 190 545-581

[13] Cockburn B, Shu C-W. The Runge-Kutta local projection $P^{1}$-discontinuous-Galerkin finite element method for scalar conservation laws. RAIRO Model. Math. Anal. Numer.. 1991, 25 3 337-361

[14] Cockburn B, Shu C-W. The Runge-Kutta discontinuous Galerkin method for conservation laws V: multidimensional systems. J. Comput. Phys.. 1998, 141 2 199-224

[15] Brezzi, F., Manzini, G., Marini, D., Pietra, P., Russo, A.:Discontinuous finite elements for diffusion problems, Atti Convegno in onore di F. Brioschi (Milano 1997), Istituto Lombardo. Accademia di Scienze e Lettere, pp. 197–217 (1999)

[16] Brezzi F, Manzini G, Marini D, Pietra P, Russo A. Discontinuous Galerkin approximations for elliptic problems. Numer. Methods Partial Differ. Equ.. 2000, 16 365-378

[17] Castillo P, Cockburn B, Perugia I . An a priori error analysis of the local discontinuous Galerkin method for elliptic problems. SIAM J. Numer. Anal.. 2000, 38 1676-1706

[18] Cockburn B, Dawson C. Whiteman JR. Some extensions of the local discontinuous Galerkin method for convection-diffusion equations in multidimensions. Proceedings of the Conference on the Mathematics of Finite Elements and Applications: MAFELAP X. 2000 NewYork: Elsevier. 225-238

[19] Cockburn B, Kanschat G, Perugia I . Superconvergence of the local discontinuous Galerkin method for elliptic problems on Cartesian grids. SIAM J. Numer. Anal.. 2001, 39 1 264-285

[20] Arnold, D.N.: An interior penalty finite element method with discontinuous elements. Ph.D. thesis, The University of Chicago, Chicago, IL (1979)

[21] Arnold DN. An interior penalty finite element method with discontinuous elements. SIAM J. Numer. Anal.. 1982, 19 742-760

[22] Baker GA. Finite element methods for elliptic equations using nonconforming elements. Math. Comput.. 1977, 31 45-59

[23] Douglas J Jr, Dupont T. Interior Penalty Procedures for Elliptic and Parabolic Galerkin Methods, Computing methods in applied sciences. 1976 Berlin: Springer

[24] Wheeler MF. An elliptic collocation-finite element method with interior penalties. SIAM J. Numer. Anal.. 1978, 15 152-161

[25] Arnold DN, Brezzi F, Cockburn B, Marini LD. Unified analysis of Discontinuous Galerkin methods for elliptic problems. SIAM J. Numer. Anal.. 2002, 39 1749-1779

[26] Auricchio F, da Veiga LB, Buffa A . A fully locking-free isogeometric approach for plane linear elasticity problems: a stream function formulation. Comp. Methods Appl. Mech. Eng.. 2007, 197 1 160-172

[27] Bazilevs Y, Hughes TJR. NURBS-based isogeometric analysis for the computation of flows about rotating components. Comput. Mech.. 2008, 43 1 143-150

[28] Brunero, F.: Discontinuous Galerkin Methods for isogeometric analysis. master’s thesis, Università degli Studi di Milano (2012)

[29] Evans JA, Hughes TJR. Isogeometric divergence-conforming B-splines for the steady Navier-Stokes equations. Math. Models Methods Appl. Sci.. 2013, 23 08 1421-1478

[30] Nguyen VP, Kerfriden P, Brino M . Nitsche’s method for two and three dimensional NURBS patch coupling. Comput. Mech.. 2014, 53 1163-1182

[31] Imran Özcan, Ali: A Nitsche-based mortar method for non-matching multi-patch discretizations in the Finite Cell Method. Master’s thesis, Technische Universität München, Germany (2012)

[32] Ruess M, Schillinger D . Weak coupling for isogeometric analysis of non-matching and trimmed multi-patch geometries. Comput. Methods Appl. Mech. Eng.. 2014, 269 46-71

[33] Dziuk G, Elliot CM. Finite element methods for surface PDEs. Acta Numerica. 2013, 22 289-396

[34] Langer, U., Moore, S.E.: Discontinuous Galerkin Isogeometric Analysis of Elliptic PDEs on Surfaces, preprint http://arxiv.org/pdf/1402.1185, 2014. arXiv preprint arXiv:1402.1185 (2014)

[35] Xu GL, Zhang Q. Geometric Partial Differential Equation Methods in Computational Geometry. 2013 Beijing: Science Press

[36] Bonito A, Cascón JM, Morin P, Nochetto RH, Gianazza U. Brezzi F, Franzone PC, Gilardi G. AFEM for geometric PDE: the Laplace–Beltrami operator. Analysis and Numerics of Partial Differential Equations. 2013 Milan: Springer. 257-306

[37] Gould PL. Introduction to Linear Elasticity. 1999 Berlin: Springer