Analytical Solution of a Forced Vibrating Tube Transporting Fluid

Moussa Tembely

doi:10.20944/preprints202401.1962.v2

Preprint

Article

Analytical Solution of a Forced Vibrating Tube Transporting Fluid

This version is not peer-reviewed.

Moussa Tembely^*

This version is not peer-reviewed.

Downloads

180

Views

Comments

Submitted:

09 March 2024

Posted:

11 March 2024

Read the latest preprint version here

Abstract

In this paper is proposed an analytical technique which can be used to obtain the forced response of a cantilevered tube conveying fluid. Considering the pipe subjected to an arbitrary harmonic force either distributed or concentrated, an analytical solution is found using Green’s function method. The solution obtained of the closed form satisfies the differential equations in a classical sense. The presented method giving the exact solutions is more precise than the classical eigenfunction expansion or Galerkin’s method with no need to the eigenfunctions or eigenvalues.

Keywords:

Tube conveying fluid

;

forced vibration

;

Analytical solution

Subject:

Engineering - Mechanical Engineering

1. Introduction

Pipes or tubes carrying fluid play a critical role in a wide range of industrial and engineering applications, including fuel pipelines, inkjet print-head nozzles, fuel-carrying tubes in space shuttle engines, and heat exchangers in nuclear power plants; thus, forced tube conveying is clearly desired [1,2]. It is also worth noting that feedback control for beams might be used for forced tube vibration [3].

From a theoretical standpoint, the fluid/structure interaction associated with the boundary moving condition is undoubtedly one of the most interesting and hardest problems in dynamical systems. For nearly 50 years, these issues have been studied both theoretically and empirically. The challenges investigated varied from changes in the natural frequencies of the pipe system caused by flowing fluid, to vibrations caused by turbulent flows, and instabilities such as flutter instability (via a Hopf bifurcation) which occurs when the flow velocity reaches a critical velocity [4]. Following [5], there are two aspects of forced vibration of pipes conveying fluid: the first is the physical aspect, which sheds more light on the dynamics of the system, and the second is related to the analytical techniques that can be used to obtain the forced response of such systems, which is the one used in this work, though the result could also be extended to the first aspect.

This paper develops an improved analytical strategy based on Green's function method to generate closed form solutions of tube carrying fluid that have continuous derivatives and satisfy the differential equations in a classical sense. Furthermore, the method does not involve solving the free vibration problem in order to obtain the eigenvalues and related eigenfunctions, which are required when utilizing classical Galerskin's method.

2. Tube Conducting Fluid General Equation of Motion

The system under consideration is a straight cantilever tube fixed to the x-axis. The fluid with an average velocity of U enters the tube at the fixed end and departs at the free end. A force is provided to the tube, allowing it to move transversely. Furthermore, the force that induces tube motion is assumed to be harmonic. Unlike the unconstrained pipe dynamical stability achieved by many others [5,6], the goal here is considerably different in the sense that we want to know the precise deflection of the tube, with the harmonic force considered either distributed or concentrated.

Considering the effects of dissipation and damping, the general equation of motion for a tube conveying fluid can be formulated based on the following assumptions: (i) the tube has a uniform annular cross-section; (ii) the length of the tube is significantly greater than its diameter; (iii) the effects of rotary inertia and shear deformation are neglected; (iv) the tube's center line is inextensible; (v) the tube is elastic and initially straight; and (vi) cross-sections of the tube remain plane and perpendicular to the axis of the beam, in accordance with Bernoulli-Euler beam theory. The equation can be expressed as follows [1,7]:

\frac{\partial^{4} η}{\partial x^{4}} + u^{2} \frac{\partial^{2} η}{\partial^{2} x} + 2 β^{1 / 2} u \frac{\partial^{2} η}{\partial t \partial x} + σ \frac{\partial η}{\partial t} + \frac{\partial^{2} η}{\partial t^{2}} = F (x, t)

(1)

The resultant equation is made dimensionless by employing for variables,

x \to \frac{x}{L}, η \to \frac{η}{L}, t \to {(\frac{E I}{m_{t} + M_{f}})}^{1 / 2} \frac{t}{L^{2}};

(2)

and for constants,

u = L U (\frac{M_{f}}{E I}), β = \frac{M_{f}}{m_{t} + M_{f}}, σ = \frac{c L^{2}}{{[E I (m_{t} + M_{f})]}^{1 / 2}}

(3)

In the following section, we will solve the dimensionless form of the equation for a cantilever conveying fluid, which is subject to a harmonic forcing function.

3. Exact Analytical Solution

We propose a new approach based on Green’s function method to solve equation (1). Assuming the applied force is harmonic, we could formulate the problem in a complex form. Let’s consider a tube subjected to an arbitrary harmonic force, such that it is governed by an equation of the form,

F (x, t) = f (x) e^{i ω t}

(4)

where

i

is the imaginary unit.

For general purpose, the force could be either distributed or concentrated at a position

x_{p}

thus,

f (x) = {\begin{matrix} F_{p} δ (x - x_{p}) \\ f (x) \end{matrix},

(5)

therefore, the solution of the problem is sought as:

w (x, t) = V (x) e^{i ω t}

(6)

where

V (x)

is the unknown function to be determined.

As consequence of (4) and (6) in complex form, the solution of the problem is the real part of

w (x, t)

as we will detail later. By replacing (4) in (1) we obtain,

V'''' (x) + a V'' (x) + i b V' (x) - c V (x) = f (x)

(7)

where we set

a = u^{2}

b = 2 ω β^{1 / 2}

and

c = ω^{2} - 1 - i σ ω

where prime (’) denotes for derivative with respect to x.

The solution of (7) could be expressed using green function as,

V (x) = \int_{0}^{L} f (x^{'}) G (x, x^{'}) d x^{'}

(8)

where we can determine the Green's function

by solving the following equation:

{V^{‴}}^{'} (x) + a V^{″} (x) + i b V^{'} (x) - c V (x) = δ (x - x^{'})

(9)

We solve this equation by adopting the Laplace transform and its inverse, deducing therefore,

\begin{array}{l} V (x, x^{'}) = H (x - x^{'}) \sum_{α_{i}} \frac{e^{α_{i} (x - x^{'})}}{4 α_{i}^{3} + 2 α_{i} a + i b} \\ + V (0) \sum_{α_{i}} \frac{α_{i}^{3} e^{α_{i} x}}{4 α_{i}^{3} + 2 α_{i} a + i b} + V^{'} (0) \sum_{α_{i}} \frac{α_{i}^{2} e^{α_{i} x}}{4 α_{i}^{3} + 2 α_{i} a + i b} \\ + \sum_{α_{i}} \frac{e^{α_{i} x}}{4 α_{i}^{3} + 2 α_{i} a + i b} (V^{‴} (0) + a V^{'} (0) + b i V (0)) \\ + \sum_{α_{i}} \frac{α_{i} e^{α_{i} x}}{4 α_{i}^{3} + 2 α_{i} a + i b} (V^{″} (0) + a V (0)) \end{array}

(10)

where

H (x)

is the Heaviside or step function and

α_{i}, i = 1...4

is the solution of the following polynomial equation:

Z^{4} + a Z^{2} + i b Z - c = 0

(11)

This quartic equation can be analytically solved using Ferrari's method.

Let the function

ψ (x)

be defined as:

ψ (x) = \sum_{α_{i}} \frac{e^{α_{i} x}}{4 α_{i}^{3} + 2 α_{i} a + i b}

(12)

We then can recast (10) in the following form,

V (x, x^{'}) = H (x - x^{'}) ψ (x - x_{p}) + ψ^{'} (x) V^{″} (0) + ψ (x) V^{‴} (0)

(13)

where we use the boundary conditions at

x = 0

with,

V (0) = V' (0) = 0

The last step is now to find the other boundary conditions namely

V'' (0)

and

V''' (0)

present in (13). For this purpose, we will use the boundary conditions at the free end of the cantilevered tube; from the bending moment at

x = L

V^{″} (L) = 0 = ψ^{″} (L - x^{'}) + ψ^{‴} (L) V^{″} (0) + ψ^{″} (L) V^{‴} (0)

(14)

and from the shear force at

x = L

V^{‴} (L) = 0 = ψ^{‴} (L - x^{'}) + ψ^{(4)} (L) V^{″} (0) + ψ^{‴} (L) V^{‴} (0)

(15)

Solving the previous equations (14-15) leads to determine the unknown boundary conditions at

x = 0

as:

V^{″} (0) = \frac{ψ^{″} (L) ψ^{‴} (L - x^{'}) - ψ^{‴} (L) ψ^{″} (L - x^{'})}{ψ^{‴} {(L)}^{2} - ψ^{″} (L) ψ^{(4)} (L)}

(16)

V^{‴} (0) = \frac{ψ^{(4)} (L) ψ^{″} (L - x^{'}) - ψ^{‴} (L) ψ^{‴} (L - x^{'})}{ψ^{‴} {(L)}^{2} - ψ^{″} (L) ψ^{(4)} (L)}

(17)

Thus, we deduce the Green's function of the cantilevered tube conveying fluid as:

\begin{array}{l} G (x, x^{'}) = H (x - x^{'}) ψ (x - x^{'}) + \frac{ψ^{″} (L) ψ^{‴} (L - x^{'}) - ψ^{‴} (L) ψ^{″} (L - x^{'})}{ψ^{‴} {(L)}^{2} - ψ^{″} (L) ψ^{(4)} (L)} ψ^{″} (x) \\ + \frac{ψ^{(4)} (L) ψ^{″} (L - x^{'}) - ψ^{‴} (L) ψ^{‴} (L - x^{'})}{ψ^{‴} {(L)}^{2} - ψ^{″} (L) ψ^{(4)} (L)} ψ (x) \end{array}

(18)

Finally, the dynamics response of a forced cantilevered tube conducting fluid subjected to an arbitrary harmonic force expresses as:

w (x) = e^{i ω t} \int_{0}^{L} f (x^{'}) G (x, x^{'}) d x^{'}

(19)

f (x)

could be taken either concentrated or distributed (5).

In the case where

f (x) = F_{p} δ (x - x_{p})

, we could establish the analytical solution of the tube conveying fluid submitted to a pointwise harmonic excitation as:

η (x, t) = ℜ [w (x, t)]

where

w (x, t)

is given by,

\begin{array}{l} w (x, t) = F_{p} \exp (i ω) {H (x - x_{p}) ψ (x - x_{p}) \\ + \frac{ψ^{″} (L) ψ^{‴} (L - x_{p}) - ψ^{‴} (L) ψ^{″} (L - x_{p})}{ψ^{‴} (L) ψ^{(4)} (L) - ψ^{″} (L) ψ^{(5)} (L)} ψ^{″} (x) \\ + \frac{ψ^{(5)} (L) ψ^{″} (L - x_{p}) - ψ^{(4)} (L) ψ^{‴} (L - x_{p})}{ψ^{‴} (L) ψ^{(4)} (L) - ψ^{″} (L) ψ^{(5)} (L)} ψ (x)} \end{array}

(20)

Thus, the solution finally is given by:

v (x, t) = ℜ [w (x, t)] \cos (ω t) - ℑ [w (x, t)] \sin (ω t)

(21)

4. Conclusion

Considering a general equation of motion for a tube carrying fluid subjected to harmonic excitation, an analytical technique based on Green's function method has been developed to find a solution for a tube carrying fluid. The derived exact solution is more accurate than the conventional series expansion or Galerkin's method, and it does not require knowledge of the free vibration problem's eigenfunctions or eigenvalues. The proposed method might also be extended to include moment solicitation and other configurations such as multi-span tubes, boundary conditions, and even tube dynamical analysis.

References

M. P. Paidoussis, Fluid-Structure Interactions: Slender Structures and Axial Flow, 2nd edition. Amsterdam; Boston MA: Academic Press, 2014. [Google Scholar]
Tembely, M.; Lecot, C.; Soucemarianadin, A. Prediction and evolution of drop-size distribution for a new ultrasonic atomizer. Appl. Therm. Eng. 2011, 31, 656–667. [Google Scholar] [CrossRef]
Lin, Y.-H.; Chu, C.-L. ACTIVE FLUTTER CONTROL OF A CANTILEVER TUBE CONVEYING FLUID USING PIEZOELECTRIC ACTUATORS. J. Sound Vib. 1996, 196, 97–105. [Google Scholar] [CrossRef]
Housner, G.W. Bending Vibrations of a Pipe Line Containing Flowing Fluid. J. Appl. Mech. 2021, 19, 205–208. [Google Scholar] [CrossRef]
Weaver, D.S.; Unny, T.E. On the Dynamic Stability of Fluid-Conveying Pipes. J. Appl. Mech. 1973, 40, 48–52. [Google Scholar] [CrossRef]
Païdoussis, M.; Semler, C.; Wadham-Gagnon, M.; Saaid, S. Dynamics of cantilevered pipes conveying fluid. Part 2: Dynamics of the system with intermediate spring support. J. Fluids Struct. 2007, 23, 569–587. [Google Scholar] [CrossRef]
Ammari, K.; Tucsnak, M. Stabilization of Bernoulli--Euler Beams by Means of a Pointwise Feedback Force. SIAM J. Control. Optim. 2000, 39, 1160–1181. [Google Scholar] [CrossRef]

Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

© 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

Downloads

180

Views

Comments

Subscription

Notify me about updates to this article or when a peer-reviewed version is published.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

Analytical Solution of a Forced Vibrating Tube Transporting Fluid

Abstract

Keywords:

Subject:

1. Introduction

2. Tube Conducting Fluid General Equation of Motion

3. Exact Analytical Solution

4. Conclusion

References

MDPI Initiatives

Important Links

Subscribe