Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “invariance”

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 631 records · Page 35

Viscoplastic Model Development to Account for Strength Differential: Application to Aged Inconel 718 at Elevated Temperature

The magnitude of yield and flow stresses in aged Inconel 718 are observed to be different in tension and compression. This phenomenon, called the Strength differential (SD), contradicts the metal plasticity axiom that the second deviatoric stress invariant alone is sufficient for representing yield and flow. Apparently, at least one of the other two stress invariants is also significant. A unified viscoplastic model was developed that is able to account for the SD effect in aged Inconel 718. Building this model involved both theory and experiments. First, a general threshold function was proposed that depends on all three stress invariants and then the flow and evolution laws were developed using a potential-based thermodynamic framework. Judiciously chosen shear and axial tests were conducted to characterize the material. Shear tests involved monotonic loading, relaxation, and creep tests with different loading rates and load levels. The axial tests were tension and compression tests that resulted in sufficiently large inelastic strains. All tests were performed at 650 C. The viscoplastic material parameters were determined by optimizing the fit to the shear tests, during which the first and the third stress invariants remained zero. The threshold surface parameters were then fit to the tension and compression test data. An experimental procedure was established to quantify the effect of each stress invariant on inelastic deformation. This requires conducting tests with nonproportional three-dimensional load paths. Validation of the model was done using biaxial tests on tubular specimens of aged Inconel 718 using proportional and nonproportional axial-torsion loading. These biaxial tests also helped to determine the most appropriate form of the threshold function; that is, how to combine the stress invariants. Of the set of trial threshold functions, the ones that incorporated the third stress invariant give the best predictions. However, inclusion of the first stress invariant does not significantly improve the model predictions. The model shows excellent predictive capability for nonproportional load paths. Additionally, it reduces to the well-known models of Mises Drucker and Drucker-Prager. The requisite experiments involve reasonably simple load paths in the axial-shear stress plane and hence can be performed on a variety of different materials: be they metallic, geological. polymeric, ceramic or granular. The general form of the threshold function allows representation of inelastic deformation in a range of materials.

Iyer, Saiganesh↗

Adiabatic theory in regions of strong field gradients

The theory for the generalized first invariant for adiabatic motion of charged particles in regions where there are large gradients in magnetic or electric fields is developed. The general condition for an invariant to exist in such regions is that the potential well in which the particle oscillates change its shape slowly as the particle drifts. It is shown how the Kruskal (1962) procedure can be applied to obtain expressions for the invariant and for drift velocities that are asymptotic in a smallness parameter epsilon. The procedure is illustrated by obtaining the invariant and drift velocities for particles traversing a perpendicular shock, and the generalized invariant is compared with the magnetic moment, and the drift orbits with the actual orbits, for a particular case. In contrast to the magnetic moment, the generalized first invariant is better for large gyroradii (large kinetic energies) than for small gyroradii. Expressions for the invariant when an electrostatic potential jump is imposed across the perpendicular shock, and when the particle traverses a rotational shear layer with a small normal component of the magnetic field are given.

Whipple, E. C.↗

A study of the topology of dissipating motions in direct numerical simulations of time-developing compressible and incompressible mixing layers

A preliminary investigation of the geometry of flow patterns in numerically simulated compressible and incompressible mixing layers was carried out using 3-D critical point methodology. Motions characterized by high rates of kinetic energy dissipation and/or high enstrophy were of particular interest. In the approach the partial derivatives of the velocity field are determined at every point in the flow. These are used to construct the invariants of the velocity gradient tensor and the rate-of-strain tensor (P, Q, R, and P(sub s), Q(sub s), R(sub s) respectively). For incompressible flow the first invariant is zero. For the conditions of the compressible simulation, the first invariant is found to be everywhere small, relative to the second and third invariants, and so in both cases the local topology at a point is mainly determined by the second and third invariants. The data at every grid point is used to construct scatter plots of Q versus R and Q(sub s) versus R(sub s). Most points map to a cluster near the origin in Q-R space. However, fine scale motions, that is motions which are characterized by velocity derivatives which scale with the square root of R(sub delta), tend to map to regions which lie far from the origin. Definite trends are observed for motions characterized by high enstrophy and/or high dissipation. The observed trends suggest that, for these motions, the second and third invariants of the velocity gradient and rate-of-strain tensors are strongly correlated. The second and third invariants of the rate-of-strain tensor are related by K(-Q(sub s))(exp 3/2), which is consistent with the above scaling of velocity derivatives. The quantity K appears to depend on Reynolds number with an upper limit K = 2(the square root of 3)/9 corresponding to locally axisymmetric flow. For both the compressible and incompressible mixing layer, regions corresponding to high rates of dissipation are found to be characterized by comparable magnitudes of R(sub ij)R(sub ij) and S(sub ij)S(sub ij). For the incompressible mixing layer, regions characterized by the highest values of enstrophy are found to have relatively low strain rates.

Chen, J. H.↗

Data-Driven Operator Theoretic Methods for Phase Space Learning and Analysis

This paper uses data-driven operator theoretic approaches to explore the global phase space of a dynamical system. In this work, we defined conditions for discovering new invariant subspaces in the state space of a dynamical system starting from an invariant subspace based on the spectral properties of the Koopman operator. When the system evolution is known locally in several invariant subspaces in the state space of a dynamical system, a phase space stitching result is derived that yields the global Koopman operator. Additionally, in the case of equivariant systems, a phase space stitching result is developed to identify the global Koopman operator using the symmetry properties between the invariant subspaces of the dynamical system and time-series data from any one of the invariant subspaces. Finally, these results are extended to topologically conjugate dynamical systems; in particular, the relation between the Koopman tuple of topologically conjugate systems is established. The proposed results are demonstrated on several second-order nonlinear dynamical systems including a bistable toggle switch. Our method elucidates a strategy for designing discovery experiments: experiment execution can be done in many steps, and models from different invariant subspaces can be combined to approximate the global Koopman operator.

42 ENGINEERING↗

The kinetic origin of the fluid helicity—A symmetry in the kinetic phase space

Helicity, a topological degree that measures the winding and linking of vortex lines, is preserved by ideal (barotropic) fluid dynamics. In the context of the Hamiltonian description, the helicity is a Casimir invariant characterizing a foliation of the associated Poisson manifold. Casimir invariants are special invariants that depend on the Poisson bracket, not on the particular choice of the Hamiltonian. The total mass (or particle number) is another Casimir invariant, whose invariance guarantees the mass (particle) conservation (independent of any specific choice of the Hamiltonian). In a kinetic description (e.g., that of the Vlasov equation), the helicity is no longer an invariant (although the total mass remains a Casimir of the Vlasov’s Poisson algebra). The implication is that some “kinetic effect” can violate the constancy of the helicity. To elucidate how the helicity constraint emerges or submerges, we examine the fluid reduction of the Vlasov system; the fluid (macroscopic) system is a “sub-algebra” of the kinetic (microscopic) Vlasov system. In the Vlasov system, the helicity can be conserved if a special helicity symmetry condition holds. To put it another way, breaking helicity symmetry induces a change in the helicity. We delineate the geometrical meaning of helicity symmetry and show that for a special class of flows (the so-called epi-two-dimensional flows), the helicity symmetry is written as ∂γ = 0 for a coordinate γ of the configuration space.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Gauge Picture of Quantum Dynamics

Although local Hamiltonians exhibit local time dynamics, this locality is not explicit in the Schrödinger picture in the sense that the wavefunction amplitudes do not obey a local equation of motion. We show that geometric locality can be achieved explicitly in the equations of motion by "gauging" the global unitary invariance of quantum mechanics into a local gauge invariance. That is, expectation values $\langle \psi|A|\psi \rangle$ are invariant under a global unitary transformation acting on the wavefunction $|\psi\rangle \to U |\psi\rangle$ and operators $A \to U A U^\dagger$, and we show that it is possible to gauge this global invariance into a local gauge invariance. To do this, we replace the wavefunction with a collection of local wavefunctions $|\psi_J\rangle$, one for each patch of space $J$. The collection of spatial patches is chosen to cover the space; e.g. we could choose the patches to be single qubits or nearest-neighbor sites on a lattice. Local wavefunctions associated with neighboring pairs of spatial patches $I$ and $J$ are related to each other by dynamical unitary transformations $U_{IJ}$. The local wavefunctions are local in the sense that their dynamics are local. That is, the equations of motion for the local wavefunctions $|\psi_J\rangle$ and connections $U_{IJ}$ are explicitly local in space and only depend on nearby Hamiltonian terms. (The local wavefunctions are many-body wavefunctions and have the same Hilbert space dimension as the usual wavefunction.) We call this picture of quantum dynamics the gauge picture since it exhibits a local gauge invariance. The local dynamics of a single spatial patch is related to the interaction picture, where the interaction Hamiltonian consists of only nearby Hamiltonian terms. We can also generalize the explicit locality to include locality in local charge and energy densities.

Physics↗

Self-organization via active exploration in robotic applications

We describe a neural network based robotic system. Unlike traditional robotic systems, our approach focussed on non-stationary problems. We indicate that self-organization capability is necessary for any system to operate successfully in a non-stationary environment. We suggest that self-organization should be based on an active exploration process. We investigated neural architectures having novelty sensitivity, selective attention, reinforcement learning, habit formation, flexible criteria categorization properties and analyzed the resulting behavior (consisting of an intelligent initiation of exploration) by computer simulations. While various computer vision researchers acknowledged recently the importance of active processes (Swain and Stricker, 1991), the proposed approaches within the new framework still suffer from a lack of self-organization (Aloimonos and Bandyopadhyay, 1987; Bajcsy, 1988). A self-organizing, neural network based robot (MAVIN) has been recently proposed (Baloch and Waxman, 1991). This robot has the capability of position, size rotation invariant pattern categorization, recognition and pavlovian conditioning. Our robot does not have initially invariant processing properties. The reason for this is the emphasis we put on active exploration. We maintain the point of view that such invariant properties emerge from an internalization of exploratory sensory-motor activity. Rather than coding the equilibria of such mental capabilities, we are seeking to capture its dynamics to understand on the one hand how the emergence of such invariances is possible and on the other hand the dynamics that lead to these invariances. The second point is crucial for an adaptive robot to acquire new invariances in non-stationary environments, as demonstrated by the inverting glass experiments of Helmholtz. We will introduce Pavlovian conditioning circuits in our future work for the precise objective of achieving the generation, coordination, and internalization of sequence of actions.

Ogmen, H.↗

The Statistical Mechanics of Ideal Homogeneous Turbulence

Plasmas, such as those found in the space environment or in plasma confinement devices, are often modeled as electrically conducting fluids. When fluids and plasmas are energetically stirred, regions of highly nonlinear, chaotic behavior known as turbulence arise. Understanding the fundamental nature of turbulence is a long-standing theoretical challenge. The present work describes a statistical theory concerning a certain class of nonlinear, finite dimensional, dynamical models of turbulence. These models arise when the partial differential equations describing incompressible, ideal (i.e., nondissipative) homogeneous fluid and magnetofluid (i.e., plasma) turbulence are Fourier transformed into a very large set of ordinary differential equations. These equations define a divergenceless flow in a high-dimensional phase space, which allows for the existence of a Liouville theorem, guaranteeing a distribution function based on constants of the motion (integral invariants). The novelty of these particular dynamical systems is that there are integral invariants other than the energy, and that some of these invariants behave like pseudoscalars under two of the discrete symmetry transformations of physics, parity, and charge conjugation. In this work the 'rugged invariants' of ideal homogeneous turbulence are shown to be the only significant scalar and pseudoscalar invariants. The discovery that pseudoscalar invariants cause symmetries of the original equations to be dynamically broken and induce a nonergodic structure on the associated phase space is the primary result presented here. Applicability of this result to dissipative turbulence is also discussed.

Shebalin, John V.↗

Non-Semisimple TQFT's and BPS q -Series

We propose and in some cases prove a precise relation between 3-manifold invariants associated with quantum groups at roots of unity and at generic q. Both types of invariants are labeled by extra data which plays an important role in the proposed relation. Bridging the two sides - which until recently were developed independently, using very different methods - opens many new avenues. In one direction, it allows to study (and perhaps even to formulate) q-series invariants labeled by spin c structures in terms of non-semisimple invariants. In the opposite direction, it offers new insights and perspectives on various elements of non-semisimple TQFT's, bringing the latter into one unifying framework with other invariants of knots and 3-manifolds that recently found realization in quantum field theory and in string theory.

3-manifold invariants↗

A Consistent AVHRR Visible Calibration Record Based on Multiple Methods Applicable for the NOAA Degrading Orbits: Validation - Part 2

Consistent cross-sensor Advanced Very High Resolution Radiometer (AVHRR) calibration coefficients are determined using desert, polar ice, and deep convective cloud (DCC) invariant Earth targets. The greatest AVHRR calibration challenge is the slow orbit degradation of the host satellite, which precesses toward a terminator orbit. This issue is solved by characterizing the invariant targets with NOAA-16 AVHRR observed radiances that have been referenced to the Aqua Moderate Resolution Imaging Spectrometer (MODIS) calibration using simultaneous nadir overpass (SNO) observations. Another benefit of the NOAA-16 invariant target-modeled reflectance method is that, because of the similarities among the AVHRR spectral response functions, a smaller spectral band adjustment factor is required than when establishing calibrations relative to a non-AVHRR reference instrument. The sensor- and band-specific calibration uncertainties, with respect to the calibration reference, are, on average, 2 percent and 3 percent for channels 1 and 2, respectively. The uncertainties are smaller for sensors that are in afternoon orbits, have longer records, and spend less time in terminator conditions. The multiple invariant targets referenced to Aqua MODIS (MITRAM) AVHRR calibration coefficients are evaluated for individual target consistency, compared against Aqua MODIS/AVHRR SNOs, and selected published calibration gains. The MITRAM and SNO relative calibration biases mostly agree to within 1 percent for channels 1 and 2, respectively. The individual invariant target and MITRAM sensor relative calibration biases are mostly consistent to within 1 percent and 2 percent for channels 1 and 2, respectively. The differences between the MITRAM and other published calibrations are mostly attributed to the reference instrument calibration differences.

Doelling, David R.↗

Fast and Scalable FFT-Based GPU-Accelerated Algorithms for Block-Triangular Toeplitz Matrices with Application to Linear Inverse Problems Governed by Autonomous Dynamical Systems

In this work, we present an efficient and scalable algorithm for performing matrix-vector multiplications (matvecs) for block Toeplitz matrices. Such matrices, which are shift-invariant with respect to their blocks, arise in the context of solving inverse problems governed by autonomous systems, and time-invariant systems in particular. In this article, we consider inverse problems that infer unknown parameters from observational data of a linear time-invariant dynamical system given in the form of partial differential equations (PDEs). Matrix-free Newton-conjugate-gradient methods are often the gold standard for solving these inverse problems, but they require numerous actions of the Hessian on a vector. Matrix-free adjoint-based Hessian matvecs require solution of a pair of linearized forward/adjoint PDE solves per Hessian action, which may be prohibitive for large-scale inverse problems. Time invariance of the forward PDE problem leads to a block Toeplitz structure of the discretized parameter-to-observable (p2o) map defining the mapping from inputs (parameters) to outputs (observables) of the PDEs. This block Toeplitz structure enables us to exploit two key properties: (1) compact storage of the p2o map and its adjoint, and (2) efficient fast Fourier transform–based Hessian matvecs. The proposed algorithm is mapped onto large multi-GPU clusters and achieves more than 80% of peak bandwidth on NVIDIA A100 GPUs. Excellent weak scaling is shown for up to 48 A100 GPUs. For the targeted problems, the implementation executes Hessian matvecs within fractions of a second, which is orders of magnitude faster than can be achieved by conventional matrix-free Hessian matvecs via forward/adjoint PDE solves.

97 MATHEMATICS AND COMPUTING↗

Properties of the new $ \mathcal{N} $ = 1 AdS 4 vacuum of maximal supergravity

The recent comprehensive numerical study of critical points of the scalar potential of four-dimensional N = 8, SO(8) gauged supergravity using Machine Learning software in [1] has led to a discovery of a new N = 1 vacuum with a triality-invariant SO(3) symmetry. Guided by the numerical data for that point, we obtain a consistent SO(3) × Z 2 -invariant truncation of the N = 8 theory to an N = 1 supergravity with three chiral multiplets. Critical points of the truncated scalar potential include both the N = 1 point as well as two new non-supersymmetric and perturbatively unstable points not found by previous searches. Studying the structure of the submanifold of SO(3) × Z 2 -invariant supergravity scalars, we find that it has a simple interpretation as a submanifold of the 14-dimensional Z $3\atop{2}$ -invariant scalar manifold (SU(1, 1)/U(1)) 7 , for which we find a rather remarkable superpotential whose structure matches the single bit error correcting (7, 4) Hamming code. This 14-dimensional scalar manifold contains approximately one quarter of the known critical points. We also show that there exists a smooth supersymmetric domain wall which interpolates between the new N = 1 AdS4 solution and the maximally supersymmetric AdS4 vacuum. Using holography, this result indicates the existence of an N = 1 RG flow from the ABJM SCFT to a new strongly interacting conformal fixed point in the IR.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Towards the Gravituhedron: new expressions for NMHV gravity amplitudes

In this paper, we present new expressions for n-point NMHV tree-level gravity amplitudes. We introduce a method of factorization diagrams which is a simple graphical representation of R-invariants in Yang-Mills theory. We define the gravity analogues which we call G-invariants, and expand the NMHV gravity amplitudes in terms of these objects. We provide explicit formulas of NMHV gravity amplitudes up to eight points in terms of G-invariants, and give the general definition for any number of points. We discuss the connection to BCFW representation, special behavior under large momentum shift, the role of momentum twistors and the intricate web of spurious poles cancelation. Because of the close connection between R-invariants and the (tree-level) Amplituhedron for Yang-Mills amplitudes, we speculate that the new expansion for gravity amplitudes should correspond to the triangulation of the putative Gravituhedron geometry.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗