Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS TENSOR”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Interaction of radio frequency waves with cylindrical density filaments: scattering and radiation pressure

The propagation of radio-frequency (RF) waves in tokamaks can be affected by filamentary structures, or blobs, that are present in the edge plasma and the scrape-off layer. The difference in the permittivity between the surrounding plasma and interior of a filament leads to reflection, refraction and diffraction of the waves. This, in turn, can affect the power flow into the core of the plasma and reduce the efficiency of heating and/or current generation. The scattering of RF waves, lower hybrid, helicon and ion cyclotron waves, by a single cylindrical filament, embedded in a background plasma, is studied using a full-wave analytical theory developed previously (Ram & Hizanidis, Phys. Plasmas , vol. 23, 2016, 022504). The theory assumes that the plasma in and around a filament is homogeneous and cold. A detailed scattering analysis reveals a variety of common features that exist among the three distinctly different RF waves. These common attributes can be inferred intuitively based on an examination of the cold plasma dispersion relation. The physical intuition is a useful step to understanding experimental observations on scattering, as well as results from simulations that include general forms of edge plasma turbulence. While a filament can affect the propagation of RF waves, the radiation force exerted by the waves can influence the filament. The force on a filament is determined using the Maxwell stress tensor. In 1905, Poynting was the first to evaluate and measure the radiation force on an interface separating two different dielectric media (Poynting, London Edinburgh Dublin Philos. Mag. J. Sci. , vol. 9, 1905, pp. 393–406). For ordinary light propagating in vacuum and incident on a glass surface, Poynting noted that the surface is ‘pulled’ towards the vacuum. In a magnetized cold plasma, there are two independent wave modes. Even if only one of these modes is excited by an RF antenna, a filament will couple power to the other mode: a consequence of electromagnetic boundary conditions. This facet of scattering has consequences on the radiation force that go beyond Poynting's seminal contribution. The direction of the force depends on the polarization of the incident wave and on the mode structure of the waves inside and in the vicinity of a filament. It can either pull the filament toward the RF source or push it away. For slow lower hybrid waves, filaments with densities greater than the ambient density are pulled in, while filaments with lower densities are pushed out, thereby enhancing the density in front of the antenna. In the case of fast helicon and ion cyclotron waves, the direction of the force depends on the plasma and wave parameters; in particular, on the ambient density. The radiation force, in all three frequency ranges, is large enough to affect the motion of a filament and could be measured experimentally. This also suggests the possibility of modifying the edge turbulence using RF waves.

Physics↗
DOE OSTI2024

Maximum Entropy Theory of Multiscale Coarse-Graining via Matching Thermodynamic Forces: Application to a Molecular Crystal (TATB)

The MSCG/FM (multiscale coarse-graining via force-matching) approach is an efficient supervised machine learning method to develop microscopically informed coarse-grained (CG) models. Here we present a theory based on the principle of maximum entropy (PME) enveloping the existing MSCG/FM approaches. This theory views the MSCG/FM method as a special case of matching the thermodynamic forces from the extended ensemble described by the set of thermodynamic (relevant) system coordinates. This set may include CG coordinates, the stress tensor, applied external fields, and so forth, and may be characterized by nonequilibrium conditions. Following the presentation of the theory, we discuss the consistent matching of both bonded and nonbonded interactions. The proposed PME formulation is used as a starting point to extend the MSCG/FM method to the constant strain ensemble, which together with the explicit matching of the bonded forces is better suited for coarse-graining anisotropic media at a submolecular resolution. The theory is demonstrated by performing the fine coarse-graining of crystalline 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a well-known insensitive molecular energetic material, which exhibits highly anisotropic mechanical properties.

1,3,5-triamino-2,4,6-trinitrobenzene↗
DOE OSTI2020

Single Crystal Plasticity Model with Deformation Twinning for the High Rate Deformation of β-HMX

Deformation twinning is an active deformation mechanism during loading along certain crystallographic directions of β-cyclotetramethylene tetranitramine (β-HMX). In this work, a finite strain thermomechanical model is extended to include twinning as a deformation mechanism in addition to plastic slip. The stress is derived from the free energy expression including a term representing the equation of state. The crystal plasticity framework is used to divide the total strain into inelastic and elastic component, where the elastic part is used in the expression for the free energy. The twin systems are treated as pseudo slip systems and the shear rate on the twin systems is evaluated in terms of the twin resistance and the projection of the stress tensor. The model parameters were calibrated against a set of plate impact experiments performed on β-HMX by Dick et al. The remaining plate impact experiments are used to evaluate the predictive capability of the model. The quality of the model fits and predictions is discussed from a physical and modeling perspective. Particularly, the role of twin modeling on results is highlighted. Finally, the model is used to explore the relationship between the propensity for twinning and crystal orientation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Anti-symmetric representation of the extended magnetohydrodynamic equations

We introduce the anti-symmetric representation of the extended magnetohydrodynamic (MHD) equations. In this representation, the use of the anti-symmetric flux operator (∇·ν+ν·∇) results in conservation theorems with discrete analogs. Inherently robust numerical applications are achieved with little effort, and conservation to machine precision is possible with simple numerical schemes. Starting from the two-fluid equations, we construct a single-fluid MHD model based on generalized center-of-mass variables for the mass (√ρ), momentum (√ρν), and pressure (√ p ). This model is shown to possess identical conservation properties to the two-fluid system, with the only restriction being the use of a single temperature. Common approximations to the Braginskii heat fluxes and to the gyroviscous stress tensor are cast into our representation for convenience. The discrete conservation properties are verified using the classic Orszag–Tang vortex problem. In addition to the favorable mass, momentum, and energy conservation properties, the time reversibility of the simulations is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Calculation of the detonation state of HN 3 with quantum accuracy

HN 3 is a unique liquid energetic material that exhibits ultrafast detonation chemistry and a transition to metallic states during detonation. We combine the Chebyshev interaction model for efficient simulation (ChIMES) many-body reactive force field and the extended-Lagrangian multiscale shock technique molecular dynamics method to calculate the detonation properties of HN 3 with the accuracy of Kohn–Sham density-functional theory. ChIMES is based on a Chebyshev polynomial expansion and can accurately reproduce density-functional theory molecular dynamics (DFT-MD) simulations for a wide range of unreactive and decomposition conditions of liquid HN 3 . We show that addition of random displacement configurations and the energies of gas-phase equilibrium products in the training set allows ChIMES to efficiently explore the complex potential energy surface. Schemes for selecting force field parameters and the inclusion of stress tensor and energy data in the training set are examined. Structural and dynamical properties and chemistry predictions for the resulting models are benchmarked against DFT-MD. We demonstrate that the inclusion of explicit four-body energy terms is necessary to capture the potential energy surface across a wide range of conditions. Our results generally retain the accuracy of DFT-MD while yielding a high degree of computational efficiency, allowing simulations to approach orders of magnitude larger time and spatial scales. The techniques and recipes for MD model creation we present allow for direct simulation of nanosecond shock compression experiments and calculation of the detonation properties of materials with the accuracy of Kohn–Sham density-functional theory.

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

Electromagnetic torque measurements in DIII-D using internal/external magnetic field decomposition

Given spatially resolved measurements of normal and tangential components of the magnetic field just outside the surface of a magnetically confined plasma, the field at the measurement location can be uniquely decomposed into contributions from the plasma and from external sources. This principle allows direct measurement of the electromagnetic torque on the plasma without knowledge of the distribution of the internal and external currents, similar to the more well-known formalism using the Maxwell stress tensor. The internal/external field decomposition also enables a mixed approach that incorporates any explicitly known current distributions (e.g., from non-axisymmetric coils). We discuss the requirements and limitations of such an approach to torque measurements. Experimental measurements of the torque evolution as a rotating tearing mode locks to the wall in the DIII-D tokamak are consistent with a simple model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Apparatus for measuring strength in biaxial compression

We report most measurements of compressive strength of ductile materials have involved Hopkinson–Kolsky bars or Taylor anvils placing samples in uniaxial compression. In these geometries, strain is limited by the tendency of the sample to petal, in analogy to necking in uniaxial tension. Estimation of strength for any other form of the stress tensor requires assuming a shape of the yield surface; because data exist only for uniaxial compression, these assumptions are untested. In an imploding spherical shell, compression is biaxial, the plastic strain may not be small, and the material behavior may be nonlinear as a result of work hardening and heating by plastic work. We outline a method of measuring the strengths of materials in biaxial compression, both quasistatically and dynamically, using the compression of thin spherical shells. We suggest surrounding the shell with an annulus filled with a mixture of H 2 and Cl 2 gases whose homogeneous ignition is initiated by a flash of blue and near-ultraviolet light. Less promising approaches are described in Appendixes A–C .

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Real-space density kernel method for Kohn–Sham density functional theory calculations at high temperature

Kohn–Sham density functional theory calculations using conventional diagonalization based methods become increasingly expensive as temperature increases due to the need to compute increasing numbers of partially occupied states. In this work, we present a density matrix based method for Kohn–Sham calculations at high temperatures that eliminates the need for diagonalization entirely, thus reducing the cost of such calculations significantly. Specifically, we develop real-space expressions for the electron density, electronic free energy, Hellmann–Feynman forces, and Hellmann–Feynman stress tensor in terms of an orthonormal auxiliary orbital basis and its density kernel transform, the density kernel being the matrix representation of the density operator in the auxiliary basis. Using Chebyshev filtering to generate the auxiliary basis, we next develop an approach akin to Clenshaw–Curtis spectral quadrature to calculate the individual columns of the density kernel based on the Fermi operator expansion in Chebyshev polynomials and employ a similar approach to evaluate band structure and entropic energy components. We implement the proposed formulation in the SPARC electronic structure code, using which we show systematic convergence of the aforementioned quantities to exact diagonalization results, and obtain significant speedups relative to conventional diagonalization based methods. Finally, we employ the new method to compute the self-diffusion coefficient and viscosity of aluminum at 116 045 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Accelerating boundary analog of a Kerr black hole

An accelerated boundary correspondence (i.e. a flat spacetime accelerating mirror trajectory) is derived for the Kerr spacetime, with a general formula that ranges from the Schwarzschild limit (zero angular momentum) to the extreme maximal spin case (yielding asymptotic uniform acceleration). Here, the beta Bogoliubov coefficients reveal the particle spectrum is a Planck distribution at late times with temperature cooler than a Schwarzschild black hole, due to the 'spring constant' analog of angular momentum. The quantum stress tensor indicates a constant emission of energy flux at late times consistent with eternal thermal equilibrium.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Local well-posedness and singularity formation in non-Newtonian compressible fluids

We investigate the initial value problem of a very general class of 3 + 1 non-Newtonian compressible fluids in which the viscous stress tensor with shear and bulk viscosity relaxes to its Navier–Stokes values. These fluids correspond to the non-relativistic limit of well-known Israel–Stewart-like theories used in the relativistic fluid dynamic simulations of high-energy nuclear and astrophysical systems. After establishing the local well-posedness of the Cauchy problem, we show for the first time in the literature that there exists a large class of initial data for which the corresponding evolution breaks down in finite time due to the formation of singularities. This implies that a large class of non-Newtonian fluids do not have finite solutions defined at all times.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Quantum critical collapse abhors a naked singularity

Classical critical collapse yields naked singularities from smooth initial data, challenging cosmic censorship, and shaping the spectrum of primordial black holes. We show that one-loop vacuum polarization near the threshold qualitatively changes this outcome by dressing the singularity with a horizon within a controlled semiclassical regime. In analytically tractable Einstein-scalar critical spacetimes, a one-loop $s$-wave treatment linearized around self-similar backgrounds shows that regularity uniquely selects an asymptotically Minkowskian, vacuum-polarization state. Its renormalized stress tensor carries a universal quantum growing mode that competes with the classical unstable mode, shifts the critical point, and generates a trapped surface along with a finite mass gap at the new threshold, thereby enforcing horizon formation even under arbitrary fine-tuning. In primordial collapse, the threshold shift enters exponentially into the formation fraction, while the mass gap truncates the low-mass tail, suggesting potentially important consequences for the predicted mass spectrum. Furthermore, these results provide a self-consistent semiclassical treatment of critical collapse and yield sharp predictions within the one-loop, near-critical, linearized regime.

Anomalies↗
DOE OSTI2026

Quantum properties of non-Dirichlet boundary conditions in gravity

The Euclidean path integral for gravity is enriched by the addition of boundaries, which provide useful probes of thermodynamic properties. Common boundary conditions include Dirichlet conditions on the boundary induced metric; microcanonical conditions, which refers to fixing some components of the Brown-York boundary stress tensor; and conformal conditions, in which the conformal structure of the induced metric and the trace of the extrinsic curvature are fixed. Boundaries also present interesting problems of consistency. The Dirichlet problem is known, under various (and generally different) conditions, to be inconsistent with perturbative quantization of graviton fluctuations, to exhibit thermodynamic instability, or to require infinite fine-tuning in the presence of matter fluctuations. We extend some of these results to other boundary conditions. We find that similarly to the Dirichlet problem, the graviton fluctuation operator is not elliptic with microcanonical boundaries, and the nonelliptic modes correspond to “boundary-moving diffeomorphisms.” However, we argue that microcanonical factorization of path integrals—essentially, the insertion of microcanonical constraints on two-sided surfaces in the bulk—is not affected by the same issues of ellipticity. We also show that for a variety of matter field boundary conditions, matter fluctuations renormalize the gravitational bulk and boundary terms differently, so that the classical microcanonical or conformal variational problems are not preserved unless an infinite fine-tuning is performed.

Draper, Patrick [Univ. of Illinois at Urbana-Champ↗
DOE OSTI2021

Structure and energetics of the elbows in the Au(111) herringbone reconstruction

We study the structure of the threading edge dislocations, or “elbows,” which are an essential component of the well-known herringbone reconstruction of the (111) surface of Au. Previous work had shown that these dislocations can be stabilized by long-range elastic relaxations into the bulk. However, the validity of the harmonic spring model that had been used to estimate the energies of the dislocations is uncertain. To enable a more refined model of the dislocation energetics, we have imaged the atomic structure of these dislocations using scanning tunneling microscopy. We find that the harmonic spring model does not adequately reproduce the observed structure. We are able to reproduce the structure, however, with a two-dimensional Frenkel-Kontorova (FK) model that uses a pairwise Morse potential to describe the interactions between the top layer Au atoms on a rigid substrate. The parameters of the potential were obtained by fitting the energy of uniaxially compressed phases, or “stripes”, computed with density functional theory, as a function of surface Au density. Within this model, the formation of the threading dislocations remains unfavorable. However, the large forces on the substrate atoms near the threading-dislocation cores, render the assumption of a completely rigid substrate questionable. Indeed, if the FK parameters are modified to account for the relaxation of just one more atomic layer, threading dislocations can, in principle, become favorable, even without bulk elastic relaxations. Additional evidence for a small elbow energy is that our computed change in the Au(111) surface stress tensor caused by the (√ 3 × 22) reconstruction is considerably smaller than previous estimates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Impact of fragment formation on shear viscosity in the nuclear liquid-gas phase transition region

Within the improved quantum molecular dynamics (ImQMD) model we follow the evolution of nuclear matter for planar Couette flow in a periodic box. We focus on the region of liquid-gas phase transition and extract the shear viscosity coefficient from the local stress tensor, directly following the viscosity definition. Here, by switching on and off the mean field and thus inducing the phase transition, we are able to observe the impact of clumping in the phase transition region on the viscosity.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Mirror at the edge of the universe: Reflections on an accelerated boundary correspondence with de Sitter cosmology

An accelerated boundary correspondence (ABC) is solved for the de Sitter moving mirror cosmology. The beta Bogoliubov coefficients reveal the particle spectrum is a Planck distribution with temperature inversely proportional to horizon radius. The quantum stress tensor indicates a constant emission of energy flux consistent with eternal equilibrium, while the total energy carried by the particles remains finite. In conclusion, the curved spacetime transformation to flat spacetime with an accelerated boundary is illustrated, and also shown for anti-de Sitter (AdS) spacetime.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Far-from-equilibrium search for the QCD critical point

Initial conditions for relativistic heavy-ion collisions may be far from equilibrium (i.e., there are large initial contributions from the shear stress tensor and bulk pressure) but it is expected that on very short time scales the dynamics converge to a universal attractor that defines hydrodynamic behavior. Thus far, studies of this nature have only considered an idealized situation at LHC energies (high temperatures T and vanishing baryon chemical potential µ B = 0) but, in this work, we investigate for the first time how far-from-equilibrium effects may influence experimentally driven searches for the Quantum Chromodynamic critical point at RHIC. We find that the path to the critical point is heavily influenced by far from equilibrium initial conditions where viscous effects lead to dramatically different {T, µ B } trajectories through the QCD phase diagram. We compare hydrodynamic equations of motion with shear and bulk coupled together at finite µB for both DNMR and phenomenological Israel-Stewart equations of motion and discuss their influence on potential attractors at finite µ B and their corresponding {T, µ B } trajectories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Soft pions and transport near the chiral critical point

During the expansion of a heavy ion collision, the system passes close to the $\textit{O}$(4) critical point of QCD, and thus the fluctuations of the order parameter $(σ, \overrightarrow{π}$) are expected to be enhanced. In this work, our goal is to compute how these enhanced fluctuations modify the transport coefficients of QCD near the pseudocritical point. We also make a phenomenological estimate for how chiral fluctuations could effect the momentum spectrum of soft pions. We first formulate the appropriate stochastic hydrodynamic equations close to the $\textit{O}$(4) critical point. Then, working in mean field, we determine the correlation functions of the stress tensor and the currents which result from this stochastic real-time theory, and use these correlation functions to determine the scaling behavior of the transport coefficients. The hydrodynamic theory also describes the propagation of pion waves, fixing the scaling behavior of the dispersion curve of soft pions. We present scaling functions for the shear viscosity and the charge conductivities near the pseudocritical point and estimate the absolute magnitude of the critical fluctuations to these parameters and the bulk viscosity. Using the calculated pion dispersion curve, we estimate the expected critical enhancement of soft pion yields, and this estimate provides a plausible explanation for the excess seen in experiment relative to ordinary hydrodynamic computations. Our results motivate further phenomenological and numerical work on the implications of chiral symmetry on real-time properties of thermal QCD near the pseudocritical point.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Slowly rotating gravastars

We solve Einstein’s equations for slowly rotating gravitational condensate stars (gravastars) up to second order in the rotation by expanding about the spherically symmetric gravastar with de Sitter interior and Schwarzschild exterior matched at their common horizon. Requiring that the perturbations are finite on the null surface reduces the exterior geometry to that of a Kerr black hole, implying that a slowly rotating gravastar cannot be distinguished from a Kerr black hole by any measurement or observation restricted to the macroscopic spacetime exterior to the horizon. Here we determine the interior solution, the surface stress tensor, and the Komar mass and angular momentum localized on the slowly rotating horizon surface. With the interior equation of state fixed at p = –ρ, finite junction conditions on the null horizon surface necessarily lead to an interior solution with a singular core, where the perturbative expansion breaks down. Comparison to other models and implications for more rapidly rotating gravastars are briefly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗