Stretching and shearing contamination analysis for Liutex and other vortex identification methods

The newly developed vortex-identification method, Liutex, has provided a new systematic description of the local fluid rotation, which includes scalar, vector, and tensor forms. However, the advantages of Liutex over the other widely used vortex-identification methods such as Q, ∆, λ2, and λci have not been realized. These traditional methods count on shearing and stretching as a part of vortex strength. But, in the real flow, shearing and stretching do not contribute to fluid rotation. In this paper, the decomposition of the velocity gradient tensor is conducted in the Principal Coordinate for uniqueness. Then the contamination effects of stretching and shearing of the traditional methods are investigated and compared with the Liutex method in terms of mathematical analysis and numerical calculations. The results show that the Liutex method is the only method that is not affected by either stretching or shear, as it represents only the local fluid rigid rotation. These results provide supporting evidence that Liutex is the superior method over others.


Introduction
A vortex is recognized as the rotational motion of fluids. Within the last several decades, a lot of vortex identification methods have been developed to track the vortical structure in a fluid flow; however, we still lack unambiguous and universally accepted vortex identification criteria. This is the impediment that has caused a lot of confusions and misunderstandings in turbulence research [1]. In most research papers and textbooks, the vorticity tube/filament is regarded as a vortex and the magnitude of the vorticity as the local rotational strength. Many researchers widely acknowledged the concept of vortex defined as the vorticity concentration and other vorticity-based methods [2,3] as the vorticity vector was believed to offer a mathematical definition of fluid rotational motion. Zhou and Antonia (1993) [4] utilized the spatially phasedcorrelated vorticity to characterize large-scale and organized structures in the cylinder wake. However, problems arose while applying in viscous flows, particularly in turbulent flows. In turbulent viscous flows, the rigid rotation strength is minimal near the wall, where shear stress is dominant. Many researchers in the literature have supported this claim, indicating the inadequacies of vorticity-based methods. Epps (2017) [5] observed that vorticity could not distinguish a vortical region with rotational motions from a strong shear layer. Robinson (1991) [6] also uncovered the fact that the regions of strong vorticity and actual vortices are weakly related. These vorticity-based methods are the first generation (1G) of vortex identification methods as per [7], while scalar based schemes such as Δ, Q, λ 2 , λ ci , and Ω are known as second-generation (2G) criteria.
Attempting to remedy the problems of vorticity based methods, 2G vortex identification criteria such as Δ [8,9], Q [10], λ 2 [11], λ ci [12], Ω [13,14], etc., have been introduced and extensively applied in visualizing vortex structures over the past three decades. Although the 2G-schemes have been able to judge the presence of the local rotational motion to some extent, these methods need an adjustment of threshold to visualize the iso-surface plot accurately. This poses a problem as it is difficult to adjust the threshold that will define the boundary of vortical structures in a particular case [13,15]. Moreover, the Q and λ 2 methods are restrictive for incompressible flows due to their incompressibility assumption. These and other issues of 1G and 2G vortex identification methods prompted the development of a new vortex identification method, the Liutex method [16,17], also known as a third-generation vortex identification scheme.
According to Chong and Perry (1990) [8], if the velocity gradient tensor of a threedimensional flow field has a pair of complex conjugate eigenvalues and a real eigenvalue, then the instantaneous streamline pattern presents a local swirling motion.   [7] acknowledged that these streamlined patterns are the direct result of the real eigenvector of the velocity gradient tensor. In fact, this real eigenvector is the axis of fluid rotation which is also known as the direction of Liutex. Unlike 2G schemes, the Liutex method is a novel eigenvector-based method that is local, accurate, unique, and systematic. Furthermore, the systematical definition of Liutex is given as a vector which has its corresponding scalar and tensor forms [17]. The vector form of Liutex gives the direction of the local fluid rotation, while the magnitude of Liutex represents the rotational strength of a fluid rotation. Also, the Liutex vector is Galilean invariant [18][19][20]. According to Liu et al. [7], the Liutex method, Liutex-Omega method [16], Liutex Core Line method, and other Liutex-based methods [21][22][23] are the third generation (3G) of vortex identification methods. The Liutex method is discussed briefly in section 2.6.
After the introduction of the Liutex method in 2018, many doors of investigations of fluid rotations and related topics have opened. Based on Liutex definition, Gao and Liu proposed a velocity gradient tensor decomposition [24] in the rotated coordinate system, which was somewhat lengthy and cumbersome. So, later in the same year, they gave explicit expressions for the Rortex/Liutex tensor and velocity gradient tensor decomposition [25] without coordinate rotation. Recently, Yu and Liu proposed a new unique definition of Principal Coordinates and conducted a unique Principal Tensor decomposition under the Principal coordinates [26]. In this paper, the concept of Principal Coordinates and Principal Tensor decomposition is revisited briefly in section 3. Also, the issue of shear and stretching contamination of 1G and several 2G vortex identification methods is investigated in terms of mathematical analysis and numerical calculations. Our investigation involves calculating the velocity gradient tensor at a point inside the numerical domain (flow field). Then, a unique Principal Tensor decomposition is performed in the Principal Coordinates. In other words, the velocity gradient tensor is uniquely decomposed into its rigid rotation part, stretching part, and shearing part in the Principal Coordinates. From our investigation, it is found that 1G and several 2G vortex identification methods misinterpret shearing and stretching as a constituent of vortical structure, which is not true. Also, these methods have dimensional problems at different physical levels, which are discussed in section 4 in detail.
The organization of this paper is as follows: In Section 2, some traditional vortex identification methods are revisited including vorticity tube/filaments, Δ, Q, λ 2 and λ ci methods. Section 3 describes the Liutex-based Principal Tensor Decomposition in the unique Principal Coordinates. Based on the unique Principal Tensor Decomposition in the Principal Coordinates, the theoretical contamination analysis of the vorticity, Δ, Q, λ 2 , λ ci and Liutex schemes are examined and presented in Section 4. Numerical experiments and their results can be found in Sections 5 and 6, respectively, and some concluding remarks are made in the last section.
2 Review of 1G and 2G vortex identification methods

Vorticity-based method
Since Helmholtz introduced the concept of vorticity tube/filament in 1858 [2], many researchers have believed that vortices consist of small vorticity tubes which are called vortex filaments, and the magnitude of vorticity gives the vortex strength. The vorticity vector is mathematically defined as the curl of the velocity. i.e.
Helmholtz also defined vortex lines and vortex filaments based on the vorticity vector. According to Helmholtz, vortex lines are drawn through the fluid mass so that their direction at every point coincides with the direction of the momentary axis of rotation of the water particles lying on it. He also defined vortex filament as the portions of the fluid mass cut out from it by way of constructing corresponding vortex lines through all points of circumference of an infinitely small surface element. Another vorticity-based method to define a vortical structure is vortex-tube. Lamb [27] defines vortex-tube as: "If through every point of a small, closed curve we draw the corresponding vortex-line, we mark out a tube known as vortex-tube."

Δ method
The Δ method defines a vortex to be the region where the velocity gradient tensor ∇ v ! has a pair of complex conjugate eigenvalues and a real eigenvalue [8]. If λ 1 , λ 2 and λ 3 are the eigenvalues of the 3 ×3 matrix of ∇ v ! ; then the characteristic equation can be written as where, I 1 , I 2 , and I 3 are the first, second, and third principal invariants of the characteristic eq. (2) and given by The discriminant of the characteristic eq. (2) is given by whereQ For incompressible flow, the first principal invariant I 1 = 0, which results in eq. (6) be- there exists one real and two conjugate complex eigenvalues. The latter means that the point is inside a vortex region. Although the Δ-method can capture the vortex region successfully, it is susceptible to the threshold value, which is man-made and arbitrary in general.

Q method
Proposed by Hunt et al. (1988) [10], the Q method is one of the most popular methods used to visualize the vortex structure. According to Chakraborty et al. (2005) [28], Q is also the second principal invariant of the velocity gradient, and the value of Q can be calculated as half of the difference of squares of the Frobenius norm of the vorticity tensor and strain-rate tensor. i.e., where A and B are the symmetric and anti-symmetric part of the velocity gradient tensor, respectively. ð8Þ The Q method considers that a vortex occurs in the region where Q > 0. Using the Q method, it is easy to track the vortical structure by iso-surface plotting. However, Q is scalar-valued, and a proper threshold, which in general does not exist, is required to visualize the vortex region. Also, there exists an inconsistency between the Δ method and Q method as ðR 2 Þ 2 is always positive, and even if Q < 0, it is still possible that Δ could be positive.

λ ci criterion
The λ ci criterion [12,28] uses the imaginary part of the complex eigenvalues of the velocity gradient tensor to visualize the vortex structure. It is based on the idea that the local time-frozen streamlines exhibit a rotational flow pattern when ∇ v ! has a pair of complex conjugate eigenvalues. In this case, the tensor transformation of ∇ v ! is given by where λ r is the real eigenvalue with the corresponding real eigenvector v ! r and the pair of complex conjugate eigenvalues are λ cr ± iλ ci with corresponding eigenvectors v ! cr ± i v ! ci : In this case, in the local curvilinear system (c 1 , c 2 , c 3 ) spanned by the eigenvector the instantaneous streamlines exhibit a spiral motion. The equations of such streamlines can be written as: Despite being used by many researchers to track the vortical structure, the strength of this swirling motion may not be quantified accurately by λ ci . This is largely because the tensor transformation given by (10) is a similar transformation, but not an orthog- In a non-orthogonal system, vorticity is not invariant ( j ∇ Â v ! |≠2λ ci Þ: As well-known that vorticity is invariant in any orthogonal transformation system and therefore, λ ci cannot be used as the strength of fluid rotation. We cannot use a coordinate transformation which changes the vorticity, and the new tensor is different from the original tensor, or ∇V ! is not ∇ v ! ; making it clear that transformation used in λ ci method is not Galilean invariant. On the other hand, the Liutex transformation is an orthogonal transformation that preserves the invariance of vorticity.

λ 2 criterion
The λ 2 -criterion is calculated based on the observation that in a vortical region, pressure tends to be the lowest on the axis of a swirling motion of fluid particles. This happens because the centrifugal force is balanced by the pressure force (the cyclostrophic balance). This method is valid only in a steady inviscid and incompressible flow [11]. However, this assumption fails to accurately identify vortices under compressible, unsteady, and viscous conditions. By neglecting these unsteady and viscous effects, the symmetric part S of the gradient of the incompressible Navier-Stokes equation can be expressed as: where p is the pressure and eq. (14) is a representation of the pressure Hessian matrix.
To capture the region of local pressure minimum in a plane perpendicular to the vortex core line, Jeong & Hussain defined the vortex core as a connected region with two positive eigenvalues of the pressure Hessian matrix, i.e., a connected region with two negative eigenvalues of the symmetric tensor S. If λ S1 , λ S2 & λ S3 are three real eigenvalues of the symmetric tensor S and when setting them in order in such a way that λ S1 ≥ λ S2 ≥ λ S3 , there must be λ S2 < 0 as two eigenvalues are negative, which confirms the existence of a vortex. In general, λ S2 cannot be expressed in terms of eigenvalues of the velocity gradient tensor; however, in some special cases when eigenvectors are orthonormal, λ S2 can be exclusively determined by eigenvalues of the velocity gradient tensor. Vortex structure can be visualized as iso-surface by selecting a proper threshold of λ S2 [7]. The relation between the eigenvalues of the symmetric tensor A 2 + B 2 and second invariant Q is given by: While the Q-criterion measures the excess of vorticity rate over the strain rate magnitude in all directions, the λ 2 criterion looks for this excess only on a specific plane [11]. For comparative studies of different schemes, another article, "on the relationship between vortex identifications methods" by Chakraborty & Balachandar [28] is also recommended. There is yet another article by Liu et al., "Third generation of vortex identification methods", which reviews the second generation of vortex identification methods. These 2G methods are based on the eigenvalues of ∇ v ! . Afterward, we are presented with the novel third generation of vortex identification methods [7].

Liutex
Definition 1: Liutex is a vector defined as R ! ¼ R r ! , where R is the magnitude of Liutex, and r ! is the direction of Liutex.
According to Ref. [29], r ! is the normalized real eigenvector of the velocity gradient tensor such that ω ! • r ! > 0 and the explicit formula of R is Dr. Liu classified the Liutex method as the third generation of vortex identification method, while the Q method, λ ci criterion and λ 2 criterion are all classified as secondgeneration methods since they are all scalar and eigenvalue related. Liutex is a vector which overcomes the drawbacks of the scalar methods, e.g., it can give the direction of rotation and represents the rotational strength accurately. Liutex is finally found as a mathematical definition of fluid rigid rotation or vortex.

Principal coordinates and principal decomposition
where λ r and λ cr are the real eigenvalue and real part of the conjugate complex eigenvalue pair of the velocity gradient tensor, respectively, for rotation points.
3. ∂U ∂Y < 0 and j ∂U ∂Y j ≤j ∂V ∂X j. 4. The rotational angle of the X-Y coordinates around the Z-axis must be smaller than 90°or −90°< θ ≤ 90°T he following two theorems are valid in the Principal Coordinate.
Proof: For an arbitrary velocity gradient tensor ∇v, there exists a special orthogonal rotation matrix Q r that aligns the Z-axis of the new frame XYZ with the Liutex direction r ! after rotation.
Where xyz is the original frame and XYZ is the new frame after Q r -rotation. Then, we apply a second rotation P r around the new Z-axis to get the Principal Coordinates.
Note that if X-and Y-axes are rotated by additional 180°, the new coordinates are still satisfied with the principal coordinate definition. However, since we limit the Protation angle to be less than 90°or −90°< θ ≤ 90°. Under such a limit, the Principal Coordinate is unique.
Definition 3: The Principal Decomposition is the decomposition of velocity gradient tensor under the Principal Coordinates i.e.
Here A represents the rotation part, B represents the shear part, and C represents the stretching part.
Mathematically, the popular Cauchy-Stokes decomposition is correct, but the physical meaning is dubious. For instance, according to the Cauchy-Stokes decomposition, the velocity gradient tensor 0 2y . This indicates that shear flow has rotation. But, in general, the shear flow does not represent fluid rotation. Also, the Cauchy-Stokes decomposition depends on the choice of coordinates, and hence it is not Galilean invariant. Yu et al. [26] used the Cauchy-Stokes decomposition to decompose the velocity gradient tensor into the rotation, deformation, and stretching matrices and found out that the norm of rotation, symmetric deformation, and stretching matrices are all different, except for 2-norm, under different coordinates. Then, a question arose. Under which coordinates, does the Cauchy-Stokes tensor decomposition give the right stretching (compression) and deformation? In other words, is it possible to decompose a velocity gradient tensor into the rotation, deformation, and stretching matrices without being affected by a coordinate change? Yu et.al [26] presume that there must be a unique coordinate that can give a unique tensor decomposition for stretching (compression), deformation, shear, and rotation. They call these coordinates as the "Principal Coordinates" and tensor decomposition as the "Principal Tensor Decomposition". The Principal Decomposition, being Galilean invariant, correctly and uniquely decomposes the velocity gradient tensor into its rotation part, shear part, and stretching part. This decomposition expresses that the shear and stretching contamination analysis can only be uniquely done in the Principal Coordinates as decomposition in Principal Coordinates is not affected by coordinate changes. The detailed information about the Principal Coordinates and Principal Tensor Decomposition can be found in reference [26].

Theoretical contamination analysis
Kolar et al [32] numerically investigated the stretching response of several popular vortex identification methods and found that Liutex does not respond to uniaxial stretching coupled with an inevitable uniform radial contraction for incompressible flow. In this paper, we give mathematical expression that shows the effect of strectching and shearing on several vortex identification methods. Also, these effects are presented numerically.The velocity gradient tensor, ∇V nates, which shows the relation between Liutex and λ ci is Apparently, λ ci is contaminated by ε (shear). The vortex structure drawn by Liutex and λ ci are quite different and detailed analysis can be found in [17]. λ ci considers shear deformation ε as a part of the fluid rotation and thus may not represent the pure/rigid rotational strength of fluid motion correctly. A detailed explanation with mathematical expression about how shear affects λ ci is given later in chapter 4.3. λ ci is also characterized as a scalar-valued criterion [7]. The velocity gradient tensor in the Principal Coordinate is: Then the Principal Decomposition is given by A is the rotation part, B is the shear part and C is the stretching part. For the first time in the literature, we are defining the several vortex identification schemes in the Principal Coordinate, which will allow us to analyze the contamination effects on these traditional vortex identification methods.

Contamination of vorticity
Plugging in the gradient components from (40) into the eq. (1) gives us the vorticity vector as and its magnitude is From the eqs. (42) and (43), it can be concluded that a vorticity vector does not only represent rotation but also claims shearing and stretching components to be a part of the vortical structure, which is clearly contaminated by shears.

Contamination of Q method
First, we perform the Cauchy-Stokes decomposition on the Principal Coordinate. Then, the scalar magnitude of the Q method can be calculated based on the definition provided by eq. (7).
In the expression of Q above, there are not only R (the magnitude of rotation) but also ε, λ cr and λ r , which are components of either the shear part or stretching part. Therefore, the value of Q is undoubtedly contaminated by shear and stretching. Also, Q contains an R 2 term, indicating dimensional distortion of a fluid rotation.

Contamination of λ ci criterion
The characteristic equation of velocity gradient tensor given by (40) is Thus, the eigenvalues are.
Since rotation is orthogonal, eigenvalues are the same as the original velocity gradient tensor, Thus, Despite being dimensionally correct, the expression of λ ci has ε, which is a component of the shear part and thus is contaminated by shear.

Contamination of Δ method
In section 4.2, the three roots of the characteristic equation are Plug their values into (2),(3), and (4) Then, the expression for Δ can be written as The expression of Δ includes ε, λ r and λ cr , indicating the contamination of Δ method by shear and stretching. Moreover, Δ contains the sixth power of R, which confirms that it cannot accurately denote rotation strength.
Next, we investigate the contamination of these schemes with the help of a numerical example.

Vortex example: Burgers vortex
In this test case, Burgers vortex is examined to justify the comparison of the effects of shearing and stretching/compression on different criteria. The Burgers vortex is an exact steady solution of the Navier-Stokes equation that can be used to model fine scales of turbulence [30]. The Burgers vortex forms when an inward radial flow concentrates and spins around the symmetric axis and the flow moves out in both directions along the zaxis [31]. The velocity gradient tensor is obtained from the Burgers vortex.
The velocity components of Burgers vortex in the cylindrical coordinate system are given by: where Γ is the circulation that characterizes the intensity of the vortex and ξ is another constant representing the axisymmetric strain rate, ν the kinematic viscosity, and r is the distance of the chosen point from the centerline in the Burgers vortex. For post-processing, the velocity components are converted into the Cartesians coordinate system given by: The existence of the vortex structure is highly contingent on the selection of parameters. For the proper vortex structure visualization, we take dimensionless quantity ξ = 0.042, Γ = 1.45 and ν = 0.01. The calculation domain is limited to 50 ×20 × 20 grid points with a step size of 0.5 due to the limited computational memory. The vortex strength is strong in the core and becomes weak when moving away from the center. The streamlines of the Burgers vortex exhibit a spiral pattern around the vortex rotation axis line which is demonstrated in Fig. 1 a and b, which show that flows enter from a radial direction and stretches outward spinning around the axis.

Numerical contamination analysis
The velocity gradient tensor in the Principal Coordinate is: where u, v and w are the three components of the velocity along the x, y, and z directions in Cartesian Coordinate. Since our purpose is to demonstrate the shear and stretching contaminations of different criteria numerically, we add shear and stretching components separately and calculate how different criteria respond.

Adding shear components
The matrix corresponding to shear is in the form of: The subscript "a" refers to the "adding" of shearing components.
The new velocity gradient tensor in the Principal Coordinates after adding shear is Under the Principal Coordinates, the local rotation axis is the Z-axis, so, ξ a (value of ∂w ∂x component) and η a (value of ∂w ∂y component) are not in the rotation plane. Thus, these components will not influence the rotation strength. However, ε a is in the rotation plane, and thus affecting rotation strength. If ε a is in the direction of rotation, the Liutex magnitude remains unchanged, but if ε a is in the reverse direction of fluid rotation or negative and |ε a | > ε, it may change the Liutex magnitude. This change occurred not because Liutex was contaminated by shear but because j 1 2 R þ ε þ ε a j < R 2 which causes the strength of Liutex reduced according to the Liutex definition that the magnitude of rotation strength is minfj ∂u ∂y j θ ; j ∂v ∂x j θ g where j ∂u ∂y j θ and j ∂v ∂x j θ are the minimum absolute values of ∂u ∂y and ∂v ∂x respectively. Here we assume the rotation is anti-clockwise along the Z-axis. Therefore, when shear is in the opposite direction of rotation, minf j ∂u ∂y j θ ; j ∂v ∂x j θ g may change, affecting the Liutex magnitude. Following Proposition 1 gives the accurate measure about when shear may change the Liutex magnitude. Proof: If ε a ≥ − ε, then j 1 2 R þ ε þ ε a j ≥ j − 1 2 Rj, which shows that the Liutex magnitude is still R which is defined as twice the angular speed of the rigid rotation.
Remark: Physically, ε a < 0 means adding a shear against the rotation direction and ε a < − ε indicates that shear is strong to the extent that it can affect rotation strength. So, keeping the rotation strength the same, the shear matrix we add should satisfy the following requirement: ε a ≥ − ε. In other words, Proposition 1 just shows in some condition, external shear disturbance cannot change Liutex, but in other conditions (ε a < − ε), the external shear disturbance may. However, in any case, for a certain velocity gradient tensor, Liutex is always defined to represent the rigid rotation including the direction and strength. The definition has no shear and no stretching and cannot be considered as contaminated by shear or stretching.

Adding stretching components
The matrix corresponding to stretching is of the form The subscript "a" means "adding" the stretching components.
The new velocity gradient tensor in the Principal Coordinate after adding stretching is Proof: Based on theorem 1, the rotation strength of ∇V 2 is still R. Therefore, to keep the rotation strength unchanged, the added ∇V stretching should satisfy α a = β a .
Remark: Proposition 1 and Proposition 2 just give some conditions under which Liutex may or may not change the magnitude and direction. They do not mean Liutex could be contaminated by shear or stretching for a certain velocity gradient tensor. Liutex is always a rigid rotation vector and never contaminated by any shear and stretching. On the other hand, look eqs. (43), (45), (48) and (54), all vorticity, Q, λ ci , Δ are contaminated by shear or stretching or both where R is defined as twice the angular speed of rigid rotation of fluids. The external shear may change Liutex under some conditions which means the external shear can change the rotation strength but does not mean Liutex is contaminated by shear.

Stretching contamination analysis
The following procedure was implemented to depict the reaction of the vortex identification schemes over the change in stretching effects. Firstly, we select a point in Burgers vortex whose velocity gradient tensor is: Then for convenience, the stretching matrix ∇V stretching is added and the values of Q, Δ, λ 2 , λ ci , and Liutex are recorded where ∇V stretching is given by (note that it satisfies the requirement of adding stretching): The sum of diagonal elements of ∇V stretching is equal to zero, so it satisfies the continuity equation of incompressible flow. The stretching effect given by ∇V stretching is increased repeatedly, step by step, and the values of every scheme are recorded. The results are presented in the following graph where the x-axis represents the relative stretching rate (the ratio of stretching component and vorticity magnitude), and the y-axis gives the corresponding values of different vortex identification methods.
The change in stretching components results in a change in Q, Δ, and λ 2 . On the other hand, λ ci and L are not changed, indicating that λ ci and L are not contaminated by the stretching effect, as shown in Fig. 2. Further observation of Fig. 2 reveals that the increase in the stretching effect on Q decreases its value after a certain point. The Q values ultimately become negative, indicating a non-vortical structure. However, other criteria indicate there is a vortex structure. So, it can be concluded that Q may conflict with other criteria and violate the Δ-rule. Theoretically, this is understandable as in Eq. (6), Δ could be positive even if Q is negative provided that the square of the third variant is large enough. The computational results shown in Fig. 2 coincide with the theoretical analysis given by Eqs. (45), (48) and (54).

Shear contamination analysis
A similar procedure to the stretching effect is implemented for the graphic representation of the shearing effect on the vortex identification schemes. Again, for our convenience, the shearing matrix ∇V shear is defined as: ∇V shear is added to ∇V ! B a few times, and then the corresponding values of the different criteria are recorded. The results are presented in the following graph, where the x-axis represents the relative shear rate, and the y-axis gives the corresponding values of different vortex identification methods.
The relative shear rate is the ratio of the shearing component over vorticity magnitude. More precisely, it is the ratio of shear over vorticity at any point in our numerical domain. The reason we are dividing shear by vorticity is to refrain shear from getting too large. Fig. 3 indicates that Q, Δ, λ 2 , and λ ci are all affected by the change in the shearing component, while Liutex has remained unaffected by shear. It can be concluded from Fig. 2 and Fig. 3 that Q, Δ, λ 2 , and λ ci are affected by  either shear, stretching, or both at different levels, whereas the Liutex method is unaffected by both.
Since the values of Δ are very small, Δ looks to be consistently zero because Δ has sixth power of R which is small; however, if we use zoom in, it can be observed that the Δ values are increasing. The Δ -plots are given in Fig. 4 below:

Conclusion
As the theoretical relation between Liutex and several 2G vortex identification methods are given, and a numerical example of Burger vortex is implemented, the following conclusions can be summarized.

1) Mathematical relations between Liutex and the 2G vortex identification methods,
such as Q, Δ, λ 2 and λ ci methods are derived mathematically, which clearly show the second generation methods are severely contaminated by shearing, stretching, or by both. On the other hand, Liutex is the exact mathematical definition of fluid rotation or vortex. 2) The velocity gradient tensor in the Principal Coordinate can be uniquely decomposed into its rigid rotational part, stretching/compression part, and shearing part. This decomposition is known as the Liutex-based decomposition.
3) The Q method conflicts with Δ, λ 2 , and λ ci methods, meaning that at some points, Q may mistreat the vortex point as a non-rotational point. Evidence of this can be seen in the graphical representations (see Fig. 2). 4) According to Eqs. (43), (45), (48), (54), we can have the following table (Table 1) as Liutex magnitude is defined as twice the angular speed of local rigid rotation. 5) The external shear may change Liutex under some conditions, especially when large shear force is applied in the opposite direction of rotation, which means the external shear can change the rotation strength but does not mean Liutex is  contaminated by shear, which can be seen clearly in the formula of the Liutex magnitude. 6) According to Table 1, Q, Δ, and λ 2 are stretch sensitive, while Q, Δ, λ 2 , and λ ci are affected by the shearing factor. The numerical results from the Burgers vortex coincide with the theoretical analysis. However, Liutex remained the same and is not affected. Therefore, we can conclude Q, λ 2 , and λ ci methods may not appropriately represent vortex strength, while Liutex is the only one that correctly represents fluid rotation or vortex.