СМЕШАННАЯ КРАЕВАЯ ЗАДАЧА ДЛЯ ОБОБЩЕННОГО УРАВНЕНИЯ БУССИНЕСКА

Раздел: Научно-практическая работа студентов педагогического вуза

Журнал: Научно-практическая деятельность студентов педагогического вуза

13 июня 2012 г.

Авторы: Маланьина И. Г.

И. Г. Маланьина

Научный руководитель: Фураев Валерий Зиновьевич.

СМЕШАННАЯ КРАЕВАЯ ЗАДАЧА ДЛЯ ОБОБЩЕННОГО УРАВНЕНИЯ БУССИНЕСКА

Рассмотрим начально-краевую задачу для уравнения

с неоднородными граничными условиями

где

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

                           (3)

Введем функцию

               

И записываем уравнение  (1) в виде

                                      (4)

Вводя функцию Грина   для оператора    с нулевыми условиями на концах отрезка [0,1],  из (4) получим:

где

Если начальная функция  непрерывна и неотрицательна на отрезке [0,1] то задача (5), (3) однозначно разрешима в целом, т. е. на отрезке [0,T] для любого T>0 [1. С. 2014 - 2015]. 

Таким образом, данная смешанная начально-краевая задача сводится к задаче Коши для интегро-дифференциального уравнения (5) с начальным условием (3).

         Решение данной задачи может быть найдено методом Эйлера, по формуле:

где

Для численного решения задачи интеграл в правой части по пространственной переменной y вычисляется по формуле прямоугольников [2. С. 24 - 26]:

         Рассмотрим случай непроницаемости левой границы области фильтрации, т.е. при d=0 при этом уравнение (5) примет вид:

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

Приведем пример решения смешанной краевой задачи для обобщенного уравнения Буссинеска с использованием пакета MathCad:

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

 

 

 

 

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

Рис. 1

Здесь жирные сплошные линии – начальная (при t=0) и конечная (по истечении большого промежутка времени t) формы свободной поверхности фильтрующейся жидкости соответственно, остальные – ее промежуточные формы. Из рисунка видно, что начальная форма поверхности с течением времени сходится к константе, которая равна заданному уровню жидкости на правой границе  области фильтрации. В данном примере H=6.

Исходя из анализа расчетов и графических представлений движения жидкости со свободной поверхностью в случае непроницаемости левой границы области фильтрации и заданном уровне жидкости на правой границе заметили, что решение интегро – дифференциального уравнения (5) равномерно сходится по x к константе, значение которой равно уровню жидкости H на правой границе области фильтрации.

Библиографический список

1. Фураев В. З. О разрешимости в целом первой краевой задачи для обобщенного уравнения Буссинеска. // Дифференциальные уравнения. 1983. т.19. № 11. С. 2014 – 2015.

2. Фураев В. З., Чухно С. А. К вопросу о численном интегрировании обобщенного уравнения Буссинеска// Краевые задачи и математическое моделирование: сборник материалов VII Всероссийской конференции 4-5 декабря 2004 г. Новокузнецк: Изд-во КемГУ, 2004. С. 24-26.

PDF