Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

Discontinuous growth mechanisms of mode II crack under high-speed impact conditions

A recoverable plate impact testing technology has been used for studying the growth mechanisms of mode II crack. The results show that interactions of microcracks ahead of a crack tip cause the crack growth unsteadily. Failure mode transitions of materials were observed. Based on the observations, a discontinuous crack growth model was established. Analysis shows that the shear crack grows unsteady as the growth speed is between the Rayleigh wave speed c(R) and the shear wave speed c(s); however, when the growth speed approaches root 2c(s), the crack grows steadily. The transient microcrack growth makes the main crack speed to jump from subsonic to intersonic and the steady growth of all the sub-cracks leads the main crack to grow stably at an intersonic speed.

郭永怀文集

Aerodynamic characteristics research on wide-speed range waverider configuration

Waverider generated from a given flow field has a high lift-to-drag ratio because of attached bow shock on leading edge. However, leading edge blunt and off-design condition can make bow shock off leading edge and have unfavorable influence on aerodynamic characteristics. So these two problems have always been concerned as important engineering science issues by aeronautical engineering scientists. In this paper, through respectively using low speed and high speed waverider design principles, a wide-speed rang vehicle is designed, which can level takeoff and accelerate to hypersonic speed for cruise. In addition, sharp leading edge is blunted to alleviated aeroheating. Theoretical study and wind tunnel test show that this vehicle has good aerodynamic performance in wide-speed range of subsonic, transonic, supersonic and hypersonic speeds.

General-domain compressible Navier-Stokes solvers exhibiting quasi-unconditional stability and high-order accuracy in space and time

This thesis presents a new class of solvers for the subsonic compressible Navier-Stokes equations in general two- and three-dimensional spatial domains. The proposed methodology incorporates: 1) A novel linear-cost implicit solver based on use of higher-order backward differentiation formulae (BDF) and the alternating direction implicit approach (ADI); 2) A fast explicit solver; 3) Dispersionless spectral spatial discretizations; and 4) A domain decomposition strategy that negotiates the interactions between the implicit and explicit domains. In particular, the implicit methodology is quasi-unconditionally stable (it does not suffer from CFL constraints for adequately resolved flows), and it can deliver orders of time accuracy between two and six in the presence of general boundary conditions. In fact this thesis presents, for the first time in the literature, high-order time-convergence curves for Navier-Stokes solvers based on the ADI strategy---previous ADI solvers for the Navier-Stokes equations have not demonstrated orders of temporal accuracy higher than one. An extended discussion is presented in this thesis which places on a solid theoretical basis the observed quasi-unconditional stability of the methods of orders two through six. The performance of the proposed solvers is favorable. For example, a two-dimensional rough-surface configuration including boundary layer effects at Reynolds number equal to one million and Mach number 0.85 (with a well-resolved boundary layer, run up to a sufficiently long time that single vortices travel the entire spatial extent of the domain, and with spatial mesh sizes near the wall of the order of one hundred-thousandth the length of the domain) was successfully tackled in a relatively short (approximately thirty-hour) single-core run; for such discretizations an explicit solver would require truly prohibitive computing times. As demonstrated via a variety of numerical experiments in two- and three-dimensions, further, the proposed multi-domain parallel implicit-explicit implementations exhibit high-order convergence in space and time, useful stability properties, limited dispersion, and high parallel efficiency.

Particle fluid mechanics in shear flows, acoustic waves, and shock waves

Three different categories of flow problems of a fluid containing small particles are being considered here. They are: (i) a fluid containing small, non-reacting particles (Parts I and II); (ii) a fluid containing reacting particles (Parts III and IV); and (iii) a fluid containing particles of two distinct sizes with collisions between two groups of particles (Part V).

Part I

A numerical solution is obtained for a fluid containing small particles flowing over an infinite disc rotating at a constant angular velocity. It is a boundary layer type flow, and the boundary layer thickness for the mixture is estimated. For large Reynolds number, the solution suggests the boundary layer approximation of a fluid-particle mixture by assuming W = W_p. The error introduced is consistent with the Prandtl’s boundary layer approximation. Outside the boundary layer, the flow field has to satisfy the “inviscid equation” in which the viscous stress terms are absent while the drag force between the particle cloud and the fluid is still important. Increase of particle concentration reduces the boundary layer thickness and the amount of mixture being transported outwardly is reduced. A new parameter, β = 1/Ω τ_v, is introduced which is also proportional to μ. The secondary flow of the particle cloud depends very much on β. For small values of β, the particle cloud velocity attains its maximum value on the surface of the disc, and for infinitely large values of β, both the radial and axial particle velocity components vanish on the surface of the disc.

Part II

The “inviscid” equation for a gas-particle mixture is linearized to describe the flow over a wavy wall. Corresponding to the Prandtl-Glauert equation for pure gas, a fourth order partial differential equation in terms of the velocity potential ϕ is obtained for the mixture. The solution is obtained for the flow over a periodic wavy wall. For equilibrium flows where λ_v and λ_T approach zero and frozen flows in which λ_v and λ_T become infinitely large, the flow problem is basically similar to that obtained by Ackeret for a pure gas. For finite values of λ_v and λ_T, all quantities except v are not in phase with the wavy wall. Thus the drag coefficient C_D is present even in the subsonic case, and similarly, all quantities decay exponentially for supersonic flows. The phase shift and the attenuation factor increase for increasing particle concentration.

Part III

Using the boundary layer approximation, the initial development of the combustion zone between the laminar mixing of two parallel streams of oxidizing agent and small, solid, combustible particles suspended in an inert gas is investigated. For the special case when the two streams are moving at the same speed, a Green’s function exists for the differential equations describing first order gas temperature and oxidizer concentration. Solutions in terms of error functions and exponential integrals are obtained. Reactions occur within a relatively thin region of the order of λ_D. Thus, it seems advantageous in the general study of two-dimensional laminar flame problems to introduce a chemical boundary layer of thickness λ_D within which reactions take place. Outside this chemical boundary layer, the flow field corresponds to the ordinary fluid dynamics without chemical reaction.

Part IV

The shock wave structure in a condensing medium of small liquid droplets suspended in a homogeneous gas-vapor mixture consists of the conventional compressive wave followed by a relaxation region in which the particle cloud and gas mixture attain momentum and thermal equilibrium. Immediately following the compressive wave, the partial pressure corresponding to the vapor concentration in the gas mixture is higher than the vapor pressure of the liquid droplets and condensation sets in. Farther downstream of the shock, evaporation appears when the particle temperature is raised by the hot surrounding gas mixture. The thickness of the condensation region depends very much on the latent heat. For relatively high latent heat, the condensation zone is small compared with Ʌ_D.

For solid particles suspended initially in an inert gas, the relaxation zone immediately following the compression wave consists of a region where the particle temperature is first being raised to its melting point. When the particles are totally melted as the particle temperature is further increased, evaporation of the particles also plays a role.

The equilibrium condition downstream of the shock can be calculated and is independent of the model of the particle-gas mixture interaction.

Part V

For a gas containing particles of two distinct sizes and satisfying certain conditions, momentum transfer due to collisions between the two groups of particles can be taken into consideration using the classical elastic spherical ball model. Both in the relatively simple problem of normal shock wave and the perturbation solutions for the nozzle flow, the transfer of momentum due to collisions which decreases the velocity difference between the two groups of particles is clearly demonstrated. The difference in temperature as compared with the collisionless case is quite negligible.

Transonic wind tunnel for investigation of turbomachinery blades

The design of a two-stream wind tunnel was undertaken to allow the simulation and study of certain features of the flow field around the blades of high-speed axial-flow turbomachineries. The mixing of the two parallel streams with designed Mach numbers respectively equal to 1.4 and 0.7 will simulate the transonic Mach number distribution generally obtained along the tips of the first stage blades in large bypass-fan engines.

The GALCIT hypersonic compressor plant will be used as an air supply for the wind tunnel, and consequently the calculations contained in the first chapter are derived from the characteristics and the performance of this plant.

The transonic part of the nozzle is computed by using a method developed by K. O. Friedrichs. This method consists essentially of expanding the coordinates and the characteristics of the flow in power series. The development begins with prescribing, more or less arbitrarily, a Mach number distribution along the centerline of the nozzle. This method has been programmed for an IBM 360 computer to define the wall contour of the nozzle.

A further computation is carried out to correct the contour for boundary layer buildup. This boundary layer analysis included geometry, pressure gradient, and Mach number effects. The subsonic nozzle is calculated {including boundary layer buildup) by using the same computer programs. Finally, the mixing zone downstream of the splitter plate was investigated to prescribe the wall contour correction necessary to ensure a constant-pressure test section.

Advancements in jet turbulence and noise modeling: accurate one-way solutions and empirical evaluation of the nonlinear forcing of wavepackets

Jet noise reduction is an important goal within both commercial and military aviation. Although large-scale numerical simulations are now able to simultaneously compute turbulent jets and their radiated sound, lost-cost, physically-motivated models are needed to guide noise-reduction efforts. A particularly promising modeling approach centers around certain large-scale coherent structures, called wavepackets, that are observed in jets and their radiated sound. The typical approach to modeling wavepackets is to approximate them as linear modal solutions of the Euler or Navier-Stokes equations linearized about the long-time mean of the turbulent flow field. The near-field wavepackets obtained from these models show compelling agreement with those educed from experimental and simulation data for both subsonic and supersonic jets, but the acoustic radiation is severely under-predicted in the subsonic case. This thesis contributes to two aspects of these models. First, two new solution methods are developed that can be used to efficiently compute wavepackets and their acoustic radiation, reducing the computational cost of the model by more than an order of magnitude. The new techniques are spatial integration methods and constitute a well-posed, convergent alternative to the frequently used parabolized stability equations. Using concepts related to well-posed boundary conditions, the methods are formulated for general hyperbolic equations and thus have potential applications in many fields of physics and engineering. Second, the nonlinear and stochastic forcing of wavepackets is investigated with the goal of identifying and characterizing the missing dynamics responsible for the under-prediction of acoustic radiation by linear wavepacket models for subsonic jets. Specifically, we use ensembles of large-eddy-simulation flow and force data along with two data decomposition techniques to educe the actual nonlinear forcing experienced by wavepackets in a Mach 0.9 turbulent jet. Modes with high energy are extracted using proper orthogonal decomposition, while high gain modes are identified using a novel technique called empirical resolvent-mode decomposition. In contrast to the flow and acoustic fields, the forcing field is characterized by a lack of energetic coherent structures. Furthermore, the structures that do exist are largely uncorrelated with the acoustic field. Instead, the forces that most efficiently excite an acoustic response appear to take the form of random turbulent fluctuations, implying that direct feedback from nonlinear interactions amongst wavepackets is not an essential noise source mechanism. This suggests that the essential ingredients of sound generation in high Reynolds number jets are contained within the linearized Navier-Stokes operator rather than in the nonlinear forcing terms, a conclusion that has important implications for jet noise modeling.

Multimode radiation from an unflanged, semi-infinite circular duct with uniform flow.

Multimode sound radiation from an unflanged, semi-infinite, rigid-walled circular duct with uniform subsonic mean flow everywhere is investigated theoretically. The multimode directivity depends on the amplitude and directivity function of each individual cut-on mode. The amplitude of each mode is expressed as a function of cut-on ratio for a uniform distribution of incoherent monopoles, a uniform distribution of incoherent axial dipoles, and for equal power per mode. The directivity function of each mode is obtained by applying a Lorentz transformation to the zero-flow directivity function, which is given by a Wiener-Hopf solution. This exact numerical result is compared to an analytic solution, valid in the high-frequency limit, for multimode directivity with uniform flow. The high-frequency asymptotic solution is derived assuming total transmission of power at the open end of the duct, and gives the multimode directivity function with flow in the forward arc for a general family of mode amplitude distribution functions. At high frequencies the agreement between the exact and asymptotic solutions is shown to be excellent.

Performance of a Boundary Layer Ingesting (BLI) propulsion system

This paper presents an assessment of the performance of an embedded propulsion system in the presence of distortion associated with boundary layer ingestion. For fan pressure ratios of interest for civil transports, the benefits of boundary layer ingestion are shown to be very sensitive to the magnitude of fan and duct losses. The distortion transfer across the fan, basically the comparison of the stagnation pressure non-uniformity downstream of the fan to that upstream of the fan, has a major role in determining the impact of boundary layer ingestion on overall fuel burn. This, in turn, puts requirements on the fidelity with which one needs to assess the distortion transfer, and thus the type of models that need to be used in such assessment. For the three-dimensional distortions associated with fuselage boundary layers ingested into a subsonic diffusing inlet, it is found that boundary layer ingestion can provide decreases in fuel burn of several per cent. It is also shown that a promising avenue for mitigating the risks (aerodynamic as well as aeromechanical) in boundary layer ingestion is to mix out the flow before it reaches the engine face.

Micro-vortex generator flow control for supersonic engine inlets

To investigate the flow control potential of micro-vortex generators for supersonic mixed-compression inlets, a basic model experiment has been designed which combines a normal shock wave with a subsonic diffuser. The diffuser is formed by a simple expansion corner, with a divergence angle of 6 degrees. The diffuser entry Mach numbers were M=1.3 and M=1.5 and a number of shock locations relative to the corner position were tested. Flow control was applied in the form of counter-rotating micro-vanes with heights of approximately 20% of boundary layer thickness. Furthermore, corner fences where employed to reduce sidewall effects. It was found that micro-vortex generators were able to significantly reduce the extent of flow separation under all conditions, but could not eliminate it altogether. Corner fences also demonstrated potential for improving the flow in rectangular cross section channels and the combination of corner fences with micro-vortex generators was found to give the greatest benefits. At M=1.3 the combination of corner fences and micro-vanes placed close to the diffuser entry could prevent separation for a wide range of conditions. At the higher diffuser entry Mach number the benefits of flow control were less significant although a reduction of separation size and an improved pressure recovery was observed. It is thought that micro-vortex generators can have significant flow control potential if they are placed close to the expected separation onset and when the adverse pressure gradient is not too far above the incipient separation level. The significant beneficial effects of corner fences warrant a more comprehensive further investigation. It is thought that the control methods suggested here are capable of reducing the bleed requirement on an inlet, which could provide significant performance advantages.

Computational fluid-structure interaction methods for simulation of inflatable aerodynamic decelerators

Inflatable aerodynamic decelerators have potential advantages for planetary re-entry in robotic and human exploration missions. It is theorized that volume-mass characteristics of these decelerators are superior to those of common supersonic/subsonic parachutes and after deployment they may suffer no instabilities at high Mach numbers. A high fidelity computational fluid-structure interaction model is employed to investigate the behavior of tension cone inflatable aeroshells at supersonic speeds up to Mach 2.0. The computational framework targets the large displacements regime encountered during the inflation of the decelerator using fast level set techniques to incorporate boundary conditions of the moving structure. The preliminary results indicate large but steady aeroshell displacement with rich dynamics, including buckling of the inflatable torus that maintains the decelerator open under normal operational conditions, owing to interactions with the turbulent wake. Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

High frequency multimode radiation from ducts with flow

Multimode sound radiation from hard-walled semi-infinite ducts with uniform subsonic flow is investigated theoretically. An analytic expression, valid in the high frequency limit, is derived for the multimode directivity function in the forward arc for a general family of mode distribution functions. The multimode directivity depends on the amplitude and directivity function of each individual mode. The amplitude of each mode is expressed as a function of cut-off ratio for a uniform distribution of incoherent monopoles, a uniform distribution of incoherent axial dipoles and for equal power per mode. The modes' directivity functions are obtained analytically by applying a Lorentz transformation to the zero flow solution. The analytic formula for the multimode directivity with flow is derived assuming total transmission of power at the open-end of the duct. This formula is compared to the exact numerical result for an unflanged duct, computed utilizing a Wiener-Hopf solution. The agreement is shown to be excellent. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

VGs for a normal SBLI with a downstream diffuser

New types of vortex generators for boundary layer control were investigated experimentally in a flow field which contains a Mach 1.4 normal Shockwave followed by a subsonic diffuser. A parametric study of device height and distance upstream of the normal shock was undertaken with two novel devices: ramped-vanes and split-ramps. Flowfield diagnostics included high-speed Schlieren, oil flow visualization, and pitot-static pressure measurements. A number of flowfield parameters including flow separation, pressure recovery, centerline incompressible boundary layer shape factor, and shock stability were analyzed and compared to the baseline. All configurations tested yielded an elimination of centerline flow separation with the presence of the vortex generators. However, the devices also tended to increase the three-dimensionality of the flow with increased side-wall interaction. When located 25δo upstream of the normal shock, the largest ramped-vane device (whose height was about 0.75 the incoming uncontrolled boundary layer thickness, δo) yielded the smallest centerline incompressible shape factor and the least streamwise oscillations of the normal shock. However, additional studies are needed to better understand the corner interaction effects, which are substantial. © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Experimental investigation into parameters governing corner interactions for transonic shock wave/ boundary layer interactions

Experiments were conducted investigating the interaction between a normal shock wave and a corner boundary layer in a constant area rectangular duct. Active corner suction and passive blowing were applied to manipulate the natural corner flows developing in the working section of the Cambridge University supersonic wind tunnel. In addition robust vane micro-vortex generators were applied to the corners of the working section. Experiments were conducted at Mach numbers of M∞=1.4 and 1.5. Flow visualisation was carried out through schlieren and surface oil flow, while static pressures were recorded via floor tappings. The results indicate that an interplay occurs between the corner flow and the centre line flow. It is believed that corner flow separation acts to induce a shock bifurcation, which in turn leads to a smearing of the adverse pressure gradient elsewhere. In addition the blockage effect from the corners was seen to result in a reacceleration of the subsonic post-shock flow. As a result manipulation of the corner regions allows a separated or attached centre line flow to be observed at the same Mach number. Copyright © 2010 by Babinsky, Burton, Bruce.