Search NASA⌕ Search

SEARCH · Search NASA

Results for “fractional diffusion”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A Public-Private-Academic Partnership to Advance Solar Power Forecasting

Executive Summary NOAA is making major contributions to the solar forecasting project in three areas. First, it is improving its forecasts of solar irradiance, clouds, and aerosols in its numerical weather prediction models. Second, it is providing advanced satellite products for DOE's FOA awardees to use in their forecast systems. Third, it is using high-quality ground-based measurements from SURFRAD and ISIS stations to verify and validate forecast model output. This reports covers results from all three areas for the period May 1, 2014 - April 30, 2015. Modeling In its modeling effort, NOAA continues work to improve the skill of solar forecasts from the Earth System Research Lab (ESRL) research versions of the 13-km Rapid Refresh (RAP) and the 3-km High-Resolution Rapid Refresh (HRRR) models, which are in turn transitioned into operations at the National Centers for Environmental Prediction (NCEP). A major milestone was achieved in September 2014 with the initial operational implementation of the HRRR at NCEP. In the ESRL research versions of the models, testing and development, in both real-time runs and retrospective experiments, is guided by an extensive in-house verification system. Early in the SFIP project, we developed the capability to verify our model forecasts against the high-quality surface radiation measurements from the SURFRAD and ISIS networks. This highlighted some shortcomings with the RAP and HRRR forecasts of incoming shortwave radiation. Most of our effort during Phase 1 of SFIP was focused on addressing these problems with a variety of model system improvements. The RAP and HRRR models during the warm season of 2014 had a noticeable warm and dry bias in near-surface conditions over most of the central and eastern United States, and our new SURFRAD/ISIS verification revealed that there was also a large excess of incoming global horizontal irradiance in the models. We hypothesized that a lack of cloud cover (particularly low-level cloud cover) in the models was resulting in too much heating of the land surface. This, in turn, caused unrealistically strong surface heat fluxes and turbulent mixing in the planetary boundary layer (PBL), which further reduced the already deficient cloud cover. We addressed these issues with a combination of data assimilation system modifications and model physics improvements. Many of our data assimilation changes were made with a view towards improving the near-term representation of clouds and precipitation. One of these changes involved better accounting for regions of weak reflectivity in the RAP cloud / hydrometeor assimilation system, in order to improve the representation of light precipitation in the RAP initial conditions and provide more realistic initial cloud cover. Additional modifications more accurately accounted for radar beam blockage and data gaps (particularly in the western United States), which improves shorter lead times forecasts of clouds and precipitation. We have also tested the assimilation of new data sources within the RAP and the HRRR, including radar radial velocity data and surface mesonet observations. Within the HRRR, we have tested the cycling of the 3-km land surface fields to allow a higher-resolution treatment of land surface processes. In terms of model physics development for SFIP, we have implemented a shallow cumulus scheme within the RAP, and have made numerous improvements to the Mellor-Yamada-Nakanishi-Niino (MYNN) PBL scheme to address insufficient low-level cloud cover in the models. We have conducted tests incorporating the radiation effects of (parameterized) boundary-layer clouds within the modified MYNN PBL scheme (independent of the convective schemes). The Grell-Freitas-Olson shallow cumulus scheme has also been tested within the 3-km HRRR. Finally, we have also modified the RUC land surface model (LSM) treatment of the vegetation wilting point, reducing it to increase evapotranspiration and increase cloud cover in the boundary layer. All of these changes work in tandem to significantly improve the model forecasts of cloud cover, incoming shortwave radiation, and near-surface temperature and moisture. Satellite The role of NOAA/NESDIS in the Solar Forecasting Improvement Project is to provide Advanced Satellite Products (ASPs) for the two forecasting teams at NCAR and IBM. The ASPs are cloud, surface, and atmosphere products derived from geostationary satellite imagery at the highest possible spatial and temporal resolution - such quantities as cloud mask, cloud probability, cloud transmission, cloud top height, cloud top temperature, cloud effective particle size, etc. Ancillary data, such as elevation and numerical weather prediction fields are provided in the files at the same resolution as well. There are at this time 147 different variables in the ASP output, including quality flags and processing information. The main goals for Year 1 of the project were to implement an Advanced Satellite Products system for the use of the IBM and NCAR teams, begin validation, and make any needed changes based on feedback from the teams. ASP files are being produced every GOES Imager acquisition, which occur on a 15-30 minute schedule. Processing is done on a dedicated computer, with a turn-around time of 8-21 minutes from image acquisition to results available on ftp. Several helpful visualizations of the data are also created for users on web pages. Users have been provided with a document titled "User's Guide for 1km Cloud Products Derived from GOES Imager Data using CLAVR-x", which discusses the basics of the source imagery, the process by which it is turned into Advanced Satellite Products, and considerations users should make when using the data. Validation of selected variables from the older 4km version of the products was also included. Future work will concentrate on validation of the 1km products and improving the turn-around time, product variety, and product quality as needed. Ground Observations In the ground-based measurement effort, NOAA's main objectives are to provide high quality radiation products for validation and verification of short-term to day-ahead solar forecasts. More specifically for the three year project, our goals include (1) Maintaining and providing data from our 7 SURFRAD and 7 ISIS; (2) Update ISIS radiation measurements from 3 min to 1 min data: (3) Purchase and install new pyrheliometers for direct solar irradiance measurements at the 7 SURFRAD sites; (4) Building, testing, and deploying two mobile SURFRAD stations at two utility plants in collaboration with DOE sponsored partners, and includes ongoing maintenance and processing of the data at the mobile sites; (5) Upgrading the data acquisition and communications at 7 SURFRAD sites and 7 ISIS sites; (6) Providing radiation data at the 7 SURFRAD sites in near real-time; (7) Develop and provide aerosol optical depth and cloud images and cloud fraction at our two mobile sites; (8) Provide data recovery rates each year; (9) Provide temporally and spatially averaged radiation products for comparison to HRRR and RAP solar forecasts and advanced satellite products; (10) Provide a data-set for analysis of conversion of direct and diffuse to sloped surfaces; (11) and as time permits develop and provide spectral solar irradiance, cloud optical depth and spectral albedo from the mobile sites. Milestones this year include working with the DOE sponsored teams to find locations to deploy two mobile SURFRAD stations. One existing unit was deployed at a 30MW PV facility in the San Luis Valley in collaboration with Xcel and the NCAR team in August, 2014. The second unit was built and tested at our facilities in Boulder, CO and deployed near Green Mountain Power's Education Center in Rutland, VT in collaboration with Green Mountain Power and the IBM Team in October, 2014. Data processing was implemented and the radiation data from these two mobile sites have been made available on our ftp server in near real-time. We also are providing images and cloud fraction from the TSI cameras for these two mobile sites on our ftp site. Another milestone was upgrading our data acquisition and communication systems at 7 SURFRAD and 7 ISIS sites. We accelerated our schedule for these upgrades to provide timely radiation products. These upgrades allow more reliable and near-real time radiation data delivery to the DOE sponsored teams to meet their goals. Lastly, we changed the data rate at the ISIS sites from 3 min to 1 min.

14 SOLAR ENERGY↗
DOE OSTI2025

Behaviour of cadmium isotopes in sulfidic waters and sediments of the Black Sea: Implications for global cadmium cycling and the application of cadmium isotopes as a paleo-oceanographic proxy

The cadmium isotope system has found use as a tracer for biological productivity, redox and organic carbon burial in the oceans. There are, however, very few observational constraints on Cd isotopic behaviour in modern sulfidic marine conditions, limiting our understanding of the modern Cd cycle, and our ability to use Cd isotope measurements of ancient sedimentary deposits as a paleoceanographic proxy. Here we study the behaviour of dissolved Cd in the water column of the Black Sea and its incorporation into sediments on the basin floor. The isotopic composition of dissolved Cd in the upper ∼50 m of the Black Sea water column is controlled by a combination of biological uptake and regeneration along with mixing of river water and Mediterranean seawater. Cadmium declines to <2 % of its peak subsurface concentration within the zone of nitrate-reduction but above the sulfide chemocline. The isotopic composition of dissolved Cd below the sulfide chemocline evolves in a manner that is consistent with experimental Cd sulfide fractionation factors, providing field-based support for these earlier studies. In contrast, the stability of the dissolved Cd isotopic signature above the chemocline, despite a large reduction in dissolved Cd concentrations, is explained by diffusion towards the Cd-sulfide sink below the chemocline. Sediments accumulating in the deep Black Sea below the sulfide chemocline have isotopic compositions that are similar to dissolved Cd in its near-surface aqueous concentration maximum. Sediments accumulating at shallower depths have lighter isotopic compositions that are likely affected by non-sulfidic burial pathways mediated by macrofaunal mixing and diagenesis of shallow burial phases. Our study shows that the isotopic composition of bulk sediments accumulating under sulfidic conditions is similar to dissolved Cd in the upper water column, supporting its use as a paleo-chemical proxy for ancient seawater as long as basin-scale controls on open ocean seawater sources are considered. We show that Cd burial in the deep Black Sea, predominantly as CdS, accounts for ∼0.5–15 % of total annual Cd burial in the global ocean.

Biogeochemical cycles↗
DOE OSTI2026

Ion Transport and Crystal Rotation in Plastic Crystal Electrolytes Under Applied Electric Fields

Organic ionic plastic crystal electrolytes, containing a plastic crystal and lithium salt, offer a potential balance between mechanical and electrochemical properties for solid state lithium-ion battery electrolytes. These electrolytes contain multiple mobile ionic species (three or four), resulting in complex transport mechanisms which have not yet been established. Plastic crystals are defined by long-range positional order and short-range rotational disorder. It is therefore necessary to quantify changes in the local crystal structure of the electrolyte as current flows through it. Herein, we examine the electrochemical properties of pyrrolidinium-based plastic crystal electrolytes containing lithium salt and zwitterion additives, including measurements of current fraction and limiting current. We obtain species-specific insight into electrolyte transport using pulsed-field gradient nuclear magnetic resonance spectroscopy and find that, while the zwitterion additive increases ionic conductivity, it decreases lithium diffusivity with respect to other ionic components. With operando spatiotemporally resolved wide-angle X-ray scattering we observe location-specific crystal rotations due to the passage of ionic current. In conclusion, we posit that reducing energy dissipation due to rotation is essential for using plastic crystal electrolytes in practical applications.

Yap, Kyra M. K. [University of California, Berkele↗
DOE OSTI2023

Morphological analysis of the polarized synchrotron emission with WMAP and Planck

The bright polarized synchrotron emission, away from the Galactic plane, originates mostly from filamentary structures. We implement a filament finder algorithm which allows the detection of bright elongated structures in polarized intensity maps. We analyse the sky at 23 and 30 GHz as observed respectively by WMAP and Planck. We identify 19 filaments, 13 of which have been previously observed. For each filament, we study the polarization fraction, finding values typically larger than for the areas outside the filaments, excluding the Galactic plane, and a fraction of about 30% is reached in two filaments. We study the polarization spectral indices of the filaments, and find a spectral index consistent with the values found in previous analysis (about -3.1) for more diffuse regions. Decomposing the polarization signals into the E and B families, we find that most of the filaments are detected in P E , but not in P B . We then focus on understanding the statistical properties of the diffuse regions of the synchrotron emission at 23 GHz. Using Minkowski functionals and tensors, we analyse the non-Gaussianity and statistical isotropy of the polarized intensity maps. For a sky coverage corresponding to 80% of the fainter emission, and on scales smaller than 6 degrees (ℓ > 30), the deviations from Gaussianity and isotropy are significantly higher than 3σ. The level of deviation decreases for smaller scales, however, it remains significantly high for the lowest analised scale (~ 1.5°). When 60% sky coverage is analysed, we find that the deviations never exceed 3σ. Finally, we present a simple data-driven model to generate non-Gaussian and anisotropic simulations of the synchrotron polarized emission. The simulations are fitted in order to match the spectral and statistical properties of the faintest 80% sky coverage of the data maps.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Tomographic Dynamics and Scale-Dependent Viscosity in 2D Electron Systems

Fermi gases in two dimensions display collective dynamics originating from head-on collisions, a collinear carrier scattering process that dominates angular relaxation at not-too-high temperatures T<< TF. In this regime, a large family of excitations emerges, with an odd-parity angular structure of momentum distribution and exceptionally long lifetimes. This leads to “tomographic” dynamics: fast 1D spatial diffusion along the unchanging velocity direction accompanied by a slow angular dynamics that gradually randomizes velocity orientation. Here, the tomographic regime features an unusual hierarchy of timescales and scale-dependent transport coefficients with nontrivial fractional scaling dimensions, leading to fractional-power current flow profiles and unusual conductance scaling versus sample width.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1977

A survey of interstellar molecular hydrogen. I

Data from the Copernicus satellite's ultraviolet telescope were used to survey column densities of atomic and molecular hydrogen from a large sample of early-type stars; these data have bearing on an eventual understanding of diffuse and dense interstellar clouds. Column densities are derived by fitting damping profiles to the observed spectra, most of which exhibit strong damping lines in the lower rotational levels surveyed. Plots of dust column density, fractional abundance of molecular hydrogen, and the logarithm of fractional abundance versus total gas column density are given for many of the stars; stars with abnormally large or small hydrogen column densities, as well as some distant stars at high galactic latitudes, are considered. Equilibrium and nonequilibrium theories accounting for the abundance of interstellar hydrogen are compared, and support is found in the data for an account which balances hydrogen formation on interstellar grains with destruction through photodissociation. Overall averages for atomic and molecular hydrogen levels in the galactic plane are also calculated.

Savage, B. D.↗
NASA NTRS1979

Spatial distribution and broad-band spectral characteristics of the diffuse X-ray background, 0.1 - 1.0 keV

Preliminary maps covering more than 85 percent of the sky are presented for three energy bands: the B band, the C band, and the M band. The study was undertaken to find evidence that most of the diffuse X-ray background at energies less than 1 keV is local to the galaxy and that it is most probably due to thermal radiation from a low density plasma which fills a substantial fraction of interstellar space. A preliminary analysis of the data is provided including a report that most of the B and C band flux has a common origin, probably in a 10 to the 6th power K region surrounding the Sun, and that most of the M band flux does not originate from the same material.

Mccammon, D.↗
NASA NTRS1983

Possible contributions of supernova remnants to the soft X-ray diffuse background (0.1 - 1 keV)

Almost all of the B band (0. 10-0.19 keV) and C band (0.15-0.28 keV) X-rays probably originate in a hot region surrounding the sun, which Cox and Anderson modeled as a supernova remnant. This same region may account for a significant fraction of the M band (0.5-1 keV) X-rays if the nonequilibrium models of Cox and Anderson are applicable. A population of distant SNR similar to the local region, with center-to-center spacing of about 300 pc, could provide enough galactic M band emission to fill in the dip in the count rate in the galactic plane that would otherwise be present due to absorption of both the extra galactic power law flux and any large-scale-height stellar (or galactic halo) emission.

Sanders, W. T.↗
NASA NTRS1985

Buoyancy effects on smoldering combustion

The effect of buoyancy on the rate of spread of a concurrent smolder reaction through a porous combustible material is investigated theoretically and experimentally. In the experiments, buoyant forces are controlled by varying the density difference, and the smolder rate spread through porous alpha cellulose (0.83 void fraction) is measured as a function of the ambient air pressure. The smolder velocity is found to increase with the ambient pressure; extinction occurs when the buoyancy forces cannot overcome the drag forces, indicating that diffusion by itself cannot support the spread of a smolder reaction. Theoretical predictions are found to be in good qualitative agreement with the experimental results.

Dosanjh, S.↗
NASA NTRS1994

Microgravity foam structure and rheology

Our long-range objective is to establish the fundamental interrelationship between the microscopic structure and dynamics of foams and their macroscopic stability and rheology. Foam structure and dynamics are to be measured directly and noninvasively through the use and development of novel multiple light scattering techniques such as diffusing-wave spectroscopy (DWS). Foam rheology is to be measured in a custom rheometer which allows simultaneous optical access for multiple light drainage of liquid from in between gas bubbles as the liquid:gas volume fraction in increased towards the rigidity-loss transition.

Durian, Douglas J.↗
NASA NTRS2001

Characteristics of Non-Premixed Turbulent Flames in Microgravity

This project is concerned with the characteristics of turbulent hydrocarbon (primarily propane) gas-jet diffusion flames in microgravity. A microgravity environment provides the opportunity to study the structure of turbulent diffusion flames under momentum-dominated conditions (large Froude number) at moderate Reynolds number which is a combination not achievable in normal gravity. This paper summarizes progress made since the last workshop. Primarily, the features of flame radiation from microgravity turbulent jet diffusion flames in a reduced gravity environment are described. Tests were conducted for non-premixed, nitrogen diluted propane flames burning in quiescent air in the NASA Glenn 5.18 Second Zero Gravity Facility. Measured flame radiation from wedge-shaped, axial slices of the flame are compared for microgravity and normal gravity flames. Results from numerical computations of the flame using a k-e model for the turbulence are also presented to show the effects of flame radiation on the thermal field. Flame radiation is an important quantity that is impacted by buoyancy as has been shown in previous studies by the authors and also by Urban et al. It was found that jet diffusion flames burning under microgravity conditions have significantly higher radiative loss (about five to seven times higher) compared to their normal gravity counterparts because of larger flame size in microgravity and larger convective heat loss fraction from the flame in normal gravity. These studies, however, were confined to laminar flames. For the case of turbulent flames, the flame radiation is a function of time and both the time-averaged and time-dependent components are of interest. In this paper, attention is focused primarily on the time-averaged level of the radiation but the turbulent structure of the flame is also assessed from considerations of the radiation power spectra.

Hegde, U.↗
DOE OSTI2014

WE-G-BRD-01: Diffusion Weighted MRI for Response Assessment of Inoperable Lung Tumors for Patients Undergoing SBRT Treatment

Purpose: To investigate early changes in tumor Apparent Diffusion Coefficients derived from diffusion weighted (DW)-MRI of lung cancer patients undergoing SBRT, as a possible early predictor of treatment response. Methods: DW-MRI scans were performed in this prospective phase I IRB-approved study of inoperable lung tumors at various time-points during the course of SBRT treatments. Axial DW scan using multi b-values ranging from 0–1000 s/mm{sup 2} were acquired in treatment position on a 3T Philips MR scanner during simulation, one hour after the first fraction (8 Gy), after a total of 5 fractions (40 Gy) and 4 weeks after SBRT delivery. A monoexponential model based on a least square fit from all b values was performed on a pixel-by-pixel basis and ADC was calculated. GTVs drawn on 4DCT for planning were mapped on the T2w MRI (acquired at exhale) after deformable registration. These volumes were then mapped on DWI scan for ADC calculation after rigid registration between the anatomical scan and diffusion scan. T2w scan on followup time points were deformably registered to the pretreatment T2 scan. Results: The first two patients in this study were analyzed. Median ADC values were 1.48, 1.48, 1.62 and 1.83 (10{sup −3}×) mm{sup 2}/s at pretreatment, after 8 Gy, after 40 Gy and 4 weeks posttreatment for the first patient and 1.57, 1.53, 1.66 and 1.72 (10{sup −3}×) mm{sup 2}/s for the second patient. ADC increased more significantly after 4 weeks of treatment rather than immediately post treatment, implying that late ADC value may be a better predictor of tumor response for SBRT treatment. The fraction of tumor pixels at high ADC values increased at 4 weeks post treatment. Conclusion: The observed increase in ADC values before the end of radiotherapy may be a surrogate for tumor response, but further patient accrual will be necessary to determine its value.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2014

Adaptive two-regime method: Application to front propagation

The Adaptive Two-Regime Method (ATRM) is developed for hybrid (multiscale) stochastic simulation of reaction-diffusion problems. It efficiently couples detailed Brownian dynamics simulations with coarser lattice-based models. The ATRM is a generalization of the previously developed Two-Regime Method [Flegg et al., J. R. Soc., Interface 9, 859 (2012)] to multiscale problems which require a dynamic selection of regions where detailed Brownian dynamics simulation is used. Typical applications include a front propagation or spatio-temporal oscillations. In this paper, the ATRM is used for an in-depth study of front propagation in a stochastic reaction-diffusion system which has its mean-field model given in terms of the Fisher equation [R. Fisher, Ann. Eugen. 7, 355 (1937)]. It exhibits a travelling reaction front which is sensitive to stochastic fluctuations at the leading edge of the wavefront. Previous studies into stochastic effects on the Fisher wave propagation speed have focused on lattice-based models, but there has been limited progress using off-lattice (Brownian dynamics) models, which suffer due to their high computational cost, particularly at the high molecular numbers that are necessary to approach the Fisher mean-field model. By modelling only the wavefront itself with the off-lattice model, it is shown that the ATRM leads to the same Fisher wave results as purely off-lattice models, but at a fraction of the computational cost. The error analysis of the ATRM is also presented for a morphogen gradient model.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2014

Solar Wind Implantation into Lunar Regolith: Hydrogen Retention in a Surface with Defects

Solar wind protons are implanted directly into the top 100 nm of the lunar near-surface region, but can either quickly diffuse out of the surface or be retained, depending upon surface temperature and the activation energy, U, associated with the implantation site. In this work, we explore the distribution of activation energies upon implantation and the associated hydrogen-retention times; this for comparison with recent observation of OH on the lunar surface. We apply a Monte Carlo approach: for simulated solar wind protons at a given local time, we assume a distribution of U values with a central peak, U(sub c) and width, U(sub w), and derive the fraction retained for long periods in the near-surface. We find that surfaces characterized by a distribution with predominantly large values of U (greater than 1 eV) like that expected at defect sites will retain implanted H (to likely form OH). Surfaces with the distribution predominantly at small values of U (less than 0.2 eV) will quickly diffuse away implanted H. However, surfaces with a large portion of activation energies between 0.3 eV less than U less than 0.9 eV will tend to be H-retentive in cool conditions but transform into H-emissive surfaces when warmed (as when the surface rotates into local noon). These mid-range activation energies give rise to a diurnal effect with diffusive loss of H at noontime.

regolith↗
DOE OSTI2014

Dynamics of Block Copolymer Nanocomposites

A detailed study of the dynamics of cadmium sulfide nanoparticles suspended in polystyrene homopolymer matrices was carried out using X-ray photon correlation spectroscopy for temperatures between 120 and 180 °C. For low molecular weight polystyrene homopolymers, the observed dynamics show a crossover from diffusive to hyper-diffusive behavior with decreasing temperatures. For higher molecular weight polystyrene, the nanoparticle dynamics appear hyper-diffusive at all temperatures studied. The relaxation time and characteristic velocity determined from the measured hyper-diffusive dynamics reveal that the activation energy and underlying forces determined are on the order of 2.14 × 10 -19 J and 87 pN, respectively. We also carried out a detailed X-ray scattering study of the static and dynamic behavior of a styrene– isoprene diblock copolymer melt with a styrene volume fraction of 0.3468. At 115 and 120 °C, we observe splitting of the principal Bragg peak, which we attribute to phase coexistence of hexagonal cylindrical and cubic double- gyroid structure. In the disordered phase, above 130 °C, we have characterized the dynamics of composition fluctuations via X-ray photon correlation spectroscopy. Near the peak of the static structure factor, these fluctuations show stretched-exponential relaxations, characterized by a stretching exponent of about 0.36 for a range of temperatures immediately above the MST. The corresponding characteristic relaxation times vary exponentially with temperature, changing by a factor of 2 for each 2 °C change in temperature. At low wavevectors, the measured relaxations are diffusive with relaxation times that change by a factor of 2 for each 8 °C change in temperature.

36 MATERIALS SCIENCE↗
NASA NTRS1979

Extensions to the classical calculation of the effect of mutual shadowing in diffuse reflection

The classical method for accounting for the mutual shadowing among closely packed particles in multiple scattering calculations is extended in the following ways. (1) By modeling the particle distribution by a Poisson process with a varying density parameter, a 'Van der Waals' type approximation allows extension to a greater fractional volume density, D. In this case it is only required that D squared be much less than 1 instead of D being much less than 1. (2) In the case that the particle distribution is not uniform the classical calculation may be weighted by the pair correlation function of the distribution. (3) The use of the Markov chain formalism for radiative transfer allows inclusion of the effect of shadowing for two orders of scattering. For conditions such as might apply in Saturn's rings, the inclusion of this effect makes less than a 0.1% difference in the calculated phase curves, compared to previous calculations which have included shadowing only in the first scattering. The latter are thus shown to be quite accurate.

Esposito, L. W.↗
NASA NTRS2005

The nakhlite meteorites: Augite-rich igneous rocks from Mars

The seven nakhlite meteorites are augite-rich igneous rocks that formed in flows or shallow intrusions of basaltic magma on Mars. They consist of euhedral to subhedral crystals of augite and olivine (to 1 cm long) in fine-grained mesostases. The augite crystals have homogeneous cores of Mg' = 63% and rims that are normally zoned to iron enrichment. The core-rim zoning is cut by iron-enriched zones along fractures and is replaced locally by ferroan low-Ca pyroxene. The core compositions of the olivines vary inversely with the steepness of their rim zoning - sharp rim zoning goes with the most magnesian cores (Mg' = 42%), homogeneous olivines are the most ferroan. The olivine and augite crystals contain multiphase inclusions representing trapped magma. Among the olivine and augite crystals is mesostasis, composed principally of plagioclase and/or glass, with euhedra of titanomagnetite and many minor minerals. Olivine and mesostasis glass are partially replaced by veinlets and patches of iddingsite, a mixture of smectite clays, iron oxy-hydroxides and carbonate minerals. In the mesostasis are rare patches of a salt alteration assemblage: halite, siderite, and anhydrite/ gypsum. The nakhlites are little shocked, but have been affected chemically and biologically by their residence on Earth. Differences among the chemical compositions of the nakhlites can be ascribed mostly to different proportions of augite, olivine, and mesostasis. Compared to common basalts, they are rich in Ca, strongly depleted in Al, and enriched in magmaphile (incompatible) elements, including the LREE. Nakhlites contain little pre-terrestrial organic matter. Oxygen isotope ratios are not terrestrial, and are different in anhydrous silicates and in iddingsite. The alteration assemblages all have heavy oxygen and heavy carbon, while D/H values are extreme and scattered. Igneous sulfur had a solar-system isotopic ratio, but in most minerals was altered to higher and lower values. High precision analyses show mass-independent fractionations of S isotopes. Nitrogen and noble gases are complex and represent three components: two mantle sources (Chas-E and Chas-S), and fractionated Martian atmosphere. The nakhlites are igneous cumulate rocks, formed from basaltic magma at approx.1.3 Ga, containing excess crystals over what would form from pure magma. After accumulation of their augite and olivine crystals, they were affected (to various degrees) by crystallization of the magma, element diffusion among minerals and magma, chemical reactions among minerals and magma, magma movement among the crystals, and post-igneous chemical equilibration. The extent of these modifications varies, from least to greatest, in the order: MIL03346, NWA817, Y000593, Nakhla = Governador Valadares, Lafayette, and NWA998. Chemical, isotopic, and chronologic data confirm that the nakhlites formed on Mars, most likely in thick lava flows or shallow intrusions. Their crystallization ages, referenced to crater count chronologies for Mars, suggest that the nakhlites formed on the large volcanic constructs of Tharsis, Elysium, or Syrtis Major. The nakhlites were suffused with liquid water, probably at approx.620 ma. This water dissolved olivine and mesostasis glass, and deposited iddingsite and salt minerals in their places. The nakhlites were ejected from Mars at approx.10.75Ma by an asteroid impact and fell to Earth within the last 10,000 years. Although the nakhlites are enriched in incompatible elements, their source mantle was strongly depleted. This depletion event was ancient, as the nakhlites source mantle was fractionated while short-lived radionuclides (e.g., t(sub 1/2 = 9 my) were still active. This differentiation event may have been core formation coupled with a magma ocean, as is inferred for the moon.

Treiman, Allan H.↗
DOE OSTI2020

Gas Transport in the Insensitive High Explosive PBX 9502

Small-scale laboratory experiments were performed to analyze gas transport mechanics in PBX 9502, an insensitive high explosive (IHE) composition. Two independent investigations are reported here. First, gas permeametry and dynamic pycnometry techniques were used to measure the molecular flow (Knudsen) coefficient and the internally interconnected void fraction of pristine and thermally damaged samples in two uniaxial pressing orientations. The permeability of PBX 9502 was found to be unmeasurably low (of the same order of magnitude as PTFE) with gas transport being diffusion-dominated. Secondly, a pressure vessel experiment was developed to measure quasi-static and dynamic gas generation as the explosive was heated to self-ignition (cookoff). The gas generation results and the permeability/diffusivity findings provide evidence that PBX 9502 remains impermeable until seconds prior to self-ignition. At ignition, internal void-pressure drives macro-scale cracking and the sample becomes uniformly incorporated in the following deconsolidating deflagration. These results are discussed within the context of previous observations of pressure-dependent cookoff behavior of PBX 9502 and provide a more complete description of thermal damage evolution in this explosive composition.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗