On the Extension of a Harmonic Balance Method for the Simulation of Compressor Casing Treatments

Author(s):  
Christian Voigt ◽  
Graham Ashcroft

In recent years both linear and nonlinear frequency domain methods have become increasingly popular in the simulation and investigation of time-periodic flows in turbomachinery. In this work the extension of an alternating frequency/time domain Harmonic Balance method to support arbitrary inter-domain block interfaces, with possibly different frames of reference, is described in detail. The approach outlined is based on the time-domain, area-based interpolation algorithm originally developed for the investigation of casing treatments. In this paper, it is shown that by solving the domain coupling problem in the time-domain it is possible to accurately and efficiently capture the flow physics of such complex, nonlinear problems as blade tip interaction with casing treatments in transonic compressors. To demonstrate and verify the basic algorithm the advection of a simple entropy disturbance in a subsonic duct flow is first computed. Secondly, unsteady flow due to rotor-stator interaction in a transonic compressor stage is simulated and the data compared with reference numerical methods. Finally, to validate the method a single stage transonic axial compressor with casing treatments is simulated and the results are compared with previously published time-domain simulations as well as experimental data based on particle image velocimetry measurements in the blade tip region.

Author(s):  
Laura Junge ◽  
Graham Ashcroft ◽  
Peter Jeschke ◽  
Christian Frey

Due to the relative motion between adjacent blade rows the aerodynamic flow fields within turbomachinery are normally dominated by deterministic, periodic phenomena. In the numerical simulation of such unsteady flows (nonlinear) frequency-domain methods are therefore attractive as they are capable of fully exploiting the given spatial and temporal periodicity, as well as capturing or modelling flow nonlinearity. Central to the efficiency and accuracy of such frequency-domain methods is the selection of the frequencies and the circumferential modes to be resolved in simulations. Whilst trivial in the context of the simulation of a single compressor- or turbine-stage, the choice of solution modes becomes substantially more involved in multi-stage configurations. In this work the importance of mode scattering, in the context of the unsteady aerodynamic field, is investigated and quantified. It is shown that scattered modes can substantially impact the unsteady flow field and are essential for the accurate modelling of wake propagation within multistage configurations. Furthermore, an iterative approach is outlined, based on the spectral analysis of the circumferential modes at the interfaces between blade rows, to identify the dominant solution modes that should be resolved in the adjacent blade row. To demonstrate the importance of mode scattering and validate the approach for their identification the unsteady blade row interaction within a 4.5 stage axial compressor is computed using both the harmonic balance method and, based on a full annulus midspan simulation, a time-domain method. Through the inclusion of scattered modes it is shown that the solution quality of the harmonic balance results is comparable to that of the nonlinear time-domain simulation.


1985 ◽  
Vol 52 (4) ◽  
pp. 958-964 ◽  
Author(s):  
C. Pierre ◽  
A. A. Ferri ◽  
E. H. Dowell

A multi-harmonic, frequency domain analysis of dry friction damped systems is presented which uses an incremental harmonic balance method. When compared with time domain solution methods, it is found that the incremental harmonic balance method can yield very accurate results with some advantages over the time domain methods. Both one and two degree-of-freedom systems are studied.


Author(s):  
Pengcheng Du ◽  
Fangfei Ning

Time periodic unsteady flows are often encountered in turbomachinery. Simulating such flows using conventional time marching approach is very time-consuming and hence expensive. To handle this problem, several Fourier-based reduced order models have been developed recently. Among these, the time-domain harmonic balance method solves the governing equations purely in the time domain and there is also no need for the turbulence model to be linearized, making it easy to be implemented in an existing RANS code. Thus, the time-domain harmonic balance method was chosen and incorporated into an in-house Navier-Stokes flow solver. Several test cases were performed for the validations of the developed code. They cover standard unsteady test cases such as the low speed vortex shedding cylinder flow and the Sajben transonic diffuser under periodically oscillating back pressure. Further, two different practical turbomachinery unsteady flows were considered. One is a transonic fan under circumferential inlet distortion and the other is the rotor-stator interactions in a single stage compressor. The results illustrate the capability of the harmonic balance method in capturing the dominant nonlinear effects. The number of harmonics should be retained in the harmonic balance method is depend on the strength of the nonlinear unsteady effects and differs from case to case. With appropriate number of harmonics retained, it can resolve the unsteady flow field satisfactory, meanwhile, reducing the computational time significantly. In a word, the harmonic balance method promise to be an effective way to simulate time periodic unsteady flows.


Author(s):  
Zhenxia Chai ◽  
Wei Liu ◽  
Xiaoliang Yang ◽  
Xu Liu

The harmonic balance method is an efficient frequency-domain approach for computing periodically unsteady flows. Applying harmonic balance method to the rapid calculation of forced periodic motions, an efficient prediction method for dynamic derivatives is established based on the Etkin unsteady aerodynamic model. The method is firstly validated by a standard model of the hypersonic missile hyper ballistic shape, and the harmonic balance method results show great consistency with both time-domain results and experimental data, consolidating the efficiency of our numerical algorithm. The method is subsequently applied for an air-breathing hypersonic vehicle to investigate the capability of the harmonic balance method in modeling unsteady flows around complex geometric configurations. Comparisons with the results of the time-domain method demonstrate that the harmonic balance method is capable of resolving such complicated flowfield with both accuracy and efficiency. To be specific, for the air-breathing hypersonic vehicle, the maximum difference of the dynamic derivatives calculated by the harmonic balance method and the time-domain method is only 6.56%, and the harmonic balance method shows at least an order of magnitude efficiency over the general time-domain method in the current study.


2014 ◽  
Vol 2014 ◽  
pp. 1-8
Author(s):  
Hongliang Yao ◽  
Qian Zhao ◽  
Qi Xu ◽  
Bangchun Wen

The efficiency and accuracy of common time and frequency domain methods that are used to simulate the response of a rotor system with malfunctions are compared and analyzed. The Newmark method and the incremental harmonic balance method are selected as typical representatives of time and frequency domain methods, respectively. To improve the simulation efficiency, the fixed interface component mode synthesis approach is combined with the Newmark method and the receptance approach is combined with the incremental harmonic balance method. Numerical simulations are performed for rotor systems with single and double frequency excitations. The inherent characteristic that determines the efficiency of the two methods is analyzed. The results of the analysis indicated that frequency domain methods are suitable single and double frequency excitation rotor systems, whereas time domain methods are more suitable for multifrequency excitation rotor systems.


Author(s):  
Christian Frey ◽  
Graham Ashcroft ◽  
Hans-Peter Kersken

This paper compares various approaches to simulate unsteady blade row interactions in turbomachinery. Unsteady simulations of turbomachinery flows have gained importance over the last years since increasing computing power allows the user to consider 3D unsteady flows for industrially relevant configurations. Furthermore, for turbomachinery flows, the last two decades have seen considerable efforts in developing adequate CFD methods which exploit the rotational symmetries of blade rows and are therefore up to several orders of magnitude more efficient than the standard unsteady approach for full wheel configurations. This paper focusses on the harmonic balance method which has been developed recently by the authors. The system of equations as well as the iterative solver are formulated in the frequency domain. The aim of this paper is to compare the harmonic balance method with the time-linearized as well as the non-linear unsteady approach. For the latter the unsteady flow fields in a fan stage are compared to reference results obtained with a highly resolved unsteady simulation. Moreover the amplitudes of the acoustic modes which are due to the rotor stator interaction are compared to measurement data available for this fan stage. The harmonic balance results for different sets of harmonics in the blade rows are used to explain the minor discrepancies between the time-linearized and unsteady results published by the authors in previous publications. The results show that the differences are primarily due to the neglection of the two-way coupling in the time-linearized simulations.


AIAA Journal ◽  
2014 ◽  
Vol 52 (1) ◽  
pp. 62-71 ◽  
Author(s):  
Frédéric Sicot ◽  
Adrien Gomar ◽  
Guillaume Dufour ◽  
Alain Dugeai

Author(s):  
Jong-Yun Yoon ◽  
Hwan-Sik Yoon

This paper presents the nonlinear frequency response of a multistage clutch damper system in the framework of the harmonic balance method. For the numerical analysis, a multistage clutch damper with multiple nonlinearities is modeled as a single degree-of-freedom torsional system subjected to sinusoidal excitations. The nonlinearities include piecewise-linear stiffness, hysteresis, and preload all with asymmetric transition angles. Then, the nonlinear frequency response of the system is numerically obtained by applying the Newton–Raphson method to a system equation formulated by using the harmonic balance method. The resulting nonlinear frequency response is then compared with that obtained by direct numerical simulation of the system in the time domain. Using the simulation results, the stability characteristics and existence of quasi-harmonic response of the system are investigated. Also, the effect of stiffness values on the dynamic performance of the system is examined.


Sign in / Sign up

Export Citation Format

Share Document