EDP Sciences logo
Submit your paper
Search
Menu
All issues Volume 23 (2022) Mechanics & Industry, 23 (2022) 3 Full HTML
Open Access
Issue
Mechanics & Industry
Volume 23, 2022
Article Number 3
Number of page(s) 20
DOI https://doi.org/10.1051/meca/2022001
Published online 07 February 2022

© T. Daniel et al., Published by EDP Sciences 2022

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

Numerical simulation is vastly used in the industry when designing a new mechanical part. Uncertainty quantification is used to compute the influence of a poorly known or uncontrolled parameter, like dispersion within manufacturing tolerances. Procedures that rely on the Monte Carlo method require solving the corresponding problem for many values of the parameter. With modeling and simulation progresses, meshes are getting larger and models more complex, leading to increased duration for the corresponding computations. For these reasons, many methods have been proposed to replace these reference models with fast approximations.

In this work, we consider an industrial model of an elastoviscoplastic high-pressure turbine blade subjected to thermal, centrifugal and pressure loadings, for the quantification of the uncertainty of dual quantities (such as the accumulated plastic strain and the stress tensor), generated by the uncertainty of the temperature loading. Computing the fatigue lifetime of one such blade requires simulating its behavior until the stabilization of the mechanical response, which last several weeks using Abaqus [1] because of the size of the mesh, the complexity of the constitutive equations, and the number of loading cycles in the transient regime. With such a computation time, uncertainty quantification with the Monte Carlo method is unaffordable. In addition, such simulations are too time-consuming to be integrated in design iterations, which limits them to the final verification steps, while the design process still relies on simplified models. Accelerating these complex simulations is a key challenge while maintaining a satisfying accuracy, as it would provide useful numerical tools to improve design processes and quantify the effect of the uncertainties on the environment of the system.

Simulations are accelerated using a dictionary of reduced order models, with a classifier able to select which local reduced order model to be used for a new temperature loading. The framework is called ROM-net [2]. A dataset of 200 solutions is computed in a Finite Element approximation space of dimension in the order of the million, for various instances of the temperature field loading, in parallel in 7 days and 9 hours on 48 cores. These solutions are computed over 11 time steps in the first cycle, using a scalable Adaptive MultiPreconditioned FETI (AMPFETI) solver [3] in Z-set finite-element software [4]. The dataset is partitioned into two clusters using a k-medoids algorithm with a Reduced Order Model (ROM)-oriented dissimilarity measure in 5 minutes; the corresponding local ROMs, using Proper Orthogonal Decomposition (POD) [5,6] and Empirical Cubature Method (ECM) [7], are trained in 2 hours and 30 minutes. An automatic reduced model recommendation procedure, allowing to decide which local ROM to be used for a new temperature loading, is trained in the form of a logistic regression classifier in 16 minutes. A contribution of this work is the use of a meta-model to reconstruct the dual quantities of interest over the complete mesh from their values on the reduced integration points, in the form of a multi-task Lasso, which takes 1 hour to train for 14 dual fields. The uncertainties on dual quantities of interest, such as the accumulated plastic strain and the stress tensor, are quantified by using our trained ROM-net on 1008 Monte Carlo draws of the temperature loading field in 2 hours and 48 minutes, which corresponds to a speedup greater than 600 with respect to our highly optimized domain decomposition AMPFETI solver. Expected values for the Von Mises stress and the accumulated plastic strain have 0.99-confidence intervals width of respectively 1.66% and 2.84% (relative to the corresponding prediction for the expected value). As a verification stage, 20 reference solutions are computed on new temperature loadings, and dual quantities of interest are predicted with relative accuracy in the order of 1–2%, while the location of the maximum value is perfectly predicted.

In Section 2, a summary of the results and the methodology constructed in our team’s previous work on dictionary-based ROM-net is provided. Details are given on the reduced order modeling strategy and the concept of ROM-net is recalled. References on alternative methods from the literature are also given. Section 3 contains a description of the industrial dataset, the hypotheses of the model, and the objective of the present study. The proposed workflow for uncertainty quantification is then applied on this industrial configuration in Section 4. Finally, conclusions are drawn in Section 5.

2 Description of the methodology

This sectionprovides elements from the literature and from our previous works to deal with uncertainty quantification in an industrial context using physical reduced order modeling.

2.1 Reduced order modeling of nonlinear parametrized partial differential equations

Consider a nonlinear physical problem described by the following parametrized differential equation: (1)

where u is the primal variable belonging to a Hilbert space , x denotes the parameters of the problem, and is an operator involving a differential operator and operators for initial conditions and/or boundary conditions. Here, the time is considered as a parameter and is included in the definition of x. The solution manifold is defined by . Classically, the problem is written on a geometrical support discretized by a mesh, and the solution is sought in a finite dimensional space, e.g. the finite element space : this defines the High Fidelity Model (HFM). Here, in the target application, is in the order of the million.

Model order reduction [9,8] is a discipline in numerical analysis consisting in replacing a computationally expensive high fidelity model by a fast ROM to calculate approximate solutions of the considered physical problem (1). A ROM can be either a data-driven metamodel (or surrogate model) calibrated with a regression algorithm, or a physics-based model obtained by numerical methods such as the Proper Generalized Decomposition [10,11], the Reduced Basis Method (RBM) [12,13], and the POD-Galerkin method [5,6], among others. It is generally used for parametrized equations whose solution must be known for different points in the parameter space. As in machine learning, a model order reduction procedure starts by a training phase (or offline stage) where the ROM is built from some training data. The ROM is then used on test data in an exploitation phase (or online stage). In the training phase, high-fidelity solutions, called snapshots, are computed with the HFM for different points of the parameter space to get a sampled representation of the solution manifold. The model order reduction algorithm analyzes these snapshots to learn how the solution is affected by parameter variations. Given the cost of computing snapshots in the training phase, a ROM is profitable only if it is extensively used in the exploitation phase.

RBM and POD-Galerkin are also called projection-based since they consist in applying the Galerkin method on a reduced order basis, which can be costly for certain parameter dependencies and nonlinear problems. In such cases, a second reduction stage is necessary. Hyper-reduction was initially the name of a method proposed in [14] in 2005, but this term has been extended to refer to all the methods proposing a second reduction stage. Hyper-reduction methods include the Empirical Interpolation Method (EIM, [15]), the Missing Point Estimation (MPE, [16]), the A Priori Hyper-Reduction (APHR, [14]), the Best Point Interpolation Method (BPIM, [17]), the Discrete Empirical Interpolation Method (DEIM, [18]), the Gauss-Newton with Approximated Tensors (GNAT, [19]), the Energy-Conserving Sampling and Weighting (ECSW, [20]), the Empirical Cubature Method (ECM, [7]), and the Linear Program Empirical Quadrature Procedure LPEQP, [21]). Hyper-reduction techniques implicitly assume that employed discretization technique and the equations modeling constitutive laws evaluated at a given point ξ ∈ Ω only involve variables evaluated at ξ.

We use the reduced order modeling framework developed in our previous work [22] for elastoviscoplastic structural mechanics, whose training phase is decomposed in three steps:

  • Data generation: snapshots u(xn), 1 ≤ n ≤ Ns, Ns being the number of available snapshots, are computed with the high-fidelity model and provide information about how the physical system reacts to changes of the parameter x. In our case, the finite-element solver Z-set [4] is used.

  • Data compression: a Reduced-Order Basis (ROB) is constructed by looking for a hidden low-rank structure in the snapshots. We apply the snapshot-POD, which consists in (i) computing the snapshot correlation matrix Cn,m = ∫Ωu(xn) ⋅ u(xm), 1 ≤ n, m ≤ Ns, (ii) retaining the eigenvalue/eigenvector pairs of C associated to the highest eigenvalues: (ξi, λi), 1 ≤ i ≤ N, and (iii) recombining them with the snapshots to create the ROB: , 1 ≤ i ≤ N. When the quantities of interest are dual variables, which is the case in the industrial application considered here, ROBs are also constructed for the corresponding fields.

  • Operator compression: this step contains the operations that guarantee the efficiency of the reduced order model in the exploitation phase. For nonlinear problems, applying the Galerkin method in the online phase on the ROB requires to compute integrals over the mesh. In our case, we use the ECM [7] to replace these costly integrals by reduced quadrature schemes trained over the snapshots at our disposal, hence tailored for our considered problem.

Then, the online stage consists in assembling a variational formulation of equation (1) on the ROB, using a reduced quadrature scheme, and solving it using a Newton algorithm. Since the constitutive laws are only computed at the locations of the integration points of the reduced quadrature scheme, the dual quantities have to be reconstructed over the complete mesh. This can be done using the Gappy-POD and the ROBs of the corresponding dual fields. See [22] for more details on the reduced order modeling framework we used. In Section 4.4, we propose to replace this last reconstruction stage by a meta-model.

2.2 Dictionary of reduced order models

This section introduces the concept of ROM-net as a dictionary of ROMs, and provides elements on its implementation.

2.2.1 Nonreducible problems

A ROM is an approximation of our considered HFM. We define the speedup of the ROM as the ratio of the computation time of the HFM to the computation time of the ROM. Consider a set of snapshots generated using the high-fidelity model over a sampling of the parameter domain. The parametrized problem is said nonreducible when applying a linear data compression over this set of snapshots leads to a ROB containing too many vectors for the online problem to feature an interesting speedup. Formally, this happens when the Kolmogorov N-width decreases too slowly with respect to N, where N is the cardinality of the ROB, (2)

with the Grassmannian being the set of all N-dimensional subspaces of and the subspace spanned by the considered ROB. Qualitatively, the solution manifold covers too many independent directions to be embedded in a low-dimensional subspace. To address this issue, several techniques have been developed:

  • Problem-specific methods tackle the difficulties of some specific physics problems that are known to be nonreducible, such as advection-dominated problems which have been largely investigated, for instance in [23–25].

  • Online-adaptive model reduction methods update the ROM in the exploitation phase by collecting new information online asexplained in [26], in order to limit extrapolation errors when solving the parametrized governing equations in a region of the parameter space that was not explored in the training phase. The ROM can be updated for example by querying the high-fidelity model when necessary for basis enrichment [14,27–30].

  • ROM interpolation methods [31–42] use interpolation techniques on Grassmann manifolds or matrix manifolds to adapt the ROM to the parameters considered in the exploitation phase by interpolating between two precomputed ROMs.

  • Dictionaries of basis vector candidates enable building a parameter-adapted ROM in the exploitation phase by selecting a few basis vectors. This technique is presented in [43,44] for the Reduced Basis method.

  • Dictionaries of ROMs rely on the construction of several local ROMs adapted to different regions of the solution manifold. These local ROMs can be obtained by partitioning the time interval [45,46], the parameter space [28,46–51], or the solution space [2,51–57].

  • Nonlinear manifold ROM methods [58–60] learn a nonlinear embedding and project the governing equations onto the corresponding approximation manifold, by means of a nonlinear function mapping a low-dimensional latent space to the solution space.

Our framework focuses on dictionaries of ROMs, where the solution manifold is partitioned to get a collection of subsets that can be covered by a dictionary of low-dimensional subspaces, enabling the use of local linear ROMs. If is a partition of , then: (3)

2.2.2 Dictionary-based ROM-nets

We introduced the concepts of ROM-net and dictionary-based ROM-net in [2], where rigorous definitions can be found. Suppose we dispose of an already computed dictionary of ROMs for the parametrized problem (1), where each element of the dictionary is a ROM that can approximate the problem on a subset of the solution manifold . A dictionary-based ROM-net is a machine learning algorithm trained to assign the parameter to the ROM of the dictionary leading to the most accurate reduced prediction. This assignment, called model recommendation in [55], is a classification task, see Figure 1.

The dictionary of ROMs is constructed in a clustering stage, during which snapshots are regrouped depending on their respective proximity on , in the sense of a particular dissimilarity measure we introduced in [61]. The dissimilarity between two parameter values , denoted by δ(x, x′), involves the sine of the principal angles between subspaces associated to the solutions of the HFM , see [61, Definition 4.10]. Applying a k-medoids clustering algorithm on the solution manifold using the dissimilarity δ leads to an optimal partitioning for a dictionary of local ROMs, in a sense introduced in [61, Property 4.13]. We refer to the remaining of [61] for the description of a practical efficiency criterion of the dictionary-based ROM-net, which enables to decide, before the computationally costly steps of the workflow, if a dictionary of ROMs is preferable to one global ROM, and how to calibrate the various hyperparameters of the ROM-net.

Remark 1.

Importance of the classification. One could argue that the classification step can be replaced by choosing the cluster k for which the dissimilarity measure between the parameter x and the cluster medoid is the smallest. However, we recall that the computation of the dissimilarity measure requires solving the HFM at the parameter value x, which wouldrender the complete model reduction framework useless. Hence, the classification step enables to bypass this HFM solve and directly recommend the appropriate local ROM.

As briefly mentioned in Section 2.2.1, local ROMs can be constructed by partitioning the parameter space [46,47], in which case the classification step is not required: the cluster affectation is made by computing distances directly in the parameter space. In other cases, partitioning in the solution space can be considered without requiring a classification step [52,57]. Consider a time-dependent problems where the initial condition is not a parameter of the problem, and suppose an efficient computation of the clustering distance in the solution space based on the reduced solution at the previous time-step. Then, local basis affectation and switching is possible without requiring classification. Efficient computation of inner-product-based dissimilarities (like L2 and sine) usually require an affine decomposition of the operators with respect to the parameters, which is not the case in our problem. Besides, base switching with complex constitutive laws would require reconstructing the complete material state in different reduced integration points sets, which could become inefficient.

The training of the classifier can be difficult when working with physical fields: simulations are costly, data are in high dimension and classical data augmentation techniques for images cannot be applied. Hence, we can consider replacing the HFM by an intermediate-fidelity solver for generating the data needed for the training of the classifier, by considering coarser meshes and fewer time steps. We point out that the HFM should be used at the end for generating the data required in the training of the local ROMs. We propose in [62] improvements for the training of the classifier in our context by developing a fast variant of the mRMR [63] feature selection algorithm, and new class-conserving transformations of our data, acting like a data augmentation procedure. In this work, there is no intermediate-fidelity solver: we use the same model for generating the data used in all the training tasks: clustering, classification, and the offline stage of the local reduced-order models.

thumbnail Fig. 1

Exploitation phase of a dictionary-based ROM-net. K local ROMs, combined with a classifier for automatic ROM recommendation, are used to predict the quantity of interest Z(x).

2.3 Uncertainty quantification

The parameter is modeled by a random variable. The uncertainty quantification consists in a Monte Carlo procedure where values of the parameter are drawn from the distribution of this random variable, and using the trained dictionary-based ROM-net to select a local ROM and predicting the corresponding quantity of interest. Statistical estimators of quantities depending on the solution of our physical problem can then be efficiently computed. Our objective is to apply the presented methodology to a real industrial case: quantifying the uncertainties on dual quantities of interest generated by the uncertainty of the temperature loading, in a high-pressure (HP) turbine blade elastoviscoplastic cyclic mechanical computation.

3 Industrial context

This section presents the industrial test case of interest. It consists in predicting the mechanical behavior of a HP turbine blade in an aircraft engine with uncertainties on the thermal loading. The industrial context and the models for the mechanical behavior and the thermal loading are presented, with a particular emphasis on the assumptions that have been made.

For confidentiality reasons, mesh sizes and numerical values corresponding to the industrial dataset are not given, and the provided figures and plots do not contain any color map or physical numerical value. Reproducible data are available on request for the numerical example proposed in [2]. The accuracy of the predictions made by our methodology are given in the form of relative errors.

3.1 Thermomechanical fatigue of high-pressure turbine blades

High-pressure turbine blades are critical parts in an aircraft engine. Located downstream of the combustion chamber, they are subjected toextreme thermomechanical loadings resulting from the combination of centrifugal forces, pressure loads, and hot turbulent fluid flows whose temperatures are higher than the material’s melting point. The repeating thermomechanical loading over time progressively damages the blades and leads to crack initiation under thermomechanical fatigue. Predicting the fatigue lifetime is crucial not only for safety reasons, but also for ecological issues, since reducing fuel consumption and improving the engine’s efficiency requires increasing the temperature of the gases leaving the combustion chamber.

High-pressure turbine blades are made of monocrystalline nickel-based superalloys that have good mechanical properties at high temperatures. To reduce the temperature inside this material, the blades contain cooling channels where flows relatively fresh air coming from the compressor. In addition, the blade’s outer surface is protected by a thin thermal barrier coating. In spite of these advanced cooling technologies, the rotor blades undergo centrifugal forces at high temperatures, causing inelastic strains. Under this cyclic thermomechanical loading repeated over the flights, the structure has a viscoplastic behavior and reaches a viscoplastic stabilized response, where the dissipated energy per cycle still has a nonzero value. This is called plastic shakedown, and leads tolow-cycle fatigue. At cruise flight, the persistent centrifugal force applied at high temperature induces progressive (or time-dependent) inelastic deformations: this phenomenon is called creep. In addition, the difference between gas pressures on the extrados and the intrados of the blade generates bending effects. Environmental factors may also locally modify the chemical composition of the material, leading to its oxidation. As oxidized parts are more brittle, they facilitate crack initiation and growth. Thermal fatigue resulting from temperature gradients is another life-limiting factor. Temperature gradients make colder parts of the structure prevent the thermal expansion of hotter parts, creating thermal stresses. Due to their higher temperatures, the hot parts are more viscous and have a lower yield stress, which make them prone to develop inelastic strains in compression. When the temperature cools down after landing, tensile residual stresses appear in parts which were compressed at high temperatures and favor crack nucleation. Given the complex temperature field resulting from the internal cooling channels and the turbulent gas flow, thermal fatigue has a strong influence on the turbine blade’s lifetime. In particular, during transient regimes such as take-off, an important temperature gradient appears between the leading edge and the trailing edge of the blade, since the latter has a low thermal inertia due to its small thickness and thus warms up faster.

In short, the behavior of a high pressure turbine blade results from a complex interaction between low-cycle fatigue, thermal fatigue, creep, and oxidation. Due to the cost and the complexity of experiments on parts of an aircraft engine, numerical simulations play a major role in the design of high-pressure turbine blades and their fatigue lifetime assessment. All this knowledge have been learned by scientist and engineers during last decades. In the proposed approach to machine learning for model order reduction, all this knowledge is preserved in local ROMs. It is even more than that, the uncertainty propagation comes to complete this valuable traditional knowledge. We do not expect from artificial intelligence to learn everything in our modeling process.

3.2 Industrial dataset and objectives

3.2.1 Industrial problem

Figure 2 gives the geometry and the finite-element mesh of a real high-pressure turbine blade. The mesh is made of quadratic tetrahedral elements, and contains a number of nodes in the order of the million. The elasto-viscoplastic mechanical behavior is described by a crystal plasticity model presented in Section A. As explained above, Monte Carlo simulations using a commercial software as Abaqus are unaffordable. With the help of domain decomposition methods, the computation time can be reduced by solving equilibrium equations in parallel on different subdomains of the geometry. Using the implementation of the Adaptive MultiPreconditioned FETI solver [3] in Z-set finite-element software [4], the simulation of one single loading cycle of the HP turbine blade with 48 subdomains takes approximately 53 minutes.

thumbnail Fig. 2

High-pressure turbine blade geometry and mesh (micro-perforations are not modeled).

3.2.2 Objectives

The objective is to use a ROM-net to quantify uncertainties on the mechanical behavior of the high-pressure turbine blade, givenuncertainties on the thermal loading. The reduction of the computation time should enable Monte Carlo simulations for uncertainty quantification. In particular, we are not interested in predicting the state of the structure after a large number of flight-representative loading cycles. Only one cycle is simulated. Cyclic extrapolation of the behavior of a high-pressure turbine blade has been studied in [22] and is out of the scope of the present work.

3.3 Modeling assumptions

3.3.1 Weak thermomechanical coupling

It is assumed that the heat produced or dissipated by mechanical phenomena has negligible effects in comparison with thermal conduction, which enables avoiding strongly coupled thermomechanical simulations and running thermal and mechanical simulations separately instead. Under a weak thermomechanical coupling, the first step consists in solving the heat equation to determine the temperature field and its evolution over time. The temperature field history defines the thermal loading and is used to compute thermal strains and temperature-dependent material parameters for the mechanical constitutive laws. Once the thermal loading is known, the temperature-dependent mechanical problem must be solved in order to predict the mechanical response of the structure.

3.3.2 Cyclic thermomechanical loading

The thermomechanical loading applied to the high-pressure turbine blade during its whole life is modeled as a cyclic loading, with one cycle being equivalent to one flight. The rotation speed of the turbine’s rotor is proportional to a periodic function of time ω(t) whose evolution over one period (or cycle, see Figure 3) is representative of one flight with its three main regimes, namely take-off, cruise, and landing. The period (or duration of one cycle) is denoted by tc. The rotation speed between flights k and k + 1 is zero, which means that ω(ktc) = 0 for any integer k. The rotation speed ω(t) is scaled so that its maximum is 1.

Let denote the solid body representing the high-pressure turbine blade, with ∂Ω denoting its outer surface. Let ∂Ωp ⊂ ∂Ω be the surface corresponding to the intrados and extrados. The thermal loading is defined as: (4)

where T0 = 293 K and Tmax is the temperature field obtained when the rotation speed reaches its maximum. This field Tmax is obtained either by an aerothermal simulation or by a stochastic model, as explained later. Similarly, the pressure load applied on ∂Ωp reads: (5)

where is the atmospheric pressure at sea level, and where is the pressure field obtained when the rotation speed reaches its maximum. The clamping of the blade’s fir-tree foot on the rotor disk is modeled by displacements boundary conditions that are not detailed here.

thumbnail Fig. 3

Function ω(t) defining onecycle for the rotation speed.

3.3.3 Geometric details and thermal barrier coating

Small geometric details of the structure have been removed to simplify the geometry. Nonetheless, the main cooling channels are considered. The effects of the thermal barrier coating (TBC) have been integrated in aerothermal simulations, but the TBC is not considered in the mechanical simulation although its damage locally increases the temperature in the nickel-based superalloy and thus affects the fatigue resistance of the structure. Additional centrifugal effects due to the TBC are not taken into account.

3.3.4 Influential factors

The predicted mechanical response of the structure depends on many different factors. Below is a nonexhaustive list of influential factors that are possible sources of uncertainties in the numerical simulation:

  • Thermal loading: The viscoplastic behavior of the nickel-based superalloy is very sensitive to the temperature field and its gradients. However, the temperature field is not accurately known because of the impossibility of validating numerical predictions experimentally. Indeed, temperature-sensitive paints are accurate to within 50 K only, and they do not capture a real surface temperature field since they measure the maximum temperature reached locally during the experiment.

  • Crystal orientation: Because of the complexity of the manufacturing process of monocrystalline blades, the orientation of the crystal is not perfectly controlled. As the superalloy has anisotropic mechanical properties, defaults in crystal orientation highly affect the location of damaged zones in the structure.

  • Mechanical loading: The centrifugal forces are well known because they are related to the rotation speed that is easy to measure. On the contrary, pressure loads are uncertain because of the turbulent nature of the incoming fluid flow. However, the effects of pressure loads uncertainties on the mechanical response are less significant than those of the thermal loading and crystal orientation uncertainties.

  • Constitutive laws: Uncertainties on the choice of the constitutive model, the relevance of the modeling assumptions, and the values of the calibrated parameters involved in the constitutive equations also influence the results of the numerical simulations.

For simplification purposes, the only source of uncertainty that is considered in this work is the thermal loading. The equations of the mechanical problem are then seen as parametrized equations, where the parameter is the temperature field Tmax (see Eq. (4)) obtained when the rotation speed reaches its maximum value. The dimension of the parameter space is then the number of nodes in the finite-element mesh. The mechanical loading is assumed to be deterministic. With the crystal orientation, the constitutive laws and their parameters (or coefficients), they are considered as known data describing the context of the studyand given by experts.

3.4 Stochastic model for the thermal loading

A stochastic model is required to take into account the uncertainties on the thermal loading. Given the definition of the thermal loading in equation (4), we only need to model uncertainties in space through the field Tmax obtained when the rotation speed reaches its maximum value. The random temperature fields must satisfy some constraints: they must satisfy the heat equation, and they must not take values out of the interval [0 K; Tmelt], where Tmelt is the melting point of the superalloy. These random fields are obtained by adding random fluctuations to a reference temperature field, see Figure 4. The reference field and comes from aerothermal simulations run with the software Ansys Fluent1. The data-generating distribution is defined as a Gaussian mixture model made of two Gaussian distributions with the same covariance function but with distinct means, and with a prior probability of 0.5 for each Gaussian distribution. The Gaussian distributions are obtained by taking the four first eigenfunctions of the covariance function (see Karhunen-Loève expansion [64]), with a standard deviation of 15 K. Therefore, realizations of the random temperature field read: (6)

where Tref is the reference field, δT0 is a temperature perturbation at the trailing edge whose maximum value is 50 K, are fluctuation modes, ϒ0 is a random variable following the Bernoulli distribution with parameter 0.5, and are independent and identically distributed random variables following the standard normal distribution . The variable ϒ0 is also independent of the other variables ϒi. The different fields involved in equation (6) can be visualized in Figure 4. Equation (6) defines a mixture distribution with two Gaussian distributions whose means are Tref and Tref + δT0. We voluntarily define this mixture distribution with δT0 adding 50 K in a critical zone of the turbine blade in order to check that our cluster analysis can successfully detect two relevant clusters, i.e. one for fields obtained with ϒ0 (θ) = 0 and one for fields obtained with ϒ0(θ) = 1. Indeed, the temperatureperturbation δT0 is expected to significantly modify the mechanical response of the high-pressure turbine blade. All the fields satisfy the steady heat equation like Tref, which ensures that the random fields always satisfy the heat equation under the assumption of a linear thermal behavior. For nonlinear thermal behaviors, equation (6) would define surface temperature fields that would be used as Dirichlet boundary conditions for the computation of bulk temperature fields. The assumption of a linear thermal behavior is adopted here to avoid solving the heat equation for every realization of the random temperature field.

Let us now give more details about the construction of the fluctuation modes . First, surface fluctuation modes are computed on the boundary ∂Ω using the method given in [65] for the construction of random fields on a curved surface. The correlation function is defined as a function of the geodesic distance dG along the surface ∂Ω: (7)

where is a correlation length. Geodesic distances are computed using the algorithm described in [66,67] and implemented in the Python library gdist2. A covariance matrix is built by evaluating the correlation function on pairs of nodes of the outer surface of the finite-element mesh, and multiplying the correlation by the constant variance. The four surface modes are then obtained by finding the four eigenvectors corresponding to the largest eigenvalues of the covariance matrix. The steady heat equation with Dirichlet boundary conditions is solved for each of these surface modes to derive the 3D fluctuation modes, using Z-set [4] finite-element solver. The Python library BasicTools3 developed by SafranTech is used to read the finite-element mesh and write the temperature fields in a format that can be used for simulations on Z-set.

thumbnail Fig. 4

Reference temperature field (on the left), temperature perturbation at the trailing edge (field 0 = δT0), and fluctuation modes (fields 1 to 4). The fluctuations in the fourth mode are located inside the blade, in the cooling channels.

4 ROM-net based uncertainty quantification applied to an industrial high-pressure turbine blade

This section develops the different stages of the ROM-net for the industrial test case presented in the previous section. Given our budget of 200 high-fidelity simulations, a dictionary containing two local ROMs is constructed using our clustering procedure. A logistic regression classifier is trained for automatic model recommendation using information identified by feature selection, followed by an alternative to the Gappy-POD for full-field reconstruction of dual quantities. Then, the results of the uncertainty quantification procedure are presented. We advise the reader to refer to Section 4.6, which contains an illustration of the proposed workflow, while reading Sections 4.1-4.5. Finally the accuracy of the ROM-net is validated using simulations for new temperature loadings.

4.1 Design of numerical experiments

Given the computational cost of high-fidelity mechanical simulations of the high-pressure turbine blade, the training data are sampled from the stochastic model for the thermal loading using a design of experiments (DoE). Our computational budget corresponds to 200 high-fidelity simulations, so a database of 200 temperature fields must be built. This database includes two separate datasets coming from two independent DoEs:

  • The first dataset is built from a Maximum Projection LHS design (MaxProj LHS DoE [68]) and contains 80 points. This dataset will be used for the construction of the dictionary of local ROMs via clustering. The MaxProj LHS DoE has good space-filling properties on projections onto subspaces of any dimension.

  • The second dataset is built from a Sobol’ sequence (Sobol’ DoE) of 120 points. Using a suboptimal DoE method ensures that this second dataset is different and independent from the first one. The lower quality of this dataset with respect to the first one is compensated by its larger population. This dataset will be used for learning tasks requiring more training examples than the construction of the local ROMs, namely the classification task for automatic model recommendation, and the training of cluster-specific surrogate models for the reconstruction of full fields from hyper-reduced predictions on a reduced-integration domain. These surrogate models (Gappy surrogates) replace the Gappy-POD [69] method that is commonly used in hyper-reduced simulations to retrieve dual variables on the whole mesh.

These DoEs are built with the platform Lagun4. The fact that these two datasets come from two separate DoEs is beneficial: as each of them is supposed to have good space-filling properties, they are both representative of the possible thermal loading and can therefore be used to define a training set and a test set for a given learning task. For instance, the classifier trained on the Sobol’ DoE can be tested on the MaxProj LHS DoE. The local ROMs built from snapshots belonging to the MaxProj LHS DoE can make predictions on the Sobol’ DoE that will be used for the training of the Gappy surrogates, which is relevant since the Gappy surrogates are supposed to analyze ROM predictions on new unseen data in the exploitation phase.

Drawing random temperature fields as defined in equation (6) requires sampling data from the random variables , where ϒ0 follows the Bernoulli distribution with parameter 0.5 and the variables ϒi for i ∈ [ [1;4] ] are independent standard normal variables and independent of ϒ0. Both DoE methods (Maximum Projection LHS and Sobol’ sequence) generate point clouds with a uniform distribution in the unit hypercube. Figures 5 and 6 show the projections onto 2-dimensional subspaces of the 5D point clouds used to build our datasets. The marginal distributions are plotted to check that they well approximate the uniform distribution. These point clouds, considered as samples of a random vector (χ0, χ1, χ2, χ3, χ4) following the uniform distribution on the unit hypercube, are transformed into realizations of the random vector (ϒ0, ϒ1, ϒ2, ϒ3, ϒ4) using the following transformations: (8)

where F−1 is the inverse of the cumulative distribution function of the standard normal distribution. The resulting samples define the MaxProj dataset and the Sobol’ dataset of random temperature fields, using equation (6). Each temperature field defines a thermal loading, using equation (4). The 200 corresponding mechanical problems are solved for one loading cycle with the finite-element software Z-set [4] with the domain decomposition method described in [3], with 48 subdomains. The average computation time for one simulation is 53 minutes.

thumbnail Fig. 5

Visualization of the MaxProj LHS DoE. The marginal distributions are represented on the diagonal. The 5D DoE is projected on 2D subspaces for visualization purposes, in order to check space-filling properties in 2D.
thumbnail Fig. 6

Visualization of the Sobol’ DoE. The marginal distributions are represented on the diagonal. The 5D DoE is projected on 2D subspaces for visualization purposes, in order to check space-filling properties in 2D.

4.2 ROM dictionary construction

4.2.1 Clustering

The 80 simulations associated to the MaxProj dataset are used as clustering data. Loading all the simulation data and computing the pairwise ROM-oriented dissimilarities takes about 5 minutes. The ROM-oriented dissimilarity defined in [61, Definition 3.11], and mentioned in Section 2.2.2, is computed with n = 1, i.e. each simulation is represented by one field. The dataset is partitioned into two clusters using our implementation of PAM [70,71] k-medoids algorithm, with 10 different random initializations for the medoids. The clustering results can be visualized using Multidimensional Scaling (MDS) [72]. MDS is an information visualization method which consists in finding a low-dimensional dataset Z0 whose matrix of Euclidean distances d(Z0) is an approximation of the true dissimilarity matrix δ. To that end, a cost function called stress function is minimized with respect to Z : (9)

This minimization problem is solved with the algorithm Scaling by MAjorizing a COmplicated Function (SMACOF, [73]) implemented in Scikit-Learn [74]. Figure 7 show the clusters on the MDS representations with the goal-oriented variants of the ROM-oriented dissimilarity measure, applied on the accumulated plastic strain (see the Appendix for details on the constitutive model). The clustering results is compared with the expected clusters corresponding to ϒ0 = 0 and ϒ0 = 1, the latter corresponds to the perturbation δT0 being activated. The obtained clusters almost correspond to the expected ones, with only 4 points with wrong labels out of 80, which quantifies the ability of the ROM-oriented dissimilarity measure on the accumulated plasticity to infer the correct value of ϒ0. The medoids of the two clusters are given in Figure 8. Cluster 0 contains temperature fields for which ϒ0 = 1, while cluster 1 contains fields for which ϒ0 = 0. It can be observed that the quantity of interest clearly differs from one cluster to the other, while the differences are hardly visible on the displacement field. The displacement field combines deformations associated to different phenomena (thermal expansion, elastic strains, viscoplastic strains) that are not necessarily related to damage in the structure, which could explain why the quantity of interest seems to be more appropriate for clustering in this example.

thumbnail Fig. 7

MDS representation of the clustering results using the ROM-oriented dissimilarity measure on the quantity of interest (goal-oriented variant). On the left, the colors correspond to the expected clusters. On the right, the colors correspond to the clusters identified by the clustering algorithm. The positions of the labels 0 and 1 coincide with the positions of the clusters’ medoids. The MDS relative error ς(Z0 ; δ)∕ς(0;δ) is 12%.
thumbnail Fig. 8

The 3 fields on the left correspond to the medoid of cluster 0, and those onthe right correspond to the medoid of cluster 1. The fields in the first and the third columns show the differences between the medoids’ temperature fields and the reference temperature field Tref (the scale istruncated for the first field). The second and the fourth columns show the displacement magnitude field (top) and the quantity of interest (bottom).

4.2.2 Construction of local ROMs

The simulations used for the clustering procedure can directly provide snapshots for the construction of the local ROMs. To control the duration of their training, only 20 simulations are selected to provide snapshots for the each local ROMs, which represents half of the clusters’ populations. These simulations are selected in a maximin greedy approach starting from the medoid (see [62, Algorithm 2, Stage 2] for a example of maximin selection). Figure 9 shows which simulations have been selected for the construction of the local ROMs.

The local ROMs are built following the methodology described in Section 2.1. The snapshot-POD and the ECM are done in parallel with shared memory on 24 cores. The tolerance for the snapshot-POD is set to 10−8 for the displacement field, and to 10−4 for dual variables (the quantity of interest and the six components of the stress tensor). The POD bases for the dual variables will be used for their reconstruction with the Gappy surrogates. The tolerance for the ECM is set to 5 × 10−4. Locals ROMs each contain 18 displacement modes and between 8 and 13 modes for stress components. The first and second local ROMs contain respectively 10 and 12 modes for the quantity of interest . The ECM selects 506 (resp. 510) integration points for the reduced-integration domain of ROM 0 (resp. 1). Building one local ROM takes approximately 2 hours and 30 minutes.

thumbnail Fig. 9

MDSrepresentation of the clustering results. Orange points represent the snapshots selected for cluster 0, while the light blue points represent the snapshots selected for cluster 1. For each cluster, the snapshots are selected by a maximin procedure starting from the medoid.

4.3 Automatic model recommendation

In this section, a classifier is trained for the automatic model recommendation task. The 120 temperature fields coming from the Sobol’ dataset are used as training data for the classifier. Their labels are determined by finding their closest medoid in terms of the ROM-oriented dissimilarity measure. Hence, for each temperature field of the Sobol’ dataset, twodissimilarities are computed: one with the medoid of the first cluster, and one with the medoid of the second cluster. Once trained, the classifier can be evaluated on the 80 labelled temperature fields of the MaxProj dataset.

4.3.1 Feature selection

Each temperature field is discretized on the finite-element mesh, which contains in the order of the million nodes. To reduce the dimension of the input space and facilitate the training phase of the classifier, we apply the geostatistical mRMR feature selection algorithm described in [62, Algorithm 1] on data from the Sobol’ dataset. First, 800 pairs of nodes are selected in the mesh, which takes 18 seconds. The 800 corresponding redundancy terms are computed with Scikit-Learn [74] in less than 3 seconds. Figure 10 plots the values of these redundancy terms versus the Euclidean distance between the nodes. We observe that the correlation between the redundancy mutual information terms and the distance between the nodes is poor, with a lot of noise. This can be due to the fact that the random temperature fields have been built using Gaussian random fields on the outer surface with an isotropic correlation function depending on the geodesic distance along the surface rather than the Euclidean distance. Since the turbine blade is a relatively thin structure, two nodes, one on the intrados and another one on the extrados, can be close to each other in the Euclidean distance, but with totally uncorrelated temperature fluctuations because of the large geodesic distance separating them. On the contrary, two points on the same side of the turbine blade can have correlated temperature variations while being separated by a Euclidean distance in the order of the blade’s thickness. The length of the mutual information’s high-variance regime seems to correspond to the blade’s chord, which supports this explanation. The thinness of the turbine blade induces anisotropy in the correlation function of the bulk Gaussian random field defining the thermal loading, which implies an anisotropic behavior of the mutual information according to [62, Property 1]. The use of a local temperature perturbation δT0 in conjunction with fluctuation modes having larger length scales may also partially explain the large variance of redundancy terms. Nonetheless,it remains clear that redundancy terms are smaller as for large distances. This trend is captured by a kriging metamodel (Gaussian process regression) trained with Scikit-Learn in a few seconds, with a sum-kernel involving the Matérn kernel with parameter 5∕2 (to get a continuous and twice differentiable metamodel) and length scale 1, and a whitekernel to estimate the noise level of the signal. The curve of the metamodel is given in Figure 10. Then, for each node of the finite-element mesh, the mutual information with the label variable is computed. The computations of these relevance terms (in the order of the million terms) are distributed between 280 cores, which gives a total computation time of 15 minutes. Among these features, 5, 986 features are preselected by discarding those with a relevance mutual information lower than 0.05. The geostatistical mRMR selects 11 features in 42 seconds. The corresponding nodes in the finite-element mesh can be visualized in Figure 10.

Remark 2.

The metamodel for redundancy terms could be improved by defining it as a function of the precomputed geodesic distances along the outer surface rather than the Euclidean distances. Each finite-element node would be associated to its nearest neighbor on the outer surface before computing the approximate mutual information from geodesic distances.

thumbnail Fig. 10

Feature selection results. The kriging metamodel for redundancy terms is represented by the red curve and built from 800 true redundancy terms (blue points). The elements containing the selected nodes are represented in the turbine blade geometry.

4.3.2 Classification

The classifier is trained on the Sobol’ dataset, using the values of the temperature fields at the 11 nodes identified by the feature selection algorithm. The classifier is a logistic regression [75–77] with elastic net regularization [78] implemented in Scikit-Learn. The two hyperparameters involved in the elastic net regularization are calibrated using 5-fold cross-validation, giving a value of 0.001 for the inverse of the regularization strength, and 0.4 for the weight of the L1 penalty term (and thus 0.6 for the L2 penalty term). Due to the L1 penalty term, the classifier only uses 5 features among the 11 input features. The classifier’s accuracy, evaluated on the MaxProj dataset to use new unseen data, reaches 98.75%. The confusion matrix indicates that 100% of the test examples belonging to class 0 have been correctly labeled, and that 2.38% of the test examples belonging to class 1 have been misclassified. Table 1 summarizes the values of precision, recall and F1-score on test data.

Table 1

Classification results.

4.4 Surrogate model for Gappy reconstruction

When using hyper-reduction, the ROM calls the constitutive equations solver only at the integration points belonging to the reduced-integration domain. It is recalled that the ECM selected 506 (resp. 510) integration points for the reduced-integration domain of ROM 0 (resp. 1), and that the finite-element mesh initially contains a number of integration points in the order of the million. Therefore, after a reduced simulation, dual variables defined at integration points are known only at integration points of the reduced-integration domain. To retrieve the full field, the Gappy-POD [69] finds the coefficients in the POD basis that minimize the squared error between thereconstructed field and the ROM predictions on the reduced-integration domain. This minimization problem defines the POD coefficients as a linear function of the predicted values on the reduced-integration domain. Although these coefficients are optimal in the least squares sense, they can be biased by the errors made by the ROM. To alleviate this problem, we propose to replace the common Gappy-POD procedure by a metamodel or Gappy surrogate. The inputs and the outputs of the Gappy surrogate are the same as for the Gappy-POD: the input is a vector containing the values of a dual variable on the reduced-integration domain, and the output is a vector containing the optimal coefficients in the POD basis. One Gappy surrogate must be built for each dual variable of interest: in our case, 7 surrogate models per cluster are required, namely one for the quantity of interest and one for every component of the Cauchy stress tensor.

The training data for these Gappy surrogates are obtained by running reduced simulations with the local ROMs, using the thermal loadings of the Sobol’ dataset. Indeed, the two local ROMs have been built on the MaxProj dataset, therefore thermal loadings of the Sobol’ dataset can play the role of test data for the ROMs. For each thermal loading in the Sobol’ dataset, the true high-fidelity solution is already known since it has been computed to provide training data for the classifier. In addition, the exact labels for these thermal loadings are known, which means that we know which local ROM to choose for each thermal loading of the Sobol’ dataset. Given ROM predictions on the reduced-integration domain, the optimal coefficients in the POD basis are given by the projections of the true prediction made by the high-fidelity model (the finite-element model) onto the POD modes. This provides the true outputs for the Sobol’ dataset, which can then be used as a training set for the Gappy surrogates.

Given the high-dimensionality of the input data (there are more than 500 integration points in the reduced-integration domains) with respect to the number of training examples (120 examples), a multi-task Lasso metamodel is used. The hyperparameter controlling the regularization strength is optimized by 5-fold cross-validation. Training the 14 Gappy surrogates (7 for each cluster) takes 1 hour. The Gappy surrogates select between 8% and 18% of the integration points in the reduced-integration domains, due to the L1 regularization. The mean cross-validated coefficients of determination are 0.9637 (resp. 0.8935) for the quantity of interest for cluster 0 (resp. cluster 1), and range from 0.9404 to 0.9938 for stress components. These satisfying results mean that it is not required to train a kriging metamodel with the variables selected by Lasso to get nonlinear Gappy surrogates. The Gappy surrogates are then linear, just as the Gappy-POD.

The accuracy gains provided by the Gappy surrogates with respect to classical Gappy-POD on the present industrial case is investigated in Figure 11. Here, 24 high-fidelity simulation in the first cluster are computed (with ϒ0 = 0) as reference, and Gappy surrogates (using two meta models Lasso and ElasticNet) and classical Gappy-POD are computed using ROM 0. For both variants of meta models and both quantities of interest: accumulated plastic strain and the component 33 of the stress tensor, Gappy surrogates provides more accurate predictions than the classical Gappy-POD.

Remark 3.

In this strategy, the local ROMs solve the equations of the mechanical problem, which enables using linear surrogate models to reconstruct dual variables. Using surrogate models from scratch instead of local ROMs would have been more difficult, given the nonlinearities of this mechanical problem and the lack of training data for regression. In addition, such surrogate models would require a parametrization of the input temperature fields, whereas the local ROMs use the exact values of the temperature fields on the RID without assuming any model for the thermal loading.

The dictionary-based ROM-net used for mechanical simulations of the high-pressure turbine blade is made of a dictionary of two local hyper-reduced order models and a logistic regression classifier. The classifier analyzes the values of the input temperaturefield at 11 nodes only, identified by our feature selection strategy. For a given thermal loading in the exploitation phase, after the reduced simulation with the local ROM recommended by the classifier, linear cluster-specific Gappy surrogates reconstruct the full dual fields (quantity of interest and stress components) from their predicted values on the reduced-integration domain.

thumbnail Fig. 11

Mean over the complete mesh of dual quantities of interest: accumulated plastic strain (left) and component 33 of the stress tensor σ33 (right), plotted as points where the x-coordinate is the reference value, and the y-coordinate is the considered reduced prediction.

4.5 Uncertainty quantification results

Once trained, the ROM-net can be applied for the quantification of uncertainties on the mechanical behavior of the HP turbine blade resulting from the uncertainties on the thermal loading. Since the ROM-net online operations can be performed sequentially on one single core, 24 cores are used in order to compute the solution for 24 thermal loadings at once. This way, 42 batches of 24 Monte Carlo simulations are run in 2 hours and 48 minutes. The 1008 thermal loadings used for this study are generated by randomly sampling points from the uniform distribution on the 5D unit hypercube and applying the transformation given in equation (8).

The expected values of and are estimated with the empirical means , where Zi are the corresponding samples. The variances of and are computed using the unbiased sample variance . The Central Limit Theorem gives asymptotic confidence intervals for the expected values: for all α ∈]0;1[, (10)

where ϕr denotes the quantile of order r of the standard normal distribution , and In is an asymptotic confidence interval with confidence level 1 − α for the expectation μ: . The widths of the confidence intervals are expressed as a percentage of the estimated value for the expectations in Table 2.

The probability density functions of the quantities of interest can be estimated using Gaussian kernel density estimation (see Section 6.6.1. of [79]). Figure 12 gives the histograms and estimated distributions for and . The shapes of these distributions highly depend on the assumptions made for the stochastic thermal loading. As observed in Figure A.1, the stress field is highly sensitive to temperature gradients, which may explain why the distribution of the Von Mises stress is bimodal.

Table 2

Widths of the confidence intervals (CI) for the expectations, expressed as percentages of the estimated expectations.

thumbnail Fig. 12

Histograms and probability density functions of the quantities of interest (left) and (right).

4.6 Workflow

Figure 13 provides an illustration of the workflow and the computational time of each step presented in Sections 4.1-4.5.

thumbnail Fig. 13

Workflow for the ROM-net methology applied to the considered industrial setting in Sections 4.1-4.5.

4.7 Verification

For verification purposes, the accuracy of the ROM-net is evaluated on 20 Monte Carlo simulations with 20 new thermal loadings. These thermal loadings are generated by randomly sampling points from the uniform distribution on the 5D unit hypercube, and applying the transformation given in equation (8). The reduced simulations are run on single cores. The total computation time for generating a new thermal loading on the fly, selecting the corresponding reduced model, running one reduced simulation and reconstructing the quantities of interest is 4 minutes on average. As a comparison, one single high-fidelity simulation with Z-set [4] with 48 subdomains takes 53 minutes, which implies that the ROM-net computes 13.25 times faster. However, one high-fidelity simulation requires 48 cores for domain decomposition, whereas the ROM-net works on one single core. Hence, using 48 cores to run 48 reduced simulations in parallel, 636 reduced simulations can be computed in 53 minutes with the ROM-net, while the high-fidelity model only runs one simulation. In addition to the acceleration of numerical simulations, energy consumption is reduced by a factor of 636 in the exploitation phase. In spite of the fast development of high-performance computing, numerical methods computing approximate solutions at reduced computational resources and time are particularly important for many-query problems such as uncertainty quantification, where the intensive use of computational resources is a major concern. Model order reduction and ROM-nets play a prominentrole toward green numerical simulations [80]. Of course, the number of simulations in the exploitation phase must be large enough to compensate the efforts made in the training phase, like in any machine learning or model order reduction problem.

Figures 14 and 15 show the results for two simulations belonging to cluster 0 and cluster 1 respectively.These figures give the difference between the current temperature field as the reference one, i.e. the field T −Tref, and the resulting variations of the quantity of interest predicted by the ROM-net and the high-fidelity model, i.e. and . The signs and the positions of the variations of the quantity of interest seem to be quite well predicted by the ROM-net.

Let us introduce a zone of interest Ω′ defined by all of the integration points at which is higher than for the thermal loading defined by Tref + δT0. This zone of interest contains 209 integration points. The values of the variables and σeq averaged over Ω′ are denoted by and . Table 3 gives different indicators quantifying the errors made by the ROM-net: the L2 relative errors on the whole domain Ω and on the zone of interest Ω′, the L∞ relative errors on Ω and Ω′, the relative errors on and , and the errors on the locations of the pointswhere the fields and σeq reach their maxima. All the relative errors remain in the order of 1% or 2%, which validates the methodology. In addition, the ROM-net perfectly predicts the position of the critical points at which and σeq reach their maxima. Figure 16 shows errors on the quantities of interest.

thumbnail Fig. 14

Comparison between high-fidelity predictions (middle column) and ROM-net’s predictions (right-hand column). The field on the left represents the difference between the current temperature field (belonging to cluster 0) and the reference one. The other fields correspond to the increments of the quantity of interest with respectto its reference state obtained with the reference temperature field.
thumbnail Fig. 15

Comparison between high-fidelity predictions (middle column) and ROM-net’s predictions (right-hand column). The field on the left represents the difference between the current temperature field (belonging to cluster 1) and the reference one. The other fields correspond to the increments of the quantity of interest with respectto its reference state obtained with the reference temperature field.
Table 3

Error indicators for the evaluation of the ROM-net on 20 new thermal loadings.

thumbnail Fig. 16

Errors on the quantity of interest . The red (resp. blue) color is used for zones where the quantity of interest is overestimated (resp. underestimated).

5 Conclusion

In this work, we used a dictionary-based ROM-net to successfully quantify the uncertainties on dual quantities of interest of an elastoviscoplastic high-pressure turbine blade, generated by the uncertainty of its temperature loading. This validates the methodology on large models with highly nonlinear behaviors. An outlook of this work would be to consider nonparametrized geometrical variability, which is of paramount interest when considering the design of mechanical parts and the uncertainty quantification of their manufacturing processes.

Acknowledgements

Aknowledgments The authors wish to thank Sébastien Da Veiga and Clément Bénard for fruitful discussions on uncertainty quantification.

Funding: This research was partially funded by the French Fonds Unique Interministériel (MOR_DICUS).

Appendix Mechanical constitutive model

It is assumed that the mechanical behavior of the high-pressure turbine blade can be described in the framework of the infinitesimal strain theory. The mechanical response of the structure during the first loading cycle is described by the following equilibrium equations and boundary conditions: (A.1)

where u(ξ, t) is the displacement field (primal variable), σ(ξ, t) is the symmetric second-order Cauchy stress tensor, fC(ξ, t) is the local volumic centrifugal force, u∂Ω(ξ, t) are the imposed displacements, and n(ξ, t) is the outward-pointing normal vector to the outer surface ∂Ω. The relationbetween the stress tensor and the displacement field is described by constitutive laws modeling the mechanical behavior of the monocrystalline nickel-based superalloy. At high temperatures, this material has an elasto-viscoplastic behavior that can be described in the crystal plasticity framework [81,82] to model inelastic strains generated by the motion of dislocations5 in different slip systems of the crystal. The strain tensor ε is defined as the symmetric part of the displacement gradient (with respect to ξ): (A.2)

The stress tensor is obtained from the elastic strain tensor using Hooke’s law: (A.3)

where εp is the tensor of inelastic strains and 1 is the identity second-order tensor. The fourth-order tensor C is the stiffness tensor. Given the face-centered cubic crystal structure of the superalloy, the stiffness tensor is anisotropic but has only three independentcoefficients. The thermal expansion of crystals with cubic symmetry is isotropic, which explains why the thermal expansion coefficient α is the same in all directions. The time evolution of hidden variables such as inelastic strains are described by ordinary differential equations that must be solved at every integration point of the finite-element mesh. The inelastic strain rate can be decomposedinto contributions of dislocations motions in 12 octahedral slip systems and 6 cubic slip systems: (A.4)

where (resp. ) is the shear strain rate in the s-th octahedral (resp. cubic) slip system. The tensor (resp. ) is the orientation tensor of the s-th octahedral (resp. cubic) slip system, defined by the normal (resp. ) to the slip plane and the slip direction (resp. ). The shear strain rates are given by a hyperbolic viscoplastic flow rule: (A.5)

where , and are materialparameters. Similar equations are satisfied in cubic slip systems. The resolved shear stresses are given by Schmid’s law: (A.6)

Again, similar equations are valid for cubic slip systems. The stress variables , , and describe hardening phenoma, i.e. the evolution of the shape of the elastic domain within which no dissipative phenoma occur. The back-stresses (and ) are the solutions of an ordinary differential equation modeling kinematic hardening with static recovery: (A.7)

Isotropic hardening is modeled by the following equations: (A.8)

with . All the constitutive equations given in this section are true for all ξ ∈ Ω and for all t ∈ [0;tc], and are solved at every integration point of the finite-element mesh. All the coefficients involved in these equations depend on the local value of the temperature field. The problem is thus seen as a system of partial differential equations and ordinary differential equations parametrized by the thermal loading. The standard procedure for the computation of a fatigue lifetime with an uncoupled damage model consists in solving the mechanical problem for a large number of cycles until the stabilization of the mechanical response (in the case of plastic shakedown). Then, a damage field can be computed in a post-processing step and can be linked to a fatigue lifetime. For high-pressure turbine blades, fatigue models generally consider interaction effects with oxidation and creep, like in [83,84]. In this work, no fatigue lifetime is computed since we only solve the problem for the very first cycle. Instead, our quantity of interest is a strain indicator that partially describes the damage state of the material. This quantity of interest corresponds to the accumulated plastic strain in octahedral slip systems at the end of the first cycle, which reads: (A.9)

with: (A.10)

For simplicity, among other possible measures of the stress tensor, we consider as quantity of interest the von Mises equivalent stress field, defined as: (A.11)

Therefore, the variables considered for the evaluation of the ROM-net and for uncertainty quantification are the accumulated plastic strain in octahedral slip systems at the end of the first cycle, and the von Mises stress σeq obtained when the rotation speed reaches its maximum value. These variables can be visualized in Figure A.1 for a reference thermal loading.

thumbnail Fig. A.1

On the left: von Mises stress field σeq obtained when the rotation speed reaches its maximum value. On the right: accumulated plastic strain in octahedral slip systems at the end of the first cycle. Note: the foot of the high-pressure turbine blade has an elastic behavior, while the rest of the blade has a viscoplastic behavior described by a crystal plasticity model.

References

  1. Dassault Systèmes. Abaqus Unified FEA. https://www.3ds.com/fr/produits-et-services/simulia/produits/abaqus/, 1978-present [Google Scholar]
  2. T. Daniel, F. Casenave, N. Akkari, D. Ryckelynck, Model order reduction assisted by deep neural networks (ROM-net), Adv. Model. Simul. Eng. Sci. 7, 1–27 (2020) [CrossRef] [Google Scholar]
  3. C. Bovet, A. Parret-Freaud, N. Spillane, P. Gosselet, Adaptive multipreconditioned FETI: scalability results and robustness assessment, Comput. Struct. 193, 1–20 (2017) [CrossRef] [Google Scholar]
  4. Mines ParisTech and ONERA the French aerospace lab. Z-set: nonlinear material & structure analysis suite. http://www.zset-software.com, 1981-present [Google Scholar]
  5. L. Cordier, M. Bergmann, Proper Orthogonal Decomposition: an overview, in Lecture series 2002-04, 2003-03 and 2008-01 on post-processing of experimental and numerical data, Von Karman Institute for Fluid Dynamics, 2008, page 46 pages. VKI (2008) [Google Scholar]
  6. C. Rowley, T. Colonius, R. Murray, Model reduction for compressible flow using POD and Galerkin projection, Physica D 189, 115–129 (2003) [Google Scholar]
  7. J.A. Hernandez, M.A. Caicedo, A. Ferrer, Dimensional hyper-reduction of nonlinear finite element models via empirical cubature, Comput. Methods Appl. Mech. Eng. 313, 687–722 (2017) [CrossRef] [Google Scholar]
  8. W. Keiper, A. Milde, S. Volkwein, Reduced-Order Modeling (ROM) for Simulation and Optimization: Powerful Algorithms as Key Enablers for Scientific Computing (Springer International Publishing, 2018) [CrossRef] [Google Scholar]
  9. A. Quarteroni, G. Rozza, Reduced Order Methods for Modeling and Computational Reduction (Springer Publishing Company, Incorporated, 2013) [Google Scholar]
  10. F. Chinesta, E. Cueto, PGD-Based Modeling of Materials, Structures and Processes (Springer, 2014) [CrossRef] [Google Scholar]
  11. F. Chinesta, P. Ladeveze, E. Cueto, A short review on model order reduction based on Proper Generalized Decomposition, Arch. Comput. Methods Eng. 18, 395–404 (2011) [CrossRef] [Google Scholar]
  12. C. Prud’homme, D. Rovas, K. Veroy, L. Machiels, Y. Maday, A. Patera, G. Turinici, Reliable real-time solution of parametrized partial differential equations: reduced-basis output bound methods, J. Fluids Eng. 124, 70 (2002) [CrossRef] [Google Scholar]
  13. G. Rozza, D. Huynh, A. Patera, Reduced basis approximation and a posteriori error estimation for affinely parametrized elliptic coercive partial differential equations, Arch. Comput. Methods Eng. 15, 1–47 (2007) [CrossRef] [Google Scholar]
  14. D. Ryckelynck, A priori hyperreduction method: an adaptive approach, J. Comput. Phys. 202, 346–366 (2005) [CrossRef] [Google Scholar]
  15. M. Barrault, Y. Maday, N.C. Nguyen, A.T. Patera, An empirical interpolation method: application to efficient reduced-basis discretization of partial differential equations, Comptes Rendus Math. 339, 666–672 (2004) [Google Scholar]
  16. P. Astrid, S. Weiland, K. Willcox, T. Backx, Missing point estimation in models described by proper orthogonal decomposition, Proc. IEEE Conf. Decis. Control 53, 1767–1772 (2005) [Google Scholar]
  17. N.C. Nguyen, A.T. Patera, J. Peraire, A best points interpolation method for efficient approximation of parametrized functions, Int. J. Numer. Methods Eng. 73, 521–543 (2008) [CrossRef] [Google Scholar]
  18. S. Chaturantabut, D. Sorensen, Discrete empirical interpolation for nonlinear model reduction, in Proceedings of the 48th IEEE Conference on Decision and Control, 2009 held jointly with the 2009 28th Chinese Control Conference, CDC/CCC 2009, 2010, pp. 4316–4321 [Google Scholar]
  19. K. Carlberg, C. Farhat, J. Cortial, D. Amsallem, The GNAT method for nonlinear model reduction: effective implementation and application to computational fluid dynamics and turbulent flows, J. Comput. Phys. 242, 623–647 (2013) [CrossRef] [MathSciNet] [Google Scholar]
  20. C. Farhat, P. Avery, T. Chapman, J. Cortial, Dimensional reduction of nonlinear finite element dynamic models with finite rotations and energy-based mesh sampling and weighting for computational efficiency, Int. J. Numer. Methods Eng. 98, 625–662 (2014) [CrossRef] [Google Scholar]
  21. M. Yano, A.T. Patera, An LP empirical quadrature procedure for reduced basis treatment of parametrized nonlinear PDEs, Comput. Methods Appl. Mech. Eng. 344, 1104–1123 (2019) [CrossRef] [Google Scholar]
  22. F. Casenave, N. Akkari, F. Bordeu, C. Rey, D. Ryckelynck, A nonintrusive distributed reduced-order modeling framework for nonlinear structural mechanics – application to elastoviscoplastic computations, Int. J. Numer. Methods Eng. 121, 32–53 (2020) [CrossRef] [Google Scholar]
  23. N. Cagniart, Y. Maday, B. Stamm, Model order reduction for problems with large convection effects, in B. Chetvershkin, W. Fitzgibbon, Y. Kuznetsov, P. Neittaanmäki, J. Periaux, O. Pironneau (Eds.), Contributions to Partial Differential Equations and Applications. Computational Methods in Applied Sciences, 2019 [Google Scholar]
  24. A. Iollo, D. Lombardi, Advection modes by optimal mass transfer, Phys. Rev. E 89, 022923 (2014) [CrossRef] [PubMed] [Google Scholar]
  25. J. Reiss, P. Schulze, J. Sesterhenn, V. Mehrmann, The shifted proper orthogonal decomposition: a mode decomposition for multiple transport phenomena, SIAM J. Sci.Comput. 40, A1322–A1344 (2018) [CrossRef] [Google Scholar]
  26. R. Zimmermann, B. Peherstorfer, K. Willcox, Geometric subspace updates with applications to online adaptive nonlinear model reduction, SIAM J. Matrix Anal. Appl. 39, 234–261 (2018) [CrossRef] [MathSciNet] [Google Scholar]
  27. F. Casenave, N. Akkari, An error indicator-based adaptive reduced order model for nonlinear structural mechanics - application to high-pressure turbine blades, Math. Comput. Appl. 24, 41 (2019) [Google Scholar]
  28. W. He, P. Avery, C. Farhat, In-situ adaptive reduction of nonlinear multiscale structural dynamics models 121, 4971–4988 (2020) [Google Scholar]
  29. T. Kim, D.L. James, Skipping steps in deformable simulation with online model reduction, ACM Trans. Graph. 28, 1–9 (2009) [Google Scholar]
  30. M. Ohlberger, F. Schindler, Error control for the localized reduced basis multiscale method with adaptive on-line enrichment, SIAM J. Sci. Comput. 37, A2865–A2895 (2015) [CrossRef] [Google Scholar]
  31. D. Amsallem, J. Cortial, C. Farhat, Towards real-time computational-fluid-dynamics-based aeroelastic computations using a database of reduced-order information, AIAA J. 48, 2029–2037 (2010) [CrossRef] [Google Scholar]
  32. D. Amsallem, C. Farhat, Interpolation method for adapting reduced-order models and application to aeroelasticity, AIAA J. 46, 1803–1813 (2008) [CrossRef] [Google Scholar]
  33. D. Amsallem, C. Farhat, An online method for interpolating linear parametric reduced-order models, SIAM J. Sci. Comput. 33, 2169–2198 (2011) [CrossRef] [MathSciNet] [Google Scholar]
  34. Y. Choi, G. Boncoraglio, S. Anderson, D. Amsallem, C. Farhat, Gradient-based constrained optimization using a database of linear reduced-order models, J. Comput. Phys. 423, 109787 (2020) [CrossRef] [MathSciNet] [Google Scholar]
  35. T. Lieu, C. Farhat, Adaptation of POD-based aeroelastic ROMs for varying Mach number and angle of attack: Application to a complete F-16 configuration, AIAA Paper 2005-7666 (2005) [Google Scholar]
  36. T. Lieu, C. Farhat, Adaptation of aeroelastic reduced-order models and application to an F-16 configuration, AIAA J. 45, 1244–1257 (2007) [CrossRef] [Google Scholar]
  37. T. Lieu, C. Farhat, M. Lesoinne, POD-based aeroelastic analysis of a complete F-16 configuration: ROM adaptation and demonstration, AIAA Paper 2005-2295 (2005) [Google Scholar]
  38. T. Lieu, C. Farhat, M. Lesoinne, Reduced-order fluid/structure modeling of a complete aircraft configuration, Comput. Methods Appl. Mech. Eng. 195, 5730–5742 (2006) [CrossRef] [Google Scholar]
  39. T. Lieu, M. Lesoinne, Parameter adaptation of reduced order models for three-dimensional flutter analysis, AIAA Paper 2004-0888 (2004) [Google Scholar]
  40. R. Mosquera, A. El Hamidi, A. Hamdouni, A. Falaize, Generalization of the Neville-Aitken Interpolation Algorithm on Grassmann Manifolds : Applications to Reduced Order Model, Int. J. Nume. Methods Fluids 93, 2421–2442 (2021) [CrossRef] [Google Scholar]
  41. R. Mosquera, A. Hamdouni, A. El Hamidi, C. Allery, POD basis interpolation via Inverse Distance Weighting on Grassmann manifolds, Discr.Continu. Dyn. Syst. Ser. S 12, 1743–1759 (2018) [Google Scholar]
  42. D. Amsallem, R. Tezaur, C. Farhat, Real-time solution of linear computational problems using databases of parametric reduced-order models with arbitrary underlying meshes, J. Comput. Phys. 326, 373–397 (2016) [CrossRef] [MathSciNet] [Google Scholar]
  43. S. Kaulmann, B. Haasdonk, Online greedy reduced basis construction using dictionaries. VI International Conference on Adaptive Modeling and Simulation (ADMOS 2013), 2013, pp. 365–376 [Google Scholar]
  44. Y. Maday, B. Stamm, Locally adaptive greedy approximations for anisotropic parameter reduced basis spaces, SIAM J. Sci. Comput. 35, A2417–A2441 (2013) [CrossRef] [Google Scholar]
  45. M. Dihlmann, M. Drohmann, B. Haasdonk, Model reduction of parametrized evolution problems using the reduced basis method with adaptive time partitioning, Proc. of ADMOS 2011, 64 (2011) [Google Scholar]
  46. M. Drohmann, B. Haasdonk, M. Ohlberger, Adaptive reduced basis methods for nonlinear convection-diffusion equations, Finite Volumes for Complex Applications VI Problems & Perspectives, Springer, 2011, pp. 369–377 [CrossRef] [Google Scholar]
  47. J. Eftang, A. Patera, E. Ronquist, An “hp” certified reduced basis method for parametrized elliptic partial differential equations, SIAM J. Sci. Comput. 32, 3170–3200 (2010) [CrossRef] [MathSciNet] [Google Scholar]
  48. B. Haasdonk, M. Dihlmann, M. Ohlberger, A training set and multiple bases generation approach for parametrized model reduction based on adaptive grids in parameter space, Math. Comput. Model. Dyn. Syst. 17, 423–442 (2011) [CrossRef] [MathSciNet] [Google Scholar]
  49. M.G. Kapteyn, D.J. Knezevic, K.E. Willcox, Toward predictive digital twins via component-based reduced-order models and interpretable machine learning. AIAA Scitech 2020 Forum, 2020, 0418 [Google Scholar]
  50. M.G. Kapteyn, K.E. Willcox, From physics-based models to predictive digital twins via interpretable machine learning. arXiv preprint, arXiv:2004.11356 (2020) [Google Scholar]
  51. B. Peherstorfer, D. Butnaru, K. Willcox, H.J. Bungartz, Localized discrete empirical interpolation method, SIAM J. Sci. Comput. 36, A168–A192 (2014) [CrossRef] [Google Scholar]
  52. D. Amsallem, M. Zahr, C. Farhat, Nonlinear model order reduction based on local reduced-order bases, Int. J. Numer. Methods Eng. 92, 891–916 (2012) [CrossRef] [Google Scholar]
  53. D. Amsallem, M. Zahr, K. Washabaugh, Fast local reduced basis updates for the efficient reduction of nonlinear systems with hyper-reduction, Adv. Comput. Math. 41, 1187–1230 (2015) [CrossRef] [MathSciNet] [Google Scholar]
  54. S. Grimberg, C. Farhat, R. Tezaur, C. Bou-Mosleh, Mesh sampling and weighting for the hyperreduction of nonlinear Petrov-Galerkin reduced-order models with local reduced-order bases, Int. J. Numer. Methods Eng. 122, 1846–1874 (2021) [CrossRef] [Google Scholar]
  55. F. Nguyen, S.M. Barhli, D.P. Muñoz, D. Ryckelynck, Computer vision with error estimation for reduced order modeling of macroscopic mechanical tests, Complexity 2018, 3791543 (2018) [Google Scholar]
  56. D. Ryckelynck, T. Goessel, F. Nguyen, Mechanical dissimilarity of defects in welded joints via Grassmann manifold and machine learning, Compt. Rend. Mécaniq. 348, 911–935 (2020) [Google Scholar]
  57. K. Washabaugh, D. Amsallem, M. Zahr, C. Farhat, Nonlinear model reduction for CFD problems using local reduced order bases, 42nd AIAA Fluid Dynamics Conference International Journal for Numerical Methods in Fluids, 2012, 2686 [Google Scholar]
  58. Y. Kim, Y. Choi, D. Widemann, T. Zohdi, A fast and accurate physics-informed neural network reduced order model with shallow masked autoencoder, Elsevier, 2021, u7 110841 [Google Scholar]
  59. K. Lee, K. Carlberg, Deep conservation: a latent-dynamics model for exact satisfaction of physical conservation laws, arXiv preprint, arXiv:1909.09754 (2019) [Google Scholar]
  60. K. Lee, K.T. Carlberg, Model reduction of dynamical systems on nonlinear manifolds using deep convolutional autoencoders, J. Comput. Phys. 404, 108973 (2020) [CrossRef] [MathSciNet] [Google Scholar]
  61. T. Daniel, F. Casenave, N. Akkari, A. Ketata, D. Ryckelynck, Physics-informed cluster analysis and a priori efficiency criterion for the construction of local reduced-order bases, 2021, https://arxiv.org/pdf/2103.13683.pdf [Google Scholar]
  62. T. Daniel, F. Casenave, N. Akkari, D. Ryckelynck, Data augmentation and feature selection for automatic model recommendation in computational physics, Math. Comput. Appl. 26, 17 (2021) [Google Scholar]
  63. H. Peng, F. Long, C. Ding, Feature selection based on mutual information: criteria of max-dependency, max-relevance, and min-redundancy, IEEE Trans. Pattern Anal. Mach. Intell. 27, 1226–1238 (2005) [CrossRef] [PubMed] [Google Scholar]
  64. B.N. Khoromskij, A. Litvinenko, H.G. Matthies, Application of hierarchical matrices for computing the Karhunen-Loève expansion, Computing 84, 49–67 (2009) [CrossRef] [MathSciNet] [Google Scholar]
  65. C. Scarth et al., Random field simulation over curved surfaces: applications to computational structural mechanics, Comput. Methods Appl. Mech. Eng. 345, 283–301 (2019) [CrossRef] [Google Scholar]
  66. V. Surazhsky, T. Surazhsky, D. Kirsanov, S.J. Gortler, H. Hoppe, Fast exact and approximate geodesics on meshes, ACM Trans. Graph. 24, 553–560 (2005) [CrossRef] [Google Scholar]
  67. J.S.B. Mitchell, D.M. Mount, C.H. Papadimitriou, The discrete geodesic problem, SIAM J, Comput. 16, 647–668 (1987) [CrossRef] [MathSciNet] [Google Scholar]
  68. V.R. Joseph, E. Gul, S. Ba, Maximum projection designs for computer experiments, Biometrika 102, 371–380 (2015) [CrossRef] [MathSciNet] [Google Scholar]
  69. R. Everson, L. Sirovich, Karhunen–Loève procedure for gappy data, J. Opt. Soc. Am. A 12, 1657–1664 (1995) [CrossRef] [Google Scholar]
  70. L. Kaufman, P.J.R. Leonard Kaufman, P.J. Rousseeuw, Finding Groups in Data: An Introduction to Cluster Analysis, Wiley, 1990 [Google Scholar]
  71. L. Kaufmann, P. Rousseeuw, Clustering by means of medoids, Data Analysis based on the L1-Norm and Related Methods, 1987, pp. 405–416 [Google Scholar]
  72. I. Borg, P. Groenen, Modern Multidimensional Scaling: Theory and Applications, 2nd edn. Springer Science & Business Media, 2005 [Google Scholar]
  73. J. de Leeuw Applications of convex analysis to multidimensional scaling, in J.R. Barra, F. Brodeau, G. Romier, B. van Cutsem (Eds.), Recent Developments in Statistics, 1977, pp. 133–145 [Google Scholar]
  74. F. Pedregosa, G. Varoquaux, A. Gramfort, et al. Scikit-learn: machine learning in Python, J. Mach. Learn. Res. 12, 2825–2830 (2011) [MathSciNet] [Google Scholar]
  75. J. Berkson, Application of the logistic function to bio-assay, J. Am. Stat. Assoc. 39, 357–365 (1944) [Google Scholar]
  76. D.R. Cox, The regression analysis of binary sequences, J. Royal Stat. Soc. Ser. B 20, 215–242 (1958) [Google Scholar]
  77. D.R. Cox, Some procedures connected with the logistic qualitative response curve, Research papers in probability and statistics, 1996, pp. 55–71 [Google Scholar]
  78. H. Zou, T. Hastie, Regularization and variable selection via the elastic net, J. Royal Stat. Soc. Ser B 67, 301–320 (2005) [CrossRef] [MathSciNet] [Google Scholar]
  79. T. Hastie, R. Tibshirani, J.H. Friedman, The Elements of Statistical Learning: Data Mining, Inference, and Prediction, Springer series in statistics 2nd edn, Springer, 2009 [Google Scholar]
  80. D. Ryckelynck, D.M. Benziane, A. Musienko, G. Cailletaud, Toward “green” mechanical simulations in materials science, Eur. J. Comput. Mech. 19, 365–388 (2010) [CrossRef] [Google Scholar]
  81. R. J. Asaro, Crystal plasticity, J. Appl. Mech. 50, 921–934 (1983) [CrossRef] [Google Scholar]
  82. L. Meric, P. Poubanne, G. Cailletaud, Single crystal modeling for structural calculations: Part 1 - model presentation, J. Eng. Mater. Technol. 113, 162–170 (1991) [CrossRef] [Google Scholar]
  83. F. Gallerneau, Etude et modelisation de l’endommagement d’un superalliage monocristallin revetu pour aube de turbine, PhD thesis - Mines ParisTech, 1995 [Google Scholar]
  84. F. Gallerneau, D. Nouailhas, J.L. Chaboche, A fatigue damage model including interaction effects with oxidation and creep damages, in Fatigue’ 96: Proceedings of the Sixth International Fatigue Congress, 1996, pp. 861–866 [Google Scholar]
5

Linear defects in the crystal structure.

Cite this article as: T. Daniel, F. Casenave, N. Akkari, D. Ryckelynck, C. Rey, Uncertainty quantification for industrial numerical simulation using dictionaries of reduced order models, Mechanics & Industry 23, 3 (2022)

All Tables

Table 1

Classification results.

In the text
Table 2

Widths of the confidence intervals (CI) for the expectations, expressed as percentages of the estimated expectations.

In the text
Table 3

Error indicators for the evaluation of the ROM-net on 20 new thermal loadings.

In the text

All Figures

thumbnail Fig. 1

Exploitation phase of a dictionary-based ROM-net. K local ROMs, combined with a classifier for automatic ROM recommendation, are used to predict the quantity of interest Z(x).
In the text
thumbnail Fig. 2

High-pressure turbine blade geometry and mesh (micro-perforations are not modeled).
In the text
thumbnail Fig. 3

Function ω(t) defining onecycle for the rotation speed.
In the text
thumbnail Fig. 4

Reference temperature field (on the left), temperature perturbation at the trailing edge (field 0 = δT0), and fluctuation modes (fields 1 to 4). The fluctuations in the fourth mode are located inside the blade, in the cooling channels.
In the text
thumbnail Fig. 5

Visualization of the MaxProj LHS DoE. The marginal distributions are represented on the diagonal. The 5D DoE is projected on 2D subspaces for visualization purposes, in order to check space-filling properties in 2D.
In the text
thumbnail Fig. 6

Visualization of the Sobol’ DoE. The marginal distributions are represented on the diagonal. The 5D DoE is projected on 2D subspaces for visualization purposes, in order to check space-filling properties in 2D.
In the text
thumbnail Fig. 7

MDS representation of the clustering results using the ROM-oriented dissimilarity measure on the quantity of interest (goal-oriented variant). On the left, the colors correspond to the expected clusters. On the right, the colors correspond to the clusters identified by the clustering algorithm. The positions of the labels 0 and 1 coincide with the positions of the clusters’ medoids. The MDS relative error ς(Z0 ; δ)∕ς(0;δ) is 12%.
In the text
thumbnail Fig. 8

The 3 fields on the left correspond to the medoid of cluster 0, and those onthe right correspond to the medoid of cluster 1. The fields in the first and the third columns show the differences between the medoids’ temperature fields and the reference temperature field Tref (the scale istruncated for the first field). The second and the fourth columns show the displacement magnitude field (top) and the quantity of interest (bottom).
In the text
thumbnail Fig. 9

MDSrepresentation of the clustering results. Orange points represent the snapshots selected for cluster 0, while the light blue points represent the snapshots selected for cluster 1. For each cluster, the snapshots are selected by a maximin procedure starting from the medoid.
In the text
thumbnail Fig. 10

Feature selection results. The kriging metamodel for redundancy terms is represented by the red curve and built from 800 true redundancy terms (blue points). The elements containing the selected nodes are represented in the turbine blade geometry.
In the text
thumbnail Fig. 11

Mean over the complete mesh of dual quantities of interest: accumulated plastic strain (left) and component 33 of the stress tensor σ33 (right), plotted as points where the x-coordinate is the reference value, and the y-coordinate is the considered reduced prediction.
In the text
thumbnail Fig. 12

Histograms and probability density functions of the quantities of interest (left) and (right).
In the text
thumbnail Fig. 13

Workflow for the ROM-net methology applied to the considered industrial setting in Sections 4.1-4.5.
In the text
thumbnail Fig. 14

Comparison between high-fidelity predictions (middle column) and ROM-net’s predictions (right-hand column). The field on the left represents the difference between the current temperature field (belonging to cluster 0) and the reference one. The other fields correspond to the increments of the quantity of interest with respectto its reference state obtained with the reference temperature field.
In the text
thumbnail Fig. 15

Comparison between high-fidelity predictions (middle column) and ROM-net’s predictions (right-hand column). The field on the left represents the difference between the current temperature field (belonging to cluster 1) and the reference one. The other fields correspond to the increments of the quantity of interest with respectto its reference state obtained with the reference temperature field.
In the text
thumbnail Fig. 16

Errors on the quantity of interest . The red (resp. blue) color is used for zones where the quantity of interest is overestimated (resp. underestimated).
In the text
thumbnail Fig. A.1

On the left: von Mises stress field σeq obtained when the rotation speed reaches its maximum value. On the right: accumulated plastic strain in octahedral slip systems at the end of the first cycle. Note: the foot of the high-pressure turbine blade has an elastic behavior, while the rest of the blade has a viscoplastic behavior described by a crystal plasticity model.
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.