Strain-rate tensor

A two-dimensional flow that, at the highlighted point, has only a strain rate component, with no mean velocity or rotational component.

In continuum mechanics, the strain-rate tensor or rate-of-strain tensor is a physical quantity that describes the rate of change of the strain (i.e., the relative deformation) of a material in the neighborhood of a certain point, at a certain moment of time. It can be defined as the derivative of the strain tensor with respect to time, or as the symmetric component of the Jacobian matrix (derivative with respect to position) of the flow velocity. In fluid mechanics it also can be described as the velocity gradient, a measure of how the velocity of a fluid changes between different points within the fluid. Though the term can refer to a velocity profile (variation in velocity across layers of flow in a pipe), it is often used to mean the gradient of a flow's velocity with respect to its coordinates. The concept has implications in a variety of areas of physics and engineering, including magnetohydrodynamics, mining and water treatment.

The strain rate tensor is a purely kinematic concept that describes the macroscopic motion of the material. Therefore, it does not depend on the nature of the material, or on the forces and stresses that may be acting on it; and it applies to any continuous medium, whether solid, liquid or gas.

On the other hand, for any fluid except superfluids, any gradual change in its deformation (i.e. a non-zero strain rate tensor) gives rise to viscous forces in its interior, due to friction between adjacent fluid elements, that tend to oppose that change. At any point in the fluid, these stresses can be described by a viscous stress tensor that is, almost always, completely determined by the strain rate tensor and by certain intrinsic properties of the fluid at that point. Viscous stress also occur in solids, in addition to the elastic stress observed in static deformation; when it is too large to be ignored, the material is said to be viscoelastic.

Dimensional analysis

By performing dimensional analysis, the dimensions of velocity gradient can be determined. The dimensions of velocity are ${𝖫}^{\mathrm{𝟣}}{𝖳}^{\mathsf{-}\mathrm{𝟣}}$, and the dimensions of distance are ${𝖫}^{\mathrm{𝟣}}$. Since the velocity gradient can be expressed as $\frac{\mathrm{\Delta }\text{velocity}}{\mathrm{\Delta }\text{distance}}$. Therefore, the velocity gradient has the same dimensions as this ratio, i.e., ${𝖳}^{\mathsf{-}\mathrm{𝟣}}$.

In continuum mechanics

In 3 dimensions, the gradient $\mathrm{\nabla }𝐯$ of the velocity $𝐯$ is a second-order tensor which can be expressed as the matrix $𝐋$: $$𝐋=\mathrm{\nabla }𝐯=\left[\begin{array}{ccc}\frac{\partial v_x}{\partial x}& \frac{\partial v_y}{\partial x}& \frac{\partial v_z}{\partial x}\\ \frac{\partial v_x}{\partial y}& \frac{\partial v_y}{\partial y}& \frac{\partial v_z}{\partial y}\\ \frac{\partial v_x}{\partial z}& \frac{\partial v_y}{\partial z}& \frac{\partial v_z}{\partial z}\end{array}\right]$$ $𝐋$ can be decomposed into the sum of a symmetric matrix $\mathbf{\text{E}}$ and a skew-symmetric matrix $\mathbf{\text{W}}$ as follows $$\begin{array}{cc}𝐄& =\frac{1}{2}(𝐋+{𝐋}^{\mathsf{\text{T}}})\\ 𝐖& =\frac{1}{2}(𝐋-{𝐋}^{\mathsf{\text{T}}})\end{array}$$ $\mathbf{\text{E}}$ is called the strain rate tensor and describes the rate of stretching and shearing. $\mathbf{\text{W}}$ is called the spin tensor and describes the rate of rotation.

Relationship between shear stress and the velocity field

Sir Isaac Newton proposed that shear stress is directly proportional to the velocity gradient: $$\tau =\mu \frac{\partial u}{\partial y}.$$

The constant of proportionality, $\mu$, is called the dynamic viscosity.

Formal definition

Consider a material body, solid or fluid, that is flowing and/or moving in space. Let Template:Math be the velocity field within the body; that is, a smooth function from Template:Math such that Template:Math is the macroscopic velocity of the material that is passing through the point Template:Math at time Template:Mvar.

The velocity Template:Math at a point displaced from Template:Math by a small vector Template:Math can be written as a Taylor series: $$𝐯(𝐩+𝐫,t)=𝐯(𝐩,t)+(\mathrm{\nabla }𝐯)(𝐩,t)(𝐫)+\text{higher order terms},$$ where Template:Math the gradient of the velocity field, understood as a linear map that takes a displacement vector Template:Math to the corresponding change in the velocity.

In an arbitrary reference frame, Template:Math is related to the Jacobian matrix of the field, namely in 3 dimensions it is the 3 × 3 matrix $${\left(\mathrm{\nabla }𝐯\right)}^{\mathrm{T}}=\left[\begin{array}{ccc}{\partial }_1{v}_1& {\partial }_2{v}_1& {\partial }_3{v}_1\\ {\partial }_1{v}_2& {\partial }_2{v}_2& {\partial }_3{v}_2\\ {\partial }_1{v}_3& {\partial }_2{v}_3& {\partial }_3{v}_3\end{array}\right]=𝐉.$$ where Template:Mvar is the component of Template:Math parallel to axis Template:Mvar and Template:Math denotes the partial derivative of a function Template:Mvar with respect to the space coordinate Template:Mvar. Note that Template:Math is a function of Template:Math and Template:Mvar.

In this coordinate system, the Taylor approximation for the velocity near Template:Math is $$v_i(𝐩+𝐫,t)=v_i(𝐩,t)+\sum _j{J}_{ij}(𝐩,t)r_j=v_i(𝐩,t)+\sum _j{\partial }_j{v}_i(𝐩,t)r_j;$$ or simply $$𝐯(𝐩+𝐫,t)=𝐯(𝐩,t)+𝐉(𝐩,t)𝐫$$

if Template:Math and Template:Math are viewed as 3 × 1 matrices.

Symmetric and antisymmetric parts

Any matrix can be decomposed into the sum of a symmetric matrix and an antisymmetric matrix. Applying this to the Jacobian matrix with symmetric and antisymmetric components Template:Math and Template:Math respectively: $$\begin{array}{cccc}𝐄& =\frac{1}{2}(𝐉+{𝐉}^{\mathsf{\text{T}}})& 𝐑& =\frac{1}{2}(𝐉-{𝐉}^{\mathsf{\text{T}}})\\ E_{ij}& =\frac{1}{2}({\partial }_j{v}_i+{\partial }_i{v}_j)& R_{ij}& =\frac{1}{2}({\partial }_j{v}_i-{\partial }_i{v}_j)\end{array}$$

This decomposition is independent of coordinate system, and so has physical significance. Then the velocity field may be approximated as $$𝐯(𝐩+𝐫,t)\approx 𝐯(𝐩,t)+𝐄(𝐩,t)(𝐫)+𝐑(𝐩,t)(𝐫),$$ that is, $$\begin{array}{cc}{v}_i(𝐩+𝐫,t)& =v_i(𝐩,t)+\sum _j{E}_{ij}(𝐩,t)r_j+\sum _j{R}_{ij}(𝐩,t)r_j\\ & =v_i(𝐩,t)+\frac{1}{2}\sum _j({\partial }_j{v}_i(𝐩,t)+{\partial }_i{v}_j(𝐩,t))r_j+\frac{1}{2}\sum _j({\partial }_j{v}_i(𝐩,t)-{\partial }_i{v}_j(𝐩,t))r_j\end{array}$$

The antisymmetric term Template:Math represents a rigid-like rotation of the fluid about the point Template:Math. Its angular velocity Template:Math is $$\overrightarrow{\omega }=\frac{1}{2}\mathrm{\nabla }\times 𝐯=\frac{1}{2}\left[\begin{array}{c}{\partial }_2{v}_3-{\partial }_3{v}_2\\ {\partial }_3{v}_1-{\partial }_1{v}_3\\ {\partial }_1{v}_2-{\partial }_2{v}_1\end{array}\right].$$

The product Template:Math is called the vorticity of the vector field. A rigid rotation does not change the relative positions of the fluid elements, so the antisymmetric term Template:Math of the velocity gradient does not contribute to the rate of change of the deformation. The actual strain rate is therefore described by the symmetric Template:Math term, which is the strain rate tensor.

Shear rate and compression rate

Template:Multiple image The symmetric term Template:Math (the rate-of-strain tensor) can be broken down further as the sum of a scalar times the unit tensor, that represents a gradual isotropic expansion or contraction; and a traceless symmetric tensor which represents a gradual shearing deformation, with no change in volume: $$𝐄(𝐩,t)(𝐫)=𝐒(𝐩,t)(𝐫)+𝐃(𝐩,t)(𝐫).$$

That is, $$E_{ij}=\underset{\text{rate-of-expansion tensor }{S}_{ij}}{\underbrace{\frac{1}{3}\left(\sum _k{\partial }_k{v}_k\right){\delta }_{ij}}}+\underset{\text{rate-of-shear tensor }{D}_{ij}}{\underbrace{\stackrel{E_{ij}}{\overbrace{\frac{1}{2}({\partial }_i{v}_j+{\partial }_j{v}_i)}}-S_{ij}}},$$

Here Template:Math is the unit tensor, such that Template:Mvar is 1 if Template:Math and 0 if Template:Math. This decomposition is independent of the choice of coordinate system, and is therefore physically significant.

The trace of the expansion rate tensor is the divergence of the velocity field: $$\mathrm{\nabla }\cdot 𝐯={\partial }_1{v}_1+{\partial }_2{v}_2+{\partial }_3{v}_3;$$ which is the rate at which the volume of a fixed amount of fluid increases at that point.

The shear rate tensor is represented by a symmetric 3 × 3 matrix, and describes a flow that combines compression and expansion flows along three orthogonal axes, such that there is no change in volume. This type of flow occurs, for example, when a rubber strip is stretched by pulling at the ends, or when honey falls from a spoon as a smooth unbroken stream.

For a two-dimensional flow, the divergence of Template:Math has only two terms and quantifies the change in area rather than volume. The factor 1/3 in the expansion rate term should be replaced by Template:Sfrac in that case.

Examples

The study of velocity gradients is useful in analysing path dependent materials and in the subsequent study of stresses and strains; e.g., Plastic deformation of metals. The near-wall velocity gradient of the unburned reactants flowing from a tube is a key parameter for characterising flame stability.Template:Rp The velocity gradient of a plasma can define conditions for the solutions to fundamental equations in magnetohydrodynamics.

Fluid in a pipe

Consider the velocity field of a fluid flowing through a pipe. The layer of fluid in contact with the pipe tends to be at rest with respect to the pipe. This is called the no slip condition. If the velocity difference between fluid layers at the centre of the pipe and at the sides of the pipe is sufficiently small, then the fluid flow is observed in the form of continuous layers. This type of flow is called laminar flow.

The flow velocity difference between adjacent layers can be measured in terms of a velocity gradient, given by $\mathrm{\Delta }u/\mathrm{\Delta }y$. Where $\mathrm{\Delta }u$ is the difference in flow velocity between the two layers and $\mathrm{\Delta }y$ is the distance between the layers.

See also

• Stress tensor (disambiguation)
• Template:Slink, the spatial and material velocity gradient from continuum mechanics