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 289 records · Page 16

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.↗

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.↗

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 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.↗

Our Sun IV: The Standard Model and Helioseismology: Consequences of Uncertainties in Input Physics and in Observed Solar Parameters

Helioseismic frequency observations provide an extremely accurate window into the solar interior; frequencies from the Michaelson Doppler Imager (MDI) on the Solar and Heliospheric Observatory (SOHO) spacecraft, enable the adiabatic sound speed and adiabatic index to be inferred with an accuracy of a few parts in 10(exp 4) and the density with an accuracy of a few parts in 10(exp 3). This has become a Serious challenge to theoretical models of the Sun. Therefore, we have undertaken a self-consistent, systematic study of the sources of uncertainties in the standard solar models. We found that the largest effect on the interior structure arises from the observational uncertainties in the photospheric abundances of the elements, which affect the sound speed profile at the level of 3 parts in 10(exp 3). The estimated 4% uncertainty in the OPAL opacities could lead to effects of 1 part in 10(exp 3); the approximately 5%, uncertainty in the basic pp nuclear reaction rate would have a similar effect, as would uncertainties of approximately 15% in the diffusion constants for the gravitational settling of helium. The approximately 50% uncertainties in diffusion constants for the heavier elements would have nearly as large an effect. Different observational methods for determining the solar radius yield results differing by as much as 7 parts in 10(exp 4); we found that this leads to uncertainties of a few parts in 10(exp 3) in the sound speed int the solar convective envelope, but has negligible effect on the interior. Our reference standard solar model yielded a convective envelope position of 0.7135 solar radius, in excellent agreement with the observed value of 0.713 +/- 0.001 solar radius and was significantly affected only by Z/X, the pp rate, and the uncertainties in helium diffusion constants. Our reference model also yielded envelope helium abundance of 0.2424, in good agreement with the approximate range of 0.24 to 0.25 inferred from helioseismic observations; only extreme Z/X values yielded envelope helium abundance outside this range. We found that other current uncertainties, namely, in the solar age and luminosity, in nuclear rates other than the pp reaction, in the low-temperature molecular opacities, and in the low-density equation of state, have no significant effect on the quantities that can be inferred from helioseismic observations. The predicted pre-main-sequence lithium depletion is uncertain by a factor of 2. The predicted neutrino capture rate is uncertain by approximately 30% for the Cl-27 experiment and by approximately 3% for Ga-71 experiments, while the B-8 neutrino flux is uncertain by approximately 30%.

Boothroyd, Arnold I.↗

The multiple nesting of mesoscale, submesoscale, and nonhydrostatic microscale numerical models for a case of mesocyclogenesis and severe storm development

A fundamental approach involving the integration of two different numerical models over several different spatial mesh lengths is considered. Model I is a hydrostatic, incompressible, and adiabatic primitive (Newtonian) equation set with no boundary layer forcing, i.e., surface heat or moisture flux or external frictional stresses. Typically this model is initialized from conventional radiosonde data and run for a 12-15 hour forecast period over a 40 km horizontal mesh. If phase I of the tornado developmental processes develops, the data is 'frozen' in space and time, interpolated to a 20 km mesh, and integrated for a shorter time period (3 hours) with a shorter time step. Model II is a nonhydrostatic compressible 18-level complement of model I. If model II indicates the development of the intensification of phase V dynamical processes, the data is frozen and interpolated to finer mesh lengths until the tornadic circulation system is approximated

Kaplan, M. L.↗

Transformation toughened ceramics for the heavy duty diesel engine technology program, phase 2

The objective of this program is to develop an insulating structural ceramic for application in a heavy duty adiabatic diesel engine. The approach is to employ transformation toughening (TT) by additions of zirconia-hafnia solid solution (ZHSS). The feasibility of using ZHSS as a toughening agent in mullite and alumina has been demonstrated in Phase 1 of this work. Based on Phase 1 results, a decision was made to concentrate the Phase 2 effort on process optimization of the TT mullite. A strong factor in that decision was the low thermal conductivity and high thermal shock resistance of the mullite. Results of the Phase 2 effort indicate that optimum toughening of mullite by additions of ZHSS is difficult to achieve due to apparent sensitivity to morphology. The 48 ksi room temperature modulus-of-rupture (MOR) achieved in selected specimens is approximately 50% of the original strength target. The MOR deteriorated to 34 ksi at 800 C.

Musikant, S.↗

Van Allen Probes Observations of Magnetic Field Dipolarization and Its Associated O+ Flux Variations in the Inner Magnetosphere at L 6.6

We investigate the magnetic field dipolarization in the inner magnetosphere and its associated ion flux variations, using the magnetic field and energetic ion flux data acquired by the Van Allen Probes. From a study of 74 events that appeared at L= 4.5-6.6 between 1 October 2012 and 31 October 2013, we reveal the following characteristics of the dipolarization in the inner magnetosphere: (1) its time scale is approximately 5 min; (2) it is accompanied by strong magnetic fluctuations that have a dominant frequency close to the O+ gyrofrequency; (3) ion fluxes at 20-50 keV are simultaneously enhanced with larger magnitudes for O+ than for H+; (4) after a few minutes of the dipolarization, the flux enhancement at 0.1-5 keV appears with a clear energy-dispersion signature only for O+; and (5) the energy-dispersed O+ flux enhancement appears in directions parallel or antiparallel to the magnetic field. From these characteristics, we discuss possible mechanisms that can provide selective acceleration to O+ ions at > 20 keV. We conclude that O+ ions at L= 5.4-6.6 undergo nonadiabatic local acceleration caused by oscillating electric field associated with the magnetic fluctuations and/or adiabatic convective transport from the plasma sheet to the inner magnetosphere by the impulsive electric field. At L= 4.5-5.4, however, only the former acceleration is plausible. We also conclude that the field-aligned energy-dispersed O+ ions at 0.1-5 keV originate from the ionosphere and are extracted nearly simultaneously to the onset of the dipolarization.

Magnetosphere↗

A Magnetoresistive Heat Switch for the Continuous ADR

In compensated elemental metals at low temperature, a several Tesla field can suppress electronic heat conduction so thoroughly that heat is effectively carried by phonons alone. In approximately one mm diameter single crystal samples with impurity concentrations low enough that electron conduction is limited by surface scattering, the ratio of zerofield to high-field thermal conductivity can exceed ten thousand. We have used this phenomenon to build a compact, solid-state heat switch with no moving parts and no enclosed fluids. The time scale for switching states is limited by time scale for charging the magnet that supplies the controlling field. Our design and fabrication techniques overcome the difficulties associated with manufacturing and assembling parts from single crystal tungsten. A clear disadvantage of the magnetoresistive switch is the mass and complexity of the magnet system for the controlling field. We have discovered a technique of minimizing this mass and complexity, applicable to the continuous adiabatic demagnetization refrigerator.

Canavan, E. R.↗

Boundary Layer Transition on Slender Cones in Conventional and Low Disturbance Mach 6 Wind Tunnels

An experimental investigation was conducted on a 5-degree half-angle cone and a 5-degree half-angle flared cone in a conventional Mach 6 wind tunnel to examine the effects of facility noise on boundary layer transition. The influence of tunnel noise was inferred by comparing transition onset locations determined from the present test to that previously obtained in a Mach 6 low disturbance quiet tunnel. Together, the two sets of experiments are believed to represent the first direct comparison of transition onset between a conventional and a low disturbance wind tunnel using a common test model and transition detection technique. In the present conventional hypersonic tunnel experiment, separate measurements of heat transfer and adiabatic wall temperatures were obtained on the conical models at small angles of attack over a range of Reynolds numbers, which resulted in laminar, transitional, and turbulent flow. Smooth model turbulent heating distributions are compared to that obtained with transition forced via discrete surface roughness. The model nosetip radius was varied to examine the effects of bluntness on transition onset. Despite wall to total temperature differences between the transient heating measurements and the adiabatic wall temperature measurement, the two methods for determining sharp cone transition onset generally yielded equivalent locations. In the 'noisy' mode of the hypersonic low disturbance tunnel, transition onset occurred earlier than that measured in the conventional hypersonic tunnel, suggesting higher levels of freestream acoustic radiation relative to the conventional tunnel. At comparable freestream conditions, the transition onset Reynolds number under low disturbance conditions was a factor of 1.3 greater than that measured on flared cone in the LaRC conventional hypersonic tunnel and a factor of 1.6 greater that the flared cone run in the low disturbance tunnel run 'noisy'. Navier-Stokes mean flow computations and linear stability analysis were conducted to assess the experimental results and have indicated N factors associated with sharp flared cone transition onset to be approximately a factor of 2 lower than that inferred from the corresponding low disturbance tunnel measurements.

Horvath, Thomas J.↗

Determination of Particle Sources for a Geotail Distribution Function Observed on May 23, 1995

On May 23, 1995, the Comprehensive Plasma Instrumentation (CPI) onboard the Geotail spacecraft observed a complex and structured ion distribution function near the magnetotail midplane at x approximately -10 R(sub E). On the same day, the Wind spacecraft observed a very high density (approximately 40/cubic cm) solar wind and an interplanetary magnetic field (IMF) that was predominantly northward but had several southward turnings. We have inferred the sources of the ions in this distribution function by following approximately 90,000 ion trajectories backward in time using time-dependent electric and magnetic fields obtained from a global MHD (magnetohydrodynamic) simulation. Wind data were used as input for the MHD model. We found that three sources contributed to this distribution: the ionosphere, the plasma mantle which had near-Earth and distant tail components, and the low latitude boundary layer (LLBL). Moreover, distinct structures in the low energy part of the distribution function were found to be associated with individual sources. Structures near 0 deg pitch angle were made up of either ionospheric or plasma mantle ions, while structures near 90 deg pitch angle were dominated by ions from the LLBL source. Particles that underwent nonadiabatic acceleration were numerous in the higher energy part of the ion distribution function, whereas ionospheric and LLBL ions were mostly adiabatic. A large proportion of the near-Earth mantle ions underwent adiabatic acceleration, while most of the distant mantle ions experienced nonadiabatic acceleration.

Ashour-Abdalla, M.↗