Підвищення бульбашки в в’язкому середовищі: Форма, Швидкість та Нестабільність Шляху

Buoyancy versus drag

A gas bubble released in a viscous liquid rises because it is less dense than the liquid around it - the same buoyancy that floats a boat - but it does not accelerate forever. As it speeds up, viscous drag from the surrounding fluid grows with velocity until it exactly balances the buoyant force, and the bubble settles into a constant terminal velocity. Everything interesting about how a bubble rises - whether it stays a smooth sphere, flattens into a cap, or wobbles chaotically - comes down to how a handful of dimensionless numbers weigh buoyancy, viscosity and surface tension against each other.

Стійкий опір та межа для малих бульбашок

Для невеликої бульбашки, яка рухається достатньо повільно, щоб інерція була незначною порівняно з в'язкістю (низький число Рейну), сила опору відповідає закону Стокса: F_опір = 6*pi*μ*r*v, де μ - динамічна в’язкість рідини, а r - радіус бульбашки. Встановлюючи стійкий опір рівним чистому плаваючому силі (вага витісненої рідини мінус власна вага бульбашки, яка незначна), отримуємо класичний результат термічного прискорення:

F_плавучість = (4/3)*pi*r^3 * Δρ * g F_опір = 6*pi*μ*r*v_термічного (Стійкий опір, дійсний для Re < ~1) v_термічного = (2/9) * Δρ * g * r^2 / μ // розв’язуючи F_плавучість = F_опір Це однакова формула, з інвертованим знаком, яка описує щільну сферу, що повільно опускається через в'язку рідину, і добре працює для дуже малих бульбашок - кілька десятків міліметрів або менше у воді, де поверхневий натяг достатньо сильний, щоб підтримувати бульбашку майже ідеально сферичною протягом її підйому, тому опір насправді є простим стійким опором твердого об’єкта. Більші бульбашки швидко відхиляються від цього зображення, як через тертя, яке ігнорується закон Стокса, так і через те, що форма бульбашки більше не сферична.

F_buoyancy = (4/3)*pi*r^3 * delta_rho * g
F_drag     = 6*pi*mu*r*v_terminal            (Stokes drag, valid for Re < ~1)
v_terminal = (2/9) * delta_rho * g * r^2 / mu   // solving F_buoyancy = F_drag

Два числа, що визначають форму

Число Еотвоса, Bo = дельта_ро * g * d^2 / сигма, порівнює плавучий виштовхнутий тиск, який намагається деформувати бульбашку, з поверхневою напруженістю сигма, яка намагається утримати її сферичною. Маленьке Bo (домінування поверхневої напруженості) підтримує сферичну форму бульбашки; велике Bo (домінування плавучості, зазвичай більші бульбашки) дозволяє їй вирівнятися в обтічну еліпсоїдну форму або, на межі, у форму «сферичного коронного пара» з плоским або навіть вигнутим донизу основами.

Число Мортона, Mo = g*му^4*дельта_ро / (ро^2*сигма^3), є незвичайною серед безрозмірних груп, оскільки залежить лише від властивостей рідини, а не від розміру або швидкості бульбашки - це фіксоване число для даної суміші рідкого та газоподібного середовища (вода-повітря, гліцерин-повітря тощо), яке вимірює відносну в'язкість і щільність рідин до їх поверхневої напруженості. Разом числа Еотвоса та Мортона (або еквівалентно триплетом Рейнольдса-Еотвоса-Мортона) визначають форму бульбашки та режим її підйому майже повністю, і дослідники малюють їх на стандартній довідковій діаграмі - діаграмі Грейс, яка розділяє цей параметричний простір на сферичні, еліпсоїдні, колискові та коронні області, дозволяючи інженеру визначити очікувану форму бульбашки та її кінцеву швидкість лише за властивостями рідини та розміром бульбашки, без запуску симуляції.

Path instability: why bubbles wobble and zigzag

Beyond a critical bubble size (and correspondingly Reynolds number, typically several hundred), a rising bubble’s path stops being a straight vertical line. Small asymmetric perturbations in the wake behind the bubble – vortex shedding, much like the vortex street behind a cylinder in cross-flow – couple back to the bubble’s motion and shape, producing a rectilinear-to-zigzag or rectilinear-to-spiral transition.

This path instability is a genuinely three-dimensional, unsteady phenomenon: the bubble’s wake sheds vorticity asymmetrically, pushing the bubble sideways, which changes the wake, which pushes it again, in a self-sustaining oscillation that also correlates with shape oscillations (wobbling) for bubbles that are not perfectly rigid.

Why bubble rise matters beyond the lab

Bubble dynamics governs gas-liquid mass transfer in bubble-column chemical reactors, aeration in wastewater treatment, boiling heat transfer in power plants, CO2 outgassing from carbonated drinks, and the degassing of magma in volcanic conduits, where large, rapidly rising gas bubbles are implicated in explosive eruption dynamics. In nearly all of these, the practical quantity engineers care about - the interfacial area available for mass or heat transfer, or the residence time a bubble spends in a reactor - depends directly on the terminal velocity and shape regime the Eotvos-Morton-Reynolds framework predicts.

Frequently asked questions

Спробуйте наживо

Усе, що вище, працює прямо у вашому браузері — відкрийте Rising Bubble і змінюйте параметри під час роботи. Нічого не встановлюється, нічого не завантажується на сервер, уся модель живе в одній вкладці.