3.1. Трехмерные уравнения (неразрывности и импульса – Навье-Стокса)
Уравнение Навье-Стокса для несжимаемой жидкости имеет вид:
3.2. Одномерные уравнения
Сила трения в расчете на единицу длины трубки и единицу сечения fr = - τL/s
Параболический профиль, течение Пуазейля τL = 8πηu
Турбулентный поток τL = ρKtu2
Уравнение неразрывности
∂S/∂t + ∂(uS)/ ∂x = 0
Закон трубки
S = S0 + ∂S/∂P * dP
3.2. Уравнения для линейного, невязкого варианта
Уравнение движения
Уравнение неразрывности ∂S/∂t = - S0 ∂u/ ∂x
3.3. Вывод формулы Моэнса-Кортевега
Уравнение неразрывности + закон трубки
(∂S/∂P)0 ∂P/∂t + S0 ∂u/∂x = 0 t
баланс импульса
∂u/∂t = - 1/ ρ ∂P/∂x x
|
|
Получаем волновое уравнение
∂2P/∂t2 = S0 /(ρ (∂S/∂P)0) ∂2P/∂x2
с = (S0 /ρ * 1/(∂S/∂P)0)1/2 = (Eh/ ρd) 1/2
3.4. Рассматривали следующую упрощенную ситуацию:
1) Предположение о линейности. Есть нелинейные эффекты:
- выбросили квадрат скорости (конвективные ускорения), это возможно, если они <<∂u/∂t, для этого u/c << 1. реально u/c~0.4
- нелинейность S-P. dP/P ~ 0.2. Искажает форму, но мало влияет на скорость.
2) Нулевая вязкость: если α >> 1, то можно пренебрегать вязкостью. Аорта человека α ≈ 20, α = (d/2)*(ω/ ν)1/2
т.е. можно.
3) Стенка чисто упругая
4) Жидкость несжимаемая
3.5. Роль отражений в точках ветвления, импеданс согласованных сосудов
P = zQ
z = ρc/S
Pотр/Pпад = (Z0-1 - (Z1-1 + Z2-1))/(Z0-1 + (Z1-1 + Z2-1)) = R
3.6. Нестационарное течение в очень длинной трубке
P=P0*sin ωt
Стационарное течение – параболический профиль (течение Пуазейля).
Нестационарное P=P0*sin ωt
α = (d/2)*(ω/ ν)1/2
J. R. Womersley. Method for the calculation of velocity, rate of flow and viscous drag in arteries when the pressure gradient is known // J Physiol. 1955 Mar 28; 127(3): 553–563.
Каро и соавт., 1981.