Documents published in Scipedia

• A Proper Generalized Decomposition (PGD) approach to crack propagation in brittle materials: with application to random field material properties

  H. Garikapati, S. Zlotnik, C. Verhoosel, E. Brummelen, P. Diez

  Computational Mechanics (2019). (Preprint) pp. 1-23

  
  

  Abstract

  Understanding the failure of brittle heterogeneous materials is essential in many applications. Heterogeneities in material properties are frequently modeled through random fields, which typically induces the need to solve finite element problems for a large number of realizations. In this context, we make use of reduced order modeling to solve these problems at an affordable computational cost. This paper proposes a reduced order modeling framework to predict crack propagation in brittle materials with random heterogeneities. The framework is based on a combination of the Proper Generalized Decomposition (PGD) method with Griffith’s global energy criterion. The PGD framework provides an explicit parametric solution for the physical response of the system. We illustrate that a non-intrusive sampling-based technique can be applied as a postprocessing operation on the explicit solution provided by PGD.We first validate the framework using a global energy approach on a deterministic two-dimensional linear elastic fracture mechanics benchmark. Subsequently, we apply the reduced order modeling approach to a stochastic fracture propagation problem.

  Abstract
  Understanding the failure of brittle heterogeneous materials is essential in many applications. Heterogeneities in material properties are frequently modeled through random [...]

  • 12
  •  read

• A new 3D equilibrated residual method improving accuracy and efficiency of flux‐free error estimates

  N. Parés, P. Diez

  Comput. Methods Appl. Mech. Engrg., (2019). (Preprint) Vol. 120, pp. 391-432

  
  

  Abstract

  The paper presents a novel strategy providing fully computable upper bounds for the energy norm of the error in the context of three‐dimensional linear finite element approximations of the reaction‐diffusion equation. The upper bounds are guaranteed regardless the size of the finite element mesh and the given data, and all the constants involved are fully computable. The upper bound property holds if the shape of the domain is polyhedral and the Dirichlet boundary conditions are piecewise‐linear. The new approach is an extension of the flux‐free methodology introduced by Parés and Díez in the paper “A new equilibrated residual method improving accuracy and efficiency of flux‐free error estimates”, which introduces a guaranteed, low‐cost, and efficient flux‐free method substantially reducing the computational cost of obtaining guaranteed bounds using flux‐free methods while retaining the good quality of the bounds. Besides extending the 2D methodology, specific new modifications are introduced to further reduce the computational cost in the three‐dimensional setting. The presented methodology also provides a new strategy to obtain equilibrated boundary tractions, which improves the quality of standard techniques while having a similar computational cost.

  Abstract
  The paper presents a novel strategy providing fully computable upper bounds for the energy norm of the error in the context of three‐dimensional linear finite element approximations [...]

  • 1
  •  read

• Sampling-based stochastic analysis of the PKN model for hydraulic fracturing

  H. Garikapati, C. Verhoosel, E. Brummelen, S. Zlotnik, P. Diez

  Computational Geosciences (2019). (Preprint) Vol. 23 (1), pp. 81-105

  
  

  Abstract

  Hydraulic fracturing processes are surrounded by uncertainty, as available data is typically scant. In this work, we present a sampling-based stochastic analysis of the hydraulic fracturing process by considering various system parameters to be random. Our analysis is based on the Perkins-Kern-Nordgren (PKN) model for hydraulic fracturing. This baseline model enables computation of high fidelity solutions, which avoids pollution of our stochastic results by inaccuracies in the deterministic solution procedure. In order to obtain the desired degree of accuracy of the computed solution, we supplement the employed time-dependent moving-mesh finite element method with two new enhancements: (i) global conservation of volume is enforced through a Lagrange multiplier; (ii) the weakly singular behavior of the solution at the fracture tip is resolved by supplementing the solution space with a tip enrichment function. This tip enrichment function enables the computation of the tip speed directly from its associated solution coefficient. A novel incremental-iterative solution procedure based on a backward-Euler time-integrator with sub-iterations is employed to solve the PKN model. Direct Monte-Carlo sampling is performed based on random variable and random field input parameters. The presented stochastic results quantify the dependence of the fracture evolution process—in particular the fracture length and fracture opening—on variations in the elastic properties and leak-off coefficient of the formation, and the height of the fracture.

  Abstract
  Hydraulic fracturing processes are surrounded by uncertainty, as available data is typically scant. In this work, we present a sampling-based stochastic analysis of the hydraulic [...]

  • 4
  •  read

• Algebraic and parametric solvers for the power flow problem: towards real-time and accuracy-guaranteed simulation of electric systems

  R. García-Blanco, P. Diez, D. Borzacchiello, F. Chinesta

  Archives of Comp. Meths. Engng. (2018). (Preprint) Vol. 25 (4), pp. 1003-1026

  
  

  Abstract

  The power flow model performs the analysis of electric distribution and transmission systems. With this statement at hand, in this work we present a summary of those solvers for the power flow equations, in both algebraic and parametric version. The application of the Alternating Search Direction method to the power flow problem is also detailed. This results in a family of iterative solvers that combined with Proper Generalized Decomposition technique allows to solve the parametric version of the equations. Once the solution is computed using this strategy, analyzing the network state or solving optimization problems, with inclusion of generation in real-time, becomes a straightforward procedure since the parametric solution is available. Complementing this approach, an error strategy is implemented at each step of the iterative solver. Thus, error indicators are used as an stopping criteria controlling the accuracy of the approximation during the construction process. The application of these methods to the model IEEE 57-bus network is taken as a numerical illustration.

  Abstract
  The power flow model performs the analysis of electric distribution and transmission systems. With this statement at hand, in this work we present a summary of those solvers [...]

  • 3
  •  read

• Algebraic PGD for tensor separation and compression: An algorithmic approach

  P. Diez, S. Zlotnik, A. García-González, A. Huerta

  Comptes Rendus Mécanique (2018). (Preprint) Vol. 346 (7), pp. 501-514

  
  

  Abstract

  Proper Generalized Decomposition (PGD) is devised as a computational method to solve high-dimensional boundary value problems (where many dimensions are associated with the space of parameters defining the problem). The PGD philosophy consists in providing a separated representation of the multidimensional solution using a greedy approach combined with an alternated directions scheme to obtain the successive rank-one terms. This paper presents an algorithmic approach to high-dimensional tensor separation based on solving the Least Squares approximation in a separable format of multidimensional tensor using PGD. This strategy is usually embedded in a standard PGD code in order to compress the solution (reduce the number of terms and optimize the available storage capacity), but it stands also as an alternative and highly competitive method for tensor separation.

  Abstract
  Proper Generalized Decomposition (PGD) is devised as a computational method to solve high-dimensional boundary value problems (where many dimensions are associated [...]

  • 4
  •  read

• Numerical estimation of the bearing capacity of resistance spot welds in martensitic boron steels using a J-integral fracture criterion

  D. Dorribo, L. Greve, P. Diez, I. Arias, X. Larráyoz-Izcara

  Theoretical and Applied Fracture Mechanics (2018). (Preprint) Vol. 96, pp. 497-508

  
  

  Abstract

  Predicting the bearing capacity of resistance spot welds (RSW) during vehicle crash tests has become a crucial task for the automotive industry, since the recent introduction of advanced high strength steels (AHSS) such as martensitic boron steels (e.g. 22MnB5). The spot weld joints of these steels exhibit relatively low bearing strengths, compared to those of more ductile high strength steels. Currently, the bearing capacity of spot weld joints is characterized through extensive experimental campaigns. In this article, a model for quantification of the bearing capacity of RSW using a finite-element J-integral fracture criterion is presented. The model takes into account geometric and mechanical features of the spot weld, namely the weld diameter and the mechanical properties distribution resulting from the welding process. An experimental loading test campaign is carried out for calibration and validation purposes, considering multiple sheet thickness combinations, loading angles and weld sizes. Experimental observations of the failed spot welds and preliminary simulations show that failure is caused mostly by stress concentration around the sharp weld notch. Consequently, the J-integral obtained from detailed finite element simulations is used to asses the stress/strain concentration along the first crack advance direction predicted by the acoustic tensor. The computed J-integral values are compared to the material toughness to obtain the joint’s maximum force. The resulting simulated and experimental bearing capacities show a good agreement for all tested configurations.

  Abstract
  Predicting the bearing capacity of resistance spot welds (RSW) during vehicle crash tests has become a crucial task for the automotive industry, since the recent introduction [...]

  • 2
  •  read

• A non-intrusive proper generalized decomposition scheme with application in biomechanics

  X. Zou, M. Conti, P. Diez, F. Auricchio

  Int. J. Numer. Meth. Engng. (2018). (Preprint) Vol. 113 (2), pp. 230-251

  
  

  Abstract

  Proper generalized decomposition (PGD) is often used for multi-query and fast-response simulations. It is a powerful tool alleviating the curse of dimensionality affecting multi-parametric partial differential equations. Most implementations of PGD are intrusive extensions based on in-house developed finite element (FE) solvers. In this work, we propose a non-intrusive PGD scheme using off-the-shelf FE codes (such as certified commercial software) as an external solver. The scheme is implemented and monitored by in-house flow-control codes. A typical implementation is provided with downloadable codes. Moreover, a novel parametric separation strategy for the PGD resolution is presented. The parametric space is split into two- or three-dimensional subspaces, to allow PGD technique solving problems with constrained parametric spaces, achieving higher convergence ratio. Numerical examples are provided. In particular, a practical example in biomechanics is included, with potential application to patient-specific simulation.

  Abstract
  Proper generalized decomposition (PGD) is often used for multi-query and fast-response simulations. It is a powerful tool alleviating the curse of dimensionality affecting [...]

  • 19
  •  read

• Explicit parametric solutions of lattice structures with proper generalized decomposition (PGD)

  A. Sibileau, A. García-González, F. Auricchio, S. Morganti, P. Diez

  Computational Mechanics (2018). (Preprint) Vol. 62 (4), pp. 871-891

  
  

  Abstract

  Architectured materials (or metamaterials) are constituted by a unit-cell with a complex structural design repeated periodically forming a bulk material with emergent mechanical properties. One may obtain specific macro-scale (or bulk) properties in the resulting architectured material by properly designing the unit-cell. Typically, this is stated as an optimal design problem in which the parameters describing the shape and mechanical properties of the unit-cell are selected in order to produce the desired bulk characteristics. This is especially pertinent due to the ease manufacturing of these complex structures with 3D printers. The proper generalized decomposition provides explicit parametic solutions of parametric PDEs. Here, the same ideas are used to obtain parametric solutions of the algebraic equations arising from lattice structural models. Once the explicit parametric solution is available, the optimal design problem is a simple post-process. The same strategy is applied in the numerical illustrations, first to a unit-cell (and then homogenized with periodicity conditions), and in a second phase to the complete structure of a lattice material specimen.

  Abstract
  Architectured materials (or metamaterials) are constituted by a unit-cell with a complex structural design repeated periodically forming a bulk material with emergent mechanical [...]

  • 7
  •  read

• Numerical Modelling of Multi-Phase Multi-Component. Reactive Transport in the Earth’s interior

  B. Oliveira, J. Afonso, S. Zlotnik, P. Diez

  Geophysical Journal International (2018). (Preprint) Vol. 212 (1), pp. 345-388

  
  

  Abstract

  We present a conceptual and numerical approach to model processes in the Earth’s interior that involve multiple phases that simultaneously interact thermally, mechanically and chemically. The approach is truly multiphase in the sense that each dynamic phase is explicitly modelled with an individual set of mass, momentum, energy and chemical mass balance equations coupled via interfacial interaction terms. It is also truly multicomponent in the sense that the compositions of the system and its constituent phases are expressed by a full set of fundamental chemical components (e.g. SiO2, Al2O3, MgO, etc.) rather than proxies. These chemical components evolve, react with and partition into different phases according to an internally consistent thermodynamic model. We combine concepts from Ensemble Averaging and Classical Irreversible Thermodynamics to obtain sets of macroscopic balance equations that describe the evolution of systems governed by multiphase multicomponent reactive transport (MPMCRT). Equilibrium mineral assemblages, their compositions and physical properties, and closure relations for the balance equations are obtained via a ‘dynamic’ Gibbs free-energy minimization procedure (i.e. minimizations are performed on-the-fly as needed by the simulation). Surface tension and surface energy contributions to the dynamics and energetics of the system are taken into account. We show how complex rheologies, that is, visco-elasto-plastic, and/or different interfacial models can be incorporated into our MPMCRT ensemble-averaged formulation. The resulting model provides a reliable platform to study the dynamics and nonlinear feedbacks of MPMCRT systems of different nature and scales, as well as to make realistic comparisons with both geophysical and geochemical data sets. Several numerical examples are presented to illustrate the benefits and limitations of the model.

  Abstract
  We present a conceptual and numerical approach to model processes in the Earth’s interior that involve multiple phases that simultaneously interact thermally, mechanically [...]

  • 8
  •  read

• Generalized parametric solutions in Stokes flow

  P. Diez, S. Zlotnik, A. Huerta

  Comput. Methods Appl. Mech. Engrg. (2017). (Preprint) Vol. 326, pp. 223-240

  
  

  Abstract

  Design optimization and uncertainty quantification, among other applications of industrial interest, require fast or multiple queries of some parametric model. The Proper Generalized Decomposition (PGD) provides a separable solution, a computational vademecum explicitly dependent on the parameters, efficiently computed with a greedy algorithm combined with an alternated directions scheme and compactly stored. This strategy has been successfully employed in many problems in computational mechanics. The application to problems with saddle point structure raises some difficulties requiring further attention. This article proposes a PGD formulation of the Stokes problem. Various possibilities of the separated forms of the PGD solutions are discussed and analyzed, selecting the more viable option. The efficacy of the proposed methodology is demonstrated in numerical examples for both Stokes and Brinkman models.

  Abstract
  Design optimization and uncertainty quantification, among other applications of industrial interest, require fast or multiple queries of some parametric model. The Proper [...]

  • 5
  •  read