It is well known that fan/compressor blade flutter stability increases with the increase of reduced frequency.Less well-known is that the least stable inter blade phase angle(IBPA)increases with the drop of reduced fr...It is well known that fan/compressor blade flutter stability increases with the increase of reduced frequency.Less well-known is that the least stable inter blade phase angle(IBPA)increases with the drop of reduced frequency.However,it is quite striking that little can be found in the open literature about the mechanism to the observations.In this paper,a numerical investigation is carried out to uncover the mechanism of the effect of reduced frequency on flutter stability and the least stable IBPA.The NASA rotor 67 has been used as the test vehicle with its first bending and torsion modes being considered.The time domain harmonic balance method together with the influence coefficient method is used to obtain the worksum-IBPA curves for all cases.It is found that:1)the deterioration of flutter stability with the decrease of reduced frequency is dictated by the dominant decrease of aerodamping due to a blade own vibration;2)the increase of the least stable IBPA with the decrease of reduced frequency arises largely from the increase of the least stable IBPA of the aerodamping from the nearest blade on a blade pressure side.展开更多
This paper presents a thorough study of the effect of the Constant Eddy Viscosity(CEV)assumption on the optimization of a discrete adjoint-based design optimization system.First,the algorithms of the adjoint methods w...This paper presents a thorough study of the effect of the Constant Eddy Viscosity(CEV)assumption on the optimization of a discrete adjoint-based design optimization system.First,the algorithms of the adjoint methods with and without the CEV assumption are presented,followed by a discussion of the two methods’solution stability.Second,the sensitivity accuracy,adjoint solution stability,and Root Mean Square(RMS)residual convergence rates at both design and offdesign operating points are compared between the CEV and full viscosity adjoint methods in detail.Finally,a multi-point steady aerodynamic and a multi-objective unsteady aerodynamic and aeroelastic coupled design optimizations are performed to study the impact of the CEV assumption on optimization.Two gradient-based optimizers,the Sequential Least-Square Quadratic Programming(SLSQP)method and Steepest Descent Method(SDM)are respectively used to draw a firm conclusion.The results from the transonic NASA Rotor 67 show that the CEV assumption can deteriorate RMS residual convergence rates and even lead to solution instability,especially at a near stall point.Compared with the steady cases,the effect of the CEV assumption on unsteady sensitivity accuracy is much stronger.Nevertheless,the CEV adjoint solver is still capable of achieving optimization goals to some extent,particularly if the flow under consideration is benign.展开更多
Manufactured blades are inevitably different from their design intent,which leads to a deviation of the performance from the intended value.To quantify the associated performance uncertainty,many approaches have been ...Manufactured blades are inevitably different from their design intent,which leads to a deviation of the performance from the intended value.To quantify the associated performance uncertainty,many approaches have been developed.The traditional Monte Carlo method based on a Computational Fluid Dynamics solver(MC-CFD)for a three-dimensional compressor is prohibitively expensive.Existing alternatives to the MC-CFD,such as surrogate models and secondorder derivatives based on the adjoint method,can greatly reduce the computational cost.Nevertheless,they will encounter’the curse of dimensionality’except for the linear model based on the adjoint gradient(called MC-adj-linear).However,the MC-adj-linear model neglects the nonlinearity of the performance function.In this work,an improved method is proposed to circumvent the lowaccuracy problem of the MC-adj-linear without incurring the high cost of other alternative models.The method is applied to the study of the aerodynamic performance of an annular transonic compressor cascade,subject to prescribed geometric variability with industrial relevance.It is found that the proposed method achieves a significant accuracy improvement over the MC-adj-linear with low computational cost,showing the great potential for fast uncertainty quantification.展开更多
A Reduced Order Model(ROM)based analysis method for turbomachinery cascade coupled mode flutter is presented in this paper.The unsteady aerodynamic model is established by a system identification technique combined wi...A Reduced Order Model(ROM)based analysis method for turbomachinery cascade coupled mode flutter is presented in this paper.The unsteady aerodynamic model is established by a system identification technique combined with a set of Aerodynamic Influence Coefficients(AIC).Subsequently,the aerodynamic model is encoded into the state space and then coupled with the structural dynamic equations,resulting in a ROM of the cascade aeroelasticity.The cascade flutter can be determined by solving the eigenvalues of the ROM.Bending-torsional coupled mode flutter analysis for the Standard Configuration Eleven(SC11)cascade is used to validate the proposed method.展开更多
This paper revisits the Space-Time Gradient(STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspecti...This paper revisits the Space-Time Gradient(STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspective. The STG method requires reordering of blade passages according to their relative clocking positions with respect to blades of an adjacent blade row. As the space-clocking is linked to an equivalent time-clocking, the passage reordering can be performed according to the alternative time-clocking. With the time-clocking perspective, unsteady flow solutions from different passages of the same blade row are mapped to flow solutions of the same passage at different time instants or phase angles. Accordingly, the time derivative of the unsteady flow equation is discretized in time directly, which is more natural than transforming the time derivative to a spatial one as with the original STG method. To improve the solution accuracy, a ninth order difference scheme has been investigated for discretizing the time derivative. To achieve a stable solution for the high order scheme, the implicit solution method of Lower-Upper Symmetric GaussSeidel/Gauss-Seidel(LU-SGS/GS) has been employed. The NASA Stage 35 and its blade-countreduced variant are used to demonstrate the validity of the time-clocking based passage reordering and the advantages of the high order difference scheme for the STG method. Results from an existing harmonic balance flow solver are also provided to contrast the two methods in terms of solution stability and computational cost.展开更多
基金National Science and Technology Major Project(2017-II-0009-0023)National Nature Science Foundation of China(51976172)。
文摘It is well known that fan/compressor blade flutter stability increases with the increase of reduced frequency.Less well-known is that the least stable inter blade phase angle(IBPA)increases with the drop of reduced frequency.However,it is quite striking that little can be found in the open literature about the mechanism to the observations.In this paper,a numerical investigation is carried out to uncover the mechanism of the effect of reduced frequency on flutter stability and the least stable IBPA.The NASA rotor 67 has been used as the test vehicle with its first bending and torsion modes being considered.The time domain harmonic balance method together with the influence coefficient method is used to obtain the worksum-IBPA curves for all cases.It is found that:1)the deterioration of flutter stability with the decrease of reduced frequency is dictated by the dominant decrease of aerodamping due to a blade own vibration;2)the increase of the least stable IBPA with the decrease of reduced frequency arises largely from the increase of the least stable IBPA of the aerodamping from the nearest blade on a blade pressure side.
基金supported by the National Science and Technology Major Project,China(No.2017-II-0009-0023)China’s 111 project(No.B17037)sponsored by Innovation Foundation for Doctor Dissertation of Northwestern Polytechnical University,China.
文摘This paper presents a thorough study of the effect of the Constant Eddy Viscosity(CEV)assumption on the optimization of a discrete adjoint-based design optimization system.First,the algorithms of the adjoint methods with and without the CEV assumption are presented,followed by a discussion of the two methods’solution stability.Second,the sensitivity accuracy,adjoint solution stability,and Root Mean Square(RMS)residual convergence rates at both design and offdesign operating points are compared between the CEV and full viscosity adjoint methods in detail.Finally,a multi-point steady aerodynamic and a multi-objective unsteady aerodynamic and aeroelastic coupled design optimizations are performed to study the impact of the CEV assumption on optimization.Two gradient-based optimizers,the Sequential Least-Square Quadratic Programming(SLSQP)method and Steepest Descent Method(SDM)are respectively used to draw a firm conclusion.The results from the transonic NASA Rotor 67 show that the CEV assumption can deteriorate RMS residual convergence rates and even lead to solution instability,especially at a near stall point.Compared with the steady cases,the effect of the CEV assumption on unsteady sensitivity accuracy is much stronger.Nevertheless,the CEV adjoint solver is still capable of achieving optimization goals to some extent,particularly if the flow under consideration is benign.
基金funded by the National Natural Science Foundation of China(No.52006177)National Science and Technology Major Project,China(No.2017-II-0009-0023)。
文摘Manufactured blades are inevitably different from their design intent,which leads to a deviation of the performance from the intended value.To quantify the associated performance uncertainty,many approaches have been developed.The traditional Monte Carlo method based on a Computational Fluid Dynamics solver(MC-CFD)for a three-dimensional compressor is prohibitively expensive.Existing alternatives to the MC-CFD,such as surrogate models and secondorder derivatives based on the adjoint method,can greatly reduce the computational cost.Nevertheless,they will encounter’the curse of dimensionality’except for the linear model based on the adjoint gradient(called MC-adj-linear).However,the MC-adj-linear model neglects the nonlinearity of the performance function.In this work,an improved method is proposed to circumvent the lowaccuracy problem of the MC-adj-linear without incurring the high cost of other alternative models.The method is applied to the study of the aerodynamic performance of an annular transonic compressor cascade,subject to prescribed geometric variability with industrial relevance.It is found that the proposed method achieves a significant accuracy improvement over the MC-adj-linear with low computational cost,showing the great potential for fast uncertainty quantification.
基金supported by the National Science and Technology Major Project, China (No. 2017-II-0009-0023)the Aeronautical Science Foundation of China(No. 2020Z039053004)the Fundamental Research Funds for the Central Universities, China (No. 3102019OQD701)
文摘A Reduced Order Model(ROM)based analysis method for turbomachinery cascade coupled mode flutter is presented in this paper.The unsteady aerodynamic model is established by a system identification technique combined with a set of Aerodynamic Influence Coefficients(AIC).Subsequently,the aerodynamic model is encoded into the state space and then coupled with the structural dynamic equations,resulting in a ROM of the cascade aeroelasticity.The cascade flutter can be determined by solving the eigenvalues of the ROM.Bending-torsional coupled mode flutter analysis for the Standard Configuration Eleven(SC11)cascade is used to validate the proposed method.
基金co-supported by the National Natural Science Foundation of China(No.51976172)the National Science and Technology Major Project of China(No.2017-Ⅱ-0009-0023)。
文摘This paper revisits the Space-Time Gradient(STG) method which was developed for efficient analysis of unsteady flows due to rotor–stator interaction and presents the method from an alternative time-clocking perspective. The STG method requires reordering of blade passages according to their relative clocking positions with respect to blades of an adjacent blade row. As the space-clocking is linked to an equivalent time-clocking, the passage reordering can be performed according to the alternative time-clocking. With the time-clocking perspective, unsteady flow solutions from different passages of the same blade row are mapped to flow solutions of the same passage at different time instants or phase angles. Accordingly, the time derivative of the unsteady flow equation is discretized in time directly, which is more natural than transforming the time derivative to a spatial one as with the original STG method. To improve the solution accuracy, a ninth order difference scheme has been investigated for discretizing the time derivative. To achieve a stable solution for the high order scheme, the implicit solution method of Lower-Upper Symmetric GaussSeidel/Gauss-Seidel(LU-SGS/GS) has been employed. The NASA Stage 35 and its blade-countreduced variant are used to demonstrate the validity of the time-clocking based passage reordering and the advantages of the high order difference scheme for the STG method. Results from an existing harmonic balance flow solver are also provided to contrast the two methods in terms of solution stability and computational cost.