Next Article in Journal
Application of a Real-Time Visualization Method of AUVs in Underwater Visual Localization
Next Article in Special Issue
Slag as an Inventory Material for Heat Storage in a Concentrated Solar Tower Power Plant: Design Studies and Systematic Comparative Assessment
Previous Article in Journal
Image Recovery with Data Missing in the Presence of Salt-and-Pepper Noise
Previous Article in Special Issue
Nucleation Triggering of Highly Undercooled Xylitol Using an Air Lift Reactor for Seasonal Thermal Energy Storage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamics of a Partially Confined, Vertical Upward-Fluid-Conveying, Slender Cantilever Pipe with Reverse External Flow

1
State Key Laboratory of Geomechanics and Geotechnical Engineering, Institute of Rock and Soil Mechanics, Chinese Academy of Sciences, Wuhan 430071, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
State Key Laboratory of Coal Mine Disaster and Control, Chongqing University, Chongqing 400044, China
*
Authors to whom correspondence should be addressed.
Appl. Sci. 2019, 9(7), 1425; https://doi.org/10.3390/app9071425
Submission received: 18 February 2019 / Revised: 28 March 2019 / Accepted: 2 April 2019 / Published: 4 April 2019
(This article belongs to the Special Issue Clean Energy and Fuel (Hydrogen) Storage)

Abstract

:
A linear theoretical model is established for the dynamics of a hanging vertical cantilevered pipe which is subjected concurrently to internal and reverse external axial flows. Such pipe systems may have instability by flutter (amplified oscillations) or static divergence (buckling). The pipe system under consideration is a slender flexible cantilevered pipe hanging concentrically within an inflexible external pipe of larger diameter. From the clamped end to the free end, fluid is injected through the annular passage between the external pipe and the cantilevered pipe. When exiting the annular passage, the fluid discharges in the counter direction along the cantilevered pipe. The inflexible external pipe has a variable length and it can cover a portion of the length of the cantilevered pipe. This pipe system has been applied in the solution mining and in the salt cavern underground energy storage industry. The planar motion equation of the system is solved by means of a Galerkin method, and Euler–Bernoulli beam eigenfunctions are used as comparison functions. Calculations are conducted to quantify the effects of different confinement conditions (i.e., the radial confinement degree of the annular passage and the confined-flow length) on the cantilevered pipe stability, for a long leaching-tubing-like system. For a long system, an increase in the radial confinement degree of the annular passage and the confined-flow length gives rise to a series of flutter and divergence. Additionally, the effect of the cantilevered pipe length is studied. Increasing the cantilevered pipe length results in an increase of the critical flow velocity while a decrease of the associated critical frequency. For a long enough system, the critical frequency almost disappears.

1. Introduction

This paper focuses on flow-induced vibration of a slender flexible cantilevered pipe which is hanging concentrically inside a shorter inflexible external pipe, thus constituting an annular passage between the cantilevered pipe and the shorter external pipe, as shown in Figure 1. The entire pipe system is submerged in incompressible fluid and is placed within a closed cavity, with the upper portion outside the cavity. As shown in Figure 1a, from the fixed end to the free end of the shorter inflexible external pipe, fluid is injected as external flow into the closed cavity through the annular passage, and flows out as a free jet into stagnant fluid. In consideration of conservation of mass, initially stagnant fluid in the closed cavity is then forced through the cantilevered pipe, as a bounded reverse internal flow, discharging from the closed cavity upwards along the cantilevered pipe. Comparing the pipe configuration shown in Figure 1a with that given in Figure 1b, it is clear that the fluid flowing directions are opposite from each other. Contrary to the direction of fluid circulation illustrated in Figure 1a, fluid [1] is injected as internal flow into the closed cavity through the cantilever pipe and discharges through the annular passage, as shown in Figure 1b.
The pipe configuration shown in Figure 1 has been applied in solution mining and in the salt cavern underground energy storage industry. Figure 2 shows the schematic view of pipe configurations applied in the salt cavern underground energy storage. Figure 2a,b correspond to Figure 1a,b, respectively. In fact, the pipe failure or excessive bending caused by flow-induced vibrations is a well-known issue in the salt cavern storage industry [2,3]. For example, excessive bending of the inner tubing in a Chinese salt cavern storage is shown in Figure 3. This paper focuses on the pipe model shown in Figure 1a, and the motivation is desired to potentially improve the design of slender cantilever pipes in order to avoid issues cause by flow-induced vibration.
For several decades, a few researchers have been done the work on conveying-fluid pipes subjected simultaneously to both internal and external axial flows. Most of the work done involves concurrent flows, i.e., internal and external axial flows in the same direction. In [4], Cesari and Curioni investigated the buckling instability of pipes subjected to internal and external axial fluid flow. In [5], Hannoyer et al. investigated the dynamics and stability of both clamped–clamped cylindrical tubular beams conveying fluid, which simultaneously are subjected to independent axial external flows. They also studied the dynamics of cantilevered pipes fitted with a tapered nozzle at the free end. They provided theoretical and experimental results to support their model. In [6], Païdoussis and Besancon studied various aspects of the dynamic characteristic of clusters of tubular pipes which were subjected to internal flows and concurrently surrounded by bounded outer axial flows. To establish general characteristics of free motions, they obtained the eigenfrequencies of the system and studied their evolution by increasing either external or internal flows. In [7], Wang and Bloom established a mathematical model to investigate the vibration and stability of a submerged and inclined concentric pipe system which subjected to internal and external flow, and obtained the resonant frequencies of the system. In [8], Païdoussis et al. developed a theoretical model to study the vibration of a hanging tubular cantilever which was centrally inside a cylindrical container, with fluid flowing within the cantilever, and discharging from the free end. The configuration thus resembles that of a drill-string with a floating fluid-powered drill-bit. Of particular interest is the study in [1]. Moditis and Païdoussis preliminarily discussed the dynamics of the pipe configuration shown in Figure 1b. Furthermore, they developed a theoretical model and carried out corresponding experiments.
Because the length of the cantilevered pipe system in practical applications is on the order of one kilometer, how the system behavior evolves is of interest as the length increases. Theoretical and experimental studies [9,10] regarding long hanging vertical pipes conveying fluid have shown that both the critical flow velocity and associated frequency for instability tend to be asymptotic toward different limiting values with increasing the length of the pipes. The literature related to aspirating cantilevered pipes was believed to be useful for analyzing the dynamics of a vertical hanging pipe shown in Figure 1a. For example, in [11,12], the authors theoretically and experimentally considered the dynamic stability of a submerged cantilever pipe aspirating fluid, which could be applied in deep ocean mining.
Several key dimensions of the system under consideration are defined in Figure 1a. The x-axis lies in the undeflected centerline of the internal cantilever pipe, while the z-axis lies along the lateral direction and is perpendicular to the x-axis. L is the length of the cantilever pipe, while Ls is the length of the shorter inflexible external pipe (i.e. the length of the annular passage). Ui is the internal flow velocity (area average flow velocity) inside the cantilevered pipe, and Uo is the external flow velocity (area average flow velocity) in the annular passage with reference to the cantilever pipe.
Firstly, the present paper studies the linear idealized system shown in Figure 1a. We present a derivation of a theoretical model as well as a method of solving it are discussed. Secondly, we give theoretical results for slender leaching-tubing-like systems which are associated with the industrial application of salt cavern storage. We also obtained some unexpected and interesting results by numerical simulation. Finally, we present general conclusions.

2. Problem Formulation

2.1. Derivation of the Motion Equation of Theoretical Model

The derivation of the linear theoretical model is carried out as follows. Take a small element of length δx of the cantilevered pipe into consideration, as shown in Figure 4a, under the action of fluid-related and structural forces and moments. Force balances in the x- and z-directions renders, respectively
F T x x ( Q w x ) + M p g F i t + F e t ( F i n + F e n ) w x = 0
Q x + x ( F T w x ) + F i n + F e n + ( F e t F i t ) w x M p 2 w t 2 = 0
where w is the lateral deflection; FT is the axial tension; Mp is the mass per unit length of the cantilever pipe; Q is the transverse shear force in the cantilevered pipe; g is the gravitational acceleration; Fin and Fit are the normal and tangential hydrodynamic forces because of the internal flow Ui, respectively; Fen and Fet are the normal and tangential hydrodynamic forces because of the external flow, respectively.
As approximated by the Euler–Bernoulli (E–B) beam theory, one obtains the following relationship
Q = x ( E I 2 w x 2 )
where EI is the flexural rigidity. Internal dissipation in the cantilevered pipe material is relatively small compared to dissipation by the surrounding fluid so it can be neglected [13]. Substituting Equation (3) into Equations (1) and (2), and neglecting second-order terms, one finds the following expressions in the x- and z-directions, respectively
F T x ( F i t + F i n w x ) + ( F e t F e n w x ) + M p g = 0
E I 4 w x 4 x ( F T w x ) ( F i n F i t w x ) ( F e n + F e t w x ) + M p 2 w t 2 = 0
The Fin and Fit (hydrodynamic forces because of internal flow) are obtained by a force balance on an element δx of internal fluid (Figure 5b) in the manner of Païdoussis [14,15]. The Fin′ and Fit′ in Figure 5b are a pair of interaction forces with the Fin and Fit in Figure 4a, Fin′ = −Fin and Fit′ = −Fit. The resulting expressions are written as below in the x- and z-directions, respectively,
F i t + F i n w x = M f g + A f p i x
and
F i n F i t w x = M f ( 2 w t 2 2 U i 2 w x t + U i 2 2 w x 2 ) A f x ( p i w x )
where Mf is the mass of internal fluid per unit length; pi is the internal pressure in the cantilevered pipe; Af is the cross-sectional area of the internal flow, and Mf = ρf Af, in which ρf is the fluid density.
Substituting Equations (6) and (7) in (4) and (5) renders the following equations in the x- and z-directions, respectively
F T x + ( M f g A f p i x ) + ( F e t F e n w x ) + M p g = 0
E I 4 w x 4 x ( F T w x ) + [ M f ( 2 w t 2 2 U i 2 w x t + U i 2 2 w x 2 ) + A f x ( p i w x ) ] ( F e n + F e t w x ) + M p 2 w t 2 = 0
The field of flow outside the cantilevered pipe can be simplified because the cantilever pipe displacements w are considered to be small according to E–B beam theory. More specifically, the external flow is modelled as the superposition of a perturbation potential because of the cantilever pipe vibrations in the mean axial flow [15,16]. Hence, the effects of fluid viscosity are considered to be different from this potential flow, and the viscosity-related forces are added separately to the system. When linear investigations of fluid-structure interaction issues are carried out, this is a common approach. There are more details in [5,15].
The resultant force of all external hydrodynamic forces acting on the cantilevered pipe (Figure 4a) can be broken down into Fen and Fet in the directions perpendicular and tangential to the cantilever pipe centerline, respectively. Projection of the external hydrodynamic force in the x- and z-directions yields the following force-balance equations
F e t F e n w x = F L F p x ,
and
F e n + F e t w x = ( F A + F N ) + F p z + F L w x .
In Equations (10) and (11), FA is the lateral inviscid hydrodynamic force; Fpx and Fpz are the resultant forces due to external pressure in the x and z directions, respectively; FN and FL are the fluid frictional viscous forces perpendicular and tangential to the cantilever pipe centerline, respectively; and high-order terms have been omitted.
The added (hydrodynamic) mass Ma is the mass per unit length for the annular (external) fluid associated with motions of the cantilevered pipe. Based on the perturbation potential, for a thin boundary layer, the added (hydrodynamic) mass Ma = χρfAo [17,18], where Ao = πDo2/4 is the external cross-sectional area of the cantilevered pipe. The parameter χ quantifies the effect of radial confinement degree (or narrowness) of the annular flow passage as
χ = ( D c h / D o ) 2 + 1 ( D c h / D o ) 2 1 , lim D c h ( χ ) = lim D c h ( D c h / D o ) 2 + 1 ( D c h / D o ) 2 1 = 1 .
In Equation (12), Dch is the inner diameter of the passage (the inflexible external pipe) and Do is the outer diameter of the cantilever pipe as shown in Figure 1. Regarding an unconfined pipe, Dch → ∞ and as a result Ma is equal to ρfAo, i.e., the displaced mass of the fluid per unit length. Hence, considering the variation of confinement along the cantilevered pipe, the added (hydrodynamic) mass may be written as
M a = [ χ + ( 1 χ ) H ( x L s ) ] ρ f A o
where χ pertains to the confined portion of the cantilevered pipe, 0 ≤ xLs, and H(xLs) is the Heaviside step function.
As is shown in Figure 1a, both the external flow velocity Uo and the confinement degree vary along the length of the cantilevered pipe. For the sake of simplicity, it is assumed that Uo is zero along the unconfined part, Ls < xL, and immediately attains the value of Uo along the confined portion of the cantilevered pipe, 0 ≤ x < Ls. The above assumptions can also be found in the literature [1]. Thus, the lateral inviscid hydrodynamic force FA can be written as
F A = ( t + U o x U o H ( x L s ) x ) { [ χ + ( 1 χ ) H ( x L s ) ] ρ f A o ( w t + U o w x U o H ( x L s ) w x ) } ,
where the initial expression [16] and later [17] have been altered to explicate the spatial variation of both the external flow velocity and the added (hydrodynamic) mass Equation (13). After various simplifications, the formula given in Equation (14) was arrived at
F A = ( 1 χ ) ρ f A o H ( x L s ) 2 w t 2 + ρ f A o χ 2 w t 2 2 A o U o ρ f χ H ( x L s ) 2 w x t + 2 A o U o ρ f χ 2 w x t A o U o 2 ρ f χ H ( x L s ) 2 w x 2 + A o U o 2 ρ f χ 2 w x 2
The fluid frictional viscous force FL tangential to the cantilevered pipe centerline is
F L = 1 2 C f ρ f D o U o 2 [ 1 H ( x L s ) ] ,
where Cf is the frictional damping coefficient, with a value of 0.0125 to give an acceptable estimate of the fluid frictional viscous force FL [8].
Equation (16) was modified from the literature [17] to take the spatial variation of the mean external flow velocity into account. Likewise, in terms of the literature [17], the fluid frictional viscous force FN perpendicular to the cantilevered pipe centerline is
F N = 1 2 C f ρ f D o U o [ 1 H ( x L s ) ] { w t + [ 1 H ( x L s ) ] U o w x } + k w t ,
where the variation of the external flow velocity Uo has been taken into account; k is a viscous damping coefficient which only relates to lateral motions of the cantilever pipe, without average external flow Uo. Based on two-dimensional flow analysis, expressions used for k are derived in [13,18]. For a thin boundary layer, the coefficient k is calculated via the following expression
k = 2 2 S 1 + γ ¯ 3 ( 1 γ ¯ 2 ) 2 ρ f A o { Ω }
In (10), S = R{Ω}Do2/(4ν) represents the oscillatory Reynolds number (i.e., Stokes number), where ν and R{Ω} are the kinematic viscosity of the fluid and the circular frequency of oscillation, respectively; γ ¯ = Do/Dch is a measure of the radial confinement degree of the annular passage. Naturally, based on Figure 1a, the value of k varies along the cantilevered pipe. This variation of k is explained by letting
F N = 1 2 C f ρ f D o U o { [ 1 H ( x L s ) ] w t + [ 1 H ( x L s ) ] U o w x } + k u [ 1 + γ ¯ 3 ( 1 γ ¯ 2 ) 2 + H ( x L s ) ( 1 1 + γ ¯ 3 ( 1 γ ¯ 2 ) 2 ) ] w t .
k u = 2 2 S ρ f A o { Ω }
In (19) and (20), ku represents the friction coefficient applicable to the cantilevered pipe without confined portion, Ls < xL.
Finally, the forces (Fpx and Fpz) caused by mean tensioning (induced by gravity) and pressurization are considered, respectively. Fpx and Fpz are cleverly derived in [15,17] and they are expressed as
F p x = x ( A o p o ) + A o p o x
and
F p z = A o x ( p o w x )
in the x- and z-directions respectively, where po represents the fluid pressure outside the cantilevered pipe.
For the purpose of analyzing the model shown in Figure 1a, the outlet effects of the flow exiting from the annular passage have been reduced to a tinily thin region at x = Ls. For 0 ≤ x < Ls, there is a pressure loss due to the friction of the flowing fluid inside the annulus, while there is a pressure increasing due to the gravity. For Ls < xL, there is a purely hydrostatic pressure distribution along the unconfined portion of the cantilever pipe. For 0 ≤ x < Ls, an annular fluid element (i.e., the external flow of length δx) (Figure 5a) is considered, and a force balance per unit length is obtained as
A c h p o x F f + A c h ρ f g = 0 ,
where Ach = π(Dch2Do2)/4 represents the cross-sectional flow area of the annular passage; Ff represents the total wall-friction force acting on the fluid element.
It is assumed that an equal wall shear stress acts on both annular passage surfaces, the total wall-friction force Ff can be expressed as
F f S t o t = F L S o
in which Stot = π(Dch + Do) is the total wetted area of the annular flow per unit length, and So is the outside wetted pipe perimeter.
Combining Equations (23) and (24), multiplied by (Ao/Ach), the resultant equation is rearranged as
A o p o x = F L ( D o D h ) + A o ρ f g ,
where Dh = 4Ach/Stot = (DchDo) is the hydraulic diameter of the external flow in the annular passage. Integration of Equation (25), where for the confined portion external fluid pressure at x = 0 is the reference pressure (po|x = 0 = 0), with respect to x gives
p o ( x ) = [ F L A o ( D o D h ) + ρ f g ] x ,
for 0 ≤ x < Ls.
For the unconfined part, i.e., Ls < xL, the external fluid pressure distribution is hydrostatic
p o x = ρ f g ,
which upon integration results in
p o ( x ) = ρ f g x + C 1 ,
in which C1 is an integration constant determined in the following description.
It is assumed that x1 = Ls is the axial location just inside the annular passage, and x2 = Ls+ is the location just outside the annular passage. Based on Bernoulli’s equation, an energy balance of the fluid is obtained from x1 to x2
p o | x 2 = p o | x 1 + 1 2 ρ f U o 2 ρ f g h e , h e = K 1 U o 2 2
where the quantity he, with K1 = 1, is the head-loss associated with sudden enlargement of the external flow in the annular passage into the surrounding fluid at x = Ls+ [19]. Combination Equations (28) and (29) yields
C 1 = [ F L A o ( D o D h ) L s ] + 1 2 ρ f U o 2 ρ f g h e ,
in which ||FL|| = 0.5ρfCfDoUo2.
Through the above analysis, the resulting fluid pressure gradient and fluid pressure distribution over the whole cantilevered pipe, for 0 ≤ xL, are, respectively
p o x = F L A o ( D o D h ) [ 1 H ( x L s ) ] + ρ f g + ( 1 2 ρ f U o 2 ρ f g h e ) δ D ( x L s ) ,
and
p o = F L A o ( D o D h ) x + F L A o ( D o D h ) ( x L s ) H ( x L s ) + ρ f g x + ( 1 2 ρ f U o 2 ρ f g h e ) H ( x L s ) ,
where δD(xLs) is the Dirac delta function.
Substituting Equations (16) and (21) in (10), and subsequently substituting (10) into (8), the force balance on the cantilevered pipe is obtained in the x-direction, namely
x ( F T A f p i + A o p o ) + M p g + ρ f A f g A o p o x + 1 2 C f ρ f D o U o 2 [ 1 H ( x L s ) ] = 0
Substituting Equations (15), (16), (19), and (22) in (11), and subsequently substituting the result in (1), the motion equation is obtained in the z-direction,
E I 4 w x 4 x [ ( F T A f p i + A o p o ) w x ] + M p 2 w t 2 + [ M f ( 2 w t 2 2 U i 2 w x t + U i 2 2 w x 2 ) ] + ( 1 χ ) ρ f A o H ( x L s ) 2 w t 2 + ρ f A o χ 2 w t 2 2 A o U o ρ f χ H ( x L s ) 2 w x t + 2 A o U o ρ f χ 2 w x t A o U o 2 ρ f χ H ( x L s ) 2 w x 2 + A o U o 2 ρ f χ 2 w x 2 + 1 2 C f ρ f D o U o [ 1 H ( x L s ) ] w t + k u [ 1 + γ ¯ 3 ( 1 γ ¯ 2 ) 2 + H ( x L s ) ( 1 1 + γ ¯ 3 ( 1 γ ¯ 2 ) 2 ) ] w t = 0 .
In Equation (34), the tensioning and pressurization term (FTpiAf + Aopo) is determined by integrating Equation (33) from x to L, which yields
( F T A f p i + A o p o ) = 1 2 C f ρ f D o U o 2 ( D o D h + 1 ) ( L s x ) [ 1 H ( x L s ) ] A o ( 1 2 ρ f U o 2 ρ f g h e ) [ 1 H ( x L s ) ] + ( M p g + ρ f A f g A o ρ f g ) ( L x ) + ( F T A f p i + A o p o ) | L .
In Equation (35) |L represents an evaluation of the term in parentheses at x = L.
Hence, the final form of the motion equation is obtained in the z-direction by substituting Equation (35) in (34), as follows
E I 4 w x 4 + { ( M p + ρ f A f ρ f A o ) g + 1 2 C f ρ f D o U o 2 ( D o D h + 1 ) [ 1 H ( x L s ) ] A o ( 1 2 ρ f U o 2 ρ f g h e ) δ D ( x L s ) } w x + { ( M p ρ f A f + ρ f A o ) g ( L x ) 1 2 C f ρ f D o U o 2 ( D o D h + 1 ) ( L s x ) [ 1 H ( x L s ) ] + A o ( 1 2 ρ f U o 2 ρ f g h e ) [ 1 H ( x L s ) ] ( T A f p i + A o p o ) | L } 2 w x 2 + M p 2 w t 2 + ρ f A f 2 w t 2 2 U i ρ f A f 2 w x t + ρ f A f U i 2 2 w x 2 + ρ f A o χ U o 2 [ 1 H ( x L s ) ] 2 w x 2 + 2 A o U o ρ f χ [ 1 H ( x L s ) ] 2 w x t + ( 1 χ ) ρ f A o H ( x L s ) 2 w t 2 + ρ f A o χ 2 w t 2 + 1 2 C f ρ f D o U o [ 1 H ( x L s ) ] w t + k u { 1 + [ 1 H ( x L s ) ] ( 1 + γ ¯ 3 ( 1 γ ¯ 2 ) 2 1 ) } w t = 0
According to the Bernoulli equation, the pressures Po|L and Pi|L are related to the energy balance of the fluid at x = L
p i | L = p o | L 1 2 ρ f U i 2 ρ f g h a ,     h a = K 2 U i 2 2 .
In Equation (37), the quantity ha is the headloss due to the stagnant fluid entering the cantilevered pipe at x = L and acquiring internal flow velocity Ui, calculated with 0.8 ≤ K2 ≤ 0.9, independent of flow velocity [19]. Consequently, evaluation of Equation (32) at x = L yields Po|L
p o | L = 1 2 A o C f ρ f D o U o 2 ( D o D h ) L + ρ f g L + 1 2 ρ f U o 2 ρ f g h a .
Finally, the internal and external flow velocities are related through conservation of mass, i.e., UiAf = UoAch, which after some transformation is expressed as
U o = U i D i 2 D c h 2 D o 2 .

2.2. Boundary Conditions

The motion Equation (36) is subjected to the boundary conditions
w ( 0 , t ) = w x | x = 0 = 0 , 2 w x 2 | x = L = 3 w x 3 | x = L = 0 .

2.3. Dimensionless Motion Equation and Boundary Conditions

The motion Equation (36) is non-dimensionalized using the following dimensionless quantities
ξ = x L ,              τ = [ E I M p + ρ f A f + ρ f A o ] 1 2 t L 2 ,              η = w L , u i = ( ρ f A f E I ) 1 2 L U i ,              u o = ( ρ f A o E I ) 1 2 L U o ,              β o = ρ f A o M p + ρ f A f + ρ f A o , β i = ρ f A f M p + ρ f A f + ρ f A o ,              γ = ( M p + ρ f A f ρ f A o ) g L 3 E I ,              Γ = T | L L 2 E I , Π iL = p i | L A f L 2 E I ,              Π oL = p o | L A o L 2 E I ,              c f = 4 C f π , κ u = k u L 2 [ E I ( M p + ρ f A f + ρ f A o ) ] 1 2 ,              ε = L D o ,              h = D o D h , α = D i D o ,              α c h = D c h D o ,              r a n n = L s L .
Thus, the resulting dimensionless form of Equation (36) is
4 η ξ 4 + { γ + 1 2 c f ε u o 2 ( h + 1 ) [ 1 H ( ξ r a n n ) ] 1 2 u o 2 ( 1 K 1 ) δ D ( ξ r a n n ) } η ξ { γ ( 1 ξ ) + 1 2 c f ε u o 2 ( h + 1 ) ( r a n n ξ ) [ 1 H ( ξ r a n n ) ] 1 2 u o 2 ( 1 K 1 ) [ 1 H ( ξ r a n n ) ] + ( Γ Π iL + Π oL ) } 2 η ξ 2 + { u i 2 + χ u o 2 [ 1 H ( ξ r a n n ) ] } 2 η ξ 2 + { 1 + ( χ 1 ) β o [ 1 H ( ξ r a n n ) ] } 2 η τ 2 + 2 ( χ u o β o 1 2 [ 1 H ( ξ r a n n ) ] u i β i 1 2 ) 2 w ξ τ + 1 2 c f ε u o β o 1 2 [ 1 H ( ξ r a n n ) ] η τ + κ u { 1 + [ 1 H ( ξ r a n n ) ] ( 1 + α c h 3 ( 1 α c h 2 ) 2 1 ) } η τ = 0 ,
with corresponding dimensionless boundary conditions
η ( 0 , τ ) = η ξ | ξ = 0 = 0 , 2 η ξ 2 | ξ = 1 = 3 η ξ 3 | ξ = 1 = 0 .
Additionally, the dimensionless (complex) frequency of vibration, ω, is related to the dimensional one, Ω, by
ω = [ M p + ρ f A f + ρ f A o E I ] 1 2 L 2 Ω .
It is noted that the external and internal fluid pressures in dimensionless form at the free end have the relationship
Π iL = α 2 Π oL 1 2 u i 2 A f ρ f g h a ( L 2 E I ) ,
where
Π oL = 1 2 c f h r a n n ε u o 2 + 1 2 u o 2 ( 1 K 1 ) + A o ρ f g L 3 E I ,
with ha given in Equation (37), and K1 as shown in Equation (29).
Finally, the dimensionless flow velocities are expressed as
u o = α α c h 2 1 u i ,
where α and αch are defined in (41).

3. Solution of Equations by a Galerkin Method

The governing Equation (42) of lateral vibration of the cantilever pipe contains the Heaviside step function, which is a discontinuous function, and a conventional Galerkin method is used here to discretize this system into an ordinary differential equation. The cantilever pipe is divided into N units. Consequently, solution of Equation (42) was effected as follows. Let an approximate solution be
η ( ξ , τ ) η ˜ ( ξ , τ ) = j = 1 N Φ j ( ξ ) q j ( τ ) .
In (48), qj(τ) is the generalized coordinate of the system and written as
q j ( τ ) =   s j e i ω j τ = s j e Im ( ω j ) τ e i Re ( ω j ) τ .
The Φj(ξ) are the appropriate comparison functions, in this case the normalized beam eigenfunctions for the fixed-free E–B beam, and satisfy the geometric and natural boundary conditions of the problem. The qj(τ) are the generalized coordinates of the system. In Equation (49), ωj are the eigenfrequencies; sj is a complex amplitude; and Im(ωj) and Re(ωj) are the imaginary and real part of ωj, respectively. It is assumed that the solution of Equation (42) is separable in terms of the dimensionless spatial variable ξ and the dimensionless time τ.
Equation (48) denotes a truncated series, in which N is finite positive integer and it represents the number of the appropriate comparison functions. For the sake of compactness of the notation, the following integral expressions are defined as
a i j ( m , n ) m n Φ i Φ j d ξ ,     b i j ( m , n ) m n Φ i ( Φ j ξ ) d ξ , c i j ( m , n ) m n Φ i ( 2 Φ j ξ 2 ) d ξ ,     d i j ( m , n ) m n ξ Φ i ( 2 Φ j ξ 2 ) d ξ ,
where i, j are indices corresponding to the relevant quantity to be integrated; m and n of (m, n) are the lower/upper bounds of integration.
Substituting Equations (48) and (50) into (42), simultaneous equations of the form are obtained by pre-multiplying by Φi and integrating from ξ = 0 to ξ = 1
M d 2 d τ 2 ( { q 1 q 2 q 3 } ) + C d d τ ( { q 1 q 2 q 3 } ) + K ( { q 1 q 2 q 3 } ) = 0 .
In Equation (51), M, C, and K matrices are mass matrix, damping matrix and stiffness matrix, respectively. The components of the M, C, and K matrices are
M i j = a i j ( 0 , 1 ) β o ( 1 χ ) a i j ( 0 , r a n n ) ,
C i j = 2 u i β i 1 2 b i j ( 0 , 1 ) + 2 χ u o β o 1 2 b i j ( 0 , r a n n ) + 1 2 c f ε u o β o 1 2 a i j ( 0 , r a n n ) + κ u a i j ( 0 , 1 ) + κ u ( 1 + α c h 3 ( 1 α c h 2 ) 2 1 ) a i j ( 0 , r a n n ) ,
K i j = λ j 4 a i j ( 0 , 1 ) + γ b i j ( 0 , 1 ) 1 2 c f ε u o 2 ( h 1 ) b i j ( 0 , r a n n ) 1 2 u o 2 ( 1 K 1 ) ( Φ i | ξ = r a n n Φ j ξ | ξ = r a n n ) ( Γ Π iL + Π oL ) c i j ( 0 , 1 ) γ ( c i j ( 0 , 1 ) d i j ( 0 , 1 ) )     + 1 2 c f ε u o 2 ( h 1 ) ( r a n n c i j ( 0 , r a n n ) d i j ( 0 , r a n n ) ) + 1 2 u o 2 ( 1 K 1 ) c i j ( 0 , 1 ) + u i 2 c i j ( 0 , 1 ) + χ u o 2 c i j ( 0 , r a n n ) ;
where λj is the jth dimensionless eigenvalue of the fixed-free E–B beam.
Solutions of Equation (51) constitute approximate solutions of Equation (42) with boundary conditions in Equation (43). In order to obtain a non-trivial solution of Equation (51), it is required that the determinant of coefficient matrix is equal to zero. This corresponds to a generalized eigenvalue problem, and the stability of the system can be determined by the generalized eigenvalue (i.e., eigenfrequency) ω of matrix E which consists of M, C, and K, as
E = [ M 1 C M 1 K I 0 ] .
For a given internal dimensionless flow velocity ui, the stability of the cantilever pipe is determined by Im(ω). When Im(ω) < 0, the qj(τ) (shown in Equation (49)) grows exponentially in time, so Im(ω) < 0 indicates instability. Consequently, for Im(ω) < 0 the cantilevered pipe becomes unstable, and for Im(ω) = 0, divergence or buckling of the cantilever pipe happens. Divergence represents a static loss of stability, vulgarly known as buckling and, in the nonlinear dynamics milieu, as a static pitchfork bifurcation. In the following passages, any discussion of the system stability involves only the stability of the internal cantilevered pipe.

4. Theoretical Analysis for Slender, Leaching-Tubing-Like Systems

The geometry shown in Figure 1 has been applied in the salt cavern underground energy storage industry, called leaching-tubing systems. Salt caverns serving as underground energy storage are generally located at a depth between 500 and 2000 m [3], and storage volumes of salt caverns range from 5 × 104 m3 to 2 × 105 m3 [20,21]. The typical depths of a salt cavern height and the salt cavern ceiling are about 400 m and 600 m, respectively, learned from the industry survey described in [2]. Hence common lengths of the leaching tubing are on the order of one kilometer.
Based on information from industrial applications, theoretical analyses for slender leaching-tubing-like systems are carried out in this section. Sample dimensional and associated dimensionless parameters for the leaching tubing are given in Table 1. In the salt cavern storage industry, the geometric parameters αch, rann and L shown in Figure 1 can be varied by operators, as required by the engineering practice. Therefore, it is physically meaningful to discuss the effects of the geometric parameters αch, rann, and L on the pipe system. The effects of varying αch, rann and L are discussed from Section 4.1 to Section 4.3.

4.1. Effect of the Radial Confinement αch

Here, it is interesting to see the effect of different radial confinement values αch on the leaching-tubing-like system behavior. According to Equation (47), decreasing αch corresponds to an increase in radial confinement and to a related increase in the dimensionless annular flow velocity uo.
The leaching-tubing-like systems behavior with the variation of αch from 1.10 to 20 is summarized in Table 2, based on rann = 0.50, and ε = 1124.9 (L = 200 m) and the remaining parameters shown in Table 1. Theoretical results of both ucr and ωcr are shown in Figure 6. ucr is a dimensionless critical instability flow velocity and indicates the onset of instability of the internal pipe. For example, Figure 7 shows the vibration and stability of the system with αch = 1.20, for the first three modes. From Figure 7, we can see that when ucr = 4.34 the first predicted instability of the system occurs, and the associated frequency ωcr = 0, which indicates the onset of divergence. Naturally, these results presented here are related to this particular parameters selection, and only indicate what the leaching-tubing-like systems behavior could be.
As shown in Table 2, when αch ≈ 20, i.e., the leaching-tubing-like system is effectively unconfined, flutter instability occurs in the second mode. For 1.35 ≤ αch < 20, the decrease in αch results in flutter instability in which the modes are sequentially lowered. However, for 1.32 ≤ αch < 1.35, flutter instability occurs in the third mode. A further decrease of αch below αch < 1.32 can still result in flutter instability in which the modes are sequentially lowered, until first-mode divergence (i.e., a static instability) arises for 1.10 ≤ αch < 1.27. The vanishing of Re(ω) indicates the beginning of divergence, as clearly seen in Figure 6.
With regard to the dimensionless critical flow velocity ucr, and with the very confined case (αch = 1.10) as a reference, the increase of αch leads initially to the rapid increase of ucr. This is followed by a steep and significant decrease of ucr as αch is increased continuously. Moreover, when αch > 1.48, ucr is less than 1.00. However, with a further increase of αch above αch > 2.00 approximately, ucr remains almost unaffected. As mentioned above, when 1.10 ≤ αch < 1.27, a first-mode divergence arises, so ωcr is equal to 0. An increase of αch to 1.27 ≤ αch < 1.35 causes a significant increase of ωcr. However, when αch > 1.35, decreased confinement causes a very steep reduction in ωcr. When 1.35 ≤ αch < 4.46, ωcr remains almost unaffected with the increase of αch. When 4.46 ≤ αch < 20, the system instability state changes from first-order flutter to second-order flutter, ωcr steeply increases and remains almost unaffected with the increase of αch.

4.2. Effect of the Confinement Length rann

The effect of the confinement length rann (i.e., confined-flow length) on the slender leaching-tubing-like system is studied next. A larger value of the confinement length rann corresponds to a larger part of the cantilevered pipe being subjected to annular flow. In the schematic view of the theoretical model shown in Figure 1, the inflexible external pipes cannot be longer than the internal cantilevered pipe. Consequently, the studied values of rann are less than 1 and the studied maximum value of rann is 0.9.
The leaching-tubing-like systems behavior with the variation of rann from 0 to 0.90 is summarized in Table 3, respectively based on αch = 1.676, ε = 562.5 (L = 100 m) and αch = 1.676, ε = 1124.9 (L = 200 m) and the remaining parameters shown in Table 1. Theoretical results of both ucr and ωcr are shown in Figure 8. The value of αch = 1.676 is used, due to the typical degree of confinement in the real application. Here, two different leaching-tubing lengths (i.e., L = 100 m and L = 200 m) are used to show the qualitative difference in system behavior.
The slender leaching-tubing-like system behavior is discussed by considering rann values increased from rann = 0. For the shorter, (L = 100 m) leaching-tubing-like system initially loses stability by the first-mode flutter. A further increase of rann results in the second-mode flutter. For the longer, (L = 200 m) leaching-tubing-like system, an increasing rann leads initially to loss of stability by the second-mode flutter. A further increase of rann leads to first-mode flutter and second-mode flutter.
For the shorter leaching-tubing length considered, in terms of the dimensionless critical flow velocity ucr, an increase of rann from 0 to 0.50 results in a slight stabilization. ucr increases nonlinearly with a further increase of rann from 0.50 to 0.72. When rann increases from 0.72 to 0.90, ucr begins to decrease and reaches a minimum at rann = 0.78, after which ucr continues to increase. As shown in Figure 8a, ωcr is almost unaffected by rann during the first-mode flutter phase and the second-mode flutter phase. However, when the unstable state changes from the first-mode flutter phase to the second-mode flutter phase, ωcr increases steeply. For the longer leaching-tubing length considered, the dimensionless critical flow velocity ucr remains almost unaffected with the increase of rann from 0 to 0.30. Then, ucr increases nonlinearly with a further increase of rann from 0.30 to 0.69. When rann increases from 0.69 to 0.90, ucr begins to decrease and reaches a minimum at rann = 0.81, after which ucr continues to increase. As shown in Figure 8b, ωcr is almost unaffected by rann during the first-mode flutter phase and the second-mode flutter phase. However, when the unstable state changes from the second-mode flutter phase to the first-mode flutter phase, ωcr steeply decreases. When the unstable state returns from the first-mode flutter back to the second-mode flutter, ωcr steeply increases. Hence, for both leaching-tubing lengths considered, the curve of ucr has a U shape, and the curve of ωcr is stepped.

4.3. Effect of the Cantilevered Pipe Length L

Dimensional length L of the cantilevered pipe has been used in the nondimensionalization of both the flow velocity ui and the frequency ω. Consequently, the dimensionless parameters given in Equation (41) are unsuitable for studying the dependence of the dynamics on the system length. Another set of dimensionless parameters have been defined in [9] to study the dynamic characteristics of slender hanging pipes conveying fluid and in [22] to study the dynamic characteristics of cylinders in axial flow. In the present paper, no new dimensionless quantities have been defined, but rather the dimensional critical flow velocity Ucr and associated frequency Ωcr have been used to carried out a brief analysis. Notwithstanding the loss of generality, this approach facilitates illustrating the variation of the Ucr and Ωcr with increasing the pipe system length L. Based on diameters, mass, and stiffness typical of a leaching-tubing, as shown in Table 1, numerical results are obtained to predict dynamical behaviors of the leaching-tubing-like systems with increasing lengths L for various values of rann.
As shown in Table 4, the leaching-tubing-like systems behavior with the variation of L from 5 m to 500 m is summarized, based on rann = 0.125, 0.25, 0.50, and 0.85. When 5 m ≤ L ≤ 500 m, with increasing of the length, the leaching-tubing loses stability by first-mode flutter, for rann = 0.125, 0.25, and 0.50, while the leaching-tubing loses stability by second-mode flutter, for rann = 0.85. Theoretical results of both Ucr and Ωcr are shown in Figure 9 and Figure 10. When the length L is below 40 m, an increasing of L leads initially to a nearly negligible stabilization. After that, Ucr increases significantly with a further increase of L from 40 m to 500 m. Ucr increases approximately linearly, for rann = 0.125, 0.25, and 0.50, while Ucr increases approximately nonlinearly, for rann = 0.85. For all studied values of rann, with increasing the pipe system length L, the critical frequency Ωcr decrease. Especially when L < 50 m, Ωcr decreases sharply. Moreover, for a long enough system, the frequency Ωcr almost vanishes, which indicates a divergence is about to begin.

5. Conclusions

A linear theoretical model has been formulated for a hanging vertical cantilevered pipe which is subjected concurrently to two dependent axial flows. It should be pointed out that after the external flow exits from the outside annular passage, fluid is conveyed upwards as the internal flow in the cantilever pipe. The motion equation of the system is solved by means of a Galerkin method, and eigenfunctions of Euler–Bernoulli beam are used as comparison functions. Theoretical predictions were obtained for a long leaching-tubing-like system with parameters related to the salt cavern energy storage.
A theoretical study into the effect of radial confinement (i.e., the degree of radial confinement of the annular passage) has shown that, the system loses stability in flutter and even in divergence regardless of length, as the radial confinement degree and consequently the annular flow velocity is increased. When αch is greater than a certain value, ucr decreases as αch increases, even to a small value. Considering that the flow velocity in industrial applications is relatively large, the value of αch should not be too large. Additionally, only when values of αch exceed a certain parameter-dependent threshold does the radial confinement has no appreciable effect on the pipe system.
The confinement length rann (i.e., confined-flow length), was shown to have a destabilizing or stabilizing effect on the pipe system, depending on value ranges of rann. The effect is a clear and significant stabilization, especially when the annular flow extends over most portion of the cantilevered pipe.
Furthermore, investigations have been conducted in dimensional terms on the effect of the pipe system length L, for the leaching-tubing-like system. The pipe system length L was shown to have a stabilizing effect on the pipe system. Increasing the pipe system length L results in an increase of dimensional critical flow velocity Ucr, while it results in a decrease of the associated dimensional frequency Ωcr. For a long enough system, the frequency Ωcr almost vanishes, which indicates a divergence is about to begin.
The next step is to conduct experiments to check the validity of the linear theoretical model, based on the device described in [23], and further research results are deferred to another paper.

Author Contributions

X.G. conceived the study, derived formulas, interpreted the results and wrote the paper under the guidance of Y.L. and C.Y. X.C. and N.Z. helped in the numerical calculation. X.S. and H.M. reviewed the study plan, edited the manuscript, and corrected the grammar mistakes. H.Y. participated in data analysis.

Funding

This research was funded by National Natural Science Foundation of China grant number 51874273, 51874274, 51774266 and by National Key Research and Development Program of China grant number 2018YFC0808401. And The APC was funded by 51874274.

Acknowledgments

The authors would like to be sincerely grateful to Jaak J.K. Daemen from University of Nevada, USA, for his English help and thoughtful review of this paper. The authors would like to gratefully acknowledge the financial support from National Natural Science Foundation of China (nos. 51874273, 51874274, 51774266), and the National Key Research and Development Program of China (grant no. 2018YFC0808401).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Moditis, K.; Païdoussis, M.P.; Ratigan, J. Dynamics of a partially confined, discharging, cantilever pipe with reverse external flow. J. Fluids Struct. 2016, 63, 120–139. [Google Scholar] [CrossRef]
  2. Ratigan, J.L. Brine string integrity and model evaluation. In Proceedings of the Processes of SMRI Fall Meeting, Galveston, TX, USA, 13–14 October 2008; pp. 273–293. [Google Scholar]
  3. Li, Y.; Yang, C.; Qu, D.; Yang, C.; Shi, X. Preliminary study of dynamic characteristics of tubing string for solution mining of oil/gas storage salt caverns. Rock Soil Mech. 2012, 33, 681–686. [Google Scholar]
  4. Cesari, F.; Curioni, S. Buckling instability in tubes subject to internal and external axial fluid flow. In Proceedings of the 4th Conference on Dimensioning, Hungarian Academy of Science, Budapest, Hungary, October 1971; pp. 301–311. [Google Scholar]
  5. Hannoyer, M.; Païdoussis, M.P. Instabilities of tubular beams simultaneously subjected to internal and external axial flows. J. Mech. Des. 1978, 100, 328–336. [Google Scholar] [CrossRef]
  6. Paidoussis, M.P.; Besancon, P. Dynamics of arrays of cylinders with internal and external axial flow. J. Sound Vib. 1981, 76, 361–379. [Google Scholar] [CrossRef]
  7. Wang, X.; Bloom, F. Dynamics of a submerged and inclined concentric pipe system with internal and external flows. J. Fluids Struct. 1999, 13, 443–460. [Google Scholar] [CrossRef]
  8. Païdoussis, M.P.; Luu, T.; Prabhakar, S. Dynamics of a long tubular cantilever conveying fluid downwards, which then flows upwards around the cantilever as a confined annular flow. J. Fluids Struct. 2008, 24, 111–128. [Google Scholar] [CrossRef]
  9. Doaré, O.; De Langre, E. The flow-induced instability of long hanging pipes. Eur. J. Mech. A. Solids 2002, 21, 857–867. [Google Scholar] [CrossRef]
  10. Lemaitre, C.; Hémon, P.; De Langre, E. Instability of a long ribbon hanging in axial air flow. J. Fluids Struct. 2005, 20, 913–925. [Google Scholar] [CrossRef]
  11. Kuiper, G.; Metrikine, A. Dynamic stability of a submerged, free-hanging riser conveying fluid. J. Sound Vib. 2005, 280, 1051–1065. [Google Scholar] [CrossRef]
  12. Kuiper, G.; Metrikine, A. Experimental investigation of dynamic stability of a cantilever pipe aspirating fluid. J. Fluids Struct. 2008, 24, 541–558. [Google Scholar] [CrossRef]
  13. Chen, S.; Wambsganss, M.T.; Jendrzejczyk, J. Added mass and damping of a vibrating rod in confined viscous fluids. J. Appl. Mech. 1976, 43, 325–329. [Google Scholar] [CrossRef]
  14. Païdoussis, M.P. Fluid-Structure Interactions: Slender Structures and Axial Flow, 2nd ed.; Academic Press: Oxford, UK, 2014; Volume 1. [Google Scholar]
  15. Païdoussis, M.P. Fluid-Structure Interactions: Slender Structures and Axial Flow, 2nd ed.; Academic Press: Oxford, UK, 2016; Volume 2. [Google Scholar]
  16. Lighthill, M. Note on the swimming of slender fish. J. Fluid Mech. 1960, 9, 305–317. [Google Scholar] [CrossRef]
  17. Paidoussis, M.P. Dynamics of cylindrical structures subjected to axial flow. J. Sound Vib. 1973, 29, 365–385. [Google Scholar] [CrossRef]
  18. Sinyavskii, V.; Fedotovskii, V.; Kukhtin, A. Oscillation of a cylinder in a viscous liquid. Sov. Appl. Mech. 1980, 16, 46–50. [Google Scholar] [CrossRef]
  19. Brater, E.F.; King, H.W.; Lindell, J.E.; Wei, C.Y. Handbook of Hydraulics for the Solution of Hydraulic Engineering Problems; McGraw-Hill: Boston, MA, USA, 1996. [Google Scholar]
  20. Liu, W.; Muhammad, N.; Chen, J.; Spiers, C.; Peach, C.; Deyi, J.; Li, Y. Investigation on the permeability characteristics of bedded salt rocks and the tightness of natural gas caverns in such formations. J. Nat. Gas Sci. Eng. 2016, 35, 468–482. [Google Scholar] [CrossRef]
  21. Liu, W.; Chen, J.; Jiang, D.; Shi, X.; Li, Y.; Daemen, J.K.; Yang, C. Tightness and suitability evaluation of abandoned salt caverns served as hydrocarbon energies storage under adverse geological conditions (AGC). Appl. Energy 2016, 178, 703–720. [Google Scholar]
  22. De Langre, E.; Paidoussis, M.; Doaré, O.; Modarres-Sadeghi, Y. Flutter of long flexible cylinders in axial flow. J. Fluid Mech. 2007, 571, 371–389. [Google Scholar] [CrossRef]
  23. Ge, X.; Li, Y.; Shi, X.; Chen, X.; Ma, H.; Yang, C.; Shu, C.; Liu, Y. Experimental device for the study of liquid–solid coupled flutter instability of salt cavern leaching tubing. J. Nat. Gas Sci. Eng. 2019, in press. [Google Scholar] [CrossRef]
Figure 1. Schematic view of the pipe configuration. A hanging flexible cantilevered pipe conveying fluid is partially confined inside a shorter inflexible external pipe. (a) the system considered in this paper; (b) the system considered in [1].
Figure 1. Schematic view of the pipe configuration. A hanging flexible cantilevered pipe conveying fluid is partially confined inside a shorter inflexible external pipe. (a) the system considered in this paper; (b) the system considered in [1].
Applsci 09 01425 g001
Figure 2. Schematic view of the salt cavern underground energy storage application. Fresh water or light brine is injected into the inner tubing, and high-concentration brine is expelled through the annular space in between the inner and outer tubings. This process is known as direct circulation, as distinguished from reverse circulation. Direct and reverse circulations are alternately employed until the salt cavern achieves the design shape and dimensions.
Figure 2. Schematic view of the salt cavern underground energy storage application. Fresh water or light brine is injected into the inner tubing, and high-concentration brine is expelled through the annular space in between the inner and outer tubings. This process is known as direct circulation, as distinguished from reverse circulation. Direct and reverse circulations are alternately employed until the salt cavern achieves the design shape and dimensions.
Applsci 09 01425 g002
Figure 3. Excessive bending of the inner tubing at a Chinese salt cavern storage [3].
Figure 3. Excessive bending of the inner tubing at a Chinese salt cavern storage [3].
Applsci 09 01425 g003
Figure 4. (a) Forces acting on an element of length δx of the cantilevered pipe. (b) Forces due to the outside fluid acting on an element δx of the cantilevered pipe.
Figure 4. (a) Forces acting on an element of length δx of the cantilevered pipe. (b) Forces due to the outside fluid acting on an element δx of the cantilevered pipe.
Applsci 09 01425 g004
Figure 5. (a) Forces acting on an element of length δx of the annular flow. (b) Forces acting on an element δx of the internal flow.
Figure 5. (a) Forces acting on an element of length δx of the annular flow. (b) Forces acting on an element δx of the internal flow.
Applsci 09 01425 g005
Figure 6. Leaching-tubing-like systems with rann = 0.50 and ε = 1124.9 (L = 200 m), theoretical variation of ucr and ωcr with the radial confinement αch. The results are based on the instability of the first prediction. A zero frequency (i.e., ωcr = 0) indicates divergence.
Figure 6. Leaching-tubing-like systems with rann = 0.50 and ε = 1124.9 (L = 200 m), theoretical variation of ucr and ωcr with the radial confinement αch. The results are based on the instability of the first prediction. A zero frequency (i.e., ωcr = 0) indicates divergence.
Applsci 09 01425 g006
Figure 7. Argand diagram of the dimensionless complex eigenfrequencies of the leaching-tubing-like system, as a function of the nondimensional internal flow velocity ui for αch = 1.20. Blue solid circles represent the nondimensional internal flow velocity ui. Red solid triangles represent the critical flow velocity when flutter instability occurs in the corresponding mode.
Figure 7. Argand diagram of the dimensionless complex eigenfrequencies of the leaching-tubing-like system, as a function of the nondimensional internal flow velocity ui for αch = 1.20. Blue solid circles represent the nondimensional internal flow velocity ui. Red solid triangles represent the critical flow velocity when flutter instability occurs in the corresponding mode.
Applsci 09 01425 g007
Figure 8. Leaching-tubing-like systems, theoretical variation of ucr and ωcr with the confinement length rann. The results are based on the instability of the first prediction. (a) Leaching-tubing with αch = 1.676 and ε = 562.5 (L = 100 m). (b) Leaching-tubing with αch = 1.676 and ε = 1124.9 (L = 200 m).
Figure 8. Leaching-tubing-like systems, theoretical variation of ucr and ωcr with the confinement length rann. The results are based on the instability of the first prediction. (a) Leaching-tubing with αch = 1.676 and ε = 562.5 (L = 100 m). (b) Leaching-tubing with αch = 1.676 and ε = 1124.9 (L = 200 m).
Applsci 09 01425 g008
Figure 9. Leaching-tubing-like systems, theoretical variation of Ucr and Ωcr with the increase of the pipe system length L. The results are based on the instability of the first prediction. (a) Leaching-tubing with αch = 1.676 and rann = 0.125. (b) Leaching-tubing with αch = 1.676 and rann = 0.25.
Figure 9. Leaching-tubing-like systems, theoretical variation of Ucr and Ωcr with the increase of the pipe system length L. The results are based on the instability of the first prediction. (a) Leaching-tubing with αch = 1.676 and rann = 0.125. (b) Leaching-tubing with αch = 1.676 and rann = 0.25.
Applsci 09 01425 g009
Figure 10. Leaching-tubing-like systems, theoretical variation of Ucr and Ωcr with the increase of the pipe system length L. The results are based on the instability of the first prediction. (a) Leaching-tubing with αch = 1.676 and rann = 0.50. (b) Leaching-tubing with αch = 1.676 and rann = 0.85.
Figure 10. Leaching-tubing-like systems, theoretical variation of Ucr and Ωcr with the increase of the pipe system length L. The results are based on the instability of the first prediction. (a) Leaching-tubing with αch = 1.676 and rann = 0.50. (b) Leaching-tubing with αch = 1.676 and rann = 0.85.
Applsci 09 01425 g010
Table 1. Sample properties of real-life leaching tubings and associated dimensionless parameters [1].
Table 1. Sample properties of real-life leaching tubings and associated dimensionless parameters [1].
Dimensional
parameters
Di (m)Do (m)Dch (m)L (m)Ls (m)EI (N·m2)Mp (kg/m)
0.1590.17780.298128310853.47 × 10638.7
Dimensionless
parameters
ααchεβiβohγrann
0.8971.67672160.2390.2971.4792.019 × 1050.85
Table 2. Summary of the leaching-tubing-like systems behavior with the vibration of αch from 1.10 to 20.
Table 2. Summary of the leaching-tubing-like systems behavior with the vibration of αch from 1.10 to 20.
Behavior a
D1F2F3F1F2
αch1.101.201.271.321.354.4620
F and D represent two kinds of unstable modes, respectively, i.e., flutter instability and divergence. As an example, D1 denotes first-mode divergence and F1 denotes first-mode flutter. a [rann = 0.50, ε = 1124.9 (L = 200 m)].
Table 3. Summary of the leaching-tubing-like systems behavior with the vibration of rann from 0 to 0.90.
Table 3. Summary of the leaching-tubing-like systems behavior with the vibration of rann from 0 to 0.90.
Behavior a
F1F2
rann00.150.300.450.600.720.750.800.850.90
Behavior b
F2F1F2
rann00.150.190.450.600.690.750.800.850.90
F represents a kind of unstable mode, i.e., flutter instability. As an example, F1 denotes first-mode flutter. a [αch = 1.676, ε = 562.5 (L = 100 m)]. b [αch = 1.676, ε = 1124.9 (L = 200 m)].
Table 4. Summary of the leaching-tubing-like systems behavior with increasing L from 5 to 500 m.
Table 4. Summary of the leaching-tubing-like systems behavior with increasing L from 5 to 500 m.
Behavior a, b, c
F1
L (m)5150300450500
Behavior dF2
L (m)5150300450500
F represents a kind of unstable mode, i.e., flutter instability. As an example, F1 denotes first-mode flutter. a [αch = 1.676, rann = 0.125], b [αch = 1.676, rann = 0.25]. c [αch = 1.676, rann = 0.50], d [αch = 1.676, rann = 0.85].

Share and Cite

MDPI and ACS Style

Ge, X.; Li, Y.; Chen, X.; Shi, X.; Ma, H.; Yin, H.; Zhang, N.; Yang, C. Dynamics of a Partially Confined, Vertical Upward-Fluid-Conveying, Slender Cantilever Pipe with Reverse External Flow. Appl. Sci. 2019, 9, 1425. https://doi.org/10.3390/app9071425

AMA Style

Ge X, Li Y, Chen X, Shi X, Ma H, Yin H, Zhang N, Yang C. Dynamics of a Partially Confined, Vertical Upward-Fluid-Conveying, Slender Cantilever Pipe with Reverse External Flow. Applied Sciences. 2019; 9(7):1425. https://doi.org/10.3390/app9071425

Chicago/Turabian Style

Ge, Xinbo, Yinping Li, Xiangsheng Chen, Xilin Shi, Hongling Ma, Hongwu Yin, Nan Zhang, and Chunhe Yang. 2019. "Dynamics of a Partially Confined, Vertical Upward-Fluid-Conveying, Slender Cantilever Pipe with Reverse External Flow" Applied Sciences 9, no. 7: 1425. https://doi.org/10.3390/app9071425

APA Style

Ge, X., Li, Y., Chen, X., Shi, X., Ma, H., Yin, H., Zhang, N., & Yang, C. (2019). Dynamics of a Partially Confined, Vertical Upward-Fluid-Conveying, Slender Cantilever Pipe with Reverse External Flow. Applied Sciences, 9(7), 1425. https://doi.org/10.3390/app9071425

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