Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userAnalysis of the finite element method for the Laplace--Beltrami equation on surfaces with regions of high curvature using graded meshes
86
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We derive error estimates for the piecewise linear finite element approximation of the Laplace--Beltrami operator on a bounded, orientable, $C^3$, surface without boundary on general shape regular meshes. As an application, we consider a problem where the domain is split into two regions: one which has relatively high curvature and one that has low curvature. Using a graded mesh we prove error estimates that do not depend on the curvature on the high curvature region. Numerical experiments are provided.
rate research
Read More
We develop a finite element method for the Laplace-Beltrami operator on a surface with boundary and nonhomogeneous Dirichlet boundary conditions. The method is based on a triangulation of the surface and the boundary conditions are enforced weakly using Nitsches method. We prove optimal order a priori error estimates for piecewise continuous polynomials of order $k geq 1$ in the energy and $L^2$ norms that take the approximation of the surface and the boundary into account.
An isogeometric approach for solving the Laplace-Beltrami equation on a two-dimensional manifold embedded in three-dimensional space using a Galerkin method based on Catmull-Clark subdivision surfaces is presented and assessed. The scalar-valued Laplace-Beltrami equation requires only C0 continuity and is adopted to elucidate key features and properties of the isogeometric method using Catmull-Clark subdivision surfaces. Catmull-Clark subdivision bases are used to discretise both the geometry and the physical field. A fitting method generates control meshes to approximate any given geometry with Catmull-Clark subdivision surfaces. The performance of the Catmull-Clark subdivision method is compared to the conventional finite element method. Subdivision surfaces without extraordinary vertices show the optimal convergence rate. However, extraordinary vertices introduce error, which decreases the convergence rate. A comparative study shows the effect of the number and valences of the extraordinary vertices on accuracy and convergence. An adaptive quadrature scheme is shown to reduce the error.
Elliptic partial differential equations on surfaces play an essential role in geometry, relativity theory, phase transitions, materials science, image processing, and other applications. They are typically governed by the Laplace-Beltrami operator. We present and analyze approximations by Surface Finite Element Methods (SFEM) of the Laplace-Beltrami eigenvalue problem. As for SFEM for source problems, spectral approximation is challenged by two sources of errors: the geometric consistency error due to the approximation of the surface and the Galerkin error corresponding to finite element resolution of eigenfunctions. We show that these two error sources interact for eigenfunction approximations as for the source problem. The situation is different for eigenvalues, where a novel situation occurs for the geometric consistency error: The degree of the geometric error depends on the choice of interpolation points used to construct the approximate surface. Thus the geometric consistency term can sometimes be made to converge faster than in the eigenfunction case through a judicious choice of interpolation points.
In this work, we present an adaptive unfitted finite element scheme that combines the aggregated finite element method with parallel adaptive mesh refinement. We introduce a novel scalable distributed-memory implementation of the resulting scheme on locally-adapted Cartesian forest-of-trees meshes. We propose a two-step algorithm to construct the finite element space at hand by means of a discrete extension operator that carefully mixes aggregation constraints of problematic degrees of freedom, which get rid of the small cut cell problem, and standard hanging degree of freedom constraints, which ensure trace continuity on non-conforming meshes. Following this approach, we derive a finite element space that can be expressed as the original one plus well-defined linear constraints. Moreover, it requires minimum parallelization effort, using standard functionality available in existing large-scale finite element codes. Numerical experiments demonstrate its optimal mesh adaptation capability, robustness to cut location and parallel efficiency, on classical Poisson $hp$-adaptivity benchmarks. Our work opens the path to functional and geometrical error-driven dynamic mesh adaptation with the aggregated finite element method in large-scale realistic scenarios. Likewise, it can offer guidance for bridging other scalable unfitted methods and parallel adaptive mesh refinement.
In this paper we develop an evolution of the $C^1$ virtual elements of minimal degree for the approximation of the Cahn-Hilliard equation. The proposed method has the advantage of being conforming in $H^2$ and making use of a very simple set of degrees of freedom, namely 3 degrees of freedom per vertex of the mesh. Moreover, although the present method is new also on triangles, it can make use of general polygonal meshes. As a theoretical and practical support, we prove the convergence of the semi-discrete scheme and investigate the performance of the fully discrete scheme through a set of numerical tests.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
