Все выпуски

Изучение физиологических и патофизиологических процессов, связанных с системой кровообращения, является на сегодняшний день актуальной темой многих исследований. В данной работе рассматривается ряд подходов к математическому моделированию кровотока, основанных на пространственном осреднении и/или использующих стационарное приближение. Обсуждаются допущения и предположения, ограничивающие область применения моделей такого рода. Приводятся наиболее распространенные математические постановки задач и кратко описываются методы их численного решения. В первой части обсуждаются модели, основанные на полном пространственном осреднении и/или использующие стационарное приближение. Один из наиболее распространенных на сегодняшний день подходов состоит в проведении аналогий между течением вязкой несжимаемой жидкости в эластичных трубках и электрическим током в цепи. Такие модели используются не только сами по себе, но и как способ постановки граничных условий в моделях, учитывающих одномерную или трехмерную пространственную зависимость переменных. Динамические, полностью осредненные по пространству модели позволяют описывать динамику кровотока на достаточно больших временных интервалах, равных длительности десятков сердечных циклов и более. Далее рассмотрены стационарные модели основанные как на полностью осредненном, так и на двухмерном подходе. Такие модели могут быть использованы для моделирования кровотока в микроциркуляторном русле. Во второй части обсуждаются модели, основанные на одномерном осреднении параметров кровотока. Преимущество данного подхода также состоит в невысоких, по сравнению с трехмерным моделированием, требованиях к вычислительным ресурсам и возможности охвата всех достаточно крупных кровеносных сосудов в организме. Модели данного типа позволяют рассчитывать параметры кровотока в каждом сосуде сосудистой сети, включенной в модель. Структура и параметры такой сети могут быть заданы как на основе данных литературы, так и с помощью методов сегментации медицинских данных. Основными и весьма существенными предположениями при выводе одномерных уравнений из уравнений Навье – Стокса с помощью асимптотического анализа или их интегрирования по объему являются радиальная симметрия течения и постоянство формы профиля скорости в поперечном сечении. Существующие в настоящее время работы, посвященные валидации одномерных моделей, их сравнению между собой и с данными клинических исследований, позволяют говорить об успешности данного подхода и подтверждают возможность его использования в медицинской практике. Одномерные модели позволяют описывать такие динамические явления, как распространение пульсовой волны и звуки Короткова. В этом приближении могут быть учтены такие факторы, как действие на кровоток силы тяжести, действие на стенки сосудов силы сжатия мышц, регуляторные и ауторегуляторные эффекты.

The study of the physiological and pathophysiological processes in the cardiovascular system is one of the important contemporary issues, which is addressed in many works. In this work, several approaches to the mathematical modelling of the blood flow are considered. They are based on the spatial order reduction and/or use a steady-state approach. Attention is paid to the discussion of the assumptions and suggestions, which are limiting the scope of such models. Some typical mathematical formulations are considered together with the brief review of their numerical implementation. In the first part, we discuss the models, which are based on the full spatial order reduction and/or use a steady-state approach. One of the most popular approaches exploits the analogy between the flow of the viscous fluid in the elastic tubes and the current in the electrical circuit. Such models can be used as an individual tool. They also used for the formulation of the boundary conditions in the models using one dimensional (1D) and three dimensional (3D) spatial coordinates. The use of the dynamical compartment models allows describing haemodynamics over an extended period (by order of tens of cardiac cycles and more). Then, the steady-state models are considered. They may use either total spatial reduction or two dimensional (2D) spatial coordinates. This approach is used for simulation the blood flow in the region of microcirculation. In the second part, we discuss the models, which are based on the spatial order reduction to the 1D coordinate. The models of this type require relatively small computational power relative to the 3D models. Within the scope of this approach, it is also possible to include all large vessels of the organism. The 1D models allow simulation of the haemodynamic parameters in every vessel, which is included in the model network. The structure and the parameters of such a network can be set according to the literature data. It also exists methods of medical data segmentation. The 1D models may be derived from the 3D Navier – Stokes equations either by asymptotic analysis or by integrating them over a volume. The major assumptions are symmetric flow and constant shape of the velocity profile over a cross-section. These assumptions are somewhat restrictive and arguable. Some of the current works paying attention to the 1D model’s validation, to the comparing different 1D models and the comparing 1D models with clinical data. The obtained results reveal acceptable accuracy. It allows concluding, that the 1D approach can be used in medical applications. 1D models allow describing several dynamical processes, such as pulse wave propagation, Korotkov’s tones. Some physiological conditions may be included in the 1D models: gravity force, muscles contraction force, regulation and autoregulation.

Тромбоз воротной вены (ТВВ) является значимым осложнением как в претрансплантационном, так и в послеоперационном периоде трансплантации печени. Многофакторная этиология ТВВ, парадоксальное состояние гемостаза при циррозе печени и ограниченная информативность стандартных коагулологических тестов определяют необходимость разработки формализованных моделей для оценки риска тромбообразования.

Цель работы — на основе математической модели системы свертывания крови исследовать влияние гидродинамических условий и уровней факторов коагуляции на возможность формирования тромба в воротной вене.

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

Показан пороговый характер генерации тромбина в зависимости от скорости кровотока. При превышении критической скорости фаза инициации не переходит в фазу амплификации, что соответствует физиологическим условиям, препятствующим тромбообразованию. Продемонстрировано, что при пониженной концентрации фибриногена (характерной для дисфункции печени) порог скорости, выше которого тромб не формируется, возрастает, что означает повышенную чувствительность системы к стазу. Дефицит протеина C в моделируемых условиях оказывает минимальное влияние на динамику тромбообразования. Результаты моделирования качественно согласуются с клиническими данными по распределению пациентов с ТВВ по уровню фибриногена ($n$ = 932, НИИ им. Склифосовского).

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

Portal vein thrombosis (PVT) is a significant complication during both the pre-transplant and postoperative periods of liver transplantation. The multifactorial etiology of PVT, the paradoxical hemostatic state in liver cirrhosis and the limited usefulness of standard coagulation tests highlight the necessity of formalized models to assess thrombosis risk.

Objective: Using a mathematical model of the blood coagulation system, investigate the influence of hemodynamic conditions and coagulation factor levels on the likelihood of thrombus formation in the portal vein.

The portal vein is modelled as a flow-through reactor with rapid convective mixing. The mathematical model is based on Panteleev et al. (2010) detailed kinetic scheme, incorporating equations for the extrinsic pathway of coagulation activation, positive and negative feedback loops, and inhibition of active factors. The resulting system of ordinary differential equations was integrated using a one-stage Rosenbrock method with complex coefficients.

Thrombin generation was shown to exhibit threshold dependence on blood flow velocity. Above a critical velocity, the initiation phase does not transition to the amplification phase. This corresponds to physiological conditions that prevent thrombus formation. We demonstrated that, at reduced fibrinogen concentrations characteristic of hepatic dysfunction, the critical velocity threshold above which thrombus formation is suppressed increases. This indicates the system’s heightened susceptibility to stasis. Protein C deficiency had minimal effect on thrombogenesis dynamics under the modeled conditions. The modeling results qualitatively agree with clinical data on fibrinogen distribution in PVT patients ($n$ = 932, Sklifosovsky Research Institute).

В данной статье рассмотрены особенности численных решений некоторых задач для кноидальных волн, которые являются периодическими решениями классического уравнения Кортевега – де Фриза типа бегущей волны. Точные решения, описывающие эти волны, получены путемс ведения автоволновым приближением уравнения Кортевега – де Фриза к обыкновенным дифференциальным уравнениям сначала третьего, затем второго и, наконец, первого порядков. Обращение к числовому примеру показывает, что полученные такимо бразом обыкновенные дифференциальные уравнения не являются равносильными. Сформулированная и доказанная в настоящей статье теорема и замечание к ней показывают, что множество решений уравнения третьего порядка самое широкое и в качестве подмножеств включает в себя множества решений уравнений первого и второго порядков, которые в свою очередь равносильными не являются. Полученное автоволновым приближением обыкновенное дифференциальное уравнение первого порядка является источником для нахождения точных формул для описания кноидальной волны (периодического решения) и солитона (уединенной волны). Несмотря на это, с вычислительной точки зрения это уравнение является самым неудобным. Для этого уравнения не выполняется условие Липшица по искомой функции в окрестности постоянных решений. Отсюда теорема о существовании и единственности решения задачи Коши для обыкновенного дифференциального уравнения первого порядка не является справедливой. В частности, в стационарных точках нарушается единственность решения задачи Коши. Поэтому для обыкновенного дифференциального уравнения первого порядка, полученного из уравнения Кортевега – де Фриза, и в случае кноидальной волны, и в случае солитона задачу Коши нельзя ставить в точках экстремума. Начальное условие может быть поставлено лишь в точке убывания или роста, а отрезок численного решения необходимо выбрать так, чтобы он лежал между соседними точками экстремума. Но для уравнений второго и третьего порядков начальные условия можно ставить как в точках убывания или роста, так и в точках экстремума. При этом отрезок для численного решения сильно расширяется и наблюдается периодичность. Для решений этих обыкновенных уравнений изучаются постановки задач Коши, проводится сравнение полученных результатов с точными решениями и между собой. Показана численная реализация перерождения кноидальной волны в солитон. Результаты статьи имеют гемодинамическую интерпретацию пульсационного течения кровотока в цилиндрическом кровеносном сосуде, состоящем из упругих колец.

This article discusses the features of the numerical solutions of some problems for cnoidal waves, which are periodic solutions of the classical Korteweg – de Vries equation of the traveling wave type. Exact solutions describing these waves were obtained by communicating the autowave approximation of the Korteweg – de Vries equation to ordinary functions of the third, second, and finally, first orders. Referring to a numerical example shows that in this way ordinary differential equations are not equivalent. The theorem formulated and proved in this article and the remark to it include the set of solutions of the first and second order, which, in their ordinal, are not equivalent. The ordinary differential equation of the first order obtained by the autowave approximation for the description of a cnoidal wave (a periodic solution) and a soliton (a solitary wave). Despite this, from a computational point of view, this equation is the most inconvenient. For this equation, the Lipschitz condition for the sought-for function is not satisfied in the neighborhood of constant solutions. Hence, the existence theorem and the unique solutions of the Cauchy problem for an ordinary differential equation of the first order are not valid. In particular, the uniqueness of the solution to the Cauchy problem is violated at stationary points. Therefore, for an ordinary differential equation of the first order, obtained from the Korteweg – de Vries equation, both in the case of a cnoidal wave and in the case of a soliton, the Cauchy problem cannot be posed at the extremum points. The first condition can be a set position between adjacent extremum points. But for the second, third and third orders, the initial conditions can be set at the growth points and at the extremum points. In this case, the segment for the numerical solution greatly expands and periodicity is observed. For the solutions of these ordinary solutions, the statements of the Cauchy problems are studied, and the results are compared with exact solutions and with each other. A numerical realization of the transformation of a cnoidal wave into a soliton is shown. The results of the article have a hemodynamic interpretation of the pulsating blood flow in a cylindrical blood vessel consisting of elastic rings.

В данной работе описана методика генерации синтетической базы данных стенозов, состоящей из 1620 записей. Каждая запись представляет собой результаты численного эксперимента по моделированию трехмерного течения вязкой несжимаемой жидкости через трубку с переменным сечением: перепад давлений, средний поток, усредненная по сечению скорость кровотока на входе в трубку, максимальная степень сужения стеноза, длина стеноза, асимметрия стеноза, радиус трубки, число Рейнольдса. База данных валидировалась путем сравнения с другими моделями (с эластичными стенками) и стендовыми экспериментами и показала отклонение перепадов давлений не выше 4%. База данных синтетических стенозов использовалась для обучения физически информированной нейронной сети для быстрой оценки перепада давления по четырем ключевым входным показателям: число Рейнольдса, длина стеноза, степень стеноза, степень асимметрии стеноза. Физическая информированность достигалась за счет введения штрафов в функцию потерь за отсутствие положительного перепада давления и за отсутствие монотонности перепада давления по входным параметрам. Физически информированная нейронная сеть показала более высокую точность на гемодинамических значимых стенозах при тестировании на валидационной выборке и на новых стенозах, не представленных в базе данных. Средняя относительная ошибка на стенозах длиной в 8 радиусов здорового сосуда составила 6% для физически информированной сети и 13% для классической нейронной сети. Ошибки на коротких стенозах длиной в 4 радиуса оказались почти одинаковы: 9,5% для физически информированной сети и 10% для классической нейронной сети. Разработанный метод функциональной оценки гемодинамической значимости стенозов может использоваться как самостоятельный инструмент для клинической оценки стенозов и как компонент сетевых моделей кровотока. Наибольшую актуальность подход приобретает при моделировании многососудистых поражений, которые доминируют в клинической практике. Ключевое преимущество метода заключается в физической корректности результатов и точности, сопоставимой с классическим моделированием, при значительно меньших вычислительных затратах.

This paper describes a method for generating a synthetic database of stenoses, consisting of 1620 entries. Each entry represents the results of a numerical experiment simulating the three-dimensional flow of a viscous incompressible fluid through a tube with a variable cross-section: pressure drop, mean flow rate, cross-sectionally averaged inlet blood flow velocity, maximum stenosis severity, stenosis length, stenosis asymmetry, tube radius, and Reynolds number. The database was validated by comparison with other models (with elastic walls) and bench experiments, showing a deviation in pressure drops of no more than 4%. The synthetic stenosis database was used to train a physics-informed neural network for the rapid estimation of pressure drop based on four key input parameters: Reynolds number, stenosis length, stenosis severity, and stenosis asymmetry coefficient. The physics-informed aspect was achieved by introducing penalties into the loss function for the absence of a positive pressure drop and for the lack of monotonicity of the pressure drop with respect to the input parameters. The physics-informed neural network demonstrated higher accuracy on hemodynamically significant stenoses when tested on a validation set and on new stenoses not represented in the database. The mean relative error for stenoses with a length of 8 healthy vessel radii was 6% for the physics-informed network and 13% for a classical neural network. The errors for short stenoses with a length of 4 radii were nearly identical: 9.5% for the physics-informed network and 10% for the classical neural network. The developed method for the functional assessment of the hemodynamic significance of stenoses can be used both as a standalone tool for clinical stenosis evaluation and as a component of network blood flow models. The approach becomes most relevant when modeling multi-vessel disease, which is predominant in clinical practice. The key advantage of the method lies in the physical correctness of the results and accuracy comparable to classical modeling, but with significantly lower computational costs.

Разработана численная модель распространения возмущений скорости тока крови и давления по кровеносному сосуду с тромбом, расположенным в точке венозной бифуркации, и их влияния на динамику тромба. Модель построена в акустическом (линейном) приближении. Результаты расчетов позволят определить условия возникновения резонансных колебаний тромба, которые могут привести к его отрыву и тромбоэмболии.

In this paper a numerical model for blood flow through a venous bifurcation with an obliterating clot is investigated. We studied propagation of perturbations of blood flow velocity and perturbations of pressure inside the vein. The model is built in acoustic (linear) approximation. Computational results reveal conditions for clot resonance oscillation, which can cause its detachment and thromboembolism.

В настоящей работе рассматривается математическая модель динамики клеточной ткани. В первой части дается вывод модели, основные положения и постановка задачи. Во второй части итоговая система исследуется численно и приводятся результаты моделирования. Постулируется, что клеточная ткань есть трехфазная среда, которая состоит из твердого скелета (представляющего собой внеклеточный матрикс), клеток и внеклеточной жидкости. Ко всему прочему учитывается наличие питательных веществ в ткани. В основу модели положены уравнения сохранения массы с учетом обмена масс, уравнения сохранения импульса для каждой фазы, а также уравнение диффузии для питательных веществ. В уравнении, описывающем клеточную фазу, также учитывается слагаемое, описывающее химическое воздействие на ткань, которое называется хемотаксисом — движением клеток, вызванным градиентом концентрации химических веществ. Исходная система уравнений сводится к системе трех уравнений для нахождения пористости, насыщенности клеток и концентрации питательных веществ. Данные уравнения дополняются начальными и краевыми условиями. В одномерном случае в начальный момент времени задается распределение пористости, концентрации клеточной фазы и питательных веществ. На левой границе задана постоянная концентрация питательных веществ, что соответствует, например, поступлению кислорода из сосуда, а также поток концентрации клеток на ней равен нулю. На правой границе рассматриваются два типа условий: первое — условие непроницаемости правой границы, второе — условие постоянной концентрации клеточной фазы и нулевой поток концентрации питательных веществ. В обоих случаях условия для матрикса и внеклеточной жидкости одинаковы, предполагается наличие источника питательных веществ (кровеносного сосуда) на левой границе области моделирования. В результате моделирования было выявлено, что хемотаксис оказывает значительное влияние на рост ткани. При отсутствии хемотаксиса зона уплотнения распространяется на всю область моделирования, но при увеличении влияния хемотаксиса на ткань образуется область деградации, в которой концентрация клеток становится ниже начальной.

In this paper, a mathematical model of cellular tissue dynamics is considered. The first part gives the conclusion of the model, the main provisions and the formulation of the problem. In the second part, the final system is investigated numerically and the simulation results are presented. It is postulated that cellular tissue is a three-phase medium that consists of a solid skeleton (which is an extracellular matrix), cells and extracellular fluid. In addition, the presence of nutrients in the tissue is taken into account. The model is based on the equations of conservation of mass, taking into account mass exchange, the equations of conservation of momentum for each phase, as well as the diffusion equation for nutrients. The equation describing the cellular phase also takes into account the term describing the chemical effect on the tissue, which is called chemotaxis — the movement of cells caused by a gradient in the concentration of chemicals. The initial system of equations is reduced to a system of three equations for finding porosity, cell saturation and nutrient concentration. These equations are supplemented by initial and boundary conditions. In the one-dimensional case, the distribution of porosity, concentration of the cell phase and nutrients is set at the initial moment of time. A constant concentration of nutrients is set on the left border, which corresponds, for example, to the supply of oxygen from the vessel, as well as the flow of cell concentration on it is zero. Two types of conditions are considered at the right boundary: the first is the condition of impermeability of the right boundary, the second is the condition of constant concentration of the cell phase and zero flow of nutrient concentration. In both cases, the conditions for the matrix and extracellular fluid are the same, it is assumed that there is a source of nutrients (blood vessel) on the left border of the modeling area. As a result of modeling, it was revealed that chemotaxis has a significant effect on tissue growth. In the absence of chemotaxis, the compaction zone extends to the entire modeling area, but with an increase in the effect of chemotaxis on the tissue, a degradation area is formed in which the concentration of cells becomes lower than the initial one.

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

Проводились расчеты в области, представляющей фрагмент КТ изображения воротной вены конкретного пациента, не содержащий сосудистых бифуркаций. В качестве граничных условий для потока использована интерполяция данных доплерографии пациентов, опубликованных в печати.

Расчеты проводились с использованием программного комплекса индустриальной гидродинамики FlowVision. Проведено сравнение расчетов при использовании модели идеальной жидкости и модели Куемады вязкопластического течения. Расчеты проводились для разных значений коэффициента сопротивления, соответствующих физиологической норме, и при повышенном значении коэффициента сопротивления.

При нормальном значении коэффициента сопротивления течение в воротной вене характеризуется сильным конвективным перемешиванием.

В результате расчетов получено, что при использовании модели Куемады течение в воротной вене стратифицировано. Характер стратификации зависит от гематокрита. При нормальном значении коэффициента сопротивления при уменьшении скорости формируется пластическое ядро течения. При повышенном значении коэффициента сопротивления при инверсии потока также формируется ядро течения. При инверсии потока пластическое ядро продолжает движение в прямом направлении, в то время как пристеночные слои жидкости начинают двигаться в обратном направлении.

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

A mathematical model of fluid movement in the portal vein is considered. The fundamental circumstance determining the whole variety of portal hemodynamic phenomena is the absence of a valve apparatus. The flow direction is solely a function of the pressure gradient, and, therefore, it is fundamentally reversible when the boundary conditions of the system change.

Calculations were performed in the area representing a CT image fragment of the portal vein of a particular patient, which does not contain vascular bifurcations. The interpolation of patient Dopplerography data published in the press was used as boundary conditions for the flow. The calculations were carried out using the FlowVision industrial hydrodynamics software package. A comparison of calculations using the ideal fluid model and the Kuemada model of viscoplastic flow is carried out.

Calculations were performed for different values of the resistance coefficient corresponding to the physiological norm and with an increased value of the resistance coefficient.

At a normal value of the resistance coefficient, the flow in the portal vein is characterized by strong convective mixing.

As a result of calculations, it was found that when using the Kuemada model, the flow in the portal vein is stratified. The nature of the stratification depends on the hematocrit. At a normal value of the resistance coefficient, a plastic core of the flow is formed as the velocity decreases. With an increased value of the resistance coefficient during flow inversion, the flow core is also formed. When the flow is inverted, the plastic core continues to move in the forward direction, while the wall layers of the liquid begin to move in the opposite direction.

It is noted that in order to obtain correct modeling results, it is necessary to refine the blood composition of the portal vein and, possibly, refine the rheological model. Such information can be obtained both from comparing simulation data with clinical data, and by laboratory examination of portal vein blood.

В работе рассмотрены модели формирования тромбоцитарного тромба в потоке плазмы крови в цилиндрическом сосуде, основанные на уравнении типа «адвекция–диффузия» и уравнении Фоккера–Планка. Приведено сравнение результатов расчетов на основе этих моделей. Рассмотренные модели демонстрируют качественно схожее поведение на начальном этапе формирования тромба. При детальном исследовании возникновения крупных сгустков необходимо уточнение моделей.

The paper considers models of platelet thrombus formation in blood plasma flow in a cylindrical vessel, based on the “advection–diffusion” equation and the Fokker–Planck equation. The comparison of the results of calculations based on these models is given. Considered models show qualitatively similar behavior at the initial stage of thrombus formation. А detailed investigation of large clots requires models’ improvement.

Целью проведенного исследования была разработка математической модели пульсирующего течения крови по участку аорты, включающему восходящий отдел, дугу аорты с ее ответвлениями и верхнюю часть нисходящего отдела. Поскольку при прохождении пульсовой волны деформации этой наиболее твердой части аорты малы, то при построении механической модели ее стенки считались абсолютно твердыми. В статье приводится описание внутренней структуры крови и ряда внутриструктурных эффектов. Этот анализ показывает, что кровь, которая по существу является суспензией, можно рассматривать только как неньютоновскую жидкость. Кроме того, кровь можно считать жидкостью только в кровеносных сосудах, диаметр которых намного больше характерного размера клеток крови и их агрегатных образований. В качестве неньютоновской жидкости была выбрана вязкая жидкость со степенным законом связи напряжения со скоростью деформации. Этот закон позволяет описывать поведение не только жидкостей, но и суспензий. При постановке граничного условия на входе в аорту, отражающего пульсирующий характер течения крови, было решено не ограничиваться заданием совокупного потока крови, который не дает представления о пространственном распределении скорости по поперечному сечению. В связи с этим было предложено моделировать огибающую поверхность этого пространственного распределения частью параболоида вращения с фиксированным радиусом основания и высотой, которая меняется во времени от нуля до максимального значения скорости. Для граничного условия на стенке сосуда предлагается использовать условие полупроскальзывания. Это связано с тем, что клетки крови, в силу своих электрохимических свойств, не прилипают к внутреннему слою сосуда. На внешних концах аорты и ее ответвлений задавалась величина давления. Для выполнения вычислений была построена геометрическая модель рассматриваемой части аорты с ответвлениями, на которую была нанесена тетраэдальная сетка с общим числом элементов 9810. Вычисления производились методом конечных элементов с шагом по времени 0.01 с с использованием пакета ABAQUS. В результате было получено распределение скоростей и давления на каждом шаге по времени. В областях ветвления сосудов было обнаружено вре́менное наличие вихрей и обратных течений. Они зарождались через 0.47 с от начала пульсового цикла и исчезали спустя 0.14 с.

The purpose of research was to develop a mathematical model for pulsating blood flow in the part of aorta with their branches. Since the deformation of this most solid part of the aorta is small during the passage of the pulse wave, the blood vessels were considered as non-deformable curved cylinders. The article describes the internal structure of blood and some internal structural effects. This analysis shows that the blood, which is essentially a suspension, can only be regarded as a non-Newtonian fluid. In addition, the blood can be considered as a liquid only in the blood vessels, diameter of which is much higher than the characteristic size of blood cells and their aggregate formations. As a non-Newtonian fluid the viscous liquid with the power law of the relationship of stress with shift velocity was chosen. This law can describe the behaviour not only of liquids but also dispersions. When setting the boundary conditions at the entrance into aorta, reflecting the pulsating nature of the flow of blood, it was decided not to restrict the assignment of the total blood flow, which makes no assumptions about the spatial velocity distribution in a cross section. In this regard, it was proposed to model the surface envelope of this spatial distribution by a part of a paraboloid of rotation with a fixed base radius and height, which varies in time from zero to maximum speed value. The special attention was paid to the interaction of blood with the walls of the vessels. Having regard to the nature of this interaction, the so-called semi-slip condition was formulated as the boundary condition. At the outer ends of the aorta and its branches the amounts of pressure were given. To perform calculations the tetrahedral computer network for geometric model of the aorta with branches has been built. The total number of meshes is 9810. The calculations were performed with use of the software package ABACUS, which has also powerful tools for creating geometry of the model and visualization of calculations. The result is a distribution of velocities and pressure at each time step. In areas of branching vessels was discovered temporary presence of eddies and reverse currents. They were born via 0.47 s from the beginning of the pulse cycle and disappeared after 0.14 s.

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

An approximate mathematical model of blood flow in an axisymmetric blood vessel is studied. Such a vessel is understood as an infinitely long circular cylinder, the walls of which consist of elastic rings. Blood is considered as an incompressible fluid flowing in this cylinder. Increased pressure causes radially symmetrical stretching of the elastic rings. Following J. Lamb, the rings are located close to each other so that liquid does not flow between them. To mentally realize this, it is enough to assume that the rings are covered with an impenetrable film that does not have elastic properties. Only rings have elasticity. The considered model of blood flow in a blood vessel consists of three equations: the continuity equation, the law of conservation of momentum and the equation of state. An approximate procedure for reducing the equations under consideration to the Korteweg – de Vries (KdV) equation is considered, which was not fully considered by J. Lamb, only to establish the dependence of the coefficients of the KdV equation on the physical parameters of the considered model of incompressible fluid flow in an axisymmetric vessel. From the KdV equation, by a standard transition to traveling waves, ODEs of the third, second and first orders are obtained, respectively. Depending on the different cases of arrangement of the three stationary solutions of the first-order ODE, a cnoidal wave and a soliton are standardly obtained. The main attention is paid to an unbounded periodic solution, which we call a degenerate cnoidal wave. Mathematically, cnoidal waves are described by elliptic integrals with parameters defining amplitudes and periods. Soliton and degenerate cnoidal wave are described by elementary functions. The hemodynamic meaning of these types of decisions is indicated. Due to the fact that the sets of solutions to first-, second- and third-order ODEs do not coincide, it has been established that the Cauchy problem for second- and third-order ODEs can be specified at all points, and for first-order ODEs only at points of growth or decrease. The Cauchy problem for a first-order ODE cannot be specified at extremum points due to the violation of the Lipschitz condition. The degeneration of the cnoidal wave into a degenerate cnoidal wave, which can lead to rupture of the vessel walls, is numerically illustrated. The table below describes two modes of approach of a cnoidal wave to a degenerate cnoidal wave.