CC BY
63
УДК 514.7
РЕШЕНИЯ УРАВНЕНИЯ БЮРГЕРСА С ПЕРИОДИЧЕСКИМ ВОЗМУЩЕНИЕМ НА ГРАНИЦЕ
A.B. САМОХИН
Изучена асимптотика решений уравнения Бюргерса с начальными/граничными данными на конечном интервале с периодическим возмущением на границе. Уравнение описывает вязкую среду и первоначальный постоянный профиль переходит в бегущую волну с убывающей амплитудой. При малых значениях вязкости асимптотический профиль имеет пилообразный профиль с периодическими разрывами производной, похожий на известное решение Фэя на полупрямой.
Ключевые слова: уравнение Бюргерса, начально-граничная задача на отрезке, пилообразные решения.
1. ВВЕДЕНИЕ
Уравнение Бюргерса
ut (X, t) = e2uxx (x, t) - 2 • u (x, t )ux (x, t). (1)
описывает колебания вязкой среды. Вязкость гасит колебания (за исключением инвариантных относительно какой-либо подалгебры симметрий).
Начально-граничные данные на конечном интервале для уравнения Бюргерса имеют вид:
u(x, 0) = f (x), u(a, t) = l(t), u(b, t) = r(t), x e [a, b]. (2)
Нас интересует поведение решений в случае периодического возмущения на границе следующего вида:
u(x, 0) = a, u(0, t) = a + b sin(&t), u(L, t) = a, x e [0, L].
Для полубесконечного интервала x e [0, вопрос о виде таких решений хорошо
изучен. При значительной вязкости колебания экспоненциально затухают при (пространственном) удалении от источника возмущений. Однако во многих задачах нелинейной акустики возмущение u периодично по x, и нелинейные эффекты концентрируются вблизи источника и содержат там регулярно разнесённые в пространстве разрывы. Таким образом, при незначительной вязкости доминируют нелинейные эффекты и возникает убывающий пилообразный профиль, который постепенно расплывается (разрывы превращаются в скачки с увеличивающейся толщиной), превращаясь в затухающую волну - иногда на значительном удалении от источника. Достаточно адекватная асимптотика для таких решений предложена Фэем [].
Задача для конечного интервала, которая рассматривается в настоящей статье, приводит к аналогичным результатам, хотя решения имеют некоторые характерные особенности.
Работа является продожением [6], [7]. Численные расчёты проводились с использованием пакета Maple PDEtools.
2. ПИЛООБРАЗНЫЕ РЕШЕНИЯ УРАВНЕНИЯ БЮРГЕРСА
Как известно, для уравнения переноса общего вида ( x e R )
W + f (w) wx = 0
(3)
момент начала градиентной катастрофы может быть найден следующим образом. Пусть w = р( x) - начальный профиль. Решение задачи (17) может быть представлено в параметрической форме w = р(%), x = £ + F(£)t, где F(£) = f (р(4)) . Характеристики x = <% + F(£)t будут пересекаться в случае р'(^) < 0, что приводит к многозначности w (опрокидыванию волны или градиентной катастрофе). Если неравенство выполняется на бесконечном интервале, то существует минимальное значение времени tc, при котором
возникает градиентная катастрофа. Можно определить tc при помощи формулы tc = —1 / F'(^c)
где | F) |= max | F'(£) | на интервале [a, b], для которого F'(£) < 0.
Мы продемонстрируем, что градиентные катастрофы наследуются уравнениями типа Бюргерса на конечном интервале при некоторых начальных данных в случае слабой диссипации, добавленной к модели (17); см. [1-4].
Для полубесконечного интервала и периодического возмущения в точке Х0 = 0 вида:
u(0, t) = u0 + a sin( cot)
асимптотика решения имеет вид:
_ a ^ sin(w#) U = R ti sh[n(1 + X )/2R]'
здесь R - число Рейнольдса, 1^=c(t — x / u0) и X - безразмерная координата, Fay, [4],. График
решения для ut(x,t) = £2uxx(x,t) — 2 • u(x,t)ux(x, t), u(0,t) = 1 + 2sin(2^t) представлен на рис. 1.
Можно заметить два эффекта, вносимых слабой вязкостью среды (диссипацией). Во-первых,
амплитуда решения достаточно быстро убывает, и решение стремится к константе u0 . Во-вторых,
эффекты нелинейности (пилообразность, т.е. периодические разрывы производной) концентрируются вблизи источника возмущений и быстро сходят на нет за счёт всё той же вязкости.
JC
Рис. 1. £ = 0.2, Г = 30
В более сложной ситуации конечного интервала появляются некоторые новые эффекты. В частности, асимптотическое значение (или его среднее по времени, в том случае, когда стабилизация не достигается) будет отличаться от величины исходного возмущения и0 на правом конце интервала.
3. РАЗВИТИЕ ПИЛООБРАЗНЫХ РЕШЕНИЙ УРАВНЕНИЯ БЮРГЕРСА
В этом разделе обсуждается процесс развития пилообразных решений уравнения Бюргерса
Ч =£2ихх - 2иих (4)
со следующими граничными условиями на интервале [0, Ь]:
и ($, 0) = 1 + И\п(2жг), — и (Г, Ь) = 0, и (0, х)1 = 1.
дх [0,Ь]
(5)
В рассмотренном примере £ = 0.5, Ь = 75.
На рис. 2-5 показан процесс установления периодического (повремени) профиля, начиная от стационарного. Выбранная вязкость довольно высока и подобие пилообразности возникает только вблизи источника возмущения на левом конце интервала. После (пространственного) затухания периодического возмущения процесс протекает в форме движения ударной волны постоянной высоты А вправо: как показано ниже, её величина равна
1 С
среднему значению — I и (0, ?.
Т •
Рис. 2Л = 0.5
1
Ил
Г \
---1-'—
2,0 1.8 1.6 1.4 11 1.0 0.80.6 0.4 02
10 20 30 40 50 60 70
х
Рис. 3. t = 5
10 20 30 40 $0 60 70
X
Рис.4. / = 10
Г
1
О 10 20 30 40 50 60 70 х
Рис. 5. / = 15
Сравним эти графики с тем, что получается при значительно меньшей диссипации. При £ = 0.05 и прежних начально граничных условиях, профиль соответствующего решения представлен на рис. 6. Убывание амплитуды идет с меньшей скоростью, но пилообразный характер графика сохраняется вплоть до правого конца интервала. Фронт распространения возмущения представляет собой ударную волну.
/[/[АААЛЛЛЛАЛ^ЛЛ
10 2(1 30 40 50 60
Рис. 6. ^ = 25. £ = 0.05
10 20 30 40 50 60 70
Рис. 1.1 = 25, е = 0.2
На рис. 7 видно, что передний фронт зубьев пилообразного решения постепенно расширяется от почти вертикального вблизи источника возмущения, перед тем как перейти в гладкий профиль на достаточном удалении; в этом случае
1(г, 0) = 2 + 3вт(2П), —и($,Ь) = 0, и(0, х)1 = 2, £ = 0.15.
дх [0,Ь]
4. ОЦЕНКИ СКОРОСТИ РАСПРОСТРАНЕНИЯ И ВЫСОТЫ УДАРНОЙ ВОЛНЫ
Уравнение Бюргерса может быть записано в форме закона сохранения
(и\ = (£2их - и2)х,
поэтому контурный интеграл по границе прямоугольника ^ = {0 < X < Ь,0 < 1 < Т} равен нулю:
ф[и •Жх + (е2их - и2) • Ж] = 0.
дЕ
Принимая во внимание начально-граничные условия, можно переписать это следующим образом:
Ь 1 0 0
Iи(х, 0)Жх +1 (е2их (Ь, 1) - и2 (Ь, 1)) & +1и(х, Т)Жх +1 (е2их (0,1) - и2(0,1)) (И = 0.
Вследствие диссипации при достаточно больших Ь и T решение на правом конце интервала приобретает постоянное значение А, которое совпадает со средним значением. Меняя порядок интегрирования и учитывая, что Ux (Ь,1) = 0 по условию, получим:
Ь T Ь T
| u (x,0)dX +1 (е2их (Ь,1) - u 2(Ь,1))ё1 -1 u (x,T)dx -1 £4(0,1) - и 2 (0,1)) & =
0 0 0 0
T
= u0Ь - Л2Т - ЛЬ -1(£2UX(0,1) - u2(0,1))ё1 = 0,
ИЛИ
и0,
0Ь - Л2Т - ЛЬ -1(£2их (0,1) - и2 (0,1)) & = 0. (6)
0
Поделив на Т, получим квадратное уравнение для определения А:
Т
-Л2 - лЬ - Т|£2их(0,1) - и2(0,1))ё1 + ч0 Ь = 0. (7)
Займёмся оставшимся интегралом, который представляет собой среднее значение подынтегральной функции. Поскольку и(0, /) = и0 + а ), их (0, /) = а]008(]), то
1 Т £ Т 1 Т
Т |(^и^с (0,1) - и2(0,1)) ё1 = Т^а^] ооб])Ж - ^ {(и0 + а Бт]))2 Ж.
Поскольку среднее по периоду для синуса и косинуса равны нулю, при больщих Т с большой точностью получим, что искомый интеграл равен:
а2 т а2 Т а2 а2
и--[Бт2(]= -ч02--[(1 - ооб(2]))& ~ -и02--Т = -и---
т - 2Т - 2Т 2
Квадратное уравнение (27) приобретёт теперь вид:
2 Т и Т
*2 2 а . Ь Ь ^
Л + ч02 +--Л— + Ч0 — = 0. (8)
Отсюда Л = (-к + ^/к2 + 4ч- + 2а2 + 4ч0к ) /2 , где к = Ь / Т .
Отметим, что при Т ^^ получим Л = ^Ч^ + , а ПРИ Ь ^^ получим, в согласии с
результатами для полубесконечного интервала, Л = Ч0.
Для рассчитанных примеров и(0, /) = 1 + 2 ) ^ Л ~ V12 + 22 / 2 ~ 1.7 при Т ^^ .
Перейдём к определению скорости распространения возмущения (или скорости ударной волны). Известно, что для нелинейных волн скорость определяется амплитудой начального возмущения, в рассматриваемом случае, высотой первого горба (см. рис. 2). Эта скорость ~ 3 сохраняется, как это видно из рис. 4 - 6 и далее.
Для аналитического определения высоты первого горба можно также воспользоваться тем, что его форма хорошо описывается решением Р. Хохлова [4], [5]:
а
и =-
1 + X
th
kRB_ 1 + X
a
здесь R =-г - число Рейнольдса, $ = CO(t — x / u0) и X - безразмерная координата
2азе
5. ЗАКЛЮЧЕНИЕ
Численные расчёты позволяют не только обнаруживать новые эффекты, возникающие при решении уравнения Бюргерса на конечном интервале, но и экспериментально проверять асимптотические оценки, полученные аналитически. Отметим, что численное моделирование функций с разрывной производной является непростым делом, поскольку вблизи разрывов стандартные методы теряют устойчивость. Потеря устойчивости приводит к мульти-осцилляциям и потере точности, не говоря уже о потере ясности. С этой проблемой удалось справиться за счёт адаптивной длины шага по пространственной координате (уменьшая шаг в 10-20 раз по сравнению со значениями по умолчанию).
REFERENCES
1. Dubrovin B., Elaeva M. On critical behavior in nonlinear evolutionary PDEs with small viscosity // ArXiv: 1301.7216v1math-ph., 30.01.2013, 16 p.
2. Dubrovin B., Grava T. and Clein C. Numerical study of breakup in generalized Korteweg de Vries and Kawahara equations // Siam J. Appl. Math, 71: 4 (2011), pp. 983-1008.
3. Dubrovin B. On Hamiltonian Perturbations of Hyperbolic Systems of Conservation Laws, II: Universality of Critical Behaviour // Comm. Math. Phys., 267 (2006), pp. 117-139.
4. Fay R.D. J.Acoust. Soc. Am., Proc., 3, 1931, pp. 222-241.
5. Rudenko O.V. Nonlinear sawtooth-shaped waves // UFN, 9 (1995), pp. 1011-1035 (in Russian).
6. Samokhin A., Gradient catastrophes for a generalized Burgers equation on a finite interval // Geometry and Physics, Elsevier, the Netherlands, 85 (November 2014), pp. 177-184
SOLUTIONS TO THE BURGERS EQUATION WITH PERIODIC PERTURBATIONS ON BOUNDARY
Samokhin A.V.
The asymptotic behavior of solutions of the Burgers equation with initial value - boundary problem on a finite interval with periodic boundary conditions is studied. The equation describes a dissipative medium, so a constant initial profile will evolve to a travelling-wave solution. Its asymptotic limit is periodic 'sawtooth' solution with periodical breaks of derivative, similar to the Fay solution on a half-line.
Keywords: Burgers equation, initial value - boundary problem, finite interval, sawtooth solutions.
Сведения об авторе
Самохин Алексей Васильевич, 1947 г.р., окончил МГУ им. М.В. Ломоносова (1971), доктор технических наук, профессор кафедры высшей математики МГТУ ГА, автор 40 научных работ, область научных интересов — уравнения математической физики, симметрии, законы сохранения.
