stringtranslate.com

Вязкость

Вязкость жидкости является мерой ее сопротивления деформации при заданной скорости . [1] Для жидкостей это соответствует неформальному понятию «густота»: например, сироп имеет более высокую вязкость, чем вода . [2] Вязкость определяется с научной точки зрения как сила, умноженная на время, разделенное на площадь. Таким образом, единицей измерения в системе СИ являются ньютоны-секунды на квадратный метр или паскаль-секунды. [1]

Вязкость количественно определяет силу внутреннего трения между соседними слоями жидкости, находящимися в относительном движении. [1] Например, когда вязкая жидкость проталкивается через трубку, она течет быстрее вблизи оси трубки, чем вблизи ее стенок. [3] Эксперименты показывают, что для поддержания потока необходимо некоторое напряжение (например, разница давлений между двумя концами трубки). Это связано с тем, что для преодоления трения между слоями жидкости, находящимися в относительном движении, требуется сила. Для трубки с постоянной скоростью потока сила компенсирующей силы пропорциональна вязкости жидкости.

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

Нулевая вязкость (отсутствие сопротивления сдвиговому напряжению ) наблюдается только при очень низких температурах в сверхтекучих средах ; в противном случае второй закон термодинамики требует, чтобы все жидкости имели положительную вязкость. [4] [5] Жидкость, имеющая нулевую вязкость (невязкая), называется идеальной или невязкой .

Этимология

Слово «вязкость» происходит от латинского viscumомела »). Вискум также назывался вязким клеем , полученным из ягод омелы. [6]

Определения

Динамическая вязкость

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

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

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

В течении Куэтта жидкость зажата между двумя бесконечно большими пластинами: одна неподвижна, а другая движется параллельно с постоянной скоростью (см. иллюстрацию справа). Если скорость верхней пластины достаточно мала (чтобы избежать турбулентности), то в установившемся режиме частицы жидкости движутся параллельно ей, а их скорость меняется от низа к верху. [7] Каждый слой жидкости движется быстрее, чем тот, что находится непосредственно под ним, и трение между ними порождает силу, сопротивляющуюся их относительному движению. В частности, жидкость прикладывает к верхней пластине силу в направлении, противоположном ее движению, и равную, но противоположную силу к нижней пластине. Поэтому требуется внешняя сила, чтобы поддерживать движение верхней пластины с постоянной скоростью.

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

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

давление , умноженное на время , энергия на единицу объема, умноженная на время.

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

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

Использование греческой буквы мю ( ) для обозначения динамической вязкости (иногда также называемой абсолютной вязкостью ) распространено среди инженеров -механиков и химиков , а также математиков и физиков. [8] [9] [10] Однако греческая буква эта ( ) также используется химиками, физиками и ИЮПАК . [11] Вязкость иногда также называют сдвиговой вязкостью . Однако, по крайней мере, один автор не рекомендует использовать эту терминологию, отмечая, что она может встречаться не только в сдвиговых потоках, но и в несдвиговых потоках. [12]

Кинематическая вязкость

В гидродинамике иногда более целесообразно работать с точки зрения кинематической вязкости (иногда также называемой коэффициентом диффузии импульса ), определяемой как отношение динамической вязкости ( μ ) к плотности жидкости ( ρ ). Обычно его обозначают греческой буквой nu ( ν ):

и имеет размеры , что приводит к единицам СИ и производным единицам :

удельная энергия , умноженная на время , энергия на единицу массы, умноженная на время.

Общее определение

В самых общих чертах вязкие напряжения в жидкости определяются как напряжения, возникающие в результате относительной скорости различных частиц жидкости. Таким образом, вязкие напряжения должны зависеть от пространственных градиентов скорости потока. Если градиенты скорости малы, то вязкие напряжения в первом приближении зависят только от первых производных скорости. [13] (Для ньютоновских жидкостей это также линейная зависимость.) В декартовых координатах общее соотношение тогда можно записать как

где – тензор вязкости, который отображает тензор градиента скорости на тензор вязкого напряжения . [14] Поскольку индексы в этом выражении могут варьироваться от 1 до 3, всего существует 81 «коэффициент вязкости». Однако если предположить, что тензор ранга вязкости 4 изотропен , эти 81 коэффициент сводятся к трем независимым параметрам , , :

и, кроме того, предполагается, что никакие силы вязкости не могут возникнуть, когда жидкость подвергается простому вращению твердого тела, таким образом , остаются только два независимых параметра. [13] Наиболее обычное разложение происходит через стандартную (скалярную) вязкость и объемную вязкость, такую ​​что и . В векторной записи это выглядит так:

где – единичный тензор, а крестик обозначает транспонирование . [12] [15] Это уравнение можно рассматривать как обобщенную форму закона вязкости Ньютона.

Объемная вязкость (также называемая объемной вязкостью) отражает тип внутреннего трения, которое препятствует сжатию или расширению жидкости без сдвига. Знание часто не требуется в задачах гидродинамики. Например, несжимаемая жидкость удовлетворяет условиям , поэтому член «содержащий» выпадает. Более того, для газов часто предполагается, что он пренебрежимо мал, поскольку он находится в одноатомном идеальном газе . [12] Одной из ситуаций, в которой может быть важно, является расчет потерь энергии в звуковых и ударных волнах , описываемых законом затухания звука Стокса , поскольку эти явления связаны с быстрым расширением и сжатием.

Определяющие уравнения вязкости не являются фундаментальными законами природы, поэтому их полезность, а также методы измерения или расчета вязкости должны устанавливаться отдельными средствами. Потенциальная проблема заключается в том, что вязкость, в принципе, зависит от полного микроскопического состояния жидкости, которое включает в себя положение и импульс каждой частицы в системе. [16] Такая подробная информация обычно недоступна в реалистичных системах. Однако при определенных условиях большая часть этой информации может оказаться незначительной. В частности, для ньютоновских жидкостей вблизи равновесия и вдали от границ (объемное состояние) вязкость зависит только от пространственно-временных макроскопических полей (таких как температура и плотность), определяющих локальное равновесие. [16] [17]

Тем не менее, вязкость все еще может иметь существенную зависимость от нескольких свойств системы, таких как температура, давление, а также амплитуда и частота любого внешнего воздействия. Таким образом, точные измерения вязкости определяются только в отношении конкретного состояния жидкости. [18] Чтобы стандартизировать сравнения экспериментов и теоретических моделей, данные о вязкости иногда экстраполируются на идеальные предельные случаи, такие как предел нулевого сдвига или (для газов) предел нулевой плотности .

Импульсный транспорт

Теория переноса предлагает альтернативную интерпретацию вязкости с точки зрения переноса импульса: вязкость — это свойство материала, которое характеризует перенос импульса внутри жидкости, точно так же, как теплопроводность характеризует перенос тепла , а коэффициент диффузии (массовой) характеризует перенос массы. [19] Эта точка зрения неявно заложена в законе вязкости Ньютона , потому что напряжение сдвига имеет единицы, эквивалентные потоку импульса , т.е. импульс в единицу времени на единицу площади. Таким образом, его можно интерпретировать как задание потока импульса в направлении от одного слоя жидкости к другому. Согласно закону вязкости Ньютона, этот поток импульса происходит поперек градиента скорости, и величина соответствующего потока импульса определяется вязкостью.

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

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

Ньютоновские и неньютоновские жидкости

Вязкость, наклон каждой линии, варьируется в зависимости от материала.

Закон вязкости Ньютона не является фундаментальным законом природы, а скорее определяющим уравнением (подобно закону Гука , закону Фика и закону Ома ), которое служит для определения вязкости . Его форма обусловлена ​​экспериментами, которые показывают, что для широкого диапазона жидкостей он не зависит от скорости деформации. Такие жидкости называются ньютоновскими . Газы , воду и многие обычные жидкости можно считать ньютоновскими в обычных условиях и контекстах. Однако существует множество неньютоновских жидкостей , поведение которых значительно отклоняется от такого поведения. Например:

Коэффициент Трутона представляет собой отношение вязкости при растяжении к вязкости при сдвиге . Для ньютоновской жидкости коэффициент Трутона равен 3. [21] [22] Жидкости, разжижающиеся при сдвиге, очень часто, но ошибочно, называют тиксотропными. [23]

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

В твердых телах

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

Однако многие жидкости (включая воду) кратковременно реагируют как упругие твердые тела при внезапном напряжении. И наоборот, многие «твердые тела» (даже гранит ) будут течь как жидкости, хотя и очень медленно, даже при сколь угодно малом напряжении. [24] Такие материалы лучше всего охарактеризовать как вязкоупругие , то есть обладающие как эластичностью (реакция на деформацию), так и вязкостью (реакция на скорость деформации).

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

В геологии земные материалы, демонстрирующие вязкую деформацию, по крайней мере на три порядка превышающую их упругую деформацию, иногда называют рейдами . [25]

Измерение

Вязкость измеряют различными типами вискозиметров и реометров . Реометр используется для жидкостей, вязкость которых невозможно определить по одному значению и поэтому требуется установка и измерение большего количества параметров, чем в случае вискозиметра. Тщательный контроль температуры жидкости необходим для получения точных измерений, особенно для таких материалов, как смазочные материалы, вязкость которых может удвоиться при изменении всего на 5 °C. [26]

Для некоторых жидкостей вязкость постоянна в широком диапазоне скоростей сдвига ( ньютоновские жидкости ). Жидкости без постоянной вязкости ( неньютоновские жидкости ) не могут быть описаны одним числом. Неньютоновские жидкости демонстрируют множество различных корреляций между напряжением сдвига и скоростью сдвига.

Одним из наиболее распространенных приборов для измерения кинематической вязкости является стеклянный капиллярный вискозиметр.

В промышленности по нанесению покрытий вязкость можно измерять с помощью чашки, в которой измеряется время истечения . Существует несколько типов чашек, таких как чашка Зана и чашка вязкости Форда , причем использование каждого типа варьируется в основном в зависимости от отрасли.

Вискозиметр Штормера , также используемый в покрытиях, использует вращение в зависимости от нагрузки для определения вязкости. Вязкость указывается в единицах Кребса (КУ), которые являются уникальными для вискозиметров Штормера.

Вибрационные вискозиметры также можно использовать для измерения вязкости. Резонансные или вибрационные вискозиметры работают за счет создания поперечных волн в жидкости. В этом методе датчик погружается в жидкость и заставляет его резонировать на определенной частоте. Когда поверхность датчика проходит через жидкость, энергия теряется из-за ее вязкости. Затем эта рассеиваемая энергия измеряется и преобразуется в показания вязкости. Более высокая вязкость приводит к большим потерям энергии. [ нужна цитата ]

Вязкость при растяжении можно измерить с помощью различных реометров , создающих напряжение при растяжении .

Объемную вязкость можно измерить акустическим реометром .

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

Нановязкость (вязкость, измеряемая нанозондами) можно измерить с помощью флуоресцентной корреляционной спектроскопии . [27]

Единицы

Единицей динамической вязкости в системе СИ является ньютон -секунда на квадратный метр (Н·с/м 2 ), также часто выражаемая в эквивалентных формах паскаль - секунда (Па·с), килограмм на метр в секунду (кг·м -1) . ·с -1 ) и Пуазейля (Pl). Единицей СГС является пуаз (П, или г·см -1 ·с -1 = 0,1 Па·с), [28] названный в честь Жана Леонара Мари Пуазейля . Обычно его выражают, особенно в стандартах ASTM , в сантипуазах (сП). Сантипуаз удобен тем, что вязкость воды при 20 °C составляет около 1 сП, а один сантипуаз равен миллипаскалю-секунде системы СИ (мПа·с).

Единицей кинематической вязкости в системе СИ является квадратный метр в секунду (м 2 /с), тогда как единицей кинематической вязкости в СГС является стокс ( Ст, или см 2 ·с -1 = 0,0001 м 2 ·с -1 ), названный в честь Сэр Джордж Гэбриел Стоукс . [29] В США слово «стоук» иногда употребляется в единственном числе. Вместо этого часто используют дробные сантистоксы (сСт), 1 сСт = 1 мм 2 · с -1 = 10  -6 м  2 · с -1 . 1 сСт равен 1 сП, разделенному на 1000 кг/м^3, что близко к плотности воды. Кинематическая вязкость воды при 20 °С составляет около 1 сСт.

Наиболее часто используемые системы обычных или имперских единиц США — это британская гравитационная (BG) и английская инженерная (EE). В системе BG динамическая вязкость измеряется в фунтах -секундах на квадратный фут (фунт·с/фут 2 ), а в системе EE – в фунтах-силах -секундах на квадратный фут (фунт-сила-секунда на квадратный фут (фунт-сила-секунда на квадратный фут 2 ). . Фунт и фунт-сила эквивалентны; две системы различаются только тем, как определяются сила и масса. В системе BG фунт является основной единицей, из которой единица массы ( снаряд ) определяется вторым законом Ньютона , тогда как в системе EE единицами силы и массы (фунт-сила и фунт-масса соответственно) являются определяется независимо посредством Второго закона с использованием константы пропорциональности g c .

Кинематическая вязкость измеряется в квадратных футах в секунду (фут 2 /с) как в системах BG, так и в системах EE.

К нестандартным единицам измерения относится рейн , британская единица динамической вязкости. [30] В автомобильной промышленности индекс вязкости используется для описания изменения вязкости с температурой.

Обратная величина вязкости — это текучесть , обычно обозначаемая или , в зависимости от используемого соглашения, измеряемая в обратном пуазе (P -1 или см · с · г -1 ), иногда называемом rhe . Текучесть редко используется в инженерной практике. [ нужна цитата ]

Одно время нефтяная промышленность полагалась на измерение кинематической вязкости с помощью вискозиметра Сейболта и выражение кинематической вязкости в единицах универсальных секунд Сейболта (SUS). [31] Иногда используются другие сокращения, такие как SSU ( универсальные секунды Saybolt ) или SUV ( универсальная вязкость Saybolt ). Кинематическая вязкость в сантистоксах может быть преобразована из SUS в соответствии с арифметикой и справочной таблицей, приведенной в ASTM D 2161.

Молекулярное происхождение

Перенос импульса в газах осуществляется дискретными столкновениями молекул, а в жидкостях — силами притяжения, которые связывают молекулы близко друг к другу. [19] Из-за этого динамическая вязкость жидкостей обычно намного выше, чем у газов. Кроме того, вязкость имеет тенденцию увеличиваться с температурой в газах и уменьшаться с температурой в жидкостях.

Выше критической точки жидкость-газ жидкая и газовая фазы заменяются одной сверхкритической фазой . В этом режиме механизмы переноса импульса интерполируют поведение жидкости и газа. Например, вдоль сверхкритической изобары (поверхности постоянного давления) кинематическая вязкость уменьшается при низкой температуре и увеличивается при высокой температуре с минимумом между ними. [32] [33] Грубая оценка минимального значения равна

где – постоянная Планка , – масса электрона , – молекулярная масса. [33]

Однако вообще вязкость системы в деталях зависит от того, как взаимодействуют молекулы, составляющие систему, и для нее не существует простых, но правильных формул. Простейшими точными выражениями являются соотношения Грина-Кубо для линейной сдвиговой вязкости или выражения временной корреляционной функции, полученные Эвансом и Морриссом в 1988 году. [34] Хотя каждое из этих выражений является точным, расчет вязкости плотной жидкости с использованием этих соотношений в настоящее время требует использования компьютерного моделирования молекулярной динамики . Несколько большего прогресса можно добиться для разбавленного газа, поскольку элементарные предположения о том, как молекулы газа движутся и взаимодействуют, приводят к базовому пониманию молекулярной природы вязкости. Более сложные методы лечения могут быть построены путем систематического уточнения уравнений движения молекул газа. Примером такой трактовки является теория Чепмена-Энскога , которая выводит выражения для вязкости разбавленного газа из уравнения Больцмана . [17]

Чистые газы

Вязкость газов возникает главным образом из-за молекулярной диффузии , которая передает импульс между слоями потока. Элементарный расчет для разбавленного газа по температуре и плотности дает

где – постоянная Больцмана , молекулярная масса и числовая константа порядка . Величина , средняя длина свободного пробега , измеряет среднее расстояние, которое молекула проходит между столкновениями. Даже без априорного знания , это выражение имеет нетривиальные последствия. В частности, поскольку он обычно обратно пропорционален плотности и увеличивается с температурой, он сам должен увеличиваться с температурой и не зависеть от плотности при фиксированной температуре. Фактически, оба этих предсказания сохраняются при более сложных методах лечения и точно описывают экспериментальные наблюдения. Напротив, вязкость жидкости обычно уменьшается с температурой. [19] [35]

Для жестких упругих сфер диаметр можно вычислить, дав

В этом случае не зависит от температуры, т. е. . Однако для более сложных молекулярных моделей она зависит от температуры нетривиальным образом, и простые кинетические аргументы, используемые здесь, недостаточны. Более фундаментально, понятие средней длины свободного пробега становится неточным для частиц, которые взаимодействуют в конечном диапазоне, что ограничивает полезность концепции для описания реальных газов. [36]

Теория Чепмена – Энскога

Методика, разработанная Сиднеем Чепменом и Дэвидом Энскогом в начале 1900-х годов, позволяет более точно рассчитать . [17] Оно основано на уравнении Больцмана , которое обеспечивает статистическое описание разбавленного газа с точки зрения межмолекулярных взаимодействий. [37] Этот метод позволяет точно рассчитать молекулярные модели, которые более реалистичны, чем жесткие упругие сферы, например, модели, включающие межмолекулярное притяжение. Это необходимо для воспроизведения правильной температурной зависимости , которая, как показывают эксперименты, увеличивается быстрее, чем тенденция, предсказанная для жестких упругих сфер. [19] Действительно, анализ Чепмена-Энскога показывает, что предсказанную температурную зависимость можно настроить, варьируя параметры в различных молекулярных моделях. Простым примером является модель Сазерленда [а] , которая описывает твердые упругие сферы со слабым взаимным притяжением. В таком случае силу притяжения можно рассматривать пертурбативно , что приводит к простому выражению для :

где не зависит от температуры и определяется только параметрами межмолекулярного притяжения. Для связи с экспериментом удобно переписать в виде

где вязкость при температуре . [38] Если известно из экспериментов при и хотя бы еще одной температуре, то ее можно рассчитать. Полученные таким образом выражения для ряда простых газов качественно точны. Немного более сложные модели, такие как потенциал Леннарда-Джонса или более гибкий потенциал Ми , могут обеспечить лучшее согласие с экспериментами, но только за счет более непрозрачной зависимости от температуры. Еще одним преимуществом этих более сложных потенциалов взаимодействия является то, что их можно использовать для разработки точных моделей для широкого спектра свойств с использованием одних и тех же потенциальных параметров. В ситуациях, когда доступно мало экспериментальных данных, это позволяет получить параметры модели путем подгонки к таким свойствам, как равновесие пар-жидкость чистой жидкости , прежде чем использовать полученные таким образом параметры для прогнозирования интересующей вязкости с разумной точностью.

В некоторых системах необходимо отказаться от предположения о сферической симметрии , как это имеет место для паров с высокополярными молекулами типа H 2 O. В этих случаях анализ Чепмена–Энскога значительно сложнее. [39] [40]

Объемная вязкость

В кинетически-молекулярной картине ненулевая объемная вязкость возникает в газах всякий раз, когда существуют существенные релаксационные временные масштабы, управляющие обменом энергией между поступательной энергией молекул и их внутренней энергией, например вращательной и колебательной . Таким образом, объемная вязкость соответствует одноатомному идеальному газу, в котором внутренняя энергия молекул незначительна, но отлична от нуля для такого газа, как углекислый газ , молекулы которого обладают как вращательной, так и колебательной энергией. [41] [42]

Чистые жидкости

Видео, показывающее три жидкости с разной вязкостью
Эксперимент, показывающий поведение вязкой жидкости с синей краской для наглядности.

В отличие от газов, в жидкостях не существует простой, но точной картины молекулярного происхождения вязкости.

На простейшем уровне описания относительному движению соседних слоев в жидкости противодействуют прежде всего молекулярные силы притяжения, действующие через границу слоя. На этой картине можно (и правильно) ожидать, что вязкость уменьшится с повышением температуры. Это связано с тем, что повышение температуры увеличивает случайное тепловое движение молекул, что облегчает им преодоление притягивающих взаимодействий. [43]

Основываясь на этой визуализации, можно построить простую теорию по аналогии с дискретной структурой твердого тела: группы молекул в жидкости визуализируются как образующие «клетки», которые окружают и заключают в себе отдельные молекулы. [44] Эти клетки могут быть заняты или незаняты, и более сильное молекулярное притяжение соответствует более сильным клеткам. Из-за случайного теплового движения молекула «прыгает» между клетками со скоростью, которая обратно пропорциональна силе молекулярного притяжения. В равновесии эти «прыжки» не смещены ни в каком направлении. С другой стороны, чтобы два соседних слоя могли двигаться относительно друг друга, «прыжки» должны быть смещены в направлении относительного движения. Силу, необходимую для поддержания этого направленного движения, можно оценить для заданной скорости сдвига, что приводит к

где – постоянная Авогадро , – постоянная Планка , – объем моля жидкости , – нормальная температура кипения . Этот результат имеет тот же вид, что и известное эмпирическое соотношение

где и – константы, соответствующие данным. [44] [45] С другой стороны, некоторые авторы выражают осторожность в отношении этой модели. При использовании уравнения ( 1 ) могут возникнуть ошибки до 30% по сравнению с уравнением подгонки ( 2 ) к экспериментальным данным. [44] Более фундаментально, физические предположения, лежащие в основе уравнения ( 1 ), подверглись критике. [46] Также утверждалось, что экспоненциальная зависимость в уравнении ( 1 ) не обязательно описывает экспериментальные наблюдения более точно, чем более простые, неэкспоненциальные выражения. [47] [48]

В свете этих недостатков разработка менее специальной модели представляет практический интерес. Отказ от простоты в пользу точности позволяет написать строгие выражения для вязкости, исходя из фундаментальных уравнений движения молекул. Классическим примером такого подхода является теория Ирвинга–Кирквуда. [49] С другой стороны, такие выражения задаются как средние значения по многочастичным корреляционным функциям и поэтому их трудно применять на практике.

В целом, выражения, полученные эмпирическим путем (основанные на существующих измерениях вязкости), кажутся единственным неизменно надежным средством расчета вязкости жидкостей. [50]

Смеси и смеси

Газовые смеси

Ту же молекулярно-кинетическую картину однокомпонентного газа можно применить и к газовой смеси. Например, в подходе Чепмена-Энскога вязкость бинарной смеси газов может быть записана через вязкости отдельных компонентов , их соответствующие объемные доли и межмолекулярные взаимодействия. [17]

Что касается однокомпонентного газа, то зависимость от параметров межмолекулярных взаимодействий проявляется через различные столкновительные интегралы, которые не могут быть выражены в замкнутой форме . Чтобы получить полезные выражения, которые разумно соответствуют экспериментальным данным, интегралы столкновений могут быть рассчитаны численно или на основе корреляций. [51] В некоторых случаях интегралы столкновений рассматриваются как подходящие параметры и подгоняются непосредственно к экспериментальным данным. [52] Это распространенный подход при разработке справочных уравнений для вязкости газовой фазы. Примером такой процедуры является подход Сазерленда для однокомпонентного газа, обсуждавшийся выше.

Было показано , что для газовых смесей, состоящих из простых молекул, пересмотренная теория Энскога точно отражает зависимость вязкости как от плотности, так и от температуры в широком диапазоне условий. [53] [51]

Смеси жидкостей

Что касается чистых жидкостей, вязкость смеси жидкостей трудно предсказать с помощью молекулярных принципов. Один из методов — расширить представленную выше теорию молекулярной «клетки» на случай чистой жидкости. Это можно сделать с разной степенью сложности. Одним из выражений, полученных в результате такого анализа, является уравнение Ледерера – Регирса для бинарной смеси:

где – эмпирический параметр, и – соответствующие мольные доли и вязкости составляющих жидкостей. [54]

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

Растворы и суспензии

Водные растворы

В зависимости от растворенного вещества и диапазона концентраций водный раствор электролита может иметь большую или меньшую вязкость по сравнению с чистой водой при той же температуре и давлении. Например, 20%-ный солевой раствор ( хлорид натрия ) имеет вязкость более чем в 1,5 раза больше, чем чистая вода, тогда как 20%-ный раствор йодида калия имеет вязкость примерно в 0,91 раза больше, чем чистая вода.

Идеализированная модель разбавленных электролитических растворов приводит к следующему предсказанию вязкости раствора: [55]

где – вязкость растворителя, – концентрация, – положительная константа, которая зависит как от свойств растворителя, так и от свойств растворенного вещества. Однако это выражение справедливо только для очень разбавленных растворов, имеющих концентрацию менее 0,1 моль/л. [56] Для более высоких концентраций необходимы дополнительные члены, которые учитывают молекулярные корреляции более высокого порядка:

где и соответствуют данным. В частности, отрицательное значение способно объяснить снижение вязкости, наблюдаемое в некоторых растворах. Ниже приведены расчетные значения этих констант для хлорида натрия и йодида калия при температуре 25 °C (моль = моль , L = литр ). [55]

Подвески

В суспензии твердых частиц (например, сфер микронного размера, взвешенных в масле) эффективная вязкость может быть определена через компоненты напряжения и деформации, которые усреднены по объему, большому по сравнению с расстоянием между взвешенными частицами, но малому относительно до макроскопических размеров. [57] Такие суспензии обычно демонстрируют неньютоновское поведение. Однако для разбавленных систем в установившихся потоках поведение является ньютоновским, и выражения для могут быть получены непосредственно из динамики частиц. В очень разбавленной системе с объемной долей взаимодействием между взвешенными частицами можно пренебречь. В таком случае можно явно рассчитать поле течения вокруг каждой частицы независимо и объединить результаты, чтобы получить . Для сфер это приводит к формуле эффективной вязкости Эйнштейна:

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

где коэффициент может зависеть от формы частиц (например, сферы, стержни, диски). [58] Однако экспериментальное определение точного значения затруднено: даже предсказание для сфер не было окончательно подтверждено, а различные эксперименты находили значения в диапазоне . Этот недостаток объясняется сложностью контроля условий эксперимента. [59]

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

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

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

Аморфные материалы

Общие кривые вязкости стекла [61]

В пределах высоких и низких температур вязкое течение в аморфных материалах (например, в стеклах и расплавах) [62] [63] [64] имеет форму Аррениуса :

где Q — соответствующая энергия активации , выраженная в молекулярных параметрах; Т – температура; R — молярная газовая постоянная ; и A приблизительно постоянна. Энергия активации Q принимает различное значение в зависимости от того, рассматривается ли верхний или нижний температурный предел: она меняется от высокого значения Q H при низких температурах (в стеклообразном состоянии) до низкого значения Q L при высоких температурах (в жидкое состояние).

Десятый логарифм зависимости вязкости от температуры для B 2 O 3 , показывающий два режима

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

где , , , – все константы, обеспечивает хорошее соответствие экспериментальным данным во всем диапазоне температур, одновременно приводя к правильной форме Аррениуса в низком и высоком температурном пределах. Это выражение можно мотивировать различными теоретическими моделями аморфных материалов на атомном уровне. [63]

Двухэкспоненциальное уравнение для вязкости может быть получено в рамках модели выталкивания Дайра для переохлажденных жидкостей, где энергетический барьер Аррениуса идентифицируется как высокочастотный модуль сдвига , умноженный на характерный объем выталкивания. [65] [66] При задании температурной зависимости модуля сдвига через тепловое расширение и отталкивающую часть межмолекулярного потенциала получается другое двухэкспоненциальное уравнение: [67]

где обозначает высокочастотный модуль сдвига материала, оцененный при температуре, равной температуре стеклования , - так называемый объем выталкивания, т.е. это характерный объем группы атомов, участвующих в событии выталкивания, при котором атом /молекула вырывается из клетки ближайших соседей, обычно порядка объема, занимаемого несколькими атомами. Кроме того, – коэффициент теплового расширения материала, является параметром, измеряющим крутизну степенного подъема восходящего склона первого пика функции радиального распределения , и количественно связан с отталкивающей частью межатомного потенциала . [67] Наконец, обозначает константу Больцмана .

Вихревая вязкость

При изучении турбулентности в жидкостях общепринятая практическая стратегия состоит в том, чтобы игнорировать мелкомасштабные вихри (или завихрения ) в движении и рассчитывать крупномасштабное движение с эффективной вязкостью, называемой «вихревой вязкостью», которая характеризует перенос и диссипация энергии в потоке меньшего масштаба (см. Моделирование больших вихрей ). [68] [69] В отличие от вязкости самой жидкости, которая согласно второму закону термодинамики должна быть положительной , вихревая вязкость может быть отрицательной. [70] [71]

Прогноз

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

Для простейших жидкостей, таких как разбавленные одноатомные газы и их смеси, квантово-механические вычисления ab initio могут точно предсказать вязкость с точки зрения фундаментальных атомных констант, т.е. без ссылки на существующие измерения вязкости. [72] Для частного случая разбавленного гелия неопределенности в расчетной вязкости ab initio на два порядка меньше, чем неопределенности в экспериментальных значениях. [73]

Для немного более сложных жидкостей и смесей при умеренных плотностях (т.е. докритических плотностях ) пересмотренная теория Энскога может использоваться для прогнозирования вязкости с некоторой точностью. [51] Пересмотренная теория Энскога является прогнозирующей в том смысле, что прогнозы вязкости могут быть получены с использованием параметров, соответствующих другим термодинамическим свойствам чистой жидкости или транспортным свойствам , что не требует априорных экспериментальных измерений вязкости.

Для большинства жидкостей высокоточные расчеты из первых принципов невозможны. Скорее, теоретические или эмпирические выражения должны соответствовать существующим измерениям вязкости. Если такое выражение подходит для высокоточных данных в широком диапазоне температур и давлений, то оно называется «эталонной корреляцией» для этой жидкости. Справочные корреляции были опубликованы для многих чистых жидкостей; несколькими примерами являются вода , углекислый газ , аммиак , бензол и ксенон . [74] [75] [76] [77] [78] Многие из них охватывают диапазоны температур и давлений, охватывающие газовые, жидкие и сверхкритические фазы.

Программное обеспечение для теплофизического моделирования часто опирается на эталонные корреляции для прогнозирования вязкости при заданных пользователем температуре и давлении. Эти корреляции могут быть собственными . Примерами являются REFPROP [79] (собственный) и CoolProp [80] (с открытым исходным кодом).

Вязкость также можно рассчитать с помощью формул, которые выражают ее через статистику траекторий отдельных частиц. Эти формулы включают соотношения Грина-Кубо для линейной сдвиговой вязкости и выражения временной корреляционной функции , полученные Эвансом и Морриссом в 1988 году. [81] [34] Преимущество этих выражений состоит в том, что они формально точны и справедливы для общих систем. . Недостатком является то, что они требуют детального знания траекторий частиц, доступного только в дорогостоящих с точки зрения вычислений симуляциях, таких как молекулярная динамика . Также необходима точная модель межчастичных взаимодействий, которую может быть сложно получить для сложных молекул. [82]

Выбранные вещества

В эксперименте с падением смолы , проведенном в Университете Квинсленда , смола медленно капала через воронку с 1927 года, со скоростью одна капля примерно каждые десять лет. Таким образом, вязкость пека была определена примерно в 230 миллиардов ((в 2,3 × 10 11 ) раз больше, чем у воды. [83]

Наблюдаемые значения вязкости различаются на несколько порядков даже для обычных веществ (см. таблицу порядков величин ниже). Например, вязкость 70%-ного раствора сахарозы (сахара) более чем в 400 раз выше вязкости воды и в 26 000 раз выше вязкости воздуха. [84] Более того, по оценкам, вязкость смолы в 230 миллиардов раз выше вязкости воды. [83]

Вода

Динамическая вязкость воды составляет около 0,89 мПа·с при комнатной температуре (25 °С) . Вязкость как функцию температуры в Кельвинах можно оценить с помощью полуэмпирического уравнения Фогеля-Фульчера-Таммана :

где A = 0,02939 мПа·с, B = 507,88 К и C = 149,3 К. [85] Экспериментально определенные значения вязкости также приведены в таблице ниже. Полезными являются значения при 20 °C: здесь динамическая вязкость составляет около 1 сП, а кинематическая вязкость — около 1 сСт.

Воздух

При стандартных атмосферных условиях (25 °C и давлении 1 бар) динамическая вязкость воздуха составляет 18,5 мкПа·с, что примерно в 50 раз меньше вязкости воды при той же температуре. За исключением очень высокого давления, вязкость воздуха зависит главным образом от температуры. Среди множества возможных приближенных формул температурной зависимости (см. Температурная зависимость вязкости ) есть одна: [86]

что соответствует точности в диапазоне от -20 °C до 400 °C. Чтобы эта формула была действительна, температура должна быть указана в кельвинах ; тогда соответствует вязкости в Па·с.

Мед поливают

Другие распространенные вещества

Оценки порядка величины

В следующей таблице показан диапазон значений вязкости, наблюдаемый в обычных веществах. Если не указано иное, предполагается температура 25 °C и давление 1 атмосфера.

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

Смотрите также

Рекомендации

Сноски

  1. ^ Последующее обсуждение основано на Chapman & Cowling 1970, стр. 232–237.
  2. ^ ab Эти материалы в высшей степени неньютоновские .

Цитаты

  1. ^ abc «Вязкость». Британская энциклопедия . 26 июня 2023 г. Проверено 4 августа 2023 г.
  2. ^ Растём вместе с наукой. Маршалл Кавендиш . 2006. с. 1928. ISBN 9780761475217.
  3. ^ Э. Дейл Мартин (1961). Исследование ламинарного течения сжимаемой вязкой жидкости в трубе, ускоренного осевой объемной силой, с применением к магнитогазодинамике. НАСА . п. 7.
  4. ^ Балеску 1975, стр. 428–429.
  5. ^ Ландау и Лифшиц 1987.
  6. ^ Харпер, Дуглас (nd). «вязкий (прил.)». Интернет-словарь этимологии . Архивировано из оригинала 1 мая 2019 года . Проверено 19 сентября 2019 г.
  7. ^ Мьюис и Вагнер 2012, с. 19.
  8. ^ Стритер, Уайли и Бедфорд 1998.
  9. ^ Холман 2002.
  10. ^ Incropera et al. 2007.
  11. ^ Нич и др. 1997.
  12. ^ abc Bird, Stewart & Lightfoot 2007, стр. 19.
  13. ^ аб Ландау и Лифшиц 1987, стр. 44–45.
  14. ^ Берд, Стюарт и Лайтфут 2007, стр. 18: В этом источнике используется альтернативное соглашение о знаках, которое здесь изменено на противоположное.
  15. ^ Ландау и Лифшиц 1987, с. 45.
  16. ^ аб Балеску 1975.
  17. ^ abcd Чепмен и Коулинг 1970.
  18. ^ Миллат 1996.
  19. ^ abcde Bird, Stewart & Lightfoot 2007.
  20. ^ Шредер 1999.
  21. ^ Ружанска и др. 2014, стр. 47–55.
  22. ^ Траутон 1906, стр. 426–440.
  23. ^ Мьюис и Вагнер, 2012, стр. 228–230.
  24. ^ Кумагай, Сасадзима и Ито 1978, стр. 157–161.
  25. ^ Шерер, Парденек и Святек 1988, с. 14.
  26. ^ Ханнан 2007.
  27. ^ Квапишевска и др. 2020.
  28. ^ Макнот и Уилкинсон 1997, уравновешенность.
  29. ^ Гилленбок 2018, с. 213.
  30. ^ «Какая единица называется рейн?». Размеры.com . Проверено 23 декабря 2023 г.
  31. ^ ASTM D2161: Стандартная практика преобразования кинематической вязкости в универсальную вязкость по Сейболту или вязкость по Сейболту-Фуролу , ASTM , 2005, стр. 1
  32. ^ Траченко и Бражкин 2020.
  33. ^ аб Траченко и Бражкин 2021.
  34. ^ аб Эванс и Моррис 1988.
  35. ^ аб Беллак, Мортессан и Батруни 2004.
  36. ^ Чепмен и Коулинг 1970, стр. 103.
  37. ^ Черчиньяни 1975.
  38. ^ Сазерленд 1893, стр. 507–531.
  39. ^ Берд, Стюарт и Лайтфут 2007, стр. 25–27.
  40. ^ Чепмен и Коулинг 1970, стр. 235–237.
  41. ^ Чепмен и Коулинг 1970, стр. 197, 214–216.
  42. ^ Крамер 2012, с. 066102-2.
  43. ^ Рид и Шервуд 1958, с. 202.
  44. ^ abc Bird, Stewart & Lightfoot 2007, стр. 29–31.
  45. ^ Рид и Шервуд 1958, стр. 203–204.
  46. ^ Хильдебранд 1977.
  47. ^ Хильдебранд 1977, с. 37.
  48. ^ Эгельстафф 1992, с. 264.
  49. ^ Ирвинг и Кирквуд 1949, стр. 817–829.
  50. ^ Рид и Шервуд 1958, стр. 206–209.
  51. ^ abc Джервелл, Вегард Г.; Вильгельмсен, Ойвинд (08 июня 2023 г.). «Пересмотренная теория Энскога для жидкостей Ми: прогнозирование коэффициентов диффузии, коэффициентов термодиффузии, вязкости и теплопроводности». Журнал химической физики . 158 (22). Бибкод : 2023JChPh.158v4101J. дои : 10.1063/5.0149865. ISSN  0021-9606. PMID  37290070. S2CID  259119498.
  52. ^ Леммон, EW; Якобсен, RT (1 января 2004 г.). «Уравнения вязкости и теплопроводности для азота, кислорода, аргона и воздуха». Международный журнал теплофизики . 25 (1): 21–69. doi :10.1023/B:IJOT.0000022327.04529.f3. ISSN  1572-9567.
  53. ^ Лопес де Аро, М.; Коэн, EGD; Кинкейд, Дж. М. (1 марта 1983 г.). «Теория Энскога для многокомпонентных смесей. I. Теория линейного переноса». Журнал химической физики . 78 (5): 2746–2759. дои : 10.1063/1.444985. ISSN  0021-9606.
  54. ^ аб Жмудь 2014, с. 22.
  55. ^ аб Вишванат и др. 2007.
  56. ^ Абдулагатов, Зейналова и Азизов 2006, стр. 75–88.
  57. ^ ab Bird, Stewart & Lightfoot 2007, стр. 31–33.
  58. ^ Берд, Стюарт и Лайтфут 2007, стр. 32.
  59. ^ ab Mueller, Llewellin & Mader 2009, стр. 1201–1228.
  60. ^ Берд, Стюарт и Лайтфут 2007, стр. 33.
  61. ^ Флюгель 2007.
  62. ^ Доремус 2002, стр. 7619–7629.
  63. ^ аб Оджован, Трэвис и Хэнд 2007, стр. 415107.
  64. ^ Оджован и Ли 2004, стр. 3803–3810.
  65. ^ Дайр, Олсен и Кристенсен 1996, с. 2171.
  66. ^ Хекшер и Дайр 2015.
  67. ^ аб Крауссер, Самвер и Закконе, 2015, стр. 13762.
  68. ^ Берд, Стюарт и Лайтфут 2007, стр. 163.
  69. ^ Лесье 2012, стр. 2–.
  70. ^ Сивашинский и Яхот 1985, с. 1040.
  71. ^ Се и Левченко 2019, с. 045434.
  72. ^ Шарипов и Бенитес 2020.
  73. ^ Роуленд, Аль Гафри и май 2020 г.
  74. ^ Хубер и др. 2009.
  75. ^ Лазеке и Музный 2017.
  76. ^ Моногениду, Ассаэль и Хубер 2018.
  77. ^ Авгери и др. 2014.
  78. ^ Веллиаду и др. 2021.
  79. ^ "Рефпроп". НИСТ . Nist.gov. 18 апреля 2013 г. Архивировано из оригинала 9 февраля 2022 г. Проверено 15 февраля 2022 г.
  80. ^ Белл и др. 2014.
  81. ^ Эванс и Моррис 2007.
  82. ^ Магинн и др. 2019.
  83. ^ abcd Edgeworth, Dalton & Parnell 1984, стр. 198–200.
  84. ^ abcde Rumble 2018.
  85. ^ Вишванатх и Натараджан 1989, стр. 714–715.
  86. ^ tec-science (25 марта 2020 г.). «Вязкость жидкостей и газов». техническая наука . Архивировано из оригинала 19 апреля 2020 г. Проверено 7 мая 2020 г.
  87. ^ Abcd Fellows 2009.
  88. ^ Янниотис, Скальци и Карабурниоти 2006, стр. 372–377.
  89. ^ аб Коочеки и др. 2009, стр. 596–602.
  90. ^ ab Citerne, Carreau & Moan 2001, стр. 86–96.
  91. ^ Кестин, Халифа и Уэйкхэм 1977.
  92. ^ Ассаэль и др. 2018.
  93. ^ Кестин, Ро и Уэйкхэм 1972.
  94. ^ Розенсон, Маккормик и Урец 1996.
  95. ^ Чжао и др. 2021.
  96. ^ Сагдеев и др. 2019.
  97. ^ Уолцер, Хендель и Баумгарднер и др.

Источники

Внешние ссылки