Model of Airflow Field on the Deck for Shipborne Helicopter Flight Dynamics Analysis

For the research of helicopter/ship dynamic interface,the method of combining steady flow and stochastic flow is adopted to establish a flow field model applied to the flight dynamics analysis of shipborne helicopter.The steady flow is calculated by computational fluid dynamics(CFD)method,while the stochastic flow is composed of the compensation velocity derived from ship motion and turbulence above the deck.The accuracy of the proposed flow field model is verified by comparing the helicopter response in the proposed flow field with the results calculated by US Army′s Military Specifications(MIL)model which is commonly used in engineering.Meanwhile,it also shows the proposed flow field model is more appliable to flight dynamics analysis of shipborne helicopter.On that the basis,ship deck flow field is simulated at different sea conditions by adjusting the wind speed on the deck,mother ship movement and shipboard turbulence,etc.And helicopter angular rate response is calculated.The results show that the difference of dynamic stability between helicopter′s roll and pitch leads to the facts that the influence of above factors on the helicopter′s roll angular rate response is greater than that of pitch angular rate,that the frequency and amplitude of mother ship roll motion are much greater than those of pitch motion,and that the disturbance caused by roll motion on the air has greater influence on the helicopter response.The shipboard turbulence is the main disturbance factor that influences helicopter flight stability and its intensity determines the amplitudes of angular rate response.     

DEVELOPMENT OF HYBRID RANS / LES CODE AND ITS VALIDATION

A computational code based on the hybrid RANSLES approach is developed. The hybrid approach com- bines the delayed detached-eddy simulation（DDES） with an improved RANS-LES hybrid model aiming at wall modeling in LES（WMLES）. In the code, the convective flux is solved using the fourth-order skew-symmetric scheme so as to diminish the negative effect of numerical dissipation. The Spalart-Allmaras（S-A） model is applied as a suhgrid scale（SGS） model. To validate the developed code, homogeneous isotropic turbulence and turbulent channel flow are simulated and the results are compared with experimental data and DNS results. The results of the isotropic turbulence show that the fourth-order skew-symmetric scheme is adequate enough and the model works well coupling with the convective scheme. The results of the turbulent channel flow agree well with the DNS data, the predicted velocity profiles at Reynolds number from 178 to 2 700 match well with the Reichardt＇s law, and the organized vortical structures are well captured.     

Drag Reducing and Increasing Mechanism on Triangular Riblet Surface

Drag reducing and increasing mechanism on riblet surface has been studied through computational fluid dynamics(CFD).Drag reduction is achieved through the optimization of riblet geometry which would affect flow structure inside riblet grooves.Force and flow structure on riblet surface are analyzed and compared with those of smooth surface based on the k-εturbulence model.Drag reducing and increasing mechanism is proved to be related to microvortexes induced inside riblets which lead to Reynolds shear stress reduction significantly and is considered to be the dominant factor resulting in wall friction reduction.Simulation results also show that the pressure drag generating from the deviation of static pressure on the front and rear ends of riblets occurs and grows exponentially with Mach number,which can cause drag increasing.Furthermore,near-wall vortical structures,Reynolds shear stress and static pressure on riblet surfaces are also analyzed in detail.     

Influence of Different Velocities on Muzzle Flow Field

A two?dimensional axisymmetric numerical simulation was successfully carried out on the muzzle flow field of a 300 mm?caliber counter?mass propelling gun. Based on the FLUENT software,using the finite volume method(FVM)and the realizable k?ε turbulence model,we adopted the holistic movement of a partitioned mesh processing method coupled with the intermediate ballistic model and the six degree?of?freedom model(6?DOF). We compared the flow field characteristics at the velocity of 1 730.4,978.3,and 323.4 m/s. The results indicate that the pressure of the hypersonic initial flow field is much higher than that of the subsonic and supersonic initial flow fields. In the case of the subsonic(323.4 m/s)flow field,the tiny disturbance spreads throughout the whole domain. But in the cases of the supersonic(978.3 m/s) and the hypersonic(1 730.4 m/s) flow fields,it cannot spread to the upstream disturbance source,and the disturbance domain of the supersonic flow field is wider than that of the hypersonic. It is noted that the subsonic flow field has a rounded shock wave before the projectile. However,in the supersonic and hypersonic flow fields,a shear layer is formed which begins from the head of the projectile and extends outward from the side of the projectile. Then a multi?layer shock wave is formed composed of coronal shock waves,bottom shock waves,reflected shock waves,and Mach disk.     

HIGH TEMPERATURE RHEOLOGICAL BEHAVIOR OF FSJ JOINTED REGION FOR 7022ALUMINUM ALLOY

The Hopkinson pressure bar tests for base metal and friction stir jointing（FSJ） jointed region of 7022 a- luminum alloy are carried out at different temperatures and strain rates. The temperature is 30--400 ℃ and the strain rate is 1 200-5 000 s 1. High strain rate for base metal and FSJ jointed region of 7022 aluminum alloy are studied. The corresponding stress-strain curves are obtained. The results show that the flow stresses of base metal and FSJ jointed region of 7022 aluminum alloy decline with the increase of temperature and increase with the in- crease of strain rate. Furthermore, the constitutive equation for base metal and FSJ jointed region of 7022 alumi- num alloy at high temperature and high strain rate is obtained based on Johnson-Cook model.     

SIMULATION OF BLEED AIR BEHAVIOR DURING AIRCRAFT IN FLIGHT BASED ON FLOWMASTER

Bleed air system is one of the most important components of air management system(AMS).It acts as transfer pipes responsible for air supply at high temperature and pressure.The thermal and flow performance of the bleed air system is a key issue for the design of AMS since the characteristics of air source have a great influence on the anti-ice system,the environmental control system and other downstream system in need of high temperature pressurized air.Based on the one-dimensional lumped parameter technology,a computer analysis model of bleed air system is developed in order to analyze the thermal and flow behaviors of the nodal points in the pipeline network.The simulation are performed with a given flight assignment using the analysis model,and the results verify that the system meets the design requirements.     

Effects of Thermo-Mechanical Loads on Aeroelastic Instabilities of Metallic and Composite Panels

Panel flutter phenomena can be strongly affected by thermal loads,and so a refined aeroelastic model is presented.Higher-order shell theories are used as structural models.The aerodynamic forces are described using the Piston theory.The temperature is considered uniform over the thickness of the panel.The aero-thermo-elastic model is derived in the framework of the Carrera unified formulation(CUF),therefore the matrices are expressed in a compact form using the″fundamental nuclei″.Composite and sandwich structures are considered and different boundary conditions are taken into account.The effects of the thermal load on the aeroelastic behavior are investigated.     

MAC Models in Thermoelasticity

There are two types of singularities in the linear thermoelasticity.The first one arises in the field of stresses if a force is applied to one point of the body.This singularity is physical and should be accepted.The second type of singularities is nonphysical and they arise in the fields of displacements and temperatures.There exist the nonlocal theories and gradient theories which have the goal to introduce the finite stresses instead of the infinite ones.The MAC model of the thermoelasticity is created to avoid the nonphysical singularities and it accepts the infinite stresses.MAC is the method of additional conditions,which allows introducing the new model to use the classical model,plus additional condition of the physical nonsingularity and/or condition of the good behavior of the solutions at infinity.The MAC Green′s functions for the heat conduction and for the elasticity could be introduced using the differential MAC models.The infinite and finite bodies are considered.The principle of superposition is applied to obtain the integral equations to solve the boundary value problems.The strength criteria based on finite stresses could be changed in this model because the infinite stresses are allowed.The strength criteria based on deformations are applicable.Classification of MAC models is given.     

Quadcopter UAV Modeling and Automatic Flight Control Design

The mathematical model of quadcopter-unmanned aerial vehicle(UAV)is derived by using two approaches:One is the Newton-Euler approach which is formulated using classical mechanics;and other is the Euler-Lagrange approach which describes the model in terms of kinetic(translational and rotational)and potential energy.The proposed quadcopter′s non-linear model is incorporated with aero-dynamical forces generated by air resistance,which helps aircraft to exhibits more realistic behavior while hovering.Based on the obtained model,the suitable control strategy is developed,under which two effective flight control systems are developed.Each control system is created by cascading the proportional-derivative(PD)and T-S fuzzy controllers that are equipped with six and twelve feedback signals individually respectively to ensure better tracking,stabilization,and response.Both proposed flight control designs are then implemented with the quadcopter model respectively and multitudinous simulations are conducted using MATLAB/Simulink to analyze the tracking performance of the quadcopter model at various reference inputs and trajectories.     

Problem of Circular Hole in Thermopiezoelectric Media with Semi-permeable Thermal Boundary Condition

The semi-permeable boundary condition is proposed to discuss the influence of the thermal conductivity acting on the stress and heat flow around the hole. Based on the Stroh formalism, the closed form solutions are de- rived, the stress and heat flow around the hole are discussed. The results show that the thermal boundary condi- tion has significant influence on the hoop stress and heat flow around the hole. The hoop stress decreases dramatic- ally with the increasement of the thermal conduction coefficient.     

SIMULATION OF LARGE CIVIL AIRCRAFT LANDING RESPONSE UNDER CONSIDERATIONS OF FLEXIBLE AIRFRAME

To evaluate the landing response of the large civil aircraft in the conceptual design phase , a method for simulating aircraft landing is given.The model for the shock absorber is investigated.The flexible airframe model is established using finite element model ( FEM ) to analyze its modes.Then , the whole aircraft model with flexible airframe is made for the multibody simulation.Tail-down , two-point , three-point and sideslip landing scenarios are studied.The influence on the landing performance considering mode superposition of the flexible airframe is analyzed.Both longitudinal and spanwise positions of the main landing gear are changed to research the influence on the landing performance.Results show that the method is feasible.The shock absorber axial force of the main landing gear with the flexible airframe is smaller than that of rigid airframe.The number of mode superposition and the position of main landing gear can influence the landing response.     

WAKE GEOMETRY CALCULATIONS FOR TILT-ROTOR USING VISCOUS VORTEX METHOD

A tilt-rotor unsteady flow analytical method has been developed based upon viscous vortex-particle meth- od. In this method, the vorticity field is divided into small assembled vortex particles. Vortex motion and diffusion are obtained by solving the velocity-vorticity-formed incompressible Navier-Stokes equations using a grid-free La- grangian simulation method. Generation of the newly vortex particles is calculated by using the Weissinger-L lifting surface model. Furthermore, in order to significantly improve computational efficiency, a fast multiple method （FMM） is introduced into the calculation of induced velocity and its gradient. Finally, the joint vertical experimen- tal （JVX） tilt-rotor is taken as numerical examples to analyze. The wake geometry and downwash are investigated for both hover and airplane modes. The proposed method for tilt-rotor flow analysis is verified by comparing its re- sults with those available measured data. Comparison indicates that the current method can accurately capture the complicated tilt-rotor wake variation and be suitable for aerodynamic interaction simulation in complex environ- ments. Additionally, the aerodynamic interactional characteristics of dual-rotor wake are discussed in different ro- tor distance. Results show that there are significant differences on interactional characteristics between hover mode and airplane mode.     

Thermodynamic Modeling and Simulation of Air System Control Device

This paper aims to obtain the thermodynamic characteristics of the air system control device sealing part in different compressor bleed air and ambient temperature. On the basis of considering the main factors affecting the heat exchange process and simplifying the physical model of the air system control device,the thermodynamic model of air system control device is established based on the basic theory of laminar flow heat transfer and heat conduction theory.Then the piston motion characteristics and the thermodynamic characteristics of the air system control device seal are simulated. The simulation results show that the valve actuation dynamic time of piston is about 0.13 s in the actual working conditions,and the temperature effect on the dynamic response of the piston rod is only 5 ms when the inlet air temperature at 300 ℃ and 370 ℃. The maximum temperature of the air system control device sealing part is not more than 290 ℃ under long time working condition of compressor air entraining. The highest temperature of the sealing part can reach up to 340 ℃ when the air flow temperature reaches the limit temperature of 370 ℃,and the longest duration working temperature limit is not more than 14 s. Therefore,the selection of control device sealing material should consider the work characteristic of instantaneous temperature limit.     

Construction of Crack Perturbation Model and Forward Semi-analytical Model of Attached Eddy Current Sensor

Fatigue damage monitoring is critical metallic structure health monitoring of aircraft.The sensor should be high sensitive,easy to be integrated into structure and well adaptable for poor working conditions.Therefore,an attached eddy current sensor with flexible plane is put forward and its characteristics are analyzed.By extracting material′s conductivity as the crack features,forward semi-analytical model is established and parameter optimizations are carried out.Crack perturbation model of attached eddy current sensor is constructed,and perturbation voltages of sensing channels under three-dimension structural crack are obtained.To verify the sensor′s performance,monitoring experiment on crack extension is conducted under condition of 3 MHz frequency.The validation experimental results show that perturbation model of 2A12-T4 aluminum alloy agrees well with experiment results,and perturbation model errors of four sensing channels are within 25%.The attached eddy current sensor is capable of testing the crack nondestructively and measuring the crack extension quantitatively with the accuracy of 1mm.     

NUMERICAL STUDY ON EFFECT OF OPERATING TABLE PROTECTED BY HORIZONTAL LAMINAR FLOW SCREEN

Transmission of airborne bacteria is the main factor causing surgical site infection （SSI）, which is harm- ful to patients＇ health and even lives. Numerical study is conducted on the effect of the operating table protected by horizontal laminar flow screen. Discrete phase model （DPM） is used. Numerical simulation is carried out to evalu- ate particle trajectories with the Lagrange approach. As a result, the protecting effect of horizontal laminar flow screen is established, and the protecting parameters of the air velocity supplied by the screen and the protecting dis- tance are optimized. The optimized air velocity supplied by the screen should be at 0.4--0.6 m/s. And the protec- ting distance should be less than 1.3 m. This work provides references for the study on the depuration of operating table or room.     

DYNAMIC RESPONSE ANALYSIS OF CARRIER-BASED AIRCRAFT DURING LANDING

In view of the complexity of landing on the deck of aircraft carrier, a systematic model, composed of six- degree-of-freedom mathematic model of carrier-based aircraft, four-degree-of-freedom model of landing gears and six-degree-of-freedom mathematic model of carrier, is established in the Matlab-Simulink environment, with damping function of landing gears and dynamic characteristics of tires being considered. The model, where the car- rier movement is introduced, is applicable for any abnormal landing condition. Moreover, the equations of motion and relevant parameter are also derived. The dynamic response of aircraft is calculated via the variable step-size Runge-Kuta algorithm. The effect of attitude angles of aircraft and carrier movement during the process of landing is illustrated in details. The analytical results can provide some reference for carrier-based aircraft design and main- tenance.     

Drag Reduction Effect for Hypersonic Lifting-Body Vehicle with Counterflowing Jet

This study takes the novel approach of using a counterflowing jet positioned on the nose of a lifting-body vehicle to explore its drag reduction effect at a range of angles of attack.Numerical studies are conducted at a freestream Mach number of 8 in standard atmospheric conditions corresponding to the altitude of 40 km.The effects of jet pressure ratio and flying angles of attack on drag reduction of the model are systematically investigated.Considering the reverse thrust generated from the counterflowing jet,the drag on the nose at hypersonic speeds could be reduced up to 66%.The maximum lift-to-drag ratio of the model is obtained at 6°;meanwhile,the counterflowing jet produces a drag reduction of 8.8%for the whole model.In addition to the nose,the counterflowing jet influences the drag by increasing the pressure drag of the model and reducing the skin friction drag of the first cone within 8°.The results show that the potential of the counterflowing jet as a means of active flow control for drag reduction is significant in the engineering application on hypersonic lifting-body vehicles.     

Analysis on Potential Conflict Frequency of Intersected Air Routes in Terminal Airspace Design

In order to obtain accurate conflict risks in terminal airspace design, the concept and calculation model of potential conflict frequency for intersected routes are proposed. Conflict frequency is represented by the product of horizontal conflict frequency and vertical conflict probability. The horizontal conflict frequency is derived from the probability density distribution of conflicts in a period of time. Based on the recorded radar trajectory data, the concept and model of ROUTE distance are proposed, and the probability density function of aircraft height at a specified ROUTE distance is deduced by kernel density estimation. Furthermore, vertical conflict probability and its horizontal distribution are achieved. Examples of three intersected arrival and departure route design schemes are studied. Compared with scheme 1, the conflict frequency values of the other two improved schemes decrease to 53 % and 24%, respectively. The results show that the model can quantify potential conflict frequency of intersec ted routes.     

Experiment on Boiling Heat Transfer of Refrigerant R134a in Mini-channels

The flow boiling heat transfer characteristics of refrigerant R134 a flowing inside two different kinds of minichannels are investigated. One channel is multi-port extruded with the hydraulic diameter of 0.63 mm,and the other one is rectangular with offset fins and a hydraulic diameter of 1.28 mm. The experiments are performed with a mass flow rate between 68 and 630 kg/(m~2·s),a heat flux between 9 and 64 kW/m~2,and a saturation pressure between 0.24 and 0.63 MPa,under the constant heat flux heating mode. It is found that the effect of mass flow rate on boiling heat transfer is related to heat flux,and that with the increase of heat flux,the effect can only be efficient in higher vapor quality region. The effects of heat flux and saturation pressure on boiling heat transfer are related to a threshold vapor quality,and the value will gradually decrease with the increase of heat flux or saturation pressure. Based on these analyses,a new correlation is proposed to predict the boiling heat transfer coefficient of refrigerant R134 a in the mini-channels under the experimental conditions.