Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fluid Structure Interaction”

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,045 records · Page 58

Materials for Ultra‐Coherent, Mobile, Electron‐Spin Qubits

This research project has had the goal of gaining a better understanding of the physics of electrons bound to the surface of superfluid helium from both experimental and theoretical perspectives. It has particularly been aimed at two areas which had not been well studied: the relaxation and decoherence of the spin of the electrons on the helium surface and how the properties of underlying metallic layers affect the behavior of the electrons when the helium covering the metal is thin. This work is motivated in part by interest in using the spin of these electrons as a quantum bit, or qubit. Low levels of decoherence are advantageous for qubits, and moving the electrons, as one might do in a quantum processor, will be easiest if thin helium films can be employed. It had been suggested that spin decoherence should be very weak for electrons bound to superfluid He, but before this work there have been no quantitative studies of spin relaxation and decoherence. It is especially important to know how moving the electrons across the helium surface would affect their spin coherence. Calculations performed as part of this project show that the Rashba effective magnetic field, the mechanism which limits the spin coherence of mobile electrons in silicon-based devices (an actively pursued qubit technology), is exceptionally weak for electrons bound to helium. This project has identified other decoherence mechanisms which are stronger, but still weak compared to analogous silicon-based structures. Calculated spin coherence times for mobile electrons approach one day, as compared to microseconds in silicon. With coherence times of this magnitude, the spin qubit errors on helium will be completely dominated by errors in the quantum gates. In related work, the possibility of using an artificial spin-orbit interaction (a gradient magnetic field) for quantum operations on the electrons spins was considered. The calculations show that a moderate gradient field, small enough to be generated by a narrow superconducting wire, will enable high-fidelity quantum operations on electrons held in lithographically-defined quantum dots by driving them with a microwave electric field. The spin and motional coherence of the electrons is sufficient to allow high-fidelity 2-qubit quantum operations between electrons in neighboring quantum dots. As an outgrowth of experiments aiming to measure electron spin coherence it was discovered that very high densities of electrons can be stably supported on thin helium films coating ultra-smooth amorphous metallic layers. The measured densities are high enough that the electron system has almost certainly transitioned from an ordered array of electrons, known as a Wigner crystal (ordered by the electrons’ mutual repulsion), to a quantum fluid known as a Fermi liquid. This transition has been a subject of intense interest for over 40 years, since the electron Wigner crystal was first observed with electrons bound to superfluid helium, but it has never been unambiguously observed. Experiments are still underway in these new structures to definitively determine whether true quantum melting of the Wigner crystal has been demonstrated. This work has also catalyzed the development of a new approach for measuring the transport of electrons across very thin helium films, as will be needed for some of the quantum computing applications. The high electron density experiments as well as experiments with electrons bound in quantum dots have led to new techniques which may enable spin coherence measurements.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Designing CO 2 -Responsive Multifunctional Nanoscale Fluids with Tunable Hydrogel Behavior for Subsurface Energy Recovery

The ability to direct subsurface fluid flow or achieve a high level of control on permeability in subsurface environments necessitates the development of tunable novel subsurface fluids. These multifunctional fluids should have the potential to form hydrogels to divert flow or enhance fracture networks at elevated pressures, ability to transport proppants, undergo reversible transformations from gel-like to fluid-like in response to chemical- or pressure-based perturbations, and have improved CO 2 carrying capacity for enhancing the miscibility and flowability of oil and gas. In this paper, we discuss the development of multifunctional nanofluids constructed from silica (SiO 2 ) nanoparticles and poly(allylamine) (PAA), amine bearing polymer chains with high affinity for CO 2 . A 2-fold increase in CO 2 absorption in SiO 2 -PAA nanofluids compared to the pure polymer was noted. Upon CO 2 absorption, the weakly interacting polymeric chains around the nanoparticles formed relatively compact hydrogels. Time-resolved ultrasmall-angle X-ray scattering/small-angle X-ray scattering measurements showed the transition from swollen branched polymers to Gaussian coils with an increased exposure to CO 2 . Further, CO 2 -induced hydrogel formation in aqueous fluids bearing 1 wt % SiO 2 -PAA nanofluids occurred at room temperature, unlike in fluids bearing 1 wt % PAA. These observations point to the feasibility of forming hydrogels at lower temperatures and pressures using novel nanofluids as opposed to using the pure polymer. As a result, the ability to tune the structures and morphologies of these fluids expands the potential applications of these fluids in a wide range of subsurface environments.

03 NATURAL GAS↗

Opening the reheating box in multifield inflation

The robustness of multi-field inflation to the physics of reheating is investigated. In order to carry out this study, reheating is described in detail by means of a formalism which tracks the evolution of scalar fields and perfect fluids in interaction (the inflatons and their decay products). This framework is then used to establish the general equations of motion of the background and perturbative quantities controlling the evolution of the system during reheating. Next, these equations are solved exactly by means of a new numerical code. Moreover, new analytical techniques, allowing us to interpret and approximate these solutions, are developed. As an illustration of a physical prediction that could be affected by the micro-physics of reheating, the amplitude of non-adiabatic perturbations in double inflation is considered. It is found that ignoring the fine-structure of reheating, as usually done in the standard approach, can lead to differences as big as ~ 50%, while our semi-analytic estimates can reduce this error to ~ 10%. Finally, we conclude that, in multi-field inflation, tracking the perturbations through the details of the reheating process is important and, to achieve good precision, requires the use of numerical calculations.

79 ASTRONOMY AND ASTROPHYSICS↗

Aeroelastic Calculations of Quiet High- Speed Fan Performed

An advanced high-speed fan was recently designed under a cooperative effort between the NASA Glenn Research Center and Honeywell Engines & Systems. The principal design goals were to improve performance and to reduce fan noise at takeoff. Scale models of the Quiet High-Speed Fan were tested for operability, performance, and acoustics. During testing, the fan showed significantly improved noise characteristics, but a self-excited aeroelastic vibration known as flutter was encountered in the operating range. Flutter calculations were carried out for the Quiet High-Speed Fan using a three-dimensional, unsteady aerodynamic, Reynolds-averaged Navier-Stokes turbomachinery code named "TURBO." The TURBO code can accurately model the viscous flow effects that can play an important role in various aeroelastic problems such as flutter with flow separation, flutter at high loading conditions near the stall line (stall flutter), and flutter in the presence of shock and boundary-layer interaction. Initially, calculations were performed with no blade vibrations. These calculations were at a constant rotational speed and a varying mass flow rate. The mass flow rate was varied by changing the backpressure at the exit boundary of the computational domain. These initial steady calculations were followed by aeroelastic calculations in which the blades were prescribed to vibrate harmonically in a natural mode, at a natural frequency, and with a fixed interblade phase angle between adjacent blades. The AE-prep preprocessor was used to interpolate the in-vacuum mode shapes from the structural dynamics mesh onto the computational fluid dynamics mesh and to smoothly propagate the grid deformations from the blade surface to the interior points of the grid. The aeroelastic calculations provided the unsteady aerodynamic forces on the blade surface due to blade vibrations. These forces were vector multiplied with the structural dynamic mode shape to calculate the work done on the blade during one vibration period, then this result was converted to an aerodynamic damping. Flutter occurs when the aerodynamic damping becomes negative, if structural damping is ignored. The results of these aeroelastic calculations are summarized in this plot of aerodynamic damping versus mass flow rate at a constant rotational speed. As the backpressure is increased, the mass flow rate through the fan decreases and the fan operating point moves towards the stall line. The aeroelastic calculations showed that the aerodynamic damping decreases as the stall line is approached, as observed during testing. In addition, the aeroelastic calculations with the TURBO code correctly predicted the aeroelastic parameters: the most unstable vibration mode and interblade phase angle, as observed during testing. The Quiet High-Speed Fan demonstrated significant noise reductions during testing, but flutter imposed limits on its operating range. The accurate calculation of the aeroelastic characteristics using the TURBO code is a significant step toward eliminating flutter from the operating range and toward realizing the benefits of reduced fan noise. The aeroelastic calculations described here were performed under a grant by University of Toledo researchers in collaboration with Glenn's researchers.

Bakhle, Milind A.↗

Bottlebrush Polymers at Liquid Interfaces: Assembly Dynamics, Mechanical Properties, and All-Liquid Printed Constructs

Bottlebrush polymer surfactants (BPSs), formed by the interfacial interactions between bottlebrush polymers (BPs) with poly(acrylic acid) side chains dissolved in an aqueous phase and amine-functionalized ligands dissolved in the oil phase, assemble and bind strongly to the fluid-fluid interface. The ratio between N BB (backbone degree of polymerization) and N SC (side chain degree of polymerization) defines the initial assembly kinetics, interface packing efficiency, and stress relaxation. The equilibrium interfacial tension (γ) increases when N BB < N SC , but decreases when N BB $\gg$ N SC , correlating to a pronounced change in the effective shape of the BPs from being spherical to worm-like structures. The apparent surface coverage (ASC), i.e., the interfacial packing efficiency, decreases as N BB increases. The dripping-to-jetting transition of an injected polymer solution, as well as fluorescence recovery after photobleaching experiments, revealed faster initial assembly kinetics for BPs with higher N BB . Furhter, Euler buckling of BPS assemblies with different N BB values was used to characterize the stress relaxation behavior and bending modulus. The stress relaxation behavior was directly related to the ASC, reflecting the strong influence of macromolecular shape on packing efficiency. The bending modulus of BPSs decreases for N BB < N SC , but increased when N BB $\gg$ N SC , showing the effect of molecular architecture and multisite anchoring. All-liquid printed constructs with lower N BB BPs yielded more stable structured liquids, underscoring the importance of macromolecular packing efficiency at fluid interfaces. Overall, this work elucidates fundamental relationships between nanoscopic structures and macroscopic properties associated with various bottlebrush polymer architectures, which translate to the stabilization of all-fluidic printed constructs.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On oscillatory flow over topography in a rotating fluid

Quasi-geostrophic beta-plane motion of a homogeneous liquid over topography is considered for situations in which there is a time-periodic forcing of a zonal current. Such an oscillatory current generates a topographic Rossby wave response that has a complicated, but periodic, temporal structure. The linear solution grows resonances at all integer values of the beta-parameter. The nonlinear analysis demonstrates that for weak friction and forcing, the resonances are bent and multiple equilibria of the time-dependent Rossby wave states are possible in certain parameter ranges. While the basic forced flow in the absence of topography has no time-mean, the nonlinear amplitude equations show that a mean retrograde (westward) Eulerian zonal flow is generated in the interactions of the forced flow with the topography.

Hart, J. E.↗

Fluid Aspects of Solar Wind Disturbances Driven by Coronal Mass Ejections

Transient disturbances in the solar wind initiated by coronal eruptions have been modeled for many years, beginning with the self-similar analytical models of Parker and Simon and Axford. The first numerical computer code (one-dimensional, gas dynamic) to study disturbance propagation in the solar wind was developed in the late 1960s, and a variety of other codes ranging from simple one-dimensional gas dynamic codes through three-dimensional gas dynamic and magnetohydrodynamic codes have been developed in subsequent years. For the most part, these codes have been applied to the problem of disturbances driven by fast CMEs propagating into a structureless solar wind. Pizzo provided an excellent summary of the level of understanding achieved from such simulation studies through about 1984, and other reviews have subsequently become available. More recently, some attention has been focused on disturbances generated by slow CMEs, on disturbances driven by CMEs having high internal pressures, and disturbance propagation effects associated with a structured ambient solar wind. Our purpose here is to provide a brief tutorial on fluid aspects of solar wind disturbances derived from numerical gas dynamic simulations. For the most part we illustrate disturbance evolution by propagating idealized perturbations, mimicking different types of CMEs, into a structureless solar wind using a simple one-dimensional, adiabatic (except at shocks), gas dynamic code. The simulations begin outside the critical point where the solar wind becomes supersonic and thus do not address questions of how the CMEs themselves are initiated. Limited to one dimension (the radial direction), the simulation code predicts too strong an interaction between newly ejected solar material and the ambient wind because it neglects azimuthal and meridional motions of the plasma that help relieve pressure stresses. Moreover, the code ignores magnetic forces and thus also underestimates the speed with which pressure disturbances propagate in the wind.

Gosling, J. T.↗

Facile polymerization in a bicellar template to produce polymer nano-rings

We report a well-defined discoidal bicelle composed of three lipids, specifically zwitterionic long-chain 1,2-dipalmitoyl phosphocholine (DPPC) and short-chain 1,2-dihexanoyl phosphocholine (DHPC) doped with anionic 1,2-dipalmitoyl phosphoglycerol (DPPG) provides a generalized template for the synthesis of hydrophobic polymer nano-rings. The lipid molar ratio of DPPC/DHPC/DPPG is 0.71/0.25/0.04. The detailed investigation and discussion were based on styrene but tested on three other vinyl monomers. The structure of nano-rings is identified through the detailed analysis of small angle X-ray/neutron scattering (SAXS and SANS) data and transmission electron micrographs (TEM), supported by the differential scanning calorimetric (DSC) data before and after polymerization. The investigation covers samples with a styrene-to-lipid ratio ranged varied from 1:50 to 1:10.The styrene monomers are initially located at both the discoidal planar (long-chain lipid rich) and rim (short-chain lipid rich) regions. During polymerization, they migrate to the more fluid rim region-section. The formation mechanism involves the interplay of hydrophobic interaction, mismatched miscibility of polystyrene between the ordered and disordered phases, and crystallinity of the long lipid acyl chains. This facile synthesis is proven applicable for several hydrophobic monomers. The well-defined nano-rings greatly enhance the interfacial area and have the potential to be the building blocks for functional materials, if monomers are incorporated with desirable functions, for future applications. Crosslinked polymer nano-rings are produced in well-defined discoidal bicelles. The enhanced immiscibility in the planar region drives the hydrophobic polymers to the fluidic rim. The platform can be generalized for synthesis of other hydrophobic polymer nano-rings.

36 MATERIALS SCIENCE↗

Multi‐scale simulation of rock compaction through breakage models with microstructure evolution

Abstract. Regional subsidence due to fluid depletion includes the interaction among multiple physical processes. Specifically, rock compaction is governed by coupled hydro-mechanical feedbacks involving fluid flow, effective stress change and pore collapse. Although poroelastic models are often used to explain the delay between depletion and subsidence, recent evidence indicates that inelastic effects could alter the rock microstructure, thus exacerbating coupling effects. Here, a constitutive law built within the framework of Breakage Mechanics is proposed to account for the inherent connection between rock microstructure, hydraulic conductivity, and pore compaction. Furthermore, it is embedded into a 1-D hydromechanical coupled finite element analysis (FEA) to explore the effects of micro-structure rearrangement on the development of reservoir compaction. Numerical examples with the proposed model are compared with simulations under constant hydraulic conductivity to illustrate the model capability to capture the non-linear processes of reservoir compaction induced by fluid depletion.

Buscarnera, Giuseppe↗

Extension of HCDstruct for Transonic Aeroservoelastic Analysis of Unconventional Aircraft Concepts

A substantial effort has been made to implement an enhanced aerodynamic modeling capability in the Higher-fidelity Conceptual Design and structural optimization tool. This additional capability is needed for a rapid, physics-based method of modeling advanced aircraft concepts at risk of structural failure due to dynamic aeroelastic instabilities. To adequately predict these instabilities, in particular for transonic applications, a generalized aerodynamic matching algorithm was implemented to correct the doublet-lattice model available in Nastran using solution data from a priori computational fluid dynamics anal- ysis. This new capability is demonstrated for two tube-and-wing aircraft configurations, including a Boeing 737-200 for implementation validation and the NASA D8 as a first use case. Results validate the current implementation of the aerodynamic matching utility and demonstrate the importance of using such a method for aircraft configurations featuring fuselage-wing aerodynamic interaction.

Quinlan, Jesse R.↗

An Oculometric Standard to Assess the Performance of the Ocular System for Long Duration Spaceflight

Changes in the human brain, due to spaceflight, have been a challenge to the space program since its inception; a challenge that is becoming more acute as extended International Space Station (ISS) missions become routine, a return to crewed lunar exploration is about to begin, and a multi-year crewed trip to Mars is being planned. Determining the extent to which neurophysiological adaptations may adversely impact performance in operational tasks, assessing the full-time course of these changes, and identifying/mitigating any potential long-term health consequences are crucial steps required to enable safe crew-autonomous, long-duration, deep-space missions. Vision is the predominant perceptual sense used to guide cognition and motor control in humans. Thus, it is critical for the success of any future, long-duration mission to understand how long-term exposure to microgravity and to other stressors of spaceflight and their interactions, impact visual function. This is an especially challenging task given the small and disparate samples from which we have been and will be able to collect human performance data. NASA has recently been carefully tracking spaceflight-induced ophthalmic changes, driven in large part by the fluid shifts related to microgravity. In particular, a recent study of Optical Coherence Tomography (OCT) postflight measures revealed that more than two-thirds of US crew members experience significant increases in retinal thickness after 6month or longer ISS missions. Space Associated Neuro-ocular Syndrome (SANS) also includes visual acuity decrements and other ocular structural changes that could adversely impact in-flight performance as missions become longer. However, the impact of spaceflight on the visual system is not limited to the retina. Recent comparisons of pre-and post-flight brain images have revealed structural changes throughout the brain that implicate visual, visuomotor, and visual-cognitive pathways, consistent with observed functional impacts (e.g., decreased speed/accuracy/timing of fine goal-oriented movements). Current limitations on inflight testing make it very difficult to determine when and under what conditions SANS and other disruptions of human neurological subsystems arise. Currently, one cannot anticipate the time course and extent of impairment and recovery. This study addresses this knowledge gap by creating an integrated framework of the relevant existing literature on visual and oculomotor function during spaceflight with an eye towards complementing current structural measures of visuomotor impairment with new oculometric standard measures. Using eye-movement based visual function assessment and diagnostics would provide a valuable enhancement of crew health and performance monitoring in support of deepspace exploration.

Vision↗

Morphology of shocked lateral outflows in colliding hydrodynamic flows

Supersonic interacting flows occurring in phenomena, such as protostellar jets, give rise to strong shocks and have been demonstrated in several laboratory experiments. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in three dimensions. Here, we introduce variations in the flow parameters of density, velocity, and cross-sectional radius of the colliding flows in order to study the propagation and conical shape of the bow shock formed by collisions between two, not necessarily symmetric, hypersonic flows. We find that the motion of the interaction region is driven by imbalances in ram pressure between the two flows, while the conical structure of the bow shock is a result of shocked lateral outflows being deflected from the horizontal when the flows are of differing cross sections.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Hyporheic Exchange Driven by Submerged Rigid Vegetation: A Modeling Study

Flow-vegetation interaction affects fluid flow hydraulics and associated material transport in river corridors. Concomitant changes in pressure within the flow field due to the presence of vegetation may act as a driver for the formation of hyporheic flow across the sediment-water interface. This potentially important process, however, has yet to be studied. In order to investigate vegetation-induced hyporheic exchange, a series of numerical models of interlinked surface-subsurface flow modified by plant stems was conducted. Periodically staggered plant stem arrays on a flat sediment bed were considered within a coupled multiphysics computational fluid dynamics approach. Plants were idealized as rigid cylinders and arranged in different streamwise and spanwise spacing distances. Each vegetation array was then subjected to a broad range of flow Reynolds Numbers (Re). The results showed that hyporheic flow occurs in all conditions with the presence of vegetation. The vegetation-induced hyporheic flux is found to be a function of Re via a power law. The flux increases with interstem space until the space reaches the distance that rigid stems no longer affect the flow structures in the vicinity of each other. Larger intervegetation distances lead to a larger hyporheic zone. A direct comparison with bedforminduced hyporheic flow showed that vegetation can induce higher hyporheic flux through relatively shallower exchange zones. The results of all the simulations were synthesized into predictive models for hyporheic flux, bulk residence time and exchange depth based on drag coefficient, vegetation density, and Reynolds Number.

54 ENVIRONMENTAL SCIENCES↗

3D Simulations and MLT. I. Renzini’s Critique

Renzini wrote an influential critique of "overshooting" in mixing-length theory (MLT), as used in stellar evolution codes, and concluded that three-dimensional fluid dynamical simulations were needed. Such simulations are now well tested. Implicit large eddy simulations connect large-scale stellar flow to a turbulent cascade at the grid scale, and allow the simulation of turbulent boundary layers, with essentially no assumptions regarding flow except the number of computational cells. Buoyant driving balances turbulent dissipation for weak stratification, as in MLT, but with the dissipation length replacing the mixing length. The turbulent kinetic energy in our computational domain shows steady pulses after 30 turnovers, with no discernible diminution; these are caused by the necessary lag in turbulent dissipation behind acceleration. Interactions between coherent turbulent structures give multi-modal behavior, which drives intermittency and fluctuations. These cause mixing, which may justify use of the instability criterion of Schwarzschild rather than the Ledoux. Chaotic shear flow of turning material at convective boundaries causes instabilities that generate waves and sculpt the composition gradients and boundary layer structures. The flow is not anelastic; wave generation is necessary at boundaries. A self-consistent approach to boundary layers can remove the need for ad hoc procedures of "convective overshooting" and "semi-convection." In Paper II, we quantify the adequacy of our numerical resolution in a novel way, determine the length scale of dissipation—the "mixing length"—without astronomical calibration, quantify agreement with the four-fifths law of Kolmogorov for weak stratification, and deal with strong stratification.

79 ASTRONOMY AND ASTROPHYSICS↗

Neutrino interactions in the late universe

The cosmic neutrino background is both a dramatic prediction of the hot Big Bang and a compelling target for current and future observations. The impact of relativistic neutrinos in the early universe has been observed at high significance in a number of cosmological probes. In addition, the non-zero mass of neutrinos alters the growth of structure at late times, and this signature is a target for a number of upcoming surveys. These measurements are sensitive to the physics of the neutrino and could be used to probe physics beyond the standard model in the neutrino sector. We explore an intriguing possibility where light right-handed neutrinos are coupled to all, or a fraction of, the dark matter through a mediator. In a wide range of parameter space, this interaction only becomes important at late times and is uniquely probed by late-time cosmological observables. Due to this coupling, the dark matter and neutrinos behave as a single fluid with a non-trivial sound speed, leading to a suppression of power on small scales. In current and near-term cosmological surveys, this signature is equivalent to an increase in the sum of the neutrino masses. Given current limits, we show that at most 0.5% of the dark matter could be coupled to neutrinos in this way.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Origins of Clustering of Metalate–Extractant Complexes in Liquid–Liquid Extraction

Effective and energy-efficient separation of precious and rare metals is very important for a variety of advanced technologies. Liquid-liquid extraction (LLE) is a relatively less energy intensive separation technique, widely used in separation of lanthanides, actinides, and platinum group metals (PGMs). In LLE, the distribution of an ion between an aqueous phase and an organic phase is determined by enthalpic (coordination interactions) and entropic (fluid reorganization) contributions. The molecular scale details of these contributions are not well understood. Preferential extraction of an ion from the aqueous phase is usually correlated with the resulting fluid organization in the organic phase, as the longer-range organization increases with metal loading. However, it is difficult to determine the extent to which organic phase fluid organization causes, or is caused by, metal loading. In this study, we demonstrate that two systems with the same metal loading may impart very different organic phase organizations and investigate the underlying molecular scale mechanism. Small-angle X-ray scattering shows that the structure of a quaternary ammonium extractant solution in toluene is affected differently by the extraction of two metalates (octahedral PtCl 6 2- and square-planar PtCl 6 2- ), although both are completely transferred into the organic phase. The aggregates formed by the metalate-extractant complexes (approximated as reverse micelles) exhibit a more long-range order (clustering) with PtCl 6 2- compared to that with PtCl 6 2- . Vibrational sum frequency generation spectroscopy and complementary atomistic molecular dynamics simulations on model Langmuir monolayers indicate that the two metalates affect the interfacial hydration structures differently. Furthermore, the interfacial hydration is correlated with water extraction into the organic phase. Overall, these results support a strong relationship between the organic phase organizational structure and the different local hydration present within the aggregates of metalate-extractant complexes, which is independent of metalate concentration.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NON-BAROTROPIC LINEAR ROSSBY WAVE INSTABILITY IN THREE-DIMENSIONAL DISKS

Astrophysical disks with localized radial structure, such as protoplanetary disks containing dead zones or gaps due to disk-planet interaction, may be subject to the non-axisymmetric Rossby wave instability (RWI) that leads to vortex formation. The linear instability has recently been demonstrated in three-dimensional (3D) barotropic disks. It is the purpose of this study to generalize the 3D linear problem to include an energy equation, thereby accounting for baroclinity in three dimensions. Linear stability calculations are presented for radially structured, vertically stratified, geometrically thin disks with non-uniform entropy distribution in both directions. Polytropic equilibria are considered but adiabatic perturbations assumed. The unperturbed disk has a localized radial density bump, making it susceptible to the RWI. The linearized fluid equations are solved numerically as a partial differential equation eigenvalue problem. Emphasis on the ease of method implementation is given. It is found that when the polytropic index is fixed and adiabatic index increased, non-uniform entropy has negligible effect on the RWI growth rate, but pressure and density perturbation magnitudes near a pressure enhancement increase away from the midplane. The associated meridional flow is also qualitatively changed from homentropic calculations. Meridional vortical motion is identified in the nonhomentropic linear solution, as well as in a nonlinear global hydrodynamic simulation of the RWI in an initially isothermal disk evolved adiabatically. Numerical results suggest that buoyancy forces play an important role in the internal flow of Rossby vortices.

79 ASTRONOMY AND ASTROPHYSICS↗

Reaction enhancement of initially distant scalars by Lagrangian coherent structures

Turbulent fluid flows have long been recognized as a superior means of diluting initial concentrations of scalars due to rapid stirring. Conversely, experiments have shown that the structures responsible for this rapid dilution can also aggregate initially distant reactive scalars and thereby greatly enhance reaction rates. Indeed, chaotic flows not only enhance dilution by shearing and stretching but also organize initially distant scalars along transiently attracting regions in the flow. To show the robustness of this phenomenon, a hierarchical set of three numerical flows is used: the periodic wake downstream of a stationary cylinder, a chaotic double gyre flow, and a chaotic, aperiodic flow consisting of interacting Taylor vortices. We demonstrate that Lagrangian coherent structures (LCS), as identified by ridges in finite time Lyapunov exponents, are directly responsible for this coalescence of reactive scalar filaments. When highly concentrated filaments coalesce, reaction rates can be orders of magnitude greater than would be predicted in a well-mixed system. This is further supported by an idealized, analytical model that was developed to quantify the competing effects of scalar dilution and coalescence. Chaotic flows, known for their ability to efficiently dilute scalars, therefore have the competing effect of organizing initially distant scalars along the LCS at timescales shorter than that required for dilution, resulting in reaction enhancement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗