Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical dispersion”

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

Multi-Scale Integrated Monitoring System for Enhancing Methane Emission Detection, Quantification & Prediction

This report details the progress and findings of a comprehensive study on reviewing existing solutions, identifying technology gaps, and formulating an “all-in-one” integrated strategy for developing the next-generation multiscale methane monitoring and modeling platform, conducted under grant number DE-FE0032292. Co-led by Dr. David Ebert, Dr. Binbin Weng, and Dr. Chenghao Wang at the University of Oklahoma, the project’s goal was to develop an integrated approach for building this engineering platform to detect, quantify, and mitigate methane emissions across various temporal scale, spatial scales, and sectors. The planning grant study began with an extensive review of various methane sensing and monitoring technologies and systems, surveying over 100 technology providers globally. This review revealed the prevalence of optical methods over chemical methods in commercially available sensors, with Non-Dispersive Infrared (NDIR), Tunable Diode Laser Absorption Spectroscopy (TDLAS), and Optical Gas Imaging (OGI) cameras being the most prevalent options. A trend towards more advanced optical techniques was observed, driven by increased regulatory focus and technological advancements. The technical evaluation of these sensing technologies provided crucial insights into their capabilities and limitations. The study examined emerging technologies such as Differential Absorption LiDAR (DIAL), which show promise for high-precision and long-range detection. The team then investigated the features and application bandwidth of various sensing platforms, including handheld, fixed/stationary, mobile, aerials, and spaceborne monitors. Pilot field studies were conducted to assess the capabilities of solutions for different emission scenarios. Field work with sensor deployments was conducted at three distinct site types: an oil & gas industry site, a cattle ranching operation, and a waste processing facility. The team also conducted a thorough review of methane flux inverse modeling approaches, focused on physically based methods. These approaches were categorized into simple, intermediate, and advanced methods. A realtime WRF-GHG (Weather Research and Forecasting-Greenhouse Gas) modeling system was developed and applied, incorporating multiple data sources to guide field experiments and inform methane plume detection. The project identified and analyzed numerous categories of methane data sources, including satellite measurements, ground-based sensors, and inventory databases. Key platforms examined include EDGAR, EPA GHGI, NASA TROPOMI, Carbon Mapper, and Climate TRACE, among others. The team proposed an architecture for a comprehensive methane monitoring platform. This system incorporates multi-source data acquisition, advanced data processing and assimilation, interactive visualization tools, and analytical capabilities for emissions forecasting and scenario analysis. The proposed platform aims to provide a user-friendly interface catering to various stakeholders, from researchers to policymakers. The architecture includes sophisticated data ingestion methods, a centralized data warehouse, and advanced analytical tools for data fusion and interpretation. To ensure the relevance and effectiveness of the proposed system, a comprehensive survey was conducted to gather stakeholder input on system requirements. Key findings include a strong need for integrating various data types and formats, a preference for real-time data updates and advanced visualization tools, and a demand for user-friendly interfaces catering to different expertise levels.

03 NATURAL GAS↗
DOE OSTI2024

Surface-active antibiotic production as a multifunctional adaptation for postfire microorganisms

Abstract Wildfires affect soils in multiple ways, leading to numerous challenges for colonizing microorganisms. Although it is thought that fire-adapted microorganisms lie at the forefront of postfire ecosystem recovery, the specific strategies that these organisms use to thrive in burned soils remain largely unknown. Through bioactivity screening of bacterial isolates from burned soils, we discovered that several Paraburkholderia spp. isolates produced a set of unusual rhamnolipid surfactants with a natural methyl ester modification. These rhamnolipid methyl esters (RLMEs) exhibited enhanced antimicrobial activity against other postfire microbial isolates, including pyrophilous Pyronema fungi and Amycolatopsis bacteria, compared to the typical rhamnolipids made by organisms such as Pseudomonas spp. RLMEs also showed enhanced surfactant properties and facilitated bacterial motility on agar surfaces. In vitro assays further demonstrated that RLMEs improved aqueous solubilization of polycyclic aromatic hydrocarbons, which are potential carbon sources found in char. Identification of the rhamnolipid biosynthesis genes in the postfire isolate, Paraburkholderia kirstenboschensis str. F3, led to the discovery of rhlM, whose gene product is responsible for the unique methylation of rhamnolipid substrates. RhlM is the first characterized bacterial representative of a large class of integral membrane methyltransferases that are widespread in bacteria. These results indicate multiple roles for RLMEs in the postfire lifestyle of Paraburkholderia isolates, including enhanced dispersal, solubilization of potential nutrients, and inhibition of competitors. Our findings shed new light on the chemical adaptations that bacteria employ to navigate, grow, and outcompete other soil community members in postfire environments.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Growth of Jupiter: Formation in disks of gas and solids and evolution to the present epoch

The formation of Jupiter is modeled via core-nucleated accretion, and the planet's evolution is simulated up to the present epoch. Throughout the phases when the planet acquires most of its heavy-element content, the calculation of solids' accretion accounts for interactions with an evolving disk of planetesimals. The phase of growth from an embryo of a few hundred kilometers in radius until the time when the accretion of gas overtakes solids' accretion was presented by D'Angelo et al., and the same numerical methods are applied here. Those calculations followed the formation for about 4 × 10 5 years, until the epoch when the heavy-element and hydrogen/helium masses were $M_Z ≈ 7.3$ and $M_{XY} ≈ 0.15$ Earth's masses ($M_⊕$), respectively, and $\dot {M} _{XY} ≈ $$\dot {M}_Z$. In this work, the calculation is continued through the phase when $M_{XY}$ grows to equal $M_Z$ , at which age, about 2.4 × 10 6 years, the total mass of the planet is $M_p ≈ 20 M_⊕$. About 9 × 10 5 years later, $M_p$ is approximately $60 M_⊕$ and $M_Z ≈ 16 M_⊕$, three-quarters of which are delivered by planetesimals larger than 10 km in radius. Around this epoch, the contraction of the envelope dictates gas accretion rates a few times 10 –3 $M_⊕$ per year, initiating the regime of disk-limited accretion, whereby the planet can accrete all the gas provided by the disk, and its evolution is therefore tied to disk's evolution. Growth is continued by constructing simplified models of protosolar accretion disks that evolve through viscous diffusion, winds, and accretion on the planet. Jupiter's formation ends after ≈ 3.4–4.2 Myr, depending on the applied disk viscosity parameter, when nebula gas disperses. The young Jupiter is 4.5–5.5 times as voluminous as it is presently and thousands of times as luminous, ~10 –5 $L_⊙$. The heavy-element mass is ≈ 20 $M_⊕$. The evolution proceeds through the cooling and contraction phase, in isolation except for solar irradiation. After 4570 Myr, the age of the solar system, radius and luminosity of the planet are within 10$\%$ of current values, accounting also for uncertainties in the power absorbed from the Sun. During formation, and soon thereafter, the planet exhibits features, e.g., luminosity and effective temperature, that may probe aspects of the latter stages of formation, if observable. These possibly distinctive features, however, seem to disappear within a few tens of Myr.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Probing the hydrolytic degradation of UF 4 in humid air

This manuscript describes the chemical transformations that occur during hydrolysis of uranium tetrafluoride (UF 4 ) due to its storage in humid air (85% and 50% relative humidity) at ambient temperatures. This hydrolysis was previously reported to proceed slowly or not at all (depending on the percent relative humidity); however, previous reports relied primarily on X-ray diffraction methods to probe uranium speciation. Here, we employ a battery of physiochemical probing techniques to explore potential hydrolysis, including Raman spectroscopy, powder X-ray diffraction, 19 F nuclear magnetic resonance spectroscopy, scanning electron microscopy, and focused ion beam microscopy with energy-dispersive X-ray spectroscopy. Of these, only Raman spectroscopy proved to be particularly useful at observing chemical changes to UF 4 . It was found that anhydrous UF 4 slightly oxidizes over the course of thirteen days to Schoepite-like uranium complexes and possibly UO 3 . In contrast, UF 4 exposed to 50% relative humidity slightly decomposes into UO 2 F 2 , Schoepite-like uranium complexes, and possibly a high order uranium oxide that eluded chemical assignment (U x O y ). Despite the rich chemical speciation observed in our Raman spectroscopy measurements, X-ray diffraction and 19 F NMR measurements on the same material showed no changes. Microscopy measurements suggest that the observed reactions between UF 4 and water occur primarily on the surface of UF 4 particulates via a method that is visually similar to surface corrosion of metals. Therefore, we postulate that NMR spectroscopy and X-ray diffraction, which are well-suited for bulk analysis, are less suited than Raman spectroscopy to observe the surface-based reactions that occur to UF 4 when exposed to humid air. Considering the importance of UF 4 in the production of nuclear fuel and weapons, the results presented herein are widely applicable to numerous nuclear science fields where uranium detection and speciation in humid environments is of value, including nuclear nonproliferation and nuclear forensics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Generation of selective single-mode guided waves by d 36 type piezoelectric wafer

In general, mixed modes and dispersion of guided wave occur at any excitation frequency, while single modes are needed for effective nondestructive evaluation of structures. In this work, we present an approach to generate a selective single-mode guided wave in plate-like structures by exploiting the unique directionality of the d 36 type piezoelectric wafer and the symmetry of fundamental guided wave modes. Specifically, we devise a unified fundamental shear horizontal (SH 0 ) wave/fundamental antisymmetric mode (A 0 ) wave directional transducer pair by attaching the d 36 type piezoelectric wafers back to back on the opposite sides of the plate to provide selective, directional single-mode guided waves for actuation and sensing. Originally, the d 36 -type piezoelectric wafer, poled and cut from the lead magnesium niobate-lead titanate crystal, produces a mixture of the in-plane SH 0 (symmetric) mode and out-of-plane A 0 (asymmetric) and S 0 (symmetric) modes. By applying the in-phase or out-of-phase applied electrical field to the coupled d 36 -type piezoelectric wafer pair, our devised approach generates selective single modes, SH 0 or A 0 , respectively, to propagate in the structure. We describe the theoretical development of our approach and conduct both numerical simulations and laboratory experiments for validations.

42 ENGINEERING↗
DOE OSTI2026

Influence of Nanoconfinement and Elevated Temperatures on Geocolloid-Facilitated Transport of Energy-related Contaminants

The main objective of this project is to obtain a better understanding of the transport behavior and interactions of geocolloids in the presence of energy- related contaminants under bulk and nanoconfined conditions. In that, we first established solid protocols to synthesize and fabricate numerous types of geocolloids with a. various sizes ranging from 40 – 800 nm and b. different hydrophile-lipophile balance (HLB) ratios in the range of 25:75 to 75:25. Direct force measurements with geocolloids having different degrees of surface coverages (i.e. HLB to realistically mimic adsorption of energy related contaminants) were conducted using the Surface Forces Apparatus (SFA) over a distance regime starting from 8μm all the way down to molecular contact. Repulsive forces were observed on approach starting from > 3μm, followed by an exponential increase of which magnitude appears to be larger than a decay length obtained from Derjaguin–Landau–Verwey–Overbeek (DLVO) theory in pure water. When the geocolloids were confined in salted water, the magnitude of onset of repulsion was varied as a function of salinity in solution, which can significantly alter a purely repulsive screened electrostatic (coulombic) interaction arising from, among geocolloids as well as between geocolloids and geosurfaces. The viscosity and flow characteristics of geocolloidal suspensions at different degrees of confinement were also investigated where we identified highly discontinuous rheological behaviors below a critical nanoconfinement level. We anticipate that the knowledge gained through this study will enable the scientists and researchers to better assess transport and fate behaviors of geocolloidal dispersions that can carry energy-related contaminants under realistically emulated geosystem.

03 NATURAL GAS↗
DOE OSTI2022

Modeling multiple scattering transient of an ultrashort laser pulse by spherical particles

The multiple scattering of an ultrashort laser pulse by a turbid dispersive medium (namely a cloud of bubbles in water) is investigated by means of Monte Carlo simulations. The theory of Gouesbet and Gréhan (2000) is used to derive an energetic model of the scattering transient. It is shown that the spreading and extinction of the pulse can be decoupled from the transient of scattering, which allows to describe each phenomenon individually. The transient of scattering is modeled with the Lorenz-Mie Theory and thus is also valid for a relative refractive index lower than one, contrary to the Debye series expansion which does not converge close to the critical angle. This is made possible after the introduction of a new physical object, the Scattering Impulse Response Function (SIRF) which allows to detect the different modes of scattering transient, in time and direction. Here the present approach is more generic, as it enables to simulate clouds of air bubbles in water, which was not possible previously. Two different approaches are proposed within the Monte Carlo framework. The first is a pure Monte Carlo approach where the delay due to the scattering is randomly drawn at each event, while the second is based on the transport of the whole scattering signal. They are both embedded in the Monte Carlo code Scatter3D. Both models produce equivalent trends and are validated against published numerical results. They are then applied to the multiple scattering of ultra short pulse by a cloud of bubble in water in the forward direction. The pulse spread due to the propagation in water is computed for a wide range of traveled distances and pulse durations, and the optimal pulse duration is given to minimize the pulse spread at a given distance. The main result is that the scattered photons exit the turbid medium earlier than the ballistic photons and produce a double peak related to the refraction in the bubble. This demonstrates the possibility to develop new diagnostics to characterize dynamic bubbly flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Bandwidth extension and conversion efficiency improvements beyond phase matching limitations using cavity-enhanced OPCPA

The conversion efficiency and phase matching bandwidth of ultrafast optical parametric amplification (OPA) are constrained by the dispersion and nonlinear coefficient of the employed crystal as well as pulse shaping effects. In our work we show that an enhancement cavity resonant with the pump seeded at the full repetition rate of the pump laser can automatically reshape the small-signal gain in optical parametric chirped-pulse amplification (OPCPA) to achieve close-to-optimal operation. This new method termed cavity-enhanced OPCPA or C-OPCPA significantly increases both the gain bandwidth and the conversion efficiency, in addition to boosting gain for high-repetition-rate amplification. The goal in C-OPCPA is to arrive at a condition of impedance matching at all temporal coordinates, such that, in the absence of linear losses, all the incident pump power is dissipated in the nonlinear loss element, i.e., converted to signal and idler. The use of a low finesse enhancement cavity resonant with a low average power (<1W) and a high repetition rate (78MHz) pump source is shown to achieve more than 50% conversion efficiency into signal and idler from the coupled pump in an optical parametric process, whereas an equivalent amount of pump power in a single-pass configuration leads to negligible conversion. Additionally, the gain bandwidth is extended by a factor of 3-4 beyond the phase-matching limit. Our empirical observations are corroborated by a numerical analysis of depletion optimizing the single-pass case, which assesses the underlying impedance matching that is responsible for the observed performance improvements.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Assessing complexity and dynamics in epidemics: geographical barriers and facilitators of foot-and-mouth disease dissemination

Introduction: Physical and non-physical processes that occur in nature may influence biological processes, such as dissemination of infectious diseases. However, such processes may be hard to detect when they are complex systems. Because complexity is a dynamic and non-linear interaction among numerous elements and structural levels in which specific effects are not necessarily linked to any one specific element, cause-effect connections are rarely or poorly observed. Methods: To test this hypothesis, the complex and dynamic properties of geo-biological data were explored with high-resolution epidemiological data collected in the 2001 Uruguayan foot-and-mouth disease (FMD) epizootic that mainly affected cattle. County-level data on cases, farm density, road density, river density, and the ratio of road (or river) length/county perimeter were analyzed with an open-ended procedure that identified geographical clustering in the first 11 epidemic weeks. Two questions were asked: (i) do geo-referenced epidemiologic data display complex properties? and (ii) can such properties facilitate or prevent disease dissemination? Results: Emergent patterns were detected when complex data structures were analyzed, which were not observed when variables were assessed individually. Complex properties–including data circularity–were demonstrated. The emergent patterns helped identify 11 counties as ‘disseminators’ or ‘facilitators’ (F) and 264 counties as ‘barriers’ (B) of epidemic spread. In the early epidemic phase, F and B counties differed in terms of road density and FMD case density. Focusing on non-biological, geographical data, a second analysis indicated that complex relationships may identify B-like counties even before epidemics occur. Discussion: Geographical barriers and/or promoters of disease dispersal may precede the introduction of emerging pathogens. If corroborated, the analysis of geo-referenced complexity may support anticipatory epidemiological policies.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Magnetic-Confinement Fusion—Plasma Theory: Tokamak Magnetohydrodynamic Equilibrium and Stability

Magnetohydrodynamics (MHD) provides a useful model to describe the crucial plasma macroscopic equilibrium and stability behaviors in toroidal tokamak devices by considering the plasma as a conducting fluid interacting with a surrounding confining electromagnetic field. MHD is the most basic plasma model, incorporating most large-scale phenomena, including plasma equilibrium and all major instabilities. MHD equations are obtained by taking moments of the Boltzmann equations for different plasma species. They provide a set of comprehensive physics constrains to compute and optimize the equilibrium plasma shape and pressure and current profiles that are critical to its stability and performance. In the ideal case, the equations have special properties that lead to efficient numerical calculation schemes, the most important of which is the ideal MHD energy principle for linear stability against small departures from equilibrium. In a tokamak plasma, equilibrium pressure is mostly destabilizing for MHD modes, whereas equilibrium current is also often a major driving force. Plasma resistivity creates new freedom for a MHD instability to grow, but there are also cases where the plasma resistivity plays a stabilizing role. Equilibrium toroidal flow and/or flow shear can affect MHD instabilities. Principal MHD instabilities include the internal kink mode, sawtooth, fishbone, external kink, resistive wall mode, resistive interchange, tearing and neoclassical tearing modes (NTMs), locked modes, toroidal Alfven eigenmodes (TAEs), and edge localized modes (ELMs). Fast-growing MHD instabilities can lead to an abrupt plasma disruption and termination that can potentially damage the device plasma facing components (PFCs) and in-vessel structures. Furthermore, an important MHD application is to develop robust techniques to mitigate and control MHD instabilities.

Dispersive pellet injection↗
DOE OSTI2023

Machine learning methods for particle stress development in suspension Poiseuille flows

Numerical simulations are used to study the dynamics of a developing suspension Poiseuille flow with monodispersed and bidispersed neutrally buoyant particles in a planar channel, and machine learning is applied to learn the evolving stresses of the developing suspension. The particle stresses and pressure develop on a slower time scale than the volume fraction, indicating that once the particles reach a steady volume fraction profile, they rearrange to minimize the contact pressure on each particle. Here we consider how the stress development leads to particle migration, time scales for stress development, and present a new physics-informed Galerkin neural network that allows for learning the particle stresses when direct measurements are not possible. The particle fluxes are compared with the Suspension Balance Model with good agreement. We show that when stress measurements are possible, the MOR-physics operator learning method can also capture the particle stresses.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗
DOE OSTI2023

Machine-to-machine variability of roughness and corrosion in additively manufactured 316L stainless steel

Numerous studies on the corrosion response of metal AM have been conducted. Nonetheless, the specimens being tested are commonly ground or polished to remove the outer as-built surface. If metal AM is truly going to be employed as a transformative technology that can produce complex shapes that do not require traditional machining, then the material needs to be evaluated in the as-built state. The reality is that AM alloys have shown significant inconsistencies regarding as-built surface texture, topology, and residual stress. One metric that has shown significant unpredictability in the literature is the susceptibility to localized corrosion of typically passive alloys, such as stainless steel (SS). There are a large number of studies that have attempted to understand the corrosion response of metal AM materials, but these studies are typically performed on materials that have been printed on a single machine and often mechanically polished to a smooth finish. This study compares the corrosion response of as-built AM, laser-beam powder bed fusion (LB-PBF), 316L SS parts that have been fabricated on five different machines. The majority of this work focused on understanding the susceptibility to localized corrosion of AM metals with respect to machine-dependent variables, namely surface roughness, and build angle. Surface roughness data was collected using scanning white light triangulation, laser scanning confocal microscopy, and coherence scanning interferometry. The results show that there is significant variability (p two-tail < 0.05) in the susceptibility to local corrosion initiation of LB-PBF 316L SS samples built on different machines. Surface oxides were probed with electron dispersive spectroscopy and revealed that variations in local corrosion susceptibility likely arise from differences in the stability of the passive film caused by chemical segregation, unique microstructure, and tortuous roughness features at the as-built surfaces. The variability of roughness and corrosion properties from test samples printed on different machines was corroborated by property measurements performed at five different testing sites, proving reproducibility of the data. Most importantly, this study shows that if the as-built surface layer is removed through grinding or electropolishing the machine-to-machine variation observed in the corrosion susceptibility is reduced.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Understanding the effect of density functional choice and van der Waals treatment on predicting the binding configuration, loading, and stability of amine-grafted metal organic frameworks

Metal organic frameworks (MOFs) are crystalline, three-dimensional structures with high surface areas and tunable porosities. Made from metal nodes connected by organic linkers, the exact properties of a given MOF are determined by node and linker choice. MOFs hold promise for numerous applications, including gas capture and storage. M2(4,4′-dioxidobiphenyl-3,3′-dicarboxylate)—henceforth simply M2(dobpdc), with M = Mg, Mn, Fe, Co, Ni, Cu, or Zn—is regarded as one of the most promising structures for CO2 capture applications. Further modification of the MOF with diamines or tetramines can significantly boost gas species selectivity, a necessity for the ultra-dilute CO2 concentrations in the direct-air capture of CO2. There are countless potential diamines and tetramines, paving the way for a vast number of potential sorbents to be probed for CO2 adsorption properties. The number of amines and their configuration in the MOF pore are key drivers of CO2 adsorption capacity and kinetics, and so a validation of computational prediction of these quantities is required to suitably use computational methods in the discovery and screening of amine-functionalized sorbents. In this work, we study the predictive accuracy of density functional theory and related calculations on amine loading and configuration for one diamine and two tetramines. In particular, we explore the Perdew–Burke–Ernzerhof (PBE) functional and its formulation for solids (PBEsol) with and without the Grimme-D2 and Grimme-D3 pairwise corrections (PBE+D2/3 and PBEsol+D2/3), two revised PBE functionals with the Grimme-D2 and Grimme-D3 pairwise corrections (RPBE+D2/3 and revPBE+D2/3), and the nonlocal van der Waals correlation (vdW-DF2) functional. We also investigate a universal graph deep learning interatomic potential’s (M3GNet) predictive accuracy for loading and configuration. These results allow us to identify a useful screening procedure for configuration prediction that has a coarse component for quick evaluation and a higher accuracy component for detailed analysis. Our general observation is that the neural network-based potential can be used as a high-level and rapid screening tool, whereas PBEsol+D3 gives a completely qualitatively predictive picture across all systems studied, and can thus be used for high accuracy motif predictions. We close by briefly exploring the predictions of relative thermal stability for the different functionals and dispersion corrections.

Chemistry↗
DOE OSTI2021

Spatial Decorrelation of Young Stars and Dense Gas as a Probe of the Star Formation–Feedback Cycle in Galaxies

The spatial decorrelation of dense molecular gas and young stars observed on ≲1 kpc scales in nearby galaxies indicates rapid dispersal of star-forming regions by stellar feedback. In this work, we explore the sensitivity of this decorrelation to different processes controlling the structure of the interstellar medium, the abundance of molecular gas, star formation, and feedback in a suite of simulations of an isolated dwarf galaxy with structural properties similar to NGC 300 that self-consistently model radiative transfer and molecular chemistry. Our fiducial simulation reproduces the magnitude of decorrelation and its scale dependence measured in NGC 300, and we show that this agreement is due to different aspects of feedback, including H 2 dissociation, gas heating by the locally variable UV field, early mechanical feedback, and supernovae. In particular, early radiative and mechanical feedback affects the correlation on ≲100 pc scales, while supernovae play a significant role on ≳ 100 pc scales. The correlation is also sensitive to the choice of the local star formation efficiency per free fall time, ϵ ff , which provides a strong observational constraint on ϵ ff when the global star formation rate is independent of its value. Finally, we explicitly show that the degree of correlation between the peaks of molecular gas and star formation density is directly related to the distribution of the lifetimes of star-forming regions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Influence of Soil Properties on Sea-Breeze Circulations in the Southeast U.S.

Sea-breeze circulations (SBCs) are common weather phenomena at and near coastal regions. They form because of a thermal gradient between the land surface at the coast and the sea surface. In a mid-day regime, a “thermal low” generated at the warm coast will lead to rising air motion, creating a wind shift coming from the sea near the surface displacing the coastal air. A “return flow” moving back towards the sea is generated by upper-level divergence because of the rising motion from the thermal low. SBCs propagate and serve as a method of urban pollutant dispersion in the Los Angeles region of California and are constrained to the coast due to the topography of surrounding mountains serving as a boundary for further inland propagation. Within the northeast U.S. SBCs are seen in the warm season but tend to remain coastally bound due to Coriolis distortion over long distances. Within the southeast U.S. (SEUS), the paradigmatic example of SBCs occurs over the Florida peninsula, where thunderstorms form on a nearly daily occurrence due to the convergence of SBCs from the east and west sides of the peninsula. However, there are further examples of sea-breezes in the SEUS that warrant study. Within the region bordering the SEUS and the Mid-Atlantic, just east of the southern Appalachian Mountains, warm-season SBCs form at the coast of Georgia and the Carolinas. Relatively flat topography ~150-200km inland allows for mostly unimpeded inland SBC propagation. Through visual analysis, Viner et al. catalogued several SBCs that propagated as far inland as the Central Savannah River Area surrounding Augusta, Georgia. Wermter et al. found that while the land-sea thermal gradient at the coast can influence coastal SBC genesis, the inland-coastal thermal gradient over the land is the primary influencer on the speed and depth of inland propagation of SBCs in this region. Additionally, soil moisture itself is a known correlative factor to sea-breeze formation, as it influences the soil temperature and the thermal gradient needed for SBC genesis and inland penetration. Physick determined that higher latent heat fluxes associated with wetter soil dampen the land-sea thermal gradient and suppress the formation of a SBC. Physick determined that higher latent heat fluxes associated with wetter soil dampen the land-sea thermal gradient and suppress the formation of a SBC. Conversely, drier soil enhances the thermal gradient and promotes SBC formation. However, while there is an inverse relationship between soil moisture and SBCs, higher soil moisture can actually promote more convective rainfall following a SBC if it does not significantly impact the thermal gradient. While the relationship between soil moisture and SBC formation has been conceptually explored and modeled numerically, there is a research gap in observed connections. The Soil Moisture Active Passive (SMAP) satellite mission has been operational since 2015 and has been used to create high-resolution re-analytical Level 4 (L4) datasets of soil moisture and soil temperature at different soil depths: the surface (0-5cm) and rootzone (0-1m). The surface soil temperature effectively acts as the “skin temperature” of the surface at these levels, and a spatial map of the land-sea as well as the coastal-inland thermal gradients can be represented. SMAP data are also assimilated in some atmospheric models such at the High Resolution Rapid Refresh (HRRR) mesoscale model. We propose leveraging the use of SMAP products to fill in spatial gaps left by weather and mesonet stations within the SEUS region, as well as assessing the effectiveness of utilizing SMAP products towards SBC forecasting in both deterministic and machine learning (ML) models.

58 GEOSCIENCES↗
DOE OSTI2020

Mesoscale to Microscale Simulations over Complex Terrain with the Immersed Boundary Method in the Weather Research and Forecasting Model

Improvements to the Weather Research and Forecasting (WRF) Model are made to enable multiscale simulations over highly complex terrain with dynamically downscaled boundary conditions from the mesoscale to the microscale. Over steep terrain, the WRF Model develops numerical errors that are due to grid deformation of the terrain-following coordinates. An alternative coordinate system, the immersed boundary method (IBM), has been implemented into WRF, allowing for simulations over highly complex terrain; however, the new coordinate system precluded nesting within mesoscale simulations using WRF’s native terrain-following coordinates. Here, the immersed boundary method and WRF’s grid-nesting framework are modified to seamlessly work together. This improved framework for the first time allows for large-eddy simulation over complex (urban) terrain with IBM to be nested within a typical mesoscale WRF simulation. Simulations of the Joint Urban 2003 field campaign in Oklahoma City, Oklahoma, are performed using a multiscale five-domain nested configuration, spanning horizontal grid resolutions from 6 km to 2 m. These are compared with microscale-only simulations with idealized lateral boundary conditions and with observations of wind speed/direction and SF 6 concentrations from a controlled release from intensive observation period 3. The multiscale simulation, which is configured independent of local observations, shows similar model skill predicting wind speed/direction and improved skill predicting SF 6 concentrations when compared with the idealized simulations, which require use of observations to set mean flow conditions. Use of this improved multiscale framework shows promise for enabling large-eddy simulation over highly complex terrain with dynamically downscaled boundary conditions from mesoscale models.

42 ENGINEERING↗
DOE OSTI2020

Tracking Radioactive Isotopes in HVAC and Application for Hot Cell Analyses - 20155

In recent years, the surge in the number of isotope production facilities under design has increased the need to analyze isotope migration in facility Heating, Ventilation, and Air Conditioning (HVAC) during normal operations and accident scenarios. These facilities are used to produce isotopes for medical, security, and industrial applications and as such are subject to license and regulatory requirements 10CFR20, 10CFR30, 10CFR50, and 10CFR70. GOTHIC, a general-purpose thermal-hydraulics software package, includes the ability to model isotopic tracers, radioactive decay and isotope migration as well as HEPA and charcoal filters for isotope retention. A GOTHIC model was developed by Zachry Nuclear Engineering (ZNE) to examine the effect of negative room pressure, HEPA filtration, and HVAC fluctuations on radiation areas, hot cells and gloveboxes. Radioactive tracers were used to simulate the concentration of spills within contaminated areas, track the migration of isotopes of interest, and determine the isotopic retention and buildup on facility HEPA filters. A variety of isotopes with concern to dose (e.g., Kr-85, Sr-90, I-131, etc.), including their decay and progeny, are included in the analysis. Negative pressures are maintained in the regions of interest by a representative central HVAC system equipped with a volumetric fan that exhausts to the environment after a series of isolation valves and HEPA filters. GOTHIC is an industry trusted tool for providing engineering solutions for a variety of applications, including fission product tracking, aerosol and particulate transport and ventilation assessments. The software provides an integrated analysis environment that includes a graphical user interface (GUI) for constructing analysis models, a numerical solver that includes parallel processing capabilities and a post-processor for evaluating simulation results. It solves the conservation equations for mass, momentum and energy for multicomponent, multi-phase flow in lumped parameter and multi-dimensional geometries (1, 2, or full 3D), including the effects of turbulence, diffusion and buoyancy. It has been developed and maintained under a Quality Assurance program in compliance with the requirements of 10CFR50 Appendix B and applicable portions of ASME NQA-1 since 1995. GOTHIC has been used for assessing both forced and natural convection conditions for a wide range of applications, including: - Tracking concentration of hazardous gases and chemicals for habitability and safety assessments - Determining ventilation and filtration requirements and optimizing location and arrangement of these systems - Room heat-up, including diverse and Flexible coping strategies for Extended Loss of AC Power (FLEX/ELAP) - Equipment Qualification (EQ) A distinctive feature of GOTHIC is the ability to track many different fields/substances in a simulation, including user defined tracer elements, in the liquid, vapor and droplet fields as well as surfaces and filters. This capability allows GOTHIC to model fission product transport and release or the removal of particulates or harmful toxins from exhaust gases using a spray scrubber or other types of filtration systems. GOTHIC also includes models for engineered equipment, such as fans, filters, charcoal filters, dryers/demisters, dampers, etc. The aerosols and other filtered material are removed or accumulated in these components. The range of aerosol and radiological applications that GOTHIC has been used for includes: - Source Term: Primary Coolant (Equilibrium) Activity; Non-Water Coolant Source/Leakage. - Conditions for Iodine Re-evolution: Sump/Suppression Pool Conditions and pH; - RWST Conditions and pH. - Isotope Removal Mechanisms: Containment Sprayed and Unsprayed Region Mixing; Charcoal Filter Heating due to Iodine decay. - Radionuclide Transport and Decay: Post-LOCA Release in containment; Transport between connected Compartments and vent systems; Groundwater transport of radionuclides. - Non-Newtonian Fluid modeling for sludge, waste tanks, etc. ADAMS ML071581053 (titled 'Best Practice Guidelines for the use of CFD in Nuclear Reactor Safety Applications') poses guidelines for applying single phase CFD codes in nuclear reactor safety problems and GOTHIC is listed as a 'tool for 3D flows' and 'dispersal and deposition of radionuclides.' The Nuclear Quality Assurance (NQA) pedigree of GOTHIC is an important aspect for applications in the nuclear industry. The fundamental tracer models (convective transport, molecular and turbulent diffusion, removal mechanisms, etc.) have been verified using analytical solutions and validated against applicable separate effects tests. Also, GOTHIC has been benchmarked to many integrated effects tests, including Phebus FP (Fission Product). GOTHIC gives good agreement for the buildup and decay of fission products in Phebus Test 3. The model developed by ZNE demonstrates GOTHIC's applicability and acceptability for use in analyzing the migration and retention of radioactive isotopes and their progeny in normal operation and accident scenario analyses for isotope production facilities, hot cells, and gloveboxes. The tracer activities calculated by GOTHIC can then be used in downstream radiation transport and shielding codes like RADTRAD-NAI{sup C}, MCNP{sup R}, and MicroShield{sup R} to determine on-site and of-site doses. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗