We consider a flow of incompressible Newtonian fluid through
a pipe with helical shape. We suppose that the flow is governed by the prescribed
pressure drop between pipe's ends. Such model has relevance to some important
engineering applications. Under small data assumption, we prove the existence and
uniqueness of the weak solution to the corresponding Navier-Stokes system with pressure
boundary condition. The proof is based on the contraction method.
1. Introduction
Engineering practice requires extensive knowledge of flow through curved pipes. Helically coiled pipes are well-known types of curved pipes which have been used in wide variety of applications. They cover a large number of devices such as pipelines, air conditioners, refrigeration systems, central heating radiators, and chemical reactors. Therefore, numerous researchers have studied the fluid flow in helical pipes with circular cross-section both theoretically and experimentally. Let us just mention some of them. Wang [1] proposed a nonorthogonal coordinate system to investigate the effects of curvature and torsion on the low-Reynolds number flow in a helical pipe. The introduction of an orthogonal system of coordinates along a spatial curve allowed Germano [2, 3] to explore in a simpler way the effects of pipe's geometry on a helical pipe flow. Yamamoto et al. [4] first studied experimentally the effects of torsion and curvature on the flow characteristics in a helical tube. The results obtained by the experiments were then compared with those obtained from the model of Yamamoto et al. in [5]. Hüttl and Friedrich [6, 7] applied the second-order finite volume method for solving the incompressible Navier-Stokes equations to study the turbulent flow in helically coiled tubes. More numerical simulations on helical pipe flow can be found in the works of Yamamoto et al. [8] and Wang and Andrews [9].
In this paper, we study the stationary flow of incompressible Newtonian fluid through a helical pipe with prescribed pressures at its ends. We suppose that the pipe's thickness and the helix step have the same small order 𝒪(ε),0<ε≪1, while the diameter of the helix is larger, of order 𝒪(1) (see Figure 1). Such assumptions cover a large variety of realistic coiled pipes and they appear naturally in many devices as, for instance, the Liebig cooler, cooling channel in the nozzle of a rocket engine, particle separators used in the mineral processing industry, and so forth.
Helical pipe.
It is well-known that the flow of incompressible viscous Newtonian fluid is described by Navier-Stokes system of nonlinear PDEs. Their nonlinearity makes most boundary-value problems difficult or impossible to solve and thus the concept of weak solutions is introduced. When the velocity is prescribed on the whole boundary (Dirichlet condition), the existence of the weak solution can be proved by constructing approximate solutions via Galerkin method (see, e.g., Temam [10], although the main ideas goes back to Leray [11]). However, applications (as this one in the present paper) often give rise to problems where it is natural to prescribe the value of pressure on some part of the boundary. In case of (linear) Stokes system with pressure boundary condition, the existence and uniqueness of weak solution is obtained in Conca et al. [12]. For the Navier-Stokes system with boundary conditions involving pressure we only have some partial results and the full proof of existence is still unknown. The main difficulty lies in the fact that the integral
∫Ω(u∇)uu=12∫Γp|u|2(u·n)
no longer equals to zero due to pressure boundary condition prescribed on Γp⊂∂Ω, (in our case Γp=Σ0ε∪Σℓε. see Figure 1) Here u denotes the velocity of the fluid which occupies a bounded domain Ω⊂R3, while n stands for the unit outward normal on Γp. As a consequence, the nonlinear term in corresponding variational formulation does not vanish (and we do not know how to control it) causing the absence of the energy equality. Described technical difficulty can be elegantly overcome by prescribing the so-called dynamic (Bernoulli) pressure p+(1/2)|u|2 as proposed by Conca et al. [12] (see also Łukaszewicz [13] for nonstationary flow). Nevertheless, it should be mentioned that there is no physical justification for prescribing such boundary condition in the case of viscous fluid. In the case of ideal fluid the interpretation of the dynamic pressure is given by the Bernoulli law.
Another possibility is to restrict to the case of small boundary data as in Heywood et al. [14] or Marušić-Paloka [15] (see also Marušić-Paloka [16]) and we follow such approach. By doing that, we manage to control the inertial term in Navier-Stokes equations and to obtain existence and uniqueness of the solution for the physically realistic situation (without prescribing the dynamic pressure).
The paper is organized as follows. First we describe pipe's geometry using parametrization rε(x1)=(x1,acos(x1/ε),asin(x1/ε)) of its central curve. Here the small parameter ε stands for the distance between two coils of the helix. We assume that the pipe has (constant) circular cross-section of size ε and use Frenet basis attached to a helix to formally define our domain. Then we state the boundary-value problem describing the flow of a Newtonian fluid inside of the pipe and introduce the corresponding variational formulation. The main result is formulated in Theorem 3.3 where we precisely establish the assumption on the prescribed pressure drop under which the weak solution exists and then prove its uniqueness. The proof is based on the auxiliary result on sharp Sobolev constants (Lemma 3.1) and the contraction method.
We end this introduction by giving few more bibliographic remarks. The asymptotic behavior of the fluid flow through a helical pipe was investigated in Marušić-Paloka and Pažanin [17, 18]. Using the techniques from Marušić-Paloka [19], enabling the treatment of the curved geometry, the asymptotic approximation of the solution is built and rigorously justified by proving the error estimate in terms of the small parameter ε. Last but not least, let us mention that in Marušić-Paloka and Pažanin [20] the nonisothermal flow of Newtonian fluid through a general curved pipe has been considered. In such flow temperature changes cannot be neglected so the Navier-Stokes equations are coupled with heat conducting equation. A simplified model showing explicitly the effects of pipe's geometry is derived via rigorous asymptotic analysis with respect to the pipe's thickness.
2. Setting of the Problem2.1. Pipe's Geometry
We start by defining the helix whose parametrization has the form
rε(x1)=(x1,acosx1ε,asinx1ε),x1∈[0,ℓ],
serving to define the center curve of the pipe. It should be noted that the curve is not parameterized by its arc length (i.e., we do not use the natural parametrization) since in that case interval length ℓ would depend on the small parameter ε which is inconvenient for further analysis. At each point of the helix we compute the Frenet basis as follows
tε(x1)=rε′(x1)|rε′(x1)|=1a2+ε2(ε,-asinx1ε,acosx1ε),nε(x1)=tε′(x1)|tε′(x1)|=(0,-cosx1ε,-sinx1ε),bε(x1)=tε(x1)×nε(x1)=1a2+ε2(a,εsinx1ε,-εcosx1ε).
It can be easily verified that the curvature (flexion) κ and the torsion τ are constant and given by
κ=aa2+ε2,τ=εa2+ε2.
For small parameter ε>0 and a unit circle B=B(0,1)⊂R2, we introduce a thin straight pipe with circular cross-section:
Sε={x=(x1,x2,x3)∈R3:x1∈(0,ℓ),x'=(x2,x3)∈εB}=(0,ℓ)×εB.
We define the mapping Φε:Sε→R3 by
Φε(x)=rε(x1)+x2nε(x1)+x3bε(x1),
and we put
Ωε=Φε(Sε).
Such domain is our thin pipe with helical shape filled with a viscous incompressible fluid. Finally, the lateral boundary of the pipe and its ends are denoted by
Γε=Φε((0,ℓ)×ε∂B),Σiε=Φε({i}×εB),i=0,ℓ.
2.2. The Equations
As mentioned before, the fluid inside of the pipe is assumed to be Newtonian so the velocity uε and pressure pε satisfy the following Navier-Stokes equations:
-μΔuε+(uε∇)uε+∇pε=0inΩε,divuε=0inΩε.
Here μ>0 and we have added the subscript ε in our notation in order to stress the dependence of the solution on the small parameter. The above system must be completed by the boundary conditions:
uε=0onΓε,uε×tε=0,onΣiε,pε=qionΣiε,i=0,ℓ.
The fluid flow is governed by the pressure drop between pipe's ends so in (2.11) we prescribe constant pressures q0 and qℓ, qℓ<q0. In addition to the value of pressure, we need to prescribe something more on the boundary in order to assure that the problem is well possed. Thus, we take the tangential velocity to be zero on Σiε, while we keep the classical no-slip condition for the velocity on Γε. Imposing that the tangential component of the velocity equals to zero is not a serious restriction since the only part that counts is the normal part, due to the Saint-Venant principle for thin domains (see, e.g., Marušić-Paloka [19]). Let us mention that instead of these conditions, one could prescribe the whole normal stress, including the viscous part, on the boundary (see, e.g., Heywood et al. [14]). Such situation can be treated using the same method, with a slight change of functional space in the corresponding variational formulation.
3. Existence and Uniqueness of the Solution3.1. Variational Formulation
Let us introduce the following natural functional space:
Vε={v∈H1(Ωε)3:divv=0inΩε,v=0onΓε,v×tε=0onΣiε,i=0,ℓ}.
The space Vε is equipped with the norm
∥v∥Vε=∥∇v∥L2(Ωε),
which is equivalent to the H1(Ωε)3-norm due to Poincaré's inequality. Now we can write the variational formulation of our problem (2.8)–(2.11):
finduε∈Vεsuchthatμ∫Ωε∇uε∇v+∫Ωε(uε∇)uεv+q0∫Σ0ε(v·tε)+qℓ∫Σℓε(v·tε)=0∀v∈Vε.
Notice that we have eliminated the pressure pε. Furthermore, it should be observed that replacing qi by qi+C(C=const.)i=0,ℓ changes nothing so the relevant quantity is the pressure drop q0-qℓ. As a consequence, the pressure pε is determined only up to an additive constant. The equivalence between the variational (3.3) and differential (2.8)–(2.11) formulation, in case of smooth solutions, is discussed in the usual way (see, e.g., Conca et al. [12]).
3.2. The Main Result
The goal of this paper is to prove the existence and uniqueness of the solution for the variational problem (3.3), under small data assumption. In order to accomplish that, we need the following auxiliary result.
Lemma 3.1.
There exist constants C1,C2>0, independent of ε, such that
∥φ∥L2(Ωε)≤C1ε∥∇φ∥L2(Ωε),∥φ∥L4(Ωε)≤C2ε1/4∥∇φ∥L2(Ωε),
for all φ∈H1(Ωε) satisfying φ=0 on Γε.
Remark 3.2.
It is well known that Poincaré’s constant depends on the geometry of the domain Ωε. Inequality (3.4) gives the precise dependence of that constant on the small parameter ε. Inequality (3.5) enable us to estimate the inertial term in the variational formulation (3.3).
Proof of Lemma 3.1.
Let φ∈H1(Ωε) be such that φ=0 on Γε. We introduce new functions
Θ(x)=φ(z),Ψ(x1,y')=Θ(x1,εy'),
with z=Φε(x). It is clear that Ψ∈H1(S), where S=(0,ℓ)×B and Ψ=0 on Γ=(0,ℓ)×∂B. Now we extend Ψ by zero on a strip-like domain Sa,b=(0,ℓ)×(a,b)2, where (a,b) is chosen such that B⊂(a,b)2. A simple integration formula yields
Ψ(x1,y')=∫ay2∂Ψ∂y2(x1,ξ,y3)dξ.
Using the Cauchy-Schwarz inequality, one can easily obtain
∫Sa,b|Ψ(x1,y')|2dx1dy'≤(b-a)22∫Sa,b|∂Ψ∂y2(x1,ξ,y3)|2dx1dξdy3.
Since Ψ=0 on Sa,b∖S, it follows that
∥Ψ∥L2(S)≤C∥∂Ψ∂y2∥L2(S),
with C independent of ε. By a simple change of variables (x'=εy'), we get
∫Sε|Θ|2dx=ε2∫S|Ψ|2dx1dy',∫Sε|∂Θ∂x2|2dx=∫S|∂Ψ∂y2|2dx1dy'.
Direct calculation gives det∇Φε=gε, where
gε=1+(x2-a)2ε2-x22a2+ε2,
implying the asymptotic behavior det∇Φε=gε=(a/ε)1+𝒪(ε). In view of that, by simple change of variables, we obtain
∫Ωε|φ|2dz=∫Sε|Θ|2gεdx≤Cε∫Sε|Θ|2dx,∫Sε|∂Θε∂x2|2dx≤∫Sε|∇Θ|2gε1gεdx≤Cε∫Ωε|∇φ|2dz.
Now from (3.9)–(3.12) we deduce
∫Ωε|φ|2dz≤Cε2∫Ωε|∇φ|2dz,
implying (3.4). Finally, using the interpolation inequality, the embedding H1(Ωε)↪L6(Ωε), and the estimate (3.4) we obtain at once
∥φ∥L4(Ωε)≤∥φ∥L2(Ωε)1/4∥φ∥L6(Ωε)3/4≤Cε1/4∥∇φ∥L2(Ωε).
Our main result can be stated as follows.
Theorem 3.3.
Assume that the pressure drop and the helix step are such that
ε2(q0-qℓ)<2μ2ℓ9C1C22.
Then the problem (3.3) admits at least one solution uε∈Vε. Moreover, such solution is unique in the ball:
Bε={v∈Vε:∥∇v∥L2(Ωε)≤μ3C22ε},
where C1,C2>0 are the constants from Lemma 3.1.
Proof.
The idea of the proof is to introduce the mapping T:Bε→H1(Ωε)3 defined by T(w)=u, where u is the solution of the variational problem:
findu∈Vεsuchthatμ∫Ωε∇u∇v+∫Ωε(w∇)uv+q0∫Σ0ε(v·tε)+qℓ∫Σℓε(v·tε)=0,∀v∈Vε.
(Although the operatorT depends onε we drop the indexε for the sake of notational simplicity.)The first step is to show that mapping T is well defined. For that purpose, let a:Vε×Vε→R be the following bilinear form
a(u,v)=μ∫Ωε∇u∇v+∫Ωε(w∇)uv.
Using estimate (3.5), we get
a(u,u)≥(μ-C22ε∥∇w∥L2(Ωε))∥∇u∥L2(Ωε)2.
Since w∈Bε, we conclude
a(u,u)≥2μ3∥∇u∥L2(Ωε)2,∀u∈Vε,
implying that the form a(·,·) is elliptic on Vε. Analogously, we can easily obtain that
|a(u,v)|≤4μ3∥∇u∥L2(Ωε)∥∇v∥L2(Ωε).
Thus, due to Lax-Milgram Lemma, problem (3.17) has a unique solution and T is well defined on Bε.
The next step is to prove that T(Bε)⊂Bε. Applying (3.4), we obtain
|q0∫Σ0ε(v·tε)+qℓ∫Σℓε(v·tε)|=|∫Ωεdiv((q0+qℓ-q0ℓx1)v)|≤q0-qℓℓC1εε∥∇v∥L2(Ωε).
From (3.17) and (3.20), it follows that
∥∇T(w)∥L2(Ωε)=∥∇u∥L2(Ωε)≤32μq0-qℓℓC1εε.
Taking into account the condition (3.15), we obtain
∥∇T(w)∥L2(Ωε)<μ3C22ε,
proving that T(Bε)⊂Bε.
Now it remains to prove that T is a contraction which will enable us to use Banach fixed point theorem. Let w1,w2∈Bε, w1≠w2 be such that T(w1)=u1, T(w2)=u2, with u1 and u2 being the solutions of problem (3.17) for w1 and w2, respectively. Subtracting the corresponding equations it follows that
μ∫Ωε∇(T(w1)-T(w2))∇v+∫Ωε((w1∇)T(w1)-(w2∇)T(w2))v=0,∀v∈Vε.
Setting v=u1-u2, we obtain
μ∥∇(T(w1)-T(w2))∥L2(Ωε)2=-∫Ωε((w1-w2)∇)T(w1)(T(w1)-T(w2))-∫Ωε(w2∇)(T(w1)-T(w2))(T(w1)-T(w2)).
For the first integral, we have
|∫Ωε((w1-w2)∇)T(w1)(T(w1)-T(w2))|≤C22ε∥∇(w1-w2)∥L2(Ωε)∥∇T(w1)∥L2(Ωε)∥∇(T(w1)-T(w2))∥L2(Ωε)≤μ3∥∇(w1-w2)∥L2(Ωε)∥∇(T(w1)-T(w2))∥L2(Ωε).
Here we used the estimate (3.5) and the fact that T(w1)∈Bε. For the second integral, we proceed in the similar way and obtain
|∫Ωε(w2∇)(T(w1)-T(w2))(T(w1)-T(w2))|≤μ3∥∇(T(w1)-T(w2))∥L2(Ωε)2,
leading to
∥∇(T(w1)-T(w2))∥L2(Ωε)≤12∥∇(w1-w2)∥L2(Ωε),
implying that T is a contraction. Now Banach fixed point theorem provides the existence and uniqueness of uε finishing the proof.
Remark 3.4.
Regarding the uniqueness of the solution, it should be observed that we managed to prove only that the weak solution uε is unique in some ball Bε around zero, under small data assumption (3.15). It means that we cannot exclude the existence of some other solutions with large norm. The existence and uniqueness (up to an additive constant) of the pressure p∈L2(Ωε) satisfying the governing equations in the sense of the distributions can be obtained in the standard way using De Rham theorem (see, e.g., Temam [10]).
Acknowledgment
This research was supported by the Ministry of Science, Education and Sports, Republic of Croatia, Grant 037-0372787-2797.
WangC. Y.On the low-reynolds-number flow in a helical pipe1981108185194EID2-s2.0-0019584472GermanoM.On the effect of torsion on a helical pipe flow198212518EID2-s2.0-0020250386GermanoM.The dean equations extended to a helical pipe flow1989203289305EID2-s2.0-0024680708YamamotoK.AkitaT.IkeuchiH.KitaY.Experimental study of the flow in a helical circular tube1995164237249EID2-s2.0-0029360919YamamotoK.YanaseS.JiangR.Stability of the flow in a helical tube1998223153170MR1606161ZBL1051.76560EID2-s2.0-0042734108HüttlT. J.FriedrichR.Influence of curvature and torsion on turbulent flow in helically coiled pipes200021334535310.1016/S0142-727X(00)00019-9EID2-s2.0-0343962734HüttlT. J.FriedrichR.Direct numerical simulation of turbulent flows in curved and helically coiled pipes200130559160510.1016/S0045-7930(01)00008-1EID2-s2.0-0035370695YamamotoK.AribowoA.HayamizuY.HiroseT.KawaharaK.Visualization of the flow in a helical pipe200230425126710.1016/S0169-5983(02)00043-6EID2-s2.0-0036542289WangJ.-W.AndrewsJ. R. G.Numerical simulation of flow in helical ducts199541510711080EID2-s2.0-0029158547TemamR.19772Amsterdam, The NetherlandsNorth-Hollandx+500Studies in Mathematics and Its ApplicationsMR0609732LerayJ.Etude de diverses équations intégrales non linéaires et du quelques
problemes que pose l'hydrodynamique193312182ConcaC.MuratF.PironneauO.The Stokes and Navier-Stokes equations with boundary conditions involving the pressure1994202279318MR1308419ZBL0826.35093ŁukaszewiczG.On the Navier-Stokes equations in time dependent domains and with boundary conditions involving the pressure199743529550MR1484601ZBL0905.35062HeywoodJ. G.RannacherR.TurekS.Artificial boundaries and flux and pressure conditions for the incompressible Navier-Stokes equations1996225325352EID2-s2.0-0029777695MR138084410.1002/(SICI)1097-0363(19960315)22:5<325::AID-FLD307>3.0.CO;2-YZBL0863.76016Marušić-PalokaE.Solvability of the Navier-Stokes system with L2 boundary data2000413365375MR1739399ZBL0952.35090Marušić-PalokaE.Rigorous justification of the Kirchhoff law for junction of thin pipes filled with viscous fluid20033315166MR1976050ZBL1114.76320Marušić-PalokaE.PažaninI.Effective flow of a viscous liquid through a helical pipe20043321297397810.1016/j.crme.2004.09.005EID2-s2.0-9644301021Marušić-PalokaE.PažaninI.Fluid flow through a helical pipe20075818199MR2293105ZBL1115.35099Marušić-PalokaE.The effects of flexion and torsion on a fluid flow through a curved pipe200144324527210.1007/s00245-001-0021-yEID2-s2.0-18844468308MR1851740ZBL1007.35060Marušić-PalokaE.emarusic@math.hrPažaninI.pazanin@math.hrOn the effects of curved geometry on heat conduction through a distorted pipeNonlinear Analysis: Real World Applications. In pressEID2-s2.0-6054910788610.1016/j.nonrwa.2008.09.016