Loading [MathJax]/jax/output/SVG/jax.js
Theory article

Cumulative STF coefficients evaluation and validation

  • Received: 22 October 2024 Revised: 02 December 2024 Accepted: 10 December 2024 Published: 20 December 2024
  • Decomposition of geophysical functions in ranks on degrees of components of a single position vector with coefficients in the form of the indexes of tensors, symmetric and traceless on any couple (symmetric and trace free [STF] tensors or deviators), is applied along with decomposition on surface harmonics (scalar, vector, and tensor). The article considers the problem of deviator decomposition of a function having the special form of a series of degrees of components of a unit radius vector. The algorithm evaluation of STF coefficients using known values of series coefficients is under consideration. Taking into account that often only the first several of these coefficients are used, the author created and presented a table with several coefficient formulas for reference and validation. The STF-formalism is mainly used for the representation of radiative gravity fields and gravitational waves in general relativity; however, it can also be applied in mathematical physics to represent spherical harmonics, including fluid dynamics in Earth's outer core and seismic wave analysis.

    Citation: Pasynok Sergey. Cumulative STF coefficients evaluation and validation[J]. Metascience in Aerospace, 2024, 1(4): 371-378. doi: 10.3934/mina.2024017

    Related Papers:

    [1] Yuan Hu, Zhihao Jiang, Xintai Yuan, Xifan Hua, Wei Liu . Isometric mapping algorithm based GNSS-R sea ice detection. Metascience in Aerospace, 2024, 1(1): 38-52. doi: 10.3934/mina.2024002
    [2] Sarker Ashraful Islam, Farhana Kabir Esheta, Md Mahir Shahriar, Dewan Hasan Ahmed . Numerical study of aerodynamic drag reduction of a circular cylinder with an inbuilt nozzle. Metascience in Aerospace, 2024, 1(4): 379-400. doi: 10.3934/mina.2024018
    [3] Shorob Alam Bhuiyan, Ikram Hossain, Redwan Hossain, Md. Sakib Ibn Mobarak Abir, Dewan Hasan Ahmed . Effect of a bioinspired upstream extended surface profile on flow characteristics and a drag coefficient of a circular cylinder. Metascience in Aerospace, 2024, 1(2): 130-158. doi: 10.3934/mina.2024006
    [4] Victor Cánepa, Pablo Servidia . Design and evaluation of INS/GNSS loose and tight coupling applied to launch vehicle integrated navigation. Metascience in Aerospace, 2024, 1(1): 66-109. doi: 10.3934/mina.2024004
    [5] Wenchao Cai, Yadong Zhou . Preventing aircraft from wildlife strikes using trajectory planning based on the enhanced artificial potential field approach. Metascience in Aerospace, 2024, 1(2): 219-245. doi: 10.3934/mina.2024010
    [6] Igor V. Bezmenov . Fast algorithm for cleaning highly noisy measurement data from outliers, based on the search for the optimal solution with the minimum number of rejected measurement data. Metascience in Aerospace, 2024, 1(1): 110-129. doi: 10.3934/mina.2024005
    [7] Weizhen Tang, Jie Dai, Zhousheng Huang . Comprehensive assessment of air pollutant and noise emissions at airports across different altitudes. Metascience in Aerospace, 2024, 1(3): 292-308. doi: 10.3934/mina.2024013
    [8] Dulce M Graciano, Fernando Z Sierra-Espinosa, Juan C García . Numerical simulation of vortex-induced vibration in a bladeless turbine: Effects of separation distance between tandem harvesters. Metascience in Aerospace, 2024, 1(3): 309-328. doi: 10.3934/mina.2024014
    [9] Jubilee Prasad Rao, Sai Naveen Chimata . Machine learning-based surrogates for eVTOL performance prediction and design optimization. Metascience in Aerospace, 2024, 1(3): 246-267. doi: 10.3934/mina.2024011
    [10] James H. Page, Lorenzo Dorbolò, Marco Pretto, Alessandro Zanon, Pietro Giannattasio, Michele De Gennaro . Sensitivity analysis of mixed analysis-synthesis flight profile reconstruction. Metascience in Aerospace, 2024, 1(4): 401-415. doi: 10.3934/mina.2024019
  • Decomposition of geophysical functions in ranks on degrees of components of a single position vector with coefficients in the form of the indexes of tensors, symmetric and traceless on any couple (symmetric and trace free [STF] tensors or deviators), is applied along with decomposition on surface harmonics (scalar, vector, and tensor). The article considers the problem of deviator decomposition of a function having the special form of a series of degrees of components of a unit radius vector. The algorithm evaluation of STF coefficients using known values of series coefficients is under consideration. Taking into account that often only the first several of these coefficients are used, the author created and presented a table with several coefficient formulas for reference and validation. The STF-formalism is mainly used for the representation of radiative gravity fields and gravitational waves in general relativity; however, it can also be applied in mathematical physics to represent spherical harmonics, including fluid dynamics in Earth's outer core and seismic wave analysis.



    Decomposition of geophysical functions in ranks on degrees of components of a single position vector with coefficients in the form of the indexes of tensors, symmetric and traceless on any couple (symmetric and trace free [STF] tensors or deviators), is applied along with decomposition on surface harmonics (scalar, vector, tensor). The term deviator comes from the elasticity theory [1]. From the point of view of group theory, the deviator part of a tensor is the first member of decomposition of tensor on irreducible representations of rotation group SO (dim).

    In group theory, it has been proven that deviator decompositions are equivalent to decomposition on surface harmonics (scalar, vector, and tensor), which makes it possible to use such decompositions to solve mathematical physics equations. At the same time, in some cases (e.g., boundary conditions), there is a need to represent the sum of the inner products of symmetric tensors ˜FM with multi components of a unit radial vector:

    F=Nm=1˜FMnM (1)

    in the form of the sum of the inner products of deviators ˆFM:

    F=Nm=1ˆFMnM, (2)

    where nMni1...niM is the multi component of a unit radial vector in Damour's notation [2]. The multi index in Thorn notation: TIlTi1i2...il [3, formulae (1.6a)] will be used in section 2.

    This task can be solved directly by evaluating the ˆFM coefficient with the integration of sum (1) multiplied by the STF-basis tensor of rank M by a unit sphere. In this paper, an easier route is considered.

    Let us consider, for example, case N = 2:

    F=2m=0˜FMnM=˜Fijninj+˜Fini+˜F0.

    The Einstein summation notation is used whenever we have an expression with a repeated index (multi index). Thus, we implicitly know to sum over that index (indexes of multi index) from 1 to dim, where dim is notation of the dimension of the space. The dimension of space (dim) can be also written as δkk=1+1+...+1dim=dim.

    One can obtain the following for STF part of tensor with rank 2: ˜Fij=ˆFij+1dimδijFkk. Substituting this into an expression for F and introducing similar ones, one can get:

    F=2m=0˜FMnM=(ˆFij+1dimδijFkk)ninj+˜Fini+˜F0=ˆFijninj+ˆFini+ˆF0,

    where deviator coefficients are equal:

    ˆFij=˜Fij1dimδijFkk,ˆFk=˜Fk,ˆF0=˜F0+1dimFkk.

    In this simple case, formulas for finding deviator coefficients from symmetric coefficients are very easy, and it is not needed in some special evaluations. However, when the number of sum members N is growing, the difficulties increase because the fragments of the members of the senior ranks are "poured" into the members of the lower ranks and the more members there are in total, the more complex they are to enter there.

    The algorithm of estimation of coefficients ˆFK of sum (2) using knowing values of symmetric coefficients ˜FJ of sum (1) is under consideration. Taking into account that often only the first several of these coefficients are used, the author created and presented a table with several coefficient formulas for reference and validation.

    The algebraic algorithm for cumulative STF coefficients estimation was proposed and proved by author in [4] and is represented here in Supplementary in the English translation. It can be formulated as the following theorem about cumulative STF coefficients.

    Theorem about cumulative STF coefficients. Sum (1) can be represented as exactly equal to its sum (2) with the deviator coefficients estimated according to formulas:

    ˆF2J=[N2]s=j˜a(2s,sj)˜F<2J>KSJKSJ, ˆF2J+1=[N12]s=j˜a(2s+1,sj)˜F<2J+1>KSJKSJ, (3)

    where ˜a(n,l) are cumulative STF coefficients and calculated by values coefficients a(n,l) of STF part of tensor by following the recurrent algorithm.

    1) for l = 0, 1, 2 cumulative coefficients are evaluated by formulas:

    ˜a(n,0)=1,˜a(n,1)=a(n,1),˜a(n,2)=a(n,1)a(n2,1)a(n,2). (4)

    2) for l > 2:

    – evaluated a(n,l2,1) by formula:

    a(n,l2,1)=a(n,1)a(n2,l1)a(n,l) (5)

    – for k = 1, …, l-2 evaluated:

    a(n,l2k,k+1)=a(n,l1k,k)a(n2(k+1),l1k)˜a(n,k+1) (6)

    Last coefficient a(n,0,l1) is equal to ˜a(n,l).

    Remark. At the same time, the coefficients ˜a(n,k+1) for k = 1, …, l-2 have to be evaluated on previous steps.

    The formulae for the coefficients a(n,l) of STF part of tensor for 3-dimentional space one can be found in [1, formulae (2.2c)]. For multidimensional space, it can be found in [4, formulae (4)].

    Taking into account that often only the first several of these coefficients are used, the algebraic algorithm for cumulative STF coefficients estimation was used by author for the evaluation the first several cumulative STF coefficients (see Supplementary). For reference, author prepared the following tables.

    The following procedure was used for the validation of formulas for cumulative STF coefficients presented at Table 1. At first, the STF coefficients of sum (2) must be expressed in terms of sum (1) for given sum limit N, according to formulae (3) with cumulative STF coefficient from Table 1. After that, these expressions must be inserted in sum (2) and similar ones have to be introduced. If the sum (1) will be obtained, then the validation is assumed to be successful. Validation for N = 2 (see example in Introduction) obviously is successful.

    Table 1.  The cumulative STF coefficients ˜a(n,l) for n from 1 to 7.
    n L
    0 1 2 3
    0 1
    1 1
    2 1 1dim
    3 1 3(dim+2)
    4 1 6(4+dim) 3(2+dim)dim
    5 1 10(6+dim) 15(4+dim)(2+dim)
    6 1 15(8+dim) 45(6+dim)(4+dim) 15(4+dim)(2+dim)dim
    7 1 21(10+dim) 105(8+dim)(6+dim) 105(6+dim)(4+dim)(2+dim)

     | Show Table
    DownLoad: CSV

    Formulas (3) for N = 3 can be written as:

    ˆF0=˜F0+˜a(2,1)˜Fss=˜F0+1dim˜Fss,ˆFi=˜Fi+˜a(3,1)˜Fiss=˜Fi+3(2+dim)˜Fiss,ˆFij=˜Fij=˜Fij1dimδij˜Fss,ˆFijk=˜Fijk=˜Fijk3(2+dim)13(δij˜Fssk+δik˜Fssj+δjk˜Fssi).

    Now, inserting these in (2) and introducing similar ones, one can obtain:

    F=˜F0+{1dim˜Fss}+˜Fini+[3(2+dim)˜Fissni]+˜Fijninj+{1dimninjδij˜Fss}+˜Fijkninjnk++[3(2+dim)13(δij˜Fkss+δik˜Fjss+δjk˜Fiss)ninjnk]=˜F0+˜Fini+˜Fijninj+˜Fijkninjnk.

    Thus, validation is successful. The shrinking members are enclosed in big curly and square brackets.

    Formulas (3) for N = 4 can be written as:

    ˆF0=˜F0+˜a(2,1)˜Fss+˜a(4,2)˜Fssmm=˜F0+1δkk˜Fss+{3(2+dim)dim˜Fssmm},
    ˆFi=˜Fi+˜a(3,1)˜Fiss=˜Fi+3(2+dim)˜Fiss,
    ˆFij=˜Fij+˜a(4,1)˜Fijss=˜Fij1dimδij˜Fss+[6(4+dim)˜Fijss]+{6(4+dim)1dimδij˜Fmmss},
    ˆFijk=˜Fijk=˜Fijk3(2+dim)13(δij˜Fssk+δik˜Fssj+δjk˜Fssi),
    ˆFijkl=˜Fijkl=˜Fijkl+[6(4+dim)16(δij˜Fsskl+δik˜Fssjl+δil˜Fssjk+δjk˜Fssil+δjl˜Fssik+δlk˜Fssij)]++{3(4+dim)(2+dim)13(δijδkl+δikδjl+δilδjk)˜Fssmm}.

    As shown above, the members which not enclosed in big brackets, give us the first 4 members of sum (1), excluding one with coefficient ˜Fijkl. Thus, only the new members with coefficients enclosed in big brackets and ˜Fijkl have to be checked. The members of sum (2) of these must give the result ˜Fijklninjnknk. After evaluating these members and introducing similar ones, one can obtain ˜Fijklninjnknk. So, validation is successful.

    Remark. In paper [4], the example for N = 4 has some mistakes because of a big rush before publication; however, all formulas here are correct.

    Similar to the previous case, it is enough to check only for new odd-ranked members:

    ˆFi=˜Fi+˜a(3,1)˜Fiss+˜a(5,2)˜Fisskk=˜Fi+3(2+dim)˜Fiss+{15(4+dim)(2+dim)˜Fisskk},
    ˆFijk=˜Fijk+˜a(5,1)˜Fijkss=˜Fijk3(2+dim)13(δij˜Fssk+δik˜Fssj+δjk˜Fssi)++[10(6+dim)˜Fijkss]+{10(6+dim)3(2+dim)13(δij˜Fsswwk+δik˜Fsswwj+δjk˜Fwwssi)},
    ˆFijklm=˜Fijklm=˜Fijklm[10(6+dim)110(δij˜Fssklm+δik˜Fssjlm+δil˜Fssjkm+δjk˜Fssilm+δjl˜Fssikm+δlk˜Fssijm+δim˜Fssklj+δkm˜Fssijl+δlm˜Fssijk+δjm˜Fssilk)]++{15(6+dim)(4+dim)115(δijδkl˜Fsswwm+δijδkm˜Fsswwl+δijδlm˜Fsswwk+δikδjl˜Fsswwm+δikδjm˜Fsswwl+δikδlm˜Fsswwj+δilδjk˜Fsswwm++δilδjm˜Fsswwk+δilδkm˜Fsswwj+δimδjk˜Fsswwl+δimδjl˜Fsswwk+δimδkl˜Fsswwj+δjkδlm˜Fsswwi+δjlδkm˜Fsswwi+δjmδkl˜Fsswwi)}

    Other members are the same as in previous case. After inserting the coefficients in (2) and introducing similar ones, one can obtain ˜Fijklmninjnknknm for the sum members with coefficients enclosed in big brackets and ˜Fijklm. So, validation is successful.

    Similar to the previous case, it is enough to check only for new even-ranked members:

    ˆF0=˜F0+˜a(2,1)˜Fss+˜a(4,2)˜Fssmm+˜a(6,3)˜Fssmm=˜F0+1dim˜Fss+3(2+dim)dim˜Fssmm++{15(4+dim)(2+dim)dim˜Fssmmrr},
    ˆFij=˜Fij+˜a(4,1)˜Fijss+˜a(6,2)˜Fijssrr=˜Fij1dimδij˜Fss+6(4+dim)(˜Fijss1dimδij˜Fmmss)++[45(6+dim)(4+dim)˜Fijssrr]+{45(6+dim)(4+dim)1dimδij˜Fmmssrr},
    ˆFijkl=˜Fijkl+˜a(6,1)˜Fijklss=˜Fijkl6(4+dim)16(δij˜Fsskl+δik˜Fssjl+δil˜Fssjk+δjk˜Fssil+δjl˜Fssik+δlk˜Fssij)++3(4+dim)(2+dim)13(δijδkl+δikδjl+δilδjk)˜Fssmm+(15(8+dim)˜Fijklss)+[15(8+dim)6(4+dim)16(δij˜Fssklrr++δik˜Fssjlrr+δil˜Fssrrjk+δjk˜Frrssil+δjl˜Frrssik+δlk˜Frrssij)]+{15(8+dim)3(4+dim)(2+dim)13(δijδkl+δikδjl+δilδjk)˜Fssmmrr},
    ˆFijklmp=˜Fijklmp=˜Fijklmp+(15(8+dim)115(δij˜Fssklmp++δjm˜Fssilk))+[45(8+dim)(6+dim)145(δijδkl˜Fsswwmp+)]++{15(8+dim)(6+dim)(4+dim)115(δijδklδmp+)˜Fssmmrr}.

    Other members are the same as in previous case. The ellipsis denotes similar terms, the form of which is understandable from context. After inserting the coefficients in (2) and introducing the similar ones, one can obtain ˜Fijklmpninjnknknmnp for sum of members with coefficients enclosed in big brackets and ˜Fijklmp. So, validation is successful.

    Similar to the previous case, it is enough to check only for new odd-ranked members:

    ˆFi=˜Fi+˜a(3,1)˜Fiss+˜a(5,2)˜Fisskk+˜a(7,3)˜Fisskkrr=˜Fi+3(2+dim)˜Fiss+15(4+dim)(2+dim)˜Fisskk++{105(6+dim)(4+dim)(2+dim)˜Fisskkrr},
    ˆFijk=˜Fijk+˜a(5,1)˜Fijkss+˜a(7,2)˜Fijkssrr=˜Fijk3(2+dim)13(δij˜Fssk+δik˜Fssj+δjk˜Fssi)++10(6+dim)(˜Fijkss3(2+dim)13(δij˜Fsswwk+δik˜Fsswwj+δjk˜Fwwssi))++[105(8+dim)(6+dim)˜Fijkssrr]+{105(8+dim)(6+dim)3(2+dim)13(δij˜Fsswwrrk+δik˜Fsswwrrj+δjk˜Fwwssrri)},
    ˆFijklm=˜Fijklm+˜a(7,1)˜Fijklmss=˜Fijklm10(6+dim)110(δij˜Fssklm+...+δjm˜Fssilk)+15(6+dim)(4+dim)115(δijδkl˜Fsswwm++...+δjmδkl˜Fsswwi)+(21(10+dim)˜Fijklmrr)+[21(10+dim)10(6+dim)110(δij˜Fsswwklm++δjm˜Fsswwilk)]+{15(6+dim)(4+dim)115(δijδkl˜Fsswwrrm++δjmδkl˜Fsswwi))},
    ˆFijklmps=˜Fijklmps=˜Fijklmps+(21(10+dim)121(δij˜Fwwklmps++δjm˜Fwwilks))++[105(10+dim)(8+dim)1105(δijδkl˜Frrwwmps+)]+{105(10+dim)(8+dim)(6+dim)1105(δijδklδmp˜Fwwmmrrs+)},

    Other members are the same as in previous case. The ellipsis denotes similar terms, the form of which is understandable from context. After inserting coefficients in (2) and introducing similar ones, one can obtain ˜Fijklmpsninjnknknmnpns for sum of members with coefficients enclosed in big brackets and ˜Fijklmps. So, validation is successful.

    The cumulative STF coefficients formulas were obtained based on algebraic algorithm [4] for space with finite dimension. The formulas were checked analytically, and validation was successful. The formulas of Table 1 can be used for easier transition from sum in form (1) to sum in form (2) compared with directly using the theorem about cumulative STF coefficients algorithm. In conclusion, it should be noted that STF-formalism can be used for the solution of all mathematical physics tasks that include spherical harmonics. For example, it was used by author for calculating the added masses tensor in the solution of the Earth's solid inner core motion [5].

    The author declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    The author would like to thank of Journal Metascience in Aerospace Editorial office for possibility of this paper publication.

    The authors declare no conflict of interest.



    [1] Kim NH (2015) Introduction to nonlinear finite element analysis, Springer Science+ Business Media, New York. https://doi.org/10.1007/978-1-4419-1746-1
    [2] Damour T, Iyer BR (1991) Multipole analysis for electromagnetism and linearized gravity with irreducible Cartesian tensors. J Phys Rev D 43: 3259. https://doi.org/10.1103/PhysRevD.43.3259 doi: 10.1103/PhysRevD.43.3259
    [3] Thorne KS (1980) Multipole expansions of gravitational radiation. Rev Mod Phys 52: 299–338. https://doi.org/10.1103/revmodphys.52.299 doi: 10.1103/revmodphys.52.299
    [4] Pasynok SL (2018) The algebraic algorithm of decomposition on deviators of functions in the form of the sum of terms with the symmetric coefficients, Izvestiya GAO RAN 225: 267,272, ISSN 0367-7966 (in russian).
    [5] Pasynok SL (1999) Uchyot tenzora prisoedinyonnyh mass v zadache o dvizhenii tvyordogo yadra Zemli, Vestnik Moskovskogo universiteta. Seriya 3. Fizika. Astronomiya. 1: 46–49.
  • mina-01-04-017-s001.pdf
  • Reader Comments
  • © 2024 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(333) PDF downloads(10) Cited by(0)

Figures and Tables

Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog