НОУ ИНТУИТ | Численные методы решения уравнений в частных производных. Лекция 9: Применение вариационных принципов для построения разностных схем

Учитесь и получайте официальные документы БЕСПЛАТНО. Вы можете поддержать наш проект.

Регистрация Вход

Твой путь к знаниям!

Опубликован: 25.10.2007 | Уровень: профессионал | Доступ: платный

|

Вам нравится? Нравится 20 студентам

| Поделиться |

Поддержать программу

Проинтегрируем уравнение

$\frac{{\partial}u}{{\partial}t} + {div}{\mathbf{W}} = 0$

по элементарной ячейке разностной сетки. Имеем:

${S_{ij} \frac{{du_{ij}}}{dt} + (W_{\eta ij} \Gamma_{\xi ij} - W_{\xi i + 1j} \Gamma_{\eta i + 1j} - W_{\eta ij + 1} \Gamma_{\xi ij + 1} + W_{\xi ij} \Gamma_{\eta ij} ) = 0, }$

( 9.4)

где $\Gamma$ — длины соответствующих ребер, S_ij — площадь элементарной ячейки.

Так как координаты всех вершин выпуклого четырехугольника известны, то поиск длин, площадей и углов — элементарная геометрическая задача.

Уравнение (9.4}) — дискретный аналог уравнения (9.2). Если возможно определить все потоки в моменты времени tⁿ, t^{n + 1}, а после применить аппроксимацию (9.2) по времени с какими - либо весами, то будет построена разностная схема для расчета температуры.

Учтем, что

$2u \frac{{\partial}u}{{\partial}t} = - 2u{div}{\mathbf{W}}.$

Построим дискретный аналог (9.3):

$\begin{gather*} F_{h} ({\mathbf{W}}^{h} ) = \sum\limits_{i, j \in \omega_{h}} \left[{S_{ij} \cdot \left(\sum\limits_{p, l = 0}^1{\frac{{W_{i + p, j + l}^2}}{{k_{ij}}}}\right) +}\right. \\ {\left. \begin{array}{l} {} \\ {} \\ {} \\ \end{array}2u_{ij}^{({\sigma})} (W_{\eta ij} \Gamma_{\xi ij} - W_{\xi i + 1j} \Gamma_{\eta i + 1j} - W_{\eta ij + 1} \Gamma_{\xi ij + 1} + W_{\xi ij} \Gamma_{\eta ij}) \right].}\end{gather*}$

( 9.5)

Скалярные квадраты, входящие в первое слагаемое дискретного аналога функционала, выражаются через контравариантные проекции следующим образом (рис. 9.7):

$\begin{gather*} W_{ij} = \frac{1}{{{\sin}^2 \varphi_1}}(W_{\eta ij}^2 + W_{\xi ij}^2 + 2W_{\eta ij} W_{\xi ij} \cos \varphi_1 ), \\ W_{i + 1j} = \frac{1}{{{\sin}^2 \varphi_2}}(W_{\eta ij}^2 + W_{\xi i + 1j}^2 - 2W_{\eta ij} W_{\xi i + 1j} \cos \varphi_2 ) \end{gather*}$

и т.д.

Рис. 9.7.

Знак " + " или " - " определяется по правилу: компоненту потока приписывается знак " + ", если проекция потока сонаправлена с внешней нормалью, а знак " - " — если противонаправлена. Таким образом, для углов ячейки $\varphi_1$ и $\varphi_3$ получим знак " + " в последнем слагаемом (проекции одинаковых знаков), а для углов $\varphi_2$ и $\varphi_4$ — знак " - " (проекции теплового потока в произведении разных знаков).

Для получения явной схемы положим в (9.5) вес верхнего слоя по времени $\sigma = 0$ и дифференцируем (9.5) по всем $W_{\xi ij}, W_{\eta ij}$ . Приравнивая производные нулю, получим схему для определения потоков, затем из (9.4) ищем все $u_{ij}^{{n} + 1}$ .

Для построения неявной схемы в (9.4) считаем $\sigma = 1$ , а вместо (9.5) пишем следующую дискретизацию:

$S_{ij} \frac{{u_{ij}^{{n} + 1} - u_{ij}^{n}}}{\tau} = - \sum\limits_{l, p = 0}^1 {(- 1)^{l}(W_\xi )_{i + l, j}^{{n} + 1} \Gamma_{\eta i + l, j} + (- 1)^{p} (W_\eta )_{i, j + p}^{{n} + 1} \Gamma_{\xi i, j + p} .}$

Выражая отсюда неизвестное пока значение $u_{ij}^{{n} + 1}$ в (9.5), получим выражение, зависящее от $u_{ij}^{n}, \{W_\xi ^{{n} + 1} \}, \{W_\eta ^{{n} + 1} \}$ , причем $F_{h} ({\mathbf{W}}^{h} )$ есть сумма квадратов контравариантных проекций.

Дифференцируя, получим линейную систему уравнений для определения потоков. Можно показать, что матрица системы будет обладать следующими свойствами:

${\mathbf{A}} = {\mathbf{A}}^* > 0;$
$\mathbf{A}$ имеет ленточную структуру;
$\mathbf{A}$ является разреженной.

Можно применить эффективные итерационные методы решения системы.

Доказано, что неявная схема будет безусловно устойчивой, а явная — условно устойчивой.

Метод легко обобщается на случай k = k(x, y, u), если уравнение не вырождается. Кроме того, метод может быть обобщен и на случай других граничных условий (не обязательно отсутствия потоков). В этом случае в функционал (9.3) добавляются соответствующие интегралы по границам, а в (9.5) — суммы по поверхностям.

Подробнее об этих схемах можно прочитать в [19.2].

Дальше >>

Численные методы решения уравнений в частных производных

Численные методы решения уравнений в частных производных

Применение вариационных принципов для построения разностных схем

Вопросы и ответы

Студенты

Авторизоваться

Численные методы решения уравнений в частных производных

Численные методы решения уравнений в частных производных

Применение вариационных принципов для построения разностных схем

Вопросы и ответы

Студенты