CC BY
25
Владикавказский математический журнал 2017, Том 19, Выпуск 1, С. 50-58
УДК 519.633
РАЗНОСТНЫЕ СХЕМЫ ДЛЯ УРАВНЕНИЯ ВЛАГОПЕРЕНОСА АЛЛЕРА - ЛЫКОВА С НЕЛОКАЛЬНЫМ УСЛОВИЕМ
М. М. Лафишева, М. А. Керефов, Р. В. Дышекова
Работа посвящена построению разностных схем для уравнения влагопереноса Адлера — Лыкова. Рассмотрена задача с нелокальными граничными условиями типа В. А. Стеклова. Установлен факт сходимости разностной схемы со скоростью 0(Н + т). Проведены численные расчеты с использованием метода окаймления.
Ключевые слова: уравнение влагопереноса, нелокальные условия, разностная схема, априорная оценка, сходимость, метод окаймления.
Вопросы тепло-влагопереноса в почвах являются фундаментальными при решении многих задач гидрологии, агрофизики, гляциологии, экологии, строительной физики и других областей науки. Сложное взаимодействие потоков тепла в почво-грунтах и снежном покрове обусловливает протекание процессов инфильтрации, миграции и морозного пучения, испарения и транспирации, метаморфизма и снеготаяния.
Вопросы теплового и водного режима корнеобитаемого слоя почвы, процессов испарения и транспирации имеют важное значения для сельского хозяйства. Данные процессы определяют условия перезимовки и произрастания сельскохозяйственных культур. Велика роль миграции и инфильтрации влаги в формировании продуктивных запасов влаги на сельскохозяйственных полях.
Исследователи все свое внимание концентрируют на возможности отражения в характере исходных уравнений специфических особенностей изучаемых массивов, их структуры, физических свойств, протекающих в них процессов [2] и т. д. Если уравнение переноса влаги
дШ _ д дЪ дх
где Ш — влажность в долях единицы, х — глубина, Ъ — время, О — коэффициент диф-фузивности, предполагает бесконечную скорость распространения возмущения, то уравнение А. В. Лыкова
дШ Л д2Ш д __1_ ___
дЪ дЪ2 дх
учитывает конечную его скорость.
В то же время весьма существенно введение дополнительного слагаемого даже когда оно мало. Особенно роль последнего становится заметной в процессах, предполагающих быстрые колебания влажности на границах исследуемого образца почвы.
D
dW
дх
(1)
D-
.dW дх
(2)
© 2017 Лафишева М. М., Керефов М. А., Дышекова Р. В.
Правильное истолкование того факта, когда и при каких условиях происходит движение влаги в прямом и обратном направлениях, возможно на основе нового модифицированного уравнения диффузии или уравнения Аллера [8]:
~дГ
д_
дх
'Ш „д2 Ш В—--\-А-
дх
дtдx
В литературе мы находим всевозможные постановки задач для таких уравнений. Одним из таких классов качественно новых задач являются нелокальные задачи для дифференциальных уравнений в частных производных.
Нелокальными называют такие задачи, в которых вместо, или вместе с граничным условием ставятся условия, связывающие значения решения (и, возможно, его производных) во внутренних точках области или в точках границы и каких-либо внутренних точках. Подобные задачи возникают при математическом моделировании процессов различной природы, например, влагопереноса, теплопроводности, при изучении задач математической биологии, задач управления и других.
1. Постановка задачи. Априорная оценка
В замкнутой области <т = {(х,£) : 0 ^ х ^ I, 0 ^ £ ^ Т} рассмотрим задачу
и(0,£) = Аи(1,£), (4)
щ*(0,£) = Апж(М), (5)
и(х, 0) = По(х), (6)
и(х, 0) = П1(х), (7)
где 0 < Со ^ к(х,£) ^ С1, |к;| ^ С2, р, А — положительные постоянные.
Нелокальные условия типа (4)-(5) рассматривались еще В. А. Стекловым [6]. Краевые задачи с нелокальным условием по времени изучались в ряде работ А. И. Кожанова, М. М. Шханукова — Лафишева и др.
В предположении существования достаточно гладкого решения задачи (3)-(7), получим для него априорную оценку. Для чего умножим уравнение (3) скалярно на и;:
(и;, и;) + р(ии, и;) = ((ких)х, и;) + А(ихх;,щ) + щ), (8)
I
где (и, $) = / и&йх, (и, и) = ||и||0-о
Преобразуем, с учетом граничных условий, интегралы, входящие в (8)
( и ; , и ; ) = | и ; | о2 ,
(ии,щ) = ~ |КИо ,
I
а,,, .-,) = ^ -1
о
I I
= их(1, Ь) (к(1, Ь) — Х2к(0, ¿)) — J ки2(1х + ^ ^ кгь2(1х
(/,«*)< 72\\Л\о + 72\Ы\20,
(иххь^г) = ихг(1,Ъ)Щ(1,Ъ) - ихг(0,Ъ)Щ(0,Ъ) - \\ихл\\2 = ихг(1,Ъ)щ(1,Ъ)(1 - X2) - ||и^||2 . Подставляя в (8) полученные соотношения, получим
I
2
+ 21 + 1ки-ах + А Ии^Ио
о
< ихл(1,Ъ) щ(1,Ъ)(1 - X2) + их(1,Ъ) щ(1,Ъ) (к(1, Ъ) - Х2к(0, Ъ))
I
+ ^ / Ьи2(1х + ^ " 1 1 "-■ "2
2„,„о+2Ы1о- (9)
о
Положим к(1,Ъ) = к(0,Ъ), X = 1. Тогда неравенство (9) примет вид
(10)
(1 1 [ (
+ ^ / ку2<1х + 2А ^ Ьи1(1х +
оо
или, с учетом |кг| ^ С2, имеем
I
И^Ио + ^ / ки1Лх + 1М1о < С2 ||пж||2 + \\fWl .
о
Проинтегрируем (10) по т от 0 до t
г г г
2 2 2 2 2 Р |к||0 + Со ||их||0 + 2А ||их*|0 1т ^ С2 / ||их|0 1т+Со ||иоУ0 + /
о о о
Перепишем неравенство (11) в виде
г
11их112 < Иг 111их||2 1т + ^(Ъ),
где ^(Ъ) = И2 ^^ ||/1|0 1т + ||иг(х)|0 + ||и0(х)||Иг, И2 — положительные постоянные.
Воспользуемся леммой Грануолла [1]. Тогда из неравенства (11) получим априорную оценку для решения исходной дифференциальной задачи (3)-(7)
0 (1т+р Ни ||2. (и)
||и*|2 + ||их|2 + / ||их^|0 (1т < И П ||/1|21т + ||их(х)|0 + ||и0(х)||01 . (12) 0 \0
Из оценки (12) следует единственность и устойчивость решения дифференциальной задачи по входным данным.
г
г
2. Разностная схема. Дискретный аналог априорной оценки
Введем в замкнутой области <т сетку
= {х» = ¿Л : I = 0,1,... , N Л = 1/Ж}, ¿т = {¿3 = ¿т : ^ =0,1,..., ^о, т = Т/^'о},
¿5Лт = ¿5Л х ¿5т с шагами Л = т = Т/^0.
А=1
Уо + руй = Л + (1 - - а2)у + ^2у) + Ау о + (13)
; хх ;
у(0,£) = у(М), (14)
Ух,о = Ух,и, (15)
у(х, 0) = ио(х), (16)
У;(х, 0) = и1 (х), (17)
где
У - У У - 2У + у „ ,-+1 з ^ 3-1 к I \ 1
уо = —, ш =-о-' У = У > У = У ' У = У > АУ = УаУ^х> а» = Ь-1/2,
t 2т т2
/ Ч 1 / У»+1 - У» У» - У»-1 \ ,3
= 7 -Г--Ог—Г— , ч> = п
л v л л
Порядок аппроксимации разностной схемы (13)—(17) 0(Л + т).
Перепишем уравнение (13) при а = а1 = а2 и к(х,£) = 1. Примем во внимание также очевидное равенство [4]
ау + (1 - 2а)у + ау = у + ат2уй. Теперь разностное уравнение (13) примет вид
Уо + РУй = Лу + ат2 Луй + АЛу; + (18)
Для получения разностного аналога априорной оценки умножим скалярно уравне-Уо;
(ус,уо) + р((Е - ат2Л) Уй,Уо) = (Лу.уо) + А(Луа,уо) + (<,уо), (19)
где а = а/р, 11и|2 = (и, и) (и, V) = ^и^Н.
Используем в дальнейших преобразованиях следующее равенство [4]
(уй, Уо) =0.5 (||уйЦо);.
Рассмотрим слагаемые равенства (19):
(Е - ат2Л) уг;,уо) = (уй,Уа) - ат2(Луй, у?) = 0.5(||у*||о)* - 5т20.5(|Ух);,
(Луо^о) = (у о,Уо) = -(у 0,У о] + у о Уо - У о Уо
V ; ;/ хх! 1 х ; х ; х х^о
Уо х;
<, Уо ;
2 1 о + 2
Имеет место формула [4]
Подставив полученное в равенство (19) и учитывая выражение (20), будем иметь
+£ (\\vt\lfy + (\ы\ъ\ +1 (ь + ыъ)
о 2 8
|ЫЮ) +
или
р( НУйИо +
Уо х;
о + *
Р°т2 - (Ш\о)г + \ (\\Ух + Ух] 1о)( <
Введем обозначение
2 о.
(20)
(21)
Е3 =р\Ыо+ [^~\)г2Ш\о + \{\\Ух + Ухт
Тогда неравенство (21) примет вид
или
Е^1 _ Е3 < т|Ы|о
Просуммируем последнее выражение по /от 1 до
Е3+1 < Е1 + ^ \\<
'II2
т.
3'=1
При а ^ 1/4 из неравенства (21) получим
(}\\ш\ 1о + С111Iо + \(\\У* + < Р||г/|Но + С1И1о + ^(Ы + ^]1о) + I]
12 , 1 (\\ ,1 , „П 1|2^
I 3'\\2 < \\от-
3'=1
Преобразуем слагаемое правой части последнего неравенства
Ы + У°х]\1 = НиМТ + 2ио,ж]|о < ^1Кх]1о + К,х]1о-
Итак, получена априорная оценка
№||о + НУхй]|§ + ||Ух + Ух] 1о < Щ ||и1 |Ю + ||и1,х]|§ + ||ио,х]|о + I] Н<"Ют
3'=1
Из оценки (22) следует устойчивость и сходимость разностной схемы (13)-(17).
2
о
1
2
2
о
2
;
3. Алгоритм численного решения задачи. Метод окаймления
Пусть о\ = о2 = 1/2, к(х,Ъ) = 1. Тогда разностное уравнение (13) примет вид
А*у- - + ВУ+ = -ъ,
(23)
где
„з-1
К =
2т
" -9. '
А_ 2т
1 О 1 I 1
Уг+1 ~ Щ + Щ-1 з
/I2
а А
А: = Д: = — +
2Л2 2тЛ2'
Запишем граничные условия (14), (15) в виде
1 р а А
У0 - Ум = 2уи - У1 - Ум-1 = 0.
Дополним задачу начальными условиями
У0 = и0,
У = ти1 + и0.
(24)
(25)
(26) (27)
Задачу (23)-(27) будем решать методом окаймления [7, 9], а не методом прогонки, так как матрица системы не трехдиагональна. Представим задачу (23)-(27) в виде
Ам • Хм = ,
(28)
где
/ 1 0 .0 -1 \
А -с В .. .0 0
Ам = 0 А -с . . .0 0
V 0 -1 0 .. . -1 2 /
Матрицу системы (28) Ам перепишем в виде окаймленной матрицы
Ам =
Ам -1 им
им амм
Здесь
/ 1 0 ..... .0
А -с В . . .0
0 А -С .. . 0 , Км
V 0 0 ..... . -с)
Ам-1 =
Обозначим Хм-1 решение усеченной системы
Ам-1 • Хм~1 = Км-
Км -1 /м
1.
(29)
Аналогично Qм-1 есть решение той же системы, но с другой правой частью
Ам-1 • Qм-1 = -им-1- (30)
Решать системы (29) и (30) можно методом прогонки. А затем, зная Хм-ь Qм-1, мы
Хм
Хм = ( , \ + -тХм-г / длг-1 ^ _ (31)
V 0 ) амм + им Qм-1 V 1 /
Итак, разностная задача (23) (27) решается с помощью метода окаймления на каждом временном слое, последовательно, начиная со второго.
Приведем результаты расчетов (в среде МаЙаЬ), при различных входных данных.
Рис. 1. Вычисления проведены со следующими входными данными: u0 = sin x, ui = cosx, A =1 p = 1 f = 0 I = 2n.
Рис. 2. Вычисления проведены со следующими входными данными: uo = sin x, ui = cosx, A =1 p = 1 f = 5 l = 2n.
Рис. 3. Вычисления проведены со следующими входными данными: u0 = sin x, ui = cos x, A = 1 p = 1 f = —1 I = 2n.
Для проверки правильности работы алгоритма сравним точное решение задачи (3) (7) с разностным решением задачи (13) (17) (рис. 4).
Рис. 4.
Возьмем функцию u = sin(x + t) (положив A = 1, р = 1), которая удовлетворяет граничным условиям при l = 2п. Тогда
uo = sin x, ui = cos x, f (x,t) = 2cos(x + t).
При этом max |y — u| = 0.4001.
Литература
1. Ладыженская О. А. Краевые задачи математической физики.—М.: Наука, 1973.—702 с.
2. Нернин С. В., Юзефонич Г. И., Янгарбер В. А. Математические методы прогнозирования водного режима // Материалы объединенной сессии ВАСХНИЛ и АН УзССР.—Ташкент: ФАН, 1967.— С. 279 293.
3. Самарский А. А., Гулин А. В. Устойчивость разностных схем.—М.: Наука, 1973.—415 с.
4. Самарский А. А. Теория разностных схем.—М.: Наука, 1983.—616 с.
5. Соколенке Э. А., Делон В. М., Зелинченко Е. Н., Канокин А. А. Моделирование и управление водно-солевым режимом почв.—Алма-Ата: Наука КазССР, 1976.—180 с.
6. Стеклов В. А. Основные задачи математической физики.—М.: Наука, 1983.—432 с.
7. Фаддеев Д. К., Фаддеева В. И. Вычислительные методы линейной алгебры.—М: Физматгиз, 1960.-656 с.
8. Чудновский А. Ф. Теплофизика почв.—М.: Наука, 1976.—353 с.
9. Шхануков М. X. Локальные и нелокальные краевые задачи для уравнений третьего порядка: Дис.... докт. физ.-мат. наук.—Нальчик, 1995.—225 с.
Статья поступила 2 июня 2016 г.
Лафишева Мадина Мухамедовна
Кабардино-Балкарский государственный университет им. X. М. Бербекова, доцент кафедры ИМОАС РОССИЯ, 360004, г. Нальчик, ул. Чернышевского, 173 E-mail: [email protected]
Керефов Марат Асланбиевич
Кабардино-Балкарский государственный университет им. X. М. Бербекова, доцент кафедры ИМОАС РОССИЯ, 360004, г. Нальчик, ул. Чернышевского, 173 E-mail: keref ovOmail. ru
Дышекова Pameta Владимировна
Кабардино-Балкарский государственный университет
им. X. М. Бербекова, магистр 2-го года обучения кафедры ИМОАС
РОССИЯ, 360004, г. Нальчик, ул. Чернышевского, 173
E-mail: rometa.dikinovaOmail.ru
DIFFERENCE SCHEMES FOR THE ALLER-LYKOV MOISTURE TRANSFER EQUATIONS WITH A NONLOCAL CONDITION
Lafisheva M. M., Kerefov M. A., Dyshekova R. V.
Questions of warm-moisture transfer in the soil are fundamental in solving of various problems of hydrology, agrophysics, ecology and others. Aller-Lykov equation obtained by introducing additional terms in the moisture transfer equation, which take into account the rapid fluctuations of humidity on the boundaries of the test sample of the soil and the final velocity of the perturbation. The paper deals with a boundary value problem for the Aller-Lykov moisture transfer equation with the first type Steklov conditions. A priori estimate for the solution of the differential problem is obtained by the method of energy inequalities, which implies the stability of its solution. Three-level scheme is built. A priori estimate for the solution of the difference problem is obtained. The fact of the convergence of a difference scheme with a rate of O(h + t) is set. The features of the application of the bordering method to the numerical solution of the difference problem are considered. Numerical experiments are conducted, the results of which are attached.
Key words: moisture transfer equation, nonlocal conditions, difference scheme, a priori estimate, convergence, bordering method.
