Численное исследование течения газа в тонких профилированных зазорах с вибрацией, объединенных в Т-образную структуру

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

Предлагаемый способ создания направленного течения газа [1] связан с тем, что профилирование, т. е. придание определенной формы зазору, способствует генерации в нем распределения давлений, отражающего форму профиля. Этот эффект рассматривался авторами при исследовании интегральных характеристик плоского слоя с вибрацией [2, 3], где ряд выводов были получены на основе асимптотического анализа. Однако, более полное представление об интегральных характеристиках слоя можно получить, основываясь на результатах численного решения модельных уравнений, определяющих течение газа в слое. В частности, численное моделирование позволяет получить результаты, достоверно отражающие работу слоя как в диапазоне больших частот, используемых для формирования эффекта «сдавливания слоя», так и в диапазоне малых частот. Именно в диапазоне малых частот можно ожидать переход количественных изменений характеристик слоя в качественные, например, появление расхода слоя, достаточного для промышленного использования в области приборостроения.

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

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

чивающих слой поверхностей в каждом из слоев, например, синфазное или противофазное.

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

Профилированный капилляр выполнен в элементе 1 в виде отверстия длиной L с входным Г]

и выходным г2 радиусами.

Профилированный плоский слой образован профилированной поверхностью диска 2 с глубиной профиля § и плоской поверхностью элемента 1. Поверхности элементов 1 и 2, формирующих плоский слой, расположены на расстоянии z друг от друга. Элемент 1 совершает вынужденные колебания обычно на одной из резонансных частот с амплитудой h_v. Боковые поверхности элемента 1 совершают колебания с амплитудой в k_h раз меньше амплитуды И„ где k_h - безразмерный коэффициент, близкий к коэффициенту Пуассона материала элемента 1.

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

Функции зазоров, образованных смазываемыми поверхностями элементов 1 и 2, определяются видом (1) для профилированного плоского слоя в полярной системе координат, и видом (2) для профилированного капилляра в цилиндрической системе координат:

hpi

=z + 5cos

Л" Х-Г] Я 7?-^ 2_у

+ Д,соз(т),

А« = Л + tg(/)j + ^kh К ^cos(r),

где х - радиальная координата (riу - продольная координата (0п - показатель степени; r=vr - безразмерное время; v - циклическая частота колебаний элемента 1; / - время.

В (1) второе слагаемое определяет профиль зазора, для описания которого используется тригонометрическая функция косинуса со степенной зависимостью от показателя крутизны профиля п (0<и<оо), что позволяет универсально описывать осесимметричный профиль рабочей поверхности элемента 2 любой формы. Третьи слагаемые в (1) и (2) определяют вибрационную составляющую зазора. В (2) второе слагаемое определяет вид профиля капилляра.

Для получения универсальных моделей вместо размерных величин вводятся безразмерные величины, так функция (1) относится к максимальному значению зазора (S+z), функция (2) к радиусу входного отверстия и, радиальная координата к радиусу R и продольная координата к длине L. Вводятся следующие безразмерные величины: r=x/R - безразмерная радиальная координата (к<г<1), где K=rVR; Z=z/(S+z); J=5/(<5+z); H_Ph=h_vl(6+zy, w=y!L - безразмерная продольная координата (0K=tg(y)-L/ry, H_lsr=k_h-h_vlr\.

С учетом введенных безразмерных величин функции (1) и (2) принимают вид

Н_р1 = Н_р10 + H_plT cos(r), H_ts = H_ts0 + H_tST cos(r),

TT         * n\ ' A /4- где Hpl0=Acos I y-^-

^ + Z, H_ts0

= 1 + Kw - независящие от времени составляющие функций

зазоров.

Для описания течения газа используем обычные допущения теории смазки. Полагаем, что смазочное вещество (газ) может рассматриваться как сплошная среда, динамический коэффициент вязкости не зависит от давления, режим течения в смазочном слое является ламинарным, тепловой режим изотермическим, силами инерции можно пренебречь по сравнению с силами вязкости, давление не зависит от поперечной зазору координаты, характерный продольный размер слоя много больше поперечного. В этих условиях для описания ламинарного изотермического течения газа с числами Рейнольдса порядка 7?е~1 в тонких газовых слоях с вибрацией можно использовать уравнение Рейнольдса [4], при этом математическую модель задачи получаем, решая его совместно с уравнениями неразрывности и состояния. В результате получаем уравнение, описывающее распределение давления в слое для вспомогательной функции W-PH, где Р - давление в слое. Функция ¥ обладает рядом полезных по сравнению с давлением Р свойств и, в частности, существенно большей гладкостью, так как направления изменения входящих в нее сомножителей всегда противоположны. Данный факт давно подмечен в теории смазки и использу-

Расчет и конструирование

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

Тогда, следуя работе [5], получаем следующие модельные уравнения в полярной системе координат (г, 6) для осесимметричного течения газа в плоском слое и капилляре в цилиндрической системе координат (w, в):

d dr

8 dw

tsO

2 dr ^pl dr

Я д^^{2 11} ts ^{и 1} ts

2 8w

= rAvpl

ir/2 ⁹Н₍$   _ „   .

^TtS ~       ^(sO^V

OW

д^р, 8т ’

8*Р ^{v J} ts

•is ’

8т

,     12uv/?2

ГДе Л-^-^—V

и Avls = —v ™ ■ - безразмерные частотные параметры сжимаемости профилированных плоского слоя и капилляра; д - динамическая вязкость, Ра - атмосферное давление.

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

Система граничных и начальных условий для уравнений (3) и (4) при независимом рассмотрении плоского слоя и капилляра отражает их естественный характер и может быть определена в виде

Р_р^г=к)=Н_р1(г=к), *^/(г=1)=Яд(г=1), P_pi(.T=O>=y_plo=H_plo, ^(w=0)=^Xw=0), ^=1)=Я_й(гг=1), ^(т=0)=^₀=Я_й0, ВД= Рр^ТС), Ч^У Р,_$(т+2тг).

При рассмотрении Т-образной структуры краевые условия (5) и (6) для общей точки, принадлежащей плоскому слою и капилляру с координатами г=к и w=0, заменяются условиями равенства давлений в общей точке и условием равенства расходов в плоском слое и капилляре.

Краевые задачи (3)-(7) будут использоваться для получения численного решения.

Численное решение. Для линеаризации задачи воспользуемся итерационной схемой с линеаризацией по методу Ньютона [6]

^№„)=^№)-ЦУ,).

где пространственно-временные дифференциальные операторы имеют вид

dr

r

pi

2 dr

* pl dr

- гA_{v t} ------для плоского слоя, 8т

LQP_t;) = — Я 8w

tsO

Н бТ1

11 ts и т ts

2 8w

,2 dH_ts ts OW

8*Р

- H_ts0A_vts---— - для капилляра, 8т

% - начальное приближение для искомых функций, Р_п и TVi - соответственно п и п+ 1 итерационное приближение, L' - производная Фреше в виде следующих линейных операторов:

Ah^

о r 6r

plT plO

Qr ZVplO Ypl

9Нр1 8г

Л SXFPi -rAvpl—--, 8т

L' Q?\ = — H

tsd

5^

ts* tsO OW

-2^о^

SH_ts 8w

Л ¹⁵

tsQ^^vts

ОТ

В результате, после преобразований, уравнение (8) принимает вид уравнения (9) для плоского слоя и (10) для капилляра:

дг     8т6г

8      8 „    8_

—--Е,_ =—

6w    дт6w

где по аналогии с уравнением теплопроводности введены следующие функции: - функции источника d^/o^z,,          9Hd1 HDl d^ln _ dHDl r piu pin ^Vp ^p pi___pl pin t u/^ ____Pl pl dr Pio pl" 9r 2 dr pin dr

Фр1^Н.

Фк=Н1хй

Г tsO tsn            IT/ ts ~         z'TtsO-rtsn dw

9Hls Htsd^n

dw

2 dw

,1 9H_ts tsn

OW

- функции потока

W_pl=r Я

pi

dr

^pl^pln^X ^

= Н,$о

н„

d’EtspEtsn+x               dHls

^^tsO-^tSn+X /-V OW             OW

- функции энергии Е_Р1 =гД,_р1У_р1п+1, Et$ ⁼ H_ts0 A_vls ^$„+1.

Дивергентная форма записи уравнений (9) и (10) позволяет, воспользовавшись формулой Грина, представить их в интегральной форме:

^Epldr-WpldT = -^—OpidrdT,

D                G иг с               red

Ь E_lsdw - W,_sdT = -f[— Ф,.dwdт,

J « iiQw ^ts

где D и I- замкнутые кусочно-гладкие контуры, ограничивающие на плоскостях (г, т) и (w, т) области G и J соответственно.

В соответствии с идеей интегро-интерполяционного метода [7], построим консервативные схемы на основе интегральных уравнений (17) и (18) для чего введем в рассматриваемых областях сеточное разбиение

Рис. 2. Ячейка сеточной области

^_Я={ (г„ Ту), lU;2, Ту); r^r^kr, Г^^ГркШД.Г, i= 1, 2,... N_P[+1; t₇-i= т₇+Дт,у=1, 2...}.

w_ser^ (w„ Ту), (w_1±i/2, Ту); w_/±i=Wj±Aw, w_tol/2=w,±l/2Aw, z=l, 2,... №+l; Tj^\= Tj +At,j=1, 2...}, где N_pk N_ts - количество узлов сетки по координате для плоского слоя и капилляра соответственно.

Сеточные функции источника и потока для плоского слоя и капилляра будем относить к полуцелым точкам (ц^хп, tj) и (w_;±i/2, Ту), а функции энергии, давления и зазора к целым узлам сетки. Интегральные уравнения (18) и (19)

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

уравнения (9) и (10):

/>/1+1/2

где Ж>$ —Г,&                       +(1-гг)ф>₁,₂.

а - параметр, позволяющий варь-

ировать интерполяцию сеточных функций (0<ц<1), причем при (7=0,5 получаем второй порядок аппроксимации по времени.

Путем интегрирования по отрезкам г,_]<г<г( функций (11), (13) и w,_i

Расчет и конструирование

Очевидно, что уравнения (21) и (22) линейны, так как имеющиеся нелинейности, следуя методу решения, отнесены в правую часть уравнений и определены или на предыдущем j-ом временном слое или на предыдущей итерации и и поэтому известны.

Краевые условия (5), (6) при независимом рассмотрении плоского слоя и капилляра на7+1/2 временном слое аппроксимируем следующим образом

= Н-“. ^.Л^, = Н^,(23)

(^„Д” = НЙ', (^, „.,^„ = НЙ,„,.

Начальные условия (7) принимают вид

^Д=Нр10, М-Н1$о.(25)

Разностные выражения расходов, записанных в полуцелых точках, определяются в виде

PaPa\5 + Z)        PQPj\       , .

где л / =     ----—, а₍₅ = - —-■--коэффициенты расхода, величина которых зависит от вве-

24//            31p.L денного ранее масштаба безразмерных величин, ра - плотность воздуха при давлении Ра.

Для численного моделирования совместной работы последовательно соединенных в

Т-образную структуру плоского слоя и капилляра используются ранее полученные трехточечные

Некрасов С.Г., Пашнина Н.А. Численное исследование течения газа в тонких профилированных __________________________________зазорах с вибрацией, объединенных в Т-образную структуру итерационные разностные уравнения (21) и (22). Вводится общая сетка для рассмотрения Т-образной структуры с общим количеством узлов i=N_ls+N_pi+1.

При сращении двух моделей, построенных для тонких профилированных зазоров различной конфигурации, вводится в рассмотрение общая точка (точка сращения) i=N_ts+1, принадлежащая одновременно сеточной области плоского слоя и капилляра. Особенности в схему вносятся дополнительными краевыми условиями, записанными для точки сращения. Для плоского слоя это точка с координатой г,=к, а для капилляра w_;=0. Дополнительные краевые условия

• -    условие равенства давлений P^r^^^^w^O), в точке сращения

^ = ^,где/=^+1;(28)

• -    условие равенства расходов A^X^⁼K+l/2Ar)=Ms(w;=0+l/2Aw), записанное на основе выражений массовых расходов плоского слоя и капилляра [3], с линеаризацией по методу Ньютона

^kK                 ■(29)

Интегрируя функции Ф_рЬ W_pi по отрезку r,в, W_ts по отрезку w,(z=jV_b+1) и подставляя результат в (29) с учетом (28), получим разностное уравнение, используемое в качестве дополнительного для формирования решения в точке сращения:

• =    "(^4О - <т№ + (^„ХО - ^/-и + (^„^^(1 - аН. +

^И»^+а - -*дJ- te+б - *к!

А -               Г,+1+ГД/ n _ (^о),_]     0), + (^й0^-д [

^/-1       Д^          2        К*Й1 + ^tsi-l ^*\“tSl ^tSi-ljP с =^М^А-^м „^-«л

Htsl Дг £ Kts

• + ^^ (^lli^L [_(_{я + я} )₊₂(_я __я л

L \                 —1 / \ lol        lol 1 /J

Уравнение (30) является связующим для уравнений (22) и (21) при переходе от расчета функции ^ для капилляра к функции Ч^ для плоского слоя, причем при проведении расчетов следует соблюдать условие гладкости функций, входящих в уравнение (30) в точке сращения, что обеспечивает устойчивость разностной схемы, что было подтверждено, в том числе и экспериментально. Трехточечный вид конечно-разностного уравнения (30) позволяет проводить решение методом прогонки [8] на каждом временном слое.

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

^X.n,*!АнХ,*^ ■            ■                     <³¹>

Основным критерием качества построенных схем является их сходимость, которая следует из аппроксимации и устойчивости. Погрешность аппроксимации [7] зависит от шага сетки и порядка аппроксимации, причем с увеличением порядка аппроксимации и с уменьшением шага сетки точность решения увеличивается. Для рассматриваемых дифференциальных уравнений выбран второй порядок аппроксимации по пространству 0(Ar2), 0(Aw2) и по времени О(А^), что обеспечено записью разностных уравнений с использованием центральных разностей в полуцелых точках сетки, а так же использованием схемы Кранка-Никольсона при (7=0,5. Шаги сетки по координате Ar, Aw выбираются из условия размещения 7-10 узлов сетки в краевых

Расчет и конструирование

зонах, в которых наблюдаются наибольшие градиенты давления. Величина зон определяется отношением 1/7д и.

Качественную оценку сходимости численной задачи к модельной можно получить через оценку сходимости линеаризованной численной модели к прямому конечноразностному аналогу нелинейной модельной задачи. Для этого введем в рассмотрение конечноразностную невязку по формулам

Значения невязок численных схем плоского слоя и капилляра для случая п=10, к=0,024, Z=0,67, /1=0,33, ки =0,35, ^=0,083 и H_Ph=0,5 при различных значениях частотного параметра представлены на рис. 3 и 4. Шаг сетки по времени при этом равен шагу сетки по координате. Увеличение числа узлов сетки по времени на период в 10 раз по отношению к максимальному числу узлов на изображенных здесь графиках приводит к изменению наиболее чувствительных к точности расчета величин расходов в плоском слое и капилляре не более чем на 0,0006 процента. Полученное значение сравнимо с точностью расчетов, что позволяет минимизировать время решения задачи.

Рис. 3. Значение невязок численной схемы плоского слоя: 1 -Д„=1; 2-Д,=10; 3-Д„=100

Рис. 4. Значение невязок численной схемы капилляра: 1 - Ду=1; 2 - Д„=10; 3 - Ду=100

Понятие устойчивости включает в себя условия существования, единственности и непрерывности решения при различных начальных условиях. В [7] показано, что разностные схемы, по типу аналогичные построенным нами, но с постоянными коэффициентами, являются абсолютно устойчивыми при <7=1 (неявная схема) и о=0,5. Исследования устойчивости для уравнений подобного типа с переменными коэффициентами, что и характерно для задач с переменным по координате зазором, показывают, что возможна неустойчивость при наличии зон конфузорного и диффузорного типа. Однако в данных задачах функция зазора монотонна и полученные разностные схемы также устойчивы.

Сходимость численной схемы складывается из итерационной (характеризует близость решений для и и n+1 итерационного приближения в у-й момент времени) и временной (характеризует близость решений в моменты времени у=0 и j=2it для п+1 итерационного приближения) сходимости. При реализации разностных схем итерационная сходимость оценивается косвенно где £ - малая заданная величина порядка КГ8.

Решение задач проводится на установление переходного процесса во времени, для чего необходимо проверять на каждом периоде с точностью £] > £ выполнение условий периодичности (7). Выполнение этого условия не всегда обеспечивает временную сходимость схем. Экспериментальная проверка путем увеличения времени решения задачи без ограничения точности показала, что медленной сходимости численной задачи здесь не существует. Сетка с количеством узлов 100 (по координате) и 100 (по времени на период) является достаточной для сходимости при любых значениях частотных параметров Лv.

Результаты численных расчетов и параметрическая оптимизация. Решения, полученные при реализации разностных схем (21), (22) и (30), проиллюстрируем эпюрами средних за период избыточных давлений для плоского слоя (рис. 5) и капилляра (рис. 6) в сравнении с аналогичными эпюрами, полученными при аналитическом исследовании, проведенном с использованием асимптотического анализа [3].

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

Проведем параметрическую оптимизацию рассматриваемой Т-образной структуры, исходя из критерия максимума расхода и сравнительный анализ численных и аналитических (на основе асимптотического метода [9]) характеристик плоского слоя и Т-образной структуры.

независимом рассмотрении плоского слоя и капилляра 4,5,6 - при последовательном соединении плоского слоя и капилляра в Т-образную структуру 1,4- результаты асимптотического анализа; 2, 5 - результаты численного решения для ЛУ=100; 3, 6 - результаты численного решения для Лу=1

Рассмотрим влияние на интегральные характеристики устройства параметра крутизны профиля п. Для поиска максимума расхода значение параметра п будем варьировать в диапазоне от 0 до 10 при прочих постоянных безразмерных параметах к=0,024, 7=0,67, 7=0,33, кь =0,35, ^=0,083 иИ,_/г=0,5.

Очевидно, что максимум безразмерного расхода в плоском слое при значении параметра и=0,17, полученный на основе ассимптотического анализа (кривая 1 на рис. 6, а), существенно смещен относительно аналогичных экстремумов при различных значениях частотного параметра (кривые 2-7 на рис. 6, а), полученных при численном решении. Из рис. 6, а видно, что при уменьшении частоты вибрации для обеспечения максимума расхода требуется более пологая форма профилированной поверхности плоского слоя, характеризуемая большей величиной параметра крутизны профиля, причем максимум становится не столь выраженным. Так

Расчет и конструирование

например, для A_vp^X оптимальное значение крутизны профиля и=3, что значительно отличается от результата, полученного на основе ассимптотического анализа.

О 2,5        5       7,5 И               0        2,5        5        7,5 и

• а)                                                                    б)

Рис. 6. Зависимость безразмерного расхода от параметра п: а - плоского слоя; б - Т-образной структуры. 1 - результат ассимптотического решения; результаты численного решения: 2 -Avpi™ 100; 3 Avpi= 5; 4 — Avpi= 3; 5 ■ Avpi™ 1; 6 Avpi™ 0,7; 7 — Avpi= 0,5

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

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

Иная картина наблюдается при рассмотрении Т-образной структуры. Из графиков на рис. 6, б следует, что об оптимальном значении параметра и можно говорить только при значениях частотного параметра более 100 и, кроме того, структура обладает явной чувствительностью к частоте вибрации рабочих поверхностей зазоров, причем максимального значения расхода можно добиться приблизительно при A_vpf=\.

б - Avpi=1; 1 - к= 0,024; 2 - к= 0,034; 3 - к= 0,05; 4 - к= 0,07; 5 - к= 0,1; 6 - к= 0,12

В Т-образной структуре интегральные характеристики в значительной мере зависят от параметров капилляра, например, от входного радиуса Г] или соответствующей ему безразмерной величины к, поэтому зависимости расхода от крутизны профиля п на рис. 7 отличаются от приведенных выше.

Например, на рис. 6, б показаны зависимости, полученные для наименьшего из рассматриваемых значений входного радиуса капилляра к=0,024. Но как показано на рис. 7, а для значения

Некрасов С.Г., Пашкина Н.А. Численное исследование течения газа в тонких профилированных __________________________________зазорах с вибрацией, объединенных в Т-образную структуру частотного параметра A_vpf= 100, при увеличении значения параметра к расход значительно увеличивается, причем четко определяется оптимальное значение параметра и=1, при котором достигается максимальный расход, во всем диапазоне изменения параметра к.

Иная картина наблюдается для варианта, представленного на рис. 7, б, где построены зависимости расхода от крутизны профиля для значения частотного параметра A_vpf= 1. Здесь четко просматривается зависимость расхода от входного радиуса капилляра и своего максимума расход достигает при значении к=0,07, при этом дальнейшее увеличение параметра к приводит к уменьшению расхода. На фоне явной зависимости расхода от входного радиуса капилляра наблюдается максмум расхода по параметру крутизны профиля и. Так для к=0,034 оптимальное значение крутизны профиля и=10, для к=0,05 - и=5, для к=0,07 - и=3, для к=0,12 - и=2.

Из проведенного исследования, очевидно, что для Т-образной структуры геометрические параметры значительно влияют на интегральные характеристики. В частности, следует обратить внимание на значение частотного параметра A_v в каждой конкретной задаче, поскольку от него в значительной мере будет зависеть поиск оптимальных параметров зазоров.

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

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

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