The hierarchical basis multigrid method for convection-diffusion equations.
The paper deals with the application of a fast algorithm for the solution of finite-difference systems for boundary-value problems on a standard domain (e.g. on a rectangle) to the solution of a boundary-value problem on a domain of general shape contained in the standard domain. A simple iterative procedure is suggested for the determination of fictitious right-hand sides for the system on the standard domain so that its solution is the desired one. Under the assumptions that are usual for matrices...
Output least squares stability for the diffusion coefficient in an elliptic equation in dimension two is analyzed. This guarantees Lipschitz stability of the solution of the least squares formulation with respect to perturbations in the data independently of their attainability. The analysis shows the influence of the flow direction on the parameter to be estimated. A scale analysis for multi-scale resolution of the unknown parameter is provided.
Output least squares stability for the diffusion coefficient in an elliptic equation in dimension two is analyzed. This guarantees Lipschitz stability of the solution of the least squares formulation with respect to perturbations in the data independently of their attainability. The analysis shows the influence of the flow direction on the parameter to be estimated. A scale analysis for multi-scale resolution of the unknown parameter is provided.
For a general class of atomistic-to-continuum coupling methods, coupling multi-body interatomic potentials with a P1-finite element discretisation of Cauchy–Born nonlinear elasticity, this paper adresses the question whether patch test consistency (or, absence of ghost forces) implies a first-order error estimate. In two dimensions it is shown that this is indeed true under the following additional technical assumptions: (i) an energy consistency condition, (ii) locality of the interface correction,...
For a general class of atomistic-to-continuum coupling methods, coupling multi-body interatomic potentials with a P1-finite element discretisation of Cauchy–Born nonlinear elasticity, this paper adresses the question whether patch test consistency (or, absence of ghost forces) implies a first-order error estimate. In two dimensions it is shown that this is indeed true under the following additional technical assumptions: (i) an energy consistency condition, (ii) locality of the interface correction,...
The projection-algebraic approach of the Calogero type for discrete approximations of linear and nonlinear differential operator equations in Banach spaces is studied. The solution convergence and realizability properties of the related approximating schemes are analyzed. For the limiting-dense approximating scheme of linear differential operator equations a new convergence theorem is stated. In the case of nonlinear differential operator equations the effective convergence conditions for the approximated...
We estimate the constant in the strengthened Cauchy-Bunyakowski-Schwarz inequality for hierarchical bilinear finite element spaces and elliptic partial differential equations with coefficients corresponding to anisotropy (orthotropy). It is shown that there is a nontrivial universal estimate, which does not depend on anisotropy. Moreover, this estimate is sharp and the same as for hierarchical linear finite element spaces.
A variant of the Total Overlapping Schwarz (TOS) method has been introduced in [Ben Belgacem et al., C. R. Acad. Sci., Sér. 1 Math. 336 (2003) 277–282] as an iterative algorithm to approximate the absorbing boundary condition, in unbounded domains. That same method turns to be an efficient tool to make numerical zooms in regions of a particular interest. The TOS method enjoys, then, the ability to compute small structures one wants to capture and the reliability to obtain the behavior of the solution...
A variant of the Total Overlapping Schwarz (TOS) method has been introduced in [Ben Belgacem et al., C. R. Acad. Sci., Sér. 1 Math.336 (2003) 277–282] as an iterative algorithm to approximate the absorbing boundary condition, in unbounded domains. That same method turns to be an efficient tool to make numerical zooms in regions of a particular interest. The TOS method enjoys, then, the ability to compute small structures one wants to capture and the reliability to obtain the...
The numerical solution of the elliptic Monge-Ampère Partial Differential Equation has been a subject of increasing interest recently [Glowinski, in 6th International Congress on Industrial and Applied Mathematics, ICIAM 07, Invited Lectures (2009) 155–192; Oliker and Prussner, Numer. Math.54 (1988) 271–293; Oberman, Discrete Contin. Dyn. Syst. Ser. B10 (2008) 221–238; Dean and Glowinski, in Partial differential equations, Comput. Methods Appl. Sci. 16 (2008) 43–63; Glowinski et al., Japan...
In this article, we present a new two-level stabilized nonconforming finite elements method for the two dimensional Stokes problem. This method is based on a local Gauss integration technique and the mixed nonconforming finite element of the pair (nonconforming linear element for the velocity, conforming linear element for the pressure). The two-level stabilized finite element method involves solving a small stabilized Stokes problem on a coarse mesh with mesh size and a large stabilized Stokes...
We derive an optimal lower bound of the interpolation error for linear finite elements on a bounded two-dimensional domain. Using the supercloseness between the linear interpolant of the true solution of an elliptic problem and its finite element solution on uniform partitions, we further obtain two-sided a priori bounds of the discretization error by means of the interpolation error. Two-sided bounds for bilinear finite elements are given as well. Numerical tests illustrate our theoretical analysis....
We derive an optimal lower bound of the interpolation error for linear finite elements on a bounded two-dimensional domain. Using the supercloseness between the linear interpolant of the true solution of an elliptic problem and its finite element solution on uniform partitions, we further obtain two-sided a priori bounds of the discretization error by means of the interpolation error. Two-sided bounds for bilinear finite elements are given as well. Numerical tests illustrate our theoretical analysis. ...
The paper concerns the solution of partial differential equations of evolution type by the finite difference method. The author discusses the general assumptions on the original equation as well as its discretization, which guarantee that the difference scheme is unconditionally stable, i.e. stable without any stability condition for the time-step. A new notion of the -acceptability of the integration formula is introduced and examples of such formulas are given. The results can be applied to ordinary...