In view of engineering application, it is practicable to decompose the aerodynamics into three components： the static aerodynamics, the aerodynamic increment due to steady rotations, and the aerodynamic increment due...In view of engineering application, it is practicable to decompose the aerodynamics into three components： the static aerodynamics, the aerodynamic increment due to steady rotations, and the aerodynamic increment due to unsteady separated and vortical flow. The first and the second components can be presented in conventional forms, while the third is described using a one-order differential equation and a radial-basis-function （RBF） network. For an aircraft configuration, the mathematical models of 6- component aerodynamic coefficients are set up from the wind tunnel test data of pitch, yaw, roll, and coupled yawroll large-amplitude oscillations. The flight dynamics of an aircraft is studied by the bifurcation analysis technique in the case of quasi-steady aerodynamics and unsteady aerodynam- ics, respectively. The results show that： （1） unsteady aerodynamics has no effect upon the existence of trim points, but affects their stability; （2） unsteady aerodynamics has great effects upon the existence, stability, and amplitudes of periodic solutions; and （3） unsteady aerodynamics changes the stable regions of trim points obviously. Furthermore, the dynamic responses of the aircraft to elevator deflections are inspected. It is shown that the unsteady aerodynamics is beneficial to dynamic stability for the present aircraft. Finally, the effects of unsteady aerodynamics on the post-stall maneuverability展开更多
This paper investigates the effect of the angle-of-attack(AOA)on the windward boundary-layer stability over a blunt cone with a nose radius of 5 mm.We consider a free-stream Mach number of 6 and a unit Reynolds number...This paper investigates the effect of the angle-of-attack(AOA)on the windward boundary-layer stability over a blunt cone with a nose radius of 5 mm.We consider a free-stream Mach number of 6 and a unit Reynolds number of 4.0×10^(7) m^(−1) and conduct both wind-tunnel experiments and stability analyses for AOAs ranging from 2°-10° at intervals of 2°.The results suggest that,as the AOA increases,the pressure gradient across the spanwise range becomes more pronounced,and the windward-side boundary layer becomes thinner.Using bi-global stability analysis,two unstable three-dimensional modes(varicose and sinuous)are identified in the windward boundary layer at various AOAs.The most unstable mode is the varicose V1 mode,in which the amplitude peak is initially close to the windward centerline and gradually shifts to the centerline downstream.Hence,the primary unstable disturbance exhibits a“V-like”distribution along the streamwise direction,which is likely to cause the V-shaped transition front observed in the wind-tunnel experiments.The eN method based on bi-global analysis is used to predict the transition location along the centerline on the windward region of the cone.The results indicate that,as the AOA increases,the transition location shifts forward,in line with our experimental results.Moreover,linear stability theory accurately predicts the eigenfunction and growth rate of the V1 mode obtained from bi-global analysis.This indicates that linear stability theory can be used to predict transitions in the windward boundary layer of a blunt cone at large AOAs under the conditions studied in this paper.展开更多
Abstract Accurate aerodynamic models are the basis of flight simulation and control law design. Mathematically modeling unsteady aerodynamics at high angles of attack bears great difficulties in model structure determ...Abstract Accurate aerodynamic models are the basis of flight simulation and control law design. Mathematically modeling unsteady aerodynamics at high angles of attack bears great difficulties in model structure determination and parameter estimation due to little understanding of the flow mechanism. Support vector machines （SVMs） based on statistical learning theory provide a novel tool for nonlinear system modeling. The work presented here examines the feasibility of applying SVMs to high angle.-of-attack unsteady aerodynamic modeling field. Mainly, after a review of SVMs, several issues associated with unsteady aerodynamic modeling by use of SVMs are discussed in detail, such as sele, ction of input variables, selection of output variables and determination of SVM parameters. The least squares SVM （LS-SVM） models are set up from certain dynamic wind tunnel test data of a delta wing and an aircraft configuration, and then used to predict the aerodynamic responses in other tests. The predictions are in good agreement with the test data, which indicates the satisfving learning and generalization performance of LS-SVMs.展开更多
Aerodynamic modeling and parameter estimation from quick accesses recorder （QAR） data is an important technical way to analyze the effects of highland weather conditions upon aerodynamic characteristics of airplane....Aerodynamic modeling and parameter estimation from quick accesses recorder （QAR） data is an important technical way to analyze the effects of highland weather conditions upon aerodynamic characteristics of airplane. It is also an essential content of flight accident analysis. The related techniques are developed in the present paper, including the geometric method for angle of attack and sideslip angle estimation, the extended Kalman filter associated with modified Bryson-Frazier smoother （EKF-MBF） method for aerodynamic coefficient identification, the radial basis function （RBF） neural network method for aerodynamic mod- eling, and the Delta method for stability/control derivative estimation. As an application example, the QAR data of a civil air- plane approaching a high-altitude airport are processed and the aerodynamic coefficient and derivative estimates are obtained. The estimation results are reasonable, which shows that the developed techniques are feasible. The causes for the distribution of aerodynamic derivative estimates are analyzed. Accordingly, several measures to improve estimation accuracy are put forward.展开更多
Mitigation of sonic boom to an acceptable stage is a key point for the next generation of supersonic transports. Meanwhile, designing a supersonic aircraft with an ideal ground signature is always the focus of researc...Mitigation of sonic boom to an acceptable stage is a key point for the next generation of supersonic transports. Meanwhile, designing a supersonic aircraft with an ideal ground signature is always the focus of research on sonic boom reduction. This paper presents an inverse design approach to optimize the near-field signature of an aircraft, making it close to the shaped ideal ground signature after the propagation in the atmosphere. Using the Proper Orthogonal Decomposition(POD) method, a guessed input of augmented Burgers equation is inversely achieved. By multiple POD iterations, the guessed ground signatures successively approach the target ground signature until the convergence criteria is reached. Finally, the corresponding equivalent area distribution is calculated from the optimal near-field signature through the classical Whitham F-function theory. To validate this method, an optimization example of Lockheed Martin 1021 is demonstrated. The modified configuration has a fully shaped ground signature and achieves a drop of perceived loudness by 7.94 PLdB. This improvement is achieved via shaping the original near-field signature into wiggles and damping it by atmospheric attenuation. At last, a nonphysical ground signature is set as the target to test the robustness of this inverse design method and shows that this method is robust enough for various inputs.展开更多
Abstract Based on the Reynolds-averaged Navier--Stokes （RANS） equations and structured grid technology, the calibration and validation of Y-Reo transition model is preformed with fifth-order weighted compact nonline...Abstract Based on the Reynolds-averaged Navier--Stokes （RANS） equations and structured grid technology, the calibration and validation of Y-Reo transition model is preformed with fifth-order weighted compact nonlinear scheme （WCNS）, and the purpose of the present work is to improve the numerical accuracy for aerodynamic characteristics simulation of low-speed flow with transition model on the basis of high-order numerical method study. Firstly, the empirical correlation functions involved in the Y-Reo transition model are modified and calibrated with experimental data of turbulent flat plates. Then, the grid convergence is studied on NLR-7301 two-element airfoil with the modified empirical correlation. At last, the modified empirical correlation is validated with NLR-7301 two-element airfoil and high-lift trapezoidal wing from transition location, velocity pro- file in boundary layer, surface pressure coefficient and aerodynamic characteristics. The numerical results illustrate that the numerical accuracy of transition length and skin friction behind transition location are improved with modified empirical correlation function, and obviously increases the numerical accuracy of aerodynamic characteristics prediction for typical transport configurations in low-speed range.展开更多
A CFD-based Numerical Virtual Flight(NVF)simulator is presented,which integrates an unsteady flow solver on moving hybrid grids,a Rigid-Body Dynamics(RBD)solver and a module of the Flight Control System(FCS).A techni...A CFD-based Numerical Virtual Flight(NVF)simulator is presented,which integrates an unsteady flow solver on moving hybrid grids,a Rigid-Body Dynamics(RBD)solver and a module of the Flight Control System(FCS).A technique of dynamic hybrid grids is developed to control the active control surfaces with body morphing,with a technique of parallel unstructured dynamic overlapping grids generating proper moving grids over the deflecting control surfaces(e.g.the afterbody rudders of a missile).For the flow/kinematic coupled problems,the 6 Degree-Of-Freedom(DOF)equations are solved by an explicit or implicit method coupled with the URANS CFD solver.The module of the control law is explicitly coupled into the NVF simulator and then improved by the simulation of the pitching maneuver process of a maneuverable missile model.A nonlinear dynamic inversion method is then implemented to design the control law for the pitching process of the maneuverable missile model.Simulations and analysis of the pitching maneuver process are carried out by the NVF simulator to improve the flight control law.Higher control response performance is obtained by adjusting the gain factors and adding an integrator into the control loop.展开更多
This paper numerically studies the aerodynamic performance of a bird-like bionic flapping wing.The geometry and kinematics are designed based on a seagull wing,in which flapping,folding,swaying,and twisting are consid...This paper numerically studies the aerodynamic performance of a bird-like bionic flapping wing.The geometry and kinematics are designed based on a seagull wing,in which flapping,folding,swaying,and twisting are considered.An in-house unsteady flow solver based on hybrid moving grids.is adopted for unsteady flow simulations.We focus on two main issues in this study,i.e.,the influence of the proportion of down-stroke and the effect of span-wise twisting.Numerical results show that the proportion of downstroke is closely related to the efficiency of the flapping process.The preferable proportion is about 0.7 by using the present geometry and kinematic model,which is very close to the observed data.Another finding is that the drag and the power consumption can be greatly reduced by the proper span-wise twisting.Two cases with different reduced frequencies are simulated and compared with each other.The numerical results show that the power consumption reduces by more than 20%,and the drag coefficient reduces by more than 60% through a proper twisting motion for both cases.The flow mechanism is mainly due to controlling of unsteady flow separation by adjusting the local effective angle of attack.These conclusions will be helpful for the high-performance micro air vehicle (MAV) design.展开更多
To understand fundamental problems in hypersonic laminar-turbulent boundary layer transition for three-dimensional complex vehicles,a new standard model with typical lifting-body features has been proposed,named as hy...To understand fundamental problems in hypersonic laminar-turbulent boundary layer transition for three-dimensional complex vehicles,a new standard model with typical lifting-body features has been proposed,named as hypersonic transition research vehicle(HyTRV).The configuration of HyTRV is fully analytical,and details of the design process are discussed in this study.The transition characteristics for HyTRV are investigated using three combined methods,i.e.,theoretical analyses,numerical simulations,and wind tunnel experiments.Results show that the fully analytic parameterization design of HyTRV can satisfy the model simplification requirements from both numerical simulations and wind tunnel experiments.Meanwhile,the flow field of HyTRV reveals typical transition mechanisms in six relatively separated regions,including the streamwise vortex instability,crossflow instability,secondary instability,and attachment-line instability.Therefore,the proposed HyTRV model is valuable for fundamental researches in hypersonic boundary layer transition.展开更多
It is common for aircraft to encounter atmospheric turbulence in flight tests.Turbulence is usually modeled as stochastic process noise in the flight dynamics equations.In this paper,parameter estimation of nonlinear ...It is common for aircraft to encounter atmospheric turbulence in flight tests.Turbulence is usually modeled as stochastic process noise in the flight dynamics equations.In this paper,parameter estimation of nonlinear dynamic system with both process and measurement noise was studied,and a practical filter error method was proposed.The linearized Kalman filter of first-order approximation was used for state estimation,in which the filter gain,along with the system parameters and the initial states,constituted the parameter vector to be estimated.The unknown parameters and measurement noise covariance were estimated alternately by a relaxation iteration method,and the sensitivities of observations to unknown parameters were calculated by finite difference approximation.Some practical aspects of the method application were discussed.The proposed filter error method was validated by the flight simulation data of a research aircraft.Then,the method was applied to the flight tests of a subscale aircraft,and the aerodynamic stability and control derivatives were estimated.All the estimation results were compared with the results of the output error method to demonstrate the effectiveness of the approach.It is shown that the filter error method is superior to the output error method for flight tests in atmospheric turbulence.展开更多
Experiment design method is a key to construct a highly reliable surrogate model for numerical optimization in large-scale project. Within the method, the experimental design criterion directly affects the accuracy of...Experiment design method is a key to construct a highly reliable surrogate model for numerical optimization in large-scale project. Within the method, the experimental design criterion directly affects the accuracy of the surrogate model and the optimization efficient. According to the shortcomings of the traditional experimental design, an improved adaptive sampling method is proposed in this paper. The surrogate model is firstly constructed by basic sparse samples. Then the supplementary sampling position is detected according to the specified criteria, which introduces the energy function and curvature sampling criteria based on radial basis function （RBF） network. Sampling detection criteria considers both the uniformity of sample distribution and the description of hypersurface curvature so as to significantly improve the prediction accuracy of the surrogate model with much less samples. For the surrogate model constructed with sparse samples, the sample uniformity is an important factor to the interpolation accuracy in the initial stage of adaptive sam- pling and surrogate model training. Along with the improvement of uniformity, the curvature description of objective function surface gradually becomes more important. In consideration of these issues, crowdness enhance function and root mean square error （RMSE） feedback function are introduced in C criterion expression. Thus, a new sampling method called RMSE and crowd- ness enhance （RCE） adaptive sampling is established. The validity of RCE adaptive sampling method is studied through typical test function firstly and then the airfoil/wing aerodynamic opti- mization design problem, which has high-dimensional design space. The results show that RCE adaptive sampling method not only reduces the requirement for the number of samples, but also effectively improves the prediction accuracy of the surrogate model, which has a broad prospects for applications.展开更多
The performance of an optimized aerodynamic shape is further improved by a second-step optimization using the design knowledge discovered by a data mining technique based on Proper Orthogonal Decomposition(POD) in the...The performance of an optimized aerodynamic shape is further improved by a second-step optimization using the design knowledge discovered by a data mining technique based on Proper Orthogonal Decomposition(POD) in the present study. Data generated in the first-step optimization by using evolution algorithms is saved as the source data, among which the superior data with improved objectives and maintained constraints is chosen. Only the geometry components of the superior data are picked out and used for constructing the snapshots of POD. Geometry characteristics of the superior data illustrated by POD bases are the design knowledge, by which the second-step optimization can be rapidly achieved. The optimization methods are demonstrated by redesigning a transonic compressor rotor blade, NASA Rotor 37, in the study to maximize the peak adiabatic efficiency, while maintaining the total pressure ratio and mass flow rate.Firstly, the blade is redesigned by using a particle swarm optimization method, and the adiabatic efficiency is increased by 1.29%. Then, the second-step optimization is performed by using the design knowledge, and a 0.25% gain on the adiabatic efficiency is obtained. The results are presented and addressed in detail, demonstrating that geometry variations significantly change the pattern and strength of the shock wave in the blade passage. The former reduces the separation loss,while the latter reduces the shock loss, and both favor an increase of the adiabatic efficiency.展开更多
An efficient MPI/OpenMP hybrid parallel Radial Basis Function (RBF) strategy for both continuous and discontinuous large-scale mesh deformation is proposed to reduce the computational cost and memory consumption.Unlik...An efficient MPI/OpenMP hybrid parallel Radial Basis Function (RBF) strategy for both continuous and discontinuous large-scale mesh deformation is proposed to reduce the computational cost and memory consumption.Unlike the conventional parallel methods in which all processors use the same surface displacement and implement the same operation,the present method employs different surface points sets and influence radius for each volume point movement,accompanied with efficient geometry searching strategy.The deformed surface points,also called Control Points (CPs),are stored in each processor.The displacement of spatial points is interpolated by using only 20-50 nearest control points,and the local influence radius is set to 5-20 times the maximum displacement of control points.To shorten the searching time for the nearest control point clouds,an Alternating Digital Tree (ADT) algorithm for 3D complex geometry is designed based on an iterative bisection technique.Besides,an MPI/OpenMP hybrid parallel approach is developed to reduce the memory cost in each High-Performance Computing (HPC) node for large-scale applications.Three 3D cases,including the ONERA-M6 wing and a commercial transport airplane standard model with up to 2.5 billion hybrid elements,are used to test the present mesh deformation method.The robustness and high parallel efficiency are demonstrated by a wing deflection case with a maximum bending angle of 450 and more than 80% parallel efficiency with 1024 MPI processors.In addition,the availability for both continuous and discontinuous surface deformation is verified by interpolating the projecting displacement with opposite directions surface points to the spatial points.展开更多
Circulation Control(CC) realizes rudderless flight control by driving compressed air jet to generate a virtual rudder surface, which significantly improves low detectability. The layout plan of combined control rudder...Circulation Control(CC) realizes rudderless flight control by driving compressed air jet to generate a virtual rudder surface, which significantly improves low detectability. The layout plan of combined control rudder surface is proposed based on the tailless flying wing aircraft. The closed-loop jet actuator system and stepless rudder surface switching control strategy are used to quantitatively study the control characteristics of circulation actuator for pitch and roll attitude through 3-DOF virtual flight test in a wind tunnel with a powered model at wind speed of 40 m/s. The results show that the combined use of circulation actuators can achieve bidirectional continuous and stable control of the aircraft’s pitch and roll attitude, with the maximum pitch rate of 12.3(°)/s and the maximum roll rate of 21.5(°)/s;the response time of attitude angular rate varying with the jet pressure ratio is less than 0.02 s, which can satisfy the control response requirements of aircraft motion stability for the control system;the jet rudder surface has a strong moment control ability, and the pitch moment of the jet elevator with a pressure ratio of 1.28 is the same as that of the mechanical elevator with 28° rudder deflection, which can expand the flight control boundary.展开更多
The interaction of an impinging oblique shock wave with an angle of 30°and a supersonic turbulent boundary layer at Ma_(∞)=2.9 and Re_(θ)=2400 over a wavy-wall is investigated through direct numerical simulatio...The interaction of an impinging oblique shock wave with an angle of 30°and a supersonic turbulent boundary layer at Ma_(∞)=2.9 and Re_(θ)=2400 over a wavy-wall is investigated through direct numerical simulation and compared with the interaction on a flat-plate under the same flow conditions.A sinusoidal wave with amplitude to wavelength ratio of 0.26 moves in the streamwise direction and is uniformly distributed across the spanwise direction.The influences of the wavy-wall on the interaction,including the characterization of the flow field,the skin-friction,pressure and the budget of turbulence kinetic energy,are systematically studied.The region of separation grows slightly and decomposes into four bubbles.Local peaks of skin-friction are observed at the rear part of the interaction region.The low-frequency shock motion can be seen in the wall pressure spectra.Analyses of the turbulence kinetic energy budget indicate that both diffusion and transport significantly increase near the crests,balanced by an amplified dissipation in the near-wall region.Proper orthogonal decomposition analyses show that the most energetic structures are associated with the separated shock and the shear layer over the bubbles.Only the bubbles in the first two troughs are dominated by a low-frequency enlargement or shrinkage.展开更多
Flying condition with angle of attack is inevitable in a hypersonic vehicle,and it may influence the thermal protection system(TPS)performance of opposing jet and its combinations.A 3D Navier–Stokes equation and shea...Flying condition with angle of attack is inevitable in a hypersonic vehicle,and it may influence the thermal protection system(TPS)performance of opposing jet and its combinations.A 3D Navier–Stokes equation and shear stress transfer(SST)k-ωmodel with compressible correction are employed to simulate the angle of attack characteristics of a blunt body with opposing jet and platelet transpiration TPS.The flowfield and heat flux transfer for angles of attack 0°,3°,6°with jet pressure ratio PR=0.1 and 0°,6°,12°with PR=0.2 are obtained.Numerical results show that the flowfield is no longer symmetrical with the effect of the angle of attack.The flowfield and heat transfer in windward and leeward performed adversely.The recompression shock wave in windward is strengthened,which increases local temperature and strengthens heat transfer.The opposing jet fails in thermal protection when the angle of attack reaches critical value;however,the critical angle of attack can be promoted by increasing PR.Finally,the transpiration gas can strengthen the cooling efficiency of windward,thereby,increasing the critical angle of attack.展开更多
The reattached boundary layer in the interaction of an oblique shock wave with a flatplate turbulent boundary layer at Mach number 2.25 is studied by means of Direct Numerical Simulation(DNS).The numerical results are...The reattached boundary layer in the interaction of an oblique shock wave with a flatplate turbulent boundary layer at Mach number 2.25 is studied by means of Direct Numerical Simulation(DNS).The numerical results are carefully compared with available experimental and DNS data in terms of turbulence statistics,wall pressure and skin friction.The coherent vortex structures are significantly enhanced due to the shock interaction,and the reattached boundary layer is characterized by large-scale structures in the outer region.The space-time correlation of fluctuating wall shear stress and streamwise velocity fluctuation reveals that the structural inclination angle exhibits a gradual decrease during the recovery process.The scale interactions are analyzed by using a twopoint amplitude modulation correlation.A possible mechanism is proposed to account for the strong amplitude modulation in the downstream region.Moreover,the mean skin-friction is decomposed to understand the physically informed contributions.Unlike the upstream Turbulent Boundary Layer(TBL),the contribution associated with the Turbulence Kinetic Energy(TKE)production is greatly amplified,while the spatial growth contribution induced by the pressure gradient largely inhibits skin-friction generation.Based on bidimensional empirical mode decomposition,the turbulence kinetic energy production contribution is further split into different terms with specific spanwise length scales.展开更多
The simulation of hypersonic flows with fully unstructured(tetrahedral)grids has severe problems with respect to the prediction of stagnation region heating,due to the random face orientation without alignment to the ...The simulation of hypersonic flows with fully unstructured(tetrahedral)grids has severe problems with respect to the prediction of stagnation region heating,due to the random face orientation without alignment to the bow shock.To improve the accuracy of aero-heating predictions,three multi-dimensional approaches on unstructured grids are coupled in our Reynolds-averaged Navier-Stokes(RANS)solver,including multi-dimensional upwind flux reconstruction(MUP),multi-dimensional limiter(MLP-u2)and multi-dimensional gradient reconstruction(MLR).The coupled multi-dimensional RANS solver is validated by several typical verification and validation(V&V)cases,including hypersonic flows over a cylinder,a blunt biconic,and a double-ellipsoid,with commonly used prism/tetrahedral hybrid grids.Finally,the coupled multi-dimensional solver is applied to simulating the heat flux distribution over a 3D engineering configuration,i.e.a Hermes-like space shuttle model.The obtained numerical results are compared with experimental data.The predicted results demonstrate that the coupled multi-dimensional approach has a good prediction capability on aerodynamic heating over a wide range of complex engineering configurations.展开更多
Based on improved multi-objective particle swarm optimization（MOPSO） algorithm with principal component analysis（PCA） methodology, an efficient high-dimension multiobjective optimization method is proposed, which,...Based on improved multi-objective particle swarm optimization（MOPSO） algorithm with principal component analysis（PCA） methodology, an efficient high-dimension multiobjective optimization method is proposed, which, as the purpose of this paper, aims to improve the convergence of Pareto front in multi-objective optimization design. The mathematical efficiency,the physical reasonableness and the reliability in dealing with redundant objectives of PCA are verified by typical DTLZ5 test function and multi-objective correlation analysis of supercritical airfoil,and the proposed method is integrated into aircraft multi-disciplinary design（AMDEsign） platform, which contains aerodynamics, stealth and structure weight analysis and optimization module.Then the proposed method is used for the multi-point integrated aerodynamic optimization of a wide-body passenger aircraft, in which the redundant objectives identified by PCA are transformed to optimization constraints, and several design methods are compared. The design results illustrate that the strategy used in this paper is sufficient and multi-point design requirements of the passenger aircraft are reached. The visualization level of non-dominant Pareto set is improved by effectively reducing the dimension without losing the primary feature of the problem.展开更多
Well-designed airfoil is very important for high-performance rotor.This paper developed an efficient multi-objective and multi-constraint optimization design system for rotor airfoils based on RANS analysis,and verifi...Well-designed airfoil is very important for high-performance rotor.This paper developed an efficient multi-objective and multi-constraint optimization design system for rotor airfoils based on RANS analysis,and verified the performance of the optimized airfoil.Using CRA09-A as the baseline rotor airfoil,the CRA09-B optimized rotor airfoil was designed successfully.Combined with the foundation of high-precision rotor airfoil stationary test technology,the CRA09-B and CRA09-A rotor airfoils were tested in the S3 MA high-speed wind tunnel of ONERA.In order to correct the aerodynamic data,a single parameter linear wall pressure method is used to consider the tunnel effects.The results indicate that multi-objective and multi-constraint optimization design method developed in this study is reliable,and that CRA09-B optimized airfoil provides better stationary performance than CRA09-A airfoil in terms of maximum lift coefficient and lift-to-drag ratio.展开更多
文摘In view of engineering application, it is practicable to decompose the aerodynamics into three components： the static aerodynamics, the aerodynamic increment due to steady rotations, and the aerodynamic increment due to unsteady separated and vortical flow. The first and the second components can be presented in conventional forms, while the third is described using a one-order differential equation and a radial-basis-function （RBF） network. For an aircraft configuration, the mathematical models of 6- component aerodynamic coefficients are set up from the wind tunnel test data of pitch, yaw, roll, and coupled yawroll large-amplitude oscillations. The flight dynamics of an aircraft is studied by the bifurcation analysis technique in the case of quasi-steady aerodynamics and unsteady aerodynam- ics, respectively. The results show that： （1） unsteady aerodynamics has no effect upon the existence of trim points, but affects their stability; （2） unsteady aerodynamics has great effects upon the existence, stability, and amplitudes of periodic solutions; and （3） unsteady aerodynamics changes the stable regions of trim points obviously. Furthermore, the dynamic responses of the aircraft to elevator deflections are inspected. It is shown that the unsteady aerodynamics is beneficial to dynamic stability for the present aircraft. Finally, the effects of unsteady aerodynamics on the post-stall maneuverability
基金supported by the National Natural Science Foundation of China(Grant No.92052301)the National Numerical Wind-tunnel Project(NNW).
文摘This paper investigates the effect of the angle-of-attack(AOA)on the windward boundary-layer stability over a blunt cone with a nose radius of 5 mm.We consider a free-stream Mach number of 6 and a unit Reynolds number of 4.0×10^(7) m^(−1) and conduct both wind-tunnel experiments and stability analyses for AOAs ranging from 2°-10° at intervals of 2°.The results suggest that,as the AOA increases,the pressure gradient across the spanwise range becomes more pronounced,and the windward-side boundary layer becomes thinner.Using bi-global stability analysis,two unstable three-dimensional modes(varicose and sinuous)are identified in the windward boundary layer at various AOAs.The most unstable mode is the varicose V1 mode,in which the amplitude peak is initially close to the windward centerline and gradually shifts to the centerline downstream.Hence,the primary unstable disturbance exhibits a“V-like”distribution along the streamwise direction,which is likely to cause the V-shaped transition front observed in the wind-tunnel experiments.The eN method based on bi-global analysis is used to predict the transition location along the centerline on the windward region of the cone.The results indicate that,as the AOA increases,the transition location shifts forward,in line with our experimental results.Moreover,linear stability theory accurately predicts the eigenfunction and growth rate of the V1 mode obtained from bi-global analysis.This indicates that linear stability theory can be used to predict transitions in the windward boundary layer of a blunt cone at large AOAs under the conditions studied in this paper.
文摘Abstract Accurate aerodynamic models are the basis of flight simulation and control law design. Mathematically modeling unsteady aerodynamics at high angles of attack bears great difficulties in model structure determination and parameter estimation due to little understanding of the flow mechanism. Support vector machines （SVMs） based on statistical learning theory provide a novel tool for nonlinear system modeling. The work presented here examines the feasibility of applying SVMs to high angle.-of-attack unsteady aerodynamic modeling field. Mainly, after a review of SVMs, several issues associated with unsteady aerodynamic modeling by use of SVMs are discussed in detail, such as sele, ction of input variables, selection of output variables and determination of SVM parameters. The least squares SVM （LS-SVM） models are set up from certain dynamic wind tunnel test data of a delta wing and an aircraft configuration, and then used to predict the aerodynamic responses in other tests. The predictions are in good agreement with the test data, which indicates the satisfving learning and generalization performance of LS-SVMs.
基金National Natural Science Foundation of China(60832012)
文摘Aerodynamic modeling and parameter estimation from quick accesses recorder （QAR） data is an important technical way to analyze the effects of highland weather conditions upon aerodynamic characteristics of airplane. It is also an essential content of flight accident analysis. The related techniques are developed in the present paper, including the geometric method for angle of attack and sideslip angle estimation, the extended Kalman filter associated with modified Bryson-Frazier smoother （EKF-MBF） method for aerodynamic coefficient identification, the radial basis function （RBF） neural network method for aerodynamic mod- eling, and the Delta method for stability/control derivative estimation. As an application example, the QAR data of a civil air- plane approaching a high-altitude airport are processed and the aerodynamic coefficient and derivative estimates are obtained. The estimation results are reasonable, which shows that the developed techniques are feasible. The causes for the distribution of aerodynamic derivative estimates are analyzed. Accordingly, several measures to improve estimation accuracy are put forward.
文摘Mitigation of sonic boom to an acceptable stage is a key point for the next generation of supersonic transports. Meanwhile, designing a supersonic aircraft with an ideal ground signature is always the focus of research on sonic boom reduction. This paper presents an inverse design approach to optimize the near-field signature of an aircraft, making it close to the shaped ideal ground signature after the propagation in the atmosphere. Using the Proper Orthogonal Decomposition(POD) method, a guessed input of augmented Burgers equation is inversely achieved. By multiple POD iterations, the guessed ground signatures successively approach the target ground signature until the convergence criteria is reached. Finally, the corresponding equivalent area distribution is calculated from the optimal near-field signature through the classical Whitham F-function theory. To validate this method, an optimization example of Lockheed Martin 1021 is demonstrated. The modified configuration has a fully shaped ground signature and achieves a drop of perceived loudness by 7.94 PLdB. This improvement is achieved via shaping the original near-field signature into wiggles and damping it by atmospheric attenuation. At last, a nonphysical ground signature is set as the target to test the robustness of this inverse design method and shows that this method is robust enough for various inputs.
基金supported by the National Basic Research Program of China(No.2014CB744803)
文摘Abstract Based on the Reynolds-averaged Navier--Stokes （RANS） equations and structured grid technology, the calibration and validation of Y-Reo transition model is preformed with fifth-order weighted compact nonlinear scheme （WCNS）, and the purpose of the present work is to improve the numerical accuracy for aerodynamic characteristics simulation of low-speed flow with transition model on the basis of high-order numerical method study. Firstly, the empirical correlation functions involved in the Y-Reo transition model are modified and calibrated with experimental data of turbulent flat plates. Then, the grid convergence is studied on NLR-7301 two-element airfoil with the modified empirical correlation. At last, the modified empirical correlation is validated with NLR-7301 two-element airfoil and high-lift trapezoidal wing from transition location, velocity pro- file in boundary layer, surface pressure coefficient and aerodynamic characteristics. The numerical results illustrate that the numerical accuracy of transition length and skin friction behind transition location are improved with modified empirical correlation function, and obviously increases the numerical accuracy of aerodynamic characteristics prediction for typical transport configurations in low-speed range.
基金supported partially by National Key Research and Development Program (No. 2016YFB0200701)National Natural Science Foundation of China (Nos. 11532016 and 11672324)
文摘A CFD-based Numerical Virtual Flight(NVF)simulator is presented,which integrates an unsteady flow solver on moving hybrid grids,a Rigid-Body Dynamics(RBD)solver and a module of the Flight Control System(FCS).A technique of dynamic hybrid grids is developed to control the active control surfaces with body morphing,with a technique of parallel unstructured dynamic overlapping grids generating proper moving grids over the deflecting control surfaces(e.g.the afterbody rudders of a missile).For the flow/kinematic coupled problems,the 6 Degree-Of-Freedom(DOF)equations are solved by an explicit or implicit method coupled with the URANS CFD solver.The module of the control law is explicitly coupled into the NVF simulator and then improved by the simulation of the pitching maneuver process of a maneuverable missile model.A nonlinear dynamic inversion method is then implemented to design the control law for the pitching process of the maneuverable missile model.Simulations and analysis of the pitching maneuver process are carried out by the NVF simulator to improve the flight control law.Higher control response performance is obtained by adjusting the gain factors and adding an integrator into the control loop.
基金Project supported by the National Key Research and Development Program(No.2016YFB0200700)the National Natural Science Foundation of China(Nos.11532016 and 11672324)
文摘This paper numerically studies the aerodynamic performance of a bird-like bionic flapping wing.The geometry and kinematics are designed based on a seagull wing,in which flapping,folding,swaying,and twisting are considered.An in-house unsteady flow solver based on hybrid moving grids.is adopted for unsteady flow simulations.We focus on two main issues in this study,i.e.,the influence of the proportion of down-stroke and the effect of span-wise twisting.Numerical results show that the proportion of downstroke is closely related to the efficiency of the flapping process.The preferable proportion is about 0.7 by using the present geometry and kinematic model,which is very close to the observed data.Another finding is that the drag and the power consumption can be greatly reduced by the proper span-wise twisting.Two cases with different reduced frequencies are simulated and compared with each other.The numerical results show that the power consumption reduces by more than 20%,and the drag coefficient reduces by more than 60% through a proper twisting motion for both cases.The flow mechanism is mainly due to controlling of unsteady flow separation by adjusting the local effective angle of attack.These conclusions will be helpful for the high-performance micro air vehicle (MAV) design.
基金This work was supported by the National Natural Science Foundation of China(Grant 11702315,92052301)the National Key Research and Development Program of China(Grant 2016YFA0401200).
文摘To understand fundamental problems in hypersonic laminar-turbulent boundary layer transition for three-dimensional complex vehicles,a new standard model with typical lifting-body features has been proposed,named as hypersonic transition research vehicle(HyTRV).The configuration of HyTRV is fully analytical,and details of the design process are discussed in this study.The transition characteristics for HyTRV are investigated using three combined methods,i.e.,theoretical analyses,numerical simulations,and wind tunnel experiments.Results show that the fully analytic parameterization design of HyTRV can satisfy the model simplification requirements from both numerical simulations and wind tunnel experiments.Meanwhile,the flow field of HyTRV reveals typical transition mechanisms in six relatively separated regions,including the streamwise vortex instability,crossflow instability,secondary instability,and attachment-line instability.Therefore,the proposed HyTRV model is valuable for fundamental researches in hypersonic boundary layer transition.
基金supported by the National Natural Science Foundation of China(No.11802325)。
文摘It is common for aircraft to encounter atmospheric turbulence in flight tests.Turbulence is usually modeled as stochastic process noise in the flight dynamics equations.In this paper,parameter estimation of nonlinear dynamic system with both process and measurement noise was studied,and a practical filter error method was proposed.The linearized Kalman filter of first-order approximation was used for state estimation,in which the filter gain,along with the system parameters and the initial states,constituted the parameter vector to be estimated.The unknown parameters and measurement noise covariance were estimated alternately by a relaxation iteration method,and the sensitivities of observations to unknown parameters were calculated by finite difference approximation.Some practical aspects of the method application were discussed.The proposed filter error method was validated by the flight simulation data of a research aircraft.Then,the method was applied to the flight tests of a subscale aircraft,and the aerodynamic stability and control derivatives were estimated.All the estimation results were compared with the results of the output error method to demonstrate the effectiveness of the approach.It is shown that the filter error method is superior to the output error method for flight tests in atmospheric turbulence.
基金co-supported by the National Natural Science Foundation of China (Nos. 11402288 and 11372254)
文摘Experiment design method is a key to construct a highly reliable surrogate model for numerical optimization in large-scale project. Within the method, the experimental design criterion directly affects the accuracy of the surrogate model and the optimization efficient. According to the shortcomings of the traditional experimental design, an improved adaptive sampling method is proposed in this paper. The surrogate model is firstly constructed by basic sparse samples. Then the supplementary sampling position is detected according to the specified criteria, which introduces the energy function and curvature sampling criteria based on radial basis function （RBF） network. Sampling detection criteria considers both the uniformity of sample distribution and the description of hypersurface curvature so as to significantly improve the prediction accuracy of the surrogate model with much less samples. For the surrogate model constructed with sparse samples, the sample uniformity is an important factor to the interpolation accuracy in the initial stage of adaptive sam- pling and surrogate model training. Along with the improvement of uniformity, the curvature description of objective function surface gradually becomes more important. In consideration of these issues, crowdness enhance function and root mean square error （RMSE） feedback function are introduced in C criterion expression. Thus, a new sampling method called RMSE and crowd- ness enhance （RCE） adaptive sampling is established. The validity of RCE adaptive sampling method is studied through typical test function firstly and then the airfoil/wing aerodynamic opti- mization design problem, which has high-dimensional design space. The results show that RCE adaptive sampling method not only reduces the requirement for the number of samples, but also effectively improves the prediction accuracy of the surrogate model, which has a broad prospects for applications.
基金supported by the National Natural Science Foundation of China(Nos.51676003,51206003,and 11702305)
文摘The performance of an optimized aerodynamic shape is further improved by a second-step optimization using the design knowledge discovered by a data mining technique based on Proper Orthogonal Decomposition(POD) in the present study. Data generated in the first-step optimization by using evolution algorithms is saved as the source data, among which the superior data with improved objectives and maintained constraints is chosen. Only the geometry components of the superior data are picked out and used for constructing the snapshots of POD. Geometry characteristics of the superior data illustrated by POD bases are the design knowledge, by which the second-step optimization can be rapidly achieved. The optimization methods are demonstrated by redesigning a transonic compressor rotor blade, NASA Rotor 37, in the study to maximize the peak adiabatic efficiency, while maintaining the total pressure ratio and mass flow rate.Firstly, the blade is redesigned by using a particle swarm optimization method, and the adiabatic efficiency is increased by 1.29%. Then, the second-step optimization is performed by using the design knowledge, and a 0.25% gain on the adiabatic efficiency is obtained. The results are presented and addressed in detail, demonstrating that geometry variations significantly change the pattern and strength of the shock wave in the blade passage. The former reduces the separation loss,while the latter reduces the shock loss, and both favor an increase of the adiabatic efficiency.
基金supported by the National Key Research and Development Program of China (No.2016YFB0200701)the National Natural Science Foundation of China (Nos. 11532016 and 91530325)
文摘An efficient MPI/OpenMP hybrid parallel Radial Basis Function (RBF) strategy for both continuous and discontinuous large-scale mesh deformation is proposed to reduce the computational cost and memory consumption.Unlike the conventional parallel methods in which all processors use the same surface displacement and implement the same operation,the present method employs different surface points sets and influence radius for each volume point movement,accompanied with efficient geometry searching strategy.The deformed surface points,also called Control Points (CPs),are stored in each processor.The displacement of spatial points is interpolated by using only 20-50 nearest control points,and the local influence radius is set to 5-20 times the maximum displacement of control points.To shorten the searching time for the nearest control point clouds,an Alternating Digital Tree (ADT) algorithm for 3D complex geometry is designed based on an iterative bisection technique.Besides,an MPI/OpenMP hybrid parallel approach is developed to reduce the memory cost in each High-Performance Computing (HPC) node for large-scale applications.Three 3D cases,including the ONERA-M6 wing and a commercial transport airplane standard model with up to 2.5 billion hybrid elements,are used to test the present mesh deformation method.The robustness and high parallel efficiency are demonstrated by a wing deflection case with a maximum bending angle of 450 and more than 80% parallel efficiency with 1024 MPI processors.In addition,the availability for both continuous and discontinuous surface deformation is verified by interpolating the projecting displacement with opposite directions surface points to the spatial points.
基金supported by the Equipment Pre-research Common Technology Project,China(No.41406010101).
文摘Circulation Control(CC) realizes rudderless flight control by driving compressed air jet to generate a virtual rudder surface, which significantly improves low detectability. The layout plan of combined control rudder surface is proposed based on the tailless flying wing aircraft. The closed-loop jet actuator system and stepless rudder surface switching control strategy are used to quantitatively study the control characteristics of circulation actuator for pitch and roll attitude through 3-DOF virtual flight test in a wind tunnel with a powered model at wind speed of 40 m/s. The results show that the combined use of circulation actuators can achieve bidirectional continuous and stable control of the aircraft’s pitch and roll attitude, with the maximum pitch rate of 12.3(°)/s and the maximum roll rate of 21.5(°)/s;the response time of attitude angular rate varying with the jet pressure ratio is less than 0.02 s, which can satisfy the control response requirements of aircraft motion stability for the control system;the jet rudder surface has a strong moment control ability, and the pitch moment of the jet elevator with a pressure ratio of 1.28 is the same as that of the mechanical elevator with 28° rudder deflection, which can expand the flight control boundary.
基金co-supported by the National Key Research and Development Program of China(Nos.2019YFA0405300 and 2016YFA0401200)the National Natural Science Foundation of China(Nos.11972356 and 91852203)+1 种基金LHD Youth Innovation Fund(No.LHD2019CX04)National Numerical Wind Tunnel Project。
文摘The interaction of an impinging oblique shock wave with an angle of 30°and a supersonic turbulent boundary layer at Ma_(∞)=2.9 and Re_(θ)=2400 over a wavy-wall is investigated through direct numerical simulation and compared with the interaction on a flat-plate under the same flow conditions.A sinusoidal wave with amplitude to wavelength ratio of 0.26 moves in the streamwise direction and is uniformly distributed across the spanwise direction.The influences of the wavy-wall on the interaction,including the characterization of the flow field,the skin-friction,pressure and the budget of turbulence kinetic energy,are systematically studied.The region of separation grows slightly and decomposes into four bubbles.Local peaks of skin-friction are observed at the rear part of the interaction region.The low-frequency shock motion can be seen in the wall pressure spectra.Analyses of the turbulence kinetic energy budget indicate that both diffusion and transport significantly increase near the crests,balanced by an amplified dissipation in the near-wall region.Proper orthogonal decomposition analyses show that the most energetic structures are associated with the separated shock and the shear layer over the bubbles.Only the bubbles in the first two troughs are dominated by a low-frequency enlargement or shrinkage.
文摘Flying condition with angle of attack is inevitable in a hypersonic vehicle,and it may influence the thermal protection system(TPS)performance of opposing jet and its combinations.A 3D Navier–Stokes equation and shear stress transfer(SST)k-ωmodel with compressible correction are employed to simulate the angle of attack characteristics of a blunt body with opposing jet and platelet transpiration TPS.The flowfield and heat flux transfer for angles of attack 0°,3°,6°with jet pressure ratio PR=0.1 and 0°,6°,12°with PR=0.2 are obtained.Numerical results show that the flowfield is no longer symmetrical with the effect of the angle of attack.The flowfield and heat transfer in windward and leeward performed adversely.The recompression shock wave in windward is strengthened,which increases local temperature and strengthens heat transfer.The opposing jet fails in thermal protection when the angle of attack reaches critical value;however,the critical angle of attack can be promoted by increasing PR.Finally,the transpiration gas can strengthen the cooling efficiency of windward,thereby,increasing the critical angle of attack.
基金co-supported by the National Key R&D Program of China (No. 2019YFA0405300)the National Natural Science Foundation of China (Nos. 11972356 and 91852203)
文摘The reattached boundary layer in the interaction of an oblique shock wave with a flatplate turbulent boundary layer at Mach number 2.25 is studied by means of Direct Numerical Simulation(DNS).The numerical results are carefully compared with available experimental and DNS data in terms of turbulence statistics,wall pressure and skin friction.The coherent vortex structures are significantly enhanced due to the shock interaction,and the reattached boundary layer is characterized by large-scale structures in the outer region.The space-time correlation of fluctuating wall shear stress and streamwise velocity fluctuation reveals that the structural inclination angle exhibits a gradual decrease during the recovery process.The scale interactions are analyzed by using a twopoint amplitude modulation correlation.A possible mechanism is proposed to account for the strong amplitude modulation in the downstream region.Moreover,the mean skin-friction is decomposed to understand the physically informed contributions.Unlike the upstream Turbulent Boundary Layer(TBL),the contribution associated with the Turbulence Kinetic Energy(TKE)production is greatly amplified,while the spatial growth contribution induced by the pressure gradient largely inhibits skin-friction generation.Based on bidimensional empirical mode decomposition,the turbulence kinetic energy production contribution is further split into different terms with specific spanwise length scales.
基金the National Key Research&Development Program of China(2016YFB020071)the National Natural Science Foundation of China(Grants 11532016 and 11702315).
文摘The simulation of hypersonic flows with fully unstructured(tetrahedral)grids has severe problems with respect to the prediction of stagnation region heating,due to the random face orientation without alignment to the bow shock.To improve the accuracy of aero-heating predictions,three multi-dimensional approaches on unstructured grids are coupled in our Reynolds-averaged Navier-Stokes(RANS)solver,including multi-dimensional upwind flux reconstruction(MUP),multi-dimensional limiter(MLP-u2)and multi-dimensional gradient reconstruction(MLR).The coupled multi-dimensional RANS solver is validated by several typical verification and validation(V&V)cases,including hypersonic flows over a cylinder,a blunt biconic,and a double-ellipsoid,with commonly used prism/tetrahedral hybrid grids.Finally,the coupled multi-dimensional solver is applied to simulating the heat flux distribution over a 3D engineering configuration,i.e.a Hermes-like space shuttle model.The obtained numerical results are compared with experimental data.The predicted results demonstrate that the coupled multi-dimensional approach has a good prediction capability on aerodynamic heating over a wide range of complex engineering configurations.
基金supported by the National Natural Science Foundation of China (No.11402288)
文摘Based on improved multi-objective particle swarm optimization（MOPSO） algorithm with principal component analysis（PCA） methodology, an efficient high-dimension multiobjective optimization method is proposed, which, as the purpose of this paper, aims to improve the convergence of Pareto front in multi-objective optimization design. The mathematical efficiency,the physical reasonableness and the reliability in dealing with redundant objectives of PCA are verified by typical DTLZ5 test function and multi-objective correlation analysis of supercritical airfoil,and the proposed method is integrated into aircraft multi-disciplinary design（AMDEsign） platform, which contains aerodynamics, stealth and structure weight analysis and optimization module.Then the proposed method is used for the multi-point integrated aerodynamic optimization of a wide-body passenger aircraft, in which the redundant objectives identified by PCA are transformed to optimization constraints, and several design methods are compared. The design results illustrate that the strategy used in this paper is sufficient and multi-point design requirements of the passenger aircraft are reached. The visualization level of non-dominant Pareto set is improved by effectively reducing the dimension without losing the primary feature of the problem.
基金supported by the National Natural Science Foundation of China(No.11902335)。
文摘Well-designed airfoil is very important for high-performance rotor.This paper developed an efficient multi-objective and multi-constraint optimization design system for rotor airfoils based on RANS analysis,and verified the performance of the optimized airfoil.Using CRA09-A as the baseline rotor airfoil,the CRA09-B optimized rotor airfoil was designed successfully.Combined with the foundation of high-precision rotor airfoil stationary test technology,the CRA09-B and CRA09-A rotor airfoils were tested in the S3 MA high-speed wind tunnel of ONERA.In order to correct the aerodynamic data,a single parameter linear wall pressure method is used to consider the tunnel effects.The results indicate that multi-objective and multi-constraint optimization design method developed in this study is reliable,and that CRA09-B optimized airfoil provides better stationary performance than CRA09-A airfoil in terms of maximum lift coefficient and lift-to-drag ratio.
