Научная статья на тему 'Влияние элементов конструкции пульсирующего детонационного двигателя на тяговые характеристики'

Влияние элементов конструкции пульсирующего детонационного двигателя на тяговые характеристики Текст научной статьи по специальности «Физика»

CC BY
167
87
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ДЕТОНАЦИОННАЯ ВОЛНА / УДАРНАЯ ВОЛНА / ДЕТОНАЦИОННАЯ КАМЕРА / ДИФФУЗОР / ЭЖЕКТОР / ФРОНТ ВОЛНЫ / РЕЖИМ ЧЕПМЕНА-ЖУГЕ / A DETONATION WAVE / A SHOCK WAVE / A DETONATION CHAMBER / A DIFFUSOR / AN EJECTOR / WAVE FRONT / REGIME CHEPMENA-ZHUGE

Аннотация научной статьи по физике, автор научной работы — Поршнев Владимир Александрович

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

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

Похожие темы научных работ по физике , автор научной работы — Поршнев Владимир Александрович

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

AGENCY OF MEMBERS OF THE PULSE DETONATION ENGINE ON THRUST CHARACTERISTICS

In paper agency of members (depth of the room, a diffusor and an ejector), rendered on tractive characteristics of the pulse detonation engine (PDE) is observed. The essence of computing experiment is stated, its results are resulted. Besides the generalised results of fire tests of demonstrator PDE are resulted.

Текст научной работы на тему «Влияние элементов конструкции пульсирующего детонационного двигателя на тяговые характеристики»

УДК 629.7.036.038

ВЛИЯНИЕ ЭЛЕМЕНТОВ КОНСТРУКЦИИ ПУЛЬСИРУЮЩЕГО ДЕТОНАЦИОННОГО ДВИГАТЕЛЯ НА ТЯГОВЫЕ ХАРАКТЕРИСТИКИ

© 2009 В. А. Поршнев ОАО «Конструкторское бюро «Электроприбор»

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

Детонационная волна, ударная волна, детонационная камера, диффузор, эжектор, фронт волны, режим Чепмена-Жуге.

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

Длина цилиндрической камеры 11, диаметр ё}, длина диффузорного участка /0, угол полураствора диффузора а. Эжектор характеризуют длины /2, /3,, диаметр ё2,, угол полураствора р. В камере находится взрывчатая смесь с давлением р}, плотностью р1. Скорость смеси \1=0. Область камеры, где находится невозмущенная взрывчатая смесь, далее называется областью 1.

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

Ось ОХ направлена вдоль оси камеры, а ее начало совпадает с дном камеры. У дна камеры в момент времени t=0 инициируется детонационная волна (рис.1), которая распространяется по взрывчатой смеси в режиме Чепмена-Жуге. После того как в момент времени t¡ фронт волны детонации закончит движение по взрывчатой смеси, по воздуху от сечения х=/] пойдет ударная волна. Подвижная граница, разделяющая продукты детонации и воздух, будет представлять собой контактный разрыв, на котором плотность, температура и химический состав будут ис-

пытывать разрыв. Ударная волна, распространяясь по диффузорному участку, придет в сечения х=/]+/0, и далее будет происходить вытекание воздуха из диффузора. Когда контактный разрыв достигнет сечения х=/}+/0+/2, начнется вытекание продуктов детонации из детонационного двигателя.

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

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

Второй этап работы двигателя начинается с момента t¡ , когда в сечении х=/} возникает ударная волна Б, которая движется по воздуху и приводит его в движение, и контактный разрыв С, который разделяет продукты детонации и воздух, прошедший через ударную волну. Контактный разрыв ведет себя как поршень. Он поджимает воздух и создает ударную волну и в то же время в продуктах детонации создает волну разрежения с фронтом Б2, движущимся в отрицательном направлении оси ОХ (рис. 1,ф.

Движение продуктов детонации в камере и диффузоре можно описать системой уравнений газовой динамики идеального газа в двумерной постановке (зависимость от двух пространственных переменных и от времени) [1].

ц- ^=о

£■*=*1

{■*=&

Ь. 3 2 1 1 ^ 1

¥ 0

і і

□ і

F.

!>Н*

этап

ЭТЭП

О&ласти:

О - невозмущенный воздух Б диффузоре и эжекторе:

1- невозмушенная взрывчатая смесьвкамере:

2-простая волна разрежения в п р оду ктая дето н а ни и ;

3 — область покоя в продуктах детонации;

4-область разрежения в п р оду ктах дето н а ни и;

5- область воздуха, прош еди_его

через ударную волну.

Фронты:

□ - ФйШ: волны детонации:

задний фронт волны разрежения;

5-фшшударной волны:

Г;- фронт волны разрежения движущейся б отрицательном направлении оси ОХ;

С - .ДШШ контактного разрыва ме;нду воздухом и продуктами детонации.

Время:

1:=0- начало детонации:

t= Ъ - завершение детонации;

11=1:;. - выхода ударной волны из диффузора;

ї=1>-выхода контактного разрыва из двигателя;

ї=1:"- прихода волны разрежения всечениех=0;

при котором на выходе из диффузора скорость п р о дукто в дето Н а І1И и обратится в нуль.

Рис. 1. Режимы и этапы работы детонационного двигателя:

I - этап распространения фронта детонации по камере;

II - этап течения продуктов детонации и воздуха после завершения детонации

Выпишем эти системы для случаев плоской и осевой симметрии.

Плоская симметрия: др д (ри) д (pv)

= 0,

дt дx ду

д (pu) д (p + pu2) д (рш)

дt дx ду

д (pv) д (puv) д (p + pv2)

= 0,

= 0,

дt дx ду

д (pe) д ((є + p )u) д (( є + p ^)

- + ■

дt дx

P = (к -1) Pe,

ду

= 0,

(1)

E = є +

2

Здесь t, x, у - независимые переменные; р - плотность;

и, v - компоненты вектора скорости вдоль осей х и у соответственно;

е - внутренняя энергия единицы массы газа;

р - давление;

Е - полная энергия единицы массы газа;

к - отношение удельных теплоёмкостей.

Осевая симметрия: др + д(ри) + д(ру) = _ ру & дх дг г

д(ри) д( р + ри2) д(риу) риу

дл

дх

ду

д(pv) + д(puv) д(р + pv )

pv

дґ

дх

ду

д^) + д((е + р)и) + д((е + р)v) (е + p)v

(2)

дґ дх

р = (к -1)^

2 2 и + V

2

ду

Е = е +

Здесь t, х, г - независимые переменные.

Остальные параметры - те же, что и в плоском случае.

Граничные условия

Решения систем (1), (2) должны удовлетворять граничным условиям: на дне камеры, на стенках камеры, на стенках диффузора, на эжекторе и в выходном сечении двигателя. На поверхности детонационного

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

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

Начальные условия

Состояние продуктов детонации и воздуха в момент t=t¡ , т.е. в конце I этапа, служит начальным условием для II этапа. Начальные условия имеют следующий вид: у=У](х), р=р}(х), р=р\(х) при t=t1 , 0 < х < /1; У=У0, р=р0, р=р0 при t=t1, /1 < х < /1 + /0.

Здесь у}(х), р}(х), р}(х) - функции распределения скоростей, плотностей и давлений в трубе в момент прихода фронта детонации Б в сечение х=/}. Эти функции определяются по известному аналитическому решению при t=t1 [3]. Протяжённость волны разрежения в момент t=t1:

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

х„

в. р.

аа -

У 2 - 1 2

1 ґ1-

В момент ґ=ґі определяется распределение V на промежутке хв.р. < х < 11. В области

0 < х < хв.р. v=0. Параметры p и р определяются через V.

P = Pd

Р = Рс!

1 +

1 +

а,,

2

а.

2

У2-1

2У2

У2-1

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

величина к = 2^р~’ Рі = р0, скорость контактного разрыва ис и скорость ударной волны

и3.

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

V

а

Г

Г

2

Г

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

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

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

/1 < х < /1 +10 и на участке эжектора, который составляет с детонационной камерой единое целое.

Теоретическое значение суммарного размерного импульса давления 1С (за промежуток времени 0 < t < tk ) складывается из донного размерного импульса 1}, дополнительного размерного импульса 12 за счет диффузорного участка и дополнительного размерного импульса 13 за счет эжекторного участка:

1с =11+12 +1 з, (3)

tk ё\!2

;: 11 = I |Щ(У)(р(0, ^t) _ р0^Уё; (4)

где:

0 0

h ¡іЧ

I2 =J J (p>(xX40-peD(xX40)

0 l

X

(5)

xFPD(x, y(x)) • tga• cXct

tk h+lo

I3 =J J (pB(xy(x)’t)-Рэ(xy(x)’t))x

0 h+lo-h . (6)

xFP3 (x, y(x)) • tgb • dxdt.

Здесь FP](y) =2ny; 0 < y < d1 /2 ;

FPD(x,y(x))- периметр сечения канала на диффузорном участке;

FP^x,y(x))- периметр сечения канала на эжекторном участке в точке x;

FPd (x) = p(d1 + 2(x - ¡1) • tga) ;

¡і £ x £ /j + ¡0.

FP3 (x) = P(d2 + 2(A + ¡0 - x) • tgP) ,

¡і + ¡0 — ¡3 £ x £ + ¡0 .

рб (х> У(х) tX рвБ (х> У(х) t) - давление на внутренней и внешней стенках диффузора;

рэ (х, у( х), t), р1 (х, у( х), t) - давление на внутренней и внешней стенках эжектора.

Внутренние интегралы в (4), (5) и (6) вычисляются по методу трапеций [2], внешние интегралы вычисляются по методу левых прямоугольников [2].

Численный метод крупных частиц (метод Давыдова) решения краевых задач

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

В отчёте [3] приведено решение подобных задач при отсутствии эжектора как функции от одной пространственной переменной х и от времени 1. Этого достаточно, чтобы исследовать течение на этапе 1 аналитически, а на этапе 2 - численно с помощью метода Лакса-Вендроффа. Учёт двух измерений (х, у) по пространству в настоящем исследовании заставляет применять более сложные системы уравнений и более мощные численные методы для их решения на II этапе. К числу таких методов относится метод крупных частиц (метод Давыдова). Расчётная область разбивается на простые геометрические фигуры (прямоугольники, треугольники и т.д.), прилегающие друг к другу. Газ, попавший в каждую такую фигуру, объявляется крупной частицей. Решение получается продвижением небольшими шагами по времени. Основной идеей этого метода является идея расщепления «по физическим факторам», когда сначала «работают» силы давления, затем происходит перетекание газа из одной крупной частицы в другую, и в конце происходит объединение результатов двух первых этапов в исходные уравнения, но уже на новом временном шаге.

Две системы (1), (2) будут численно интегрироваться, используя метод Давыдова (крупных частиц) [1,4,5], модифицированный на случай произвольной треугольной сетки [6].

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

Расчётное поле с помощью алгоритма триангуляции разбивается на множество треугольников. Достоинство треугольной сетки - в её быстрой адаптации к произвольным границам.

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

Расчёт одного временного шага разбивается на три этапа. На первом, эйлеровом этапе, пренебрегают дивергентными членами. Укороченные дифференциальные уравнения Эйлерова этапа дискретизируются на треугольной сетке:

ди

P

P

P

д ґ дv дґ дЕ

+ = 0,

+

дх д р ду

= 0,

д ри д pv _ + ^— + ^— = 0.

дґ дх

ду

В частности, члены с давлением можно заменить через контурные интегралы (применяются теорема Грина и теорема о среднем [7]):

д= | рСу /1 хСу ,

дх

др

д д

ду і д

дри + дpv

[£ риСу - pvdx] /1 хСу

дх дУ а д

Контурные интегралы при численных расчётах обходятся по сторонам треугольника против часовой стрелки:

3

д Е(р}+1 + р})(У;+1 _ У;)

д р » ___________________________

д х

д р » ду~ 3

Е (хі+1 + хі)(у;+1 - у; )

І=1 3

Е (рі+1 + рі )(І - хі )

Е (хі+1 + хі)(у;+1 - у;)

і=1

д( ) 3( ) (р;+1иі+1 + ри)(уі+1 -уІ)-

д(ри) + д(pv)»І=1

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

дх

ду

3

(х,+1 + хі ) ■

і+1

і=1

- ( рі+^і+1+р^і)(хІ+1- хі)}

•(У;+1 _ У; )

Здесь: р4=рГ, х4=хГ, У4=УГ; ри р2, р3 - значения давления в трёх соседних треугольниках; (хг,уг), (х2,У2), (хз,Уз) - координаты центров этих треугольников.

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

После аппроксимации производных с давлением величины с «волнами» определяются обычным образом, как это принято в классическом методе Давыдова [4,5]:

ип = и _ 4 (д р/ д х)П,

р

У = У _ р (д р/д У)П, (7)

р

Е = Е _ р [(д( ри)1 дх)п + (д( ру)1 д У)П ].

На лагранжевом этапе учитывают то, что было опущено на эйлеровом этапе, то есть члены, описывающие перенос массы, импульса, энергии.

Выберем произвольную сторону треугольника, например, отрезок (х1; у1) _ (х2, у2).

Выпишем формулу для потока массы через выбранную границу:

ДМпК1Л = Pnv

где ДМ К

+ и

ї(норм) ¿1 2(нОрм)

2

ДґЩх^ уД(х2, у2)), (8)

поток массы через сторону

(х1, У1) _ (х2. У2). рр - плотность в текущем треугольнике, если поток вытекает из него, и плотность в соседнем через указанную сторону

треугольнике, если поток втекает;

и

і (норм)

нормальная составляющая вектора скорости

~ п

для текущего треугольника; и^ 2(норм) - нормальная составляющая вектора скорости для соседнего треугольника; А t - шаг по времени;

Щхх, УД( х2, У2)) = д/ (х _ х2)2 + (У! _ У2)2 -

длина стороны (х1, у1) _ (х2, у2).

На заключительном этапе

3

р"" = р" + Еш"і /Ї

і=1

Г р" +'£ш"В; / )/ р"

і=1

=(~; р"+’Е,АМ"~; / ї)/р"

і=1

3

Е"+1 = (Е" р" + Е ЬМ"Епр / )/ р"+1

и"+1 = (

„"+1

(9)

і=1

где - площадь текущего (і-го) треугольника.

Нижний индекс р у переменных

и.

V

р

ка,

Е"

р , равен номеру текущего треугольни-

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

На теле ставится условие непротека-

ния:

и.

и

1 гр.( норм) И вирт.(норм) •

(10)

На внешних границах ставятся обычные условия «мягкого» сопряжения (равенство нулю производных от интегрируемых функций по нормали к границе).

Таким образом, формулы (7), (8), (9), (10) представляют метод Давыдова для случая произвольной треугольной сетки.

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

На рис.2 и в табл.1 представлены результаты вычислений превышения донного давления над давлением окружающей среды во времени (0<<4). Время завершения решения краевой задачи tk. Угол а =15°; длина диффузорного участка /0=0.050 м; длина цилиндрического участка /у =0.450 м; длина эжекторного участка /3=0.400 м; угол в=15°, к2=120.

Рис. 2. Превышение донного давления Таблица 1-Превышение донного давления

Время, мс Давление, бар Время, мс Давление, бар

0.0 12.5 0.8 0.31

0.1 11.9 0.9 0.19

0.2 11.7 1.0 0.25

0.3 8.21 1.1 0.21

0.4 4.11 1.2 0.16

0.5 2.00 1.3 0.14

0.6 1.01 1.4 0.12

0.7 0.47 1.5 0.10

Интегральные импульсы (н*сек)

Рис. 3. Размерные интегральные импульсы

На рис. 3 представлены зависимости суммарного интегрального импульса 1С от времени: 1-без эжектора; 2-с эжектором; 3-отдельные компоненты суммарного интегрального импульса: 3.1 - донный импульс, 3.2 - диффу-зорный импульс, 3.3-эжекторный импульс.

Эжектор повышает отношение превышения донного давления над давлением окружающей среды на 19%, причём в начальные

"

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

На рис.4 представлены поля давления в детонационном двигателе в фиксированные моменты времени t=0•10tk, 0.254, 1.004 . Как видно из приведенных результатов, поток существенно перестраивается при изменении времени. Если в начальные моменты времени (t=0Atk, 0.254) масса газа выбрасывается через эжектор, то в последующие моменты (>0.54) воздух из окружающей среды засасывается через эжектор, что создаёт дополнительный положительный суммарный импульс 13 эжектора (для нашего варианта

расчёта 19% при t = 4).

Физические закономерности, отмеченные в результате расчётов

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

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

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

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

На рис. 4 изображены поля давления р/р0 в три момента времени ^=0.104 0.254, 1.004). Первые два момента иллюстрируют отрицательное действие эжектора, когда избыточное давление выдавливает газ наружу. На третьем рисунке газ начинает засасываться.

Рис. 4. Поля давления в детонационном двигателе в различные моменты времени

Определение влияния элементов конструкции на тяговые характеристики ПДД в процессе огневых испытаний

Огневые испытания проводились на частотах следования детонационных импульсов 20...30 Гц. В процессе испытаний исследованы несколько вариантов конструкции демонстратора ПДД:

- с цилиндрической ДК внутренним диаметром ё = 25 мм и длиной ЬдК = 320, 400, 500 мм;

- с цилиндрической ДК внутренним диаметром ё = 25 мм и длиной ЬдК = 320, 400, 500

мм, диффузором с углом раскрытия а = 30, 40, 45°;

- с цилиндрической ДК внутренним диаметром ё = 25 мм и длиной ЬдК = 320, 400, 500 мм, диффузором с углом раскрытия а = 30, 40, 45°, эжектором диаметром ёэж = 56 мм и длиной Ьэж = 85, 170, 255, 340 мм.

Всего проведено 20 серий огневых экспериментов по 3-5 в каждой серии. В статье приведены результаты исследований, проводимых на рабочей частоте следования детонационных импульсов 30 Гц при постоянном расходе для всех серий.

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

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

1 у =

P

т 0 + т г

где то + тг - массовый расход

рабочей смеси.

Кроме этого определялась скорость распространения детонационной волны по длине

ДК.

Скорость детонационной волны определена как отношение расстояния к времени прохождения ДВ между датчиками давления 014МТ и составляет Б = (1547,0 - 1739,1) м/с. (Расстояние между датчиками составляет 0,2;

0,28; 0,3 м, время прохождения измеряется при помощи цифрового осциллографа АСК-3107).

Основные результаты исследований представлены в виде диаграммы на рис. 5.

7G0

Удельный

ГТ-БОО

(Н с)/кг

5G0

40G

3G0

2G0

Л Л Л

5 5

Рис. 5. Диаграмма результатов исследований

J 5 5

Конструкция ДК

5

В пределах проведенных исследований наибольший прирост удельного импульса тяги при работе демонстратора ПДД в земных статических условиях дают увеличение длины цилиндрической части ДК до 20 калибров и длины эжектора до 6-7 калибров.

Удельный импульс тяги демонстратора ПДПД (среднее его значение на измеряемом промежутке времени) достигает необходимых в сравнении с ТРД значений кН ■ с

при условии, что длина ци-

кг/с. В этом случае

Iy = 0,6....0,8

кг

линдрической части ДК составляет 20 калибров, длина эжектора - 6,07 калибра, угол раскрытия диффузора а=30°, при частоте следования детонационных импульсов £=30Гц и суммарном массовом расходе ком-

понентов 0,0143

_ ^ кН ■ с

I = 0,646-----.

кг

Наиболее эффективным является диффузор с углом раскрытия не более 30°.

Исследования проводятся в рамках гранта МО РФ.

Библиографический список

1. Белоцерковский, О.М. Численное моделирование в механике сплошных сред / О.М. Белоцерковский. -М.: Наука. - 1984.

2. Гутер, Р.С. Элементы численного анализа и математической обработки результатов опыта / Р.С. Гутер, Б.В. Овчинский. - М.: Наука. - 1970.

3. Исследование и разработка пульсирующих детонационных двигателей для перспективных образцов вооружения и военной техники. Шифр «Почин» (II этап): науч.-техн. отчет ОАО «КБ «Электроприбор». Вып. по Гос. Контр. № 5144 от 31.05.2004. 2005 г.

4. Математическая энциклопедия. -Т.3. М.: Сов. энцикл., 1982.

5. Математический энциклопедический словарь. - М.: Сов. энцикл., 1988.

6. Шевырёв, С.П. Разностные схемы метода Давыдова на произвольной сетке. - Математика. Механика: сб. науч. трудов. / С.П. Шевырёв. - Саратов: Изд-во Сарат. ун-та, 2005. -С. 205-209.

7. Фихтенгольц, Г.М. Курс дифференциального и интегрального исчисления / Г.М. Фихтенгольц.- М.: Наука. - 1970. Т.3.

References

1. Belotserkovsky O.M. Chislennoe modelling in a mechanics of continua. М: Nauka. The main edition of the physical and mathematical literature, 1984.

2. Guter R. S, Ovchinsky B.V. Element of a numerical analysis and mathematical machining of results of experience. - M: Nauka. The main edition of the physical and mathematical literature, 1970.

3. The Technological report of «Electropri-bor» Design Bureau. «Research and working out of pulsing detonation engine for perspectiv samples of armament and the military engineering», the code «Pochin» (II stage). It is executed on the state contract № 5144 from 31.05.2004. 2005 year.

4. The mathematical encyclopaedia. Volume

3. M: the Soviet encyclopaedia, 1982.

5. The mathematical encyclopaedic dictionary. M: the Soviet encyclopaedia, 1988.

6. Shevyrev S.P. Difference of the circuit design of a method of Davidov on any grid. - Mathematics. Mechanics. The collector of proceedings. Publishing house Saratov Univercity, 2005, page 205-209.

7. Fihtengolts G.M. Kurs differential and integral calculus. Volume III. M: Nauka. The main edition of the physical and mathematical literature, 1970.

AGENCY OF MEMBERS OF THE PULSE DETONATION ENGINE ON THRUST CHARACTERISTICS

© 2009 V. A. Porshnev «Electropribor» Design Bureau, Saratov

In paper agency of members (depth of the room, a diffusor and an ejector), rendered on tractive characteristics of the pulse detonation engine (PDE) is observed. The essence of computing experiment is stated, its results are resulted. Besides the generalised results of fire tests of demonstrator PDE are resulted.

A detonation wave, a shock wave, a detonation chamber, a diffusor, an ejector, wave front, regime Chepmena-Zhuge

Информация об авторах

Поршнев Владимир Александрович, кандидат технических наук, начальник научноисследовательского центра ОАО «КБ Электроприбор», г. Саратов, 410065, 2-й Красноармейский тупик, д. 3. E-mail: [email protected]. Область научных интересов: детонация, детонационный двигатель, конструкция, рабочий процесс, статические и динамические характеристики, испытания.

Porshnev Vladimir Alexandrovich, candidate of technical science, the chief of the research centre, «Electropribor» Design Bureau, 410065, Saratov, 2-i Krasnoarmeiskii Typik, 3. E-mail: [email protected]. Area of research: Detonation, a detonation engine, design, working process, static and dynamic characteristics, tests.

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