Пред.
След.
Макеты страниц
Распознанный текст, спецсимволы и формулы могут содержать ошибки, поэтому с корректным вариантом рекомендуем ознакомиться на отсканированных изображениях учебника выше Также, советуем воспользоваться поиском по сайту, мы уверены, что вы сможете найти больше информации по нужной Вам тематике Составные разложения, полученные в п. 4.1.1-4.1.7, являются частным случаем разложений вида где $y$-зависимая переменная, $\varepsilon$-малый параметр, $x$-внешняя переменная, $\zeta$-внутренняя переменная. Составное разложение может рассматриваться как сумма двух членов $F(x ; \varepsilon)=y^{\circ}$ и $G(\zeta ; \varepsilon)=y^{i}-\left(y^{o}\right)^{i}$, т. е. Вместо того чтобы определять внешнее и внутреннее разложения, сращивать их и затем строить составное разложение, Бромберг [1956] и Вишик и Люстерник [1957] предположили, что решение имеет вид (4.2.1б) и пригодно всюду. Следовательно, оно удовлетворяет всем граничным условиям. Взяв внешний предел от (4.2.1б) получим Эта функция должна удовлетворять исходному дифференциальному уравнению, выраженному через внешнюю переменную. Аналогично, функция должна удовлетворять исходному дифференциальному уравнению, записанному через внутреннюю переменную. Чтобы найти приближенное решение, $F$ и $G$ раскладывают по $\varepsilon$ и для каждого уровня приближения получают уравнения и краевые условия. Этот метод применил Чудов [1966] для вязкого обтекания плоской пластины. Вариант метода Бромберга заново открыл О’Малли [1971]. Другой метод составных разложений ранее был предложен Латта [1951]. В соответствии с этим методом предполагалось, что решение также имеет вид (4.2.1б), но $G$ является функцией внешней переменной и внутренней переменной $\zeta$, которая имела более общий вид $g(x) / \delta(\varepsilon)$, а не $x / \delta(\varepsilon)$, причем функция $g$ определялась в результате анализа. Кроме того, Латта исследовал внутреннее разложение и искал специальные функции, которые могут быть использованы для представления $G(x, \zeta ; \varepsilon)$. Ниже мы проиллюстрируем оба метода, применив их к частным примерам. 4.2.1. Уравнение второго порядка с постоянными коэффициентами Как показано в п. 4.1.1, прямое разложение становится непригодным вблизи точки $x=0$, а чтобы описать поведение $y$ в области неравномерности, вводилось внутреннее разложение, использующее преобразование растяжения $\zeta=x \varepsilon^{-1}$. Было показано, что внутреннее разложение содержит функцию $e^{-\zeta}=e^{-x / \varepsilon}$. Поскольку при дифференцировании функция $e^{-x / \varepsilon}$ выражается через саму себя, то нет других специальных функций, необходимых для представления составного разложения. Поэтому Латта пред- положил, что $y$ имеет равномерно пригодное разложение вида Подставляя (4.2.3) в (4.2.2а) и (4.2.2б) и приравнивая к нулю коэффициенты при $\varepsilon^{n}$ и $\varepsilon^{n} e^{-x / \varepsilon}$ для всех $n$, получим уравнения для $f_{n}$ и $h_{n}$. Уравнения для $n=0,1$ и 2 имеют вид Краевые условия имеют вид Здесь мы пренебрегли экспоненциально малыми членами $e^{-1 / \varepsilon} h_{n}(1)$. Подставив (4.2.9) в (4.2.5) и решив полученные уравнения при краевых условиях (4.2.8), получим Подставляя решение первого порядка в (4.2.6) и решая полученные уравнения при краевых условиях (4.2.8), получим Используя полученные решения, найдем разложение (4.2.3). Имеем Легко проверить, что внешнее разложение (предел при $\varepsilon \rightarrow 0$ и фиксированном $x$ ) первых двух членов этого разложения дается выражением (4.1.49), а внутреннее разложение ( $\varepsilon \rightarrow 0$ при фиксированном $\zeta=x / \varepsilon$ ) дается выражением (4.1.76). Таким обра- зом, метод составных разложений дает равномерно пригодное разложение непосредственно без определения внешнего и внутреннего разложений, сращивания их и построения затем составного разложения. Найдем теперь разложение для $y$ с использованием метода Бромберга и Вишика и Люстерника. Предположим, что причем функцией $G(\zeta ; \varepsilon)$ вне внутренней области можно пренебречь (Бромберг [1956]), т. е. $G(\zeta ; \varepsilon) \rightarrow 0$ при $\zeta \rightarrow \infty$, поэтому Поскольку $x=\varepsilon \zeta$, то Так как предполагается, что функцией $G(\zeta ; \varepsilon)$ вне пограничного слоя можно пренебречь, то $F(x ; \varepsilon)$ удовлетворяет граничному условию $y(1)=\beta$. Следовательно, Граничному условию $y(0)=\alpha$ должна удовлетворять функция $F+G$, т. е. Чтобы найти уравнения для $F_{n}$, подставим (4.2.14) в (4.2.2а) и приравняем коэффициенты при равных степенях $\varepsilon$, предполагая $x$ фиксированным. Получим Чтобы определить уравнение для $G_{n}$, выразим сначала (4.2.2а) через внутреннюю переменную $\zeta$. Имеем Подставляя (4.2.15) в (4.2.20), приравнивая коэффициенты при равных степенях $\varepsilon$ и считая $\zeta$ фиксированным, получим Решение уравнения (4.2.18) при условии (4.2.16) имеет вид Следовательно, из (4.2.17) $G_{0}(0)=\alpha-\beta e$, и, таким образом, решение уравнения (4.2.21), стремящееся к нулю при $\zeta \longrightarrow \infty$, имеет вид Решение уравнения (4.2.19) при краевом условии (4.2.16) в случае $n=1$ имеет вид что вместе с условием (4.2.17) дает $G_{1}(0)=-\beta$. Подставляя выражения для $F_{0}$ и $G_{0}$ в (4.2.22), получим Решение этого уравнения, подчиненное условию $G_{1}(0)=-\beta e$ и стремящееся к нулю при $\zeta \rightarrow \infty$, имеет вид Переходя ко второму порядку, получим, что решение уравнения (4.2.19) при условии (4.2.16) имеет вид Поэтому из (4.2.17) получим $G_{2}(0)=-5 \beta e / 2$, и уравнение (4.2.23) примет вид Решение этого уравнения при условии $G_{2}(0)=-5 \beta e / 2$, стремящееся к нулю при $\zeta \rightarrow \infty$, имеет вид Таким образом, первые два члена результирующего равномерно пригодного разложения в точности совпадают с решением (4.1.81), полученным ранее с помощью метода сращивания асимптотических разложений. 4.2.2. Уравнение второго порядка с переменными коэффициентами В качестве второго примера рассмотрим следующую задачу, которая является частным случаем задачи, рассмотренной в II. 4.1.3: Поскольку коэффициент при $y^{\prime}$ положителен, неравномерность будет иметь место в окрестности $x=0$. Чтобы описать поведение $y$ в области неоднородности, необходимо ввести преобразование растяжения $\zeta=x \varepsilon^{-1}$ и внутреннее разложение описывать с помощью функции $e^{-\xi}=e^{-x / \varepsilon}$. Поскольку задача содержит переменные коэффициенты, то $y$ имеет равномерно пригодное разложение вида где функция $g(x)$, которая определитя при анализе, эквивалентна $x$ при $x \rightarrow 0$. Подставляя (4.2.32) в (4.2.30) и (4.2.31) и приравнивая нулю коэффициенты при $\varepsilon^{n}$ и $\varepsilon^{n} e^{-g(x) / \varepsilon}$ при всех $n$, получим уравнения для определения $g, f_{n}$ и $h_{n}$. Первые три уравнения и краевые условия имеют вид Чтобы существовало нетривиальное решение для $h_{0}$, в силу (4.2.33) требуется, чтобы Первый случай (4.2.37) приводит к $g=$ const и должен быть отброішен, так как $g(x) / x \rightarrow 1$ при $x \rightarrow 0$. Следовательно, Решение уравнения (4.2.34), подчиненное условию $f_{0}(1)=\beta$, имеет вид Подставив (4.2.38) в (4.2.35) и решив полученное уравнение при условии (4.2.36), получим Поэтому Рассмотрим далее применение второго варианта метода составных разложений в этой задаче. В этом случае можно применить (4.2.13)-(4.2.17). Подставляя (4.2.14) в (4.2.30), прирав- нивая коэффициенты при равных степенях $\varepsilon$ к нулю и полагая $x$ фиксированным, получим Чтобы получить уравнение для $G_{n}$, выразим (4.2.30) через внутреннюю переменную $\zeta=x / \varepsilon$. Имеем Подставляя (4.2.15) в это уравнение, приравнивая коэффициенты при равных степенях $\varepsilon$ и считая $\zeta$ фиксированным, получим Решение уравнения (4.2.42), подчиненное условию (4.2.16), имеет вид Это выражение вместе с (4.2.17) дает $G_{0}(0)=\alpha-3 \beta$. Следовательно, решение уравнения (4.2.45), стремящееся к нулю при $\zeta \longrightarrow \infty$, имеет вид Подставляя выражение для $F_{0}$ в $(4.2 .43)$ и решая полученное уравнение при условии (4.2.16), придем к Тогда (4.2.17) дает $G_{1}(0)=-16 \beta / 3$, а (4.2.46) принимает вид Решение этого уравнения, стремящееся к нулю при $\zeta \longrightarrow \infty$, имеет вид Поэтому 4.2.3. Краевая задача с начальными условиями для уравнения теплопроводности В качестве третьего примера применения метода составных разложений рассмотрим краевую задачу с начальными условиями для уравнения теплопроводности, поставленную Келлером [1968]. Предположим, что температура $u(x, t ; \varepsilon)$ зависит от одной пространственной переменной $x$, которая изменяется от 0 до $b(\varepsilon t)$, где $b$-известная функция, а $\varepsilon$-малый параметр. Таким образом, $b$-слабо меняющаяся функция $t$. Математически задача записывается в виде Заменив переменную $t$ на $\tau=\varepsilon t$, уравнение (4.2.51) и краевые условия (4.2.52) перепишем в виде Поскольку $\varepsilon$ умножается на $u_{\tau}$, то прямое разложение метода возмущений при малом $\varepsilon$ и фиксированном $\tau$ не может, вообще говоря, удовлетворить начальному условию (4.2.53) и неравномерно вблизи $\tau=0$ всюду, за исключением окрестности концов. Чтобы описать поведение функции $u$ в окрестности $t=0$, необходимо применить преобразование растяжения $t=\tau / \varepsilon$. Как подтвердится ниже, функция, описывающая поведение в этой области, имеет вид $\exp [-g(\tau) / \varepsilon]$, где $g(\tau) / \tau \rightarrow 1$ при $\tau \rightarrow 0$. Поэтому предположим, что равномерно пригодное асимптотическое разложение для $и$ имеет вид Подставляя это разложение в (4.2.53)-(4.2.55) и приравнивая к нулю коэффициенты при $\varepsilon^{n}$ и $\varepsilon^{n} e^{-g(\tau) / \varepsilon}$ для всех $n$, получим и при $n \geqslant 1$ где через $g^{\prime}$ обозначено $d g / d \tau$. Решение задачи (4.2.57) имеет вид Поскольку краевые условия для $h_{0}$ однородны, то уравнение для $h_{0}$ имеет нетривиальное решение, только если $g^{\prime}$ равно одному из собственных значений Соответствующие нормированные собственные функции имеют вид Следовательно, где $a_{0}$ – неизвестная пока функция, которая определится при исследовании уравнения для $h_{1}$. При известном $h_{0}$ уравнение (4.2.61) в случае $n=1$ примет вид Предположим, что $h_{1}$ может быть разложено по собственным функциям $\chi_{k}$, т. е. Подставляя (4.2.68) в (4.2.67) и используя тот факт, что $\chi_{s, x x}=$ $=-g_{s}^{\prime} \chi_{s}$, получим Если мы теперь умножим это уравнение на $\chi_{k}$ и проинтегрируем от $x=0$ до $b(\tau)$, то правая часть обратится в нуль, так как $\chi_{k}$ ортогональна $\chi_{s}$ при $k Это является условием разрешимости задачи (4.2.67). Поскольку Так как $\chi_{k}[b(\tau), \tau]=0$, то Следовательно, (4.2.70) и (4.2.71) приводят к Поэтому решение нулевого порядка имеет вид где $a_{k}$-постоянная, определяемая равенством 4.2.4. Ограничения метода составных разложений При попытке применить метод Латты к нелинейным задачам могут возникнуть осложнения. Кроме того, могут возникнуть трудности, если для описания поведения рассматриваемой функции во внутренней области необходимо использовать большое число специальных функций. Несмотря на эти ограничения, этот метод является огправной точкой для развития метода многих масштабов, описанного в гл. 6. Модифицированный метод составных разложений Бромберга, Вишика и Люстерника преодолевает эти осложнения, как будет показано на примере применения этого метода к нелинейному уравнению описывающему одномерную задачу о космическом корабле ЗемляЛуна, которая изучалась в п. 4.1.7 с помощью метода сращивания асимптотических разложений. Предположим, что составное разложение имеет вид где $\xi=(1-x) / \mu-$ внутренняя переменная, найденная в п. 4.1.7, и $G_{n} \rightarrow 0$ при $\xi \rightarrow \infty$. Начальное условие $t(0)=0$ дает Из (4.2.78) имеем Это разложение, будучи подставленным в (4.2.77), дает Решения этих уравнений, подчиненные условиям (4.2.79), имеют вид Из (4.2.78) и (4.2.82) следует Чтобы определить $G_{0}$ и $G_{1}$, перейдем в уравнении (4.2.77) к внутренней переменной $\xi$. Имеем Подставляя (4.2.83) в это уравнение и приравнивая коэффициенты при равных степенях $\mu$, получим Решение уравнения (4.2.85), стремящееся к нулю при $\xi \rightarrow \infty$, есть $G_{0}=0$, в то время как решение уравнения (4.2.86), стремящееся к нулю при $\xi \rightarrow \infty$, имеет вид Подставляя в (4.2.78) выражения для $F_{0}$ и $F_{1}$ из (4.2.82) и используя найденные значения $G_{0}$ и $G_{1}$, получим разложение, в точности совпадающее с (4.1.202), которое было получено с помощью метода сращивания асимптотических разложений. Упражнения 4.1. Рассмотреть задачу 4.2. Определить разложения второго порядка (трехчленные разложения) для задачи 4.3. Определить равномерно пригодные разложения второго порядка для задач 4.4. Определить равномерно пригодные разложения первого порядка для задачи 4.5. Рассмотреть задачу 4.6. Определить равномерно пригодные разложения первого порядка (двучленные равномерно пригодные разложения) для задачи (2.2.28)-(2.2.30), используя оба варианта MCP. 4.7. Показать, что МСАР не может быть использован для получения равномерно пригодного разложения для 4.8. Рассмотреть задачу, определяемую (2.2.28) при $a \leqslant r \leqslant b$ с граничными условиями Определить двучленные равномерно пригодные разложения, используя а) МСАР и б) МСР. 4.9. Колебания балки с жестко закрепленными концами описываются уравнением Определить разложение первого порядка при малых $\varepsilon$ для $u$ и $\lambda$. 4.10. Теплопередача в одномерном стационарном потоке без диссипации описывается краевой задачей (Xанкс [1971]) 4.11. Определить равномерно пригодное разложение первого порядка для уравнения 4.12. Определить одночленные разложения для решения задач 4.13. Определить одночленные разложения для 4.14. Определить разложение первого порядка для задачи из предыдущего упражнения, если $a(x)$ имеет простой нуль $\mu$, где $0 \leqslant \mu \leqslant 1$. 4.15. Определить равномерно пригодные разложения первого порядка для 4.16. Ламинарный поток в канале с пористыми стенками в случае всасывання приводит к задаче (Праудмен [1960]; Террил и Шреста [1965]) 4.17. Рассмотреть задачу из предыдущего упражнения в случае $a<0$. 4.18. Рассмотреть задачу 4.19. Рассмотреть задачу 4.20. Рассмотреть задачу 4.21. Используя МСР, получить равномерное разложение при больших $R$ для задачи (2.1.37)-(2.1.40), описывающей обтекание сферы.
|
1 |
Оглавление
|