The possibility of taking into account unsteady flow effects if performing turbomachinery shape optimization is attractive to accurately address inherently time dependent design problems. The harmonic balance method is an efficient solution for computational fluid dynamics problems of turbomachinery characterized by quasi-periodic flows. If applied in combination with adjoint methods, it enables the possibility to deal with unsteady fluid-dynamic design in a cost effective manner, opening the way towards multi-disciplinary applications. This paper presents the development of a novel fully-turbulent discrete adjoint based on the time domain harmonic balance method and its application to the constrained fluid dynamic optimization of an axial turbine stage. As opposed to previous works, the proposed method does not require any assumption on the turbulent eddy viscosity and on the set of input frequencies. The results show that the method provides accurate gradients, if compared with second order finite differences, and significant deviation with respect to the sensitivity computed with the constant eddy viscosity approximation. The application of the method to the fluid-dynamic shape optimization of the exemplary stage leads to improve the total-to-static efficiency of 0.8%. The efficiency increase is found to be higher than that obtained by means of a steady state optimization method.

This paper documents a fully turbulent discrete ad joint method for three-dimensional multistage turbomachinery design. The method is based on a duality preserving algorithm and is implemented in the open-source computational fluid dynamics tool SU2. The SU2 Reynolds-averaged Navier–Stokes solver is first extended to treat three dimensional steady turbomachinery flow using a conservative formulation of the mixing-plane coupled to non reflective boundary conditions. The numerical features of the flow solver are automatically inherited by the discrete ad joint solver, ensuring the same convergence rate of the primal solver. The flow solver is then validated against experimental data available for three turbine configurations, namely, a one-and-half axial turbine stage, a transonic radial turbine coupled to a downstream diffuser, and a supersonic mini–organic Rankine cycle radial turbine operating with a fluid made by a heavy molecule. Finally the ad joint-based optimization framework is applied to the concurrent shape optimization of three rows of the axial turbine, demonstrating the advantages deriving from adopting multi row automated design methods in the context of turbomachinery design.

Fast and accurate computation of thermo-physical properties is essential in computationally expensive simulations involving fluid flows that significantly depart from the ideal gas or ideal liquid behavior. A look-up table algorithm based on unstructured grids is proposed and applied to non-ideal compressible fluid dynamics simulations. The algorithm grants the possibility of a fully automated generation of the tabulated thermodynamic region for any boundary and to use mesh refinement. Results show that the proposed algorithm leads to a computational cost reduction up to one order of magnitude, while retaining the same accuracy level compared to simulations based on more complex equation of state. Furthermore, a comparison of the LuT algorithm with a uniformly spaced quadrilateral tabulation method resulted in similar performance and accuracy.

The capabilities of the open-source SU2 software suite for the numerical simulation of viscous flows over unstructured grid are extended to non-ideal compressible-fluid dynamics (NICFD). A built-in thermodynamic library is incorporated to account for the non-ideal thermodynamic characteristics of fluid flows evolving in the close proximity of the liquid-vapour saturation curve and critical point. The numerical methods, namely the Approximate Riemann Solvers (ARS), viscous fluxes and boundary conditions are generalised to non-ideal fluid properties. Quantities of interest for turbomachinery cascades, as loss coefficients and flow angles, can be automatically determined and used for design optimization. A variety of test cases are carried out to assess the performance of the solver. At first, numerical methods are verified against analytical solution of reference NICFD test cases, including steady shock reflection and unsteady shock tube. Then, non-ideal gas effects in planar nozzles and past turbine cascades, typically encountered in Organic Rankine Cycle applications, are investigated and debated. The obtained results demonstrate that SU2 is highly suited for the analysis and the automatic design of internal flow devices operating in the non-ideal compressible-fluid regime.

High temperature Organic Rankine Cycles power systems of low power capacity, i.e. 3-50 kW_e, are receiving recognition for distributed and mobile energy generation applications. For this type of power plants, it is customary to adopt a radial-turbine as prime mover, essentially for their ability to cope with very large volumetric flow ratio with limited fluid-dynamic penalty. To date, the design of such turbines is based on design guidelines and loss models developed mainly for turbo-chargers, subsequently adapted by means of non-validated computational fluid-dynamic calculations. In the attempt to provide data sets for CFD validation and calibration of loss models, a mini-ORC radial inflow turbine delivering 10kW of mechanical power will be realized and tested in the Propulsion and Power Laboratory of TU-Delft. The fluid dynamic design and characterization of the machine is detailed in this paper. According to available models, the results indicate that the optimal layout of mini-ORC turbines can substantially differ from that of radial-inflow turbines utilized in more traditional applications, strengthening the need of experimental campaigns to support the conception of new design practices.