Размещено на http://www.allbest.ru/

Размещено на http://www.allbest.ru/

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

Автореферат

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

Общая характеристика работы

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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. I Всероссийская конференция молодых ученых «Математическое моделирование фрактальных процессов, родственные проблемы анализа и информатики». Кабардино-Балкарская республика, пос. Терскол, 2010.

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

Публикации. Основное содержание диссертации опубликовано в 23 работах, из которых 6 - статьи в научных рецензируемых журналах из перечня ВАК

Структура и объем диссертации. Диссертация состоит из введения, четырех глав, заключения и списка литературы из 129 наименований. Объем работы 120 стр. в том числе 30 рисунков и 2 таблицы.

Основное содержание работы

диффузионный конвективный температуропроводность стефан

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

Первая глава посвящена изложению современного состояния математической теории теплопроводности. В основе анализа лежит классификация уравнений теплопроводности с позиций принципов неравновесной термодинамики - принципов локального равновесия и локального неравновесия. Кратко изложены существующие подходы в теории теплопроводности. Отмечается, что традиционный подход учета нелокальных свойств с помощью интегрального оператора имеет принципиальные трудности при практическом применении. Новым этапом развития математических основ теории теплопроводности стала теория нелокальных уравнений теплопроводности на основе математического аппарата интегродифференцирования дробного порядка - дробного исчисления [Oldham Keith B., Spanier Jerome. 1974; Самко С.Г., Килбас А.А., Маричев О.И. 1987; Нахушев А.М. 2003; Псху А.В. 2005; Потапов А.А. 2005; Нахушева В.А. 2006; Сербина Л.И. 2007; Учайкин В.В. 2008; Бабенко Ю.И. 2009].

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

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

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

(1)

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

,

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

,

.

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

, (2)

где - функция Миттаг-Леффлера,

- начальное условие для температуры. Решение (2) в частном случае совпадает с решением работы [Mainardi F., Luchko Y., Pagnini G. 2001]. Дается анализ полученного решения. Выяснено влияние учета нелокальности по времени и координате на распределение температуры, как при отсутствии, так и при наличии конвекции. Установлено, что характер распределения температуры в начальные моменты времени при учете нелокальности по времени определяется некоторым характерным временем.

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

(3)

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

. (4)

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

. (5)

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

. (6)

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

Рис. 1. Распределение безразмерной температуры по координате для различных значений параметров б, в; D=1. а) без конвекции V=0; б) с конвекцией V=1

Показано, что существует характерное время, определяющее характер влияния учета нелокальности на распределение температуры: при при . Характер распределения температуры по координате, определяемый решениями (5) и (6) зависит от соотношения и . На рисунке 2а вблизи начала координат при уменьшении значения параметра в величина температуры увеличивается при и уменьшается при . Однако распределение температуры вблизи начала координат для случаев и качественно отличается. В области значений времени как видно на рисунке 2б, зависимость распределения температуры вблизи начала координат от учета нелокальности по координате и времени меняется наоборот: при учете нелокальности по координате величина температуры вблизи начала координат уменьшается.

а) б)

Рис. 2. Зависимости безразмерной температуры от

а) , б)

Такой характер зависимости распределения температуры вблизи начала координат от соотношения между и объясняется природой сингулярности функций (5) и (6) при , когда . Для решения (5) ; для решения (6) .

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

, (7)

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

[Тихонов А.Н., Самарский А.А. 1972, Фарлоу С. 1985, Лыков А.В. 1967]

Рассматривая случай и , получим

,

где - интеграл вероятностей. Если же , то имеем новое решение:

.

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

, (8)

которое при совпадает с известным решением

.



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

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

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

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

Характерные параметры равны:

,

Температуропроводность

; .

Рис. 4. Расчетные значения температуропроводности и параметров нелокальности б, в для ф =3420 с. Сплошная кривая: зависимость D(о); пунктирная - б(о), точечная - в(о)

Как следует из рисунка 4, значение температуропроводности у торца образца больше чем в объеме, которое оказывается равным м2/с. При этом видно, что параметры нелокальности б и в также меняются.

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

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

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

, (9)

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

.

Здесь - граничное значение температуры. На практике часто возникают задачи, когда граничное условие задается в виде . Представляя интеграл (9) в виде

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

,

,

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

.

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

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

,

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

.

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

, (10)

где . Рассмотрен случай, когда , где - частота изменения температуры. Окончательно (10) принимает вид

. (11)

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

.

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



Рис. 5. Экспериментальные и расчетные данные распределения температуры

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

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

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

, (12)

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

Рис. 6 Расчетные значения температуропроводности (D(о) - сплошная кривая) и параметра нелокальности по времени (б (о) - точечная кривая) по данным режимных наблюдений (Избербаш 10.05.1959).

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

Из рисунка 6 видно, что значение температуропроводности удовлетворительно согласуется с данными других авторов. Однако, в отличие от традиционного расчета, в данном случае мы рассчитываем зависимость температуропроводности и от координаты. Температуропроводность вычисляется по формуле . В нашем случае параметр , для характерного параметра времени, поскольку рассматриваются годичные колебания, равен году. Таким образом м2/с. В соответствии с расчетами получается, что температуропроводность меняется в пределах (м2/с), что соответствует результатам других авторов [Новиков С.В. 2009]. Кроме того, температуропроводность в данном случае увеличивается с глубиной. Как видно на рис. 6. с глубиной меняется параметр нелокальности по времени б. Это соответствует тому, что свойства почвы изменяются вглубь.

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

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

, (13)

. (14)

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

, ,

при,

при .

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

для,

для .

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

(15)

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

,

.

Граничные условия:

, ,

при ,

при,, (16)

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

), ,

, ,

.

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

, (17)

где уравнение для определения дается выражением

.

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

Таким образом, при учете эффектов памяти и пространственных корреляций имеет место обобщенный закон изменения координаты области фазового перехода, определяемый соотношением . Величина оказывается зависящей от двух параметров: . При б=1 и в=2 уравнение (17) совпадает с классическим уравнением. В случае и в=2 уравнение (17) совпадает с уравнением работы [Liu Junyi, Xu Mingyu. 2009].

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

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

,, ,

, уравнение для определения принимает вид

Численные значения приведены в виде таблицы для нескольких значений б, в. В таблице по строке меняются значения , по столбцу меняются значения с шагом 0.05. Исходя из этих значений и определяя для фиксированных и , можно определить зависимость координаты от б и в. На рисунке 7а приведены результаты расчета зависимости от значений б для некоторых фиксированных значений в. Для значения в=1.8 эта зависимость имеет монотонно возрастающий характер при уменьшении значения параметра б. Для значения в=1.05 эта зависимость становится немонотонной и имеет пороговый характер. На рисунке 7б приведены результаты расчета зависимости от значений в для некоторых фиксированных значений б. Как следует из рисунка, характер зависимости от значений в определяется значением параметра б. Для значения б=0.7 эта зависимость имеет монотонно убывающий характер при уменьшении значения параметра в. Для значений б=1 эта зависимость становится немонотонной и имеет пороговый характер. Общая зависимость от параметров б, в приводится на рисунке 8. Как видно, вблизи области значений б=1,в=1 значения функции аномально растут.

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

а) от параметра б; б) от параметра в.

Таблица значений

Аномальный рост значений функции наблюдается в области значений параметров б=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. - C. 278 - 282.

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

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

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

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

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

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

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

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

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

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.

Размещено на Allbest.ru

...