Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interactive model”

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

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

At least 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Uncertainty characterization in a coupled human-natural system: Modeling agricultural adaptation in the Great Lakes Region

The Great Lakes Region's water quality and ecological health are threatened by the export of nutrients from agricultural lands, which causes eutrophication, hypoxia, and destructive algal blooms. The intensification of hydrologic cycles brought about by climate change is expected to exacerbate nutrient loading in the region, and, at the same time, agricultural adaptation to changing conditions is also expected to affect loading through shifting amounts and timing of fertilization. Quantifying these future effects and their interactions necessitates modeling both the human and natural processes as a coupled system, by pairing land use and agricultural management with hydrologic modeling. At the same time, compounding uncertainties arising from the complex interactions in both systems significantly limit our predictive understanding of the region's impacts. This study utilizes the Soil and Water Assessment Tool (SWAT), developed for simulating the impact of various farmer decisions on watershed functions in Western Lake Erie watersheds, and an under-development agent-based model (ABM) for agricultural management decisions. The aim of this study is to use global sensitivity analysis on the coupled ABM and SWAT models to quantify how uncertainty in both models interactively affects nutrient loading. To do so, we will conduct Sobol sensitivity analysis experiments at different levels of coupling assumptions to quantify how various uncertain factors (e.g., soil moisture and crop choice) and their interactions affect our estimates of nutrient loading. The results of this analysis will allow us to quantify how complex interactions and dependencies between both systems amplify the effect of uncertainties. Insights gained from this study will have broader implications for modeling the adaptive co-evolution of human and natural systems under climate change and can inform effective management of nutrient loading in the Great Lakes Region.

Climate Change↗
DOE OSTI2023

Sensitivity of Dislocation-GB interactions to simulation setups in atomistic models

Dislocation-grain boundary (GB) interactions play an integral role in strengthening of crystalline materials, and are dependent on external loading conditions and atomic arrangements at the grain boundary. While molecular dynamics (MD) simulations can provide critical insights into relationships between these parameters and the exact interaction, the outcomes are sensitive to the computational setup. Here, in this work, we explore the effect of computational setup (system size, GB orientation and boundary conditions) on stress-controlled dislocation-grain boundary (DGB) reactions using MD simulations for a well-studied copper symmetric $\langle$1 1 2$\rangle$ tilt boundary. The study demonstrates that the DGB reactions (pinning, absorption, transmission, etc.) are sensitive to the system size and the orientation of the GB to the applied stress. Additionally, the current work shows that there is a critical system size for the dislocation-grain boundary reaction stresses to converge, one which was found to be much larger than the setups in comparable MD studies. This was attributed to the specific stress states and the influence of surfaces, fixed or free, as these computational setups are modified. Finally, a stress-controlled setup with minimal model artifacts is examined, and the effect of coupled Schmid and non-Schmid stress components on the dislocation-GB reactions are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Nonlinear solution of classical three-wave interaction via finite-dimensional quantum model

The quantum three-wave interaction, the lowest-order nonlinear interaction in plasma physics, describes energy–momentum transfer between three resonant waves in the quantum regime. We describe how it may also act as a finite-degree-of-freedom approximation to the classical three-wave interaction in certain circumstances. By promoting the field variables to operators, we quantize the classical system, show that the quantum system has more free parameters than the classical system and explain how these parameters may be selected to optimize either initial or long-term correspondence. We then numerically compare the long-time quantum–classical correspondence far from the fixed point dynamics. We discuss the Poincaré recurrence of the system and the mitigation of quantum scrambling.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A pion-argon cross section measurement in the ProtoDUNE-SP experiment with cosmogenic muon

Neutrinos are tiny mysterious fundamental particles with small cross sections. Through neutrino physics, scientists across the world are trying to answer many intriguing questions about nature such as the dominance of matter over antimatter, CP violation in the lepton sector, number of supernovas in the early universe, etc. Detection of neutrinos requires massive particle detectors and intense neutrino beam owing to their small cross section. Deep Underground Neutrino Experiment (DUNE) is a next-generation neutrino experiment that is planned to start taking data beginning in 2026. DUNE will consist of 4 massive detectors, the first of which will be using single-phase liquid argon time projection chamber (LArTPC) technology. The ProtoDUNE-SP experiment is a prototype of the DUNE built at the CERN neutrino platform and uses the same detector technology that will be used in DUNE first module. The ProtoDUNE-SP experiment collected months of test beam and cosmic ray data beginning in September 2018. It was built to provide a testbed for the installation of detector parts for DUNE, showing long-term stability of the detector, understanding detector response for different test beam particles (including protons, pions, electrons, kaons, muons), and measurement of hadron-argon cross sections. When a particle passes through LArTPC electron-ion pairs are produced. To reconstruct the position and energy of a particle passing through the medium knowledge of ionization electron drift velocity is essential. The electron drift velocity is distorted by an excess positive charge built up in the detector, known as space charge. This study discusses a novel technique for measuring the ionization electron drift velocity using cosmic-ray muons. The technique uses tracks that travel the entire drift distance of the TPC for drift velocity determination. Secondly, the study discusses a method for converting the charge deposited into energy. The method is carried out in two step s. In th e first step detector response for energetic cosmic ray muons crossing the entire the TPC is used to make the charge deposition uniform throughout the TPC, and in the second step stopping cosmic-ray muons are used for determining the energy scale. Finally, the study discusses a pion-argon cross section measurement based on reweighting of Monte Carlo simulations using J. Calcutt's Geant4Reweight framework. Neutrinos cannot be directly detected; they are identified based on the interaction products. Pions are a common interaction product in a neutrino interaction. For precise modeling of neutrino event generators, it is essential to understand the pion-argon interaction. Pion-argon cross section measurement serves as an important input for neutrino interaction models. The results of the pion-argon total reaction cross section using the Geant4 reweighting technique are found to be in good agreement with Geant4 predictions. The many studies carried out in the ProtoDUNE-SP experi ment wil l be useful for current and future neutrino experiments using LArTPC technology including ICARUS, MicroBooNE, DUNE

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Modeling of aqueous species interaction energies prior to nucleation in cement-based gel systems

Arguably the most ubiquitous construction material in modern civilization, concrete is enabling the development of megacities around the globe together with increasing living standards in developing nations. However, the exact formation mechanisms of the strength-giving calcium-rich gels remain a topic of debate. Using density functional modeling, we simulate the fundamental solution-based building blocks of cement hydrates (calcium ions and silicate and aluminate monomers) and their propensity to form pair-wise complexes with bonding environments characteristic of those found in calcium-silicate-hydrate, calcium-alumino-silicate-hydrate and sodium-containing calcium-alumino-silicate-hydrate gels, as assessed from Gibbs free energies of chemical reactions. By accurately simulating the high pH pore solution chemistry in Portland cements and related systems, along with discrete solvation of the species, we hypothesize potential early age formation routes of the gels and discuss limitations and future work associated with this approach.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Analytical model for the motion and interaction of two-dimensional active nematic defects

Here, we develop an approximate, analytical model for the velocity of defects in active nematics by combining recent results for the velocity of topological defects in nematic liquid crystals with the flow field generated from individual defects in active nematics. Importantly, our model takes into account the long-range interactions between defects that result from the flows they produce as well as the orientational coupling between defects inherent in nematics. Our work complements previous studies of active nematic defect motion by introducing a linear approximation that allows us to treat defect interactions as two-body interactions and incorporates the hydrodynamic screening length as a tuning parameter. We show that the model can analytically predict bound states between two +1/2 winding number defects, effective attraction between two –1/2 defects, and the scaling of a critical unbinding length between ±1/2 defects with activity. The model also gives predictions for the trajectories of defects, such as the scattering of +1/2 defects by –1/2 defects at a critical impact parameter that depends on activity. In the presence of circular confinement, the model predicts a braiding motion for three +1/2 defects that was recently seen in experiments, as well as stable and ergodic trajectories for four or more defects.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Constraints on long-ranged interactions between dark matter and the Standard Model

Dark matter's existence is known thanks to its gravitational interaction with Standard Model particles, but it remains unknown whether this is the only force present between them. While many searches for such new interactions with dark matter focus on short-range, contact-like interactions, it is also possible that there exist weak, long-ranged forces between dark matter and the Standard Model. In this work, we present two types of constraints on such new interactions. First, we consider constraints arising from the fact that such a force would also induce long range interactions between Standard Model particles themselves, as well as between dark matter particles themselves. Combining the constraints on these individual forces generally sets the strongest constraints available on new Standard Model-dark matter interactions. Second, we consider the possibility of constraining new long-ranged interactions between dark matter and the Standard Model using the effects of dynamical friction in ultrafaint dwarf galaxies, especially Segue I. Such new interactions would accelerate the transfer of kinetic energy from stars to their surrounding dark matter, slowly reducing their orbits; the present-day stellar half-light radius of Segue I therefore allows us to exclude new forces which would have reduced stars' orbital radii below this scale by now.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

A multirate mass transfer model to represent the interaction of multicomponent biogeochemical processes between surface water and hyporheic zones (SWAT-MRMT-R 1.0)

Abstract. Surface water quality along river corridors can be modulated by hyporheic zones (HZs) that are ubiquitous and biogeochemically active. Watershed management practices often ignore the potentially important role of HZs as a natural reactor. To investigate the effect of hydrological exchange and biogeochemical processes on the fate of nutrients in surface water and HZs, a novel model, SWAT-MRMT-R, was developed coupling the Soil and Water Assessment Tool (SWAT) watershed model and the reaction module from a flow and reactive transport code (PFLOTRAN). SWAT-MRMT-R simulates concurrent nonlinear multicomponent biogeochemical reactions in both the channel water and its surrounding HZs, connecting the channel water and HZs through hyporheic exchanges using multirate mass transfer (MRMT) representation. Within the model, HZs are conceptualized as transient storage zones with distinguished exchange rates and residence times. The biogeochemical processes within HZs are different from those in the channel water. Hyporheic exchanges are modeled as multiple first-order mass transfers between the channel water and HZs. As a numerical example, SWAT-MRMT-R is applied to the Hanford Reach of the Columbia River, a large river in the United States, focusing on nitrate dynamics in the channel water. Major nitrate contaminants entering the Hanford Reach include those from the legacy waste, irrigation return flows (irrigation water that is not consumed by crops and runs off as point sources to the stream), and groundwater seepage resulting from irrigated agriculture. A two-step reaction sequence for denitrification and an aerobic respiration reaction is assumed to represent the biogeochemical transformations taking place within the HZs. The spatially variable hyporheic exchange rates and residence times in this example are estimated with the basin-scale Networks with EXchange and Subsurface Storage (NEXSS) model. Our simulation results show that (1), given a residence time distribution, how the exchange fluxes to HZs are approximated when using MRMT can significantly change the amount of nitrate consumption in HZs through denitrification and (2) source locations of nitrate have a different impact on surface water quality due to the spatially variable hyporheic exchanges.

58 GEOSCIENCES↗
DOE OSTI2025

Electromagnetic moments of 215,217 Bi: Probing shell evolution beyond N = 126

The nuclear properties of bismuth isotopes (Z = 83) , with just one valence proton above the closed spherical shell at (Z = 82) , are expected to be governed by a single unpaired proton. However, already in semimagic 209 Bi (Z = 83, N = 126) , , the magnetic moment (μ) strongly deviates from the single-particle Schmidt value. A near linear decrease in μ with the increase of N after the N = 126 magic number was observed up to N = 130 . In order to test whether this trend is kept at N > 130 and to reveal the underlying mechanisms, an investigation of 215,217 Bi (N = 132, 134) has been undertaken. The magnetic dipole and electric quadrupole moments of the I π = 9/2 - nuclear ground states in these isotopes have been measured for the first time using the in-source resonance-ionization spectroscopy technique at ISOLDE (CERN). It has been shown that the linearly decreasing trend of μ( 209,211,213 Bi g ) is broken in 215,217 Bi with a nearly constant value of μ observed. Experimental data have been compared to calculations in the framework of the configuration-interaction shell model with the monopole-based universal V MU +LS interaction. The peculiarities in the behavior of μ(Bi, 9/2 - ) with increasing neutron number are explained as being due to the shell evolution, change of the neutron orbitals occupancies and strong configuration mixing beyond N = 130 . Also, the difference in the μ trends for bismuth (Z = 83) and astatine (Z = 85) isotopes with N > 126 are reproduced by the shell-model calculations. It is shown that monopole interaction plays noticeable role in the description of the peculiarities of the μ behaviour. Additionally, the extension of the application of the V MU interaction to the μ isotopic trends for heavy nuclei is important for further study of the capabilities of this promising version of the shell-model calculations.

Dipole magnetic moments↗
DOE OSTI2023

Phase-field modeling of dislocation–interstitial interactions

The mechanical behavior of body-centered cubic (BCC) materials can be dramatically affected by the presence of interstitial solute atoms. Here we present a new phase-field dislocation dynamics formulation to include the diffusion of interstitials. Short-range interactions are accounted for by a concentration-dependent lattice energy, and long-range interactions are accounted for by modifications to the elastic energy. The interstitial diffusion law introduces gradients that require methods for minimizing Gibbs oscillations, which is done via a modified Green’s function. The formulation is general to any solute-solvent system and is applied here to Nb-O as a model system, whose interstitial parameters are obtained from ab initio calculations. The effect of O on the core structures of Nb edge and screw dislocations is calculated. The diffusion of O to form interstitial atmospheres around dislocation cores is simulated, as well as the critical stresses required for dislocations to break away or cross slip from these atmospheres. Future applications of the method to simulate complex interstitial embrittlement mechanisms are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A detailed study of interpretability of deep neural network based top taggers

Abstract Recent developments in the methods of explainable artificial intelligence (XAI) allow researchers to explore the inner workings of deep neural networks (DNNs), revealing crucial information about input–output relationships and realizing how data connects with machine learning models. In this paper we explore interpretability of DNN models designed to identify jets coming from top quark decay in high energy proton–proton collisions at the Large Hadron Collider. We review a subset of existing top tagger models and explore different quantitative methods to identify which features play the most important roles in identifying the top jets. We also investigate how and why feature importance varies across different XAI metrics, how correlations among features impact their explainability, and how latent space representations encode information as well as correlate with physically meaningful quantities. Our studies uncover some major pitfalls of existing XAI methods and illustrate how they can be overcome to obtain consistent and meaningful interpretation of these models. We additionally illustrate the activity of hidden layers as neural activation pattern diagrams and demonstrate how they can be used to understand how DNNs relay information across the layers and how this understanding can help to make such models significantly simpler by allowing effective model reoptimization and hyperparameter tuning. These studies not only facilitate a methodological approach to interpreting models but also unveil new insights about what these models learn. Incorporating these observations into augmented model design, we propose the particle flow interaction network model and demonstrate how interpretability-inspired model augmentation can improve top tagging performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Evidence for triaxial shape coexistence in 74 Ge

The deformation properties of the low-lying states in 74 Ge have been investigated using multistep Coulomb excitation. The measurements were carried out with the advanced 𝛾-ray tracking array, GRETINA, and the CHICO2 particle detector. A comprehensive set of 𝐸⁢2 transition and diagonal matrix elements was deduced following an analysis with the semiclassical coupled-channels code GOSIA. The data were compared with results of calculations carried out within the framework of the generalized triaxial rotor model as well as with the configuration interaction shell model and the symmetric rotor model. Results from calculations with covariant density functional theory were used to construct a five-dimensional collective Hamiltonian for further comparisons with the data. Collectively, the calculations provide an accurate reproduction of the experimental matrix elements and further support an understanding in terms of the coexistence of two axially asymmetric shapes. In conclusion, this leads to an overall interpretation of the underlying structure of 74 Ge requiring triaxiality, as is also the case in the neighboring even-mass Ge isotopes.

59 ≤ A ≤ 89↗
DOE OSTI2025

Drought adaptation index (DAI) based on BLUP as a selection approach for drought-resilient switchgrass germplasm

This study introduces a Drought Adaptation Index (DAI), derived from Best Linear Unbiased Prediction (BLUP), as a method to assess drought resilience in switchgrass (Panicum virgatum L.). A panel of 404 genotypes was evaluated under drought-stressed (CV) and well-watered (UC) conditions over four consecutive years (2019–2022). BLUP-estimated biomass yields were used to calculate the DAI, which enabled classification of genotypes into four adaptation groups: very well-adapted, well-adapted, adapted, and unadapted. The DAI was compared with conventional drought tolerance indices, including the Stress Susceptibility Index (SSI), Stress Tolerance Index (STI), Geometric Mean Productivity (GMP), and Yield Stability Index (YSI). Correlation analyses demonstrated strong agreement between DAI and these indices, supporting its validity and consistency. Biplot analyses using the Genotype plus Genotype-by-Environment Interaction (GGE) and Additive Main Effects and Multiplicative Interaction (AMMI) models revealed significant genotype-by-environment interactions (GEI) and identified J222.A, J463.A, and J295.A. A as high-performing genotypes, with J222.A exhibiting greater yield stability across treatments and years. Additionally, DAI isoline curves provided a graphical representation of differential genotype performance under drought and control conditions. These visualizations aided in distinguishing genotypes with stable and superior biomass yield across contrasting environments. Overall, the BLUP-based DAI is a robust and practical selection tool that improves the accuracy of identifying drought-resilient, high-yielding switchgrass genotypes. Its integration into breeding programs offers a comprehensive framework for improving biomass productivity and stress adaptation under variable climatic conditions. The application of DAI supports the development of climate-resilient cultivars and contributes to sustainable bioenergy and forage production systems.

BLUP↗
DOE OSTI2020

Reconciling Contrasting Relationships Between Relative Dispersion and Volume-Mean Radius of Cloud Droplet Size Distributions

Cloud droplet spectral relative dispersion is critical to parameterizations of cloud radiative properties, warm-rain initiation, and aerosol-cloud interactions in models; however, there is no consistent relationship between relative dispersion and volume-mean radius in literature, which hinders improving relative dispersion parameterization and calls for physical explanation. Here we show, by analyzing aircraft observations of cumulus clouds during Routine AAF [Atmospheric Radiation Measurement (ARM) Aerial Facility] Clouds with Low Optical Water Depths (CLOWD) Optical Radiative Observations, that the correlation between relative dispersion and volume-mean radius changes from positive to negative as volume-mean radius increases. With the new observation, we postulate that the sign of the correlation is determined by whether or not condensation (evaporation) occurs simultaneously with significant new activation (deactivation). The hypothesis is validated by simulations of both an adiabatic cloud parcel model and a parcel model accounting for entrainment-mixing. A new quantity, first bin strength, is introduced to quantify this new observation. Theoretical analysis of truncated gamma and modified gamma size distributions further supports the hypothesis and reconciles the contrasting relationships between relative dispersion and volume-mean radius, including the results in polluted fog observations. The results could shed new light on the so-called “twilight zone” between cloudy and cloud-free air, which in turn affects evaluation of aerosol-cloud interactions and retrieval of aerosol optical depth.

Lu, Chunsong↗
DOE OSTI2024

Integrating Tide‐Driven Wetland Soil Redox and Biogeochemical Interactions Into a Land Surface Model

Abstract Redox processes, aqueous and solid‐phase chemistry, and pH dynamics are key drivers of subsurface biogeochemical cycling and methanogenesis in terrestrial and wetland ecosystems but are typically not included in terrestrial carbon cycle models. These omissions may introduce errors when simulating systems where redox interactions and pH fluctuations are important, such as wetlands where saturation of soils can produce anoxic conditions and coastal systems where sulfate inputs from seawater can influence biogeochemistry. Integrating cycling of redox‐sensitive elements could therefore allow models to better represent key elements of carbon cycling and greenhouse gas production. We describe a model framework that couples the Energy Exascale Earth System Model (E3SM) Land Model (ELM) with PFLOTRAN biogeochemistry, allowing geochemical processes and redox interactions to be integrated with land surface model simulations. We implemented a reaction network including aerobic decomposition, fermentation, sulfate reduction, sulfide oxidation, methanogenesis, and methanotrophy as well as pH dynamics along with iron oxide and iron sulfide mineral precipitation and dissolution. We simulated biogeochemical cycling in tidal wetlands subject to either saltwater or freshwater inputs driven by tidal hydrological dynamics. In simulations with saltwater tidal inputs, sulfate reduction led to accumulation of sulfide, higher dissolved inorganic carbon concentrations, lower dissolved organic carbon concentrations, and lower methane emissions than simulations with freshwater tidal inputs. Model simulations compared well with measured porewater concentrations and surface gas emissions from coastal wetlands in the Northeastern United States. These results demonstrate how simulating geochemical reaction networks can improve land surface model simulations of subsurface biogeochemistry and carbon cycling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A crystal plasticity model with an atomistically informed description of grain boundary sliding for improved predictions of deformation fields

Crystal plasticity (CP) is a powerful meso-scale technique for deformation modeling in polycrystalline materials. Crystal plasticity models typically do not include an explicit description of grain boundary (GB) sliding, which could lead to inaccurate predictions of strain distributions, especially in the vicinity of GBs. In the present study, a CP model is developed that includes GB sliding as an additional deformation mechanism and models the interaction between slip and GB sliding. Atomistic simulations are used to formulate the constitutive model for GB sliding and its interaction with incident slip. The deformation fields and structure of the GBs are obtained directly from experimental characterization and are faithfully reproduced as bicrystal systems for molecular dynamics simulations. The GBs are modeled as random-type (non-coincidence, asymmetrical) boundaries as observed from experimental data. The underlying atomistic-scale mechanism for pure sliding of random GBs was found to be analogous to fluid-flow. The interaction between incident slip and a sliding GB caused local increases in stress concentration, which further led to a momentary increase in the local sliding rate. The displacement profiles at the sliding GBs computed from the CP model with sliding accommodation are 38% more accurate than the baseline CP model as quantified by the mean squared error between the simulations and experiment. Here these results help improve the sophistication and accuracy of deformation modeling by including physics-based descriptions for GB sliding and its interaction with slip, eventually leading to more reliable predictions of micromechanical quantities.

36 MATERIALS SCIENCE↗