CC BY
48
(Выпуск 1
^ВЕСТНИК
~ гасудАравЕНншУнйв^сшА
х^рЛ0РСК0т*1 очного флота имени адмирала с о. Макарова
МОРСКИЕ И ВНУТРЕННИЕ ВОДНЫЕ ПУТИ, ГИДРОТЕХНИЧЕСКИЕ СООРУЖЕНИЯ И ПОРТЫ
УДК 532.546 Э. Н. Береславский.
д-р физ.-мат. наук, профессор, Санкт-Петербургский государственный университет гражданской авиации;
Б. В. Пестерев.
аспирант,
Санкт-Петербургский государственный университет гражданской авиации
МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ ТЕЧЕНИЙ ИЗ КОТЛОВАНОВ И ВОДОХРАНИЛИЩ
MATHEMATICAL MODELING OF FLOW THROUGH TRENCHES AND WATER RESERVOIRS
В гидродинамической постановке рассматривается плоская установившаяся фильтрация несжимаемой жидкости по закону Дарси из строительных котлованов, огражденных шпунтами Жуковского, через однородный и изотропный слой грунта, подстилаемый хорошо проницаемым напорным водоносным горизонтом, на кровле которого содержится непроницаемый участок. Для изучения инфильтрации на свободную поверхность грунтовых вод формулируется смешанная краевая многопараметрическая задача теории аналитических функций, которая решается с помощью метода Полубариновой-Кочиной и способов конформного отображения областей специального вида, характерных для задач подземной гидромеханики. Рассматриваются предельные случаи течения.
The filtration from the pits, which are fenced rabbets of Zhukovsky through a soil layer consider. At the bottom of the soil is highly permeable pressure aquifer with nonpermeable site. Mixed multiparametric boundary value problem of the theory of analytic functions is formulated to study the infiltration of the free surface. The problem is solved using the Polubarinova-Cochina’s method. The limiting cases are considered. They associated with the lack of one of the factors which characterize the simulated process.
Ключевые слова: фильтрация, грунтовые воды, шпунт, область комплексной скорости, конформные отображения.
Key words: filtration, groundwater, dam, groove, velocity hodograph, conformal mappings.
И
ССЛЕДОВАНИЯ по изучению фильтрационных течений из строительных котлованов и водохранилищ, огражденных симметричными шпунтами Жуковского (см., например, обзоры [1; 2; 3, с. 585-648; 4]). При этом в одних случаях предполагалось, что водопроницаемый слой грунта имеет неограниченную мощность, в других нижележащий сильнопроницаемый напорный пласт моделировался одной или двумя дренами в форме горизонтальной щели Жуковского [5, с. 297-332]. В ряде работ изучалась свободная фильтрация, то есть течение без подпора, а в отдельных случаях — напорная, то есть наличием свободной поверхности пренебрегали. Во всех упомянутых исследованиях учет инфильтрации не производился. Кроме того, для решения задач авторы использовали разные приемы: применялась функция Жуковского и способ Ведерникова-Павловского, сводящие дело к конформному отображению прямолинейных многоугольников с последующим использованием формулы Кристоффеля-Шварца.
В отличие от названных исследований, ниже, как непосредственное продолжение и развитие предыдущих результатов автора [6, с. 529-533], изучается задача о течении жидкости из котлова-
нов через грунтовой массив конечной мощности, подстилаемый хорошо проницаемым напорным водоносным горизонтом, содержащим на своей кровле водонепроницаемый участок, при наличии инфильтрации на свободную поверхность. Исследуется наиболее общий случай движения, при котором на обоих водопроницаемых участках границы области фильтрации расход принимает экстремальные значения и точка нулевой скорости потока выходит на шпунт (что, по-видимому, до сих пор не встречалось в литературе). Отмечаются предельные случаи течения, связанные с отсутствием как подпора, непроницаемого включения или инфильтрации, так и случай вырождения котлованов в полубесконечную слева полосу затопления, исследованный ранее В. В. Ведерниковым [7, с. 619-622]. Приводятся результаты для схемы, возникающей при отсутствии критических точек в том случае, когда скорость обтекания на конце шпунта оказывается конечной; полученное решение является некоторым аналогом классической задачи Жуковского [5].
Для решения смешанной краевой многопараметрической задачи теории аналитических функций используется метод Полубариновой-Кочиной [1-4], а также разработанные для областей специального вида [8] способы конформного отображения круговых многоугольников [9, с. 296301; 10, с. 590-594], которые весьма характерны для задач подземной гидромеханики. Учет характерных особенностей движения позволяет представить решения через специальные, а в ряде случаев элементарные функции, что делает их применение простым и удобным.
На основе построенных точных аналитических зависимостей и посредством численных расчетов проводится гидродинамический анализ влияния всех физических параметров схемы на картину явления и отмечаются некоторые особенности разрабатываемых моделей.
1. Основная модель. Постановка задачи. Рассматривается плоское установившееся течение из котлована (водохранилища) ЛЛ шириной 21, огражденного симметричными шпунтами Жуковского одинаковой длины S, через водопроницаемый слой грунта мощности Т с нижележащим водоносным горизонтом EE, содержащим подземные или артезианские воды, напор в котором имеет постоянное значение Н0 (рис. 1). На кровле этого пласта расположен водонепроницаемый участок D'D, моделируемый горизонтальным отрезком длины 2L. В силу симметрии картины движения ограничимся изучением правой половины области фильтрации ABCDEGR.
Грунтовые воды, обтекая шпунт ARG под действием разности напоров в котловане и нижележащем хорошо проницаемом водоносном горизонте, поднимаются за ним на некоторую высоту RG и, преодолевая точку M нулевой скорости на шпунте, образуют свободную поверхность GE, на которую поступают инфильтрационные воды с интенсивностью s (0 < s < 1), отнесенной к коэффициенту фильтрации грунта к = const. Задача состоит в определении положения кривой депрессии GE и, стало быть, высоты RG поднятия грунтовой воды за шпунтом, то есть величины S - d.
Будем полагать, что движение грунтовых вод подчиняется закону Дарси с известным коэффициентом фильтрации к и происходит в однородном и изотропном грунте, который считается несжимаемым. Скорость обтекания на острие шпунта vR принимается бесконечной (до п. 5), глубина воды Н в котловане остается неизменной во времени.
С математической точки зрения задача состоит в нахождении комплексного потенциала течения ш = ф + /ф (ф — потенциал скорости; ф — функция тока) как аналитической внутри области фильтрации функции комплексной координаты z = x + iy при следующих граничных условиях:
АВ:у = 0, ф = -Н; ВС :х = 0, \|/ = 0; CD:у = -Т,\|/ = 0; AG: х = /,ф = Q\ DE:у = -Т,ф = -Н0; GE:q> = -y-T,\\i = Q + s(x-1),
(1)
где Q — искомый фильтрационный расход из котлована. Исследование осуществляется в терминах приведенных величин ши z, связанных с одноименными фактическими величинами щф и z^i посредством равенств ш = шф/кТ, z = zJT.
Выпуск 1
Выпуск 1
li
Рис. 1
2. Построение решения краевой задачи. Обратимся к области комплексной скорости w, соответствующей краевым условиям (1). Эта область, являющаяся круговым многоугольником с тремя разрезами, вершины N1 и N2 двух из которых соответствуют экстремумам функции тока на непроницаемых участках AB и DE, принадлежит классу многоугольников в полярных сетках [8] и совпадает с таковой для случая, рассмотренного ранее [6; 11, с. 130-138] (рис. 2). Однако, в отличие от [6; 11], в данном случае на границе области движения появляется дополнительная угловая точка — точка B; общее число особых точек становится равным девяти, что сильно осложняет задачу.
Для решения краевой задачи используем метод Полубариновой-Кочиной, который основан на применении аналитической теории линейных дифференциальных уравнений класса Фукса [1-4]. Вводится прямоугольник вспомогательной параметрической переменной плоскости
г (0 < Reх < 1/2, 0 < Imx < 0,5р, p(k) = K'/K, K = K(k'), k' = s[\, где K(k) — полный эллиптический интеграл первого рода при модуле k [27,28] при соответствии точек %е = 0, zQ = 0,5, %с = 0,5(1 + /р), %D = 0,5/р) и функции
d со
Z
dz
dQ
(2)
Определяя характеристические показатели функций (2) около регулярных особых точек [1-4] и принимая во внимание то, что функция w = d®/dz = Q/Z имеет прежний вид [6; 11], придем к зависимостям:
Q = -^N^,Z = iN\Q, %±(х) = (1 + ^(х)±(1-^)У2(х),
А(х) А(х)
Y12 (х) = $q3 ('c)^i (4 ± i‘y)9 2 (т - *а)Э2 (х + г’Р) ехр(+глх), (3)
Д(х) = sd(2Kt, *)^[l -(l-£'2Z2)sn2 (2Kz, *)] [l - (l - k'2B2) sn2 (2 Kz, it)],
Ve = th7t(0,5p + p-a-y). (4)
Здесь sn(M, k) — эллиптическая функция Якоби (синус) при модуле k; SQ(x), ^1(х) и S2(x) — тета-функции с параметром q = ехр(-лр), который однозначно связан с модулем k; N> 0 — масштабная постоянная моделирования; A = sn(2Ka, k'), B = sn(2Kb, k'), a и b — неизвестные ординаты точек A и B области x. В представлениях (3) постоянные конформного отображения a, Р иу, которые связаны соотношением (4), подчинены условиям:
Выпуск 1
Выпуск 1
0 < а < r < Р < m < a < b < 0,5 p; 0 < у < 0,5 p,
(5)
регламентирующим положение на границе области течения точек нулевой скорости M и острия шпунта R, а также N1 и N2; m и r — неизвестные ординаты точек M и R в плоскости т.
Можно проверить, что функции (3) удовлетворяют условия (1), переформулированные в терминах функций d<s>/dx и dz/dx, таким образом являющиеся параметрическим решением исходной краевой задачи.
Запись соотношений (3) для разных участков границы области т с последующим интегрированием по всему контуру вспомогательной области приводит к замыканию области движения и тем самым служит контролем вычислений.
В результате получаем выражения для задаваемых и искомых геометрических и фильтрационных характеристик модели:
0,5
0,5р
J ^RAdt — J XAedt — h j XCDdt — L, J
yBCdt -
0,5p 0,5
J j &CDdt = H — H0,
ь о
(6)
0,5
d = T — H0 — J ФEGdt о
(7)
и координат точек кривой депрессии EG:
0,5 0,5
xeg(м) — l + j XEGdt, уEG(и) = -d+\ YEGdt, 0 < u < 0,5.
(8)
Контролем счета являются другие выражения для величин d, L и фильтрационного расхода Q:
0,5 0,5 0,5р
d = Т — Н0 — J YEGdt, L = l — j XEGdt — J XDEdt,
ooo (9)
0,5p 0,5
Q= j *YDEdt — e j XEGdt.
о 0
В формулах (6)-(9) подынтегральные функции — выражения правых частей равенств (3) на соответствующих участках контура области т.
3. Анализ численных результатов для основной фильтрационной модели. Представления (3), (6)-(9) содержат семь неизвестных постоянных: ординаты a, b, r прообразов точек A, B, R в плоскости т, параметры конформного отображения а, Р, у, удовлетворяющие соотношению (4) и неравенствам (5), а также модуль k (0 < k < 1) и постоянную моделирования N. Для их определения при заданных S, l, L, H и T служит система уравнений (6), наряду с которыми используются соотношения:
w 1(0,5 + г>) = 0,
0,5
j(®
0,5р
EG +
)dt + joGAdt+ J ФдСЛ = 0.
(10)
Первое из этих соотношений означает, что скорость на конце шпунта обращается в бесконечность, а второе непосредственно вытекает из рассмотрения граничных условий (1). После определения неизвестных постоянных находятся искомые величины d и Q по формулам (7) и, наконец, по формулам (8) рассчитываются координаты точек свободной поверхности EG.
На рис. 1 изображена картина течения, рассчитанная при s = 0,6, T = 7, S = 3, H0 = 3, L = 15, H = 7, l = 10 (базовые значения). В табл. 1 и 2 (варьируется в допустимом диапазоне один из указанных параметров, а остальные фиксируются базовыми значениями) приведены результаты расчетов влияния определяющих физических параметров s, T, S, H0, L, H и l на глубину d (отрицательные величины d означают, что свободная поверхность поднимается выше оси абсцисс) и расход Q.
Таблица 1
£ d Q T d Q S d Q Н0 d Q
0,5 2,651 0,182 6,5 1,349 0,234 2,0 1,726 0,635 2,0 3,155 0,038
0,6 1,804 0,394 7,0 1,804 0,394 2,5 1,745 0563 4,0 0,441 0,769
0,7 0,586 0,457 7,5 2,299 0,457 3,5 1,844 0255 5,0 -0,93 1,159
0,8 -1,195 0,094 8,0 2,745 0,627 4,0 1,873 0,129 6,0 -2,35 1,815
Таблица 2
L d Q Н d Q l d Q
14 2,555 0,086 4,0 0,650 1,815 10,0 1,804 0,394
16 1116 0,599 5,0 1,070 1,159 10,3 2,018 0,320
17 0,413 0,888 6,0 1,441 0,769 10,6 2,229 0,253
18 -0,256 1,084 8,0 2,155 0,038 11,0 2,478 0,230
Анализ данных таблиц и графиков позволяет сделать следующие выводы.
Прежде всего обращает на себя внимание одинаковый качественный характер зависимостей величин d и Q от параметров T и l, S и Н и в то же время противоположное поведение искомых характеристик при изменении параметров S и H с одной стороны йот L и Н0 — с другой.
Увеличение интенсивности инфильтрации, ширины непроницаемого включения и напора в нижележащем пласте и уменьшение мощности слоя, длины шпунта, напора воды в котловане и его ширины приводят к уменьшению глубины d, то есть к увеличению ординаты точки G выхода кривой депрессии из-под шпунта. Так, согласно данным табл. 1 и 2 увеличению параметров s, S, Н0 и l в 1,6, 2,0, 2,0 и 1,1 раза сопутствует изменение величины d в 2,2, 1,1, 1,3 и 1,6 раза соответственно. Однако наибольшее влияние на глубину d оказывает непроницаемый участок: данные табл. 2 показывают, что при возрастании ширины L всего на 28 % глубина dувеличивается почти в 10 раз.
При е = 0,8, Н0 = 5 и 6 и L = 18, то есть при достаточно больших значениях параметров е, Н0 и L, свободная поверхность поднимается выше оси абсцисс, при этом величина d принимает отрицательные значения. Если ввести безразмерную величину h(d) = (S - d)/S, h(S) = 0, характеризующую относительную высоту поднятия грунтовой воды за шпунтом, то для указанных значений параметров е, Н0 и L получаем h(-1,1945) = 1,3981, h(-0,9297) = 1,3099, h(-2,3500) = 1,7833 и h(-0,2560) = 1,0853 соответственно, причем величина h растет с увеличением ординаты точки G выхода кривой депрессии из-под шпунта.
Зависимости глубины d от параметров T, Н0, L и Н близки к линейным.
Что касается расхода, то с увеличением ширины непроницаемого включения величина Q также возрастает: из табл. 2 следует, что увеличение параметра L в 1,28 раза влечет увеличение расхода более чем в 12 раз. Таким образом, выявляется значительное подпирающее воздействие непроницаемого участка по отношению к фильтрации из котлована.
Подобное поведение расхода наглядно прослеживается как с увеличением мощности слоя, ширины котлована и напора в нижележащем горизонте, так и с уменьшением длины шпунта и напора в бьефе. Из табл. 2 вытекает, что уменьшению параметраН всего в 2 раза сопутствует увеличение расхода Q почти в 48 раз, что свидетельствует о наибольшем влиянии на расход напора воды в котловане
4. Предельные случаи. 1. Случай Н0 = 0. Остановимся прежде всего на случае отсутствия подпора, то есть движении при Н0 = 0. Решение для этого предельного случая получается из зависимостей (3), (6)-(9) при у = у, = 0. При таком значении параметра у круговой разрез EG области w, трансформируясь, вырождается в правую полуокружность (штриховая линия на рис. 2), и таким образом исходная область превращается в круговой шестиугольник, у которого выпадает правая часть полукруга |w - 0,5/(1 + е)| < 0,5(1 - е). В плоскости течения z при у = у, кривая депрес-
Выпуск 1
Выпуск 1
сии выполаживается в точке E, с которой сливается ее точка перегиба F, и выходит на кровлю нижележащего горизонта под прямым углом.
2. Случай l = да. Рассмотрим случай, когда ширина котлована неограниченно возрастает. Если совершить преобразование z' = z + l, переведя точку A' в начало координат, и зафиксировать все физические параметры модели, то по мере возрастания ширины котлована l постоянная конформного отображения b ^ b* = 0,5 р. В пределе при l = да параметры b = b*, B = 1. В плоскости движения z точки B и C сливаются на бесконечности, так что область фильтрации становится полубесконечной слева полосой затопления. Из (3), (6)-(9) при b = b* вытекают результаты, полученные ранее [6; 11].
3. Случай L = 0. Если отсутствует непроницаемый участок на кровле нижележащего водоносного горизонта, то последний на всем протяжении становится хорошо проницаемым. При слиянии точек C и D в области комплексной скорости w ее левая полуплоскость отсекается, круговой разрез EG переходит в правую полуплоскость, а исходная область трансформируется в круговой треугольник (рис. 3). В плоскости движения z точка D, сливаясь с точкой C, выходит на ось ординат, а прямоугольник плоскости т преобразуется в полуполосу 0 < Rex < 0,5, 0 < Imx < да, поскольку параметр р =K'/K = да, K = 0,5 л.
Решение для этого предельного случая получается из формул (3), (6)-(9), если в них положить k = 0 и учесть, что при этом эллиптические функции вырождаются в тригонометрические, а тета-функции обрываются на своих первых членах или константах:
sin2m'cosxA(x)
Z=iAfsi° 2(1
cosxA(x)
Д(т) = sj(a- sin2 x)(b - sin2 x)(c - sin2 x), e = tg2m'ctg2r'.
(11)
Здесь /w' = arcsin 4m, r' = arcsin 4~r , m', r', c' (0 < m' < r' < a' < b' < c' < 0,5л) — прообразы точек M, R, C на оси абсцисс плоскости т. Этот случай подробно описан ранее [12, с. 1032-1040], где дан анализ влияния всех физических параметров модели.
4. Случай е = 0. В рамках задачи L = 0 остановимся на случае отсутствия инфильтрации. Учитывая связь между параметрами т' и г', отмеченную в (11), видим, что в случае, когда е = 0, решение задачи вытекает из зависимостей (11) при т' = 0, то есть когда в плоскости т точки кривой G и E сливаются в начале координат с точкой M нулевой скорости.
Таким образом, получается решение задачи, впервые рассмотренной В. В. Ведерниковым [7], только другим способом.
5. Случай конечной величины скорости обтекания на конце шпунта. Анализ задачи Жуковского. В рамках краевой задачи (1) рассмотрим случай, когда скорость обтекания на конце шпунта vR, 0 < vR < е конечна и функция тока на водопроницаемых участках AB и DE не имеет экстремумов. Тогда в области комплексной скорости w исчезают оба вертикальных разреза, левая полуплоскость отсекается, как и ранее в случае L = 0, однако, в отличие от последнего, участок MR переносится в первый квадрант (штриховая линия на рис. 3). В результате исходная область превращается в круговой пятиугольник.
Параметрическое решение задачи имеет прежний вид (3) с заменой интегралов Т12(т) и постоянных конформного отображения а и р на следующие:
Решение аналогичной задачи для случая отсутствия подпора вытекает из представлений (3), (12) при у = у,.
Анализ проведенных численных результатов показывает, что в случае vR < да сохраняется качественный характер зависимостей фильтрационного расхода от физических параметров схемы, свойственный случаю vR = да. Так, например, наблюдается тот же, что и ранее, характер поведения расхода от величин T и l с одной стороны и противоположный характер от параметров S и H — с другой. Существенное влияние на расход Q по-прежнему оказывают инфильтрация, непроницаемое включение и мощность слоя.
На рис. 4 изображена картина движения, рассчитанная при s = 0,5, T = 6, S = 3, L = 16,2, H = 3,
Обращает на себя внимание то обстоятельство, что для всех расчетных вариантов оказывается d = S и, следовательно, величина h(d) = h(S) = 0. Это означает, что в плоскости течения точка G выхода кривой депрессии из-под шпунта сливается с точкой R его острия; из рассмотрения области комплексной скорости w (рис. 3) следует, что в таком случае скорость на конце шпунта равна интенсивности инфильтрации: vR = 8, 0 < 8 < 1.
Если совершить преобразование т' = 0,5 + /р'т, переводящее прямоугольник вспомогательной переменной т в ему подобный с параметром р' = 1/р = KIK', то соответствующее основной фильтрационной схеме неравенство на параметры (5) примет вид
\2 ("0 = V(T)9i0 ± *у) ехр(+гат), а = Р = 0,5(1 + гр).
(12)
l = 15.
0 < b' < a' < r' < 0,5,
где b', a', r' — абсциссы прообразов точек B, A, R в плоскости т.
(13)
\\
i
Э
©
н
В
A
10
15
18 Л-
S,d
-a-
Puc. 4
Выпуск 1
(Выпуск 1
^ВЕСТНИК
-- “7 ГОСУДАРСТВЕННСИТ) университета
^JMACKOm И РЕЧНОГО ОЛОГЛ ИМЕНИ АДМИРАЛА С О, МАКАРОВА
Расчеты показывают, что для любого значения интенсивности инфильтрации s (0 < s < 1) соотношение d = S выполняется лишь для единственного значения величины г' — ее предельного значения г*, когда на плоскости т' сливаются точки G и R: r = г* = 0,5. Все другие допустимые значения r < г* приводят к несогласованности с реальной картиной течения — соотношению d > S, то есть к отрыву потока.
Подобный результат в предельном для рассматриваемой модели случае, когда водопроницаемый слой грунта имеет неограниченную мощность, отсутствует непроницаемый участок и инфильтрация, то есть при T = да (к' = 0, к = 1), L = 0 (b' = b* = 0) ие = 0 (т' = 0), был впервые получен в свое время Н. Е. Жуковским [5]. Решение для этого предельного случая получается из зависимостей (3), (12), если в них положить K = да, K = л/2, к' = 0, к = 1, b' = 0, q' = 0 и учесть, что при этом эллиптические функции вырождаются в гиперболические, а тета-функции, которые на этот раз характеризуются параметром q' = 0, обрываются на своих первых членах или константах.
Таким образом, в предельном случае изучаемой схемы получается решение задачи Жуковского, только иным путем.
Список литературы
1. Полубаринова-Кочина П. Я. Теория движения грунтовых вод / П. Я. Полубаринова-Кочи-на. — М.: Гостехиздат, 1952. — 676 с.; 2-е изд. — М.: Наука, 1977. — 664 с.
2. Развитие исследований по теории фильтрации в СССР (1917-1967). — М.: Наука, 1969. —
545 с.
3. Михайлов Г. К. Движение жидкостей и газов в пористых средах / Г. К. Михайлов, В. Н. Николаевский // Механика в СССР за 50 лет. — М.: Наука, 1970. — Т. 2.
4. Кочина П. Я. Избранные тр. Гидродинамика и теория фильтрации / П. Я. Кочина. — М.: Наука, 1991. — 351 с.
5. Жуковский Н. Е. Просачивание воды через плотины / Н. Е. Жуковский // Жуков -ский Н. Е. Собр. соч. — М.: Гостехиздат, 1950. — Т. 7.
6. Береславский Э. Н. О некоторых гидродинамических моделях, связанных с задачей Жуковского об обтекании шпунта / Э. Н. Береславский // Докл. РАН. — 2013. — Т. 448, № 5.
7. Ведерников В. В. Фильтрация при наличии дренирующего или водоносного слоя / В. В. Ведерников // Докл. АН СССР. — 1949. — Т. 69, № 5.
8. Кочина П. Я. Аналитическая теория линейных дифференциальных уравнений класса Фукса и некоторые задачи подземной гидромеханики / П. Я. Кочина, Э. Н. Береславский, Н. Н. Кочина. — М.: Ин-т проблем механики РАН, 1996. — Ч. 1. — 122 с. — (Препринт № 567).
9. Береславский Э. Н. О дифференциальных уравнениях класса Фукса, связанных с конформным отображением круговых многоугольников в полярных сетках / Э. Н. Береславский // Дифференциальные уравнения. — 1997. — Т. 33, № 3.
10. Береславский Э. Н. О некоторых дифференциальных уравнениях класса Фукса, встречающихся в задачах механики жидкостей и газов / Э. Н. Береславский // Дифференциальные уравнения. — 2012. — Т. 48, № 4.
11. Береславский Э. Н. О некоторых гидродинамических схемах, связанных с обтеканием шпунта Жуковского / Э. Н. Береславский, Е. В. Пестерев // Вестник СПбГУ. Сер. 1. «Математика, механика, астрономия». — 2013. — Вып. 1.
12. Береславский Э. Н. О фильтрации жидкости из котлована, огражденного шпунтами / Э. Н. Береславский // ИФЖ. — 2013. — Т. 86, № 5.
13. Аравин В. И. Теория движения жидкостей и газов в недеформируемой пористой среде /
В. И. Аравин, С. Н. Нумеров. — М.: Гостехиздат, 1953. — 616 с.
