Development of a Harmonic Balance Method-Based Numerical Strategy for Blade-Tip/Casing Interactions: Application to NASA Rotor 37

Author(s):  
Yann Colaïtis ◽  
Alain Batailly

Abstract In this study, a frequency-domain approach based on the harmonic balance method coupled to a predictor-corrector continuation algorithm is implemented for the qualitative analysis of blade-tip/casing contacts in aircraft engines. Unilateral contact and dry friction are taken into account through a regularized penalty law. To enhance the robustness of the methodology, particular attention is paid to the mitigation of the Gibbs phenomenon. To this end, the employed Alternating Frequency/Time scheme features a Lanczos σ-approximation so that spurious oscillations of the computed nonlinear contact forces become negligible. This approach is applied in combination with a model reduction technique on an industrial compressor blade: NASA rotor 37. In order to assess the influence of both the contact law regularization and the Lanczos σ-approximation, obtained results are thoroughly compared to an existing time integration-based numerical strategy relying on a Lagrange multiplier-based approach for contact treatment and that was previously confronted to experimental results. Presented results underline the very good agreement between the proposed methodology and the reference time integration numerical strategy. The proposed developments thus complement existing results on blade-tip/casing contact adding a much needed qualitative understanding of the interaction and an accurate assessment of the contact stiffening phenomenon.

2021 ◽  
Author(s):  
Yann Colaïtis ◽  
Alain Batailly

Abstract In this study, a frequency-domain approach based on the harmonic balance method coupled to a predictor-corrector continuation algorithm is implemented for the qualitative analysis of blade-tip/casing contacts in aircraft engines. Unilateral contact and dry friction are taken into account through a regularized penalty law. To enhance the robustness of the methodology, particular attention is paid to the mitigation of the Gibbs phenomenon. To this end, the employed Alternating Frequency/Time scheme features a Lanczos σ-approximation so that spurious oscillations of the computed nonlinear contact forces become negligible. This approach is applied in combination with a model reduction technique on an industrial compressor blade: NASA rotor 37. In order to assess the influence of both the contact law regularization and the Lanczos σ-approximation, obtained results are thoroughly compared to an existing time integration-based numerical strategy relying on a Lagrange multiplier-based approach for contact treatment and that was previously confronted to experimental results. Presented results underline the very good agreement between the proposed methodology and the reference time integration numerical strategy. The proposed developments thus complement existing results on blade-tip/casing contact adding a much needed qualitative understanding of the interaction and an accurate assessment of the contact stiffening phenomenon.


Author(s):  
Loi¨c Salles ◽  
Laurent Blanc ◽  
Fabrice Thouverez ◽  
Alexander M. Gouskov ◽  
Pierrick Jean

Contact interfaces with dry friction are frequently used in turbomachinery. Dry friction damping produced by the sliding surfaces of these interfaces reduces the amplitude of bladed-disk vibration. The relative displacements at these interfaces lead to fretting-wear which reduces the average life expectancy of the structure. Frequency response functions are calculated numerically by using the multi-Harmonic Balance Method (mHBM). The Dynamic Lagrangian Frequency-Time method is used to calculate contact forces in the frequency domain. A new strategy for solving non-linear systems based on dual time stepping is applied. This method is faster than using Newton solvers. It was used successfully for solving Nonlinear CFD equations in the frequency domain. This new approach allows identifying the steady state of worn systems by integrating wear rate equations a on dual time scale. The dual time equations are integrated by an implicit scheme. Of the different orders tested, the first order scheme provided the best results.


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):  
Mayank Sharma ◽  
Nathan A. Wukie ◽  
Matteo Ugolotti ◽  
Mark G. Turner

The Harmonic Balance method is well suited for analyzing unsteadiness in turbomachinery flows comprised of a few dominant frequencies. A harmonic condition is imposed on the temporal derivatives through a Fourier transform operation. The solution is then reinterpreted as a time-domain problem, where several instances of time (lying within the largest period) are solved for simultaneously with the enforcement of the time-harmonic condition providing coupling between time instances. A discontinuous Galerkin discretization is used together with overset grids to provide higher-order spatial accuracy and flexibility in representing complex geometry. In this work, the discontinuous Galerkin infrastructure is extended for unsteady problems with a Harmonic Balance method and a Diagonally Implicit Runge-Kutta time-integrator. Verification results are presented for both time integration approaches in addition to results for a turbine blade with unsteadiness driven by a prescribed unsteady inlet boundary condition. Comparisons of results from the Harmonic Balance and Diagonally Implicit Runge-Kutta approaches are very close, with some small discrepancies that require further investigation. Significantly, rapid convergence from the Newton solver is obtained for the Harmonic Balance approach applied to the Euler equations for the turbine blade problem. Solutions converged by 8–10 orders of magnitude are obtained in between 5 and 16 Newton steps.


Author(s):  
Ender Cigeroglu ◽  
Ning An ◽  
Chia-Hsiang Menq

In this paper, a forced response prediction method for the analysis of constrained and unconstrained structures coupled through frictional contacts is presented. This type of frictional contact problem arises in vibration damping of turbine blades, in which dampers and blades constitute the unconstrained and constrained structures, respectively. The model of the unconstrained/free structure includes six rigid body modes and several elastic modes, the number of which depends on the excitation frequency. In other words, the motion of the free structure is not artificially constrained. When modeling the contact surfaces between the constrained and free structure, discrete contact points along with contact stiffnesses are distributed on the friction interfaces. At each contact point, contact stiffness is determined and employed in order to take into account the effects of higher frequency modes that are omitted in the dynamic analysis. Depending on the normal force acting on the contact interfaces, quasistatic contact analysis is initially employed to determine the contact area as well as the initial preload or gap at each contact point due to the normal load. A friction model is employed to determine the three-dimensional nonlinear contact forces, and the relationship between the contact forces and the relative motion is utilized by the harmonic balance method. As the relative motion is expressed as a modal superposition, the unknown variables, and thus the resulting nonlinear algebraic equations in the harmonic balance method, are in proportion to the number of modes employed. Therefore the number of contact points used is irrelevant. The developed method is applied to a bladed-disk system with wedge dampers where the dampers constitute the unconstrained structure, and the effects of normal load on the rigid body motion of the damper are investigated. It is shown that the effect of rotational motion is significant, particularly for the in-phase vibration modes. Moreover, the effect of partial slip in the forced response analysis and the effect of the number of harmonics employed by the harmonic balance method are examined. Finally, the prediction for a test case is compared with the test data to verify the developed method.


Author(s):  
Loïc Salles ◽  
Laurent Blanc ◽  
Fabrice Thouverez ◽  
Alexander M. Gouskov ◽  
Pierrick Jean

Contact interfaces with dry friction are frequently used in turbomachinery. Dry friction damping produced by the sliding surfaces of these interfaces reduces the amplitude of bladed-disk vibration. The relative displacements at these interfaces lead to fretting-wear which reduces the average life expectancy of the structure. Frequency response functions are calculated numerically by using the multi-harmonic balance method (mHBM). The dynamic Lagrangian frequency-time method is used to calculate contact forces in the frequency domain. A new strategy for solving nonlinear systems based on dual time stepping is applied. This method is faster than using Newton solvers. It was used successfully for solving Nonlinear CFD equations in the frequency domain. This new approach allows identifying the steady state of worn systems by integrating wear rate equations a on dual time scale. The dual time equations are integrated by an implicit scheme. Of the different orders tested, the first order scheme provided the best results.


Sign in / Sign up

Export Citation Format

Share Document