Published: 13 June 2022

Dynamical stability analysis for a micro shell subjected to swirling annular flow including the coupling effect of small size and fluid rotation

Cheng Chen1
Xinlei Chen2
Wenbo Ning3
Junzhang Wu4
Yuehong Liu5
1, 2, 3, 4Jiangsu Key Laboratory of Advanced Manufacturing Technology, HuaiYin Institute of Technology, Huai’an, China
5Huaian Weili Pump Technology Co., Ltd, Huai’an, China
Corresponding Author:
Cheng Chen
Views 27
Reads 7
Downloads 396

Abstract

This paper is concerned with stability analysis for an outer thin-walled shell conveying a swirling flowing fluid in the annular space between the inner shell and outer shell with micro-dimensions. It describes the development of theoretical model in the framework of the modified coupled stress theory along with the Donnell shell theory. The theoretical model uses the Hamilton’s principle to derive the governing motion equations of the shell. The fluid is assumed to be incompressible and irrotational and the hydrodynamic pressure is considered with the version of the potential flow theory. The zero-level contour method is utilized to display the complex dynamical behaviors of micro shell subjected to swirling annular flow. The effects of the material intrinsic parameter, the fluid rotation on stabilities of the fluid microscale shell system are discussed fully. The results elucidate that the coupling effect of the fluid rotation and the material intrinsic parameter on such system.

1. Introduction

As a fundamental structural element, macro-scale shells/tubes containing flowing fluid are encountered in many engineering applications. A remarkable number of scientists and engineers have made great efforts on different aspects of the topic. Especially, dynamical behaviors of coaxial cylindrical shells conveying axial flow have been considered [1, 2], and later the further work of shells conveying helical annular flow was performed [3].

Owing to low power, super mechanical and electrical properties, micro/nano structures such as shells/tubes have become more important in the past fifty years. The classical continuum mechanics cannot predict size-dependent responses of the microstructures. To explain the size dependent mechanical behaviors of such microstructures, some outstanding researchers have developed the nonlocal elasticity theory/couple stress theory /strain gradient theory, in which the modified couple stress theory [4] has been used widely because of its symmetric form and the involvement of only one material length scale parameter. Based on the Differential Quadrature Method (DQM), the dynamical responses of fluid-conveying microtubes were investigated by Wang [5] from the modified coupled stress theory. The threshold values of instability are higher because the modified couple stress theory considered the size effects. The nonlinear studies of fluid-conveying micro pipes was studied by the literature [6]. It was found that the material intrinsic parameter causes micro pipes more stable. And complex vibrations could be observed when the micropipe boundary meets certain conditions. Based on the nonlocal strain gradient theory, Farajpour et al. [7] extended this work. They also found that the nonlinear dynamical behaviors of the nanotubes including chaos and period-1, 3 motions could occur. Considering the fluid viscosity and density, the literature [8] performed the investigations of the nonlinear vibration and instability of fluid-conveying nanotubes. Further works on the wave propagation characteristics of fluid-filled nanotubes were performed [9]. Furthermore, the dynamical responses of micro shells/tubes is often subjected to thermal loads, electromagnetic loads. Arani et al. [10] studied the effects of the magical field on the stability of fluid-conveying micro pipes. For small-scale structures including micro/nano-shells, the reports on the mechanics of such structures were presented by Farajpour et al. [11] and Ghayesh and Farajpour [12]. For the micro/nano-system, the material intrinsic parameter must be considered. The material intrinsic parameter is not unchangeable, which was confirmed by Khorshidi [13].

The present paper performs numerical calculation to the results by means of the software package MATLAB. It will be shown that the fluid rotation effect, size effect on the dynamical behaviors of the shell. These studies will provide useful information in design and utilization of microfluidic structures.

2. Theoretical model

2.1. The definition of the problem

Fig. 1 displays the two coaxially micro shells under consideration with same materials properties by elastic modulus E, Poission’ ratio υ, and density ρs. The length, radius and wall thickness of the inner shell are denoted by L, rihi. And corresponding quantities of the outer shell are L, ro and ho. In a cylindrical coordinate system O;r,θ,x, the origin O is set at the middle surface of the inner shell. x, θ, and r denote in the axial, circumferential, and radial directions, respectively.

The density of an incompressible helical fluid is ρf. The axial mean flow velocity of the annular flow is Uo and the angular mean flow velocity is Ωo.

We raised three assumptions in this paper: (1) the shell motion is small; (2) the fluid is irrotational and isothermal; (3) the gap between the two shells is much smaller than the length L.

Fig. 1Sketch of the coaxial micro shells under swirling annular flow

Sketch of the coaxial micro shells under swirling annular flow

2.2. Hydrodynamic pressure

Due to the small shell motion, the inviscid fluid forces induced by shell motions are obtained from the potential flow theory. From the impermeability boundary conditions on the shell wall surfaces, the inviscid fluid forces p (the perturbation pressure) on the inner shell or the outer shell are given [14], because the fluid cannot flow on the interface between the shell and the fluid:

1
pi=-ρfΓt+Vorθ+Uox2InλriKn'λro-In'λroKnλriwi,
2
po=-ρfΓt+Vorθ+Uox2In'λriKnλro-InλroKn'λriwo,

in which Γ=In'(λri)Kn'(λro)-In'(λro)Kn'(λri); the mean tangential velocity Vo=r×Ωo, where In and Kn are the nth-order the first and second modified Bessel functions; the subscripts i and o are the inner shell and the outer shell, respectively; t is the time; the prime indicates the differentiation; w is the r-direction displacement; λ pertains to the order of axial vibration mode.

2.3. The governing equations of motion

Compared with the traditional couple stress theory, the modified couple stress theory proposed by Yang et al. [4] has been using widely due to its compaction and symmetry. Based on the modified couple stress theory and the Hamilton principle, the shell motions are given:

3
Nxx+1rNθxθ+122Yxzrxθ+122Yθzr2θ2=ρs2ut2,
4
Nxθx+1rNθθ-122Yxzx2-122Yθzrxθ=ρs2vt2,
5
2Mxx2+1r2Mxθxθ+1r2Mθxxθ-Nθr-1r2Yxxθ+1r2Yθxθ-122Yxθr2θ2+122Yxθx2=ρs2wt2+pi/o,

where Nij is the internal forces; Mij is the moment resultants; γij is the couple forces; Tij is the couple moments; i, j=x, θ, r.

The shell motion equations are given:

6
Ζui/ovi/owi/o=00pi/o,

where u, v and w are the x-, θ-, and r- directions displaces of a shell, respectively; Ζ is the linear matrix operator.

There exists non-zero solutions from the Eq. (6), we arrive at:

7
Ζ=0.

We use the zero-level contour method to determine the dispersion relation of ω--U-oV-o(V-o/U-o is given) by the software package MATLAB. From the feature points of the frequency-fluid velocity curve, we can find the instability velocities and correspond modes.

2.4. Validation of the present method

After using the zero-level contour method, the calculated dimensionless results are tabulated in Table 1, and are compared with results of Zhou and Wang [15].

From the comparison study, the maximum error is 6.0 %. It can be seen that the obtained critical flow velocities in present paper are good agreement with those in published literature.

Table 1Comparisons of critical axial velocities instability by divergence for a simply supported micro shell

Modified couple stress theory
Class theory
Uo (m/s)
84.6 [15]
85.8 (present)
14.8 [15]
13.9 (present)
Mode (m, n)
(1, 0)
(3, 1)
(1, 2)
(1, 1)

3. Results and discussion

In this paper, the case of elastic outer and perfectly stiff inner shells is considered. We use the simply supported boundary conditions to analyze the dynamical behavior of micro shell. The boundary conditions are given in the form of displaces.

We adopt the following geometry parameters:

L=100 μm, ro=10 μm, k=ro-ri/ri, ri=9.0901 μm (k=0.1),
h=hi=ho=0.1 μm.

The material properties in the calculations are chosen as: E= 1.44×1011 Pa, υ= 0.38, ρs= 1000 kgm-3, ρf= 1000 kgm-3, l=h= 0.1 μm.

3.1. Effects of small size

For the epoxy micro shell in present paper, the definition of the characteristic length is given. l=bh/3(1-υ), where l is the material intrinsic parameter; bh is a higher order bending parameter.

The material intrinsic parameter l indeed affects the results. To illustrate the effect of the material intrinsic parameter on the critical velocity, we obtain different l after defining and calculating bh, as listed in Table 2. Using the same data, the obtained critical velocities by divergence are summarized in Table 3.

Table 2Different values of the material intrinsic parameter l with different bh

bh
0
0.14
0.50
1.00
l
0
0.10
0.37
0.74
l'=l/ri
0
0.0161
0.0403
0.0807

Table 3Comparisons of dimensionless critical axial velocities by divergence with different l

l'
0
0.0161
0.0403
0.0807
U-0
0.00498
0.0060
0.0151
0.0158
modes (m, n)
(1, 1)
(1, 2)
(1, 2)
(1, 2)

It is clear that, the material intrinsic parameter l exerts a very significant influence on the critical axial velocity, which goes up with the increase of l. When keeping the length-radius ratio, thickness-radius ratio and gap ratio constant, the relative radius of micro shell ro'=ro/roi=ro/10 changes from 1 to 50.

Fig. 2 shows the variation of the critical velocities with different relative radii. We found that the material intrinsic parameter effect gradually diminishes when the radius of micro-shell becomes larger enough. To be concluded, the size effect cannot be considered for macro structures.

Fig. 2Variation of dimensionless critical axial flow velocities by divergence with different relative radii

Variation of dimensionless critical axial flow velocities by divergence with different relative radii

Fig. 3Effects of fluid rotation on the instability of the micro shell

Effects of fluid rotation  on the instability of the micro shell

3.2. Effects of fluid rotation

Fig. 3 depicts the effect of the fluid rotation on the instability behavior of the micro shell. It can be seen that with the increase of the rotating speed of fluid, the critical axial velocity by divergence evidently reduces to zero, which indicating the compound form of divergence and flutter changing into the flutter. We obtain the dimensionless critical angular velocity (0.00086) at the transition point by solving Eq. (7). It should be pointed out that the influence of the angular flow velocity for micro shells is similar to that for macro-shells.

The dimensionless total flow velocity is defined as Vt=U-o2+V-o2.

The effect of the fluid rotation on the instability of shell base on the modified coupled stress theory and the classical continuum mechanics theory is illustrated in Fig. 4. From the three curves in the figure, it is seen that the critical total fluid velocities decrease when the ratio of the tangential velocity to the axial velocity increases, and thus the fluid ration lowers the stability of the shell. Furthermore, there appears to be a nonuniform decrease of the total critical flow velocity with the increase of the angular flow velocity. It is also seen that the two curves (for l= 0, or 0.1) present more obvious downward trend while the third curve (for l= 17.6) is declined slowly with increasing of the tangential velocity. To be concluded that, the macro-shell can be affected easily the influences of fluid rotation and the micro shell considering the size effect might be more stable when the material intrinsic parameter l becomes larger.

Fig. 4Effects of the ratio of the tangential velocity to the axial velocity on the critical total flow velocity

Effects of the ratio of the tangential velocity to the axial velocity on the critical total flow velocity

Fig. 5Coupling effect of the fluid rotation and size on the system instability

Coupling effect of the fluid rotation  and size on the system instability

3.3. Coupling effects of fluid rotation and small size

We define the dimensionless physical parameter Rlv=l/Vt, which represents the coupling effects of the fluid rotation and small size.

Fig. 5 displays the coupling effects of the fluid rotation and the material intrinsic parameter under different dimensionless parameters Rlv. It is seen that the effects of the fluid rotation or the material intrinsic parameter on the stability of the system are not the same. It is also seen that for a big enough Rlv, the material intrinsic parameter changes have greater influence on the system instability than the fluid rotation. The coupling effects might be considered in industrial applications.

4. Conclusions

This paper is concerned with the size-dependent stability analysis of the outer cylindrical shell under the annular flow with two flow velocity components based on the modified couple stress theory by the aid of the Hamilton’s principle. The obtained results indicate that, the fluid rotation lowers the stability of fluid-shell system. Furthermore, when the fluid speed value of rotation increases a certain value, the form of instability could be changed. The obtained results indicate that, the material intrinsic parameter effect results in the fluid-micro shell system more stable. The coupling effect of the small size and fluid rotation on the system stability are evaluated by introducing the new physical coefficient.

References

  • M. P. Païdoussis, S. P. Chan, and A. K. Misra, “Dynamics and stability of coaxial cylindrical shells containing flowing fluid,” Journal of Sound and Vibration, Vol. 97, No. 2, pp. 201–235, 1984, https://doi.org/10.1016/0022-460x(84)90319-5
  • M. P. Paidoussis, A. K. Misra, and S. P. Chan, “Dynamics and stability of coaxial cylindrical shells conveying viscous fluid,” Journal of Applied Mechanics, Vol. 52, No. 2, pp. 389–396, Jun. 1985, https://doi.org/10.1115/1.3169059
  • T. L. C. Chen and C. W. Bert, “Dynamic Stability of Isotropic or Composite-Material Cylindrical Shells Containing Swirling Fluid Flow,” Journal of Applied Mechanics, Vol. 44, No. 1, pp. 112–116, Mar. 1977, https://doi.org/10.1115/1.3423973
  • F. Yang, A. C. M. Chong, D. C. C. Lam, and P. Tong, “Couple stress based strain gradient theory for elasticity,” International Journal of Solids and Structures, Vol. 39, No. 10, pp. 2731–2743, May 2002, https://doi.org/10.1016/s0020-7683(02)00152-x
  • L. Wang, “Size-dependent vibration characteristics of fluid-conveying microtubes,” Journal of Fluids and Structures, Vol. 26, No. 4, pp. 675–684, May 2010, https://doi.org/10.1016/j.jfluidstructs.2010.02.005
  • K. Hu, Y. K. Wang, H. L. Dai, L. Wang, and Q. Qian, “Nonlinear and chaotic vibrations of cantilevered micropipes conveying fluid based on modified couple stress theory,” International Journal of Engineering Science, Vol. 105, pp. 93–107, Aug. 2016, https://doi.org/10.1016/j.ijengsci.2016.04.014
  • A. Farajpour, H. Farokhi, and M. H. Ghayesh, “Chaotic motion analysis of fluid-conveying viscoelastic nanotubes,” European Journal of Mechanics – A/Solids, Vol. 74, pp. 281–296, Mar. 2019, https://doi.org/10.1016/j.euromechsol.2018.11.012
  • A. Ghorbanpour Arani, M. R. Bagheri, R. Kolahchi, and Z. Khoddami Maraghi, “Nonlinear vibration and instability of fluid-conveying DWBNNT embedded in a visco-Pasternak medium using modified couple stress theory,” Journal of Mechanical Science and Technology, Vol. 27, No. 9, pp. 2645–2658, Sep. 2013, https://doi.org/10.1007/s12206-013-0709-3
  • Y. Yang, J. Wang, and Y. Yu, “Wave propagation in fluid-filled single-walled carbon nanotube based on the nonlocal strain gradient theory,” Acta Mechanica Solida Sinica, Vol. 31, No. 4, pp. 484–492, Aug. 2018, https://doi.org/10.1007/s10338-018-0035-5
  • A. Ghorbanpour Arani, S. Amir, P. Dashti, and M. Yousefi, “Flow-induced vibration of double bonded visco-CNTs under magnetic fields considering surface effect,” Computational Materials Science, Vol. 86, pp. 144–154, Apr. 2014, https://doi.org/10.1016/j.commatsci.2014.01.047
  • A. Farajpour, M. H. Ghayesh, and H. Farokhi, “A review on the mechanics of nanostructures,” International Journal of Engineering Science, Vol. 133, pp. 231–263, Dec. 2018, https://doi.org/10.1016/j.ijengsci.2018.09.006
  • M. H. Ghayesh and A. Farajpour, “A review on the mechanics of functionally graded nanoscale and microscale structures,” International Journal of Engineering Science, Vol. 137, pp. 8–36, Apr. 2019, https://doi.org/10.1016/j.ijengsci.2018.12.001
  • M. Akbarzadeh Khorshidi, “The material length scale parameter used in couple stress theories is not a material constant,” International Journal of Engineering Science, Vol. 133, pp. 15–25, Dec. 2018, https://doi.org/10.1016/j.ijengsci.2018.08.005
  • W.-B. Ning, R. Gao, Y. Li, and S. Wang, “Size-dependent stability analysis of a functionally graded cylindrical shell subjected to swirling annular flow including the fluid viscosity,” Ocean Engineering, Vol. 224, p. 108755, Mar. 2021, https://doi.org/10.1016/j.oceaneng.2021.108755
  • X. Zhou and L. Wang, “Vibration and stability of micro-scale cylindrical shells conveying fluid based on modified couple stress theory,” Micro and Nano Letters, Vol. 7, No. 7, pp. 679–684, 2012, https://doi.org/10.1049/mnl.2012.0184

About this article

Received
14 May 2022
Accepted
31 May 2022
Published
13 June 2022
SUBJECTS
Flow induced structural vibrations
Keywords
micro shell
coupling effect
size effect
annular flow
modified coupled stress
Acknowledgements

The authors would like to gratefully acknowledge the support provided by the Huai’an natural science research plan project (joint special project) (No. HABL202103). The authors also thank the reviewers for their valuable suggestions.