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

УДК 532.546

И. К. Волянская, А. С. Шашков, А. Я. Шпилевой

МОДЕЛИРОВАНИЕ ФИЛЬТРАЦИОННЫХ ТЕЧЕНИЙ В ОБЛАСТИ, ОГРАНИЧЕННОЙ СТОРОНАМИ ПРЯМОУГОЛЬНИКА

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

The article uses the method of spatial point reflection for the definition of complex potentials of filtrational flows in an area limited by rectangle. The article considers possible conditions and generates solutions in the form of infinite series. For point source the series are presented by means of the Jacobi theta-function or the Weierstrass sigmafunction.

Ключевые слова: комплексные потенциалы, фильтрационное течение, точечный источник, тета-функция Якоби, сигма-функция Вей-ерштрасса.

Keywords: complex potentials, filtrational flows, point source, the Jacoby theta-function, the Weierstrass sigma-function.

1. Для исследования фильтрационных течений жидкости, подчиняющихся закону Дарси, используется комплексный потенциал W(z) = q>(x, y) + ip(x, y), с помощью которого можно определить поле скоростей и поле давлений [1].

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

а б

Рис. 1. Область фильтрации с однородными граничными условиями

A y 1 y y

b_ /////////////// ь.

Вестник Россииск ого государственного ун • z0 0 a Г77Т7Т7Т777Т777 иверсит гета им. и. Канта. zuu * z0 0 a .Вып. 4. С. 12-1/.

Jy

ь >

I I

.WWWWWWVs,

б

а

Рис. 2. Область фильтрации с различными граничными условиями

Аналитическим выражением метода отражений особых точек являются теоремы о прямой и окружности [1; 2]. Используем только теорему о прямой.

Пусть границей области фильтрации является ось у, а областью фильтрации полуплоскость Ие z<0. Полагаем, что особые точки функции /(х) находятся в области Ие z <0. Если граница х=0 непроницаема (+) или заполнена свободной жидкостью (-), то теорема о прямой имеет вид

Ш

Если границей раздела является ось х, то аргументом второго сла-

гаемого выражения (1) вместо -г будет z.

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

3. Пусть область фильтрации окружена свободной жидкостью (рис. 1, б) [3]. Применяя метод отражений особых точек, получаем:

ГО _ _

Ш

к ,ш=-ю

где Т = Тк ш = 2ка + 2тЫ .

Пусть в области фильтрации находится точечный сток (скважина). В этом случае

/

2п

Записывая выражение (2) для функции (3) и преобразуя результат, получаем:

Ш

2п ст

где o(z) — сигма-функция Вейерштрасса. Аналогичный результат получен другим методом в работе [4] для электрического поля заряда, расположенного внутри прямоугольника.

В работе [5] получено решение этой задачи с помощью тета-функции Якоби

W

2п 3ü

z - z0 + a z + z0 + a z - z0 + a z + z0 + a

где v =------0---,v2 =-----0----,v3 =----0---,v4 =---------0-.

2ib 2ib 2ib 2ib

4. Пусть область фильтрации ограничена непроницаемыми границами (1, а). Зададим в области фильтрации непроницаемые источник и сток:

f (z) = -Q-^n(z - z0) - -Q-ln(z - z¡) .

2n 2n

Решение задачи для этого случая имеет вид [2]

W

2n cr

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

f (z) = -Q-ln(z - z0).

2п

Решение задачи находим в два этапа.

Сначала находим комплексный потенциал точечного источника в области между двумя непроницаемыми стенками y=0 и y=b [2]:

п(z - z0) п(z - z0) , .

W

2ib 2ib

Если границами являются прямые x=0 и x=a и область фильтрации

ограничена свободной жидкостью, то для функции f ( z) решение

имеет вид

W

k=-w

Записывая выражение (8) для функции (7), получаем комплексный потенциал источника в заданной области (рис. 2, а):

. П . П _

„ sin

W

П П _

k=-” sin

2bi 2bi 0

Выражение (9) преобразуем, пользуясь определением тета-функции Якоби:

W H]l[l-q2mle2nv\-q2mle-2nv\ Ha = fl

m=1 m=1

В результате получаем:

ж (г) = о_ 1п , ?=е-ь,

2п 50

аргументы VI, V2, vз, V4 такие же, как в выражении (5).

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

Ж(2) = £(-1)к+” [7 (2 + Т) + 7 (-2 + Т) + ¡(2 + Т) + 7 (-г + Т)],

Рис. 3. Положение особых точек: * — источник, * — сток

Решение (12), записанное для точечного источника, представим в виде

1п(2 - 20 - Т) + 1п(2 - 20 - Т) + 1п(2 + 20 - Т)+ '

+ 1п( + 20 - Т) - 1п( - 2. - Т) - 1п( - 2[ - Т)-

- 1п( + 2[ - Т) - 1п( + 2. - Т) - 1п( - 22 -Т )-

- 1 2 - Т) - 1п( + 22 - Т) - 1п( + 22 -Т)+

- ( + - Т) - 1п( - 2з - Т) + 1п( + 2, -Т)+

+ 1п

где

Ж (2) = £

г1 = г0 - ТЬ1, 2[ = 20 - ТЫ, г2 = г0 - 2а, г3 = г0 - 2а - 2Ь/,

= 20 - 2а - 2Ы, Т = 4ка + 4шЬ1.

Воспользуемся определением сигма-функции Вейерштрасса:

(2) =2 п|.- ті ехР

( _ 2 2 2

----1------;

Т 2Т'

, к, да = ±1,±2,.

Каждое из слагаемых выражения (13) преобразуем с помощью формулы (14); тогда решение задачи принимает вид

Ж(2) = О- 1п[ст(2 - 20)ст(2 - 20)ст(2 + 20)ст(2 + 20)ст(2 - 2,)ст(2 - 2, ) * 2п

* СТ(2 + 2,)ст(2 + 2,)/ст(2 - 2.)ст(2 + 2.)ст(2 - 2[' )ст(2 + 2[ )ст(2 - 22 ) * *СТ(2 + 22)ст(2 - 22)ст(2 + 2)] .

к.т=-да

к.т=-да

7. Определим комплексный потенциал для точечного источника в области, указанной на рисунке 2, г. Задача решается аналогично предыдущей.

Общее решение задачи, полученное методом отражений для произвольной аналитической функции / (2), имеет вид

Ш(£) = £

/

Решение (16) записываем для точечного источника и преобразуем к вицу

№

2п ст

где 21 = 20 + 2а.

8. Определим комплексный потенциал точечного источника для области, указанной на рисунке 2, в. Используя метод отражении особых точек для аналитической функции / (2), получаем общее решение в виде

№(2) = £

к,т=-«

/(2 -7) - /(-2 +7) -/2 -7) + /(-2+7) - /(2+Т-2а) +

+/

где Т имеет такой же вид, как в п. 6, 7. Решение (18) записываем для функции / (2) = 1п(2 - 20) и результат преобразуем к виду

№

2п а

где 21 = 20 + 2а.

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

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

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

Список литературы

1. Голубева О. В. Курс механики сплошных сред. М., 1972.

2. Зайцев А.А, Шпилевой А. Я. Теория стационарных физических полей в кусочно-однородных средах. Калининград, 2001.

3. Шпилевой А. Я. Моделирование фильтрационных течений жидкости в кусочно-однородных средах методом изображения особых точек // Материалы всероссийской научно-практической конференции: «Вклад земляков-орловцев в развитии и становление российской науки, культуры и образования». Орел, 2003. Т. 3. С. 134—136.

4. Лаврентьев М. А., Шабат Б. В. Методы теории функций комплексного переменного. М., 1965.

5. Севастьянова Н. В., Шпилевой А. Я. Решение задач теории фильтрации, выраженные через тета-функцию Якоби для областей в виде прямоугольника и равностороннего треугольника // Труды международной конференции: «Современные методы физико-математических наук». Орел, 2006. Т. 2. С. 67—70.

Об авторах

А. Я. Шпилевой — канд. физ.- мат. наук, доц., РГУ им. И. Канта.

И. К. Волянская — студ., РГУ им. И. Канта.

А. С. Шашков — учащийся малого физ.-тех. факультета, РГУ им. И. Канта и лицея № 18.

Authors

A. Shpilevoy — Dr., IKSUR.

I. Volyanskaya — student, IKSUR.

A. Shashkov — student, IKSUR.