Научная статья на тему 'Исследование закрученной струи газа, истекающей из вихревого струйного захватного устройства на плоскую преграду'

Исследование закрученной струи газа, истекающей из вихревого струйного захватного устройства на плоскую преграду Текст научной статьи по специальности «Физика»

CC BY
71
14
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
вихревое захватное устройство / струя газ / преграда / подъемная сила. / vortex gripping device / gas stream / barrier / carrying power.

Аннотация научной статьи по физике, автор научной работы — О. В. Конищева, Е. В. Брюховецкая, И. В. Кудрявцев

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по физике , автор научной работы — О. В. Конищева, Е. В. Брюховецкая, И. В. Кудрявцев

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Researches of the of Gas Expiring From the Vortex Jet Gripping Device on the Flat Barrier

The interaction of the twisted stream of gas expiring from the vortex jet gripping device on a plane barrier is considered. The developed mathematical model allows defining a pattern of a current of a flow in the device camera, and also in a zone of contact with the retained object. Rarefaction in the field of capture and influence of different parameters on its value is defined. Value of forces of suction of a stream of compressed air allowing to retain a plane disk is received.

Текст научной работы на тему «Исследование закрученной струи газа, истекающей из вихревого струйного захватного устройства на плоскую преграду»

Journal of Siberian Federal University. Engineering & Technologies, 2019, 12(2), 171-181

yflK 621.86.062

Researches of the of Gas Expiring From the Vortex Jet Gripping Device on the Flat Barrier

Olga V. Konischeva, Elena V. Bryuhoveckaya and Ilya V. Kudryavcev*

Siberian Federal University 79 Svobodny, Krasnoyarsk, 660041, Russia

Received 23.11.2016, received in revised form 01.06.2018, accepted 15.01.2019

The interaction of the twisted stream of gas expiring from the vortex jet gripping device on a plane barrier is considered. The developed mathematical model allows defining a pattern of a current of a flow in the device camera, and also in a zone of contact with the retained object. Rarefaction in the field of capture and influence of different parameters on its value is defined. Value offorces of suction of a stream of compressed air allowing to retain a plane disk is received.

Keywords: vortex gripping device, gas stream, barrier, carrying power.

Citation: Konischeva O.V., Bryuhoveckaya E.V., Kudryavcev I.V. Researches of the of gas expiring from the vortex jet gripping device on the flat barrier, J. Sib. Fed. Univ. Eng. technol., 2019, 12(2), 171-181. DOI: 10.17516/1999-494X-0126.

Исследование закрученной струи газа, истекающей из вихревого струйного захватного устройства на плоскую преграду

О.В. Конищева, Е.В. Брюховецкая, И.В. Кудрявцев

Сибирский федеральный университет Россия, 660041, Красноярск, пр. Свободный, 79

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

© Siberian Federal University. All rights reserved

This work is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License (CC BY-NC 4.0). Corresponding author E-mail address: olgakon13@mail.ru

*

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

Ключевые слова: вихревое захватное устройство, струя газ, преграда, подъемная сила.

Введение

Вихревые струйные захватные устройства (ВСЗУ), принцип действия которых основан на использовании энергии закрученной струи газа или жидкости, широко применяются в различных отраслях промышленности [1, 2]. В первую очередь, это области топочной техники и энерготехнологии.

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

В классических вихревых устройствах воздух выходит через центральное отверстие диаметром, меньшим диаметра вихревой камеры, при этом, как известно [3, 4], тангенциальная составляющая скорости в центральной части изменяется по закону вращения твердого тела (v9 /г = const). На определенном радиусе rmax тангенциальная скорость достигает своего максимального значения, а затем уменьшается по закону квазипотенциального течения (v9 rn = const), где п может принимать значения, отличные от единицы.

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

Известно множество работ по экспериментальному и теоретическому определению v9max и радиуса rmax, [3-6]. Вихревые устройства обычно имеют оптимальный диаметр выходного отверстия, при котором v9 max и разрежение на оси будут наибольшими. При уменьшении этого диаметра уменьшаются v9 max и разрежение, которое затем совсем исчезает, а при увеличении диаметра отверстия также происходит уменьшение данных величин, однако область разрежения увеличивается. При этом v9 max, уменьшаясь, сдвигается к периферии.

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

боковую кольцевую щель высотой h, как и в предлагаемой нами новой конструкции ВСЗУ [7-9].

Объекты и методики исследований

Предлагаемое устройство ВСЗУ (рис. 1) состоит из цилиндрической камеры 1 диаметром dk, плавно переходящей в развитой торец (насадок) 2, имеющий диаметр Б. В верхней части имеется кожух 3, создающий замкнутую полость вокруг камеры для подачи воздуха в четыре отверстия 4 диаметром dc.

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

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

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

В нижней части устройства радиальная составляющая V, становится соизмеримой с тангенциальной V,,, а осевая составляющая vz, наоборот, становится незначительной в том случае,

Рис. 1. Вихревое струйное захватное устройство Fig. 1. Vortex jet gripping device

если на пути истекающего из камеры воздуха устанавливается преграда, каковой и является пластина 7.

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

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

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

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

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

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

Выход воздуха ыэыекыервфериыныющель эк^е^ы^ы^С^^тэ]^^! рваличынию °тережения зо счет повышения радиальной скорости, которая вызывает эжекцию в центральной части, и градиент радиальной скоростивыIЗывaнвпoкыжeтыедэвтeлия, т.к. ecлипpенсOкыоьoceвойсостявлхю-щей скорости,то

д р ^ д Vг д

—=р—-р^ ^Т+^—г ■ (1)

r

д r r ' д r д z

Из этой зависимости видно, что давление уменьшается сильнее, если градиент радиальной скорости будет положительным, т.е. радиальная скорость будет увеличиваться при удалении от центра. Это возможно в том случае, если Vг на входе в кольцевой зазор будет равна скорости звука, тогда дальнейшее увеличениеплощади кольцевойщели вызовет увеличениескорости и падение давления, как всоплеЛаваля[10].

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

Таким образом, подъемная (притягивающая) сила F, удерживающая пластину 7, будет возникать за счет:

- разрежения в центре вихря;

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

- разрежения в области торца из-за значительного повышения составляющей скорости

v„

- сил вязкостного трения.

Течение газа, как было сказано выше, существенно трехмерно, но все параметры течения периодически (через каждые 90°) повторяются, если двигаться по окружности. В связи с этим предложено брать расчетную сетку в виде 1/4 реального устройства (рис. 2).

Течение газа в полости устройства считаем установившимся. Используем ортогональную систему координат. Ось z направим вдоль геометрической оси устройства, а оси x и y будем располагать в плоскости пластины. Искомыми функциями служат три компонента вектора скорости u, v, w, являющихся проекциями вектора скорости на оси координат x, y, z соответственно, и давление р.

Для учета турбулентности будем использовать так называемую к-е-модель турбулентности [11]. В соответствии с этой моделью к искомым функциям добавляются еще две переменные: кинетическая энергия турбулентности к и скорость диссипации турбулентности е. Газовой средой является воздух при температуре 20 °С. Наличие каких-либо массовых сил не учитывается.

В рассматриваемой нами задаче число Рейнольдса может достигать значений 105-106, что свид тельс вует о урбулентном режиме течения.

Для расчета течениявВСЗУ используем уравнения Рейнольдса [12]:

д(р") + d(pv) +d(pw)= 0 (2)

dx dy dz

Рис. 2. Расчетная область и направления осей декартовых координат Fig. 2. Settlement area and directions of axes of the Cartesian coordinates

- 175 -

d(puu} дх ) ) d(puv) ду + ¿((puw) Э z д p —— + д х Эх^ ох д х 5xZ + оу д у 5xZ uz +- dz

d(pvu) д х д у | ; рр д(рш=) dd z dp д у дтш +- д х ду +- д z

wu) d(pwv) d(pwrw) _ dp d+Twx ++Twy ddx: дх dy dz dz dx dy dz

X

wz

где p - плотность воздуха; т■ ■ т вязкие напряжения; puHи\- - турбулентные напряжения; т- -

дО 1 J б]

турбулентные потоки импульс a; i,j = 1, 2, 3, для удоб-твт в в:ы^р>^>з:;ен^я:х (6) и (2'), приняты обозначения и1 = и, и2 = v, и3 = w, x1 = x, x2 = y, x3 = z.Например,т12 = тиу. Вязкиенапряжхних обипделяютс= по сохтношсгтЕ^-н^^ю!

Tp=-puiu-,

; = ^еЮ-;

ОХ:

т ,,= Ц

fduu du: ^ —— + —-

дх,-

V J

(8)

duf

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

Турбулентные потоки импульса обычш представляют т виде прдазведенитврыииентов соответствующих осредненхых нт э Ыфептинный турТухонтнi>iaB поэффициентперев

носа:

t дЭи; du;y 2,

Tj = Ц-иГ—ь. р. —--g,

J T0 ydx{ Нх/ 3 4

(9)

к

где цeff-эффектавнаюе турбтле нтнти вяз-ость; - - кинетическая энергия турбулентности; Ъу = 0, если / да j, Ъу = 1, если / =/

Вырржение для эффективной турфю^ии-й вятдости хледующее:

(10)

где ц - коэффициент динамической вязкости; k - кинетическая эиергия турбулентности; е -скоростьдиссипации кинетической энергии; С= - эмэиричеэктя константа.

Уравнения, которые определяют кинетическую энергию турбулентности и скорость ее диссипации для стандартной &-е-модели,имеютвид:

д {ри-) + d{pvk) р—-) д

дх ду dz 5 х

( f д-Л

дх

_д_ ду

е-

y к дУ

д_ "dz

Vef- дк

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

д z

+ рр ps,

дрю) d(ve) ppwe)_ d ] V-eff де

-т —*—s + -

дх ду dz

дх

с C Тху

ду

(11)

(12)

д z

Veff д+s

е is ^ рL

тСх-роС( р — тС3 —. - к рк

Эмпирические констаитыС^О.О^а^^О, ае = 1.3, Cx = 1.43, C= = 1.92, (= = Одля стандартной k-в-модели взяtoi изра^тыЛауидертО^для модифицированы^ £ s-модееиЧеоа-Кима [13] ак = 0.8, ав = 1.15. ^ = 1.10, С2=1.9,С3 = 0.2= Скорость генерации турбулентности G такова:

Для исследования использовалось устройство (рис. 1), имеющее камеру диаметром ^ = 70 тт, диаметр сопел dc = 4 тт, длина сопел 1С = 55 тт, высота камеры Н =40 тт, диаметр насадка Б = 200 тт. Зазор к менялсяот0.1 до 2.5тт.

Давление сжатого воздуха на входе в ВСЗУ изменялось от 0.2 до 3 а! При этом расход Q изменялся от 8.5 до26.5 т3/к

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

На рис. 3-7 в качестве примера представлены только некоторые результаты, полученные теоретически и экспериментально. На рис. 3 и 4 показаны картины скоростей движения воздуха в горизонтальных сечениях зазора между нижней поверхностью устройства и плоскостью пластины 7. Расстояние 2 от пластины до рассматриваемого сечения, а также величина зазора к указаны в подрисуночной надписи.

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

У д x . '

(13)

Результаты и обсуждение

Рис. 3. Поле скоростей в сечениио = 45.8тт,е = 0.6тт Fig. 3. The field of speeds in the section at z = 45.8 mm, h = 0.6 mm

- 177 -

Рис. 5. Поле скоростей в сечении z = 0.45 mm, h = 0.6 mm Fig. 5. The field of speeds in the section at z = 0.45 mm, h = 0.6 mm

По результатам видно, что радиальная составляющая скорости мала. В области зазора (рис. 5) увеличивается радиальная составляющая скорости, которая все более возрастает при удалении от оси.

На рис. 6 представлены диаграммы зависимости давления в области зазора от радиуса r при разных величинах этого зазора h. Расход составлял 21 m3/h.

При увеличении зазора происходит рост разрежения вблизи центральной оси. При дальнейшем росте до h = 2 mm наблюдается незначительное падение разрежения, однако кривая становится более ровной, без скачков в области, где r ~ 0.5 mm.

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

На рис. 7 представлено изменение перепада давления от радиуса при зазоре h = 0.6 mm. Точками отмечены экспериментальные значения, пунктиром - расчет с помощью стандартной £-£-модели, сплошной линией - расчет с помощью модифицированной модели Чена.

Рис. 4. Поле скоростей в сечении z = 12.75 mm, h = 0.6 mm Fig. 4. The field of speedsin thesectionat z =12.75mm, h = 0.6 mm

P

100500-50-100-150-200-250-300-350-

~1—I—I--г

2___^

M

>

у

"1—I—I--1—г

0,00 0,02 0,04 0,06 0,08 0,10 0,12

Рис. 6. Зависимо стипо.маодадавления (mmH20) от радиуса (т)ариразном загора й(тт): 1 - 0.6; 2 - 1.2; 3 - 1.5; 4-Т.Т

Fig. 6. Dependence oPpreshBoediiTucepae (mmH20) ho radiur )nUa0adiffnron) pap oh(t (mm): 0 - 0.6; 2 - 1.2; 3 - 1.5; 4 - 2.0

100500-50-100-150-200-250-

/ / i / \\ V • •

/ / / / / / ' V

/ / / / /

/ / у ✓ / •

А

•У

0,00

0,02

0,04 0,0

0,08

0,10

Рис. 7. Зависимость перепада давления (mm H2O) от радиуса (m) при зазоре h = 0.6 mm (• - эксперимент; сплошная линия -модельЧеиа; пунктир - стандартная к - s модель)

Fig. 7. Dependence of pressure difference (mm H2O) on the radius (m) at h gap = 0.6 mm (• - an experiment; the continuous line - tT; iMTdel Chena; dotted ))ие - standasper-д model)

е

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

Заключение

Резукедакы исснедований показывают, что подъемная (удерживающая) сила повышается по параболическому оаипну при уве личонои давлиния сжатого воздота. а киже не-

- Л7б —

посредственно зависит от геометрических параметров, в частности зазора h. Влияние других параметров пренебрежимо мало и в работе не рассматривалось.

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

Наибольшая подъемная сила для данного устройства при давлении на входе 2.8 at и расходе 29.4 m3/h составила 96 N, в этом случае удерживаемый с грузом диск не контактировал с устройством, т.е. зазор устанавливался произвольно.

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

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

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

[1] Красов В.Б. Вихревой поток в автоматизации технологических процессов, Автоматизация технологических процессов, 1981, 6, 93-104. [Krasov V.B. Vortex stream in automation of technological processes, Automation of technologies and production, 1981, 6, 93-104 (in Russian)]

[2] Меркулов А.П. Вихревой поток и его применение в технике. М.: Машиностроение, 1997. 185 с. [Merkulov A.P. Vortex stream and its application in the machines, M.: Mashinostroenie, 1997. 185 p. (in Russian)]

[3] Гольдштик М.А. Вихревые потоки, Новосибирск: Наука, 1981. 365 с. [Goldshtik M.A. Vortex streams, Novosibirsk: Nauka, 1981. 365 p. (in Russian)]

[4] Волчков Э.П., Дворников Н.А., Терехов В.И. К расчету закрученной газовой завесы в цилиндрическом канале, Журнал прикладной механики и технической физики, 1986, 4, 59-68. [Volchkov EH. P., Dvornikov N.A., Terekhov V.I. To calculation of the twirled gas veil in the cylindrical channel, Journal of Applied Mechanics and Technical Physics, Novosibirsk: Nauka, 1986, 4, 59-68 (in Russian)]

[5] Пиролишвили Ш.А., Поляев В.М., Сергеев М.Н. Вихревой эффект. Эксперимент, теория, технические решения, М.: УНПЦ Энергомаш, 2000. 412 с. [Pirolishvili Sh.A., Polyaev V.M., Sergeev M.N. Vortex effect. Experiment, theory, technical solutions, M.: UNPC EHnergomash, 2000. 412 p. (in Russian)]

[6] Черныш Н.К. Теория и расчет идеальных вихревых устройств, М.: Медисонт, 2010. 370 с. [Chernysh N.K. Theory and calculation of ideal vortex devices, M.: Medisont, 2010. 370 p. (in Russian)]

[7] Конищева О.В., Конищев В.М. Вакуумное грузозахватное устройство / Патент 2114782 // БИ. 1998. № 19. [Konishcheva O.V., Konishchev V.M. Vacuum load gripping device, (in Russian)]

[8] Конищева О.В., Синенко Е.Г., Брюховецкая Е. . Захватные устройства для автоматической загрузки технологического оборудования, Технология машиностроения, 2012, 11, 22-24.

[Konishcheva O.V., Sinenko E.G., Bryuhoveckaya E.V. Gripping devices for automatic loading of processing equipment, 2012, 11, 22-24 (in Russian)]

[9] Конищева О.В., Брюховецкая Е.В. Расчет подъемной силы вихревого струйного захватного устройства (ВСЗУ) Современные технологии. Системный анализ. Моделирование, 2012, 4, 57-60. [Konishcheva O.V., Bryuhoveckaya E.V. Calculation of carrying power of the vortex jet gripping device (VJGD), Modern technologies. System analysis. Modeling, 2012, 4, 57-60 (in Russian)]

[10] Ландау Л.Д., Лившиц Е.М. Гидродинамика, М.: Наука, 1988. 203 с. [Landau L.D., Livshic E.M. Hydrodynamics, M.: Nauka, 1988. 203 p. (in Russian)]

[11] Launder B.E., Spalding D.B. The numerical computation of turbulent flow computer methods, Comput. Methods Appl. Mech. & Eng, 1974, 3, 269-289.

[12] Лойцянский Л.Г. Механика жидкости и газа, М.: Наука, 1987, 365 с. [Lojcyanskij L.G. Mechanics of liquid and gas, M.: Nauka, 1987, 365 p. (in Russian)]

[13] Chen Y.S. Computation of viscous incompressible flows, AIAA Paper, 1988, 88, 1241-1268.

i Надоели баннеры? Вы всегда можете отключить рекламу.