книги2 / 10-1
.pdfЗдесь L = π (2/β)1/2, а функции a(t), b(t) удовлетворяют системе дифференциальных уравнений
|
a˙ = αa + β a2 |
+ |
b2 |
, |
b˙ = αb + |
3 |
β ab . |
(3) |
||
|
|
|
|
|
||||||
|
|
2 |
2 |
|||||||
Семейство |
решений |
w(x, t) |
полностью |
описывается: |
a(0), b(0) — начальными данными решений системы (3) и координатой x0.
Будем считать, что решение u(x, t) обладает компактным носителем, начальный размер которого r меньше L . Тогда найдутся такие: точка x(1)0 и значения параметров a1(0), b1 > 0, |a1| b1, для которых имеет место неравенство u(1)(x, 0) u(x, 0). Выбором па-
раметров x(1)0 , a1(0), b1(0) среди совокупности всех допустимых для них значений можно добиться, чтобы функция u(1)(x, 0) аппроксимировала функцию u(x, 0) снизу наиболее оптимальным образом.
Точно также построим эталонное решение u(2)(x, t) с набором па-
раметров a2(0), b2(0), x(2)0 таким образом, чтобы выполнялось неравенство u(2)(x, 0) u(x, 0), но его носитель не превосходил L . Вы-
брав параметры aj(0), bj(0), x(0j), мы, тем самым, зафиксировали ре-
шения aj(t), bj(t) , j {1, 2} динамической системы a(t), b(t) , для которых выбранные значения являются начальными данными. В ре-
зультате, аппроксимируемое точное решение уравнения (1) подчинено, в силу Теорем 2.2, неравенствам u(1)(x, t) u(x, t) u(2)(x, t).
Пусть это решение обладает обострением режима с временем обострения t < ∞. Заметим, что при α 0 любое эталонное решение обладает обострнием режима. Если же α < 0, и, в этом случае, параметры a−(0), b−(0) могут быть выбраны так, что эталонное решение u(1)(x, t) также обладает обострением режима с неко-
торым временем обострения t(1), то, в силу указанного неравенства, t(1) t . Кроме того, обострением режима обладает эталонное решение u2(x, t) с временем обострения t(2), которое удовлетворяет неравенству t(2) t .
В силу неравенства u(1)(x, t) u(x, t) u(2)(x, t), зависящий от времени размер носителя r(t) решения u(x, t) также подчинен нера-
венствам r1(t) r(t) r+(t) L , где rj(t) — размеры носителей эталонных решений j {1, 2}. Поэтому, имеет место r(t ) L . Из
Теоремы 3 следует, что предельные значения всех носителей совпадают и равны L . Таким образом, мы получаем возможность оценивать время обострения и размер области локализации произвольного
81
решения уравнения (1) с компактным носителем, не превосходящим
L .
Система (3) интегрируется. Это позволяет построить эталонные
решения uj(x, t), j {1, 2}, у которых b(t) > |a(t)|. На основе этих решений находятся главные члены асимптотик решений при
u(x, t) → ∞. Эти асимптотики характеризуются временем обострения
t = α−1 ln |
1 + βb(0) |
1 − b(0) |
|
|
− |
× |
|
|
|
|
|
|
|||||||||
|
|
α |
|
|
|
a(0) |
|
2 |
|
3/2 |
|
|
|
|
|
|
|
||||
|
|
|
a(0) |
2 |
|
|
1/2 |
|
|
a(0) |
|
|
a(0) |
2 |
|
1/2 |
|
||||
× harcsin h1 − |
|
|
i |
+ |
|
|
h1 |
− |
|
|
|
i |
|
i, |
|||||||
b(0) |
b(0) |
b(0) |
|
выраженным через начальные данные.
Литература
1.Маслов В.П. Математическое моделирование процессов тепломассопереноса. Эволюция диссипативных структур / В.П. Маслов, В.Г. Данилов, К.А. Волосов.// — М. : Наука, 1969. — 352 с.
2.Потемкина Е.В. Режимы с обострениями в задаче Коши для неоднородного уравнения теплопроводности/ Е.В. Потемкина // Успехи мат. наук. — 1996. — Т. 51, № 6. — C. 221–222.
АСИМПТОТИКИ ДЛИННЫХ НЕЛИНЕЙНЫХ БЕРЕГОВЫХ ВОЛН ИХ СВЯЗЬ С БИЛЛИАРДАМИ
СПОЛУЖЕСТКИМИ СТЕНКАМИ1
М.М.Вотякова, Д.С. Миненков
(Москва, ИПМех РАН; МГУ, механико-математический факультет) votiakova.mm@phystech.edu, minenkov.ds@gmail.com
Под береговыми волнами мы понимаем периодические или близкие к периодическим по времени гравитационные волны на воде в бассейне глубины D(x), x = (x1, x2), локализованные в окрестности береговой линии Γ0 = {D(x) = 0}. В двух конкретных примерах мы строим отвечающие береговым волнам асимптотические решения системы нелинейных уравнений мелкой воды в виде параметрически заданных функций, определяемых через асимптотики линеаризованной системы (см. [1]), которые в свою очередь связаны с асимп-
тотическими собственными функциями оператора ˆ − .
L = gD(x)
1 Работа выполнена при финансовой поддержке гранта РНФ 22-71-10106 в МГУ имени М.В.Ломоносова
© Вотякова М.М., Миненков Д.С., 2024
82
Область определения оператора — гладкие функции ξ(x) в области Ω = {x : D(x) > 0} с конечной энергией: |ξ|x Γ0 < ∞. Также обсуждается связь построенных асимптотик с классическими (почти интегрируемыми) "биллиардами с полужесткими стенками”.
Литература
1. Dobrokhotov S.Y. Asymptotic Solutions of the Cauchy Problem for the Nonlinear Shallow Water Equations in a Basin with a Gently Sloping Beach / S.,Y. Dobrokhotov, D.S. Minenkov, V.E. Nazaikinskii // Russ. J. Math. Phys. — 2022. — vol. 29, — p. 28–36.
ЗАДАЧИ ВОЗНИКАЮЩИЕ ПРИ ИССЛЕДОВАНИИ НЕКОТОРОГО КЛАССА НЕЛИНЕЙНЫХ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
И МЕТОДЫ ИХ РЕШЕНИЯ М.В. Гасанов, В.Н. Орлов (Москва, НИУ МГСУ)
GasanovMV@mgsu.ru, Orlovvn@mgsu.ru
При описании процессов и явлений очень часто используются дифференциальные уравнения. В случае линейных дифференциальных уравнений, существует классическая теория для их решения, но когда при описаний явлений и процессов используются нелинейные дифференциальные уравнения, то возникают сложности, связанные с наличием подвижных особых точек, которые являются условием неразрешимости таких уравнений, в общем случае, в квадратурах. На данный момент, для решения нелинейных уравнений существует всего два варианта: либо разрешимость в квадратурах, с помощью использования специальной замены переменной (преобразование Шварца, преобразование групп Ли и т. п.), как показано в работах [1], [2], либо аналитический приближенный метод решения [3 – 5], основанный на решении следующих математических задач:
1.Теорема существования и единственности решения в области аналитичности и окрестности подвижной особой точки
2.Построение структуры аналитического приближенного решения в области аналитичности и окрестности подвижной особой точки.
3.Влияние возмущения начальных данных на структуру аналитического приближенного решения в области аналитично-
©Гасанов М.В., Орлов В.Н., 2024
83
сти. Влияние погрешности приближенного значения подвижной особой точки на структуру аналитического приближенного решения в окрестности подвижной особой точки
4.Расширение границ области применения аналитического приближенного решения в окрестности приближенного значения подвижной особой точки
5.Необходимые, необходимые и достаточные условия существования подвижной особой точки
6.Алгоритм нахождения подвижной особой точки с заданной точностью
Рассматривается задача Коши |
|
y′′′ = y2(z) + r(z). |
(1) |
y′(z00) = y01, |
(2) |
y(z ) = y , |
|
y′′(z0) = y2.
Теорема 1. Потребуем выполнение следующих условий:
1.Пусть z – подвижная особая точка решения задачи Коши (1) –
(2)
2.Функция r (z) голоморфна в области |z − z| < ρ1,
тогда решение задачи Коши (1) – (2) имеет вид:
y (z) = (z − z)−3 |
|
∞ |
Cn (z − z)n |
(3) |
|||||||||||
X |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
||
в области |
|
|z − z| < ρ2, |
|
|
|
|
|||||||||
где |
|
|
|
|
(4) |
||||||||||
ρ2 = min ρ1, |
|
|
1 |
|
, |
|
|||||||||
|
√6 |
|
|
||||||||||||
(|y0 |
M |
+ 1 |
|
||||||||||||
| , |y1| , |y2 |
| , |
n |
( |
|
|
n! |
|
)) |
|
||||||
|
|
|
|
|
|
r(n)(z ) |
|
|
|
||||||
M = max |
|
|
sup |
|
|
|
|
|
|
|
|
, n = 0, 1, 2, . . . |
|
||
|
|
|
|
|
|
|
|
|
|
|
|
Для определения критериев существования воспользуемся технологией регуляризации подвижной особой точки, сделав замену y (z) = w(1z) .
84
Инверсная задача Коши имеет вид: |
|
w′′′w2 = −6w′3 + 6ww′w′′ − w2 − r (z) w4 |
(5) |
w′ (z00) = w01, |
(6) |
w (z ) = w , |
|
w′′ (z0) = w2.
Решение данной задачи Коши, с учетом замены, имеет вид:
∞
w (z) = (z − z)3 XDn (z − z)n .
0
Критерии существования в комплексной и вещественной области отличаются. В случае комплексной области используется переход к фазовым пространствам.
Далее будут сформулированы критерии для случая вещественной области. В формулировках теорем будем использовать независимую переменную x.
Теорема 2. (точечный критерий существования подвижных особых точек) Точка x является подвижной особой точкой, кратности 3, функции y(x), решение задачи Коши (1)–(2), тогда и только тогда когда для функции x (w), x (w) — есть обратная функция к решению инверсной задачи Коши (5)–(6), выполняются условия:
x (0) = x , x′ (0) = 0, x′′ (0) = 0, x′′′ (0) = |
|
1 |
. |
(7) |
|
10 |
|||||
|
|
|
Теорема 3. (интервальный критерий существования подвижной особой точки) x является подвижной особой точкой решения задачи (1)–(2) тогда и только тогда, когда существует некоторая окрестность подвижной особой точки [x1; x2], x [x1; x2], для которой функция w (x) являлась бы непрерывной, и выполнялось условие:
w(x1) · w(x2) < 0.
Литература
1. Carillo S. Schwarzian derivative, Painleve XXV–Ermakov equation, and Backlund transformations,/ S. Carillo, A. Chichurin, G. Filipuk, F. Zullo // Math. Nachr. — 2023. — P. 1–19.
https://doi.org/10.1002/mana.202200180
2. Carillo S. A short note on the Painleve XXV–Ermakov equation / S. Carillo, A. Chichurin, G. Filipuk, F.
85
Zullo // Applied Mathematics Letters, Vol. 131, — 2022.
https://doi.org/10.1016/j.aml.2022.108064
3. Orlov V. Moving Singular Points and the Van der Pol Equation, as Well as the Uniqueness of Its Solution. / V. Orlov // Mathematics —
2023. https://doi.org/10.3390/math11040873
4.Pchelova A. Construction of approximate solutions for a class of first-order nonlinear differential equations in the analyticity region, Herald of the Bauman Moscow State Tech. Univ. / Pchelova // Nat. Sci., — 2016. — P. 3–15. DOI: 10.18698/1812-3368-2016-3-3-15
5.Gasanov M. A Study of a Mathematical Model with a Movable
Singular Point in a Fourth-Order Nonlinear Differential |
Equation/ |
M. Gasanov, A. Gulkanov // Rus. J. Nonlin. Dyn., |
— 2023. |
https://doi.org/10.20537/nd230904 |
|
ОПЕРАТОРНЫЙ ПОДХОД К ПЕРИОДИЧЕСКИМ РЕШЕНИЯМ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ
НА ГРУППАХ ЛИ1
Ю.Е. Гликлих (Воронеж, ВГУ) yeg@math.vsu.ru
К настоящему веремени множество процессов в технике и инженерном деле описывается математически в терминах обыкновенных дифференциальных уравнений на нелинейных многообразиях. И очень часто возникает задача о нахождении периодических решений указанных уравнений. Эта задача особенно сложна, если правая часть уравнения только непрерывна и не удовлетворяет, например, условию Липшица, т.е. задача Коши для уравнения не имеет единственного решения. В этом случае для уравнений в линейном пространстве для исследования периодических решений используется метод интегральных операторов. Напомним, что обыкновенное дифференциальное уравнение в векторном пространстве может быть преобразовано в эквивалентное интегральное уравнение. Например, задача Коши x˙ = f t, x(t) , x(0) = x0 в Rn эквивалентна интеграль-
ному уравнению x(t) = x0 + R0t f τ, x(τ) dτ.
Однако классические интегральные операторы на многообразиях не ковариантны, т.е., зависят от выбора карты. Ранее нами были построены так называемые интегральные операторы с римановым параллельным переносом, которые ковариантны, что позволило
1 Работа выполнена при финансовой поддержке РНФ (проект № 24-21-00004). © Гликлих Ю.Е., 2024
86
использовать их для исследования дифференциальных уравнений с непрерывной правой частью, но которые были не применимы в задаче о периодических решениях. Здесь мы строим новый тип указанных операторов на группах Ли такой, что их неподвижные точки являются периодическими решениями дифференциальных уравнений с непрерывными периодическими правыми частями. Отметим, что эти операторы действуют в пространстве (многообразии) C1- кривых и только из вторые итерации являются вполне непрерывными операторами. Мы описываем конструкцию топологического индекса для этих операторов и не основе использования направляющих функций приводим пример, в котором индекс оператора не равен нулю откуда следует существование периодического решения.
ЗАДАЧА ХОДЖА-ГЕЛЬМГОЛЬЦА И УСЛОВИЕ
ПРИЛИПАНИЯ ВО ВНЕШНИХ ОБЛАСТЯХ А.В. Горшков (Москва, МГУ) alexey.gorshkov.msu@gmail.com
Задачей Ходжа-Гельмгольца называется задача восстановления соленоидального векторного поля по вихревой функции. Рассмотрим внешнюю область Ω R2, которая является дополнением односвязной области с кусочно-гладкой границей. С условием прилипания и заданным потоком на бесконечности внешняя задача имеет вид
div v(x) = 0, |
(1) |
curl v(x) = w(x), |
(2) |
v(x) = 0, x ∂Ω, |
(3) |
v(x) → (v∞, 0), |x| → ∞. |
(4) |
Здесь v = (v1, v2) - двумерное векторное поле, curl v(x) = ∂x1 v2 −
∂x2 v1.
С условием прилипания эта задача вообще говоря не разрешима. Для её разрешимости требуется ряд дополнительных условий на вихревую функцию.
В работе Квартапэлли, Уальс-Гриза[1] выведены проекционные условия на вихревую функцию, заключающиеся в ортогональности ротора гармоническим функциям. Векторное поле представляется как орто-градиент v = ψ функции тока ψ, где x = (−∂x2 , ∂x1 ).
© Горшков А.В., 2024
87
Тогда в ограниченных областях с нулевыми граничными услови-
ями ψ(x) = ∂ψ(x) = 0, x ∂Ω из соотношения ∆ψ = curlv согласно
∂n
формуле Грина следует ортогональность curlv гармоническим функциям. Это и есть условие прилипания, выраженное в терминах только вихревой функции.
Оказывается, этот принцип ортогональности может быть распространен на внешние области[2]. Аналогичное условие ортогональности для внешних областей имеет следующий вид:
2π |
ZΩ Φ(z)k dx = |
(iv∞, |
k = 1.{ |
} |
̸ |
(5) |
||
1 |
|
w(x) |
0, k |
N |
0 |
|
, k = 1, |
|
Здесь Φ — это отображение Римана области Ω на единичный круг B1. Для функции ротора будем рассматривать весовое пространство
Z
L2,N (Ω) = {f(x) : f(·) 2L2,N (Ω) = |f(x)|2(1 + |x|2)N dx < ∞}.
Ω
Теорема Если w L2,N (Ω) с N > 1 удовлетворяет (5), то тогда существует единственное решение задачи (1)-(4), удовлетворяющее оценке
v(·) − v∞ H1(Ω) C w L2,N (Ω).
Как было показано Финном[4], для стационарных течений величина v(·) − v∞ имеет бесконечную L2 норму. Однако из (5) следует, что среднее ротора равняется нулю, и, следовательно, согласно формуле Стокса, отсутствует циркуляция на бесконечности, и энергия становится конечной.
Условие ортогональности (5) позволяет найти граничное условие для ротора curl v(t, x) нестационарной системы Стокса. Пусть w(t, x) = curl v(t, Φ−1(x)) есть вихревая функция потока Стокса, заданная во внешности круга. Для нее граничное условие прилипание в терминах коэффициентов Фурье wk будет иметь вид
|
∂wk(t, r) |
r=r0 |
+ |k|wk(t, r0) = 0, k Z. |
r0 ∂r |
|||
|
|
|
|
|
|
|
|
Используя эти граничные условия в работе автора [2] была найдена явная формула решения системы Стокса обтекания кругового цилиндра в вихревой форме. В работе [3] найдено схожее граничное условие для ротора, соответствующее условию прилипания, для системы Навье-Стокса.
88
Литература
1.Quartapelle L. Projection conditions on the vorticity in viscous incompressible flows. / Quartapelle L., Valz-Gris, F. // International Journal for Numerical Methods in Fluids. — 1981. — 1(2), — pp. 129– 144.
2.Gorshkov A.V. No-Slip Boundary Condition for Vorticity Equation in 2D Exterior Domain / Gorshkov A.V. // J. of Math. Fluid Mech. — 2023. — 25:47.
3.Gorshkov A.V. Associated Weber-Orr Transform, Biot-Savart Law and Explicit Form of the Solution of 2D Stokes System in Exterior of the Disc. / Gorshkov A.V. // J. Math. Fluid Mech. — 2019. 21:3.
4.Finn R. An energy theorem for viscous fluid motions./ Finn R. // Arch. Rational Mech. Anal. — 1960. 6 — pp. 371–381.
ТУННЕЛЬНЫЕ ПО МАСЛОВУ АСИМПТОТИКИ ДЛЯ ЗАДАЧ С НЕПРЕРЫВНЫМИ И ДИСКРЕТНЫМИ
АРГУМЕНТАМИ В.Г. Данилов (Москва, НИУ ВШЭ)
vgdanilov@mail.ru
Туннельные асимптотики это асимптотические решения имеющие вид
Uas = exp(−S/ε)φ(x, t, ε), S const,
ε → +0 – малый параметр. Естественное приложение для такой техники – параболические задачи различного происхождения. Принципиальное отличие асимптотических решений такого вила от обычных ВКБ решений состоит в том, что exp(−S/ε) играет роль "второго"малого параметра. Это проявляется в том, что оценка max |U − Uas| = O(εN ) – фиксированное число >0, вполне хорошая для асимптотических ВКБ решений, оказывается малосодержатель-
ной в области S δ > 0, где Uas = O(εM ) для всех M > 0, то есть асимптотика меньше разности между асимптотическим и точ-
ным решением. Далее, в отличие от ВКБ решений, к туннельным задачам с трудом применимы техники, основанные на преобразовании Фурье. Это затрудняет построение фундаментальных решений и разрушает конструкцию канонического оператора Маслова, широко используемого для гиперболических задач.
© Данилов В.Г., 2024
89
В докладе я расскажу о малоизвестной (моей) формуле для асимптотических фундаментальных решений задач Коши в пространстве и на сетке, о возможных способах построения глобальных асимптотических решений (включая случай с каустиками). Также будет рассказано о применении этой техники к задачам случайного блуждания.
Литература
1.Danilov V. G. A representation of the delta-function via creation operators and gaussian exponentials, and multiplicative fundamental solution asymptotics for some parabolic pseudodifferential-equations / V. G. Danilov // Russian Journal of Mathematical Physics. — 1995. — Vol. 3. — no. 1. — P. 25–40.
2.Danilov V. G. Nonsmooth non-oscillating exponential-type asymptotics for linear parabolic PDE / V. G. Danilov // SIAM Journal on Mathematical Analysis. — 2017. — Vol. 49. — no. 5. — P. 3550–3572.
3.Данилов В.Г. Асимптотика фундаментальных решений параболических задач / В. Г. Данилов // Математические заметки. — 2024. — Т. 115. — В. 2.
ДВУМЕРНЫЕ ИНТЕГРАЛЬНЫЕ ОПЕРАТОРЫ С ПОДВИЖНЫМИ И ФИКСИРОВАННЫМИ
ОСОБЕННОСТЯМИ ПО ОГРАНИЧЕННОЙ ОБЛАСТИ Г. Джангибеков., Г. Козиев (Душанбе, ИM НАНТ и МУТПТ) gulkhoja@list.ru, gulnazar88@mail.ru
Пусть D – конечная односвязная область комплексной плоскости, ограниченная простой замкнутой кривой Ляпунова Γ и содержащая внутри точку z = 0; D = D Γ; a(z), b(z), c(z), d(z), e(z) - непрерывные в D функции; K1(z, ζ), K2(z, ζ) - измеримые ограничен-
ные функции, имеющие пределы limz,ζ→0 Kj(z, ζ) = Kj(0, 0); h1(σ) и h2(σ) - измеримые на всей комплексной плоскости функции, причем
ZZ
|hj(σ)||σ|−βdsσ < ∞, j = 1, 2,
|σ|<∞
где β – некоторое число из интервала (0, 2), dsσ – элемент плоской меры Лебега; пусть, наконец, B(z, ζ) – керн-функция Бергмана области D, представимая в виде (см., например, [1], стр. 252, 258)
B(z, ζ) = |
1 |
|
ω′(z) |
ω′(ζ) |
|
|
, |
|
|
|
|
|
|
|
|||
|
π (1 − ω(z)ω(ζ))2 |
© Джангибеков Г.,Козиев Г., 2024
90