The paper examines the dynamic stall characteristics of a finite wing with an aspect ratio of eight in order to explore the 3D effects on flow topology,aerodynamic characteristics,and pitching damping.Firstly,CFD meth...The paper examines the dynamic stall characteristics of a finite wing with an aspect ratio of eight in order to explore the 3D effects on flow topology,aerodynamic characteristics,and pitching damping.Firstly,CFD methods are developed to calculate the aerodynamic characteristics of wings.The URANS equations are solved using a finite volume method,and the two-equation k-ωshear stress transport(SST)turbulence model is employed to account for viscosity effects.Secondly,the CFD methods are used to simulate the aerodynamic characteristics of both a static,rectangular wing and a pitching,tapered wing to verify their effectiveness and accuracy.The numerical results show good agreement with experimental data.Subsequently,the static and dynamic characteristics of the finite wing are computed and discussed.The results reveal significant 3D flow structures during both static and dynamic stalls,including wing tip vortices,arch vortices,Ω-type vortices,and ring vortices.These phenomena lead to differences in the aerodynamic characteristics of the finite wing compared with a 2D airfoil.Specifically,the finite wing has a smaller lift slope during attached-flow stages,higher stall angles,and more gradual stall behavior.Flow separation initially occurs in the middle spanwise section and gradually spreads to both ends.Regarding aerodynamic damping,the inboard sections mainly generate unstable loading.Furthermore,sections experiencing light stall have a higher tendency to produce negative damping compared with sections experiencing deep dynamic stall.展开更多
To reduce the high computational cost of the uncertainty analysis, a procedure is proposed for the aerodynamic optimization under uncertainties, in which the surrogate model is used to simplify the computation of the ...To reduce the high computational cost of the uncertainty analysis, a procedure is proposed for the aerodynamic optimization under uncertainties, in which the surrogate model is used to simplify the computation of the uncertainty analysis. The surrogate model is constructed by using the Latin Hypercube design and the Kriging model. The random parameters are used to account for the small manufacturing errors and the variations of operating conditions. Based on the surrogate model, an uncertainty analysis approach, called the Monte Carlo simulation, is used to compute the mean value and the variance of the predicated performance. The robust optimization for aerodynamic design is formulated, and solved by the genetic algorithm. And then, an airfoil optimization problem is used to test the proposed procedure. Results show that the optimal solutions obtained from the uncertainty-based optimization formulation are less sensitive to uncertainties. And the design constraints are still satisfied under the uncertainties.展开更多
The roughness effect based on the wall function method is introduced into the numerical simulation of the rime ice accretion and the resulting effect on the aerodynamic performance of the airfoil. Incorporating the tw...The roughness effect based on the wall function method is introduced into the numerical simulation of the rime ice accretion and the resulting effect on the aerodynamic performance of the airfoil. Incorporating the two-phase model of air/super-cooled droplets in the Eulerian coordinate system, this paper presents the simulation of the rime ice accretion on the NACA 0012 airfoil. The predicted rime ice shape is compared with those results of measurements and simulations by other icing codes. Also the resulting effects of rime ice on airfoil aerodynamic performance are discussed. Results indicate that the rime ice accretion leads to the loss of the maximum lift coefficient by 26%, the decrease of the stall angle by about 3° and the considerable increase of the drag coefficient.展开更多
This paper describes an improved model of curved vortex element on the circular arc (CVEC) for rotor wake analysis.As the key of the paper,two approximate formulas are derived by the series of limited terms to replace...This paper describes an improved model of curved vortex element on the circular arc (CVEC) for rotor wake analysis.As the key of the paper,two approximate formulas are derived by the series of limited terms to replace the Legendre incomplete elliptical integrals from the Biot-Savart integration,and the analytical solution of the induced velocity for the CVEC is obtained, which is more efficient in the complex rotor free wake calculation. Furthermore,the approximate formulas with the chosen factors are selected to avoid sigularity and give finite result of the induced velocity on the Vortex line,and an equivalent viscous vortex core radius might be evaluated.As examples, the induced velocity calculations on the vortex ring and two turns of a skew vortex helix are performed, and the comparisons between the circular-arc vortex element and the conventional straightline vortex element (SLVE) are given.It is shown that this curved vortex element model is advantageous over the SLVE model and is suitable for the rotor wake analysis.展开更多
The muhi-body analysis of the aeroelastic stability of the tiltrotor aircraft is presented. Muhi-body dynamic differential equations are combined with the equations of the unsteady dynamic inflow model to establish th...The muhi-body analysis of the aeroelastic stability of the tiltrotor aircraft is presented. Muhi-body dynamic differential equations are combined with the equations of the unsteady dynamic inflow model to establish the complete unsteadily aeroelastic coupling analytical model of the tiltrotor. The stability of the tiltrotor in the helicopter mode is analyzed aiming at a semi span soft-inplane tihrotor model with an elastic wing. Parametric effects of the lag stiffness of blades and the flight speed are analyzed. Numerical simulations demonstrate that the multibody analytical model can analyze the aeroelastic stability of the tiltrotor aircraft in the helicopter mode.展开更多
To calculate the aerodynamics of flapping-wing micro air vehicle（MAV） with the high efficiency and the engineering-oriented accuracy,an improved unsteady vortex lattice method （UVLM） for MAV is proposed. The metho...To calculate the aerodynamics of flapping-wing micro air vehicle（MAV） with the high efficiency and the engineering-oriented accuracy,an improved unsteady vortex lattice method （UVLM） for MAV is proposed. The method considers the influence of instantaneous wing deforming in flapping,as well as the induced drag,additionally models the stretching and the dissipation of vortex rings,and can present the aerodynamics status on the wing surface. An implementation of the method is developed. Moreover,the results and the efficiency of the proposed method are verified by CFD methods. Considering the less time cost of UVLM,for application of UVLM in the MAV optimization,the influence of wake vortex ignoring time saving and precision is studied. Results show that saving in CPU time with wake vortex ignoring the appropriate distance is considerable while the precision is not significantly reduced. It indicates the potential value of UVLM in the optimization of MAV design.展开更多
The aeromechanical st ability for the coupled rotor/fuselage system of helicopters in forward flight i s investigated. The periodic time-varying equations of motion are developed thr ough building a new 24DOF coupled ...The aeromechanical st ability for the coupled rotor/fuselage system of helicopters in forward flight i s investigated. The periodic time-varying equations of motion are developed thr ough building a new 24DOF coupled rigid/elastic blended element based on the fle xible multibody system theory in this paper. It accounts for the effects of prec one, sweep, and the moderately large elastic deflections on the blade and elasti city of shaft and fuselage of the helicopter. The dynamic coupling between the r igid motion of blades about the flap, lag and pitch hinges of articulated rotor and moderately large elastic deflections are included. There is no restriction o n the rotation amplitudes of flap, lag and pitch in the formulation. The stabili ty of periodic solution is studied using the Floquet theory. The transition matr ix is calculated by the Newmark integration method. The aeromechanical stability of a new helicopter is studied. The results show that it is stable in the given forward flight. But the instability arises with the decrease of the bending and torsion stiffness of the shaft.展开更多
We propose an efficient and robust algorithm to solve the steady Euler equa- tions on unstructured grids.The new algorithm is a Newton-iteration method in which each iteration step is a linear multigrid method using b...We propose an efficient and robust algorithm to solve the steady Euler equa- tions on unstructured grids.The new algorithm is a Newton-iteration method in which each iteration step is a linear multigrid method using block lower-upper symmetric Gauss-Seidel(LU-SGS)iteration as its smoother To regularize the Jacobian matrix of Newton-iteration,we adopted a local residual dependent regularization as the replace- ment of the standard time-stepping relaxation technique based on the local CFL number The proposed method can be extended to high order approximations and three spatial dimensions in a nature way.The solver was tested on a sequence of benchmark prob- lems on both quasi-uniform and local adaptive meshes.The numerical results illustrated the efficiency and robustness of our algorithm.展开更多
A comprehensive aeroelastic analysis of helicopter rotor with advanced tips based on finite element theory in space and time is coupled with an efficient rotor wake modeling. The effects of different wake models and t...A comprehensive aeroelastic analysis of helicopter rotor with advanced tips based on finite element theory in space and time is coupled with an efficient rotor wake modeling. The effects of different wake models and tip sweep angles on blade response and loads are investigated in forward flight. Each blade is assumed to undergo flap bending, lag bending and elastic twist deflections. The blade response is calculated from nonlinear periodic equations using a finite element in time scheme. For induced inflow distributions on the rotor disk, a constant vorticity contour wake model is used. Results show that blade response and loads are sensitive to wake model and tip sweep angle, and the rotor wake analysis is important for capturing the harmonics of vertical hub loads, flap bending moments and low speed aerodynamic loadings.展开更多
This paper concerns the real time control of the boundary layer on an aircraft wing. This new approach consists in heating the surface in an unsteady regime using electrically resistant strips embedded in the wing ski...This paper concerns the real time control of the boundary layer on an aircraft wing. This new approach consists in heating the surface in an unsteady regime using electrically resistant strips embedded in the wing skin. The control of the boundary layer＇s separation and transition point will provide a reduction in friction drag, and hence a reduction in fuel consumption. This new method consists in applying the required thermal power in the different strips in order to ensure the desired temperatures on the aircraft wing. We also have to determine the optimum size of these strips （length, width and distance between two strips）. This implies finding the best mathematical model corresponding to the physics enabling us to facilitate the calculation for any type of material used for the wings. Secondly, the heating being unsteady, and, as during a flight the flow conditions or the ambient temperatures vary, the thermal power needed changes and must be chosen as fast as possible in order to ensure optimal operating conditions.展开更多
In view of the reduction of hovering efficiency near high tension when a helicopter rotor hovers,a numerical simulation method of lifting rotor hovering aerodynamic characteristics based on leading edge droop is estab...In view of the reduction of hovering efficiency near high tension when a helicopter rotor hovers,a numerical simulation method of lifting rotor hovering aerodynamic characteristics based on leading edge droop is established in this paper. It is dominated by Reynolds average N-S equation in integral form. Firstly,VR-12 airfoil is taken as the research object,and the influence of leading edge droop angle on the aerodynamic characteristics of two-dimensional airfoil is studied. Secondly,the modified 7 A rotor is taken as the research object,and the effects of different leading edge droop angles at the position of blade r/R=0.75—1 on the aerodynamic characteristics in hover are explored. It is found that the leading edge droop can significantly improve the aerodynamic characteristics of two-dimensional airfoil and three-dimensional hovering rotor near high angle of attack,and can effectively inhibit the generation of stall vortex.展开更多
The purpose of increasing the aerodynamic efficiency and enhancing the supermaneuverability for the selected supersonic aircraft is presented. Aerodynamic characteristics, the surface pressure distribution and the max...The purpose of increasing the aerodynamic efficiency and enhancing the supermaneuverability for the selected supersonic aircraft is presented. Aerodynamic characteristics, the surface pressure distribution and the maximum lift are estimated for the baseline configuration for different Mach numbers and attack angles in subson- ic and supersonic potential flows, using a low-order three-dimensional panel method supported with the semi-empirical formulas of the data compendium （DATCOM）. Total nose-up and nose-down pitching moments about the center of gravity of the complete aircraft in the subsonic region depending on flight conditions and aircraft performance limitations are estimated. A software package is developed to implement the two-dimensional thrust vectoring flight control technique （pitch vectoring up and down） controlled by the advanced aerodynamic and control surface （the foreplane or the canard）. Results show that the canard with the thrust vectoring produces enough nose-down moment and can support the stabilizer at high maneuvers. The suggested surface can increase the aerodynamic efficiency （lift-to-drag ratio） of the baseline configuration by 5%-6% in subsonic and supersonic flight regimes.展开更多
基金supported by the National Natural Science Foundation of China(No.12072156)the National Key Laboratory Foundation of China(No.61422202103)the Priority Academic Program Development of Jiangsu Higher Education Institutions(PAPD)。
文摘The paper examines the dynamic stall characteristics of a finite wing with an aspect ratio of eight in order to explore the 3D effects on flow topology,aerodynamic characteristics,and pitching damping.Firstly,CFD methods are developed to calculate the aerodynamic characteristics of wings.The URANS equations are solved using a finite volume method,and the two-equation k-ωshear stress transport(SST)turbulence model is employed to account for viscosity effects.Secondly,the CFD methods are used to simulate the aerodynamic characteristics of both a static,rectangular wing and a pitching,tapered wing to verify their effectiveness and accuracy.The numerical results show good agreement with experimental data.Subsequently,the static and dynamic characteristics of the finite wing are computed and discussed.The results reveal significant 3D flow structures during both static and dynamic stalls,including wing tip vortices,arch vortices,Ω-type vortices,and ring vortices.These phenomena lead to differences in the aerodynamic characteristics of the finite wing compared with a 2D airfoil.Specifically,the finite wing has a smaller lift slope during attached-flow stages,higher stall angles,and more gradual stall behavior.Flow separation initially occurs in the middle spanwise section and gradually spreads to both ends.Regarding aerodynamic damping,the inboard sections mainly generate unstable loading.Furthermore,sections experiencing light stall have a higher tendency to produce negative damping compared with sections experiencing deep dynamic stall.
文摘To reduce the high computational cost of the uncertainty analysis, a procedure is proposed for the aerodynamic optimization under uncertainties, in which the surrogate model is used to simplify the computation of the uncertainty analysis. The surrogate model is constructed by using the Latin Hypercube design and the Kriging model. The random parameters are used to account for the small manufacturing errors and the variations of operating conditions. Based on the surrogate model, an uncertainty analysis approach, called the Monte Carlo simulation, is used to compute the mean value and the variance of the predicated performance. The robust optimization for aerodynamic design is formulated, and solved by the genetic algorithm. And then, an airfoil optimization problem is used to test the proposed procedure. Results show that the optimal solutions obtained from the uncertainty-based optimization formulation are less sensitive to uncertainties. And the design constraints are still satisfied under the uncertainties.
文摘The roughness effect based on the wall function method is introduced into the numerical simulation of the rime ice accretion and the resulting effect on the aerodynamic performance of the airfoil. Incorporating the two-phase model of air/super-cooled droplets in the Eulerian coordinate system, this paper presents the simulation of the rime ice accretion on the NACA 0012 airfoil. The predicted rime ice shape is compared with those results of measurements and simulations by other icing codes. Also the resulting effects of rime ice on airfoil aerodynamic performance are discussed. Results indicate that the rime ice accretion leads to the loss of the maximum lift coefficient by 26%, the decrease of the stall angle by about 3° and the considerable increase of the drag coefficient.
文摘This paper describes an improved model of curved vortex element on the circular arc (CVEC) for rotor wake analysis.As the key of the paper,two approximate formulas are derived by the series of limited terms to replace the Legendre incomplete elliptical integrals from the Biot-Savart integration,and the analytical solution of the induced velocity for the CVEC is obtained, which is more efficient in the complex rotor free wake calculation. Furthermore,the approximate formulas with the chosen factors are selected to avoid sigularity and give finite result of the induced velocity on the Vortex line,and an equivalent viscous vortex core radius might be evaluated.As examples, the induced velocity calculations on the vortex ring and two turns of a skew vortex helix are performed, and the comparisons between the circular-arc vortex element and the conventional straightline vortex element (SLVE) are given.It is shown that this curved vortex element model is advantageous over the SLVE model and is suitable for the rotor wake analysis.
文摘The muhi-body analysis of the aeroelastic stability of the tiltrotor aircraft is presented. Muhi-body dynamic differential equations are combined with the equations of the unsteady dynamic inflow model to establish the complete unsteadily aeroelastic coupling analytical model of the tiltrotor. The stability of the tiltrotor in the helicopter mode is analyzed aiming at a semi span soft-inplane tihrotor model with an elastic wing. Parametric effects of the lag stiffness of blades and the flight speed are analyzed. Numerical simulations demonstrate that the multibody analytical model can analyze the aeroelastic stability of the tiltrotor aircraft in the helicopter mode.
基金Supported by the Aviation Science Foundation of China (2007ZA56001)the National Natural Science Foundation of China(50865009)~~
文摘To calculate the aerodynamics of flapping-wing micro air vehicle（MAV） with the high efficiency and the engineering-oriented accuracy,an improved unsteady vortex lattice method （UVLM） for MAV is proposed. The method considers the influence of instantaneous wing deforming in flapping,as well as the induced drag,additionally models the stretching and the dissipation of vortex rings,and can present the aerodynamics status on the wing surface. An implementation of the method is developed. Moreover,the results and the efficiency of the proposed method are verified by CFD methods. Considering the less time cost of UVLM,for application of UVLM in the MAV optimization,the influence of wake vortex ignoring time saving and precision is studied. Results show that saving in CPU time with wake vortex ignoring the appropriate distance is considerable while the precision is not significantly reduced. It indicates the potential value of UVLM in the optimization of MAV design.
文摘The aeromechanical st ability for the coupled rotor/fuselage system of helicopters in forward flight i s investigated. The periodic time-varying equations of motion are developed thr ough building a new 24DOF coupled rigid/elastic blended element based on the fle xible multibody system theory in this paper. It accounts for the effects of prec one, sweep, and the moderately large elastic deflections on the blade and elasti city of shaft and fuselage of the helicopter. The dynamic coupling between the r igid motion of blades about the flap, lag and pitch hinges of articulated rotor and moderately large elastic deflections are included. There is no restriction o n the rotation amplitudes of flap, lag and pitch in the formulation. The stabili ty of periodic solution is studied using the Floquet theory. The transition matr ix is calculated by the Newmark integration method. The aeromechanical stability of a new helicopter is studied. The results show that it is stable in the given forward flight. But the instability arises with the decrease of the bending and torsion stiffness of the shaft.
文摘We propose an efficient and robust algorithm to solve the steady Euler equa- tions on unstructured grids.The new algorithm is a Newton-iteration method in which each iteration step is a linear multigrid method using block lower-upper symmetric Gauss-Seidel(LU-SGS)iteration as its smoother To regularize the Jacobian matrix of Newton-iteration,we adopted a local residual dependent regularization as the replace- ment of the standard time-stepping relaxation technique based on the local CFL number The proposed method can be extended to high order approximations and three spatial dimensions in a nature way.The solver was tested on a sequence of benchmark prob- lems on both quasi-uniform and local adaptive meshes.The numerical results illustrated the efficiency and robustness of our algorithm.
文摘A comprehensive aeroelastic analysis of helicopter rotor with advanced tips based on finite element theory in space and time is coupled with an efficient rotor wake modeling. The effects of different wake models and tip sweep angles on blade response and loads are investigated in forward flight. Each blade is assumed to undergo flap bending, lag bending and elastic twist deflections. The blade response is calculated from nonlinear periodic equations using a finite element in time scheme. For induced inflow distributions on the rotor disk, a constant vorticity contour wake model is used. Results show that blade response and loads are sensitive to wake model and tip sweep angle, and the rotor wake analysis is important for capturing the harmonics of vertical hub loads, flap bending moments and low speed aerodynamic loadings.
文摘This paper concerns the real time control of the boundary layer on an aircraft wing. This new approach consists in heating the surface in an unsteady regime using electrically resistant strips embedded in the wing skin. The control of the boundary layer＇s separation and transition point will provide a reduction in friction drag, and hence a reduction in fuel consumption. This new method consists in applying the required thermal power in the different strips in order to ensure the desired temperatures on the aircraft wing. We also have to determine the optimum size of these strips （length, width and distance between two strips）. This implies finding the best mathematical model corresponding to the physics enabling us to facilitate the calculation for any type of material used for the wings. Secondly, the heating being unsteady, and, as during a flight the flow conditions or the ambient temperatures vary, the thermal power needed changes and must be chosen as fast as possible in order to ensure optimal operating conditions.
基金supported by the National Natural Science Foundation of China(No.11972190)the Aeronautical Science Foundation of China(No. 20185752)
文摘In view of the reduction of hovering efficiency near high tension when a helicopter rotor hovers,a numerical simulation method of lifting rotor hovering aerodynamic characteristics based on leading edge droop is established in this paper. It is dominated by Reynolds average N-S equation in integral form. Firstly,VR-12 airfoil is taken as the research object,and the influence of leading edge droop angle on the aerodynamic characteristics of two-dimensional airfoil is studied. Secondly,the modified 7 A rotor is taken as the research object,and the effects of different leading edge droop angles at the position of blade r/R=0.75—1 on the aerodynamic characteristics in hover are explored. It is found that the leading edge droop can significantly improve the aerodynamic characteristics of two-dimensional airfoil and three-dimensional hovering rotor near high angle of attack,and can effectively inhibit the generation of stall vortex.
文摘The purpose of increasing the aerodynamic efficiency and enhancing the supermaneuverability for the selected supersonic aircraft is presented. Aerodynamic characteristics, the surface pressure distribution and the maximum lift are estimated for the baseline configuration for different Mach numbers and attack angles in subson- ic and supersonic potential flows, using a low-order three-dimensional panel method supported with the semi-empirical formulas of the data compendium （DATCOM）. Total nose-up and nose-down pitching moments about the center of gravity of the complete aircraft in the subsonic region depending on flight conditions and aircraft performance limitations are estimated. A software package is developed to implement the two-dimensional thrust vectoring flight control technique （pitch vectoring up and down） controlled by the advanced aerodynamic and control surface （the foreplane or the canard）. Results show that the canard with the thrust vectoring produces enough nose-down moment and can support the stabilizer at high maneuvers. The suggested surface can increase the aerodynamic efficiency （lift-to-drag ratio） of the baseline configuration by 5%-6% in subsonic and supersonic flight regimes.
