Renormalization and energy conservation for axisymmetric fluid flows

We study vanishing viscosity solutions to the axisymmetric Euler equations without swirl with (relative) vorticity in Lp\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$L^p$$\end{document} with p>1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p>1$$\end{document}. We show that these solutions satisfy the corresponding vorticity equations in the sense of renormalized solutions. Moreover, we show that the kinetic energy is preserved provided that p>3/2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p>3/2$$\end{document} and the vorticity is nonnegative and has finite second moments.


Introduction
For axisymmetric incompressible flows without swirl, the (originally three-dimensional) Navier-Stokes and Euler equations can be reduced to two-dimensional Communicated  Institut für Analysis und Numerik, Westfälische Wilhelms-Universität Münster, Münster, Germany mathematical models which are obtained by assuming a cylindrical symmetry for both the physical space variables and the velocity components. Despite this simplification, such flows are still able to describe interesting physical phenomena like the motion and interaction of toroidal vortex rings. On the mathematical level, even though two-dimensional, the (vaguely defined) degree of difficulty of analyzing solution properties lies somewhere between that of the two-dimensional planar equations and the full three-dimensional model. Indeed, as we shall see later on, axisymmetric flows 1 do still feature vortex stretching and some of the standard global estimates have an unambiguous three-dimensional character. On the other hand, many of the features of the Biot-Savart kernel are typically two-dimensional even though some helpful symmetry properties are lost.
In the present work, our aim is to study renormalization and energy conservation of solutions to the Euler equations that are obtained as vanishing viscosity solutions from the axisymmetric Navier-Stokes equations. Here, renormalization is to be understood in the sense of DiPerna and Lions [24], that is, a solution is called renormalized if the chain rule of differentiation applies in a suitable way. We are particularly interested into solutions whose vorticity is merely L p integrable in a sense that will be made precise later.
The analogous (though in some parts technically much simpler) studies for the two-dimensional planar equations have been conducted quite recently: As long as the vorticity is L p integrable with exponent p ≥ 2, DiPerna and Lions's theory for transport equations (combined with Calderón-Zygmund theory) ensures that the vorticity is a renormalized solution of the corresponding vorticity equation [39]. This fact is true regardless of the construction of the solution. If p ∈ (1, 2), renormalization properties are proved in [18] for vanishing viscosity solutions. The argument in this work relies on a duality argument and exploits the DiPerna-Lions theory. This theory, however, does not apply to the p = 1 case, in which the associated velocity gradient is a singular integral of an L 1 function. Instead, a stability-based theory for continuity equations proposed in [42,43] can be suitably generalized in order to handle this situation and to extend the results from [18,39] to the limiting case p = 1; see [17].
Conservation of kinetic energy for vanishing viscosity solutions with L p vorticity, p > 1, is established in [14] for the planar two-dimensional setting (on the torus). The corresponding three-dimensional problem gained much attention in recent years, particularly in connection with Onsager's conjecture [40], which states that the threshold Hölder regularity for the validity of energy conservation is the exponent 1/3. Energy conservation for larger Hölder exponents was proved in [16], see also [26] for partial results and [13] for improvements. In particular, in the last paper the authors show conservation of energy for velocities in the Besov space B 1/3 3,s , for 1 ≤ s < ∞, which contains W 1 3 ,3 for all s ≥ 3. Note that by fractional Sobolev inequalities in R 3 the W 1/3+,3 regularity holds for any vorticity in L p with p > 9/5. More recently, the sharpness of the Hölder exponent was proved in [32], building up on the theory developed in [7][8][9]20,21].
Before discussing our precise findings and the relevant earlier results for the axisymmetric equations, we shall introduce the mathematical model. The Euler equations for an ideal fluid in R 3 are given by the system where u = u(t, x) ∈ R 3 is the fluid velocity and p = p(t, x) ∈ R is the pressure. In this formulation, the (constant) fluid density is set to 1. Whenever the fluid has locally finite kinetic energy, which will be the case in the regularity framework considered in this paper, the Euler equations can be interpreted in the sense of distributions.
Definition 1 Let T > 0 and u 0 ∈ L 2 loc (R 3 ) 3 be given. A vector field u ∈ L 2 loc ((0, T )× R 3 ) 3 is called a distributional solution to the Euler equations (1), (2) We restrict ourself to the case of axisymmetric solutions without swirl. That is, if (r , θ, z) are the cylindrical coordinates of a point x ∈ R 3 , i.e., x = (r cos θ, r sin θ, z) T , we shall assume that u = u(t, r , z), and u = u r e r + u z e z , where e r and e z are the unit vectors in radial and vertical directions, which form together with the angular unit vector e θ a basis of R 3 , We remark that u θ = u · e θ is the swirl direction of the flow, that we assume to vanish identically. Under these hypotheses on the velocity field, the vorticity vector is unidirectional, ∇ x × u = (∂ z u r − ∂ r u z )e θ , and we write ω = ∂ z u r − ∂ r u z . A direct computation reveals that this quantity, that we will call vorticity from here on, satisfies the continuity equation on the half-space H = (r , z) ∈ R 2 : r > 0 . We remark that ω is thus a conserved quantity, because the no-penetration boundary condition u r = 0 on ∂H comes along with the symmetry assumptions. However, opposed to the situation for the twodimensional planar Euler equations, the vorticity is not transported by the flow, as the divergence-free condition (2) becomes in cylindrical coordinates. Indeed, the continuity equation can be rewritten as a damped transport equation, where the damping term on the right-hand side describes the phenomenon of vortex stretching, 1 r u r ωe θ = (∇ x × u) · ∇ x u. What is transported instead is the relative vorticity ξ = ω/r , We remark that the flow is entirely determined by the (relative) vorticity, as the associated velocity field can be reconstructed with the help of the Biot-Savart law in R 3 , A transformation into cylindrical coordinates and an analysis of the axisymmetric Biot-Savart law can be found, for instance, in [29]. Thanks to this relation, we may thus study (5), (6) instead of (1), (2). Working with the vorticity formulation has certain advantages: At least on a formal level, it is readily seen that the vorticity equation (5) preserves any L p norm, if ξ 0 is the initial relative vorticity. 2 This observation is crucial, for instance, in order to prove uniqueness in the case of bounded vorticity fields [19]. The drawback of working with (5) is that there is no direct way of giving a meaning to the transport term in low integrability settings (opposed to the momentum equation (1)). For instance, it is not obvious to us, how to extend common symmetrization techniques that allow for an alternative formulation of the transport nonlinearity in the planar two-dimensional setting, see, e.g., [6,22,48]. Whenever the product uξ is locally integrable, we can interpret the transport equation (5) in the sense of distributions.
Definition 2 Let T > 0 and p, q ∈ (1, ∞) be given with 1 is called a distributional solution to the transport equation (5) with initial datum ξ 0 if ξ is axisymmetric and Notice that our formulation relies on the identity u r ∂ r ξ = u 1 ∂ 1 ξ +u 2 ∂ 2 ξ that allow us to switch between Cartesian and cylindrical coordinates. Moroever the definition provides a distributional formulation of the continuity equation (3) in which ω is replaced by r ξ .
Simple scaling arguments show that the local integrability of the product uξ can be expected to hold true only if p ≥ 4/3. For this insight, it is crucial to observe that the Sobolev inequality is valid as in the planar two-dimensional setting, cf. [29,Proposition 2.3]. For vorticity fields with smaller integrability exponents, we propose the notion of renormalized solutions.
is called a renormalized solution to the transport equation (5) with initial datum ξ 0 if ξ is axisymmetric, ξ(t, r , z) and and any bounded β ∈ C 1 (R) vanishing near zero.
We remark that the notion of renormalized solutions implies the conservation of the L p integral of vorticity in the sense of (7) via a standard approximation argument. Moreover, it is shown in [3,24] that renormalized solutions are transported by the Lagrangian flow of the vector field u as in the smooth situation. We will further comment on this in Sect. 2 below. The relation between Lagrangian transport and the partial differential equations (3) and (5) was thoroughly reviewed in [4].
In the present paper, we study solutions to the vorticity equation (5) in the case where the initial (relative) vorticity can be unbounded, more precisely, for some p ∈ (1, ∞). We are thus outside of the class of functions in which uniqueness is known to hold [2,19]. On the positive side, existence of distributional solutions to the Euler equations (1), (2) was proved in [33] for initial vorticities satisfying (9) and under the additional assumption that the initial kinetic energy is finite, u 0 ∈ L 2 (R 3 ) 3 . (Notice that local L 2 bounds on the initial velocity can be deduced from the integrability assumptions on the vorticity via Sobolev embeddings, cf. (8).) For larger integrability exponents and (near) vortex sheet initial data, (crucial insights on) existence results were previously obtained in [11,12,23,34,35,41,44,47]. To the best of our knowledge, renormalized solutions (Definition 3) have not been considered in the context of the axisymmetric Euler equations. We are particularly interested into solutions that are obtained as the vanishing viscosity limit from the Navier-Stokes equations, which are, in fact, physically meaningful approximations to the Euler equation. Hence, for any viscosity constant ν > 0, we consider solutions (u ν , p ν ) to the Navier-Stokes equations We furthermore impose fixed initial conditions, u ν (0) = u 0 and shall assume that u ν is axisymmetric, that is, u ν = u ν (t, r , z) and u ν = (u ν ) r e r + (u ν ) z e z . Instead of working with the momentum equation (10), will mostly study its vorticity formulation, which is a viscous version of (3) (or (5)), see (22) (or (23)) below. It was shown in [29] that under the assumption (9) on the initial data, which implies that Starting from this solution to the Navier-Stokes equations, our first result addresses compactness and convergence to the Euler equations.
Theorem 1 (Compactness and convergence to Euler) Let u ν be the unique solution to the Navier-Stokes equations (10), (11) with initial datum u 0 ∈ L 2 loc (R 3 ) such that the associated relative vorticity ξ 0 belongs to L 1 ∩ L p (R 3 ) for some p > 1. and Moreover, u is a distributional solution to the Euler equations (1), (2) and ω = r ξ is the corresponding vorticity that is (in a distributional sense) related to u by the Biot-Savart law (6).
The vanishing viscosity limit was studied for finite energy solutions with mollified initial datum satisfying the bound (9) in [33]. The novelty in the above result is the kinetic energy may be unbounded. For earlier and related convergence results for non-classical solutions, we refer to [1,31,35,47,49] and references therein. Our next statement concerns the renormalization property of the relative vorticity.
Theorem 2 (Renormalization) Let u and ξ be the velocity field and relative vorticity, respectively, from Theorem 1. Then ξ is a renormalized solution to the transport equation (5) with velocity u. In particular, it holds that and ξ is transported by the regular Lagrangian flow of u in R 3 .
To the best of our knowledge, in this result, renormalized solutions to the axisymmetric Euler equations are considered for the first time. We recall from the above discussion that for p ∈ (1, 4/3), the interpretation of the transport equation (5) as a distributional solution does not apply anymore as the transport nonlinearity is no longer integrable. In particular, while for p ≥ 4/3 our result implies that distributional and renormalized solutions coincide, in the low integrability range, we show the existence of renormalized solutions. We also recall that for p ≥ 2, the result in Theorem 2 is already covered in DiPerna and Lions's original paper [24]. In Sect. 2, we recall the theory from [24] and explain what we mean by ξ being transported by a flow. For a precise definition of regular Lagrangian flows, we refer to [3,4].
Our final result addresses the conservation of the kinetic energy.
Theorem 3 Let p ≥ 3 2 . Suppose that the fluid has finite kinetic energy, u 0 ∈ L 2 (R 3 ) 3 , and that ω 0 is nonnegative and has finite impulse, Then the kinetic energy is preserved, In order to show conservation of energy, the growth of vorticity at infinity has to be suitably controlled. Here, we choose a growth condition that is natural as it can be interpreted as the control of the fluid impulse. Notice that the latter is conserved by the evolution, cf. Lemma 8. This is in principle not required by our method of proving Theorem 3, and any estimate of the form r 2 ω(t) L 1 (H) r 2 ω 0 L 1 (H) would be sufficient. It is, however, not clear to us whether such an estimate holds true under our integrability assumptions apart from the special case considered in Lemma 8, that is, for nonnegative (or nonpositive) vorticity fields. Also, if higher order moments could be controlled, our method shows that the value of p could be lowered (at least up to p > 6 5 ). See, for instance, [12] for similar results in the setting with p > 3 (and general solutions). We leave this issue for future research and consider the simplest case here.
From the result in Theorem 3, it follows that we are outside of the range in which Kolmogorov's celebrated K41 theory of three-dimensional turbulence applies, since, similar to the case of planar two-dimensional turbulence, there cannot be anomalous diffusion.
From here on, we will simplify the notation by writing The advective derivatives f · ∇ and f · ∇ x are to be interpreted correspondingly.
The remainder of the article is organized as follows: In Sect. 2 we recall the parts of the DiPerna-Lions theory for transport equations and explain how the results apply to the setting under consideration. In Sect. 3 we provide estimates for the velocity field that are essentially based on the Biot-Savart law. Section 4 contains global estimates for the axisymmetric Navier-Stokes equations, while the proof of Theorems 1, 2 and 3 are given in Sects. 5, 6 and 7, respectively. This work, finally, contains an appendix in which a helpful interpolation estimate is provided.

Renormalized solutions for linear transport equations
In this section, we shall briefly recall DiPerna and Lions's theory for linear transport equations [24] in the general setting of transport equations in R 3 , thus neglecting the assumption of axisymmetry for a moment. We are particularly interested into wellposedness and renormalization properties of the vorticity equation (5), which we shall now treat as a (linear) passive scalar equation for some scalar quantity θ and a velocity field u that does not depend on θ . Notice that working in cylindrical coordinates would at this point become problematic as the cylindrical divergence of the velocity field u in our fluid dynamics problem might in general be unbounded opposed to the Cartesian divergence, which vanishes identically. In order to apply the DiPerna-Lions theory, in which that boundedness is a crucial assumption, it is therefore advantageous to go back to the Cartesian formulation and rewrite (12) as Both formulations are indeed equivalent if u and θ are axisymmetric, because then If, in addition, u is Sobolev regular, as is the case for the axisymmetric Euler equations under the integrability assumption (7) on the vorticity, the theory in [24] applies. We summarize some of the main results, not aiming for the most general assumptions.

Theorem 4 ([24]
) Let T > 0 and p ∈ (1, ∞) be given and (i) There exists a unique renormalized solution θ ∈ L ∞ ((0, T ); L p (R 3 )) of the transport equation (12) with initial datum θ 0 . (ii) This solution is stable under approximation in the following sense: Let {θ k 0 } k∈N be a sequence that approximates θ 0 in L p (R 3 ) and {u k } k∈N a sequence that approximates u in L 1 ((0, T ); , then distributional solutions are renormalized solutions and vice versa.
It has been proved in [3,24] that renormalized solutions are in fact transported by the (regular) Lagrangian flow of the vector field u, and this feature carries over to the cylindrical setting. Hence, it holds that θ(t, φ(t, x)) = θ 0 (x), where φ satisfies a suitably generalized formulation of the ordinary differential equation In terms of the vorticity, the transport identity can be rewritten as ω(t, φ(t, x)) = ω 0 (x)φ r (t, x)/r , and thus, r /φ r (t, x) is the Jacobian. See also [4] for a review of the connection between the Lagrangian and Eulerian descriptions of transport by nonsmooth velocity fields.
Following [17,18], our strategy for proving that vanishing viscosity solutions to the axisymmetric Euler equations are renormalized solutions relies on duality arguments both in the viscous and in the inviscid setting. In the latter, we quote a suitable duality theorem from DiPerna and Lions's original work.

Estimates on the velocity field
In this section, we provide some estimates on the velocity field that turn out to be helpful in the subsequent analysis. We continue denoting by ω and ξ the vorticity and relative vorticity, respectively, of a given (steady) axisymmetric velocity field u, that is, ω = ∂ z u r − ∂ r u z and ξ = ω/r independently from the Euler or Navier-Stokes background. In particular, any of the following estimates are consequences of the explicit definitions or follow from suitable properties of the Biot-Savart kernel.
We start by verifying that the velocity field generated by (relative) vorticities in the class (9) satisfies the growth condition needed to apply DiPerna and Lions's theory recalled in Theorem 4.

Lemma 2
Let ω ∈ L 1 (H). Then it holds that Proof It is proved in [29] that the axisymmetric Biot-Savart kernel satisfies similar decay estimates as the planar two-dimensional one, namely, if G is obtained from restricting the three-dimensional Biot-Savart kernel to the axisymmetric setting, so that cf. [29, Eq. (2.11)]. We now denote by G 1 the restriction of G to the unit half ball B 1 (0) and set G 2 = G − G 1 , and decompose u = u 1 + u 2 accordingly. Then, on the one hand, by Young's convolution inequality, it holds Thus, in view of the trivial bound u 1 /|x| L 1 (R 3 ) u 1 L 1 (H) , we verify the growth condition for u 1 . On the other hand, since χ B 1 (0) c 1 |·| is uniformly bounded, we deduce a uniform control on u 2 . This establishes (16).
Our second result is a fairly standard identity for the enstrophy, that is, the (square of the) L 2 norm of the velocity gradient.

Lemma 3 It holds that
We provide the argument for this standard identity for the convenience of the reader.
Proof From the definition of the vorticity, we infer that We have to identify the third term on the right-hand side: It holds that Indeed, using the no-penetration boundary condition u r = 0 on ∂H together with the incompressibility condition (4), a multiple integration by parts reveals on the one hand that On the other hand, it holds that It remains to notice that to conclude the statement of the lemma.
In the following lemma, we provide a maximal regularity estimate for the velocity gradient in terms of the relative vorticity. Our proof relies on the classical theories by Calderón, Zygmund and Muckenhoupt.
Proof We note that in view of the Biot-Savart law (6), the velocity gradient can be represented as a singular integral of convolution type, It is well-known that Calderón-Zygmund theory guarantees that for any p ∈ (1, ∞). Our goal is to produce a weighted version of this estimate, namely which is nothing but (18). We are thus led to the theory of Muckenhoupt weights: If p ∈ (1, ∞) and m is in the class of Muckenhoupt weights A p then the weighted-maximal regularity estimate (18) holds. We recall that a nonnegative locally integrable functions m is in the class A p if there exists a constant C > 0 such that for all balls B in R 3 the condition holds with q ∈ (1, ∞) with 1 p + 1 q = 1. This well known result that can be found, for instance, in [25] (Theorem 7.11, Chapter 7).
We thus have to show that m = m(r ) = r − p satisfies (19) for p ∈ (1, 2). For this, consider a ball in R 3 with radius R and centered in a generic point X = (X 1 , X 2 , X 3 ) ∈ R 3 , i.e., B = B R (X ). We denote by d the distance of X to the z-axis, that is, d = X 2 1 + X 2 2 . We split our argumentation into the two cases when d ≥ 2R (far field) and d < 2R (near field).
Let us first consider the case where d ≥ 2R. Notice that we have d − R ≤ for any x ∈ B by the triangle inequality, and thus For Making use of the fact that d+R d−R ≤ 3 for all d ≥ 2R, we deduce that We now turn to the case where d < 2R. We first observe that x 2 1 + x 2 2 < d + R and |x 3 − X 3 | < R for all x ∈ B, and we may thus bound the integral over the ball by an integral over the cylinder. Making relative transformations in cylindrical coordinates, we then have the estimates Taking the product and using that d+R R ≤ 3 for all d < 2R, we conclude that Hence, in either cases, we proved (19) and, thus, the proof is over.
We conclude this section with an estimate on the velocity field in large annuli.
Proof In the following, we will use the notation B R for both the ball in R 3 , that is B R (0) ⊂ R 3 , and the half ball in H, that is B R (0) ∩ H ⊂ H. It should be clear from the situation, which one is considered. We start by using the Sobolev embedding in two dimensions and find To estimate the gradient term in (20), we write the velocity field with the help of a vector stream function, u = ∇ x × ψ, where − x ψ = ωe θ . Then recalling (17), we find that Thanks to local Calderon-Zygmund estimates (see, e.g., Theorem 9.11 in [30]), we have a control on the term on right-hand side in terms of the vorticity, Notice that ψ = ψ θ e θ by construction. Moreover, as the stream function is unique up to additive constants, we may, without loss of generality, assume that ψ θ has zero average in the annulus B 2α R \B β R 2 ⊂ R 3 , so that an application of the Poincaré inequality implies Hence, because we can write u = ∂ z ψ θ e r + 1 r ψ θ + ∂ r ψ θ e z , we have proved that A combination of the previous estimates yields Plugging this estimate into (20), we arrive at With regard to the first term in (21), we notice that by Jensen's inequality, it holds that and, thus, in view of our assumptions on r , the first term in (21) vanishes as R → ∞.
For the second term, we appeal to Hölder's inequality, where we use the fact that r < 2 p 2− p , which holds true because 2 p 2− p > 6 p−2 3− p for any p < 2. We can easily smuggle in some weights to the effect that It remains to observe that the exponent on R is negative by our assumption on r . This shows the convergence of the second term in (21). The proof is over.

Global estimates for the axisymmetric Navier-Stokes equations
In this section, we provide some global estimates for solutions to the Navier-Stokes equations that will turn out to be helpful later on. We start by rewriting the momentum equation (10) in terms of the vorticity ω ν = ∂ z u r ν − ∂ r u z ν and the relative vorticity ξ ν = ω ν /r . The evolution equation for the vorticity is given by and is equipped with homogeneous Dirichlet conditions on the boundary of the halfspace, i.e. ω ν = 0 on ∂H. It follows that the vorticity equation is conservative, as expected, because r −1 ∂ r ω ν − r −2 ω ν = ∂ r (r −1 ω ν ). The relative vorticity satisfies the nonconservative equation which is supplemented with homogeneous Neumann boundary conditions, ∂ r ξ ν = 0 on ∂H. We will mostly work with the latter equation. For initial data ξ ν (0) = ξ 0 in (9), well-posedness for either formulation can be inferred from the theory developed by Gallay and Šverák [29]. In the following, ω ν will always be the unique mild solution to the vorticity equation (22) in the class C([0, T ); L 1 (H)) ∩ C((0, T ); L ∞ (H)) and ξ ν = ω ν /r . We start by recalling some useful properties which can be found in various references. Yet, we provide their short proofs for the convenience of the reader. Our first concern is an L p estimate.

Lemma 6 It holds that
Proof We can perform a quite formal computation as solutions can be assumed to be smooth by standard approximation procedures. A direct calculation yields d dt where we made use of the no-penetration boundary conditions on the velocity field u, i.e., u r = 0 at ∂H, to eliminate the advection term. The Cartesian Laplacian as can be seen by an integration by parts. Another integration by parts reveals that the first order term is nonpositive and can thus be dropped, A combination of the previous estimates yields d dt and an integration in time yields the desired estimate (24).
Our next estimate quantifies integrability improving features of the advectiondiffusion equation (23) by suitably extending the estimates on the L p norm established in the previous lemma to any q ∈ [p, ∞).

Lemma 7 For any q ∈ [p, ∞], it holds that
Proof Our proof is a small modification of the argument of Feng and Šverák in [27,Lemma 3.8], where the case p = 1 is considered. We define E q (t) = ξ ν (t) for some q ∈ [p, ∞) and claim that Let us postpone the proof of this estimate a bit and explain first how it implies (26). Notice that, by interpolation of Lebesgue spaces, it is enough to show (26) for exponents q = 2 k p with k ∈ N 0 and q = ∞. We first treat the case for finite exponents, which will be achieved by induction. We start by observing that the base case k = 0 is settled in Lemma 6 above. The induction step from k to k + 1 is based on estimate (27). We setq = 2 k and q = 2 k+1 = 2q. Plugging (26) withq = q 2 into (27), we find Integrating in time yields where we have used that q > p. Notice that all constants can be chosen uniformly in q. We have thus proved (26) for q = 2q, which settles the case where q = 2 k p.
If q = ∞, we may now simply take the limit in (26) and use the convegence of the Lebesgue norms, · L ∞ = lim q→∞ · L q . It remains to provide the argument for (27). We start by recalling that Notice that the constants in the estimate can be chosen independently of q as q > 1, and can thus be dropped. We estimate the right-hand-side with the help the 3D Nash |ξ ν | q dx 5 3 , which can be rewritten as (27).
We also note that the fluid impulse is conserved along the viscous flow.

Lemma 8 Suppose that r
This identity can be seen in several ways, see, for instance [29,Lemma 6.4] for a proof that is based on the symmetry properties of the Biot-Savart kernel and applies to our regularity setting. We omit the proof and remark only that whenever u is an axisymmetric vector field and ω the associated scalar vorticity. The conservation of momentum follows immediately from the Euler equations (1), (2).
The last global estimate concerns the energy balance law, for which we assume that the initial kinetic energy is bounded.
It is a classical result by Leray that for any divergence-free initial datum u 0 in L 2 (R 3 ), there exists a weak solution to the Navier-Stokes equations (10), (11) satisfying the energy inequality cf. [38]. Whether there is an energy equality (28) for such solutions is an important open problem. There are various conditions available in the literature under which an equality can be established, most notably, Serrin's condition u ∈ L q ((0, T ); L p (R d )) with d p + 2 q ≤ 1 or Shinbrot's criterion 2 p + 2 q ≤ 1 and p ≥ 4, cf. [45,46]. We refer to [15] for an extension of the previous results to a larger class of function spaces and to [5] for a recent improvement based on assumptions on the gradient of the velocity.
It is not difficult to see that we can construct mild solutions in the setting of [29] that satisfy the inequality (29), and thus, thanks to the uniqueness in that setting, our solutions do as well. We remark that in [10] Buckmaster and Vicol construct weak solutions for the three-dimensional Navier for which the energy inequality is not automatically achieved. Unfortunately, it is not obvious how to check Serrin's or Shinbrot's integrability conditions to ensure an energy equality in the axisymmetric setting. The problem is the appearance of weights as, for instance, in (18) and in suitable Sobolev inequalities. For this reason, we provide a proof of (28) that is tailored to our needs but still mimics the original arguments in [45,46].
Proof By interpolation between Lebesgue spaces, we may without loss of generality assume that p ∈ (1, 2). Thanks to the well-posedness result in [29], we may suppose that (29) holds true in our setting. In particular, we deduce In addition, thanks to the L p bound on the vorticity in Lemma 6 and the weighted maximal regularity estimate in Lemma 4, it holds that By standard density arguments, we may thus find a sequence {u δ ν } δ↓0 of axisymmetric divergence-free functions in C ∞ c ((0, T ) × R 3 ) 3 which satisfy (30) and (31) and that converges towards u ν in L 2 ((0, T ); H 1 (R 3 ) 3 ), staying bounded in all the spaces in which u ν is contained. We furthermore denote by η ε a standard mollifier on R. Because is an admissible test function in the definition of distributional solution of the Navier-Stokes equations, we find that In a first step, we send δ to zero with ε > 0 fixed. The convergence is obvious for all but the nonlinear term. It is enough to show that the nonlinear term vanishes when u δ ν is replaced by v δ = u δ ν − u ν . Performing an integration by parts, we can throw the derivative on one of the u ν (t, x). Hölder's inequality then yields where by * we denote the convolution-type operation between η ε and v δ . We now have to make use of the interpolation inequality in Lemma 15 in the appendix and notice that |∇u| ≤ |∇ x u| for any axisymmetric velocity field u. We find that where λ = 3 p−3 7 p−6 . Using Hölder's and Young's convolution inequality, we then infer that Notice that from the definition of v δ , the triangle inequality and estimate (31) applied to u δ ν , it follows that 1 is bounded for all δ. From (30) and (31) and the assumptions on v δ , we deduce that the right-hand side in the above estimate is vanishing as δ → 0. Passing to the limit in the weak formulation of the Navier-Stokes equations above thus yields Notice that the term that involved the time derivative on η ε dropped out by imposing that η ε is an even function. We finally send ε to zero and may thus choose ε < T from here on. Notice first that thanks to standard convergence properties of the mollifier. For the convergence of the end-point integrals, we make use of the fact that our solutions are continuous in time with respect to the weak topology in L 2 (R 3 ), see, e.g., [46,Corollary 3.2]. Because η ε is chosen even, Lebesgue's convergence theorem then yields It remains to argue that the nonlinear term is vanishing. Notice first that for any δ if u δ ν is defined as above. This identity carries over to the limit δ → 0 as can be seen by using the same kind of estimates that we used above in order to control the nonlinear term. We may thus rewrite the nonlinear term above as and, by applying the same kind of estimates again, we observe that this term vanishes as ε → 0 by the convergence properties of the mollifier.

Vanishing viscosity limit. Proof of Theorem 1
In this section, we turn to the proof of Theorem 1. The compactness argument is based on the a priori estimate (24) on the relative vorticity and local estimates on the velocity field. The latter are provided by the following two lemmas.

Lemma 10
For any R > 0, p * ∈ (1, p] ∩ (1, 2) and q * ∈ (2, 2 p 2− p ) ∩ (2, ∞), there exists a constant C(R) such that Proof By standard interpolation between Lebesgue spaces, we may without loss of generality assume that p = p * < 2. The bound on the gradient is an immediate consequence of the maximal regularity estimate in Lemma 4, where B R = B R (0) ∩ H denotes the open half ball of radius R centered at 0 in the half-space H.
In order to deduce an estimate on the velocity field itself, we first invoke the Poincaré estimate for mean-zero functions, formula (17) and the previous bound to observe that uniformly in time. It remains to bound the L 1 norm of u. For this purpose, we make use of the decay behavior of the Biot-Savart kernel. In [29], the authors show that the decay of the axisymmetric Biot-Savart kernel is identical (in scaling) to that of the planar Biot-Savart kernel, that is, if we rewrite (6) For the near-field, we use Fubini's theorem, Young's convolution estimate and Lemma 6 to deduce For the far-field, we simply observe that the kernel is bounded below, and thus Plugging the previous bounds into (33) and recalling the already established gradient bounds yields a uniform in time control on u ν W 1, p (B R ) . From this, we deduce the estimate in L q * in (32) via standard two-dimensional Sobolev embedding.

Lemma 11
For any R > 0, it holds that The proof of this estimate is fairly standard. We sketch the argument for the convenience of the reader.
Proof Let F be a divergence-free vector field in W 1,∞ 0 (B R (0)) 3 . Then and a similar bound holds for the dissipation term −ν u ν . The statement thus follows directly from the momentum equation and Lemma 10.
We are now in the position to prove the compactness result.

Proof of Theorem 1
Notice that the norms in the statement of Lemma 10 can be replaced by the corresponding norms on the three-dimensional balls B R (0), because dx ≤ 2π Rd(r , z) in the domain of integration. Therefore, thanks to Lemmas 10 and 11, the sequence of velocity fields {u ν } ν↓0 satisfies the hypotheses of the Aubin-Lions 0)) and the latter in is continuously embedded in W −1,1 σ (B R (0)). Therefore, for any R > 0, there exists a subsequence that converges strongly in C([0, T ]; L 2 (B R (0))). By applying a diagonal sequence argument, this convergence carries over to the space C([0, T ]; L 2 (K )) for any compact K in R 3 . Hence, there exists a subsequence (not relabelled) and a vector field It is readily checked that u is a distributional solution to the Euler equations (1), (2). Moreover, from the a priori estimate on the relative vorticity in Lemma 6, we deduce that there exists a function ξ ∈ L ∞ ((0, T ); L p (R 3 )) such that, upon taking a further subsequence, We finally notice that the velocity field u and the vorticity ω = r ξ are related by the Biot-Savart law that holds true in the sense of distributions.

Renormalization. Proof of Theorem 2
In this section, we provide the argument for the renormalization property of the relative vorticity obtained as the vanishing viscosity solution of the Navier-Stokes equations in Theorem 1. Our approach is based on the duality formula in Lemma 1 established in [24] and follows closely the argumentation from [17,18]. By interpolation of Lebesques spaces, we may without loss of generality assume that p < 2 in (9).
We now show a compactness result for a backwards advection-diffusion equation, that is, as we will see, dual to the vorticity formulation (23) of the Navier-Stokes equations.
We remark that renormalized solutions to advection-diffusion equations have been considered, for instance, in [24,28,37].

Proof
We start with an a priori estimate. A direct computation reveals that d dt where we have used the Dirichlet boundary conditions on f ν , and thus d dt f ν L q ≥ − χ L q . Via integration and by our choice of the final datum, we deduce that Hence, there exists a subsequence {ν k } k∈N that can be chosen as a subsequence of the one found in Theorem 1 and anf ∈ L ∞ ((0, T ); L q (R 3 )) such that Since at the same time by the virtue of Theorem 1 and q ≥ 2, we find in the limit thatf solves the backward advection equation (15) in the sense of distributions. From Theorem 1 and Lemma 4 we have u ∈ L 1 ((0, T ); W 1, p loc (R 3 )) for 1 < p < 2 and thus, by Theorem 4, the solution is renormalized and thus unique; hencef = f . In particular, the convergence result holds true for the subsequence from Theorem 1.
We finally turn to the proof of the renormalization property.

Proof of Theorem 2
be an arbitrarily fixed axisymmetric function and f ν a solution to the backwards advection-diffusion equation considered in Lemma 12. From the statement of the lemma, it follows that { f ν k } k∈N converges to f weakly-in L ∞ ((0, T ); L q (R 3 )). Moreover, there exists a positive number s such that for any radius R, the sequence of time derivatives {∂ t f ν k } k∈N is bounded in L 1 ((0, T ); H −s (B R (0))), and thus, invoking an Arzelà-Ascoli-type argument, we conclude that the convergence of { f ν k } k∈N can be upgraded to hold in C([0, T ], L q weak (R 3 )), that is, where 1/q + 1/q = 1. Upon a standard approximation argument, f ν can be considered as a test function in the distributional formulation of the vorticity formulation (23) of the Navier-Stokes equations. Thus As a consequence of Theorem 1, Lemma 12 and (36), we can pass to the limit in this identity and find On the other hand, because u satisfies the general assumptions of Theorem 4, see Lemma 2, there exists a unique renormalized solutionξ ∈ L ∞ ((0, T ); L p (R 3 )) to the transport equation (5) with u being the given solution to the Euler equations and with initial datum ξ 0 . Of course,ξ is axisymmetric. By Lemma 1, we then find that and thus, Because χ was an arbitrarily fixed smooth axisymmetric function and ξ andξ are both axisymmetric, we infer thatξ = ξ almost everywhere, and thus, ξ coincides almost everywhere with the renormalized solutionξ .

Energy conservation. Proof of Theorem 3
We now prove Theorem 3. Throughout this section, we thus suppose that ω ν is nonnegative and has finite impulse. One of the main ingredients of the proof is the convergence of the kinetic energy that is established in the following lemma.

Lemma 13
Let p > 11 9 . Let {ν k } k∈N be the subsequence found in Theorem 1. Then it holds that for any t ∈ [0, T ].
Proof We have already seen in Theorem 1 that u ν k converges to u strongly in C(0, T ; L 2 loc (R 3 )). We have to turn this result into a global convergence result. In fact, it is enough to show that Indeed, if (37) holds true, given ε > 0, we can find a radius R ≥ 1 such that Moreover, thanks to the strong convergence in B 2R (0), we have that Combining both estimates, we find that for |h| sufficiently small.
By Riesz' compactness criterion, the latter result together with (37) and the standard energy estimate (29) imply strong convergence in L 2 (R 3 ) for ant t ∈ [0, T ]. We now give the argument for (37). For notational convenience, we write u and ν instead of u ν k and ν k . We consider a smooth cut-off function η R that is 1 in B R = B R (0) and 0 outside B 2R = B 2R (0). Testing the Navier-Stokes equations with (1 − η R ) 2 u and integrating by parts yields d dt The error term in (40) is quite easily estimated. Indeed, using the Cauchy-Schwarz inequality together with the elementary inequality 2ab ≤ εa 2 + 1 ε b 2 , we can absorb the gradient term in (40) in the second term in (38) and we are left with an error term of the form ν R 2 u 2 L 2 (R 3 ) . In view of the energy inequality for the Navier-Stokes equations, this term is obviously vanishing as R → ∞ uniformly in t.
As a next step, we address the first error term in (39). Using the properties of the cut-off function, this term is bounded as follows: z). (41) Here, we have used the same notation for both the ball in R 3 and the half ball in H. It should be clear from the situation, which one is considered. The expression on the right-hand side vanishes as R → ∞ as a consequence of Lemmas 5,6, 8, and 9. We finally turn to the term that involves the pressure, that is, the second term in (39). We choose r ∈ 3, 6 p− p 3−r and assume for convenience that p ≤ 5 3 , which can be achieved by interpolation between Lebesgue spaces. Using the properties of the cut-off function and Hölder's inequality, we observe that Again, the velocity term vanishes by the virtue of Lemmas 5,6, 8, and 9, however, this time we need to extract a small amount of more information on the rate of decay, namely which is the leading order term for large R, because p ≤ 5 3 ≤ 3 2 . It is thus enough to consider the pressure term and to show that lim R→∞ R −1− p(r +6)−3r −2 2( p−1)(r −1) For this, we recall that p solves the Poisson equation − x p = ∇ 2 x : u ⊗ u, and thus, we have that p = i j ∂ x i ∂ x j G * (u i u j ), where G is the Newtonian potential in R 3 , i.e., G(x) = 1 4π 1 |x| . Let us write f = G * (u i u j ), so that | p| |∇ 2 x f |. The localized Calderón-Zygmund estimates (see, e.g., Theorem 9.11 in [30]) yield We observe that 2 pr − 13 p − 9 2( p − 1)r < r p − 8 p + r + 4 4( p − 1)r ⇐⇒ r < 18 p − 14 3 p − 1 , and the latter condition can be achieved by our choice of r 's because 18 p−14 3 p−1 > 3 for any p > 11 9 . To prove (45), we use the kernel representation of f and estimate and decompose Note that in this step, the assumption r > 3 is crucial. Plugging in the bound in (42) .
With these preparations, we are now in the position to prove Theorem 3. Our short proof is strongly inspired by [14].

Proof of Theorem 3
In order to prove conservation of energy, we choose a subsequence as in Theorem 1, which we will not relabel for notational convenience, and recall the energy identity in Lemma 9, which we rewrite as Thanks to Lemmas 3, 8 and 7, we observe that 1 νs