Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation models”

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,405 records · Page 78

Stellar feedback-regulated black hole growth: driving factors from nuclear to halo scales

ABSTRACT Several recent simulations of galaxy formation predict two main phases of supermassive black hole (BH) accretion: an early, highly intermittent phase (during which BHs are undermassive relative to local scaling relations), followed by a phase of accelerated growth. We investigate physical factors that drive the transition in BH accretion in cosmological zoom-in simulations from the FIRE project, ranging from dwarf galaxies to galaxies sufficiently massive to host luminous quasars. The simulations model multichannel stellar feedback, but neglect AGN feedback. We show that multiple physical properties, including halo mass, galaxy stellar mass, and depth of the central gravitational potential correlate with accelerated BH fuelling: constant thresholds in these properties are typically crossed within ∼0.1 Hubble time of accelerated BH fuelling. Black hole masses increase sharply when the stellar surface density in the inner 1 kpc crosses a threshold $\Sigma^\star _{1\,\rm kpc}\approx 10^{9.5} \, {\rm M_{\odot }}\,{\rm kpc}^{-2}$, a characteristic value above which gravity prevents stellar feedback from ejecting gas, and similar to the value above which galaxies are observed to quench. We further show that accelerated BH growth correlates with the emergence of long-lived thin gas discs, as well as with virialization of the inner circumgalactic medium. The halo mass Mhalo ∼ 1012 M⊙ and stellar mass M* ∼ 1010.5 M⊙ at which BH growth accelerates correspond to ∼L⋆ galaxies. The fact that stellar feedback becomes inefficient at ejecting gas from the nucleus above this mass scale may play an important role in explaining why AGN feedback appears to be most important in galaxies above L⋆.

Astronomy & Astrophysics↗

Measurements of the z > 5 Lyman-α forest flux autocorrelation functions from the extended XQR-30 data set

We present the first observational measurements of the Lyman-α (Ly α) forest flux autocorrelation functions in ten redshift bins from 5.1 ≤ z ≤ 6.0. We use a sample of 35 quasar sightlines at z > 5.7 from the extended XQR-30 data set; these data have signal-to-noise ratios of >20 per spectral pixel. We carefully account for systematic errors in continuum reconstruction, instrumentation, and contamination by damped Ly α systems. With these measurements, we introduce software tools to generate autocorrelation function measurements from any simulation. Our measurements of the smallest bin of the autocorrelation function increase with redshift when normalizing by the mean flux, $\langle{F}\rangle$. This increase may come from decreasing $\langle{F}\rangle$ or increasing mean free path of hydrogen-ionizing photons, λmfp. Recent work has shown that the autocorrelation function from simulations at z > 5 is sensitive to λmfp, a quantity that contains vital information on the ending of reionization. For an initial comparison, we show our autocorrelation measurements with simulation models for recently measured λmfp values and find good agreements. Further work in modelling and understanding the covariance matrices of the data is necessary to get robust measurements of λmfp from this data.

79 ASTRONOMY AND ASTROPHYSICS↗

Investigating event-shape methods in the search for the chiral magnetic effect in relativistic heavy ion collisions

The chiral magnetic effect (CME) is a phenomenon in which electric charge is separated by a strong magnetic field from local domains of chirality imbalance and parity violation in quantum chromodynamics. The CME-sensitive observable, the charge-dependent three-point azimuthal correlator Δ⁢𝛾 , is contaminated by a major physics background proportional to the particle's elliptic flow anisotropy 𝑣 2 . Event-shape engineering (ESE) binning events in dynamical fluctuations of 𝑣 2 and event-shape selection (ESS) binning events in statistical fluctuations of 𝑣 2 are two methods to search for the CME by projecting Δ⁢𝛾 to the measured anisotropy 𝑣 2 = 0 intercept. Here, we conduct a systematic study of these two methods using physics models as well as toy model simulations. It is observed that the ESE method fulfills the general premise of measuring the CME but is statistically hungry. It is found that the intercept from the ESS method depends on the details of the event content, such as the mixtures of background-contributing sources, because of statistical fluctuations of intertwining variables used in the method, and is thus not practically useful to measure the CME.

Relativistic heavy-ion collisions↗

Two-dimensional multifunctional materials from endohedral fullerenes

A multifunctional two-dimensional (2D) material is theoretically predicted based on systematic ab initio calculations and model simulations for the honeycomb lattice of endohedral fullerene W@C 28 molecules. It has structural bistability, ferroelectricity, multiple magnetic phases, and excellent valley characters and can be easily functionalized by the proximity effect with magnetic isolators such as MnTiO 3 . Furthermore, we may also manipulate the valley Hall and spin transport properties by selectively switching a few W@C 28 molecules to the metastable phase. These conclusions contribute to integrating different functions in a single 2D material for technological innovations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Back-to-back relative-excess observable to identify the chiral magnetic effect

Background: The chiral magnetic effect (CME) is extensively studied in heavy-ion collisions at RHIC and LHC. In the commonly used reaction-plane (RP)-dependent, charge-dependent azimuthal correlator ($Δ_γ$), both the close and back-to-back pairs are included. Many backgrounds contribute to the close pairs (e.g., resonance decays, jet correlations), whereas the back-to-back pairs are relatively free of those backgrounds. Purpose: In order to reduce those backgrounds, we propose a new observable, which only focuses on the back-to-back pairs, namely, the relative back-to-back opposite-sign over same-sign pair excess ($r_{BB}$) as a function of the pair azimuthal orientation with respect to the RP ($φ_{BB}$). Methods: We use analytical calculations and toy model simulations to demonstrate the sensitivity of $r_{BB}(φ_{BB})$ to the CME and its insensitivity to backgrounds. Results: With finite CME, the $φ_{BB}$ distribution of $r_{BB}$ shows a clear characteristic modulation. Its sensitivity to background is significantly reduced compared to the previous $Δ_γ$ observable. The simulation results are consistent with our analytical calculations. Conclusions: In this work, our studies demonstrate that the $r_{BB}(φ_{BB})$ observable is sensitive to the CME signal and rather insensitive to the resonance backgrounds.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Photo response of 164 Dy

Background: Little data is available for the pygmy dipole resonance (PDR) in axially deformed nuclei. Photon-scattering experiments are complicated by high level densities in the PDR region and the small energy difference of transitions to the ground state and to excited states. Purpose: In this work, we report on an experimental study of the low-energy dipole strength distribution of the well-deformed nucleus 164 Dy between 4.0–7.7 MeV. Methods: The low-lying photoresponse of 164 Dy has been investigated using the method of nuclear resonance fluorescence using a quasimonochromatic linearly polarized $\textit{γ}$-ray beam in the energy range of 4.0–7.7 MeV in steps of 0.2 MeV. Results: For excitation energies between 4 MeV and 5 MeV, sufficiently low level densities allow for the identification of individual states, including level energies, reduced transition widths and branching ratios. Energy-averaged mean decay branching ratios, mean population ratios and partial absorption cross sections were determined above 5 MeV up to the neutron-separation threshold at 7.7 MeV. A Lorentzian-shaped enhancement of the partial photo absorption cross section followed by decays back to the ground-state band is found at 6.10(5) MeV with a width of 0.77(23) MeV. A comparison with results from complementary measurements is performed using the framework of the statistical model. Conclusions: The experimental results for the mean population ratios deviate systematically from the statistical model simulation by 30(6)%. However, they are in agreement within one standard deviation of the simulation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Angular momentum effects in fission

Background: The role of angular momentum in fission has long been discussed but the observable effects are difficult to quantify. Purpose: Here we discuss a variety of effects associated with angular momentum in fission and present quantitative illustrations. Methods: We employ the fission simulation model freya , which is well suited for this purpose because it obeys all conservation laws, including linear and angular momentum conservation at each step of the process. We first discuss the implementation of angular momentum in freya and then assess particular observables, including various correlated observables. We also study potential effects of neutron-induced fission of the low-lying isomeric state of U 235 relative to the ground state. Results: The fluctuations inherent in the fission process ensure that the spin of the initial compound nucleus has only a small influence on the fragment spins, which are therefore nearly uncorrelated. There is a marked correlation between the spin magnitude of the fission fragments and the photon multiplicity. We also consider the dynamical anisotropy caused by the rotation of an evaporating fragment and study especially the distribution of the projected neutron-neutron opening angles, showing that while it is dominated by the effect of the evaporation recoils, it is possible to extract the signal of the dynamical anisotropy by means of a Fourier decomposition. Finally, we note that the use of an isomeric target, U 235 m ( n th ,f), may enhance the symmetric yields and can thus result in higher neutron multiplicities for low total fragment kinetic energies. Conclusions: While the initial angular momentum of the fissioning nucleus tends to have little effect on the observables, those of the produced fragments influence the emitted neutrons and photons in a significant and correlated manner which may be exploited experimentally to elucidate the fission process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Ratios of collective flow observables in high-energy isobar collisions are insensitive to final-state interactions

The ratios of bulk observables, such as harmonic flow v 2 and v 3 , between high-energy 96 Ru + 96 Ru and 96 Zr + 96 Zr collisions were recently argued to be a clean probe of the nuclear structure differences between 96 Ru and 96 Zr . Using a transport model simulation of isobar collisions, we quantify this claim from the dependence of the ratios v 2,Ru /v 2, Zr and v 3,Ru /v 3,Zr on various final-state effects, such as the shear viscosity, hadronization, and hadronic cascade. Although the v 2 and v 3 change by more than 50% when varying the final-state effects, the ratios are unchanged. In addition, these ratios are independent of the transverse momentum p T and hadron species, despite of up to a factor of 2 change in v n . Further, the ratio of mean transverse momentum $\langle$p T $\rangle$ is found to be controlled by the nuclear skin and nuclear radius but is only slightly impacted by the final-state effects. Therefore, these isobar ratios serve as a clean probe of the initial condition of the quark-gluon plasma, which, in turn, is controlled by the collective structure of the colliding nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Neural network emulation of flow in heavy-ion collisions at intermediate energies

Applications of new techniques in machine learning are speeding up progress in research in various fields. In this work, we construct and evaluate a deep neural network (DNN) to be used within a Bayesian statistical framework as a faster and more reliable alternative to the Gaussian process (GP) emulator of an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions at intermediate beam energies. We found strong evidence of the DNN being able to emulate the IBUU simulator's prediction on the strengths of protons' directed and elliptical flow very efficiently even with small training datasets and with accuracy about ten times higher than the GP. Here, limitations of our present work and future improvements are also discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Evolution of in-medium baryon-baryon scattering cross sections and stiffness of dense nuclear matter from Bayesian analyses of FOPI proton-flow excitation functions

Within a Bayesian statistical framework using a Gaussian process (GP) emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions with momentum-independent Skyrme interactions, we infer from the proton directed and elliptical flow in mid-central Au + Au reactions, measured at beam energies from 150 to 1200 MeV/nucleon by the FOPI Collaboration, the posterior probability distribution functions (PDFs) of the in-medium baryon-baryon scattering cross section (BBSCS) modification factor 𝑋 (with respect to their free-space values) and the stiffness parameter 𝐾 of dense nuclear matter. Here, we find that the most probable value of 𝑋 evolves from around 0.7 to 1.0 as the beam energy 𝐸 beam /𝐴 increases. On the other hand, the posterior PDF⁡(𝐾) may have dual peaks having roughly the same height or extended shoulders at high 𝐾 values. More quantitatively, the posterior PDF⁡(𝐾) changes from having a major peak around 220 MeV, characterizing a soft EOS in the reaction at 𝐸 beam /𝐴 = 150 MeV, to one that peaks around 320 MeV, indicating a stiff EOS in the reactions at 𝐸 beam /𝐴 higher than about 600 MeV. The transition from soft to stiff happens in mid-central Au + Au reactions at beam energies around 250 MeV/nucleon in which 𝐾 = 220 MeV and 𝐾 = 320 MeV are approximately equally probable. Altogether, the FOPI proton-flow excitation function data indicate a gradual hardening of hot and dense nuclear matter as its density and temperature increase in reactions with higher beam energies.

nuclear matter↗

Learning efficient erasure protocols for an underdamped memory

Here we apply evolutionary reinforcement learning to a simulation model to identify efficient time-dependent erasure protocols for a physical realization of a 1-bit memory using an underdamped mechanical cantilever. We show that these protocols, when applied to the cantilever in the laboratory, are considerably more efficient than our best hand-designed protocols. The learned protocols allow reliable high-speed erasure by minimizing the heating of the memory during its operation. More generally, the combination of methods used here opens the door to the rational design of efficient protocols for various physics applications.

74 ATOMIC AND MOLECULAR PHYSICS↗

Coexistence of Surface Ferromagnetism and a Gapless Topological State in MnBi 2 Te 4

Surface magnetism and its correlation with the electronic structure are critical to understanding the topological surface state in the intrinsic magnetic topological insulator MnBi 2 Te 4 . Here, using static and time resolved angle-resolved photoemission spectroscopy (ARPES), we find a significant ARPES intensity change together with a gap opening on a Rashba-like conduction band. Comparison with a model simulation strongly indicates that the surface magnetism on cleaved MnBi 2 Te 4 is the same as its bulk state. Finally, the inability of surface ferromagnetism to open a gap in the topological surface state uncovers the novel complexity of MnBi 2 Te 4 that may be responsible for the low quantum anomalous Hall temperature of exfoliated MnBi 2 Te 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Sparsity-Independent Lyapunov Exponent in the Sachdev-Ye-Kitaev Model

The saturation of a recently proposed universal bound on the Lyapunov exponent has been conjectured to signal the existence of a gravity dual. This saturation occurs in the low-temperature limit of the dense Sachdev-Ye-Kitaev (SYK) model, N Majorana fermions with q body ( q > 2 ) infinite-range interactions. We calculate certain out-of-time-order correlators (OTOCs) for N ≤ 64 fermions for a highly sparse SYK model and find no significant dependence of the Lyapunov exponent on sparsity up to near the percolation limit where the Hamiltonian breaks up into blocks. This provides strong support to the saturation of the Lyapunov exponent in the low-temperature limit of the sparse SYK. A key ingredient to reaching N = 64 is the development of a novel quantum spin model simulation library that implements highly optimized matrix-free Krylov subspace methods on graphical processing units. This leads to a significantly lower simulation time as well as vastly reduced memory usage over previous approaches, while using modest computational resources. Strong sparsity-driven statistical fluctuations require both the use of a much larger number of disorder realizations with respect to the dense limit and a careful finite size scaling analysis. The saturation of the bound in the sparse SYK points to the existence of a gravity analog that would enlarge substantially the number of field theories with this feature. Published by the American Physical Society 2024

Physics↗

Characterizing the Oscillatory Properties of Bulk Electric Systems

This paper presents a process for characterizing the oscillatory dynamics of a large bulk power system. As a demonstration, the process is applied to the Western Interconnection of North America. Several complementary analysis approaches, both new and existing, are employed to provide a comprehensive understanding of the oscillatory properties of the system. Established modal analysis techniques based on ringdown and mode-meter algorithms are utilized. In addition, we derive and apply methods based on spectral correlation analysis to identify modal frequencies, distinguish between modes that are closely spaced in frequency, and determine locations at which the modes are observable. Critical interarea modes are identified and characterized using actual-system synchrophasor measurements taken over several years of operation in concert with industry-standard simulation models. This includes 145 hours of PMU data and two planning base cases.

42 ENGINEERING↗

Synchronverter-based Control of Multi-Port Autonomous Reconfigurable Solar Plants (MARS)

Control of integrated photovoltaic (PV) plants with energy storage systems (ESSs) has become an important research and development topic in recent times. In this context, a Multi-port Autonomous Reconfigurable Solar (MARS) plant that integrates PV and ESS to alternating current transmission grid and high-voltage direct current (HVdc) link is studied in this paper. With penetration of power electronic based resources in the grid, the grid’s capability to recover from frequency or voltage disturbances are reduced. Therefore, one of the vital objectives of any new grid integrated power electronic resource is to provide advanced control functions like voltage and frequency support to the grid during disturbances. In this research work, a detailed implementation of a synchronverter-based control algorithm of MARS is presented. The proposed control algorithm and the MARS control architecture are evaluated through simulations on PSCAD/EMTDC simulation platform to showcase the performance in different operating conditions. In addition, they are evaluated in Opal-RT offline simulation models which can also be used to perform control-hardware-in-the-loop (cHIL) tests.

Marthi, Phani Ratna Vanamali↗

Real-time Simulation Framework for Hardware-in-the-Loop Testing of Multi-port Autonomous Reconfigurable Solar Power Plant (MARS)

Multiport autonomous reconfigurable solar power plant (MARS) has been proposed for integrated development of photovoltaic (PV) and energy storage system (ESS) that can connect to high-voltage direct current (HVdc) and alternating current (ac) transmission grid. To de-risk the development of this complex integrated system that consists of hundreds to thousands of power electronics modules, a controller hardware-in-the-loop (cHIL) test setup will be extremely beneficial. The cHIL testing can be used for evaluation of modules as well as the hierarchical control system in MARS. With the unique configuration of power electronics modules in MARS, it becomes necessary to develop custom-designed real-time simulation models in the cHIL setup in absence of off-the-shelf models. In this paper, high-fidelity dynamic model of MARS, control algorithms at the lower level, and required communication algorithms are developed and optimized for real-time performance in the cHIL setup. Real-time experimental results from the cHIL are provided.

Dong, Zerui↗

A High-Fidelity Cyber-Physical Testbed-Based Benchmarking Dataset For Testing Operational Technology Specific Intrusion Detection Systems

Quality datasets serve a critical purpose in cyber security research. Data is needed to understand system behavior and develop security controls to protect critical systems. However, for critical infrastructure operational environments there is a lack of available datasets to study because of the high cost and specialized capabilities necessary to generate them. This paper documents the development of a dataset of high fidelity hardware in the loop laboratory simulated models of electric and natural gas distribution systems with real cyber attack test cases. A deep dive discussion for the experimental setup and controls for generating the data is provided along with observations from using the data in evaluating intrusion detection approaches.

Ashok, Aditya↗

Using Image Processing Techniques to Identify and Quantify Spatiotemporal Carbon Cycle Extremes

Rising atmospheric carbon dioxide due to human activities through fossil fuel emissions and land use changes have increased climate extremes such as heat waves and droughts that have led to and are expected to increase the occurrence of carbon cycle extremes. Carbon cycle extremes represent large anomalies in the carbon cycle that are associated with gains or losses in carbon uptake. Carbon cycle extremes could be continuous in space and time and cross political boundaries. Here, we present a methodology to identify large spatiotemporal extremes (STEs) in the terrestrial carbon cycle using image processing tools for feature detection. We characterized the STE events based on neighborhood structures that are three-dimensional adjacency matrices for the detection of spatiotemporal manifolds of carbon cycle extremes. We found that the area affected and carbon loss during negative carbon cycle extremes were consistent with continuous neighborhood structures. In the gross primary production data we used, 100 carbon cycle STEs accounted for more than 75% of all the negative carbon cycle extremes. This paper presents a comparative analysis of the magnitude of carbon cycle STEs and attribution of those STEs to climate drivers as a function of neighborhood structures for two observational datasets and an Earth system model simulation.

Sharma, Bharat↗