%  LaTeX support: latex@mdpi.com 
%  For support, please attach all files needed for compiling as well as the log file, and specify your operating system, LaTeX version, and LaTeX editor.

%=================================================================
\documentclass[fractalfract,article,submit,pdftex,oneauthor]{Definitions/mdpi} 

%--------------------
% Class Options:
%--------------------
%----------
% journal
%----------
% Choose between the following MDPI journals:
% acoustics, actuators, addictions, admsci, adolescents, aerobiology, aerospace, agriculture, agriengineering, agrochemicals, agronomy, ai, air, algorithms, allergies, alloys, analytica, analytics, anatomia, animals, antibiotics, antibodies, antioxidants, applbiosci, appliedchem, appliedmath, applmech, applmicrobiol, applnano, applsci, aquacj, architecture, arm, arthropoda, arts, asc, asi, astronomy, atmosphere, atoms, audiolres, automation, axioms, bacteria, batteries, bdcc, behavsci, beverages, biochem, bioengineering, biologics, biology, biomass, biomechanics, biomed, biomedicines, biomedinformatics, biomimetics, biomolecules, biophysica, biosensors, biotech, birds, bloods, blsf, brainsci, breath, buildings, businesses, cancers, carbon, cardiogenetics, catalysts, cells, ceramics, challenges, chemengineering, chemistry, chemosensors, chemproc, children, chips, cimb, civileng, cleantechnol, climate, clinpract, clockssleep, cmd, coasts, coatings, colloids, colorants, commodities, compounds, computation, computers, condensedmatter, conservation, constrmater, cosmetics, covid, crops, cryptography, crystals, csmf, ctn, curroncol, cyber, dairy, data, ddc, dentistry, dermato, dermatopathology, designs, devices, diabetology, diagnostics, dietetics, digital, disabilities, diseases, diversity, dna, drones, dynamics, earth, ebj, ecologies, econometrics, economies, education, ejihpe, electricity, electrochem, electronicmat, electronics, encyclopedia, endocrines, energies, eng, engproc, entomology, entropy, environments, environsciproc, epidemiologia, epigenomes, est, fermentation, fibers, fintech, fire, fishes, fluids, foods, forecasting, forensicsci, forests, foundations, fractalfract, fuels, future, futureinternet, futurepharmacol, futurephys, futuretransp, galaxies, games, gases, gastroent, gastrointestdisord, gels, genealogy, genes, geographies, geohazards, geomatics, geosciences, geotechnics, geriatrics, grasses, gucdd, hazardousmatters, healthcare, hearts, hemato, hematolrep, heritage, higheredu, highthroughput, histories, horticulturae, hospitals, humanities, humans, hydrobiology, hydrogen, hydrology, hygiene, idr, ijerph, ijfs, ijgi, ijms, ijns, ijpb, ijtm, ijtpp, ime, immuno, informatics, information, infrastructures, inorganics, insects, instruments, inventions, iot, j, jal, jcdd, jcm, jcp, jcs, jcto, jdb, jeta, jfb, jfmk, jimaging, jintelligence, jlpea, jmmp, jmp, jmse, jne, jnt, jof, joitmc, jor, journalmedia, jox, jpm, jrfm, jsan, jtaer, jvd, jzbg, kidneydial, kinasesphosphatases, knowledge, land, languages, laws, life, liquids, literature, livers, logics, logistics, lubricants, lymphatics, machines, macromol, magnetism, magnetochemistry, make, marinedrugs, materials, materproc, mathematics, mca, measurements, medicina, medicines, medsci, membranes, merits, metabolites, metals, meteorology, methane, metrology, micro, microarrays, microbiolres, micromachines, microorganisms, microplastics, minerals, mining, modelling, molbank, molecules, mps, msf, mti, muscles, nanoenergyadv, nanomanufacturing,\gdef\@continuouspages{yes}} nanomaterials, ncrna, ndt, network, neuroglia, neurolint, neurosci, nitrogen, notspecified, %%nri, nursrep, nutraceuticals, nutrients, obesities, oceans, ohbm, onco, %oncopathology, optics, oral, organics, organoids, osteology, oxygen, parasites, parasitologia, particles, pathogens, pathophysiology, pediatrrep, pharmaceuticals, pharmaceutics, pharmacoepidemiology,\gdef\@ISSN{2813-0618}\gdef\@continuous pharmacy, philosophies, photochem, photonics, phycology, physchem, physics, physiologia, plants, plasma, platforms, pollutants, polymers, polysaccharides, poultry, powders, preprints, proceedings, processes, prosthesis, proteomes, psf, psych, psychiatryint, psychoactives, publications, quantumrep, quaternary, qubs, radiation, reactions, receptors, recycling, regeneration, religions, remotesensing, reports, reprodmed, resources, rheumato, risks, robotics, ruminants, safety, sci, scipharm, sclerosis, seeds, sensors, separations, sexes, signals, sinusitis, skins, smartcities, sna, societies, socsci, software, soilsystems, solar, solids, spectroscj, sports, standards, stats, std, stresses, surfaces, surgeries, suschem, sustainability, symmetry, synbio, systems, targets, taxonomy, technologies, telecom, test, textiles, thalassrep, thermo, tomography, tourismhosp, toxics, toxins, transplantology, transportation, traumacare, traumas, tropicalmed, universe, urbansci, uro, vaccines, vehicles, venereology, vetsci, vibration, virtualworlds, viruses, vision, waste, water, wem, wevj, wind, women, world, youth, zoonoticdis 
% For posting an early version of this manuscript as a preprint, you may use "preprints" as the journal. Changing "submit" to "accept" before posting will remove line numbers.

%---------
% article
%---------
% The default type of manuscript is "article", but can be replaced by: 
% abstract, addendum, article, book, bookreview, briefreport, casereport, comment, commentary, communication, conferenceproceedings, correction, conferencereport, entry, expressionofconcern, extendedabstract, datadescriptor, editorial, essay, erratum, hypothesis, interestingimage, obituary, opinion, projectreport, reply, retraction, review, perspective, protocol, shortnote, studyprotocol, systematicreview, supfile, technicalnote, viewpoint, guidelines, registeredreport, tutorial
% supfile = supplementary materials

%----------
% submit
%----------
% The class option "submit" will be changed to "accept" by the Editorial Office when the paper is accepted. This will only make changes to the frontpage (e.g., the logo of the journal will get visible), the headings, and the copyright information. Also, line numbering will be removed. Journal info and pagination for accepted papers will also be assigned by the Editorial Office.

%------------------
% moreauthors
%------------------
% If there is only one author the class option oneauthor should be used. Otherwise use the class option moreauthors.

%---------
% pdftex
%---------
% The option pdftex is for use with pdfLaTeX. Remove "pdftex" for (1) compiling with LaTeX & dvi2pdf (if eps figures are used) or for (2) compiling with XeLaTeX.

%=================================================================
% MDPI internal commands - do not modify
\firstpage{1} 
\makeatletter 
\setcounter{page}{\@firstpage} 
\makeatother
\pubvolume{1}
\issuenum{1}
\articlenumber{0}
\pubyear{2023}
\copyrightyear{2023}
%\externaleditor{Academic Editor: Firstname Lastname}
\datereceived{ } 
\daterevised{ } % Comment out if no revised date
\dateaccepted{ } 
\datepublished{ } 
%\datecorrected{} % For corrected papers: "Corrected: XXX" date in the original paper.
%\dateretracted{} % For corrected papers: "Retracted: XXX" date in the original paper.
\hreflink{https://doi.org/} % If needed use \linebreak
%\doinum{}
%\pdfoutput=1 % Uncommented for upload to arXiv.org

%=================================================================
% Add packages and commands here. The following packages are loaded in our class file: fontenc, inputenc, calc, indentfirst, fancyhdr, graphicx, epstopdf, lastpage, ifthen, float, amsmath, amssymb, lineno, setspace, enumitem, mathpazo, booktabs, titlesec, etoolbox, tabto, xcolor, colortbl, soul, multirow, microtype, tikz, totcount, changepage, attrib, upgreek, array, tabularx, pbox, ragged2e, tocloft, marginnote, marginfix, enotez, amsthm, natbib, hyperref, cleveref, scrextend, url, geometry, newfloat, caption, draftwatermark, seqsplit
% cleveref: load \crefname definitions after \begin{document}

%=================================================================
% Please use the following mathematics environments: Theorem, Lemma, Corollary, Proposition, Characterization, Property, Problem, Example, ExamplesandDefinitions, Hypothesis, Remark, Definition, Notation, Assumption
%% For proofs, please use the proof environment (the amsthm package is loaded by the MDPI class).
% \usepackage{amsmath} 
%GENERAL
\DeclareMathOperator{\sign}{sign}
\newcommand{\Mathematica}{\textit{Mathematica\textsuperscript{\resizebox{!}{0.8ex}{\textregistered}}}}
\def\8{\infty}
\def\oh{\frac{1}{2}}
\def\ot{\frac{1}{3}}
\def\oq{\frac{1}{4}}
\def\tt{\frac{2}{3}}
\def\ft{\frac{4}{3}}
\def\tq{\frac{3}{4}}
\def\ofi{\frac{1}{5}}
\def\twfi{\frac{2}{5}}
\def\thfi{\frac{3}{5}}
\def\sifi{\frac{6}{5}}

% \def\partial{\partial}
% \def\i{\imath\,}
\newcommand*{\I}{\imath}%
\def\alp{\alpha}
\def\bet{\beta}
\def\gam{\gamma}
\def\lam{\lambda}
\def\eps{\epsilon}
\def\ome{\omega}
\def\del{\delta}
\def\dal{\partial_{\alpha}}
\def\dbe{\partial_{\beta}}
\def\dga{\partial_{\gamma}}
\def\dla{\partial_{\lambda}}
\def\dmu{\partial_{\mu}}
\def\dnu{\partial_{\nu}}
\def\const{\textit{ const }}
\def\undertext#1{\vtop{\hbox{#1}\kern 1pt \hrule}}
\def\ra{\rightarrow}
\def\Ra{\Rightarrow}
\def\lra{\longrightarrow}
\def\MO#1{O\left(#1\right)}
\def\abs#1{\left| #1\right|}
\def\tens#1{\lVert #1\rVert}
\def\VEV#1{\left\langle #1\right\rangle}
\def\VVEV#1{\left\langle\left\langle #1\right\rangle\right\rangle}
\def\tr{\hbox{tr}\,}
\def\diag#1{\hbox{diag}\left(#1\right)}
\def\dd#1{\frac{d}{d#1}}
\def\dbyd#1#2{\frac{d#1}{d#2}}
\def\pp#1{\frac{\partial}{\partial#1}}
\def\pbyp#1#2{\frac{\partial#1}{\partial#2}}
\def\ff#1{\frac{\delta}{\delta#1}}
\def\fbyf#1#2{\frac{\delta#1}{\delta#2}}
\def\pd#1{\partial_{#1}}
\def\br{\\ \nonumber & &}
\def\inv#1{\frac{1}{#1}}
% \def\begin{equation}{\begin{equation}}
% \def\end{equation}{\end{equation}}
\def\bea{\begin{eqnarray} & &}
\def\eea{\end{eqnarray}}
\def\ct#1{\cite{#1}}
\def\rf#1{(\ref{#1})}
\let\oldexp\exp
\renewcommand{\exp}[1]{\oldexp\left(#1\right)}
\def\INT#1#2{\int_{#1}^{#2}}
\def\INT#1#2{\int_{#1}^{#2}}
\def\IINT#1{\int_{#1}\int_{#1}}
\def \Eabg{e_{\alpha\beta\gamma}} 
\def \E#1{e_{#1}}
\def \De#1{\delta_{#1}}
%LOCAL
\def\PBR#1{\left[#1\right]}
\def\JAC#1{\frac{\partial\left(#1\right)}{\partial\left(\rho_1,
\rho_2,\rho_3\right)}}
\def\Mint#1{\int \frac{d^{3}#1}{\left(2\pi\right)^{3}}}
\def\C{\oint_C d }
\def\Y{\Psi_C(\gamma)}
\def\Po#1{\int \frac{d^3 \rho}{4 \pi}
\frac{\rho_{#1}}{ \rho^3} }
\def\et{{\mathcal E}}
\def\CL{Clebsch}
\def\SD{Schwinger-Dyson}
\def\PB{Poisson brackets}
\def\KO{Kolmogorov}
\def\RE{\textbf{Rey}}
\def\NS{Navier-Stokes}
\def\KH{Kelvin-Helmholtz}
\def\BS{Biot-Savart}
\def\LI {Liouville}
\def\BL {Beltrami}
\def \BT{Burgers-Townsend}
\def \BRE{Birkhoff-Rott equation}
\def\SYM {symplectomorphisms}
\def\GBF {$\mbox{GBF}$}
\def\PF {\mbox{pf }}
\def\DS {discontinuity surface}
\def\CVS{\textit{CVS}}
\def\CFE{Clebsch Flow Equation}
\def\vi{viscosity}
\def\tu{turbulent}
\def\ve{velocity}
\def\ci{circulation}
\def\pr{probability}
\def\dis{dissipation}
\def\fl{fluctuation}
\def\fu{function}
\def\dt{\partial_t\,}
\def\val{v_{\alpha}}
\def\vbe{v_{\beta}}
\def\vga{v_{\gamma}}
\def\vla{v_{\lambda}}
\def\ral{r_{\alpha}}
\def\rbe{r_{\beta}}
\def\rga{r_{\gamma}}
\def\rla{r_{\lambda}}
\def\rmu{r_{\mu}}
\def\rnu{r_{\nu}}
\def\oal{\omega_{\alpha}}
\def\obe{\omega_{\beta}}
\def\oga{\omega_{\gamma}}
\def\ola{\omega_{\lambda}}
\def\omu{\omega_{\mu}}
\def\onu{\omega_{\nu}}
\newcommand{\Mod}[1]{\ (\mathrm{mod}\ #1)}

\def\Xint#1{\mathchoice
   {\XXint\displaystyle\textstyle{#1}}%
   {\XXint\textstyle\scriptstyle{#1}}%
   {\XXint\scriptstyle\scriptscriptstyle{#1}}%
   {\XXint\scriptscriptstyle\scriptscriptstyle{#1}}%
   \!\int}
\def\XXint#1#2#3{{\setbox0=\hbox{$#1{#2#3}{\int}$}
     \vcenter{\hbox{$#2#3$}}\kern-.5\wd0}}
\def\ddashint{\Xint=}
\def\dashint{\Xint-}
\def\pint{\Xint/}
% \DeclarePairedDelimiter\ceil{\lceil}{\rceil}
% \DeclarePairedDelimiter\floor{\lfloor}{\rfloor}
\newcommand{\image}[2][\empty]{\begin{center}\begin{adjustbox}{max size={\textwidth}{\textheight}}\includegraphics[#1]{#2}\end{adjustbox}\end{center}}
\newcommand{\bZ}{\mathbb{Z}}
\newcommand{\bS}{\mathbb{S}}
\newcommand{\bR}{\mathbb{R}}
\newcommand{\bD}{\mathbb{D}}
\newcommand{\bT}{\mathbb{T}}
\DeclareMathOperator*{\argmax}{arg\,max}
\DeclareMathOperator*{\argmin}{arg\,min}
\DeclareMathOperator{\sech}{sech}
\DeclareMathOperator{\csch}{csch}
\DeclareMathOperator{\arcsec}{arcsec}
\DeclareMathOperator{\arccot}{arcCot}
\DeclareMathOperator{\arccsc}{arcCsc}
\DeclareMathOperator{\arccosh}{arcCosh}
\DeclareMathOperator{\arcsinh}{arcsinh}
\DeclareMathOperator{\arctanh}{arctanh}
\DeclareMathOperator{\arcsech}{arcsech}
\DeclareMathOperator{\arccsch}{arcCsch}
\DeclareMathOperator{\arccoth}{arcCoth} 
\DeclareMathOperator{\erf}{erf}
\DeclareMathOperator{\Li}{Li} 
\renewcommand{\Re}{\textbf{Re }}
\renewcommand{\Im}{\textbf{Im }}

\newcommand{\pct}[2]{
\begin{figure}
    % \centering
    \includegraphics[width=\textwidth]{#1}
    \caption{#2}
    \label{fig::#1}
\end{figure}
}

%=================================================================
% Full title of the paper (Capitalized)
\Title{TOPOLOGICAL VORTEXES, ASYMPTOTIC FREEDOM AND MULTIFRACTALS}

% MDPI internal command: Title for citation in the left column
\TitleCitation{Title}

% Author Orchid ID: enter ID or remove command
\newcommand{\orcidauthorA}{0000-0003-2987-0897} % Add \orcidA{} behind the author's name
%\newcommand{\orcidauthorB}{0000-0000-0000-000X} % Add \orcidB{} behind the author's name

% Authors, for the paper (add full first names)
\Author{Alexander Migdal $^{1,\dagger,\ddagger}$\orcidA{}}

%\longauthorlist{yes}

% MDPI internal command: Authors, for metadata in PDF
\AuthorNames{Alexander Migdal}

% MDPI internal command: Authors, for citation in the left column
\AuthorCitation{Migdal, A.~A.}
% If this is a Chicago style journal: Lastname, Firstname, Firstname Lastname, and Firstname Lastname.

% Affiliations / Addresses (Add [1] after \address if there is only one affiliation.)
\address{
$^{1}$ \quad Department of Physics, New York University Abu Dhabi,%Department and Organization
            Saadiyat Island, 
            Abu Dhabi,
            PO Box 129188, 
           Abu Dhabi,
            United Arab Emirates; am10485@nyu.edu}

% Contact information of the corresponding author
\corres{Correspondence: sasha.migdal@gmail.com}

% Current address and/or shared authorship
%\firstnote{Current address: Affiliation 3.} 
%\secondnote{These authors contributed equally to this work.}
% The commands \thirdnote{} till \eighthnote{} are available for further notes

%\simplesumm{} % Simple summary

%\conference{} % An extended version of a conference paper

% Abstract (Do not insert blank lines, i.e. \\) 
\abstract{We study the Kelvinons: monopole ring solutions to the Euler equations, regularized as the Burgers vortex in the viscous core.
 There is finite anomalous dissipation in the inviscid limit. However,  in the anomalous Hamiltonian, some terms are growing as logarithms of Reynolds number; these terms come from the core of the Burgers vortex. In our theory, the turbulent multifractal phenomenon is similar to asymptotic freedom in QCD, with these logarithmic terms summed up by an RG equation. The small effective coupling does not imply small velocity; on the contrary, velocity is large compared to its fluctuations, which opens the way for a quantitative theory. 
 In the leading order in the perturbation theory in this effective coupling constant, we compute running multifractal dimensions for high moments of velocity circulation, in good agreement with the data for quantum Turbulence and available data for classical Turbulence. The logarithmic dependence of fractal dimensions on the loop size comes from the running coupling in anomalous dimensions. This slow logarithmic drift of fractal dimensions would be barely observable at Reynolds numbers achievable at modern DNS.}

% Keywords
\keyword{Turbulence, Multifractals, Anomalous dissipation, Fixed point, Velocity circulation, Burgers vortex, asymptotic freedom} 

% The fields PACS, MSC, and JEL may be left empty or commented out if not applicable
%\PACS{J0101}
%\MSC{}
%\JEL{}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Only for the journal Diversity
%\LSID{\url{http://}}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Only for the journal Applied Sciences
%\featuredapplication{Authors are encouraged to provide a concise description of the specific application or a potential application of the work. This section is not mandatory.}
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Only for the journal Data
%\dataset{DOI number or link to the deposited data set if the data set is published separately. If the data set shall be published as a supplement to this paper, this field will be filled by the journal editors. In this case, please submit the data set as a supplement.}
%\datasetlicense{License under which the data set is made available (CC0, CC-BY, CC-BY-SA, CC-BY-NC, etc.)}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Only for the journal Toxins
%\keycontribution{The breakthroughs or highlights of the manuscript. Authors can write one or two sentences to describe the most important part of the paper.}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Only for the journal Encyclopedia
%\encyclopediadef{For entry manuscripts only: please provide a brief overview of the entry title instead of an abstract.}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Only for the journal Advances in Respiratory Medicine
%\addhighlights{yes}
%\renewcommand{\addhighlights}{%

%\noindent This is an obligatory section in "Advances in Respiratory Medicine", whose goal is to increase the discoverability and readability of the article via search engines and other scholars. Highlights should not be a copy of the abstract, but a simple text allowing the reader to quickly and simplified find out what the article is about and what can be cited from it. Each of these parts should be devoted up to 2~bullet points.\vspace{3pt}\\
%\textbf{What are the main findings?}
% \begin{itemize}[labelsep=2.5mm,topsep=-3pt]
% \item First bullet.
% \item Second bullet.
% \end{itemize}\vspace{3pt}
%\textbf{What is the implication of the main finding?}
% \begin{itemize}[labelsep=2.5mm,topsep=-3pt]
% \item First bullet.
% \item Second bullet.
% \end{itemize}
%}

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

\begin{document}
\setlength{\headheight}{20.0pt}
\setcounter{section}{-1} %% Remove this when working on the template.

\section{Introduction}

In 1948 Burgers discovered his vortex solution to the \NS{} equations \cite{BURGERS1948}. This exact stationary solution with cylindrical geometry had an anomalous dissipation (a finite dissipation at vanishing viscosity).

It should have been a breakthrough in the theory of Turbulence. Instead, a phenomenological K41 theory by Kolmogorov and Obukhov \cite{Kolm41} dominated the turbulence studies for the next 80 years, while the Burgers' discovery was almost forgotten.  

While qualitatively describing some important turbulence features, the K41 scaling laws led the microscopic theory to a dead end.

The main question: how is the inviscid \NS{} theory different from the Euler theory was unanswered, nor were the other nascent questions like what is the microscopic mechanism of the spontaneous stochasticity?

In the last few years, there has been some progress in understanding the role of the singular topological solutions of the Euler equations, with Burgers solutions resolving singularities. This progress was inspired by the geometric approach to Turbulence initiated in the 80-ties and 90-ties \cite{M88, AM89, M93, TSVS}, leading to the area law prediction \cite{M93}.

This prediction was verified recently in DNS \cite{S19, S21}, which triggered a flux of new studies in the geometric theory of Turbulence.

This work was recently revised and summarized in a review paper \cite{M23PR}.

Judging by the first responses to this review, this theory is too radical to swallow for the turbulent community. Such a  theory needs to be explained and advanced further.

There are many subtleties to be discussed, some minor corrections to be made, and some questions to be answered. Also, there are more data to compare with the Kelvinon theory and the Loop equation solutions.

\section{Trying to bend the Burgers cylindrical vortex into a torus}

Let us start with the Burgers vortex solution of the \NS{} equation:
\begin{subequations}
\begin{eqnarray}\label{BurgersVortex}
    && \vec v =\left\{-a x- g(r) y, -b y + g(r) x, c z \right\};\\
    && \vec \omega = \left\{0,0,\frac{c  \Gamma_B  e^{-\frac{c r^2}{4 \nu }}}{8 \pi  \nu }\right\};\\
    && g(r) = \frac{\left(1-e^{-\frac{c  r^2}{4 \nu }}\right)\Gamma_B  }{2 \pi  r^2};\\
    && r = \sqrt{x^2 + y^2};\\
    && a = b = c/2 ;
\end{eqnarray}
\end{subequations}

This solution describes an infinite cylinder with vorticity decaying as a Gaussian with the width
\begin{eqnarray}
    w = \sqrt{2 \nu/c};
\end{eqnarray}
This constant $c$ represents the strain $\hat S$ in the direction of the symmetry axis.
\begin{eqnarray}
   && c = \vec t \cdot \hat S \cdot \vec t;\\
   && \hat t = \{0,0,1\}; \\
   \label{strain}
   && \hat S = \tens{ S_{\alpha\beta} } =  \frac{\pd\beta \val + \pd\alpha \vbe}{2} = \diag{-a,-b,c}
\end{eqnarray}

The parameter $\Gamma_B$  represents the circulation in the closed loop surrounding this axis far from its center. By the Stokes theorem, this circulation reduces to an area integral over the surface bounded by this loop, which tends to $\Gamma_B$ in the limit when the radius is much larger than the width $w$ of the vorticity core.

This relation can be directly verified by integrating Burgers velocity around the contour and neglecting exponential terms in $g(r)$.

It is also straightforward to compute anomalous dissipation
\begin{eqnarray}
    \et  =  \nu \int d^3r \vec \omega^2 =  \frac{L c \Gamma_B^2}{8 \pi}
\end{eqnarray}
where $L \ra \infty$ is the length of the cylinder.


Let us try to compactify the Burgers vortex by bending a cylinder into the torus.
Consider a circular vortex line with a radius $R$ much larger than the thickness $w$ of the viscous core.

Locally, at the distances from the core $r \sim w$, the torus is equivalent to the cylinder, up to higher order terms in $w/R$. 

At larger distances, the vorticity is equivalent to the delta function in the cross-section plane, and velocity is purely potential.
It adds up from  a linear term and a singular term
\begin{eqnarray}
&&\vec \omega(\vec r) \ra \Gamma_B \vec t \delta(x) \delta(y);\\
   &&\vec v(\vec r) \ra \left(-a x,-b y, c z\right) + \frac{\Gamma_B}{2 \pi} \frac{ \left(-y,x,0\right)}{x^2 + y^2};
\end{eqnarray}

An important property of the Burgers solution is that the strain tensor \eqref{strain} is constant and has no singularities in the inviscid limit, unlike the rotational part of velocity.
This constant strain represents the local value of the nonsingular Euler strain outside the vortex core, continued inside the tube..

The invariant formula for the inviscid limit of vorticity would be the line integral (with $l$ being the length of the line)
\begin{eqnarray}
    &&\vec \omega(\vec r) = \oint_C d \vec C(l) \Gamma_B(l) \delta^3(\vec r - \vec C(l));\\
    && |\vec C'(l)| =1
\end{eqnarray}
When the point $\vec r$ approaches some point $\vec C(l_0)$ at the loop, the loop integration cancels the delta function, and one recovers the Burgers delta function in the $x y$ plane.
\begin{eqnarray}
    &&\vec \omega(\vec r) \ra \vec C'(l_0) \Gamma_B(l_0) \delta(x) \delta(y);\\
    &&\vec C'(l_0) = (0,0,1);
\end{eqnarray}

In the linear vicinity of a point $\vec C(l_0)$ at this loop, we use the Burgers solution with $z = l-l_0$ to find the derivatives
\begin{subequations}\label{velocityOnLoop}
\begin{eqnarray}
    &&\partial_l \Gamma_B =0;\\
    &&\partial_l\vec v\left(\vec C(l)\right) = \hat S\left(\vec C(l)\right) \cdot \vec C'(l) ;
\end{eqnarray}
\end{subequations}
The velocity field here is taken at the loop $C$, which is a center of the Burgers core.

There is an important boundary condition for the strain tensor at the loop
\begin{eqnarray}\label{EigenvalueCondition}
    && \hat S \left(\vec C(l)\right) \cdot \vec C'(l) = c(l) \vec C'(l);\label{EgenEq}
\end{eqnarray}
In general, the eigenvalue $c(l)$ depends on the point $l$ at the loop, as the strain could be a function of coordinates and the loop changes direction.
This eigenvalue must be positive.

This boundary condition is an analog of the \CVS{} conditions \cite{M23PR} (confined vortex surface) relating the vortex sheet shape and the boundary value of the strain.
In the case of the vortex sheet, this relation imposed restrictions on a surface shape, as the Euler strain did not have any free parameters to adjust. 

In the case of the loop, these restrictions can be treated as extra boundary conditions on the more complex Euler flow outside the tube, as we shall see in the next Sections.

These restrictions are quite strong. In the case of a curved loop we are looking for, the constant strain cannot have the loop tangent vector $\vec C' (l)$ as its eigenvector. 

Unlike the Burgers formulas, these relations \eqref{velocityOnLoop} are parametric invariant and do not depend on the coordinate frame.
This invariance makes them correct generalizations of Burger's solution to an arbitrary smooth loop.

These equations can be readily integrated\footnote{A rookie mistake would be to forget the $l$ dependence of the tangent vector $\vec C'(l)$ and the parameter $c$ and write  a non-periodic formula $\vec v \stackrel{?}{=} c\vec C'(l) l$ for the velocity field.}
\begin{eqnarray}
    &&\Gamma_B(l) = const;\\
    && \vec v\left(\vec r\in C\right)= \int_C^{\vec r} \hat S \left(\vec r'\right) \cdot d \vec r'
\end{eqnarray}

We found a parametric invariant vector integral of the strain along the loop. Another way to derive this formula is as follows:
\begin{eqnarray}
    &&\vec v(\vec r \in C) = \int_C^{\vec r} d \vec r'\cdot \vec \nabla \otimes  \vec v\left(\vec r'\right) =\nonumber\\
    && \int_C^{\vec r} d \vec r'\cdot \hat S(\vec r') + \int_C^{\vec r} d \vec r' \times \vec \omega(\vec r')
\end{eqnarray}
The second term here vanishes because the vorticity at the loop $\vec \omega(\vec r')$ is aligned with its tangent vector $d \vec r'$.

From this relation, by integrating by parts, one can derive the following exact relation for the circulation over the large loop $C$
\begin{eqnarray}
    \Gamma_C = \oint_C d \ral \val = - \oint_C \ral d \val = - \oint_C \ral   S_{\alpha\beta} d\rbe 
\end{eqnarray}
In virtue of the eigenvalue equation, this is also equivalent to
\begin{eqnarray}
    \Gamma_C = -\oint_C d l  \vec C'(l) \cdot \vec C(l) c(l)
\end{eqnarray}

The generalization of Burger's anomalous dissipation is also straightforward \cite{M23PR}:
\begin{eqnarray}\label{AnomalousDissipation}
    \et = \frac{\Gamma_B^2}{8\pi} \oint_C d l \vec C'(l)\cdot \hat S \cdot \vec C'(l) = \frac{\Gamma_B^2}{8\pi} \oint_C d l c(l)
\end{eqnarray}
The problem we are facing  is related to the strain. As we have seen, it cannot be a constant tensor (its highest eigenvector must be equal to the local direction of the loop everywhere).

The generic Euler (singular) velocity field, corresponding to the above singular vorticity line, reads
\begin{eqnarray}\label{SingularLineVelocity}
    \vec v(\vec r) \stackrel{?}{=} \vec \nabla \Phi - \frac{\Gamma_B}{4 \pi} \vec \nabla \times \oint_C\frac{d \vec r'}{\abs{\vec r - \vec r'}}
\end{eqnarray}

The harmonic potential $\Phi(\vec r)$ must be such that the strain at the loop
\begin{eqnarray}
    S_{\alpha\beta}\left(\vec C(l)\right) \stackrel{?}{=} \dal \dbe \Phi\left(\vec C(l)\right) 
\end{eqnarray}
has the local tangent vector $\vec C' (l)$ as its main eigenvector at every point on the loop. 

This main eigenvalue $c$ must be positive.
Two lower eigenvalues $-a,-b$ do not have to be equal, as it was revealed by a further study of the Burgers vortex. As was shown by Moffat et. al. \cite{MKO94}, the vortex solution for a general non-axisymmetric strain tends to the symmetric Burgers solution in the turbulent limit $|\Gamma_B| \gg \nu$. 

We are only interested in that limit, so we can skip the requirement of equal lower eigenvalues.

So, is this it? A weak Euler flow regularized by a Burger vortex core?
Not so fast.

This solution would not be valid for an arbitrary smooth loop $C$, because, in general, the eigenvalue condition will not hold.

We are unaware of any theorems that would prove the existence of the harmonic potential with the prescribed main eigenvector of its Hessian $S_{\alpha\beta}$ on a closed loop $C$ in space.
Apparently, such a harmonic potential does not exist for an arbitrary (or even smooth) loop.



\section{Matching principle and anomalies in the Euler Hamiltonian}

The stationary Euler solution we are looking for must minimize the Hamiltonian
\begin{eqnarray}
    H = \int_{\bR^3} \frac{\vec v^2}{2}
\end{eqnarray}

In our case of the singular velocity field, we have to split this energy integral into two parts: inside and outside a thin tube $\mathcal T$ surrounding the loop $C$. The radius $R$ of a local cross-section of the boundary $\partial \mathcal T$ of this tube must be much larger than Burger's thickness $w$ but much smaller than the local curvature radius $R_C$ of the loop $C$.

Under these conditions, the inside of the tube is described by a cylindrical Burgers vortex, while the outside is some Euler flow in remaining space $\mathcal G = \bR^3 \setminus \mathcal T$. This remaining space has the topology of the full torus, the same as $\mathcal T$.
We understand the space $\bR^3$ is compactified as a sphere $\bS^3$ by including the infinity.

At the surface of the tube and its vicinity, for the whole range $ w \ll R \ll R_C$, there must be a match of the inside Burgers flow up to higher order corrections in $w/R$ with the (yet unspecified) Euler flow up to $R/R_C$ corrections.

This is the matching principle, which we suggested first for the vortex sheets, and then for the vortex lines (see \cite{M23PR} and references to earlier work within).

Thus, the Euler Hamiltonian can be written as the sum of two terms
\begin{eqnarray}
    H = \int_{\mathcal G} \frac{\vec v_E^2}{2}+ \int_{\mathcal T} \frac{\vec v_B^2}{2}
\end{eqnarray}

The last term is an anomaly; it was calculated in \cite{M23PR}. Up to negligible power corrections in $R/R_C, w/R$
\begin{eqnarray}
    &&\int_{\mathcal T} \frac{\vec v_B^2}{2} = \frac{\Gamma_B^2}{8\pi}\oint_C d l \left(\gamma +\log \frac{c(l) R^2}{8 \nu} \right);\\
    && c(l) = \vec C'(l) \cdot \hat S\left(\vec C(l)\right) \cdot \vec C'(l) 
\end{eqnarray}

The term with $\log R$
\begin{eqnarray}
    \frac{\Gamma_B^2}{4\pi} \oint_C d l \log R
\end{eqnarray}
is canceled by a similar term coming from the Euler field. 
This cancellation is a consequence of the matching conditions, as discussed in \cite{M23PR}. Here are these calculations. 
Taking the derivative in $R$, we find the contribution of the surface $\partial \mathcal T$
\begin{eqnarray}
   && \vec v_E(\vec r) \ra \frac{\Gamma_B}{2 \pi} \frac{ \left(-y,x,0\right)}{x^2 + y^2};\\
   && \partial_R \int_{\mathcal G} \frac{\vec v_E^2}{2}  = -\int_{\partial \mathcal T} \frac{\vec v_E^2}{2} \ra -\frac{\Gamma_B^2}{8 \pi^2}\int_C d l  \frac{2 \pi}{R} = -\frac{\Gamma_B^2}{4 \pi } \oint_C d l \frac{1}{R}
\end{eqnarray}
which adds up to zero with the derivative of the above $\log R$ term.
Therefore, the sum of these two terms in the Hamiltonian does not depend on $R$ ( up to neglected power corrections).

As a result of this independence, we can take a limit $R \ra 0$ in the regularized Euler part and finally find the anomalous Hamiltonian
\begin{subequations}\label{AnomalousHamiltonian}
    \begin{eqnarray}
   && H = H_E + \frac{\Gamma_B^2}{8\pi} \oint_C d l \left(\gamma +\log \frac{c(l) |C|^2}{8 \nu} \right);\\
   && H_E = \lim_{R \ra 0}\left(\int_{\mathcal G} \frac{\vec v_E^2}{2}+ \frac{\Gamma_B^2}{8\pi} |C| \log\frac{R^2}{|C|^2}\right);
\end{eqnarray}
\end{subequations}

This Hamiltonian has the common flaw of all Euler flows (except topological ones): it continuously scales down to zero $\vec v_E =0, \Gamma_B =0$ by changing the scale of the velocity field.

Therefore, the zero velocity provides the absolute minimum of the Hamiltonian unless there are some topological restrictions on a flow, preventing continuous scaling to zero within a topological class.



\section{Topological Euler Flow}

The way to build topological solutions to the Euler equations was discovered in the eighties \cite{KM80, L81}.
These are spherical \CL{} or Faddeev variables, as we suggested in \cite{M23PR}.
These variables are elements of 2-sphere geometrically, 3D rigid rotators in the Hamiltonian dynamics, or the $\bS^2$ sigma model in the statistical field theory language.

The vorticity is locally parametrized as
\begin{eqnarray}\label{CLS2}
    &&\vec \omega = \oh Z e_{a b c}  S_a  \vec \nabla S_b \times \vec\nabla S_c;\\
    &&S_1^2 + S_2^2 + S_3^2= 1;
\end{eqnarray}
 where $Z$ is some parameter with the dimension of viscosity, staying finite when $\nu \ra 0$. 

 One could use various coordinates on the sphere. In particular, there are canonical coordinates
 \begin{eqnarray}
  &&\phi_1 =Z S_3;\\
  &&\phi_2 = \arg{ (S_1 + \I S_2)};
 \end{eqnarray}

 With these coordinates, the vorticity becomes parametrized as with the ordinary \CL{} variables
 \begin{eqnarray}
     \vec \omega = \vec \nabla \phi_2 \times \vec \nabla \phi_1
 \end{eqnarray}
 except these two variables $\phi_1, \phi_2$ vary on a rectangle $-Z \le \phi_1 \le Z, -\pi \le \phi_2 \le \pi$ rather than the whole plane $R_2$.

The Euler equations are then equivalent to passive convection of the \CL{} field by the velocity field (modulo gauge transformations, as we argue in \cite{M23PR}):
\begin{eqnarray}\label{CLEq}
     &&\partial_t \phi_a = -\vec v \cdot \vec \nabla \phi_a\\
     &&\vec v =\left(\phi_2 \vec \nabla \phi_1\right)^\perp; 
\end{eqnarray}
Here $V^\perp$ denotes projection to the transverse direction in Fourier space, or:

\begin{equation}
    V^\perp_\alpha(r) = V_\alpha(r) + \dal \dbe \int d^3 r' \frac{V_\beta(r')}{4 \pi |r-r'|}
\end{equation}

As we can see from this representation of velocity, it has a gap 
\begin{eqnarray}
    \Delta \vec v(\mathcal S) = 2 \pi n \vec \nabla \phi_1(\mathcal S)
\end{eqnarray}
at the surface where the phase $\phi_2$ has the gap $\Delta \phi_2 = 2 \pi n$.

This surface is bounded by a singular line $C$ where the angular velocity in the transverse plane diverges as $1/r$, preserving a finite circulation for an infinitesimal dual loop $\delta C$ encircling $C$ (Fig.\ref{fig::DeltaC}).

The velocity gap must vanish at the edge of the surface, which requires the boundary condition
\begin{eqnarray}
    \vec \nabla S_3(C) =0;
\end{eqnarray}

The velocity itself must be directed along the edge, i.e., parallel to the tangent vector $\vec C' (l)$ of the loop
\begin{eqnarray}
    \vec v(\vec C(l)) \parallel \vec C'(l)
\end{eqnarray}

 \pct{DeltaC}{The dual loop (red) encircling the monopole ring (blue). The Burgers vortex resolves the singular vortex line. }

This tube cross-section's boundary $\delta C$ is mapped on the loop $\gamma$ on $\bS^2$, covered $n $ times.
Topologically, a circle is mapped on a circle with homotopy $\pi_1(\mathbb S_1) = \bZ$.

The circulation $\oint_{\delta C} \val d \ral$ is related to the vorticity flux through the tube's cross-section, which equates it to the Burgers parameter $\Gamma_B$.

The circulation around a loop on a 2-sphere equals the area of one of two complementary spherical caps $\Omega_\pm;\; \partial \Omega_\pm = \gamma$, depending upon the orientation of $\gamma$, i.e., the sign of the winding number $n$ (see Fig.\ref{fig::OmegaPM}).

\pct{OmegaPM}{The regions at $\Omega_\pm \in \bS^2$, with opposite orientations of the boundary loop $\partial \Omega_\pm = \gamma$. The areas $|\Omega_\pm|$ add up to $4 \pi$.}

In our coordinates, this loop $\gamma$ is some horizontal circle $S_3 = \cos\lambda =\const$.

\begin{eqnarray}\label{GammaB}
    &&\Gamma_B = \oint_{\delta C} \val d \ral  = \int_{\Omega_\pm} d \phi_2 \wedge d \phi_1=  2 \pi n Z(1 - \sign{n}  \cos\lambda);
\end{eqnarray}

The boundary value of $S_3(C)= \cos(\lambda)$ remains as a constant parameter of our \CL{} field, to be determined later from the minimization of the Hamiltonian.

There is, of course, a possibility to get zero circulation (and, therefore, zero velocity, in case $Z=0$ or $\cos \lambda = \sign n$).

Let us now turn to the circulation $\Gamma_C$ around the original loop. In the same way, as with the Burgers circulation $\Gamma_B$, this circulation can be written as a vorticity flux through some Stokes surface bounded by the loop and passing through the Euler region.

A small part of this surface will pass through the Burgers tube. As the vorticity inside the Burgers tube is directed towards its axis,  we can choose this surface to pass this tube in the local tangent plane to the loop. The normal to this local tangent plane $ \vec n \parallel \vec C' (l) \times \vec C" (l)$ is orthogonal to the direction of the vorticity $\vec \omega \parallel \vec C' (l)$.

Therefore, the flux is determined solely by the vorticity in the Euler region, which means that the singular-line solution \eqref{SingularLineVelocity} with potential flow outside the loop does not provide a finite circulation $\Gamma_C$.

This requirement is the ultimate reason for the topological solution.

There is no contradiction at this level with the \CL{} field, but there is an interesting relation based on the Stokes theorem.

Let us compute the Euler flux through the discontinuity surface on the upper side $\mathcal S^+_C$. 
There is only a tangent discontinuity of velocity at the surface $S_C\setminus C$,  coming from the delta function in vorticity. The boundary values are
 \begin{eqnarray}
     &&\vec \omega \ra 2 \pi n \delta(z) \vec \sigma \times \vec \nabla \phi_1;\\
     \label{velGap}
     && \vec v^+ - \vec v^- =  2 \pi n \vec \nabla \phi_1;\\
     && S_3^+  = S_3^-;\\
     && \phi_2^+ =\phi_2^- - 2 \pi n;\\
     && \omega^+_n =  \omega^-_n = - \vec \sigma \cdot \vec \nabla \phi_2^\pm \times\vec \nabla \phi_1^\pm
 \end{eqnarray}
 Here $\vec \sigma $ is the local normal vector to the surface, and $z$ is the normal coordinate.

The flux through the surface is
\begin{eqnarray}\label{GammaC}
    &&\Gamma_C =  \int_{S^+_C} d \vec \sigma \cdot \vec \omega^+ =\int_{\Omega_\pm} d \phi_1 \wedge d \phi_2 = 2\pi m Z(1 - \sign{m} \cos\lambda);
\end{eqnarray}
where $m \in \bZ$ is a winding number for $\phi_2$ around the loop $C$.

These two winding numbers $n, m$ are consistent with the boundary condition \cite{M23PR} at the surface $\partial \mathcal T$ of the infinitesimal tube $\mathcal T$
\begin{subequations}
\begin{eqnarray}
\label{Phi2Map}
    &&\phi_2\left(\vec C(l) + \vec \xi\right) = m \alpha + n \beta;\\
    && \vec \xi =\eps \left( \vec n \cos\beta + \vec \sigma \sin\beta \right);\\
    &&\vec n = \frac{\vec C''}{|\vec C''|};\\
    &&\vec \sigma = \vec C' \times\vec n;\\
    && \alpha = \frac{2 \pi l}{|C|};\\
    && |C| = \oint_C d l;
\end{eqnarray}
\end{subequations}
Here, $\eps \ra 0$ is the radius of the tube.

Let us briefly discuss the topology of the Kelvinon.
The equation \eqref{Phi2Map} implements the mapping of a torus on a circle with homotopy $ \pi_1(\mathbb T^2) \cong \bZ \times \bZ$, which corresponds to a pair of integer winding numbers $n, m \in \bZ$.

What about mapping the 3D space by the \CL{} field $S_a(\vec r)$?

The Hopf mapping of the compactified 3-space to a 2-sphere,  $\bS^3 \mapsto \bS^2$, was already implemented by a spherical \CL{} field in \cite{KM80, L81}, but our Kelvinon is different.

The Kelvinon implements a mapping of the compactified 3-space $\mathcal G \cong\bS^3$ \textit{without a monopole ring} $\mathcal T \cong \bT^3$ onto the spherical cap $\Omega_\pm\cong \bD^2$ rather than the full sphere $\bS^2$.

Topologically,
\begin{eqnarray}
    \left(\mathcal G \cong  \bS^3 \setminus \bT^3 \cong \bT^3 \right)\mapsto \bD^2;
\end{eqnarray}
As we have seen, this mapping is described by two winding numbers. Depending upon the signs of these winding numbers, the \CL{} field maps physical space on one of the two complementary caps on a sphere separated by a circle $S_3 = \cos\lambda$.

The Appendix presents a family of smooth \CL{} fields with desired properties, including winding numbers $n,m$ and a decrease of vorticity at infinity. We do not compute velocity for these examples, as it would require a solution of a nontrivial Neumann problem on a minimal surface bounded by an arbitrary smooth loop $C$.

The vorticity is shown as a vector flow in 3D in Fig.\ref{fig::OmegaPlot},\ref{fig::CylindricLayerOmega},\ref{fig::CutOutOmega}. 
The \Mathematica code producing these plots is published on Wolfram Cloud \cite{MB18}.

\pct{OmegaPlot}{The stream plot of the 3D vorticity field from the topological example of the Kelvinon field is in Appendix.}
\pct{CylindricLayerOmega}{A toroidal layer far from the core of the stream plot of the 3D vorticity field from the topological example of the Kelvinon field in the Appendix.}
\pct{CutOutOmega}{A toroidal layer close to the core of the stream plot of the 3D vorticity field from the topological example of the Kelvinon field in the Appendix.}


Comparing the expression \eqref{GammaC} to Burger's contribution to the line integral of velocity at the center of the core of the tube, we get a self-consistency relation between parameters of the Euler flow
\begin{equation}\label{SelfConsistencyGammaC}
    \oint_C \ral   S_{\alpha\beta} d \rbe =-2\pi m Z(1 - \sign{m} \cos\lambda);
\end{equation}

This relation being linear, the normalization factor $Z$ can be canceled so that this is a restriction on a Faddeev vector field $S_a(\vec r)$.
The eigenvalue equation \eqref{EgenEq} also restricts the field $S_a(\vec r)$, leaving the normalization factor arbitrary.

As suggested in \cite{M23PR}, this factor will be determined from the energy balance between the incoming energy flow and anomalous dissipation.

There is some advance we have recently made in this part of the Kelvinon theory; we present it in the next Section.

\section{Energy balance revisited}

For every stationary solution of the \NS{} equation, the time derivative of any functional of the velocity field, including the Euler Hamiltonian, must vanish.

Our matching principle suggests using this \NS{} relation to fix the remaining free parameter $Z$ of the Euler flow in the \CL{} variables.

The energy dissipation is localized in the Burgers vortex core. We know it in an inviscid limit as a functional of the local Euler strain at the loop \eqref{AnomalousDissipation}. At fixed \CL{} field $S_a(\vec r)$, this dissipation is proportional to $Z^3$ (two powers of $Z$ from circulation $\Gamma_B$ and one from the strain.)

The matching energy pumping into this region can be written in many forms. In the review paper \cite{M23PR}, we used the energy flux through the boundary of the volume $V$, which flux we then estimated as $\eps V$.

There is a more direct approach, using notorious external random forces $\vec f(\vec r)$, which we take as a random uniform vector inside this volume.

This point needs clarification. Usual forces are also a function of time, correlated by a $\delta(t-t')f(\vec r - \vec r')$ with some slow function of space distance $\vec r - \vec r'$. Then the time averaging is assumed, leading to some expressions for the equal time velocity correlations.

In our approach, we replace time averages with ensemble averages. We have an ensemble of stationary flows, each with its uniform force $\vec f$ drawn from a Gaussian distribution. 
Our force is a space-independent Gaussian random vector with delta variance
\begin{eqnarray}
    \VEV{f_\alpha f_\beta} = \sigma \delta_{\alpha\beta}
\end{eqnarray}

This ensemble averaging is equivalent to time averaging with delta-correlated Gaussian force (such force is a different sample of a Gaussian vector at different times).

The energy pumping created by such a force inside a volume $V$ before Gaussian averaging is simply the work made by this force over the net momentum of the fluid inside this volume
\label{epsEq}
\begin{eqnarray}
    &&\et =  \int_V \delta \vec v \cdot \vec f;
\end{eqnarray}

This perturbed velocity $\delta \vec v$, in turn, depends on the force through the \NS{} equation.
We shall assume this relation to be linear, and later we justify that assumption in the inviscid limit utilizing asymptotic freedom.

In the linear approximation, we have to solve linearized Euler equations for $\delta \vec v, \delta p$
\begin{subequations}\label{LinEq}
    \begin{eqnarray}
    &&\hat G\cdot \delta \vec v  +  \vec \nabla \delta p = \vec f;\\
    && \hat G_{\alpha\beta} = (\dbe \val) \hat I + \delta_{\alpha \beta} \vga \pd \gamma \;\\
    && \vec \nabla^2 \delta p = - \vec \nabla \cdot  \hat G\cdot \delta \vec v ;
\end{eqnarray}
\end{subequations}

Symbolically, we can write the result
\begin{eqnarray}\label{Qmatrix}
    &&\et = \vec f \cdot \hat Q\cdot \vec f;\\
    && \hat Q = \IINT V d^3 \vec r d^3 \vec r'\left(\hat G - \frac{1}{\vec \nabla^2} \vec \nabla \otimes \vec \nabla \cdot \hat G\right)^{-1}(\vec r, \vec r')
\end{eqnarray}
Let us count the factors of $Z$ here, assuming $\vec f$ to be $\MO 1$. As we have shown in \cite{M23PR} (see also the next Section), in the inviscid limit, there is asymptotic freedom:
$Z$ grows as a power of the logarithm of the effective Reynolds number 
\begin{eqnarray}
   && Z \sim \left(\log \RE\right)^\ot \ra \infty;\\
   && \RE = \frac{\VEV{c}_C |C|^2}{ \nu}
\end{eqnarray} 
Then, we have $\hat G \sim Z, \hat Q \sim 1/Z$ in the \eqref{Qmatrix}.


This makes the correction $\delta \vec v \sim 1/Z$ much smaller than $\vec v \sim Z$
\begin{equation}
    \delta \vec v /\vec v \sim Z^{-2} \sim \left(\log \RE\right)^{-\tt}
\end{equation}


Now we have justified linear approximation for the perturbation of the Euler equation by an external pumping force.

In general, asymptotic freedom in Turbulence makes the fluctuations of the velocity field around Kelvinon go to zero as a power of the logarithm of the Reynolds number. 

The same thing happens in QCD with fluctuations of the gluon field around the instanton: these fluctuations logarithmically die out compared to the instanton field.

As the background velocity fluctuation $\vec v_0 =\VEV{\delta \vec v}_V$ is a linear function of $\vec f$, the Gaussian distribution of $\vec f$ is equivalent to the Gaussian distribution of $\vec v_0$, which was considered in \cite{M23PR}. The difference is that we now have a microscopic equation \eqref{LinEq}, which allows us to compute this background velocity field once the base Kelvinon flow is known.

Let us turn back to the energy balance. Naturally, only the total energy of the \NS{} flow is stationary. There are some contributions $\et'_d, \et'_p$ from the remainder of the fluid  both to the energy pumping $\et_p$ and to the energy dissipation $\et_d$ in the energy balance equations
\begin{eqnarray}
     \et'_p +\vec f \cdot \hat Q \cdot \vec f - \sigma \tr \tilde Q=   \et'_d + \frac{\Gamma_B^2}{8 \pi}\oint_C d l c(l);\\
\end{eqnarray}
where we subtracted the mean over the random forces from the $\tilde Q$ term to satisfy the mean energy balance. The outside volume energy flow components $\et'_d, \et'_p$ are treated as constants rather than random numbers (self-averaging of the random forces acting in the remaining infinite volume.)

We have to factor our $Z$ and solve the resulting equation (with variables $\tilde X $ corresponding to  $X$ with $Z=1$)
\begin{subequations}
    \begin{eqnarray}
    &&Z^3 A = B + \frac{1}{Z}\left(\vec f \cdot \tilde Q \cdot \vec f  - \sigma \tr \tilde Q\right);\\
    \label{Aconst}
    && A = \frac{\tilde \Gamma_B^2}{8 \pi}\oint_C d l \tilde c(l);\\
    && B =  \et'_p - \et'_d ;
\end{eqnarray}
\end{subequations}


The unknown parameter $B$ here can be estimated as excessive energy pumped into our volume to be dissipated inside the singular vorticity tube.
This parameter is proportional to the missing volume in each of the energy flows $\et'_d, \et'_p$
\begin{eqnarray}
    B \approx \eps V[C] \propto \eps |C|^3
\end{eqnarray}
where $\eps$ is a \KO{} energy flow per unit volume and $V[C]$ is the volume occupied by our soliton around the singular loop $C$. 

The last estimate $V \propto |C|^3$ implies that the loop has only one scale, which can be taken as its length $C$. 

There are two scales for a large, almost flat loop with small normal deviations: the minimal area $|S[C]|$ and the width $\Delta$ of vorticity field distribution around this surface.
In this case $V[C] \sim |S[C]| |\Delta|$.

We postpone the discussion of large flat loops to the next section.

The solution of this quartic equation at small $A $ equals to
\begin{subequations}
    \begin{eqnarray}
   && Z \ra \left(\frac{B}{A}\right)^\ot  + \vec f \cdot \hat M\cdot \vec f - \sigma \tr \hat M;\\
   && \hat M = \frac{\tilde Q}{3 B} ;\\
   && A \sim \frac{1}{\log \RE} \ra 0
\end{eqnarray}
\end{subequations}


\section{Asymptotic freedom revisited}

Now we can revisit and extend the analysis of \cite{M23PR} of asymptotic freedom (logarithmic decrease of running coupling constant of the Kelvinon theory with the local Reynolds number).

Let us elaborate on the relation \eqref{Aconst}. Using the basic formula \eqref{GammaB}, we find
\begin{subequations}\label{TildeConds}
    \begin{eqnarray}
    &&A = g\oint_C d l \tilde c(l) ;\\
    && g = \frac{ \pi n^2 (1- \sign{n}  \cos\lambda)^2 }{2} ;\\
    \label{S3Cond}
    && S_3(\vec C(l)) = \cos\lambda ;\ \forall l;\\
    \label{streainCond}
    && \tilde S_{\alpha\beta}(\vec C(l)) C'_\beta(l) = \tilde c(l) C'_\alpha(l);\\
    && \oint \ral d\rbe \tilde S_{\alpha\beta} = -2 \pi m ( 1 - \sign{m} \cos\lambda);\\
    && \tilde S_{\alpha\beta}(\vec r) = \frac{\dal \tilde \vbe + \dbe \tilde \val}{2};\\
    && \tilde \val = \left(\phi_2 \dal S_3 \right)_\perp;
\end{eqnarray}
\end{subequations}
The variables with tilde $\tilde S,\tilde v,\dots$ correspond to the normalized flow, with $Z \Ra 1$.

Factoring out the $Z-$ dependence in the anomalous Euler Hamiltonian \eqref{AnomalousHamiltonian}, we find
\begin{subequations}\label{TildeHamiltonian}
\begin{eqnarray}
     && H = Z^2 \left(\tilde H + g |C| \log\frac{Z}{\nu}\right) ;\\
   && \tilde H = g\oint_C d l \left(\gamma +\log \frac{\tilde c(l) |C|^2 }{8 }\right) + \lim_{R \ra 0}\left(\int_{\mathcal P} \frac{\tilde v_E^2}{2}+ g |C| \log\frac{R^2}{|C|^2}\right);\\
   && Z = \left(\frac{\eps V[C]}{\oint_C d l \tilde c(l) }\right)^\ot g^{-\ot}
\end{eqnarray}
\end{subequations}

We note that this factor $g$ depends on only one dynamical variable: the boundary value $S_3(C) = \cos\lambda$.
The \CL{} fields $\phi_2, S_3 $ outside the loop $C$ satisfy the passive advection equations \eqref{CLEq}.
The boundary value of the strain $\tilde S$ satisfies the boundary conditions \eqref{TildeConds}.

The parameter $g$ must minimize the Hamiltonian, with fixed loop $C$ and fixed \KO{} constant $\eps$.

Varying the Hamiltonian \eqref{TildeHamiltonian} with respect to $g$, we find the transcendental equation for $g,Z$
\begin{eqnarray}
&& g = \frac{g_0}{-1 +\log \frac{Z}{\nu }};\\
&& g_0 = \frac{2 \tilde H}{|C|};\\
&& \left(-1 +\log \frac{Z}{\nu }\right)Z^3 = \eps R^4[C];\\
&& R^4[C] =\frac{g_0 V[C] }{\oint_C d l \tilde c(l) }
\end{eqnarray}

This $R[C]$ is a functional of the loop $C$. By dimensional counting, it scales as $|C|$
\begin{eqnarray}
    R[C] =   |C|f\left[ \frac{C}{|C|}\right]
\end{eqnarray}
This estimate assumes that there are no other scales in the solution. As mentioned above, for the smooth, almost flat loop $C$, there may be an effective scale in the normal direction to the minimal surface.

Neglecting such a case (it will be considered in the next Section), we can take $R=R[C]$ as a definition of the loop scale and write down the "renormalization group" equation for effective coupling $g_R =\frac{g}{g_0}$
\begin{eqnarray}
   && \dbyd{g_R}{\log R} = \beta(g_R);\\
   && \beta(g_R) = - \frac{4g_R^2}{3 + g_R};
\end{eqnarray}

This beta function corresponds to the transcendental equation for $g_R$
\begin{eqnarray}\label{eqforgR}
    \frac{3}{g_R} - \log g_R = -3 + 4 \log \left(\frac{\eps^\oq R}{ \nu^\tq}\right)
\end{eqnarray}


\pct{RunningCoupling}{The beta function $\beta(g)$ and the running coupling constant $g_R(\log R)$}

The beta function and the solution for $g_R(\log R)$ are plotted in Fig.\ref{fig::RunningCoupling}. There are no roots of this beta function; it monotonously goes to $-\infty$. 

The roots of the beta function would correspond to the fixed point of RG, leading to scaling laws with anomalous dimensions. \textbf{That} would justify the traditional multifractal scaling laws. However, in our theory, the beta function is calculable in explicit form and does not have a root.

Therefore, we have asymptotic freedom instead of the usual multifractal scaling laws. As discussed below, the difference is hard to observe in DNS and even harder in real experiments on Earth.

Now, with asymptotic freedom, we can justify the linearization of the Euler equation in the presence of fixed external random force: this is the leading term of expansion in our running coupling constant.

In the same leading approximation, the $Z$ factor grows as a power of $\log \RE$
\begin{eqnarray}
     Z \ra \left(\frac{\eps}{3} R^4\log \left(\frac{\eps R^4}{ \nu^3}\right)\right)^\ot 
\end{eqnarray}

\section{Multifractals}

For a purely historic reason, the multifractals are assumed to be pure scaling laws for the moments of the velocity differences, or circulation, in our case, as a function of the scale $R$.

Parisi and Frisch \cite{FP85}  borrowed this idea from Conformal Field Theory, where the moments of similar structure functions are described by scaling laws with anomalous dimensions. 

However, this was only half of the story, and the other half was somehow dropped.

In the conformal field theory of critical phenomena in statistical mechanics, the anomalous dimensions are universal functions of the normal dimension of the fluctuating variable  (power of circulation in our case).

At the same time, in asymptotically free theories, like  QCD, the anomalous dimensions depend on the so-called running coupling constant, which tends to zero inversely proportional to the logarithm of scale.

So, the multifractal as a critical phenomenon has some precedent in statistical and quantum field theory. However, it is a dynamic question of whether the effective coupling is a universal number (a simple root of the beta function at finite coupling constant in the case of conformal theory) or it is running to zero like in asymptotically free theories (the double root of the beta function at zero coupling constant).

We can now answer that question for the Kelvinon theory in favor of asymptotic freedom\cite{M23PR}.

Our asymptotic freedom does not correspond to small velocity. On the contrary, velocity is large while its fluctuations are small; the same happens with the gluon field around instanton in QCD.

The unexpected phenomenon in Turbulence is that the powers of the Reynolds number do not appear in our expansion. In the inviscid limit, the positive powers of viscosity become negligible, and the convergence to this limit slows down to the inverse powers of the logarithm of the viscosity.

This expansion in the zeroth approximation involves the nontrivial Kelvinon field $S_a(\vec r) \in \bS^2$, mapping the physical space (with added infinity and removed loop $C$) onto the spherical cap $\Omega_\pm \in \bS^2$.

There are two possibilities for the orientation $\sign n = \pm 1$, corresponding to an upper or lower cap on a sphere. These two possibilities are equivalent as they correspond to the reflection $S_3 \Ra -S_3$. In the asymptotically free limit, the boundary value $S_3(C) = \cos\lambda$ tends to the North or South pole, depending upon the sign of $n$.
\begin{eqnarray}
   \cos \lambda =  \sign n + \MO{(\log \RE)^{-\ot}}
\end{eqnarray}

This asymptotic formula will be compatible with our self-consistency relation \eqref{SelfConsistencyGammaC} for the circulation $\Gamma_C$ only if $\sign m = - \sign n$.
Then both sides of this equation are $\MO 1$ when the running coupling constant $g_R$ goes to zero.
\begin{eqnarray}
     && \oint_C \ral \tilde S_{\alpha\beta} d \rbe =-2\pi m (1 + \sign{n} \cos\lambda) \ra -4 \pi m ;
\end{eqnarray}
In this case, with the $\vec f$ correction in the linear approximation:
\begin{subequations}
    \begin{eqnarray}
    && \Gamma_C \ra 4 \pi m Z \ra  4 \pi m \eps^\ot R^\ft g_R^{-\ot} + 4 \pi m \left(\vec f \cdot \hat M\cdot \vec f - \sigma \tr \hat M\right)
\end{eqnarray}
\end{subequations}

This expression for the velocity circulation has the same structure as the one in \cite{M23PR}
\begin{subequations}
\begin{eqnarray}
    &&\Gamma_C = \tau + \vec \xi \cdot\hat q \cdot \vec \xi ; \\
    &&\vec \xi  \sim \mathcal N(0,1);\\
    && \tau = 4 \pi m \eps^\ot R^\ft g_R^{-\ot} - 4 \pi m \sigma \tr \hat M;  \\
    && \hat q = 4 \pi m \sigma \hat M;
\end{eqnarray}
\end{subequations}

The circulation shift $\tau$ shows the logarithmic growth with the scale $R$ through the running coupling $g_R^{-\ot} \sim (\log R)^\ot$, but the parabolic term $\vec \xi \cdot\hat q \cdot \xi$ does not depend on the running coupling constant.

According to \cite{M23PR}, this corresponds to the circulation PDF
\begin{eqnarray}
    W(\Gamma) \propto \frac{1}{\sqrt{|\Gamma -\tau|}} \exp{ -\frac{|\Gamma - \tau|}{q_0}}
\end{eqnarray}
where $q_0>0$ is the leading eigenvalue of the $3\times 3$ matrix $\hat q$.

The effective fractal dimension for higher moments, corresponding to the saddle point in the integral for the moments
\begin{eqnarray}
    &&\VEV{\Gamma^p} \propto \int d \Gamma  \Gamma^p W(\Gamma) \ra \textit{saddle point};\\
    && \lambda(p, \log R) \equiv \pbyp{\log \VEV{\Gamma^p}}{\log R} \ra  p \pbyp{\log q_0}{\log R} + \left(1 + \frac{1}{2 p}\right) \pbyp{(\tau/q_0)}{\log R}; \text{ if } p > p_c;
\end{eqnarray}

If there is no internal scale in a Kelvinon, the parameter $R$ is proportional to the size of the loop $|C|$, and the area inside the loop scales as $R^2$.
In that region, we have a K41 formula for $\lambda$, up to higher corrections in asymptotically free coupling constant $g_R$
\begin{eqnarray}
    \lambda(p, \log R) \ra  p \pbyp{\log \tau}{\log R} \ra \frac{4 p}{3}\text{ if } p < p_c;
\end{eqnarray}
As it was shown in \cite{M23PR} from an alternative dynamical theory (loop equations), the asymptotic dependence of the circulation PDF in the true inertial range $|\Gamma| \gg \nu, \eps^\ot|R|^\ft \gg \nu  $ must tend to the scaling law $ \Gamma \sim \sqrt{|S_{min}[C]|}$, where $|S_{min}[C]|$ is the minimal area of the surface bounded by $C$.

In our context, this means that $R^{\ft}$ is to be replaced by  $\sqrt{|S_{min}|}$. 

As for the next correction in $\lambda(p, \log \sqrt{|S_{min}|})$ at large $p$, it depends upon the dimensionless ratio $\tau/q_0$.

Asymptotic freedom tells us this ratio goes to infinity as $g_R^{-\ot}$. Differentiating that and using the RG equation, we find
\begin{eqnarray}
   && \frac{\tau}{q_0} \propto  g_R^{-\ot};\\
   && \pbyp{(\tau/q_0)}{\log R} \ra  g_R^\tt \left(\const{}  + \MO{g_R}\right)
\end{eqnarray}
Thus, our prediction for the asymptotic behavior of the running index $\lambda(p, \log \sqrt{|S_{min}|})$ is
\begin{eqnarray}\label{lambdaPred}
    &&\lambda(p > p_c,\log \sqrt{|S_{min}|}) \ra  p + \const \left(1 + \frac{1}{2 p}\right) g_R^{\tt};\\
\end{eqnarray}
where $g_R$ is related to $\log \sqrt{|S_{min}|}$ by a transcendental equation \eqref{eqforgR}.
\pct{LambdaClassical}{DNS data \cite{S21} for classical circulation fractal dimensions $\lambda(p)$ fitted against our asymptotic formula \eqref{lambdaPred} with $\const{} g_R^{\tt} = 1.359$ (blue) and K41 line $\frac{4 p}{3}$ (green).}
\pct{LambdaQuantum}{DNS data  \cite{QuantumCirculation21} for quantum circulation fractal dimensions $\lambda(p)$ fitted against our asymptotic formula \eqref{lambdaPred} with $\const{} g_R^{\tt}= 1.442$ (blue) and K41 line $\frac{4 p}{3}$ (green).}

We do not have enough data to compare this formula for varying $\log \sqrt{|S_{min}|}$. However, we have the data obtained for the fractal dimension of classical \cite{S19, S21} and quantum \cite{QuantumCirculation21} circulation in turbulent flows.

These authors used a square loop $C$ with a variable side $a$ and fitted their data to constant values of $\lambda(p)$ for each $p$ over the inertial range of $\log a$, neglecting possible systematic deviations from that fit. 

In the future, it would be interesting to fit that data by our formula with the running coupling constant $g_R(\log a)$.

Here are the DNS data from these papers, fitted by our formula with constant $g_R$ Fig.\ref{fig::LambdaClassical}, Fig.\ref{fig::LambdaQuantum}.
The more detailed analysis of underlying data may reveal this slow dependence of $\log a$.

We are not trying to undermine the important numerical work \cite{S19, S21, QuantumCirculation21} and successful phenomenological models \cite{SY21, SY22} explaining these data by approximating the chain of the Hopf equation for velocity moments.  

The phenomenological theory cannot distinguish between the constant fractal dimensions and those running with logarithms of scale; this is a task for a microscopic theory.

The microscopic theory so far says nothing about the transient region $p \sim 5$ where both asymptotic laws break. Only the phenomenological theory \cite{SY21, SY22} accurately describes the fractal dimensions (for velocity differences only !) in the whole domain of small and large $p$.

\section{Conclusions}

Here is what we added to the previous results published in the review paper \cite{M23PR}.
\begin{itemize}
    \item We clarified the topology of the Kelvinon. Its boundary value $\phi_2(\partial \mathcal T)$ at the surface of the infinitesimal tube $\mathcal T$ surrounding the singular line $C$ maps a torus on a circle, which mapping is described by two integer winding numbers related to velocity circulations around two cycles of the torus.
    \item The 3D field $S_a(\vec r)$ maps the compactified 3-space without the infinitesimal tube $\mathcal T$ onto one of the two caps on a sphere $\bS^2$ separated by the circle $\gamma: S_3 = \const{}$.
    \item We modified the energy balance analysis of \cite{M23PR} using conventional random forces and expanding the energy pumping into the Kelvinon in series in the running coupling constant $g_R \sim 1/\log R$. This approach gives us a microscopic definition and corrects the $\log R$ dependence of the phenomenological parameters in the circulation PDF tails \cite{M23PR}.
    \item We studied the vorticity field in the topological family of Kelvinons and presented 3D stream plots of vorticity Fig.\ref{fig::OmegaPlot},\ref{fig::CylindricLayerOmega},\ref{fig::CutOutOmega}. 
    \item We found the self-consistency conditions for the Kelvinon field \eqref{SelfConsistencyGammaC}, \eqref{EigenvalueCondition} from the matching conditions with the Burgers vortex, overlooked in the \cite{M23PR}.
    \item Using these conditions, we removed the ambiguity in relative signs of the winding numbers $n, m$: they must have opposite signs.
    \item We computed and compared predictions of the Kelvinon theory for the fractal dimension $\lambda(p, \log r)$ in the leading perturbation expansion in $1/\log r$ with the DNS data (Fig.\ref{fig::LambdaClassical}, \ref{fig::LambdaQuantum}) with a good fit except for the transient region $4 < p < 7$.
\end{itemize}

In conclusion, we suggested a microscopic quantitative approach to the turbulence problem, and we presented some predictions for the multifractal indexes, modified by powers of the logarithm of the scale.

Here is what still needs to be elaborated and clarified.
\begin{itemize}
    \item The self-consistency conditions \eqref{SelfConsistencyGammaC}, \eqref{EigenvalueCondition} need to be investigated further. Presumably, the boundary values of the \CL{} field on each side of the discontinuity surface provide the set of free parameters needed to satisfy these self-consistency conditions.
    \item The notion of the region occupied by Kelvinon needs to be clarified and defined unambiguously. With correct definition, observable results should not depend upon the shape of the boundary of this region, and its volume should be a well-defined functional of the loop $C$.
    \item Higher correction in perturbation expansion in the running coupling constant $g_R$ need to be computed; fractal dimensions should become universal functions of the logarithm of scale without any phenomenological parameters to fit the DNS data.
\end{itemize}
\newpage
\section*{Acknowledgments}

I benefited from discussions of this theory with my friends and colleagues: Konstantin Khanin, Camillo De Lellis, Nader Masmoudi, Sasha Polyakov, Karim Shariff, Katepalli Sreenivasan, Grigory Volovik, Pavel Wiegmann and Victor Yakhot. 

I am also grateful to Kartik Iyer and Juan Ignacio Polanco for providing me with some data from their DNS of classical and quantum Turbulence.

This research was supported by a Simons Foundation award ID $686282$ at NYU Abu Dhabi.

\section*{Data Availability}
Data sharing does not apply to this article, as no new data were created or analyzed in this study.
% \bibliography{aipnum}

\bibliography{bibliography}
\appendix
\setcounter{section}{-1}
\section{Topological family of the Kelvinon fields}\label{ExampleClebsch}
   

Let us present an explicit example of the \CL{} field with the required topology, which could serve as initial data for the Hamiltonian minimization by relaxation.

We introduce a surface of the minimal area bounded by our loop $C$
\begin{eqnarray}
   S_{min}(C) = \argmin_{S: \partial S = C} \int_S  d S 
\end{eqnarray}
For every point $\vec r \in R_3$, there is the nearest point $\vec r_1$ at the minimal surface $S_{min}(C)$. 
\begin{eqnarray}
   \vec r_1 = \argmin_{ \vec r' \in S_{min}(C)} (\vec r - \vec r')^2;
\end{eqnarray}

   
For this point $\vec r_1 $ at the surface there is also a nearest point $\vec r_0$ at its edge $C$, minimizing the geodesic distance  $d( a,  b)$ along the surface from $\vec r_1$ to the edge
\begin{eqnarray}
   &&s_0 = \argmin_{s} d(\vec r_1, \vec C(s));\\
   && \vec r_0 = \vec C(s_0);
\end{eqnarray}

Let us also introduce the local frame with vectors $\vec t(s), \vec n(s), \vec \sigma(s)$ at the loop:
\begin{subequations}
\begin{eqnarray}
&& \vec C'(s)^2 =1 ;\\
&&\vec t(s) = \vec C'(s);\\
&& \vec n(s) = \frac{ \vec C''(s)}{|\vec C''(s)|}\\
&&\vec \sigma(s) =\vec t(s) \times \vec n(s);
\end{eqnarray}
\end{subequations}\label{frame}

   

Our field is then defined as follows:
\begin{subequations}\label{alphabetarho}
\begin{eqnarray}
&& \alpha = \frac{2 \pi s_0}{ \oint | d \vec C|};\\
&& \beta = \arg \left((\vec r - \vec r_0)\cdot \left(\vec n(s_0) + \imath\vec \sigma(s_0)\right) \right);\\
&& \rho = \sqrt{(\vec r - \vec r_1)^2 + d(\vec r_1, \vec r_0)^2} ;\\
\label{thetaEq}
&& \theta =f(\rho^2);\\
&& f(0) = \lambda;\;f(\infty) = \pi;\; f'(\rho^2) >0;\\
&& \phi =  m \alpha +  n \beta;\\
&&\vec S =  \left(\sin\theta \cos\phi, \sin\theta \sin\phi, \cos\theta\right) ;
\end{eqnarray}
\end{subequations}

In the variational solution for the Kelvinon flow, one may use this Anzatz and optimize the smooth monotonous function $f(\rho^2)$ to reach the minimum of the Euler Hamiltonian at fixed winding numbers $n, m$.

In this example, the \CL{} field maps the physical space $\mathcal G = \bR^3 \setminus\mathcal T$ onto the disk on a 2-sphere $\Omega_+ : S_3 < \cos \lambda$.
With the function $f(\rho^2)$  monotonously \textit{decreasing} from $f(0) =\lambda $ to $f(\infty) = 0$, this Kelvinon would map the physical space to the complementary region $\Omega_-$.

When the point $\vec r $ approaches the nearest point $\vec r_1$ at the minimal surface, the difference $ \vec \eta = \vec r - \vec r_1$ is normal to this surface.

When the point $\vec r_1$ approaches the nearest point $\vec r_0$ at the edge $C$ of the surface, the geodesic becomes a straight line in $R_3$, tangent to the surface and $\rho$ becomes Euclidean distance to the loop
\begin{eqnarray}
  &&d(\vec r_1, \vec r_0) \ra |\vec r_1 - \vec r_0|;\\
   &&\rho^2 \ra  |\vec r - \vec r_1|^2 + |\vec r_1 - \vec r_0|^2 = |\vec r - \vec r_0|^2 
\end{eqnarray}

   
Note that all variables $s_0, \vec r_0, \vec r_1, \alpha, \beta,\rho,\theta, \phi$ depend on $\vec r \in R_3$ through the minimization of the distance to the surface and the loop. By construction, $\rho= |\vec r - \vec r_0|$ away from the surface, when $\vec r_1 = \vec r_0, d(\vec r_1, \vec r_0) =0$.


The Euler angles $\theta, \phi$ for the \CL{} field take the boundary values at the loop :
\begin{eqnarray}
&& \phi(\vec r \ra C) \ra  m \alpha + n \beta;\\
&&\theta(\vec r \ra C) =  \lambda + \MO{(\vec r-C)^2};
\end{eqnarray}
and $\theta(\infty) = 0 \textit{  or }\pi$. 

Assuming the decay $f(\rho^2) \ra f(\infty) + \const{}/\rho^2$, one may estimate the decay rate of $\vec \nabla \cos\theta \sim 1/|\vec r|^3, \vec \nabla \phi\sim 1/|\vec r|$, which corresponds to vorticity decaying as $1/|\vec r|^4$.
This decay is sufficient for the convergence of the enstrophy integral at infinity.

The \BS{} integral for velocity corresponding to such vorticity would decay as $1/\vec r^2$ or faster, sufficient for convergence of the Euler Hamiltonian at infinity.
   
Let us move $\vec r$ along the normal from the surface at $\vec r_1$. Our parametrization of $\theta$ does not change in the first order in normal shift $\vec \eta = \vec r - \vec r_1$, as $\rho^2$ has only quadratic terms in $\vec \eta$.

 We conclude that the normal derivative of the \CL{} field $\phi_1 = Z(1+ \cos\theta)$ vanishes
\begin{eqnarray}
\partial_n \phi_1=0;
\end{eqnarray}

 We requested vanishing normal velocity at the discontinuity surface for this surface to be stationary. 
 
 
 In terms of the \CL{} parametrization, the normal velocity would vanish provided
 \begin{eqnarray}\label{dnS3}
\partial_n \Phi = \left(\vec \nabla \times \Psi\right)_n; \; \vec r \in S_C\setminus C
\end{eqnarray}

   
As for the angular field $\phi$ in \eqref{alphabetarho}, its normal derivative does not vanish in the general case. The angle $\alpha$ does not change when the point $\vec r$ moves in normal direction  $\vec N(\vec r_1)$ from the surface projection $\vec r_1$ by infinitesimal shift $\vec \eta = \eps \vec N(\vec r_1) $. However, another angle $\beta$ changes in the linear order in $\eps$ as
\begin{eqnarray}
   \vec N(\vec r_1) \cdot \left(\vec n(s_0) + \imath\vec \sigma(s_0)\right) \neq 0
\end{eqnarray}

\end{document}
