A posteriori error analysis of an augmented mixed-primal formulation for the stationary Boussinesq model

In an earlier work of us, a new mixed finite element scheme was developed for the Boussinesq model describing natural convection. Our methodology consisted of a fixed-point strategy for the variational problem that resulted after introducing a modified pseudostress tensor and the normal component of the temperature gradient as auxiliary unknowns in the corresponding Navier-Stokes and advection-diffusion equations defining the model, respectively, along with the incorporation of parameterized redundant Galerkin terms. The well-posedness of both the continuous and discrete settings, the convergence of the associated Galerkin scheme, as well as a priori error estimates of optimal order were stated there. In this work we complement the numerical analysis of our aforementioned augmented mixed-primal method by carrying out a corresponding a posteriori error estimation in two and three dimensions. Standard arguments relying on duality techniques, and suitable Helmholtz decompositions are used to derive a global error indicator and to show its reliability. A globally efficiency property with respect to the natural norm is further proved via usual localization techniques of bubble functions. Finally, an adaptive algorithm based on a reliable, fully local and computable a posteriori error estimator induced by the aforementioned one is proposed, and its performance and effectiveness are illustrated through a few numerical examples in two dimensions.     

Formulation and analysis of fully-mixed methods for stress-assisted diffusion problems

This paper is devoted to the mathematical and numerical analysis of a mixed-mixed PDE system describing the stress-assisted diffusion of a solute into an elastic material. The equations of elastostatics are written in mixed form using stress, rotation and displacements, whereas the diffusion equation is also set in a mixed three-field form, solving for the solute concentration, for its gradient, and for the diffusive flux. This setting simplifies the treatment of the nonlinearity in the stress-assisted diffusion term. The analysis of existence and uniqueness of weak solutions to the coupled problem follows as combination of Schauder and Banach fixed-point theorems together with the Babu?ka–Brezzi and Lax–Milgram theories. Concerning numerical discretization, we propose two families of finite element methods, based on either PEERS or Arnold–Falk–Winther elements for elasticity, and a Raviart–Thomas and piecewise polynomial triplet approximating the mixed diffusion equation. We prove the well-posedness of the discrete problems, and derive optimal error bounds using a Strang inequality. We further confirm the accuracy and performance of our methods through computational tests.     

A posteriori error estimates of two-grid finite volume element methods for nonlinear elliptic problems

In this paper, we study the a posteriori error estimates of two-grid finite volume element method for second-order nonlinear elliptic equations. We derive the residual-based a posteriori error estimator and prove the computable upper and lower bounds on the error in $H^1$-norm. The a posteriori error estimator can be used to assess the accuracy of the two-grid finite volume element solutions in practical applications. Numerical examples are provided to illustrate the performance of the proposed estimator.     

A posteriori error estimates of finite element method for the time-dependent Darcy problem in an axisymmetric domain

We consider the time dependent Darcy problem in a three-dimensional axisymmetric domain and, by writing the Fourier expansion of its solution with respect to the angular variable, we observe that each Fourier coefficient satisfies a system of equations on the meridian domain. We propose a discretization of these equations in the case of general solution. This discretization relies on a backward Euler’s scheme for the time variable and finite elements for the space variables. We prove a posteriori error estimates that allow for an efficient adaptivity strategy both for the time steps and the meshes. Computations for an example with a known solution are presented which support the a posteriori error estimate.     

A priori and a posteriori error analyses of an HDG method for the Brinkman problem

In this paper we introduce and analyze a hybridizable discontinuous Galerkin (HDG) method for the linear Brinkman model of porous media flow in two and three dimensions and with non-homogeneous Dirichlet boundary conditions. We consider a fully-mixed formulation in which the main unknowns are given by the pseudostress, the velocity and the trace of the velocity, whereas the pressure is easily recovered through a simple postprocessing. We show that the corresponding continuous and discrete schemes are well-posed. In particular, we use the projection-based error analysis in order to derive a priori error estimates. Furthermore, we develop a reliable and efficient residual-based a posteriori error estimator, and propose the associated adaptive algorithm for our HDG approximation. Finally, several numerical results illustrating the performance of the method, confirming the theoretical properties of the estimator and showing the expected behavior of the adaptive refinements are presented.     

A posteriori error estimates for nonconforming approximations of Steklov eigenvalue problems

This paper deals with a posteriori error estimators for the non conforming Crouzeix-Raviart finite element approximations of the Steklov eigenvalue problem. First, we define an error estimator of the residual type which can be computed locally from the approximate eigenpair and we prove the equivalence between this estimator and the broken energy norm of the error with constants independent of the corresponding eigenvalue. Next, we prove that edge residuals dominate the volumetric part of the residual and that the volumetric part of the residual terms dominate the normal component of the jumps of the discrete fluxes across interior edges. Finally, based on these results, we introduce two simpler equivalent error estimators. The analysis shows that these a posteriori error estimates are optimal up to higher order terms and that may be used for the design of adaptive algorithms.     

A posteriori error analysis of an augmented mixed formulation in linear elasticity with mixed and Dirichlet boundary conditions

We develop a residual-based a posteriori error analysis for the augmented mixed methods introduced in [13], [14] for the problem of linear elasticity in the plane. We prove that the proposed a posteriori error estimators are both reliable and efficient. Numerical experiments confirm these theoretical properties and illustrate the ability of the corresponding adaptive algorithms to localize the singularities and large stress regions of the solutions.     

A robust a posteriori error estimate for the Fortin-Soulie finite-element method

The subject of a posteriori error estimation is widely studied, and a variety of such error estimates have been used for elasticity problems in recent years. Of particular interest is the work carried out in ^{1. and 2.}. In this paper, we derive a new a posteriori error estimator for the quadratic nonconforming Fortin-Soulie element for the error in an energy-like norm. Then, we illustrate the new error bound by presenting some numerical examples, and show an example of a sequence of adaptively refined meshes.     

Simple a posteriori error estimators for the <Emphasis Type="Italic">h</Emphasis>-version of the boundary element method

The h-h/2-strategy is one well-known technique for the a posteriori error estimation for Galerkin discretizations of energy minimization problems. One considers to estimate the error , where is a Galerkin solution with respect to a mesh and is a Galerkin solution with respect to the mesh obtained from a uniform refinement of . This error estimator is always efficient and observed to be also reliable in practice. However, for boundary element methods, the energy norm is non-local and thus the error estimator η does not provide information for a local mesh-refinement. We consider Symm’s integral equation of the first kind, where the energy space is the negative-order Sobolev space . Recent localization techniques allow to replace the energy norm in this case by some weighted L ²-norm. Then, this very basic error estimation strategy is also applicable to steer an h-adaptive algorithm. Numerical experiments in 2D and 3D show that the proposed method works well in practice. A short conclusion is concerned with other integral equations, e.g., the hypersingular case with energy space and , respectively, or a transmission problem. Dedicated to Professor Ernst P. Stephan on the occasion of his 60th birthday.     

A priori and a posteriori error analysis of an augmented mixed finite element method for incompressible fluid flows

In this paper we extend recent results on the a priori and a posteriori error analysis of an augmented mixed finite element method for the linear elasticity problem, to the case of incompressible fluid flows with symmetric stress tensor. Similarly as before, the present approach is based on the introduction of the Galerkin least-squares type terms arising from the constitutive and equilibrium equations, and from the relations defining the pressure in terms of the stress tensor and the rotation in terms of the displacement, all of them multiplied by stabilization parameters. We show that these parameters can be suitably chosen so that the resulting augmented variational formulation is defined by a strongly coercive bilinear form, whence the associated Galerkin scheme becomes well-posed for any choice of finite element subspaces. Next, we present a reliable and efficient residual-based a posteriori error estimator for the augmented mixed finite element scheme. Finally, several numerical results confirming the theoretical properties of this estimator, and illustrating the capability of the corresponding adaptive algorithm to localize the singularities and the large stress regions of the solution, are reported.     

An adjustment for edge effects using an augmented neighborhood model in the spatial auto-logistic model

Parameter estimation in the spatial auto-regressive models has difficulty due to the edge sites which have unobserved neighborhood sites. Some ad hoc remedies suggested in the literature are the free boundary condition, the toroidal boundary condition or estimation using only internal data sites. However, parameter estimates are often sensitive to assumptions on the unobserved neighborhood sites and all the above assumptions have some apparent shortcomings such as systematic bias or inflated variance. In this paper, we propose a new way to incorporate the edge sites by introducing an augmented random neighborhood, denoted by the augmented neighborhood model, which represents the entire external field. To estimate the model, we derive the EM procedures for the maximum pseudo-likelihood estimator and the maximum likelihood estimator. Several simulation studies show that the random external field provides better performance of the maximum pseudo-likelihood estimator and the maximum likelihood estimator than other assumptions on the edge sites. As an example, we apply the random external field to modeling the distribution of Plantago lanceolata in Kansas.     

A priori and a posteriori error estimates for the method of lumped masses for parabolic optimal control problems

In this paper, we examine the method of lumped masses for the approximation of convex optimal control problems governed by linear parabolic equations, where the lumped mass method is used for the discretization of the state equation. We derive some a priori and a posteriori error estimates for both the state and control approximations with control constraints of obstacle type. Numerical experiments are given to show the efficiency and reliability of the lumped mass method.     

A finite element formulation for dynamic parameter identification of robot manipulators

In this paper a finite element based approach is described for the automatic generation of models suitable for dynamic parameter identification. The method involves a nonlinear finite element formulation in which both links and joints are considered as specific finite elements [6, 7]. Since the identification procedure considers rigid-link robot models, the inertial properties of the link elements are described using a lumped mass formulation. The parameters to be identified are masses, first-order moments and inertial tensor components of the links. The equations of motion are written in a form which is linear in the dynamic parameters. This formulation is obtained by employing Jourdain’s principle of virtual power. The parameters are estimated using a linear least squares technique. Singular value decomposition of the regression matrix is used to find the minimum parameter set. Simulation results obtained from the 6 DOF PUMA 560 robot based on the estimated parameters show that the method yields accurate responses.     

Superconvergence in a DPG method for an ultra-weak formulation

In this work we study a DPG method for an ultra-weak variational formulation of a reaction–diffusion problem. We improve existing a priori convergence results by sharpening an approximation result for the numerical flux. By duality arguments we show that higher convergence rates for the scalar field variable are obtained if the polynomial order of the corresponding approximation space is increased by one. Furthermore, we introduce a simple elementwise postprocessing of the solution and prove superconvergence. Numerical experiments indicate that the obtained results are valid beyond the underlying model problem.     

An adaptive mesh refinement of quadrilateral finite element meshes based upon an a posteriori error estimation of quantities of interest: modal response

Modal analysis is commonly performed in a vehicle development process to assess dynamic responses of structure designs. This paper presents an adaptive quadrilateral refinement process for modal analysis of elastic shells based upon a posteriori error estimation in natural frequencies. The process provides engineers with an estimation of their modal analysis quality and an effective adaptive refinement tool for quadrilateral meshes. The effectiveness of the process is demonstrated on the eigenvalue analyses of two numerical examples, a shock tower cap and a roof structure. It shows that the solution error in the frequency of interest is effectively reduced through the adaptive refinement process, and the resulting frequency of interest converges to the solution of a very fine model.     

Two guaranteed equilibrated error estimators for Harmonic formulations in eddy current problems

In this paper a guaranteed equilibrated error estimator is developed for the 3D harmonic magnetodynamic problem of Maxwell’s system. This system is recasted in the classical $\mathbf{A}?\phi$ potential formulation and solved by the Finite Element method. The error estimator is built starting from the $\mathbf{A}?\phi$ numerical solution by a local flux reconstruction technique. Its equivalence with the error in the energy norm is established. A comparison of this estimator with an equilibrated error estimator already developed through a complementary problem points out the advantages and drawbacks of these two estimators. In particular, an analytical benchmark test illustrates the obtained theoretical results and a physical benchmark test shows the efficiency of these two estimators.     

Numerical simulations of the Boussinesq equation by He's variational iteration method

In this paper, we present a reliable algorithm to study the well-known model of nonlinear dispersive waves proposed by Boussinesq. We solve the Cauchy problem of Boussinesq equation using variational iteration method (VIM). The numerical results of this method are compared with the exact solution of an artificial model to show the efficiency of the method. The approximate solutions show that VIM is a powerful mathematical tool for solving nonlinear problems.     

Boussinesq systems in two space dimensions over a variable bottom for the generation and propagation of tsunami waves

Considered here are Boussinesq systems of equations of surface water wave theory over a variable bottom. A simplified such Boussinesq system is derived and solved numerically by the standard Galerkin-finite element method. We study by numerical means the generation of tsunami waves due to bottom deformation and we compare the results with analytical solutions of the linearized Euler equations. Moreover, we study tsunami wave propagation in the case of the Java 2006 event, comparing the results of the Boussinesq model with those produced by the finite-difference code MOST, that solves the shallow water wave equations.     

An adaptive mesh refinement of quadrilateral finite element meshes based upon a posteriori error estimation of quantities of interest: linear static response

A posteriori error estimation in finite element analysis serves as an important guide to the meshing tool in an adaptive refinement process. However, the traditional posteriori error estimates, which are often defined in the energy or energy-type norms over the entire domain, provide users insufficient information regarding the accuracy of specific quantities in the solution. This paper describes an adaptive quadrilateral refinement process with a goal-oriented error estimation, in which a posteriori error is estimated with respect to the specified quantity of interest. A highlight of this paper is the demonstration of tools described in the paper used in a practical industrial environment. The performance of this process is demonstrated on several practical problems where the comparison is with the adaptive process based on the traditional error estimation.     

A Posteriori Error Estimates of Recovery Type for Distributed Convex Optimal Control Problems

In this paper, we derive a posteriori error estimates of recovery type, and present the superconvergence analysis for the finite element approximation of distributed convex optimal control problems. We provide a posteriori error estimates of recovery type for both the control and the state approximation, which are generally equivalent. Under some stronger assumptions, they are further shown to be asymptotically exact. Such estimates, which are apparently not available in the literature, can be used to construct adaptive finite element approximation schemes and as a reliability bound for the control problems. Numerical results demonstrating our theoretical results are also presented in this paper.