• Home
• E-submission
• Sitemap
J. Ocean Eng. Technol. Search

CLOSE

 J. Ocean Eng. Technol. > Volume 36(2); 2022 > Article Choi: Analysis of Steady Vortex Rings Using Contour Dynamics Method for Fluid Velocity

### 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^ŽĢ
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)
ŽēŽĢ={╬®ŽāinsideŌłéD0outsideŌłéD
, 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ŽĆŌł½Ōł½Ōł½DŌłćŌĆ▓2uŌåÆ(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^xŌĆēinsideŌĆēŌłéD╬┤-functionŌĆēbehaviorŌĆēacrossŌĆēŌłéD
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ŌåÆ=uŌåÆC+uŌåÆJ
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)
uŌåÆC=╬®Ōł«c[G(sŌĆ▓){-(x-xŌĆ▓)cos╬ĖŌĆ▓+ŽāŌĆ▓sin╬ĖŌĆ▓}-H(sŌĆ▓)Žāsin╬ĖŌĆ▓]e^xdsŌĆ▓
##### (9)
uŌåÆj=╬®Ōł«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)
uŌåÆŌŗģe^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 (x╠ä╬ō ) and impulse center (x╠ä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=╬®╬ōŌł«cxuŌåÆŌŗģe^nds=╬®╬ō(Ōł«cuxŽāxdŽā-Ōł«cuxŽāxdx)
##### (23)
x┬»P=ŽĆŽüf╬®PŌł½Ōł½ŌłæŽā3xdxdŽā
##### (24)
U=dxP┬»dt=ŽĆŽüf╬®PŌł«cŽā3xuŌåÆŌŗģe^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╬▒2Ōł«cŽāŌĆ▓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,ŽāŌĆēonŌĆēc
Non-dimensionalized cross-sectional area of the core, circulation, and impulse are as follows (Choi, 2020):
##### (31)
AŌłæ=ŽĆ╬▒2=12Ōł½02ŽĆq2dŽå
##### (32)
╬ō=1╬▒2Ōł½Ōł½ŌłæŽādxdŽā=1╬▒2Ōł½02ŽĆ(12q2+13q3cosŽå)dŽå
##### (33)
P=ŽĆ╬▒2Ōł½Ōł½ŌłæŽā3dxdŽā=ŽĆ╬▒2Ōł½02ŽĆ(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╬▒2Ōł«cŽā6dx-12Ōł«cŽā2Žłdx+18Ōł«cŽā4uŽādŽā+18Ōł«cŽā4uxdx]
However, the stream function value satisfies the following relation on contour c (Norbury, 1973; Choi, 2020).
##### (35)
╬║=Žł(x,Žā)-12UŽā2=constŌĆēŌĆēforŌĆēŌĆē(x,Žā)ŌĆēonŌĆēc
, 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╬▒2Ōł«cŽā6dx+18Ōł«cŽā4uŽādŽā+18Ōł«cŽā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ŽånŌĆēforŌĆēn=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ŽāŌĆ▓+ŽānŌł½ciH(sŌĆ▓)dxŌĆ▓]
##### (45)
╬┤uŽā,i(xn,Žān)=1╬▒2Ōł½ciŽāŌĆ▓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+1ŌłÆxi, 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╬Š+ŽānlxŌł½01Hd╬Š]
##### (48)
╬┤uŽā,i(xn,Žān)=1╬▒2lŽāŌł½01(Žāi+╬ŠlŽā)Hd╬Š
When nŌēĀi 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╬Š+ŽānŌł½01Hd╬Š]
##### (50)
╬┤uŽā,n(xn,Žān)=1╬▒2lŽā[ŽānŌł½01Hd╬Š+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╬Š),whereŌĆēl=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ŌĆēmŌēź1
dq/dŽå 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).

### 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

##### Fig.┬Ā1.
Core of vortex ring and cylindrical polar coordinates ##### Fig.┬Ā2.
Meridional cross section of the core of vortex ring ##### Fig.┬Ā3.
Discretization of contour (Choi, 2020) ##### Fig.┬Ā4.
Contour shapes for various values of ╬▒ ##### Fig.┬Ā5.
Translation velocity (U), ╬║, and core volume (V╬Ż) ##### Fig.┬Ā6.
Circulation (╬ō), vortical impulse (P), and kinetic energy (T) ### 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 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 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 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 . ### 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)

BROWSE ARTICLES
ARTICLE CATEGORY

Browse all articles >

PUBLICATION ETHICS
FOR CONTRIBUTORS
Editorial Office
President Office BD Rm. 1302, 13 Jungang-daero 180beon-gil, Dong-gu, Busan 48821, Republic of Korea
Tel: +82-51-759-0656    Fax: +82-51-759-0656    E-mail: ksoehj@ksoe.or.kr     