Shock-Induced Near-Wall Two-Phase Flow Structure Over a Micron-Sized Particles Bed

The present paper studies numerical modelling of near-wall two-phase flows induced by a normal shock wave moving at a constant speed, over a micronsized particles bed. In this two-fluid model, the possibility of particle trajectory intersection is considered and a full Lagrangian formulation of the dispersed phase is introduced. The finiteness of the Reynolds and Mach numbers of the flow around a particle as well as the fineness of the particle sizes are taken into account in describing the interactions between the carrier- and dispersed- phases. For the small mass-loading ratio case, the numerical simulation of flow structure of the two phases is implemented and the profiles of the particle number density are obtained under the constant-flux condition on the wall. The effects of the shock Mach number and the particle size and material density on particle entrainment motion are discussed in detail.The obtained results indicate that interphase non-equilibrium in the velocity and temperature is a common feature for this type of flows and a local particle accumulation zone may form near the envelope of the particle trajectory family.

Assessment of the contribution of surface roughness to airframe noise

The generation of sound by turbulent boundary-layer flow at low Mach number over a rough wall is investigated by applying a theoretical model that describes the scattering of the turbulence near field into sound by roughness elements. Attention is focused on the numerical method to approximately quantify the absolute level of far-field radiated roughness noise. Models for the source statistics are obtained by scaling smooth-wall data by the increased skin friction velocity and boundary-layer thickness for a rough surface. Numerical integration is performed to determine the roughness noise, and it reproduces the spectral characteristics of the available empirical formula and experimental data. Experiments are conducted to measure the radiated sound from two rough plates in an open jet The measured noise spectra of the rough plates are above that of a smooth plate in 1-2.5 kHz frequency and exhibit reasonable agreement with the predicted level. Estimates of the roughness noise for a Boeing 757 sized aircraft wing with idealized levels of surface roughness show that hi the high-frequency region the sound radiated from surface roughness may exceed that from the trailing edge, and higher overall sound pressure levels are observed for the roughness noise. The trailing edge noise is also enhanced by surface roughness somewhat A parametric study indicates that roughness height and roughness density significantly affect the roughness noise with roughness height having the dominant effect The roughness noise directivity varies with different levels of surface roughness. Copyright © 2007 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Effects of surface roughness on airframe noise

The generation of sound by turbulent boundary layer flow at low Mach number over a rough wall is investigated by applying the theoretical model which describes the scattering of the turbulence near field into sound by roughness elements. Attention is focused on the numerical method to approximately quantify the absolute level of the roughness noise radiated to far field. Empirical models for the source statistics are obtained by scaling smooth-wall data through increased skin friction velocity and boundary layer thickness for the rough surface. Numerical integration is performed to determine the roughness noise, and it reproduces the spectral characteristics of the available empirical formula and experimental data. Experiments are conducted to measure the radiated sound from two rough plates in an open jet by four 1/2'' free-field condenser microphones. The measured noise spectra of the rough plates are above that of a smooth plate in 1-2.5 kHz frequency and exhibits encouraging agreement with the predicted spectra. Also, a phased microphone array is utilized to localize the sound source, and it confirms that the rough plates generate higher source strengthes in this frequency range. A parametric study illustrates that the roughness height and roughness density significantly affect the far-field radiated roughness noise with the roughness height having the dominant effect. The estimates of the roughness noise for a Boeing 757 sized aircraft wing show that in high frequency region the sound radiated from surface roughness may exceed that from the trailing edge, and higher overall sound pressure levels for the roughness noise are also observed.

Numerical Investigation of Toroidal Shock Wave Focusing in a Cylindrical Chamber

In this paper, focusing of a toroidal shock wave propagating from an annular shock tube into a cylindrical chamber was investigated numerically with the dispersion controlled dissipation (DCD) scheme. The first case for an incident Mach number of 1.5 was conducted and compared with experiments for validation. Then, several cases were calculated for higher incident Mach numbers varying from 2.0 to 5.0, and complicated flow structures were observed. The numerical study was mainly focused on two aspects: focusing process and flow structures. The process, including diffraction, focusing, and reflection, is displayed to reveal the focusing mechanism, and the flow structures at different incident. Mach numbers are used to demonstrate shock reflection styles and focusing characteristics.

Modal scattering at an impedance transition in a lined flow duct

An explicit Wiener-Hopf solution is derived to describe the scattering of duct modes at a hard-soft wall impedance transition in a circular duct with uniform mean flow. Specifically, we have a circular duct r = 1, - ∞ < x < ∞ with mean flow Mach number M > 0 and a hard wall along x < 0 and a wall of impedance Z along x > 0. A minimum edge condition at x = 0 requires a continuous wall streamline r = 1 + h(x, t), no more singular than h = Ο(x1/2) for x ↓ 0. A mode, incident from x < 0, scatters at x = 0 into a series of reflected modes and a series of transmitted modes. Of particular interest is the role of a possible instability along the lined wall in combination with the edge singularity. If one of the "upstream" running modes is to be interpreted as a downstream-running instability, we have an extra degree of freedom in the Wiener-Hopf analysis that can be resolved by application of some form of Kutta condition at x = 0, for example a more stringent edge condition where h = Ο(x3/2) at the downstream side. The question of the instability requires an investigation of the modes in the complex frequency plane and therefore depends on the chosen impedance model, since Z = Z (ω) is essentially frequency dependent. The usual causality condition by Briggs and Bers appears to be not applicable here because it requires a temporal growth rate bounded for all real axial wave numbers. The alternative Crighton-Leppington criterion, however, is applicable and confirms that the suspected mode is usually unstable. In general, the effect of this Kutta condition is significant, but it is particularly large for the plane wave at low frequencies and should therefore be easily measurable. For ω → 0, the modulus fends to |R001| → (1 + M)/(1 -M) without and to 1 with Kutta condition, while the end correction tends to ∞ without and to a finite value with Kutta condition. This is exactly the same behaviour as found for reflection at a pipe exit with flow, irrespective if this is uniform or jet flow.

DNS of Compressible Turbulent Boundary Layer Over a Blunt Wedge

Direct numerical simulation of spatially evolving compressible boundary layer over a blunt wedge is performed in this paper. The free-stream Mach number is 6 and the disturbance source produced by wall blowing and suction is located downstream of the sound-speed point. Statistics are studied and compared with the results in incompressible flat-plate boundary layer. The mean pressure gradient effects on the vortex structure are studied.

The Interaction Between Shock Waves and Solid Spheres Arrays in a Shock Tube

When a shock wave interacts with a group of solid spheres, non-linear aerodynamic behaviors come into effect. The complicated wave reflections such as the Mach reflection occur in. the wave propagation process. The wave interactions with vortices behind each sphere's wake cause fluctuation in the pressure profiles of shock waves. This paper reports an experimental study for the aerodynamic processes involved in the interaction between shock waves and solid spheres. A schlieren photography was applied to visualize the various shock waves passing through solid spheres. Pressure measurements were performed along different downstream positions. The experiments were conducted in both rectangular and circular shock tubes. The data with respect to the effect of the sphere array, size, interval distance, incident Mach number, etc., on the shock wave attenuation were obtained.

Behaviour analysis of numerical solutions and simulation of coherent structures in compressible mixing layers

For understanding the correctness of simulations the behaviour of numerical solutions is analysed, Tn order to improve the accuracy of solutions three methods are presented. The method with GVC (group velocity control) is used to simulate coherent structures in compressible mixing layers. The effect of initial conditions for the mixing layer with convective Mach number 0.8 on coherent structures is discussed. For the given initial conditions two types of coherent structures in the mixing layer are obtained.

Microdamage Evolution, Energy Dissipation and Their Trans-Scale Effects on Macroscopic Failure

The process of damage evolution concerns various scales, from micro- to macroscopic. How to characterize the trans-scale nature of the process is on the challenging frontiers of solid mechanics. In this paper, a closed trans-scale formulation of damage evolution based on statistical microdamage mechanics is presented. As a case study, the damage evolution in spallation is analyzed with the formulation. Scaling of the formulation reveals that the following dimensionless numbers: reduced Mach number M, damage number S, stress wave Fourier number P, intrinsic Deborah number D*, and the imposed Deborah number De*, govern the whole process of deformation and damage evolution. The evaluation of P and the estimation of temperature increase show that the energy equation can be ignored as the first approximation in the case of spallation. Hence, apart from the two conventional macroscopic parameters: the reduced Mach number M and damage number S, the damage evolution in spallation is mainly governed by two microdamage-relevant parameters: the Deborah numbers D* and De*. Higher nucleation and growth rates of microdamage accelerate damage evolution, and result in higher damage in the target plate. In addition, the mere variation in nucleation rate does not change the spatial distribution of damage or form localized rupture, while the increase of microdamage growth rate localizes the damage distribution in the target plate, which can be characterized by the imposed Deborah number De*.

Normal shock boundary layer control with various vortex generator geometries

Various vortex generators which include ramp, split-ramp and a new hybrid concept "ramped-vane" are investigated under normal shock conditions with a diffuser at Mach number of 1.3. The dimensions of the computational domain were designed using Reynolds Average Navier-Stokes studies to be representative of the flow in an external-compression supersonic inlet. Using this flow geometry, various vortex generator concepts were studied with Implicit Large Eddy Simulation. In general, the ramped-vane provided increased vorticity compared to the other devices and reduced the separation length downstream of the device centerline. In addition, the size, edge gap and streamwise position respect to the shock were studied for the ramped-vane and it was found that a height of about half the boundary thickness and a large trailing edge gap yielded a fully attached flow downstream of the device. This ramped-vane also provided the largest reduction in the turbulent kinetic energy and pressure fluctuations. Additional benefits include negligible drag while the reductions in boundary layer displacement thickness and shape factor were seen compared to other devices. © 2010 by Sang Lee.

Normal shock boundary layer control with various vortex generator geometries

Various vortex generators which include ramp, split-ramp and a new hybrid concept "ramped-vane" are investigated under normal shock conditions with a diffuser at Mach number of 1.3. The dimensions of the computational domain were designed using Reynolds Average Navier-Stokes studies to be representative of the flow in an external-compression supersonic inlet. Using this flow geometry, various vortex generator concepts were studied with Implicit Large Eddy Simulation. In general, the ramped-vane provided increased vorticity compared to the other devices and reduced the separation length downstream of the device centerline. In addition, the size, edge gap and streamwise position respect to the shock were studied for the ramped-vane and it was found that a height of about half the boundary thickness and a large trailing edge gap yielded a fully attached flow downstream of the device. This ramped-vane also provided the largest reduction in the turbulent kinetic energy and pressure fluctuations. Additional benefits include negligible drag while the reductions in boundary layer displacement thickness and shape factor were seen compared to other devices. © 2011 Elsevier Ltd.

Direct Numerical Simulation Of Hypersonic Boundary-Layer Transition Over A Blunt Cone

Direct numerical simulation of transition How over a blunt cone with a freestream Mach number of 6, Reynolds number of 10,000 based on the nose radius, and a 1-deg angle of attack is performed by using a seventh-order weighted essentially nonoscillatory scheme for the convection terms of the Navier-Stokes equations, together with an eighth-order central finite difference scheme for the viscous terms. The wall blow-and-suction perturbations, including random perturbation and multifrequency perturbation, are used to trigger the transition. The maximum amplitude of the wall-normal velocity disturbance is set to 1% of the freestream velocity. The obtained transition locations on the cone surface agree well with each other far both cases. Transition onset is located at about 500 times the nose radius in the leeward section and 750 times the nose radius in the windward section. The frequency spectrum of velocity and pressure fluctuations at different streamwise locations are analyzed and compared with the linear stability theory. The second-mode disturbance wave is deemed to be the dominating disturbance because the growth rate of the second mode is much higher than the first mode. The reason why transition in the leeward section occurs earlier than that in the windward section is analyzed. It is not because of higher local growth rate of disturbance waves in the leeward section, but because the growth start location of the dominating second-mode wave in the leeward section is much earlier than that in the windward section.

Dns Of Compressible Turbulent Boundary Layer Around A Sharp Cone

Direct numerical simulation of the turbulent boundary layer over a sharp cone with 20 degrees cone angle (or 10 degrees half-cone angle) is performed by using the mixed seventh-order up-wind biased finite difference scheme and sixth-order central difference scheme. The free stream Mach number is 0.7 and free stream unit Reynolds number is 250000/inch. The characteristics of transition and turbulence of the sharp cone boundary layer are compared with those of the flat plate boundary layer. Statistics of fully developed turbulent flow agree well with the experimental and theoretical data for the turbulent flat-plate boundary layer flow. The near wall streak-like structure is shown and the average space between streaks (normalized by the local wall unit) keeps approximately invariable at different streamwise locations. The turbulent energy equation in the cylindrical coordinate is given and turbulent energy budget is studied. The computed results show that the effect of circumferential curvature on turbulence characteristics is not obvious.

Analysis and application of ellipticity of stability equations on fluid mechanics

By using characteristic analysis of the linear and nonlinear parabolic stability equations (PSE), PSE of primitive disturbance variables are proved to be parabolic intotal. By using sub-characteristic analysis of PSE, the linear PSE are proved to be elliptical and hyperbolic-parabolic for velocity U, in subsonic and supersonic, respectively; the nonlinear PSE are proved to be elliptical and hyperbolic-parabolic for relocity U + u in subsonic and supersonic, respectively. The methods are gained that the remained ellipticity is removed from the PSE by characteristic and sub-characteristic theories, the results for the linear PSE are consistent with the known results, and the influence of the Mach number is also given out. At the same time, the methods of removing the remained ellipticity are further obtained from the nonlinear PSE.

Numerical simulation of Richtmyer-Meshkov instability

The compressible Navier-Stokes equations discretized with a fourth order accurate compact finite difference scheme with group velocity control are used to simulate the Richtmyer-Meshkov (R-M) instability problem produced by cylindrical shock-cylindrical material interface with shock Mach number Ms = 1.2 and density ratio 1:20 (interior density/outer density). Effect of shock refraction, reflection, interaction of the reflected shock with the material interface, and effect of initial perturbation modes on R-M instability are investigated numerically. It is noted that the shock refraction is a main physical mechanism of the initial phase changing of the material surface. The multiple interactions of the reflected shock from the origin with the interface and the R-M instability near the material interface are the reason for formation of the spike-bubble structures. Different viscosities lead to different spike-bubble structure characteristics. The vortex pairing phenomenon is found in the initial double mode simulation. The mode interaction is the main factor of small structures production near the interface.