Publications

  1. ResearchGate
    Hypervolume Scalarization for Shape Optimization to Improve Reliability and Cost of Ceramic Components

    Authors: Camilla Hahn, Kathrin Klamroth, Johanna Schultes and Michael Stiglmayr
    Abstract: In engineering applications one often has to trade-off among several objectives as, for example, the mechanical stability of a component, its efficiency, its weight and its cost. We consider a biobjective shape optimization problem maximizing the mechanical stability of a ceramic component under tensile load while minimizing its volume. Stability is thereby modeled using a Weibull-type formulation of the probability of failure under external loads. The PDE formulation of the mechanical state equation is discretized by a finite element method on a regular grid. To solve the discretized biobjective shape optimization problem we suggest a hypervolume scalarization, with which also unsupported efficient read more

    solutions can be determined without adding constraints to the problem formulation. We investigate the relation of the hypervolume scalarization to the weighted sum scalarization and to direct multiobjective descent methods. Since gradient information can be efficiently obtained by solving the adjoint equation, the scalarized problem can be solved by a gradient ascent algorithm. We evaluate our approach on a 2 D test case representing a straight joint under tensile load.

  2. arXiv:2004.10549
    An Analytical Study in Multi Physics and Multi Criteria Shape Optimization

    Authors: Hanno Gottschalk and Marco Reese
    Abstract: A simple multi-physical system for the potential flow of a fluid through a shroud in which a mechanical component, like a turbine vane, is placed, is modeled mathematically. We then consider a multi criteria shape optimization problem, when the shape of the component is allowed to vary under a certain set of 2nd order Hölder continuous differentiable transformations of a baseline shape with boundary of the same continuity class. As objective functions, we consider a simple loss model for the fluid dynamical efficiency and the probability of failure of the component due to repeated application of loads that stem from the fluid’s static pressure. read more

    For this multi-physical system, it is shown that under certain conditions the Pareto front is maximal in the sense that the Pareto front of the feasible set coincides with Pareto front of its closure. We also show that the set of all optimal forms with respect to scalarization techniques deforms continuously (in the Hausdorff metric) with respect to preference parameters.

  3. arXiv:2004.10820
    Tracing Locally Pareto Optimal Points by Numerical Integration

    Authors: Matthias Bolten, Onur Tanil Doganay, Hanno Gottschalk and Kathrin Klamroth
    Abstract: We suggest a novel approach for the efficient and reliable approximation of the Pareto front of sufficiently smooth unconstrained bi-criteria optimization problems. Optimality conditions formulated for weighted sum scalarizations of the problem yield a description of (parts of) the Pareto front as a parametric curve, parameterized by the scalarization parameter (i.e., the weight in the weighted sum scalarization). Its sensitivity w.r.t. parameter variations can be described by an ordinary differential equation (ODE). Starting from an arbitrary initial Pareto optimal solution, the Pareto front can then be traced by numerical integration. read more

    We provide an error analysis based on Lipschitz properties and suggest an explicit Runge-Kutta method for the numerical solution of the ODE. The method is validated on bi-criteria convex quadratic programming problems for which the exact solution is explicitly known, and numerically tested on complex bi-criteria shape optimization problems involving finite element discretizations of the state equation

  4. ResearchGate
    Towards Multidisciplinary Turbine Blade Tolerance Design Assessment using Adjoint Methods

    Authors: Matthias Bolten, Onur Tanil Doganay, Hanno Gottschalk and Kathrin Klamroth
    Abstract: This paper shows how to use discrete CFD and FEM ad-joint surface sensitivities to derive objective-based tolerances for turbine blades, instead of relying on geometric tolerances. For this purpose a multidisciplinary adjoint evaluation tool chain is introduced to quantify the effect of real manufacturing imperfections on aerodynamic efficiency and probabilistic low cycle fatigue life time. Before the adjoint method is applied, a numerical validation of the CFD and FEM adjoint gradients is performed using 102 heavy duty turbine vane scans. The results show that the absolute error for adjoint CFD gradients is below 0.5%, read more

    while the FEM life time gradient absolute errors are below 5%. The adjoint assessment tool chain further reduces the computational cost by around 85% for the investigated test case compared to non-linear methods. Through the application of the presented tool chain, the definition of specified objective-based tolerances becomes available as a design assessment tool and allows to improve overall turbine efficiency and the accuracy of life time prediction.

  5. ResearchGate
    GivEn – Shape Optimization for Gas Turbines in Volatile Energy Networks (Preprint)

    Authors: Backhaus, J. and Bolten, M. and Doganay, O.T., and Ehrhardt, M. and Engel, B. and Frey, Ch. and Gottschalk, H. and Günther, M. and Hahn, C., Jaschke, J., and Jaksch, P. and Klamroth, K. and Liefke, A. and Luft, D. Mäde, L. and Marciniak, V. and Reese, M. and Schultes, J. and Schulz, V. and Schmitz, S. and Steiner, J. and Stiglmayr, M.
    Abstract: This paper describes the project GivEn that develops a novel multicrite-ria optimization process for gas turbine blades and vanes using modern ”adjoint”shape optimization algorithms. read more

    Given the many start and shut-down processes ofgas power plants in volatile energy grids, besides optimizing gas turbine geometriesfor efficiency, the durability understood as minimization of the probability of fail-ure is a design objective of increasing importance. We also describe the underlyingcoupling structure of the multiphysical simulations and use modern, gradient based multicriteria optimization procedures to enhance the exploration of Pareto-optimalsolutions.

  6. mdpi.com
    Probabilistic Modeling of Slip System-Based Shear Stresses and Fatigue Behavior of Coarse-Grained Ni-Base Superalloy Considering Local Grain Anisotropy and Grain Orientation

    Authors: Benedikt Engel 1, Lucas Mäde, Philipp Lion, Nadine Moch, Hanno Gottschalk and Tilmann Beck
    Abstract: New probabilistic lifetime approaches for coarse grained Ni-base superalloys supplement current deterministic gas turbine component design philosophies; in order to reduce safety factors and push design limits. The models are based on statistical distributions of parameters, which determine the fatigue behavior under high temperature conditions. In the following paper, Low Cycle Fatigue (LCF) test data of several material batches of polycrystalline Ni-base superalloy René80 with different grain sizes and orientation distribution (random and textured) is presented and evaluated. read more

    The textured batch, i.e., with preferential grain orientation, showed higher LCF life. Three approaches to probabilistic crack initiation life modeling are presented. One is based on Weibull distributed crack initiation life while the other two approaches are based on probabilistic Schmid factors. In order to create a realistic Schmid factor distribution, polycrystalline finite element models of the specimens were generated using Voronoi tessellations and the local mechanical behavior investigated in dependence of different grain sizes and statistically distributed grain orientations. All models were first calibrated with test data of the material with random grain orientation and then used to predict the LCF life of the material with preferential grain orientation. By considering the local multiaxiality and resulting inhomogeneous shear stress distributions, as well as grain interaction through polycrystalline Finite Element Analysis (FEA) simulation, the best consistencies between predicted and observed crack initiation lives could be achieved.

  7. arXiv:1905.07566
    Gradient Based Biobjective Shape Optimization to Improve Reliability and Cost of Ceramic Components

    Authors: Onur T. Doganay, Camilla Hahn, Hanno Gottschalk, Kathrin Klamroth, Johanna Schultes, Michael Stiglmayr
    Abstract: We consider the simultaneous optimization of the reliability and the cost of a ceramic component in a biobjective PDE constrained shape optimization problem. A probabilistic Weibull-type model is used to assess the probability of failure of the component under tensile load, while the cost is assumed to be proportional to the volume of the component. Two different gradient-based optimization methods are suggested and compared at 2D test cases. read more

    The numerical implementation is based on a first discretize then optimize strategy and benefits from efficient gradient computations using adjoint equations. The resulting approximations of the Pareto front nicely exhibit the trade-off between reliability and cost and give rise to innovative shapes that compromise between these conflicting objectives.

  8. arXiv:1904.08650
    Efficient Techniques for Shape Optimization with Variational Inequalities using Adjoints

    Authors: Daniel Luft, Volker H. Schulz, Kathrin Welker
    Abstract: In general, standard necessary optimality conditions cannot be formulated in a straightforward manner for semi-smooth shape optimization problems. In this paper, we consider shape optimization problems constrained by variational inequalities of the first kind, so-called obstacle-type problems. Under appropriate assumptions, we prove existence of adjoints for regularized problems and convergence…read more

    to adjoints of the unregularized problem. Moreover, we derive shape derivatives for the regularized problem and prove convergence to a limit object. Based on this analysis, an efficient optimization algorithm is devised and tested numerically.

  9. arXiv:1705.05776
    Shape optimization to decrease failure probability

    Authors: Matthias Bolten, Hanno Gottschalk, Camilla Hahn, Mohamed Saadi
    Abstract: Ceramic is a material frequently used in industry because of its favorable properties. Common approaches in shape optimization for ceramic structures aim to minimize the tensile stress acting on the component, as it is the main driver for failure. In contrast to this, we follow a more natural approach by minimizing the component’s probability of failure under a given tensile load. Since the…read more

    fundamental work of Weibull, the probabilistic description of the strength of ceramics is standard and has been widely applied. Here, for the first time, the resulting failure probabilities are used as objective functions in PDE constrained shape optimization. To minimize the probability of failure, we choose a gradient based method combined with a first discretize then optimize approach. For discretization finite elements are used. Using the Lagrangian formalism, the shape gradient via the adjoint equation is calculated at low computational cost. The implementation is verified by comparison of it with a finite difference method applied to a minimal 2d example. Furthermore, we construct shape flows towards an optimal / improved shape in the case of a simple beam and a bended joint.

  10. arXiv:1806.04389
    Shape gradients for the failure probability of a mechanical component under cyclical loading

    Authors: Hanno Gottschalk, Mohamed Saadi
    Abstract: This work provides a numerical calculation of shape gradients of failure probabilities for mechanical components using a first discretize, then adjoint approach. While deterministic life prediction models for failure mechanisms are not (shape) differentiable, this changes in the case of probabilistic life prediction. The probabilistic, or reliability based, approach thus opens the way for…read more

    efficient adjoint methods in the design for mechanical integrity. In this work we propose, implement and validate a method for the numerical calculation of the shape gradients of failure probabilities for the failure mechanism low cycle fatigue (LCF), which applies to polycrystalline metal. Numerical examples range from a bended rod to a complex geometry from a turbo charger in 3D.

  11. arXiv:1803.01216
    Adjoint Method to Calculate Shape Gradients of Failure Probabilaties for Turbomachinery Components

    Authors: Hanno Gottschalk, Mohamed Saadi, Onur Tanil Doganay, Kathrin Klamroth, Sebastian Schmitz
    Abstract: In the optimization of turbomachinery components, shape sensitivities for fluid dynamical objective functions have been used for a long time. As peak stress is not a differential functional of the shape, such highly efficient procedures so far have been missing for objective functionals that stem from mechanical integrity. This changes, if deterministic lifing criteria are replaced by…read more

    probabilistic criteria, which have been introduced recently to the field of low cycle fatigue (LCF). Here we present a finite element (FEA) based first discretize, then adjoin approach to the calculation of shape gradients (sen- sitivities) for the failure probability with regard to probabilistic LCF and apply it to simple and complex geometries, as e.g. a blisk geometry. We review the computation of failure probabilities with a FEA postprocessor and sketch the computation of the relevant quantities for the adjoint method. We demonstrate high accuracy and computational efficiency of the adjoint method compared to finite difference schemes. We discuss implementation details for rotating components with cyclic boundary conditions. Finally, we shortly comment on future development steps and on potential applications in multi criteria optimization.

  12. arXiv:1702.05759
    Probabilistic LCF Risk Evaluation of a Turbine Vane by Combined Size Effect and Notch Support Modeling

    Authors: Lucas Mäde, Sebastian Schmitz, Hanno Gottschalk, Tilman Beck
    Abstract: A probabilistic risk assessment for low cycle fatigue (LCF) based on the so-called size effect has been applied on gas-turbine design in recent years. In contrast, notch support modeling for LCF which intends to consider the change in stress below the surface of critical LCF regions is known and applied for decades. Turbomachinery components often show sharp stress gradients and very localized…read more

    critical regions for LCF crack initiations so that a life prediction should also consider notch and size effects. The basic concept of a combined probabilistic model that includes both, size effect and notch support, is presented. In many cases it can improve LCF life predictions significantly, in particular compared to \textit{E-N} curve predictions of standard specimens where no notch support and size effect is considered. Here, an application of such a combined model is shown for a turbine vane.

  13. Conference: ECCM-ECFD 2018, At Glasgow
    Using Adjoint CFD to Quantify the Impact of Manufacturing Variations on a Heavy Duty Turbine Vane
    Authors: Alexander Liefke, Vincent Marciniak, Uwe Janoske, Hanno Gottschalk
    Abstract: Turbine efficiency is one of the main design criteria for heavy duty gas turbines. After the design, margin adaption factors are applied on the baseline to predict the impact of manufacturing variations (MV). These margins are normally based on testbed experience. A more detailed knowledge of the impact of MV, prior to testing, would therefore improve the margin prediction accuracy and could benefit in product cost and global efficiency. …read more


    For turbomachines the impact of MV can be quantified with a Monte Carlo (MC) simulation in combination with steady non-linear CFD calculations e.g. RANS. The drawback of this approach is the large number of RANS computations needed to quantify the impact of MV, which is prohibitive for a daily use in an industrial context. Assuming that the MV are small enough, the adjoint CFD method, which linearizes the governing equations, can be an alternative to the RANS evaluations. This kind of approach has been successfully used for compressors and turbines.

    The first part of this paper presents a systematic approach to evaluate a hand-derived and an algorithmic-derived version of the discrete adjoint CFD solver TRACE. To do so, the ERCOFTAC axial flow turbine known as Aachen Turbine has been selected. For the adjoint version comparison a NACA-like parametrization is applied to compare and validate the adjoint-generated with finite difference gradients.

    In the second part the adjoint-based method is applied to an industrial turbine vane to quantify the impact of MV. For this case real MV have been measured using 102 optical blade scans. The scans are used to generate the corresponding deformed geometries for which an adjoint and a RANS simulation are computed. The comparison between each computation demonstrates that the impact of realistic MV can be handled by the adjoint approach.