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).
This paper presents a mathematical model for photo-excited carrier decay in a semiconductor. Due to the carrier trapping states and recombination centers in the bandgap, the carrier decay process is defined by the system of nonlinear differential equations. The system of nonlinear ordinary differential equations is approximated by linearized backward Euler scheme. Some a priori estimates of the discrete solution are obtained and the convergence of the linearized backward Euler method is proved....
Parallel replica dynamics is a method for accelerating the computation of processes characterized by a sequence of infrequent events. In this work, the processes are governed by the overdamped Langevin equation. Such processes spend much of their time about the minima of the underlying potential, occasionally transitioning into different basins of attraction. The essential idea of parallel replica dynamics is that the exit distribution from a given well for a single process can be approximated by...
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 variables where 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 -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 -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...
This paper discusses some conceptional questions of the numerical simulation of viscous incompressible flow which are related to the presence of boundaries.
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...