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

Страницы работы

Содержание работы

ЗАДАНИЕ 1 (МФ, 3 курс, весна)

ЧИСЛЕННОЕ РЕШЕНИЕ КРАЕВОЙ ЗАДАЧИ

ДЛЯ  УРАВНЕНИЯ ТЕПЛОПРОВОДНОСТИ

Плоский слой  (a < x < b)  с начальным распределением температур u(x,0)=u0(x) (a < x <b) в течение некоторого конечного времени (0 < t <= T) подвергается внешнему тепловому воздействию при x=a:

a)  первого рода  u(x,t) =w(t);

b)  либо второго рода   ux(x,t)=-q(x,t);                                                                (1)

c)  либо третьего рода   ux(x,t)=α (u(x,t)- w(t))

и  при x=b:

d)  первого рода  u(x,t) =v(t);

e)  либо второго рода   ux(x,t)=r(x,t);                                                                  (2)

f)  либо третьего рода   ux(x,t)= - β (u(x,t)- v(t)).

В (1)–(2)  функции w(t), v(t) суть температуры среды вблизи границы слоя, функции q(x,t), r(x,t) имеют смысл внешних потоков тепла к границам слоя,  а положительные постоянные α и β  означают коэффициенты теплоотдачи, влияющие на теплосъем с поверхности (например, чем сильнее ветер, тем больше коэффициент теплоотдачи).

В области изменения переменных (a < x < b, 0 < t < T) температура u(x,t) удовлетворяет уравнению теплопроводности

u(x,t) t = u(x,t) xx + f(x,t),                                                                      (3)

где f – плотность внутренних источников тепла.

Уравнение (3) аппроксимируется на равномерной сетке  {x i = a + i h, i=0, …, N; h= (b–a)/N} неявной двухслойной разностной схемой

(u n+1 – u n)/τ = Λ (σ u n+1 + (1– σ) u n) + g n,                                          (4)

где u n = u (t) – решение на временном слое t = n τ, Λ– трехточечный разностный аналог операции двойного дифференцирования по пространственной переменной

(Λz) i = (z i+1 - 2 z i + z i-1) / h2.

Схема (4) имеет различную точность в зависимости от того, какое значение имеет вес схемы σ, как аппроксимируется правая часть и граничные условия.

1.  Схема с погрешностью O(τ + h 2).

σ - любое число из отрезка (0, 1). При σ = 0 схема явная. В противном случае схема неявная. В этом варианте рекомендуется взять σ = 1, т.е. чисто неявную схему.

g n = f(n τ)  с точностью до  величин O(τ).  Первые производные в граничных условиях при этом должны аппроксимироваться с точностью O(τ +h 2).

2.  Схема с погрешностью O(τ 2 + h 2).

σ =0.5.

g n = f(n τ +τ/2)   или g n = (f n+1 + f n)/2 с точностью до  величин O(τ 2).  Первые производные в граничных условиях при этом должны аппроксимироваться также с точностью O(τ 2 + h 2).

3.  Схема с погрешностью O(τ 2 + h 4).

σ =1/2 – h 2/(12 τ), g n = [E + h 2 / 12 Λ]  f(t+τ/2) или g n = f(t +τ/2) + h 2 / 12 fxx(t +τ/2) с точностью до  величин O(τ 2 + h 2) Первые производные в граничных условиях при этом должны аппроксимироваться с точностью O(τ 2 + h 4).

Об аппроксимации граничных условий см. Приложение 1.

Задание состоит в сравнении методов 1, 2, 3 при одном из сочетаний вариантов граничных условий (a, b, c) и (d, e, f) (в индивидуальных заданиях разные сочетания) на тестовых задачах.

Тестовую задачу можно создать стандартным образом, исходя из "придуманного" решения u(x,t). По нему устанавливается правая часть f(x,t) и начальное распределение u0(x), а после фиксации области, т.е. задания параметров a, b, T (при условиях третьего рода следует задать также коэффициенты теплоотдачи), устанавливаются зависимости от времени температур внешней среды и падающих тепловых потоков). Решая затем численно краевую задачу с установленными входными данными, в случае сходимости схемы должны получить приближенные решения, близкие к точным. Не следует брать в качестве решений самые простые функции (например, полиномы невысоких степеней), чтобы погрешность аппроксимации не была равна нулю тождественно.

Похожие материалы

Информация о работе

Тип:
Задания на контрольные работы
Размер файла:
45 Kb
Скачали:
0