ИНСТРУМЕНТ C-31 · v1.0.0 · БЕСПЛАТНО · БЕЗ РЕГИСТРАЦИИ

C-31 · ТЕПЛООБМЕН & ГИДРОДИНАМИКА · СЕРДЕЧНО-СОСУДИСТАЯ СИСТЕМА

Пульс в артерии

Как кровь движется по артерии за один сердечный цикл: профиль скорости, касательное напряжение на стенке и моменты, когда оно меняет направление. Точное решение Уомерсли для аорты, артерии меньшего размера или любого введённого вами сосуда.

БЕСПЛАТНО · БЕЗ РЕГИСТРАЦИИ v1.0.0 · проверено по 10 опубликованным значениям
ИСХОДНЫЕ ДАННЫЕточное решение
1 · Течение за сердечный цикл
2 · Сосуд и кровь

По умолчанию: восходящая аорта, смоделированная De Tullio et al. (2009): 27 mm, 69,3 удара в минуту (866 ms), 0,81 m/s на пике. Попробуйте артерию поменьше или более частое сердцебиение и посмотрите, как изменится профиль.

3 · Воспроизведение
СКОРОСТЬ ПО СЕЧЕНИЮ СОСУДА

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

ЗА ОДИН ЦИКЛрасход · касательное напряжение на стенке
ЧИСЛО УОМЕРСЛИ20.0α = R·√(ω/ν)
СКОРОСТЬ НА ОСИ / СРЕДНЯЯ В ПИКЕ1.21 ×Пуазейль: 2
ПИК КАСАТЕЛЬНОГО НАПРЯЖЕНИЯ НА СТЕНКЕ6.40 Pa0.83 Pa если стационарно
ИНДЕКС КОЛЕБАТЕЛЬНОГО СДВИГА0.4370 в одну сторону · 0,5 вперёд и назад

Что показывает этот цикл

  • Число Уомерсли 20.0: инерция управляет ядром потока. Профиль плоский почти по всему сосуду, а вязкость действует лишь в слое толщиной около 0.95 mm у стенки.
  • На пике расхода скорость на оси в 1.21× больше средней; течение Пуазейля дало бы 2.
  • Пиковое касательное напряжение на стенке 6.40 Pa против 0.83 Pa, если считать течение стационарным в каждый момент: разница +674 %.
  • Касательное напряжение на стенке опережает поток: на 42° на частоте сердечных сокращений (0° для стационарного течения, до 45° при больших числах Уомерсли).
  • Осреднённое по времени касательное напряжение на стенке 0.15 Pa; индекс колебательного сдвига 0.437 (0 = всегда в одну сторону, 0,5 = назад столько же, сколько вперёд).
  • В течение 70 % цикла кровь у стенки течёт назад, тогда как ядро потока ещё движется вперёд.
ПРОВЕРЕНО ПО ОПУБЛИКОВАННЫМ ЗНАЧЕНИЯМ · ЗАПУСКАЕТСЯ В НАШЕМ НАБОРЕ ТЕСТОВ
ЧтоОпубликованоЭтот калькулятор
Waite & Fine §7.7, гармоника 1: 1 − 2J₁/(ΛJ₀)0.5783 · 0.32920.5783 · 0.3292согласуется
Waite & Fine §7.7, гармоника 2: 1 − 2J₁/(ΛJ₀)0.7088 · 0.24670.7088 · 0.2467согласуется
De Tullio et al. 2009: средний расход аортальной волны5 L/min5.03 L/minсогласуется
Waite & Fine §7.7: стационарная скорость на оси0.35525 m/s0.35536 m/sсогласуется
Womersley 1955, бедренная артерия собаки: M′10 при α = 4.720.74360.7437согласуется
Womersley 1955, бедренная артерия собаки: M′10 при α = 5.780.78390.7839согласуется
Womersley 1955, бедренная артерия собаки: M′10 при α = 6.670.80960.8096согласуется
Womersley 1955, бедренная артерия собаки: M′10 при α = 7.460.82780.8278согласуется
Womersley 1955, бедренная артерия собаки: M′10 при α = 8.170.84160.8414согласуется
Womersley 1955, бедренная артерия собаки: M′10 при α = 3.340.66510.6554см. примечание
  • Значение Womersley для α = 3,34 напечатано как 0,6651; та же формула и scipy дают 0,6554, тогда как остальные пять значений совпадают до 0,0002. Мы считаем это опечаткой.
  • Малое α: парабола стационарного течения, касательное напряжение на стенке синфазно с потоком. Большое α: плоское ядро, касательное напряжение на стенке опережает поток на 45°. Оба предела входят в набор тестов.
  • Только для обучения и сравнения. Не для диагностики и лечения.
  • Упругие стенки, ветвления и развивающееся течение за клапаном требуют полного моделирования, для которого мы и создали Tessera. Подробнее о Tessera →
Бесплатно с учётной записью · Войти

КАК ЭТО РАССЧИТЫВАЕТСЯ · С ВАШИМИ ЧИСЛАМИ

Шаг за шагом, с вашими числами

  1. 01
    Число Уомерслиα = R·√(ω/ν) = 13.50 mm × √(2π / 0.866 s / 3.302 mm²/s) = 20.01Отношение радиуса сосуда к глубине, на которую за один цикл проникает влияние трения о стенку.
  2. 02
    Вязкий слой у стенкиδ = √(2ν/ω) = 0.954 mmСлой Стокса основной гармоники; у высших гармоник он тоньше.
  3. 03
    Среднее течение (Пуазейль)ū₀ = 0.1450 m/s · Q₀ = 4.981 L/min · u₀(r) = 2ū₀(1 − r²/R²)Постоянная часть цикла.
  4. 04
    Каждая гармоника (Womersley 1955)u_n(r) = Re{ C_n [1 − J₀(Λ_n r/R) / J₀(Λ_n)] e^(inωt) }, Λ_n = α√n · i^(3/2)Гармоник: 15; для первой α₁ = 20.01 и F = 1 − 2J₁(Λ)/(ΛJ₀(Λ)) = 0.9318∠4.19°.
  5. 05
    Касательное напряжение на стенкеτ_w = 4μū₀/R − Re Σ μ C_n (Λ_n/R) J₁(Λ_n)/J₀(Λ_n) e^(inωt)Пик 6.40 Pa; квазистационарная оценка 4μū/R дала бы 0.83 Pa.
  6. 06
    Индекс колебательного сдвигаOSI = ½·(1 − |∫τ dt| / ∫|τ| dt) = 0.437Ku et al. (1985): 0, если касательное напряжение на стенке никогда не меняет направление, и 0,5, если оно направлено назад так же долго, как и вперёд.
  7. 07
    Число РейнольдсаRe = ρ ū d / μ: 6 526 at peak 1 186 meanРешение предполагает ламинарное течение; оно не указывает, где реальный цикл становится турбулентным.

Источники: Womersley, J. Physiol. 127 (1955); Waite & Fine, Applied Biofluid Mechanics (McGraw-Hill 2007), гл. 1 и 7; De Tullio et al., J. Fluid Mech. 622 (2009); Ku et al., Arteriosclerosis 5 (1985).

ЗА ПРЕДЕЛАМИ ЖЁСТКОЙ ТРУБЫ

Упругие артерии, ветвления, стенты и течение за клапаном: мы моделируем их методами CFD и расчёта взаимодействия жидкости и конструкции.

Отправить задачу →

Другие калькуляторы