Расчет критической скорости флюида, протекающего в однослойной углеродной нанотрубке в полимерной матрице Текст научной статьи по специальности «Физика»

Лолов Д.С., Лилкова-Маркова Св.В. Расчет критической скорости флюида, протекающего в однослойной углеродной нанотрубке в полимерной матрице // Вестник Пермского национального исследовательского политехнического университета. Механика. 2019. № 4. С. 114-119. DOI: 10.15593/perm.mech/2019.4.11

Lolov D.S., Lilkova-Markova Sv.V. Determination of the critical velocity of the fluid flowing in a single-walled carbon nanotubes, embedded in a polymer matrix. PNRPUMechanics Bulletin, 2019, no. 4, pp. 114-119. DOI: 10.15593/perm.mech/2019.4.11

ВЕСТНИК ПНИПУ. МЕХАНИКА

№ 4, 2019 PNRPU MECHANICS BULLETIN

http://vestnik.pstu.ru/mechanics/about/inf/

DOI: 10.15593/perm.mech/2019.4.11 УДК 539.3

РАСЧЕТ КРИТИЧЕСКОЙ СКОРОСТИ ФЛЮИДА, ПРОТЕКАЮЩЕГО В ОДНОСЛОЙНОЙ УГЛЕРОДНОЙ НАНОТРУБКЕ В ПОЛИМЕРНОЙ МАТРИЦЕ

Д.С. Лолов, Св.В. Лилкова-Маркова

Университет архитектуры, строительства и геодезии, София, Болгария

О СТАТЬЕ

АННОТАЦИЯ

Получена: 11 марта 2019 г. Принята: 05 сентября 2019 г. Опубликована: 30 декабря 2019 г.

Ключевые слова:

нанотруба, динамическая устойчивость, критическая скорость, упругое основание, дифференциально-квадратичный метод.

Начиная с 1990-х годов прошлого века нанотрубки широко используют в нанофизике, на-нобиологии и наномеханике для транспортировки жидкости в виде наноконтейнеров - для хранения газа и для других целей. Они представляют собой полую цилиндрическую структуру диаметром от десятых до нескольких десятков нанометров и длиной от одного микрометра до нескольких сантиметров. Нанотрубки обладают высокой электропроводностью и превосходящей сталь прочностью. Рассмотрение проблемы взаимодействия «жидкость-нанотрубка» на наноуровне сопряжено со значительными трудностями и является дорогостоящим. Это основные причины расчетного исследования динамической устойчивости нанотрубок, транспортирующих жидкость, с использованием модели упругой балки Эйлера или Тимошенко.

В данной статье рассматривается динамическая устойчивость монослойной углеродной нанотрубки, встроенной в полимерную матрицу. Динамика и устойчивость исследовались на основе модели балки Эйлера и применения обобщенного дифференциально-квадратичного метода.

Исследуемая трубка встроена в полимерную матрицу и имеет шарнирное опирание. Для изучения влияния окружающей упругой среды (например, полимера) на устойчивость трубы вводится эластичное основание Пастернака. Представлено дифференциальное уравнение, описывающее поперечные колебания нанотрубки, встроенной в полимерную матрицу. Введены безразмерные параметры. Для дискретизации использована схема Чебышева-Гаусса-Лобато. Коэффициенты рассчитаны с помощью интерполяционных функций Лагранжа. Записана система однородных уравнений в матричной форме.

Исследуется влияние отношения масс (отношения массы жидкости к общей массе жидкости и трубки) на величину критической скорости жидкости (скорость, при которой система теряет устойчивость) при различной изгибной жесткости основания Пастернака. Полученные результаты представлены в графической форме. Сделаны выводы об устойчивости системы. Установлено понижение критической скорости при увеличении массового отношения (отношения массы жидкости к сумме масс жидкости и трубки на единицу длины).

©ПНИПУ

© Лолов Димитр Сенков - к.т.н., доц., e-mail: dlolov@yahoo.com, : 0000-0002-8138-0265 Лилкова-Маркова Светлана Велкова - к.т.н., проф., зав. каф., e-mail: lilkova_fhe@uacg.bg : 0000-0003-0582-8176.

Dimitar S. Lolov - CSc in Technical Sciences, Associate Professor, e-mail: dlolov@yahoo.com, Svetlana V. Lilkova-Markova - CSc in Technical Sciences, Professor, Head of Department, e-mail: lilkova_fhe@uacg.bg, : 0000-0003-0582-8176.

Эта статья доступна в соответствии с условиями лицензии Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)

This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0)

DETERMINATION OF THE CRITICAL VELOCITY OF THE FLUID FLOWING

IN A SINGLE-WALLED CARBON NANOTUBES EMBEDDED IN A POLYMER MATRIX

D.S. Lolov, Sv.V. Lilkova-Markova

University of Architecture, Civil Engineering and Geodesy, Sofia, Bulgaria

ARTICLE INFO ABSTRACT

Receded 11 March 2019 Since 90's carbonic nanotubes are broadly used in nanophysics, nanobiology and

Accepted: 05 September 2019 nanomechanics in nanofluidic devices, nanocontainers for gas storage and nanopipes conveying

Published: 30 December 2019 fluid. They have a perfect hollow cylindrical geometry and superior mechanical strength. The

--flowing fluid can be water, oil, dynamic flow of methane, ethane and ethylene molecules. The

Keywords. problem of the fluid-structure interaction could be considered in the case of nanoscale. However,

nanotube, dynamic stability, the experiments at the nanoscale are difficult and expensive. That is why the continuum elastic

critical stability, elastic foundation, models have been used to study the fluid-structure interaction. The carbon nanotubes are con-

differential quadrature method. sidered with Euler- and Timoshenko-beam models.

In this paper the dynamic stability of a single-walled carbon nanotube is investigated on the basis of the Euler-beam model and with the employment of the Generalized Differential Quadrature Method. The tube under investigation is assumed hinged at its both ends and is embedded in a polymer matrix. To study the influence of the surrounding elastic medium (for example, a polymer) on the stability of the pipe, an elastic base of Pasternak is introduced. A differential equation is presented that describes the transverse vibrations of a nanotube embedded in a polymer matrix. Dimensionless parameters are introduced. The scheme of Chebycheva-Gauss-Lobato is used for sampling. The coefficients are calculated using Lagrange interpolation functions. A system of homogeneous equations is written in the matrix form.

The obtained numerical results are for flowing fluids with different densities. In order to study the effect of the surrounding elastic medium (such as polymer) on the stability of the pipe the Pasternak elastic foundation is introduced. The critical velocities of each type of fluid are determined for different stiffnesses of this matrix. A decrease in the critical speed with the increasing mass ratio has been established.

Введение

Углеродные нанотрубки применяются в различных областях промышленности. Механические свойства этих трубок являются исключительными. Исследования по углеродным нанотрубкам проводятся с момента их открытия в 1991 году.

Ряд статей посвящен динамической устойчивости нанотрубок с текущим флюидом. В [5] показано влияние массы движущегося флюида на собственные частоты и критическую скорость одностенной углеродной нанотрубки. В [1] эти колебания изучаются в случае, когда нанотрубка лежит на упругом основании Винклера.

При исследовании собственных колебаний и устойчивости труб с протекающим флюидом часто применяется дифференциально-квадратичный метод. В [3] с использованием дифференциально-квадратичного метода решены разные задачи и показана его эффективность. В [4] приведены различные схемы дифференциально-квадратичного метода при исследовании колебаний балок для разных граничных условий.

В [2] представлена техника решения некоторых дифференциальных уравнений в частных производных. В [9] анализируется влияние плотности транспортируемой жидкости, длины трубы и размеров ее поперечного сечения на критическую скорость потока жидкости

©PNRPU

в двухслойной нанотрубке. В [10] исследуется динамика труб с флюидом на основаниях Винклера и Пастернака. Труды с [11] по [30] посвящены исследованиям колебаний различных нанотрубок с протекающим флюидом.

В [6] исследована консольная труба с протекающим флюидом. Применен обобщенный дифференциально-квадратичный метод и проведены параметрические исследования устойчивости трубы. В настоящей работе дифференциально-квадратичный метод применяется для определения критической скорости протекающего флюида в однослойной углеродной нанотрубке, лежащей на упругом основании Пастернака. Рассмотрены различные параметры этого упругого основания.

1. Методология

Однослойная углеродная нанотрубка исследуется на динамическую устойчивость с использованием континуальной модели балки Эйлера. Эта модель (рассмотренная в [8]) чрезвычайно проста и надежна в исследованиях такого типа. Эксперименты в области устойчивости углеродных нанотрубок демонстрируют эффективность модели Эйлера [7].

Дифференциальное уравнение, описывающее поперечные колебания нанотрубки с длиной Ь, изгибной

Lolov D.S., Lilkova-Markova SvV. /PNRPU Mechanics Bulletin 4 (2019) 114-119

жесткостью Е1, по которой течет флюид со скоростью V, имеет следующий вид:

д4м , 2 , \ д2м

Е1—+КУ - к +

+ 2т у

дх4 д2 м дх дг

■ + (т/ + тР )

д2 м

Иё

+ кы = 0,

(1)

где т/ и т - масса флюида на единицу длины трубы и масса трубы на единицу длины соответственно; х -осевая координата вдоль оси трубы; г - время; м (х, г) -поперечные перемещения оси трубы; к и к являются параметрами упругого основания Пастернака, моделирующего полимерную матрицу.

Для нанотрубки с шарнирными опорами на обоих концах применяются следующие граничные условия:

ы (0,г) = ы (Ь,г) = 0;

"д2 м ( х, г )" "д2 м (х, г )"

дх2 х= 0 дх2

= 0.

(2)

Решение дифференциального уравнения (1) ищется следующим образом:

м (х, г) = Г (х) еш г,

(3)

где ш - комплексная круговая частота системы.

После подставления выражения (3) в (1) и преобразования получается дифференциальное уравнение

^и4 Ж , 2 , г

Е1 ——— + (т}У - К ) , +

ё х

иЖ

ё х

и ГУ / \ О

+ 2т,Кш-+ (т, + т )ш2Г + к Г = 0. (4)

и х х '

Для удобства введены следующие безразмерные параметры:

е х Г

% = —; т =—; ь ь

\т, \т, + тР

■■ У^А—; П = ш ¿Л—--;

Е1

Р = -

т.

-; т = -

Е1

т + тр

1т+

_ к Ь - к Ь

Х„ = ^ ; £ = КЬ. (5)

Е1

Е1

Уравнение (4), записанное в безразмерной форме, принимает следующий вид:

и4 т

4 +(«2-кг) 0 + 2ф и П ^ + (^2 + к„ )т = 0. (6)

Уравнение (6) можно преобразовать в систему алгебраических уравнений дифференциально-квадратичным методом, использованным в [4], [6]. Согласно этому методу значения производных функции т(%) представляются в виде суммы умноженных с коэффи-

и т

циентами значений функции в каждой точке оси

трубы после дискретизации.

ип т(%)

и %п

=ЕР(Г)тв), ■=1,

(7)

1=1

где т - количество точек, которые разделяют нанот-рубку на участки; Р(п - коэффициенты для п-й производной в г-й точке оси трубы.

В настоящей статье ось трубы дискретизируется с использованием схемы Чебышева-Гаусса-Лобато, согласно которой

1

1 I ■ -1 1 - -л

т-1

, = 1,...,т .

(8)

Коэффициенты рассчитываются с помощью интерполяционных функций Лагранжа [6]:

ь(1)(%,)

Рм=-

Р" -О ь(1)Ы'

,,1 = 1,...,т,, Ф 1,

р(1 )=-ЕрС), ■, 1 = 1,..., т,

1=и *>

Р(; )= п

(п-1Ь(1)__\Л

Р(П-1)Р(1)-

(п-1) Л

1 У

(9) (10)

(11)

1 = 1,...,т;, Ф ]; п = 2,..,(т -1),

рг=- Е рг, ,, 1 = 1,...,т; п = 2,...,(т-1). (12)

]=1,1 ф,

В формуле (9) первая производная интерполяционных полиномов Лагранжа в каждой точке с ординатой вычисляется по формуле

т

Ь(1)(%к)=П(%к), к = 1,...,т . (13)

I=и ф к

Показано, что применение интерполяционных полиномов Лагранжа при дискретизации в схемах Чебы-шева-Гаусса-Лобато обеспечивает конвергенцию решаемой задачи. Увеличение количества точек приводит к уменьшению ошибки в решении [6].

Применением формул (7), (8), (9), (10), (11), (12) и (13) уравнение (6) представляется следующей системой алгебраических уравнений в точках 1= 3,4,...,(т-2) оси трубы:

т т

Е р(4Ч% 1 )+(и2 - кГ р(Т)л(% ])+

1=1 1=1

т _

р(1)т (% 1)+(^2 + К ) т (%,) = 0. (14)

1=1

Граничные условия (2) записываются следующим образом:

т т

) = т(%т) = 0; ЕЕр(1)т(%1 ) = Ер!т'т(%1 ) = 0. (15)

1

1

т

т

г

+

Уравнение (14), записанное в матричной форме, дает следующую форму:

+ (и2 -кг)Б(2) + 20и%/рБ(1))8 +

+ (о2 + кг) /8, = 0, (16)

В (16)

Б (4)=р(4); Б (2)=р(2); Б (1)=р(1),

V 1V ' V 1V ' V 1V '

' = 3,4,...,(т-2), V = 1,...,т являются матрицами размерностью (т - 4) х т. Векторы

8={^),...,^т)}Г , (18)

8, ={^3 ) ,...,^-2 )}Г . (19)

Четыре граничных условия записываются в матричной форме:

Къ 8 = 0, (20)

где матрица К имеет размерность (4 хт) и следующий вид:

Кь =

1 0 0 .. 0

Р1? Р? .. Р(2) .. Р1т

0 0 0. .. 1

1(2) 'т 1 Р(2) Рт 2 Р(2). Рт 3 .. Р(2) тт

(21)

Уравнения (16) и (20) могут быть записаны в системе т алгебраических уравнений с т неизвестными:

К

Б4>+'-2

(и2 - к ) Б(2)+ 20иЛ/рБ'

8 +

(о2 + К) / 8,

= 0.

(22)

После преобразований уравнение (22) принимает следующий вид:

Кьь Кь, Б(4) Б(4)

Здесь

-(и2 -кг)

0 0

Б(2) Б(2)

-20и

00

Б(1) Б(1)

\

8Ь /о — \ 0 0 8Ъ 0

+ (о2 + к ) =

8 ь V "/ 0 / 8 ь 0

(23)

Кьь =

10 0 0

в(2) в(2) в(2) в(2

Р11 Р12 Р1 (т-1) Р1,

00

0

Р(2) Р(2) Р(2) Р

рт 1 рт 2 рт (т-1) Р

(24)

Кьа =

0

Р12з' 0

Рт 3

0

р124} 0

(2) о(2) Рт 4

... 0

Р(2) ... Р1 (т-2)

... 0

... Р(2)

т (т-2)

(25)

Б(4) =

Б(4) =

р341) р342 Р(4) Р3 (т-1) Р(4) р3 т

Р(4? р442 Р(4) Р4 (т-1) Р(4) р4 т

(4) (т-3) 1 Р(4) Р(т-3) 2 Р(4) Р(т-3) (т-1) р(4) Р(т-3)

(4) (т-2) 1 р(4) Р(т-2) 2 р(4) Р(т-2) (т-1) р(4) Р(т-2)

Р33

1Ч(4)

Р344 Р444

р(4)

3 (т-2)

(4)

4 (т-2)

Р(4)

Р(т-2) 3 Р(т-2) 4

Р(,4)

Р(4)

(т-2) (т-2) Т

(27)

8ь={л(^0Л&), (28)

8, = { ) ... Л (¡^-3)) Л (V2)) }Т . (29)

Учитывается, что

8ь =-К-1Кьа 8а.

(30)

В результате получается система с (т - 4) уравнениями и (т - 4) неизвестными:

г м БУ Б1' КььКь<з )

8а +

(Б$?-К- 8, +(«2 -к^( +20«^« -Б0ЫК-Кьа) 8, + (о2 + К)/8, = 0. (31)

Чтобы уравнение (31) имело ненулевое решение, определитель перед вектором 8ё должен быть нулевым. Этот определитель является функцией неизвестной критической скорости и и неизвестной круговой частоты 0 = Яе 0+'1т0. В случае когда Яе 0< 0, система устойчива. Если Яе 0> 0, система теряет свою стабильность. Если Яе 0 = 0, система находится в состоянии безразличного равновесия.

Следующая вычислительная процедура используется для определения критической скорости: задается конкретное значение и , и из условия, что определитель перед вектором 8Л в уравнении (31) равен нулю, определяется круговая частота О, на основе которой система оценивается как стабильная.

2. Численные исследования

Для иллюстрации предложенного метода исследована устойчивость однослойной углеродной нанотрубки с внешним радиусом поперечного сечения Я = 55 нм и

толщиной И = 10 нм. Длина трубки Ь = 2000 нм. Трубка имеет шарнирное опирание. Исследованы случаи жидкостей с плотностями в диапазоне р = [700; 1600] кг/м3. Определена критическая скорость флюида Усг для параметров упругого основания Пастернака к^ = 100 кПа;

43

0

+

X

тт

кг = 100 Нм/м и = 100 кПа; кг = 200 Нм/м. Полученные результаты представлены на рисунке. Параметры ^ и кг соответствуют механическим характеристикам полимерной матрицы, в которую вставлена иследованная нанотрубка.

Исследования показали, что при т > 14 представленные результаты для разных значений т довольно похожи. Наблюдаемые отклонения < 4 %.

Выводы

Из графика на рисунке видно, что снижение жесткости кг упругого основания оказывает стабилизирующее воздействие на систему. В пределах исследуемого интервала Р более высокие критические скорости флюида получены для упругого основания с меньшей жесткостью кг. В обоих исследованных упругих осно-

ваниях видно снижение критической скорости при увеличении массового отношения Р .

Рис. Зависимость критической скорости флюида Vcr от массового отношения Р

Fig. Dependence of the critical fluid velocity Vcr and the mass ratio Р

Библиографический список

1. Belhadi A., Boukhalfa A., Belalia S.A. Free vibration modeling of single-walled carbon nanotubes using the differential quadrature method // Mathematical Modeling of Engineering Problems. - 2017. - Vol. 4. - No. 1. - С. 33-37.

2. Bellman R., Kashef B.G., Casti J. Differential quadrature: a technique for the rapid solution of nonlinear partial differential equations // Journal of Computational Physics. - 1972 - Vol. 10. -Р. 40-52.

3. Bert C.W., Malik M. Differential quadrature method in computational mechanics: a review // Applied Mechanics Reviews -1996 - Vol. 49. - Р. 1-27.

4. Ng C.H.W., Zhao Y. Xiang, G.W. We. On the accuracy and stability of a variety of differential quadrature formulations for the vibration analysis of beams // International Journal of Engineering and Applied Sciences (IJEAS) - 2009. - Vol. 1. -Iss. 4. - Р. 1-25.

5. Free vibration analysis of fluid conveying single-walled carbon nanotubes / C.D. Reddy, C. Lu, S. Rajendran, K.M. Liew // Applied Physics Letters. - 2007. - Vol. 90 - P. 133.

6. Critical flow speeds of pipes conveying fluid using Generalized Quadrature Method / F. Tornabene, A. Marzani, E. Viola, I. Elisha-koff // Adv. Theor. Appl. Mech. - 2010. - No. 3. - С. 121-138.

7. Yoon J., Ru C., Mioduchowski A. Vibration and instability of carbon nanotubes conveying fluid // Composites Science and Technology - 2005 - No. 65. - С. 1326-1336.

8. Лилкова-Маркова Св.В., Лолов Д.С. Устойчивост на тръби с протичащ флуид. - София: АВС Техника, 2016. - 124 с.

9. Lolov D.S., Lilkova-Markova Sv.V. Dynamic stability of double-walled carbon nanotubes // Journal of the Serbian Society for Computational Mechanics. - 2018. - Vol. 12. - No. 1 - С. 1-8.

10. Lilkova-Markova Sv.V., Lolov D.S. Free and forced vibrations of a pipe on Winkler and Pasternak elastic foundations // Beams and Frames on Elastic Foundation 3. - Ostrava, 2010. -P. 337-343.

11. Hanasaki I., Nakatani A. Water flow through carbon nano-tube junctions as molecular convergent nozzles // Nanotechnology. -2006. - Vol. 17 (11) - С. 2794-2804.

12. Anomalously immobilized water: a new water phase induced by confinement in nanotubes / R.J. Mashl, S. Joseph, N.R. Aluru, E. Jakobsson // Nano Lett. - 2003 - Vol. 3 - P. 589-592.

13. Sokhan V.P., Nicholson D., Quirke N. Fluid flow in nanopores: an examination of hydrodynamic boundary conditions // J. Chem. Phys. - 2001 - Vol. 115 (8). - P. 3878-3887.

14. Molecular dynamics simulation of contact angle of water droplet in carbon nanotubes / T. Werder, J.H. Walther, R.L. Jaffe, T. Halicioglu, F. Noca, P. Koumoutsakos // Nano Lett. - 2001. -Vol. 1. - C. 697-702.

15. Rapid transport of gases in CNTs / A. Skoulidas, D.M. Ackerman, K.J. Johnson, D.S. Sholl // Phys. Rev. Lett. -2002. - Vol. 89 (6). - P. 185-901.

16. Dynamics of fluid flow inside CNTs / R.E. Tuzun, D.W. Noid, B.G. Sumpter, R.C. Merkle // Nanotechnology. -1996. - Vol. 7. - P. 241-246.

17. Mahan G.D. Oscillations of a thin hollow cylinder: carbon nanotubes // Phys. Rev. - 2002. - B 65. - P. 235.

18. Popov V.N., Doren V.E. Van, Balkanski M. Elastic properties of single-walled carbon nanotubes // Phys. Rev. - 2000. -B 61. - P. 3078-3084.

19. Gogotsi Y. In situ multiphase fluid experiments in hydrothermal CNTs // Appl. Phys. Lett. - 2001. - Vol. 79. - P. 1021-1023.

20. Pozhar L.A., Structure and dynamics of nanofluids: theory and simulations to calculate viscosity // Phys. Rev. - 2000. -E 61 (2). - P. 1432.

21. Yoon J., Ru C.Q., Mioduchowski A. Sound wave propagation in multiwall carbon nanotubes // J. Appl. Phys. - 2003. -Vol. 93 (8). - P. 4801.

22. Supple S., Quirke N. Rapid imbibition of fluids in CNTs // Phys. Rev. Lett. - 2003. - Vol. 90. - P. 214.

23. Mao Z., Sinnott S.B. A computational study of molecular diffusion and dynamics flow through CNTs // J. Chem. Phys. -2000. - B 104. - C. 4618-4624.

24. Yoon J., Ru C.Q., Mioduchowski A. Noncoaxial resonance of an isolated multi-walled carbon nanotubes // Phys. Rev. -2002. - B 66.

25. Yoon J., Ru C.Q., Mioduchowski A. Flow-induced flutter instability of cantilever CNTs // Int. J. Solids Struct. - 2006. -Vol. 43. - P. 3337-3349.

26. Yoon J., Ru C.Q., Mioduchowski A. Vibration and instability of CNTs conveying fluid // Compos. Sci. Technol. - 2005. -Vol. 65 - P. 1326-1336.

27. Wang X.Y., Wang X., Sheng G.G. The coupling vibration of fluid-filled carbon nanotubes // J. Phys. D: Appl. Phys. -2007. - Vol. 40 - P. 2563-2572.

References

1. Belhadi A., Boukhalfa A., Belalia S.A. Free vibration modeling of single-walled carbon nanotubes using the differential quadrature method. Mathematical modeling of engineering problems, 2017, vol. 4, no.1, pp. 33-37.

2. Bellman R., Kashef B.G., Casti J. Differential quadrature: a technique for the rapid solution of nonlinear partial differential equations. Journal of Computational Physics, 1972, vol. 10, pp. 40-52.

3. Bert C.W., Malik M. Differential quadrature method in computational mechanics: a review. Applied Mechanics Reviews. 1996, vol. 49, pp. 1-27.

4. Ng C.H.W., Zhao, Y. Xiang, G. W. We. On the accuracy and stability of a variety of differential quadrature formulations for the vibration analysis of beams. International journal of engineering and applied scinces (IJEAS), 2009, vol.1, is. 4, pp.1-25.

5. Reddy, C.D., Lu C., S. Rajendran S., Liew , K.M. Free vibration analysis of fluid conveying single-walled carbon nano-tubes. Applied Physics Letters, 2007, no. 90, pp.133.

6. Tornabene, F., Marzani, A., Viola, E., Elishakoff I. Critical flow speeds of pipes conveying fluid using Generalized Quadrature Method. Adv. Theor. Appl. Mech., 2010, no.3, pp. 121-138.

7. Yoon, J., Ru, C., Mioduchowski, A. Vibration and instability of carbon nanotubes conveying fluid. Composites Science and Technology, 2005, no.65, pp. 1326-1336.

8. Lilkova-Markova Sv. V., Lolov D.S. Stability of pipes conveying fluid. Sofia, ABC Technika, 2016, 124 p.

9. Lolov D.S., Lilkova-Markova Sv. V. Dynamic stability of double-walled carbon nanotubes. Journal of the Serbian Society for Computational Mechanics. 2018, vol.12, no 1, pp.1-8.

10. Lilkova-Markova Sv. V., Lolov D.S. Free and forced vibrations of a pipe on Winkler and Pasternak elastic foundations. Beams and frames on elastic foundation 3, Ostrava, 2010, pp. 337-343.

11. Hanasaki I., Nakatani A., Water flow through carbon nanotube junctions as molecular convergent nozzles. Nanotech-nology, 2006, vol. 17, no. 11 pp. 2794-2804.

12. Mashl R.J., Joseph S., Aluru N.R., Jakobsson E., Anomalously immobilized water: a new water phase induced by confinement in nanotubes. Nano Lett, 2003, no. 3, pp. 589-592.

13. Sokhan V.P., Nicholson D., Quirke N., Fluid flow in nanopores: an examination of hydrodynamic boundary conditions. J. Chem. Phys., 2001, vol. 115, no. 8, pp. 3878-3887.

14. Werder T., Walther J.H, Jaffe R.L., Halicioglu T., Noca F., Koumoutsakos P., Molecular dynamics simulation of contact angle of water droplet in carbon nanotubes. Nano Lett., 2001, no. 1, pp. 697-702.

28. Dong K., Wang X. Wave propagation in carbon nanotubes under shear deformation // Nanotechnology. - 2006. - Vol. 17. -P. 2773-2782.

29. Dong K., Wang X., Sheng G.G. Wave dispersion characteristics in fluid-filled carbon nanotubes embedded in an elastic medium // Model. Simul. Mater. Sci. Eng. - 2007. - Vol. 15 - P. 427-439.

30. Anomalous potential barrier of double-wall carbon nanotube / R. Saito, R. Matsuo, T. Kimura G. Dresselhaus, M.S. Dresselhaus // Chem. Phys. Lett. - 2001. - Vol. 348. - P. 187-193.

15. Skoulidas A., Ackerman D.M., Johnson K.J., Sholl D.S., Rapid transport of gases in CNTs. Phys. Rev. Lett., 2002, vol. 89, no. 6, pp. 185-901.

16. Tuzun R.E., Noid D.W., Sumpter B.G., Merkle R.C. Dynamics of fluid flow inside CNTs Nanotechnology, 1996, no. 7, pp. 241-246.

17. Mahan G.D., Oscillations of a thin hollow cylinder: carbon nanotubes. Phys. Rev., 2002, B 65, pp. 235.

18. Popov V.N., Van Doren V.E., Balkanski M., Elastic properties of single-walled carbon nanotubes. Phys. Rev., 2000, B 61, pp. 3078-3084.

19. Gogotsi Y., In situ multiphase fluid experiments in hydrothermal CNTs. Appl. Phys. Lett., 2001, vol. 79, pp. 1021-1023.

20. Pozhar L.A., Structure and dynamics of nanofluids: theory and simulations to calculate viscosity. Phys. Rev., 2000, E 61 (2), pp. 1432.

21. Yoon J., Ru C.Q., Mioduchowski A., Sound wave propagation in multiwall carbon nanotubes. J. Appl. Phys., 2003, vol. 93, no. 8, pp. 4801.

22. Supple S., Quirke N., Rapid imbibition of fluids in CNTs. Phys. Rev. Lett., 2003, vol. 90, pp. 214.

23. Mao Z., Sinnott S.B., A computational study of molecular diffusion and dynamics flow through CNTs. J. Chem. Phys., 2000, B 104, pp. 4618-4624.

24. Yoon J., Ru C.Q., Mioduchowski A., Noncoaxial resonance of an isolated multi-walled carbon nanotubes. Phys. Rev., 2002, B 66.

25. Yoon J., Ru C.Q., Mioduchowski A, Flow-induced flutter instability of cantilever CNTs, Int. J. Solids Struct., 20066 no. 43, pp. 3337-3349.

26. Yoon J., Ru C.Q., Mioduchowski A, Vibration and instability of CNTs conveying fluid. Compos. Sci. Technol, 2005, no. 65, pp. 1326-1336.

27. Wang X.Y., Wang X., Sheng G.G., The coupling vibration of fluid-filled carbon nanotubes. J. Phys. D: Appl. Phys, 2007, no. 40, pp.2563-2572.

28. Dong K., Wang X., Wave propagation in carbon nanotubes under shear deformation. Nanotechnology, 2006, no. 17, pp. 2773-2782.

29. Dong K., Wang X., Sheng G.G., Wave dispersion characteristics in fluid-filled carbon nanotubes embedded in an elastic medium. Model. Simul. Mater. Sci. Eng., 2007, 15, pp. 427-439.

30. Saito R., Matsuo R., Kimura T., Dresselhaus G., Dresselhaus M.S., Anomalous potential barrier of double-wall carbon nanotube. Chem. Phys. Lett., 2001, no. 348, pp.187-193.