Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation analytics”

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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Gravitational waves from global cosmic strings and cosmic archaeology

Global cosmic strings are predicted in many motivated extensions to the Standard Model of particle physics, with close connections to axion dark matter physics. Recent studies suggest that, although subdominant relative to Goldstone emission, gravitational wave (GW) signals from global strings can be detectable with current and planned GW detectors such as LIGO, LISA, DECIGO/BBO, ET/CE and AEDGE/AION, as well as pulsar timing arrays such as PPTA, NANOGrav and SKA. This work is an extensive, updated study on GWs from a global cosmic string network, taking into account of the most recent developments related to the subject. The main analysis is based on the analytical Velocity-dependent One-Scale (VOS) model calibrated with recent simulation results, which provides a generic protocol for such calculations with details given. We also demonstrate how the GW signal can be influenced with variations to the baseline model: this includes considering the uncertainties of model parameters and the potential deviation from the conventional VOS model prediction (i.e. the scaling behavior) as suggested by some of the recent simulation results. Furthermore, we investigated in detail the effect of a non-standard cosmology (e.g. early matter domination or kination) or new particle species on the GW signals from global strings. We demonstrate that the frequency spectrum of GW background from global cosmic strings can be used to probe the cosmic history prior to the Big Bang nucleosynthesis (BBN) (i.e. the primordial dark age) up to a temperature of T ~ 10 8 GeV.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

An integrated numerical model for coupled poro-hydro-mechanics and fracture propagation using embedded meshes

Integrated models for fluid-driven fracture propagation and general multiphase flow in porous media are valuable to the study and engineering of several systems, including hydraulic fracturing, underground disposal of waste, and geohazard mitigation across such applications. This work extends the coupled model multiphase flow and poromechanical model of Ren et al. (2018) to admit fracture propagation (FP). The coupled XFEM-EDFM scheme utilizes a separate fracture mesh that is embedded on a static background mesh. The onset and dynamics of fracture propagation are governed by the equivalent stress intensity factor (SIF) criterion. A domain-integral method (J integral) is applied to compute this information. Additionally, an adaptive time-marching scheme is proposed to rapidly restrict and grow temporal resolution to match the underlying time-scales. The proposed model is verified with analytical solutions, and shows the capability to accurately and adaptively co-simulate fluid transport and deformation as well as the propagation of multiple fractures.

42 ENGINEERING↗
DOE OSTI2022

Frhodo: A program for simulating chemical kinetic measurements and optimizing kinetic mechanisms

Chemical kinetic simulations are frequently used to extract kinetic and mechanistic information from experimental data by fitting simulations to the data. For simple systems of consecutive single channel reactions this is a relatively straightforward task. However, as the complexity of the system increases and particularly for systems involving competitive multi-channel reactions manually optimizing a mechanism to obtain the best fits to a range of experimental data quickly becomes a time-intensive and challenging process, especially when rate coefficients have pressure and temperature dependencies. There is considerable potential for exploiting automated optimization methods to rapidly screen reaction mechanisms and optimize the fit by adjusting rate coefficients within well-defined constraints. A new chemical kinetics simulation program, Frhodo, has been developed for this purpose. The program allows either manual simulation of individual experiments or fully automated optimization against a range of experimental data by adjusting user selected reactions within user defined constraints. Frhodo incorporates machine learning-based optimization routines that work by either minimizing standardized residuals or through Bayesian parameter estimation. Frhodo's optimization capabilities are demonstrated by reexamining two previously studied systems of relevance to combustion. Dissociation of diacetyl followed by recombination of methyl radicals is an example of a sequence of single channel, consecutive reactions. Pyrolysis of 2-methyl furan exemplifies multi-channel, unimolecular reactions that are often encountered and have temperature and pressure dependencies in the rate coefficients of the channels and branching ratios between them. For both systems the optimization routines resulted in solutions similar to the original studies. Finally, comments are made on the strengths and limitations of the approach.

0D simulation↗
DOE OSTI2023

Quantum Simulation of Polarized Light-Induced Electron Transfer with a Trapped-Ion Qutrit System

Electron transfer within and between molecules is crucial in chemistry, biochemistry, and energy science. Here this study describes a quantum simulation method that explores the influence of light polarization on electron transfer between two molecules. By implementing precise and coherent control among the quantum states of trapped atomic ions, we can induce quantum dynamics that mimic the electron-transfer dynamics in molecules. We use three-level systems (qutrits), rather than traditional two-level systems (qubits), to enhance the simulation efficiency and realize high-fidelity simulations of electron-transfer dynamics. We treat the quantum interference between the electron coupling pathways from a donor with two degenerate excited states to an acceptor and analyze the transfer efficiency. We also examine the potential error sources that enter the quantum simulations. The trapped-ion systems have favorable scalings with system size compared to those of classical computers, promising access to richer electron-transfer simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Hydride Accessibility and Reactivity in the Configurational and Stoichiometric Space of β-Ga 2 O 3 for CO 2 Hydrogenation

Understanding how surface species evolve under reaction conditions is essential for improving catalyst design for efficient CO 2 hydrogenation. This work combines systematic DFT calculations with grand canonical sampling to investigate the stability and reactivity of Ga–H species on β-Ga 2 O 3 across a range of reaction conditions. Initial DFT studies reveal that when Ga–H species are present, they facilitate formate formation via a low-barrier pathway, largely independent of the surface termination or hydrogen site. However, grand canonical sampling shows that under a broad range of reaction conditions─especially at high oxygen chemical potentials associated with high water content─Ga–H species are thermodynamically inaccessible. Furthermore, adsorbed water molecules can block reactive sites, inhibiting CO 2 activation even when hydrides are present. These findings suggest that the lack of accessible hydride species, rather than their intrinsic reactivity, could contribute to reduced catalytic performance of β-Ga 2 O 3 under more oxidizing, high-conversion conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Plasma heating and current sheet structure in anti-parallel magnetic reconnection

A theoretical model and an analytic theory of current sheet structure are presented for understanding anti-parallel driven magnetic reconnection in 2-1/2 dimension in collisionless plasmas. The theoretical model provides formulation to compute the current sheet y-profiles by specifying the profiles of electron and ion flow velocities V ex (x,y) and V ix (x,y). The current sheet solutions depend on the plasma density n in , merging magnetic field B 0 , ion velocity v i , and electron velocity ve in the upstream and the S evz =V ez /V dz parameter where V ez is the electron velocity accelerated by the reconnection electric field E z in the electron orbit meandering region, V dz ≃cE y /B x is the $\overrightarrow{E}$ × $\overrightarrow{B}$ drift velocity as electrons enter the orbit meandering region, B x is the merging magnetic field, and E y is the electrostatic electric field. With simplifying assumptions on the y-profiles of V ex and V ix , we have also developed an analytic theory of the current sheet structure. Analytic expressions for the anomalous resistivity, the electrostatic potential drop, and the maximum E y amplitude E max are obtained. The analytic results agree reasonably well with both the particle-in-cell simulation results and the numerical solutions of the theoretical model. The ions energy gain due to the potential drop is $∝B^2_0/n_{in}$. The electron energy gain is $∝(B^2_0 /8πn_{in})S_{evz}$. The $B^2_0/n_{in}$ scaling of the average ion and electron energy gains are consistent with laboratory experiments and space plasma observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

SMBH seeds from dissipative dark matter

The existence of supermassive black holes (SMBHs) with masses greater than ~ 10 9 M ⊙ at high redshift ( z ≳ 7) is difficult to accommodate in standard astrophysical scenarios. We study the possibility that (nearly) totally dissipative self-interacting dark matter (tdSIDM)–in rare, high density dark matter fluctuations in the early Universe — produces SMBH seeds through catastrophic collapse. We use a semi-analytic model, tested and calibrated by a series of N-body simulations of isolated dark matter halos, to compute the collapse criteria and timescale of tdSIDM halos, where dark matter loses nearly all of its kinetic energy in a single collision in the center-of-momentum frame. Applying this model to halo merger trees, we empirically assign SMBH seeds to halos and trace the formation and evolution history of SMBHs. We make predictions for the quasar luminosity function, the M BH -σ v * relation, and cosmic SMBH mass density at high redshift and compare them to observations. We find that a dissipative dark matter interaction cross-section of σ/ m ~ 0.05 cm 2 /g is sufficient to produce the SMBHs observed in the early Universe while remaining consistent with ordinary SMBHs in the late Universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Widen the resonance at ultra-high energies: novel probes of neutrino self-interactions in the high-mass regime

Neutrino self-interactions beyond the Standard Model are well motivated by the nonzero masses of neutrinos, which are the only known particles guaranteed to have new physics. Meanwhile, cosmic messengers, especially neutrinos, play a central role in probing new physics, as they provide experimental conditions far beyond the reach of laboratories and serve as the link between laboratory fundamental-physics discoveries and their roles in the Universe, where many new physics motivations originate.In this work, we propose a novel probe of neutrino self-interactions through ultrahigh-energy neutrinos scattering off the cosmic neutrino background when the lightest neutrino species remains relativistic today.This allows us to “Widen the Resonance” of such scattering [1].In addition, we also provide a semi-analytic framework for cosmogenic ultrahigh-energy neutrino production, avoiding computationally intensive simulations and yielding results precise enough for beyond-the-standard-model studies.The widened resonance enables future ultrahigh-energy neutrino telescopes,in particular GRAND, to probe mediator masses from MeV to GeV, reaching couplings down to g ∼ 10$^{-3}$ — up to two orders of magnitude beyond current bounds.Our results enhance the discovery potential of neutrino self-interactions in the high-mass regime, potentially offering crucial insights into the connections between the neutrino sector and dark sector.

Machado, Pedro A.N. [Fermilab] (ORCID:000000029118↗
DOE OSTI2025

Energetic particle physics: Chapter 7 of the special issue: on the path to tokamak burning plasma operation

We review the physics of energetic particles (EPs) in magnetically confined burning fusion plasmas with focus on advances since the last update of the ITER Physics Basis (Fasoli et al 2007 Nucl. Fusion 47 S264). Topics include basic EP physics, EP generation, diagnostics of EPs and instabilities, the interaction of EPs and thermal plasma instabilities, EP-driven instabilities, energetic particle modes (EPMs), and turbulence, linear and nonlinear stability and simulation of EP-driven instabilities and EPMs, 3D effects, scenario optimization strategies based on EP phase-space control, EPs in reduced field scenarios in ITER before DT, and the physics of runaway electrons. We describe the simulation and modeling of EPs in fusion plasmas, including instability drive and damping as well as EP transport, with a range of approaches from first-principles to reduced models, including gyrokinetic simulations, kinetic-MHD models, gyrofluid models, reduced models, and semi-analytical approaches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Heavy flavor jet substructure at lepton colliders

We provide a detailed analysis of event-shape observables, namely the energy correlation function and jet angularity, for heavy-flavor jets produced in electron-positron collisions, focusing on quantum chromodynamics (QCD) interactions. Using modern jet substructure techniques, we investigate the dead-cone effect, where QCD radiation is suppressed around a heavy quark within an angle proportional to its mass. Our analysis achieves next-to-leading logarithmic accuracy, combined with partial fixed-order contributions, to improve the description of the transition near the dead-cone threshold. To ensure a comprehensive perspective, we compare our analytical results with predictions from the pythia, herwig, and sherpa Monte Carlo simulations at past and future lepton colliders.

Perturbative QCD↗
DOE OSTI2021

Magnetic Field Amplification by a Nonlinear Electron Streaming Instability

Magnetic field amplification by relativistic streaming plasma instabilities is central to a wide variety of high-energy astrophysical environments as well as to laboratory scenarios associated with intense lasers and electron beams. We report on a new secondary nonlinear instability which arises for relativistic dilute electron beams after the saturation of the linear Weibel instability. This instability grows due to the transverse magnetic pressure associated with the beam current filaments, which cannot be quickly neutralized due to the inertia of background ions. We show that it can amplify the magnetic field strength and spatial scale by orders of magnitude, leading to large-scale plasma cavities with strong magnetic field and to very efficient conversion of the beam kinetic energy into magnetic energy. The instability growth rate, saturation level, and scale length are derived analytically and shown to be in good agreement with fully-kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

A Data-Driven Approach to Nation-Scale Building Energy Modeling

In 2019, 125 million U.S. residential and commercial buildings consumed $412 billion in energy bills. These buildings currently consume 40% of the nation's primary energy, 73% of electricity, 80% of energy during peak electric grid use, and responsible for 39% of greenhouse gas emissions [14]. Urban-scale building energy modeling has grown significantly in the past decade, allowing individual campuses or communities of buildings to be modeled, simulated, and cost-effective solutions for intelligent management to be identified and implemented. While traditionally limited to individual counties and usually less than 2,000 buildings, the Automatic Building Energy Modeling (AutoBEM) soft-ware suite has been developed to process unconventional, nation-scale data sources to generate unique OpenStudio and EnergyPlus models of each building. Through the use of High Performance Computing (HPC) resources, every U.S. building has been simulated. This paper showcases the data layout, node partitioning, algorithmic approaches, and analytic results that were used to create, share, and analyze 124.4 million U.S. building models.

Berres, Andy↗
DOE OSTI2020

Progress on Understanding Rayleigh–Taylor Flow and Mixing Using Synergy Between Simulation, Modeling, and Experiment

Simultaneous advances in numerical methods and computing, theoretical techniques, and experimental diagnostics have all led independently to better understanding of Rayleigh–Taylor (RT) instability, turbulence, and mixing. In particular, experiments have provided significant motivation for many simulation and modeling studies, as well as validation data. Numerical simulations have also provided data that is not currently measurable or very difficult to measure accurately in RT unstable flows. Thus, simulations have also motivated new measurements in this class of buoyancy-driven flows. This overview discusses simulation and modeling studies synergistic with experiments and examples of how experiments have motivated simulations and models of RT instability, flow, and mixing. First, a brief summary of measured experimental and calculated simulation quantities, of experimental approaches, and of issues and challenges in the simulation and modeling of RT experiments is presented. Implicit large-eddy, direct numerical, and large-eddy simulations validated using RT experimental data are then discussed. This is followed by a discussion of modeling using analytical, modal, buoyancy–drag, and turbulent transport models of RT mixing experiments. The discussion will focus on three-dimensional RT mixing arising from multimode perturbations. Finally, this focused review concludes with a perspective on future simulation, modeling, and experimental directions for further research. Research in simulation and modeling of RT unstable flows, coupled with experiments, has made significant progress over the past several decades. This overview serves as an opportunity to both discuss progress and to stimulate future research on simulation and modeling of this unique class of hydrodynamically unstable turbulent flows.

42 ENGINEERING↗
DOE OSTI2022

maestro

MÆSTRO stands for Multi-fidelity Adaptive Ensemble Stochastic Trust Region Optimization and it is a plug n play derivate fee stochastic optimization solver. The problem being considered in MÆSTRO involves fitting Monte Carlo simulations that describe complex phenomena to experiments. This is done by finding parameters of the resource intensive and noisy simulation that yield the least squares objective function value to the noisy experimental data. This problem is solved using a stochastic trust-region optimization algorithm where in each iteration, a local approximation of the simulation signal and of the simulation noise is constructed over data, which is obtained by running the simulation at strategically placed design points within the trust-region around the current iterate. Then the simulation components of the objective are replaced by their approximations and this analytical and closed-form optimization problem is solved to find the next iterate within the trust-region. Then the trust region is moved and the iterations continue until a satisfactory convergence criteria is met.

KRISHNAMOORTHY, MOHAN↗
DOE OSTI2020

Benchmark for Fuel Shuffling and Depletion for Pebble-Bed Reactors

Pebble bed reactors have specific operational characteristics when their fuel-cycle and fueling operations are considered. They are specifically distinguished by other type of nuclear reactor designs by their online fuel recycling scheme, where the fuel elements that have not yet reached discharge burnup can be reloaded and recycled continuously during normal operation. The fuel in a pebble bed reactor is not stationary and stochastically moves through the core once or several times during its lifetime, which allows them to operate without requiring a large excess reactivity hold for the burnup. However, this characteristic of pebble bed reactors introduces challenges in simulation, as each pebble can take many different trajectories through the core, its composition depends on the details of the irradiation history that is unique to its aggregated path through the core. For predicting the safety performance characteristics, such as source term, maximum fuel temperatures and fuel failure rates, etc., it is important to accurately incorporate the movement of pebbles through the core during their lifetime in a multi-physics simulation together with other phenomena. The equilibrium core analysis for pebble bed reactors are performed with multi-physics tools including fuel depletion in a multi pass reload coupled to the fuel movement. Currently, there are only a few legacy multi-physics simulation tools that can implement the pebble flow characteristics and perform equilibrium core analysis for pebble bed reactors. However, there are development efforts on-going under Department of Energy's Nuclear Energy Advanced Modelling and Simulation program and also in private industry for including these capabilities into their modelling and simulation tools. Any new development in the modelling and simulation tools needs to be validated by using tools such as experiments, analytical solutions or code-to-code benchmarks. In this work, a code-to-code benchmark for the equilibrium core analysis capability of pebble bed reactors was developed. Multiple cases were identified to capture different fuel cycle strategies that can be used in PBRs. The results of each case are presented in terms of overall equilibrium core characteristics: the discharge burnup; spatial burnup distribution; spatial isotopic distributions; axial and radial neutron flux distributions and power history of fuel elements per pass through the core for both a prototypical pebble bed High Temperature Gas-cooled Reactor and a prototypical pebble bed Fluoride-salt cooled High temperature Reactor.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Compare low-gain 1D and 3D models toward a storage ring-based X-ray free-electron-laser oscillator

Fourth-generation diffraction-limited synchrotron light sources (DL-SLS) present a compelling platform for advanced light source applications due to their ultra-low emittance (around tens of picometers), sub-micrometer beam stability, and high repetition rates (greater than 1 MHz). Despite these advantages, fundamental limitations, such as long bunch durations (tens of picoseconds), large relative energy spread (~10⁻³), and low peak current (<300 A), have hindered their effectiveness in free-electron laser (FEL) applications. Using a low-gain three-dimensional (3D) model based on NSLS-II upgrade lattice parameters, we demonstrate that a 6% single-pass gain is achievable with an optimized transverse gradient undulator (TGU). Concurrently, recent advances in X-ray optics—specifically, a four-mirror bow-tie resonator utilizing diamond mirrors and beryllium compound refractive lenses (CRLs)—enable total round-trip losses below 4% and an output coupling efficiency near 1%, rendering a Storage Ring based X-ray FEL oscillator (XFELO) a realistic possibility. In the 3D framework, the TGU effectively mitigates the detrimental effects of large energy spread by spatially dispersing the beam, leading to enhanced gain relative to conventional one-dimensional (1D) models. We establish both analytical and numerical connections between the 1D and 3D descriptions, leveraging GENESIS simulations to map equivalent gain conditions. While 1D resonance is defined in the energy domain, the 3D TGU system shifts the resonance dependence to transverse gradient and dispersion parameters. Once optimized, this framework allows predictive modeling using 1D formulas, which can then be validated through full 3D simulations.

43 PARTICLE ACCELERATORS↗