Система алгоритмов и программа для осреднения параметров неравномерных воздушных и газовых потоков при обработке результатов испытаний гтд и его узлов Текст научной статьи по специальности «Механика и машиностроение»

МАШИНОСТРОЕНИЕ • ТЕПЛОВЫЕ, ЭЛЕКТРОРАКЕТНЫЕ ДВИГАТЕЛИ И ЭНЕРГОУСТАНОВКИ ЛА

УДК 629.7.036

В. М. КОФМАН

СИСТЕМА АЛГОРИТМОВ И ПРОГРАММА ДЛЯ ОСРЕДНЕНИЯ ПАРАМЕТРОВ НЕРАВНОМЕРНЫХ ВОЗДУШНЫХ И ГАЗОВЫХ ПОТОКОВ ПРИ ОБРАБОТКЕ РЕЗУЛЬТАТОВ ИСПЫТАНИЙ ГТД И ЕГО УЗЛОВ

Описаны разработанная система алгоритмов и программа для осреднения параметров неравномерных воздушных и газовых потоков по семи, наиболее часто применяемым в практике испытаний ГТД и его узлов, способам. Программа может входить в качестве модуля в системы автоматизированных испытаний ГТД и позволяет в зависимости от решаемой задачи использовать наиболее обоснованный способ осреднения. Программа апробирована при проведении численных экспериментов и при обработке результатов измерений в сечениях проточной части ГТД параметров неравномерных воздушных и газовых потоков. ГТД; неравномерный поток; способы осреднения

ВВЕДЕНИЕ

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

Несмотря на большое количество исследований, посвященных методам осреднения неравномерных воздушных и газовых потоков, (например, [1, 2, 3]), и установления в них целесообразности использования в тех или иных физических задачах конкретного способа осреднения, в практике работы некоторых отечественных и зарубежных двигательных ОКБ, еще довольно часто, независимо от решаемой задачи, используются различные, в том числе и упрощенные способы осреднения (например, по площади). В большинстве случаев при неравномерном поле температуры принимается допущение о постоянстве по сечению удельной теплоемкости рабочего тела. Применение упрощенных способов осреднения и принимае-

Контактная информация: (347) 273-08-44

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

На современном этапе дальнейшее совершенствование ГТД как сложной технической системы требует в процессе его доводки достоверной оценки изменения величин показателей эффективности работы узлов двигателя (коэффициентов полезного действия, восстановления полного давления и др.) в пределах 2...1% и менее. Поэтому задача осреднения параметров при обработке результатов испытаний ГТД и его узлов и при оценке показателей эффективности работы узлов является важной.

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

В настоящей статье приведено описание алгоритмов и созданной на их основе программы (комплекса подпрограмм) для осреднения

неравномерных воздушных и газовых потоков. Программа, ранее частично описанная в [4], к настоящему времени модернизирована с целью ее методического усовершенствования и выполнения осреднения по большему количеству способов. Программа составлена на языке программирования Compaq Visual Fortran Version 6.6. Разработанная программа носит как научный, так и прикладной характер. Основа разработанной программы позволяет проводить ее наращивание (путем ввода в программу других способов осреднения) или упрощение (путем отключения подпрограмм и способов, носящих научно-исследовательский характер).

2. ОПИСАНИЕ АЛГОРИТМОВ И ПРОГРАММЫ

2.1. Алгоритм расчета параметров газа в струйках и расчета интегральных характеристик неравномерного действительного потока

Переходя к описанию алгоритмов и программы и не повторяя основных положений теории осреднения, изложенных в известной литературе, следует отметить, что одним из допущений, принятых при их разработке является одномерная аппроксимация течения элементарной газовой струйки и соответствующая форма записи уравнений интегральных характеристик потока. То есть предполагается, что в сечении измерения параметров практически отсутствует закрутка потока в окружном направлении, а угол между вектором скорости потока и осью канала мал. Принимаемое допущение, как правило, соответствует реальным условиям работы узлов ГТД, так как оно применяется для выбираемых экспериментатором сечений, например, во входном и выходном сечениях узла, и предполагает, что в области течения между этими сечениями закрутка потока и угол наклона вектора скорости к оси канала могут быть значительными [5]. Разработанные алгоритмы и программа позволяют осуществлять осреднение параметров неравномерных воздушных и газовых потоков в цилиндрических и кольцевых каналах по 12 способам: по площади, по массовому расходу, с сохранением в неравномерном действительном и каноническом осредненном потоках одинаковыми трех интегральных характеристик потока: G, I, S; G, I, Ф; G, I*, S; G, I*, Ф; G, Ф, S; Ф, S, I; E, G, I; E, G, Ф; E, G, S; G, I*, E, где G -массовый расход, Ф - поток импульса, E - поток кинетической энергии, I* - поток полного теплосодержания, I - поток теплосодержания,

- поток энтропии. Осреднение по каждому из перечисленных способов осреднения производится по соответствующей подпрограмме. Учет теплофизических свойств рабочего тела осуществляется на уровне элементарных струек в области предполагаемой установки термоприемников с помощью метода и подпрограмм расчета термодинамических параметров воздуха и продуктов сгорания углеводородных топлив, разработанных Л. Н. Дружининым, Л. И. Швецом, Н. С. Малининой [6]. При расчетах предполагается, что коэффициент избытка воздуха в элементарных струйках (и соответственно, газовая постоянная) имеют одинаковое значение.

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

Статическая температура Т (при известных

* т* \

Ри р , и Т ) определяется с помощью подпрограммы ТР1 [6] из уравнения, полученного при совместном решении 1-го закона термодинамики, уравнения энтропии и уравнения состояния идеального газа для случая изоэнтропического процесса расширения

Ti dT

JCp(T) — = R ln

'T* T

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

Локальные скорости в струйках потока рассчитываются с помощью уравнения энергии по

формуле

W =42(Ai* -Ai,),

в которой приращения энтальпии определяются по формулам

а; = і срёг,

То

Ті

а/ = і срёт

То

с помощью подпрограммы БЇ [6], где Т0 - условная начальная температура газа (при х > 0

Т0 = 273 К, при х = 0 Т0 = 180 К, где х = —, а -

a

коэффициент избытка воздуха).

*

Плотность газа в струйках потока рассчитывается из уравнения состояния идеального газа

р, =

КТ

В первом блоке программы (подпрограмме OSROS) также рассчитываются параметры

Ц р(г, ф)Ж2 (г, ф)ММф т„ = Ц р(г, $)гёгё$—-----------------— ;

Л р(г, ф) Ж (г, ф)ММф

и = -

Я р(г , ф)Ж (г, ф)тйтйф Я Р(г , ф)га?гаф

где 5 = (К1 < К < Яы , ф1 < ф < фм) - область на плоскости в полярной системе координат; К = (К1,Я2,..КМ), ф = (ф1,ф2, фм) - массивы радиусов и углов, на которых осуществляет-

* гг*

ся измерение параметров р1, р1 , Т в элементарных струйках.

Параметр тг, называемый коэффициентом поля [2], характеризует неравномерность скорости потока в рассматриваемом сечении. В случае равномерного поля параметров тг = 1, неравномерного - тг > 1,0. Параметр и - среднемассовая скорость потока газа, рассчитываемая для равномерного потока по уравнению =

и = Р^ '

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

Во втором блоке программы (подпрограмма ROSGP) определяются обладающие свойством аддитивности интегральные характеристики действительного неравномерного потока:

° = Я р(г , ф)Ж (г, ф)га?га?ф,

Ф = Л (р(г, ф) + р(г, ф)Ж 2(г, ф)^^ф,

5=к л 1п

Т (г, ф)

к-1

р(г, ф)

р(г, ф) Ж (г, ф)гёгёф,

I = Лср (Т)Т(г, ф)р(г, ф)Ж(г, ф)ЫЫф ,

I* = Л ср (Т*)Т*(Я, ф)р(К, ф)Ж(К, ф)ЯdЯdф ,

ф К

Е = — Лр(г,ф)Ж3(г,ф)МЫф .

2 5

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

Т + Т*

т ___ , ,

ср ' ~ 2 '

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

, =4-1( ¥. ( к),+! + ¥. ( к)

2

\(К+1 - К),

где уДг) - подынтегральная функция при заданной .-й величине угловой координаты.

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

Фм

Ф1

.=М-1 ( 5. (ф) + ^(ф) \ )

1------- (ф.+1 - ф. ).

2

Программа позволяет выполнять операцию осреднения параметров при их измерении на 20 радиусах кольцевого или цилиндрического канала при 720 угловых положениях термопар и приемников полного и статического давления при их установке (или перемещении) с равномерным и неравномерным шагом.

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

5

К

I =1

5Т

к

бам с сохранением в осредненном потоке О, I*,

5; О, I*, Ф; О, I, 5; О, I, Ф; О, I*, Е.

2.2. Алгоритм осреднения по площади

В подпрограмме, предназначенной для осреднения параметров по площади, сначала рассчитываются средние значения давлений рр,

* Л7-Г*

рЕ и температуры торможения Тр

Pf = -1 Яp(r, 9)rdrafy ,

F S

Pf = т Я p*(r, 9)rdrafy,

F S

TF = — JJr*(r, §)rdrd§ . F ?

Далее, по известным рр, рр и Тр, с помощью подпрограммы ТР1 [6], с помощью уравнения

tf dT

JCp — = R In

'T* T

f Pf_ л *

Pf

рассчитывается средняя статическая температура Тр.

С помощью подпрограммы DI [6] рассчитываются приращения энтальпий

DiF = \CpdT,

—о

TF

DiF = J CpdT

и определяется средняя скорость

Wf =V2(DiF -DiF),

скорость звука a

a =

C„

RTf ,

где Cp = f(TF). Величина газовой постоянной Rr, зависящей от состава газа, здесь и далее определяется с помощью подпрограммы CPS [6].

Критическая скорость рассчитывается по уравнению

акр =.

C

КгТкр ,

* 1

где Ср = У(Ткр), а Ткр (при известных Т р и х = —)

а

определяется с помощью подпрограммы ТШ [6] из уравнения

* 1

i = i + —

2

С,

полученного из уравнения энергии и уравнения для определения критической скорости. Далее в данной подпрограмме по величинам осреднен-ных параметров Тр, рр, Жр рассчитываются число Мр, приведенная скорость Хр, плотность по уравнению состояния идеального газа

PF

Pf

R—f

массовый расход

Ор = р рЖрр

и другие интегральные характеристики I, I*, 5, Ф, Е канонического осредненного потока, которые при неравномерном поле параметров могут быть отличны от интегральных характеристик неравномерного действительного потока. Степень этого отличия характеризует «погрешность» определения средних параметров, обусловленную применением, в данном случае -способа осреднения по площади.

2.3. Алгоритм осреднения по массовому расходу

В подпрограмме, предназначенной для осреднения параметров по массовому расходу, сначала (при известном массовом расходе, рассчитанном ранее путем интегрирования по струйкам потока) рассчитываются средние по расходу значения давлений рО, р О и температуры Т*О

ро = Т!Яр(г, ф)р(г, ф)Ж(г, ф)гdгdф ,

О 5

ро =1Яр*(г,ф)р(г,ф)Ж(г,ф)гdгdф ,

О 5

Т* = — ЦТ*(г,ф)р(г,ф)Ж(г,ф)^^ф .

О 5

Далее порядок и формулы для расчета параметров ТО, А/ О, А/О, ЖО, МО, ХО, такие же, как и при осреднении по площади.

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

F

T

о

скорость (X = 0,2...0,72), величины потока импульса Ф и потока кинетической энергии Е в каноническом осредненном потоке практически равны соответствующим характеристикам действительного неравномерного потока. Например, для неравномерного воздушного потока, имеющего тг = 1,082 и X = 0,72, при осреднении по массовому расходу погрешности расчета вышеуказанных характеристик составляют соответственно 5Ф = 0,1% и 5Е = 0,27% [7]. При этом величина средней скорости, рассчитанная из уравнения энергии, превышает на 8,5% величину среднемассовой скорости, рассчитанную независимо от применяемого способа осреднения по приведенному выше интегральному соотношению. При расчете же средней скорости канонического потока из уравнения расхода, которая при этом меньше среднемассовой скорости на 0,9%, величины погрешностей расчета интегральных характеристик Ф и Е осредненного потока значительно возрастают и составляют 5Ф = -3,53% и 5Е = = -16,4%.

2.4. Основные положения алгоритмов осреднения по способам с сохранением в осредненном потоке О, I, I*, S, Ф, Е

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

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

—* —

О = ¥р-^ = тр ^р,

V _ Я

I = Оср (Т )Т ,

I * = Оср (Т *)Т *,

5 = ЯгО 1п

к к

Тк-1 Т к—1

= ЯтО 1п —;

р р

__ к +1 -

Ф = ОЖ + рЕ =--------Оакр г(1),

2к кр

Е = 2

2

и определяются для выбранного способа осреднения соответствующие средние параметры

р, Т, Ж или р*, Т *, 1. Далее из трех оставшихся уравнений определяются интегральные характеристики канонического осредненного потока, величины которых при данном способе осреднения не сохраняются равными интегральным характеристикам неравномерного действительного потока. Исключением является способ осреднения с сохранением Ф, 5, I, где из-за наличия четырех неизвестных в трех уравнениях предполагается, что величина О известна.

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

Т =

ср

т+т *

2

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

При осреднении по способам с сохранением

О, I*, 5 и О, I*, Ф уравнения массового расхода и потока импульса записываются, соответственно, с использованием газодинамических функций д(к) и z(X) [2]

тр д(1)Е

О =

к +1 -

Ф =------Оакр z(1).

2к кр

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

—*

жения потока Т , рассчитанной из уравнения потока полного теплосодержания.

2.5. Алгоритм осреднения по способу с сохранением в осредненном потоке О, I, S

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

По известным величинам потока энтальпии

I и массового расхода О неравномерного дейст-

вительного потока, из уравнения потока теплосодержания

I = Оср (Т )Т

методом последовательных приближений рассчитываются средняя статическая температура

Т . Далее и из уравнения потока энтропии определяется среднее статическое давление р ^

р = ^—,

У

где у

(

У = exp

S

v RG у

Величина средней скорости рассчитывается из уравнения массового расхода

Ж=.

рр

Затем рассчитываются «местные» (не сохраняемые при данном способе осреднения) интегральные характеристики канонического осредненного потока - поток импульса и поток кинетической энергии

Ф = ОЖ + рр,

Е = 2.

2

Далее производится обращение к вспомогательной подпрограмме PDT, в которой рассчитываются следующие параметры:

• скорость звука

a =

R T

где С р = / (Т),

• число м в осредненном потоке

— Ж

м = —,

а

• приращения энтальпий (с помощью подпрограммы DI [6])

__ Т

м = Ц с^т ,

-----2

—* W —

Ai =--------+ Ai,

2

и их разность

Средняя температура торможения Т для этого способа осреднения рассчитывается с помощью подпрограммы Т1 [6] методом последовательных приближений с помощью уравнения

Т *

М - ЦСр(Т)dТ = 0.

Т

Как показали расчеты, величина средней

—*

температуры торможения Т , полученная по этому способу осреднения, несколько отличается от температуры торможения, рассчитанной из уравнения потока полного теплосодержания. С помощью подпрограммы ТШ [6] по из-

—*

вестной Т и составу смеси газа рассчитывается критическая температура Ткр Ткр из уравнения

'Г* Т

Т 1 кр

Мкр = Ц Ср (Т)dТ - Ц Ср (Т)dТ =

Т0 Т0

= Ср (Ткр) КгТкр

Ср (ТКр) - Кг 2 .

Далее рассчитываются критическая и приведенная скорости

% =.

pV кр)

—RгTкp

Rг кр

- W

1 =W.

акр

Затем с помощью подпрограммы PIT [6] из уравнения изоэнтропического процесса по известным средним температурам T* и T рассчитывается соответствующее отношение средних давлений

П=^ = exp

p

TJ cp (T )dT

R J

T

Ai = Ai - Ai.

и среднее полное давление

р = П р .

Окончательно, по величине средней

—*

температуры торможения Т , рассчитывается величина «местного» потока полного теплосодержания I* (не сохраняемого при способе осреднения по О, I, 5).

2.6. Алгоритм осреднения по способу с сохранением в осредненном потоке О, I, Ф

В подпрограмме, предназначенной для осреднения параметров по способу с сохранением

О, I, Ф, сначала по известным величинам пото-

*

T

0

*

ка теплосодержания I и массового расхода О неравномерного действительного потока, из уравнения потока теплосодержания

I = Оср (Т )Т

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

оЖ 2 - ФЖ+ОКг Т = 0

рассчитывается величина средней скорости Ж .

Величина среднего статического давления рассчитывается из уравнения массового расхода

- ОКг Т р = ^~ .

Жр

Затем рассчитывается «местный» (не сохраняемый при данном способе осреднения) поток кинетической энергии. Далее производится обращение к подпрограмме PDTS для

расчета скорости звука, числа М , приращений энтальпий, их разности, средней температура

—*

торможения Т , критической скорости, приведенной скорости 1 и среднего полного давле-

— *

ния р . Окончательно в подпрограмме PDTS рассчитывается величина «местного» потока полного теплосодержания по величине средней

—*

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

5 = Кг 01п

к

Тк-1

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

Т =

ср

Т + Т

2

2.7. Алгоритм осреднения по способу с сохранением в осредненном потоке О, I*, S

В подпрограмме, предназначенной для осреднения параметров с сохранением О, I*, 5, сначала по известным величинам потока полного теплосодержания I* массового расхода О

неравномерного действительного потока из уравнения потока полного теплосодержания

I * = 0ср (Т *)Т *

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

—*

Т . Из уравнения потока энтропии определяется среднее полное давление

__*

р =

Т

к-1

где

(

У = ехр

Кг0

Из уравнения массового расхода

—* —

О = тр д(Х)р

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

Ж = Чр.

Окончательно определяются «местные» (не сохраняемые при данном способе осреднения) поток импульса, поток кинетической энергии и поток теплосодержания. «Местный» поток теплосодержания рассчитывается по величине средней статической температуры полученной при этом способе осреднения.

2.8. Алгоритм осреднения по способу с сохранением в осредненном потоке О, I*, Ф

В подпрограмме, предназначенной для осреднения параметров с сохранением О, I*, Ф, сначала по известным величинам потока полной энтальпии I* и массового расхода О неравномерного действительного потока из уравнения потока полного теплосодержания

I * = 0ср (Т *)Т *

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

—*

Т и далее величина критической скорости.

Из уравнения потока импульса

к

*

*

к +1 -

Ф =------0акр z(1)

2к кр

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

z(1) = 1 + 1/1

величины средней приведенной скорости 1 и газодинамических функций д(1), л(1), т(1).

Из уравнения массового расхода газа

тр д(Х)р

тр д(Х)р

О=

рассчитывается величина среднего полного

—*

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

Затем рассчитываются скорость звука в потоке, среднее число м и средняя скорость ос-редненного потока

Ж = Чр.

Окончательно определяются «местные» (не сохраняемые при данном способе осреднения) поток энтропии, поток кинетической энергии и поток теплосодержания.

2.9. Алгоритм осреднения по способу с сохранением в осредненном потоке О, I*, Е

В подпрограмме, предназначенной для осреднения параметров с сохранением О, I*, Е, сначала по известным величинам потока полной энтальпии I* и массового расхода О неравномерного действительного потока, из уравнения потока полного теплосодержания

I * = 0ср (Т *)Т *

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

—*

Т и далее величина критической скорости.

Из уравнения потока кинетической энергии

Е = 2

2

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

Из уравнения массового расхода газа

О =

рассчитывается величина среднего полного

—*

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

Окончательно определяются «местные» (не сохраняемые при данном способе осреднения) поток энтропии, поток теплосодержания, поток импульса.

При этом поток теплосодержания рассчитывается по величине средней статической температуре полученной при данном способе осреднения, а показатель адиабаты в уравнении «местного» потока энтропии

5 = Кг01п

определяется по величине средней температуры

Т

к-1

Т =

ср

Т + Т

2

3. АПРОБИРОВАНИЕ ПРОГРАММЫ

Апробирование программы осуществлялось при осреднении неравномерных полей параметров воздуха на входе в ГТД [7, 8, 9], на входе и выходе вентилятора ГТД [10], на выходе из камеры сгорания, на выходе из турбины ГТД. Тривиальной проверкой программы является расчет при равномерном поле давлений и температур в потоке. В этом случае, независимо от способа осреднения, величины интегральных характеристик и средних параметров потока идентичны.

С помощью разработанной программы была осуществлена описанная ниже сравнительная оценка средних величин полного давления неравномерного потока воздуха во входном сечении ГТД при осреднении параметров по площади и с сохранением массового расхода, потока теплосодержания и потока импульса (О, I, Ф). Необходимость такой оценки обусловлена следующим. Как известно, в процессе стендовых испытаний ГТД при приведении тяги к стандартным атмосферным условиям в случае неравномерного поля полного давления

к

*

*

на входе в ГТД величину измеренной тяги относят к величине полного давления воздуха во входном сечении осредненного, как правило, по площади. Кроме этого, при имитации на стенде высотно-скоростных условий при испытании ГТД с присоединенным трубопроводом на входе имитируемое число м полета рассчитывается по отношению осредненного по площади полного давления воздуха во входном сечении двигателя к статическому давлению в термобарокамере. Очевидно, что от используемого способа осреднения полного давления воздуха во входном сечении ГТД (в данном случае осреднения по площади) зависят величины приведенной тяги двигателя и имитируемого числа м полета.

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

I, Ф. Однако соответствующих количественных оценок к настоящему времени не приведено. С целью выявления влияния способа осреднения на величину среднего полного давления воздуха во входном сечении ГТД были проведены необходимые расчеты, результаты которых приводятся ниже.

Оценка отличия величин полного давления воздуха при его осреднении по площади и с сохранением массового расхода, потока теплосодержания и потока импульса осуществлялась путем сравнительных численных расчетов при различной неравномерности условного исследуемого потока (коэффициент поля тг изменялся от 1,0 до 1,45) при различном уровне заторможенной температуры воздуха (Т = = 288^500 К), поле которой было равномерным, и различном числе М на входе в двигатель (М = = 0,1^0,8). Обобщение результатов расчета производилось с использованием двух критериев: коэффициента поля тг и числа М. Результаты обобщения сравнительных расчетов приведены на рисунке, из рассмотрения которого следует, что отличие полного давления, осред-ненного по площади, от полного давления, ос-редненного с сохранением О, I, Ф, возрастает при увеличении числа М в потоке воздуха на входе в двигатель и степени неравномерности поля его скоростей. При величинах коэффициента поля тг =1,4 и М = 0,8 отличие достигает 7%. Давление, осредненное по площади, всегда меньше давления, осредненного по способу с сохранением О, I, Ф.

*

ротФ )

Рис. Зависимость отношения полного давления, осредненного по площади, к полному давлению, осредненному с сохранением О, I, Ф, от коэффициента поля тг и числа М

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

ВЫВОДЫ

Разработана система алгоритмов и программа осреднения параметров неравномерных воздушных и газовых потоков. Программа позволяет выполнять осреднение по нескольким способам, наиболее часто применяемым в практике доводки ГТД и его узлов. Программа апробирована при осреднении параметров неравномерных потоков на входе в ГТД, на входе и выходе из вентилятора ГТД, на выходе из камеры сгорания, на выходе из турбины и может быть использована как модуль в системах автоматизированных испытаний ГТД.

СПИСОК ЛИТЕРАТУРЫ

1. Седов Л. И. Методы подобия и размерности в механике. М: Наука, 1967. 428 с.

2. Абрамович Г. Н. Прикладная газовая динамика. М.: Наука, 1976. 888 с.

3. Самойлович Г. С. Гидрогазодинамика. М.: Машиностроение, 1990. 384 с.

4. Кофман В. М., Тагирова Р. А. Система алгоритмов и программ для осреднения параметров неравномерных газовых потоков при обработке результатов испытаний ГТД и его узлов // Вопросы авиационной науки и техники. Сер. «Авиационное двигателестроение»: науч.-техн. сб. ЦИАМ. 1992. Вып. 4. С. 13-17.

5. Крокко Л. Одномерное рассмотрение газовой динамики установившихся течений // Основы газовой динамики под ред. Г. Эммонса. М.: Иностр. лит., 1963. С. 64-324.

6. Дружинин Л. Н., Швец Л. И., Малинина Н. С. Метод и подпрограмма расчета термодинамических параметров воздуха и продуктов сгорания углеводородных топлив // Двигатели авиационные и газотурбинные: Руководящий техн. материал авиационной техники. 1983. 1667-83. 68 с.

7. Кофман В. М. Сравнительный анализ способов осреднения при обработке параметров неравномерного воздушного потока на входе в ГТД // Вестник УГАТУ. 2009. Т.12, № 2(31). С. 35-42.

8. Кофман В. М., Ахметов Р. Г. Экспериментальное определение потерь и неравномерности полного давления в радиально-осевом входном устройстве ВГТД // Вопросы авиационной науки и техники. Сер. «Авиационное двигателестроение»: на-уч.-техн. сб. ЦИАМ. 1992. Вып.4. С. 5-9.

9. Кофман В. М. Методы оценки показателей эффективности работы узлов ГТД по результатам измерений параметров неравномерных воздушных и

газовых потоков // Наука и технологии: Тр. XXVIII Российск. шк. М.: РАН, 2008. Т. 1. С. 108-120.

10. Кофман В. М. Определение КПД вентилятора по результатам измерения параметров неравномерных воздушных потоков // Полет. 2009. № 9. С. 38-47.

ОБ АВТОРЕ

Кофман Вячеслав Моисеевич, доц. каф. авиац. теплотехники и теплоэнергетики. Дипл. инж.-мех. по авиац. двигателям (УГАТУ, 1970). Канд. техн. наук (УГАТУ, 1983).

Иссл. в обл. идентификации и мат. моделир. ГТД.