Renormalization and energy conservation for axisymmetric fluid flows

We study vanishing viscosity solutions to the axisymmetric Euler equations with (relative) vorticity in $L^p$ with $p>1$. 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$ 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 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 [37]. 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 [40,41] can be suitably generalized in order to handle this situation and to extend the results from [37,18] 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 [38], 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 [25] for partial results and [13] for improvements, while the sharpness of this exponent was proved in [30], building up on the theory developed in [20,21,8,7,9]. We note that Hölder-1/3 regularity is guaranteed for any vorticity field in L p with p ≥ 9 2 thanks to three-dimensional Sobolev embeddings.
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) if 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 , e r = (cos θ, sin θ, 0) T , e θ = (− sin θ, cos θ, 0) T , e z = (0, 0, 1) T .
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 A transformation into cylindrical coordinates and an analysis of the axisymmetric Biot-Savart law can be found, for instance, in [28]. 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., [22,46,6]. 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 2 We caution the reader that throughout the manuscript, we carefully distinguish between the Lebesgue spaces on the full three-dimensional space, L p (R 3 ), and those on the two-dimensional half-space L p (H). Notice also that the three-dimensional Lebesgue measure reduces to the weighted measure 2πrd(r, z) on H when restricted to axisymmetric configurations as in (7). In particular, Notice that 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. [28,Proposition 2.3]. For vorticity fields with smaller integrability exponents, we propose the notion of renormalized solutions.
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 [24,3] 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 Section 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 is unbounded, more precisely, for some p ∈ (1, ∞). We are thus outside of the class of functions in which uniqueness is known to hold [19,2]. On the positive side, existence of distributional solutions to the Euler equations (1), (2) was proved in [31] 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 [45,23,39,12,11,33,32].
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 (18) (or (19)) below. It was shown in [28] that under the assumption (9) on the initial data, which implies that ω 0 ∈ L 1 (H), there exists a unique global (mild) solution ω ν ∈ C 0 ([0, ∞) × L 1 (H)) ∩ C 0 ((0, ∞) × L ∞ (H)) to the viscous vorticity equation.
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 [31]. 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 [45,33,29,47,1] 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 section 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 6. 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 6, 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 were 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 ∇ = ∂r ∂z , with the interpretation that The advective derivatives f · ∇ and f · ∇ x are to be interpreted correspondingly.
The remainder of the article is organized as follows: In Section 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 Section 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 Sections 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] and explain how it applies to the situation at hand. We are particularly interested into well-posedness and renormalization properties of the transport 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 θ. Yet, we have in mind that u has its origin in the fluid dynamics problem that is considered in the main part of this paper. We shall thus continue assuming that the flow is incompressible in the sense of (4), and that u r = 0 on ∂H, which is ensured in the nonlinear setting by the Biot-Savart law (6), see also the discussion in [28]. Working in cylindrical coordinates becomes at this point problematic as the cylindrical divergence of the velocity field u 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 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, again formulated for the axisymmetric setting, and 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 ); W 1,1 loc (H) 3 ) and such that r −1 ∂ r (ru r k ) + ∂ z u z k = 0. Let θ k denote the corresponding renormalized solution. Then , then distributional solutions are renormalized solutions and vice versa.
It has been proved in [24,3] 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, r, z)) = θ 0 (r, z), 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, r, z)) = ω 0 (r, z)φ r (t, r, z)/r, and thus, r/φ r (t, r, z) is the Jacobian. See also [4] for a review of the connection between the Lagrangian and Eulerian descriptions of transport by non-smooth velocity fields.
It follows a discussion of the validity of the growth condition (13) in the Euler setting.

Remark 1.
In the nonlinear setting in which u can be reconstructed from θ = ξ with the help of the Biot-Savart kernel (6), the growth condition (13) it automatically fullfilled provided that ξ ∈ L ∞ ((0, T ); L 1 (R 3 )) as it is assumed in this paper. Indeed, it is proved in [28] 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. [28,Eq. (2.11)]. We now denote by G 1 the restriction of G to the unit ball B 1 (r, z) 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, On the other hand, Using that ω L 1 (H) = 1 2π ξ L 1 (R 3 ) , we deduce (13). Following [18,17], 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.
Our first result is a fairly standard identity for the enstrophy, that is, the (square of the) L 2 norm of the velocity gradient.

Lemma 2. 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, ∇ x u = K * (ωe θ ), where |K(x)| ∼ 1 |x| 3 . 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 with m = m(r) = 1 r p and r = x 2 1 + x 2 2 , which is nothing but (16). 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 (16) holds. Here, A p is the set of nonnegative locally integrable weight functions satisfying for a universal constant C > 0 and all balls B in R 3 , and q ∈ (1, ∞) with 1 p + 1 q = 1. This well known result that can be found, for instance, in the book of Stein [44].
We thus have to show that m = m(r) = r −p satisfies (17) 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 ≤ x 1 2 + x 2 2 ≤ d + R for any x ∈ B by the triangle inequality, and thus 1 ( For m(x) = ( 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 (17) and, thus, 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 we 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 L 1 (R 3 ) ∩ L p (R 3 ), cf. (9), well-posedness for either formulation can be inferred from the theory developed by Gallay andŠverák [28]. In the following, ω ν will always be the unique mild solution to the vorticity equation (18) 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 4. 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 where we made use of the no-penetration boundary conditions on the velocity field u 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 (20).
Our next estimate quantifies integrability improving features of the advectiondiffusion equation (19) by suitably extending the estimates on the L p norm established in the previous lemma to any q ∈ [p, ∞).

Lemma 5.
For any q ∈ [p, ∞], it holds that Proof. Our proof is a small modification of the argument of Feng andŠverák in [26,Lemma 3.8], where the case p = 1 is considered. We define E q (t) = ξ ν (t) q L q (R 3 ) for some q ∈ [p, ∞) and claim that Let us postpone the proof of this estimate a bit and explain first how it implies (22). Notice that, by interpolation of Lebesgue spaces, it is enough to show (22) 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 4 above. The induction step from k to k + 1 is based on estimate (23). We setq = 2 k and q = 2 k+1 = 2q. Plugging (22) withq = q 2 into (23), 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 (22) for q = 2q, which settles the case where q = 2 k p.
If q = ∞, we may now simply take the limit in (22) and use the convegence of the Lebesgue norms, · L ∞ = lim q→∞ · L q .
It remains to provide the argument for (23). 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 inequality , which can be rewritten as (23).
We also note that the fluid impulse is conserved along the viscous flow. This identity can be seen in several ways, see, for instance [28, 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 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.
for all t > 0.
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. [36]. Whether there is an energy equality (24) 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. [42,43]. 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 [28] that satisfy the inequality (25), 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 (16) and in suitable Sobolev inequalities. For this reason, we provide a proof of (24) that is tailored to our needs but still mimics the original arguments in [42,43].
Proof. Thanks to the well-posedness result in [28], we may suppose that (25) holds true in our setting. In particular, we deduce In addition, thanks to the L p bound on the vorticity in Lemma 4 and the weighted maximal regularity estimate in Lemma 3, it holds that By standard density arguments, we may thus find a sequence {u δ ν } δ↓0 of divergencefree functions in C ∞ c ((0, T ); C ∞ c (R 3 ) 3 ) that converges towards u ν in L 2 ((0, T ); H 1 0 (R 3 ) 3 ) and stays 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 fist 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 13 in the appendix and notice that |∇u| ≤ |∇ x u| for any axisymmetric velocity field u. We find that Using Hölder's and Young's convolution inequality, we then infer that From (26) and (27) 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., [43,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 (20) on the relative vorticity and local estimates on the velocity field. The latter are provided by the following two lemmas.
Lemma 8. For any R > 0 and any p * ∈ (1, p] ∩ (1, 2), there exists a constant C(R) such that where B R (0) is the ball in R 3 centered in the origin.
Proof. By standard interpolation between Lebesgue spaces, we may without loss of generality assume that p = p * . The bound on the gradient is an immediate consequence of the maximal regularity estimate in Lemma 3 and formula (15), In order to deduce an estimate on the velocity field itself, we first invoke the Poincaré estimate for mean-zero functions 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 [28], 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) as u ν (r, z) = H K(r, z,r,z)ω ν (r,z) d(r,z), then the kernel K obeys the estimate |K(r, z,r,z)| 1 |r −r| + |z −z| .
We thus have and write For the near-field, we use Fubini's theorem, Young's convolution estimate and Lemma 4 to deduce For the far-field, we simply observe that the kernel is bounded below, and thus Plugging the previous bounds into (29) yields (28) as desired.
Lemma 9. For any R > 0, it holds that where W −1,1 σ (B R (0)) 3 is the Banach space that is dual to the space of divergence-free vector fields in W 1,∞ 0 (B R (0)) 3 .
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 8.
We are now in the position to prove the compactness result.
Proof of Theorem 1. Thanks to Lemmas 8 and 9, the sequence of velocity fields {u ν } ν↓0 satisfies the hypotheses of the Aubin-Lions Lemma, and thus, 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 u ∈ C([0, T ]; L 2 loc (R 3 ) 3 ) such that 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 4, 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 [18,17].
In a first step, we show a compactness result for a backwards advection-diffusion equation, that is, as we will see, dual to the vorticity formulation (19) of the Navier-Stokes equations.
We remark that renormalized solutions to advection-diffusion equations have been considered, for instance, in [24,35,27].
Proof. We start with an a priori estimate. A direct computation reveals that and thus, by a Gronwall argument and our choice of the final datum, 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 and q ≥ 2, we find in the limit thatf solves the backward advection equation (14), and thus,f = f by uniqueness. 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. Let χ ∈ C ∞ c ((0, T )×H) be given and f ν a solution to the backwards advection-diffusion equation considered in Lemma 10. 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 )) for any q ∈ (2, ∞). By using the advection-diffusion equation, this convergence can be upgraded to hold in C([0, T ], L q weak (R 3 )) for any q ∈ (2, ∞), that is, where 1/q + 1/q = 1. Indeed, it is not difficult to obtain this result for smooth test functions and the full statement is obtained by standard approximation procedures. Upon a standard approximation argument, f ν can be considered as a test function in the distributional formulation of the vorticity formulation (19) of the Navier-Stokes equations. Thus As a consequence of Theorem 1, Lemma 10 and (30), we can pass to the limit in this identity and find On the other hand, because u satisfies the assumptions of Theorem 4, see also Remark 1, there exists a unique distributional 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 . By Lemma 1, we then find that and thus, Because χ was arbitrarily fixed, we infer thatξ = ξ almost everywhere, and thus, ξ coincide 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. Moreover, we assume that p > 3 2 as in the assumption of the theorem. Notice that by interpolation between Lebesgue spaces, we may always suppose that p ∈ 3 2 , 2 , which we will do from here on. One of the main ingredients of the proof is the convergence of the kinetic energy that is established in the following lemma.
Lemma 11. Let {ν k } k∈N be the subsequence found in Theorem 1. Then it holds that lim 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 (31) holds true, given ε > 0, we can find a radius R ≥ 1 such that sup k u ν k (t, · + h) L 2 (R 3 \B 2R (0)) ≤ ε for any |h| ≤ 1.
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 (31) and the standard energy estimate (25) imply strong convergence in L 2 (R 3 ) for ant t ∈ [0, T ]. We now give the argument for (31). 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 (34) 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 (34) in the second term in (32) 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 (33). Using the properties of the cut-off function, this term is bounded as follows: Here, we have used the same notation for both the ball in R 3 and the one in H. It should be clear from the situation, which one is considered. We now use the Sobolev embedding in two dimensions and find Moreover, thanks to local Calderon-Zygmund estimates (which we can perform on the level of the standard Biot-Savart kernel in R 3 ), we find that With regard to the first term in (36), we notice that by Jensen's inequality and the energy inequality, it holds that and, thus, the first term in (36) vanishes as R → ∞, uniformly in t. For the second term in (36), we appeal to Hölder's inequality, .
We can easily smuggle in some weights, It remains to observe that the exponent on R is negative because p > 3 2 > 4 3 . Combining the above estimates, we conclude that uniformly in t, and thus, in view of (35), the first term in (33) vanishes.
We finally turn to the term that involves the pressure, that is, the second term in (33). Using the properties of the cut-off function and Hölder's inequality, we observe that In view of (37), it is enough to show that the pressure term is bounded, in the sense that 1 uniformly in R, ν and t. To establish this estimate, we recall that p solves the Poisson equation −∆p = ∇ 2 : u ⊗ u, and thus, we have that where G is the Newtonian potential in R 3 , G(z) = 1 4π 1 |z| . Let us write f = G * (u i u j ). The localized Calderón-Zygmund estimates yield The first term is controlled thanks to (37). (Notice that the exact value of the radii in (37) is not of importance but the scaling in R.) Now, since for any m ∈ 3 2 , ∞ , it holds that by Hölder's inequality, and R 1 3 − 3 m → 0 as R → ∞ for m < 9, it suffices to show that ∇f L m (R 3 ) is bounded for some m ∈ 3 2 , 9 . Notice first that the standard Sobolev inequality in R 3 yields as long as s = 3m 3+m ∈ [1, 3). By our choice of m, we have to restrict the range of admissible s to the interval 1, 9 4 . Now we use the maximal regularity properties of the Laplacian, in the sense that In order to estimate the velocity field in L 2s , we use the Sobolev inequality in three dimensions and Calderón-Zygmund estimates for the gradient of the Biot-Savart kernel, .
Let us now write q = 6s 2s+3 and notice that q ∈ 6 5 , 27 15 by our choice of s. Recall from our line of proof that we have to show the boundedness of ω L q (R 3 ) for some value of q in the above interval. This is achieved via interpolation of the estimates in Lemmas 4 and 6. Indeed, setting θ = q−1 p−1 for some q ∈ 6 5 , 3 2 whose explicit value we will specify in a moment, we have that where for the second to last inequality we used Hölder estimate with exponents 1/θ and 1 1−θ and in the last identity we used the definition of θ. We may now determine q by requiring that q(p−1) p−q = 2, which yields q = 2p p+1 . It remains to notice that this choice of q is admissible because 2p p+1 ∈ 6 5 , 4 3 for any p ∈ 3 2 , 2 which we may assume here as explained in the introduction to this section. By using the a priori bounds in Lemmas 4 and 6 we conclude that (38) holds uniformly in ν and t. We have thus established Lemma 11.
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 7, which we rewrite as Thanks to Lemmas 2, 6 and 5, we observe that 1 νs

2p
, and thus, the energy identity implies that because p > 3 2 . Sending ν to zero, we conclude that lim ν→∞ u ν (t) L 2 (R 3 ) = u 0 L 2 (R 3 ) , and the statement of the theorem follows upon applying Lemma 11, in which the convergence of the kinetic energy is established.

Appendix: Two auxiliary inequalities
We conclude this paper with two auxiliary inequalities, that are weighted versions of standard Sobolev and interpolation inequalities.
Then H |f | t r α d(r, z) This estimate is proved, for instance, in [34]. We recall the argument for completeness.
We set γ = α t and let g ∈ C ∞ c (H) be defined by g(r, z) = f (2r, z) and A = [R, 2R]×R and B = [R, 4R] × R be two subsets of H for some R > 0 fixed. By Hölder's inequality, we then have that We now use the embeddingẆ 1,1 ⊂ L ∞ , that holds true in one space dimension, in each variable. On the one hand, using the embedding in r (in form of the fundamental theorem of calculus), we have On the other hand, it holds that where we have used the triangle inequality and a rescaling argument in the last inequality. Combining the previous three estimates, we find that |∇f |r γ d(r, z).
Our next goal is the Hardy-type inequality A |f − g|r γ d(r, z) B |∇f |r γ+1 d(r, z), which holds true provided that γ > −1, and thus α + t > 0. It can be established as follows. Using the fundamental theorem again, we observe that which implies (40) because the prefactor R γ+1 can be smuggled into the integrand.
Towards the weighted Sobolev inequality with s = 1, we set A k = [2 k , 2 k+1 ] × R and B k = [2 k , 2 k+2 ] × R and estimate with the help of the triangle inequality H |f − g| t r α d(r, z) Interpolation between Lebesgue spaces and an application of (39) and (40)  for some universal constant C. It remains to apply the triangle inequality and a change of variables to the effect that