805 research outputs found
A Nitsche Finite Element Approach for Elliptic Problems with Discontinuous Dirichlet Boundary Conditions
We present a numerical approximation method for linear diffusion-reaction
problems with possibly discontinuous Dirichlet boundary conditions. The
solution of such problems can be represented as a linear combination of
explicitly known singular functions as well as of an -regular part. The
latter part is expressed in terms of an elliptic problem with regularized
Dirichlet boundary conditions, and can be approximated by means of a Nitsche
finite element approach. The discrete solution of the original problem is then
defined by adding the singular part of the exact solution to the Nitsche
approximation. In this way, the discrete solution can be shown to converge of
second order with respect to the mesh size
A note on the penalty parameter in Nitsche's method for unfitted boundary value problems
Nitsche's method is a popular approach to implement Dirichlet-type boundary
conditions in situations where a strong imposition is either inconvenient or
simply not feasible. The method is widely applied in the context of unfitted
finite element methods. From the classical (symmetric) Nitsche's method it is
well-known that the stabilization parameter in the method has to be chosen
sufficiently large to obtain unique solvability of discrete systems. In this
short note we discuss an often used strategy to set the stabilization parameter
and describe a possible problem that can arise from this. We show that in
specific situations error bounds can deteriorate and give examples of
computations where Nitsche's method yields large and even diverging
discretization errors
The diffuse Nitsche method: Dirichlet constraints on phase-field boundaries
We explore diffuse formulations of Nitsche's method for consistently imposing Dirichlet boundary conditions on phase-field approximations of sharp domains. Leveraging the properties of the phase-field gradient, we derive the variational formulation of the diffuse Nitsche method by transferring all integrals associated with the Dirichlet boundary from a geometrically sharp surface format in the standard Nitsche method to a geometrically diffuse volumetric format. We also derive conditions for the stability of the discrete system and formulate a diffuse local eigenvalue problem, from which the stabilization parameter can be estimated automatically in each element. We advertise metastable phase-field solutions of the Allen-Cahn problem for transferring complex imaging data into diffuse geometric models. In particular, we discuss the use of mixed meshes, that is, an adaptively refined mesh for the phase-field in the diffuse boundary region and a uniform mesh for the representation of the physics-based solution fields. We illustrate accuracy and convergence properties of the diffuse Nitsche method and demonstrate its advantages over diffuse penalty-type methods. In the context of imaging based analysis, we show that the diffuse Nitsche method achieves the same accuracy as the standard Nitsche method with sharp surfaces, if the inherent length scales, i.e., the interface width of the phase-field, the voxel spacing and the mesh size, are properly related. We demonstrate the flexibility of the new method by analyzing stresses in a human vertebral body
A Nitsche-based domain decomposition method for hypersingular integral equations
We introduce and analyze a Nitsche-based domain decomposition method for the
solution of hypersingular integral equations. This method allows for
discretizations with non-matching grids without the necessity of a Lagrangian
multiplier, as opposed to the traditional mortar method. We prove its almost
quasi-optimal convergence and underline the theory by a numerical experiment.Comment: 21 pages, 5 figure
- …