Development of flapping wing aerial vehicle (FWAV) has been of interest in the aerospace community with ongoing research into unsteady and low Reynolds number aerodynamics based on the vortex lattice method. Most of t...Development of flapping wing aerial vehicle (FWAV) has been of interest in the aerospace community with ongoing research into unsteady and low Reynolds number aerodynamics based on the vortex lattice method. Most of the previous research has been about pitching and plunging motion of the FWAV. With pitching and flapping motion of FMAV, people usually study it by experiment, and little work has been done by numerical calculation. In this paper, three-dimension unsteady vortex lattice method is applied to study the lift and thrust of FWAV with pitching and flapping motion. The results show that: 1) Lift is mainly produced during down stroke, however, thrust is produced during both down stroke and upstroke. The lift and thrust produced during down stroke are much more than that produced during upstroke. 2) Lift and thrust increase with the increase of flapping frequency;3) Thrust increases with the increase of flapping amplitude, but the lift decreases with the increase of flapping amplitude;4) Lift and thrust increase with the increase of mean pitching angle, but the effect on lift is much more than on thrust. This research is helpful to understand the flight mechanism of birds, thus improving the design of FWAV simulating birds.展开更多
Numerical simulation of wing stall of a blended flying wing configuration at transonic speed was conducted using both delayed detached eddy simulation（DDES） and unsteady Reynolds-averaged Navier-Stokes（URANS） equa...Numerical simulation of wing stall of a blended flying wing configuration at transonic speed was conducted using both delayed detached eddy simulation（DDES） and unsteady Reynolds-averaged Navier-Stokes（URANS） equations methods based on the shear stress transport（SST） turbulence model for a free-stream Mach number 0.9 and a Reynolds number 9.6 × 10. A joint time step/grid density study is performed based on power spectrum density（PSD） analysis of the frequency content of forces or moments, and medium mesh and the normalized time scale0.010 were suggested for this simulation. The simulation results show that the DDES methods perform more precisely than the URANS method and the aerodynamic coefficient results from DDES method compare very well with the experiment data. The angle of attack of nonlinear vortex lift and abrupt wing stall of DDES results compare well with the experimental data. The flow structure of the DDES computation shows that the wing stall is caused mainly by the leeward vortex breakdown which occurred at x/x= 0.6 at angle of attack of 14°. The DDES methods show advantage in the simulation problem with separation flow. The computed result shows that a shock/vortex interaction is responsible for the wing stall caused by the vortex breakdown. The balance of the vortex strength and axial flow, and the shock strength, is examined to provide an explanation of the sensitivity of the breakdown location. Wing body thickness has a great influence on shock and shock/vortex interactions, which can make a significant difference to the vortex breakdown behavior and stall characteristic of the blended flying wing configuration.展开更多
The objective of the paper is to compute the optimal burn-out conditions and control requirements that would result in maximum down-range/cross-range performance of a waverider type hypersonic boost-glide(HBG) vehicle...The objective of the paper is to compute the optimal burn-out conditions and control requirements that would result in maximum down-range/cross-range performance of a waverider type hypersonic boost-glide(HBG) vehicle within the medium and intermediate ranges,and compare its performance with the performances of wing-body and lifting-body vehicles vis-a-vis the g-load and the integrated heat load experienced by vehicles for the medium-sized launch vehicle under study.Trajectory optimization studies were carried out by considering the heat rate and dynamic pressure constraints.The trajectory optimization problem is modeled as a nonlinear,multiphase,constraint optimal control problem and is solved using a hp-adaptive pseudospectral method.Detail modeling aspects of mass,aerodynamics and aerothermodynamics for the launch and glide vehicles have been discussed.It was found that the optimal burn-out angles for waverider and wing-body configurations are approximately 5° and 14.8°,respectively,for maximum down-range performance under the constraint heat rate environment.The down-range and cross-range performance of HBG waverider configuration is nearly 1.3 and 2 times that of wing-body configuration respectively.The integrated heat load experienced by the HBG waverider was found to be approximately an order of magnitude higher than that of a lifting-body configuration and 5 times that of a wing-body configuration.The footprints and corresponding heat loads and control requirements for the three types of glide vehicles are discussed for the medium range launch vehicle under consideration.展开更多
The aerodynamic forces and flow structure of a model insect wing is studied by solving the Navier-Stokes equations numerically.After an initial start from rest,the wing is made to execute an azimuthal rotation(sweepin...The aerodynamic forces and flow structure of a model insect wing is studied by solving the Navier-Stokes equations numerically.After an initial start from rest,the wing is made to execute an azimuthal rotation(sweeping)at a large angle of attack and constant angular velocity.The Reynolds number(Re)considered in the present note is 480(Re is based on the mean chord length of the wing and the speed at 60% wing length from the wing root).During the constant-speed sweeping motion,the stall is absent and large and approximately constant lift and drag coefficients can be maintained.The mechanism for the absence of the stall or the maintenance of large aerodynamic force coefficients is as follows.Soon after the initial start,a vortex ring,which consists of the leading-edge vortex(LEV),the starting vortex,and the two wing-tip vortices,is formed in the wake of the wing.During the subsequent motion of the wing,a base-to-tip spanwise flow converts the vorticity in the LEV to the wing tip and the LEV keeps an approximately constant strength.This prevents the LEV from shedding.As a result, the size of the vortex ring increases approximately linearly with time,resulting in an approximately constant time rate of the first moment of vorticity,or approximately constant lift and drag coefficients. The variation of the relative velocity along the wing span causes a pressure gradient along the wing- span.The base-to-tip spanwise flow is mainly maintained by the pressure-gradient force.展开更多
Sail is the core part of autonomous sailboat and wing sail is a new type of sail. Wing sail generates not only propulsion but also lateral force and heeling moment. The latter two will affect the navigation status and...Sail is the core part of autonomous sailboat and wing sail is a new type of sail. Wing sail generates not only propulsion but also lateral force and heeling moment. The latter two will affect the navigation status and bring resistance. Double sail can effectively reduce the center of wind pressure and heeling moment. In order to study the effect of distance between two sails, airfoil and attack angle on the total lift coefficient of double sail propulsion system, pressure coefficient distribution and lift coefficient calculation model have been established based on vortex panel method. By using the basic finite solution, the fluid dynamic forces on the two-dimensional sails are computed.The results show that, the distance in the range of 0 to 1 time chord length, when using the same airfoil in the fore and aft sail, the total lift coefficient of the double sail increases with the increase of distance, finally reaches a stable value in the range of one to three times chord length. Lift coefficients of thicker airfoils are more sensitive to the change of distance. The thicker the airfoil, the longer distance is required of the total lift coefficient toward stable.When different airfoils are adopted in fore and aft sail, the total lift coefficient increases with the increase of the thickness of aft sail. The smaller the thickness difference is, the more sensitive to the distance change the lift coefficient is. The thinner the fore sail is, the lower the influence will be on the lift coefficient of aft sail.展开更多
文摘Development of flapping wing aerial vehicle (FWAV) has been of interest in the aerospace community with ongoing research into unsteady and low Reynolds number aerodynamics based on the vortex lattice method. Most of the previous research has been about pitching and plunging motion of the FWAV. With pitching and flapping motion of FMAV, people usually study it by experiment, and little work has been done by numerical calculation. In this paper, three-dimension unsteady vortex lattice method is applied to study the lift and thrust of FWAV with pitching and flapping motion. The results show that: 1) Lift is mainly produced during down stroke, however, thrust is produced during both down stroke and upstroke. The lift and thrust produced during down stroke are much more than that produced during upstroke. 2) Lift and thrust increase with the increase of flapping frequency;3) Thrust increases with the increase of flapping amplitude, but the lift decreases with the increase of flapping amplitude;4) Lift and thrust increase with the increase of mean pitching angle, but the effect on lift is much more than on thrust. This research is helpful to understand the flight mechanism of birds, thus improving the design of FWAV simulating birds.
基金supported by the National Natural Science Foundation of China (No. 11372337)
文摘Numerical simulation of wing stall of a blended flying wing configuration at transonic speed was conducted using both delayed detached eddy simulation（DDES） and unsteady Reynolds-averaged Navier-Stokes（URANS） equations methods based on the shear stress transport（SST） turbulence model for a free-stream Mach number 0.9 and a Reynolds number 9.6 × 10. A joint time step/grid density study is performed based on power spectrum density（PSD） analysis of the frequency content of forces or moments, and medium mesh and the normalized time scale0.010 were suggested for this simulation. The simulation results show that the DDES methods perform more precisely than the URANS method and the aerodynamic coefficient results from DDES method compare very well with the experiment data. The angle of attack of nonlinear vortex lift and abrupt wing stall of DDES results compare well with the experimental data. The flow structure of the DDES computation shows that the wing stall is caused mainly by the leeward vortex breakdown which occurred at x/x= 0.6 at angle of attack of 14°. The DDES methods show advantage in the simulation problem with separation flow. The computed result shows that a shock/vortex interaction is responsible for the wing stall caused by the vortex breakdown. The balance of the vortex strength and axial flow, and the shock strength, is examined to provide an explanation of the sensitivity of the breakdown location. Wing body thickness has a great influence on shock and shock/vortex interactions, which can make a significant difference to the vortex breakdown behavior and stall characteristic of the blended flying wing configuration.
基金the Chinese Scholarship Council for supporting the research
文摘The objective of the paper is to compute the optimal burn-out conditions and control requirements that would result in maximum down-range/cross-range performance of a waverider type hypersonic boost-glide(HBG) vehicle within the medium and intermediate ranges,and compare its performance with the performances of wing-body and lifting-body vehicles vis-a-vis the g-load and the integrated heat load experienced by vehicles for the medium-sized launch vehicle under study.Trajectory optimization studies were carried out by considering the heat rate and dynamic pressure constraints.The trajectory optimization problem is modeled as a nonlinear,multiphase,constraint optimal control problem and is solved using a hp-adaptive pseudospectral method.Detail modeling aspects of mass,aerodynamics and aerothermodynamics for the launch and glide vehicles have been discussed.It was found that the optimal burn-out angles for waverider and wing-body configurations are approximately 5° and 14.8°,respectively,for maximum down-range performance under the constraint heat rate environment.The down-range and cross-range performance of HBG waverider configuration is nearly 1.3 and 2 times that of wing-body configuration respectively.The integrated heat load experienced by the HBG waverider was found to be approximately an order of magnitude higher than that of a lifting-body configuration and 5 times that of a wing-body configuration.The footprints and corresponding heat loads and control requirements for the three types of glide vehicles are discussed for the medium range launch vehicle under consideration.
基金The project supported by the National Natural Science Foundation of China(10232010)
文摘The aerodynamic forces and flow structure of a model insect wing is studied by solving the Navier-Stokes equations numerically.After an initial start from rest,the wing is made to execute an azimuthal rotation(sweeping)at a large angle of attack and constant angular velocity.The Reynolds number(Re)considered in the present note is 480(Re is based on the mean chord length of the wing and the speed at 60% wing length from the wing root).During the constant-speed sweeping motion,the stall is absent and large and approximately constant lift and drag coefficients can be maintained.The mechanism for the absence of the stall or the maintenance of large aerodynamic force coefficients is as follows.Soon after the initial start,a vortex ring,which consists of the leading-edge vortex(LEV),the starting vortex,and the two wing-tip vortices,is formed in the wake of the wing.During the subsequent motion of the wing,a base-to-tip spanwise flow converts the vorticity in the LEV to the wing tip and the LEV keeps an approximately constant strength.This prevents the LEV from shedding.As a result, the size of the vortex ring increases approximately linearly with time,resulting in an approximately constant time rate of the first moment of vorticity,or approximately constant lift and drag coefficients. The variation of the relative velocity along the wing span causes a pressure gradient along the wing- span.The base-to-tip spanwise flow is mainly maintained by the pressure-gradient force.
基金financially supported by the JIANG Xinsong Innovation Fund(Grant No.Y8F7010701)
文摘Sail is the core part of autonomous sailboat and wing sail is a new type of sail. Wing sail generates not only propulsion but also lateral force and heeling moment. The latter two will affect the navigation status and bring resistance. Double sail can effectively reduce the center of wind pressure and heeling moment. In order to study the effect of distance between two sails, airfoil and attack angle on the total lift coefficient of double sail propulsion system, pressure coefficient distribution and lift coefficient calculation model have been established based on vortex panel method. By using the basic finite solution, the fluid dynamic forces on the two-dimensional sails are computed.The results show that, the distance in the range of 0 to 1 time chord length, when using the same airfoil in the fore and aft sail, the total lift coefficient of the double sail increases with the increase of distance, finally reaches a stable value in the range of one to three times chord length. Lift coefficients of thicker airfoils are more sensitive to the change of distance. The thicker the airfoil, the longer distance is required of the total lift coefficient toward stable.When different airfoils are adopted in fore and aft sail, the total lift coefficient increases with the increase of the thickness of aft sail. The smaller the thickness difference is, the more sensitive to the distance change the lift coefficient is. The thinner the fore sail is, the lower the influence will be on the lift coefficient of aft sail.
