Implicitly enriched Galerkin methods for numerical solutions of fourth‐order partial differential equations containing singularities. Issue 6 (2nd May 2018)

• Record Type: Journal Article
• Authors: Kim, Sinae; Oh, Hae‐Soo; Palta, Birce; Kim, Hyunju
• Is Part Of: Numerical methods for partial differential equations. Volume 34:Issue 6(2018)
• Page Start: 2079
• Page End: 2112
• Publication Date: 2018-05-02
• Subjects: B‐spline (NURBS) -- biharmonic equation -- enriched isogeometric analysis(finite element method) -- Galerkin method with B‐spline basis functions -- Partition of Unity function with flat top; Differential equations, Partial -- Numerical solutions -- Periodicals; 515.353
• DOI: 10.1002/num.22277
• Languages: English
• ISSNs: 0749-159X

Abstract

Abstract : Highlights are the following: For any integer n ≥ 0, we construct C n ‐continuous partition of unity (PU) functions with flat‐top from B‐spline functions to have numerical solutions of fourth‐order equations with singularities. B‐spline functions are modified to satisfy clamped boundary conditions. To handle singularity arising in fourth‐order elliptic differential equations, these modified B‐spline functions are enriched either by introducing enrichment basis functions implicitly through particular geometric mappings or by adding singular basis functions explicitly. To show the effectiveness of the proposed implicit enrichment methods (mapping method), the accuracy, the number of degrees of freedom (DOF), and matrix condition numbers are computed and compared in the h ‐refinement, the p ‐refinement, and the k ‐refinement of the approximation space of B‐spline basis functions. Using Partition of unity (PU) functions with flat‐top, B‐spline functions are modified to satisfy boundary conditions of the fourth‐order equations. Since the standard isogeometric analysis (IGA) as well as the conventional FEM have limitations in handling fourth‐order differential equations containing singularities, we consider two enrichment methods (explicit and implicit) in the framework of the p ‐, the k, and the h ‐refinements of IGA. We demonstrate that both enrichment methods yield good approximate solutions, but explicit enrichment methods give large (almost singular) matrix condition numbers and face integrating singular functions. Because of these limitations of external enrichment methods, we extensively investigate implicit enrichment methods (mapping methods) that virtually convert fourth‐order elliptic problems with singularities to problems with no influence of the singularities. Effectiveness of the proposed mapping method extensively tested to one‐dimensional fourth‐order equation with singularities. The implicit enrichment (mapping) method is extended to the two‐dimensional cases and test it to fourth‐order partial differential equations on cracked domains.