Закон Фарадея или как магнит застревает в медной трубе

Магнит в медной трубе
Изображение взято с сайта «Популярная механика»

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

И если вы дочитали до этого места и не испугались, добро пожаловать под кат — дальше будет интереснее!
Сам я давно подумывал над тем, чтобы хорошенько разобраться в этом вопросе. И вот недавно зашёл разговор с коллегой по работе. Его ребёнку задали сделать научную демонстрацию в школе, на что папа раздобыл кусок медной трубы и неодим-железо-борный магнит. Ребёнок разобрался, произвёл демонстрацию опыта перед классом, дал пояснения, но ни класс ни учитель особо не впечатлились. На конкурсе научных опытов победил вулкан (!) из соды и лимонной кислоты =) Мы с коллегой прикинули на словах и поняли, что дело ясное, что дело тёмное. Да и в литературе не особо много написано по данной тематике. Поговорили и разошлись. Ну, а я всё-таки решил попробовать продраться сквозь дебри. В этой статье пишу, что у меня получилось.

Описание эксперимента



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

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

Итак, в опыте мы наблюдаем, как постоянный магнит движется внутри полой медной трубы с постоянной скоростью. Зафиксируем произвольную точку в теле медной трубки и мысленно проведем поперечное сечение. Через данное сечение медной трубы проходит магнитный поток, содаваемый постоянным магнитом. Из-за того, что магнит движется вдоль трубы, в сечении проводника возникает переменный магнитный поток, то ли нарастающий, то ли убывающий в зависимости от того, приближается или отдаляется магнит от точки, где мы мысленно провели сечение. Переменный магнитный поток, согласно уравнениям Максвелла, порождает вихревое электрическое поле, вообще говоря, во всём пространстве. Однако, только там, где есть проводник, это электрическое поле приводит в движение свободные заряды, находящиеся в проводнике — возникает круговой электрический ток, который создает уже своё собственное магнитное поле и взаимодействует с магнитным полем движущегося постоянного магнита. Проще говоря, круговой электрический ток создает магнитное поле того же знака, что и постоянный магнит, и на магнит действует некая диссипативная сила, а если конкретно — сила трения. Читатель может справедливо задать вопрос: «Трение чего обо что?» Трение возникает между магнитным полем диполя и проводником. Да, это трение не механическое. Вернее сказать, тела не соприкасаются. Ну и пусть! Трение всё равно есть!

В целом, на словах всё выглядит более или менее складно, а можно ли это описать на языке математики? Приступим…

Математическое описание


Постоянный магнит в медной трубе

Перво-наперво, нам понадобится математическая модель постоянного магнита. На мой взгляд, будет удобно представить постоянный магнит как магнитный диполь.

$\vec{B}=\frac{\mu_0}{4\pi}\left(\frac{3(\vec{p}_m\cdot\vec{r})\vec{r}}{r^5}-\frac{\vec{p}_m}{r^3}\right)$


Здесь приняты обозначения $\vec{r}=(r,z)$ — радиус-вектор из центра диполя в точку наблюдения, $\vec{p}_m$ — вектор дипольного момента.

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

$B_z(r,z) = \frac{\mu_0\,p_m}{4\pi}\frac{2z^2 - r^2}{\left(r^2 + z^2\right)^\frac{5}{2}}$


Теперь запишем выражение для магнитного потока через площадь, охватываемую окружностью радиуса $r$ на расстоянии $z$ от диполя.

$\Phi(r,z) = \int_0^{2\pi}\int_0^r B_z(r',z)\,r'\,dr'\,d\varphi = 2\pi\int_0^r \frac{\mu_0\,p_m}{4\pi}\frac{2z^2 - r'^2}{\left(r'^2 + z^2\right)^\frac{5}{2}}\,r'\,dr'$


Вы не поверите, но этот интеграл берётся. Не буду утомлять. В ответе получается очень красиво

$\Phi(r,z) = \frac{\mu_0\,p_m}{2}\frac{r^2}{\left(r^2 + z^2\right)^\frac{3}{2}}$


Из-за того, что диполь движется вдоль оси $z$ со скоростью $v$, нужно также сделать стандартную подстановку $\Phi(r, z)\rightarrow \Phi(r, z - v t)$
Похоже, пора призвать на помощь одно из великих уравнений Максвелла, а именно, то самое уравнение, которое описывает закон Фарадея:

Изменение потока магнитной индукции, проходящего через незамкнутую поверхность $S$, взятое с обратным знаком, пропорционально циркуляции электрического поля на замкнутом контуре $L$, который является границей поверхности $S$

$\oint_L \vec{E}\,dl = -\frac{\partial}{\partial t}\int_S \vec{B}\,d\vec{s}$


Или, что то же самое,

$2\pi r E_\varphi = -\frac{\partial}{\partial t}\Phi(r,z - v t)$


Здесь мы воспользовались аксиальной симметрией задачи по отношению к оси $z$, а также учли, что индуцированное электрическое поле имеет только азимутальную компоненту $\vec{E} = E_\varphi \vec{e}_\varphi$.
Отсюда можно найти азимутальную компоненту электрического поля, индуцированного магнитом.

$E_\varphi(r,z) = -\frac{1}{2\pi r}\frac{\partial}{\partial t}\Phi(r, z - v t)=-\frac{3\mu_0\,p_m}{4\pi}\frac{rv(z-v t)}{\left(r^2 + (z-v t)^2\right)^\frac{5}{2}}$


Теперь, когда у нас есть выражение для электрического поля, можно вспомнить и о трубе. Как показано на рисунке выше, внутренний радиус трубы равен $a$, а внешний — $b$. Материал трубы — медь. В данный момент нам будет нужна только электрическая проводимость меди. Обозначим проводимость за $\sigma$.
Электрическое поле внутри проводника вызывает электрический ток. Поэтому можем записать закон Ома в дифференциальной форме

$\vec{j} = \sigma \vec{E}$


Электрический ток, в свою очередь вызывает омические потери внутри проводника. Иными словами, энергия рассеивается внутри проводника и переходит в форму тепла, строго говоря, в нашем случае во всём объёме проводника.
Объёмная плотность мощности омических потерь по определению равна

$w = \vec{j}\cdot\vec{E} = \sigma E^2$


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

$\frac{dE_p}{dt} = P$

$-mg\dot{z} = \int_V w\,dV$

$mgv = \int_{-\infty}^{\infty}\int_0^{2\pi}\int_a^b \sigma E^2\,r\,dr\,d\varphi\,dz$


Здесь необходимо заметить, что потенциальная энергия в координатах, изображенных на рисунке выше будет равна $E_p = - mgz$, а чтобы найти полную мощность омических потерь, следует проинтегрировать $w$ по всему объёму проводника. Длину трубы считаем бесконечной. Это не так далеко от истины, если учесть, что в опыте из видеоролика диаметр магнитика много меньше длины трубы.

Последний тройной интеграл выглядит очень сложным. И так оно и есть! Но, во-первых, интегрирование по азимутальному углу $\varphi$ можно заменить просто домножением на $2\pi$ в силу аксиальной симметрии задачи. Во-вторых, порядок интегрирования в данном конкретном интеграле можно изменить и сначала проинтегрировать по $z$, а уж потом по $r$. В-третьих, при интегрировании по $z$ по бесконечным пределам можно смело отбросить слагаемое $-vt$. Оставшийся интеграл берется машиной.

$\int_{-\infty}^\infty \frac{z^2\,dz}{\left(r^2 + z^2\right)^5} = \frac{5\pi}{128r^7}$


В итоге получается ответ для полной мощности омических потерь

$P = \frac{15}{1024}\,\mu_0^2\,p_m^2\,\sigma\left(\frac{1}{a^3} - \frac{1}{b^3}\right)v^2 = k\,v^2$


Здесь после второго знака равенства мы обозначили коэффициент трения

$k = \frac{15}{1024}\,\mu_0^2\,p_m^2\,\sigma\left(\frac{1}{a^3} - \frac{1}{b^3}\right)$


Отметим что, коэффициент трения $k$ зависит только от намагниченности магнита $p_m$, свойств материала проводника $\sigma$ и геометрических размеров трубы $a$ и $b$ — то есть зависит исключительно от параметров магнита и трубы и не зависит от, например, скорости или времени. Это хороший знак для нас и маленький зачётик в копилку найденных формул! Отсюда же становится понятно, почему для демонстрации опыта выбрана именно медная труба, а не, скажем, стальная. Трение зависит от проводимости линейно $\sigma$, а у стали проводимость меньше на порядок.

А что если труба сделана из сверхпроводника?

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

Теперь можно записать

$mgv=kv^2\\mg = kv$


И внезапно (!), перед нами третий закон Ньютона! Сила действия равна силе противодействия. Можем найти установившуюся скорость движения магнита

$v_s = \frac{mg}{k}$

Уравнение движения


Настал черёд уравнения движения. С помощью второго закона Ньютона его будет записать очень просто

$ma=mg - kv\\ m\ddot{z} + k\dot{z} = mg$


Решать уравнение для $z(t)$ неинтересно, потому что ну просто координата меняется с постоянной скоростью. Гораздо полезнее знать, как быстро стабилизируется падение, чему равна установившаяся скорость падения. В общем, надо решать это уравнение для скорости

$\dot{v} + \frac{k}{m} v = g$


А решение будет такое

$v(t) = v_0\,e^{-\alpha t} + v_s\left(1 - e^{-\alpha t}\right)$


Здесь $\alpha = k/m$ — коэффициент затухания. Характерное время выхода на установившийся режим падения — $\tau = \alpha^{-1}$. Начальная скорость — $v_0$, установившаяся скорость — $v_s$.

А вообще, это уравнение парашютиста. Вот, наверное, почему статья Популярной Механики называется «Магнитный парашют».

Численный эксперимент


А теперь будет то, ради чего всё это затевалось. Навели тут, понимаешь, теорию. А на что она способна? Вдруг это всего лишь как тень на плетень? Или вообще не работает…

Начинаем снова с магнита. На сайтах интернет магазинов нашёл. Продаются магниты NdFeB, остаточная намагниченность $B_r = 10000 .. 13000$ Гс. Однако, это не означает, что поле непосредственно рядом с магнитом такое же большое. Когда-то доводилось измерять прибором. Поля прямо на поверхности таких магнитов обычно оказываются меньше и составляют 4000 — 5000 Гс. А масса такого магнитика будет равна примерно $m = 8$ г.

Видео у нас из MIT, стало быть, американское. Попробую угадать размеры их демонстрационной установки в дюймах (они же в дюймах любят всё измерять). Размер магнитика похож на $d = 1/2$ дюйма в диаметре. Это из тех какие есть в продаже. Размер медной трубы в длину похож на $l = 12$ дюймов (1 фут), а внутренний и внешний диаметры трубы, скорее всего, $2a = 5/8$ дюйма, $2b = 2$ дюйма.

С геометрией, вроде разобрались. Теперь физические свойства. Проводимость меди $59.5\times 10^6$ См/м. Ну и должен признаться, что как я ни пытался, ничего не смог найти на тему: как связаны остаточная намагниченность и магнитный момент постоянного магнита. Поэтому, держа в уме те измерения, которые когда-то делал, я просто подобрал такое значение $p_m$, чтобы поле на поверхности магнита было где-то между 4500 и 5000 Гс. Получилось $p_m \approx 0.6$ А·м². Чтобы убедиться, что поля похожи на настоящие, построил график магнитной индукции от радиальной координаты на расстоянии половины диаметра от центра шарика.

z-компонента магнитного поля на поверхности постоянного магнита
z-компонента магнитного поля рядом с поверхностью постоянного магнита

Теперь воспользуемся решением уравнения движения, чтобы построить график скорости магнита. Для всех выбранных выше параметров коэффициент трения получается равным $k = 1.032$ Н/(м/с), установившаяся скорость — $v_s = 7.6$ см/с — как раз примерно 3 дюйма в секунду! На видео шарик проходит через трубу длиной в 12 дюймов примерно за 4 секунды.

ЭТО ЗАЧОТ!

Знаю, что правильно «зачёт» писать через «ё», но в данном случае правильнее будет через «о» ;-)


А мы продолжаем. Рассеиваемая мощность оказывается равной примерно $P \approx 6$ мВт, а характерное время выхода на установившийся режим — $\tau \approx 8$ мс. Ниже построены графики $v(t)$ для двух разных начальных скоростей: нулевой, и $v_0 = 15$ см/с.

График решения уравнения движения магнитика в медной трубе
График решения уравнения движения магнитика в медной трубе

Вот такое вот получилось исследование. Надеюсь, что было интересно. Оставляйте ваши комментарии. Постараюсь ответить всем. Если вам понравилась статья, поддержите автора лайком или плюсиком в карму. Спасибо, что прочитали.

Литература


  1. Джексон, Дж. Классическая электродинамика: Пер. с англ. Мир, 1965.
  2. Ландау, Л. Д., & Лифшиц, Е. М. (1941). Теория поля. Москва; Ленинград: Государственное издательство технико-теоретической литературы.
  3. Сивухин, Д. В. «Общий курс физики. Том 3. Электричество.» Москва, издательство «Наука», главная редакция физико-математической литературы (1977).
  4. Яворский, Б. М., and А.А. Детлаф. «Справочник по физике.» (1990).

© Habrahabr.ru