Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixed Layer Depth”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS1998

Monitoring of Volcanogenic CO2-Induced Tree Kills with AVIRIS Image Data at Mammoth Mountain, California

Elevated cold CO2 emissions from the flank of Mammoth Mountain volcano on the southwest rim of the Long Valley Caldera, eastern California, have been the cause of over 100 acres of dead trees in that area since 1990. The source of the CO2 gas is thought to be from one or more magmatic intrusion(s) beneath Mammoth Mountain and is probably related to a period of seismic unrest that began in 1989. The gas rises to the surface probably from depths of a few kilometers, along faults and fracture zones. The gas is at ambient temperature and diffuses from the soil rather than discharging from distinct vents. Typically, soil gas concentrations in tree-kill areas range from 10% to over 90% CO2 by volume, as compared to normal background of < 1% in healthy forest. The gas composition is predominantly CO2 mixed with air (sulfur gases are not elevated), and C and He isotopic ratios are consistent with a magmatic origin for the gas. The total CO2 emission has been estimated at 1200 tons/day, comparable to the emissions at Kilauea. Some of the dead trees are as old as 250 years, suggesting that similar anomalous gas discharge has not occurred over the previous few hundred years. The delta C-13/12 ratio in the Mammoth Mountain CO2 emission averages about -4.5 (PDB standard). This is consistent with a mantle source for the carbon. However, the large volume of the emission suggests that not all of the CO2 is necessarily being generated from the 1989 intrusion. The voluminous gas could be leaking from a vapor-rich zone, capped by an impermeable layer, that was supplied CO2 from degassing of many small magma bodies that intruded beneath the mountain over a period of decades or centuries. Earthquakes in 1989 could have fractured the capping layer and provided pathways for the escape Of CO2 to the surface. Alternatively, some of the CO2 could be derived from contact metamorphism of carbonate rocks intruded by magma. Carbonate-bearing Paleozoic roof pendents crop out in close proximity to Mammoth Mountain. It is possible that similar rocks could occur at depth beneath Mammoth Mountain, and could have contributed CO2 from thermal decomposition caused from recent intrusions. We hope to determine the C-13/12 ratio of a suite of samples to demonstrate if the carbonate rocks could be the source of at least part of the 1990-97 CO2 emission. To better understand the behavior of the CO2 gas, we have used hyperspectral imagery data of Mammoth Mountain acquired from the Airborne Visual/Infrared Imaging Spectrometer (AVIRIS) to map out areas of dead trees. The areas of tree kill have increased in size from about 50 acres in 1994 to about 100 acres in 1997. Tree kill is the major surface manifestation of the carbon dioxide flux at Mammoth Mountain, is widely dispersed, and has been cursorily mapped by regular field mapping techniques in the area. Initial investigations using airborne digital imagery from the Thematic Mapper Simulator (NS001) and AVIRIS instruments have shown extremely encouraging results for complete delineation of the vegetation anomalies. The most successful maps (when compared with ground truth) were developed using AVIRIS data with spectral angle mapper and matched filter algorithms with a data set that was reduced to maximum variance via the minimum noise fraction transformation. The result of this work is a series of maps that show the tree kill areas occurring in an halo-pattern surrounding the base of Mammoth Mountain. We are applying these same techniques to earlier AVIRIS images of Mammoth Mountain to examine the progression of the tree kill areas over time. Temporal maps of the tree kill areas may assist in constructing a picture of the structure beneath Mammoth Mountain.

Hausback, Brian P.↗
NASA NTRS1992

Photosynthetic characteristics and estimated growth rates indicate grazing is the proximate control of primary production in the equatorial Pacific

Macronutrients persist in the surface layer of the equatorial Pacific because the production of phytoplankton is limited; the nature of this limitation has yet to be resolved. Measurements of photosynthesis as a function of irradiance (P-I) provide information on the control of primary productivity, a question of great biogeochemical importance. Accordingly, P-I was measured in the equatorial Pacific along 150 deg W, during February-March 1988. Diel variability of P-I showed a pattern consistent with nocturnal vertical mixing in the upper 20 m followed by diurnal stratification, causing photoinhibition near the surface at midday. Otherwise, the distribution of photosynthetic parameters with depth and the stability of P-I during simulated in situ incubations over 2 days demonstrated that photoadaptation was nearly complete at the time of sampling: photoadaptation had not been effectively countered by upwelling or vertical mixing. Measurements of P-I and chlorophyll during manipulations of trace elements showed that simple precautions to minimize contamination were sufficient to obtain valid rate measurements and that the specific growth rates of phytoplankton were fairly high in situ, a minimum of 0.6/d. Diel variability of beam attenuation also indicated high specific growth rates of phytoplankton and a strong coupling of production with grazing. It appears that grazing is the proximate control on the standing crop of phytoplankton. Nonetheless, the supply of a trace nutrient such as iron might ultimately regulate productivity by influencing species composition and food-web structure.

Cullen, John J.↗
NASA NTRS2016

Elevated Aerosol Layers and Their Radiative Impact over Kanpur During Monsoon Onset Period

Accurate information about aerosol vertical distribution is needed to reduce uncertainties in aerosol radiative forcing and its effect on atmospheric dynamics. The present study deals with synergistic analyses of aerosol vertical distribution and aerosol optical depth (AOD) with meteorological variables using multisatellite and ground-based remote sensors over Kanpur in central Indo-Gangetic Plain (IGP). Micro-Pulse Lidar Network-derived aerosol vertical extinction (sigma) profiles are analyzed to quantify the interannual and daytime variations during monsoon onset period (May-June) for 2009-2011. The mean aerosol profile is broadly categorized into two layers viz., a surface layer (SL) extending up to 1.5 km (where sigma decreased exponentially with height) and an elevated aerosol layer (EAL) extending between 1.5 and 5.5 km. The increase in total columnar aerosol loading is associated with relatively higher increase in contribution from EAL loading than that from SL. The mean contributions of EALs are about 60%, 51%, and 50% to total columnar AOD during 2009, 2010, and 2011, respectively. We observe distinct parabolic EALs during early morning and late evening but uniformly mixed EALs during midday. The interannual and daytime variations of EALs are mainly influenced by long-range transport and convective capacity of the local emissions, respectively. Radiative flux analysis shows that clear-sky incoming solar radiation at surface is reduced with increase in AOD, which indicates significant cooling at surface. Collocated analysis of atmospheric temperature and aerosol loading reveals that increase in AOD not only resulted in surface dimming but also reduced the temperature (approximately 2-3 C) of lower troposphere (below 3 km altitude). Radiative transfer simulations indicate that the reduction of incoming solar radiation at surface is mainly due to increased absorption by EALs (with increase in total AOD). The observed cooling in lower troposphere in high aerosol loading scenario could be understood as a dynamical feedback of EAL-induced stratification of lower troposphere. Further, the observed radiative effect of EALs increases the stability of the lower troposphere, which could modulate the large-scale atmospheric dynamics during monsoon onset period. These findings encourage follow-up studies on the implication of EALs to the Indian summer monsoon dynamics using numerical models.

aerosol vertical distribution↗
DOE OSTI2022

Sound velocity and compressibility of melts along the hedenbergite (CaFeSi2O6)-diopside (CaMgSi2O6) join at high pressure : Implications for stability and seismic signature of Fe-rich melts in the mantle

Iron-rich silicate melts play an important role in the magmatic history of the Earth and the 16 Moon. However, their elastic properties at high pressures, especially the sound velocity, are poorly 17 understood. Here we determined the ultrasonic sound velocity for the first time of a hedenbergite 18 (Hd, CaFeSi2O6) melt and a melt mixture of 50 mol% hedenbergite + 50 mol% diopside (Hd50Di50) 19 at high pressure and temperature conditions up to 6 GPa and 2329 K, using high-pressure ultrasonic 20 technique combined with synchrotron radiation in a multi-anvil apparatus. Our results show that 21 Fe can significantly reduce the sound velocity while increase the density of silicate melts. 22 Comparing the Di, Hd, and Hd50Di50 melts, we find that the sound velocity does not mix linearly 23 2 for melts in the Hd-Di join, whereas the density for Hd-Di melts at high pressures can be well-24 described by linear mixing. Combined with melt geometry and melt compaction models, we 25 applied our results to study the stability and seismic signature of Fe-rich silicate melts in the Earth’s 26 upper mantle. For the low-velocity zone (LVZ) in mantle asthenosphere, although the degree of 27 seismic velocity reduction can be explained by the presence of a small amount of partial melt 28 distributed in film/band geometry along grain boundaries, Fe-rich melts formed at this depth are 29 unlikely to be gravitationally stable, but may be dynamically unextractable if the melt supply is 30 continuous. For the low-velocity layer (LVL) above the mantle transition zone, the presence of Fe-31 rich melts (with FeO>~10 wt%) distributed in textural equilibrium with the ambient mantle is a 32 plausible explanation.

diopside↗
DOE OSTI2025

Remote-sensing detectability of airborne Arctic dust

Remote-sensing (RS)-based estimates of Arctic dust are oftentimes overestimated due to a failure in separating out the dust contribution from that of spatially homogeneous clouds or low-altitude cloud-like plumes. A variety of illustrations are given with a particular emphasis on questionable claims of using brightness temperature differences (BTDs) as a signature indicator of Arctic dust transported from mid-latitude deserts or generated by local Arctic sources. While there is little dispute about the presence of both Asian and local dust across the Arctic, the direct RS detectability of airborne dust, as ascribed to satellite (MODIS and AVHRR) measurements of significantly negative brightness temperature differences at 11 and 12 µm (BTD 11–12 ), has been misrepresented in certain cases. While it is difficult to account for all examples of strongly negative BTD 11–12 values in the Arctic, it is unlikely that airborne dust plays a significant role. One much more likely contributor would be water clouds in the Arctic inversion layer. The RS detectability of the impact of Arctic dust (notably due to Arctic dust from local sources) can, however, be of significance. Sustained dust deposition can substantially decrease (visible to shortwave IR) snow and ice reflectance albedo (pan-chromatic reflectance) and the signal measured by satellite sensors. Significantly negative BTD 11–12 values would, however, only represent a limited area near the drainage basin sources according to our event-level case studies. The enhanced ice-nucleating particle (INP) role of local Arctic dust can, for example, induce significant changes in the properties of low-level mixed-phase clouds (cloud optical depth changes <~ 1) that can readily be detected by active and passive RS instruments. It is critical that the distinction between the RS detectability of airborne Arctic dust versus the RS detectability of the impacts of that dust be understood if we are to appropriately parameterize, for example, the radiative forcing influence of dust in this climate-sensitive region.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1992

Chemical differentiation of a convecting planetary interior: Consequences for a one-plate planet such as Venus

Chemically depleted mantle forming a buoyant, refractory layer at the top of the mantle can have important implications for the evolution of the interior and surface. On Venus, the large apparent depths of compensation for surface topographic features might be explained if surface topography were supported by variations in the thickness of a 100-200 km thick chemically buoyant mantle layer or by partial melting in the mantle at the base of such a layer. Long volcanic flows seen on the surface may be explained by deep melting that generates low-viscosity MgO-rich magmas. The presence of a shallow refractory mantle layer may also explain the lack of volcanism associated with rifting. As the depleted layer thickens and cools, it becomes denser than the convecting interior and the portion of it that is hot enough to flow can mix with the convecting mantle. Time dependence of the thickness of a depleted layer may create episodic resurfacing events as needed to explain the observed distribution of impact craters on the venusian surface. We consider a planetary structure consisting of a crust, depleted mantle layer, and a thermally and chemically well-mixed convecting mantle. The thermal evolution of the convecting spherical planetary interior is calculated using energy conservation: the time rate of change of thermal energy in the interior is equated to the difference in the rate of radioactive heat production and the rate of heat transfer across the thermal boundary layer. Heat transfer across the thermal boundary layer is parameterized using a standard Nusselt number-Rayleigh number relationship. The radioactive heat production decreases with time corresponding to decay times for the U, Th, and K. The planetary interior cools by the advection of hot mantle at temperature T interior into the thermal boundary layer where it cools conductively. The crust and depleted mantle layers do not convect in our model so that a linear conductive equilibrium temperature distribution is assumed. The rate of melt production is calculated as the product of the volume flux of mantle into the thermal boundary layer and the degree of melting that this mantle undergoes. The volume flux of mantle into the thermal boundary layer is simply the heat flux divided by amount of heat lost in cooling mantle to the average temperature in the thermal boundary layer. The degree of melting is calculated as the temperature difference above the solidus, divided by the latent heat of melting. A maximum degree of melting is prescribed corresponding to the maximum amount of basaltic melt that the mantle can initially generate. As the crust thickens, the pressure at the base of the crust becomes high enough and the temperature remains low enough for basalt to transform to dense eclogite.

Parmentier, E. M.↗
DOE OSTI2021

Identifying dominant environmental predictors of freshwater wetland methane fluxes across diurnal to seasonal time scales

While wetlands are the largest natural source of methane (CH 4 ) to the atmosphere, they represent a large source of uncertainty in the global CH 4 budget due to the complex biogeochemical controls on CH 4 dynamics. Here we present, to our knowledge, the first multi-site synthesis of how predictors of CH 4 fluxes (FCH4) in freshwater wetlands vary across wetland types at diel, multiday (synoptic), and seasonal time scales. In this work we used several statistical approaches (correlation analysis, generalized additive modeling, mutual information, and random forests) in a wavelet-based multi-resolution framework to assess the importance of environmental predictors, nonlinearities and lags on FCH4 across 23 eddy covariance sites. Seasonally, soil and air temperature were dominant predictors of FCH4 at sites with smaller seasonal variation in water table depth (WTD). In contrast, WTD was the dominant predictor for wetlands with smaller variations in temperature (e.g., seasonal tropical/subtropical wetlands). Changes in seasonal FCH4 lagged fluctuations in WTD by ~17 ± 11 days, and lagged air and soil temperature by median values of 8 ± 16 and 5 ± 15 days, respectively. Temperature and WTD were also dominant predictors at the multiday scale. Atmospheric pressure (PA) was another important multiday scale predictor for peat-dominated sites, with drops in PA coinciding with synchronous releases of CH 4 . At the diel scale, synchronous relationships with latent heat flux and vapor pressure deficit suggest that physical processes controlling evaporation and boundary layer mixing exert similar controls on CH 4 volatilization, and suggest the influence of pressurized ventilation in aerenchymatous vegetation. In addition, 1- to 4-h lagged relationships with ecosystem photosynthesis indicate recent carbon substrates, such as root exudates, may also control FCH4. By addressing issues of scale, asynchrony, and nonlinearity, this work improves understanding of the predictors and timing of wetland FCH4 that can inform future studies and models, and help constrain wetland CH 4 emissions.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS1993

Spectral evidence of size dependent space weathering processes on asteroid surfaces

Most compositional characterizations of the minor planets are derived from analysis of visible and near-infrared reflectance spectra. However, such spectra are derived from light which has only interacted with a very thin surface layer. Although regolith processes are assumed to mix all near-surface lithologic units into this layer, it has been proposed that space weathering processes can alter this surface layer to obscure the spectral signature of the bedrock lithology. It has been proposed that these spectral alteration processes are much less pronounced on asteroid surfaces than on the lunar surface, but the possibility of major spectral alteration of asteroidal optical surfaces has been invoked to reconcile S-asteroids with ordinary chondrites. The reflectance spectra of a large subset of the S-asteroid population have been analyzed in a systematic investigation of the mineralogical diversity within the S-class. In this sample, absorption band depth is a strong function of asteroid diameter. The S-asteroid band depths are relatively constant for objects larger than 100 km and increase linearly by factor of two toward smaller sizes (approximately 40 km). Although the S-asteroid surface materials includes a diverse variety of silicate assemblages, ranging from dunites to basalts, all compositional subtypes of the S-asteroids conform to this trend. The A-, R-, and V-type asteroids which are primarily silicate assemblages (as opposed to the metal-silicate mixtures of most S-asteroids) follow a parallel but displaced trend. Some sort of textural or regolith equilibrium appears to have been attained in the optical surfaces of asteroids larger than about 100 km diameter but not on bodies below this size. The relationships between absorption band depth, spectral slope, surface albedo and body size provide an intriguing insight into the nature of the optical surfaces of the S-asteroids and space weathering on these objects.

Gaffey, M. J.↗
NASA NTRS2009

Evaluation of the OMI Cloud Pressures Derived from Rotational Raman Scattering by Comparisons with other Satellite Data and Radiative Transfer Simulations

In this paper we examine differences between cloud pressures retrieved from the Ozone Monitoring Instrument (OMI) using the ultraviolet rotational Raman scattering (RRS) algorithm and those from the thermal infrared (IR) Aqua/MODIS. Several cloud data sets are currently being used in OMI trace gas retrieval algorithms including climatologies based on IR measurements and simultaneous cloud parameters derived from OMI. From a validation perspective, it is important to understand the OMI retrieved cloud parameters and how they differ with those derived from the IR. To this end, we perform radiative transfer calculations to simulate the effects of different geophysical conditions on the OMI RRS cloud pressure retrievals. We also quantify errors related to the use of the Mixed Lambert-Equivalent Reflectivity (MLER) concept as currently implemented of the OMI algorithms. Using properties from the Cloudsat radar and MODIS, we show that radiative transfer calculations support the following: (1) The MLER model is adequate for single-layer optically thick, geometrically thin clouds, but can produce significant errors in estimated cloud pressure for optically thin clouds. (2) In a two-layer cloud, the RRS algorithm may retrieve a cloud pressure that is either between the two cloud decks or even beneath the top of the lower cloud deck because of scattering between the cloud layers; the retrieved pressure depends upon the viewing geometry and the optical depth of the upper cloud deck. (3) Absorbing aerosol in and above a cloud can produce significant errors in the retrieved cloud pressure. (4) The retrieved RRS effective pressure for a deep convective cloud will be significantly higher than the physical cloud top pressure derived with thermal IR.

Vasilkov, Alexander↗
NASA NTRS1991

Convection in three dimensions with surface plates - Generation of toroidal flow

This work presents numerical calculations of mantle convection that incorporate some of the basic observational constraints imposed by plate tectonics. The model is three-dimensional and includes surface plates; it allows plate velocity to change dynamically according to the forces which result from convection. It is shown that plates are an effective means of introducing a toroidal component into the flow field. After initial transients the plate motion is nearly parallel to transform faults and in the direction that tends to minimize the toroidal flow field. The toroidal field decays with depth from its value at the surface; the poloidal field is relatively constant throughout the layer but falls off slightly at the top and bottom boundaries. Layered viscosity increasing with depth causes the toroidal field to decay more rapidly, effectively confining it to the upper, low-viscosity layer. The effect of viscosity layering on the poloidal field is relatively small, which is attributed to its generation by temperature variations distributed throughout the system. The generation of toroidal flow by surface plates would seem to account for the observed nearly equal energy of toroidal and poloidal fields of plate motions on the earth. A low-viscosity region in the upper mantle will cause the toroidal flow to decay significantly before reaching the lower mantle. The resulting concentration of toroidal flow in the upper mantle may result in more thorough mixing there and account for some of the geochemical and isotopic differences proposed to exist between the upper and lower mantles.

Gable, Carl W.↗
NASA NTRS1980

A model of early lunar differentiation

The presence of anorthosites and the scarcity of complementary mafic rocks in the lunar crust suggest a body of melt with depth of the order of the depth of the crust itself (about 60 km). This depth is consistent with geophysical and geochemical models which call for a depth of melting of the order of 200-500 km. The major and minor element compositions of the ancient pristine rocks require complex igneous processes if the moon accreted homogeneously and has chondritic interelement ratios among the refractory elements. Anorthosites crystallized from magmas that had assimilated plagioclase and a light REE-enriched component. KREEP as well as most of the troctolites and norites formed from relatively primitive magmas that were mixed with a LIL-enriched, Sc-Ti depleted component. Assimilation and mixing processes inferred for pristine rocks are consistent with crystallization from a convecting magma ocean which initially developed a floating anorthosite crust with a plagioclase-saturated boundary layer beneath it, transitional in temperature and composition to the main body of the magma beneath it.

Longhi, J.↗
NASA NTRS1989

Synthesis of H2 in dirty ice mantles by fast ion energy loss: New experimental results increase the relevance of this mechanism

Recent experimental results support the importance of H2 production in molecular clouds by cosmic ray bombardment of the mantles of grains. The formation of molecules different from those originally present in the irradiated layer can be explained by the production of molecular fragments induced by the release of energy if the impinging fast particle. One way of considering the process is in terms of a transiently hot cylinder, initially about 50 A in diameter, that exists around the track of an individual fast ion. Since ice has a relatively low thermal conductivity, energy lost by the ion in the ice layers remains confined around the track for time long enough to be thermalized. The hot cylinder increases in diameter and decreases in temperature on a time scale of 10(exp -11) to 10(exp -10) sec. Molecular fragments that are formed in this high temperature region acquire enough mobility to recombine with different partners, forming new molecules. A Monte Carlo simulation of the interaction between cosmic rays and grain mantles, at various depths in the core of a spherical molecular cloud, was performed. The simulation was continued until 40,000 ions had hit each grain of the type and size chosen. During the performed experiments thin icy films made of H2O and CD4 mixed in the gas phase and deposited on a cold finger at 9 K were irradiated with 1.5 MeV helium beams. Among synthesized molecules were found H2, HD, and D2.

Pirranello, Valero↗
NASA NTRS2004

Lunar and Planetary Science XXXV: Mars Polar Processes: Land and Sky

The titles in this section include: 1) Morphological Evidence for the Large-scale Evolution of Martian North Polar Troughs? 2) New Views of the Martian Polar Regions: The Latest Results from the Mars Odyssey THEMIS Investigation; 3) Observations at the Chasma Australe Re-Entrant, South Polar Region, Mars; 4) Stratigraphic Details of Uppermost Units within South Polar Layered Deposits on Mars; 5) Detailed Geologic Analysis of Part of the South Polar Layered Deposits, Planum Australe, Mars: Part II; 6) SHARAD: Radar Volume Scattering and the Polar Layered Terrains on Mars; 7) Arabia and Memnonia Equatorial Regions with High Content of Water: Data from HEND/Odyssey; 8) Alaskan Permafrost Analogs of Martian Small Valley Networks, Thermokarst, Terrain Softening, Terraces, and Volcanic Craters; 9) Periglacial Processes in the Southern Hemisphere of Mars; 10) Thermal Modeling of Possible Surface Water Ice Deposits in Juventae Chasma; 11) Modeling Depth to Ground Ice on Mars; 12) Ground Ice at the Phoenix Landing Site: A Preflight Assessment; 13) Possible Temperature-related Differences in Slope Angle between the North and South Walls of Coprates Chasma, Mars; 14) Seasonal Variations within Korolev Crater, Mars; 15) South Polar Ar Enhancement as a Tracer for Southern Winter Horizontal Meridional Mixing; 16) An Investigation of the Correlation of Water-Ice and Dust Retrievals via the MGS TES Data Set; 17) Observations of Plasma Waves near Mars and Their Implications for Atmospheric Loss.

Source record↗
DOE OSTI2026

Gauging nexus between topological and fracton phases

Coupled layer constructions are a valuable tool for capturing the universal properties of certain interacting quantum phases of matter in terms of the simpler data that characterizes the underlying layers. In the study of fracton phases, the X-Cube model in 3+1D can be realized via such a construction by starting with a stack of 2+1D Toric Codes and turning on a coupling which condenses a composite "particle-string" object. In a recent work [Phys. Rev. B 112, 125124 (2025)], we have demonstrated that in fact, the particle-string can be viewed as a symmetry defect of a topological 1-form symmetry. In this paper, we study the result of gauging this symmetry in depth. We unveil a rich gauging web relating the X-Cube model to symmetry protected topological (SPT) phases protected by a mix of subsystem and higher-form symmetries, subsystem symmetry fractionalization in the 3+1D Toric Code, and non-trivial extensions of topological symmetries by subsystem symmetries. Here, our work emphasizes the importance of topological symmetries in non-topological, geometric phases of matter.

Anyons↗
NASA NTRS2005

Simulations of Evaporating Multicomponent Fuel Drops

A paper presents additional information on the subject matter of Model of Mixing Layer With Multicomponent Evaporating Drops (NPO-30505), NASA Tech Briefs, Vol. 28, No. 3 (March 2004), page 55. To recapitulate: A mathematical model of a three-dimensional mixing layer laden with evaporating fuel drops composed of many chemical species has been derived. The model is used to perform direct numerical simulations in continuing studies directed toward understanding the behaviors of sprays of liquid petroleum fuels in furnaces, industrial combustors, and engines. The model includes governing equations formulated in an Eulerian and a Lagrangian reference frame for the gas and drops, respectively, and incorporates a concept of continuous thermodynamics, according to which the chemical composition of a fuel is described by use of a distribution function. In this investigation, the distribution function depends solely on the species molar weight. The present paper reiterates the description of the model and discusses further in-depth analysis of the previous results as well as results of additional numerical simulations assessing the effect of the mass loading. The paper reiterates the conclusions reported in the cited previous article, and states some new conclusions. Some new conclusions are: 1. The slower evaporation and the evaporation/ condensation process for multicomponent-fuel drops resulted in a reduced drop-size polydispersity compared to their single-component counterpart. 2. The inhomogeneity in the spatial distribution of the species in the layer increases with the initial mass loading. 3. As evaporation becomes faster, the assumed invariant form of the molecular- weight distribution during evaporation becomes inaccurate.

Bellan, Josette↗
NASA NTRS1997

Planetary Geophysics and Tectonics

Research supported by grant NAGW-1928 has addressed a variety of problems related to planetary evolution. One important focus has been on questions related to the role of chemical buoyancy in planetary evolution with application to both Venus and the Moon. We have developed a model for the evolution of the Moon (Hess and Parmentier, 1995) in which dense, highly radioactive, late stage magma ocean cumulates sink forming a core. This core heats the overlying, chemically layered mantle giving rise to a heated, chemically well-mixed layer that thickens with time. This Mixed layer eventually becomes hot enough and thick enough that its top begins to melt at a pressure low enough that melt is buoyant, thus creating mare basalts from a high pressure source of the correct composition and at an appropriate time in lunar evolution. In work completed during the last year, numerical experiments on convection in a chemically stably stratified fluid layer heated from below have been completed. These results show us how to calculate the evolution of a mixed layer in the Moon, depending on the heat production in the ilmenite- cumulate core and the chemical stratification of the overlying mantle. Chemical stratification of the mantle after its initial differentiation is would trap heat in the deep interior and prevent the rapid rise of plumes with accompanying volcanism. This trapping of heat in the interior can explain the thickness of the lunar lithosphere as a function of time as well as the magmatic evolution. We show that heat transported to the base of the lithosphere at a rate determined by current estimates of radioactivity in the Moon would not satisfy constraints on elastic lithosphere thickness from tectonic feature associated with basin loading. Trapping heat at depth by a chemically stratified mantle may also explain the absence of global compressional features on the surface that previous models predict for an initially hot lunar interior. For Venus, we developed a model in which the chemical buoyancy of crust and a depleted mantle layer stabilizes the lithosphere for long periods of time and provides a mechanism of episodic planetary evolution (Parmentier and Hess, 1992). Continued thickening of a residual depleted mantle layer eventually suppresses pressure release melting and the creation of depleted mantle. Continued cooling then allows the lithosphere to become heavier than the underlying hotter, undepleted mantle. This repeated instability can occur on time scales appropriate for episodic global resurfacing on Venus. We have also examined the role of the gabbro-eclogite phase transformation on crust and lithosphere stability and as a mechanism of crustal recycling in the absence of plate tectonics. Our work thus far concentrates on the scale of instability that would occur due to cooling or crustal thickening associated with horizontal shortening. Whether repeated overturn can explain the evolution of Venus depends in part on whether sufficient heat transfer can occur between overturns and on constraints provided by understanding observed surface features and evolution.

Parmentier, E. M.↗
DOE OSTI2021

The diurnal cycle of entrainment and detrainment in LES of the Southern Ocean driven by observed surface fluxes and waves

Empirical rules for both entrainment and detrainment are developed from LES of the Southern Ocean boundary layer when the turbulence, stratification, and shear cannot be assumed to be in equilibrium with diurnal variability in surface flux and wave (Stokes drift) forcing. A major consequence is the failure of downgradient eddy viscosity, which becomes more serious with Stokes drift and is overcome by relating the angle between the stress and shear vectors to the orientations of Lagrangian shear to the surface and of local Eulerian shear over 5 m. Thus, the momentum flux can be parameterized as a stress magnitude and this empirical direction. In addition, the response of a deep boundary layer to sufficiently strong diurnal heating includes boundary layer collapse and the subsequent growth of a morning boundary layer, whose depth is empirically related to the time history of the forcing, as are both morning detrainment and afternoon entrainment into weak diurnal stratification. Below the boundary layer, detrainment rules give the maximum buoyancy flux and its depth, as well a specific stress direction. Another rule relates both afternoon and nighttime entrainment depth and buoyancy flux to surface layer turbulent kinetic energy production integrals. These empirical relationships are combined with rules for boundary layer transport to formulate two parameterizations; one based on eddy diffusivity and viscosity profiles and another on flux profiles of buoyancy and of stress magnitude. Evaluations against LES fluxes show the flux profiles to be more representative of the diurnal cycle, especially with Stokes drift.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Long-Term Cyclability of NCM-811 at High Voltages in Lithium-Ion Batteries: an In-Depth Diagnostic Study

High-nickel layered oxides, such as LiNi 0.8 Co 0.1 Mn 0.1 O 2 (NCM-811), offer higher energy density than their low-nickel counterparts at a given voltage and are gaining major traction in automotive lithium-ion batteries for electric vehicles. Besides high-Ni content, higher charging voltages above 4.3 V vs Li + /Li boost the energy and represent another focus in battery development. Here, we investigate the long-term cyclability of NCM-811 in graphite pouch cells over 1000 deep cycles between 2.5–4.2, 2.5–4.4, and 2.5–4.5 V through a suite of sensitive characterization techniques. The NCM-811 full cells show severely deteriorated cyclability with higher charging voltages, from 78% at 4.2 V to 52 and 32% at 4.4 and 4.5 V, respectively. At 4.2 V, minor parasitic electrolyte oxidation as well as cathode Li/Ni mixing and cracking are revealed after cycling, while transition-metal dissolution and surface reconstruction into rock-salt NiO are virtually undetectable. At 4.4 and 4.5 V, transition-metal dissolution and crossover to the anode become much more pronounced and a primary contributor to the capacity fade, while significant surface NiO formation causes substantial voltage polarization, which is less noticeable at 4.2 V. Meanwhile, more severe electrolyte oxidation, Li/Ni mixing, and particle pulverization exacerbate the voltage and capacity fade. Finally, these results outline distinct challenges for stable high-Ni layered oxide cathodes in high-voltage Li-ion batteries.

25 ENERGY STORAGE↗