Quality assessment and control of finite element solutions

Status and some recent developments in the techniques for assessing the reliability of finite element solutions are summarized. Discussion focuses on a number of aspects including: the major types of errors in the finite element solutions; techniques used for a posteriori error estimation and the reliability of these estimators; the feedback and adaptive strategies for improving the finite element solutions; and postprocessing approaches used for improving the accuracy of stresses and other important engineering data. Also, future directions for research needed to make error estimation and adaptive improvement practical are identified.     

Error estimates for shooting methods in two-point boundary value problems for second-order equations

This paper is concerned with a procedure for estimating the global discretization error arising when a boundary value problem for a system of second order differential equations is solved by the simple shooting method, without transforming the original problem in an equivalent first order problem. Expressions of the global discretization error are derived for both linear and nonlinear boundary value problems, which reduce the error estimation for a boundary value problem to that for an initial value problem of same dimension. The procedure extends to second order equations a technique for global error estimation given elsewhere for first order equations. As a practical result the accuracy of the estimates for a second order problem is increased compared with the estimates for the equivalent first order problem.     

Asymptotic expansions of finite element solutions to Robin problems in H 3 and their application in extrapolation cascadic multigrid method

For the Poisson equation with Robin boundary conditions,by using a few techniques such as orthogonal expansion(M-type),separation of the main part and the finite element projection,we prove for the first time that the asymptotic error expansions of bilinear finite element have the accuracy of O(h3)for u∈H3.Based on the obtained asymptotic error expansions for linear finite elements,extrapolation cascadic multigrid method(EXCMG)can be used to solve Robin problems effectively.Furthermore,by virtue of Richardson not only the accuracy of the approximation is improved,but also a posteriori error estimation is obtained.Finally,some numerical experiments that confirm the theoretical analysis are presented.     

双相滞热传导方程的有限元分析

本文考虑一类具有广泛应用背景的双相滞热传导方程混合边界问题．建立了其有限元和交替方向有限元的两种数值逼近格式．利用微分方程的先验估计理论与技巧,作出了数值解的L^2—范数估计结果．基于一系列的误差估计,也研究了两种逼近格式数值的稳定性和收敛性。     

A posteriori finite element error estimators for indefinite elliptic boundary value problems ∗

A general construction technique is presented for a posteriori error estimators of finite element solutions of elliptic boundary value problems that satisfy a Gång inequality. The estimators are obtained by an element–by–element solution of ‘weak residual’ with or without considering element boundary residuals. There is no order restriction on the finite element spaces used for the approximate solution or the error estimation; that is, the design of the estimators is applicable in connection with either one of the h–p–, or hp– formulations of the finite element method. Under suitable assumptions it is shown that the estimators are bounded by constant multiples of the true error in a suitable norm. Some numerical results are given to demonstrate the effectiveness and efficiency of the approach.     

Bandwidth choice for local polynomial estimation of smooth boundaries

Local polynomial methods hold considerable promise for boundary estimation, where they offer unmatched flexibility and adaptivity. Most rival techniques provide only a single order of approximation; local polynomial approaches allow any order desired. Their more conventional rivals, for example high-order kernel methods in the context of regression, do not have attractive versions in the case of boundary estimation. However, the adoption of local polynomial methods for boundary estimation is inhibited by lack of knowledge about their properties, in particular about the manner in which they are influenced by bandwidth; and by the absence of techniques for empirical bandwidth choice. In the present paper we detail the way in which bandwidth selection determines mean squared error of local polynomial boundary estimators, showing that it is substantially more complex than in regression settings. For example, asymptotic formulae for bias and variance contributions to mean squared error no longer decompose into monotone functions of bandwidth. Nevertheless, once these properties are understood, relatively simple empirical bandwidth selection methods can be developed. We suggest a new approach to both local and global bandwidth choice, and describe its properties.     

Estimates of modeling errors for linear elliptic problems

We derive new estimates of modeling errors for linear elliptic boundary value problems with periodic coefficients solved by homogenization method. Our approach is based on the concept of functional a posteriori error estimation. The estimates are obtained for the energy norm and use solely the global flux of the non-oscillatory solution of the homogenized model and solution of a boundary value problem on the cell of periodicity. (© 2014 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

On a boundary value problem for delay differential equations of population dynamics and Chebyshev approximation

In this paper, a boundary value problem for delay differential equations of population dynamics is considered. We obtain approximate solutions by using Chebyshev polynomial series and Newton–Raphson's procedure and give the error estimation. The method of the error estimation has been obtained in an existence theorem proved by a part of the authors. We carry out some numerical experiments by a computer language MATLAB.     

Mesh selection for stiff two-point boundary value problems

This paper is concerned with the mesh selection algorithm of COLSYS, a well known collocation code for solving systems of boundary value problems. COLSYS was originally designed to solve non-stiff and mildly stiff problems only. In this paper we show that its performance for solving extremely stiff problems can be considerably improved by modifying its error estimation and mesh selection algorithms. Numerical examples indicate the superiority of the modified algorithm.Dedicated to John Butcher on the occasion of his sixtieth birthday     

A new error estimation approach based on fuzzy logic system for H-R adaptive boundary element method

Conventional adaptive boundary element method cannot be universally applied to solve many more problems than the subject it discussed, and different error estimation formulas need to be designed for varied problems. This paper put forward a new error analysis method based on the fuzzy logic system, which is able to make error estimation effectively using human expert experience, and solve the two classical elasticity problems in conjunction with the H-R adaptive boundary element method. Numerical examples have illustrated the effectiveness, superiority and potential of a fuzzy logic approach in the adaptive boundary element method.     

一类非线性向量微分方程无穷边值问题的奇摄动

讨论一类非线性向量微分方程无穷边值问题的奇摄动,虽然此类边值问题在有限区间上曾被广泛地研究过,但在无限区间上还未曾采用对角化的方法去研究它.在适当的条件下,该文采用新的方法去证明解的存在性和任意阶的一致有效的渐近展开式,同时也给出误差估计.     

$C^1$ error estimation on the boundary for an exterior Neumann problem in $\mathbb R^3$

Summary.   In this paper we establish a error estimation on the boundary for the solution of an exterior Neumann problem in . To solve this problem we consider an integral representation which depends from the solution of a boundary integral equation. We use a full piecewise linear discretisation which on one hand leads to a simple numerical algorithm but on the other hand the error analysis becomes more difficult due to the singularity of the integral kernel. We construct a particular approximation for the solution of the boundary integral equation, for the solution of the Neumann problem and its gradient on the boundary and estimate their error. Received May 11, 1998 / Revised version received July 7, 1999 / Published online August 24, 2000     

Boundary output feedback stabilization for a class of coupled parabolic systems

In this paper, the boundary output feedback stabilization problem is addressed for a class of coupled nonlinear parabolic systems. An output feedback controller is presented by introducing a Luenberger‐type observer based on the measured outputs. To determine observer gains, a backstepping transform is introduced by choosing a suitable target system with nonlinearity. Furthermore, based on the state observer, a backstepping boundary control scheme is presented. With rigorous analysis, it is proved that the states of nonlinear closed‐loop system including state estimation and estimation error of plant system are locally exponentially stable in the L²norm. Finally, a numerical example is proposed to illustrate the effectiveness of the presented scheme.     

A posteriori error majorants of the modeling errors for elliptic homogenization problems

In this paper, we derive new two-sided a posteriori estimates of the modeling errors for linear elliptic boundary value problems with periodic coefficients solved by homogenization. Our approach is based on the concept of functional a posteriori error estimation. The estimates are obtained for the energy norm and use solely the global flux of the non-oscillatory solution of the homogenized model and solution of a boundary value problem on the cell of periodicity.     

Time step selection for the numerical solution of boundary value problems for parabolic equations

An algorithm is proposed for selecting a time step for the numerical solution of boundary value problems for parabolic equations. The solution is found by applying unconditionally stable implicit schemes, while the time step is selected using the solution produced by an explicit scheme. Explicit computational formulas are based on truncation error estimation at a new time level. Numerical results for a model parabolic boundary value problem are presented, which demonstrate the performance of the time step selection algorithm.     

屈曲特征值问题的边界元方法及收敛性分析

讨论了屈曲特征值问题的定解条件,建立了相应的具约束的积分方程组及带Lagrange乘子的边界变分方程,给出了解的存在唯一性定理。建立了相应的边界元方法并讨论了近似解的误差估计。文末给出了数值算例。     

Free-edge stress analysis of general composite laminates under extension, torsion and bending

In this study, based on the reduced form of elasticity displacement field for a long laminate, an analytical method is established to exactly obtain the interlaminar stresses near the free edges of generally laminated composite plates subjects to extension, torsion, and bending. The constant parameters being in the displacement field, which describe the global deformation of a laminate, are appropriately calculated by using the improved first-order shear deformation theory. Reddy’s layerwise theory is subsequently employed for analytical and numerical examinations of the boundary layer stresses within arbitrary laminated composite plates. Various numerical results are developed for the interlaminar normal and shear stresses along the interfaces and through the thickness of laminates near the free edges. Finally the effects of end conditions of laminates and geometric parameters on the boundary-layer stress are studied.     

Solution for a problem of linear plane elasticity with mixed boundary conditions on an ellipse by the method of boundary integrals

A numerical boundary integral scheme is proposed for the solution of the system of field equations of plane, linear elasticity in stresses for homogeneous, isotropic media in the domain bounded by an ellipse under mixed boundary conditions. The stresses are prescribed on one half of the ellipse, while the displacements are given on the other half. The method relies on previous analytical work within the Boundary Integral Method [1], [2].The considered problem with mixed boundary conditions is replaced by two subproblems with homogeneous boundary conditions, one of each type, having a common solution. The equations are reduced to a system of boundary integral equations, which is then discretized in the usual way and the problem at this stage is reduced to the solution of a rectangular linear system of algebraic equations. The unknowns in this system of equations are the boundary values of four harmonic functions which define the full elastic solution inside the domain, and the unknown boundary values of stresses or displacements on proper parts of the boundary.On the basis of the obtained results, it is inferred that the tangential stress component on the fixed part of the boundary has a singularity at each of the two separation points, thought to be of logarithmic type. A tentative form for the singular solution is proposed to calculate the full solution in bulk directly from the given boundary conditions using the well-known Boundary Collocation Method. It is shown that this addition substantially decreases the error in satisfying the boundary conditions on some interval not containing the singular points.The obtained results are discussed and boundary curves for unknown functions are provided, as well as three-dimensional plots for quantities of practical interest. The efficiency of the used numerical schemes is discussed, in what concerns the number of boundary nodes needed to calculate the approximate solution.     

A Dual Coupled Method for Boundary Value Problems of PDE with Coefficients of Small Period

1.IntroductionThemechanicalperformanceanalysisofthestructuresmadeofwovncompositematerialandperiodicallyperforatedmaterialisoftenencounteredinthemodernengi-neeringanalysis.Sincethiskindofcompositematerialhasperiodicallybasicconfigu-rations,thestaticanalysisofthestructuresmadefromthiscompositematerialleadstotheboundaryvalueproblemofellipticPDEwithperiodiccoefficients,forexample,theequllibriumproblemofwovenmembraneundertraverseloadingscanbeexpressedintheboundaryvalueproblemoftwodimensiontwoorder…     

Error estimates for a discretized Galerkin method for a boundary integral equation in two dimensions

We present a priori and a posteriori estimates for the error between the Galerkin and a discretized Galerkin method for the boundary integral equation for the single layer potential on the square plate. Using piecewise constant finite elements on a rectangular mesh we study the error coming from numerical integration. The crucial point of our analysis is the estimation of some error constants, and we demonstrate that this is necessary if our methods are to be used. After the determination of these constants we are in the position to prove invertibility and quasioptimal convergence results for our numerical scheme, if the chosen numerical integration formulas are sufficiently precise. © 1992 John Wiley & Sons, Inc.