# Analysis of Steady Vortex Rings Using Contour Dynamics Method for Fluid Velocity

## Article information

J. Ocean Eng. Technol. 2022;36(2):108-114
Publication date (electronic) : 2022 March 31
doi : https://doi.org/10.26748/KSOE.2021.070
1Professor, School of Naval Architecture and Ocean Engineering, University of Ulsan, Ulsan, Korea
Corresponding author Yoon-Rak Choi: +82-52-259-2158, yrchoi@ulsan.ac.kr
Received 2021 September 23; Revised 2022 January 6; Accepted 2022 February 12.

## Abstract

Most studies on the shape of the steady vortex ring have been based on the Stokes stream function approach. In this study, the velocity approach is introduced as a trial approach. A contour dynamics method for fluid velocity is used to analyze the Norbury–Fraenkel family of vortex rings. Analytic integration is performed over the logarithmic-singular segment. A system of nonlinear equations for the discretized shape of the vortex core is formulated using the material boundary condition of the core. An additional condition for the velocities of the vortical and impulse centers is introduced to complete the system of equations. Numerical solutions are successfully obtained for the system of nonlinear equations using the iterative scheme. Specifically, the evaluation of the kinetic energy in terms of line integrals is examined closely. The results of the proposed method are compared with those of the stream function approaches. The results show good agreement, and thereby, confirm the validity of the proposed method.

## 1. Introduction

A water-jet can be used as one of the propulsion systems for ships and marine life. When a jet is injected to obtain thrust, a vortex ring is formed at a nozzle and then propagated downstream (Krueger et al., 2008). Furthermore, a vortex ring is generated due to volcanic eruption or nuclear explosion (Akhmetov, 2009).

The flow of a vortex ring is formulated with the Helmholtz vorticity equation in inviscid and incompressible fluids (Batchelor, 1967). A steady vortex ring was first reported by Helmholtz (1867) who examined a vortex ring of an small circular cross section, while a spherical vortex was first analyzed by Hill (1894). Norbury (1973) analyzed a vortex ring in a steady state for general circumstances, which is referred to as the Norbury–Fraenkel family (N-F family) of vortex rings.

A dynamic analysis is required for analyzing the instability due to the disturbance or interaction between vortex rings. A contour dynamics (CD) method for fluid velocity is used for analyzing the complex evolution of the contour of a vortex core. The CD method is a two-dimensional or axisymmetric flow analysis method due to the isolated vorticity in an inviscid, incompressible, and irrotational flow field (Pullin, 1992; Smith et al., 2018). The CD method can drastically reduce the burden of computations because the computation is performed in the form of line integrals on the boundary contour of the vorticity region. The fluid velocity on the contour is calculated using the CD method and then applied with time integrals to estimate the dynamic changes in the shape of the vortex core. Zabusky et al. (1979) introduced the CD method in dynamic analysis of two-dimensional vortex patches. Various examples of dynamic analysis for three- dimensional axisymmetric vortex rings are provided in the study by Shariff et al. (1989).

In this study, the CD method was applied to the analysis of the N-F family of vortex rings which are flows in steady state. Choi (2020) combined the CD method for a stream function (Shariff et al., 1989) and the direct shape-calculation method, and thus obtained results that were superior that those reported by Norbury (1973) wherein surface integrals and Fourier analysis were used. As a follow-up study to Choi (2020), in this study, we analyzed the N-F family of vortex rings using the CD method for fluid velocity examined in studies by Shariff et al. (1989) and Shariff et al. (2008). A stream function has been mostly used for analyzing a vortex ring in a steady state (Batchelor, 1967; Fraenkel, 1970; Fraenkel, 1972; Norbury 1973). In this study, we examined whether the CD method for fluid velocity, which is used in dynamic analysis, can also be applied to the analysis of a vortex ring in a steady state.

By applying the boundary conditions of a vortex core to the fluid velocity determined via contour integral, we proposed nonlinear simultaneous equations for the forward speed of a vortex ring and core shape nodal points. The shape nodal points and forward speed were determined via iterative calculations by applying an additional conditional equation, which specified that the speed of vortex center and impulse center were the same. The characteristics of core shapes and vortex rings proved that the analytical method used in this study is valid based on the comparison with previous results.

## 2. Problem Formulation

The fluid is assumed as inviscid and incompressible. In Fig. 1, the core of a vortex ring is analyzed, which is a rotational flow present in an irrotational and infinite flow field. When cylindrical polar coordinates, as shown in the figure, are introduced for analyzing an axisymmetric flow for the x-axis, the fluid velocity vector u⃗ and vorticity vector ω are as follows:

(1) u=uxe^x+uσe^σ
(2) ω=×u=ωϕe^ϕ Core of vortex ring and cylindrical polar coordinates

For an axisymmetric flow, the Helmholtz vorticity equation can be expressed as follows (Batchelor, 1967):

(3) DDt(ωϕσ)=0

Based on Eq. (3), vorticity ωϕ is proportional to σ inside the core, and it is expressed in the equation below. Whereas it is an irrotational flow outside the core.

(4) ωϕ={ΩσinsideD0outsideD

, where ∂D denotes the core boundary, and Ω is a constant. The aforementioned formulation is also reported in Choi (2020).

The fluid velocity vector u⃗ can be calculated using the Green’s 2nd identiy (Shariff et al., 1989; Shariff et al., 2008) as follows:

(5) u(x)=-14π∫∫∫D2u(x)|x-x|dV=14π∫∫∫D×ω(x)|x-x|dV

However, the curl of vorticity vector in Eq. (5) exhibits the following behavior due to the jump in vorticity across the core boundary ∂D in Eq. (4).

(6) ×ω(x)={2Ωe^xinsideDδ-functionbehavioracrossD

Therefore, u⃗ in Eq. (5) can be expressed as the sum of the integrated value inside the core u⃗C and integrated value in the small region crossing the boundary u⃗J as follows:

(7) u=uC+uJ

According to Shariff et al. (1989) and Shariff et al.(2008), u⃗C and u⃗J are determined based on the following equations, which are expressed as contour integrals.

(8) uC=Ωc[G(s){-(x-x)cosθ+σsinθ}-H(s)σsinθ]e^xds
(9) uj=Ωcσ[-G(s)sinθe^x+H(s)cosθe^σ]ds
(10) G(s)=σ4π02π1Δdϕ=12πσσkK(k)
(11) H(s)=σ4π02πcosϕΔdϕ=12πσσ[(2k-k)K(k)-2kE(k)]Δ=(x-x)2+σ2+σ2-2σσcosϕ
, where s, θ, and path c are defined in Fig. 2 (Norbury, 1973; Choi, 2020). Point C in the figure is the mid-point between point A and point B. Furthermore, K(k) and E(k) denote complete elliptic integrals of the first and second kinds, respectively (Gradshteyn and Ryzhik, 2000), and modulus is defined as follows:
(12) k=2σσ(x-x)2+(σ+σ)2
ux and uσ are calculated using Eqs. (7)(9) as follows:
(13) ux=-Ωc(x-x)G(s)dσ+ΩσcH(s)dx
(14) uσ=ΩcσH(s)dσ

When (x, σ)∈c in Eqs. (13) and (14), the movement of a vortex ring can be analyzed. In this study, the N-F family of vortex rings, which are in a steady state, were analyzed. In this case, a vortex ring advances forward in the positive x-direction at a constant speed U. Given that the core boundary c exhibits the characteristics of a material boundary, it should satisfy the following impermeable boundary condition.

(15) ue^n=uxcosθ+uσsinθ=Ucosθfor (x,σ)on c

To define the size of a core, Norbury (1973) expressed the cross-sectional area AΣ of the core using the ring radius L in Fig. 2 and non-dimensionalized mean core radius α as follows:

(16) A=πL2α2

Circulation (Γ), vortical impulse (Γ), and kinetic energy (T) of a inviscid vortex ring are invariant and expressed as follows (Lamb, 1932):

(17) Γ=∫∫ωϕdxdσ=Ω∫∫σdxdσ
(18) P=πρf∫∫ωϕσ2dxdσ=πρfΩ∫∫σ3dxdσ
(19) T=πρf∫∫ψωϕdxdσ=πρfΩ∫∫ψσdxdσ
, where ρf denotes fluid density, and ψ denotes Stokes’ stream function.
(20) ux=1σψσ,uσ=-1σψx

Furthermore, forward speed U can be derived by determining vorticity center (Γ ) and impulse center (P) using Eq. (17) and (18), respectively, and then taking their time derivatives (Lamb, 1932). By considering circulation and impulse as invariants and using the Reynolds transport theorem, U can be expressed as follows:

(21) x¯Γ=ΩΓ∫∫σxdxdσ
(22) U=dx¯Γdt=ΩΓcxue^nds=ΩΓ(cuxσxdσ-cuxσxdx)
(23) x¯P=πρfΩP∫∫σ3xdxdσ
(24) U=dxP¯dt=πρfΩPcσ3xue^nds=πρfΩP(cuxσ3xdσ-cuσσ3xdx)

In this study, physical quantities were non-dimensionalized as the methods in studies by Norbury (1973) and Choi (2020).

(25) [x,σ,s,A,V]=[Lx˜,Lσ˜,Ls˜,L2A˜,L3V˜]
(26) [t,ωϕ,U,ψ]=[t˜/(ΩLα2),ΩLα2ωϕ˜,ΩL2α2U˜,ΩL4α2ψ˜]
(27) [Γ,P,T]=[ΩL3α2Γ˜,ρfΩL5α2P˜,ρfΩ2L7α4T˜]
, where VΣ denotes the core volume, and symbol (∼) denotes non-dimensionalized physical quantities. For the convenience of expression, symbol (∼) is omitted henceforth when expressing non-dimensionalized physical quantities.

Eqs. (13) and (14), which are fluid velocity components, can be non-dimensionalized as follows:

(28) ux=1α2[-c(x-x)G(s)dσ+σcH(s)dx]
(29) uσ=1α2cσH(s)dσ

For expressing contour shapes, radius q(φ) from point C and the parameter angel φ are introduced as shown in Fig. 2 (Norbury, 1973; Choi, 2020).

(30) x(φ)=-q(φ)sinφ,σ(φ)=1+q(φ)cosφfor(x,σonc

Non-dimensionalized cross-sectional area of the core, circulation, and impulse are as follows (Choi, 2020):

(31) A=πα2=1202πq2dφ
(32) Γ=1α2∫∫σdxdσ=1α202π(12q2+13q3cosφ)dφ
(33) P=πα2∫∫σ3dxdσ=πα202π(12q2+q3cosφ+34q4cos2φ+15q5cos3φ)dφ

For a given value of α, in this study, we aim to ensure that the fluid velocity components calculated from Eqs. (28) and (29) satisfy the boundary conditions in Eq. (15) and to find a contour shape whose cross-sectional area satisfies Eq. (31).

Pozrikidis (1986) demonstrated that kinetic energy in Eq. (19) can be calculated using line integral as shown below.

(34) T=πα2∫∫ψσdxdσ=πα2[148α2cσ6dx-12cσ2ψdx+18cσ4uσdσ+18cσ4uxdx]

However, the stream function value satisfies the following relation on contour c (Norbury, 1973; Choi, 2020).

(35) κ=ψ(x,σ)-12Uσ2=const for (x,σ)onc
, where constant κ can be calculated based on the method proposed by Choi (2020) using the deduced contour shape. By substituting Eq. (35) into Eq. (34) and considering the calculation equations of Γ and P, kinetic energy can be expressed as follows:
(36) T=πκΓ+UP+πα2[148α2cσ6dx+18cσ4uσdσ+18cσ4uxdx]

The kinetic energy was computed using a two-dimensional integral in Eq. (34) because the fluid velocity cannot be determined if the CD method for a stream function is used (Choi, 2020). By using the CD method for the fluid velocity proposed in this study, we can substantially reduce the number of computations because contour integral in Eq. (36) is used in computations.

U values deduced from vortical center and impulse center are as follows:

(37) U=dx¯Γdt=1α2Γ[cuxσxdσ-cuσσxdx]
(38) U=dx¯Pdt=πα2P[cuxσ3xdσ-cuσσ3xdx]

## 3. Numerical Analysis Method

Considering the shape symmetry, the shape is discretized using nodal points and segments, as shown in Fig. 3, to determine the contour shape (Choi, 2020).

Based on the definition of point C (Fig. 2) and symmetry, the following relations are established.

(39) φ1=0,φN+1=π,qN+1=q1,x1=xN+1=0,σ1=1+q1,σN+1=1-q1
(40) φ2N-n+2=2π-φn,q2N-n+2=qn,x2N-n+2=-xn=-qnsinφn,σ2N-n+2=σn=1+qncosφnforn=2~N

If N number of q1 to qN, which satisfy Eqs. (15) and (31), can be determined for given α, then the discretized shape of contour can be determined. Furthermore, the forward speed U is an unknown value that should be determined. Therefore, the total number of unknowns is N + 1. Given that Eq. (15) is automatically satisfied if the field point is (x1,σ1) or (xN+1,σN+ 1), the number of conditional equations that can be applied using this equation is N − 1. The following conditional equation determined from the relation between Eqs. (37) and (38) can be used; the total number of conditional equations is N + 1 when Eq. (31) is used, and therefore, N + 1 unknowns can be calculated.

(41) 1Γ[cuxσxdσ-cuσσxdx]=πP[cuxσ3xdσ-cuxσ3xdx]

To determine the solution for N + 1 nonlinear system of equations, an initial shape was assumed and an iterative method was applied (Choi, 2020). In this study, Broyden’s method was applied as the iterative method (Press et al., 1992). Furthermore, angle was divided into equal intervals.

Contour integrals in Eq. (28) and Eq. (29) are summed as segment integrals as follows:

(42) ux(xn,σn)=i=12Nδux,i(xn,σn)
(43) uσ(xn,σn)=i=12Nδuσ,i(xn,σn)
(44) δux,i(xn,σn)=1α2[-ci(xn-x)G(s)dσ+σnciH(s)dx]
(45) δuσ,i(xn,σn)=1α2ciσH(s)dσ
, where (xn,σn) is the field point, and ci denotes the i-th segment. In this study, segments were simplified as linear segments for taking integrals (Shariff et al., 1989; Shariff et al., 2008; Choi, 2020).
(46) x=xi+ξlx,σ+σi+ξlσ,

where lx=xi+1xi, lσ=σi+1σi, 0 ≤ ξ ≤ 1

Eq. (46) can be substituted into Eqs. (44) and (45) as follows:

(47) δux,i(xn,σn)=1α2[-lσ01(xn-xi-ξlx)Gdξ+σnlx01Hdξ]
(48) δuσ,i(xn,σn)=1α2lσ01(σi+ξlσ)Hdξ

When ni and n≠(i+1), general numerical integrals can be considered because integrands of Eqs. (47) and (48) show a non-singular behavior. In this study, numerical integrals were performed using three-point Gauss–Legendre quadrature (Choi, 2020).

However, when n = i or n =(i+1), integrals should be performed by considering logarithmic singularity generated when ξ = 0 or ξ = 1, respectively. When n = i, Eqs. (47) and (48) are expressed as follows.

(49) δux,n(xn,σn)=1α2lx[lσ01ξGdξ+σn01Hdξ]
(50) δuσ,n(xn,σn)=1α2lσ[σn01Hdξ+lσ01ξHdξ]

Based on the methods proposed by Shariff et al. (1989), Shariff et al. (2008), and Choi (2020), H and G were asymptotically expanded at ξ = 0 to take integrals analytically.

(51) (HorG)=ln(8σnlξ)j=0JMajξj+j=0JMbjξj+O(ξJM+l1nξ),wherel=lx2+lσ2

In Eq. (51), the integration values calculated by setting JM = 5 are presented in Appendix. When δux,n(xn +1,σn+1) and δuσ,n(xn+1,σn+1), the signs of lx and lσ are reversed in Eqs. (49)(50) and Eqs. (A1)(A3); when σn+1 is used instead of σn, the signs of the integrals are reversed.

A derivative of a contour curve is required to find cosθ and sinθ, which are the components of a normal vector in the boundary condition of Eq. (15).

(52) cosθ=(dq/dφ)cosφ+x(dq/dφ)2+q2,sinθ=(dq/dφ)sinφ-1+σ(dq/dφ)2+q2

A Fourier cosine series was used in this study to determine the slope of a continuous contour curve from the discretized contour shape.

(53) q(φ)m=0Nqc,mcosmφ,
(54) qc,0=1π0πq(φ)dφ12πi=1N(qi+qi+1)(φi+1-φi)
(55) qc,m=2π0πq(φ)cosmφdφ1πi=1N(qicosmφi+qi+1cosmφi+1)(φi+1-φi)for m1
dq/ is approximated from Eq. (53) as follows:
(56) dqdφ-m=1Nmqc,msinmφ

The solution does not converge or cannot be determined even when numerical differentiation was performed by using the nearby nodal points or derivatives of a cubic spline curve to determine the derivative of the contour curve.

## 4. Analysis Results

The initial guess for the iterative method was set similar to the method in a study by Choi (2020). When α≤0.95, a circle with a radius of α was set as the initial guess for counter; for 0.95 <α≤1.0, a circle with the radius of 0.95 was used as the initial guess for contour. Here, the initial guess of U was set to 0.5. When α >1.0, the solution for α = 1.0, 1.1, 1.2, 1.3, 1.39 was sequentially calculated to set the initial guess value.

The number of segments used for the analysis was set identical to that in the study by Choi (2020), or 2N = 120.

The core shape of several α values is illustrated in Fig. 4. An interpolated value based on the Fourier analysis was used for the shape between nodal points. Specifically, half of the results of this study and those of the study by Choi (2020) were illustrated considering symmetry. The core shape was almost identical to the result in the study by Choi (2020), in which the CD method for a stream function was used. Thus, using the CD method for the fluid velocity is also appropriate for analyzing the vortex ring in a steady state. Furthermore, the excellence of the CD method based on the direct method is exhibited when compared to the results of the study by Norbury (1973) (Choi, 2020).

Shariff et al. (1989) and Shariff et al. (2008) conducted a dynamic analysis by applying the shape results of Norbury (1973) as the initial condition when α = 0.6 to verify the convergence of the dynamic analysis results based on the number of segments. Regarding the analysis results based on 200, 400, 800, and 1,200 segments, the shape of a steady state was maintained even after time had passed when the number of segments was 1,200. Conversely, in this study, an analytical method was proposed to determine the shape of the N-F family vortex ring via iterative methods based on the initially estimated shape, and the results were superior to those reported by Norbury (1973).

If only computation convenience is considered, the method proposed by Choi (2020), which uses a stream function to directly express shapes, is superior than the method proposed in present study, which estimates shapes based on two velocity components and shape slope. However, the CD method for fluid velocity should be introduced to expand the proposed method for an unsteady fluid analysis.

Fig. 5 illustrates the analysis results for forward speed U, κ in Eq. (35), and the core volume VΣ. Furthermore, constant κ was determined at each nodal point by applying the deduced contour to the method proposed by Choi (2020). The constant value must be identical in principle, but a small numerical error is observed at each nodal point. In this study, an arithmetic mean of κ of each nodal point is used. The results of this study correspond to those of the studies by Choi (2020) and Norbury (1973), in which stream functions were used.

The analysis results of circulation, vortical impulse, and kinetic energy are illustrated in Fig. 6. Given that Norbury (1973) and Choi (2020) only used stream functions when calculating kinetic energy, a two-dimensional integral in Eq. (34) should be performed. However, the number of computations can be drastically reduced because kinetic energy is calculated using the line integral of Eq. (36) based on the fluid velocity. The kinetic energy result obtained in this study corresponds to the result obtained by Choi (2020), but the values are slightly greater than the result reported by Norbury (1973). Circulation (Γ), vortical impulse (P), and kinetic energy (T)

## 5. Conclusion

In this study, we examined whether the CD method for fluid velocity, which is used for dynamic analysis, can also be applied for the analysis of a vortex ring in a steady state. When compared to conventional analysis results based on the Stokes’ stream function, the CD method for fluid velocity is applicable for analyzing the N-F family of vortex rings.

The speed of the vortical center and impulse center corresponds to the forward speed of a vortex ring, and the solution can be obtained by additionally applying a conditional equation, which sets the two speeds as identical. A normal vector can be obtained at the discretized nodal point of a shape via Fourier analysis. The normal vector was eventually used to determine the converged iterative calculation results. The degree of numerical integration was heightened by analytically performing the integral of logarithmic singularity.

In future studies, the stability of the N-F family of vortex rings due to a small disturbance can be examined based on the proposed method.

## Notes

Conflict of Interest

## References

Akhmetov DG. 2009. Vortex Rings New York, USA: Springer.
Batchelor GK. 1967. An Introduction to Fluid Dynamics Cambridge, UK: Cambridge University Press.
Choi Y.-R. 2020;Analysis of Steady Vortex Rings Uising Contour Dynamics Method for the Stream Function. Journal of Ocean Engineering and Technology 34(2):89–96. https://doi.org/10.26748/KSOE.2020.010 .
Fraenkel LE. 1970;On Steady Vortex Rings of Small Cross-Section in an Ideal Fluid. Proceedings of the Royal Society A 316(1524):29–62. https://doi.org/10.1098/rspa.1970.0065 .
Fraenkel LE. 1972;Examples of Steady Vortex Rings of Small Cross-Section in an Ideal Fluid. Journal of Fluid Mechanics 51(1):119–135. https://doi.org/10.1017/S0022112072001107 .
Gradshteyn IS, Ryzhik IM. 2000. Table of Integrals, Series, and Products 6th edth ed. San Diego, USA: Academic Press.
Helmholtz H. 1867;On Integrals of the Hydrodynamical Equations, Which Express Vortex-Motion. Philosophical Magazine and Journal of Science (4th series) 33(226):485–512. https://doi.org/10.1080/14786446708639824 .
Hill MJM. 1894;On a Spherical Vortex. Philosophical Transactions of the Royal Society A 185:213–245. https://doi.org/10.1098/rsta.1894.0006 .
Lamb H. 1932. Hydrodynamics 6th edth ed. Cambridge, UK: Cambridge University Press.
Krueger PS, Moslemi AA, Nichols JT, Bartol IK, Stewart WJ. 2008;Vortex Rings in Bio-inspired and Biological Jet Propulsion. Advances in Science and Technology 58:237–246. https://doi.org/10.4028/www.scientific.net/AST.58.237 .
Norbury J. 1973;A Family of Steady Vortex Rings. Journal of Fluid Mechanics 57(3):417–431. https://doi.org/10.1017/S0022112073001266 .
Pozrikidis C. 1986;The Nonlinear Instability of Hill’s Vortex. Journal of Fluid Mechanics 168:337–367. https://doi.org/10.1017/S002211208600040X .
Press WH, Teukolsky SA, Vetterling WT, Flannery BP. 1992. Numerical Recipes in Fortran 77 2nd edth ed. Cambridge, UK: Cambridge University Press.
Pullin DI. 1992;Contour Dynamics Method. Annual Review of Fluid Mechanics 24:89–115. https://doi.org/10.1146/annurev.fl.24.010192.000513 .
Shariff K, Leonard A, Ferziger JH. 1989. Dynamics of a Class of Vortex Rings (NASA Technical Memorandum 102257) Moffett Field, USA: NASA.
Shariff K, Leonard A, Ferziger JH. 2008;A Contour Dynamics Algorithm for Axisymmetric Flow. Journal of Computational Physics 227(10):9044–9062. https://doi.org/10.1016/j.jcp.2007.10.005 .
Smith SGL, Chang C, Chu T, Blyth M, Hattori Y, Salman H. 2018;Generalized Contour Dynmics: A Review. Regular and Chatic Dynamics 23(5):507–518. https://doi.org/10.1134/S1560354718050027 .
Zabusky NJ, Hughes MH, Roberts KV. 1979;Contour Dynamics for the Euler Equations in Two Dimensions. Journal of Computational Physics 30:96–106. https://doi.org/10.1016/0021-9991(79)90089-5 .

## Appendices

### Appendix

The integral values at the logarithmic-singular segment based on the asymptotic expansion of Eq. (51) are expressed as follows.

(A1) 01Hdξ1122880πσn5(225l4lσ-180l4σn-600l2lσ3+864l2lσ2σn-1440l2lσσn2+3840l2σn3+280lσ5-480lσ4σn+960lσ3σn2-2560lσ2σn3+15360lσσn4+61440σn5)ln(8σnl)-1921600πσn5(2025l4lσ-1125l4σn-7650l2lσ3+9504l2lσ2σn-11700l2lσσn2+5270lσ5-8400lσ4σn+15000lσ3σn2-32000lσ2σn3+57600lσσn4+460800σn5)
(A2) 01ξHdξ1430080πσn5(675l4lσ-525l4σn-1800l2lσ3+2520l2lσ2σn-4032l2lσσn2+10080l2σn3+840lσ5-1440lσ4σn+2688lσ3σn2-6720lσ2σn3+35840lσσn4+107520σn5)ln(8σnl)-190316800πσn5(173475l4lσ-95550l4σn-651600l2lσ3+793800l2lσ2σn-959616l2lσσn2+176400l2σn3+446880lσ5-695800lσ4σn+120422lσ3σn2-2469600lσ2σn3+5017600lσσn4+3386800σn5)
(A3) 01ξHdξ1430080πσn5(-405l4lσ+315l4σn+600l2lσ3-840l2lσ2σn+1344l2lσσn2-3360l2σn3+840lσ5-1440lσ4σn+2688lσ3σn2-6720lσ2σn3+35840lσσn4+107520σn5)ln(8σnl)-190316800πσn5(115425l4lσ-66150l4σn-301200l2lσ3+382200l2lσ2σn-508032l2lσσn2+529200l2σn3+94080lσ5+107800lσ4σn-75246lσ3σn2-352800lσ2σn3+10035200lσσn4+11289600σn5)