Теория моделей кавитации
В
уравнении 2.169 переноса пара R
e и R
c - это характеристики источника переноса массы, связанные с ростом и схлопыванием пузырьков пара в кавитационных потоках. Эти условия и переменные учитываются для расчета обмена масс между паровой и жидкой фазами во время процесса кавитации. R
e и R
c моделируются в уравнении
Рэлея - Плессета (Rayleigh-Plesset), описывающем рост одного пузырька пара в жидкости.
Перенос массы "жидкость-пар"
Чтобы получить выражение для результирующей скорости фазового перехода в кавитации, рассматривайте несжимаемый жидкий пар как два фазовых потока с нулевой скоростью скольжения, где неконденсирующийся газ не рассматривается. Если ввести R для представления результирующей скорости переноса массы из жидкости в пар, то уравнения объемных долей жидкости и пара и уравнение непрерывности полной массы можно представить следующим образом:
• Жидкая фаза
уравнение 2.179
• Фаза пара
уравнение 2.180
• Смесь (непрерывность общей массы)
уравнение 2.181
В двухфазной системе "жидкость-пар" плотность смеси ρ выражается через объемную долю пара и плотности фаз:
уравнение 2.182
Поскольку предполагается, что обе плотности - жидкости и пара - постоянны (несжимаемы), отношение между градиентом скорости и долей объема пара выводится из
уравнения 2.181 и
уравнения 2.182:
уравнение 2.183
Комбинация уравнений
2.179 и
уравнение 2.183 дает следующее выражение для характеристики источника общей массы:
уравнение 2.184
Если подставить
уравнение 2.184 в
уравнение 2.180, уравнение для объемной доли пара можно переписать в общем виде:
уравнение 2.185
Если применить соотношение между массовой долей и долей объема из
уравнения 2.173 для пара, получится
уравнение 2.185 в терминах массовой доли пара:
уравнение 2.186
Из
уравнения 2.185 и
уравнения 2.186 ясно, что при условии нулевой скорости скольжения между фазами жидкости и пара кавитацию можно моделировать как однофазный поток с дополнительным транспортным уравнением для массовой доли пара или как
эйлеровский поток многофазной смеси с переносом массы между фазами жидкости и пара. Если не учитывать эффект диффузии и различия в фазовой скорости, эти два подхода математически идентичны.
Creo Flow Analysis использует однофазный метод для моделирования кавитационных потоков.
Описание динамики пузырьков
Для большинства природных ситуаций и технических систем в жидкости для возникновения кавитации должно существовать соответствующее число ядер, таких как пузырьки, неконденсируемые газы и т. д. Соответственно для моделирования процесса кавитации необходим, прежде всего, правильный расчет роста и схлопывания пузырьков. Если предположить, что в жидкости, движущейся в потоке, существует нулевая скорость скольжения между пузырьками жидкости и пара, уравнение динамики пузырьков выводится из обобщенного уравнения Релея - Плессета, описывающего рост газового пузырька в жидкости:
уравнение 2.187
где
RB | радиус пузырька |
pB | давление в пузырьке (предполагается, что это давление пара при данной температуре жидкости в отсутствие других газов) |
p | давление в жидкости, окружающей пузырек |
σ | коэффициент поверхностного натяжения между жидкостью и паром |
Это уравнение получается из механического баланса (нет тепловых барьеров для роста пузырьков).
Если пренебречь производной по времени второго порядка (для низкочастотных колебаний), членом вязкого трения и силой поверхностного натяжения, мы получим упрощенную форму
уравнения 2.187, которая действительна для асимптотического состояния:
уравнение 2.188
Это упрощенное уравнение
Релея - Плессета описывает физический подход, позволяющий ввести эффекты динамики пузырьков в модели кавитации. Радиус пузырьков может увеличиваться или уменьшаться в зависимости от знака выражений (p
B - p): пузырек растет, если p < p
B, и сжимается, если p > p
B. Поэтому
уравнение 2.188 переписывается следующим образом:
уравнение 2.189
Если η - это плотность числа пузырьков в жидкости (число пузырьков на единицу объема) и все пузырьки пара являются совершенными сферами с одинаковым радиусом RB, то объемная доля фазы пара имеет вид:
уравнение 2.190
Предполагается, что пузырьки пара не могут создаваться или разрушаться в жидкости, но они могут расти (испарение) и схлопываться (конденсация) во время процесса кавитации. В
уравнении 2.190 плотность числа пузырьков пара (η) остается постоянной, но радиус пузырьков (R
B) увеличивается или уменьшается. Тогда производная по времени доли объема пара может рассчитываться следующим образом:
уравнение 2.191
Подстановка
уравнения 2.191 в
уравнение 2.186 дает транспортное уравнение, управляющее массовой долей паровой фазы:
уравнение 2.192
Применение
уравнения 2.189 дает для результирующей скорости переноса массы на единицу объема между жидкостью и паром следующий вид:
уравнение 2.193
где
уравнение 2.193 показывает, что в кавитации скорость переноса массы в единице объема (R) является функцией (пропорциональна) плотности паровой и жидкой фаз, а также обратно пропорциональна плотности смеси. Поскольку
уравнение 2.192 получено непосредственно из непрерывности фазы и массы смеси, оно является точным и должно точно представлять перенос массы между фазами жидкости и пара в кавитации. С введением динамики пузырьков
уравнение 2.193 использует аналогичный метод для моделирования двух противоположных (и физически различающихся) процессов переноса массы: из жидкости в пар (рост или испарение пузырьков) и из пара в жидкость (схлопывание или конденсация пузырьков). Для транспортного
уравнения 2.192 массовой доли пара рост пузырьков является характеристикой источника, в то время как схлопывание пузырьков обрабатывается как характеристика теплопоглощения.
В практических моделях кавитации локальное давление далекого поля p обычно считается совпадающим с давлением в центре ячейки. Давление пузырьков pB равно давлению насыщенного пара (psat, свойство материала) в отсутствие растворенных газов, движения масс и вязкого трения, pB = psat
Сравнивая
уравнение 2.192 и
уравнение 2.193 с общим уравнением 2.169 массовой доли пара, получаем следующие характеристики источника R
e и R
c:
• Если локальное давление потока ниже давления насыщенного пара p<psat, возникает только испарение, так что:
уравнение 2.194
• Если локальное давление потока выше давления насыщенного пара p>psat, существует только конденсация:
уравнение 2.195
Уравнение 2.194 и
уравнение 2.195 формируют основу почти всех доступных механических двухфазных моделей кавитации. В
Creo Flow Analysis использован метод моделирования
Сингхала и др..
Поглощение/растворение и выделение газа
Неконденсирующиеся газы часто присутствуют в рабочей жидкости и могут оказывать значительное влияние на кавитацию. Иногда неконденсирующийся газ свободно перемещается с потоком, а также растворяется в жидкости или выделяется из жидкости, естественным образом стремящейся к динамическому равновесию концентраций масс между фазами жидкости и газа. Поглощение или растворение газа и выделение его в жидкости - это также явления переноса массы жидкости-газа, управляемое перепадами и градиентами концентраций массы. Чтобы моделировать кавитационные потоки, необходимо также учитывать воздействие неконденсирующегося газа и возможный перенос массы между жидкостью и газом в потоке смеси.
Если предположить, что в двухфазовом потоке "жидкость-газ" неконденсируемый газ, такой как воздух или кислород, существует и в фазе жидкости (растворенный газ), и в фазе газа (свободный газ), транспортные уравнения для массовой доли газа в каждой фазе будут следующими:
• Свободный газ (фаза газа)
уравнение 2.196
• Растворенный газ (фаза жидкости)
уравнение 2.197
где
fg,f | массовые доли свободного газа |
fg,d | массовые доли свободного газа и растворенного газа, |
Sg,f, Sg,d | внешние или определяемые пользователем источники. |
Dg,f | диффузионная способность свободного газа и растворенного газа |
Dg,d | диффузионная способность растворенного газа |
Если массовая доля неконденсирующегося газа предварительно описана как fg, получаем:
fg=fg,f
или его перемещение в пространстве и времени получается из решения следующего уравнения:
уравнение 2.199
В
уравнении 2.196 и
уравнении 2.197 источник R
d→f указывает уровень выделения растворенного газа, а R
f→d указывает уровень поглощения или растворения свободного газа.
Если эти две фазы соприкасаются, наблюдается тенденция к перемещению свободного газа f и растворенного газа d из одной фазы в другую, чтобы достигнуть динамического равновесия между этими двумя фазами. Модели равновесия предполагают, что волюметрические скорости переносов массы зависят от градиентов или различий в концентрациях масс:
уравнение 2.200
уравнение 2.201
где
AI | область на границе раздела жидкости-газа |
kf,d(=kd,f) | объемный коэффициент волюметрического переноса массы |
ρg,d(=ρfg,d) | локальные концентрации массы растворенного газа |
ρg,f(=ρfg,f) | локальные концентрации массы свободного газа |
| равновесные концентрации массы растворенного газа в содержащих его фазах |
| равновесные концентрации массы свободного газа в содержащих его фазах |
Отметим, что в
![](../../../simulate/cfd/images/cfd-rds.5.57.25.jpg)
используется единица измерения 1/с как показатель эффективности переноса массы. Таким образом,
уравнение 2.200 и
уравнение 2.201 также имеют следующий вид:
уравнение 2.202
уравнение 2.203
Обычно
![](../../../simulate/cfd/images/cfd-rds.5.57.28.jpg)
и
![](../../../simulate/cfd/images/cfd-rds.5.57.29.jpg)
не совпадают (нарушение непрерывности). Между этими двумя концентрациями существует четко определенная кривая равновесия, которая зависит от температуры, давления и состава смеси. Эта кривая, как правило, является монотонной и нелинейной и часто выражается как квазилинейное уравнение с коэффициентом
![](../../../simulate/cfd/images/cfd-rds.5.57.30.jpg)
уравнение 2.204
где K
f,d обычно решается с помощью физических законов или эмпирических корреляций. Один общий метод заключается в выполнении закона
Генри, обеспечивающего обобщенную равновесную взаимосвязь. Он утверждает, что для жидкой смеси, соприкасающейся с газовой фазой, частичное давление свободного газа ρ
g,f равно произведению равновесной мольной доли растворенного газа в жидкой фазе
![](../../../simulate/cfd/images/cfd-rds.5.57.32.jpg)
и константы Генри, H
x:
уравнение 2.205
Если для фазы свободного газа выполняется закон идеального газа, то закон Дальтона для частичного давления дает следующее уравнение:
уравнение 2.206
уравнение 2.207
Таким образом, константа Генри имеет единицу измерения давления и может рассматриваться в качестве эталонного давления. Если предположить, что для идеальной жидкой смеси в контакте с газом константа Генри равна давлению насыщенного пара p
sat, то
уравнение 2.207 также имеет следующий вид:
уравнение 2.208
где p
d,equil,ref - ссылочное давление для растворенной равновесной массовой доли. Тогда
уравнение 2.204 можно переписать следующим образом:
уравнение 2.209
В
уравнении 2.202 и
уравнении 2.203 переменные
![](../../../simulate/cfd/images/cfd-rds.5.57.38.jpg)
и
![](../../../simulate/cfd/images/cfd-rds.5.57.39.jpg)
представляют различия/градиенты концентрации масс, которые определяют поглощение/растворение и выделение неконденсирующихся газов. Это показывает, что для переноса неконденсирующегося газа через две фазы требуется отклонение концентрации массы от ее равновесного состояния. Направление переноса массы может способствовать приближению системы к состоянию равновесия в зависимости от локальных концентраций масс и их равновесия в обеих фазах. Из модели равновесия можно получить следующий результат:
◦ Поглощение/растворение газа в жидкости - имеет место перенос массы из газовой фазы (свободный газ) в жидкую фазу (растворенный газ). Модель равновесия предполагает, что свободный газ в газовой фазе находится в состоянии равновесия:
![](../../../simulate/cfd/images/cfd-rds.5.57.40.jpg)
. Используя
уравнение 2.204 и
уравнение 2.209, получаем следующий результат:
уравнение 2.210
уравнение 2.211
уравнение 2.212
◦ Выделение газа из жидкости - имеет место перенос массы из жидкой фазы (растворенный газ) в газовую фазу (свободный газ). В этом процессе модель равновесия предполагает, что растворенный газ в жидкой фазе всегда находится в состоянии равновесия:
![](../../../simulate/cfd/images/cfd-rds.5.57.44.jpg)
. Из
уравнения 2.202,
уравнения 2.203,
уравнения 2.204 и
уравнения 2.209 скорости переноса массы имеют следующий вид:
уравнение 2.213
уравнение 2.214
Отметим, что при растворении части неконденсирующегося газа в жидкости свободно расширяющийся газ является единственной частью, остающейся в газовой фазе, fg,f. Следовательно, плотность смеси вычисляется следующим образом:
уравнение 2.215
и объемная доля свободного газа имеет следующий вид:
уравнение 2.216