Large eddy simulation of shock boundary layer interaction control using micro-vortex generators

The performance of supersonic engine inlets and external aerodynamic surfaces can be critically affected by shock wave / boundary layer interactions (SBLIs), whose severe adverse pressure gradients can cause boundary layer separation. Currently such problems are avoided primarily through the use of boundary layer bleed/suction which can be a source of significant performance degradation. This study investigates a novel type of flow control device called micro-vortex generators (µVGs) which may offer similar control benefits without the bleed penalties. µVGs have the ability to alter the near-wall structure of compressible turbulent boundary layers to provide increased mixing of high speed fluid which improves the boundary layer health when subjected to flow disturbance. Due to their small size,µVGs are embedded in the boundary layer which provide reduced drag compared to the traditional vortex generators while they are cost-effective, physically robust and do not require a power source. To examine the potential of µVGs, a detailed experimental and computational study of micro-ramps in a supersonic boundary layer at Mach 3 subjected to an oblique shock was undertaken. The experiments employed a flat plate boundary layer with an impinging oblique shock with downstream total pressure measurements. The moderate Reynolds number of 3,800 based on displacement thickness allowed the computations to use Large Eddy Simulations without the subgrid stress model (LES-nSGS). The LES predictions indicated that the shock changes the structure of the turbulent eddies and the primary vortices generated from the micro-ramp. Furthermore, they generally reproduced the experimentally obtained mean velocity profiles, unlike similarly-resolved RANS computations. The experiments and the LES results indicate that the micro-ramps, whose height is h≈0.5δ, can significantly reduce boundary layer thickness and improve downstream boundary layer health as measured by the incompressible shape factor, H. Regions directly behind the ramp centerline tended to have increased boundary layer thickness indicating the significant three-dimensionality of the flow field. Compared to baseline sizes, smaller micro-ramps yielded improved total pressure recovery. Moving the smaller ramps closer to the shock interaction also reduced the displacement thickness and the separated area. This effect is attributed to decreased wave drag and the closer proximity of the vortex pairs to the wall. In the second part of the study, various types of µVGs are investigated including micro-ramps and micro-vanes. The results showed that vortices generated from µVGs can partially eliminate shock induced flow separation and can continue to entrain high momentum flux for boundary layer recovery downstream. The micro-ramps resulted in thinner downstream displacement thickness in comparison to the micro-vanes. However, the strength of the streamwise vorticity for the micro-ramps decayed faster due to dissipation especially after the shock interaction. In addition, the close spanwise distance between each vortex for the ramp geometry causes the vortex cores to move upwards from the wall due to induced upwash effects. Micro-vanes, on the other hand, yielded an increased spanwise spacing of the streamwise vortices at the point of formation. This resulted in streamwise vortices staying closer to the wall with less circulation decay, and the reduction in overall flow separation is attributed to these effects. Two hybrid concepts, named “thick-vane” and “split-ramp”, were also studied where the former is a vane with side supports and the latter has a uniform spacing along the centerline of the baseline ramp. These geometries behaved similar to the micro-vanes in terms of the streamwise vorticity and the ability to reduce flow separation, but are more physically robust than the thin vanes. Next, Mach number effect on flow past the micro-ramps (h~0.5δ) are examined in a supersonic boundary layer at M=1.4, 2.2 and 3.0, but with no shock waves present. The LES results indicate that micro-ramps have a greater impact at lower Mach number near the device but its influence decays faster than that for the higher Mach number cases. This may be due to the additional dissipation caused by the primary vortices with smaller effective diameter at the lower Mach number such that their coherency is easily lost causing the streamwise vorticity and the turbulent kinetic energy to decay quickly. The normal distance between the vortex core and the wall had similar growth indicating weak correlation with the Mach number; however, the spanwise distance between the two counter-rotating cores further increases with lower Mach number. Finally, various µVGs which include micro-ramp, split-ramp and a new hybrid concept “ramped-vane” are investigated under normal shock conditions at Mach number of 1.3. In particular, the ramped-vane was studied extensively by varying its size, interior spacing of the device and streamwise position respect to the shock. The ramped-vane provided increased vorticity compared to the micro-ramp and the split-ramp. This significantly reduced the separation length downstream of the device centerline where a larger ramped-vane with increased trailing edge gap yielded a fully attached flow at the centerline of separation region. The results from coarse-resolution LES studies show that the larger ramped-vane provided the most reductions in the turbulent kinetic energy and pressure fluctuation compared to other devices downstream of the shock. Additional benefits include negligible drag while the reductions in displacement thickness and shape factor were seen compared to other devices. Increased wall shear stress and pressure recovery were found with the larger ramped-vane in the baseline resolution LES studies which also gave decreased amplitudes of the pressure fluctuations downstream of the shock.

Investigating wetting characteristics on microstructured surfaces for superhydrophobicity and metal microcasting

The engineering of liquid behavior on surfaces is important for infrastructure, transportation, manufacturing, and sensing. Surfaces can be rendered superhydrophobic by microstructuring, and superhydrophobic devices could lead to practical corrosion inhibition, self-cleaning, fluid flow control, and surface drag reduction. To more fully understand how liquid interacts with microstructured surfaces, this dissertation introduces a direct method for determining droplet solid-liquid-vapor interfacial geometry on microstructured surfaces. The technique performs metrology on molten metal droplets deposited onto microstructured surfaces and then frozen. Unlike other techniques, this visualization technique can be used on large areas of curved and opaque microstructured surfaces to determine contact line. This dissertation also presents measurements and models for how curvature and flexing of microstructured polymers affects hydrophobicity. Increasing curvature of microstructured surfaces leads to decreased slide angle for liquid droplets suspended on the surface asperities. For a surface with regularly spaced asperities, as curvature becomes more positive, droplets suspended on the tops of asperities are suspended on fewer asperities. Curvature affects superhydrophobicity because microscopic curvature changes solid-liquid interaction, pitch is altered, and curvature changes the shape of the three phase contact line. This dissertation presents a model of droplet interactions with curved microstructured surfaces that can be used to design microstructure geometries that maintain the suspension of a droplet when curved surfaces are covered with microstructured polymers. Controlling droplet dynamics could improve microfluidic devices and the shedding of liquids from expensive equipment, preventing corrosion and detrimental performance. This dissertation demonstrates redirection of dynamic droplet spray with anisotropic microstructures. Superhydrophobic microstructured surfaces can be economically fabricated using metal embossing masters, so this dissertation describes casting-based microfabrication of metal microstructures and nanostructures. Low melting temperature metal was cast into flexible silicone molds which were themselves cast from microfabricated silicon templates. The flexibility of the silicone mold permits casting of curved surfaces, which this dissertation demonstrates by fabricating a cylindrical metal roller with microstructures. The metal microstructures can be in turn used as a reusable molding tool. This dissertation also describes an industrial investment casting process to produce aluminum molds having integrated microstructures. Unlike conventional micromolding tools, the aluminum mold was large and had complex curved surfaces. The aluminum was cast into curved microstructured ceramic molds which were themselves cast from curved microstructured rubber. Many structures were successfully cast into the aluminum with excellent replication fidelity, including circular, square, and triangular holes. This dissertation demonstrates molding of large, curved surfaces having surface microstructures using the aluminum mold. This work contributes a more full understanding of the phenomenon of superhydrophobicity and techniques for the economic fabrication of superhydrophobic microstructures.

Simulations of droplet dispersion in the wakes of cylinders and icing tunnel spray bars

Simulations of droplet dispersion behind cylinder wakes and downstream of icing tunnel spray bars were conducted. In both cases, a range of droplet sizes were investigated numerically with a Lagrangian particle trajectory approach while the turbulent air flow was investigated with a hybrid Reynolds-Averaged Navier-Stokes/Large-Eddy Simulations approach scheme. In the first study, droplets were injected downstream of a cylinder at sub-critical conditions (i.e. with laminar boundary layer separation). A stochastic continuous random walk (CRW) turbulence model was used to capture the effects of sub-grid turbulence. Small inertia droplets (characterized by small Stokes numbers) were affected by both the large-scale and small-scale vortex structures and closely followed the air flow, while exhibiting a dispersion consistent with that of a scalar flow field. Droplets with intermediate Stokes numbers were centrifuged by the vortices to the outer edges of the wake, yielding an increased dispersion. Large Stokes number droplets were found to be less responsive to the vortex structures and exhibited the least dispersion. Particle concentration was also correlated with vorticity distribution which yielded preferential bias effects as a function of different particle sizes. This trend was qualitatively similar to results seen in homogenous isotropic turbulence, though the influence of particle inertia was less pronounced for the cylinder wake case. A similar study was completed for droplet dispersion within the Icing Research Tunnel (IRT) at the NASA Glenn Research Center, where it is important to obtain a nearly uniform liquid water content (LWC) distribution in the test section (to recreate atmospheric icing conditions).. For this goal, droplets are diffused by the mean and turbulent flow generated from the nozzle air jets, from the upstream spray bars, and from the vertical strut wakes. To understand the influence of these three components, a set of simulations was conducted with a sequential inclusion of these components. Firstly, a jet in an otherwise quiescent airflow was simulated to capture the impact of the air jet on flow turbulence and droplet distribution, and the predictions compared well with experimental results. The effects of the spray bar wake and vertical strut wake were then included with two more simulation conditions, for which it was found that the air jets were the primary driving force for droplet dispersion, i.e. that the spray bar and vertical strut wake effects were secondary.

Evaluation of Adoptive T Cell Therapy and Oncolytic Virotherapy for Treatment of Brain Tumors

Adoptive immunotherapy and oncolytic virotherapy are two promising strategies for treating primary and metastatic malignant brain tumors. We demonstrate the ability of adoptively transferred tumor-specific T cells to rapidly mediate the clearance of established brain tumors in several mouse models. Similar to the clinical situation, tumor recurrences are frequent and result from immune editing of tumors. T cells can eliminate antigen-expressing tumor cells but are not effective against antigen loss variant (ALV) cancer cells that multiply and repopulate a tumor. We show that the level of tumor antigen present affects the success of adoptive T cell therapy. When high levels of antigen are present, tumor stromal cells such as microglia and macrophages present tumor peptide on their surface. As a result, T cells directly eliminate cancer cells and cross-presenting stromal cells and indirectly eliminate ALV cells. We were able to show the first direct evidence of tumor antigen cross-presentation by CD11b+ stromal cells in the brain using soluble, high-affinity T cell receptor monomers. Strategies that target brain tumor stroma or increase antigen shedding from tumor cells leading to increased crosspresentation by stromal cells may improve the clinical success of T cell adoptive therapies. We evaluated one potential strategy to complement adoptive T cell therapy by characterizing the oncolytic effects of myxoma virus (MYXV) in a syngeneic mouse brain tumor model of metastatic melanoma. MYXV is a rabbit poxvirus with strict species tropism for European rabbits. MYXV can also infect mouse and human cancer cell lines due to signaling defects in innate antiviral mechanisms and hyperphosphorylation of Akt. MYXV kills B16.SIY melanoma cells in vitro, and intratumoral injection of virus leads to robust, selective and transient infection of the tumor. We observed that virus treatment recruits innate immune cells iii to the tumor, induces TNFα and IFNβ production in the brain, and results in limited oncolytic effects in vivo. To overcome this, we evaluated the safety and efficacy of co-administering 2C T cells, MYXV, and neutralizing antibodies against IFNβ. Mice that received the triple combination therapy survived significantly longer with no apparent side effects, but eventually relapsed. Based on these findings, methods to enhance viral replication in the tumor and limit immune clearance of the virus will be pursued. We conclude that myxoma virus should be further explored as a vector for transient delivery of therapeutic genes to a tumor to enhance T cell responses.

Superconducting Transition-Edge Sensor Physics

Despite record-setting performance demonstrated by superconducting Transition Edge Sensors (TESs) and growing utilization of the technology, a theoretical model of the physics governing TES devices superconducting phase transition has proven elusive. Earlier attempts to describe TESs assumed them to be uniform superconductors. Sadleir et al. 2010 shows that TESs are weak links and that the superconducting order parameter strength has significant spatial variation. Measurements are presented of the temperature T and magnetic field B dependence of the critical current Ic measured over 7 orders of magnitude on square Mo/Au bilayers ranging in length from 8 to 290 microns. We find our measurements have a natural explanation in terms of a spatially varying order parameter that is enhanced in proximity to the higher transition temperature superconducting leads (the longitudinal proximity effect) and suppressed in proximity to the added normal metal structures (the lateral inverse proximity effect). These in-plane proximity effects and scaling relations are observed over unprecedentedly long lengths (in excess of 1000 times the mean free path) and explained in terms of a Ginzburg-Landau model. Our low temperature Ic(B) measurements are found to agree with a general derivation of a superconducting strip with an edge or geometric barrier to vortex entry and we also derive two conditions that lead to Ic rectification. At high temperatures the Ic(B) exhibits distinct Josephson effect behavior over long length scales and following functional dependences not previously reported. We also investigate how film stress changes the transition, explain some transition features in terms of a nonequilibrium superconductivity effect, and show that our measurements of the resistive transition are not consistent with a percolating resistor network model.