Сферически симметричные решения нелинейного уравнения Гельмгольца

peregoudov
Сообщений: 620
Зарегистрирован: 29 дек 2015, 13:17

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение peregoudov » 14 апр 2021, 21:45

Утащил с dxdy, где ничего путного не написали и даже немного наговорили глупостей. А задачка мне показалась интересной.

Найти сферически-симметричные решения уравнения

$$ \Delta u+u^n=0. $$

$\Delta$ --- трехмерный лапласиан. Пусть n будет натуральным и нечетным. В линейном случае n=1, например, решениями являются $\frac{\sin r}r$ и $\frac{\cos r}r$. Будут ли решения при других n похожи на эти?

Аватар пользователя
Ian
Сообщений: 960
Зарегистрирован: 18 янв 2016, 19:42

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение Ian » 15 апр 2021, 17:52

[math]
[math]
Именно при [math] имеется однородное решение [math], проверяемое подстановкой. Про другие решения не думал

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение zykov » 15 апр 2021, 22:43

При $n=3$ это решение вырождается в $u(r)=0$.

Аватар пользователя
Ian
Сообщений: 960
Зарегистрирован: 18 янв 2016, 19:42

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение Ian » 16 апр 2021, 11:35

[math]
[math] случай с n=3 наиболее мутный

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение zykov » 16 апр 2021, 12:58

На википедии такое есть: Lane–Emden equation.

peregoudov
Сообщений: 620
Зарегистрирован: 29 дек 2015, 13:17

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение peregoudov » 16 апр 2021, 15:43

Ну вот и на dxdy наибольшие терки были по поводу n=3... Википедия --- это не спортивно. На мой взгляд, тут довольно элементарными методами можно исследовать качественное поведение решений при любых n. Первый шаг Ian уже фактически сделал.

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение zykov » 17 апр 2021, 00:22

На вики написано, что аналитические решения есть только для $n=0,1,5$.
Там видимо имеется ввиду граничное условие $u(0)=u_0$ (конечное), $u'(0)=0$.
(Без потери общности при $n\neq 1$ можно положить $u_0=1$ просто изменяя масштаб $u$ и $r$.)
Решения Ian в нуле стремятся к бесконечности при $n>1$, так что не подходят.

Так что, можно уточнить, что имеется ввиду?
peregoudov писал(а):Source of the post Найти сферически-симметричные решения уравнения

Аналитических там больше нет.
В виде ряда легко представить. Нечётные производные равны нулю, $u''(0)=-u_0^n/3$ и дальше производные 4ой, 6ой и т.д. степеней в нуле находятся из того же уравнения.

Аватар пользователя
Ian
Сообщений: 960
Зарегистрирован: 18 янв 2016, 19:42

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение Ian » 17 апр 2021, 21:37

Замена независимой переменной [math] упрощает уравнение. А что узнать-то хотели

[math]
Последний раз редактировалось Ian 18 апр 2021, 05:22, всего редактировалось 1 раз.

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение zykov » 18 апр 2021, 01:49

Вроде при $n=3$ при больших $r$ оно похоже на $$\frac{\sin(2^{1/4} \sqrt r)}{2^{1/4} \sqrt r}$$.
В нуле он $$1-\frac{{{r}^{2}}}{6}+\frac{{{r}^{4}}}{40}-\frac{19 {{r}^{6}}}{5040}+\frac{619 {{r}^{8}}}{1088640}-\frac{17117 {{r}^{10}}}{199584000}+\frac{1208293 {{r}^{12}}}{93405312000}+\mbox{...}$$.

Аватар пользователя
Ian
Сообщений: 960
Зарегистрирован: 18 янв 2016, 19:42

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение Ian » 18 апр 2021, 08:37

Кажется. кроме n=3 они осциллирующие
https://www.wolframalpha.com/input/?i=s ... %29%29%5E5
А при n=3 осциллируют не все
https://www.wolframalpha.com/input/?i=s ... %29%29%5E3


peregoudov
Сообщений: 620
Зарегистрирован: 29 дек 2015, 13:17

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение peregoudov » 19 апр 2021, 16:51

Я думаю, следует прежде всего ставить задачу качественного исследование решений. Если этого не сделано, то и разложение построить не получится. А качественно исследовать мы хорошо умеем автономные уравнения...

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение zykov » 21 апр 2021, 08:07

zykov писал(а):Source of the post Вроде при $n=3$ при больших $r$ оно похоже на $$\frac{\sin(2^{1/4} \sqrt r)}{2^{1/4} \sqrt r}$$.

Если сделать замену переменных $t = 2^{1/4} \sqrt r$ и $u(r) = \frac{f(t)}{t}$, то уравнение будет $f''+f'/t+2f^3-f/t^2=0$.
Если предположить, что $f(t)$ - что-то вроде синуса, амплитуда которого не слишком быстро убывает (чтобы $f^3$ падал не слишком быстро), то можно оставить только $f''+2f^3=0$.
Тут решение конечно не синус (из-за куба), но на него похоже.
Тоже что-то периодическое с постоянной амплитудой.
(Нелинейные колебания на пружине, где возвращающая сила пропорциональна кубу сдвига. Потенциальная яма 4-ой степени.)

Нестационарная сила линейного отталкивания $f/t^2$ даёт малый вклад. Возможно как-то несильно и повлияет на амплитуду.
Для линейной по скорости силы трения $-f'/t$ коэффициент трения линейно падает по времени.
За один период она отъест $\overline{(f')^2} \frac{T}{t}$ энергии.
Т.е. по идее энергия должна падать как $t^{-1}$ (тут правда в отличии от линейного осцилятора $\overline{(f')^2}$ может быть не пропорционален энергии - нужно проверить).
Тогда амплитуда будет падать как $t^{-1/4} \sim r^{-1/8}$.
Т.е. для $u$ амплитуда будет падать не как $r^{-1/2}$, а как $r^{-5/8}$.
Ну и частота осциляций будет менятся, т.к. для нелинейного осцилятора период колебаний зависит от амплитуды.

peregoudov
Сообщений: 620
Зарегистрирован: 29 дек 2015, 13:17

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение peregoudov » 21 апр 2021, 09:10

На самом деле тут нужна другая замена переменных. Ее идея состоит в рассмотрении того масштабного преобразования, которое оставляет уравнение инвариантным: $r\to\lambda r$, $u\to\lambda^{-2/(n-1)}u$. Если перейти к переменным $t=\ln r$, $v=r^{2/(n-1)}u$, то в новых переменных преобразование принимает вид $t\to t+\ln\lambda$, $v\to v$. То есть в новых переменных уравнение будет инвариантным относительно сдвига $t$, другими словами, автономным. В явном виде

$$ v''+\frac{n-5}{n-1}v'-\frac{2(n-3)}{(n-1)^2}v+v^n=0. $$

Отсюда сразу видно, чем выделен случай n=5: получается нелинейный маятник с потенциалом в форме "жопы Лифшица". Решение Ian'а $\sim r^{-1/2}$ соответствует одному из минимумов, а указанное в википедии $\sim(1+r^2)^{-1/2}$ --- сепаратрисе.
Последний раз редактировалось peregoudov 21 апр 2021, 09:45, всего редактировалось 2 раз.

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение zykov » 21 апр 2021, 09:28

peregoudov писал(а):Source of the post На самом деле тут нужна другая замена переменных.

Какая замена нужна - зависит от целей...

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение zykov » 21 апр 2021, 11:46

zykov писал(а):Source of the post (тут правда в отличии от линейного осцилятора $\overline{(f')^2}$ может быть не пропорционален энергии - нужно проверить).

Есть подозрение, что для такого осцилятора $\overline{(f')^2}$ пропорционально энергии в степени $3/4$. Тогда энергия падает до 0 за конечное время (хотя это время и растёт экспоненциально).
Тогда правда слагаемое $f/t^2$ может помешать занулению.

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение zykov » 21 апр 2021, 14:28

Да, если вернутся к исходному уравнению $u''+2u'/r+u^3=0$, то это тот же осцилятор в поле 4-ой степени с линейным по скорости трением (коэффициент трения падает как "-1" степень со временем).
Форма осциляций само собой негармоническая.
Амплитуда падает как $c - k\ln r$, т.е. падает до 0 за конечное время (экспоненциальное по начальной амплитуде). По мере падения амплитуды период колебаний растёт как "-1" степень амплитуды. Т.е. возможно нуля оно и не достигнет, если период колебаний уйдёт в бесконечность.

Всё же пока не понятно, как оно приближается к нулю на бесконечности при этом загадочном $n=3$...

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение zykov » 21 апр 2021, 14:40

Тут вот в чём может быть дело.
Сам результат $\dot E \sim -E^{3/4}$ основан на усреднении потери энергии на трение по одному периоду.
Для этого сам период должен быть малым, чтобы за один период терялась малая доля энергии.
Но если при малой амплитуде период становится большим, то это уже не верно.
Т.е. нужно анализировать детальнее при малой амплитуде.

peregoudov
Сообщений: 620
Зарегистрирован: 29 дек 2015, 13:17

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение peregoudov » 21 апр 2021, 14:41

Я, конечно, не претендую на точную оценку при $r\to\infty$, но грубую сделать нетрудно. При n=3 уравнение в моих переменных принимает вид

$$ v''-v'+v^3=0. $$

При этом $u=v/r$. Функцией Ляпунова для такого уравнения, очевидно, будет $L=\frac{v^{\prime2}}2+\frac{v^4}4$. Имеем $L'=v^{\prime2}\leq 2L$, а потому $L\leq c_1e^{2t}$. Значит $|v|\leq(4L)^{1/4}\leq c_2e^{t/2}=c_2\sqrt r$. То есть при больших r имеем $|u|\leq c_2r^{-1/2}$.

Гораздо труднее оценить поведение при $r\to0$.

zykov
Сообщений: 1393
Зарегистрирован: 06 янв 2016, 17:41

Сферически симметричные решения нелинейного уравнения Гельмгольца

Сообщение zykov » 21 апр 2021, 16:01

Тоже решил посмотреть в Ваших переменных.
Тут тот же нелинейный осцилятор в яме 4-ой степени. Только теперь "антитрение" линейное по скорости с постоянным коэффициентом трения.
Это "антитрение" будет разгонять осцилятор. Энергия должна расти как $t^4$, значит амплитуда и частота будут расти как $t$.
Значит приближение с усреднением по одному периоду должно работать хорошо.

Тогда получается $u=\frac{\ln r}{r} p(k \cdot \ln^2 r)$, где $p(t)$ - это единичные колебания осцилятора $p''+p^3=0$.
Последний раз редактировалось zykov 21 апр 2021, 16:31, всего редактировалось 1 раз.


Вернуться в «Математика»

Кто сейчас на форуме

Количество пользователей, которые сейчас просматривают этот форум: нет зарегистрированных пользователей и 5 гостей