Search NASA⌕ Search

SEARCH · Search NASA

Results for “spatial prediction”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Observation of the molecular response to light upon photoexcitation

Abstract When a molecule interacts with light, its electrons can absorb energy from the electromagnetic field by rapidly rearranging their positions. This constitutes the first step of photochemical and photophysical processes that include primary events in human vision and photosynthesis. Here, we report the direct measurement of the initial redistribution of electron density when the molecule 1,3-cyclohexadiene (CHD) is optically excited. Our experiments exploit the intense, ultrashort hard x-ray pulses of the Linac Coherent Light Source (LCLS) to map the change in electron density using ultrafast x-ray scattering. The nature of the excited electronic state is identified with excellent spatial resolution and in good agreement with theoretical predictions. The excited state electron density distributions are thus amenable to direct experimental observation.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

A renewable double plasma mirror for Petawatt-class lasers

Exceptional pulse contrast can be critical for ultraintense laser experiments, particularly when using solid density targets, and their use is becoming widespread. However, current plasma mirror technology is becoming inadequate for the new generation of high repetition rate, high power lasers now available. We describe a novel double plasma mirror configuration based on renewable, free standing, ultrathin liquid crystal films tested at the BELLA Petawatt Laser Center. Although operating at a repetition rate of several shots per minute, this system can be scaled to a high repetition rate exceeding 1 Hz and represents an important step towards enabling sustained, continuous operation of plasma mirrors. We demonstrate an improvement of two to three orders of magnitude in contrast and a total throughput of 80%. We present the first measurements of a beam reflected from a single or double plasma mirror system using a wavefront sensor, showing a well preserved wavefront and spatial mode. Finally, we introduce a model that predicts the total throughput through this double plasma mirror. This is the first model that accurately predicts the peak reflectivity of a plasma mirror when given the laser temporal profile.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Topological Circular Dichroism in Chiral Multifold Semimetals

Uncovering the physical contents of the nontrivial topology of quantum states is a critical problem in condensed matter physics. Here, we study the topological circular dichroism in chiral semimetals using linear response theory and first-principles calculations. We show that, when the low-energy spectrum respects emergent SO(3) rotational symmetry, topological circular dichroism is forbidden for Weyl fermions, and thus is unique to chiral multifold fermions. This is a result of the selection rule that is imposed by the emergent symmetry under the combination of particle-hole conjugation and spatial inversion. Using first-principles calculations, we predict that topological circular dichroism occurs in CoSi for photon energy below about 0.2 eV. In conclusion, our Letter demonstrates the existence of a response property of unconventional fermions that is fundamentally different from the response of Dirac and Weyl fermions, motivating further study to uncover other unique responses.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Estimating CO 2 fluxes through integrating spatial and temporal input layers via deep learning algorithms

Background Accurate estimation of net ecosystem exchange of CO 2 fluxes (Fc) is essential for understanding carbon cycle processes and assessing ecosystem carbon budgets. However, conventional modeling approaches often emphasize temporal dynamics while overlooking the pronounced spatial heterogeneity within the footprint of eddy covariance (EC) towers, potentially limiting predictive accuracy and interpretability of Fc estimates. To address this challenge, we developed a spatiotemporal model that integrates high-resolution footprint-weighted spatial information with sequential environmental drivers. Results The integrated model combines a deeper graph convolutional network to characterize fine-scale spatial variability within EC footprints and a gated recurrent unit network to capture temporal dependencies in biophysical conditions. Using multi-year flux tower observations, remote sensing vegetation indices and footprint modeling, we evaluate the proposed method across three land cover types. This spatiotemporal model consistently outperforms temporal-only and spatial-only baselines, achieving the highest overall accuracy (R 2 = 0.9569) and the lowest RMSE (1.8128 μmol m −2 s −1 ) and MAE (1.1939 μmol m −2 s −1 ). Performance gains are particularly evident in ecosystems with strong vegetation heterogeneity, where spatial structure substantially modulates Fc variability. Conclusions This study demonstrates the importance of joint modeling spatial heterogeneity and temporal dynamics for improving Fc estimation and provides a robust method for advancing footprint-based Fc estimates across diverse ecosystems, supporting refined assessments of terrestrial carbon fluxes, and enhancing scientific foundations for carbon studies.

CO2 flux estimate↗
DOE OSTI2022

Accurate and Rapid Forecasts for Geologic Carbon Storage via Learning-Based Inversion-Free Prediction

Carbon capture and storage (CCS) is one approach being studied by the U.S. Department of Energy to help mitigate global warming. The process involves capturing CO 2 emissions from industrial sources and permanently storing them in deep geologic formations (storage reservoirs). However, CCS projects generally target “green field sites,” where there is often little characterization data and therefore large uncertainty about the petrophysical properties and other geologic attributes of the storage reservoir. Consequently, ensemble-based approaches are often used to forecast multiple realizations prior to CO 2 injection to visualize a range of potential outcomes. In addition, monitoring data during injection operations are used to update the pre-injection forecasts and thereby improve agreement between forecasted and observed behavior. Thus, a system for generating accurate, timely forecasts of pressure buildup and CO 2 movement and distribution within the storage reservoir and for updating those forecasts via monitoring measurements becomes crucial. This study proposes a learning-based prediction method that can accurately and rapidly forecast spatial distribution of CO 2 concentration and pressure with uncertainty quantification without relying on traditional inverse modeling. The machine learning techniques include dimension reduction, multivariate data analysis, and Bayesian learning. The outcome is expected to provide CO 2 storage site operators with an effective tool for timely and informative decision making based on limited simulation and monitoring data.

58 GEOSCIENCES↗
DOE OSTI2022

3D Hydrodynamics of Pre-supernova Outbursts in Convective Red Supergiant Envelopes

Eruptive mass loss likely produces the energetic outbursts observed from some massive stars before they become core-collapse supernovae (SNe). The resulting dense circumstellar medium may also cause the subsequent SNe to be observed as Type IIn events. The leading hypothesis of the cause of these outbursts is the response of the envelope of the red supergiant (RSG) progenitor to energy deposition in the months to years prior to collapse. Early theoretical studies of this phenomenon were limited to 1D, leaving the 3D convective RSG structure unaddressed. Using FLASH's hydrodynamic capabilities, we explore the 3D outcomes by constructing convective RSG envelope models and depositing energies less than the envelope binding energies on timescales shorter than the envelope dynamical time deep within them. We confirm the 1D prediction of an outward-moving acoustic pulse steepening into a shock, unbinding the outermost parts of the envelope. However, we find that the initial 2–4 km s -1 convective motions seed the intrinsic convective instability associated with the high-entropy material deep in the envelope, enabling gas from deep within the envelope to escape and increasing the amount of ejected mass compared to an initially "quiescent" envelope. The 3D models reveal a rich density structure, with column densities varying by ≈10× along different lines of sight. Our work highlights that the 3D convective nature of RSG envelopes impacts our ability to reliably predict the outburst dynamics, the amount, and the spatial distribution of the ejected mass associated with deep energy deposition.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

DESI Observations of the Andromeda Galaxy: Revealing the Immigration History of Our Nearest Neighbor

Abstract We present Dark Energy Spectroscopic Instrument (DESI) observations of the inner halo of M31, which reveal the kinematics of a recent merger—a galactic immigration event—in exquisite detail. Of the 11,416 sources studied in 3.75 hr of on-sky exposure time, 7438 are M31 sources with well-measured radial velocities. The observations reveal intricate coherent kinematic structure in the positions and velocities of individual stars: streams, wedges, and chevrons. While hints of coherent structures have been previously detected in M31, this is the first time they have been seen with such detail and clarity in a galaxy beyond the Milky Way. We find clear kinematic evidence for shell structures in the Giant Stellar Stream, the Northeast Shelf, and Western Shelf regions. The kinematics are remarkably similar to the predictions of dynamical models constructed to explain the spatial morphology of the inner halo. The results are consistent with the interpretation that much of the substructure in the inner halo of M31 is produced by a single galactic immigration event 1–2 Gyr ago. Significant numbers of metal-rich stars ([Fe/H] > − 0.5) are present in all of the detected substructures, suggesting that the immigrating galaxy had an extended star formation history. We also investigate the ability of the shells and Giant Stellar Stream to constrain the gravitational potential of M31, and estimate the mass within a projected radius of 125 kpc to be log 10 M NFW ( < 125 kpc ) / M ⊙ = 11.80 − 0.10 + 0.12 . The results herald a new era in our ability to study stars on a galactic scale and the immigration histories of galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A staged deep learning approach to spatial refinement in 3D temporal atmospheric transport

High-resolution spatiotemporal simulations effectively capture the complexities of atmospheric plume dispersion in complex terrain. However, their high computational cost makes them impractical for applications requiring rapid responses or iterative processes, such as optimization, uncertainty quantification, or inverse modeling. To address this challenge, this work introduces the Dual-Stage Temporal Three-dimensional UNet Super-resolution (DST3D-UNet-SR) model, a highly efficient deep learning model for plume dispersion predictions. DST3D-UNet-SR is composed of two sequential modules: the temporal module (TM), which predicts the transient evolution of a plume in complex terrain from low-resolution temporal data, and the spatial refinement module (SRM), which subsequently enhances the spatial resolution of the TM predictions. We train DST3D-UNet-SR using a comprehensive dataset derived from high-resolution large eddy simulations (LES) of plume transport. We propose the DST3D-UNet-SR model to significantly accelerate LES of three-dimensional (3D) plume dispersion by three orders of magnitude. Additionally, the model demonstrates the ability to dynamically adapt to evolving conditions through the incorporation of new observational data, substantially improving prediction accuracy in high-concentration regions near the source.

3D temporal sequences↗
DOE OSTI2024

Development of spatially and temporally resolved electron density measurements for the assessment of radiation hydrodynamics simulations of laboratory X-ray photoionized plasmas

The photoionized plasma gas cell experiment is an established platform we use to make at-parameter (ξ > > 1 ergs cm s –1 ) measurements of plasma properties with application to high-energy astrophysical systems. We model the experiments with 1D radiation hydrodynamics simulations using the HELIOS-CR code to inform our understanding and assist in the interpretation of results. The simulations predict that the bulk of the plasma is in a quasi-uniform and hydrodynamically unperturbed state throughout the duration of the experiment. To evaluate this prediction, we introduced a photonic Doppler velocimetry (PDV) diagnostic to measure spatially and temporally resolved plasma electron density. The initial measurements were successful but had limitations that made model-data comparisons challenging. To address this, we re-designed the gas cell PDV diagnostic and doubled the number of measurement locations to sample across two thirds of the depth of the cell. We also present a comparison of the results from the upgraded PDV diagnostic to the HELIOS-CR simulations for the first time. As a result, the experimental data confirms the prediction of an unperturbed region in the bulk of the plasma but reveals discrepancies in the time evolution and spatial distribution of the simulated electron density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Physics-Informed Deep Learning for Reconstruction of Spatial Missing Climate Information in the Antarctic

Understanding the influence of the Antarctic on the global climate is crucial for the prediction of global warming. However, due to very few observation sites, it is difficult to reconstruct the rational spatial pattern by filling in the missing values from the limited site observations. To tackle this challenge, regional spatial gap-filling methods, such as Kriging and inverse distance weighted (IDW), are regularly used in geoscience. Nevertheless, the reconstructing credibility of these methods is undesirable when the spatial structure has massive missing pieces. Inspired by image inpainting, we propose a novel deep learning method that demonstrates a good effect by embedding the physics-aware initialization of deep learning methods for rapid learning and capturing the spatial dependence for the high-fidelity imputation of missing areas. We create the benchmark dataset that artificially masks the Antarctic region with ratios of 30%, 50% and 70%. The reconstructing monthly mean surface temperature using the deep learning image inpainting method RFR (Recurrent Feature Reasoning) exhibits an average of 63% and 71% improvement of accuracy over Kriging and IDW under different missing rates. With regard to wind speed, there are still 36% and 50% improvements. In particular, the achieved improvement is even better for the larger missing ratio, such as under the 70% missing rate, where the accuracy of RFR is 68% and 74% higher than Kriging and IDW for temperature and also 38% and 46% higher for wind speed. In addition, the PI-RFR (Physics-Informed Recurrent Feature Reasoning) method we proposed is initialized using the spatial pattern data simulated by the numerical climate model instead of the unified average. Compared with RFR, PI-RFR has an average accuracy improvement of 10% for temperature and 9% for wind speed. When applied to reconstruct the spatial pattern based on the Antarctic site observations, where the missing rate is over 90%, the proposed method exhibits more spatial characteristics than Kriging and IDW.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Machine learning enables national assessment of wind plant controls with implications for land use

Summary As deployment of wind energy continues to expand, computationally efficient tools for predicting wind plant performance over a wide range of layout designs, technology innovations, and spatial locations are increasingly important for policy and investment decisions. We demonstrate two approaches to training a surrogate model to predict annual energy production (AEP) of parameterized wind plant layouts: one using a Gaussian process (GP) and the other using a fully convolutional neural network (FCNN). We leverage the powerful FCNN architecture by encoding wind plant design parameters and output response surface as an image. The FCNN produces more accurate results than the GP with mean absolute errors equivalent to 1% and 1.9% of plant rated power, respectively, although the GP performs well under limited training data and provides useful uncertainty information. We also evaluate a surrogate model for wake steering, enabling a nationwide assessment of the impact of plant control strategies and plant layout decisions. Across two million locations, we find that wake steering strategies boost AEP with relative gains upwards of 3%. Gains are most pronounced at sites without a dominant wind direction and where layout optimization is less fruitful. Additionally, we perform a nationwide sensitivity analysis showing that wake steering can mitigate wake losses from higher density plant layouts. Our results suggest that regions which have not been previously viable for wind deployment due to moderate wind resources are especially well enhanced by wake steering strategies that could help overcome land constraints and inflexible layout options, potentially identifying new deployment opportunities.

17 WIND ENERGY↗
DOE OSTI2025

Comparison of spatial dynamics and point kinetics approaches in multiphysics modeling of the molten salt reactor experiment

In this work, we present validation test results of fully coupled neutronics and thermal-hydraulics models of the Molten Salt Reactor Experiment (MSRE) against experimental data of the zero power pump transients and the natural circulation tests at low power. To capture the strong coupling between neutronics and thermal-hydraulics due to fuel circulation, and to account for the delayed neutron precursor (DNP) distribution, the porous media thermal-hydraulics solver Pronghorn was fully coupled to the spatial neutron dynamics code Griffin, which solves the neutron diffusion equation, and to the 0-D point kinetics solver Squirrel, using a 2-D homogenized representation of the MSRE. The validation test results show very good agreement with experimental data for both point kinetics and spatial dynamics simulations, capturing the strong feedback effect and DNP losses in the MSRE. The 0-D code Squirrel accurately predicted the time-dependent behavior in the MSRE given the steady-state spatial dynamics solution of Griffin.

42 - ENGINEERING↗
DOE OSTI2022

Deep Learning Provides Substantial Improvements to County‐Level Fire Weather Forecasting Over the Western United States

Abstract The recent wildfires in the western United States during 2018 and 2020 caused record‐breaking fire damage and casualties. Despite remarkable advances in fire modeling and weather forecasting, it remains challenging to anticipate catastrophic wildfire events and associated damage. One key missing component is a fire weather prediction system with sufficiently long lead time capable of providing useful regional details. Here, we develop a hybrid prediction model of wildfire danger called CFS with super resolution (CFS‐SR) as a proof of concept to fill that void. The CFS‐SR model is constructed by integrating the Climate Forecast System version 2 with a deep learning (DL) technique from Single Image Super Resolution, a method widely used in enhancing image resolution. We show that for the 2018–2019 fire season, the CFS‐SR model significantly improves accuracy in forecasting fire weather at lead times of up to 7 days with an enhanced spatial resolution up to 4 km. This level of high resolution provides county‐level fire weather forecast, making it more practical for allocating resources to mitigate wildfire danger. Our study demonstrates that a proper combination of ensemble climate predictions with DL techniques can boost predictability at finer spatial scales, increasing the utility of fire weather forecasts for practical applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Understanding Hailstone Temporal Variability and Contributing Factors over the U.S. Southern Great Plains

Abstract Hailstones are a natural hazard that pose a significant threat to property and are responsible for significant economic losses each year in the United States. Detailed understanding of their characteristics is essential to mitigate their impact. Identifying the dynamic and physical factors contributing to hail formation and hailstone sizes is of great importance to weather and climate prediction and policymakers. In this study, we have analyzed the temporal and spatial variabilities of severe hail occurrences over the U.S. southern Great Plains (SGP) states from 2004 to 2016 using two hail datasets: hail reports from the Storm Prediction Center and the newly developed radar-retrieved maximum expected size of hail (MESH). It is found that severe and significant severe hail occurrences have a considerable year-to-year temporal variability in the SGP region. The interannual variabilities have a strong correspondence with sea surface temperature anomalies over the northern Gulf of Mexico and there is no outlier. The year 2016 is identified as an outlier for the correlations with both El Niño–Southern Oscillation (ENSO) and aerosol loading. The correlations with ENSO and aerosol loading are not statistically robust to inclusion of the outlier 2016. Statistical analysis without the outlier 2016 shows that 1) aerosols that may be mainly from northern Mexico have the largest correlation with hail interannual variability among the three factors and 2) meteorological covariation does not significantly contribute to the high correlation. These analyses warrant further investigations of aerosol impacts on hail occurrence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Predicting Fracture Network Development in Crystalline Rocks

The geometric properties of fractures influence whether they propagate, arrest, or coalesce with other fractures. Thus, quantifying the relationship between fracture network characteristics may help predict fracture network development, and perhaps precursors to catastrophic failure. To constrain the relationship and predictability of fracture characteristics, we deform eight one centimeter tall rock cores under triaxial compression while acquiring in situ X-ray tomograms. The tomograms reveal precise measurements of the fracture network characteristics above the spatial resolution of 6.5 µm. We develop machine learning models to predict the value of each characteristic using the other characteristics, and excluding the macroscopic stress or strain imposed on the rock. The models predict fracture development more accurately in the experiments performed on granite and monzonite, than the experiments on marble. Fracture network development may be more predictable in these igneous rocks because their microstructure is more mechanically homogeneous than the marble, producing more systematic fracture development that is not strongly impeded by grain contacts and cleavage planes. The varying performance of the models suggest that fracture volume, length, and aperture are the most predictable of the characteristics, while fracture orientation is the least predictable. Orientation does not correlate with length, as suggested by the idea that the orientation evolves with increasing differential stress and thus fracture length. This difference between the observed and expected relationship between orientation and length highlights the influence of mechanical heterogeneities and local stress perturbations on fracture growth as fractures propagate, link, and coalesce.

58 GEOSCIENCES↗
DOE OSTI2025

Convergent Manufacturing of Large-Scale Components for Nuclear Applications, via Additive Manufacturing and Powder Metallurgy Hot Isostatic Pressing

Powder metallurgy (PM)–hot isostatic pressing (PM-HIP) has long been recognized as a powerful route for producing fully dense, near net shape metallic components. By consolidating powders under high temperature and pressure, HIP provides isotropic properties, uniform microstructures, and scalability to complex geometries that are vital for sectors such as aerospace, energy, and nuclear power. Yet despite these advantages, the technology has remained constrained by costly trial and error canister fabrication, limitations of conventional forging, and incomplete knowledge about how the canister design influences final part properties. Additive manufacturing (AM), by contrast, thrives on design freedom and geometric flexibility but struggles with speed, scalability, and cost when applied to very large structures. The research presented in this report investigated how a convergent manufacturing approach, combining AM with PM-HIP, can merge the strengths of both technologies, leveraging AM’s flexibility for canister design and HIP’s consolidation capability to deliver reliable, large, and complex parts. The work progressed through three case studies that built on one another in scale and complexity. Small cylindrical canisters fabricated by conventional methods, laser powder bed fusion, and directed energy deposition were filled with stainless steel powders and subjected to HIP. The resulting parts demonstrated near-full density and mechanical properties on par with wrought stainless steel, showing for the first time that AM canisters can be a direct substitute for conventional ones without sacrificing quality. The next step involved a medium-scale, noncentrosymmetric T-valve, which is an enclosed, multibranch geometry that tested the limits of AM + PM-HIP integration. The T-valve achieved predictable shrinkage and uniform densification, confirming feasibility for enclosed designs. However, this study also revealed oxide inclusions and interfacial challenges at the AM + HIP boundary, underscoring the critical importance of controlling interface chemistry and employing robust, in situ strategies, such as melt pool monitoring and thermal monitoring, coupled with nondestructive evaluation techniques such as x-ray computed tomography. Finally, the effort culminated in fabricating a large-scale impeller weighing nearly 2000 lb and spanning 5 ft in diameter. Produced via multirobot wire arc AM and hot isostatic pressed to near-full density, the impeller validated industrial-scale feasibility. Predictive models closely matched experimental shrinkage, tensile properties were spatially uniform across the component, and the AM + PM-HIP interface proved mechanically sound despite the presence of oxide-decorated prior particle boundaries. This large-scale demonstration is a major milestone, showing that hybrid AM + PM‑HIP can reliably deliver components at reactor-relevant scales. Collectively, these studies charted a logical pathway: small-scale work built scientific confidence, medium-scale work highlighted opportunities and challenges, and large-scale work proved industrial impact. The overarching conclusion of this report is that AM + PM-HIP should not be seen as a replacement for forging but as a complementary pathway that provides the US with flexibility, resilience, and new options for manufacturing nuclear-grade components. Looking ahead, several directions emerge as critical to sustaining progress. Predictive modeling must become faster, more accessible, and more accurate, with digital twins and machine learning reducing reliance on trial and error. Powders and alloys must be optimized for HIP, with improved cleanliness, reduced oxides, and tailored chemistries that enhance creep, fatigue, and irradiation resistance. Interfaces between AM and HIP regions must be better engineered through coatings, machining strategies, and surface treatments to mitigate oxide formation and ensure reliable bonding to explore opportunities for HIP of targeted compositional parts, as well as multimaterial HIP cladding applications. Monitoring and nondestructive evaluation need to expand, incorporating multimodal sensors, x-ray computed tomography, and real-time data integration through platforms such as Pelican. At the same time, the pathway to industrial adoption requires techno-economic analysis, machinability studies, and qualification frameworks aligned with industry and regulatory standards. Finally, workforce and academic engagement must be strengthened. Programs that train technicians and engineers for US Navy and US Department of Energy manufacturing challenges should be paired with academic partnerships to support fundamental research, with open sharing of non-export-controlled data to accelerate innovation and build the next generation of experts. In conclusion, this report demonstrates that hybrid AM + PM-HIP is scientifically viable and strategically important. By combining the design agility of AM with the consolidation strength of HIP and embedding modeling, monitoring, and workforce development, this approach provided a transformative new capability for US manufacturing. The path forward is clear: hybrid AM + PM-HIP is not just a promising research direction but is also potentially an industrially relevant pathway that can reshape how nuclear-grade components are designed, qualified, and deployed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Simulated plant-mediated oxygen input has strong impacts on fine-scale porewater biogeochemistry and weak impacts on integrated methane fluxes in coastal wetlands

Methane (CH 4 ) emissions from wetland ecosystems are controlled by redox conditions in the soil, which are currently underrepresented in Earth system models. Plant-mediated radial oxygen loss (ROL) can increase soil O 2 availability, affect local redox conditions, and cause heterogeneous distribution of redox-sensitive chemical species at the root scale, which would affect CH 4 emissions integrated over larger scales. In this study, we used a subsurface geochemical simulator (PFLOTRAN) to quantify the effects of incorporating either spatially homogeneous ROL or more complex heterogeneous ROL on model predictions of porewater solute concentration depth profiles (dissolved organic carbon, methane, sulfate, sulfide) and column integrated CH 4 fluxes for a tidal coastal wetland. From the heterogeneous ROL simulation, we obtained 18% higher column averaged CH 4 concentration at the rooting zone but 5% lower total CH 4 flux compared to simulations of the homogeneous ROL or without ROL. This difference is because lower CH 4 concentrations occurred in the same rhizosphere volume that was directly connected with plant-mediated transport of CH 4 from the rooting zone to the atmosphere. Sensitivity analysis indicated that the impacts of heterogeneous ROL on model predictions of porewater oxygen and sulfide concentrations will be more important under conditions of higher ROL fluxes or more heterogeneous root distribution (lower root densities). Despite the small impact on predicted CH 4 emissions, the simulated ROL drastically reduced porewater concentrations of sulfide, an effective phytotoxin, indicating that incorporating ROL combined with sulfur cycling into ecosystem models could potentially improve predictions of plant productivity in coastal wetland ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Predicting mesoscale spectral thermal conductivity using advanced deterministic phonon transport techniques

We present a review, demonstration and simulation of phonon transport for the purposes of predicting materials performance at the mesoscale. We focus primarily on the development and implementation of a unified methodology to enable predictive heat transport. We report on the current state of the art as it pertains to deterministic phonon transport methodologies, discussing various topics concerning phonons. In application, we focus on the Self-Adjoint Angular Flux (SAAF) formulation of the Boltzmann transport equation for phonons, and develop the spatial, angular, and material property discretization required to accurately simulate the predictive physics of heat transport in dielectrics. We discuss thermal interfacial resistance and present our formulation of the diffuse mismatch model for simulating phonon interactions at internal boundaries. We have recently developed a deterministic, spectral phonon transport method for predicting effective thermal conductivity ($\kappa_{\textrm{eff}}$), using Bose-Einstein source terms coupled through an average material temperature. In this method, we introduce a closure term to the phonon transport system which acts as a redistribution function for the total energy of the system, and serves as a guide for the amount of non-equilibrium behavior occurring in the system. This method predicts thermal conductivity and equilibrium temperature distributions in homogeneous and heterogeneous materials using data generated by ab initio density functional theory methods. We employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. The equations associated with this method are solved via a modification of traditional source iteration. We compare the performance of source iteration applied to an existing uncoupled, traditional SAAF method to our new method and comment on the iterative performance of each. We observe ballistic and diffusive phonon scattering as acoustic thickness of the domain changes, and are able to make comparisons between the accuracy and efficiency of both methods.

36 MATERIALS SCIENCE↗