Уравнение Лыкова учитывает, в отличие от уравнения Аллера, конечную скорость распространения возмущения. Оно выведено из принципов необратимой теплодинамики на основе анализа релаксационных процессов, а не исходя из молекулярных представлений капиллярности. В работе [4] для уравнения влагопереноса Лыкова с краевыми условиями первого рода получены априорные оценки в дифференциальной и разностной трактовках. Получены экономичные факторизованные схемы в многомерной области, устойчивость которых доказана с использованием общей теории устойчивости схем. В работе [6] проведено исследование локально-одномерных схем для уравнения теплопроводности с сосредоточенной теплоемкостью. В работе [3] выведена априорная оценка решения уравнения диффузии дробного порядка с сосредоточенной теплоемкостью, построена схема и получена соответствующая априорная оценка в равномерной метрике. В работе [2] был рассмотрен случай многомерной задачи для уравнения диффузии дробного порядка, когда на границе области по каждому направлению помещена сосредоточенная теплоемкость некоторой величины.
Постановка задачи. Устойчивость
В области рассматривается задача для уравнения влагопереноса с переменными коэффициентами с сосредоточенной теплоемкостью вида:
(1)
(2)
(3)
Для коэффициентов требуем выполнения условий ,
,
,
,
,
,
,
. Здесь
– класс функций, непрерывных вместе со своими частными производными третьего порядка по x и первого порядка по t в
.
Доказательство устойчивости решения задачи (1)–(3) проводим с помощью метода энергетических неравенств [5]. Для чего умножаем уравнение (1) скалярно на и приходим к энергетическому тождеству:
(4)
где скалярное произведение , норма
. Для первого и второго скалярных произведений в (4) имеем:
, (5)
(6)
Для третьего скалярного произведения в (4) применяем интегрирование по частям. Учитывая условия (2), будем иметь:
. (7)
Используя равенство , можем записать
. Тогда для интеграла в (7) будем иметь:
.
Теперь (7) с учётом последнего можем переписать:
. (8)
С учётом условий на k(x, t) и kt(x, t) для интегралов в правой части (8) запишем:
,
. (9)
Для четвёртого интеграла в (4), учитывая условие на q(x, t), можем записать:
. (10)
Для последнего скалярного произведения в тождестве (4) также запишем:
. (11)
С учётом (5)–(11) и условий на , используя лемму ([1], стр.152), из (4) будем иметь
.
Учтём условия ,
и, пренебрегая положительными слагаемыми, отчего неравенство только усилится, последнее неравенство перепишем:
. (12)
где m1, m2, m3, m4 – положительные постоянные величины.
Неравенство (12) проинтегрируем теперь по τ в пределах от 0 до и, учитывая начальные условия (3), имеем интегралы:
,
,
,
,
,
,
.
Учитывая вышеизложенное, перепишем неравенство (12):
или с учётом того, что , перепишем неравенство
, (13)
где M1, M2 – положительные величины. Введем далее обозначение:
.
Тогда неравенство (13), пренебрегая положительными слагаемыми ,
в левой части (отчего неравенство только усилится) перепишем в виде
. (14)
Введем обозначение: , тогда
. С учетом принятых обозначений неравенство (14) можем переписать:
. (15)
Применяя известную лемму [1], из неравенства (15) запишем: . С учётом последнего и вышесделанных обозначений имеем:
.
Из (13), учитывая последнее, получаем неравенство
.
Поскольку , окончательно получим неравенство
. (16)
Здесь положительная величина M(t) зависит от коэффициентов уравнения и размеров области QT. Из априорной оценки (16) заключаем устойчивость решения задачи (1)–(3) по входным данным.
Разностная схема
В строим сетку
, где
, h – шаг сетки, N – число разбиений;
, h – шаг сетки, j0 – число разбиений. Через
обозначается значение сеточной функции y в узле (xi, tj), определенной на сетке
[5]. Дифференциальному уравнению (1) ставится в соответствие параметрическое разностное уравнение вида
, (17)
σ1, σ2 – вещественные параметры, ,
,
,
,
,
,
,
. Разностное уравнение (17) аппроксимирует уравнение (1) с точностью
. Следуя [5], получим аппроксимацию краевых условий (2) с точностью
. Для ku' при x = 0 запишем:
.
Тогда для краевого условия (2) при x = 0 имеем
(18)
Здесь ,
,
. Также для – ku' в (2) при x = l запишем:
.
Тогда аппроксимируется второе граничное условие (2) при x = l в виде
(19)
Здесь ,
.
Первому начальному условию (3) ставим в соответствие разностное условие вида
. (20)
Второе начальное условие в (3) выражением аппроксимируется только с точностью O(τ). Получим аппроксимацию второго порядка по τ. Используя формулу Тейлора, с учётом
запишем:
. (21)
Выразив из уравнения (1) и учитывая,
и перенося
влево и разделив на τ, из (21) получим
,
откуда
. (22)
Итак, получена аппроксимация начального условия с точностью O(τ2).
Таким образом, разностная схема, аппроксимирующая краевую задачу (1)–(3) с точностью имеет вид
(23)
(24)
(25)
(26)
где .
Для реализации на ЭВМ (23)–(26) приводится к системе уравнений с трёхдиагональной матрицей, которая решается методом прогонки на каждом временном слое.
Библиографическая ссылка
Нахушева Ф.М., Куважукова Р.Х., Максидова З.Т. КРАЕВАЯ ЗАДАЧА ДЛЯ УРАВНЕНИЯ ВЛАГОПЕРЕНОСА С СОСРЕДОТОЧЕННОЙ ТЕПЛОЁМКОСТЬЮ // Международный журнал прикладных и фундаментальных исследований. – 2017. – № 7-1. – С. 48-53;URL: https://applied-research.ru/ru/article/view?id=11689 (дата обращения: 16.08.2024).