Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “FREE OSCILLATION”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Dynamic Nucleation of Supercooled Melts and Measurement of the Surface Tension and Viscosity

We investigate the phenomenon of acoustic pressure-induced nucleation by using a novel approach involving the large amplitude resonant radial oscillations and collapse of a single bubble intentionally injected into a supercooled liquid. Using a combination of previously developed and proven techniques, the bubble is suspended in a fluid host by an ultrasonic field which supplies both the levitation capability as well as the forcing of the radial oscillations. We observe the effects of an increase in pressure (due to bubble collapse) in a region no larger than 100 microns within the supercooled melt to rigorously probe the hypothesis of pressure-induced nucleation of the solid phase. The use of single bubbles operating in narrow temporal and spatial scales will allow the direct and unambiguous correlation between the origin and location of the generation of the disturbance and the location and timing of the nucleation event. In a companion research effort, we are developing novel techniques for the non-contact measurements of the surface tension and viscosity of highly viscous supercooled liquids. Currently used non-invasive methods of surface tension measurement for the case of undercooled liquids generally rely of the quantitative determination of the resonance frequencies of drop shape oscillations, of the dynamics of surface capillary waves, or of the velocity of streaming flows. These methods become quickly ineffective when the liquid viscosity rises to a significant value. An alternate and accurate method which would be applicable to liquids of significant viscosity is therefore needed. We plan to develop such a capability by measuring the equilibrium shape of levitated undercooled melt droplets as they undergo solid-body rotation. The experimental measurement of the characteristic point of transition (bifurcation point) between axisymmetric and two-lobed shapes will be used to calculate the surface tension of the liquid. Such an approach has already been validated through the experimental verification of numerical modeling results. The experimental approach involves levitation, melting, and solidification of undercooled droplets using a hybrid ultrasonic-electrostatic technique in both a gaseous as well as a vacuum environment. A shape relaxation method will be investigated in order to derive a reliable method to measure the viscosity of undercooled melts. The analysis of the monotonic relaxation to equilibrium shape of a drastically deformed and super-critically damped free drop has been used to derive interfacial tension of immiscible liquid combinations where one of the component has high viscosity. A standard approach uses the initial elongation of a droplet through shear flows, but an equivalent method could involve the initial deformation of a drop levitated in a gas by ultrasonic radiation pressure, electric stresses, or even solid body rotation. The dynamic behavior of the free drop relaxing back to equilibrium shape will be modeled, and its characteristic time dependence should provide a quantitative means to evaluate the liquid viscosity.

Trinh, E. H.↗
NASA NTRS2005

Dynamic Weakening (Extinction) of Simple Hydrocarbon-air Counterflow Diffusion Flames by Oscillatory Inflows

This study of laminar non-premixed HC-air flames used an Oscillatory-input Opposed Jet Burner (OOJB) system developed from a previously well-characterized 7.2-mm Pyrex-nozzle OJB system. Over 600 dynamic Flame Strength (FS) measurements were obtained on unanchored (free-floating) laminar Counterflow Diffusion Flames (CFDF's). Flames were stabilized using plug inflows having steady-plus-sinusoidal axial velocities of varied magnitude, frequency, f, up to 1600 Hz, and phase angle from 0 (most data) to 360 degrees. Dynamic FS is defined as the maximum average air input velocity (U(sub air), at nozzle exit) a CFDF can sustain before strain-induced extinction occurs due to prescribed oscillatory peak-to-peak velocity inputs superimposed on steady inputs. Initially, dynamic flame extinction data were obtained at low f, and were supported by 25-120 Hz Hot-Wire cold-flow velocity data at nozzle exits. Later, expanded extinction data were supported by 4-1600 Hz Probe Microphone (PM) pk/pk P data at nozzle exits. The PM data were first obtained without flows, and later with cold stagnating flows, which better represent speaker-diaphragm dynamics during runs. The PM approach enabled characterizations of Dynamic Flame Weakening (DFW) of CFDF's from 8 to 1600 Hz. DFW was defined as % decrease in FS per Pascal of pk/pk P oscillation, namely, DFW = - 100 d(U(sub air) / U(sub air),0Hz) / d(pkpk P). The linear normalization with respect to acoustic pressure magnitude (and steady state (SS) FS) led to a DFW unaffected by strong internal resonances. For the C2H4/N2-air system, from 8 to 20 Hz, DFW is constant at 8.52 plus or minus 0.20 (% weakening)/Pa. This reflects a quasi-steady flame response to an acoustically induced dU(sub air)/dP. Also, it is surprisingly independent of C2H4/N2 mole fraction due to normalization by SS FS. From 20 to approximately 150 Hz, the C2H4/N2 air-flames weakened progressively less, with an inflection at approximately 70 Hz, and became asymptotically insensitive (DFW approximately 0) at approximately 300 Hz, which continued to 1600 Hz. The DFW of CH4-air flames followed a similar pattern, but showed much greater weakening than C2H4/N2-air flames; i.e., the quasi-steady DFW (8 to approximately 15 Hz) was 44.3 %/Pa, or approximately 5x larger, even though the 0 Hz (SS) FS was only 3.0 x smaller. The quasi-steady DFW's of C3H8-air and C2H6-air were intermediate at 34.8 and 20.9 %Pa, respectively. The DFW profiles of all four fuels, at various frequencies, correlated well but non-linearly with respective SS FS's. Notably, the DFW profile for C3H8 air fell more rapidly in the range greater than 15 to 60 Hz, compared with the 1- and 2-carbon fuels. This may indicate a shift in chemical kinetics, and/or O2 transport to a flame that moved closer to the fuel-side. In conclusion, Dynamic Flame Weakening limits appear significant and unique for each fuel, and correlate closely, but non-linearly, with Steady-State Flame Strengths at any given frequency. For reasons unknown, the dynamic flames didn't weaken more at intermediate frequencies (e.g., at 20-50 Hz) than they did at low frequencies (less than 15 Hz), where quasi-steady weakening appears to dominate. Quasi-steady flame weakening ostensibly represents a transient input strain rate maximum that just exceeds the steady-state strain-rate-limited extinction limit for a few cycles. Clearly, further detailed mechanistic understanding is needed in the fall-off region.

Pellett, G.↗
DOE OSTI2020

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗
NASA NTRSDate not supplied

Aero-Propulsive Damping Characterization for eVTOL Aircraft Using Free Motion Wind-Tunnel Testing

This paper describes an electric vertical takeoff and landing (eVTOL) aircraft system identification method applied using three degree-of-freedom (3DOF) free motion wind-tunnel testing. The approach, similar to flight-test system identification, allows for efficient mathematical model development of the aero-propulsive moments applied on an eVTOL vehicle, including aerodynamic damping effects. The approach is demonstrated using a subscale tiltrotor eVTOL aircraft mounted on a new 3DOF wind-tunnel apparatus. To execute the test, a model-based 3DOF control system is designed to track attitude commands and transition the aircraft based on the freestream dynamic pressure. While the flight controller is active, orthogonal phase-optimized multisine inputs are injected into the attitude command and control effector command signals to enable collection of informative data for model identification. Aero-propulsive models are then identified at several reference conditions in the transition flight envelope using the equation-error method in the frequency domain. The identified models are shown to have a good fit to the modeling data and good prediction capability of data not used for model identification. The method yields aerodynamic damping estimates using less wind-tunnel test time compared to traditional forced oscillation experiments and supplements static wind-tunnel testing to produce a comprehensive transition aero-propulsive model suitable for use in flight dynamics simulations.

RAVEN↗
NASA NTRS2012

Optical Phase Recovery and Locking in a PPM Laser Communication Link

Free-space optical communication holds great promise for future space missions requiring high data rates. For data communication in deep space, the current architecture employs pulse position modulation (PPM). In this scheme, the light is transmitted and detected as pulses within an array of time slots. While the PPM method is efficient for data transmission, the phase of the laser light is not utilized. The phase coherence of a PPM optical signal has been investigated with the goal of developing a new laser communication and ranging scheme that utilizes optical coherence within the established PPM architecture and photon-counting detection (PCD). Experimental measurements of a PPM modulated optical signal were conducted, and modeling code was developed to generate random PPM signals and simulate spectra via FFT (Fast Fourier Transform) analysis. The experimental results show very good agreement with the simulations and confirm that coherence is preserved despite modulation with high extinction ratios and very low duty cycles. A real-time technique has been developed to recover the phase information through the mixing of a PPM signal with a frequency-shifted local oscillator (LO). This mixed signal is amplified, filtered, and integrated to generate a voltage proportional to the phase of the modulated signal. By choosing an appropriate time constant for integration, one can maintain a phase lock despite long dark times between consecutive pulses with low duty cycle. A proof-of-principle demonstration was first achieved with an RF-based PPM signal and test setup. With the same principle method, an optical carrier within a PPM modulated laser beam could also be tracked and recovered. A reference laser was phase-locked to an independent pulsed laser signal with low-duty-cycle pseudo-random PPM codes. In this way, the drifting carrier frequency in the primary laser source is tracked via its phase change in the mixed beat note, while the corresponding voltage feedback maintains the phase lock between the two laser sources. The novelty and key significance of this work is that the carrier phase information can be harnessed within an optical communication link based on PPM-PCD architecture. This technology development could lead to quantum-limited efficient performance within the communication link itself, as well as enable high-resolution optical tracking capabilities for planetary science and spacecraft navigation.

Aveline, David C.↗
DOE OSTI2022

Radiation-Hardened GaN HEMT and Cell Design, Modeling, and Fabrication for Nuclear Instrumentation Applications

Recent advances in nuclear power generation technologies show great promise for efficient, safe generation of carbon-free power in the near future. Consequently, many new generation reactor designs are being pursued by government and industry groups, including both fission- and fusion-based technologies. Radiation- and temperature-tolerant sensor technology continues to advance, yet development of the electronics/instrumentation technologies suitable for these environments have fallen behind. Industry has adequately addressed harsh environment electronics needs for low earth orbit satellites, but the much more extreme conditions associated with electronics placed near or in a reactor core remain unaddressed. This project investigates the use of gallium nitride (GaN) circuit technology to address the unique needs for sensor interface electronics and communications in reactor environments. This report provides a summary of the first year’s project activities related to developing and optimizing devices and circuits for this GaN high electronic mobility transistor (HEMT) process. Activities reported include multiple analog and digital circuit designs simulated using custom Verilog-A models generated from measured GaN devices fabricated in the target process. Digital circuit designs included fundamental logic circuits such as an inverter, NAND, NOR, AND and OR gates, as well as a 5-stage ring oscillator. Each design was simulated using an open-source SPICE simulator, and an integrated circuit layout was produced of each design for use in future fabrication. Analog circuit designs focused on the fundamental building blocks to be used to construct sensor interface and communications circuits. These included current mirrors, matched differential pairs, and a Gilbert cell mixer design, each with an associated integrated circuit layout. Layout tools used for these designs included two open-source packages—KLayout and Magic—which were customized for specific use with the OSU GaN process layers. Design rules and preliminary device extraction capabilities were built for the Magic tool to enable device extraction and subsequent netlist generation for SPICE simulation. Finally, future research directions for year two are summarized. When the year-two directions are implemented, the project will be in position to advance the state of the art in electronics for near- or in-reactor sensor interfacing and communications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1998

An Assessment of Molten Metal Detachment Hazards During Electron Beam Welding in Space

The safety issue has been raised with regards to potential molten metal detachments from the weld pool and cold filler wire during electron beam welding in space. This investigation was undertaken to evaluate if molten metal could detach and come in contact with astronauts and burn through the fabric of the astronauts' Extravehicular Mobility Unit (EMU) during electron beam welding in space. Molten metal detachments from either the weld/cut substrate or weld wire could present harm to a astronaut if the detachment was to burn through the fabric of the EMU. Theoretical models were developed to predict the possibility and size of the molten metal detachment hazards during the electron beam welding exercises at Low Earth Orbit (LEO). The primary molten metal detachment concerns were those cases of molten metal separation from the metal surface due to metal cutting, weld pool splashing, entrainment and release of molten metal due to filler wire snap-out from the weld puddle, and molten metal accumulation and release from the end of the weld wire. Some possible ways of obtaining molten metal drop detachments would include an impulse force, or bump, to the weld sample, cut surface, or filler wire. Theoretical models were developed for these detachment concerns from principles of impact and kinetic energies, surface tension, drop geometry, surface energies, and particle dynamics. The surface tension represents the force opposing the liquid metal drop from detaching whereas the weight of the liquid metal droplet represents a force that is tending to detach the molten metal drop. Theoretical calculations have indicated that only a small amount of energy is required to detach a liquid metal drop; however, much of the energy of an impact is absorbed in the sample or weld plate before it reaches the metal drop on the cut edge or surface. The tendency for detachment is directly proportional to the weld pool radius and metal density and inversely proportional to the surface tension of the liquid metal. For a detachment the initial kinetic energy of the weld pool with respect to the plate has to exceed the energy to form the extra surface required for the detachment of the pool. The difficulty is in transferring the energy from the point of impact through the plate and sample to the cut edge. It is likely that not all of the kinetic energy is available for detaching the pool; some may be sequestered in weld pool oscillations. The coefficient of restitution for the collision will be lower than one if irreversible deformation, for example plastic flow deformation, takes place during the collision. Thus determining the amount of energy from an impact that actually reaches the molten metal droplet is critical. Various molten metal detachment scenarios were tested experimentally in an enclosed vacuum chamber using the Ukrainian Universal Hand Tool, an electron beam welder designed for space welding. The experimental testing was performed in a 4 ft. X 4 ft. vacuum chamber at Marshall Space Flight Center, evacuated to vacuum levels of at least 50 microTorr, and also some welding garment material was utilized to observe the effect of the molten metal detachments on the material. A "carillon" apparatus consisting of four pendulum hammer strikers, each weighing approximately 3.65 lbs, raised to predetermined specific heights was used to apply an impact force to the weld sample/plate during electron beam welding and cutting exercises. The strikers were released by switching on an electric motor to rotate a pin holding wires retaining the strikers at desired heights. The specimens were suspended so as to be free to respond to the blows with a sudden velocity increment. The specimens were mounted on a hinged plate for minimizing effective mass with the option to fasten it down so as to raise its effective mass closer to that anticipated for an actual space welding scenario. Measurements were made of the impact energy and the horizontal fling distances of the detached metal drops. It was not particularly easy to generate the detachments for this experiment. This document presents the details of the theoretical modeling effort and a summary of the experimental effort to measure molten metal drop detachments from terrestrial electron beam welding in the enclosed vacuum chamber. The results of the experimental effort have shown that molten metal detachments can occur from the sample/weld plate only if a sufficiently large impact force is applied to the weld plate. A "weld pool detachment parameter" was determined to indicate whether detachment would occur. Detachment can be either full or partial (dripping), Partial detachment means that the weld pool detached from one side of the liquid-solid boundary so as to leave a hole at the puddle site but remained attached over part of the liquid-solid boundary and dripped down the plate with no fully detached material detected. Full detachment, however, does not necessarily mean that the whole pool fully detached; in some cases only a smaller portion of the pool detached, the remainder dripping down the plate. The weld pool detachment parameter according to theory and according to the empirical data allows a determination of whether full detachments might occur. Theoretical calculations indicated titanium alloy would be the most difficult from which to detach molten metal droplets followed by stainless steel and then by aluminum. The experimental results were for the most part consistent with the theoretical analysis and predictions. The above theory is applicable to other situations as desired for assessing the potential for molten metal detachments.

Fragomeni, James M.↗
DOE OSTI2020

Extended Lagrangian Born–Oppenheimer molecular dynamics using a Krylov subspace approximation

It is shown how the electronic equations of motion in extended Lagrangian Born–Oppenheimer molecular dynamics simulations can be integrated using low-rank approximations of the inverse Jacobian kernel. This kernel determines the metric tensor in the harmonic oscillator extension of the Lagrangian that drives the evolution of the electronic degrees of freedom. The proposed kernel approximation is derived from a pseudoinverse of a low-rank estimate of the Jacobian, which is expressed in terms of a generalized set of directional derivatives with directions that are given from a Krylov subspace approximation. The approach allows a tunable and adaptive approximation that can take advantage of efficient preconditioning techniques. The proposed kernel approximation for the integration of the electronic equations of motion makes it possible to apply extended Lagrangian first-principles molecular dynamics simulations to a broader range of problems, including reactive chemical systems with numerically sensitive and unsteady charge solutions. This can be achieved without requiring exact full calculations of the inverse Jacobian kernel in each time step or relying on iterative non-linear self-consistent field optimization of the electronic ground state prior to the force evaluations as in regular direct Born–Oppenheimer molecular dynamics. We note the low-rank approximation of the Jacobian is directly related to Broyden’s class of quasi-Newton algorithms and Jacobian-free Newton–Krylov methods and provides a complementary formulation for the solution of nonlinear systems of equations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

NUMERICAL CALCULATION OF LOSSES OF TRAPPED VORTICES UNDER STRONG RF MEISSNER CURRENT AND DC SUPERHEATING FIELD IN TYPE II SUPERCONDUCTORS

Research on the vortex dynamics and enhancing of superheating field in superconductors has attracted much attention in accelerator physics community to develop next-generation high-performance accelerator cavities. However, the extreme dynamics of curvilinear elastic vortices driven by very strong currents close to the depairing limit or superheating field of a superconductor with a nanostructured surface has not been well understood. We calculated the superheating field Hsh and critical momentum kc characterizing the wavelength of the instability ?m of the Meissner state to flux penetration by solving numerically the Ginzburg?Landau equations. A bulk superconductor, superconductor with the inhomogeneous surface disorder (S-S), and multilayered surface (S-I-S) have been thoroughly investigated in this work. Our result showed that S-S and S-I-S structures can enhance the superheating field well above their clean limit. In this work extensive numerical simulation of the power dissipated by an oscillating vortex segment driven by the surface ac Meissner currents was performed. Our simulations take into account the nonlinear vortex line tension, vortex mass, Bardeen?Stephen viscous vortex drag applicable at low fields, and nonlinear Larkin-Ovchinnikov (LO) viscous drag coefficient ?(v) at high fields and pinning force. We showed that the LO decrease of ?(v) with the vortex velocity v could radically change the field dependence of the surface resistance Ri(H) caused by trapped vortices. At low frequencies Ri(H) exhibits a conventional increase with H. However, as frequency increases, the surface resistance becomes a nonmonotonic function of H which decreases with H at higher fields irrespective of the pinning distribution. Overheating can mask the descending field dependence of Ri(H) as frequency increases. Our numerical simulations also show that the LO effect can cause a vortex bending instability at high field amplitudes and frequencies, giving rise to the formation of dynamic kinks along with the vortex when a vortex is pinned strongly to one end. Nonlinear losses of trapped vortices in thick films under high-amplitude RF fields as functions of frequency, mean free path and pinning characteristics have been calculated.

Pathirana, Walive↗
NASA NTRS1986

Analytical and numerical studies of the thermocapillary flow in a uniformly floating zone

The microgravity environment of an orbiting vehicle permits crystal growth experiments in the presence of greatly reduced buoyant convection in the liquid melt. Crystals grown in ground-based laboratories do not achieve their potential properties because of dopant variations caused by flow in the melt. The floating zone crystal growing system is widely used to produce crystals of silicon and other materials. However, in this system the temperature gradient on the free sidewall surface of the melt is the source of a thermocapillary flow which does not disappear in the low-gravity environment. The idea of using a uniform rotation of the floating zone system to confine the thermocapillary flow to the melt sidewall leaving the interior of the melt passive is examined. A cylinder of fluid with an axial temperature gradient imposed on the cylindrical sidewall is considered. A half zone and the linearized, axisymmetric flow in the absence of crystal growth is examined. Rotation is found to confine the linear thermocapillary flow. A simplified model is extended to a full zone and both linear and nonlinear thermocapillary flows are studied theoretically. Analytical and numerical methods are used for the linear flows and numerical methods for the nonlinear flows. It was found that the linear flows in the full zone have more complicated and thicker boundary layer structures than in the half zone, and that these flows are also confined by the rotation. However, for the simplified model considered and for realistic values for silicon, the thermocapillary flow is not linear. The fully nonlinear flow is strong and unsteady (a weak oscillation is present) and it penetrates the interior. Some non-rotating flow results are also presented. Since silicon as a large value of thermal conductivity, one would expect the temperature fields to be determined by conduction alone. This is true for the linear and weakly nonlinear flows, but for the stronger nonlinear flow the results show that temperature advection is also important. Uniform rotation may still be a means of confining the flow and the results obtained define the procedure to be used to examine this hypothesis.

Fowlis, W. W.↗
DOE OSTI2025

DESI 2024 VI: cosmological constraints from the measurements of baryon acoustic oscillations

We present cosmological results from the measurement of baryon acoustic oscillations (BAO) in galaxy, quasar and Lyman-α forest tracers from the first year of observations from the Dark Energy Spectroscopic Instrument (DESI), to be released in the DESI Data Release 1. DESI BAO provide robust measurements of the transverse comoving distance and Hubble rate, or their combination, relative to the sound horizon, in seven redshift bins from over 6 million extragalactic objects in the redshift range 0.1 < z < 4.2. To mitigate confirmation bias, a blind analysis was implemented to measure the BAO scales. DESI BAO data alone are consistent with the standard flat ΛCDM cosmological model with a matter density Ω m =0.295±0.015. Paired with a baryon density prior from Big Bang Nucleosynthesis and the robustly measured acoustic angular scale from the cosmic microwave background (CMB), DESI requires H 0 =(68.52±0.62) km s -1 Mpc -1 . In conjunction with CMB anisotropies from Planck and CMB lensing data from Planck and ACT, we find Ω m =0.307± 0.005 and H 0 =(67.97±0.38) km s -1 Mpc -1 . Extending the baseline model with a constant dark energy equation of state parameter w, DESI BAO alone requirew=-0.99 +0.15 -0.13 . In models with a time-varying dark energy equation of state parametrised by w 0 and w a , combinations of DESI with CMB or with type Ia supernovae (SN Ia) individually prefer w 0 > -1 and w a < 0. This preference is 2.6σ for the DESI+CMB combination, and persists or grows when SN Ia are added in, giving results discrepant with the ΛCDM model at the 2.5σ, 3.5σ or 3.9σ levels for the addition of the Pantheon+, Union3, or DES-SN5YR supernova datasets respectively. For the flat ΛCDM model with the sum of neutrino mass ∑ m ν free, combining the DESI and CMB data yields an upper limit ∑ m ν < 0.072 (0.113) eV at 95% confidence for a ∑ m ν > 0 (∑ m ν > 0.059) eV prior. These neutrino-mass constraints are substantially relaxed if the background dynamics are allowed to deviate from flat ΛCDM.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Nonlinear dynamics of drops and bubbles and chaotic phenomena

Nonlinear phenomena associated with the dynamics of free drops and bubbles are investigated analytically, numerically and experimentally. Although newly developed levitation and measurement techniques have been implemented, the full experimental validation of theoretical predictions has been hindered by interfering artifacts associated with levitation in the Earth gravitational field. The low gravity environment of orbital space flight has been shown to provide a more quiescent environment which can be utilized to better match the idealized theoretical conditions. The research effort described in this paper is a closely coupled collaboration between predictive and guiding theoretical activities and a unique experimental program involving the ultrasonic and electrostatic levitation of single droplets and bubbles. The goal is to develop and to validate methods based on nonlinear dynamics for the understanding of the large amplitude oscillatory response of single drops and bubbles to both isotropic and asymmetric pressure stimuli. The first specific area on interest has been the resonant coupling between volume and shape oscillatory modes isolated gas or vapor bubbles in a liquid host. The result of multiple time-scale asymptotic treatment, combined with domain perturbation and bifurcation methods, has been the prediction of resonant and near-resonant coupling between volume and shape modes leading to stable as well as chaotic oscillations. Experimental investigations of the large amplitude shape oscillation modes of centimeter-size single bubbles trapped in water at 1 G and under reduced hydrostatic pressure, have suggested the possibility of a low gravity experiment to study the direct coupling between these low frequency shape modes and the volume pulsation, sound-radiating mode. The second subject of interest has involved numerical modeling, using the boundary integral method, of the large amplitude shape oscillations of charged and uncharged drops in the presence of a static or time-varying electric field. Theoretically predicted non linearity in the resonant frequency of the fundamental quadrupole mode has been verified by the accompanying experimental studies. Additional phenomena such as hysteresis in the frequency response of ultrasoncially levitated droplets in the presence of a time varying electric field, and mode coupling in the oscillations of ultrasonically modulated droplets, have also been uncovered. One of the results of this ground-based research has been the identification and characterization of phenomena strictly associated with the influence of the gravitational field. This has also allowed us to identify the specific requirements for potential microgravity investigations yielding new information not obtainable on Earth.

Trinh, Eugene H.↗
DOE OSTI2019

Microscopic subsurface characterization of layered magnetic materials using magnetic resonance force microscopy

This project resulted in the following accomplishments: 1) The invention of scanned probe ferromagnetic resonance microscopy. Upon implementation and refinement this approach was used to measure static and dynamic magnetic resonance properties of ferromagnets with an unprecedented combination of spectral and spatial resolution. 2) Demonstration and study of spin pumping in a pristine ferromagnet across a boundary defined entirely by a magnetic field gradient, and so free of materials discontinuities. 3) First demonstrated highly efficient dynamic spin injection from Y 3 Fe 5 O 12 (YIG) into NiO, an antiferromagnetic (AF) insulator and robust spin propagation in AF NiO over distances up to 100-nm thickness; both spin injection and spin transport are mediated by the AF spin correlations. 4) Demonstration of exceptionally high efficiency conversion of spin current to charge current achieved in spin pumping from the dynamic magnetization of a ferromagnetic insulator into a variety of adjacent metal films. 5) Seminal study of the nature of spin-orbit coupling in d-electron metals that revealed the dominant sensitivity of the d-orbital moment to the d-electron count and revealing the important roles of both d-electron configuration and atomic number in the physics of the spin Hall effect. 6) Demonstration of observation of FMR by means of optically-detected NV diamond relaxometry that provided a highly sensitive, spatially resolved FMR detection method and revealed spin pumping across 100-nm length scales by means of dipolar coupling. 7) Using magnetic field localization, developed for FMR imaging, demonstrated a novel approach to discretizing spin wave modes to improve the characteristics of auto-oscillators driven by pure spin currents produced by the spin Hall effect. 8) We demonstrated large spin pumping signals from propagating high wavevector spinwaves in thin film YIG find that the spin pumping efficiency was insensitive to wavevector thus emphasizing the promise of thin YIG films for spinwave-based spintronics by allowing for predictable down-scaling of device dimensions and enabling a broad range of operational frequencies.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2003

Linear Back-Drive Differentials

Linear back-drive differentials have been proposed as alternatives to conventional gear differentials for applications in which there is only limited rotational motion (e.g., oscillation). The finite nature of the rotation makes it possible to optimize a linear back-drive differential in ways that would not be possible for gear differentials or other differentials that are required to be capable of unlimited rotation. As a result, relative to gear differentials, linear back-drive differentials could be more compact and less massive, could contain fewer complex parts, and could be less sensitive to variations in the viscosities of lubricants. Linear back-drive differentials would operate according to established principles of power ball screws and linear-motion drives, but would utilize these principles in an innovative way. One major characteristic of such mechanisms that would be exploited in linear back-drive differentials is the possibility of designing them to drive or back-drive with similar efficiency and energy input: in other words, such a mechanism can be designed so that a rotating screw can drive a nut linearly or the linear motion of the nut can cause the screw to rotate. A linear back-drive differential (see figure) would include two collinear shafts connected to two parts that are intended to engage in limited opposing rotations. The linear back-drive differential would also include a nut that would be free to translate along its axis but not to rotate. The inner surface of the nut would be right-hand threaded at one end and left-hand threaded at the opposite end to engage corresponding right- and left-handed threads on the shafts. A rotation and torque introduced into the system via one shaft would drive the nut in linear motion. The nut, in turn, would back-drive the other shaft, creating a reaction torque. Balls would reduce friction, making it possible for the shaft/nut coupling on each side to operate with 90 percent efficiency.

Waydo, Peter↗
NASA NTRS2020

The Role of Oceanic Processes in the Initiation of Boreal Winter Intraseasonal Oscillations Over the Indian Ocean

Observational analyses and ocean general circulation model experiments were performed to understand the influence of oceanic processes on warm sea surface temperature anomalies (SSTAs) prior to convection initiation of boreal winter intraseasonal oscillations (ISOs), including the Madden‐Julian Oscillation (MJO), in the tropical Indian Ocean. We found 39 strong ISOs that passed over the Indian Ocean Warm Pool (WP) region during the November–April season of the 2001–2012 period. Seventeen of them (44%) initiated over the Seychelles‐Chagos Thermocline Ridge (SCTR) before propagating eastward to the WP and Maritime Continent. Including only global‐scale MJOs, 71% (24%) of the remaining MJO events initiated from the SCTR (WP). Four (seven) SCTR (WP) ISO events were preceded by SSTAs that were strongly influenced by wind stress‐driven oceanic processes; all four (2/7) SCTR (WP) events were MJOs. Composite analyses of the four oceanic process influenced SCTR MJO events showed that wind stress‐driven oceanic processes contributed about 30% (between 0.1 and 0.2 °C) of the preconvection warming in the SCTR, which coincides with the Intertropical Convergence Zone with mean SST > 28 °C. Reduced upwelling and entrainment played a more important role than horizontal advection in elevating the +SSTAs. Case studies revealed that the two ocean process‐influenced primary MJO SCTR events were also associated with oceanic equatorial Rossby waves. Furthermore, results of a linear atmospheric mixed layer model indicated that boundary layer processes related to SSTAs contributed approximately half of the total convergence during MJO initiation, relative to the effect of the free troposphere, for the Rossby wave‐influenced events.

MJO Initiation↗
NASA NTRSDate not supplied

Global Precipitation Means and Variations with the New Version of GPCP

Knowledge of global precipitation means, patterns and variations is essential for understanding the global water cycle. The observation-based global analysis of the Global Precipitation Climatology Project (GPCP) has been a key input to many such studies, including those related to means of the global (and regional) water and energy cycles. A new version (Version 3.1) of the GPCP monthly analysis is now available (1983-2019), with finer spatial resolution (0.5º latitude/longitude), updated satellite algorithms, latest gauge analysis over land from the Global Precipitation Climatology Center (GPCC), and with ocean climatologies adjusted using information from the Tropical Rainfall Measuring Mission (TRMM), the Global Precipitation Measurement (GPM) mission and CloudSat. The presentation will give key findings from the new analysis, compare with the previous version, link to studies of the water cycle and compare to CMIP6 climate model results. The GPCP Monthly V3.1 analysis uses a Tropical Composite Climatology (TCC) using 22 years (1998-2019) of TRMM and GPM-based surface precipitation estimates from passive microwave, radar and combined passive microwave and radar observations to adjust the long-term mean values in the tropics over ocean. At higher latitudes over ocean climatological values of merged CloudSat and GPM combined passive microwave and radar surface precipitation estimates are used. The result of the new algorithms and procedures is an ocean mean value 60º N to 60º S of 3.2 mm/d, an increase of 7% from the older V2.3. Over land the gauge analysis (from GPCC in Germany), combined with satellite estimates results in a small decrease (~ 0.5%) from the previous version, giving a total global precipitation mean of 2.81 mm/d for V3.1, an increase of 4.5%. Variations from inter-annual to trend scales over most of the ocean are driven by the Colorado State University (CSU) Goddard Profiling (GPROF) algorithm applied to SSMI/SSMIS satellite data. The new GPCP version retains a near zero trend of global precipitation, with significant positive trends in the deep tropics along the Pacific ITCZ and elsewhere, countered by middle latitude decreases, very similar to the previous version of the GPCP Monthly analysis. The pattern of trends is similar to that of AMIP climate model results, driven by observed SSTs for the period in question, but very different from results for “free-running” CMIP historical ensembles, likely due in part to the relatively short comparison period and effects of inter-decadal variations in the GPCP and AMIP results that are not in the “history” models. Interannual variations are also nearly the same in the new version, but with finer detail, although ENSO variations over the ocean have slightly larger amplitude than before, an effect likely related to the change in ocean satellite algorithm used in Version 3. For precipitation intensity (percentiles) at the monthly scale, especially in the tropics, GPCP shows a positive trend for the upper one third of the percentiles (Pct ≥ 70th) and a much weaker positive trend for the lowest percentiles (Pct ≤ 10th), while negative trends appear for the middle one-half percentiles (20th-65th). AMIP results agree with those from the GPCP in terms of the sign of the changes/trends for high and intermediate percentiles. The CMIP historical results also agree in the sign of the trends, but the trends are weaker. Comparisons with other types of CMIP historical simulations including the GHG-only, aerosol-only, and nature-only runs suggest that the observed changes/trends in precipitation amount and intensity during the GPCP period are dominated by a combination of the effects of the Pacific Decadal Oscillation (PDO) and anthropogenic GHG-related surface warming.

Robert Adler↗
NASA NTRS1999

Theory of Dust Voids in Plasmas

Dusty plasmas in a gas discharge often feature a stable void, i.e., a dust-free region inside the dust cloud. This occurs under conditions relevant to both plasma processing discharges and plasma crystal experiments. The void results from a balance of the electrostatic and ion drag forces on a dust particle. The ion drag force is driven by a flow of ions outward from an ionization source and toward the surrounding dust cloud, which has a negative space charge. In equilibrium the force balance for dust particles requires that the boundary with the dust cloud be sharp, provided that the particles are cold and monodispersive. Numerical solutions of the one-dimensional nonlinear fluid equations are carried out including dust charging and dust-neutral collisions, but not ion-neutral collisions. The regions of parameter space that allow stable void equilibria are identified. There is a minimum ionization rate that can sustain a void. Spatial profiles of plasma parameters in the void are reported. In the absence of ion-neutral collisions, the ion flow enters the dust cloud's edge at Mach number M = 1. Phase diagrams for expanding or contracting voids reveal a stationary point corresponding to a single stable equilibrium void size, provided the ionization rate is constant. Large voids contract and small voids expand until they attain this stationary void size. On the other hand, if the ionization rate is not constant, the void size can oscillate. Results are compared to recent laboratory and microgravity experiments.

Goree, J.↗