SEARCH · Engineering Papers

Results for “Spatial dependence”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

DOE OSTI2021

Depth-resolved cathodoluminescence and surface photovoltage spectroscopies of gallium vacancies in β-Ga 2 O 3 with neutron irradiation and forming gas anneals

The gallium vacancy is one of the dominant native point defects in β-Ga 2 O 3 , one that, together with its complexes, can have a major effect on free carrier densities and transport in this wide bandgap semiconductor. We used a combination of depth-resolved cathodoluminescence spectroscopy and surface photovoltage spectroscopy to identify the optical and energy-level properties of these defects as well as how their defect densities and spatial distributions vary with neutron irradiation and temperature-dependent-forming gas anneals. These studies reveal optical signatures that align closely with theoretical energy-level predictions. Likewise, our optical techniques reveal variations in these defect densities that are consistent with hydrogen passivation of gallium vacancies as a function of temperature and depth from the free Ga 2 O 3 surface. These techniques can help guide the understanding and control of dominant native point defects in Ga 2 O 3 .

42 ENGINEERING↗

DOE OSTI2022

Fast Solution of Fully Implicit Runge--Kutta and Discontinuous Galerkin in Time for Numerical PDEs, Part II: Nonlinearities and DAEs

Fully implicit Runge--Kutta (IRK) methods have many desirable accuracy and stability properties as time integration schemes, but high-order IRK methods are not commonly used in practice with large-scale numerical PDEs because of the difficulty of solving the stage equations. This paper introduces a theoretical and algorithmic framework for solving the nonlinear equations that arise from IRK methods (and discontinuous Galerkin discretizations in time) applied to nonlinear numerical PDEs, including PDEs with algebraic constraints. Several new linearizations of the nonlinear IRK equations are developed, offering faster and more robust convergence than the often-considered simplified Newton, as well as an effective preconditioner for the true Jacobian if exact Newton iterations are desired. Inverting these linearizations requires solving a set of block 2 x 2 systems. Under quite general assumptions, it is proven that the preconditioned 2 x 2 operator's condition number is bounded by a small constant close to one, independent of the spatial discretization, spatial mesh, and time step, and with only weak dependence on the number of stages or integration accuracy. Moreover, the new method is built using the same preconditioners needed for backward Euler-type time stepping schemes, so can be readily added to existing codes. The new methods are applied to several challenging fluid flow problems, including the compressible Euler and Navier--Stokes equations, and the vorticity-streamfunction formulation of the incompressible Euler and Navier--Stokes equations. Up to 10th-order accuracy is demonstrated using Gauss IRK, while in all cases fourth-order Gauss IRK requires roughly half the number of preconditioner applications as required by standard Singly diagonally implicit Runge--Kutta methods.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

LLNL FY22 Aging and Lifetimes Exit Criteria (Milestone 8118, GC#5)

This report introduces a novel characterization tool that can rapidly collect in-situ time-dependent oxide thickness and hydride growth data together with spatial data using while-light interferometry. White-light interferometry is a non-destructive and non-contact optical surface profiling technique for materials characterization that offers excellent lateral and vertical resolution. A depleted uranium sample was used to assess the ability of white-light interferometry to provide quantitative corrosion reaction rates at metal surfaces for by hydriding and oxidation and compare this technique to existing characterization tools for capturing early-stage hydride and oxide growth kinetics.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Modeling Carbon-Based Nanomaterials (CNMs) and Derived Composites and Devices

A review of different modeling techniques, specifically in the framework of carbon-based nanomaterials (CNMs, including nanoparticles such as graphene and carbon nanotubes—CNTs) and the composites and devices that can be derived from them, is presented. The article emphasizes that the overall performance of these materials depends on mechanisms that operate across different time and spatial scales, requiring tailored approaches based on the material type, size, internal structure/configuration, and the specific properties of interest. Far from attempting to cover the entire spectrum of models, this review examines a wide range of analysis and simulation techniques, highlighting their potential use, some of their weaknesses and strengths, and presenting the latest developments and some application examples. In this way, it is shown how modeling can provide key information for tailoring or designing new materials for specific components or systems or to obtain certain functionalities. At the same time, it is revealed to be an area constantly undergoing development and improvement, as evidenced by the progress made by various of these techniques and the new modeling approaches that have emerged in recent years.

Chiminelli, Agustίn (ORCID:0000000230312107)↗

DOE OSTI2021

Microstructure-dependent rate theory model of defect segregation and phase stability in irradiated polycrystalline LiAlO 2

We report gamma lithium aluminate (LiAlO 2 ) is a breeder material for tritium and is one of key components in a tritium-producing burnable absorber rod (TPBAR). Dissolution and precipitation of second phases such as LiAl 5 O 8 and voids are observed in irradiated LiAlO 2 . Such microstructure changes cause the degradation of thermomechanical properties of LiAlO 2 and affect tritium retention and release kinetics, and hence, the TPBAR performance. In this work, a microstructure-dependent model of radiation-induced segregation (RIS) has been developed for investigating the accumulation of species and phase stability in polycrystalline LiAlO 2 structures under irradiation. Three sublattices (i.e. [Li, Al, V] I [O, V o ] II [Li i , Al i , O i , V i ] III ), and concentrations of six diffusive species (i.e. Li; vacancy of Li or Al at [Li, Al, V] I sublattice, O vacancy at [O, V o ] II sublattice, and Li, Al and O interstitials at [Li i , Al i , O i , V i ] III interstitial sublattices; are used to describe spatial and temporal distributions of defects and chemistry. Microstructure-dependent thermodynamic and kinetic properties including the generation, reaction, and chemical potentials of defects and defect mobility are taken into account in the model. The parametric studies demonstrated the capability of the developed RIS model to assess the effect of thermodynamic and kinetic properties of defects on the segregation and depletion of species in polycrystalline structures and to explain the phase stability observed in irradiated LiAlO 2 samples. The developed RIS model will be extended to study the precipitation of LiAl 5 O 8 and voids and tritium retention by integrating the phase-field method.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Spatially Patterned Architectures to Modulate CO 2 Reduction Cascade Catalysis Kinetics

Electrochemical CO 2 reduction using renewable sources of electrical energy holds promise for converting CO 2 into fuels and chemicals. The complex interactions among chemical/electrochemical reactions and mass transport make it difficult to analyze the effect of an individual process on electrode performance based only on experimental methods. Here, we developed a generalized steady-state simulation to describe an electrode surface in which sequential cascade catalysts are patterned in a periodic trench design. If appropriately constructed, this trench geometry is hypothesized to be able to yield a higher net current density for a CO 2 reduction (CO 2 R) cascade reaction. We have used realistic experimental reaction kinetics to investigate the role of trench geometry in mass transport, local microenvironments, and selectivity for a model CO 2 R cascade reaction. The model considers local concentration gradients of bicarbonate species at quasi-equilibrium and catalytic surface reactions based on concentration-dependent Butler–Volmer kinetics. Our results suggest that varying the spatial distribution of active sites plays a significant role in facilitating effective mass transport between active sites, modulating selectivity for the cascade reaction, and enhancing the yield of desirable cascade products. Moreover, we observe that this trench geometry significantly alters the cascade reaction rate by affecting the local pH, which can cause inadvertent depletion of available aqueous CO 2 to limit the CO 2 R cascade kinetics and modest suppression of the hydrogen evolution reaction (HER). The results highlight the trade-offs between mass transport, pH, and reaction kinetics that become apparent only when considering the coupled physics of all processes at the electrode surface. Here, this model can thus serve as a primary tool to build more selective and efficient patterned architectures for the CO 2 R cascade catalysis.

CO2 reduction↗

DOE OSTI2020

Massively parallel transport sweeps on meshes with cyclic dependencies

When solving the first-order form of the linear Boltzmann equation, a common misconception is that the matrix-free computational method of “sweeping the mesh”, used in conjunction with the Discrete Ordinates method, is too complex or does not scale well enough to be implemented in modern high performance computing codes. This has led to considerable efforts in the development of matrix-based methods that are computationally expensive and is partly driven by the requirements placed on modern spatial discretizations. In particular, modern transport codes are required to support higher order elements, a concept that invariably adds a lot of complexity to sweeps because of the introduction of cyclic dependencies with curved mesh cells. In this article we will present a comprehensive implementation of sweeping, to a piecewise-linear DFEM spatial discretization with particular focus on handling cyclic dependencies and possible extensions to higher order spatial discretizations. We find that these methods are implemented in a new C++ simulation framework called Chi-Tech (). We present some typical simulation results with some performance aspects that one can expect during real world simulations, we also present a scaling study to >100k processes where Chi-Tech maintains greater than 80% efficiency solving a total of 87.7 trillion angular flux unknowns for a 116 group simulation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2024

Quantifying Dissolution Dynamics in Porous Media Using a Spatial Flow Focusing Profile

Abstract The diverse range of patterns in porous media formed by dissolution processes depends on the relative magnitude of flow, transport, and chemical reactions at pore surfaces. However, distinguishing between regimes often relies solely on qualitative, visual comparisons of emergent structures. Here, we propose a quantitative measure capable of identifying different regimes using the concept of the spatial flow focusing profile, which segments the medium into cross sections along the flow direction to calculate the flow focusing index for each section. We employ this measure in numerical simulations of a dissolving porous medium using a pore network model. We obtain a morphological phase diagram of dissolution patterns, which we characterize using the flow focusing profile. In particular, we demonstrate that analyzing the temporal changes in the profile allows one to quantitatively distinguish between wormholing and channeling. The transition between them is shown to be affected by the heterogeneity of the system.

58 GEOSCIENCES↗

DOE OSTI2020

Steric effects in light-induced solvent proton abstraction

The significance of solvent structural factors in the excited-state proton transfer (ESPT) reactions of Schiff bases with alcohols is reported here. In this work, we use the super photobase FR0-SB and a series of primary, secondary, and tertiary alcohol solvents to illustrate the steric issues associated with solvent to photobase proton transfer. Steady-state and time-resolved fluorescence data show that ESPT occurs readily for primary alcohols, with a probability proportional to the relative –OH concentration. For secondary alcohols, ESPT is greatly diminished, consistent with the barrier heights obtained using quantum chemistry calculations. ESPT is not observed in the tertiary alcohol. We explain ESPT using a model involving an intermediate hydrogen-bonded complex where the proton is “shared” by the Schiff base and the alcohol. The formation of this complex depends on the ability of the alcohol solvent to achieve spatial proximity to and alignment with the FR0-SB* imine lone pair stabilized by the solvent environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Temporal and Spatial Evolution of Non-Elastic Strain Accumulation in Stanstead Granite During Brittle Creep

Understanding the long-term behavior of brittle rocks requires fundamental consideration of time-dependent strain evolution and brittle creep processes. Previous studies have evaluated sub-critical crack growth during time-dependent deformation and damage evolution in brittle rocks; however, there is an incomplete knowledge of how damage evolves spatially and temporally within the body of intact rocks, where distributed regions of damage interact and coalesce during creep. This paper presents laboratory research focusing on evaluating brittle creep damage processes in Stanstead granite (SG) using 2-dimensional digital image correlation (2D-DIC). In the laboratory, the prismatic SG specimens were loaded beyond an estimated Crack Damage stress threshold (CD) level and then maintained a constant stress to initiate the creep process. DIC was used to characterize full-field spatiotemporal strain evolution, which was then interpreted in the context of local regions of “damage”, determined according to a strain-based criterion. Here, a method was proposed for identifying “existing” and “new” damage regions over specified intervals during the test, followed by spatial clustering of these regions to assess their spatiotemporal evolution. The clustering analysis results demonstrated the extension of existing damage regions was the main damage process during brittle creep, which is consistent with existing models of sub-critical crack growth. In addition, temporal analysis of tensile and shear strains on a point-by-point basis revealed both new damage formation and the strain concentration within existing damaged regions significantly contribute to overall specimen strain during primary creep. In contrast, during secondary creep, increases in specimen deformation are influenced by the accumulation of strains within already damaged regions.

58 GEOSCIENCES↗

DOE OSTI2025

Dependent scattering and fractal microstructure determine the transparency of aerogel monoliths

This study reveals how dependent scattering and microstructure significantly affect electromagnetic wave propagation through aerogel monoliths, contributing to their transparency. Light scattering by particle ensembles is considered “dependent” when the scattering properties rely not only on particle size and optical constants but also on their spatial distribution, typically occurring when the average interparticle distance is small in comparison with the wavelength of incident radiation. Addressing dependent scattering requires solving Maxwell’s equations for complex heterogeneous structures, which is computationally demanding and usually limited to sample thicknesses on the same scale as the wavelength. This study combines computer-generated ambigel microstructures of fractal aggregates of polydisperse nanoparticles and the radiative transfer with reciprocal transaction method to predict the transmittance of thick ambigel slabs. Transmittance measurements of ambiently dried aerogel monoliths (ambigels) with porosities from about 50% to 90% closely matched the predicted values for their digital twins. However, ignoring dependent scattering or particle aggregation led to inaccurate predictions. This study validated the computational framework, and its findings offer insights for designing photonic metamaterials and analyzing their interactions with electromagnetic waves.

Yalcin, Refet A. (ORCID:0000000339973494)↗

DOE OSTI2023

Limited Role of Absolute Humidity in Intraurban Heat Variability

Abstract Monitoring and understanding the variability of heat within cities is important for urban planning and public health, and the number of studies measuring intraurban temperature variability is growing. Recognizing that the physiological effects of heat depend on humidity as well as temperature, measurement campaigns have included measurements of relative humidity alongside temperature. However, the role the spatial structure in humidity, independent from temperature, plays in intraurban heat variability is unknown. Here we use summer temperature and humidity from networks of stationary sensors in multiple cities in the United States to show spatial variations in the absolute humidity within these cities are weak. This variability in absolute humidity plays an insignificant role in the spatial variability of the heat index and humidity index (humidex), and the spatial variability of the heat metrics is dominated by temperature variability. Thus, results from previous studies that considered only intraurban variability in temperature will carry over to intraurban heat variability. Also, this suggests increases in humidity from green infrastructure interventions designed to reduce temperature will be minimal. In addition, a network of sensors that only measures temperature is sufficient to quantify the spatial variability of heat across these cities when combined with humidity measured at a single location, allowing for lower-cost heat monitoring networks. Significance Statement Monitoring the variability of heat within cities is important for urban planning and public health. While the physiological effects of heat depend on temperature and humidity, it is shown that there are only weak spatial variations in the absolute humidity within nine U.S. cities, and the spatial variability of heat metrics is dominated by temperature variability. This suggests increases in humidity will be minimal resulting from green infrastructure interventions designed to reduce temperature. It also means a network of sensors that only measure temperature is sufficient to quantify the spatial variability of heat across these cities when combined with humidity measured at a single location.

Meteorology & Atmospheric Sciences↗

DOE OSTI2024

Flat Spectra of Energetic Particles in Interplanetary Shock Precursors

The observed energy spectra of accelerated particles at interplanetary shocks often do not match the diffusive shock acceleration (DSA) theory predictions. In some cases, the particle flux forms a plateau over a wide range of energies, extending upstream of the shock for up to seven flux e-folds before submerging into the background spectrum. Remarkably, at and downstream of the shock we have studied in detail, the flux falls off in energy as ϵ -1 , consistent with the DSA prediction for a strong shock. The upstream plateau suggests a particle transport mechanism different from those traditionally employed in DSA models. We show that a standard (linear) DSA solution based on a widely accepted diffusive particle transport with an underlying resonant wave–particle interaction is inconsistent with the plateau in the particle flux. To resolve this contradiction, we modify the DSA theory in two ways. First, we include a dependence of the particle diffusivity κ on the particle flux F (nonlinear particle transport). Second, we invoke short-scale magnetic perturbations that are self-consistently generated by, but not resonant with, accelerated particles. They lead to the particle diffusivity increasing with the particle energy as ∝ϵ 3/2 that simultaneously decreases with the particle flux as 1/F. The combination of these two trends results in the flat spectrum upstream. We speculate that nonmonotonic spatial variations of the upstream spectrum, apart from being time-dependent, may also result from non-DSA acceleration mechanisms at work upstream, such as stochastic Fermi or magnetic pumping acceleration.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2025

Inhomogeneous energy injection in the 21-cm power spectrum: Sensitivity to dark matter decay

The 21-cm signal provides a novel avenue to measure the thermal state of the Universe during cosmic dawn and reionization (redshifts z ~ 5–30), and thus to probe energy injection from decaying or annihilating dark matter (DM). These DM processes are inherently inhomogeneous: both decay and annihilation are density-dependent, and furthermore, the fraction of injected energy that is deposited at each point depends on the gas ionization and density, leading to further anisotropies in absorption and propagation. In this work, we develop a new framework for modeling the impact of spatially inhomogeneous energy injection and deposition during cosmic dawn, accounting for ionization and baryon density dependence, as well as the attenuation of propagating photons. We showcase how this first completely inhomogeneous treatment affects the predicted 21-cm power spectrum in the presence of exotic sources of energy injection, and forecast the constraints that upcoming HERA measurements of the 21-cm power spectrum will set on DM decays to photons and to electron/positron pairs. These projected constraints considerably surpass those derived from CMB and Lyman-a measurements, and for decays to electron/positron pairs they exceed all existing constraints in the sub-GeV mass range, reaching lifetimes of ~10 28 s. Our analysis demonstrates the unprecedented sensitivity of 21-cm cosmology to exotic sources of energy injection during the cosmic dark ages. Our code, dm21cm, includes all these effects and is publicly available in an accompanying release.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2021

Competition within low-density bacterial populations as an unexpected factor regulating carbon decomposition in bulk soil

Bacterial decomposition of organic matter in soils is generally believed to be mainly controlled by the access bacteria have to organic substrate. The influence of bacterial traits on this control has, however, received little attention. Using the substrate-dependent Monod growth model, we develop a bioreactive transport model to screen the interactive impacts of spatial dispersion and bacterial traits on mineralization. Bacterial traits primarily involved in the bacterial response to the substrate concentration, such as the maximum specific uptake rate and efficiency, the adaptation time of the uptake rate and the initial population density, are considered. We compare the model results with two sets of previously performed cm-scale soil-core experiments in which the mineralization of the pesticide 2,4-D was measured under well-controlled initial distributions and transport conditions. Bacterial dispersion away from the initial substrate location induced a significant increase in 2,4-D mineralization. Additionally, it reveals an increase of specific uptake rates at lower bacterial densities, more than compensating the decrease of specific uptake rates caused by substrate dilution. This regulation of bacterial activities by density, caused by the local depletion of substrate by competing bacteria, becomes dominant for bacteria with an efficient uptake of substrate at low substrate concentrations (a common feature of oligotrophs). Such oligotrophs, commonly found in soils, compete with each other for substrate even at remarkably low population densities. The ratio-dependent Contois growth model, which includes a density regulation in the expression of the uptake efficiency, is more accurate and convenient to calibrate than the substrate-dependent Monod model, at least under these conditions. In view of their strong interactions, bioreactive and transport processes cannot be handled independently but should be integrated, in particular when reactive processes of interest are carried out by oligotrophs.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

A solution framework for linear PDE-constrained mixed-integer problems

Abstract We present a general numerical solution method for control problems with state variables defined by a linear PDE over a finite set of binary or continuous control variables. We show empirically that a naive approach that applies a numerical discretization scheme to the PDEs to derive constraints for a mixed-integer linear program (MILP) leads to systems that are too large to be solved with state-of-the-art solvers for MILPs, especially if we desire an accurate approximation of the state variables. Our framework comprises two techniques to mitigate the rise of computation times with increasing discretization level: First, the linear system is solved for a basis of the control space in a preprocessing step. Second, certain constraints are just imposed on demand via the IBM ILOG CPLEX feature of a lazy constraint callback. These techniques are compared with an approach where the relations obtained by the discretization of the continuous constraints are directly included in the MILP. We demonstrate our approach on two examples: modeling of the spread of wildfire and the mitigation of water contamination. In both examples the computational results demonstrate that the solution time is significantly reduced by our methods. In particular, the dependence of the computation time on the size of the spatial discretization of the PDE is significantly reduced.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

Imputation of urban environmental sensor data using gated attention bidirectional long short-term memory (GA-BiLSTM): methods, performance, and implications

Urban environmental monitoring networks frequently encounter significant data gaps due to sensor malfunctions, environmental disturbances, and communication failures. Reliable approaches to address these gaps are essential for ensuring the continuity and quality of environmental data streams. In this study, we developed a gated attention bidirectional long short-term memory (GA-BiLSTM) model to impute missing data in a dense urban monitoring network. Using observations from the CROCUS network in Chicago, we evaluated GA-BiLSTM against widely used approaches (XGBoost and K-nearest neighbors) under scenarios of both short-term intermittent gaps and prolonged outages. GA-BiLSTM consistently outperformed comparative methods, particularly during extended outages of up to ten days, demonstrating its ability to capture spatiotemporal dependencies across sensor nodes. Beyond performance metrics, feature importance and spatial network analyses highlighted the unexpected but critical predictive role of peripheral rural nodes, underlining their strategic value for maintaining robust urban monitoring systems. These results emphasize that advanced imputation methods can substantially improve the reliability of environmental monitoring networks and support more resilient data infrastructures for urban sustainability.

Data imputation↗

DOE OSTI2022

Calibrating hypersonic turbulence flow models with the HIFiRE-1 experiment using data-driven machine-learned models.

In this paper we study the efficacy of combining machine-learning methods with projection-based model reduction techniques for creating data-driven surrogate models of computationally expensive, high-fidelity physics models. Such surrogate models are essential for many-query applications e.g., engineering design optimization and parameter estimation, where it is necessary to invoke the high-fidelity model sequentially, many times. Surrogate models are usually constructed for individual scalar quantities. However there are scenarios where a spatially varying field needs to be modeled as a function of the model’s input parameters. Here we develop a method to do so, using projections to represent spatial variability while a machine-learned model captures the dependence of the model’s response on the inputs. The method is demonstrated on modeling the heat flux and pressure on the surface of the HIFiRE-1 geometry in a Mach 7.16 turbulent flow. The surrogate model is then used to perform Bayesian estimation of freestream conditions and parameters of the SST (Shear Stress Transport) turbulence model embedded in the high-fidelity (Reynolds-Averaged Navier–Stokes) flow simulator, using shock-tunnel data. The paper provides the first-ever Bayesian calibration of a turbulence model for complex hypersonic turbulent flows. We find that the primary issues in estimating the SST model parameters are the limited information content of the heat flux and pressure measurements and the large model-form error encountered in a certain part of the flow.

42 ENGINEERING↗