SEARCH · Engineering Papers

Results for “gradient sampling”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

1,346 records · Page 75

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Induction and repair of HZE induced cytogenetic damage

Wistar rats were exposed to high-mass, high energy (HZE) 56Fe particles (1000 GeV/AMU) using the Alternating Gradient Synchrotron (AGS). The animals were sacrificed at 1-5 hours or after a 30-day recovery period. The frequency of micronuclei in the tracheal and the deep lung epithelial cells were evaluated. The relative effectiveness of 56Fe, for the induction of initial chromosome damage in the form of micronuclei, was compared to damage produced in the same biological system exposed to other types of high and low-LET radiation. It was demonstrated that for animals sacrificed at short times after exposure, the tracheal and lung epithelial cells, the 56Fe particles were 3.3 and 1.3 times as effective as 60Co in production of micronuclei, respectively. The effectiveness was also compared to that for exposure to inhaled radon. With this comparison, the 56Fe exposure of the tracheal epithelial cells and the lung epithelial cells were only 0.18 and 0.20 times as effective as radon in the production of the initial cytogenetic damage. It was suggested that the low relative effectiveness was related to potential for 'wasted energy' from the core of the 56Fe particles. When the animals were sacrificed after 30 days, the slopes of the dose-response relationships, which reflect the remaining level of damage, decreased by a factor of 10 for both the tracheal and lung epithelial cells. In both cases, the slope of the dose-response lines were no longer significantly different from zero, and the r2 values were very high. Lung epithelial cells, isolated from the animals sacrificed hours after exposure, were maintained in culture, and the micronuclei frequency evaluated after 4 and 6 subcultures. These cells were harvested at 24 and 36 days after the exposure. There was no dose-response detected in these cultures and no signs of genomic instability at either sample time.

Non-NASA Center↗

DOE OSTI2025

Sharp detection of low-dimensional structure in probability measures via dimensional logarithmic Sobolev inequalities

Identifying low-dimensional structure in high-dimensional probability measures is an essential pre-processing step for efficient sampling. To identify this structure, we approximate the target measure as a perturbation of an arbitrary reference measure along a few directions in $\mathbb{R}^{d}$. These directions are determined by minimizing an upper bound on the Kullback–Leibler (KL) divergence between the target and its approximation. Our contribution improves upon previous works by leveraging dimensional logarithmic Sobolev inequalities to refine the bound on the KL divergence. These inequalities lead to a uniformly tighter bound on the KL divergence, thereby enhancing the identification of the most significant perturbation directions. In particular, when the target and reference are both Gaussian, minimizing the resulting bound is equivalent to minimizing the KL divergence. We further demonstrate the applicability of this analysis to the squared Hellinger distance, where analogous reasoning shows that the dimensional Poincaré inequality offers improved bounds.

Bayesian inference↗

NASA NTRS1999

Simulating the Formation of Lunar Floor-Fracture Craters Using Elastoviscoplastic Relaxation

Lunar floor-fracture craters formed during the height of mare basalt emplacement. Due to a general temporal and spatial relation with the maria, these craters, numbering some 200, may be diagnostic of the thermal structure of the crust during this time. As the name suggests, these craters exhibit brittle failure, generally limited to the central floor region. That, and a shallower depth than fresh lunar craters, has led to two main theories as to their formation: laccolith emplacement under the crater and viscous relaxation. The implications of each model for the state of the Moon's crust during this time are quite different, so the viability of each model must be checked. Laccolith emplacement has been treated elsewhere. However, previous attempts to study the relaxation of the craters have assumed only a uniform, Newtonian viscous response of the near surface to the topographic driving forces, and simply postulated that the fractures resulted from tensile stresses associated with floor uplift. Here, we use a more sophisticated rheological model that includes not only non-Newtonian viscous behavior (i.e., the viscosity is stress-dependent), but also incorporates elastic behavior and a plastic component to the rheology to directly simulate the formation of the floor fractures. The results of our simulations show that while elastoviscoplastic relaxation is potentially viable for larger floor-fracture craters, it is not viable for craters with diameters < or = 60 km, the size of the majority of floor-fracture craters. We employ the finite element method, a numerical technique well suited for boundary-value problems, via the commercially available MARC software package. To test the viability of topographic relaxation, our goal is to prepare the simulations as to maximize the amount of relaxation. We take advantage of the natural axisymmetry of craters, simulating one radial plane. Initial shapes are based on data for fresh craters from Pike. To simplify implementation, a fourth order polynomial is used for the basin, while a third order inverse function is used for the rim. This form closely approximates the long-wavelength behavior of complex craters, while ignoring higher-frequency topography, save the rim. This approximation is appropriate because crater relaxation is strongly controlled by long-wavelength topography. Loading is accomplished assuming a uniform gravity field (1.62 m/s-square) and a uniform density of 2900 kg/cubic m. The initial stress state is set to be hydrostatic, with an additional pressure term to account for any overlying topography. This additional pressure term is tapered exponentially with depth. While the simulations quickly settle on a preferred stress state, and while the final solution is fairly insensitive to the choice of the e-folding depth of the taper, selecting an e-folding depth close to the diameter of the crater sets the initial stress state near the preferred state. We assume a diuranally averaged surface temperature of -20 C, and allow temperature to increase with depth at a rate of 50 K//km. Assuming a thermal conductivity of 2 W/in/K, this gradient translates to a heat flow of 100mW/square m, an extremely high value for the Moon. Temperature, of course, will not increase without bound. To maximize relaxation, we allow our temperature profile to increase linearly until it reaches the solidus (assumed to be 1200C) at a depth of 24.4 km, at which point it is kept constant. The presence of melt will drop the bulk viscosity; however, we have no rheological control for partial melts. Therefore, we make no attempt to simulate this situation. Elastoviscoplastic rheological model. In general, geologic materials can behave in three main ways: elastically, viscously (via solid-state creep), and brittly (plasticity is a continuum approach to simulate this phenomenon). We combine these three deformation mechanisms in an extended Maxwell solid, where the total strain can be broken down into a simple summation of the elastic, creep, and plastic strains. In relaxation phenomena in general, the system takes advantage of any means possible to eliminate deviatoric stresses by relaxing away the topography. Previous analyses have only modeled the viscous response. Comparatively, the elastic response in our model can augment the relaxation, to a point. This effect decreases as the elastic response becomes stiffer; indeed, in the limit of infinite elastic Young's modulus (and with no plasticity), the solution converges on the purely viscous solution. Igneous rocks common to the lunar near-surface have Young's modulii in the range of 10-100 GPa. To maximize relaxation, we use a Young's modulus of 10 GPa. (There is negligible sensitivity to the other elastic modulus, the Poisson's ratio; we use 0.25.) For the viscous response, we use a flow law for steady-state creep in thoroughly dried Columbia diabase, because the high plagioclase (about 70 vol%) and orthopyroxene (about 17 vol%) content is similar to the composition of the lunar highland crust as described by remote sensing and sample studies: noritic anorthosite. This flow law is highly non-Newtonian, i.e., the viscosity is highly stress dependent. That, and the variability with temperature, stands in strong contrast to previous examinations of lunar floor-fracture crater relaxation. To model discrete, brittle faulting, we assume "Byerlee's rule," a standard geodynamical technique. We implement this "rule" with an-angle of internal friction of about 40 deg, and a higher-than-normal cohesion of about 3.2 MPa (to approximate the breaking of unfractured rock). The actual behavior of geologic materials is more complex than in our rheological model, so the uncertainties in the plasticity do not represent the state-of-the-art error. Additional information is contained in the original.

Dombard, A. J.↗

NASA NTRS2012

Paleo-Environment and C-14 Dating: The Key to the Depositional Age of the Tha Chang and Related Sand Pits, Northeastern Thailand

Tha Chang sand pits, Nakhon Ratchasima Province and many other sand pits in the area adjacent to the Mun River are characterized by their fluviatile environment in association with mass wasting deposits, along the paleo-river channel and the flood plain of the Mun River. Sediments of these deposits are characterized by clasts of various rock types especially the resistant ones with frequent big tree trunks, logs and wood fragments in different sizes and various stages of transformation from moldering stage to lignification and petrification. Widespread pyritization of the lower horizon suggests strongly reducing environment during burial. The Tha Chang deposits have been received much attention from geoscientists especially paleontologist communities, as they contain fragments of some distinct vertebrate species such as Stegadon sp., hominoid primate, rhinoceros Aceratherium and others. Based on the associated mammal fauna and hominoid fossils, the late Miocene ( 9 - 6 Ma) was given for the time of deposition of this sand and gravel unit. Some other reports believed that sediments and materials of these sand and gravel quarries (pits) were deposited by high-energy flood pulses contemporaneous with the tektites forming event during mid-Pleistocene at c. 0.8 Ma. Interpretation from Palynostratigraphical study suggested that the lower horizon of Tha Chang sand pit was deposited during Pliocene/Pleistocene period and the upper horizons are Pleistoncene/Holocene. It is crystal clear that all the fluviatile sediments including tektites and almost all fossil fragments being deposited in these sand pits were, likely a multiple times reworked materials. Only some old bamboo trees, some old crowling trees and fossils grasses observed on the old river bank are considered in situ. C-14 dating of 5 old wood specimens from Tha Chang Sand Pits, 15 old wood specimens from Chumpuang Sand Pits and one sample of old pottery from a Chumpuang Sand Pit were carried out in the NSF- Arizona AMS Laboratory. Although, there is no sharp boundary between the unconsolidated sedimentary horizons in the pits, C-14 ages obtained from the Tha Chang vary from 34,340 BP at the middle horizon (approx 10 m below ground zero) to >49,900 BP at the lower horizon with unknown basal formation (highly pyritized zone approx 20 - 25 m below ground zero). The ages for the Chumpuang vary from 41,700 BP, >45,900 BP and >49,900 BP from the upper most to the lower most of a broad horizon (approx 8 m to approx 12 m below ground zero). The C-14 age of the pottery collected from layer approximately 5 m below ground zero is 2,514 BP. The nature of fluviatile together with occasional mass wasting characteristics of all sand pits studies suggest the relatively faster depositional rate of the lower horizon which involved more flooding and mass wasting deposits than those of the upper horizons. The apparent of some mixing of the wood ages may indicate reworking and lag deposits nature of the area. The depositional rate of the upper most sand and soil horizon (5 m thick) is approximately 1 m per 500 years which mean both erosion and deposition had played a significant role during that time period. In term of the true age of the formation, we argue that since most of the materials deposited are reworked materials, all ages obtained from fossil fragments could not be the age of sand and gravel formation. Furthermore, the maximum age of all the tektite bearing horizons cannot be older than 0.8 Ma. The oldest C-14 age of 49,900 BP is interpreted as the minimum age of the Tha Chang and related sand pits formation when geomorphology of the area was a lot more hilly and much higher gradient than that of the present day.

Putthapiban, P.↗

DOE OSTI2024

Greenhouse gas and dissolved oxygen production and consumption rates associated with Regier et al. (2023)

Tidal inundation along the coastal terrestrial-aquatic interface controls soil and sediment biogeochemistry and gas dynamics. Although a rich literature exist on studies of the influence of tidal waters on the biogeochemistry of coastal ecosystem soils, few studies have experimentally addressed the reverse question: How do soils (or sediments) from different coastal ecosystems influence the biogeochemistry of the tidal waters that inundate them? We conducted short-term microcosm laboratory experiments where seawater was amended with sediments and soils collected across regional gradients of inundation exposure (i.e., frequently to rarely inundated) and measured changes in dissolved oxygen and greenhouse gas concentrations to calculate gas consumption or production rates occurring during seawater exposure to terrestrial materials. This data package contains dissolved oxygen and greenhouse gas data collected during incubation of soils and sediments collected at 18 sites, which were used in the publication Regier et al. (2023) entitled “Coastal inundation regime moderates the short-term effects of sediment and soil additions on seawater oxygen and greenhouse gas dynamics: a microcosm experiment” which is published in Frontiers in Marine Science (DOI: https://doi.org/10.3389/fmars.2023.1308590).---Acknowledging EXCHANGE: General Support and Data Product UseWe ask that users of EXCHANGE data add the following acknowledgement when publishing data in scholarly articles and data repositories:"This research is based on work supported by COMPASS-FME, a multi-institutional project supported by the U.S. Department of Energy, Office of Science, Biological and Environmental Research as part of the Environmental System Science Program."

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

Expedition UT-GOM2-2 Site H

Pressure and conventional cores were collected at Site H of the Walker Ridge Protracted Area Block 313 in the Terrebonne Basin, deepwater Gulf of America (Gulf of Mexico) during the University of Texas (UT) Deepwater Hydrate Coring Expedition (UT-GOM2-2). Pressure and conventional cores were collected continuously to a depth of 155.1 meters below the seafloor (mbsf). At deeper depths, cores were taken periodically from hydrate-bearing sands and their bounding muds to a total depth of 861.3 mbsf. 162.6 m of conventional core and 54.8 m of pressure core were recovered. Twelve temperature measurements were made between 27.1 and 144.5 mbsf to determine the geothermal gradient. At the seafloor, more than 4 m of sandy silt of unknown origin was encountered. Beneath this sand, to a depth of about 200 mbsf, the section was composed of interbedded mud and biogenic carbonate ooze. The ooze correlated to low density and high porosity intervals observed in the previously acquired logging while drilling (LWD) data and as measured. These ooze intervals also correspond to lighter sediment color, increased Ca content based on X-ray florescence (XRF) core scanning, and increased calcareous nannofossil abundance. Calcareous nannofossil biostratigraphy constrains the entire record to the Pleistocene (< 0.91 million years), with a pronounced increase in sedimentation rate with depth. Below 200 mbsf, the section was predominantly composed of mud with two thicker, hydrate-bearing coarse-grained intervals, which are commonly known as the Blue and Orange sands. The dissolved gas concentration was quantified from pressure cores. In the shallow section, dissolved methane concentration increased below the sulfate-methane transition zone (SMTZ) and reaches saturation (the limit of solubility for methane) at 147 mbsf. Gas expansion was very common in conventional and depressurized pressure (conventionalized) cores below the SMTZ. At deeper depths, the methane concentration within muds bounding the Blue and Orange reservoirs was generally found to be less than saturation. The dissolved and hydrate gas composition is consistent with a microbial source, containing greater than 99.99% methane and only trace concentrations of ethane, propane, and butane. The methane to ethane ratio (C1 /C2 ) and the methane to ethane plus propane (C1 /(C2 +C3 )) decrease with depth down to at least 678 mbsf, mainly driven by the increase in ethane with depth. It is unclear if this trend continues through the Orange sand interval. The δ13C isotopic signature of methane ranges between -69.9 and -78.5 ‰ relative to the Vienna Pee Dee Belemnite (VPDB) standard. Pressure core recovery in sandy intervals was poor. However, pressure core logs of the Orange sand show intervals of low density and high velocity, which are indicative of high hydrate saturation. One pressure core was degassed and the average hydrate saturation in the core was determined to be 24%. One core from within the Orange sand was composed of interbedded graded sandy silt and mud. The sandy silts from this core are composed of mainly quartz and feldspar with some lithics. Most of the recovered pressure core samples are maintained at near in-situ pressure and temperature (within the hydrate stability field) at the University of Texas Pressure Core Center awaiting analysis. In the shallow section, samples will be used to determine the flux of organic carbon through the basin system, find the rate at which that carbon was consumed, and understand the microbial population responsible for these processes. In the deeper section, samples from in and around the hydrate reservoirs will be used to determine the petrophysical properties of the reservoir and bounding seals in these systems.

03 NATURAL GAS↗

DOE OSTI2025

Data from a throughfall exclusion experiment: Fine root dynamics, morphology, chemistry, and AMF colonization across four lowland Panamanian forests

Fine roots regulate forest nutrient, carbon, and water cycling, yet their variation within and among tropical forests remains under-characterized. We quantified root productivity, disappearance, and stocks to 1 m using minirhizotron imaging, and we measured morphology, elemental composition [root carbon (C), root nitrogen (N), root phosphorus (P)], and arbuscular mycorrhizal fungi (AMF) colonization to 20 cm using ingrowth cores and sequential coring. Sampling took place in four distinct lowland Panamanian forests (32 plots; 8 per forest) from 2018 through 2022 under control and throughfall-exclusion (drought) treatments in the Panama Rainforest Changes with Experimental Drying (PARCHED) experiment.The dataset is presented as an Excel workbook with six tabs. The first tab is the data dictionary. Tab S1 contains ingrowth-core production and mortality, morphology and soil moisture. Tab S2 contains sequential-coring standing stocks with associated morphology and soil moisture. Tab S3 contains minirhizotron row data records to 1 m depth, including per-frame root length and diameter, normalized length metrics, and session timing. Tab S4 contains AMF colonization. Tab S5 contains fine-root chemistry at 0–10 cm, reporting %P, %C, %N, and C:N for samples collected via ingrowth cores and sequential-coring standing stocks. CSV mirrors for each tab are provided, and a KML file supplies coordinates for all 32 plots.Key variables span live and dead fine-root biomass (and coarse fractions where applicable), specific root length (SRL) and area (SRA), diameter, root tissue density (RTD), soil moisture, AMF colonization, root %N, %C, %P, and C:N, along with minirhizotron root length and diameter. Depth, season, treatment, and plot/site identifiers are included to support cross-tab integration and analysis from 0–100 cm (minirhizotron) and 0–20 cm (cores).Units are reported in-column and missing values are coded as NA. No special software is required to open or use the files (Excel, CSV, and KML compatible).

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Foam Fracturing Study for Stimulation Development of Enhanced Geothermal Systems

The large thermal gradients and high subsurface temperatures of the western region of the U.S. hold great potential for the implementation of enhanced geothermal systems (EGS). The development of these potential EGS resources requires stimulation of the reservoir to enhance permeability and it has been widely reported that a substantial amount of water will be required should conventional hydraulic stimulation be used. This presents a huge challenge and a high risk to the geothermal development because the water stress1 in these areas is already high or extremely high. The use of foam, a gas/liquid mixture predominantly composed of gas, in fracturing is considered and explored in this project as a potential approach to address water concerns with hydraulic stimulation in the development of EGS. This project, led by Oak Ridge National Laboratory (ORNL) in collaboration with Temple University, was awarded in an open lab call in 2018, and was part of the DOE GTO waterless stimulation initiative. The goal of the project was to demonstrate the feasibility of foam fracturing for EGS development through two primary tasks: Task 1: Laboratory study of the effectiveness of foam fracturing for representative geological materials, including cyclic pressurization using foam (led by ORNL) and Task 2: High temperature foam material selection and characterization (led by Temple University). In FY19, ORNL finished the critical review on serval issues associated with foam fracturing and the implementation of the proposed tasks in a lab study (Wang, et al., 2019), and completed the foam fracturing testing using cement as a model material (Wang, et al, 2020a). The work at ORNL was geared up to develop a brand-new foam testing system in FY20. The purchase of main components for the new system was finished in the first half of the FY20. The assembly of the foam testing system and foam fracturing testing were completed in the second half of the FY20 (Wang, et al., 2021a). Task 1 required the development of a test system which can be used to perform hydraulic fracturing of geological specimens with both water and foamed liquids at pressure up to 6,000 psi (41.4 MPa). The system possesses several capabilities that conventional injection systems lack for hydraulic fracturing. In addition to its ability to generate foam with controlled quality, it is capable of cycling pressure levels between specified values up to frequencies of 50 Hz. The latter capability was developed to evaluate the hypothesis that cyclic loading of samples would produce enhanced fracturing. The system consists of two sections: one for foam generation and another for foam injection. The foam is generated through separate control and pressurization of liquid and gas phases with controlled flow rates. The injection section is equipped with a low-flow Coriolis flow that monitors the density of foam to ensure the injection is in the range of target foam quality2. Experimental results of foam fracturing are reported for cylindrical granite specimens using water and aqueous N2 foam as the fracturing fluids. All experiments were performed for unconfined conditions. The effects of injection mode (i.e., monotonic vs cyclic pressurization) on breakdown pressure and failure response sample were investigated using water alone as a fracturing fluid and foams with a range of compositions. It was found that in the case of monotonic injection, the breakdown pressure of granite specimens tended to be slightly higher when fracturing with foam. Additionally, with a foam quality of 90%, the water use can be reduced by 50 to 84%, depending on hole size. On the other hand, it was observed that the breakdown pressure can be brought down to 70% of the monotonic breakdown pressure by using low cycle fatigue. Finally, discussions are presented regarding injectivity and water use reduction.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

Baltimore Ecosystem Study: Soil atmosphere fluxes of carbon dioxide, nitrous oxide and methane, 1998 - ongoing

The Baltimore Ecosystem Study (BES) established a network of long-term permanent biogeochemical study plots in 1998. These plots provide long-term data on vegetation, soil and hydrologic processes in the key ecosystem types within the urban ecosystem. The network of study plots includes forest plots (upland and riparian), chosen to represent the range of forest conditions in the area and grass plots (to represent home lawns). Plots are instrumented with lysimeters (drainage and tension) to sample soil solution chemistry, time domain reflectometry probes to measure soil moisture, dataloggers to measure and record soil temperature, and trace gas flux chambers to measure the flux of carbon dioxide, nitrous oxide and methane from soil to the atmosphere. Measurements of in situ nitrogen mineralization, nitrification and denitrification were made at approximately monthly intervals from Fall 1998 - Fall 2000. Detailed vegetation characterization (all layers) was done in summer 1998 and 2015. Data from these plots has been published in Groffman et al. (2006, 2009), Groffman and Pouyat (2009), Savva et al. (2010), Costa and Groffman (2013), Duncan et al. (2013), Waters et al. (2014), Ni and Groffman (2018), Templeton et al. (2019). Literature Cited Costa, K.H. and P.M. Groffman. 2013. Factors regulating net methane flux in urban forests and grasslands. Soil Science Society of America Journal 77:850 - 855. Duncan, J. M., L. E. Band, and P. M. Groffman. 2013. Towards closing the watershed nitrogen budget: Spatial and temporal scaling of denitrification. Journal of Geophysical Research Biogeosciences 118:1-5; DOI: 10.1002/jgrg.20090 Groffman PM, Pouyat RV, Cadenasso ML, Zipperer WC, Szlavecz K, Yesilonis IC,. Band LE and Brush GS. 2006. Land use context and natural soil controls on plant community composition and soil nitrogen and carbon dynamics in urban and rural forests. Forest Ecology and Management 236:177-192. Groffman, P.M., C.O. Williams, R.V. Pouyat, L.E. Band and I.C. Yesilonis. 2009. Nitrate leaching and nitrous oxide flux in urban forests and grasslands. Journal of Environmental Quality 38:1848-1860. Groffman, P.M. and R.V. Pouyat. 2009. Methane uptake in urban forests and lawns. Environmental Science and Technology 43:5229-5235. DOI: 10.1021/es803720h. Ni, X. and P.M. Groffman. 2018. Declines in methane uptake in forest soils. Proceedings of the National Academies of Science of the United States of America 115:8587-8590. Savva, Y., K. Szlavecz, R. V. Pouyat, P. M. Groffman, and G. Heisler. 2010. Effects of land use and vegetation cover on soil temperature in an urban ecosystem. Soil Science Society of America Journal 74:469-480. Templeton, L., M.L. Cadenasso, J. Sullivan, M. Neel and P.M. Groffman. 2019. Changes in vegetation structure and composition of urban and rural forest patches in Baltimore from 1998 to 2015. Forest Ecology and Management. In press. Waters, E.R., J.L. Morse, N.D. Bettez and P.M. Groffman. 2014. Differential carbon and nitrogen controls of denitrification in riparian zones and streams along an urban to exurban gradient. Journal of Environmental Quality 43:955–963.

Groffman, Peter M↗

NASA NTRS1999

Containerless, Low-Gravity Undercooling of Ti-Ce Alloys in the MSFC Drop Tube

Previous tests of the classical nucleation theory as applied to liquid-liquid gap miscibility systems found a discrepancy between experiment and theory in the ability to undercool one of the liquids before the L1-L2 separation occurs. To model the initial separation process in a two-phase liquid mixture, different theoretical approaches, such as free-energy gradient and density gradient theories, have been put forth. If there is a large enough interaction between the critical liquid and the crucible, both models predict a wetting temperature (T(sub w)) above which the minority liquid perfectly wets and layers the crucible interface, but only on one side of the immiscibility dome. Materials with compositions on the other side of the dome will have simple surface adsorption by the minority liquid before bulk separation occurs when the coexistence (i.e., binoidal) line in reached. If the interaction between the critical liquid and the crucible were to decrease, T(sub w) would increase, eventually approaching the critical consolute temperature (T(sub cc)). If this situation occurs, then there could be large regions of the miscibility gap where non-perfect wetting conditions prevail resulting in droplets of L1 liquid at the surface having a non-zero contact angle. The resulting bulk structure will then depend on what happens on the surface and the subsequent processing conditions. In the past several decades, many experiments in space have been performed on liquid metal binary immiscible systems for the purpose of determining the effects that different crucibles may have on the wetting and separation process of the liquids. Potard performed experiments that showed different crucible materials could cause the majority phase to preferentially wet the container and thus produce a dispersed microstructure of the minority phase. Several other studies have been performed on immiscibles in a semi-container environment using an emulsion technique. Only one previous study was performed using completely containerless processing of immiscible metals and the results of that investigation are similar to some of the emulsion studies. In all the studies, surface wetting was attributed as the cause for the similar microstructures or the asymmetry in the ability to undercool the liquid below the binoidal on one side of the immiscibility dome. By removing the container completely from the separation process, it was proposed that the loss of the crucible/liquid interaction would produce a large shift in T(sub w) and thus change the wetting characteristics at the surface. By investigating various compositions across the miscibility gap, a change in the type and amount of liquid wetting at the surface of a containerless droplet should change the surface nucleating behavior of the droplet - whether it be the liquid-liquid wetting or the liquid-to-solid transition. Undercooling of the liquid into the metastable region should produce significant differences in the separation process and the microstructure upon solidification. In this study, we attempt to measure these transitions by monitoring the temperature of the sample by optical pyrometry. Microstructural analysis will be made to correlate with the degree of undercooling and the separation mechanisms involved.

Robinson, M. B.↗

DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2025

Abiotic influences on continuous conifer forest structure across a subalpine watershed

Understanding the abiotic drivers of high-elevation forest physiognomy is essential for forecasting how mountain ecosystems will respond to emerging environmental pressures. Most prior studies of these relationships have relied on small samples of the full landscape, resulting in limited power to detect dominant covariates and their interactions. Here we report the first evaluation of abiotic influences on a complement of accurate, wall-to-wall estimates of conifer forest structure and composition at the watershed scale. In a subalpine conifer domain in the Colorado Rocky Mountains (USA), we developed a novel method for deriving stand structure metrics from waveform LiDAR data, which showed high fidelity with field inventory. We quantified the relationships between structural and compositional metrics and climate, topographic, edaphic, and geologic factors. Our results showed that peak snow water equivalent (SWE), snow disappearance rate, and elevation explained most of the variation in forest structure. The highest stand density, basal area, maximum canopy height, and quadratic mean diameter occurred in sites with SWE around one standard deviation below mean, but with long snow residence times. Stand density decreased linearly with elevation, while other metrics peaked between 3000 m and 3200 m. Substrate properties had weaker influence. Continuous mapping of through-canopy forest structure enabled our novel findings of the dominant role of snowpack in explaining structural and compositional variation, and of elevation thresholds. Our reproducible approach facilitates assessment of forest-topoclimate relationships in other conifer-dominated landscapes and improves understanding of the baseline patterns controlling forest structure, which is needed for predicting long-term ecological change.

forest ecology↗

NASA NTRS1999

Ultrasound Thermal Field Imaging of Opaque Fluids

We have initiated an experimental program to develop an ultrasound system for non-intrusively imaging the thermal field in opaque fluids under an externally imposed temperature gradient. Many industrial processes involve opaque fluids, such as molten metals, semiconductors, and polymers, often in situations in which thermal gradients are important. For example, one may wish to understand semiconductor crystal growth dynamics in a Bridgman apparatus. Destructive testing of the crystal after the process is completed gives only indirect information about the fluid dynamics of the formation process. Knowledge of the coupled thermal and velocity fields during the growth process is then essential. Most techniques for non-intrusive velocity and temperature measurement in fluids are optical in nature, and hence the fluids studied must be transparent. In some cases (for example, LDV (laser Doppler velocimetry) and PIV (particle imaging velocimetry)) the velocities of small neutrally buoyant seed particles suspended in the fluid, are measured. Without particle seeding one can use the variation of the index of refraction of the fluid with temperature to visualize, through interferometric, Schlieren or shadowgraph techniques, the thermal field. The thermal field in turn gives a picture of the pattern existing in the fluid. If the object of study is opaque, non-optical techniques must be used. In this project we focus on the use of ultrasound, which propagates easily through opaque liquids and solids. To date ultrasound measurements have almost exclusively relied on the detection of sound scattered from density discontinuities inside the opaque material of interest. In most cases it has been used to visualize structural properties, but more recently the ultrasound Doppler velocimeter has become available. As in the optical case, it relies on seed particles that scatter Doppler shifted sound back to the detector. Doppler ultrasound techniques are, however, not useful for studying convective fluid flow in crystal growth, because particle seeding is unacceptable and flow velocities are typically too low to be resolved, and may be even lower in microgravity conditions where buoyancy forces are negligible. We will investigate a different use of ultrasound to probe the flows of opaque fluids non-intrusively and without the use of seed particles: our goal is to ultrasonically visualize the thermal field of opaque fluids with relatively high spatial resolution. The proposed technique relies upon the variation of sound speed with temperature of the fluid. A high frequency ultra-sound pulse passing through a fluid-filled chamber will traverse the chamber in a time determined by the relevant chamber dimension and the temperature of the fluid through which the pulse passes. With high time-resolution instrumentation that compares the excitation signal with the received pulse we can detect the influence of the fluid temperature on the pulse travel time. This is effectively an interferometric system, which in its optical form is an extremely sensitive approach to measuring thermal fields in fluids. Moreover, the temperature dependence of sound velocity in liquid metals is comparable to the temperature dependence of the speed of light required for accurate interferometric thermal images in transparent fluids. With an array of transducers scanned electronically a map of the thermal field over the chamber could be produced. An alternative approach would be to use the ultrasound analog of the shadowgraph technique. In the optical version, collimated light passes through the fluid, where it is focused or defocused locally by temperature field induced variations of the index of refraction. The resulting image reveals the thermal field through the spatial pattern of light intensity variations. By analogy, an ultrasound plane wave traversing an opaque fluid sample would be also locally focused or defocused depending on the speed of sound variations, giving rise to spatial variations in sound intensity that will reveal the thermal field pattern. These approaches could be applied to any situation in which temperature differences are expected to occur, and will rapidly provide information about the flow that simply cannot be obtained by any current intrusive or non-intrusive diagnostic technique. As materials processing in microgravity matures it will become increasingly important to have available simple and versatile diagnostic tools, such as we will develop, for studying the flows of opaque fluids under thermal forcing.

Andereck, C. David↗

DOE OSTI2025

Why the electron chooses one of two symmetry-related paths in the type II bacterial photosynthetic reaction centers

A series of electron and proton transfer reactions are carried out in photosynthetic reaction centers that convert the energy of light into high energy reduced products and add to the proton gradient. All photosynthetic reaction centers (photosystems I, II, and different bacterial reaction centers (bRCs)) have their charge separating cofactors arranged with two, c 2 -symmetric paths from a primary donor to a terminal acceptor. In Type II reaction centers (PSII and bRCs), electrons utilize only one active branch. In bRCs the electron transfer occurs in the A branch from the excited state of the bacteriochlorophyll (BChl) dimer P 860 through BChl A , bacteriopheophytin A (BPh A ), to ubiquinone Q A , and then to the B-branch Q B . B-branch BChl B and BPh B do not participate. A similar path is found in PSII. These proteins challenge our understanding of how proteins can turn on or turn off the activity of bound ligands. The shift of the electrochemical midpoint potential (E m ) of each cofactor by the protein environment is calculated using continuum electrostatics within the Multi-Conformation Continuum Electrostatic (MCCE) program. In bRCs from Rhodobacter sphaeroides, the electrostatic potentials favor reduction of the active branch cofactors. The residues that contribute to a more positive potential on the A branch are identified. These are often distant from the cofactors. The large difference of the potential at the competing BChls may direct the electron transfer to the A branch. Residues TyrM210 and PheL181 are at symmetrical positions on the A and B branches, near P 860 , between BChl A and BChl B . Mutated bRCs, that swap the Tyr and Phe (YF mutant), have been shown to support some electron transfer to the B-branch cofactors. Here, the calculated E m s with the YF swap show more positive B-branch potentials, which help explain why there is now electron transfer along the B branch.

Bacterial Reaction Center↗