Модели вязкости смесей
Сдвиговая вязкость (или, короче говоря, вязкость) жидкости — это свойство материала, которое описывает трение между внутренними соседними поверхностями (или слоями) жидкости, текущими с разными скоростями жидкости. Это трение представляет собой эффект (линейного) обмена импульсом, вызванный молекулами, имеющими достаточную энергию для перемещения (или «перепрыгивания») между этими слоями жидкости из-за колебаний их движения. Вязкость — это не константа материала, а свойство материала, которое зависит от температуры, давления, состава смеси жидкостей, локальных изменений скорости. Эта функциональная зависимость описывается математической моделью вязкости, называемой определяющим уравнением , которая обычно намного сложнее, чем определяющее уравнение сдвиговой вязкости. Одной из таких усложняющих особенностей является связь между моделью вязкости для чистой жидкости и моделью для смеси жидкостей, которая называется правилами смешивания. Когда ученые и инженеры используют новые аргументы или теории для разработки новой модели вязкости, вместо улучшения существующей модели это может привести к созданию первой модели в новом классе моделей. В этой статье будут показаны одна или две репрезентативные модели для разных классов моделей вязкости, а именно:
- Элементарная кинетическая теория и простые эмпирические модели [1] [2] [3] - вязкость для разбавленного газа с почти сферическими молекулами
- Силовая серия [2] [3] - самый простой подход после разбавленного газа
- Уравнение аналогии состояний [3] между PVT и T П
- Соответствующее состояние [2] [3] модель — масштабирование переменной по ее значению в критической точке
- Теория силы трения [3] - аналогия внутренней скользящей поверхности с выдвижной коробкой на наклонной поверхности
- Много- и однопараметрическая версия теории сил трения
- Аналогия с переходным состоянием - молекулярная энергия, необходимая для того, чтобы втиснуться в вакансию, аналогична молекулам, сцепляющимся друг с другом в химической реакции.
- Теория свободного объема [3] - молекулярная энергия, необходимая для перехода на свободное место на соседней поверхности
- Теория значимой структуры [3] - основано на концепции Айринга о жидкости как о смеси твердого и газообразного поведения/свойств.
Отдельные вклады этих направлений развития представлены в следующих разделах. Это означает, что некоторые известные вклады в области исследований и разработок не включены. Например, не отображается ли метод группового вклада , примененный к модели сдвиговой вязкости. Несмотря на то, что это важный метод, считается, что он является методом параметризации выбранной модели вязкости, а не самой моделью вязкости.
Микроскопическое или молекулярное происхождение жидкостей означает, что коэффициенты переноса, такие как вязкость, могут быть рассчитаны с помощью временных корреляций , которые справедливы как для газов, так и для жидкостей, но это трудоемкие компьютерные вычисления. Другой подход - уравнение Больцмана , которое описывает статистическое поведение термодинамической системы, не находящейся в состоянии равновесия. Его можно использовать для определения того, как изменяются физические величины, такие как тепловая энергия и импульс, когда жидкость находится в транспортировке, но это требует интенсивного компьютерного моделирования.
Из уравнения Больцмана можно также аналитически вывести (аналитические) математические модели свойств, характерных для жидкостей, таких как вязкость, теплопроводность и электропроводность (путем рассмотрения носителей заряда в материале как газа). См. также уравнение конвекции-диффузии . Математика настолько сложна для полярных и несферических молекул, что очень трудно получить практические модели вязкости. Поэтому чисто теоретический подход в оставшейся части статьи будет опущен, за исключением некоторых упоминаний, связанных с разбавленным газом и существенной структурной теорией.
Использование, определение и зависимость
[ редактировать ]Классическое уравнение Навье-Стокса представляет собой уравнение баланса плотности импульса изотропной, сжимающей и вязкой жидкости, которое используется в механике жидкости в целом и в динамике жидкости в частности:
В правой части находится (расхождение) тензора полных напряжений. который состоит из тензора давления и диссипативный (или вязкий, или девиаторный) тензор напряжений . Диссипативное напряжение представляет собой тензор напряжений сжатия (член № 2) и тензор касательных напряжений (термин №3). Самый правый термин - гравитационная сила, которая представляет собой вклад объемной силы, и - массовая плотность, а это скорость жидкости.
Для жидкостей пространственная или эйларова форма основных уравнений предпочтительнее материальной или лагранжевой формы, а концепция градиента скорости предпочтительнее эквивалентной концепции тензора скорости деформации . Таким образом, предположения Стокса для широкого класса жидкостей гласят, что для изотропной жидкости напряжения сжатия и сдвига пропорциональны градиентам их скорости: и соответственно, и назвал этот класс жидкостей ньютоновскими жидкостями . Классическое определяющее уравнение объемной вязкости и сдвиговая вязкость соответственно:
Классический «градиент» скорости сжатия представляет собой диагональный тензор, который описывает сжимающий (или расширяющийся) поток или ослабляющие звуковые волны :
Классический градиент скорости сдвига Коши представляет собой симметричный и бесследовый тензор, который описывает чистый сдвиговый поток (где чистый означает исключение нормального оттока, что в математических терминах означает бесследовую матрицу), например, вокруг крыла, гребного винта, корпуса корабля или, например, в реке. труба или жила с изгибами и граничной обшивкой или без них:
где симметричная матрица градиента с ненулевым следом равна
Насколько объемная вязкость влияет на характеристики потока, например, в дроссельном потоке , таком как сужающееся-расширяющееся сопло или поток в клапане , не очень хорошо известно, но сдвиговая вязкость, безусловно, является наиболее используемым коэффициентом вязкости. Объемная вязкость теперь будет оставлена, а остальная часть статьи будет посвящена сдвиговой вязкости.
Еще одним применением моделей сдвиговой вязкости является закон Дарси для многофазного потока .
- где а = вода, нефть, газ
и и – абсолютная и относительная проницаемость соответственно. Эти 3 (векторных) уравнения моделируют течение воды, нефти и природного газа в подземных резервуарах нефти и газа в пористых породах. Хотя изменения давления велики, жидкие фазы будут медленно течь через пласт из-за ограничения потока, вызванного пористой породой.
Приведенное выше определение основано на движении жидкости под действием сдвига , которое в наиболее общей форме моделируется тензором сдвигового напряжения и тензором градиента скорости. Однако гидродинамика сдвигового потока очень хорошо иллюстрируется простым течением Куэтта . В этой экспериментальной схеме напряжение сдвига и градиент скорости сдвига (где сейчас ) принимает простой вид:
Включение этих упрощений дает нам определяющее уравнение, которое можно использовать для интерпретации экспериментальных измерений:
где - площади движущейся и неподвижной пластины, – пространственная координата, нормальная к пластинам. В этой экспериментальной установке значение силы выбирается первым. Тогда максимальная скорость измеряется, и, наконец, оба значения вводятся в уравнение для расчета вязкости. Это дает одно значение вязкости выбранной жидкости. Если выбрано другое значение силы, будет измерена другая максимальная скорость. Это приведет к другому значению вязкости, если жидкость является неньютоновской жидкостью, такой как краска, но это даст то же значение вязкости для ньютоновской жидкости, такой как вода, нефтяное масло или газ. Если другой параметр, например температура, , изменяется и эксперимент повторяется с той же силой, будет рассчитано новое значение вязкости как для неньютоновских, так и для ньютоновских жидкостей. Подавляющее большинство свойств материалов изменяются в зависимости от температуры, в том числе и вязкости. Вязкость также является функцией давления и, конечно же, самого материала. Для смеси жидкостей это означает, что сдвиговая вязкость также будет меняться в зависимости от состава жидкости . Чтобы отобразить вязкость как функцию всех этих переменных, требуется большая последовательность экспериментов, которая генерирует еще больший набор чисел, называемых измеренными данными, наблюдаемыми данными или наблюдениями . До или одновременно с экспериментами создается модель свойств материала (или краткая модель материала), предлагаемая для описания или объяснения наблюдений. Эта математическая модель называется определяющим уравнением сдвиговой вязкости. Обычно это явная функция, содержащая некоторые эмпирические параметры, которые корректируются так, чтобы соответствовать наблюдениям настолько хорошо, насколько способна математическая функция.
Для ньютоновской жидкости основное уравнение сдвиговой вязкости обычно является функцией температуры , давления и состава жидкости :
где – состав жидкой фазы с молфракцией для жидкого компонента i, и и – газовая фаза и общий состав флюида соответственно. Для неньютоновской жидкости (в смысле обобщенной ньютоновской жидкости ) материальное уравнение сдвиговой вязкости также является функцией градиента скорости сдвига:
Существование градиента скорости в функциональном соотношении для неньютоновских жидкостей говорит о том, что вязкость обычно не является уравнением состояния, поэтому термин «конституционное уравнение» обычно будет использоваться для уравнений (или функций) вязкости. Свободные переменные в двух приведенных выше уравнениях также указывают на то, что конкретные определяющие уравнения для сдвиговой вязкости будут сильно отличаться от простого определяющего уравнения для сдвиговой вязкости, которое показано ниже. Оставшаяся часть этой статьи покажет, что это, безусловно, правда. Поэтому от неньютоновских жидкостей придется отказаться, и оставшаяся часть статьи будет посвящена ньютоновским жидкостям.
Предел разбавленного газа и масштабируемые переменные
[ редактировать ]Элементарная кинетическая теория
[ редактировать ]В учебниках по элементарной кинетике теория [1] можно найти результаты моделирования разбавленных газов, которые получили широкое распространение. Вывод кинетической модели сдвиговой вязкости обычно начинается с рассмотрения течения Куэтта , в котором две параллельные пластины разделены слоем газа. Этот неравновесный поток накладывается на Максвелла – Больцмана равновесное распределение молекулярных движений .
Позволять быть сечением столкновения одной молекулы с другой. Плотность числа определяется как количество молекул на (обширный) объем . Сечение столкновения на объем или плотность сечения столкновения равна , и это связано со средней длиной свободного пробега к
Объединение кинетических уравнений молекулярного движения с определяющим уравнением сдвиговой вязкости дает хорошо известное уравнение сдвиговой вязкости для разбавленных газов:
где
где — постоянная Больцмана , – постоянная Авогадро , газовая постоянная , это молярная масса и это молекулярная масса . В приведенном выше уравнении предполагается, что плотность газа мала (т.е. давление низкое), поэтому в переменной указан нулевой индекс. . Это означает, что кинетическая поступательная энергия доминирует над вращательной и колебательной энергией молекул. Уравнение вязкости, представленное выше, дополнительно предполагает, что существует только один тип молекул газа и что молекулы газа представляют собой идеально эластичные частицы с твердым ядром сферической формы. Это предположение о том, что частицы подобны бильярдным шарам с радиусом , означает, что сечение столкновения одной молекулы можно оценить по формуле
Но молекулы не являются твердыми частицами. Для достаточно сферической молекулы потенциал взаимодействия больше похож на потенциал Леннарда-Джонса или даже больше на потенциал Морса . Оба имеют отрицательную часть, которая притягивает другую молекулу с расстояний, намного превышающих радиус твердого ядра, и, таким образом, моделирует силы Ван-дер-Ваальса. Положительная часть моделирует силы отталкивания, когда электронные облака двух молекул перекрываются. Поэтому радиус нулевого потенциала взаимодействия подходит для оценки (или определения) сечения столкновения в кинетической теории газа, а r-параметр (конф. ) поэтому называется кинетическим радиусом . D-параметр (где ) называется кинетическим диаметром .
Макроскопическое сечение столкновения часто связан с критическим молярным объемом и часто без дополнительных доказательств или подтверждающих аргументов, путем
где — это параметр формы молекулы, который принимается в качестве эмпирического параметра настройки, и включается чистая числовая часть, чтобы сделать окончательную формулу вязкости более подходящей для практического использования. Вставив эту интерпретацию и использование пониженной температуры , дает
откуда следует, что эмпирический параметр безразмерен, и это и имеют одинаковые единицы. Параметр масштабный параметр, включающий газовую постоянную и критический молярный объем , и он использовался для масштабирования вязкости. В этой статье параметр масштабирования вязкости часто будет обозначаться как которые включают один или несколько параметров , , помимо критической температуры и молярная масса . Неполные параметры масштабирования, такие как параметр выше, где газовая постоянная поглощается эмпирической константой, часто встречается на практике. В этом случае уравнение вязкости принимает вид
где эмпирический параметр не является безразмерной, и предлагаемая модель вязкости для плотной жидкости не будет безразмерной, если — общий коэффициент масштабирования. Обратите внимание, что
Подстановка критической температуры в уравнение для разбавленной вязкости дает
Значения параметров по умолчанию и должны быть достаточно универсальными ценностями, хотя зависит от системы единиц. Однако критический молярный объем в параметрах масштабирования и трудно получить из экспериментальных измерений, и это является существенным недостатком. Общее уравнение состояния реального газа обычно записывается в виде
где критический коэффициент сжимаемости , отражающее объемное отклонение реальных газов от идеального, также нелегко получить из лабораторных экспериментов. Однако критическое давление и критическая температура более доступны из измерений. Следует добавить, что критическая вязкость также не всегда может быть получена экспериментально.
Уехара и Ватсон (1944) [4] предложил поглотить универсальное среднее значение (и газовая постоянная ) в значение параметра настройки по умолчанию как практическое решение трудностей получения экспериментальных значений и/или . Тогда модель вязкости для разбавленного газа имеет вид
Подставляя критическую температуру в приведенную выше формулу, критическая вязкость рассчитывается как
На основе среднего критического коэффициента сжимаемости и измерили значения критической вязкости 60 различных типов молекул, Уехара и Уотсон (1944). [4] определил среднее значение быть
Кубическое уравнение состояния (EOS) — очень популярные уравнения, которые достаточно точны для большинства промышленных вычислений как в парожидкостном равновесии, так и в молярном объеме. Их самым слабым местом, пожалуй, является молярный объем в жидкой области и в критической области.Принимая кубическую EOS, молярный объем твердого ядра может быть рассчитана на основе ограничения точки поворота в критической точке. Это дает
где константа — универсальная константа, специфичная для выбранного варианта кубической ЭОС. Это говорит о том, что использование , и пренебрегая вариациями компонентов жидкости На практике это эквивалентно утверждению, что макроскопическое сечение столкновения пропорционально молярному объему твердого ядра, а не критическому молярному объему.
В жидкой смеси, такой как нефтяной газ или нефть, существует множество типов молекул, и внутри этой смеси существуют семейства типов молекул (т.е. группы компонентов жидкости). Самая простая группа — н-алканы, представляющие собой длинные цепочки CH 2 -элементов. Чем больше CH 2 -элементов, или атомов углерода, тем длиннее молекула. Таким образом, критическая вязкость и критические термодинамические свойства н-алканов демонстрируют тенденцию или функциональное поведение, если их построить в зависимости от молекулярной массы или числа атомов углерода в молекуле (т.е. числа углерода). Параметры в уравнениях для таких свойств, как вязкость, обычно также демонстрируют такое же поведение. Это означает, что
Это говорит о том, что параметр масштабирования сам по себе не является истинным или полным коэффициентом масштабирования, если только все компоненты жидкости не имеют достаточно схожей (и предпочтительно сферической) формы.
Наиболее важным результатом этого кинетического вывода является, возможно, не формула вязкости, а полуэмпирический результат. параметр который широко используется в промышленности и прикладных научных кругах в качестве коэффициента масштабирования (сдвиговой) вязкости. В литературе часто сообщается об обратном параметре и обозначается как .
Вклад вязкости разбавленного газа в общую вязкостьжидкости будет важно только при прогнозировании вязкости паров при низких давлениях или вязкости плотных жидкостей.при высоких температурах. Показанная выше модель вязкости для разбавленного газа широко используется в промышленности и прикладных научных кругах. Поэтому многие исследователи не указывают модель вязкости разбавленного газа, когда предлагают модель общей вязкости, а оставляют пользователю право выбора и включения вклада разбавленного газа. Некоторые исследователи не включают отдельный термин в модель разбавленного газа, а предлагают модель общей вязкости газа, охватывающую весь исследованный ими диапазон давлений и температур.
В этом разделе нашими основными макроскопическими переменными и параметрами и их единицами являются температура. [К], давление [бар], молярная масса [г/моль], вязкость газа низкой плотности (низкого давления или разбавленного газа) [мкП]. Однако в отрасли принято использовать другую единицу измерения вязкости жидкости и газа высокой плотности. [КП].
Кинетическая теория
[ редактировать ]На основе уравнения Больцмана Чепмен и Энског вывели модель вязкости разбавленного газа .
где (абсолютное значение) энергетической глубины потенциальной ямы (см., например, потенциал взаимодействия Леннарда-Джонса ). Термин называется интегралом столкновений и возникает как общая функция температуры, которую должен указать пользователь, а это непростая задача. Это иллюстрирует ситуацию с молекулярным или статистическим подходом: (аналитическая) математика становится невероятно сложной для полярных и несферических молекул, что очень затрудняет создание практических моделей вязкости, основанных на статистическом подходе. Поэтому чисто статистический подход в оставшейся части статьи будет опущен.
Эмпирическая корреляция
[ редактировать ]Зеберг-Миккельсен (2001) [3] предложил эмпирические модели вязкости газа для достаточно сферических молекул, которые представлены в разделе, посвященном теории силы трения и ее моделям для разбавленных газов и простых легких газов. Эти простые эмпирические корреляции показывают, что эмпирические методы конкурируют со статистическим подходом в отношении моделей вязкости газа для простых жидкостей (простых молекул).
Кинетическая теория с эмпирическим расширением
[ редактировать ]Модель вязкости газа Chung et alios (1988). [5] представляет собой комбинацию Чепмена -Энскога (1964).кинетическая теория вязкости для разбавленных газов и эмпирическое выражение Нойфельд и другие (1972) [6] для уменьшенного интеграла столкновений, но расширенного эмпирического для многоатомных, полярных и водородныхСвязующие жидкости в широком диапазоне температур. Эта модель вязкости иллюстрирует успешное сочетание кинетической теории и эмпиризма и отображается в разделе Теория существенных структур и ее модель газоподобного вклада в общую вязкость жидкости.
Функции тренда и масштабирование
[ редактировать ]В разделе моделей, основанных на элементарной кинетической теории, обсуждалось несколько вариантов масштабирования уравнения вязкости, и они показаны ниже для компонента жидкости i, в качестве услуги читателю.
Зеберг-Миккельсен (2001) [3] предложил эмпирическую корреляцию параметр для н-алканов, который
Критический молярный объем компонента i связано с критической мольной плотностью и критическая мольная концентрация по уравнению . Из приведенного выше уравнения для отсюда следует, что
где — коэффициент сжимаемости компонента i, который часто используется в качестве альтернативы . Путем установления функции тренда для параметра для гомологического ряда, групп или семейств молекул значения параметров для неизвестных компонентов жидкости в гомологичной группе можно найти путем интерполяции и экстраполяции, а значения параметров можно легко сгенерировать заново при необходимости. Использование функций тренда для параметров гомологичных групп молекул значительно повысило полезность уравнений вязкости (и термодинамических УС) для смесей жидкостей, таких как нефтяной газ и нефть. [2]
Уехара и Ватсон (1944) [4] предложил корреляцию критической вязкости (для компонента жидкости i) н-алканов с использованием их среднего параметра и классический параметр масштабирования с доминированием давления :
Зеберг-Миккельсен (2001) [3] предложил эмпирическую корреляцию для параметра критической вязкости η ci для н-алканов, которая равна
Единичные уравнения для двух определяющих уравнений, приведенных Зебергом-Миккельсеном (2001), имеют вид
Включение критической температуры в три уравнения вязкости из элементарной кинетической теории дает уравнения с тремя параметрами.
Три уравнения вязкости теперь объединяются в одно уравнение вязкости.
потому что безразмерное масштабирование для всего уравнения вязкости используется . Стандартное рассуждение о безразмерности выглядит следующим образом: создание безразмерных переменных (с индексом D) путем масштабирования дает
Утверждение о безразмерности дает
Сечение столкновения и критический молярный объем, к которым трудно получить экспериментальный доступ, избегают или обходят. С другой стороны, критическая вязкость появилась как новый параметр, и критическая вязкость так же труднодоступна экспериментально, как и два других параметра. К счастью, лучшие уравнения вязкости стали настолько точными, что оправдывают расчеты в критической точке, особенно если уравнение согласовано с окружающими экспериментальными точками данных.
Классические правила смешивания
[ редактировать ]Классические правила смешивания газа
[ редактировать ]Вилке (1950) [7] вывел правило смешивания на основе кинетической теории газа
Правило смешивания Уилке способно описать правильное поведение вязкости газовых смесей, демонстрируя нелинейное и немонотонное поведение или демонстрируя характерную форму выступа, когда вязкость построена в зависимости от массовой плотности при критической температуре, для смесей, содержащих молекулы очень разных размеры. Из-за своей сложности он не получил широкого распространения. Вместо этого было использовано несколько более простое правило смешивания, предложенное Хернингом и Зипперером (1936): [8] оказывается пригодным для газов смесей углеводородов.
Классические правила смешивания жидкости
[ редактировать ]Классический Аррениус (1887). [9] Правило смешивания жидких смесей:
где - вязкость жидкой смеси, - вязкость (уравнение) для компонента жидкости i при движении в виде чистой жидкости, и – мольфракция компонента i в жидкой смеси.
Грюнберг-Ниссан (1949 г.) [10] Правило смешивания расширяет правило Аррениуса до
где — эмпирические коэффициенты бинарного взаимодействия, специальные для теории Грюнберга-Ниссана. Бинарные коэффициенты взаимодействия широко используются в кубических EOS, где они часто используются в качестве параметров настройки, особенно если компонент j является неопределенным компонентом (т.е. имеет неопределенные значения параметров).
Катти Чаудхари (1964) [11] правило смешивания
где - парциальный мольный объем компонента i, а представляет собой молярный объем жидкой фазы и получается из расчета парожидкостного равновесия (VLE) или EOS для однофазной жидкости.
Модификация правила смешивания Катти-Чаудри:
где – избыточная энергия активации вязкого течения, – энергия, характерная для межмолекулярных взаимодействий между компонентом i и компонентом j и, следовательно, ответственная за избыточную энергию активации вязкого течения. Это правило смешивания теоретически обосновано представлением Эйринга о вязкости чистой жидкости согласно Гласстону и др. (1941). [12] Количество было получено из выражения временной корреляции для сдвиговой вязкости Цванцигом (1965). [13]
Силовая серия
[ редактировать ]Очень часто просто выбирают известную корреляцию для вязкости разреженного газа. и вычитает этот вклад из общей вязкости, измеряемой в лаборатории. Это дает термин остаточной вязкости, часто обозначаемый , который представляет вклад плотной жидкости, .
Таким образом, вязкость плотной жидкости определяется как вязкость, превышающая вязкость разбавленного газа. Этот метод часто используется при разработке математических моделей как чисто эмпирических корреляций, так и моделей, имеющих теоретическое обоснование. Вклад вязкости разбавленного газа становится важным при приближении к пределу нулевой плотности (т.е. пределу нулевого давления). Также очень часто вязкость плотной жидкости масштабируют по критической вязкости или по оценке критической вязкости, которая является характерной точкой далеко в области плотной жидкости. Простейшая модель вязкости плотной жидкости представляет собой (усеченный) степенной ряд приведенной мольной плотности или давления. Джосси и др. (1962) [14] представил такую модель, основанную на уменьшенной мольной плотности, но наиболее распространенной ее формой является версия, предложенная Лоренц и др. (1964) [15] который отображается ниже.
Затем функция LBC разлагается в (усеченный) степенной ряд с эмпирическими коэффициентами, как показано ниже.
Таким образом, окончательное уравнение вязкости имеет вид
Список местной номенклатуры:
- : мольная плотность [моль/см 3 ]
- : пониженная мольная плотность [1]
- : молярная масса [г/моль]
- : критическое давление [атм]
- : температура [К]
- : критическая температура [К]
- : критический молярный объем [см 3 /моль]
- : вязкость [сП]
Смесь
[ редактировать ]Формула для который был выбран LBC, отображается в разделе «Вклад разбавленного газа».
Правила смешивания
[ редактировать ]Индекс C7+ относится к совокупности молекул углеводородов в пластовой жидкости с нефтью и/или газом, которые имеют 7 или более атомов углерода в молекуле. Критический объем фракции С7+ имеет единицу фута. 3 /фунт моль, и он рассчитывается по формуле
где – удельный вес фракции С7+.
Молярная масса (или молекулярная масса) обычно не включается в формулу EOS, но обычно входит в характеристику параметров EOS.
ЭОС
[ редактировать ]Из уравнения состояния рассчитывается мольный объем пластовой жидкости (смеси).
Молярный объем преобразуется в мольную плотность (также называемая мольной концентрацией и обозначаемая ), а затем масштабируется для уменьшения мольной плотности .
Вклад разбавленного газа
[ редактировать ]Корреляция для вязкости разбавленного газа смеси взята из работы Хернинга и Зипперера. (1936) [8] и есть
Соотношение вязкости разбавленного газа отдельных компонентов взято из Стиль и Тодос (1961) [16] и есть
где
Соответствующий государственный принцип
[ редактировать ]Принцип соответствующих состояний (принцип CS или CSP) был впервые сформулирован Ван дер Ваальсом и гласит, что две жидкости (индекс a и z) группы (например, жидкости неполярных молекул) имеют примерно одинаковый приведенный молярный объем. (или уменьшенный коэффициент сжимаемости) при сравнении при одинаковой пониженной температуре и пониженном давлении. В математическом плане это
Когда общий принцип CS, описанный выше, применяется к вязкости, он гласит:
Обратите внимание, что принцип CS изначально был сформулирован для состояний равновесия, но теперь он применяется к свойству переноса — вязкости, и это говорит нам о том, что для вязкости может потребоваться другая формула CS.
Чтобы увеличить скорость вычислений вязкости на основе теории CS, что важно, например, при моделировании композиционного коллектора, сохраняя при этом точность метода CS, Pedersen et al. (1984, 1987, 1989) [17] [18] [2] предложил метод CS, который использует простую (или традиционную) формулу CS при расчете приведенной массовой плотности, которая используется в константах вращательной связи (показана в разделах ниже), и более сложную формулу CS, включающую константы вращательной связи, в других местах. .
Смесь
[ редактировать ]Простой принцип соответствующего состояния расширен за счет включения коэффициента вращательной связи. как предложили Тэм и Габбинс (1970). [19] Эталонной жидкостью является метан, ему присвоен индекс z.
Правила смешивания
[ редактировать ]Условия взаимодействия критической температуры и критического объема:
Параметр обычно неопределенна или недоступна. Поэтому желательно избегать этого параметра. Замена с общим средним параметром для всех компонентов дает
Приведенное выше выражение для теперь подставляется в уравнение для . Это дает следующее правило смешивания
Аналогично устанавливают правило смешивания для критического давления смеси.
Правило смешивания молекулярной массы гораздо проще, но оно не совсем интуитивно понятно. Это эмпирическая комбинация более интуитивных формул с массовым взвешиванием. и мольное взвешивание .
Параметр вращательной связи смеси равен
Эталонная жидкость
[ редактировать ]Точность конечной вязкости метода CS требует очень точного прогноза плотности эталонной жидкости. Поэтому молярный объем эталонной жидкости метана рассчитывается с помощью специального EOS, а Бенедикт-Уэбб-Рубин (1940) [20] вариант уравнения состояния, предложенный Маккарти (1974), [21] и сокращенно BWRM, рекомендован Pedersen et al. (1987) для этой цели. Это означает, что массовая плотность флюида в ячейке сетки модели пласта может быть рассчитана, например, с помощью кубического EOS или с помощью входной таблицы с неизвестным установлением. Поэтому, чтобы избежать итеративных вычислений, эталонная (массовая) плотность, используемая в параметрах вращательной муфты, рассчитывается с использованием более простого принципа соответствующего состояния, который гласит, что
Молярный объем используется для расчета массовой концентрации, которая называется (массовой) плотностью, а затем масштабируется до приведенной плотности, которая равна обратной величине приведенного молярного объема, поскольку существует только один компонент (тип молекулы). В математическом плане это
Формула для параметра вращательной связи смеси показана далее, а параметр вращательной связи для эталонной жидкости (метана) равен
Массовая плотность метана, используемая в формулах вязкости, основана на расширенном соответствующем состоянии, показанном в начале этой главы о CS-методах. Используя BWRM EOS, молярный объем эталонной жидкости рассчитывается как
Опять же, молярный объем используется для расчета массовой концентрации или массовой плотности, но эталонная жидкость представляет собой однокомпонентную жидкость, и приведенная плотность не зависит от относительной молярной массы. В математическом плане это
Эффект изменения состава, например, жидкой фазы, связан с масштабными коэффициентами вязкости, температуры и давления, и это является соответствующим принципом состояния.
Эталонная корреляция вязкости Pedersen et al. (1987) [18] является
Формулы для , , взяты из Hanley et al. (1975). [22]
Вклад разбавленного газа
Температурно-зависимый коэффициент вклада первой плотности равен
Термин плотная жидкость
где показательная функция записывается как и как . Молярный объем эталонной жидкости метана, который используется для расчета массовой плотности в приведенных выше формулах вязкости, рассчитывается при пониженной температуре, пропорциональной приведенной температуре смеси. Из-за высоких критических температур более тяжелых молекул углеводородов пониженная температура более тяжелых пластовых нефтей (т.е. смесей) может давать передаваемую пониженную температуру метана, близкую к температуре замерзания метана. Это проиллюстрировано на примере двух довольно тяжелых молекул углеводородов в таблице ниже. Выбранные температуры представляют собой типичную температуру нефтяного или газового пласта, эталонную температуру Международного стандарта метрических условий для природного газа (и аналогичных жидкостей) и температуру замерзания метана ( ).
Педерсен и др. (1987) добавили четвертый член, который корректирует эталонную формулу вязкости при низких пониженных температурах. Температурные функции и являются весовыми коэффициентами. Срок их коррекции
Уравнение аналогии состояний
[ редактировать ]Филлипс (1912) [23] построенная температура в зависимости от вязкости для разных изобар пропана и заметил сходство этих изобарических кривых с классическими изотермическими кривыми поверхность. Позже Литтл и Кеннеди (1968) [24] разработал первую модель вязкости, основанную на аналогии между и с использованием Ван-дер-Ваальсовой ЭОС. Ван-дер-Ваальсова EOS была первой кубической EOS, но кубическая EOS с годами была улучшена и теперь составляет широко используемый класс EOS. Таким образом, Го и др. (1997) [25] разработали две новые модели-аналоги для вязкости на основе PR EOS (Пэн и Робинсон, 1976) и PRPT EOS (Патель и Теджа, 1982). [26] соответственно. В следующем году Т.-М. Го (1998) [27] [3] немного модифицировал модель вязкости, основанную на PR, и именно эта версия будет представлена ниже как представитель моделей вязкости, аналогичных EOS.
PR EOS отображается в следующей строке.
Уравнение вязкости Го (1998) отображается в следующей строке.
Чтобы подготовиться к правилам смешивания, уравнение вязкости переписывается для одного компонента жидкости i.
Подробности о том, как составные элементы уравнения связаны с основными параметрами и переменными, показаны ниже.
Смесь
[ редактировать ]Правила смешивания
[ редактировать ]Теория силы трения
[ редактировать ]Многопараметрическая теория сил трения
[ редактировать ]Многопараметрическая версия теории силы трения (короткая теория FF и модель FF), также называемая теорией трения (короткая F-теория), была разработана Quiñones-Cisneros et al. (2000, 2001a, 2001b и Z 2001, 2004, 2006), [28] [29] [30] [3] [31] [32] и его основные элементы с использованием некоторых хорошо известных кубических EOS показаны ниже.
Распространенным методом моделирования является принятие модели вязкости для разбавленного газа ( ), а затем построить модель вязкости плотной жидкости . Теория FF утверждает, что для жидкости, находящейся в сдвиговом движении, напряжение сдвига (т.е. силу сопротивления), действующую между двумя движущимися слоями, можно разделить на термин вызванные столкновениями разреженных газов, и термин вызванное трением в плотной жидкости.
Вязкость разбавленного газа (т. е. поведение предельной вязкости, когда давление, нормальное напряжение, стремится к нулю) и вязкость плотной жидкости (остаточная вязкость) можно рассчитать по формуле
где ду/ди — локальный градиент скорости, ортогональный направлению потока. Таким образом
Основная идея QZS (2000) заключается в том, что внутренние поверхности в потоке Куэтта действуют как (или аналогичны) механические плиты с силами трения, действующими на каждую поверхность, когда они скользят мимо друг друга. Согласно закону трения Амонтона-Кулона в классической механике соотношение кинетической силы трения и нормальная сила дается
где известен как кинетический коэффициент трения, A — площадь внутренней поверхности потока, напряжение сдвига и нормальное напряжение (или давление ) между соседними слоями в потоке Куэтта .
Теория QZS FF утверждает, что когда жидкость подвергается сдвиговому движению, межмолекулярные силы притяжения и отталкивания будут способствовать усилению или уменьшению механических свойств жидкости. Термин напряжения сдвига при трении Таким образом, можно считать, что плотная жидкость состоит из притягивающего вклада трения и сдвига. и вклад отталкивающего трения и сдвига . Вставка этого дает нам
Известные кубические уравнения состояний ( SRK, PR и PRSV EOS ), в общем виде можно записать как
Пара параметров (u,w)=(1,0) дает EOS SRK , а (u,w)=(2,-1) дает EOS PR и EOS PRSV, поскольку они различаются только температурой и составом. зависимый параметр/функция а. Входными переменными в нашем случае являются давление (P), температура (T), а для смесей также состав жидкости, который может быть однофазным (или полным) составом. , состав пара (газа) или жидкий (в нашем примере масляный) состав . Выход представляет собой мольный объем фазы (V). Поскольку кубическое EOS не является идеальным, молярный объем является более неопределенным, чем значения давления и температуры.
ЭОС состоит из двух частей, связанных с силами Ван-дер-Ваальса, или взаимодействиями, возникающими в статических электрических полях сталкивающихся частей/пятен двух (или более) сталкивающихся молекул. Отталкивающая часть ЭОС обычно моделируется как твердое ядро молекул, отсюда и обозначение (P h ), а притягивающая часть (Pa ) основана на притягивающем взаимодействии между молекулами (по конф. силе Ван-дер-Ваальса ). Таким образом, EOS можно записать как
Предположим, что молярный объем (V) известен из расчетов EOS и предыдущих расчетов парожидкостного равновесия (VLE) для смесей. Тогда две функции и могут быть использованы, и ожидается, что эти функции будут более точными и надежными, чем сам молярный объем (V). Эти функции
Поэтому теория трения предполагает, что остаточное напряжение притяжения и остаточное напряжение отталкивания являются функциями члена давления притяжения и член отталкивающего давления , соответственно.
Первая попытка — это, конечно, попробовать линейную функцию в терминах/функциях давления.
Все Коэффициенты, как правило, являются функциями температуры и состава и называются функциями трения. Оказалось, что для достижения высокой точности в широком диапазоне давлений и температур необходим член второго порядка даже для неполярных типов молекул, таких как углеводородные жидкости в нефтяных и газовых пластах, чтобы достичь высокой точности при очень высоких значениях. высокие давления. Испытание с предположительно сложной трехкомпонентной смесью типов неполярных молекул требовало мощности третьего порядка для достижения высокой точности при самых экстремальных сверхкритических давлениях.
В этой статье основное внимание будет уделено версии второго порядка, но по мере возможности будет включаться термин третьего порядка, чтобы показать полный набор формул. В качестве введения в обозначение смеси приведенное выше уравнение повторяется для компонента i в смеси.
Единичные уравнения для центральных переменных в многопараметрической FF-модели:
Функции трения
[ редактировать ]Функции трения для i-го компонента жидкости в 5-параметрической модели для чистых молекул н-алкана представлены ниже.
Функции трения для жидкого компонента i в 7- и 8-параметрических моделях представлены ниже.
Эмпирические константы в функциях трения называются константами трения. Константы трения для некоторых н-алканов в пятипараметрической модели с использованием SRK и PRSV EOS (и, следовательно, PR EOS) представлены в таблицах ниже. Константы трения для некоторых н-алканов в 7-параметрической модели с использованием PRSV EOS также представлены в таблице ниже. Константа для трех компонентов жидкости представлены ниже в последней таблице данной серии таблиц.
Смесь
[ редактировать ]В однофазных областях мольный объем жидкой смеси определяется входными переменными, такими как давление (P), температура (T) и (общий) состав жидкости. . В двухфазной газожидкостной области расчет парожидкостного равновесия (VLE) разделяет жидкость на паровую (газовую) фазу с составом и мольфракция смеси фаз n g и жидкая фаза (в нашем примере нефть) состава и молфракция смеси фаз n o . Для жидкой фазы, паровой фазы и однофазной жидкости связь с переменными VLE и EOS равна
В симуляторе композиционного пласта давление рассчитывается динамически для каждой ячейки сетки и каждого временного шага. Это дает динамическое давление для пара и жидкости (масла) или однофазной жидкости. Предполагая нулевое капиллярное давление между углеводородной жидкостью (нефтью) и газом, код программного обеспечения симулятора выдаст единое динамическое давление. которое относится как к паровой смеси, так и к жидкой (масляной) смеси. В этом случае код программного обеспечения симулятора резервуара может использовать
или
Модель трения для вязкости смеси:
Член кубической степени необходим только тогда, когда в смесь включены молекулы с достаточно жесткой двумерной структурой или пользователю требуется очень высокая точность при чрезвычайно высоких давлениях. Стандартная модель включает только линейные и квадратичные члены в функциях давления.
Правила смешивания
[ редактировать ]где эмпирическая весовая доля равна
Рекомендуемые значения для являются
- показал лучшую производительность для SRK EOS
- показал лучшую производительность для PRSV EOS
Эти значения установлены для бинарных смесей н-алканов с использованием 5-параметрической модели вязкости и, по-видимому, также используются для 7- и 8-параметрических моделей. Мотивация этого весового параметра , и, таким образом, -параметр заключается в том, что в асимметричных смесях, таких как CH 4 - C 10 H 12 , самый легкий компонент имеет тенденцию уменьшать вязкость смеси более чем линейно, если построить график зависимости от мольфракции легкого компонента (или тяжелого компонента).
Коэффициенты трения некоторых выбранных компонентов жидкости представлены в таблицах ниже для 5,7 и 8-параметрических моделей. Для удобства в таблицы также включены критические вязкости.
.
Однопараметрическая теория силы трения
[ редактировать ]Однопараметрическая версия теории силы трения (теория FF1 и модель FF1) была разработана Quiñones-Cisneros et al. (2000, 2001а, 2001б и З 2001, 2004), [28] [29] [30] [3] [31] и его основные элементы с использованием некоторых хорошо известных кубических EOS показаны ниже.
Первым шагом является определение приведенной вязкости плотной жидкости (или фрикционной) вязкости для чистой (т.е. однокомпонентной) жидкости путем деления на критическую вязкость. То же самое касается вязкости разбавленного газа.
Второй шаг – замена функций притяжения и отталкивания приведенными функциями давления. Это, конечно, повлияет и на функции трения. Поэтому вводятся новые функции трения. Их называют приведенными функциями трения, и они носят более универсальный характер. Приведенная фрикционная вязкость равна
Возвращаясь к неприводимой фрикционной вязкости и перефразируя формулу, получим
Критическая вязкость измеряется редко, и попыток предсказать ее по формулам немного. Для чистой жидкости или компонента i в смеси жидкостей для оценки критической вязкости часто используется формула кинетической теории.
где является константой, а критический молярный объем V ci считается пропорциональным сечению столкновения. Критический молярный объем V ci является значительно более неопределенным, чем параметры P ci и T ci . Чтобы избавиться от V ci , критический коэффициент сжимаемости Z ci часто заменяют универсальным средним значением. Это дает
где является константой. Основываясь на среднем критическом коэффициенте сжимаемости Z c = 0,275 и измеренных значениях критической вязкости 60 различных типов молекул, Уехара и Уотсон (1944) [4] определил среднее значение K p как
Зеберг-Миккельсен (2001) предложил эмпирическую корреляцию V ci с параметрами н-алканов, которая равна
где . Из приведенного выше уравнения и определения коэффициента сжимаемости следует, что
Зеберг-Миккельсен (2001) также предложил эмпирическую корреляцию для η ci с параметрами н-алканов, которая равна
Единичные уравнения для двух определяющих уравнений, приведенных Зебергом-Миккельсеном (2001), имеют вид
Следующим шагом является разделение формул на формулы для четко определенных компонентов (обозначенных индексом d) с учетом критической вязкости и формулы для неопределенных компонентов (обозначенных индексом u), где критическая вязкость оценивается с помощью и универсальная константа который будет рассматриваться как параметр настройки текущей смеси. Вязкость плотной жидкости (для жидкого компонента i в смеси) записывается как
Формулы теории трения затем связаны с четко определенными и неопределенными компонентами жидкости. Результат
Однако для получения характеристической критической вязкости тяжелых псевдокомпонентов можно использовать следующую модификацию выражения Уехары и Уотсона (1944) для критической вязкости. Тогда фрикционная (или остаточная) вязкость записывается как
Единичные уравнения: и и .
Функции пониженного трения
[ редактировать ]Единичное уравнение является .
Однопараметрическая модель была разработана на основе однокомпонентных жидкостей в ряду от метана до н-октадекана (от C 1 H 4 до C 18 H 38 ). Эмпирические параметры приведенных выше функций приведенного трения рассматриваются как универсальные константы и перечислены в следующей таблице. Для удобства в таблицы для моделей с 5 и 7 параметрами включены критические вязкости, представленные далее.
.
Смесь
[ редактировать ]Вязкость смеси определяется выражением
Вязкость смеси четко определенных компонентов определяется выражением
Функция вязкости смеси неопределенных компонентов определяется выражением
Вязкость смеси можно настроить в соответствии с измеренными данными вязкости путем оптимизации (регрессии) параметра .
где коэффициенты трения смеси получаются по формулам (I.7.45)–(I.7.47) и и являются членами давления притяжения и отталкивания смеси.
Правила смешивания
[ редактировать ]Правила смешивания четко определенных компонентов таковы:
QZS рекомендует отказаться от термина «разбавленный газ» для неопределенных компонентов жидкости, которыми обычно являются более тяжелые (углеводородные) компоненты. Формула сохранена здесь для единообразия. Правила смешивания неопределенных компонентов таковы:
Предел разбавления газа
[ редактировать ]Зеберг-Миккельсен (2001) [3] предложил эмпирическую модель вязкости разбавленного газа для довольно сферических молекул следующим образом.
или
Единичные уравнения для вязкости и температуры:
Второй член является поправочным членом для высоких температур. Обратите внимание, что большинство параметры отрицательные.
.
Легкие газы
[ редактировать ]Зеберг-Миккельсен (2001) предложил FF-модель вязкости легких газов следующим образом:
Функции трения для легких газов просты.
ФФ-модель для легких газов справедлива для низких, нормальных, критических и сверхкритических условий для этих газов. Хотя рекомендуется использовать модель FF для вязкости разбавленного газа, с хорошими результатами можно также использовать любую точную модель вязкости для разбавленного газа.
Единичные уравнения для вязкости и температуры:
.
Аналогия переходного состояния
[ редактировать ]Эта статья началась с вязкости смесей с отображения уравнений для разбавленного газа, основанных на элементарной кинетической теории, теории твердого ядра (кинетической), и перешла к избранным теориям (и моделям), направленным на моделирование вязкости для плотных газов, плотных жидкостей и сверхкритических жидкостей. Многие или большинство из этих теорий основаны на философии того, как газы ведут себя с молекулами, летающими вокруг, сталкивающимися с другими молекулами и обменивающимися (линейными) импульсами и, таким образом, создающими вязкость. Когда жидкость стала жидкой, модели начали отклоняться от измерений, поскольку небольшая ошибка в расчете молярного объема по EOS связана с большим изменением давления и наоборот, а, следовательно, и вязкости. Теперь статья подошла к другому концу, где теории (или модели) основаны на философии того, как жидкость ведет себя и приводит к вязкости. Поскольку молекулы в жидкости расположены гораздо ближе друг к другу, можно задаться вопросом, как часто молекула на одной скользящей поверхности жидкости находит свободный объем на соседней скользящей поверхности, достаточно большой для того, чтобы молекула могла в него прыгнуть. Это можно перефразировать так: когда молекула в своих колебательных движениях обладает достаточной энергией, чтобы втиснуться в небольшой открытый объем на соседней поверхности скольжения, подобно молекуле, которая сталкивается с другой молекулой и замыкается в ней в результате химической реакции, и, таким образом, создает новое соединение, как смоделировано в теория переходного состояния (теория ТС и модель ТС).
Теория свободного объема
[ редактировать ]Теория свободного объема (короткая теория FV и модель FV) берет свое начало от Дулитл (1951) [33] который предположил, что вязкость связана с долей свободного объема аналогично уравнению Аррениуса . Модель вязкости Дулитла (1951)
где молярный объем и – молярный объем твердого ядра.
Однако там не было особой активности в отношении теории FV до тех пор, пока Аллал и др. (1996, 2001а) [34] [35] предложил связь между долей свободного объема и параметрами (и/или переменными) на молекулярном уровне жидкости (также называемой микроструктурой жидкости). Модель 1996 года стала началом периода высокой исследовательской активности, когда были выдвинуты различные модели. Сохранившаяся модель была представлена Allal et al. (2001б), [36] и эта модель будет отображаться ниже.
Модель вязкости состоит из вклада разбавленного газа. (или ) и вклад плотной жидкости (или вклад плотного состояния или ).
Аллал и др. (2001б) [36] показали, что вклад плотной жидкости в вязкость может быть связан с коэффициентом трения. поверхности скольжения жидкости и Дуллиенс (1963). [37] показал, что коэффициент самодиффузии связано с коэффициентом трения внутренней поверхности жидкости. Эти два отношения показаны здесь:
Устранив коэффициент трения , Бонед и др. (2004) [38] выражал характерную длину как
Правая часть соответствует так называемому инварианту Дюллиена, выведенному Дюллиеном (1963, 1972). [37] [39] В результате характерная длина интерпретируется как среднее расстояние передачи импульса молекуле, которая войдет в участок свободного объема и столкнется с соседней молекулой.
Коэффициент трения смоделировано Allal et alios (2001b). [36] как
Доля свободного объема теперь связана с энергией E соотношением
где - это полная энергия, которую молекула должна использовать, чтобы диффундировать в свободный объем, и связан с работой (или энергией), необходимой для образования или расширения свободного объема, доступного для диффузии молекулы. Энергия - это энергия барьера, которую молекула должна преодолеть, чтобы диффундировать, и она смоделирована как пропорциональная массовой плотности, чтобы улучшить соответствие измеренных данных вязкости. Обратите внимание, что чувствительный термин в знаменателе модели Дулитла (1951) исчезло, что сделало модель вязкости Аллала и других (2001b) более устойчивой к численным расчетам молярного объема жидкости с помощью несовершенного EOS. Предэкспоненциальный множитель A теперь является функцией и становится
Модель вязкости, предложенная Аллалом и др. (2001b). [36] таким образом
Отступление: коэффициент самодиффузии Boned et al. (2004) [38] становится
Список местной номенклатуры:
- параметр, характеризующий перекрытие свободного объема или параметр эмпирической настройки [1]
- молярный объем твердого ядра [м 3 /моль]
- полная энергия, которую молекула должна использовать для диффузии [Дж/моль]
- барьерная энергия, которую молекула должна преодолеть для диффузии [Дж/моль]
- среднее расстояние передачи импульса для молекулы, которая передает линейный импульс (конф. радиус твердого ядра) и/или угловой момент (конф. радиус вращения) [Å]
- длина диссипации к энергии E [Å]
- составной параметр, характеризующий вязкость [Å]
- молярная масса, конф. молекулярная масса [кг/моль]
- постоянная Авогадроса
- давление [МПа]
- газовая постоянная R = 8,31451 [К·Дж/моль]
- молярный объем [м 3 /моль]
- характеристический параметр или эмпирический параметр настройки [1]
- вязкость [Па]
- массовая плотность [кг/м 3 ]
- коэффициент трения молекулы, связанный с подвижностью молекулы [1]
- коэффициент трения при нулевой массовой плотности, т.е. для разбавленной системы/предел низкого давления [1]
Смесь
[ редактировать ]Вязкость смеси
Вязкость разбавленного газа взято из Chung et al. (1988). [40] который отображается в разделе теории СС. Вклад плотной жидкости в вязкость в теории FV равен
где – три характерных параметра расчета вязкости жидкости. Для смесей жидкостей эти три параметра рассчитываются с использованием правил смешивания. Если коэффициент самодиффузии включен в основные уравнения, возможно, через уравнение диффузии, использование четырех характеристических параметров (т.е. использование L p и L d вместо L c ) даст согласованную модель потока, но исследования потока, которые включают уравнение диффузии принадлежит к небольшому классу специальных исследований.
Единицей вязкости является [Пас], тогда как все остальные единицы измерения сохраняются в единицах СИ.
Правила смешивания
[ редактировать ]В конце периода интенсивных исследований Аллал и др. (2001с) [41] и Кане (2001) [42] предложил два разных набора правил смешивания, и, согласноАлмаси (2015) [43] В литературе не существует единого мнения относительно того, какие правила смешивания являются лучшими. Поэтому Алмаси (2015) рекомендовал классические линейные моль-взвешенные правила смешивания, которые показаны ниже для смеси N флюидных компонентов.
Три характеристических параметра вязкости обычно устанавливаются путем оптимизации формулы вязкости на основе измеренных данных вязкости для чистых жидкостей (т. е. однокомпонентных жидкостей).
Функции тренда
[ редактировать ]Три характеристических параметра вязкости обычно устанавливаются путем оптимизации формулы вязкости на основе измеренных данных вязкости для чистых жидкостей (т. е. однокомпонентных жидкостей). Данные по этим параметрам затем можно хранить в базах данных вместе с данными и информацией о других химических и физических свойствах материалов. Это происходит чаще, если использование уравнения становится широко распространенным. Молекулы углеводородов — это огромная группа молекул, состоящая из нескольких подгрупп, которые сами содержат молекулы одинаковой базовой структуры, но разной длины. Алканы — простейшие из этих групп. Свойство материала молекул в такой группе обычно отображается как функция при построении графика зависимости от другого свойства материала. Затем на основе физических/химических знаний, опыта и интуиции выбирается математическая функция, а эмпирические параметры (т.е. константы) функции определяются путем аппроксимации кривой. Такая функция называется трендом или функцией тренда, а группа типов молекул — гомологическим рядом. . Лловелл и др. (2013а, 2013б) [44] [45] предлагаемые функции тренда для трех параметров FV для алканов.Оливейра и др. (2014) [46] предложенные функции тренда для параметров FV для метиловых эфиров жирных кислот (FAME) и этиловых эфиров жирных кислот (FAEE), включая соединения с количеством ненасыщенных связей до трех, которые показаны ниже.
Молярная масса M [г/моль] (или молекулярная масса/вес), связанная с параметрами, используемыми в процессе аппроксимации кривой (где , , и являются эмпирическими параметрами) соответствует числам атомов углерода в диапазоне 8-24 и 8-20 для FAME и FAEE соответственно.
Теория значимой структуры
[ редактировать ]Модели вязкости, основанные на существенной теории структуры, это обозначение происходит от Айринг, [47] [48] (краткая теория СС и модель СС) в первые два десятилетия 2000-х годов превратились в эстафету развития. Это начинается с Масиас-Салинас и др. (2003), [49] продолжился при значительном вкладе со стороны Круз-Рейес и др. (2005), [50] за которым следует третий этап развития Масиас-Салинас и др. (2013), [51] чья модель здесь представлена. Теории СС имеют три основных предположения:
- Жидкость во многих аспектах ведет себя аналогично твердому телу, например, в чувствительной зависимости между молярным объемом (или массовой плотностью) и давлением; положение и расстояние между молекулами подобны квазирешетке с «ожиженными вакансиями» молекулярного размера, распределенными случайным образом по всей квазирешетке. Предполагается, что вакансии имеют молекулярный размер и свободно перемещаются по квазирешетчатой структуре.
- Вязкость жидкости рассчитывается из двух компонентов: газообразного и твердого, и оба вклада содержат все типы молекул, встречающиеся в жидкой фазе. Говорят, что молекула, которая перепрыгивает с одной поверхности скольжения на свободное место на соседней поверхности, ведет себя подобно газу. Говорят, что молекула, которая остается на своем месте на поверхности скольжения в течение некоторого времени, ведет себя как твердое тело.
- Столкновения между молекулами из соседних слоев эквивалентны прыжкам молекул на свободные места, и эти события в рамках моделирования вязкости аналогичны химическим реакциям между сталкивающимися молекулами в теории ТС.
Доля газообразных молекул и твердые молекулы являются
где – мольный объем рассматриваемой фазы, - молярный объем твердоподобных молекул и – молярный объем твердого ядра. Вязкость жидкости представляет собой смесь этих двух классов молекул.
Газоподобный вклад
[ редактировать ]Вклад газоподобной вязкости взят из модели вязкости Чунг и др. (1984, 1988), [52] [5] который основан на книге Чепмена-Энскога (1964).кинетическая теория вязкости для разбавленных газов и эмпирическое выражение Нойфельд и др. (1972) [6] для уменьшенного интеграла столкновений, но расширенного эмпирического для многоатомных, полярных и водородныхСвязующие жидкости в широком диапазоне температур. Модель вязкости Chung et al. (1988):
где
Список местной номенклатуры:
- : коэффициент молекулярной формы и полярности разбавленных газов [1]
- : молярная масса, конф. молекулярная масса [г/моль]
- : температура [К]
- : критическая температура [К]
- : молярный критический объем [см 3 /моль]
- : вклад газоподобной вязкости [мкП]
- : поправочный коэффициент для эффектов водородных связей [1]
- : уменьшенный дипольный момент [1]
- : уменьшенный интеграл столкновений [1]
- : ацентрический фактор [1]
Солидный вклад
[ редактировать ]В 2000-х годах разработка вклада твердоподобной вязкости началась с Масиас-Салинас и др. (2003) [49] который использовал уравнение Эйринга в теории TS как аналог вклада вязкости твердого тела и как обобщение первой модели экспоненциальной вязкости жидкости, предложенной Рейнольдсом ( 1886). [53] Уравнение Эйринга моделирует необратимые химические реакции при постоянном давлении, поэтому в уравнении используется Гиббса энергия активации : , чтобы смоделировать энергию переходного состояния, которую система использует для перемещения вещества (т. е. отдельных молекул) из начального состояния в конечное состояние (т. е. новое соединение). В течении Куэтта система перемещает вещество с одной поверхности скольжения на другую за счет флуктуации внутренней энергии, а также, вероятно, за счет давления и градиента давления. Кроме того, влияние давления на вязкость несколько иное для систем в диапазоне среднего давления, чем для систем в диапазоне очень высокого давления. Круз-Рейес и др. (2005) [50] использует энергию Гельмгольца (F = U-TS = G-PV) в качестве потенциала в экспоненциальной функции. Это дает
Круз-Рейес и др. (2005). [50] утверждает, что энергия активации Гиббса отрицательно пропорциональна внутренней энергии испарения (и, таким образом, рассчитывается в точке на кривой замерзания), но Масиас-Салинас и др. (2013) [51] меняет это на остаточную внутреннюю энергию, , при общем давлении и температуре системы. В качестве альтернативы можно использовать большой потенциал ( = U-TS-G = -PV, иногда называемая энергией или потенциалом Ландау) в экспоненциальной функции и утверждаем, что поток Куэтта не является однородной системой , поэтому необходимо добавить член с остаточной внутренней энергией. Оба аргумента дают предлагаемый твердый вклад, который
Предэкспоненциальный фактор принимается как
Частота прыжка молекулы, которая перепрыгивает из исходного положения на вакантное место, , ставится в зависимость от количества вакансий, и давление с целью расширения применимости в гораздо более широких диапазонах температуры и давления, чем это можно было бы сделать при постоянной скачкообразной частоте. Окончательная модель частоты скачков:
Постоянной проблемой для моделей вязкости является расчет молярного объема жидкости для заданного давления с использованием несовершенного EOS. Это требует введения некоторых эмпирических параметров. Использование регулируемых коэффициентов пропорциональности как для остаточной внутренней энергии, так и для Z-фактора является естественным выбором. Чувствительность значений P и Vb для жидкостей делает естественным введение эмпирического показателя степени (степени) в безразмерный Z-фактор. Эмпирическая сила оказывается очень эффективной в области высокого давления (высокого Z-фактора). Вклад твердоподобной вязкости, предложенный Масиасом-Салинасом и др. (2013). [51] тогда
Список местной номенклатуры:
- : молярный объем твердого ядра жидкой фазы [см 3 /моль]
- : давление [бар]
- : температура [К]
- : молярный объем жидкой фазы [см 3 /моль]
- : объемная доля j-подобного вклада j=g,s [1]
- : коэффициент сжимаемости (Z-фактор) [1]
- : коэффициент пропорциональности [1]
- : регулируемые параметры i=0,1 [1]
- : вязкость жидкой фазы [мкПа·с]
- : вклад твердоподобной вязкости [мкПа·с]
- : регулируемые параметры i=0,1 [с −1 ] и [бар −1 с −1 ]
- : энергия активации жидкости [Дж/моль]
- : остаточная внутренняя энергия жидкости [Дж/моль]
Смесь
[ редактировать ]Чтобы прояснить приведенные выше математические утверждения, вклад твердого тела для смеси жидкостей показан более подробно ниже.
Правила смешивания
[ редактировать ]Переменные а все параметры EOS для смеси жидкостей берутся из EOS (конф. W) и правил смешивания, используемых EOS (конф. Q). Более подробная информация об этом представлена ниже.
Жидкость с n молями в однофазной области, где общий состав жидкости равен [мольфракции]:
Газовая фаза, нг моль в двухфазной области, где состав газа [мольфракции]:
Жидкая фаза n л моль в двухфазной области, где состав жидкости [мольфракции]:
где
Поскольку почти все исходные данные для этой модели вязкости предоставляются EOS и расчетами равновесия, эту модель SS (или модель TS) для вязкости должно быть очень просто использовать для смесей жидкостей. Модель вязкости также имеет некоторые эмпирические параметры, которые можно использовать в качестве параметров настройки для компенсации несовершенства моделей EOS и обеспечения высокой точности также для смесей жидкостей.
См. также
[ редактировать ]- Теория столкновений
- Коэффициент сжимаемости
- Концентрация
- Преобразование единиц
- Подгонка кривой
- Банк данных Дортмунда
- Идеальный газ
- Закон идеального газа
- Кинетическая теория газов
- Список термодинамических свойств
- Молярный объем
- Обвинение
- настоящий газ
- Закон о реальном газе
- Статистическая механика
- Статистическая физика
- Таблица термодинамических уравнений
- Транспортные явления
- Оценка тренда
- ВТПР
Ссылки
[ редактировать ]- ^ Jump up to: а б Сирс, ФРВ; Сэлинджер, Г.Л. (1975). «10». Термодинамика, кинетическая теория и статистическая термодинамика (3-е изд.). Ридинг, Массачусетс, США: Addison-Wesley Publishing Company, Inc., стр. 286–291. ISBN 978-0201068948 .
- ^ Jump up to: а б с д и Педерсен, Канзас; Фреденслунд, Аа .; Томассен, П. (1989). Свойства нефти и природных газов . Том. 1989. Издательство Gulf Publishing Company, Хьюстон. стр. 1–252. ISBN 9780872015883 .
- ^ Jump up to: а б с д и ж г час я дж к л м н Зеберг-Миккельсен, СК (2001). «Исследование вязкости углеводородных жидкостей в пластовых условиях – моделирование и измерения». доктор философии Диссертация в Техническом университете Дании. Кафедра химической инженерии . Июнь (2001 г.): 1–271. ISBN 9788790142742 .
- ^ Jump up to: а б с д Уехара, ОА; Уотсон, К.М. (1944). «Универсальная корреляция вязкости». Национальные нефтяные новости . 39 (октябрь): Р-714–Р-722.
- ^ Jump up to: а б Чунг, Т.-Х.; Аджлан, М.; Ли, LL; Старлинг, К.Э. (1988). «Обобщенная многопараметрическая корреляция свойств переноса неполярных и полярных жидкостей». Индийский англ. хим. Рез . 27 (4): 671–679. дои : 10.1021/ie00076a024 .
- ^ Jump up to: а б Нойфельд, PD; Янзен, Арканзас; Азиз, Р.А. (1972). «Эмпирические уравнения для расчета 16 интегралов транспортных столкновений Ω (l, s) * для потенциала Леннарда-Джонса (12-6)». Журнал химической физики . 57 (3): 1100–1102. Бибкод : 1972ЖЧФ..57.1100Н . дои : 10.1063/1.1678363 .
- ^ Вилке, ЧР (1950). «Уравнение вязкости газовых смесей». Журнал химической физики . 18 (1950): 517–519. Бибкод : 1950JChPh..18..517W . дои : 10.1063/1.1747673 .
- ^ Jump up to: а б Хернинг, Ф.; Зипперер, Л. (1936). «Немецкий: Вклад в расчет вязкости технических газовых смесей по значениям вязкости отдельных компонентов; английский: Расчет вязкости технических газовых смесей по вязкости отдельных газов». Газо-водяной отсек . 79 (1936): 49–54 и 69–73.
- ^ Аррениус, С. (1887). «О внутреннем трении разбавленных водных растворов». З. Физ. Хим . 1 (1887): 2855–298.
- ^ Грюнберг, Л.; Ниссан, АХ (1949). «Закон смеси вязкости». Природа . 164 (1949): 799–800. Бибкод : 1949Natur.164..799G . дои : 10.1038/164799b0 . ПМИД 15395375 . S2CID 4096294 .
- ^ Катти, ПК; Чаудри, ММ (1964). «Вязкость бинарных смесей бензилацетата с диоксаном, анилином и м-крезолом» . Журнал химических и инженерных данных . 9 (1964): 442–443. дои : 10.1021/je60022a047 . hdl : 2027/mdp.39015086509596 .
- ^ Гласстоун, С.; Лейдлер, К.Дж.; Айринг, Х. (1941). Теория скоростных процессов, кинетика химических реакций, вязкости, диффузии и электрохимических явлений . Макгроу Хилл, Нью-Йорк.
- ^ Цванциг, Р. (1965). «Корреляционные функции времени и коэффициенты переноса в статистической механике». Ежегодный обзор физической химии . 16 (1965): 67–102. Бибкод : 1965ARPC...16...67Z . дои : 10.1146/annurev.pc.16.100165.000435 .
- ^ Джосси, Дж.А.; Стил, Л.И.; Тодос, Г. (1961). «Вязкость чистых веществ в плотных газообразных и жидких фазах». Журнал Айше . 8 (1962): 59–63. дои : 10.1002/aic.690080116 .
- ^ Лоренц Дж.; Брей, Б.Г.; Кларк, CR (1964). «Расчет вязкости пластовых жидкостей по их составам» . Журнал нефтяных технологий . Октябрь (1964): 1171–1176. дои : 10.2118/915-ПА .
- ^ Стил, Л.И.; Тодос, Г. (1961). «Вязкость неполярных газов при нормальном давлении». Айше Дж . 7 (1961): 611–615. дои : 10.1002/aic.690070416 .
- ^ Педерсен, Канзас; Фреденслунд, Аа .; Томассен, П. (1984). «Вязкость сырой нефти». хим. англ. Наука . 39 (1984): 1011–1016. Бибкод : 1984ЧЭнС..39.1011П . дои : 10.1016/0009-2509(84)87009-8 .
- ^ Jump up to: а б Педерсен, Канзас; Фреденслунд, Аа. (1987). «Улучшенная модель соответствующих состояний для прогнозирования вязкости и теплопроводности нефти и газа». хим. англ. Наука . 42 (1987): 182–186. Бибкод : 1987ЧЭнС..42..182П . дои : 10.1016/0009-2509(87)80225-7 .
- ^ Тэм, MJ; Габбинс, К.Э. (1970). «Принцип соответствия транспортных свойств плотных жидкостей неполярным многоатомным жидкостям». Индийский англ. хим. Фундам . 9 (1975): 63–70. дои : 10.1021/i160033a010 .
- ^ Бенедикт, В.; Уэбб, Великобритания; Рубин, Л.К. (1940). «Эмпирическое уравнение термодинамических свойств легких углеводородов и их смесей. I. Метан, этан. Пропан и бутан». Дж. Хим. Физ . 8 (4): 334–345. Бибкод : 1940ЖЧФ...8..334Б . дои : 10.1063/1.1750658 .
- ^ Маккарти, Р.Д. (1974). «Модифицированное уравнение состояния Бенедикта-Уэбба-Рубина для метана с использованием последних экспериментальных данных». Криогеника . 14 (5): 276–280. Бибкод : 1974Крио...14..276М . дои : 10.1016/0011-2275(74)90228-8 .
- ^ Хэнли, HJM; Маккарти, доктор медицинских наук; Хейнс, WM (1975). «Уравнение для коэффициентов вязкости и теплопроводности метана». Криогеника . 15 (1975): 413–417. Бибкод : 1975Крио...15..413H . дои : 10.1016/0011-2275(75)90010-7 .
- ^ Филлипс, П. (1912). «Вязкость углекислого газа». Труды Лондонского королевского общества . 87А (1912): 48–61. Бибкод : 1912РСПСА..87...48П . дои : 10.1098/rspa.1912.0058 . ISSN 0950-1207 .
- ^ Литтл, Дж. Э.; Кеннеди, ХТ (1968). «Корреляция вязкости углеводородных систем с давлением, температурой и составом» . Журнал Общества инженеров-нефтяников . 8 июня (2): 157–162. дои : 10.2118/1589-ПА .
- ^ Го, X.-Q.; Ван, Л.-С.; Ронг, С.-Х.; Го, Т.-М. (1968). «Модель вязкости, основанная на уравнениях состояния углеводородных жидкостей и газов». Журнал Общества инженеров-нефтяников . 139 (1997): 405–421. дои : 10.1016/S0378-3812(97)00156-8 .
- ^ Патель, Северная Каролина; Теджа, А.С. (1982). «Новое кубическое уравнение состояния жидкостей и смесей жидкостей». Химико-техническая наука . 37 (1982): 463–473. Бибкод : 1982ЧЭнС..37..463П . дои : 10.1016/0009-2509(82)80099-7 .
- ^ Го, X.-Q. (2001). «Частные сообщения с К. К. Зеберг-Миккельсеном». доктор философии Диссертация в Техническом университете Дании. Кафедра химической инженерии . Июнь (2001 г.): 1–271. ISBN 9788790142742 .
- ^ Jump up to: а б Киньонес-Сиснерос, ЮВ; Зеберг-Миккельсен, СК; Стенби, Э.Х. (2000). «Теория трения (f-теория) для моделирования вязкости». Жидкостно-фазовые равновесия . 169 (2000): 249–276. дои : 10.1016/S0378-3812(00)00310-1 .
- ^ Jump up to: а б Киньонес-Сиснерос, ЮВ; Зеберг-Миккельсен, СК; Стенби, Э.Х. (2001a). «Однопараметрические модели теории трения для вязкости». Жидкостно-фазовые равновесия . 178 (2001а): 1–16. дои : 10.1016/S0378-3812(00)00474-X .
- ^ Jump up to: а б Киньонес-Сиснерос, ЮВ; Зеберг-Миккельсен, СК; Стенби, Э.Х. (2001b). «Теория трения для моделирования вязкости: распространение на системы сырой нефти». Жидкостно-фазовые равновесия . 56 (2001б): 7007–7015. Бибкод : 2001ЧЭнС..56.7007Q . дои : 10.1016/S0009-2509(01)00335-9 .
- ^ Jump up to: а б Киньонес-Сиснерос, ЮВ; Дальберг, А.; Стенби, Э.Х. (2004). «PVT-характеристика, моделирование вязкости и прогноз сырой нефти». Журнал «Нефтяная наука и технология» . 22 (9–10): 1309–1325. дои : 10.1081/LFT-200034092 . S2CID 96591609 .
- ^ Киньонес-Сиснерос, ЮВ; Дейтерс, Великобритания (2006). «Обобщение теории трения для моделирования вязкости». Журнал физической химии Б. 110 (25): 12820–12834. дои : 10.1021/jp0618577 . ПМИД 16800618 .
- ^ Дулитл, АК (1951). «Исследования ньютоновского течения. II - Зависимость вязкости жидкостей от свободного пространства». Журнал прикладной физики . 22 (12): 1471–1475. Бибкод : 1951JAP....22.1471D . дои : 10.1063/1.1699894 .
- ^ Аллал, А.; Монтфорд, Япония; Марин, Г. (1996). «Молекулярная реология: расчет вязкоупругих свойств по микроструктуре полимеров». Материалы XII Международного конгресса по реологии под редакцией Айта Кади А., Дили Дж. М., Джеймса Д. Ф. и Уильямса М. К. из канадской группы реологии . 317 . ISBN 9782980510908 .
- ^ Аллал, А.; Моха-Ушан, М.; Бонед, К. (2001a). «Новая модель свободного объема для динамической вязкости и плотности плотных жидкостей в зависимости от давления и температуры». Физика и химия жидкостей . 39 : 1–30. дои : 10.1080/00319100108030323 . S2CID 95517075 .
- ^ Jump up to: а б с д Аллал, А.; Бонед, К.; Байлаук, А. (2001b). «Модель вязкости в свободном объеме для жидкостей в плотном и газообразном состояниях». Физ. Преподобный Е. 64 (1): 1203–. Бибкод : 2001PhRvE..64a1203A . дои : 10.1103/PhysRevE.64.011203 . ПМИД 11461236 .
- ^ Jump up to: а б Даллиен, ФАЛ (1963). «Новая связь между вязкостью и коэффициентами диффузии, основанная на теории диффузии Ламма». Пер. Фарадей Соц . 59 : 856–868. дои : 10.1039/TF9635900856 .
- ^ Jump up to: а б Бонед, К.; Аллал, А.; Байлаук, А.; Зеберг-Миккельсен, СК; Бессьер, Д.; Киньонес-Сиснерос, SE (2004). «Одновременное моделирование коэффициента самодиффузии и динамической вязкости при высоком давлении в свободном объеме» (PDF) . Физ. Преподобный Е. 69 (3): 1–6. Бибкод : 2004PhRvE..69c1203B . дои : 10.1103/PhysRevE.69.031203 . ПМИД 15089276 . S2CID 9314966 .
- ^ Даллиен, ФАЛ (1972). «Прогнозирующие уравнения самодиффузии в жидкостях: другой подход». Айше Дж . 18 (1): 62–70. дои : 10.1002/aic.690180113 .
- ^ Чунг1988
- ^ Аллал, А.; Бонед, К.; Дауге, П. (2001c). «Новая модель свободного объема для определения динамической вязкости плотных жидкостей в зависимости от давления и температуры. Расширение модели прогнозирования для не очень ассоциативных смесей». Физика и химия жидкостей . 39 (5): 607–624. дои : 10.1080/00319100108030681 . S2CID 98468228 .
- ^ Кане, X. (2001). «Динамическая вязкость и плотность при высоких давлениях бинарных и тройных смесей тяжелых и легких углеводородов». Докторская диссертация, Университет По, По, Франция .
- ^ Алмаси, М. (2015). «Температурная зависимость и влияние длины цепи на плотность и вязкость бинарных смесей нитробензола и 2-спиртов». Журнал молекулярных жидкостей . 209 : 346–351. дои : 10.1016/j.molliq.2015.05.045 .
- ^ Лловелл, Ф.; Маркос, РМ; Вега, Л.Ф. (2013а). «Теория свободного объема в сочетании с soft-SAFT для расчета вязкости: сравнение с данными молекулярного моделирования и экспериментов». Дж. Физ. хим. Б. 117 (27): 8159–8171. дои : 10.1021/jp401307t . ПМИД 23789584 .
- ^ Лловелл, Ф.; Маркос, РМ; Вега, Л.Ф. (2013б). «Транспортные свойства смесей по теории мягкого SAFT + свободного объема: применение к смесям н-алканов и гидрофторуглеродов». Дж. Физ. хим. Б. 117 (17): 5195–5205. дои : 10.1021/jp401754r . ПМИД 23566079 .
- ^ Оливейра, МБ; Фрейтас, СВД; Лловелл, Ф.; Вега, ЛФ; Коутиньо, Япония (2014). «Разработка простых и переносимых молекулярных моделей производства биодизеля с использованием уравнения состояния soft-SAFT». Химико-технологические исследования и проектирование . 92 (12): 2898–2911. дои : 10.1016/j.cherd.2014.02.025 .
- ^ Айринг, Х.; Ри, Т.; Хираи, Н. (1958). «Значительные структуры в жидком состоянии» . Учеб. Натл. акад. наук. США . 44 (7): 683–691. Бибкод : 1958PNAS...44..683E . дои : 10.1073/pnas.44.7.683 . ПМК 528643 . ПМИД 16590259 .
- ^ Хендерсон, Д. (2010). «Генри Айринг: квантовая химия, статистическая механика, теория жидкостей и теория существенной структуры». Бык. Хист. Хим . 35 (2).
- ^ Jump up to: а б Масиас-Салинас, Р.; Гарсиа-Санчес, Ф.; Эрнандес Гардуса, О. (2003). «Модель вязкости для чистых жидкостей на основе теории Айринга и кубического EoS». Айше Дж . 49 (3): 799–804. дои : 10.1002/aic.690490324 .
- ^ Jump up to: а б с Круз-Рейес, Г.; Луна-Барсенас, Г.; Альварадо, JFJ; Санчес, IC; Масиас-Салинас, Р. (2005). «Одновременная корреляция вязкостей насыщенных газов и жидкостей с использованием теории значимой структуры». Индийский англ. хим. Рез . 44 (6): 1960–1966. дои : 10.1021/ie049070v .
- ^ Jump up to: а б с Масиас-Салинас, Р.; Акино-Оливос, Массачусетс; Гарсиа-Санчес, Ф. (2013). «Моделирование вязкости пластовых жидкостей в широком диапазоне температур и давлений». Химико-технологические операции . 32 : 1573. дои : 10.3303/CET1332263 . ISBN 978-88-95608-23-5 . ISSN 1974-9791 .
- ^ Чунг, Т.-Х.; Ли, LL; Старлинг, К.Э. (1984). «Применение кинетических теорий газа и многопараметрической корреляции для прогнозирования вязкости и теплопроводности разбавленного газа». Основы промышленной и инженерной химии . 23 (1): 8–13. дои : 10.1021/i100013a002 .
- ^ Рейнольдс, О. (1886). «О теории смазки и ее применении к экспериментам г-на Бошампа Тауэра, включая экспериментальное определение вязкости оливкового масла». Фил. Пер. Р. Сок. Лонд . 177 : 157–234. дои : 10.1098/rstl.1886.0005 . S2CID 110829869 .