В неидеальной динамике жидкости , то течение пуазёйль , также известный как закон Хаген-Пуазейль , закон Пуазейля или уравнение Пуазейля , является физическим законом , который дает падение давления в качестве несжимаемой и ньютоновской жидкости в ламинарном потоке протекающее через длинный цилиндрическую трубу постоянного сечения. Его можно успешно применять для потока воздуха в альвеолах легких , потока через соломинку для питья или через иглу для подкожных инъекций . Он был экспериментально получен независимоПуазёйль в 1838 году [1] и Gotthilf Генрих Людвиг Хаген , [2] и опубликованы Пуазейля в 1840-41 и 1846. [1] Теоретическое обоснование закона Пуазейля дал Джордж Стокс в 1845. [3]
Предположения уравнения таковы, что жидкость несжимаемая и ньютоновская ; поток ламинарный через трубу постоянного круглого поперечного сечения , которая является , по существу , больше , чем его диаметр; и нет ускорения жидкости в трубе. Для скоростей и диаметров трубы выше порогового значения фактический поток жидкости не ламинарный, а турбулентный , что приводит к большим перепадам давления, чем рассчитывается по уравнению Хагена – Пуазейля.
Уравнение Пуазейля описывает падение давления из- за вязкости жидкости; В жидкости все еще могут возникать другие типы перепадов давления (см. Демонстрацию здесь). [4] Так , например, давление , необходимое для приведения в действие вязкой жидкости вверх против силы тяжести будет содержать как то, что при необходимости в законе Пуазейля плюс , что при необходимости в уравнении Бернулли , таким образом, что любая точка в потоке будет иметь большее давление , чем ноль ( в противном случае никакого потока бы не было).
Другой пример, когда кровь течет в более узком сужение , его скорость будет больше , чем в большем диаметре (за счетом непрерывности в объемной скорости потока ), и его давление будет ниже , чем в большем диаметре [4] ( в связи с уравнением Бернулли ). Однако вязкость крови вызовет дополнительное падение давления вдоль направления потока, которое пропорционально пройденной длине [4] (согласно закону Пуазейля). Оба эффекта влияют на фактическое падение давления.
Уравнение
В стандартных обозначениях кинетики жидкости: [5] [6] [7]
где:
- Δ p - разница давлений между двумя концами,
- L - длина трубы,
- μ - динамическая вязкость ,
- Q - объемный расход ,
- R - радиус трубы ,
- А - поперечное сечение трубы.
Уравнение не выполняется вблизи входа в трубу. [8] : 3
Уравнение не работает в пределе низкой вязкости, широкой и / или короткой трубы. Низкая вязкость или широкая труба могут привести к турбулентному потоку, что потребует использования более сложных моделей, таких как уравнение Дарси – Вайсбаха . Отношение длины к радиусу трубы должно быть больше одной сорок восьмой числа Рейнольдса, чтобы закон Хагена – Пуазейля действовал. [9] Если труба слишком короткая, уравнение Хагена – Пуазейля может привести к нефизически высоким расходам; поток ограничен принципом Бернулли , при менее ограничительных условиях,
потому что невозможно иметь давление ниже нуля (абсолютное) (не путать с манометрическим давлением ) в потоке несжимаемой жидкости.
Связь с уравнением Дарси – Вайсбаха.
Обычно поток Хагена – Пуазейля подразумевает не только соотношение для перепада давления, приведенное выше, но также полное решение для ламинарного профиля потока, который является параболическим. Однако результат для падения давления может быть расширен до турбулентного потока путем определения эффективной турбулентной вязкости в случае турбулентного потока, даже если профиль потока в турбулентном потоке, строго говоря, не является фактически параболическим. В обоих случаях, ламинарном или турбулентном, падение давления связано с напряжением на стенке, которое определяет так называемый коэффициент трения. Напряжение стенки может быть определено феноменологически с помощью уравнения Дарси – Вайсбаха в области гидравлики , учитывая соотношение для коэффициента трения через число Рейнольдса. В случае ламинарного течения для круглого поперечного сечения:
где Re - число Рейнольдса , ρ - плотность жидкости, а v - средняя скорость потока, которая составляет половину максимальной скорости потока в случае ламинарного потока. Более полезным оказывается определение числа Рейнольдса в терминах средней скорости потока, потому что эта величина остается хорошо определенной даже в случае турбулентного потока, тогда как максимальная скорость потока может не быть, или, в любом случае, может быть трудно сделать вывод. . В такой форме закон аппроксимирует коэффициент трения Дарси , то коэффициент потерь энергии (головки) , коэффициент потерь на трение или Дарси (трение) фактор Л в ламинарном потоке при очень низких скоростях в цилиндрической трубе. Теоретический вывод несколько иной формы закона был независимо сделан Видманом в 1856 г. и Нейманом и Э. Хагенбахом в 1858 г. (1859, 1860). Хагенбах был первым, кто назвал этот закон законом Пуазейля.
Закон также очень важен в гемореологии и гемодинамике , обеих областях физиологии . [10]
Позднее, в 1891 году, закон Пуазейля был распространен на турбулентный поток Л. Р. Уилберфорсом на основе работы Хагенбаха.
Вывод
Уравнение Хагена – Пуазейля может быть получено из уравнений Навье – Стокса . Ламинарный поток через трубу равномерной (круговой) поперечное сечение , известен как поток Хагена-Пуазейля. Уравнения, описывающие течение Хагена – Пуазейля, могут быть получены непосредственно из уравнений импульса Навье – Стокса в трехмерных цилиндрических координатах. сделав следующий набор предположений:
- Течение устойчивое ( ).
- Радиальная и азимутальная составляющие скорости жидкости равны нулю ( ).
- Течение осесимметричное ( ).
- Поток полностью развит (). Однако здесь это можно доказать с помощью сохранения массы и сделанных выше предположений.
Тогда угловое уравнение в уравнениях импульса и уравнение неразрывности выполняются тождественно. Уравнение радиального импульса сводится к, т.е. давление является функцией осевой координаты Только. Для краткости используйте вместо . Уравнение осевого импульса сводится к
где - динамическая вязкость жидкости. В приведенном выше уравнении левая часть является только функцией а член в правой части является только функцией , подразумевая, что оба члена должны быть одной и той же константой. Оценить эту константу несложно. Если принять длину трубы равной и обозначим перепад давления между двумя концами трубы через (высокое давление минус низкое давление), то константа просто определено таким образом, что положительный. Решение
С должен быть конечным в , . Граничное условие отсутствия проскальзывания на стенке трубы требует, чтобы в (радиус трубы), что дает Таким образом, мы имеем, наконец, следующий параболический профиль скорости :
Максимальная скорость достигается на средней линии трубы (), . Среднюю скорость можно получить интегрированием по поперечному сечению трубы ,
Легко измеряемая величина в экспериментах - это объемный расход. . Его перестановка дает уравнение Хагена – Пуазейля
Тщательный вывод, начиная с первых принципов Альтернативный метод вывода уравнения Хагена – Пуазейля, хотя и более длинный, чем прямое использование уравнений Навье – Стокса , заключается в следующем. Поток жидкости по трубе
Предположим, что жидкость имеет ламинарный поток . Ламинарный поток в круглой трубе предполагает наличие группы круглых слоев (пластин) жидкости, скорость каждого из которых определяется только их радиальным расстоянием от центра трубы. Также предположим, что центр движется быстрее всего, в то время как жидкость, касающаяся стенок трубки, неподвижна (из -за условия прилипания ).Чтобы определить движение жидкости, необходимо знать все силы, действующие на каждую пластинку:
- Сила давления, толкающая жидкость через трубку, представляет собой изменение давления, умноженное на площадь: F = - A Δ p . Эта сила направлена в направлении движения жидкости. Отрицательный знак происходит от обычного способа определения Δ p = p end - p top <0 .
- Эффекты вязкости будут вытягиваться из более быстрой пластинки непосредственно ближе к центру трубки.
- Эффекты вязкости будут перемещаться от более медленной пластинки непосредственно к стенкам трубки.
Вязкость
Когда два соприкасающихся друг с другом слоя жидкости движутся с разной скоростью, между ними возникает сила сдвига . Эта сила пропорциональна к зоне контакта А , с градиентом скорости , перпендикулярной к направлению потокаΔ v x/Δ y, и константа пропорциональности (вязкость) и определяется выражением
Здесь стоит отрицательный знак, потому что нас интересует более быстро движущаяся жидкость (вверху на рисунке), которая замедляется более медленной жидкостью (внизу на рисунке). Согласно третьему закону движения Ньютона сила, действующая на более медленную жидкость, равна силе, действующей на более быструю жидкость, и противоположна ей (без отрицательного знака). Это уравнение предполагает, что площадь контакта настолько велика, что мы можем игнорировать любые эффекты от краев и что жидкости ведут себя как ньютоновские жидкости .
Более быстрая пластинка
Предположим, что мы вычисляем силу, действующую на пластину радиуса r . Из приведенного выше уравнения нам нужно знать площадь контакта и градиент скорости . Думайте о пластине как о кольце радиуса r , толщины dr и длины Δ x . Площадь контакта между пластиной и более быстрой пластиной - это просто площадь внутренней части цилиндра: A = 2π r Δ x . Мы еще не знаем точной формы скорости жидкости в трубке, но мы знаем (из нашего предположения выше), что она зависит от радиуса. Следовательно, градиент скорости - это изменение скорости по отношению к изменению радиуса на пересечении этих двух пластин. Это пересечение находится в радиусе r . Итак, учитывая, что эта сила будет положительной по отношению к движению жидкости (но производная скорости отрицательна), окончательная форма уравнения станет
где вертикальная черта и индекс r после производной указывают на то, что ее следует брать при радиусе r .
Более медленная пластинка
Теперь давайте найдем силу сопротивления более медленной пластинки. Нам нужно рассчитать те же значения, что и для силы от более быстрой пластины. В этом случае площадь контакта будет r + dr вместо r . Кроме того, мы должны помнить, что эта сила противодействует направлению движения жидкости и поэтому будет отрицательной (и что производная скорости отрицательна).
Собираем все вместе
Чтобы найти решение для потока ламинарного слоя через трубку, нам нужно сделать одно последнее предположение. В трубе нет ускорения жидкости, и, согласно первому закону Ньютона , нет чистой силы. Если нет чистой силы, мы можем сложить все силы вместе, чтобы получить ноль.
или же
Во-первых, чтобы все происходило в одной точке, используйте первые два члена разложения градиента скорости в ряд Тейлора :
Выражение справедливо для всех пластинок. Группируем одинаковые члены и опускаем вертикальную черту, поскольку предполагается, что все производные находятся на радиусе r ,
Наконец, представьте это выражение в форме дифференциального уравнения , отбросив член, квадратичный по dr .
Приведенное выше уравнение совпадает с уравнением, полученным из уравнений Навье – Стокса, и вывод отсюда следует, как и ранее.
Запуск потока Пуазейля в трубе
Когда постоянный градиент давления применяется между двумя концами длинной трубы, поток не сразу приобретает профиль Пуазейля, скорее, он развивается во времени и достигает профиля Пуазейля в установившемся состоянии. Уравнения Навье – Стокса сводятся к
с начальными и граничными условиями,
Распределение скорости дается формулой
где - функция Бесселя первого рода нулевого порядка и положительные корни этой функции и - функция Бесселя первого рода порядка. В виде, Раствор Пуазейля восстанавливается. [11]
Течение Пуазейля в кольцевом сечении
Если - радиусы внутреннего цилиндра и - радиусы внешнего цилиндра с приложенным градиентом давления между двумя концами. , распределение скорости и объемный поток через кольцевую трубу равны
Когда , исходная проблема устранена. [12]
Течение Пуазейля в трубе с колеблющимся градиентом давления
Поток по трубам с колеблющимся градиентом давления находит применение в кровотоке через крупные артерии. [13] [14] [15] [16] Накладываемый градиент давления определяется как
где , а также константы и это частота. Поле скорости определяется выражением
где
где а также - функции Кельвина и.
Плоский поток Пуазейля
Плоский поток Пуазейля - это поток, создаваемый между двумя бесконечно длинными параллельными пластинами, разделенными расстоянием. с постоянным градиентом давления наносится по направлению потока. Поток по существу однонаправленный из-за бесконечной длины. Уравнения Навье – Стокса сводятся к
с противоскользящим покрытием на обеих стенах
Следовательно, распределение скорости и объемный расход на единицу длины равны
Поток Пуазейля через некруглые поперечные сечения
Джозеф Буссинеск вывел профиль скорости и объемный расход в 1868 году для прямоугольного канала и труб равностороннего треугольного сечения и для эллиптического сечения. [17] Джозеф Праудман вывел то же самое для равнобедренных треугольников в 1914 году. [18] Пусть - постоянный градиент давления, действующий в направлении, параллельном движению.
Скорость и объемный расход в прямоугольном канале высотой и ширина находятся
Скорость и объемный расход трубки с равносторонним треугольным поперечным сечением длины стороны находятся
Скорость и объемный расход в прямоугольном равнобедренном треугольнике находятся
Распределение скорости для труб эллиптического сечения с полуосью а также это [11]
Здесь, когда , Поток Пуазейля для круглой трубы восстанавливается, а когда , восстанавливается плоское течение Пуазейля . Более явные решения с поперечными сечениями, такими как секции в форме улитки, секции, имеющие форму круга с надрезом, следующего за полукругом, кольцевые секции между гомофокальными эллипсами, кольцевые секции между неконцентрическими окружностями, также доступны, как описано Ratip Berker . [19] [20]
Течение Пуазейля через произвольное сечение
Поток через произвольное сечение удовлетворяет условию, что на стенах. Основное уравнение сводится к [21]
Если мы введем новую зависимую переменную как
то легко видеть, что проблема сводится к интегрированию уравнения Лапласа
удовлетворяющий условию
на стене.
Уравнение Пуазейля для идеального изотермического газа
Для сжимаемой жидкости в трубке объемный расход (но не массовый расход) и осевая скорость непостоянны вдоль трубы. Расход обычно выражается при давлении на выходе. Когда жидкость сжимается или расширяется, работа выполняется, и жидкость нагревается или охлаждается. Это означает, что скорость потока зависит от теплопередачи к жидкости и от жидкости. Для идеального газа в изотермическом случае, когда температуре жидкости позволяют уравновеситься с окружающей средой, можно получить приблизительное соотношение для перепада давления. [22] Используя уравнение состояния идеального газа для процесса с постоянной температурой, соотношениеможет быть получен. На коротком участке трубы газ, текущий по трубе, можно считать несжимаемым, так что закон Пуазейля можно использовать локально,
Здесь мы предположили, что локальный градиент давления не слишком велик, чтобы иметь какие-либо эффекты сжимаемости. Хотя локально мы игнорировали эффекты изменения давления из-за изменения плотности, на больших расстояниях эти эффекты учитываются. С не зависит от давления, приведенное выше уравнение можно проинтегрировать по длине дать
Следовательно, объемный расход на выходе из трубы определяется выражением
Это уравнение можно рассматривать как закон Пуазейля с дополнительным поправочным коэффициентом. п 1 + п 2/2 п 2 выражая среднее давление относительно давления на выходе.
Аналогия электрических схем
Первоначально считалось, что электричество - это своего рода жидкость. Эта гидравлическая аналогия до сих пор концептуально полезна для понимания схем. Эта аналогия также используется для изучения частотной характеристики жидкостно-механических сетей с использованием схемотехнических инструментов, и в этом случае жидкостная сеть называется гидравлической схемой . Закон Пуазейля соответствует закону Ома для электрических цепей, V = IR . Поскольку результирующая сила, действующая на жидкость, равна, где S = π r 2 , т.е. Δ F = π r 2 Δ P , то из закона Пуазейля следует, что
- .
Для электрических цепей пусть n будет концентрацией свободных заряженных частиц (в м −3 ) и пусть q * будет зарядом каждой частицы (в кулонах ). (Для электронов q * = e =1,6 × 10 −19 Кл .) Тогда nQ - количество частиц в объеме Q , а nQq * - их полный заряд. Это зарядкоторый проходит через поперечное сечение в единицу времени, то есть ток I . Следовательно, I = nQq * . Следовательно, Q = я/nq *, а также
Но Δ F = Eq , где q - полный заряд в объеме трубки. Объем трубки равен π r 2 L , поэтому количество заряженных частиц в этом объеме равно n π r 2 L , а их общий заряд равенПоскольку напряжение V = EL , то
Это в точности закон Ома, где сопротивление R = V/я описывается формулой
- .
Отсюда следует, что сопротивление R пропорционально длине резистора L , что верно. Однако также следует, что сопротивление R обратно пропорционально четвертой степени радиуса r , то есть сопротивление R обратно пропорционально второй степени площади поперечного сечения S = π r 2 резистора, которая отличается от электрическая формула. Электрическое соотношение для сопротивления:
где ρ - удельное сопротивление; т.е. сопротивление R обратно пропорционально площади поперечного сечения S резистора. [23] Причина, по которой закон Пуазейля приводит к другой формуле для сопротивления R, заключается в разнице между потоком жидкости и электрическим током. Электронный газ является невязким , так что его скорость не зависит от расстояния до стенки проводника. Сопротивление возникает из-за взаимодействия между текущими электронами и атомами проводника. Следовательно, закон Пуазейля и гидравлическая аналогия применимы к электричеству только в определенных пределах. И закон Ома, и закон Пуазейля иллюстрируют явления переноса .
Медицинские применения - внутривенный доступ и доставка жидкости
Уравнение Хагена – Пуазейля полезно для определения сосудистого сопротивления и, следовательно, скорости потока внутривенных (IV) жидкостей, что может быть достигнуто с использованием периферических и центральных канюль различных размеров . Уравнение утверждает, что скорость потока пропорциональна радиусу в четвертой степени, а это означает, что небольшое увеличение внутреннего диаметра канюли приводит к значительному увеличению скорости потока жидкостей для внутривенного вливания. Радиус канюли для внутривенных вливаний обычно измеряется в "калибровке", которая обратно пропорциональна [ сомнительно ] радиусу. Периферические канюли для внутривенных вливаний обычно выпускаются (от больших до малых) 14G, 16G, 18G, 20G, 22G, 26G. Например, поток канюли 14G обычно примерно в два раза [ необходима цитата ] , чем у канюли 16G, и в десять раз [ требуется ссылка ] , чем у канюли 20G. В нем также указано, что поток обратно пропорционален длине, а это означает, что более длинные линии имеют более низкие скорости потока. Это важно помнить, поскольку в экстренных случаях многие врачи предпочитают более короткие и большие катетеры по сравнению с более длинными и узкими катетерами. В то время как менее клинического значения, увеличенное изменение давления (Δ р ) - например, путь повышения давления в мешке жидкости, сжимая сумку, или висят мешок выше ( по отношению к уровню канюли) - может быть использовано для ускорения скорость потока. Также полезно понимать, что вязкие жидкости будут течь медленнее (например, при переливании крови ).
Смотрите также
- Поток Куэтта
- Закон Дарси
- Пульс
- Волна
- Гидравлический контур
Заметки
Цитированные ссылки
- ^ a b Sutera, Salvatore P .; Скалак, Ричард (1993). «История закона Пуазейля» . Ежегодный обзор гидромеханики . 25 : 1–19. Bibcode : 1993AnRFM..25 .... 1S . DOI : 10.1146 / annurev.fl.25.010193.000245 .
- ↑ István Szabó, ;; Geschichte der Mechanischen Prinzipien und ihrer wichtigsten Anwendungen , Базель: Birkhäuser Verlag, 1979.
- ^ Стокс, GG (1845). О теориях внутреннего трения движущихся жидкостей, равновесия и движения упругих тел. Труды Кембриджского философского общества, 8, 287–341.
- ^ а б в «Давление» . hyperphysics.phy-astr.gsu.edu . Проверено 15 декабря 2019 .
- ^ Кирби, Би Джей (2010). Микро- и наномасштабная механика жидкости: транспорт в микрофлюидных устройствах . Издательство Кембриджского университета. ISBN 978-0-521-11903-0.
- ^ Брус, Х. (2007). Теоретическая микрофлюидика .
- ^ Пфицнер, Дж. (1976). «Пуазей и его закон» (PDF) . Анестезия . 31 (2): 273–275. DOI : 10.1111 / j.1365-2044.1976.tb11804.x . PMID 779509 . S2CID 40607063 . Архивировано из оригинального (PDF) 10 августа 2017 года.
- ^ Фогель, Стивен (1981). Жизнь в движущихся жидкостях: физическая биология потока . PWS Kent Publishers. ISBN 0871507498.
- ^ tec-science (02.04.2020). «Энергетический анализ закона Хагена – Пуазейля» . тек-наука . Проверено 7 мая 2020 .
- ^ Детерминанты сопротивления кровеносных сосудов .
- ^ a b Бэтчелор, Джордж Кейт. Введение в гидродинамику . Издательство Кембриджского университета, 2000.
- ^ Розенхед, Луи, изд. Ламинарные пограничные слои. Кларендон Пресс, 1963.
- ^ Сексл, Т. (1930). "Über den von EG Richardson entdeckten 'Annulareffekt ' ". Zeitschrift für Physik . 61 (5–6): 349–362. DOI : 10.1007 / BF01340631 .
- ^ Ламбосси П. (1952). «Колебания заставляют несжимаемую жидкость и вязкость в жесткой и горизонтальной трубке. Расчет силы взрыва». Helv. Physica Acta . 25 : 371–386.
- ^ Уомерсли, младший (1955). «Метод расчета скорости, скорости потока и вязкого сопротивления в артериях, когда известен градиент давления» . Журнал физиологии . 127 (3): 553–563. DOI : 10.1113 / jphysiol.1955.sp005276 . PMC 1365740 . PMID 14368548 .
- ^ Учида, С. (1956). «Пульсирующий вязкий поток, наложенный на установившееся ламинарное движение несжимаемой жидкости в круглой трубе». Zeitschrift für angewandte Mathematik und Physik . 7 (5): 403–422. DOI : 10.1007 / BF01606327 .
- ^ Буссинеск, Джозеф (1868). "Воспоминания о влиянии жидкостей на регулирующие движения жидкостей". J. Math. Pures Appl . 13 (2): 377–424.
- ^ Праудмен, Дж. (1914). «Заметки о движении вязкой жидкости в каналах» . Лондонский, Эдинбургский и Дублинский философский журнал и научный журнал . 28 (163): 30–36. DOI : 10.1080 / 14786440708635179 .
- ^ Беркер, Р. (1963). "Интеграция условий движения несжимаемой жидкости". Handbuch der Physik . 3 . С. 1–384.
- ^ Дразин, Филип Г .; Райли, Норман (2006). Уравнения Навье – Стокса: классификация потоков и точные решения . № 334. Издательство Кембриджского университета. ISBN 9780521681629.
- ^ Кёрл, Сэмюэл Ньюби ; Дэвис, HJ (1971). Современная гидродинамика . Том 1, Несжимаемый поток. Ван Ностранд Рейнхольд.
|volume=
есть дополнительный текст ( справка ) - ^ Ландау, ЛД ; Лифшиц, Э.М. (1987). Механика жидкости . Pergamon Press. п. 55, проблема 6. ISBN 0-08-033933-6.
- ^ Fütterer, C .; и другие. (2004). «Система впрыска и регулирования потока для микроканалов». Лаборатория на чипе . 4 (4): 351–356. DOI : 10.1039 / B316729A . PMID 15269803 .
Рекомендации
- Сутера, ИП; Скалак Р. (1993). «История закона Пуазейля» . Ежегодный обзор гидромеханики . 25 : 1–19. Bibcode : 1993AnRFM..25 .... 1S . DOI : 10.1146 / annurev.fl.25.010193.000245 ..
- Пфицнер, Дж (1976). «Пуазей и его закон». Анестезия . 31 (2) (опубликовано в марте 1976 г.). С. 273–5. DOI : 10.1111 / j.1365-2044.1976.tb11804.x . PMID 779509 ..
- Беннетт, Колорадо; Майерс, Дж. Э. (1962). Импульс, тепло и массообмен . Макгроу-Хилл..
Внешние ссылки
- Закон Пуазейля для степенной неньютоновской жидкости
- Закон Пуазейля в слегка сужающейся трубке
- Калькулятор уравнения Хагена – Пуазейля