To improve the cruise flight performance of aircraft, two new configurations of plasma actuators(grid-type and super-dense array) were investigated to reduce the turbulent skin friction drag of a low-speed airfoil. Th...To improve the cruise flight performance of aircraft, two new configurations of plasma actuators(grid-type and super-dense array) were investigated to reduce the turbulent skin friction drag of a low-speed airfoil. The induced jet characteristics of the two actuators in quiescent air were diagnosed with high-speed particle image velocimetry(PIV), and their drag reduction efficiencies were examined under different operating conditions in a wind tunnel. The results showed that the grid-type plasma actuator was capable of producing a wall-normal jet array(peak magnitude: 1.07 m/s) similar to that generated in a micro-blowing technique, while the superdense array plasma actuator created a wavy wall-parallel jet(magnitude: 0.94 m/s) due to the discrete spanwise electrostatic forces. Under a comparable electrical power consumption level,the super-dense array plasma actuator array significantly outperformed the grid-type configuration,reducing the total airfoil friction drag by approximately 22% at a free-stream velocity of 20 m/s.The magnitude of drag reduction was proportional to the dimensionless jet velocity ratio(r), and a threshold r = 0.014 existed under which little impact on airfoil drag could be discerned.展开更多
A numerical method has been used to analyze the flow field related to a NACA 0015 airfoil with and without a flap and assess the influence of the flap height and angle on the surface pressure coefficient,lift coeffici...A numerical method has been used to analyze the flow field related to a NACA 0015 airfoil with and without a flap and assess the influence of the flap height and angle on the surface pressure coefficient,lift coefficient,and drag coefficient.The numerical results demonstrate that the flap can effectively improve the lift coefficient of the airfoil;however,at small attack angles,its influence is significantly reduced.When the angle of attack exceeds the critical stall angle and the flap height is 1.5%of the chord length,the influence of the flap becomes very evident.As the flap height increases,the starting point of the separation vortex gradually moves forward and generates a larger wake vortex.Optimal aerodynamic characteristics are obtained for 1.5%(of the chord length)flap height and a 45°flap angle;in this case,the separation vortex is effectively reduced.展开更多
The current work aims at employing a gradient descent algorithm for optimizing the thrust of a flapping wing. An in-house solver has been employed, along with mesh movement methodologies to capture the dynamics of flo...The current work aims at employing a gradient descent algorithm for optimizing the thrust of a flapping wing. An in-house solver has been employed, along with mesh movement methodologies to capture the dynamics of flow around the airfoil. An efficient framework for implementing the coupled solver and optimization in a multicore environment has been implemented for the generation of optimized solutionsmaximizing thrust performance & computational speed.展开更多
Aerodynamic force and flow structures of two airfoils in a tandem configuration in flapping motions axe studied, by solving the Navier-Stokes equations in moving overset grids. Three typical phase differences between ...Aerodynamic force and flow structures of two airfoils in a tandem configuration in flapping motions axe studied, by solving the Navier-Stokes equations in moving overset grids. Three typical phase differences between the fore- and aft-airfoil flapping cycles are considered. It is shown that: (1) in the case of no interaction (single airfoil), the time average of the vertical force coefficient over the downstroke is 2.74, which is about 3 times as large as the maximum steady-state lift coefficient of a dragonfly wing; the time average of the horizontal force coefficient is 1.97, which is also large. The reasons for the large force coefficients are the acceleration at the beginning of a stroke, the delayed stall and the 'pitching-up' motion near the end of the stroke. (2) In the cases of two-airfoils, the time-variations of the force and moment coefficients on each airfoil are broadly similar to that of the single airfoil in that the vertical force is mainly produced in downstroke and the horizontal force in upstroke, but very large differences exist due to the interaction. (3) For in-phase stroking, the major differences caused by the interaction are that the vertical force on FA in downstroke is increased and the horizontal force on FA in upstroke decreased. As a result, the magnitude of the resultant force is almost unchanged but it inclines less forward. (4) For counter stroking, the major differences are that the vertical force on AA in downstroke and the horizontal force on FA in upstroke are decreased. As a result, the magnitude of the resultant force is decreased by about 20 percent but its direction is almost unchanged. (5) For 90 degrees -phase-difference stroking, the major differences axe that the vertical force on AA in downstroke and the horizontal force on FA in upstroke axe decreased greatly and the horizontal force on AA in upstroke increased. As a result, the magnitude of the resultant force is decreased by about 28% and it inclines more forward. (6) Among the three cases of phase angles, inphase flapping produces the largest vertical force (also the largest resultant force); the 90 degrees -phase-difference flapping results in the largest horizontal force, but the smallest resultant force.展开更多
In the present work,a parametric numerical study is conducted in order to assess the effect of airfoil cambering on the aerodynamic performance of rigid heaving airfoils.The incompressible Navier-Stokes equations are ...In the present work,a parametric numerical study is conducted in order to assess the effect of airfoil cambering on the aerodynamic performance of rigid heaving airfoils.The incompressible Navier-Stokes equations are solved in their velocity-pressure formulation using a second-order accurate in space and time finite-difference scheme.To tackle the problem of moving boundaries,the governing equations are solved on overlapping structured grids.The numerical simulations are performed at a Reynolds number of Re=1100 and at different values of Strouhal number and reduced frequency.The results obtained show that the airfoil cambering geometric parameter has a strong influence on the average lift coefficient,while it has a smaller impact on the average thrust coefficient and propulsive efficiency of heaving airfoils.展开更多
Simulations have been done to assess the lift, thrust and propulsive efficiency of different types of non-symmetrical airfoils under different flapping configurations. The variables involved are reduced frequency, Str...Simulations have been done to assess the lift, thrust and propulsive efficiency of different types of non-symmetrical airfoils under different flapping configurations. The variables involved are reduced frequency, Strouhal number, pitch amplitude and phase angle. In order to analyze the variables more efficiently, the design of experiments using the response surface methodology is applied. Results show that both the variables and shape of the airfoil have a profound effect on the lift, thrust, and efficiency. By using non- symmetrical airfoils, average lift coefficient as high as 2.23 can be obtained. The average thrust coefficient and efficiency also reach high values of 2.53 and 0.61, respectively. The lift production is highly dependent on the airfoil's shape while thrust production is influenced more heavily by the variables. Efficiency falls somewhere in between. Two-factor interac- tions are found to exist among the variables. This shows that it is not sufficient to analyze each variable individually. Vorticity diagrams are analyzed to explain the results obtained. Overall, the S1020 airfoil is able to provide relatively good efficiency and at the same time generate high thrust and lift force. These results aid in the design of a better ornithopter's wing.展开更多
The objective of this project is to improve the performance of the efficiency, thrust and lift of flapping wings in tandem arrangement. This research investigates the effect of the arrangement of the airfoils in tande...The objective of this project is to improve the performance of the efficiency, thrust and lift of flapping wings in tandem arrangement. This research investigates the effect of the arrangement of the airfoils in tandem on the performance of the airfoils by varying the phase difference and distance between the airfoils. Three flapping configurations from an earlier phase of a research which gives high efficiency, thrust and lift are used in the tandem simulation. It is found all the different flapping configurations show improvement in the efficiency, thrust or lift when the distance between the two airfoils and the phase angle between the heaving positions of the two airfoils are optimal. The average thrust coefficient of the tandem arrangement managed to attain more than twice that of the single one (4.84 vs. 2.05). On the other hand, the average lift coefficient of the tandem arrangement also increased to 4.59, as compared to the original single airfoil value of 3.04. All these results obtained will aid in the design of a better ornithopter with tandem wing arrangement.展开更多
Trailing edge serrations(TESs)are capable of noticeably suppressing the turbulent trailing edge noise induced by rotating wind turbine blades and become an integral part of a blade.However,the challenges involved in t...Trailing edge serrations(TESs)are capable of noticeably suppressing the turbulent trailing edge noise induced by rotating wind turbine blades and become an integral part of a blade.However,the challenges involved in the dimensional design of serration height 2 h,wavelengthλand flap angleϕare yet to be dealt with in a satisfactory manner.To address the problem,a general model for simulating the effects of serrations on the hydrodynamic and aeroacoustic performance is proposed due to its ease of use and relatively low requirements for user input.The solid serrations are replicated by momentum sources calculated by its aerodynamic forces.Then,a case relevant to wind turbine airfoil is examined,a hybrid improved delay detached eddy simulation(IDDES)method coupled with FW-H integration is deployed to obtain the flow features and far-field sound pressure level.It is found that the modeling method could reproduce the flow field and noise as serrated airfoil.展开更多
A bi-objective optimization problem for flapping airfoils is solved to maximize the time-averaged thrust coefficient and the propulsive efficiency. Design variables include the plunging amplitude, the pitching amplitu...A bi-objective optimization problem for flapping airfoils is solved to maximize the time-averaged thrust coefficient and the propulsive efficiency. Design variables include the plunging amplitude, the pitching amplitude and the phase shift angle. A well defined Kriging model is used to substitute the time-consuming high fidelity model, and a multi-objective genetic algorithm is employed as the search algorithm. The optimization results show that the propulsive efficiency can be improved by reducing the plunging amplitude and the phase shift angle in a proper way. The results of global sensitivity analysis using the Sobol’s method show that both of the time-averaged thrust coefficient and the propulsive efficiency are most sensitive to the plunging amplitude, and second most sensitive to the pitching amplitude. It is also observed that the phase shift angle has an un-negligible influence on the propulsive efficiency, and has little effect on the time-averaged thrust coefficient.展开更多
The paper is to integrate aerodynamic and aero-acoustic optimizatiom design of high lift devices,especially for two-element airfoils with slat.Aerodynamic analysis on flow field utilizes a high-order,high-resolution s...The paper is to integrate aerodynamic and aero-acoustic optimizatiom design of high lift devices,especially for two-element airfoils with slat.Aerodynamic analysis on flow field utilizes a high-order,high-resolution spatial differential method for large eddy simulation(LES),which can guarantee accuracy and efficiency.The aeroacoustic analysis for noise level is calculated with Ffowcs Williams-Hawkings(FW-H)integration formula.Fidelity of calculation is verified by standard models.Method of streamline-based Euler simulation(MSES)is used to obtain the aerodynamic characters.Based on the confirmation of numerical methods,detailed research has been conducted for the leading edge slat on multi-element airfoils.Various slot parameter influences on noise are analyzed.The results of the slot optimization parameters can be used in multi-element airfoil design.展开更多
In this paper, a new unsteady aerodynamic design method is presented based on the Navier-Stokes equations and a continuous adjoint approach. A basic framework of time-accurate unsteady airfoil optimization which adopt...In this paper, a new unsteady aerodynamic design method is presented based on the Navier-Stokes equations and a continuous adjoint approach. A basic framework of time-accurate unsteady airfoil optimization which adopts time-averaged aerodynamic coefficients as objective functions is presented. The time-accurate continuous adjoint equation and its boundary conditions are derived. The flow field and the adjoint equation are simulated numerically by the finite volume method (FVM). Feasibility and accuracy of the approach are perfectly validated by the design optimization results of the plunging NACA0012 airfoil.展开更多
A hybrid Cartesian grid/gridless method is developed for calculating viscous flows over multi-element airfoils.The method adopts an unstructured Cartesian grid to cover most areas of the computational domain and leave...A hybrid Cartesian grid/gridless method is developed for calculating viscous flows over multi-element airfoils.The method adopts an unstructured Cartesian grid to cover most areas of the computational domain and leaves only small region adjacent to the aerodynamic bodies to be filled with the cloud of points used in the gridless methods,which results in a better combination of the computational efficiency of the Cartesian grid and the flexibility of the gridless method in handling complex geometries.The clouds of points in the local gridless region are implemented in an anisotropic way according to the features of the thin boundary layer of the viscous flows over the airfoils,and the clouds of points at the vicinity of the interface between the grid and the gridless regions are also controlled by using an adaptive refinement technique during the generation of the unstructured Cartesian grid.An implementation of the resulting hybrid method is presented for solving two-dimensional compressible Navier-Stokes(NS)equations.The simulations of the viscous flows over a RAE2822airfoil or a two-element airfoil are successfully carried out,and the obtained results agree well with the available experimental data.展开更多
In this paper, on the basis of the variational principles developed the finite element method (FEM) is employed for numerical solution of the inverse pro blem of 2 D unsteady compressible flow around oscillating a...In this paper, on the basis of the variational principles developed the finite element method (FEM) is employed for numerical solution of the inverse pro blem of 2 D unsteady compressible flow around oscillating airfoils by incorporating the non reflecting far field boundary conditions and a new unsteady Kutta condition. All unknown boundary (airfoil contour) and discontinuities(shocks and free trailing vortex sheets) are determined via the functional variation with variable domain and artificial density concept. For the numerical realization of the variable domain variation, a special finite element with self adjusting nodes is also suggested herein. The numerical results show that the present method is effective for the design of unsteady airfoil.展开更多
This paper presents a numerical prediction of the unsteady flow field around oscillating airfoils at high angles of attack by solving unsteady Reynolds-averaged Navier-Stokes equations with SST turbulence model in ord...This paper presents a numerical prediction of the unsteady flow field around oscillating airfoils at high angles of attack by solving unsteady Reynolds-averaged Navier-Stokes equations with SST turbulence model in order to simulate the effects of wind tunnel model vibrations on the aerodynamic properties of airfoils,especially high-aspect-ratio wings in a wind tunnel.The effects of the phase lagging between different modes of oscillations,i.e.,the airfoil plunging oscillation mode,the pitching oscillation mode,and the forward-backward oscillation mode,are also studied.It is shown that the vibrations (oscillations) of airfoils can cause the unsteady shedding of large-size separated vortex to precede the stationary stall incidence,hence lead to a stall onset at some earlier (lower) incidence than that in the steady sense.The different phase lagging has different effect on the flow field.When the pitching oscillation mode has small phase lagging behind the plunging oscillation mode,the effect of vibrations is large.Besides,if the amplitude of the oscillations is large enough,and the different modes of vibrations match or combine appropriately,the unsteady stall may occur 2° earlier in angle of attack than the case where airfoils keep stationary.展开更多
In the present work, the recent developed Unsteady Double Wake Model, USDWM, is used to simulate separated flows past a wind turbine airfoil at high angles of attack. The solver is basically an unsteady two-dimensiona...In the present work, the recent developed Unsteady Double Wake Model, USDWM, is used to simulate separated flows past a wind turbine airfoil at high angles of attack. The solver is basically an unsteady two-dimensional panel method which uses the unsteady double wake technique to model flow separation and its dynamics. In this paper, the calculated integral forces have been successfully validated against wind tunnel measurements for the FFA-W3-211 airfoil. Furthermore, the computed highly unsteady flow field is analyzed in detail for a set of angles of attack ranging from light to deep stall conditions.展开更多
The performances of turbine blades have a significant impact on the energy conversion efficiency of vertical solar power plants.In the present study,such a relationship is assessed by considering two kinds of airfoil ...The performances of turbine blades have a significant impact on the energy conversion efficiency of vertical solar power plants.In the present study,such a relationship is assessed by considering two kinds of airfoil blades,designed by using the Wilson theory.In particular,numerical simulations are conducted using the SST K−ω model and assuming a wind speed of 3–6 m/s and seven or eight blades.The two airfoils are the NACA63121(with a larger chord length)and the AMES63212;It is shown that the torsion angle of the former is smaller,and its wind drag ratio is larger;furthermore,the resistance is increased by about 66.3%on average.Within the scope of the study,the results show that the NACA63212 airfoil is better than the AMES63212 airfoil in terms of power,with an average improvement of about 2.8%.The simulation results have a certain guiding significance for selecting turbine blade airfoils and improving turbine efficiency.展开更多
An efficient data-driven approach for predicting steady airfoil flows is proposed based on the Fourier neural operator(FNO),which is a new framework of neural networks.Theoretical reasons and experimental results are ...An efficient data-driven approach for predicting steady airfoil flows is proposed based on the Fourier neural operator(FNO),which is a new framework of neural networks.Theoretical reasons and experimental results are provided to support the necessity and effectiveness of the improvements made to the FNO,which involve using an additional branch neural operator to approximate the contribution of boundary conditions to steady solutions.The proposed approach runs several orders of magnitude faster than the traditional numerical methods.The predictions for flows around airfoils and ellipses demonstrate the superior accuracy and impressive speed of this novel approach.Furthermore,the property of zero-shot super-resolution enables the proposed approach to overcome the limitations of predicting airfoil flows with Cartesian grids,thereby improving the accuracy in the near-wall region.There is no doubt that the unprecedented speed and accuracy in forecasting steady airfoil flows have massive benefits for airfoil design and optimization.展开更多
The plenum chamber of a heat setting machine is a key structure for distributing hot air to different air channels.Its outlet velocity uniformity directly determines the heating uniformity of textiles,significantly af...The plenum chamber of a heat setting machine is a key structure for distributing hot air to different air channels.Its outlet velocity uniformity directly determines the heating uniformity of textiles,significantly affecting the heat setting performance.In a traditional heat setting machine,the outlet airflow maldistribution of the plenum chamber still exists.In this study,a novel plenum chamber with an airfoil baffle was established to improve the uniformity of the velocity distribution at the outlet in a heat setting machine.The structural influence of the plenum chamber on the velocity distribution was investigated using a computational fluid dynamics program.It was found that a chamber with a smaller outlet partition thickness had a better outlet velocity uniformity.The structural optimization of the plenum chamber was conducted using the particle swarm optimization algorithm.The outlet partition thickness,the transverse distance and the longitudinal distance of the optimized plenum chamber were 20,686.2 and 274.6 mm,respectively.Experiments were carried out.The experimental and simulated results showed that the optimized plenum chamber with an airfoil baffle could improve the outlet velocity uniformity.The air outlet velocity uniformity index of the optimized plenum chamber with an airfoil baffle was 4.75%higher than that of the plenum chamber without an airfoil baffle and 5.98%higher than that of the conventional chamber with a square baffle in a commercial heat setting machine.展开更多
Numerical simulations are conducted to improve the energy acquisition efficiency of H-type vertical axis wind turbines through the optimization of the related blade airfoil aerodynamic performance.The Bézier curve...Numerical simulations are conducted to improve the energy acquisition efficiency of H-type vertical axis wind turbines through the optimization of the related blade airfoil aerodynamic performance.The Bézier curve is initi-ally used tofit the curve profile of a NACA2412 airfoil,and the moving asymptote algorithm is then exploited to optimize the design of the considered H-type vertical-axis wind-turbine blade airfoil for a certain attack angle.The results show that the maximum lift coefficient of the optimized airfoil is 8.33%higher than that of the original airfoil.The maximum lift-to-drag ratio of the optimized airfoil exceeds the maximum lift-to-drag ratio of the ori-ginal airfoil by 11.22%.Moreover,the power coefficient is increased by 12.19%and the torque coefficient of the wind turbine is significantly improved.展开更多
High-speed rotor rotation under the low-density condition creates a special low-Reynolds compressible flow around the rotor blade airfoil where the compressibility effect on the laminar separated shear layer occurs. H...High-speed rotor rotation under the low-density condition creates a special low-Reynolds compressible flow around the rotor blade airfoil where the compressibility effect on the laminar separated shear layer occurs. However, the compressibility effect and shock wave generation associated with the increase in the Mach number (M) and the trend change due to their interference have not been clarified. The purpose is to clear the compressibility effect and its impact of shock wave generation on the flow field and aerodynamics. Therefore, we perform a two-dimensional unsteady calculation by Computational fluid dynamics (CFD) analysis using the CLF5605 airfoil used in the Mars helicopter Ingenuity, which succeeded in its first flight on Mars. The calculation conditions are set to the Reynolds number (Re) at 75% rotor span in hovering (Re = 15,400), and the Mach number was varied from incompressible (M = 0.2) to transonic (M = 1.2). The compressible fluid dynamics solver FaSTAR developed by the Japan aerospace exploration agency (JAXA) is used, and calculations are performed under multiple conditions in which the Mach number and angle of attack (α) are swept. The results show that a flow field is similar to that in the Earth’s atmosphere above M = 1.0, such as bow shock at the leading edge, whereas multiple λ-type shock waves are observed over the separated shear layer above α = 3° at M = 0.80. However, no significant difference is found in the C<sub>p</sub> distribution around the airfoil between M = 0.6 and M = 0.8. From the results, it is found that multiple λ-type shock waves have no significant effect on the airfoil surface pressure distribution, the separated shear layer effect is dominant in the surface pressure change and aerodynamic characteristics.展开更多
基金supported by National Natural Science Foundation of China (Nos.12002384, U2341277,and 52025064)Foundation Strengthening Program (No.2021JJ-0786)。
文摘To improve the cruise flight performance of aircraft, two new configurations of plasma actuators(grid-type and super-dense array) were investigated to reduce the turbulent skin friction drag of a low-speed airfoil. The induced jet characteristics of the two actuators in quiescent air were diagnosed with high-speed particle image velocimetry(PIV), and their drag reduction efficiencies were examined under different operating conditions in a wind tunnel. The results showed that the grid-type plasma actuator was capable of producing a wall-normal jet array(peak magnitude: 1.07 m/s) similar to that generated in a micro-blowing technique, while the superdense array plasma actuator created a wavy wall-parallel jet(magnitude: 0.94 m/s) due to the discrete spanwise electrostatic forces. Under a comparable electrical power consumption level,the super-dense array plasma actuator array significantly outperformed the grid-type configuration,reducing the total airfoil friction drag by approximately 22% at a free-stream velocity of 20 m/s.The magnitude of drag reduction was proportional to the dimensionless jet velocity ratio(r), and a threshold r = 0.014 existed under which little impact on airfoil drag could be discerned.
基金supported by the National Natural Science Foundation Project(Grant Numbers 51966018 and 51466015)the Key Research&Development Program of Xinjiang(Grant Number 2022B01003).
文摘A numerical method has been used to analyze the flow field related to a NACA 0015 airfoil with and without a flap and assess the influence of the flap height and angle on the surface pressure coefficient,lift coefficient,and drag coefficient.The numerical results demonstrate that the flap can effectively improve the lift coefficient of the airfoil;however,at small attack angles,its influence is significantly reduced.When the angle of attack exceeds the critical stall angle and the flap height is 1.5%of the chord length,the influence of the flap becomes very evident.As the flap height increases,the starting point of the separation vortex gradually moves forward and generates a larger wake vortex.Optimal aerodynamic characteristics are obtained for 1.5%(of the chord length)flap height and a 45°flap angle;in this case,the separation vortex is effectively reduced.
文摘The current work aims at employing a gradient descent algorithm for optimizing the thrust of a flapping wing. An in-house solver has been employed, along with mesh movement methodologies to capture the dynamics of flow around the airfoil. An efficient framework for implementing the coupled solver and optimization in a multicore environment has been implemented for the generation of optimized solutionsmaximizing thrust performance & computational speed.
文摘Aerodynamic force and flow structures of two airfoils in a tandem configuration in flapping motions axe studied, by solving the Navier-Stokes equations in moving overset grids. Three typical phase differences between the fore- and aft-airfoil flapping cycles are considered. It is shown that: (1) in the case of no interaction (single airfoil), the time average of the vertical force coefficient over the downstroke is 2.74, which is about 3 times as large as the maximum steady-state lift coefficient of a dragonfly wing; the time average of the horizontal force coefficient is 1.97, which is also large. The reasons for the large force coefficients are the acceleration at the beginning of a stroke, the delayed stall and the 'pitching-up' motion near the end of the stroke. (2) In the cases of two-airfoils, the time-variations of the force and moment coefficients on each airfoil are broadly similar to that of the single airfoil in that the vertical force is mainly produced in downstroke and the horizontal force in upstroke, but very large differences exist due to the interaction. (3) For in-phase stroking, the major differences caused by the interaction are that the vertical force on FA in downstroke is increased and the horizontal force on FA in upstroke decreased. As a result, the magnitude of the resultant force is almost unchanged but it inclines less forward. (4) For counter stroking, the major differences are that the vertical force on AA in downstroke and the horizontal force on FA in upstroke are decreased. As a result, the magnitude of the resultant force is decreased by about 20 percent but its direction is almost unchanged. (5) For 90 degrees -phase-difference stroking, the major differences axe that the vertical force on AA in downstroke and the horizontal force on FA in upstroke axe decreased greatly and the horizontal force on AA in upstroke increased. As a result, the magnitude of the resultant force is decreased by about 28% and it inclines more forward. (6) Among the three cases of phase angles, inphase flapping produces the largest vertical force (also the largest resultant force); the 90 degrees -phase-difference flapping results in the largest horizontal force, but the smallest resultant force.
基金supported by Maric Curie actions EST project FLUBIO(Grant:MEST-CT-2005-020228)support of the HPC-Europa++ project(Project number:211437)support of the European Community-Research Infrastructure Action of the FP7
文摘In the present work,a parametric numerical study is conducted in order to assess the effect of airfoil cambering on the aerodynamic performance of rigid heaving airfoils.The incompressible Navier-Stokes equations are solved in their velocity-pressure formulation using a second-order accurate in space and time finite-difference scheme.To tackle the problem of moving boundaries,the governing equations are solved on overlapping structured grids.The numerical simulations are performed at a Reynolds number of Re=1100 and at different values of Strouhal number and reduced frequency.The results obtained show that the airfoil cambering geometric parameter has a strong influence on the average lift coefficient,while it has a smaller impact on the average thrust coefficient and propulsive efficiency of heaving airfoils.
文摘Simulations have been done to assess the lift, thrust and propulsive efficiency of different types of non-symmetrical airfoils under different flapping configurations. The variables involved are reduced frequency, Strouhal number, pitch amplitude and phase angle. In order to analyze the variables more efficiently, the design of experiments using the response surface methodology is applied. Results show that both the variables and shape of the airfoil have a profound effect on the lift, thrust, and efficiency. By using non- symmetrical airfoils, average lift coefficient as high as 2.23 can be obtained. The average thrust coefficient and efficiency also reach high values of 2.53 and 0.61, respectively. The lift production is highly dependent on the airfoil's shape while thrust production is influenced more heavily by the variables. Efficiency falls somewhere in between. Two-factor interac- tions are found to exist among the variables. This shows that it is not sufficient to analyze each variable individually. Vorticity diagrams are analyzed to explain the results obtained. Overall, the S1020 airfoil is able to provide relatively good efficiency and at the same time generate high thrust and lift force. These results aid in the design of a better ornithopter's wing.
文摘The objective of this project is to improve the performance of the efficiency, thrust and lift of flapping wings in tandem arrangement. This research investigates the effect of the arrangement of the airfoils in tandem on the performance of the airfoils by varying the phase difference and distance between the airfoils. Three flapping configurations from an earlier phase of a research which gives high efficiency, thrust and lift are used in the tandem simulation. It is found all the different flapping configurations show improvement in the efficiency, thrust or lift when the distance between the two airfoils and the phase angle between the heaving positions of the two airfoils are optimal. The average thrust coefficient of the tandem arrangement managed to attain more than twice that of the single one (4.84 vs. 2.05). On the other hand, the average lift coefficient of the tandem arrangement also increased to 4.59, as compared to the original single airfoil value of 3.04. All these results obtained will aid in the design of a better ornithopter with tandem wing arrangement.
基金This work was supported by the National Natural Science Foundation of China(Grant No.51736008)“Transformational Technologies for Clean Energy and Demonstration”,Strategic Priority Research Program of the Chinese Academy of Sciences(Grant No.XDA21050303).
文摘Trailing edge serrations(TESs)are capable of noticeably suppressing the turbulent trailing edge noise induced by rotating wind turbine blades and become an integral part of a blade.However,the challenges involved in the dimensional design of serration height 2 h,wavelengthλand flap angleϕare yet to be dealt with in a satisfactory manner.To address the problem,a general model for simulating the effects of serrations on the hydrodynamic and aeroacoustic performance is proposed due to its ease of use and relatively low requirements for user input.The solid serrations are replicated by momentum sources calculated by its aerodynamic forces.Then,a case relevant to wind turbine airfoil is examined,a hybrid improved delay detached eddy simulation(IDDES)method coupled with FW-H integration is deployed to obtain the flow features and far-field sound pressure level.It is found that the modeling method could reproduce the flow field and noise as serrated airfoil.
基金Supported by the National Science Foundation for Post-doctoral Scientists of China (20090460216 )the National Defense Fundamental Research Foundation of China(B222006060)
文摘A bi-objective optimization problem for flapping airfoils is solved to maximize the time-averaged thrust coefficient and the propulsive efficiency. Design variables include the plunging amplitude, the pitching amplitude and the phase shift angle. A well defined Kriging model is used to substitute the time-consuming high fidelity model, and a multi-objective genetic algorithm is employed as the search algorithm. The optimization results show that the propulsive efficiency can be improved by reducing the plunging amplitude and the phase shift angle in a proper way. The results of global sensitivity analysis using the Sobol’s method show that both of the time-averaged thrust coefficient and the propulsive efficiency are most sensitive to the plunging amplitude, and second most sensitive to the pitching amplitude. It is also observed that the phase shift angle has an un-negligible influence on the propulsive efficiency, and has little effect on the time-averaged thrust coefficient.
文摘The paper is to integrate aerodynamic and aero-acoustic optimizatiom design of high lift devices,especially for two-element airfoils with slat.Aerodynamic analysis on flow field utilizes a high-order,high-resolution spatial differential method for large eddy simulation(LES),which can guarantee accuracy and efficiency.The aeroacoustic analysis for noise level is calculated with Ffowcs Williams-Hawkings(FW-H)integration formula.Fidelity of calculation is verified by standard models.Method of streamline-based Euler simulation(MSES)is used to obtain the aerodynamic characters.Based on the confirmation of numerical methods,detailed research has been conducted for the leading edge slat on multi-element airfoils.Various slot parameter influences on noise are analyzed.The results of the slot optimization parameters can be used in multi-element airfoil design.
基金supported by the Shanghai Municipal Natural Science Foundation(No.13ZR1415700)
文摘In this paper, a new unsteady aerodynamic design method is presented based on the Navier-Stokes equations and a continuous adjoint approach. A basic framework of time-accurate unsteady airfoil optimization which adopts time-averaged aerodynamic coefficients as objective functions is presented. The time-accurate continuous adjoint equation and its boundary conditions are derived. The flow field and the adjoint equation are simulated numerically by the finite volume method (FVM). Feasibility and accuracy of the approach are perfectly validated by the design optimization results of the plunging NACA0012 airfoil.
基金Supported by the National Natural Science Foundation of China(11172134)the Funding of Jiangsu Innovation Program for Graduate Education(CXZZ110192)the Priority Academic Program Development of Jiangsu Higher Education Institutions
文摘A hybrid Cartesian grid/gridless method is developed for calculating viscous flows over multi-element airfoils.The method adopts an unstructured Cartesian grid to cover most areas of the computational domain and leaves only small region adjacent to the aerodynamic bodies to be filled with the cloud of points used in the gridless methods,which results in a better combination of the computational efficiency of the Cartesian grid and the flexibility of the gridless method in handling complex geometries.The clouds of points in the local gridless region are implemented in an anisotropic way according to the features of the thin boundary layer of the viscous flows over the airfoils,and the clouds of points at the vicinity of the interface between the grid and the gridless regions are also controlled by using an adaptive refinement technique during the generation of the unstructured Cartesian grid.An implementation of the resulting hybrid method is presented for solving two-dimensional compressible Navier-Stokes(NS)equations.The simulations of the viscous flows over a RAE2822airfoil or a two-element airfoil are successfully carried out,and the obtained results agree well with the available experimental data.
文摘In this paper, on the basis of the variational principles developed the finite element method (FEM) is employed for numerical solution of the inverse pro blem of 2 D unsteady compressible flow around oscillating airfoils by incorporating the non reflecting far field boundary conditions and a new unsteady Kutta condition. All unknown boundary (airfoil contour) and discontinuities(shocks and free trailing vortex sheets) are determined via the functional variation with variable domain and artificial density concept. For the numerical realization of the variable domain variation, a special finite element with self adjusting nodes is also suggested herein. The numerical results show that the present method is effective for the design of unsteady airfoil.
基金Supported by National Natural Science Foundation of China (Contract/Grant Number:10572120,10872171)
文摘This paper presents a numerical prediction of the unsteady flow field around oscillating airfoils at high angles of attack by solving unsteady Reynolds-averaged Navier-Stokes equations with SST turbulence model in order to simulate the effects of wind tunnel model vibrations on the aerodynamic properties of airfoils,especially high-aspect-ratio wings in a wind tunnel.The effects of the phase lagging between different modes of oscillations,i.e.,the airfoil plunging oscillation mode,the pitching oscillation mode,and the forward-backward oscillation mode,are also studied.It is shown that the vibrations (oscillations) of airfoils can cause the unsteady shedding of large-size separated vortex to precede the stationary stall incidence,hence lead to a stall onset at some earlier (lower) incidence than that in the steady sense.The different phase lagging has different effect on the flow field.When the pitching oscillation mode has small phase lagging behind the plunging oscillation mode,the effect of vibrations is large.Besides,if the amplitude of the oscillations is large enough,and the different modes of vibrations match or combine appropriately,the unsteady stall may occur 2° earlier in angle of attack than the case where airfoils keep stationary.
基金The authors would like to acknowledge the support from the Danish Council for Strategic Research for the project Off Wind China(Sagsnr.0603-00506B)the European Union’s Seventh Programme for research,technological development and demonstration for the project“AVATAR:AdVanced Aerodynamic Tools for large Rotors”(FP7-ENERGY-2013-1/no.608396).
文摘In the present work, the recent developed Unsteady Double Wake Model, USDWM, is used to simulate separated flows past a wind turbine airfoil at high angles of attack. The solver is basically an unsteady two-dimensional panel method which uses the unsteady double wake technique to model flow separation and its dynamics. In this paper, the calculated integral forces have been successfully validated against wind tunnel measurements for the FFA-W3-211 airfoil. Furthermore, the computed highly unsteady flow field is analyzed in detail for a set of angles of attack ranging from light to deep stall conditions.
文摘The performances of turbine blades have a significant impact on the energy conversion efficiency of vertical solar power plants.In the present study,such a relationship is assessed by considering two kinds of airfoil blades,designed by using the Wilson theory.In particular,numerical simulations are conducted using the SST K−ω model and assuming a wind speed of 3–6 m/s and seven or eight blades.The two airfoils are the NACA63121(with a larger chord length)and the AMES63212;It is shown that the torsion angle of the former is smaller,and its wind drag ratio is larger;furthermore,the resistance is increased by about 66.3%on average.Within the scope of the study,the results show that the NACA63212 airfoil is better than the AMES63212 airfoil in terms of power,with an average improvement of about 2.8%.The simulation results have a certain guiding significance for selecting turbine blade airfoils and improving turbine efficiency.
文摘An efficient data-driven approach for predicting steady airfoil flows is proposed based on the Fourier neural operator(FNO),which is a new framework of neural networks.Theoretical reasons and experimental results are provided to support the necessity and effectiveness of the improvements made to the FNO,which involve using an additional branch neural operator to approximate the contribution of boundary conditions to steady solutions.The proposed approach runs several orders of magnitude faster than the traditional numerical methods.The predictions for flows around airfoils and ellipses demonstrate the superior accuracy and impressive speed of this novel approach.Furthermore,the property of zero-shot super-resolution enables the proposed approach to overcome the limitations of predicting airfoil flows with Cartesian grids,thereby improving the accuracy in the near-wall region.There is no doubt that the unprecedented speed and accuracy in forecasting steady airfoil flows have massive benefits for airfoil design and optimization.
基金National Natural Science Foundation of China(No.62173307)the Key R&D Projects of Science and Technology Department of Zhejiang Province,China(Nos.2023C01158,2022C01065 and 2022C01188)the Fundamental Research Funds of Zhejiang Sci-Tech University,China(No.22242298-Y)。
文摘The plenum chamber of a heat setting machine is a key structure for distributing hot air to different air channels.Its outlet velocity uniformity directly determines the heating uniformity of textiles,significantly affecting the heat setting performance.In a traditional heat setting machine,the outlet airflow maldistribution of the plenum chamber still exists.In this study,a novel plenum chamber with an airfoil baffle was established to improve the uniformity of the velocity distribution at the outlet in a heat setting machine.The structural influence of the plenum chamber on the velocity distribution was investigated using a computational fluid dynamics program.It was found that a chamber with a smaller outlet partition thickness had a better outlet velocity uniformity.The structural optimization of the plenum chamber was conducted using the particle swarm optimization algorithm.The outlet partition thickness,the transverse distance and the longitudinal distance of the optimized plenum chamber were 20,686.2 and 274.6 mm,respectively.Experiments were carried out.The experimental and simulated results showed that the optimized plenum chamber with an airfoil baffle could improve the outlet velocity uniformity.The air outlet velocity uniformity index of the optimized plenum chamber with an airfoil baffle was 4.75%higher than that of the plenum chamber without an airfoil baffle and 5.98%higher than that of the conventional chamber with a square baffle in a commercial heat setting machine.
基金This study was supported by the following research funding.Natural Science Foundation of Anhui Province,China,Grant Number 1908085ME166Research on the Key Technology of Multipole Grain Sampling and Inspection Equipment Based on Machine Vision,Anhui Provincial Grain Machinery Rural Development Collaborative Technology Service Center,Grant Number GXXT-2022-077+3 种基金Research on the Preparation Process and Application of Biochar Made of Bamboo,Science and Technology Bureau of Chuzhou City,Grant Number 2022ZN014The Development and Industrialization of Fruit Sorting Equipment,Science and Technology Bureau of Chuzhou City,Grant Number 2022ZN016Natural Science Major Project of Anhui Provincial Education Department,Anhui Provincial Education Department,Grant Number 2022AH040238Key Scientific Research Project of Anhui Provincial Education Department,Anhui Provincial Education Department,Grant Number KJ2021A0877.
文摘Numerical simulations are conducted to improve the energy acquisition efficiency of H-type vertical axis wind turbines through the optimization of the related blade airfoil aerodynamic performance.The Bézier curve is initi-ally used tofit the curve profile of a NACA2412 airfoil,and the moving asymptote algorithm is then exploited to optimize the design of the considered H-type vertical-axis wind-turbine blade airfoil for a certain attack angle.The results show that the maximum lift coefficient of the optimized airfoil is 8.33%higher than that of the original airfoil.The maximum lift-to-drag ratio of the optimized airfoil exceeds the maximum lift-to-drag ratio of the ori-ginal airfoil by 11.22%.Moreover,the power coefficient is increased by 12.19%and the torque coefficient of the wind turbine is significantly improved.
文摘High-speed rotor rotation under the low-density condition creates a special low-Reynolds compressible flow around the rotor blade airfoil where the compressibility effect on the laminar separated shear layer occurs. However, the compressibility effect and shock wave generation associated with the increase in the Mach number (M) and the trend change due to their interference have not been clarified. The purpose is to clear the compressibility effect and its impact of shock wave generation on the flow field and aerodynamics. Therefore, we perform a two-dimensional unsteady calculation by Computational fluid dynamics (CFD) analysis using the CLF5605 airfoil used in the Mars helicopter Ingenuity, which succeeded in its first flight on Mars. The calculation conditions are set to the Reynolds number (Re) at 75% rotor span in hovering (Re = 15,400), and the Mach number was varied from incompressible (M = 0.2) to transonic (M = 1.2). The compressible fluid dynamics solver FaSTAR developed by the Japan aerospace exploration agency (JAXA) is used, and calculations are performed under multiple conditions in which the Mach number and angle of attack (α) are swept. The results show that a flow field is similar to that in the Earth’s atmosphere above M = 1.0, such as bow shock at the leading edge, whereas multiple λ-type shock waves are observed over the separated shear layer above α = 3° at M = 0.80. However, no significant difference is found in the C<sub>p</sub> distribution around the airfoil between M = 0.6 and M = 0.8. From the results, it is found that multiple λ-type shock waves have no significant effect on the airfoil surface pressure distribution, the separated shear layer effect is dominant in the surface pressure change and aerodynamic characteristics.