Сетевое издание
Современные проблемы науки и образования
ISSN 2070-7428
"Перечень" ВАК
ИФ РИНЦ = 1,006

РАЗНОСТНАЯ СХЕМА ДЛЯ УРАВНЕНИЯ ДИФФУЗИИ ДРОБНОГО ПОРЯДКА С СОСРЕДОТОЧЕННОЙ ТЕПЛОЁМКОСТЬЮ

Нахушева Ф.М. 1 Кудаева Ф.Х. 1 Кайгермазов А.А. 1 Кармоков М.М. 1
1 ФГБОУ ВПО «Кабардино-Балкарский государственный университет им. Х.М. Бербекова»
Необходимость изучения краевых задач для дифференциальных уравнений дробного порядка связана с тем, что многие проблемы теории фильтрации жидкости во фрактальной среде, фильтрации жидкости в трещиноватой среде с фрактальной геометрией трещин приводят к дифференциальным уравнениям дробного порядка. Дробные производные применяются при описании физических процессов стохастического переноса, изучении деформационно-прочностных свойств полимерных материалов. Задачи, когда на границе области помещена сосредоточенная теплоемкость некоторой величины. В предыдущей работе построена локально-одномерная схема для многомерного уравнения теплопроводности с сосредоточенной теплоемкостью. Рассмотрен случай многомерной задачи и построена локально-одномерная схема для уравнения диффузии дробного порядка, когда на границах области по каждому направлению помещена сосредоточенная теплоемкость некоторой величины.
априорная оценка
устойчивость разностной схемы
разностная схема
сосредоточенная теплоёмкость
производная дробного порядка
уравнение диффузии
краевая задача
1. Ладыженская О.А. Краевые задачи математической физики. – М.: Наука, – 1973.
2. Нахушева Ф.М., Водахова В.А., Кудаева Ф.Х., Абаева З.В. Локально-одномерная разностная схема для уравнения диффузии дробного порядка с сосредоточенной теплоемкостью // Современные проблемы науки и образования, М.: Издание РАЕ. – 2015. – № 2.
3. Нахушева Ф.М., Шхануков-Лафишев М.Х. Об устойчивости локально-одномерной схемы для уравнения диффузии дробного порядка в многомерной области // Доклады Адыгской (Черкесской) Международной академии наук. – 1997. –Т. 1. – № 2. – С. 23.
4. Самарский А.А. Об одной задаче распространения тепла // Избранные труды А.А. Самарского. – М.: МАКС Пресс, –2003. –С. 1–22.
5. Самарский А.А. Теория разностных схем. – М.: Наука, –1973.
6. Самко С.Г., Килбас А.А., Маричев О.И. Интегралы и производные дробного порядка и некоторые их приложения. – Минск: Наука и техника, –1997.
7. Шхануков-Лафишев М.Х., Нахушева Ф.М., Лафишева М.М., Мамбетова М.М. Локально-одномерная схема для уравнения теплопроводности с сосредоточенной теплоемкостью // Владикавказский математический журнал, Владикавказ. –2013. –Т. 15. – Вып. 4. –С. 59–65.
8. Шхануков-Лафишев М.Х., Нахушева Ф.М. Локально-одномерная схема для уравнения диффузии дробного порядка в р-мерном параллелепипеде // Известия Кабардино-Балкарского научного центра РАН. – 1999. – № 2. – С. 35.
9. Шхануков М.Х. О сходимости разностных схем для дифференциальных уравнений с дробной производной. // Докл. РАН. – 1996. –Т. 348, — С. 746–748.

Необходимость изучения краевых задач для дифференциальных уравнений дробного порядка связана с тем, что многие проблемы теории фильтрации жидкости во фрактальной среде, фильтрации жидкости в трещиноватой среде с фрактальной геометрией трещин приводят к дифференциальным уравнениям дробного порядка. Дробные производные применяются при описании физических процессов стохастического переноса, изучении деформационно-прочностных свойств полимерных материалов. Численным методам решения уравнения диффузии дробного порядка в многомерных областях посвящены работы [3, 8]. Задачи, когда на границе области помещена сосредоточенная теплоемкость некоторой величины, рассмотрены в [4]. В работе [7] построена локально-одномерная схема для многомерного уравнения теплопроводности с сосредоточенной теплоемкостью. В [2] рассмотрен случай многомерной задачи и построена локально-одномерная схема для уравнения диффузии дробного порядка, когда на границах области по каждому направлению помещена сосредоточенная теплоемкость некоторой величины.

1. ПОСТАНОВКА ЗАДАЧИ. УСТОЙЧИВОСТЬ

В области рассматривается краевая задача для обобщенного уравнения диффузии с сосредоточенной теплоемкостью вида:

, (1)

(2)

, , (3)

где коэффициенты удовлетворяют условиям ; ; ; , – регуляризованная производная дробного порядка [6].

Для доказательства устойчивости решения задачи (1)-(3) будем пользоваться методом энергетических неравенств. После умножения уравнения (1) скалярно на получаем энергетическое тождество:

, (4)

где , норма .

Для второго скалярного произведения в правой части (4), интегрируя по частям с учетом краевых условий (2), будем иметь:

. (5)

Используя лемму 1 ([1], стр. 152) и учитывая условия на коэффициенты и , из (5) получим:

.

Для третьего интеграла в энергетическом тождестве (4) с учетом условия запишем:

.

Последнее скалярное произведение в тождестве (4) оценим с помощью -неравенства. Для будем иметь:

.

С учетом полученных неравенств из (4) имеем:

. (6)

Интегрируя (6) по от 0 до , получаем:

.

Пренебрегая положительным интегралом, получаем неравенство:

. (7)

Обозначим в последнем неравенстве и потребуем, чтобы . Будем иметь:

. (8)

Применяя лемму 1 для нестационарных задач ([1], стр. 152), из последнего получим:

, (9)

где – положительная величина, зависящая от коэффициентов уравнения и размеров области . Из априорной оценки (9) следует устойчивость решения задачи (1)-(3) по входным данным задачи, а также его единственность.

2. РАЗНОСТНАЯ СХЕМА

В замкнутой области строится сетка , где – шаг сетки по переменной , – шаг сетки по переменной , , ; – число разбиений по переменной ; – число разбиений по переменной .

Для задачи (1)-(3) построена разностная схема:

(10) , (11)

, (12)

, , (13)

Здесь , , , , , , , , ,

– разностный аналог дробной производной , и в классе достаточно гладких функций справедливо равенство [9].

3. УСТОЙЧИВОСТЬ РАЗНОСТНОЙ СХЕМЫ

Для доказательства устойчивости разностной схемы (10)-(13) используем принцип максимума. Для этого схему приводим к каноническому виду [5]:

(14)

где – связная сетка; – окрестность узла , не содержащая самого узла . Для коэффициентов (14) должны выполняться условия:

(15)

Расписав (10) в индексной форме с учетом в точке получим:

. (16)

Сравнивая (16) с (14), видим:

, , .

Так как по условию , то

, , (17)

Здесь .

В точке , расписывая (11) в индексной форме, будем иметь:

. (18)

Ввиду условий ,, имеем:

, , . (19)

Аналогично, в точке для граничного условия (12) имеем:

. (20)

С учетом условий ,, имеем:

, , . (21)

Из неравенств (17), (19), (21) на основании теоремы 3 [5] для задачи (10)-(13) верна оценка:

, (22)

из которой следует устойчивость решения разностной схемы (10)-(13).

Рецензенты:

Шхануков-Лафишев М. Х., д.ф.-м.н., профессор ФГБУН Институт информатики и проблем регионального управления Кабардино-Балкарского научного центра РАН, г. Нальчик;

Ашабоков Б.А., д.ф.-м.н., профессор Высокогорного Геофизического Института, г. Нальчик.


Библиографическая ссылка

Нахушева Ф.М., Кудаева Ф.Х., Кайгермазов А.А., Кармоков М.М. РАЗНОСТНАЯ СХЕМА ДЛЯ УРАВНЕНИЯ ДИФФУЗИИ ДРОБНОГО ПОРЯДКА С СОСРЕДОТОЧЕННОЙ ТЕПЛОЁМКОСТЬЮ // Современные проблемы науки и образования. – 2015. – № 2-2. ;
URL: https://science-education.ru/ru/article/view?id=22638 (дата обращения: 28.03.2024).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1,674