Свешников, Боголюбов, Кравцов. Лекции по математической физике (1125176), страница 49
Текст из файла (страница 49)
Оно называется функцией Римана. Зная функцию Римана, легко построить решение задачи (11.12) — (11.14). Из формул (11.17) — (11.20) с учетом (11.12) получается 306 Непосредственной проверкой легко устанавливается, что для операторов 7. и К выполняется равенство (11.15), где ф(А)о(А)+~р(В)о(В)+~ ( )~( ) 1 ! 2 в (11.21» — — 1 сг' с1х т(',7 с(у 2 .) АВ причем последний интеграл в формуле (11.2!) легко вычисляется, поскольку функции о, ер и ф известны. 3 а м е ч а н и е.
В начале этого параграфа мы отметили, что любая характеристика уравнения (1!.!2) должна пересекать кривую С не более одного раза. Действительно, если характеристика пересекает кривую С в двух точках А и Мь то значение и(М,) не может быть задано произвольно, а определяется по формуле и(А) о(А)+и(В,) о(В,) у и (М,)— 2 + в, — 1 '1е(хс(у — — ( Удх+б1с(у 12,) в Л с начальным значением, заданным на дуге АВ, и функцией 1(х, у), заданной в области х)ы— криволинейном треугольнике М,В,А (рнс. 7.13). Рис. 7.13 5 12. НЕЛИНЕЙНЫЕ УРАВНЕНИЯ ю 1. Простейшие уравнения и метод характеристик Рассматривая простейшие уравнения колебаний с постоянным коэффициентом а иа бесконечной прямой: и„=пи„„хе= йх, Г)О, (12.1) мы установили, что его общее решение представляется в виде суперпозиции прямой и обратной бегущих волн: и(х, 1) =1',(х — а!)+)а(х+а(), (12.2) где 1, и )т — пРоизвольные дважды непРеРывно диффеРенцируемые функции своих аргументов.
Как легко проверить, пря- маЯ ио1=1,(х — а1) и обРатнаЯ и1'>=Гх(х+а1) волны Удовлетворяют следующим уравнениям в частных производных первого порядка: ю Этот параграф написан профессором С. А. Габовым Более подробно о рассматриваемых уравнениях смл Габон С. А. Введение в теорию нелинейных волн Изд-во МГУ, 1988. 307 (12. 4) или — !(и,+ии„)1'Я) =О. Так как и(х, 1) — решение уравнения (!2 6), то данное соотношение удовлетворяется при любой дифференцируемой функции !($). Отсюда следует, что решение исходной задачи Коши (12.6), (12.7) определяется из неявного уравнения и (х, !) = и, (х — и (х, !) ! ).
(12. 9 им> — аи!и = О, Уравнения (12.3) и (12.4) часто называются уравнениями переноса, поскольку их решения (функции 1,(х — аг) и 1,(х+ат)) переносят соответственно начальный профиль 7, (х) или !2(х) вдоль оси х с постоянной скоростью а в прямом или обратном .направлении. Легко показать, что если решение уравнения переноса является дважды непрерывно дифференцируемой функцией своего аргумента, то оно удовлетворяет и уравнению ко.лебаний. При этом не происходит искажения профиля волны во времени, т.
е. отсутствует эффект дисперсии волны, что связано с постоянством скорости распространения, в свою очередь определяющимся постоянством параметров среды. Однако во многих физических задачах приходится учитывать изменение свойств среды под действием распространяющихся в ней волн. Это приводит к зависимости скорости распространения от решения, что связано с необходимостью рассматривать квазилинейное уравнение переноса вида и, +,'Г (и) и, = О. (12.5) В этом параграфе мы рассмотрим свойства решения простейшего квазнлинейного уравнения переноса вида и,+ии„=О, хе=к', !)О. (!2.6) Поставим для него задачу Коши, задав начальное условие и (х, 0) =- и, (х) .
(12.7) По аналогии с линейным уравнением переноса с постоянным коэффициентом будем искать решение уравнения (!2.6) в виде и (х, г) = 7 (х — и (х, !) 1) (12. 8) н попытаемся найти функцию 7 такую, чтобы удовлетворялись .как уравнение (12.6), так и начальное условие (12.7). Вычисляя частные производные и„и иь получим и„= (! — и„1) )' ($), и,=( — и — и,г) 7'я), $=х — и!. Подстановка этих выражений в уравнение (!2.6) дает ! — и — и!) )' (~) + ( и — и иД 1' $) = 0 звв Мы не будем заниматься анализом соотношения (12.9), а попытаемся получить необходимые в дальнейшем результаты другим способом — методом характеристик.
Рассмотрим на плоскости (х, 1) кривую, определяемую дифференциальным уравнением Ш ах 1 и(х, 1) (12. 10) илн й !=а 3 Ф ч 'а Рис. 7.!5 Рис 7.14 Проведем на рисунке горизонтальную прямую (=х, и опреде- 309 — =и(х, г), йх ш и называемую характеристикой уравнения (12.б).
Пусть х=х(1) — решение уравнения (12.10), тогда ах — и (х ((), () = и, + и„— = и, + ии„! „=„! ~! = О, ш Ж т. е. и(х(1),1) есть константа на кривой х=х(1) и, следова- тельно, как вытекает из (12.10), х=х(1) представляет собой прямую линию на плоскости (х, 1) с наклоном и(х(1), 1) =и(х(0), О) =иа(х(0) ), определяемым начальной функцией иа(Ц, ~=х(0).
Для этой прямой можно записать уравнение (12. 11) ис (41 Чаким образом, мы получаем однопараметрнческое семейство прямых, зависящих от параметра ~ и обладающих тем свойст- вом, что решение и(х, () уравнения (!2.6) на этих прямых ока- зывается постоянным. Это позволяет по начальной функции иа(я) определить функцию и(х,1) в любой момент времени й Покажем, как это можно сделать практически. Пусть начальная функция ис имеет вид, изображенный на рис. 7.14, н равна нулю вне интервала (а, б). Выберем неко- торую точку ахен!и, Ь) и построим соответствующую ей харак- теристику Гг„.
х=-ах+(их($х) с углом наклона 1д<р=1(исДх) на рис. 7.15. Всюду на этой характеристике и(г =иа(ах). 1и ! 1 ! а(Х, ~ 2. Обобщенное решение. Условия на разрыве Для того чтобы исключить неоднозначные решения, оказывается необходимым расширить понятие решения уравнения (12.6) и вместо непрерывно днфференцируемых решений рассматривать разрывные. При этом, разумеется, необ- з!о лим точку пересечения этой прямой с характеристикой Г! .
Пусть эта точка имеет координаты (х,,1,). Изобразим теперь на плоскости (х, и) точку с координатами (хх, и, Я„) ) . Проведем аналогичные построения для различных точек ~~~[а, 61, получим совокупность точек (хм из($,)), образующую некоторую линию на плоскости (х, и), представляющую график решения и(х,1) уравнения (12.6) для момента времени 1=1, (рнс. 7.16). С помощью этой процедуры мы можем вычислить функцию и(х,1) для любых моментов времени. Отметим, что поскольку скорость переноса начального значения и,Я~) вдоль соответствующей характеристики 1 11 1 зависит от этого решения, то ! начальный профиль и,(х) с и(х,й/ и(х,г,) течением времени искажает- 1 ся — имеет место явление ! дисперсии бегущей волны.
0 Обратимся вновь к рис.7.!5 и заметим, что, начиная с неРяс 7.Щ которого момента времени 1= =1, т. е. при всех 1>1, характеристики, отвечающие различным ~„начинают пересекаться, что приводит к тому, что при 1>1 профиль решения и(х,1) оказывается неоднозначным, т, е. одному значению х могут отвечать два и более значения функции и(х,1) (см. рис. 7.16). Таким образом, здесь мы сталкиваемся с явлением, носящим название опрокидывания волн.
Понять причину этого опрокидывания легко, если обратиться к представлению решения (!2.6), (12.7) в виде (12.9). Из этого представления следует, что чем выше амплитуда точки, тем с большей скоростью она распространяется. Поэтому точки вершины волны обгоняют в своем движении точки подошвы волны. Заметим теперь, что возникновение неоднозначного профиля решения оказывается, как правило, противоречащим сути физической модели, описываемой уравнением (12.6), согласно которому и(х, 1) — однозначная функция. Например, если мы рассматриваем волны в сплошных средах: в одной точке физические параметры не могут иметь различные значения. Выходом из сложившейся ситуации является расширение понятия решения задачи (12.6), (12.7). з(олимп придать новый смысл выражению «функция и(х, () удовлетворяет уравнению (12.6)».
Естественным здесь является введение обобщенных решений так, как они вводятся в теории обобщенных функций. Напомним о п р е д е л е н и е. Будем говорить, что функция (((х, 1) удовлетворяет уравнению (12.6) в обобщенном смысле, если для любого прямоугольника П,(=((х, !): х)<х<х„0<11< <(<Г,) и любой бесконечно дифференцируемой и финитной в П,( функции ()) (х, !) справедливо интегральное тождество ~ и)р(+ — ((»(р„~~ ((х ((( = О, ! 2 (12.!2) пх! Из (!2.12) следует, что если обобщенное решение (!2.6) непрерывно дифференцируемо, то оно удовлетворяет (12.6) в обычном смысле, поскольку, интегрируя (!2.12) по частям, по- лучим ~ (и(+ ии,) ()) ((х ((( = О. пм ( и)р(+ — 'иф„1 ((х(((= ('!)рсоа(п"!) и( — ) +)рсоа(п х) " 1'~ ((з, „1 2 "1 .1! 2 т(«( (1) (и+)~1 (и)г(+ — и'))),(((х(((= — ( ))))соз(п"!)игь)+(рсоа(п х) ) ~ (Ь.
2 ),11 2 и(" 5 ы Здесь учтено, что функция и(х, () в П,'(' (/г=1, 2) удовлетворяет уравнению (12.6), и приняты следующие обозначения: соз(п" !) и соз(л х) — направляющие косинусы вектора нор- з(! Отсюда в силу произвольности П„, и функции (()(х, !) будем иметь (12.6). Тем самым данное выше определение обобщенного решения действительно является расширенным понятием классического решения. Применим введенное понятие обобщенного решения для установления простейших свойств разрывных решений уравнения (!2.6).
Пусть и(х, () — разрывное решение, имеющее единственный разрыв на кривой, описываемой уравнением х=з(1) и представляющей собой на плоскости (х, !) линию 5=((х, !): хс з(()Е Рассмотрим произвольный прямоугольник П„(, содержащий линию 5. Предположим также, что в частях Пм и П„, (1) 12) на которые линия 5 делит прямоугольник П (, решение и(х, !) непрерывно дифференцируемо, т. е. удовлетворяет (! 2.6) в обычном классическом смысле (рис.
7.!7). Преобразуем формулу (12.12) отдельно в частях П,( и П„(, (1) (2) Интегрируя по частям, получим мали п=(соз(п 1), соз(п" х)), изображенного на рис. 7.17; индексы «+» и « — » обозначают предельные значения функции и(х, 1) на кривой 5 при стремлении к ней «справа» и «слева» от кривой. Складывая полученные формулы, имеем ф '[соз(п 1)[и]+сов(п х) [ — ~~ с(з=О, (12.13) 2 где [и]=ие — и —. Учитывая произвольность функции лр(х,)), из равенства (12.13) легко выводим, что «2 1 соз(п 1) [и]+сов(п"х) ~ — ]! =О. (12.14) ()о ]5) Формула (12.!5) связывает скорость о„=з(1) распространения разрыва и значения решения ие и и- слева и справа от разрыва. Она называется формулой Гюгонио — Ренкина, нли формулой условий на разрыве.
Знание этой формулы позволяет определить скорость распространения разрыва по значениям и= решения на разрыве, но, однако, не дает ответа на вопрос о положении разрыва х=з(1). Мы не будем останавливаться на детальном решении этого вопроса, а укажем (без обоснования) лишь способ построения разрыва. Основная идея, на которой основано это по. строение, опирается на следующий факт. Уравнение (!2.6), записанное в виде (12.16) имеет вид одномерного уравнения неразрывности или закона сохранения.
![](https://s.studizba.com/z.php?f=/templates/si/i/noavatar.png&w=100&h=100&t=1"){kind=avatar}
