Experimental and model study of swirling fluid flow in a converging channel as a simulation of blood flow in the heart and aorta

Мұқаба

Дәйексөз келтіру

Толық мәтін

Ашық рұқсат Ашық рұқсат
Рұқсат жабық Рұқсат берілді
Рұқсат жабық Тек жазылушылар үшін

Аннотация

The study of swirling flows in channels corresponding to the static approximation of flow channels of the heart and major vessels with a longitudinal-radial profile zR2 = const and a concave streamlined surface at the beginning of the longitudinal coordinate has been carried out. A comparative analysis of the flow structure in channel configurations zRN = const, where N = –1; 1; 2; 3, in the absence and presence of a concave surface was carried out. The numerical modelling was compared with the results of hydrodynamic experiments on the flow characteristics and the shape of the flow lines. The numerical model was used to determine the velocity structure, viscous friction losses and shear stresses. Numerical modelling of steady-state flows for channels without a concave surface showed that in the channel zR2 = const there is a stable vortex flow structure with the lowest viscous friction losses. The presence of a concave surface of sufficient size significantly reduces viscous friction losses and shear stresses in both steady state and pulsed modes.

Толық мәтін

Рұқсат жабық

Авторлар туралы

Ya. Zharkov

A. N. Bakulev National Medical Research Center for Cardiovascular Surgery

Email: agorodkov@bk.ru
Ресей, Moscow

Sh. Zhorzholiani

A. N. Bakulev National Medical Research Center for Cardiovascular Surgery

Email: agorodkov@bk.ru
Ресей, Moscow

А. Sergeev

A. N. Bakulev National Medical Research Center for Cardiovascular Surgery

Email: agorodkov@bk.ru
Ресей, Moscow

A. Agafonov

A. N. Bakulev National Medical Research Center for Cardiovascular Surgery

Email: agorodkov@bk.ru
Ресей, Moscow

A. Gorodkov

A. N. Bakulev National Medical Research Center for Cardiovascular Surgery

Хат алмасуға жауапты Автор.
Email: agorodkov@bk.ru
Ресей, Moscow

L. Bockeria

A. N. Bakulev National Medical Research Center for Cardiovascular Surgery

Email: agorodkov@bk.ru

Academician

Ресей, Moscow

Әдебиет тізімі

  1. Frazin L.J., et al. Confirmation and initial documentation of thoracic and abdominal aortic helical flow. An ultrasound study //ASAIO Journal (American Society for Artificial Internal Organs: 1992). 1996. Т. 42. №. 6. С. 951–956.
  2. Gorodkov A., Dobrova N.B., Kuzmina N.B., et al. Anatomical structures determining blood flow in the heart left ventricle / // Journal of Materials Science: Materials in Medicine. 1996. V. 7. № 3. P. 153–160.
  3. Жоржолиани Ш.Т., Миронов А.А., Талыгин Е.А. и др. Анализ динамической геометрической конфигурации проточного канала аорты с позиций смерчевой самоорганизации потока крови / // Бюллетень экспериментальной биологии и медицины. 2017. Т. 164. № 10. С. 519–524.
  4. Кикнадзе Г. И. и др. О структуре потока в левом желудочке сердца и аорте с применением точных решений нестационарных уравнений гидродинамики и морфометрических исследований // Докл. АН. 1996. Т. 351. №. 1. С. 119.
  5. Bockeria L.A., Gorodkov A.Y., Kiknadze G.I., Gachechiladze I.A. Application of Tornado-flow fundamental hydrodynamic theory to the study of blood flow in the heart – Further development of Tornado-like jet technology // ASME 2011 International Mechanical Engineering Congress and Exposition, IMECE 2011, Denver, CO, 11–17 ноября 2011 г. V. 2. Denver, CO, 2011. P. 287–296.
  6. Талыгин Е.А., Зазыбо Н.А., Жоржолиани Ш.Т. и др. Количественная оценка состояния внутрисердечного потока крови по динамической анатомии левого желудочка сердца на основании точных решений нестационарных уравнений гидродинамики для класса смерчеобразных потоков вязкой жидкости // Успехи физиологических наук. 2016. Т. 47, № 1. С. 48–68.
  7. Zhorzholiani S.T., Talygin E.A., Krasheninnikov S.V., et al. Elasticity Change along the Aorta is a Mechanism for Supporting the Physiological Self-organization of Tornado-like Blood Flow // Human Physiology. 2018. V. 44. № 5. P. 532–540.
  8. Г.И. Кикнадзе, Ю.К. Краснов. Эволюция смерчеобразных течений вязкой жидкости // Докл. АН СССР, 290:6 (1986). 1315–1319.
  9. Burgers J.M. A mathematical model illustrating the theory of turbulence // Advances in applied mechanics. 1948. Т. 1. С. 171–199.
  10. Кикнадзе Г.И., Талыгин Е.А., Городков А.Ю. Патент № 2691705 C1 Российская Федерация, МПК F15D 1/00. Способ отсасывания пограничного слоя сплошной среды с поверхности тела и устройство для его реализации : № 2018119493 : заявл. 28.05.2018 : опубл. 17.06.2019
  11. Пиралишвили Ш.А., Писаревский А.С. Численное моделирование вихревого эффекта в несжимаемой жидкости // Известия Российской академии наук. Механика жидкости и газа. 2013. №. 3. С. 138–147.
  12. Гайфуллин А.М., Жвик В.В. Взаимодействие двух противоположно закрученных затопленных струй // Известия РАН. Механика жидкости и газа. 2019. № 3. С. 48–57.
  13. Платонов Д.В., Минаков А.В., Дектерев А.А., Сентябов А. В. Численное моделирование пространственных течений с закруткой потока // Компьютерные исследования и моделирование. 2013. Т. 5. № 4. С. 635–648.
  14. Высотина В.Г. Структура закрученного течения в осесимметричных каналах // Матем. моделирование. 2003. Т. 15. № 1. С. 69–77.
  15. Shipkowitz T., et al. Numerical study on the effect of steady axial flow development in the human aorta on local shear stresses in abdominal aortic branches //Journal of biomechanics. 1998. Т. 31. №. 11. С. 995–1007.
  16. Wilcox D.C. Formulation of the kw turbulence model revisited // AIAA journal. 2008. Т. 46. №. 11. С. 2823–2838.
  17. Versteeg H.K., Malalasekera W. An introduction to computational fluid dynamics: the finite volume method. – Pearson education, 2007.
  18. Митрофанова, О.В. Гидродинамика и теплообмен закрученных потоков в каналах с завихрителями. // Теплофизика высоких температур. 2003. Т. 41. № 4. С. 587–633.
  19. Brown C.H., et al. Morphological, biochemical, and functional changes in human platelets subjected to shear stress // The Journal of laboratory and clinical medicine. 1975. Т. 86. – №. 3. С. 462–471.
  20. Sutera S.P. Flow-induced trauma to blood cells // Circulation research. 1977. Т. 41. №. 1. С. 2–8.

Қосымша файлдар

Қосымша файлдар
Әрекет
1. JATS XML
Жүктеу
2. Fig. 1. Geometric approximation of the shape of the left atrium (a) and left ventricle (b) using a statically shaped channel with a longitudinal-radial profile zR2 = const with a concave surface.

Жүктеу (110KB)
Метадеректер
3. Fig. 2. Scheme of the experimental setup (a) and geometric configuration of the numerical model (b, c): a) Ht – height of the water column in a water tank, He – channel height, D – cut-off diameter, d – channel outlet diameter; b) Pin – specified static pressure at the inlet to the channel; (c) α is the angle of fluid velocity relative to the normal of the source surface. Physical modeling for various geometric configurations of the channel was carried out in water, the pressure drop in the channel was created using a column of water in the reservoir, and the swirl of the liquid was organized using a rotating blade. The flow rate was recorded at the output of the experimental channel; the flow lines were visualized by introducing dye into the channel using two needles: measuring and control. In Fig. 2 b, c shows that the numerical simulation was carried out under conditions identical to the experimental ones. The static pressure boundary condition at the channel inlet was varied depending on the modeling requirements. The creation of the azimuthal component of the fluid velocity was carried out by parameterizing the angle α between the velocity vector at the channel entrance to the normal of the input surface.

Жүктеу (299KB)
Метадеректер
4. Fig. 3. Dependence of flow characteristics on the height of the liquid column in the tank for channels of various shapes: Cone – dependence of flow rate for a conical-shaped channel. N = 1, 2, 3 – flow rate dependence for channels zRN = const with the specified exponent. It can be seen that the flow characteristics for channels zRN = const have little difference. The flow characteristics of a conical channel are 4% lower.

Жүктеу (131KB)
Метадеректер
5. Fig. 4. Photographs of streamlines in the absence (a) and in the presence (b) of swirl for the channel zR2 = const. It can be seen that in the absence of the azimuthal component of the flow velocity, the shape of the fluid flow line follows the shape of the channel, while in its presence the shape has the form of a converging spiral.

Жүктеу (126KB)
Метадеректер
6. Fig. 5. Experimental dependence of the number of turns of the line on the height of the measuring needle in the channel slit (a) and photographs of streamlines with the camera positioned above the channel (b-d): The control needle is located on the right, the measuring needle is on the left; a jet of dye from above visualizes the channel boundary. From Fig. 5 it can be seen that the number of turns of the spiral current line increases with increasing height of the measuring needle. According to Fig. 5d shows that visual tracking of the number of turns in the extreme position of the measuring needle is difficult.

Жүктеу (206KB)
Метадеректер
7. Fig. 6. Flow characteristics during numerical and physical modeling for different channel shapes (a), flow velocity structure for conical channels (b) and zR2 = const shape (c): Cone – flow dependence for a conical channel. N = 1, 2, 3 – flow rate dependence for channels zRN = const with the specified exponent. It can be seen that the consumption characteristics are in good agreement. According to Fig. 6b, c it is clear that the decrease in flow rate for channels zRN = const is associated with high hydrodynamic resistance in the zone of maximum narrowing of the channel.

Жүктеу (200KB)
Метадеректер
8. Fig. 7. Dependence of the number of turns of the streamline on the height of the needle during physical and numerical experiments (a), three-dimensional visualization of numerically determined streamlines for heights of 3, 5 and 7 mm (b). It can be seen that the number of turns of spiral current lines in numerical and physical modeling has a satisfactory agreement.

Жүктеу (162KB)
Метадеректер
9. Fig. 8. Cut lines in the axial projection of the numerical model (a), distribution of azimuthal velocity in the radial projection for channels zRN = const at N = 1; 2 and 3 (b, c, d): Blue line – cut at the zero height of the channel, green curve – cut at a height of –3 cm, green curve – cut at a height of –6 cm. It can be seen that a stable vortex structure is formed only in zR2 = const.

Жүктеу (315KB)
Метадеректер
10. Fig. 9. Integral viscous friction losses for channels zRN = const of different orders of degree (a), distribution of viscous friction losses in channels at N = 1; 2 and 3 (b, c, d) It can be seen that in a channel with exponent N = 1, a vortex structure is not formed. In channels with an exponent N = 2, 3, the flow is of a vortex nature, while in the channel N = 2 the loss values in the region of vortex formation are lower than in the case of N = 3.

Жүктеу (196KB)
Метадеректер
11. Fig. 10. Integral of viscous forces along the streamlined surface of the channels zRN = const with an exponent from 1 to 3 (a), distribution of viscous forces along the streamlined surfaces of the channels (b–d). It can be seen that in a channel with an exponent N = 2 there is no clearly defined region of a local maximum of viscous forces along the channel surface, and the integral of viscous forces is minimal among all the cases considered.

Жүктеу (156KB)
Метадеректер
12. Fig. 11. An example of a geometric configuration of a channel with an exponent N = 2 with a concave surface with parameters of the overall surface height h = 5 mm and the ratio of the overall radial size of the concave surface to the entrance radius of the channel k = 0.5.

Жүктеу (38KB)
Метадеректер
13. Fig. 12. Distribution of losses due to viscous friction (a) and structure of azimuthal velocities (b) for different sizes of the concave surface. It can be seen that with an increase in the additional volume created by the concave surface, the flow becomes vortex with the expansion of the vortex formation region, which helps reduce losses in this zone.

Жүктеу (468KB)
Метадеректер
14. Fig. 13. Specific losses due to viscous friction for various overall dimensions of the concave surface: k is the ratio of the overall radius of the concave surface to the maximum radius of the channel, h is the overall height of the concave surface. It can be seen that with an increase in the additional volume created by the concave surface, the integral losses due to viscous friction decrease. At maximum volume, viscous friction losses become lower than in the absence of a concave surface.

Жүктеу (82KB)
Метадеректер
15. Fig. 14. Calculated oscillogram of the static pressure pulse at the entrance to the channel during dynamic modeling, simulating the systole of the cardiac cycle.

Жүктеу (129KB)
Метадеректер
16. Fig. 15. Values of maximum shear stresses arising in a second-order channel depending on time in the presence and absence of a concave surface of maximum volume: solid line - maximum shear stresses in a channel without a concave surface, broken line - maximum shear stresses in a channel in its presence. Red lines mark the time points to represent the shear stress structure. It can be seen that the presence of a concave surface of maximum volume contributes to a twofold reduction in the maximum shear stresses in the channel.

Жүктеу (118KB)
Метадеректер
17. Fig. 16. Values of maximum shear stresses arising in a second-order channel depending on time in the presence and absence of a concave surface of maximum volume: solid line - maximum shear stresses in a channel without a concave surface, broken line - maximum shear stresses in a channel in its presence. Red lines mark the time points to represent the shear stress structure. It can be seen that the presence of a concave surface of maximum volume contributes to a twofold reduction in the maximum shear stresses in the channel.

Жүктеу (276KB)
Метадеректер
18. Fig. 17. Maximum values of time integrals of shear stresses in channels in the presence and absence of a concave surface (a), structure of the integral of the action of shear stresses in these channels (b, c). It can be seen that the presence of a concave surface makes it possible to reduce the effect of shear stresses by half in the near-wall zone, and also significantly reduces their effect in the area of formation of vortex motion.

Жүктеу (207KB)
Метадеректер

© Russian Academy of Sciences, 2024

Согласие на обработку персональных данных с помощью сервиса «Яндекс.Метрика»

1. Я (далее – «Пользователь» или «Субъект персональных данных»), осуществляя использование сайта https://journals.rcsi.science/ (далее – «Сайт»), подтверждая свою полную дееспособность даю согласие на обработку персональных данных с использованием средств автоматизации Оператору - федеральному государственному бюджетному учреждению «Российский центр научной информации» (РЦНИ), далее – «Оператор», расположенному по адресу: 119991, г. Москва, Ленинский просп., д.32А, со следующими условиями.

2. Категории обрабатываемых данных: файлы «cookies» (куки-файлы). Файлы «cookie» – это небольшой текстовый файл, который веб-сервер может хранить в браузере Пользователя. Данные файлы веб-сервер загружает на устройство Пользователя при посещении им Сайта. При каждом следующем посещении Пользователем Сайта «cookie» файлы отправляются на Сайт Оператора. Данные файлы позволяют Сайту распознавать устройство Пользователя. Содержимое такого файла может как относиться, так и не относиться к персональным данным, в зависимости от того, содержит ли такой файл персональные данные или содержит обезличенные технические данные.

3. Цель обработки персональных данных: анализ пользовательской активности с помощью сервиса «Яндекс.Метрика».

4. Категории субъектов персональных данных: все Пользователи Сайта, которые дали согласие на обработку файлов «cookie».

5. Способы обработки: сбор, запись, систематизация, накопление, хранение, уточнение (обновление, изменение), извлечение, использование, передача (доступ, предоставление), блокирование, удаление, уничтожение персональных данных.

6. Срок обработки и хранения: до получения от Субъекта персональных данных требования о прекращении обработки/отзыва согласия.

7. Способ отзыва: заявление об отзыве в письменном виде путём его направления на адрес электронной почты Оператора: info@rcsi.science или путем письменного обращения по юридическому адресу: 119991, г. Москва, Ленинский просп., д.32А

8. Субъект персональных данных вправе запретить своему оборудованию прием этих данных или ограничить прием этих данных. При отказе от получения таких данных или при ограничении приема данных некоторые функции Сайта могут работать некорректно. Субъект персональных данных обязуется сам настроить свое оборудование таким способом, чтобы оно обеспечивало адекватный его желаниям режим работы и уровень защиты данных файлов «cookie», Оператор не предоставляет технологических и правовых консультаций на темы подобного характера.

9. Порядок уничтожения персональных данных при достижении цели их обработки или при наступлении иных законных оснований определяется Оператором в соответствии с законодательством Российской Федерации.

10. Я согласен/согласна квалифицировать в качестве своей простой электронной подписи под настоящим Согласием и под Политикой обработки персональных данных выполнение мною следующего действия на сайте: https://journals.rcsi.science/ нажатие мною на интерфейсе с текстом: «Сайт использует сервис «Яндекс.Метрика» (который использует файлы «cookie») на элемент с текстом «Принять и продолжить».