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 199 records · Page 11

Accuracy Enhancement of Nuclear Power Plant Simulators Utilizing High Accuracy Simulation Predictions

More recently, reactor core simulators for core designs associated with commercial nuclear power plants that utilize what is believed to be higher fidelity models have been developed. Features such as neutronics models that utilize transport equation solvers with fine spatial meshes and many energy-groups, thermal-hydraulic models that utilize sub-channel solvers with fine spatial mesh and capable of treating a wide range of fluid conditions, and fuel-coolant chemistry interaction models capable of treating CRUD deposition are to be found in these higher fidelity core simulators. These reactor core simulators require access to higher performance computers, characterized by many processors, cores and large memory. So associated with utilization of these simulators is access to high performance computers and ability to accommodate in one’s workflow longer execution times. By contrast, currently used core simulators by the nuclear industry can execute on engineering workstations and have execution times of seconds to minutes. The desirability for having short execution times is not only desired for support of time critical tasks but supports the mental process of decision making by engineers. The goal of the work reported upon here has the objective of retaining the fidelity of higher fidelity models while retaining the ability to utilize engineering workstations. Beyond the core simulator goal, additional goals of this work include incorporating the just described core simulator capability into a Nuclear Steam Supply System (NSSS) simulator, and to incorporate the resulting capability into an environment supportive of design and operational decision making associated with nuclear power stations. The model selected for the core neutronics model is the NESTLE code, for the core thermal-hydraulic model is the CTF code utilizing coarse mesh, and for the NSSS model is the RELAP5-3D code. WSC’s proprietary 3KEYMASTERTM platform is being used to provide software coupling, user interface, visualization, and reporting. The NESTLE core neutronics simulator was first integrated with the CTF core thermal-hydraulic simulator using CTF developed communication commands which are also used for CTF to communicate with RELAP5-3D under WSC’s proprietary 3KEYMASTERTM platform. To assure NESTLE prediction consistency with higher fidelity core neutronic simulators, buffer codes have been created to automatically generate from output files written by the VERA core simulator the NESTLE nodal neutronic parameter’ library, geometry, and pin-power reconstruction input files, thereby avoiding a number of challenges associated with utilizing lattice physics codes and providing consistency with VERA predictions. To treat absorber rod effects a multi-set library is utilized, where a set refers to a specific absorber rod fully inserted pattern. A coarse spatial mesh CTF model was developed with features added that support using CTF as envisioned in the engineering quality simulator. A hybrid meshing approach was implemented to allow for automated construction of models with mixed levels of refinement. Specifically, a core model could resolve some assemblies at a nodal level (4 subchannels per assembly) and others at a pin-resolution (one subchannel per coolant subchannel in the assembly). The intention is that this will allow for better resolution of limiting conditions such as DNBR and PCT, which are based on local rod and subchannel conditions. Further development was done of features that enhance the capabilities for the envisioned engineering quality simulator that has been developed, but now for RELAP-3D. The RELAP5-3D code development includes ability to model more than 999 components and the addition of the cross-channels turbulence mixing model and the void drift model that are implemented in CTF, aiming to achieve closer prediction agreement of the two codes for transient simulations, specifically, more accurate matches of the overall mass, momentum, and energy exchanges of both the liquid and gas phases between the neighboring core assemblies. Graphics were also developed for the Instructor Station for this project under WSC’s proprietary 3KEYMASTERTM platform to facilitate design and operational decision making.

42 ENGINEERING↗

System Analysis Modeling and Intermodal Transportation for Commercial Spent Nuclear Fuel

The United States Department of Energy (DOE) Office of Nuclear Energy is applying knowledge and understanding in the areas of systems level engineering, analysis, and decision making to better inform the waste management pathways for U.S. spent nuclear fuel (SNF). Currently there are 93 operating and 23 shutdown commercial nuclear reactors in the United States. SNF at most of these locations is being stored in spent fuel pools, dry storage, or both [1, 2]. This paper initially goes over the basics of the agent-based simulation tool known as the Next Generation System Analysis Model (NGSAM) that has the capability to model interaction and movements of individual components or groups like casks, railcars etc. The next section covers some of the methods implemented in the Java Transportation Operations Model (JTOM) to facilitate the movement of assemblies, casks, railcars etc. and the various intermodal transfer options. This section also covers the current cask loading and intermodal transfer times implemented in NGSAM. The last section goes over the proposed values for intermodal transportation, and cask transfer/loading operations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

High‐Speed Embedded Ink Writing of Anatomic‐Size Organ Constructs

Embedded ink writing (EIW) is an emerging 3D printing technique that fabricates complex 3D structures from various biomaterial inks but is limited to a printing speed of ∼10 mm s −1 due to suboptimal rheological properties of particulate-dominated yield-stress fluids when used as liquid baths. In this work, a particle-hydrogel interactive system to design advanced baths with enhanced yield stress and extended thixotropic response time for realizing high-speed EIW is developed. In this system, the interactions between particle additive and three representative polymeric hydrogels enable the resulting nanocomposites to demonstrate different rheological behaviors. Accordingly, the interaction models for the nanocomposites are established, which are subsequently validated by macroscale rheological measurements and advanced microstructure characterization techniques. Filament formation mechanisms in the particle-hydrogel interactive baths are comprehensively investigated at high printing speeds. To demonstrate the effectiveness of the proposed high-speed EIW method, an anatomic-size human kidney construct is successfully printed at 110 mm s −1 , which only takes ∼4 h. This work breaks the printing speed barrier in current EIW and propels the maximum printing speed by at least 10 times, providing an efficient and promising solution for organ reconstruction in the future.

36 MATERIALS SCIENCE↗

Modeling Plant Nutrient Acquisition Strategies Alters Projections of Carbon and Nitrogen Dynamics in Bioenergy Agroecosystems

Plant strategies to acquire nutrients from limited environments help shape ecosystem carbon (C) and nitrogen (N) cycling and response to environmental change. The effects of plant strategies on ecosystem dynamics are largely uncharacterized in bioenergy agroecosystems, where the impacts could determine bioenergy's ability to meet its sustainability goals of storing C and reducing N loss. We used FUN-BioCROP (Fixation and Uptake of Nitrogen-Bioenergy Carbon, Rhizosphere, Organisms and Protection), a plant–microbe interaction model of coupled plant nutrient uptake and soil organic matter decomposition, to simulate the effects of nutrient acquisition strategies on soil microbial activity and ecosystem nutrient cycling in bioenergy feedstocks miscanthus (Miscanthus × giganteus) and sorghum (Sorghum bicolor (L.) Moench). We examined the model's ability to reproduce the relative effects of belowground nutrient uptake on microbial activity using a reanalysis of empirical data showing that miscanthus root exudation provoked a larger soil microbial response than sorghum. From baseline model simulations, we found that the ability of miscanthus to retranslocate N resulted in higher N uptake at a lower C cost than the sorghum/soybean rotation and that soil C and N pools increased under perennial (miscanthus) and decreased under annual (sorghum/soybean) cultivation. The model also predicted that greater root exudation increased soil C accumulation, highlighting the role of roots in forming stable soil C. Overall, the baseline model was unable to reproduce field observations of miscanthus root exudation stimulating microbial activity more than sorghum. To improve the model, we updated the soil microbial parameters in miscanthus to have faster decomposition, a higher C/N ratio, and greater carbon use efficiency. These changes improved the simulated soil microbial response to miscanthus root exudation, supporting the hypothesis that miscanthus soils foster a microbial community that is more responsive to root exudation than that of sorghum.

09 BIOMASS FUELS↗

Cluster-Graph Fingerprinting: A Framework for Quantitative Analysis of Machine-Learned Interatomic Model Training and Simulation Data

Machine-learned interatomic models represent a significant advancement in simulation methods, extending the predictive ability of first-principles methods to previously inaccessible length and time scales. However, the data-driven nature of these models can lead to difficult-to-detect errors that can compromise prediction accuracy. To address this challenge, we introduce a novel fingerprinting approach based on the Chebyshev Interaction Model for Efficient Simulation (ChIMES) ML-IAM graph-based descriptor. Our strategy enables efficient and statistically rigorous analysis of system configurations used in ML-IAM training and those generated by their application, e.g., in molecular dynamics simulations. We demonstrate that these fingerprints can effectively assess novelty of a configuration relative to an existing data set and determine dissimilarity among individual configurations, which are two key tasks in workflows for active learning-based ML-IAM training, data set curation, and on-the-fly uncertainty quantification.

36 MATERIALS SCIENCE↗

Muon Neutrino Charged Current Single Pion Production on Various Targets in the MINERvA Detector

Measurements of neutrino oscillation parameters at accelerator-based experiments require precise knowledge of neutrino interaction cross sections. A significant fraction of neutrino interactions observed in these experiments produce single charged pions. While single charged pion production has been measured on hydrocarbon and water, few precise measurements exist for heavier targets, and there are no high statistics measurements that directly compare pion production on different targets. This thesis measures differential cross sections of muon neutrino charged current single charged pion production with no other mesons in the final state in a neutrino beam with a mean energy of $\approx 6$ GeV. The pion production is studied on targets of scintillator, carbon, water, iron, and lead at MINERvA, a neutrino cross section experiment located in the NuMI beamline at Fermilab. These measured cross sections are differential in $\mu^-$ and $\pi^+$ kinematics. A strong suppression of low $Q^2 $ and an enhancement of low $T_\pi$ are observed relative to commonly used neutrino interaction models in both the light and heavy nuclear targets. Incorporation of these observed effects in future neutrino oscillation experiments will improve their ability to measure energy of neutrino interactions, and therefore measure the parameters governing neutrino oscillations.

Bercellie, Aaron↗

Muon neutrino charged current single pion production on various targets in the MINERvA detector

Measurements of neutrino oscillation parameters at accelerator-based experiments require precise knowledge of neutrino interaction cross sections. A significant fraction of neutrino interactions observed in these experiments produce single charged pions. While single charged pion production has been measured on hydrocarbon and water, few precise measurements exist for heavier targets, and there are no high statistics measurements that directly compare pion production on different targets. This thesis measures differential cross sections of muon neutrino charged current single charged pion production with no other mesons in the final state in a neutrinobeam with a mean energy of ≈ 6 GeV. The pion production is studied on targets of scintillator, carbon, water, iron, and lead at MINERvA, a neutrino cross section experiment located in the NuMI beamline at Fermilab. These measured cross sections are differential in μ - and π + kinematics. A strong suppression of low Q 2 and an enhancement of low π are observed relative to commonly used neutrino interaction models in both the light and heavy nuclear targets. Incorporation of these observed effects in future neutrino oscillation experiments will improve their ability to measure energy of neutrino interactions, and therefore measure the parameters governing neutrino oscillations.

Bercellie, Aaron↗

Application of Artificial Intelligence/Machine Learning to Operations Research

This report examines the transformative impact of Artificial Intelligence (AI) and Machine Learning (ML) on operations research, private industry, and government sectors, highlighting their applications in automating processes, enhancing decision-making, and optimizing complex systems. AI/ML technologies have revolutionized industries through predictive maintenance, supply chain optimization, and autonomous systems, while also advancing public safety and defense operations. However, challenges such as data integrity, model transparency, and the need for human oversight persist, particularly in high-consequence environments. The report emphasizes the critical role of explainable AI (XAI) and human-computer interaction models like Human-in-the-Loop (HITL) and Human-on-the-Loop (HOTL) in fostering trust and accountability. Balancing automation with ethical responsibility and transparency is essential for the continued successful integration of AI/ML into operational and strategic decision-making frameworks.

97 MATHEMATICS AND COMPUTING↗

Effects of EMIC Wave-Driven Proton Precipitation on the Ionosphere

In this report we investigate the proton precipitation caused by electromagnetic ion cyclotron (EMIC) waves and its impact on the ionosphere with the self-consistent ring current-atmosphere interactions model (RAM-SCBE) coupled with an ionospheric particle transport model Global Airglow (GLOW). The EMIC wave diffusion process causes significant precipitation of tens of keV protons, particularly in the afternoon to midnight sector. These precipitating energetic protons further impact the ionosphere-thermosphere and contribute to the ionization in the E/F regions. The integrated auroral conductance is significantly enhanced in the dusk-to-midnight sector. Although the EMIC waves do not directly interact with the ring current electrons, after the EMIC wave scattering included in the model, remarkable changes are found in the global distribution of precipitating electron flux. This means that the addition of proton precipitation driven by EMIC waves results in feedback effects on the ring current electron dynamics through the circulation system. Validation is also conducted by comparing the simulated precipitation flux and ionospheric electron density with observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Methodology of the Extraction of Multi-Differential Cross Sections of Charged-Current $\nu_\mu$-Argon Interactions in MicroBooNE Using Wire-Cell $\nu_\mu$CC Selection

Neutrino physics experiments rely on accurate nuclear-interaction models that are in part guided by experimental observations. Multi-differential cross section measurements are particularly valuable for exploring the underlying physics described in these models, making them a significant step forward in the field. This note presents methodology in preparation to unfold neutrino flux-averaged double and triple-differential cross sections of the inclusive muon neutrino charged-current interaction on Argon. This work builds upon the existing framework and methodology presented in recent energy dependent cross section measurements at MicroBooNE using the Wire-Cell tomographic event reconstruction. The signal definition, choice of binning, and handling of estimated detector uncertainties are discussed, and the MicroBooNE simulation model is validated over the multi-dimensional phase space of muon momentum, muon polar angle, and visible hadronic energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Magnetic structure and exchange interactions in the Heisenberg pyrochlore antiferromagnet Gd 2 Pt 2 O 7

In this report the Heisenberg pyrochlore antiferromagnet Gd 2 Pt 2 O 7 is one of a series of gadolinium pyrochlore compounds with a variety of B-site cations. Despite the expected simplicity of a spin-only Gd 3+ Heisenberg interaction model, the gadolinium pyrochlore series exhibits various complex magnetic ground states at low temperature. Gd 2 Pt 2 O 7 displays the highest temperature magnetic order of the series with $T_N$=1.6 K, which has been attributed to enhanced superexchange pathways facilitated by empty $5d$ $e_g$ Pt orbitals. In this study, we use various neutron scattering techniques on an isotopically enriched polycrystalline 160 Gd 2 Pt 2 O 7 sample to examine the magnetic structure and spin-wave excitation spectrum below $T_N$ in order to extract the dominant exchange interactions. We find that the ground-state magnetic structure is the Palmer-Chalker state previously seen in Gd 2 Sn 2 O 7 with an associated gapped excitation spectrum consistent with enhanced exchange interactions between further near-neighbor Gd 3+ ions. We confirm this exchange model with analysis of the magnetic diffuse scattering in the paramagnetic regime using polarized neutrons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Progress toward double-differential cross-section measurements of single charged-pion production in charged current muon neutrino interactions on argon with SBND

Resonant neutrino–nucleus processes constitute a significant portion of neutrino interactions in the few-GeV energy region. Charged-current interactions with a muon and one single charged pion in the final state (CC1$\pi$) are primarily sensitive to resonant processes, while also receiving contributions from non-resonant processes and being strongly affected by nuclear effects and final-state interactions. A precise understanding of this channel is essential for improving neutrino interaction models and reducing systematic uncertainties in oscillation measurements. The Short-Baseline Near Detector (SBND) at Fermilab has collected the largest neutrino–argon dataset to date, providing an excellent opportunity to study CC1$\pi$ interactions on argon. This analysis aims to improve existing single-differential cross-section measurements by reporting results in a variety of kinematic variables and to perform the first double-differential CC1$\pi$ cross-section measurement on argon, using novel techniques for pion energy reconstruction.

Pelegrina-Gutiérrez, Luis [Granada U.] (ORCID:0000↗

The Spectacular Ultraviolet Flash from the Peculiar Type Ia Supernova 2019yvq

Early observations of Type Ia supernovae (SNe Ia) provide essential clues for understanding the progenitor system that gave rise to the terminal thermonuclear explosion. Here, we present exquisite observations of SN 2019yvq, the second observed SN Ia, after iPTF 14atg, to display an early flash of emission in the ultraviolet (UV) and optical. Our analysis finds that SN 2019yvq was unusual, even when ignoring the initial flash, in that it was moderately underluminous for a SN Ia ($M_g \approx -18.5$ mag at peak) yet featured very high absorption velocities ($v \approx 15,0000$ km s -1 for Si II λ6355 at peak). We find that many of the observational features of SN 2019yvq, aside from the flash, can be explained if the explosive yield of radioactive 56 Ni is relatively low (we measure $M_{^{56}\mathrm{Ni}} = 0.31 \pm 0.05\,M_\odot$) and it and other iron-group elements are concentrated in the innermost layers of the ejecta. To explain both the UV/optical flash and peak properties of SN 2019yvq we consider four different models: interaction between the SN ejecta and a nondegenerate companion, extended clumps of 56 Ni in the outer ejecta, a double-detonation explosion, and the violent merger of two white dwarfs. Each of these models has shortcomings when compared to the observations; it is clear additional tuning is required to better match SN 2019yvq. In closing, we predict that the nebular spectra of SN 2019yvq will feature either H or He emission, if the ejecta collided with a companion, strong [Ca II] emission, if it was a double detonation, or narrow [O I] emission, if it was due to a violent merger.

79 ASTRONOMY AND ASTROPHYSICS↗

Extracting the number of short-range correlated nucleon pairs from inclusive electron scattering data

The extraction of the relative abundances of short-range correlated (SRC) nucleon pairs from inclusive electron scattering is studied using the generalized contact formalism (GCF) with several nuclear interaction models. GCF calculations can reproduce the observed scaling of the cross-section ratios for nuclei relative to deuterium at high xB and large Q 2 , a 2 =(σ A /A)/(σ d /2). In the nonrelativistic instant-form formulation, the calculation is very sensitive to the model parameters and only reproduces the data using parameters that are inconsistent with ab initio many-body calculations. Using a light-cone GCF formulation significantly decreases this sensitivity and improves the agreement with ab initio calculations. The ratio of similar mass isotopes, such as 40 Ca and 48 Ca, should be sensitive to the nuclear asymmetry dependence of SRCs, but is found to also be sensitive to low-energy nuclear structure. Thus the empirical association of SRC pair abundances with the measured a2 values is only accurate to about 20%. Finally, improving this will require cross-section calculations that reproduce the data while properly accounting for both nuclear structure and relativistic effects.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Probabilistic modeling approach for tracking each collision and breakup of computational parcels

A modeling approach for binary collision and particle breakup is presented for a stochastic particle model in which a computational parcel represents a number of similar physical particles. Presented in this report includes sampling from a distribution of discrete data, parcel number reduction, probabilistic calculation of collision rates, and various particle-particle interactions. Modeling of elastic scattering, coalescence and agglomeration, and breakup are discussed for solid-solid, solid-liquid, and liquid-liquid particle collisions. Numerical procedures are also summarized.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Scalable simulation of coupled adsorption and transport of methane in confined complex porous media with density preconditioning

The growing significance of shales and tight formations in the transition to less carbon-intensive and clean energy drives the research endeavor to understand the physics of gas flow within these systems. However, shales are composed of massively heterogeneous physical and chemical features. Most nano-sized pores connect to millimeter-scale fractures, leading to multiscale transport. These nano-scale pore throats demonstrate non-classical flow behavior, such as non-negligible slip velocities and adsorbed gas layers at the boundary. As a result, classical computational fluid dynamics models do not capture the physics. In this work, we develop a coupling scheme for the multiple-relaxation-time (MRT) lattice Boltzmann (LB) method that integrates the Peng-Robinson equation of state into a pseudo-potential interaction model to capture the physics of methane flow in irregular networks of channels that represent nano-scale porous media. We use atomistic simulations to calibrate and validate our model in slit nano-channels. We propose a preconditioning scheme to initialize the coupled transport and adsorption simulation of methane in complex porous media. The results of this implementation of LB agree with Direct Simulation Monte Carlo (DSMC) and Molecular Dynamics (MD) simulations. We then scale up the LB implementation through vectorization and indirect addressing. We parallelize it using Message Passing Interface (MPI) and OpenMP frameworks to simulate transport and adsorption in complex media with a million lattices. Additionally, we analyze the differences between coupled and transport-only simulations in two case studies and show that considering phase behavior, i.e., adsorption, can significantly change the flow behavior. This work constitutes an important step towards bridging the gap between molecular flow and system-scale behavior of complex disordered porous media.

42 ENGINEERING↗

A reactive molecular dynamics model for uranium/hydrogen containing systems

Uranium-based materials are valuable assets in the energy, medical, and military industries. However, understanding their sensitivity to hydrogen embrittlement is particularly challenging due to the toxicity of uranium and the computationally expensive nature of quantum-based methods generally required to study such processes. In this regard, we have developed a Chebyshev Interaction Model for Efficient Simulation (ChIMES) that can be employed to compute energies and forces of U and UH3 bulk structures with vacancies and hydrogen interstitials with accuracy similar to that of Density Functional Theory (DFT) while yielding linear scaling and orders of magnitude improvement in computational efficiency. Here, we show that the bulk structural parameters, uranium and hydrogen vacancy formation energies, and diffusion barriers predicted by the ChIMES potential are in strong agreement with the reference DFT data. We then use ChIMES to conduct molecular dynamics simulations of the temperature-dependent diffusion of a hydrogen interstitial and determine the corresponding diffusion activation energy. Our model has particular significance in studies of actinides and other high-Z materials, where there is a strong need for computationally efficient methods to bridge length and time scales between experiments and quantum theory.

36 MATERIALS SCIENCE↗

Continuum contribution to charged-current absorption of low-energy $ν_e$ on $^{40}$Ar

Accurate modeling of the absorption of tens-of-MeV $ν_e$ on $^{40}$Ar is needed to enable measurements of astrophysical neutrinos using large liquid argon time projection chamber (LArTPC) detectors, such as those planned for the Deep Underground Neutrino Experiment (DUNE). We revisit the MARLEY neutrino interaction model used in present estimates of DUNE sensitivity to supernova and solar neutrino signals. Multiple theoretical refinements are pursued, especially in the unbound continuum region of nuclear excitation energy. Inclusive charged-current neutrino-argon cross sections are calculated using a hybrid strategy. Nuclear transitions to unbound states are treated using a Hartree-Fock Continuum Random Phase Approximation (HF-CRPA) model, including forbidden contributions. Allowed transitions to low-lying discrete levels are also included using indirect measurements and approximate corrections for the momentum transfer dependence. Exclusive predictions are obtained by coupling these calculations with a statistical nuclear de-excitation model. The impact on observables of interest for DUNE and similar experiments is examined in terms of both total and differential cross sections. Our refined calculations predict a lower allowed portion of the cross section relative to the prior MARLEY model. At neutrino energies appreciably below 100 MeV, the inclusion of forbidden transitions does not fully compensate for the loss of allowed strength. For a representative neutrino burst from a galactic core-collapse supernova, our results suggest that MARLEY 1.2.0 overestimates the event yield in a DUNE-like detector by approximately 20%. However, because this overestimation is more severe at backwards angles, use of the charged-current $ν_e$-$^{40}$Ar reaction for supernova pointing may be more feasible than previously expected.

Gardiner, Steven [Fermilab]↗