Finite volume distance field solution applied to medial axis transform

Accurate and efficient computation of the nearest wall distance d (or level set) is important for many areas of computational science/engineering. Differential equation-based distance/ level set algorithms, such as the hyperbolic-natured Eikonal equation, have demonstrated valuable computational efficiency. Here, in the context, as an 'auxiliary' equation to the main flow equations, the Eikonal equation is solved efficiently with two different finite volume approaches (the cell vertex and cell-centered). Application of the distance solution is studied for various geometries. Moreover, a procedure using the differential field to obtain the medial axis transform (MAT) for different geometries is presented. The latter provides a skeleton representation of geometric models that has many useful analysis properties. As an alternative approach to the pure geometric methods (e.g. the Voronoi approach), the current d-MAT procedure bypasses many difficulties that are usually encountered by pure geometric methods, especially in three dimensional space. It is also shown that the d-MAT approach provides the potential to sculpt/control the MAT form for specialized solution purposes. Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Unsteady ejectors: The effect of driver jet mark-space ratio

mark Unsteady ejectors can be driven by a wide range of driver jets. These vary from pulse detonation engines, which typically have a long gap between each slug of fluid exiting the detonation tube (mark-space ratios in the range 0.1-0.2) to the exit of a pulsejet where the mean mass flow rate leads to a much shorter gap between slugs (mark-space ratios in the range 2-3). The aim of this paper is to investigate the effect of mark-space ratio on the thrust augmentation of an unsteady ejector. Experimental testing was undertaken using a driver jet with a sinusoidal exit velocity profile. The mean value, amplitude and frequency of the velocity profile could be changed allowing the length to diameter ratio of the fluid slugs L/D and the mark-space ratio (the ratio of slug length to the spacing between slugs) L/S to be varied. The setup allowed L/S of the jet to vary from 0.8 to 2.3, while the L/D ratio of the slugs could take any values between 3.5 and 7.5. This paper shows that as the mark-space ratio of the driver jet is increased the thrust augmentation drops. Across the range of mark-space ratios tested, there is shown to be a drop in thrust augmentation of 0.1. The physical cause of this reduction in thrust augmentation is shown to be a decrease in the percentage time over which the ejector entrains ambient fluid. This is the direct result ofthe space between consecutive slugs in the driver jet decreasing. The one dimensional model reported in Heffer et al. [1] is extended to include the effect of varying L/S and is shown to accurately capture the experimentally measured behavior ofthe ejector. Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Once again on the importance of propagation effects for jet noise modelling

Several turbulent jet noise models starting from the classical Lighthill acoustic analogy to state-of-the art models are considered. No attempt is made to present any complete overview of jet noise theories. Instead, the aim is to emphasise the importance of sound generation and meanflow effects for the understanding and prediction of jet noise. For a recent acoustic analogy model, the consequences of jet flow simplification on the predicted sound spectra shape and the effective noise source location in the jet are discussed. © 2010 by the American Institute of Aeronautics and Astronautics, Inc.

Boundary-layer suction system design for laminar-flying-wing aircraft

The present study aims to provide insight into the parameters affecting practical laminar-flow-control suction power requirements for a commercial laminar-flying-wing transport aircraft. It is shown that there is a minimum power requirement independent of the suction system design, associated with the stagnation pressure loss in the boundary layer. This requirement increases with aerofoil section thickness, but depends only weakly on Mach number and (for a thick, lightly loaded laminar flying wing) lift coefficient. Deviation from the optimal suction distribution, due to a practical chamber-based architecture, is found to have very little effect on the overall suction coefficient; hence, to a good approximation, the power penalty is given by the product of the optimal suction flow rate coefficient and the average skin pressure drop. In the spanwise direction, through suitable choice of chamber depth, the pressure drop due to frictional and inertial effects may be rendered negligible. Finally, if there are fewer pumps than chambers, the average pressure drop from the aerofoil surface to the pump collector ducts, rather than to the chambers, determines the power penalty. For the representative laminar-flying-wing aircraft parameters considered here, the minimum power associated with boundary-layer losses alone contributes some 80-90% of the total power requirement. © 2011 by the American Institute of Aeronautics and Astronautics, Inc.

Corner effect and asymmetry in transonic channel flows

An experimental and numerical investigation into transonic shock/boundary-layer interactions in rectangular ducts has been performed. Experiments have shown that flow development in the corners of transonic shock/boundary-layer interactions in confined channels can have a significant impact on the entire flowfield. As shock strength is increased from M∞ = 1:3 to 1.5, the flowfield becomes very slightly asymmetrical. The interaction of corner flows with one another is thought to be a potential cause of this asymmetry. Thus, factors that govern the size of corner interactions (such as interaction strength) and their proximity to one another (such as tunnel aspect ratio) can affect flow symmetry. The results of the computational study show reasonable agreement with experiments, although simulations with particular turbulence models predict highly asymmetrical solutions for flows that were predominantly symmetrical in experiments. These discrepancies are attributed to the tendency of numerical schemes to overprediction corner-interaction size, and this also accounts for why computational fluid dynamics predicts the onset of asymmetry at lower shock strengths than in experiments. The findings of this study highlight the importance of making informed decisions about imposing artificial constraints on symmetry and boundary conditions for internal transonic flows. Future effort into modeling corner flows accurately is required. Copyright © 2011 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Monte Carlo optimization strategies for air-traffic control

The safety of the flights, and in particular conflict resolution for separation assurance, is one of the main tasks of Air Traffic Control. Conflict resolution requires decision making in the face of the considerable levels of uncertainty inherent in the motion of aircraft. We present a Monte Carlo framework for conflict resolution which allows one to take into account such levels of uncertainty through the use of a stochastic simulator. A simulation example inspired by current air traffic control practice illustrates the proposed conflict resolution strategy. Copyright © 2005 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Small-diameter membrane reflector wrinkling

This study concerns the wrinkling performance of thin membranes for use as novel reflectors in space-based telescopes. We introduce small-scale experiments for inducing and interrogating wrinkling patterns in at membranes, and we capture these details computationally by performing a range of finite element analysis. The overall aim is to assess the sophistication of modelling, to verify the feasibility of a small-diameter reector concept proposed in accompanying work. © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

Robust airfoil design with respect to boundary layer transition

In order to disign an airfoil of which maximum lift coefficient (CL max) is not sensitive to location of forced top boundary layer transition. Taking maximizing mean value of CL max and minimizing standard deviation as biobjective, leading edge radius, manximum thickness and its location, maximum camber and its location as deterministic design variables, location of forced top boundary layer transition as stochastic variable, XFOIL as deterministic CFD solver, non-intrusive polynomial chaos as substitute of Monte Carlo method, we completed a robust airfoil design problem. Results demonstrate performance of initial airfoil is enhanced through reducing standard deviation of CL max. Besides, we also know maximum thickness has the most dominating effect on mean value of CL max, location of maximum thickness has the most dominating effect on standard deviation of CL max, maximum camber has a little effect on both mean value and standard deviation, and maximum camber is the only element of which increase can lead increase of mean value and standard deviation at the same time. Copyright © 2009 by the American Institute of Aeronautics and Astronautics, Inc.

Combustion rumble prediction with integrated computational-fluid-dynamics/ low-order-model methods

The pressure oscillation within combustion chambers of aeroengines and industrial gas turbines is a major technical challenge to the development of high-performance and low-emission propulsion systems. In this paper, an approach integrating computational fluid dynamics and one-dimensional linear stability analysis is developed to predict the modes of oscillation in a combustor and their frequencies and growth rates. Linear acoustic theory was used to describe the acoustic waves propagating upstream and downstream of the combustion zone, which enables the computational fluid dynamics calculation to be efficiently concentrated on the combustion zone. A combustion oscillation was found to occur with its predicted frequency in agreement with experimental measurements. Furthermore, results from the computational fluid dynamics calculation provide the flame transfer function to describe unsteady heat release rate. Departures from ideal one-dimensional flows are described by shape factors. Combined with this information, low-order models can work out the possible oscillation modes and their initial growth rates. The approach developed here can be used in more general situations for the analysis of combustion oscillations. Copyright © 2012 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Experimental study into the flow physics of three-dimensional shock control bumps

This paper describes a fundamental experimental study of the flow structure around a single three-dimensional (3D) transonic shock control bump (SCB) mounted on a flat surface in a wind tunnel. Tests have been carried out with a Mach 1.3 normal shock wave located at a number of streamwise positions relative to the SCB. Details of the flow have been studied using the experimental techniques of schlieren photography, surface oil flow visualization, pressure sensitive paint, and laser Doppler anemometry. The results of the work build on the findings of previous researchers and shed new light on the flow physics of 3D SCBs. It is found that spanwise pressure gradients across the SCB ramp and the shape of the SCB sides affect the magnitude and uniformity of flow turning generated by the bump, which can impact on the spanwise propagation of the quasi-two-dimensional (2D) shock structure produced by a 3DSCB. At the bump crest, vortices can form if the pressure on the crest is significantly lower than at either side of the bump. The trajectories of these vortices, which are relatively weak, are strongly influenced by any spanwise pressure gradients across the bump tail. Asignificant difference between 2D and 3D SCBs highlighted by the study is the impact of spanwise pressure gradients on 3D SCB performance. The magnitude of these spanwise pressure gradients is determined largely by SCB geometry and shock position. Copyright © 2011 by the American Institute of Aeronautics and Astronautics, Inc.

Control of a shock-wave/boundary-layer interaction and subsequent subsonic diffuser using a combination of vortex generators and bleed

An experimental investigation has been undertaken in which vortex generators (VGs) have been employed to inhibit boundary-layer separation produced by the combined adversepressure- gradient of a terminal shock-wave and subsonic diffuser. This setup has been developed as part of a program to produce a more inlet relevant flow-field using a small-scale wind tunnel than previous studies. The resulting flow is dominated by large-scale separation, and as such, is thought to be a good test-bed for flow control. In this investigation, VGs have been added to determine their potential for shock-induced separation mitigation. In line with previous studies, it was observed that the application of VGs alone was not able to significantly alleviate separation overall, because enlarged corner separations was observed. Only when control of the corner separations using corner bleed was employed alongside centre-span control using VGs was a significant improvement in both wall pressure recovery (6% increase) and stagnation pressure recovery (2.4% increase) observed. Copyright © 2012 by the American Institute of Aeronautics and Astronautics, Inc.

Large-eddy simulations and measurements of a small-scale high-speed coflowing jet

Measurements and predictions are made of a short-cowl coflowing jet with a bypass ratio of 8:1. The Reynolds number is 300,000, and the inlet Mach numbers are representative of aeroengine conditions. The low Reynolds number of the measurements makes the case well suited to the assessment of large-eddy-simulation-related strategies. The nozzle concentricity is carefully controlled to deal with the emerging metastability issues of jets with coflow. Measurements of mean quantities and turbulence statistics are made using both laser Doppler anemometry and particle image velocimetry. The simulations are completed on 6× 106, 12× 106, and 50 × 106 cell meshes. To overcome near-wall modeling problems, a hybrid large-eddy-simulation-Reynolds-averaged-Navier-Stokesrelated method is used. The near-wall Reynolds-averaged-Navier-Stokes layer is helpful in preventing nonphysical separation from the nozzle wall.Copyright © 2010 by the American Institute of Aeronautics and Astronautics, Inc.