Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “system simulation”

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 1,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Energy response characterization of current-mode MeV Cherenkov aerogel detectors

Current-mode Cherenkov aerogel detectors have been used to measure time resolved spectra of MeV radiographic x-ray sources and have been considered to replace gas Cherenkov systems. The simulations of the Cherenkov response of aerogel detectors have uncertainties in how to model the complex structure of aerogels. Data from three sites—the Idaho electron beam, the Los Alamos Site electron beam and microtron sources at Los Alamos—have been used to calibrate the relative energy response of 15 aerogels, varying in density from 10 to 430 mg/cm 3 , from incident electrons with energies from 1 to 12 MeV. Detector simulations need to introduce density variation in the aerogels to match the collected data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Atlantic hurricane surge response to geoengineering

Devastating Atlantic hurricanes are relatively rare events. However their intensity and frequency in a warming world may rapidly increase by a factor of 2-7 for each degree of increase in mean global temperature. Geoengineering by stratospheric sulphate aerosol injection cools the tropics relative to the polar regions, including the hurricane main development region in the Atlantic, suggesting that geoengineering may be an effective method of controlling hurricanes. We examine this hypothesis using 8 Earth System Model simulations of climate under the GeoMIP G3 and G4 schemes that use stratospheric aerosols to reduce the radiative forcing under the RCP4.5 scenario. Global mean temperature increases are greatly ameliorated by geoengineering, and tropical temperature increases are at most half of those in RCP4.5, but sulphate injection would have to double between 2020 and 2070 to balance RCP 4.5 to nearly 10 Tg SO2 yr-1, with consequent implications for damage to stratospheric ozone. We project changes in storm frequencies using a temperature-dependent Generalized Extreme Value statistical model calibrated by historical storm surges from 1923 and observed temperatures. The numbers of storm surge events as big as the one that caused the 2005 Katrina hurricane are reduced by about 50% compared with no geoengineering, but this is only marginally statistically significant. Furthermore, when sea level rise differences at 2070 between RCP4.5 and geoengineering are factored in to coastal flood risk, we find that expected flood levels are reduced by about 40 cm for 5 year events and perhaps halved for 50 year surges.

58 GEOSCIENCES↗
DOE OSTI2021

Entropic formation of a thermodynamically stable colloidal quasicrystal with negligible phason strain

Quasicrystals have been discovered in a variety of materials ranging from metals to polymers. Yet, why and how they form is incompletely understood. In situ transmission electron microscopy of alloy quasicrystal formation in metals suggests an error-and-repair mechanism, whereby quasiperiodic crystals grow imperfectly with phason strain present, and only perfect themselves later into a high-quality quasicrystal with negligible phason strain. The growth mechanism has not been investigated for other types of quasicrystals, such as dendrimeric, polymeric, or colloidal quasicrystals. Soft-matter quasicrystals typically result from entropic, rather than energetic, interactions, and are not usually grown (either in laboratories or in silico) into large-volume quasicrystals. Consequently, it is unknown whether soft-matter quasicrystals form with the high degree of structural quality found in metal alloy quasicrystals. Here, we investigate the entropically driven growth of colloidal dodecagonal quasicrystals (DQCs) via computer simulation of systems of hard tetrahedra, which are simple models for anisotropic colloidal particles that form a quasicrystal. Using a pattern recognition algorithm applied to particle trajectories during DQC growth, we analyze phason strain to follow the evolution of quasiperiodic order. As in alloys, we observe high structural quality; DQCs with low phason strain crystallize directly from the melt and only require minimal further reduction of phason strain. Finally, we also observe transformation from a denser approximant to the DQC via continuous phason strain relaxation. Our results demonstrate that soft-matter quasicrystals dominated by entropy can be thermodynamically stable and grown with high structural quality––just like their alloy quasicrystal counterparts.

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

Numerical Simulation and Experimental Comparison of System Analysis Module 1D Mixing Model for Cold Shock Transients in the Gallium Thermal-Hydraulic Mixing Facility

Abstract Liquid metals are being investigated as coolants in many advanced reactor designs because of their high thermal conductivity and effectiveness at high temperatures. However, they often pose challenges to reactor operation and safety because of the complex thermal mixing and stratification in the plenum of pool-type reactor designs. The advanced system analysis code System Analysis Module (SAM) currently under development at Argonne National Laboratory aims to develop and implement thermal mixing models to accurately capture these complex thermal fluid behaviors. In this study, the SAM thermal mixing model was compared against experimental data from the Gallium Thermal-Hydraulic Experiment facility, a scaled liquid metal test facility that uses gallium as a surrogate fluid to investigate the stratification and thermal mixing of low-Prandtl-number fluids in the upper plenum of a liquid metal-cooled reactor. Two cold shock transient cases were used: one with stable stratified flow (Ri = 32) and one with stronger thermal mixing (Ri = 0.5). The resultant temperatures were then compared with the experimental temperatures over the entire plenum to assess the ability of the mixing models to capture the thermal behavior and to better correspond mixing parameters to various flow scenarios. Generally, the zero-dimensional mixing model was more capable of capturing the bulk temperature of the component modeled assuming that an accurate mass flow rate was provided, but it was inherently unable to capture thermal gradients in space. The one-dimensional mixing model was capable of capturing that the thermal gradients provided accurate selection of the mixing coefficients. Further, the temperature at the outlet junction was compared over time for each of the mixing models with the recorded experimental temperature. The implemented mixing models demonstrated the ability to effectively capture the overall thermal behavior for stronger mixing scenarios but struggled with more stably stratified flows. It was found that a system analysis code's covering of the entire range of different operating conditions still remains a challenging task, and it is suggested that further model and closure improvements are necessary to accurately capture complex thermal mixing and stratification phenomena.

stratification↗
DOE OSTI2014

Seven new carbon-enhanced metal-poor RR Lyrae stars

We report estimated carbon-abundance ratios, [C/Fe], for seven newly discovered carbon-enhanced metal-poor (CEMP) RR Lyrae stars. These are well-studied RRab stars that had previously been selected as CEMP candidates based on low-resolution spectra. For this pilot study, we observed eight of these CEMP RR Lyrae candidates with the Wide Field Spectrograph on the ANU 2.3 m telescope. Prior to this study, only two CEMP RR Lyrae stars had been discovered: TY Gru and SDSS J1707+58. We compare our abundances to new theoretical models of the evolution of low-mass stars in binary systems. These simulations evolve the secondary stars, post accretion from an asymptotic giant-branch (AGB) donor, all the way to the RR Lyrae stage. The abundances of CEMP RR Lyrae stars can be used as direct probes of the nature of the donor star, such as its mass, and the amount of material accreted onto the secondary. We find that the majority of the sample of CEMP RR Lyrae stars is consistent with AGB donor masses of around 1.5-2.0 M {sub ☉} and accretion masses of a few hundredths of a solar mass. Future high-resolution studies of these newly discovered CEMP RR Lyrae stars will help disentangle the effects of the proposed mixing processes that occur in such objects.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Reduced-order autodifferentiable ensemble Kalman filters

This paper introduces a computational framework to reconstruct and forecast a partially observed state that evolves according to an unknown or expensive-to-simulate dynamical system. Our reduced-order autodifferentiable ensemble Kalman filters (ROAD-EnKFs) learn a latent low-dimensional surrogate model for the dynamics and a decoder that maps from the latent space to the state space. The learned dynamics and decoder are then used within an EnKF to reconstruct and forecast the state. Numerical experiments show that if the state dynamics exhibit a hidden low-dimensional structure, ROAD-EnKFs achieve higher accuracy at lower computational cost compared to existing methods. If such structure is not expressed in the latent state dynamics, ROAD-EnKFs achieve similar accuracy at lower cost, making them a promising approach for surrogate state reconstruction and forecasting.

Mathematics↗
DOE OSTI2025

Spin–orbit torque in a three-fold-symmetric bilayer and its effect on magnetization dynamics

Field-free switching of perpendicular magnetization has previously been observed in an epitaxial L1 1 -ordered CoPt/CuPt bilayer and attributed to spin–orbit torque (SOT) arising from the crystallographic 3m point group of the interface. Using a first-principles nonequilibrium Green's function formalism combined with the Anderson disorder model, we calculate the angular dependence of the SOT in a CoPt/CuPt bilayer and find that the magnitude of the 3m SOT is about 20% of the conventional dampinglike SOT. We further study the magnetization dynamics in perpendicularly magnetized films in the presence of 3m SOT and Dzyaloshinskii–Moriya interaction, using the equations of motion for domain wall dynamics and micromagnetic simulations. For systems where strong interfacial DMI results in the Néel character of domain walls, we find that a very large current density is required to achieve deterministic switching, because reorientation of the magnetization inside the domain wall is necessary to induce the switching asymmetry. For thicker films with relatively weak interfacial DMI and the Bloch character of domain walls, deterministic switching is possible at much smaller current densities, which may explain the recent experimental findings.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Simulations of a FIR Oscillator with Large Slippage parameter at Jefferson Lab for FIR/UV pump-probe experiments

We previously proposed a dual FEL configuration on the UV Demo FEL at Jefferson Lab that would allow simultaneous lasing at FIR and UV wavelengths. The FIR source would be an FEL oscillator with a short wiggler providing diffraction-limited pulses with pulse energy exceeding 50 microJoules, using the exhaust beam from a UVFEL as the input electron beam. Since the UV FEL requires very short pulses, the input to the FIR FEL is extremely short compared to a slippage length and the usual Slowly Varying Envelope Approximation (SVEA) does not apply. We use a non-SVEA code to simulate this system both with a small energy spread (UV laser off) and with large energy spread (UV laser on).

43 PARTICLE ACCELERATORS↗
DOE OSTI2014

Fish Passage though Hydropower Turbines: Simulating Blade Strike using the Discrete Element Method

mong the hazardous hydraulic conditions affecting anadromous and resident fish during their passage though turbine flows, two are believed to cause considerable injury and mortality: collision on moving blades and decompression. Several methods are currently available to evaluate these stressors in installed turbines, i.e. using live fish or autonomous sensor devices, and in reduced-scale physical models, i.e. registering collisions from plastic beads. However, a priori estimates with computational modeling approaches applied early in the process of turbine design can facilitate the development of fish-friendly turbines. In the present study, we evaluated the frequency of blade strike and nadir pressure environment by modeling potential fish trajectories with the Discrete Element Method (DEM) applied to fish-like composite particles. In the DEM approach, particles are subjected to realistic hydraulic conditions simulated with computational fluid dynamics (CFD), and particle-structure interactions—representing fish collisions with turbine blades—are explicitly recorded and accounted for in the calculation of particle trajectories. We conducted transient CFD simulations by setting the runner in motion and allowing for better turbulence resolution, a modeling improvement over the conventional practice of simulating the system in steady state which was also done here. While both schemes yielded comparable bulk hydraulic performance, transient conditions exhibited a visual improvement in describing flow variability. We released streamtraces (steady flow solution) and DEM particles (transient solution) at the same location from where sensor fish (SF) have been released in field studies of the modeled turbine unit. The streamtrace-based results showed a better agreement with SF data than the DEM-based nadir pressures did because the former accounted for the turbulent dispersion at the intake but the latter did not. However, the DEM-based strike frequency is more representative of blade-strike probability than the steady solution is, mainly because DEM particles accounted for the full fish length, thus resolving (instead of modeling) the collision event.

CFD↗
DOE OSTI2022

A robust solver for wavefunction-based density functional theory calculations

A new iterative solver is proposed to efficiently calculate the ground state electronic structure in density functional theory calculations. This algorithm is particularly useful for simulating physical systems considered difficult to converge by standard solvers, in particular metallic systems. Here, the effectiveness of the proposed algorithm is demonstrated on various applications.

42 ENGINEERING↗
DOE OSTI2023

Electron-scale current sheets and energy dissipation in 3D kinetic-scale plasma turbulence with low electron beta

ABSTRACT Three-dimensional kinetic-scale turbulence is studied numerically in the regime where electrons are strongly magnetized (the ratio of plasma species pressure to magnetic pressure is βe = 0.1 for electrons and βi = 1 for ions). Such a regime is relevant in the vicinity of the solar corona, the Earth’s magnetosheath, and other astrophysical systems. The simulations, performed using the fluid-kinetic spectral plasma solver (sps) code, demonstrate that the turbulent cascade in such regimes can reach scales smaller than the electron inertial scale, and results in the formation of electron-scale current sheets (ESCS). Statistical analysis of the geometrical properties of the detected ESCS is performed using an algorithm based on the medial axis transform. A typical half-thickness of the current sheets is found to be on the order of electron inertial length or below, while their half-length falls between the electron and ion inertial length. The pressure–strain interaction, used as a measure of energy dissipation, exhibits high intermittency, with the majority of the total energy exchange occurring in current structures occupying approximately 20 per cent of the total volume. Some of the current sheets corresponding to the largest pressure–strain interaction are found to be associated with Alfvénic electron jets and magnetic configurations typical of reconnection. These reconnection candidates represent about 1 per cent of all the current sheets identified.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Physics prospects with MeV neutrino-argon charged current interactions using enhanced photon detection in future LArTPCs

Here, we investigate MeV-scale electron neutrino charged current interactions in a liquid argon time projection chamber (LArTPC) equipped with an enhanced photon detection system. Using simulations of deposited energy in charge and light calorimetry, we explore the potential for dual calorimetric neutrino energy reconstruction. We found energy reconstruction based on light-only calorimetry has a better resolution than combined charge and light calorimetry when hadrons are produced in these events. Meanwhile, enhanced light detection offers improved nanosecond timing resolution and broad optical coverage, enabling neutron tagging and identification of delayed low-energy gamma emissions. These advancements open new avenues in low-energy neutrino physics in next-generation LArTPCs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Systematically localizable operators for quantum simulations of quantum field theories

Correlations and measures of entanglement in ground-state wave functions of relativistic quantum field theories are spatially localized over length scales set by the mass of the field. We utilize this localization to design digital quantum circuits for preparing the ground states of lattice scalar quantum field theories. Controlled rotations that are exponentially localized in their position-space extent are found to provide exponentially convergent wave function fidelity. These angles scale with the correlation between sites and the classical two-point correlation function, as opposed to the more localized mutual information or the hyperlocalized (1×1)-site negativity. We anticipate that further investigations will uncover quantum circuit designs with controlled rotations dictated by the measures of entanglement. This work is expected to impact quantum simulations of systems of importance to nuclear physics, high-energy physics, and basic energy sciences research.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Collective charge excitations and the metal-insulator transition in the square lattice Hubbard-Coulomb model

Here in this article, we discuss the nontrivial collective charge excitations (plasmons) of the extended square lattice Hubbard model. Using a fully nonperturbative approach, we employ the hybrid Monte Carlo algorithm to simulate the system at half-filling. A modified Backus-Gilbert method is introduced to obtain the spectral functions via numerical analytic continuation. We directly compute the single-particle density of states which demonstrates the formation of Hubbard bands in the strongly correlated phase. The momentum-resolved charge susceptibility also is computed on the basis of the Euclidean charge-density-density correlator. In agreement with previous extended dynamical mean-field theory studies, we find that, at high strength of the electron-electron interaction, the plasmon dispersion develops two branches.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2018

Nonadiabatic quantum molecular dynamics with detailed balance

In this paper, we present an approach for carrying out nonadiabatic molecular dynamics simulations of systems in which nonadiabatic transitions arise from the coupling between the classical atomic motions and a quasicontinuum of electronic quantum states. Such conditions occur in many research areas, including chemistry at metal surfaces, radiation damage of materials, and warm-dense-matter physics. The classical atomic motions are governed by stochastic Langevin-like equations, while the quantum electron dynamics is described by a master equation for the populations of the electronic states. These working equations are obtained from a first-principles derivation. Remarkably, unlike the widely used Ehrenfest and surface-hopping methods, the approach naturally satisfies the principle of detailed balance at equilibrium and therefore can describe the evolution to thermal equilibrium from an arbitrary initial state. Lastly, a practical algorithm is cast in the form of the widely used fewest-switches surface-hopping algorithm but with switching probabilities that are not specified ad hoc like in the standard algorithm but are instead derived.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Lattice model of heavy-light three-body system

We present a study of a $1+1$ dimensional heavy-light three-body system in finite volume. The heavy-light system is simulated by a coupled-channel $\phi^4$ type lattice model, and both ground state and excited states of multiparticle energy spectra are measured on various lattices. The lattice simulation data analysis is performed based on variational approach.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗