Журналы
НЕФТЕГАЗОВЫЕ ТЕХНОЛОГИИ И ЭКОЛОГИЧЕСКАЯ БЕЗОПАСНОСТЬ Серия: Морская техника и технология Серия: РЫБНОЕ ХОЗЯЙСТВО Серия: УПРАВЛЕНИЕ, ВЫЧИСЛИТЕЛЬНАЯ ТЕХНИКА И ИНФОРМАТИКА Серия: ЭКОНОМИКА

Вестник АГТУ (АРХИВНЫЙ)
История Ссылки Новости Контакты
Отправить рукопись
Главная / Журналы / Вестник Астраханского государственного технического университета. Серия: Управление, вычислительная техника и информатика / Номер 4 / ЧИСЛЕННОЕ РЕШЕНИЕ НЕСТАЦИОНАРНОГО ДРОБНОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ В ЗАДАЧАХ МОДЕЛИРОВАНИЯ РАСПРОСТРАНЕНИЯ ТОКСИЧНЫХ ВЕЩЕСТВ В ПОДЗЕМНЫХ ВОДАХ
ЧИСЛЕННОЕ РЕШЕНИЕ НЕСТАЦИОНАРНОГО ДРОБНОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ В ЗАДАЧАХ МОДЕЛИРОВАНИЯ РАСПРОСТРАНЕНИЯ ТОКСИЧНЫХ ВЕЩЕСТВ В ПОДЗЕМНЫХ ВОДАХ
Отправить рукопись Скачать PDF
Текст
Цитировать
Цитирований:

ЧИСЛЕННОЕ РЕШЕНИЕ НЕСТАЦИОНАРНОГО ДРОБНОГО ДИФФЕРЕНЦИАЛЬНОГО УРАВНЕНИЯ В ЗАДАЧАХ МОДЕЛИРОВАНИЯ РАСПРОСТРАНЕНИЯ ТОКСИЧНЫХ ВЕЩЕСТВ В ПОДЗЕМНЫХ ВОДАХ
УДК 504.054
Афанасьева Александра Алексеевна 1
Швецова-Шиловская Татьяна Николаевна 2
Иванов Дмитрий Евгеньевич 3
Назаренко Денис Игоревич 4
Казарезова Елена Викторовна 5
Авторы:
1.  Государственный научно-исследовательский институт органической химии и технологии

2.  Государственный научно-исследовательский институт органической химии и технологии

3.  Государственный научно-исследовательский институт органической химии и технологии

Россия
4.  Государственный научно-исследовательский институт органической химии и технологии

5.  Государственный научно-исследовательский институт органической химии и технологии

Тип:
Статья
DOI:
https://doi.org/10.24143/2072-9502-2019-4-70-80
Страницы:
с 70 по 80
Статус:
Опубликован
Получено:
05.10.2019
Одобрено:
25.10.2019
Опубликовано:
25.10.2019
Классификаторы:
УДК 504.054
Язык материала:
русский
Ключевые слова:
дробная производная, производная по времени, дробное дифференциальное уравнение, разностная схема, подземные воды
Аннотация и ключевые слова
Аннотация (русский):
В настоящее время во многих областях науки для моделирования различных процессов широко применяется теория дробного исчисления. Дифференциальные уравнения с производными дробного порядка используются при моделировании миграции токсичных веществ в пористых неоднородных средах и позволяют более корректно описать поведение токсичных веществ на больших расстояниях от источника. Аналитическое решение дифференциальных уравнений с производными дробного порядка зачастую оказывается очень сложным или даже невозможным. Предложен численный метод решения дифференциальных уравнений в частных производных дробного порядка по времени для описания миграции токсичных веществ в подземных водах. Для численного решения нестационарного дробного дифференциального уравнения разработана неявная разностная схема, являющаяся аналогом известной неявной разностной схемы Кранка - Николсона. Система разностных уравнений представлена в матричном виде. Решение задачи сводится к многократному решению трехдиагональной системы линейных алгебраических уравнений методом прогонки. Представлены результаты оценки распространения токсичного вещества в подземных водах на основе численного метода для модельных примеров. Выполнено сравнение концентраций вещества, полученных на основе аналитического и численного решения нестационарного одномерного дробного дифференциального уравнения. Результаты, полученные с помощью предлагаемого метода и на основании известного аналитического решения дробного дифференциального уравнения, достаточно хорошо согласуются между собой. Относительная ошибка составляет в среднем 9 %. В отличие от известного аналитического решения разработанный численный метод может использоваться при моделировании распространения токсичных веществ в подземных водах с учетом их биодеградации.

Ключевые слова:
дробная производная, производная по времени, дробное дифференциальное уравнение, разностная схема, подземные воды
Текст
Текст произведения (PDF): Читать Скачать
В настоящее время многих исследователей интересует использование математических моделей с дробными производными в силу точности моделей при описании объектов в различных областях науки. Причинами возросшего интереса к дробным производным и дробным дифференциальным уравнениям являются следующие: 1) определение дробной производной нелокально и может быть использовано для математического моделирования сред с памятью; 2) дробные производные и дифференциальные уравнения дробных порядков точнее описывают рассматриваемые явления во многих задачах физики, биологии, геологии; 3) дробными дифференциальными уравнениями можно описывать немарковские процессы. Повышенный интерес к дробному исчислению, наблюдаемый в последнее время, позволяет говорить о нем как о новой области науки [1-10]. Особенный интерес к дробным производным проявляют гидрогеологи в связи с проблемой обоснования безопасности хранения опасных химических и радиоактивных отходов в геологических формациях [1, 2]. В соответствии с имеющимися экспериментальными данными изменение концентрации вещества на больших расстояниях от источника не может быть описано законами классической диффузии, поскольку подчиняется степенному закону убывания, а не экспоненциальному, как это считалось ранее. Данное явление связывают с тем, что область распространения подземных вод представляет собой неоднородную гетерогенную среду, которая включает различные по литологическому составу (по водопроницаемости) породы [1, 2]. Как указано в [1, 2], в таких случаях имеет место так называемая аномальная диффузия, которую подразделяют на субдиффузию и супердиффузию. Субдиффузия протекает намного медленнее классической диффузии, поскольку диффузант захватывается ловушками (дефектами, адсорбционно- или химически активными центрами), уходит в боковые, тупиковые пути и на некоторое время или навсегда выводится из миграционного процесса. Наоборот, супердиффузия осуществляется со скоростями, существенно превышающими классическую диффузию. Этот режим наблюдается, если в системе есть облегченные пути (например, трещины) или присутствуют процессы случайной или направленной адвекции (увлечение диффузанта потоками флюидов). Эти явления невозможно корректно описать классическим уравнением турбулентной диффузии. С учетом этого обстоятельства для более точного описания данных процессов необходимо использовать уравнения диффузии в частных производных дробных порядков [1, 2, 10]. Чтобы учесть в математической модели структуру среды или изменения внешнего поля со временем, нужно заменить обычный оператор производной оператором производной вариационного порядка. Дифференциальные уравнения в частных производных дробных порядков подразделяют на два вида: дифференциальные уравнения с дробной производной по пространству и дробной производной по времени. Чаще всего получить решение для таких уравнений в явном виде не удается. В настоящее время известны аналитические решения дробных дифференциальных уравнений для ограниченного круга задач. Например, в работе [11] приведено аналитическое решение дробного дифференциального уравнения, которое описывает распространение вещества в подземных водах без учета его возможных химических превращений. Поэтому для решения дробных дифференциальных уравнений в частных производных очень эффективными являются приближенные методы [3]. Несмотря на развитый математический аппарат дробного интегрирования и дифференцирования, вопросы численного решения уравнений с дробными производными практически не нашли отражения в существующей литературе [1, 2, 12-15]. В данной работе представлены результаты исследований по разработке численного метода и разностной схемы для решения одномерного дифференциального уравнения с дробной производной по времени. Наибольшее распространение в практике гидрогеологических расчетов получил метод конечных разностей. Однако данный метод обладает некоторыми недостатками [4, 5]. Избавиться от недостатков данного метода можно, используя различные подходы к построению разностной схемы [7]. В данной статье рассматривается уравнение параболического типа, описывающее процесс диффузии, и один из подходов к построению разностных схем для такого типа уравнений - схема Кранка - Николсона, называемая также схемой с весом 1/2. Устойчивость и сходимость решения конечно-разностных уравнений на основе схемы Кранка - Николсона доказана во многих работах [1, 2]. Применение данной схемы позволяет без существенных вычислительных затрат повысить (по сравнению с другими подходами) порядок аппроксимации разностной схемы [13-19]. Таким образом, наиболее целесообразной является разработка неявной разностной схемы для решения дробного дифференциального уравнения на основе известной схемы Кранка - Николсона. Неявная разностная схема для численного решения нестационарного одномерного дробного дифференциального уравнения Рассмотрим общее одномерное уравнение для моделирования миграции токсичных веществ в подземных водах с начальными и граничными условиями вида (1) начальные условия: c(x, 0) = c0, 0 < x < L; граничные условия: где c - концентрация токсичного вещества, мг/м3; c0- начальная концентрация токсичного вещества, мг/м3; D - коэффициент дисперсии в направлении оси х, м2/с; λ - константа скорости превращения токсичного вещества, 1/с; ν - скорость движения воды в порах почвы в направлении оси x, м/с; t - время, с; L - расстояние, на котором концентрация принимает значение с0, м; α - порядок производной. В 1967 г. М. Капуто ввел новое и полезное определение дробной производной, которое сегодня широко известно. Приведем это определение для дробной производной по времени порядка α [20, 21]: (2) где Г (.) - Гамма-функция. Разработаем неявную разностную схему, являющуюся аналогом схемы Кранка - Николсона, для численного решения дробного дифференциального уравнения. Схема Кранка - Николсона может быть записана в виде (3) В качестве сетки возьмем совокупность точек , координаты которых ; ; ; ; ; ; . На данной сетке будем строить разностную аппроксимацию дифференциального уравнения (1), используя шеститочечный шаблон: Ш . Следовательно, схема Кранка - Николсона для модели (1) и шеститочечного шаблона имеет вид (4) (5) (6) Приближенное решение задачи (1) получаем, заменяя значение исходной непрерывной функции на функцию, определяемую в узлах сетки (k - номер узла по времени и l - номер узла по пространству): и означает численную аппроксимацию точного решения (или обозначает численное приближение к точному решению). Аналог неявной разностной схемы Кранка - Николсона для численного решения дробного дифференциального уравнения (1) с учетом выражения для дробной производной Капуто по времени (2) можно сформулировать следующим образом: (7) Выполним последовательно преобразования правой части уравнения (7): (8) Запишем уравнение (1) с учетом (4)-(6), (8): (9) где ; . Проведем упрощающие преобразования. Примем, что коэффициент перед выражением в левой части равен единице, и введем следующие обозначения: (10) После преобразований (8)-(10) система разностных уравнений для первоначальной краевой задачи с дробной производной (1) имеет вид (11) Отметим, что при k = 0 слагаемое с суммой исчезает, что упрощает реализацию метода при проведении вычислений. Выполним преобразование уравнения (11): или (12) Коэффициенты в узлах сетки вычисляются по следующим формулам: Систему разностных уравнений (12) с учетом преобразований можно представить в следующем виде: или Первоначальная краевая задача с дробными производными (2) может быть представлена в матричном виде: для k = 0; (13) для , (14) где ; начальные условия и граничные условия: Матрицы А и В имеют следующий вид: ; 1 Решение задачи (13), (14) в конечном счете сводится к многократному решению системы линейных алгебраических уравнений (СЛАУ) с трехдиагональной матрицей. Для решения таких СЛАУ известны специальные методы, например метод прогонки и метод циклической редукции [9]. Наиболее простым и эффективным методом является метод прогонки, который привлекателен своей экономичностью и простотой реализации в последовательном режиме [8]. Именно поэтому метод прогонки использовался в данном случае для решения системы линейных алгебраических уравнений. Разработанный численный метод решения дробного дифференциального уравнения реализован в программном комплексе [22] для решения задачи распространения токсичных веществ в условиях аномальной диффузии. Алгоритм реализации численного метода представлен на рис. 1. Рис. 1. Алгоритм реализации разработанного численного метода Данный алгоритм предусматривает проверку корректности исходных и расчетных данных на всех этапах вычислений. Новый численный метод реализован на языке программирования Microsoft Visual Basic 6.5, в качестве системы управления базой данных использовался Microsoft Access. Оценка распространения токсичного вещества в подземных водах на основе численного метода В данном разделе приведены результаты численного эксперимента, иллюстрирующие вычислительные свойства нового метода. Выполнено сравнение концентраций вещества, полученных на основе решения нестационарного одномерного дробного дифференциального уравнения разработанным методом и на основании подхода, представленного в работе [11], для двух значений порядка производной по времени - α. Результаты приведены для модельного примера, в котором химические превращения токсичного вещества не рассматривались. Параметры расчетной сетки: h = 0,10; τ = 0,1. Число расчетных точек по осям - 100. На рис. 2 представлены графики зависимости концентрации вещества от расстояния, полученные с помощью разработанного численного метода и на основании аналитического решения, приведенного в работе [11], при α = 0,2. Рис. 2. Зависимость концентрации токсичного вещества от расстояния (α = 0,2) На рис. 3 представлены графики зависимости концентрации вещества от расстояния, полученные с помощью разработанного численного метода и на основании аналитического решения, приведенного в работе [11], при α = 0,6. Рис. 3. Зависимость концентрации токсичного вещества от расстояния (α = 0,6) Из рис. 2 и 3 видно, что графики зависимостей концентрации токсичного вещества от расстояния изменяются симбатно. Результаты, полученные с помощью предлагаемого метода и известного точного решения дробного дифференциального уравнения, достаточно хорошо согласуются между собой. Относительная ошибка составляет в среднем 9 %. В отличие от известного аналитического решения разработанный численный метод может использоваться при моделировании распространения токсичных веществ в подземных водах с учетом их биодеградации. Заключение В результате проведенных исследований разработан новый численный метод решения дифференциального уравнения с производной дробного порядка по времени, а также неявная разностная схема, являющаяся аналогом схемы Кранка - Николсона. Система разностных уравнений представлена в матричном виде. Решение задачи сводится к многократному решению трехдиагональной системы линейных алгебраических уравнений методом прогонки. Представлены результаты оценки распространения токсичных веществ в подземных водах на основе численного метода при решении модельных задач. Полученные результаты позволяют рекомендовать разработанный метод численного решения нестационарного одномерного дробного дифференциального уравнения, основанный на использовании линейных интегро-дифференциальных уравнений достаточно общей структуры, для использования на практике при моделировании процессов тепло- и массопереноса в сильно неоднородных средах.
Список литературы

1. Головизнин В. М., Киселев В. П., Короткин И. А., Юрков Ю. И. Некоторые особенности вычислительных алгоритмов для уравнений дробной диффузии. М.: Изд-во Ин-та проблем безопасного развития атомной энергетики РАН, 2002. 57 с.

2. Головизнин В. М., Киселев В. П., Короткин И. А. Численные методы решения уравнений дробной диффузии в одномерном случае. М.: Изд-во Ин-та проблем безопасного развития атомной энергетики РАН, 2002. 35 с.

3. Меркулова Н. Н., Михайлов М. Д. Методы приближенных вычислений: учеб. пособие. Томск: Изд-во ТГУ, 2011. Ч. III. 184 с.

4. Самарский А. А., Гулин А. В. Численные методы. М.: Наука, 1989. 432 с.

5. Вержбицкий В. М. Численные методы. М.: Высш. шк., 2001. 382 с.

6. Байков В. А., Жибер А. В. Уравнения математической физики. М.; Ижевск: Изд-во Ин-та компьютерных исследований, 2003. 252 с.

7. Самарский А. А., Николаев Е. С. Методы решения сеточных уравнений. М.: Наука, 1987. 130 с.

8. Тихонов А. Н., Самарский А. А. Уравнения математической физики. М.: Наука, 1977. 735 с.

9. Кустикова В. Д. Дифференциальные уравнения в частных производных. Н. Новгород: Изд-во НГУ им. Н. И. Лобачевского, 2011. 54 с.

10. Хенди Ахмед Санд Абделазиз. Численные алгоритмы решения дробных дифференциальных уравнений с запаздыванием: дис. … канд. физ.-мат. наук. Екатеринбург, 2017. 124 с.

11. Shukla H. S., Tamsir M., Srivastava V. K., Kumar J. Approximate Analytical Solution of Time-fractional order Cauchy-Reaction Diffusion equation // Computer Modeling in Engineering & Sciences. 2014. V. 103. N. 1. Р. 1-17.

12. Самко С. Г., Килбас А. А., Маричев О. И. Интегралы и производные дробного порядка и некоторые их приложения. Минск: Наука и техника, 1987. 687 с.

13. Штумпф C. А., Бахтин М. А. Математические методы компьютерных технологий в научных исследованиях: учеб. пособие. СПб.: Изд-во НИУ ИТМО, 2012. 148 с.

14. Fadugba S. E., Nwozo C. R. Crank-Nicolson Finite Difference Method for the Valuation of Options // The Pacific Journal of Science and Technology. 2013. V. 14. N. 2. Р. 136-145.

15. Wani S. S., Thakar S. H. Crank-Nicolson Type Method for Burgers Equation // International Journal of Applied Physics and Mathematics. 2013. V. 3. N. 5. Р. 324-328.

16. Kilbas A. A., Trujillo J. J. Differential equations of fractional order: methods, results and problems // Applicable Analysis. 2001. V. 78. N. 1. Р. 153-192.

17. Atangana1 А., Kilicman А. On the Generalized Mass Transport Equation to the Concept of Variable Fractional Derivative // Mathematical Problems in Engineering. 2014. Р. 1-9.

18. Ерохин С. В. Математическое моделирование напряженно-деформированного состояния вязкоупругих тел с использованием методов дробного исчисления: дис. … канд. техн. наук. М., 2016. 130 с.

19. Бэгли Р. Л., Торвик П. Дж. Дифференциальное исчисление, основанное на производных дробного порядка - новый подход к расчету конструкций с вязкоупругим демпфированием // Аэрокосмическая техника. 1984. Т. 2. № 2. С. 84-91.

20. Caputo M. Linear model of dissipation whose Q is almost frequency independent // Geophis. J. R. Astr. Soc. 1967. V. 13. P. 529-539.

21. Caputo M., Mcinardi F. Linear models of dissipation in an elastic solids // Riv. Nuovo Cimento. 1971. V. 1. Р. 161-196.

22. Свидетельство о государственной регистрации программы для ЭВМ № 2019614483. Программное обеспечение для моделирования распространения загрязняющих веществ в подземных водах / Афанасьева А. А., Назаренко Д. И., Швецова-Шиловская Т. Н., Казарезова Е. В.; зарег. 26.03.2019; опубл. 05.04.2019.

Контент доступен под лицензией Creative Commons Attribution 3.0 Unported

© "Редакция научных журналов"
Поддержка Powered by Editorum,  Инструкции
Войти или Создать
* Забыли пароль?