Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 1,189 records · Page 66
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Transport coefficients for magnetic-field evolution in inviscid magnetohydrodynamics

The magnetized resistivity and electrothermal tensors when substituted into the induction equation lead to electrothermal magnetic field generation, resistive magnetic diffusion, and magnetic field advection due to resistivity gradients, temperature gradients and currents. Here, the advection terms driven by temperature gradient and current have cross-field components (perpendicular to both the magnetic field and the driving term) that depend on significantly modified versions of Braginskii’s transport coefficients. The improved fits to Braginskii’s coefficients given by Epperlein and Haines and Ji and Held give physically incorrect results for cross-field advection at small Hall parameters (product of cyclotron frequency and collision time). The errors in Epperlein and Haines’ fits are particularly severe, giving increasing advection velocities below a Hall parameter of one when they should decrease linearly to zero. Epperlein and Haines’ fits can also give erroneous advection terms due to variations in effective atomic number. The only serious error in Braginskii’s fits is an overestimate in advection due to perpendicular resistivity. New fits for the cross-field advection terms are obtained from a direct numerical solution of the Fokker-Planck equation and Ji and Held’s higher order expansion approach that are continuous functions of the effective atomic number.

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

DECaNT: Simulation tool for diffusion of excitons in carbon nanotube films

In this work, we present the numerical tool DECaNT (Diffusion of Excitons in Carbon NanoTubes) that simulates exciton transport in thin films of carbon nanotubes. Through a mesh of nanotubes generated using the Bullet Physics C++ library, excitons move according to an ensemble Monte Carlo algorithm, with the scattering rates that account for tube chirality, orientation, and distance. We calculate the diffusion tensor from the position–position correlation functions and analyze its anisotropy and dependence on the film composition, morphology, and defect density.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Reduced scaling formulation of CASPT2 analytical gradients using the supporting subspace method

We present a reduced scaling and exact reformulation of state specific complete active space second-order perturbation (CASPT2) analytical gradients in terms of the MP2 and Fock derivatives using the supporting subspace method. This work follows naturally from the supporting subspace formulation of the CASPT2 energy in terms of the MP2 energy using dressed orbitals and Fock builds. For a given active space configuration, the terms corresponding to the MP2-gradient can be evaluated with O(N5) operations, while the rest of the calculations can be computed with O(N3) operations using Fock builds, Fock gradients, and linear algebra. When tensor-hyper-contraction is applied simultaneously, the computational cost can be further reduced to O(N4) for a fixed active space size. The new formulation enables efficient implementation of CASPT2 analytical gradients by leveraging the existing graphical processing unit (GPU)-based MP2 and Fock routines. We present benchmark results that demonstrate the accuracy and performance of the new method. Example applications of the new method in ab initio molecular dynamics simulation and constrained geometry optimization are given.

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

The ultrafast Kerr effect in anisotropic and dispersive media

The ultrafast optical Kerr effect (OKE) is widely used to investigate the structural dynamics and interactions of liquids, solutions, and solids by observing their intrinsic nonlinear temporal responses through nearly collinear four-wave mixing. Non-degenerate mixing schemes allow for background free detection and can provide information on the interplay between a material’s internal degrees of freedom. Here, we show a source of temporal dynamics in the OKE signal that is not reflective of the internal degrees of freedom but arises from a group index and momentum mismatch. It is observed in two-color experiments on condensed media with sizable spectral dispersion, a common property near an optical resonance. In particular, birefringence in crystalline solids is able to entirely change the character of the OKE signal via the off-diagonal tensor elements of the nonlinear susceptibility. We develop a detailed description of the phase-mismatched ultrafast OKE and show how to extract quantitative information on the spectrally resolved birefringence and group index from time-resolved experiments in one and two dimensions.

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

Gyrokinetic and extended-MHD simulations of a flow shear stabilized Z-pinch experiment

Axisymmetric (m = 0) gyrokinetic and extended-MHD simulations of a sheared-flow Z-pinch plasma are performed with the high-order finite volume code COGENT. The present gyrokinetic model solves the long-wavelength limit of the gyrokinetic equation for both ion and electron species coupled to the electrostatic gyro-Poisson equation for the electrostatic potential. The extended-MHD model is electromagnetic and includes the effects of the gyro-viscous pressure tensor, diamagnetic electron and ion heat fluxes, and generalized Ohm's law. Additionally, a prominent feature of this work is that the radial profiles for the plasma density and temperature are taken from the fusion Z-pinch experiment (FuZE), and the magnetic field profile is obtained as a solution of the MHD force balance equation. Such an approach allows to address realistic plasma parameters and provide insights into the current and planned experiments. In particular, it is demonstrated that the radial profiles play an important role in stabilization, as the embedded guiding center (E×B ) drift has a strong radial shear, which can contribute to the Z-pinch stabilization even in the absence of the fluid flow shear. The results of simulations for the FuZE plasma parameters show a decrease in the linear growth rate with an increase in the flow shear; however, full stabilization in the linear regime is not observed even for large (comparable to the Alfvén velocity) radial variations of the axial flow. Nonlinear stability properties of the FuZE plasmas are also studied, and it is found that profile broadening can have a pronounced stabilizing effect in the nonlinear regime.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Extended Lagrangian Born–Oppenheimer molecular dynamics for orbital-free density-functional theory and polarizable charge equilibration models

We report extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) is formulated for orbital-free Hohenberg–Kohn density-functional theory and for charge equilibration and polarizable force-field models that can be derived from the same orbital-free framework. The purpose is to introduce the most recent features of orbital-based XL-BOMD to molecular dynamics simulations based on charge equilibration and polarizable force-field models. These features include a metric tensor generalization of the extended harmonic potential, preconditioners, and the ability to use only a single Coulomb summation to determine the fully equilibrated charges and the interatomic forces in each time step for the shadow Born–Oppenheimer potential energy surface. The orbital-free formulation has a charge-dependent, short-range energy term that is separate from long-range Coulomb interactions. This enables local parameterizations of the short-range energy term, while the long-range electrostatic interactions can be treated separately. The theory is illustrated for molecular dynamics simulations of an atomistic system described by a charge equilibration model with periodic boundary conditions. The system of linear equations that determines the equilibrated charges and the forces is diagonal, and only a single Ewald summation is needed in each time step. The simulations exhibit the same features in accuracy, convergence, and stability as are expected from orbital-based XL-BOMD.

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 OSTI2021

Cartesian message passing neural networks for directional properties: Fast and transferable atomic multipoles

The message passing neural network (MPNN) framework is a promising tool for modeling atomic properties but is, until recently, incompatible with directional properties, such as Cartesian tensors. We propose a modified Cartesian MPNN (CMPNN) suitable for predicting atom-centered multipoles, an essential component of ab initio force fields. The efficacy of this model is demonstrated on a newly developed dataset consisting of 46 623 chemical structures and corresponding high-quality atomic multipoles, which was deposited into the publicly available Molecular Sciences Software Institute QCArchive server. We show that the CMPNN accurately predicts atom-centered charges, dipoles, and quadrupoles and that errors in the predicted atomic multipoles have a negligible effect on multipole–multipole electrostatic energies. The CMPNN is accurate enough to model conformational dependencies of a molecule’s electronic structure. This opens up the possibility of recomputing atomic multipoles on the fly throughout a simulation in which they might exhibit strong conformational dependence.

Glick, Zachary L. (ORCID:0000000309002849)↗
DOE OSTI2021

Measurement of Local Magnetic Fields in Actinide Tetrafluorides

Fluorine-19 magnetic shielding tensors have been measured in a series of actinide tetrafluorides (AnF4) by solid state nuclear magnetic resonance (NMR) spectroscopy. Tetravalent actinide centers with 0 to 8 valence electrons can form tetrafluorides with the same unit cell, making these compounds an attractive choice for a systematic study of the variation in electronic structure across the 5f row of the Periodic Table. Pronounced deviations from predictions based on localized valence electron models have been detected by these experiments, which suggests that this approach may be used as a quantitative probe of electronic correlations.

Walter, Eric D.↗
DOE OSTI2021

A classification of invertible phases of bosonic quantum lattice systems in one dimension

We study invertible states of 1D bosonic quantum lattice systems. Here we show that every invertible 1D state is in a trivial phase: after tensoring with some unentangled ancillas, it can be disentangled by a fuzzy analog of a finite-depth quantum circuit. If an invertible state has symmetries, it may be impossible to disentangle it in a way that preserves the symmetries, even after adding unentagled ancillas. We show that in the case of a finite unitary symmetry G , the only obstruction is an index valued in degree-2 cohomology of G . We show that two invertible G -invariant states are in the same phase if and only if their indices coincide.

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

Wave-driven torques to drive current and rotation

In the classic Landau damping initial value problem, where a planar electrostatic wave transfers energy and momentum to resonant electrons, a recoil reaction occurs in the nonresonant particles to ensure momentum conservation. To explain how net current can be driven in spite of this conservation, the literature often appeals to mechanisms that transfer this nonresonant recoil momentum to ions, which carry negligible current. However, this explanation does not allow the transport of net charge across magnetic field lines, precluding E X B rotation drive. In this study, we show that in steady state, this picture of current drive is incomplete. Using a simple Fresnel model of the plasma, we show that for lower hybrid waves, the electromagnetic energy flux (Poynting vector) and momentum flux (Maxwell stress tensor) associated with the evanescent vacuum wave become the Minkowski energy flux and momentum flux in the plasma and are ultimately transferred to resonant particles. Thus, the torque delivered to the resonant particles is ultimately supplied by the electromagnetic torque from the antenna, allowing the nonresonant recoil response to vanish and rotation to be driven. We present a warm fluid model that explains how this momentum conservation works out locally, via a Reynolds stress that does not appear in the one-dimensional initial value problem. This model is the simplest that can capture both the nonresonant recoil reaction in the initial-value problem, and the absence of a nonresonant recoil in the steady-state boundary value problem, thus forbidding rotation drive in the former while allowing it in the latter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Engineering local strain for single-atom nuclear acoustic resonance in silicon

Mechanical strain plays a key role in the physics and operation of nanoscale semiconductor systems, including quantum dots and single-dopant devices. In this report we describe the design of a nanoelectronic device, where a single nuclear spin is coherently controlled via nuclear acoustic resonance (NAR) through the local application of dynamical strain. The strain drives spin transitions by modulating the nuclear quadrupole interaction. We adopt an AlN piezoelectric actuator compatible with standard silicon metal–oxide–semiconductor processing and optimize the device layout to maximize the NAR drive. We predict NAR Rabi frequencies of order 200 Hz for a single 123Sb nucleus in a wide region of the device. Spin transitions driven directly by electric fields are suppressed in the center of the device, allowing the observation of pure NAR. Using electric field gradient-elastic tensors calculated by the density-functional theory, we extend our predictions to other high-spin group-V donors in silicon and to the isoelectronic 73Ge atom.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Symmetry-specific orientational order parameters for complex structures

A comprehensive framework of characterizing complex self-assembled structures with a set of orientational order parameters is presented. It is especially relevant in the context of using anisotropic building blocks with various symmetries. Two classes of tensor order parameters are associated with polyhedral nematic and bond orientational order. For the latter, a variation of classical bond order parameters (BOPs) is introduced, which takes advantage of the symmetry of constituent particles and/or expected crystalline phases. These symmetrized BOPs can be averaged over an entire system or assigned locally to an individual bond. By combining that with the bond percolation procedure, one is able to identify coherent domains within a self-assembled structure. As a demonstration of the proposed framework, we apply it to a simulated hybrid system that combines isotropic and patchy particles with octahedral symmetry. In conclusion, not only does the methodology allow one to identify individual crystalline domains but also it detects coherent clusters of a peculiar compact amorphous structure that is not space-filling and lacks any long-range order.

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

An incompressible Eulerian method for fluid–structure interaction with mixed soft and rigid solids

Here, we present a general simulation approach for incompressible fluid–structure interactions in a fully Eulerian framework using the reference map technique. The approach is suitable for modeling one or more rigid or finitely deformable objects or soft objects with rigid components interacting with the fluid and with each other. It is also extended to control the kinematics of structures in fluids. The model is based on our previous Eulerian fluid–soft solver [Rycroft et al., “Reference map technique for incompressible fluid–structure interaction,” J. Fluid Mech. 898, A9 (2020)] and generalized to rigid structures by constraining the deformation-rate tensor in a projection framework. Several numerical examples are presented to illustrate the capability of the method.

42 ENGINEERING↗
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 OSTI2022

Polar liquids at charged interfaces: A dipolar shell theory

The structure of polar liquids and electrolytic solutions, such as water and aqueous electrolytes, at interfaces underlies numerous phenomena in physics, chemistry, biology, and engineering. In this work, we develop a continuum theory that captures the essential features of dielectric screening by polar liquids at charged interfaces, including decaying spatial oscillations in charge and mass, starting from the molecular properties of the solvent. The theory predicts an anisotropic dielectric tensor of interfacial polar liquids previously studied in molecular dynamics simulations. We explore the effect of the interfacial polar liquid properties on the capacitance of the electrode/electrolyte interface and on hydration forces between two plane-parallel polarized surfaces. In the linear response approximation, we obtain simple formulas for the characteristic decay lengths of molecular and ionic profiles at the interface.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗