SEARCH · Engineering Papers

Results for “Physical quantities”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Drude weights in one-dimensional systems with a single defect

Ballistic transport of a quantum system can be characterized by Drude weight, which quantifies the response of the system to a uniform electric field in the infinitely long timescale. The Drude weight is often discussed in terms of the Kohn formula, which gives the Drude weight by the derivative of the energy eigenvalue of a finite-size system with the periodic boundary condition in terms of the Aharonov-Bohm flux. Recently, the Kohn formula is generalized to nonlinear responses. However, the nonlinear Drude weight determined by the Kohn formula often diverges in the thermodynamic limit. In order to elucidate the issue, in this work we examine a simple example of a one-dimensional tight-binding model in the presence of a single defect at zero temperature. We find that its linear and nonlinear Drude weights given by the Kohn formula (i) depend on the Aharonov-Bohm flux and (ii) diverge proportionally to a power of the system size. Here, we argue that the problem can be attributed to different order of limits. The Drude weight according to the Kohn formula (“Kohn-Drude weight”) indicates the response of a finite-size system to an adiabatic insertion of the Aharonov-Bohm flux. While it is a well-defined physical quantity for a finite-size system, its thermodynamic limit does not always describe the ballistic transport of the bulk. The latter should be rather characterized by a “bulk Drude weight” defined by taking the thermodynamic limit first before the zero-frequency limit. While the potential issue of the order of limits has been sometimes discussed within the linear response, the discrepancy between the two limits is amplified in nonlinear Drude weights. We demonstrate the importance of the low-energy excitations of O(1/L), which are excluded from the Kohn-Drude weight, in regularizing the bulk Drude weight.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2023

Transport across interfaces in symmetric orbifolds

We examine how conformal boundaries encode energy transport coefficients — namely transmission and reflection probabilities — of corresponding conformal interfaces in symmetric orbifold theories. These constitute a large class of irrational theories and are closely related to holographic setups. Our central goal is to compare such coefficients at the orbifold point (a field theory calculation) against their values when the orbifold is highly deformed (a gravity calculation) — an approach akin to past AdS/CFT-guided comparisons of physical quantities at strong versus weak coupling. At the orbifold point, we find that the (weighted-average) transport coefficients are simply averages of coefficients in the underlying seed theory. We then focus on the symmetric orbifold of the T 4 sigma model interface CFT dual to type IIB supergravity on the 3d Janus solution. We compare the holographic transmission coefficient, which was found by [1], to that of the orbifold point. We find that the profile of the transmission coefficient substantially increases with the coupling, in contrast to boundary entropy. We also present some related ideas about twisted-sector data encoded by boundary states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Probabilistic neural networks for improved analyses with phenomenological R -matrix

Here we present a method for measurement analyses based on probabilistic deep neural networks that provide several advantages over conventional analyses with phenomenological models. These include predicting physical quantities directly from data, the rapid generation of statistically robust uncertainties, and the ability to bypass some parameters that may induce ambiguities and complications in data analysis. As deep learning methods make predictions through “black boxes,” the uncertainty quantification is typically challenging. We use a probabilistic framework that provides thorough uncertainty quantification and is straightforward to follow in practice. With the network architecture based on the Transformer, we demonstrate the current method for predicting nuclear resonance parameters from scattering data using the phenomenological R-matrix model.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Probing Postmeasurement Entanglement without Postselection

We study the problem of observing quantum collective phenomena emerging from large numbers of measurements. These phenomena are difficult to observe in conventional experiments because, in order to distinguish the effects of measurement from dephasing, it is necessary to postselect on sets of measurement outcomes with Born probabilities that are exponentially small in the number of measurements performed. An unconventional approach, which avoids this exponential “postselection problem”, is to construct cross-correlations between experimental data and the results of simulations on classical computers. However, these cross-correlations generally have no definite relation to physical quantities. We first show how to incorporate classical shadows into this framework, thereby allowing for the construction of quantum information-theoretic cross-correlations. We then identify cross-correlations that both upper and lower bound the measurement-averaged von Neumann entanglement entropy, as well as cross-correlations that lower bound the measurement-averaged purity and entanglement negativity. These bounds show that experiments can be performed to constrain postmeasurement entanglement without the need for postselection. To illustrate our technique, we consider how it could be used to observe the measurement-induced entanglement transition in Haar-random quantum circuits. We use exact numerical calculations as proxies for quantum simulations and, to highlight the fundamental limitations of classical memory, we construct cross-correlations with tensor-network calculations at finite bond dimension. Our results reveal a signature of measurement-induced criticality that can be observed using a quantum simulator in polynomial time and with polynomial classical memory. Published by the American Physical Society 2024

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1999

Parallel Grid Manipulations in Earth Science Calculations

The National Aeronautics and Space Administration (NASA) Data Assimilation Office (DAO) at the Goddard Space Flight Center is moving its data assimilation system to massively parallel computing platforms. This parallel implementation of GEOS DAS will be used in the DAO's normal activities, which include reanalysis of data, and operational support for flight missions. Key components of GEOS DAS, including the gridpoint-based general circulation model and a data analysis system, are currently being parallelized. The parallelization of GEOS DAS is also one of the HPCC Grand Challenge Projects. The GEOS-DAS software employs several distinct grids. Some examples are: an observation grid- an unstructured grid of points at which observed or measured physical quantities from instruments or satellites are associated- a highly-structured latitude-longitude grid of points spanning the earth at given latitude-longitude coordinates at which prognostic quantities are determined, and a computational lat-lon grid in which the pole has been moved to a different location to avoid computational instabilities. Each of these grids has a different structure and number of constituent points. In spite of that, there are numerous interactions between the grids, e.g., values on one grid must be interpolated to another, or, in other cases, grids need to be redistributed on the underlying parallel platform. The DAO has designed a parallel integrated library for grid manipulations (PILGRIM) to support the needed grid interactions with maximum efficiency. It offers a flexible interface to generate new grids, define transformations between grids and apply them. Basic communication is currently MPI, however the interfaces defined here could conceivably be implemented with other message-passing libraries, e.g., Cray SHMEM, or with shared-memory constructs. The library is written in Fortran 90. First performance results indicate that even difficult problems, such as above-mentioned pole rotation- a sparse interpolation with little data locality between the physical lat-lon grid and a pole rotated computational grid- can be solved efficiently and at the GFlop/s rates needed to solve tomorrow's high resolution earth science models. In the subsequent presentation we will discuss the design and implementation of PILGRIM as well as a number of the problems it is required to solve. Some conclusions will be drawn about the potential performance of the overall earth science models on the supercomputer platforms foreseen for these problems.

Sawyer, W.↗

DOE OSTI2021

Design and engineering challenges of a multi-energy hard x-ray camera for long-pulse profile measurements at WEST tokamak

The WEST tokamak has recently been prepared for long pulse operation with a water-cooled full first wall. Heating is provided by radiofrequency systems, including Lower Hybrid Current Drive (LHCD). The Princeton Plasma Physics Laboratory has developed a multi-energy hard x-ray camera for profile measurements at WEST. The camera is based on a new generation of 2D pixel array detectors that allow the energy threshold to be independently set for each pixel. The diagnostic will provide spatial, temporal and energy resolved measurements of the hard x-ray emission from the full plasma cross-section, investigating several physical quantities such as the electron temperature from continuum emission, the fast electron tail density produced by radiofrequency current drive and runaway electrons, as well as characteristic tungsten x-rays due to beam-target emission at the edge. Overall, this work describes the engineering challenges that the WEST long pulse scenario poses for this diagnostic and how the design addresses and solves them. Vacuum, thermal stress and heat transfer calculations are presented and discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Stanford Thermal Earth Model for the Conterminous United States

Provided here are various forms of the Stanford Thermal Earth Model, as well as the data and methods used for its creation. The predictions produced by this model were visualized in two-dimensional spatial maps across the modeled depths (0-7 km) for the conterminous United States. The thermal earth model is made available as an application programming interface (API) and as feature layers on ArcGIS, which are both provided via links below. A data-driven spatial interpolation algorithm based on physics-informed graph neural networks was used to develop these national temperature-at-depth maps. The model satisfied the three-dimensional heat conduction law by predicting subsurface temperature, surface heat flow, and rock thermal conductivity. Many physical quantities, including bottomhole temperature, depth, geographic coordinates, elevation, sediment thickness, magnetic anomaly, gravity anomaly, gamma-ray flux of radioactive elements, seismicity, and electric conductivity were used as model inputs. Surface heat flow, temperature, and thermal conductivity predictions were constructed for depths of 0-7 km at an interval of 1 km with spatial resolution of 18 km2 per grid cell. The model showed superior temperature, surface heat flow and thermal conductivity mean absolute errors of 4.8C, 8.1 mW/m2 and 0.07 W/(C-m), respectively..

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Numerical minimization of neoclassical poloidal viscosity for supersonic equilibria in tokamak geometry

Extensive experimental evidence has shown that the presence of poloidal flow in tokamaks can dramatically improve transport properties. However, theory indicates that poloidal flows are damped by poloidal viscosity, thus necessitating external drivers, such as neutral beam injection or radio frequency heating. In this work, ideal magnetohydrodynamic equilibria are calculated via the FORTRAN code FLOW [Guazzotto et al., Phys. Plasmas 11, 604 (2004)] and a postprocessor is used to estimate the neoclassical poloidal viscosity. The equilibrium inputs, which correspond to intuitive physical quantities, are then numerically optimized to reduce a viscosity figure of merit. We present supersonic equilibria in tokamak geometry with minimized neoclassical poloidal viscosities for various velocity free function inputs, plasma aspect ratios, and collisionality regimes. Benchmarks are made against an analytic theory as well as a classical expression of poloidal viscosity. Numerical confirmation of the analytic theory is obtained in the high aspect ratio and high collisionality limit. Good agreement is also seen near the plasma core and edge, with discrepancies arising in the intermediate region. Outside of these limits, rotation input function profiles are found that provide ∼order of magnitude improvements over the analytic theory, with additional progress being made toward predictions for tokamak-relevant equilibria.

Physics↗

DOE OSTI2024

Quantum Scales of Galaxies from Self-interacting Ultralight Dark Matter

We derive the characteristic scales for physical quantities of dwarf galaxies, such as mass, size, acceleration, and angular momentum, within the self-interacting ultralight dark matter (ULDM) model. Due to the small mass of ULDM, even minor self-interactions can drastically alter these scales in the Thomas-Fermi limit. We suggest that these characteristic scales are connected to mysteries of observed galaxies. Oscillation of ULDM field can explain the current cosmological density of dark matter. Many cosmological constraints suggest that the energy scale $\tilde{m}$ for self-interacting ULDM is typically of the order $10~eV$, whereas the mass $m$ for the non-interacting case is around $10^{-21}~eV$. Self-interacting ULDM provides the better explanation for cosmological observations than the non-interacting case.

Astrophysics of Galaxies (astro-ph.GA)↗

DOE OSTI2025

Do We Know How to Model Reionization?

I compare the power spectra of the radiation fields from two recent sets of fully-coupled simulations that model cosmic reionization: “Cosmic Reionization On Computers” (CROC) and “Thesan”. While both simulations have similar power spectra of the radiation sources, the power spectra of the photoionization rate are significantly different at the same values of cosmic time or the same values of the mean neutral hydrogen fraction. However, the power spectra of the photoionization rate can be matched at large scales for the two simulations when the matching snapshots are allowed to vary independently. I.e., on large scales, the clustering of the radiation field in two simulations evolves similarly, but the exact timing of this evolution is different in different simulations and is not parameterized by an easily interpretable physical quantity like the mean neutral fraction or the mean free path. On small scales, large differences are present and remain partially unexplained. Both CROC and Thesan use the Variable Eddington Tensor approximation for modeling radiative transfer, but adopt different closure relations (optically thin OTVET versus M1). The role of this key difference is tested by using smaller simulations with a new cosmological simulation code that implements both closure relations in a controlled environment (the same hydro, cooling, and gravity solvers and the star formation recipe). In these controlled tests, both the M1 closure and the OTVET ansatz follow the expected behavior from a simple analytical approximation, demonstrating that the differences in the 2-point function of the radiation field induced by the choice of the Eddington tensor are not dominant.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

HPC-driven computational reproducibility in numerical relativity codes: a use case study with IllinoisGRMHD

Abstract Reproducibility of results is a cornerstone of the scientific method. Scientific computing encounters two challenges when aiming for this goal. Firstly, reproducibility should not depend on details of the runtime environment, such as the compiler version or computing environment, so results are verifiable by third-parties. Secondly, different versions of software code executed in the same runtime environment should produceconsistent numerical results for physical quantities. In this manuscript, we test the feasibility of reproducing scientific results obtained using theIllinoisGRMHDcode that is part of an open-source community software for simulation in relativistic astrophysics, theEinstein Toolkit. We verify that numerical results of simulating a single isolated neutron star withIllinoisGRMHDcan be reproduced, and compare them to results reported by the code authors in 2015. We use two different supercomputers: Expanse at SDSC, and Stampede2 at TACC. By compiling the source code archived along with the paper on both Expanse and Stampede2, we find thatIllinoisGRMHDreproduces results published in its announcement paper up to errors comparable to round-off level changes in initial data parameters. We also verify that a current version ofIllinoisGRMHDreproduces these results once we account for bug fixes which have occurred since the original publication.

Astronomy & Astrophysics↗

DOE OSTI2020

Frustration and Atomic Ordering in a Monolayer Semiconductor Alloy

Frustrated interactions can lead to short-range ordering arising from incompatible interactions of fundamental physical quantities with the underlying lattice. The simplest example is the triangular lattice of spins with antiferromagnetic interactions, where the nearest-neighbor spin-spin interactions cannot simultaneously be energy minimized. In this work, we show that engineering frustrated interactions is a possible route for controlling structural and electronic phenomena in semiconductor alloys. Using aberration-corrected scanning transmission electron microscopy in conjunction with density functional theory calculations, we demonstrate atomic ordering in a two-dimensional semiconductor alloy as a result of the competition between geometrical constraints and nearest-neighbor interactions. Statistical analyses uncover the presence of short-range ordering in the lattice. In addition, we show how the induced ordering can be used as another degree of freedom to considerably modify the band gap of monolayer semiconductor alloys.

2-dimensional systems↗

DOE OSTI2024

The inviscid incompressible limit of Kelvin–Helmholtz instability for plasmas

The Kelvin–Helmholtz Instability (KHI) is an interface instability that develops between two fluids or plasmas flowing with a common shear layer. KHI occurs in astrophysical jets, solar atmosphere, solar flows, cometary tails, planetary magnetospheres. Two applications of interest, encompassing both space and fusion applications, drive this study: KHI formation at the outer flanks of the Earth’s magnetosphere and KHI growth from non-uniform laser heating in magnetized direct-drive implosion experiments. Here, we study 2D KHI with or without a magnetic field parallel to the flow. We use both the GAMERA code, which solves the compressible Euler equations, and the STRATOSPEC code, which solves the Navier-Stokes equations under the Boussinesq approximation, coupled with the magnetic field dynamics. GAMERA is a global three-dimensional MHD code with high-order reconstruction in arbitrary nonorthogonal curvilinear coordinates, which is developed for a large range of astrophysical applications. STRATOSPEC is a three-dimensional pseudo-spectral code with an accuracy of infinite order (no numerical diffusion). Magnetized KHI is a canonical case for benchmarking hydrocode simulations with extended MHD options. An objective is to assess whether or not, and under which conditions, the incompressibility hypothesis allows to describe a dynamic compressible system. For comparing both codes, we reach the inviscid incompressible regime, by decreasing the Mach number in GAMERA, and viscosity and diffusion in STRATOSPEC. Here, we specifically investigate both single-mode and multi-mode initial perturbations, either with or without magnetic field parallel to the flow. The method relies on comparisons of the density fields, 1D profiles of physical quantities averaged along the flow direction, and scale-by-scale spectral densities. We also address the triggering, formation and damping of filamentary structures under varying Mach number or Atwood number, with or without a parallel magnetic field. Comparisons show very satisfactory results between the two codes. The vortices dynamics is well reproduced, along with the breaking or damping of small-scale structures. We end with the extraction of growth rates of magnetized KHI from the compressible regime to the incompressible limit in the linear regime assessing the effects of compressibility under increasing magnetic field. The observed differences between the two codes are explained either from diffusion or non-Boussinesq effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2021

Probing phonon softening in ferroelectrics by scanning probe microwave spectroscopy

Microwave measurements have recently been successfully applied to measure ferroelectric materials on the nanoscale, including detection of polarization switching and ferroelectric domain walls. In this study we discuss the question of whether scanning probe microscopy operating at microwave frequency can identify the changes associated with the soft phonon dynamics in a ferroic. The analytical expressions for the electric potential, complex impedance, and dielectric losses are derived and analyzed, since these physical quantities are linked to experimentally measurable properties of the ferroic. As a ferroic we consider virtual or proper ferroelectric with an optic phonon mode that softens at a Curie point. We also consider a decay mechanism linked to the conductance of the ferroic, thus manifesting itself as the dielectric loss in the material. Our key finding is that the influence of the soft phonon dispersion on the surface potential distribution, complex impedance, and dielectric losses are evidently strong in the vicinity (~10–30 K) of the Curie temperature. Additionally, we quantified how the spatial distribution and frequency spectra of the complex impedance and the dielectric losses react on the dynamics of the soft phonons near the Curie point. These results set the stage for characterization of polar phase transitions with nanoscale microwave measurements, providing a complementary approach to well established electromechanical measurements for fundamental understanding of ferroelectric properties as well as their applications in telecommunication and computing.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Coupling of CTF and TRACE for Modeling of Transients

This report documents the improvements that have been made to the capabilities for coupling CTF to systems codes-specifically, the US Nuclear Regulatory Commission (NRC) TRACE code. An initial systems coupling capability had been set up previously using a nonoverlapping domain approach with the codes exchanging data at the core boundaries. The present work adds a new approach using overlapping domains, in which the system code models the core as well. A new input format has been added to allow the user to specify the physical quantities to be exchanged and their location in the system model, which gives the flexibility of applying one-way or two-way coupling between the codes using the desired data exchanges. In addition to applying thermal-hydraulics (T/H) boundary condition (BC) values obtained from TRACE, a capability was added to allow CTF to apply flow resistance feedback to TRACE to match the CTF core pressure drop. Support was added for executing parallel CTF models within the CTF-systems coupling. The system coupling capability was successfully applied to a parallel MSLB transient, demonstrating that both the one-way and two-way coupling behaved as expected and provided substantial improvements to numerical stability and routine compared to the previous nonoverlapping domain coupling. An initial capability was also developed for performing restart calculations in CTF which will be used in the future for restarting CTF-TRACE simulations at specific points in the transient simulation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Stellarator optimization for good magnetic surfaces at the same time as quasisymmetry

A method is demonstrated to optimize a stellarator's geometry to eliminate magnetic islands and achieve other desired physics properties at the same time. For many physics quantities that have been used in stellarator optimization, including quasisymmetry, neoclassical transport, and magnetohydrodynamic stability, it is convenient to use a magnetic equilibrium representation that assures the existence of magnetic surfaces. However, this representation hides the possible presence of magnetic islands, which are typically undesirable. To include both surface-based objectives and island widths in a single optimization, two fixed-boundary equilibrium calculations are run at each iteration of the optimization: one that enforces the existence of magnetic surfaces (the Variational Moments Equilibrium Code) [S. P. Hirshman and J. C. Whitson, Phys. Fluids 26, 3553 (1983)] and one that does not (the Stepped Pressure Equilibrium Code) [Hudson et al., Phys. Plasmas 19, 112502 (2012)]. By penalizing the island residues in the objective function, the two magnetic field representations are brought into agreement during the optimization. Here, an example is presented in which, particularly on the surface where quasisymmetry was targeted, quasisymmetry is achieved more accurately than in previously published examples.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Duality defect in a deformed transverse-field Ising model

Physical quantities with long lifetimes have both theoretical significance in the study of quantum many-body systems and practical implications for quantum technologies. In this manuscript, we investigate the roles played by topological defects in the construction of quasiconserved quantities, using as a prototypical example the Kramers-Wannier duality defect in a deformed one-dimensional quantum transverse-field Ising model. We construct the duality defect Hamiltonian in three different ways: half-chain Kramers-Wannier transformation, utilization of techniques in the Ising fusion category, and defect-modified weak integrability breaking deformation. The third method is also applicable for the study of generic integrable defects under weak integrability breaking deformations. We also work out the deformation of defect-modified higher charges in the model and study their slower decay behavior. Furthermore, we consider the corresponding duality defect twisted deformed Floquet transverse-field Ising model and investigate the stability of the isolated zero mode associated with the duality defect in the integrable Floquet Ising model, under such weak integrability breaking deformation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Uncertainty quantification for Multiphase-CFD simulations of bubbly flows: a machine learning-based Bayesian approach supported by high-resolution experiments

In this paper, we developed a machine learning-based Bayesian approach to inversely quantify and reduce the uncertainties of multiphase computational fluid dynamics (MCFD) simulations for bubbly flows. The proposed approach is supported by high-resolution two-phase flow measurements, including those by double-sensor conductivity probes, high-speed imaging, and particle image velocimetry. Local distributions of key physical quantities of interest (QoIs), including the void fraction and phasic velocities, are obtained to support the Bayesian inference. In the process, the epistemic uncertainties of the closure relations are inversely quantified while the aleatory uncertainties from stochastic fluctuations of the system are evaluated based on experimental uncertainty analysis. The combined uncertainties are then propagated through the MCFD solver to obtain uncertainties of the QoIs, based on which probability-boxes are constructed for validation. The proposed approach relies on three machine learning methods: feedforward neural networks and principal component analysis for surrogate modeling, and Gaussian processes for model form uncertainty modeling. The whole process is implemented within the framework of an open-source deep learning library PyTorch with graphics processing unit (GPU) acceleration, thus ensuring the efficiency of the computation. The results demonstrate that with the support of high-resolution data, the uncertainties of MCFD simulations can be significantly reduced. The proposed approach has the potential for other applications that involve numerical models with empirical parameters.

42 ENGINEERING↗