Физика · Гидродинамика

Почему мяч закручивается?

Эффект Магнуса, уравнение Бернулли и парадокс Даламбера

Штрафной обходит стенку и заходит под штангу, ни с чем не столкнувшись. Выводим уравнение Бернулли, разбираем обтекание цилиндра и парадокс Даламбера, получаем силу Магнуса и считаем настоящий удар.

  • Арцгольд Максим Дмитриевичфизик, главный редактор, руководитель "Академии"
  • Ребрикова Мария Максимовнавыпускающий редактор

Листай вниз

01постановка задачи

Поворот мяча в полёте

Штрафной с двадцати пяти метров. Мяч идёт правее ворот, мимо стенки, а на последних метрах сворачивает внутрь и заходит под штангу. Ни с чем он при этом не сталкивается.

Поперечной силе взяться вроде бы неоткуда, ведь тяжесть тянет вниз, а сопротивление воздуха назад. Но мяч ещё и вращается, и это меняет картину обтекания. С одной стороны воздух проходит мимо мяча быстрее, с другой медленнее, а там, где течение быстрее, давление ниже. Разность давлений и уводит мяч в сторону.

Дальше выведем уравнение Бернулли из закона сохранения энергии, разберём обтекание цилиндра и парадокс Даламбера, введём циркуляцию Γ\Gamma и получим силу F=ρvΓF=\rho v\Gamma. В конце посчитаем настоящий штрафной.

02гидродинамика

Уравнение Бернулли

Возьмём трубку тока — воображаемую трубку, стенки которой сотканы из линий тока, так что жидкость сквозь них не течёт. Пусть за время dtdt через вход прошёл объём dVdV. Жидкость несжимаема, поэтому через выход прошёл ровно такой же объём.

Работа сил давления на входе и выходе идёт на изменение кинетической и потенциальной энергии:

(p1−p2) dV=ρ dV v222−ρ dV v122+ρ dV g(h2−h1).(p_1-p_2)\,dV=\frac{\rho\,dV\,v_2^{2}}{2}-\frac{\rho\,dV\,v_1^{2}}{2}+\rho\,dV\,g(h_2-h_1).

Сократим на объёме и соберём слагаемые по сечениям. Получится уравнение Бернулли:

p+ρv22+ρgh=const вдоль линии тока.p+\frac{\rho v^{2}}{2}+\rho g h=\text{const вдоль линии тока.}

это закон сохранения энергии, записанный на единицу объёма жидкости

Когда им можно пользоваться

Вывод опирался на три допущения, и каждое ограничивает применимость. Течение стационарно, то есть картина не меняется со временем. Жидкость несжимаема, что для воды верно всегда, а для воздуха пока скорость много меньше скорости звука. Вязкости нет, и часть энергии не уходит в тепло.

Постоянная при этом своя для каждой линии тока. Поэтому сравнивать давления в двух точках по Бернулли можно, только если через них проходит одна и та же линия.

Сузь горло. Частица в нём разгоняется, а столбик манометра над ним опускается, ведь где быстрее, там давление ниже.

Из неразрывности следует, что в узком месте скорость выше, ведь через любое сечение за секунду проходит один и тот же объём. Значит, там ниже и давление. На этом работают распылитель, водоструйный насос и расходомер Вентури, который меряет расход по перепаду давления.

v1S1=v2S2⟹Δp=ρ2(v22−v12).v_1 S_1=v_2 S_2\quad\Longrightarrow\quad \Delta p=\frac{\rho}{2}\left(v_2^{2}-v_1^{2}\right).
03обтекание

Линии тока вокруг тела

Поставим в поток цилиндр. Если жидкость идеальная, задачу можно решить точно, ведь течение описывается функцией тока, линии уровня которой и есть линии тока. Для цилиндра радиуса aa в потоке UU функция тока равна

ψ=Uy(1−a2x2+y2).\psi=U y\left(1-\frac{a^{2}}{x^{2}+y^{2}}\right).

Скорость получается дифференцированием: u=∂ψ/∂yu=\partial\psi/\partial y, v=−∂ψ/∂xv=-\partial\psi/\partial x. На самой поверхности радиальная составляющая обращается в ноль — жидкость сквозь тело не течёт, — а касательная равна

q(θ)=−2Usin⁡θ.q(\theta)=-2U\sin\theta.

В точках θ=0\theta=0 и θ=180∘\theta=180^\circ скорость равна нулю, это и есть точки торможения, спереди и сзади. На боках, при θ=±90∘\theta=\pm90^\circ, скорость вдвое больше скорости набегающего потока. Там, по Бернулли, глубокое разрежение.

Добавь циркуляцию. Сверху линии тока сгущаются, там воздух быстрее и давление ниже, и на цилиндр действует направленная вверх сила Магнуса.

Пока циркуляции нет, картина симметрична и сверху вниз, и спереди назад. Симметрия сверху вниз означает, что поперечной силы нет. А симметрия спереди назад приводит к выводу, который расходится с опытом, и им мы займёмся в следующем разделе.

04парадокс

Парадокс Даламбера

Посчитаем силу сопротивления цилиндра в идеальной жидкости. Давление на поверхности даёт Бернулли:

p(θ)=p∞+ρ2(U2−4U2sin⁡2θ).p(\theta)=p_\infty+\frac{\rho}{2}\left(U^{2}-4U^{2}\sin^{2}\theta\right).

Оно зависит только от sin⁡2θ\sin^{2}\theta, то есть одинаково спереди и сзади. Проинтегрируем проекцию силы давления на направление потока по всей окружности и получим ноль.

Значит, ни цилиндр, ни шар, ни любое другое тело в идеальной жидкости не испытывает силы сопротивления. Это и есть парадокс Даламбера, сформулированный в 1752 году.

Вывод расходится с опытом, ведь грести веслом тяжело, а велосипедист на скорости тратит на воздух больше половины мощности. Значит, какое-то из допущений не работает, и это отсутствие вязкости.

Пограничный слой и отрыв

Вязкость воздуха мала, и в основной части потока ею действительно можно пренебречь. Но у самой поверхности скорость обязана обращаться в ноль, потому что молекулы прилипают к телу. Значит, есть тонкий слой, в котором скорость быстро падает от UU до нуля. Его называют пограничным, и толщина его невелика:

δ∼νLU=LRe.\delta\sim\sqrt{\frac{\nu L}{U}}=\frac{L}{\sqrt{\mathrm{Re}}}.

За серединой тела давление начинает расти, ведь жидкость, разогнавшаяся на боку, должна затормозить у кормы. Внешнему потоку энергии на это хватает, а обеднённому пограничному слою — нет. Слой останавливается, разворачивается и отрывается от поверхности.

За телом образуется след с пониженным давлением, симметрия спереди назад ломается, и появляется настоящее сопротивление. Поэтому у обтекаемых тел длинный сужающийся хвост, который отодвигает отрыв назад.

вязкость мала, но именно она решает, где поток оторвётся

Тот же пограничный слой понадобится нам через раздел. Вращающееся тело увлекает его за собой, и так у мяча появляется циркуляция.

05циркуляция

Циркуляция и подъёмная сила

Циркуляция — это интеграл скорости по замкнутому контуру вокруг тела:

Γ=∮v⃗⋅dl⃗.\Gamma=\oint \vec v\cdot d\vec l.

К решению для цилиндра её можно добавить, не нарушив ни одного уравнения. Течение остаётся безвихревым всюду вне тела, а сквозь поверхность по-прежнему ничего не течёт. Функция тока получает логарифмическую добавку:

ψ=Uy(1−a2r2)−Γ2πln⁡ra.\psi=U y\left(1-\frac{a^{2}}{r^{2}}\right)-\frac{\Gamma}{2\pi}\ln\frac{r}{a}.

А скорость на поверхности — постоянное слагаемое:

q(θ)=−2Usin⁡θ+Γ2πa.q(\theta)=-2U\sin\theta+\frac{\Gamma}{2\pi a}.

Этим и исчерпывается несимметрия. С одной стороны циркуляция складывается со скоростью обтекания, с другой вычитается. Точки торможения съезжают и стоят теперь там, где sin⁡θ=Γ/(4πaU)\sin\theta=\Gamma/(4\pi a U). Когда циркуляция превышает 4πaU4\pi a U, точки сливаются и вовсе сходят с поверхности.

Считаем силу

Подставим скорость в Бернулли и проинтегрируем давление по окружности. Члены с sin⁡3θ\sin^{3}\theta и просто sin⁡θ\sin\theta за полный оборот дают ноль, выживает только перекрёстное слагаемое:

Fy=−∮psin⁡θ a dθ=−ρUΓ.F_y=-\oint p\sin\theta\,a\,d\theta=-\rho U\Gamma.

Подъёмная сила равна произведению плотности, скорости потока и циркуляции.

Это теорема Жуковского. Сила перпендикулярна потоку, а размер тела и вязкость в неё не входят, ведь вся геометрия задачи спрятана в циркуляции.

Знак читается так. Если поверхность тела движется по потоку с той стороны, где мы ждём подъёмную силу, там воздух разгоняется, давление падает, и туда тело и потянет.

Откуда циркуляция у мяча

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

Γ∼2πa⋅ωa=2πa2ω.\Gamma\sim 2\pi a\cdot\omega a=2\pi a^{2}\omega.

Подставим в теорему Жуковского и получим F∼2πρa2ωUF\sim 2\pi\rho a^{2}\omega U, так что сила растёт и с вращением, и со скоростью. Для шара строгой формулы такого вида нет, потому что след турбулентный, и на практике пользуются коэффициентом подъёмной силы, измеренным в трубе:

F=CL ρU22 πa2,CL≈0,6 ωaU.F=C_L\,\frac{\rho U^{2}}{2}\,\pi a^{2},\qquad C_L\approx0{,}6\,\frac{\omega a}{U}.

Коэффициент насыщается около 0,30{,}3 при сильном вращении. Дальше формула перестаёт работать, потому что вступает в дело перестройка следа и зависимость становится нелинейной.

06практика

Футбол, теннис, бейсбол

Оценим, на сколько уводит футбольный мяч. Радиус a=0,11a=0{,}11 м, масса 0,430{,}43 кг, скорость 2525 м/с, вращение 88 оборотов в секунду, то есть ω≈50\omega\approx50 рад/с.

S=ωaU=50⋅0,1125≈0,22,CL≈0,13.S=\frac{\omega a}{U}=\frac{50\cdot0{,}11}{25}\approx0{,}22,\qquad C_L\approx0{,}13.
F=CLρU22πa2=0,13⋅1,2⋅6252⋅0,038≈1,9 Н.F=C_L\frac{\rho U^{2}}{2}\pi a^{2}=0{,}13\cdot\frac{1{,}2\cdot625}{2}\cdot0{,}038\approx1{,}9\ \text{Н}.

Ускорение получается a⊥=F/m≈4,4a_\perp=F/m\approx4{,}4 м/с². За секунду полёта это даёт a⊥t2/2≈2a_\perp t^{2}/2\approx2 метра вбок — ровно столько, чтобы обойти стенку и вернуться в створ.

———

Покрути сцену: пунктир — тот же удар без вращения, розовая кривая — с вращением, и сбоку видно, что увод лежит в горизонтальной плоскости. Поставь вращение в ноль, и линии совпадут.

Боковое ускорение почти постоянно, а смещение набирается как a⊥t2/2a_\perp t^{2}/2. Поэтому увод растёт как квадрат времени полёта, и на длинном штрафном подкрутка заметна, а на коротком пасе почти нет.

Пробей. Закрученный мяч обходит стенку снаружи и сворачивает в створ, а пунктиром показан тот же удар без вращения, он уходит мимо.

Другие виды спорта

В теннисе верхнее вращение направляет силу Магнуса вниз, мяч падает круче параболы, и бить можно сильнее, не рискуя вылететь за линию. Нижнее вращение действует наоборот, оно удерживает мяч в воздухе, и резаный удар летит дольше и ниже отскакивает.

В бейсболе подача с верхним вращением падает резче обычной, а с нижним кажется поднимающейся, хотя падает просто медленнее. Есть и подача почти без вращения, у неё неустойчивый след, и мяч виляет непредсказуемо.

В гольфе нижнее вращение задают специально, потому что сила Магнуса направлена вверх и заметно увеличивает дальность. Ямочки на поверхности работают на то же самое, они делают пограничный слой турбулентным, отодвигают отрыв и уменьшают сопротивление.

Подъёмная сила крыла

Крыло не вращается, но циркуляция у него есть, и создаёт её форма профиля вместе с условием схода потока с острой задней кромки. Подъёмная сила считается той же теоремой Жуковского, и для оценок обычно записывается через коэффициент:

L=CL ρU22 S.L=C_L\,\frac{\rho U^{2}}{2}\,S.

Отсюда следуют два практических вывода. Подъёмная сила растёт как квадрат скорости, поэтому при удвоении скорости она вырастает вчетверо. И есть минимальная скорость, ниже которой крыло не удержит вес, её и называют скоростью сваливания.

07задачи

Проверь себя

Посчитай, впиши ответ и нажми «Проверить». Разбор открывается кнопкой под задачей. Плотность воды 10001000 кг/м³, воздуха 1,21{,}2 кг/м³, g=9,8g=9{,}8 м/с². Всё нужное для решения собрано в трёх формулах.

v1S1=v2S2,p+ρv22+ρgh=const,F=ρUΓ.v_1S_1=v_2S_2,\qquad p+\frac{\rho v^{2}}{2}+\rho g h=\text{const},\qquad F=\rho U\Gamma.
1
По трубе диаметром 44 см вода течёт со скоростью 22 м/с. Какой станет скорость в сужении диаметром 22 см?
v₂ =м/с
Разбор

Решение

Через оба сечения за секунду проходит один и тот же объём:

v2=v1S1S2=v1(d1d2)2=2⋅4=8 м/с.v_2=v_1\frac{S_1}{S_2}=v_1\left(\frac{d_1}{d_2}\right)^{2}=2\cdot 4=8\ \text{м/с}.

Площадь зависит от квадрата диаметра, поэтому вдвое уже — вчетверо быстрее.

Ответ: 8 м/с
2
Для той же трубы найди перепад давления между широким и узким сечением. Ответ в килопаскалях.
Δp =кПа
Разбор

Решение

Δp=ρ2(v22−v12)=10002(64−4)=30 000 Па.\Delta p=\frac{\rho}{2}\left(v_2^{2}-v_1^{2}\right)=\frac{1000}{2}(64-4)=30\,000\ \text{Па}.

Тридцать килопаскалей — почти треть атмосферы. Если сужение сделать ещё сильнее, давление опустится ниже давления насыщенного пара, и вода вскипит прямо в трубе. Это и есть кавитация.

Ответ: 30 кПа
3
Цилиндр обтекается воздухом со скоростью 1010 м/с, циркуляция вокруг него 33 м²/с. Найди подъёмную силу на метр длины.
F =Н/м
Разбор

Решение

F=ρUΓ=1,2⋅10⋅3=36 Н/м.F=\rho U\Gamma=1{,}2\cdot 10\cdot 3=36\ \text{Н/м}.

Ни радиус цилиндра, ни вязкость в формулу не входят. Радиус прячется внутри циркуляции, именно он определяет, какую Γ\Gamma удастся раскрутить.

Ответ: 36 Н/м
4
Теннисный мяч бьют с верхним вращением, то есть верх мяча движется вперёд по полёту. Что будет с траекторией?
Разбор

Решение

Перейдём в систему мяча, где набегающий воздух идёт назад. Верх мяча при верхнем вращении движется вперёд, навстречу этому воздуху, и тормозит его. Низ движется назад, вместе с потоком, и разгоняет его.

Значит, снизу скорость обтекания больше, а давление меньше. Сила направлена вниз, и мяч падает круче параболы. Поэтому верхним вращением и бьют, ведь так можно разогнать мяч сильнее, не вылетая за линию.

Ответ: упадёт круче, чем по параболе
5
В боковой стенке сосуда сделали отверстие на глубине 22 м под уровнем воды. С какой скоростью бьёт струя?
v =м/с
Разбор

Решение

Записываем Бернулли вдоль линии тока от свободной поверхности до отверстия. И там, и там давление атмосферное, скорость на поверхности пренебрежимо мала:

ρgh=ρv22⟹v=2gh=2⋅9,8⋅2≈6,3 м/с.\rho g h=\frac{\rho v^{2}}{2}\quad\Longrightarrow\quad v=\sqrt{2gh}=\sqrt{2\cdot9{,}8\cdot2}\approx6{,}3\ \text{м/с}.

Это формула Торричелли, и она совпадает со скоростью свободного падения с той же высоты. Плотность сократилась, так что ртуть била бы так же.

Ответ: ≈ 6,3 м/с
6
Во сколько раз вырастет подъёмная сила крыла при удвоении скорости, если угол атаки тот же?
Разбор

Решение

L=CLρU22S∝U2.L=C_L\frac{\rho U^{2}}{2}S\propto U^{2}.

Вчетверо. Отсюда и понятие скорости сваливания, ниже неё подъёмной силы не хватает, чтобы удержать вес самолёта, как ни задирай нос.

Ответ: вчетверо

Поперечную силу создаёт несимметрия обтекания, а её задаёт вращение тела.

В Компендиуме разбираем физические задачи, в которых привычная интуиция требует аккуратной проверки.

залетай в Академию→