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.

Comparison of les to LDA and PIV measurements of a small scale high speed coflowing jet

Measurements and predictions are made of a short cowl co-flowing jet with a bypass ratio of 8:1. The Reynolds number for computations and measurements are matched at 300,000 and the Mach numbers representative of realistic jet conditions with core and co flow velocities of 240m/s and 216m/s respectively. The low Reynolds number of the measurements makes the case well suited to the assessment of large eddy resolving computational strategies. Also, the nozzle concentricity was carefully controlled to deal with the emerging metastability issues of jets with coflow. Measurements of mean quantities and turbulence statistics are made using both two dimensional coincident LDA and PIV systems. The computational simulations are completed on a modest 12×106 mesh. The simulation is now being run on a 50×106 mesh using hybrid RANSNLES (Numerical Large Eddy Simulation). Close to the nozzle wall a k-l RANS model is used. For an axisymmetric jet, comparison is made between simulations which use NLES, RANSNLES and also a simple imposed velocity profile where the nozzle is not modeled. The use of a near wall RANS model is shown to be beneficial. When compared with the measurements the NLES results are encouraging. Copyright © 2008 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

High-order DES and zonal LES of a labyrinth seal

Hybrid numerical large eddy simulation (NLES) and detached eddy simulation (DES) methods are assessed on a labyrinth seal geometry. A high sixth order discretization scheme is used and is validated using a test case of a two dimensional vortex. The hybrid approach adopts a new blending function and along with DES is initially validated using a simple cavity flow. The NLES method is also validated outside of RANS zones. It is found that there is very little resolved turbulence in the cavity for the DES simulation. For the labyrinth seal calculations the DES approach is problematic giving virtually no resolved turbulence content. It is seen that over the tooth tips the extent of the LES region is small and is likely to be a strong contributor to excessive flow damping in these regions. On the other hand the zonal Hamilton-Jacobi approach did not suffer from this trait. In both cases the meshes used are considered to be hybrid RANS-LES adequate. Fortunately (or perhaps unfortunately) the DES profiles are in agreement with the time mean experimental measurements. It is concluded that for an inexperienced CFD practitioner this could have wider implications particularly if transient results such as unsteady loading are desired. Copyright © 2012 by the American Institute of Aeronautics and Astronautics, Inc.

Shock wave/boundary-layer interaction control using a combination of vortex generators and bleed

To investigate whether vortex generators can be an effective form of passive flow control an experimental investigation has been conducted in a small-scale wind tunnel. With specific emphasis on supersonic inlet applications flow separation was initiated using a combined terminal shock wave and subsonic diffuser: a configuration that has been developed as a part of a program to produce a more inlet-relevant flowfield in a small-scale wind tunnel than previous studies. When flow control was initially introduced little overall flow improvement was obtained as the losses tended to be redistributed instead of removed. It became apparent that there existed a strong coupling between the center-span flow and the corner flows. As a consequence, only when flow control was applied to both the corner flows and center-span flow was a significant flow improvement obtained. When corner suction and center-span vortex generators were employed in tandem separation was much reduced and wall-pressure and stagnation pressure were notably improved. As a result, when applied appropriately, it is thought that vortex generators do have the potential to reduce the dependence on boundary-layer bleed for the purpose of separation suppression. Copyright © 2012 by Neil Titchener and Holger Babinsky. Published by the American Institute of Aeronautics and Astronautics, Inc.

Differential geometry and design coupling in MDO

The notion of coupling within a design, particularly within the context of Multidisciplinary Design Optimization (MDO), is much used but ill-defined. There are many different ways of measuring design coupling, but these measures vary in both their conceptions of what design coupling is and how such coupling may be calculated. Within the differential geometry framework which we have previously developed for MDO systems, we put forth our own design coupling metric for consideration. Our metric is not commensurate with similar types of coupling metrics, but we show that it both provides a helpful geo- metric interpretation of coupling (and uncoupledness in particular) and exhibits greater generality and potential for analysis than those similar metrics. Furthermore, we discuss how the metric might be profitably extended to time-varying problems and show how the metric's measure of coupling can be applied to multi-objective optimization problems (in unconstrained optimization and in MDO). © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Geometric perspectives on MDO and MDO architectures

Two main perspectives have been developed within the Multidisciplinary Design Optimization (MDO) literature for classifying and comparing MDO architectures: a numerical point of view and a formulation/data flow point of view. Although significant work has been done here, these perspectives have not provided much in the way of a priori information or predictive power about architecture performance. In this report, we outline a new perspective, called the geometric perspective, which we believe will be able to provide such predictive power. Using tools from differential geometry, we take several prominent architectures and describe mathematically how each constructs the space through which it moves. We then consider how the architecture moves through the space which it has constructed. Taken together, these investigations show how each architecture relates to the original feasible design manifold, how the architectures relate to each other, and how each architecture deals with the design coupling inherent to the original system. This in turn lays the groundwork for further theoretical comparisons between and analyses of MDO architectures and their behaviour using tools and techniques derived from differential geometry. © 2012 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Enhancing and developing the practical optimisation capabilities and intelligence of automatic design software

In the modern engineering design cycle the use of computational tools becomes a neces- sity. The complexity of the engineering systems under consideration for design increases dramatically as the demands for advanced and innovative design concepts and engineering products is expanding. At the same time the advancements in the available technology in terms of computational resources and power, as well as the intelligence of the design software, accommodate these demands and make them a viable approach towards the chal- lenge of real-world engineering problems. This class of design optimisation problems is by nature multi-disciplinary. In the present work we establish enhanced optimisation capabil- ities within the Nimrod/O tool for massively distributed execution of computational tasks through cluster and computational grid resources, and develop the potential to combine and benefit from all the possible available technological advancements, both software and hardware. We develop the interface between a Free Form Deformation geometry manage- ment in-house code with the 2D airfoil aerodynamic efficiency evaluation tool XFoil, and the well established multi-objective heuristic optimisation algorithm NSGA-II. A simple airfoil design problem has been defined to demonstrate the functionality of the design sys- tem, but also to accommodate a framework for future developments and testing with other state-of-the-art optimisation algorithms such as the Multi-Objective Genetic Algorithm (MOGA) and the Multi-Objective Tabu Search (MOTS) techniques. Ultimately, heav- ily computationally expensive industrial design cases can be realised within the presented framework that could not be investigated before. © 2012 by the authors. Published by the American Institute of Aeronautics and Astronautics, Inc.

Average corotation of line segments near a point and vortex identification

An easy-to-interpret kinematic quantity measuring the average corotation of material line segments near a point is introduced and applied to vortex identification. At a given point, the vector of average corotation of line segments is defined as the average of the instantaneous local rigid-body rotation over "all planar cross sections" passing through the examined point. The vortex-identification method based on average corotation is a one-parameter, region-type local method sensitive to the axial stretching rate as well as to the inner configuration of the velocity gradient tensor. The method is derived from a well-defined interpretation of the local flow kinematics to determine the "plane of swirling" and is also applicable to compressible and variable-density flows. Practical application to direct numerical simulation datasets includes a hairpin vortex of boundary-layer transition, the reconnection process of two Burgers vortices, a flow around an inclined flat plate, and a flow around a revolving insect wing. The results agree well with some popular local methods and perform better in regions of strong shearing. Copyright © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Multispacecraft refueling optimization considering the J2 perturbation and window constraints

The optimization of a near-circular low-Earth-orbit multispacecraft refueling problem is studied. The refueling sequence, service time, and orbital transfer time are used as design variables, whereas the mean mission completion time and mean propellant consumed by orbital maneuvers are used as design objectives. The J2 term of the Earth's nonspherical gravity perturbation and the constraints of rendezvous time windows are taken into account. A hybridencoding genetic algorithm, which uses normal fitness assignment to find the minimum mean propellant-cost solution and fitness assignment based on the concept of Pareto-optimality to find multi-objective optimal solutions, is presented. The proposed approach is demonstrated for a typical multispacecraft refueling problem. The results show that the proposed approach is effective, and that the J2 perturbation and the time-window constraints have considerable influences on the optimization results. For the problems in which the J2 perturbation is not accounted for, the optimal refueling order can be simply determined as a sequential order or as the order only based on orbitalplane differences. In contrast, for the problems that do consider the J2 perturbation, the optimal solutions obtained have a variety of refueling orders and use the drift of nodes effectively to reduce the propellant cost for eliminating orbital-plane differences. © 2013 by the American Institute of Aeronautics and Astronautics, Inc. All rights reserved.

Joint experimental and numerical approach to three-dimensional shock control bump research

Previous studies of transonic shock control bumps have often been either numerical or experimental. Comparisons between the two have been hampered by the limitations of either approach. The present work aims to bridge the gap between computational fluid dynamics and experiment by planning a joint approach from the outset. This enables high-quality validation data to be produced and ensures that the conclusions of either aspect of the study are directly relevant to the application. Experiments conducted with bumps mounted on the floor of a blowdown tunnel were modified to include an additional postshock adverse pressure gradient through the use of a diffuser as well as introducing boundary-layer suction ahead of the test section to enable the in-flow boundary layer to be manipulated. This has the advantage of being an inexpensive and highly repeatable method. Computations were performed on a standard airfoil model, with the flight conditions as free parameters. The experimental and computational setups were then tuned to produce baseline conditions that agree well, enabling confidence that the experimental conclusions are relevant. The methods are then applied to two different shock control bumps: a smoothly contoured bump, representative of previous studies, and a novel extended geometry featuring a continuously widening tail, which spans the wind-tunnel width at the rear of the bump. Comparison between the computational and experimental results for the contour bump showed good agreement both with respect to the flow structures and quantitative analysis of the boundary-layer parameters. It was seen that combining the experimental and numerical data could provide valuable insight into the flow physics, which would not generally be possible for a one-sided approach. The experiments and computational fluid dynamics were also seen to agree well for the extended bump geometry, providing evidence that, even though thebumpinteracts directly with the wind-tunnel walls, it was still possible to observe the key flow physics. The joint approach is thus suitable even for wider bump geometries. Copyright © 2013 by S. P. Colliss, H. Babinsky, K. Nubler, and T. Lutz. Published by the American Institute of Aeronautics and Astronautics, Inc.

The absorption of sound by helmholtz resonators with and without flow

Helmholtz resonators are commonly used as absorbers of incident acoustic power. Theoretical and experimental investigations have been performed in the four cases of no mean flow, grazing mean flow, bias mean flow and a combination of grazing and bias mean flows. In the absence of a mean flow, the absorption coefficient (deflned as the proportion of incident energy absorbed) is a non-linear function of the acoustic pressure and high incident acoustic pressures are required before the absorption becomes signiflcant. In contrast, when there is a mean flow present, either grazing or bias, the absorption is linear and thus absorption coefficient is independent of the magnitude of the acoustic pressure, and absorption is obtained over a wider range of frequencies. Non-linear effects are only discernible very close to resonance and at very-high amplitude. With grazing mean flow, there is the undesirable effect that sound can be generated over a range of frequencies due to the interaction between the unsteadily shed vorticity waves and the downstream edge of the aperture. This production is not observed when there is a bias flow because here the vorticity is shed all around the rim of the aperture and swept away by the mean flow. When there is both a grazing mean flow and a mean bias flow, we flnd that only a small amount of bias mean flow, compared with grazing mean flow, is required to destroy the production of acoustic energy. © 2002 by the author(s). Published by the American Institute of Aeronautics and Astronautics, Inc.

APSK coded modulation schemes for nonlinear satellite channels with high power and spectral efficiency

A new class of 16-ary Amplitude Phase Shift Keying (APSK) coded modulations deemed double-ring PSK modulations best suited for (satellite) nonlinear channels is proposed. Constellation parameters optimization has been based on geometric and information-theoretic considerations. Furthermore, pre- and post-compensation techniques to reduce the nonlinearity impact have been examined. Digital timing clock and carrier phase have been derived and analyzed for a Turbo coded version of the same new modulation scheme. Finally, the performance of state-of the art Turbo coded modulation for this new 16-ary digital modulation has been investigated and compared to the known TCM schemes. It is shown that for the same coding scheme, double-ring APSK modulation outperforms classical 16-QAM and 16-PSK over a typical satellite nonlinear channel due to its intrinsic robustness against the High Power Amplifier (HPA) nonlinear characteristics. The new modulation is shown to be power- and spectrally-efficient, with interesting applications to satellite communications. © 2002 by the American Institute of Aeronautics and Astronautics, Inc.

PHYSICAL AGING OF HIGH-PERFORMANCE THERMOPLASTICS - ENTHALPY RELAXATION IN PEK-C AND PES-C

The developments of physical aging in phenolphthalein poly(aryl-ether-ketone) (PEK-C) and poly(aryl-ether-sulfone) (PES-C) with time at two aging temperatures up to 20 K below their respective glass transition temperatures (T-g = 495 and 520 K) have been studied using differential scanning calorimetry (DSC). Substantial relaxation within the aging course of several hours were observed by detecting T-g decreasing during physical aging process at the two aging temperatures. The relaxation processes of both polymers are extremely nonlinear and self-retarding. The time dependencies of their enthalpies during the initial stages of annealing were approximately modeled using the Narayanaswamy-Tool model. The structure relaxation parameters obtained from this fitting were used to predict the possibility of physical aging occurring at their respective using temperatures. (C) 1995 John Wiley and Sons, Inc.

MECHANICAL-PROPERTIES OF PHENOLPHTHALEIN POLYETHER KETONE - YIELD STRESS, YOUNGS MODULUS, AND FRACTURE-TOUGHNESS

A series of tensile and three-point bending studies was conducted at various temperatures and loading rates using phenolphthalein polyether ketone (PEK-C). Yield stress, Young's modulus, fracture toughness, and crack opening displacement data were obtained for various conditions. In general, both yield stress and Young's modulus increase with decreasing temperature. However, the relationships between fracture toughness, loading rate, and temperature are very complex. This behavior is due to the simultaneous intersection of viscoelasticity and localized plastic deformation. The increased yield stress is the main factor contributing to the reduction in fracture toughness and crack opening displacement. The relationship between fracture toughness and yield stress are discussed. (C) 1995 John Wiley and Sons, Inc.

THERMAL AND MECHANICAL-PROPERTIES OF PHENOLPHTHALEIN POLYETHERSULFONE POLY(PHENYLENE SULFIDE) BLENDS

The thermal and mechanical properties of phenolphthalein polyethersulfone/poly(phenylene sulfide) (PES-C/PPS) blends were studied using a differential scanning calorimeter, a dynamic mechanical analyzer, and mechanical characterization. The morphologies of fracture surfaces were observed by scanning electron microscopy. The blends are multiphase systems with strong interaction between the two phases. It is of interest that, although the strength and ductility of PPS are lower than those of PES-C, the addition of PPS can improve markedly the impact strength of PES-C without changing its higher strength. The PPS can also act as a flow aid for PES-C. (C) 1995 John Wiley and Sons, Inc.