Invariants of tensors

Last updated

In mathematics, in the fields of multilinear algebra and representation theory, the principal invariants of the second rank tensor are the coefficients of the characteristic polynomial [1]

Contents

,

where is the identity operator and are the roots of the polynomial and the eigenvalues of .

More broadly,any scalar-valued function is an invariant of if and only if for all orthogonal . This means that a formula expressing an invariant in terms of components, , will give the same result for all Cartesian bases. For example, even though individual diagonal components of will change with a change in basis, the sum of diagonal components will not change.

Properties

The principal invariants do not change with rotations of the coordinate system (they are objective, or in more modern terminology, satisfy the principle of material frame-indifference) and any function of the principal invariants is also objective.

Calculation of the invariants of rank two tensors

In a majority of engineering applications, the principal invariants of (rank two) tensors of dimension three are sought, such as those for the right Cauchy-Green deformation tensor which has the eigenvalues , , and . Where , , and are the principal stretches, i.e. the eigenvalues of .

Principal invariants

For such tensors, the principal invariants are given by:

For symmetric tensors, these definitions are reduced. [2]

The correspondence between the principal invariants and the characteristic polynomial of a tensor, in tandem with the Cayley–Hamilton theorem reveals that

where is the second-order identity tensor.

Main invariants

In addition to the principal invariants listed above, it is also possible to introduce the notion of main invariants [3] [4]

which are functions of the principal invariants above. These are the coefficients of the characteristic polynomial of the deviator , such that it is traceless. The separation of a tensor into a component that is a multiple of the identity and a traceless component is standard in hydrodynamics, where the former is called isotropic, providing the modified pressure, and the latter is called deviatoric, providing shear effects.

Mixed invariants

Furthermore, mixed invariants between pairs of rank two tensors may also be defined. [4]

Calculation of the invariants of order two tensors of higher dimension

These may be extracted by evaluating the characteristic polynomial directly, using the Faddeev-LeVerrier algorithm for example.

Calculation of the invariants of higher order tensors

The invariants of rank three, four, and higher order tensors may also be determined. [5]

Engineering applications

A scalar function that depends entirely on the principal invariants of a tensor is objective, i.e., independent of rotations of the coordinate system. This property is commonly used in formulating closed-form expressions for the strain energy density, or Helmholtz free energy, of a nonlinear material possessing isotropic symmetry. [6]

This technique was first introduced into isotropic turbulence by Howard P. Robertson in 1940 where he was able to derive Kármán–Howarth equation from the invariant principle. [7] George Batchelor and Subrahmanyan Chandrasekhar exploited this technique and developed an extended treatment for axisymmetric turbulence. [8] [9] [10]

Invariants of non-symmetric tensors

A real tensor in 3D (i.e., one with a 3x3 component matrix) has as many as six independent invariants, three being the invariants of its symmetric part and three characterizing the orientation of the axial vector of the skew-symmetric part relative to the principal directions of the symmetric part. For example, if the Cartesian components of are

the first step would be to evaluate the axial vector associated with the skew-symmetric part. Specifically, the axial vector has components

The next step finds the principal values of the symmetric part of . Even though the eigenvalues of a real non-symmetric tensor might be complex, the eigenvalues of its symmetric part will always be real and therefore can be ordered from largest to smallest. The corresponding orthonormal principal basis directions can be assigned senses to ensure that the axial vector points within the first octant. With respect to that special basis, the components of are

The first three invariants of are the diagonal components of this matrix: (equal to the ordered principal values of the tensor's symmetric part). The remaining three invariants are the axial vector's components in this basis: . Note: the magnitude of the axial vector, , is the sole invariant of the skew part of , whereas these distinct three invariants characterize (in a sense) "alignment" between the symmetric and skew parts of . Incidentally, it is a myth that a tensor is positive definite if its eigenvalues are positive. Instead, it is positive definite if and only if the eigenvalues of its symmetric part are positive.

See also

Related Research Articles

In mathematics, particularly in linear algebra, a skew-symmetricmatrix is a square matrix whose transpose equals its negative. That is, it satisfies the condition

In linear algebra, the characteristic polynomial of a square matrix is a polynomial which is invariant under matrix similarity and has the eigenvalues as roots. It has the determinant and the trace of the matrix among its coefficients. The characteristic polynomial of an endomorphism of a finite-dimensional vector space is the characteristic polynomial of the matrix of that endomorphism over any base. The characteristic equation, also known as the determinantal equation, is the equation obtained by equating the characteristic polynomial to zero.

In linear algebra, a square matrix  is called diagonalizable or non-defective if it is similar to a diagonal matrix. That is, if there exists an invertible matrix  and a diagonal matrix such that . This is equivalent to . This property exists for any linear map: for a finite-dimensional vector space , a linear map  is called diagonalizable if there exists an ordered basis of  consisting of eigenvectors of . These definitions are equivalent: if  has a matrix representation as above, then the column vectors of  form a basis consisting of eigenvectors of , and the diagonal entries of  are the corresponding eigenvalues of ; with respect to this eigenvector basis,  is represented by .

<span class="mw-page-title-main">Jordan normal form</span> Form of a matrix indicating its eigenvalues and their algebraic multiplicities

In linear algebra, a Jordan normal form, also known as a Jordan canonical form, is an upper triangular matrix of a particular form called a Jordan matrix representing a linear operator on a finite-dimensional vector space with respect to some basis. Such a matrix has each non-zero off-diagonal entry equal to 1, immediately above the main diagonal, and with identical diagonal entries to the left and below them.

In mathematics, the Hessian matrix, Hessian or Hesse matrix is a square matrix of second-order partial derivatives of a scalar-valued function, or scalar field. It describes the local curvature of a function of many variables. The Hessian matrix was developed in the 19th century by the German mathematician Ludwig Otto Hesse and later named after him. Hesse originally used the term "functional determinants". The Hessian is sometimes denoted by H or, ambiguously, by ∇2.

In mathematics, a Casimir element is a distinguished element of the center of the universal enveloping algebra of a Lie algebra. A prototypical example is the squared angular momentum operator, which is a Casimir element of the three-dimensional rotation group.

In numerical analysis, one of the most important problems is designing efficient and stable algorithms for finding the eigenvalues of a matrix. These eigenvalue algorithms may also find eigenvectors.

In linear algebra, a generalized eigenvector of an matrix is a vector which satisfies certain criteria which are more relaxed than those for an (ordinary) eigenvector.

In linear algebra, an eigenvector or characteristic vector is a vector that has its direction unchanged by a given linear transformation. More precisely, an eigenvector, , of a linear transformation, , is scaled by a constant factor, , when the linear transformation is applied to it: . It is often important to know these vectors in linear algebra. The corresponding eigenvalue, characteristic value, or characteristic root is the multiplying factor .

<span class="mw-page-title-main">Cauchy stress tensor</span> Representation of mechanical stress at every point within a deformed 3D object

In continuum mechanics, the Cauchy stress tensor, also called true stress tensor or simply stress tensor, completely defines the state of stress at a point inside a material in the deformed state, placement, or configuration. The second order tensor consists of nine components and relates a unit-length direction vector e to the traction vectorT(e) across an imaginary surface perpendicular to e:

Linear dynamical systems are dynamical systems whose evolution functions are linear. While dynamical systems, in general, do not have closed-form solutions, linear dynamical systems can be solved exactly, and they have a rich set of mathematical properties. Linear systems can also be used to understand the qualitative behavior of general dynamical systems, by calculating the equilibrium points of the system and approximating it as a linear system around each such point.

In mathematics, in particular in algebra, polarization is a technique for expressing a homogeneous polynomial in a simpler fashion by adjoining more variables. Specifically, given a homogeneous polynomial, polarization produces a unique symmetric multilinear form from which the original polynomial can be recovered by evaluating along a certain diagonal.

In mathematics, the Bianchi classification provides a list of all real 3-dimensional Lie algebras. The classification contains 11 classes, 9 of which contain a single Lie algebra and two of which contain a continuum-sized family of Lie algebras. The classification is important in geometry and physics, because the associated Lie groups serve as symmetry groups of 3-dimensional Riemannian manifolds. It is named for Luigi Bianchi, who worked it out in 1898.

In linear algebra, eigendecomposition is the factorization of a matrix into a canonical form, whereby the matrix is represented in terms of its eigenvalues and eigenvectors. Only diagonalizable matrices can be factorized in this way. When the matrix being factorized is a normal or real symmetric matrix, the decomposition is called "spectral decomposition", derived from the spectral theorem.

In statistics, principal component regression (PCR) is a regression analysis technique that is based on principal component analysis (PCA). More specifically, PCR is used for estimating the unknown regression coefficients in a standard linear regression model.

<span class="mw-page-title-main">Yeoh hyperelastic model</span> Phenomenological model of elastic materials

The Yeoh hyperelastic material model is a phenomenological model for the deformation of nearly incompressible, nonlinear elastic materials such as rubber. The model is based on Ronald Rivlin's observation that the elastic properties of rubber may be described using a strain energy density function which is a power series in the strain invariants of the Cauchy-Green deformation tensors. The Yeoh model for incompressible rubber is a function only of . For compressible rubbers, a dependence on is added on. Since a polynomial form of the strain energy density function is used but all the three invariants of the left Cauchy-Green deformation tensor are not, the Yeoh model is also called the reduced polynomial model.

In mathematics, the Hilbert–Mumford criterion, introduced by David Hilbert and David Mumford, characterizes the semistable and stable points of a group action on a vector space in terms of eigenvalues of 1-parameter subgroups.

In mathematics, particularly in linear algebra and applications, matrix analysis is the study of matrices and their algebraic properties. Some particular topics out of many include; operations defined on matrices, functions of matrices, and the eigenvalues of matrices.

The Lambda2 method, or Lambda2 vortex criterion, is a vortex core line detection algorithm that can adequately identify vortices from a three-dimensional fluid velocity field. The Lambda2 method is Galilean invariant, which means it produces the same results when a uniform velocity field is added to the existing velocity field or when the field is translated.

The Batchelor–Chandrasekhar equation is the evolution equation for the scalar functions, defining the two-point velocity correlation tensor of a homogeneous axisymmetric turbulence, named after George Batchelor and Subrahmanyan Chandrasekhar. They developed the theory of homogeneous axisymmetric turbulence based on Howard P. Robertson's work on isotropic turbulence using an invariant principle. This equation is an extension of Kármán–Howarth equation from isotropic to axisymmetric turbulence.

References

  1. Spencer, A. J. M. (1980). Continuum Mechanics. Longman. ISBN   0-582-44282-6.
  2. Kelly, PA. "Lecture Notes: An introduction to Solid Mechanics" (PDF). Retrieved 27 May 2018.
  3. Kindlmann, G. "Tensor Invariants and their Gradients" (PDF). Retrieved 24 Jan 2019.
  4. 1 2 Schröder, Jörg; Neff, Patrizio (2010). Poly-, Quasi- and Rank-One Convexity in Applied Mechanics. Springer.
  5. Betten, J. (1987). "Irreducible Invariants of Fourth-Order Tensors". Mathematical Modelling. 8: 29–33. doi: 10.1016/0270-0255(87)90535-5 .
  6. Ogden, R. W. (1984). Non-Linear Elastic Deformations. Dover.
  7. Robertson, H. P. (1940). "The Invariant Theory of Isotropic Turbulence". Mathematical Proceedings of the Cambridge Philosophical Society. 36 (2). Cambridge University Press: 209–223. Bibcode:1940PCPS...36..209R. doi:10.1017/S0305004100017199. S2CID   122767772.
  8. Batchelor, G. K. (1946). "The Theory of Axisymmetric Turbulence". Proc. R. Soc. Lond. A. 186 (1007): 480–502. Bibcode:1946RSPSA.186..480B. doi: 10.1098/rspa.1946.0060 .
  9. Chandrasekhar, S. (1950). "The Theory of Axisymmetric Turbulence". Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences. 242 (855): 557–577. Bibcode:1950RSPTA.242..557C. doi:10.1098/rsta.1950.0010. S2CID   123358727.
  10. Chandrasekhar, S. (1950). "The Decay of Axisymmetric Turbulence". Proc. R. Soc. A. 203 (1074): 358–364. Bibcode:1950RSPSA.203..358C. doi:10.1098/rspa.1950.0143. S2CID   121178989.