Finite element approximation of spectral problems with Neumann boundary conditions on curved domains

This paper deals with the finite element approximation of the spectral problem for the Laplace equation with Neumann boundary conditions on a curved nonconvex domain . Convergence and optimal order error estimates are proved for standard piecewise linear continuous elements on a discrete polygonal domain in the framework of the abstract spectral approximation theory.

     

Finite element approximation of spectral acoustic problems on curved domains

Summary This paper deals with the finite element approximation of the displacement formulation of the spectral acoustic problem on a curved non convex two-dimensional domain . Convergence and error estimates are proved for Raviart-Thomas elements on a discrete polygonal domain _h in the framework of the abstract spectral approximation theory. Similar results have been previously proved only for polygonal domains. Numerical tests confirming the theoretical results are reported.Mathematics Subject Classification (2000):65N25, 65N30, 70J30Supported by FONDECYT 2000114 (Chile)Supported by FONDAP in Applied Mathematics (Chile)     

Finite element approximation to a contact problem for a nonlinear thermoviscoelastic beam

In this paper we propose and analyze a finite element method to the solution of a quasi-static contact problem between a nonlinear beam and a rigid obstacle. Error estimates and energy decay are obtained and some numerical simulations described.     

Stabilized dual-mixed method for the problem of linear elasticity with mixed boundary conditions

We extend the applicability of the augmented dual-mixed method introduced recently in Gatica (2007), Gatica et al. (2009) to the problem of linear elasticity with mixed boundary conditions. The method is based on the Hellinger–Reissner principle and the symmetry of the stress tensor is imposed in a weak sense. The Neumann boundary condition is prescribed in the finite element space. Then, suitable Galerkin least-squares type terms are added in order to obtain an augmented variational formulation which is coercive in the whole space. This allows to use any finite element subspaces to approximate the displacement, the Cauchy stress tensor and the rotation.     

Error estimates for the finite element approximation of linear elastic equations in an unbounded domain

In this paper we present error estimates for the finite element approximation of linear elastic equations in an unbounded domain. The finite element approximation is formulated on a bounded computational domain using a nonlocal approximate artificial boundary condition or a local one. In fact there are a family of nonlocal approximate boundary conditions with increasing accuracy (and computational cost) and a family of local ones for a given artificial boundary. Our error estimates show how the errors of the finite element approximations depend on the mesh size, the terms used in the approximate artificial boundary condition, and the location of the artificial boundary. A numerical example for Navier equations outside a circle in the plane is presented. Numerical results demonstrate the performance of our error estimates.

     

Finite element approximation of -surfaces

In this paper a piecewise linear finite element approximation of -surfaces, or surfaces with constant mean curvature, spanned by a given Jordan curve in is considered. It is proved that the finite element -surfaces converge to the exact -surfaces under the condition that the Jordan curve is rectifiable. Several numerical examples are given.

     

Boundary element approximation for Maxwell's eigenvalue problem

We introduce a new method for computing eigenvalues of the Maxwell operator with boundary finite elements. On bounded domains with piecewise constant material coefficients, the Maxwell solution for fixed wave number can be represented by boundary integrals, which allows to reduce the eigenvalue problem to a nonlinear problem for determining the wave number along with boundary and interface traces. A Galerkin discretization yields a smooth nonlinear matrix eigenvalue problem that is solved by Newton's method or, alternatively, the contour integral method. Several numerical results including an application to the band structure computation of a photonic crystal illustrate the efficiency of this approach. Copyright © 2013 John Wiley & Sons, Ltd.     

A note on the approximation of Dirichlet boundary control problems for the wave equation on curved domains

We consider an exact boundary control problem for the wave equation with given initial and terminal data and Dirichlet boundary control. The aim is to steer the state of the system that is defined on a given domain to a position of rest in finite time. The optimal control that is obtained as the solution of the problem depends on the data that define the problem, in particular on the domain. Often for the numerical solution of the control problem, this given domain is replaced by a polygon. This is the motivation to study the convergence of the optimal controls for the polygon to the optimal controls for the given domain. To study the convergence, the values of the optimal controls that are defined on the boundaries of the approximating polygons are mapped in the normal directions of the polygon to control functions defined on the boundary of the original domain. This map has already been used by Bramble and King, Deckelnick, Guenther and Hinze and by Casas and Sokolowski. Using this map, we can show the strong convergence of the transformed controls as the polygons approach the given domain. An essential tool to obtain the convergence is a regularization term in the objective functions to increase the regularity of the state.     

Finite element theory on curved domains with applications to discontinuous Galerkin finite element methods

In this paper we provide key estimates used in the stability and error analysis of discontinuous Galerkin finite element methods (DGFEMs) on domains with curved boundaries. In particular, we review trace estimates, inverse estimates, discrete Poincaré–Friedrichs' inequalities, and optimal interpolation estimates in noninteger Hilbert–Sobolev norms, that are well known in the case of polytopal domains. We also prove curvature bounds for curved simplices, which does not seem to be present in the existing literature, even in the polytopal setting, since polytopal domains have piecewise zero curvature. We demonstrate the value of these estimates, by analyzing the IPDG method for the Poisson problem, introduced by Douglas and Dupont, and by analyzing a variant of the hp-DGFEM for the biharmonic problem introduced by Mozolevski and Süli. In both cases we prove stability estimates and optimal a priori error estimates. Numerical results are provided, validating the proven error estimates.     

Finite element formulation of thermal boundary layers

In this article we discuss the finite element discretization of the two-dimensional, incompressible, and turbulent boundary layers. The formulation of the momentum equation is essentially due to Baker and Soliman [1] with some modifications.The versatility and the accuracy of the method is established by considering several test cases. The predictions are satisfactory and compare favorably with alternative numerical techniques.     

Finite element approximation to a contact problem in linear thermoelasticity

A finite element approximation to the solution of a one-dimensional linear thermoelastic problem with unilateral contact of the Signorini type and heat flux is proposed. An error bound is derived and some numerical experiments are performed.

     

Analysis of a finite element method for pressure/potential formulation of elastoacoustic spectral problems

A finite element method to approximate the vibration modes of a structure enclosing an acoustic fluid is analyzed. The fluid is described by using simultaneously pressure and displacement potential variables, whereas displacement variables are used for the solid. A mathematical analysis of the continuous spectral problem is given. The problem is discretized on a simplicial mesh by using piecewise constant elements for the pressure and continuous piecewise linear finite elements for the other fields. Error estimates are settled for approximate eigenvalues and eigenfrequencies. Finally, implementation issues are discussed.

     

Error estimates of the finite element method for the exterior Helmholtz problem with a modified DtN boundary condition

A priori error estimates in the H¹- and L²-norms are established for the finite element method applied to the exterior Helmholtz problem, with modified Dirichlet-to-Neumann (MDtN) boundary condition. The error estimates include the effect of truncation of the MDtN boundary condition as well as that of discretization of the finite element method. The error estimate in the L²-norm is sharper than that obtained by the author [D. Koyama, Error estimates of the DtN finite element method for the exterior Helmholtz problem, J. Comput. Appl. Math. 200 (1) (2007) 21-31] for the truncated DtN boundary condition.     

Finite element approximation of the Neumann eigenvalue problem in domains with multiple cracks

^** Email: belhach{at}poncelet.univ-metz.fr^*** Email: bucur{at}math.univ-metz.fr^**** Email: jmse{at}math.univ-metz.fr We study the Neumann–Laplacian eigenvalue problem in domainswith multiple cracks. We derive a mixed variational formulationwhich holds on the whole geometric domain (including the cracks)and implements efficient finite-element discretizations forthe computation of eigenvalues. Optimal error estimates aregiven and several numerical examples are presented, confirmingthe efficiency of the method. As applications, we numericallyinvestigate the behaviour of the low eigenvalues in domainswith a large number of cracks.     

Ciarlet–Raviart mixed finite element approximation for an optimal control problem governed by the first bi-harmonic equation

The Ciarlet–Raviart mixed finite element approximation is constructed to solve the constrained optimal control problem governed by the first bi-harmonic equation. The optimality conditions consisting of the state and the co-state equations is derived. Also, the a priori error estimates are analyzed. In the analysis of the a priori error estimates, the improved convergent rate of the higher order than existed results is proved. Some numerical experiments are performed to confirm the theoretical analysis for the a priori error estimate.     

Mortar spectral element discretization of the stokes problem in axisymmetric domains

The Stokes problem in a tri‐dimensional axisymmetric domain results into a countable family of two‐dimensional problems when using the Fourier coefficients with respect to the angular variable. Relying on this dimension reduction, we propose and study a mortar spectral element discretization of the problem. Numerical experiments confirm the efficiency of this method. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 44–73, 2014     

The structure of the boundary element matrix for the three-dimensional Dirichlet problem in elasticity

The algebraic properties of the matrix arising for the three-dimensional Dirichlet problem for Lamé equations in a rotational domain by the boundary element method are considered. The use of the special basis leads to a matrix having a block structure with sparse blocks. The possible strategies for the efficient solution of the above problem are discussed. © 1998 John Wiley & Sons, Ltd.     

Legendre-Laguerre coupled spectral element methods for second- and fourth-order equations on the half line

Some Legendre spectral element/Laguerre spectral coupled methods are proposed to numerically solve second- and fourth-order equations on the half line. The proposed methods are based on splitting the infinite domain into two parts, then using the Legendre spectral element method in the finite subdomain and Laguerre method in the infinite subdomain. C⁰ or C¹-continuity, according to the problem under consideration, is imposed to couple the two methods. Rigorous error analysis is carried out to establish the convergence of the method. More importantly, an efficient computational process is introduced to solve the discrete system. Several numerical examples are provided to confirm the theoretical results and the efficiency of the method.