composition operators on hardy spaces on lavrentiev domains in this note composition operators on bergman spaces of a simply connected domain are studied characterizing boundedness and compactness of such operators whenever the domain is lavrentiev

a unified approach to a posteriori error estimation using element residual methods this paper deals with the problem of obtaining numerical estimates of the accuracy of approximations to solutions of elliptic partial differential equations it is shown that by solving appropriate local residual type problems one can obtain upper bounds on the error in the energy norm moreover in the special case of adaptiveh p finite element analysis the estimator will also give a realistic estimate of the error a key feature of this is the development of a systematic approach to the determination of boundary conditions for the local problems the work extends and combines several existing methods to the case of fullh p finite element approximation on possibly irregular meshes with elements of non uniform degree as a special case the analysis proves a conjecture made by bank and weiser some a posteriori error estimators for elliptic partial differential equations math comput 44 283 301 1985

a simple error estimator and adaptive procedure for practical engineerng analysis a new error estimator is presented which is not only reasonably accurate but whose evaluation is computationally so simple that it can be readily implemented in existing finite element codes the estimator allows the global energy norm error to be well estmated and alos gives a good evaluation of local errors it can thus be combined with a full adaptive process of refinement or more simply provide guidance for mesh redesign which allows the user to obtain a desired accuracy with one or two trials when combined with an automatic mesh generator a very efficient guidance process to analysis is avaiable estimates other than the energy norm have successfully been applied giving for instance a predetermined accuracy of stresses

a posteriori error estimation of h p finite element approximations of frictional contact problems abstract dynamic and static frictional contact problems are described using the normal compliance law on the contact boundary dynamic problems are recast into quasistatic problems by time discretization an a posteriori error estimator is developed for a nonlinear elliptic equation of corresponding static or quasistatic problems the a posteriori error estimator is applied to a frictionless case and extended to frictional contact problems also an adaptive strategy is introduced and h p finite element meshes are obtained through a procedure based on a priori and a posteriori error estimations numerical examples are given to support the theoretical results

a posteriori pointwise error computation for 2 d transport equations based on the variational multiscale method abstract this article presents a general framework to estimate the pointwise error of linear partial differential equations the error estimator is based on the variational multiscale theory in which the error is decomposed in two components according to the nature of the residuals element interior residuals and inter element jumps the relationship between the residuals coarse scales and the error components fine scales is established yielding to a very simple model in particular the pointwise error is modeled as a linear combination of bubble functions and green s functions if residual free bubbles and the classical green s function are employed the technology leads to an exact explicit method for the pointwise error if bubble functions and free space green s functions are employed then a local projection problem must be solved within each element and a global boundary integral equation must be solved on the domain boundary as a consequence this gives a model for the so called fine scale green s functions the numerical error is studied for the standard galerkin and supg methods with application to the heat equation the reaction diffusion equation and the convection diffusion equation numerical results show that stabilized methods minimize the propagation of pollution errors which stay mostly locally

recovering lower bounds of the error by postprocessing implicit residual a posteriori error estimates classical residual type error estimators approximate the error flux around the elements and yield upper bounds of the exact or reference error lower bounds of the error are also needed in goal oriented adaptivity and for bounds on functional outputs this work introduces a simple and cheap strategy to recover a lower bound estimate from standard upper bound estimates this lower bound may also be used to assess the effectivity of the former estimate and to improve it

toward a universal h p adaptive finite element strategy part 2 a posteriori error estimation abstract in this second paper of the series on h p adaptive finite element methods we present five different techniques for a posteriori error estimation of h p methods for linear elliptic boundary value problems the techniques are referred to as the residual estimation method the duality method the subdomain residual method a method based on interpolation theory and a post processing method the present study examines the theory and behavior of such techniques applied specifically to h p adaptive finite element methods several numerical examples illustrate the efficiency of the methods

adaptive mesh refinement in strain localization problems an adaptive meshing method tailored to problems of strain localization is given the adaption strategy consists of equi distributing the variation of the velocity field over the elements of the mesh a heuristic justification for the use of variations as indicators is advanced and possible connections with interpolation error bounds are discussed meshes are constructed by delaunay triangulation it is shown how the hu washizu principle determines a consistent transfer operator for the state variables examples of application are given which demonstrate the versatility of the method

enrichment and coupling of the finite element and meshless methods a mixed hierarchical approximation based on finite elements and meshless methods is presented two cases are considered the first one couples regions where finite elements or meshless methods are used to interpolate continuity and consistency is preserved the second one enriches a finite element mesh with particles thus there is no need to remesh in adaptive refinement processes in both cases the same formulation is used convergence is studied and examples are shown copyright 2000 john wiley sons ltd

analysis of the efficiency of an a posteriori error estimator for linear triangular finite elements this paper adresses the problem of determining upper and lower bounds for the effectivity index on the a posteriori estimate of the error in the finite element method these bounds are given explicitly for a certain concrete estimator for linear elements and unstructured triangular meshes they depend strongly on the geometry of the triangles and relatively weakly on the smoothness of the solution an example shows that the bounds are not over pessimistic in babuska plank and rodriguez 4 detailed numerical experimentation is given

the multiscale approach to error estimation and adaptivity abstract a new explicit a posteriori error estimator is investigated which emanates from the variational multiscale theory the error estimator uses an approximation of the green s function which propagates the error according to the dual problem the technique predicts the error in the l2 norm from element residuals both interior and inter element jumps therefore leading to a very economical procedure finally the methodology is applied to fluid flow transport showing that for convection dominated flows the efficiency index is independent of the diffusion coefficient

a posteriori error estimation for standard finite element analysis a new residual type estimator based on projections of the error on subspaces of locally supported functions is presented the estimator is defined by a standard element by element refinement first an approximation of the energy norm of the error is obtained solving local problems with homogeneous dirichlet boundary conditions a later enrichment of the estimation is performed by adding the contributions of projections on a new family of subspaces this estimate is a lower bound of the measure of the actual error the estimator does not need to approximate local boundary conditions for the error equation therefore computation of flux jumps is not necessary moreover the estimator can be applied in mixed meshes containing elements of different shapes and its implementation in a standard finite element code is straightforward the presented results show the effectiveness of the estimator approximating both the distribution and the global measure of the error as well as its usefulness in adaptive procedures

a unified approach to remeshing strategies for finite element h adaptivity in h adaptivity a remeshing strategy is needed to compute the distribution of required element size using the estimated error distribution several authors have introduced different remeshing strategies yielding very different results in this work these methods are included in a unified framework emphasizing the role of the underlying hypotheses moreover an objective tool to evaluate the accuracy of the resulting finite element solution is presented thus a new remeshing strategy is introduced to optimize the accuracy of the adapted solutions the different remeshing strategies are compared with well known numerical examples

adaptive fem analysis of selected elastic visco plastic problems abstract elastic visco plastic deformations modeled by the bodner partom constitutive equations s r bodner y partom constitutive equations elastic viscoplastic strain hardening materials j appl mech 42 1975 385 389 have been evaluated by the adaptive finite element method j m bass j t oden adaptive finite element methods for a class of evolution problems in viscoplasticity int j engrg sci 25 623 653 1987 l demkowicz j t oden w rachowicz o hardy toward a universal hp adaptive finite element strategy part 1 constrained approximation and data structure comput methods appl mech engrg 77 1989 79 112 we use the explicit residual error estimate i babuska w c rheinboldt error estimates for adaptive finite element computations int j numer methods engrg 12 1978 1597 1615 supplemented with a special term accounting for the fact that the right hand side of the equations we use is only approximated since it depends on the inelastic strains which are computed only at the gauss integration points we also address the problems of adaptation strategy and the method of solution transfer to a new mesh results of recent tests confirm efficiency and robustness of the approach

new advances on a posteriori error on constitutive relation in f e analysis this paper deals with error estimators based on residuals on the constitutive relation which have been developed for the past 20 years at cachan especially in the case of elasticity these error estimators can become poor in situations where the material and the mesh are strongly anisotropic the crucial point herein is the construction of equilibrated stress fields from the computed finite element solution therefore in this paper the approach and in particular this last point are re examined several important modifications are introduced in the concepts and basic techniques hence a new generation of error estimators which seem to be extremely robust upon initial applications has been generated

a posteriori estimation and adaptive control of the pollution error in the h version of the finite element method abstract we studied the pollution error in the h version of the finite element method and its effect on the quality of the local error indicators resp the quality of the derivatives recovered by local postprocessing in the interior of the mesh here we show that it is possible to construct a posteriori estimates of the pollution error in a patch of elements by employing the local error indicators over the entire mesh we also give an adaptive algorithm for the local control of the pollution error in a patch of elements of interest

adapting broyden method to handle linear constraints imposed via lagrange multipliers various non linear equation solvers are adapted to handle linear constraints via the lagrange multiplier technique this adaptation process turns out to be quite straightforward for newton raphson methods and rank two quasi newton methods bfgs and dfp but rather more involved for broyden method in fact two broyden methods can be obtained the standard one and a modified one better adapted to the lagrange multiplier environment some numerical examples are used to assess the relative performance of the various adapted solvers these tests illustrate the superiority of the modified broyden method over the standard one

bounds of functional outputs for parabolic problems part i exact bounds of the discontinuous galerkin time discretization classical implicit residual type error estimators require using an underlying spatial finer mesh to compute bounds for some quantity of interest consequently the bounds obtained are only guaranteed asymptotically that is with respect to the reference solution computed with the fine mesh exact bounds that is bounds guaranteed with respect to the exact solution are needed to properly certify the accuracy of the results especially if the meshes are coarse the paper introduces a procedure to compute strict upper and lower bounds of the error in linear functional outputs of parabolic problems in this first part the bounds account for the error associated with the spatial discretization the error coming from the time marching scheme is therefore assumed to be negligible in front of the spatial error the time discretization is performed using the discontinuous galerkin method both for the primal and adjoint problems in the error estimation procedure equilibrated fluxes at interelement edges are calculated using hybridization techniques

some a posteriori error estimators for elliptic partial differential equations we present three new a posteriori error estimators in the energy norm for finite element solutions to elliptic partial differential equations the estimators are based on solving local neumann problems in each element the estimators differ in how they enforce consistency of the neumann problems we prove that as the mesh size decreases under suitable assumptions two of the error estimators approach upper bounds on the norm of the true error and all three error estimators are within multiplicative constants of the norm of the true error we present numerical results in which one of the error estimators appears to converge to the norm of the true error

a posteriori error estimation and adaptivity for linear elasticity using the reciprocal theorem the reciprocal theorem of betti and rayleigh or in other notions duality arguments are used to derive error estimators for the finite element approximation of various quantities including local variables like single displacements and stresses the proposed error estimator is evaluated solving the set of equations for an additional right hand side simply applying the energy norm error estimators two times furthermore a general h adaptive algorithm is introduced which allows us to optimize meshes with respect to different user specified variables the efficiency of the current approach is demonstrated for plate and shell examples

the problem of the selection of an a posteriori error indicator based on smoothening techniques abstract this paper addresses the problem of assessing the quality of an a posteriori error estimate of a finite element solution an error estimate based on local l projections is analyzed in the case of translation invariant meshes it is shown that for general meshes this technique does not lead to an asymptotically exact estimator the problem is analyzed in detail in the one dimensional setting it is shown that an asymptotically exact estimator is not the optimal one when the solution is not sufficiently smooth an optimal estimator for adaptively constructed meshes is given finally a general mathematical framework for the quality assessment of estimators is introduced

a posteriori error estimation the relationship between different procedures to address the issue of accurate and efficient estimation of the discretization error of the finite element solution the relationship between different types of error estimators is investigated we found that some of the widely used error estimators are in fact related although they were originally derived from different methodologies the recognition of the relation between these error estimators provides a different viewpoint on the theory of a posteriori error estimation and a guidance in the development of the error estimators

subdomain based flux free a posteriori error estimators a new residual type flux free error estimator is presented it estimates upper and lower bounds of the error in energy norm the proposed approach precludes the main drawbacks of standard residual type estimators circumvents the need of flux equilibration and results in a simple implementation that uses standard resources available in finite element codes this is specially interesting for 3d applications where the implementation of this technique is as simple as in 2d recall that on the contrary the complexity of the flux equilibration techniques increases drastically in the 3d case the bounds for the energy norm of the error are used to produce upper and lower bounds of linear functional outputs representing quantities of engineering interest the presented estimators demonstrate their efficiency in numerical tests producing sharp estimates both for the energy and the quantities of interest

error estimation and adaptivity for discontinuous failure discontinuous failure is simulated via the introduction of a geometrical discontinuity the cohesive zone is modelled via the use of the partition of unity property of the finite element interpolation by this approach a crack can pass through the elements without any restriction to the underlying mesh despite such a feature it has been confirmed that a sufficiently fine mesh discretization still needs to be ensured in order to obtain a correct crack path and mechanical response the p adaptive scheme driven by error in an energy norm measure or a goal oriented measure has been examined due to its implementational simplicity the results have shown that if considering only increasing the polynomial degree the p approach can greatly improve the results copyright 2008 john wiley sons ltd

error estimation and adaptivity for nonlocal damage models abstract nonlocal damage models are typically used to model failure of quasi brittle materials due to brittleness the choice of a particular model or set of parameters can have a crucial influence on the structural response to assess this influence it is essential to keep finite element discretization errors under control if not the effect of these errors on the result of a computation could be erroneously interpreted from a constitutive viewpoint to ensure the quality of the fe solution an adaptive strategy based on error estimation is proposed here it is based on the combination of a residual type error estimator and quadrilateral h remeshing another important consequence of brittleness is that it leads to structural responses of the snap through or snap back type this requires the use of arc length control with a definition of the arc parameter that accounts for the localized nature of quasi brittle failure all these aspects are discussed for two particular nonlocal damage models mazars and modified von mises and for two tests the brazilian tensile splitting test and the single edge notched beam test for the latter test the capability of the mazars model to capture the curved crack pattern observed in experiments a topic of debate in the literature is confirmed

pollution error in the h version of the finite element method and the local quality of the recovered derivatives abstract in this paper we address the quality of the solution derivatives which are recovered from finite element solutions by local averaging schemes as an example we consider the zienkiewicz zhu superconvergent patch recovery scheme the zz spr scheme and we study its accuracy in the interior of the mesh for finite element approximations of solutions of laplace s equation in polygonal domains we will demonstrate the following 1 in general the accuracy of the derivatives recovered by the zz spr or any other local averaging scheme may not be higher than the accuracy of the derivatives computed directly from the finite element solution 2 if the mesh is globally adaptive i e it is nearly equilibrated in the energy norm then we can practically always gain in accuracy by employing the recovered derivatives instead of the derivatives computed directly from the finite element solution 3 it is possible to guarantee that the recovered solution derivatives have higher accuracy than the derivatives computed directly from the finite element solution in any patch of elements of interest by employing a mesh which is adaptive only with respect to the patch of interest i e it is nearly equilibrated in a weighted energy norm 4 in practice we are often interested in obtaining highly accurate derivatives or heat fluxes stresses etc only in a few critical regions which are identified by a preliminary analysis a grid which is adaptive only with respect to the critical regions of interest may be much more economical for this purpose because it may achieve the desired accuracy by employing substantially fewer degrees of freedom than a globally adaptive grid which achieves comparable accuracy in the critical regions

a variational multiscale a posteriori error estimation method for mixed form of nearly incompressible elasticity this paper presents an error estimation framework for a mixed displacement pressure finite element method for nearly incompressible elasticity the proposed method is based on variational multiscale vms concepts wherein the displacement field is decomposed into coarse scales that can be resolved by a given finite element mesh and fine scales that are beyond the resolution capacity of the mesh variational projection of fine scales onto the coarse scale space via variational embedding of the fine scale solution into the coarse scale formulation leads to the stabilized method with two major attributes first it is free of volumetric locking and second it accommodates arbitrary combinations of interpolation functions for the displacement and pressure fields this vms based stabilized method is equipped with naturally derived error estimators and offers various options for numerical computation of the error specifically two error estimators are explored the first method employs an element based strategy and a representation of error via a fine scale error equation defined over element interiors which is evaluated by a direct post solution evaluation this quantity when combined with the global pollution error results in a simple explicit error estimator the second method involves solving the fine scale error equation through localization to overlapping patches spread across the domain thereby leading to an implicit calculation of the local error this implicit calculation when combined with the global pollution error results in an implicit error estimator the performance of the stabilized method and the error estimators is investigated through numerical convergence tests conducted for two model problems on uniform and distorted meshes the sharpness and robustness of the estimators is shown to be consistent across the test cases employed

a posteriori estimation of the error in the recovered derivatives of the finite element solution this work addresses the accuracy of the solution derivatives which are recovered by local averaging of the finite element solution the main results of the study are 1 the error in the locally averaged derivatives e g the derivatives which are recovered by the zienkiewicz zhu superconvergent patch recovery zz spr or other similar local recoveries can be more than the error in the derivatives computed directly from the finite element solution especially in the case of unsmooth solutions and or coarse meshes 2 in order to determine which solution derivatives should be relied upon the locally averaged ones or the ones computed directly from the finite element solution one must be able to estimate their errors it is shown that one can obtain indicators of the error in the derivatives recovered by the zz spr by employing an additional local averaging of the recovered derivatives recycling of the zz spr or by comparing the derivatives computed by the zz spr with the derivatives obtained using a different local averaging which takes into account the character of the exact solution harmonic averaging

error estimation including pollution assessment for nonlinear finite element analysis a residual type error estimator for nonlinear finite element analysis is introduced this error estimator solves local problems avoiding both the computation of the flux jumps and the associated flux splitting procedure pollution errors are taken into account by a feedback strategy that is an error estimate based on local computations is used as the input of the pollution analysis this estimator is used in the frame of an adaptive procedure numerical examples show that the estimator is able to drive adaptive procedures leading to likely good solutions moreover one of the examples demonstrates that adaptive procedures are essential for complex highly nonlinear mechanical problems because they may discover secondary collapse mechanisms

the background of error estimation and adaptivity in finite element computations the paper summarizes the present state of the art concerning error estimation and adaptive remeshing recent work shows that methods based on recovery result on most accurate error estimators

the discrete origin of fetd newmark late time instability and a correction scheme the finite element time domain fetd method is commonly used for transient simulation of electromagnetic wave phenomena most practitioners consider fetd when time integrated using the newmark beta method to be unconditionally stable when s 0 25 unlike the finite difference time domain fdtd courant criterion fetd newmark has no limiting timestep above which the method exhibits exponential growth however herein the stability properties of fetd newmark will be rigorously investigated by deducing the jordan canonical form of the fetd newmark amplification matrix and it will be demonstrated that the method does exhibit linear growth for certain field configurations these modes are none other than the pure gradient fields associated with late time instability though many practical simulations are of short duration and will never observe a linearly growing gradient solution it can be problematic for simulations which require long time periods to be integrated a correction scheme for eliminating this late time instability shall be suggested and numerical results will verify its performance