留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

Schur Forms and Normal-Nilpotent Decompositions

LI Zhen

LI Zhen. Schur Forms and Normal-Nilpotent Decompositions[J]. Applied Mathematics and Mechanics, 2024, 45(9): 1200-1211. doi: 10.21656/1000-0887.450129
Citation: LI Zhen. Schur Forms and Normal-Nilpotent Decompositions[J]. Applied Mathematics and Mechanics, 2024, 45(9): 1200-1211. doi: 10.21656/1000-0887.450129
李震. Schur形式与正规-幂零分解[J]. 应用数学和力学, 2024, 45(9): 1200-1211. doi: 10.21656/1000-0887.450129
引用本文: 李震. Schur形式与正规-幂零分解[J]. 应用数学和力学, 2024, 45(9): 1200-1211. doi: 10.21656/1000-0887.450129

Schur Forms and Normal-Nilpotent Decompositions

1 A matrix is said to be quasiorthogonal if its columns are mutually orthogonal and so are its rows. But the columns and rows are not required to be normalized to unit length. The accurate name for this class of matrices should be orthogonal, unfortunately, which has been widely accepted for matrices that should have been called orthonormal matrices.
(Recommended by WU Chuijie, M. AMM Editorial Board)

Schur形式与正规-幂零分解

doi: 10.21656/1000-0887.450129
详细信息
    作者简介:
  • 中图分类号: O302; O35

  • 摘要: 实的和复的Schur形式近年来受到流体力学界(特别是与旋涡和湍流相关)越来越多的关注。几个速度梯度张量分解(例如三元运动分解TDM和正规-幂零分解NND)被提出用于分析流体微元的局部运动。然而,由于Schur形式存在不同类型和非唯一性,以及NND有多种可能定义,一些混淆广泛传播并正在对研究造成危害。该工作旨在清除这种混淆。为此,复的和实的Schur形式由很基本的知识构造性地推导出来,其非唯一性被特别加以考虑,唯一性条件被提出。在对正规性和幂零性加以一般讨论后,一个复NND和几个实NND以及正规-非正规分解被构造出来,并简要地比较了复的和实的分解。在这些基础上,几个混淆点得到澄清,例如NND与TDM的差异以及复的和实的NND之间的内在鸿沟。此外,笔者提议将复本征值情况下实的块Schur形式及其对应的NND拓展到实本征值情况,不过其合理性有待进一步研究。
  • In the thesis of Li[1], a normal-nilpotent decomposition (NND) of real square matrices is proposed and applied to the velocity gradient tensor of a fluid element for analyzing vortex criteria (also see ref. [2]). The decomposition can be stated as that every 3-by-3 real matrix is a sum of a normal matrix and a nilpotent matrix. As a special application, the velocity gradient tensor $ {\boldsymbol{A}}=\nabla {\boldsymbol{v}} $ can be decomposed as

    $$ {\boldsymbol{A}} = {\boldsymbol{N}} + {\boldsymbol{S}}, $$ (1)

    where $ {\boldsymbol{N}} $ is a normal tensor and $ {\boldsymbol{S}} $ is a nilpotent tensor. Then $ {\boldsymbol{N}} $ can be further decomposed into three components and results in a quadruple decomposition of $ {\boldsymbol{A}} $ as described in ref. [1-2]. The decomposition shows itself as a powerful framework for analyzing local motions of fluid elements and vortex criteria.

    The NND proposed in ref. [1-2] is derived from a canonical form of real matrices discovered earlier by Murnaghan, Wintner[3], which is a counterpart over the field of real numbers to the famous Schur form for complex matrices[4]. In that paper, Isaai Schur generalized the previous results restricted to orthogonal[5] and Hermitian[6] matrices. The canonical form of Murnaghan, Wintner[3] is often referred to as the real Schur form. It is known that both the complex and real Schur forms of a given matrix are not unique. Actually, the NND is based on a special real Schur form among the possible variants. However, ref. [1-2] didn’t present a clarification of which variant of real Schur forms has been chosen and how to maintain uniqueness and consistency throughout the flow fields.

    Recently, real and complex Schur forms as well as NND have received increasing attention. For example, Keylock[7-8] employed the complex Schur form of velocity gradient tensor in their statistical model of fine structures in turbulence. Das, Girimaji[9] analyzed fine structures of turbulence in the framework of NND. Hoffman[10] used NND to analyze energy stability in incompressible turbulence. Zhu[11-12] derived a real Schur flow (RSF) as the compressible Taylor-Proudman limit in magnetohydrodynamics, for which the velocity gradient tensor takes a globally uniform real Schur form. Kronborg et al. [13] applied NND to analyse the shear in blood flows. Arun, Colonius[14] analyze the collision of vortex rings with NND. Just to mention a few.

    Apart from the direct application of Schur forms and NND as mentioned above, there are also studies of fluid motions that are related to or similar in some sense to NND. For example, Kolář[15-16] proposed a triple decomposition of motion (TDM) of the velocity gradient tensor, which is very similar to the quadruple NND proposed in ref. [1-2]. Unfortunately, several works, such as ref. [9-10, 13-14], misidentified the tools they used as TDM, when in fact they were NND.

    Kronborg, Hoffman[17] made an effort to clarify the relationship between TDM and the Schur forms. But they overlooked the gap between the two; especially, they misidentified NND as TDM. Also, they didn’t reveal the intrinsic gap between the NNDs derived from real and complex Schur forms properly, and treated it merely as an issue of algorithmic convenience. Besides, they failed to characterize the non-uniqueness of these forms completely and to present a satisfactory standardization procedure, although they claimed so.

    Considering the different types and non-uniqueness of Schur forms as well as the various possible definitions of NNDs, a clarification of their relationships is in demand, which is the aim of the current article. We hope that this article will help eliminate the confusion that has existed in previous research and clear the way for further investigations and applications of Schur forms and NNDs in fluid mechanics and other areas. The readers may also have an interest in the geometric or kinematical interpretations of Schur forms, which is out of the scope of this article. We refer to Zou, Xu, Tang[18], which provides a nice discussion on streamline patterns of real Schur forms.

    The structure of the article is as follows. Since NND is based on a real Schur form, we first describe the Schur forms. In section 1, the general complex Schur form for a 3-by-3 real matrix is derived with a discussion of its non-uniqueness. In section 2, the general real Schur form for a 3-by-3 real matrix is derived with a discussion of its non-uniqueness. A condition of uniqueness and the special real Schur form determined by this condition are proposed. In section 3, a brief discussion of normal and nilpotent matrices are presented firstly. Based on the general complex Schur form, a construction of complex NND is presented. Based on the special real Schur form proposed in section 2, several types of real NNDs are proposed. Besides, we also proposed two normal-nonnormal decompositions. Then, the distinction between the NNDs and the TDM proposed by Kolář[15-16] is demonstrated. In section 4, we demonstrate that it is generally impossible to convert a complex NND into a real NND through unitary transformations, which implies that the two types of NNDs are nonequivalent.

    For later use, let’s first restate a classical result about matrices. Let $ {\boldsymbol{A}}\in \mathbb{C}^{n\times n} $ be a complex $ n $-by-$ n $ matrix. A theorem of Schur[4] says that there is a unitary matrix $ {\boldsymbol{U}} $ such that $ \tilde{{\boldsymbol{A}}}= {\boldsymbol{U}} {\boldsymbol{A}} {\boldsymbol{U}}^\dagger $ is upper triangular, i.e., $ \tilde{A}_{ij}=0 $ for all $ i>j $, and the diagonal entries $ \tilde{A}_{ii} $’s are eigenvalues of $ {\boldsymbol{A}} $. Denote the $ i $-th row of $ {\boldsymbol{U}} $ as $ {\boldsymbol{u}}_i $, then

    $$ {\boldsymbol{A}} = {\boldsymbol{U}}^\dagger \tilde{{\boldsymbol{A}}} {\boldsymbol{U}} = \sum\limits_{i,j:i\le j} \tilde{A}_{ij} {\boldsymbol{u}}_i^\dagger {\boldsymbol{u}}_j . $$ (2)

    The eq.(2) is sometimes stated as a theorem and referred to as the representation theorem of matrices.

    Theorem 1 Every matrix $ {\boldsymbol{A}}\in \mathbb{C}^{n\times n} $ has a representation as eq.(2) in some orthonormal basis $ {\boldsymbol{u}}_i $.

    Indeed, it is an extension of the representation theorem of self-adjoint matrices, which says that every self-adjoint matrix has a diagonal representation in certain orthonormal basis. An immediate consequence of the theorem is

    Corollary 1 In the representation formula eq.(2), $ {\boldsymbol{u}}_n $ is a unit eigenvector of $ {\boldsymbol{A}} $.

    Since $ {\boldsymbol{u}}_i $ forms an orthonormal basis of $ \mathbb{C}^n $, it is easy to see that $ {\boldsymbol{u}}_n {\boldsymbol{A}}= \tilde{A}_{nn} {\boldsymbol{u}}_n $, i.e., $ {\boldsymbol{u}}_n $ must be a unit eigenvector of $ {\boldsymbol{A}} $. Without specification, we always refer to left eigenvectors.

    Let $ {\boldsymbol{A}}\in \mathbb{R}^{3\times 3} $ be a real 3-by-3 matrix. It can be the component matrix of a second-order tensor in the basis aligned with the Cartesian coordinates $ x,y,z $. Since the characteristic polynomial of $ {\boldsymbol{A}} $ is cubic, $ {\boldsymbol{A}} $ always has a real eigenvalue $ \lambda_3=\lambda_{ {\rm{r}}} $. Let $ {\boldsymbol{e}}_{ {\rm{r}}} $ be a unit eigenvector associated with $ \lambda_{ {\rm{r}}} $. If $ {\boldsymbol{e}}_{ {\rm{r}}} $ is required to be real, then it has two possible choices, which are in opposite directions. If it is allowed to be complex, then it can only be determined up to an arbitrary unitary factor $ {\rm{e}}^{ {\rm{i}}\phi_3} $. Notice that we can always require that $ {\boldsymbol{e}}_{ {\rm{r}}} $ is a real vector. But we will not make this choice until necessary.

    Fix $ {\boldsymbol{e}}_{ {\rm{r}}} $, let $ {\boldsymbol{u}}, {\boldsymbol{v}}\in \mathbb{C}^3 $ such that

    $${\boldsymbol{U}}:=\left(\begin{matrix}{\boldsymbol{u}} \\ {\boldsymbol{v}} \\ {\boldsymbol{e}}_{{\rm{r}}}\end{matrix}\right)$$ (3)

    is a unitary matrix. There is

    $$ \tilde{{\boldsymbol{A}}} := {\boldsymbol{U}} {\boldsymbol{A}} {\boldsymbol{U}}^\dagger = \left(\begin{matrix} {\tilde{{\boldsymbol{A}}}_2 }& { {\boldsymbol{c}}^\dagger }\\ { {\bf{0}}} & {\lambda_{ {\rm{r}}} } \end{matrix}\right) $$ (4)

    with submatrices (i.e., blocks)

    $$\tilde{{\boldsymbol{A}}}_2 := \left(\begin{matrix}{{\boldsymbol{u}} {\boldsymbol{A}} {\boldsymbol{u}}^\dagger} & {{\boldsymbol{u}} {\boldsymbol{A}} {\boldsymbol{v}}^\dagger }\\ { {\boldsymbol{v}} {\boldsymbol{A}} {\boldsymbol{u}}^\dagger} & {{\boldsymbol{v}} {\boldsymbol{A}} {\boldsymbol{v}}^\dagger } \end{matrix}\right),\; {\boldsymbol{c}}^\dagger := \left(\begin{matrix} { {\boldsymbol{u}} {\boldsymbol{A}} {\boldsymbol{e}}_{ {\rm{r}}}^\dagger }\\ { {\boldsymbol{v}} {\boldsymbol{A}} {\boldsymbol{e}}_{ {\rm{r}}}^\dagger } \end{matrix}\right), $$ (5)

    where $ ^\dagger $ stands for adjoint (i.e., conjugate transpose). The choices of $ {\boldsymbol{u}}, {\boldsymbol{v}} $ are not unique. If $ {\boldsymbol{V}}_2\in \mathbb{C}^{2\times2} $ is a unitary matrix, then $ {\boldsymbol{V}} {\boldsymbol{U}} $ with

    $$ {\boldsymbol{V}}:= \left(\begin{matrix} { {\boldsymbol{V}}_2 }& { {\bf{0}}} \\ { {\bf{0}} }&{ 1 } \end{matrix}\right) $$ (6)

    will do the equivalent job as $ {\boldsymbol{U}} $ does, i.e., eliminate the first and second elements of the third row of $ {\boldsymbol{A}} $. Conversely, any unitary transformation on $ \mathbb{C}^3 $ that leaves $ {\boldsymbol{e}}_{ {\rm{r}}} $ unchanged must have the form of $ {\boldsymbol{V}} $.

    Besides of $ \lambda_3=\lambda_{ {\rm{r}}} $, $ {\boldsymbol{A}} $ has another two eigenvalues $ \lambda_1,\lambda_2\in \mathbb{C} $, which could be real or a pair of complex conjugate numbers. They are also eigenvalues of $ \tilde{{\boldsymbol{A}}}_2 $. If $ \lambda_1,\lambda_2\notin \mathbb{R} $, then $ {\boldsymbol{A}} $ can not be transformed into an upper triangular matrix by real orthogonal transformations. With unitary transformations, however, $ {\boldsymbol{A}} $ can be converted to an upper triangular form.

    Let

    $$ {\boldsymbol{V}}_2 = \left(\begin{matrix} {\boldsymbol{a}} \\ {\boldsymbol{b}} \end{matrix}\right) \; ,$$ (7)

    with $ {\boldsymbol{a}}, {\boldsymbol{b}}\in \mathbb{C}^2 $ and

    $$\tilde{\tilde{{\boldsymbol{A}}}}_2 := {\boldsymbol{V}}_2 \tilde{{\boldsymbol{A}}}_2 {\boldsymbol{V}}_2^\dagger = \left(\begin{matrix} {{\boldsymbol{a}} \tilde{{\boldsymbol{A}}}_2 {\boldsymbol{a}}^\dagger }& {{\boldsymbol{a}} \tilde{{\boldsymbol{A}}}_2 {\boldsymbol{b}}^\dagger} \\{{\boldsymbol{b}} \tilde{{\boldsymbol{A}}}_2 {\boldsymbol{a}}^\dagger }& { {\boldsymbol{b}} \tilde{{\boldsymbol{A}}}_2 {\boldsymbol{b}}^\dagger } \end{matrix}\right). $$ (8)

    If $ \tilde{\tilde{{\boldsymbol{A}}}}_2 $ is upper triangular, then $ {\boldsymbol{b}} $ must be a unit eigenvector of $ \tilde{{\boldsymbol{A}}}_2 $. Notice that $ \tilde{{\boldsymbol{A}}}_2 $ is just the restriction of $ {\boldsymbol{A}} $ in the subspace orthogonal to $ {\boldsymbol{e}}_{ {\rm{r}}} $, which is independent of the choice of $ {\boldsymbol{u}}, {\boldsymbol{v}} $. Therefore, $ {\boldsymbol{b}} $ only depends on $ {\boldsymbol{A}} $. Again, $ {\boldsymbol{b}} $ can only be determined up to an arbitrary unitary factor $ {\rm{e}}^{ {\rm{i}} \phi_2} $. $ {\boldsymbol{a}} $ is a unit vector orthogonal to $ {\boldsymbol{b}} $, i.e., $ {\boldsymbol{a}} {\boldsymbol{a}}^\dagger = 1 $ and $ {\boldsymbol{a}} {\boldsymbol{b}}^\dagger = 0 $, hence $ {\boldsymbol{a}} $ is uniquely determined up to an arbitrary unitary factor $ {\rm{e}}^{ {\rm{i}} \phi_1} $.

    Being an upper triangular matrix, $ {\boldsymbol{a}} \tilde{{\boldsymbol{A}}}_2 {\boldsymbol{a}}^\dagger $ and $ {\boldsymbol{b}} \tilde{{\boldsymbol{A}}}_2 {\boldsymbol{b}}^\dagger $ must be eigenvalues of $ \tilde{\tilde{{\boldsymbol{A}}}}_2 $ and $ \tilde{{\boldsymbol{A}}}_2 $. But their order is undetermined. Let’s denote $ {\boldsymbol{a}} \tilde{{\boldsymbol{A}}}_2 {\boldsymbol{a}}^\dagger = \lambda_1 $ and $ {\boldsymbol{b}} \tilde{{\boldsymbol{A}}}_2 {\boldsymbol{b}}^\dagger = \lambda_2 $. There is

    $$\tilde{\tilde{{\boldsymbol{A}}}} := {\boldsymbol{V}} \tilde{{\boldsymbol{A}}} {\boldsymbol{V}}^\dagger = \left(\begin{matrix} { \tilde{\tilde{{\boldsymbol{A}}}}_2 }& { {\boldsymbol{V}}_2 {\boldsymbol{c}}^\dagger }\\ { {\bf{0}} }& {\lambda_{ {\rm{r}}} } \end{matrix}\right) = \left(\begin{matrix} { \lambda_1 }& {{\boldsymbol{a}} \tilde{{\boldsymbol{A}}}_2 {\boldsymbol{b}}^\dagger} & { {\boldsymbol{a}} {\boldsymbol{c}}^\dagger }\\ { 0} & {\lambda_2} & { {\boldsymbol{b}} {\boldsymbol{c}}^\dagger} \\ { 0 }&{ 0} & {\lambda_{\rm{r}} } \end{matrix}\right), $$ (9)

    which is a complex Schur form of $ {\boldsymbol{A}} $. Notice that the complex Schur form is not unique. Apart from permutations of diagonal entries, there are also undetermined unitary factors. For all $ \phi_1,\phi_2,\phi_3\in[0,2{\text{π}}) $, let

    $$ {\boldsymbol{P}} := {\rm{diag}}\,({{\rm{e}}^{ {\rm{i}} \phi_1},{\rm{e}}^{ {\rm{i}} \phi_2},{\rm{e}}^{ {\rm{i}} \phi_3}}).$$ (10)

    Then $ {\boldsymbol{P}} {\boldsymbol{V}} {\boldsymbol{U}} $ will do the equivalent job as $ {\boldsymbol{V}} {\boldsymbol{U}} $ to transform $ {\boldsymbol{A}} $ into its complex Schur form. Hence the general complex Schur form of $ {\boldsymbol{A}} $ can be written as

    $$ \tilde{\tilde{{\boldsymbol{A}}}} = \left(\begin{matrix}{ \lambda_1 } &{ a_{12} {\rm{e}}^{ {\rm{i}} (\phi_1-\phi_2)} } &{ a_{13} {\rm{e}}^{ {\rm{i}} (\phi_1-\phi_3)} }\\ { 0 } &{ \lambda_2 } &{ a_{23} {\rm{e}}^{ {\rm{i}} (\phi_2-\phi_3)} }\\ { 0 } &{ 0 } &{ \lambda_{ {\rm{r}}} } \end{matrix}\right). $$ (11)

    Notice that only the differences between $ \phi_1,\phi_2,\phi_3 $ enter the results. Let $ \phi_{13}=\phi_1-\phi_3,\,\phi_{23}=\phi_2-\phi_3 $, we obtain

    $$ \tilde{\tilde{{\boldsymbol{A}}}} = \left(\begin{matrix} { \lambda_1 } &{ a_{12} {\rm{e}}^{ {\rm{i}} (\phi_{13}-\phi_{23})} } &{ a_{13} {\rm{e}}^{ {\rm{i}} \phi_{13}} }\\ { 0 } &{ \lambda_2 } &{ a_{23} {\rm{e}}^{ {\rm{i}} \phi_{23}} }\\ { 0 } &{ 0 } &{ \lambda_3} \end{matrix}\right). $$ (12)

    Since each of $ {\boldsymbol{a}}, {\boldsymbol{b}}, {\boldsymbol{c}} $ can be determined up to an arbitrary unitary factor, we can require that any two of $ a_{12},a_{13},a_{23} $ are real. Let’s assume that $ a_{13},a_{23}\in \mathbb{R} $. But in general, $ a_{12} $ is complex. Assume that $ a_{12}=\mu+ {\rm{i}}\nu $ with $ \mu,\nu\in \mathbb{R} $.

    If $ \lambda_1,\lambda_2\in \mathbb{R} $, then the $ {\boldsymbol{U}}, {\boldsymbol{V}} $ can be chosen as special orthogonal matrices, hence $ a_{12}\in \mathbb{R} $. Furthermore, we can set $ \phi_{13}=m{\text{π}},\phi_{23}=n{\text{π}} $, then

    $$ \tilde{\tilde{{\boldsymbol{A}}}} = \left(\begin{matrix} { \lambda_1 } &{ a_{12} (-1)^{m-n} } &{ a_{13} (-1)^m }\\ { 0 } &{ \lambda_2 } &{ a_{23} (-1)^n }\\ { 0 } &{ 0 } &{ \lambda_3 } \end{matrix}\right). $$ (13)

    This is the real Schur form of $ {\boldsymbol{A}} $ when $ \lambda_1,\lambda_2\in \mathbb{R} $, which is a special case of eq.(12).

    As mentioned at the end of section 1, When $ {\boldsymbol{A}}\in \mathbb{R}^{3\times3} $ has three real eigenvalues, its general real Schur form is a particular case of its general complex Schur form, given in eq.(13). When $ {\boldsymbol{A}} $ has a pair of complex conjugate eigenvalues, it can not be transformed into an upper triangular matrix in the field of real numbers. In this case, we can only have a block upper triangular matrix in the form of $ \tilde{{\boldsymbol{A}}} $ obtained in section 1. Especially, we can require that $ \tilde{{\boldsymbol{A}}} $ is real by restricting $ {\boldsymbol{U}} $ to orthogonal matrices. It gives as the general real Schur form of $ {\boldsymbol{A}} $ with a pair of complex conjugate eigenvalues. But this leaves us one degree of freedom, which is the rotation in the subspace normal to $ {\boldsymbol{e}}_{ {\rm{r}}} $. Therefore, an orthogonal transform can be applied to $ \tilde{{\boldsymbol{A}}} $ to meet one more requirement. Let

    $$ {\boldsymbol{Q}} = \left(\begin{matrix} { \cos\,\varrho } &{ \sin\,\varrho } &{ 0 }\\ { -\sin\,\varrho } &{ \cos\,\varrho } &{ 0 }\\ { 0 } &{ 0 } &{ 1 } \end{matrix}\right), $$ (14)

    then

    $$ \hat{{\boldsymbol{A}}} := {\boldsymbol{Q}} \tilde{{\boldsymbol{A}}} {\boldsymbol{Q}}^ {\rm{T}}= $$ (15)
    $$\qquad \left(\begin{matrix} {\begin{array}{*{20}{c}}{\tilde{A}_{11} \cos ^2 \varrho+\tilde{A}_{22} \sin\, \varrho+}\\{(\tilde{A}_{12}+\tilde{A}_{21}) \cos \,\varrho \sin \,\varrho}\end{array}} & {\begin{array}{*{20}{c}}{\tilde{A}_{12} \cos ^2 \varrho-\tilde{A}_{21} \sin ^2 \varrho+}\\{(\tilde{A}_{22}-\tilde{A}_{11}) \cos \,\varrho \sin \,\varrho}\end{array}} & {c_1 \cos \,\varrho+c_2 \sin \,\varrho} \\ {\begin{array}{*{20}{c}}{ \tilde{A}_{21} \cos ^2 \varrho-\tilde{A}_{12} \sin ^2 \varrho+}\\{ (\tilde{A}_{22}-\tilde{A}_{11}) \cos \,\varrho \sin \,\varrho} \end{array}}& {\begin{array}{*{20}{c}}{\tilde{A}_{22} \cos ^2 \varrho+\tilde{A}_{11} \sin ^2 \varrho-}\\{(\tilde{A}_{12}+\tilde{A}_{21}) \cos \,\varrho \sin \,\varrho} \end{array}}& {{c_2 \cos \,\varrho-c_1 \sin \,\varrho }}\\ {0} & {0} & {\lambda_{{\rm{r}}}}\end{matrix}\right).$$ (16)

    We require that $ \hat{A}_{11}= \hat{A}_{22} $, hence

    $$ ( \tilde{A}_{11}- \tilde{A}_{22})({\cos^2\varrho - \sin^2\varrho}) + 2( \tilde{A}_{12}+ \tilde{A}_{21})\cos\,\varrho\sin\,\varrho = 0, $$ (17)

    which implies that

    $$ \tan\,(2\varrho) = \frac{ \tilde{A}_{22}- \tilde{A}_{11}}{ \tilde{A}_{12}+ \tilde{A}_{21}}. $$ (18)

    Since the period of $ \tan $ is $ {\text{π}} $, there are four values of $ \varrho\in[0,2{\text{π}}) $ that satisfy the equation. They determine two lines that are orthogonal to each other. We can also obtain

    $$ \cos^2\varrho = \frac{1}{2}\left({1\pm\frac{ \tilde{A}_{12}+ \tilde{A}_{21}}{| \tilde{{\boldsymbol{A}}}_2|}}\right),\;\; \sin\,\varrho\cos\,\varrho = \pm\frac{ \tilde{A}_{22}- \tilde{A}_{11}}{2| \tilde{{\boldsymbol{A}}}_2|}, $$ (19)

    where

    $$ | \tilde{{\boldsymbol{A}}}_2|:=\sqrt{\|{ \tilde{{\boldsymbol{A}}}_2}\|_F^2-2\det \tilde{{\boldsymbol{A}}}_2}=\sqrt{({ \tilde{A}_{12}+ \tilde{A}_{21}})^2 + ({ \tilde{A}_{22}- \tilde{A}_{11}})^2} $$ (20)

    is invariant under unitary transformations, hence the components of $ \tilde{{\boldsymbol{A}}} $ in its expression can be replaced by those of $ \hat{{\boldsymbol{A}}} $. With these values of $ \varrho $, we have

    $$ \hat{A}_{11} = \hat{A}_{22} = \frac{ \tilde{A}_{11}+ \tilde{A}_{22}}{2},\;\; \hat{A}_{12} = \frac{ \tilde{A}_{12}- \tilde{A}_{21}\pm| \tilde{{\boldsymbol{A}}}_2|}{2}, $$ (21)
    $$ \hat{A}_{21} = \frac{ \tilde{A}_{21}- \tilde{A}_{12}\pm| \tilde{{\boldsymbol{A}}}_2|}{2}. $$ (22)

    Therefore,

    $$ \hat{{\boldsymbol{A}}} = \left(\begin{matrix} { \dfrac{ \tilde{A}_{11}+ \tilde{A}_{22}}{2} } &{ \dfrac{ \tilde{A}_{12}- \tilde{A}_{21}\pm| \tilde{{\boldsymbol{A}}}_2|}{2} } &{ \hat{A}_{13} }\\ { \dfrac{ \tilde{A}_{21}- \tilde{A}_{12}\pm| \tilde{{\boldsymbol{A}}}_2|}{2} } &{ \dfrac{ \tilde{A}_{11}+ \tilde{A}_{22}}{2} } &{ \hat{A}_{23}} \\ { 0 } &{ 0 } &{ \lambda_{ {\rm{r}}} } \end{matrix}\right). $$ (23)

    Denote

    $$ \chi := \hat{A}_{11} = \hat{A}_{22},\;\; \omega_3:= \hat{A}_{12}- \hat{A}_{21},\;\; \gamma := \hat{A}_{12} + \hat{A}_{21} = \pm| \hat{{\boldsymbol{A}}}_2|. $$ (24)

    The 2-by-2 diagonal block of $ \hat{{\boldsymbol{A}}} $ can be written as

    $$ \hat{{\boldsymbol{A}}}_2 := \left(\begin{matrix} { \chi} &{ \dfrac{\gamma + \omega_3}{2}} \\ { \dfrac{\gamma - \omega_3}{2}} & {\chi } \end{matrix}\right) . $$ (25)

    We have mentioned that $ {\boldsymbol{e}}_{ {\rm{r}}} $ has two possible choices, which are in opposite directions. If $ \omega_3<0 $, rotate the basis about the $ \hat{x} $-axis for an angle of $ {\text{π}} $. The operation change the signs of $ \omega_3 $ and $ \gamma $ simultaneously. Therefore, the choice of $ {\boldsymbol{e}}_{ {\rm{r}}} $ is uniquely determined by requiring $ \omega_3>0 $. We have also mentioned that $ \varrho $ has four possible choices, which determine two lines that are orthogonal to each other. If $ \gamma<0 $, rotate the basis about the $ \hat{z} $-axis for an angle of $ {\text{π}}/2 $. This operation changes the sign of $ \gamma $ while keeping the sign of $ \omega_3 $. Therefore, the choice of $ \varrho $ is reduced to two opposite directions by requiring $ \gamma>0 $. Denote

    $$ \hat{A}_{13} := -\beta, \;\; \hat{A}_{23} := \alpha. $$ (26)

    The two choices of $ \varrho $ lead to different signs of $ \alpha $ and $ \beta $ simultaneously. Therefore, the choice of $ \varrho $ can be uniquely determined by assigning a sign for $ \alpha $ or $ \beta $. We require that $ \alpha>0 $. According to the discussions above, in the generic case that all of $ \omega_3,\gamma,\alpha $ are nonzero,

    $$ \hat{{\boldsymbol{A}}} = \left(\begin{matrix} { \chi } &{ \dfrac{\gamma + \omega_3}{2} } &{ -\beta }\\ { \dfrac{\gamma - \omega_3}{2} } &{ \chi } &{ \alpha }\\ { 0 } &{ 0 } &{ \lambda_{ {\rm{r}}} } \end{matrix}\right) $$ (27)

    can be uniquely determined by requiring

    $$ \omega_3>0,\; \gamma>0,\; \alpha>0. $$ (28)

    Of course, this uniqueness condition is not the only choice. For example, one can require $ \omega_3<0 $ instead. But that only affects the appearance of the resulting standard matrix form and is merely a matter of convention. If some of $ \omega_3,\gamma,\alpha $ is zero, then the form of $ \hat{{\boldsymbol{A}}} $ can be nonunique. We should mention that Kronborg, Hoffman[17] also noticed the nonuniqueness of Schur forms, but their consideration is restricted to a permutation of basis vectors, which is not essential to the forms and the construction of NNDs.

    When $ {\boldsymbol{A}} $ has a pair of complex conjugate eigenvalues $ \lambda_{1,2} = \lambda_{ {\rm{c}} {\rm{r}}}\pm {\rm{i}}\lambda_{ {\rm{c}} {\rm{i}}} $, eq.(27) is a special form of the real Schur forms of $ {\boldsymbol{A}} $. In this case, there is

    $$ \lambda_{ {\rm{c}} {\rm{r}}} = \chi,\; \lambda_{ {\rm{c}} {\rm{i}}} =\frac{1}{2} \sqrt{\omega_3^2-\gamma^2}. $$ (29)

    We should emphasize that when $ {\boldsymbol{A}} $ has three real eigenvalues, the form of $ \hat{{\boldsymbol{A}}} $ is still uniquely determined by the above procedure. In other words, $ \hat{{\boldsymbol{A}}} $ is a unique canonical form of $ {\boldsymbol{A}} $ even when $ {\boldsymbol{A}} $ has three real eigenvalues. But in this case, it is not a real Schur form of $ {\boldsymbol{A}} $. In this case,

    $$ \lambda_{1,2} = \chi \pm \frac{1}{2} \sqrt{\gamma^2-\omega_3^2}. $$ (30)

    Remark 1 In eq.(27), notice that

    $$ 4 \hat{A}_{12} \hat{A}_{21} = \gamma^2 - \omega_3^2 $$ (31)

    is the discriminant of the characteristic polynomial of $ \hat{{\boldsymbol{A}}}_2 $. It can be used to determine whether $ {\boldsymbol{A}} $ has a pair of complex conjugate eigenvalues.

    Remark 2 Since $ \hat{{\boldsymbol{A}}}=( {\boldsymbol{Q}} {\boldsymbol{U}}) {\boldsymbol{A}}( {\boldsymbol{Q}} {\boldsymbol{U}})^{\rm{T}} $, $ {\boldsymbol{Q}} $ is a rotation matrix and $ {\boldsymbol{U}} $ can be chosen as a rotation matrix, therefore, the six nonzero components $ \chi,\omega_3,\lambda_{ {\rm{r}}},\alpha,\beta,\gamma $ of $ \hat{{\boldsymbol{A}}} $ are rotational invariants of $ {\boldsymbol{A}} $.

    A matrix $ {\boldsymbol{A}}\in \mathbb{C}^{n\times n} $ is said to be normal if $ {\boldsymbol{A}} {\boldsymbol{A}}^\dagger= {\boldsymbol{A}}^\dagger {\boldsymbol{A}} $. $ {\boldsymbol{A}} $ is normal if its complex Schur form is diagonal or its real Schur form is block diagonal. $ {\boldsymbol{A}} $ is said to be nilpotent if $ {\boldsymbol{A}}^n=0 $. $ {\boldsymbol{A}} $ is nilpotent if it only has zero eigenvalues, if its Schur form is strictly upper triangular. Thus we can conclude that if $ {\boldsymbol{A}} $ is both normal and nilpotent, then $ {\boldsymbol{A}}={\boldsymbol{0}} $. Therefore, it is very natural to decompose $ {\boldsymbol{A}} $ as a sum of a normal matrix and a nilpotent matrix by splitting the diagonal and off-diagonal parts of its complex Schur forms or splitting block diagonal and off-diagonal parts of its real Schur forms. However, the splitting approach is not the only way to realize such a decomposition. Actually, there may be more than one such decomposition for each Schur form of $ {\boldsymbol{A}} $. Furthermore, considering the non-uniqueness of Schur forms, based on each variant of complex or real Schur forms of $ {\boldsymbol{A}} $, such a decomposition can be constructed. All of them can be called normal-nilpotent decompositions (NNDs).

    There is another set of necessary and sufficient conditions for normal matrices, which may serve as guidance when constructing NNDs. Let

    $$ {\boldsymbol{S}} = \frac{ {\boldsymbol{A}}+ {\boldsymbol{A}}^\dagger}{2},\; {\boldsymbol{W}} = \frac{ {\boldsymbol{A}}- {\boldsymbol{A}}^\dagger}{2}. $$ (32)

    Then $ {\boldsymbol{A}} = {\boldsymbol{S}} + {\boldsymbol{W}} $ and $ {\boldsymbol{S}}^\dagger= {\boldsymbol{S}}, {\boldsymbol{W}}^\dagger=- {\boldsymbol{W}} $. It is evident that $ {\boldsymbol{A}} $ is normal if $ {\boldsymbol{S}} {\boldsymbol{W}}= {\boldsymbol{W}} {\boldsymbol{S}} $. In the orthonormal basis composed of eigenvectors of $ {\boldsymbol{S}} $, suppose

    $$ {\boldsymbol{S}} = \left(\begin{matrix} {s_1 } &{ 0 } &{ 0 }\\{0 } &{ s_2 } &{ 0 }\\{0 } &{ 0 } &{ s_3 } \end{matrix}\right),\; {\boldsymbol{W}} = \left(\begin{matrix} { 0 } &{ w_3 } &{ -w_2 }\\{-w_3 } &{ 0 } &{ w_1 }\\ { w_2 } &{ -w_1 } &{ 0 } \end{matrix}\right). $$ (33)

    Then

    $$ {\boldsymbol{S}} {\boldsymbol{W}}- {\boldsymbol{W}} {\boldsymbol{S}} = \left(\begin{matrix} {0} &{ w_3(s_1-s_2)} & w_2(s_3-s_1) \\{w_3(s_1-s_2)} & {0} & {w_1(s_2-s_3) }\\{w_2(s_3-s_1) }&{ w_1(s_2-s_3)} &{ 0} \end{matrix}\right), $$ (34)

    which leads to the following

    Theorem 2 The matrix $ {\boldsymbol{A}}\in \mathbb{C}^{3\times3} $ is normal only in one of the following three cases:

    1) All $ s_1,s_2,s_3 $ are distinct, $ w_1=w_2=w_3=0 $. In this case, $ {\boldsymbol{A}} $ is symmetric.

    2) $ s_i=s_j\ne s_k $ for distinct $ i,j,k $, $ w_i=w_j=0 $. In this case, $ {\boldsymbol{A}} $ is quasiorthogonal. 1

    3) $ s_1=s_2=s_3 $. That is $ {\boldsymbol{S}}=( {\rm{tr}}\, {\boldsymbol{A}}/3) {\boldsymbol{I}} $.

    Notice that the validity of these conditions does not depend on the diagonal form of $ {\boldsymbol{S}} $ since normality and nilpotency are independent of unitary transformations. But the $ s_i $’s and $ w_i $’s refer to the matrix entries in the basis of eigenvectors of $ {\boldsymbol{S}} $. This also justifies the procedure of constructing NND in the basis of Schur forms since they are generalizations of the basis of eigenvectors of self-adjoint matrices, as stated by the representation theorem eq.(2).

    Let $ \tilde{\tilde{{\boldsymbol{\varLambda}}}} $ be the diagonal part of $ \tilde{\tilde{{\boldsymbol{A}}}} $ given in eq.(12) and $ \tilde{\tilde{{\boldsymbol{\varDelta}}}} $ be the off-diagonal part. That is

    $$ \tilde{\tilde{{\boldsymbol{\varLambda}}}} = \left(\begin{matrix} {\lambda_1 } &{ 0 } &{ 0 }\\ { 0 } &{ \lambda_2 } &{ 0 }\\{0 } &{ 0 } &{ \lambda_3 } \end{matrix}\right),\; \tilde{\tilde{{\boldsymbol{\varDelta}}}} = \left(\begin{matrix}{0 } &{ a_{12} {\rm{e}}^{ {\rm{i}} (\phi_{13}-\phi_{23})} } &{ a_{13} {\rm{e}}^{ {\rm{i}} \phi_{13}} }\\{0 } &{ 0 } &{ a_{23} {\rm{e}}^{ {\rm{i}} \phi_{23}}} \\{0 } &{ 0 } &{ 0 }\end{matrix}\right). $$ (35)

    Then we have $ \tilde{\tilde{{\boldsymbol{A}}}} = \tilde{\tilde{{\boldsymbol{\varLambda}}}} + \tilde{\tilde{{\boldsymbol{\varDelta}}}} $. Obviously, $ \,\tilde{\tilde{{\boldsymbol{\varLambda}}}} $ is normal and $ \tilde{\tilde{{\boldsymbol{\varDelta}}}} $ is nilpotent. In the generic case, all $ \lambda_1,\lambda_2,\lambda_3 $ are distinct from each other. In the same sprite of the case 1 of theorem 2, the normal part should contains no off-diagonal entries. In particular, we can not add an anti-symmetric matrix to $ \,\tilde{\tilde{{\boldsymbol{\varLambda}}}} $. Besides, we can not add a symmetric matrix to $ \,\tilde{\tilde{{\boldsymbol{\varLambda}}}} $ either, otherwise it will violate the nilpotency of $ \tilde{\tilde{{\boldsymbol{\varDelta}}}} $. Hence the above construction of NND is the only natural one in this case. Of course, the decomposition is complex, and $ \tilde{\tilde{{\boldsymbol{\varDelta}}}} $ is not unique. When $ {\boldsymbol{A}} $ has three real eigenvalues, $ \,\tilde{\tilde{{\boldsymbol{\varLambda}}}} $ is real and we can make $ \tilde{\tilde{{\boldsymbol{\varDelta}}}} $ real too, as shown by eq.(13). When $ {\boldsymbol{A}} $ has a pair of complex conjugate eigenvalues, $\, \tilde{\tilde{{\boldsymbol{\varLambda}}}} $ can be converted to a real matrix through a unitary transformation, but $ \tilde{\tilde{{\boldsymbol{\varDelta}}}} $ can not be converted to a real matrix by the same unitary transformation. The details are presented in section 4.

    The NND used by Keylock[7-8] is of this type. It has several advantages, such as

    1) The Schur form has a simple upper triangular form with eigenvalues on its diagonal.

    2) There is only one natural way of decomposing the Schur form into the sum of a normal part and a nilpotent part.

    3) The normal part and the nilpotent part have no overlap, resulting in a clean decomposition.

    4) The cases of real and complex conjugate eigenvalues are treated in a unified way.

    As a mathematical model of physical quantities such as the velocity gradient tensor, the complex NND is adequate. However, the physical interpretation of the complex matrices is not clear yet. In contrast, real NNDs can have a much clearer physical meaning, as demonstrated in the literature.

    In Li[1], the real NND of $ {\boldsymbol{A}} $ with three eigenvalues is a special case of the complex NND as constructed above. When $ {\boldsymbol{A}} $ has a pair of complex conjugate eigenvalues, the real NND is constructed based on the special real Schur form eq.(27). However, in both the real and complex eigenvalue cases, those are not the only ways to construct real NNDs.

    Let’s consider the case of complex eigenvalues first. Since the special real Schur form eq.(27) has two identical diagonal entries, the normal part of $ \hat{{\boldsymbol{A}}} $ corresponds to the case 2 of theorem 2. In particular, the normal part of $ \hat{{\boldsymbol{A}}} $ can only have off-diagonal entries in the positions of 12 and 21. We have three natural ways to decompose the off-diagonal part of the submatrix $ \hat{{\boldsymbol{A}}}_2 $, which are:

    1) Quasiorthogonal-nilpotent;

    2) Symmetric-nilpotent;

    3) Symmetric-antisymmetric.

    Each type of decomposition above leads to a type of decomposition of $ \hat{{\boldsymbol{A}}} $. Since both symmetric and antisymmetric matrices are normal, the first two ways of decomposition belong to NND and the third does not.

    Denote

    $$ \phi := \frac{\omega_3 + \gamma}{2},\; \; \psi := \frac{\omega_3 - \gamma}{2}. $$ (36)

    The decomposition

    $$ \hat{{\boldsymbol{A}}} = \hat{{\boldsymbol{N}}} + \hat{{\boldsymbol{\varGamma}}} $$ (37)

    can be realized in the following ways.

    3.3.1   Quasiorthogonal-Nilpotent Decomposition

    The normal component $ \hat{{\boldsymbol{N}}} $ contains the diagonal of $ \hat{{\boldsymbol{A}}} $ and the anti-symmetric part of $ \hat{{\boldsymbol{A}}}_2 $, hence is orthogonal. There are two natural ways to separate an anti-symmetric part from $ \hat{{\boldsymbol{A}}}_2 $, which leads to two variants of NND: the one minimizes $ \| \hat{{\boldsymbol{N}}}\|_F $, the other maximizes $ \| \hat{{\boldsymbol{N}}}\|_F $. The $ \| \hat{{\boldsymbol{\varGamma}}}\|_F $ in the two approaches are the same. In both cases, $ \hat{{\boldsymbol{N}}} $ is normal and $ \hat{{\boldsymbol{\varGamma}}} $ is nilpotent.

    (a) Minimize the normal part:

    $$ \hat{{\boldsymbol{N}}} := \left(\begin{matrix}{ \chi } &{ \psi } &{ 0} \\{-\psi } &{ \chi } &{ 0 }\\{0 } &{ 0 } &{ \lambda_{ {\rm{r}}}}\end{matrix}\right),\; \; \hat{{\boldsymbol{\varGamma}}} := \left(\begin{matrix} {0 } &{ \gamma } &{ -\beta }\\{0 } &{ 0 } &{ \alpha }\\{0 } &{ 0 } &{ 0}\end{matrix}\right). $$ (38)

    In this case, $ \| \hat{{\boldsymbol{N}}}\|_F $ is minimized because $ |\gamma-\omega_3|<\gamma+\omega_3 $. This variant of NND is the one proposed in ref. [1].

    (b) Maximize the normal part:

    $$ \hat{{\boldsymbol{N}}} := \left(\begin{matrix} { \chi } &{ \phi } &{ 0 }\\ { -\phi } &{ \chi } &{ 0 }\\ { 0 } &{ 0 } &{ \lambda_{ {\rm{r}}} } \end{matrix}\right),\; \; \hat{{\boldsymbol{\varGamma}}} := \left(\begin{matrix} { 0 } &{ 0 } &{ -\beta} \\ { \gamma } &{ 0 } &{ \alpha }\\ { 0 } &{ 0 } &{ 0 } \end{matrix}\right). $$ (39)

    In this case, $ \| \hat{{\boldsymbol{N}}}\|_F $ is maximized because the off-diagonal entry of $ \hat{{\boldsymbol{A}}}_2 $ with a larger absolute value is kept in $ \hat{{\boldsymbol{N}}} $.

    3.3.2   Symmetric-Nilpotent Decomposition

    The normal component $ \hat{{\boldsymbol{N}}} $ contains the diagonal of $ \hat{{\boldsymbol{A}}} $ and the symmetric part of $ \hat{{\boldsymbol{A}}}_2 $. Similarly, there are two natural ways to separate a symmetric part from $ \hat{{\boldsymbol{A}}}_2 $, which leads to two variants of NND:

    (a) Minimize the normal part:

    $$ \hat{{\boldsymbol{N}}} := \left(\begin{matrix} { \chi } &{ -\psi } &{ 0} \\ { -\psi } &{ \chi } &{ 0 }\\ { 0 } &{ 0 } &{ \lambda_{ {\rm{r}}} } \end{matrix}\right),\; \hat{{\boldsymbol{\varGamma}}} := \left(\begin{matrix}{0 } &{ \omega_3 } &{ -\beta }\\ { 0 } &{ 0 } &{ \alpha }\\ { 0 } &{ 0 } &{ 0} \end{matrix}\right). $$ (40)

    (b) Maximize the normal part:

    $$ \hat{{\boldsymbol{N}}} := \left(\begin{matrix}{\chi } &{ \phi } &{ 0} \\{\phi } &{ \chi } &{ 0} \\{0 } &{ 0 } &{ \lambda_{ {\rm{r}}} } \end{matrix}\right),\; \hat{{\boldsymbol{\varGamma}}} := \left(\begin{matrix}{0 } &{ 0 } &{ -\beta }\\{-\omega_3 } &{ 0 } &{ \alpha }\\{0 } &{ 0 } &{ 0 } \end{matrix}\right). $$ (41)
    3.3.3   Normal-Nonnormal Decompositions

    The $ \hat{{\boldsymbol{N}}} $ contains the diagonal of $ \hat{{\boldsymbol{A}}} $ and the symmetric part or antisymmetric part of $ \hat{{\boldsymbol{A}}}_2 $. In both cases, $ \hat{{\boldsymbol{N}}} $ is normal but $ \hat{{\boldsymbol{\varGamma}}} $ is generally not nilpotent. Hence the resulting decompositions of $ \hat{{\boldsymbol{A}}} $ are not NNDs.

    (a) Symmetric-nonnormal:

    $$ \hat{{\boldsymbol{N}}} := \left(\begin{matrix} {\chi } &{ \dfrac{\gamma}{2} } &{ 0} \\{\dfrac{\gamma}{2} } &{ \chi } &{ 0 }\\{0 } &{ 0 } &{ \lambda_{ {\rm{r}}}}\end{matrix}\right),\; \hat{{\boldsymbol{\varGamma}}} := \left(\begin{matrix}{0 } &{ \dfrac{\omega_3}{2} } &{ -\beta }\\{-\dfrac{\omega_3}{2} } &{ 0 } &{ \alpha }\\{0 } &{ 0 } &{ 0}\end{matrix}\right). $$ (42)

    (b) Quasiorthogonal-nonnormal:

    $$ \hat{{\boldsymbol{N}}} := \left(\begin{matrix} {\chi } &{ \dfrac{\omega_3}{2} } &{ 0} \\{-\dfrac{\omega_3}{2} } &{ \chi } &{ 0 }\\{0 } &{ 0 } &{ \lambda_{ {\rm{r}}} }\end{matrix}\right),\; \hat{{\boldsymbol{\varGamma}}} := \left(\begin{matrix} {0 } &{ \dfrac{\gamma}{2} } &{ -\beta} \\{\dfrac{\gamma}{2} } &{ 0 } &{ \alpha }\\{0 } &{ 0 } &{ 0}\end{matrix}\right). $$ (43)

    Remark 3 When $ {\boldsymbol{A}} $ has three real eigenvalues, eq.(27) is not a real Schur form of $ {\boldsymbol{A}} $, but it is still a uniquely determined canonical form of $ {\boldsymbol{A}} $. In this case,

    $$ \chi = \frac{\lambda_1 + \lambda_2}{2} $$ (44)

    and $ |\gamma|\ge|\omega_3| $, hence $ \psi\le0 $. The decompositions constructed above based on the canonical form eq.(27) extends to this case.

    Remark 4 Kolář[15-16] defined a triple decomposition of motion (TDM) earlier than NND of Li[1], but the two decompositions have very similar appearances, which has caused some confusion. Here we point out several distinctions between the two:

    1) The purely asymmetric tensor in ref. [15-16], defined by the vanishing of the product of each pair of off-diagonal entries, is generally not nilpotent. For example, the matrix

    $$ \left(\begin{matrix}{ 0 } &{ \gamma } &{ 0} \\{0 } &{ 0 } &{ \alpha }\\{\beta } &{ 0 } &{ 0}\end{matrix}\right) $$ (45)

    is purely asymmetric according to ref. [15-16], but its determinant is nonzero, indicating nonzero eigenvalues. Hence it is not nilpotent. Of course, the purely asymmetric tensor can be nilpotent in several special cases, especially in the 2D case. But the concept behind it is different from the nilpotent tensor. Actually, a condition of nilpotency can be phrased as a given $ n $-by-$ n $ matrix is nilpotent if among any $ n $ entries that occupy $ n $ rows and $ n $ columns, there is a zero.

    2) TDM maximizes the norm of the purely asymmetric tensor through changing of basis, while the NNDs do not maximize the nilpotent tensor, neither do they choose basis according to optimization principles.

    3) Regardless of the difference in their choice of basis, TDM and NND produce different results. For example, suppose the basic reference frame (BRF) of ref. [15-16] is the same as the basis of the special real Schur form eq.(27). When $ {\boldsymbol{A}} $ has complex conjugate eigenvalues, its TDM coincides with eq.(38); but when $ {\boldsymbol{A}} $ has three real eigenvalues, its TDM coincides with eq.(40) instead. Of course, we can define a new type of NND by combining eq.(38) and eq.(40) in each case to include the results of TDM, but it is somewhat artificial and we don’t know whether this is always the case.

    In this section, we show that it is generally impossible to transform a complex NND into a real NND. Consider the general complex Schur form eq.(12). Its complex NND is just the splitting of the diagonal and off-diagonal parts, as given by eq.(35). Let’s find a unitary transformation that converts it to a real NND. Since $ \lambda_3=\lambda_{ {\rm{r}}} $ is real, the required unitary transformation should leave $ {\boldsymbol{e}}_{ {\rm{r}}} $ unchanged, hence is essentially a unitary transformation in the subspace normal to $ {\boldsymbol{e}}_{ {\rm{r}}} $.

    Denote

    $$ {\boldsymbol{\varLambda}}_2 := {\rm{diag}}\,\left({\lambda_1,\lambda_2}\right), $$ (46)

    where $ \lambda_{1,2}=\lambda_{ {\rm{c}} {\rm{r}}}\pm {\rm{i}}\lambda_{ {\rm{c}} {\rm{i}}} $ is a pair of complex conjugate numbers. We are seeking a unitary transformation that converts $ {\boldsymbol{\varLambda}}_2 $ into a real matrix. Let

    $$ {\boldsymbol{R}}_2 := \frac{1}{\sqrt{2}} \left(\begin{matrix}{- {\rm{i}} } &{ 1} \\ { 1 } &{ - {\rm{i}} } \end{matrix}\right). $$ (47)

    It can be verified that $ {\boldsymbol{R}}_2 $ is unitary and

    $$ {\boldsymbol{N}}_2 := {\boldsymbol{R}}_2 {\boldsymbol{\varLambda}}_2 {\boldsymbol{R}}_2^\dagger = \left(\begin{matrix} { \lambda_{ {\rm{c}} {\rm{r}}} } &{ \lambda_{ {\rm{c}} {\rm{i}}}} \\{-\lambda_{ {\rm{c}} {\rm{i}}} } &{ \lambda_{ {\rm{c}} {\rm{r}}} } \end{matrix}\right), $$ (48)

    which is real. Of course, to transform $ {\boldsymbol{\varLambda}}_2 $ into a real matrix, the choice of $ {\boldsymbol{R}}_2 $ is not unique. For every real orthogonal matrix

    $$ {\boldsymbol{Q}}_2 = \left(\begin{matrix}{\cos\,\theta } &{ \sin\,\theta }\\{-\sin\,\theta } &{ \cos\,\theta}\end{matrix}\right), $$ (49)

    $ {\boldsymbol{Q}}_2 {\boldsymbol{R}}_2 $ will also do the job. Let

    $$ {\boldsymbol{R}} := \left(\begin{matrix}{{\boldsymbol{R}}_2 } &{ {\bf{0}} }\\{{\bf{0}} } &{ 1}\end{matrix}\right),\; {\boldsymbol{Q}} := \left(\begin{matrix} { {\boldsymbol{Q}}_2 } &{ {\bf{0}}} \\ { {\bf{0}} } &{ 1 } \end{matrix}\right) \;,$$ (50)

    and $ \tilde{\tilde{{\boldsymbol{A}}}} $ be the general complex Schur form eq.(12), then

    $$ \tilde{\tilde{\tilde{{\boldsymbol{A}}}}} := ( {\boldsymbol{Q}} {\boldsymbol{R}})\tilde{\tilde{{\boldsymbol{A}}}}( {\boldsymbol{Q}} {\boldsymbol{R}})^\dagger = \left(\begin{matrix} { {\boldsymbol{M}} } &{ {\boldsymbol{d}}^\dagger} \\ { {\bf{0}} } &{ \lambda_{ {\rm{r}}} } \end{matrix}\right), $$ (51)

    where

    $$ {\boldsymbol{M}} := \left(\begin{matrix} {\lambda_{ {\rm{c}} {\rm{r}}} - \dfrac{ {\rm{i}}}{2}a_{12}{\rm{e}}^{ {\rm{i}}(2\theta+\phi_{13}-\phi_{23})} }&{ \lambda_{ {\rm{c}} {\rm{i}}} + \dfrac{1}{2}a_{12}{\rm{e}}^{ {\rm{i}}(2\theta+\phi_{13}-\phi_{23})}} \\{-\lambda_{ {\rm{c}} {\rm{i}}} + \dfrac{1}{2}a_{12}{\rm{e}}^{ {\rm{i}}(2\theta+\phi_{13}-\phi_{23})}} & {\lambda_{ {\rm{c}} {\rm{r}}} + \dfrac{ {\rm{i}}}{2}a_{12}{\rm{e}}^{ {\rm{i}}(2\theta+\phi_{13}-\phi_{23})} } \end{matrix}\right) $$ (52)

    and

    $$ {\boldsymbol{d}} := \left({ \frac{ {\rm{i}}}{\sqrt{2}}a_{13}{\rm{e}}^{- {\rm{i}}(\theta+\phi_{13})} + \frac{1}{\sqrt{2}}a_{23}{\rm{e}}^{ {\rm{i}}(\theta-\phi_{23})}, \frac{1}{\sqrt{2}}a_{13}{\rm{e}}^{- {\rm{i}}(\theta+\phi_{13})} + \frac{ {\rm{i}}}{\sqrt{2}}a_{23}{\rm{e}}^{ {\rm{i}}(\theta-\phi_{23})}}\right). $$ (53)

    That’s all we can do. Now we need to determine $ \phi_{13},\phi_{23} $ and $ \theta $ such that $ \tilde{\tilde{\tilde{{\boldsymbol{A}}}}} $ is real. We rewrite $ {\boldsymbol{M}} = {\boldsymbol{N}}_2 + \xi {\boldsymbol{L}}_2 $, where $ \xi:=\dfrac{1}{2}a_{12}{\rm{e}}^{ {\rm{i}}(2\theta+\phi_{13}-\phi_{23})} $ and

    $$ {\boldsymbol{L}}_2 := \left(\begin{matrix} { - {\rm{i}} }& {1} \\ { 1 }& { {\rm{i}}} \end{matrix}\right). $$ (54)

    However, whatever $ \phi_{13},\;\phi_{23} $ and $ \theta $ are chosen, $ {\boldsymbol{M}} $ can never be real, so can $ \tilde{\tilde{\tilde{{\boldsymbol{A}}}}} $. Therefore, a complex NND can not be converted into a real NND by unitary transformations in general. They are intrinsically distinct decompositions.

    The major contributions of this article are:

    1) The general complex and real Schur forms are derived constructively, with special consideration for their non-uniqueness. Conditions of uniqueness are found to obtain standardized Schur forms.

    2) Conditions for normality and nilpotency are presented. Based on the general complex Schur form, a complex NND is constructed. Based on a special real Schur form, several real NND as well as normal-nonnormal decompositions are constructed. A comparison of advantages and disadvantages between complex and real NNDs is presented.

    3) Several confusing points are explained, including the distinction between NND and the TDM proposed by Kolář[15-16], the intrinsic gap between complex and real NNDs.

    We hope that the objective of the article, i.e., clarifying the confusion about Schur forms and NNDs existing in research, has been achieved.

    Now the open question is, provided so many versions of decompositions, which is the right one to use? Typically, there are two widely adopted usages:

    1) Use complex Schur form and complex NND throughout.

    2) Use the special block real Schur form eq.(27) and correspongding NND when $ {\boldsymbol{A}} $ has complex conjugate eigenvalues; use the real Schur form eq.(13) and the corresponding NND when $ {\boldsymbol{A}} $ has three real eigenvalues.

    But now one has an additional option: use the special block real Schur form eq.(27) and the corresponding NND throughout, i.e., for both real and complex eigenvalue cases. Of course, the above description didn’t take into account the different variants of real NNDs as proposed in section 3, which provides even more options. These options will serve different needs.

    The author would like to thank professors WU Jiezhi and YAO Jie for inspiring communication. But the author is solely responsible for the content and point of view of this communication.

  • [1] LI Z. Theoretical study on the definition of vortex[D]. Beijing: Tsinghua University, 2010.
    [2] LI Z, ZHANG X, HE F. Evaluation of vortex criteria by virtue of the quadruple decomposition of velocity gradient tensor[J]. Acta Physica Sinica, 2014, 63(5): 054704. doi: 10.7498/aps.63.054704
    [3] MURNAGHAN F D, WINTNER A. A canonical form for real matrices under orthogonal transformations[J]. Proceedings of the National Academy of Sciences, 1931, 17(7): 417-420. doi: 10.1073/pnas.17.7.417
    [4] SCHUR I. Über die charakteristischen wurzeln einer linearen substitution mit einer anwendung auf die theorie der integralgleichungen[J]. Mathematische Annalen, 1909, 66(4): 488-510. doi: 10.1007/BF01450045
    [5] STICKELBERGER L. Ueber Reelle Orthogonale Substitutionen[M]. Bern: Schweizerische Nationalbibliothek, 1877.
    [6] LÉON A M. Sur l’ Hermitien[J]. Rendicontidel Circolo Matematico di Palermo (1884—1940), 1907, 16: 104-128.
    [7] KEYLOCK C J. Synthetic velocity gradient tensors and the identification of statistically significant aspects of the structure of turbulence[J]. Physical Review Fluids, 2017, 2(8): 084607. doi: 10.1103/PhysRevFluids.2.084607
    [8] KEYLOCK C J. The Schur decomposition of the velocity gradient tensor for turbulent flows[J]. Journal of Fluid Mechanics, 2018, 848: 876-905. doi: 10.1017/jfm.2018.344
    [9] DAS R, GIRIMAJI S S. Revisiting turbulence small-scale behavior using velocity gradient triple decomposition[J]. New Journal of Physics, 2020, 22(6): 063015. doi: 10.1088/1367-2630/ab8ab2
    [10] HOFFMAN J. Energy stability analysis of turbulent incompressible flowbased on the triple decomposition of the velocity gradient tensor[J]. Physics of Fluids, 2021, 33(8): 081707. doi: 10.1063/5.0060584
    [11] ZHU J Z. Thermodynamic and vortic structures of real Schur flows[J]. Journal of Mathematical Physics, 2021, 62(8): 083101. doi: 10.1063/5.0052296
    [12] ZHU J Z. Compressible helical turbulence: fastened-structure geometry and statistics[J]. Physics of Plasmas, 2021, 28(3): 032302. doi: 10.1063/5.0031108
    [13] KRONBORG J, SVELANDER F, ERIKSSON-LIDBRINK S, et al. Computational analysis of flow structures in turbulent ventricular blood flow associated with mitral valve intervention[J]. Frontiers in Physiology, 2022, 13: 806534. doi: 10.3389/fphys.2022.806534
    [14] ARUN R, COLONIUS T. Velocity gradient analysis of a head-on vortexring collision[J]. Journal of Fluid Mechanics, 2024, 982: A16. doi: 10.1017/jfm.2024.90
    [15] KOLÁŘ V. 2D Velocity-field analysis using triple decomposition of motion[C]//Proceedings of the Fifteenth Australasian Fluid Mechanics Conference. Sydney, Australia, 2004.
    [16] KOLÁŘ V. Vortex identification: new requirements and limitations[J]. International Journal of Heat and Fluid Flow, 2007, 28(4): 638-652. doi: 10.1016/j.ijheatfluidflow.2007.03.004
    [17] KRONBORG J, HOFFMAN J. The triple decomposition of the velocity gradient tensor as a standardized real Schur form[J]. Physics of Fluids, 2023, 35: 031703. doi: 10.1063/5.0138180
    [18] ZOU W, XU X Y, TANG C X. Spiral streamline pattern around a critical point: its dual directivity and effective characterization by right eigen representation[J]. Physics of Fluids, 2021, 33(6): 067102. doi: 10.1063/5.0050555
  • 加载中
计量
  • 文章访问数:  214
  • HTML全文浏览量:  99
  • PDF下载量:  40
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-05-08
  • 修回日期:  2024-07-03
  • 刊出日期:  2024-09-01

目录

/

返回文章
返回