CC BY
9
4. Zelenetskaya E. P. Model of mass transfer processes in the cascade of centrifugal extractors / E P Zelenetskaya // 2016 IOP Conf. Ser.: Mater. Sci. Eng. 135 012053 doi:10.1088/1757-899X/135/1/012053.
5. Abramzon A. A. Emulsii [Emulsions] / A. A. Abramzon; trans. with English; edited by A. A. Abramzova. - M.: Khimiya. 1972. - 443 p. [in Russian]
6. Planovskiy A. N. Protsessy i apparaty khimicheskoy i neftekhimicheskoy tekhnologii [Processes and apparatuses of chemical and petrochemical technology] / A. N. Planovskiy. P.I. Nikolayev // 3rd edition., revised and additional. - M.: Khimiya. 1987. - 496 p. [in Russian]
7. Kleyton V. Emulsii. Ikh teoriya i tekhnicheskiye primeneniya: spravochnoye rukovodstvo po obshchim voprosam issledovaniya emulsiy i ikh ispolzovaniya v tekhnike [Emulsions. Their theory and technical applications: reference manual on general issues of emulsion research and their use in engineering] / V. Kleyton; trans. with English A. Pleteneva and other. -Moskow: Izdatelstvo inostrannoy literatury. 1950. - 680 p. [in Russian]
8. Zimon A.D. Kolloidnaya khimiya [Colloid chemistry]: uchebnik dlya vuzov / A. D. Zimon, N. F. Leshchenko. - M.: Agar. 2003. - 319 p. [in Russian]
9. Novikov M. A. Strukturnyye osobennosti prirodnykh vodoneftyanykh emulsify [Structural features of natural water-oil emulsions] / M. A. Novikov. - M.:. 2007. - 85 p. [in Russian]
10. Kasatkin A. G. Osnovnyye protsessy i apparaty khimicheskoy tekhnologii [The main processes and apparatus of chemical technology] / A. G. Kasatkin. - M.: Alians. 2004. - 753 p. [in Russian]
11. Planovskiy A. N. Protsessy i apparaty khimicheskoy tekhnologii [Processes and apparatuses of chemical technology]: Uchebnik dlya vuzov / A. N. Planovskiy. V. M. Ramm. S. Z. Kagan. - M.: Goskhimizdat. 1962. — 841 p. [in Russian]
12. Analiticheskaya khimiya urana [Analytical chemistry of uranium] / edition by D. I. Ryabchikova. M. M. Senyavina. -M.: Izd. AN SSSR. 1962. - 432 p. [in Russian]
13. Goryunov A. G. Matematicheskoye modelirovaniye protsessa ekstraktsii urana kak obyekta upravleniya [Mathematical modeling of uranium extraction as a control object]: uch. posobiye / A. G. Goryunov, S. N. Liventsov, A. A. Lysenok and others // Tomsk: Izd-vo Tomskogo politekhnicheskogo universiteta. 2008. - 143 p. [in Russian]
DOI: https://doi.org/10.23670/IRJ.2018.71.028 Телятников И.С.
ORCID: 0000-0001-8500-2133, Кандидат физико-математических наук, Федеральный исследовательский центр «Южный научный центр Российской академии наук», Ростов-на-Дону, Россия
Исследование выполнено в рамках реализации ГЗ ЮНЦ РАН, проект № 01201354241 и при поддержке РФФИ и
администрации Краснодарского края (проект 16-41-230184) ОБ ОДНОЙ МОДЕЛИ ВЗАИМОДЕЙСТВИЯ ЛИТОСФЕРНЫХ СТРУКТУР
Аннотация
Получил развитие эффективный подход к исследованию напряженно-деформированного состояния литосферных структур, моделируемых системой пластин Кирхгофа на упругом слое, вблизи разломов. Реализован новый метод решения задачи, соответствующей модели прямолинейного разлома, в плоской постановке, проведены расчеты и проанализированы полученные результаты при варьировании физико-механических характеристик пластин и условий их контакта. Использование изложенного подхода позволит упростить алгоритм определения смещений покрытия и контактных напряжений, а также сделать выводы о влиянии типа разлома (граничных условий на стыке пластин) на характер волнового процесса в моделируемой геофизической среде, в частности изменение формы сигнала при прохождении через разлом, применимые для изучения характеристик разломов в верхней части коры Земли.
Ключевые слова: деформируемое основание, разлом, составное покрытие, вибрационная нагрузка, метод факторизации.
Telyatnikov I.S.
ORCID: 0000-0001-8500-2133, PhD in Physics and Mathematics, Federal Research Centre «Southern Scientific Centre of Russian Academy of sciences», Rostov-on-Don, Russia This work was carried out within the framework of the implementation of the State Order in SSC ofRAS, project No.
01201354241 and supported by the Russian Fundamental Research Fund and the administration of the Krasnodar Krai
(project 16-41-230184) ON A MODEL OF INTERACTION OF LITHOSPHERIC STRUCTURES
Abstract
The paper presents the effective approach to the study of the stress-strain state of lithospheric structures simulated by the Kirchhoff plate system on the elastic layer near the splits. A new method for solving the problem corresponding to the rectilinear fault model was implemented in a flat setting, calculations were performed, and the obtained results were analyzed with varying physical and mechanical characteristics of the plates and the conditions for their contact. The use of this approach will simplify the algorithm for determining the displacements of the coating and contact stresses, and also draw conclusions about the effect of the fault type (boundary conditions at the plate interface) on the nature of the wave process in the simulated geophysical environment, in particular, the change in the shape of the signal when passing through a split, characteristics of splits in the upper part of the Earth's crust.
Keywords: deformable foundation, fault, composite coating, vibration load, factorization method.
В рамках основанных на применении концепции блоковой делимости литосферы моделей [1], [2] в настоящее время получены интересные результаты, характеризующие зарождения и развития сейсмических событий. Зачастую эпицентры землетрясений приурочены к границам блоков различных рангов [3]. Возможность
возникновения землетрясения, вызванного концентрацией сейсмогенерирующих напряжений в области сближения и контакта литосферных структур и названного «стартовым», исследована в работе [4]. При этом в исследовании краевых задач для систем уравнений с частными производными (УЧП) использован метод блочного элемента, алгоритм которого изложен в [5], [6].
В работе развит эффективный подход к исследованию напряженно-деформированного состояния литосферных структур, моделируемых пластинами Кирхгофа на упругом основании, вблизи разломов. Реализован метод решения задачи, описывающей одну из моделей для случая прямолинейного разлома, в плоской постановке. Скалярный случай подобной задачи был рассмотрен автором в [7], при этом использован метод решения, связанный с преобразованием дифференциального оператора задачи. Предлагаемый в данной работе подход, основанный на упрощенном методе блочного элемента, является более универсальным в сравнении с [7] и одновременно обеспечивает получение более простых представлений решения задачи, чем топологический метод в [6, С. 16-18]. В результате численной реализации разработанного алгоритма проведены расчеты и проанализированы полученные результаты при варьировании физико-механических характеристик пластин и условий их контакта.
Рассматривается задача о вибрации двух занимающих полуплоскости контактирующих пластин на поверхности упругого основания под действием поверхностной сосредоточенной нагрузки ^(х) = ^(х — х0)ехр(— ).
Декартова система координат связана со срединной плоскостью покрытия, ось Ох3 направлена по нормали вверх. В
плоском случае установившийся процесс вертикальных колебаний пластин с частотой СО на упругом слое, жестко сцепленном с недеформируемым основанием, может быть описан с помощью линеаризованных уравнений [8, С. 27-28]
Я±(дх1)и±(х) — ^±,5^±(х) = Ь±(х), ±х >0, (1)
где и+ (х ) и (х ) - амплитуды перемещений и контактных напряжений на границе покрытия и подложки;
£-*±,; ь(х,)=— м±(х.); „=& ^= V -
коэффициенты Пуассона, р+ - плотности, - толщины пластин, - модули Юнга их материалов. Здесь и далее
использованы обозначения работы [7, С. 3], временной множитель ехр (—iat) опущен.
Общий вид граничных условий на стыке частей покрытия, описанных в [8, С. 30-31], может быть представлен в следующей формулировке [7, С. 5]:
(ах )и+ (о+о)+^ (ах К (о—о)=/к, к=. (2)
Вид дифференциальных операторов ^ (5х) и ^ определяется типом контактного взаимодействия пластин. В рассматриваемой задаче соотношения для упругой подложки имеют вид [9, С. 26]
<Ю
и(х)= |к(х1 —)Я(#)> к(х ) = — |К(а)ехр(Ча^х)¿а >
—ад 7
где и (х) и я(х ) - амплитуды перемещений и напряжений на поверхности упругого слоя соответственно; функция Грина К (а) упругого слоя с защемленной нижней гранью, содержащая в качестве параметров характеристики подложки: V - коэффициент Пуассона, р - плотность, Ц - модуль сдвига, И -толщину. Расположение в комплексной плоскости контура 7 определяется принципом предельного поглощения [10, С. 101108]. В трансформантах Фурье последнее соотношение можно записать
и (а)=к (а) о (а), ' (3)
где и (а) = V (а) и ( х ), О (а) = V (а) ё ( х ) - Фурье-образы функций и ( х ) и я (х) соответственно.
Условия сопряжения подложки с пластинами покрытия принимают вид: и+ (х )= и (х ), (х ) = Я (х )' ±х > 0. Или в образах Фурье
и(а)=и+ (а)+и_(а), о(а) = о (а)+о (а). (4)
Общие решения уравнений (1), ограниченные на бесконечности в своих областях определения, найденные с учетом принципа предельного поглощения, имеют представление
и (х) = С±1еТЧЛ + С±2е±щл + V 1 (х)[л. ' {ах)(Е±0±(а,) + В (а,))], ±х > 0, где С+ j = 1,2 - произвольные константы; > 0 - вещественные положительные корни уравнений Я+ (а) = 0;
справедливы представления
В (а ) = ехр(шхх{0 ); В (а ) = 0. Тогда для Я
Я± (а ) = £±,3 (а — д± )(а— щ± )(а + д± )(а + щ±). Применив с учетом принципа предельного поглощения к выражениям для и+ (х ) преобразование Фурье по соответствующим полупрямым, получим
и± (а)=я— [о (а)+В (а)]+К (а), (5)
w±(4)= ±ic* ■ ±lC- -
«1 ±щ± «1 ± Ч 2q± (ql+q2±)s±
g±,s G±(+q±) ig±sG±(+iff±) ax + q± ax + i q±
| B±(±q±) | L8±(+i<7±)
a ea.
ах + q± ах+ щ±
Из последнего соотношения можно выразить преобразования Фурье амплитуд контактных напряжений на нижних поверхностях пластин
а± )- Мч) + ¥± ц)
(6)
y a )=-
К ( ÜC, ^ R [G±(±q±) iG±(±iq±)
S±,5
—± 2
a ± iq± a ± q±
V "a
«1 + í/ «1 + lí/
2q± ( q2 + q2 )s± | B±(±q±) | ü?±(±ig±)
Из уравнения (1), соотношений для Фурье-образов перемещений и напряжений (3), а также условий контакта покрытия с подложкой (4) после ряда преобразований приходим к функциональному уравнению
К+ (a )Ц+ (a ) = К (a )G- (a1) + F(a1), (7)
где к (a) = ± (к (a) - е±,5R-), f = в+ (a) R-+(a)+Д. (a) R-1+(a).
В отличие от изложенного в [6] подхода предлагаемый в данной работе метод не требует преобразования оператора задачи и приводит к уравнению относительно интегральных характеристик контактных напряжений. Построенное функциональное уравнение (7) является нагруженным уравнением Винера - Хопфа в связи с наличием неизвестных G+ (±q+), G+ (±iq+), также подлежащих определению. Для этого в выражение для G+ (a)
последовательно подставляются вместо a¡ значения iq+ и q+ , а в выражение для G_ (a ) - iq_ и q . Таким образом,
получается система алгебраических уравнений относительно G+ (±q+ ) и G+ (± iq+ ), разрешив которую, остается
подставить последние в выражения для трансформант Фурье перемещений (5).
Предлагаемый в данной работе алгоритм позволяет воспользоваться интегральным методом факторизации при решении (7), так как факторизуемые функции удовлетворяют условиям теорем, сформулированных в [10, С. 113-114]. Для факторизации по a¡ относительно выбранного согласно принципу предельного поглощения контура a
функция К (a) представляется в виде К = S(a)n(a) [9, С. 61]. В качестве S используется функция S(a) = dj Jaf + d2 , где dx = lim axK(a). При этом параметр d2 выбирается достаточно большим [10, С.116-
|a
121], при численной реализации положено d2 = 10. Поведение S(a) совпадает с поведением функции K на бесконечности, на вещественной оси у S (a) нет особых точек. Особенности функции n(a) = S 1 (a) К (a)
совпадают с особенностями K, кроме того, lim n(a ) = 1 и она может быть приближена рациональной функцией
u i ™ V i /
a \-1
п(а)«п*(а) = П(а2 -21 )(а2 -р1 ) ' ' р2 - соответственно вещественные нули и полюса функции к:.
к=1
Функции К+ в (7) имеют такое же поведение на бесконечности, как и К. Их полюса - нули и полюса К. Тогда их
можно представить как к = , ^ п (а )' < = 11т а,К+ (а)'
К± п—п±(а1) 1 а.и» 1 ± V ^
у а + 1
-1 -1 N N -1
п± (а) - п* = (а2+q2) (а2 - ч1) П (а2 - ^ )П (а2 - р2 ) . Для каждого значения частоты с функции
ͱ=(a1 + q2^~2 ^ "Ь
к=\ к=1
П*, П* могут быть аппроксимированы с помощью полиномов Бернштейна [10, С. 124].
±
После подстановки найденных интегральных характеристик напряжений соотношения (5) будут содержать в качестве неизвестных только С+ ., У = 1,2, определяемые из граничных условий в области стыковки пластин (2). Для
этого к выражениям (а ) применяется обратное преобразование Фурье для определения приближенных значений
перемещений, которые следует подставить в граничные условия на стыке составляющих покрытия для получения амплитуд смещений поверхности.
Описанный подход позволяет исследовать зависимость характера прохождения сигнала от вида граничных условий в области контакта пластин, а также от свойств покрытия и основания. Результаты вычислительных экспериментов хорошо согласуются с полученными на основе подхода работы [7], что является свидетельством их
1
достоверности, в то же время очевидное преимущество нового метода - отсутствие дополнительного преобразования дифференциального оператора задачи для последующего применения метода Винера - Хопфа.
Численная реализация разработанных алгоритмов позволила проанализировать результаты для некоторых задач при варьировании физико-механических характеристик покрытия и упругой подложки. На рисунках 1,2 представлены
результаты модельных расчетов плоской задачи для безразмерной частоты О = 2, которая определялась как
(О = р+Оа2ц^ (р+ - плотность, /и+ - модуль сдвига правой пластины, а - характерный линейный размер), при
этом выбраны безразмерные значения плотностей пластин и подложки р = р = р = 1, модуля сдвига упругого слоя
р = 0,67, коэффициенты Пуассона у+ = у_ = V = 0,125, к = 2. Вертикальная сосредоточенная нагрузка приложена
в точке х0 = 5. На графиках по оси ординат отложена величина амплитуды вертикального смещения (линия красного
цвета соответствует вещественной части комплексных амплитуд, зеленого - мнимой), по оси абсцисс - координата точки поверхности системы. Расчеты проведены для трех соотношений жесткостей пластин.
Рисунок 1 иллюстрирует амплитуды смещений поверхности слоя с составным покрытием для условия жесткого сцепления его частей
^ = и .
+ \х1 = 0 - \х1 = 0
ди^
Б
д2 и,
дх2
= Б
д2 и
дхг
Б
ди
X! = 0
д3 и,
дхг
х, = О
дх3
= Б
д3 и
дх3
где Б =
кк
12 (1 -V.2)
] = 1,2.
Как видно из рис.1, с увеличение жесткости левой пластины амплитуды ее смещений уменьшаются. Тот же
результат наблюдается для других типов условий контакта.
Рис. 1 - Комплексные амплитуды колебаний при жестком сцеплении пластин
На рис. 2 приведены результаты расчетов для случая свободно смещающихся пластин, т.е. когда в точке их контакта заданы условия равенства нулю изгибающих моментов и поперечных сил
д и
дх,2
= о ди
дх3
= 0.
Конфигурация гармонического сигнала изменяется (рис.2) при переходе через начало координат, т.е. точку стыковки элементов покрытия. Смещения терпят разрыв (скачок), амплитуда смещения менее жесткой пластины больше.
2
х =0
х, =0
Рис. 2 - Комплексные амплитуды колебаний при свободном смещении пластин на стыке
Полученные результаты позволили выявить характер распространения гармонического сигнала в рассматриваемой структуре для однотипных и разнотипных элементов покрытия и различных условиях их контакта. Непосредственно на разломе определены конфигурации прохождения сигнала, которые могут служить для идентификации типа взаимодействия берегов разлома (свободное смещение, сцепление, контакт с трением и т. д.).
Развитие механико-математических методов идентификации разломов как сейсмогенерирующих зон в системе блоковой делимости земной коры на примере рассматриваемых моделей литосферных плит в виде системы контактирующих пластин Кирхгофа на поверхности деформируемой трехмерной подожки позволяет для криволинейных участков границ пластин строить кусочно-линейные аппроксимации, используя для каждого прямолинейного участка описанный в работе упрощенный подход.
Список литературы / References
1. Садовский М. А. Естественная кусковатость горной породы / М. А. Садовский // Доклады АН СССР. - 1979. -Т. 247, № 4. - С. 829-831.
2. Садовский М.А. Блоковая тектоника литосферы / М. А. Садовский, Л. И. Красный // Доклады АН СССР. - 1986.
- Т. 287, № 6. - С. 1451-1454.
3. Oparin V.N. Gutenberg velocity section of the earth and its possible geomechanical explanation. I. Zonal disintegration and the hierarchical series of geoblocks / V. N. Oparin, M. V. Kurlenya // Journal of Mining Science. - 1994. - V. 30(2). - P. 97-108.
4. Babeshko V.A. Properties of «started» earthquakes / V.A. Babeshko, O.M. Babeshko, O.V. Evdokimova / Doklady Physics. - 2016. - V. 61, Is. 4. - P. 188-191.
5. Babeshko V.A. The differential factorization method for a block structure / V.A. Babeshko // Doklady Physics. - 2009.
- V. 54, Is. 1. - P. 25-28.
6. Block element method for body, localizations and resonances / Babeshko V.A. // Экологический вестник научных центров Черноморского экономического сотрудничества. - 2014. - № 2. - С. 13-19.
7. Колесников М.Н. К методам исследования разломов в условиях вибрационных воздействий / М.Н. Колесников, И.С. Телятников // Научный журнал КубГАУ. - 2016. - № 121(07). - C. 647-659. http://ej.kubagro.ru/2016/07/pdf/33.pdf.
8. Вольмир А.С. Нелинейная динамика пластинок и оболочек / А.С. Вольмир. - М.: Наука, 1972. - 432 с.
9. Ворович И.И. Динамика массивных тел и резонансные явления в деформируемых средах / И.И. Ворович, В.А. Бабешко, О.Д. Пряхина. - М.: Научный мир, 1999. - 246 с.
10. Ворович И.И. Динамические смешанные задачи теории упругости для неклассических областей / И.И. Ворович, В.А. Бабешко. - М.: Наука, 1979. - 319 с.
Список литературы на английском языке / References in English
1. Sadovskiy M.A. Estestvennaja kuskovatost' gornoj porody [Natural lumpiness of rocks] // Doklady AS USSR [Proceedings of the USSR Academy of Sciences]. - 1979. - V. 247, № 4. - P. 829-831. [in Russian]
2. Sadovskiy M.A. Blokovaja tektonika litosfery [Block tectonics of lithosphere] / Sadovskiy M.A., Krasnyi L.I. // Doklady AS USSR [Proceedings of the USSR Academy of Sciences]. - 1986. - V. 287, № 6. P. 1451-1454. [in Russian]
3. Oparin V.N. Gutenberg velocity section of the earth and its possible geomechanical explanation. I. Zonal disintegration and the hierarchical series of geoblocks / V. N. Oparin, M. V. Kurlenya // Journal of Mining Science. - 1994. - V. 30(2). - P. 97-108.
4. Babeshko V.A. Properties of «started» earthquakes / V.A. Babeshko, O.M. Babeshko, O.V. Evdokimova / Doklady Physics. - 2016. - V. 61, Is. 4. - P. 188-191.
5. Babeshko V.A. The differential factorization method for a block structure / V.A. Babeshko // Doklady Physics. - 2009.
- V. 54, Is. 1. - P. 25-28.
6. Block element method for body, localizations and resonances / Babeshko V.A. // EHkologicheskij vestnik nauchnyh centrov CHernomorskogo ehkonomicheskogo sotrudnichestva [Ecological Bulletin of Research Centers of the Black Sea Economic Cooperation]. - 2014. - № 2. - C. 13-19.
7. Kolesnikov M.N. K metodam issledovanija razlomov v uslovijah vibracionnyh vozdejstvij [To the research methods of faults under the vibration impacts] / M.N. Kolesnikov, I.S. Telyatnikov // Scientific Journal of KubSAU. - 2016. - V. 121(07).
- P. 647-659. http://ej.kubagro.ru/2016/07/pdf/33.pdf. [in Russian]
8. Volmir A.S. Nelinejnaja dinamika plastinok i obolochek [Nonlinear Dynamics of Plates and Shells] / A.S. Volmir. -M.: Nauka, 1972. - 432 p. [in Russian]
9. Vorovich I.I. Dinamika massivnyh tel i rezonansnye javlenija v deformiruemyh sredah [Dynamic Behavior of Massive Bodies and Resonance Phenomena in Deformable Media] / I.I. Vorovich, V.A. Babeshko, O.D. Pryakhina. - M.: Nauchn. Mir, 1999. - 246 p. [in Russian]
10. Vorovich I.I. Dinamicheskie smeshannye zadachi teorii uprugosti dlja neklassicheskih oblastej [Dynamic mixed problem of elasticity theory for nonclassical domains]. / I.I. Vorovich, V.A. Babeshko. - M.: Nauka. 1979. - 319 p. [in Russian]
DOI: https://doi.org/10.23670/IRJ.2018.71.027 Шерыхалина Н.М.1, Зарипов А.А.2
1 ORCID: 0000-0002-2808-1311, Доктор технических наук, 2 ORCID: 0000-0002-6541-4680, Аспирант, 1,2 Уфимский государственный авиационный технический университет, Уфа, Россия Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (код проекта 1707-00356)
МОДЕЛИРОВАНИЕ ТЕЧЕНИЙ ЖИДКОСТИ В ПРОНИЦАЕМЫХ СЛОЯХ ГРУНТА С ПОМОЩЬЮ
РЕШЕНИЯ ЗАДАЧ ХЕЛЕ-ШОУ
Аннотация
Решается задача подземной гидродинамики для моделирования утечек воды при создании гидросооружений типа водохранилищ. Сжимаемостью жидкости и грунта пренебрегается, что позволяет в качестве модели использовать задачу Хеле-Шоу. Разработан численно-аналитический метод решения задачи с применением теории функций комплексного переменного. Рассмотрены два режима изменения давления: скачкообразное и линейное по времени с ограничением. Решены задачи с различной конфигурацией границ.
Ключевые слова: подземная гидродинамика; несжимаемая жидкость; плоская задача; комплексный потенциал.
Sherykhalina N.M.1, Zaripov A.A.2
1 ORCID: 0000-0002-2808-1311, PhD in Engineering,
2 ORCID: 0000-0002-6541-4680, PhD student, 1,2 Ufa State Aviation Technical University, Ufa, Russia This work was supported by the Russian Foundation for Basic Research (project code 17-07-00356) MODELING OF LIQUID FLOWS IN PERVIOUS GROUND LAYERS WITH USE OF HELE-SHAW SOLUTION
Abstract
The problem of underground hydrodynamics is solved for water leaks modeling during the construction of hydroconstructions of reservoir type. The compressibility of the liquid and soil is neglected, which allows us to use the Hele-Shaw problem as a model. A numerical-analytical method for solving the problem using the theory of functions of a complex variable is developed. Two modes of pressure change are considered: jumplike and linear in time with constraint. Problems with different border configurations were solved.
Keywords: underground hydrodynamics; incompressible fluid; two-dimensional problem; complex potential.
Рассмотрим плоскую задачу о течении жидкости в водопроницаемом слое грунта, обусловленном наличием втолба жидкости в резервуаре или водохранилище. До подъема уровня водохранилища проницаемый слой считается сухим, и давление в нем равно атмосферному из-за наличия близко расположенного оврага или канала (рис. 1).
Для границ с очень сложной криволинейной формой можно предложить метод декомпозиции задачи (рис. 2). Для этого область течения условно разрезается на несколько подобластей непроницаемыми пластинами. Для уменьшения погрешности, связанной со спрямлением линий тока можно использовать наложение отдельных подобластей, т. е. соседние подобласти могут пересекаться. Ширина расчетной области - L.
