CC BY
651
Изв. вузов «ПНД», т. 20, № 1, 2012
УДК 53.081.5(075.2)
ДВЕ ЛЕКЦИИ. АНАЛИЗ РАЗМЕРНОСТЕЙ ИЛИ РАЙСКАЯ ЖИЗНЬ В ФИЗИКЕ
Д.И. Трубецков
В двух лекциях, прочитанных на школе для старшеклассников «Нелинейные дни в Саратове для молодых», изложены основные определения и положения теории размерности и подобия, а также некоторые специальные приемы решения конкретных задач, в частности, использование я-теоремы и дополнения Хантли. Приведены решения задач из различных областей физики.
Ключевые слова: Размерность, подобие, модели.
Введение
Лекции посвящены, как следует из названия, анализу размерностей и его применению в физике. Вторая часть названия заимствована из статьи В. Вайскопфа «Современная физика в элементарном изложении», который писал:
«...Я буду вычислять только порядки физических величин. Используемый в расчетах знак равенства поэтому будет означать вовсе не "равно", а "такого порядка величины". Численными множителями, например, 2п и т.д. я буду пренебрегать. Райская жизнь!».
В его слова включается метод оценок (в большей степени) и анализ размерностей, о чем мы и будем говорить.
Тому же Вайскопфу принадлежит статья об удивительной эффективности математики в физике (замечу, что не только в физике). Мы привыкли к чуду абстрактных чисел, например, к их сложению: когда числа абстрактные, то все равно, что складывается, фрукты или галоши. Но ситуация меняется, если вы знаете, к какому объекту относится число. Вспомним детские стишки: «И вышло у меня в ответе два землекопа и две трети». В данной конкретной ситуации число 2| недопустимо, хотя в арифметике - допустимо.
Естественный вопрос: «Если мы имеем дело с размерными величинами, то что делать? И вообще, есть ли какая-нибудь наука на сей счет»? Да, есть, и ее называют анализом размерности или теорией размерности. Наука небольшая, но ее и возможные ошибки из-за ее неумелого использования должен знать, по-моему, любой школьник, изучающий естественные науки.
О предсказаниях и моделях
В основе теории любого процесса или явления лежит модель, справедливая лишь в определенных пределах. В каких именно, мы узнаем чаще всего тогда, когда она перестает описывать реальность. Нужно либо поправлять модель, либо искать новую. В этот период ломки модели и важна теория размерности. Чем бы вы ни занимались - теорией или экспериментом, важно не только поставить задачу (создать теоретическую модель или придумать постановку эксперимента), но и попробовать представить себе, что должно получиться в результате вашей работы. Вот несколько мнений на эту тему.
Английский физик-экспериментатор Г. Липсон пишет: «Находятся люди, которые считают, что у ученого не должно быть никакой предвзятой идеи относительно исхода его эксперимента; ученый - говорят они - должен быть совершенно объективным. Это вздор. Настоящий ученый почти всегда ставит эксперименты с целью проверить ожидаемые результаты. Он испытывает удовольствие, обнаружив то, что ожидал, и разочарование, когда результаты эксперимента не совпадают с ожидаемыми. Если же ученый ничего не ожидает получить, то он не может быть вполне уверен в значимости своих результатов».
А вот мнение выдающегося ученого и педагога Джона Уилера - учителя Ричарда Фейнмана, которое получило название «правило Уилера»: «Никогда не начинай вычислений, пока не знаешь ответа. Каждому вычислению предпосылай оценочный расчет: привлеки простые физические соображения (симметрию! инвариантность! сохранение!) до того, как начинать подробный вывод; продумай возможные ответы на каждую загадку. Будь смелее: ведь никому нет дела до того, что именно ты предположил. Поэтому делай предположения быстро, интуитивно. Удачные предположения укрепляют эту интуицию. Ошибочные предположения дают полезную встряску».
Повторюсь, что тому, кому придется серьезно заниматься физикой, нужно прежде всего научиться делать оценки, которые часто подсказывают и сам путь более точного решения задачи. Именно это помогает делать анализ размерностей -самый распространенный метод физических оценок. По мнению одного из самых выдающихся физиков ХХ столетия Энрико Ферми, «в физике нет места для путаных мыслей... действительно понимающие природу того или иного явления должны получать основные законы из соображений размерности».
А что значит - понимать природу того или иного явления? Это означает -создать его модель.
Модель (от лат. modulus - мера, образец, норма) - в широком смысле в науке - аналог, «заместитель» оригинала (фрагмента действительности), который при определенных условиях воспроизводит интересующие исследователя свойства оригинала. Модель - карикатура на изучаемое явление, но зато одну и ту же модель можно использовать для описания широкого круга явлений.
Иногда модели рассматривают как «картины», соотносящиеся с чем-то.
Академик Н.Н. Моисеев считал, что модель можно рассматривать как специальную форму кодирования информации, причем модель кодирует ту информацию,
которую люди раньше не знали. В этом заключена предсказательная способность модели.
Создание модели - главный момент исследований с помощью анализа размерностей.
1. Выражение единиц измерения произвольной физической величины через единицы измерения величин, принятых за основные, называется размерностью.
2. Размерность любой физической величины может быть только произведением возведенных в степень величин, принятых за основные.
3. Размерности обеих частей равенства, выражающего некоторую физическую закономерность, должны быть одинаковы.
4. Безразмерные комплексы размерных величин, представляющие собой произведение различных степеней этих величин, называются критериями подобия; их обычно обозначают как П = idem.
В этих достаточно прозрачных определениях неясным остается одно: какие величины принять за основные. Дадим ответ на этот вопрос.
Когда-то не считалось неудобным, что разные длины измерялись по-разному. Вот, к примеру, старорусские меры длины:
И только во времена Великой французской революции была разработана и принята метрическая система мер. Первый эталон метра - доля длины парижского меридиана, равного 40 ■ 106 м.1 Измерение меридиана, необходимое для создания метрической системы, велось с восьмидесятых годов восемнадцатого века и длилось много лет. Самой яркой фигурой в этой измерительной эпопее был, несомненно, французский ученый и человек удивительной судьбы Доминик Франсуа Араго. В замечательной книге Даниила Александровича Гранина «Сад камней»2 есть «Повесть об одном ученом и одном императоре», посвященная Араго. Начинается она так: «Имя Араго хранилось в моей памяти со школьных лет. Щетина железных опилок вздрагивала, ершилась вокруг проводника... Стрелка намагничивалась внутри соленоида. Красивые, похожие на фокусы опыты, описанные во всех учебниках...».
Думаю, что и вы помните эти опыты.
Но в жизни Араго кроме науки были переодевания, побеги, рабство, пираты... По Гранину, Араго - «...молодой француз, немножко д'Артаньян, немножко хитроумный Одиссей, влюбчивый, неприступный, любознательный и легкомысленный...», благополучно выходящий из самых невероятных ситуаций. Я не могу и не хочу пересказы-
1В настоящее время метр есть длина, равная 1650763.73 длины волны (в вакууме) излучения, соответствующего переходу между уровнями 2рю и 5(15 атома криптона-86. До 1960 года за метр было принято расстояние между двумя штрихами, нанесенными на платиноиридиевом стержне, который находился в Международном бюро мер и весов в Севре, близ Парижа.
2Гранин Д.А. Сад камней. М.: Современник, 1972, с. 125-177.
Основные определения теории размерностей
1 вершок -
4.445 см, 0.7112 м, 2.1336 м, 1.0668 км, 10.954 км2.
1 аршин - 16 вершков 1 сажень - 3 аршина -1 верста - 500 саженей 1 десятина -
вать повесть. Прочитайте ее. Еще только одна цитата из повести Гранина, имеющая отношение к нашей теме.
«У меридиана были свои законы. Участок, измеренный Лакайлем, давно слился с отрезком, измеренным Кассини, а тот - с отрезком, измеренным Мешеном, теперь к ним добавится дуга, измеренная Араго.
На меридиане не пишут имен, на нем не будут обозначены три года жизни Араго, его скитания, свист пуль, чума, там будут лишь градусы, минуты, секунды».
Наряду с единицей длины была введена и единица измерения времени - се-кунда3 - доля звездного (сидерического) года, который равен 3.1557 ■ 107 с. Метра и секунды оказалось достаточно, чтобы измерять то, что называют пространством-временем. Следовательно, кинематика по сути дела совпадает с геометрией пространства-времени. Напомню, что в школьных учебниках кинематику определяют как раздел физики, в котором изучается механическое движение тел без выяснения причин, вызвавших это движение. Если есть взаимодействие тел, то это уже динамика, и для изучения ее законов необходимо ввести понятие силы, а также связать силу, действующую на тело, с ускорением, которое приобретается телом, или связать силу тяготения с массой, которая порождает гравитационное поле.
Появилась «двуликая масса»: в одних построениях она была инертной, в других - гравитационной. Однако Галилей своими опытами доказал, что все тела падают с одинаковым ускорением. Масса оказалась только одна: «двуликость» исчезла, поскольку гравитационная и инертная масса совпали. (Впрочем, мы позднее вернемся к «двуликости», используя ее при решении одной из задач.) В этом совпадении состоит великий принцип эквивалентности. Массу стали определять как некоторую величину, характеризующую реакцию тела на заданную силу, как меру инертности. Был создан рукотворный эталон, равный 1 кг.
Так появилась система единиц длина-масса-время (ЬМТ). Указанные размерности независимые (из них нельзя составить безразмерную комбинацию) и именно их в анализе размерностей чаще всего принимают за основные. Исходя из физических законов, можно было все новые и новые единицы сводить к основным. Действительно, обозначая размерность любой величины квадратными скобками, получим, например, для скорости [V] = ЬТ-1, для ускорения [а] = ЬТ-2, для силы = ЬМТ-2, для энергии [Е] = МЬ2Т-2, исходя из закона Кулона для заряда [д] = Ь3/2М1/2Т-1 и т.д.
Как указывает Я. Смородинский4, в возможности создания эталонов - единственных на весь мир - скрыто глубокое удивительное свойство природы: при установлении эталонов длины, массы и времени человек использовал важное свойство нашего физического мира - свойство однородности пространства по расстояниям, по времени и по массе. Впрочем, последнее не совсем верно даже на Земле, поскольку 1 кг ведет себя по-разному на разных широтах.
Давайте теперь, следуя Я. Смородинскому, введем в изложение некоторую детективную интригу. Пусть случилась катастрофа: Парижская лаборатория мер и весов погибла. А нам нужно определить период малых колебаний маятника. Что делать? Вот здесь и поможет анализ размерностей, поэтому обратимся к определениям,
3Секунда - единица измерения времени, равная 1/31556925.9747 доли тропического года для момента 1900 г., января 0, в 12 часов эфемеридного времени (31 декабря 1899 г., полдень).
4Смородинский Я.А. Подобие в природе и фундаментальные постоянные // Наука и жизнь, 1977,
№ 8, с. 79-85.
которые мы сформулировали раньше. Уже из этих утверждений следует рецепт решения задач методом размерности.
• Выпишите группу N физических величин, между которыми, как вам кажется, есть какая-то взаимосвязь; это - самый важный и ответственный этап -этап создания модели, на котором нельзя не учесть главного в описываемом явлении, но нельзя ввести и лишние величины.
• Не забудьте рядом с ними поставить их размерности, выраженные через K<N размерностей величин, принятых за основные.
• Попытайтесь составить из выписанных величин безразмерные произведения; некоторые величины при этом, возможно, придется возвести в какие-то степени. Часто не так уж просто составить такие произведения без использования формальных приемов, поэтому можно воспользоваться определением 2 (это основная формула метода размерности), составив матрицу размерностей -таблицу, в которой по горизонтали расположены выбранные вами для решения задачи величины (а 1, а2 ... ап), а по вертикали - величины, выбранные как основные; в клеточках матрицы ставятся показатели степеней, в которые нужно возвести основные величины, чтобы получить размерности (а1, а2 ... ап).
• Если N — К = 1, то безразмерное произведение будет единственным и, приравняв его безразмерной константе, вы получите искомую закономерность. Вернемся к задаче о маятнике и на примере ее решения используем описанный
выше рецепт.
Рассмотрим для простоты математический маятник - тяжелую материальную точку массой т, подвешенную на нерастяжимой и невесомой нити длиной I. Поскольку маятник движется к положению равновесия под действием силы тяжести, его период т может зависеть от ускорения силы тяжести д, от массы т маятника и его длины I. Ограничим совокупность возможных движений маятника условием, что они плоские. Будем считать, что период колебаний маятника не зависит от амплитуды колебаний, то есть ограничимся малыми колебаниями около положения равновесия. Пренебрежем затуханием колебаний (в число определяющих величин не войдут вязкость и температура воздуха), не будем также учитывать ускорение точки подвеса маятника вместе с Землей (в число определяющих величин не войдет скорость вращения Земли). Тогда для периода колебаний можно записать формулу
т = / (1,т,д). (1)
Выберем в качестве основной систему ЬМТ. Поскольку N = 4, а К = 3, то N — К = 1 и безразмерное произведение будет единственным. В данном случае легко сообразить, что безразмерная комбинация имеет вид: П = л/т2/(7/д) или т = Пу/1/д. Однако решим задачу и формальным образом. Составим матрицу размерностей, исходя из того, что [т] = Т1; [I] = Ь1; [т] = М1; [д] = Ь1Т-2. Тогда
т 1т д
Ь 0 10 1
М 0 0 1 0
Т 10 0 —2.
Используя определение 2 и формулу (1), имеем
т = Шатвду, (2)
где а, в и у - не известные нам показатели степеней. На основе положения 3, из формулы (2) получим
Ь0М0Т1 = ЬаМ вГ'Т-2"'. (3)
Приравнивая показатели степеней в формуле (3) при одинаковых основаниях, приходим к системе уравнений
а + у = 0, 1 = -2у, в = 0,
из которой следует, что
11
у = —, а = -. г 2' 2
Таким образом,
Согласно определению 4, мы получили критерий подобия
П = ^ или П = П2 = ^. (4)
Дадим еще одно определение.
Явления называются подобными, если они отличаются только численными значениями определяющих параметров и притом так, что для них соответствующие безразмерные величины П¿(г = 1,2, ...m) совпадают.
Применительно к нашей задаче о маятнике закон подобия говорит, что, если мы увеличим длину маятника в n раз, то отношение длины нового маятника к квадрату его периода остается неизменным, поскольку li/xf = 12/т2 = g/П, а g = const. Таким образом, если существует хоть один маятник, для которого известна величина l/т2, то мы можем найти периоды маятников любой длины.
В большинстве случаев, прежде чем приступить к изготовлению какого-либо дорогостоящего и крупного сооружения, например, корабля или самолета, для получения его характеристик в предстоящих условиях работы прибегают к испытаниям на моделях - к моделированию. При этом надо знать, как пересчитать результаты опыта на модели на натуру. Если этого не знать, моделирование бесполезно. Для рационального моделирования основным является понятие подобных явлений и, следовательно, установление критериев подобия.
Разумеется, закон подобия, как и любые другие законы природы, должен иметь ограниченную область применения. Очевидно, что маятник, длина которого сравнима с радиусом R Земли, будет качаться не так, как маятник обычных размеров l, потому что натяжение, обусловленное силой тяжести, будет существенно меняться вдоль его длины.
Таким образом, закон подобия для маятника будет справедлив только в том случае, если l ^ R. Следовательно, геометрического подобия в гравитационном поле нет.
Мы уже упоминали, что Галилей открыл закон подобия для свободно падающих тел (а мы его использовали для маятника, считая g = const). Но по иронии судьбы он же отметил, что некоторые законы не остаются неизменными при изменении масштабов. Вот его размышления над прочностью костей. Животное вдвое
большей длины, ширины и высоты должно весить в 8 раз больше. Однако его вдвое более широкие кости имеют лишь вчетверо большее поперечное сечение, а значит, способны выдерживать лишь вчетверо больший вес. Это означает, что для выдерживания полного веса тела поперечный размер костей должен был бы увеличиться более чем в 2 раза. Такое отклонение от простого подобия вводит естественный масштаб в строение тела как сухопутных, так и водных животных: при некоторых предсказуемых в общих чертах размерах диаметр костей начинает увеличиваться быстрее, чем остальные части тела, разрушая тем самым подобие, но сохраняя в живых животных (на эту тему есть интересная книга: Дж.Б.С. Холдейн «Как важно быть нужного размера»5).
Перейдем теперь к решению разных задач.
Доказательство теоремы Пифагора. Историки науки считают, что развитию математических способностей Эйнштейна способствовал его дядя Якоб - инженер по образованию. Он давал мальчику математические задачи, и тот испытывал удовольствие от их решения. В 1891 году Эйнштейн приобрел, по его словам, «священную книгу по геометрии» - геометрию Эвклида. Изучая ее, он почувствовал, что некоторые доказательства в книге неоправданно сложны. Вот одно из возможных доказательств теоремы Пифагора, сделанное одиннадцатилетним Эйнштейном. Площадь прямоугольного треугольника Б определяется величиной его гипотенузы с и, для определенности, меньшим из острых углов ф, то есть Б = /(с, ф). Очевидно, что анализ размерности дает 5 = с2ф(ф).
Определение гравитационного радиуса Земли. Начнем с, казалось бы, абстрактной задачи на технику решения с помощью анализа размерностей.
Попробуем составить из гравитационной постоянной О, входящей в закон всемирного тяготения ^ = О(тхт2)/т2 (считаю, что вы его знаете), массы т и скорости света с величину, имеющую размерность длины. Обозначая ее К, имеем
Примеры различных задач, решаемых методами размерности и подобия
Из рисунка видно, что сумма площадей треугольников Б1 и Б2 равна Б. Тогда
с2ф(ф) = а2ф(ф) + Ь2ф(ф) и с2 = а2 + Ь2.
с
R = f (G, m, c) ^ R = ПGamecY.
(5)
Составим матрицу размерностей в системе LMT
R G m c
L 13 0 1
M 0-11 0
T 0 -2 0 -1,
5Haldane J.B.S. On being the right size. L.: Oxford University Press, 1928.
из которой с учетом формулы (5) следует:
Ь1М0Т0 = Ь3аМ-а Т-2аМ вЬ Т
Тогда 1 = 3а+у, 0 = —а+в, 0 = —2а—у и, следовательно, а = в = 1, У = -2. Таким образом, К = П(Ст/с2). Если т = МЗ - массе Земли, то Кгр = П(СМЗ/с2) - гравитационный радиус Земли. Важное замечание: величину П, которая присутствует во всех задачах, решаемых методом размерности, в рамках метода определить нельзя (вспомните Вайскопфа: «Я буду вычислять только порядки физических величин...»). Эта величина может быть определена либо из строгой теории, либо из эксперимента. В нашей задаче П = 2, поэтому Кгр Земли получается примерно равным 0.4 см, Кгр Солнца - около 3 км, Кгр пульсаров - порядка размеров Земли.
Для тела массой т, лежащего на поверхности Земли, тд = С(тМЗ/т2), поэтому д = Кгр с2/(2г2), то есть гравитационное поле само характеризуется постоянной Кгр, имеющей размерность длины. Вот почему в гравитационном поле нет геометрического подобия, и короткий и длинный маятники колеблются по-разному.
Задача для будущих президентов. В студенческой или школьной аудитории я всегда задаю вопрос: «Кто хочет стать президентом страны?» И всегда находятся претенденты. Недавно появилась любопытная книга Ричарда А. Мюллера «Физика для будущих президентов»6, в которой автор пытается ответить на вопрос: «Зачем президенту знать физику?». Ответ в том, что именно с физикой связано множество проблем, с которыми ему придется столкнуться (терроризм, энергетические проблемы, ядерные технологии, космос, глобальное потепление), и еще потому, что физика в наши дни становится главной наукой - наукой, открывающей человечеству двери в будущее. Естественно, что одна из проблем - ядерное оружие.
Поэтому - вот одна важная задача, которую с помощью анализа размерностей может решить президент, даже недолюбливающий физику. Решим задачу о сильном точечном взрыве, под которым понимают взрыв заряда, имеющего малый размер и массу, но обладающего большой энергией. Именно к таким взрывам и относится ядерный взрыв. Эту задачу решали почти одновременно Дж. Тэйлор в Англии, фон Нейман в США и академик Л.И. Седов в СССР. Вот изящное и простое решение, которое предложил Л.И. Седов.
Пусть в течение очень короткого времени в атмосфере выделилась настолько большая энергия Е, что при дальнейшем распространении ударной волны можно пренебречь атмосферным давлением за ударной волной. Найдем закон движения ударной волны, то есть ее расстояние г от центра взрыва в момент времени 1 Масштабы явления должны зависеть от мгновенно выделившейся энергии Е и инерционных свойств атмосферы - плотности р. Матрица размерности в системе ЬМТ имеет вид
Г Е р г
Ь 1 2 —3 0
М 0 1 1 0
Т 0 2 0 1
6Мюллер Р.А. Физика для будущих президентов. М.: АСТ: Астрель: Полиграфиздат, 2011, 411 с.
Будем искать зависимость г = /(Е, р,£). Тогда, согласно основной формуле теории размерностей (определению 2), г = ПЕ арв^. Используя матрицу размерностей, находим
Ь1М0Т0 = (Ь2МТ-2)а(МЬ-3)в Т*,
откуда 1 = 2а — 3в, 0 = а + в, 0 = —2а + у и а = 1/5, в = —1/5, у = 2/5. Таким образом, закон движения ударной волны имеет вид
г = ¿2/5. (6)
Соотношение (6) оказалось весьма полезным, поскольку экономило усилия и средства экспериментаторов. Действительно, достаточно произвести единственный взрыв с известными Е и р, измерить, до какого г дойдет ударная волна за время - и можно определить безразмерную постоянную П. После чего соотношение (6) можно использовать для определения энергии взрыва Е, если во время взрыва измерять г(£) и р. Разумеется, постоянная П (критерий подобия), должна быть одинаковой для всех взрывов, то есть взрывы должны быть однотипными.
В логарифмических переменных (5/2) ^ г, ^ £ экспериментальные точки должны лечь на прямую
55
2^=2 ^(ПЕ1/5р"1/5) + ^^, (7)
имеющую наклон, равный единице. Это подтвердил Тейлор, обработавший кинофильм о распространении огненного шара, снятый во время американских ядерных испытаний Дж. Маком. Из более строгих расчетов следует, что П ~ 1. Зная это, по экспериментальной зависимости радиуса фронта ударной волны от времени можно определить энергию взрыва. Публикация Тейлором этой величины (она оказалась примерно равной 1014 Дж) вызвала в свое время, по его словам, немалое смущение в американских правительственных кругах, поскольку эта цифра считалась весьма секретной, хотя фильм Мака секретным не был.
Проблема пульсации звезд и п-теорема. Пока у нас задачи решались легко, поскольку N — К = 1. А что делать, если N — К > 1? В теории размерности и подобия есть п-теорема, впервые сформулированная и доказанная, по-видимому Бэкингемом. Ее можно изложить следующим образом.
Пусть существует физическая закономерность а = /(а1, а2, ... , а^+1, ... ап), выраженная в виде зависимости некоторой, вообще говоря, размерной величины от размерных же определяющих ее параметров.
Эта зависимость может быть представлена в виде П = /(П1, П2,... Пп-&) некоторой безразмерной величины от безразмерных комбинаций определяющих параметров. Количество этих безразмерных комбинаций меньше общего числа определяющих параметров на число размерных определяющих параметров с независимыми размерностями (параметров, принятых за основные). Рассмотрим пример использования теоремы в астрофизике.
Для классификации звезд используют интенсивность их блеска. В 1784 году Джон Гудрайтис установил, что звезда в созвездии Цефея, которой в астрономических каталогах присвоена буква 8 (8-Цефея), периодически изменяет блеск с
периодом 5.366341 суток (порядка 130 часов). Ее диаметр изменяется в пределах от 40 до 45 диаметров Солнца. 8-Цефея дала имя целому классу звезд - цефеид. Цефеиды - переменные звезды-сверхгиганты, периодически изменяющие свою светимость и яркость. Причиной переменности служат колебания радиуса и температуры фотосферы (пульсации), возникающие из-за нарушения баланса между силами гравитации, стремящимися сжать звезду, и силами горячего газа и излучения, стремящимися звезду расширить. Существование связи между периодом и средней светимостью звезды установила в 1908 году Генриетта Суон Левит. В 1925-1926 годах, применив полученные Левит результаты к цефеидам, обнаруженным в туманности Андромеды, Эдвин Пауэлл Хаббл смог определить расстояние до нее и установить, что она расположена далеко за границами Млечного пути и является самостоятельной галактикой.
Как сформулировать задачу о пульсациях и какие величины выбрать в качестве определяющих?
Предположим, что пульсации звезд - их сжатия и разрежения - имеют строго выраженный радиальный характер. Амплитуду пульсаций будем считать малой. В качестве определяющих величин, следуя книге Э.А. Дибая и С.А. Каплана «Размерности и подобие астрофизических величин»7, выберем массу звезды М, светимость С и радиус П. Светимость - одна из важнейших единиц астрофизики. Под светимостью понимают количество энергии, излучаемой небесным телом за единицу времени, поэтому в системе ЬМТ [С] = Ь2МТ-3. Заметим, что М, С и П имеют независимые размерности, то есть из этих трех величин нельзя составить безразмерный комплекс. Кроме указанных определяющих величин следует добавить гравитационную постоянную Я, поскольку сжатие звезды при пульсации вызвано полем тяготения звезды. Наконец, при малых амплитудах пульсации определяются только периодом V. При этом отброшены такие «тонкости», например, как зависимость периода пульсаций от структуры звезды и ее химического состава, роль ионизации в поверхностных слоях, конвекция и т.п. Все эти факторы, конечно, влияют на период пульсаций, но это уже более высокое приближение в решаемой задаче.
Запишем на этот раз искомую физическую зависимость в неявной форме
/(М, С, П, Я, V)=0. (8)
Число независимых размерностей равно трем. Согласно п-теореме, должны быть два критерия подобия, П1 и П2, а уравнение (8) можно переписать следующим образом:
/1(П1, П2) =0.
Составим матрицу размерностей в системе ЬМТ.
М С п Я V
ь 0 2 1 3 0
М 1 1 0 —1 0
Т 0 3 0 2 1.
7Дибай Э.А., Каплан С.А. Размерности и подобие астрофизических величин. М.: Наука, 1976, 399 с.
Используя матрицу размерности и основную формулу размерности, найдем два безразмерных комплекса, в один из которых входит Р, а в другой - Я.
[V] = 1 т = 1
[С]а [М]в[П\1 ¿2а+уМа+рТ-За
2а + у = 0, а + в = 0, 1 + 3а = 0, а = -1/3, в = 1/3, у = 2/3. Таким образом,
П = {СГ)1/3 (9)
1 (мк2/г2у/з-
Есть ли какой-либо физический смысл у этого критерия подобия? Оказывается, что удобнее интерпретировать величину П^, которая есть отношение энергии СР, излученной звездой за один период колебаний, к величине, характеризующей механическую энергию колебаний звезды МП2/Р 2.
Аналогичным образом будем отыскивать и второй критерий подобия, который обозначим пока П^. Итак,
[Я] = 1 ь3м~1т~2 = х
[С]а1 [М]в1 [Я.] VI ' ¿2ах+ух М«1+в1 Т-За1
2а1 + У1 = 3, а! + в1 = -1, 3а1 =2; а! = -2/3, в1 = -5/3, л = 5/3 и, следовательно,
дс-2/3
11п
2 м~5' 3тг5/3'
Образуем критерий подобия
' р 2 п^ = п2п ¡ = д_ш_1п3.
Тогда
П2 = РЯ1/2М1/2П-3/2. (10)
Критерий П2 определяет период колебаний звезды, которые являются механическими. Вот почему мы используем вместо П2 критерий П2, который не содержит светимость звезды. Для средней плотности звезды можно написать соотношение р = М/((4/3)лЯ3). С учетом этого выражения из формулы (10) получим
Р(р)1/2 = "2 . (П)
7(4/3)^
Если считать правую часть формулы (11) постоянной, то (11) выражает хорошо известное в астрофизике соотношение «период-плотность» для пульсирующих звезд. А в чем физический смысл безразмерного комплекса П2?
Вспомним формулу Гюйгенса для периода колебаний математического маятника (формула (4) при П = 2п) и применим ее к пульсациям звезд, считая, что «длина» маятника I равна радиусу звезды Я. Кроме того, положим ускорение силы тяжести равным его значению на поверхности звезды, а именно д = ЯМЯ-2. Тогда
Г = 2лх1- = 2л\
д V МЯ
или с учетом выражения для р
V
Р(Р)1/2
(12)
Легко видеть, что соотношения (12) и (11) совпадают, если П2 = 2п. Откуда такая близость численных значений критериев подобия и их приближенное постоянство для объектов столь разных масштабов - маятника и звезды? Все дело в физической общности явлений, приводящих к малым пульсациям звезд и малым колебаниям маятника. И в том и в другом случае сущность одна и та же - механическое движение в поле тяжести, соответствующее небольшим отклонениям от положения равновесия.
Конечно, во многом сравниваемые явления различны. Они подобны только по одному критерию подобия (такое подобие явлений называется ограниченным). Более того, критерий подобия П1 не может считаться постоянным, он разный для разных звезд и всегда много меньше единицы.
Повторю еще раз, что численные значения критериев подобия из анализа размерности нельзя найти. Но очень часто они не отличаются от единицы. Почему? Одно из объяснений состоит в следующем: если анализируемая система находится в равновесном состоянии, то это значит, что два противоположно действующих эффекта (существенных для данной задачи) примерно компенсируют друг друга.
У формул размерности есть некоторые неудовлетворительные особенности. Например, угловая скорость и частота имеют одинаковую размерность. Формула размерности скорости изменения момента количества движения в системе ЬМТ имеет вид Ь2МТ-2, что совпадает с формулой размерности для энергии. Подобные примеры можно найти в книге Г. Хантли «Анализ размерностей»8.
Уже из двух приведенных примеров следует определенная двусмысленность системы единиц ЬМТ. Таким образом, если каждой физической величине можно приписать лишь одну вполне определенную формулу размерности, то обратное утверждение не верно: существуют некоторые формулы размерности, справедливые для более чем одной физической величины. Последнее приводит к неполному решению задачи, когда число независимых физических величин превышает число основных величин на два и более, поскольку в этом случае число неизвестных показателей степени превышает число уравнений, связывающих их друг с другом. В этом случае приходится выражать некоторые показатели через другие. Следовательно, было бы разумно увеличить число основных размерностей.
Рассмотрим следующую простую задачу. Пуля выпущена с начальной скоростью г0 в горизонтальном направлении на высоте Н от земной поверхности. Определить дальность К горизонтального полета пули. Очевидно, что К = /(г0,Н,д).
8Хантли Г. Анализ размерностей. М.: Мир, 1970, 175 с.
Дополнение Хантли
Тогда К = Пг^НУ; в системе ЬМТ имеем Ь = (ЬТ-1)аЬвЬ^Ти 1 = а + в + у; 0 = —а — 2у. Из двух последних уравнений у = -а/2, в = 1 — а/2 и, следовательно,
Я = Ш(-|=)<\ (13)
то есть решение неполное.
Дополнение Хантли состоит во введении векторных единиц длины с различающимися индексами, соответствующими декартовой системе координат, то есть Ьх, Ьу, Ьг. Этим способом система ЬМТ превратилась в систему ЬХЬУЬхМТ, в которой пять основных единиц вместо трех. В нашей задаче пуля обладает равномерной горизонтальной скоростью ЬХТ-1 и равномерным вертикальным ускорением ЬгТ-2. Используемый прием не только увеличивает число независимых основных размерностей, что дает возможность найти неизвестные показатели степени (в нашем случае а), но и устраняет путаницу в формулах размерности, когда одна и та же формула справедлива для двух и более физических величин.
Очевидно, что формулы площади ЬхЬу, ЬхЬ, ЬуЬ^ лучше, чем
Ь2.
Формулу, например, для давления Ь-1МТ-2 лучше преобразовать, вспомнив, что давление равно силе, действующей на единицу площади, то есть (МЬТ-2)/(ЬХЬу). Примеры можно продолжить.
Решим задачу о пуле заново, используя дополнение Хантли и считая, что [К] = Ьх, [го] = ЬхТ-1, [д] = Ь*Т-2, [Н] = Ь. Тогда К = П ■ г£Нвд^ и Ьх = (ЬХТ-1)аЬв(ЬгТ-2)^, что дает 1 = а, 0 = —а — 2у, 0 = в + У и, следовательно, у = —1/2, в = 1/2. Таким образом,
га
П = Пу0\ (14)
V д
Подставляя а = 1 в (13), также получим формулу (14).
С введением векторной величины длины все разумно и понятно, а можно ли разделить понятие массы? Ведь мы ранее, восхищаясь Галилеем, провозгласили принцип эквивалентности - масса одна. Но я обещал вернуться к этому вопросу. Следуя Хантли, решим такую задачу.
Определить массовый расход вязкой жидкости, протекающей через трубу круглого поперечного сечения радиусом г и длиной I. Массовый расход т - это масса жидкости, протекающей в единицу времени, то есть [т] = МТ-1. Так как течение жидкости поддерживается за счет разности давлений на концах трубы Р1 — Р2, то изменение давления на единицу длины (градиент давления) Р = (Р1 — Р2)/I и [Р] = МТ-2/Ь2. Вязкость жидкости п (динамическую) применительно к потоку, параллельному неподвижной поверхности, определяют как силу, действующую на единицу площади (которая перпендикулярна поверхности потока), и деленную на градиент скорости, который имеет место в направлении, перпендикулярном поверхности. Таким образом, в системе ЬМТ [п] = Ь-1МТ-1. Тогда уравнение задачи имеет вид: т = /(Р, р, п, г), где р - плотность жидкости, и, следовательно, в системе ЬМТ определяющих величин пять при трех основных.
Хантли предлагает вопреки общей теории относительности различать массу как количество вещества Ы^ и как меру инерции Ыг, то есть считать, что [ш] = ЫЦТ-1, [Р] = Ь-2ЫгТ-2, [р] = Ы,Ь-3, [п] = [Ь-1ЫгТ-1], [г] = Ь. Используя эти соотношения и формулу ш = ПРарвп^т8, получим
Ь°Ы^Т-1Ы° = (Ь-2ЫгТ-2)а (Ь-3Ы^ )в(Ь-1ЫгТ-1)* Ь8
и 0 = -2а - 3в - у + 8, в = 1, -1 = -2а - у, 0 = а + у. Откуда а = 1, у = -1, 8 = 4, и окончательно приходим к закону Пуазейля
Грт4
тп = П-^—. (15)
п
Задачу можно однозначно решить и в системе ЬМТ, если вместо ш ввести объемный расход жидкости V. Получится
Тт4
У = П-—. (16)
п
Если умножить (16) на р, то получится формула (15).
Задача Рэлея. Почему небо голубое? Тиндаль первым наблюдал, что белый свет при рассеянии становится синеватым и высказал мысль, что голубой цвет неба связан с рассеянием солнечного света на частичках пыли, которые всегда есть в достаточном количестве в атмосфере Земли. Не будь рассеяния света, небо было бы совершенно черным, и на нем ярко выделялись бы звезды и другие светила даже днем. Таким и видят небо космонавты в своих полетах. Но небо голубое... Вслед за Тиндалем Рэлей считал, что рассеяние света атмосферой объясняется наличием в ней взвешенных капелек жидкости, мелких пылинок и твердых частичек, размеры которых малы по сравнению с длиной световой волны. Он показал, что количественный результат можно получить, основываясь на анализе размерностей в сочетании с известными законами оптики. Решение Рэлея весьма поучительно и с точки зрения использования анализа размерностей. Приведем его. Рэлей формулирует задачу следующим образом.
Пусть частица с линейным размером I рассеивает солнечный свет с длиной волны X и амплитудой А. Амплитуда рассеянной волны уменьшается с увеличением расстояния от частицы. Пусть она равна Б на расстоянии т от частицы. Требуется определить зависимость Б от остальных переменных величин. Вся информация о сформулированной задаче собрана в таблице (используется система единиц ЬМТ).
Таблица
Физические величины и их размерности (к задаче Рэлея)
Физическая величина Обозначение Формула размерности
Амплитуда волны рассеянного света Ь
Амплитуда падающей волны А Ь
Линейный размер частицы 1 Ь
Расстояние от частицы г Ь
Длина волны света X Ь
По правилам теории размерностей выразим 5 как произведение остальных переменных, возведенных в ту или иную степень. Тогда получим, что
S = ПАа1МХ0. (17)
Из формулы (17) получается весьма необычное уравнение для размерностей
[L] = [Lja[L]e[L]Y[Ljô,
из которого следует, что
а + в + у + Ô = 1. (18)
Если теперь не дополнить анализ размерностей физическими представлениями, то задачу не решить. Вспомним, что амплитуда волны рассеянного света пропорциональна амплитуде падающего света, поэтому а = 1; амплитуда волны рассеянного света обратно пропорциональна расстоянию от частицы, то есть у = -1. Следовательно, возвращаясь к соотношениям (17) и (18), находим
A
0 = 1- 6, S = n~lpX1~p
r
или
Далее Рэлей отмечает: «Судя по динамике явления, N (отношение амплитуд волн падающего и рассеянного света) изменяется пропорционально V (объем рассеивающей частицы)». Таким образом, следует считать, что N = S/A ~ Iе ~ V ~ I3 и в = 3. Окончательно из формулы (19) имеем
S = —.
rX2
По определению, интенсивность I рассеянного света пропорциональна квадрату его амплитуды S, поэтому
I ~ ^ ~ «Л (20)
где ю - круговая частота световой волны. Если принять, что Xкрасная
~ 1.5ХСиняя, то 1синяя ^ 51красная, то есть синий цвет
рассеивается намного сильнее, чем красный.
Хантли заключает изложение решения Рэлея следующими словами: «Читатель может заметить, что в этом примере в большей мере использованы физическая интуиция и знание законов физики, чем анализ размерностей. С этим приходится согласиться, но верно также и то, что использование обоих источников привело простым и изящным образом к интересному результату, который невозможно получить методами элементарного анализа».
Результат Рэлея легко объясняет голубой цвет неба. Действительно, из-за присутствия атмосферы значительная часть солнечного света рассеивается в стороны. Причем, согласно формуле (20), рассеянная часть излучения тем больше, чем короче длина волны. Это значит, что рассеянный свет богат короткими волнами. Отсюда и голубой цвет неба, поскольку максимум интенсивности в рассеянном свете попадает в голубую область спектра.
Все было бы хорошо в этих объяснениях, если бы они были еще и правильными. Ошибка в исходном допущении о загрязнении атмосферы. Ведь по мере поднятия над Землей в воздухе уменьшается содержание посторонних частиц и пыли: воздух чище над океанами, чем в больших городах, он чище в Антарктиде, чем в Москве. Но тогда насыщенность рассеянного света в этих местах должна была бы уменьшаться. Однако научные наблюдения, например, в высокогорных обсерваториях показали, что все обстоит как раз наоборот: чем чище воздух, чем меньше в нем посторонних частиц, тем ярче небесная лазурь, тем богаче излучение неба синими лучами. При этом формула (20) остается справедливой и в этих случаях (закон был количественно подтвержден измерениями интенсивности рассеянного света).
Исходя из подобных фактов, Рэлей позднее понял, что рассеяние вызывается не посторонними частицами, а самими молекулами воздуха. Такое рассеяние света стали называть рэлеевским или молекулярным рассеянием.
Правда, на этом история ответа на вопрос «Почему небо голубое?» не заканчивается, хотя формула (20) сохраняется. Но это - другая история9.
На этом я заканчиваю лекции, и в заключение приведу цитату из книги Р.А. Мюллера: «Для понимания физики современной жизни важно изменить те свои представления, которые не соответствуют истине. Часто цитируют высказывание, приписываемое Марку Твену: "Неведение не приносит большинству людей никаких неудобств. Проблема в их обширных знаниях, которые таковыми не являются"». Думаю, что эти лекции пополнили ваши знания.
Тем, кого заинтересовала теория размерности и подобия, могу порекомендовать уже упомянутую мою книгу «Введение в синергетику» (естественно, что в лекцию вошло многое из нее), книгу Хантли, книгу Дибая и Каплана, а также следующие издания:
1. Брук Ю.М., Стасенко А.Л. Как физики делают оценки - метод размерности и порядки физических величин // О современной физике - учителю. М.: Знание, 1975. C. 54-131.
2. Седов Л.И. Методы подобия и размерности в механике. М.: Наука, 1977.
3. Бриджмен П. Анализ размерностей. Ижевск: НИЦ «Регулярная и хаотическая динамика», 2001, 148 с.
4. Журнал «Квант» за разные годы.
Саратовский государственный Поступила в редакцию 24.12.2011
университет им. Н.Г. Чернышевского
9См., например: Трубецков Д.И. Введение в синергетику. Колебания и волны. Изд. 4-е. М.: Книжный дом «ЛИБРОКОМ», 2012, глава 5.
TWO LECTURES. DIMENSIONAL ANALYSIS OR PARADISE IN PHYSICS
D.I. Trubetskov
In the two lectures read at school for senior pupils «Nonlinear days in Saratov for youth» the basic definitions of the theory of dimension and similarity are stated, and also some special ways for concrete problems decision, in particular, use Pi-theorem and additions of Huntley are presented. The solutions of problems from different areas of physics are shown.
Keywords: Dimension, similarity, models.
Трубецков Дмитрий Иванович - родился в Саратове (1938). Окончил физический факультет Саратовского государственного университета (1960). Защитил диссертации на соискание ученой степени кандидата (1965) и доктора физико-математических наук в СГУ (1978) в области радиофизики. Заведующий кафедрой электроники, колебаний и волн факультета нелинейных процессов СГУ, профессор, член-корреспондент Российской академии наук, заслуженный деятель науки РФ, лауреат премии Президента РФ в области образования. Научный руководитель Лицея прикладных наук и факультета нелинейных процессов СГУ. Область научных интересов: вакуумная электроника и микроэлектроника сверхвысоких частот, теория колебаний и волн, нелинейная динамика, история науки. Автор более двадцати учебных пособий и монографий, а также более двухсот статей в периодической печати.
410012 Саратов, ул. Астраханская, 83
Саратовский государственный университет им. Н.Г. Чернышевского E-mail: TrubetskovDI@nonlin.sgu.ru
