Численные методы решения прямых и обратных задач для уравнения диффузии дробного порядка по времени тема диссертации и автореферата по ВАК РФ 05.13.18, кандидат физико-математических наук Иващенко, Дмитрий Сергеевич

• Иващенко, Дмитрий Сергеевич
• кандидат физико-математических наук
• 2008

Оглавление диссертации кандидат физико-математических наук Иващенко, Дмитрий Сергеевич

Введение

1. Аналитические методы решения прямых и обратных задач для уравнения диффузии дробного порядка по времени с постоянным коэффициентом

1.1. Определения и основные свойства интегродифференциальных операторов дробного порядка.

1.2. Уравнение диффузии дробного порядка по времени

1.2.1. Случайное блуждание при непрерывном времени

1.2.2. Вывод уравнения диффузии дробного порядка по времени

1.3. Аналитический подход к решению прямых и обратных задач для уравнения диффузии дробного порядка по времени.

1.3.1. Решение краевой задачи без начальных условий с периодическим источником для уравнения диффузии дробного порядка по времени.

1.3.2. Фундаментальное решение уравнения диффузии дробного порядка по времени.

1.3.3. Решение первой краевой задачи для уравнения диффузии дробного порядка по времени.

1.3.4. Связь между фундаментальным решением уравнения диффузии дробного порядка по времени и функцией Грина первой краевой задачи.

1.3.5. Решение краевой задачи с однородными граничными условиями для уравнения диффузии дробного порядка по времени методом Фурье.

1.3.6. Постановка и решение обратных задач для уравнения диффузии дробного порядка по времени.

1.4. Выводы.

2. Численные методы решения прямых и обратных задач для уравнения диффузии дробного порядка по времени с постоянным коэффициентом 54 2.1. Метод конечных разностей численного решения краевых задач для уравнения диффузии дробного порядка по времени с постоянным коэффициентом.

2.1.1. Разностная схема для уравнения диффузии дробного порядка по времени с постоянным коэффициентом

2.1.2. Обобщенный метод прогонки решения первой краевой задачи для уравнения диффузии дробного порядка по времени с постоянным коэффициентом.

2.1.3. Получение оценок решения разностной краевой задачи

2.1.4. Исследование устойчивости разностной схемы для уравнения диффузии дробного порядка по времени

2.1.5. Определение порядка аппроксимации разностной схемы для уравнения диффузии дробного порядка но времени

2.2. Метод Монте-Карло численного решения краевых задач для уравнения диффузии дробного порядка по времени с постоянным коэффициентом.

2.2.1. Связь дробно-устойчивых распределений с фундаментальным решением уравнения с дробной производной по времени.

2.2.2. Классическая модель случайного блуждания.

2.2.3. Построение дискретной модели случайного блуждания для диффузии дробного порядка по времени

2.2.4. Разложение случайного блуждания на диффузионную и дисперсионную составляющие

2.3. Результаты вычислительных экспериментов.

2.3.1. Численное решение краевой задачи с однородными граничными условиями.

2.3.2. Решение уравнения диффузии дробного порядка по времени методом Монте-Карло.

2.4. Выводы.

3. Численные методы решения прямых и обратных задач для уравнения диффузии дробного порядка по времени с переменным коэффициентом

3.1. Метод конечных разностей численного решения уравнения диффузии дробного порядка по времени с переменным коэффициентом

3.1.1. Интегро-интерполяционный метод построения разностной схемы для уравнения диффузии дробного порядка с переменным коэффициентом.

3.1.2. Обобщенный метод прогонки численного решения краевой задачи для разностного уравнения диффузии дробного порядка по времени с переменным коэффициентом.

3.1.3. Получение оценок решения разностной краевой задачи

3.1.4. Исследование устойчивости разностной схемы для уравнения диффузии дробного порядка по времени с переменным коэффициентом

3.2. Метод минимизации функционала невязки решения обратных задач для уравнения диффузии дробного порядка по времени с переменным коэффициентом.

3.2.1. Постановка прямых и обратных задач.

3.2.2. Свойства решения разностной прямой задачи.

3.2.3. Построение и основные свойства функционала невязки

3.2.4. Метод Левенберга—Марквардта минимизации функционала невязки.

3.2.5. Метод секущих минимизации функционала невязки

3.2.6. Метод Флстчера—Ривса минимизации функционала невязки.

3.3. Результаты вычислительных экспериментов.

3.3.1. Выбор критериев останова алгоритмов в методе минимизации функционала невязки.

3.3.2. Исследование свойств функционала невязки.

3.3.3. Численное решение обратных задач для уравнения диффузии дробного порядка по времени.

3.3.4. Случай постоянного коэфициента

3.4. Выводы.

4. Эволюционные методы решения обратных задач для уравнения диффузии дробного порядка по времени с переменным коэффициентом

4.1. Применение генетических алгоритмов для минимизации функционала невязки.

4.2. Применение алгоритма поиска по шаблону для улучшения приближенного решения

4.3. Результаты вычислительных экспериментов.

4.3.1. Численное решение обратных задач для уравнения диффузии дробного порядка по времени.

4.3.2. Сравнительный анализ эффективности классических и эволюционных алгоритмов

4.3.3. Пример последовательного использования ГА и метода Левенберга—Марквардта.

4.4. Выводы.

Введение диссертации (часть автореферата) на тему «Численные методы решения прямых и обратных задач для уравнения диффузии дробного порядка по времени»

В последнее время значительный интерес представляют физические исследования диффузионных процессов аномальной природы, отклоняющихся от классической гауссовской диффузии, которые встречаются во множестве физических систем (см. [77, 78] и цитированную там литературу). Для аномального диффузионного процесса характерно в первую очередь то, что зависимость среднеквадратического смещения от времени отклоняется от «нормального» линейного закона (x2(t)) = 2 Kit. Здесь К\ и Ка — «обычный» и обобщенный коэффициенты диффузии размерности см2 с-1 и см2с~а соответственно (при а = 1 имеем Г(2) = 1). Показатель аномальной диффузии а ф 1 определяет, будет ли процесс классифицирован как субдиффузионпый (дисперсионный, медленный) при 0 < а < 1 или супердиффузионный (ускоренный, быстрый) при а > 1. Обычно рассматривается область 1 < а < 2, где а = 2 — баллистический предел, описываемый волновым уравнением.

Соотношение (0.1) описывает так называемые «странные» процессы переноса в нелинейных динамических системах [21], то есть негаус-совые процессы, допускающие корреляции на сколь угодно больших пространственно-временных масштабах. Так, например, броуновское движение, для которого функция плотности вероятности имеет вид

1 f х1 \

Рвм(х,Ь) = . - ехр y/MCtf V ±Kyt) ' допускает обобщение с помощью модели случайного блуждания при непрерывном времени (Continuous Time Random Walk — CTRW) на случай субдиффузии или переноса с дисперсией, где

PsD{x,t) ~ С1г°/2Г(1"а)/(2~а) ехр [-c2£2/(2-a> при £ = |z|/ia/2, 0 < a < 1.

Рассматриваемым моделям могут быть поставлены в соответствие дифференциальные уравнения дробного порядка, из которых видно, что данные процессы являются сильно нелокальными и характеризуются широкими корреляциями во времени и/или пространстве, представимыми в виде медленно убывающих по степенному закону ядер в соответствующих интегро-дифференциальных уравнениях.

В работе [78] подробно рассматриваются примеры систем из различных областей (геофизики, геологии, физики, химии, астрономии, экономики), в которых были обнаружены процессы аномальной природы (см. также цитированную там литературу). В частности, с точки зрения аномальной диффузии было предложено объяснение водородного эффекта в структуре силиконовых электродов при некоторых электрохимических условиях, так же, как и в контексте нелинейного электрофореза. Методы дробного исчисления применялись при анализе процессов аномальной диффузии, обнаруженных в аморфных электроактивиых материалах. Аномальная диффузия катионов была обнаружена при изучении механизмов роста образцов наплавленного оксида молибдена, кинетики переноса электронов в пленках из поли-3,4-этилепдиокситиофена, а также атомных процессов переноса и химических реакций в диэлектрических пленках с высокой диэлектрической постоянной. Дробная динамика (fractional dynamics) может лежать в основе статистики объединенной функции плотности вероятности «скорость—координата» движения частицы в турбулентном поле. Обобщение закона Ричардсона с позиций дробного исчисления было предложено для описания переноса воды в ненасыщенном грунте.

Функция плотности вероятности для диффузии дробного порядка описывает самодиффузиониые профили на перколирующих фрактальных структурах. С этих же позиций рассматривается 1//а-шум (т. е. сигнал или процесс со степенной спектральной плотностью, пропорциональной 1 /fa, где / — частота) и соответствующее промежуточное состояние при моделировании молекулярной динамики замерзания воды. Замечено также, что для процессов «старения» в неупорядоченных системах, равно как и в динамических системах, характерна «нелокальность» во времени, подчиняющаяся степенному закону вида (0.1).

Процессы вида (0.1) обнаружены, в частности, в электроактивных материалах [53], лабораторной плазме, турбулентных жидкостях выше некоторого порогового значения числа Рейнольдса, замагпиченных вихревых потоках (см. [21] и указанные там ссылки), в средах с временной дисперсией, таких как биологические ткани и материалы, имеющие самоподобную «архитектуру». Известны многочисленные примеры физических структур, обладающих самоподобием на уровне, промежуточном между микро- и макроуровнями (на мезоуровне), которые стали объектами пристального внимания со стороны физиков-экспериментаторов [23]. Оказалось, что материалы, мезоскопическая структура которых обладает свойством масштабной инвариантности всего 5—8 порядков, имеют уникальные физические свойства, являющиеся результатом их внутренней самоподобпой «архитектуры» [21]. В связи с этим актуальной становится проблема построения адекватных математических моделей таких сред.

С макроскопической точки зрения, диффузионный процесс описывается уравнением диффузии ut(x,t) — Kiuxx(x,t), (0.2) где u(x,t) представляет собой плотность вероятности обнаружить частицу в точке х в момент времени t. С микроскопической точки зрения, диффузия представляет собой марковский процесс, в котором микроскопические частицы выполняют случайные «прыжки» конечной длины и конечной дисперсии. С другой стороны, если рассматривается пемарковский процесс, в котором «прыжки» частиц, выбираются из распределения с длинным временным хвостом ¿~а1, то диффузионный процесс является аномальным (см., например, работу G. Rangarajan и М. Ding [96]). Как показано в работах М. Giona и Е. Roman [57], R. Metzler и J. Klafter [77, 78], W. Schneider [98], W. Wyss [104], R. Gorenflo и F. Mainardi [61], А.И. Саиче-ва и С.Г. Уткина [38, 39], функция плотности вероятности u(x,t), которая описывает движение частиц в случае аномальной диффузии, удовлетвориет уравнению с дробной производной вида

Уравнения вида (0.3) часто встречаются в современной литературе при описании различных аномальных диффузионных процессов [77, 78, 108], в частности процессов, протекающих в сильно пористых (фрактальных) средах [86, 102]. Показатель а может интерпретироваться как характеристика фрактальной размерности «активного» времени [21], в котором реальные блуждания частиц выглядят как случайный процесс; интервал активного времени пропорционален ta.

Важный класс объектов с самоподобной структурой образуют множества, описывающие геометрию перколяции (или протекания), то есть случайного распространения жидкости через среду, причем абстрактные понятия «жидкость» и «среда» могут быть интерпретированы в соответствии с физическим смыслом рассматриваемой задачи. Диффузионные процессы на перколирующих фрактальных структурах существенно не гауссовы и не согласуются с традиционными представлениями о процессах переноса как о случайном броуновском движении частиц в среде. Уравнения вида (0.3) учитывают эффекты памяти и нелокальное™, выходящие далеко за пределы традиционной гауссовой статистики (0.2).

Примеры физических систем, которые описываются в терминах дробного исчисления, приведены в работах Jl. М. Зеленого и А. В. Мило-ванова [21], A. Le Mehaute [69], F. Mainardi [73], Р. Р. Нигматулли-на [87, 88, 89, 90], А. И. Олемского [93], R. Hilfer [62], I. Sokolov [100]. Однако в соответствующей литературе часто отсутствует достоверная информация о результатах проверок адекватности получаемых моделей в ходе реальных экспериментов. В диссертации предлагается подход, использующий предположение, что процесс аномальной диффузии описывается уравнением с дробной производной по времени вида (0.3), однако порядок дифференцирования не известен, то есть ставится задача нахождения вида уравнения аномальной диффузии по реально измеряемым данным. Такие задачи не являются классическими коэффициентными обратными задачами для уравнений в частных производных, вид которых заранее известен [36, 37].

В книгах А.Н. Тихонова и A.A. Самарского [45] и А. Зоммерфелъда [22] рассматривалась известная задача Фурье о распространении температурных волн в почве и было показано, что получить информацию о строении среды можно, измерив ее реакцию на воздействие температурной волны. С другой стороны, в ряде оптических экспериментов по исследованию распространения света в биологических средах было замечено, что источник света с периодически модулированной интенсивностью порождает волну плотности энергии. Эта волна, распространяясь в оптически плотной среде, имеет сферический волновой фронт, форма которого зависит от функции плотности среды. Этому явлению было дано название «диффузионная волна» [74]. Таким образом, по реакции среды на воздействие диффузионной волны молено определить основные параметры среды. Применение данного подхода на мезоуровне позволяет по данным измерений провести проверку адекватности предлагаемой математической модели и определить вид соответствующего дифференциального уравнения.

Актуальность темы диссертации определяется возрастающим интересом к исследованиям в области мезоскопического моделирования процессов переноса. Трудности, связанные с применением аналитических методов решения прямых задач для неоднородных сред па мезоуровне, приводят к необходимости построения эффективных численных методов, чему способствуют и высокие темпы развития компьютерной техники. Разработанные в диссертации методы численного решения обратных задач позволяют па основе реальных данных построить математическую модель аномального диффузионного процесса.

Цель работы заключается в построении точных решений обратных задач для уравнения диффузии дробного порядка по времени с постоянным коэффициентом, а также разработке численных методов решения прямых и обратных задач для уравнения диффузии дробного порядка по времени с переменным коэффициентом.

В рамках указанной цели были поставлены следующие задачи:

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

2. Разработать метод статистического моделирования (Монте-Карло) для численного решения начально-краевых задач для уравнения диффузии дробного порядка по времени с постоянным коэффициентом.

3. На основании классической теории разностных схем построить разностную схему с весами и предложить обобщенный метод прогонки численного решения начально-краевых задач для дифференциального уравнения с дробной производной по времени с постоянным и переменным коэффициентом.

4. Получить условия устойчивости разностных схем и оценки решений разностных краевых задач для дифференциального уравнения с дробной производной по времени с постоянным и переменным коэффициентом.

5. Разработать модифицированные методы минимизации функционала невязки для численного решения обратных задач, состоящих в определении (переменного) обобщенного коэффициента диффузии и показателя аномальной диффузии, и провести их сравнительный анализ.

Методика исследования. При решении поставленных задач использовались методы математической физики, дифференциальных уравнений в частных производных, конечных разностей, статистического моделирования, теории вероятностей, дробного, операционного и вариационного исчислений, а также классические и эволюционные методы безусловной оптимизации.

Научная новизна работы состоит в следующем:

1. Впервые получено решение краевой задачи без начальных условий с периодическим источником методом разделения переменных и построено интегральное преобразование, связывающее решение краевой задачи без начальных условий с периодическим источником для уравнения диффузии дробного порядка по времени с решением аналогичной задачи для уравнения параболического типа.

2. Предложены новые постановки обратных задач для уравнения диффузии дробного порядка по времени, заключающихся в восстановлении обобщенного коэффициента диффузии и дробного показателя диффузии. В некоторых простых случаях решения обратных задач представлены в виде точных аналитических формул, а в общем случае предложен алгоритм численного решения.

3. На основании существующей дискретной модели случайного блуждания для уравнения диффузии дробного порядка по времени разработан метод Монте-Карло численного решения первой краевой задачи. Впервые в явном виде получено разложение случайного блуждания на диффузионную и дисперсионную составляющие.

4. Впервые для уравнения диффузии дробного порядка по времени с переменным коэффициентом с помощью интегро-интерполяционного метода построена разностная схема с весами и разработан обобщенный метод прогонки численного решения первой краевой задачи и получены равномерные оценки решения. Исследована устойчивость разностной схемы с весами для уравнения диффузии дробного порядка по времени с переменным коэффициентом и условие устойчивости получено в явном виде.

5. Впервые предложены постановки обратных задач для уравнения диффузии дробного порядка по времени, заключающихся в восстановлении обобщенного коэффициента диффузии и дробного показателя диффузии, и разработан комплекс программ и реализованы оптимизационные методы ньютоновского типа, сопряженных градиентов, а также эволюционные алгоритмы для их решения.

Теоретическое значение работы заключается в том, что предложенные в ней постановки обратных задач являются новыми, а разработанные для их решения методы представляют собой обобщение существующих подходов па случай уравнений с дробной производной по времени.

Практическая ценность работы состоит в том, что результаты исследований, проведенных в диссертации, позволят на практике по данным измерений строить математические модели сред с неизвестными характеристиками.

Достоверность и обоснованность результатов подтверждена сравнением результатов решения прямых задач теории аномальной диффузии сеточными методами, методом Монте-Карло и аналитическими методами.

На защиту выносятся следующие положения:

1. Разработаны методы аналитического решения обратных задач для уравнения диффузии дробного порядка по времени с постоянным коэффициентом, заключающихся в восстановлении его параметров: коэффициента и порядка временной производной.

2. Разработай алгоритм статистического моделирования (метод Монте-Карло) для решения прямых задач теории аномальной диффузии в однородных средах.

3. Разработан метод конечных разностей численного решения краевых задач для уравнения диффузии дробного порядка по времени как с постоянным, так и с переменным коэффициентом.

4. С помощью классической техники получены условия устойчивости разностных схем с весами для уравнения диффузии дробного порядка по времени с постоянным и переменным коэффициентом, а также оценки решений краевых задач для данного уравнения.

5. Представлены новые постановки обратных задач для уравнения диффузии дробного порядка по времени с переменным коэффициентом, заключающихся в восстановлении его параметров: коэффициента и порядка временной производной. Для их решения разработаны и реализованы в рамках программного комплекса ШРР1/ШВ1 оптимизационные алгоритмы ньютоновского типа, представляющие собой модификации методов Левенберга—Марквардта, секущих и Флетчера—Ривса, а также гибридные эволюционные алгоритмы.

Апробация работы. Основные положения диссертации и отдельные ее результаты докладывались и обсуждались в рамках семинаров, проводимых на кафедре Высшей Математики НГТУ; на семинаре академика

B.Н. Монахова в Институте гидродинамики СО РАН; на семинарах чл.-корр. В.Г. Романова и проф. A.M. Блохина в Институте математики им.

C.J1. Соболева СО РАН, а также на следующих конференциях:

1. Региональная научная конференция «Наука. Техника. Инновации». Новосибирск, 2001, 2002 гг.

2. Всероссийская научная конференция молодых ученых «Наука. Технологии. Инновации». Новосибирск, 2003, 2004, 2005, 2006 гг.

3. Korea—Russia International Symposium on Science and Technology. Novosibirsk, 2002; Tomsk, 2004; Novosibirsk, 2005.

4. Международная конференция «Дифференциальные уравнения, теория функций и приложения», посвященная 100-летию со дня рождения академика Ильи Несторовича Векуа. Новосибирск, 2007 г.

5. Всероссийская конференция но вычислительной математике КВМ-2007. Новосибирск, 2007 г.

6. Международная конференция «Обратные и некорректные задачи математической физики», посвященной 75-летию академика М. М. Лаврентьева. Новосибирск, 2007 г.

Публикации. По теме диссертации опубликовано 15 работ, в том числе 1 работа в журнале из перечня ВАК.

Структура работы. Диссертация состоит из введения, четырех глав, заключения, приложения и списка литературы из 109 наименований. Общий объем диссертации составляет 187 страниц, в том числе основной текст 163 страницы.

Заключение диссертации по теме «Математическое моделирование, численные методы и комплексы программ», Иващенко, Дмитрий Сергеевич

4.4. Выводы

В данной главе разработаны эволюционные алгоритмы решения обратных задач для уравнения диффузии дробного порядка по времени с переменным коэффициентом.

Впервые для решения таких задач предложен гибридный эволюционный алгоритм, сочетающий масштабирование целевого функционала с применением функции селекции и отбор родительских особей по методу рулетки. Для уточнения полученного приближенного решения предложено использовать один из методов прямого поиска — алгоритм поиска по шаблону.

Установлено, что использование генетических алгоритмов позволяет получить решение обратной задачи в тех случаях, когда применение классических оптимизационных методов не дает результата.

Проведен сравнительный анализ эффективности ГА и алгоритмов ньютоновского типа. Показано, что основными недостатками ГА являются мед-лепная сходимость и большие затраты машинного времени.

Показано, как эволюционные алгоритмы могут быть использованы в качестве инструмента для построения начального приближения для классических оптимизационных методов.

Заключение

На основе проведенных в диссертационной работе исследований получены следующие теоретические и практические результаты.

1. Впервые построено решение краевой задачи без начальных условий с периодическим источником для уравнения диффузии дробного порядка по времени с постоянным коэффициентом. Построено интегральное преобразование, связывающее решение краевой задачи без начальных условий с периодическим источником для уравнения диффузии дробного порядка по времени с решением аналогичной задачи для уравнения параболического типа.

2. Построены точные аналитические решения обратных задач для уравнения диффузии дробного порядка по времени с постоянным коэффициентом, заключающихся в восстановлении обобщенного коэффициента диффузии и показателя аномальной диффузии, а также одновременном определении этих параметров.

3. Разработан модифицированный метод Монте-Карло для численного решения прямых задач теории аномальной диффузии, в основу которого положены дискретные модели случайного блуждания для уравнений с дробной производной по времени. Впервые в явном виде получено разложение случайного блуждания на диффузионную и дисперсионную составляющие.

4. Разработан метод конечных разностей численного решения краевых задач для уравнения диффузии дробного порядка по времени с постоянным и переменным коэффициентом. Впервые аналитически получены условия устойчивости соответствующих разностных схем с весами, оценки решений разностных краевых задач и определены порядки аппроксимации на сетке в случаях явной и чисто неявной схем.

5. Разработаны модифицированные оптимизационные алгоритмы ньютоновского типа решения обратных задач для уравнения диффузии дробного порядка по времени с переменным коэффициентом, реализующие метод минимизации функционала невязки. Установлено, что топография поверхности, задаваемой функционалом невязки, имеет овражную структуру и включает области локального минимума. В случае постоянного коэффициента данные обратной задачи моделируются методом Монте-Карло.

6. Разработаны гибридные эволюционные алгоритмы решения обратных задач для уравнения диффузии дробного порядка по времени с переменным коэффициентом, сочетающие масштабирование целевого функционала с селекцией по методу рулетки. Установлено, что в ряде случаев они более эффективны, чем классические методы, а также могут быть использованы в качестве инструмента построения начального приближения для последних.

7. Для численного моделирования процесса аномальной диффузии и решения обратных задач разработан программный комплекс ШРИАИБ!, в рамках которого на основе оптимизационного подхода разработаны численные методы решения обратных задач восстановления обобщенного коэффициента диффузии и показателя аномальной диффузии неоднородной среды, а также задачи одновременного определения этих параметров.

Список литературы диссертационного исследования кандидат физико-математических наук Иващенко, Дмитрий Сергеевич, 2008 год

1. Бахвалов Н. С., Жидков Н. П., Кобельков Г. М. Численные методы. М.: Наука, 1987.

2. Бейтмен Г., Эрдейи А. Высшие трансцендентные функции. — М.: Наука, 1973.

3. Бондаренко А. Н., Иващепко Д. С. Задачи неразрушающего контроля фрактальной среды // Тез. докл. региональной науч. конф. «Наука, Техника, Инновации» — Новосибирск, 2001, 4.1, С. 107—108.

4. Бондаренко А. Н., Иващенко Д. С. Волновые процессы в средах с временной дисперсией // Тез. докл. региональной науч. конф. «Наука, Техника, Инновации» Новосибирск, 2002, 4.1, С. 200-202.

5. Бондаренко А. Н., Иващенко Д. С. Численное решение обратных задач аномальной диффузии // Тез. докл. всероссийской науч. конф. «Наука, Технологии, Инновации» — Новосибирск, 2003, 4.1, С. 223— 224.

6. Бондаренко А. Н., Иващенко Д. С. Численные алгоритмы в обратных задачах восстановления параметров среды на мезоуровне // Тез. докл. всероссийской науч. конф. «Наука, Технологии, Инновации» — Новосибирск, 2004, 4.1, С. 213-214.

7. Бондаренко А. Н., Иващенко Д. С. Метод Монте-Карло в прямых задачах теории аномальной диффузии для неоднородных сред // Тез. докл. всероссийской науч. конф. «Наука, Технологии, Инновации» — Новосибирск, 2005, 4.1, С. 277-278.

8. Бондаренко А. Н., Иващенко Д. С. Численные методы решения обратных задач для уравнения дробной диффузии с гладким коэффициентом // Тез. докл. всероссийской науч. конф. «Наука, Технологии, Инновации» Новосибирск, 2006 4.1, С. 182-183.

9. Бондаренко А. Н., Иващенко Д. С. Численные алгоритмы решения обратных задач аномальной диффузии // Сб. науч. тр. НГТУ. — 2003. №4(34), С. 59-64.

10. Бондаренко А. Н., Иващенко Д. С. Восстановление параметров слоистой среды методом минимизации функционала невязки // Сб. науч. тр. НГТУ. 2004. - №3(37), С. 21-26.

11. Бондаренко А. Н., Иващенко Д. С. Исследование функционала невязки в задачах мониторинга слоистых сред // Сб. науч. тр. НГТУ. — 2004. №4(38), С. 9-14.

12. Бондаренко А. Н., Иващенко Д. С. Оптимизация вычислений в рамках пакета программ «Численное решение обратных задач аномальной диффузии» // Сб. науч. тр. НГТУ. 2004. - №3(37), С. 27-32.

13. Бондаренко А. Н., Иващенко Д. С., Селезнев В. А. Диффузионные волны в средах с остаточной памятью// Науч. вестник НГТУ. — Новосибирск: НГТУ, 2002. №1(12). - С 45-55.

14. Бугров Я. С. Дробные разностные операторы и классы функций // Труды Мат. ин-та АН СССР, 1985, Т. 172, С. 60-70.

15. Васильев Ф. П. Численные методы решения экстремальных задач. — М.: Наука, 1980.

16. Владимиров В. С. Уравнения математической физики. — М.: Наука, 1988.

17. Глушак А. В. О периодических решениях абстрактных дифференциальных уравнений с дробной производной // Вестник ВГУ, Серия физика, математика, 2003, N2 1, С. 96—98.

18. Дэннис Дж., Шнабель Р. Численные методы безусловной минимизации и решения систем нелинейных уравнений. Пер. с англ. — М.: Мир, 1988.

19. Желудев В. А. Производные дробного порядка и численное решение одного класса уравнений в свертках // Диф. уравнения, 1982, Т. 18,11, С. 1950-1960.

20. Зеленый JI. М., Милованов А. В. Фрактальная топология и странная кинетика: от теории перколяции к проблемам космической электродинамики // Успехи физических наук, 2004, Т. 174, № 8, С. 809—852.

21. Зоммерфельд А. Дифференциальные уравнения в частных производных физики. — М.: Изд. ин. лит, 1950.

22. Зосимов В. В., Лямшев Л. М. Фракталы в волновых процессах // Успехи физических наук, 1995, Т. 165, № 4, С. 361-402.

23. Кабанихин С. И. Проекционно-разностные методы определения коэффициентов гиперболических уравнений. — Новосибирск: Наука. Сиб. отд., 1988.

24. Кабанихин С. И., Искаков К. Т. Оптимизационные методы решения коэффициентных обратных задач. — Новосибирск, 2001.

25. Калиткии H.H. Численные методы. — М.: Наука, 1978.

26. Карманов В.Г. Математическое программирование. — М.: Наука, 1975.

27. Кочубей А.Н. Задача Коши для эволюционных уравнений дробного порядка // Дифференциальные уравнения, 1989, Т. 25, № 8, С. 1359— 1368.

28. Кочубей А.Н. Диффузия дробного порядка // Дифференциальные уравнения, 1990, Т. 26, № 4, С. 660-670.

29. Ермаков С.М., Михайлов Г.А. Статистическое моделирование. — М.: Наука, 1982.

30. Моисеев Н. Н., Иванилов Ю. П., Столярова Е. М. Методы оптимизации. — М.: Наука, 1978.

31. Нахушев А. М. Дробное исчисление и его применение. — М.: Физмат-лит, 2003.

32. Рутковская Д., Пилиньский М., Рутковский Л. Нейронные сети, генетические алгоритмы и нечеткие системы. — М.: Горячая линия — Телеком, 2007.

33. Полак Э. Численные методы оптимизации. — М.: Мир, 1974.

34. Рабинович М. И. Введение в теорию колебаний и волн. — М.: Наука, 1984.

35. Романов В. Г. Обратные задачи математической физики. — М.: Наука, 1984.

36. Романов В. Г., Кабанихин С. И. Обратные задачи геоэлектрики. — М.: Наука, 1991.

37. Саичев А. И., Уткин С. Г. Модели дробной диффузии / В сб.: Актуальные проблемы статистической радиофизики. — Нижний Новгород, 2002, Т.1, №1, С. 5-43.

38. Саичев А. И., Уткин С. Г. Асимптотические законы супердиффузии // Журнал технической физики, 2003, Т.73, №7, С. 1—6.

39. Самарский А. А. Теория разностных схем. — М.: Наука, 1977.

40. Самарский А. А., Гулин А. В. Устойчивость разностных схем. — М.: Наука, 1973.

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

42. Стренг Г. Линейная алгебра и ее приложения. — М.: Мир, 1980.

43. Тихонов А. Н., Арсении В. Я. Методы решения некорректных задач.1. М.: Наука, 1979.

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

45. Феллер В. Введение в теорию вероятностей и ее приложения. — М.: Мир, 1984.

46. Уилкинсон Дж. Алгебраическая проблема собственных значений. — М.: Наука, 1970.

47. Учайкин В. В. К теории аномальной диффузии частиц с конечной скоростью свободного движения // Теор. и Мат. Физ., 1998, Т. 115, №1, С. 154—161.

48. Учайкин В. В. Автомодельная аномальная диффузия и устойчивые законы // Успехи физических наук, 2003, Т. 173, №8, С. 847—876.

49. Чукбар К. В. Стохастический перенос и дробные производные // ЖЭТФ, 1995, Т. 108, Вып. 5(11), С. 1875-1884.

50. Arridge S. R., Schweiger М. A General Framework for Iterative Reconstrution Algorithms in Optical Tomography, Using a Finite Element Method // Pre-Print, 1999.

51. Berger J M and Mandelbrot В В 1963 IBM Journal of Research and Development 7 224

52. Bisquert J. Fractional diffusion in the multiple-trapping regime and revision of the equivalence with the continuous-time random walk // Phys. Rev. Lett., 2003, Vol. 91, №1.

53. Carpinteri A and Mainardi F (eds) Fractals and Fractional Calculus in Continuum Mechanics. — Wien and New York, Springer Verlag, 1997.

54. Giona M., Roman E. Fractional diffusion equation on fractals: one-dimensional case and asymptotic behavior // J. Phys. A: Math. Gen.,1992, Vol. 25, P. 2093-2105.

55. Gorenflo R., Mainardi F., Moretti D., Pagnini G., Paradisi P. Discrete random walk models for space-time fractional diffusion // Chemical Physics, 2002, Vol. 284, P. 521-544.

56. Gorenflo R., Mainardi F., Moretti D., Paradisi P. Time fractional diffusion: a discrete random walk approach // Nonlinear Dynamics, 2002, Vol. 29, P. 129-143.

57. Gorenflo R., Vivoli A., Mainardi F. Discrete and continuous random walk models for space-time fractional diffusion // Nonlinear Dynamics, 2004, Vol 38, P. 101-116.

58. Gorenflo R., Mainardi F. Simply and multiply scaled diffusion limits for continuous time random walks // Journal of Physics: Conference series, 2005, Vol. 7, P. 1-16.

59. Hilfer R. Applications of Fractional Calculus in Physics. — Singapore, World Scientific, 2000.

60. Hilfer R., Anton L. // Phys. Rev. E, 1995, Vol.51, P. R848.

61. Hughes В. H. Random Walks and Random Invironments. — Oxford, Clareton Press, 1995.

62. Ismail A. Lecture 18: Non-Markovian Diffusion Equations / In: M.Z. Bazant 18.366 Random Walks and Diffusion. Lecture Notes for Spring 2003. — http://www-math.mit.edu/18.366/

63. Klafter J., Blumen A., Shlesinger М. F. // Phys. Rev. A, 1987, Vol. 35, P. 3081.

64. V. Kolokoltsov, V. Korolev and V. Uchaikin, Fractional stable distributions. J. Math. Sci. 105 (2001), 2569-2576.

65. Le Mchaute A. Fractal Geometries: Theory and Applications. — Boca Raton: CRC Press, 1991.

66. Levenberg K. A method for solution of certain problems in least squares // Quart. Appl. Math., 1944, Vol. 2, P. 164-168.

67. Liu F., Shen S., Anh V., Turner I. Analysis of a discrete non-Markovian random walk approximation for the time fractional diffusion equation // ANZIAM J., 2005, Vol. 46(E), P. C488-C504.

68. Liu F., Zhuang P., Anh V., Turner I. A fractional-order implicit difference approximation for the space-time fractional diffusion equation // ANZIAM J., 2006, Vol. 47(EMAC2005), P. C48-C68.

69. Mainardi F. On the initial-value problem for the fractional diffusion-wave equation // Waves and Stability in Continuous Media, World Scientific, Singapore, 1994, P. 246-251.

70. Mainardi F. Fractional diffusive waves in viscoelastic solids // Nonlinear Waves in Solids, 1995, P. 93-97, ASME/AMR.

71. Mainardi F. The fundamental solutions for the fractional diffusion-wave equation // Appl. Math. Lett., 1996, Vol. 9, No. 6, P. 23-28.

72. Marquardt D. An algorithm for least-squares estimation of nonlinear parameters // SIAM J. Appl. Math., 1963, Vol. 11, P. 431-441.

73. Metzler R., Klafter J. The random walk's guide to anomalous diffusion: a fractional dynamics approach // Phys. Rep., 2000, Vol. 339, P. 1—77.

74. Metzler R., Klafter J. The restaurant at the end of the random walk: recent developments in the description of anomalous transport by fractional dynamics// J. Phys. A: Math. Gen., 2004, Vol. 37, P. R161— R208.

75. Metzler R., Klafter J., Sokolov I. Anomalous transport in external fields: Continuous time random walks and fractional diffusion equations extended // Phys. Rev. E, 1998, Vol. 58. P. 1621-1633.

76. Miller K. S., Ross B. An introduction to the fractional calculus and fractional differential equations. — John Wiley & Sons. Inc., New York, 1993.

77. Montroll E. W., Weiss G. H. // J. Math. Phys., 1965, Vol. 6, P. 167.

78. Montroll E. W., Weiss G. H. // J. Math. Phys., 1969, Vol. 10, P. 753.

79. Montroll E. W., Scher H. // J. Stat. Phys., 1973, Vol. 9, P. 101.

80. Montroll E. W., Scher H. // J. Stat. Phys., 1973, Vol. 34, P. 129.

81. Moré J. J. The Levenberg—Marquardt algorithm: implementation and theory — Numerical Analysis, Lecture Notes in Math., Springer—Verlag, Berlin, 1977, Vol. 60, P. 105-116.

82. Nigmatullin R. R. // Phys. Stat. Sol. (b), 1986, Vol. 133, № 1, P. 425430.

83. Nigmatullin R. R. On the theory of relaxation for systems with Remnant memory // Phys. Stat. Sol. (b), 1984, Vol. 124, P. 389-393.

84. Nigmatullin R. R. The realization of the generalized transfer equation in a medium with fractal geometry // Phys. Stat. Sol. (b), 1984, Vol. 123, P. 534-540.

85. Nigmatullin R. R. To the theoretical explanation of the ■ Universal Response // Phys. Stat. Sol. (b), 1984, Vol. 123, P. 739-745.

86. Nigmatullin R. R., Le Mehaute A. To the nature of irreversibility in linear systems // Magnetic Resonance in Solids, 2004, Vol. 6, No. 1, P. 165—179.

87. Oldham K., Spanier J. Fractional Calculus. — London, New York: Academic Press, 1973.

88. Oldham K., Spanier J. Fractional Calculus and its applications // Bui. Inst. Politehn. I asi. Sec. 1., 1978. Vol. 24, No. 3-4, P. 29-34.

89. Olemskoi A. I. Fractional-differential equations of motion // Pre-Print, 1999.

90. Podlubny I. Fractional Differential Equations. — CA: Academic, San Diego, 1999.

91. Podlubny I. The Laplace Transform method for Linear Differential Equations of the Fractional Order // Pre-Print, 1994.

92. Rangarajan G., Ding M. Anomalous diffusion and the first passage time problem // Phys. Rev. E, 2000, Vol. 62.

93. Scher H, Shlesinger M. F., Bcndlcr J. T. // Phys. Today, 1991, Vol. 44, P. 26.

94. Schneider W. R., Wyss W. Fractional diffusion and wave equations // J.Math.Phys., 1989, Vol. 30, P. 134-144.

95. Shlesinger M. F. // Annual Reviews of Physical Chemistry, 1988, Vol. 39, P. 269.

96. Sokolov I. M., Klafter J., Blumen A. // Phys. Today, 2002, Vol. 55, P. 48.

97. Stanislavsky A. A. Fractional dynamics from the ordinary Langevin equation // Phys. Rev. E, 2003, Vol. 30, P. 87-102.

98. Wegner F., Grossman S. // Zeitschr. Phys. B, 1985, Vol. 59, № 2, P. 197-206.

99. Weiss G. H. Aspects and Applications of the Random Walk. — North-Holland, Amsterdam, 1994.

100. Wyss W. Fractional diffusion equation //J. Math. Phys., 1986, Vol. 27, P. 2782-2785.

101. Wyss M. M., Wyss W. Evolution, its fractional extension and generalization // Pre-Print, 1999.

102. Yuste S, B., Acedo L. An explicit finite difference method and a new von Neumann-type stability analysis for fractional diffusion equations // SIAM J. Numer. Anal., 2005, Vol. 42, P. 1862-1874.

103. Yuste S. B. Weighted average finite difference methods for fractional diffusion equations //J. Comp. Phys., 2006, Vol. 216, P. 264-274.

104. Zaslavsky G. M. Chaos, fractional kinetics, and anomalous transport // Physics Reports, 2002, Vol. 371, P. 461-580.

105. Zauderer E. Partial Differential Equations of Applied Mathematics. — Wiley, New York, 1989.