Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial process”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Machine Learning and Economic Models to Enable Risk-Informed Condition Based Maintenance of a Nuclear Plant Asset

The primary objective of this research is to address challenges in the implementation of risk-informed, condition-based predictive maintenance (PdM), which reduces operating costs while still maintaining the safety and reliability of commercial nuclear power plants (NPPs). To achieve the objective, risk models are being developed by taking advantage of advancements in data analytics, deep learning, machine learning (ML), and artificial intelligence (AI). The notable outcomes presented in the report include ? Development of a ML models using heterogeneous plant process and vibration data collected at different spatial and temporal resolutions from the Salem?s CWS to diagnose a circulating water pump (CWP) failure based on salient fault signatures. The developed diagnostic models are extendable to other faults associated with CWPs and CWP motors given associated fault signatures. ? Development of a natural language processing (NLP) technique to automatically classify the WO data into different categories. The developed NLP technique was validated on independent WO data. This automates the tedious and time-consuming activity of mining and classifying WOs by subject matter experts. ? Estimation of mean time between downtime (i.e., time duration between time instances when 1 or more CWPs are not available) and developed an approach to establish reliability of CWS components using unstructured WO data along with CWS plant process data. ? Formulation of economic model based on Markov chain models. The parameters of associated with the transition rate between different states of Markov chain models were estimated using WO data. The economic model formulation and discussion captures both time-independent and time-dependent parameter variation, leading to risk-informed decision-making.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

BioRT-Flux-PIHM v1.0: a biogeochemical reactive transport model at the watershed scale

Watersheds are the fundamental Earth surface functioning units that connect the land to aquatic systems. Many watershed-scale models represent hydrological processes but not biogeochemical reactive transport processes. This has limited our capability to understand and predict solute export, water chemistry and quality, and Earth system response to changing climate and anthropogenic conditions. Here we present a recently developed BioRT-Flux-PIHM (BioRT hereafter) v1.0, a watershed-scale biogeochemical reactive transport model. The model augments the previously developed RT-Flux-PIHM that integrates land-surface interactions, surface hydrology, and abiotic geochemical reactions. It enables the simulation of (1) shallow and deep-water partitioning to represent surface runoff, shallow soil water, and deeper groundwater and of (2) biotic processes including plant uptake, soil respiration, and nutrient transformation. The reactive transport part of the code has been verified against the widely used reactive transport code CrunchTope. BioRT-Flux-PIHM v1.0 has recently been applied in multiple watersheds under diverse climate, vegetation, and geological conditions. This paper briefly introduces the governing equations and model structure with a focus on new aspects of the model. It also showcases one hydrology example that simulates shallow and deep-water interactions and two biogeochemical examples relevant to nitrate and dissolved organic carbon (DOC). These examples are illustrated in two simulation modes of complexity. One is the spatially lumped mode (i.e., two land cells connected by one river segment) that focuses on processes and average behavior of a watershed. Another is the spatially distributed mode (i.e., hundreds of cells) that includes details of topography, land cover, and soil properties. Whereas the spatially lumped mode represents averaged properties and processes and temporal variations, the spatially distributed mode can be used to understand the impacts of spatial structure and identify hot spots of biogeochemical reactions. The model can be used to mechanistically understand coupled hydrological and biogeochemical processes under gradients of climate, vegetation, geology, and land use conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Development of the tangent linear and adjoint models of the global online chemical transport model MPAS-CO 2 v7.3

We describe the development of the tangent linear (TL) and adjoint models of the Model for Prediction Across Scales (MPAS)-CO 2 transport model, which is a global online chemical transport model developed upon the non-hydrostatic Model for Prediction Across Scales – Atmosphere (MPAS-A). The primary goal is to make the model system a valuable research tool for investigating atmospheric carbon transport and inverse modeling. First, we develop the TL code, encompassing all CO 2 transport processes within the MPAS-CO 2 forward model. Then, we construct the adjoint model using a combined strategy involving re-calculation and storage of the essential meteorological variables needed for CO 2 transport. This strategy allows the adjoint model to undertake a long-period integration with moderate memory demands. To ensure accuracy, the TL and adjoint models undergo vigorous verifications through a series of standard tests. The adjoint model, through backward-in-time integration, calculates the sensitivity of atmospheric CO 2 observations to surface CO 2 fluxes and the initial atmospheric CO 2 mixing ratio. To demonstrate the utility of the newly developed adjoint model, we conduct simulations for two types of atmospheric CO 2 observations, namely the tower-based in situ CO 2 mixing ratio and satellite-derived column-averaged CO 2 mixing ratio (X CO 2 ). A comparison between the sensitivity to surface flux calculated by the MPAS-CO 2 adjoint model with its counterpart from CarbonTracker–Lagrange (CT-L) reveals a spatial agreement but notable magnitude differences. These differences, particularly evident for X CO 2 , might be attributed to the two model systems' differences in the simulation configuration, spatial resolution, and treatment of vertical mixing processes. Moreover, this comparison highlights the substantial loss of information in the atmospheric CO 2 observations due to CT-L's spatial domain limitation. Furthermore, the adjoint sensitivity analysis demonstrates that the sensitivities to both surface flux and initial CO 2 conditions spread out throughout the entire Northern Hemisphere within a month. MPAS-CO 2 forward, TL, and adjoint models stand out for their calculation efficiency and variable-resolution capability, making them competitive in computational cost. In conclusion, the successful development of the MPAS-CO 2 TL and adjoint models, and their integration into the MPAS-CO 2 system, establish the possibility of using MPAS's unique features in atmospheric CO 2 transport sensitivity studies and in inverse modeling with advanced methods such as variational data assimilation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Using aircraft measurements to characterize subgrid-scale variability of aerosol properties near the Atmospheric Radiation Measurement Southern Great Plains site

Abstract. Complex distributions of aerosol properties evolve in space and time as a function of emissions, new particle formation, coagulation, condensational growth, chemical transformation, phase changes, turbulent mixing and transport, removal processes, and ambient meteorological conditions. The ability of chemical transport models to represent the multi-scale processes affecting the life cycle of aerosols depends on their spatial resolution since aerosol properties are assumed to be constant within a grid cell. Subgrid-scale-dependent processes that affect aerosol populations could have a significant impact on the formation of particles, their growth to cloud condensation nuclei (CCN) sizes, aerosol–cloud interactions, dry deposition and rainout and hence their burdens, lifetimes, and radiative forcing. To address this issue, we characterize subgrid-scale variability in terms of measured aerosol number, size, composition, hygroscopicity, and CCN concentrations made by repeated aircraft flight paths over the Atmospheric Radiation Measurement (ARM) program's Southern Great Plains (SGP) site during the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystem (HI-SCALE) campaign. Subgrid variability is quantified in terms of both normalized frequency distributions and percentage difference percentiles using grid spacings of 3, 9, 27, and 81 km that represent those typically used by cloud-system-resolving models as well as the current and next-generation climate models. Even though the SGP site is a rural location, surprisingly large horizontal gradients in aerosol properties were frequently observed. For example, 90 % of the 3, 9, and 27 km cell mean organic matter concentrations differed from the 81 km cell around the SGP site by as much as ∼ 46 %, large spatial variability in aerosol number concentrations and size distributions were found during new particle formation events, and consequently 90 % of the 3, 9, and 27 km cell mean CCN number concentrations differed from the 81 km cell mean by as much as ∼ 38 %. The spatial variability varied seasonally for some aerosol properties, with some having larger spatial variability during the spring and others having larger variability during the late summer. While measurements at a single surface site cannot reflect the surrounding variability of aerosol properties at a given time, aircraft measurements that are averaged within an 81 km cell were found to be similar to many, but not all, aerosol properties measured at the ground SGP site. This analysis suggests that it is reasonable to directly compare most ground SGP site aerosol measurements with coarse global climate model predictions. In addition, the variability quantified by the aircraft can be used as an uncertainty range when comparing the surface point measurements with model predictions that use coarse grid spacings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Derivation and assessment of regional electricity generation emission factors in the USA

Electricity production is one of the largest sources of environmental emissions - especially greenhouse gases (GHGs) - in the USA. Emission factors (EFs) vary from region to region, which requires the use of spatially relevant EF data for electricity production while performing life cycle assessments (LCAs). Uncertainty information, which is sought by LCA practitioners, is rarely supplied with available life cycle inventories (LCIs). To address these challenges, we present a method for collecting data from different sources for electricity generation and environmental emissions; discuss the challenges involved in agglomerating such data; provide relevant suggestions and solutions to merge the information; and calculate EFs for electricity generation processes from various fuel sources for different spatial regions and spatial resolutions. The EFs from the US 2016 Electricity Life Cycle Inventory (eLCI) are analyzed and explored in this study. We also explore the method of uncertainty information derivation for the EFs. We explore the EFs from different technologies across Emissions & Generation Resource Integrated Database (eGRID) regions in the USA. We find that for certain eGRID regions, the same electricity production technology may have worse emissions. This may be a result of the age of the plants in the region, the quality of fuel used, or other underlying factors. Region-wise life cycle impact assessment (LCIA) ISO 14040 impacts for total generation mix activities provide an overview of the total sustainability profile of electricity production in a particular region, rather than only global warming potential (GWP). We also find that, for different LCIA impacts, several eGRID regions are consistently worse than the US average LCIA impact for every unit of electricity generated. Here this work describes the development of an electricity production LCI at different spatial resolutions by combining and harmonizing information from several databases. The inventory consists of emissions, fuel inputs, and electricity and steam outputs from different electricity production technologies located across various regions of the USA. This LCI for electricity production in the USA will prove to be an enormous resource for all LCA researchers - considering the detailed sources of the information and the breadth of emissions covered by it.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Projecting Life-Cycle Environmental Impacts of Corn Production in the U.S. Midwest under Future Climate Scenarios using machine learning approach

Emission factors and process-based mechanism models are popular approaches used to estimate life-cycle environmental impacts. However, emission factors are often incapable of describing spatial and temporal heterogeneity of agricultural emissions, whereas process-based mechanism models, capable of capturing the heterogeneity, tend to be very complicated, and time-consuming. Efficient predictions of life-cycle environmental impacts from agricultural production are lacking. This study develops a rapid predictive model to quantify life-cycle global warming (GW), eutrophication (EU) and acidification (AD) impacts of corn production using a novel machine learning approach. We used the boosted regression tree (BRT) model to estimate current and future life-cycle environmental impacts of corn production in 12 U.S. Midwest at county scale under four Intergovernmental Panel on Climate Change (IPCC)’s climate scenarios (RCP2.6. 4.5, 6.0, and 8.5) for years 2009-2100. Results show that future life-cycle GW, EU and AD impacts of corn production will increase under all four scenarios, with the worst environmental impacts shown under the RCP 8.5 scenario. Significant spatial variation were estimated to be 3.6, 29.9 and 1.8-folds under the RCP4.5 scenario, and 6.6, 152.9 and 1.6-folds under the RCP8.5 scenario, respectively. Findings from this study demonstrate the importance of considering county-level life-cycle environmental impacts and the influences of climate change in future years to aid in developing potential adaptation and mitigation programs and policies.

Machine learning, climate change, projections, env↗
DOE OSTI2024

A prototype scintillator real‐time beam monitor for ultra‐high dose rate radiotherapy

Background: FLASH Radiotherapy (RT) is an emergent cancer RT modality where an entire therapeutic dose is delivered at more than 1000 times higher dose rate than conventional RT. For clinical trials to be conducted safely, a precise and fast beam monitor that can generate out-of-tolerance beam interrupts is required. This paper describes the overall concept and provides results from a prototype ultra-fast, scintillator-based beam monitor for both proton and electron beam FLASH applications. Purpose: A FLASH Beam Scintillator Monitor (FBSM) is being developed that employs a novel proprietary scintillator material. The FBSM has capabilities that conventional RT detector technologies are unable to simultaneously provide: (1) large area coverage; (2) a low mass profile; (3) a linear response over a broad dynamic range; (4) radiation hardness; (5) real-time analysis to provide an IEC-compliant fast beam-interrupt signal based on true two-dimensional beam imaging, radiation dosimetry and excellent spatial resolution. Methods: The FBSM uses a proprietary low mass, less than 0.5 mm water equivalent, non-hygroscopic, radiation tolerant scintillator material (designated HM: hybrid material) that is viewed by high frame rate CMOS cameras. Folded optics using mirrors enable a thin monitor profile of ∼10 cm. A field programmable gate array (FPGA) data acquisition system generates real-time analysis on a time scale appropriate to the FLASH RT beam modality: 100–1000 Hz for pulsed electrons and 10–20 kHz for quasi-continuous scanning proton pencil beams. An ion beam monitor served as the initial development platform for this work and was tested in low energy heavy-ion beams ( 86 Kr +26 and protons). A prototype FBSM was fabricated and then tested in various radiation beams that included FLASH level dose per pulse electron beams, and a hospital RT clinic with electron beams. Results: Results presented in this report include image quality, response linearity, radiation hardness, spatial resolution, and real-time data processing. Furthermore, the HM scintillator was found to be highly radiation damage resistant. It exhibited a small 0.025%/kGy signal decrease from a 216 kGy cumulative dose resulting from continuous exposure for 15 min at a FLASH compatible dose rate of 237 Gy/s. Measurements of the signal amplitude versus beam fluence demonstrate linear response of the FBSM at FLASH compatible dose rates of >40 Gy/s. Comparison with commercial Gafchromic film indicates that the FBSM produces a high resolution 2D beam image and can reproduce a nearly identical beam profile, including primary beam tails. The spatial resolution was measured at 35–40 µm. Tests of the firmware beta version show successful operation at 20 000 Hz frame rate or 50 µs/frame, where the real-time analysis of the beam parameters is achieved in less than 1 µs. Conclusions: The FBSM is designed to provide real-time beam profile monitoring over a large active area without significantly degrading the beam quality. A prototype device has been staged in particle beams at currents of single particles up to FLASH level dose rates, using both continuous ion beams and pulsed electron beams. Using a novel scintillator, beam profiling has been demonstrated for currents extending from single particles to 10 nA currents. Radiation damage is minimal and even under FLASH conditions would require ≥50 kGy of accumulated exposure in a single spot to result in a 1% decrease in signal output. Beam imaging is comparable to radiochromic films, and provides immediate images without hours of processing. Real-time data processing, taking less than 50 µs (combined data transfer and analysis times), has been implemented in firmware for 20 kHz frame rates for continuous proton beams.

2D beam imaging↗
DOE OSTI2021

Chemical vapor infiltration of additively manufactured preforms: Pore-resolved simulations and experimental validation

The densification of additively manufactured porous preforms by chemical vapor infiltration (CVI) is studied using pore-resolved simulations and experiments. Experimentally, 3D printed silicon carbide (SiC) preforms are subject to CVI synthesis using methyltrichlorosilane (MTS) precursor to obtain high purity SiC/SiC composites. Optical images of the cross sections of the processed preforms are analyzed to obtain the spatial porosity distribution. The numerical method is based on a level set formulation to capture the spatial distribution and time evolution of the pore scale microstructural characteristics. The coupled transport and kinetic effects are represented using a dimensionless Thiele modulus. Simulations are initialized using representative synthetic preform geometries comprising of packed particles based on the size distribution of the powder used for 3D printing. The simulation results are validated against the experimental observations in terms of total density and the distribution of residual porosity. The densification characteristics, porosity classification, concentration profiles, and structure functions are analyzed as functions of processing temperature and Thiele modulus.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Investigating the Effects of Mixing Dynamics on Twin-Screw Granule Quality Attributes via the Development of a Physics-Based Process Map

Twin-screw granulation (TSG) is an emerging continuous wet granulation technique that has not been widely applied in the industry due to a poor mechanistic understanding of the process. This study focuses on improving this mechanistic understanding by analyzing the effects of the mixing dynamics on the granule quality attributes (PSD, content uniformity, and microstructure). Mixing is an important dynamic process that simultaneously occurs along with the granulation rate mechanisms during the wet granulation process. An improved mechanistic understanding was achieved by identifying and quantifying the physically relevant intermediate parameters that affect the mixing dynamics in TSG, and then their effects on the granule attributes were analyzed by investigating their effects on the granulation rate mechanisms. The fill level, granule liquid saturation, extent of nucleation, and powder wettability were found to be the key physically relevant intermediate parameters that affect the mixing inside the twin-screw granulator. An improved geometrical model for the fill level was developed and validated against existing experimental data. Finally, a process map was developed to depict the effects of mixing on the temporal and spatial evolution of the materials inside the twin-screw granulator. This process map illustrates the mechanism of nucleation and the growth of the granules based on the fundamental material properties of the primary powders (solubility and wettability), liquid binders (viscosity), and mixing dynamics present in the system. Furthermore, it was shown that the process map can be used to predict the granule product quality based on the granule growth mechanism.

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

Performance of Compact Pulsed Thermal Imaging System for In-Service Applications. Pulsed thermal tomography nondestructive examination of additively manufactured reactor materials and components

Additive manufacturing (AM) is an emerging method for cost-efficient fabrication of complex topology nuclear reactor parts from high-strength corrosion resistance alloys, such as stainless steel and Inconel. AM of metallic structures for nuclear energy applications is currently based on laser powder bed fusion (LPBF) process, which has the capability of melting metallic powder and net shaping the structures with relatively high precision. Some of the challenges with using LPBF method for nuclear manufacturing include the possibility of introducing pores into metallic structures. Integrity of AM structures needs to be evaluated nondestructively because material flaws could lead to premature failures due to creep in high temperature nuclear reactor environment. Currently, there exist limited capabilities to evaluate actual AM structures nondestructively. Pulsed Thermography (PT) imaging provides a capability for non-destructive evaluation (NDE) of sub-surface defects in arbitrary size structures. The PT method is based on recording material surface temperature transients with infrared (IR) camera following thermal pulse delivered on material surface with flash light. The PT method has advantages for NDE of actual AM structures because the method involves one-sided non-contact measurements and fast processing of large sample areas captured in one image. The data cube of PT measurements consists of surface temperature taken at sequential time intervals T(x,y,t). Material defects can be detected either by analyzing the thermograms T(x,y,t) data cube, or by using thermal tomography (TT) algorithm to obtain 3D spatial reconstruction of thermal effusivity e(x,y,z). To reduce the cost and enable in-service NDE in spatially constrained environment, it is highly desirable to develop PT with compact and inexpensive IR camera. Following initial qualification of an AM component for deployment in a nuclear reactor, a compact PT system can also be used for in-service nondestructive evaluation (NDE) applications. However, data cube obtained with PT based on compact IR camera suffers from strong thermal noises and loss of features due to relatively low sampling rate. In this report we describe two unsupervised machine learning (ML) algorithms for enhancement of PT images obtained with compact IR camera. In one approach, we introduce Sparse Coding Discrete Cosine Transform (SC/DCT) algorithm to remove additive white Gaussian noise (AWGN) from spatial thermal effusivity reconstructions. In another approach we introduce a Spatial Temporal Denoised Thermal Source Separation (STDTSS) ML algorithm to process thermograms. The STDTSS algorithm consists of spatial and temporal denoising using Gaussian and Savitzky–Golay filtering, followed by the matrix decomposition using Principal Component Analysis (PCA), and Independent Component Analysis (ICA) to automatically detect flaws. In the work described in this report, we constructed a compact PT system using a relatively small and low-cost FLIR A65 camera, consisting on uncooled microbolometer detector. Performance of SC/DCT algorithm was demonstrated on enhancing TT images of Inconel 718 AM plate. Performance of the STDTSS methods was investigated using thermography data obtained from imaging stainless steel 316L specimens produced with LPBF method with imprinted calibrated porosity defects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Optimizing quantification of MK6240 tau PET in unimpaired older adults

Accurate measurement of Alzheimer's disease (AD) pathology in older adults without significant clinical impairment is critical to assessing intervention strategies aimed at slowing AD-related cognitive decline. The U.S. Study to Protect Brain Health Through Lifestyle Intervention to Reduce Risk (POINTER) is a 2-year randomized controlled trial to evaluate the effect of multicomponent risk reduction strategies in older adults (60-79 years) who are cognitively unimpaired but at increased risk for cognitive decline/dementia due to factors such as cardiovascular disease and family history. The POINTER Imaging ancillary study is collecting tau-PET ([ 18 F]MK6240), beta-amyloid (Aβ)-PET ([ 18 F]florbetaben [FBB]) and MRI data to evaluate neuroimaging biomarkers of AD and cerebrovascular pathophysiology in this at-risk sample. Here 481 participants (70.0±5.0; 66% F) with baseline MK6240, FBB and structural MRI scans were included. PET scans were coregistered to the structural MRI which was used to create FreeSurfer-defined reference regions and target regions of interest (ROIs). We also created off-target signal (OTS) ROIs to examine the magnitude and distribution of MK6240 OTS across the brain as well as relationships between OTS and age, sex, and race. OTS was unimodally distributed, highly correlated across OTS ROIs and related to younger age and sex but not race. Aiming to identify an optimal processing approach for MK6240 that would reduce the influence of OTS, we compared our previously validated MRI-guided standard PET processing and 6 alternative approaches. The alternate approaches included combinations of reference region erosion and meningeal OTS masking before spatial smoothing as well as partial volume correction. To compare processing approaches we examined relationships between target ROIs (entorhinal cortex (ERC), hippocampus or a temporal meta-ROI (MetaROI)) SUVR and age, sex, race, Aβ and a general cognitive status measure, the Modified Telephone Interview for Cognitive Status (TICSm). Overall, the processing approaches performed similarly, and none showed a meaningful improvement over standard processing. Across processing approaches we observed previously reported relationships with MK6240 target ROIs including positive associations with age, an Aβ+> Aβ- effect and negative associations with cognition. In sum, we demonstrated that different methods for minimizing effects of OTS, which is highly correlated across the brain within subject, produced no substantive change in our performance metrics. This is likely because OTS contaminates both reference and target regions and this contamination largely cancels out in SUVR data. Caution should be used when efforts to reduce OTS focus on target or reference regions in isolation as this may exacerbate OTS contamination in SUVR data.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Mapping Glacier Basal Sliding Applying Machine Learning

During the RESOLVE project (“High-resolution imaging in subsurface geophysics: development of a multi-instrument platform for interdisciplinary research”), continuous surface displacement and seismic array observations were obtained on Glacier d’Argentière in the French Alps for 35 days in May 2018. The data set is used to perform a detailed study of targeted processes within the highly dynamic cryospheric environment. In particular, the physical processes controlling glacial basal motion are poorly understood and remain challenging to observe directly. Especially in the Alpine region for temperate based glaciers where the ice rapidly responds to changing climatic conditions and thus, processes are strongly intermittent in time and heterogeneous in space. Spatially dense seismic and Global Positioning System (GPS) measurements are analyzed applying machine learning to gain insight into the processes controlling glacial motions of Glacier d’Argentière. Using multiple bandpass-filtered copies of the continuous seismic waveforms, we compute energy-based features, develop a matched field beamforming catalog and include meteorological observations. Features describing the data are analyzed with a gradient boosting decision tree model to directly estimate the GPS displacements from the seismic noise. We posit that features of the seismic noise provide direct access to the dominant parameters that drive displacement on the highly variable and unsteady surface of the glacier. The machine learning model infers daily fluctuations and longer term trends. The results show on-ice displacement rates are strongly modulated by activity at the base of the glacier. The techniques presented provide a new approach to study glacial basal sliding and discover its full complexity.

58 GEOSCIENCES↗
DOE OSTI2022

Characterizing soil water content variability across spatial scales from optimized high-resolution distributed temperature sensing technique

Fiber-optic Distributed Temperature Sensing, when combined with the Single-probe Heat-pulse technique can measure soil moisture (θ) across spatial scales. The key limitation of this system is in obtaining the relationship between soil thermal conductivity (λ) and θ for a specific field. Using the Department of Energy Atmospheric Radiation Measurement (ARM) site, this study tested a new methodology to account for the spatial variability in the λ-θ relationship using a Gaussian processes model. The resulting accurate θ measurements (RMSE = 0.03 m 3 m –3 ) were used to characterize the spatial variability of θ across scales and to develop an empirical equation that can correct for the changes in the θ spatial variability observed at different spatial resolutions. In addition, the number of required samples to accurately characterize θ and its variability over scales ranging from 5 m and 350 m were estimated. Finally, these findings provide key information to scale soil moisture from centimeters to hundreds of meters for process understanding.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

System and method for monitoring atomic absorption during a surface modification process

A processing system monitors and/or controls a surface modification process occurring on a substrate within a processing chamber. An optical processing module having a light emission submodule to output a generated light signal and an optical detection submodule to detect a resultant light signal, is connected via fiber optic cables to light illuminating and light receiving components located within the chamber. A processor determines an amount of atomic absorption by an atomic element encountered by a probing beam passing between the illuminating and receiving components, based on the intensity of the generated light signal, the intensity of the received light signal and optionally the spontaneous emission of the atomic element in the absence of illumination by a probing beam. Based on the determined amount, the system derives a plurality of parameters of the modified substrate, their spatial and temporal uniformity, and information about process conditions in the processing chamber.

Atanasoff, George↗
DOE OSTI2022

A New Model for Simulating the Imbibition of a Wetting-Phase Fluid in a Matrix-Fracture Dual Connectivity System

The imbibition experiment is an effective approach for measuring petrophysical properties of porous media, with many such experiments performed over the past decade. Quite some empirical, analytical, and numerical models have been developed to simulate spontaneous imbibition of the wetting phase fluid into porous media, but limitations still exist. In previous studies, the imbibition process has been considered to give a piston-like displacement or the porous medium modeled as multiply-sized pores linked with bonds; both approaches fail to yield comprehensive results due to their neglect of the presence of irregular fractures or nonuniform flow paths through the matrix. By building a numerical model for simulating laboratory-scale experimental data, we performed imbibition tests on several fractured Barnett Shale samples having fractures either parallel ( P ) or transverse ( T ) to the bedding plane and used MATLAB to build a new numerical model by combining the imbibition process in fractures and the matrix using concepts from percolation theory. The experimental data show that the rocks with P -direction fractures have a more steady increase of imbibition rates than the case of T -direction one. As the shale matrix with low pore connectivity hampers the upward water movement, the imbibition rate of shales with T -direction fractures will decrease suddenly after the bottom layer in contact with water is saturated during the initial period. This wetting phase movement (WPM) model can simulate 3D porous media with 2D fractures. The rate of imbibition by fractured porous media is associated with physical parameters such as porosity and fracture distribution (e.g., the number and angle of fractures). Using Monte Carlo methods, we examined fracture parameters and predicted elapsed time and cumulative water imbibition, for the Barnett Shale samples. The results show that the rate of imbibed water mass is sensitive to the number of fractures directly connected to water source, and the connectivity between two neighboring grid cells is a key parameter for the wetting-front progression. The findings of this study can help to better understand the imbibition process with multiple influencing processes and factors in fractured-matrix rocks. Although the experiments, data simulation, and prediction results are based only on Barnett Shale samples, the model is readily applicable to imbibition tests of other fractured rocks to show the spatial and temporal behavior during a dynamic imbibition process that are not easily captured experimentally.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Soil moisture and hydrology projections of the permafrost region – a model intercomparison

Abstract. This study investigates and compares soil moisture and hydrology projections of broadly used land models with permafrost processes and highlights the causes and impacts of permafrost zone soil moisture projections. Climate models project warmer temperatures and increases in precipitation (P) which will intensify evapotranspiration (ET) and runoff in land models. However, this study shows that most models project a long-term drying of the surface soil (0–20 cm) for the permafrost region despite increases in the net air–surface water flux (P-ET). Drying is generally explained by infiltration of moisture to deeper soil layers as the active layer deepens or permafrost thaws completely. Although most models agree on drying, the projections vary strongly in magnitude and spatial pattern. Land models tend to agree with decadal runoff trends but underestimate runoff volume when compared to gauge data across the major Arctic river basins, potentially indicating model structural limitations. Coordinated efforts to address the ongoing challenges presented in this study will help reduce uncertainty in our capability to predict the future Arctic hydrological state and associated land–atmosphere biogeochemical processes across spatial and temporal scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Noise suppression for MHD characterization with electron cyclotron emission imaging 1D technique

Significant noise suppression for magnetohydrodynamics (MHD) mode characterization in the spatial and spectral domain is achieved by processing two-dimensional (2D) electron cyclotron emission imaging (ECEI) data with a one-dimensional (1D) ECEI technique using a short time window ($1\;{\text{ms}}$). The technique is applied to detect toroidal Alfven eigenmodes (TAEs) in the temporal spectrum and fit their radial envelope using the data from the DIII-D tokamak W-band 2D ECEI system. Using the data length (time window) of only 1 ${\text{ms}}$, the 1D ECEI can clearly detect the TAEs (~100 ${\text{kHz}}$) on the spectrum, while similar spectrum quality requires ~10 ${\text{ms}}$ data length with the cross power spectrum between two midplane ECEI channels. The 1D ECEI technique also effectively avoids biased fitting when resolving the fine structure of the TAE's radial envelope. The radially spatial resolution of 1D ECEI is constrained by the finite ECE radiation volume of the ECEI receiver. Here, with forward radiation modeling, we find the DIII-D ECEI system can sensitively measure the even parity MHD activities, for which the mode width is >15 mm, and tearing modes (odd parity MHD activities), for which the island full width is >30 mm.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗