Search NASA⌕ Search

SEARCH · Search NASA

Results for “vortex structures”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,531 records · Page 85

Breaking Kelvin-Helmholtz waves and cloud-top entrainment as revealed by K-band Doppler radar

Radars have occasionally detected breaking Kelvin-Helmholtz (KH) waves under clear-air conditions in the atmospheric boundary layer and in the free troposphere. However, very few direct measurements of such waves within clouds have previously been reported and those have not clearly documented wave breaking. In this article, we present some of the most detailed and striking radar observations to date of breaking KH waves within clouds and at cloud top and discuss their relevance to the issue of cloud-top entrainment, which is believed to be important in convective and stratiform clouds. Aircraft observations reported by Stith suggest that vortex-like circulations near cloud top are an entrainment mechanism in cumuliform clouds. Laboratory and modeling studies have examined possibility that KH instability may be responsible for mixing at cloud top, but direct observations have not yet been presented. Preliminary analyses shown here may help fill this gap. The data presented in this paper were obtained during two field projects in 1991 that included observations from the NOAA Wave Propagation Laboratory's K-band Doppler radar (wavelength = 8.7 mm) and special rawinsonde ascents. The sensitivity (-30 dBZ at 10 km range), fine spatial resolution (375-m pulse length and 0.5 degrees beamwidth), velocity measurement precision (5-10 cm s-1), scanning capability, and relative immunity to ground clutter make it sensitive to non-precipitating and weakly precipitating clouds, and make it an excellent instrument to study gravity waves in clouds. In particular, the narrow beam width and short pulse length create scattering volumes that are cylinders 37.5 m long and 45 m (90 m) in diameter at 5 km (10 km) range. These characteristics allow the radar to resolve the detailed structure in breaking KH waves such as have been seen in photographic cloud images.

Martner, Brooks E.↗

Vortex-Body Interactions: A Critical Assessment. Coupled Gap-Wake Instabilities/Turbulence: A Source of Noise

This program has involved, first of all, a critical state-of-the-art assessment of vortex-body interactions. Then, efforts were focused on experimental investigation on coupled-wake instabilities and turbulence occurring in a two-cylinder system. An extensive review was undertaken on the effect of incident vortices on various types of bodies. These incident vortices have a length scale of the same order of magnitude as the scale of the body. The body can take on various forms, including, for example, a circular cylinder, a blade or a wing. The classes of vortex-body interaction that were critically assessed include: (1) Periodic distortion of the incident (primary) vortex and shedding of secondary vorticity from the surface of the body. (2) Modulated vortex distortion and shedding at a leading-edge or surface due to incidence of a complex system of vortices. (3) Vortex distortion and shedding in presence of body oscillation. (4) Three-dimensional vortex interaction and shedding. For all of these classes of vortex-body interaction, quantitative topologies of the vorticity distributions and streamline patterns were found to be central to a unified description of mechanisms of vortex distortion and shedding. In most cases, it was possible to define relationships between vortex interactions and unsteady loading at the body surface. This phase of the program was an experimental investigation of a two-cylinder system, which simulated a central aspect of a four-wheel bogie on a large-scale commercial aircraft. The overall aim of this experimental research program was to determine the crucial elements of the unsteadiness in the gap and near-wake regions as a function of time using cinema-based techniques. During the research program, various image evaluation techniques were employed. They involved assessment of instantaneous velocity fields, streamline topology and patterns of vorticity. Experiments were performed in a large-scale water channel using a high-resolution version of digital particle image velocimetry. The program has focused on acquisition of images of velocity and vorticity for varying gap widths between the two-cylinder system. As a result of analysis of a relatively large number of images, it is demonstrated that low frequency instabilities can occur in the gap region between the cylinder. These low frequency instabilities are hypothesized to influence the near-wake structure of the entire two-cylinder system. The nature of the unstable shear layers in the gap region involves generation of small-scale Kelvin-Helmholtz instabilities. These unsteady shear layers then impinge upon the upper and lower surfaces of the cylinders, thereby influencing both the unsteady structure and the time-averaged patterns of the near-wake. Initial efforts have focused on characterization of the patterns of instantaneous and averaged streamlines using topological concepts. The end result of this investigation is a series of documented instantaneous images. They will serve as a basis for various types of post-processing, which will lead to a fuller understanding of the instantaneous and time-averaged unstable-turbulent fields in the gap region and downstream of the two-cylinder system. This further assessment is the focus of a subsequent program.

Rockwell, Donald↗

On the nonlinear development of the most unstable Goertler vortex mode

The nonlinear development of the most unstable Gortler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Gortler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Gortler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Gortler equations valid for O(1) Gortler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Denier, James P.↗

On the nonlinear development of the most unstable Goertler vortex mode TLSP - Final report

The nonlinear development of the most unstable Goertler vortex mode in boundary layer flows over curved walls is investigated. The most unstable Goertler mode is confined to a viscous wall layer of thickness O(G -1/5) and has spanwise wavelength O(G 11/5); it is, of course, most relevant to flow situations where the Goertler number G is much greater than 1. The nonlinear equations covering the evolution of this mode over an O(G -3/5) streamwise lengthscale are derived and are found to be of a fully nonparallel nature. The solution of these equations is achieved by making use of the numerical scheme used by Hall (1988) for the numerical solution of the nonlinear Goertler equations valid for O(1) Goertler numbers. Thus, the spanwise dependence of the flow is described by a Fourier expansion, whereas the streamwise and normal variations of the flow are dealt with by employing a suitable finite difference discretization of the governing equations. Our calculations demonstrate that, given a suitable initial disturbance, after a brief interval of decay, the energy in all the higher harmonics grows until a singularity is encountered at some downstream position. The structure of the flowfield as this singularity is approached suggests that the singularity is responsible for the vortices, which are initially confined to the thin viscous wall layer, moving away from the wall and into the core of the boundary layer.

Denier, James P.↗

Experiments with Turbulent Jets at Mach Number 0.9

A systematic investigation of the structure of turbulent jets before their interaction with shock or expansion waves was undertaken during the last year. In particular compressibility and density effects in circular jets issuing in still air were investigated experimentally. Jets with nitrogen, helium, and krypton gases at 0.3, 0.6, and 0.9 Mach numbers were investigated in detail. Particle Image Velocimetry technique was developed, tested, and used to obtain qualitative information of the two-dimensional velocity field on a plane inside the flow field, which was illuminated by a laser sheet. The motion of particles was recorded by a CCD camera, which was appropriately triggered to capture two images within a fraction of a microsecond. Statistical averaging of the data at each location reduced the large amount of acquired data. It was found that the spreading rate of the jets was reduced with increased Mach numbers or increased density ratio. It was also found that decay rates of centerline Mach numbers are higher in gases with reduced density ratio. Mach number fluctuations appear to decrease with increasing Mach number of the flow. It has been proposed that the reason for this behavior is the reduction of vortex stretching activities with increased Mach number.

Agui, Juan↗

Dynamic, Adaptive, Systems and Materials: Complex, Simple and Emergent Behaviors

This program has been funded by DoE/BES for twenty years. It has moved into and out of various subjects as it has developed, but it has retained its focus on complexity and complex systems. The project has evolved in the following way: Self-Assembly and Biomimetic Self-Assembly: All self-assembling systems depend upon a minimum of two types of interaction: a repulsion and an attraction. For the familiar molecular systems, attractive interactions are typically hydrogen bonds and electrostatic interactions. Repulsive interactions include steric effects, hydrophobic effects (in biological systems), and charge-charge repulsion. We have expanded this repertoire to include surface interactions, magnetic interactions, and others. I list these systems in the order in which we have explored them: i) A key emphasis in current work is in understanding how the movement of ions in a magnetic field (the Lorentz effect) interacts with catalytic systems. We have demonstrated that an acceleration in rate of reduction of CO 2 to CO can be accomplished by applying an external magnetic field. This acceleration is largely due to the application of the Lorentz effect on mass transport at the catalyst’s surface. ii) We have also extensively explored the influence of electrostatics, as exhibited in self-assembling systems, by tribocharging. iii) Another key system involves surface tension effects; examples include interactions between heavy particles floating at a liquid-air interface, and interacting by changes in surface area; interactions of bubbles and bubble rafts, behaviors of bubble trains in microfluidic networks, and behaviors of microorganisms in constraining environments. iv) This work has intentionally de-emphasized biological systems; but it does include some work on protein-ligand interactions and interactions among microorganisms. v) We have also explored applications of some of these effects, these explorations include bubble rafts as diffraction gratings, exploration of the structures that can be obtained by tribocharging and uses of these structures in exploring nucleation and melting of crystals. vi) Although not a major focus of this work, several other topics have emerged and offer opportunities for future work. These include the behavior of bubble trains and bubble rafts in microfluidic systems. A particularly interesting example is the formation of bubble trains that repeat in the alteration of large and small bubbles according to rules we do not presently understand, but are uniquely large-period oscillating systems. These systems offer a new route into understanding the instabilities of the type represented by oscillations. vii) We have also begun exploratory projects on magnetic levitation (especially to determine molecular density), and information storage (in molecules). Magnetic Levitation: Self-assembling and biomimetic systems require both attraction and repulsion. We have used electrostatics (tribocharging), interfacial free-energies (surface tension and related forces) and others. Potential uses include reconfigurable diffraction gratings and liquid lenses; exploration of mechanisms in tribocharging; tunneling in EGaIn junctions; and bubble trains (especially in micro-fluidic systems). Examples of systems representing these topics is included in the following papers: Complexity: Disks rotating at a water-air interface; Benard-Marangoni effects; Vortex-Crystals from spinning magnetic disks (Marangoni effects); EGaIn Electrode to study quantum tunneling; Self-Assembly of electrostatically-charged metallic spheres (electrets); Dynamically reconfigurable lens; Using computational designs of ligands for enzymes; Electrostatic self-assembly by tribocharging; Monodisperse bubble trains in microchannel systems; Inverted dripping faucet; Flames; Printing of micro-organisms to regenerate the “ink” of printing device; Using micro-organisms to move loads (“microoxen”); Motion of bacterial swarms near surfaces; Making monodisperse particles in microfluidic systems; Coding/decoding of information stored in droplet trains in microfluidic networks; Magnetic levitation; and Information storage. i) Tribocharging. The change in focus of this work on electrets from the fundamentals of charging to applications of these materials in studying self-assembly using electrostatic interactions. ii) Bubbles in Microchannels. The realization that systems of bubbles in microchannels represented a major opportunity to study complexity in a very tractable system, and the development of a semi-quantitative theory of this subject. iii) Flames. The growth of “flames” remains an exploratory subject for the research, although their currently relatively little active work involving it ongoing. iv) Systems with Microorganisms. The removal of work in biological systems from this project. Based on work supported in this program, we now have a significant project on the development of microfluidic tools for studying C. elegans (a nematode), but this work was not appropriate for a program focused on complexity, and we developed separate support for it. (It is, however, an example of successful seeding of a new area by BES.) The work on electrets has gone through a period in which a part of the program was the subject of a MURI; the focus of this work was to develop materials that did not charge electrostatically on friction or contact. The MURI is now over, and the work on dynamic self-assembly (supported by BES) is the major focus. “Flames” has also enjoyed synergistic support, in terms of a project supported by DARPA on flame suppression (in the absence of extinguishing agents, using acoustic and electrostatic interactions). This work was helpful in understanding some of the basics of flames, but is entirely distinct from the BES focus in complexity. A growing interest is in the Lorentz effect. The Lorentz effect is the force exerted on charged particles (electrons, ions, charged molecules) when they move through a perpendicular magnetic field. The Lorentz effect is almost ubiquitous in modern technology: examples of applications include electric motors, dynamos, cathode ray tubes, many batteries, and most systems that control electrical currents with magnetic forces. We have begun to explore the Lorentz effect in electrochemical systems and heterogeneous catalytic systems involving charged organic species and inorganic ions. This work is still at an early stage, but initial studies that Lorentz effects can be large when ions move through magnetic fields, or magnetic fields move in the presence of ions.

36 MATERIALS SCIENCE↗

Tuning superconductivity and spin-vortex instabilities in CaKFe 4 As 4 through in-plane antisymmetric strains

Lattice strains of appropriate symmetry have served as an excellent tool to explore the interaction of superconductivity in the iron-based superconductors with orthorhombic-nematic and stripe spin-density-wave (SSDW) order. In this Letter, we contribute to a broader understanding of the coupling of strain to superconductivity and competing normal-state orders by studying CaKFe 4 As 4 under large, in-plane strains of B 1 g and B 2 g symmetry. In contrast to the majority of iron-based superconductors, pure CaKFe 4 As 4 exhibits superconductivity with a relatively high transition temperature of T c ∼ 35 K in proximity of a noncollinear, tetragonal, hedgehog spin-vortex crystal (SVC) order. Through experiments and calculations, we demonstrate an anisotropic in-plane strain response of T c and the favored SVC configuration in CaKFe 4 As 4 . This supports a scenario, in which the change in spin fluctuations dominates the strain response of superconducting T c . Overall, by suggesting moderate B 2 g strains as an effective parameter to change the stability of SVC and SSDW, we outline a pathway to a unified phase diagram of iron-based superconductivity. Published by the American Physical Society 2024

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Electrical, Structural, and Thermal Properties of Ferrite/Superconductor (Ni{sub 0.5}Zn{sub 0.5}Fe{sub 2}O{sub 4}){sub x}/YBa{sub 2}Cu{sub 3}O{sub 7-δ}) Nanocomposite Materials

The superconductor YBa{sub 2}Cu{sub 3}O{sub 7-δ} (Y123) is prepared by the solid-state route, while the nanosized Ni{sub 0.5}Zn{sub 0.5}Fe{sub 2}O{sub 4} (Ni-Zn ferrite) is synthesized via citric acid sol-gel technique. The composite material of the type (Ni{sub 0.5}Zn{sub 0.5}Fe{sub 2}O{sub 4}){sub x}/YBa{sub 2}Cu{sub 3}O{sub 7-δ}, where x = 0.00, 0.03, 0.10, and 0.50 wt%, is prepared by a modified solid-state method to study the effect of doping Ni-Zn nanosized ferrite on the superconducting, structure, and thermal properties of Y123. The prepared samples show Meissner effect which confirmed the presence of superconducting phase at liquid nitrogen temperature (77 K). The prepared samples have been characterized by X-ray diffraction (XRD), scanning electron microscopy (SEM), energy dispersion of X-ray (EDX), infra-red (IR) spectroscopy, thermogravimetric analysis (TGA), vibrating sample magnetometer (VSM), DC resistance by four-probe, and critical current density measurements. It was found that by increasing Ni{sub 0.5}Zn{sub 0.5}Fe{sub 2}O{sub 4} nanoparticles in the Y123 compound, the unit cell volume of the orthorhombic phase decreases, but the transition temperature T{sub c} and critical current density J{sub c} decrease for low doping (x = 0.03 wt%) then increase (x = 0.10 and 0.50 wt%) for Ni{sub 0.5}Zn{sub 0.5}Fe{sub 2}O{sub 4} doped in Y123. The transport properties of the (Ni{sub 0.5}Zn{sub 0.5}Fe{sub 2}O{sub 4}){sub x}/YBa{sub 2}Cu{sub 3}O{sub 7}-{sub δ} composite as electrical resistance R, T{sub c}, and J{sub c} are interpreted according to Bean’s critical state model and charge-vortex interaction combined with vortex pinning in Y123 high-T{sub c} mixed-state type-II superconductor. Highlights: The critical current density J{sub c} decreased for low doping of nanosized Ni-Zn ferrite into Y123 with ratio of 0.03 wt% and increased by increasing doping ratio to 0.50 wt%. The physical properties are improved by the coexistence of the superconductor Y123 and nanosized Ni-Zn ferrites in nanocomposite material.

77 NANOSCIENCE AND NANOTECHNOLOGY↗

Structural, magnetic, and superconducting properties of pulsed-laser-deposition-grown La 1.85 Sr 0.15 CuO 4 / La 2 / 3 Ca 1 / 3 MnO 3 superlattices on (001)-oriented LaSrAlO 4 substrates

Epitaxial La 1.85 Sr 0.15 CuO 4 / La 2 / 3 Ca 1 / 3 MnO 3 (LSCO/LCMO) superlattices (SL) on (001)- oriented LaSrAlO 4 substrates have been grown with pulsed laser deposition (PLD) technique. Their structural, magnetic and superconducting properties have been determined with in-situ reflection high energy electron diffraction (RHEED), x-ray diffraction, specular neutron reflectometry, scanning transmission electron microscopy (STEM), electric transport, and magnetization measurements. We find that despite the large mismatch between the in-plane lattice parameters of LSCO (a = 0.3779 nm) and LCMO (a = 0.387 nm) these superlattices can be grown epitaxially and with a high crystalline quality. While the first LSCO layer remains clamped to the LSAO substrate, a sizeable strain relaxation occurs already in the first LCMO layer. The following LSCO and LCMO layers adopt a nearly balanced state in which the tensile and compressive strain effects yield alternating in-plane lattice parameters with an almost constant average value. No major defects are observed in the LSCO layers, while a significant number of vertical antiphase boundaries are found in the LCMO layers. The LSCO layers remain superconducting with a relatively high superconducting onset temperature of T c onset ≈ 36 K. The macroscopic superconducting response is also evident in the magnetization data due to a weak diamagnetic signal below 10 K for H ∥ ab and a sizeable paramagnetic shift for H ∥ c that can be explained in terms of a vortex-pinning-induced flux compression. The LCMO layers maintain a strongly ferromagnetic state with a Curie temperature of T Curie ≈ 190 K and a large low-temperature saturation moment of about 3.5 (1) μ B . These results suggest that the LSCO/LCMO superlattices can be used to study the interaction between the antagonistic ferromagnetic and superconducting orders and, in combination with previous studies on YBCO/LCMO superlattices, may allow one to identify the relevant mechanisms.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Numerical Modeling of a Vortex Stabilized Arcjet

Arcjet thrusters are being actively considered for use in Earth orbit maneuvering applications. Experimental studies are currently the chief means of determining an optimal thruster configuration. Earlier numerical studies have failed to include all of the effects found in typical arcjets including complex geometries, viscosity, and swirling flow. Arcjet geometries are large area ratio converging nozzles with centerbodies in the subsonic portion of the nozzle. The nozzle walls serve as the anode while the centerbody functions as the cathode. Viscous effects are important because the Reynolds number, based on the throat radius, is typically less than 1,000. Experimental studies have shown that a swirl or circumferential velocity component stabilizes a constricted arc. This dissertation describes the equations governing flow through a constricted arcjet thruster. An assumption that the flowfield is in local thermodynamic equilibrium leads to a single fluid plasma temperature model. An order of magnitude analysis reveals the governing fluid mechanics equations are uncoupled from the electromagnetic field equations. A numerical method is developed to solve the governing fluid mechanics equations, the Thin Layer Navier-Stokes equations. A coordinate transformation is employed in deriving the governing equations to simplify the application of boundary conditions in complex geometries. An axisymmetric formulation is employed to include the swirl velocity component as well as the axial and radial velocity components. The numerical method is an implicit finite-volume technique and allows for large time steps to reach a converged steady-state solution. The inviscid fluxes are flux-split, and Gauss-Seidel line relaxation is used to accelerate convergence. Converging-diverging nozzles with exit-to-throat area ratios up to 100:1 and annular nozzles were examined. Quantities examined included Mach number and static wall pressure distributions, and oblique shock structures. As the level of swirl and viscosity in the flowfield increased the mass flow rate and thrust decreased. The technique was used to predict the flow through a typical arcjet thruster geometry. Results indicate swirl and viscosity play an important role in the complex geometry of an arcjet.

Pawlas, Gary E.↗

Dynamics of an inertially collapsing gas bubble between two parallel, rigid walls

The collapse of cavitation bubbles in channel flows can give rise to structural damage along neighbouring walls. Although the collapse of a bubble near a single wall has been studied extensively, less is known about bubble collapse between two walls, e.g. as in a channel. We conduct highly resolved, direct simulations of the Navier–Stokes equations to investigate the bubble dynamics and pressures produced by the collapse of a bubble between two parallel rigid walls. We examine the dependence of the dynamics and pressures on the initial bubble location, confinement and driving pressure. For a fixed initial stand-off distance, as the channel width increases the bubble volume, migration distance and re-entrant jet speed approach their single-wall counterparts. We obtain an expression for the minimum channel width at which the confinement does not affect the bubble dynamics depending on the driving pressure difference and initial stand-off distance. For a fixed channel width, varying stand-off distance reduced the maximum wall pressures in the channel relative to the single wall; the trend was consistent for three different driving pressures. Two different jetting behaviours are seen when the bubble is centred in the channel, depending on the channel width. Under significant confinement, wall-parallel re-entrant jets impinge upon each other and further intensify the collapse of the vortex ring.

Rodriguez, Jr, Mauro (ORCID:0000000305450265)↗

Flow Physics at a Large Rotor Tip Gap in a Multi-Stage Axial Compressor

The flow physics in a large rotor tip gap in a 1.5-stage axial compressor is investigated in the current study. The flow structure in the rotor tip region is complex with several dominant vortical structures of opposite rotation, resulting in inhomogeneous and highly anisotropic turbulence. Earlier measurements show that eddy viscosity is negative over large parts of the tip region and eddy viscosity varies among stress/strain components. The present study aims to understand how the complex nature of rotor tip leakage flow affects compressor performance when the tip gap size is greater than 4–5% of the rotor span, which is typical of advanced small core engines. Unsteady Reynolds-averaged Navier-Stokes (URANS) and Large Eddy Simulation (LES) techniques are applied to study flow physics in a large rotor tip gap (5.5% of rotor span) in a 1.5-stage axial compressor. Calculated flow fields from the two different approaches are compared with available measurements and examined in detail. LES calculates the pressure rise in the present compressor fairly well, while URANS with a standard two-equation turbulence closure underpredicts the pressure rise by 15-20% of the measured values. The current study shows that URANS with the current turbulence closure produces much higher all-positive eddy viscosity in the tip-gap region compared to measurements and LES. The distribution of eddy viscosity in the URANS simulation is also wrong. Consequently, the flow in the tip region is highly damped with significantly larger blockage generation, which results in the tip leakage vortex (TLV) staying closer to the blade suction side compared to the measurement. When the TLV stays closer to the blade, both flow turning and the pressure rise across the compressor are reduced compared to the measurements. It appears that this effect is amplified by a large rotor tip gap.

tip gap↗

Electron Tunneling and X-Ray Photoelectron Spectroscopy Studies of the Superconducting Properties of Nitrogen-Doped Niobium Resonator Cavities

Scanning tunneling microscopy and spectroscopy and X-ray photoelectron spectroscopy (XPS) have been used to investigate the differences between the surface electronic structures and chemical structure of typically prepared and N-doped Nb cutouts from superconducting radio frequency (SRF) cavities. The goal of this work is to get insights into the fundamental physics and materials mechanisms behind the striking decrease of the surface resistance with the radiofrequency magnetic field, which has been observed by many groups on N-doped Nb cavities. In particular, we address the effects of N-doping on the superconducting properties at the surface of SRF cavities. XPS measurements reveal a significantly more oxidized Nb 3d states on the N-doped Nb surfaces which is confirmed by tunneling spectroscopy measurements. Analysis of the tunneling spectra in the framework of a recent model (A. Gurevich and T. Kubo Phys. Rev. B 96, 184515 (2017)) shows that the N-doping greatly reduces local distribution of superconducting properties on the surface, causes a significant shrinkage of the metallic suboxide and changes the contact resistance between the metallic suboxide and the bulk niobium toward an optimum value, resulting in a lower surface resistance. Combination of these factors enables one to effectively tune the density of states at the surface and to reveal the decrease of the surface resistance with the radio frequency field which follows from the BCS theory. A slightly reduced average gap and a smaller coherence length, revealed by tunneling spectroscopy in the vortex cores, have been found in the N-doped Nb samples compared to typically prepared Nb samples, indicating a stronger impurity scattering caused by nitrogen doping in a moderately disordered material.

Lechner, Eric M.↗

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗

Multiresolution With Super-Compact Wavelets

The solution data computed from large scale simulations are sometimes too big for main memory, for local disks, and possibly even for a remote storage disk, creating tremendous processing time as well as technical difficulties in analyzing the data. The excessive storage demands a corresponding huge penalty in I/O time, rendering time and transmission time between different computer systems. In this paper, a multiresolution scheme is proposed to compress field simulation or experimental data without much loss of important information in the representation. Originally, the wavelet based multiresolution scheme was introduced in image processing, for the purposes of data compression and feature extraction. Unlike photographic image data which has rather simple settings, computational field simulation data needs more careful treatment in applying the multiresolution technique. While the image data sits on a regular spaced grid, the simulation data usually resides on a structured curvilinear grid or unstructured grid. In addition to the irregularity in grid spacing, the other difficulty is that the solutions consist of vectors instead of scalar values. The data characteristics demand more restrictive conditions. In general, the photographic images have very little inherent smoothness with discontinuities almost everywhere. On the other hand, the numerical solutions have smoothness almost everywhere and discontinuities in local areas (shock, vortices, and shear layers). The wavelet bases should be amenable to the solution of the problem at hand and applicable to constraints such as numerical accuracy and boundary conditions. In choosing a suitable wavelet basis for simulation data among a variety of wavelet families, the supercompact wavelets designed by Beam and Warming provide one of the most effective multiresolution schemes. Supercompact multi-wavelets retain the compactness of Haar wavelets, are piecewise polynomial and orthogonal, and can have arbitrary order of approximation. The advantages of the multiresolution algorithm are that no special treatment is required at the boundaries of the interval, and that the application to functions which are only piecewise continuous (internal boundaries) can be efficiently implemented. In this presentation, Beam's supercompact wavelets are generalized to higher dimensions using multidimensional scaling and wavelet functions rather than alternating the directions as in the 1D version. As a demonstration of actual 3D data compression, supercompact wavelet transforms are applied to a 3D data set for wing tip vortex flow solutions (2.5 million grid points). It is shown that high data compression ratio can be achieved (around 50:1 ratio) in both vector and scalar data set.

Lee, Dohyung↗

Progress Toward Generation of a Navier-Stokes Database for a Harrier in Ground Effect

The Harrier YAV-8B aircraft is capable of vertical and short-field take-off and landing (V/STOL) by directing its four exhaust nozzles toward the ground, or conventional flight by rotating its nozzles into a horizontal position. The British Royal Air Force and the United States Marine Corps have used this aircraft for more than 30 years to provide a quick reaction time for troop support, and reduce the need for long runways. The success of this powered-lift (PL) vehicle has also prompted the more recent design of the Joint Strike Fighter (JSF). However there are significant safety issues that must be addressed when operating a PL vehicle in close proximity to the ground. Hot Gas Ingestion (HGI) by the inlets can result in a rapid loss of powered lift; and high-speed jet flows along the ground plane can induce low pressures underneath the vehicle, causing a 'suck-down' effect. Under these conditions, departure from controlled flight may occur. Moreover, unsteady ground vortices and jet fountains can affect the aircraft,s controllability and its proximity to ground troops. The viscous, time-dependent flow fields of PL vehicles are difficult to accurately and efficiently predict using Computational Fluid Dynamics (CFD). A number of researchers have used the time-dependent Reynolds-averaged Navier-Stokes (RANS) equations to compute flows for single and multiple jets in a cross-flow. A few have added some geometric complexity to the problem by computing flows for jet-augmented delta wings near a ground plane. Smith et.al. computed for the first time a single RANS solution about a simplified Harrier. This geometry included a fuselage, wing, leading edge root extension (LERX), inlets, and exhaust nozzles. All of these investigations cite two practical problems with computing these flows: 1) the need for improved solution accuracy; and, 2) the need for faster solution methods. We view the need for faster solution methods as key to improving the solution accuracy and making this class of computation more routine. One can hardly refine grids, explore the use of advanced turbulence models, and generate databases when it takes weeks of dedicated computer time for a single solution. Chaderjian, Ahmad, Pandya, and Murman have focused on reducing the time-to-solution for this very difficult and complex problem through process automation and exploitation of parallel computing. They began with the Harrier geometry reported, and added a deflected wing flap and empennage for greater realism. To date more than 80 solutions have been carried out. This paper will describe this process and progress made in reducing the time required to generate a simple longitudinal force and moment database for a Harrier in ground effect. It shows a typical snap-shot from an unsteady streakline animation, where fluid particles are colored by temperature. The ground vortex and a jet-fountain vortex are highlighted. It also shows a similar streakline image, where HGI occurs due to the vehicle in close proximity to the ground. It is show the mean lift coefficient as a function of angle of attack and height. The angle of attack range was 4 deg less than or = alpha less than or = 10 deg with an increment of 1 degree, and the height range was 10 ft less than or = h less than or = 30ft with an increment of 5 feet. This 35 solution database was extended to over 2500 cases using a monotone cubic-spline interpolation procedure. The suck-down effect (reduction of lift near the ground) is highlighted in the figure. The "cushion effect," the conventional reduction of lift as the vehicle moves out of ground effect, is also indicated. All 35 RANS solutions were obtained using 952 Silicon Graphics Origin 2000 and 3000 processors in dedicated mode for one week. Typically, 112 processors were assigned to each case. Some other cases used fewer processors to utilize all available CPUS. The final paper will report on the automation of the solution process, including: grid generation, job monitoring, solution completion criteria, and post processing. Moreover, improvements in parallel efficiency for a dual time-step algorithm for the RANS equations will also be presented. Results will be discussed in detail using unsteady streakline flow visualization to correlate unsteady flow structures with dominant aerodynamic frequencies. The stability derivatives, CL, and CL, will also be presented.

Chaderjian, Neal M.↗

A comprehensive review of nonlinear oscillators in hydrokinetic energy harnessing using flow-induced vibrations

A comprehensive review of hydrokinetic energy converters based on alternating lift technology (ALT) is provided. Emphasis is on nonlinear oscillators based on Flow Induced Vibration (FIV) or Oscillation (FIO). Due to strong coupling in Fluid-Structure Interaction (FSI), and in order to maximize the hydrokinetic harnessed energy, design of nonlinear oscillators and analysis by model tests or computational fluid dynamics dominates this area. Research confirmed that the nonlinear oscillator can harvest energy from a stochastic excitation modeled by a generic wide spectrum, and overcome the most severe oscillator limitations: specifically, the need for continuous frequency tuning due to the narrow bandwidth response, and low efficiency outside the narrow bandwidth oscillator response. This review covers the following aspects of nonlinear oscillators in ALT converters: (1) Geometric changes in oscillator cross-section; e.g., circular, square, rectangular, or trilateral shapes. (2) Passive turbulence control of FIV/FIO. (3) Position based nonlinear stiffness. (4) Multi-cylinder synergistic FIV/FIO. (5) Mechanically linked oscillators. (6) Velocity-based, nonlinear, adaptive harnessing damping.

42 ENGINEERING↗

The Turbulent/Non-Turbulent Interface Bounding a Far-Wake

The velocity fields of a turbulent wake behind a flat plate obtained from the direct numerical simulations of Moser et al. are used to study the structure of the flow in the intermittent zone where there are, alternately, regions of fully turbulent flow and non-turbulent velocity fluctuations either side of a thin randomly moving interface. Comparisons are made with a wake that is 'forced' by amplifying initial velocity fluctuations. There is also a random temperature field T in the flow; T varies between constant values of 0.0 and 1.0 on the sides of the wake. The value of the Reynolds number based on the centreplane mean velocity defect and halfwidth b of the wake is Re approx. = 2000. It is found that the thickness of the continuous interface is about equal to 0.07b, whereas the amplitude of fluctuations of the instantaneous interface displacement y(sub I)(t) is an order of magnitude larger, being about 0.5b. This explains why the mean statistics of vorticity in the intermittent zone can be calculated in terms of the probability distribution of y(sub I) and the instantaneous discontinuity in vorticity across the interface. When plotted as functions of y - y(sub I), the conditional mean velocity (U) and temperature (T) profiles show sharp jumps Delta(U) and Delta(T) at the interface adjacent to a thick zone where (U) and (T) vary much more slowly. Statistics for the vorticity and velocity variances, available in such detail only from DNS data, show how streamwise and spanwise components of vorticity are generated by vortex stretching in the bulges of the interface. Flow fields around the interface, analyzed in terms of the local streamline pattern, confirm previous results that the advancement of the vortical interface into the irrotational flow is driven by large-scale eddy motion. It is argued that because this is an inviscid mechanism the entrainment process is not sensitive to the value of Re, and that small-scale nibbling only plays a subsidiary role. While mean Reynolds stresses decrease gradually in the intermittent zone, conditional stresses are found to decrease sharply towards zero at the interface. Using one-point turbulence models applied to either unconditional or conditional statistics for the turbulent region and then averaged, the entrainment rate E(sub b) would, if calculated exactly, be zero. But if computed with standard computational methods, E(sub b) would be non-zero because of numerical diffusion. It is concluded that the current practice in statistical models of approximating entrainment by a diffusion process is computationally arbitrary and physically incorrect. An analysis shows how E(sub b) is related to Delta(U) and the jump in shear stress at the interface, and correspondingly to Delta(T) and the heat flux.

Bisset, David K.↗