Нелокальные математические модели процессов переноса в водоносных природных системах с фрактальной структурой

Сербина, Людмила Ивановна

Математическое моделирование, численные методы и комплексы программ

автореферат диссертации по информатике, вычислительной технике и управлению, 05.13.18, диссертация на тему:Нелокальные математические модели процессов переноса в водоносных природных системах с фрактальной структурой

доктора физико-математических наук: Сербина, Людмила Ивановна
город: Таганрог
год: 2004
специальность ВАК РФ: 05.13.18

Диссертация по информатике, вычислительной технике и управлению на тему «Нелокальные математические модели процессов переноса в водоносных природных системах с фрактальной структурой»

Автореферат диссертации по теме "Нелокальные математические модели процессов переноса в водоносных природных системах с фрактальной структурой"

На правах рукописи

Сербина Людмила Ивановна

Нелокальные математические модели процессов переноса в водоносных природных системах с фрактальной структурой

05.13.18 - «Математическое моделирование, численные методы и комплексы программ»

Автореферат диссертации на соискание ученой степени доктора физико-математических наук

Таганрог - 2004

Работа выполнена в Научно-исследовательском институте прикладной математики и автоматизации Кабардино-Балкарского научного центра Российской академии наук

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

профессор Нахушев Адам Маремович

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

профессор Вабищевич Петр Николаевич;

доктор физико-математических наук, профессор Елеев Валерий Абдурахманович;

доктор физико-математических наук, профессор Сухинов Александр Иванович

Ведущая организация - Институт математики им. С В. Соболева Сибирского отделения Российской академии наук

Зашита состоится ____0..&__ 2004 г. в «-/-<£-» часов на заседа-

нии диссертационного совета Д 212.259.03 при Таганрогском государственном радиотехническом университете по адресу: 347928, г. Таганрог, пер. Некрасовский 44, корпус Д, ауд. Д-406.

С диссертацией можно ознакомиться в научной библиотеке Таганрогского государственного радиотехнического университета по адресу: ул. Чехова 22. корпус А.

Автореферат разослан __^^ 2004 г.

Ученый секретарь диссертационного совета, д.т.н., проф. А.Н. Целых

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

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

П. Я. Полубаринова-Кочина, С. Ф. Аверьянов, М. Г. Андерсен, Т. П. Берт, Я. Бэр, Л. Дакштейн, Д. 3. Заславски, С. Ирмей, Р. Дж. Ханке и другие в своих исследованиях обращали внимание на важность выяснения физико-математических основ процессов переноса в природных водоносных системах и создания их математической модели для решения задач гидрогеологического прогнозирования, защиты подземных вод от загрязнения и борьбы с засолением орошаемых земель.

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

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

¥№?Шй£|№№1|Я9Г

БИБЛИОТ8КА

¡?жиз

|1Л1

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

Исследование по теме докторской диссертации проводилось в рамках плана научно-исследовательских работ Научно-исследовательского института прикладной математики и автоматизации Кабардино-Балкарского научного центра Российской академии наук по следующим двум направлениям фундаментальных исследований:

— «Нелокальные дифференциальные операторы основных и смешанных типов и их применение к фундаментальным исследованиям в различных областях знаний» (№ГР 01.20.00 12841);

— «Развитие дробного исчисления и анализа на фракталах для разработки математических моделей физико-биологических процессов и сред с фрактальной структурой» (№ГР 01.20.00 12845);

а также проекта № 00-01-00311 - «Исследование класса задаваемых дифференциальными операторами дробного порядка математических моделей тепломассопереноса в средах с фрактальной структурой», поддержанного в 2000, 2001, 2002 годах Российским фондом фундаментальных исследований.

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

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

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

1. Создана нелокальная математическая модель движения грунтовых

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

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

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

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

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

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

Эти шесть научно обоснованных результатов выносятся на публичную защиту, последний из них включен в «Отчет о деятельности Российской академии наук» в 2002 году как один из основных результатов в области математического моделирования.

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

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

Апробация работы. Основные положения и выводы диссертации были предметом систематического обсуждения на еженедельных заседаниях научно-исследовательского семинара по современному анализу и информатике Института прикладной математики и автоматизации Кабардино-Балкарского научного центра Российской академии наук и прошли апробацию на следующих научных мероприятиях:

1. Научно-техническая конференция профессорско-преподавательского состава СевКавГТУ (Ставрополь, 1992 г.);

2. Первая Международная конференция "Нелокальные краевые задачи и родственные проблемы математической биологии, информатики и физики" (Нальчик, 1996 г.);

3. Третья Международная конференция "Математика, компьютер, образование" (Москва, 1996 г.);

4. Международная конференция "Спектральная теория дифференциальных операторов и смежные вопросы" (Стерлитамак, 1998 г.);

5. Международная конференция "Воздействие интенсивных потоков на вещество" (Терскол, 1999 г.);

6. Третья региональная научно-техническая конференция "Вузовская наука- Северо-Кавказскому региону" (Ставрополь, 1999 г.);

7. Научно-техническая конференция профессорско-преподавательского состава СевКавГТУ (Ставрополь, 1999 г.);

8. Международная конференция "Аналитические методы анализа и дифференциальных уравнений. AMADE" (Минск, 1999 г.);

9. Воронежская весенняя математическая школа "Понтрягинские чтения" "Современные методы в теории краевых задач" (Воронеж, 1999 г.);

10. Третий Всероссийский симпозиум "Математическое моделирование и компьютерные технологии", "Математическое моделирование и вычислительный эксперимент в естественных и гуманитарных науках" (Кисловодск, 1999 г.);

11. Девятый Международный симпозиум "Методы дискретных особенностей в задачах математической физики" (Орел, 2000 г.);

12. Четвертый Всероссийский симпозиум "Математическое моделирование и компьютерные технологии", "Математическое моделирование и вычислительный эксперимент в естественных и гуманитарных науках" (Кисловодск, 2000 г.);

13. Всероссийская научная конференция "Математическое моделирование в научных исследованиях" (Ставрополь, 2000 г.);

14. Четвертый Сибирский конгресс по прикладной и индустриальной математике (ИНПРИМ) (Новосибирск, 2000 г.);

15. Международная конференция "Аналитические методы анализа и дифференциальных уравнений. AMADE" (Минск, 2001 г.);

16. Международная конференция "Новые подходы к решению дифференциальных уравнений" (Дрогобыч, 2001 г.);

17. Воронежская весенняя математическая школа "Понтрягинские чтения" "Современные методы в теории краевых задач" (Воронеж, 2001 г.);

18. Международная конференция "Математическое моделирование, статистика и информатика в современном управлении экономикой" (Самара, 2001 г.);

19. Вторая Международная конференция "Нелокальные краевые за-, дачи и родственные проблемы математической биологии, информатики и физики" (Нальчик, 2001 г.);

20. Десятая Международная конференция "Математика. Экономика. Образование" (Ростов-на-Дону, 2002 г.);

21. Четвертая Международная конференция по неравновесным процессам в соплах и струях (Санкт-Петербург, 2002 г.);

22. Российско-Узбекский симпозиум "Уравнения смешанного типа и родственные проблемы анализа и информатики" (Нальчик-Эльбрус, 2003 г.);

23. Международная конференция "Аналитические методы анализа и дифференциальных уравнений. AMADE" (Минск, 2003 г.);

24. Десятая Всероссийская школа-семинар "Современные проблемы математического моделирования" (Ростов-на-Дону, 2003 г.).

Публикации. По теме диссертации опубликованы 33 научных работы [1]—[33], в том числе одна монография [25] объемом 144 страницы, которая рекомендована к изданию Ученым советом Научно-исследовательского института прикладной математики и автоматизации и Президиумом Кабардино-Балкарского научного центра Российской академии наук.

Структура и объем работы. Диссертация состоит из введения, пяти глав, разбитых на параграфы, заключения и списка литературы, включающего 101 наименование. Всего 23 параграфа: 1.1.-1.4; 2.1-2.8; 3.1-3.3; 4.1-4.3; 5.1-5.5. Объем работы - 183 стр, включая иллюстрации.

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

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

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

верхностью г — Н(х, у, г) в области П евклидовой плоскости точек (х, у) с абсциссой х и ординатой у в момент времени г от начального г = 0 до расчетного г = Т, а также разработке и анализу методов их линеаризации.

В этой главе, как и во всей диссертации, существенно используются элементы дробного исчисления и в первую очередь операторы и дЦ. Оператор .[Уд, есть оператор дробного интегро-дифференцирования порядка |а| в смысле Римана-Лиувилля:

где Г(х) - гамма-функция Эйлера, [а] - целая часть числа а, которая удовлетворяет неравенству [а] < а < [а] + 1. Если п — 1 < а < л = 1,2,..., то по определению величина у, г) =

_ ту>-пдаФ,У,т)

г порядка а :или производной по М. СарШ:о.

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

В параграфе 1.1 на основе интерпретации почвогрунта как сил-пористой среды с распределенными параметрами, структура пор ье рой образует множество с фрактальной размерностью <1], и гипоте том, что за промежуток времени происходит изменение уровня г товых вод к(х,у,Ь) в точке (х,у) в момент времени £ на 8Ь. =

Л, = д£Л(:с,у,т),

получено нагруженное обобщенное уравнение Буссинеска

дЦ1[сгк(х, у, т)] = Ь(к - hJ к) + Ю0 +

Здесь ц - число из интервала ]0,1[; о — водоотдача или недостаток: сьпцения грунта;

г,т .х д (. _дк\ д Л тЭН\

ка

- дифференциальное выражение; гу0 = ~гг (Но — И) - скорость прос

Л/о

вания через водоупор, к - коэффициент фильрации; ко — коэффицис фильтрации слабопроницаемого водоупора, который значительно м-ше к; Но - напор в нижележащем водоносном пласте; Мо = мощность водоупора; го = ъи(х, у, ¿) - величина инфильтрации (или 5 ность между инфильтрацией и испарением) на единицу площади в ницу времени.

Для широких классов пористых сред можно приближенно положи

Ц = <1}.

Уравнение (2) эквивалентно уравнению

у, г)] = ^¿(1; Л2) - £(Ао; Л) - с^кк + Д

с Со = к0/(кМ0), Ь = + ксоЩ. Поскольку

<1)

(2)

тта-

<3)

/1-и

дк(х, у, г)

К '

то уравнение (2) при д —»■ 1 совпадает с обобщенным (по терминолоз П.Я. Полубариновой-Кочиной) уравнением Буссинеска

даИ

~дГ

д („дК\ д (Т,дК\ „ г

-СО

^г^де К = k(h — ho) - проводимость пласта, /0 = сок(Но — hо).

Пусть а, к, Мо - постоянные величины, водоупор - горизонтален и = 0. Тогда из (3) получаем уравнение

ad'oth(x, у,т)= - Ah2 - соkh + fu (5)

жгде

дх2 ду2

-оператор Лапласа.

Только при (1 —>• 1 равнение (5) переходит в известное в классиче--<зкой теории движения грунтовых вод уравнение

adt = 2Ah2 + (б)

В отличие от уравнений (4) и (6) нелинейные уравнения (3) и (5) .являются нелокальными и относятся к классу нагруженных уравнений и частных производных параболического типа.

Пусть решение h(x, у, t) уравнения (2) удовлетворяет условию Коши

h(x,y,0) = ip(x,y), (х,у)€П, (7)

где (р(х,у) - уровень грунтовой воды в начальный момент времени.

В этом же параграфе доказана эквивалентность уравнения (2) уравнению

ЩЦх, у, т) - <р(х, у)] = L{h, h) - L(hо, h) + w0 + w. (8)

Целью исследований, проводимых в параграфах 1.2-1.4, является разработка методов линеаризации существенно нелинейного уравнения (2). В этих параграфах предложены и обоснованы три способа линеаризации нагруженного обобщенного уравнения Буссинеска.

Суть первого способа линеаризации, изложенного в § 1.2, состоит в следующем. Разность К/к = h — hа представляет собой мощность водоносной части пласта. Пусть

¿W = j^j J[Мж> У> *) ~ М*> y))dxdy (9)

— среднее значение этой разности по области П с площадью |П|. В дифференциальном выражении L{h—ho\h) функцию h—ho в каждый момент

времени t заменим ее средним значением 5(Ь) и нелинейное уравнение (2) аппроксимируем нагруженным (линейным относительно h при заданном значении 8) уравнением

д&[<тН(х, У-* Т)1 = Л) - - Ло) + /о,

(10)

точное решение которого принимается за приближенное его решение.

Уравнение (10) и нелокальное для него условие (9) в случае уравнения (5) записываются в виде

доМх, у, г) = аДЛ - &(Л - Н0) + е, = 0 <*<Т,

(И) (12)

где а = кб/о называется коэффициентом уровнепроводности, Ь = с^к/сг, £ = ш/ст, Л(П;£) - среднее значение функции h(x, у, $ по области $7.

При линеаризации уравнения (4) с постоянными а, к и CQ ПО методу Буссинеска полагают (П.Я. Полубаринова-Кочина, В.Т. Пряжинская, В.Н. Эмих) функцию 8 = не зависящей от времени. В силу (7) и (12) это означает, что для уравнения (11) должно задаваться нелокальное условие типа условия А.А. Самарского для одномерных уравнений параболического типа:

Когда процесс движения грунтовых вод существенно нестационарен, задача определения коэффициента уровнепроводности а = a(t) экспериментальным или теоретическим путем становится весьма затруднительной и дорогостоящей. В этом случае в диссертации предлагается математическая модель, в основе которой лежат следующие уравнения:

где - единичный вектор внешней нормали к дО., сЬ - элемент границы дО..

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

|П|(Л-йо)а

/ИЖ>У) - ко{х,у)]<1х<1у' п

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

с?Л. , „,. . А, dh , dw

которое при ц —> 1 принимает вид

d2h tilMA. dh ,dw

(16)

(17)

Эффективность третьего способа демонстрируется на примере следующей математической модели одномерной задачи равномерного подъема уровня воды в каналах: h(x, 0) = 0, Л(0, i) =at,a = const > 0,

dh _ к d2h2 dt ~ 2a dx2'

h — h(x,t).

(18)

Когда инфильтрация w не зависит от времени, водоупор непроницаем и подъем среднего уровня грунтовых вод происходит равномерно (¿'(£) = а), из уравнений (16) и (17) соответственно получаем

где

При этом соответствующее уравнению (18) линеаризованное уравнение Буссинеска принимает такой же вид, что и одномерное волновое уравнение

&2h _

dt2 dx2'

(21)

Волновые уравнения (20) и (21) учитывают конечность скорости распространения передней границы возмущенной области движения. Линейная математическая модель (21) движения грунтовых вод сохраняет важное свойство нелинейного уравнения параболического типа (18), обнаруженное в известных работах Я.Б. Зельдовича и А.С. Компанейца, Г.И. Баренблатта и М.И. Вишика в начале 50-х годов.

Уравнение (19) относится к классу уравнений дробной диффузии (фильтрации), исследованной в основополагающей работе А.Н. Кочубея, опубликованной в 1990 г. в журнале «Дифференциальные уравнения» (т. 26, № 24).

Представительная библиография работ, посвященных уравнениям дробной диффузии, содержится в монографии A.M. Нахушева «Элементы дробного исчисления и их применение» (Нальчик: изд-во КБНЦ РАН, 2000).

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

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

Параграф 2.2 посвящен анализу математических моделей динамики грунтовых вод, в основе которых лежит уравнение Лаврентьева-Бицадзе с нулевым начальным условием. Здесь на базе теоремы Коши-Ковалевской и принципа экстремума Зарембы-Жиро обнаружены экстремальные свойства модели и ее разрешимости в малом. Обращено внимание, что нелинейность процесса проявляется через некорректность по Адамару задачи Коши для уравнения Лапласа, и эффективным выходом из этой ситуации может стать не замена уравнения Лаврентьева-Бицадзе уравнением Трикоми, а его замена уравнением вида

к,

(22)

где

Н{ц) - функция Хевисайда, Л = ~ уровень грунтовой воды в

точке £ в момент времени t.

Уравнение (22) в случае, когда fi(t) меняется по степенному закону: fi(t) = — f.Csign^. — t), где До = const > О, ¡3 = const > О, U~

критическое время, допускает следующую запись:

дШЦ~у)и Q2U

(23)

где х = ^yjo/knv, y = t-t„ и(х, у) = Цх^/Щц/а, у +1,).

Уравнение (23) входит в класс уравнений смешанного параболо-гиперболического типа, исследованного Т.Д. Джураевым, В.А. Елее-вым, A.M. Нахушевым, К.Б. Сабитовым и их учениками.

В § 2.3 исследуется модель (23) в специальной прямоугольной области на плоскости точек (х, у) при нулевых начальных распределениях, однородном условии на депрессионной кривой dx/dy = {—у)^2 и нелокальном краевом условии

где - флуктурирующая

сила.

Основным научным результатом этого параграфа является Теорема 2.3.1. Пусть и{х, 0) €С2]0, г[ПС[0, г], щ(х, 0) еС^О, г[П£[0, г], А=Д/(2/^4)^71=оГ(2А)Г(1-А)(2-4А)1-ад/[Г(А)Г(2-2Д)]. Тогда

В § 2.4 разработан класс линейных математических моделей динамики грунтовых вод с горизонтальным водоупором, учитывающих явление последействия, когда интенсивность фильтрации 1[с — Л£ — {к/а){1 и скорость расхода в слое меняются по линейным законам:

где Ai(f), Аг(0 - заданные функции из С[0,/], а с и т - заданные положительные числа.

В случае, когда / = kv(h — На)/Мо и рассматриваются задачи долгосрочного гидрогеологического прогнозирования при естественном или искусственном орошении больших площадей, правая часть (25) заменя-ечся величиной, пропорциональной — í|"'sign (t, — t), где t, —время,

когда расход грунтовой воды в слое 0 <1 достигает максимального значения, а затем падает до значения, не нарушающего экологию зоны орошения. Это условие позволяет аппроксимировать уравнение Бусси-неска уравнением смешанного типа

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

Важным результатом параграфов 2.5 и 2.6 является развитие одного способа поиска нелокального краевого условия для уравнения Буссине-ска (18), предложенного A.M. Нахушевым в 1982 году (Дифференциальные уравнения. Т. 18, № 1). В основе предлагаемого в §2.6 алгоритма поиска нелокального краевого условия для нагруженного уравнения Буссинеска в случае горизонтального водоупора лежит уравнение

д&5(т) + b(t)5{t) - e(t) = a(t)Ah, a(t) = k6{t)/a

(27)

с начальным условием

¿(0) = S0 = ~ J <p(x, y) dxdy.

(28)

Принцип экстремума и функция Грина С(х,у;£,т)) для уравнения Лапласа Аи = 0 позволяют из (27) получить нелинейное дифференциальное уравнение дробного порядка

d^(r)+B(t)8(t)-S2(t)/(G)=e(t),

(29)

где

N( - внутренняя нормаль ЗП в точке £ = £(s) + ^¡(s), s - длина кривой di2, отсчитываемая от фиксированной точки в положительном направлении.

Основной результат § 2.6 при В — const сформулирован в виде следующей теоремы.

Теорема 2.6.1. Приближенное решение Sn(t) задачи Коти (28) для уравнения (29) можно построить по итерационной схеме

где

ou j

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

При ц —> 1 уравнение (29) переходит в уравнение Риккати и его приближенное решение можно найти одним из численных методов, например, методом Рунге-Кутта. Для широких классов граничных и начальных условий средний уровень грунтовой воды меняется по логистическому закону и можно ввести понятие емкости водоносного пласта.

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

где и = и(х, t) - влажность в точке х в момент времени t\ D(u) и k(u) — коэффициенты диффузитивности и влагопроводности при влажности и.

Моделирование баланса почвенных вод на основе различных вариантов уравнения (30) был объектом исследования многих авторов: С.Ф. Аверьянов, Я. Бэр, Д.З. Заславский, С. Ирмей, СВ. Нерпин, А.Ф. Чудновский, Р. Дж. Ханке и др.

Численной, как классической, так и суммарной аппроксимации уравнений тепломассопереноса параболического типа вида (4) и (30) посвящено значительное число работ, среди которых следует отметить работы П.Н. Вабищевича, А.Я. Горбачевского, Л.А. Крукиера, A.M. Наху-шева, А.И. Сухинова и М.Х. Шханукова.

В §2.7 анонсируется, что: для широких классов почв, включающих почвы типа Гарднера, имеет место логистическое уравнение типа Ферхюльста-Пирла

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

д2и

ди

71) = аа— - 6„—, 0 < а < 1,

(32)

где а„ = А//х - коэффициент фрактальной диффузий, - фактическая скорость движения влаги в порах грунта.

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

дачи Коши

г) = тЦ), «х(о,¿) = , о<t<т

(33)

для уравнения (32). Суть метода состоит в замене уравнения (32) на каждом шаге нагруженным

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

д2и ди 1

и{х„7]) Хь — Х и(х,+1,Т)) Х,+1~Х

, » = 0,1,...,п

и эффективном использовании элементов дробного исчисления.

Отметим, что решение -Р(и) = А///, Ц Ф 0 уравнения (31) является устойчивым по Ляпунову. Если в уравнении (30) £*(ы) = Х/ц и скорость движения влаги под действием гравитационных сил постоянна, то и к условию (33) надо присоединить условие

и(х, 0) = <р(х), 0 < х < г,

(34)

задающее глубинный ход влажности в начальный момент времени. Смешанная задача (33), (34) не является корректной по Адамару для уравнения (32) при а = 1.

В § 2.8 в качестве базового уравнения математической модели одномерного движения почвенного раствора в почвогрунтах, интерпретиру-

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

где и = и(х, t) - концентрация почвенного раствора в точке х в момент времени t > 0; а - коэффициент фрактальной диффузии, 6 - скорость конвекции,а €]0,1[, ß - коэффициент растворения соли; и» - предельная концентрация насыщения.

Сходимость разностных схем для дифференциальных уравнений в частных производных с дробной производной впервые исследована М.Х. Шхануковым (ДАН, 1996, т. 348, №6).

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

Интересные исследования в этом направлении провели К.В. Oldhain, J. Spanier, А.Р. Бечелова, М.А. Керефов, A.M. Нахушев, М.З. Худалов.

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

Задача 2.8.1. Найти регулярное в любой точке х €]0,г[ и для любого момента времени t > 0 решение и = u(x,t) уравнения (35), ограниченное при f —> 0 и удовлетворяющее условию Коши: u(0, t) = r(t), u*(0,f) = *(i),0<i<T.

В данном случае - минерализация почвенного раствора, а -"поток" концентрации на поверхности почвогрунта х = 0.

Здесь следует отметить глубоко содержательные исследования К.М. Магомедова, Р.П. Мейланова, A.M. Нахушева, В.А. Нахушевой, P.P. Нигматулина, Р.Ш. Нигматулина, посвященные реализации концепции фрактала в теории фильтрации в средах со сложной геометрической структурой порового пространства.

Весьма важные исследования по применению операции дробного интегро-дифференцирования постоянного и переменного порядков проведены в работах В.Л. Кобелева, О.Л. Кобелевой, Я.Л. Кобелева и Е.П. Романова.

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

В §3.1 в области ü„ = {{х,у) : 0 < х < г, —пг < у < ß} рассматри-

вается уравнение Лаврентьева-Бицадзе д2и д2и

sisay'd^ + df = 0, 0<ж<г = //\^ (36)

со следующими нелокальными и локальными начально-краевыми условиями:

Здесь: п - целое число, у = t—t„ х = и(х,у) = h(x\/c,y+t,); предполагается, что горизонтальный водоупор непроницаем; инфильтрация на единицу площади в единицу времени пренебрежительно мала (см. уравнение (26); объект исследования -

Задача 3.1.1. Найти решение и(х,у) уравнения (36), удовлетворяющее условиям (37), (38) и (39), где д(у), kn(x), h(x), Ф0(у), Фг(у) - заданные функции, ц(у) € С1 [0,/?], ФоЫ £ С[-тп, 0], Фг(у) 6 С[-гп,0\, h(x) и hu(x) е С[0, г] .

Доказательство теорем единственности и существования решения задачи 3.1.1 существенно опирается на следующую лемму.

Лемма 3.1.1. Пусть существует решение и(х,у) задачи 3.1.1. Тогда для любой точки х € [0, г] »

и(х, 0) = ]Г{а4Фг(х - кг) + Фо( х -(к- 1)г)] -*=i

-а*+1[Фг(-х - (к - 1)г) + Ф0(х - *г)]} + Н„(х), (40)

где

_/1, А; = 1 (mod 2); . . _ Г -hn (г - х), n=l(mod2); °к ~ \ 0, к = 0 (mod 2); л"~ \ hn (х), н = 0 (mod 2).

Теорема 3.1.1 единственности решения задачи 3.1.1 доказана с помощью леммы 3.1.1 с использованием свойств гармонических функций. Лемма 3.1.1 позволяет свести вопрос о существовании решения задачи 3.1.1 в эллиптической части Í2+ смешанной области Í7„ к следующей новой нелокальной задаче для уравнения Лапласа

Д и = 0. (41)

Задача 3.1.2. Найти регулярное в области и непрерывное в замыкании i2+ решение и(х, у) уравнения (41), которое удовлетворяет условию гладкости их е С(\х\ < г, 0 < у < ¡3), краевым условиям и(г,у) = Ф,(у), их(г,у) - ит (0, у) - {у) для всех у е]0,г[ и u(x,(i) = h{x), и(х,0) = т(х) для всех х € [0,г].

Справедлива

Теорема 3.1.2. Пусть г(х) € С^О.г] Л С2]0, r[, h{х) € C'fO.r],

т'(О) — /(г) = fj,'(0). Тогда существует решение задачи 3.1.2.

Фурье решения задачи 3.1.1 по схеме, предложенной Н.И. Ионкиным и Е.И. Моисеевым для уравнений параболического типа (Дифференциальные уравнения, 1979. Т. 15, № 27).

В §3.2 для уравнения (36) исследована эталонная видоизмененная начально-краевая задача 3.2.1, которая отличается от задачи 3.1.1 тем, что условие и(г, у) = Ф,(у) в (38) заменено условием hr(r,y) = Фг(у), где Уг{у) - закон изменения уклона потока грунтовых вод в точке х = г.

Задачи 3.1.1, 3.2.2 существенно отличаются от задачи Дирихле для уравнения смешанного типа в прямоугольных областях, исследованной в работах Н.Н. Вахания, A.M. Нахушева и М.М. Хачева.

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

и начально-краевых условий: и(х,0) = т(ж); иу(х, 0) = 0 при т ф 0, w„(x, 0) = А с/г при т = 0,

Уравнение (42) - следствие уравнения (23) при у < 0. Основная цель четвертой главы - разработка линейных математических моделей эволюции малых возмущений в каналах с пористыми и проницаемыми стенками.

В §4.1 произведены выбор и анализ базовых уравнений модели рас-

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

Пусть с - скорость звука в жидкости, х - коэффициент пьезо-проводности, т - пористость, к — высота канала, а2 = 2л/Тст/к",

Главное научное достижение этого параграфа - вывод базового уравнения

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

В § 4.2 для нелокального волнового уравнения (43) доказаны теоремы единственности и существования решения следующей смешанной задачи.

Задача 4.2.1. В области Б— {(х,*) : 0<х<г, 0<<<Т} найти регулярное решение и = и(х, уравнения (43), удовлетворяющее условиям и(0,г) = 0, и(г,г) = 0 € [О, Т]; и(х,0) = <р(х), щ{х, 0) = Ф(х) Уж е]0,г[; обладающее тем свойством, что и((г, 4) 6 ¿[0, Т] для любого х > 0, 1р(х) и Ф(х) - заданные функции, причем <р(х) € С2[0, г], Ф(х) 6 С1 [0, г].

Основным результатом § 4.2 являются

Теорема 4.2.1. Задача 4.2.1 не имеет более одного решения.

Теорема 4.2.2. Задача 4.2.1 имеет решение.

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

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

Здесь ро - плотность, отнесенная к невозмущенному состоянию, х0 -характерная точка среды, где реализуемо измерение давления в заданные моменты времени. '

Гипотеза (44) приводит к существенно новому уравнению математи-

ческой модели фильтрации с условием Коши:

(45)

(46)

тогда

и(ю, *) = МЫ 0 - с2 а2 ! {г- т))-" х

-чГ'х

X Е^ [-<^(4 - т?)1-"; 1 -а] КЫ ч)<1т)- (47)

Формула (47) - основной научный результат § 4.3. Она позволяет свести задачу Коши (46) для уравнения (45) к этой же задаче для неоднородного одномерного волнового уравнения, является базовой при разработке эффективных аналитических и численных методов ее решения.

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

Объектом исследования в § 5.1 является математическая модель процесса трансформации полей температуры и влажности при стационарных условиях.

Метеорологический режим искусственно или естественно орошаемой земной поверхности По евклидова пространства Е3 точек (х, у, г) существенно отличается от режима окружающей неорошаемой территории. Влажность на поверхности По и ее температура зависят от влажности и+ и ее температуры и+ в приземном слое 0 < г < к атмосферы.

При определенных, принадлежащих Д.П. Лайхтману (Физика пограничного слоя атмосферы, 1970) предположениях, когда характеристики метеорологического режима не зависят от времени, температура и влажность могут быть определены как

решение следующей расщепленной системы дифференциальных уравне-

ний с частными производными второго порядка параболического типа:

где к+ (г) и т+ (г) - коэффициент турбулентности и скорость, соответственно, на высоте z от поверхности fi0-

Пусть v~ = v~ (z) — влажность, а и~ = и~ (г) - температура на неорошаемой территории Oq , окружающей область fí0, или на самом массиве fio до орошения. Тогда функции и (z„ х) = и+ (z„ х) — и~ (z) и v (л, х) = v+ (z, х) — v~ (г) удовлетворяют уравнению

(50)

и граничному условию

В случае, когда к < 50 — 100 м, функции тп+ {£) и к+ (г) можно по Лайхтману аппроксимировать степенными функциями и уравнение (50) переписать в виде

ТП1

/ z\' ди _ ^ д / 2 \ \ZiJ дх idz /

1-е

ди dz'

(52)

где mi и к \ — скорость ветра и турбулентность на характерной высоте Z\ от поверхности земли, е - положительное число меньше 1, S = const > 0.

Уравнение (52) относится к классу уравнений в частных производных параболического типа с вырождением порядка на линии 2 = 0.

В § 5.2 доказана пропорциональность турбулентного потока дробной производной от удельной влажности на деятельной поверхности по следующей схеме:

кгх (z\e ,_„ , , 5 — е + 1

£ =

riliz\

,r/=c(f)'

<?-»е2 = 1, р =

уравнение (52) эквивалентно уравнению Жевре

2. Пусть

N(x) = -k

f z\ du \zij dz

u(0,x) - pa

- турбулентный поток, Т(х) = и(0,х) - распределение субстанции и (z, х) на деятельной п о в „, т(£) = UI „.э

гда

3. Пусть и(г, х) - регулярное при х > О, г > 0 решение уравнения (52), удовлетворяющее условию (51) и обладающее тем свойством, что существуют Итди/дг] = Итди/д£ = т^), где 1/(£) и т(£)

суммируемы на сегменте О < £ < Для любого £0 > 0; г (0) = 0,

1ш1 и (г, х) = 0 для любой фиксированной точки х > 0- Тогда

D^r(t), а =

р+ 2

(54)

При выводе связи (54) существенно использовано фундаментальное решение уравнения (53), впервые построенное Жевре. Ее можно получить и из формулы Келлера-Лайхтмана, которая выведено ими при дополнительном и весьма жестком предположении, что

(zV"du А

lim - — = 0.

i-too \Z-iJ OZ

Пусть теперь N.

+(х) = lim ( — )

dv ^ dv

JV+(z) и T'+{x) € Ц0,х0] для любого Хо > 0; Т+(0) = 0, v - v(z, x) - удельная влажность, удовлетворяющая уравнению (52) при х > 0, Z > 0 lim v(z, х) = 0 я любой фиксированной точки х > 0. В этом

г—юо

же параграфе установлено, что

N+(x) = kseDlT+(t),

(55)

где

(ЛЛ1'24 [тЛчТ(1-д)

\яъ) и; Г (в) •

В § 5.3 впервые дано представление турбулентного потока, удельной влажности и температуры на деятельной поверхности через функцию Миттаг-Леффлера; в частности, строго доказано, что при граничных условиях:

с заданными характеристиками А,-, ^ = 1,2,..., 6 приземного слоя имеет место представление

где /и /з и А однозначно определяются через А,-, Еч (г) = Ещ (г; 1) -функция Миттаг-Леффлера.

Вывод формул вида (56) реализован методом редукции задачи к обыкновенным дифференциальным уравнениям дробного порядка.

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

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

где - постоянные положительные величины, зависящие от па-

раметров граничных условий.

Г+ (х) = А {Г (х) + Да, А3ЛГ (х) + (х) + А 5Т (х) = А6

Т(х) = /,^(Ах')-/1/А,

(56)

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

1. Экстремальные свойства линейных математических моделей процесса движения грунтовых вод, в основе которых лежит весьма важное в газовой динамике околозвуковых течений уравнение Лаврентьева-Бицадзе; утверждение о том, что нелинейность этого процесса проявляется через некорректность по Адамару задачи Коши в локальной постановке для уравнения Лапласа.

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

3. Установление непосредственной связи уравнения Лайхтмана с уравнением Жевре и дано представление турбулентного потока через функцию Миттаг-Леффлера.

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

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

Список публикаций по теме диссертации:

1. Сербина Л. И. Об одной линейной математической модели динамики грунтовых вод // Доклады Адыгской (Черкесской) международной академии наук. - 1997. Т.2, №2. - С. 6-9.

2. Сербина Л.И. Об одной краевой задаче для уравнения Лаврентьева-Бицадзе с нелокальным условием Самарского // Доклады Адыгской (Черкесской) международной академии наук. — 1998. Т.З, № 2. - С. 6-9.

3. Сербина Л.И. Об одной проблеме для линеаризованного уравнения Буссинеска с нелокальным условием Самарского. Сб. науч. трудов Международной конференции "Спектральная теория дифференциальных операторов и смежные вопросы", посвященной юбилею академика Ильина В.А. Стерлитамак-Республика Башкортостан. - 1998. - С. 5758.

4. Сербина Л.И. Математическая модель эволюции малых возмущений в полубесконечных каналах с пористыми и проницаемыми стенками // Доклады Адыгской (Черкесской) международной академии наук. - 1999. Т.4, № 1. - С. 60-63.

5. Сербина Л.И. Об одной задаче определения скорости фильтрации в средах с фрактальной геометрией // Материалы III региональной научно-технической конференции "Вузовская наука - СевероКавказскому региону". Ставрополь. — 1999. С. 3-4.

6. Сербина Л.И. Об одном волновом уравнении с дробной производной, моделирующем динамику малых возмущений. Тезисы докладов Международной конференции " Аналитические методы анализа и дифференциальных уравнений. AMADE-99". Минск. - 1999. - С. 203-205.

7. Сербина Л.И. Об одной математической модели эволюции малых возмущений в полубесконечных каналах с пористыми и проницаемыми стенками. Сб. науч. трудов III Всероссийского симпозиума "Математическое моделирование и компьютерные технологии", т. 2. Кисловодск -

1999. - С. 34-35.

8. Сербина Л.И. Качественный анализ одной математической модели турбулентной диффузии. Тезисы докладов Воронежской весенней математической школы «Понтрягинские чтения-Х», 1999 г. С. 224.

9. Сербина Л.И. Об одной математической модели динамики взаимодействия деятельной поверхности почвы с приземным слоем атмосферы // Доклады Адыгской (Черкесской) международной академии наук. -

2000. Т.5, № 1. - С. 85-90.

10. Сербина Л.И. Об одной математической модели процесса трансформации полей температуры и влажности при стационарных услови-

ях. Труды IX Международного симпозиума "Методы дискретных особенностей в задачах математической физики". Орел. - 2000. - С. 410-420.

11. Сербина Л.И. Об одной проблеме для линеаризованного уравнения Буссинеска с нелокальным условием Самарского. Отчет НИИ ПМА КБНЦ РАН по научно-исследовательским, опытно-конструкторским работам за 1995-2000 гг. по теме "Исследование структурных и качественных свойств решений локальных и нелокальных краевых задач для широких классов уравнений и систем основных типов и их приложений". № гос. регистрации 01.950004494, №рег. в ВИНИТИ 02.2 00 104695, кн.1.

- С. 54-59.

12. Сербина Л.И. Математические модели динамики микрометеорологического режима, описываемые дифференциальными операторами дробного порядка. Тезисы докладов IV Сибирского конгресса по прикладной и индустриальной математике (ИНПРИМ-2000). Новосибирск.

- 2000. 4.1. - С. 79-80.

13. Сербина Л.И. Об одной математической модели динамики атмосферных процессов. Сб. науч. трудов IV Всероссийского симпозиума "Математическое моделирование и компьютерные технологии", т. 2. Кисловодск. - 2000. - С. 34-35.

14. Сербина Л.И. Решение качественно новых краевых задач для уравнения Лаврентьева-Бицадзе. Труды международной конференции "Математическое моделирование, статистика и информатика в современном управлении экономикой". Самара. - 2001. - С. 193-194.

15. Сербина Л.И. Об одной математической модели динамики солей в почвах и алгоритме ее компьютерной реализации // Доклады Адыгской (Черкесской) международной академии наук. - 2001. Т.5, № 2. - С. 72-76.

16. Сербина Л.И. Об одной качественно новой краевой задаче для уравнения Лаврентьева-Бицадзе // Труды Института математики НАН Беларуси. - 2001. Т.10. - С. 142-146.

17. Сербина Л.И. Об одной математической модели эволюции малых возмущений в полубесконечных каналах с пористыми и проницаемыми стенками. Труды международной конференции " Математическое моделирование, статистика и информатика в современном управлении экономикой". Самара. - 2001. - С. 194-195.

18. Сербина Л.И. О разрешимости одного волнового уравнения с дробной производной. Тезисы докладов Международной конференции "Новые подходы к решению дифференциальных уравнений". Дрогобы-чи - Украина. - 2001. - С. 134.

19. Сербина Л.И. Алгоритм релизации математической модели соле-

вого режима в почвогрунтах. Тезисы докладов II Международной конференции "Нелокальные краевые задачи и родственные проблемы математической биологии, информатики и физики". Нальчик. — 2001. -С. 37-38.

20. Сербина Л.И. Об одной краевой задаче для математической модели переноса во фрактальной среде. Тезисы докладов Воронежской весенней математической школы "Понтрягинские чтения"-ХП "Современные методы в теории краевых задач". Воронеж. - 2001. - С. 141 - 142.

21. Сербина Л.И. Об одной математической модели безнапорного движения грунтовых вод // Известия КБНЦ РАН. 2002. № 1(8). С. 94-102.

22. Сербина Л.И. Об одной модели динамики микрометеорологического режима при орошении // Математическое моделирование. - 2002. Т.14, №4.-С. 67-79.

23. Сербина Л.И. Об одной проблеме для линеаризованного уравнения Буссинеска с нелокальным условием Самарского // Дифференц. уравнения. - 2002. Т. 38, № 8. - С. 1113-1119.

24. Сербина Л.И. Математическая модель расчета влагосолепереноса во фрактальных средах. Тезисы докладов X Международной конференции «Математика. Экономика. Образование». Ростов-на-Дону, 2002. С. 57-58.

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

26. Сербина Л.И. Анализ математической модели динамики грунтовых вод // Известия Вузов. Северо-Кавказский регион. Естественные науки. - 2002. № 2. - С. 19-21.

27. Сербина Л.И. Об одном волновом уравнении с дробной производной, моделирующем динамику малых возмущений. Тезисы докладов IV международной конференции по неравновесным процессам в соплах и струях. Санкт-Петербург. - 2002. - С. 389-391.

28. Сербина Л.И. О нелокальной начально-краевой задаче для модельного уравнения Лаврентьева-Бицадзе. Материалы Международного Российско-Узбекского симпозиума "Уравнения смешанного типа и родственные проблемы анализа и информатики". Нальчик-Эльбрус. 2003. С. 149-151.

29. Сербина Л.И. Об одной математической модели переноса субстанции во фрактальных средах // Математическое моделирование. 2003. Т. 15, №9. С. 17-28.

30. Сербина Л.И. Математическая модель динамики почвенной вла-

ги во фрактальных средах // Вестник Дагестанского научного центра. 2003. № 14. С. 16-21.

31. Сербина Л.И. Нелокальная начально-краевая задача для модельного уравнения нестационарной фильтрации//Доклады Адыгской (Черкесской) международной академии наук. - 2003. Т. 6, № 2. - С. 125-130.

32. Сербина Л.И. Решение одной начально-краевой задачи теории фильтрации с нелокальным краевым условием // Вестник СамГТУ. Серия: Физико-математические науки. - 2003. Вып. 19. - С. 16-21.

33. Сербина Л.И. Об одной математической модели влагопереноса в пористых средах с фрактальной структурой. Тезисы докладов Международной конференции " Аналитические методы анализа и дифференциальных уравнений. АМАББ-2003". Минск. - 2003. - С. 160-161.

Изд. лиц. серия ИД № 00502 Подписано к печати 10.03.04 г.

Формат 60x84. 1/16 Усл. печ. л. - 1,9. Уч.-изд. л. - 1,5.

Бумага офсетная. Печать офсетная. Заказ 920 Тираж 100 экз. Северо-Кавказский государственный технический университет 355029 г. Ставрополь пр. Кулакова, 2

Отпечатано в типографии СевКавГТУ Издательство Северо-кавказского государственного технического университета

»-67 10

Оглавление автор диссертации — доктора физико-математических наук Сербина, Людмила Ивановна

Введение

1. Основные уравнения моделей движения грунтовых вод

§1.1. Обобщенное уравнение Буссинеска с дробной производной по времени

§ 1.2. Первый способ линеаризации нагруженного уравнения Буссинеска

§ 1.3. Второй способ линеаризации нагруженного уравнения Буссинеска

§ 1.4. Третий способ линеаризации нагруженного уравнения Буссинеска

2. Линейные одномерные математические модели движения грунтовых вод и почвенной влаги

§2.1. Анализ математической модели одномерного движения грунтовых вод, основанной на волновом уравнении.

§ 2.2. Анализ математической модели динамики грунтовых вод, основанной на уравнении Лаврентьева-Бицадзе с нулевым начальным условием.

§ 2.3.Математическая модель движения грунтовых вод, основанная на уравнении смешанного параболо- гиперболического типа с нулевым начальным условием

§2.4.06 одном классе математических моделей динамики грунтовых вод с горизонтальным водоупором

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

§ 2.6 Алгоритм поиска нелокального краевого условия для нагруженного уравнения Буссинеска в случае горизонтального водоупора

§ 2.7.06 одной математической модели движения почвенной влаги и алгоритме ее компьютерной реализации.

§ 2.8 Алгоритм реализации математической модели движения почвенного раствора

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

§ 3.1.Эталонная начально-краевая задача для смешанного типа уравнения одномерного движения грунтовых вод с горизонтальным водоупором

§ 3.2. Видоизмененная эталонная начально-краевая задача для уравнения

Лаврентьева-Бицадзе.

§ 3.3.Об алгоритме долгосрочного прогноза динамики грунтовых вод.

4. Математическая модель эволюции малых возмущений в каналах с пористыми и проницаемыми стенками

§ 4.1.Выбор и анализ базовых уравнений.

§ 4.2.Смешанная задача для нелокального волнового уравнения с оператором дробного дифференцирования в младшем члене.

§ 4.3.Модификация уравнения модели фильтрации, учитывающая явления последействия .•.

5. Математическая модель динамики микрометеорологического режима при орошении

§ 5.1.Математическая модель процесса трансформации полей температуры и влажности при стационарных условиях.

§ 5.2. Пропорциональность турбулентного потока дробной производной от удельной влажности на деятельной поверхности.

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

§ 5.4.Математические модели водопотребления и нормы орошения.

§ 5.5.Качественный и сравнительный анализ математической модели динамики микрометеорологического режима при орошении и формул Лайхтмана

Введение 2004 год, диссертация по информатике, вычислительной технике и управлению, Сербина, Людмила Ивановна

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

П. Я. Полубаринова-Кочина [1]; С.Ф. Аверьянов [2]; М.Г. Андерсен, Т.П. Берт [3], Я. Бэр, Д. Заславски, С. Ирмей [3]; JI. Дакштейн (см. [3. С. 681]); Р.Дж, Ханке (см. [3. С. 27]) и другие в своих исследованиях обращали внимание на важность выяснения физико-математических основ процессов переноса в природных водоносных системах и создания их математической модели для решения задач гидрогеологического прогнозирования, защиты подземных вод от загрязнения и борьбы с засолением орошаемых земель.

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

Исследование по теме докторской диссертации проводилось в рамках плана научно-исследовательских работ Института прикладной математики и автоматизации КБНЦ РАН по следующим двум направлениям фундаментальных исследований:

- «Нелокальные дифференциальные операторы основных и смешанных типов и их применение к фундаментальным исследованиям в различных областях знаний» (№ГР 01.20.00 12841);

- «Развитие дробного исчисления и анализа на фракталах для разработки математических моделей физико-биологических процессов и сред с фрактальной структурой» (№ГР 0'1.20.00 12845); а также проекта № 00-01-00311 - «Исследование класса задаваемых дифференциальными операторами дробного порядка математических моделей тепломассопереноса в средах с фрактальной структурой», поддержанного в 2000, 2001, 2002 годах Российским фондом фундаментальных исследований.

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

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

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

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

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

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

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

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

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

Основные результаты диссертации изложены в работах автора [56-87], в том числе в монографии [74] объемом 144 страницы.

Первая глава посвящена выводу основных уравнений моделей неустановившегося движения грунтовых вод в безнапорном пласте со сла-боизменяющейся поверхностью водоупора z = h0(x,y) и свободной поверхностью z = h(x,y,t) в области Q евклидовой плоскости точек (х, у) с абсциссой х и ординатой у в момент времени t от начального t = 0 до расчетного t = Т, а также разработке и анализу методов их линеаризации.

В этой главе, как и во всей диссертации, существенно используются элементы дробного исчисления и в первую очередь операторы и 8qV Оператор D%t есть оператор дробного интегро-дифференцирования порядка |а| в смысле Римана-Лиувилля:

D%tu{x, у,т)= < t

J f u(x,y,r)dy

-a) J (t- ' a < 0,

Г (-a) J (t — r)c о u(x,y,t), a = 0, d[al+1 a-fal-l

-Dqi и(х,у,т), a> 0; где F(x) - гамма-функция Эйлера, [a] - целая часть числа а, которая удовлетворяет неравенству [а] < а < [а]+1. Если п—1 < а < п = 1,2,., то по определению величина и она называется регуляризованной производной по t порядка а или производной по М. Caputo [5].

В параграфе 1.1 на основе интерпретации почвогрунта как сильно пористой среды с распределенными параметрами, структура пор которой образует множество с фрактальной размерностью cf/, и гипотезы о том, что за промежуток времени At происходит изменение уровня грунтовых вод на 5h = htAt, где

Ы = д&к(х,у,т), (1) получено нагруженное обобщенное уравнение Буссинеска d%t[ah{x, у, г)] = L(h - Л0; h) + w0 + w. (2)

Здесь fi - число из интервала ]0,1[; а - водоотдача или недостаток насыщения грунта; dh\ ду)

- дифференциальное выражение;

- скорость просачивания через водоупор, к - коэффициент фильрации; к0 - коэффициент фильтрации слабопроницаемого водоупора, который значительно меньше h, HQ - напор в нижележащем водоносном пласте; Mq = Мо(х,у) - мощность водоупора; w = w(x,y,t) - величина инфильтрации (или разность между инфильтрацией и испарением) на единицу площади в единицу времени.

Для широких классов пористых сред можно приближенно положить li = df.

Уравнение (2) эквивалентно уравнению dt[ah(x, у, т)\ = 1; h2) - L(hQ; h) - c^kh + Д (3)

С Co = уfi = w + kcoHQ.

KMq

Поскольку limд&и(х, y, t) = /1-Я ut то уравнение (2) при ц —> 1 совпадает с обобщенным (по терминологии

П.Я. Полубариновой-Кочиной) уравнением Буссинеска dah д (ТМ\ д (ТМ\ „ , . где К = k{h — hQ) - проводимость пласта, /0 = cok(HQ — ho).

Пусть а, к, М0 - постоянные величины, водоупор - горизонтален и ho = 0. Тогда из (3) получаем уравнение к adSth(x, у, т) = -Ah2 - cokh 4- /ь (5) д2 д2 где Д = ——- + -——■ - оператор Лапласа. ох2 оу1

Только при fx 1 уравнение (5) переходит в известное в классической теории движения грунтовых вод уравнение а^ = ^Ah2 - cokh + fx. (6)

В отличие от уравнений (4) и (б) нелинейные уравнения (3) и (5) являются нелокальными и относятся к классу нагруженных уравнений в частных производных параболического типа.

Пусть решение h(x,y,t) уравнения (2) удовлетворяет условию Коши h(x,y,0) = (р(х,у), (х,у)еП, (7) где (р(х,у) - уровень грунтовой воды в начальный момент времени.

В этом же параграфе доказана эквивалентность уравнения (2) уравнению

Dj£t[h(x,y,T) - (р(х,у)] = L(h,h) - L(hQ,h) + wQ + ги. (8)

Суть первого способа линеаризации, изложенного в § 1.2, состоит в следующем. Разность = h - h0 к представляет собой мощность водоносной части пласта. Пусть = Щ / у' ~ y^dxdy ® п

- среднее значение этой разности по области Q с площадью ]0| . В дифференциальном выражении L(h — ho',h) функцию h — hQ в каждый момент времени t заменим ее средним значением £(£) и нелинейное уравнение (2) аппроксимируем нагруженным (линейным относительно h при заданном значении 8) уравнением d&[<rh(x, у, г)] = Ц6; h) - cok{h - h0) + /о, (10) точное решение которого принимается за приближенное его решение.

Уравнение (10) и нелокальное для него условие (9) в случае уравнения (5) записываются в виде

• 35Л(®, у, т) = aAh - b(h - tf0) + е, (11) h(Q;t) = S(t), 0<t<T, (12) kS с0к где а = — называется коэффициентом уровнепроводности, b = -, а а fjj — е = —, h(Q]t) - среднее значение функции h{x,y,t) по области Q. о

При линеаризации уравнения (4) с постоянными сг, к и Cq по методу Буссинеска полагают (П.Я. Полубаринова-Кочина, В.Т. Пряжинская, В.Н. Эмих [6]) функцию 5 = S(t) не зависящей от времени. В силу (7) и (12) это означает, что для уравнения (11) должно задаваться нелокальное условие типа условия А.А. Самарского для одномерных уравнений параболического типа: h(Q; t) = Tp{Q,\0).

ЗДт) = [N(t) - b(t)]5(t) + -МП; t), (14), т дП где it - единичный вектор внешней нормали к dQ, ds - элемент границы дО,.

-„,() =шк*-*»)'. f[<p(x, у) - hQ(x, у)] dxdy п

2) нагруженными уравнениями гиперболического и смешанного типов. Примерами таких уравнений является, уравнение tdh , . . А , dh , dw . aD%t— = kS (t)Ah - о,— + к—, (16) которое при ц -> 1 принимает вид d2h . г//«\ л 1 dh кты-сй-ш+к-. (17)

Эффективность третьего способа демонстрируется на примере следующей математической модели одномерной задачи равномерного подъема уровня воды в каналах: dh к d2h2 h{x, 0) = 0, h{О, t) = at, а = const > 0. (19)

Когда инфильтрация w не зависит от времени, водоупор. непроницаем и подъем .среднего уровня грунтовых вод происходит равномерно ($'(£) = а), из уравнений (16) и (17) соответственно получаем irjj^M)=c^Aht (20) ' • (21) otJc где с1+/1 = —. При этом соответствующее уравнению (18) линеаризованное уравнение Буссинеска принимает такой же вид, что и одномерное волновое уравнение

2d2h

W =са? (22)

Волновые уравнения (21) и (22) учитывают конечность скорости распространения передней границы возмущенной области движения. Линейная математическая модель (22) движения грунтовых вод сохраняет важное свойство нелинейного уравнения параболического типа (18), обнаруженное в известных работах Я.Б. Зельдовича и А.С. Компанейца, Г.И. Баренблатта и М.И. Вишика в начале 50-х годов.

Уравнение (20) относится к классу уравнений дробной диффузии (фильтрации), исследованной в основополагающей работе А.Н. Кочубея

7].

Представительная библиография работ, посвященных уравнениям дробной диффузии, содержится в монографии A.M. Нахушева [8].

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

Н (д) - функция Хевисайда.

Уравнение (23) в случае, когда j.i(t) меняется по степенному закону: fi(t) = fio\t - t^sign (£« - t), где fiQ = const > 0, /3 = const > 0, t* - критическое время, допускает следующую запись:

Ql+H(-y) Q2 \уф, 0 <х <г, (24, где х y = t-U, и(х, y) = h У + t*

Уравнение (24) входит в класс уравнений смешанного парабологиперболического типа, исследованного Т.Д. Джураевым [9], В.А. Елеевым [10], A.M. Нахушевым [11], К.Б. Сабитовым [12] и их учениками.

В § 2.3 исследуется модель (24) в специальной прямоугольной области на плоскости точек (х, у) при нулевых начальных распределениях, одноdx родном условии на депрессионной кривой —— = (—yf!2 и нелокальном краевом условии dy г д Г J u(x,y)dx = <Рп{у), -U<y<T„ (25) где

I 771

Vn{y) = [ld0\yfsigny + f(y + Q] J—, f(t) - флуктурирующая сила.

Основным научным результатом этого параграфа является Теорема 2.3.1. Пусть и(х, 0) 6 С2]0, г[ПС[0, г], иу(х, 0) <ЕС1]0, г[ПЬ[0, г], А =^/(2/3 + 4), Ti = Г(2/51)Г(1-А)(2-4Д1)1-2^/[Г(А)Г(2-2А)]. Тогда иу(х) = 7lD^u(t,0).

В §2.4 разработан класс линейных математических моделей динамики грунтовых вод с горизонтальным водоупором, учитывающих явление последействия, когда интенсивность фильтрации a ~ dt ад£ \ д£) и скорость расхода в слое 0 < £ < I меняются по линейным законам:

- = АхЮЛК, 0 + Аа(0 D£h(Z,T), 7 I У л К, t) d£ = сГ, (26) о где Ai(£), А2(£) - заданные функции из С[0,1], а с и т - заданные положительные числа. В случае, когда ko(h-HQ) Мп и рассматриваются задачи долгосрочного гидрогеологического прогнозирования при естественном и искусственном орошении больших площадей, правая часть (26) заменяется величиной, пропорциональной £|msign (£, — t), где t* - время, когда расход грунтовой воды в слое О < £ < I достигает максимального значения, а затем падает до значения, не нарушающего экологию зоны орошения. Это условие позволяет аппроксимировать уравнение Буссинеска уравнением смешанного типа d2h , ,т . , х d2h кп dh 1 dw .

- = ф. - trswlt. - t)-- — - + --. (2-0

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

Важным результатом параграфов 2.5 и 2.6 является развитие одного способа поиска нелокального краевого условия для уравнения Буссинеска (18), предложенного А. М. Нахушевым в 1982 году [13]. В основе предлагаемого в § 2.6 алгоритма поиска нелокального краевого условия для нагруженного уравнения Буссинеска в случае горизонтального водоугюра лежит уравнение адт) + b(t)S(t) - e(t) = a(t)Ah, a{t) = ^ (28) 7 с начальным условием

S(0) = 5° = щ/ у) dxdy• (29) n

Принцип экстремума и функция Грина G(x, у, rj) для уравнения Лапласа А и = 0 позволяют из (28) получить нелинейное дифференциальное уравнение дробного порядка

В&(Т) + В№*)-Щ = «*),■ (30) где

B(t) = b(t) + (<?>, (С) = щ1 dx dy J G(x, у; e, T))dZdr],

Q Q

W = dxdyJ -Щ

П dSl

Nq - внутренняя нормаль dQ в точке £ = €{s)+irj(s), s - длина кривой отсчитываемая от фиксированной точки в положительном направлении.

Основной результат § 2.6 при В = const сформулирован в виде следующей теоремы.

Теорема 2.6.1. Приближенное решение 5n(t) задачи Коши (29) д.ия уравнения (30) можно построить по итерационной схеме

Fn-i(t) = щПйХ-Лт) + ЩМт) + ОД, t

Sn{t) = Fn! - В J(t- ry-'Ey^-Bit - г)"; ^]Fn.x{r) dr, о где

OO j

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

При fi —у 1 уравнение (30) переходит в уравнение Риккати и его приближенное решение можно найти одним из численных методов, например, методом Рунге-Кутта. Для широких классов граничных и начальных условий средний уровень грунтовой воды меняется по логистическому закону и можно ввести понятие емкости водоносного пласта.

Параграф 2.7 посвящен линейной математической модели движения почвенной влаги и алгоритму ее компьютерной реализации. Принцип локальности в случае одномерного движения влаги в почвенном слое 0 < х < г приводит к следующему локальному нелинейному уравнению параболического типа: ди д dt дх

D(u- k{u)

31) где и = u{x,t) - влажность в точке х в момент времени t; D(u) и k(u) -коэффициенты диффузитивности и влагопроводности при влажности и.

Моделирование баланса почвенных вод на основе различных вариантов уравнения (31) был объектом исследования многих авторов: С.Ф. Аверьянов [2]; Я. Бэр, Д. Заславски, С. Ирмей [4], С.В. Нерпин [14], А.Ф. Чуд-новский [15], Р. Дж. Ханке (см. [3. С. 27]) и др.

Численной, как классической, так и суммарной аппроксимации уравнений тепломассопереноса параболического типа вида (4) и (30) посвящено значительное число работ, среди которых следует отметить работы П.Н. Вабищевича [16]; П.Н. Вабищевича, А.Я. Горбачевского [17]; Л.А. Кру-киера, И.В. Шевченко [18]; A.M. Нахушева, М.Х. Шханукова [19] и А.И. Сухинова [20], [21].

В §2.7 анонсируется, что: для широких классов почв, включающих почвы типа Гарднера, имеет место логистическое уравнение типа Ферхюльста-Пирла [А - fiD(u)]D(u), \ф = const > 0; (32) du процесс движения влаги происходит в субдиффузионном режиме, имеет фрактальную во времени природу, подобен модели макроскопического движения при стохастическом переносе и с определенной точностью может быть описан уравнением д2и ди

ДМ^) = аа—-Ьа—, 0 < а < 1, (33) А где аа =--коэффициент фрактальной диффузии, оа - фактическая скорость движения влаги в порах грунта.

Главный результат § 2.7 - алгоритм высокой точности, реализованный в виде численно-аналитического метода, позволяющий найти в классе функций, ограниченных при t 0, приближенное решение задачи Коши u(0, t) = r{t), их( 0, t) = v(t), 0 <t<T (34) для уравнения (33). Суть метода состоит в замене уравнения (33) на каждом шаге hi = Xi+1 — яг», 0 < ятц '< Xi < . < хп = г нагруженным уравнением фрактальной диффузии u(Xi,7]) Xi~X u(xi+urj) xi+i-x и эффективном использовании элементов дробного исчисления. Отметим, что решение = -, А4 уравнения (32) является устойчивым по Ляпунову. Если в уравнении (31) D(u) = А/д и скорость ъ dK{u) ' du д2и . ди „, 1 , « = 0,1,.,п движения влаги под действием гравитационных сил постоянна, то и к условию (34) надо присоединить условие и(х, 0) = <р(х), 0 < х < г, (35) задающее глубинный ход влажности в начальный момент времени. Смешанная задача (34), (35) не является корректной по Адамару для уравнения (33) при а = 1.

В § 2.8 в качестве базового уравнения математической модели одномерного движения почвенного раствора в почвогрунтнх, интерпретируемых как среды с фрактальной структурой, предлагается уравнение д2и ди

Dgu{x, г,) = а— - Ъ— + /3(и. - и), 0<х< г, (36) где и = и{х, t) - концентрация почвенного раствора в точке х в момент времени t > 0; а - коэффициент фрактальной диффузии, 6 - скорость конвекции, а (Е]0,1[, f3 - коэффициент растворения соли; и* - предельная: концентрация насыщения.

Сходимость разностных схем для дифференциальных уравнений в частных производных с дробной производной впервые исследована М.Х., Шхануковым [21].

Интересные исследования в этом направлении провели К.В. Oldham, J. Spanier [23]; А.Р. Бечелова [24]; М.А. Керефов [25]; A.M. Нахушев [26]; М.З. Худалов [27].

Задача 2.8.1. Найти регулярное в любой точке х €]0, г[ и для любого момента времени t > 0 решение и = и(х, t) уравнения (36), ограниченное при t —» 0 и удовлетворяющее условию Коши: u(0,t) = r(t), ^^ = Ф(0, 0 <t<T. x-Q

В данном случае r(t) - минерализация почвенного раствора, a 4f(t) -"поток" концентрации на поверхности почвогрунта х = 0.

Здесь следует отметить весьма важные исследования К.М. Магомедо-ва [28], Р.П. Мейланова [29], [30], A.M. Нахушева [31], В.А. Нахушевой [32], P.P. Нигматулина [33], Р.Ш. Нигматулина [34], посвященные реализации концепции фрактала в теории фильтрации в средах со сложной геометрической структурой порового пространства.

Глубоко содержательные исследования по применению операции дробного интегро-дифференцирования постоянного и переменного порядков проведены в работах B.JI. Кобелева, Я.Л. Кобелевой, Л.Я. Кобелева и Е.П. Романова [35].

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

В §3.1 в области Г2П = {{х,у) : 0 < х < г, —пт < у < /3} рассматривается уравнение Лаврентьева-Бицадзе д2и(х,у) д2и(х,у) I

•-^- + ^- = 0, 0 < х < г = (37) со следующими нелокальными и локальными начально-краевыми условиями: г

Jи{х, y)dx = ii(y), 0 < у < {3 = Т — (38) о

ЦО,у) = ФоЫ, и{г,у) = Фг{у), -и<у<р\ (39) и(х,/3) = h(x), и(х, -U) = hn(x), 0 < х < г. (40)

Здесь: п - целое число, у = t — х = и(х,у) = h(xy/c,y + £»);. vе предполагается, что горизонтальный водоупор непроницаем; инфильтрация на единицу площади в единицу времени пренебрежительно мала (см. уравнение (26); объект исследования

Задача 3.1.1. Найти решение и(х,у) уравнения (37), удовлетворяющее условиям (38), (39) и (40), где ц(у), hn(x), h(x), Ф0(у), Фг(у) ~ заданные функции, fi(y) е С1[0,Р], Ф0{у) £ С[—гп, 0], Фr(y) <Е С[-гп,/3], h{x) и hn(x) <Е C[0,r] .

Лемма 3.1.1. Пусть существует решение и(х, у) задачи 3.1.1. Тогда для любой точки х £ [0, г] п и(ж,0) = ^2{ак[Фг{х - кг) + Ф0(-яг -{к- 1)г)] -к=1

-ам[Фг(-х -(к- 1)г) + Ф0(х - Ат)]} + Я„(я), (41) где и ( \ - / ~hn (г ~ ' 71 ~ 1 (mod 2)5 ПпКХ)- \ hn{x), п = 0 (mod 2).

Г 1, к = 1 (mod 2); 0, к = 0(mod 2);

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

Аи = 0. (42)

Задача 3.1.2. Найти регулярное в области Q+ и непрерывное в замыкании Г2+ решение и(х, у) уравнения (42), которое удовлетворяет условию гладкости jIL ec(|x|<r, о<у</3), дх краевым условиям: и(г,у) = Ф,М, . = *(»), V у е]0,г[; х=г дх х=0 u(x,(3) = h(x), и(х,0) = т(х), V хв [0,7-].

Справедлива

Теорема 3.1.2. Пусть т(х) £ С1[0,г]Г\С2]0,г[, h(x) £ С^О.г], ц(у) £ с2[0,(3], Фг(у) € СЧО,0\ и т(г) = Фг(0), Л(г) = Фг(/?), т'(0) -т'(г) = рЩ.

Тогда существует решение задачи 3.1.2.

В этом параграфе разработан простой, но вместе с тем эффективный метод редукции нелокальной задачи 3.1.2 к двум локальным краевым задачам для гармонических функций, позволяющий построить конструктивную блок-схему поиска решения задач 3.1.2 и 3.1.1; для модельного случая, когда г = (3 = 1, Фг(?/) = 0, h{x) = 0 реализован и метод Фурье решения задачи 3.1.1 по схеме, предложенной Н.И. Ионкиным и Е.И. Моисеевым для уравнений параболического типа [36]. x=r %(v),

В §3.2 для уравнения (37) исследована эталонная видоизмененная начально-краевая задача 3.2.1, которая отличается от задачи 3.1.1 тем, что условие и(г,у) = Фг(у) в (39) заменено условием dh(x,y) дх где Фг(у) - закон изменения уклона потока грунтовых вод в точке х — т.

Задачи 3.1.1, 3.2.2 существенно отличаются от задачи Дирихле для уравнения смешанного типа в прямоугольных областях, исследованной в работах Н.Н. Вахания [37], A.M. Нахушева [38] и М.М. Хачева [39].

Параграф 3.3 посвящен разработке алгоритма долгосрочного прогноза динамики грунтовых вод на базе уравнения Геллерстедта и начально-краевых условий: u(a;,0) = т{х)\

Ас диI • л , ди =0 при тф 0, — оу\у=о ду при т = 0; у=о г 7 J u{x, y)dx = \сут, ^ ^ = ф(у), Л = о

Уравнение (43) - следствие уравнения (24) при у < 0.

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

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

Пусть с - скорость звука в жидкости, х - коэффициент пьезопровод-ности, т - пористость, h - высота канала, о2 =

2 2 \/>ст h

2Г>3/2 , 2 д2

3/2 = azDQ't + — - cz-dt2 дх2

Главное научное достижение этого параграфа - вывод базового уравнения

Lzl2[u+{x,t)} = 0

44) математической модели распространения волн давления р — u+{x,t) в каналах с плоскопараллельными стенками, заполненных жидкостью или газом.

В § 4.2 для нелокального волнового уравнения (44) доказаны теоремы единственности и существования решения следующей смешанной задачи.

Задача 4.2.1. В области D = {(x,t) : 0 < х < г, 0 < t < Т} найти регулярное решение и = u(x,t) уравнения (44), удовлетворяющее условиям: u{0,t) = 0, u(r,t) = 0, VtG[0,T]; и(х,0) = ф), = Vse]0,r[; обладающее тем свойством, что е Цо ,т] для любого х > 0, ср(х) и - заданные функции, причем ср(х) £ с2[о,г], Ффесчо,»-]:

Основным результатом § 4.2 являются:

Теорема 4.2.1. Задача 4.2.1 не имеет более одного решения. Теорема 4.2.2. Задача 4.2.1 имеет решение.

В § 4.3 рассмотрен случай, когда скорость фильтрации v+ через стенки канала не зависит от точки х и меняется по следующему закону: h v+{x, t) = А2-—D%tu+(x0, t), 0 < xQ < г. (45)

2po

Здесь po - плотность, отнесенная к невозмущенному состоянию, Xq - характерная точка среды, где реализуемо измерение давления u+(xo,t) в заданные моменты времени.

Гипотеза (45) приводит к существенно новому уравнению математической модели фильтрации с условием Коши: ~ c2£ + (сА)2 t) = 0, (46) u(x,0) = 0, ut(x,0) = i/(x), (47) где v(x) - заданная функция из класса С1 [0,г].

Пусть x+ct

K(x,t) = ± J i/KK, x—ct тогда t u(®0, 0 = t) - с2л2 j(t- 7])~a X x [—c2A2(£ - t7)1e; 1-a] (48)

Формула (48) - основной научный результат § 4.3. Она позволяет свести задачу Коши (47) для уравнения (46) к этой же задаче для неоднородного одномерного волнового уравнения, является базовой при разработке эффективных аналитических и численных методов ее решения.

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

Объектом исследования в §5.1 является математическая модель процесса трансформации полей температуры и влажности при стационарных условиях.

Метеорологический режим искусственно или естественно орошаемой земной поверхности евклидова пространства R3 точек (х, у,z) существенно отличается от режима окружающей неорошаемой территории. Влажность на поверхности Г20 и ее температура зависят от влажности v+ и ее температуры и+ в приземном слое 0 < z < h атмосферы.

При определенных, принадлежащих Д.П. Лайхтману [40] предположениях, когда характеристики метеорологического режима не зависят от времени, температура и+ = u+(z,x) и влажность v+ = v+(z,x) могут быть определены как решение следующей расщепленной системы дифференциальных уравнений с частными производными второго порядка параболического типа: где к+ (z) и т+ (г) - коэффициент турбулентности и скорость, соответственно, на высоте z от поверхности П0

Пусть v~ = v~ {z) - влажность, a u~ = u~ (z) - температура на неорошаемой территории fig > окружающей область Q0, или на самом массиве fio до орошения. Тогда функции и (z, х) = и+ (z, х) — и~ (z) и v (z, х) = v+ (z, х) — v~ (z) удовлетворяют уравнению

51) и граничному условию ф,0) = 0. (52)

В случае, когда h < 50 — 100 м, функции т+ (z) и к+ (z) можно по Лайхтману аппроксимировать степенными функциями и уравнение (51) переписать в виде mi z \ ди 1 д f zV с ди хZi) дх Kldz\zi) dz' где mi и к\ - скорость ветра и турбулентность на характерной высоте zx от поверхности земли, е - положительное число меньше 1, <5 = const > 0.

Уравнение (53) относится к классу уравнений в частных производных параболического типа с вырождением порядка на линии z = 0.

В §5.2 доказана пропорциональность турбулентного потока дробной производной от удельной влажности на деятельной поверхности по следующей схеме: 1. Пусть д

1-е

53) = kix

2 ' rrilZl z2~pe2 = 1, р =

6-e + l

Тогда уравнение (53) эквивалентно уравнению Жевре дU d2U ^ di ~~ drf '

54) где

U = и&г,) = и L

1/е rrh4 с 1 1 £ к,

N(x) = -кх

2. Пусть z \ 1-е ди ZiJ dz

- турбулентный поток, Т (х) = и(0, х) - распределение субстанции и (z, х) на деятельной поверхности, ди

2=0 дг) г,=0 г(0 - и т]=0'

Тогда ко

3. Пусть u(z,x) - регулярное при х > 0, z > 0 решение уравнения (53), удовлетворяющее условию (52) и обладающее тем свойством, что существуют dU . dU lim— = z/(£), lim— = r 0, v^o дт] rt-*о где и (£) и т (£) суммируемы на сегменте 0 < £ < £0 для любого £0 > 0; т (0) = 0, lim и (z, х) = 0 для любой фиксированной точки х > 0. Тогда

Z—¥00 r(l-tt)

Г(а) + 2'

При выводе связи (55) существенно использовано фундаментальное решение уравнения (54), впервые построенное Жевре. Ее можно получить и из формулы Келлера-Лайхтмана, которая выведено ими при дополнительном и весьма жестком предположении, что

1-е

КО = « = ГТ-^- (55)

Пусть теперь

N+(x) = lim (-)1 £ Т+{х) = lim ^

-*о дх' iV+(x) и ТЦя) € Ь[0,я0] для любого х0 > 0; Т+(0) = 0, v = удельная влажность, удовлетворяющая уравнению (53) при х > 0, z > 0 и lim v(z,x) = 0 для любой фиксированной точки х > 0. В этом же

Z—УОО

параграфе установлено, что

N+(x) = kseDlT+ (f), (56) где

1-2 q к. = ki

4 5 + е +1' е \qzj \кг; Г (9) ' В §5.3 впервые дано представление турбулентного потока, удельной влажности и температуры на деятельной поверхности через функцию Миттаг-Леффлера; в частности, строго доказано, что при граничных условиях:

Т+ (х) = АгТ {х) + Аа, A3N- (®) + A4TV+ (я?) + А5Т (®) = Аб с заданными характеристиками А^, j = 1,2,., 6 приземного слоя имеет место представление где /ь /з и А однозначно определяются через Xj, Eq(z) = Ei/q (z; 1) -функция Миттаг-Деффлера.

Вывод формул вида (57) реализован методом редукции задачи к обыкновенным дифференциальным уравнениям дробного порядка.

В последнем параграфе 5.5 главы 5 проведен качественный и сравнительный анализ математической модели динамики микрометеорологического режима при орошении и формул Лайхтмана, а также соответствующих им интегро-дифференциальных уравнений. Найдены формулы, уточняющие формулы для вычисления функции Xq(x) — v(0,x). Показано, что она удовлетворяет дифференциальному уравнению дробного порядка и представима в виде где Xi и Х2 ~ постоянные положительные величины, зависящие от параметров граничных условий.

Заключение содержит выводы из проведенного научного исследования и вытекающие из него рекомендации, изложение в компактной форме сущности основных научных достижений автора, а также вспомогательные научные результаты, которые могут представлять самостоятельный интерес и в других областях знаний, например, в математической биологии, газовой динамике смешанных течений, при разработке математического и программного обеспечения САПР оросительных систем и систем магистральных водопроводов для сельскохозяйственного водообес-печения. В нем отмечается, что основные положения и выводы диссертации прошли апробацию на 24 научных мероприятиях, в том числе на 14 международных научных конференциях и симпозиумах, на III и IV Всероссийских симпозиумах " Математическое моделирование и компьютерные технологии", "Математическое моделирование и вычислительный эксперимент в естественных и гуманитарных науках"; на Всерос

T(x) = f3Eq(\x«

57) сийской научной конференции "Математическое моделирование в научных исследованиях"; Четвертом Сибирском конгрессе по прикладной и индустриальной математике; на десятой Всероссийской школе-семинаре " Современные проблемы математического моделирования"; на Воронежской весенней математической школе "Понтрягинские чтения", "Современные методы в теории краевых задач" в 1999 и 2001 годах.

Заключение диссертация на тему "Нелокальные математические модели процессов переноса в водоносных природных системах с фрактальной структурой"

Заключение

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

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

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

2. Проведен структурный и качественный анализ математических моделей движения грунтовых вод и почвенной влаги, учитывающих явления последействия и основанных на волновом уравнении, уравнении смешанного эллшггико-гиперболического и параболо-гиперболического типов, уравнении фрактальной диффузии. Предложены конструктивные и высокой точности алгоритмы их компьютерной реализации, учитывающие дискретные эквиваленты необходимых краевых и внутреннекраевых условий типа фундаментального принципа среднего для волнового уравнения [42. С. 165] и уравнения Лапласа, когда приближенное значение гармонической функции в произвольном узле прямоугольной сетки должно совпадать со средним арифметическим значением, принимаемых ее в четырех "соседних" узлах [88. С. 323].

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

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

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

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

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

1. Экстремальные свойства линейных математических моделей процесса движения грунтовых BOflt в основе которых лежит весьма важное в газовой динамике околозвуковых течений уравнение Лаврентьева-Бицадзе [89]; утверждение о том, что нелинейность этого процесса проявляется через некорректность по Адамару задачи Коши в локальной постановке для уравнения Лапласа.

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

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

Основные положения и выводы диссертации были предметом систематического обсуждения на еженедельных заседаниях научно-исследовательского семинара по современному анализу и информатике Института прикладной математики и автоматизации Кабардино-Балкарского научного центра Российской академии наук и прошли апробацию на следующих научных и научно-технических мероприятиях:

1. Научно-техническая конференция профессорско-преподавательского состава СевКавГТУ (Ставрополь, 1992 г.);

3. Третья Международная конференция "Математика, компьютер, образование" (Москва, 1996 г.);

5. Международная конференция "Воздействие интенсивных потоков на вещество" (Терскол, 1999 г.);

6. Третья региональная научно-техническая конференция "Вузовская наука - Северо-Кавказскому региону" (Ставрополь, 1999 г.);

7. Научно-техническая конференция профессорско-преподавательского состава СевКавГТУ (Ставрополь, 1999 г.);

8. Международная конференция "Аналитические методы анализа и дифференциальных уравнений. AMADE" (Минск, 1999 г.);

10. Третий Всероссийский симпозиум "Математическое моделироваг ние и компьютерные технологии", "Математическое моделирование и вычислительный эксперимент в естественных и гуманитарных науках"

Кисловодск, 1999 г.);

14. Четвертый Сибирский конгресс по прикладной и индустриальной математике (ДОГОРИМ) (Новосибирск, 2000 г.);

15. Международная конференция "Аналитические методы анализа и дифференциальных уравнений. AMADE" (Минск, 2001 г.);

16. Международная конференция "Новые подходы к решению дифференциальных уравнений" (Дрогобычи, 2001 г.);

19. Вторая Международная конференция "Нелокальные краевые задачи и родственные проблемы математической биологии, информатики и физики" (Нальчик, 2001 г.);

20. Десятая Международная конференция "Математика. Экономика. Образование" (Ростов-на-Дону, 2002 г.);

21. Четвертая Международная конференция по неравновесным процессам в соплах и струях (Санкт-Петербург, 2002 г.);

23. Международная конференция "Аналитические методы анализа и дифференциальных уравнений. AMADE" (Минск, 2003 г.);

Формула (5.2.11), устанавливающая пропорциональность турбулентного потока дробной производной от удельной влажности на деятельной поверхности земли, прошла особую апробацию: она вошла в монографию A.M. Нахушева [8. С. 265].

В заключение отметим монографии [9], [90] -[95] и работы [96]—[101], посвященные краевым задачам для уравнений параболо-гиперболического типа и некоторым их приложениям к задачам распространения электрических колебаний в составных линиях и теплопереноса в системе, составленной из двух стержней; развитию метода уравнений в частных производных смешанного типа как одного из конструктивных способов математического моделирования физических процессов и явлений в составных (смешанных) средах. "

Библиография Сербина, Людмила Ивановна, диссертация по теме Математическое моделирование, численные методы и комплексы программ

1. Полубаринова-Кочина П.Я. Теория движения грунтовых вод.— М.: Наука, 1977. 664 с.

2. Аверьянов С.Ф. Борьба с засолением орошаемых земель. — М.: Колос, 1978. 288 с.

3. Гидрогеологическое прогнозирование.— Пер. с англ.: Под редакцией М.Г. Андерсена и Т.П. Берта. М.: Мир, 1988. — 736 с.

4. Бэр Я., Заславски Д., Ирмей С. Физико-математические основы фильтрации воды: Пер. с англ. — М.: Мир, 1971. 452 с.

5. Caputo М. Elasticita е Dissipazione.— Zanichelli, Bologna, 1969 (in Italian).

6. Полубаринова-Кочина П.Я., Пряжинская В.Г., Эмих В.Н. Математические методы в вопросах орошения. — М.: Наука, 1969. 414 с.

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

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

9. Джураев Т.Д., Cony ев А., Мамажанов М. Краевые задачи для уравнений параболо-гиперболического типа. — Ташкент: Фан, 1986. -220 с.

10. Елеев В.А. Краевые задачи для уравнений смешанного гиперболо-параболического типа. Автореферат диссертации на соискание ученой степени д.ф.-м.н. Киев, 1995.

11. Нахушев A.M. К теории краевых задач для уравнения второго порядка смешанного гиперболо-параболического типа // Дифференциальные уравнения. 1978. Т. 14, №1.

12. Сабитов К.Б. Некоторые вопросы качественной и спектральной теории уравнений смешанного типа. Автореферат диссертации на соискание ученой степени д.ф.-м.н. Киев, 1992.

13. Нахушев A.M. // Дифференц. уравнения. 1982. Т.18, № 1. -С. 72 - 81.

14. Нерпин С.В., Чудновский А.Ф. Энерго- и массообмен в системе растение-почвагвоздух. JL: Гидрометеоиздат, 1975. 356 с.

15. Чудновский А.Ф. Теплофизика почв. М.: Наука, 1976. - 352 с.

16. Вабищевич П.Н. Метод фиктивных областей в задачах математической физики. — М.: Изд-во МГУ, 1991. —156 с.

17. Вабищевич П.Н., Горбачевский А.Я. Математическое моделирование массопереноса и кристаллизации в пористой среде // Математическое моделирование. 1995. Т.7, №9. - С. 15-22.

18. Крукиер Л.А., Шевченко И.В. Сравнение моделей гравитационного режима течения грунтовых вод // Математическое моделирование. 2002. Т. 14, №2. - С. 51-60.

19. Сухинов А.И. Локально-двумерные схемы для решения многомерных параболических уравнений на вычислительных системах матричного типа // Известия вузов. Математика. 1984. № 11. - С. 45-53.

20. Сухинов А.И. Локально-двумерные схемы для аппроксимации трехмерного уравнения теплопроводности в тороидальных координатах // Известия вузов. Математика. 1996. № 3. - С. 58-67.

21. Шхануков М.Х. О сходимости разностных схем для дифференциальных уравнений с дробной производной // ДАН. 1996. Т. 348, № 6. -С. 746-748.

22. Oldham Keith В., Spanier Jerome. The Fractional Calculus (Theory and Applications of Differentiation and Integration to Arbitrary Order). Academic Press, New York and London, 1974. 233 p.

23. Керефов М.А. Краевые задачи для модифицированного уравнения влагопереноса с дробной по времени производной. Автореферат диссертации на соискание ученой степени к.ф.-м.н. Нальчик, 2000.

24. Нахушев A.M. О непрерывных дифференциальных уравнениях и их разностных аналогах // ДАН СССР. 1988. Т. 300, №4. С. 769-799.

25. Худалов М.З. Разностные методы решения краевых задач для уравнений параболического типа с дробной производной. Автореферат диссертации на соискание ученой степени к.ф.-м.н. Нальчик, 2003.

26. Магомедов К.М. Теоретические основы геотермии. — М.: Наука, 2001.- 277 с.

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

28. Мейланов Р.П. Обобщенные уравнения одномерной фильтрации с дифференцированиями дробной степени // Инженерно-физический журнал. 2001. Т.74, № 2. - С.34-37.

29. Нахушев A.M. Об уравнениях состояния непрерывных одномерных систем и их приложениях. — Нальчик: Логос, 1995. 59 с.

30. Нахушева В.А. Некоторые классы дифференциальных уравнений математических моделей нелокальных физических процессов.— Нальчик: Изд-во КБНЦ РАН, 2002. 100 с.

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

32. Нигматуллин РЖ, Белавин Б.А. // Труды КАИ. 1964. Т. 82. -С. 58-61.

33. Кобеле в B.JI., Романов Е.П., Кобелев Я.Л., Кобелев Л.Я. Недебаев-ская релаксация и диффузия в фрактальном пространстве // ДАН. 1998. Т. 361, № 6. С. 755-758.

34. Ионкин Н.И., Моисеев Е.И. // Дифференц. уравнения. 1979. Т.15, №27.-0.1284-1295.

35. Бахания Н.Н. Об одной особой задаче для уравнения смешанного типа // Труды АН ГрузССР. 1963. Т.З. С. 69-80.

36. Нахушев A.M. Критирий единственности решения задачи Дирихле для уравнений смешанного типа в цилиндрической области // Дифференц. уравнения. 1970. Т.6, № 1. - С. 191-196.

37. Хачев М.М. Первая краевая задача для линейных уравнений смешанного типа. — Нальчик: Эльбрус, 1998. 168 с.

38. Лайхтман Д.П. Физика пограничного слоя атмосферы. — Л.: Ги-дрометеоиздат, 1970. 241 с.

39. Соболев С.Л. Локально-неравновесные модели процессов переноса // Успехи физических наук. 1997. Т. 167, №10. С. 1096-1106.

40. Нахушев A.M. Уравнения математической биологии. — М.: Высшая школа, 1995. 301 с.

41. Лебедев Я.Я. Специальные функции и их приложения.—М.: Гос. изд-во физ.-мат. литературы, 1963. 358 с.

42. Курант Р. Уравнения с частными производными. — М.: Мир, 1964. -830 с.

43. Бицадзе А.В. Некоторые классы уравнений в частных производных.1. М.: Наука, 1981. 448 с.

44. Самарский А. А. Введение в численные методы. — М.: Наука, 1997.- 240 с.

45. Самарский А.А. // Дифференц. уравнения. 1980. Т.16, № 11. -С. 1925-1935.

46. Бицадзе А.В. Уравнения смешанного типа. — М.: Изд-во АН СССР, 1959. 155 с.

47. Джрбаишн М.М. Интегральные преобразования и представления функций в комплексной области. — М.: Наука, 1966. 672 с.

48. Лаврентьев М.М., Савельев Л.Я. Теория операторов и некорректные задачи. — Новосибирск: Институт математики, 1999.

49. Ионкин Н.И. // Дифференц. уравнения. -1977. Т.8, № 2. С.174-294.

50. Хлесткина Н.М., Шагалов В.Ш. // Прикладная математика и техническая физика. 1996. Т.37, № 5. - С.82-92.

51. Чарный И.А. Неустановившееся движение реальной жидкости в трубах. — М.: Недра, 1975. 296 с.

52. Нахушев A.M., Нахушева В.А., Сербина Л.И. О некоторых прикладных аспектах дробного исчисления. Тезисы докладов международной конференции "Воздействие интенсивных потоков на вещество". Тер-скол. 1999. - С. 37-38.

53. Deuvres de Maurice Gevrey. Editions du Centre National de la Recherche Scientifique 15, quai Anatole-Prance-Paris-VII. Paris, 1970. 575 p.

54. Сербина Л.И. Об одной линейной математической модели динамики грунтовых вод // Доклады Адыгской (Черкесской) международной академии наук. 1997. Т.2, №2. - С. 19-22.

55. Сербина Л.И. Об одной краевой задаче для уравнения Лаврентьева-Бицадзе с нелокальным условием Самарского // Доклады Адыгской (Черкесской) международной академии наук. 1998. Т.З, №2. -С. 6-9.

56. Сербина Л.И. Математическая модель эволюции малых возмущений в полубесконечных каналах с пористыми и проницаемыми стенками // Доклады Адыгской (Черкесской) международной академии наук. 1999. Т.4, №1. - С. 60-63.

57. Сербина Л.И. Об одном волновом уравнении с дробной производной, моделирующем динамику малых возмущений. Тезисы докладов Международной конференции "Аналитические методы анализа и дифференциальных уравнений. AMADE-99". Минск. 1999. -С. 203-205.

58. Сербина Л.И. Качественный анализ одной математической модели турбулентной диффузии. Тезисы докладов Воронежской весенней математической школы «Понтрягинские чтения-Х», 1999 г. С. 224.

59. Сербина Л.И. Об одной математической модели динамики взаимодействия деятельной поверхности почвы с приземным слоем атмосферы // Доклады Адыгской (Черкесской) международной академии наук. 2000. Т.5, №1. - С. 85-90.

60. Сербина Л.И. Об одной математической модели процесса трансформации полей температуры и влажности при стационарных условиях. Труды IX Международного симпозиума "Методы дискретных особенностей в задачах математической физики". Орел. 2000. -С. 410 - 420.

61. Сербина Л.И. Об одной математической модели динамики солей впочвах и алгоритме ее компьютерной реализации // Доклады Адыгской (Черкесской) международной академии наук. 2001. Т.5, №2. -С. 72-76.

62. Сербина Л.И. Об одной качественно новой краевой задаче для уравнения Лаврентьева-Бицадзе // Труды Института математики НАН Беларуси. 2001. Т.10. - С. 142-146.

63. Сербина Л.И. О разрешимости одного волнового уравнения с дробной производной. Тезисы докладов Международной конференции "Новые подходы к решению дифференциальных уравнений". Дро-гобычи Украина. - 2001. - С. 134.

64. Сербина Л.И. Решение качественно новых краевых задач для уравнения Лаврентьева-Бицадзе. Труды Международной конференции "Математическое моделирование, статистика и информатика в современном управлении экономикой". Самара. 2001. - С. 193-194.

65. Сербина Л.И. Об одной краевой задаче для математической модели переноса во фрактальной среде. Тезисы докладов Воронежскойвесенней математической пшолы "Понтрягинские чтения"-XII "Современные методы в теории краевых задач". Воронеж. 2001. -С. 141 - 142.

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

67. Сербина Л.И. Анализ математической модели динамики грунтовых вод // Известия вузов. Северо-Кавказский регион. Естественные науки. 2002. № 2. - С. 19-21.

68. Сербина Л.И. Об одной математической модели безнапорного движения грунтовых вод // Известия КБНЦ РАН. 2002. № 1(8). С. 94-102.

69. Сербина Л.И. Об одной проблеме для линеаризованного уравнения Буссинеска с нелокальным условием Самарского // Дифференц. уравнения. 2002. Т.38, N°. 8. - С. 1113-1119.

70. Сербина Л.И. Математическая модель расчета влагосолепереноса во фрактальных средах. Тезисы докладов X Международной конференции «Математика. Экономика. Образование». Ростов-на-Дону, 2002.

71. Сербина Л.И. Об одном волновом уравнении с дробной производной, моделирующем динамику малых возмущений. Тезисы докладов IV международной конференции по неравновесным процессам в соплах и струях. Санкт-Петербург. 2002. - С. 389-391.

72. Сербина Л.И. Об одной модели динамики микрометеорологического режима при орошении // Математическое моделирование. 2002.1. Т.14, №4. С. 67-79.

73. Сербина Л.И. Математическая модель динамики почвенной влаги во фрактальных средах // Вестник Дагестанского научного центра. 2003. №14. С. 19-24.

74. Сербина Л.И. Нелокальная начально-краевая задача для модельного уравнения нестационарной фильтрации // Доклады Адыгской (Черкесской) международной академии наук. 2003. Т. 6, №2. - С. 125130.

75. Сербина Л.И. Об одной математической модели переноса субстанции во фрактальных средах // Математическое моделирование. 2003. Т. 15, №9. С. 17-28.

76. Сербина Л.И., Решение одной начально-краевой задачи теории фильтрации с нелокальным краевым условием // Вестник СамГТУ. Серия: Физико-математические науки. 2003. Вып. 19. С. 16-21.

77. Сербина Л.И. Об одной математической модели влагопереноса в пористых средах с фрактальной структурой // Тезисы докладов международной конференции "AMADE"-2003. 4-9 сент. 2003. Минск. С. 160-161.

78. Трикоми Ф. Лекции по уравнениям в частных производных. М.: Иностранная литература, 1957. - 443 с.

79. Франкль Ф.И. Избранные труды по газовой динамике. М.: Наука, 1973.- 711 с.

80. Калъменов Т.Ш. Краевые задачи для линейных уравнений в частных производных гиперболического типа. Шымкент: Гылым, 1993. -328 с.

81. Кожанов А.И. Краевые задачи для уравнений математической физики нечетного порядка. Новосибирск: Новосиб. университет, 1990. -132 с.

82. Кожанов А. И. Задача сопряжения для одного класса уравнений составного типа переменного направления. В сб. научн. работ " Неклассические уравнения математической физики". Новосибирск: изд-во Института математики, 2002. С. 97-109.

83. Моисеев Е.И. Уравнения смешанного типа со спектральным параметром. М.: изд-во МГУ, 1988. 150 с.

84. Салахитдинов М.С., Уринов А.К. Краевые задачи для уравнений смешанного типа со спектральным параметром. Ташкент: Фан, 1997. 165 с.

85. Смирнов М.М. Уравнения смешанного типа. М.: Наука, 1970. 296 с.

86. Бжихатлов Х.Г., Нахушев A.M. // ДАН СССР. 1968. Т.183, №2. С. 261-264.

87. Елеев В.А. // Дифференц. уравнения. 1994. Т.ЗО, №2. С. 230-237.

88. Елеев В.А. If Украинский математический журнал. 1995. Т.47, № 1. С. 20-29.

89. Елеев В.А. // Украинский математический журнал. 1995. Т.47, № 12. С. 1636-1652.

90. Капустин Н.Ю. Задача для параболо-гиперболического уравнения с нелокальным условием склеивания нормальных производных на линии изменения типа // Доклады АН СССР. 1989. Т. 305, №1. -С. 31-33.

91. Моисеев Е.И., Капустин Н.Ю. // Дифференц. уравнения. 2003. Т.39, №5. С. 656-662.

Похожие работы

Информатика, вычислительная техника и управление
05.13.00