Over the past years, nonlinear frequency-domain methods have become a state-of-the-art technique for the numerical simulation of unsteady flow fields within multistage turbomachinery as they are capable of fully exploiting the given spatial and temporal periodicities, as well as modeling flow nonlinearities in a computationally efficient manner. Despite this success, it still remains a significant challenge to capture nonlinear interaction effects within the context of configurations with multiple fundamental frequencies. If all frequencies are integer multiples of a common fundamental frequency, the interval spanned by the sampling points typically resolves the period of the common base frequency. For configurations in which the common frequency is very low in relation to the frequencies of primary interest, many sampling points are required to resolve the highest harmonic of the common fundamental frequency and the method becomes inefficient. In addition, when a problem can no longer be described by harmonic perturbations that are integer multiples of one fundamental frequency, as it may occur in two-shaft configurations or when simulating the nonlinear interaction in the context of forced response or flutter, then the standard discrete Fourier transform is no longer suitable, and the basic harmonic balance (HB) method requires extension. One possible approach is to use almost periodic Fourier transforms with equidistant or nonequidistant time sampling. However, the definition of suitable sampling points that lead to well-conditioned Fourier transform matrices and small aliasing errors is an intricate issue and far from straightforward. To overcome the issues regarding multifrequency problems described earlier, a new harmonic balance approach based on multidimensional Fourier transforms in time is presented. The basic idea of the approach is that instead of defining common sampling points in a common time period, separate time domains, one for each base frequency, are spanned, and the sampling points are computed equidistantly within each base frequency's period. Since the sampling domain is now extended to a multidimensional time-domain, all time instant combinations covering the whole multidimensional domain are computed as the Cartesian product of the sampling points on the axes. In a similar fashion, the frequency-domain is extended to a multidimensional frequency-domain by the Cartesian product of the harmonics of each base frequency, so that every point defined by the Cartesian product is an integer linear combination of the occurring frequencies. In this way, the proposed method is capable of fully integrating the nonlinear coupling effects between higher harmonics of different fundamental frequencies by using multidimensional discrete Fourier transforms within the harmonic balance solution procedure. The aim of this paper is to introduce the multidimensional harmonic balance method in detail and demonstrate the capability of the approach to simultaneously capture unsteady disturbances with arbitrary excitation frequencies. Therefore, the well-established aeroelasticity testcase standard configuration 10 in the presence of an artificial inflow disturbance, which mimics an upstream blade wake, is investigated. The crucial aspect of the proposed testcase is that a small ratio of the frequency of the inflow disturbance and the blades vibration frequency is chosen. To demonstrate the advantages of the newly proposed multidimensional harmonic balance approach, the results are compared to unsteady simulations in the time-domain and to state-of-the-art frequency-domain methods based on one-dimensional discrete Fourier transforms.

Flutter stability is a dominant design constraint of modern turbines. Thus, flutter-tolerant designs are currently explored, where the resulting vibrations remain within acceptable bounds. In particular, friction damping has the potential to yield limit cycle oscillations (LCOs) in the presence of a flutter instability. To predict such LCOs, it is the current practice to model the aerodynamic forces in terms of aerodynamic influence coefficients for a linearized structural model with fixed oscillation frequency. This approach neglects that both the nonlinear contact interactions and the aerodynamic stiffness cause a change in the deflection shape and the frequency of the LCO. This, in turn, may have a substantial effect on the aerodynamic damping. The goal of this paper is to assess the importance of these neglected interactions. To this end, a state-of-the-art aero-elastic model of a low pressure turbine blade row is considered, undergoing nonlinear frictional contact interactions in the tip shroud interfaces. The LCOs are computed with a fully coupled harmonic balance method, which iteratively computes the Fourier coefficients of structural deformation and conservative flow variables, as well as the a priori unknown frequency. The coupled algorithm was found to provide excellent computational robustness and efficiency. Moreover, a refinement of the conventional energy method is developed and assessed, which accounts for both the nonlinear contact boundary conditions and the linearized aerodynamic influence. It is found that the conventional energy method may not predict a limit cycle oscillation at all while the novel approach presented here can.

The impact of the unsteadiness in the considered turbulence quantities on the numerical prediction of the aeroelastic behavior of a low-pressure turbine (LPT) rotor blade is evaluated by means of a numerical study. In this context, one of the main objectives of this work is to compare different nonlinear harmonic balance (HB) approaches—one neglecting and one considering the unsteadiness in the employed turbulence models—with a conventional nonlinear solver of the unsteady Reynolds-averaged Navier–Stokes (URANS) equations in the time domain. In order to avoid unphysical oscillations in the turbulence quantities caused by the Gibbs phenomenon in the chosen HB approach, a filter method based on the Lanczos filter is developed. The developed filter method is applied in the course of the HB simulations considering the unsteadiness in the underlying turbulence model. Furthermore, the impact of its application on the solution of the flow field and on the unsteady surface pressure of the rotor blade, in particular, is discussed in the context of this work.

This paper describes the development and initial application of an adjoint harmonic balance (HB) solver. The HB method is a numerical method formulated in the frequency domain which is particularly suitable for the simulation of periodic unsteady flow phenomena in turbomachinery. Successful applications of this method include unsteady aerodynamics as well as aeroacoustics and aeroelasticity. Here, we focus on forced response due to the interaction of neighboring blade rows. In the simulation-based design and optimization of turbomachinery components, it is often helpful to be able to compute not only the objective values—e.g., performance data of a component—themselves but also their sensitivities with respect to variations of the geometry. An efficient way to compute such sensitivities for a large number of geometric changes is the application of the adjoint method. While this is frequently used in the context of steady computational fluid dynamics (CFD), it becomes prohibitively expensive for unsteady simulations in the time domain. For unsteady methods in the frequency domain, the use of adjoint solvers is feasible but still challenging. The present approach employs the reverse mode of algorithmic differentiation (AD) to construct a discrete adjoint of an existing HB solver in the framework of an industrially applied CFD code. The paper discusses implementational issues as well as the performance of the adjoint solver, in particular regarding memory requirements. The presented method is applied to compute the sensitivities of aeroelastic objectives with respect to geometric changes in a turbine stage.