$$This paper presents the numerical analysis for a variational formulation of rate-independent phase transformations in elastic solids due to Mielke et al. The new model itself suggests an implicit time-discretization which is combined with the finite element method in space. A priori error estimates are established for the quasioptimal spatial approximation of the stress field within one time-step. A posteriori error estimates motivate an adaptive mesh-refining algorithm for efficient discretization....

This paper presents the numerical analysis for a variational formulation of rate-independent phase transformations in elastic solids due to Mielke et al. The new model itself suggests an implicit time-discretization which is combined with the finite element method in space. A priori error estimates are established for the quasioptimal spatial approximation of the stress field within one time-step. A posteriori error estimates motivate an adaptive mesh-refining algorithm for efficient...

Numerical analysis of a model Stokes interface problem with the homogeneous Dirichlet boundary condition is considered. The interface condition is interpreted as an additional singular force field to the Stokes equations using the characteristic function. The finite element method is applied after introducing a regularization of the singular source term. Consequently, the error is divided into the regularization and discretization parts which are studied separately. As a result, error estimates...

This paper is concerned with the dual formulation of the interface problem consisting of a linear partial differential equation with variable coefficients in some bounded Lipschitz domain Ω in ${ℝ}^n$ (n ≥ 2) and the Laplace equation with some radiation condition in the unbounded exterior domain Ωc:= ${ℝ}^n\setminus \overline{\Omega }$. The two problems are coupled by transmission and Signorini contact conditions on the interface Γ = ∂Ω. The exterior part of the interface problem is rewritten using a Neumann to Dirichlet mapping (NtD)...

This paper is concerned with the dual formulation of the interface problem consisting of a linear partial differential equation with variable coefficients in some bounded Lipschitz domain Ω in ${ℝ}^n$ (n ≥ 2) and the Laplace equation with some radiation condition in the unbounded exterior domain Ωc := ${ℝ}^n\setminus \overline{\Omega }$. The two problems are coupled by transmission and Signorini contact conditions on the interface Γ = ∂Ω. The exterior part of the interface problem is rewritten using a Neumann to Dirichlet mapping...

This paper deals with the numerical approximation of mild solutions of elliptic-parabolic equations, relying on the existence results of Bénilan and Wittbold (1996). We introduce a new and simple algorithm based on Halpern’s iteration for nonexpansive operators (Bauschke, 1996; Halpern, 1967; Lions, 1977), which is shown to be convergent in the degenerate case, and compare it with existing schemes (Jäger and Kačur, 1995; Kačur, 1999).

This paper deals with the numerical approximation of mild solutions of elliptic-parabolic equations, relying on the existence results of Bénilan and Wittbold (1996). We introduce a new and simple algorithm based on Halpern's iteration for nonexpansive operators (Bauschke, 1996; Halpern, 1967; Lions, 1977), which is shown to be convergent in the degenerate case, and compare it with existing schemes (Jäger and Kačur, 1995; Kačur, 1999).

In this paper, we carry out the numerical analysis of a distributed optimal control problem governed by a quasilinear elliptic equation of non-monotone type. The goal is to prove the strong convergence of the discretization of the problem by finite elements. The main issue is to get error estimates for the discretization of the state equation. One of the difficulties in this analysis is that, in spite of the partial differential equation has a unique solution for any given control, the uniqueness...

In this paper, we carry out the numerical analysis of a distributed optimal control problem governed by a quasilinear elliptic equation of non-monotone type. The goal is to prove the strong convergence of the discretization of the problem by finite elements. The main issue is to get error estimates for the discretization of the state equation. One of the difficulties in this analysis is that, in spite of the partial differential equation has a unique solution for any given control, the uniqueness...

Many problems in quantum chemistry deal with the computation of fundamental or excited states of molecules and lead to the resolution of eigenvalue problems. One of the major difficulties in these computations lies in the very large dimension of the systems to be solved. Indeed these eigenfunctions depend on $3n$ variables where $n$ stands for the number of particles (electrons and/or nucleari) in the molecule. In order to diminish the size of the systems to be solved, the chemists have proposed many...

Many problems in quantum chemistry deal with the computation of fundamental or excited states of molecules and lead to the resolution of eigenvalue problems. One of the major difficulties in these computations lies in the very large dimension of the systems to be solved. Indeed these eigenfunctions depend on 3n variables where n stands for the number of particles (electrons and/or nucleari) in the molecule. In order to diminish the size of the systems to be solved, the chemists have proposed...

The present paper deals with solving the general biharmonic problem by the finite element method using curved triangular finit $C^1$-elements introduced by Ženíšek. The effect of numerical integration is analysed in the case of mixed boundary conditions and sufficient conditions for the uniform $V_{Oh}$-ellipticity are found.

The meshless element-free Galerkin method is developed for numerical analysis of hyperbolic initial-boundary value problems. In this method, only scattered nodes are required in the domain. Computational formulae of the method are analyzed in detail. Error estimates and convergence are also derived theoretically and verified numerically. Numerical examples validate the performance and efficiency of the method.

The Method of Fundamental Solutions (MFS) is a boundary-type meshless method for the solution of certain elliptic boundary value problems. In this work, we investigate the properties of the matrices that arise when the MFS is applied to the Dirichlet problem for Laplace’s equation in a disk. In particular, we study the behaviour of the eigenvalues of these matrices and the cases in which they vanish. Based on this, we propose a modified efficient numerical algorithm for the solution of the problem...

The Method of Fundamental Solutions (MFS) is a boundary-type meshless method for the solution of certain elliptic boundary value problems. In this work, we investigate the properties of the matrices that arise when the MFS is applied to the Dirichlet problem for Laplace's equation in a disk. In particular, we study the behaviour of the eigenvalues of these matrices and the cases in which they vanish. Based on this, we propose a modified efficient numerical algorithm for the solution of the problem...

In this article, we provide a priorierror estimates for the spectral and pseudospectral Fourier (also called planewave) discretizations of the periodic Thomas-Fermi-von Weizsäcker (TFW) model and for the spectral discretization of the periodic Kohn-Sham model, within the local density approximation (LDA). These models allow to compute approximations of the electronic ground state energy and density of molecular systems in the condensed phase. The TFW model is strictly convex with respect to the...

In this article, we provide a priori error estimates for the spectral and pseudospectral Fourier (also called planewave) discretizations of the periodic Thomas-Fermi-von Weizsäcker (TFW) model and for the spectral discretization of the periodic Kohn-Sham model, within the local density approximation (LDA). These models allow to compute approximations of the electronic ground state energy and density of molecular systems in the condensed phase. The TFW model is strictly convex with respect to the...