Next Article in Journal
Diauxie Studies in Biogas Production from Gelatin and Adaptation of the Modified Gompertz Model: Two-Phase Gompertz Model
Previous Article in Journal
Predicting Implicit User Preferences with Multimodal Feature Fusion for Similar User Recommendation in Social Media
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of Coriolis Force on Vibration of Annulus Pipe

by
Gian Maria Santi
*,
Daniela Francia
and
Francesco Cesari
Department of Industrial Engineering—DIN, University of Bologna, Viale del Risorgimento 2, 40100 Bologna, Italy
*
Author to whom correspondence should be addressed.
Appl. Sci. 2021, 11(3), 1058; https://doi.org/10.3390/app11031058
Submission received: 23 December 2020 / Revised: 19 January 2021 / Accepted: 21 January 2021 / Published: 25 January 2021
(This article belongs to the Section Mechanical Engineering)

Abstract

:
Annulus pipe conveying fluids have many practical applications, such as hydraulic control lines and aircraft fuel lines. In some applications, these tubes are exposed to high speeds. Normally, this leads to a vibration effect which may be of a catastrophic nature. The phenomenon is not only driven by the centrifugal forces, but an important role is played also by the Coriolis forces. Many theoretical approaches exist for a simple configuration or a complex three-dimensional configuration. Finite element models are tested. This paper provides a numerical technique for solving the dynamics of annulus pipe conveying fluid by means of the mono-dimensional Finite Element Method (FEM). In particular, this paper presents a numerical solution to the equations governing a fluid conveying pipeline segment, where a Coriolis force effect is taken into consideration both for fix and hinge constraint.

1. Introduction and Background

Fluid conveying pipes [1] are attractive form an engineering point of view because the instability effects can be a dangerous structural weakness. These systems interact with axial flows in flexible conduits with prevalent applications in the oil and gas production industries. Examples of these practical applications can be found in [2]. However, the examples found in the literature are not only limited to the field of engineering but cut across other areas of human endeavor, such as the study of pulmonary and urinary tract systems or even hemodynamics within human physiology. A lot of research interest is now focused on models for studying the stability of certain classes of dynamical systems, finding novel numerical and analytical methods for solving such problems. Thus, for linear dynamics for axial flows along slender structures, the pipe conveying fluid is regarded as the main paradigm. Direct applications in technical fields can be found in the behavior of aspirating pipes for ocean mining and Liquefied Natural Gas (LNG) in situ production, as well as for the offshore mining of methane liquid-crystal deposits and carbon sequestration. Flow-induced vibrations and instabilities are seeking to be resolved by simple, fast, and robust models.
Many researchers developed different methodologies to study the problem [3,4,5]. In particular, the Coriolis force became important in the analysis of stability, due to its role in the energy equation. Studies can be found in [6,7,8]. Helped by the experiments described in [9,10] it was possible to assume that the contribution of both the Coriolis force and the centrifugal force lead to flutter effects that can be dangerous for the strength of the structure. In [11], Eriksson et al. investigated the impact of the Coriolis force on the long distance wake behind wind farms by using Large Eddy Simulations (LES) combined with a Forced Boundary Layer (FBL) technique. The results indicate that FBL can be used for studies of long distance wakes without including a Coriolis correction, but efforts need to be taken to use a wind shear with a correct mean wind veer. The subject of piping vibration has attracted a lot of attention from various researchers in recent times, due to vast applications. In [12] Shaik et al. described a vibration analysis and mathematical model using Euler–Bernoulli and Hamilton’s energy expressions for fluid conveying a welded galvanized iron pipe with a clamped-clamped boundary condition. Many other attempts concerning the fluid conveying pipe vibrations were made and different analytical methods were applied and compared to numerical solutions by software. Avinash B. Kokare et al. [13], studied the vibrational characteristics of pipe conveying fluid and FE simulation to evaluate velocity and pressure distribution in a single phase fluid flow. Gongfa Li et al. [14] obtained the element standard equation for natural frequencies by the Lagrangian interpolation function, the first order Hermite interpolation function, and the Ritz method. Wentao Xiao et al. [15] used the Lagrangian interpolation function, the first order Hermit interpolation function, and the Ritz method to obtain the element standard equation for the nonlinear vibration response and then integrated a global matrix equation, obtaining the response of a conveying fluid pipe with the New Mark method and Matlab. Muhsin J. Jweeg et al. [16] compared results obtained by experimental verification with the aid of Smart Materials with the results performed by using analytical solution for equation of motion and also, with the results performed by using ANSYS Software. The field of nanostructures is currently turning to vibration issues and their formulation as well. Carbon nanotubes (CNTs) have been extensively used in numerous areas, due to their excellent mechanical properties. As an important mechanical property, the vibration characteristics of a CNT have become another hot topic for research in recent years. In [17], the spectral element method (SEM) is applied to analyze the dynamic characteristics of fluid conveying single-walled carbon nanotubes (SWCNTs).

2. Objectives

This paper proposes a solution to the pipe conveying fluid with Coriolis forces using a simple one-dimensional Finite Element Method with a Hermitian shape function. In particular, the modal frequencies and the critical frequencies are analyzed by changing the velocity field of the fluid. An imaginary solution is taken into account for the problem, as shown in the following. Mono-dimensional elements make possible the solution of simple numerical problems for a complex phenomena, avoiding large Fluid Structure Simulations (FSI) with difficult boundary conditions and numerical issues. In [18], an extensive literature review on the state-of-the-art numerical models in 1D, the Fluid-Structure Interaction is proposed. In the presented paper, two systems are studied: a double fixed pipe and a hinged pipe to its boundary. The first system is mostly theoretical, while the hinged pipe is more representative of the true environment that can occur in a real application. Engineering solutions are presented for the latter problem, but many other fields can be researched, such as pulmonary and urinary tract systems in the human endeavor area, highlighting the importance of simple and fast numerical methods for the analysis.

3. Methodology

In the modeling of the mechanics of fluid conveying pipes, the Coriolis force was assigned to have the role of energy absorption that counters the centrifugal effect that normally arises in free motions, affecting the stability in conservative and nonconservative systems, as reported in [1,6,7,8,19,20,21,22,23,24,25,26]. For example, in [27], sewer condition prediction models are developed to provide a framework to forecast future conditions of pipes and to schedule inspection frequencies. Furthermore, in [28], the finite-element method was effectively used to model the soil–pipe interaction for five full-scale laboratory tests conducted on a steel pipe. Such models can be used for the analysis of flexible pipe embedment design for layered embedment conditions.

3.1. Analytical Approach

From the studies reported in [29,30,31], the analytical expression for the motion equation can be written. Let us start by assuming that the curvature of a pipe in which a fluid flows with a velocity v depends on the transversal displacement w caused by its weight, the external forces and the internal forces. Let us assume that m is the total mass of fluid and pipe and ρ is the mass of fluid per length span. If the pipe is immersed in the fluid, it is possible to sum the mass of the external fluid to the mass of the pipe, and it can be calculated by multiplying the density of the fluid itself times the area of the considered section. The forces acting on the element d x of the pipe are:
  • m w ¨ = inertia force due to the vertical acceleration of the pipe;
  • E J w = elastic forces due to Euler–Bernoulli’s theory.
E is the Young module of the material, J is the inertia moment of the geometry, and the derivatives over time and space of the displacement w are defined as follows:
w ˙ = d w d t , w = d w d x , w = d 2 w d x 2 , ω = w ˙ = d 2 w d t d x
When a fluid is flowing in a pipe, and it is subjected to the Coriolis acceleration [32] through the mechanical introduction of the apparent rotation into the pipe, the amount of deflecting force generated by the Coriolis inertial effect will be a function of the mass flow rate of the fluid. If a pipe is rotated around a point while a liquid is flowing through it, that fluid will generate an inertial force acting on the pipe that will be at right angles to the direction of the flow.
Figure 1 shows the entire methodology, and Figure 2 shows a particle d m traveling at a velocity v inside a tube. The tube is rotating around a fixed point P, and the particle is at a distance of one radius r from the fixed point. The particle moves with an angular velocity ω under two components of acceleration, a centripetal acceleration directed toward P, and a Coriolis acceleration acting at right angles
  • a r (centripetal) = ω 2 r
  • a t (Coriolis) = 2 ω v
Moreover, when the pipe deflects, a centripetal force acts on an element d x of mass d m that tends to move the pipe back to the equilibrium position. If R is the curvature radius, the result is:
d F = d m v 2 R d m v 2 w d x
Substituting d m = ρ d x
d F = ρ v 2 w d x
The sign ‘-’ appears because the force d F with positive curvature points to the opposite direction with respect to the displacement w. Because of that, the forces imposed by the fluid to the pipe per unit length neglecting a r are:
F t = 2 ρ v w ˙
F = ρ v 2 w
and the equation of motion becomes:
E J 4 w x 4 + m 2 w t 2 + ρ v 2 2 w x 2 + 2 ρ v 2 w t x = 0
Before reducing Equation (1) to its dimensional form, it is possible to write down the integration form that represents the potential energy equation without considering Neumann boundary conditions
Π = 1 2 0 L E J 2 w x 2 2 + ρ v 2 2 w x 2 + 2 ρ v 2 w t x + m 2 w t 2 w d x + W
where W is the work of external loads. Interpolating the displacement field using the Hermitian shape function, separating the variables
w ( x , t ) = Δ T ( t ) N ( x ) Δ ( t ) = ( w 1 , ϕ 1 , w 2 , ϕ 2 )
where Δ ( t ) is the vector of displacements and rotation referred to each node of the finite element. In particular, the considered element is a mono-dimensional element with two nodes where the subscripts 1 and 2 indicate the quantities related to the first or the second node. Applying the principle of virtual displacement [33]:
δ Π = 0
substituting Equation (3) in Equation (2) and imposing Equation (4) lead to
δ Δ T 0 L m Δ ¨ T N N T + E J ( N ) T + 2 ρ v ( Δ ˙ T N ) N T + ρ v 2 Δ T N N T d x = δ Δ T R
where δ indicates the variation of the quantities described in the vector Δ . Reordering the members in order to form matrices result in:
M Δ ¨ + C Δ ˙ + E J ( K K n l ) Δ = R
M, K, C, and K n l are matrices defined in the following section and R is the residual matrix.

3.2. FEM Approach

In Equation (6) M, K, C, and K n l are defined as:
M = 0 L m N N T d x , K = 0 L E J ( N ) T d x , C = 0 L 2 ρ v N N T d x , K n l = 0 L ρ v 2 N N T d x
In order to use the FEM theory, it is important to integrate all the matrices in the local reference system, where the shape functions are defined as follow:
N 1 = 1 4 ( 2 3 ξ + ξ 3 ) N 2 = L 8 ( 1 ξ ξ 2 + ξ 3 ) N 3 = 1 4 ( 2 + 3 ξ ξ 3 ) N 4 = L 8 ( 1 ξ + ξ 2 + ξ 3 )

3.2.1. Fixed Pipe

Let us now consider a pipe fixed to its end nodes and defined by two finite elements. These two elements have the same length L and the total amount of degrees of freedom is two, as shown in Figure 3.
The motion equation becomes
m 26 L 35 0 0 2 L 3 105 w 2 ¨ ϕ 2 ¨ + 2 ρ v 0 L 5 L 5 0 w 2 ˙ ϕ 2 ˙ + E J 24 L 0 0 8 L ρ v 2 12 5 L 0 0 4 L 15 w 2 ϕ 2 = F 2 M 2
The system is now ready to be solved by determining the values of natural frequencies while the velocity changes. In particular, the load vector is set equal to zero while the displacement w and the rotation ϕ are defined as an exponential function in the complex domain ( C ) multiplied to a unknown constant (W for displacement and Φ for rotation):
w 2 = W e i ω t w 2 ˙ = t w 2 w 2 ¨ = t , t w 2 ϕ 2 = Φ e i ω t ϕ 2 ˙ = t ϕ 2 ϕ 2 ¨ = t , t ϕ 2
Substituting Equation (9) in Equation (8) and dividing everything by e i ω t results in:
2 ( 420 E J W + 42 L 2 v 2 W ρ + L 4 ω ( 7 i v ρ Φ + 13 m W ω ) ) 35 L 3 8 E J Φ L 2 105 L ( 14 v 2 ρ Φ 21 i v W ρ ω + L 2 m Φ ω 2 ) = 0
that can be transformed in the following homogeneous system collecting the coefficients of w and ϕ
A 11 A 12 A 21 A 22 w ϕ = 0
excluding the trivial solution of the Equation (10), it is possible to find the eigenvalues substituting ω 2 with λ and solving
D e t [ A f i x ] = 0
where
A f i x = A 11 A 12 A 21 A 22
leading to
λ 1 = 1 26 L 8 m 2 7 ( 840 E J L 4 m + 32 L 6 m v 2 ρ 21 L 6 v 2 ρ 2 + ( L 8 ( 518400 E J 2 m 2 720 E J L 2 m v 2 ( 40 m 49 ρ ) ρ + L 4 v 4 ρ 4 ( 400 m 2 1344 m ρ + 441 ρ 2 ) ) ) 1 / 2 ) λ 2 = 1 26 L 8 m 2 7 ( 840 E J L 4 m 32 L 6 m v 2 ρ + 21 L 6 v 2 ρ 2 + ( L 8 ( 518400 E J 2 m 2 720 E J L 2 m v 2 ( 40 m 49 ρ ) ρ + L 4 v 4 ρ 4 ( 400 m 2 1344 m ρ + 441 ρ 2 ) ) ) 1 / 2 )
The corresponding eigenvector associated to these eigenvalues describes the vibration modes of the structure. For a better representation, it is worth to transform the characteristic polynomial in its dimensionless form, introducing
θ = L 2 v 2 ρ E J
η = ρ m
so the power of the frequency is
ω 2 = E J L 4 m ( 226.2 + θ ( 8.165 + 5.654 η ) ± 0.0385 ( 2.54 e 7 + 10 6 θ ( 1.410 + 1.729 η ) + θ 2 ( 19600 65856 η + 21609 η 2 ) 1 / 2 ) )
and finally:
Ψ = L 4 m E J ω = *

3.2.2. Hinged Pipe

Let us now consider a pipe hinged to its end nodes and defined with two finite elements. These two elements have the same length L and the total amount of degrees of freedom is 4, as shown in Figure 4.
In this case, the elements of the Equation (6) are defined for four degrees of freedom derived from the imposed boundary condition and result in:
K = E L L 3 4 L 2 6 L 2 L 2 0 6 L 12 6 L 0 2 L 2 6 L 4 L 2 0 0 0 0 0 + E L L 3 0 0 0 0 0 12 6 L 6 L 0 6 L 4 L 2 2 L 2 0 6 L 2 L 2 4 L 2
M = m L 420 4 L 2 13 L 3 L 2 0 13 L 156 22 L 0 3 L 2 22 L 4 L 2 0 0 0 0 0 + m L 420 0 0 0 0 0 156 22 L 13 L 0 22 L 4 L 2 3 L 2 0 13 L 3 L 2 4 L 2
C = 2 ρ v 60 0 6 L L 2 0 6 L 30 6 L 0 L 2 6 L 0 0 0 0 0 0 + 2 ρ v 60 0 0 0 0 0 30 6 L 6 L 0 6 L 0 L 2 0 6 L L 2 0
K n l = ρ v 2 2 L 15 1 10 L 30 0 1 10 6 5 L 1 10 0 L 30 1 10 2 L 15 0 0 0 0 0 + ρ v 2 0 0 0 0 0 6 5 L 1 10 1 10 0 1 10 2 L 15 L 30 0 1 10 L 30 2 L 15
Following the procedure of the previous case, the unknown quantities are expressed as follows:
w 2 = W e i ω t w 2 ˙ = t w 2 w 2 ¨ = t , t w 2 ϕ 1 = Φ e i ω t ϕ 1 ˙ = t ϕ 1 ϕ 1 ¨ = t , t ϕ 1 ϕ 2 = Φ e i ω t ϕ 2 ˙ = t ϕ 2 ϕ 2 ¨ = t , t ϕ 2 ϕ 3 = Φ e i ω t ϕ 3 ˙ = t ϕ 3 ϕ 3 ¨ = t , t ϕ 3
Substituting Equation (15), Equation (14) in Equation (6), dividing by e i ω t and collecting the coefficient of ω and ϕ , it is possible to build the homogeneous system
A 11 A 12 A 13 A 14 A 21 A 22 A 23 A 24 A 31 A 32 A 33 A 34 A 41 A 42 A 43 A 44 w ϕ 1 ϕ 2 ϕ 3 = 0
Again, excluding the trivial solution, the natural frequencies are found solving
D e t [ A h i n g e ] = 0
where
A h i n g e = A 11 A 12 A 13 A 14 A 21 A 22 A 23 A 24 A 31 A 32 A 33 A 34 A 41 A 42 A 43 A 44

4. Results and Discussions

The results are presented as a comparison between theoretical and numerical solutions for simple applications. In particular, a fixed pipe and a hinged pipe are taken into account as case studies. In both cases, the natural frequencies, due to the geometry only and the vibration modes relative to the inner fluid velocity, are calculated. A steel pipe filled with water is considered for both applications.
Table 1 shows the data for the following problems. L is the total length, A is the section area, ρ is the fluid density, D is the diameter, d is the thickness, E is the Young module of the steel, and m 1 , m 2 are the two masses of the fluid and the pipe.

4.1. Fixed Pipe

Considering the pipe only ( η = θ = 0 because the mass of the fluid ρ = 0 ),
Equation (13) lets us calculate the two natural frequencies:
Ψ 1 = 5.684 Ψ 2 = 20.494
against the theoretical values:
Ψ 1 t = 5.593 Ψ 2 t = 15.417
In particular, the critical loads are calculated solving
D e t [ E J 24 L 3 0 0 8 L + ρ v 2 12 5 L 0 0 4 L 15 ] = 0
thus
θ = 10 θ = 30
Figure 5 shows the decrease of the frequencies increasing the velocities of the flow underlying the critical values for the pipe instability (dotted lines).
Solving Equation (18) for the data reported in Table 1, lead to these four critical velocities of the pipe:
v = 219.481 v = 126.717 v = 126.717 v = 219.481
The results for the void tube with no velocity are described solving Equation (11) imposing ρ = 0 , v = 0 and m = m 2 . The frequencies are:
ω = 223.653 ω = 62.0302 ω = 62.0302 ω = 223.653
Following the same path but considering a tube filled with water with zero velocity ( ρ = 1000 A , v = 0 and m = m 1 ) results in
ω = 151.244 ω = 41.9475 ω = 41.9475 ω = 151.244
Finally, Figure 6 shows the curves of critical frequencies ω with respect to the fluid velocity v for the first quadrant of the graph.

4.2. Hinged Pipe

Considering the pipe only and solving Equation (17), it is possible to evaluate the frequencies:
ω = 2.47714 E J L 4 m ω = 2.47714 E J L 4 m ω = 27.5349 E J L 4 m ω = 27.5349 E J L 4 m      ω = 10.9545 E J L 2 m ω = 10.9545 E J L 2 m      ω = 50.1996 E J L 2 m ω = 50.1996 E J L 2 m
and the theoretical value is:
ω t = π 2 4 L 2 E J ρ
Introducing m = m 1 for the fluid mass, it is possible to calculate the frequencies, due to the water velocity. Figure 7 shows the result for the first frequency.
Figure 8 shows all the four frequencies associated with the finite element solution of the hinged structure.
In this specific application, a comparison with the theoretical solution is possible since and the analytical solution is known:
w ( x , t ) = n A 2 n 1 S i n ( ( 2 n 1 ) π x L ) S i n ( ω i t ) + m A 2 m S i n ( 2 m π x L ) C o s ( ω i ) t
Equation (19) is substituted in Equation (1). The result can be rewritten by grouping the trigonometric terms and, in particular, it is possible to expand with a Fourier series the cosines terms:
C o s ( 1 + 2 n ) π x L = m 2 n 1 ( 2 n 1 ) 2 ( 2 m ) 2
C o s 2 m π x L = m 2 m ( 2 m ) 2 ( 2 n 1 ) 2
Setting the similar trigonometric terms equal to zero and deleting all the modes except for the first two, the following condition can be written from the homogeneous system
32 ρ v ω 3 L E J π 4 L 4 π 2 v 2 ρ L 2 m ω 2 = 3 L 8 ρ v ω 16 E J π 4 L 4 4 π 2 v 2 ρ L 2 m ω 2
Equation (20) combined with data in Table 1 defines the critical velocity
v c r = π E J ρ L 2
Figure 9 shows the comparison between the theoretical result expressed in Equation (20) and the numerical result of for the first frequency. The two results match, showing the good representation offered by a simple discretization of two elements in the finite element model.
In particular, Figure 10 shoes the error for the first frequency calculated as
E r r = ω t ω n ω t
where w t is the theoretical result of the frequency and w n is the numerical result.
The error increases with the velocity as expected, suggesting the use of more elements for high velocities.

5. Conclusions

The analysis presented in this paper shows the importance of FEM in complex phenomena. In particular, a steel tube filled with a water under Coriolis forces is studied using two element in two different configurations. The results show how the FEM solution matches the analytical one using the Hermitian shape function especially for the first frequency. In particular, the elements used in this paper are mono-dimensional, since the problem is axial-symmetrical, leading to the most efficient discretization, considering the complexity and the number of elements. However, the results are not directly comparable with other method of application because similar examples of the effect of the Coriolis force on the vibration of the annulus pipe solved by FEM are lacking in the literature, it is proved that the first frequencies of the pipe for both configurations of the filled-with-water fixed pipe and hinged pipe are stackable for theoretical results as well as for the numerical results. Moreover, the extreme simplicity of the mesh leads to a fast solution for a complex mathematical problem, avoiding heavy three-dimensional meshes and difficult boundary impositions.

6. Recommendations for Future Research

A further study should be oriented to the stress distribution, as suggested in [34] in order to investigate mechanical aspects. The use of this method could be very apt to conjure approximate closed form solutions also for nonlinear problems. Further studies and applications could be aimed also at solving the problem applied to cylindrical shell conveying fluids, whose dynamic behavior is of practical interest in the field of power plants or oil pipelines. Cylindrical shells are essential structural elements in offshore structures, submarines, and airspace crafts. They are often subjected to combined compressive stress and external pressure, and therefore must be designed to meet strength requirements [35]. As discussed in [36], these structures are often stiffened by frames or ribs. Calculation of the sound scattering properties of stiffened circular shapes is difficult. Much research has focused on the acoustic radiation from a stiffened infinite shell with a simple shape but, when the shape is more complicated, such as for a submarine or an aircraft, the numerical methods are not appropriate, while the finite element method seems justified with the finite element characterizing the stiffeners. This phenomenon could be applied to modern tools such as Virtual or Augmented Reality. In particular, Augmented Reality (AR) is a computer technology where the perception of the user is enhanced by the seamless blending between a realistic environment and computer-generated virtual objects coexisting in the same space. The resulting mixture supplements reality, rather than replacing it. The possibility of interacting with external information could be very useful in simulations difficult to be interpreted by common users. As suggested in [37,38,39,40,41,42,43,44,45,46,47], Augmented Reality applications could have many interesting advantages and could be a viable tool in many fields now investigated by Industry 4.0.

Author Contributions

Conceptualization F.C.; methodology and writing G.M.S.; writing and validation D.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

FEMFinite Element Method
LNGLiquefied Natural Gas
LESLarge Eddy Simulations
CTNCarbon NanoTubes
SEMSpectral Element Method
SWCNTSingle-Walled Carbon NanoTubes
FSIFluid Structure Interaction
ARAugmented Reality

References

  1. Ibrahim, R.A. Overview of Mechanics of Pipes Conveying Fluids—Part I: Fundamental Studies. J. Press. Vessel. Technol. 2010, 132. [Google Scholar] [CrossRef]
  2. Paidoussis, M. Fluid-Structure Interactions: Slender Structures and Axial Flow; Academic Press: Cambridge, MA, USA, 2004; Volume 2. [Google Scholar]
  3. Lee, U. Spectral element modelling and analysis of a pipeline conveying internal unsteady fluid. J. Fluids Struct. 2006, 22, 273–292. [Google Scholar] [CrossRef]
  4. Dodds, H.L.; Runyan, H.L. Effect of High Velocity Fluid Flow on the Bending Vibrations and Static Divergence of a Simply Supported Pipe; National Aeronautics and Space Administration: Washington, DC, USA, 1965.
  5. Yamaguchi, R.; Tanaka, G.; Liu, H.; Hayase, T. Fluid Vibration Induced in T-Junction with Double Side Branches. World J. Mech. 2016, 6, 169–179. [Google Scholar] [CrossRef] [Green Version]
  6. Kuiper, G.; Metrikine, A. On stability of a clamped-pinned pipe conveying fluid. Heron 2004, 49, 211–232. [Google Scholar]
  7. Chellapilla, K.R.; Simha, H. Vibrations of Fluid-Conveying Pipes Resting on Two-parameter Foundation. Open Acoust. J. 2008, 1, 24–33. [Google Scholar] [CrossRef] [Green Version]
  8. Leklong, J.; Chucheepsakul, S.; Kaewunruen, S. Dynamic Responses of Marine Risers/Pipes Transporting Fluid Subject to Top End Excitations; Faculty of Engineering, UOW: Wollongong, Australia, 2008. [Google Scholar]
  9. Murai, M.; Yamamoto, M. An Experimental Analysis of the Internal Flow Effects on Marine Risers. Proc. MARTEC 2010, 1, 159–165. [Google Scholar]
  10. Marakala, N.; Kuttan, A.K.K.; Kadoli, R. Experimental and Theoretical Investigation of Combined Effect of Fluid and Thermal Induced Vibration on Vertical Thin Slender Tube. Iosr J. Mech. Civ. Eng. 2013, 63–68. [Google Scholar]
  11. Eriksson, O.; Breton, S.P.; Nilsson, K.; Ivanell, S. Impact of Wind Veer and the Coriolis Force for an Idealized Farm to Farm Interaction Case. Appl. Sci. 2019, 9, 922. [Google Scholar] [CrossRef] [Green Version]
  12. Shaik, I.; Uddien, S.; Arkanti, K.; Sutar, S. Numerical Analysis of Clamped Fluid Conveying Pipe. HAL 2017. [Google Scholar] [CrossRef]
  13. Kokare, A.B.; Paward, P.M. Vibration Characteristics of Pipe Conveying Fluid. Int. J. Sci. Eng. Technol. Res. 2015, 4, 10951–10954. [Google Scholar]
  14. Li, G.; Liu, J.; Jiang, G.; Kong, J.; Xie, L.; Xiao, W.; Zhang, Y.; Cheng, F. The nonlinear vibration analysis of the fluid conveying pipe based on finite element method. Comput. Model. New Technol. 2014, 18, 19–24. [Google Scholar]
  15. Li, G.; Xiao, W.; Jiang, G.; Liu, J. Finite element analysis of fluid conveying pipe line of nonlinear vibration response. Comput. Model. New Technol. 2014, 18, 37–41. [Google Scholar]
  16. Jweeg, M.J.; Ntayeesh, T.J. Dynamic analysis of pipes conveying fluid using analytical, numerical and experimental verification with the aid of smart materials. Int. J. Sci. Res. 2015, 4, 1594–1605. [Google Scholar]
  17. Yi, X.; Li, B.; Wang, Z. Vibration Analysis of Fluid Conveying Carbon Nanotubes Based on Nonlocal Timoshenko Beam Theory by Spectral Element Method. Nanomaterials 2019, 9, 1780. [Google Scholar] [CrossRef] [Green Version]
  18. Ferras, D.; Manso, P.; Schleiss, A.; Covas, D. One-Dimensional Fluid-Structure Interaction Models in Pressurized Fluid-Filled Pipes: A Review. Appl. Sci. 2018, 8, 1844. [Google Scholar] [CrossRef] [Green Version]
  19. Elishakoff, I. Controversy Associated with the So-Called “Follower Forces”: Critical Overview. Appl. Mech. Rev. 2005, 58, 117. [Google Scholar] [CrossRef]
  20. Öz, H.R.; Boyaci, H. Transverse vibrations of tensioned pipes conveying fluid with time-dependent velocity. J. Sound Vib. 2000, 236, 259–276. [Google Scholar] [CrossRef]
  21. Szmidt, T.; Przybyłowicz, P. Critical flow velocity in a pipe with electromagnetic actuators. J. Theor. Appl. Mech. 2013, 51, 487–496. [Google Scholar]
  22. Askarian, A.; Abtahi, H.; Haddadpour, H. Dynamic instability of cantilevered composite pipe conveying flow with an end nozzle. In Proceedings of the 21st International Congress on Sound and Vibration 2014, ICSV 2014, Beijing, China, 13–17 July 2014; Volume 4, pp. 3564–3571. [Google Scholar]
  23. Ntayeesh, T.; Alhelli, A. Free Vibration Characteristics of Elastically Supported Pipe Conveying Fluid. J. Eng. Sci. 2013, 16, 9–19. [Google Scholar]
  24. Guo, Q.; Zhang, L.; Xiao, L. Damage analysis of the vehicle’s pipe conveying fluid induced by complex random-shock loads. J. Press. Equip. Syst. 2006, 4, 100–103. [Google Scholar]
  25. Zhang, Y.L.; Gorman, D.G.; Reese, J.M.; Horacek, J. Observations on the vibration of axially tensioned elastomeric pipes conveying fluid. J. Mech. Eng. Sci. 2000, 214, 423–434. [Google Scholar] [CrossRef] [Green Version]
  26. Modarres-Sadeghi, Y.; Païdoussis, M. Nonlinear dynamics of extensible fluid conveying pipes, supported at both ends. J. Fluids Struct. 2009, 25, 535–543. [Google Scholar] [CrossRef]
  27. Sharma, J.; Najafi, M.; Marshall, D.; Kaushal, V.; Hatami, M. Development of a Model for Estimation of Buried Large-Diameter Thin-Walled Steel Pipe Deflection due to External Loads. J. Pipeline Syst. Eng. Pract. 2019, 10. [Google Scholar] [CrossRef] [Green Version]
  28. Malek Mohammadi, M.; Najafi, M.; Kermanshachi, S.; Kaushal, V.; Serajiantehrani, R. Factors Influencing the Condition of Sewer Pipes: A State-of-the-Art Review. J. Pipeline Syst. Eng. Pract. 2020, 11. [Google Scholar] [CrossRef]
  29. Olunloyo, V.O.S.; Osheku, C.A.; Olayiwola, P.S. A Note on an Analytic Solution for an Incompressible Fluid Conveying Pipeline System. Adv. Acoust. Vib. 2017, 2017, 20. [Google Scholar] [CrossRef]
  30. Paidoussis, M.; Issid, N. Dynamic stability of pipes conveying fluid. J. Sound Vib. 1974, 33, 267–294. [Google Scholar] [CrossRef]
  31. Semler, C.; Li, G.; Paidoussis, M. The Non-linear Equations of Motion of Pipes Conveying Fluid. J. Sound Vib. 1994, 169, 577–599. [Google Scholar] [CrossRef]
  32. Enz, S.; Thomsen, J.J. Predicting phase shift effects for vibrating fluid conveying pipes due to Coriolis forces and fluid pulsation. J. Sound Vib. 2011, 330, 5096–5113. [Google Scholar] [CrossRef]
  33. Bathe, K. Finite Element Procedures in Engineering Analysis; Prentice Hall Inc.: Upper Saddle River, NJ, USA, 2014. [Google Scholar]
  34. Zachwieja, J. Stress analysis of vibrating pipelines. Aip Conf. Proc. 2017, 1822, 020017. [Google Scholar] [CrossRef] [Green Version]
  35. Bai, Y.; Jin, W.L. Ultimate Strength of Cylindrical Shells; Elsevier: Oxford, UK, 2016; pp. 353–365. [Google Scholar] [CrossRef]
  36. Debus, J.C. 4-Finite element analysis of flexural vibration of orthogonally stiffened cylindrical shells with ATILA. In Applications of ATILA FEM Software to Smart Materials; Uchino, K., Debus, J.C., Eds.; Woodhead Publishing Series in Electronic and Optical Materials; Woodhead Publishing: New Delhi, India, 2013; pp. 69–93. [Google Scholar] [CrossRef]
  37. Ceruti, A.; Marzocca, P.; Liverani, A.; Bil, C. Maintenance in aeronautics in an Industry 4.0 context: The role of Augmented Reality and Additive Manufacturing. J. Comput. Des. Eng. 2019, 6, 516–526. [Google Scholar] [CrossRef]
  38. Osti, F.; Ceruti, A.; Liverani, A.; Caligiana, G. Semi-automatic Design for Disassembly Strategy Planning: An Augmented Reality Approach. In Proceedings of the 27th International Conference on Flexible Automation and Intelligent Manufacturing, FAIM2017, Modena, Italy, 27–30 June 2017; Volume 11, pp. 1481–1488. [Google Scholar] [CrossRef]
  39. Frizziero, L.; Liverani, A.; Caligiana, G.; Donnici, G.; Chinaglia, L. Design for disassembly (DfD) and augmented reality (AR): Case study applied to a gearbox. Machines 2019, 7, 29. [Google Scholar] [CrossRef] [Green Version]
  40. De Marchi, L.; Ceruti, A.; Marzani, A.; Liverani, A. Augmented Reality to Support On-Field Post-Impact Maintenance Operations on Thin Structures. J. Sens. 2013, 2013, 619570. [Google Scholar] [CrossRef] [Green Version]
  41. Liverani, A.; Leali, F.; Pellicciari, M. Real-time 3D features reconstruction through monocular vision. Int. J. Interact. Des. Manuf. 2010, 4, 103–112. [Google Scholar] [CrossRef]
  42. Liverani, A.; Kuester, F.; Hamann, B. Towards interactive finite element analysis of shell structures in virtual reality. In Proceedings of the 1999 IEEE International Conference on Information Visualization (Cat. No. PR00210), London, UK, 14–16 July 1999; pp. 340–346. [Google Scholar] [CrossRef] [Green Version]
  43. Francia, D.; Seminerio, D.; Caligiana, G.; Frizziero, L.; Liverani, A.; Donnici, G. Virtual Design for Assembly Improving the Product Design of a Two-Way Relief Valve. In Proceedings of the International Conference on Design, Simulation, Manufacturing: The Innovation Exchange, Kharkiv, Ukraine, 9–12 June 2020; pp. 304–314. [Google Scholar] [CrossRef]
  44. Osti, F.; Santi, G.C.G. Real Time Shadow Mapping for Augmented Reality Photorealistic Rendering. Appl. Sci. 2019, 9, 2225. [Google Scholar] [CrossRef] [Green Version]
  45. Francia, D.; Ponti, S.; Frizziero, L.; Liverani, A. Virtual Mechanical Product Disassembly Sequences Based on Disassembly Order Graphs and Time Measurement Units. Appl. Sci. 2019, 9, 3638. [Google Scholar] [CrossRef] [Green Version]
  46. Amicis, R.; Ceruti, A.; Francia, D.; Frizziero, L.; Simões, B. Augmented Reality for virtual user manual. Int. J. Interact. Des. Manuf. 2018, 12. [Google Scholar] [CrossRef]
  47. Frizziero, L.; Donnici, G.; Caligiana, G.; Liverani, A.; Francia, D. Project of Inventive Ideas Through a TRIZ Study Applied to the Analysis of an Innovative Urban Transport Means. Int. J. Manuf. Mater. Mech. Eng. 2018, 8, 35–62. [Google Scholar] [CrossRef]
Figure 1. Methodology flow chart.
Figure 1. Methodology flow chart.
Applsci 11 01058 g001
Figure 2. Fluid moving through a pipe.
Figure 2. Fluid moving through a pipe.
Applsci 11 01058 g002
Figure 3. Element representation of the pipe fixed to its boundary nodes.
Figure 3. Element representation of the pipe fixed to its boundary nodes.
Applsci 11 01058 g003
Figure 4. Element representation of the pipe hinged to its boundary nodes.
Figure 4. Element representation of the pipe hinged to its boundary nodes.
Applsci 11 01058 g004
Figure 5. Frequencies with respect to the velocity.
Figure 5. Frequencies with respect to the velocity.
Applsci 11 01058 g005
Figure 6. Frequencies respect to the velocity.
Figure 6. Frequencies respect to the velocity.
Applsci 11 01058 g006
Figure 7. Frequencies with respect to the velocity.
Figure 7. Frequencies with respect to the velocity.
Applsci 11 01058 g007
Figure 8. Frequencies with respect to the velocity.
Figure 8. Frequencies with respect to the velocity.
Applsci 11 01058 g008
Figure 9. Frequencies with respect to the velocity.
Figure 9. Frequencies with respect to the velocity.
Applsci 11 01058 g009
Figure 10. Error between theoretical and numerical result for the first frequency.
Figure 10. Error between theoretical and numerical result for the first frequency.
Applsci 11 01058 g010
Table 1. Input data for a steel pipe filled with water.
Table 1. Input data for a steel pipe filled with water.
L = 4 m D = 0.1 m d = 0.095 m
A = π d 2 4 m 2 E = 2 · 10 11 Pa J = π ( D 2 d 2 ) 64
ρ = 1000 π d 2 4 m 1 = π ( D 2 d 2 ) 4 7800 + 1000 A m 2 = π ( D 2 d 2 ) 4 7800
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Santi, G.M.; Francia, D.; Cesari, F. Effect of Coriolis Force on Vibration of Annulus Pipe. Appl. Sci. 2021, 11, 1058. https://doi.org/10.3390/app11031058

AMA Style

Santi GM, Francia D, Cesari F. Effect of Coriolis Force on Vibration of Annulus Pipe. Applied Sciences. 2021; 11(3):1058. https://doi.org/10.3390/app11031058

Chicago/Turabian Style

Santi, Gian Maria, Daniela Francia, and Francesco Cesari. 2021. "Effect of Coriolis Force on Vibration of Annulus Pipe" Applied Sciences 11, no. 3: 1058. https://doi.org/10.3390/app11031058

APA Style

Santi, G. M., Francia, D., & Cesari, F. (2021). Effect of Coriolis Force on Vibration of Annulus Pipe. Applied Sciences, 11(3), 1058. https://doi.org/10.3390/app11031058

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