|
3.2. Пример на построение тестовой задачи.
Рассмотрим одномерное дифференциальное уравнение параболического типа,
для решения которого требуется подобрать устойчивую разностную схему:
| |
\(\displaystyle \frac{\partial u}{\partial t} = 10 \frac{\partial^2 u}{\partial x^2} + x^2, \qquad x \in [0, 1], \quad t \in [0, 1].\) |
(15.13) |
Зададим тестовую функцию в виде:
| |
\(\displaystyle \widetilde{u} = t x^2.\) |
(15.14) |
Определим выражения для производных от тестовой функции (15.14), которые входят в исходное уравнение (15.13):
| |
\(\displaystyle \frac{\partial \widetilde{u}}{\partial t} = x^2, \qquad \frac{\partial \widetilde{u}}{\partial x} = 2 t x, \qquad \frac{\partial^2 \widetilde{u}}{\partial x^2} = 2 t.\) |
(15.15) |
Представим дифференциальное уравнение, истинным решением которого должна быть тестовая функция (15.14), в следующем виде:
| |
\(\displaystyle \frac{\partial \widetilde{u}}{\partial t} = 10 \frac{\partial^2 \widetilde{u}}{\partial x^2} + x^2 + \varphi(t, x).\) |
|
Подставляя выражения (15.15) в данное уравнение, определим вид неизвестной функции \(\displaystyle \varphi(t, x)\), входящей в состав свободного члена:
| |
\(\displaystyle x^2 = 20 t + x^2 + \varphi(t, x) \quad \Longrightarrow \quad \varphi(t, x) = -20 t.\) |
|
Таким образом, мы построили дифференциальное уравнение, которое отличается
от исходного уравнения (15.13) только видом свободного члена и для которого, в то же время, тестовая
функция (15.14) является истинным решением:
| |
\(\displaystyle \frac{\partial \widetilde{u}}{\partial t} = 10 \frac{\partial^2 \widetilde{u}}{\partial x^2} + x^2 - 20 t.\) |
(15.6) |
Однако для того чтобы его решить с использованием численных методов, требуется задать начальное и граничное условия. Зададим эти условия с помощью истинного решения уравнения (15.16), то есть с помощью
тестовой функции (15.14):
| |
\(\displaystyle \widetilde{u}\,(t = 0, x) = 0; \qquad \widetilde{u}\,(t, x = 0) = 0; \qquad \widetilde{u}\,(t, x = 1) = t.\) |
|
Теперь мы можем записать для уравнения (15.16) какую-либо разностную схему, решить уравнение (15.16)
с её помощью, сравнить полученные численные значения со значениями тестовой функции (15.14) в тех же
точках разностной сетки и на основании выражения (15.12) сделать вывод об устойчивости использовавшейся
разностной схемы. Поскольку уравнения (15.13) и (15.16) отличаются только видом свободного члена,
в состав которого в обоих случаях не входит искомая функция, то одинаковые разностные схемы для этих
уравнений будут обладать одинаковым типом устойчивости.
|