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 163 records · Page 9

Models of the Cartwheel ring galaxy: Spokes and starbursts

Recent observations of this famous ring galaxy, including optical and near-infrared CCD surface photometry, and VLA radio continuum and 21 cm line mapping (Higdon 1992b, in prep.), have inspired a renewed modeling effort. Toomre's (1978, in The Large-scale Structure of the Universe, eds. Longair and Einasto) series of restricted three-body simulations demonstrated how the multiple rings could be produced in a nearly head-on galaxy collision. New models with a halo-dominated potential based on the 21 cm rotation curve are able to reproduce such details as the spacing between rings, ring widths, offset of the nucleus, and several kinematical features, thus providing strong support for the collisional theory. The new observations have shown there are little or no old stars in Cartwheel; it may consist almost entirely of gas and stars produced as a result of compression in the ring wave. To model this process Smooth Particle Hydrodynamics (SPH) simulations of the Cartwheel disk have been performed. Fixed gravitational potentials were used to represent the Cartwheel and a roughly 30 percent mass collision partner. The interaction dynamics was treated as in the usual restricted three-body approximation, and the effects of local self-gravity between disk particles were calculated. We are particularly interested in testing the theory that enhanced star formation in waves is the result of gravitational instability in the compressed region (see e.g. Kennicutt 1989, ApJ 344, 685). The gas surface density in a number of simulations was initialized to a value slightly below the threshold for local gravitational instability throughout most of the disk. The first ring wave produces relatively modest compressions (a factor of order a few), triggering instability in a narrow range of wavelengths. Self-gravity in the disk is calculated over a comparable range of scales. Simulations were run with isothermal, adiabatic, and adiabatic with radiative cooling characterized by a relatively short timescale. The isothermal approximation is good except in the vicinity of the strong second (inner) ring, and several snapshots from one case are shown in the figure below. Flocculent spiral segments are present before the collision, and these are compressed into dense knots in the ring wave. These knots are likely to be sites of vigorous star formation. In the strong rarefaction behind the outer ring most of the knots are radially stretched and sheared, giving rise to spoke-like features. A few dense knots are evidently very tightly bound, because they retain their coherence and are stretched relatively little through the rarefaction. This is in accord with evidence for continuing star formation in some spokes (Marcum, Appleton and Higdon 1992). The number and spacing of spokes is a direct function of the scale of the gravitational instability in the disk. Thus, the gravitational instability theory, together with the hypothesis that massive stars are only formed in dense knots of gas, can account for most of the distinct morphology of the Cartwheel.

Struck-Marcell, Curtis↗

Approximate Solution for Choked Flow in Gas Seal Pads

Previous analyses of high pressure seals have considered adiabatic flow with friction but neglected effects of seal rotation. Most of this work analyzed a one-dimensional flow field. This works well to calculate stiffness and leakage of full circular seals, either face seals or annular ring seals. However, it cannot provide accurate results for a rectangular seal pad with its strongly two-dimensional flow field and its reliance on hydrodynamic forces to maintain a full fluid film. On the other hand, solutions of Reynolds lubrication equation have been obtained for the two-dimensional flow in a seal pad. But these solutions do not account for choking which occurs at high seal pressure ratios, nor do they consider the pressure loss that occurs in the entrance region of the flow field. The aim of the present work is to build on the Reynolds equation solution by use of an approximate choked flow analysis. This will account for the pressure losses in the flow entrance region, ensure that fluid velocities remain subsonic, and enable fluid inertial effects within the pad film to be accounted for. Results show that, in general, fluid inertia acts to decrease pad film load capacity and leakage, and increase film stiffness.

Fleming, David P.↗

Adiabatic evolution of trapped particle pitch angle distributions during a storm main phase

Equatorial pitch angle distributions of over 200-keV ions obtained near L = 4 by Explorer 45 during a storm main phase show the development of marked minima at 90 deg pitch angle in direct association with a greater-than-a-factor-of-2 reduction in the equatorial magnetic field magnitude. The evolution of these pitch angle distributions can be quantitatively explained as a result of conservation of the first two adiabatic invariants provided the magnitude of the magnetic field decrease was approximately constant along field lines.

Lyons, L. R.↗

A two-fluid approximation for calculating the cosmic microwave background anisotropies

We present a simplified treatment for calculating the cosmic microwave background anisotropy power spectrum in adiabatic models. It consists of solving for the evolution of a two-fluid model until the epoch of recombination and then integrating over the sources to obtain the cosmic microwave background (CMB) anisotropy power spectrum. The approximation is useful both for a physical understanding of CMB anisotropies as well as for a quantitative analysis of cosmological models. Comparison with exact calculations shows that the accuracy is typically 10%-20% over a large range of angles and cosmological models, including those with curvature and cosmological constant. Using this approximation we investigate the dependence of the CMB anisotropy on the cosmological parameters. We identify six dimensionless parameters that uniquely determine the anisotropy power spectrum within our approximation. CMB experiments on different angular scales could in principle provide information on all these parameters. In particular, mapping of the Doppler peaks would allow an independent determination of baryon mass density, matter mass density, and the Hubble constant.

Seljak, Uros↗

Approximation of nearly-periodic symplectic maps via structure-preserving neural networks

A continuous-time dynamical system with parameter ε is nearly-periodic if all its trajectories are periodic with nowhere-vanishing angular frequency as ε approaches 0. Nearly-periodic maps are discrete-time analogues of nearly-periodic systems, defined as parameter-dependent diffeomorphisms that limit to rotations along a circle action, and they admit formal U(1) symmetries to all orders when the limiting rotation is non-resonant. For Hamiltonian nearly-periodic maps on exact presymplectic manifolds, the formal U(1) symmetry gives rise to a discrete-time adiabatic invariant. In this paper, we construct a novel structure-preserving neural network to approximate nearly-periodic symplectic maps. This neural network architecture, which we call symplectic gyroceptron, ensures that the resulting surrogate map is nearly-periodic and symplectic, and that it gives rise to a discrete-time adiabatic invariant and a long-time stability. This new structure-preserving neural network provides a promising architecture for surrogate modeling of non-dissipative dynamical systems that automatically steps over short timescales without introducing spurious instabilities.

97 MATHEMATICS AND COMPUTING↗

A New Analytical Expression for Estimating the Adiabatic Temperature Rise in Lithium-Ion Batteries During High-Power Pulses

A coupled, thermal-electrochemical model is used to explain why the Joule heating assumption (i.e., I 2 R) does not provide a good representation of the temperature rise during high power pulses in lithium-ion batteries, even in cases where the reversible heat generation can be neglected. The poor agreement occurs because the internal resistance changes during the pulse due to the opposing effects of mass transport in the electrolyte, which raises the resistance, and heat generation (temperature rise), which lowers the resistance. These insights are used to propose a new analytical expression for predicting the temperature rise during adiabatic pulses with limited experimental and physical data. The expression accounts for thermal effects using a Taylor series expansion of an Arrhenius-type equation. It accounts for transport effects using an approximate solution to the one-dimensional diffusion equation. The new expression is shown to accurately estimate the simulated, adiabatic temperature rise across a range of loadings (1 to 4 mAh cm -2 ) and C-rates (1C to 10C) for cells containing LiNi 0.5 Mn 0.3 Co 0.2 O 2 positive electrodes and graphite negative electrodes. In conclusion, it is also shown to accurately estimate the experimental, adiabatic temperature rise measured for cells with LiFePO 4 positive electrodes and mesocarbon negative electrodes with minimal changes to the fitted parameters.

25 ENERGY STORAGE↗

Thermally generated magnetic fields in laser-driven compressions and explosions

The evolution of thermally generated magnetic fields in a plasma undergoing a nearly spherically symmetric adiabatic compression or expansion is calculated. The analysis is applied to obtain approximate results for the development of magnetic fields in laser-driven compression and explosion of a pellet of nuclear fuel. Localized sources, such as those occurring at composition boundaries in structured pellets or at shock fronts, give stronger fields than those deriving from smoothly distributed asymmetries. Although these fields may approach 10 million G in the late stages of compression, this is not expected to present difficulties for the compression process. Assuming ignition of a nuclear explosion occurs, the sources become much stronger, and values of approximately 10 billion G are obtained at tamper boundaries assuming a 20% departure from spherical symmetry during the explosion.

Tidman, D. A.↗

Monte-Carlo model of pitch-angle scattering in solar cosmic ray events

Monte Carlo simulations of the propagation of solar cosmic-rays in interplanetary space are reported, including the effects of pitch-angle scattering and adiabatic focusing. Time intensity profiles agree well with the corresponding spatial diffusion approximation for models close to those used by Gombosi and Owens (1980). Monte Carlo simulations of the same model problem yield results that disagree with those of Gombosi and Owens, and the most likely possible source of error in the Monte Carlo simulation is considered to be in the region near the 90 deg pitch angle where the finite step size introduces certain inconsistencies.

Palmer, I. D.↗

Degradation of performance in ICF implosions due to Rayleigh–Taylor instabilities: A Hamiltonian perspective

The Rayleigh–Taylor instability (RTI) is an ubiquitous phenomenon that occurs in inertial-confinement-fusion (ICF) implosions and is recognized as an important limiting factor of ICF performance. To analytically understand the RTI dynamics and its impact on ICF capsule implosions, we develop a first-principle variational theory that describes an imploding spherical shell undergoing RTI. The model is based on a thin-shell approximation and includes the dynamical coupling between the imploding spherical shell and an adiabatically compressed fluid within its interior. Using a quasilinear analysis, we study the degradation trends of key ICF performance metrics (e.g., stagnation pressure, residual kinetic energy, and areal density) as functions of initial RTI parameters (e.g., the initial amplitude and Legendre mode), as well as the 1D implosion characteristics (e.g., the convergence ratio). We compare analytical results from the theory against nonlinear results obtained by numerically integrating the governing equations of this reduced model. Our findings emphasize the need to incorporate polar flows in the calculation of residual kinetic energy and demonstrate that higher convergence ratios in ICF implosions lead to significantly greater degradation of key performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Atmospheric lapse rate regimes and their parameterization

Lapse rates, moist adiabatic lapse rates and the critical lapse rate for baroclinic adjustment are calculated and compared for the mean annual, January and July states in the Northern Hemisphere. The small seasonal changes in hemispheric mean lapse rates can be attributed to counteracting seasonal changes in the strength of moist convection and baroclinic eddies. A major conclusion is that the vertical temperature structure may be well approximated by a radiative-convective equilibrium model with two critical lapse rates, viz. the moist adiabatic lapse rate and the critical lapse rate for baroclinic adjustment.

Stone, P. H.↗

Approximations for the study of drift boundaries in the magnetosphere

An approximate relation between energy and equatorial radial distance for a particle of given equatorial pitch angle moving adiabatically in a dipole magnetic field is used to describe the Alfven layer for particles of arbitrary pitch angle. The convection electric field is assumed to be uniform in the equatorial plane of the magnetosphere, and magnetic field lines are taken to be equipotentials in the region traversed by the magnetospheric particles being studied. Approximate solutions in the high- and low-energy limits are given in a form directly applicable to the interpretation of measurements at fixed particle energy and pitch angle. These results are particularly useful for electrons, since the approximations provide lower bounds to the exact solutions. The results are used to interpret recent Explorer 45 electron measurements of Williams et al. (1974).

Kivelson, M. G.↗

Investigation into the transport properties of planetary interiors through inelastic X-ray scattering experiments and quantum molecular dynamics (Final Technical Report)

This was a two-year research project for the period 08/15/2018 - 08/14/2020 conducted at the University of Nevada, Reno using Pronghorn, a newly built high-performance cluster located at the University. The short-term project goal was to provide computational support for the experimental campaigns at LCLS through data analysis of previous and upcoming scheduled experiments, along with performing state-of-the-art atomistic simulations of dense plasmas. The longer-term goal was to create a new computational high energy density physics group located at the University of Nevada, Reno. Guided by results from past and future high-resolution scattering experiments, we aimed to research and develop non-equilibrium and non-adiabatic atomistic simulations of dense plasmas that went beyond the Born-Oppenheimer approximation

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Pitch Angle Distributions of Solar Energetic Particles and the Local Scattering Properties of the Interplanetary Magnetic Field

An approximate solution of the Fokker-Planck equation containing pitch angle scattering and adiabatic focusing is discussed. For modest focusing effects the omnidirectional density obeys an ordinary diffusion equation with a modified diffusion coefficient. The anisotropic part of the distribution function is properly normalized and split into an even and an odd part. The even part is determined by the ratio between the scattering mean free path and the focusing length and by the degree of polarization of the magnetic field fluctuations. The odd part is determined by the deviation of the pitch angle scattering from isotropic scattering. The method supplies a powerful tool to obtain the local characteristics of pitch angle scattering. It is insensitive to long lasting solar injections and to moderate radial variations of the mean free path. The method is applied to solar particle events observed on Helio-1 and -2.

Beeck, J.↗

Near equality of ion phase space densities at earth, Jupiter, and Saturn

Energetic-ion phase-space density profiles are strikingly similar in the inner magnetospheres of earth, Jupiter, and Saturn for ions of first adiabatic invariant near 100 MeV/G and small mirror latitudes. Losses occur inside L approximately equal to 7 for Jupiter and Saturn and inside L approximately equal to 5 at earth. At these L values there exist steep plasma-density gradients at all three planets, associated with the Io plasma torus at Jupiter, the Rhea-Dione-Tethys torus at Saturn, and the plasmasphere at earth. Measurements of ion flux-tube contents at Jupiter and Saturn by the low-energy charged-particle experiment show that these are similar (for O ions at L = 5-9) to those at earth (for protons at L = 2-6). Furthermore, the thermal-ion flux-tube contents from Voyager plasma-science data at Jupiter and Saturn are also very nearly equal, and again similar to those at earth, differing by less than a factor of 3 at the respective L values. The near equality of energetic and thermal ion flux-tube contents at earth, Jupiter, and Saturn suggests the possibility of strong physical analogies in the interaction between plasma and energetic particles at the plasma tori/plasma sheets of Jupiter and Saturn and the plasmasphere of earth.

Cheng, A. F.↗

Adiabatic invariants in stellar dynamics. 2: Gravitational shocking

A new theory of gravitational shocking based on time-dependent perturbation theory shows that the changes in energy and angular momentum due to a slowly varying disturbance are not exponentially small for stellar dynamical systems in general. It predicts significant shock heating by slowly varying perturbations previously thought to be negligible according to the adiabatic criterion. The theory extends the scenarios traditionally computed only with the impulse approximation and is applicable to a wide class of disturbances. The approach is applied specifically to the problem of disk shocking of star clusters.

Weinberg, Martin D.↗

Consequences of magnetotail ion dynamics

The trajectories of a large ensemble of particles are calculated in a modified Tsyganenko magnetic field model with a uniform cross-tail electric field. The model magnetotail can be divided into several distinct dynamical regimes of ion motion. Near Earth, where the field lines are dipolar the adiabatic formalism is adequate. In the mid-tail and distant tail, guiding-center theory breaks down and must be replaced by a quasi-adiabatic formalism. There is an important transition region between the adiabatic and quasi-adiabatic regions where ion trajectories become more complicated and no simple analytical description holds. This wall region is characterized by rapid ion acceleration and a major loss of particles to the dusk flank. The moments of the ion distribution function are constructued from the ion trajectories, including density, temperature, and pressure in the x-z and x-y planes. In the noon-midnight meridian plane, parameters are relatively constant except near the Earth, while the x-y plots show strong gradients across the magnetotail. Magnetotail plasma convects earthward, drifts toward dusk, and is squeezed out of the tail in the near-Earth region. A thin current sheet forms in the quasi-adiabatic region, and the pressure tensor has significant off-diagonal terms at its edges. These terms are the result of quasi-adiabatic ion trajectories which lead to azimuthally asymmetric distribution functions capable of maintaining approximate stress balance across the current sheet. Simplified analytical descriptions provide further physical insight into ion dynamics that are observed.

Ashour-Abdalla, Maha↗

Constraints on Mantle Plume Melting Conditions in the Martian Mantle Based on Improved Melting Phase Relationships of Olivine-Phyric Shergottite Yamato 980459

Martian meteorite Yamato 980459 (hereafter Y98) is an olivine-phyric shergottite that has been interpreted as closely approximating a martian mantle melt [1-4], making it an important constraint on adiabatic decompression melting models. It has long been recognized that low pressure melting of the Y98 composition occurs at extremely high temperatures relative to martian basalts (1430 degC at 1 bar), which caused great difficulties in a previous attempt to explain Y98 magma generation via a mantle plume model [2]. However, previous studies of the phase diagram were limited to pressures of 2 GPa and less [2, 5], whereas decompression melting in the present-day martian mantle occurs at pressures of 3-7 GPa, with the shallow boundary of the melt production zone occurring just below the base of the thermal lithosphere [6]. Recent experimental work has now extended our knowledge of the Y98 melting phase relationships to 8 GPa. In light of this improved petrological knowledge, we are therefore reassessing the constraints that Y98 imposes on melting conditions in martian mantle plumes. Two recently discovered olivine- phyric shergottites, Northwest Africa (NWA) 5789 and NWA 6234, may also be primary melts from the martian mantle [7, 8]. However, these latter meteorites have not been the subject of detailed experimental petrology studies, so we focus here on Y98.

Kiefer, Walter S.↗

Collisionless tearing in a field-reversed sheet pinch assuming nonparallel propagation

The problem of collisionless linear tearing is examined assuming a wave vector with a component normal to the equilibrium field. The geometry is defined and the general form of the linear dispersion equation is calculated. The linear theory results when k is parallel to B are reviewed, and Ampere's law is calculated for the external adiabatic region when k times B does not equal zero, using two-fluid theory. A solution is obtained for the approximate form of the perturbed currents and vector potential assuming quasi-parallel k. The resonant current contributions within the singular layer are calculated, obtaining an estimate of the dispersion equation. The form of the adiabatic currents within the singular layer is calculated, showing that an x-z current system persists even in the limit k perpendicular to B goes to zero. Finally, the perturbed vector potential solutions across the singular layer are matched to obtain the shape of the complete eigenfunction.

Quest, K. B.↗