Aspects of quality control of wind profiler measurements in complex topography

It is well known in the scientific community that some remote sensing instruments assume that sample volumes present homogeneous conditions within a defined meteorological profile. At complex topographic sites and under extreme meteorological conditions, this assumption may be fallible depending on the site, and it is more likely to fail in the lower layers of the atmosphere. This piece of work tests the homogeneity of the wind field over a boundary layer wind profiler radar located in complex terrain on the coast under different meteorological conditions. The results reveal the qualitative importance of being aware of deviations in this homogeneity assumption and evaluate its effect on the final product. Patterns of behavior in data have been identified in order to simplify the analysis of the complex signal registered. The quality information obtained from the homogeneity study under different meteorological conditions provides useful indicators for the best alternatives the system can offer to build wind profiles. Finally, the results are also to be considered in order to integrate them in a quality algorithm implemented at the product level.

Boundary-layer transition on a slender cone in hypervelocity flow with real gas effects

The laminar to turbulent transition process in boundary layer flows in thermochemical nonequilibrium at high enthalpy is measured and characterized. Experiments are performed in the T5 Hypervelocity Reflected Shock Tunnel at Caltech, using a 1 m length 5-degree half angle axisymmetric cone instrumented with 80 fast-response annular thermocouples, complemented by boundary layer stability computations using the STABL software suite. A new mixing tank is added to the shock tube fill apparatus for premixed freestream gas experiments, and a new cleaning procedure results in more consistent transition measurements. Transition location is nondimensionalized using a scaling with the boundary layer thickness, which is correlated with the acoustic properties of the boundary layer, and compared with parabolized stability equation (PSE) analysis. In these nondimensionalized terms, transition delay with increasing CO₂ concentration is observed: tests in 100% and 50% CO₂, by mass, transition up to 25% and 15% later, respectively, than air experiments. These results are consistent with previous work indicating that CO₂ molecules at elevated temperatures absorb acoustic instabilities in the MHz range, which is the expected frequency of the Mack second-mode instability at these conditions, and also consistent with predictions from PSE analysis. A strong unit Reynolds number effect is observed, which is believed to arise from tunnel noise. N_Tr for air from 5.4 to 13.2 is computed, substantially higher than previously reported for noisy facilities. Time- and spatially-resolved heat transfer traces are used to track the propagation of turbulent spots, and convection rates at 90%, 76%, and 63% of the boundary layer edge velocity, respectively, are observed for the leading edge, centroid, and trailing edge of the spots. A model constructed with these spot propagation parameters is used to infer spot generation rates from measured transition onset to completion distance. Finally, a novel method to control transition location with boundary layer gas injection is investigated. An appropriate porous-metal injector section for the cone is designed and fabricated, and the efficacy of injected CO₂ for delaying transition is gauged at various mass flow rates, and compared with both no injection and chemically inert argon injection cases. While CO₂ injection seems to delay transition, and argon injection seems to promote it, the experimental results are inconclusive and matching computations do not predict a reduction in N factor from any CO₂ injection condition computed.

Efficient methods for stochastic optimal control

The Hamilton Jacobi Bellman (HJB) equation is central to stochastic optimal control (SOC) theory, yielding the optimal solution to general problems specified by known dynamics and a specified cost functional. Given the assumption of quadratic cost on the control input, it is well known that the HJB reduces to a particular partial differential equation (PDE). While powerful, this reduction is not commonly used as the PDE is of second order, is nonlinear, and examples exist where the problem may not have a solution in a classical sense. Furthermore, each state of the system appears as another dimension of the PDE, giving rise to the curse of dimensionality. Since the number of degrees of freedom required to solve the optimal control problem grows exponentially with dimension, the problem becomes intractable for systems with all but modest dimension.

In the last decade researchers have found that under certain, fairly non-restrictive structural assumptions, the HJB may be transformed into a linear PDE, with an interesting analogue in the discretized domain of Markov Decision Processes (MDP). The work presented in this thesis uses the linearity of this particular form of the HJB PDE to push the computational boundaries of stochastic optimal control.

This is done by crafting together previously disjoint lines of research in computation. The first of these is the use of Sum of Squares (SOS) techniques for synthesis of control policies. A candidate polynomial with variable coefficients is proposed as the solution to the stochastic optimal control problem. An SOS relaxation is then taken to the partial differential constraints, leading to a hierarchy of semidefinite relaxations with improving sub-optimality gap. The resulting approximate solutions are shown to be guaranteed over- and under-approximations for the optimal value function. It is shown that these results extend to arbitrary parabolic and elliptic PDEs, yielding a novel method for Uncertainty Quantification (UQ) of systems governed by partial differential constraints. Domain decomposition techniques are also made available, allowing for such problems to be solved via parallelization and low-order polynomials.

The optimization-based SOS technique is then contrasted with the Separated Representation (SR) approach from the applied mathematics community. The technique allows for systems of equations to be solved through a low-rank decomposition that results in algorithms that scale linearly with dimensionality. Its application in stochastic optimal control allows for previously uncomputable problems to be solved quickly, scaling to such complex systems as the Quadcopter and VTOL aircraft. This technique may be combined with the SOS approach, yielding not only a numerical technique, but also an analytical one that allows for entirely new classes of systems to be studied and for stability properties to be guaranteed.

The analysis of the linear HJB is completed by the study of its implications in application. It is shown that the HJB and a popular technique in robotics, the use of navigation functions, sit on opposite ends of a spectrum of optimization problems, upon which tradeoffs may be made in problem complexity. Analytical solutions to the HJB in these settings are available in simplified domains, yielding guidance towards optimality for approximation schemes. Finally, the use of HJB equations in temporal multi-task planning problems is investigated. It is demonstrated that such problems are reducible to a sequence of SOC problems linked via boundary conditions. The linearity of the PDE allows us to pre-compute control policy primitives and then compose them, at essentially zero cost, to satisfy a complex temporal logic specification.

Stochastic optimal control

H. J. Kushner has obtained the differential equation satisfied by the optimal feedback control law for a stochastic control system in which the plant dynamics and observations are perturbed by independent additive Gaussian white noise processes. However, the differentiation includes the first and second functional derivatives and, except for a restricted set of systems, is too complex to solve with present techniques.

This investigation studies the optimal control law for the open loop system and incorporates it in a sub-optimal feedback control law. This suboptimal control law's performance is at least as good as that of the optimal control function and satisfies a differential equation involving only the first functional derivative. The solution of this equation is equivalent to solving two two-point boundary valued integro-partial differential equations. An approximate solution has advantages over the conventional approximate solution of Kushner's equation.

As a result of this study, well known results of deterministic optimal control are deduced from the analysis of optimal open loop control.

Non-Linear Scale Interactions in a Forced Turbulent Boundary Layer

This thesis explores the dynamics of scale interactions in a turbulent boundary layer through a forcing-response type experimental study. An emphasis is placed on the analysis of triadic wavenumber interactions since the governing Navier-Stokes equations for the flow necessitate a direct coupling between triadically consist scales. Two sets of experiments were performed in which deterministic disturbances were introduced into the flow using a spatially-impulsive dynamic wall perturbation. Hotwire anemometry was employed to measure the downstream turbulent velocity and study the flow response to the external forcing. In the first set of experiments, which were based on a recent investigation of dynamic forcing effects in a turbulent boundary layer, a 2D (spanwise constant) spatio-temporal normal mode was excited in the flow; the streamwise length and time scales of the synthetic mode roughly correspond to the very-large-scale-motions (VLSM) found naturally in canonical flows. Correlation studies between the large- and small-scale velocity signals reveal an alteration of the natural phase relations between scales by the synthetic mode. In particular, a strong phase-locking or organizing effect is seen on directly coupled small-scales through triadic interactions. Having characterized the bulk influence of a single energetic mode on the flow dynamics, a second set of experiments aimed at isolating specific triadic interactions was performed. Two distinct 2D large-scale normal modes were excited in the flow, and the response at the corresponding sum and difference wavenumbers was isolated from the turbulent signals. Results from this experiment serve as an unique demonstration of direct non-linear interactions in a fully turbulent wall-bounded flow, and allow for examination of phase relationships involving specific interacting scales. A direct connection is also made to the Navier-Stokes resolvent operator framework developed in recent literature. Results and analysis from the present work offer insights into the dynamical structure of wall turbulence, and have interesting implications for design of practical turbulence manipulation or control strategies.

Aspects of quality control of wind profiler measurements in complex topography

It is well known in the scientific community that some remote sensing instruments assume that sample volumes present homogeneous conditions within a defined meteorological profile. At complex topographic sites and under extreme meteorological conditions, this assumption may be fallible depending on the site, and it is more likely to fail in the lower layers of the atmosphere. This piece of work tests the homogeneity of the wind field over a boundary layer wind profiler radar located in complex terrain on the coast under different meteorological conditions. The results reveal the qualitative importance of being aware of deviations in this homogeneity assumption and evaluate its effect on the final product. Patterns of behavior in data have been identified in order to simplify the analysis of the complex signal registered. The quality information obtained from the homogeneity study under different meteorological conditions provides useful indicators for the best alternatives the system can offer to build wind profiles. Finally, the results are also to be considered in order to integrate them in a quality algorithm implemented at the product level.

Three-dimensional SBLI control for transonic airfoils

The application of shock control to transonic airfoils and wings has been demonstrated widely to have the potential to reduce wave drag. Most of the suggested control devices are two-dimensional, that is they are of uniform geometry in spanwise direction. Examples of such techniques include contour bumps and passive control. Recently it has been observed that a spanwise array of discrete three-dimensional controls can have similar benefits but also offer advantages in terms of installation complexity and drag. This paper describes research carried out in Cambridge into various three-dimensional devices, such as slots, grooves and bumps. In all cases the control device is applied to the interaction of a normal shock wave (M=1.3) with a turbulent boundary layer. Theoretical considerations are proposed to determine how such fundamental experiments can provide estimates of control performance on a transonic wing. The potential of each class of three-dimensional device for wave drag reduction on airfoils is discussed and surface bumps in particular are identified as offering potential drag savings for typical transonic wing applications under cruise conditions.

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.

Behaviour of unsteady transonic shock / boundary layer interactions with three-dimensional effects

An experimental investigation into the response of transonic SBLIs to periodic down-stream pressure perturbations in a parallel walled duct has been conducted. Tests have been carried out with a shock strength of M ∞ = 1.5 for pressure perturbation frequencies in the range 16-90 Hz. Analysis of the steady interaction at M∞ = 1.5 has also been made. The principle measurement techniques were high speed schlieren photography and laser Doppler anemometry. The structure of the steady SBLI was found to be highly three-dimensional, with large corner flows and sidewall SBLIs. These aspects are thought to influence the upstream transmission of pressure information through the interaction by affecting the post-shock flow field, including the extent of regions of secondary supersonic flow. At low frequency, the dynamics of shock motion can be predicted using an inviscid analytical model. At increased frequencies, viscous effects become significant and the shock exhibits unexpected dynamic behaviour, due to a phase lag between the upstream transmission of pressure information in the core flow and in the viscous boundary layers. Flow control in the form of micro-vane vortex generators was found to have a small impact on shock dynamics, due to the effect it had on the post-shock flow field outside the viscous boundary layer region. The relationship between inviscid and viscous effects is developed and potential destabilising mechanisms for SBLIs in practical applications are suggested. Copyright © 2009 by Paul Bruce and Holger Babinsky.

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.

Non-reflective boundary conditions for non-uniform mean velocity inlet and outlet flows

The numerical solution of problems in unbounded physical space requires a truncation of the computational domain to a reasonable size. As a result, the conditions on the artificial boundaries are generally unknown. Assumptions like constant pressure or velocities are only valid in the far field and lead to spurious reflections if applied on the boundaries of the truncated domain. A number of attempts have been made over the past decades to design conditions that prevent such reflections. One approach is based on characteristics. The standard analysis assumes a spatially uniform mean flow field but this is often impractical. In the present paper we show how to extend the formulation to the more general case of a non-uniform mean velocity field. A number of test cases are provided and our results compare favourably with other boundary conditions. In principle the present approach can be extended to include non-uniformities in all variables.

Can fundamental shock-wave/boundary-layer interaction research be relevant to inlet aerodynamics?

The important influence of shock waves on supersonic inlet performance has led to much time and effort being expended in the area of shock wave/boundary layer interaction research (SWBLI) and SWBLI control. In this short review, the impact of SWBLIs on supersonic inlet aerodynamic research is discussed and is contrasted with fundamental SWBLI research. Inlet research focussed on internal flow performance is reviewed, based on the salient results, conclusions, and the limitations of such work. The role of fundamental SWBLI research in relation to supersonic inlet research is considered, and the possible positive impact of improving the link between fundamental SWBLI research and inlet design is considered. A simple flow-field is discussed which is thought to be able to simulate at least some more of the flow physics found in a typical inlet. A brief review of real inlet parameters is then given to help determine appropriate fundamental experimental parameters such as incoming Mach number, incoming boundary-layer thickness and subsonic difiuser angle. Copyright © 2012 by N. Titchener, H. Babinsky, and E. Loth.

Vortical structures on three-dimensional shock control bumps

The vortical wake structure produced by a three-dimensional shock control bump (SCB) is thought to be useful for controlling transonic buffet on airfoils. However, at present the vorticity produced is relatively weak and the production mechanism is not well understood. Using a combined experimental and computational approach, a preliminary investigation on the wake vorticity for different bump geometries has been carried out. The structure of the wake for on and off-design conditions are considered, and the effects on the downstream boundary layer demonstrated. Three main vortical structures are observed: a primary vortex pair, weak inter-bump vortices and shear flow in the lambda-shock region. The effect of pressure gradients on vortex strength is examined and it is found that spanwise pressure gradients on the front section of the bump are the most significant parameter influencing vortex strength. © 2013 by S.P. Colliss et al.

Canonical normal shock wave/boundary-layer interaction flows relevant to external compression inlets

The normal shock wave/boundary-layer interaction is important to the operation and performance of a supersonic inlet, and the normal shock wave/boundary-layer interaction is particularly prominent in external compression inlets. To improve understanding of such interactions, it is helpful to make use of fundamental flows that capture the main elements of inlets, without resorting to the level of complexity and system integration associated with full-geometry inlets. In this paper, several fundamental flowfield configurations have been considered as possible test cases to represent the normal shock wave/boundary-layer interaction aspects found in typical external compression inlets, and it was found that the spillage diffuser more closely retains the basic flow features of an external compression inlet than the other configurations. In particular, this flowfield allows the normal shock Mach number as well as the amount and rate of subsonic diffusion to all be held approximately constant and independent of the application of flow control. In addition, a survey of several external compression inlets was conducted to quantify the flow and geometric parameters of the spillage diffuser relevant to actual inlets. The results indicated that such a flow may be especially relevant if the terminal Mach number is about 1.3 to 1.4, the confinement parameter is around 10%, and the width is around twice or three times the height. In addition, the area expansion downstream of the shock should be limited to the conservative side of incipient stall based on incompressible diffusers. Copyright © 2013 by the authors.

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.