Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “initial state”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Chiral effects and Joule heating in hot and dense matter

Initial states of dense matter with nonzero electron chiral imbalance could potentially give rise to strong magnetic fields through chiral plasma instability. Previous work indicated that unless chiral chemical potential is as large as the electron vector chemical potential, the growth of magnetic fields due to the instability is washed out by chirality flipping rate enabled by electron mass. We reexamine this claim in a broader range of parameters and find that at higher temperatures the hierarchy is reversed supporting a growing magnetic field for an initial electron chiral chemical potential much smaller than the electron vector chemical potential. Further, we identify a qualitatively new effect relevant for magnetized hot and dense medium where chiral magnetic effect (CME) sourced by density fluctuation acts as a powerful source of Joule heating. Remarkably, even modest chiral chemical potentials (keV) in such environment can deposit energy densities set by the QCD scale in a relatively short time of the order of a few milliseconds or seconds. We speculate how this mechanism makes CME-driven Joule heating a potentially critical ingredient in the dynamics of turbulent density fluctuation of supernovae and neutron star mergers.

Particle astrophysics↗
DOE OSTI2026

Direct estimation of the density of states for fermionic systems

Simulating time evolution is one of the most natural applications of quantum computers and is thus one of the most promising prospects for achieving practical quantum advantage. Here, we develop quantum algorithms to extract thermodynamic properties by estimating the density of states (DOS), which is a central object in quantum statistical mechanics. We introduce several key innovations that significantly improve the practicality and extend the generality of previous techniques. First, our approach allows one to estimate the DOS only for a specific subspace of the full Hilbert space. This is crucial for fermionic systems, since both canonical and grand canonical ensemble thermal equilibrium properties depend on subspaces of fixed number. Second, in our approach, by time evolving very simple, random initial states, such as randomly chosen computational basis states, we can exactly recover the DOS on average. Third, due to circuit-depth limitations, we only reconstruct the DOS up to a convolution with a Gaussian window—thus all imperfections that shift the energy levels by less than the width of the convolution window will not significantly affect the estimated DOS. For these reasons, we find the approach is a promising candidate for early quantum advantage as even short-time, noisy dynamics can yield a semiquantitative reconstruction of the DOS (convolution with a broad Gaussian window), while early fault-tolerant devices will likely enable higher-resolution DOS reconstruction through longer time evolutions. We demonstrate the practicality of our approach in representative Fermi-Hubbard and spin models and indeed find that our approach is highly robust against algorithmic errors in the time evolution and against gate noise. We further demonstrate that our approach is compatible with noisy intermediate-scale quantum (NISQ) computing NISQ-friendly variational techniques, introducing and leveraging a technique for variational time evolution.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Efficient state preparation for the Schwinger model with a theta term

We present a comparison of different quantum state preparation algorithms and their overall efficiency for the Schwinger model with a theta term. While adiabatic state preparation is proved to be effective, in practice it leads to large gate counts to prepare the ground state. The quantum approximate optimization algorithm (QAOA) provides excellent results while keeping the counts small by design, at the cost of an expensive classical minimization process. We introduce a “blocked” modification of the Schwinger Hamiltonian to be used in the QAOA that further decreases the length of the algorithms as the size of the problem is increased. The rodeo algorithm (RA) provides a powerful tool to efficiently prepare any eigenstate of the Hamiltonian, as long as its overlap with the initial guess is large enough. We obtain the best results when combining the blocked QAOA ansatz and the RA, as this provides an excellent initial state with a relatively short algorithm without the need to perform any classical steps for large problem sizes. Published by the American Physical Society 2025

Bazavov, Alexei (ORCID:0000000321411901)↗
DOE OSTI2020

Disentangling contributions to small-system collectivity via scans of light nucleus-nucleus collisions

The observation of multiparticle azimuthal correlations in high-energy small-system collisions has led to intense debate on its origin and the possible coexistence from two competing theoretical scenarios: one based on initial-state intrinsic momentum anisotropy, and the other based on final-state interaction model (FSM) collective response to the collision geometry. To complement the previous scan of asymmetric collision systems ( p +Au, d +Au, and He+Au), we propose a scan of small symmetric collision systems at the Brookhaven National Laboratory Relativistic Heavy Ion Collider, such as C+C, O+O, Al+Al, and Ar+Ar$\sqrt{s_\text{NN}}$ = 0.2 TeV , to provide further insights in disentangling contributions from these two scenarios. These symmetric small systems have the advantage of providing a better controlled initial geometry dominated by the average shape of the overlap region as opposed to fluctuation-driven geometries in asymmetric systems. In this work, a transport model is employed to investigate the expected geometry response in the FSM scenario. Different trends of elliptic flow with increasing charge particle multiplicity are observed between symmetric and asymmetric systems, whereas triangular flow appears to show a similar behavior. Furthermore, a comparison of O+O collisions at $\sqrt{s_\text{NN}}$ = 0.2 TeV and at $\sqrt{s_\text{NN}}$ = 2.76 ₋ 7 TeV as proposed at the CERN Large Hadron Collider, provides a unique opportunity to disentangle the collision geometry effects at the nucleon level from those arising from subnucleon fluctuations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Discrete Time-Crystalline Order Enabled by Quantum Many-Body Scars: Entanglement Steering via Periodic Driving

The control of many-body quantum dynamics in complex systems is a key challenge in the quest to reliably produce and manipulate large-scale quantum entangled states. Recently, quench experiments in Rydberg atom arrays [Bluvstein et al. Science 371, 1355 (2021)] demonstrated that coherent revivals associated with quantum many-body scars can be stabilized by periodic driving, generating stable subharmonic responses over a wide parameter regime. We analyze a simple, related model where these phenomena originate from spatiotemporal ordering in an effective Floquet unitary, corresponding to discrete time-crystalline behavior in a prethermal regime. Unlike conventional discrete time crystals, the subharmonic response exists only for Néel-like initial states, associated with quantum scars. We predict robustness to perturbations and identify emergent timescales that could be observed in future experiments. Our results suggest a route to controlling entanglement in interacting quantum systems by combining periodic driving with many-body scars.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS2014

A Self-Stabilizing Hybrid-Fault Tolerant Synchronization Protocol

In this report we present a strategy for solving the Byzantine general problem for self-stabilizing a fully connected network from an arbitrary state and in the presence of any number of faults with various severities including any number of arbitrary (Byzantine) faulty nodes. Our solution applies to realizable systems, while allowing for differences in the network elements, provided that the number of arbitrary faults is not more than a third of the network size. The only constraint on the behavior of a node is that the interactions with other nodes are restricted to defined links and interfaces. Our solution does not rely on assumptions about the initial state of the system and no central clock nor centrally generated signal, pulse, or message is used. Nodes are anonymous, i.e., they do not have unique identities. We also present a mechanical verification of a proposed protocol. A bounded model of the protocol is verified using the Symbolic Model Verifier (SMV). The model checking effort is focused on verifying correctness of the bounded model of the protocol as well as confirming claims of determinism and linear convergence with respect to the self-stabilization period. We believe that our proposed solution solves the general case of the clock synchronization problem.

Malekpour, Mahyar R.↗
NASA NTRS2009

Examining the Impacts of High-Resolution Land Surface Initialization on Model Predictions of Convection in the Southeastern U.S.

One of the most challenging weather forecast problems in the southeastern U.S. is daily summertime pulse convection. During the summer, atmospheric flow and forcing are generally weak in this region; thus, convection typically initiates in response to local forcing along sea/lake breezes, and other discontinuities often related to horizontal gradients in surface heating rates. Numerical simulations of pulse convection usually have low skill, even in local predictions at high resolution, due to the inherent chaotic nature of these precipitation systems. Forecast errors can arise from assumptions within physics parameterizations, model resolution limitations, as well as uncertainties in both the initial state of the atmosphere and land surface variables such as soil moisture and temperature. For this study, it is hypothesized that high-resolution, consistent representations of surface properties such as soil moisture and temperature, ground fluxes, and vegetation are necessary to better simulate the interactions between the land surface and atmosphere, and ultimately improve predictions of local circulations and summertime pulse convection. The NASA Short-term Prediction Research and Transition (SPORT) Center has been conducting studies to examine the impacts of high-resolution land surface initialization data generated by offline simulations of the NASA Land Informatiot~ System (LIS) on subsequent numerical forecasts using the Weather Research and Forecasting (WRF) model (Case et al. 2008, to appear in the Journal of Hydrometeorology). Case et al. presents improvements to simulated sea breezes and surface verification statistics over Florida by initializing WRF with land surface variables from an offline LIS spin-up run, conducted on the exact WRF domain and resolution. The current project extends the previous work over Florida, focusing on selected case studies of typical pulse convection over the southeastern U.S., with an emphasis on improving local short-term WRF simulations over the Mobile, AL and Miami, FL NWS county warning areas. Future efforts may involve examining the impacts of assimilating remotely-sensed soil moisture data, and/or introducing weekly greenness vegetation fraction composites (as opposed to monthly climatologies) into ol'fline NASA LIS runs. Based on positive impacts, the offline LIS runs could be transitioned into an operational mode, providing land surface initialization data to NWS forecast offices in real time.

Case, Jonathan L.↗
NASA NTRS1993

Diffusion effects on the helium abundance of the solar transition region and corona

The diffusion of helium in the solar transition region is studied by solving the mass and momentum conservation equations for a hydrogen-helium plasma given a representative temperature profile. Steady state solutions show that two distinct atmospheres may result. In cases where the thermal force on alpha-particles is balanced by the partial pressure gradient force, helium is the dominant coronal species. On the other hand, if it is the frictional force between protons and alpha-particles which balances the thermal force on alpha-particles then hydrogen is the major coronal component. In order to explore which of these solutions are attainable within reasonable time scales, the time-dependent equations are solved, starting from an initial state with a uniform helium abundance of 10 percent. The atmosphere as a whole is close to hydrostatic equilibrium, but due the thermal forces the individual elements are not. This force inbalance leads to a differential flow between species. It is found that this differential flow leads to a significant enhancement of the coronal helium abundance. Even for the relatively shallow temperature gradient used the helium abundance in the lower corona increases to 30 percent over a 24 hr period.

Hansteen, Viggo H.↗
DOE OSTI2021

O 16 O 16 collisions at energies available at the BNL Relativistic Heavy Ion Collider and at the CERN Large Hadron Collider comparing α clustering versus substructure

Collisions of . light and heavy nuclei in relativistic heavy-ion collisions have been shown to be sensitive to nuclear structure. With a proposed 16 O 16 O run at the CERN Large Hadron Collider (LHC) and at the BNL Relativistic Heavy Ion Collider (RHIC) we study the potential for finding α clustering in 16 O. Here we use the state-of-the-art iEBE-VISHNU package with 16 O nucleonic configurations from ab initio nuclear lattice simulations. This setup was tuned using a Bayesian analysis on p Pb and PbPb systems. We find that the 16 O 16 O system always begins far from equilibrium and that at LHC and RHIC it approaches the regime of hydrodynamic applicability only at very late times. Finally, by taking ratios of flow harmonics we are able to find measurable differences between α-clustering, nucleonic, and subnucleonic degrees of freedom in the initial state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1985

A Comparison of the Bounded Derivative and the Normal Mode Initialization Methods Using Real Data

Browning et al. (1980) proposed an initialization method called the bounded derivative method (BDI). They used analytical data to test the new method. Kasahara (1982) theoretically demonstrated the equivalence between BDI and the well known nonlinear normal mode initialization method (NMI). The purposes of this study are the extension of the application of BDI to real data and comparison with NMI. The unbalanced initial state (UBD) is data of January, 1979 OOZ which were interpolated from the adjacent sigma levels of the GLAS GCM to the 300 mb surface. The global barotropic model described by Takacs and Balgovind (1983) is used. Orographic forcing is explicitly included in the model. Many comparisons are performed between various quantities. However, we only present a comparison of the time evolution at two grid points A(50 S, 90 E) and B(10 S, 20 E) which represent low and middle latitude locations. To facilitate a more complete comparison an initialization experiment based on the classical balance equation (CBE) was also included.

Semazzi, F. H. M.↗
NASA NTRS2003

The MAP Spacecraft Angular State Estimation After Sensor Failure

This work describes two algorithms for computing the angular rate and attitude in case of a gyro and a Star Tracker failure in the Microwave Anisotropy Probe (MAP) satellite, which was placed in the L2 parking point from where it collects data to determine the origin of the universe. The nature of the problem is described, two algorithms are suggested, an observability study is carried out and real MAP data are used to determine the merit of the algorithms. It is shown that one of the algorithms yields a good estimate of the rates but not of the attitude whereas the other algorithm yields a good estimate of the rate as well as two of the three attitude angles. The estimation of the third angle depends on the initial state estimate. There is a contradiction between this result and the outcome of the observability analysis. An explanation of this contradiction is given in the paper. Although this work treats a particular spacecraft, the conclusions have a far reaching consequence.

Bar-Itzhack, Itzhack Y.↗
NASA NTRS1981

Quasi-linear blocks forced by orography in a hemispheric, quasi-geostrophic barotropic model

Stationary linear perturbation responses to Northern Hemisphere orography are calculated in a quasi-geostrophic barotropic model in solid-body rotation. The stationary mountain torque induced by these perturbations is then used to construct graphical solutions to the steady-state wave, mean-flow interaction problem. It is shown that multiple solutions exist in the system and are near either the forcing equilibrium of the zonal forcing or near the resonance points in the system. Some of these near-resonance solutions have block-like configurations with a confluence zone upstream from a large-amplitude structure consisting of a high at high latitudes and a low at low latitudes. These block-like configurations are shown to be near stable solutions of the system. Time-dependent calculations show that the initial state and the zonal forcing equilibrium are important in determining the long-term time evolution of the system.

Roads, J. O.↗
DOE OSTI2025

Global stellarator coil optimization with quadratic constraints and objectives

Most present stellarator designs are produced by costly two-stage optimization: the first for an optimized equilibrium, and the second for a coil design reproducing its magnetic configuration. Few proxies for coil complexity and forces exist at the equilibrium stage. Rapid initial state finding for both stages is a topic of active research. Most present convex coil optimization codes use the least square winding surface method by Merkel (NESCOIL), with recent improvements in conditioning, regularization, sparsity, and physics objectives. While elegant, the method is limited to modeling the norms of linear functions in coil current. We present QUADCOIL, a global coil optimization method that targets combinations of linear and quadratic functions of the current. It can directly constrain and/or minimize a wide range of physics objectives unavailable in NESCOIL and REGCOIL, including the Lorentz force, magnetic energy, curvature, field-current alignment, and the maximum density of a dipole array. QUADCOIL requires no initial guess and runs nearly $10$ 2 x faster than filament optimization. Integrating it in the equilibrium optimization stage can potentially exclude equilibria with difficult-to-design coils, without significantly increasing the computation time per iteration. QUADCOIL finds the exact, global minimum in a large parameter space when possible, and otherwise finds a well-performing approximate global minimum. It supports most regularization techniques developed for NESCOIL and REGCOIL. We demonstrate QUADCOIL’s effectiveness in coil topology control, minimizing non-convex penalties, and predicting filament coil complexity with three numerical examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Comparison of Interface Growth Models Applied to Rayleigh–Taylor and Richtmyer–Meshkov Instabilities

Sophisticated numerical models that contain fluid interfaces rely upon interface evolution models to approximate the transition to turbulence near interfaces, in the presence of Rayleigh–Taylor (RTI) or Richtmyer–Meshkov (RMI) instability. Semi-analytical models have been developed in recent decades to predict the interface growth from an initial state into the nonlinear regime. Two of these models are considered in this study. They are the Goncharov and the z-models. Both of these interface models have strengths and weaknesses, which are examined here. Both of them have been implemented in the xRAGE compressible flow solver, which models fluid interfaces. The flow solver provides the fluids acceleration as a body force to the interface model. The purpose of such interface model is to evolve the early times of interface position as a subgrid model within a compressible flow simulation in order to then initialize a turbulence model. In this work, the interface models are assessed and compared for their evolution of RTI and RMI. The z-model performed better than the Goncharov model for all cases that were explored.

42 ENGINEERING↗
DOE OSTI2024

Dark Matter searches with photons at the LHC

We unveil blind spot regions in dark matter (DM) direct detection (DMDD), for weakly interacting massive particles with a mass around a few hundred GeV that may reveal interesting photon signals at the LHC. We explore a scenario where the DM primarily originates from the singlet sector within the Z 3 -symmetric Next-to-Minimal Supersymmetric Standard Model (NMSSM). A novel DMDD spin-independent blind spot condition is revealed for singlino-dominated DM, in cases where the mass parameters of the higgsino and the singlino-dominated lightest supersymmetric particle (LSP) exhibit opposite relative signs (i.e., κ < 0), emphasizing the role of nearby bino and higgsino-like states in tempering the singlino-dominated LSP. Additionally, proximate bino and/or higgsino states can act as co-annihilation partner(s) for singlino-dominated DM, ensuring agreement with the observed relic abundance of DM. Remarkably, in scenarios involving singlino-higgsino co-annihilation, higgsino-like neutralinos can distinctly favor radiative decay modes into the singlino-dominated LSP and a photon, as opposed to decays into leptons/hadrons. In exploring this region of parameter space within the singlino-higgsino compressed scenario, we study the signal associated with at least one relatively soft photon alongside a lepton, accompanied by substantial missing transverse energy (E T ) and a hard initial state radiation jet at the LHC. In the context of singlino-bino co-annihilation, the bino state, as the next-to-LSP, exhibits significant radiative decay into a soft photon and the LSP, enabling the possible exploration at the LHC through the triggering of this soft photon alongside large E T and relatively hard leptons/jets resulting from the decay of heavier higgsino-like states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A Bayesian Approach for In-Situ Stress Prediction and Uncertainty Quantification for Subsurface Engineering

Many subsurface engineering applications require accurate knowledge of the in-situ state of stress for their safe design and operation. Existing methods to meet this need primarily include field measurements for estimating one or more of the principal stresses from a borehole, or optimization methods for constructing a 3D geomechanical model in terms of geophysical measurements. These methods, however, often contain considerable uncertainty in estimating the state of stress. Here, in this paper, we build on a Bayesian approach to quantify uncertainty in stress estimations for subsurface engineering applications. This approach can provide an estimate of the 3D distribution of stress throughout the volume of interest and provide an estimate of the uncertainty arising from the stress measurement, the rheology parameters, and a paucity of measurements. The value of this approach is demonstrated using stress measurements from the In Salah carbon storage site, which was one of the world’s first industrial carbon capture and storage projects. This demonstration shows the application of this Bayesian approach for estimating the initial state of stress for In Salah and quantifying the uncertainty in the estimated stress. Also, an assessment of a maximum injection pressure to prevent geomechanical risks from CO 2 injection pressures is provided in terms of the probability distribution of the minimum principal stress quantified by the approach. With the In Salah case study, this paper demonstrates that using the Bayesian approach can provide additional insights for site explorations and/or project operations to make informed-site decisions for subsurface engineering applications.

58 GEOSCIENCES↗
DOE OSTI2021

One-directional polarization transport in electron/nuclear spin chains with loss and gain

Understanding the joint dynamics of electron and nuclear spins is central to core concepts in solid-state magnetic resonance such as spin-lattice relaxation and dynamic nuclear polarization, but a generalization that capitalizes on competing polarization loss and gain channels is still lacking. Here, we study theoretically the non-Hermitian dynamics of hybrid electron/nuclear spin systems in the simultaneous presence of electron spin pumping and spin-lattice relaxation. Focusing on periodic, one-dimensional chains, we find that, by adjusting the electron spin pumping to a critical level, it is possible to steer the flow of nuclear polarization to create site-dependent distributions where either end of the array polarizes in opposite ways, irrespective of the initial state. By contrast, we show that ring-like patterns—where the limit nuclear polarization is uniform—exhibit a non-decaying, externally driven nuclear spin current. Interestingly, cyclic magnetic field modulation can render these processes largely robust to defects in the chain, a response featuring some interesting similarities—and differences—with recent findings in other non-Hermitian physical platforms.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Search for the Z γ decay mode of new high-mass resonances in pp collisions at s = 13 TeV with the ATLAS detector

This letter presents a search for narrow, high-mass resonances in the Zγ final state with the Z boson decaying into a pair of electrons or muons. The $\sqrt{s}$ =13 TeV pp collision data were recorded by the ATLAS detector at the CERN Large Hadron Collider and have an integrated luminosity of 140 fb -1 . The data are found to be in agreement with the Standard Model background expectation. Upper limits are set on the resonance production cross section times the decay branching ratio into Zγ. For spin-0 resonances produced via gluon-gluon fusion, the observed limits at 95% confidence level vary between 65.5 fb and 0.6 fb, while for spin-2 resonances produced via gluon-gluon fusion (or quark-antiquark initial states) limits vary between 77.4 (76.1) fb and 0.6 (0.5) fb, for the mass range from 220 GeV to 3400 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗