Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Adiabatic approximation”

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 415 records · Page 23

Solar modulation of galactic antiprotons

Galactic antiproton data of current interest lie in an energy regime heavily influenced by solar modulation. Correcting for it needs to be done more carefully than it has been in the past. The well-known force-field analytic approximation of the spherically-symmetric, steady-state, cosmic-ray transport equation is applied in order to account for modulation down to at least 100 MeV. A sample solution which applies to the currently available antiproton data set (1979-80), and can be used to accurately modulate any possible interstellar antiproton spectrum, is given. The solution is easily adapted for comparison to future measurements. It also shows that boosting the low-energy (less than 600 MeV) side of the interstellar antiproton spectrum will not affect the low-energy spectrum at 1 AU, due to strong adiabatic deceleration during that time.

Perko, J. S.↗

Large-angle cosmic microwave background anisotropies in an open universe

If the universe is open, scales larger than the curvature scale may be probed by observation of large-angle fluctuations in the cosmic microwave background (CMB). We consider primordial adiabatic perturbations and discuss power spectra that are power laws in volume, wavelength, and eigenvalue of the Laplace operator. Such spectra may have arisen if, for example, the universe underwent a period of `frustated' inflation. The resulting large-angle anisotropies of the CMB are computed. The amplitude generally increases as Omega is decreased but decreases as h is increased. Interestingly enough, for all three Ansaetze, anisotropies on angular scales larger than the curvature scale are suppressed relative to the anisotropies on scales smaller than the curvature scale, but cosmic variance makes discrimination between various models difficult. Models with 0.2 approximately less than Omega h approximately less than 0.3 appear compatible with CMB fluctuations detected by Cosmic Background Explorer Satellite (COBE) and the Tenerife experiment and with the amplitude and spectrum of fluctuations of galaxy counts in the APM, CfA, and 1.2 Jy IRAS surveys. COBE normalization for these models yields sigma(sub 8) approximately = 0.5 - 0.7. Models with smaller values of Omega h when normalized to COBE require bias factors in excess of 2 to be compatible with the observed galaxy counts on the 8/h Mpc scale. Requiring that the age of the universe exceed 10 Gyr implies that Omega approximately greater than 0.25, while requiring that from the last-scattering term in the Sachs-Wolfe formula, large-angle anisotropies come primarily from the decay of potential fluctuations at z approximately less than 1/Omega. Thus, if the universe is open, COBE has been detecting temperature fluctuations produced at moderate redshift rather than at z approximately 1300.

Kamionkowski, Marc↗

Further observational support for the lossy radial diffusion model of the inner Jovian magnetosphere

A mathematical model describing radial diffusion, with violation of the third adiabatic invariant and with local losses, is applied to the Pioneer 10 and 11 observations of omnidirectional integral intensities of electrons with energy exceeding 21 MeV. Local losses outside L = 3 are much stronger than synchrotron losses but weaker than the losses one would expect for strong pitch angle diffusion. If radial diffusion is driven by ionospheric winds, as was suggested by Brice and McDonough (1973), it is found that the theoretical solution which best fits the data requires that the radial diffusion coefficient equals 3 x 10 to the -9th L-cubed/s and the lifetime against local losses is approximately a few million seconds.

Goertz, C. K.↗

Radio science with Voyager 1 at Jupiter - Preliminary profiles of the atmosphere and ionosphere

A preliminary profile of the atmosphere of Jupiter in the South Equatorial Belt shows (1) the tropopause occurring at a pressure level of 100 millibars and temperature of about 113 K, (2) a higher warm inversion layer at about the 35-millibar level, and (3) a lower-altitude constant lapse rate matching the adiabatic value of about 2 K/km, with the temperature reaching 150 K at the 600-millibar level. Preliminary afternoon and predawn ionospheric profiles at 12 deg and near the equator, respectively, have topside plasma scale heights of 590 km changing to 960 km above an altitude of 3500 km for the dayside, and about 960 km at all measured heights above the peak for the nightside. The higher value of scale height corresponds to a plasma temperature of 1100 K under the assumption of a plasma of protons and electrons in ambipolar diffusive equilibrium. The peak electron concentration in the upper ionosphere is approximately 200,000/cu cm for the dayside and about a factor of 10 less for the nightside. These peaks occur at altitudes of 1600 and 2300 km, respectively. Continuing analyses are expected to extend and refine these results, and to be used to investigate other regions and phenomena.

Eshleman, V. R.↗

Cirrus Parcel Model Comparison Project

The Cirrus Parcel Model Comparison (CPMC) is a project of the GEWEX Cloud System Study Working Group on Cirrus Cloud Systems (GCSS WG2). The primary goal of this project is to identify cirrus model sensitivities to the state of our knowledge of nucleation and microphysics. Furthermore, the common ground of the findings may provide guidelines for models with simpler cirrus microphysics modules. We focus on the nucleation regimes of the warm (parcel starting at -40 C and 340 hPa) and cold (-60 C and 170 hPa) cases studied in the GCSS WG2 Idealized Cirrus Model Comparison Project. Nucleation and ice crystal growth were forced through an externally imposed rate of lift and consequent adiabatic cooling. The background haze particles are assumed to be lognormally-distributed H2SO4 particles. Only the homogeneous nucleation mode is allowed to form ice crystals in the HN-ONLY runs; all nucleation modes are switched on in the ALL-MODE runs. Participants were asked to run the HN-lambda-fixed runs by setting lambda = 2 (lambda is further discussed in section 2) or tailoring the nucleation rate calculation in agreement with lambda = 2 (exp 1). The depth of parcel lift (800 m) was set to assure that parcels underwent complete transition through the nucleation regime to a stage of approximate equilibrium between ice mass growth and vapor supplied by the specified updrafts.

Lin, Ruei-Fong↗

Charged particle motions in the distended magnetospheres of Jupiter and Saturn

Charged particle motion in the guiding center approximation is analyzed for models of the Jovian and Saturnian magnetospheric magnetic fields based on Voyager magnetometer observations. Field lines are traced and exhibit the distention which arises from azimuthally circulating magnetospheric currents. The spatial dependencies of the guiding center bounce period and azimuthal drift rate are investigated for the model fields. Non-dipolar effects in the gradient-curvature drift rate are most important at the equator and affect particles with all mirror latitudes. The effect is a factor of 10-15 for Jupiter with its strong magnetodisc current and 1-2 for Saturn with its more moderate ring current. Limits of adiabaticity, where particle gyroradii become comparable with magnetic scale lengths, are discussed and are shown to occur at quite modest kinetic energies for protons and heavier ions.

Birmingham, T. J.↗

Nonguiding Center Motion and Substorm Effects in the Magnetotail

Thick and thin models of the middle magnetotail were developed using a consistent orbit tracing technique. It was found that currents carried near the equator by groups of ions with anisotropic distribution functions are not well approximated by the guiding center expressions. The guiding center equations fail primarily because the calculated pressure tensor is not magnetic field aligned. The pressure tensor becomes field aligned as one moves away from the equator, but here there is a small region in which the guiding center equations remain inadequate because the two perpendicular components of the pressure tensor are unequal. The significance of nonguiding center motion to substorm processes then was examined. One mechanism that may disrupt a thin cross-tail current sheet involves field changes that cause ions to begin following chaotic orbits. The lowest-altitude chaotic region, characterized by an adiabaticity parameter kappa approx. equal to 0.8, is especially important. The average cross-tail particle drift is slow, and we were unable to generate a thin current sheet using such ions. Therefore, any process that tends to create a thin current sheet in a region with kappa approaching 0.8 may cause the cross-tail current to get so low that it becomes insufficient to support the lobes. A different limit may be important in resonant orbit regions of a thin current sheet because particles reach a maximum cross-tail drift velocity. If the number of ions per unit length decreases as the tail is stretched, this part of the plasma sheet also may become unable to carry the cross-tail current needed to support the lobes. Thin sheets are needed for both resonant and chaotic orbit mechanisms because the distribution function must be highly structured. A description of current continuity is included to show how field aligned currents can evolve during the transition from a two-dimensional (2-D) to a 3-D configuration.

Kaufmann, Richard L.↗

Exponential acceleration of macroscopic quantum tunneling in a Floquet Ising model

The exponential suppression of macroscopic quantum tunneling (MQT) in the number of elements to be reconfigured is an essential element of broken symmetry phases. This suppression is also a core bottleneck in quantum algorithms, such as traversing an energy landscape in optimization, and adiabatic state preparation more generally. In this work, we demonstrate exponential acceleration of MQT through Floquet engineering with the application of a uniform, high frequency transverse drive field. Using the ferromagnetic phase of the transverse field Ising model in one and two dimensions as a prototypical example, we identify three phenomenological regimes as a function of drive strength. For weak drives, the system exhibits exponentially decaying tunneling rates but robust magnetic order; in the crossover regime at intermediate drive strength, we find polynomial decay of tunnelling alongside vanishing magnetic order; and at very strong drive strengths both the Rabi frequency and time-averaged magnetic order are approximately constant with increasing system size. We support these claims with extensive full wavefunction and tensor network numerical simulations, and theoretical analysis. An experimental test of these results presents a technologically important and novel scientific question accessible on NISQ-era quantum computers.

Grattan, George↗

Elemental Fractionation During Rapid Accretion of the Moon Triggered by a Giant Impact

Recently, Ida et al. made an N-body simulation of lunar accretion from a protolunar disk formed by a giant impact. One of their important conclusions is that the accretion time of the Moon is as short as one month. Such rapid accretion is a necessary consequence of the high surface density of a lunar mass disk accreting just beyond the Roche limit (about 3Re); the Safronov accretion time (a few days) is even shorter. The energy of accretion always exceeds the gravitational binding energy of newly arriving matter. Hence, without an energy sink, the accreting body is thermally unstable. For the Earth and other planets, radiation acts as the sink. However, in such a short accretion time, the Moon cannot radiate the accretional energy. Even radiating at a silicate cloudtop temperature of roughly 2000 K, it would take more than 100 yr to radiatively cool the Moon. The plausible alternative heat sinks are heat capacity, latent heat of vaporization, and thermal escape of the gas to space (i.e., hydrodynamic blowoff). The latter becomes plausible for the Moon because the scale height at 2000 K (about 300 km) is a significant fraction of the lunar radius. The early stages of lunar (or "lunatesimal") growth release relatively little energy and can occur simply by heating the material, especially if the accreting material is originally cold. However, the material is unlikely to be cold, because the disk itself is hot and cooling time is long, while the lunar accretion time iss very short. Therefore, the moon is likely to accrete condensed material just after it condenses. Accordingly, the newly accreted material will be on the verge of vaporization and will have very little heat capacity to spare. The immediate heat sink is the latent heat of vaporization. Most of the vapor will escape from the moon, because the thermal energy in the gas can be used to drive escape. However, vaporization is generally incomplete. the latent heat of vaporization exceeds the energy of accretion. Viewed globally, the accretional energy is about half the energy required to vaporize the entire Moon. Thus to first approximation, half of the Moon-forming material can be vaporized and lost during accretion. During this process, we would expect preferential loss of relatively volatile elements. Escape will retard the rate of accretion. To test these ideas, we computed detailed models of the thermal state of the Moon during accretion. We pay special attention to the structure of the silicate atmosphere and its loss rate by calculating the chemical species at equilibrium. We used the PHEQ program which includes 12 elements (H,O,C,Mg,Si,Fe,Ca, Al, Na,Ti, and N.) and 272 compounds (including ionic compounds). Because of the large heats of vaporization and ionization, the adiabatic atmosphere is nearly isothermal and massive escape is expected. The pressure of the atmosphere is determined by the balance between vaporization of a accreting material and escape. If the accretion time is one month, a 0.3 bar atmosphere is expected. Elemental fractionation depends strongly on the temperature of the accreting material. The initial temperature of the material can be estimated from the condition of gravitational instability in the protolunar disk. As shown by Ida et al, accretion starts when gravitational instability occurs when more than 99% of the material condenses. At this point, all of Ca, Al, Si, Mg, and Fe, and 95% of Na (probably K also), are in condensed phases. If the moon is formed from the accretion of such material, volatile elements such as Na, and K are retained by the moon only early in accretion. At later times, K and Na are lost and a fraction of the MG, Si and Fe is lost. However, refractory elements such as Ca and Al are retained and so achieve a mild degree (factor 2) of superabundance.

Abe, Y.↗

Transient mixing behavior of a supercritical fluid injected into supercritical and subcritical environments

Transient, near-field mixing of a supercritical fluid injected into supercritical and subcritical environments is investigated experimentally. Fuel (propane) from a 100 μm diameter Bosch diesel injector is injected into a high-pressure test rig with a steady flow of heated, low-velocity nitrogen to control the ambient pressure and temperature. The fuel supply pressure is approximately 190 bar, allowing short injections lasting 3–4 ms to be performed in quick succession to acquire statistically significant datasets to represent the unsteady fuel–air mixing. Two test cases with the same ambient density and injector inlet pressure and temperature were chosen in an attempt to isolate effects of mixture phase behavior. Microscopic rainbow Schlieren deflectometry (RSD) is used for visualization and quantitative data acquisition. In the case of the supercritical ambient, the injector inlet and exit and ambient conditions are all supercritical. The RSD analysis reveals a mixing process similar to the single phase fluid mixing, though condensation is also observed locally, and can be explained by the adiabatic mixing theory presented in the paper. For the subcritical ambient case, the injector inlet is supercritical, but the injector exit and ambient conditions are both subcritical. The results show a rapid expansion in the near-field promoting mixing between propane and nitrogen. At downstream locations, significant condensation is also observed as predicted by the adiabatic mixing theory. As a result, the supercritical ambient case exhibited longer and narrower jet penetration compared to the subcritical ambient case.

02 PETROLEUM↗

Self-Consistent Simulation of the Magnetotail

Kinetic effects may be critically important at various stages of magnetotail evolution prior to and after the break up of substorms. The quasi-adiabaticity of ion dynamics results in the appearance of pronounced off diagonal terms in the pressure tensor (an effect that cannot be reproduced using the standard MHD (magnetohydrodynamic) approach). Only the appearance of these terms could support the formation of thin current sheets that recent experimental data have repeatedly indicated. We present in this paper the results of a self-consistent large-scale kinetic simulation of the Earth's magnetotail. We launch ions from the plasma mantle and self-consistently adjust the background magnetic field by computing the perturbations magnetic field from the particle current. We find that a thin current sheet with a half-thickness of approximately 0.3 R(sub E) quickly develops in the simulation. Also, the magnetotail evolves into a periodic state where the location of the X-line in the tail oscillates between x approximately 40 R(sub E) and x approximately 70 R(sub E) on a 5-minute time scale. Evidence of this periodic motion is also seen in the auroral precipitation profile of ions.

Peroomian, Vahe↗

Interior structure of Saturn inferred from Pioneer 11 gravity data

The structure of Saturn is studied via a fourth-order theory for rotating planets and equations of state for the envelope which depend parametrically on the helium abundance, on the starting temperature for the adiabat, and on adopted forms of the pressure-density curve in the region of transition from molecular to metallic hydrogen. Models are constrained by the values of J2 and J4 obtained from the Pioneer-Saturn celestial mechanics experiment. Equations of state are tested by computing Jupiter models, which can now be subjected to a more stringent comparison with observed zonal harmonics. It is found that Saturn has a low-density hydrogen-helium envelope with no evidence for enhancement of H2O, CH4, or other abundant compounds. Such compounds are presumably located near the core. The helium mass abundance for Saturn's envelope appears to be in the range of approximately 0.12 to 0.19, but this result is very model-dependent. The helium abundance in the envelope of Jupiter is apparently very similar to that of Saturn.

Hubbard, W. B.↗

An approximate Riemann solver for hypervelocity flows

We describe an approximate Riemann solver for the computation of hypervelocity flows in which there are strong shocks and viscous interactions. The scheme has three stages, the first of which computes the intermediate states assuming isentropic waves. A second stage, based on the strong shock relations, may then be invoked if the pressure jump across either wave is large. The third stage interpolates the interface state from the two initial states and the intermediate states. The solver is used as part of a finite-volume code and is demonstrated on two test cases. The first is a high Mach number flow over a sphere while the second is a flow over a slender cone with an adiabatic boundary layer. In both cases the solver performs well.

Jacobs, Peter A.↗

Simulation of proton radiation belt formation during the March 24, 1991 SSC

The rapid formation of a new proton radiation belt at L approximately = 2.5 following the March 24, 1991 Storm Sudden Commencement (SSC) observed at the Combined Release and Radiation Effects Satellite (CRRES) satellite is modeled using a relativistic guiding center test particle code. The SSC is modeled by a bipolar electric field and associated compression and relaxation in the magnetic field, superimposed on a dipole magnetic field. The source population consists of both solar and trapped inner zone protons. The simulations show that while both populations contribute to drift echoes in the 20-80 MeV range, primary conditions is from the solar protons. Proton acceleration by the SSC differs from relativistic electron acceleration in that different source populations contribute and nonrelativistic conservation of the first adiabatic invariation leads to greater energization of protons for a given decrease in L. Model drift echoes and flux distribution in L at the time of injection compare well with CRRES observations.

Hudson, M. K.↗

A 10 Kelvin 3 Tesla Magnet for Space Flight ADR Systems

Many future space flight missions are expected to use adiabatic demagnetization refrigerators (ADRs) to reach detector operating temperatures well below one Kelvin. The goal is to operate each ADR with a mechanical cooler as its heat sink, thus avoiding the use of liquid cryogens. Although mechanical coolers are being developed to operate at temperatures of 6 Kelvin and below, there is a large efficiency cost associated with operating them at the bottom of their temperature range. For the multi-stage ADR system being developed at Goddard Space Flight Center, the goal is to operate with a 10 Kelvin mechanical cooler heat sink. With currently available paramagnetic materials, the highest temperature ADR stage in such a system will require a magnetic field of approximately three Tesla. Thus the goal is to develop a small, lightweight three Tesla superconducting magnet for operation at 10 Kelvin. It is important that this magnet have a low current/field ratio. Because traditional NbTi magnets do not operate safely above about six Kelvin, a magnet with a higher Tc is required. The primary focus has been on Nb3Sn magnets. Since standard Nb3Sn wire must be coated with thick insulation, wound on a magnet mandrel and then reacted, standard Nb,Sn magnets are quite heavy and require high currents Superconducting Systems developed a Nb3Sn wire which can be drawn down to small diameter, reacted, coated with thin insulation and then wound on a small diameter coil form. By using this smaller wire and operating closer to the wire s critical current, it should be possible to reduce the mass and operating current of 10 Kelvin magnets. Using this "react-then-wind" technology, Superconducting Systems has produced prototype 10 Kelvin magnets. This paper describes the development and testing of these magnets and discusses the outlook for including 10 Kelvin magnets on space-flight missions.

Tuttle, Jim↗

Approximate Quantum Codes From Long Wormholes

We discuss families of approximate quantum error correcting codes which arise as the nearly-degenerate ground states of certain quantum many-body Hamiltonians composed of non-commuting terms. For exact codes, the conditions for error correction can be formulated in terms of the vanishing of a two-sided mutual information in a low-temperature thermofield double state. We consider a notion of distance for approximate codes obtained by demanding that this mutual information instead be small, and we evaluate this mutual information for the SYK model and for a family of low-rank SYK models. After an extrapolation to nearly zero temperature, we find that both kinds of models produce fermionic codes with constant rate as the number, N , of fermions goes to infinity. For SYK, the distance scales as N 1 / 2 , and for low-rank SYK, the distance can be arbitrarily close to linear scaling, e.g. N .99 , while maintaining a constant rate. We also consider an analog of the no low-energy trivial states property which we dub the no low-energy adiabatically accessible states property and show that these models do have low-energy states that can be prepared adiabatically in a time that does not scale with system size N . We discuss a holographic model of these codes in which the large code distance is a consequence of the emergence of a long wormhole geometry in a simple model of quantum gravity.

Physics↗

Exact Factorization Adventures: A Promising Approach for Non-Bound States

Modeling the dynamics of non-bound states in molecules requires an accurate description of how electronic motion affects nuclear motion and vice-versa. The exact factorization (XF) approach offers a unique perspective, in that it provides potentials that act on the nuclear subsystem or electronic subsystem, which contain the effects of the coupling to the other subsystem in an exact way. We briefly review the various applications of the XF idea in different realms, and how features of these potentials aid in the interpretation of two different laser-driven dissociation mechanisms. We present a detailed study of the different ways the coupling terms in recently-developed XF-based mixed quantum-classical approximations are evaluated, where either truly coupled trajectories, or auxiliary trajectories that mimic the coupling are used, and discuss their effect in both a surface-hopping framework as well as the rigorously-derived coupled-trajectory mixed quantum-classical approach.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Approximate Riemann solver for hypervelocity flows

We describe an approximate Riemann solver for the computation of hypervelocity flows in which there are strong shocks and viscous interactions. The scheme has three stages, the first of which computes the intermediate states assuming isentropic waves. A second stage, based on the strong shock relations, may then be invoked if the pressure jump across either wave is large. The third stage interpolates the interface state from the two initial states and the intermediate states. The solver is used as part of a finite-volume code and is demonstrated on two test cases. The first is a high Mach number flow over a sphere while the second is a flow over a flow over a slender cone with an adiabatic boundary layer. In both cases the solver performs well.

Jacobs, P. A.↗