The Analysis of Nonlinear Vibrations of Top-Tensioned Cantilever Pipes Conveying Pressurized Steady Two-Phase Flow under Thermal Loading

This paper studied the nonlinear vibrations of top-tensioned cantilevered pipes conveying pressurized steady two-phase flow under thermal loading. The coupled axial and transverse governing partial differential equations of motion of the system were derived based on Hamilton’s mechanics, with the centerline assumed to be extensible. Using the multiple-scale perturbation technique, natural frequencies, mode shapes, and first order approximate solutions of the steady-state response of the pipes were obtained. The multiple-scale assessment reveals that at some frequencies the system is uncoupled, while at some frequencies a 1:2 coupling exists between the axial and the transverse frequencies of the pipe. Nonlinear frequencies versus the amplitude displacement of the cantilever pipe, conveying two-phase flow at super-critical mixture velocity for the uncoupled scenario, exhibit a nonlinear hardening behavior; an increment in the void fractions of the two-phase flow results in a reduction in the pipe’s transverse vibration frequencies and the coupled amplitude of the system. However, increases in the temperature difference, pressure, and the presence of top tension were observed to increase the pipe’s transverse vibration frequencies without a significant change in the coupled amplitude of the system.


Introduction
Two-phase flow is a common flow phenomenon in various industrial pipes: in nuclear heat exchangers, pipes in process plants, thermal plants, subsea oil and gas explorations, and many more. However, in spite of the vast occurrences of two-phase flow in pipes, most of the existing publications on the flow-induced vibrations of pipes conveying fluids focus on the fluidelastic instability of pipes conveying single-phase flow. Miwa et al. [1] did an in-depth review of the extent of existing work on two-phase flow-induced vibrations, stating that there exist very few results on the instability behavior of pipes due to internal two-phase flow. In the review, it was explained that internal two-phase flow-induced vibration can be initiated by various hydrodynamic phenomena; depending on the geometrical configurations of the flow channels and operating conditions, gas-liquid two-phase flow may create vibrations with various modes of amplitude and frequency. Monette and Pettigrew [2] presented an excellent experimental work on the fluidelastic instability of flexible tubes subjected to two-phase flow which might be one of the premier papers on the dynamics of pipes conveying two-phase flow and also reveals the relationship between the void fraction and the linear dynamics of the pipe for a two-phase liquid-gas flow. The early studies on the nonlinear dynamics of cantilevered pipes conveying single-phase flow either studied only the transverse displacement of the pipe or considered the longitudinal displacement using inextensible centerline assumptions [3][4][5]. However, the pioneering work by Ghayesh et al. [6] studied the nonlinear dynamics of a cantilevered extensible pipe conveying fluid, with equations of motion of the coupled transverse and longitudinal displacements derived using the Lagrange equations for a system containing non-material volumes, and highlighted that, conversely to inextensible pipe, an extensible pipe elongates in the axial direction as the flow velocity increases. Luczko et al. [7] highlighted that the dynamic behavior of continuous systems, such as beams, moving tapes, or pipes with flowing fluid is governed by non-linear partial differential equations with appropriate boundary and initial conditions. As highlighted by Païdoussis [8], the nonlinear problems of pipes conveying fluids cannot be resolved analytically, but recourse needs to be taken to adopt specialized analytical methods like perturbation techniques, numerical time difference methods, or a combined analytical-numerical method. The direct Lagrangian discretization method (Galerkin method) to convert the partial differential equations (PDEs) to ordinary differential equations (ODEs) and then resolve the resulting ODEs using numerical techniques has been adopted by some publications, namely, Modarres and Païdoussis [9], Wang et al. [10], Sinir [11], Ritto et al. [12], Chen et al. [13], Tian-Zhi et al. [14]. The usage of analytical methods like perturbation techniques is highly common with researchers working on nonlinear problems, such as Nayfeh [15], Nayfeh [16], Kesimli et al. [17], and Oz and Boyaci [18], where the solutions were sought using an asymptotic expansion or by perturbing the original set of equations in terms of a small parameter which is either present in the equation or introduced artificially. Some researchers working on the fluidelastic instability of pipes conveying fluids have adopted this technique to resolve the nonlinear dynamics of the pipes, including the works of Enz [19] on a simple supported straight pipe using perturbation analysis with a multiple time-scaled method and comparison with measurements made by Coriolis flowmeters; the study by Xiao-Dong et al. [20] on the dynamic stability simply-supported viscoelastic pipe in transverse vibration for conveying pulsating fluid; the study on the transverse vibrations of tension pipes conveying fluid with time-dependent velocity using the multiple-scale perturbation technique by Oz and Pakdemirli [21]; and the study on the analysis of nonlinear vibrations of a pipe conveying an ideal fluid by Sinir and Demir [22].
Most of these existing publications on the nonlinear dynamics of a cantilevered pipe conveying fluid were focused on single-phase flow, resolving the governing equations using various methods as highlighted in the review of literature. However, profound among authors who adopted the perturbation approach is the study of the uncoupled problem, solving only the transverse vibrations independent of the axial vibration. To the best of our knowledge, a perturbation approach for the resolution of the coupled nonlinear dynamics of a top-tensioned cantilevered pipe conveying pressurized two-phase flow under thermal loading is still a reserved topic with a possibly intriguing outcome. In this paper, the governing equation of motion for the nonlinear coupled axial and transverse vibrations of a cantilever pipe conveying two-phase flow is derived and resolved by imposing the multiple-scale perturbation technique directly to the system equations (direct-perturbation method).

Problem Formulation and Modelling
Consider a system comprising a cantilever cylindrical pipe of length (L), having a cross-sectional area (A), mass per unit length (m p ) and flexural rigidity (EI), conveying multiphase flow, flowing parallel to the pipe's centerline. The centerline axis of the pipe in its undeformed state is assumed to overlap with the Y-axis, and the cylinder is assumed to vibrate in the (Y, X) plane (see Figure 1). To derive the system's governing equations of motion, the following basic assumptions were made for the cylinder and the fluid: (i) the mean flow velocity is constant; (ii) the cylinder is slender, so that the Euler-Bernoulli beam theory is applicable; (iii) although the deflections of the cylinder may be large, the strains are small; and (iv) the cylinder centerline is extensible. The centerline of the cantilever pipe is assumed to be extensible to account for possible expansion due to the high temperature of the fluid content. The strain expressions and the geometric relation of the centerline of an extensible pipe are derived as expressed by Semler et al. [3].

Derivation of the Equation of Motion
The equations of motion are derived using the energy method. The energy method is based on Hamilton's principle, which is defined as the variations of the time derivative of the Lagrangian. This can be mathematically expressed as where is the number of phases in the fluid, which will be 2 for a two-phase flow; is the mass of the phases in the fluid; is the flow velocity of the phases in the fluid; and ℒ is the Lagrangian operator expressed in Equation (2): where and are the kinetic and potential energies of the pipe, respectively, and and are the kinetic and potential energies associated with the conveyed fluid, respectively.

Kinetic Energy
The total kinetic energy of the system is the summation of the kinetic energy of the pipe and the kinetic energies of the phases/components of the flowing fluid. The velocity vector of the pipe's centerline is expressed as Therefore, the kinetic energy of the pipe is expressed as The centerline of the cantilever pipe is assumed to be extensible to account for possible expansion due to the high temperature of the fluid content. The strain expressions and the geometric relation of the centerline of an extensible pipe are derived as expressed by Semler et al. [3].

Derivation of the Equation of Motion
The equations of motion are derived using the energy method. The energy method is based on Hamilton's principle, which is defined as the variations of the time derivative of the Lagrangian. This can be mathematically expressed as where n is the number of phases in the fluid, which will be 2 for a two-phase flow; M j is the mass of the phases in the fluid; U j is the flow velocity of the phases in the fluid; and L is the Lagrangian operator expressed in Equation (2): where T p and V p are the kinetic and potential energies of the pipe, respectively, and T f and V f are the kinetic and potential energies associated with the conveyed fluid, respectively. The following expressions hold: .

Kinetic Energy
The total kinetic energy of the system is the summation of the kinetic energy of the pipe and the kinetic energies of the phases/components of the flowing fluid. The velocity vector of the pipe's centerline is expressed as Therefore, the kinetic energy of the pipe is expressed as As illustrated by Semler, the axial elongation of the pipe is complemented by a lateral contraction, due to the Poisson effect. This will impact the flow velocity of the fluid phases/components such that where U 0 j and U j are the flow velocities before and after elongation, the subscript (j) is used to identify the various phases/components of the conveyed fluid, ( ) is the axial strain, and (a) relates to the Poisson ratio (v) as a = 1 − 2v; for an extreme case, v = 0.5 and a becomes zero, as in Ghayesh et al. [6]. The flow velocity relative to the centerline axis of the pipe is expressed as Therefore, the Kinetic energy of the conveyed fluid is expressed as

Potential Energy
Semler highlighted that the potential energy is a result of the elastic deformation of the pipe and the effect of gravity, and the deformation from elastic behavior of the pipe can be linked to the strain energy. This is expressed as This is clearly the combination of the axial strain effect and the bending strain effect where (E) denotes the Young's modulus, (I) denotes the pipe moment of inertia, (A) denotes the cross-sectional area, ( ) is the axial strain, and (k) is the curvature term as expressed by Semler.
The thermal effect can be introduced by considering the linear strain tensor as a sum of the strain contributions from the mechanical stress and the thermal effect. Semler further decomposed the axial strain into a steady strain component due to externally applied tension (T 0 ), a pressure force component (P = pA), and an oscillatory strain component due to the oscillations of the pipe. These can be expressed as while the stress-contributing strain component is as expressed by Semler as We note that the gradient of the transverse displacement of the pipe is far greater than the gradient of the longitudinal displacement ∂v ∂x > ∂u ∂x . Also, the thermal contributing strain component can be expressed in terms of the thermal expansivity (α) and the difference in temperatures (∆T) as ij ∆ = (−α∆T). (11) Substituting Equations (10) and (11) into Equation (9), and then substituting Equation (9) into Equation (8), the resulting expression is With the reference plane in the same direction as the gravitational acceleration, the effect of gravity can be expressed as The variations of the time derivative of the algebraic sum of the kinetic energy and the potential energy of the systems gives the expression on the left-hand side term of Hamilton's equation.

Non-Conservative Work Done
As detailed by Semler for a single-phase flow, the right-hand side term of Hamilton's equation can be expressed for a multiphase flow as Physically, this implies a non-classical boundary condition at the free end for a discharging cantilever pipe. Therefore, a force is imposed at the free end if the velocity of the exiting fluid is not tangential to the pipe. However, this study assumes that the exiting flow remains tangential to the pipe at the free end, and therefore the classical boundary condition holds at the free end.

Equation of Motion for Multiphase Flow
The equation of motion for an extensible cantilever pipe conveying pressurized unsteady multiphase flow under thermal loading can be expressed as The associated boundary conditions are

Dimensionless Equation of Motion for Two-Phase Flow
The equation of motion may be rendered dimensionless to make the analysis of the system more robust and not constrained to one specific system by introducing the following non-dimensional quantities: The dimensionless governing equation can be reduced to that of a two-phase as follows: ..
For a steady flow, velocities are not changing with time, therefore We define the following notation: The governing equation for a cantilevered pipe conveying steady two-phase flow can be expressed as ..
. The dimensionless boundary conditions are In these equations, u and v are, respectively, the dimensionless displacements in the longitudinal and transverse direction; (U j ) is the flow velocities of the constituent phases/components used in the parametric studies of the dynamics of the system; (β j ) is the same mass ratio as in single-phase flows as derived by Semler and Paidoussis [8]; (Ψ j ) is another mass ratio which is unique to multiphase flow relating the fluid mass independent of the mass of the pipe; (γ) is the gravity term; and (Π 0 , Π 1 , Π 2 ) represent the Tension term, Flexibility term, and the Pressurization term, respectively.

Empirical Gas-Liquid Two-Phase Flow Model
Adopting the Chisholm empirical relations as presented in [23], the void fraction is expressed as Volume of gas Volume of gas + Volume of Liquid .
The Slip Ratio is expressed as where x is the vapor density; and ρ l , ρ g are the densities of the liquid and gas phases, respectively. Mixture Velocity is expressed as where U g and U l are the superficial flow velocities. Individual Velocities for the two phases can be expressed as For various void fractions (0.3, 0.4, and 0.5) and a series of mixture velocities, the corresponding slip ratio and individual velocities are estimated and used for numerical calculations.

Method of Solution
Exact solutions of nonlinear equations are almost always not available; an approximate solution will be sought by utilizing the multiple time-scale perturbation technique. This approach is applied directly to the partial differential Equations (22) and (23), given that the common method of discretizing the equations first and then applying the perturbation method yields less accurate results for finite mode truncations and higher order perturbation schemes [15][16][17][18].
Adopting perturbation techniques, it is necessary to decide the terms to be considered small or weak. This study considers the contributions of the nonlinear terms, gradient term and gravity term to be small compared to the linear terms. ..
We seek an approximate solution for u and v in the form where ε is a small dimensionless measure of the amplitude of u and v, used as a book-keeping parameter. Two time scales are needed: T 0 = t and T 1 = εt. Then, the time derivatives are where D n = ∂ ∂T n . Substituting Equations (32)-(35) into Equations (30) and (31) and equating the coefficients of (ε) to zero and one, respectively, we get the following: The order zero problems for both the axial and transverse vibration of the cantilever pipe have the form of an undamped and unforced flow-induced vibration problem; this will be used to estimate the linear natural frequencies and mode shapes, while the order one problem will be solved to obtain the amplitude of the nonlinear response of the pipe.

Linear Analysis
The leading order equations present a set of linear equations which relates the flow velocity-generated forces (Coriolis and Centrifugal forces) to the stiffness of the pipe without neglecting the mass ratios. This set of linear equations exhibits the form of an eigenvalue problem, which upon resolution will produce the natural frequency and mode shapes, and also predict the stability of the system.

Natural Frequencies and Modal Functions
Estimation of the natural frequencies and modal function is an order zero problem that can be determined by resolving Equations (36) and (38).
The homogeneous solution of the leading order Equations (36) and (38) can be expressed as where (CC) is the complex conjugate, φ(x) n and η(x) n are the complex modal functions for the axial and transverse vibrations for each mode (n), and ω n and λ n are the eigenvalues for the axial and transverse vibrations for each mode (n).
The eigenvalues are complex values with a complex conjugate pair of solutions which can be expressed as The real parts of the eigenvalues are associated with the natural frequency of oscillation, and the imaginary parts with the damping.

Solution to Axial Vibration Problem
Substituting Equation (44) into Equation (42) gives a quadratic relation of the form We solve the quadratic Equation (46) for the wave numbers (k j ) as a function of the eigenvalue (ω n ): In order to obtain the eigenvalue, Equation (51) is substituted into the boundary conditions in Equation (25): and In matrix form: For a non-trivial solution, the determinant of (D) must varnish: We substitute Equations (47) and (48) into Equation (52) and solve for the eigenvalue: where a = C21U 1 2 Equation (53) is the pipe's axial vibration eigenvalue. Solving Equations (49) and (50) gives the constants G1 and G2. Therefore, the modal function for the axial vibration of the pipe is expressed as Substituting Equation (44) into Equation (40) yields It can be observed from Equation (55) that the real part is the natural frequency and the imaginary part is the amplitude. However, as the mixture velocity is varied, a critical value of the mixture velocity is attained when the imaginary parts of any of the eigenvalues (ω n ) take on a negative value which will cause the axial displacement (u) to grow exponentially in time and linearly signifies the onset of the system's oscillatory instability.
This gives four algebraic equations which can be expressed in matrix form as For a non-trivial solution, the determinant of (G) must varnish, that is In order to find modal solutions of (λ n ), Equations (57) and (58) must be solved simultaneously; this can be solved numerically using a nonlinear numerical routine.
The mode function of the transverse vibration corresponding to the nth eigenvalue is expressed as It can be observed from Equation (60) that the real part is the natural frequency and the imaginary part is the amplitude. However, as the mixture velocity is varied, a critical value of the mixture velocity is attained when the imaginary parts of any of the eigenvalues (λ n ) take on a negative value which will cause the transverse displacement (v) to grow exponentially in time and linearly signifies the onset of the system's flutter instability.

Nonlinear Analysis
Linear analysis of the system is sufficient for the prediction of the critical velocity at which instability will occur, but cannot predict the post-buckling behavior. Linear theory has shown that buckling amplitudes will grow unboundedly with time after the critical velocity, however, as amplitudes grow, the effect of nonlinearities comes into play. This predictably limits the growth to some finite value [24][25][26].

Nonlinear Axial and Transverse Vibration Problem
We solve the nonlinear axial and transverse vibration problem by seeking an approximate solution for u and v of the form The zero order solutions produce the undamped and uncoupled linear solution of the axial and transverse vibrations, respectively: where X and Y are unknown complex-valued functions of the slow time scale T 1 , (CC) is the complex conjugate, φ(x) and η(x) are the modal functions for the axial and transverse vibrations, and ω = Re(ω) and λ = Re(λ). (the real parts of the complex frequencies) are the natural frequencies for the axial and transverse vibrations. Substituting Equations (63) and (64) into Equations (37) and (39) gives where where CC and NST denote complex conjugates and non-secular terms, respectively.
The next task is to determine the requirements for X(T1) and Y(T1) that permit the solutions of u 1 . and v 1 . to be independent of secular terms. However, examining Equations (65) and (66), it can be observed that two scenarios exist: ω = 2λ and ω = 2λ.

When ω Is Far from 2λ
If ω is far from 2λ, then none of the coupled nonlinear terms will generate secular terms, therefore resulting to an uncoupled response.
Equations (65) and (66) will have bounded solutions only if the solvability condition holds. The solvability condition demands that the coefficients of exp(iωT 0 ) and exp(iλT 0 )-that is, X(T1) and Y(T1)-should satisfy the following relation: Q3· The inner product for complex functions on [0, 1] is defined as Equations (67) and (68) can be cast as where P = where P is a complex number such that We substitute Equation (73) into Equation (71) and sort the outcome into real and imaginary parts, resolving for αy and βy.

When ω Is Close to 2λ
In order to examine the coupled nonlinear dynamics of the system, which is the scenario when ω = 2λ, a detuning parameter σ is introduced.
The two Equations (65) and (66) will have bounded solutions only if the solvability condition holds. The solvability condition demands that the coefficient of epe(iωT 0 ) and exp(iλT 0 ) vanishes.
That is, X(T1) and Y(T1) should satisfy the following relation: With the inner product defined for complex functions on [0, 1] as expressed in Equation (77), the Equations (87) and (88) can be cast as To determine X(T1) and Y(T1), the solution of Equations (89) and (90) is expressed in polar form: We substitute the polar expressions into the solvability condition and separate the real and imaginary parts. The following set of modulation equations is formed: The first order nonlinear solutions can be obtained by solving for αx(T1) and αy(T1) completely. We define the following notation: CS = cos(ψ1) and SS = sin(ψ1). Now αy(T1) can be obtained by resolving the quartic frequency amplitude relation: where The solution of the quartic Equation (96) will produce four roots of αy(T1): However, the acceptable solution for αy(T1) is the root of the quartic Equation (96) that is real and positive [27][28][29]. The expression for estimating αx(T1) is defined in terms of αy(T1) as Considering the nth values of (T1), αy(T1), βx(T1), and βy(T1) corresponding to the nth modal functions and the nth natural frequencies, the nth solution of the coupled problem is expressed as v(x, t) n = αy(T1) n η(x) n cos(λ n T0 + βy(T1) n ) + O(ε).

Numerical Results
This section presents the numerical solutions of the nonlinear dynamics of a cantilever pipe conveying a steady pressurized air/water two-phase flow, summary of the pipe and flow parameters are presented in Table 1. The axial linear natural frequencies for different flow velocities are estimated analytically from Equation (53), while the transverse linear natural frequencies are estimated by solving Equations (57) and (58) simultaneously with a numeric code written in Matlab. Considering a simple system with β = 0.2 and Π1 = 100, Π0 = Π2 = 0, a = α∆T = 0 for a single-phase flow through the pipe, the natural frequencies Argand diagram plots for the axial and transverse vibrations are presented in Figures 2 and 3.  Considering a simple system with β = 0.2 and Π1 = 100, Π0 = Π2 = 0, a = αΔT = 0 for a single-phase flow through the pipe, the natural frequencies Argand diagram plots for the axial and transverse vibrations are presented in Figures 2 and 3.    The Argand diagram of the axial vibrations shows that, as the fluid mixture velocity tends towards the critical velocity, all paths move towards the origin of the Argand diagram; for the transverse vibration it can be seen that, as the mixture velocity attains higher values, the ( ) in the second mode of the system starts to diminish and, in time, becomes negative; Therefore, a Hopf The Argand diagram of the axial vibrations shows that, as the fluid mixture velocity tends towards the critical velocity, all paths move towards the origin of the Argand diagram; for the transverse vibration it can be seen that, as the mixture velocity attains higher values, the Im(λ) in the second mode of the system starts to diminish and, in time, becomes negative; Therefore, a Hopf bifurcation occurs at an approximate dimensionless velocity of 5.65, which is the critical velocity at which the systems becomes transversely unstable, summary of the results are presented in Table 2. The nonlinear behavior of the pipe will be examined for the scenario when there is coupling of the axial and transverse vibration and for the scenario when both are uncoupled. The uncoupled nonlinear transverse frequency presented in Equation (80) demonstrates a cubic nonlinear dependence of the uncoupled transverse frequencies on the amplitudes. The quartic expression presented in Equation (96) relates the frequency and the detuning parameter of the coupled axial and transverse vibration with the amplitude. These expressions will be used to plot the frequency response curves for both the uncoupled and coupled scenarios and the nonlinear behavior of the pipe as it conveys two-phase flow at a supercritical mixture velocity-"flow velocity larger than the critical velocity".

Effects of Two-Phase Flow on the Dynamic Behavior of the Pipe
The effect of two-phase flow is studied by considering the nonlinear response of the cantilever pipe as the void fraction of the two-phase flow changes. Similar to the single-phase flow, the Argand diagram of the eigen-frequencies is used to find the critical velocities of the two-phase flow for the various void fractions (0.3, 0.4, and 0.5); with the corresponding slip ratios estimated from the Chisholm empirical relations presented in Equations (26) to (29), summary of the flow parameters are presented in Table 3.
Considering a supercritical mixture velocity of 15 and book-keeping parameter (ε) of 0.1, the nonlinear frequency amplitude variations and coupled frequency response of the second mode, which is the stability-determining mode in the linear sense, are plotted for various void fractions.
It can be seen in Figure 4 that, as the void fraction increases, the natural frequency reduces. A nonlinear hardening behavior is observed in the dynamic response of the pipe for all the examined void fractions. As seen in Figure 5, as the detuning parameter is increasing, bifurcation is observed for the various void fractions examined. Also, it can be observed that the amplitude of the coupled transverse vibration is reducing as the void fraction increases.

Effects of Temperature Difference on the Dynamic Behavior
The effect of temperature differences on a cantilever pipe conveying two-phase flow is studied by considering the nonlinear response of the cantilever pipe as the temperature difference of the two-phase flow with void fraction of 0.3 changes. Similar to the single-phase flow, the Argand diagram of the eigen-frequencies is used to find the critical velocities of the two-phase flow for the various temperature differences (0, 40, and 50); with the corresponding slip ratios estimated from the Chisholm empirical relations presented in Equations (26)- (29), summary of the flow parameters are presented in Table 4.

Parameter Void Fraction Thermal Expansivity
Critical Mixture Velocity

Effects of Temperature Difference on the Dynamic Behavior
The effect of temperature differences on a cantilever pipe conveying two-phase flow is studied by considering the nonlinear response of the cantilever pipe as the temperature difference of the two-phase flow with void fraction of 0.3 changes. Similar to the single-phase flow, the Argand diagram of the eigen-frequencies is used to find the critical velocities of the two-phase flow for the various temperature differences (0, 40, and 50); with the corresponding slip ratios estimated from the Chisholm empirical relations presented in Equations (26)- (29), summary of the flow parameters are presented in Table 4.

Parameter Void Fraction Thermal Expansivity
Critical Mixture Velocity

Effects of Temperature Difference on the Dynamic Behavior
The effect of temperature differences on a cantilever pipe conveying two-phase flow is studied by considering the nonlinear response of the cantilever pipe as the temperature difference of the two-phase flow with void fraction of 0.3 changes. Similar to the single-phase flow, the Argand diagram of the eigen-frequencies is used to find the critical velocities of the two-phase flow for the various temperature differences (0, 40, and 50); with the corresponding slip ratios estimated from the Chisholm empirical relations presented in Equations (26)- (29), summary of the flow parameters are presented in Table 4. Considering a supercritical mixture velocity of 15 and book-keeping parameter (ε) of 0.1, the nonlinear frequency amplitude variations and coupled frequency response of the second mode, which is the stability-determining mode in the linear sense, is plotted for various temperature differences.
It can be seen in Figure 6 that, as the temperature difference increases, the natural frequency increases. A nonlinear hardening behavior is observed in the dynamic response of the pipe for all the examined temperature differences. As seen in Figure 7, as the detuning parameter is increasing, bifurcation is observed for the various temperature differences examined. Also, it can be observed that the increase in temperature difference did not create significant changes in the amplitude of the coupled transverse vibration. It can be seen in Figure 6 that, as the temperature difference increases, the natural frequency increases. A nonlinear hardening behavior is observed in the dynamic response of the pipe for all the examined temperature differences. As seen in Figure 7, as the detuning parameter is increasing, bifurcation is observed for the various temperature differences examined. Also, it can be observed that the increase in temperature difference did not create significant changes in the amplitude of the coupled transverse vibration.

Effects of Flow Pressure on the Dynamic Behavior
The effect of pressurization on a cantilever pipe conveying two-phase flow is studied by considering the nonlinear response of the cantilever pipe as the pressurization of the two-phase flow with void fraction of 0.3 changes. Similar to the single-phase flow, the Argand diagram of the eigen-frequencies is used to find the critical velocities of the two-phase flow for the various pressures (0, 5, and 10); with the corresponding slip ratios estimated from the Chisholm empirical relations presented in Equations (26)-(29), summary of the flow parameters are presented in Table 5. It can be seen in Figure 6 that, as the temperature difference increases, the natural frequency increases. A nonlinear hardening behavior is observed in the dynamic response of the pipe for all the examined temperature differences. As seen in Figure 7, as the detuning parameter is increasing, bifurcation is observed for the various temperature differences examined. Also, it can be observed that the increase in temperature difference did not create significant changes in the amplitude of the coupled transverse vibration.

Effects of Flow Pressure on the Dynamic Behavior
The effect of pressurization on a cantilever pipe conveying two-phase flow is studied by considering the nonlinear response of the cantilever pipe as the pressurization of the two-phase flow with void fraction of 0.3 changes. Similar to the single-phase flow, the Argand diagram of the eigen-frequencies is used to find the critical velocities of the two-phase flow for the various pressures (0, 5, and 10); with the corresponding slip ratios estimated from the Chisholm empirical relations presented in Equations (26)-(29), summary of the flow parameters are presented in Table 5.

Effects of Flow Pressure on the Dynamic Behavior
The effect of pressurization on a cantilever pipe conveying two-phase flow is studied by considering the nonlinear response of the cantilever pipe as the pressurization of the two-phase flow with void fraction of 0.3 changes. Similar to the single-phase flow, the Argand diagram of the eigen-frequencies is used to find the critical velocities of the two-phase flow for the various pressures (0, 5, and 10); with the corresponding slip ratios estimated from the Chisholm empirical relations presented in Equations (26)- (29), summary of the flow parameters are presented in Table 5. Considering a supercritical mixture velocity of 15 and book-keeping parameter (ε) of 0.1, the nonlinear frequency amplitude variations and coupled frequency response of the second mode, which is the stability-determining mode in the linear sense, is plotted for various pressures.
It can be seen in Figure 8 that, as the pressure increases, the natural frequency increases. A nonlinear hardening behavior is observed in the dynamic response of the pipe for all the examined pressures. As seen in Figure 9, as the detuning parameter is increasing, bifurcation is observed for the various pressures examined. Also, it can be observed that increasing the pressure did not create significant changes in the amplitude of the coupled transverse vibration. Considering a supercritical mixture velocity of 15 and book-keeping parameter (ε) of 0.1, the nonlinear frequency amplitude variations and coupled frequency response of the second mode, which is the stability-determining mode in the linear sense, is plotted for various pressures.
It can be seen in Figure 8 that, as the pressure increases, the natural frequency increases. A nonlinear hardening behavior is observed in the dynamic response of the pipe for all the examined pressures. As seen in Figure 9, as the detuning parameter is increasing, bifurcation is observed for the various pressures examined. Also, it can be observed that increasing the pressure did not create significant changes in the amplitude of the coupled transverse vibration.  The effect of top tension on a cantilever pipe conveying two-phase flow is studied by considering the nonlinear response of the cantilever pipe for a situation with no top tension, a Considering a supercritical mixture velocity of 15 and book-keeping parameter (ε) of 0.1, the nonlinear frequency amplitude variations and coupled frequency response of the second mode, which is the stability-determining mode in the linear sense, is plotted for various pressures.
It can be seen in Figure 8 that, as the pressure increases, the natural frequency increases. A nonlinear hardening behavior is observed in the dynamic response of the pipe for all the examined pressures. As seen in Figure 9, as the detuning parameter is increasing, bifurcation is observed for the various pressures examined. Also, it can be observed that increasing the pressure did not create significant changes in the amplitude of the coupled transverse vibration.  The effect of top tension on a cantilever pipe conveying two-phase flow is studied by considering the nonlinear response of the cantilever pipe for a situation with no top tension, a tensioning value of 5, and a compressing value of 5, with a two-phase flow and void fraction of 0.3.

Effects of Top Tension on the Dynamic Behavior
The effect of top tension on a cantilever pipe conveying two-phase flow is studied by considering the nonlinear response of the cantilever pipe for a situation with no top tension, a tensioning value of 5, and a compressing value of 5, with a two-phase flow and void fraction of 0.3. Similar to the single-phase flow, the Argand diagram of the eigen-frequencies is used to find the critical velocities of the two-phase flow for the various tensioning values, and the corresponding slip ratios are estimated from the Chisholm empirical relations presented in Equations (26)- (29), summary of the simulated parameters are presented in Table 6.
Considering a supercritical mixture velocity of 15 and book-keeping parameter (ε) of 0.1, the nonlinear frequency amplitude variations and coupled frequency response of the second mode, which is the stability-determining mode in the linear sense, is plotted for the top tensions. It can be seen in Figure 10 that the tensioning top load reduces the natural frequency increase, while a compressing top load increases the natural frequency. A nonlinear hardening behavior is observed in the dynamic response of the pipe for all the examined cases. As seen in Figure 11, as the detuning parameter is increasing, bifurcation is observed for the various pressures examined. Also, it can be observed that the top tensions did not create significant changes in the amplitude of the coupled transverse vibration.  It can be seen in Figure 10 that the tensioning top load reduces the natural frequency increase, while a compressing top load increases the natural frequency. A nonlinear hardening behavior is observed in the dynamic response of the pipe for all the examined cases. As seen in Figure 11, as the detuning parameter is increasing, bifurcation is observed for the various pressures examined. Also, it can be observed that the top tensions did not create significant changes in the amplitude of the coupled transverse vibration.

Time History and Phase Plots
Equations (83)

Time History and Phase Plots
Equations (83) The uncoupled response of the transverse vibrations-as shown in Figures 12-14 for the various void fractions-looks similar; they all reveal that the uncoupled transverse vibrations exhibit an oscillation that converges to a limit cycle with time. The initial amplitudes are greater than the amplitudes of the limit cycles; hence, a positive initial damping is observed and the amplitudes decay until they attain the limit cycles, as shown in the phase plot with a set of concentric circles inside the phase trajectories. Contrary to this, the uncoupled axial vibration-as shown in Figures 15-17-is observed to exhibit uniform periodic oscillations, which traces out as a closed orbit in the phase plots. The amplitudes of the displacements are observed to reduce as the void fraction increases.
As a result of the coupling between the axial and the transverse vibration, both the coupled transverse and axial solutions-as shown in Figures 18-23-exhibit uniform periodic oscillations, which trace out as a closed orbit in the phase plots. The amplitudes of the coupled displacements are observed to reduce as the void fraction increases. exhibit an oscillation that converges to a limit cycle with time. The initial amplitudes are greater than the amplitudes of the limit cycles; hence, a positive initial damping is observed and the amplitudes decay until they attain the limit cycles, as shown in the phase plot with a set of concentric circles inside the phase trajectories. Contrary to this, the uncoupled axial vibration-as shown in Figures  15-17-is observed to exhibit uniform periodic oscillations, which traces out as a closed orbit in the phase plots. The amplitudes of the displacements are observed to reduce as the void fraction increases.
As a result of the coupling between the axial and the transverse vibration, both the coupled transverse and axial solutions-as shown in Figures 18-23-exhibit uniform periodic oscillations, which trace out as a closed orbit in the phase plots. The amplitudes of the coupled displacements are observed to reduce as the void fraction increases.

Conclusions
This study examines the dynamic behavior of a cantilever pipe conveying two-phase flow. Taking into consideration the extensible theory, nonlinear equations of motion and boundary conditions were obtained using Hamilton's principle. The equations were made to be non-dimensional so as to remove the dependence on geometric and dimensional parameters. Using the multiple-scale perturbation technique, natural frequencies, mode shapes, and first order approximate solutions of the steady state response of the pipes were obtained. From the multiple-scale assessment, it was observed that a 1:2 coupling exists between the axial and the transverse vibrations of the pipe.
The critical flow mixture velocities for various void fractions were obtained from the Argand diagram plot of the eigen-frequencies; it was observed that the critical velocity increases as the void fraction increases. The investigation on the uncoupled nonlinear dynamic behavior of the pipe as it conveys two-phase flow at a super-critical mixture velocity reveals that the system exhibits a nonlinear hardening behavior. As a result of the dynamic analysis, it has been observed that for a two-phase flow, an increase in the void fraction reduces the natural frequency and the coupled amplitude of the system. Also, an increase in temperature difference, an increase in pressure, and the presence of top tension were observed to increase the natural frequencies without a significant alteration in the coupled amplitude of the system, while compression load at the top was observed to reduce the natural frequencies without a significant change in the coupled amplitude of the system.

Conclusions
This study examines the dynamic behavior of a cantilever pipe conveying two-phase flow. Taking into consideration the extensible theory, nonlinear equations of motion and boundary conditions were obtained using Hamilton's principle. The equations were made to be non-dimensional so as to remove the dependence on geometric and dimensional parameters. Using the multiple-scale perturbation technique, natural frequencies, mode shapes, and first order approximate solutions of the steady state response of the pipes were obtained. From the multiple-scale assessment, it was observed that a 1:2 coupling exists between the axial and the transverse vibrations of the pipe.
The critical flow mixture velocities for various void fractions were obtained from the Argand diagram plot of the eigen-frequencies; it was observed that the critical velocity increases as the void fraction increases. The investigation on the uncoupled nonlinear dynamic behavior of the pipe as it conveys two-phase flow at a super-critical mixture velocity reveals that the system exhibits a nonlinear hardening behavior. As a result of the dynamic analysis, it has been observed that for a two-phase flow, an increase in the void fraction reduces the natural frequency and the coupled amplitude of the system. Also, an increase in temperature difference, an increase in pressure, and the presence of top tension were observed to increase the natural frequencies without a significant alteration in the coupled amplitude of the system, while compression load at the top was observed to reduce the natural frequencies without a significant change in the coupled amplitude of the system. Author Contributions: Adeshina Adegoke and Ayo Oyediran conceived and formulated the model of the system; Adeshina Adegoke performed the resolution of the problem; Adeshina Adegoke and Ayo Oyediran analyzed the results; Ayo Oyediran contributed materials/analysis tools; Adeshina Adegoke wrote the paper.

Conflicts of Interest:
The authors declare no conflict of interest.