Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mean current vector”

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.

29 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Pseudogauge dependence of the spin polarization and of the axial vortical effect

The mean spin polarization vector of spin 1/2 particles in a relativistic fluid at local thermal equilibrium (LTE) in different pseudogauges (PGs), i.e., with different choices for the decomposition of orbital and spin angular momentum, is obtained. The spin polarization obtained in the canonical PG differs from the one obtained in the Belinfante PG. It is found that this difference can not be written by replacing the thermal vorticity with the spin potential in the usual polarization formula. Other PG choices affect the contribution of thermal shear to the spin polarization. So, the choice of a PG affects the predictions for the polarization measured in heavy-ion collisions. In general, it is shown that the Wigner function of a noninteracting Dirac field at LTE is affected by the PG transformations. Explicit expressions of the mean axial current are also calculated and it is found that even the axial vortical effect conductivity depends on the PG.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Incommensurate Spin Fluctuations in the Spin-Triplet Superconductor Candidate UTe 2

Spin-triplet superconductors are of extensive current interest because they can host topological state and Majorana fermions important for quantum computation. The uranium-based heavy-fermion superconductor UTe 2 has been argued as a spin-triplet superconductor similar to UGe2, URhGe, and UCoGe, where the superconducting phase is near (or coexists with) a ferromagnetic (FM) instability and spin-triplet electron pairing is driven by FM spin fluctuations. Here we use neutron scattering to show that, although UTe 2 exhibits no static magnetic order down to 0.3 K, its magnetism in the [0,K,L] plane is dominated by incommensurate spin fluctuations near an antiferromagnetic ordering wave vector and extends to at least 2.6 meV. We are able to understand the dominant incommensurate spin fluctuations of UTe 2 in terms of its electronic structure calculated using a combined density-functional and dynamic mean-field theory.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Magnetic anisotropy in single-crystalline antiferromagnetic Mn 2 Au

Multiple recent studies have identified the metallic antiferromagnet Mn 2 ⁢Au to be a candidate for spintronic applications due to apparent in-plane anisotropy, preserved magnetic properties above room temperature, and current-induced Néel vector switching. Crystal growth is complicated by the fact that Mn 2 ⁢Au melts incongruently. We present a bismuth flux method to grow millimeter-scale bulk single crystals of Mn 2 ⁢Au in order to examine the intrinsic anisotropic electrical and magnetic properties. Flux quenching experiments reveal that the Mn 2 ⁢Au crystals precipitate below 550°⁢C, about 100⁢°⁢C below the decomposition temperature of Mn 2 ⁢Au. Bulk Mn 2 ⁢Au crystals have a room-temperature resistivity of 16–19 µ⁢Ωcm and a residual resistivity ratio of 41. Mn 2 ⁢Au crystals have a dimensionless susceptibility on the order of 10 –4 (SI units), comparable to calculated and experimental reports on powder samples. Single-crystal neutron diffraction confirms the in-plane magnetic structure. The tetragonal symmetry of Mn 2 ⁢Au constrains the ab-plane magnetic susceptibility to be constant, meaning that χ 100 =χ 110 in the low-field limit, below any spin-flop transition. We find that three measured magnetic susceptibilities χ 100 , χ 110 , and χ 001 are the same order of magnitude and agree with the calculated prediction, meaning the low-field susceptibility of Mn 2 ⁢Au is quite isotropic, despite clear differences in ab-plane and ac-plane magnetocrystalline anisotropy. Mn 2 ⁢Au is calculated to have an extremely high in-plane spin-flop field above 30 T, which is much larger than that of another in-plane antiferromagnet, Fe 2 ⁢As (less than 1 T). Finally, the subtle anisotropy of intrinsic susceptibilities may lead to dominating effects from shape, crystalline texture, strain, and defects in devices that attempt spin readout in Mn 2⁢ Au.

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

Huge ensembles – Part 1: Design of ensemble weather forecasts using spherical Fourier neural operators

Abstract. Simulating low-likelihood high-impact extreme weather events in a warming world is a significant and challenging task for current ensemble forecasting systems. While these systems presently use up to 100 members, larger ensembles could enrich the sampling of internal variability. They may capture the long tails associated with climate hazards better than traditional ensemble sizes. Due to computational constraints, it is infeasible to generate huge ensembles (comprised of 1000–10 000 members) with traditional, physics-based numerical models. In this two-part paper, we replace traditional numerical simulations with machine learning (ML) to generate hindcasts of huge ensembles. In Part 1, we construct an ensemble weather forecasting system based on spherical Fourier neural operators (SFNOs), and we discuss important design decisions for constructing such an ensemble. The ensemble represents model uncertainty through perturbed-parameter techniques, and it represents initial condition uncertainty through bred vectors, which sample the fastest-growing modes of the forecast. Using the European Centre for Medium-Range Weather Forecasts Integrated Forecasting System (IFS) as a baseline, we develop an evaluation pipeline composed of mean, spectral, and extreme diagnostics. With large-scale, distributed SFNOs with 1.1 billion learned parameters, we achieve calibrated probabilistic forecasts. As the trajectories of the individual members diverge, the ML ensemble mean spectra degrade with lead time, consistent with physical expectations. However, the individual ensemble members' spectra stay constant with lead time. Therefore, these members simulate realistic weather states during the rollout, and the ML ensemble passes a crucial spectral test in the literature. The IFS and ML ensembles have similar extreme forecast indices, and we show that the ML extreme weather forecasts are reliable and discriminating. These diagnostics ensure that the ensemble can reliably simulate the time evolution of the atmosphere, including low-likelihood high-impact extremes. In Part 2, we generate a huge ensemble initialized each day in summer 2023, and we characterize the simulations of extremes.

Mahesh, Ankur↗
DOE OSTI2025

Superfluid weight of strongly inhomogeneous superconductors

In this study, we obtain the superfluid weight of strongly inhomogeneous superconductors. We find that, in general, the response of the superconductor’s pairing potential to the perturbing vector potential is essential to obtain the correct value of the superfluid weight. We consider two exemplary cases: the case when strong inhomogeneities in the pairing potential are induced by a periodic potential, and the case when superconducting vortices are induced by an external magnetic field. For both cases we show that the correction to the superfluid weight due to the response of the paring potential to the perturbing vector potential can be significant, it must be included to obtain quantitatively correct results, and that for the case when vortices are present the expression of the superfluid weight that does not include such correction returns qualitatively wrong results. The results are directly relevant for superconducting states in moiré multilayer where the almost inevitable spatial variations of the magic angle and the proximity of superconductivity to anisotropic phases induce inhomogeneities in the pairing potential, and to the efforts to optimally pin superconducting vortices to achieve large superconducting critical currents.

Bogoliubov-de Gennes equations↗
DOE OSTI2021

A compute-bound formulation of Galerkin model reduction for linear time-invariant dynamical systems

This work aims to advance computational methods for projection-based reduced-order models (ROMs) of linear time-invariant (LTI) dynamical systems. For such systems, current practice relies on ROM formulations expressing the state as a rank-1 tensor (i.e., a vector), leading to computational kernels that are memory bandwidth bound and, therefore, ill-suited for scalable performance on modern architectures. This weakness can be particularly limiting when tackling many-query studies, where one needs to run a large number of simulations. This work introduces a reformulation, called rank-2 Galerkin, of the Galerkin ROM for LTI dynamical systems which converts the nature of the ROM problem from memory bandwidth to compute bound. We present the details of the formulation and its implementation, and demonstrate its utility through numerical experiments using, as a test case, the simulation of elastic seismic shear waves in an axisymmetric domain. We quantify and analyze performance and scaling results for varying numbers of threads and problem sizes. In conclusion, we present an end-to-end demonstration of using the rank-2 Galerkin ROM for a Monte Carlo sampling study. We show that the rank-2 Galerkin ROM is one order of magnitude more efficient than the rank-1 Galerkin ROM (the current practice) and about 970 times more efficient than the full-order model, while maintaining accuracy in both the mean and statistics of the field.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Transport Properties of Magnetized High-Energy Density Plasma (Final Project Report)

This report summarizes results of project DE-SC0016159 "Transport Properties of Magnetized High-Energy-Density Plasma,'' which ran from 7/15/2016 - 7/14/2021. The primary objectives of the work, as stated in the original proposal, were to develop a theory that describes transport coefficients in magnetized high energy density plasmas and to test the theory with molecular dynamics (MD) simulations. The science challenge is that strong ion coupling, strong magnetization of electrons, and partial degeneracy of electrons, cause the system to be in a regime that is not well described by current theory. The particular processes that were to be investigated include: ion stopping power, diffusion, electron-ion temperature relaxation, thermal conduction and viscosity. All of the primary research objectives were accomplished during the course of this work. In addition, some unexpected results were discovered along the way that led to new and productive research directions. This work resulted in 17 publications in well-respected peer-reviewed journals (primarily Physics of Plasmas and Physical Review E), and 2 more are being prepared for publication. A few of these were selected as Editor's Picks and one was published as a Rapid Communication. Highlights of the research results include: Transport phase space: The first identification of the parameter space that defines fundamental transport regimes in terms of the Coulomb coupling and magnetization parameters. Mean force kinetic theory: Systematic derivation of a kinetic theory for strongly coupled plasmas based on a new expansion parameter of the BBGKY hierarchy. This provided the derivation of a theory that we had previously posed phenomenologically, and also revealed a new term that captures the equation of state properties at all coupling strengths. Transverse friction force: Discovery that the friction force on a test charge in a strongly magnetized plasma includes a component that is perpendicular to the motion of the test charge in the plane defined by the velocity and magnetic field vectors. Gyrofriction force: Discovery that the friction force on a test charge in a plasma that is both strongly magnetized and strongly coupled includes a component of the force in the direction of the Lorentz force. Kinetic theory for warm dense matter: Extension of the mean force kinetic theory to include partial electron degeneracy, so that it applies to warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Seismic moment tensor classification using elliptical distribution functions on the hypersphere

Discrimination of underground explosions from naturally occurring earthquakes and other anthropogenic sources is one of the fundamental challenges of nuclear explosion monitoring. In an operational setting, the number of events that can be thoroughly investigated by analysts is limited by available resources. The capability to rapidly screen out events that can be robustly identified as not being explosions is, therefore, of great potential benefit. Nevertheless, possible mis-classification of explosions as earthquakes currently limits the use of screening methods for verification of test-ban treaties. Moment tensors provide a physics-based classification tool for the characterization of different seismic sources and have enabled the advent of new techniques for discriminating between earthquakes and explosions. Following normalization and projection of their six-degree vectors onto the hypersphere, existing screening approaches use spherically symmetric metrics to determine whether any new moment tensor may have been an explosion. Here, we show that populations of moment tensors for both earthquakes and explosions are anisotropically distributed on the hypersphere. Distributions possessing elliptical symmetry, such as the scaled von Mises–Fisher distribution, therefore provide a better description of these populations than the existing spherically symmetric models. We describe a method that uses these elliptical distributions in combination with a Bayesian classifier to achieve successful classification rates of 99 per cent for explosions and 98 per cent for earthquakes using existing catalogues of events from the western United States. The 1983 May 5 Crowdie underground nuclear test and 2018 July 20 DAG-1 deep-borehole chemical explosion are the only two explosions out of 140 that are incorrectly classified. Application of the method to the 2006–2017 nuclear tests in the Democratic People’s Republic of Korea yields 100 per cent identification rates and we provide a simple routine MTid for general usage. The approach provides a means to rapidly assess the likelihood of an event being an explosion and can be built into monitoring workflows that rely on simultaneously assessing multiple different discrimination metrics.

58 GEOSCIENCES↗
DOE OSTI2023

A Study of the Properties and Dynamics of the Disk of Satellites in a Milky-Way-like Galaxy System

The dynamics of the satellite systems of Milky-Way-like galaxies offer a useful means by which to study the galaxy formation process in the cosmological context. It has been suggested that the currently observed anisotropic distribution of the satellites in such galaxy systems is inconsistent with the concordance ΛCDM cosmology model on the galactic scale if the observed satellites are random samples of the dark matter (DM) sub-halos that are nearly isotropically distributed around the central galaxy. In this study, we present original high-resolution zoom-in studies of central galaxies and satellite systems based upon initial conditions for the DM distribution from the Aquarius simulations but with substantial high-resolution baryon physics added. We find that the galaxy most like the Milky Way in this study does indeed contain a disk of satellites (DOS). Although one galaxy DOS system does not answer the question of how common such disks are, it does allow the opportunity to explore the properties and dynamics of the DOS system. Our investigation centers on the spatial arrangement (distances, angles, etc.) of satellites in this Milky-Way-like galaxy system with a specific emphasis on identifying and analyzing the disk-like structure along with its dynamical and morphological properties. Among the conclusions from this study, we find that the satellites and DM sub-halos in the galaxy simulations are anisotropically distributed. The dynamical properties of the satellites, however, indicate that the direction of the angular momentum vector of the whole satellite system is different from the normal direction of the fitted DOS and from the normal direction of the velocity dispersion of the system. Hence, the fitted DOS appears to be comprised of infalling sub-halos and is not a rotationally supported system.

Zhao, Xinghai↗
DOE OSTI2020

Incorporating the Impacts of Climate Change on Hydrology in a Performance Assessment Model - 20403

The evidence of climate change is increasingly well-documented and impacts should be incorporated in performance assessment studies. The current climate literature provides both observational evidence and climate model projections of climate trends and/or climate change in the late 20. and early 21. centuries for North America and the northeast United States. Probabilistic modeling is a core requirement for quantifying uncertainty and evaluating its impacts. Not evaluating future climate states in a performance assessment because of the existence of uncertainty is contradictory to good modeling practices - the most uncertain issues and parameters require the most attention in effective probabilistic modeling. Excluding climate change limits development of modeling information that could aid in effective decision making. In this work we develop methods to use the output from hydrologic models and analysis of historical aerial imagery to quantify and implement the impacts of climate change on hydrologic processes at a nuclear waste site in West Valley, New York. Specifically, we used the HELP (Hydraulic Performance of Landfill Performance) model to characterize key hydrologic processes under both current and future climate conditions to assess the impacts of changing climate on hydrology. A suite of previous climatic models were reviewed and synthesized to produce a cohesive representation of the current state of knowledge of the impact of climate change on important model inputs such as precipitation. Output from the HELP simulations was coupled to the GoldSim model that was used to develop the Probabilistic Performance Assessment (PPA) approach through the application of a novel 'nearest neighbor' technique. First, several thousand realizations were generated from the HELP model using a Latin Hypercube experimental design to ensure adequate coverage of the parameter space of explanatory variables used to drive HELP. For example, porosity is a physical parameter that is used as an input to both HELP and the GoldSim PA model. We then ran sensitivity analysis (SA) algorithms on the output of HELP for each of the responses of interest. For each predictor, each time we build an SA model we get a different value for the sensitivity index (SI). From the collection of all the SA models, the average was computed among all of the SIs to represent the predictor within the context of the nearest neighbor approach. That is, we conduct SA on each HELP outcome for each scenario. This gives us parameter sensitivity indices for the outcomes. We average the parameter sensitivity indices across the outcomes to get the average SI for a scenario. For each realization that is generated from the Goldsim PA model, Goldsim generates random values for physical/empirical parameters that HELP uses as well. For each vector of physical/empirical parameters that Goldsim generates, the vector from the 5,000 HELP runs that is most 'similar' to the Goldsim vector is computed using the nearest neighbor approach. In this context 'similar' means minimization of the SA-weighted sum of the absolute differences among the 5,000 values computed for this statistic, where each value corresponds to a different HELP realization. In order to account for the impacts of climate change, this process was repeated using the spatially downscaled future climate projections. For each of the key parameters of interest, it was assumed that a linear change depicted the relationship between the values for the present day and those for 2100. In this way, the climatically-driven changes in key parameters used to inform the GoldSim model are quantified and incorporated into the PA model output for the future. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2025

Spacetime pq theory for AC and DC electric power systems

The 50/60 Hz alternating current (AC) electric power has been the standard and most flexible energy source powering our modern societies for one and a half centuries since the war of the currents: AC versus direct current (DC). A reactive power concept that was introduced at the beginning of the AC power was very useful for circuit/system analysis, design, control, optimization, and ultimately for more efficient and stable generation, transmission, distribution, and consumption. The initial reactive power theory was based on single-phase sinusoidal AC power to capture inductive and capacitive power that yields to net-zero average power over one fundamental cycle. Soon it was expanded to non-sinusoidal AC power and finally to instantaneous three-phase AC power. However, these reactive power theories remain separate and limited to special cases and have never been consolidated and made valid to all cases. Today, more widespread adoption of power electronics and renewable energy is bringing back DC power into the electric grids. The reactive power concept has never been applied to DC power systems. There is no reactive power in DC power systems according to the existing reactive power theories. Do DC power systems really have no reactive power? Capacitors and inductors are widely used in DC just like in AC power systems. Are they not reactive power components? Why are they different from their AC counterparts? Furthermore, are batteries active or reactive power components? What about active devices like power converters (or inverters) with AC (or DC) on one side and DC (or AC) on the other? Do they generate or consume reactive power? Finally, what about AC and DC hybrid power systems? How to define reactive power in such a complex power system that has a multitude of loads, buses, and sources? Is there reactive power between any two loads, any two buses, or any two sources in a power system and what is the total reactive power in such a complex power system as a whole? As the motivation and goal of this paper to answer the above basic questions, to unify the existing AC reactive power theories and to ultimately provide theoretical and insightful guidance for system analysis, design, control, efficiency, optimization, and operation of complex power systems, a concept of spacetime (both spatial and temporal) active and reactive power (pq) theory—the spatiotemporal aspect of active and reactive power—is developed for both AC and DC power systems. The theoretical definitions and physical meanings of the spacetime reactive power will be developed, and real applications and thought experiments/cases/exercises will be explored and discussed. The developed mathematics to define the active (or real) and reactive (or imaginary) power— p and q respectively by dot (scalar) and cross (vector) products of multi-dimension spacetime vectors and time-space mapping principle/law can have some fundamental implications as well.

24 POWER TRANSMISSION AND DISTRIBUTION↗