Documents published in Scipedia

• Numerical modelling of tectonic plates subduction using X-FEM

  S. Zlotnik, P. Diez, M. Fernández, J. Vergés

  Comput. Methods Appl. Mech. Engrg., (2007). Vol. 196 (41-44), pp. 4283-4293

  
  

  Abstract

  The numerical modelling of plate subduction requires solving a coupled thermo-mechanical highly-nonlinear transient problem. The mechanical description of the phenomenon results in a multiphase quasi-static Stokes flow, where the inertia terms are neglected. The transient thermal problem is dominated by the advection term. Here, the representation and evolution of the different phases are described using level sets. The phase tracking is carried out transporting the level set along with the material, using a pure advective model. The gradient discontinuities induced by the viscosity jump across the interface are resolved numerically by enriching the solution using a partition of unity method in a eXtended Finite Element Method (X-FEM) context. These numerical tools are used to simulate plate subduction with different parameters and to derive useful correlations between relevant geophysical factors.

  Abstract
  The numerical modelling of plate subduction requires solving a coupled thermo-mechanical highly-nonlinear transient problem. The mechanical description of the phenomenon results [...]

  • 4
  •  read

• Equilibrated Patch Recovery error estimates: simple and accurate upper bounds of the error

  P. Diez, J. Ródenas, O. Zienkiewicz

  Int. J. Numer. Meth. Engng. (2007). Vol. 69 (10), pp. 2075-2098

  
  

  Abstract

  This paper introduces a new recovery‐type error estimator ensuring local equilibrium and yielding a guaranteed upper bound of the error. The upper bound property requires the recovered solution to be both statically equilibrated and continuous. The equilibrium is obtained locally (patch‐by‐patch) and the continuity is enforced by a postprocessing based on the partition of the unity concept. This postprocess is expected to preserve the features of the locally equilibrated stress field. Nevertheless, the postprocess phase modifies the equilibrium, which is no longer exactly fulfilled. A new methodology is introduced that yields upper bound estimates by taking into account this lack of equilibrium. This requires computing the ℒ︁2 norm of the error or relating it with the energy norm. The guaranteed upper bounds are obtained by using a pessimistic bound of the error ℒ︁2 norm, derived from an eigenvalue problem. Nevertheless, these bounds are not sharp. An additional strategy based on a more accurate assessment of the error ℒ︁2 norm is introduced, providing sharp estimates, which are practical upper bounds as it is demonstrated in the numerical tests.

  Abstract
  This paper introduces a new recovery‐type error estimator ensuring local equilibrium and yielding a guaranteed upper bound of the error. The upper bound property requires [...]

  • 3
  •  read

• Goal-oriented error estimation for transient parabolic problems

  P. Diez, G. Calderón

  Computational Mechanics (2007). Vol. 39 (5), pp. 631-646

  
  

  Abstract

  This work focuses on controlling the error and adapting the discretization in the context of parabolic problems. In order to obtain a sound mathematical framework, the time domain is discretized using a Discontinuous Galerkin (DG) approach. This allows to formulate the time stepping procedure in a variational format. The error is measured in the basis of an output of interest of the solution, defined by a linear functional. A dual problem, associated with this linear output is introduced. The dual problem has to be solved backward in time. An error representation is introduced, based on the weak residual of the primal error applied to the dual solution. Two different alternatives are studied to estimate the error in the dual solution: 1) recovery based error estimators and 2) implicit residual type estimators. Once the error assessment is performed implicitly in the dual problem, the obtained estimate is plugged into the primal residual to obtain the error in the quantity of interest. The implementation of the estimator is drastically simplified by using the weak version of the residual instead of the strong version used in previous works. Thus, the output error is assessed using a mixed technique, explicit for the primal problem and implicit for the dual. In the framework of adaptive computations of transient problems, this approach is very attractive because it allows using first the implicit scheme for the dual problem and then integrating the primal problem, estimating the error explicitly and eventually adapting the space-time grid. Thus, at every time step of the time marching scheme, the estimate of the dual error is injected into the primal residual (explicit estimate for the primal problem).

  Abstract
  This work focuses on controlling the error and adapting the discretization in the context of parabolic problems. In order to obtain a sound mathematical framework, the time [...]

  • 0
  •  read

• Remeshing criteria and proper error representations for goal oriented h-adaptivity

  P. Diez, G. Calderón

  Comput. Methods Appl. Mech. Engrg., (2007). Vol. 196 (4-6), pp. 719-733

  
  

  Abstract

  A key ingredient in h-adaptivity pertains to the transformation of output data from a given error estimator into input data, usually in the form of an element-size distribution, that needs to be supplied to a mesh generator. This paper analyzes the different possibilities of de¯ning remeshing criteria in the framework of goal oriented adaptivity. In standard energy norm driven adaptivity, the optimal mesh is clearly obtained if the local error distribution is uniform. The goal-oriented paradigm introduces new di±culties associated with the different possibilities for the spatial error representation and the signs of the local error contributions. A nodal error representation is introduced in order to improve the communication with the automatic mesh generation tool, precluding the transfer of information from elements to nodes. Numerical experiments demonstrate the ability of the introduced remeshing strategies to drive e±cient adaptive procedures and to control the error in quantities of interest. The results of the numerical tests fit the expected properties of the different remeshing strategies.

  Abstract
  A key ingredient in h-adaptivity pertains to the transformation of output data from a given error estimator into input data, usually in the form of an element-size distribution, [...]

  • 2
  •  read

• Subdomain-based flux-free a posteriori error estimators

  N. Parés, P. Diez, A. Huerta

  Comput. Methods Appl. Mech. Engrg., (2006). Vol. 195 (4-6), pp. 297-323

  
  

  Abstract

  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 .nite 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.

  Abstract
  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 [...]

  • 0
  •  read

• Accurate upper and lower error bounds by solving flux-free local problems in “stars”

  P. Diez, N. Parés, A. Huerta

  European Journal of Finite Elements (2004). Vol. 13 (5-7), pp. 497-507

  
  

  Abstract

  Two implicit residual type estimators yielding upper bounds of the error are presented which do not require flux equilibration. One of them is based on the ideas introduced in [MAC 00, CAR 99, MOR 03, PRU 02]. The new approach introduced here is based on using the estimated error function rather than the estimated error norms. Once the upper bounds are computed, also lower bounds for the error are obtained with little supplementary effort.

  Abstract
  Two implicit residual type estimators yielding upper bounds of the error are presented which do not require flux equilibration. One of them is based on the ideas introduced [...]

  • 3
  •  read

• Convergence of finite elements enriched with meshless methods

  S. Fernández-Méndez, P. Diez, A. Huerta

  Numerische Mathematik (2003). Vol. 96 (1), pp. 43-59

  
  

  Abstract

  A combined hierarchical approximation based on finite elements and mesh-less methods is proposed and studied. Finite Elements are enriched adding hierarchical shape functions based on a particle distribution. Convergence results are presented and proved.

  Abstract
  A combined hierarchical approximation based on finite elements and mesh-less methods is proposed and studied. Finite Elements are enriched adding hierarchical shape functions [...]

  • 2
  •  read

• A note on the convergence of the secant method for simple and multiple roots

  P. Diez

  Applied Numerical Letters (2003). Vol. 16 (8), pp. 1211-1215

  
  

  Abstract

  The secant method is one of the most popular methods for root finding. Standard text books in numerical analysis state that the secant method is super linear: the rate of convergence is set by the gold number. Nevertheless, this property holds only for simple roots. If the multiplicity of the root is larger than one, the convergence of the secant method becomes linear. This communication includes a detailed analysis of the secant method when it is used to approximate multiple roots. Thus, a proof of the linear convergence is shown. Moreover, the values of the corresponding asymptotic convergence factors are determined and are found to be also related with the golden ratio.

  Abstract
  The secant method is one of the most popular methods for root finding. Standard text books in numerical analysis state that the secant method is super linear: the rate of [...]

  • 8
  •  read

• Goal-oriented adaptivity for shell structures: Error assessment and remeshig criteria

  P. Diez, I. Morata, A. Huerta

  European Journal of Finite Elements (2003). Vol. 12 (6), pp. 691-715

  
  

  Abstract

  The reliable computation of shell structures requires a tool to assess and control the quality of the finite element solution. For practical purposes, the quality of the numerical solution must be measured using a quantity of engineering interest rather than in the standard energy norm. However, the assessment of the error in an output of interest is based on a standard energy norm error estimator. The standard error estimator has to be applied to both the original problem (primal) and a dual problem related with the selected engineering quantity. In shells with assumed-strain models, the combination of primal and dual error estimation is performed differently than in the continuum mechanics case. Moreover, a part from the goal-oriented error estimator, the adaptive process requires a remeshing criterion. This work introduces an specific remeshing criterion for goal-oriented adaptivity and its particularization to the context of shell elements.

  Abstract
  The reliable computation of shell structures requires a tool to assess and control the quality of the finite element solution. For practical purposes, the quality of the numerical [...]

  • 3
  •  read

• Recovering lower bounds of the error by postprocessing implicit residual a posteriori error estimates

  P. Diez, N. Parés, A. Huerta

  Int. J. Numer. Meth. Engng. (2003). Vol. 56 (10), pp. 1465-1488

  
  

  Abstract

  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.

  Abstract
  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 [...]

  • 4
  •  read