Уравнение Рэлея – Плессета

редактировать
Уравнение Рэлея – Плессета часто применяется для изучения кавитации пузырьков, которые здесь образуются позади

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

R d 2 R dt 2 + 3 2 (d R dt) 2 + 4 ν LR d R dt + 2 γ ρ LR + Δ P (t) ρ L = 0 {\ displaystyle R {\ frac {d ^ {2} R} {dt ^ {2}}} + {\ frac {3} {2}} \ left ({\ frac {dR} {dt}} \ right) ^ {2 } + {\ frac {4 \ nu _ {L}} {R}} {\ frac {dR} {dt}} + {\ frac {2 \ gamma} {\ rho _ {L} R}} + {\ frac {\ Delta P (t)} {\ rho _ {L}}} = 0}{\ displaystyle R {\ frac {d ^ {2} R} {dt ^ {2}}} + {\ frac {3} {2}} \ left ({\ frac {dR} {dt}} \ справа) ^ {2} + {\ frac {4 \ nu _ {L}} {R}} {\ frac {dR} {dt}} + {\ frac {2 \ gamma} {\ rho _ {L} R }} + {\ frac {\ Delta P (t)} {\ rho _ {L}}} = 0}

, где

ρ L {\ displaystyle \ rho _ {L}}\ rho _ {L} - плотность окружающей жидкости, принимаемая постоянной
R (t) {\ displaystyle R (t)}R (t) - радиус пузыря
ν L {\ displaystyle \ nu _ {L}}\nu_L- кинематическая вязкость окружающей жидкости, принимаемая постоянной
γ {\ displaystyle \ gamma}\ gamma - поверхностное натяжение границы раздела пузырь-жидкость
Δ P (t) = P ∞ (t) - PB (t) {\ displaystyle \ Delta P (t) = P _ {\ infty} (t) -P_ {B} (t)}{\ displaystyle \ Delta P (t) = P _ {\ infty} (t) -P_ {B} (t)} , в котором PB (t) {\ displaystyle P_ {B} (t)}{\ displaystyle P_ {B} (t)} - давление внутри пузыря, считается однородным и P ∞ (t) {\ displaystyle P _ {\ infty} (t)}P _ {\ infty} (t) - это внешнее давление бесконечно далеко от пузыря

при условии, что PB (t) {\ displaystyle P_ {B} (t) }P_ {B} (t) известен и P ∞ (t) {\ displaystyle P _ {\ infty} (t)}P _ {\ infty} (t) задан, уравнение Рэлея – Плессета можно использовать для решения изменяющийся во времени радиус пузырька R (t) {\ displaystyle R (t)}R (t) .

Уравнение Рэлея – Плессета выводится из уравнений Навье – Стокса в предположении сферическая симметрия.

Содержание
  • 1 История
  • 2 Вывод
    • 2.1 Сохранение массы
    • 2.2 Сохранение импульса
    • 2.3 Граничные условия
  • 3 Решения
  • 4 Ссылки
История

Без учета поверхностного натяжения и вязкости уравнение было впервые выведено У. Г. Безант в своей книге 1859 г. постановка задачи сформулирована так: Бесконечная масса однородной несжимаемой жидкости, на которую не действуют никакие силы, находится в состоянии покоя, а сферическая часть жидкости внезапно аннигилирует; требуется найти мгновенное изменение давления в любой точке массы и время, в течение которого полость будет заполнена, при этом давление на бесконечном расстоянии должно оставаться постоянным (фактически, Безант приписывает проблему Кембриджу Сенатские проблемы 1847 г.). Пренебрегая колебаниями давления внутри пузырька, Безант предсказал, что время, необходимое для заполнения полости, составит

t = a (6 ρ p ∞) 1/2 ∫ 0 1 z 4 dz 1 - z 6 = a (π ρ 6 п ∞) 1/2 Γ (5/6) Γ (4/3) ≈ 0,91468 a (ρ p ∞) 1/2 {\ displaystyle {\ begin {align} t = a \ left ({\ frac {6 \ rho} {p _ {\ infty}}} \ right) ^ {1/2} \ int _ {0} ^ {1} {\ frac {z ^ {4} \, dz} {\ sqrt {1-z ^ {6}}}} \\ = a \ left ({\ frac {\ pi \ rho} {6p _ {\ infty}}} \ right) ^ {1/2} {\ frac {\ Gamma (5/6)} {\ Gamma (4/3)}} \\ \ приблизительно 0,91468a \ left ({\ frac {\ rho} {p _ {\ infty}}} \ right) ^ {1/2} \ end {выровнено }}}{\ displaystyle {\ begin {al igned} t = a \ left ({\ frac {6 \ rho} {p _ {\ infty}}} \ right) ^ {1/2} \ int _ {0} ^ {1} {\ frac {z ^ { 4} \, dz} {\ sqrt {1-z ^ {6}}}} \\ = a \ left ({\ frac {\ pi \ rho} {6p _ {\ infty}}} \ right) ^ { 1/2} {\ frac {\ Gamma (5/6)} {\ Gamma (4/3)}} \\ \ приблизительно 0,91468a \ left ({\ frac {\ rho} {p _ {\ infty}} } \ right) ^ {1/2} \ end {align}}}

где интегрирование было выполнено лордом Рэлеем в 1917 году, который вывел уравнение из баланса энергии. Рэлей также осознал, что предположение о постоянном давлении внутри полости станет неверным по мере уменьшения радиуса, и он показывает, используя закон Бойля, если радиус полости уменьшается в 4 1/3 раз. {\ displaystyle 4 ^ {1/3}}{\ displaystyle 4 ^ {1/3}} , то давление у границы полости становится больше, чем давление окружающей среды. Уравнение было впервые применено к движущимся кавитационным пузырькам Милтоном С. Плессетом в 1949 году путем включения эффектов поверхностного натяжения.

Вывод
Численное интегрирование RP eq.. включая поверхностное натяжение и вязкость. Первоначально в состоянии покоя при атмосферном давлении с R0 = 50 мкм пузырек, подвергающийся колебательному давлению на его собственной частоте, расширяется, а затем схлопывается. Численное интегрирование RP eq. включая поверхностное натяжение и вязкость. Первоначально находящийся в состоянии покоя при атмосферном давлении с R0 = 50 мкм пузырь, подверженный перепаду давления, расширяется, а затем схлопывается.

Уравнение Рэлея – Плессета может быть полностью выведено из первых принципов с использованием радиуса пузырька как динамический параметр. Рассмотрим сферический пузырь с зависящим от времени радиусом R (t) {\ displaystyle R (t)}R (t) , где t {\ displaystyle t}t - время. Предположим, что пузырек содержит однородно распределенный пар / газ с постоянной температурой TB (t) {\ displaystyle T_ {B} (t)}T_{B}(t)и давлением PB (t) {\ стиль отображения P_ {B} (t)}P_ {B} (t) . За пределами пузыря находится бесконечная область жидкости с постоянной плотностью ρ L {\ displaystyle \ rho _ {L}}\ rho _ {L} и динамической вязкостью мкл {\ displaystyle \ mu _ {L}}\ mu _ {L} . Пусть температура и давление вдали от пузыря будут T ∞ {\ displaystyle T _ {\ infty}}T _ {\ infty} и P ∞ (t) {\ displaystyle P _ {\ infty} (t) }P _ {\ infty} (t) . Температура T ∞ {\ displaystyle T _ {\ infty}}T _ {\ infty} считается постоянной. На радиальном расстоянии r {\ displaystyle r}r от центра пузыря меняющиеся свойства жидкости: давление P (r, t) {\ displaystyle P (r, t) }P (r, t) , температура T (r, t) {\ displaystyle T (r, t)}T(r,t)и скорость в радиальном направлении наружу u (r, t) {\ displaystyle u (r, t)}u (r, t) . Обратите внимание, что эти свойства жидкости определяются только за пределами пузырька, для r ≥ R (t) {\ displaystyle r \ geq R (t)}r \ geq R (t) .

Сохранение массы

По сохранению массы, закон обратных квадратов требует, чтобы скорость в радиальном направлении наружу u (r, t) {\ displaystyle u (r, t)}u (r, t) была обратно пропорциональна пропорционально квадрату расстояния от начала координат (центра пузыря). Следовательно, позволяя F (t) {\ displaystyle F (t)}F(t)быть некоторой функцией времени,

u (r, t) = F (t) r 2 {\ displaystyle u (r, t) = {\ frac {F (t)} {r ^ {2}}}}u (r, t) = {\ frac {F (t)} {r ^ {2}}}

В случае нулевого переноса массы через поверхность пузыря скорость на границе раздела должна быть

u (Р, t) знак равно d р dt знак равно F (t) р 2 {\ displaystyle u (R, t) = {\ frac {dR} {dt}} = {\ frac {F (t)} {R ^ { 2}}}}u (R, t) = {\ frac {dR} {dt}} = {\ frac {F (t)} {R ^ {2}}}

, что дает

F (t) = R 2 d R / dt {\ displaystyle F (t) = R ^ {2} dR / dt}F (т) = R ^ {2} dR / dt

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

dm V dt = ρ V d V dt = ρ V d (4 π R 3/3) dt = 4 π ρ VR 2 d R dt {\ displaystyle {\ frac {dm_ {V}} {dt}} = \ rho _ {V} {\ frac {dV} {dt}} = \ rho _ {V} {\ frac {d (4 \ pi R ^ { 3} / 3)} {dt}} = 4 \ pi \ rho _ {V} R ^ {2} {\ frac {dR} {dt}}}{\ frac {dm_ {V}} {dt}} = \ rho _ {V} {\ frac {dV} {dt}} = \ rho _ {V} {\ frac {d (4 \ pi R ^ {3 } / 3)} {dt}} = 4 \ pi \ rho _ {V} R ^ {2} {\ frac {dR} {dt}}

с V {\ displaystyle V}V - объем пузыря. Если u L {\ displaystyle u_ {L}}u_ {L} - скорость жидкости относительно пузыря в r = R {\ displaystyle r = R}r = R , тогда масса, попадающая в пузырек, определяется выражением

dm L dt = ρ LA u L = ρ L (4 π R 2) u L {\ displaystyle {\ frac {dm_ {L}} {dt}} = \ rho _ {L} Au_ {L} = \ rho _ {L} (4 \ pi R ^ {2}) u_ {L}}{\ frac {dm_ {L}} {dt}} = \ rho _ {L} Au_ {L} = \ rho _ {L} (4 \ pi R ^ {2}) u_ {L}

с A {\ displaystyle A}A - площадь поверхности пузыря. Теперь, сохраняя массу dmv / dt = dm L / dt {\ displaystyle dm_ {v} / dt = dm_ {L} / dt}dm_{v}/dt=dm_{L}/dt, поэтому u L = (ρ V / ρ L) d R / dt {\ displaystyle u_ {L} = (\ rho _ {V} / \ rho _ {L}) dR / dt}u_ {L} = (\ rho _ {V} / \ rho _ {L}) dR / dt . Следовательно,

u (R, t) = d R dt - u L = d R dt - ρ V ρ L d R dt = (1 - ρ V ρ L) d R dt {\ displaystyle u (R, t) = {\ frac {dR} {dt}} - u_ {L} = {\ frac {dR} {dt}} - {\ frac {\ rho _ {V}} {\ rho _ {L}}} {\ frac {dR} {dt}} = \ left (1 - {\ frac {\ rho _ {V}} {\ rho _ {L}}} \ right) {\ frac {dR} {dt}}}u (R, t) = {\ frac {dR} {dt}} -u_ {L} = {\ frac {dR} {dt}} - {\ frac {\ rho _ {V}} {\ rho _ {L}}} {\ frac {dR} {dt}} = \ left (1 - {\ frac {\ rho _ {V}} {\ rho _ {L}}} \ right) {\ frac {dR} {dt}}

Следовательно,

F (t) = (1 - ρ V ρ L) R 2 d R dt {\ displaystyle F (t) = \ left (1 - {\ frac {\ rho _ {V}} {\ rho _ {L}}} \ right) R ^ {2} {\ frac {dR} {dt}}}F (t) = \ left (1 - {\ frac {\ rho _ {V}} {\ rho _ {L}}} \ right) R ^ {2} {\ frac {dR} {dt}}

Во многих случаях плотность жидкости намного больше, чем плотность пара, ρ L ≫ ρ V {\ displaystyle \ rho _ {L} \ gg \ rho _ {V}}\ rho _ {L} \ gg \ rho _ {V} , так что F (t) {\ displaystyle F (t)}F(t)может быть аппроксимировано исходной формой нулевого массопереноса F (t) = R 2 d R / dt {\ displaystyle F (t) = R ^ {2} dR / dt}F (т) = R ^ {2} dR / dt , так что

U (r, t) = F (t) r 2 = R 2 r 2 d R dt {\ displaystyle u (r, t) = {\ frac {F (t)} {r ^ {2}}} = {\ frac {R ^ {2}} {r ^ {2}}} {\ frac {dR} {dt}}}u (r, t) = {\ frac {F (t)} {r ^ {2}}} = {\ frac {R ^ {2}} {r ^ {2}}} {\ frac {dR} {dt}}

Сохранение импульса

Предполагая, что жидкость является Ньютоновская жидкость, несжимаемая Navi Уравнение Эр – Стокса в сферических координатах для движения в радиальном направлении дает

ρ L (∂ u ∂ t + u ∂ u ∂ r) = - ∂ P ∂ r + μ L [ 1 р 2 ∂ ∂ р (р 2 ∂ u ∂ р) - 2 ur 2] {\ displaystyle \ rho _ {L} \ left ({\ frac {\ partial u} {\ partial t}} + u {\ frac {\ partial u} {\ partial r}} \ right) = - {\ frac {\ partial P} {\ partial r}} + \ mu _ {L} \ left [{\ frac {1} {r ^ { 2}}} {\ frac {\ partial} {\ partial r}} \ left (r ^ {2} {\ frac {\ partial u} {\ partial r}} \ right) - {\ frac {2u} { r ^ {2}}} \ right]}\ rho _ {L} \ left ({\ frac {\ partial u} {\ partial t}} + u {\ frac {\ partial u} {\ partial r} } \ right) = - {\ frac {\ partial P} {\ partial r}} + \ mu _ {L} \ left [{\ frac {1} {r ^ {2}}} {\ frac {\ partial } {\ partial r}} \ left (r ^ {2} {\ frac {\ partial u} {\ partial r}} \ right) - {\ frac {2u} {r ^ {2}}} \ right]

Подставляем кинематическая вязкость ν L = μ L / ρ L {\ displaystyle \ nu _ {L} = \ mu _ {L} / \ rho _ {L}}\ nu _ {L} = \ mu _ {L} / \ rho _ { L} и перестановка дает

- 1 ρ L ∂ P ∂ r = ∂ u ∂ t + u ∂ u ∂ r - ν L [1 r 2 ∂ ∂ r ( р 2 ∂ u ∂ r) - 2 ur 2] {\ displaystyle - {\ frac {1} {\ rho _ {L}}} {\ frac {\ partial P} {\ partial r}} = {\ frac { \ partial u} {\ partial t}} + u {\ frac {\ partial u} {\ partial r}} - \ nu _ {L} \ left [{\ frac {1} {r ^ {2}}} {\ frac {\ partial} {\ partial r}} \ left (r ^ {2} {\ frac {\ partial u} {\ partial r}} \ right) - {\ frac {2u} {r ^ {2 }}} \ right]}- {\ frac {1} {\ rho _ {L}}} {\ frac {\ partial P} {\ partial r}} = {\ frac {\ partial u} {\ partial t} } + u {\ frac {\ partial u} {\ partial r}} - \ nu _ {L} \ left [{\ frac {1} {r ^ {2}}} {\ frac {\ partial} {\ частичный r}} \ left (r ^ {2} {\ frac {\ partial u} {\ partial r}} \ right) - {\ frac {2u} {r ^ {2}}} \ right]

при этом sub сшивание u (r, t) {\ displaystyle u (r, t)}u (r, t) из сохранения массы дает

- 1 ρ L ∂ P ∂ r = 2 R r 2 (d R dt) 2 + R 2 r 2 d 2 R dt 2 - 2 R 4 r 5 (d R dt) 2 = 1 r 2 (2 R (d R dt) 2 + R 2 d 2 R dt 2) - 2 R 4 р 5 (d R dt) 2 {\ displaystyle - {\ frac {1} {\ rho _ {L}}} {\ frac {\ partial P} {\ partial r}} = {\ frac {2R} {r ^ {2}}} \ left ({\ frac {dR} {dt}} \ right) ^ {2} + {\ frac {R ^ {2}} {r ^ {2}}} {\ frac {d ^ {2} R} {dt ^ {2}}} - {\ frac {2R ^ {4}} {r ^ {5}}} \ left ({\ frac {dR} {dt}} \ right) ^ {2} = {\ frac {1} {r ^ {2}}} \ left (2R \ left ({\ frac {dR} {dt}} \ right) ^ {2} + R ^ {2} {\ frac {d ^ {2} R} {dt ^ {2}}} \ right) - {\ frac {2R ^ {4}} {r ^ {5}}} \ left ({\ frac {dR} {dt }} \ right) ^ {2}}- {\ frac {1} {\ rho _ {L}}} {\ frac {\ partial P} {\ partial r}} = {\ frac {2R} {r ^ {2}}} \ left ({\ frac {dR} {dt}} \ right) ^ {2} + {\ frac {R ^ {2}} { r ^ {2}}} {\ frac {d ^ {2} R} {dt ^ {2}}} - {\ frac {2R ^ {4}} {r ^ {5}}} \ left ({\ frac {dR} {dt}} \ right) ^ {2} = {\ frac {1} {r ^ {2}}} \ left (2R \ left ({\ frac {dR} {dt}} \ right) ^ {2} + R ^ {2} {\ frac {d ^ {2} R} {dt ^ {2}}} \ right) - {\ frac {2R ^ {4}} {r ^ {5}} } \ left ({\ frac {dR} {dt}} \ right) ^ {2}

Обратите внимание, что вязкие члены отменяются во время замены. Разделение переменных и интегрирование от границы пузырька r = R {\ displaystyle r = R}r = R к r → ∞ {\ displaystyle r \ rightarrow \ infty}r \ rightarrow \ infty дает

- 1 ρ L ∫ P (R) P (∞) d P = ∫ R ∞ [1 r 2 (2 R (d R dt) 2 + R 2 d 2 R dt 2) - 2 R 4 r 5 (d R dt) 2] dr {\ displaystyle - {\ frac {1} {\ rho _ {L}}} \ int _ {P (R)} ^ {P (\ infty)} dP = \ int _ {R} ^ {\ infty} \ left [{\ frac {1} {r ^ {2}}} \ left (2R \ left ({\ frac {dR} {dt}} \ right) ^ {2} + R ^ {2} {\ frac { d ^ {2} R} {dt ^ {2}}} \ right) - {\ frac {2R ^ {4}} {r ^ {5}}} \ left ({\ frac {dR} {dt}} \ right) ^ {2} \ right] dr}{\ displaystyle - {\ frac { 1} {\ rho _ {L}}} \ int _ {P (R)} ^ {P (\ infty)} dP = \ int _ {R} ^ {\ infty} \ left [{\ frac {1} {r ^ {2}}} \ left (2R \ left ({\ frac {dR} {dt}} \ right) ^ {2} + R ^ {2} {\ frac {d ^ {2} R} { dt ^ {2}}} \ right) - {\ frac {2R ^ {4}} {r ^ {5}}} \ left ({\ frac {dR} {dt}} \ right) ^ {2} \ справа] dr}
P (R) - P ∞ ρ L = [- 1 r (2 R (d R dt) 2 + R 2 d 2 R dt 2) + R 4 2 р 4 (d R dt) 2] R ∞ = R d 2 R dt 2 + 3 2 (d R dt) 2 {\ displaystyle {{\ frac {P (R) -P _ {\ infty}} {\ rho _ {L}}} = \ left [- {\ frac {1} {r}} \ left (2R \ left ({\ frac {dR} {dt}} \ right) ^ {2} + R ^ {2 } {\ frac {d ^ {2} R} {dt ^ {2}}} \ right) + {\ frac {R ^ {4}} {2r ^ {4}}} \ left ({\ frac {dR } {dt}} \ right) ^ {2} \ right] _ {R} ^ {\ infty} = R {\ frac {d ^ {2} R} {dt ^ {2}}} + {\ frac { 3} {2}} \ left ({\ frac {dR} {dt}} \ right) ^ {2}}}{{\ frac {P (R) -P _ {\ infty}} {\ rho _ {L}}} = \ left [- {\ frac {1} {r}} \ left (2R \ left ({\ frac {dR} {dt}} \ right) ^ {2} + R ^ {2} {\ frac {d ^ {2} R} {dt ^ {2}}} \ right) + { \ frac {R ^ {4}} {2r ^ {4}}} \ left ({\ frac {dR} {dt}} \ right) ^ {2} \ right] _ {R} ^ {\ infty} = R {\ frac {d ^ {2} R} {dt ^ {2}}} + {\ frac {3} {2}} \ left ({\ frac {dR} {dt}} \ right) ^ {2}}

Граничные условия

Пусть σ rr {\ displaystyle \ sigma _ {rr}}\ sigma _ {{rr}} быть нормальным напряжением в жидкости, которое направлено радиально наружу от центра пузыря. В сферических координатах для жидкости с постоянной плотностью и постоянной вязкостью

σ rr = - P + 2 μ L ∂ u ∂ r {\ displaystyle \ sigma _ {rr} = - P + 2 \ mu _ {L} {\ frac {\ partial u} {\ partial r}}}\ sigma _ {{rr}} = - P + 2 \ mu _ {L} {\ frac {\ partial u} {\ partial r}}

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

σ rr (R) + PB - 2 γ R = - P (R) + 2 μ L ∂ u ∂ r | r = R + PB - 2 γ R = - P (R) + 2 μ L ∂ ∂ r (R 2 r 2 d R dt) r = R + PB - 2 γ R = - P (R) - 4 μ LR d R dt + PB - 2 γ R {\ displaystyle {\ begin {align} \ sigma _ {rr} (R) + P_ {B} - {\ frac {2 \ gamma} {R}} = - P ( R) + \ left.2 \ mu _ {L} {\ frac {\ partial u} {\ partial r}} \ right | _ {r = R} + P_ {B} - {\ frac {2 \ gamma} {R}} \\ = - P (R) +2 \ mu _ {L} {\ frac {\ partial} {\ partial r}} \ left ({\ frac {R ^ {2}} {r ^ {2}}} {\ frac {dR} {dt}} \ right) _ {r = R} + P_ {B} - {\ frac {2 \ gamma} {R}} \\ = - P (R) - {\ frac {4 \ mu _ {L}} {R}} {\ frac {dR} {dt}} + P_ {B} - {\ frac {2 \ gamma} {R}} \\\ конец {выравнивается}}}{\ displaystyle {\ begin {align} \ sigma _ {rr} (R) + P_ {B} - {\ frac {2 \ gamma} {R}} = - P (R) + \ left.2 \ mu _ {L} {\ frac {\ partial u} {\ partial r}} \ right | _ {r = R} + P_ {B} - {\ frac {2 \ gamma} {R}} \ \ = - P (R) +2 \ mu _ {L} {\ frac {\ partial} {\ partial r}} \ left ({\ frac {R ^ {2}} {r ^ {2}}} {\ frac {dR} {dt}} \ right) _ {r = R} + P_ {B} - {\ frac {2 \ gamma} {R}} \\ = - P (R) - {\ frac {4 \ mu _ {L}} {R}} {\ frac {dR} {dt}} + P_ {B} - {\ frac {2 \ gamma} {R}} \\\ конец {выровнено}}}

, где γ {\ displaystyle \ gamma}\ gamma - это поверхностное натяжение. Если нет массопереноса через границу, то эта сила на единицу площади должна быть равна нулю, поэтому

P (R) = PB - 4 μ LR d R dt - 2 γ R {\ displaystyle P (R) = P_ {B} - {\ frac {4 \ mu _ {L}} {R}} {\ frac {dR} {dt}} - {\ frac {2 \ gamma} {R}}}{\ displaystyle P (R) = P_ {B} - {\ frac {4 \ mu _ {L}} {R}} {\ frac {dR} {d t}} - {\ frac {2 \ gamma} {R}}}

и поэтому результат сохранения импульса становится

P (R) - P ∞ ρ L = PB - P ∞ ρ L - 4 μ L ρ LR d R dt - 2 γ ρ LR = R d 2 R dt 2 + 3 2 (d R dt) 2 {\ displaystyle {\ frac {P (R) -P _ {\ infty}} {\ rho _ {L}}} = {\ frac {P_ {B} -P _ {\ infty}} {\ rho _ {L}}} - {\ frac {4 \ mu _ {L}} {\ rho _ {L} R}} {\ frac {dR} {dt}} - {\ frac {2 \ gamma} {\ rho _ {L} R}} = R {\ frac {d ^ {2} R} {dt ^ {2}}} + {\ frac {3} {2}} \ left ({\ frac {dR} { dt}} \ right) ^ {2}}{\ displaystyle {\ frac {P (R) -P _ {\ infty}} {\ rho _ {L}}} = {\ frac {P_ {B} -P _ {\ infty}} {\ rho _ {L}}} - {\ frac {4 \ mu _ {L}} {\ rho _ {L} R}} {\ frac {dR} {dt}} - {\ frac {2 \ gamma} {\ rho _ {L} R}} = R {\ frac {d ^ {2} R} {dt ^ {2}}} + {\ frac {3 } {2}} \ left ({\ frac {dR} {dt}} \ right) ^ {2}}

, переставляя и позволяя ν L = μ L / ρ L {\ displaystyle \ nu _ {L} = \ mu _ {L} / \ rho _ { L}}\ nu _ {L} = \ mu _ {L} / \ rho _ { L} дает уравнение Рэлея – Плессета

PB (t) - P ∞ (t) ρ L = R d 2 R dt 2 + 3 2 (d R dt) 2 + 4 ν LR d R dt + 2 γ ρ LR {\ displaystyle {\ frac {P_ {B} (t) -P _ {\ infty} (t)} {\ rho _ {L}}} = R {\ frac {d ^ { 2} R} {dt ^ {2}}} + {\ frac {3} {2}} \ left ({ \ frac {dR} {dt}} \ right) ^ {2} + {\ frac {4 \ nu _ {L}} {R}} {\ frac {dR} {dt}} + {\ frac {2 \ gamma} {\ rho _ {L} R}}}{\ displaystyle {\ frac {P_ {B} (t) -P _ {\ infty} (t) } {\ rho _ {L}}} = R {\ frac {d ^ {2} R} {dt ^ {2}}} + {\ frac {3} {2}} \ left ({\ frac {dR } {dt}} \ right) ^ {2} + {\ frac {4 \ nu _ {L}} {R}} {\ frac {dR} {dt}} + {\ frac {2 \ gamma} {\ rho _ {L} R}}}

Используя точечную нотацию для представления производных по времени, уравнение Рэлея – Плессета можно более кратко записать как

PB (t) - P ∞ (t) ρ L знак равно RR ¨ + 3 2 (R ˙) 2 + 4 ν LR ˙ R + 2 γ ρ LR {\ displaystyle {\ frac {P_ {B} (t) -P _ {\ infty } (t)} {\ rho _ {L}}} = R {\ ddot {R}} + {\ frac {3} {2}} ({\ dot {R}}) ^ {2} + {\ frac {4 \ nu _ {L} {\ dot {R}}} {R}} + {\ frac {2 \ gamma} {\ rho _ {L} R}}}{\ displaystyle {\ frac {P_ {B} (t) -P _ {\ infty} (t)} {\ rho _ {L}}} = R {\ ddot {R}} + {\ frac { 3} {2}} ({\ dot {R}}) ^ {2} + {\ frac {4 \ nu _ {L} {\ dot {R}}} {R}} + {\ frac {2 \ gamma} {\ rho _ {L} R }}}
Решения

Недавно аналитические решения в замкнутой форме были найдены для уравнения Рэлея – Плессета как для пустого, так и для газового пузыря и были обобщены на N-мерный случай. Также был изучен случай, когда поверхностное натяжение присутствует из-за эффектов капиллярности.

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

В статическом случае уравнение Рэлея – Плессета упрощается, давая уравнение Юнга-Лапласа :

PB - P ∞ = 2 γ R {\ displaystyle P_ {B} -P _ {\ infty} = { \ frac {2 \ gamma} {R}}}{\ displaystyle P_ {B} -P _ {\ infty} = {\ frac { 2 \ gamma} {R}}}

Когда рассматриваются только бесконечно малые периодические изменения радиуса пузыря и давления, уравнение RP также дает выражение собственной частоты колебаний пузыря.

Ссылки
Последняя правка сделана 2021-06-03 09:25:33
Содержание доступно по лицензии CC BY-SA 3.0 (если не указано иное).
Обратная связь: support@alphapedia.ru
Соглашение
О проекте