Reconnection of Infinitely Thin Antiparallel Vortices and Coherent Structures

One of the characteristic features of turbulent flows is the emergence of many vortices which interact, deform, and intersect, generating a chaotic movement. The evolution of a pair of vortices, e.g., condensation trails of a plane, can be considered as a basic element of a turbulent flow. This simple example nevertheless demonstrates very rich behavior which still lacks a complete explanation. In particular, after the reconnection of the vortices some coherent structures with the shape of a horseshoe emerge. They have a high level of complexity generated by the interaction of waves running from the reconnection region. These structures also turn to be very reminiscent to the ones obtained from the localized induction approximation applied to a polygonal vortex. It can be considered as an evidence that a pair of vortices creates a corner singularity during the reconnection. In this work we focus on a study of the reconnection phenomena and the emerged structures. In order to do it we present a new model based on the approximation of an infinitely thin vortex, which allows us to focus on the chaotic movement of the vortex center line. The main advantage of the developed model consists in the ability to go beyond the reconnection time and to see the coherent structures. It is also possible to define the reconnection time by analyzing the fluid impulse.


Introduction
One possible way of transition from laminar to turbulent flow is the interaction between vortices.This is a fascinating process characterized by the reconnection phenomena when the topology of vortices changes producing a cascade of smaller structures which can also reconnect.Simple examples of such processes are the interaction of a pair of vortices like a reconnection of aircraft condensation trials, considered by Crow in [3], or the collision of vortex rings [30,37].In both cases the vortices firstly undergo long wave deformation and then reconnect generating a series of smaller rings or eye-shaped structures as those in the figure 1.This sudden change of the flow structure is quite impressive and still not completely understood.It is also interesting that there is a surprising similarity between the statistical behavior of vortex filaments in turbulent flows for quantum and classical fluids [8,33].Therefore, an explanation of this phenomena may be very useful for the understanding of turbulence.
Vortices are regions in the fluid were a circular movement of particles happens.These regions can be recognized by a high magnitude of the vorticity field which is the curl of the velocity.However, the vorticity field moves with the flow obeying the Navier-Stokes equations, so the vortices are also moving and deforming as if they were objects embedded into the flow.The fact that the vortices are part of a flow and also moved by it leads to that the detection and extraction of vortices from the solutions of the Navier-Stokes equations is quite challenging, and it is even more complicated to follow their evolution [23].
We develop a new model of vortex interaction based on the approximation of infinitely thin vortex embedded into a non-viscous fluid.This model allows us to analyze the behavior of vortices before and after the reconnection moment.After the reconnection time the model predicts the emergence of structures very reminiscent to those obtained from the localized induction approximation (LIA) applied to an isolated eye-shaped vortex.The solution of LIA is related to Riemann's non-differentiable function (RNDF, see formula (43)).The ability of the developed model to pass beyond the reconnection time and see the coherent structures allows to establish a relation between a classical mathematical object such as RNDF and the turbulent flows.
In the simplest cases of the reconnection it is enough to consider a pair of vortex tubes, that is, cylindrical regions where the vorticity lines are parallel, and the vorticity magnitude is decaying far from the cylinder central line.Even in this case we have two different processes: deformation of the vortex core and deformation of the central line.The first one leads to the creation of the helical Kelvin waves seen in many experiments [34], numerical simulations [7], and described in detail in [32].The presence of the Kelvin waves leads to a less clean reconnection process what makes more challenging to understand the phenomena.Their amplitude, however, depends on the size of the vortex core so in this article we focus on infinitely thin vortices in order to avoid these waves.We will refer to this approximation as the vortex filament approximation.In our model the asymmetric helical structures of Kelvin waves before the reconnection disappear, but the symmetric ones persist.It can be considered as an indicator that these waves have different nature: the ones before are due to the deformation of the core, but the ones after are due to the emerged singularity.
Many attempts to describe the vortex reconnection were done in the last 50 years.It started with the pioneer paper of Crow [3] where the evolution of aircraft condensation trails is studied.A pair of initially straight vortices undergoes slow, nearly planar sinusoidal instability.However, the planes containing the deformed vortices are different and incline one to another by approximately π/4.The amplitude of the deformation grows and after some time the vortices reconnect forming a train of vortex rings.The stability analysis is performed in [3] under the assumption that the vortices are infinitely thin, and the viscosity can be neglected.Suppose that X i (s, t) ⊂ R 3 are curves defining central lines of both vortices.We can find the velocity v of the flow in any point x outside the vortices using the Biot-Savart integral: Here t is time, s is the parameter of the curve, i ∈ {1, 2}, Γ i is the vortex strength, the symbol ∧ defines the vector product, the domain is supposed to be R 3 , and the vortices are infinite.The integrals (1) after the introduction of a cutoff can be also applied to a point on the vortex providing the velocity ∂ ∂t X i (s, t).In [3] a linear perturbation analysis is applied to a pair of initially antiparallel vortices finding that the most unstable mode is the long symmetric wave, called Crow wave.There are also short and asymmetric waves seen in experiments and predicted by the model, but their growing rate is much slower, and they probably emerge only under certain atmospheric conditions.The wave length predicted by the Crow model for aircraft condensation trails equals 8.6b where b is the initial distance between vortices (aircraft wingspan).This result slightly exceeds the wavelength obtained in experiments [18] and numerical simulations [12].It happens because the finite size core enhance the growth of shorter waves especially in the case of intense turbulence.
Even though the Crow model predicts instability it is hard to use it for the numerical simulation or any further analysis.A series of simpler models was proposed in several papers by Klein, Majda, and Damodaran [25], [26].Since there is only long wave deformation we can choose an orthonormal basis e 1 , e 2 , e 3 and suppose that the vortices are nearly parallel to e 3 : where δ 1, and X i (s, t) is always orthogonal to e 3 , so it can be considered as a 2D vector.The allowed wavelength is proportional to δ, but they are still long comparing to the initial distance between vortices that is proportional to δ 2 .Under this assumption the Biot-Savart integral can be approximated up to the leading order in δ: where α i is a constant which depends on vortex core, Γ i is the vortex strength, and There are two terms in (3): the first one is the local self-induction that is the velocity of the flow generated by the vortex itself, whereas the second one is the velocity produced by the external flow generated by other vortices.The reconnection after finite time for this model is proven in [39].The equations (3) a much easier to use for numerical simulation.However, due to the singularity in the second term it is impossible to go beyond the reconnection time.Furthermore, it is not clear if the assumption (2) holds true when the amplitude of Crow waves is large and consequently the distance between vortices is small.Another possible approach to study the reconnection phenomena is to consider the Navier-Stokes equations.In this case we do not have any problems related to the singularity and can include all the details such as compressibility, viscosity, and core deformation.The simulation, however, requires more computations and is more difficult.The direct numerical simulation (DNS) of the incompressible Navier-Stokes equations is done in [6].The visualisation of vortices with the λ 2 criterion [13] shows the flattering of the vortex core near the reconnection region and the formation of threads between vortex rings.These threads are stretching and may reconnect again if the Reynolds number is big enough.The presence of viscosity makes it difficult to follow vortex lines due to the dissipation.Therefore, the topology of the vortices is defined by surfaces not by the curves and is much more complicated to analyse.In [23] the vortex lines are extracted from the solution of the Navier-Stokes equations and classificated into ones that reconnect, threads, and the additional vortex rings which emerges in the reconnection zone.It again demonstrates that these phenomena contains a lot of different effects.Another attempt is done in [14] where the reconnection processes is divided in 3 stages: (i) the vortex cores flattering and stretching in the reconnection region; (ii) cutting and reconnection of the inner vortex lines that leads to formation of bridges; (iii) formation of threads from the rest vortex lines where the energy is dissipating through a turbulent cascade.More details are given in the review article [15].The evolution of threads in the last stage is very complicated and chaotic for high Reynolds number, that is for thin vortices or almost non-viscous fluids.Thus, if we consider the Navier-Stokes equation for small Reynolds number then we have to deal with the deformation of the core and the dissipation due to viscosity which does not allow us to extract main coherent structures.On the other hand for high Reynolds number we face the chaotic behavior of the threads after reconnection which is also difficult to filter.We can also note that the coherent structures emerge not only in the reconnection region but also far from it as a result of the interaction of running waves.These structures are clearer, have a distinctive horseshoe shape, and they are the focus of our attention.In particular, we highlight in this paper that they have a behavior very reminiscent to the evolution of the corner vortex under LIA.
We show in figure 1 the main features that emerge due to the reconnection process.We perform the solution of the Navier-Stokes equations using the large eddy simulation (LES) in OpenFOAM software.In figure 1a the symmetric large length Crow waves emerge.Further, in figure 1b we can notice small length asymmetric Kelvin waves in the region close to the reconnection.And finally after the reconnection we can see in figure 1c one horseshoe structure in black rectangle and the bridge going to another symmetric one.In [6,14,15,21] the bridge is called thread whereas the horseshoes are called bridges but in this article we will follow the introduced terminology.In this paper we are not interested in the vortex reconnection cascade [22] but in the evolution of the horseshoe structures and in the complexity of the interaction of the waves that emerge.The extraction of these structures and the definition of the reconnection time is quite challenging due to the finite thickness of the vortices in the Navier-Stokes simulation.
The reconnection of vortices also happens in superfluids [8,5].Even though, there are many differences between classical fluids and superfluids some features of the turbulent regime have similarities, e.g.classical turbulence flows have a filamentary structure [33].Therefore, the results in superfluid reconnection are also important for us.Instead of classical vortices, in superfluids there are quantum vortices that are topological defects where the density tends to zero.Usually the quantum vortices are studied using Gross-Pitaevsky equations [2] or doing vortex filament approximation [19] since the quantum vortices are infinitely thin.There is still a problem on how to model the reconnection and the change of the topology of the vortices.In [19] it is done by a heuristic way of    The configuration with the horseshoe and helical waves is very reminiscent to the evolution of a vortex filament that moves according to LIA and at the initial time is given by two half-lines that meet at point (the corner) with an angle θ [28,11,35].We will call this vortex the corner vortex.Its evolution is depicted in figure 2. The curve is self-similar and has constant curvature c(s, t) = c 0 / √ 2t, and torsion τ (s, t) = s/t, where s is the arclength of the curve, t is time, c 0 and θ satisfy the relation sin θ 2 = e −πc 2 0 /2 (see [28]).It is possible to see that the corner is turned into a circular horseshoe quite similar to the one we can see in figure 1c.The parameter c 0 and the angle θ are estimated from the quantum vortex reconnection experiment in the superfluid helium 4 He studied in [5].The obtained data is consistent with the analytical results from [28], so one would expect the local induction approximation to work reasonably well, at least for quantum vortices.One of the goals of this work is to generalize the model by adding the interaction term and to establish a relation between the self-similar behavior of an infinitely thin corner vortex and the reconnection of vortices.
In this article we present a new model based on the vortex filament approximation but with the interaction term similar to the one in [26].This model is able to go beyond the reconnection moment in a natural way, without any heuristic, and generate coherent structures.It seems to behave as the one showed in the figure 2 that we have just explained.The reconnection time can not be found from the configuration of vortices, so we can not change the topology without introducing an error.However, if we focus on an integral quantity the reconnection time becomes much more clear.We study the fluid impulse that is an integral of a cross product of position and vorticity calculated around the reconnection region.Initially this integral changes monotonously but at some time starts to oscillate.We consider this time as the reconnection time.Furthermore, the behavior of the fluid impulse after reconnection is very reminiscent to the RNDF.A similar effect for a polygonal vortex was discussed in [11], and we can consider this as an evidence, that the vortices at the reconnection seem to create a corner similar to the one described in [28].
The paper has the following structure.In the section 2 the derivation of the equations from the Biot-Savart law are shown.Next, in the section 3 we discuss some properties of the new model and its relation to the previous ones.In the section 4 we describe the numerical method we use to solve the derived b e 1 e 3 e 2 X(s, t) Y(s, t) Figure 3: Initial configuration of vortices equations.Finally, in the section 5 we present results of the numerical simulation, and in the section 6 we compare the behavior after the reconnection with the one of an isolated vortex which has the shape of an eye.It can be used as an approximation of the vortex which emerges after the reconnection (see figure 1c) and also as analog of a curvilinear polygon with only two corners.In the section 7 and we make conclusions and discuss possible directions of further research.The appendix A is devoted to an alternative formulation and possible simplification, the appendix B contains the numerical algorithm, and in the appendix C the evolution of the eye-shaped vortex is described.

Derivation
The velocity of the flux produced by a pair of infinitely thin antiparallel vortices is given by the Biot-Savart law: where Γ is the circulation, X(s) and Y(s) are curves in R 3 defining central lines of both vortices, the symbol ∧ defines the vector product.The vortices are moving by the flow generated by them, therefore X(s, t) and Y(s, t) are also functions of time.For a point x belonging to the vortex X(s, t) the first integral represents the velocity due to the local self-induction v lia , and the second integral is the velocity of the external flow v ext produced by the vortex Y(s, t).Then, we can decompose the velocity of the vortex filament into a sum: Let us call the components of vector X(s, t) = x 1 (s, t) x 2 (s, t) x 3 (s, t) T .
The initial configuration of the vortices is depicted in figure 3. Now we make the following assumptions: 1.The vortices are symmetric respect to the plane x 1 = 0 hence we can reduce our problem to only one unknown curve X(s, t) obtaining the second one by: Y(s, t) = −x 1 (s, t) x 2 (s, t) x 3 (s, t) T . (5) 2. For the velocity of the external flow we will use Rosenhead regularization [20]: where r c is a parameter related to the vortex core, s is the parameter of the curve, and q is the parameter of the curve used inside the integral.This regularization prevents the singularity in the interaction term and can be understood as a viscosity effect during the merging of the cores of the vortices.For the self-induction part we do not use the regularization because there is no core merging.
3. The second vortex can be linearized near any point s that is in the interval q ∈ (s ) for some parameter α, furthermore the tails can be neglected in the Biot-Savart law.The length of this interval is almost proportional to the distance between vortices.When the second vortex is far, a large part of it makes a relevant contribution.On the other hand, when it is close the situation is similar to LIA because the second integral in ( 4) is close to singular, thus we can consider only a small leading piece of it.
4. Either x 1 (s, t) or ∂x1(s,t) ∂s are small so the product x 1 (s, t) ∂x1(s,t) ∂s is neglectable.It means that a piece of vortex can be oriented in e 1 direction only when it is close to the second vortex.

Self-induction
Let us start from the first term in (4) which corresponds to the self-induction.Here we follow the standard derivation of the vortex filament equation for binormal flow [24] using LIA.Fix a time moment and select a point X(s).We choose a perturbation ρξ(s) = ρ (N(s) cos(θ) + B(s) sin(θ)) with small ρ, normal N(s), binormal B(s), and some angle θ.The self-induced velocity of the vortex is found as the limit of Biot-Savart integral when ρ goes to zero.Since the expression in the integral becomes singular we can estimate it using only a neighborhood of the point s with a cut-off L lia where we can expand X(q) in the Taylor series up to the third order of |q − s|: subindeces here designate corresponding derivatives.We have also used that |ξ(s)| = 1 and that it is orthogonal to X s (s).The first term in (7) represents rotation around the vortex central line without its alteration whereas the second one is the movement of the central line that gives the self-induced velocity: In the limit ρ/L lia → 0 we obtain the local induction approximation (or the binormal flow): Assume that the cut-off is inverse to the modulus of tangential vector, that is L lia = Llia /|X s (s)|, and introduce the first parameter of our model: We will see later that the modulus of the tangential vector is growing when the time is close to the reconnection moment.However, due to the regularization in the interaction term we can bound it from above with a power of r c .The parameter ε represents the strength of the interaction between vortices.

External flow
The second integral in (4) after applying symmetry assumption, Rosenhead regularization (6), and fixing the time reads dq.
(10) where component y 1 (s) of vector Y(s) is given by (5).We can apply the assumption 3 obtaining: Neglecting the term x 1 (s)x 1,s (s) due to the assumption 4 we obtain that the external velocity is if we choose α large enough.
Main equations Summing up ( 8) with (11) and rescaling the time with [t] = 8πε/Γ we obtain the main equations describing the evolution of a pair of symmetric vortices due to self-induction and interaction: In (12) we have also used that The equations should be equipped with an initial condition X 0 (s) = X(s, 0) which is usually a small perturbation of a line be 1 + se 3 shifted from the origin by a positive value b in e 1 direction and oriented in e 3 , and the boundary conditions which we will suppose periodic on an interval s ∈ (0, S).
The parameter ε represents the strength of the vortex interaction when compared with the self-induction.The bigger it is, the faster the reconnection happens.A more detailed relation between ε and the velocity of the vortices is shown in section 3.In that section we also show how ε influences on the vortex stretching.The parameter r c is necessary in order to avoid the singularity in the interaction term when x 1 = 0.However, it has a physical meaning of viscosity.According to [24] the radius of the vortex core r c ∼ √ νt, where ν is viscosity and t is time.Since we are mainly interested in non-viscous reconnection the parameter r c should be as small as possible.Even though the presence of r c does not allow to see a sharp corner we still can see the effect and complexity of wave interaction at later times.

Some properties of the derived equations
Relation to previous models The Klein-Majda system of equations [26] for a pair of nearly parallel counter-rotating vortices in the symmetric case can be obtained from the equations (12).Indeed, in the considered case, taking into account that the vortices are nearly parallel to e 3 and including regularization, the Klein-Majda system reads: Parameter σ here depends on the structure of the vortex core, X and Y here are 2-dimensional vectors, the third component x 3 (s, t) = y 3 (s, t) = s is known and ignored in the system.We can reduce the number of equations using symmetry: where t = Γ 4π σt is the rescaled time and ε = 1/σ.Now we will consider the equation (12) supposing that the vortex central line is given in the shape where X (2) • e 3 = 0 and ξa 1.It means that similarly to (14) we have only 2 unknowns, and that the vortices may deviate from the straight line in only a long wave shape, comparing with the distance between vortices.In the Klein-Majda paper [25] these waves are called short-waves.However, if we compare them with the distance between vortices they are long.Plugging (15) into (12) and cancelling the amplitude a we obtain: We suppose that the left-hand side and the first two terms of the right-hand side are of the same order whereas the rest is smaller, thus multiplying by a we get: that is equivalent to (14).
Crow instability Linear stability analysis of ( 12) predicts a long-wave instability described by Crow in [3].Suppose that initially the vortices are parallel to e 3 , as depicted in figure 3, and add a perturbation: where b is a half of the initial distance between vortices, δ 1.Without the perturbation the pair of vortices will move in e 2 direction with velocity that is proportional to ε and almost inverse to the distance between vortices.This result coincide with many previous researches [3], with experiments, and with numerical simulation using the Navier-Stokes equations.Moreover, here we get another physical meaning of the parameter ε: the bigger it is the faster the pair of vortices moves in the e 2 direction.Now let us find the frequencies of perturbations for which this straight line solution is not stable.Plugging ( 16) into (12) and keeping only linear terms respect to δ we obtain We have an eigenvalue problem and the perturbed solution ( 16) is unstable if at least one eigenvalue of ( 18) has positive real part.It happens for the following frequencies ω and wavelengths λ = 2π/ω: Since r c represents the radius of the vortex core and is always smaller than the initial distance between vortices the square root in the expression ( 19) is always real.These waves are long and called Crow waves since they were firstly described in [3].
It is important to note that since we consider the equations ( 12) on the interval s ∈ (0, S) with periodic boundary conditions we have to be sure that the Crow waves (19) fit in this interval, that is S ≥ λ.It is also interesting that the velocity (17) depends on ε/b whereas the wavelength (19) depends on √ ε/b therefore we can not reduce number of parameters and consider only the ratio.Further we will see that ε affects not only on the speed of the reconnection but also on the angle the vortices make at that moment.
The modulus of the tangential vector The vortex filament equations preserve the modulus of the tangential vector T = X s .However, when we have the interaction term as in (12) it is not true anymore.Nevertheless, we can derive a closed expression for the modulus |T|.In order to do it we take a derivative of (12) respect to s and calculate the inner product with T: We can also calculate the inner product of X t with e 1 simplifying the expression of the right-hand side of (20): Combining ( 20) with (21) and integrating respect to time we find the expression for the modulus of the tangential vector: where L 0 (s) is a function which does not depend on time and is given by the initial conditions.The modulus T is growing in the reconnection region (that is x 1 goes to 0) and even tends to a singularity when r c tends to 0. It can be understood as a vortex stretching phenomenon.
Self-similar solution When r c goes to zero the equations ( 12) have selfsimilar solutions.Let us define η = s/ √ t and plug X(s, t) = √ tG(η) into (12) assuming that r c = 0: It is easy to see that after multiplying by √ t we get an ODE for G(η): In order to extract the highest derivative we can calculate the cross product of (23) with G (η): , and it can be expressed by lower derivatives similarly to the previous paragraph.Thus, the final equation is: The equation ( 24) should be equipped with two initial conditions: G(0) and G (0).It is not clear which initial conditions we have to impose for the reconnection problem.
The self-similar solution for the model ( 3) is studied in [39].It appears that if a singularity is introduced at the beginning it will persist for the infinite time.Thus, the self-similar reconnection in the model (3) will never have a clear horseshoe in difference with the corner vortex studied in [28].One of the reasons can be that in the model (3) the LIA term has a linear approximation.In the equation (24) we include this term in the complete nonlinear form.
Behavior close to the reconnection point Using formula (22) we can bound T 1 / T before the reconnection moment from bellow.Assume that in the interval s ∈ [s 0 , s 1 ] the component x 1 (s) is growing monotonically (so T 1 (s) ≥ 0) for value m = x 1 (s 0 ) to M = x 1 (s 1 ).This assumption is correct before reconnection, but the numerical simulation shows that it does not hold true after it since the helical waves emerge.Subject to the proposed assumption we can write the following estimation: sup s∈(s0,s1) where L 0 C([s1,s0]) = sup s∈(s0,s1) |L 0 (s)|, and 2 F 1 (a, b, c, d) is the hypergeometric function, the modulus of the tangential vector is given by (22).Assume now for simplicity that r c = 0, the reconnection happens at s 0 = 0, and designate s 1 = s.Then, the estimate (25) reads sup q∈(0,s) If initially the vortices were oriented into x 3 direction and separated by value 2b the norm L 0 C([0,s]) = b ε .The furthest point between vortices corresponds to s = λ/2, where λ is the wavelength of Crow waves given by (19).The value of x 1 (s) in this point is not smaller than b, so we can use it to estimate the ratio This bound is not optimal and is very far from it.Nevertheless, we can expect that the first component of the tangent vector will grow when we increase ε tending it to 1.The numerical experiments in section 5 show that it becomes almost parallel to e 1 , so the shape of the vortex after reconnection is very close to a horseshoe.
We can also use the formula (26) to estimate the maximal possible value of x 1 .Indeed, the left hand side can not be bigger than 1, so assuming that L 0 C([0,s]) = b ε and the maximum for x 1 (s) is achieved at s = λ/2 we obtain; .

Numerical method
In this section we describe the numerical method we use to solve the system (12).The main problem for numerical stability of the method is related to the interaction term which grows when the vortices are close to each other.We consider a simpler case of Klein-Majda equations ( 14) to derive possible restrictions for the numerical method.They come from the relation between the time step τ , the spatial discretization step h, and the regularization parameter r c .Even though we use a more advanced Runge-Kutta-Felhberg method and the equations ( 12) have higher nonlinearity the derived restrictions hold true in a qualitative way.
Necessary stability conditions We derive the necessary stability conditions for a simpler case of Klein-Majda equations (14).Consider the following semiimplicit numerical scheme: n and y n are approximation of first and second components respectively of the solution X(s n , t k ), h and τ are discretization steps for the parameter s and time respectively.Assume now that there is a high-frequency but small amplitude numerical error δ (k) cos(ωn) in the second component and let us analyse how it will grow on the next time step.Plugging the perturbed values x x cos(ωn) and y y cos(ωn) into (28), (29) we get for the following expression for the linear approximation of the error: where The necessary stability condition requires the eigenvalues of the error transformation matrix be not bigger than 1 by modulus.From (30), (31) we obtain the equation of the eigenvalues ν: The product of the eigenvalues is always 1, therefore the necessary stability condition is satisfied if and only if the roots are complex, that is: implying This condition can be resolved providing a constraint for τ : for any ω and x (k) n .The condition has the strongest form when ζ achieves it maximal value (ζ = 4, see (32)) and x (k) n = 0, so the condition for τ reads: Note that for smaller r c we have to use smaller time step.Furthermore, the stability of the scheme is lost when r c tends to 0 what corresponds to the emergence of a singularity in the interaction term.The formula (36) implies that for say a two times smaller regularization parameter r c we have to use a two times smaller time step τ .However, this relation does not hold true for a higher order scheme.Indeed, in that case we will have a higher derivative of the interaction term respect to x.That is the terms h 3 /r 3 c and further will be presented in the constraint for τ .The stability can be obtained by choosing h proportional to r c so all terms in the Taylor expansion of the interaction term will be bounded.However, this choice leads to a very fast growth of computations making it very hard to perform the simulation for small r c .

Description of the numerical scheme
The main challenge in the numerical solution of the equations ( 12) is that at the reconnection moment the behavior of the interaction term is close to singular.This time period should be passed with very small time step which is not needed when the vortices are far from each other.Therefore, we use an adaptive time step tecnique: an embedded 5th Runge-Kutta method in time with 8th order finite difference discretization in the filament parameter s.The 8th order scheme gives the best results of those we have tried.On one hand, the spectral method that has a higher order requires a higher order time scheme.On the other hand, a lower order spatial discretization does not provide sufficient accuracy.We have also studied the possibility of use of implicit methods, such as [36] but these methods suffer the same requirement of the small time step at the reconnection moment.Besides, we are interested in the multifractal behavior of the trajectories of the vortex filament points therefore we need data with very high discretization in time thus such advantage of implicit method as large time step can not be used.
In the implementation of the method we follow the book of J. Butcher [16].The idea of the embedded Runge-Kutta method consists in realization of two Runge-Kutta schemes on the same points one of order p and another one of order p + 1.The difference between the outputs of these methods on each step is used for the error estimation which should have the decay τ 5 .If the error is bigger then we decrease the time step until the accuracy test is not passed.The explicit k-steps Runge-Kutta scheme for an ODE x (t) = f (x, t) at step n is given by where τ is the time step, α i , β il , and c i are the coefficients of the scheme.Usually the coefficients are given in the Butcher table: α β c T .In order to add the accuracy test we have to add another vector of coefficients ĉ for the embedded method.We are using Runge-Kutta-Felhberg method with the coefficients obtained in [4].The adaptive time step allows to decrease the time step when it is necessary.In our case when the reconnection happens the interaction term is very close to singular and therefore, a much smaller time step comparing with the rest of the simulation has to be used.
To make the solution more stable we use the idea of [10] and resolve the equations for X and T at the same time adding also the arclength correction according to (22).The new equations read Expression (22) for the modulus of the tangential vector allows us to avoid the calculation of the derivative: Here L 0 (s) and its derivative are given as initial condition.Following [10] we also perform the correction of the tangential vector modulus after each interaction using formula (22).The method is explained in detail in the appendix B.
Even though there are two connected unknowns x 1 and T 1 , the solution of the system (39)-( 40) provides the correct result T 1 = x 1,s .One can wonder if it is possible to reduce the number of unknowns and what is the minimal number of independent functions which describe the evolution of the vortex reconnection.This is studied in the appendix A using the Frenet frame.It turns out that the reconnection of vortices can be described in terms of just two functions: x 1 (s, t) and its derivative with respect to time x 1,t (s, t).However, the equations in this case contain derivatives respect to s up to the 4th order.Thus, it is much more complicated for the numerical solution.

Results
We use the method described in the section 4 to solve the equations (12).There are a few things we are mainly interested in: (i) the emergence of Crow waves and their length; (ii) the influence of r c to the solution; (iii) the formation of the horseshoe structure and the direction of the tangential vector in that region.

Crow waves
In the first test we start from a random perturbation and check the formation of Crow waves.We start from a small perturbation of a straight vortex separated by b = 0.11 and consider the evolution following the equations (12) selecting ε = 0.05 and r c = 0.025.The results are depicted in figure 4. We can see that at time t = 2 we have almost sinusoidal waves.The wavelength is around π that is very close to one predicted by formula (19) for the given values of ε, b, and r c .
Influence of the regularization parameter r c The behavior of the vortices far from the reconnection is not dependent of r c .However, at times when the vortices are close to each other the regularization parameter starts to play a crucial role.Furthermore, the evolution of vortices after reconnection also depends on the parameter r c .We may expect that when the regularization parameter tends to zero the shape of vortices will be less smooth.Therefore, the behavior of vortices after reconnection will resemble the behavior of a corner vortex.These expectations are confirmed by the numerical solution, see figures 5,6,7.
We have consider a pair of symmetric vortices with initial conditions where b = √ ε/2, δ = b/20, in the interval (0, 2π) discretized with 6000 nodes.The boundary conditions are periodic.The parameter b is selected using formula (19) in such a way that there is exactly one Crow wave in the considered interval.If we decrease b the reconnection may happen in multiple points thus  The vortices start to touch each other at time t = 1.01, figure 5.The influence of r c can be notice only close to the reconnection region and the shape of vortices is sharper for small r c .It is not completely clear what we can call "the reconnection moment": the first touch or the moment when the horseshoe emerges.Both these moments are dependent on r c .However, the second one has a more complicated dependence since the smallest size of the horseshoe is dictated by r c : the smaller r c , the smaller the horseshoe structure will be.This effect is demonstrated in figure 6 where the configuration of vortices at time t = 1.025 is depicted.We can see that even though for all r c we have a contact the horseshoe appears only for small values of r c .Furthermore, only for small values of the regularization parameter we can see the helical waves at this time moment.The fact that we can not define the reconnection moment does not allow us to change the topology of vortices that leads to the formation of the bridge between the horseshoes.The bridge is growing and represents a source of numerical difficulties and possible instabilities at later times.The bigger bridge at time t = 1.05 can be seen in figure 7.There we can also see the horseshoe and the helical waves for large values of r c .It is interesting to note that the difference between solutions for r c = 1.25 • 10 −3 and r c = 3.125 • 10 −4 is almost neglectable.We can expect that there is a convergence when r c tends to zero.
The analysis of the influence of r c to the solution shows an important phenomenon: the reconnection moment can not be determined so the "real corner" between vortices can not be seen.We can see only the consequences, such as the horseshoe starting from some minimal size related to r c .This effect also appears in the experiments [5].Since the regularization parameter have the physical meaning of the radius of vortex core we can expect that for thick vortices we will never see the singularity and the shape after the reconnection will be closer to a vortex ring without any helical waves, since the minimal size of the horseshoe is big and close to the length of Crow waves.When the vortices get thinner the shape after reconnection becomes more complicated.However,    we still do not see the singularity.One of the challenges related to this phenomenon is the incapacity to perform the reconnection, that is to say the change topology, because we do not know when we have to do it.On the one hand if we reconnect the vortices when they touch each other, the corner they create is not the one that generates the horseshoe structure.On the other hand, if we wait until the horseshoe structure emerges we always find an artifact on its tip related to the bridge.
Influece of the parameter ε The results obtained in section 3 predicts that the horseshoe will be closer to circular when the value of ε increases (formula ( 27)).We use the same initial condition (42) as before but now the initial distance b = √ ε/2, and the perturbation δ = b/20 both depend on the parameter ε.It is necessary for two things: firstly, we have to use such initial distance, so there is at least one Crow wave in the interval, secondly the change of perturbation amplitude allows to achieve the reconnection almost at the same time.The regularization parameter is set to r c = 0.05 and the interval (0, 2π) is discretized with 6000 nodes.
The configuration of the vortices at different time moments are depicted in figures 8,9, and 10.In the first figure the reconnection has not happened yet.However, we can see that the vortices with bigger values of ε move faster in x 2 direction and also the amplitude of Crow wave is bigger due to the bigger initial distance b.In the figure 9 the moment when the horseshoe structure emerges is depicted, but there are still no helical waves.For the smaller value of ε the horseshoe structure is not planar and has a sharper tip.The configuration with helical waves and the horseshoe structure is shown in figure 10.For the large value of ε the horseshoe is almost planar and the vortices look very similar to the shape we can see in evolution of a corner vortex in figure 2.
Another important question one may ask is the behavior of the vortices around the reconnection.In particular, is there a cusp?Or can we expect a smooth horseshoe?In order to do it we can analyze the components of the tangential vector T. The bigger value of the first component T 1 corresponds to a smoother horseshoe.Moreover, when T is parallel to e 1 the cusp disappears.A very rough analytical result on the ratio T 1 /|T| is presented in the section 3. Using numerical simulations we can see that the value of T 1 /|T| may be very close to 1 for some ε.In figure 11a it is possible to see that before the reconnection the projection of normalized vector T to the direction e 1 is small   because it is mostly oriented in the e 3 direction.However, in figure 11b when the reconnection already happens we can see that there is a bridge, where T 1 is close to 0, and a jump almost of size 1 at tip of the horseshoe.Moreover, the bigger is ε the bigger is the maximum of T 1 /|T|.This coincides with the analytical result from the section 3.In figure 11c we can see that the bridge is growing and the horseshoe persists.

Comparison with the eye-shaped vortex
We can compare the evolution of the reconnected vortices with an isolated vortex which has an eye-shape and deforms obeying the local induction approximation.These not planar vortex can be considered as an approximation of the shape after reconnection (see figure 1c) or as a "two-corner curvilinear polygon".We will see that the evolution of the eye-shaped vortex has similarities with both the reconnection process and the evolution of polygonal vortex.It means that even though we can not see the corner at the reconnection time, the further evolution of the vortices has a similar structure to the corner vortex (figure 2) under LIA.
It is known [11,27] that a polygonal vortex with M corners has a periodic behavior with period T = 2π/M 2 whereas the trajectory of the corner tends to a modification of Riemann's non-differential function (RNDF): when M tends to infinity, i 2 = −1.RNDF is a periodic multifractal that has a peculiar behavior at points corresponding to rational multiples of the period [40].It can be seen in figure 12.We pick a point t * = π/4 and multiply it by different rational numbers.It is easy to see that most of them corresponds to local minima and maxima of the absolute value of (43) and to corners of the graph in the complex plane.However, points that correspond to the rational multiples with odd denominator fall into a cusp singularity of the absolute value of RNDF seen as a spirale on the complex plane.It was shown in [11] that the polygonal vortices have a similar behavior.Furthermore, at each rational multiple p/q of the period the shape of the corresponding vortex is also polygonal though not necessary planar.Moreover, the angle between two adjacent sides tends to π when q tends to infinity.This dependence on q is more visible in the H δ process studied in [31].
The similar effect can be seen for the eye-shaped vortex which can be considered as a curvilinear polygon with two corners.The eye-shaped vortex is more similar to the configuration we have after the reconnection than a polygonal one therefore the comparison with it is interesting to us.In appendix C it is shown that the evolution of this vortex is quasi-periodic.Moreover, the Fourier analysis of the trajectory that starts in a corner shows that the dominating frequencies are still the squares like in the polygonal case.
The fluid impulse of the reconnected vortices We define the fluid impulse around the corner as where the corner is located at s = 0 and the interval l is 20% of the whole perimeter of the vortex.The fluid impulse is important for us because in the case of the reconnection we can not extract the singularity point hence we can not find its trajectory, whereas the formula (44) is always applicable.It depends however on the position of the origin X 0 which in the case of the eyeshaped vortex we define as X 0 (t) = X(0, t).We will be mainly interested in the oscillations and multi-fractal behavior for which the choice of origin does not have any influence.Despite this similarity to the RNDF, the real structure of the fluid impulse of the eye-shaped vortex is much more complicated, as can be seen in figure 13.The curve is still planar since F 1 (t) = 0 for all time, the rational multiples of a local maxima t * = 0.10848 also corresponds to the local maxima, minima, and singular points but the scaling of self-similar structures are deformed and there is no real period.The dependence on the size of the maxima with respect to the size of the denominator is also presented.
In the case of the reconnection we can not specify the reconnection point and follow its trajectory.However, we can calculate the fluid impulse (44) around the reconnection region.The results for different r c and ε = 0.05 are depicted in figure 14.Analysing the fluid impulse we can detect a sudden change from monotone to oscillating behavior at time t ≈ 1.01.We call this moment the reconnection time.Note that from the configuration of the vortices we could not to define this time, thus analysis of integral quantities such as the fluid impulse is beneficial for understanding of the reconnection phenomena.Moreover, one can note that the behavior after the reconnection time is quite reminiscent to the one of the eye-shaped vortex though the period is different.The smaller is the regularization parameter r c , the more details we can see in the fluid impulse.
The extraction of squares however is not possible for the reconnection problem due to the noise generated by the bridge.Analysis of distribution of minima, maxima, and singular points using wavelet transform have also faced problems related to the noise.Maxima for the considered types of signals corresponds to singular points that can be studied using the multifractal analysis [1,17].The main idea of this approach is to construct the singularity spectrum D(h), that is the function for which each Hölder exponent h yields the Hausdorff dimension of the set of points where the function has this exponent.For RNDF the singularity spectrum is known: D(h) = 4h−2 for h ∈ [0.5, 0.75], D(h) = 0 for h = 1.5 and D(h) = −∞ otherwise [29].We have tried to use the p-leaders method [9] to estimate the spectrum of RNDF and the fluid impulse.Unfortunately these signals are quite delicate.Therefore, even in the case of RNDF the approximation of D(h) is not very accurate (especially for exponents corresponding to D(h) = 0).For the fluid impulse of the polygonal or reconnected vortices the situation is much worse even though there is a tendency that h corresponding to maximal D(h) is decreasing getting closer to 0.75.
Since the multifractal analysis failed, in this paper we perform only a qualitative analysis of similarity between the fluid impulse of the eye-shaped vortex Figure 13: Fluid impulse of an eye-shaped vortex for b = 0.4, θ = π/6 and its modulus.Some maxima a located in times 3/11t * , 1/3t * , 3/7t * , 3/5t * , t * , and the bigger is the denominator, the smaller is the peak at that (similar to what happens in [31]).and the one of the reconnection problem.We choose a point t * related to a local maxima of the reconnection fluid impulse and check if the rational multiples of this point also fall in maxima, minima, or singular points.The results are depicted in figure 15 for ε = 0.03, b = 0.22, r c = 6.25 • 10 −4 .We can see that the rational points mostly fall into local minima and maxima at least for a short time after reconnection.However later we can see that the self-similar structure of the fluid impulse vector is completely lost.In order to obtain a cleaner structure we have to focus on a small time interval after the reconnection and choose a small regularization parameter r c .This is quite challenging, due to the stability condition (36) of the method.Thus, a new approach is required to find more similarities between the reconnection process and the eye-shaped vortex.
The vortex separation rate The scaling law which determines how the distance between vortices changes during the reconnection was studied in multiple works [2,5,14].It is not completely clear if there are different laws before and after the reconnection or if this law varies for quantum and classical fluids.Nevertheless, there are multiple evidences that the separation rate δ(t) of the vortices after the reconnection should be of the scale √ t − t rec where t rec is the reconnection time.This rate can be observed in experiments [5] and also coincide with the analytical result for the corner vortex [28].
In figure 16 the separation rates before and after the reconnection are depicted for ε = 0.05 different values of r c .The black dashed lines corresponds to the scale |t − t rec |.We can see that before and after the reconnection the separation rate is very close to the square root law.The agreement is better for small values of r c that correspond to a case of smaller viscosity.We can also   16b).This result can be considered as another evidence that the behavior of the vortices after the reconnection resembles the evolution of the corner vortex even though we can not see the corner at the reconnection time due to the presence of the regularization parameter r c .

Conclusions
Even though it is not entirely clear what happens in the reconnection time the further evolution of vortices contains coherent structures reminiscent to the ones generated by a polygonal vortex [11,27].In order to extract and analyze these structures a new model (12) describing the interaction of a pair of antiparallel vortices is developed.The introduction of the regularization parameter r c allows us to go beyond the reconnection time.Moreover, the model provides a closed expression (22) for the length of the tangent vector of the vortex central line.This length is proportional to the circulation and, according to our result, it is increasing, when the distance between vortices tends to 0. It can be considered as the vortex stretching phenomenon.The model also predicts the Crow waves firstly described in [3] and the formation of coherent structures.These structures are a pair of horseshoes in the spirit of [19,28] connected by a bridge which is artificial but can not be removed due to restrictions of the model.Nevertheless, we have shown analytically in section 3, that the non smoothness of the horseshoe due to the bridge can be bounded.Finally, the model predicts the square root timescale for the separation rate after the reconnection.
We performed a numerical simulation for different values of the vortex interaction parameter ε and the regularization parameter r c .There is a difference with the evolution of the corner vortex filament shown in figure 2 and the reconnection.In the first case the horseshoe emerges immediately and can be infinitely small, whereas for our model the smallest possible size of the horseshoe is finite and dictated by the parameter r c .This effect makes impossible to determine the reconnection time, so we can speak only about the reconnection interval which is getting smaller when r c tends to zero.However, the reconnection time can be defined better if an integral quantity instead of the configuration of the vortices is considered.We focus on the fluid impulse that is an integral of the cross product between the position and the vorticity calculated in the reconnection region.It is possible to see a sudden change of the behavior of this quantity from monotone to oscillatory happening at the reconnection time (see figure 14).Note also that the smaller r c the more complicated the behavior after the reconnection is.The oscillations look quite reminiscent to the ones of the fluid impulse of the polygonal vortex [11] that it tends to Riemann's non-differential function (RNDF, formula (43)) when the number of polygon sides tends to infinity.It can be considered as an evidence that the antiparallel vortices indeed form a corner at the reconnection.However, a further research is necessary for estimation of the noise produced by the bridge, studying the influence of the interaction, and calculation of the corner angle.
A possible way to find more similarities between the fluid impulse of reconnecting vortices and RNDF is to improve the numerical method.In particular the condition (36) is very restrictive because we have to reduce both space and time discretization to be able to solve problems with small r c .Probably this restriction can be surpassed by choosing right orders in space and time discretizations.Another way is to apply filtering to the already obtained fluid impulse or try more advanced methods of the analysis.We have tried to apply methods of multifractal analysis such as the p-leader method [9].This approach is based on studying the distribution of the singular points of the function and the calculation of the singularity spectrum which can be considered as a fingerprint of the signal.Similar singularity spectrum means similar multifractal properties of the signals.Since, the RNDF is very delicate for such methods, and the fluid impulse of the reconnecting vortices have a lot of noise, our analysis could not arrive to any conclusion.
Advances in the study of the vortex reconnection may have a huge impact to the understanding of turbulence and its structure.
where L = |X s |.On the other hand we can write that X s = LT and take derivative with respect to t: Since for X(s, t) the order of differentiation does not matter we can find the equation for L t : and also expressions for coefficients λ 1 and λ 2 : We can continue this process and obtain the equations for κ 1,t and κ 2,t : What do we need to close the system?In the expression for u and v we use x 1 , n 1 , and n 2 whose equations can be obtained from (54),(52), and (53) respectively: Here we got rid of λ 1 and λ 2 using (60),(61).Expressions for ω and λ 3 are still missing.There are two ways to find the first quantity.First one is to use equations ( 47) and (48): This is a functional equation, and it does not include λ 3 .However, in order the denominator to be nonzero we have to require n 2 1 + n 2 2 = 0 that is T = e 1 .Another way is to take a derivative of (47) respect to t and a derivative of (52) respect to q making them equal, that gives a differential equation for ω t : In order to find λ 3 we have to make assumptions about our frame.Suppose that n 2 = 0 and use the equation (66): Then the direct expressions for N 1 and N 2 are: To make the first vector not zero the tangential vector should not be oriented in the e 1 direction.Numerical experiments show that if r c > 0 this is true before the reconnection and a long time after it.The second vector corresponds to the direction of the interaction term in the system (12).This choice of frame vectors has another advantage that we can express κ 1 using L, x 1 and their derivatives.Indeed, calculating inner product of (46) with N 1 and taking into account that T q • T = 0 we obtain: Using this we can reduce the system of equations to only 3 unknowns: We will show later that the equation ( 72) can be resolved analytically.Besides, similarly to κ 1 we can find κ 2 , and see that it is proportional to the velocity x 1,t .Thus, the behavior of the system of antiparallel vortices is governed by 2 quantities: the distance between vortices x 1 and the velocity of their approximation represented by κ 2 .
Even though the choice of the interaction frame allows us to reduce the number of unknowns it is not very useful in numerical simulation since the 4th derivative of x 1 respect to s is required: equation (73) includes the second derivative of v, that according to (56) depends on κ 1 , that is proportional to the second derivative of x 1 due to (70).In practice, it is better to use the initial formulation (12).

21:
if error < threshold then T ← T + ∆T

C The eye-shaped vortex
The initial configuration of the eye-shaped vortex is given by The evolution of the eye-shape vortex (74) with θ = π/6 and b = 0.4 is shown in the figure 17.It is possible to see that the movement is quasi-periodic with period T = 3.55 since at that time we see that the vortex has again the eyeshape with the same orientation (blue line) but a slightly different parameters than at the initial time (red dashed line).We can also see that at a half-period time the vortex also has an eye-shape but is rotated (green line) similarly to the polygonal vortex [11].At a rational fraction p/q of the quasi-period T we can also see a non planar curve with q or 2q corners for even and odd values of q respectively.This behavior also coincides with the one of a polygon.
In figure 18a the trajectory of the corner X(0, t) and the fluid impulse around the corner are depicted.Analysing the Fourier coefficients (figures 18b and 18c) we can see that for both the trajectory and the fluid impulse the frequencies corresponding to squares of integers are dominating similarly to what happens in the case of regular polygons [11].It makes the behavior similar to the RNDF (43).
(a) Crow waves (time t = 5) (b) Kelvin waves (time t = 6) (c) Horseshoes (one in the black rectangle and the symmetric one) and bridge between them (time t = 7)

Figure 2 :
Figure 2: Self-similar evolution of a corner vortex from different points of view

Figure 4 :
Figure 4: Formation of Crow waves for ε = 0.05 and r c = 0.0025.500 nodes are used in the discretization.In the front view picture we show only right vortex, in the top view the vortices are shifted by the mean value in x 2 direction.

Figure 12 :
Figure 12: Riemann's non-differential function (43) on the complex plane and its absolute value

Figure 14 :
Figure 14: Fluid impulse (44) for the reconnection problem with different values of r c , ε = 0.05, and b = √ ε.The blue line corresponds to normalized and shifted fluid impulse of the eye-shaped vortex with θ = π/6 and b = 0.4.At time t ≈ 1.01 the behavior suddenly changes from monotone growth to oscillation.We can use it as the definition of the reconnection time.

Figure 16 :
Figure 16: The minimal distance between vortices before and after the reconnection for ε = 0.05 and different values of r c

Figure 18 :
Figure 18: Trajectory, fluid impulse, and their Fourier coefficients for the eyeshaped vortex with θ = π/6 and b = 0.4, the lines connect frequencies corresponding to squares of integers.