Next Article in Journal
A Novel Class of Separation Axioms, Compactness, and Continuity via C-Open Sets
Next Article in Special Issue
On a Generalized Gagliardo–Nirenberg Inequality with Radial Symmetry and Decaying Potentials
Previous Article in Journal
Composing Diverse Ensembles of Convolutional Neural Networks by Penalization
Previous Article in Special Issue
Superiorization with a Projected Subgradient Algorithm on the Solution Sets of Common Fixed Point Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Global Regular Axially Symmetric Solutions to the Navier–Stokes Equations: Part 1

by
Wojciech M. Zaja̧czkowski
1,2
1
Institute of Mathematics, Polish Academy of Sciences, Śniadeckich 8, 00-656 Warsaw, Poland
2
Institute of Mathematics and Cryptology, Cybernetics Faculty, Military University of Technology, S. Kaliskiego 2, 00-908 Warsaw, Poland
Mathematics 2023, 11(23), 4731; https://doi.org/10.3390/math11234731
Submission received: 28 September 2023 / Revised: 10 November 2023 / Accepted: 13 November 2023 / Published: 22 November 2023
(This article belongs to the Special Issue Variational Problems and Applications, 2nd Edition)

Abstract

:
The axially symmetric solutions to the Navier–Stokes equations are considered in a bounded cylinder Ω R 3 with the axis of symmetry. S 1 is the boundary of the cylinder parallel to the axis of symmetry and S 2 is perpendicular to it. We have two parts of S 2 . For simplicity, we assume the periodic boundary conditions on S 2 . On S 1 , we impose the vanishing of the normal component of velocity, the angular component of velocity, and the angular component of vorticity. We prove the existence of global regular solutions. To prove this, it is necessary that the coordinate of velocity along the axis of symmetry vanishes on it. We have to emphasize that the technique of weighted spaces applied to the stream function plays a crucial role in the proof of global regular axially symmetric solutions. The weighted spaces used are such that the stream function divided by the radius must vanish on the axis of symmetry. Currently, we do not know how to relax this restriction. In part 2 of this topic, the periodic boundary conditions on S 2 are replaced by the conditions that both the normal component of velocity and the angular component of vorticity must vanish. Moreover, it is assumed that the normal derivative of the angular component of velocity also vanishes on S 2 . A transformation from part 1 to part 2 is not trivial because it needs new boundary value problems, so new estimates must be derived.

1. Introduction

The regularity problem for axially symmetric solutions to the Navier–Stokes equations has a long history. However, there are only two results where the global regular axially symmetric solutions are proved, assuming the vanishing of the angular component of velocity (see papers [1] by O.A. Ladyzhenskaya and [2] by M.R. Ukhovskii and V.I. Yudovich).
Other results (see the papers cited in [3,4,5,6,7]) describe the existence of global regular axially symmetric solutions imposing different Serrin-type conditions. The conditions are such that certain coordinates, either of velocity or of derivatives of velocity or vorticity, belong to L q ( 0 T ; L p ( Ω ) ) spaces for appropriately chosen parameters p and q.
This paper closely aligns with the results presented by O.A. Ladyzhenskaya, M.R. Ukhovskii, and V.I. Yudovich, as the vanishing of the stream function divided by the radius implies the existence of global regular axially symmetric solutions. The aim of this paper is to provide a proof of the global estimate (24).
The estimate can imply any global regularity of solutions to problem (6), assuming appropriate regularity of data.
We must emphasize that the methods and proofs presented in this paper are completely new. The proofs and results in Section 3, Section 5, and Section 6 are original.
Before the formal introduction starts, we outline the main steps of the proof of Theorem 1. The main difficulty in the regularity theory of the Navier–Stokes equations lies in handling the nonlinear terms. We need to transform them in such a way that they can be absorbed by the main linear terms. In this paper, we consider problems (17)–(20) for functions Φ and Γ defined by (16). Applying the energy method, we derive inequality (111) with a strongly nonlinear term denoted by I 3 .
The main task of this paper is to estimate I 3 by quantities that can be absorbed by the terms from the l.h.s. of (111).
I 3 is estimated in (122). Using notation (132), we derive from (111) and (122) the inequality (see (134))
(∗)
X 2 ϕ 1 X 2 δ + ϕ ( data ) ,
where ϕ 1 depends on | v φ | d , , Ω t , | v φ | , Ω t and d > 3 , δ > 0 .
For δ > 0 , the Young inequality can be applied, so (133) holds. The existence of such a positive δ follows from inequality (173), which can be written in the following form:
(∗∗)
| Φ | 2 , Ω t 2 ϕ 2 Γ 1 , 2 , Ω t + ϕ ( data ) ,
where ϕ 2 depends on | v φ | , Ω t .
Inequality ( ) implies the existence of the positive δ . For δ = 0 , we were not able to apply the Young inequality in ( ) , so we could not prove Theorem 1.
Hence, ( ) is probably the most important inequality in this paper. It is a totally new result.
In the next step, we eliminate | v φ | d , , Ω t , d = 12 . To show this, we need to delve into the proof of Lemma 13. To derive (141) from (140), we need the estimate
(∗∗∗)
Ω ψ 1 2 r 6 d x c Γ 1 , Ω 2 .
The Hardy inequality implies that ( ) does not hold for 6 but holds for any number less than 6. It is denoted by 6 . Inequality ( ) follows from (202).
Then, we derive (145). Using (133) in (145) yields the following inequality:
(∗∗∗∗)
| v φ | 12 , , Ω t 6 c | v φ | 12 , , Ω t 4 ε θ + ϕ ( data ) .
To apply the Young inequality in ( ) , we require that 6 > 4 ε θ . It is shown in Remark 4 that the inequality holds. We need 6 to be close to 6. Then, ( ) implies (137).
Moreover, to prove (137), we need the existence of such solutions to problem (6), where v φ is not small. The existence of such solutions is proven in Appendix A. Hence, for such local solutions, we prove global estimate (24). Once we have (24), we can extend the local solution incrementally over time.
Finally, we can easily derive estimate (152) because | v φ | , . Ω t appears in (137) with arbitrarily small power.
Using estimates (137) and (152) in (133) implies (24) and proves Theorem 1.
In this paper, we prove the existence of global regular axially symmetric solutions to the Navier–Stokes equations in a cylindrical domain Ω R 3 :
Ω = { x R 3 : x 1 2 + x 2 2 < R 2 , | x 3 | < a } ,
where a and R are the given positive numbers. We denote by x = ( x 1 , x 2 , x 3 ) the Cartesian coordinates. It is assumed that the x 3 -axis is the axis of symmetry of Ω .
Moreover,
S 1 = { x R 3 : x 1 2 + x 2 2 = R , x 3 ( a , a ) } , S 2 = S 2 ( a ) S 2 ( a ) and S 2 ( a 0 ) = { x R 3 : x 1 2 + x 2 2 < R , x 3 = a 0 { a , a } } ,
where S 1 is parallel to the axis of symmetry and S 2 ( a 0 ) is perpendicular to it. S 2 ( a 0 ) meets the axis of symmetry at a 0 .
To describe the considered problem, we introduce cylindrical coordinates r, φ , and z by the relations
x 1 = r cos φ , x 2 = r sin φ , x 3 = z .
The following orthonormal system
e ¯ r = ( cos φ , sin φ , 0 ) , e ¯ φ = ( sin φ , cos φ , 0 ) , e ¯ z = ( 0 , 0 , 1 )
is connected to the cylindrical coordinates.
Any vector, u, for the axially symmetric motions can be decomposed as follows:
u = u r ( r , z , t ) e ¯ r + u φ ( r , z , t ) e ¯ φ + u z ( r , z , t ) e ¯ z ,
where u r , u φ , and u z are cylindrical coordinates of u.
Therefore, velocity v and vorticity ω = rot v are decomposed in the form
v = v r ( r , z , t ) e ¯ r + v φ ( r , z , t ) e ¯ φ + v z ( r , z , t ) e ¯ z
and
ω = ω r ( r , z , t ) e ¯ r + ω φ ( r , z , t ) e ¯ φ + ω z ( r , z , t ) e ¯ z .
The paper is devoted to the proof of global regular axially symmetric solutions to the problem
v , t + v · v ν Δ v + p = f in Ω T = Ω × ( 0 , T ) , div v = 0 in Ω T , v satisfies periodic boundary conditions on S 2 T = S 2 × ( 0 , T ) , v · n ¯ | S 1 = 0 , ω φ | S 1 = 0 , v φ | S 1 = 0 on S 1 T = S 1 × ( 0 , T ) , v | t = 0 = v ( 0 ) in Ω ,
where v = ( v 1 ( x , t ) , v 2 ( x , t ) , v 3 ( x , t ) ) R 3 is the velocity of the fluid, p = p ( x , t ) R is the pressure, f = ( f 1 ( x , t ) ,   f 2 ( x , t ) ,   f 3 ( x , t ) ) R 3 is the external force field, and ν > 0 is the constant viscosity coefficient.
Expressing problem (6) in the cylindrical coordinates of velocity yields
v r , t + v · v r v φ 2 r ν Δ v r + ν v r r 2 = p , r + f r , v φ , t + v · v φ + v r r v φ ν Δ v φ + ν v φ r 2 = f φ , v z , t + v · v z ν Δ v z = p , z + f z , ( r v r ) , r + ( r v z ) , z = 0 v r | S 1 = 0 , v φ | S 1 = 0 , v r , z v z , r | S 1 = 0 , v r | t = 0 = v r ( 0 ) , v φ | t = 0 = v φ ( 0 ) , v z | t = 0 = v z ( 0 ) ,
where we have the periodic boundary conditions on S 2 and
v · = ( v r e ¯ r + v z e ¯ z ) · = v r r + v z z , Δ u = 1 r ( r u , r ) , r + u , z z .
Formulating problem (6) in terms of the cylindrical coordinates of vorticity implies
ω r , t + v · ω r ν Δ ω r + ν ω r r 2 = ω r v r , r + ω z v r , z + F r , ω φ , t + v · ω φ v r r ω φ ν Δ ω φ + ν ω φ r 2 = 2 r v φ v φ , z + F φ , ω z , t + v · ω z ν Δ ω z = ω r v z , r + ω z v z , z + F z , ω r | t = 0 = ω r ( 0 ) , ω φ | t = 0 = ω φ ( 0 ) , ω z | t = 0 = ω ( 0 )
and we have boundary conditions (7)5 on S 1 and the periodic boundary conditions on S 2 , where F = rot f and
F = F r ( r , z , t ) e ¯ r + F φ ( r , z , t ) e ¯ φ + F z ( r , z , t ) e ¯ z .
The function
u = r v φ
is called swirl. It is a solution to the problem
u , t + v · u ν Δ u + 2 ν r u , r = r f φ f 0 , u | S 1 = 0 and u satisfies periodic boundary conditions on S 2 , u | t = 0 = u ( 0 ) .
The cylindrical components of vorticity can be described in terms of the cylindrical components of the velocity and swirl in the following form
ω r = v φ , z = 1 r u , z , ω φ = v r , z v z , r , ω z = 1 r ( r v φ ) , r = v φ , r + v φ r = 1 r u , r .
Equation ( 7 ) 4 implies the existence of the stream function ψ , which is a solution to the problem
Δ ψ + ψ r 2 = ω φ , ψ | S 1 = 0 , ψ satisfies periodic boundary conditions on S 2 .
Moreover, cylindrical components of velocity can be expressed in terms of the stream function in the following way
v r = ψ , z , v z = 1 r ( r ψ ) , r = ψ , r + ψ r , v r , r = ψ , z r , v r , z = ψ , z z , v z , z = ψ , r z + ψ , z r , v z , r = ψ , r r + 1 r ψ , r ψ r 2 .
Introduce the pair
( Φ , Γ ) = ( ω r / r , ω φ / r ) .
Formula (6) from [8] implies that quantities (16) satisfy the following equations
Φ , t + v · Φ ν Δ + 2 r r Φ ( ω r r + ω z z ) v r r = F r / r F ¯ r
and
Γ , t + v · Γ ν Δ + 2 r r Γ + 2 v φ r Φ = F φ / r F ¯ φ .
We add the following initial and boundary conditions to solutions of (17) and (18)
Φ | S 1 = 0 , Γ | S 1 = 0 , Φ , Γ satisfy the periodic boundary conditions on S 2 ,
Φ | t = 0 = Φ ( 0 ) , Γ | t = 0 = Γ ( 0 ) .
Next, we express cylindrical coordinates of velocity in terms of ψ 1 = ψ / r
v r = r ψ 1 , z , v z = ( r ψ 1 ) , r + ψ 1 = r ψ 1 , r + 2 ψ 1 , v r , r = ψ 1 , z r ψ 1 , r z , v r , z = r ψ 1 , z z , v z , z = r ψ 1 , r z + 2 ψ 1 , z , v z , r = 3 ψ 1 , r + r ψ 1 , r r .
The aim of this paper is to prove the existence of global regular axially symmetric solutions to problem (6). For this purpose, we have to find a global estimate that guarantees the existence of global regular solutions.
Function ψ 1 is a solution to the problem
Δ ψ 1 2 r ψ 1 , r = ω 1 in Ω = ( 0 , R ) × ( a , a ) , ψ 1 | r = R = 0 , ψ 1 satisfies the periodic boundary conditions on S 2 ,
where
ω 1 = ω φ / r .
We have ω 1 = Γ .
To state the main result, we first introduce the necessary assumptions.
Assumption 1.
Assume that the following quantities are finite:
D 1 = f L 2 ( Ω t ) + v ( 0 ) L 2 ( Ω ) , D 2 = f 0 L , 1 ( Ω t ) + u ( 0 ) L ( Ω ) , f 0 = r f φ , u = r v φ , D 3 2 = D 1 2 D 2 2 + u , z ( 0 ) L 2 ( Ω ) 2 + f 0 L 2 ( Ω t ) 2 , D 4 2 = D 1 2 ( 1 + D 2 ) + u , r ( 0 ) L 2 ( Ω ) 2 + f 0 L 2 ( Ω t ) 2 + f 0 L 2 ( 0 , t ; L 4 / 3 ( S 1 ) ) ,
where D 1 and D 2 are introduced in (46) and (52), respectively, and D 3 and D 4 are introduced in (159) and (160), respectively. Let
D 5 = D 2 ( D 1 + D 2 + D 3 ) , D 6 = D 2 1 ε 0 D 3 ,
where ε 0 is an arbitrarily small positive number. Moreover,
D 7 = F r L 2 ( 0 , t ; L 6 / 5 ( Ω ) ) 2 + F z L 2 ( 0 , t ; L 6 / 5 ( Ω ) ) 2 + ω r ( 0 ) L 2 ( Ω ) 2 + ω z ( 0 ) L 2 ( Ω ) 2
is defined in Lemma 16.
Next,
D 8 = ϕ ( D 2 ) ( F ¯ r L 2 ( 0 , t ; L 6 / 5 ( Ω ) ) 2 + F ¯ φ L 2 ( 0 , t ; L 6 / 5 ( Ω ) ) 2 ) + Φ ( 0 ) L 2 ( Ω ) 2 + Γ ( 0 ) L 2 ( Ω ) 2 ,
where F ¯ r = F r / r , F ¯ φ = F φ / r , Φ = ω r r , Γ = ω φ r and D 8 appear in (111).
In Lemma 13, the following quantity is defined
D 9 ( 12 ) = 12 f φ L 12 ( 0 , t ; L 36 / 25 ( Ω ) ) + v φ ( 0 ) L 12 ( Ω ) .
Finally, in Lemma 14, we introduce the quantity
D 10 = f φ / r L 1 ( 0 , t ; L ( Ω ) ) + v φ ( 0 ) L ( Ω ) .
The main result is as follows:
Theorem 1.
Assume that Assumption 1 holds. Then, there exists an increasing positive function ϕ, such that
Φ V ( Ω t ) + Γ V ( Ω t ) ϕ ( D 1 , , D 10 ) .
Remark 1.
Estimate (24) implies any regularity of solutions to problem (6), assuming sufficient regularity of data.
To prove (24), we require that ψ 1 and v z vanish on the axis of symmetry.
Proof of Theorem 1. 
Inequality (113) in the form
d d t | Φ | 2 , Ω 2 + | Φ | 2 , Ω 2 I + Ω F ¯ r Φ d x
is the first step of the proof of (24), where Φ = v φ , z r , F ¯ r = F r r and
I Ω | v φ r v r r Φ , z | d x + Ω | v φ z v r r Φ , r | d x I 1 + I 2 .
Our aim is to estimate I 1 and I 2 by a product of norms Φ V ( Ω t ) , Γ V ( Ω t ) .
Since the L -estimate of swirl r v φ is bounded by D 2 (see Lemma 2) and v r r = ψ 1 , z , we obtain the estimates
I 1 D 2 | ψ 1 , r z r | 2 , Ω | Φ , z | 2 , Ω , I 2 D 2 | ψ 1 , z z r | 2 , Ω | Φ , r | 2 , Ω .
To examine estimate (26), we recall that ψ 1 is a solution to problem (22).
In Lemma 4, we prove the existence of weak solutions to problem (22) and derive the estimate (56)
ψ 1 1 , Ω c | w 1 | 6 / 5 , Ω .
In Section 3, we increase the regularity of weak solutions by deriving estimates for higher derivatives.
From (82), we have
| ψ 1 , r z r | 2 , Ω c | Γ , z | 2 , Ω .
The estimate holds for the weak solutions to problem (22) because [9] yields the expansion of ψ 1 near the axis of symmetry
ψ 1 = a 1 ( z , t ) + a 2 ( z , t ) r 2 + a 3 ( z , t ) r 4 +
Hence, ψ 1 , r = 2 a 2 ( z , t ) r and the norm | ψ 1 , r z r | 2 , Ω can be finite.
To estimate I 2 , we need
| ψ 1 , z z r | 2 , Ω c | Γ , z | 2 , Ω .
The estimate holds for such a class of regularized weak solutions to problem (22), where
ψ 1 | r = 0 = 0 .
This means that in expansion (29), we have a 1 ( z , t ) = 0 .
The existence of solutions to problem (22) (also see (61)), satisfying restriction (31) and estimate (30), follows from the theory developed by Kondratiev (see [10]) for elliptic boundary value problems in domains with cones in weighted Sobolev spaces.
In this paper, the existence is proven in Lemmas 8 and 17. From [10], it also follows that we can prove the existence of different solutions to problem (22) belonging to different weighted Sobolev spaces.
The difference between such two distinct solutions is equal to the expression derived from the Cauchy theorem for complex functions related to contour integration. This is described in more detail in [11].
Restriction (31) means that we have to work with a very restricted class of weak solutions to (22). This also means that v z must vanish on the axis of symmetry.
Using estimates (28) and (30) in (25) yields
d d t | Φ | 2 , Ω 2 + | Φ | 2 , Ω 2 c D 2 | Γ , z | 2 , Ω | Φ | 2 , Ω + Ω F ¯ r Φ d x .
We have to emphasize that we are not able to prove estimate (24) without restriction (31).
Now, we integrate (120) with respect to time. Then, we obtain
| Γ | 2 , Ω 2 + Γ 1 , 2 , Ω t 2 2 | Ω t v φ r Φ Γ d x d t | + c | F ¯ φ | 6 / 5 , 2 , Ω t 2 + c | Γ ( 0 ) | 2 , Ω 2 .
Integrating (32) with respect to time and adding to (33) yields
Φ V ( Ω t ) 2 + Γ V ( Ω t ) 2 c ( D 2 ) | Ω t v φ r Φ Γ d x d t | + c ( D 2 ) ( | F ¯ r | 6 / 5 , 2 , Ω t 2 + | F ¯ φ | 6 / 5 , 2 , Ω t 2 ) + c ( D 2 ) ( | Φ ( 0 ) | 2 , Ω 2 + | Γ ( 0 ) | 2 , Ω 2 ) .
Now, we have to estimate the first term on the r.h.s. of (34).
Introducing the quantity (see (132))
X ( t ) = Φ V ( Ω t ) + Γ V ( Ω t )
and recalling that constant D 8 is introduced in Assumption 1, inequality (34) takes the form
X 2 ( t ) c ( D 2 ) | Ω t v φ r Φ Γ d x d t | + c D 8 2 ,
where the first integral is called I 3 .
Using estimate (123) and estimate of L 1 2 in the proof of Lemma 11, we obtain from (36) the inequality
X 2 ( t ) = c ( D 2 ) | v φ | d , , Ω t ε | Φ | 2 , Ω t θ | Φ | 2 , Ω t 1 θ | Γ | 2 , Ω t + c D 8 2 ,
where θ = 1 3 d ε 1 3 d ε 2 , d > 3 , ε = ε 1 + ε 2 < 1 .
To derive any estimate from (37), we use (173) in the form
| Φ | 2 , Ω 2 c ( D 5 + D 6 | v φ | , Ω t ε 0 ) Γ 1 , 2 , Ω t + c D 7 ,
where ε 0 can be assumed to be an arbitrarily small positive number and D 5 , D 6 , and D 7 are defined in Assumption 1. This is a very important estimate because the square of | Φ | 2 , Ω t depends linearly on Γ 1 , 2 , Ω t .
Using (38) in (37) yields the following (the estimate of I 3 is described in (122)):
X 2 ( t ) c | v φ | d , , Ω t ε [ c 1 ( 1 + | v φ | , Ω t 1 2 θ ε 0 ) X 1 2 θ + c 2 ] X 2 θ + c D 8 2 ,
where c 1 and c 2 depend on D 5 , D 6 , and D 7 .
Since 2 1 2 θ , 2 θ are less than 2 Lemma 12 yields the inequality
X 2 c 0 | v φ | d , , Ω t 4 ε θ ( 1 + | v φ | , Ω t 2 ε 0 ) + c 0 | v φ | d , , Ω t 2 ε θ + c D 8 2 ,
where c 0 = ϕ ( D 2 , D 5 , D 6 , D 7 ) .
Setting d = 12 and assuming that v φ is not small, we derive (137) in the form
| v φ | 12 , , Ω t 6 c | v φ | , Ω t b 0 ε 0 + ϕ ( D 2 , D 5 , D 6 , D 7 , D 8 , D 9 ) ,
where b 0 is a positive number.
The smallness of v φ , which must be excluded in the proof of (41), is described in Appendix A.
To prove (41), we have to pass from (140) to (141). Therefore, we need the following estimate:
Ω t ψ 1 2 r 6 d x d t c Γ 1 , 2 , Ω t 2 ,
where 6 < 6 and 6 are very close to 6. Moreover, 6 is such a number where (42) holds (also see Remark 8).
Estimate (42) is crucial to the proof of (141), which is very important in deriving (151).
Inequalities (151) and (152) imply the main result of this paper: estimate (24).
Replacing 6 with 6 estimate (42) takes the form
0 t Ω ψ 1 2 r 6 d x d t c 0 t Γ H 0 1 ( Ω ) 2 d t ,
where the r.h.s. cannot be estimated by Γ V ( Ω t ) .
Estimate (42) follows from Lemma 18 and imposes the following additional restrictions on ψ 1 :
ψ 1 | r = 0 = 0 , ψ 1 , r | r = 0 = 0 .
However, the theory developed in [9] implies that ψ 1 , r | r = 0 = 0 .
Exploiting (41) in (40) yields
X c ( 1 + | v φ | , Ω t d 1 ε 0 ) | v φ | , Ω t d 2 ε 0 + ϕ ( D 2 , D 5 , D 6 , D 7 , D 8 , D 9 ) ,
where d 1 and d 2 are positive finite numbers.
Finally, we find the estimate for | v φ | , Ω t (see (152)). Using (152) in (45) yields (24). This ends the proof of Theorem 1. □
The problem of regularity of axially symmetric solutions to the Navier–Stokes equations has a long history. The first regularity results in the case of vanishing swirl were independently derived by O. A. Ladyzhenskaya and Ukhovskii-Yudovich, as referenced in [1,2]. Many references concerning the nonvanishing swirl case can be found in [3].
We have to emphasize that we were able to prove Theorem 1 because the theory of weighted Sobolev spaces developed in [11] was used.

2. Notation and Auxiliary Results

First, we introduce some notations
Definition 1.
We use the following notation for Lebesgue and Sobolev spaces
u L p ( Ω ) = | u | p , Ω , u L p ( Ω t ) = | u | p , Ω t , u L p , q ( Ω t ) = u L q ( 0 , t ; L p ( Ω ) ) = | u | p , q , Ω t ,
where p , q [ 1 , ] . Next,
u H s ( Ω ) = u s , Ω , u W p s ( Ω ) = u s , p , Ω , u L q ( 0 , t ; W p k ( Ω ) ) = u k , p , q , Ω t , u k , p , p , Ω t = u k , p , Ω t ,
where s , k N { 0 } , H s ( Ω ) = W 2 s ( Ω ) .
We need energy-type space V ( Ω t ) , which is appropriate for the description of weak solutions to the Navier–Stokes equations
u V ( Ω t ) = | u | 2 , , Ω t + | u | 2 , Ω t .
We recall weighted Sobolev spaces, defined by
f H μ k ( R + ) = R + j = 0 k | r j f | 2 r 2 ( μ + j k ) r d r 1 / 2
and
f H μ k ( Ω ) = Ω | α | = 0 k | D r , z α f | 2 r 2 ( μ + | α | k ) r d r d z 1 / 2 ,
where Ω contains the axis of symmetry, D r , z α = r α 1 z α 2 , | α | = α 1 + α 2 , α i N { 0 } , i = 1 , 2 , k N { 0 } , and μ R + . Moreover, we have
H 0 0 ( Ω ) = L 2 , 0 ( Ω ) = L 2 ( Ω ) , H μ 0 ( Ω ) = L 2 , μ ( Ω )
and
f L 2 , μ ( Ω ) = | f | 2 , μ , Ω .
Lemma 1.
Let f L 2 , 1 ( Ω t ) , v ( 0 ) L 2 ( Ω ) . Then, solutions to (7) satisfy the estimate
v ( t ) L 2 ( Ω ) 2 + ν Ω t ( | v r | 2 + | v φ | 2 + | v z | 2 ) d x d t + ν Ω t v r 2 r 2 + v φ 2 r 2 d x d t 3 f L 2 , 1 ( Ω t ) 2 + 2 v ( 0 ) L 2 ( Ω ) 2 D 1 2 .
Proof. 
Multiplying ( 7 ) 1 by v r , ( 7 ) 2 by v φ , ( 7 ) 3 by v z , adding the results, and integrating over Ω yield
1 2 d d t Ω ( v r 2 + v φ 2 + v z 2 ) d x + ν Ω ( | v r | 2 + | v φ | 2 + | v z | 2 ) d x + ν Ω v r 2 r 2 + v φ 2 r 2 d x + Ω ( p , r v r + p , z v z ) d x = Ω ( f r v r + f φ v φ + f z v z ) d x .
The last term on the l.h.s. of (47) vanishes in virtue of the equation of continuity ( 7 ) 4 and the boundary conditions.
Using v 2 = v r 2 + v φ 2 + v z 2 , (47) takes the form
1 2 d d t v L 2 ( Ω ) 2 + ν Ω ( | v r | 2 + | v φ | 2 + | v z | 2 ) d x + ν Ω v r 2 r 2 + v φ 2 r 2 d x = Ω ( f r v r + f φ v φ + f z v z ) d x .
Applying the Hölder inequality to the r.h.s. of (48) yields
d d t v L 2 ( Ω ) f L 2 ( Ω ) ,
where f 2 = f r 2 + f φ 2 + f z 2 .
Integrating (49) with respect to time gives
v L 2 ( Ω ) f L 2 , 1 ( Ω t ) + v ( 0 ) L 2 ( Ω ) .
Integrating (48) with respect to time, using the Hölder inequality on the r.h.s. of (48), and exploiting (50), we obtain
1 2 v ( t ) L 2 ( Ω ) 2 + ν Ω t ( | v r | 2 + | v φ | 2 + | v z | 2 ) d x d t + ν Ω t v r 2 r 2 + v φ 2 r 2 d x d t f L 2 , 1 ( Ω t ) ( f L 2 , 1 ( Ω t ) + v ( 0 ) L 2 ( Ω ) ) + 1 2 v ( 0 ) L 2 ( Ω ) 2 .
The above inequality implies (46). This concludes the proof. □
Lemma 2.
Consider problem (12). Assume that f 0 L , 1 ( Ω t ) and u ( 0 ) L ( Ω ) . Then,
u ( t ) L ( Ω ) f 0 L , 1 ( Ω t ) + u ( 0 ) L ( Ω ) D 2 .
Proof. 
Multiplying ( 12 ) 1 by u | u | s 2 , s > 2 , integrating over Ω and by parts, we obtain
1 s d d t u L s ( Ω ) s + 4 ν ( s 1 ) s 2 | u | s / 2 L 2 ( Ω ) 2 + ν s Ω ( | u | s ) , r d r d z = Ω f 0 u | u | s 2 d x .
From [9], it follows that u | r = 0 = 0 . Moreover, using boundary conditions, (53) implies
d d t u L s ( Ω ) f 0 L s ( Ω ) .
Integrating (54) with respect to time and passing with s , we derive (52). This ends the proof. □
Lemma 3.
Let estimates (46) and (52) hold. Then,
v φ L 4 ( Ω t ) D 1 1 / 2 D 2 1 / 2 .
Proof. 
We have
Ω t | v φ | 4 d x d t = Ω t r 2 v φ 2 v φ 2 r 2 d x d t r v φ L ( Ω t ) 2 ( Ω t ) v φ 2 r 2 d x d t D 2 2 D 1 2 .
This implies (55) and concludes the proof. □
Lemma 4.
Consider problem (22). Assume that ω 1 L 6 / 5 ( Ω ) , where Ω = ( 0 , R ) × ( a , a ) . Then, there exists a weak solution to problem (22), such that ψ 1 H 1 ( Ω ) and the estimate
ψ 1 1 , Ω c | ω 1 | 6 / 5 , Ω
holds.
Proof. 
Multiplying ( 22 ) 1 by ψ 1 and using the boundary conditions, we obtain
ψ 1 1 , Ω 2 + a a ψ 1 2 | r = 0 d z = Ω ω 1 ψ 1 d x .
Applying the Hölder and Young inequalities to the r.h.s. implies (56). The Fredholm theorem gives existence. This ends the proof. □
Remark 2.
We have to emphasize that the weak solution ψ 1 of (22) does not vanish on the axis of symmetry. It also follows from [9].
From Lemma 2.4 in [8], we also have
Lemma 5.
Let f C ( ( 0 , R ) × ( a , a ) ) , f | r R = 0 . Let 1 < r 3 , 0 s r , s 2 , q r , r ( 3 s ) 3 r . Then, there exists a positive constant c = c ( s , r ) , such that
Ω | f | q r s d x 1 / q c | f | r , Ω 3 s q 3 r + 1 | f | r , Ω 3 r 3 s q ,
where f does not depend on φ.
Notation 1
(see [11]). First, we introduce the Fourier transform. Let f S ( R ) , where S ( R ) is the Schwartz space of all complex-valued rapidly decreasing and infinitely differentiable functions on R . Then, the Fourier transform of f and its inverse are defined by
f ^ ( λ ) = 1 2 π R e i λ τ f ( τ ) d τ , f ˇ ^ ( τ ) = f ( τ ) = 1 2 π R e i λ τ f ^ ( λ ) d λ
and f ^ ˇ = f ˇ ^ = f .
By H μ k ( R + ) , we denote a weighted space with the norm
u H μ k ( R + ) 2 = i = 0 k R + | r i u | 2 r 2 ( μ k + i ) r d r .
In view of transformation τ = ln r , r = e τ , d r = e τ d τ , we have the equivalence
i = 0 k R + | r i u | 2 r 2 ( μ k + i ) r d r i = 0 k R | τ i u | 2 e 2 h τ d τ
which holds for u ( τ ) = u ( ln r ) = u ( r ) , h = k 1 μ .
In view of the Fourier transform (58) and the Parseval identity, we have
+ i h + + i h j = 0 k | λ | 2 j | u ^ ( λ ) | 2 d λ = R j = 0 k | τ j u | 2 e 2 h τ d τ .

3. Estimates for the Stream Function ψ 1

In this section, we derive many estimates for ψ 1 = ψ / r , where ψ is the stream function, in terms of Γ 1 , 2 , Ω t + | Γ | 2 , , Ω t (recall that Γ = ω 1 ). Function ψ 1 was introduced by Thomas Hou in [12]. Lemma 6 is proved by applying the energy-type method. Inequalities (85) and (93) are proved by applying the technique of Kondratiev (see [10]) to problem (61). We need the inequalities to prove inequality (173), so they are very important.
To prove (85) and (93), we require that ψ 1 | r = 0 = 0 ; however, the theory developed by [9] does not imply the restriction.
Recall that ψ 1 is a solution to the problem
ψ 1 , r r ψ 1 , z z 3 r ψ 1 , r = ω 1 in Ω = ( 0 , R ) × ( a , a ) , ψ 1 | r = R = 0 , ψ 1 satisfies the periodic boundary conditions on S 2 .
Lemma 6.
For sufficiently regular solutions to (61), the following estimates hold
Ω ( ψ 1 , r r 2 + ψ 1 , r z 2 + ψ 1 , z z 2 ) d x + Ω 1 r 2 ψ 1 , r 2 d x + a a ψ 1 , z 2 | r = 0 d z + a a ψ 1 , r 2 | r = R d z c | ω 1 | 2 , Ω 2
and
Ω ( ψ 1 , z z r 2 + ψ 1 , z z z 2 ) d x + a a ψ 1 , z z 2 | r = 0 d z c | ω 1 , z | 2 , Ω 2
and
Ω ( ψ 1 , r r z 2 + ψ 1 , r z z 2 + ψ 1 , z z z 2 ) d x + a a ψ 1 , z z 2 | r = 0 d z + a a ψ 1 , r z 2 | r = R d z c | ω 1 , z | 2 , Ω 2 .
Proof. 
First, we prove (62). Multiplying ( 61 ) 1 by ψ 1 , z z and integrating over Ω yields
Ω ψ 1 , r r ψ 1 , z z d x Ω ψ 1 , z z 2 d x 3 Ω 1 r ψ 1 , r ψ 1 , z z d x = Ω ω 1 ψ 1 , z z d x .
Integrating by parts with respect to r in the first term implies
Ω ( ψ 1 , r ψ 1 , z z r ) , r d r d z + Ω ψ 1 , r ψ 1 , z z r d x + Ω ψ 1 , r ψ 1 , z z d r d z Ω ψ 1 , z z 2 d x 3 Ω ψ 1 , r ψ 1 , z z d r d z = Ω ω 1 ψ 1 , z z d x .
Continuing, we have
a a ψ 1 , r ψ 1 , z z r | r = 0 r = R d z + Ω ψ 1 , r ψ 1 , z z r d x Ω ψ 1 , z z 2 d x 2 Ω ψ 1 , r ψ 1 , z z d r d z = Ω ω 1 ψ 1 , z z d x .
The first integral in (66) vanishes because ψ 1 , r r | r = 0 = 0 , ψ 1 , z z | r = R = 0 . Integrating by parts with respect to z in the last term on the l.h.s. of (66), and using the periodic boundary conditions on S 2 , we obtain
Ω ψ 1 , r ψ 1 , z z r d x Ω ψ 1 , z z 2 d x + 2 Ω ψ 1 , r z ψ 1 , z d r d z = Ω ω 1 ψ 1 , z z d x .
Integrating by parts with respect to z in the first term in (67) and using the boundary conditions on S 2 , we have
Ω ( ψ 1 , z r 2 + ψ 1 , z z 2 ) d x Ω ( ψ 1 , z 2 ) , r d r d z = Ω ω 1 ψ 1 , z z d x ,
where the last term on the l.h.s. equals
a a ψ 1 , z 2 | r = 0 r = R d z = a a ψ 1 , z 2 | r = 0 d z
because ψ 1 , z | r = R = 0 . Using this in (68) and applying the Hölder and Young inequalities to the r.h.s. of (68) yield
Ω ( ψ 1 , r z 2 + ψ 1 , z z 2 ) d x + a a ψ 1 , z 2 | r = 0 d z c | ω 1 | 2 , Ω 2 .
We multiply ( 61 ) 1 by 1 r ψ 1 , r and integrate over Ω. Then, we have
3 Ω | 1 r ψ 1 , r | 2 d x = Ω ψ 1 , r r 1 r ψ 1 , r d x Ω ψ 1 , z z 1 r ψ 1 , r d x Ω ω 1 1 r ψ 1 , r d x .
The first term on the r.h.s. of (70) equals
1 2 Ω r ψ 1 , r 2 d r d z = 1 2 a a ψ 1 , r 2 | r = 0 r = R d z = 1 2 a a ψ 1 , r 2 | r = R d z ,
because ψ 1 , r | r = 0 = 0 (see [9]). Applying the Hölder and Young inequalities to the last two terms on the r.h.s. of (70) implies
Ω | 1 r ψ 1 , r | 2 d x + 1 2 a a ψ 1 , r 2 | r = R d z c ( | ψ 1 , z z | 2 , Ω 2 + | ω 1 | 2 , Ω 2 ) .
Inequalities (69) and (71) imply the estimate
Ω ( ψ 1 , r z 2 + ψ 1 , z z 2 ) d x + Ω | 1 r ψ 1 , r | 2 d x + a a ψ 1 , z 2 | r = 0 d z + a a ψ 1 , r 2 | r = R d z c | ω 1 | 1 , Ω 2 .
From ( 61 ) 1 , we have
| ψ 1 , r r | 2 , Ω 2 | ψ 1 , z z | 2 , Ω 2 + 3 | 1 r ψ 1 , r | 2 , Ω 2 + | ω 1 | 2 , Ω 2 .
Inequalities (72) and (73) imply (62).
Now, we show (63). We differentiate ( 61 ) 1 with respect to z, multiply by ψ 1 , z z z , and integrate over Ω. Then, we obtain
Ω ψ 1 , r r z ψ 1 , z z z d x + Ω ψ 1 , z z z 2 d x + 3 Ω 1 r ψ 1 , r z ψ 1 , z z z d x = Ω ω 1 , z ψ 1 , z z z d x .
Integrating by parts with respect to z yields
Ω ψ 1 , r r z ψ 1 , z z z d x = Ω ( ψ 1 , r r z ψ 1 , z z ) , z d x Ω ψ 1 , r r z z ψ 1 , z z d x ,
where the first integral vanishes in view of periodic boundary conditions on S 2 . Integrating by parts with respect to r in the second integral in (75) gives
Ω ( ψ 1 , r z z ψ 1 , z z r ) , r d r d z + Ω ψ 1 , r z z 2 d x + Ω ψ 1 , r z z ψ 1 , z z d r d z ,
where the first integral vanishes because
ψ 1 , r z z r | r = 0 = 0 , ψ 1 , z z | r = R = 0 .
In view of the above considerations, (74) takes the form
Ω ( ψ 1 , r z z 2 + ψ 1 , z z z 2 ) d x + Ω ψ 1 , r z z ψ 1 , z z d r d z + 3 Ω ψ 1 , r z ψ 1 , z z z d r d z = Ω ω 1 , z ψ 1 , z z z d x .
Integrating by parts with respect to z in the last term on the l.h.s. of (76) and using the periodic boundary conditions on S 2 , we have
Ω ( ψ 1 , r z z 2 + ψ 1 , z z z 2 ) d x Ω r ψ 1 , z z 2 d r d z = Ω ω 1 , z ψ 1 , z z z d x .
Applying the Hölder and Young inequalities to the r.h.s. of (77) yields
Ω ( ψ 1 , r z z 2 + ψ 1 , z z z 2 ) d x + a a ψ 1 , z z 2 | r = 0 d z c | ω 1 , z | 2 , Ω 2 ,
where we used that ψ 1 , z z | r = R = 0 .
The above inequality implies (63).
Finally, we show (64). We differentiate ( 61 ) 1 with respect to z, multiply by ψ 1 , r r z , and integrate over Ω. Then, we have
Ω ψ 1 , r r z 2 d x Ω ψ 1 , z z z ψ 1 , r r z d x 3 Ω 1 r ψ 1 , r z ψ 1 , r r z d x = Ω ω 1 , z ψ 1 , r r z d x .
Integrating by parts with respect to z in the second term in (78) implies
Ω ψ 1 , z z z ψ 1 , r r z d x = Ω ψ 1 , z z ψ 1 , r r z z d x = Ω ( ψ 1 , z z ψ 1 , r z z r ) r d r d z Ω ψ 1 , r z z 2 d x Ω ψ 1 , z z ψ 1 , r z z d r d z ,
where the first term vanishes because
ψ 1 , r z z r | r = 0 = 0 , ψ 1 , z z | r = R = 0 .
Then, (78) takes the form
Ω ( ψ 1 , r r z 2 + ψ 1 , r z z 2 ) d x + Ω ψ 1 , z z ψ 1 , r z z d r d z + 3 Ω ψ 1 , r z ψ 1 , r r z d r d z = Ω ω 1 , z ψ 1 , r r z d x .
The second term in (79) equals
1 2 a a ψ 1 , z z 2 | r = 0 r = R d z = 1 2 a a ψ 1 , z z 2 | r = 0 d z
because ψ 1 , z z | r = R = 0 , and the last term on the l.h.s. of (79) has the form
3 2 Ω r ψ 1 , r z 2 d r d z = 3 2 a a ψ 1 , r z 2 | r = 0 r = R d z = 3 2 a a ψ 1 , r z 2 | r = R d z
because ψ 1 , r z | r = 0 = 0 .
Using the above expressions in (79) implies the equality
Ω ( ψ 1 , r r z 2 + ψ 1 , r z z 2 ) d x 1 2 a a ψ 1 , z z 2 | r = 0 d z + 3 2 a a ψ 1 , r z 2 | r = R d z = Ω ω 1 , z ψ 1 , r r z d x .
Applying the Hölder and Young inequalities on the r.h.s. of (80) gives
Ω ( ψ 1 , r r z 2 + ψ 1 , r z z 2 ) d x 1 2 a a ψ 1 , z z 2 | r = 0 d z + 3 2 a a ψ 1 , r z 2 | r = R d z c | ω 1 , z | 2 , Ω 2 .
Inequalities (81) and (63) imply (64). This ends the proof. □
Lemma 7.
For sufficiently regular solutions to (61), the following inequality
| 1 r ψ 1 , r z | 2 , Ω c | ω 1 , z | 2 , Ω
holds.
Proof. 
Differentiating (61) with respect to z implies
ψ 1 , r r z ψ 1 , z z z 3 r ψ 1 , r z = ω 1 , z .
From (83), we have
| 1 r ψ 1 , r z | 2 , Ω c ( | ψ 1 , r r z | 2 , Ω + | ψ 1 , z z z | 2 , Ω + | ω 1 , z | 2 , Ω ) .
Using (64) in (84) yields (82). This concludes the proof. □
Now, we estimate | ψ 1 , z z r | 2 , Ω .
Lemma 8.
Let ψ 1 be a weak solution to problem (61), which vanishes on the axis of symmetry. Then, such sufficiently regular solutions to problem (61) satisfy the following estimate:
Ω ψ 1 , z z 2 r 2 d x + Ω ψ 1 , z r r 2 + ψ 1 , z r 2 r 2 + ψ 1 , z 2 r 4 d x c | ω 1 , z | 2 , Ω 2 .
Proof. 
Differentiating (61) with respect to z yields
Δ ψ 1 , z 2 r ψ 1 , z r = ω 1 , z , ψ 1 , z | S 1 = 0 , ψ 1 , z satisfies periodic boundary conditions on S 2 .
Applying Lemma 17 (also see Lemma 3.1 from [11]) to problem (86) gives
Ω ψ 1 , z r r 2 + ψ 1 , z r 2 r 2 + ψ 1 , z 2 r 4 d x c ( | ω 1 , z | 2 , Ω 2 + | ψ 1 , z z z | 2 , Ω 2 ) c | ω 1 , z | 2 , Ω 2 ,
where (63) is used in the last inequality.
To examine solutions to (86) we use the notation
u = ψ 1 , z .
Then, (86) takes the form
Δ u 2 r u , r = ω 1 , z , u | S 1 = 0 , u satisfies periodic boundary conditions on S 2 .
We multiply ( 88 ) 1 by u r 2 , integrate over Ω, and express the Laplacian operator in cylindrical coordinates. Then, we have
Ω u , r r + 1 r u , r + u , z z u r 2 d x 2 Ω 1 r u , r u r 2 d x = Ω ω 1 , z u r 2 d x .
Integrating by parts with respect to z in the third term under the first integral, we obtain
Ω u , z 2 r 2 d x = Ω u , r r + 3 r u , r u r 2 d x + Ω ω 1 , z u r 2 d x .
Applying the Hölder and Young inequalities to the r.h.s. integrals, using u = ψ 1 , z and (87), we derive
Ω ψ 1 , z z 2 r 2 d x c Ω ψ 1 , z r r 2 + ψ 1 , z r 2 r 2 + ψ 1 , z 2 r 4 d x + c | ω 1 , z | 2 , Ω 2 .
Using (87) in (92) implies (85). This concludes the proof. □
Remark 3.
Lemma 8 is necessary in the proof of global regular axially symmetric solutions to problem (6). However, it imposes strong restrictions on the solutions to (6) because the condition ψ 1 | r = 0 = 0 implies v z | r = 0 = 0 . We do not know how to omit this restriction in the proof presented in this paper.
Lemma 9.
Let μ > 0 and ω 1 H μ 1 ( Ω ) . Then, for sufficiently smooth solutions to (61), the following estimate is valid
Ω ψ 1 , r r r 2 + ψ 1 , r r 2 r 2 + ψ 1 , r 2 r 4 r 2 μ d x c R 2 μ ω 1 1 , Ω 2 .
Proof. 
To prove the lemma, we introduce a partition of unity { ζ ( i ) ( r ) } i = 1 , 2 , such that
i = 1 2 ζ ( i ) ( r ) = 1
and
ζ ( 1 ) ( r ) = 1 r r 0 , 0 r r 0 + λ , ζ ( 2 ) ( r ) = 0 r r 0 , 1 r r 0 + λ ,
where r 0 < R and ζ ( i ) ( r ) , i = 1 , 2 , are smooth functions.
We introduce the notation
ψ 1 ( i ) = ψ 1 ζ ( i ) , ω 1 ( i ) = ω 1 ζ ( i ) , i = 1 , 2 .
Then, function (94) satisfies the equations
ψ 1 , r r ( i ) ψ 1 , z z ( i ) 3 r ψ 1 , r ( i ) = 2 ψ 1 , r ζ ˙ ( i ) ψ 1 ζ ¨ ( i ) 3 r ψ 1 ζ ˙ ( i ) + ω 1 ( i ) g ( i ) , i = 1 , 2 ,
where the dot denotes the derivative with respect to r.
First, we consider the case i = 1 . Differentiating (95) for i = 1 with respect to r yields
ψ 1 , r r r ( 1 ) ψ 1 , r z z ( 1 ) 3 r ψ 1 , r r ( 1 ) + 3 r 2 ψ 1 , r ( 1 ) = g , r ( 1 ) .
We introduce the notation
v = ψ 1 , r ( 1 ) , f = g , r ( 1 ) .
Then, (96) takes the form
v , r r v , z z 3 r v , r + 3 r v = f in Ω r 0 , v | r = r 0 = 0 , v | S 2 satisfies periodic boundary conditions ,
where Ω r 0 = { x Ω : r ( 0 , r 0 ) , z ( a , a ) } and r 0 < R .
Multiplying ( 98 ) 1 by r 2 yields
r 2 v , r r 3 r v , r + 3 v = r 2 ( f + v , z z ) g ( r , z )
or equivalently
r r ( r r v ) 2 r r v + 3 v = g ( r , z ) .
We introduce the new variable
τ = ln r , r = e τ .
Since r r = τ , Equation (99) takes the form
τ 2 v + 2 τ v + 3 v = g ( e τ , z ) g ( τ , z ) .
Applying the Fourier transform (58) to (100) gives
λ 2 v ^ + 2 i λ v ^ + 3 v ^ = g ^ .
Looking for solutions to the algebraic equation
λ 2 + 2 i λ + 3 = 0
we see that it has two solutions
λ 1 = 3 i , λ 2 = i .
For λ { 3 i , i } , we can write solutions to (101) in the form
v ^ = 1 λ 2 + 2 i λ + 3 g ^ R ( λ ) g ^ .
Since R ( λ ) does not have poles on the line Im λ = 1 μ = h , μ ( 0 , 1 ) , we can use Lemma 3.1 from [11]. Then, we obtain
+ i h + i h j = 0 2 | λ | 2 ( 2 j ) | v ^ | 2 d λ c + i h + + i h j = 0 2 | λ | 2 ( 2 j ) | R ( λ ) g ^ | 2 d λ c + i h + + i h | g ^ | 2 d λ .
By applying the Parseval identity, inequality (103) becomes
R j = 0 2 | τ j v | 2 e 2 h τ d τ c R | g | 2 e 2 h τ d τ .
Passing to variable r yields
j = 0 2 R + | r j v | 2 r 2 ( μ + j 2 ) r d r c R + | g | 2 r 2 ( μ 2 ) r d r .
Using that g = r 2 ( f + v , z z ) , we have
j = 0 2 R + | r j v | 2 r 2 ( μ + j 2 ) r d r c R + | f + v , z z | 2 r 2 μ r d r .
Recalling notation (97), we derive from (104) the inequality
j = 0 2 Ω | r j ψ 1 , r ( 1 ) | 2 r 2 ( μ + j 2 ) d x c Ω | g , r ( 1 ) | 2 r 2 μ d x + c Ω | ψ 1 , r z z | 2 r 2 μ d x .
In view of (63),
| ψ 1 , r z z | 2 , Ω c | ω 1 , z | 2 , Ω .
The first term on the r.h.s. of (105) can be estimated by
| g , r ( 1 ) | 2 , μ , Ω c ( | ψ 1 , r r | 2 , Ω + | ψ 1 , r | 2 , Ω + | ψ 1 | 2 , Ω + | ω 1 , r | 2 , Ω + | ω 1 | 2 , Ω ) .
Lemma 6 and inequalities (105), (106), and (107) imply
Ω | ψ 1 , r r r ( 1 ) | 2 + | ψ , r r ( 1 ) | 2 r 2 + | ψ , r ( 1 ) | 2 r 4 r 2 μ r d r d z + Ω | ψ 1 , r z z | 2 d x c ( | ω 1 , r | 2 , Ω 2 + | ω 1 | 2 , Ω 2 ) .
Function ψ 1 ( 2 ) is a solution to the problem
Δ ψ 1 ( 2 ) = 2 ψ 1 , r ζ ˙ ( 2 ) ψ 1 ζ ¨ ( 2 ) + 2 r ψ 1 , r ( 2 ) 3 r ψ 1 ζ ˙ ( 2 ) + ω 1 ( 2 ) in Ω ¯ r 0 , ψ 1 ( 2 ) | r = R = 0 , ψ ( 2 ) = 0 for r r 0 , ψ ( 2 ) satisfies periodic boundary conditions on S 2 ,
where Ω ¯ r 0 = { x R 3 : r 0 r R , z ( a , a ) } and the dot denotes the derivative with respect to r.
For solutions to (109), the following estimate holds
ψ 1 ( 2 ) 3 , Ω c ( ψ 1 , r 1 , Ω + ψ 1 1 , Ω + ω 1 ( 2 ) 1 , Ω ) c ω 1 1 , Ω .
From (56), (108), and (110), inequality (93) follows. This ends the proof. □

4. Estimates for Φ and Γ

Let Ω = { ( r , z ) : r ( 0 , R ) , z ( a , a ) } . Let Φ = ω r / r , Γ = ω φ / r and Φ, Γ are solutions to problems (17)–(20).
Lemma 10.
Assume that Φ ( 0 ) , Γ ( 0 ) L 2 ( Ω ) , F ¯ r , F ¯ φ L 2 ( 0 , t ; L 6 / 5 ( Ω ) ) . Let D 2 be defined by (52) and let
I 3 = Ω t | v φ r Φ Γ | d x d t < .
Then
| Φ ( t ) | 2 , Ω 2 + | Γ ( t ) | 2 , Ω 2 + ν ( Φ 1 , 2 , Ω t 2 + Γ 1 , 2 , Ω t 2 ) ϕ ( D 2 ) | Ω t v φ r Φ Γ d x d t | + ϕ ( D 2 ) ( | F ¯ r | 6 / 5 , 2 , Ω t 2 + | F ¯ φ | 6 / 5 , 2 , Ω t 2 ) + | Φ ( 0 ) | 2 , Ω 2 + | Γ ( 0 ) | 2 , Ω 2 ϕ ( D 2 ) I 3 + D 8 .
Proof. 
Multiplying (17) by Φ and integrating over Ω yields
1 2 d d t | Φ | 2 , Ω 2 + | Φ | 2 , Ω 2 a a Φ | r = 0 r = R d z = Ω ( ω r r + ω z z ) v r r Φ d x + Ω F ¯ r Φ d x .
To derive the second term on the l.h.s. of (112), we consider (17) in
Ω ¯ = { x R 3 : r < R , z ( a , a ) , φ ( 0 , 2 π ) } .
Then, by applying the Green theorem and the boundary conditions, we obtain the second term on the l.h.s. of (112) on Ω ¯ . Considering that all quantities in (112) do not depend on φ, we can omit the integration with respect to φ and obtain (112).
Using Φ | r = R = 0 and (113), we have
1 2 d d t | Φ | 2 , Ω 2 + | Φ | 2 , Ω 2 Ω ( ω r r + ω z z ) v r r Φ d x + Ω F ¯ r Φ d x Ω v φ , z r v r r + r ( r v φ ) r z v r r Φ r d r d z + Ω F ¯ r Φ d x = Ω v φ ( z r v r r Φ + r v r r z Φ ) d x + Ω r r v φ z v r r Φ d r d z Ω v φ z r v r r Φ + z v r r r Φ d x + Ω F ¯ r Φ d x = a a r v φ z v r r Φ | r = 0 r = R d z + Ω v φ r v r r z Φ z v r r r Φ d x + Ω F ¯ r Φ d x a a r v φ z v r r Φ | r = 0 r = R d z + I + Ω F ¯ r Φ d x ,
Using the periodic boundary conditions on S 2 , the boundary term vanishes because v φ | r = R = 0 , v r | r = R = 0 , Φ | r = R = 0 , and
a a r v φ z v r r Φ | r = 0 d z = 0
because [9] implies the following expansions near the axis of symmetry
v φ = a 1 ( z , t ) r + a 2 ( z , t ) r 3 + , v r = a ¯ 1 ( z , t ) r + a ¯ 2 ( z , t ) r 3 +
and Φ = v φ , z r .
Finally, I I 1 + I 2 , where
I 1 Ω | v φ r v r r Φ , z | d x , I 2 Ω | v φ z v r r Φ , r | d x .
Now, we estimate I 1 and I 2 . Recall that v r r = ψ 1 , z . Then,
I 1 Ω | v φ ψ 1 , r z Φ , z | d x = Ω | r v φ ψ 1 , r z r Φ , z | d x | r v φ | , Ω | ψ 1 , r z r | 2 , Ω | Φ , z | 2 , Ω I 1 1 .
From (52) and (82), we have (recall that Γ = ω 1 ):
I 1 1 c D 2 | Γ , z | 2 , Ω | Φ , z | 2 , Ω .
Similarly, we calculate
I 2 Ω | v φ ψ 1 , z z Φ , r | d x | r v φ | , Ω | ψ 1 , z z r | 2 , Ω | Φ , r | 2 , Ω c D 2 | Γ , z | 2 , Ω | Φ , r | 2 , Ω ,
where (85) is used.
Finally, the last term on the r.h.s. of (113) is bounded by
ε | Φ | 6 , Ω 2 + c ( 1 / ε ) | F ¯ r | 6 / 5 , Ω 2 .
Using estimates (115)–(117) in (113), assuming that ε is sufficiently small and applying the Poincaré inequality, we obtain
d d t | Φ | 2 , Ω 2 + Φ 1 , Ω 2 c D 2 | Γ , z | 2 , Ω | Φ | 2 , Ω + c | F ¯ r | 6 / 5 , Ω 2 .
Multiplying (18) by Γ, integrating over Ω, using the boundary conditions and explanations about applying the Green theorem appeared below (112), we obtain
1 2 d d t | Γ | 2 , Ω 2 + | Γ | 2 , Ω 2 a a Γ 2 | r = 0 r = R d z 2 | Ω v φ r Φ Γ d x | + Ω F ¯ φ Γ d x .
Using Γ | r = R = 0 , applying the Hölder and Young inequalities to the last term on the r.h.s. of (119) and using the Poincaré inequality, we derive
d d t | Γ | 2 , Ω 2 + Γ 1 , Ω 2 2 Ω v φ r Φ Γ d x + c | F ¯ φ | 6 / 5 , Ω 2 .
From (118) and (120), we have
d d t ( | Φ | 2 , Ω 2 + | Γ | 2 , Ω 2 ) + Φ 1 , Ω 2 + Γ 1 , Ω 2 ϕ ( D 2 ) | Ω v φ r Φ Γ d x d t | + ϕ ( D 2 ) ( | F ¯ r | 6 / 5 , Ω 2 + | F ¯ φ | 6 / 5 , Ω 2 ) ,
where ϕ is an increasing positive function. Integrating (121) with respect to time yields (111). This ends the proof. □
Lemma 11.
Let the assumptions of Lemma 16 hold.
Let v φ L ( 0 , t ; L d ( Ω ) ) , d > 3 . Let θ = 1 3 d ε 1 3 d ε 2 > 0 , ε = ε 1 + ε 2 . Let ε 0 > 0 be arbitrarily small.
Then,
I 3 c | v φ | d , , Ω t ε [ c 1 ( 1 + | v φ | , Ω t 1 2 θ ε 0 ) Γ 1 , 2 , Ω t 1 2 θ + c 2 ] | Φ | 2 , Ω t 1 θ | Γ | 2 , Ω t ,
where c 1 and c 2 , depending on D 5 , D 6 , and D 7 , are introduced in L 1 4 below.
Proof. 
We examine
I 3 = Ω t | r v φ Φ r Γ r | d x d t Ω t | r v φ | 1 ε | v φ | ε | Φ r 1 ε 1 | | Γ r 1 ε 2 | d x d t = I 3 1 ,
where ε = ε 1 + ε 2 and ε i , for i = 1 , 2 , are positive numbers.
Using (52) and applying the Hölder inequality in I 3 1 yields
I 3 1 D 2 1 ε Ω t | v φ | 2 ε | Φ r 1 ε 1 | 2 d x d t 1 / 2 | Γ r 1 ε 2 | 2 , Ω t D 2 1 ε L | Γ / r 1 ε 2 | 2 , Ω t I 3 2 .
By the Hardy inequality, we obtain
Γ r L 2 , ε 2 ( Ω t ) c Γ L 2 , ε 2 ( Ω t ) c R ε 2 | Γ | 2 , Ω t .
Now, we estimate L,
L = 0 t Ω | v φ | 2 ε | Φ r 1 ε 1 | 2 d x d t 1 / 2 0 t | v φ | 2 ε σ , Ω 2 ε Ω | Φ r 1 ε 1 | q d x 2 / q d t 1 / 2 L 1 ,
where 1 / σ + 1 / σ = 1 , q = 2 σ . Let d = 2 ε σ . Then,
σ = d d 2 ε so q = 2 d d 2 ε .
Continuing,
L 1 sup t | v φ | d , Ω ε 0 t | Φ r 1 ε 1 | q , Ω 2 d t 1 / 2 L 1 1 L 1 2 .
Now, we estimate the second factor L 1 2 .
For this purpose, we use Lemma 5 for r = 2 . Let s q = 1 ε 1 . Then, q [ 2 , 2 ( 3 s ) ] . Since s = ( 1 ε 1 ) q we have the restriction 2 q 6 2 s = 6 2 ( 1 ε 1 ) q . Then,
2 q 6 3 2 ε 1
and 6 3 2 ε 1 > 2 for any ε 1 ( 0 , 1 ) .
Hence, Lemma 5 implies
L 1 2 = 0 t | Φ r 1 ε 1 | q , Ω 2 d t 1 / 2 c 0 t | Φ | 2 , Ω 2 ( 3 s q 1 2 ) | Φ | 2 , Ω 2 ( 3 2 3 s q ) d t 1 / 2 c | Φ | 2 , Ω t 3 s q 1 2 | Φ | 2 , Ω t 3 2 3 s q L 1 3 ,
where we used that for θ = 3 s q 1 2 , 1 θ = 3 2 3 s q so the Hölder inequality can be applied.
Using (173) in L 1 3 , we have
L 1 3 c ( D 5 1 2 θ | Γ | 2 , Ω t 1 2 θ + D 6 1 2 θ | v φ | , Ω t 1 2 θ ε 0 Γ 1 , 2 , Ω t 1 2 θ + D 7 1 2 θ ) · | Φ | 2 , Ω t 1 θ [ c 1 ( 1 + | v φ | , Ω t 1 2 θ ε 0 ) Γ 1 , 2 , Ω t 1 2 θ + c 2 ] | Φ | 2 , Ω t 1 θ L 1 4 ,
where c 1 and c 2 depend on D 5 , D 6 , and D 7 .
To justify the above inequality, we have to know that the following inequalities hold:
θ = 3 s q 1 2 > 0
and
1 θ = 3 2 3 s q > 0 .
Consider (125). Using the form of q and s q , we have
3 q s q 1 2 > 0 so 3 ( d 2 ε ) 2 d ( 1 ε 1 ) 1 2 > 0 .
Hence
3 2 3 d ε 1 + ε 1 1 2 > 0 so ε 1 3 d ( ε 1 + ε 2 ) > 0 .
Therefore, the following inequality
1 3 d ε 1 3 d ε 2 > 0
holds for d > 3 , and ε 2 is sufficiently small. Moreover, (127) implies
ε 1 > 3 d d d 3 ε 2 = 3 d 3 ε 2 .
To examine (126), we calculate
3 2 3 ( d 2 ε ) 2 d + 1 ε 1 = 1 + 3 d ε ε 1 = 1 1 3 d ε 1 + 3 d ε 2 .
Since (129) must be positive, we have the restriction
1 + 3 d ε 2 > 1 3 d ε 1
Using (128) in (130) implies
1 + 3 d ε 2 > 3 d ε 2
so there is no contradiction.
Hence, we have
θ = 1 3 d ε 1 3 d ε 2 , 1 θ = 1 1 3 d ε 1 + 3 d ε 2 ,
where d > 3 .
Finally,
I 3 c | v φ | d , , Ω t ε [ c 1 ( 1 + | v φ | , Ω t 1 2 θ ε 0 ) Γ 1 , 2 , Ω t 1 2 θ + c 2 ] | Φ | 2 , Ω t 1 θ · | Γ | 2 , Ω t .
This implies (122) and ends the proof. □
We introduce the quantity
X ( t ) = Φ V ( Ω t ) + Γ V ( Ω t ) .
Lemma 12.
Let the assumptions of Lemmas 10 and 11 hold. Let θ = 1 3 d ε 1 3 d ε 2 , ε = ε 1 + ε 2 . Then
X 2 c 0 | v φ | d , , Ω t 4 ε θ ( 1 + | v φ | , Ω t 2 ε 0 ) + c 0 | v φ | d , , Ω t 2 ε θ + D 8 2 ,
where c 0 = ϕ ( D 5 , D 6 , D 7 ) .
Proof. 
In view of notation (132), inequalities (111) and (122) imply
X 2 c | v φ | d , , Ω t ε [ c 1 ( 1 + | v φ | , Ω t 1 2 θ ε 0 ) X 1 1 2 θ + c 2 X 1 θ ] X + D 8 α 1 X 2 1 2 θ + α 2 X 2 θ + D 8 2 .
Applying the Young inequality in (134) implies
X 2 c α 1 4 θ + c α 2 2 θ + D 8 2 .
This yields (133) and concludes the proof. □
Remark 4.
Consider exponents in (133). Then,
δ = 4 ε θ = 4 ε ( 1 3 d ) ε 1 3 d ε 2 , δ 0 = 2 ε ( 1 3 d ) ε 1 3 d ε 2 .
For ε 2 small, we have
δ = 4 1 3 d + ε , δ 0 = 2 1 3 d + ε 0 ,
where ε and ε 0 are positive numbers that can be chosen to be very small.
For d = 12 , it follows that
δ = 16 3 + ε , δ 0 = 8 3 + ε 0 .
This ends the remark.
Lemma 13.
Assume that ε 1 > a ε 2 , s > 1 , a = 16 + 6 3 · 6 16 , b = 2 · 6 ( 3 ε 1 ε 2 ) ( 6 · 3 16 ) ( ε 1 a ε 2 ) , and we choose 6 to be arbitrarily close to 6, and
D 9 s ( s ) = s 2 | f φ | 3 s 2 s + 1 , s , Ω t s + | v φ ( 0 ) | s , Ω s < .
Then, excluding cases where either v φ = 0 or v φ is small, we have
| v φ | 12 , , Ω t 6 c | v φ | , Ω t b ε 0 + ϕ ( D 5 , D 6 , D 7 ) + c ( D 8 + D 9 12 )
Proof. 
We multiply ( 7 ) 2 by v φ | v φ | s 2 , integrate over Ω, and exploit the relation v r r = ψ 1 , z . Then, we obtain
1 s d d t | v φ | s , Ω s + 4 ν ( s 1 ) s 2 | | v φ | s / 2 | 2 , Ω 2 = Ω ψ 1 , z | v φ | s d x + Ω f φ v φ | v φ | s 2 d x .
Integrating by parts in the first term on the r.h.s. of (138) and applying the Hölder and Young inequalities yield
| Ω ψ 1 , z | v φ | s d x | ε | z | v φ | s / 2 | 2 , Ω 2 + c ( 1 / ε ) Ω ψ 1 2 | v φ | s d x .
By the Poincaré inequality,
| | v φ | s / 2 | 2 , Ω 2 c | v φ | 3 s , Ω s
we can estimate the second term on the r.h.s. of (138) by
| f φ | 3 s 2 s + 1 , Ω | v φ | 3 s , Ω s 1 ε 1 | v φ | 3 s , Ω s + c ( 1 / ε 1 ) | f φ | 3 s 2 s + 1 , Ω s .
Using the above estimates with sufficiently small ε, ε 1 in (138) we derive the inequality
1 s d d t | v φ | s , Ω s + 1 s | | v φ | s / 2 | 2 , Ω 2 + 1 s | v φ | 3 s , Ω s c s Ω ψ 1 2 | v φ | s d x + c s | f φ | 3 s 2 s + 1 , Ω s .
In view of Lemma 2, the first term on the r.h.s. of (139) is bounded by
c s | u | , Ω t 6 Ω ψ 1 2 r 6 | v φ | s 6 d x c s D 2 6 | v φ | , Ω s 6 Ω ψ 1 2 r 6 d x ,
where 6 < 6 but 6 may be assumed arbitrarily close to 6.
Using the estimate in (139) yields
1 s d d t | v φ | s , Ω s c s D 2 6 | v φ | , Ω s 6 Ω ψ 1 2 r 6 d x + c s | f φ | 3 s 2 s + 1 , Ω s .
Integrating (140) with respect to time and using Lemma 18, we obtain
| v φ | s , Ω s c 1 s 2 D 2 6 | v φ | , Ω t s 6 Γ 1 , 2 , Ω t 2 + c s 2 | f φ | 3 s 2 s + 1 , s , Ω t s + | v φ ( 0 ) | s , Ω s c 1 s 2 D 2 6 | v φ | , Ω t s 6 Γ 1 , 2 , Ω t 2 + c D 9 s ( s ) ,
c 1 = c R 2 μ 1 + 4 ( 6 6 ) 2 .
Dividing (141) by | v φ | , Ω t s 6 implies
| | v φ | s , , Ω t | v φ | , Ω t | s 6 | v φ | s , Ω 6 c 1 s 2 D 2 6 Γ 1 , 2 , Ω t 2 + c | v φ | , Ω t s 6 D 9 s ( s ) .
The dividing by | v φ | , Ω t is justified because the following two cases are excluded from this paper:
Case 1: In the case, v φ = 0 the existence of global regular solutions to problem (6) is proven in [1,2,13].
Case 2: The existence of global regular solutions to problem (6) for v φ sufficiently small is proven in Appendix A.
Since Cases 1 and 2 are not considered in this paper, we can show the existence of positive constants c 0 and c 1 , such that
| | v φ | s , , Ω t | v φ | , Ω t | s 6 c ¯ 0
and
1 | v φ | , Ω t s 6 c ¯ 1 .
In view of (143) and (144), inequality (142) takes the form
c ¯ 0 | v φ | s , , Ω t 6 c 1 s 2 Γ 1 , 2 , Ω t 2 + c c ¯ 1 D 9 s ( s ) .
Let d = 12 . Then, θ = 1 4 ( 3 ε 1 ε 2 ) and (133) for d = 12 takes the form
X 2 c 0 | v φ | 12 , , Ω t 16 ε 3 ε 1 ε 2 ( 1 + | v φ | , Ω t 2 ε 0 ) + c 0 | v φ | 12 , , Ω t 8 ε 3 ε 1 ε 2 + D 8 .
Taking (145) for s = 12 and using (146) yield
| v φ | 12 , , Ω t 6 c 2 | v φ | 12 , , Ω t 16 ε 3 ε 1 ε 2 ( 1 + | v φ | , Ω t 2 ε 0 ) + c 2 | v φ | 12 , , Ω t 8 ε 3 ε 1 ε 2 + c D 8 + c D 9 12 ,
where C 2 = 144 c 1 c 0 c ¯ 0 .
To derive any estimate from (147), we need
16 ε 3 ε 1 ε 2 < 6
We see that (148) holds for
ε 1 > 16 + 6 3 · 6 16 ε 2 a ε 2 ,
where a > 11 .
In view of the Young inequality, (147) implies
| v φ | 12 , , Ω t 6 c | v φ | , Ω t b ε 0 + c + c ( D 8 + D 9 12 ) ,
where b = 2 · 6 ( 3 ε 1 ε 2 ) ( 6 · 3 16 ) ( ε 1 a ε 2 ) . The above inequality implies (137) and concludes the proof. □
Remark 5.
Exploiting (150) in (146) implies the inequality
X 2 c ( 1 + | v φ | , Ω t 2 ε 0 ) | v φ | , Ω t d ε 0 + ϕ ( D 5 , D 7 , D 8 , D 9 ) ,
where d = 16 b ε 3 ε 1 ε 2 and X is introduced in (132).
To prove Theorem 1, we need an estimate for | v φ | , Ω t . For this purpose, we need the result
Lemma 14.
Assume that quantities D 2 , D 5 , D 7 , D 8 , and D 9 are bounded. Assume that f φ / r L 1 ( 0 , t ; L ( Ω ) ) , v φ ( 0 ) L ( Ω ) .
Then, there exists an increasing positive function ϕ, such that
| v φ | , Ω t ϕ ( D 2 , D 5 , D 7 , D 8 , D 9 , f φ / r L 1 ( 0 , t ; L ( Ω ) ) , | v φ ( 0 ) | , Ω ) .
Proof. 
Recall Equation ( 7 ) 2 for v φ
v φ , t + v · v φ ν Δ v φ 1 r 2 v φ = ψ 1 , z v φ + f φ ,
where v r r = ψ 1 , z .
Multiplying (153) by v φ | v φ | s 2 and integrating over Ω yields
1 s d d t | v φ | s , Ω s + 4 ν ( s 1 ) s 2 | | v φ | s / 2 | 2 , Ω 2 + ν Ω | v φ | s r 2 d x = Ω ψ 1 , z v φ 2 | v φ | s 2 d x + Ω f φ v φ | v φ | s 2 d x .
The first term on the r.h.s. of (154) is bounded by
Ω | ψ , z | | v φ | s / 2 | v φ | s / 2 r d x ε Ω | v φ | s r 2 d x + c ( 1 / ε ) Ω ψ , z 2 | v φ | s d x ,
where the second integral is bounded by
| r v φ | , Ω 2 Ω | ψ 1 , z | 2 | v φ | s 2 d x D 2 2 | ψ 1 , z | s , Ω 2 | v φ | s , Ω s 2 .
The second term on the r.h.s. of (154) is estimated by
Ω | f φ | | v φ | s 1 d x = Ω | f φ r | r | v φ | s 1 d x | r v φ | , Ω Ω | f φ r | | v φ | s 2 d x D 2 | f φ r | s / 2 , Ω | v φ | s , Ω s 2 .
Using the above estimates in (154) and assuming that ε is sufficiently small, we obtain the inequality
1 s d d t | v φ | s , Ω s D 2 2 | ψ 1 , z | s , Ω 2 | v φ | s , Ω s 2 + | f φ r | s / 2 , Ω | v φ | s , Ω s 2 .
Simplifying, we have
1 2 d d t | v φ | s , Ω 2 D 2 2 | ψ 1 , z | s , Ω 2 + | f φ r | s / 2 , Ω .
Integrating with respect to time and passing with s , we derive
| v φ ( t ) | , Ω 2 D 2 2 0 t | ψ 1 , z | , Ω 2 d t + 0 t | f φ r | , Ω d t + | v φ ( 0 ) | , Ω 2 .
Since 0 t | ψ 1 , z | , Ω 2 d t X 2 , we can apply (151). Then (155) takes the form
| v φ | , Ω t 2 D 2 2 ( 1 + | v φ | , Ω t 2 ε 0 ) | v φ | , Ω t 96 ε ε 1 11 ε 2 ε 0 + D 2 ϕ ( D 5 , D 7 , D 8 , D 9 ) + D 2 2 0 t | f φ r | , Ω d t + | v φ ( 0 ) | , Ω 2 .
Hence, for ε 0 sufficiently small, we derive (152). This ends the proof. □
Remark 6.
Inequalities (151) and (152) imply
X ϕ ( D 2 , D 5 , D 7 , D 8 , D 9 , | f φ / r | , 1 , Ω t , | v φ ( 0 ) | , Ω )
The above inequality proves Theorem 1.

5. Estimates for the Swirl

Applying the energy method and using the estimate for the weak solution (see Lemma 1) and L -estimate for swirl (see Lemma 2), we derive the estimate
u L ( 0 , t ; H 1 ( Ω ) ) + u L 2 ( 0 , T ; H 2 ( Ω ) ) ϕ ( data ) .
This is a new result. It is necessary in the proof of (173).
In this section, we find estimates for solutions to the problem
u , t + v · u ν Δ u + 2 ν u , r r = r f φ f 0 in Ω t , u | S 1 = 0 , in Ω t , u | S 2 periodic boundary conditions , u | t = 0 = u ( 0 ) in Ω .
Lemma 15.
Assume that D 1 and D 2 are described by (46) and (52), respectively. Let u , z ( 0 ) , u , r ( 0 ) L 2 ( Ω ) , f 0 L 2 ( Ω t ) .
Then, solutions to (158) satisfy the estimates
| u , z ( t ) | 2 , Ω 2 + ν | u , z | 2 , Ω t 2 c ( D 1 2 D 2 2 + | u , z ( 0 ) | 2 , Ω 2 + | f 0 | 2 , Ω t 2 ) c D 3 2 ,
| u , r ( t ) | 2 , Ω 2 + ν ( | u , r r | 2 , Ω t 2 + | u , r z | 2 , Ω t 2 ) c D 1 2 ( 1 + D 2 2 ) + | u , r ( 0 ) | 2 , Ω 2 + | f 0 | 2 , Ω t 2 + | f 0 | 4 / 3 , 2 , S 1 t 2 c D 4 2 .
Proof. 
We differentiate (158) with respect to z, multiply by u , z , and integrate over Ω. To apply the Green theorem, we have to consider problem (158) in domain Ω ¯ = { x R 3 : r < R , z ( a , a ) , φ ( 0 , 2 π ) } . Then, we obtain
1 2 d d t | u , z | 2 , Ω ¯ 2 ν Ω ¯ div ( u , z u , z ) d x ¯ + ν Ω ¯ | u , z | 2 d x ¯ = Ω ¯ v , z · u · u , z d x ¯ + Ω ¯ f 0 , z u , z d x ¯ ,
where d x ¯ = d x d φ . The second term on the l.h.s. implies a boundary term, which vanishes due to boundary conditions. Since all functions in (161) do not depend on φ, any integral with respect to φ can be omitted.
Integrating by parts with respect to z in the term from the r.h.s. of (161) and using the boundary conditions on S 2 , we obtain
1 2 d d t | u , z | 2 , Ω 2 + ν | u , z | 2 , Ω 2 ν 2 a a u , z 2 | r = 0 r = R d z = Ω v , z z · u u d x + Ω v , z · u , z · u d x Ω f 0 u , z z d x .
The last term on the l.h.s. of (162) vanishes and the first term on the r.h.s. equals
1 2 Ω v , z z · u 2 d x = 1 2 S 1 v , z z · n ¯ u 2 d S 1 = 0 .
Applying the Hölder and Young inequalities to the other terms from the r.h.s. of (162) yields
d d t | u , z | 2 , Ω 2 + ν | u , z | 2 , Ω 2 c | u | , Ω 2 | v , z | 2 , Ω 2 + c | f 0 | 2 , Ω 2 .
Integrating (163) with respect to time gives
| u , z ( t ) | 2 , Ω 2 + ν | u , z | 2 , Ω 2 c | u | , Ω t 2 | v , z | 2 , Ω t 2 + | u , z ( 0 ) | 2 , Ω 2 + c | f 0 | 2 , Ω t 2 c D 1 2 D 2 2 + | u , z ( 0 ) | 2 , Ω 2 + c | f 0 | 2 , Ω t 2 .
The above inequality implies (159).
Differentiating (158) with respect to r gives
u , r t + v · u , r + v , r · u ν ( Δ u ) , r + 2 ν r u , r r 2 ν r 2 u , r = f 0 , r .
Multiplying (165) by u , r and integrating over Ω yields
1 2 d d t | u , r | 2 , Ω 2 + Ω v , r · u u , r d x ν Ω ( Δ u ) , r u , r d x + 2 ν Ω 1 r u , r r u , r d x 2 ν Ω u , r 2 r 2 d x = Ω f 0 , r u , r d x .
Now, we examine the particular terms in (166). The second term equals
Ω v , r · u u , r r d r d z = Ω ( v r , r r u + v z , r z u ) u , r r d r d z = Ω ( r v r , r u , r + r v z , r u , z ) u , r d r d z = Ω ( r v r , r u , r u , r + r v z , r u , r u , z ) d r d z = Ω [ ( r v r , r u , r ) , r + ( r v z , r u , r ) , z ] u d r d z I ,
where we use r v r , r u , r u | r = 0 = 0 (see [9]). Continuing, we write I in the following form:
I = Ω [ ( r v r , r ) , r + ( r v z , r ) , z ] u , r u d r d z Ω [ r v r , r u , r r + r v z , r u , r z ] u d r d z I 1 + I 2 .
To estimate I 1 , we calculate
I 1 1 = ( r v r , r ) , r + ( r v z , r ) , z = r v r , r r + v r , r + r v z , r z .
Since v = v r e ¯ r + v z e ¯ z is divergence free, we have
v r , r + v z , z + v r r = 0 .
Since Equation (167) is satisfied identically in Ω, we can differentiate (167) with respect to r. Then, we have
v r , r r + v z , z r + v r , r r v r r 2 = 0 .
Hence,
I 1 1 = v r r .
Then, I 1 equals
I 1 = Ω v r r u , r u d r d z .
Therefore,
| 0 t I 1 d t | | v r r | 2 , Ω t | u , r r | 2 , Ω t | u | , Ω t .
Next,
| I 2 | ε ( | u , r r | 2 , Ω 2 + | u , r z | 2 , Ω 2 ) + c ( 1 / ε ) | u | , Ω 2 ( | v r , r | 2 , Ω 2 + | v z , r | 2 , Ω 2 ) .
The third integral in (166) equals
J = ν Ω ( Δ u ) , r u , r d x = ν Ω u , r r r + 1 r u , r , r + u , r z z u , r r d r d z = ν Ω u , r r + 1 r u , r u , r r , r d r d z + ν Ω u , r r ( u , r r ) , r d r d z + ν Ω 1 r u , r ( u , r r ) , r d r d z + Ω u , r z 2 d x = ν a a u , r r + 1 r u , r u , r r | r = 0 r = R d z + ν Ω ( u , r r 2 + u , r z 2 ) d x + ν Ω u , r 2 r 2 d x + 2 ν Ω u , r r u , r d r d z ,
where the last term equals
ν Ω ( u , r 2 ) , r d r d z = ν a a u , r 2 | r = 0 r = R d z = ν a a u , r 2 | r = R d z
because u , r | r = 0 = ( v φ + v φ , r r ) | r = 0 = 0 .
To examine the boundary term in J, we recall the expansion of v φ near the axis of symmetry (see [9])
v φ = a 1 ( z , t ) r + a 2 ( z , t ) r 3 + ,
so
u = a 1 ( z , t ) r 2 + a 2 ( z , t ) r 4 +
Then, u , r r + 1 r u , r u , r r | r = 0 = 0 , and we have to emphasize that all calculations in this paper are performed for sufficiently regular solutions.
Therefore, the boundary term in J equals
J 1 = ν a a u , r r + 1 r u , r u , r r | r = R d z .
Projecting ( 158 ) 1 on S 1 yields
ν u , r r + 1 r u , r + 2 ν u , r r = f 0 on S 1 .
Hence,
u , r r | S 1 = u , r r 1 ν f 0 | S 1 .
Using the expression in J 1 gives
J 1 = 2 ν a a u , r 2 | r = R d z + a a f 0 u , r r | r = R d z .
The fourth term in (166) equals (169).
Using the above estimates and expressions in (166) yields
1 2 d d t | u , r | 2 , Ω 2 + ν Ω ( u , r r 2 + u , r z 2 ) d x + ν Ω u , r 2 r 2 d x 2 ν Ω u , r 2 r 2 d x Ω | v r r u , r u | d r d z + ε ( | u , r r | 2 , Ω 2 + | u , r z | 2 , Ω 2 ) + c ( 1 / ε ) | u | , Ω 2 ( | v r , r | 2 , Ω 2 + | v z , r | 2 , Ω 2 ) + c ( 1 / ε ) | f 0 | 2 , Ω 2 + | a a f 0 u , r r | r = R d z .
Integrating (170) with respect to time and assuming that ε is sufficiently small, we obtain
| u , r ( t ) | 2 , Ω 2 + ν ( | u , r r | 2 , Ω t 2 + | u , r z | 2 , Ω t 2 ) ν | u , r r | 2 , Ω t 2 + c | v r r | 2 , Ω t | u , r r | 2 , Ω t | u | , Ω t + c | u | , Ω t 2 ( | v r , r | 2 , Ω t 2 + | v z , r | 2 , Ω t 2 ) + c | f 0 | 2 , Ω t 2 + | u , r ( 0 ) | 2 , Ω 2 + ν 0 t a a u , r 2 | r = R d x d t + | 0 t a a f 0 u , r r | r = R d x d t | .
Using
Ω t | u , r r | 2 d x d t Ω t | v φ , r | 2 + v φ 2 r 2 d x d t c D 1 2
and
0 t a a u , r 2 | r = R d x d t ε | u , r | 2 , Ω t 2 + c ( 1 / ε ) | u , r | 2 , Ω t 2 ,
| 0 t a a f 0 u , r | r = R d x d t ε 1 | u , r | 4 , 2 , S 1 t 2 + c ( 1 / ε 1 ) | f 0 | 4 / 3 , 2 , S 1 t 2 ε 1 ( | u , r r | 2 , Ω t 2 + | u , r z | 2 , Ω t 2 ) + c ( 1 / ε 1 ) | f 0 | 4 / 3 , 2 , S 1 t 2
and Lemmas 1, 2, we have
| u , r ( t ) | 2 , Ω 2 + ν ( | u , r r | 2 , Ω t 2 + | u , r z | 2 , Ω t 2 ) c ( D 1 2 + D 1 2 D 2 + D 1 2 D 2 2 ) + c | f 0 | 2 , Ω t 2 + c | f 0 | 4 / 3 , 2 , S 1 t 2 + | u , r ( 0 ) | 2 , Ω 2 .
This inequality implies (160) and concludes the proof. □

6. Estimates for ω r , ω z

Inequality (173) is the most important inequality in this paper. To prove it, we need results from Section 3 and Section 5 as well as Lemma 2. By the energy method, we derive (174), where the first term on the r.h.s. is nonlinear. The aim of the proof of Lemma 16 is to show that
(∗)
J c | u | , Ω t ( D 1 + u L 2 ( 0 , t ; H 2 ( Ω ) ) ) Γ L 2 ( 0 , t ; H 1 ( Ω ) )
To show ( ) , we replace ω r , ω z in J by derivatives of u described by (13) and express components of velocities v r and v z by derivatives of ψ using (15). We perform such calculations in J, where we are able to extract the norm | u | , Ω t . Then, J becomes bilinear. Then, estimates (46), (159), and (160) imply ( ) by the Hölder inequality. Hence, J is bounded by the quantity, which is linear with respect to the norm of Γ. This implies that the main Theorem 1 can be proved.
Lemma 16.
Assume that D 5 = D 2 ( D 1 + D 3 + D 4 ) and D 6 = D 2 1 ε 0 D 3 , where D 1 , and D 2 are introduced in (46) and (52), and D 3 and D 4 are introduced in (159) and (160), respectively. Let
D 7 = | F r | 6 / 5 , 2 , Ω t 2 + | F z | 6 / 5 , 2 , Ω t 2 + | ω r ( 0 ) | 2 , Ω 2 + | ω z ( 0 ) | 2 , Ω 2 + | f φ | 2 , S 1 t ( D 3 + D 4 ) < .
Let ε 0 be an arbitrarily small positive number and let v φ L ( Ω t ) .
Let Γ L 2 ( 0 , t ; H 1 ( Ω ) ) .
Then,
ω r V ( Ω t ) 2 + ω z V ( Ω t ) 2 + | Φ | 2 , Ω t 2 c D 5 | Γ , z | 2 , Ω t + c D 6 | v φ | , Ω t ε 0 Γ 1 , 2 , Ω t + c D 7 .
Proof. 
Multiplying ( 9 ) 1 by ω r , ( 9 ) 3 by ω z , integrating over Ω t , and adding yield
1 2 ( | ω r ( t ) | 2 , Ω 2 + | ω z ( t ) | 2 , Ω 2 ) + ν ( | ω r | 2 , Ω t 2 + | ω z | 2 , Ω t 2 ) + ν | ω r r | 2 , Ω t 2 ν S 1 t n ¯ · ω r ω r d S 1 d t ν S 1 t n ¯ · ω z ω z d S 1 d t = Ω t [ v r , r ω r 2 + v z , z ω z 2 + ( v r , z + v z , r ) ω r ω z ] d x d t + Ω t ( F r ω r + F z ω z ) d x d t + 1 2 ( | ω r ( 0 ) | 2 , Ω 2 + | ω z ( 0 ) | 2 , Ω 2 ) J + Ω t ( F r ω r + F z ω z ) d x d t + 1 2 ( | ω r ( 0 ) | 2 , Ω 2 + | ω z ( 0 ) | 2 , Ω 2 ) .
Since ω r = v φ , z and v φ | r = R = 0 , we obtain
S 1 n ¯ · ω r ω r d S 1 = 0 .
Using ω z = v φ , r + v φ r , we derive
ν S 1 t n ¯ · ω z ω z d S 1 d t = ν 0 t a a r v φ , r + v φ r v φ , r + v φ r | r = R R d z d t I 1 .
Since v φ | r = R = 0 I 1 takes the form
I 1 = ν 0 t a a v φ , r r + v φ , r r v φ , r | r = R R d z d t .
Projecting ( 7 ) 2 on S 1 yields
ν v φ , r r + 1 r v φ , r = f φ on S 1 .
Hence,
I 1 = R 0 t a a f φ v φ , r | r = R d z d t = 0 t a a f φ u , r 1 R u | r = R d z d t .
Using (13) and (21) in J implies
J = Ω t [ 1 r 2 u , z 2 ( ψ 1 , z + r ψ 1 , r z ) + 1 r u , r 2 ( r ψ 1 , z r + 2 ψ 1 , z ) 1 r 2 u , r u , z ( r ψ 1 , z z + 3 ψ 1 , r + r ψ 1 , r r ) ] d x d t J 1 + J 2 + J 3 .
We integrate by parts in J 1 and use the boundary conditions on S 2 . Then, we have
J 1 = Ω t 1 r 2 u u , z z ( ψ 1 , z + r ψ 1 , r z ) d x d t + Ω t 1 r 2 u u , z ( ψ 1 , z z + r ψ 1 , r z z ) d x d t .
Now, we estimate the particular terms in J 1 ,
J 11 = | Ω t u u , z z 1 r ψ 1 , r z d x d t | | u | , Ω t | u , z z | 2 , Ω t | 1 r ψ 1 , r z | 2 , Ω t , J 12 = | Ω t u u , z r ψ 1 , r z z d x d t | | u | , Ω t | u , z r | 2 , Ω t | ψ 1 , r z z | 2 , Ω t , J 13 = | Ω t u u , z r ψ 1 , z z r d x d t | | u | , Ω t | u , z r | 2 , Ω t | ψ 1 , z z r | 2 , Ω t , J 14 = | Ω t 1 r 2 u u , z z ψ 1 , z d x d t | = | Ω t u u , z z ψ 1 , z r 2 d x d t | | u | , Ω t | u , z z | 2 , Ω t | ψ 1 , z r 2 | 2 , Ω t ,
where integration by parts can be performed in view of periodic boundary conditions on S 2 .
Next, we consider J 2 ,
J 2 = Ω t 1 r 2 u , r 2 ( r ψ 1 , z r + 2 ψ 1 , z ) r d r d z d t = Ω t 1 r u , r 2 ( r ψ 1 , z r + 2 ψ 1 , z ) d r d z d t = 0 t a a 1 r u u , r ( r ψ 1 , z r + 2 ψ 1 , z ) | r = 0 r = R d z d t Ω t u u , r r 1 r ψ 1 , z r + 2 r 2 ψ 1 , z d x d t Ω t u u , r ψ 1 , z r r 2 r 2 ψ 1 , z + 2 r ψ 1 , z r d r d z d t ,
where the boundary term for r = R vanishes because u | r = R = 0 . To examine the boundary term at r = 0 , we recall from [9] the expressions near the axis of symmetry
u = a 1 ( z , t ) r 2 + a 2 ( z , t ) r 4 + ,
so
u , r = 2 a 1 ( z , t ) r + 4 a 2 ( z , t ) r 3 +
Then,
1 r u u , r ( r ψ 1 , z r + 2 ψ 1 , z ) c r 2 ( r ψ 1 , z r + 2 ψ 1 , z ) .
The above expression vanishes for r = 0 because ψ 1 , z is bounded near the axis of symmetry.
Now, we estimate the particular terms in J 2 ,
J 21 = | Ω t u u , r r 1 r ψ 1 , r z d x d t | | u | , Ω t | u , r r | 2 , Ω t | 1 r ψ 1 , z r | 2 , Ω t , J 22 = | Ω t u u , r r 1 r 2 ψ 1 , z d x d t | | u | , Ω t | u , r r | 2 , Ω t | 1 r 2 ψ 1 , z | 2 , Ω t , J 23 = | Ω t u u , r r ψ 1 , z r r d x d t | | u | , Ω t | u , r r | 2 , Ω t | ψ 1 , z r r | 2 , Ω t , J 24 = | Ω t u u , r r 1 r 2 ψ 1 , z d x d t | | u | , Ω t | u , r r | 2 , Ω t | 1 r 2 ψ 1 , z | 2 , Ω t , J 25 = | Ω t u u , r r 1 r ψ 1 , z r d x d t | | u | , Ω t | u , r r | 2 , Ω t | 1 r ψ 1 , z r | 2 , Ω t .
Finally, we examine J 3 . Integrating by parts with respect to z, and using the periodic boundary conditions on S 2 , we have
J 3 = Ω t u 1 r 2 u , r z ( r ψ 1 , z z + 3 ψ 1 , r + r ψ 1 , r r ) d x d t + Ω t u 1 r 2 u , r ( r ψ 1 , z z z + 3 ψ 1 , r z + r ψ 1 , r r z ) d x d t .
Now, we estimate the particular terms in J 3 ,
J 31 = | Ω t u u , r z 1 r ψ 1 , z z d x d t | | u | , Ω t | u , r z | 2 , Ω t | ψ 1 , z z r | 2 , Ω t , J 32 = | Ω t u 1 r 2 u , r z ψ 1 , r d x d t | = | Ω t u r ε 0 u , r z ψ 1 , r r 2 ε 0 d x d t | | u | , Ω t 1 ε 0 | v φ | , Ω t ε 0 | u , r z | 2 , Ω t | ψ 1 , r r 2 | L 2 ( 0 , t ; L 2 , ε 0 ( Ω ) ) ,
where ε 0 > 0 can be as small as we want. Thus,
J 33 = | Ω t u r ε 0 u , r z 1 r 1 ε 0 ψ 1 , r r d x d t | | u | , Ω t 1 ε 0 | v φ | , Ω t ε 0 | u , r z | 2 , Ω t | ψ 1 , r r r 1 ε 0 | 2 , Ω t , J 34 = | Ω t u u , r r ψ 1 , z z z d x d t | | u | , Ω t | u , r r | 2 , Ω t | ψ 1 , z z z | 2 , Ω t , J 35 = | Ω t u u , r r 1 r ψ 1 , r z d x d t | | u | , Ω t | u , r r | 2 , Ω t | ψ 1 , r z r | 2 , Ω t , J 36 = | Ω t u u , r r ψ 1 , r r z d x d t | | u | , Ω t | u , r r | 2 , Ω t | ψ 1 , r r z | 2 , Ω t .
Summarizing the above estimates, we obtain
| J | c | u | , Ω t ( | u , z z | 2 , Ω t + | u , z r | 2 , Ω t + | u , r r | 2 , Ω t ) · · | 1 r ψ 1 , r z | 2 , Ω t + | 1 r ψ 1 , z z | 2 , Ω t + | 1 r 2 ψ 1 , z | 2 , Ω t + | u , r r | 2 , Ω t + | u , z r | 2 , Ω t | ψ 1 , r z z | 2 , Ω t + | ψ 1 , z r r | 2 , Ω t + | ψ 1 , z z z | 2 , Ω t + | 1 r ψ 1 , z z | 2 , Ω t + | 1 r ψ 1 , z r | 2 , Ω t + | 1 r 2 ψ 1 , z | 2 , Ω t ) + c | u | , Ω t 1 ε 0 | v φ | , Ω t ε 0 | u , r z | 2 , Ω t | ψ 1 , r r r | L 2 ( 0 , t ; L 2 , ε 0 ( Ω ) ) + | ψ 1 , r r 2 | L 2 ( 0 , t ; L 2 , ε 0 ( Ω ) ) .
Using (52), (159), (160), and the estimates from (46)
| u , r r | 2 , Ω t | v φ r | 2 , Ω t + | v φ , r | 2 , Ω t c D 1 , | u , z r | 2 , Ω t | v φ , z | 2 , Ω t c D 1
we obtain the following estimate for J,
| J | c [ D 2 ( D 3 + D 4 ) + D 1 D 2 ] | ψ 1 , r r z | 2 , Ω t + | ψ 1 , r z z | 2 , Ω t + | ψ 1 , z z z | 2 , Ω t + | 1 r ψ 1 , r z | 2 , Ω t + | 1 r ψ 1 , z z | 2 , Ω t + | 1 r 2 ψ 1 , z | 2 , Ω t ) + c D 2 1 ε 0 D 3 | v φ | , Ω t ε 0 | 1 r ψ 1 , r r | L 2 ( 0 , t ; L 2 , ε 0 ( Ω ) ) + | 1 r 2 ψ 1 , r | L 2 ( 0 , t ; L 2 , ε 0 ( Ω ) ) J .
From (64), we have the following (recall that ω 1 = Γ ):
| ψ 1 , r r z | 2 , Ω t + | ψ 1 , r z z | 2 , Ω t + | ψ 1 , z z z | 2 , Ω t c | Γ , z | 2 , Ω t .
estimates (82) and (85) imply
| ψ 1 , r z r | 2 , Ω t + | ψ 1 , z z r | 2 , Ω t + | ψ 1 , z r 2 | 2 , Ω t c | Γ , z | 2 , Ω t .
Finally, (93) yields
| 1 r ψ 1 , r r | L 2 ( 0 , t ; L 2 , ε 0 ( Ω ) ) + | 1 r 2 ψ 1 , r | L 2 ( 0 , t ; L 2 , ε 0 ( Ω ) ) c R ε 0 Γ 1 , 2 , Ω t
Recall that (77) is valid for ψ 1 | r = 0 = 0 .
This restriction implies that v z | r = 0 = 0 , so it is a strong restriction on solutions proved in this paper.
Using (176)–(178) in J yields
J c D 2 ( D 1 + D 3 + D 4 ) | Γ , z | 2 , Ω t + c D 2 1 ε 0 D 3 | v φ | , Ω t ε 0 Γ 1 , 2 , Ω t .
In view of Lemma 15, the term I 1 introduced in (175) is bounded by
I c | f φ | 2 , S 1 t u 2 , 2 , Ω t c | f φ | 2 , S 1 t ( D 3 + D 4 ) .
Using the estimates in (174), we obtain
ω r V ( Ω t ) 2 + ω z V ( Ω t ) 2 + | Φ | 2 , Ω t 2 c D 2 ( D 1 + D 3 + D 4 ) | Γ , z | 2 , Ω t + c D 2 1 ε 0 D 3 | v φ | , Ω t ε 0 Γ 1 , 2 , Ω t + c ( | F r | 6 / 5 , 2 , Ω t 2 + | F z | 6 / 5 , 2 , Ω t 2 ) + c ( | ω r ( 0 ) | 2 , Ω 2 + | ω z ( 0 ) | 2 , Ω 2 ) + c | f φ | 2 , S 1 t ( D 3 + D 4 ) ,
where we used
| Ω t ( F r ω r + F z ω z ) d x d t | ε ( | ω r | 6 , Ω 2 + | ω z | 6 , Ω 2 ) + c ( 1 / ε ) ( | F r | 6 / 5 , Ω 2 + | F z | 6 / 5 , Ω 2 ) .
Hence, (179) implies (173) and concludes the proof. □

7. Estimates for the Stream Function in Weighted Sobolev Spaces

Recall that the stream function ψ 1 is a solution to problem (22). To increase the regularity of weak solutions, we need appropriate estimates for ψ 1 , assuming sufficient regularity of the vorticity ω 1 .
Remark 7.
In Lemma 4, the existence of weak solutions to problem (22) satisfying estimate (56) is proven. Inequality (62) implies that the weak solution belongs to H 2 ( Ω ) and the estimate holds
ψ 1 2 , Ω c | ω 1 | 2 , Ω .
Assuming that ω 1 , z L 2 ( Ω ) , estimates (63) and (64) increase the regularity of ψ 1 , such that ψ 1 , z H 2 ( Ω ) and the estimate holds
ψ 1 , z H 2 ( Ω ) c ( | ω 1 , z | 2 , Ω + | ω 1 | 2 , Ω ) .
Estimate (181) is derived using the technique of the energy method. However, this method is not sufficiently robust to derive an estimate for | ψ 1 , r r r | 2 , Ω .
Moreover, estimate (181) is not sufficient to prove estimate (24) of Theorem 1. To prove Theorem 1 we need estimates (85) and (93). To prove the estimates, we need the theory of weighted Sobolev spaces developed by Kondratiev [10], which is used to examine elliptic boundary value problems in domains with cones.
Unfortunately, estimates (85) and (93) hold for weak solutions where ψ 1 vanishes on the axis of symmetry. This implies that the v z coordinate of velocity must also vanish on the axis of symmetry. Therefore, Theorem 1 is applicable to a smaller class than the class of weak solutions. This indicates that the regularity problem for axially symmetric solutions to the Navier–Stokes equations is only partially solved.
Now, we show the existence of solutions to problem (22) in weighted Sobolev spaces.
Lemma 17.
Assume that ψ 1 is a solution to (61). Assume that ω 1 , z , ω 1 L 2 ( Ω ) .
Then,
Ω ψ 1 , z r r 2 + 1 r 2 ψ 1 , z r 2 + 1 r 4 ψ 1 , z 2 d x + Ω ψ 1 , z z z 2 d x c Ω ( | ω 1 , z | 2 + | ω 1 | 2 ) d x .
Proof. 
To prove the lemma, we need weighted Sobolev spaces defined by Fourier transform (58) and introduced in (59) and (60). Therefore, to examine problem (22) in weighted Sobolev spaces, we have to derive estimates with respect to r and z, separately. To derive an estimate with respect to r, we have to examine solutions to (22) independently, as well in a neighborhood of the axis of symmetry, such as in a neighborhood located in a positive distance from it. To perform such considerations, we treat z as a parameter and introduce a partition of unity { ζ ( 1 ) ( r ) , ζ ( 2 ) ( r ) } , such that
i = 1 2 ζ ( i ) ( r ) = 1
and
ζ ( 1 ) ( r ) = 1 for r r 0 0 for r 2 r 0 , ζ ( 2 ) ( r ) = 0 for r r 0 1 for r 2 r 0 ,
where 0 < r 0 is fixed in such a way that 2 r 0 < R .
Let ψ 1 ( i ) = ψ 1 ζ ( i ) , ω 1 ( i ) = ω 1 ζ ( i ) and ζ ˙ ( i ) = d d r ζ ( i ) , ζ ¨ ( i ) = d 2 d r 2 ζ ( i ) , i = 1 , 2 . Moreover, functions ζ ( 1 ) and ζ ( 2 ) are smooth.
Then, we obtain from (22) the following two problems
Δ ψ 1 ( 1 ) 2 r ψ 1 ( 1 ) = ω 1 ( 1 ) 2 ψ 1 , r ζ ˙ ( 1 ) ψ 1 ζ ¨ ( 1 ) 2 r ψ 1 ζ ˙ ( 1 ) in Ω ( 1 ) , ψ 1 ( 1 ) satisfies periodic boundary conditions on S 2 ( 1 ) ,
where
Ω ( 1 ) = { ( r , z ) : r > 0 , z ( a , a ) } , S 2 ( 1 ) = { ( r , z ) : r > 0 , z { a , a } }
and
Δ ψ 1 ( 2 ) 2 r ψ 1 ( 2 ) = ω 1 ( 2 ) 2 ψ 1 , r ζ ˙ ( 2 ) ψ 1 ζ ¨ ( 2 ) 2 r ψ 1 ζ ˙ ( 2 ) in Ω ( 2 ) ψ 1 ( 2 ) = 0 on S 1 , ψ 1 ( 2 ) satisfies periodic boundary conditions on S 2 ( 2 ) ,
where
Ω ( 2 ) = { ( r , z ) : r 0 < r < R , z ( a , a ) } , S 2 ( 2 ) = { ( r , z ) : r 0 < r < R , z { a , a } } .
We temporarily simplify the notation using
u = ψ 1 ( 1 ) , w = ψ 1 ( 2 ) , f = ω 1 ( 1 ) 2 ψ 1 , r ζ ˙ ( 1 ) ψ 1 ζ ¨ ( 1 ) 2 r ψ 1 ζ ˙ ( 1 ) , b = ω 1 ( 2 ) 2 ψ 1 , r ζ ˙ ( 2 ) ψ 1 ζ ¨ ( 2 ) 2 r ψ 1 ζ ˙ ( 2 ) .
Then, (183) and (184) become
Δ u 2 r u , r = f in Ω ( 1 ) , u satisfies periodic boundary conditions on S 2 ( 1 )
and
Δ w 2 r w , r = b in Ω ( 2 ) , w = 0 on S 1 , w satisfies periodic boundary conditions on S 2 ( 2 ) .
Next, we use the simplified notation
p b c periodic boundary conditions .
First, we consider problem (186). We rewrite it in the form
u , r r 3 r u , r = f + u , z z in Ω ( 1 ) , u satisfies p b c on S 2 1 .
For a fixed z ( a , a ) , we treat (189) as
u , r r 3 r u , r = f + u , z z in R + .
Multiplying (190) by r 2 yields
r 2 u , r r 3 r u , r = r 2 ( f + u , z z ) g ( r , z )
or equivalently
r r ( r r u ) 2 r r u = g ( r , z ) .
Introduce the new variable
τ = ln r , r = e τ .
Since r r = τ , we see that (191) takes the form
τ 2 u + 2 τ u = g ( e τ , z ) = g ( τ , z ) .
Utilizing the Fourier transform (58) to (193), we have
λ 2 u ^ + 2 i λ u ^ = g ^ .
For λ { 0 , 2 i } we have
u ^ = 1 λ ( λ + 2 i ) g ^ R ( λ ) g ^ .
We introduce the quantity
h ( k , μ ) = 1 + k μ .
Consider the case where k = 0 and μ = 0 . Then, h ( 0 , 0 ) = 1 . Theorem 1.1 from Section 1 in [10] (also see Lemma 3.1 from [11]) yields the following:
Let f + u , z z L 2 ( R + ) , and R ( λ ) does not have poles on the line I m λ = 1 .
Then, we have
+ i h ( 0 , 0 ) + + i h ( 0 , 0 ) j = 0 2 | λ | 2 ( 2 j ) | u ^ | 2 d λ c + i h ( 0 , 0 ) + + i h ( 0 , 0 ) | g ^ | 2 d λ .
Using (60) and h ( 0 , 0 ) = 1 , we obtain
R j = 0 2 | τ j u | 2 e 2 τ d τ c R | g | 2 e 2 τ d τ .
Passing to variable r yields
R + | u , r r | 2 + 1 r 2 | u , r | 2 + 1 r 4 | u | 2 r d r c R + | f + u , z z | 2 r d r .
Using notation (185) and the estimate for the weak solutions, we obtain from (197) the inequality
R + supp ζ ( 1 ) ψ 1 , r r 2 + 1 r 2 ψ 1 , r 2 + 1 r 4 ψ 1 2 r d r c R + supp ζ ( 1 ) | ω 1 | 2 r d r + c R + supp ζ ( 1 ) ( | ψ 1 , r | 2 + | ψ 1 | 2 ) r d r + c R + supp ζ ( 1 ) | ψ 1 , z z | 2 r d r .
For solutions to (187), we have the estimate
w H 2 ( ( 0 , R ) supp ζ ( 2 ) ) c b L 2 ( ( 0 , R ) supp ζ ( 2 ) ) + c w , r L 2 ( ( 0 , R ) supp ζ ( 2 ) ) .
In view of notation (185) we obtain
( 0 , R ) supp ζ ( 2 ) ( | ψ 1 , r r | 2 + | ψ 1 , r | 2 + | ψ 1 | 2 ) r d r c ( 0 , R ) supp ζ ( 2 ) ( | ω 1 | 2 + | ψ 1 , z z | 2 + | ψ 1 , r | 2 + | ψ 1 | 2 ) r d r .
Adding (198) and (200), integrating the result with respect to z, and using (56) yield
Ω ψ 1 , r r 2 + 1 r 2 ψ 1 , r 2 + 1 r 4 ψ 1 2 d x c Ω ( | ω 1 | 2 + | ψ 1 , z z | 2 ) d x .
Replacing ψ 1 with ψ 1 , z , and ω 1 with ω 1 , z , we obtain estimate (182) from (201) and (63). This ends the proof. □
Lemma 18.
Assume that ψ 1 is a solution to (61). Assume that μ ( 0 , 1 ) , ω 1 H 1 ( Ω ) , Ω = ( 0 , R ) × ( a , a ) .
Then,
Ω ψ 1 , r r r 2 + 1 r 2 ψ 1 , r r 2 + 1 r 4 ψ 1 , r 2 + 1 r 6 ψ 1 2 r 2 μ d x + ψ 1 H 2 ( Ω ) 2 + Ω ( ψ 1 , z r r 2 + ψ 1 , z z r 2 + ψ 1 , z z z 2 ) d x c 1 + 1 μ 2 ω 1 H 1 ( Ω ) 2 .
Proof. 
Recall the partition of unity introduced in the proof of Lemma 17. Also recall the local problems (183) and (184), as well as notation (185). Then, we can examine problems (186) and (187). First, we examine problem (186).
Applying the Mellin transform, any solution to (190) can be expressed in the form (194).
In this case, we introduce the following quantity:
h ( 1 , μ ) = 2 μ .
Since operator R ( λ ) does not have poles on the line Im λ = h ( 1 , μ ) , we have (see Theorem 1.1 from Section 1 in [10])
+ i h ( 1 , μ ) + + i h ( 1 , μ ) j = 0 3 | λ | 2 ( 3 j ) | u ^ | 2 d λ c + i h ( 1 , μ ) + + i h ( 1 , μ ) j = 0 1 | λ | 2 ( 1 j ) | g ^ | 2 d λ .
Using (60) for h ( 1 , μ ) = 2 μ , we obtain
R j = 0 3 | τ j u | 2 e 2 ( 3 j ) τ d τ c R j = 0 1 | τ j g | 2 e 2 ( 1 j ) τ d τ .
In view of equivalence (59), inequality (205) takes the form
R + | u r r r | 2 + 1 r 2 | u r r | 2 + 1 r 4 | u r | 2 + 1 r 6 | u | 2 r 2 μ r d r c R + | ( f + u z z ) , r | 2 r 2 μ r d r + c R + | f + u z z | 2 r 2 μ 2 r d r ,
where z ( a , a ) and μ ( 0 , 1 ) .
Integrating (206) with respect to z, and exploiting notation (185), yield
a a d z R + supp ζ ( 1 ) ψ 1 , r r r 2 + 1 r 2 ψ 1 , r r 2 + 1 r 4 ψ 1 , r 2 + 1 r 6 ψ 1 2 r 2 μ r d r c a a d z R + supp ζ ( 1 ) ( | r ( ω 1 + ψ 1 , z z ) | 2 + | ω 1 + ψ 1 , z z | r 2 ) r 2 μ r d r .
For solutions to problem (187) and notation (185), we obtain
a a d z ψ 1 H μ 3 ( R + supp ζ ( 2 ) 2 c a a d z ( ω 1 H 1 ( R + supp ζ ( 2 ) ) 2 + ψ 1 , z z H 1 ( R + supp ζ ( 2 ) ) 2 ) .
From (207) and (208), as well as the Hardy inequality (see [14] (Chapter 1, Section 2.16)):
R + | ω 1 + ψ 1 , z z | 2 r 2 μ 2 r d r 1 μ 2 R + | ( ω 1 + ψ 1 , z z ) , r | 2 r 2 μ r d r
we obtain
Ω ψ 1 , r r r 2 + 1 r 2 ψ 1 , r r 2 + 1 r 4 ψ 1 , r 2 + 1 r 6 ψ 1 2 r 2 μ d x c 1 + 1 μ 2 ω 1 H 1 ( Ω ) 2 + Ω ( ψ 1 , z z r 2 + ψ 1 , z z z 2 ) d x .
Using estimates (56), (62), (63) in (210) implies (202) and ends the proof. □
Remark 8.
Since μ > 0 , the Hardy inequality (209) does not need ω 1 + ψ 1 , z z | r = 0 = 0 .

8. Conclusions

The main result of this paper is the proof of (24). Since Γ = ω φ / r , we obtain from (24) the estimate
ω φ V ( Ω t ) ϕ ( data ) ,
where we used the fact that r < R and R is finite. This means that (211) does not hold for the Cauchy problem.
Using problem (14) and relation (15), we obtain
v L ( 0 , t ; L 6 ( Ω ) ) c ψ L ( 0 , t ; H 2 ( Ω ) ) c ω φ L ( 0 , t ; L 2 ( Ω ) ) ϕ ( data ) ,
where v = ( v r , v z ) .
Consider the Stokes problem, which follows from (6)
v t ν Δ v + p = v · v + f in Ω T , div v = 0 in Ω T , v · n ¯ | S 1 = 0 , ω φ | S 1 = 0 , v φ | S 1 = 0 on S 1 T , v satisfies periodic boundary conditions on S 2 T , v | t = 0 = v ( 0 ) in Ω .
Using (212) and the energy estimate (46), we have
v · v L 2 ( 0 , t ; L 3 / 2 ( Ω ) ) ϕ ( data ) .
Assuming more regularity on data than in the proof of Theorem 1, and using the result from [15], we obtain the following estimate for solutions to (213)
v W 3 / 2 , 2 2 , 1 ( Ω t ) ϕ ( data ) .
Anisotropic weighted spaces can be found in [16].
The above inequality implies
v L 5 / 2 ( Ω t ) ϕ ( data )
thus, the increase in regularity holds because in (46) we have v L 2 ( Ω t ) ϕ ( data ) .
Continuing the considerations, we derive the estimate
v W 2 2 , 1 ( Ω t ) + p L 2 ( Ω t ) ϕ ( data ) .
The existence of solutions to problem (6) follows from appropriately choosing a fixed point theorem.
W p , q 2 , 1 ( Ω t ) is the Sobolev space with mixed norm. We have
u W p , q 2 , 1 ( Ω t ) = u L q ( 0 , t ; L p ( Ω ) ) + x 2 u L q ( 0 , t ; L p ( Ω ) ) + t u L q ( 0 , t ; L p ( Ω ) ) .

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the author ([email protected]).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Existence of Regular Solutions to (1) for Small Data

Recall the quantities
u 1 = v φ r , ω 1 = ω φ r , ψ 1 = ψ r , f 1 = f φ r , F 1 = F φ r .
In view of [12] system (1) is equivalent to the following one
u 1 , t + v · u 1 ν Δ u 1 + 2 r u 1 , r = 2 u 1 ψ 1 , z + f 1 , ω 1 , t + v · ω 1 ν Δ ω 1 + 2 r ω 1 , r = 2 u 1 u 1 , z + F 1 , Δ ψ 1 2 r ψ 1 , r = ω 1 , periodic boundary conditions on S 2 , u 1 | t = 0 = u 1 ( 0 ) . ω 1 | t = 0 = ω 1 ( 0 ) .
Functions u 1 , ω 1 , ψ 1 have compact support with respect to variable r.
We multiply (A2) 1 by u 1 | u 1 | 2 , integrate over Ω, and use boundary conditions to yield the following:
d d t | u 1 | 4 , Ω 4 + ν | u 1 | 4 , Ω 4 c | ω 1 | 2 , Ω 2 | u 1 | 4 , Ω 4 + c | f 1 | 4 , Ω 4 .
We multiply (A2) 2 by ω 1 , integrate over Ω, and exploit boundary conditions. Then, we have
d d t | ω 1 | 2 , Ω 2 + ν | ω 1 | 2 , Ω 2 c | u 1 | 4 , Ω 4 + c | F 1 | 2 , Ω 2 .
We introduce the quantity
X ( t ) = | u 1 ( t ) | 4 , Ω 4 + | ω 1 ( t ) | 2 , Ω 2 .
Then, (A3) and (A4) imply
d d t X + ν X c 0 X 2 + G ( t ) ,
where
G ( t ) = c ( | f 1 ( t ) | 4 , Ω 4 + | F 1 ( t ) | 2 , Ω 2 ) .
We consider (A6) on the time interval ( 0 , T ) . We assume that for t ( 0 , T ) , the following inequality holds
G ( t ) c 0 k 0 .
Then, (A6) takes the form
d d t X c 0 k 0 1 k 0 X 2 + 1 .
Let X = α X . Then,
d d t X c 0 k 0 α α 2 k 0 X 2 + 1 .
Setting α 2 = k 0 yields
d d t X c 0 k 0 ( X 2 + 1 ) .
Integrating (A10) with respect to time implies
arctg X ( t ) arctg X ( 0 ) c 0 k 0 t .
Hence,
X ( t ) tg ( c 0 k 0 t ) + X ( 0 ) 1 X ( 0 ) tg ( c 0 k 0 t ) .
Recalling that X ( 0 ) = X ( 0 ) k 0 and setting y = c 0 k 0 t , we obtain
X ( t ) ( tg y + X ( 0 ) k 0 ) k 0 1 X ( 0 ) c 0 t tg y y β ( T ) ,
where t T . Hence, for large T, (A11) holds for sufficiently small X ( 0 ) and k 0 .
Consider (A6) in the interval ( 0 , T ) . Using (A11), we can write (A6) in the form
d d t X + ν X G ( t ) ,
where ν = ν c 0 β .
Integrating (A12) with respect to time yields
X ( t ) e 0 ν t 0 t G ( t ) e ν t d t + e ν 8 t X ( 0 ) .
Setting t = T implies
X ( T ) 0 T G ( t ) d t + e ν T X ( 0 ) .
For sufficiently small X ( 0 ) and k 0 , the time interval ( 0 , T ) can be chosen to be large. Then, (A14) can imply that
X ( T ) X ( 0 ) .
Therefore, the previous considerations can be performed for any time interval ( k T , ( k + 1 ) T ) , k N .

References

  1. Ladyzhenskaya, O.A. Uniquae global solvability of the three-dimensional Cauchy problem for the Navier–Stokes equations in the presence of axial symmetry. Zap. Nauchn. Sem. LOMI 1968, 7, 155–177. (In Russian) [Google Scholar]
  2. Ukhovskii, M.R.; Yudovich, V.I. Axially symmetric flows of ideal and viscous fluids filling the whole space. J. Appl. Math. Mech. 1968, 32, 52–61. [Google Scholar]
  3. Nowakowski, B.; Zaja̧czkowski, W.M. Global regular axially-symmetric solutions to the Navier-Stokes equations with small swirl. J. Math. Fluid Mech. 2023, 25, 73. [Google Scholar] [CrossRef]
  4. Liu, Y. Solving the axisymmetric Navier-Stokes equations in critical spaces (I): The case with small swirl component. J. Diff. Equ. 2022, 314, 287–315. [Google Scholar] [CrossRef]
  5. Lei, Z.; Zhang, Q.S. Criticality of the axially symmetric Navier-Stokes equations. Pacific J. Math. 2017, 289, 169–187. [Google Scholar] [CrossRef]
  6. Zhang, P.; Zhang, T. Global axisymmetric solutions to three-dimensional Navier-Stokes system. Int. Math. Res. Not. 2014, 2014, 610–642. [Google Scholar] [CrossRef]
  7. Wei, D. Regularity criterion to the axially symmetric Navier-Stokes equations. J. Math. Anal. Appl. 2016, 435, 402–413. [Google Scholar] [CrossRef]
  8. Chen, H.; Fang, D.; Zhang, T. Regularity of 3d axisymmetric Navier-Stokes equations. Disc. Cont. Dyn. Syst. 2017, 37, 1923–1939. [Google Scholar] [CrossRef]
  9. Liu, J.G.; Wang, W.C. Characterization and regularity for axisymmetric solenoidal vector fields with application to Navier-Stokes equations. SIAM J. Math. Anal. 2009, 41, 1825–1850. [Google Scholar] [CrossRef]
  10. Kondrat’ev, V. Boundary value problems for elliptic equations in domains with conical or angular points. Trudy Moskov. Mat. Obshch. 1967, 16, 209–292. (In Russian); English Translation Trans. Mosc. Math. Soc. 1967, 16, 227–313. [Google Scholar]
  11. Nowakowski, B.; Zaja̧czkowski, W.M. On weighted estimates for the stream function of axially-symmetric solutions to the Navier-Stokes equations in a bounded cylinder. arXiv 2022, arXiv:2210.15729. [Google Scholar]
  12. Hou, T.Y.; Li, C. Dynamic stability of the three-dimensional axisymmetric Navier-Stokes equations with swirl. Comm. Pure Appl. Math. 2008, 61, 661–697. [Google Scholar] [CrossRef]
  13. Ladyzhenskaya, O.A. Mathematical Theory of Viscous Incompressible Flow, Nauka, Moscow 1970, 2nd English ed.; Silverman, R.A., Chu, J., Eds.; Mathematics and Its Applications; Gordon and Breach, Science Publishers: New York, NY, USA, 1970; Volume 2, p. xviii+224. (In Russian) [Google Scholar]
  14. Besov, O.V.; Il’in, V.P.; Nikolskii, S.M. Integral Representations of Functions and Imbedding Theorems, Nauka, Moscow (1975); English Translation; Scripta Series in Mathematics; V. H. Winston: New York, NY, USA, 1978; Volume I. (In Russian) [Google Scholar]
  15. Maremonti, P.; Solonnikov, V.A. On the estimates of evolution Stokes problem in anisotropic Sobolev spaces with mixed norm. Zap. Nauchn. Sem. LOMI 1994, 223, 124–150. [Google Scholar]
  16. Ladyzhenskaya, O.A.; Solonnikov, V.A.; Uraltseva, N.N. Linear and Quasilinear Equations of Parabolic Type, Nauka, Moscow (1967); American Mathematical Society: Providence, RI, USA, 1968. (In Russian) [Google Scholar]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Zaja̧czkowski, W.M. Global Regular Axially Symmetric Solutions to the Navier–Stokes Equations: Part 1. Mathematics 2023, 11, 4731. https://doi.org/10.3390/math11234731

AMA Style

Zaja̧czkowski WM. Global Regular Axially Symmetric Solutions to the Navier–Stokes Equations: Part 1. Mathematics. 2023; 11(23):4731. https://doi.org/10.3390/math11234731

Chicago/Turabian Style

Zaja̧czkowski, Wojciech M. 2023. "Global Regular Axially Symmetric Solutions to the Navier–Stokes Equations: Part 1" Mathematics 11, no. 23: 4731. https://doi.org/10.3390/math11234731

APA Style

Zaja̧czkowski, W. M. (2023). Global Regular Axially Symmetric Solutions to the Navier–Stokes Equations: Part 1. Mathematics, 11(23), 4731. https://doi.org/10.3390/math11234731

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop