Numerical simulation of transonic compressor under circumferential inlet distortion and rotor/stator interference using harmonic balance method

2018 ◽  
Vol 32 (12n13) ◽  
pp. 1840021
Author(s):  
Ziwei Wang ◽  
Xiong Jiang ◽  
Ti Chen ◽  
Yan Hao ◽  
Min Qiu

Simulating the unsteady flow of compressor under circumferential inlet distortion and rotor/stator interference would need full-annulus grid with a dual time method. This process is time consuming and needs a large amount of computational resources. Harmonic balance method simulates the unsteady flow in compressor on single passage grid with a series of steady simulations. This will largely increase the computational efficiency in comparison with the dual time method. However, most simulations with harmonic balance method are conducted on the flow under either circumferential inlet distortion or rotor/stator interference. Based on an in-house CFD code, the harmonic balance method is applied in the simulation of flow in the NASA Stage 35 under both circumferential inlet distortion and rotor/stator interference. As the unsteady flow is influenced by two different unsteady disturbances, it leads to the computational instability. The instability can be avoided by coupling the harmonic balance method with an optimizing algorithm. The computational result of harmonic balance method is compared with the result of full-annulus simulation. It denotes that, the harmonic balance method simulates the flow under circumferential inlet distortion and rotor/stator interference as precise as the full-annulus simulation with a speed-up of about 8 times.

Author(s):  
Javier Crespo ◽  
Roque Corral ◽  
Jesus Pueblas

An implicit harmonic balance method for modeling the unsteady non-linear periodic flow about vibrating airfoils in turbomachinery is presented. As departing point, an implicit edge-based three-dimensional Reynolds Averaged Navier-Stokes equations solver for unstructured grids that runs both on central processing units (CPUs) and graphics processing units (GPUs) is used. The harmonic balance method performs a spectral discretization of the time derivatives and marches in pseudo-time a new system of equations where the unknowns are the variables at different time samples. The application of the method to vibrating airfoils is discussed. It is shown that a time spectral scheme may achieve the same temporal accuracy at a much lower computational cost than a Backward Finite Difference method at the expense of using more memory. The performance of the implicit solver has been assessed with several application examples. A speed-up factor of 10 is obtained between the spectral and finite difference version of the code whereas and an additional speed-up factor of 10 is obtained when the code is ported to GPUs, totalizing a speed factor of 100. The performance of the solver in GPUs has been assessed using the 10th standard aeroelastic configuration and a transonic compressor.


2020 ◽  
Vol 32 (12) ◽  
pp. 126103
Author(s):  
Di Zhou ◽  
Zhiliang Lu ◽  
Tongqing Guo ◽  
Guoping Chen

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.


AIAA Journal ◽  
2012 ◽  
Vol 50 (3) ◽  
pp. 741-745 ◽  
Author(s):  
Dong-Kyun Im ◽  
Jang Hyuk Kwon ◽  
Soo Hyung Park

Author(s):  
Chad H. Custer ◽  
Jonathan M. Weiss ◽  
Venkataramanan Subramanian ◽  
William S. Clark ◽  
Kenneth C. Hall

The harmonic balance method implemented within STAR-CCM+ is a mixed frequency/time domain computational fluid dynamic technique, which enables the efficient calculation of time-periodic flows. The unsteady solution is stored at a small number of fixed time levels over one temporal period of the unsteady flow in a single blade passage in each blade row; thus the solution is periodic by construction. The individual time levels are coupled to one another through a spectral operator representing the time derivative term in the Navier-Stokes equation, and at the boundaries of the computational domain through the application of periodic and nonreflecting boundary conditions. The blade rows are connected to one another via a small number of fluid dynamic spinning modes characterized by nodal diameter and frequency. This periodic solution is driven to the correct solution using conventional (steady) CFD acceleration techniques, and thus is computationally efficient. Upon convergence, the time level solutions are Fourier transformed to obtain spatially varying Fourier coefficients of the flow variables. We find that a small number of time levels (or, equivalently, Fourier coefficients) are adequate to model even strongly nonlinear flows. Consequently, the method provides an unsteady solution at a computational cost significantly lower than traditional unsteady time marching methods. The implementation of this nonlinear harmonic balance method within STAR-CCM+ allows for the simulation of multiple blade rows. This capability is demonstrated and validated using a 1.5 stage cold flow axial turbine developed by the University of Aachen. Results produced using the harmonic balance method are compared to conventional time domain simulations using STAR-CCM+, and are also compared to published experimental data. It is shown that the harmonic balance method is able to accurately model the unsteady flow structures at a computational cost significantly lower than unsteady time domain simulation.


2015 ◽  
Vol 138 (3) ◽  
Author(s):  
Javier Crespo ◽  
Roque Corral ◽  
Jesus Pueblas

An implicit harmonic balance (HB) method for modeling the unsteady nonlinear periodic flow about vibrating airfoils in turbomachinery is presented. An implicit edge-based three-dimensional Reynolds-averaged Navier–Stokes equations (RANS) solver for unstructured grids, which runs both on central processing units (CPUs) and graphics processing units (GPUs), is used. The HB method performs a spectral discretization of the time derivatives and marches in pseudotime, a new system of equations where the unknowns are the variables at different time samples. The application of the method to vibrating airfoils is discussed. It is shown that a time-spectral scheme may achieve the same temporal accuracy at a much lower computational cost than a backward finite-difference method at the expense of using more memory. The performance of the implicit solver has been assessed with several application examples. A speed-up factor of 10 is obtained between the spectral and finite-difference version of the code, whereas an additional speed-up factor of 10 is obtained when the code is ported to GPUs, totalizing a speed factor of 100. The performance of the solver in GPUs has been assessed using the tenth standard aeroelastic configuration and a transonic compressor.


Sign in / Sign up

Export Citation Format

Share Document