Search

Goal oriented error estimation and adaptive procedures are essential for the accurate and efficient evaluation of finite element numerical simulations that involve complex domains. By locally improving the approximation qual- ity, for example, by using the extended finite element method (XFEM), we can solve expensive problems which could result intractable otherwise. Here, we present an error estimation technique for enriched finite element approxi- mations that is based on an equilibrated recovery technique, which considers the stress intensity factor as the quantity of interest. The locally equilibrated superconvergent patch recovery is used to obtain enhanced stress fields for the primal and dual problems defined to evaluate the error estimate.

Some industrial processes are modelled by parametric partial differential equations. Integrating computational modelling and data assimilation into the control process requires obtaining a solution of the numerical model at the characteristic frequency of the process (real-time). This paper introduces a computational strategy allowing to efficiently exploit measurements of those industrial processes, providing the solution of the model at the required frequency. This is particularly interesting in the framework of control algorithms that rely on a model involving a set of parameters. For instance, the curing process of a composite material is modelled as a thermo-mechanical problem whose corresponding parameters describe the thermal and mechanical behaviours. In this context, the information available (measurements) is used to update the parameters of the model and to produce new values of the control variables (data assimilation). The methodology presented here is devised to ensure the possibility of having a response in real-time of the problem and therefore the capability of integrating it in the control scheme. The Proper Generalized Decomposition is used to describe the solution in the multi-parametric space. The real-time data assimilation requires a further simplification of the solution representation that better fits the data (reconstructed solution) and it provides an implicit parameter identification. Moreover, the analysis of the assimilated data sensibility with respect to the points where the measurements are taken suggests a criterion to locate of the sensors., Peer Reviewed, Postprint (author’s final draft)

Significant research effort has been devoted to produce one-sided error estimates for Finite Element Analyses, in particular to provide upper bounds of the actual error. Typically, this has been achieved using residual-type estimates. One of the most popular and simpler (in terms of implementation) techniques used in commercial codes is the recovery-based error estimator. This technique produces accurate estimations of the exact error but is not designed to naturally produce upper bounds of the error in energy norm. Some attempts to remedy this situation provide bounds depending on unknown constants. Here, a new step towards obtaining error bounds from the recovery-based estimates is proposed. The idea is (1) to use a locally equilibrated recovery technique to obtain an accurate estimation of the exact error, (2) to add an explicit-type error bound of the lack of equilibrium of the recovered stresses in order to guarantee a bound of the actual error and (3) to efficiently and accurately evaluate the constants appearing in the bounding expressions, thus providing asymptotic bounds. The numerical tests with h-adaptive refinement process show that the bounding property holds even for coarse meshes, providing upper bounds in practical applications. (C) 2015 Elsevier B.V. All rights reserved., Peer Reviewed, Postprint (author’s final draft)

Goal-oriented error estimates (GOEE) have become popular tools to quantify and control the local error in quantities of interest (QoI), which are often more pertinent than local errors in energy for design purposes (e.g. the mean stress or mean displacement in a particular area, the stress intensity factor for fracture problems). These GOEE are one of the key unsolved problems of advanced engineering applications in, for example, the aerospace industry. This work presents a simple recovery-based error estimation technique for QoIs whose main characteristic is the use of an enhanced version of the Superconvergent Patch Recovery (SPR) technique previously used for error estimation in the energy norm. This enhanced SPR technique is used to recover both the primal and dual solutions. It provides a nearly statically admissible stress field that results in accurate estimations of the local contributions to the discretisation error in the QoI and, therefore, in an accurate estimation of this magnitude. This approach leads to a technique with a reasonable computational cost that could easily be implemented into already available finite element codes, or as an independent postprocessing tool. © 2013 Springer-Verlag Berlin Heidelberg.

Over the last ten years there has been an increase on the use of goal-oriented error estimates aimed to quantify the local error on a (non)linear quantity of interest (QoI) that might result relevant for design purposes. Residual-based error estimators have been used recursively to obtain upper and lower bounds of the error in quantities of interest for finite element approximations. In this paper, we present a recovery technique for 2D linear elasticity problems, based on the superconvergent patch recovery (SPR), which provides recovered displacement and stress fields that are then utilised to evaluate practical upper and lower error bounds in QoI.

Traditionally, recovery based error estimators in linear elasticity have considered the evaluation of an enhanced stress field obtained from the raw Finite Element (FE) stress solution. Instead of that, one can also obtain a recovered displacement field from the FE displacements. Herein, we describe a superconvergent patch recovery of the displacement field which considers the local fulfilment of boundary and internal equilibrium equations, Dirichlet constraints and, for singular problems, the splitting of the displacement and stress fields into singular and smooth parts. Numerical tests using problems with known analytical solution have been carried out to validate the proposed technique for error estimation and error bounding in energy norm and quantities of interest.

In this work, an experimental technique is presented for the acoustic characterization of materials and devices. The developed technique is based on the two source method with simultaneous excitation, thus reducing the time requirements in comparison with the standard two source technique. A decomposition technique is applied to the pressure measurements associated with each microphone, that include noise of two sources acting simultaneously, and the input signals used to excite the two sources. The procedure allows the calculation of the transfer functions between each microphone and each input separately. Once the individual transfer functions have been obtained, it is possible to characterize an acoustic device by means of its four-pole matrix as well as its transmission loss. The technique has been validated by comparison with theoretical results and experimental methods of earlier studies, for several test cases including a straight duct and a reactive muffler., QC 20140910

This paper introduces a recovery-type error estimator yielding upper bounds of the error in energy norm for linear elastic fracture mechanics problems solved using the extended finite element method (XFEM). The paper can be considered as an extension and enhancement of a previous work in which the upper bounds of the error were developed in a FEM framework. The upper bound property requires the recovered solution to be equilibrated and continuous. The proposed technique consists of using a recovery technique, especially adapted to the XFEM framework that yields equilibrium at a local level (patch by patch). Then a postprocess based on the partition of unity concept is used to obtain continuity. The result is a very accurate but only nearly-statically admissible recovered stress field, with small equilibrium defaults introduced by the postprocess. Sharp upper bounds are obtained using a new methodology accounting for the equilibrium defaults, as demonstrated by the numerical tests., Postprint (published version)