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 91 records · Page 5

Conceptual Design of Integrated Energy Systems with Market Interaction Surrogate Models

Most integrated energy system (IES) optimization frameworks employ the price-taker approximation, which ignores important interactions with market and can result in overestimated economic values. In this work, we pro-pose a machine learning surrogate-assisted optimization framework to quantify the IES/market interactions and thus go beyond price taker. We use time series clustering to generate representative IES operation profiles for the IES optimization problem and use machine learning surrogate models to predict the IES/market interaction. We quantify the accuracy of the time series clustering and surrogate models in a case study to optimally retrofit a nuclear power plant with polymer electrolyte membrane electrolyzer to co-produce electricity and hydrogen.

Chen, Xinhe↗

Conceptual Design of Integrated Energy Systems with Market Interaction Surrogate Models

Most integrated energy system (IES) optimization frameworks employ the price-taker approximation, which ignores important interactions with the market and can result in overestimated economic values. In this work, we propose a machine learning surrogate-assisted optimization framework to quantify IES/market interactions and thus go beyond price-taker. We use time series clustering to generate representative IES operation profiles for the optimization problem and use machine learning surrogate models to predict the IES/market interaction. We quantify the accuracy of the time series clustering and surrogate models in a case study to optimally retrofit a nuclear power plant with a polymer electrolyte membrane electrolyzer to co-produce electricity and hydrogen.

Chen, Xinhe↗

Characteristics of Fluid‐Solid Interaction Constitutive Models Within Poroelastodynamics at Higher Strain‐Rates and Large Deformations Implemented in 1D

The large deformation, mixed formulation, finite element (FE) modeling approach presented in Irwin et al. 2024 is extended herein to include improved constitutive models for representing dynamic solid-fluid interactions at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger overpressure magnitudes (𝒪⁡(1⁢0 2 )⁢kPa) within a biphasic soft porous material using Theory of Porous Media (TPM) at finite strain. Specifically, these constitutive modeling improvements are the following: (i) a more physically robust constitutive model for pore fluid seepage velocity via inclusion of pore fluid viscous stress, and (ii) a modified deformation-dependent-permeability model and updated hyperelastic constitutive model better suited for handling larger volumetric compressions and extensions. The novelty of the present work is mainly the contribution (i): inclusion of pore fluid viscous stress at higher strain-rate and large deformations, which requires 𝐶 1 continuity in the weak formulation, accomplished by employing Hermite cubic interpolation functions within a mixed nonlinear poromechanical finite element formulation. In (ii), the model is updated to weakly enforce solid phase incompressibility, such that this assumption is not violated numerically, which provides improved numerical stability for achieving larger overpressure magnitudes on 𝒪⁡(1⁢0 2 ) kPa, which were not achievable with the previous Kozeny–Carman model in Irwin et al. 2024. Also in (ii), the volumetric part of the solid skeleton free energy function is modified to ensure proper bounds on the solid skeleton Jacobian of deformation 𝐽 s related to incompressibility of the solid phase. Uniaxial strain, unidirectional flow examples at higher strain rates (𝒪⁢(1⁢0 2 −1⁢0 3 )⁢s −1 ) and larger deformations (up to 0.2 (or 20%) nominal axial strain) demonstrate the improved physical representation—and numerical stability—of these constitutive model improvements.

42 ENGINEERING↗

Development of an Interactive Valuation Model for Distributed Energy Resource Value: Cooperative Research and Development (Final Report)

The Oklahoma Office of the Secretary of Energy and Environment (OSEE) seeks to address perceived challenges and barriers to the deployment of distributed energy generation resources (DERs), including solar and other technologies, in the state of Oklahoma. The Scope of Work (SOW) associated with CRD-18-762 covers two tasks: (1) Interactive Valuation Model Requirements Document and (2) Consumer Protection Assessment. Follow-on work was completed under a separate CRADA (CRD-19-00788).

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Lifetime measurements in 102 Mo interpreted in the interacting boson model and the X(5) symmetry

Lifetimes of low-lying excited states in 102 Mo populated in the two-neutron transfer reaction 100 Mo ⁢( 18 O, 16 O)⁢ 102 Mo were measured using the recoil distance Doppler shift method at the IFIN-HH Tandem accelerator. Lifetimes of the $2^+_1$, $0^+_2$, $4^+_1$, $2^+_2$, $2^+_3$, $3^+_1$, $6^+_1$, ($0^+_3$), $4^+_2$, ($3^−_1$), and ($5^−_1$) states were obtained. The deduced electromagnetic transition strengths have been compared to calculations performed in the interacting boson model framework including models representing the U(5) and X(5) symmetries. It is found that 102 Mo lies between the U(5) limit and the X(5) critical point symmetry.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Majorana parameters of the interacting boson model of nuclear structure and their implication for 0 ν β β decay

Here, the well-known spherical-deformed-transitional nucleus and potential 0$\textit{νββ}$ emitter 150 Nd and its daughter 150 Sm were investigated in nuclear resonance fluorescence experiments using quasimonoenergetic, linearly polarized γ-ray beams. For both nuclei transitions from the 1 + scissors mode to the $0^+_2$ and $2^+_2$ states were observed for the first time and their respective $\textit{M}$ 1 transition strengths were determined. Through a systematic investigation, a sensitivity of these transition strengths to the three Majorana parameters of the interacting boson model-2 (IBM-2) was established. In combination with the novel experimental data, this poses strong constraints to the Majorana parameters in improved IBM-2 representations of both nuclei. A subsequent recalculation of the nuclear matrix elements (NMEs) for the 150 Nd → 150 Sm 0$\textit{νββ}$ decay in the IBM-2 with these improved representations results in $M^{(0νββ)}_{\text{IBM-2}}[0^+_1]$ = 3.35 for the NME for 0$\textit{νββ}$ decay into the ground state of 150 Sm and $M^{(0νββ)}_{\text{IBM-2}}[0^+_2]$ = 1.30 for 0$\textit{νββ}$ decay to its $0^+_2$ state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Developing a simulator-based satellite dataset for using machine learning techniques to derive aerosol-cloud-precipitation interactions in models and observations in a consistent framework

Aerosol-cloud-precipitation interactions (ACPI) remain a major uncertainty in understanding the Earth’s radiation budget and water cycle (including extremes). After decades of active research, various observationally based metrics have been developed to constrain ACPI in Earth System Models (ESMs), but direct comparison of model and data estimates can confound scientific understanding because limitations and uncertainties in sampling and retrieval procedures may combine with model deficiencies in process representations of ACPI to obstruct understanding. Furthermore, conventional ACPI metrics often vary from one regime to another, and the ACPI process representation in ESMs is also typically derived based on only a limited area/regime (even though the parameterization applies globally). To bridge the gap between models and data and to correctly describe ACPI in all regimes, we propose to construct a new CALIPSO-CloudSat merged dataset that is produced by the same algorithms used in satellite simulators in ESMs, and to use machine learning techniques to derive new ACPI metrics that can be accurately estimated by satellites and can provide meaningful constraints on cloud microphysical process representations in ESMs. The dataset will include measured and retrieved variables for aerosol, cloud, and precipitation from CALIPSO and CloudSat, and environmental variables from meteorological reanalysis. The data will be used to train a neural network to construct the ACPI metrics as a function of environmental conditions. The new ACPI formula will be used to constrain the ACPI in the Energy Exascale Earth System Model (E3SM), and to augment/reformulate the ACPI process representation in the E3SM to improve the simulation of the evolution of the atmosphere under different environmental conditions.

54 ENVIRONMENTAL SCIENCES↗

Modeling Thermal Interactions between Buildings in an Urban Context

Thermal interactions through longwave radiation exchange between buildings, especially in a dense urban environment, can strongly influence a building’s energy use and environmental impact. However, these interactions are either neglected or oversimplified in urban building energy modeling. We developed a new feature in EnergyPlus to explicitly consider this term in the surface heat balance calculations and developed an algorithm to batch calculating the surrounding surfaces’ view factors using a ray-tracing technique. We conducted a case study with a district in the Chicago downtown area to evaluate the longwave radiant heat exchange effects between urban buildings. Results show that the impact of the longwave radiant effects on annual energy use ranges from 0.1% to 3.3% increase for cooling and 0.3% to 3.6% decrease for heating, varying among individual buildings. At the district level, the total energy demand increases by 1.39% for cooling and decreases 0.45% for heating. We also observe the longwave radiation can increase the exterior surface temperature by up to 10 °C for certain exterior surfaces. These findings justify a detailed and accurate way to consider the thermal interactions between buildings in an urban context to inform urban planning and design.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A machine learning-based interaction force model for non-spherical and irregular particles in low Reynolds number incompressible flows

In this study, interaction force of non-spherical particles in low Reynolds number gas-solid flow is investigated by neural network approaches. An artificial neural network (ANN) model is developed to correlate the non-spherical particle shape and the flow conditions with the interaction force. To define the particle shape, spherical harmonic expansion is applied. Furthermore, variational autoencoder model is then used to extract latent geometric features. The latent vector is utilized as an input with the Reynolds number for the ANN. The interaction force data, which is used as output data of the ANN, is obtained by particle resolved direct numerical simulation for 5200 non-spherical particles. The proposed model enables unsupervised extraction for non-spherical particle shapes and accurate predictions on the interaction force without heavy computation. This study provides the model that can explain complicated shapes of particles and be applied to a large scale, computational fluid dynamics simulation.

01 COAL, LIGNITE, AND PEAT↗

SD-pair shell model: Vibrational and rotational limits in the interacting boson–fermion model for like-nucleon system

Typical features of the Bose–Fermi symmetries associated with [Formula: see text] and [Formula: see text] limits in the interacting boson model are described in the SD-pair shell model (SDPSM) framework for like-nucleon system. It is found that the limiting spectra associated with [Formula: see text] (vibrational) and [Formula: see text] (rotational) in the interacting boson–fermion model can be well reproduced in the SDPSM framework. It is shown that coupling of collective degree of freedom with the single-particle degree of freedom in the pairing interaction almost has no effect on the spectrum of odd-A system comparing to the adjacent even–even system, while the coupling of the collective degree of freedom with the single-particle degree of freedom in the quadrupole–quadrupole interaction results in the level splitting of several excited levels.

Physics↗

A new method for diagnosing effective radiative forcing from aerosol–cloud interactions in climate models

Aerosol–cloud interactions (ACIs) are a leading source of uncertainty in estimates of the historical effective radiative forcing (ERF). One reason for this uncertainty is the difficulty in estimating the ERF from aerosol–cloud interactions (ERFaci) in climate models, which typically requires multiple calls to the radiation code. Most commonly used methods also cannot disentangle the contributions from different processes to ERFaci. Here, we develop a new, computationally efficient method for estimating the shortwave (SW) ERFaci from liquid clouds using histograms of monthly averaged cloud fraction partitioned by cloud droplet effective radius (r e ) and liquid water path (LWP). Multiplying the histograms with SW cloud radiative kernels gives the total SW ERFaci from liquid clouds, which can be decomposed into contributions from the Twomey effect, LWP adjustments, and cloud fraction (CF) adjustments. We test the method with data from five CMIP6-era models, using the Moderate Resolution Imaging Spectroradiometer (MODIS) satellite instrument simulator to generate the histograms. Our method gives similar total SW ERFaci estimates to other established methods in regions of prevalent liquid cloud and indicates that the Twomey effect, LWP adjustments, and CF adjustments have contributed -0.34 ± 0.23, -0.22 ± 0.13, and -0.09 ± 0.11 W m -2 , respectively, to the effective radiative forcing of the climate since 1850 in the ensemble mean (95 % confidence). These results demonstrate that widespread adoption of a MODIS r e –LWP joint histogram diagnostic would allow the SW ERFaci and its components to be quickly and accurately diagnosed from climate model outputs, a crucial step for reducing uncertainty in the historical ERF.

54 ENVIRONMENTAL SCIENCES↗

Modeling Market Interactions

This slide deck provides a high-level overview of wholesale electricity market design model development activities occurring at NREL. The motivation for this work, as well as the broader NREL modeling landscape in which it belongs, is first presented. Then the various model formulations and features are described.

agent-based simulation↗

Topological phase diagram of generalized Su-Schrieffer-Heeger models with interactions

We investigate interacting Su-Schrieffer-Heeger (SSH) chains with two- and three-site unit cells using density matrix renormalization group (DMRG) simulations. By selecting appropriate filling fractions and sweeping across interaction strength 𝐽 𝑧 and dimerization 𝛿, we map out their phase diagrams and identify transition lines via entanglement entropy and magnetization measurements. In the two-site model, we observe the emergence of an interaction-induced antiferromagnetic intermediate phase between the topologically trivial and nontrivial regimes, as well as a critical region at negative 𝐽 𝑧 with suppressed magnetization and finite-size scaling of entanglement entropy. In contrast, the three-site model lacks an intermediate phase and exhibits asymmetric edge localization and antiferromagnetic ordering in both positive and negative 𝐽 𝑧 regimes. We further examine the response of edge states to Ising perturbations. In the two-site model, zero-energy edge modes are topologically protected and remain robust up to a finite interaction strength. However, in the three-site model, where the edge states reside at finite energy, this protection breaks down. Despite this, the edge-localized nature of these states survives in the form of polarized modes whose spatial profiles reflect the noninteracting limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Temperature dependence of magnetic anisotropy and magnetoelasticity from classical spin-lattice calculations

Here we present a classical molecular-spin dynamics (MSD) methodology that enables accurate computations of the temperature dependence of the magnetocrystalline anisotropy as well as magnetoelastic properties of magnetic materials. The nonmagnetic interactions are accounted for by a spectral neighbor analysis potential (SNAP) machine-learned interatomic potential, whereas the magnetoelastic contributions are accounted for using a combination of an extended Heisenberg Hamiltonian and a Néel pair interaction model, representing both the exchange interaction and spin-orbit-coupling effects, respectively. All magnetoelastic potential components are parameterized using a combination of first-principles and experimental data. Our framework is applied to the α phase of iron. Initial testing of our MSD model is done using a 0 K parametrization of the Néel interaction model. After this, we examine how individual Néel parameters impact the $B$ 1 and $B$ 2 magnetostrictive coefficients using a moment-independent δ sensitivity analysis. The results from this study are then used to initialize a genetic algorithm optimization which explores the Néel parameter phase space and tries to minimize the error in the B 1 and B 2 magnetostrictive coefficients in the range of 0–1200 K. Our results show that while both the 0 K and genetic algorithm optimized parametrization provide good experimental agreement for $B$ 1 and $B$ 2 , only the genetic algorithm optimized results can capture the second peak in the $B$ 1 magnetostrictive coefficient which occurs near approximately 800 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Vertex finding in neutrino-nucleus interaction: a model architecture comparison

We compare different neural network architectures for machine learning algorithms designed to identify the neutrino interaction vertex position in the MINERvA detector. The architectures developed and optimized by hand are compared with the architectures developed in an automated way using the package “Multi-node Evolutionary Neural Networks for Deep Learning” (MENNDL), developed at Oak Ridge National Laboratory. While the domain-expert hand-tuned network was the best performer, the differences were negligible and the auto-generated networks performed as well. There is always a trade-off between human, and computer resources for network optimization and this work suggests that automated optimization, assuming resources are available, provides a compelling way to save significant expert time.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗

Tensor force role in β decays analyzed within the Gogny-interaction shell model

The half-life of the famous C 14 β decay is anomalously long, with different mechanisms: the tensor force, cross-shell mixing, and three-body forces, proposed to explain the cancellations that lead to a small transition matrix element. In this study, we revisit and analyze the role of the tensor force for the β decay of C 14 as well as of neighboring isotopes. We add a tensor force to the Gogny interaction, and derive an effective Hamiltonian for shell-model calculations. The calculations were carried out in a p – s d model space to investigate cross-shell effects. Furthermore, we decompose the wave functions according to the total orbital angular momentum L in order to analyze the effects of the tensor force and cross-shell mixing. The inclusion of the tensor force significantly improves the shell-model calculations of the β -decay properties of carbon isotopes. In particular, the anomalously slow β decay of C 14 can be explained by the isospin T = 0 part of the tensor force, which changes the components of N 14 with the orbital angular momentum L = 0 , 1 , and results in a dramatic suppression of the Gamow-Teller transition strength. At the same time, the description of other nearby β decays are improved. Decomposition of wave function into L components illuminates how the tensor force modifies nuclear wave functions, in particular suppression of β -decay matrix elements. Cross-shell mixing also has a visible impact on the β -decay strength. Inclusion of the tensor force does not seem to significantly change, however, binding energies of the nuclei within the phenomenological interaction.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Measurement of neutron production in atmospheric neutrino interactions at Super-Kamiokande

We present measurements of total neutron production from atmospheric neutrino interactions in water, analyzed as a function of electron-equivalent visible energy over a range of 30 MeV to 10 GeV. These results are based on 4,270 days of data collected by Super-Kamiokande, including 564 days with 0.011 wt% gadolinium added to enhance neutron detection. Neutron signal selection is based on a neural network trained on simulation, with its performance validated using an Am/Be neutron point source. The measurements are compared to predictions from neutrino event generators combined with various hadron-nucleus interaction models, which include an intranuclear cascade model and a nuclear deexcitation model. We observe significant variations in the predictions depending on the choice of hadron-nucleus interaction model. We discuss key factors that contribute to describing our data, such as in-medium effects in the intranuclear cascade and the accuracy of statistical evaporation modeling.

Cherenkov detectors↗