Monitoring a PGD solver for parametric power flow problems with goal-oriented error assessment

This is the peer reviewed version of the following article: [García-Blanco, R., Borzacchiello, D., Chinesta, F., and Diez, P. (2017) Monitoring a PGD solver for parametric power flow problems with goal-oriented error assessment. Int. J. Numer. Meth. Engng, 111: 529–552. doi: 10.1002/nme.5470], which has been published in final form at http://onlinelibrary.wiley.com/doi/10.1002/nme.5470/full. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.The parametric analysis of electric grids requires carrying out a large number of Power Flow computations. The different parameters describe loading conditions and grid properties. In this framework, the Proper Generalized Decomposition (PGD) provides a numerical solution explicitly accounting for the parametric dependence. Once the PGD solution is available, exploring the multidimensional parametric space is computationally inexpensive. The aim of this paper is to provide tools to monitor the error associated with this significant computational gain and to guarantee the quality of the PGD solution. In this case, the PGD algorithm consists in three nested loops that correspond to 1) iterating algebraic solver, 2) number of terms in the separable greedy expansion and 3) the alternated directions for each term. In the proposed approach, the three loops are controlled by stopping criteria based on residual goal-oriented error estimates. This allows one for using only the computational resources necessary to achieve the accuracy prescribed by the end- user. The paper discusses how to compute the goal-oriented error estimates. This requires linearizing the error equation and the Quantity of Interest to derive an efficient error representation based on an adjoint problem. The efficiency of the proposed approach is demonstrated on benchmark problems.Peer ReviewedPostprint (author's final draft

Proper generalized decomposition solution of the parameterized Helmholtz problem: application to inverse geophysical problems

The identification of the geological structure from seismic data is formulated as an inverse problem. The properties and the shape of the rock formations in the subsoil are described by material and geometric parameters, which are taken as input data for a predictive model. Here, the model is based on the Helmholtz equation, describing the acoustic response of the system for a given wave length. Thus, the inverse problem consists in identifying the values of these parameters such that the output of the model agrees the best with observations. This optimization algorithm requires multiple queries to the model with different values of the parameters. Reduced order models are especially well suited to significantly reduce the computational overhead of the multiple evaluations of the model. In particular, the proper generalized decomposition produces a solution explicitly stating the parametric dependence, where the parameters play the same role as the physical coordinates. A proper generalized decomposition solver is devised to inexpensively explore the parametric space along the iterative process. This exploration of the parametric space is in fact seen as a post‐process of the generalized solution. The approach adopted demonstrates its viability when tested in two illustrative examples

Generalized parametric solutions in Stokes flow

The final publication is available at Springer via https://doi.org/10.1016/j.cma.2017.07.016Design 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.Peer ReviewedPostprint (author's final draft

Static and dynamic global stiffness analysis for automotive pre-design

In order to be worldwide competitive, the automotive industry is constantly challenged to produce higher quality vehicles in the shortest time possible and with the minimum costs of production. Most of the problems with new products derive from poor quality design processes, which often leads to undesired issues in a stage where changes are extremely expensive. During the preliminary design phase, designers have to deal with complex parametric problems where material and geometric characteristics of the car components are unknown. Any change in these parameters might significantly affect the global behaviour of the car. A target which is very sensitive to small variations of the parameters is the noise and vibration response of the vehicle (NVH study), which strictly depends on its global static and dynamic stiffness. In order to find the optimal solution, a lot of configurations exploring all the possible parametric combinations need to be tested. The current state of the art in the automotive design context is still based on standard numerical simulations, which are computationally very expensive when applied to this kind of multidimensional problems. As a consequence, a limited number of configurations is usually analysed, leading to suboptimal products. An alternative is represented by reduced order method (ROM) techniques, which are based on the idea that the essential behaviour of complex systems can be accurately described by simplified low-order models.This thesis proposes a novel extension of the proper generalized decomposi-tion (PGD) method to optimize the design process of a car structure with respect to its global static and dynamic stiffness properties. In particular, the PGD method is coupled with the inertia relief (IR) technique and the inverse power method (IPM) to solve, respectively, the parametric static and dynamic stiffness analysis of an unconstrained car structure and extract its noise and vibrations properties. A main advantage is that, unlike many other ROM methods, the proposed approach does not require any pre-processing phase to collect prior knowledge of the solution. Moreover, the PGD solution is computed with only one offline computation and presents an explicit dependency on the introduced design variables. This allows to compute the solutions at a negligible computational cost and therefore opens the door to fast optimisation studies and real-time visualisations of the results in a pre-defined range of parameters. A novel algebraic approach is also proposed which allows to involve both material and com-plex geometric parameters, such that shape optimisation studies can be performed. In addition, the method is developed in a nonintrusive format, such that an interaction with commercial software is possible, which makes it particularly interesting for industrial applications. Finally, in order to support the designers in the decision-making process, a graphical interface app is developed which allows to visualise in real-time how changes in the design variables affect pre-defined quantities of interest

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

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

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

This is the peer reviewed version of the following article: Zou, X., Conti, M., Diez, P., Auricchio, F. A non-intrusive proper generalized decomposition scheme with application in biomechanics. "International journal for numerical methods in engineering", 7 Setembre 2017., which has been published in final form at http://dx.doi.org/10.1002/nme.5610. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Self-Archiving.Proper generalized decomposition (PGD) is often used for multiquery and fast-response simulations. It is a powerful tool alleviating the curse of dimensionality affecting multiparametric partial differential equations. Most implementations of PGD are intrusive extensions based on in-house developed FE solvers. In this work, we propose a nonintrusive 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.Peer ReviewedPostprint (author's final draft

Virtual, Digital and Hybrid Twins: A New Paradigm in Data-Based Engineering and Engineered Data

Engineering is evolving in the same way than society is doing. Nowadays, data is acquiring a prominence never imagined. In the past, in the domain of materials, processes and structures, testing machines allowed extract data that served in turn to calibrate state-of-the-art models. Some calibration procedures were even integrated within these testing machines. Thus, once the model had been calibrated, computer simulation takes place. However, data can offer much more than a simple state-of-the-art model calibration, and not only from its simple statistical analysis, but from the modeling and simulation viewpoints. This gives rise to the the family of so-called twins: the virtual, the digital and the hybrid twins. Moreover, as discussed in the present paper, not only data serve to enrich physically-based models. These could allow us to perform a tremendous leap forward, by replacing big-data-based habits by the incipient smart-data paradigm