Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “surface complexation model”

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 1,135 records · Page 63
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Transformation of an Embedded Five-Membered Ring in Polycyclic Aromatic Hydrocarbons via the Hydrogen-Abstraction–Acetylene-Addition Mechanism: A Theoretical Study

Five-membered rings are constituents of many polycyclic aromatic hydrocarbons (PAHs) and their presence on the edges of large PAHs has been repeatedly observed experimentally. However, modern kinetic combustion models often do not consider the growth of PAHs through the transformation of the five-membered rings. Here, in connection with the above, we carried out a theoretical study of the mechanism of Hydrogen-Abstraction-aCetylene-Addition (HACA) transformation of an embedded five-membered ring on the armchair PAH edge to a six-membered ring, considering cyclopenta[d,e,f]phenanthrene (4,5-methylenephenanthrene) as a prototype system for this process. The potential energy surface for the reactions of cyclopenta[d,e,f]phenanthrenyl radicals produced by direct H abstractions from cyclopenta[d,e,f]phenanthrene with acetylene has been compiled at the G3(MP2,CC)//B3LYP/6-311G(d,p) level of theory including zero-point vibrational energy corrections. The computed energies and molecular parameters were then used to solve the Rice-Ramsperger-Kassel-Marcus Master Equation (RRKM-ME) in order to calculate reaction rate at various pressures and temperatures, which were fitted to modified Arrhenius equation for further kinetic modeling. The results show that the HACA transformation of the embedded five-membered ring to a six-membered ring is possible, albeit slow. The most viable reaction mechanism involves the R2 + C 2 H 2 reaction, where the acetylene molecules adds to a σ-radical in the six-membered ring adjacent to the five-membered ring via a low entrance barrier. The predominant product of R2 + C 2 H 2 is predicted to be 3 ethynyl 4H cyclopenta[def]phenanthrene Pr5 via immediate H elimination from the initial addition complex. Next, Pr5 undergoes H-assisted isomerization to 4aH-pentaleno[4,3,2,1-cdef]phenanthrene Pr4, and the latter adds an H atom eventually forming the 1-pyrenylmethyl radical Pr3: R2 + C 2 H 2 ⇌ 3 ethynyl 4H cyclopenta[def]phenanthrene (Pr5) + H ⇌ 4aH-pentaleno[4,3,2,1-cdef]phenanthrene (Pr4) + H Pr5 + H ⇌ Pr4 + H Pr4 + H → 1-pyrenylmethyl (Pr3) This HACA sequence may be competitive with the methyl radical addition to the R1 radical formed by H abstraction from the CH 2 group in the five-membered ring of cyclopenta[d,e,f]phenanthrene, which provides a pathway to pyrene following two H atom losses. Relative contributions of the two mechanisms of the five- to six-membered ring transformation would strongly depend on the branching ratios of the R1 and R2 radicals produced by the H abstractions and the available concentration of C 2 H 2 vs. CH 3 and hence differ in different flames.

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

Variability in Ice Nucleating Particles Across Greater Houston Texas

The concentration and cloud-forming potential of a region's ice nucleating particle (INP) population have uncertain impacts on deep convective clouds. Specifically, ice nucleating particles (INPs) may affect various cloud properties related to the formation, lifetime, and precipitation of deep convective clouds. As part of the U.S. Department of Energy's TRacking Aerosol and Convection interaction ExpeRiment (TRACER) campaign, researchers from Texas A&M University deployed three Davis Rotating-drum Universal-size-cut Monitoring (DRUM) samplers throughout Greater Houston, Texas from June through September 2022. Ambient particles, collected at the surface with the DRUM samplers in four aerodynamic diameter size ranges (>3, 3–1.2, 1.2–0.34, and 0.34–0.15 μm), were analyzed in offline cold-stage ice nucleation experiments. The INP population in Greater Houston is complex, varying by site and day, but can be generalized by a weak to moderately efficient mode of INPs at −24°C and an efficient mode at −15°C. Analysis reveals that supermicron particles are largely responsible for ice nucleation warmer than −20°C across the region while submicron particles dominate at temperatures colder than −20°C. Additionally, significant spatial diversity in the INP population was observed, with differences in mean nucleation temperature between sites for nearly every size cut. Although INP concentrations were typically ∼0.08 L −1 at −20°C throughout the campaign, a notable region-wide increase in INP concentration for particles freezing at temperatures warmer than −20°C occurred from mid-August to mid-September. This comprehensive characterization of Greater Houston's INP population, including spatial, temporal, and particle size variations, can help constrain ice microphysics parameterizations in weather and climate models.

Thompson, Seth A. [Texas A & M Univ., College Stat↗
DOE OSTI2021

Which Combinations of Environmental Conditions and Microphysical Parameter Values Produce a Given Orographic Precipitation Distribution?

This study applies an idealized modeling framework, alongside a Bayesian Markov chain Monte Carlo (MCMC) algorithm, to explore which combinations of upstream environmental conditions and cloud microphysical parameter values can produce a particular precipitation distribution over an idealized two-dimensional, bell-shaped mountain. Simulations focus on orographic precipitation produced when an atmospheric river interacts with topography. MCMC-based analysis reveals that different combinations of parameter values produce a similar precipitation distribution, with the most influential parameters being relative humidity (RH), horizontal wind speed (U), surface potential temperature (θ sfc ), and the snow fall speed coefficient (A s ). RH, U, and A s exhibit interdependence: changes in one or more of these factors can be mitigated by compensating changes in the other(s) to produce similar orographic precipitation rates. The results also indicate that the parameter sensitivities and relationships can vary for spatial subregions and given different environmental conditions. In particular, high θ sfc values are more likely to produce the target precipitation rate and spatial distribution, and thus the ensemble of simulations shows a preference for liquid precipitation at the surface. The results presented here highlight the complexity of orographic precipitation controls, and have implications for flood and water management, observational efforts, and climate change.

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

Host-Directed, Bioelectronic Immunomodulation for Protection Against Emerging Pathogens

Acute care of patients with severe infections often relies on systemic administration of pharmaceuticals and monitoring of complex physiological symptoms to identify immune system dysfunction, which can lead to increased mortality. Furthermore, determining disease-specific treatment plans often leads to a delay in patient care. To address this, we proposed an immune modulation system that electrically detects and responds to a patient’s immune system status, creating an agnostic means of treating illness and infection. Two pieces of hardware were developed for this task: a minimally-invasive sensor and a vagus nerve stimulator. Stimulation of the vagus nerve is known to modulate the immune system. The sensor is a microfabricated, silicon-based microneedle array capable of interfacing with interstitial fluid to detect small molecules such as inflammatory proteins (cytokines) and pharmaceuticals (vancomycin). Process optimization to manufacture the needles refined the silicon etch process, creating needle patches long enough to penetrate skin and reach interstitial fluid. The needles were tested for mechanical strength and stability, and did not shatter when inserted into skin models. The needles are coated with a thin film metal, turning them into electrodes for electrochemical sensing of our target molecules. We hybridized aptamers to the surface of the electrode to act as the sensing layer and were able to detect changes in the conformation of the aptamer electrochemically in the presence of the target molecule. The stimulator was a cuff electrode that encircled the vagus nerve. Rodent studies were conducted in which rodents were exposed to an inflammatory event and vagus nerve stimulation (VNS) was applied. It was demonstrated that optimized electrical stimulation of the vagus nerve created measurably different levels of cytokines in blood samples, and certain cytokines released during the inflammatory event were either upregulated or downregulated. In sum, this project successfully developed new platforms and technologies that can, with further development, enable better temporal insight into biomarker changes in the body, letting healthcare providers know of possible immune system dysfunction before they are detected physiologically. We also demonstrated the value of VNS and its possible use in treating immune system response to inflammation and illness.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

The interplay between surface facet and reconstruction on isopropanol conversion over SrTiO 3 nanocrystals

Strontium titanate (SrTiO 3 ) is an extensively investigated perovskite for various applications due to its optical, electrical and chemical properties. To gain an in-depth understanding of the active sites involved in heterogeneous catalysis over the broadly used SrTiO 3 (STO), we studied a model reaction, isopropanol conversion, on three differently shape-controlled nanocrystals: cube, truncated cube and dodecahedra. SEM, XRD and XPS confirmed the morphology, phase and composition of STO shapes. Low energy ion scattering (LEIS) revealed the occurrence of surface reconstruction over STO shapes during O 2 pretreatment at different temperatures. Based on the catalytic activities, scanning transmission electron microscopy images and density functional theory calculations, the step sites on STO derived from surface reconstruction were proposed to be the active sites for isopropanol conversion. This was further confirmed by steady state isotopic kinetic analysis (SSITKA) which demonstrated similar intrinsic turnover frequencies (TOFs) for the differently reconstructed STO shapes. It is concluded that the crystal facets impose an indirect effect on the catalysis of STO via controlling the degrees of surface reconstruction: the less stable STO (1 1 0) facet (dodecahedra) leads to more step sites after reconstruction and hence higher overall reaction rate than the more stable (1 0 0) facet (cube). This work highlights the important interplay between the crystal facet and surface reconstruction in controlling the nature and density of active sties and thus catalysis over complex oxides.

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

Fast redox switches lead to rapid transformation of goethite in humid tropical soils: A Mössbauer spectroscopy study

Abstract Humid tropical forest soils experience frequent rainfall, which limits oxygen diffusion and creates redox heterogeneity in upland soils. In this study we gauged the effect of short‐term anoxic conditions on Fe mineralogy of relatively Fe‐ and C‐rich surface soils (C/Fe mole ratio ∼5) from a humid tropical forest in Puerto Rico. Soils subjected to 4‐d oxic/anoxic oscillation were characterized by selective chemical extractions, Mössbauer spectroscopy (MBS), and X‐ray diffraction. Chemical extraction data suggested that rapidly switching redox conditions had subtle effects on bulk Fe mineralogy. Mössbauer, on the other hand, indicated that (a) the soil Fe is a mixture of goethites of varying characteristics with minor contributions from ferrihydrite (<5%) and Fe(III)‐organic matter (OM) phases (∼10%), and (b) anoxic conditions rapidly transformed all forms of goethite to relatively stable Fe oxides. Such fast changes in goethite features could be due to rapid depletion and sorption of bio‐released structural Al or sorbed OM onto residual soil components. The rapid temporal changes in MBS parameters and corresponding pore water nominal oxidation state of C values suggest that Fe‐C transformations in these upland tropical soils are complex and intricately coupled. A comprehensive understanding of the fate of Fe, Al, and OM (Fe‐organic moieties) during redox switches and concurrent changes in pore water chemistry is critical for development of robust transport models in humid tropical soils, which are subject to episodic low‐redox events.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A novel active optimization approach for rapid and efficient design space exploration using ensemble machine learning

In this work, a novel design optimization technique based on active learning, which involves dynamic exploration and exploitation of the design space of interest using an ensemble of machine learning algorithms, is presented. In this approach, a hybrid methodology incorporating an explorative weak learner (regularized basis function model) which fits high-level information about the response surface, and an exploitative strong learner (based on committee machine) that fits finer details around promising regions identified by the weak learner, is employed. For each design iteration, an aristocratic approach is used to select a set of nominees, where points that meet a threshold merit value as predicted by the weak learner are selected to be evaluated using expensive function evaluation. In addition to these points, the global optimum as predicted by the strong learner is also evaluated to enable rapid convergence to the actual global optimum once the most promising region has been identified by the optimizer. This methodology is first tested by applying it to the optimization of a two-dimensional multi-modal surface. The performance of the new active learning approach is compared with traditional global optimization methods, namely micro-genetic algorithm (pGA) and particle swarm optimization (PSO). It is demonstrated that the new optimizer is able to reach the global optimum much faster with a significantly fewer number of function evaluations. Subsequently, the new optimizer is also applied to a complex internal combustion (IC) engine combustion optimization case with nine control parameters related to fuel injection, initial thermodynamic conditions, and in-cylinder flow. It is again found that the new approach significantly lowers the number of function evaluations that are needed to reach the optimum design configuration (by up to 80%) when compared to particle swarm and genetic algorithm-based optimization techniques.

Machine learning, active learning, design optimiza↗
DOE OSTI2025

Experiments to validate Thermodynamics and Transport models of Strongly Coupled Dusty Plasma Matter (Final Technical Report for DE-SC0023416)

The goal of this two-year grant is to provide access to the PI to dusty plasma experimental facilities at the DOE-funded Collaborative Research Facility Magnetized Plasma Research Laboratory, Auburn University to become a user of that facility, to obtain experimental data to support another ongoing grant DE-SC0021146 (an Early Career Award to the PI that is focused on modeling of dusty plasma thermodynamics and transport processes), to generate experimental data for funding proposals, and to provide exposure to University of Memphis students to advanced experimental techniques. The following technical accomplishments were made: 1. Development of a novel Bidirectional Electrode Control Arms Assembly (BECAA) for producing perfect 2D grain layers for complex plasma experimentation. BECAA uses movable electrode arms to tilt or move the electrode in a RF discharge from outside the chamber, allowing for the manipulation of grain clouds without needing to change the plasma parameters or gas pressure. This work addresses a longstanding gap in the literature for a method to produce clusters of selectable number of grains and that are perfectly two dimensional as opposed to being only quasi-2D. 2. Experimental investigation of the structural properties of finite-N clusters with N=2 to 50. Individual particle behavior in clusters could vary from grain to grain and this study measured systematically produced clusters for two different grain sizes. Analysis (funded by another grant DE-SC0021146) is currently underway to quantify the differences between grains that are found on the surface vs. the interior of clusters, the shell structure, and the decay of correlations in position, velocity, and kinetic energy. 3. An experimental method to measure the structural entropy of clusters was developed by observing the self-induced structural transitions between various possible arrangements. In a series of heating and cooling cycles, the number of times each possible arrangement was attained was experimentally observed and used to compute the probability of existence of that arrangement, and subsequently the configurational entropy of the cluster. Analysis (funded by another grant DE-SC0021146) is currently underway to produce the entropy of clusters as a function of the number of grains and use the same to compute thermodynamic state variables for 2D complex plasma/grain clusters. 4. A preliminary experimental study of multibody collisions between N grains (N=2 – 10) was conducted. The clusters were produced using the BECAA technique and velocities were imparted to the grains using manipulation laser pulses. Analysis (funded by another grant DE-SC0021146) is currently underway to develop a theoretical framework to describe multibody collisions analogous to classical two-body interactions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Improving energy yield in photovoltaic modules with photonic structures

When installed outdoors, Si solar cells typically operate 20 – 30 K above ambient conditions. These elevated operating temperatures lead to diminished performance, reducing efficiency by ~0.4%/K for crystalline Si cells and decreasing module lifetime. Since much of the elevated temperature of the modules derives from parasitic absorption of sub-bandgap sunlight by the contacts, encapsulants, or other materials, reflectors that remove the sub-bandgap radiation from the module before it is absorbed promise to reduce the operating temperature. However, there are many competing strategies to reject this light: selective reflectors could be integrated into the outer module glass, at the textured Si-encapsulant interface, or on the rear. Furthermore, many of these approaches are complex and expensive to fabricate. The issue of heating can be even more critical in bifacial modules: high fractions of rear irradiance can lead to higher operating temperatures. Prior to this project, relatively little was known about the optical strategies for temperature reduction in bifacial modules. The goal of this project was to quantitatively assess optical strategies for improving the energy yield via a combination of improved anti-reflection and sub-bandgap light rejection. We studied two complementary concepts: the performance of ideal structures to determine the upper limits to performance, and low-complexity structures made from standard materials as cost-effective implementations. We also studied the performance of these reflectors in different types of modules, including Al BSF, PERC, and bifacial. Here we build from our prior accomplishments in developing mirror optimization techniques, and make use of simulation methods that accurately capture the optical properties of the modules throughout the visible and near-infrared spectrum. We connect these optical properties to energy yield under realistic, outdoor conditions. We also experimentally fabricated mirrors, incorporated them into mini-modules, and performed outdoor testing of performance at NREL. We found definitively that spectrally-selective mirrors on the outside of the module cover glass offer the best performance, both under ideal conditions and for realistic, low- complexity mirrors. Mirrors at the Si cell/encapsulant interface are limited by the multiple reflections needed to remove sub-bandgap light. Mirrors at the cell surface and mirrors on the rear also cannot remove all the sub-bandgap light, since parasitic absorption also occurs in the cover glass, encapsulant, and front contacts, depending on the module configuration. Also, cell texturing places higher demands on any cell-surface mirror because multiple interactions with the mirror are required to reject sub bandgap light which compounds losses. For bifacial modules, mirrors on the back do provide additional cooling and energy yield improvements, and the mirrors can have the same structure. In the monofacial case up to 3.3K of cooling is possible for mirrors on the front, compared to 2.2K at the cell surface and 1.2K at the rear. For bifacial modules, up to 2.4K of cooling is possible when mirrors are included at both the front and back. Our realistic designs utilize 4 – 6 layers of common photovoltaic materials such as SiO 2 , TiO 2 , and SiN x , making them more cost effective. Finally, we used our models to study thermal management in bifacial modules with different back lamination materials, finding that glass-glass laminated modules operated hotter than equivalent glass-polymer backsheet modules, due solely to the increased absorption of light incident from the rear of the module. Improving the energy yield by 2%, which is achievable with these realistic designs, would lead to LCOE reductions of $0.03/kWh by 2030. These benefits make photovoltaic panels more accessible to the public and may reduce associated maintenance costs.

14 SOLAR ENERGY↗
DOE OSTI2022

Electrically controlling and optically observing the membrane potential of supported lipid bilayers

Supported lipid bilayers are a well-developed model system for the study of membranes and their associated proteins, such as membrane channels, enzymes, and receptors. These versatile model membranes can be made from various components, ranging from simple synthetic phospholipids to complex mixtures of constituents, mimicking the cell membrane with its relevant physiochemical and molecular phenomena. In addition, the high stability of supported lipid bilayers allows for their study via a wide array of experimental probes. In this work, we describe a platform for supported lipid bilayers that is accessible both electrically and optically, and demonstrate direct optical observation of the transmembrane potential of supported lipid bilayers. We show that the polarization of the supported membrane can be electrically controlled and optically probed using voltage-sensitive dyes. Membrane polarization dynamics is understood through electrochemical impedance spectroscopy and the analysis of an equivalent electrical circuit model. In addition, we describe the effect of the conducting electrode layer on the fluorescence of the optical probe through metal-induced energy transfer, and show that while this energy transfer has an adverse effect on the voltage sensitivity of the fluorescent probe, its strong distance dependency allows for axial localization of fluorescent emitters with ultrahigh accuracy. We conclude with a discussion on possible applications of this platform for the study of voltage-dependent membrane proteins and other processes in membrane biology and surface science.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Extracting single fiber transverse and shear moduli from off-axis misalignment fiber tensile testing

Small diameter (<100 μm) fibers (e.g. carbon fibers, Kevlar, and fiberglass) and wires (e.g. ultrafine copper and aluminum wires) are frequently used in many different engineering applications, such as for light weighting structures, electromagnetic shielding for aircraft/infrastructure/EVs, vibration damping, biological sensors, aerospace electrical devices, and electric windings just to name a few. Due to the manufacturing process, the fibers and wires are pulled and stretched to produce a preferential alignment. Therefore, thin fibers and wires typically display different properties along the length of the fiber as opposed to their cross section and many fibers/wires are considered transversely isotropic. The axial properties of fibers/wires can be ascertained via tensile testing of single-filaments or fiber tows, but the radial properties require much more effort to measure. Knowing these properties is important for the accurate prediction of micromechanical models and manipulation of fibers during micromanufacturing. In this paper, a new technique was developed to determine the transverse/shear moduli and strength of a material by conducting tensile tests of the material at increasing misalignment angles from the tensile axis. Due to the transversely isotropic nature of the material, the transverse/shear moduli and strength influence the experimental results recorded by the test machine to different degrees based on the amount of misalignment in the test setup. An equation was derived to determine the influence of each of the material properties based on the misalignment angle by manipulating the stiffness matrix for transversely isotropic materials using the transformation matrices. Then, curve fitted coefficients were used to identify the material properties. Here, the strengths were similarly determined by curve fitting an off-axis Tsai-Hill failure criteria to determine the influence of transverse, shear, and tensile strengths based on the complex loading condition provided by the off-axis tensile test. Zoltek Panex 35 carbon fibers were used to demonstrate this new technique and the determined properties were then compared to those obtained from nanoindentation and from literature. Fracture surfaces provide insight into the different failure mechanisms at various misalignment angles.

36 MATERIALS SCIENCE↗
DOE OSTI2024

FIU Project 2: Environmental Remediation Science & Technology [Slides]

FIU’s research under this project involves conducting basic and applied science to fill knowledge gaps and validate potential remediation technologies for contaminated soil and groundwater and the assessment of the fate and transport of contaminants in the environment. The aim of FIU’s research is to reduce the potential for contaminant mobility or toxicity in the surface and subsurface through the development and application of state-of-the-art scientific and environmental remediation technologies at the Hanford Site, Savannah River Site (SRS), and the Waste Isolation Pilot Plant (WIPP), which is the Nation’s only mined geologic repository for permanent disposal of transuranic waste. FIU collaborates with scientists from Pacific Northwest National Laboratory (PNNL), Savannah River National Laboratory (SRNL), Savannah River Ecology Laboratory (SREL), Los Alamos National Laboratory (LANL) and the DOE Carlsbad Field Office (CBFO) in order to plan and execute research that is synergistic with the work being conducted at the sites, and that supports the resolution of critical science and engineering needs which leads to a better understanding of the long-term behavior of subsurface contaminants. The knowledge gained through this research will be used to transform experimental and modeling innovations into practical applications deployed at the sites to support EM’s primary goal of expediting the closure of major contaminated soil and groundwater sites and waste units. Collaborative relationships between FIU and the national laboratories have provided large benefits over the years to FIU, the national laboratories, the DOE complex, and the DOE EM mission. By working closely with the national laboratories, FIU’s research is not only closely aligned with the cleanup mission priorities at the DOE sites, but complements and supports ongoing work at the national laboratories for screening of new remedial technologies. This coordination and leveraging of research efforts results in time- and cost-savings, and will accelerate progress of the DOE EM environmental restoration mission.

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

Combined ambient vibration and surface displacement measurements for improved progressive failure monitoring at a toppling rock slab in Utah, USA

Seismic resonance and surface displacement measurements can be implemented in tandem to improve landslide characterization and progressive failure monitoring. Crack aperture data are frequently used in rock slope stability monitoring and often exhibit recognizable trends prior to failure, such as accelerated crack opening. Alternatively, ambient resonance data offer multiple parameters including modal frequencies, damping, and polarization that can be monitored alongside crack aperture and may respond differently to environmental forcings and complex failure evolution. We analyzed data from continuous ambient vibration monitoring and concomitant crack aperture measurements at the Courthouse Mesa instability, a large toppling sandstone slab in Utah, USA. Three years of data revealed crack aperture increases of 2–4 mm/year with no clearly detectable irreversible changes in modal parameters, including frequency. Annually, frequency and displacement varied by 29% and 19% of the mean, respectively, with average and maximum daily frequency fluctuations of 6.5% and 16%, respectively. These reversible cyclic changes were primarily temperature-driven, but annually, frequency was in-phase with temperature whereas crack aperture lagged temperature changes by ~37 days. Polarization and damping also varied seasonally but were less strongly correlated with temperature. Conceptual 3D finite element modeling demonstrated consistent frequency decreases associated with crack propagation but variable changes in crack aperture measured at a single point; i.e., crack propagation did not always result in increased crack opening but always generated a resonance frequency decrease. Taken together, our data suggest a possible thermal wedging-ratcheting mechanism at the Courthouse Mesa instability, where annual thermoelastic crack closure is impeded by debris infill but the absence of downward crack propagation during the monitoring period is evidenced by no permanent resonance frequency changes. Our study demonstrates that combined seismic resonance and crack aperture data provide an improved description of rock slope instability behavior, supporting refined characterization and monitoring of changes accompanying progressive failure.

58 GEOSCIENCES↗
DOE OSTI2023

Boundary layer moisture variability at the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic observatory during marine conditions

Abstract. Boundary layer moisture variability at the Eastern North Atlantic (ENA) site during marine conditions is examined at monthly and daily timescales using 5 years of ground-based observations and output from the European Center for Medium range Weather Forecast (ECMWF) reanalysis model. The annual cycle of the mixed-layer total water budgets is presented to estimate the relative contribution of large-scale advection, local moisture tendency, entrainment, and precipitation to balance the moistening due to surface latent heat flux on monthly timescales. When marine conditions prevail, advection of colder and dry air from the north acts as an important moisture sink (∼ 50 % of the overall budget) during fall and winter driving the seasonality of the budget. Entrainment and precipitation contribute to the drying of the boundary layer (∼ 25 % and ∼ 15 % respectively), and the local change in moisture contributes to a small residual part. On a daily temporal scale, moist and dry mesoscale columns of vapor (∼ 10 km) are analyzed during 10 selected days of precipitating stratocumulus clouds. Adjacent moist and dry columns present distinct mesoscale features that are strongly correlated with clouds and precipitation. Dry columns adjacent to moist columns have more frequent and stronger downdrafts immediately below the cloud base. Moist columns have more frequent updrafts, stronger cloud-top cooling, and higher liquid water path and precipitation compared to the dry columns. This study highlights the complex interaction between large-scale and local processes controlling the boundary layer moisture and the importance of spatial distribution of vapor to support convection and precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effects of spatial resolution on WRF v3.8.1 simulated meteorology over the central Himalaya

The sensitive ecosystem of the central Himalayan (CH) region, which is experiencing enhanced stress from anthropogenic forcing, requires adequate atmospheric observations and an improved representation of the Himalaya in the models. However, the accuracy of atmospheric models remains limited in this region due to highly complex mountainous topography. This article delineates the effects of spatial resolution on the modeled meteorology and dynamics over the CH by utilizing the Weather Research and Forecasting (WRF) model extensively evaluated against the Ganges Valley Aerosol Experiment (GVAX) observations during the summer monsoon. The WRF simulation is performed over a domain (d01) encompassing northern India at 15 km x 15 km resolution and two nests (d02 at 5 km x 5 km and d03 at 1 km x 1 km) centered over the CH, with boundary conditions from the respective parent domains. WRF simulations reveal higher variability in meteorology, e.g., relative humidity (RH = 70.3%-96.1%) and wind speed (WS = 1.1-4.2m s(-1)), compared to the ERA-Interim reanalysis (RH D 80.0%-85.0%, WS D 1.2-2.3m s(-1)) over northern India owing to the higher resolution. WRF-simulated temporal evolution of meteorological variables is found to agree with balloon-borne measurements, with stronger correlations aloft (r = 0.44-0.92) than those in the lower troposphere (r = 0.18-0.48). The model overestimates temperature (warm bias by 2.8 degrees C) and underestimates RH (dry bias by 6.4 %) at the surface in d01. Model results show a significant improvement in d03 (P = 827.6 hPa, T = 19.8 degrees C, RH = 92.3 %), closer to the GVAX observations (P = 801.4 hPa, T = 19.5 degrees C, RH = 94.7 %). Interpolating the output from the coarser domains (d01, d02) to the altitude of the station reduces the biases in pressure and temperature; however, it suppresses the diurnal variations, highlighting the importance of well-resolved terrain (d03). Temporal variations in near-surface P, T, and RH are also reproduced by WRF in d03 to an extent (r>0:5). A sensitivity simulation incorporating the feedback from the nested domain demonstrates the improvement in simulated P, T, and RH over the CH. Our study shows that the WRF model setup at finer spatial resolution can significantly reduce the biases in simulated meteorology, and such an improved representation of the CH can be adopted through domain feedback into regional-scale simulations. Interestingly, WRF simulates a dominant easterly wind component at 1 km x 1 km resolution (d03), which is missing in the coarse simulations; however, the frequency of southeasterlies remains underestimated. The model simulation implementing a highresolution (3 s) topography input (SRTM) improved the prediction of wind directions; nevertheless, further improvements are required to better reproduce the observed local-scale dynamics over the CH.

Singh, Jaydeep↗
DOE OSTI2020

Bayesian Approach to Estimation of Water Table Elevations Using Historical Rasters as Prior Information 2019 - 20430

In cases of complex but only partially known geology and a lack of spatial control in observation well locations, water table elevation estimation is very challenging. In some cases, auxiliary information, such as observations of the movement of tracers, operation of injection and extraction wells, and calibration of groundwater models against the historical elevation data, can be combined with expert judgement to estimate flow directions in areas of sparse data and to aid in the production of more reliable contour maps (and associated rasters) than could be produced by relying on sparse well elevation data alone. Given a historical sequence of these raster maps, the question arises how to automate, to the extent possible, the process of producing new raster maps to reflect data from previous times, the current data and the operation of expert judgement. One solution is to adopt a Bayesian point of view and to regard the historical well elevation data, auxiliary information and historical raster maps as prior information. The well elevations for water table wells, as well as those for injection/extraction wells and the data associated with other relevant variables, can be viewed as predictors for the raster surface. From this prior information, we can, conditional on the values of these predictors for a new time period, compute an expected value map and a standard deviation map for the new raster. These then can be taken to specify a prior predictive distribution for the pixels in the new raster map. Then we condition the pixels, corresponding to water level observation wells within the raster, on the observed values in those wells (which in general will differ from the regression estimate) for the new time period. Given the smoothness of the water table surface, we then smooth the surface of deviations from the mean surface, based on the variograms of the historical rasters, and add this smoothed surface to the regression mean surface. The error structure of the produced raster map is defined by the regression error structure and the error due to smoothing based on the estimated variograms. This methodology has been developed and is being further refined for groundwater monitoring and remediation at LANL. It is a very flexible method that can also be applied with a variety of other predictors applied to model the water level wells in the area of interest over the historical record. The smoothness of the spatial process and its possible evolution over time can then be estimated from the residuals from this regression. This can be augmented by expert hydrogeological opinion based on site topography and hydrogeology. (authors)

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

FY25 MOOSE Usability Improvements: 3D Meshing Capabilities, Initiation of Geometry Support for Monte Carlo Tools, and Enhancement of MOOSE/Workbench User Input Interactions

Usability improvements have been made to MOOSE and Workbench in FY25 to enhance usability and user workflows. Assorted enhancement have been made to MOOSE’s intrinsic meshing capabilities in order to enable more flexible and complex meshing of nuclear reactor systems, in particular for 3D applications. Mesh generators have been added to perform operations such as batch mesh generation, surface mesh generation, and creation of 3D transition layers. These mesh generation capabilities make it much easier to generate high quality non-extruded 3D meshes. Additionally, work to integrate Monte Carlo reactor physics simulations into MOOSE-based multi-physics workflows has reached another milestone with the implementation of the Constructive Solid Geometry (CSG) base framework. This framework lays the foundation for mesh generators to offer the user a generic CSG output option (as opposed to a finite element mesh). To support users, workshop on the MOOSE Reactor Module was delivered which featured hands-on examples using the NEAMS Workbench on INL’s High Performance Computing system. Recent updates to the NEAMS Workbench, WASP, and the MOOSE language server have introduced several improvements aimed at making MOOSE-based simulation setup and input management faster, more accurate, and easier to use. Key capabilities that have been added include multi-tab-stop autocompletion, visual input diagnostics, developer-directed data visualizations, upgraded ParaView integration, and Workspace-level file tracking. Together, these changes make it easier for users to build, validate, and manage complex MOOSE-based simulation models — especially those involving reusable components, included files, and datasets. The improvements are designed to save time, reduce input errors, and help users get to a successful simulation run faster, with more confidence in the results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate

Characterizing the properties of groundwater aquifers is essential for predicting aquifer response and managing groundwater resources. In this work, we develop a high-dimensional scalable Bayesian inversion framework governed by a three-dimensional quasi-static linear poroelastic model to characterize lateral permeability variations in groundwater aquifers. In this work, we determine the maximum a posteriori (MAP) point of the posterior permeability distribution from centimeter-level surface deformation measurements obtained from Interferometric Synthetic Aperture Radar (InSAR). The scalability of our method to high parameter dimension is achieved through the use of adjoint-based derivatives, inexact Newton methods to determine the MAP point, and a Mat´ern class sparse prior precision operator. Together, these guarantee that the MAP point is found at a cost, measured in number of forward/adjoint poroelasticity solves, that is independent of the parameter dimension. We apply our methodology to a test case for a municipal well in Mesquite, Nevada, in which InSAR and GPS surface deformation data are available. We solve problems with up to 320,824 state variable degrees of freedom (DOFs) and 16,896 parameter DOFs. A consistent treatment of noise level is employed so that the aquifer characterization result does not depend on the pixel spacing of surface deformation data. Our results show that the use of InSAR data significantly improves characterization of lateral aquifer heterogeneity, and the InSAR-based aquifer characterization recovers complex lateral displacement trends observed by independent daily GPS measurements.

54 ENVIRONMENTAL SCIENCES↗