#### Abstract

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 while the diameter of the helix is larger, of order (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.

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

no longer equals to zero due to pressure boundary condition prescribed on , (in our case . see Figure 1) Here denotes the velocity of the fluid which occupies a bounded domain , while stands for the unit outward normal on . 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 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 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 Problem

##### 2.1. Pipe's Geometry

We start by defining the helix whose parametrization has the form 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 It can be easily verified that the curvature (flexion) and the torsion are constant and given by For small parameter and a unit circle we introduce a thin straight pipe with circular cross-section: We define the mapping by and we put 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

##### 2.2. The Equations

As mentioned before, the fluid inside of the pipe is assumed to be Newtonian so the velocity and pressure satisfy the following Navier-Stokes equations: Here 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: The fluid flow is governed by the pressure drop between pipe's ends so in (2.11) we prescribe constant pressures and , . 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 , 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 Solution

##### 3.1. Variational Formulation

Let us introduce the following natural functional space: The space is equipped with the norm which is equivalent to the -norm due to Poincaré's inequality. Now we can write the variational formulation of our problem (2.8)–(2.11): Notice that we have eliminated the pressure . Furthermore, it should be observed that replacing by changes nothing so the relevant quantity is the pressure drop . As a consequence, the pressure 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 , independent of , such that for all satisfying 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 be such that on . We introduce new functions with . It is clear that , where and on . Now we extend by zero on a strip-like domain , where is chosen such that . A simple integration formula yields Using the Cauchy-Schwarz inequality, one can easily obtain Since on , it follows that with independent of By a simple change of variables , we get Direct calculation gives , where implying the asymptotic behavior . In view of that, by simple change of variables, we obtain Now from (3.9)–(3.12) we deduce implying (3.4). Finally, using the interpolation inequality, the embedding and the estimate (3.4) we obtain at once

Our main result can be stated as follows.

Theorem 3.3. Assume that the pressure drop and the helix step are such that Then the problem (3.3) admits at least one solution . Moreover, such solution is unique in the ball: where are the constants from Lemma 3.1.

Proof. The idea of the proof is to introduce the mapping defined by , where is the solution of the variational problem: (Although the operator depends on we drop the index for the sake of notational simplicity.)The first step is to show that mapping is well defined. For that purpose, let be the following bilinear form Using estimate (3.5), we get Since we conclude implying that the form is elliptic on . Analogously, we can easily obtain that Thus, due to Lax-Milgram Lemma, problem (3.17) has a unique solution and is well defined on .
The next step is to prove that . Applying (3.4), we obtain From (3.17) and (3.20), it follows that Taking into account the condition (3.15), we obtain proving that .
Now it remains to prove that is a contraction which will enable us to use Banach fixed point theorem. Let , be such that , , with and being the solutions of problem (3.17) for and respectively. Subtracting the corresponding equations it follows that Setting , we obtain For the first integral, we have Here we used the estimate (3.5) and the fact that . For the second integral, we proceed in the similar way and obtain leading to implying that is a contraction. Now Banach fixed point theorem provides the existence and uniqueness of 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 is unique in some ball 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 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.