Теория моделей кавитации

В уравнении 2.169 переноса пара Re и Rc - это характеристики источника переноса массы, связанные с ростом и схлопыванием пузырьков пара в кавитационных потоках. Эти условия и переменные учитываются для расчета обмена масс между паровой и жидкой фазами во время процесса кавитации. Re и Rc моделируются в уравнении Рэлея - Плессета (Rayleigh-Plesset), описывающем рост одного пузырька пара в жидкости.

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

В двухфазной системе "жидкость-пар" плотность смеси ρ выражается через объемную долю пара и плотности фаз:

уравнение 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

Это упрощенное уравнение Релея - Плессета описывает физический подход, позволяющий ввести эффекты динамики пузырьков в модели кавитации. Радиус пузырьков может увеличиваться или уменьшаться в зависимости от знака выражений (pB - p): пузырек растет, если p < pB, и сжимается, если p > pB. Поэтому уравнение 2.188 переписывается следующим образом:

уравнение 2.189

Если η - это плотность числа пузырьков в жидкости (число пузырьков на единицу объема) и все пузырьки пара являются совершенными сферами с одинаковым радиусом RB, то объемная доля фазы пара имеет вид:

уравнение 2.190

Предполагается, что пузырьки пара не могут создаваться или разрушаться в жидкости, но они могут расти (испарение) и схлопываться (конденсация) во время процесса кавитации. В уравнении 2.190 плотность числа пузырьков пара (η) остается постоянной, но радиус пузырьков (RB) увеличивается или уменьшается. Тогда производная по времени доли объема пара может рассчитываться следующим образом:

уравнение 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 массовой доли пара, получаем следующие характеристики источника Re и Rc:

Уравнение 2.194 и уравнение 2.195 формируют основу почти всех доступных механических двухфазных моделей кавитации. В Creo Flow Analysis использован метод моделирования Сингхала и др..

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

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

fg,f	массовые доли свободного газа
fg,d	массовые доли свободного газа и растворенного газа,
Sg,f, Sg,d	внешние или определяемые пользователем источники.
Dg,f	диффузионная способность свободного газа и растворенного газа
Dg,d	диффузионная способность растворенного газа

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

AI	область на границе раздела жидкости-газа
kf,d(=kd,f)	объемный коэффициент волюметрического переноса массы
ρg,d(=ρfg,d)	локальные концентрации массы растворенного газа
ρg,f(=ρfg,f)	локальные концентрации массы свободного газа
	равновесные концентрации массы растворенного газа в содержащих его фазах
	равновесные концентрации массы свободного газа в содержащих его фазах

Обычно

не совпадают (нарушение непрерывности). Между этими двумя концентрациями существует четко определенная кривая равновесия, которая зависит от температуры, давления и состава смеси. Эта кривая, как правило, является монотонной и нелинейной и часто выражается как квазилинейное уравнение с коэффициентом

где Kf,d обычно решается с помощью физических законов или эмпирических корреляций. Один общий метод заключается в выполнении закона Генри, обеспечивающего обобщенную равновесную взаимосвязь. Он утверждает, что для жидкой смеси, соприкасающейся с газовой фазой, частичное давление свободного газа ρg,f равно произведению равновесной мольной доли растворенного газа в жидкой фазе

и константы Генри, Hx:

Таким образом, константа Генри имеет единицу измерения давления и может рассматриваться в качестве эталонного давления. Если предположить, что для идеальной жидкой смеси в контакте с газом константа Генри равна давлению насыщенного пара psat, то уравнение 2.207 также имеет следующий вид:

В уравнении 2.202 и уравнении 2.203 переменные

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

◦ Поглощение/растворение газа в жидкости - имеет место перенос массы из газовой фазы (свободный газ) в жидкую фазу (растворенный газ). Модель равновесия предполагает, что свободный газ в газовой фазе находится в состоянии равновесия:

. Используя уравнение 2.204 и уравнение 2.209, получаем следующий результат:

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

. Из уравнения 2.202, уравнения 2.203, уравнения 2.204 и уравнения 2.209 скорости переноса массы имеют следующий вид:

Отметим, что при растворении части неконденсирующегося газа в жидкости свободно расширяющийся газ является единственной частью, остающейся в газовой фазе, fg,f. Следовательно, плотность смеси вычисляется следующим образом: