Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “subsurface flow”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

3D Quantification of Pore Networks and Anthropogenic Carbon Mineralization in Stacked Basalt Reservoirs

Basalt formations are promising candidates for the geologic storage of anthropogenic CO 2 due to their storage capacity, porosity, permeability, and reactive geochemical trapping ability. The Wallula Basalt Carbon Storage Pilot Project demonstrated that supercritical CO 2 injected into >800 m deep Columbia River Basalt Group stacked reservoir flow tops mineralizes to ankerite-siderite-aragonite on month-year timescales, with 60% of the 977 metric tonnes of CO 2 converted within two years. The potential impacts of mineral precipitation and consequent changes on rock porosity, pore structure, pore size, and pore size distributions have likely been underestimated hitherto. Herein, we address these knowledge gaps using X-ray Micro Computed Tomography (XMT) to evaluate the pore network architecture of sidewall cores recovered two years after CO 2 injection. In the present study, we performed a detailed quantitative analysis of the CO 2 -reacted basalt cores by XMT imaging. Reconstructed 3D images were analyzed to study the distribution and volumetric details of porosity and carbonate nodules in the cores along with the various other phases, providing insight into paragenesis and carbonate growth mechanisms, including mineralogic/chemical zonation. Finally, these findings are being used to parametrize multiphase reactive transport models to predict the fate and transport of subsurface CO 2 , enabling scale-up to commercial-scale geologic carbon storage in basalts and other reactive mafic-ultramafic formations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Automatic Time Step Control to Resolve Hydromechanically Driven Fault Reactivation, Spontaneous Nucleation, and Seismic Arrest

Abstract A physical understanding of the progression from flow‐driven (quasi‐static) poromechanical deformation to dynamic fault rupture is critical to the resilient operations of several engineering systems. These processes are bridged by a progression from fault reactivation to the spontaneous nucleation of unstable sliding. Toward addressing this challenge, novel automatic time step size control methods are developed to enable accurate and efficient simulation of these dynamics and transitions from the first principles. The controllers combine local models for discretization error and Coulomb failure conditions to automatically adjust the time step size across several orders of magnitude. The methods do not require additional empirical or theoretical input and can resolve the pre‐rupture, interseismic, and seismic periods to the allowed accuracy. The computational results reveal that the proposed methods automatically capture the onset of reactivation and nucleation for homogeneous and heterogeneous fields. Hydrodynamic and structural heterogeneity lead to disparate critical nucleation sizes compared to those predicted by theory. The results highlight its potential in predicting induced seismicity in realistic subsurface engineering systems and at practical scales.

Environmental Sciences & Ecology↗
DOE OSTI2023

Model simulations of Plum Island Ecosystems LTER low marsh site using ELM-PFLOTRAN

Model simulations using the E3SM Land Model (ELM) coupled to the PFLOTRAN reactive transport model via the Alquimia interface. The simulations were conducted for a tidal salt marsh at the Plum Island Ecosystems LTER near Rowley, Massachusetts, USA. Model simulations were forced using site-specific tidal cycles and salinity, and the simulations used a biogeochemical reaction network including aerobic decomposition, sulfate reduction, iron reduction, and methanogenesis. Model outputs include simulated carbon stocks, carbon dioxide and methane fluxes, and porewater concentrations of key solutes related to sulfur, iron, and carbon cycling. The model simulations included a saline simulation (with tidal sulfate inputs), a fresh simulation (with low salinity and low sulfate inputs), and a saline simulation with lower vegetation productivity to represent the effect of salinity on vegetation. These simulations were conducted to demonstrate that a new model framework incorporating subsurface redox and biogeochemical interactions into a land surface model could reproduce measured surface greenhouse gas fluxes and biogeochemical dynamics in tidal marsh ecosystems, and to test whether including redox interactions in a land surface model would allow the model to resolve contrasts in biogeochemical cycling and greenhouse gas production between saline and freshwater wetlands.The data package includes gzipped tar archives (which can be expanded using standard tar and gzip utilities) of model outputs from three model configurations: saline subsurface and reduced vegetation productivity related to salinity; saline subsurface with vegetation productivity not reduced; and freshwater. Also included are code for the modified E3SM model, Alquimia interface, and PFLOTRAN reactive transport simulator in gzipped tar format; plain text parameter and configuration files; python code files for visualizing model output and defining model configurations; and model output, tide and salinity forcing, and configuration files in netCDF format. See the README.md file in the data package for a detailed description of all files contained in the package. All files are in netCDF (.nc), gzipped tar archive (.tar.gz or .tgz), or text (all other files).Updated: May 13, 2024. Model output, E3SM code, PFLOTRAN input files, and python codes for visualizing results were updated to reflect changes made for the manuscript revision. The updated archive reflects the code and model output from the final accepted manuscript. Changes included updated reaction parameters reflecting improved parameterization and additional comparisons with field measurements. E3SM code changes included better support for multiple grid cells and improved flow and transport parameterization.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Integration of Large-Scale Electrical Imaging into Geological Framework Development and Refinement

Geologic framework models (GFMs) are critical to the construction of reliable simulation models of groundwater flow and contaminant transport. To support GFM development, direct information (e.g., core samples, fluid samples, hydraulic testing) tends to be sparse and separated by large distances relative to the spatial scales of aquifer heterogeneity. There are additional challenges associated with highly contaminated legacy waste sites, where drilling is particularly costly, and invasive sampling requires specialized handling and disposal of hazardous materials. At these sites in particular, non-invasive geophysical imaging can play an important role in filling spatial gaps between boreholes and reducing characterization costs by optimizing and minimizing the number of necessary boreholes. Here this paper presents a case study demonstrating the use of large-scale (> 30 km 2 ) electrical mapping to identify hydrostratigraphy and potential paleochannels at the Hanford Site, located in Washington State, USA. In two field campaigns, over 36 line-kilometers of electrical resistivity tomography (ERT) data were collected along 14 transects. ERT surveys were sited and performed to image critical aspects (e.g., paleochannels, stratigraphic contacts) of the subsurface, demonstrating a general workflow for integrating ERT with GFM development. Inconsistencies between the GFM and ERT were catalogued to provide a basis for future site characterization using complementary geophysical methods and (or) direct sampling.

58 GEOSCIENCES↗
NASA NTRS1999

Unique Aeolian Transport Mechanisms on Mars: Respective Roles of Percussive and Repercussive Grain Populations in the Sediment Load

Experiments show that when sand-size grains impact a sediment surface with energy levels commensurate for Mars, small craters are formed by the ejection of several hundred grains from the bed. The experiments were conducted with a modified crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism. Individual grains of sand could be fired at loose sand targets to observe ballistic effects unhindered by aerodynamic mobilization of the bed. Impact trajectories simulated the saltation process on dune surfaces. Impact craters were not elongated despite glancing (15 deg.) bed impact; the craters were very close to being circular. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 um-diameter grains into similar material. This behavior is explained by deposition of elastic energy in the bed by the "percussive" grain. Impact creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains to open packing and they consequently become forcefully ejected from the site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. A stress model based on repercussive bed dilatancy and interparticle adhesive forces (for smaller grains) predicts, to first order, the observed crater volumes for various impact conditions. On earth, only a few grains are mobilized by a percussive saltating grain; some grains are "knudged" along the ground, and some are partly expelled on short trajectories. These motions constitute reptation transport. On Mars, saltation and reptation become indistinct: secondary or "repercussive" trajectories have sufficient vertical impulse to create a dense saltation population of many tens or hundreds of grains for each single high-speed saltation percussion of the bed. Impact cascading will lead to near-surface distortion of the boundary layer, and choked flow formed by a dense "slurry" of sand, with the majority of grains mobilized by repercussive forces rather than by aerodynamic lift. This proceeds until a fully-matured transport layer imposes self-limitations as grain-population density constrains the free-path motion of individial grains.

Marshall, John R.↗
NASA NTRS1999

Unique Aeolian Transport Mechanisms on Mars: Respective Roles of Percussive and Repercussive Grain Populations in the Sediment Load

Experiments show that when sand-size grains impact a sediment surface with energy levels commensurate for Mars, small craters are formed by the ejection of several hundred grains from the bed. The experiments were conducted with a modified crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism. Individual grains of sand could be fired at loose sand targets to observe ballistic effects unhindered by aerodynamic mobilization of the bed. Impact trajectories simulated the saltation process on dune surfaces. Impact craters were not elongated despite glancing (15 deg.) bed impact; the craters were very close to being circular. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. This behavior is explained by deposition of elastic energy in the bed by the "percussive" grain. Impact creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains to open packing and they consequently become forcefully ejected from the site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. A stress model based on repercussive bed dilatancy and interparticle adhesive forces (for smaller grains) predicts, to first order, the observed crater volumes for various impact conditions. On earth, only a few grains are mobilized by a percussive saltating grain; some grains are "knudged" along the ground, and some are partly expelled on short trajectories. These motions constitute reptation transport. On Mars, saltation and reptation become indistinct: secondary or "repercussive" trajectories have sufficient vertical impulse to create a dense saltation population of many tens or hundreds of grains for each single high-speed saltation percussion of the bed. Impact cascading will lead to near-surface distortion of the boundary layer, and choked flow formed by a dense "slurry" of sand, with the majority of grains mobilized by repercussive forces rather than by aerodynamic lift. This proceeds until a fully-matured transport layer imposes self- limitations as grain-population density constrains the free-path motion of individual grains.

Marshall, John R.↗
NASA NTRS2012

Astrobiological Journeys to and from the South Polar Sea of Enceladus - Bidirectional Interactions with the Saturn Magnetosphere

The spectacularly cryovolcanic moon Enceladus is a major source of plasma for the Saturn magnetosphere via ionization of the ejected molecular species and ice grains. Field-aligned plasma flows from the Enceladus environment visibly impact the moon's magnetic footpoint in Saturn's polar auroral region, while water group and other ions from the moon emissions diffuse radially throughout the magnetosphere and may be the dominant source of oxygen for Titan's oxygen-poor upper atmosphere. But the moon-magnetosphere interaction is bidirectional in the sense that the moon surface is globally exposed to constant irradiation by the returning magnetospheric ions and by energetic electrons from the field-aligned and radially diffusing populations. The returning ion source operates both on global scales of the magnetosphere and locally for highly reactive species produced in the ejecta plume. All of these sources likely combine to produce a highly oxidized global surface layer. Since plasma electrons and ions are cooled by interaction with neutral gas and E-ring ice grains from Enceladus, the moon emissions have a governing effect on the seed populations of energetic particles that irradiate the surface. The proposed subsurface polar sea and transient crustal overturn in the south polar region could bring the polar surface oxidants into contact with hydrocarbons and ammonia to make oxidation product gases contributing to the cryovolcanic jets, a process first proposed by Cooper et al. (Plan. Sp. Sci., 2009). As has been previously suggested for Europa, the 'oxidants could contribute to enhanced astrobiological potential of Enceladus, perhaps even higher than for Europa where organic hydrocarbons have not yet been directly detected. Unlike Europa, Enceladus shows no sign of an oxygen-dominated exosphere that could otherwise be indicative of extreme surface and interior oxidation inhibiting the detectable survival and evolution of complex organics.

Cooper, John F.↗
NASA NTRS2016

Cryovolcanic Emplacement of Domes on Europa

Here we explore the hypothesis that certain domes on Europa may have been produced by the extrusion of viscous cryolavas. A new mathematical method for the emplacement and relaxation of viscous lava domes is presented and applied to putative cryovolcanic domes on Europa. A similarity solution approach is applied to the governing equation for fluid flow in a cylindrical geometry, and dome relaxation is explored assuming a volume of cryolava has been rapidly emplaced onto the surface. Nonphysical sin- gularities inherent in previous models for dome relaxation have been eliminated, and cryolava cooling is represented by a time-variable viscosity. We find that at the onset of relaxation, bulk kinematic viscosities may lie in the range between 10(exp 3) and 10(exp 6) sq m/s, while the actual fluid lava viscosity may be much lower. Plausible relaxation times to form the domes, which are linked to bulk cryolava rheology, are found to range from 3.6 days to 7.5 years. We find that cooling of the cryolava, while dominated by conduction through an icy skin, should not prevent fluids from advancing and relaxing to form domes within the timescales considered. Determining the range of emplacement conditions for putative cryolava domes will shed light on Europa's resurfacing history. In addition, the rheologies and compositions of erupted cryolavas have implications for subsurface cryomagma ascent and local surface stress conditions on Europa.

volcanism↗
DOE OSTI2023

Hawai‘i Play Fairway (Final Report)

Most of Hawai'i's geothermal resources are blind—their manifestations, such as hot springs and steam vents, do not appear on the ground surface because the heated water flows far below. With the exception of $K\bar{i}lauea East Rift Zone$, in most areas of Hawai'i, high lateral permeability in the first kilometer below ground surface prevents surface thermal features from developing. As a methodology for discovering these blind resources, Play Fairway Analysis (PFA) involves finding potential locations of blind hydrothermal systems and describing potential geothermal sources in rift-zone settings. Using the PFA to find Hawai'i's geothermal resources, the University of Hawai'i (UH) conducted the Hawai‘i Play Fairway Project, Hawai'i's first statewide geothermal resource assessment since 1985. Sponsored by the U.S. Department of Energy, the Hawai'i Play Fairway Project provided an updated resource assessment, a roadmap for additional exploration activities, and the identification of areas for further exploration. Benefitting from UH's core competency in earth sciences and experienced geothermal researchers, the project comprised three phases. During the first phase, the team identified, compiled, and ranked existing geologic, groundwater, and geophysical datasets relevant to subsurface heat, fluid and permeability. Using a Bayesian statistical approach, the team developed a statistical methodology to integrate these data into a resource probability map. The team evaluated the confidence in the probability value and considered development viability of areas with geothermal resources. With these analyses, the team identified 10 locations in the Hawaiian Islands for exploration activities. For the second phase, the team collected new groundwater data in 10 locations across the state and new geophysical data on $L\bar{a}na‘i, Maui$, and central Hawai'i Island and modeled topographically induced stress to better characterize subsurface permeability. Analyzing the subsurface stresses, the team evaluated the potential for fracture-induced permeability. The team inverted the MT and gravity data to produce 3D models of resistivity and density, respectively, on $L\bar{a}na‘i$, across $Haleakal\bar{a}'s$ SW rift (Maui), and surrounding Mauna Kea (Hawai‘i Island). The team developed and applied a new method for incorporating depth information about resistivity, density, and potential for fracture-induced permeability into the statistical method for computing resource probability in these three focus areas. The team incorporated the new groundwater results with the new geophysical results and the calculations of potential for fracture-induced permeability to produce updated maps of resource probability and confidence. Through combining data from the first and second phases, the team determined locations for further exploration during the third phase. For MT and gravity surveys, the team recommended $Kaua'i's$ $L\bar{i}hu'e$ $Basin$, the east rift of $Maui's$ $Haleakal\bar{a}$ volcano, and the southwest rift of Hawai'i Island's Mauna Loa volcano. The MT and gravity surveys aimed to enable improved confidence in the resource potential in these locations. For drilling deep groundwater well(s), the team recommended Southeast Mauna Kea and $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$. During the third phase, further exploration involved drilling a groundwater well in $L\bar{a}na's$ $P\bar{a}l\bar{a}wai$ $Basin$ and performing more geophysical surveys. We deepened an existing water well proximal to our target area on $L\bar{a}na'i$ due to funding constraints that precluded us from spudding a new well that would exceed 1km depth. Drilling was preceded by a number of substantial elements including: writing an Environmental Assessment and the subsequent legal process, performance of deviation logging, lowering a camera down the well, coordinating site preparation with $P\bar{u}lama$ $L\bar{a}na'i$, shipping the UH-owned rig interisland, procuring supplies, and leading 3 community meetings on $L\bar{a}na'i$. Drilling occurred 24/7 the entire month of June 2019 over which time $L\bar{a}na'i$ $Well$ 10 was deepened from 427 m to 1057 m, with continuous core collected. We measured a roughly linear temperature gradient averaging 42°C/km and a maximum bottom hole temperature of 66°C. This gradient is more than twice the background for Hawai'i and within a range of gradients measured in this depth range for some exploration wells within KERZ. We consider these results encouraging for $L\bar{a}na'i's$ resource potential and recommend following with a slim hole within $L\bar{a}na'i's$ caldera (our target zone) to ~ 2 km. Further, the positive implications such results have for the island of O‘ahu are substantial - the shield stage of O'ahu's volcanoes ended 1-2 My earlier. However, O'ahu uses more electricity than the rest of the islands combined, and the utility recently called for 500-700MW of firm, dispatchable renewable electricity on O'ahu by 2033. In Phase 3, we also collected limited new encouraging groundwater data, and updated our thoughts on the probabilities of fluid and permeability at resource depths (Pr F = 1; Pr P = mostly unconstrained). Ultimately, we advocate for using our final probability of heat, and confidence in this probability, to drive the next phase of exploration. We contend further development of geothermal in Hawai‘i will enable the state to achieve its 100% renewable policy objective and Hawai'i to transition off of fossil fuels through geothermal discovery and development. The project not only produced a large amount of data and expanded the existing knowledge of Hawai'i's geothermal resources, but also produced publications, theses, presentations, core photos, datasets, media reports, television interviews, community events, and a blog. Students and new professionals benefitted from the project's hands-on research experiences and educational opportunities and earned awards and recognition.

15 GEOTHERMAL ENERGY↗
NASA NTRS2016

Megafans-Some New Perspectives from a Global Study

A global study of megafans (greater than 100 km long) has revealed their widespread existence on all continents, with almost 200 documented, 93 in Africa where research is most thorough. The largest measures 705 km. Megafans are a major subset of "DFS" (distributive fluvial systems, a category that includes all fan-like features greater than 30 km long). 1. Many researchers now recognize megafans as different from floodplains, small coarse-grained alluvial fans, and deltas. Although smaller architectural elements in megafans are the same as those encountered in floodplains (channel, overbank, etc.), larger architectures differ because of the unconfined setting of megafans, versus the valley-confined setting of floodplains. 2. A length continuum is now documented between steep alluvial fans 10-20 km in length, and fluvial fans 30-50 km long. This implies a continuum of process from end-member alluvial fan processes (e.g. high-energy flows that emplace gravels, debris-flow units) to the relatively fine-grained channel and overbank deposits common to purely fluvial fans. Combinations of these different processes will then occur in many mid-sized fans. 3. The global distribution suggests a prima facie relationship with tectonic environment rather than climatic zones, with local controls being the slope of the formative river and the existence of a basin subsiding below the long profile of the river. But the global population has revealed that most megafans are relict. So it is possible that further research will show relationships to prior climatic regimes. 4. Megafans can have regional importance: e.g., along the east flank of the central Andes, nested megafans total approximately 750,000 km2-and 1.2m km2 if all megafans in S. America are counted. Modern megafan landscapes thus have basinal importance, orders of magnitude greater than alluvial fan bajadas. 5. Because so many aggrading basins are dominated today by DFS, it is claimed that DFS ought to be significant in the subsurface; and that existing fluvial models therefore may not apply to the majority of fluvial sedimentary units. Arguments have been raised against this view, but as modern megafan systems become better known they are rapidly being applied as a model in many fluvial basins. A small literature has arisen with apparent examples from every part of the world.

Wilkinson, M. Justin↗
NASA NTRS1985

Volcano/ground ice interactions in Elysium Planitia, Mars

The occurrence within Elysium Planitia of meltwater deposits, possible pseudocraters, collapse features within troughs, and outflow channels indicates that a layer of subsurface volatiles existed at the time of volcanic activity within this area. The pseudocraters are interpreted to be indicators of near-surface volatiles, while meltwater deposits and the degree of preservation of trough walls and floors are thought to signify greater volatile depths. A latitudinal variation in the distribution of these features indicates either that the depth to the volatile layer increased from less than about 50 m at 35 deg N to greater than 600 m at 24 deg N, or that an ice wedge that existed at 35 deg N thinned to nonexistence at 24 deg N. Braided distributary channel systems within the chaotic terrain north of Elysium Planitia show that ephemeral lakes were repeatedly created and drained at this locality. The existence of volatiles contemporaneous with volcanic activity permits a search to be made for explosively generated landforms predicted to exist by previous theoretical models. Morphological evidence for strombolian, vulcanian and plinian eruptions is lacking within western Elysium Planitia; there are no identifiable cinder cones, pyroclastic flow deposits, or mantled areas indicative of large airfall deposits at an image resolution of 50-150 m/pixel. However, the pseudocraters indicate that small-scale phreatomagmatic activity may have taken place.

Mouginis-Mark, P. J.↗
NASA NTRS1999

Aeolian Sand Transport in the Planetary Context: Respective Roles of Aerodynamic and Bed-Dilatancy Thresholds

The traditional view of aeolian sand transport generally estimates flux from the perspective of aerodynamic forces creating the airborne grain population, although it has been recognized that "reptation" causes a significant part of the total airborne flux; reptation involves both ballistic injection of grains into the air stream by the impact of saltating grains as well as the "nudging" of surface grains into a creeping motion. Whilst aerodynamic forces may initiate sand motion, it is proposed here that within a fully-matured grain cloud, flux is actually governed by two thresholds: an aerodynamic threshold, and a bed-dilatancy threshold. It is the latter which controls the reptation population, and its significance increases proportionally with transport energy. Because we only have experience with terrestrial sand transport, extrapolations of aeolian theory to Mars and Venus have adjusted only the aerodynamic factor, taking gravitational forces and atmospheric density as the prime variables in the aerodynamic equations, but neglecting reptation. The basis for our perspective on the importance of reptation and bed dilatancy is a set of experiments that were designed to simulate sand transport across the surface of a martian dune. Using a modified sporting crossbow in which a sand-impelling sabot replaced the bolt-firing mechanism, individual grains of sand were fired at loose sand targets with glancing angles typical of saltation impact; grains were projected at about 80 m/s to simulate velocities commensurate with those predicted for extreme martian aeolian conditions. The sabot impelling method permitted study of individual impacts without the masking effect of bed mobilization encountered in wind-tunnel studies. At these martian impact velocities, grains produced small craters formed by the ejection of several hundred grains from the bed. Unexpectedly, the craters were not elongated, despite glancing impact; the craters were very close to circular in planform. High-speed photography showed them to grow in both diameter and depth after the impactor had ricochetted from the crater site. The delayed response of the bed was "explosive" in nature, and created a miniature ejecta curtain spreading upward and outward for many centimeters for impact of 100-300 micron-diameter grains into similar material. Elastic energy deposited in the bed by the impacting grain creates a subsurface stress regime or "quasi-Boussinesq" compression field. Elastic recovery of the bed occurs by dilatancy; shear stresses suddenly convert the grains from closed to open packing, and grains are consequently able to eject themselves forcefully from the impact site. Random jostling of the grains causes radial homogenization of stress vectors and a resulting circular crater. There is a great temptation to draw parallels with cratering produced by meteorite impacts, but a rigorous search for common modelling ground between the two phenomena has not been conducted at this time. For every impact of an aerodynamically energized grain, there are several hundred grains ejected into the wind for the high-energy transport that might occur on Mars. Many of these grains will themselves become subject to the boundary layer's aerodynamic lift forces (their motion will not immediately die and add to the creep population), and these grains will become indistinguishable from those lifted entirely by aerodynamic forces. As each grain impacts the bed, it will eject even more grains into the flow. A cascading effect will take place, but because it must be finite in its growth, damping will occur as the number of grains set in motion causes mid-air collisions that prevent much of the impact energy from reaching the surface of the bed -thus creating a dynamic equilibrium in a high-density saltation cloud. It is apparent that for a given impact energy, the stress field permits a smaller volume of grains to convert to open packing as the size of the bed grains increases, or as the energy of the "percussive" grain decreases (by decrease in velocity or mass). Thus, the mass of the "repercussive" grain population that is ejected from the impact site becomes a function of the scale of the stress field in relation to the scale of the bed material (self-similarity being applicable if both bed size and energy are simultaneously adjusted). In other words, in a very high energy aeolian system where an aerodynamically raised grain can ballistically raise many more grains, the amount of material lifted into the wind becomes largely a function of a dilatancy threshold. If this threshold is exceeded, grains are repercussively injected into the saltation cloud. The "dilatancy threshold" may be defined in terms of the saltation percussive force required to convert the bed, through elastic response, from a closed to an open packing system. If open packing cannot be created, the grains cannot escape from the impact site, even though the elastic deformation and percussive force may be able to reorganize the grains with respect to one another. As the crossbow experiments showed, for an ever-increasing bed grain size, a point is reached when no material can be moved because the energy of the percussive grain is insufficient to dilate the relatively coarse bed. Although this seems to be stating the obvious -- that too little energy will not cause the bed to splash -- the consequences of exceeding the "splash threshold" by dilatancy are not so obvious for high-energy aeolian transport. It is noted that the force required to elastically dilate the bed has to overcome Coulombic grain attractions such as dipole-dipole coupling, dielectric, monopole, contact-induced dipole attractions, van der Waals forces, molecular monolayer capillary forces, as well as the mechanical interlocking frictional resistance of the grains. On Mars, it is predicted that the dilatancy threshold may be the prime control of grain flux. Additional information is contained in the original.

Marshall, J. R.↗
DOE OSTI2025

Autonomous Infrared and Small (Wide) Angle X-Ray Scattering (IR-S(W)AXS) Capability

Thin water films are 2-D, nanoconfined layers that form on solid surfaces exposed to humid atmospheres—environments ubiquitous across catalysis, corrosion science, soil science, and subsurface geochemistry. At relative humidity (RH) values below saturation, these films are Å–nm thick and exhibit properties that differ sharply from bulk water, including disrupted H-bonding and impeded mass transport. Owing to their high surface-to-volume ratio, dissolution of the solid can rapidly drive strong supersaturation with respect to secondary phases. Reactivity in thin water films is highly sensitive to film thickness, and critically, thickness evolves during reaction because the hygroscopicity of the interfacial system changes as ions accumulate or diminish in the film and as reaction products transform. To accurately probe and control these dynamics, a capability is needed that can measure and automatically maintain a constant water-film thickness while simultaneously monitoring solid dissolution, nucleation, and growth. This project developed an autonomous Infrared/Small Angle X-ray Scattering-Wide Angle X-ray Scattering (IR/(W)SAXS) for investigating reactivity in thin water films on solid surfaces exposed to humidified gases. The capability consists of an IR spectrometer, a (W)SAXS instrument, and a mass flow controller system for generating variably humidified gas flows to a custom reaction cell. Progress on each of the major components of the capability are detailed below.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2003

An Experimental Investigation of Unsteady Surface Pressure on an Airfoil in Turbulence

Measurements of fluctuating surface pressure were made on a NACA 0015 airfoil immersed in grid generated turbulence. The airfoil model has a 2 ft chord and spans the 6 ft Virginia Tech Stability Wind Tunnel test section. Two grids were used to investigate the effects of turbulence length scale on the surface pressure response. A large grid which produced turbulence with an integral scale 13% of the chord and a smaller grid which produced turbulence with an integral scale 1.3% of the chord. Measurements were performed at angles of attack, alpha from 0 to 20 . An array of microphones mounted subsurface was used to measure the unsteady surface pressure. The goal of this measurement was to characterize the effects of angle of attack on the inviscid response. Lift spectra calculated from pressure measurements at each angle of attack revealed two distinct interaction regions; for omega(sub r) = omega b / U(sub infinity) is less than 10 a reduction in unsteady lift of up to 7 decibels (dB) occurs while an increase occurs for omega(sub r) is greater than 10 as the angle of attack is increased. The reduction in unsteady lift at low omega(sub r) with increasing angle of attack is a result that has never before been shown either experimentally or theoretically. The source of the reduction in lift spectral level appears to be closely related to the distortion of inflow turbulence based on analysis of surface pressure spanwise correlation length scales. Furthermore, while the distortion of the inflow appears to be critical in this experiment, this effect does not seem to be significant in larger integral scale (relative to the chord) flows based on the previous experimental work of McKeough suggesting the airfoils size relative to the inflow integral scale is critical in defining how the airfoil will respond under variation of angle of attack. A prediction scheme is developed that correctly accounts for the effects of distortion when the inflow integral scale is small relative to the airfoil chord. This scheme utilizes Rapid Distortion Theory to account for the distortion of the inflow with the distortion field modeled using a circular cylinder.

Mish, Patrick F.↗
NASA NTRS2005

Lunar and Planetary Science XXXVI, Part 15

Contents include the following: Impact Metamorphism of Subsurface Organic Matter on Mars: A Potential Source for Methane and Surface Alteration. Preliminary Study of Polygonal Impact Craters in Argyre Region, Mars. Geochemistry of the Dark Veinlets in the Granitoids from the Souderfjarden Impact Structure, Finland: Preliminary Results. An Experimental Method to Estimate the Chemical Reaction Rate in Vapor Clouds: An Application to the K/T Impact. Study of the Apollo 16 Landing Site: Re-Visit as a Standard Site for the SELENE Multiband Imager. First X-Ray Observation of Lunar Farside from Hayabusa X-Ray Spectrometer. Lunar X-Ray Fluorescence Spectrometry from SELENE Lunar Polar Orbiter. Origin and Thermal History of Lithic Materials in the Begaa LL3 Chondrite. Evidence of Normal Faulting and Dike Intrusion at Valles Marineris from Pit Crater Topography. Evidence of Tharsis-Radial Dike Intrusion in Southeast Alba Patera from MOLA-based Topography of Pit Crater Chains. Are They Really Intact? Evaluation of Captured Micrometeoroid Analogs by Aerogel at the Flyby Speed of Stardust. Numerical Simulations of Impactor Penetration into Ice-Over-Water Targets. A Probable Fluid Lava Flow in the Hebes Mensa (Mars) Studied by HRSC Images. New Drill-Core Data from the Lockne Crater, Sweden: The Marine Excavation and Ejection Processes, and Post-Impact Environment. Cross-Sectional Profile of Baltis Vallis Channel on Venus: Reconstruction from Magellan SAR Brightness Data.

Source record↗
DOE OSTI2025

Permeability scaling relationships of volcanic tuff from core to field scale measurements

A recent chemical explosive test in P-Tunnel at the Nevada National Security Site, Nevada, USA, was conducted to better understand how signals propagate from explosions in the subsurface. A primary signal of interest is the migration of gases that can be used to differentiate chemical from nuclear explosions. Gas migration is highly dependent on the rock permeability which is notoriously difficult to determine experimentally in the field due to a potentially large dependence on the scale over which measurements are made. Here, we present pre-explosion permeability estimates to characterize the geologic units surrounding the recent test. Permeability measurements were made at three scales of increasing size: core samples (≈2 cm), borehole packer system tests (≈1 m), and a pre-shot cavity pressurization test (> 10 m) across ten tuff units. Permeability estimates based on core measurements showed little difference from borehole packer tests. However, permeability in most rock units calibrated from cavity pressurization tests resulted in higher permeability estimates by up to two orders of magnitude. Here, we demonstrate that the scale of the measurement significantly impacts the characterization efforts of hydraulic properties in volcanic tuff, and that local-scale measurements (< 10 m scale) do not incorporate enough heterogeneity to accurately predict field-scale flow and mass transport.

Environmental sciences↗
NASA NTRS1999

Thorium Enrichment within the Procellarum KREEP Terrane: The Record in Surface Deposits and Significance for Thermal Evolution

The nearside-farside structural and compositional asymmetry of the Moon was recognized during the early days of Apollo and the suggestion was made that the migration of mantle melts to the nearside would have been favored by early Earth-Moon orbital dynamics and nonuniform planetesimal bombardment. Recent global geochemical mapping by Lunar Prospector has provided additional data, particularly in the Th distribution, that strongly supports the notion of global, preferential melt migration, which led in part to the development of the Procellarum KREEP Terrane (PKT) [2-5]. The surface distribution of Th was then reshaped by basin-forming impacts into the PKT, especially the Imbrium impact, which was the last and largest to strike in that region. The Imbrium event probably excavated material from a partially molten zone deep in the crust and delivered Th-rich ejecta Moon-wide. A fundamentally important but poorly understood aspect of the global Th distribution is the concentration of Th in the subsurface rocks of the PKT crustal section. For example, depending on what assumptions are made, the PKT crustal section, which is about 12% of the crust and only about 1.2% of the whole Moon, may contain as much as 40% of the Moon's entire Th budget. Such a distribution of Th and related heat-producing elements would have had a profound effect on melting, mixing, and the thermal evolution of the PKT and the underlying mantle. In this abstract, we examine the compositions of terra formations within the PKT and relate them to some of the Th-bearing rock types known from the Apollo samples. It appears that the existence of the PKT may be a unifying concept for a number of petrologic and geochemical observations. From the initial Lunar Prospector gamma-ray spectrometer data(-about 5 deg. resolution) and from the preliminary low-orbit data, there appears to be a number of relatively hotter "spots" within the PKT in terms of Th concentration. Some of the hotter spots correspond to intermediate-sized craters that penetrated volcanic flows and excavated Th-rich, submare material, such as Aristarchus, Aristillus, and Kepler. Other spots, however, correspond to surficial formations that constitute mainly rough topography associated with Imbrium ejecta or circum-Imbrium ring mountains and do not necessarily imply the presence of exposed KREEP basalts. The Fra Mauro Formation south of Copernicus toward the Apollo 14 site and regions of the Alpes Formation southwest of Copernicus in the vicinity of Reinhold lie within the most prominent hotspot; here them concentration is consistent with that found in the Apollo 14 soils (about 12-13 ppm) when the proportions of Fra Mauro Formation and mare basalt are considered. The area between Copernicus and Kepler and northwest of Copernicus in terra extending to the Carpathians is similarly enriched in Th. The Apennines from Eratosthenes toward the Apollo 15 site contain elevated Th concentrations, as does the northwestern quadrant of circum-Imbrium terra, especially between (but not including) La Condamine and Plato, and in the region northwest of the Jura mountains extending southward past Mairan to the Gruithuisen-Domes region. Within the main topographic rim of Imbrium, the Apennine Bench formation south of Archimedes appears to have relatively elevated Th concentration. Comparing the map of Th distribution to a digital-elevation map derived from Clementine altimetry, it appears that most of the areas richest in Th occur where the surface is elevated relative to the majority of PKT volcanic plains. Not all rough topography within the PKT has such elevated Th, however. Based on an analysis of the 5 deg. data, and using the calibration of, the mean Th concentrations for mainly volcanic-resurfaced terrain and rugged terrain are similar (about 5.5 ppm). This occurs in part because craters that penetrated mare basalt excavated Th-rich material. Even so, there appear to be extensive areas of volcanic resurfacing that have no obvious extrinsic source of Th-rich material, suggesting that the basalts, themselves, may contain as much as 5-6 ppm Th. High-FeO concentrations (18 to >20 wt%) indicate that these are not KREEP basalts but mare basalts. If so, this is surprising because most of the Apollo-sampled mare basalts have very low-Th (typically <2 ppm). A variety of Th-rich materials occur in the sample collection, particularly in the samples from the Apollo 12,14, and 15 sites. The most abundant Th-rich rock types are the mafic impact-melt breccias, which although found at all sites, are most abundant at Apollo 14, where they dominate the rock samples and make up some 40% of the rock particles in the soil. These have Th concentrations ranging up to about 30 ppm and averaging about 18 ppm. The Apollo 14 soils contain about 13 ppm Th, reflecting the high abundance of this melt-breccia component. At the nearby Apollo 12 site, the rocks consist mainly of mare basalts, and these have low-Th concentrations, mostly <1 ppm. Among the nonmare rocks, however, a few such as complex breccia 12013 contain a variety of evolved lithologies and represent potential sources of Th-rich components ranging from 17 to 50 ppm. A single fragment of KREEP basalt from the Apollo 12 soil contains about 50 ppm Th. Despite the identification of highly Th-enriched lithologic components at the Apollo 12 site, the soils there vary linearly in composition so as to extrapolate to a moderate Th KREEP-basalt component or a composition like that of the Fra Mauro formation as reflected by Apollo 14 soil (14 ppm Th at 10 wt% FeO). The nearby Lansberg Crater (39-km diameter) is a likely candidate to have delivered submare material such as buried Fra Mauro or Alpes material to the Apollo 12 site. Additional information contained in the original.

Jolliff, B. L.↗
NASA NTRS1983

Crustal deformation, the earthquake cycle, and models of viscoelastic flow in the asthenosphere

The crustal deformation patterns associated with the earthquake cycle can depend strongly on the rheological properties of subcrustal material. Substantial deviations from the simple patterns for a uniformly elastic earth are expected when viscoelastic flow of subcrustal material is considered. The detailed description of the deformation pattern and in particular the surface displacements, displacement rates, strains, and strain rates depend on the structure and geometry of the material near the seismogenic zone. The origin of some of these differences are resolved by analyzing several different linear viscoelastic models with a common finite element computational technique. The models involve strike-slip faulting and include a thin channel asthenosphere model, a model with a varying thickness lithosphere, and a model with a viscoelastic inclusion below the brittle slip plane. The calculations reveal that the surface deformation pattern is most sensitive to the rheology of the material that lies below the slip plane in a volume whose extent is a few times the fault depth. If this material is viscoelastic, the surface deformation pattern resembles that of an elastic layer lying over a viscoelastic half-space. When the thickness or breath of the viscoelastic material is less than a few times the fault depth, then the surface deformation pattern is altered and geodetic measurements are potentially useful for studying the details of subsurface geometry and structure. Distinguishing among the various models is best accomplished by making geodetic measurements not only near the fault but out to distances equal to several times the fault depth. This is where the model differences are greatest; these differences will be most readily detected shortly after an earthquake when viscoelastic effects are most pronounced.

Cohen, S. C.↗