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 91 records · Page 5

Solar Sail Torque Model Characterization for the Near Earth Asteroid Scout Mission

Near Earth Asteroid Scout (NEA Scout) was a mission to test solar sail propulsion for orbital transfer from cislunar space to flyby and image an asteroid. Had it succeeded, one of the mission goals was to characterize the solar torque on the sail to ensure successful attitude control for the orbit transfer and imaging the asteroid. The simulation used to develop the flight attitude control software uses the generalized model for solar sails, a tensor equation of the forces and torques on sails of arbitrary shape. Rios-Reyes and Scheeres developed a general process to update the torque tensor coefficients using estimates of sail torque over a range of directions to the sun. Their process was adapted and implemented for the specific case of NEA Scout using spacecraft telemetry collected during sail characterization maneuvers in combination with simulation models and parameters. The NEA Scout maneuvers were limited to the operating range of the mission and constraints of the control hardware and allowed safe testing of each attitude before proceeding to the next. The NEA Scout reaction wheel speeds are used to measure accumulated momentum, while the Active Mass Translator (AMT) position is used to subtract out the torque from the center of mass crossed with the sail force and isolate the torque from only the sail shape. The process was tested by running attitude control simulations of the characterization maneuvers, generating simulated telemetry, estimating the solar torques, then using a least squares estimating the solar torque coefficients using least-squares and then performing a least-squares fit to the solar torque tensor coefficients. These estimated coefficients were tested by evaluating the solar torques under the same conditions as the simulated telemetry and comparing to the true simulated torques. Solar force model updates can be performed separately by observing the effect of the sail on the trajectory, and the torque model can be refined using those solar force updates. This process met the needs of the NEA Scout mission and can be adapted to characterize the solar torque for other missions with different sails.

solar sail↗

Solar Sail Torque Model Characterization for the Near Earth Asteroid Scout Mission

Near Earth Asteroid Scout (NEA Scout) was a mission to test solar sail propulsion for orbital transfer from cislunar space to flyby and image an asteroid. One of the goals of the mission was to characterize the solar torque on the sail to ensure successful attitude control for the orbit transfer and imaging the asteroid. The simulation used to develop the flight attitude control software uses the generalized model for solar sails, a tensor equation of the forces and torques on sails of arbitrary shape. Rios-Reyes and Scheeres developed a general process to update the torque tensor coefficients using estimates of sail torque over a range of directions to the sun. Their process was adapted and implemented for the specific case of NEA Scout using spacecraft telemetry collected during sail characterization maneuvers in combination with simulation models and parameters. The NEA Scout maneuvers were limited to the operating range of the mission and constraints of the control hardware and allowed safe testing of each attitude before proceeding to the next. The NEA Scout reaction wheel speeds are used to measure accumulated momentum, while the Active Mass Translator (AMT) position is used to subtract out the torque from the center of mass crossed with the sail force and isolate the torque from only the sail shape. The process was tested by running attitude control simulations of the characterization maneuvers, generating simulated telemetry, estimating the solar torques, then using a least squares estimating the solar torque coefficients using least-squares and then performing a least-squares fit to the solar torque tensor coefficients. These estimated coefficients were tested by evaluating the solar torques under the same conditions as the simulated telemetry and comparing to the true simulated torques. Solar force model updates can be performed separately by observing the effect of the sail on the trajectory, and the torque model can be refined using those solar force updates. This process met the needs of the NEA Scout mission and can be adapted to characterize the solar torque for other missions with different sails.

solar sail↗

Screening the fifth force in the Horndeski's most general scalar-tensor theories

We study how the Vainshtein mechanism operates in the most general scalar-tensor theories with second-order equations of motion. The field equations of motion, which can be also applicable to most of other screening scenarios proposed in literature, are generally derived in a spherically symmetric space-time with a matter source. In the presence of a field coupling to the Ricci scalar, we clarify conditions under which the Vainshtein mechanism is at work in a weak gravitational background. We also obtain the solutions of the field equation inside a spherically symmetric body and show how they can be connected to exterior solutions that accommodate the Vainshtein mechanism. We apply our general results to a number of concrete models such as the covariant/extended Galileons and the DBI Galileons with Gauss-Bonnet and other terms. In these models the fifth force can be suppressed to be compatible with solar-system constraints, provided that non-linear field kinetic terms coupled to the Einstein tensor do not dominate over other non-linear field self-interactions.

79 ASTRONOMY AND ASTROPHYSICS↗

Generalized and efficient algorithm for computing multipole energies and gradients based on Cartesian tensors

Accurate representation of intermolecular forces has been the central task of classical atomic simulations, known as molecular mechanics. Recent advancements in molecular mechanics models have put forward the explicit representation of permanent and/or induced electric multipole (EMP) moments. The formulas developed so far to calculate EMP interactions tend to have complicated expressions, especially in Cartesian coordinates, which can only be applied to a specific kernel potential function. For example, one needs to develop a new formula each time a new kernel function is encountered. The complication of these formalisms arises from an intriguing and yet obscured mathematical relation between the kernel functions and the gradient operators. Here, I uncover this relation via rigorous derivation and find that the formula to calculate EMP interactions is basically invariant to the potential kernel functions as long as they are of the form f(r), i.e., any Green’s function that depends on inter-particle distance. I provide an algorithm for efficient evaluation of EMP interaction energies, forces, and torques for any kernel f(r) up to any arbitrary rank of EMP moments in Cartesian coordinates. The working equations of this algorithm are essentially the same for any kernel f(r). Recently, a few recursive algorithms were proposed to calculate EMP interactions. Depending on the kernel functions, the algorithm here is about 4–16 times faster than these algorithms in terms of the required number of floating point operations and is much more memory efficient. I show that it is even faster than a theoretically ideal recursion scheme, i.e., one that requires 1 floating point multiplication and 1 addition per recursion step. This algorithm has a compact vector-based expression that is optimal for computer programming. The Cartesian nature of this algorithm makes it fit easily into modern molecular simulation packages as compared with spherical coordinate-based algorithms. A software library based on this algorithm has been implemented in C++11 and has been released.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Tests of modified gravity with dwarf galaxies

In modified gravity theories that seek to explain cosmic acceleration, dwarf galaxies in low density environments can be subject to enhanced forces. The class of scalar-tensor theories, which includes f(R) gravity, predict such a force enhancement (massive galaxies like the Milky Way can evade it through a screening mechanism that protects the interior of the galaxy from this ''fifth'' force). We study observable deviations from GR in the disks of late-type dwarf galaxies moving under gravity. The fifth-force acts on the dark matter and HI gas disk, but not on the stellar disk owing to the self-screening of main sequence stars. We find four distinct observable effects in such disk galaxies: 1. A displacement of the stellar disk from the HI disk. 2. Warping of the stellar disk along the direction of the external force. 3. Enhancement of the rotation curve measured from the HI gas compared to that of the stellar disk. 4. Asymmetry in the rotation curve of the stellar disk. We estimate that the spatial effects can be up to 1 kpc and the rotation velocity effects about 10 km/s in infalling dwarf galaxies. Such deviations are measurable: we expect that with a careful analysis of a sample of nearby dwarf galaxies one can improve astrophysical constraints on gravity theories by over three orders of magnitude, and even solar system constraints by one order of magnitude. Thus effective tests of gravity along the lines suggested by Hui, Nicolis, and Stubbs (2009) and Jain (2011) can be carried out with low-redshift galaxies, though care must be exercised in understanding possible complications from astrophysical effects.

79 ASTRONOMY AND ASTROPHYSICS↗

Configurational forces in electronic structure calculations using Kohn-Sham density functional theory

Here we derive the expressions for configurational forces in Kohn-Sham density functional theory, which correspond to the generalized variational force computed as the derivative of the Kohn-Sham energy functional with respect to the position of a material point x. These configurational forces that result from the inner variations of the Kohn-Sham energy functional provide a unified framework to compute atomic forces as well as stress tensor for geometry optimization. Importantly, owing to the variational nature of the formulation, these configurational forces inherently account for the Pulay corrections. The formulation presented in this work treats both pseudopotential and all-electron calculations in a single framework, and employs a local variational real-space formulation of Kohn-Sham density functional theory (DFT) expressed in terms of the nonorthogonal wave functions that is amenable to reduced-order scaling techniques. We demonstrate the accuracy and performance of the proposed configurational force approach on benchmark all-electron and pseudopotential calculations conducted using higher-order finite-element discretization. To this end, we examine the rates of convergence of the finite-element discretization in the computed forces and stresses for various materials systems, and, further, verify the accuracy from finite differencing the energy. Wherever applicable, we also compare the forces and stresses with those obtained from Kohn-Sham DFT calculations employing plane-wave basis (pseudopotential calculations) and Gaussian basis (all-electron calculations). Finally, we verify the accuracy of the forces on large materials systems involving a metallic aluminum nanocluster containing 666 atoms and an alkane chain containing 902 atoms, where the Kohn-Sham electronic ground state is computed using a reduced-order scaling subspace projection technique

36 MATERIALS SCIENCE↗

A priori testing of subgrid-scale models for the velocity-pressure and vorticity-velocity formulations

Subgrid-scale models for Large Eddy Simulation (LES) in both the velocity-pressure and the vorticity-velocity formulations were evaluated and compared in a priori tests using spectral Direct Numerical Simulation (DNS) databases of isotropic turbulence: 128(exp 3) DNS of forced turbulence (Re(sub(lambda))=95.8) filtered, using the sharp cutoff filter, to both 32(exp 3) and 16(exp 3) synthetic LES fields; 512(exp 3) DNS of decaying turbulence (Re(sub(Lambda))=63.5) filtered to both 64(exp 3) and 32(exp 3) LES fields. Gaussian and top-hat filters were also used with the 128(exp 3) database. Different LES models were evaluated for each formulation: eddy-viscosity models, hyper eddy-viscosity models, mixed models, and scale-similarity models. Correlations between exact versus modeled subgrid-scale quantities were measured at three levels: tensor (traceless), vector (solenoidal 'force'), and scalar (dissipation) levels, and for both cases of uniform and variable coefficient(s). Different choices for the 1/T scaling appearing in the eddy-viscosity were also evaluated. It was found that the models for the vorticity-velocity formulation produce higher correlations with the filtered DNS data than their counterpart in the velocity-pressure formulation. It was also found that the hyper eddy-viscosity model performs better than the eddy viscosity model, in both formulations.

Winckelmans, G. S.↗

Few-body physics with CLAS

The study of few-body, nuclear systems with electromagnetic probes is an essential piece of the scientific program at Jefferson Lab. Reactions using real photons and electrons (up to energies of 6 GeV) are measured using the CEBAF large acceptance spectrometer (CLAS) detector in Hall B, a nearly 4π magnetic spectrometer. We focus here on three areas. (1) Short-range correlations (SRCs) probe the high-momentum components of the nuclear wave function. Recent CLAS experiments map out their isospin character and reveal the importance of the tensor part of the nuclear force. (2) Three-body forces are an essential feature of nuclei. We will show results using real photons and 3He and 4He targets that remain largely unexplained. (3) Evidence for the transition to a quark-gluon description of nuclei has been observed with photon beams in CLAS on deuterium and 3-He targets. Alternative explanations reveal the geography of the transition is complex.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Effect of the tensor part of Skyrme interaction on the description of elastic nucleon-nucleus scattering on the basis of the optical model

A microscopic optical nucleon-nucleus potential constructed on the basis of calculating the mass operator for the single-particle Green's function with Skyrme nucleon-nucleon forces was used to study the effect of the tensor part of Skyrme forces on describing differential cross sections for elastic nucleon-nucleus scattering and the structure of nuclei within a self-consistent approach. It was shown that an increase in the tensor terms of nucleon-nucleon forces impaired the description of nucleon-nucleus scattering within the model being considered. The parameters of Skyrme forces were optimized on the basis of an analysis of cross sections for neutron-nucleus scattering, the properties of nuclear matter and the structure of nuclei being tested in doing this. This optimization led to nucleon-nucleon force versions where the tensor part was small or zero and which described satisfactorily the structure of nuclei and cross sections for elastic neutron and proton scattering on nuclei over a broad range of target mass numbers.

74 ATOMIC AND MOLECULAR PHYSICS↗

New Examples of Integrable Systems with Dissipation on the Tangent Bundles of Multidimensional Spheres

In many problems of multidimensional dynamics, systems appear whose state spaces are spheres of finite dimension. Clearly, phase spaces of such systems are tangent bundles of these spheres. In this paper, we examine nonconservative force fields in the dynamics of a multidimensional rigid body in which the system possesses a complete set of first integrals that can be expressed as finite combinations of elementary transcendental functions. We consider the case where the moment of nonconservative forces depends on the tensor of angular velocity.

97 MATHEMATICS AND COMPUTING↗

Seismic Waves, 4th order accurate

SW4 is a program for simulating seismic wave propagation on parallel computers. SW4 colves the seismic wave equations in Cartesian corrdinates. It is therefore appropriate for regional simulations, where the curvature of the earth can be neglected. SW4 implements a free surface boundary condition on a realistic topography, absorbing super-grid conditions on the far-field boundaries, and a kinematic source model consisting of point force and/or point moment tensor source terms. SW4 supports a fully 3-D heterogeneous material model that can be specified in several formats. SW4 can output synthetic seismograms in an ASCII test format, or in the SAC finary format. It can also present simulation information as GMT scripts, whixh can be used to create annotated maps. Furthermore, SW4 can output the solution as well as the material model along 2-D grid planes.

Source record↗

Toward Global Regional Seismic Moment Tensor Inversion with Three-Dimensional Earth Models for Nuclear Explosion Monitoring with Sparse Networks: Demonstration of Reciprocity for Strain Greens Tensor Database Simulation with Salvus

Seismic source characterization is an essential function of global nuclear explosion monitoring (NEM). While large events (roughly with moment magnitude, M w , greater than 5.0) can often be easily detected, located and identified with high signal-to-noise ratios at teleseismic distances (> 20°), trends in NEM research require confident source characterization at much lower magnitudes (say down to 3.0) and exploitation of sparse observations (from only a few stations) at regional distance (< 20°). Regional distance waveform inversion to characterize sources is now widely used and effective (e.g. Ford et al., 2009; Alvizuri and Tape, 2018; Alvizuri et al., 2018; Chiang et al., 2018; Ford et al., 2022). These methods obtain the magnitude, depth and seismic moment tensor, which represents the forces that excited the observed seismic waves (slip on an earthquake fault, explosion, collapse or a combination of various forces). Common to many problems in seismology, the isolation of the source 2 properties requires removal of path propagation effects that waves experience while traveling through the three-dimensional (3D) Earth (the structure exists due to different rock types, material properties, temperature and tectonic processes).

58 GEOSCIENCES↗

Geometrization of the Dirac theory of the electron

Using the concept of parallel displacement of a half vector, the Dirac equations are generally written in invariant form. The energy tensor is formed and both the macroscopic and quantum mechanic equations of motion are set up. The former have the usual form: divergence of the energy tensor equals the Lorentz force and the latter are essentially identical with those of the geodesic line.

Fock, V.↗

Seismic Waves, 4th order accurate

SW4 is a program for simulating seismic wave propagation on parallel computers. SW4 colves the seismic wave equations in Cartesian corrdinates. It is therefore appropriate for regional simulations, where the curvature of the earth can be neglected. SW4 implements a free surface boundary condition on a realistic topography, absorbing super-grid conditions on the far-field boundaries, and a kinematic source model consisting of point force and/or point moment tensor source terms. SW4 supports a fully 3-D heterogeneous material model that can be specified in several formats. SW4 can output synthetic seismograms in an ASCII test format, or in the SAC finary format. It can also present simulation information as GMT scripts, whixh can be used to create annotated maps. Furthermore, SW4 can output the solution as well as the material model along 2-D grid planes.

Source record↗

Light pressure on a solid body immersed in a liquid medium

We have solved the problem about the force with which an electromagnetic pulse in a liquid (or gaseous) medium at rest affects a solid body (also at rest) immersed in it. We have shown that under certain conditions (relating to the characteristics of the medium and the pulse shape), the formula for the force exerted per unit area of a body surface is obtained from the Landau - Lifshitz equations for static fields in the same way as, according to Pitaevskii, the field stress tensor is obtained from the static field stress tensor with the dispersion taken into account. The formula for the force acting on the wall, from which an incident quasi-monochromatic plane wave with a given intensity is reflected, differs from the corresponding formula for the case when the body is in a vacuum by the factor {+-}n{sub 1}, where n{sub 1} is the refractive index, and the upper (lower) sign corresponds to a positive (negative) group velocity of the wave in the medium. (laser applications and other topics in quantum electronics)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Tensor factorization for elucidating mechanisms of piezoresponse relaxation via dynamic Piezoresponse Force Spectroscopy

Spatially resolved time and voltage-dependent polarization dynamics in PbTiO 3 thin films is explored using dynamic piezoresponse force microscopy (D-PFM) in conjunction with interferometric displacement sensing. This approach gives rise to 4D data sets containing information on bias-dependent relaxation dynamics at each spatial location without long-range electrostatic artifacts. To interpret these data sets in the absence of defined physical models, we employ a non-negative tensor factorization method which clearly presents the data as a product of simple behaviors allowing for direct physics interpretation. Correspondingly, we perform phase-field modeling finding the existence of ‘hard’ and ‘soft’ domain wall edges. This approach can be extended to other multidimensional spectroscopies for which even exploratory data analysis leads to unsatisfactory results due to many components in the decomposition.

36 MATERIALS SCIENCE↗

Gyroscopic tensor of a magnetic soliton

The gyroscopic force acting on a ferromagnetic soliton reflects the tendency of spins to precess. Its components increase linearly with (generalized) velocities of the soliton. The proportionality coefficients form the gyroscopic tensor, a generalization of Thiele’s gyrovector. Originally introduced as an auxiliary quantity, the gyroscopic tensor turns out to be a mathematical object of fundamental importance with a long history going back to Lagrange. Here, we review the applications of the gyroscopic tensor and its historical roots.

36 MATERIALS SCIENCE↗

Energy balance in peridynamics.

The peridynamic model of solid mechanics treats internal forces within a continuum through interactions across finite distances. These forces are determined through a constitutive model that, in the case of an elastic material, permits the strain energy density at a point to depend on the collective deformation of all the material within some finite distance of it. The forces between points are evaluated from the Frechet derivative of this strain energy density with respect to the deformation map. The resulting equation of motion is an integro-differential equation written in terms of these interparticle forces, rather than the traditional stress tensor field. Recent work on peridynamics has elucidated the energy balance in the presence of these long-range forces. We have derived the appropriate analogue of stress power, called absorbed power, that leads to a satisfactory definition of internal energy. This internal energy is additive, allowing us to meaningfully define an internal energy density field in the body. An expression for the local first law of thermodynamics within peridynamics combines this mechanical component, the absorbed power, with heat transport. The global statement of the energy balance over a subregion can be expressed in a form in which the mechanical and thermal terms contain only interactions between the interior of the subregion and the exterior, in a form anticipated by Noll in 1955. The local form of this first law within peridynamics, coupled with the second law as expressed in the Clausius-Duhem inequality, is amenable to the Coleman-Noll procedure for deriving restrictions on the constitutive model for thermomechanical response. Using an idea suggested by Fried in the context of systems of discrete particles, this procedure leads to a dissipation inequality for peridynamics that has a surprising form. It also leads to a thermodynamically consistent way to treat damage within the theory, shedding light on how damage, including the nucleation and advance of cracks, should be incorporated into a constitutive model.

42 ENGINEERING↗