An official website of the United States government
Abstract In this paper we study the biharmonic equation with Navier boundary conditions in a polygonal domain. In particular, we propose a method that effectively decouples the fourth-order problem as a system of Poisson equations. Our method differs from the naive mixed method that leads to two Poisson problems but only applies to convex domains; our decomposition involves a third Poisson equation to confine the solution in the correct function space, and therefore can be used in both convex and nonconvex domains. A $C^0$ finite element algorithm is in turn proposed to solve the resulting system. In addition, we derive optimal error estimates for the numerical solution on both quasi-uniform meshes and graded meshes. Numerical test results are presented to justify the theoretical findings.
more »
« less
Geometric transformation of finite element methods: Theory and applications
We present a new technique to apply finite element methods to partial differential equations over curved domains. A change of variables along a coordinate transformation satisfying only low regularity assumptions can translate a Poisson problem over a curved physical domain to a Poisson problem over a polyhedral parametric domain. This greatly simplifies both the geometric setting and the practical implementation, at the cost of having globally rough non-trivial coefficients and data in the parametric Poisson problem. Our main result is that a recently developed broken Bramble-Hilbert lemma is key in harnessing regularity in the physical problem to prove higher-order finite element convergence rates for the parametric problem. Numerical experiments are given which confirm the predictions of our theory.
more »
« less
- Award ID(s):
- 2012857
- PAR ID:
- 10524676
- Publisher / Repository:
- Elsevier
- Date Published:
- Journal Name:
- Applied Numerical Mathematics
- Volume:
- 192
- Issue:
- C
- ISSN:
- 0168-9274
- Page Range / eLocation ID:
- 389 to 413
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We consider an unconstrained tangential Dirichlet boundary control problem for the Stokes equations with an $ L^2 $ penalty on the boundary control. The contribution of this paper is twofold. First, we obtain well-posedness and regularity results for the tangential Dirichlet control problem on a convex polygonal domain. The analysis contains new features not found in similar Dirichlet control problems for the Poisson equation; an interesting result is that the optimal control has higher local regularity on the individual edges of the domain compared to the global regularity on the entire boundary. Second, we propose and analyze a hybridizable discontinuous Galerkin (HDG) method to approximate the solution. For convex polygonal domains, our theoretical convergence rate for the control is optimal with respect to the global regularity on the entire boundary. We present numerical experiments to demonstrate the performance of the HDG method.more » « less
-
We obtain regularity results in weighted Sobolev spaces for the solution of the obstacle problem for the integral fractional Laplacian [Formula: see text] in a Lipschitz bounded domain [Formula: see text] satisfying the exterior ball condition. The weight is a power of the distance to the boundary [Formula: see text] of [Formula: see text] that accounts for the singular boundary behavior of the solution for any [Formula: see text]. These bounds then serve us as a guide in the design and analysis of a finite element scheme over graded meshes for any dimension [Formula: see text], which is optimal for [Formula: see text].more » « less
-
null (Ed.)Abstract The paper studies a higher order unfitted finite element method for the Stokes system posed on a surface in ℝ 3 . The method employs parametric P k - P k −1 finite element pairs on tetrahedral bulk mesh to discretize the Stokes system on embedded surface. Stability and optimal order convergence results are proved. The proofs include a complete quantification of geometric errors stemming from approximate parametric representation of the surface. Numerical experiments include formal convergence studies and an example of the Kelvin--Helmholtz instability problem on the unit sphere.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1016/j.apnum.2023.07.002
