Теплопроводность в дробном исчислении: приложения к нестационарным методам определения теплофизических характеристик веществ и к задаче Стефана

Шабанова, Муминат Руслановна

Теплопроводность в дробном исчислении: приложения к нестационарным методам определения теплофизических характеристик веществ и к задаче Стефана тема автореферата и диссертации по физике, 01.04.14 ВАК РФ

Шабанова, Муминат Руслановна АВТОР

кандидата технических наук УЧЕНАЯ СТЕПЕНЬ

Махачкала МЕСТО ЗАЩИТЫ

2011 ГОД ЗАЩИТЫ

01.04.14 КОД ВАК РФ

Диссертация по физике на тему «Теплопроводность в дробном исчислении: приложения к нестационарным методам определения теплофизических характеристик веществ и к задаче Стефана»

Автореферат диссертации на тему "Теплопроводность в дробном исчислении: приложения к нестационарным методам определения теплофизических характеристик веществ и к задаче Стефана"

4854886

Шабанова Муминат Руслановна

ТЕПЛОПРОВОДНОСТЬ В ДРОБНОМ ИСЧИСЛЕНИИ: ПРИЛОЖЕНИЯ К НЕСТАЦИОНАРНЫМ МЕТОДАМ ОПРЕДЕЛЕНИЯ ТЕПЛОФИЗИЧЕСКИХ ХАРАКТЕРИСТИК ВЕЩЕСТВ И К ЗАДАЧЕ СТЕФАНА

Специальность 01.04.14 - Теплофизика и теоретическая теплотехника

АВТОРЕФЕРАТ

диссертации на соискание ученой степени кандидата технических наук

2 9 ГРН 2011

Махачкала - 2011

4854886

Работа выполнена в учреждении Российской академии наук «Институт проблем геотермии Дагестанского научного центра»

Научный руководитель: доктор физико-математических наук,

профессор

Назаралиев Магомед-Шафи Ахмедович.

Официальные оппоненты: доктор технических наук, профессор

Вавилов Владимир Платонович

Ведущая организация: ГОУ ВПО «Самарский государственный

технический университет»

Защита состоится 19 октября 2011г. в 15 часов на заседании объединенного диссертационного совета ДМ 002.071.01 при учреждении Российской академии наук «Институт проблем геотермии Дагестанского научного центра РАН» по адресу: 367030, г. Махачкала, пр. И.Шамиля, 39-а, актовый зал

С диссертацией можно ознакомиться в библиотеке УРАН ДНЦ.

Автореферат разослан « , 7 » сентября 2011г.

Ученый секретарь объединенного диссертационного совета

доктор технических наук, профессор Алишаев Мухтар Гусейнович

ДМ 002.071.01 д.т.н.

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

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

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

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

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

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

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

Задачи исследований:

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

- исследовать влияние нелокальностей по времени и координате на распределение температуры при рассмотрении диффузионного и конвективного механизмов переноса тепла;

- разработать метод определения температуропроводности и параметров нелокальностей по времени и координате на основе решения нелокального уравнения теплопроводности и экспериментальных данных нестационарных методов определения распределения температуры;

- разработать математическую модель нелокального переноса тепла на основе дробного исчисления для задачи без начальных условий и его приложения к определению температурных волн в полуограниченных средах;

- на основе математического аппарата дробного исчисления разработать математическую модель задачи Стефана и приложить ее к системе вода - лед.

Объект исследований: процессы переноса тепла в сложных гетеро-фазных средах, в том числе и с фрактальной структурой, с учетом нелокальных свойств по времени и координате.

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

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

Основные научные положения, защищаемые автором.

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

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

3. Метод определения температуропроводности и параметров нелокально-стей по экспериментальным данным нестационарных методов определения распределения температуры по координате и времени и решениям нелокального уравнения теплопроводности для полупрямой.

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

Научная новизна работы.

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

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

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

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

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

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

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

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

Реализация результатов работы. Развитая в работе теория нелокальной теплопроводности в дробном исчислении и его приложения к нестационарным методам измерения теплофизических характеристик веществ могут быть использованы в академических и отраслевых институтах, работающих по соответствующей тематике. Предложенный в работе метод расчета температуропроводности принят ОАО «Геотермнефт-газ» (г. Махачкала) для обработки каротажных данных скважин. По материалам диссертационной работы читается спец. курс «Концепция фрактала и ком-

пьютерное моделирование», выполняются дипломные работы на математическом факультете Дагестанского государственного университета.

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

Апробация работы. Основные положения и выводы диссертации были обсуждены на научных семинарах УРАН «ИПГ ДНЦ» и были предметом обсуждения на следующих научных мероприятиях:

1. XIV международная конференция по химической термодинамике. С.Петербург, 2002.

2. Международная конференция «Возобновляемая энергетика: проблемы и перспективы». Махачкала, 2005.

3. Школа молодых ученых «Актуальные проблемы освоения возобновляемых энергоресурсов». Махачкала, 2006.

4. Школа молодых ученых «Нелокальные краевые задачи и проблемы современного анализа и информатики». Нальчик-Эльбрус, 2007.

5. Российская заочная конференция «Современные наукоемкие технологии». Москва, 2007.

6. Школа молодых ученых «Актуальные проблемы освоения возобновляемых энергоресурсов». Махачкала, 2008.

7. Международный Российско-Абхазский симпозиум. «Уравнения смешанного типа и родственные проблемы анализа и информатики». Нальчик-Эльбрус, 2009.

8. Всероссийская научная конференция с международным участием «Математическое моделирование и краевые задачи». Самара, 2010.

9. Международный Российско-Болгарский симпозиум «Уравнения смешанного типа и родственные проблемы анализа и информатики». Нальчик, 2010.

10. II Международная конференция «Возобновляемая энергетика: Проблемы и перспективы». Махачкала, 2010.

11.1 Всероссийская конференция молодых ученых «Математическое моделирование фрактальных процессов, родственные проблемы анализа и информатики». Кабардино-Балкарская республика, пос. Терскол, 2010. 12. II Международный Российско-Казахский симпозиум «Уравнения смешанного типа и родственные проблемы анализа и информатики». Кабардино-Балкарская республика, Нальчик, 2011.

Публикации. Основное содержание диссертации опубликовано в 23 работах, из которых 6 — статьи в научных рецензируемых журналах из перечня ВАК Структура и объем диссертации. Диссертация состоит из введения, четырех глав, заключения и списка литературы из 129 наименований. Объем работы 120 стр. в том числе 30 рисунков и 2 таблицы.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Во введении обоснована актуальность темы диссертационной работы, сформулированы цели и задачи исследований, основные научные результаты, выносимые на защиту, кратко изложено содержание работы.

Излагается современное состояние уравнений теплопроводности в производных дробного порядка по времени и координате [Mainardi F., Luchko Y, Pagnini G. 2001; Kilbas A.A., Srivastava H.M., Trujillo JJ. 2006; Podlubny I. 1999; Povstenko Y.Z. 2005, Hristov J. 2005] и отмечены нерешенные задачи.

Дается анализ современного состояния нестационарных методов определения температуропроводности [Parker W.J., Jenkins R.J., Butler С.Р., Abbot. 1961; Л.П. Филиппов. 1984; Вавилов В.П. 1991; Фокин В.М., Чернышев В.Н. 2004].

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

Здесь дробная производная по времени (производная Капуто) задается выражением:

а дробные производные по координате (производные Рисса) определены выражениями:

этм=_1_A

а? 2Г(1 - у)соз(7г(1 - у)/ 2) - sf

дрТ^,т)_ 1 а2 -г Т{?,т)

2Г(2 - /3)сов(42 - р)!2) З^2 ¿|£ - ^

Здесь < оо, т > О, Щ, г)- температура, О = Ш0! х2 - безразмерный коэффициент температуропроводности, а = Я!с р- коэффициент температуропроводности, Я - коэффициент теплопроводности, Ср - удельная изобарная теплоемкость, р - плотность вещества, 7 = ух0/¿0 - безразмерная скорость, V-скорость конвективного потока среды,

О < а < 1,0 < Р < 2,0 < у < \,т = = х/х0 - безразмерные время и ко-

ордината; ¡0,Х0- характерные время и масштаб. При рассмотрении задачи Коши для неограниченной прямой получено решение

-ОО -00

оо /*п па

где с (-¡йса) = Т1 С-11" " Функция Миттаг-Леффлера,

аД й Г(аи+1)

Г(£,0)- начальное условие для температуры. Решение (2) в частном случае V —0 совпадает с решением работы [МатагсН Б., ЬисЬко У., Ра^ш в. 2001]. Дается анализ полученного решения. Выяснено влияние учета нелокальности по времени и координате на распределение температуры, как при отсутствии, так и при наличии конвекции. Установлено, что характер распределения температуры в начальные моменты времени при учете нелокальности по времени определяется некоторым характерным временем.

На асимптотическое поведение распределения температуры влияет учет нелокальности по координате. Рассмотрены различные случаи начального распределения температуры. В случае дельта-источника, когда Т(^,0) = Т0с>(^) из (2) получаем

т) = — |йй{ехр(-й4)£аД(--1УК<)ха)+ехр(;^)£а1 (-+ШУ)та)} • (3) 271 0

В частном случае, полагая в (3) а = 1, имеем

ГТ1 со

Г(&т) = -2- \с1к со%(ц - Укгт))ехр(-Вкрт). (4)

71 о

Далее, полагая р=2 и у=1, получим решение

Г(£,г) = ^£к=ехр(-(£-Гг)2/4/)г) ■ (5)

л/4 яОт

Решение (5) при V = 0 совпадает с известным решением [Тихонов А.Н., Самарский A.A. 1972, Фарлоу С. 1985, Лыков A.B. 1967]. Если же в (4) положить у = 1 и ß 1 получим новое решение

_£1_ • (6)

Исследовано влияние учета нелокальности по координате и времени на распределение температуры. Далее, на рисунках ©(£, г) = Т(£, г) / Та - безразмерная температура. Как следует из рисунка 1а влияние нелокальности по времени (пунктирная кривая) и нелокальности по координате (точечная кривая) качественно отличаются. Нелокальность по времени влияет на распределение температуры в начальное время, а нелокальность по пространству влияет на асимптотическое поведение распределения температуры. Исследовано влияние учета нелокальности по времени и координате на распределение температуры с учетом конвективного переноса тепла (рис.16). Учет нелокальности по времени приводит к смещению максимума распределения температуры в область меньших значений координаты.

-10 -5 0 5 ¡0 - 10 -5 0 5 10

Рис.1. Распределение безразмерной температуры в(£,т)по координате для различных значений параметров а, /?; 0=1. а) без конвекции У=0; б) с конвекцией К=/.

Показано, что существует характерное время тх =4/пО, определяющее характер влияния учета нелокальное™ на распределение температуры: Т(°' Т)\Р=г < Т(°> Г)|/м ПРИ 1 < Т(°' > Т(°' Т)\р=л ПРИ т > V ХаРаКТеР распределения температуры по координате, определяемый решениями (5) и (6) зависит от соотношения Т и Тх . На рисунке 2а вблизи начала координат при

уменьшении значения параметра Р величина температуры увеличивается при т<т1и уменьшается при т>хдг- Однако распределение температуры вблизи

начала координат для случаев х < х^и х > хх качественно отличается. В области значений времени х > Хх, как видно на рисунке 26, зависимость распределения температуры вблизи начала координат от учета нелокальности по коор-

динате и времени меняется наоборот: при учете нелокальности по координате величина температуры вблизи начала координат уменьшается.

а) б)

— а-,:,-2 ---« = 1,^ = 1

ötf.r)

0.4

4 ,

Рис.2. Зависимости безразмерной температуры 0(<;,т) от а) т = 0.8<т1, б) т = 1.8 > 1Х Такой характер зависимости распределения температуры вблизи начала координат от соотношения между т и тх объясняется природой сингулярности функций (5) и (6) при \ - 0, когда т -»0. Для решения (5) Т(0, г 0)|^=2 «г",/2; для решения (6) Т(0, г «г4.

В § 2.2 рассмотрена задача для полупространства без учета конвективного члена переноса тепла. Получено общее решение

Щ, г) = - jdk {</£>(£') sin(*#) sin{Ц')Еа1 {-Dk ßr") +

л о о

+ ^)dk)dT'™ßüju(T - г') ■ Н>* V«)

(7)

о о

где ц(т) - заданные начальное и граничное условия. Решение (7) в

частных случаях совпадает с известными традиционными решениями. [Тихонов А.Н., Самарский A.A. 1972, Фарлоу С. 1985, Лыков A.B. 1967]

Рассматривая случай //(г) =Т 0 и а = 1, ß = 2, получим

ехр

4ZX

-ехр

ß+4')2

4£h

1 -er.

24Di.

где erf(x) = -fexp(~z2)dz - интеграл вероятностей. Если же а = \,ß -» 1, 0J

то имеем новое решение:

гй,х) = I _ + (Dt)T - + {Dxf j-

Полагая в (7) = TQ, /j(j) =TQ, получим следующее выражение

Щ,т) = f0 +§(Г0

известным решением

(В)

которое при

т,х) = Тй+(Т0~Т0)егД

а = 1,Р = 2

совпадает

2л/5Т

В остальных случаях (0 < а < 1,1 < Р < 2) имеем новый класс решений. Исследовано влияние учета нелокальности по времени и координате на распределение температуры. В частно-

0.5

- ■ г 1 1 - 1

/■ - /А a-Vfi-2

/ ' 1 ---а=0.1,Д=2

■■■■ д=1,£=1.2 1 1 1 . .. .

сти, на рисунке 3 приведены результаты расчета распределения температуры для полуограниченной прямой. Качественный характер влияния учета нелокальности на зависимость температуры от координаты в случае полуограниченной прямой аналогичен случаю неограниченной прямой: асимптотическое поведение температуры зависит от учета нелокальности по координате.

Рис.3 Зависимость безразмерной температуры 0(£,,т) от .

В §2.3 анализированы экспериментальные данные по определению распределения температуры для органического стекла, полученные в работе [Власов А.Б. 2004]. Метод расчета заключается в определении значений безразмерной температуропроводности и параметров нелокальности аир. Для расчета использовалась формула (8). Подставляя заданные значения Х( и экспериментально определенные значения температуры Г(^(,т()в формулу

(8) определяем значения параметров О, а, р. На рис. 4. приведены значения температуропроводности и параметров нелокальности, рассчитанные по экспериментальным данным. Характерные параметры равны:

=103с, г=///0, х0 =1(Г2м, £=х/х0. Температуропроводность а = £)-а0; а0 = х2а Нй = 10~7м2/с ■

Как следует из рисунка 4, значение температуропроводности у торца образца больше чем в объеме, которое оказывается равным а = 1.12 ■ 10-7 м2/с. При этом видно, что параметры нелокальности аир также меняются. Рис.4. Расчетные значения температуропроводности и параметров нелокальности а, у5 для т =3420с. Сплошная кривая: зависимость О(с); пунктирная -а(ф, точечная -

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

В третьей главе рассмотрена задача без начальных условий. Отмечается, что задача без начальных условий имеет широкое практическое применение. В § 3.1 кратко рассмотрена классическая задача без начальных условий. Причем, в отличие от традиционного изложения, решение классической задачи без начальных условий выведено на основе полученного во второй главе решения (7). Полагая в (7) а=1и р=2 получим,

Т&т) =

2^(7 -Г0) о

ехр

4Д(г-г0)

Л ( -ехр

2 Л

4В(т -Т0);

2£>

ю Г

¡сИс ^(1х'!и{т')к ехр(-М2 (г - г'))

Учитывая, что первый интеграл при Г0 —> -со обращается в нсль, а во вто ром интеграле, используя результат

л/п £

|Ьт(^)ехр(-т2(т - х '))<& = -

ехр(—

4 (Дт4Д(т-т')

(9)

получаем решение в виде

и%'-^—-ехр(--^-)ц(т')

Здесь ц(т) - граничное значение температуры. На практике часто возникают задачи, когда граничное условие задается в виде )Д.(т) = Т0 С08(От) . Представляя интеграл (9) в виде

£ со р2

Ь Г„-3/2 „„„с <э

Г&т) = ¡х-3'2 exp(--i-)cos(Q(T-x))

2-sjTiD $ 4Dx

и используя результаты

® е 2 п Гг\ _

\х'ъп cos(Qr) ехр(- ——)dx = л/я cosßVflI2D) exp(-£Vü/ 2D) , I 4 Dx %

?2 ? Гп _ _

ix'312 sin(ilr) exp(—-—)dx - -Jn —— sin(^Vq/2D)exp(-§Vo/ID), 0 4DX 4

получим известное решение задачи о бегущих температурных волнах (задачи без начальных условий) в виде

Щ,х) = Т0 expKVfi/2D)COS(QT -I^QJID) .

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

В § 3.2 рассматривается задача без начальных условий на основе нелокального по времени уравнения теплопроводности. В этом случае уравнение теплопроводности имеет вид

a'Tfcr) dö2Tfcr)

дта ö|2

граничное условие первого типа Г(О, г) = <р(т). Производная по времени определена соотношением

э»т(£т)= 1 а

дта Т{\-а)дт1{т-2у Получено общее решение в виде

Т(с, X) = ^ )dk • к sin(^) ffi ехр(-гсдт) п ^ , (Ю)

я * л D-k +(т)

со

где <р,,(ш) = |й/техр(-/ют)ф(х) • Рассмотрен случай, когда

—00

ср[т) = Т0 cos(Qr), где Q - частота изменения температуры. Окончательно (10) принимает вид

Щ,х) = Т0exp(-§Vüa/D■ cosfrax/4))-cosfox-tD-sin(rax/4)). (11)

При а = 1 решение (11) переходит в известное решение [Тихонов А.Н., Самарский A.A. 1972]:

Г(£,т) = Г0 ехр(- со^Ох - ^П/2о).

Как видно из (11), учет эффектов памяти приводит к перенормировке характерного масштаба затухания и времени запаздывания температурных волн. Для характерного масштаба затухания температурных волн имеем

/ = Зс0 /л/2 • со^(иа/4), где Зс0 = л/2а2 /ю - характерный масштаб затухания при отсутствии учета эффектов памяти (а=1). Для характерного времени запаздывания температурных волн имеем: { = 70 ■ -У2зт(яа/4), где

7а = ■ а2) - характерное время запаздывания температурных волн без

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

Щ, г) = Та ехр(-£л/П/2/>) • совета / 4).

При а<1 характерная длина затухания увеличивается, а время запаздывания

0.5

10 12 14 16

Щл)

уменьшается.

В § 3.3. рассмотрены некоторые прикладные аспекты задачи без начальных условий. Отметим, что, несмотря на накопленный экспериментальный материал и многочисленность теоретических построений, в изучении тепло-физических свойств сложных пористых многофазных сред, которыми являются большинство осадочных горных пород, недостаточное внимание уделяется такой важной теплофизиче-ской характеристике среды, как ее температуропроводность. На рис. 5 приведены результаты экспериментальных и расчетных данных распределения температуры в верхних слоях Земли.

Рис.5. Экспериментальные и расчетные данные распределения температуры. Предложенный в § 2.3 метод определения температуропроводности позволяет найти ее на основе данных измерений распределения температуры в поверхностных слоях Земли. Экспериментальная кривая (сплошная кривая) взята из работы [Амирханов Х.И., Суетнов В.В., Левкович P.A., Гаирбеков Х.А. 1972] и соответствует результатам режимных наблюдений на геотермической станции (г. Избербаш 10.05.1959 г). Пунктирная и точечная кривые соответствуют расчетным кривым по формуле

Г(£, т) = 70 + Т0 ехр(- ^Üa/D • cos(m / 4))- cos(qt - / D • sin(mx /4)), (12)

1.....1 i i

--——^—

* ■

• Л r

- \\ -

- \ -

- ЗКСП

---a=G5

...... a-QS

i i i 1

В©

где т среднесуточная температура поверхности земли.

0.S

На рисунке 6 приведены результаты расчета температуропроводности как функции координаты по методу, предложенному

Об

в §2.3.

- D(Ç)

......а©

Рис.6 Расчетные значения тем-nepamyponpoeoànocmu(D(Ç)-

0.4

ОД) 0.2

сплошная кривая) и параметра нелокальности по времени

0.5

1.5

(a (Q - точечная кривая) по данным режимных наблюдений

Из рисунка б видно, что значение температуропроводности удовлетворительно согласуется с данными других авторов. Однако, в отличие от традиционного расчета, в данном случае мы рассчитываем зависимость температуропроводности и от координаты. Температуропроводность вычисляется по формуле а= D-xl I/0 ■ В нашем случае параметр х0 = 10м, для характерного параметра

времени, поскольку рассматриваются годичные колебания, t0 равен году. Таким образом a = D-xl//„ = £>-3.17-10~7 м2/с. В соответствии с расчетами получается, что температуропроводность меняется в пределах 0.16-10"7 -0.3-10"7 (м2/с), что соответствует результатам других авторов [Новиков C.B. 2009]. Кроме того, температуропроводность в данном случае увеличивается с глубиной. Как видно на рис, 6. с глубиной меняется параметр нелокальности по времени а. Это соответствует тому, что свойства почвы изменяются вглубь.

Четвертая глава посвящена известной задаче Стефана [Stefan J. 1889]. Постоянный и повышенный интерес к задаче Стефана связан с тем, что она удачно сочетав математические и физические проблемы, охватывает широкий круг фундаментальных проблем математики и физики [Данилюк И.И. 1985, Мейерманов A.M. 1986]. С точки зрения математики она представляет продуктивную модель класса нелинейных задач. С точки зрения физики задача Стефана в ее классической постановке позволяет обобщить себя, охватывая особенности подвижной области фазовых переходов. Межфазная область - это особое состояние вещества, занимающее промежуточное положение между сосуществующими фазами, природа, которого до сих пор не понята до конца, особенно в условиях фазового перехода. В настоящее время как математические, так и физические аспекты задачи Стефана интенсивно развиваются.

В § 4.1 Кратко изложены результаты классической задачи Стефана. В прямой задаче Стефана определяется распределение температуры воды, льда и зависимость межфазной границы от времени 9(т) . Исходная система уравнений имеет вид

эт.&т) аг,&т) = 0> 0<£<19(г); (13)

дт: 8

я! ^ =о, ${т)<^<^ .

от б£,

С краевыми условиями

т;(о,т)=г1, г2(#,о)=/2 ,

(14)

^(£,0 Щх)

при £ = 5(г).

е^ '' " аг

В (13,14) Д(2) - й1(2./0/х2— безразмерный коэффициент температуропроводности, Д1(2) = / ¿,1(2)/;Р1(2)~ коэффициент температуропроводности,

Х0 — безразмерные время и координата, 1(),Х() - характерные время и масштаб, Q - количество тепла, выделяемое или поглощаемое в процессе таяния льда или замерзания воды. Физический смысл условия Стефана

ЯШ, X) _ ^ дТ2 (Е, х) _ ^ д9(т)

= б-

заключается в равенстве потоков теп-

' д\ ' 54 ловой энергии с учетом скрытой теплоты выделяемой или поглощаемой на границе между фазами в зависимости от направления движения фазовой границы. Решения задачи имеют вид

Тх (<?, г) = 4- / 2лД\г )/ег/(а 1для 0 < £ < 3(т),

Ш,т) = /2 •ет/(СТ/2Л/А)/0^-егАа/2^)[ег/({/2^)/ег/(<т/2^5~2)-\\

для .9(г) < £ < 00.

Причем зависимость координаты межфазной границы от времени имеет вид Э(т) = ат1/2, где константа с определяется уравнением, которое получается из условия Стефана

ехр

' £1 I 40,

Л2?2

ехр

/ а2

4£>,

л/А

ет/

1-ег/

йо.-

(15)

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

дха 1 Г(1-ос)та

дха 2 д^ Г(1 - а)та Граничные условия: Г,( 0,*) = *,, Г2(£0) = /2, Г,(#,г) = Г2(#,г) = Г3 при | =

1 2 дЕ,у дх»

где i9(r)- координата подвижной границы раздела фаз. Условие (16) обобщает известное условие Стефана и совпадает с ним, когда у — 1 и Д = 1. В остальных случаях мы имеем нетривиальное обобщение условия Стефана. Решения задачи Стефана в производных дробного порядка имеют вид:

= r, + 2| Sinz ll2P), 0 < £ < т>

1 l-^iA/o") l-Faj(D,/a>)n i z вД

2 "f^ sin(z) z

Зависимость координаты межфазной границы от времени имеет вид

Ä(x) = a(a,ß)-Ta/p , (17)

где уравнение для определения О дается выражением (t3-t,)X1KaP(D1/cß) (t2-Q\2Ka3{D2lo*) ^ rq + g/ß) F^DJ cxß) (1 -Fafi(Dt/^)) r(l-Ya/ß)

Здесь

К* Ä*) = —T~-vf^/cosO? + -/)) + t/,-rfeO)V1(V)!

ЖМуО-у)'«2 V 1 }

TT * 7

где U (x, у) - функция Ломмеля от двух аргументов.

Таким образом, при учете эффектов памяти и пространственных корреляций имеет место обобщенный закон изменения координаты области фазового перехода, определяемый соотношением Э(т) = ст(а, (3) • . Величина <Т оказывается зависящей от двух параметров: а = а(а, р). При а=1 и Р=2 уравнение

(17) совпадает с классическим уравнением. В случае 0 < а < 1 и р=2 уравнение (17) совпадает с уравнением работы [Liu Junyi, Xu Mingyu. 2009].

В §4.3 рассмотрена система вода-лед. Принципиальное отличие закона 9(т) = ст(сс,р)-та'р от ранее известных заключается в том, что 3(г) зависит от двух параметров аир. Для случая системы вода - лед проведены расчеты и анализирована зависимость координаты межфазной границы при различных а и р. Показано, что в области значений параметров а=1 и Р=1 значения координаты межфазной границы аномально растут.

Для определения значения функции Э(т) = а(а, (?) • та/р установим, как зависит а(а, Р) от параметров а, Р для случая системы вода-лед. Исходя из значений параметров температура льда =-5°С, воды /2 = 0°С, ti =0°С, Р»л, = 917.4кг/м3,?^ =2.24Вт/м- К, = ЪЪОкДжЫт, альда = 1.2-10^лг2/с, уравнение для определения <т(а,р) принимает вид ^р(1.2/0р)^27озГ(1 + а/р)с,п Fafi( 1,2/ор) ' Г(1-усс/р) Численные значения ст(а, Р) приведены в виде таблицы для нескольких значений а, р. В таблице по строке меняются значения а е (0,1), по столбцу меняются значения р е (1,2) с шагом 0.05. Исходя из этих значений и определяя а(а, Р) для фиксированных а е (0,1) и р е (1,2), можно определить зависимость координаты от а и р. На рисунке 7а приведены результаты расчета зависимости Э(т) от значений а для некоторых фиксированных значений р. Для значения р=1.8 эта зависимость имеет монотонно возрастающий характер при уменьшении значения параметра а. Для значения р=1.05 эта зависимость становится немонотонной и имеет пороговый характер. На рисунке 76 приведены результаты расчета зависимости Э(т) от значений р для некоторых фиксированных значений а. Как следует из рисунка, характер зависимости -Э(т) от значений р определяется значением параметра а. Для значения а=0.7 эта зависимость имеет монотонно убывающий характер при уменьшении значения параметра р. Для значений а=1 эта зависимость становится немонотонной и имеет пороговый характер. Общая зависимость Э(т) = а(а, Р) • Та/р от параметров а, р приводится на рисунке 8. Как видно, вблизи области значений а=1,р=1 значения функции Э(т) = ст(а,Р) • та/р аномально растут.

Рис. 7. Зависимости безразмерной координаты межфазной границы: а) от параметра а; б) от параметра р.

Таблица значений с(а, Р)

' 2272 2277 Ш 2 286 229 2294 2297 23 2303 2306 2309 2311 2314 2316 2318 232 2322 2324 2326 2327

1.ЯЯ 2Ш 21)3 1Ш 2.056 2068 2079 2С89 2.098 2.107 2.116 2.124 2.131 2.133 2.145 2151 2157 2.163 2.169 2.174

1.719 1.744 1.768 1.789 Ш 1.828 1Мб Ш т 1.893 1.И7 Ш1 13® 1943 1536 1567 1577 1987 1997 3.003

1.45! 1.457 ш ш ш Ш Ш 1642 1.663 1ЙЗ 1.702 1.719 1.736 1.732 1.768 1.782 1.796 ¡Л 1.823 1.835

1217 1254 1229 1321 1352 1381 1.409 1.435 1.45 1.484 1306 1328 1348 1368 1386 1ЯИ 1.621 1.638 163 1.669

1.009 1Л8 1.083 1.121 1.154 1.Ш 1217 Ш 1273 13 1325 13« 1372 1394 1.416 1.436 1.436 1.474 1.493 131

0.8303 0.8699 05079 05442 09789 1Ш2 1.044 1175 1.104 1.132 1.139 1.185 121 1234 1257 128 1301 1322 1342 1361

0.6786 0.7174 от от «5257 0.8593 0.8917 0.9229 оет 0582 1.01 11)37 1МЗ 1.088 1.113 1.136 1.159 1181 1202 1222

0.5515 0.5834 0.6243 от от 0.7261 0.758 0.789 01191 0.8482 08764 от от 03>36 от 1.С03 11123 Ш1 1Л73 1.095

0.4459 0.4Ш2 03139 ш 03793 0.611 0.6419 0.672 0.7014 0.73 0.7379 0.7831 0.8115 0.8372 0.8622 ОШ 09103 05333 05358 09776

03589 0302 0.4213 0.4521 04224 03122 03416 03704 03987 0.6263 0.6534 067» 0.7038 0.7311 0.7538 0.78 еш от 0.849 0.871

от 03159 03441 03723 0.4303 0.428 0.4533 0.4827 от от 03617 03872 0.6123 0.6368 ОМ» 0 6845 0.7077 0.7303 0.7525 0.7743

от 02547 0.28 03053 0331 03566 0382 0.4012 0.4323 0.4571 0.4817 03039 03298 03534 0.5766 03994 0В19 0.644 0.6636 от

от 02048 02272 025 0273 02962 03194 03427 05139 0389 0.412 0.4348 0.4575 0.4798 0302 032В 03454 03667 03877 0 6083

от 0.1641 0.1838 0204 02245 02454 02663 02877 0309 03304 03317 0373 03942 0.4132 0.4362 0.4569 0.4775 0.4978 03179 03378

0.1146 0.1311 0.1482 0.166 0.1842 02028 02217 0241 от 02799 от от 0339 03387 03783 03978 0.4172 0.4363 0.4537 0.4747

0Ш7 0.1045 0.1193 0.1347 0.1307 0.1672 0.1841 02014 02189 02367 02547 0 2728 0291 03092 03273 0 3458 0364 03822 0.4003 0.4184

0.07106 0.08302 009573 0.1091 0.1231 0.1376 0.1526 0.168 0.1838 0.1998 02161 02327 Ш4 02662 02831 03801 03171 03342 03512 03682

03057 от 03832 0.08816 от 0.113 0.1» 0.1399 0,1339 0.1684 0.1831 0.1981 02134 02288 02444 02601 02759 02918 03077 03236

0289 0293 03234 03595 058156 0.0926 0.10« 0.1162 0.1287 0.1416 0.1349 0.1684 0.1823 0.1964 02106 02251 02397 0.2544 0.2692 02841,

Аномальный рост значений функции Э(т) = ст(а, (3) ■ т"/р наблюдается в области значений параметров а=1, р=1. Это связано с двумя факторами: во - первых, в этой области зависимость Э-(т) от времени имеет линейный характер Э{х) и т, во-вторых, в данном случае системы вода-лед, значения параметра а(а,Р) в этой области увеличиваются.

Рис.8 Зависимость безразмерной координаты межфазной границы ^ от параметров

а, р для значения г =20. Таким образом, нелокальные уравнения теплопроводности, благодаря наличию двух новых параметров аир, которые являются показателями производных дробного порядка по времени и координате, значительно расширяют область применимости задачи Стефана, позволяя тем самым создать адекватные количественные модели процессов переноса тепла с учетом фазовых переходов. ЗАКЛЮЧЕНИЕ

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

Основные теоретические положения и практические результаты работы следующие:

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

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

3. Установлено, что учет нелокальности по времени и по пространству по разному влияют на распределение температуры вблизи источника тепла. На асимптотическое поведение распределения температуры влияет учет нелокальности по координате.

4. Построена модель переноса тепла для полупространства с учетом нелокальности по времени и координате.

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

6. Получено решение задачи о бегущих температурных волнах (задача без начальных условий). Исследована зависимость температуропроводности и параметра нелокальности от глубины на основе данных по измерению распределения температуры в верхних слоях земли.

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

Основное содержание диссертации опубликовано в работах:

1. Мейланов Р.П., Назаралиев М.А., Бейбалаев В.Д., Шабанова М.Р. Уравнение параболического типа с дифференцированием дробного порядка // Вестник Дагестанского научного центра РАН. - 2006. - С. 11-15. (из перечня ВАК)

2. Бейбалаев В.Д., Шабанова М.Р. Численный метод решения начально-граничной задачи для двумерного уравнения теплопроводности с производными дробного порядка // Вестник Самарского государственного технического университета. Сер. ФМН. -2010. №5(21). - С. 244 - 251. (из перечня ВАК)

3. Мейланов Р.П., Шабанова М.Р. Задача Стефана в дробном исчислении // Доклады Адыгской (Черкесской) Международной Академии наук. - 2010. Т. 12. №1. С. 53 - 56. (из перечня ВАК).

4. Мейланов Р.П., Шабанова М.Р. Задача Стефана // Нелинейный мир. -2011. Т.9. №7. - С. 477 - 481. (из перечня ВАК).

5. Мейланов Р.П., Шабанова М.Р. Особенности решений уравнения теплопроводности в производных дробного порядка // Журнал технической физики. -2011.Т.81, №7.-С. 1-6. (из перечня ВАК).

6. Алхасов А.Б., Мейланов Р.П., Шабанова М.Р. Уравнение теплопроводности в производных дробного порядка // Инженерно-физический журнал. — 2011. Т.84. №2. - С. 309 -317. (из перечня ВАК).

7. Мейланов Р.П., Рамазанова А.Э., Шабанова М.Р. Равновесная термодинамика систем с фрактальной структурой // Тезисы докладов XIV международная конференция по химической термодинамике. — С.-Петербург 2002.-С. 232-233.

8. Мейланов Р.П., Шабанова М.Р., Янполов М.С. Обобщенное уравнение Фок-кера-Планка в задачах тепломассопереноса // Материалы международной конференции «Возобновляемая энергетика: проблемы и перспективы». - Махачкала. - 2005.-С. 278-282.

9. Шабанова М.Р. Обобщенное уравнение теплопроводности в задачах тепло-переноса // Материалы Школы молодых ученых «Актуальные проблемы освоения возобновляемых энергоресурсов» Махачкала. - 2006. - С. 244 -249.

10. Мейланов Р.П., Назаралиев М.А., Бейбалаев В.Д., Шабанова М.Р. Обобщенная задача диффузии на полупрямой // Современные наукоемкие технологии. - 2007. - №8. - С. 82-84.

11. Мейланов Р.П., Шабанова М.Р. Уравнение теплопроводности для сред с фрактальной структурой // Современные наукоемкие технологии. - 2007. -№8. - С. 84 - 85.

12. Мейланов Р.П., Шабанова М.Р. Особенности теплопереноса в средах с фрактальной структурой // Материалы Школы молодых ученых «Нелокальные краевые задачи и проблемы современного анализа и информатики». - Нальчик-Эльбрус. - 2007. - С. 104 -109.

13. Мейланов Р.П., Шабанова М.Р. Уравнение теплопроводности в производных дробного порядка и приложение к задачам геотермии // Сборник научных трудов «Тепловое поле Земли и методы его изучения» - Москва РГГУ. - 2008. -С. 145 -150.

14. Шабанова М.Р. Инвариантные свойства уравнения теплопроводности с производными дробного порядка II Материалы Школы молодых ученых « Ак-

туальные проблемы освоения возобновляемых энергоресурсов». Махачкала. -2008.-С. 219-221.

15. Мейланов Р.П., Магомедов Р.М., Шабанова М.Р., Ахмедова Г.М. Задача без начальных условий для нелокального уравнения теплопроводности // Материалы Международного Российско-Абхазский симпозиума «Уравнения смешанного типа и родственные проблемы анализа и информатики». - Нальчик-Эльбрус. - 2009. - С. 162 - 165.

16. Назаралиев. М.А., Мейланов Р.П., Бейбалаев В.Д., Шабанова М.Р. Особенности фазовой траектории фрактального «брюсселятора» // Труды Всероссийской научной конференции с международным участием «Математическое моделирование и краевые задачи». — Самара. - 2010. - С. 204 -210.

17. Мейланов Р.П., Шабанова М.Р. Задача Стефана в дробном исчислении // Труды Всероссийской научной конференции с международным участием «Математическое моделирование и краевые задачи». Самара. -2010. - С. 192197.

18. Мейланов Р.П., Шабанова М.Р. Задача Стефана в средах с фрактальной структурой // Материалы Международного Российско-Болгарского симпозиума «Уравнения смешанного типа и родственные проблемы анализа и информатики». - Нальчик. - 2010. - С. 163 -165.

19. Мейланов Р.П., Шабанова М.Р. Особенности решений уравнения теплопе-реноса в производных дробного порядка // Материалы Международного Российско-Болгарского симпозиума «Уравнения смешанного типа и родственные проблемы анализа и информатики». - Нальчик— 2010. - С. 166-169.

20. Мейланов Р.П., Шабанова М.Р., Ахмедова Г.М. Задача Стефана на основе нелокального уравнения теплопроводности // Материалы II Международной конференции «Возобновляемая энергетика: Проблемы и перспективы». - Махачкала.-2010.-С. 160-164.

21. Бейбалаев В.Д., Мейланов Р.П., Назаралиев М-Ш.А., Шабанова М.Р. Численное решение нелокального уравнения теплопроводности // Материалы II Международной конференция «Возобновляемая энергетика: Проблемы и перспективы». - Махачкала . - 2010. - С. 221-225.

22. Мейланов Р.П., Шабанова М.Р. Об определении производной Рисса // Материалы I Всероссийской конференции молодых ученых «Математическое моделирование фрактальных процессов, родственные проблемы анализа и информатики». - Кабардино-Балкарская республика пос. Терскол. - 2010. -С.125 -129.

23. Шабанова М.Р. Аномальные решения нелокальной задачи Стефана // Материалы II Международного Российско-Казахского симпозиума «Уравнения смешанного типа, родственные проблемы анализа и информатики». — Нальчик.-2011.-С. 215-218.

Подписано в печать 02.09.2011г. Формат 60х84Ш6. Печать ризографная. Бумага офсетная. Гарнитура «Тайме». Усл. п. л. 1. Тираж 100 экз.

Отпечатано в типографии АЛЕФ, ИП Овчинников М.А. Тел.: +7-928-264-88-64, +7-903-477-55-64, +7-988-2000-164

Содержание диссертации автор исследовательской работы: кандидата технических наук, Шабанова, Муминат Руслановна

Введение.

ГЛАВА I. Современное состояние математической теории теплопроводности и нестационарных методов определения теплофизических характеристик веществ.

§1.1 Уравнение теплопроводности на основе принципа локального равновесия.

§ 1.2 Уравнение теплопроводности на основе принципа локального неравновия.

§ 1.3 Уравнение теплопроводности в производных дробного порядка.

§1.4 Нестационарные методы определения температуропроводности.

ГЛАВА II. Нелокальное уравнение теплопроводности-в дробном исчислении.

§2.1 Математическая модель нелокального переноса тепла для неограниченной прямой.

§2.2 Математическая модель нелокального переноса тепла на полуограниченной прямой.

§2.3 Определение температуропроводности и параметров нелокальности по экспериментальным данным нестационарных методов измерения температуры.

ГЛАВА III. Нелокальное уравнение теплопроводности для задачи без начальных условий.

§3.1 Классическая задача без начальных условий.

§3.2 Математическая модель нелокального переноса тепла для задачи без начальных условий.

§3.3 Определение температуропроводности по экспериментальным данным распределения температуры в поверхностных слоях земли.

ГЛАВА IV. Математическая модель задачи Стефана в дробном исчислении.

§4.1 Классическая задача Стефана.

§4.2 Математическая модель нелокальной задачи Стефана.

§4.3 Аномальные решения задачи Стефана в дробном исчислении.

Введение диссертация по физике, на тему "Теплопроводность в дробном исчислении: приложения к нестационарным методам определения теплофизических характеристик веществ и к задаче Стефана"

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

Фундаментальной ф* зической причиной необходимости учета нелокальных эффектов в гетерофазных системах является сложная природа спектра характерных времен релаксаций неравновесного состояния к равновесному состоянию, приводящая к медленной релаксации корреляционных связей, когда многочастичные функции распределения не распадаются на произведения од-ночастичных функций распределения. В результате нарушаются условия выполнения принципа локального равновесия, традиционные методы «сокращенного» описания в статистической физике становятся непригодными, поэтому необходимо исходить из принципа локального неравновесия. Исследование неравновесных процессов в условиях принципа локального неравновесия приводит к необходимости учета эффектов памяти (нелокальность по времени) и пространственных корреляций (нелокальность по координате) и развития принципиально новых методов анализа, основанных на применении математического аппарата интегродифференцирования дробного порядка - дробного исчисления.

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

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

Математический аппарат интегродифференцирования дробного порядка имеет давнюю историю, однако его интенсивное развитие началось с развития концепции фрактала, которая потребовала разработки математического аппарата интегродифференцирования дробного порядка как фундаментальной математической основы физики фракталов. Существенный вклад в развитие дробного исчисления и его прикладных аспектов внесли работы Нахушева A.M. [55,56], Нахушевой В.А. [57], Псху A.B. [64], Чукбара К.В. [87], Сербиной Л.И. [70,71], Потапова A.A. [61], Учайкина В.В. [77], Самко С.Г., Килбаса A.A., Маричева О.И. [69], Мейланова Р.П. [34-37], Mainardi F., Gorenflo R. [109], Mainardi F., Luchko Y., Pagnini G. [110], Kilbas A.A., Srivastava H.M., Tmjillo JJ. [103], Podlubny I. [120]. В настоящее время для решения дифференциальных уравнений дробного порядка используются как аналитические, так и численные методы [6,8,26].

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

В развитие нестационарных методов измерения теплофизических параметров и его приложений к разработке неразрушающего контроля теплофизических характеристик веществ существенный вклад внесли работы Тихонова А.Н., Самарского A.A. [76], Лыкова А.В.[30,31], Карслоу Г.С., [22,23], Карта-шова Е.М. [24], Вавилова В.П.[9,10], Власова В.В.[11-13], Филиппова Л.П., [81],

Филиппова П.И., Тимофеева A.M. [82]; Фокина В.М., Чернышова В.Н. [83], Фролова Н.М.[85] и др.

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

Основные задачи, решаемые в работе:

- получить фундаментальные решения нелокальных уравнений теплопроводности на основе математичеЬкого аппарата дробного исчисления;

- исследовать влияние нелокальности по времени и координате на распределение температуры при рассмотрении диффузионного и конвективного'механизмов переноса тепла;

- разработать метод определения температуропроводности и параметров нелокальности по времени и координате на основе решения нелокального уравнения теплопроводности и экспериментальных данных нестационарных методов определения распределения температуры; •

- на основе математического аппарата дробного исчисления разработать математическую модель задачи Стефана и приложить ее к системе вода — лед.

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

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

4. На- основе уравнения теплопроводности в производных дробного порядка по времени разработана математическая модель нелокального переноса тепла для задачи без начальных условий;

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

6: Предложена обобщенная модель задачи; Стефана, на основе нелокального уравнения теплопроводности в производных дробного порядка и на ее основе получена новая зависимость координаты межфазной границы от времени и параметров нелокальности повремени и координате;

7., Обнаружена область аномальной зависимости координаты. межфазноЙЕ границы от параметров нелокальности по времени и координате.

Апробация работы; Основные положения и выводы. диссертации были обсуждены на научных семинарах УРАН «ИПГ ДНЦ» и были предметом обсуждения на следующих научных мероприятиях:

1. XIV международная конференция по химической термодинамике. С.Петербург, 2002. , ' ;

2. Международная конференция «Возобновляемая энергетика: проблемы и перспективы». Махачкала, 2005.

3. Школа молодых ученых «Актуальные проблемы, освоения возобновляемых энергоресурсов». Махачкала, 2006.

5. Российская заочная конференция «Современные наукоемкие технологии». Москва, 2007:

6. Школа; молодых ученых «Актуальные проблемы освоениям возобновляемых энергоресурсов».Махачкала, 2008.

7. Международный Российско-Абхазский симпозиум. «Уравнения смешанного типа и родственные проблемы анализа: и информатики»; Нальчик-Эльбрус, 2009. . .

8. Всероссийская научная; конференция с международным: участием «Математическое моделирование и краевые задачи». Самара, 2010.

10. II Международная: конференция «Возобновляемая энергетика: Проблемы и перспективы»: Махачкала; 2010:

1 Г. I Всероссийская конференция молодых, ученых «Математическое моделирование фрактальных процессов; родственные проблемыанализа^и информатики». Кабардино-Балкарская республика, пос: Терскол, 2010: ;'.'.■. 12. II Международный Российско-Казахский симпозиум «Уравнения' смешанного типа и родственные проблемы; анализа и информатики». Кабардино-Балкарская республика, Нальчик, 2011. ,

Основное.содержание диссертации опубликовано в работах [2,7,8,38-53,88-90].

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

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

05Т

Э"ТЫ

Здесь - производная дробного порядка по времени (производная ер ч о у

Капуто), и - производные по координате (производные Рисса).

Для задачи Коши получено решение

- СО 00 / \

00 —00 со С"х"а где еа, (~0са) = У (-1)" - - функция Миттаг-Леффлера, т(^,0) - начальное о г(ап +1) условие для температуры. Дается анализ полученного решения. Рассмотрены различные случаи начального распределения температуры. В случае дельта -источника, когда т(#,о) = из (2) получаем

271 0

В частном случае, полагая в (3) а = 1, имеем т 1 о

7Г

-р СО

Г(4,т)= \dkcos\kt, - К&ут))ехр(-£>£рт). (4) тг »

Далее, полагая в (4) р = 2, у = 1 получим решение

-ехр

К 4£>т л/фкОт

Это решение при v = О совпадает с известным решением [76] . Если же в (4) положить у = 1 и Р —> 1, получим новое решение

7Г )2+(1)х)2 •

Выяснено влияние учета нелокальности по времени и координате на распределение температуры, как при отсутствии, так и при наличии конвекции. Установлено, что характер распределения температуры в начальные моменты времени при учете нелокальности по времени определяется соотношением между характерным временем т^ =4/7гО и временем т. На асимптотическое поведение распределения температуры влияет учет нелокальности по координате. Исследовано влияние нелокальных эффектов на распределение температуры для начального условия Щ,0) = Г0[1-#(£)], где #(£) - функция Хевисайда. Решение определяется выражением гп 00 со

Т&г) = + ОКл (- - № К )+ ехр(/&(% + (- +№)та)]. (5) о о

В случае отсутствия конвекции, когда v = 0 решение (5) принимает вид Т о 7С гр 00 ос т) = ^ | + (- ). о о

Отсюда при а = 1, (3 = 2 имеем 4

2л/От

Здесь = -|=г}ехр(-г2)Л - интеграл вероятностей. Для случая а = 1,р

71 о получим

-и

1 2 (О

71 У От У

Для данной задачи также существует характерное время хх~л[п/0, определяющее различное поведение распределения температуры в окрестности начала координат. При учете конвекции для случая а = 1 решение (5) принимает вид Т о

71 оо оо т) = ^ | | + - ¥тку ))ехр(- Вхкр ) (6) о о

Отсюда для Р = 2, у = 1 получаем решение, совпадающее с ранее известным решением [78]

Т(£т) =

2 2л]От

Тут-г'

Полагая в (6) а = 1, (3 —> 1, у = 1, получаем новый класс решений

Т(£г)->

Тп 2 "

1--—) л их 2

1Н—агсШ—7==) п 24Вт

Л<Ут•

Как и в предыдущих случаях, качественный характер влияния нелокальных эффектов на распределения температуры определяется соотношением между временем х и характерным временем тх, определяемое соотношением

• \1Хх 1 V

2 л/£> я £>

В параграфе 2.2 рассмотрена задача для полупространства без учета конвективного члена переноса тепла и, используя принцип Дюамеля, получено общее решение

2 л г\ ои ао

Щ, г) = -1 ¿к | У (<Г) ) зт(к?)ЕаЛ (-Пкрта ) о о

2П

ОО т п о - кх~р 71 о о К

7) где ЦСО - заданные граничное и начальное условия. Решение (7) в частных случаях совпадает с известными традиционными решениями. Для сравнения полученного решения с известными решениями рассмотрен случай, когда граничное условие постоянное: /л(т)=Т0. Полагая а = 1,р = 2, получено решение, совпадающее с ранее известными решениями [76,78]:

1 со г ехр

4 Dt

2 Л ехр v

Gf + <T)

4Z)r

2 V

1-ег/ Г V у V i4D )

1 х

Здесь erf (х) = —¡= [ехр{-z2)dz - интеграл вероятностей. Для случая а = 1, р -> 1 поv Я" о лучен новый класс решении

Вт Ит 1 Г0

2 Г £ VI arctgl

1--arctg

V л- ^Гуу

Асимптотическое поведен! е полученных решений меняется от экспоненциального до степенного характера. Далее рассмотрен практически важный случай, когда 1|/(£') = Т0, //(г) =Т0 и получено решение

Т& т) = Г0 + — (Т0 (8) п I к '

При а = 1, р = 2 решение (8) дает т) = Г0 + (Г0 - Т0 )ег/(—5 что совпадает с традиционным решением [31]. При а = 1,(3 = 1 имеем новое решение

ПЪ 1) = Г0+ (Т0 - Т0 .

Вт

В остальных случаях (0 < а < 1,1 < р < 2) имеем новый класс решений.

В §2.3 предлагается метод расчета температуропроводности на основе экспериментальные данных по определению распределения температуры для полимерного материала - полиметилметакрилат (ПММК, орг. стекло) полученных в работах [11 - 13]. Метод расчета заключается в определении значений безразмерной температуропроводности, параметров нелокальности а и (3. Для» расчета использовалась формула? = т0 +—(т0-т0)\с1к-—,(-£>&рта).

-' V ' ' ; • ; ■■ о- к '

Подставляя заданные значений ,, и экспериментально найденного значения температуры 7X4,,т,)в эту формулу, определяем значения параметров Л. а, р. В отличие от известных методов определения температуропроводности, где используются данные нестационарных методов -определения распределения температуры, в предлагаемом методе определяется не только значения температу-. ропроводности, но и его распределение по координате. Так, значение темпера-; туропроводности у торца образца; оказывается больше, чем в объеме, которое равно а = 1.12 -10(м1 / с). При этом параметры нелокальности ? а и (5 также меняются; Предлагаемый метод анализа позволяет получить более подробную информацию о температуропроводности, определяя и неоднородность его распределения вдоль образца.

Третья; глава посвящена задаче без начальных условий, которая имеет важные прикладные аспекты и встречается при изучении процессов теплопроводности в моменты времени, удаленные от начала процесса, когда; начальное условие не влияет на распределение температуры. В §3.1 кратко изложена классическая задача без начальных условий. На основе полученного во второй главе общего решения нелокального уравнения; теплопроводности для полупрямой получено решение рассматриваемой задачи, совпадающее с известным решением [76]:

Здесь граничное условие задается в.виде Г^т)^ ц(т),41(т) - некоторая заданная функция; Рассматривая; случай1 ц(т) = т0 соз(£2т), показано, что решение имеет вид/ ■

Щ, х) = Т() ехр(-^/Ш21)) сс^Пт - . (9)

В § 3.2 рассмотрена нелокальная по времени задача без начальных условий, когда нелокальное уравнение имеет вид дта дС где нелокальная производная по времени задается в виде 1 а ¡ЖеЪь дта Т(1-а)дт1(т-гУ ■

Получено общее решение в виде

Щ, т) = — }Л ■ квт(ф ехр(-ш>г) . и ^ ¿271 О-к + (ко) со

Здесь Ф/7 (со) = / ¿/х ехр(-/ют)ф(т) . оо

В случае, когда ^(г) = 7^ соз(0 г) 5 решение принимает вид т) = Г0 ехр(- /I) • соБ^та/4))-соз(от - IО • втСяа/4)) , (10) где О - частота изменения температуры

Если а = 1, то это решение переходит в известное решение (9). Дается анализ полученного решения.

Показано, что учет эффектов памяти приводит к перенормировке характерного масштаба затухания и времени запаздывания температурных волн. Для характерного масштаба затухания температурных волн имеем I = х0 л/2 • соБ^а/Ч), где = а / со - характерный масштаб затухания при отсутствии учета эффектов памяти (а=1). Для характерного времени запаздывания температурных волн имеем: /3 = ¿0 - л/2 зт(?га/4), где 70 = л/1/(ю • а2) - характерное время запаздывания температурных волн без учета эффектов памяти. Глубина проникновения тепла в почву зависит от периода колебаний температуры на поверхности и параметра нелокальности = Т0 ехр(-£л/Ш2£>) • соз(ла/ 4). При а<1 характерная длина затухания увеличивается, а время запаздывания уменьшается.

В § 3.3 рассмотрены приложения задачи без начальных условий к задачам геотермии. На основе предложенного во второй главе метода определены температуропроводность поверхности земли в области до нейтрального слоя, где сказывается влияние колебаний температуры на поверхности Земли. Полученные данные удовлетворительно согласуются с ранее известными данными. В отличие от ранее известных методов, когда определяется среднее значение температуропроводности, в предлагаемом нами методе определяется зависимость температуропроводности от глубины. Пределы изменения и характер изменения температуропроводности с глубиной совпадает с результатами других авторов.

Четвертая глава посвящена обобщению известной задачи Стефана на основе нелокальных уравнений теплопроводности в производных дробного порядка. Постоянный и повышенный интерес к задаче Стефана связан с тем, что она, удачно сочетав математические и физические проблемы, охватывает широкий круг фундаментальных проблем как математики, так и физики. С точки зрения математики она представляет продуктивную модель класса нелинейных задач. С точки зрения физики задача Стефана в ее классической постановке позволяет обобщить себя, охватывая особенности подвижной области фазовых переходов. Межфазная область - это особое состояние вещества, занимающее промежуточное положение между сосуществующими фазами, физическая природа которого до сих пор до конца не понята, особенно в условиях фазового перехода. В настоящее время как математические, так и физические аспекты задачи Стефана интенсивно развиваются. В §4.1 кратко изложена классическая задача Стефана в постановке работы [76] применительно к системе вода-лед. Зависимость координаты межфазной границы от времени имеет вид:

Цт) = ат2, где константу а необходимо определить из уравнения, которое выводится из условия Стефана. Новый этап развития задачи Стефана, которому посвящен § 4.2 , связан с ее распространением на системы со сложной структурой, включая и фрактальные структуры, где возникает необходимость учета нелокальных эффектов по времени и координате. Процессы нелокального переноса тепла при этом описываются на основе дифференциальных уравнений в производных дробного порядка. Постановка задачи Стефана в производных дробного порядка задается уравнениями нелокального переноса тепла, предложенными, во второй главе и имеют вид: дт° д^ Г(1 —а)т" ~ 0<£<$(т)

11)

5%(1,х) Г, (5.0) --—' *«<#<»

Краевые условия для системы уравнений (11) определяются соотношениями

Т2(£,т) = г2 при £ = т) ла^(т)

--— = = = (12)

Здесь <9(т) - координата подвижной границы раздела фаз. Решения задачи Стефана в производных дробного порядка имеют вид:

Ф Ф л ® а г< =+ ^), о < е < .

Г,ах) - ^) + Ь-Ь 1]лЁ!®.е

Здесь

Зависимость координаты межфазной границы от времени имеет вид т) = с(а,(3)-та/р , (13) где уравнение для определения ст дается выражением (¿з-ОХ^рСД/сгР) «2-Г3)\2Ка^Р2/<уУ) = Г(1 + а/(3) 1+у К., (А / ) (1 - Fa¿ (А / ар)) Г(1- уа / р)

Здесь ек ( где им (х,у) - функция Ломмеля от двух аргументов.

Принципиальное отличие закона $(т) = а(а, Р) • та/р от ранее известных заключается в том, что &(т) зависит от двух параметров а и р. В § 4.3. приведены результаты приложения нелокальной задачи Стефана к системе вода - лед. Проведены расчеты и анализирована зависимость координаты межфазной границы при различных а и (3 .

Для определения значения функции 0(т) = а(а, Р) • та/р установим, как зависит а(а. Э) от параметров а, Р для случая системы вода-лед.

Исходя из значений параметров температура льда /, =-5°С, воды (2 = 0°С, /3=0°С, р1[Аз = 9ПЛкг/м3,Х1Ьйа -2.24Вт/м• К, дльда=330цЦж/кг, альда = 1.2-1(Г6лГ/с уравнение для определения а(а,(3) принимает вид 27.03 + ^(1,2/^) ГО - уа / (5) '

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

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

Заключение диссертации по теме "Теплофизика и теоретическая теплотехника"

Заключение

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

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

В рамках дробного исчисления, где операция дробного дифференцирования является определенным сочетанием обычных операций интегрирования и дифференцирования, открывается новый подход к нелокальным уравнениям теплопроводности. Описание процессов переноса тепла на основе дифференциальных уравнений дробного порядка позволяет естественным образом учесть пространственную и временную нелокальности Дифференциальные уравнения в производных дробного порядка имеют более широкий, зависящий от трех параметров а , р и У, класс решений. В частном случае (а = 1,/? = 2,у = 1) полученные решения совпадают с ранее известными решениями. В остальных случаях мы имеем новый класс решений, асимптотическое поведение которых качественно отличается от ранге известных решений. Важно и то, что функциональный вид полученных решений определяется значениями показателей производных дробного порядка а, р и У, которые тем самым, становятся новыми параметрами теории для интерпретации экспериментальных данных.

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

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

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

2. Получено трехпараметрическое семейство решений-нелокального уравнения теплопроводности для неограниченной прямой; с учетом диффузионного; и конвективного механизмов перенора тепла., •

3. Установлено; что; учет нелокальности по времени и по пространству по разному влияют на распределение температуры вблизи источника- тепла: На асимптотическое поведение распределения температуры влияет учет нелокальности по координате.

4. Построена модель переноса тепла: для полупространства с учетом нёлокальности по времени и координате.

5. Разработан метод определения температуропроводности и параметров, нелокальности по времени и координате на основе решений, нелокального уравнения теплопроводности.

6. Получено решение задачи о бегущих температурных волнах (задача: без начальных условий):. Исследована зависимость температуропроводности и параметра нелокальности от глубины на основе данных по измерению распределения температуры в верхних слоях земли.

Список источников диссертации и автореферата по физике, кандидата технических наук, Шабанова, Муминат Руслановна, Махачкала

1. Александров Д.В., A.A. Иванов. Задача Стефана затвердевания трехкомпонентных систем при наличии движущихся областей фазового перехода // ЖЭТФ. - 2009. - Т. 135. - № 5. - С. 942-950.

2. Алхасов А.Б., Мейланов Р.П., Шабанова М.Р. Уравнение теплопроводности в производных дробного порядка / / ИФЖ. 2011. - Т.84. -№2. - С. 309-317.

3. Амирханов Х.И., Суетнов В.В., Левкович P.A., Гаирбеков Х.А. Тепловой режим осадочных пород. — Махачкала, 1972. 230 с.

4. Бабенко Ю.И. Метод дробного дифференцирования в прикладных задачах теории тепломассообмена. — СПб. НПО «Профессионал», 2009. 584 с.

5. Боголюбов H.H. Проблемы динамической теории в статистической физике. Гостехиздат, 1946. 120 с.

6. Бейбалаев В.Д. Численный метод решения математической модели теплопереноса в средах с фрактальной структурой. // Фундаментальные исследования. Москва- 2007.-№12.- С. 249-251.

7. Бейбалаев В.Д., Мейланов Р.П., Назаралиев М-Ш.А., Шабанова М.Р. Численное решение нелокального уравнения теплопроводности // II Международная конференция «Возобновляемая энергетика: Проблемы и перспективы». Махачкала . - 2010. - С. 221—225.

8. Бейбалаев В.Д., Шабанова М.Р. Численный метод решения начально-граничной задачи для двумерного уравнения теплопроводности с производными дробного порядка // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. Науки. 2010. №5(21). - С. 244-251.

9. Вавилов В.П. Тепловые методы неразрушающего контроля: Справочник. М.: Машиностроение, 1991. 240 с.

10. Ю.Вавилов В.П. Инфракрасная термография и тепловой контроль. — М.: ИД Спектр, 2009. 544 с.

11. Власов А.Б. Расчет эксплуатационных показателей надежности контактных соединений с помощью тепловизионного контроля // Электротехника, №8. — 2002. С. 30-35.

12. Власов А.Б. Исследование нестационарных тепловых процессов в диэлектрике с помощью тепловизора // Вестник МГТУ. Труды Мурманского государственного технического университета, 2003. — Т.6. №1. С. 29 -34.

13. Головизнин В.М., Киселев В.П., Короткин И.А., Юрко Ю.И. Прямые задачи неклассического переноса радионуклидов в геологических формациях // Известия РАН. Энергетика. 2004. -№4. - С. 121-130.

14. Гуров К.П. Феноменологическая термодинамика необратимых процессов. М.: Наука, 1974. 128 с.

15. Данилюк И.И. О задаче Стефана. // УМЫ.- 1985.-Т.40-№5(245).- С. 133-185.

16. Де Гроот С., Мазур П. Неравновесная термодинамика. М.: Мир, 1964. 456 с.

17. Дьярмати И. Неравновесная термодинамика. М.Мир, 1974. 404 с.

18. Карслоу Г.С. Теория теплопроводности. М. Л. ОГИЗ, 1947. - 288 с.

19. Карслоу Г., Егер Р. Теплопроводность твердых тел. М.: Наука, 1964. — 488 с

20. Карташов Е.М. Аналитические методы в теории теплопроводности твердЬ1х тел. Высшая школа, 2001. 550 с.

21. Климонтович Ю.Л. Введение в физику открытых систем. М.: Яцус^ 2002. - 284 с.

22. Кольцова Э.М., Василенко В.А., Тарасов В.В. Численные методы рещеНвд уравнений переноса во фрактальных средах // ЖФХ, 2000. Т.74. №5. q 954-956.

23. Куни Ф.М. Статистическая физика и термодинамика. -М.:Наука,1981. — 352 с.

24. Курдюмов С.П., Курленка Е.С. Тепловые структуры в среде с нелинейной теплопроводностью // Нелинейный мир. — 2005. Т.З. - С. 5-6.

25. Ландау Л.Д, Лифшиц Е.М. Гидродинамика. М.: Наука, 1986. 796 с.

26. Лыков A.B. Тепломассообмен (справочник). М. «Энергия». 1971. —560 с

27. Лыков A.B. Теория теплопроводности М.: ВШ, 1967. - 600с.

28. Мандельброт Б. Фрактальная геометрия природы М.: Наука, 2002. — 654 с

29. Мейерманов A.M. Задача Стефана. Новосибирск: Наука. 1986. - 240 с.

30. Мейланов Р.П. К теории фильтрации в пористых средах с фрактальной структурой // Письма ЖТФ. 1996. - Т.22. - № 23. - С.40-43.

31. Мейланов Р.П. Обобщенное уравнение одномерной фильтрации с дифференцированием дробного порядка // ИЖФ- 2001.-Т.74.-№2.— С.34

32. Мейланов Р.П., Янполов М.С. Особенности фазовой траектории «фрактального» осциллятора // Письма в ЖТФ 2002.-T.28.-JNbl-C.67—73

33. Мейланов Р.П., Свешникова Д.А., Шабанов О.М. Метод дифференциальных уравнений дробного порядка в описании кинетики сорбции // ЖФХ. 2000, — Т.77. №2. С. 260 - 264.

34. Мейланов Р.П., Рамазанова А.Э., Шабанова М.Р. Равновесная термодинамика систем с фрактальной структурой // XIV международная конференция по химической термодинамике. — С.-Петербург 2002 .

35. Мейланов Р.П., Шабанова М.Р., Янполов М.С. Обобщенное уравнение Фоккера-Планка в задачах тепломассопереноса // Возобновляемая энергетика: проблемы и перспективы. — Махачкала — 2005 — С. 278—282.

36. Мейланов Р.П., Назаралиев М.А., Бейбалаев В.Д., Шабанова М.Р. Уравнение параболического типа с дифференцированием дробного порядка // Вестник ДНЦ РАН. -2006.- С. 11-15.

37. Мейланов Р.П., Шабанова М.Р. Особенности теплопереноса в средах с фрактальной структурой // Школа молодых ученых « Нелокальные краевые задачи и проблемы современного анализа и информатики» Нальчик-Эльбрус. - 2007. - С. 104-109.

38. Мейланов Р.П., Назаралиев М.А., Бейбалаев В.Д., Шабанова М.Р. Обобщенная задача диффузии на полупрямой // Современные наукоемкие технологии. 2007. -№8.- С. 82-84.

39. Мейланов Р.П., Шабанова М.Р. Уравнение теплопроводности для сред с фрактальной структурой // Современные наукоемкие технологии. 2007. -№8.- С.84-85.

40. Мейланов Р.П., Шабанова М.Р. Уравнение теплопроводности в производных дробного порядка и приложение к задачам геотермии // Сборник научных трудов «Тепловое поле Земли и методы его изучения» -Москва РГГУ. 2008. - С. 145-150.

41. Мейланов Р.П., Магомедов P.M., Шабанова М.Р., Ахмедова Г.М. Задача без начальных условий для нелокального уравнения теплопроводности //

42. Международный Российско-Абхазский симпозиум. «Уравнения смешанного типа и родственные проблемы анализа и информатики». Нальчик-Эльбрус. -2009.-С. 162-165. .

43. Мейланов Р.П., Шабанова М.Р. Задача Стефана в дробном исчислении. // Труды * Всероссийской научной конференции? с международным участием « Математическое моделирование и краевые задачи» — Самара. — 2010 г. —1. С. 192-197. ;

44. Мейланов Р.П., Шабанова М;Р. Задача Стефана в средах с фрактальной структурой // Международный Российско-Болгарский симпозиум «Уравнения смешанного типа и родственные проблемы анализа и информатики» Нальчик. -201 ().- С. 163-165. .

45. Мейланов Р.П., Шабанова М.Р., Ахмедова Е.М; Задача Стефана на основе нелокального уравнения теплопроводности // Материалы II';Международной конференции «Возобновляемая энергетика: Проблемы и . перспективы». — Махачкала . 2010. - С. 160-164.

46. Мейланов Р.П., Шабанова М.Р. Задача Стефана в дробном исчислении // Доклады Адыгской (Черкесской) Международной Академии наук. 2010. Т. 12. №1. С. 53-56.

47. Мейланов Р.П., Шабагова М.Р. Задача Стефана // Нелинейный мир. -2011. Т.9. №7. С. 477 — 481.

48. Мейланов Р.ГГ., Шабанова М.Р., Особенности решений уравнения теплопереноса в производных дробного порядка // ЖТФ. 2011.— Т.81. №7 -С. 1-6.

49. Нахушев A.M. Элементы дробного исчисления и их применение. Нальчик • издательство КБНЦ РАН, 2000. 299 с.

50. Нахушев A.M. Дробное исчисление и его применение. М.: ФИЗМАТЛИТ 2003.-272 с.

51. Нахушева В.А. Дифференциальные уравнения математических моделей нелокальных процессов // М., 2006. 174 с.

52. Нигматулин P.P. Дробный интеграл и его физический смысл // ТМф. 1992 -Т.90.№3.- С. 354-368.

53. Новиков С .В. Тепловые свойства терригенных коллекторов и насыщающих флюидов // Автореферат диссертации на соискание ученой степени кандидата технических наук. Москва, 2009.

54. Олемской А.И., Флат А .Я. Использование концепции фрактала в физике конденсированной среды // Успехи физических наук. -1993.- Т. 163 №12-С. 1-50.

55. Потапов A.A. Фракталы в радиофизике и радиолокации: Топология: выборки. М.: Университетская книга, 2005. — 848 с.

56. Подстригач Я.С., Коляно Ю.М. Обобщенная термодинамика. Киев. Наукова Думка, 1976. - 310 с.

57. Пригожин И. Введение в термодинамику необратимых процессов. М.: ИЛ, 1960.- 160 с.

58. Псху A.B. Уравнения в частных производных дробного порядка. М.: Наука, 2005. 199 с.

59. Пулькина JI.C. Краевые задачи с нелокальными граничными условиями для уравнений высокого порядка // Неклассические уравнения математической физики. Сборник научных работ. Новосибирск, 2010. С . 220 -232.

60. Рутман С.С. К статье P.P. Нигматулина «Дробный1 интеграл и его физический смысл» // ТМФ. -1994. -Т. 100. №3. С. 476-478.

61. Рутман С.С. О физических интерпретациях фрактального интегрирования и дифференцирования // ТМФ. -1995. ТЛ05. №3. - С. 393-^104.

62. Самарский A.A., Галактионов В.А., Курдюмов С.П., Михайлов А.П. Режимы с обострением в задачах для квазилинейных параболических уравнений. М.: Наука, 1987. - 480 с.

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

64. Сербина Л.И. Нелокальные математические модели переноса в системах с фрактальной структурой. Нальчик: Издательство КБНЦ РАН, 2002. - 144 с.

65. Сербина Л.И. Нелокальные математические модели переноса в водоносных системах. М.: Наука. 2007. — 167 с.

66. Силин В.П. Введение в кинетическую теорию газов. М.: Наука, 1971 г. 339 с.

67. Соболев С JÍ. Локально неравновесные моделишроцессов переноса // УФН. -1997. -Т.167. № 10.-С. 1095-1106.

68. Станиславский А.А. Вероятностная интерпретация интеграла дробного ; порядка // ТМФ. 20041- Т.138. №3. -С.491-507.

69. Суегнов В В. О некоторых закономерностях изменения температур в приповерхностных участках земной коры // Труды Института геологии Дагестанского филиала АН СССР. Вопросы гидрогеологии и геотермии Дагестана. Выпуск.5.С. 42- 46. • "

70. Тихонов А.Н., Самарский А.А. Уравнения математической физики; М.: i Наука. 1972.-736 с. \

71. Учайкин В В. Метод дробных производных: Ульяновск: Издательство «Артишок», 2008; -512 с. , ' ; .

72. Фарлоу С. Уравнения в частных производных. М.: Мир, 1985. 382 с.

73. Федер Е. Фракталы. М.: Мир, 1991; 254 с.

74. Ферцигер Дж., Г. Капер. Математическая теория процессов переноса в газах. М.: Мир, 1976 . 554 с. . . : г - ;

75. Фокин В;М;, Чернышов В.И. Неразрушающий контроль теплофизических характеристик . строительных" материалов. М.: "Издательство Машиностроение-!", 2004.212 с.

76. Фридман А. Вариационные принципы и задачи со свободными границами. М.: Наука, 1990. 536 с.

77. Фролов НМ. Температурный режим гелиотермозоны. М., 1966. -156 с.

78. Череменский Г.А. Прикладная геотермия. Л., «Недра», 1977. — 224 с.

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

80. Шабанова М.Р. Инвариантные свойства уравнения теплопроводности с производными дробного порядка // Школа молодых ученых «Актуальные проблемы освоения возобновляемых энергоресурсов» Махачкала — 2008.-С. 219-221.

81. Шабанова М.Р. Обобщенное уравнение теплопроводности в задачах теплопереноса // Школа молодых ученых «Актуальные проблемы освоения возобновляемых энергоресурсов». Махачкала - 2006 — С. 244—249.

82. Шабанова М.Р. Аномальные решения нелокальной задачи Стефана // Материалы II Международного Российско-Казахского симпозиума «Уравнения смешанного типа, родственные проблемы анализа и информатики». Нальчик. — 2011. — С. 215 — 218.

83. Эккерт Э.Р., Дрейк P.M. Теория тепло — и массообмена. М.: Госэнергоиздат, 1961.-680 с.

84. Cattaneo С. Sulla conduzione del calore // Atti Sem. Mat. Fis. Univ. Modena. -1948.-2. P. 83-101.

85. Cattaneo C. Sur une forme de l'équation de la chaleur éliminant le paradoxe d'une propagation instantanée // C.R. Acad. Sci. 1958. - 247. №.4. - P. 411433.

86. Charles Tadjeran, Mark M. Meerschaert, Hana-Peter Scheffler. A second-order accurate numerical approximation for the fractional diffusion equation.- Journal of Computational Physics 213 ( 2006 ) -P. 205-213.

87. Chen P J., Gurtin M.E. On second sound in materials with memoiy // ZA. Angew. Math. Phys. 1970. - 21, No. 2. - P. 232-241

88. Day W. The thermodynamics of simple materialist with fading memory. -Berlin: Springer. 1972. 134 p.

89. Fujita Y. Integrodifferential equation which interpolates the heat equation and the wave equation 11 Osaka J. Math. 1990. - 27, No.2. - P.309-321.

90. Germant A. On Fractional Differential Equation // Philosophical Magazine. — 1938.-V. 25. P. 540-549.

91. Gorenflo R., Mainardi F., Moretti D., Paradisi P. Time fractional diffusion: a discrete random walk approach // Nonlinear Dynamics. — 2002. 29, No. 1-4. — P.129- 143.

92. Gorenflo R., Iskenderov A., Luchko Y. Mapping between solutions of fractional diffusion wave equations // Fractional Calculus Appl.Anal. - 2000. — 3, No. l.-P. 75-86.

93. Gurtin M.E., Pipkin A.C. A general theory of heat conduction with finite wave speeds // Arch. Rational Mech. Anal. 1968. - 31, No.2. - P. 113-126.

94. Hristov J. Heat-Balance integral to Fractional (hale-time) Heat Diffusion submode 1// Journal Thermal Stresses 2005. V.14, N.2. - P. 291- 316.

95. Kilbas A.A., Srivastava H.M., Trujillo J.J. Theory and Applications of Fractional Differential Equations. North-Holland Mathematics Studies. Vol. 204. Amsterdam, etc.: Elsevier. 2006. 523 p.

96. Kulish V.V., Lage Fractional-Diffusion Solutions for Transient Local Temperature and Heat Flux // J. Heat Transfer . 2000 V.122, Issue 2, P. 372377

97. Lame G., Clapeiron B.P. Memoire sur la solidification par refroidissement d'un globe solid. Ann. De Chem. Et de Phys. 1831. t. 47. - P. 250 - 256.

98. Lenzi E.K., L.R. da Silva, A.T. Silva, L.R. Evangelista, M.K. Lenzi. Some results for a fractional diffusion equation with radial symmetry in a confined region // Physical A: 2009. V387. P. 806-810.

99. Liu Junyi, Xu Mingyu. Some exact solutions to Stefan problems with fractional and Applications // Journal of Mathematical Analysis and Applications. 2009. V.351. -P. 536-542.

100. Liu Q., Liu F., Turner I., Anh V. Approximation of the Levy-Feller^advection-dispersion process by random walk and finite difference method // Journal of Computational Physics 222 .2007. P. 57-70.

101. Mainardi F and5 Gorenflo R. On Mittag-Leffler-type functions;" .in.' fractional evolution processes// J. Comput. Appl. Math. 118. 2000. P. 283 - 299.

102. Meerschacrt M.M., Tadjeran C. Finite difference approximations for two-sided space-fractional partial differentia 1 equations 1/ Applied Numerical Mathematics 56. 2006. P. 80 - 90.

103. Meilanov R.P., Sveshnikova D.A„ Shabanov O.M. Fractal Neuter of SOrption Kinetics // Journal of Physical Chemistry A. 2002. V/106. P. 11771 - 11774.

104. Moodi T.B., Tait R.J. On theraial transients with Unite wave speeds // Acta Mech. 1983. - 50, No. ./2. - P. 97-104.

105. Nigmatullin R.R. To the theoretical explanation of the «universal response» // ; Physr Stat. Sol. (b). -1984. -.123, No! 2.— P. 739-745:

106. Norwood F.R: Transient then-nal waves in the general theory of heat conduction with finite wave speeds //J. Appl. Mech. -1972. 39, No. 3. -P. 673-676. "•

107. Nunziato J.W. On heat conduction in materials with memory // Quart. Appl. Math.-1971.-29,N0.2.-P. 187-204.

108. Oldham Keith B., Spanier Jerome, fhe fractional calculus (Theory and Applications of Differentiation and Integration to Arbitrary Order). Academic Press, New York and London, 1974. 234 p.

109. Parker W.J., Jenkins R.J., Butler C.P., Abbot G.L. Flash method of determining thermal diffusivity, heat capacity and thermal conductivity. J. Appl. Physics, Sept. 1961. Vol. 32. P. 1679 — 1684.

110. Podlubny I. Fractional Differential Equations. New-York: Academic Press, 1999. -340 p.

111. Povstenko Y.Z. Fractional heat conduction equation and associated thermal stresses // Journal Thermal Stresses 2005. V.28, N.l. - P. 83 - 102.

112. Povstenko Y.Z. Thermo elasticity that uses fractional heat conduction equation // Journal of mathematical Sciences. 2008.-V.52, N2, P. 239 - 246.

113. Povstenko Y.Z. Fractional Cattaneo-Type Equations and Generalized Thermo elasticity // Journal Thermal Stresses. 2011. V.34, N.2. P. 97- 114.

114. Saichev A.I., Zaslavsky G.M. Fractional kinetic equations: solutions and applications // Chaos. 1997. 7. No. 4 - P. 753 - 764

115. Stefan J. Uber einige Probleme der Theorie der Warmeleitung .- S.B. Wien. Akad. Mat, Natur., 1889, B. 98, S. 473-484.

116. Stefan J. Uber die Diffusion von Sauren und Basen gegen einander.- S.B. Wien. Akad. Mat, Natur., 1889, Bd. 98, S. 614 634.

117. Stefan J. Uber die Theorie der Eisbildung, insbesondere über die Eisbildung in Polarmeere .- S.B. Wien. Akad. Mat, Natur., 1889, Bd. 98, S. 965-983.

118. Stefan J. Uber die Verdampfung und die Auflösung als Vorgange der Diffusion .- S.B. Wien. Akad. Mat, Natur., 1889, Bd. 98, 8. 1418-1442.

119. Xiaoyun Jiang, Mingyu Xu. The time fractional heat conduction equation in the general orthogonal curvilinear coordinate and the cylindrical coordinate systems.//Physica A: Statistical Mechanics and Applications. 2010. V. 389. Issue 17. -P. 3368-3374.