The research for partitioned algorithm of coupling of elastic pipe and cross-flow in a cylinder

Qian Bei Yue1 , Ri Zhi Dong2 , Ju Bao Liu3

1, 2, 3College Mechanical Science Engineering, Northeast Petroleum University, Hei long-jiang, Daqing, 163318, China

1Corresponding author

Journal of Vibroengineering, Vol. 19, Issue 5, 2017, p. 3754-3767. https://doi.org/10.21595/jve.2017.17731
Received 14 September 2016; received in revised form 26 December 2016; accepted 19 January 2017; published 15 August 2017

Copyright © 2017 JVE International Ltd. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Creative Commons License
Table of Contents Download PDF Acknowledgements References
Cite this article
Views 31
Reads 9
Downloads 811
WoS Core Citations 0
CrossRef Citations 0
Abstract.

In the marine flatbed, nuclear power stations and oil fields, the marine risers and heat exchangers are induced vibrations and collision, even caused failure with cross flow. This is a typical problem of bundles vibration in fluid-structure interaction. In this paper, the elastic pipe was separated into several beam elements. And the dynamical equation was adopted to describe the structure. The fluid domain of cylinder was separated into solid elements. And the computational fluid dynamics equation was adopted to describe the fluid. At the fluid-structure coupling interface, the formulas of displacement, velocity and load was derived and the corresponding convergence criterion was also derived. Thus, the partitioned coupling algorithm was established for the elastic pipe and the cross-flow in cylinder. The example showed that, comparing to the results of monolithic coupling algorithm, the error was small. Thus, it proved that the correctness of the partitioned coupling algorithm in this article. And it will provide an effective computational method for vibrations and collision of bundles in a more complex fluid domain.

Keywords: partitioned coupling algorithm, monolithic coupling algorithm, elastic pipe, fluid-structure interaction.

1. Introduction

In the marine flatbed, nuclear power stations and oil fields, the marine risers (seen as Fig. 1) and heat exchangers [1, 2] are induced vibrations and collision, even caused failure with internal flow and external flow. This is a typical problem of bundles vibration in structure-interaction coupling. It related to multiple movable interfaces and multiple physical fields and which was difficult to calculate. Most of the results were studied the tube vibrations with internal flow and the vibration law was analyzed in longitudinally. M. Ragulskis [3] adopted harmonic balance method to study the non-Newtonian fluid flow in longitudinally vibrating tube. The tube vibrations induced by internal flow and the tube was a non-deformable rigid body. So far, the numerical methods for the tube vibrations with cross flow (external flow) are monolithic coupling algorithm and partitioned coupling algorithm [4-7]. In monolithic coupling algorithm, equations are solved simultaneously for fluid and structural domain. The advantage is the high accuracy [8, 9]. The disadvantage is that because the equations are integrated for fluid and structure and the solution efficiency is very low. Especially, it is not suitable for problems of large-scale and more complex domain. In partitioned coupling algorithm, equations are solved respectively. Within each time step or iteration, the information is transferred at the coupling interface. And the coupled solution is achieved. The advantage is that the solving range is wide and easy to conduct. The disadvantage is that the difference of grid density for fluid domain and structural domain at the coupling interface and which causes the grid do not match [10, 11]. Therefore, it need data transfer at the coupling interface. And it is the first to solve the problem in partitioned coupling algorithm. Comparing to the results of monolithic coupling algorithm, the error was small. The research outcome will make commutations of vibration and collision for bundles possible in a more complex fluid domain.

Fig. 1. The tubes in engineering

The tubes in engineering

a) The marine risers in the ocean

The tubes in engineering

b) The heat exchanger tubes

2. The mechanics model of the coupling of elastic pipe and cross-flow in cylinder

The fluid and elastic pipe was selected as research object in cylinder. The coupling mechanics model of the elastic pipe was set up. As shown in Fig. 2, the model was divided into two domains: the structural domain s and the fluid domain Ωf. The assumptions were used as followed.

(1) The fluid medium is water. And it assumes the water is sticky and incompressible.

(2) The solid is assumed to be perfectly elastic body. And the elastic tube is prismatic. The cross section is always circular.

As shown in Fig. 1, the boundary of cylindrical entrance was velocity. And it was denoted by ΓDf(1). The boundary of outlet was pressure and it was denoted by ΓNf. The boundary that the fluid was flowing along the inner cylinder’s wall surface was free slip and it was denoted by ΓDf(2). The interface between the fluid and solid interaction was fluid-structure coupling boundary and it denoted by ΓFSI .

The boundary conditions for the cylinder fluid domain Ωf were:

(1)
d ˙ r f = n r f u ,           Γ D f 1 , d ˙ r f = 0 ,           Γ D f ( 2 ) , p r = n r f p 0 ,         Γ N f ,

where: nrf, d˙rf, pr were the vectors of normal unit, velocity, pressure, respectively. u was inlet velocity and p0 was outlet pressure. The superscripts of f,s were represented fluid, solid respectively.

The structural domain Ωs was the elastic and the boundary was hinged on both ends.

Fig. 2. The mechanics model of the elastic pipe and the cross-flow in cylinder

 The mechanics model of the elastic pipe and the cross-flow in cylinder

Supposed that the fluid and solid surface without sliding, the displacement and speed should meet the condition of coordination and the load should meet the condition of load balancing.

1) The condition of displacement coordination:

(2)
d   s n s + d   s n f = 0 ,

where: ds, df were the displacement vectors of solid and fluid respectively. ns, nf were the unit displacement vectors of solid and fluid domain respectively. ns=nrs+nθs+nzs,nf=nrf+nθf+nzf. nrs, nθs, nzs were unit vectors along the r, θ, z direction.

2) The condition of velocity coordination:

(3)
d ˙ s n   s + d ˙ f n   f = 0 ,

where: d˙f, d˙f were the velocity vectors of solid and fluid respectively. n s, n f were the same as the Eq. (2).

3) The condition of load balancing:

(4)
σ s n s + σ f n f = 0 ,

where: σs, σf were the stress vectors of solid and fluid respectively.

3. The monolithic algorithm of the coupling of elastic pipe and the cross-flow in cylinder

3.1. The numerical model of the coupling of elastic tube and cross-flow with monolithic algorithm

In monolithic algorithm, the fluid and were both divided by solid elements. The numerical model was shown in Fig. 3. The solid elements were adopted eight-node hexahedral element or six-node pentahedral element, as shown in Fig. 4.

Fig. 3. The numerical model of transverse fluid and elastic tube with the method of universe

 The numerical model of transverse fluid and elastic tube with the method of universe

a) The finite element of cylindrical fluid

 The numerical model of transverse fluid and elastic tube with the method of universe

b) The finite element model of elastic tube

Fig. 4. Solid element

 Solid element

a) Hexahedral element

 Solid element

b) Pentahedral element

3.2. The monolithic algorithm of the coupling of elastic tube and cross-flow

The fonctionelle of generalized variational principles was as followed in fluid-solid coupled system:

(5)
Π = Π 0 + Q c + Q L ,

where, Π0 was the fonctionelle of generalized variational principles in fluid and structural domain. Qc and QL were constraints and external work which cited by the mouthed of Lagrange multiplier in fluid-solid coupled system.

As mentioned above, the fluid and structural domain were disserted by eight-node hexahedral element or six-node pentahedral element. Then the fonctionelle of generalized variational principle was obtained as followed for discrete model of fluid-solid coupling:

(6)
Π = i = 1 N f Π f ( e ) + i = 1 N s Π s e + i = 1 M c Q c e + i = 1 M c Q L e ,

where, Nf and Ns were the numbers of discrete element in fluid and structural domain respectively. Mc was the number of nodes on the fluid-structure coupling interface.

Then the finite element equations and constraints were written with matrix form. Then, the dynamic coupling finite element equation was:

(7)
A α β 0 0 0 0 M 0 0 0 0 0 0 0 0 0 0 d ¨   f d ¨   s 0 0 + B α β γ j u ˙ α j   + D α β + C 11 c C 12 c C 13 c C α β   C 21 c C + C 22 c C 23 c 0 C 31 c C 32 c C 33 c 0 G α β   0 0 0 d ˙   f d ˙   s λ p
            + 0 0 0 0 0 K 0 0 0 0 0 0 0 0 0 0 d   f d   s 0 0 = F f   F s   a I H ,

where, Cijc was the sub-block of fluid-structure coupling matrix. The remaining letter meaning was seen in Ref. [12, 13].

4. The partitioned algorithm of the coupling of elastic pipe and the cross-flow in cylinder

The monolithic algorithm is more accurate, the calculation is easier to implement. But the fluid and elastic tube are dispersed by solid element and the model is large which leads the low efficiency in monolithic algorithm. In particular, it is difficult to solve the coupling dynamics problem of tube bundle vibration and collision induced by cross-flow. Based on this, the paper presents a partitioned algorithm for calculating the vibration of elastic tube induced by cross-flow in a cylinder.

In monolithic algorithm, if we do not consider the coupling of cross-flow and elastic tube, the Eq. (7) can be simplified as:

(8)
A α β d ¨   f + B α β γ j u ˙ α j   + D α β d ˙   f + C α β   p = F f   , G α β   d ˙   f = H ,
(9)
M d ¨   s + C d ˙   s + K d   s = F s .

The Eq. (8) is the control equation for fluid domain. The Eq. (9) is the dynamic equation for solid domain of the elastic tube. In partitioned algorithm, it is not necessary to solve the fluid-solid coupling Eq. (7). According to the set order, the Eqs. (8) and (9) are solved respectively. The information which calculated in fluid domain and solid domain are transferred to each other at fluid-solid interface (FSI) [14, 15]. Until the convergence of this time to meet the requirements, then the next moment of the calculation is carried. The final result will be obtained in repeated cycle.

4.1. The numerical model of the coupling of elastic tube and cross-flow with partitioned algorithm

The numerical model of the coupling of elastic tube and cross-flow was seen as Fig. 5 with partitioned algorithm. Considering the slender beam characteristics of the elastic tube, the elastic tube was divided into several beam elements to improve the calculation scale and efficiency in partitioned algorithm. The coupling interface between fluid domain and solid domain of the elastic tube was shown in Fig. 6.

It can be seen from Fig. 5 and Fig. 6 that, due to the mismatch between solid element and elastic tube beam element at the coupling interface, the transfer information of displacement, velocity and force was deduced and partitioned algorithm for the coupling of cross-flow and elastic tube in the cylinder is achieved.

Fig. 5. The numerical model of transverse fluid and elastic tube with partitioned method

 The numerical model of transverse fluid and elastic tube with partitioned method

a) The finite model element of fluid domain

 The numerical model of transverse fluid and elastic tube with partitioned method

b) The finite element model of elastic tube

Fig. 6. The coupling interface

 The coupling interface

4.2. The information transfer between the elastic pipe and cross-flow in cylinder

The Fig. 7 showed the transitive relationship between the elastic pipe and cross-flow in cylinder at the coupling interface.

Fig. 7. The transitive relationship between the elastic pipe and cross-flow

 The transitive relationship between the elastic pipe and cross-flow

The red arrow from the fluid domain to structural domain was represented shear stress and pressure. Namely the fluid affected the elastic pipe by pressure, shear stress. The red arrow from the structural domain to fluid domain was represented displacement and velocity speed. Namely the elastic pipe transmitted displacement and velocity to fluid domain.

4.2.1. The information transmission of displacement and velocity at fluid-solid coupling interface

The information transmission of displacement and velocity was shown in Fig. 8. The beam element named e was taken and its nodes were i and j. The vectors of nodal displacement and velocity in local coordinate system were as followed:

(10a)
d e s = [ u i x s u   i y s u   i z s θ   i x s θ i y s θ i z s u   j x s u   j y s u   j z s θ   j x s θ   j y s θ   j z s ] Τ ,
(10b)
d ˙ e s = [ u ˙ i x s u ˙   i y s u ˙   i z s θ ˙   i x s θ ˙ i y s θ ˙ i z s u ˙   j x s u ˙   j y s u ˙   j z s θ ˙   j x s θ ˙   j y s θ ˙   j z s ] Τ .

Fig. 8. Force, displacement and velocity interpolation at coupling interface

 Force, displacement and velocity interpolation at coupling interface

a) The fluid element and fluid-solid coupling interface

 Force, displacement and velocity interpolation at coupling interface

b) The information transmission

The displacement and velocity of arbitrary axis position xm on the beam element was as followed by shape function N:

(11a)
d e , m s = u x , m s u y , m s u z , m s θ x , m s θ y , m s θ z , m s Τ = N   x = x m d e s ,
(11b)
d ˙ e , m s = d ˙ x , m s d ˙ y , m s d ˙ z , m s θ ˙ x , m s θ ˙ y , m s θ ˙ z , m s Τ = N   x = x m d ˙ e s .

According to the assumption that the cross section was always circular, the displacement components of the node k at axis position xm was as followed:

(12)
u x , m f = u x , m s + r s i n θ y , m s s i n θ k - r s i n θ z , m s c o s θ k , u y , m f = u y , m s + r c o s θ x , m s + θ k - r c o s θ k , u z , m f = u z , m s + r s i n θ x , m s + θ k - r s i n θ k ,

where, r was outer radius of the elastic pipe. xm was distance from the node i. θk was circumferential angle of node k, which was shown in Fig. 8.

The displacement vector de,mf can be expressed as:

(13)
d e , m f = u x , m f u y , m f u z , m f Τ = H F S I e d e , m s = H F S I e N   x = x m d e s ,

where, HFSIe was the displacement transformation matrix from structural domain to fluid domain.

Assigned HΓFSIe=HFSIeN x=xm, then the displacement, velocity vector of the fluid units on the fluid-solid coupling interface can be expressed as:

(14a)
d e , m f = H Γ F S I e d e s ,
(14b)
d ˙ e , m f = H Γ F S I e d ˙ e s .

Thence, the displacement, velocity vector of the fluid at time t can be expressed as:

(15a)
d e , m f ( t ) = H Γ F S I e ( t ) d e s ( t ) ,
(15b)
d ˙ e , m f ( t ) = H Γ F S I e ( t ) d ˙ e s ( t ) .

When the elastic pipe was deformed, according to the Eq. 15(a) and 15(b), the displacement and velocity would be re-applied on the boundary wall in fluid domain. Therefore, the coordination of displacement and velocity was achieved at the fluid-solid coupling interface.

4.2.2. The information transmission of load at the fluid-solid coupling interface

As shown in Fig. 9, the interface loads were obtained at the outer surface of the elastic pipe in fluid domain. And the loads were decomposed into normal pressure p, tangential shear stress τrϕ and axial shear stress τrx. First, loads were integrated along the surface and simplified at the corresponding axis. Then according to the principle of virtual work, loads were equivalent to the values of nodes i, j of beam element.

Fig. 9. The interface loads at the outer surface of the elastic pipe in fluid domain

 The interface loads at the outer surface of the elastic pipe in fluid domain

a) Space

 The interface loads at the outer surface of the elastic pipe in fluid domain

b) Plane

In the beam element at x=xm, the length of element was Lxm.The forces and torques was simplified was as followed:

(16)
q F x m = 0 L x m 0 2 π τ r θ s i n θ + P c o s θ d θ d z , q F y m = 0 L x m 0 2 π τ r θ c o s θ + P s i n θ d θ d z , q F z m = 0 L x m 0 2 π τ r z d θ d z , q m x m = r 0 L x m 0 2 π τ r θ s i n θ d θ d z , q m y m = r 0 L x m 0 2 π τ r θ c o s θ d θ d z , q m z m = r 0 L x m 0 2 π τ r z d θ d z .

The Eq. (16) was represented by a vector:

(17)
q m f = [ q F x m q F y m q F z m q m x m q m y m q m z m ] T .

According to the principle of virtual work, the above-mentioned axial loads were further equivalent to values of the node on beam elements:

(18)
F e f = m = 1 n f o x F f i = m = 1 n f o x T T 0 L x m N T q m f d x ,

where, T was the converted matrix of coordinates.

All nodal loads which were equivalent to the beam element were assembled. Then, the equivalent nodal loads Ff(t) were obtained from the fluid domain to the elastic pipe. Thus, it met the conditions of load balancing at the coupling interface.

Substituted Ff(t) into the dynamic Eq. (9) of the elastic pipe, the dynamic equations of the elastic pipe was obtained as followed when considering the effect of fluid:

(19)
M d ¨   s + C d ˙   s + K d   s = F s + F f ,

where, the Ff was the vector of fluid force at the coupling interface.

4.3. The convergence criterion of coupling of elastic pipe and cross-flow in the cylinder

In the coupling, iterative step of fluid domain Ωf(t) and solid domain Ωs(t), the normalized convergence criterion at coupling interface at any one time step was as followed:

(20)
ε Γ i j k t = ξ Γ i n e w t - ξ Γ i j k t 2 ξ Γ i n e w t 2 < ε Γ i 0 ,         i = f , s ,

where, t was the current time step. jk was the iteration number at the current time step. i=f, s were represented the fluid-solid coupling interface of fluid domain and structural domain respectively.

The physical quantity of fluid-structure coupling interface was ξΓi (t)=dΓi (t)uΓi (t)FΓi (t)Τ. Where, dΓi (t),uΓi (t) and FΓi (t) were vectors of displacement, velocity and force at the coupling interface in fluid domain or structural domain. ξΓijk(t) was the physical quantity which was applied at current iteration step. ξΓinew(t) was the new physical quantity. εΓijk(t) was the normalized convergence value. εΓijk(t)=εdΓijkεuΓijkεFΓijkΤ. Where, εΓi0=εdΓi0εuΓi0εFΓi0Τ,which was the convergence tolerance.

4.4. The frame diagram of partitioned algorithm for elastic pipe and cross-flow in cylinder

The Fig. 10 was the total frame diagram of partitioned algorithm. tend was the total computational time. The iterative process was as followed.

(1) The models of fluid domain and structural domain were established. Give inlet velocity u and outlet pressure P. Give the initial displacement and initial velocity at the coupling interface.

(2) Solve the fluid domain Ωf and get the coupling interface load FΓfjk(t);

(3) According to the Eq. (18), get the coupling interface loads F(t) of structural domain;

(4) Solve structural domain Ωs, get displacement and velocity at coupled interface in structural domain.

(5) According to the Eq. (15), get the displacement dΓsjk(t) and velocity uΓsjk(t) at coupled interface in fluid domain.

(6) According to the convergence Eq. (20), if the values of displacement, velocity and load do not meet the convergence criteria, return the coupled iteration from the Eq. (3) to Eq. (6).

(7) If the values of displacement, velocity and load meet the convergence Eq. (20) at the same time, Order t=t+Δt, follow the previous coupling step Eq. (2) to (6). Otherwise, end of the iteration.

The coupling iteration loop at the current time step was seen as the Fig. 11.

In partitioned algorithm of the elastic pipe and the cross-flow in cylinder, the data of displacement, velocity and load was transferred with the time. The process of time-step cycle was seen in Fig. 12.

Fig. 10. The total frame diagram of partitioned algorithm

 The total frame diagram  of partitioned algorithm

Fig. 11. The calculation process at the current time step

The calculation process at the current time step

Fig. 12. The process of time-step cycle

 The process of time-step cycle

Seen from the Fig. 12, the results files f resf(t), f ress(t) of fluid domain and structural domain were used as initial values at time t. Then solving the fluid domain and structural domain at time t+Δt, the results files f resf(t+Δt), f ress(t+Δt) were obtained. And the FΓf (t+Δt) in fluid domain, the dΓf (t+Δt) and uΓf (t+Δt) in structural domain were taken at the same time. Thus, these physical quantities would be initial values at time t+Δt.

5. Examples and verification

5.1. Cases and calculation parameters

The calculation parameters of the elastic pipe and fluid were listed in Table.

Table 1. The calculation parameters of the elastic pipe and fluid

Length (m)
0.5
Outer diameter (m)
0.02
Inner diameter (m)
0.016
Elastic modulus (GPa)
3.5
Poison’s ratio
0.3
The density of elastic pipe (kg/m3)
2250
The density of fluid(kg/m3)
1000
Dynamic viscosity (mPa.s)
1
Inlet velocity (m/s)
2.236
outlet pressure (MPa)
0
Natural frequency of elastic pipe (Hz)
23.91
Calculation time step (s)
0.001

Fig. 13. The path I and cross section II in the model

 The path I and cross section II in the model

As shown in Fig. 1, the elastic was placed at the center in cylinder. The elastic pipe in the cylinder fluid was calculated by the monolithic and partitioned algorithm respectively. Then the vortex-induced vibration of the elastic was obtained in solid domain. And the distribution of the pressure and velocity were also obtained in fluid domain. The corresponding results were compared by the two calculation methods. The path I and cross section II of the elastic pipe was shown in Fig. 13. The contrast results of the elastic pipe in solid domain were shown in Fig. 14 and Fig. 15 and the values were listed in Table 2 along the path I. The contrast results in fluid domain were shown in Fig. 16 and Fig. 17 and the values were listed in Table 3 at cross Section 2.

5.1.1. The contrast results of elastic pipe in solid domain.

Seen from Fig. 14(a), (b) and (c), the curves of the acceleration, velocity, displacement were identical by the methods of monolithic and partitioned for the elastic pipe. From Table 2, the maximum error of the acceleration amplitude was 5.33 % at time 0.0s – 0.3s. And the maximum error of the velocity amplitude was 4.90 % for the elastic pipe. And the maximum error of the displacement amplitude was 4.20 %. All error is less than 5 %. Hence, it presented the correctness of the partitioned algorithm for the coupling of the elastic pipe and the cross-flow in cylinder.

Fig. 14. The contrast curves of acceleration, velocity, displacement for the elastic pipe

 The contrast curves of acceleration, velocity, displacement for the elastic pipe

a) Acceleration

 The contrast curves of acceleration, velocity, displacement for the elastic pipe

b) Velocity

 The contrast curves of acceleration, velocity, displacement for the elastic pipe

c) Displacement

Table 2. The amplitude list of the acceleration, velocity, displacement

Time period (s)
0-0.03
0.03-0.06
0.06-0.09
0.09-0.12
0.12-0.15
Acceleration
Amplitude
(m/s2)
Monolithic
4006
2267
1615
1101
771
Partitioned
4015
2272
1632
1118
786
Error (%)
0.22
0.22
1.05
1.54
1.94
Velocity
Amplitude
(mm/s)
Monolithic
2011
1130
790
555
395
Partitioned
2079
1152
798
561
401
Error (%)
3.38
1.95
1.01
1.08
1.52
Displacement
Amplitude
(mm)
Monolithic
7.32
5.59
4.03
2.99
2.09
Partitioned
7.90
5.79
4.23
3.05
2.19
Error (%)
7.92
3.58
4.96
2.01
4.78
Time period (s)
0.15-0.18
0.18-0.21
0.21-0.24
0.24-0.27
0.27-0.30
Acceleration
Amplitude
(m/s2)
Monolithic
514
348
251
197
150
Partitioned
521
352
258
203
158
Error (%)
1.36
1.15
2.79
3.05
5.33
Velocity
Amplitude
(mm)
Monolithic
235
186
132
102
82
Partitioned
239
189
136
107
85
Error (%)
1.70
1.61
3.03
4.90
3.66
Displacement
Amplitude
(m/s2)
Monolithic
1.86
1.19
0.70
0.60
0.39
Partitioned
1.92
1.24
0.72
0.62
0.42
Error (%)
3.23
4.20
2.86
3.33
7.69

5.1.2. The contrast results of pressure and velocity distribution in fluid domain

The distribution of pressure and velocity was shown as Fig. 15 along path I in fluid domain at t= 0.281 s.

Seen from Fig. 15, the distribution of pressure and velocity was identity by the algorithm of monolithic and partitioned coupling along path I. Of note, in 48°-132°, namely in the position of upstream face, the fluid pressure was positive, greater than atmospheric pressure and the maximum value was at the position of 90°. In 132°-210°, 330°-360°, 0°-48°, namely in the position of side face, the fluid pressures were relatively negative, less than atmospheric pressure and the maximum negative value was occurred in 174° and 6° position. In 220°-320°, namely in the position of lee face, the fluid pressures were relatively positive, greater than atmospheric pressure. But the values were smaller than values of upstream face.

Seen from the Table 3, the errors of pressure calculated by the two methods were –0.17 %, 0.04 %-0.87 %, 0.15 % respectively, in the location of 6°, 90°, 174°, 27° along path I. And the errors of velocity were –0.27 %, 1.71 %, –0.27 %, –1.35 % respectively. Therefore, the values of pressure and velocity obtained by the two methods were equivalent in the fluid domain.

Fig. 15. The pressure and velocity distribution along path I

 The pressure and velocity distribution along path I

a) Pressure

 The pressure and velocity distribution along path I

b) Velocity

Table 3. The list of the pressure and velocity along path I

Positions (°)
Pressure
Velocity
Monolithic (Pa)
Partitioned (Pa)
Error (%)
Monolithic (m/s)
Partitioned (m/s)
Error (%)
6
–4589.70
–4581.91
–0.17
–3.73
–3.72
–0.27
90
3896.77
3898.17
0.04
–0.115
–0.117
1.71
174
–4590.43
–4550.63
–0.87
–3.68
–3.67
–0.27
270
1552.97
1555.29
0.15
0.075
0.074
–1.35

Fig. 16. The pressure and velocity contours at cross section II by monolithic algorithm

 The pressure and velocity contours at cross section II by monolithic algorithm

a) Pressure

 The pressure and velocity contours at cross section II by monolithic algorithm

b) Velocity

Fig. 17. The pressure and velocity contours at cross Section 2 by partitioned algorithm

 The pressure and velocity contours at cross Section 2 by partitioned algorithm

a) Pressure

 The pressure and velocity contours at cross Section 2 by partitioned algorithm

b) Velocity

The Fig. 16 and Fig. 17 were pressure and velocity contours at cross Section 2 by monolithic and partitioned algorithm. From the contours, the distribution of pressure and velocity at cross Section 2 were basically the same.

6. Conclusions

1) The coupling mechanics model was established for the elastic pipe and the cross-flow in cylinder. The boundary conditions of the inlet flow rate, outlet pressure and wall were considered in the model. And it also met the conditions of displacement, velocity and load balancing at fluid-structure coupling interface.

2) In monolithic algorithm, the model of fluid and elastic tube is both dispersed by solid elements. The dynamic coupling finite element equation was applied to solve the dynamic characteristics of elastic pipe in solid domain and the velocity and pressure in fluid domain. Although the monolithic algorithm is easy and accurate, the model is very lager in general and it may lead a low efficiency. In particular, it is difficult to solve the coupling dynamics problem of multiple bundles vibration and collision induced by cross-flow in a more complex fluid domain. Based on this, the partitioned algorithm was established for elastic pipe and cross-flow in cylinder. In partitioned algorithm, the elastic pipe was separated into several beam elements. The cross-flow was still separated into solid elements. At the fluid-structure coupling interface, the formulas of displacement, velocity and load was derived and the corresponding convergence criterion was also derived. The iterative block diagram of partitioned coupling algorithm was designed.

3) According the above two methods, the example was calculated by the algorithm of monolithic and partitioned. The results proved that the curves of acceleration, velocity, displacement for elastic pipe were consistent by two algorithms in solid domain. And pressure and velocity distribution were also consistent by two algorithms at cross section II in fluid domain. It proved the correctness of the partitioned coupling algorithm which was established in this article.

4) The method will provide an effective computational method for vibrations and collision of multiple bundles in a more complex fluid domain. And it will have important theoretical significance in engineering.

Acknowledgements

This work was supported by the Project of National Natural Science Foundation of China (11272085).

References

  1. Tang You Gang, Fan Juan Juan, Zhang Jie In line and transverse vortex-induced vibration analysis for a circular cylinder under high Reynolds number. Journal of Vibration and Shock, Vol. 36, Issue 13, 2013, p. 88-92. [Search CrossRef]
  2. Ji Jia Dong, Ge Pei Qi, Bi Wen Bo Numerical analysis on the combined flow induced vibration response of elastic to tube bundle in heat exchange. Journal of Xi’an Jiao Tong University, Vol. 45, Issue 4, 2015, p. 24-29. [Search CrossRef]
  3. Ragulskis M., Fedaravicius A., Ragulskis K. Harmonic balance method for FEM analysis of fluid flow in a vibrating pipe. Communications in Numerical Methods in Engineering, Vol. 22, Issue 4, 2006, p. 347-356. [Search CrossRef]
  4. Heil M., Hazel Boyle A. L. J. Solvers for large-displacement fluid–structure interaction problems: segregated versus monolithic approaches. Computational Mechanics, Vol. 43, Issue 1, 2008, p. 91-101. [Search CrossRef]
  5. Yan Ke, Ge Peiqi, Bi Wenbo, et al. Vibration characteristics of fluid-structure interaction of conical spiral tube bundle. Journal of Hydrodynamics, Vol. 22, Issue 1, 2010, p. 121-128. [Search CrossRef]
  6. Antoci Carla, Gallati Mario, Sibilla Stefano Numerical simulation of fluid-structure interaction by SPH. Computers and Structures, Vol. 85, Issues 11-14, 2007, p. 879-890. [Search CrossRef]
  7. Zhang A., Suzuki K. A comparative study of numerical simulmions for fluid-structure interaction of liquid-filled tank during ship collision. Ocean Engineering, Vol. 34, Issue 5, 2007, p. 645-652. [Search CrossRef]
  8. Ali Eken, Mehmet Sahin A parallel monolithic algorithm for the numerical simulation of large-scale fluid structure interaction problems. International Journal for Numerical Methods in Fluids, Vol. 80, Issue 12, 2016, p. 687-714. [Search CrossRef]
  9. Bathe K. J., Zhang H. Finite element developments for general fluid flows with structural interactions. International Journal for Numerical Methods in Engineering, Vol. 60, Issue 1, 2004, p. 213-232. [Search CrossRef]
  10. Longatte E., Verreman V., Souli M. Time marching for simulation of fluid-structure interaction problems. Journal of Fluids and Structures, Vol. 25, Issue 1, 2009, p. 95-111. [Search CrossRef]
  11. Hu Shiliang, Lu Chuanjing, He Yousheng Fluid-structure interaction simulation of three-dimensional elastic hydrofoil in water tunnel. Applied Mathematics and Mechanics, Vol. 37, Issue 1, 2016, p. 15-26. [Search CrossRef]
  12. Fabian D., Raul G., Srinivasan N. Arbitrary Lagrangian-Eulerian method for Navier-Stokes equations with moving boundaries. Computer Methods in Applied Mechanics and Engineering, Vol. 193, Issues 45-47, 2004, p. 4819-4836. [Search CrossRef]
  13. Liu Yun He Transient Dynamic Analysis of Fluid-Structure Coupling on the Application on Theory and Hydraulic Engineering. Xi’an Jiao Tong University, 2001, p. 76-78. [Search CrossRef]
  14. Yue Qian Bei, Dong Ri Zhi, Liu Ju Bao The research of vortex-induced vibration for the elastic pipe at different locations in a limited fluid domain. Applied Mathematics and Mechanics, Vol. 37, Issue 3, 2016, p. 277-289. [Search CrossRef]
  15. Liu Ju Bao, Luo Min, Wang Lin Methodology of fluid-structure interaction of rotary slender beam pipe. Engineering Mechanics, Vol. 28, Issue 4, 2011, p. 251-256. [Search CrossRef]