This study investigates the influence of mesh resolution and turbulence model selection on the accuracy of numerical simulations for transonic flow,with particular emphasis on shock-boundary layer interaction phenomen...This study investigates the influence of mesh resolution and turbulence model selection on the accuracy of numerical simulations for transonic flow,with particular emphasis on shock-boundary layer interaction phenomena.Accurate prediction of such flows is notoriously difficult due to the sensitivity to near-wall resolution,global mesh density,and turbulence model assumptions,and this problem motivates the present work.Two solvers were employed,rhoCentralFoam(unsteady)and TSLAeroFoam(steady-state),both are compressible and density-based and implemented within the OpenFOAM framework.The investigation focuses on three different non-dimensional wall distance(y+)values of 1,2.5 and 5,each implemented with both moderate and fine mesh resolutions.Three turbulence models—Spalart-Allmaras(SA),k-ωShear Stress Transport(SST),and k- Realizable—were evaluated at M=0.74,Re=2.7×106,andα=3.19°.Results showed that while both solvers achieved good overall agreement with experimental data,particularly in terms of pressure distribution,lift coefficient,and shock location,noticeable differences still emerged.The k-ωSST model consistently delivered the most robust performance across all cases,capturing the shock position on y+≈1 meshes with deviations below0.02c compared to the experiment,and maintaining accuracy even at y+≈5.The k- Realizable model was highly sensitive to near-wall resolution,displacing shocks downstream at higher y+values,whereas Spalart-Allmaras remained broadly comparable to the k-ωSST model in predictive performance.The rhoCentralFoam solver achieved consistently better lift predictions,staying within about 2%of the experimental value on average,whereas TSLAeroFoam overpredicted it by around 4%.For transonic Reynolds-Averaged Navier-Stokes(RANS)simulations,unsteady k-ωSST with y+≈1 is recommended for maximum fidelity,whereas steady k-ωSST or SA simulations offer a practical option for quick and reasonably accurate aerodynamic predictions.展开更多
The utilization of Inlet Guide Vane (IGV) plays a key factor in affecting the instability evolution. Existing literature mainly focuses on the effect of IGV on instability inception that occurs in the rotor region. Ho...The utilization of Inlet Guide Vane (IGV) plays a key factor in affecting the instability evolution. Existing literature mainly focuses on the effect of IGV on instability inception that occurs in the rotor region. However, with the emergence of compressor instability starting from the stator region, the mechanism of various instability inceptions that occurs in different blade rows due to the change of IGV angles should be further examined. In this study, experiments were focused on three types of instability inceptions observed previously in a 1.5-stage axial flow compressor. To analyze the conversion of stall evolutions, the compressor rotating speed was set to 17 160 r/min, at which both the blade loading in the stator hub region and rotor tip region were close to the critical value before final compressor stall. Meanwhile, the dynamic test points with high-response were placed to monitor the pressures both at the stator trailing edges and rotor tips. The results indicate that the variation of reaction determines the region where initial instability occurs. Indeed, negative pre-rotation of the inlet guide vane leads to high-reaction, initiating stall disturbance from the rotor region. Positive pre-rotation results in low-reaction, initiating stall disturbance from the stator region. Furthermore, the type of instability evolution is affected by the radial loading distribution under different IGV angles. Specifically, a spike-type inception occurs at the rotor blade tip with a large angle of attack at the rotor inlet (−2°, −4° and −6°). Meanwhile, the critical total pressure ratio at the rotor tip is 1.40 near stall. As the angle of attack decreases, the stator blade loading reaches its critical boundary, with a value of approximately 1.35. At this moment, if the rotor tip maintains high blade loading similar to the stator hub, the partial surge occurs (0° and +2°);otherwise, the hub instability occurs (+4° and +6°).展开更多
The influence of Impedance Boundary Condition (IBC) on transonic compressors is investigated. A systematic input–output analytical framework is developed, which treats the nonlinearities as unknown forcing terms. The...The influence of Impedance Boundary Condition (IBC) on transonic compressors is investigated. A systematic input–output analytical framework is developed, which treats the nonlinearities as unknown forcing terms. The framework is validated through the experiments of rotating inlet distortion within a low-speed compressor. The input–output method is subsequently applied to transonic compressors, including NASA Rotor37 and Stage35, wherein impedance optimization is studied along with the exploration of its fundamental mechanisms. The IBC is employed to model the effect of Casing Treatment (CT). The optimal complex impedance values are determined through predicted results and tested across a range of circumferential modes and forcing frequencies. The IBC significantly reduces the energy and Reynolds stress gain, notably at the first-order circumferential mode and within the Rotor Rotating Frequency (RRF) range. Output modes reveal that transonic compressors with fine-tuned impedance values exhibit a more confined perturbation distribution and redistribute the perturbations compared to the uncontrolled case. Additionally, the roles of resistance and reactance are elucidated through input–output analysis, and resistance determines the energy transfer direction between flow and pressure waves and modulates the amplitude, whereas reactance modifies the phase relationships and attenuates the perturbations.展开更多
In this paper,a series of flutter simulations are carried out to investigate the effects of split drag rudder(SDR)on the transonic flutter characteristic of rigid NACA 64A010.A structural dynamic model addressing two-...In this paper,a series of flutter simulations are carried out to investigate the effects of split drag rudder(SDR)on the transonic flutter characteristic of rigid NACA 64A010.A structural dynamic model addressing two-degree-of-freedom pitch-plunge aeroelastic oscillations was coupled with the unsteady Reynolds-averaged Navier-Stokes equations to perform flutter simulation.Meanwhile,the influence mechanism of SDR on flutter boundary is explained through aerodynamic work and the correlated shock wave location.The results show that the SDR delays the shock wave shifting downstream,and the Mach number corresponding to reaching freeze region increases as the split angle increases.Therefore,the peak value of aerodynamic moment coefficient amplitude and the sharp ascent process of phase occurs at higher Mach number,which leads to the delay in the occurrence of the transonic dip.Besides,before the transonic dip of airfoil without SDR occurs,the aerodynamic moment phase of airfoil with the SDR decreases slowly due to the decrease in the speed of shock wave moving downstream.This results in an increased flutter speed when employing the SDR before the transonic dip of airfoil without SDR occurs.Meanwhile,the effects of asymmetric split angles on the transonic flutter characteristics are also investigated.Before the transonic dip of airfoil without SDR occurs,the flutter characteristic is dominated by the smaller split angle.展开更多
An investigation on the ventral diverterless high offset S-shaped inlet is carried out at Mach numbers from 0.600 to 1.534, angles of attack from -4° to 9.4°, and yaw angles from 0° to 8°. Results ...An investigation on the ventral diverterless high offset S-shaped inlet is carried out at Mach numbers from 0.600 to 1.534, angles of attack from -4° to 9.4°, and yaw angles from 0° to 8°. Results indicate： （1） a large region of low total pressure exists at the lower part of the inlet exit caused by the counter-rotating vortices in the S-shaped duct; （2） the performances of the inlet at Mach number 1.000 reach almost the highest, so the propulsion system could work efficiently in terms of aerodynamics; （3） the total pressure recovery increases slowly at first and then remains unvaried as the Mach number rises from 0.6 to 1.0, however, it does in an opposite manner in the conventional diverter-equipped S-shaped inlet; （4） the performances of the inlet are generally insensitive to angles of attack from -4° to 9.4° and yaw angles from 0° to 8° at Mach number 0.850, and angles of attack from -2° to 6° and yaw angles from 0° to 5° at Mach number 1.534.展开更多
The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic inter...The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic interaction between the blade rows. The results show that, compared to the compressor with unurved rotors, the compressor under scrutiny acquires remarkable increases in efficiency with significantly reduced amplitudes of the time-dependent fluctuation. The amplitude of the pressure fluctuation around the stator leading edge decreases at both endwalls, but increases at the mid-span in the curved rotors. The pressure fluctuation near the stator leading edge, therefore, becomes more uniform in the radial direction of this compressor. Except for the leading edge area, the pressure fluctuatinn amplitude declines remarkably in the tip region of stator surface downstream of the curved rotor, but hardly changes in the middle and at the hub.展开更多
A coupled fluid-structure method is developed for flutter analysis of blade vibrations in turbomachinery. The approach is based on the time domain solution of the fluid-structure interaction in which the aerodynamic a...A coupled fluid-structure method is developed for flutter analysis of blade vibrations in turbomachinery. The approach is based on the time domain solution of the fluid-structure interaction in which the aerodynamic and structural equations are marched simultaneously in time. The three-dimensional （3D） unsteady Reynolds average Navier-Stokes （RANS） equations are solved with a multiblock finite volume scheme on dynamic deforming grids to evaluate the aerodynamic force. Dual time-stepping technique and an efficient implicit scheme with multigrid are employed to march the solution in time. The blade vibration is modeled with an aeroelasticity model in which blade motion is computed by linear combination of responses of each mode under unsteady loads. The code is validated in prediction of the unsteady flow flutter behavior of an oscillating cascade and is applied to flutter analysis of a transonic fan at the design speed.展开更多
In transonic flow,buffet is a phenomenon of flow instability caused by shock wave/boundary layer interaction and flow separation.The phenomenon is common in transonic flow,and it has serious impact on the structural s...In transonic flow,buffet is a phenomenon of flow instability caused by shock wave/boundary layer interaction and flow separation.The phenomenon is common in transonic flow,and it has serious impact on the structural strength and fatigue life of aircraft.In this paper,three typical airfoils:the supercritical OAT15A,the high-speed symmetrical NACA64A010,and the thin,transonic/supersonic NACA64A204 are selected as the research objects.The flow fields of these airfoils under pre-buffet and buffet onset conditions are simulated by Unsteady Reynolds Averaged Navier-Stokes (URANS) method,and the mode analysis of numerical results is carried out by Dynamic Mode Decomposition (DMD).Qualitative and quantitative analysis of the shock wave motion,shock wave intensity,shock foot bubble and trailing edge separation,and pressure coefficient fluctuation were performed to attain deep insight of transonic buffet flow features of different airfoils near buffet onset conditions.The results of DMD analysis show that the energy proportion of the steady mode of these airfoils decreases dramatically when approaching the buffet onset angle of attack,while the growth rate of the primary mode increases inversely.It was found that at the onset of buffet,there exist different degrees of merging behavior between shock foot bubble and trailing edge separation during one buffet cycle,and the instability of shock wave and separation induced shear layer are closely related to the merging behavior.展开更多
A Non-Intrusive Reduced-Order Model(NIROM)based on Proper Orthogonal Decomposition(POD)has been proposed for predicting the flow fields of transonic airfoils with geometry parameters.To provide a better reduced-order ...A Non-Intrusive Reduced-Order Model(NIROM)based on Proper Orthogonal Decomposition(POD)has been proposed for predicting the flow fields of transonic airfoils with geometry parameters.To provide a better reduced-order subspace to approximate the real flow field,a domain decomposition method has been used to separate the hard-to-predict regions from the full field and POD has been adopted in the regions individually.An Artificial Neural Network(ANN)has replaced the Radial Basis Function(RBF)to interpolate the coefficients of the POD modes,aiming at improving the approximation accuracy of the NIROM for non-samples.When predicting the flow fields of transonic airfoils,the proposed NIROM has demonstrated a high performance.展开更多
Transition prediction is a hot research topic of fluid mechanics.For subsonic and transonic aerodynamic flows,e^(N) method based on Linear Stability Theory(LST)is usually adopted reliably to predict transition.In 2013...Transition prediction is a hot research topic of fluid mechanics.For subsonic and transonic aerodynamic flows,e^(N) method based on Linear Stability Theory(LST)is usually adopted reliably to predict transition.In 2013,Coder and Maughmer established a transport equation for Tollmien-Schlichting(T-S)instability so that the e^(N) method can be applied to general Reynolds-Average-Navier-Stokes(RANS)solvers conveniently.However,this equation focuses on T-S instability,and is invalid for crossflow instability induced transition which plays a crucial role in flow instability of three-dimensional boundary layers.Subsequently,a transport equation for crossflow instability was developed in 2016,which is restricted to wing-like geometries.Then,in 2019,this model was extended to arbitrarily shaped geometries based on local variables.However,there are too many tedious functions and parameters in this version,and it can only be used for incompressible flows.Hence,in this paper,after a large amount of LST analyses and parameter optimization,an improved version for subsonic and transonic boundary layers is built.The present improved model is more robust and more concise,and it can be applied widely in aeronautical flows,which has great engineering application value and significance.An extensive validation study for this improved transition model will be performed.展开更多
In this paper, a new permeable adaptive integration surface is developed in order to evaluate transonic rotor noise in accordance with FW-H_pds equations(Ffowcs Williams-Hawkings equations with penetrable data surface...In this paper, a new permeable adaptive integration surface is developed in order to evaluate transonic rotor noise in accordance with FW-H_pds equations(Ffowcs Williams-Hawkings equations with penetrable data surface). Firstly, a nonlinear near-field solution is computed on the basis of Navier-Stokes equations, which is developed on moving-embedded grid methodology.The solution calculated through the present CFD method is used as the input for acoustic calculations by FW-H_pds equations. Then, two criteria for constructing integration surfaces are established based on the analysis of the quadrupole source strength and the nonlinear characteristic.A new surface is determined adaptively by the pressure gradient or density in a given flowfield,eschewing the uncertainties associated with determining cylinder-shaped integration surfaces. For varying hover cases, transonic noises are simulated with new integration surfaces for a UH-1 model rotor. Furthermore, numerical results of the new integration surface derived from the density perturbation value conform better to experimental data than results derived from the pressure gradient.Finally, the integration surface given by jrqj being 0.1, which is an applicable criterion obtained from hover cases, is used to predict transonic rotor noise in forward flight. The computational accuracy of the new integration surface method has been validated in predicting transonic rotor noise of an AH-1 model rotor at different advance ratios.展开更多
Current research shows that the traditional shock control bump（SCB） can weaken the intensity of shock and better the transonic buffet performance. The author finds that when SCB is placed downstream of the shock, it...Current research shows that the traditional shock control bump（SCB） can weaken the intensity of shock and better the transonic buffet performance. The author finds that when SCB is placed downstream of the shock, it can decrease the adverse pressure gradient. This may prevent the shock foot separation bubble to merge with the trailing edge separation and finally improve the buffet performance. Based on RAE2822 airfoil, two types of SCB are designed according to the two different mechanisms. By using Reynolds-averaged Navier-Stokes（RANS） and unsteady Reynolds-averaged Navier-Stokes（URANS） methods to analyze the properties of RAE2822 airfoil with and without SCB, the results show that the downstream SCB can better the buffet performance under a wide range of freestream Mach number and the steady aerodynamics characteristic is similar to that of RAE2822 airfoil. The traditional SCB can only weaken the intensity of the shock under the design condition. Under the off-design conditions, the SCB does not do much to or even worsen the buffet performance. Indeed, the use of backward bump can flatten the leeward side of the airfoil, and this is similar to the mechanism that supercritical airfoil can weaken the recompression of shock wave.展开更多
To simulate the transonic atomization jet process in Laval nozzles,to test the law of droplet atomization and distribution,to find a method of supersonic atomization for dust-removing nozzles,and to improve nozzle eff...To simulate the transonic atomization jet process in Laval nozzles,to test the law of droplet atomization and distribution,to find a method of supersonic atomization for dust-removing nozzles,and to improve nozzle efficiency,the finite element method has been used in this study based on the COMSOL computational fluid dynamics module.The study results showed that the process cannot be realized alone under the two-dimensional axisymmetric,three-dimensional and three-dimensional symmetric models,but it can be calculated with the transformation dimension method,which uses the parameter equations generated from the two-dimensional axisymmetric flow field data of the three-dimensional model.The visualization of this complex process,which is difficult to measure and analyze experimentally,was realized in this study.The physical process,macro phenomena and particle distribution of supersonic atomization are analyzed in combination with this simulation.The rationality of the simulation was verified by experiments.A new method for the study of the atomization process and the exploration of its mechanism in a compressible transonic speed flow field based on the Laval nozzle has been provided,and a numerical platform for the study of supersonic atomization dust removal has been established.展开更多
Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate f...Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate factorization algorithm and internal Newton iterations. An integral boundary layer method based on the dissipation integral is used to account for viscous effects. The computational results about unsteady transonic forces on wings, bodies and control surfaces are in agreement with experimental data.展开更多
Numerical simulations are performed to study the aeroelastic responses of an elastically suspended airfoil in transonic buffet flow, by coupling the unsteady Reynolds-averaged Navier- Stokes （RANS） equations and str...Numerical simulations are performed to study the aeroelastic responses of an elastically suspended airfoil in transonic buffet flow, by coupling the unsteady Reynolds-averaged Navier- Stokes （RANS） equations and structural motion equation. The current work focuses on the char- acteristic analysis of the lock-in phenomenon. Great attentions are paid to studying the frequency range of lock-in and the effects of the three parameters, namely the structural natural frequency, mass ratio and structural damping, on lock-in characteristic of the elastic system in detail. It is found that when the structural natural frequency is close to the buffet frequency, the coupling fre- quency of the elastic system is no longer equal to the buffet frequency, but keeps the same value as the structural natural frequency. The frequency lock-in occurs and stays present until the structural nature frequency is near the double buffet frequency. It means that the lock-in presents within a broad range, of which the lower threshold is near the buffet frequency, while the upper threshold is near the double buffet frequency. Moreover, the frequency range of lock-in is affected by mass ratio and structural damping. The lower the mass ratio and structural damping are, the wider the range of lock-in will be. The upper threshold of lock-in grows with the mass ratio and structural damping decreasing, but the lower threshold always keeps the same.展开更多
To predict the flutter dynamic pressure of a wind tunnel model before flutter test,an accurate Computational Fluid Dynamics/Computational Structural Dynamics(CFD/CSD)-based flutter prediction method is proposed under ...To predict the flutter dynamic pressure of a wind tunnel model before flutter test,an accurate Computational Fluid Dynamics/Computational Structural Dynamics(CFD/CSD)-based flutter prediction method is proposed under the conditions of a 2.4 m×2.4 m transonic wind tunnel with porous wall.From the CFD simulations of the flows through an inclined hole of this wind tunnel,the Nambu's linear porous wall model between the flow rate and the differential pressure is extended to the porous wall with inclined holes,so that the porous wall can be conveniently modeled as a boundary condition.According to the flutter testing approach for the current wind tunnel,the steady CFD calculation is conducted to achieve the required inlet Mach number.A timedomain CFD/CSD method is then employed to evaluate the structural response of the experimental model,and the critical flutter point is obtained by increasing the dynamic pressure step by step at a fixed Mach number.The present method is applied to the flutter calculations for a vertical tail model and an aircraft model tested in the current transonic wind tunnel.For both models,the computed flutter characteristics agree well with the experimental results.展开更多
Transonic rudder buzz responses based on the computational fluid dynamics or computational structural dynamics(CFD/CSD)loosely method are analyzed for a tailless flying wing unmanned aerial vehicle(UAV).The Reynolds-a...Transonic rudder buzz responses based on the computational fluid dynamics or computational structural dynamics(CFD/CSD)loosely method are analyzed for a tailless flying wing unmanned aerial vehicle(UAV).The Reynolds-averaged Navier-Stokes(RANS)equations and finite element methods based on the detailed aerodynamic and structural model are established,in which the aerodynamic dynamic meshes adopt the unstructured dynamic meshes based on the combination of spring-based smoothing and local remeshing methods,and the lower-upper symmetric-Gauss-Seidel(LU-SGS)iteration and Harten-Lax-van Leer-Einfeldt-Wada(HLLEW)space discrete methods based on the shear stress transport(SST)turbulence model are used to calculate the aerodynamic force.The constraints of the rudder motions are fixed at the end of structural model of the flying wing UAV,and the structural geometric nonlinearities are also considered in the flying wing UAV with a high aspect ratio.The interfaces between structural and aerodynamic models are built with an exact match surface where load transferring is performed based on 3Dinterpolation.The flying wing UAV transonic buzz responses based on the aerodynamic structural coupling method are studied,and the rudder buzz responses and aileron,elevator and flap vibration responses caused by rudder motion are also investigated.The effects of attack,height,rotating angular frequency and Mach number under transonic conditions on the flying wing UAV rudder buzz responses are discussed.The results can be regarded as a reference for the flying wing UAV engineering vibration analysis.展开更多
Partial surge is a type of instability inception in transonic compressors and occurs in the form of axisymmetric low-frequency disturbances localized in the hub region.Previous studies illustrate that the frequency of...Partial surge is a type of instability inception in transonic compressors and occurs in the form of axisymmetric low-frequency disturbances localized in the hub region.Previous studies illustrate that the frequency of partial surge is set by the Helmholtz frequency of the entire system,which motivates to propose a hypothesis that the system response performs an important role in the formation of partial surge.For further verification,a series of experiments are conducted to explore the link between the propagating of the partial surge and the system feedback in this study.In the first case,an additional test point is set on the wall of the plenum to detect the system response.Combining the flow behaviors inside the plenum with the disturbances in the rotor tip and stator hub/tip regions,the effects of the system feedback on the occurrence of the continuous disturbances and the rotating stall cells are illustrated.In the second case,a screen is mounted at the compressor outlet to prevent positive feedback from the plenum.The experimental results demonstrate that in the absence of system feedback,it is the occurrence of spike-type stall inception that leads to the flow instability instead of that of partial surge.In addition,three flow phenomena in the second case are discussed,including the occurrence of the single pulse,the unstable process during the stall evolution and the switch of instability inception.展开更多
An inviscid base pressure model for transonic turbine blade has been presented. It has been shown that for a given back pressure the base pressure at the trailing edge, and the profile loss of a turbine blade are fixe...An inviscid base pressure model for transonic turbine blade has been presented. It has been shown that for a given back pressure the base pressure at the trailing edge, and the profile loss of a turbine blade are fixed according to the model and the base pressure can be calculated with the help of an inviscid numerical scheme. A parameteric study on the model shows that a blade profile with positive curvature downstream of the throat is advantageous for generating less loss, whilst the worst situation is when the exit flow reaches the sonic condition.展开更多
A stall inception model for transonic fan/compressors is presented in this paper. It can be shown that under some assumptions the solution of unsteady flow field consists of pressure wave which propagates upstream or ...A stall inception model for transonic fan/compressors is presented in this paper. It can be shown that under some assumptions the solution of unsteady flow field consists of pressure wave which propagates upstream or downstream, vortex wave and entropy wave convected with the mean flow speed. By further using the mode-matching technique and applying the conservation law and conditions reflecting the loss characteristics of a compressor in the inlet and outlet of the rotor or stator blade rows, a group of homogeneous equations can be obtained from which the stability equation can be derived. Based on the analysis of the unsteady phenomenon caused by casing treatments, the function of casing treatments has been modeled by a wall impedance condition which has been included in the stability model through the eigenvalues and the corresponding eigenfunctions of the system. Besides, the effect of shock waves in cascade channel on the stability prediction is also considered in the stall inception model. Finally, some numerical analysis and experimental investigation are also conducted with emphasis on the mutual comparison.展开更多
文摘This study investigates the influence of mesh resolution and turbulence model selection on the accuracy of numerical simulations for transonic flow,with particular emphasis on shock-boundary layer interaction phenomena.Accurate prediction of such flows is notoriously difficult due to the sensitivity to near-wall resolution,global mesh density,and turbulence model assumptions,and this problem motivates the present work.Two solvers were employed,rhoCentralFoam(unsteady)and TSLAeroFoam(steady-state),both are compressible and density-based and implemented within the OpenFOAM framework.The investigation focuses on three different non-dimensional wall distance(y+)values of 1,2.5 and 5,each implemented with both moderate and fine mesh resolutions.Three turbulence models—Spalart-Allmaras(SA),k-ωShear Stress Transport(SST),and k- Realizable—were evaluated at M=0.74,Re=2.7×106,andα=3.19°.Results showed that while both solvers achieved good overall agreement with experimental data,particularly in terms of pressure distribution,lift coefficient,and shock location,noticeable differences still emerged.The k-ωSST model consistently delivered the most robust performance across all cases,capturing the shock position on y+≈1 meshes with deviations below0.02c compared to the experiment,and maintaining accuracy even at y+≈5.The k- Realizable model was highly sensitive to near-wall resolution,displacing shocks downstream at higher y+values,whereas Spalart-Allmaras remained broadly comparable to the k-ωSST model in predictive performance.The rhoCentralFoam solver achieved consistently better lift predictions,staying within about 2%of the experimental value on average,whereas TSLAeroFoam overpredicted it by around 4%.For transonic Reynolds-Averaged Navier-Stokes(RANS)simulations,unsteady k-ωSST with y+≈1 is recommended for maximum fidelity,whereas steady k-ωSST or SA simulations offer a practical option for quick and reasonably accurate aerodynamic predictions.
基金support of the National Natural Science Foundation of China(No.52322603)the Science Center for Gas Turbine Project of China(Nos.P2022-B-II-004-001 and P2023-B-II-001-001)+1 种基金the Fundamental Research Funds for the Central Universities,Chinathe Beijing Nova Program of China(Nos.20220484074 and 20230484479).
文摘The utilization of Inlet Guide Vane (IGV) plays a key factor in affecting the instability evolution. Existing literature mainly focuses on the effect of IGV on instability inception that occurs in the rotor region. However, with the emergence of compressor instability starting from the stator region, the mechanism of various instability inceptions that occurs in different blade rows due to the change of IGV angles should be further examined. In this study, experiments were focused on three types of instability inceptions observed previously in a 1.5-stage axial flow compressor. To analyze the conversion of stall evolutions, the compressor rotating speed was set to 17 160 r/min, at which both the blade loading in the stator hub region and rotor tip region were close to the critical value before final compressor stall. Meanwhile, the dynamic test points with high-response were placed to monitor the pressures both at the stator trailing edges and rotor tips. The results indicate that the variation of reaction determines the region where initial instability occurs. Indeed, negative pre-rotation of the inlet guide vane leads to high-reaction, initiating stall disturbance from the rotor region. Positive pre-rotation results in low-reaction, initiating stall disturbance from the stator region. Furthermore, the type of instability evolution is affected by the radial loading distribution under different IGV angles. Specifically, a spike-type inception occurs at the rotor blade tip with a large angle of attack at the rotor inlet (−2°, −4° and −6°). Meanwhile, the critical total pressure ratio at the rotor tip is 1.40 near stall. As the angle of attack decreases, the stator blade loading reaches its critical boundary, with a value of approximately 1.35. At this moment, if the rotor tip maintains high blade loading similar to the stator hub, the partial surge occurs (0° and +2°);otherwise, the hub instability occurs (+4° and +6°).
基金co-supported by the National Natural Science Foundation of China(Nos.52325602,52306036 and 52306035)the National Science and Technology Major Project of China(No.Y2022-II-0003-0006 and Y2022-II-0002-0005)+1 种基金the project funded by China Postdoctoral Science Foundation(No.2022M720346)supported by the Key Laboratory of Pre-Research Management Centre of China(No.6142702200101).
文摘The influence of Impedance Boundary Condition (IBC) on transonic compressors is investigated. A systematic input–output analytical framework is developed, which treats the nonlinearities as unknown forcing terms. The framework is validated through the experiments of rotating inlet distortion within a low-speed compressor. The input–output method is subsequently applied to transonic compressors, including NASA Rotor37 and Stage35, wherein impedance optimization is studied along with the exploration of its fundamental mechanisms. The IBC is employed to model the effect of Casing Treatment (CT). The optimal complex impedance values are determined through predicted results and tested across a range of circumferential modes and forcing frequencies. The IBC significantly reduces the energy and Reynolds stress gain, notably at the first-order circumferential mode and within the Rotor Rotating Frequency (RRF) range. Output modes reveal that transonic compressors with fine-tuned impedance values exhibit a more confined perturbation distribution and redistribute the perturbations compared to the uncontrolled case. Additionally, the roles of resistance and reactance are elucidated through input–output analysis, and resistance determines the energy transfer direction between flow and pressure waves and modulates the amplitude, whereas reactance modifies the phase relationships and attenuates the perturbations.
文摘In this paper,a series of flutter simulations are carried out to investigate the effects of split drag rudder(SDR)on the transonic flutter characteristic of rigid NACA 64A010.A structural dynamic model addressing two-degree-of-freedom pitch-plunge aeroelastic oscillations was coupled with the unsteady Reynolds-averaged Navier-Stokes equations to perform flutter simulation.Meanwhile,the influence mechanism of SDR on flutter boundary is explained through aerodynamic work and the correlated shock wave location.The results show that the SDR delays the shock wave shifting downstream,and the Mach number corresponding to reaching freeze region increases as the split angle increases.Therefore,the peak value of aerodynamic moment coefficient amplitude and the sharp ascent process of phase occurs at higher Mach number,which leads to the delay in the occurrence of the transonic dip.Besides,before the transonic dip of airfoil without SDR occurs,the aerodynamic moment phase of airfoil with the SDR decreases slowly due to the decrease in the speed of shock wave moving downstream.This results in an increased flutter speed when employing the SDR before the transonic dip of airfoil without SDR occurs.Meanwhile,the effects of asymmetric split angles on the transonic flutter characteristics are also investigated.Before the transonic dip of airfoil without SDR occurs,the flutter characteristic is dominated by the smaller split angle.
基金National Basic Research Program of China (5130802)
文摘An investigation on the ventral diverterless high offset S-shaped inlet is carried out at Mach numbers from 0.600 to 1.534, angles of attack from -4° to 9.4°, and yaw angles from 0° to 8°. Results indicate： （1） a large region of low total pressure exists at the lower part of the inlet exit caused by the counter-rotating vortices in the S-shaped duct; （2） the performances of the inlet at Mach number 1.000 reach almost the highest, so the propulsion system could work efficiently in terms of aerodynamics; （3） the total pressure recovery increases slowly at first and then remains unvaried as the Mach number rises from 0.6 to 1.0, however, it does in an opposite manner in the conventional diverter-equipped S-shaped inlet; （4） the performances of the inlet are generally insensitive to angles of attack from -4° to 9.4° and yaw angles from 0° to 8° at Mach number 0.850, and angles of attack from -2° to 6° and yaw angles from 0° to 5° at Mach number 1.534.
基金National Natural Science Foundation of China (506460210) Chinese Specialized Research Fund for the Doctoral Program of Higher Education (20060213007)Development Program for Outstanding Young Teachers in Harbin Institute of Technology (HITQNJS.2006.046)
文摘The unsteady 3D flow fields in a single-stage transonic compressor under designed conditions are simulated numerically to investigate the effects of the curved rotors on the stage performance and the aerodynamic interaction between the blade rows. The results show that, compared to the compressor with unurved rotors, the compressor under scrutiny acquires remarkable increases in efficiency with significantly reduced amplitudes of the time-dependent fluctuation. The amplitude of the pressure fluctuation around the stator leading edge decreases at both endwalls, but increases at the mid-span in the curved rotors. The pressure fluctuation near the stator leading edge, therefore, becomes more uniform in the radial direction of this compressor. Except for the leading edge area, the pressure fluctuatinn amplitude declines remarkably in the tip region of stator surface downstream of the curved rotor, but hardly changes in the middle and at the hub.
文摘A coupled fluid-structure method is developed for flutter analysis of blade vibrations in turbomachinery. The approach is based on the time domain solution of the fluid-structure interaction in which the aerodynamic and structural equations are marched simultaneously in time. The three-dimensional （3D） unsteady Reynolds average Navier-Stokes （RANS） equations are solved with a multiblock finite volume scheme on dynamic deforming grids to evaluate the aerodynamic force. Dual time-stepping technique and an efficient implicit scheme with multigrid are employed to march the solution in time. The blade vibration is modeled with an aeroelasticity model in which blade motion is computed by linear combination of responses of each mode under unsteady loads. The code is validated in prediction of the unsteady flow flutter behavior of an oscillating cascade and is applied to flutter analysis of a transonic fan at the design speed.
基金supported by the National Natural Science Foundation of China (No. 11802009)
文摘In transonic flow,buffet is a phenomenon of flow instability caused by shock wave/boundary layer interaction and flow separation.The phenomenon is common in transonic flow,and it has serious impact on the structural strength and fatigue life of aircraft.In this paper,three typical airfoils:the supercritical OAT15A,the high-speed symmetrical NACA64A010,and the thin,transonic/supersonic NACA64A204 are selected as the research objects.The flow fields of these airfoils under pre-buffet and buffet onset conditions are simulated by Unsteady Reynolds Averaged Navier-Stokes (URANS) method,and the mode analysis of numerical results is carried out by Dynamic Mode Decomposition (DMD).Qualitative and quantitative analysis of the shock wave motion,shock wave intensity,shock foot bubble and trailing edge separation,and pressure coefficient fluctuation were performed to attain deep insight of transonic buffet flow features of different airfoils near buffet onset conditions.The results of DMD analysis show that the energy proportion of the steady mode of these airfoils decreases dramatically when approaching the buffet onset angle of attack,while the growth rate of the primary mode increases inversely.It was found that at the onset of buffet,there exist different degrees of merging behavior between shock foot bubble and trailing edge separation during one buffet cycle,and the instability of shock wave and separation induced shear layer are closely related to the merging behavior.
基金supported by the National Natural Science Foundation of China(No.11802245).
文摘A Non-Intrusive Reduced-Order Model(NIROM)based on Proper Orthogonal Decomposition(POD)has been proposed for predicting the flow fields of transonic airfoils with geometry parameters.To provide a better reduced-order subspace to approximate the real flow field,a domain decomposition method has been used to separate the hard-to-predict regions from the full field and POD has been adopted in the regions individually.An Artificial Neural Network(ANN)has replaced the Radial Basis Function(RBF)to interpolate the coefficients of the POD modes,aiming at improving the approximation accuracy of the NIROM for non-samples.When predicting the flow fields of transonic airfoils,the proposed NIROM has demonstrated a high performance.
基金supported by the National Science Foundation for Young Scholars of China(No.:11802245)。
文摘Transition prediction is a hot research topic of fluid mechanics.For subsonic and transonic aerodynamic flows,e^(N) method based on Linear Stability Theory(LST)is usually adopted reliably to predict transition.In 2013,Coder and Maughmer established a transport equation for Tollmien-Schlichting(T-S)instability so that the e^(N) method can be applied to general Reynolds-Average-Navier-Stokes(RANS)solvers conveniently.However,this equation focuses on T-S instability,and is invalid for crossflow instability induced transition which plays a crucial role in flow instability of three-dimensional boundary layers.Subsequently,a transport equation for crossflow instability was developed in 2016,which is restricted to wing-like geometries.Then,in 2019,this model was extended to arbitrarily shaped geometries based on local variables.However,there are too many tedious functions and parameters in this version,and it can only be used for incompressible flows.Hence,in this paper,after a large amount of LST analyses and parameter optimization,an improved version for subsonic and transonic boundary layers is built.The present improved model is more robust and more concise,and it can be applied widely in aeronautical flows,which has great engineering application value and significance.An extensive validation study for this improved transition model will be performed.
基金supports of the National Natural Science Foundation of China (Nos. 11272150 and 11572156)the Scientific Research Innovation Program of Jiangsu Province of China (No. KYLX15_0250)
文摘In this paper, a new permeable adaptive integration surface is developed in order to evaluate transonic rotor noise in accordance with FW-H_pds equations(Ffowcs Williams-Hawkings equations with penetrable data surface). Firstly, a nonlinear near-field solution is computed on the basis of Navier-Stokes equations, which is developed on moving-embedded grid methodology.The solution calculated through the present CFD method is used as the input for acoustic calculations by FW-H_pds equations. Then, two criteria for constructing integration surfaces are established based on the analysis of the quadrupole source strength and the nonlinear characteristic.A new surface is determined adaptively by the pressure gradient or density in a given flowfield,eschewing the uncertainties associated with determining cylinder-shaped integration surfaces. For varying hover cases, transonic noises are simulated with new integration surfaces for a UH-1 model rotor. Furthermore, numerical results of the new integration surface derived from the density perturbation value conform better to experimental data than results derived from the pressure gradient.Finally, the integration surface given by jrqj being 0.1, which is an applicable criterion obtained from hover cases, is used to predict transonic rotor noise in forward flight. The computational accuracy of the new integration surface method has been validated in predicting transonic rotor noise of an AH-1 model rotor at different advance ratios.
文摘Current research shows that the traditional shock control bump（SCB） can weaken the intensity of shock and better the transonic buffet performance. The author finds that when SCB is placed downstream of the shock, it can decrease the adverse pressure gradient. This may prevent the shock foot separation bubble to merge with the trailing edge separation and finally improve the buffet performance. Based on RAE2822 airfoil, two types of SCB are designed according to the two different mechanisms. By using Reynolds-averaged Navier-Stokes（RANS） and unsteady Reynolds-averaged Navier-Stokes（URANS） methods to analyze the properties of RAE2822 airfoil with and without SCB, the results show that the downstream SCB can better the buffet performance under a wide range of freestream Mach number and the steady aerodynamics characteristic is similar to that of RAE2822 airfoil. The traditional SCB can only weaken the intensity of the shock under the design condition. Under the off-design conditions, the SCB does not do much to or even worsen the buffet performance. Indeed, the use of backward bump can flatten the leeward side of the airfoil, and this is similar to the mechanism that supercritical airfoil can weaken the recompression of shock wave.
基金Supported by the National Natural Science Foundation of China (NO: 51704146, 51274116, 51704145).
文摘To simulate the transonic atomization jet process in Laval nozzles,to test the law of droplet atomization and distribution,to find a method of supersonic atomization for dust-removing nozzles,and to improve nozzle efficiency,the finite element method has been used in this study based on the COMSOL computational fluid dynamics module.The study results showed that the process cannot be realized alone under the two-dimensional axisymmetric,three-dimensional and three-dimensional symmetric models,but it can be calculated with the transformation dimension method,which uses the parameter equations generated from the two-dimensional axisymmetric flow field data of the three-dimensional model.The visualization of this complex process,which is difficult to measure and analyze experimentally,was realized in this study.The physical process,macro phenomena and particle distribution of supersonic atomization are analyzed in combination with this simulation.The rationality of the simulation was verified by experiments.A new method for the study of the atomization process and the exploration of its mechanism in a compressible transonic speed flow field based on the Laval nozzle has been provided,and a numerical platform for the study of supersonic atomization dust removal has been established.
基金Aeronautical Science Foundation of China (99A52007)
文摘Algebraic methods and rapid deforming techniques are used to generate three-dimensional boundary-fitted dynamic grids for assemblies. The conservative full-potential equation is solved by a time-accurate approximate factorization algorithm and internal Newton iterations. An integral boundary layer method based on the dissipation integral is used to account for viscous effects. The computational results about unsteady transonic forces on wings, bodies and control surfaces are in agreement with experimental data.
基金supported by the project of the National Natural Science Foundation of China(No.11272262)
文摘Numerical simulations are performed to study the aeroelastic responses of an elastically suspended airfoil in transonic buffet flow, by coupling the unsteady Reynolds-averaged Navier- Stokes （RANS） equations and structural motion equation. The current work focuses on the char- acteristic analysis of the lock-in phenomenon. Great attentions are paid to studying the frequency range of lock-in and the effects of the three parameters, namely the structural natural frequency, mass ratio and structural damping, on lock-in characteristic of the elastic system in detail. It is found that when the structural natural frequency is close to the buffet frequency, the coupling fre- quency of the elastic system is no longer equal to the buffet frequency, but keeps the same value as the structural natural frequency. The frequency lock-in occurs and stays present until the structural nature frequency is near the double buffet frequency. It means that the lock-in presents within a broad range, of which the lower threshold is near the buffet frequency, while the upper threshold is near the double buffet frequency. Moreover, the frequency range of lock-in is affected by mass ratio and structural damping. The lower the mass ratio and structural damping are, the wider the range of lock-in will be. The upper threshold of lock-in grows with the mass ratio and structural damping decreasing, but the lower threshold always keeps the same.
基金supported by the National Natural Science Foundation of China(No.11872212)a project funded by the Priority Academic Program Development of Jiangsu Higher Education Institutions。
文摘To predict the flutter dynamic pressure of a wind tunnel model before flutter test,an accurate Computational Fluid Dynamics/Computational Structural Dynamics(CFD/CSD)-based flutter prediction method is proposed under the conditions of a 2.4 m×2.4 m transonic wind tunnel with porous wall.From the CFD simulations of the flows through an inclined hole of this wind tunnel,the Nambu's linear porous wall model between the flow rate and the differential pressure is extended to the porous wall with inclined holes,so that the porous wall can be conveniently modeled as a boundary condition.According to the flutter testing approach for the current wind tunnel,the steady CFD calculation is conducted to achieve the required inlet Mach number.A timedomain CFD/CSD method is then employed to evaluate the structural response of the experimental model,and the critical flutter point is obtained by increasing the dynamic pressure step by step at a fixed Mach number.The present method is applied to the flutter calculations for a vertical tail model and an aircraft model tested in the current transonic wind tunnel.For both models,the computed flutter characteristics agree well with the experimental results.
基金supported by the Natural Science Foundation of China(No.61074155)the Shaanxi Provincial Natural Science Foundation of China(No.2013JM015)
文摘Transonic rudder buzz responses based on the computational fluid dynamics or computational structural dynamics(CFD/CSD)loosely method are analyzed for a tailless flying wing unmanned aerial vehicle(UAV).The Reynolds-averaged Navier-Stokes(RANS)equations and finite element methods based on the detailed aerodynamic and structural model are established,in which the aerodynamic dynamic meshes adopt the unstructured dynamic meshes based on the combination of spring-based smoothing and local remeshing methods,and the lower-upper symmetric-Gauss-Seidel(LU-SGS)iteration and Harten-Lax-van Leer-Einfeldt-Wada(HLLEW)space discrete methods based on the shear stress transport(SST)turbulence model are used to calculate the aerodynamic force.The constraints of the rudder motions are fixed at the end of structural model of the flying wing UAV,and the structural geometric nonlinearities are also considered in the flying wing UAV with a high aspect ratio.The interfaces between structural and aerodynamic models are built with an exact match surface where load transferring is performed based on 3Dinterpolation.The flying wing UAV transonic buzz responses based on the aerodynamic structural coupling method are studied,and the rudder buzz responses and aileron,elevator and flap vibration responses caused by rudder motion are also investigated.The effects of attack,height,rotating angular frequency and Mach number under transonic conditions on the flying wing UAV rudder buzz responses are discussed.The results can be regarded as a reference for the flying wing UAV engineering vibration analysis.
基金the support of National Natural Science Foundation of China(Nos.51706008,51636001,51976005 and 52006002)National Science and Technology Major Project,China(No.2017-Ⅱ-0005-0018)Aeronautics Power Foundation,China(No.6141B09050375)。
文摘Partial surge is a type of instability inception in transonic compressors and occurs in the form of axisymmetric low-frequency disturbances localized in the hub region.Previous studies illustrate that the frequency of partial surge is set by the Helmholtz frequency of the entire system,which motivates to propose a hypothesis that the system response performs an important role in the formation of partial surge.For further verification,a series of experiments are conducted to explore the link between the propagating of the partial surge and the system feedback in this study.In the first case,an additional test point is set on the wall of the plenum to detect the system response.Combining the flow behaviors inside the plenum with the disturbances in the rotor tip and stator hub/tip regions,the effects of the system feedback on the occurrence of the continuous disturbances and the rotating stall cells are illustrated.In the second case,a screen is mounted at the compressor outlet to prevent positive feedback from the plenum.The experimental results demonstrate that in the absence of system feedback,it is the occurrence of spike-type stall inception that leads to the flow instability instead of that of partial surge.In addition,three flow phenomena in the second case are discussed,including the occurrence of the single pulse,the unstable process during the stall evolution and the switch of instability inception.
文摘An inviscid base pressure model for transonic turbine blade has been presented. It has been shown that for a given back pressure the base pressure at the trailing edge, and the profile loss of a turbine blade are fixed according to the model and the base pressure can be calculated with the help of an inviscid numerical scheme. A parameteric study on the model shows that a blade profile with positive curvature downstream of the throat is advantageous for generating less loss, whilst the worst situation is when the exit flow reaches the sonic condition.
基金National Natural Science Foundation of China (50736007, 51010007)
文摘A stall inception model for transonic fan/compressors is presented in this paper. It can be shown that under some assumptions the solution of unsteady flow field consists of pressure wave which propagates upstream or downstream, vortex wave and entropy wave convected with the mean flow speed. By further using the mode-matching technique and applying the conservation law and conditions reflecting the loss characteristics of a compressor in the inlet and outlet of the rotor or stator blade rows, a group of homogeneous equations can be obtained from which the stability equation can be derived. Based on the analysis of the unsteady phenomenon caused by casing treatments, the function of casing treatments has been modeled by a wall impedance condition which has been included in the stability model through the eigenvalues and the corresponding eigenfunctions of the system. Besides, the effect of shock waves in cascade channel on the stability prediction is also considered in the stall inception model. Finally, some numerical analysis and experimental investigation are also conducted with emphasis on the mutual comparison.
