Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TENSOR FORCES”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Short-range components of nuclear forces: Experiment versus mythology

The present-day situation around the description of various (central, spin-orbit, and tensor) components of short-range nuclear forces is discussed. A traditional picture of these interactions based on the idea of one-meson exchange is contrasted against numerous results of recent experiments. As is shown in the present study, these results often deviate strongly from the predictions of traditional models. One can therefore state that such models are inapplicable to describing short-range nuclear forces and that it is necessary to go over from a traditional description to some alternative QCD-based (or QCD-motivated) picture. This means that, despite the widespread popularity of traditional concepts of short-range nuclear forces and their applicability in many particular cases, these concepts are not more than scientific myths that show their inconsistency when analyzed from the viewpoint of the modern experiment.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Journey to Time-Variable Moment Tensors through Inversion of Acoustic and Seismoacoustic Data

We explore the capability of acoustic and seismoacoustic datasets to directly resolve a complex, time-variable source consisting of a buried mechanism, represented as a moment tensor, and a spall mechanism, represented as a vertical force at the surface. Traditionally, each component of a resolved moment tensor assumes one underlying source time function, which likely fails to capture the full evolution of a dynamic source, such as an explosion followed by slip on near-source joints or development of spallation. Specifically, we expand previous work to resolve a time-variable moment tensor using single-modality and joint-modality inversion frameworks through analysis of infrasound and seismoacoustic data recorded as part of the Source Physics Experiment Phase II: Dry Alluvium Geology (DAG). We investigate the impact of including signals from seismic-to-air coupling that are local to each infrasound sensor in comparison to mainly atmosphere-propagating acoustic signals, which occur from coupling of the wavefield from the subsurface to the atmosphere directly above the source. Additionally, we assess the ability of our inversion algorithm to fit observed infrasound data using a variety of time-variable source mechanisms. First, we consider the buried moment tensor source alone, which assumes that the determined Green’s functions incorporate effects from spallation or that the impact from spallation is minimal. Second, we examine the estimated buried moment tensor and vertical surface spallation as terms that must both be resolved in the inversion. Third, we assess the ability for an estimated vertical surface spallation source to fit the acoustic data on its own. Finally, we compare results from the joint inversion of both seismic geophone and infrasound acoustic data for the buried-only source compared to buried and spallation sources. Our results are a preliminary investigation into the applications of the inversion technique to recorded datasets and show the technique has limited capabilities using acoustic data alone. Instead, this method shows promise for seismic and seismoacoustic datasets to resolve the time-variable mechanisms of a buried source.

47 OTHER INSTRUMENTATION↗
DOE OSTI2018

Efficient tree decomposition of high-rank tensors [Efficient decomposition of high-rank tensors]

Tensors are a natural way to express correlations among many physical variables, but storing tensors in a computer naively requires memory which scales exponentially in the rank of the tensor. This is not optimal, as the required memory is actually set not by the rank but by the mutual information amongst the variables in question. Representations such as the tensor tree perform near-optimally when the tree decomposition is chosen to reflect the correlation structure in question, but making such a choice is non-trivial and good heuristics remain highly context-specific. In this work I present two new algorithms for choosing efficient tree decompositions, independent of the physical context of the tensor. The first is a brute-force algorithm which in most cases produces optimal decompositions but is generally impractical for high-rank tensors, as the number of possible choices grows exponentially in rank. Here, the second is a greedy algorithm, and while it is not optimal it performs extremely well in numerical experiments while having runtime which makes it practical even for tensors of very high rank.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2015

An EQT-cDFT approach to determine thermodynamic properties of confined fluids

We present a continuum-based approach to predict the structure and thermodynamic properties of confined fluids at multiple length-scales, ranging from a few angstroms to macro-meters. The continuum approach is based on the empirical potential-based quasi-continuum theory (EQT) and classical density functional theory (cDFT). EQT is a simple and fast approach to predict inhomogeneous density and potential profiles of confined fluids. We use EQT potentials to construct a grand potential functional for cDFT. The EQT-cDFT-based grand potential can be used to predict various thermodynamic properties of confined fluids. In this work, we demonstrate the EQT-cDFT approach by simulating Lennard-Jones fluids, namely, methane and argon, confined inside slit-like channels of graphene. We show that the EQT-cDFT can accurately predict the structure and thermodynamic properties, such as density profiles, adsorption, local pressure tensor, surface tension, and solvation force, of confined fluids as compared to the molecular dynamics simulation results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2026

Beyond $\rho ^{2/3}$ scaling: microscopic origins and multimessengers of high-density nuclear symmetry energy

The nature and Equation of State (EOS) of dense neutron-rich matter are still very poorly known, while they have broad impacts on many interesting issues in both astrophysics and nuclear physics. In particular, the nuclear symmetry energy $E_{\textrm{sym}}(\rho )$ encoding the cost to make nuclear matter more neutron-rich has been the most uncertain component of the EOS of dense neutron-rich nucleonic matter. It significantly affects the radii, tidal deformations, cooling rates, and frequencies of various oscillation modes of isolated neutron stars as well as the strain amplitude and frequencies of gravitational waves from their mergers, besides its many effects on the structures of nuclei as well as the dynamics and observables of their collisions. Siemens (1970s) observed that $E_{\textrm{sym}}(\rho )$ scales as $(\rho /\rho _0)^{2/3}$ near the saturation density $\rho _0$ of nuclear matter, since both the kinetic part and the potential contribution (quadratic in momentum) exhibit this dependence. The scaling holds if: (1) the nucleon isoscalar potential is quadratic in momentum, and (2) the isovector interaction is weakly density-dependent. After examining many empirical evidences and understanding theoretical findings in the literature, we conclude that: (1) Siemens’ $\rho ^{2/3}$ scaling is robust and serves as a valuable benchmark for both nuclear theories and experiments up to $2\rho _0$ but breaks down at higher densities, (2) Experimental and theoretical findings about $E_{\textrm{sym}}(\rho )$ up to $2\rho _0$ are broadly consistent, but uncertainties remain large for its curvature $K_{\textrm{sym}}(\rho )$ and higher-order parameters, (3) Above $2\rho _0$ , uncertainties grow due to poorly constrained spin-isospin-dependent tensor and three-body forces as well as the resulting nucleon short-range correlations. Looking forward, combining signals from both observations of neutron stars and terrestrial heavy-ion reaction experiments is the most promising path to finally constraining the high-density $E_{\textrm{sym}}(\rho )$ and the EOS of supradense neutron-rich matter. Multiple examples of community efforts to further constrain the high-density $E_{\textrm{sym}}(\rho )$ using both real and mocked data of present and future high-precision observations of neutron stars, as well as heavy-ion collisions involving high-energy rare isotopes, are briefly reviewed.

Li, Bao-An [East Texas A & M Univ., Commerce, TX (↗
NASA NTRS1992

Comments on the present state of second-order closure models for incompressible flows

Second-order closure models account for history and nonlocal effects of the mean velocity gradients on the Reynolds stress tensor. Turbulent flows involving body forces or curvature, Reynolds stress relaxational effects, and counter-gradient transport are usually better described. The topics are presented in viewgraph form and include: (1) the Reynolds stress transport equation; (2) issues in second-order closure modeling; and (3) near wall models.

Speziale, Charles G.↗
NASA NTRS1997

Dibromine Monoxide, Br2O: The Rotational Spectrum and Molecular Properties

The rotational spectra of (79)Br2O, (79)BrO(81)Br, and Br2O in their ground vibrational states as well as (79)BrO(81)Br in its v (sub 2) = 1 state have been studied in selected regions between 90 and 523 GHz. Transitions involving a large range of quantum numbers, 6 less than or equal to J less than or equal to 123 and 0 less than or equal to K (sub a) less than or equal to 12, have been observed permitting precise rotational and a large set of centrifugal distortion constants to be determined. All isotopic species as well as the excited state data were fit simultaneously. Ground-state effective and average structural parameters as well as an estimate of the equilibrium structure have been derived. The quartic distortion constants were used for a calculation of the harmonic force field. The complete quadrupole tensor has been determined. Its diagonalization reveals a largely covalent BrO bond with little pi-bonding. The derived properties of Br2O are compared with those of related compounds such as Cl2O, HOBR, and HOCl.

Mueller, Holger S. P.↗
DOE OSTI2024

Numerical modelling of impact seismic sources using the stress glut theory

SUMMARY Meteorite impacts have proved to be a significant source of seismic signal on the Moon, and have now been recorded on Mars by InSight seismometers. Understanding how impacts produce seismic signal is key to the interpretation of this unique data, and to improve their identification in continuous seismic records. Here, we use the seismic Representation Theorem, and particularly the stress glut theory, to model the seismic motion resulting from impact cratering. The source is described by equivalent forces, some resulting from the impactor momentum transfer, and others from the stress glut, which represents the mechanical effect of plasticity and non linear processes in the source region. We condense these equivalent forces into a point-source with a time-varying single force and nine-component moment tensor. This analytical representation bridges the gap between the complex dynamics of crater formation, and the linear point-source representation classically used in seismology. Using the multiphysics modelling software HOSS, we develop a method to compute the stress glut of an impact, and the associated point-source from hypervelocity impact simulations. For a vertical and an oblique impact at 1000 m s−1, we show that the moment tensor presents a significant deviatoric component. Hence, the source is not an ideal isotropic explosion contrary to previous assumptions, and draws closer to a double couple for the oblique impact. The contribution of the point force to the seismic signal appears negligible. We verify this model by comparing two signals: (1) HOSS is coupled to SPECFEM3D to propagate the near-source signal elastically to remote seismic stations; (2) the point-source model derived from the stress-glut theory is used to generate displacements at the same distance. The comparison shows that the point-source model is accurately simulating the low-frequency impact seismic waveform, and its seismic moment is in trend with Lunar and Martian impact data. High-frequencies discrepancies exist, which are partly related to finite-source effects, but might be further explained by the difference in mathematical framework between classical seismology and HOSS’ numerical modelling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Quantum stress and torsion distributions in the deuteron

Stress distributions in the deuteron are related to form factors of the asymmetric energy-momentum tensor through three-dimensional Fourier transforms. There are eleven such form factors, which we calculate in an impulse approximation. We compare the obtained form factors to prior results for the six form factors that have been previously calculated. We then elaborate on the formalism for relating the form factors to internal distributions of mass, mass flux, momentum, stresses, and forces, and obtain results for all of these distributions. We obtain the principal stresses for the symmetric part of the stress tensor, and show that the antisymmetric part describes the reorientation of fermion spin by torsion stress when the nucleon moves between the S and D waves. Force distributions in the nucleons depend on the so-called nonconserved form factors through the Cauchy momentum equation, and are nonradial owing to the presence of tensor forces and spin-orbit coupling.

QCD phenomenology↗
DOE OSTI2015

Extracting the diffusion tensor from molecular dynamics simulation with Milestoning

We propose an algorithm to extract the diffusion tensor from Molecular Dynamics simulations with Milestoning. A Kramers-Moyal expansion of a discrete master equation, which is the Markovian limit of the Milestoning theory, determines the diffusion tensor. To test the algorithm, we analyze overdamped Langevin trajectories and recover a multidimensional Fokker-Planck equation. The recovery process determines the flux through a mesh and estimates local kinetic parameters. Rate coefficients are converted to the derivatives of the potential of mean force and to coordinate dependent diffusion tensor. We illustrate the computation on simple models and on an atomically detailed system—the diffusion along the backbone torsions of a solvated alanine dipeptide.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
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 OSTI2022

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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Self-interaction in the Bopp–Podolsky electrodynamics: Can the observable mass of a charged particle depend on its acceleration?

In this paper we obtain the expression for the self-force in the model with the Lagrangian containing additional terms, quadratic in Maxwell tensor derivatives (so-called Bopp–Podolsky electrodynamics). Features of this force are analyzed for various limiting cases. When a charged particle moves along straight line with a uniform acceleration, an explicit formula is found. In the framework of the considered model, an observable renormalized particle mass is shown to depend on its acceleration. This dependence allows, in principle, to extract experimentally a value of the particle bare mass. -- Highlights: •An expression for the self-force in the Bopp–Podolsky electrodynamics is given. •For a uniformly accelerated charged particle an explicit formula for the self-force is obtained. •Dependence between the observable mass of a charged particle and its acceleration is found.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Nature of Non-Adiabatic Electron-Ion Forces in Liquid Metals

We report that an accurate description of electron-ion interactions in materials is crucial for our understanding of their equilibrium and non-equilibrium properties. Here, we assess the properties of frictional forces experienced by ions in non-crystalline metallic systems, including liquid metals and warm dense plasmas, that arise from electronic excitations driven by the nuclear motion due to the presence of a continuum of low-lying electronic states. To this end, we perform detailed ab-initio calculations of the full friction tensor that characterizes the set of friction forces. The nonadiabatic electron-ion interactions introduce hydrodynamic couplings between the ionic degrees of freedom, which are sizeable between nearest neighbors. The friction tensor is generally inhomogeneous, anisotropic and non-diagonal, especially at lower densities.

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

A stable partitioned FSI algorithm for rigid bodies and incompressible flow. Part I: Model problem analysis

A stable partitioned algorithm is developed for fluid-structure interaction (FSI) problems involving viscous incompressible flow and rigid bodies. This added-mass partitioned (AMP) algorithm remains stable, without sub-iterations, for light and even zero mass rigid bodies when added-mass and viscous added-damping effects are large. The scheme is based on a generalized Robin interface condition for the fluid pressure that includes terms involving the linear acceleration and angular acceleration of the rigid body. Added mass effects are handled in the Robin condition by inclusion of a boundary integral term that depends on the pressure. Added-damping effects due to the viscous shear forces on the body are treated by inclusion of added-damping tensors that are derived through a linearization of the integrals defining the force and torque. Added-damping effects may be important at low Reynolds number, or, for example, in the case of a rotating cylinder or rotating sphere when the rotational moments of inertia are small. In this second part of a two-part series, the general formulation of the AMP scheme is presented including the form of the AMP interface conditions and added-damping tensors for general geometries. A fully second-order accurate implementation of the AMP scheme is developed in two dimensions based on a fractional-step method for the incompressible Navier-Stokes equations using finite difference methods and overlapping grids to handle the moving geometry. Here, the numerical scheme is verified on a number of difficult benchmark problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Avalanche boron fusion by laser picosecond block ignition with magnetic trapping for clean and economic reactor

Measured highly elevated gains of proton–boron (HB11) fusion (Picciotto et al. , Phys. Rev. X 4 , 031030 (2014)) confirmed the exceptional avalanche reaction process (Lalousis et al. , Laser Part. Beams 32 , 409 (2014); Hora et al. , Laser Part. Beams 33 , 607 (2015)) for the combination of the non-thermal block ignition using ultrahigh intensity laser pulses of picoseconds duration. The ultrahigh acceleration above $10^{20}~\text{cm}~\text{s}^{-2}$ for plasma blocks was theoretically and numerically predicted since 1978 (Hora, Physics of Laser Driven Plasmas (Wiley, 1981), pp. 178 and 179) and measured (Sauerbrey, Phys. Plasmas 3 , 4712 (1996)) in exact agreement (Hora et al. , Phys. Plasmas 14 , 072701 (2007)) when the dominating force was overcoming thermal processes. This is based on Maxwell’s stress tensor by the dielectric properties of plasma leading to the nonlinear (ponderomotive) force $f_{\text{NL}}$ resulting in ultra-fast expanding plasma blocks by a dielectric explosion. Combining this with measured ultrahigh magnetic fields and the avalanche process opens an option for an environmentally absolute clean and economic boron fusion power reactor. Finally, this is supported also by other experiments with very high HB11 reactions under different conditions (Labaune et al. , Nature Commun. 4 , 2506 (2013)).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗