Fig. 2.13. The modification of gas generators flames form by changing air excess coefficient or.
r - rn Lfl
i
Growth of parameter _t
Fig. 2.14 The relative influence of different pwameters on diffusion flame length by secondary burning in ram/chamber of the products of propellant primary burnuig in one nozzle gas generator SPRR 181. 2-24
I ...
Fig 2.15. Multi-nozzle and equivalent one-nozzle SPRR, having the same equivalent relative combustor chamber length and the same fuel comblistion efficiency at the same a value 181. 0 one-nozzle SPRR, multi-nozzle SPRR
e
Fig 2.16 The Generalized dependence of secondary ram combustion from relative equivalence combustor length[8]. Absolute testing date: fuels - see table 3. D, = 100 - 1100 mm, Lb = 400 - 2800 (4000) mm, N,,= 1 - 228, = 0,74 - 2,5 , To = 288 - 1150 K. 2-25
1 a Fig 2.17. The complete SPRR propellant combustion effectiveness in dependence of air excess coefficient 01 by different flight Mach Number (calculation data). M3 > M2 > MI (Tot3 > Tot2 > TOtd
Fig. 2.18 The comparison of two SPRR in flight - - - Leq’ - Leq2 = 4L,,,, [8] (calculation) 2-26
Type of mixer I Open grain I 1
a Rear device
Conical device be%= m T'J - IJ Circular device - Pylon nozzle device J r? 2-1 2-1 Mixed type 91
Fig. 2.19. Possible SPRR mixing devices. [8] 2-27
b C a I I I I I 1I
Fig. 2.20 Test-bed to SI'RR with connected pipe tesring. a - test-bed inlet device. b - engine (without intakes, and supersotuc nozzle part). c - detachable "cold" calibration device. 1. - second;iry combustion chamber. 2. - gas gener;itor with propellant grain 3. - net for tiniforrnitig the flow. 4. - moving platforiii for thrtist measitremetit 5. - device for thrust measurement
mixer water
...... --.. --.
Mixers
Fig. 2.2 1 Experimental engine of SAG type with the open grain propellant with number of different mixing devices.[8]. 2-28
Fig 2.22 The tests prototype of SA6 missile SPRR with four variants of propellant nozzle devices[81.
Fig 2.23. Combustion completeness for fig 2.21, 2.22 variants [8]. Fuel CH-I (Table 3). T,, = 4503500K. ---- open end grain without mixer 3, 8 - open grain with mixers N93, and 8 00 - grain in gas generator, Nn = 12 (variant N4) I - calculation for Ceq= 21. 2-29
U
Fig 2.24 Experimental nozzle devices of pylon and mixed types.
= 16 4
1
0 1 2 3 4 5 6 Lb IDb Fig 2.25. Secondary combustion effectiveness in dependence of relative combustor length and of gas generdtor nozzle iiiimber [8] CL = 1,3; T,, = 500K. 2-30
M Fig. 2.26. Typical flight zone of middle-range air-to-air missile.
Fuei flow rate program 5
Fig. 2.27. Integral SPRR with fixed fuel flow rate programme.
!dultisectional gas-generator ------
AI H Types of trajectories tory T
Cl Trajeclorg Ill
Fuel flow rate programs
Fig. 2.28. Integral SPRR with discrete fuel flow rate control. 2-31
control pressure central body p
flow governor I
V C = A P ; 1) = 0,5-0,8 (usually 0,2-0,3) flow governor
Fig. 2.29. Integral SPRR with flexible fuel flow control.
,km
Fig. 2.30. Typical dependence of flight range of air-to-air missile on air excess coefficient at different altitude HI < H, < ... < H,. ---- optimal value of a.
3-1
SOME PROBLEMS OF SCRAMJET PROPULSION FOR AEROSPACEPLANES 0 PART I. -- SCRAMJET: AIMS AND FEATURES by Dr A.ROUDAKOV CIAM (Central Institute of Aviation Motors) 111250 MOSCOW RUSSIA
Introduction.
This lecture topic is problems of scramjet propulsion for single-stage aerospace planes. The aerospace plane main destination is acceleration of some pay load up to orbital speed. It defbed a requirement to propulsion system to be fiugal in the wide flight speed range under conditions of propulsion mass limitation and thrust compared with flight vehicle weight. In accordance with modern estimations, real hydrogen he1 scramjet specific impulse may be in several times more than that of rocket engine in the wide flight speed range, from Mf = 5...6 up to Mf = 15 ...20. But hydrogen scramjet propulsion requires sigtuficant air ram connected with high heating of vehicle structure, large hydrogen tauks 0 (and vehicle), special forms of flying vehicle. As a result scramjet advantage as well as any other air breathing engine advantage over the modem liquid rocket engines is not doubtless. Nevertheless scramjets can provide good efficiency of aerospace plane if each scramjet element and units are developed very carehlly and have high performances. In other case scramjet may lose competition with liquid rocket engines and hll scale scramjet will be never created. Theoretical and experimental investigations, ground and flight model test have demonstrated possibility of supersonic combustion ramjet for hypersonic flights. But we must cany out serious investigations to obtain high efficiency of scramjet elements and scramjet propulsion as a whole to be sure of scramjet success in aerospace plane competition with liquid rocket engine. In this lecture I shall try to repeat nothing known from books and shall try to use only last study results known only to a small number of scientists. The authors of this paper are the CIAM scientists: Dr. L.Gogish, dr. e V.Krjutchenko, dr. N.Dulepov, dr. V. Semenov, dr. Yu. Shikhman and me. Some Features of Scramiet Some years ago many scientists and Nowadays all space objects are being specialists hoped on nuclear rocket engines that may accelerated by launch systems based on chemical give propulsive efficiency 1.5... 2 times more than rocket engines. These engines' technologes are modern chemical rocket engines, but problems of enough completed and they permit to carry out some ecology and safety compelled to give up hope to use space experiments and investigations, to create some nuclear rocket engines in near future launch systems. space communication system, to use space vehicles for some unique technology, but not high propulsive At present many scientists pay their attention economies of chemical rocket engines requires to for problems of metastable substances used as rocket create multi-stage launch systems contented many one engine propulsive. These substances have energetic used components. It leads to high launch cost, to bad performances in many times better than chemical influence on ecology, to the earth surface and the rocket propulsive, but practical realization of space pollution by lost parts of rockets. metastable propulsive rocket engme has many hard problems. Solution time of these problems is not clear Modern chemical rocket engine efficiency is now. near to its limit and excludes possibility to create future effective one-stage multi-used rocket launch Now single possibility to increase sig- systems for variable aims. nificantly launch systems efficiency and to decrease
Presented at an AGARD Lecture Series on 'Researchand Development of Randscramjets and Turboramjets in Russia: December 1993 to January 1994. 3-2 cost of tespace transportation operations is air corrected and optimized continuously on base of new breathing engines using. Air breathing engines results of computational and experimental researches. 0 modern technology level, modern materials, scientific background of gasdynaniics and technology permit to The main part of thrust of any air jet engine design projects of one-stage multi-use wing launch is the difference between nozzle jet momentum and system to accelerate space objects ,to orbit and to give intake air momentum: back that to the earth surface. This system will use all aviation technology, existent runways and all ground aviation systems. These wing launch vehicles with air breathing propulsion are named the Aerospace Planes. there are R -- thrust, Aa -- air mass flow, hp -- - propulsive mass flow, V’-- flight speed, VII -- nozzle - The main advantage of air breathing engines jet speed. Pressure intakehozzle difference influence over rocket engines is atmospheric air using as is not taken into account in this formula. oxidizer and work-substance of jet nozzle. That permits to decrease consumption of on-board As a rule mp value is only several percent froni propellant in several times. Besides of air breathing ma. Nozzle jet speed is defined by energy supplied to engines using leads to economical, atmospheric flight flow due to combustion. This energy is comparable possibility. As a result air-breathing propulsion with intake air flow energy under hypersonic flight launches system is transformed to aerospace plane conditions in case of the best fuel - hydrogen - even. capable to make some maneuver in atmosphere For example ratio of combustion energy to intake air practically without additional propellant consumption. flow energy is less than 0.5 if flight Mach number is This capability will give new important quality: to more than 11.5. Therefore nozzle jet velocity is more conie to orbit of any declination practically from any than flight velocity at several percent only, if 15. point of the earth surface and in any time moment. (See Fig. 1.1.) This simple result leads to very im- This quality is practically impossible for rocket launch portant effect the sniall losses of air momentum give system due to large additional propellant significant decrease of thrust and specific impulse. consumption. Fig. 1.2 shows influences of various losses In accordance to many experts opinion air- on scramjet specific impulse. These curves were breathing technology may permit to create aerospace computed for some real losses: plane as new kind of aircraft, new kind of launch -- stognatic pressure intake ratio in accordance three system, to open new stage of space using for wide shock wave intake; aims‘ kinds in next century beginning. -- relative nozzle momentum losses are equal to 0.03; -- nozzle intake area ratio -- 3; The main success of airbreathing propulsion -- combustion perfection -- 98; technology needed for aerospace plane is connected -- relative flow losses in combustor -- .02. with supersonic ramjet concept. The scranijet is single kind of air breathing engine capable to operate at Some small variations of scramjet part hypersonic flight speed accordant to flight Mach parameters are very significant for propulsion number -- Mf> 7. Scramjets have the widest operating performances. This significance increases if flight flight speed range: dA4fz 9 in accordance to very speed increases. strict estimations and Mf > 19 in accordance to optimistic estimates. These scramjet qualities Some ideal hypothetical scramjet Is capable to determine importance of it’s application to aerospace operate up to orbital flight speed but real losses in plane in comparison with other air breathing engines. scramjets duct lead to zero thrust at hvz 23. Real scramjet specific impulse is less than that of rocket Only engines of liquid air circle may compete engme if Mf>19. The rocket engine uses more dense with scramjet in width of operating speed range but propellant and it is capable to operate in high altitude liquid air circle engines achieve the best efficiency in to decrease flight vehicle drag force. It must be taken mutual operating with scramjets. into account to compare the rocket engine and the scranijet as aerospace plane propulsion. Practically, At present idea of supersonic flow classic scramjet is not better than rocket engme at combustion engine can be admitted as realized due to Mf>15 ... 17. many ground experiments and due to two flight tests of scramjet model. But wide and serious scientific and As a rule aerospace plane propulsion concept technology researches must be carried out, many estimation is not a siniple task. Very careful problems of gasdynamics, combustion, high computational analysis and trajectory optimization are temperature materials must be soluted for creation of needed to estimate propulsion efficiency correctly scramjet to obtain necessary aerospace plane Specific impulse or some another engme parameter is propulsion efficiency. The scramjet concepts must be not enough for real estimation of engine and 3-3
propulsion efficiency. 0 A set of additional substances was selected to Hypersonic intake and nozzle are all down find out the effect of physical processes in the scramjet surface of the plane. Intake and nozzle configurations gas flow passage on its performance and to propose influence on either scramjet propulsion performances for future detailed investigations some substances, ap or outside plane aerodynamics. Nozzle/intake area plication of which may provide actual advantages and - ratio is very important for either scramjet high flight which are valid in terms of economics and ecology. Mach nuniber performances or plain transonic Passive additions were considered: inert gases performances. For these reasons nozzle/intake area (He,Ne,Ar,Kr,Xe), nitrogen (N2), water (H20). - ratio, intake and nozzle configurations must be chosen Besides the effect of 8/9 excessive hydrogen by oxygen - on base of total vehicle optimization. change was studied in order to supply additional enthalpy into the air flow. Practically all performances of scramjet - influence on the plane configuration, sizes, dry mass It is assumed, that fuel cpmponents are heat- and other parameters. It causes very hard and ed in combustion chamber cooling system and mixed important problem in design of aerospace plane and in gas generator and injected into scramjet combustion its propulsion. It may be the most significant problem chamber through gas generator nozzles tangentially to of aerospace plane/propulsion integration. air flow.
Some examples of scranijet parameters In case of oxygen application, fuel injectors influence on propulsion performances and aerospace are converted into micro rocket engines and scranijet plane efficiencies are being suggested to your is converted into combined cycle air breathing-rocket attention. engne, i.e. rocket scramjet.
Some computational investigations of addi- tional fuel component properties effect on the real Influence of fuel injection technique, fuel and scranijet main performances were carried out on the combustion properties on scramjet performances. basis of equations describing one-dimensional imper- fect gas flow (accounting heat capacity/teniperature Scranijet for aerospace planes must provide sufficient relation, equilibrium chemical reaction, dissociation thrust and efficiency in the wide range of flight speed. and recombination reactions). The scramjets perfor- A scramjet with fixed structure will have small mass mances were computed accounting some inlet and and simple technology, so it is desirable. But its thrust nozzles total pressure losses. Besides, nozzle area decreases quickly with increase of flight speed. expansion ratio was supposed restricted and constant.
For flight Mach number more than 10... 12 required thrust may be provided with significant fuel Statement of the problem excess. Excessive part of hydrogen doesn't release chemical energy. It only increases propulsive mass Pure hydrogen (monopropellant) scramjet through the nozzle. Hydrogen is best fuel for thrust efficiency performances were compared with ' aerospace plane engines, due to its high heat capacity that of bipropellant scramjet, provided that hydrogen and cooling capacity. But usage of hydrogen as consumption in the first engine (MI)) was equal to the passive propulsive mass doesn't seem valid, because of sum of hydrogen consumption (mp) and additional very low density of liquid hydrogen that results in component consumption (rnef2)) in the second engme: necessary large tank. One may suggest to substitute it by any non reactive substance.
Really, physical properties of additional All other conditions were the same, i.e. the substance may influence on gas' expansion processes same scramjet geometry, flight conditions, efficiency in the nozzle and may lead to changes in the engine of engne units and so on. These conditions lead to performances. equality of air flow rates through both scramjets:
Furthermore, the scranijet performance can be improved by the fuel impulse utilization. The fuel will have high pressure and good enthalpy if it will be All the results were obtained on condition the used after pump for the cooling of gas flow passage. hydrogen part, reacting with air, was the same in Therefore injector fuel impulse may be compared with every case. That is, for usage of any additional non difference of exhaust and inlet impulses. The value of reactive substance, the reacting part of hydrogen was fuel impulse depends on fuel substance properties kept the same also. 3-4
correspond to boiling temperature values at the see- level pressure with the exception of water ; 0 where La,+= 34.3 -- stoichiometric mixture - nozzle exhaust area -- inlet throat area ratio ratio of air-hydrogen reaction. Fe = Fe/Ft = 30; - nozzle impulse factor Ie = 0.97. lntroducing additional fuel component mass - fuel components are heated in combustion flow ratio chamber cooling system to total temperature 1000 K, - mixed in gas generator chamber and injected into combustion chamber tangentially to air flow through gas generator nozzles; expression (1) may be written - for hydrogen-oxygen fuel components gas generator is a liquid-hydrogen-oxygen micro rocket engme and scramjet is a rocket scramjet engine; - gas generator nozzle exhaust area -- throat .' for all the considered additional substances area ratio Fge = Fge/Fgt = IO, but for oxygen. ( B,+)= BW= BO)= B -- total fuel mass -gas generator nozzle impulse factor Ige- flow ratio for the first and for the second engines. ) 0.98.
Additional oxygen consumption was selected Fig. 1.4 shows the effects of relative addi- such that hydrogen-air-oxygen mixture was sustained tional component consumption on gas generator ex- stoichiometric: haust specific impulse Ig for different components. The value d is defined as
where heoQ) -- oxygen consumption, LOH = 7.94 stoichiometric mixture ratio of oxygen-hydrogen The relationships between d and B can be react ion. obtained from previous forniulas for non reactive Expression (4) may bewritten: component
d = I - I/B
or and for oxygen
d = (I - l/B)/(l +- I/Lo)
Total fuel (propellant) mass consumption in For non reactive components the distinction the bipropellant scramjet was always kept equal to the between impulses for different components and im- hydrogen consumption in the monopropellant pulse decrease with increase of impulse value are scranijet (see eq. 1). Therefore: explained by effect of gas mixture molecule mass on nozzle exhaust impulse fg. It is known, that for fixed gas total temperature nozzle exhaust impulse depends mainly on gas molecule mass (Ig-l/(p)/'2). For this Computations of scramjet performances were reason the substitution of excessive part of hydrogen carried out with the aid of the program, computation by gas with the molecule mass higher than that of hy- procedure of which is based on the solution of one- drogen (see table 1) leads to increase of mixture dimensional gas dynamics' equations, accounting molecule mass and, consequently, to decrease of gas equilibrium chemical reactions, air and combustion generator impulse. gases dissociation. The computations were carried out in the flight Mach number range Ad := 8...20 with the The highest impulse for hydrogen-oxygen following assumptions: components is explained by energy release in chemi- - flight path corresponds to the constant air cal reaction between hydrogen and oxygen in gas ram value q = 75 kPa; generator combustion chamber. - air inlet characteristics are shown in Fig. 1.3 as flight Mach number dependent relationships for Hydrogen-oxygen gas generator nozzle inlet total pressure recovery factor Pi and mass flow exhaust specific impulse is presented in Fig. 1.5 for factor Fi = Fa/Ft (where Fa -- captured air flow area, different hydrogen total temperature Ttl. The value /3 Ff -- inlet throat area); is gas generator equivalence fuel ratio - the combustion chamber is expanding, expansion ratio FdFf = 1.2; - fuel component tank temperature values 3-5
Note, that in case of TH21000 K pure 9...16. Specific impulses corresponding to Ar,Xe are hydrogen gas generator impulse is higher than that of 0.1... 0.3 kN s/kg lower than that corresponding to Ne, the hydrogen - oxygen gas generator and in case of that makes 2...3%. But density impulse corresponding TH= 500 K vice versa. to He is 0.1... 0.5 MN s/cub.m (20 ...30%) lower than that, corresponding to Ne, inspite of the higher values Fig. 1.6 shows effect of gas generator impulse of specific impulse. Zg on pure hydrogen fuel scramjet specific impulse Zm. Dashed lines correspond to normal frrel injection (Is= That is explained by the fact that Ne density 0). and continuous lines correspond to tangential fuel is one order higher and that results in 1.7 times higher injection (lg 0). One can see, that fuel impulse has an density of bipropellant fuel with Ne compared to fuel essential influence on scramjet impulse. For example, with He. Values of density impulse for Ne,Ar,Kr and at M = 12 and B = 2.5 the increase of scramjet Xe nearly coincide. Density impulse curves for Ar and impulse is about 1.2 kNs/kg (20%). Kr lie between the curves for Ne and Xe and are omitted in Fig. 1.12. The marked features are typical Scramjet mass and density specific impulses, for bipropellant fuels with different percentage of inert obtained for neon application as additional fuel coin- gases. When total fuel flow ratio (13) increases, ponent with total bipropellant fuel consumption ratio absolute values of increments lm and lv increase too, B = 1.43; 2.5; 5, are shown in Fig. 1.7 and Fig. 1.8. and vice versa. For comparison pure hydrogen fuel scramjet impulses are given. One can see, that specific impulse (Fig. 1.7) Specific impulse - fuel composition relation- for bipropellant fuel with neon is somewhat lower in ship, chemical energy supply being constant and all the considered range of flight Mach number (M = the engme parameters being identical, is explained by 7...20). Density impulse (Fig. 1.8) is significantly properties of real combustion gases. higher, if compared with that for pure hydrogen fuel. So at M = 12 for pure hydrogen engine and B = 2.5, In reality combustion gases heat capacity specific impulse lg = 8.2 kNs/kg and density impulse increases with temperature growth and the combustion lv = 1.18 MN s/cub.m, i.e. replacement of hydrogen gases are partly dissociated. It particularly regards to excessive part with neon results in specific impulse water vapor and nitrogen. reduction by 12% and double density impulse increase. Fig. 1.13 shows relative specific heat capasity of N, H, 0, Ar - temperature curves. Relative specific Engine thrust changes exactly like specific heat of each gas is the ratio of its specific heat and impulse because comparison of mono- and that at temperature 400 K. One can see significant bipropellant fuels is made at the equal fuel flow rates. heat capacity increase of nitrogen and particularly of That is why thrust versus M relationships are omitted water with temperature raise. Argon heat capacity is here and hereinafter. constant in wide range of temperatures, argon (as well as other inert gases) is practically ideal gas. Replacement of hydrogen excessive part by Inconstancy of heat capacity results in lower kinetic helium (Fig. 1.9 and Fig. 1.10) changes scramjet per- energy increase of real gases molecules compared to formances less significantly. Specific impulses of that of ideal gases, in case of heating. Some kinetic hydrogen-helium and pure hydrogen scramjets practi- energy increase delay occurs, part of lunetic energy is cally coincide. Density impulse increase is not so realized in the nozzle at lower pressure. As a result significant due to the comparable densities of liquid effectiveness of engme thermodynamics cycle de- helium and hydrogen. grades.
If fuels with the other inert gases (Ar.Kr,Xe) Another significant feature is nozzle perfor- are used, scramjet performances are close to that of mance dependence (for fixed nozzle expansion ratio) the hydrogen-neon scramjet. on adiabatic index of combustion gas mixture. The higher adiabatic index - the higher the conversion In Fig. 1.11 and Fig. 1.12 engine perfor- factor of total gas enthalpy into exhaust jet kinetic mances for bipropellant fuel are coinpared with that energy in the nozzle. Inert inonoatomic gases are for monopropellant hydrogen fuel (total fuel flow ratio practically ideal gases, their heat capacities are prac- being 13 r. 2.5). Increments of specific impulse I = tically constant, they can not dissociate, their adiabatic I,, -- lll,rrand of density impulse I, = I, -- lyH, indexes are the highest ( = l.67). corresponding the considered gas, are shown. At last, scramjet performances depend on gas It can be seen that values of specific impulse, generator exhaust impulse (Fig. 1.6). the value of that corresponding to helium, are approximately 1.0 kN depends on propellant composition ( Fig. 1.4). sfig higher than that corresponding to neon, that makes 10... 30% in the Mach number 'range M = In case of excessive hydrogen replacement 3-6 with inert monoatomic gas, gas mixture in the nozzle propellant injection. One can see, that for the case of beconies closer to ideal gas, its adiabatic index in- tangential injection hydrogen-oxygen scramjet has the creases and one should expect specific impulse in- highest values of density impulse for the all crease. But lower inert gas heat capacity compared considered Mach numbers and the highest values of with that of hydrogen (especially for Ne and heavier specific impulse in Mach number range 9... IG at B . gases) results in higher gas mixture temperatures. 2.5. Note, that by normal injection the highest value of density impulse corresponds to neon-hydrogen fuel In consequence, non-ideality effect of water scramjet and specific impulse of hydrogen-oxygen and nitrogen, forming essential part in the nozzle scranijet is higher than that of hydrogen scranijet only exhaust flow, becomes stronger. for A4 = 9.5.
The effect of these two factors (adiabatic The specific impulse increase of hydrogen- index and heat capacity) is shown in Fig. 1.14 and oxygen scramjet is achieved by additional chemical Fig. 1.15, obtained at the condition of normal fuel energy, that is provided by hydrogen and oxygen. But injection (lg = U) into combustion chamber. One can specific impulse increases slightly, because the part of see, that the effect of the second factor is stronger than water in the nozzle gas mixture increases and tem- that of the first factor (except the case of part helium perature raise shows non-ideal properties of nitrogen application and B>2.5 ). and water. In the case of tangential fuel injection the effect of non-ideality appears lower than in case of The effect of the third factor (gas generator normal fuel injection, because tangential fuel injection impulse) increases the distinction between the leads to gas velocity increase (gas temperature impulses of nionopropellant hydrogen fuel scramjet decrease) in combustion chamber and, in conse- and of bipropellant fuel scramjet (compare Fig. I. 14 quence, the larger part of additional chemical energy and Fig. I. 11) because of the gas generator impulse is converted into exhaust jet lunetic energy. In Fig. for hydrogen fuel is higher than that for inert gas and 1.20 and Fig. 1.21 specific and density impulses of hydrogen mixture (see Fig. 1.4). hydrogen-oxygen scramjet at tangential fuel injection are compared with that of pure hydrogen scramjet for In consequence, in case of tangentially different values of R = 1.43; 2.5; 5. Note, specific injected fuel, specific impulse of bipropellant scramjet impulse of hydrogen-oxygen scranijet is higher then appears lower than that of pure hydrogen scranijet for that of hydrogen scramjet at B = 5 and hf = 6...20. any inert gases as fuel additio. But dencity inipulse of bipropellant scramjet is higher than that of hydrogen Fig. 1.22 differs from Fig. 1.20 only by scranijet for tangential (Fig. I. 12) and nornial (Fig. values of gas generator exhaust impulse /g. In Fig. 1.15) fuel injection (especially for Ne and heavier 1.20 they correspond to the hydrogen temperature at gases, Fig. I. 12). gas generator chamber inlet T = /UUUK,while in Fig. 1.22 they correspond to 7’ = 2OUUK (see Fig. 1.5). One This effect, in spite of its specific impulse can see, that inspite of higher gas generator impulse reduction, may be very useful in practice, because it for pure hydrogen scramjet (lg = 6.9 kNs/kg) specific permits to reduce the needed volume and dry mass of impulse of hydrogen-oxygen scramjet is higher than tanks. that of hydrogen scramjet in the wide flight Mach number range (for B = 2.5 at A4 = a...12). Effect of real-gas-mixture properties influ- ence on the engine performances is observed also Previously mentioned effects of gas mixture when non-monoatomic substances are used as addi- reality influence on exhaust velocity from nozzle are tional components of the fuel. Fig. 1.16 and Fig. I,17 small, but essential for scranijets at high speeds. It is (like Fig. 1.11 and Fig. 1.12) for the case of tan- known, that in high flight speeds, thrust of any air- gentially injected fuel and Fig. 1.18 and Fig. 1.19 breathing engme is defined by the small difference of (like Fig. 1.14 and Fig. 1.15) for the case of normal two large values: exhaust and inlet impulses. Even fuel injection show the results, when water, nitrogen slight change in exhaust velocity may result in sig- and oxygen were used as additional components. The nificant change of thrust and specific impulse of the results of water and nitrogen application are worse engine. It complicates the quality of investigations of compared to that of neon in the both cases. It can be scranijet for aerospace planes and requires as detailed explained by non-ideal properties of water and nitro- computations as possible. gen: inconstancy of heat capacity and low value of adiabatic index. Carried out computations could not take into account possible effect of the additional substance on Replacement of excessive hydrogen part by spatial effects in the engine passage. Besides, prob- oxygen and hydrogen in stoichiometric ratio is sig- lems of engme construction cooling in case of nificantly inore effective in the case of tangential gas bipropellant fuel application were not considered. generator gas injection than in the case of normal fuel Accounting of these effects requires more complicated 3-7
model of processes in a scramjet. As for aerospace -- inaxinnuin cross acceleration - 1.5g planes, to estimate the engme effectiveness would -- maximum trajectory angle 35 degrees. require detailed computations of the plane trajectory, particularly, tank volume effects on the aerodynamics The optimization of ascent trajectory was must be accounted. Due to these reasqns one cannot conducted with account to these and some other conclude that application of additional substances limitations on different phases of motion. The
~ instead of excessive hydrogen is absolutely advanta- aerodynamic scheme of hypothetical ASP can be geous. No doubt, the problem is worth carrying out defined as wing body without horizontal tail. The more deep investigations. bottom nose part of the fuselage forms engme compression surface, and fuselage afterbody is used Some conclusions for nozzle gases expansion.
Tangential injection of fuel gives significant For efficiency comparison of various PS - improvement of scramjet performances and must be unified model was accepted. The following mass- and applied on aerospace plane scramjet propulsion. volumetric characteristics were assumed:
Some additional fuel component effects on Mnss Charcteristics hydrogen scramjet performances show significance of small differences of propellant and combustion prod- a) Propulsion system: @ ucts properties -- specific mass of thrust loaded structure panels, It seems attractive to use oxygen as addi- 16 kg/sq. m; tional fuel component; oxygen's application allows to -- specific mass of SCRJ (SCRJ-R) cooled panels, increase as density impulse (1.5 ...2 times), as well as 20 kg/sq.m; specific impulse (by 3...10%) of engme and to reduce -- specific mass of air entry and nozzle panels, the needed tank volume and tank dry mass. 20 kg/sq.m; -- specific mass of LR, 0.0009;kg/N; - Some inert gases application is interesting -- specific mass of ATR, 0.004 k@N. also. More sure results of scrannjets modifications b) Airframe: comparison may be received by the aerospace plane trajectory optimization methods. -- fuselage structure specific mass, 23 kgkq.m; -- wing structure specific mass, 30 kg/jsq.m; -- fin structure specific mass, 35 kg/sq.m; -- hydrogen tanks specific mass, 18 kg/sq.m; Trajectory Efficiency Analysis of Oxygen Boosted -- oxygen tanks relative mass, 0.015; Scramjet. -- constant nnass (crew, equipment, undercarriage, etc.), 19 t. Computational efficiency analysis of single Geometrical Charnclerisics. a stage aerospace planes (ASP) was carried out to compare two propulsion system (PS) concepts: I-- wing span to length ratio, 0.44; -- the first PS uses classical scramjet (SCRJ) -- basic wing leading edge sweep angle, 52 deg.; operating only on hydrogen fuel; -- strake leading edge sweep angle, 83 deg.; -- the second PS uses on-board oxygen boosted -- wing profile thickness, 0.05; scramjet (SCRJ-R). -- wing fineness ratio, 1.83;
Efficiency analysis of two PS concepts -- on The wing loading and disposable fiiel basis of SCRJ and SCRJ-R -- was conducted for an capacity were determined by solving two equations: example of hypothetical single-stage ground starting ASP with start mass of 280 t. The ASP is designed to -- the equation of vehicle existence (mass balance), launch the payload on near-earth circular orbit ( H = -- the equation of matching the required and 200 kni ) with inclination angle of 51 deg. The disposable vehicle volume. latitude of starting place -- 47 deg. The additional characteristical velocity needed to orbit maneuver and Airframe modelling . was assumed for re-entry starting was taken equal to 120 ids. The comparative evaluation of alternative PS versions and ascent trajectory was determined with following design parameters. In accordance with this airframe constraints on vehicle motion parameters: modelling, the absolute sizes of vehicle (length, wing -- maximum dynamic pressure allowed - 75 kPa, surface, etc.) are changing in some limits while -- niaxiinuni longtudinal acceleration - 3g, relative parannaters (span to length ratio, relative 3-8
maximum fuselage area, sweep angle, aspect ratio, flight Mach number are shown too. etc.) remain unchanged. Thus the aerodynamic characteristics also remain constant. The results of aerodynamic characteristics computations were confirmed by flight experiments of Bouran Space Vehicle program and hypersonic Aerodynamic Characteristics of scramjet test flights. Aerospaceplane (ASP)
Aerodynamic characteristics of the ASP Propulsion System Concepts were determined by mathematical modelling method. It is known that the possibility to obtain experimental Two concepts of propulsion system are- data for hypersonic velocities range is limited by considered: aerodynamic facilities capacity. At the same time the CFD (coniputational fluid dynamics) methods used for - the first PS based on scramjet; large computational investigations connected with - the second PS based on scramjet-rocket. optimization of ASP propulsion system (PS) require They are identified as: SCRJ concept and very great computer resources. Therefore the SCRJ-R concept. engneering method published in was used for evalua- tion of aerodynamic forces and moments acting on These combined PS concepts consist on vehicle at hypersonic flight speeds. following engine types:
The friction forces, heats flux and equi- -- In the first acceleration phase (MG) air-turbo- librium wall temperature are defined on base of some rocket engine (ATR) o'perating on oxygen-hydrogen boundary layer model. This model permit to deterni fuelk used. gas parameters on viscous boundary layer on base of -- In the second phase from M=5 SCRJ or SCRJ-R Newton and Prandtl-Mayer theories, theory of tangent are used. cones and wedges and takes into account entropy effect and hypersonic viscous interaction. The -- In the third (final) phase of acceleration oxygen- corresponded algorithms and computer codes were hydrogen Liquid Rocket Engme (LE) is operating described in. with space specific impulse equal to 4576 m/s and oxygenhydrogen flow ratio equal to 7. The method based on supersonic potential theory, on slender body theory and on some empirical SCFU/LRE or (SCRJ-R)/LRE transition dependencies was used to determine approximetly flight Mach number were optimazed for highest profile drag force, polar slope at subsonic, transsonic payload mass. and moderate hypersonic flight speeds.
In Fig. 1.23... 1.26 are shown the confgu- Optimum Propulsion Operating Modes of ration and aerodynamic characteristics of ASP SCRJ and SCRJ-R considered in this paper are shown as drag force factor and lift force factor dependencies on Mach number The values of ATR effective specific and angle of attack. Maximum aerodynamic fineness impulse along the ascent trajectory dependmg . on (Kmax) and corresponding angle of attack versus Mach number and flight altitude are lisied in Table 1.
A TR Performances Table 1.
H, kni 0.00 0.00 0.33 11.80 15.50 18.50 21.60 Jsp, knds 20.59 20.1.2 21.75 21.65 21.30 20.90 21.80 T, K 288 295 380 600 880 1075 1250
optimized for all operating part of trajectiry, bu The operating modes and flight speed range limitation of SCRJ combustor choking requied to on diverse hases of erospaceplane acceleration clidis air excess at Mf < 7 without optimisation trajectory were chosen in three ways: they could be practicaly. set, or defined from the "combustion chamber choking" condition or optimized. So the operating The operating mode of SCRJ-R may be mode of ATR on the first phase of trajectory (M < 5) optimized in the whole range of Mach numbers. It was was set by choosing air excess air coefficient equal to characterized by two parameters: total air excess 1.0. The operating mode of SCRJ and SCRJ-R were coefficient and oxygen-fuel excess coefficient in 3-9
Table 2, for SCRJ-R in Table 3. In the tables the a lPsgenerator. values of total effective thrust. of specific impulse and The results of computational study of of angle of attack along the flight trajectory are gven. optimum operating modes are shown: for SCRJ - in
Scramjet Performances Table 2.
Scramjet--Rocket Performances Table 3
Mach number range for SCRJ and SCRJ-R effective As a result of optimization of SCRJ hyper- operation show that optimum transition Mach number sonic operating mode, the air excess coefficient (to LRE) makes up: for SCRJ 14 and for SCRJ-R 15.2. remains constant and is equal to 1.0 at acceleration phase from M. = 7 to M. = 10. (Table 2.) Futher with The simultaneous operation is also useful: of the speed increasing SCRJ operates on a fuel enriched LRE and ATRgg in transsonic flight range, of LRE by hydrogen. and SCRJ -- by Mach numbers 13.6... 14.0 and of LRE and SCRJ-R by Mach numbers 14.8... 15.2. Optimization of SCRJ-R operation modes gives following results: in the concluding acceleration phase (M = 12 ... 15) for providing high thrust values engine may operate on oxygen enriched propellant. Comparative Evaluation of Fliaht In moderate hypersonic speed range (6 < M Performances of Aerospaceplane with SCRJ and < 11) corresponding to constant dynamic pressure q = SCRJ-R 75 kPa the optimum chosen on payload criterion is represented by highly economical operating mode Basic aerospaceplane flight performances without on-board oxygen supply. with two PS concepts -- SCRJ and SCRJ-R -- are shown in Table 4. The results of computational optimization of 3-10
SCR.J ani -7R. R Performances Table 4.
Flight performances parameters I PS c cepts SCRJ SCRJ-R PS -- ASP the main parameters Start accelaration, ndsq. s 6.25 6.25 Start relative loading, kg/sq.m 603 620 Landing relative loading, kg/sq.m 192 203 Mass ballance parameters, tonns Pay load mass 1.4 4.6 Mass in orbit 91.1 93.8 Propulsion mass 20.1 20.9 Air fraini mass 39.0 37.9 Tank mass 9.4 9.1 Propellents mass, tonns Oxygen for ATR-gg 15.7 15.7 Hydrogen for ATR-gg 15.7 15.7 Oxygen for SCRJ--R _--- 21.3 Hydrogen for SCRJKCRJ--R 45.6 43.8 Oxygen for LR 101.5 82.5 Hydrogen for LR 14.43 11.8 Total oxygen 117.3 119.5 Total hydrogen 76.1 71.3 Geometrical parameters Vehicle length, m 65.7 64.8 Wings span, m 29.1 28.7 Vehicle valume, cubm 1564 1500 Flight performances Accent time, s 1345 1340 Accent distance, km 435 1 4428 Speed losses due to drag, nds 3367 3210
One can see from Table 4 that SCRJ thrust Some Conclusions. augmentaton (use of SCRJ-R) allows to diminish aerospaceplane volunie what is particularly important 1. A 4-5 t payload may be launched to the under conditions of high tanks. and airframe elements near-earth orbit (H=20 km, i=51 deg). by specific mass. aerospaceplane with starting mass of 280 t using So using SCRJ-R instead of SCRJ allows to oxygen boosted SCRJ, i.e. SCRJ-R. This payload mass reduce total hydrogen mass by 5 t, aerospaceplane is three times greater than that launched by ASP with volume by 64 cub.m, aerospaceplane wetted surface by not oxygen SCRJ 43 sq.m, airframe and tanks mass by 1.5 t. All these improvements ensure launching into near-earth orbit 2. Optimuni Mach number for transition to 4.5 ..St of payload, i.e. some three times greater than LRE operation mode is for SCRJ - 14 and for SCRJ-R 0 the payload launched when using SCRJ only. - 15.2. Using oxygen as additional fuel component leads to larger the effective operating diapason of SCRJ-R in comparison with not oxygen SCRJ that.
3. Using of SCRJ-R permits to increase volumetric efficiency of aerospaceplane and to reduce it's structure mass. 3-11 Biblionruth E 15. Schetz T.A. and Billig F.S. "Studies 01 Scramjet Flowfields" AlAA Paper No. 87- 2 16 I, 1987. l.-Kyp3~~epP.M. "PeaKTMBHbIe ABHlXTeJIM AJlH 6on bUI Mx C BepX3BYKOBbIX C KO POCTe fi IlOJleTa" 16. Billig F.S. "Combustion Processes in M. Muwutiocmpoetiue, 1989. Supersonic Flow" Jet Propulsion, Vol. 4, No. 3, 1988.
2. Sosounov V. "Study of Propulsion for High 17. Pouliquen M.F., Doublier M., Scherer D. Velocity Flight" ISABE, 1991. "Combined Engines for Future Launchers" AIM Paper. 88-2823, 1988. 3. Escher W.J.D., Teeter R.R., Rice E.E. "hrbreathing and Rocket Propulsion Synergsni: 18. Roudakov A., Dulepov N. et.al. "Analysis of Enabling Measures for Tomorrow's Orbital Transports" Efficiency of System with Oxygen Liquifaction and AlAA Paper, 86-1680, 1986 Accumulation for Improvement of Aerospaceplane Performances" 1AI;-91-270, 1991. 4. Bendot J. "Composite Engines for Application to a Single-Stage-to-Orbit Vehicle" NASA CR-2613, 19. Czisz P.A. "Space Transportation Systems 1975. Requirments Derived from the Propulsion Perforance" MI;-92-0858. 5. Murthy S.N.B., Czysz P.A. "Energy Analysis of High speed Vehicle Systems" AIAA Paper, 90-0089,1990. 20. Vinogradov V., Stepanov V. "Numerical and Experimental Investigation of Airfranie-Integrated Inlet 6. Hewitt F.A., Ward B.D. "Advanced for High Velocities". AIAA Paper 89-2679. 1989. Airobreathing Powerplant for Hypersonic Vehicles" ISABE 89-7107, 1989. 2 1. Pouliquen M. "Propulsion of Space Transportation System", Paris, 1991. 7. Heppenheimer. "The National Aerospace Plane" Pasha Publikations, 1989. 22. Wolf D.M., Dauni A., et.al. "Air-Borne Launching with Existing technology", lAF-93-V.3.615, 8. Czysz P.A. "Thermodynamic Spectrum of 1993 Airobreathing Propulsion" SAL' T.1'. 881203, 1988. 23. Billant C., Boueldleu E., Wagner A., 9. Saber A.J. and Chen X. "H2/Air Subsystems "Radiance: A High Staging Speed Airbreathing First Combustion Kinetics in Aerospaceplan Powerplants" Stage TSTO Launcher", IAF-93-V.3.616, 1993 lA(AF-91- 276, 1991. 24. Roudakov A., Novelli Ph., et.al., "Flight IO. Robert A. Jones, Caleman du P., Donaldson. testing an Axisymnietric scramjet - Russian Recent "From Earth to Orbitin at Single Stage" Aerospace Advances." IAF- 93-S.4.485, 1993 a America, August 1987. Vol. 25, No. 8. 11. Rudakov A., Krutchenko V. "Additional Fuel Component Application for Hydrogen Scramjet Boosting" SAE T.P. 900990, 1990.
12. Hewitt F.A.. Ward B.D. (Rolls-Royce). "Advanced Airbreathing Powerplant for Hypersonic Vehicles" ISABE 89 -- 7107, 1989.
13. Kors D. (Aerojet Propulsion Division). "Design Cosiderations for Combined Air Breathing-Rocket Propulsion systems" AIM Paper No. 90-5216,1990.
14. Bapra@TMK H.E. "CnpaBOYHMK no TCnJlO@M3M'WCKMM CBOGCTBaM M30B M XwKOCTefi'' M. @u3~umzu3,1972. 3-12
t
- t .5 AV
0.6 LO 0.4
U. 5 0.2
0 10 I5 20 2 40 r5 20 5
Fig. I. 1. Relative energy of hydrogen combustion and theoretical relative limit of flow velocity increase for scramjet.
5 10 15 20
Fig. 1.2. Influence of various losses on scramjet specific impulse along the flight path. 3-13
Pi Pi
Mf
Fig. 1.3. lnlct performances. 0. .2 .4 .6 .8 1 .0 l/BIlO
Fig. I 5. Hydrogen oxygen gas exhaust specific impulse.
20 -
3 $ 52 3 10 - E a,U 2 3 +Ha0 +Ne + Nz 0 1 I I I +Ar 5 10 15 20 Mf +de Fig. 1.6. Specific impulses of hydrogen scramjet for -- I - normal and tangential fuel injection with .2 .4 .6 .8 1.0 different fuel flow ratios. d
Fig. I 4. Gas generator cshaust specific iinpulse as function of rclative additional coniponcnt coiisuinption for diflcrent propcllant additions. 3-14
- 20 20 2 9 10 - 10 n 3 E E c-( U
0 t I 0 4 10 15 l Mf
Fig. 1.7. Specific impulses of hydrogen scramjet and Fig. 1.9. Specific impulses of hydrogen scramjet and a hydrogen-neon scramjet with different fuel hydrogen-helium scramjet with different fuel flow ratios. flow ratios.
Fig. 1.8. Density impulses of hydrogen scramjet and Fig. I. 10. Density impulses of hydrogen scramjet and hydrogen-neon scramjet with different he1 hydrogen-neon scramjet with different fuel flow ratios flow ratio. 3-15
I "2 l-. --- / n -51
U 4 -2.01 1. 2. 3 Fig. 1.1 1. Spccific impulses diITerences of bipropc lant scranijct and pure hydrogen scramjet for Fig. 1.13. Specific heat variations (divided by that a1 0 taiigeiitial fuel injcction at B = 2.5. 400 K). Mf 20
n
Uz
Fig. 1.14. Specific iiiipulses differences of bipropellaM scramjet and pure-hydrogen scramjet for noriiial fuel injection at B = 2.5. n 3: 2- .4. > U II -d" .2-
0 10 15 20 Mf
Fig. I, 12. Density inipulses differences of bipropcllillll scramjet and pure-hydrogen scranijet for nornial fuel injection at B = 2.5.
1s 0 Mf
Fig. I. 15. Density inipulses differences of bipropeilanl scranijet and pure-hydrogen scmijct for normal fuel injection at B = 2.5. 3-16
Fig. 1.16. Specific impulses differences of bipropellant Fig. 1.18. Specific impulscs differcnccs of bipropcllant scramjet and pure-hydrogen scramjet for scramjet and pure hydrogen scramjet for tangential fuel injection at B = 2.5. normal fuel injection at B = 2.5.
.8 .8 E E D D tl x 2 .(i $ n z> 3 .4 .4 U -I I > U> U II 'I> .2 .2 U Q2- Q
1 0 1 I 5 10 15 Mf Mf Fig. 1.19. Density impulses differences of bipropellant scramjet and pure-hydrogen scramjet for Fig. 1.17. Dcnsity impulscs difTcrcnccs of bipropcllanr normal fucl injection at B = 2.5. scramjet and pure hydrogen scramjet for tangential fuel injection at B = 2.5. 3-17
5 10 15 20 Mf Fig. 1.20. Specific impulses of hydrogen scramjet and hydrogcn-osygcn rockct scrainjct for hydrogen tcmpcrature 1000 K.
0 L __ -1.. - 1.---- 5 IO IS 20 Mf Fig. 1.2 1 Dcnsity impulses oE hydrogen scramjet and hydrogen-oxygen rocket scramjet for hydrogen tcinperature 1000 K.
Mf
Fig. 1.22. Specific inipulscs of hydrogcn scranijct and hydrogen-oxygen rocket scramjet for hydrogen tcnipcraturc 2000 K. 3-18
--...
--( -- 3-19
The A.S.P. .2G DRAG FACTOR .18
.16 .14 i .12
.io
.06
.06
.04
.Q2
. oc I.----- I 2 4 E 8 {Q 12 14 i6 '18 z 3
FLIGHT flACH NUflBER Fig. 1.24
The A.S.P. LIFT FACTOR .7 a
.J
.2
..1
0 I 1 I I I 2 q 6 8 10 12 14 16 18
0 FLIGHT flACH NUMBER Fig. 1. 25 3-20
The A.S.P. 1C AERODYNAMIC FINENESS FACTOR 5 MAX IMUM and a CORRESPONDING 7 ANGLE OF ATTACK CINGLE OF ATTACK 6
3
4 FINENESS
3
2
.!
0 I 1 I 1 I 1 0 2 4 6 8 13 12 14 !6 18 20
FLIGHT flACH NUMBER Fig. 1. 26 4-1
SCRAMJET CFD METHODS AND ANALYSIS 0 PART 1. SCRAMJET CFD METHODS
NUMERICAL SIMULATION OF THE FXOW IN SCRAMJET DUCT
bY
V.KOPCHENOV, K.LOMKOV, L.MILLER, V.KRJUKOV, LRULEV, V.VINOGRADOV, V. STEPANOV, N.ZACHAROV, R.TAGUIROV, M.AUKIN
CIAM (Central Institute of Aviation Motors) Aviamotornaya street 2 111250 MOSCOW RUSSIA
1.1 Introduction 1.2 Forebody and inlet flow model a 1.3 Combustor CFD model structure 1.4 Nozzle and afterbody CFD model
Introduction incorporate the model of hydrogen-air combustion. Air- hydrogen chemical processes simulation is important in The computer analysis of scramjet flow became considering the flow in combustor and nozzle. Thus, it of great importance because of the limited possibilities is desirable in general case to solve the three- of ground tests and difficulties of measurements in high dimensional Reynolds-averaged Navier-Stokes speed/enthaIpy flows. This fact is a powerful stimulus in equations for multicomponent flow with adequate CFD development on the other hand. The short system of chemical reactions and turbulence model. description and examples of applications of the mathematical model for scramjet duct flow developed in This level of numerical simulation requires the CIAM are presented in this paper. usage of the most powerful modem supercomputers. We have used the computers with like personal one power Before formulating the main requirements to the in our work. Therefore, the system of sufficiently simple mathematical model it is necessary to take into account models has been developed to evaluate the influence of the following circumstances. The airbreathing the main processes, regime and design parameters on .hypersonic vehicle is characterized by high level of performances of propulsion system elements. This way a integration of airframe and propulsion system. leaves some questions open, but seems to be justified as Therefore, the correct prediction of propulsion system really possible for preliminary design evaluations. performances demands to consider the forebody pre- compression of air entering the scramjet engine. The Recently some surveys about CFD in afterbody of the vehicle is a nozzle extension. Thus, the SCRAMJET applications were published [1-6]. We mathematical model of the scramjet duct flow must describe our own experience in this line in this report. include the description of the forebody and afterbody The system of mathematical models and codes flows. developed in CIAM has the modular structure and includes such elements: forebody, inlet, combustor, At the numerical simulation of scramjet duct nozzle and afterbody. Comparatively simple, flow, it is necessary to take into account the three- engineering models are used on the first step of the dimensional effects, the real gas effects, viscosity propulsion system performance evaluations. This is the influence, turbulent mixing processes, chemical level "0" in accordance with [5]. Another models allow reactions. The complex gasdynamical structure of the to evaluate the influence of three-dimensional effects, flow and viscosity influence require to use the most real gas effects, to estimate the role of physical-chemical general hydrodynamic model - the Navier-Stokes processes. The most part of these codes is based on the equations with Reynolds averaging. The turbulence Euler and boundary layer equations and some model adequate to high supersonic Mach numbers flows procedures of viscous-inviscid interaction. This is level is required. It is necessary to take into account the "1" in accordance with [5]. Several codes are based on complex physicalchemical processes in air flow over the Reynolds-averaged parabolized Navier-Stokes the forebody and in inlet. It becomes necessary to equations and may be attributed to level "2" in
Presented at an AGARD Lecture Series on 'Researchand Development of RadScramjets and Turboramjets in Russia: December 1993 to January 1994. 4-2
classification mentioned above. Finally some codes are equilibrium in addition to usual Euler equations [8,9]. developed for the numerical solution of the full Navier- More simple method is based on the use of analytical 0 Stokes equations in 2-D case. These codes may be used approximations for the air properties [ 101. This method in local regions where strong viscous-inviscid is implemented in combination with Euler equations in interaction is sigruficant. Besides, they allow to evaluate some CIAM codes for forebody and inlet flows. the influence of effects neglected in the models of lower level. It can be done on the examples of special test Alternative method was proposed in [ll] for - problems. approximations of gas properties in equilibrium approach. In this case effective specific heat ratio It is necessary to note, that dividing on the levels (K=a2p/p where a is a frozen speed of sound, p is - is conventional in our case. Indeed, the more simple pressure and p is density), specific enthalpy h and physico-chemical model may be used with more molecular mass p are represented as functions of two complex hydrodynamic model wavier-Stokes equations successfully chosen independent variables. One of them versus Euler equations for example) and vice-versa. is the decimal logarithm of pressure, and another - the- Moreover, the combustion model of higher level may be special function of pressure and density. It is possible to used in 2-D case in comparison with 3-D case in codes obtain the values of K, h, CL as fhctions of chosen for the numerical solution of Reynolds averaged independent variables with the aid of approximation parabolized Navier-Stokes equations. procedure, based on some "node" points. These "node" points are stored as tables, in which the main The development of alternative versions of parameters K, h, p are represented in sufficient number mathematical models and the comparison of the results of points in the space of independent variables to obtained with the aid of various models are useful, provide the acceptable approximation accuracy. The because the most part of models is approximate. The proposed method provides the accuracy 1-2% in wide experimental testing of existing mathematical models is range of pressures ( from 100 Pa up to lo8 Pa ) and very important problem. This is necessary to establish temperatures (from 200" K up to 20000" K ). the credibility of CFD results used in design [7]. Two codes are used for the numerical simulation of the flow over forebody (developed by V.Kopchenov, 1.2 Forebodv and inlet flow model V.Kr~ukov,K.Lomkov and L.Miller). The first code for the computation of sub-, tran- and supersonic flows is The main requirements to the airspace plane based on time relaxation method. The region near the forebody design may be formulated. The forebody must blunt nose is calculated with the aid of this code. The have the minimal drag at fixed volume. The lower numerical method is the modified version of the well- surface of forebody must precompress the flow with known Godunovk scheme [12]. The piecewise linear minimal losses, to provide the required massflow and distributions of the main parameters in all spatial uniform, as far as possible, air flow at the inlet entrance. directions instead of piecewise constant distributions in The forebody must have blunt nose part. Bluntness of each computational cell are supposed [13]. The the nose must be chosen on the base of compromise modified minimal increments principle [14] is used to between demands to be stable to the action of high heat provide the monotonicity condition [15]. The explicit fluxes and to decrease the wave drag and to provide the higher order accuracy scheme proposed in [16] for uniform air flow at the entry of inlet. From practical hyperbolic systems of equations is realized as predictor- point of view, it seems to be justified to choose the corrector scheme in accordance with [17]. minimal acceptable bluntness from the heat fluxes analysis and then to shape a forebody in order to satisfy The modified version of shock fitting procedure, other requirements. which is based on some principles proposed for 2-D case in [18] is developed for 3-D case. The proposed Thus, the numerical model for the forebody must method provides the second order accuracy on a regular provide the opportunity to analyze the influence of blunt uniform grids and conserves approximation on arbitrary nose, the three-dimensional effects, the influence of nonuniform irregular grids. The local time step is physicochemical processes in high temperature air, and chosen for each cell to decrease the run time necessary effects of viscous-inviscid interactions. to obtain the steady state solution with the aid of the time relaxation procedure. The experience shows that The existing system of codes is based on local time step marching provides almost the twice gain numerical solution of the 3-D Euler equations. These in run time, when the transition from one flight Mach codes are modified to include real gas effects in local number regime to another is accomplished. Moreover, equilibrium approach. It is supposed that gas this method with local time step is comparable by the composition, molecular mass, heat capacity and speclfic number of time iterations with the implicit method [19], enthalpy are functions of local pressure and based on the scheme with space directions and physical temperature. To complete the problem it is necessary to processes splitting and with global time integration step. use the system of equations for the thermodynamic 4-3
The marching method is used in the regions with Fig. 1.2.4 in four cross-sections to illustrate the supersonic longitudinal velocity component. This influence of the shape change on the flow parameters method is the higher order accuracy version of the distributions. One can see the generation of the internal steady analogy of Godunov's scheme [12]. The same shock and region of high compression within the shock model [ll] of equilibrium air is incorporated. The layer connected with the deformation of the upwind side modified version of the scheme for 3-D steady of the forebody. Therefore the visible nonuniformity in - supersonic flows is based on the same principles as the parameters distribution may be obtained at the entry of scheme for the unsteady case [13, 14, 16-18]. The base inlet. version of this scheme for 3-D steady flows of perfect - gas was proposed by A.Kraiko and S.Schipin. The entropy layer which is formed on blunt nose part is also unfavorable factor for the inlet operation. It is necessary to evaluate the influence of the The conventional "total pressure" is introduced to viscous effects on the flow fields. The simplified method distinguish the entropy layer. This value is calculated - was implemented until recently. The boundary layer for locally frozen flow parameters. Such "total pressure" displacement thickness is evaluated in approximate distributions in some cross-sections are shown in Fig. manner along generatrices of the body in each cell 1.2.5. It is necessary to note, that the entropy layer nearest the wall. Then the body surface is corrected on occupies the all shock layer in cross-sections near the the boundary layer thickness and new "external" nose. But at the end of the conical part of the forebody, inviscid flow is calculated and so on. Such iterative the major part of the entropy layer flows on the leeward procedure must be repeated up to convergence in some side of the vehicle. sense. This iterative procedure was developed and realized in [20]. The integral turbulent boundary layer The calculation with the boundary layer method [21] allows to estimate the boundary layer estimation shows the following. The boundary layer displacement thickness, the momentum thickness, the displacement thickness in the last cross-section (see Fig. skin friction coefficient and the specific heat flux to the 1.2.44) constitutes approximately 5% from the shock wall. layer thickness on the upwind side of the forebody. These estimations were obtained by NZacharov. The flowfield is calculated up to the inlet entrance. It is possible to evaluate in this cross-section It is necessary to note that all these calculations the entropy layer thickness, the extent of parameters were performed on the IBM PUAT-486. The blunt nose nonuniformity, the possible inlet area, the boundary flow was simulated on the grid containing 5000 nodes layer thickness, and forces acting on the forebody and and it was required for this calculation about 2 hours of additional forces acting on the streamtube captured by run time. The run time with the marching method up to inlet. The last is necessary for the analysis of control the inlet entry constitutes approximately 5 minutes on volume forces. the same computer.
Typical forebody solutions are illustrated in The numerical simulation of the inlet flow is following figures. The flow over the sphere is calculated diliicult task. This may be explained by several reasons. for two hypersonic regimes with flight Mach numbers One of them is a wide range of operational regimes at MflO for the altitude H=10 km (Fig.1.2.1) and with flight Mach numbers variance from 0 up to 20-25. It is M~20for H=20 km (Fig.1.2.2). Results are compared necessary also to note the essential part of viscous- with the primary standard results obtained in [9]. The inviscid interaction in the flow with deceleration. This pressure and density distributions on the surface of part is most important on regimes of starting and sphere are compared. The pressure and density are unstarting. Besides that, the large nonuniformity in related accordingly to dynamic pressure and density in parameters distributions may exist at the entrance the free flow. The results are obtained on the grid because it is impossible to control actively the boundary containing the lOxlOxl0 nodes by means of the base and entropy layers. Godunov's scheme and its modified version. The distance from the sphere surface to the bow shock along But on the other hand, it is necessary to analyze the sphere radius direction is also shown as a function large number of inlet alternative schemes at the stage of of the polar angle in Fig. 1.2.1-c and 1.2.2-c. It is preliminary design. Therefore it is desirable to use some necessary to note the higher level of accuracy of simplified mathematical model of inlet flows. modified scheme in comparison with original base Nevertheless it is necessary to take into account the version scheme. main factors influencing the inlet performances at simplification of mathematical model. These main The second example corresponds to the inviscid factors are: the flow spreading over a forebody and forebody flow. The geometrical shape of the blunted additional drag concerned with this effect; the three- forebody is presented in Fig. 1.2.3. The calculations dimensional effects in the inlet flow; the real gas effects were performed for MplO at an angle of attack 3.9" for in high temperature flow; the bounw and entropy the altitude 30 km. The pressure fields are presented in layer nonuniformity at the inlet entrance; the boundary 4-4 layer influence. reduced. Relative inlet throat area At (ratio of the area of captured stream tube at the inlet throat to the inlet 0 In accordance with the forebody flow model, it is frontal area) was about 0.21 taking into account the possible to take into account the major part of the mentioned expansion. factors mentioned above by using the following inlets flow mathematical model. This mathematical model Numerical investigation of inlet flow was carried includes: the code for the 3-D calculations of steady out for Mp4-8 with the grid in the plane ZOY - supersonic flow in inlet on the base of marching containing 24*80 nodes. Calculations of integral Godunov's type scheme [22]; the real gas effects in high characteristics were made according to the conservation temperature flow are analyzed with the aid of equations written for control volume, which was defined equilibrium air model [lo]; the boundary layer by base plate, cowl, side wedges, entry and exit sections - estimations on wetted surfaces are performed assuming and upper surface limiting the freestream tube. While that the flow is attached with the aid of integral method calculating C,, (additive drag coefficient), it was [21] in accordance with procedure proposed in [20]. assumed for convenience in defining the airframe - Some examples of proposed mathematical model integrated inlet drag, that the additive drag force applications for inlet flow fields and performances includes not only the force acting on the liquid surface predictions are presented in [20,22]. but also the force acting on parts of sidewalls (marked region in Fig.1.2.6-b) wetted by the flow that does not It is necessary to note that mathematical model is enter the inlet. Averaging was made with conservation approximate and following justification is required for of the flow rate, total enthalpy and entropy in real and the recommendations have been made on the base of averaged flows. this model. But it is possible to "optimizett inlet configuration on the stage of preliminary design with The total inlet characteristics are shown in the the aid of this simple and computationally fast model Fig.1.2.7 for three-strut inlet with sweep angle 48'. One and to diminish the number of variants for the following can see a varying behavior of q=f(Mr) connected with investigations. the increased recovery losses at Mp5 because of critical flow regime in the central passage near the cowl. But it The applicability of this configuration must be justified by experimental investigation or with the .aid of is necessary to note, that the evaluation of the possible influence of the viscous-inviscid interaction on the main some more full (if possible) and more expensive mathematical models. Numerous examples of proposed inlet performances is required. Then it is necessary to method application to solve various gasdynamic use the reliable experimental data or more full, problems show that local and integral flow complex and expensive mathematical models. This fact characteristics are defined with an accuracy sufficient is confirmed by the comparison of experimental and for practical applications. Particularly, the computational pressure distributions along the base computational errors in the definition of integral plate and side wedge of considered inlet (see Fig. 1.2.8). characteristics (inlet capture ratio that is equal to ratio Unfortunately, there is no acceptable mathematical model now to evaluate the inlet performances for of areas of captured stream tube at the inlet face and inlet frontal area, mass averaged total pressure recovery regimes near starting and unstarting. The experimental investigations are the single possible method in this factor, additive drag coefficient ) do not exceed 0.5 - 3 %. case.
The calculation of the single inlet performances was made [22]. The inlet configuration is shown in Fig. 1.3 Combustor CFD model structure 1.2.6. The computational region is divided into some subregions. The computational region at zero slip angle The numerical simulation of the combustor flow is presented in Fig. 1.2.6-a. The inlet sidewall that requires to take into account: the turbulent mixing; the physicochemical processes in hydrogen-air mixtures; begins in section x=O has a sharp leading edge with the complex gasdynamic structure of the flow; the three- sweep angle x. The cutback cowl (x=L,) has also a dimensional effects. sharp leading edge with sweep angle ~~45'. The sidewall and cowl are wedges with angles 8=5.3O and It is necessary to note, that two characteristic p=lo in planes XOZ and XOY (Fig. 1.2.6-a). Cowl combustion regimes can be considered in the scramjet internal surface coordinates (upper duct wall) and base combustor. The combustion process in duct at small plate (lower wall) are specified and, in our case, supersonic combustor entrance Mach numbers correspond to planes y=h and y=O respectively. corresponding to flight Much numbers 6-7 realizes basically at subsonic velocity. This is caused by the The struts dividing the inlet into central and side pseudoshock wave system existing due to the mass-heat passages (for the half of the inlet) have the same sweep supply. The effects of interaction of shock wave system angle as side walls. To reduce the excessive flow with boundary layer are of the great importance in the compression near the cowl, the central strut height is flow regime realization. Combustion regime with the 4-5 creation of extensive subsonic regiones and recirculation this paper. Some results obtained with the aid of the zones takes place as a result. The estimations of the developed system of codes are presented in Fig.l.3.1- combustion efficiency are made in CIAM for these 1.3.3. The combustion efficiency versus longtudinal regimes on the base of experimental data generalization. distance is presented in Fig. 1.3.1 for axisymmetric flow in the cylindrical duct. The regime is characterized by At high flight Mach numbers combustion occurs following parameters in the initial cross-section: in supersonic stream, when the velocity remains hydrogen jet - Il&=2.9, T,=373" K; supersonic everywhere including combustion zone with air stream - Mi=4.8, Te=1500" K; the exception of wall boundary layers. The estimations pressures ratio in the hydrogen jet to air stream is equal of the combustion efficiency for these regimes are to 1.92. This regime corresponds to conditions at the carried out using mathematical model, which is based entry of supersonic combustor chamber for conventional on Reynolds averaged parabolized Navier-Stokes vehicle at flight Mach number 12. The calculations equations. were performed for two levels of turbulent viscosity (~~=3*10~and 2.5*103) in initial cross-section. The It seems to be justified to develop the system of turbulent viscosity is referred to the hydrogen velocity mathematical models and a number of codes to evaluate and the hydrogen jet radius in initial cross-section. It is the performances of scramjet combustor and the interesting to point out the higher level of combustion influence of the regime and design parameters on efficiency for higher level of turbulent viscosity. This combustor performances. This system must include both example demonstrates the influence of the mixing the simple engineering models, based on generalization process on the combustion efficiency. Therefore the of experimental data, and the most full models, which mixing enhancement becomes of primary importance are based on modem numerical simulation of supersonic for combustion efficiency augmentation. combustion. The influence of the heat release on the Apparently, the most general model must be gasdynamical structure of the flow in the duct is based on the 3-D Reynolds averaged Navier-Stokes illustrated in Fig.1.3.2 and 1.3.3. The calculation was equations with acceptable for high Mach numbers flows performed in the first case (Fig.1.3.2-a and 1.3.3-a) turbulence model, for multicomponent mixtures with only with turbulent mixing without combustion. In the sufficiently detailed scheme of chemical reactions. But second case, the turbulent combustion was simulated this level is inappropriate to the available computers. At with the aid of the flame sheet model (Fig.1.3.2-b and the same time, some simplifications may be used, if only 1.3.3-b). The pressure (Fig.1.3.2) and Mach number regimes with combustion in supersonic regions are fields (Fig.1.3.3) are shown. It is interesting to note, considered . Then parabolized Navier-Stokes (PNS) that the system of shocks in duct is much more distinct equations may be used for the numerical simulation. It in the case of heat release. Moreover, the longitudinal is possible to use effective marching numerical methods dimensions of "periodical" structures becomes smaller. and to realize these models on available computers. The developed mathematical model may be used It seems that the most realistic way is to develop for the comparative study of various methods of mixing the hierarchy of models including: and combustion enhancement. This model may be - engineering models based on one-dimensional implemented to evaluate the influence of the flight calculations of combustor duct flows with known curve regime parameters on the combustor efficiency. It is of combustion efficiency obtained by experimental data possible to perform the efficiency comparison for some generalization (see, for example, [23]); versions of injection system in combustor. Some CIAM - codes for 2-D calculations with turbulent mixing and experience in supersonic combustion model applications combustion simulation; is presented in the second part of this report. - codes for 3-D effects evaluations.
The second and the third levels may be based on 1.4 Nozzle and afterbodv CFD model the PNS equations. But the more detailed combustion models may be incorporated into 2-D codes. In codes for The exhaust system of airspace plane is 3-D flows, it is justified to use the more simple characterized by the following features: the integral combustion models on the first stage. For example, the scheme of the nozzle and airframe; the principle of diffusion flame sheet model seems rather attractive. But combined propulsion system; the 3-D schemes. the opportunity to estimate the applicability of this simplified model and to evaluate the influence of the In the integral scheme of exhaust system and finite rates of chemical reactions on the combustion airframe, the surfaces of the nozzle are at the same time process must be provided. the airframe surfaces. Therefore the exhaust system performances must be determined taking into account The more detailed description of the the afterbody flow. The principle of combined mathematical model is presented in the second part of propulsion system supposes that various engines or their 4-6 combinations may be used on some sections of flit t constitutes in this case more than YOfrom the ideal trajectory. In this case the exhaust system is a nozzles thrust. 0 combinations of various engines (liquid rocket, turbojet, scramjet), operating simultaneously or sequentially in The flow picture is shown in Fig.1.4.2 for time. The mathematical model must provide the another conventional regime. It is supposed that the determination of longitudinal and transverse forces, combined propulsion system operates at flight Mach moments of forces, skin friction forces and heat fluxes number Mf2 in such manner: the liquid rocket and . on the walls of the exhaust system, thrust direction and turbojet engine are operating, and massflow bypass is its location point. realized through the scramjet duct. The external flows over afterbody and nozzle cowl are also considered in - The proposed mathematical model may be used this case. Besides that, the base region exists near the - to evaluate the influence of two factors mentioned base face of afterbodyexhaust system configuration. above. This simplified mathematical model for 2-D The thrust losses constitutes in this case approximately flows is based on the following principles. The 20 %. computational region is divided into several subregions. The supersonic inviscid flow is described by the steady The Mach number contours are shown in Euler equations, which are solved numerically with the Fig. 1.4.3 for conventional regime, corresponding to aid of Godunov's type marching scheme. The version of MflO. Only scramjet engine operates in this case. The the scheme is used, which was proposed in [24]. The performances of the exhaust system were evaluated for integral method is used to describe the boundary layer the conventional regimes of Mr from 4 up to 16. The and mixing layer flows. The skin friction and wall heat exhaust system configuration presented in Fig.1.4.3 is fluxes are evaluated with the aid of the integral method considered. The thrust losses are obtained in comparison with ideal thrust as function of flight Mach number and are shown in Fig.1.4.4. The dependence of The method [25,26] was developed to determine the slope of the thrust vector on flight Mach number is the base pressure and base enthalpy behind the presented in Fig. 1.4.5. backward step, which is flowed by two different supersonic streams. The method allows to determine the It is necessary to note, that the developed base pressure and enthalpy taking into account the mathematical model is approximate and has serious influence of Mach numbers, Reynolds numbers, initial restrictions. This model is being modified now to boundary layers thicknesses in main streams. The total include the three-dimensional and real gas effects. The pressures and temperatures are different in two flows. experimental investigations are required especially for The gas is considered as perfect with Merent specific regimes with large separation zones to ver@ the heat ratios in two streams. The method permits to numerical codes. estimate the influences of mass blowing into the base region and heat flux through the wall of the base face on the base pressure and enthalpy. The model of viscous- inviscid interaction [25,26] was used to determine the parameters in the base region. The original condition in reattachment point was proposed [25].
In general case, the configuration of the exhaust system-afterbody involves the base face. Moreover, the exhaust system is considered as uncontrolled. Therefore, large base regions may arise in exhaust system flow on the acceleration part of flight trajectory, when scramjet engine does not operate, and some additional engines of combined propulsion system operate. It is known that large thrust losses may exist in this case.
The possibilities of the proposed mathematical model are illustrated in Fig.1.4.1-1.4.5. These results are obtained by RTaguirov and M.Aukin. Some estimations were performed for conventional vehicle exhaust system in the case of turbojet engine operation. The flow field is illustrated in Fig.l.4.1 for supposed regime, corresponding to flight Mach number 1.3. The bypass of the gas was realized through the duct of the non operating scramjet engine. Two base regions in the exhaust system flow take place. The thrust losses 4-1 a References 1. White M.E., Drummond J.P., Kumar A., 60nmecm CHCTeM C AByMII He3aBHCHMbIMI.I "Evolution and Application of CFD Techniques for nepeMeHHbIMH", XyPHaJ'I BbIYHCJIHTeJIbHOii Scramjet Engine Analysis", J. of Propulsion and Power, MaTeMaTHm H MaTeMaTmeCKOii @i3HKIi, T.23, V.3, N 5, 1987, ~p.423-439 N 4, 1983, c.848-859. 2. Schetz J.A., Billig F.S., Favin S., "Modular 17. POAHOHOBA.B., "MOHOMHH~cxeMa Analysis of Scramjet Flowfields", J. of Propulsion and BToporO nopsiqKa annpoKcHMaqHH rn Power, V.5, N 2, 1989, pp.172-180. CKB03HOrO PaCYeTa HepaBHOBeCHbIX TeYeHEiii", XypHm BbIYHCJIHTeJIbHOfi MaTeMaTHm . 3. Barber T.J., Cox G.B.,Jr., "Hypersonic H Vehicle Propulsion: A Computational Fluid Dynamics MaTeMaTmecKoii @H3Hm, T.27, N 4, 1987, Application Case Study", J. of Propulsion and Power, c.585-593. V.5, N 4, 1989, pp.492-501. 18. KpaiiKo A.H., MaKapoB B.E., 4. Povinelli L.A. "Advanced Computational Tmesa H.H., "K YHCJIeHHOMY nOCTpOeHMI0 Technology for Hypersonic Propulsion", ISABE 89- ~POHTOB YAapHbm BOJIH", X~PHU 7106, 1989, pp.993-1008. BbIYHCJIl€X?JIbHO~MaTeMaTHm Ei MaTeMa- 5. McClinton C., "CFD Support of NASP TEiYeCKOfi @H3EiI(zI, T.20, N 3, 1980, c.716-723. Design", AIAA 90-5249, 1990. 19. KOBeHH A.A., HHeHKO H.H., "MeToA 6. Anderson G., Kumar A., Erdos J., "Progress pacueme~mB swayax ra30~0ii amaMHm", in Hypersonic Combustion Technology with HOBOCH~HPCK,HayKa, 1981, 304 c. Computation and Experiment", AIAA 90-5254, 1990. 20. BHHOI~~~OBB.A., A~HOBB.B., 7. Mehta U., "The Aerospace Plane Design CTenaHOB B.A., "npHMeHeHHe YHCJIeHHbIX Challenge: Credible Computational Fluid Dynamic MeTOAOB K pacYeTy XapampHcTHK Results", AIM 90-5248, 1990. CBepX3BJXOBbIX H lXnep3B)'KOBbIX B03fiYXO- 8. Ea6e~~oK.H., BocKpeceHcmii r.n., 3a6OpHHKOB BPA, YYeHbIe 3an~cmWM, JIIO~HMOBA.H., PyCaHOB B.B., "npOCTpaH- T.13, N 2, 1982, c.62-68. 21. ABnyeBcmii B.C., "MeTon pacYeTa CTBeHHOe 06T€?rPRESSURE DISTRIBUTION ON SPHERE SURFACE P SX y\
05 Mf=10, H=lO km
- - - I-STORDER SCHEME - 2-ND ORDER SCHEME a 0 0 0 191
I I ,w Fi ,g.1.2.1- A. 20° 40" 60" 8 0'
DENSITY DISTRIBUTION ON SPHERE SURFACE
.5
- I- ST ORDER SCHEME 2.5 - - - 2 - ND ORDER SCHEME 0 0 0 191
0 Fig.l.2.1 - B. a 0 20° 40" 60" 8 0' 4-9
AS BOW SHOCK DEPARTURE
0.4
0.3
,432 PRESSURE DISTRIBUTION ON SPHERE SURFACE /
1 - ST ORDER SCHEME 2 - ND ORDER SCHEME 191
I I I I Fig.1.2.2 - A.
DENSITY DISTRIBUTION ON SPHERE SURFACE
-4
--- I- 2- 0 0 0 191
0 Fig. 1.2.2 - B.
I I , ,U 0 40' 80 4-11 e BOW SHOCK DEPARTW
0 0 0 191
I I ,o 0 zoo doa 60' BO
Fig.1.2.2 - C.
THE FOREBODY SCHEME
.OD D
Fig. 1.2.3. 4-12
THE PRESSURE
I ”= 10.00 ALFRn 3.90 BETA= .OO cnn 101’10 P/PH-CONST
-vi-.-
Fig.1.2.4- A. THE CROSS-SECTION X=X3 Fig.1.2.4 - B. THE CROSS-SECTION X=X4
THE PRESSURE FIELDS
M= 10.00 w 10.00 ALFA- 3.90 ALFR- 3.90 BETA= .OO BETA= .OO cnn 10rq0 cnn 10.’10 P/PH=CONST P/PH=CONST 0
--e-- +--
Fig.1.2.4 - C. THE CROSS-SECTION X=X5 Fig.1.2.4 - D. THE CROSS-SECTION X=xg 0 4-13 a THE "TOTAL"PRESSURE FlELDS
m- 10.00 WR- 3.90 BRA- .OO
m 10.'10 Pm -CONST
THE CROSS-SECTION X=X1 Fig.i.2.5 -A. Fig. 1.2.5 - B. TI3E CROSS-SECTION X=X2
THE "TOTAL" PRESSURE FIELDS
IIH- 10.00 W- 10.00 ALFA= 3.90 RLFA- 3.90 BETA- .OO BETA- .OO
CM lODClO CM 10.W Pm -CONST
~e Fig.1.2.5- C. THE CROSS-SECTION X=X3 Fig. 1.2.5 - D. THE CROSS-SECTION X=xg 4-14
THE PERFORMANCES OF THE THREE - STRUT INLET
1 - CENTER PASSAGE 16 5.5 2 SIDE I- - - P 3 - TOTAL "
Oi W m W 5 4.5 E z v) v) I W 0 3 3.5 c- y I/----- # < cI cI I I I 2.5 I I I 45 6 7 8 5 6 7 8
Il.0 1.2 r 70.01 5 >- az 1 E 0 3 n W 0 2i $0.8 0.010 ". @ 0.8 I 8 W W 0 K 3 3 W In 5 0.4 odd 0.005 E M 0.6 w 0 E 1- 0 a Q I- I- W 2 \' -I I- 0.0 Ll I I I I I I I 5 7 0.000 0.4 5 6 7 8 6 8i MACH NUMBER, M, MACH NUMBER, M, Figl.2.7. e 4-15
STATIC PRESSURE DISTRIBUTION ALONG THE BASE PLATE
CENTER PASSAGE , ---. n /‘\
E I M,=6M, = 5.1 DATA ‘8WU 3 Vl I I- 34rY // U1I - fia ,D-]A\ &I
SIDE PASSAGE 7 SIDE PASSAGE
Fig. 1.2.8 - A.
STATIC PRESSURE DISTRIBUTION ALONG THE SIDE WALL
I I 60 - @ 0;
la 9 40- I ‘;z 0’ (L 01 W b IY I 20- y/h = 0.91 W [L a rrl 0 - _------I I I .o 1.5 2.0 2.5 3.0 I
y/h = 0.54
0 I 1.5 2.0 3.0
y/h = 0.18
-A-- la O 11.b 1.5I 2.0 I I 2.5I 3.0I x’/h Fig. 1.2.8 - B. 4- 16
THE COMBUSTION EFFICIENCY FOR TWO INITIAL LEVELS OF TURBULENT VISCOSITY
I I
THE FUEL - AIR EQUIVALENCE RATIO 0 = 1.66
Fig. 1.3.1. 4- 17
c I- h 4-19
COMBINED NOZZLE SMALL SUPERSONIC FLIGHT MACH NUMBER 1
Fig. 1.4.1. MACH NUMBERS CONTOURS
COMBINED NOZZLE SMALL SUPERSONIC FLIGHT MACH NUMBER
Mf=2.0 BASE PRESSURE
Fig. 1.4.2. MACH NUMBERS CONTOURS \ 4-20
HYPERSONIC NOZZLE
Mf= 10
Fig. 1.4.3. MACH NUMBERS CONTOURS L 1 I 1 I 1 I l
NOZZLE PERFORMANCES THRUST LUSSES
A%
0 Fig. 1.4.4.
0 5 10 Yf THRUST VECTOR ANGLE
Fig. 1.4.5.
qn 0 5' IV 5-1
0 RESEARCH AND DEVELOPMENT OF RAMJETS/RARIROCKETS PART 11 INTEGRAL LIQUID FUEL RAMJETS by Prof., D.Sc. V.SOSOUNOV ClAM (Central Iiistitiite of Aviatioii Motors) 2, Avianiororn;iy;i St. I I1250 MOSCOW RUSSIA -SUMMARY LR - liquid rocker; RI - ramjet engine; 1iiregr;iI liquid fuel ramjets have some SPRR - solid propelli~~itramrocket; - specific feiit 11res i nfliie iicing on the i I' ii rrii iigenient SR - solid rocket. and coiistriiction. Aiiioiig siich features are: mostly bigger size and working diiratioii than by SPRR, rhe necessity of fuel manifolds aiid fliimeholders, ;,lid illso Will1 cooli11g svstem I. INTEGRATION OF LIQUID FUEL ii rrii iige inent i 11 combiist ion c h;i mbe I', f he RAh1.I ETS WITH BOOSTER possibility of fuel flow and iiozzle posilion 0 control driring flight on different speed/altitucle The coiicepts of' the missiles ii1tegr;ilion trajectories. In this lecture there will be disciissed with iiquid fuel rill11jeIs (LFRJ) ;in6 SPRR ;IS some iteins, concerning tlie mentioned LFRJ propulsion systems have their features, specific to feat {ires: the conditiolis of iisiiig each missile type. As - The ilitegriUiol1 of riirn combustor atid known by practice, the solid propellant ram- booster. rocket engines are the most preferable on sinall - Duct flow iiist;ibility during tlie ejection of "air-to-air" aiid "~I'oriiid-to-grouiid"missiles, the booster case. liqiiid fuel raiinjets are preferable on I;irger "air- - Fuel siipply device and effective coinbiistion to-ground" iiiid "eroillid-to-groiiiid" missiles. Tlie in ramjet. main characteristic sections of the iiitegrxl - Adaptive control of rainjet 011 different flight ranijets (solid or liquid fiiels) :ire piveii iii Fig. trajectories . 3.1. N O M EN C LATU RE The relative low gross rake-off weight ;tiid low operation times of "air-to-air" and "ground- A - i1re;i; to air" missiles predetermine using on them the d - diameter; SPRR with fixed g;is-air flow passage iis a G - inass flow rate; 'propiilsion system of tlie second missile stage. H - altitude; The cast or rainjet case-well-bonded graiii of 21 - Mach iirimber; solid propell;int booster \till1 ejectiible iiozzlr (or OMq = p,,V,,'/2 - dynarnic pressure; nozzleless solid booster) is iised to till ;it lilosl tlle r - rdnge; voliir~ieiii the SPRR. T - te lnpe r;it 11 re; t - time; Fig. 3.2 shows tlie principle scheriies 01' the V - flight velocity; integral ramjets of different class and appliciitio11. a" - angle of attack; h - dimensionless velocity coefficient; Opposed to the missiles with the SR the P - density. iiitegnition of the inissiles with ramjets is ;Issociiited with solving some probleins caused, INDICES first of all, by Iiigh strrictur;il coinplexity of separ;ite elements and components of the f - fuel; combustor and jet nozzle iis well as of the file1 and vehicle ;ictwfor coiirrol sysrenis. 0 - ii~iibie~it i~tniosphere:
The type tlir stait (iiir, grolliid ;iiid AB B R EV I AT1 ON S of subiiiiirine) defiiirs ;I booster gritin iiiass. which AB - after kiiriier; may highly differ, due to ;I significant difference ABE - air breat hi iig engi nes; in time from the missile start to the rainjet LFRJ - liquid fuel ramjet; ope r;i t ion ini t iat i o 11. LH2 - liquid hydrogen;
Presented at an AGARD Lecture Series on 'Researchand Development of RadScramjets and Turboramjets in Russia; December 1993 to January 1994. 5-2
Wlieii designing the ramjets for "air-to- The individual siispeiisioii under ii fiiselage, grouiid" siii;ill missiles (d = 350-450 mm) clue to as ;I rule, does not limit the choice of intake type. 0 the severe restrictioris III ii free voltime, they itre The missiles with veiitrdl intake :ire the most iised, where possible, the simplest solutiolis iii rational for cassette lar~iicliers of ;I drrlni type; term of the striictiire: the organization of itsing the iiose air illtake is the most reasoii;ible biirning process stiibillzed 011 21 bottolil without for missiles of siibiiiarine start. flariieholders, ii fixed frontal unit at the comblistor inlet with ;I fixed supersonic nozzle Fig. 3.3 illiistr;ites the sclieine of the (the missile of Kh-3 I type) and using, where i1itegr;il liqiiicl firel ranijet with the insetled solid possible, t1ierrii;il barriers of the cornbitstor walls propell;iiit booster. Fig. 3.3 shows also :ill possible - (of ASALM type) instead of air-film oiies (see types of intakes tliitt ;ire used 011 existing and - Fig. 3.2). develop;ible niissiles.
The other coiicept of the integration - the Fig. 3.4 shows loiigitudinal section of the - location of the solid propellant booster together integral liqiiid fuel ramjet, iii which the frontal with ii cxe in the liquid fitel ramjet combiistor is itnit and flaineliolders are unfolded after booster used for the high gross rake-off weights ;tiid high ejection and the strpersoiiic nozzle thro;it section engine operation time of "gro~ind-to-Eroiii~d"and is vi\rii\ble. "air-to-prorind" missiles. 111 this case the conip;ict 1oc;itioii of the eiigine eleinents, cornponetits aiid The possible scheines of iiitegr;iiioii of the 11 11its CoITlpl icated i 11 tlie strllct 11r;i I Ilia 1111 fiictll ri lig I'il I iijet ;I lid sol id pr~peI 1; I 11t booste I., prese 11 red in a (unfolding fuel rn;3nifolds and flameholders, Fig. 3.1, 3.2, and 3.3, to ;I considerable exteiit ;ire nozzle tlap coiitrol mec1i;iiiism) in the limited associated with the size of the vehicle and engine, volitnie of combustor is needed to fill the depeiid on, the velocity, tip to which a solid coriibmtor volrinie by inserted solid propellaiit prope I I ;i I it boost er ope rates , t I1 e st i I rt i ng booster the most fiilly. c o 11 d it i o tis, et c .
Usi i ig t I1 e i I ise rt ed so I id p rope I la I it boost e r So, [lie problem of solid propellaiit booster has its ;idv;int:iges atid rl. IIows: striicture mtcliiiig witli crirising ramjet is ;I key - independently to develop ii booster and ramjet question in designing the whole of propulsion on the speci;ilized test beds; system arid nieaiis a re;.isoii;ible comproiiiise iii - to we ;i light ramjet combustor with ;iir-filiii tlie selection of tlie contigiir;itioii ancl the cooling, the most reasonable for high operatioii vehicle. time; - to iise the s;.ime niodiile of the ranijet 011 missiles of different application with sigriificaiitly 2. INTEGRATION OF LFRJ WITH BOOSTER different flight paths and different solid AND DUCT FLOW INSTABILITY propell:int booster enpilies such as in ASM-MSS ni issi I e of cl o 11 ble ;I ppl ica t i 011 . Unsteady processes, such as ejection, inlet start - iiiistart, combustor boosterstiirt - a A wide raiige of iisiiig siich iiii engine for titist;irt, change of hiel ilijection, self - excited fliplit altittide and velocity is provided by also ;I oscillation iii iiilet or coiiibwtor, rainjet resporise stepped or smooth control of the nozzle thro;it 011 external sliock w;ives and so 011, ;ire very ;ire;> after the solid propelhiit booster ejection ;is i iriporta lit p;i rt of i 11t eg r;i I ra nijet s fii I ic t io11 i lip, well as by a wide range of fuel flow control that because they ciiii lead to the engine distructioii. is possible with liquid fuel. Nuinericril simiil;ition is the only way of The type of launchers (or the loc;ition of iiivestig;itions of this processes, as a rule, because missiles oil il currier), as well as I;iunchiiig experimental resem-ches are very expensive, need conditiolis iiitlriences directly the iritegr;ition of ;I perfect equipmelit, aiid ciiii lead to experirneiital vehicle with ;I ramjet. Both these filctors define, niodel destruction. first of all, the type of the air ilitiike ;itid its placement on :.I missile. At present. iis known, the Tlierefore, the iiiathematical models and missiles with air intakes of different types are liIlIi1ericill codes for ui1steady processes developed: nose (frontal), reiiiovable (two- siinul;itioii were worked out in CIAM. This dimeiisioiial, circiilar, segmented) and veiitr:il. iiiodels based oii I D :iiid 2D iiiisre;idy Euler Their location is deteriniiied by the type of ;I equations for rextiiig @;IS inixture. They iriclude la II tic he r . vririoiis experiinent;il ;iiid empiric;il d;ites, illid take into accoiint different constrriction peculiarities of researched ramjet. Niirnerical 5-3
codes ;ire based 011 tlie second order accriracy part, picked out in Fig. 3.10, is possible the God1111ov- Kolgati iilgol.itlim. booster ejection.
Now, let me present you some resiilts of Tlie combustor start. In Fig. 3. I I yo11 can unsteady processes numerical sirniil;ition, see some results of the next unsteady process fulfilled in CIAM. 111 Fig. 3.5 yo11 can see tlie iiumerical simulation - tlie combustor start. scheme of Liqiiid Fuel Ramjet with the ejected Mathematical model, which we use, bases on booster, used as tlie first stage of missile. All unsteady Euler eqi1;ition for reacting gas mixture. investigation, I'll tell you, were fiilfilled for srlcli We siippose, that there exist an equilibrium state type of r;irrijets. ill siibsonic combustioil chamber. We take into account experiment:il tiiites for the leiigtli of Gas oscillation in closed inlet during biirning zone, coefficient of cornbitstion missile acceleration. On tlie first part of missile completeness ;tiid so on. We find tlie dyliamic - trajectory tlie ramjet's inlet is closed by tlie loiids, tlie m;ixi~num reniperatiire in conibustor booster. If flight Mach number is greater then 1.5 and other parameters of uiisteady flow. The main - 2.0, tlien selfexcited oscillations in ramjet's inlet and very important effect, we've found, is follow. appear very often. They can lead to ramjet The initiatioii of air-fuel mixture leads to distrrictioii. We use numerical simulation to ;ippe;.ir;iiice unsteady stress waves, that move predict tlie ;ippe;irance of selfexcited oscillations, from combustor to runijet inlet atid decreases to c;ilciil;ite their iimplitiides and freqiiencies, to diiriiig short time tlie total airflow through riiinjet 0 predict tlie efficacy of different ways of this diict. In Fig. 3.1 I yoii can see the the evolutioii osci1l;rtion atnplitiide decreasing. In Fig. 3.6 yoii of air excess coefficient. Curve (1) is quasi steady see t 11 e e X;I I ii p I e of s I Ic 11 I 111 me ri c ii 1 si m I I I ;* t i o 11. approxim;ition, when totiil airflow is co1i~t;111t;illd Flight Mach iiiiriiber eq1Iiils 2.5, tlie Inaximal only fuel weight flow changes. Curve (2) undimensional amplitude is about 0.4. describes a,real 1lIlSte'iidy process, curve (3) is a stable combiistioii boiinchry. Booster ejection. Tlie next problem, we deal with, is the booster ejection - see fig 3.7. You see, that iinsteady decrease of air- Miithematical model, we use, describes booster excess coefficient cm be very strong (from steady motion atid gas dynamic flow in ramjet's duct value 1.8 to minimal value about 0.9). Therefore, toget her. Ae rodyiia mic and inertia I forces, if tlie optirnal value of air excess coefficient 011 friction, estern;il airflow around ejected booster tlie start part of missile's trajectory is about atid other effects are taken into iiccoiiiit. Yoit see a= I .2- I .3 and the stable combiistion boundary is in Fig. 3.7 tlie time evolritioii of booster velocity aboiit rx=0,9, the ;tiitoiiiatic control system mist and its location. So we caii predict the possibility provide combustor start with CL=I .8, and only of blow interaction between the booster and after that incre:ises slowly the weight file1 flow to ramjet nozzle. necess;iry valiie cx= 1.2- I .3. And diiring conibiistor restart regime, ;iiitom;itic control system work In Fig. 3.8 the time evolution of self - miist be org;iiiized by tlie same way - see Fig. 0 excited oscillation during tlie booster ejection is 3.11. presented. And in Fig. 3.9 yoti caii see the comp;irisoii of experimetital (curve I) and Self-excited oscillijtioll in combustor. In niiinerical (curve 2) siinulatioii of gas dy1iiirllic Fig. 3.12 yo11 can see the ni~mericalsimu1;ition flow, iiidiiced by the moving booster. You see 21 re<s of -self-excited oscillatioil iii combiistor, good iigreenietit between numerical and working 011 ne;-ir-stoicliiometric air-fuel mixtit re. experimeiit;il results. So, we can c:tlciilate the Tlie physical tiiltllre of this oscillation is follow. values of dynmnic loads, iipplied to the internal Lets suppose, that we increase fuel iiijection so, surface of ramjet and initiated by moving shocks tliat air-excess coefficient decreases to the value and supersonic zones. It's a very irnportanr result, Q = I .O. The11 the combustor heat 111~reit~etoo. becarise the uiisteridy aerodynamic lo;~dscaii teiir At the same moment tlie stress wave appearatices off the defense screens in combustor. in combustor atid moves to tlie ramjet inlet, decreasing tlie total airflow. Therefore tlie value All this results can be summarized as a safe of becomes less than 1.0, for example c1 = 0.8. ejection area, presented in the space of CL combustor heat decreases, tlie stress wave p;irameters: flight Mncli number, flight altitiide in T1i;in becomes weaker, the total airflow increases for Fig. 3.10. Cirrve 1 is ii boiindary of blow intefiictioti between booster and ramjet iiozzle; tlie shoit time to value Q = 1.2. But very soon ciirve 2 is the bouiidary of combitstor destruction the constant vriliie of fuel injection leads air- '.d '.d lie to high pressure in front of booster, ciirve 3 - the missile trajectory. The missile's trdjectory 5-4 excess coefficient to it's steady villue a= 1 .O, and and is equal to L, = (10-12)B, where B is an all process. described above, repeats. average distance between the axes of two adjacent flameholders. The given scheme of the Ranijet response on external shock wave combiistion process orgiiiiizatioti hiis proved itself and iitmospheric nonuniforinities. In Fig. 3.13 and at present is being used widely in turbofan you can see the result of numerical simulation of afterburners with flow mixing. Fig. 3.14 shows the shock wave motion through the ramjet. The typical ;\rrangemeiit of fuel manifold aiid flame moving shock iiicreases the total airflow throiigh st ii b i 1 i ze rs for the t II rboj e t i+ ii d t 11rbo fa ii the rilmjet duct. Therefore air excess coefficient cy afterbii riie rs. increases too. New value of (x must be less ttiaii Using the mentioned scheme for the s t ii b I e c om b tis t io ii bo 11nd a ry . combustion process orgaiuzation in tlie ramjet All this niinierical researches were fillfilled combustors led, as a rule, to sigiuficant problems in tlie ClAM and were used for developnieiit of with providing the stable operation of these different Russia 11 ramjets. combustors at low inlet air temperatiire (T,s SOOK). These problems were associated with the fact thiit the combustors of the meiitioned 3. FEATURES OF THE COMBUSTION scheme of the combiistioii process org;iiiization PROCESS ORGANIZATION IN RAMJET show a great fendelicy to the excitatioii of the regiil;ir longitiidinal gas pressure tluctuatioii in The scheme selection of the combiistioii them, particulwrly in filet reach mixtures. process organization in tlie straight-throiigli-flo~~ type conibiistors (ram combiistors) is mostly Studyuig the coinbiistion process of the defined by the flow parameters ;it ;I combustor atomized liquid fuel at low flow tempertitiires inlet such iis pressure P,, temperatiire TI aiid h - showed that a decrease in combustion process velocity coeffic ieiit . sr;ibility WiiS due mostly to a reduction iii a degree of the fuel vaporization. ahead of A possible range of change in these flameliolders, tlie deposition of filet droplets 011 parameters in modern afterbiirners and ramjet the flameholder walls arid, hence, a significant combustors is tabulate in Tdble I. enrichment of tlie gas with fuel, being in the recirciihtioii zones after tlieni 131. Table I. Combiistor type PI T, h Fig. 3. I5 gives the experiiiieiit;il k Pa K dependence of rhe air excess coefficient iii the flameholder recirculation zone on the fuel Turbojet 40-600 900- 1200 0.20-0.25 v;iporiziitioli degree z upstream the flameholder. iifterburlier It follows from the presented data that at z + 0 Tu rbo fa Ii 40-600 600-850 0.15-0.22 tlie giis composition in the recirciilation zone is afterburner with approximately 2.2 as rich as the mixture flow mixing ;ippro;iching to the stabilizer. Ramjet combristor40-600 350- I300 0.20-0.22 At high ;iir teniperiit1Ire (TI>I 100- I2OOK) It follows from the Tiible that feature of the the t'iiel self ignition IS posslble upstream booster-siistainer ramjet combiistors is ii wide flaiiieholders. The biirri-out and destriictioii of rmige of variations in a combustor inlet air flameholders in this case ;ire inevitably. temperatiire which covers a possible change of t ti e t e ni pe rii t 11 re i ii a I I exist i iig a fr e rbu rii e rs . Therefore, for the ramjet cornbustors operating at the low inlet air temperatiire CIAM, Whe 11 the t 11rbojet afterburners we re suggested the other scheme of combustion creating the combustion process organizatioti process organization, so-called scheme with tlie scheme wit ti coiisec 11t ive realjza t ion of filet separate fuel siipply, in which the processes of mixing, vaporization and combustion processes fuel mixture, vaporizatioii and buriiiiig take place along the combustor W;IS developed. In such ii simiilf;ineously. scheme of combustion process orgiiiiizarioii the fuel is supplied to flow upstream flnmeholders iii The proposed scheme of tlie coiiibiisrioii a distance siifficient to its practically full process o rga 11 i zi i f io i 1 w i t 11 se pi3 rat e fi i e I s II p p I y vaporization and mixing with gases elitering to i+llOWS to provide: tlie afterbiiriier. Under coiiditioiu typical to - a stable combustor operation at low air inlet afterburners the fuel burning length L, depends temperature (T, 5 SOOK); e niaiiily 011 mritiial ;irr;iiigement of fitel stabilizers 5-5
- ii high reliability of the combiistor oper;ition Two coinbiisrion chambers witli [lie siiiie 0 at high air inlet temperatiires (T, > 900K); relative length of the hor secrioii (L,/H = coiist) - inaxirniim possible sriibility of tlie conibirstion and approximately the siinie sclienie of the iiiter process under high altitude conditions; combiistor process control were selecred ro - the maximum simplicity of the structiire when compare. developing the combiistor with unfolding - flameholders. Fig. 3.18 and 3.19 illustrare tlie regions of the stable burning in rliese combustion chambers. However, there are some problems in It follows from the data presented in Fig. 3.18 - providing high combustion efficiency at low air that at T ,,,, = 770-820 I<, h = 0.15 and P,/Po = - temperatures and, besides, the special 0.5 the burning in afterburner is srable iii the requirements are miide on consiimption range of the air excess coefficient variation of 0.9 characteristics of fuel manifolds at such sclieme I a 5 6.7. The minimum value of pressiire in - of the combustioil process organization. afterbiirner, ;it which the fuel cornbustion is possible under rliese conditions is PJP,, = 0.27 (a Schemes of Combustors of Liquid Fuel Ramjet = 2.0) 141. Integrated with Booster Engine Decrease in afrerbiiriier iiilet temper;iture Let 11s colisider now the possible schenies to T, = 500 K leads to the facr rliar ;it P,/P,, = of the flameholder mounting in the air cooled 0.5 the cornbustion is stable only at a > I .25 and combustors of booster-siisrainer ramjet. 1 ti the case when the length of the booster engine, that the minimum pressure at this temperitiire is is placed into combustor, is less thijn the extent P,/P,, = 0.32 and correspoiids to rx = 3.0. of fuel burning zone length or equal to it, rlie flame holders a nd fuel manifolds are st ;it ioiia ry The ' rarher coriipliciited ciir-ve, boiiiiding located in the combiistor as in the rurbofiin the region of rhe stable biiriiing in the ;ifrerbiiriier afterburners. with obviously marked rniiiimiim, is traiisforined for tlie ramjet cornbiisror into a straight line above which the burning is stable and below Iii case when the length of tlie booster engine exceeds the extent of the burn-orit zone, which it is impossible (see Fig. 3.19). the iinfolding (turning) flameholders are needed to install. This is associated with the fact thnt the At the smie v;iliies of T, arid A the enlargemelit in combustor length beyond required miiiiinum gas pressure when the combustion is for the fiiel burning out leiids to the increiise in still possible is miicli the same iii both air flow to its cooling system and in this combiistion cliambers. The difference is in fact connection to a reducrion in fuel combiisrioii rht itlider the minimum pas pressure condition efficiency in the fuel rich region. the brirning iii the afterburner is possible at priictic;illy the s;iiiie v;il\ie a, whereas in rhe Fig. 3.16 presents the schematic di;igr;im of rainjet combustor it is possible in a wide range of srich a combustor. At the initial morneiit of time a variatioii. the flameholders are located close to the combustor wall as the booster engine occupies Fig. 3.20 compares the combustion practically the whole volume of the combustor. efficiency for two considered schemes of the After its jettisori the flameholders are iinfolded combiistion process organization in combiist ion (turned) itlid are set into operiting position. After chambers with close relative length of hot section that the runijet combustor is starred. Scheme of (L,,/B) 151. It follows from the diagram that it' the siicli unfolding radial flameholder is sliowii in combiisror inlet air tempemtiire lies in the range Fig. 3.17. typical to the afterbiiriiers, the conhiistion efficiency is pracrically the same in both schemes. Characteristics of Burning Process in LFRJ Combiistors Increase in flow temperatiire lip to T, = 1170 K in the ramjet combusror slightly To compare the characteristics of the fuel influeiiced the combustion efficiency valiie in it. combustion process in the modern acceleritiiy- Decrease in ramjet combustor iiiler air sustainiiig ramjet combiistor and turbof:in temperitiire to T, = 360 K with a simiilraiieoiis afterburner with flow mixing is of interest,, when reduction in fuel temperature to rl. = -50°C led to the cornbiistioii process in tlie latter is organized the fact that the maximum combustion efficiency by the traditional scheme. in it did nor exceed the value 11 = 0.68. 5-6
So, the we of the developed scheme of the atid r;ither cornp1ic:ited flight paths (see Fig. combustion process 0t.gatiizatioti with separation 3.22), over the course of which the flight illtitilde fit el sit pply allows to create the booste r-siistai ne r and speed are v;tried in i\ wide range. It reqiiires ramjet combustor operating in more wide ratige the change iri wide r;inges of the engine tlirust of inlet air temperature in comparison with ;lnd fitel consittiiption delivered to it 011 different combustors designed according to triiditioti;il segments of tliglit path. As the etipirie thritst scheme. depends 011 the ;ittnospheric conditions, they initst be tilketi itito cotisideration in controlling Features of the LFRJ Cornbitstor Start the fitel cotisitrnptioti tlirough the eiigiiie. In connection with tlie fact that LFRJ iire used, as In coticlusion it shoiild be said several it rule, on tttlrniltitled aircraft the engine ACS words ilbottt the features of the combustor must be algorithmic;tlly integrated with the starting in the moderii booster-sustainer ramjet. :iircriift flight control system to provide the required thrust on all the flight path segments. If tlie ramjet has a frilly-variable nozzle, In this case both ititer-etigine parilmeters ;itid the combitstor may be started under tlie sitlne also the piirilnieters used in the aircraft tlight conditioiis iis the afterburner of modern turbofatis control system are ittilized to cotitrol the fitel and priicticiilly by the s;ime meatis. If the riiriijet corisitmption tliroiipli tlie etigine [ 21. is equipped with the nozzle, having ;I cotist;lnt throat il rea, or two-position nozzle, the The necessity of r;itional itsing all the combustor starting rinder these conditions volumes 011 itiim;itiiied aircraft leads to the becomes a more complicated problem. advisability not only of algorithmic but also 1i;irdw;ire integr;ttion of the engine ACS with the The point is that in this case, after the flight control system under which all the booster engine jettison illid unfold of the ciilc\tliitioi1s are performed by the s;inie onboard flameholders in their operiiti~igposition, such it computer aind only the transducers arid actitators flow regime in the cornbitstor is set when the ;ire itistalled iti the enpitie. flow velocity between the flameholders becomes equal to the sonic one (A = 1 .O). The coinbitstion Target of LFRJ Cotitrol process stabilizatioll with sitch velocities of tlie flow in the combitstor, using hydrocarbori fiiels, It1 selectitig the cotitrol parameters atid is impossible. progr;tm, in addition to the provisioti of the required engine thritst on all the flight path Therefore, the special starting devices of high segments, the main task is also provision of the heat power with rather high time of operdtion are best flight characteristics of aircraft 011 a111 check used to start the rmijet combitstor. In this case flight path, narnely, tlie most effective use of after switching 011 the starting device, firstly, the sustain fuel reserve. Besides, the LFRJ control gas flow iii the combitstor is transformed to the tasks involve the litnitiitioii of the engine siibsonic one due to heat release in flow as ;I operation regimes atid iiircraft flight, providing result of the fitel portion combustion, delivered to the conditioiis for no-failure operation of rhe the combustor, in hot jet of the starting device. main engine atid aircraft compotietits. Among Only thereafter the coinbustor is ignited. s I t c 11 restrict i o t 1 ;I re : for the aircraft - the niinimiini rani determiiiing, the possibilities of aircraft ma ne wer; 4. CONTROL OF LIQUID FUEL RAMJET - the maximitin rim determining the loads on (LFRJ) aircraft cornpotients and elements; - the maxiinitin air stagtiation temperatiire Solid propell;int ramrockets have limited defining the materid strength of the aircraft possibilities of iti-flight control. 111 cornp;-lrison st ruct I I re; with thein liquid fuel ramjets can be widely for the engine controlled. Control stages of SPRR in - the provision of adequate tiiiirgitls of fuel com pa riso ti with LFRJ ;ire show ti schema t ic;i I1 y comb I I st i o I i st :i b i Iit y ; iii Fig. 3.21. - the provisioti of adeqitate margiris 01' illtitke ope rii t io ti st il bi Ii t y ( s t i~ I I tiu rgi 11) ; The automstic control system (ACS) is one - the prevetitioti of reduction in tlight Mach of the basic components of LFRJ oti nurnber below the minimitm Mach nuniber under maniifiictitring quality of which the flight which the flight becomes impossible at ;ill values performance of aircraft depends. It is necessary of fitel consuinptiori and air excess coefficient. to provide the LFRJ missile flight on different 5-7
is reqiiired to provide ;I wide raiige 01' vilri;l~i~>Ii @ Tlie LFRJ control tasks involve iilso tlie i ii tl ig lit M iic 11 1111 I n be r. provision of its reliable starting. Tlie coiirrol loops of the nozzle and illtake The nientioned restrictions of tlie engine inay be made both with uninterrupted change of operiition regime aiid vehicle flight deteriiiine tlie its geometry iilid iiS two-position OIieS - with - boiiiidaries of the latter application, areti in tlie relay cliange of the iiitiike aiid nozzle geonietry. M-H coordin;ites. Fig. 3.23 depicts the Iii the first case the aircfiift flight c1i;iracteristics is a p p I i cii t i o ii ;i rea of t li e "ii i r - t o -si I rfii c e" in i ss i I e iIS only little more t ha 11 the a irc raft c Iia riicte rist ics . an example. The left bottom of tlie ;ire;\ with relay control of tlie nozzle and intake. At .- boiind;iry of ";iir-surfilce" missile app1ic;ition is the siime tinie in tlie second cwse the gkis defined by tlie minimrirn Mach number. The left dyiiaiiiics forces inay be used to move flops of
upper part of tlie application area ' boundary is the iiozzle and illtiike tll;lt sigiiificantly decreases - determined by tlie minimum ram, surge tlie cost, dimensions and weight of tlie control boundary of the intake, combustion instability system aiid makes it possible to place some and again by surge boundary of tlie intake. Tlie additional quantity of the sustain fuel 011 board right upper part of the ijpplicatioll area boundary aiid to improve tlie aircrdft flight c'har;icteristics. is defined by niaxinially ;idmissible stagnation teniperattire of air flow and tlie right lower p;wt Using the control loops of the nozzle and of the boundary - by tlie maximally permissible int;ike is more effective, tlie larger the aircraft 0 values of air ram. sizes (the more tlie fuel reserve) and tlie wider rlie rtiiige of v;iri;ition in ;iircr;ift altitride and The botiiid;iries of the vehicle app1ic;iIioii speed. area depend on tlie parameters of tlie vehicle illld engine ;is well as tlie engine control progr;rni Fig. 3.25 gives ;IS ai1 ex;irnple the effect of atmospheric conditions. tlie nozzle control iiirroductiorl on the flight cliarxterist ics of the "ii i r-sii rface" missile. Curve The flight control program is selected so I shows tlie missile flight range - la\incli altitude that Macli number - flight altitude dependencies dependence only with the fuel coiisiimption (operating modes lines) were not beyond tlie coiitrol according to tlie program boundaries of tlie vehicle appliciitiori are;i during the whole fliglit for all possible tr;ijectories ;ind atmospheric conditions (see Fig. 3.23). In this case the conditions of no-failure operation of all Curve 2 iii Fig. 3.25 illustrates tlie similar tlie basic components of the eiigine and vehicle dependence will) tlie fiiel consuinptioii coiitrol are provided. iiiid two-position coiitrol of the nozzle according to the progr~lnl Structure ;ind Control Loops of ACS
0 Fig. 3.24 shows diagram-scheme of the ramjet ACS. Tlie control system generally involves the control loops of fuel consumption, nozzle thriist area rind intake geometry. Curve 3 in Fig. 3.2s displays tlie missile tliglit range will1 tlie file1 consiimptioil control The fuel consumption control loop ;iiid two-positioii control of the nozzle iiiid iiitake performs all the above mentioned fiiiictions. Tlie ;iccording to pr'ogriiin fuel consumption control loop is made ;is liydraiilic fuel meteriiig iiiijt controlled by electrical signals from tlie aircraft compiiter.
Tlie control loop of tlie nozzle throat area is iised in the control system to match tlie intake and eiigiiie operatioil regimes ;IS well ;IS to provide the adeqii;ite milrgilis of tlie intake It is seen that tlie flight range, depending stability iIt high thrust levels necessary on aircrxfr on ;I 1;iunch altitude, wilh tlie introduction of the climb ;iiid acceleration regimes. two-position nozzle control increases by 4048% and with tlie introdiiction of the intake control The control loop of tlie intake geometry is loop - by I3 - 1776 more. used for improving tlie engine file1 efficiency if it 5-8
To provide high flight aircraft control flexibility the LFRJ-missiles h;is bigger characteristics with sufficiently reliable operiitioii mean tr:ijectories speed then the SR missiles. 0 of ;ill the m;iiii airframe and engiiie coniponents, it is iiecessiiry very c;irefully to select the As the ratlier new propulsioii the liquid fuel coin m;i lid 1x1ra nie t ers a lid to opt i iii ize t lie e ngi iie ramjets this system Ii;is many possibilities of control progriim. The engine control progrnni is future developnieiit: fuel eiiergetic, more effective to be optimized iii close miitual trade-off with control system modes, iiew materials and light - the selection of the aircraft flight control structtire, improved exteriial gas dynamics. program. Developed for this piirpose are the charxteristics calculation program complex of Many successfiilly acting in service LFRJ- - the engirie and aircraft iIS a whole and also ;I missiles show that advaiitages. nriniber of service arid auxiliary progranis, fxilita ti iig sigii i ficii 11t ly and reduc i iig process and ACKNOWLEDGMENTS time of the control program optimizatioii respectively. The author is iridebted to ClAM colle;igues for the available iiiateri:ils and for help of its Diiriiig the optimization of the iiircraft prep;iration: Dr. M. Miklieenko (I), Dr. Sc. F. flight and engine control programs the coinmarid Slobodkiiia ;ind Dr. Egorushkiii (2), Dr. A. signals ACS parameters are selected in such a Kiidijavtsev (3), Dr. Yu. Kulikov (4) and ;also way as to provide rather wide area of aircraft Mr. V. Grachev, Mrs. L. Zhemuraiiova, Mrs. M. appliciitioii, covering all the required regimes of Sapronova for their permanent assistance. e its flight, incliidiiig the most ecoiiomic ones aiid at the same time to obtniii siich thrust-ecoiiomic REFERENCES engine cliaracteristics which permit the flight to be made ;ilong the a11 check paths (see Fig. 3.23).
The complex selection method of the program of the ramjet and vehicle flight coiitrol permits all the possibilities of the engine aiid vehicle to be used aiid so to improve the flight characteristics of vehicles with LFRJ. Besides, the developed method makes it possible to select such directions of works on the main engine and vehicle componeiits that will lead to the maximum improvement of the flight characteristics of vehicles with LFRJ. ' The calculatioiis show th;it the flight ch;i racteristics of the moderii vehicles with LFRJ may be iiiiproved by a fiictor of 2-3 ;is a result optimizing the charxteristics of the main engine and vehicle compoiieiits ;is well as the progrini of the eiigiiie non pen. C.M. WJIHXT~HKO.M., and flight control. MaUlllHOcTp~e~l~ie,1975.
CONCLUDlNG REMARKS
The integral liquid fuel rilljets are the most promisiiig propulsion for missiles of quite big size and long distance flight became of their excliisive energetic effectiveness and comp;rctness. In compnrisoii with SR missiles, which thriist valiies are permanent during the flight, the thrust of liqiiid fuel niinjets is depend 011 airflow, i.e. velocity and altitude, and 011 fiiel, i.e. control system illid 011 intake/nozzle control, and c;iii therefore change in wide range tinder required flight coiiditions.
So LFRJ-missiles are very effective 011 low- altitiide trajectories, on high iIltitiide and 011 trajectories of combined type. Because of thriisr 5-9
n*
I I 0 in din bin bex nexI a Fig. 3.1. Liquid fuel ramjet with inserted booster. Base crossections.
Kh-31 - type
.WLM - type
ASM-MSS - type
Fig. 3.2. Typical layouts of different classes missiles with integral liquid fuel rdinjets. 8 5-10
Inlets type circular nose separated underfuselage air inlet air inlets air inlet
Integration problem
0 SR ejection 0 folding flame stabilizer 0 choise of air inlet type 0 variable node
Fig. 3.3. The LFRJ and booster integration mail1 feattires.
Fig. 3.4. The view of LFRJ combustion chamber (the folding fuel manifolds and stabilizators, and variable nozzle are seen). 5-11
Fig. 3.5. Liquid fuel ramjet with separated solid fuel engine.
tt Fig. 3.6. Self-excited oscillation in ramjet's inlet, closed by solid fuel engine. 5-12
EY
X 1 YC
t -.--c t- Fig. 3.7. Numerical simulation of stage separation. Time evolution of solid fuel engine's velocity. 5-13
c = t, t ='t2
t - Fig. 3.8. The evolution of self-excited oscillatio in rdmjet inlet during the stage separation.
t = t, t=tt
lI- 1 .----.I a \ a" t
t Fig. 3.9. Experimental (1) and numerical (2) simulation of gas dynamic flo in ramjet duct, induced by the stage separation. 5-14
-1
Fig. 3.10. The stage separation area in the space of parameters: flight mach number, flight altitude
ciirve ( 1) - boundary of blow interaction between solid fuel engine and ramjet's nozzle; curve (2) - boundary of combustor destruction, due to high pressure in front of solid fuel engine; curve (3) - missile's trajectory.
I 5-15
I 8) t Fig. 3.11. Combustor start numerical sit11llli~tiol1s
a) Time evolution of air excess coefficient
curve (1) - quasisteady approximation; curve (2) - unsteady process; ciirve (3) - stable combustion boundary;
b) Fuel weight flow along the missile trajectory
curve (1) - optimal regime; curve (2) - the fuel weight flow, governed by automatic coiitrol system; curve (3) - cornbuster restart. 5-16
Q
I I I -a 1
Fig. 3.12. Numerical simulation of self-excited in combustor, using near-stoichiometric air-fuel mixture.
a) Time evolrition of air exes coefficient; b) Mass heat combustion dependance 011 air exess coefficient; c) The total airflow change during self-exited oscillation. 5-17
2.4
1.6
- t Fig. 3.13. The air excess coefficient change, indwed by the shock wave, passing throw the ramjet
Fig. 3.14. Scheme of modern tiirbolin afterbiirner. I - flameholders, 2 - nianifolds, 3 - flow inixiiig device. 5-18
0 q5
I\ I
.-.-. -. -. .-.-.-.-
Fig. 3.16. Scheme of combustor of the ramjet integrated with a booster. 3a - before booster ejection, b - after booster ejection; 1 - flameholders, 2 - booster engine. 5-19
A-A
Fig. 3.17. Scheme of an folding radial flameholder with "carburator" 1 - flameholder, 2 - fuel manifolds, 3 - "carburator"
0 I 2 3 4 5 6 7 o( Fig. 3.18. The boundaries of stable combustion region in turbofm afterburner. - T, = 770 - 820 K, h = 0.15; 0 - T, = 500 K, h = 0.14 - 0.15. '0 P"
Stable operation
n n 0 -U U
0- 0-
Fig. 3.1:). The boundaries of stable combustion region in ramjet combustor. 0 - T, = 500 K, h = 0.14 - 0.15; 0 - T, = 360 K, h = 0.15 - 0,16.
Fig. 3.20. Combustion efficiency of modern afterburners and ramjet combustors.
\\\ - afterburners, L,/B = 10.5, T, =600 - 850 K, P, = 70 - 100 kPa, h = 0.15. Ramjet combustor, L,/B = 9.5;
0 - Tt = 990 K, P, = 65 kPa, h = 0.17 - 0.22; - T, = 1170 K, P, = 60 kPa, h = 0.25; 0 - T, = 700 K, P, = 50 - 60 kPa, h = 0.2; A - T, = 360 K, P, = 300 kPa, h = 0.2. 5-21
I IT m Ez P SPRR SPRR/LFRJ LFRJ LFRJ LFRJ
0
Nozzle 0 0 0 @
Intake 0 0 0 0 0 - iincontrolable ;
@ - step-to-step control :
0 - f1iU control ; Fig 3.21. Control stages of SPRR/LFRJ. H
2 Fig. 3.22. Typical flight trajectories of air-surface missile. 5-22
i
Fig. 3.23. Air-siirfxe missile envelope with typical trajectories.
0 - start, * - finish.
MISSILE
Fig. 3.24. Automatic control system c1i;igrarri. 5-23
relative range
Fig. 3.25. Dependence of air-surface missile flight of the ramjet control complication, by different start height.
1. Fuel flow control.
2. Fuel+nozzle control.
3. Fuel+nozzle+inta ke control.
6-1 0 RESEARCH AND DEVELOPMENT OF RAMJETS/RAMROCKETS PART 111. THE STUDY OF GASEOUS HYDROGEN RAM COMBUSTORS by Prof., D. Sc. V.A. SOSOIII~OV, ClAM (Central Institute of Aviatioii Motors), 2, Aviamotorn;iya St., 11 1250 Moscow, Riissia
SUMMARY GH, - gaseoiis hydrogen; GH, T - gas (vapour) hydrogen tiirbine; The peculiarities of H, ramjet turboramjet HE - heat exchanger; high efficiency ram combustors are discussed, LA - Liquid air; which include some special qiiestioiis, such as: LH, - liquid hydrogen; - flame stabilization behind nozzle edge; P - pump; - increase of number of nozzles and ram R - rediiction gearbox; combustor length reduction; T - turbine; - separate combiistion witho~~tmixing of tlames; TF - turbofaii; - uniform spread of fuel aiid ;iir in chamber TJ - t 11 rbojr:t ; sect ion; TR - turborocket eiigiile. - the experiment on model burner series; - methods of effective H2 combustion in short I NTRO DUCTI 0 N combustor. Unique properties of hydrogen ;illow to N 0 M E NC LATURE study aiid develop a great number of new types of eiigiiies with complex thermodynamic cycles A - area; havi iig a si1 bsta lit i;i I I y bet ter c ha rac te rist ics aiid d - diameter; parameters iii comp;irison with the coiiventional G - mass flow; gas tiirbine engines. L - length; M - Mach number: These iiniclue properties of hydrogen e~liible N - number of nozzles; to significantly improve tlie paranieters and P - pressure; characteristics of t lie air-breathing engines: S - distance between nozzles; - to increase tlie specific tliriist (per I kg/s of T - temperature; air flow) by 1.5 - 2.0 times; W - speed; - to increase the specific impulse at maximal p - eqiiiva~encefrie~/air ratio; thriist by 1.2 - 1 .5 times; 6 - nozzle edge thickness; - to increase thriist/weigIit ratio by 1.5 - .2.0 11 - fuel combustion completeness; times; P - density. - to increase the maximum speed use of tlie gas tiirbine engiile by using of cooling capwiry of INDICES LH, till = 5-7.
a - air; These iidvaiitages ;ire gained iii different b - combustor, combustion; types of the LH, - ABE. At preseiit in different f - fuel; countries and compaiiies there is under way an fl - flame; extensive work for searching of tlie optimal types of the LH, - air-brei*thing engines suitable for I1 - iiozzle. the initial acceleration phase of the aerospace AB B R EV I AT1 ON system - the third generation of space vehicles. This significant theme is not a main subject of A - air; th;it report. Therefore one may confine to the AB - afterbiiriier; demonstratioil of some ABE proposed by ;.I ABE - air-breathing eiigine; number of engineers and companies and rising - combrisror; the cooling and energy capacity of the hydrogen e: - compressor; or both these properties (Fig. 4. I). Morphological F - fall; resei3rch by Zwicky I I I shows that such types of FV - flying vehicle; the LH, - ABE can be inore than 10000. Using Presented at an AGARD Lecture Series on ‘Researchand Development of RadScramjets and Turboramjetsin Russia: December 1993 to January 1994. 6-2 tlie liquid oxygen if it is abroad permits 2.STABILIZATION OF THE HYDROGEN additionally to improve thrust and weight of the FLAMES e t ii rbo roc ke t engines. The first problem is a stabilizatioii of tlie All of LH, - ABE, which are developed for hydrogen flames at the nozzles (Fig. 4.3). The high Mach numbers flight vehicle, have reheat or experiments in the open air flow showed that the ram combustor chamber (for example, nozzles must have a siifficient edge thickness to .. turborainjet engine for Zanger). The operational realize the condition of the flame stabilization in behavioiir of H2 combustion chamber is the recirculation zone (RZ) behind the edges, significantly differ from ram cornbiistors, using iiamely tlie time of air-fuel mixtiire staying in it . liquid or solid fueb. must be more than the time of burning. Obtained - in the experimeiit tlie relationslup generalized in 1. MAIN PROBLEMS correlating parameters, shows that stability of the hydrogen diffiision flame is well above than of There are considered ABE designed for methane flame. Received dependeiicies allow to subsonic and supersonic FV. The special select the required area ratio of the fuel nozzles problems for scramjet are not discussed further. and its edges thickness to obtain the stabilization of the hydrogen flame in the required operating Liquid hydrogen in talks as a resiilt of range. cooling of hot engine components or in special heat exchangers mrist be gasified and enters the 3. RELATIVE HYDROGEN FLAME LENGTH combustion cliamber as a gaseous state. There are examined the variations of the The experiments shows that the effective relative hydrogeii flame length Lll/d,, from tlie burning of tlie hydrogen in the prim~y speed ratio of the air to fuel Wa/Wl. in the opeii combustion chambers of the gas turbine engines air flow (Fig. 4.4).When the nozzle edge (ahead of the turbine) presented no special thickness is equal to 0, tlie m:iximum flame problems: it is necessary to increase the air sripply lengtli is observed following tlie theory of G.N. to tlie primary zone of the chamber and to Abraniovich 161 at W,,/W,. = I. With the finite elislire an intensive mixture of the hydrogen with nozzle edge thickness the flame length decrwses the air, for example, with the help of the swirling continually with the air speed increasing that is of these component flows in the injectors and related to more intensive hydrogen burning in the combustor domes. recirciilation zones. The rate of decreasing in the flame length growth with the relative edge A wide range of problems is posed while thickness increasing. burning tlie hydrogen in the reheat and ramjet direct-flow combustion chambers of a small One can obtili~ithe generalized relationship length (Fig. 4.2). Srich combustors are specific to of the relative flame le~igth on the complex reheated turbojet and tiirbofan engines and also parameter involving the speeds and densities ratio to high-speed ramjet, turboramjet and of air and fuel ;It the finite nozzle edge thickiiess turborocket engines with subsonic speed of the (in Fig. 4.4 below, left). air and the gas in tlie ram combi~stors. In such chambers there occiirs the diffusion combustion 4. FUEL NOZZLE NUMBER of the hydrogen flowing out of many nozzles as opposed to ‘he biirning of the homogenized It is known 121 the idea about tlie diffusioii mixtiire of the kerosene with the air in the flame reduction with fuel nozzle niiinber conventional combustion chambers. increxiiig. In this case ;it the same fuel flow the flame length L, 2s against that of a single flame * Research series of the burning of the L,, decreases with hydrogen in such cornbiistors carried out in CIAM by the atlthor together with coll;lbor;irors 12, 3, 4, 51 found out the conditioiis the hydrogen flame stabilization at the nozzles and also the possibilities of nozzle number increase to where N,, - is ;i nozzle number (Fig right). reduce the combustor length in the event of 4.4, isolated combustion without confluence of This dependencies used for designing the flames. While meeting these conditions the reheat combustion chambers. However, nozzles in and ramjet combustioil chamber iising the pylons miist be located so that the flames did nor hydrogen can be m;’lde 2 - times as short ;IS ;I 2.5 interact with each other and did not contlrient combustor rising the kerosene. into one comnion flame (Fig. 4.5). 6-3
The minimal critical distance between operation region = L,,G,./G, = 0.2 - 1.0 by 0 nozzle axes S* was determined by experiment. using the special nozzle devices with great This distance decrease with increasing of the number of nozzles i.e. by long equivaleilt relative relative air speed and the nozzle edge area. combustor chninber (2) the experimentally achieved value of 11 were 0.9 - 0.99 ill ;in 5. HYDROGEN COMBUSTION CHAMBER- absolutely short combustion cliamber. NOZZLE SYSTEMS CONCLUSION 111 Fig. 4.6 are given tlie results of the experimeiir on a model burner liaving some The results of investigation of hydrogen . lateral fuel manifolds located U1 one plane or combustion chamber give 11s tlie methodology of separately. The number of the nozzles in tlie developing a short combustion climbers manifolds is varied as well as a11 inclination of (decreasing of length more than two times) with . their axes. Spacing between tlie nozzles met the st;ible operation in wide region of fitel/;iir ratio, requirements S > S*. The separate disposition of the possibility of developing of CDF-codes and manifolds with an axial fuel feed proved to be design engine with hydrogen ram combustor. best in combustion eficiency and hydraulic resistance (burner 3). REFERENCES
In accordance with the relationship (1) oiie I. F. Zwicky Morphology of Propulsive Power. 0 may introduce the concept of the equivalent Monograplls 011 Morphological Research, N? 1, relative lenpth of multi-nozzle chimber CL. This Society for Morphological Research, Pasadena, Calif., 1962, 382 pp. length resembles the relative one of a single nozzle chamber of the same length L, and 2. V.A. Sosouilov. Some problems concerning diameter dhcq 121 in the combustion efficiency optimal ducted rocket engine with secoiid2iry brirning. Proceedings of rhe 2d ISABE, 1974.
BY means of the parameter L: oiie can compare the chambers with different nozzle number such that Li can be varied both by chamber length L, or the nozzle number N,,.
The combustioi~efficiency 11, derived from the heat combustion Hit, depends on the 0 equivalent cliamber length and fuel-air nitio p. One may select on the basis of the experimental data the type of the function f(p) at which the -* 5. V.A. Sosoilnov, Yl1.M. Annushkin, E.D. parameter Z - L, f(p) will define a unique Sverdlov, D.G. Pagy. Investigation of hydrogen diffusion flames in direct-flow combiistors. fxliion of the combustion efficiency 11. The Hydrogen Energy Progress. VI I Proceedings function f(p) has a different form at p < I (air of the 7th World Hydrogen Energy Conference. excess) and p > 1 (fuel excess). Obtained in such Vol. 3. Pergamon Press, 1988, pp. 2009-2024. a way the generalized relationship of biiriiing completeness of hydrogen (Fig. 4.6) enables to 6. G.N. Abramovich. The theory of tiirbdent jets. estimate by calculation the combustion efficiency Tr;ins. from Russian. The M.I.T. Press, characteristic at a different ntimber of tlie nozzles Cambridge, Massachusetts, 1963. or chamber length, iiamely at 1; = var (Fig. 4.6 on the right). Here is shown the comblistion efficiency 11, determined from the maximally possible heat generation, therefore at p < I: I)* = 11 (left part), and at fi > 1: q* = 11p (riglit part).
These data show the possibility to achieve very high H, combustion efficiency. In the usual 6-4
A
B
v H2 -@
The use of cooling A capacity and ciiergy 0 Fig. 4.1. New types of LH2 air-breathing engines. 1- TJ with precooling; 2 - TJ with intercooling; 3 - liquid rocket engine whit11 air liquefaction; 4 - turborocket engine (TR) with a gaseous hydrogen turbune GH2T; 5 - TJ with GH2T; 6 - TF with the GH2T; 7 - TF with the intercooling and GH,T; 8 - TR with precooling; 9 - TR with air liquefaction.
turbojet 01 combined engine ~ Lreheabt ch. Kerosene I
JP - homogenized mixture combustion 0 Ha- gaseous fuel diffusion combustion Problems - flame stabilization behind nozzles edge - ignition and flame Spread in all nozzles - increase of number of nozzles and reheat chamber length reduction - separate combusth without mixing of flames - even spread of fuel and air in chamber section Fig. 4.2. Impleineiitllig the buruuJig in reheating and ramjet combustion chamber. 6-5
1- ' Stabilization: pa Ta Pa I- Recirculation zone (RZ) Ks =aW 6Pa Ta Hydrogen: 0.8
0 Reliable stabilization o.6 of burning on edges
0 Broad range of 0.4 burning without flameout 0.2
Selection of required nozzles edge 0 thickness Fig. 4.3. Hydrogen flame stabilization.
200 '00 70 50 30 20 10
C - Fuel consumption L, A, = idem
Nozzles area EA.= idem [.caiK-\ N, - number of nozzles
Fig. 4.4. Hydrogen flame length in an open air flow. Single flames Confluent flames
0 Design curve
No confluence n
2- Confluence
0, %-- 031 0,3 0,s 1 3 5 10 30 50 100 Wf *b
Fig. 4.5. Disposition of nozzles in pylons.
qb I Burner 3 I*b 1.0 --,--.--&- -- y"- 1.o 0.8 0.8 0.6 f-- - 0.6 c Lt=L, ''.{T< Oa4 dbeq 2 A, - 0.4 i' Relative tenth of one 0.2 norrel chamber 0.2 I I I
Fig. 4.6. Experimeiital hivestigatioils of model burners. 0 7-1
CIAM EXPERIMENTAL TURBORAMJETS
by
Dr M.M.TSKHOVREBOV V .I SOLONIN . P.A X ADJARDOUZOV
CIAM (Central Institute of Aviation Motors) 2, Aviamotornaya st., 111250 MOSCOW, RUSSIA
ABSTRACT m - bypass ratio. N - RPM. In this report some results of the TRE operating process P - pressure. including those performed on specially developed Pa - ambient air pressure. demonstrator engines are considered. Essential technical PS - propulsion system. requirements and design features of the experimental TRE PRC - compressor pressure ratio. designed and developed at CIAM for investigatory study q - dynamic pressure. on a test bench are shown. Test facilities for experimental RJ - ramjet engine. investigations of TRE in "connected tube" manner and R - gas specific work. experimental work technology features are presented. SFC - specific fuel consumption. Some results of the CIAM TRE experimental study Sin - intake recovery coefficient. including engine parameters on changing to ramjet T - temperature, K. operation mode variation. working parameters in the T - engine thrust. windmilling (ranijet) mode of operation, and flow Ta - ambient air temperature. characteristics in the afterburner-ramjet combustion t - time. chamber are shown. A TRE based on TF parts TF - turbofan engine. characteristics matching on changing to ranijet operation TFRJ - turboramjet with TF as a core. mode and engine working parameter variations when TFRJs - TRE with TF as a core and going to windmilling in the ramjet mode are discussed. bypass ARCC. TJ - turbojet engine. TRE - turboramjet engine. NOMENCLATURE AND AP BREVIATIONS TRJ - TRE with TJ as a core. TRJs - TRE with TJ as a core and A - area. . bypass RJ combustion chamber. A rd - 7'RJ outer duct area. VX - flow velocity at ARCC inlet. ABE - air-breathing engine. - mass airflow. ARCC - augmented-ramjet Combustion chamber. a - air - to - fuel stoichiometry ratio SUBSCRIPT coefficient. ni a - ambient air. b=- - relative airflow C - compressor. ni+ 1 in RJ duct or TF ouler duct. col - cooling system, air. cor - corrected parameters. C - flow absolute velocity in compressor or ex - exit. turbine flowpath. f - fan. FV - flight vehicle. 6 - gas, turbine inlet. fi 11 -ratio of free stream cross section area toi in - intake. inlet. inlet area. ni - mean value. - adiabatic exponent. n - nozzle. - Mach number. oil - oil. - flow M at the ARCC inlet. r - ARCC. reheat.
Presented at an AGARD Lecture Series on 'Researchand Development of Rambcramjets and Turboramjets in Russia; December 1993 to January 1994. 7-2 rd - Ram duct. INTRODUCTION S - screen. t - total parameters. In the second half of the fifties in connection with the U - tangential component. successful launch of the first Earth satellite there was a V - flight conditions. significant growth of interest in advanced high supersonic W - windmilling mode. and hypersonic flight technology. CIAM began X - ARCC inlet. fundamental studies of combined ABE whose operation 0 - sea level static conditions. process is based essentially on evolution of the concept of * - throat: critical the bypass engine with reference to flight conditions at high supersonic speeds. At the same time the foundation - of their theory was laid.
Studies were conducted on a variety of TRE types based on TJ and TF as a core and useful fields of application . were determined. For the first time, methods of engine configuration selection were determined;which is essential in the case of the TRE for optimum matching of the parameters of gas gas turbine and ramjets. A theory of TRE operation was developed when switching from one operation mode to another, operation in ramjet mode at maximum flight speed under conditions of high kinetic heating, including the modes with turbofan (turbocompressor) windmilling. The extensive theoretical research into turboramjet performances, structures and application possibilities was carried out before experimental work began.
Bench tests of experimental full-scale turboramjets were carried out at CIAM in the 70-s [l]. For experimental work the TRE types on the basis of TJ (with common ARCC) and TF (Fig.1) were chosen.
Investigations into real operation and performance of the TRE were conducted within the programme of development, together with test bench studies of several versions of full-scale TRE demonstrators, created in production TJ's. There were studied Performances on switching for various operation modes, in the ramjet mode. with cooling of the heat stressed components; also studied were the operation conditions of rotor supports etc. For the first time tests were conducted on an experimental TRJ under simulated flight conditions at M = 3 - 4 (Tin = 300 - 600 C) over a long period (several hours) [ll.
All these investigations were intended for development of a cruise mission FV using kerosene. However, the results obtained and the Russian experience of liquid hydrogen use [3]open up possibilities for development of TRE for FV using liquid hydrogen as fuel. *) Dr., Senior Scientist, Head of Division. **) Dr.. Senior Scientist. Head of Department. **+) Dr.. Senior Scientist. Experimental demonstrator engine investigations are part of an advanced power plant design methodology created parameters. in ClAM - the national scientific centre of the engine industry. CIAM carries out scientific supervision of the work conducted in engine and accessories design bureaus and related research institutes [l]. The principle of this methodology, which assumes stage-by-st age development, 7-3
consists of.building a scientific-technical base prior to the During the TRJ heat state problems study the ARCC was a engine development phase. This principle includes equipped with a perforated screen along its length, a development of new construction materials and distribution system of manifold fuel supply and a separate technologies. design. and testing methods. design and air cooling system. This cooling system could provide development of experimental parts, components, gas variable parameters of cooling air (Wcol, Pcol, Tcol) to generators, engines with gas-dynamic, duration and cyclic the cooling duct i.e. into the space between the external testing. Although this initial phase is characterized by case and the perforated screen of ARCC. The bearing relatively low costs with respect to the total programme it supports of TRJ were reconstructed too. To provide is however very important in providing new engine high operability of engine transmission, thermo-protective , quality and efficiency (Fig.2). Application of the screens were located over support cases. Air cooling - methodology favours aircraft and engine development system with variable parameters (Wcol, Pcol, Tcol) and oil matching. system with variable oil feed were used for bearing support cooling. CIAM has the largest rigs to investigate engines and their components in altitude conditions. The unique The main goals of the experimental research were as experimental base of CIAM allows performance of follows: comprehensive tests of engines of practically any thrust and power. - complex study of operation process; - afterburning ramjet combustor operation; As an example of an up-to-date powerful facility for - engine switching modes and operation stability; 0 investigation of supersonic engines there is an altitude - windmilling and ramjet operation mode characteris- tics; facility with a multimode generator of flow disturbances at engine entry. an automated system for monitoring - passage hydraulic characteristics; experimental processes, up-to-date means of data - power output possibilities; processing (Fig.3). which enable the investigation of - ramjet combustor and nozzle cooling at M = 3.5-4.5; non-uniform flow effects on engine performance under - cooling system impact on operating parameters; conditions closest to those in flight at steady-state and - structure heat state and transmission operability. transient regimes. The facility imitates the flight dynamic conditions along the trajectory. simulating non-uniformity Tests on experimental TRE engines were conducted at the of pressure distributions and pulsations at engine entry, facility with an attached air pipeline. with the air which correspond to inlct operation at different angles of parameters behind the inlet reproduced at the engine entry attack and slip [3]. (Fig.6). The maximum entry air temperature realised at the facility corresponds to approximately M = 5. It is Some results of the TRE operating process experimental provided by heat exchanger and flame heaters (with investigations are presented below 14.5. 6.7.101. recovery of normal air composition by oxygen replenishment) and. in case of need. may be increased. Air pressure at entry reaches 900 kPa. The facility has 1. EXPERIMENTAL TURBORAMJETS AND several hundred channels for temperature and pressure SPECIALIZED TEST CELL measuring in static and djnaniic modes. The facility e dimensions allow testing of any current supersonic engine. This cell was also equipped with an injector exhaust Investigations into real operation and performance of the system. The test facility was equipped with separate TRE were conducted at CIAM as part of the programme controlled air cooling and oil systems. of development and test bench studies of several versions of full-scale TRE demonstrators. based on production TJ's The parameters of air flow at the engine inlet were and their parts. The ex pcrinicntal TRJ was developcd on changed in the range: the basis or two TJ and AB engines of the R-11-300 - Ptin = 0.1 ...02 MPa, family (Fig.4) [SI. - Ttin = 270 - 900 K.
To provide a simple solution to the problcnis of maintain The variable parameters of cooling air for ARCC cooling ability, scrvicc, preparation. manufacturing. and were changed in the range: controlling the unified rani-duct with a shut-off device with reniotcd control was uscd. The ram-duct consists of - Ttcol = 400...900 K, front and back manifolds and four connecting sidc tubcs: - PtcolPtr = 0.95 ... 1.40. replaceablc after-turbinc nozzles which provided different - Wcol/Wxin = 0.03 ...0.23. exit turbine areas wcrc used (Fig.5). 7-4
The bearing cooling air pressure was varied in range up that stable operation of the inlet and the engine is ensured to 0.5 MPa and oil pressure on entrance to bearing [41. a supports was changed in the range: The RJ starting mode corresponds to zero air flow through the ram-duct. The curves B and C ( Fig.7 ) correspond to - Poi1 = 0.16 ... 0.22 MPa. this mode for different values of turbine exit duct area (this area for curve C is bigger than one for curve B ). A The air mass flows through ram-duct and ARCC cooling certain engine nozzle throat area at the selected turbine air tube were controlled by changing the shut-off device exit corresponds to this condition. The region above the damper angles. During the experimental investigations curves B and C corresponds to the reverse flow of the gas measuredthe following parameters : part from behind the turbine to the entry of the engine. - Such reverse flow occurs if the engine nozzle throat area - - total and static pressures and total temperature at the An * is less than the value An*o corresponding to the engine inlet; curve Wrd = 0. With lower values of An* unstable - total/static pressures and total temperature at the operation of TRJ and propulsion system is possible. The compressor inlet and outlet, at the turbine outlet: points below the curves Wrd = 0 correspond to the direct - total and static pressures in the ARCC and in ARCC air flow through the ram-duct. Such flow occurs if the cooling duct: engine nozzle throat area An* is bigger than the An*o. If - total and static pressures in the ram duct; An* is increased additionally, there is an increase of air - total and static pressures and temperature at the outlet flow through the ram-duct at practically the same air flow section of both engine parts ducts: through TJ. The experimental points in the joint operation - total temperatures in both the ARCC and the cooling air mode are located nearly the curves B and C. 0 st reams: - wall temperature of perforated screen along its length: At definite values of corrected RPM the curves Wrd = 0 - cooling air flows through both the cooling duct and the cross the curve A and approach the compressor stall bearing supports: curve. But usually for TRJ the corrected RPM for the - temperatures of bearing support cases: transition to RJ-mode is located right of the curve - cooling air temperatures at entrance and exit of supports: crossing point. Therefore, an error of control at the start - oil temperatures at entrance and exit of supports: of RJ (the shut-off device prematurely open) presents a - oil flow through bearing supports: danger rather for stable operation of the inlet (Fig.8). As - engine rotational speed and some other parameters. follows from the Figure. the forestalling shut-off device opening leads to lowering of the TJ inlet corrected airflow The layout of the engine and ARCC preparation is shown with corresponding consequences. This was confirmed by in Fig.6. Engine continuous operation time at M =4 ...4.5 tests of TRJ [5]. was about 3 hours. The rational control of an engine on transition to RJ-mode process must provide not only stable operation of the 2. THE EXPERIMENTAL INVESTIGATIONS OF TRJ propulsion system but also thrust onset without significant TRANSITION TO RAMJET MODE drops.
Investigations have revealed a rather narrow operation The transition from gasturbine to ramjet mode (also interval, inside which a stable working process of the windmilling mode) is of great interest for TRJ. The first power plant with parallel operation of TJ and RJ parts part of tests were directed just to these modes [6].Some could be realized. Maintaining the TJ operating parameters results of these investigations will be discussed. The with respect to the limitations provides propulsion system rational value of flight speed for transition to ramjet mode reliability on transition to RJ mode. In transition operation depends on the design operating parameters, the TJ the profile of ARCC inlet flow velocity changes control mode and the size correlation between typical considerably. The transition to RJ-mode finisheds with TJ propulsion system duct areas. The transition Mach number windmilling, which may be used to drive engine is usually in the range 2.5 - 3.5 [8]. accessories.
The working lines on the compressor map for different In the transient mode the structure of flow at the ARCC modes of operation are shown on Fig.7. The gasturbine entry is strongly changed (Fig.9). The results of TRJ tests mode Corresponds to the curve A. The joint operation of show that ARCC inlet flow velocity increases at constant TJ and RJ is realized in the transition mode. In thcse rotation speed approximately in proportion to the increase conditions the pressure behind the turbine is grcater than in engine nozzle area. In this case a considerable that at the ram-duct entrance. If the shut-off device opens. distortion of flow takes place. For example, in the the control system of TRJ must respond in such a way TJ-mode the distortion of flow is not more than 10 - 20% (the curve A in Fig.9), but in the engine parts joint mode 7-5
the maximum flow speed can exceed the mean flow were: - ARCC combustion performance with the velocity by about 2 - 2.4 times (the curves B and C in distortion of ARCC inlet total pressure and flow Fig.9). Therefore, special measures possibly would be velocities fields; demanded for stable combustion in ARCC in the -joint operation of ARCC and its cooling system; transition mode. - heat state of ARCC structure; - heat state of bearing supports in the ramjet operation mode. Some results of experimental investigations on full-scale ARCC in the 3. THE RJ-MODE OF TRJ WITH WINDMILLING TJ ramjet mode working process and estimation of cooling system operation influence on engine performance are discussed below. .. In the RJ-mode TRJ can operate with either open or closed gas turbine part of the engine. The curve D on 4.1. Operation of ARCC cooling system in RJ modes Fig.7 corresponds to the windmilling mode of the TJ . turbocompressor. The pressure losses in TJ in windmilling The operation of TRJ is characterized by a wide range of mode depend very much on compressor corrected air flow ARCC inlet parameter change with significant distortion (the curve A in Fig.10 with the ram-duct closed). In the of inlet total pressure fields. The value of inlet total RJ-mode the TRJ pressure losses do not exceed 5 - 7% if pressure distortion depends on the operating mode and the TJ is closed (the curve C in Fig.10) or operated in the value of the corrected air flow through the engine. In the windmilling mode (the curve B in Fig.10). Application of course of the experimental investigations of TRJ the the TJ windmilling mode provides an increase of engine ARCC inlet parameters of flow were characterized by the a air flow as the additional part of the air flow is passed values: through TJ and gives some possibilities for accessories drive. However these possibilities depend on TJ Ttx = 450 ...900 K, efficiency in the windmilling mede. On Fig.11 the Mx =0.1 ... 0.25. experimental curve for TJ efficiency in windmilling with energy output is presented [4]. The radial distortion of ARCC inlet total pressure was 6...12% with the gasturbine core wind- milling and In the RJ-mode the ARCC inlet flow has relatively little 7...14% with the closed gasturbine core. Nozzle flap speed distortion. At the begining of transition mode the location and the cooling air blowing through perforated region of maximum speed is located in the central part of screen slightly influence the shape of the total pressure the ARCC. and very strong distortion takes place, but in fields in the nozzle throat section. Inlet total pressure the windmilling mode it moves to the peripheral part of distortion decreases along the ARCC and corresponding ARCC (the curve D in Fig.9). The unevenness of ARCC to estimates. in the throat nozzle section it is not more inlet flow decreases if the corrected air flow also 2...3% when combustion takes place. ARCC combustion decreases. efficiency under the conditions considered depends mainly on the air-to-fuel ratio coefficient and is in a range typical for afterburners. 4. COOLING OF THE AFTERBURNER-RAM COMBUSTION CHAMBER The length of the combustion zone significantly influences a distribution of cooling air along the perforated screen of ARCC. Combustion zone length can be characterized The ARCC and nozzle are one of the most heat loaded indirectly by the static pressure drop on the ARCC screen units of the TRE operating in the ramjet mode. To measured along the ARCC. As experiments show (Fig.13). provide reliability of these units. it is natural to use an air fuel combustion is finished in a distance of less than 2/3 . cooling system in which a part of the air is taken away of ARCC length. from the engine flowpath to the ARCC cooling duct (i.e. a space between the ARCC case and the perforated The cooling air flow variation and its distribution along screen). the ARCC are defined by very complex processes in non- uniform flow with combustion in ARCC. The hydraulic In the experimental TRJ. an air cooling system with a characteristics of the ARCC cooling duct depends on the perforated screen was used for cooling the ARCC and the discharge coefficient of the holes and duct pressure loss slot system was used for cooling the nozzle (Fig.12). coefficient. These factors variation was examined in the course of the investigation [7]. Experimental investigations of heat state and operability of high heat loaded TRJ units were carried out in the Analysis of the cooling air parameters in the ARCC ramjet mode operation with windmilling of gasturbine cooling duct has revealed that the cooling air temperature core and simulation of Mach numbers in the range of along the cooling duct changes slightly. When the total 2.5 ...4. The next heat state problem tasks investigated pressure ratio is Ptcol/Ptx < 1 the cooling air temperature 7-6 at the duct exit is 20 ...30% higher than the value of the 4.2. Heat state of ARCC duct inlet temperature (Fig.14). This temperature variation is explained by penetration of hot gas into the cooling The heat state of ARCC, which is characterized by the duct through the holes in the perforated screen part of the perforated screen temperature Ts, depends on the engine duct. However, the change of cooling air temperature in operating mode, flight conditions, parameters and flow of this case is not great. It should be noted that a small cooling air. Using the hydraulic characteristic of the air change of Ttcol along the ARCC has no practical cooling system, one can consider the correlation between influence on the hydraulic characteristics of the ARCC relative cooling air flow Wcol/Wx, mean screen cooling duct. and it is possible to use an assumption temperature Tsm and flight M number bearing in mind, about the constant value of the cooling air temperature for example, the non controlled cooling system of ARCC along the ARCC cooling duct. without conditioning of cooling air (PtcoLin = Ptx, . TtcolJn = Ttx). The two possible laws of PS control in RJ The results of experimental investigations showed that operating mode are characteristic: distribution of cooling air along the ARCC perforated screen depends on the variation of both the cooling air - with a constant value of engine inlet corrected air flow, and the ARCC gas parameters in a wide range. The which corresponds to the constant value of ARCC inlet variations of the corrected airflow in the underscreen duct Mach numbers Mx independently of flight conditions: and the gasflow in ARCC as a function of the ratio of total and static pressures in the characteristic sections - with a decreasing value of inlet engine corrected air were obtained. flow .vs flight M number, which corresponds to a drop in ARCC inlet Mach numbers: in this case the air flow The results of experimental investigations enabled through the engine is limited by maximum inlet definition of the hydraulic characteristics of the cooling productivity. duct and ARCC in joint operation. In this case the following reasons were adopted as a base : The variation of mean screen temperature and relative cooling air flow is shown in Fig.16 depending on flight - For pressurec defining characteristi cross section which Mach numbers in ramjet mode operation and the is located at approximately half length of ARCC was air-to-fuel ratio coefficient a = 1.2. adopted: in this section the combustion is mainly completed and main part of pressure losses is realized: - It follows that the relative cooling air flow changes very air total pressure in cooling duct characteristic cross little, although the change of TRJ parameters for these section is taken as equal to the mean value between the two laws of PS control is significant. At flight Mach values for inlet and exit of cooling duct: ARCC total number M = 4 the mean screen temperature is Tsm = pressure is taken as equal to the ARCC exit total pressure. 1150... 1230 K and the relative cooling air flow is 12% approximately. - one part of the cooling air from the cooling duct blows out through the perforated screen and the other part of the At flight Mach numbers M < 4 the cooling air flow in an cooling air is used for nozzle cooling. Obtained by means uncontrolled cooling system is significantly higher than of experiment and calculation. the hydraulic characteristic necessary for optimum cooling. Because of this the mean of ARCC and its cooling system is shown in Fig.15 in the screen temperature at lower flight Mach numbers is rather form of the dependence of the ratio of the corrected low. For the two laws of PS control there is very little airflow in the cooling duct to the gasflow in the ARCC, distinction between the mean screen temperatures. For on the corresponding total pressures ratio. The curves are example, at M = 4 the increase of Mx from 0.13 to 0.18 given at constant values of corrected flow density in the provides a decrease of Tsm by 80 K. In this case, the duct and ARCC. The hydraulic characteristic obtained cooling airflow decreases by 2% approx. helps to form the principles governing coolant flow. TRJ throttling at flight Mach number M = 4 by decrease The hydraulic Characteristic of the cooling system shows. of fuel supply to ARCC causes a decrease in Wcol that when the total pressure ratio is near to 1. the mass because of lower heat pressure losses. However, the flow ratio change is rather small in a wide range of both decrease of air cooling flow is lower than necessary for ARCC and cooling duct parameter variation. an optimum cooling system. In this case, bigger coolant flow also provides low screen temperature. These data enable an approach to the development of an optimally controlled cooling system for ARCC. For the The results of full-scale TRJ experimental investigations development of a controlled air cooling system with the in simulated flight conditions at Mach numbers M = possibility of cooling air flow change in different engine 2.5 ...4 have shown that with an uncontrolled ARCC air operating modes it is necessary to select the total pressure cooling system with perforated screen in non-design ratio somewhat higher than 1. flight conditions (lower flight Mach numbers or engine 7-7
throttling at M = 4) the relative cooling air flow is maximum flight speed, fuel cooling capacity significantly higher than the required value. This leads to (hydrocarbons, LH2) etc. thrust performance aggravation in these modes. Improving per- formance is possible by choosing rational cooling air parameters and using the controlled c(ooling air flow 5. BEARING SUPPORTS HEAT STATE IN RAMJET system. MODE OPERATION (WINDMILLING CORE)
The results of TRJ experimental investigations enabled development of a method to compute the engine thrust The bearings which are used in the supports of the TRJ performance characteristics taking account of the influence gas turbine core define the engine's reliability to a high - of cooling air flow mixing with the main stream. degree. Reliable and long-term operation of bearings is Performance estimatioas were made according to the two possible only when the bearing temperature is lower by types of ARCC which are used in different types of TRE 40 ..SO C than the tempering temperature of the bearing's . [9]: ARCC in engine common exit duct ("cylindrical" material [7]. ARCC) and ARCC in separate ramduct ("coaxial" ARCC). In performance definition account was taken of the Synthetic oils are used for cooling and lubrication of nonuniformity of total pres- sure and temperature in aviation engine bearings. Increased flight speed leads to ARCC and nozzle and the screen heat state. the increase of inlet engine air temperature and heat loads on units and construction elements of the engine. The heat It was considered that the cooling system was operating flow into the engine's supports increases if an oil with 0 in optimum mode, in which the necessary value of limited operating temperature is used. When flight cooling airflow is defined by the limit value of screen speeds are higher than M = 2.5 ...3, operability of temperature. The ratio of thrust and SFC values to their gasturbine engine supports should be achieved using conditional quantities, without accounting for cooling air special structural measures. flow output (fig.17) was used as a measure of the degree of influence of cooling air flow output. The performances One of the most effective methods to ensure bearing given are dependent on flight Mach number along the support operability is the use of active thermo-protection. trajectory with dynamic pressure q = SO00 kg/ni2 and at This comprises the setting up of a protective screen over air-to-fuel ratio coefficient a = 1.2 for screen materials the support case under which conditioned air from the with different temperature limits. engine cooling system is supplied. The layout of an experimental TRJ support is shown in Fig.19. All supports A significant increase of cooling air flow takes place at for experimental TRJ had autonomously controlled oil flight Mach numbers M = 4...S. At the higher flight Mach supply and thermo-protective screen under which cooling numbers the necessary value of the cooling air flow air with controlled parameters was supplied. becomes very high with significant increase of loss of performance (Fig.17). Some results of heat state investigation of experimental TRJ bearing supports in ramjet operation mode with It is possible to decrease the amount of cooling air at high windmilling of gas-turbine core are discussed below. flight Mach numbers by the use of a conditioning system, 0 which decreases the cooling air temperature. Fig.18 as an The summary value of heat flow into the oil cavity of example shows the influence of cooling air temperature supports depends on: Ttcol on both thrust and SFC for two different types of - heat of friction in bearings and seals; ARCC using the same screen material in the flight - heat flow through the details of the bearing supports; conditions M = 5. H = 25 kni. TRE with common ARCC - heat flow from hot air, which penetrates from the engine has a higher thrust value of 7...18% and lower SFC of 2% flowpath through the seals into the oil cavity. than TRE with separate ramduct because of the latten larger area of cooling surface. To provide equal values of The experimental investigation of TRJ support heat state cooling air flow output in both engine types the cooling showed that along with the TRJ transition to ramjet mode air temperature in TRE with "coaxial" ARCC must be with windmilling of gasturbine core the redistribution of lower by 100...150 K than in another type of ARCC. The relative values of each of the above mentioned heat SFC variation vs Ttcol is rather small. components takes place. It has been determined that part of the friction heat in the bearings decreases because of To ensure ARCC operability up to flight Mach numbers the decrease of both the rotor RPM and the axial forces. M = 4...4.S there is no need for a conditioning system. At Heat flow through the details of supports and heat coming higher flight Mach numbers the use of such a system is from hot air penetrating into the oil cavity components a necessity. A rational type of cooling air conditioning rises because of inlet temperature increase and aggravation system depends on duration of cruise flight at the of seal operating conditions. The total heat flow into the e supports rises as compared with gas-turbine mode 7-a
(Figzo), so the exit support oil temperature rises too and the maximum design value. The fan corrected RPM achieves its limit at certain flight Mach numbers (Fig.21). decreases by the same order. Let us consider the peculiarities of fan operation at low Autonomous oil supply to each engine support in rotational speed (Ncor/Ncor o 0, fan expended work is positive (Ptcol/Ptin) > 1.7 the cooling air flow ceases to rise and and the fan pressure ratio value is near 1. When the the support exit oil temperature also stabilises. counterpressure behind the fan reduces, the operating point is shifted to the right on the characteristic, and this A cooling air supply under screen space enables a is followed by "stretching" of the velocity triangles, significant decrease in the support exit oil temperature and decreased of angle of attack, Delta Cu, ex- pended work ensures engine transmission operability at flight Mach and PRCf. The point 2 corresponds to the case when the numbers M > 2.5 ... 3. The use of a thermo-protective Delta Cu and expended work values are equal to 0. i.e. screen, under which cooling air is supplied, ensures the fan operates in the free wind- milling mode and engine support operability up to M = 4...4.5 (Fig.21). presents hydraulic drag in the air flow (PRCf < 1). The combination of such points at various rotational speeds gives a line of fan windmilling mode operation (curve B). 6. TFRJ PERFORMANCE IN THE RJ OPERATING With further shift of the operation point down along the MODE characteristic (up to exit section choking) the transition to "turbine" operation mode takes place. At point 3 the Delta Cu and expended work values are negative, and the The operating parameters. thrust and SFC performances fan can provide some power. and control principles of TRE types based on TF core (TFRJ and TFRJs ) at higher flight M numbers (transition The example of zero work (free windmilling) performance and ramjet operation modes). including some experiments lines for fans with different design pressure ratio values with imitation of these conditions. are considered below. is shown in Fig.35. The less PRCfo, the less pressure loss and the wider the range of air flow in the fan on windmilling. At design pressure ratio values 2...2.5 and 6.1. Fan performance in the RJ operating mode (Wcor/Wcor 0) < 0.4 the total pressure decrease in windmilling fan flowpath is of the same order as that of To maintain low losses in the air intake with flight Mach the pressure losses in the separate ramjet duct of the TRE number increase the fan corrected airflow should be PI. decreased, as follows from the equality of airflows in the stream at infinity before the air intake and in the fan at In the RJ operating mode the airflow in the engine inner the area ratio fin . Ain / Af part (gasturbine core flowpath) is near to zero. Under these conditions the fan turbine becomes the consumer of power being produced by the windmil- ling fan. The turbine resistance to rotation is defined by so called ventilation losses, which arise as a result of energy feed to the airflow in the turbine flowpath and also by the turbine disc friction losses. The turbine driving work is where (Wfcor/Wfcor 0)is the relative corrected airflow in proportional to the air density and to the rotational speed the fan. This tendency is particularly clear when the area value in the powcr of about 3.7. A turbofan working in ratio becomes constant (Fig.33). typical RJ operating mode condition analysis has showed that because of low rotational speed ad- ditional pressure For example. when area ratio is equal to Ain/Af = 1.5. up decrease in the fan driving the turbine is usually small 10 certain flight M numbers. the airflow in PS is defined [91. by the fan working at the maximum RPM (fin < 1. the part A-B of the operating line). At higher flight M Experimental investigations of TFRJ models based on numbers the airflow is limited by air intake maximum production TFs with replaceable nozzles (Fig.26) and productivity (fin = 1. the B-C part).The typical area ratio external oil supply unit were performed on the special would be approximately equal to 1...2 [9]. As follows CIAM test facility on the connected tube scheme (see from Fig.23, at M = 4...5 , corresponding to such area above). ratios corrected airflow decreases to nearly 20 ...40% of a 7-9
The higher the flight Mach number and inlet pressure An example of engine parameter and performance e recovery and the lower the design fan pressure ratio, the variations in transition to RJ operating mode is given less the pressure loss in the fan on windmilling. firstly for TFRJs (solid lines on Fig.28 ... 30). As a Experiments on full-scale TFRJ models have demonstrated characteristic value the area ratio Ain/Af = 1.3 was the possibility of RJ mode performance with rather low adopted in this case. The points B and E correspond to losses with a windmilling fan (Fig.27). It is possible to the start and finish of the transition to RJ mode process; use a windmilling fan for driving engine accessories in the the point W corresponds to the TFRJ parameters just at RJ operating mode. The expediency of using the fan to the moment of going to the windmilling mode. drive the power plant accessories in the ramjet mode of operation depends on its efficiency as a "turbine". values At M = 2.7 the inlet maximum airflow capacity is of the work taken-off, thermal conditions and other reached. With further flight speed increase the engine-inlet factors. matching for maintaining finmax is provided by a corresponding decrease of the combustion chamber fuel If a windmilling fan is loaded by external energy supply. For the first time, the inner and outer nozzle consumers the operating points on the performance map throat areas An*l. and An*2 are maintained constant. As lie beneath the zero work line. However in the TFRJ a result at M > Mb the temperature ratio Ttntin and the ramjet operating mode this deflection is usually turbo-fan RPM quickly decrease. Maintaining the insignificant and it is possible to adopt the zero work line maximum ARCC flow capacity in conjunction with as an operating line for the engine RJ operating mode. decreasing the airflow in the engine gas-turbine part leads to a rapid increase in the bypass ratio and the b value, i.e. 0 The possibility of using a windmilling fan to drive the the engine operating mode approaches the RJ mode accessories in the RJ operating mode is a significant (Fig.29). As a result the overall heat supply in the engine advantage of TFRJ and TFRJs. increases because of maintaining the ARCC rating at minimum air-to-fuel ratio. However, reduction of the In some cases, particularly when the fan design pressure engine nozzle expanding ratios is found to be a prevailing ratio is rather high it is possible for TFRJ thrust factor in this case, and as a result the engine specific performance in the RJ operating mode to be improved by thrust slightly decreases. For this reason the thrust arranging the bypass duct around the fan case with a increase rate is somewhat lower when M > Mb; the SFC corresponding controlled shut-off device. increases simultaneously (the section B-W on Fig.30).
After going to windmilling mode (point W at M = 3) 6.2. TFRJ transition to RJ mode engine governing is provided by optimum variation of the nozzle throat areas. Under the flight conditions Engine transition to RJ operating mode essentially consists corresponding to Mw the b value is equal approximately of switching off the gasturbine part without engine thrust to 0.9. i.e. about 10% of the overall airflow goes through loss. For the TFRJ and TFRJs a gradual transition process the engine inner part (Fig.39). When the flight M number is typical at a certain range of flight Mach numbers. The is higher than Mw it is expedient to decrease the airflow transition process begins at Mach number Mb when in the engine gasturbine part and the latter's drag by engine airflow at maximum r:ting approaches the inlet decreasing the An*l. For this purpose it is possible , for maximum airflow capacity (point B on Fig.33) and example. to shift the inner nozzle central body. The finishes at flight M = Me when the whole (or practically nozzle area decreases, which corresponds to inlet engine whole) airflow is feeding through the outer (ranijet) matching with minimum pressure losses and the ARCC engine part (b = 1). As flight speed increases in the range inlet flow velocity limitations; the outer nozzle throat area of Mb < M < Me it is necessary to gradually decrease the being constant as is shown in Fig.29 (section W-E). fan corrected airflow by decreasing the core combustion chamber fuel supply and the throat areas of air-inlet and Along with An*l decrease the turbine expansion ratio and nozzle thus maintaining minimum intake pressure losses. work fall off. which leads to a decrease in fan RPM and In the transition proceeding the operating process in the pressure ratio (Fig.28). When the flow in the fan turbine ARCC is carried out at maximum inlet flow velocity Vx becomes sub-critical, the inner nozzle area decrease leads and minimum air-to-fuel ratio value. to a decrease in compressor turbine work along with core RPM and corrected airflow lowering. As a result the To lower the TF airflcw it is necessary to decrease the engine airflow redistributes in favour of RJ part, and with engine RPM. It is expedient to use a rather simple niethod full gas-turbine part shut-off (An*l = 0) b = 1. At this of RPM reduction by gradually decreasing the turbine moment (M = 3.4 on Fig.29.30) the transition to RJ inlet temperature until the engine goes into windmilling operation mode is being completed. The corrcs- ponding mode. In this case the nozzle throat areas or the common operating point on the fan performance map is found nozzle throat area must decrease just after cessation of the practically on the zero work line (point E on Fig.28). a combustion chamber fuel supply. 7-10
Along with inner nozzle area lowering the engine specific TFRJ thrust performance variation vs ight M number is thrust increases because of the decrease in core losses nearly the same as in the case of TFRJs. Any difference 0 and increase of the overall heat supply. For this reason in is connected with fan efficiency decrease and higher fan the range Mw-Me engine.thrust rises, and SFC varies duct pressure losses. In the RJ operation mode (M > Me) slightly. At M = Me flight conditions the TFRJs inlevengine matching is provided by a corresponding parameters differ from the "pure" RJ parameters by the nozzle area decrease, and this is accompanied by RPM value corresponding to the additional pressure losses in and engine flowpath pressure losses decreasing (dotted the windmilling fan. lines on Fig.29.30).
In the RJ operation mode (M > Me) inleuengine matching When analysing engine parameter variation on transition is provided by a corresponding decrease of the RJ part to RJ mode it was assumed that the air inlet operates at nozzle throat area An*2. The fan RPM and pressure losses maximum capacity rating. This approach, although it does are gradually decreasing (PRCf -> I), and the engine not account for possible inlet capacity limitations when M thrust performance is approaching RJ performance. < Mmax, still enables simplification of the analysis and reveals the main correlations. The taking into account of Now let us consider some peculiarities of the TFRJ real variable geometry inlet performance characteristics operating process when transition to RJ mode is carried would not change engine parameter variation principle, out. using as an example the engine with the same cycle but in the event of lower inlet capacity the transition parameters and transonic (M = 1.3) to hypersonic (M = flight M numbers range would shift to lower flight 4...5) thrust ratio as discussed above. For this condition speeds. the area ratio Ain/Af would be somewhat higher than for the TFRJ engine type (in correlation of transonic specific The thrust performances of TRE in the RJ operation mode thrust values). In this case the area ratio for TFRJ would differ from that of the "pure" RJ and are conditioned be about 1.4. as against 1.3 for TFRJs. mainly by additional pressure losses in the engine flowpath. The TRE performance aggravation in The TFRJ operating parameter variation peculiarities are comparison with "pure" RJ at hypersonic flight speeds is conditioned mainly by the mutual influence of gas and air possible also in connection with higher airflow or other flows in the common ARCC inlet. As a result of core agent expenditure for engine structure cooling. pressure ratio decrease vs flight speed. which is not compensated by fan duct outlet pressure variation. fan Some engine thrust performance aggravation is also turbine specific work decreases. Because of the fan duct possible in connection with the part of airflow energy flow influence the fan operating line on the Performance used for driving auxiliaries. In this case the use of a map deflects down from the corresponding line for TFRJs. windmilling turbofan would be possible. When power and the bypass ratio vs flight M number rises more take-off is low. this type of auxiliaries drive would be rapidly (Fig.28.29). Meanwhile corrected airflow variation expedient in spite of rather low fan efficiency when vs flight M number is approximately the same as for the running in the turbine mode (Fig31). TFRJs. CONCLUSION As a result of t.he higher Ain/Af ratio the TFRJ propulsion plant approaches the maximum inlet capacity conditions The set of problems posed in connection with the e at a higher flight M number than that of TFRJ (Fig.39). application of the new types of ABE, concern the need to and corresponding fan pressure ratio is nearer to 1. Then. increase flight speed to M = 4...6. since the transition beginning the TFRJ operation mode is nearer to the RJ mode. Lowering of the turbine inlet Broad programmes of investigation of engines for high temperature at M > Mb leads to b value increasing to speed flights were carried out in Russia from the 60s to almost 1 when TFRJ approaches the windmilling mode. the 80s. The corresponding fan operating point lies practically on In the framework of these programmes an investigation of the fan zero work line (Fig.28). Then. in the case of TFRJ advisable engine concepts and operating process the completion of transition to RJ mode practically parameters of TRE were conducted. Test programmes on coincides with cessation of the fuel supply to the experimental full-scale TRE of different types were combustion chamber. As a result. the Mc number for carried out at CIAM. The experimental TRE were TFRJ is lower than that of the TFRJ (in the example assembled from units of production gasturbine engines. considered respectively Me = 3.1 against 3.4. Fig.39). Flight conditions corresponding to Mach numbers M = 4.0 ...4.5. were simulated at the facility. A rather fast TFRJ transition to RJ mode is a consequence of the simultaneous effect of the turbine inlet temperature The results of theoretic and experimental investigations of decrease and the gas dynamic influence of fan duct flow TRE allow us to define: on gasturbine part performance. 7-11
- rational methods of transition from gasturbine to ramjet mode: - the conditions of stable operation of the propulsion systems with TRE in the transition mode; - performance characteristics and pressure losses in the windmilling mode for TRJ and TFRJ turbo-compressors; - the expediency of the application of windmilling mode in the RJ operation mode; - the performance of ARCC with distortion of ARCC inlet total pressure fields; .. - the conditions of joint operation of ARCC with an air cooling system and the hydraulic characteristics of the cooled ARCC: - - the heat state of ARCC with perforated screen and of the bearing supports; - the influence of ARCC air cooling on TRE thrust performance: - the possibilities of heat protection of bearing supports by controlled oil and cooling air supplies.
0 Solving the complicated problems posed by the development of new aviation techniques requires high expenditure. International cooperation in the building up of a scientific-technical base for adianced powerplants is one of the important factors in the successful realisation of projects for 21st century air transport.
ACKNOWLEDGEMENTS
The authors wish to acknowledge contributions from collegues during the preparation of this paper and thank the Directorate of Central Institute of Aviation Motors (CIAM) for permission to publish. 7-12
REFERENCES
2. V. A. Sosorinov, "Some Aspect.s of Hydrogen and Ot.her A1 t.e!rnc?t.ive Firel s for Ap- pl i cat.i on in Air-Rreat.hi ng Engines;" TX TSARE, Athens, Greece, Sept.emher, 1989.
5. V. A. Sosoirnov, "St.iidy of Pro- p1i1si on for High Vel oci 1;y F1 -i ght." X TSAHE, Not,t.i ng- ham, lIK, Sept-ember, 1991.
6. V. A. Sosoiinov, M .M. Tskhovre- hov, V.T.Solnnin, V.A.Pa1- ki n "The St,iidy of Experi - ment,al Trlrboramjet.s". ATAA PAPRR 92-3720, 28t.h .Tni nt. Prop11 1 sion Cnnference, .Jill y, 1992.
7. V. A. Sosoiinov, V. T .Sol nni ri, M . M .Tskhovr~h~v, P. A. Kad- jardoiizov, V. A. Pal ki n, "The s t.I I d y of Experi mpn t.al Tiirhornmjc!t.s: Heat. St.atx and Cool ing Prohl ems". ATAA PAPFR 93-1 989, 29t.h Joi nt, Prop111si on Conference, .Tirnc?, 28-30 1993, Mont.erey, CA.
8. M.M.Tskhovrebov, V.A.Pa1- kin, "Combi ned Engines for Hype rson i c F1 i ght" . TCAS- -92-3.4.3. 18-t.h TCAS Cong- ress, 21 -25 Sept,emher, 1992, Rei jing, PRC. 7-13
Without energy With energy t rans f e r r ran8 f er to ramjet part to ramjet part
I.
Fig 1 TRE types [8]
-Forecasts -I' Le sso ns 1e a r n ed" Engine design & -Propulsion requirements EEG$> developement - R&D programmes -aircmff-propulsirm process sysfem aplimisafaun
CIAM engine DB's Ts AG I systems, gasgenerators VIAM / -advanced engines ezperi- mental components ness. life time. erplda- -dual components relia- lion technofogy. ecology) IEPT -test facilities serial production -environmental probfems plants - AS institutes
new materials and
-erperimental core engines and demonstrat -aircraft-propulsim in1 e
-new generation engine
Fig 2 Scientific base for AE development 181 7-14
1 RECEIVER 2 ALTITUDE TEST CHAMBER COLD AIR PIPE 4. HOT AIR PIPE 5. INLET 6 HEAT EXCHANCER llli :I Z CONTROL PANEL FI umsumc SYSTFM
state and transient conditions-[3]
0 Complex Study of operation al 0 Engines switching modes 0 Afterburning and ramjet cornbush !ARC) ooeration 0 Ramjet operation mode characteristics (M > 2.5) 0 Windmilling mode 0 Power output 0 Passage hydraulic characteristics 0 Structure heat state 0 ARC and nozzle cooling at M - 3.5 - 4.5 0 Transmission operability Cooling system impact on T and SFc
Elg 4 ClAM experimental TRJ [1,4,5,10] 7-15
Gas turbine engine mode
fianilet cngine modc Fig 5 Experimental TRJ layout [5]
Controled sliut
~ Fig 6 Test facility schematic layout [6] PRC
Compressor corrected air flow 7-16
1.o 0.5 0.
Flow velocity
Fig 9 Row velocity at ARCC entry distortion [4.6]
Fig 8 TRJ performance at transition mode [5] 1 - stable operation 2,3 - unstable operation
0.7 I I ‘1 7 I I I II LTJl windmilling 0.6 .
0.5 corrected air flow Fig 11 Windmilling TJ with ower offtake relative efficiency [4P Fig 10 TRJ pressure losses at RJ mode (windmilling core) [6] 7-17
Cooling air -----
S: 3 scotions Fig 12 ARCC schematic layout 171
Fig 13 0.06 Static pressure variation along ARCC [7] 0.92 I-Tex /Tx=2; Mx=0,16 2-Tex /Tx=% ~p0.2 3-Tex /Tx=2,6; Mp0.22 0.88 , I I I
1.o * I,
0.8 0 0.2 0.4 0.6 0.8 'CF
Fig 14 Cooling air temperature variation along cooling duct 171 ucorcol Wcor 0,20
0.15
0.10
0.05 e I J 0 L I I I. 1.o Ptcol /Pter Fig 15 ARCC hydraulic characteristics cooling duct [7]
TsrnBK e 1000
800
600 2.5 3.0 3.5 N
Fig 16 Relative cooling airflow and mean screen temperature va flight Mach number [7]
1 - Wc0rx= const 2 - Wcorx= V= 7-19
4 nn Relative SFC
1 p2- -4/
I-
Relati I -- 0.8
0.6 sir, 4.0 4.5 M
a Fig 17 Relative thrust and SFC variation v11 flight M number [7] 1-common " cylindrical' ARCC %separate "coaxid ARCC
Relative thrust .
0.8
0.6 600 800 1000 Tt col
Fig 18 Relative thrust and SFC variation vs cooling air temperature [7] 1-common "cylindricalD ARCC %separate "coaxid ARCC 7-20 ,
Fig 19 Bearing unit shematic layout [7]
I 1 11 II
~ Heat Flows: wary tmm~- Friction 20 Bearing support Fig 1.0 0Others heat flows [7]
0.5
0
Fig 21 Influence of flight Mach number on support exit
I oil temperature [7]
2.5 3.5
Fig 22 &dative support exit oil temperature and cooling air flow in dependence of air pressure ratio [7] 7-21
c Wf c
0.8 Fig 23 Relative fan corrected airflow in dependence on free stream to fan 0.6 areas ratio and flight M number [SI 0.4
0.2
0
Fig 24 Fan stage parameters variation along the fan characteristic (low RPM region) 191
I ACU (0 N cor =const
0 qfcor
pig 25 Zero work lines for fans with different design pressure ratios [SJ 1-22
@ TFRJ @ TFRJs
Replaceable nozzle<
Fig 26 "FRJ models shematic layout 161
Fig 27 Fan pressure ratio at windmilling mode [lo] .+PRCfo= 1.8 ...2.6
PRCf
Fig 28 TFRJ fan characteristics map [6] 7-23
T
Fig 29 TFRJ parameters variation when going to RJ mode [Q] -TFRJs --- TFRJ b
Fig 30 TFRJ relative thrust and SPC variation vs flight M number [Q] -TFRJs --- TFRJ 7
1 2 3 4 M
I thrust
Fig 31 Power offtake influence on engine SFC at RJ operation mode [IO]
8-1
SCRAMJET CFD METHODS AND ANALYSIS a PART 2. SCRAMJET CFD ANALYSIS
NUMERICAL SIMULATION OF SUPERSONIC MIXING AND COMBUSTION APPLIED TO SCRAMJET COMBUSTOR
by V.KOPCHENOV, K.LOMKOV, S.ZAITSEV, 1.BORISOV
CIAM (Central Institute of Aviation Motors) Aviamotornaya street 2 111250 MOSCOW RUSSIA
2.1 Introduction 2.2 Mathematical model and numerical method 2.3 The main results of the CIAM investigations on supersonic mixing and combustion enhancement 2.4 Some notices about the jetlshock interaction and mixing enhancement 2.5 Evaluation of combustion efficiency and choice of acceptable design parameters for combustor 2.6 Estimation of the role of chemical kinetics
2.1 Introduction numerical simulation of supersonic combustion. The m&ed k-E turbulence model of Launder-Jones was It is necessary to note, that the problem of used. The equilibrium and nonequilibrium combustion supersonic combustion of hydrogen in air stream was models were realized. The influence of the investigated in Russia both experimentally and concentrations fluctuations on the reactions rates was numerically. It is not possible to give the full survey of also taken into account with the aid of the model [9]. publications about investigations concerning with numerical simulation in short description. As an The numerical investigation of supersonic example, monographs [1,2] may be pointed out. The combustion of hydrogen jet, injected tangentially along boundary layer model was used in [l]. The experience the wall into supersonic air stream, was performed in numerical simulation of supersonic combustion within the scope of two - dimensional parabolized within the scope of Reynolds averaged Navier-Stokes Navier-Stokes equations [13-151 . The system of 13 equations and their parabolized version is presented in reactions for 9 components was considered. The PI. algebraic turbulence models were used (modified Prandtl and Baldwin-Lomax). The influence of finite rates of chemical reactions and concentrations fluctuations on combustion Detailed survey of modern numerical investi- efficiency in duct was investigated in [3-51 with the aid gations of the supersonic combustion on the base of of boundary layer model. By this, the msion flame three-dimensional Reynolds averaged Navier-Stokes sheet model was used in [3], and suBciently detailed equations or their parabolized version is presented in kinetical scheme, including 20 reactions for 9 [16-191. It is necessary to note, that simplified kinetical components was considered in [5]. The influence of scheme [20] was most commonly used. The possibility duct geometry factors on combustion efficiency was also to include more detailed kinetical mechanism into this analyzed. class of mathematical models is directly dependent upon the progress in computer power. The supersonic combustion process near the wall was numerically simulated within the boundary layer The short description of the mathematical approach [6-81. The necessity to take into account the model developed in CIAM for numerical simulation of fluctuations of scalar parameters was shown. The model supersonic mixing and combustion, which permits to of unmixedness [9] was used to describe the influence of analyze the three-dimensional effects, is presented in concentrations fluctuations on the chemical kinetics. this paper. Some aspects of mixing and combustion enhancement will also be considered on the base of The system of mathematical models and CIAM experience. numerical codes using the boundary layer and parabolized Navier-Stokes equations for two- The scramjet combustor efficiency depends dimensional flows was developed in [lo-121 for the strongly upon the efficiency of the mixing process
Presented at an AGARD Lecture Series on ‘Researchand Development of RadScramjets and Turboramjets in Russia: December 1993 to January 1994. 8-2
between fuel and oxidizer. There are two reasons for mechanisms responsible for the mixing enhancement this effect. The first reason isthat the mixing effkiency may be revealed when using the above mentioned 0 of supersonic flows decreases with Mach number mixing control methods: the increase of the mixing increase because of the influence of compressibility [21- surface; the vorticity generation in the flow; the 231. The increase of the so called convective Mach influence on the stability characteristics; the influence number is accompanied by the essential mixing on the large and small-scale turbulent structures. Ths efficiency decrease in comparison with the mechanisms revealing in the specific control method is incompressible fluid at the same velocities and densities very difficult task and requires using of very sensitive ratios. This decrease observed in the experimental modem measuring techniques and complex investigations for the plane shear layer spreading rates mathematical models elaboration. At the same time it at supersonic convective Mach numbers may be near seems to be justified to use more simple measurements 70% of that value in incompressible case [24]. The and existing mathematical models to predict the similar data for axisymmetrical coflowing jets were efficiency of special mixing control methods and to obtained, for example, in [29]. evaluate the additional losses. Such evaluations for the special mixing control methods make possible the Second reason consists in the decrease of the analysis of the efficiency of the real.injection systems, residence time in the combustor as the flight Mach which include several control methods simultaneously number increases. It is obvious that the problem of as a rule. But the physical investigations of the mixing mixing and combustion enhancement is very important enhancement mechanisms are of primary importance. to create high effective scramjet combustor.
The following factors must be taken into account at estimation of possible mixing and combustion 2.2.Mathematical model and numerical method enhancement methods. It is necessary to evaluate the additional losses connected with the mixing The 3-Dparabolized Navier-Stokes equations are augmentation process and the influence of this losses on used for the numerical simulation of two supersonic the overall performances of the propulsion system such flows mixing. These equations are derived from the full as momentum or thrust [25]. It is desirable to use the 3-D Navier-Stokes equations with exception of the own momentum of hydrogen jet injected into the air second derivatives with differentiation in longitudinal stream as better as possible. This hydrogen has been direction [55]. The equations are written in Cartesian heated by using in the inlet and airframe cooling coordinates X,Y,Z and longitudinal X-axis coincides systems. This circumstance is particularly important for with the air stream direction. The multicomponent flow high flight Mach numbers and makes advantageous the of perfect gases is considered. The heat capacity of axial hydrogen injection. Moreover, the possibility to components depends upon temperature. The chemical use injection struts is doubtful at high hypersonic flight reactions are not taken into account at mixing Mach numbers [26] and then it is necessary to analyze calculations. The turbulent regimes are considered in the opportunity to apply the near wall injection system. the quasilaminar approach. In this case the turbulent viscosity is introduced. The turbulent heat and mass One can distinguish the following control transfer coefficients are expressed through turbulent methods for the mixing enhancement at present : viscosity coefficient using Prandtl and Schmidt a 1) direct influence on the hydrogen jet shape numbers. The differential turbulence model with one with the aid of the hydrogen injector nozzle design or equation for turbulent viscosity is used [56,57]. use of special heads to distort the jet [27-341; 2) the shock-jet or mixing layer interaction [35- The approximate diffusion flame sheet model is 421. Some enhancement mechanisms including the used for the numerical simulation of flows with unsteady effects may exist in this case; combustion. This approach [58,59] is based on the 3) the hydrogen strut design in such manner to assumption that the combustion is controlled only by the create intensive secondaq flows with the convective mixing of fuel and oxidizer. Additionally, it is assumed mixing mechanism [43-471; that one global infinitely fast reaction takes place with 4) the swirling of the hydrogen jet [48 - 491; formation of reaction products from fuel (hydrogen) and 5) the wall jet injection at some angle to the oxidizer (oxygen) and with heat release. The main air stream as an alternative to the slot injection in instantaneous reaction occurs at the surface of the front axial direction. The permissible angles of injection must where the diffusion fluxes of fuel and oxidizer are in be evaluated from the compromise between mixing stoichiometric relation. Only fuel is contained on one enhancement and hydrogen momentum and side of the flame sheet, only oxidizer - on the another . gasdynamical losses [50-531; The products of combustion diffuses in both sides from 6) the use of unsteady effects, including acoustic the flame front surface --I581. influence on hydrogen jet (see, for example, [54]). It is known for the flame sheet model that if 0 It must be pointed out that the following diffusivities of all components are the same, and initial 8-3
concentrations profiles and boundary conditions for The mathematical model presented here is based 0 them are self-similar, it is possible to introduce the on the assumption that the chemical processes are mixture fracture and to solve the single equation for this infinitely fast and the combustion is controlled only by value with following definition of fuel, oxidizer and mixing rate (flame sheet model). Nevertheless there is a reaction products mass fractions taking into account number of scramjet operational regimes with the their initial profiles. residence time in combustor is the same order as characteristic time scale for chemical processes. There The total enthalpy profile may also be obtained are two different cases here. The first one is the "low1' by use mixture fracture profiles with some additional flight Mach number operational regimes (Me-8). As a assumptions. It is assumed that Prandtl and Lewis rule the air flow temperature at the entrance of numbers are equal to unit, and the initial profile for combustor is not sufticiently high for self-ignition at enthalpy is piecewise constant as for the component these regimes. The combustion is possible in this case mass fraction. Also it is necessary to omit some terms only if reliable flameholders and igniters are available. from the energy equation of PNS system, which are not essential in inviscid region and in shear layers. The second case is more important for high flight Mach numbers. The residence time dramatic Thus, in the case of diffusion flame sheet model, decreasing in this case may lead to incompleteness of it is necessary to solve numerically the continuity sufficiently slow relaxation processes, because these equation for the mixture, three momentum equations, processes are controlled by third order slow reactions. and equation for the mixture fracture. Then the total On the other hand these are the processes responsible 0 enthalpy and components mass concentrations are for main part of energy release. Hence the combustion calculated on the base of mixture fracture. inefficiency controlled by chemical kinetics may be observed. In supersonic regions the marching procedure may be used in downstream direction for the numerical It is obviously seen that in both cases the flame solution of PNS equations [55]. In initial cross-section sheet model is not valid and more complex, based on the distributions of all parameters must be detailed kinetics approach, mathematical model is predetermined. The conditions on the lateral boundaries needed. To detect such kinetically controlled operational of computational region depend upon the physical regimes and to avoid the possible errors co~ectedwith problem. Some boundary conditions will be discussed the usage of the flame sheet model, the simplified later. But it is necessary to note that the boundary layers mathematical model for reacting gas flow in combustor on the channel walls with subsonic regions are not was suggested by S.Zaitcev and LBorisov. Here it is. taken into account, because it requires special consideration including global pressure iterations [55]. Let us suppose, that the flow field is one The wall boundary layer effects were neglected because dimensional, and the static pressure is constant in the mixing process is the main object of our combustor. The heat transfer from the reaction zone is consideration. neglected. The variation of the axial velocity U along the combustor is supposed to be negligible (the results of The governing equations are solved numerically calculations show that U remains nearly the same and using the method, which is based on the higher order equals U along the combustor length unlike the Mach accuracy monotone Godunov's finite volume scheme for number decreases dramatically). Hence, the equations the steady supersonic flow [60-62], generalized to three- for species concentrations and the energy conservation dimensional case by A.Kraiko and S.Schipin. This law are to be taken into account explicit predictor-corrector scheme including the modified [63] principal of minimal increments of functions [64] provides the monotonicity condition, the second order accuracy on regular uniform grids, and H=const. approximation on arbitrary nonuniform grids. The generalization of this method on the multicomponent Here, Ci is the species concentrations, Wi - flows taking into account variable heat capacities for chemical production rate, H - total enthalpy including PNS equations was performed by V.Kopchenov and species formation specific heats. The initial conditions K.Lomkov. at the entrance of combustion chamber have to be formulated to complete this problem. Two different sets The proposed method permits to use the adaptive of initial conditions were used for calculations. The first and regular grids. When the mixing layer is thin, the set corresponds to infhitely fast fuel and air mixing adaptive grid is used. This grid is matched to the with conservation of the total mass, momentum and dividing stream surface. When the shear layer becomes energy ( MIX ). The second set of initial conditions thick, the regular grid is used. In this way the accuracy corresponds to assumption that the chemical processes e of the results obtained is increased and the run time is take place at the stoichiometric mixture surface placed reduced. closely to the external boundary of jet. It was assumed 8-4 that .initial temperature and pressure here are equal to uniform. When the jet is injected at an angle to main air the air parameters, and the hydrogen concentration is flow (Fig.2.3.l-b,c), the jet parameters are given as stoichiometric (DIF). The Dimitrov's detailed kinetics uniform, including constant vertical or horizontal scheme was used in calculations [65]. It includes 30 velocity component. At the investigation of the cases C1 elementary reactions between 8 reactants: H,, O,,H, 0, and C,, (Fig.2.3.1-qd) the hydrogen jet parameters in OH, HO,, H,Oz,H20. The specific enthalpies of species initial cross-section are obtained from the calculations were approximated. The Gear-based numerical code of the inviscid supersonic nozzle flow. was used to solve the problem [66]. To investigate possibilities of the mixing There were two values of main interest in the enhancement two series of calculations were performed. calculations. The first one is the Li=U * ti - the The design Mach number for the hydrogen characteristic ignition length scale ( estimated for 5% axisymmetric nozzle was equal to 3.05 and the enlarges of the temperature from its initial value). The temperature T of the hydrogen at the nozzle exit is second one is the L,=U * t, - characteristic length scale equal to 340" K. The air flow parameters for three for energy release (estimated as the length for regimes are shown in table 2.3.1. The ratio of temperature reaches 99% of its equilibrium value). The longitudinal velocity components of hydrogen and air results of calculations are presented in section 2.6. approaches to unit for the second regime. The turbulent viscosity generation mechanism is switched off according to the turbulence model, and the mixing efficiency is determined in this case by the initial 2.3 The main results of the CIAM investigations turbulent viscosity level. Let us consider some results, on sumrsonic mixing and combustion enhancement obtained for regime RI .
The three-dimensional effects due to the The mixing efficiency is evaluated with the aid interaction of the supersonic H2jet with supersonic air of the following characteristics. The value q is flow were analyzed as possible mechanism for mixing introduced as and combustion enhancement [67]. The design of supersonic part of the nozzle for hydrogen injection is q = (1-c n)/(l-). considered as one of the possible passive control methods. The second method includes injection of Here C, is maximum value of H2 mass fraction in hydrogen jet at an angle to the main air stream. cross-section and is massflow averaged H2 mass fraction The mixing efficiency calculations were performed for four cases. In the first case, the hydrogen =jjpuc dF /jjpu dF , was injected through the axisymmetric nozzle with F F conical supersonic part (the reference case Cl). In the second case (Q, the hydrogen was injected as circular where F is the channel cross-section area. This value uniform jet at an angle 10" to the air flow in horizontal varies from zero in the initial cross-section up to one at direction, and in the third case (C,) - at an angle 10" in full mixing. This characteristic permits to determine vertical direction. In the fourth case (C,), the accurately the initial region length and the mixing nonaxisymmetric hydrogen injection nozzle is used. The efficiency on transition and main regions. To evaluate nozzle has elliptical exit cross-section with 2:l aspect the mixing efficiency in the whole region the second ratio. The walls take straight-line form in meridional characteristic is introduced planes. The elliptical and axisymmetric nozzles have the same expansion ratio.
The test problem is formulated in the following manner. It is assumed that the cascade of struts is which defines the nonuniformity of mass fraction placed in the combustor without skew angle. The Hz is distribution in the channel cross-section. The mixing injected in the air flow from the base surface of each efficiency characteristic q, plotted as a function of strut through individual nozzles. The number of nozzles streamwise coordinate X, is shown in Fig.2.3.2. The is large enough. Therefore, the influence of the wall on cwes 1-4 correspond accordingly to investigated cases the "central" jet may be neglected. Approximately, it is c,-c,. possible to divide the flow field into the elementary regions with some symmetry or periodicity conditions The mixing process itself is connected with on the elementary computational region boundaries. irreversible losses. It is necessary to analyze the The computational regions for each of considered cases influence of the mixing control method on the (Cl-C4) are shown in Fig.2.3.1-ad. The direction of additional losses. For example, the additional losses transversal injection component is indicated by pointer. origin may be concerned with shock generation in It is supposed, that the air flow in initial section is supersonic flow. The losses are estimated by the total 8-5
pressure averaged by the massflow. The value < is vortices arising is different. In our case, the longitudinal defined as the ratio of the averaged total pressure to the vortices are the result of three-dimensional fluid same parameter in the initial cross-section. The dynamic interaction behind the nozzle exit cross- dependence of losses upon mixing efficiency 8 is section. This vorticity is not generated in the nozzle presented in Fig.2.3.3. The value 0 is introduced as flow.
In order to illustrate the influence of the mixing enhancement on the combustion efficiency, the same where Do is value D in initial cross-section. It varies regime was calculated using the flame sheet model. The from zero in initial cross-section to unit at full mixing. combustion efficiency as a function of streamwise - The total pressure losses level for ideal full mixing is coordinate is shown in Fig.2.3.6 for axisymmetric designated by a marker. The results comparisons for nozzle (NI), and for two elliptical nozzles with 2: 1 (N2) these four cases confirms that the total pressure losses and 3:l (N3) aspect ratio in the exit cross-section. The - are approximately the same for identical mixing location of the flame front is presented in different efficiency, determined by 0. But it is necessary to note, cross-sections of the channel in Fig.2.3.7. The change of that identical mixing efficiency is achieved on different the flame front surface and its ''interaction'' with distances from initial cross-section (see Fig.2.3.2). computational region boundaries may explain the Hence, the mixing is enhanced for this regime without dependence of the combustion efficiency upon the noticeable additional total pressure losses. The mixing longitudinal distance, and, in particular, the reduction enhancement is also illustrated in Fig.2.3.4, where the of the "combustion rate" at some distance, and also the H2 mass fraction fields are shown for cases C, - C4 in superiority of the nozzle N3 in comparison with nozzles Merent cross-sections. N2 and NI.
The H2 jet injection at an angle 10" is Thus, it is seen, that the three-dimensional accompanied by momentum losses of hydrogen jet of effects, arising due to interaction of the supersonic 1.5%, these losses are equal to 2.2% for conical hydrogen jet from nonaxisymmetric nozzle with axisymmetric nozzle, and for elliptical nozzle - 3.5%. coflowing air stream, provide the mechanisms for the Thus, the additional momentum losses for elliptical mixing and combustion enhancement. The first nozzle in comparison with axisymmetrical nozzle are mechanism is connected with the increase of mixing equal 1.3%, and are approximately the same, as in the and combustion surfaces. The second convective cases C, and C3.Hence, providing the superior mixing mechanism is the formation of intensive secondary level, the mixing enhancement using elliptical nozzle flows with longitudinal component of vorticity. It is does not introduce visible additional losses in necessary to note, that for the mathematical model, comparison with others alternative mixing control presented in this paper, the investigated mixing and methods investigated at this paper. combustion enhancement effects are not connected with additional turbulent viscosity generation. Therefore, it is It is important to reveal the cause of mixing interesting to vellfy the role of the investigated effects enhancement. As follows from the analysis of the for regimes with weak turbulent viscosity generation numerical results, presented in Fig.2.3.4, the hydrogen and simultaneously with low initial turbulent viscosity jet strongly changes its form with increase of "mixing" level. surface in the cases G-C,. The hydrogen mass fraction fields give evidence about mixing enhancement. But for For the regime R2 the turbulent viscosity cases C,, C3, the confluence of neighboring hydrogen generation is small, because the ratio of longitudinal jets takes place at smaller distances from initial cross velocity components approaches to one, and the section than in the case of elliptical nozzle. Moreover, turbulent mixing efficiency is determined by the initial the lateral eruption of hydrogen jet is less intensive in turbulent viscosity level. The dependence of the mixing the case C3in comparison with elliptical jet. efficiency q upon the longitudinal distance for axisymmetric and elliptical nozzles is represented in The second effect, arising at the three- Fig.2.3.8 by curves 1 and 2 respectively. The initial dimensional fluid dynamic interaction of hydrogen jet dimensionless turbulent viscosity level is equal to 3 * 10" and air flow, is concerned with intensive secondary (the ratio of the turbulent viscosity to hydrogen velocity flows formation. The fields of the transverse velocity and radius of axisymmetric nozzle in the initial cross components are shown for all cases in Fig.2.3.5. It is section). The curve 1 shows, that at small initial evident, that the axial vortices are formed as a result of turbulent viscosity levels the turbulent mixing is three-dimensional interaction for cases Q-C4. These unessential. It is obvious, that in this case the elliptical vortices are observed far from the initial cross-section nozzle provides much more superior mixing and provide the convection mechanism for mixing characteristics relative to axisymmetric one. enhancement. sigmiicance of longitudinal vortices for mixing and combustion augmentation was also In order to investigate the influence of the nozzle mentioned for example in [33]. But the nature of design on combustion efficiency for regime R2, some 8-6 calculations were performed for axisymmetric NI and The ramp with angle of 7.1" produced the oblique shock elliptical N2 nozzles with two Merent initial levels of with pressure ratio about 2.5. Computational grid a turbulent viscosity (104 and 2*lO-3) and for elliptical consists of 75x50 cells in cross-section. In the base nozzle N3 with small initial turbulent viscosity level. (reference) case the hydrogen jet Mach number M, was The combustion efficiency as a function of streamwise equal to 2, and jet/air density ratio p' was equal to 0.5. coordinate is shown in Fig.2.3.9. The curves 1,2 correspond to axisymmetric nozzle with small and high The density ratios 0.125, 0.5, 1.0 and 2.0 ~ levels of turbulent viscosity. The similar results for the (corresponding to hydrogen temperatures 1120' K, elliptical nozzle N2 are represented by curves 3,4 and 280' K, 140' K and 70' K ) were studied in the first curve 5 corresponds to the nozzle N3 and low initial series for the inviscid case (Fig.2.4.1 and table 2.4.1). In - - level of turbulent viscosity. The elliptical nozzle this series M, was equal to 2. The peculiarity of the case provides superior mixing characteristics relative to with p'=l consists in the absence of the baroclinic circular nozzle for both initial levels of turbulent source term in the vorticity equation [40], because the - viscosity. In the case of elliptical nozzle N2 with small normal component of density gradient on the jet initial level of turbulent viscosity, the combustion boundary is turned off, and in unviscid case the density efficiency is higher at large distances from initial gradient as the pressure gradient is caused only by wave section, than in the case of circular nozzle with large structure if the numerical scheme dissipation is turbulent viscosity level. Taking into account the negligible. However, the patterns of density fields (see influence of the initial turbulent viscosity level on the Fig.2.4.2, 2.4.3) show, that jet shape deformation combustion efficiency for the elliptical nozzle, it is doesn't depend on p', and remains even if the density possible to assume, that the mixing and combustion gradient is parallel to the pressure gradient. surface increase has the main sigmficance for this regime. Possibly, the convective mechanism, provided In the second series, p' was fixed and equal to by the secondary flows, is less essential in any case for 0.5, but jet Mach number MI was changed: Ml=2,4,6. the higher initial level of turbulent viscosity. The Air stream Mach number was equal to 6. To provide the comparison of curves 1,3,5 shows, that the nozzle N3 identical shock slope in the jet and in the air in the last provides superior combustion efficiency. case, the air and hydrogen specific heat ratios were set equal to 1.4. The rapid interface deformation is proved At the third regime, air flow has high Mach to be reduced with the increase of jet Mach number. number, hydrogen jet is strongly underexpanded, Fig.2.4.4 shows, that jet is only compressed by the besides that the air flow velocity exceeds jet velocity. shock in vertical direction at equal Mach numbers and The comparison was made for three nozzles: NI, N2 and then retains its shape. The non-dimensional circulation N3. One can see from Fig.2.3.10 that the nozzle N3 is presented in Fig.2.4.5 for the investigated cases. It is provides the best combustion efficiency again. interesting to note, that there is no switching of circulation sign with reversing of density ratio, as it Thus, it may be noticed, that the elliptical nozzle would be supposed according to [40]. Moreover, the for hydrogen injection provides the superior combustion essential influence of the density ratio on the circulation efficiency in wide range of regimes in the duct with the level is not observed. But the level of circulation rectangular cross - section. It seems, that it is possible to diminishes with the rise of jet Mach number. So one can choose the nozzle aspect ratio, which provides the best suppose that just the Mach number difference but not level of combustion efficiency for this duct. the density ratio is responsible for the rapid jet distortion.
Some explanation may be proposed on the base 2.4 Some notices about the ietkhock interaction of the detailed pattern of density contours obtained with and mixing enhancement the grid 95x70 for the reference case (see Fig.2.4.6). After the jet and shock touching the part of shock inside The problem of hydrogen jet passage through an the jet "moves" faster than shock in air due to the Mach oblique shock wave was investigated in [39-42,68]. It numbers difference. Besides, the expansion wave arise was pointed out that a strong interface distortion and on the lower part of jet boundary and moves toward the baroclinic torque take place as a result of shock- ramp surface. Probably, these waves result in secondary impingement [39-421 except the possible jet flows with longitudinal vorticity. At some distance, compression [68]. This distortion and baroclinic torque these flows induce jet dividing on two parts and are connected with the generation of the vorticity considerable increase of the mixing surface. The jet located along the fueVair interface [39-41]. dividing effect is clearly discernible in the Fig. 2.4.3.
The goal of this study was to investigate The next series of computations was run for the numerically the process of shock-impingement itself reference regime with turned on viscosity terms to and conditions, providing mentioned effects. Free air obtain information about mixing enhancement: a) for stream conditions corresponded to M2=6, T2=2000'K. free jet without any shock, b) for the base case, when the 8-7
ramp leading edge is in the same section where jet is nozzles for hydrogen injection in each strut, and to injected (Xrmp=O), and c) for Xrmp=30. Density contour the nozzle for hydrogen injection. The injector contours for cases b), c) are shown in Fig.2.4.7, 2.4.8 design, the struts and nozzles disposition must provide and concentration fields in section x=120 - in Fig.2.4.9. the uniform hydrogen jet distribution in the combustor One can see, that in spite of different XrW the resulting cross-section. For example, it is desirable to divide the concentration fields are similar. Only the height of the combustor cross-section on the elementary square jet center is slightly larger in the case of Xmp=O. This regions of the air stream and to inject the circular fact can be explained by longer lifting action of vortices. hydrogen jet to the center of this region with the aid of As for the mixing efficiency, it is seen in the Fig.2.4.10 injection system. The uniform hydrogen jet distribution that the presence of the shock provides the superior in air stream is realized in this case. But the following mixing and the final effect is weakly dependent on the circumstances should be taken into account. The strut place of the shock impingement. system induces the drag, besides it is necessary to solve the strut cooling problem. Therefore, the number of But two others curves in Fig.2.4.10, struts is limited and is comparatively small. Thus, it is corresponding to hydrogenfair mixing in the case of necessary to design the effective injection system with a equal Mach numbers with and without shock, also fixed, as small as possible number of struts. But the demonstrate some mixing enhancement by shock. The dimension of the square elementary region is large in reduction of length, where the mixing efficiency 0.4 is the case of small number of struts. The combustion achieved, constitutes 37% for jet Mach number 6 and length at given parameters in initial cross-section is 63% for jet Mach number 2 (comparisons of lengths are defined by chemical reaction rates and by mixing made for cases with and without shock impingement). It process. If the combustion process is limited by mixing, is diflicult to divide the influence of general then the absolute flame length increases as the hydrogen compression in the shock (the decrease of the jet cross- jet radius or elementary region dimension are increased section area) and jet distortion on mixing enhancement, (at given aidfuel equivalence ratio they are correlated). but the quantitative comparison allows to suppose, that Therefore, really the absolute combustion length is large the jet distortion may provide the additional mixing in the case of uniform distribution of circular hydrogen enhancement. jets for small number of struts.
It is necessary to note that numerical simulation It seems to be useful to increase the number of of shocWjet interaction was performed in CIAM by nozzles on the strut in order to decrease the absolute V.Vasiljev and S.Zakotenko [68]. They point out that combustion length. But when the number of struts is the compression of the jet at processing through the fixed, the elementary air region becomes rectangular. shock is followed by decrease of the jet cross - section The partial confluence of hydrogen jets may cause the area. At the same time, the turbulent viscosity does not decrease of the mixing and combustion efficiency in this change its level after processing through the shock. The case. In fact, as the ratio of hydrogen and air mass flows decreased cross - section area and the unchanged level is fixed, then the distance h between nozzles is of turbulent viscosity at the shocWjet interaction provide decreased as 1/N when the number of nozzles N in strut the main mechanism of the mixing enhancement from is increased. At the same time the axisymmetric nozzle the point of authors view [68]. The secondary flows with radius r is decreased as Us,so that h/r -l/fi+O as vortical structure were observed. But from the authors of N+m. The confluence of hydrogen jets prevents from paper [68] opinion, this mechanism does not influence the effective combustion of circular hydrogen jet in mixing process. rectangular elementary region. This effect was shown in the previous section. But it is possible to use the Thus, from our point of view, the rapid jetlair potential opportunity to decrease the absolute interface distortion and intensive secondary transversal combustion length with increased number of nozzles flows are promoted by the difference of Mach numbers and decreased hydrogen jet absolute dimensions by in the hydrogen jet and air stream, but not only by the choosing the "optimal" elliptical nozzle for rectangular density gradients in the baroclinic torque effect. This elementary channel. As it was demonstrated in the distortion increases the mixing enhancement. The jet previous section, the elliptical nozzle provides the contraction may become the prevailing enhancement decrease of relative mixing and combustion lengths in mechanism at small Mach number difference. comparison with axisymmetric nozzle for the case of increased number of hydrogen nozzles with rectangular elementary channel. The gain in combustion efficiency may be obtained in this case for equivalent absolute 2.5 Evaluation of combustion efficiency and combustor length and identical number of struts in choice of acceptable desiw parameters for combustor comparison with the base injector version, providing the uniform distribution of circular hydrogen jets in the Let us consider the question about the injector cross-section. e design parameters choice. It is necessary to choose the number of struts for hydrogen injection, the number of Thus, it is possible to fix the number of struts 8-8
and to choose the 'number of elliptical nozzles on the increase the number of hydrogen injection nozzles and strut with the "optimal" axis ratio. The previous the number of struts and to decrease simultaneously a investigations showed, that it is possible to increase the absolute dimensions of jets and elementary ducts combustion efficiency with the increasing the elliptical maintaining the uniform distribution of hydrogen jets nozzle axes ratio up to value 3 and more. The over combustor cross section. The corresponding simultaneous variation of the number of nozzles in scheme, illustrating the case of two times greater order to obtain the best number for each ellipticity is number of strut and hydrogen nozzles on each strut is needed too. But such two - parameters "optimization" shown in Fig.2.5.1-b. In this case it is possible to realize procedure seems to be very expensive task, because it is more high combustion efficiency levels on the same necessary to perform direct calculations with two combustor length or to decrease required combustor parameters variation. But some practical limitations length due to reduction of absolute dimensions of jets allow to reduce the volume of calculations. If the and elementary ducts under condition that combustion elliptical nozzle aspect ratio is too large, the excessive process is controlled by mixing rather than by chemical expansion angles in the large axis direction, and kinetics. possible converging of nozzle with circular throat in the small axis direction result in high losses in comparison But under some conditions we can't select this with axisymmetric nozzle. Besides, some design way, because of increase of struts number. On the other restrictions may exist on the nozzle shape and hand the simple increase of nozzles number on strut at dimensions. For example, it may be the overall constant struts number may be ineffective due to quick dimensions restriction. confluence of jets for neighboring elementary ducts with rectangular cross-section. The situation, corresponding e Taking into account the aforementioned to doubled number of nozzles on strut at constant struts circumstances, the following design procedure seems to number is shown schematically in Fig. 2.5.1s. It is be justified. For a given number of struts the base possible to avoid negative influence of jets confluence injector version is considered, providing the uniform for rectangular elementary duct by use of elliptical distribution of axisymmetric hydrogen jets. The injector nozzles elongated in the direction of greater side of with greater number of elliptical nozzles is investigated elementary rectangular duct. The results presented as alternative version. The general hydrogen mass flow above (see, for example section 2.3) show that jet from through the elliptical nozzles is equal to that in the the elliptical nozzle divides into two parts at some axisymmetric nozzles case. The overall dimensions of distance from injection section. This provides the nozzle and strut for the alternative injector do not practically uniform distribution of jets over duct cross- exceed that for base version. The minimal cross-section section, similar to distribution with two times greater of the modified nozzle remains circular, but it has some struts number and simultaneously two times greater times lesser diameter in comparison with base number of axisymmetric nozzles on each strut. The axisymmetric nozzle. The length of nozzle supersonic latter is shown schematically in Fig. 2.5.1- d (the jets, part remains the same in both cases. The supersonic into which breaks down the jet from elliptical nozzle are part of the modified elliptical nozzle has the rectilinear shown schematically by two circles near the ellipse). generatrices. These generatrices connect throat and exit But it is necessary to select axes ratio of elliptical cross-sections. The maximal generatrix inclination nozzles providing the best dividing of hydrogen jets angle is obtained taking into account not only the from the considerations of their mixing and burning. a geometrical considerations but also the desire to provide the acceptable level of nozzle performances. This As an example some combustion efficiency configuration is considered as alternative injector estimations were performed for the model combustor version. schematically presented in Fig.2.5.2. The strut system for hydrogen injection is located in combustor within But simplified method of preliminary injector the symmetrical angle in horizontal plane, which is design may be proposed, which provides efficiency equal to 36O. The number of struts was chosen to be results equivalent to two times larger numbers of struts equal to 15. The air parameters at the combustor and axisymmetrical nozzles on each strut in comparison entrance are presented in the table 2.5.1. The Mach with base injector version. The corresponding schemes, number at the exit of hydrogen nozzle and the illustrating the suggested method are shown in parameters (pressure and temperature) are also shown Fig.2.5.1-a and 2.5.1-d. in table 2.5.1.
The base injector configuration ensures uniform Some injection schemes were tested to compare distribution of hydrogen jets over combustor cross their relative efficiency. In the first case the hydrogen is section, but absolute dimensions of jets and / or injected through the slot in the backward face of the elementary duct falling on one jet are too great to strut. The slot vertical dimension is equal to strut realize acceptable combustion efficiency levels on height. In the second case the high slot is divided into chosen combustor length with chosen number of struts several sections. The hydrogen is injected through the e (see Fig.2.5.1-a) . The most simple solution is to each section at an angle 10" to the main air flow 8-9
direction with alternation of lateral directions of number Mf < 8-9. e injection in neighboring sections. In the third case the injection is realized through the individual The characteristic energy release length scale Le axisymmetric nozzles, which are placed inside strut versus flight Mach number is presented in Fig.2.6.2. It with nozzle exit section on the strut backward face. is seen that the heat release zone length is small in both Then using the number of struts, the number of (MTX and DIF) approaches at Mfi12. The role of hydrogen jets injected from each strut, the elementary recombination processes at Mf increase (Mp12) is region per one jet was obtained in accordance with dramatically increased too and at flight Mach numbers combustor duct geometry (see Fig. 2.5.2). All Mf14-20 these processes can substantially influence estimations were performed for a single jet in the combustor operation. - elementary channel with coflowing air flow. It is obvious that such evaluation method is approximate. The obtained results indicate that the reliable But probably it is possible to compare alternative mathematical model for scramjet combustor has to be - injection schemes efficiency in such manner. developed taking into account the detailed chemical kinetics. Nevertheless, the usage of the developed Let us consider some results about the simplified mathematical model for kinetics effects combustion efficiency, which are presented in Fig.2.5.3. estimations makes it possible to apply the flame sheet It is necessary to note, that the injector with individual model in scramjet combustor flowfield calculations nozzles was chosen to provide the uniform distribution more argumentally. of hydrogen jets in the injector control cross-section. 0 Therefore the number of nozzles is equal to 20 on each strut. The curves 1,2,3 in the Fig.2.5.3 show the combustion efficiency as the function of longitudinal coordinate for test cases 1-3 (correspondingly for slot, section slot and axisymmetric nozzles injection schemes). It is necessary to note the low combustion efficiency in the cases of slot and section slot injection. The better results are obtained in the third case for injection with the aid of individual axisymmetric nozzles.
The proposed scheme of combustion enhancement due to doubled number of elliptical nozzles at the same number of struts was also tested. In this case the elliptical nozzles with aspect ratio 2:l in exit cross-section were used. The results obtained in this case are shown in Fig.2.5.3 by curve 4. The obtained results demonstrate the possibility to reduce the necessary combustor length with the aid of this injection system with doubled number of elliptical nozzles. So, if the combustion efficiency value 0.9 is reached on the distance approximately one conventional length unit for injector with elliptical nozzles, it is needed approximately 2.14 units for the injector with axisymmetric nozzles, and the distance more than 2.85 units is necessary to reach the combustion efficiency level 0.9 in the case of section slot and slot injectors.
2.6 Estimation of the role of chemical kinetics
To estimate the role of chemical kinetics at the Merent scramjet operational regimes, the set of parameters at the entrance of combustor presented in table 2.6.1 was used. The characteristic ignition length Li versus flight Mach number is presented in Fig.2.6.1. It is seen that L~is sufficiently smaller then the typical combustor length scale (-lm), and the self-ignition can be the limitation mechanism only if the flight Mach 8-10
References
1. EaeB B.K., rOJIOBkIYeB B.M., HCaKOB Typ6yJIeHTHbIX pealXpyIOlQEar CTpyji)l, QWHKa B.A., "flByhlepHbIe ~yp6yne~~~b1eTeYeHm TopeHm EI B~P~IB~,N 3, 1987, c.85-88. peaIupyIoum WOB", HOBOCH~EI~CK,HayKa, 13. KOJIeCHMKOB O.M., "PaCYeT IIOBbIlIIeH~ 1976, 264 c. ~aBJIe"H IIPEI lWPeHHEi IIJIOCKOG CBepX3ByKOBOfi 2. EaeB B.K., rOJIOBEiYeB B.H., TpeTbHKOB CTpyIl BOnOpOna B CBepX3BJXOBOM ClTJTHOM n.K., TapaHHH A.@., KOHCTaHTEZHOBCIMfi B.A., noToKe", YYeHbIe 3an~cmI&WM, T.13, N6, . RCaKOB B.A., "rOpeHEie B CBepX3BJXOBOM 1982, c.49-58. noToKe", HOBOCH~EIPCK,Haw, 1984, 302 c. 14. KOneCHEIKOB O.M., "BnmHkie 3. MelQepHKoB E.A., Ca6eJIJdlEIKOB B.A., IIJVIbCauEIfi KOHI&HTpaqEIEf Ha BOClJJIaMeHeHRe "rOpeHEIe BOnOpOna B CBePX3BJXOBOM npEICTeHHOfi CTpyEl BOnOpOna B CBepX3ByKOBOM Typ6yneHTHOM IIOTOKe B KaHaJIe IIPS ClTJTHOfi nomKe", @~13mI'opeHm EI B~P~IB~,N 1, nonaye ropmqero EI OIMC~~JIII", QEI~EIK~ 1985, c.53-57. TopeHm EI BpbIBa, N 2, 1981, c.55-64. 15. KOneCHEIKOB O.M., "BnmHMe 4. 3mOHT B.n., neBElH B.M., MeWepHKOB EiHweKqm Bonopona Ha mamoAeiic-mie E.A., Ca6enb~a~o~BA., "OCO6eHHOCTII HepaBHOBeCHOID IIOI'paHEiYHOrO CJIOH C ropesm HenepeMeluaHHbrx WOB B K~H~TI~X", BHeluHElM CBepX3BJXOBbIM nOTOKOM", YYeHbIe @EI~HK~Topeam EI B~P~IB~,N 4, 1983, c.75-78. 3an~c~1.rI.@I'M, T.18, N 6, 1987, c.101-106. 5. MeI4epHKOB E.A., Ca6eJIbHEIKOB B.A., 16. White M.E., Drummond J.P., Kumar A., "POJIb CMelIIeHHH EI KIIHeTElKEI B yMeHblIIeHHEi "Evolution and Application of CFD Techniques for 0 TelJJIOBbZqeJIeHWI npEI CBePX3BJXOBOM IOpeHEIEi Scramjet Engine Analysis", J. of Propulsion and Power, HenepeMemmHm OB B pacmpHmwmcH V.3, N 5, pp.423-439 K~H~JI~X",@~3m ropeHm EI B~PMB~,N 5, 17. Barber T.J., Cox G.B., Jr., "Hypersonic 1988, c.23-32. Vehicle Propulsion: A Computational Fluid Dynamics 6. rpoMoB B.T., nap" O.E., Application Case Study", J. of Propulsion and Power, "TaHD3H~aJIbHbIfi WyB BOnOpOna B V.5, N 4, 1989, pp.492-501. Typ6y~1eHTHbIfi CBepX3BJXOBOfi IIOI'paHEiYHbIfi 18. McClinton C., Bittner R., Kamath P., "CFD cnoii", M3~ecm AH CCCP, Mexmma Support of NASP Design", AIAA 90-5249,1990. XWJKOCTEIEI Tma, N 3, 1982, c.189-191. 19. Anderson G., Kumar A., Erdos J., "Progress in 7. rpOMOB B.T., nap= O.E., neBElH B.A., Hypersonic Combustion Technology with Computation "npHMeHeHHe MOneJIH "HeCMelIIa"0CTEI" and Experiment", AIAA 90-5254,1990. pacnem TY~~YJI~HTHO~~~PEICT~HHO~~ clpp 20. Rogers R.C., Chinitz W., "Using a Global BOnOpOna B CIIJTHOM CBePX3BJXOBOM IIOTOKe Hydrogen-Air Combustion Model in Turbulent Reacting ~o3nyxa",XIlMmecKarr QEZ~EIK~,T.3, N 8, 1984, CalCulatiOnS", AIAA J., V.21, N 4, 1983, pp.586-592. c.1190-1195. 21. Papamoschou D., Roshko A., "Observation of 8. rpoMoB B.T., napmi (I.E., JeBm B.A., Supersonic Free Shear Layers", AIM 86-0162, 1986. "Typ6yne~~~oeropeme Bonopona B 22. Bogdanoff D. W., "Compressibility Effects in ~~EICT~HHO~~clpye, mxemqeii B C~TJTHJ~ Turbulent Shear Layers", AIAA J., V01.21, N 6, 1983, CBepX3BJXOBOfi IIOTOK B03nyXa", @EI3EIKa pp. 926-927. 0 Topeam EI B3pm, N 6, 1987, c.3-9. 9. Spiegler E., Wolfshtein M., Manheimer- 23. Samimy M., Erwin D.E., Elliott G.S., That Y., "A model of unmixedness for "Compressibility and Shock Vawe Interaction Effects on turbulent reacting flows", Acta Astronautica, V.3, Free Shear Layers", AIM 89-2460, 1989. N 3-4, 1976, p. 265. 24. Gilreath H.E., Sullins G.A, "Investigation of 10. EapmoBcm6 C.M., nesm B.M., the Mixing of Parallel Supersonic Streams", ISABE 89- HWBO~C~~A.C., TypEiqeB A.M., 7060,1989. "3KCnepEIMeHTanbHO-TeOp~eCKOe EICCJIe- 25. Swithenbank J., Eames I., Chin S., Ewan B., n0Ba"e CrryIlibIX CBepX3BJXOBbIX PeallipyIo~ Yang Z., Cao J., Zhao X., "Turbulent Mixing in c~pyjil', C~OPHEIK"Typ6yne~~"e TeYeHm B Supersonic Combustion Systems", AIAA 89-0260, peampyIom IIOTOK~X", MHHCK,Hap, 1986, 1989. c.84-98. 26. Hewitt F.A., Ward B.D., "Advanced 11. EapmoBcmfi C.H., HWBOPCIM~~A.C., Airbreathing Powerplant for Hypersonic Vehicles", IlepMmoB B.A., "PacqeT ropemw ISABE 89-7107, 1989. TY~~YJI~HTHO~~CBePX3BJXOBOG HeEI306apEiYeCKOfi 27. Gutmark E., Wilson K. J., Schadow K.C., Parr CTpyEi BOAOpOAa B ClTJTHOM CBepX3BJXOBOM T.P., Hanson-Parr D.M., "Combustion Enhancement in nmKe BO~A~X~",@EI~EIK~ ropeHm EI &pbIBa, Supersonic Coaxial Flows", AIAA 89-2788, 1989. N 4, 1986, c.14-18. 28. Schadow K.C., Gutmark E., "Rewiew of 12. EapmoBcKufi C.M., IIepmHoB B.A., "0 Passive Shear Flow Control Research for Improved IIpEiMeHEiMOCTEl MOneJIEi AH@@y3MOHHOrO Subsonic and Supersonic Combustion", AIAA 89 -2786, a ropeam K pacveTy CBepX3BJXOBm 1989. 8-11
29. Schadow K.C., Gutmark E., Wilson K.J., Combustion in Scramjets", AIAA 91-2394, 1991. 0 "Passive Mixing Control in Supersonic Coaxial Jets at 47. Riggms D. W., McClinton C.R., Rogers R.C., DBerent Convective Mach Numbers", AIAA 89-0995, Brittner R.D., "A Comparative Study of Scramjet 1989. Injection Strategies for High Mach Number Flows", 30. Schadow K.C., Gutmark E., Koshigoe S., AIM 92-3287, 1992. Wilson K.J., "Combustion-Related Shear-Flow 48. Naughton J., Cattafesta Z., Settles G., 'I An Dynamics in Elliptic Supersonic Jets", AIM J., Vo1.27, Experimental Study of the Effect of Streamwise N 10,1989, pp.1347 -1360. Vorticity on Supersonic Mixing Enhancement", AIAA 31. Schadow K.C., Gutmark E., Wilson K.J., 89-2456, 1989. "Noncircular Jet Dynarmcs in Supersonic Combustion", 49. Naughton J.W., Settles G.S., "Experiments on - J. of Propulsion and Power, Vo1.5., N 5, 1989 , pp.529- the Enhancement of Compressible Mixing via 533. Streamwise Vorticity, Part 1-Optical Measurements", 32. Gutmark E., Schadow K.C., Bicker G.S., "On AIAA 92-3549, 1992. the Near Acoustic Field and Shock Structure of 50. Mays R.B., Thomas R.H., Schetz J.A., "Low Rectangular Supersonic Jets", AIAA 89-1053, 1989. Angle Injection into a Supersonic Flow", AIAA 89- 33. Gutmark E., Schadow K.C., Parr T.P., Parr 2461,1989. D.M., Wilson K.J., "Combustion Enhancement by Axial 51. Fuller E.J., Mays R.B., Thomas R.H., "Mixing Vortices", J. of Propulsion and Power, Vo1.5, N 5, 1989, Studies of Helium in Air at Mach 6", AIM 91-2268, pp. 555-559. 1991. 34. Zaman K.B.M.Q., Reeder M.F., Samimy M., S2. Fuller E.J., Walters R. W., "Navier-Stokes '3upersonic Jet Mixing Enhancement by "Delta-Tabs" Calculations for 3-D Gaseous Fuel Injection with Data 'I, AIAA 92-3548, 1992. Comparisons", AIAA 91-5072, 1991. 35. Drummond J.P., Carpenter M.M., Mukunda 53. Wood C.W., Thomas R.H., Schetz J.A., H.C., "Mixing and Combustion Enhancememt in "Effects of Oscillating Shock Impingement on the Supersonic Reacting Flows", Proceedings of the Fah Mixing of a Gaseous Jet in a Mach 3 Air Stream", International Symposium on Numerical Methods in AIAA 90-1982,1990. Engineering, Vol. 1, Springer-Verlag, 1989, pp. 175-184. 54. Kumar A., Bushnell D.M., Hussaini M.Y., "A 36. Drummond J.P., Mukunda H.S., "A Numerical Mixing Augmentation Technique for Hypervelocity Study of Mixing Enhancement in Supersonic Reacting Scramjets", J. of Propulsion and Power, Vo1.5, N 5, Flow Fields", AIAA 88-3260, 1988. 1989, pp.514-522. 37. Shau Y.R., Dolling D.S., "Experimental Study 55. AHnepcoH fl.,TaHHexm Am., nnewep of Spreading Rate Enhancement of High Mach Number P., "BbIYEicnIiTeJIbHaH WpoMeXaHEiKa EI Turbulent Shear Layers", AIM 89-2458, 1989. T~~OO~M~H",T. 1-2, MocKsa, MH~,1990. 38. Menon S., "Shock-Wave-Induced Mixing 56. Kosno~ B.E., CeKyHAoB A.H., Enhancement in Scramjet Combustors", AIAA 89-0104, CMEIpHOBa kI.n.,"MoAenH TYP6YJIeHTHOCTEi HJIII 1989. onmaHm TeYeHm B cTpye CmMaeMoro ma", 39. Marble F.E., Hendricks G.J., Zukoski E.E., kI3BeCTm AH CCCP, MexaHHKa &ll'(KOCTH H "Progress Toward Shock Enhancement of Supersonic Tma, N 6, 1986, c.38-44. 0 Combustion Processes", AIAA 87-1880, 1987. 57. K03JIOB B.E., "YlTpO~eHHbIfi MeTOA 40. Waitz I.A., Marble F.E., Zukoski E.E., YHcneHHoro pacveTa TpexMepHoii Typ6y~1eHTHOfi "Vorticity Generation by Contoured Wall Injectors", CBO6OnHO6 CTpyEI", kI3BeCTm AH CCCP, AIM 92-3550, 1992. MexaHEiKa ~KOCTMH Tma, N 1, 1989, c. 170- 41. Waitz LA., Marble F.E., Zukoski E.E., "A 173. Systematic Experimental and Computational 58. 3eJIbAOBlW B.E., Eape~6nan mi., Investigation of a Class of Contoured Wall Fuel JIEIGPOBEIY B.E., Maxsmm3e T.M., Injectors", AIAA 92-0625, 1992. "MaTeMaTEIYeCKaH mopm ropeam w mpmn", 42. Drummond J.P., "Mixing Enhancement of MOCKB~,HayKa, 1980, 478 c. Reacting Parallel Fuel Jets in a Supersonic Combustor", 59. KY3HeuOB B.P., Ca6eJIbHEiKOB B.A., AIM 91-1914,1991, "T~~~YJI~~HTHOCT~EI ropeakie", Mocma, HayKa, 43. Northam G.B., Greenberg I., Byngton C.S., 1986, 287 c. "Evaluation of Parallel Injector Configurations for 60. rOAJWOB C.K.,3a6po~a~A.B., kIBaHOB M.H., KpaiiKo A.H., ~POKOIIOB Supersonic Combustion", AIAA 89-2525, 1989. r.n.," %iCJ'IeHHOe PelueHEie MHOD3MepHbIX 3waY 44. Drummond J.P., Carpenter M.H., Riggins mOBO% AEIHaMEIKM", Mocma, HayKa, 1976, 400 D.W., Adams M.S., "Mixing Enhancement in Supersonic Combustor", AIAA 89-2794, 1989. 45. Riggins D.W., McClinton C.R., "A Computational Investigation of Flow Losses in a Supersonic Combustor", AIAA 90-2093, 1990. 46. Northam G.B., Capriotti D.P., Byington C.S., Greenberg I., "Mach 2 and Mach 3 Mixing and 8-12
MaTeMaTmeCKOfi @Ei3EiKki7 T.23, N 4, 1983, c.848-859. 62. POAHOHOBA.B., "MOHOTOHH~cxeMa moporo nopqqm annpoKcEiMaqEiEi ~nrr CKB03HOrO PaCYeTa HepaBHOBeCHbIX -Te'-EHEi~', XYpHa BbIYEiCn&fTt?JIbHofi MaTeMaTHm Ei MaTeMaTmecKoii @Ei3Hm, T.27, N 4, 1987, -I c.585-593. 63. TEylruIesa H.B., "0606we~~e MO&H@HuHpOBaHHOfi CXeMbI C.K.rOAyHOBa Ha IIpOCTpaHCTBeHHO HeOAHOpOAHbIe CeTKEI", YyeHbIe 3an~icmI.&4M, T.17, N 2, 1986, c.19- 26. 64. Konm B.n., "KOHeSHO-pa3HOCTHm CXeMa AIM PaCYeTa AByMepHbIX Pa3pbIBHbIX HecTagHoHapHm ra3onEiHaMmecmx pemeaEiii", YqeHbIe 3an~cm&41Ti7 T.6, N 1, 1975, c.9-14. 65. Dimitrow V.I., "The maximum kinetic mechanism and rate constants in the H2-02 system", React. Kinetic Catal. Lett., V.7, No.1, 1977, pp.81-86. e 66. IIonaK JI.C., rom~e~6eprM.X., JIeBEigKkifi A.A., "BbmEiCJIEITeJIbHbIe MeTOAbI B XHMEIYeCKOfi ICHHeTEiKe", MoCKBa, Hap, 1984, 280 c. 67. Kopchenov V.I., Lomkov K.E., "The Enhancement of the Mixing and Combustion Processes Applied to Scramjet Engine", AIM 92-3428, 1992. 68. BaCHnbeB B.kl., ~~KOT~HKOC.H., "06 EIHT~HCH@E~K~UIIIcMemeHm IIOTOKOB 3a KocbiM CKaYKOM ~OTHeHkUI", kI3BeCTkUI PM, MexaHEiKa ~KOCTH~i lka, N 2, 1992, c.61- 68. 8-13
TABLE 2.3.1
AIR STREAM PARAMETERS:
Re&e T2 "K M2 p2/p1 U2h
1050 2.8 0.98 0.41 R1
1950 5.5 5.88 1.10 R2
R3 2300 6.6 6.37 1.35
TABLE 2.4.1
AIR STREAM M2~6.0, T2=2000" K
series case P1P2 Ml Ti OK Ul'U2
1 1.1 2.0 2.0 70 0.25
1 1.2 1.0 2.0 140 0.35
reference 8-14
TABLE 2.5.1 I air stream parameters hydrogen parameters
1.5*105 Pa 425OK 2.6 1 1
fuel - air equivalence ratio 4=0.98
TABLE 2.6.1
8 12 16 20 Mf
4 0.91 1.66 5.0 5.0
1100 1500 1970 2300 T& OK
Pair Pa 1.8*105 1.2*105 0.8*105 0.65*105
Mair 3.05 4.8 5.5 6.6
n=PH,/Pk 0.84 1.92 5.88 6.37
2.9 2.9 2.9 2.9
THz "K 373 3 73 373 373 8-15 0 0COHPUTRT ON I ONRL REG I 0
A) INITIALCROSS-SECTION FOR CASE C, 0 0 COtlPUTRTIONRL REGION 1 I e
B) INITIALCROSS-SECTION FOR CASE ,C2
COnPUlRTIONAt REGION 0 0 0
c> INITIALCROSS-SECTION FOR CASE C, a 0 0 COtlPUTRTIONRL REGION 0
D> INITIALCROSS-SECTION FOR CASE C,
Fig.2.3.1 INITIAL CROSS-SECTIONS AND COMPUTATIONAL REGIONS FOR CASES c,-c, e 8-16
MIXING EFFICIENCY CHARACTERISTIC 9 FOR REGIME, R,
Fig. 2.3.2
THE DEPENDENCE OF TOTAL PRESSURE LOSSES UPON MIX)N(3 EFFIC1ENCY Fig -2.3.3 8-17
AXISYMMETRICAL NOZZLE HORIZONTAL INTERACTION
e
A) Hi MASS FRACTION FIELDS FOR C SE c, B) H, MASS FRACTION FIELDS .FOR CASE c,
VERTICAL INTERACTION ELLIPTICAL NOZZLE
D) Hz MASS FRACTION FIELDS FOR CASE c,
C) H, MASS FRACTION FIELDS FOR CASE c,
H, MASS FRACTION FIELDS FOR'CASES c,-c,
Fig.2.3.4 8-18
I ...... 1 I ......
I' I
...... e. .. .I I ...... tell ...... Illll 'I I1 \I \I \I
...... " * ' a'.' ......
...... n 8-19
WR
COMBUSTIONEFFICIENCY FOR NOZZLES NI, N, AND N, WITH ASPECT RATIO $1 21 AND 3:l
Fig .2.3.6
Fig. 2.3.7 8-20
a I I I I 140 160 ibo tbo $0 240
WII
MIXINO EFFICIENCY FOR REGIMER;!
Fig.2.3.8
C
COMBUSTIONEFFICIENCY FOR REGIME R,
Fig.2.3.9 8-21
0.. 4 C I RCLE 00 ---- ELLIPS 2:l -.-.- ELLIPS 3:l
(D l- .
3-.
cu.
a . I I I I I I I ( Sb 1 bo 150 2b0 250 3b0 3kO X/R
Fig.2.3.10, COMBUSTION EFFICIENCY FOR REGIME R,
l I e- 2
Fig.2.4.1 THE SCHEME OF SHOCK - JET INTERACTION IN CHANNEL 8 -22
SHOCK - JET INTERACTION L -1e.0 0 ” I\ R/RH-CONST
b
DENSITY FIELD e
tl14 -6 Rl/ft2=0.5 SERIES I
Fig.2.4.2 - A
SHOCK - JET INTERACTION
R/RH=CONSI e
DENSITY FIELD e 8-23
SHOCK - JET INTERACTION
DENSITY FIELD
n1=2 R2-6 Rl/R2=2 SERIES 1
Fig.2.4.2- C
R/RH-CONST
DENSITY FIELD
Fig.2.3.3- A 8-24
SHOCK - JET INTERACTION
DENSITY FTELD
Fig.2.4.3 - B
R/RH=JMNSI
IN- - I\ I DENSI'IY FIELD
tll=6 H2-6 RI/W.S
Fig.2.4.4 8-25
THE INFLUENCE OF THE DENSITY RATIO AND e MACH NUMBER DIFFERENCE ON THE TOTAL
THE DETAILED STRUCTURE OF THE SHOCK - JET INTERACTION 10.0 I R/RH=CONSl xs Q*0 I DENSITY FIELD a
e Hlt.2 tl24 Rl/R2=0.5 Fig .2.4.6 SHOCK - JET INTERACTION
R/RH=CONST
DENSITY FIELD 0
Fig. 2 -4.7
SHOCK - JET INTERACTION
PIRH-CONS
DENSITY FIELD 8-27
HYDROGEN MASS FRACTION FXELDS
1'9.5 .- RRTION CF 1-51 (;Rs x- 120.00 75r50 A) , x-=o
.11.s+ 15. OJ J wm aw .)a00
Fig .2.4.9
INTERACTION OF VISCOUS JET AND SHOCK
0
-- tllu2 r12..6 NO SHOW I-? ------tll=ll2=6 SHOCK
X/R 8-28
Fig.2.5.1
INJECTlON SCHEMES.
THE COMBUSTOR SCHEME
THE ELEMENTARY COMPUTATIONAL REGIONS FOR VARIOUS INJECTION SCHEMES
Fig .2.5.2 8-29 a
Fig.2.5.3 THE COMBUSTION EFFICIENCY FOR VARIOUS INJECTION SCHEMES 8-30
Li (meters)
1...1...1.-.1- 0.1 2 - ...... :.I/- -I 0.08 ...... -I
-. (...... r ‘I\ -i \ 0.04 -. (...... ~.....
-. (......
0 -...... I... Mf 8 12 16 20
Fig.2.6.1. THE CHARACTERISTIC SELF - IGNITION LENGTH V.S. FLIGHT MACH NUMBER
Le (meters)
c .... ! ... ! .. Y , . IT 3-...... Ti - i......
2- ......
- ......
1- ...... , ...... ,.-
- ...... ----+------0- ...... - - I..., Mf 8 12 16 20
Fig.2.6.2. THE CHARAclTERISTIC HEAT RELEASE LENGTH V.S. FLIGHT MACH NUMBER 9-1
SOME PROBLEMS OF SCRAMJET PROPULSION FOR AEROSPACEPLANES PART II. -- SCRAMJET: DEVELOPMENT AND TEST PROBLEMS by Dr A.ROUDAKOV CIAM (Central Institute of Aviation Motors) 111250 MOSCOW RUSSIA
Scramid development stratem. and temperature factor T=TWrr, ate presented in Fig. 2.3 in coordinates f/, h$ Wall temperature of scramjet Scramjet is considered as the base of single duct is taken to equal 1000K. One can see, that stage aerospace plane propulsion system providmg its Reynolds number is approximately constant along the acceleration from moderate supersonic (MI= 5) up to flight trajectory. When characteristic length of duct high hypersonic flight speeds (MJ= 25). element L is about IO m, real Reynolds number for all duct elements is very high, He= IO?. . I@. The feature of scramjet scheme seems to be simple and comfortable for propulsion system and Since Mj=12, temperature factor Tw values vehicle integration. The method of scramjet specific become lesser than that in liquid propellant engine parameter's computations seem to be also simple due to combustors, Tw.-.0.2. In this case, overcooling wall supersonic flow in scramjet (M > I) and configuration layer corresponds to high enthalpy core flow. Besides, simplicity of scramjet duct and its elements: inlet, flow in combustor and exhausts nozzle in characterized combustor and exhaust nozzle. by high level of turbulence generated by shock waves. The structure of these shock waves' changes along the But preliminary analyses of computation re- flight trajectory and is determined by flight Mach sults show high difficulty of specfic parameter's number Mfand engme operating regme. scramjet realization and rigorous demands to its design, operating process realization and its simulation in test The dependence of inlet throat Mach number facilities. Aerospace plane flight scheme and main MHon the flight Mach number is shown in Fig. 2.4. parameters of flow in its duct along the flight path pre- The values of Mach number at the exit of hydrogen sented in Fig. 2.1 ...2.4 illustrate it's fact. Conventional combustor with constant cross section area at air-fuel flight scheme is presented in Fig. 2.1 in coordinates: ratio a=I.l M* are presented in the same fig.2.4. It is flight altitude H (km) and flight Mach number n./J = evident, that flow Mach number in such combustor ///U. This flight scheme is limited in vertical direction remains approximately constant at given flight speed, by lines of maximum (qmm= 2bar) and minilnuin changes strongly along the flight path and reaches value (qmin=0.2bar) air ram value, and it is limited in hori- M* = 4...5 at high flight Mach number Mf 15. zontal direction by minimum (Mmin = 5) and maximum (Mmax = 22) Mach number values at which Scramjet specific parameters depend on gven scramjet specific impulse exceed specific impulse of flight path optimized with talung into account restric- ramjet and liquid propellant engine respectively. The tions concerned, in particular with ecology problems. grid of constant total enthalpy and of engine entrance The restriction due to condition of minimum distur- pressure is presented in the same Fig. 2.1. One can see, bance of atmospheric ozone layer is shown in Fig. 2.1 that total enthalpy change along the trajectory more by section-lined belt. than 20 times and reaches 20 MJIkg, that rouses large problem of scramjet design and creation. Scramjet specific impulse this shown along the
flight path corresponding to y : 0.75 at a==/./ in Fig. The lines of constant values of quasi equilibri- 2.5 with taking into account different components of um total temperature To and equilibrium total pressure losses. The losses of specific impulse are taken into Po and the engine entrance computation with the as- account by integral factors: inlet total pressure recovery sumption of thermodynamic equilibrium isotropic com- factor c=uw),combustion efficiency factor qc nozzle pression process are presented in Fig. 2.2 in the same impulse factor qn and nozzle under expansion factor n coordinates H and A4f These parameters correspond to = PJPpI. Solid line 1 shows calculated specific air parameters in high pressure cylinders that can be impulse of hydrogen scramjet with taking into account used at tests of scramjet duct and its element in test real losses of total pressure in inlet (u = @v),'losses facilities. One can see, that since 19=8total tempera- due to combustion incompleteness ( qc=O. 95) and ini- ture To>2500K and total pressure P&500bar exceed pulse losses in nozzle (pn=0.97) at nozzle expansion the possible values of these'parameters that can be re- ratio Fa=3. Solid line 5 presents specific impulse of alized stationary at earth conditions. ideal scramjet. The rest of dotted lines 2,3,4 show the The lines of constant value of the Reinolds values of coinponents of specific impulse losses, de- number calculated for unit length of main streani (Re,) scribing perfectinn nf ccramiet duct elements. Influence Presented at an AGARD Lecture Series on 'Research and Development of RadScramjets and Turboramjetsin Russia: December 1993 to January 1994. 9-2
coefficients of every component of losses Ki, or par- tems. These forks create the basic of knowledge, data, ticular logarithmic derivatives of specific impulse J on and experience for new stage of scientific and technical appropriate loss's factor I;. (Ki=(hJ)G) are shown by advance, named conventionally hypersonic technology deshed lines. It should be noted that losses in exhaust in parallel with solution of main problem (scramjet nozzle and other's scramjet duct elements have strong development). effect on engine specific impulse and this fact imposes special requirement for physical or numerical simula- The second lead is connected immediately with tion of scramjet operating process along the flight path. development and realization of scramjet starting with its pr.eliminary design and including tests of engine, that Thus, scramjet design and development de- determines lock of aerospace transport of XXI-st mand realization of operating process in its duct with century. flow parameters in the range that can be realized in test conditions stationary only at Mf< 8. Reliable data on And finally, the third lead of works is con- physical chemical conversions in flow combustion, tur- nected with creation of industrial base for tests of full bulence, three-dimensional and unsteady structure of scale units of scramjet, their development and develop flow wall flows and heat exchange and their coupling ment of engine in general, including full scale service- are absent for high Mach number at present. The ex- life tests of scranijet in real range of parameters. tension of scramjet operating domain and improvement of scramjet performances is connected with duct com- On one hand, all three leads are sufficiently plication and in particular with transition to multimode independent, and on the other hand, they interact es- ducts with variable geometry, with use of gasgenerator sentially for a long time. This fact causes the necessity like liquid propellant engine and so on. For this reason, of coordination of effects in each of mentioned leads. the problems of new materials creation and development The conception works on scramjet initiate works on of active heat protection system of structure are of leads 1' and 3. But new results in field of hypersonic important significance. technology, the time of appearance of which can't be predicted, can change essentially assumed design of Thus, in spite of apparent simplicity of scramjet. Long time high-cost design and development scramjet scheme, the development of this engine is of test facilities can not give expected results without extremely complex problem and solution is connected taking into account this fact. with working out some fundamental and applied scien- tific leads, methods of investigation and simulation of Hypersonic technology (see Fig. 2.7) includes 6 scramjet operating process scientific and methodical parts: accompaniment of general conception of scramjet, which is fully integrated with vehicle and determines -- fundamental scientific investigations; preliminary design of aerospace plane in general. -- applied scientific studies of scramjet operating process; If we take into account mutual interaction of -- experimental siniulation of engme process; considered problem elements, time and cost, which are -- numerical simulation of operating process; necessary for developments, the possibility of appear- -- new technical approaches and systems; ance of new results, which can change former concep- -- control systems and units of scranijet. tions and technical approaches, it seems to be expedient to have scheme of scientific research and design works 1. Fundamental scientific investigations (Fig. for scramjet development. This scheme is presented in 2.8) will cover physical-chemical processes in high- Fig. 2.6. enthalpy flows with shock waves, new materials creation and ecology problems. It seems to be very The whole of works, presented in Fig. 2.6 are important to obtain new theoretical and experimental divided into three different directions: data on kinetics of chemical reactions in high-enthalpy flows with shock waves, detonation and combustion on -- hypersonic technology turbulence structure in such flows and also on stability -- scramjet conceptions and control of these processes. -- test base. These data are the necessary data base for The first of these directions is the most exten- mathematical simulation of processes of propagation of sive. It is the complex of fundamental and applied sci- Qfhsion and detonation combustion in supersonic entlfic researches and design works with common pur- combustors and also for practical realization of energy pose scramjet creation. It covers field of physical- supply to supersonic flow. chemical processes in high enthalpy supersonic turbu- lent flows, their physical and numerical simulation and New material problem is the key problem for methods of measurements and diagnostics, instruments scranijet development because the use of existing and materials and appropriate technical equipment and sys- promising materials doesn't provide scramjet flight at 9-3
high flight Mach numbers n/’S ;, 16 due to restrictions 0 caused by structure thermal state. 4. Numerical simulation of scranijet operating process permits to generalize the whole experimental Ecology problem of hig.1-speed motion in at- data accumulated in test studies and interpolate they on mosphere and, in particular, in ozone layer demands real conditions of engme operating in test facility and in also in-depth study. flight.(Fig 2.1 I)
2. Applied scientific studies of scramjet oper- Numerical simulation of operating process ating process (Fig. 2.9) will cover new sections of tech- should be based on mathematical model of turbulence in nical sciences, connected with selection of design high-enthalpy hypersonic flow of reacting gas mixture scheme and parameters of scramjet inlet, combustor and with talung into account the absence of equilibrium of exhaust nozzle and with organization and control of physical-chemical conversions. Supercomputers are operating processes in them in a wide range of de- necessary for computations of three-dimensional and termining parameters and conception of reacting gas unsteady interaction of shock waves with channel mixture. boundary layer determining thermal state of structure. The problems of flow stability at combustion in bounded The development of computational model of supersonic flow in essential in connection with operating process, of methods of duct unit’s computa- combustor operating process control. tion, representing the basis for technical approaches is connected with investigations in the fields of thermo- 5. New technical approaches (Fig. 2.12). dynamics, gasdynamics, thermal physics and combus- Scramjet realization is impossible without creation of tion, which are relied on new methods of measurements new technical approaches. It should be noted the fol- of supersonic reacting gas flow parameters and on lowings: results of numerical processing and simulation of three- - the control of three-dimensional gasdynamic dimensional and unsteady phenomena of various scale. structure of supersonic flow in duct with the aim of its Wide use of super computers is necessary for stabilizing, avoiding boundary layer separation and mathematical simulation of scramjet operating process. eliminating heat flux peaks;
3. Experimental simulation of scramjet operat- - the control of combustion and mixing in ing process (Fig. 2. IO) will be carried out in ground test supersonic flow with the aim of burning length reduc- complex and in flying test laboratories. Ground test tion and combustion efficiency of hydrogen increase in complex is the most prepared in spite of it is necessary wide range of Mach numbers; to modernize it and equip it by laser methods of measurements of reacting gas parameters. This complex - the control of structure thermal state in Iiigh- includes of stationary supersonic wind tunnels, enthalpy supersonic flow by use of active cooling system hypersonic facilities of short-time action, impulse tun- operating heat flux peaks moving along the scramjet nels, ballistic facilities and special equipment for pro- duct surface; duction of high-enthalpy air flows with extremely clean composition with the aid of fire, plasma and electric - the use of hydrogen slush as coolant and fuel; heaters. - the use of compact hydrogen-air heat The problems of diffusion and detonation exchangers; combustion in supersonic flow, of control of boundary layer and wall flow at presence of shock waves and heat - the use of light-weight structures of compos- flux peaks can be considered in ground conditions at ites for controlled scramjet duct. MCI and global Reynolds number Re-107. Air chemical composition and initial turbulence levels do not 6. Scramjet systems and elements (Fig. 2.13). practically correspond by this to conditions at flight in Unique systems and elements should be made on the atmosphere almost always. The use of flying laborato- base of new technical approaches, materials, technology ries is for this reason of extremely importance. Wing diagnostic and control systems in the process of flying laboratories of ballistic flying laboratories, scramjet developing. It is possible to point out for ex- launched by ballistic rocket can realize design flight ample the followings: path of aerospace plane in atmosphere H=H(AI” both on separate parts and in a whole and to ensure the - supersonic mixers on the base of gas-genera- condition of scramjet vehicle interference without sim- tors like liquid propellant engines for thrust creation ulating real Reynolds numbers due to small dimensions and subsequent effective mixing of supersonic air flow of model scramjet duct. The orbital laboratories with with supersonic fuel jet in wide ranges of Mach num- planning descent and following return to orbit or land- bers and fuel components ratios; ing can be useful for studies of final period of scramjet operation at Mf> 16. - multiniode supersonic combustor, permitting 9-4 to realize conditions of subsonic and supersonic com- path. bustion in wide range of Mach numbers and fuel coni- ponents ratios; Flight tests of subscale models permit to in- vestigate operating process at conditions made more - active cooling system with system of diag- realistic as for duct in whole and for its elements along nostic and control of coolant injection, ensuring struc- the flight paths close to real flight path. ture thermal state at conditions of changing heat flux with high nonuniforniity. Mathematical simulation of flows including physical-chemical conversions has as objectives the Scramjet conception, that will be developed investigations of scramjet performances in whole range jointly with aerospace plane conception, can be realized of external and internal parameters, the investigation - - step by step in the form of modules, at first, in the form various parameters influence on efficiency, economy of subscale module and, then in form of full-scale and specific weight of scranijet, and the optimization of module, and in the end-in form of full scale scramjet as duct parameters for given goal function. a main engne of flying aerospace plane propulsion system (Fig. 2.14). Full-scale flight test of modules and propulsion In the process of scranijet development, its system on the whole permit to investigate all parameters elements and modules should be subjected to various of scranijet along possible flight paths. tests with the aim of providing design parameters along the flight path, reliability and service life (Fig. 2.15). 2. The advantages of various methods of in- The necessary industrial base should have available vestigations.(Fig. 2.18) e specialized test facilities for development of each element and system for module development and for full Experimental investigations in laboratory per- scale scranijet tests. Let's cite as example specialized mit to' study non stationary and non equilibrium test facility for conibustor development both in respect gasdynamic phenomena in detail including viscosity to thermal state and in respect to combustion efficiency and turbulence and mechanism of flow stabilization at and impulse. In particular, combustor impulse its separation and energy release diffusion and detona- measurements, which is connected with appreciable tion combustion in range Mc=2... 5 in core on the flow error at test facility conditions. Let's cite as the second and in wall regon. These studies permit to obtain ex- example supply system of liquid or slush hydrogen, peririiental information with the aid of optics-electronic which is connected with system of active cooling of methods of measurements. duct. The main restrictions of ground tests are connected with restricted ranges of Mach number Operating process simulation in test facility (M' <: 8) and Reynolds number (Re .-: 107) and permits to simulate phenomena on comprehensive in incomplete simulation of interference of scramjet with limited test duration. Flight tests of sub-scale models vehicle. It is necessary to use the set of methods of permit to simulate comprehensive the phenomenon for scramjet study and development to eliminate these real air composition along the whole flight path of real disadvantages: operating process simulation in the scranijet in limited range of Re and at limited test du- laboratory and in test facility, flying tests of models, ration. Flight test permit to simulate propulsion systeni- mathematical simulation and full-scale flight tests. vehicle interference and to use independent methods of measurements of thrust balance and drag, and scramjet This point is illustrated by comparison of var- specific impulse. ious methods of scramjet duct operating process in Fig. 2.16. The objects of the investigation, merits and Mathematical simulation of flows including demerits of various methods and their interaction are physical-chemical conversions permits to calculate all formulated in presented Fig. 2.17 ...2.20. parameters of problem in limitless ranges, gives neces- sary for design detailed data on flow patterns both on 1. Operating process investigation objects (Fig. the whole and in the small, permits to generalize ex- 2.17). Tests in laboratory pursue the objects of revealing perimental result obtained in laboratory, test facility and and elaborating of physical flow pattern and some flight tests and to simulate interference of propulsion characteristics of phenomena in scramjet duct elements system and plane. (inlet, combustor, exhaust nozzle) with possible siniu- lation of physical-chemical conversions in high Full-scale flight tests of modules and propul- enthalpy flow with talung into account influence of high sion system on the whole permit to simulate all real cooled walls. conditions of scranijet operation.
Simulation of operating process in test facility 3. The disadvantages of various methods of permits to investigate it at close to real conditions both investigations (Fig. 2.19) for duct in whole and for its elements in possible range a of M and Re corresponding to different parts of flight Unfortunately, each of mentioned methods 9-5
doesn't exhibit necessary completeness for comprehen- The verification and extrapolation of test facilities re- a sive investigation of scramjet operating process. Exper- sults by methods of sub-scale flight tests are also nec- imental studies in laboratory don't comprehensive sim- essary. ulation of operating process both in scramjet duct on the whole and in its elements due to impossibility to keep Flight tests of sub-scale models should use all similarity criteria, detailed duct configuration and results of experimental laboratory studies, operating initial and boundary conditions of flow. They permit to process simulation in test facilities, and mathematical simulate only some characteristic features of flows and simulation of flows. physicochemical phenomena and to obtain heterogeneous qualitative results. The verification and extrapolation of obtained results are necessary a full-size flight test. Operating process simulation in test facility is possible only in limited range of hyand Re (enthalpy Mathematical simulation of flows including and absolute flow velocity) and for limited model di- physical-chemical conversions permits to use the whole mensions. This leads to only partial simulation of date base, obtained in experiments in laboratory, in test physical-chemical composition of air, initial and facility, at sub-scale model's flight tests and at full size boundary conditions of flow in model scramjet duct. flight tests. All previous results obtained by all methods Limited nature of information obtained by measure- of study are also used at pre starting procedure of full- ments in test facility, together with mentioned above scale flight tests and at their result's analysis. disadvantages leads to necessity to use a set of facilities 0 for parameter's simulation along the flight path. Finally, it should be noted that scramjet design and development are the key problem for transorbital air Small dimensions of scramjet model duct and transport. The solution of this problem will result in limited natures of measurements methods are the dis- creation of hypersonic passenger into orbit enough advantages of flight tests of sub-scale models. space industrial opening up.
Disadvantages of mathematical simulation of In accordance with high complexity and a large flows are connected with use of semi-empirical models body of fundamental and applied problems, scramjet of turbulence and physical-chemical phenomena and problem can be classified as one of the most difficult with impossibility of exact integration of motion equa- scientific-technical problems. Scramjet development is tions for real numbers A4 and Re and gas composition. connected with parallel and interacting development of f its conception, hypersonic technology and industrial test Full-size flight tests of modules and propulsion base. system on the whole are possible only for limited number of tests. The disadvantage of this kind of tests in The realization of scramjet duct operating pro- the limited information about operating process and cess should be carried out on the base of all available extremely high cost and high risk of tests. scientific and technical means including ground, flight and orbital laboratories with comprehensive mathemat- 4. Interaction of investigation methods (Fig. ical simulation. 2.20) As to experiment in laboratory, the comparison of results obtained by different methods is necessary for The solution of considered problems, due to estimations of their reliability and measurements' errors their importance, complexity, cost, and risk, are possible and to obtain reference data, verification and ex- and expedient on the base of international collaboration trapolation of experimental results by use of all others and coordination of efforts of scientists, engineers and (mentioned above) methods are necessary. governments.
The use of results of experimental tests in lab- oratory and mathematical simulation of flows in neces- sary at operating process simulation in tests facilities. 9-6
jointly. The purpose of the second flight zst was 3 Subscale Model SCRAMJET Flight Tests define more exactly engine inside parameters at subsonic and supersonic combustion modes in during Ground test facility can not realize all HFL flight along the trajectory with hgher condition of scramjet operations. The limitations of air ram values. Stable operation of scramjet on both subsonic enthalpy, air chemical composition breaking, or supersonic combustion modes was in during more turbulence conditions breaking and some other are than 23 seconds of the second flight test. typical features of hypersonic test facility. These defects can be corrected practically by flight tests only. Hypersonic Laboratory and But real size aerospace plane scramjet flight tests are Flying Evperimental Double Mode,SCRAMJET. extremely expensive and possible as final stage of system creation. But many necessitate data may be To provide the HFL flight trajectory, close to obtain due to flight testing of subscale model scramjet. a standard flight trajectory of the vehicles with the propulsion system incorporating the scramjet, was one Subscale model scramjet flight test can not of the main causes due to which "Surface-to-Air" to use real size intake to save condition of boundary missile SA-5 was chosen the "Kholod" layer transition and limited sizes of nozzle don't as HFL booster. This missile flight trajectory is very close to permit to obtain good thrust efficiency of scramjet. All the required trajectory of the HFL "Kholod". another conditions may be very close to real size flight Proceedings from the requirements of aerodynamics, vehicle and propulsion. stability and control ability of SA-5, experimental scramjet and all on-board service modules were The first scramjet flight test could be X-15 manufactured in a form of axisymnietrical circular scramjet flight test. But plane wreck disturbed this body; their diameter did not exceed 750 mm diameter project. of the standard SA-5 modules (Fig. 2.21). Some preparation of scramjet flight test was Experimental hydrogen scramjet of axisym- beginning in Soviet Union by head company CIAM as metrical configuration installed in the in 70-ths beginning. The preliminary engineering by was HFL forebody. The casing of the section N1 is an M.M. Bondarjuk and E.S. Schetinkov have effect on intermediate section between the experimental flight tests preparing. The works were conducted scramjet and on-board tank, occupying the section N2. under leadership of R.I. Kursiner, D. A. Ogorodnikov and V.A. Sosounov. The section N 1 is washed from the outside by high temperature jet emerging from the scramjet Our model scramjet flight tests were planed annular nozzle and for this reason it is protected by 12 as part of scramjet researches. A small surface-to-air mm thermal-insulating layer. Hydrogen consumption rocket is suggested as booster. Ground tests of sinular controller (RPST-2) with actuating units and scramjet must be carried out before flight tests and regulating flaps, ignition system electronics, flow many kinds of theoretical and computational meter devices and transducers to measure the investigated were planed too. parameters in the experimental scramjet were installed within the section N 1. Hydrogen supply The axisymmetrical double mode scramjet lines, measuring tubes for pressure transducers, developed and created by Russia aviation was thermocouples and electric cables are also installed in specialists for flight tests specially. Due to strong the above section. financial limit the first unburning scramjet flights were in end of 80-ths. Flight tests of the Hypersonic On-board hydrogen tank (section N 2) has Flying Laboratory with hydrogen scramjet (HFL) shield-vacuum heat insulation with residual pressure firing have been conducted on November 28 1991 on lesser than 0.1 mm Hg. The operational pressure firing ground in Kazachstan for the first time. The within the internal vessel at hydrogen supply to the preliminary analysis has shown that experimental scramjet is 2,2 MPa. Hydrogen level monitoring scramjet had started in flight and was operating with during on-board tank filling is carried out by fuel injecting corresponding t'o subsonic and capacitance level transducer and by internal vessel supersonic combustion modes. wall temperature sensor. Due to high altitude and not enough ram Expulsion system of hydrogen supply from value stable operating duration of supersonic the on-board tank to the scramjet combustor is combustion mode was several seconds only at the first installed in section 3A. Hydrogen expulsion system flight test of scramjet. incorporates high pressure helium spherical container (37 ma), the volume of which is 42 liters, pressure The second flight test was in November 17, reducers, isolating and safety valves, and sensors for 1992. The works on the second flight test were parameters monitoring in the hydrogen supply system. conducted by Russian and French aviation specialists 9-7
High pressure nitrogen spherical container (37 MPa) combustor. At supersonic combustion mode (M$ 0 is also located in section 3A. The volume of this 5.0). part of a fuel is additionally supplied through the container is 17 liters. Nitrogen is used for functioning first collector with combustion in the first part of the of pneumatic actuated valves and hydrogen supply combustor. The injectors are located uniformity along system units in flight. Joining units of on-board tank the circle and their number is the same in all filling by hydrogen and spherical containers by helium collectors (N=42 holes). The injector nozzle diameter and nitrogen are placed in the same section; is 1,7 nini for the first row of injection, the injector electropneumatic joints to fulfill the technological and diameter is 2.1 mm for the second and the third rows control operation on maintenance position and on of injection. Hydrogen is injected into the flow at launching pad during integrated checks of on-board angle 30 degreases in the first and the second row of - systems and hydrogen filling is also installed in injection and it is injected at angle 90 degreases in the section 3A. Telemetry system units, electric power third row of injection. batteries, automatic pilot, yning indicator of scramjet combustor wall hazardous temperature and ' The walls of both the central body and the the ignition equipment of SA-5 missile control system external ring are cooled by the hydrogen used for are located in the section 3B. combustion. At the exit of the tank, hydrogen first flows along the combustor double walls and then is I-IFL measuring complex together with RTS sent to the injection valve supplying the last two (radio telemetry station) makes available recording collectors and to that supplying the first injectors' row. signals of boost rocket systems transducers, pressure If necessary, the third valve permits to increase the a transduceh, temperature transducers and indicators of hydrogen flow rate in the cooling jacket and to scramjet and HFL equipment operation from the start exhaust the extra hydrogen out of the combustion time up to 144th s of flight. HFL was equipment by chamber. 83 transducers of pressure, vibration and overload detectors, 58 temperature sensors and 46 indicators of The chart of fuel supply systeni and basic equipment operation. Engine and its elements' principles of fuel consumption control systeni design parameters were measured by 68 pressure transducers, are presented in Fig, 2.22. As operating regmes 49 thermocouples and temperature sensors and by 20 change should be ensured only by fuel consumption indicators of equipment operation. Information from variation and by fuel redistribution over fuel injector the measuring complex was transmitted to telemetry rings, basic requirements to computer-aided control station with 50 Hz sampling rate (from temperature systeni of dual mode scramjet are concerned only sensors it was transmitted with 1,5 Hz sampling rate) with fuel consumption control and change of fuel through on-board telemetry station. External trajectory injection position in accordance with flight Mach measurements were fulfilled by radar and phase number. The requirements of providing necessary direction finder complex. cooling of engine structure elements and engme stable operation at maximum heat supply regmes are also The experimental dual mode scramjet of important. axisymmetric configuration consists of the uncooled three-shock inlet, the combined cooled combustor and As it is shown in Fig. 2.22a, engine internal the exhaust nozzle with one-sided expansion. The parameters, i.e. static pressures PI, P4, P5, PHIand PI.I, design of the mentioned engine was carried out as four and total pressure Po determines fuel consumption and basic units: intake central body, combustor central temperatures of structure elements (there can be body, lips of combustor and outside shell (Fig. 2.22). several elements) determine coolant consuniption through cooling jacket. By this, inlet first cone The engine is'designed to operate in a wide pressure determines incoming flow Mach number range of flight conditions (Mf =3,5 ...6 and H=15 ...35 km) at subsonic and supersonic combustion modes experimental scramjet is made with constant duct dimensions to simpllfy engme structure and to ensure pressures PHI and Pm determine real more high reliability. performances of conibustor operating process,
The first part presents a small divergence and includes at its beginning a row of injectors followed by a flame stabilization. The second one, which is and pressures P4 and P5 characterize flow divergent, comprises at its entrance a step and two regime at duct entrance. If P4 P5, cavity. The last one has a constant area. For subsonic then flow with detached flow shock wave and decrease combustion mode, (M~3.5... 5,O) the fuel is supplied air flow mass rate occurs. through the injectors of the collectors with fuel burnout in the second and in the third part of the 9-8
scramjet is corroborated by overload factor Nx Incoming Flow Parameters HFL measurements. The difficulties of such integral a measurements are caused by rocket engine operation The main parameters of the incoming flow, till 57 s of flight and after its termination by small necessary for analysis of engine operation in flight, -- value of dual mode scramjet thrust in comparison with are the following: pressure PH , temperature TH and air drag force of the whole complex HFL-dual mode density p14, well as the angle of attack of the as scramjet. But overload sensor has recorded dual inode incoming flow. External trajectory measurements and scramjet thrust in the initial period of test at ~40...44 special meteorological measurements:are used to s (Fig. 2.25). determine these parameters (H,W)
Experimental Dual Mode SCRAMJET Parameters. . Mach number was calculated by the typical relation M=Wd(kgRTdl/2 The change of M, H and Gas-air flow parameters engme duct and ram value q,, depending on HFL flight time are scramjet inlet and combustor performances were presented at Fig. 2.23. Presented data show that the determined by use of measured values of static maximal flight Mach number value was M =5.3 and pressure and wall temperature. incoming flow the maximal altitude was H z 22.4kni. parameters, air mass flow rate, fuel and coolant consumption's. The calculations were carried out on SCRAMJET Operation in Flight one-dimensional gas dynamic equations including real Conditions. thermodynamic properties of combustion products, drag forces of struts and cavity stabilizers. friction and The changes of combustor wall temperature heat removal into duct walls. Heat fluxes qw near fuel injection positions and hydrogen temperature accumulated by wall substance was calculated with measured at cooling system exit are shown in Fig. use of unsteady heat conduction method. In 2.24. The presented data testify to active hydrogen accordance with this method, qw is proportional to combustion in dual mode scramjet combustor. The dTw /dT (qMW/dT). Heat amount removed by analysis has shown that the highest wall temperature coolant in cooling system channels' cowl and central at subsonic combustion regime T~1140KOwas body were calculated with the assumption that coolant recorded on cowl near 2 cavity stabilizer real wall. consumption on cowl constitutes 80% of coolant The maxinial values of Tw at supersonic combustion consumption on central body. regimes Tw =I400 ... 1450 KO were recorded on central body of the first part of conibustor near cavity The typical data on changes of flow pressure stabilizer rear wall at ~70...72 s. Measured value of increase ratio P=P/PH,wall temperature Tw, and total hydrogen temperature at cooling system exit was heat fluxes Qsw along experimental dual mode maximal in the end of dual mode scramjet operation scramjet duct at subsonic and supersonic combustion duration and constituted Tm=950 KO(Fig. 2.24). regimes are presented in Fig. 2.26. One can see from considering dependencies P(X) that pseudo-shock, in We can see from considering Tw and T, which air stream deceleration from supersonic velocity dependence on time that thermal state of combustor to subsonic one, accompanied by sharp pressure structure elements at engne operation at subsonic increase (Fig. 2.26a) occurs, lies at subsonic 0 combustion regimes close to steady state after was combustion regime in initial part of combustor. coolant consumption increase at 49...52 s. But after Disturbances caused by heat release in 2 and the transition to supersonic combustion regimes with second parts of combustor don't extend upstream in additional heat release in the first part of combustor, the first part of combustor farther than till 1 cavity the thermal instability increased, gradients dTw/dr stabilizer. At the third part of conibustor with constant and dTm/dr increased too, that resulted in dual mode cross section area, pressure decreases due to intensive scramjet construction break-down. The complete heat supply and reaches minimum value at combustor analysis of experimental data, taking into account exit. Unusual for such regimes' non monotonous conibustor walls thermal state, combustion efficiency change of pressure P in zone X=6.5 ...7.5 can be and pressure losses in cooling system channels, has explained apparently by combined effect on the flow shown that the first symptoms of combustor walls' of fuel jets system of the second and the third injection burnout have appeared already at ~63...64 s. The rows, cavity stabilizers, active heat supply and change appearance of regimes with detached shock wave throat are of the second part of combustor. subsequently was caused evidently by this. For this reason, the analysis of dual mode scramjet main The analysis has shown that pseudo-shock parameters was made for HFL flight time from 39 till moves downstream at flight speed increase 63 s. accompanied by decrease of relative heat supply, and static pressure in the second and the third parts of Effective hydrogen combustion in combustor combustor increases. The maximum pressure increase and producing thrust by experimental dual mode ratio in engine duct is reached values P = 30... 32 at 9-9
subsonic combustion regmes. scramjet inlet parameters determined in flight test e conditions are shown in Fig. 2.27. The first flight test The flow field in engine duct changes data are presented in the same figure for comparison. fundamentally at supersonic conibustion regimes with Reynolds numbers calculated on inconiing flow additional hydrogen supply through the first injection parameters and cowl sharp edge diameter D,=O.226 ni row in combustor begnning (see Fig. 2.26b). The varied within the limits Re=(0.5-2.1)106 and RH1.4- region with maximum pressure is located in the first 2.7)106 in the first and the second flight tests, part of conibustor and the pressure in the second and respectively. One can see a good agreement of inlet the third parts of combustor decrease by a factor parameters (pressure increase ratio Pi, Mach number 1.5... 2. The maximum pressure increase ratio in duct Mi and total pressure recovery factor 0;)obtained in is at this operation mode appreciably higher than that both tests are in good agreement. at subsonic combustion regimes and reaches values PP@5 and greater e-63 s. Some distinction between Mi and oi values, obtained in different tests, appreciable at supersonic P(X) considering in other flight moment has combustion regmes MH >5 is caused, apparently, by shown that the flow in the second and the third heat release in the first part of combustor effects on combustor parts during the initial period of engine boundary layer and separation zone in inlet throat. It operation at supersonic combustion mode is close to is possible that pressure increase on central body isobaric one. The flow field, characteristic for compression surfaces is caused by the same reason, at subsonic combustion mode , with flow deceleration least near pressure P5 points of measuring. and subsequent flow acceleration, is reached after 58- th second of flight at flight Mach number decrease in Flow Mach number M, total pressure losses tlie third part of combustor (XB7.8). Therefore, it is and hydrogen fuel combustion efficiency distributions possible to suppose that both completely supersonic along the dual mode scramjet combustor on flow along the duct and that combined with subsonic characteristic subsonic and supersonic conibustion zone in the third part of conibustor are realized at fuel regmes are illustrated in Fig. 2.28. Flow Mach supply corresponding to supersonic combustion number change in the first part of combustor at regimes. subsonic combustor regime (Fig. 2.28a) characterizes flow deceleration in pseudo-shock wave zone. The It should be noted, that the pressure in the mentioned above peculiarity of pressure distribution first part of combustor increases appreciably at engine along the length at combustor second part (see Fig. operation with detached shock wave (P = 72) and 2.26a) becomes apparent in flow acceleration up to disturbances spread upstream and give rise to average transonic velocity in section x=7.0 with two strongly pressure increase more than 2 times. pronounced subsonic flow zones at the second and the third part of combustor. Main losses of total pressure Tw(X) distributions, shown in Fig. 2.26 gve at these regimes are concentrated in the first part of additional insight into thermal loads of dual mode combustor and in the beginning of the second part in scramjet inlet and conibustor as compared with data of the regon of air flow deceleration. Combustion Fig. 2.24. One can see those heat fluxes on inlet efficiency, for an air to fuel ratio of 1.2 reaches 0.8. compression surfaces (X<4) are small and values of The most active fuel jets burning occur till ~4.7,in wall temperature Tw of central body increase during the second and in the beginning of the third part of time interval ~43.5...57.5 s not higher than by ATw combustor (Fig. 2.28a). =50 ... 100 KO. The comparison of Tw(X) and Qsw(X) shows that specific heat fluxes in the first part of The results presented in Fig. 2.28b example combustor increase appreciably at feed of fuel part of operation regme for which supersonic flow is through the first row of injection. This fact together achieved all along the engine duct. Here the main with distributions P(X) change testifies to hydrogen total losses occur in the first part of the chamber, combustion along the whole engine duct at subsonic where combustion takes place in a supersonic flow combustion regimes at 7 >52 s. with high initial velocity (about Mach 2.2) and also in the second and third flame stabilizers zone. The total As the analysis has shown, static pressure pressure decrease is about 50% in both regions. It does values P measured on inlet compression surfaces are not exceed 30% or 40% in the rest of the combustor. close to design values and experimental data obtained in ground test conditions. For this reason, calculated Non monotonous variations of Mach number experimental dependencies of inlet coefficient cp and M in the first part of combustor is explained by liquid streamline additional drag force factor Cx upon combined effect of heat release and cross section area incoming flow Mach number and angle of attack were increase. It is obvious, that the effect of heat release is used to determine engine air mass flow rate and flow prevailing in stabilizer's region due fuel burns in tlie parameters in inlet throat (section X4.2, before the circulation zone and heat release is maximum at the first injection ring). Experimental dual mode beginning of hydrogen jets combustion. 9-10
The analysis results also show that the total pressure recovery ratio cr&,,,b is constant behind the stabilizer cavity at that combustion efficiency, for the first injection; q1, is about 0.7 at the end of the first combustor part. This can lead to the conclusion that Some Conclusion. probably, for small fuel mass flow rate combustion mainly takes place in the flame stabilizer operates not During the second flight test of the duel in flame regime. The global combustion efficiency at mode scramjet in the Hypersonic flight Laboratory the combustor exit is about 0.7--0.8. "Kholod", the engine operated in the range of flight Mach nuniber running from Me3.5 to Me5.35 and Data scatter at combustion efficiency factor between altitude of 15 km and 24 km is subsonic and variation along the combustor in both examples supersonic combustion operating modes. Hydrogen presented in Fig. 2.28 is explained by approximates of combustion in a flow that was supersonic all along the one-dimensional approach to the analysis of scramjet duct was achieved for the first time in flight experimental data on pseudeshock zone and high test conditions, between Mach number of 5.1 and 5.3 flow nonuniformity near fuel jets burnout 'begmning, and with air to fuel ratio of about 1.8. including repon near the second and the second cavity stabilizer. But the results of comparison of gas For subsonic combustion regime. with air to dynamic method used in present paper with thrust fuel ratio running from 1.2 to 2.5, combustion method of q determination results, obtained at efficiency varied from 0.65 to 0.88. In supersonic experimental investigations of scramjet bed models , combustion mode, and for air to fuel ratio comprised e have corroborated the authenticity of determining between 1.7 and 2, it varied from 0.7 to 0.95. averaged flow parameters in region with moderate nonuniformity extent: at combined combustor exit and The cooling system was able to ensure in the end of its first part. normal scramjet operation without structure break- down during 23 s. Hydrogen supplying system and The generalization of the first and second engine regulation could also achieve engine starting flight tests experimental data related to combustion and it's in two different combustion modes. efficiency is presented in Fig. 2.29. The synibols indicated experimental data and the solid lines the Thus, results have demonstrated technical approximation of these data. The leading parts of feasibility of flight tests with rocket boosters and their diffusion processes in hydrogen combustion are ability to provide significant information about corroborated by typical shape of correlation's between scramjet operation. combustion efficiency and equivalence ratio p. with efficiency decrease near p = 1 and increase when p greater or low then 1. The higher values of efficiency in supersonic combustion regimes are explained by increase of fuel injectors' number and by the larger combustor length over which fuel jets mixing with air and conibustion occurs. 9-11
Biblionrafhv 10. Kuczera H., et.al., "The German Hypersonics Technology Programme - Status 1993 and Perspectives", 1.-Kyp31i~ep P.M. "PeaKTkiBHble ABHlXTeJIM IAF-93-V.4.629, 1993. &llH 6onbum CBepX3ByKOBbIX CKOPOCTefi IlOJleTa" M. Mauunocmpoenue, 1989. 11. Koelle D., "Launch Cost Assessment of Operational Winged Launch Vehicles", AIAA Paper 92- 2. Ogorodnikov D., "Englnes for Hypersonic 5021, 1992. Aircraft", '2viadvigatelestroenie-92" International Symposium, 1992 12. Vandenkerckhove J., Barrere M., "Energy Managment", ISABE - 1993. .3. Contensou P., Marguet R., Huet Ch., "Etude Theorique et Experimentale dun Statoreacteur a 13. Billig F., "Research on Supersonic Combustion mixte (domaine de vol Mach 3.5/7)" La Combustion," Jornal oJPropulsion and Power,.. . .. Vol. 9, _RecherchAerospatiale - No. 5, 1973. No. 4, 1993.
4. Hirsinger F. "Optimisation des Performances 14. Heppenheimer. "The national Aerospace dun Statoreacteur Supersonique - Etude Theorique et plane" Pasha Publications''s 1989. Experimentale", TP ONEM No. 1 106, 1992. 15. Escher W.J.D., Teeter R.R., Rice E.E. 5. Vinogradov V., et.al. "Experimental "hrbreathing and Rocket Propulsion Synergism: envestigation of 2-D Dual mode Scramjet with Hydrigen Enabling Measures for Tomorrods Orbital Transports" Fuel at Mach 4-6",AM Paper 90-5269,1990. AIAA Paper, 86-1680, 1986
6. h6ero~O., Kyp3"ep P., M ap., 16. Murthy S.N.B., Czysz P.A. "Energy Analysis "HeKOTOpbIe OCO6eHHOCTM pa6osero npouecca B of High speed Vehicle Systems" AIAA Paper, 90- KaMepe CJOPaHMR C A03BYKOBblM M . 0089,1990. CBepX3ByKOBblM IlOTOKdMM", "Teopm ABMrdTeJIefi JleTaTeJlbHblX aIlnapaTOB", AH CCCP, 1986.
7. Roudakov A., Novelli Ph., et.al., "Flight testing an hisymmetric scramjet - Russian Recent Advances." IAF- 93-S.4.485, 1993
8. Sosounov V. "Study of Propulsion for High Velocity Flight" ISABE, 1991.
9. Kleinau W., Koelle D., "Evolutionary Develoopment to a Reusable Single-Stage-to-Orbit encher", IAF-93-V. 3.6 19, 1993. FLIGHT PATH REGION, ENTHALPY CONSTANT LINES AND RAM CONSTANT LINES
H,km
60
50
40
30
20
10
0 5 10 15 20 25 Mf Fig. 2. I
LINES OF CONSTANT QUASI-EQUILIBRIUM TOTAL TEMPERATURE Tt, TOTAL PRESSURE P, AT SCRAMJET ENTRANCE AND THE FLIGHT PATH REGION
H, Kn
60
50
40
30
20
10
0 5 10 15 20 25 Fig. 2.2 9-13
LINES OF CONSTANT REYNOLDS NUMBER RE1, AND TEMPERATURE FACTOR Tw=~w/~,AT SCRAMJET ENTRANCE AND THE FLIGHT PATH REGION
H, urn
60
50
40
30
20
10
0 5 10 15 20 25 Mf Fig. 2.3
INLET THROAT MACH NUMBER M,, AND COMBUSTOR EXIT MACH NUMBER Mm DEI'ENDING OEJ FLIGHT MACH NUMBER M,
Fig. 2.4 9-14
INFLUENCES of VARIOUS LOSSES ON SCRAMJET SPECIFIC IMPULSE ALONQ THE FLIGHT PATH
3
25
20 2
10
5
5 10 15 20 ME Fig. 2.5
Scheme of Research and Design Work for Scramiet Development
/{
\
Fig. 2.6 9-15
Hypersonic TechnolopIy Development Directions
Hypersonic technology I
Experimental New technical Fundamental studies simulation of engine approaches process
Applied scientific Numerical simulation studies of engine of engine process elements process 1 Fig. 2.7
Fundamental Studies --Physical process and turbulence. --Materials Applied Scientific Studies of Engine Process --Ecology. I Experimental Simulation of Engine Process
Numerical Simulation of Engine Process
New Technical Approaches
Scramjet Systems and Elements
Fig. 2.8 9-16
Fundamental Studies --Thermodynamics.
--Thermal physics and combustion. Applied Scientific Studies of Engine Process --Gas dynamics and turbulence. Experimental Simulation of \ --Numerical simulation. Engine Process --Diagnostics and measurements Numerical Simulation of methods. Engine Process --Super computers. New Technical Approaches
Scramjet Systems and Elements
Fig. 2.9
Fundamental Studies --Ground test facility: --gas dynamics; Applied Scientific Studies of Engine Process --combustion and detonation;
Experimental Simulation of --heat exchange. Engine Process --Flying laboratories:
Numerical Simulation of --ballis tic; Engine Process --winged 3New Technical Approaches --Orbital lOFig. 2.10 9-17
Fundamental Studies --Mathematical models of turbulence in hypersonic high- enthalpy flow. Applied Scientific Studies of Engine Process --Nonqulibrium physical-chemical reactions in supersonic flow. Experimental Simulation of Engine Process --Three-dimensional and unsteady interaction of shock waves with Numerical Simulation of / boundary layer in channel. Engine Process --Detonation combustion. --Direct numerical simulation of New Technical Approaches scramjet flow. 0 --Scramjet thermal state Scramjet Systems and simulation. Elements --Operating process control.
Fig. 2.11
Fundamental Studies --Structure. --Materials. Applied Scientific Studies of a Engine Process --Technology. Experimental Simulation of --Efficiency. Engine Process --Diagnostics and control. Numerical Simulation of Engine Process --Supersonic mixers on base of liquid rocket engine. New Technical Approaches i’ --Multimode supersonic combustor.
Scramjet Systems and --Automatic control system. Elements e Fig. 2.12
9-19
Test Facilities Development
Test facilities
I
Development of scramjet Scramjet module Real-size scramjet service units and systems life tests -- Combustor: --Demonstration and 3 < Mf <8 1. Thermal state service life tests, 3 < Mf 2. Combustion efficiency <8 3. Thrust impulse -- Scramjet module and --Supply system Dlane interference
Fig. 2. I5
Methods of Scramiet Duct Operating; Process Study
Mathematical Simulation of Flows Including PhysicalChemical Conversions
Test Facility Simulation of Operating Process
Flight Tests of Sub--Scale Models
Experimental Studies in Laboratory
Real--Size Flight Tests of Modules and Propulsion System on The Whole
Fig. 2.16 9-20
Scramiet Duct Operating Process Study Methods
Experimental Studies in Test Facility Flight Tests of Sub- Mathematical Real--Size Flight Laboratory Simulation of Scale Models Simulation of Flows Tests of Modules Operating Process Including and Propulsion PhysicalChemical Conversions ivestigation Objects --Revealing and --Investigation at --Investigation at --Investigation of --Investigation of elaboration of physical conditions close to conditions close to scramjet all scramjet duct flow pattern and some real both for duct on real both for duct on performances on the parameters along characteristics of the whole and for its the whole and for its whole range of possible flight phenomena in scramjet elements in possible elements along flight external parameters. paths. duct elements (inlet, range of M and Re paths close-to real --Investigation of combustor, exhaust numbers, various parameters nozzle) with possible corresponding to mfluence on scramjet I simulation of different parts of efficiency, economy physicalchemical flight path and specific mass. conversions in high-- --Optimization of enthalpy flow including duct parameters for mfluence of high cooled given goal function. walls I
Fig. 2.17
Scramiet Duct Operating Process Study Methods
Experimental Studies in Test Facility Flight Tests of Sub--Scale Mathematical Real--Size Laboratory Simulation of Models Simulation of Flows Flight Tests of Operating Process Including Modules and PhysicalChemical Propulsion Conversions System on The
--The opportunity of detailed --Comprehensive --Comprehensive simulation --Unlimited range of Simulation of nvestigation of unsteady and simulation of of phenomena for real air parameters. dl real ion equilibrium gasdynamics phenomena in composition along the whole --Detailed data on flow Zonditions of ihenomena including flow limited ranges of flight path of real scramjet in pattern on the whole and ;cramjet nscosity and turbulence and M and Re limited range of Re number in the small. Jperation. nechanics of flow stabilization numbers and at and limited test duration --Generalization of ind energy release limited test --Simulation of propulsion experimental results, --The same for diffusion and duration system plane interference obtained in laboratory, letonation combustion process --Independent methods of test facility and flight n range M = 2...6 in flow core measurements of scramjet test. ind wall region thrust and drag balance and --Simulation of --The grater information specific impulse interference of Ibtained by optic--electronics propulsion system and neasurements plane.
Fig. 2.18 9-21
Scramiet Duct OperatinP Process Study Methods
Experimental Studies in Test Facility Simulation of Flight Tests of Mathematical Real--Size Laboratory Operating Process Sub--Scale Simulation of Flows Flight Tests of Models Including Modules and PhysicalChemical Propulsion Conversions System on The Whole iethods disadva tages --Incomplete simulation of --Limited ranges of M and --Small --Use of semi-- --Limited operating process both in Re numbers(entha1py and dimensions of empirical models of number of tests scramjet duct on the whole absolute flow velocity), model turbulence and --Limited and in its elements due to ' model dimensions, scramjet duct. physical--chemical mformation on impossibility to keep all incomplete simulation of Limited nature phenomena in operating similarity criteria, detailed physical--chemical of connection with process. duct configuration, initial compositions of flow in measurements impossibility of exact --Extremely and boundary conditions of model scramjet duct. methods. integration of motion high cost and flow. --Limited nature of equations for real Re risk of tests. --Simulation of only some mformation obtained by and M and gas characteristic features of measurements. composition. flows and phenomena. --Necessity to use a set of --Heterogeneous test facilities for simulation qualitative results. of parameters along the
flight path. I
Fig. 2.19
Scramiet Duct Operating Process Study Methods
Experimental Stydies in Test Facility Simulation Flight Tests of Sub-- Mathematical Real--Size Laboratory of Operating Process Scale Models Simulation of Flight Tests 01 Flows Including Modules and PhysicalChemical Propulsion Conversions System on The Whole :thods Interaction --Comparison of different --Results of methods 1 --Results of methods --The whole --Results of laboratory methods is and 4 are used. 1,2,4 are used. data base methods 1... 4 necessary to estimate --A set of --Verification and obtained by are used. measurements reliability and measurements methods extrapolation of methods 1,2,3,5 errors and to obtain reference is necessary to obtain results by methods is used, especially data. reference results. 2,4,5 are necessary. data of methods --Venfication and --Verification and 1,2. extrapolation of laboratory extrapolation of results methods results by methods by methods 3,4,5 are 2...5 is necessarv. necessarv.
Fig. 2.20 \-EXPERIMENTAL DMSCRAMJET
Fig. 2.21. Hypersonic flying laboratory with the experimental dual node scramjet on start position 9-23
collision b tlow rate +-,7 t enpera t ore
pbPI P4 Ps Pri fro1 tank dtnpinq W REQUIRMEYTS TO TEE EXPEEIHENTAL DHSCRAHJET CONTROL SYSTEH:
@ FUEL FLOW RATE CRANCINC AND DISTRIBUTION ON THE INJECTOR PONS Iltlltlll) IN ACCORDANCE WlTB FLICBT HACB NUHBDR
0 STABLE ENCIWE OPERATION FOR TEE STEADY AND UNSTEADY HODES (STARTI HAXIHUH BPAT ADDITION RODESl IWLXT rLON STALLI ac toators COOLIYC STRUCIURK** * 1
8 WIDK RAWCL 01 TBO OPERATIOWAL HODES BY PUP1 lLOW RATE AND FLlCIIT MCB YUHBOR FOR SUBSONIC AND SUPPRSONIC COHBUSTiON .FUEL FLOW e CBANCINC OY CONTROL PROMS AND LIHITS RATE e SUFPlClEYT ACCURACY, SIHPLICITY, RELIABILITY...
Fig. 2.22. Conception of the experimental dmscramjet control system a - functional scheme of the control system; b - scheme of the fuel feed and cooling system
1I HE, S
Fig. 2.23. HFL traiectorv rmm eters (external traiectorv measureme nts EM - (light Mach number; EH:- altitude (km); OH - dynamic pressure (kg/m3) Fig. 2.24. Hal 1 temperatures (THl2, TH20, TU4 1 and hydrogen temperature at cooling system exit (TB53) in the dmscramjet operational the + - TW12, X - TW20, o - TW24, D - TB53
NX
TIC, 5
Fig. 2.25 Axial overload factor NI t = 40 s beginning of engine operation 9-25
- a - b e D P
30 L(0 20 20 10
-. . Y JW .h TU .K 1000
8 00
6 00
-. YO0 ...
30
Fig. 2.26 Chanpe of dmscraajet parameters by envlne length (P - pressure ratlo. TW - uall temperature. QSH - heat flux: llvht symbols - coullnp. dark symbols - central body, - repulator measurements) a - subsonlc conbustlon, TRlJz43.5 S. M~z4.32. H.17.4 kn. RLFR=I.II; b - ;z;ersonlc combiistlon. TflU:57.5 s. "15.33. t!:2!.', h.11!-?9:!.3?
a H b
-... .8 .6 I .o .v .s .I .o .o
Flou paraneters In the dmscramlet conbustlon chamber tH - Mach number. SIGMCHIIHB Fig. 2.28 - a total Pressure recovery coefflclent. ETR - combustlon efficiency. modes as FIp.8) - Pi
15
3 4 5 M 10
I.._ 7 4 5 H 5 3 4 5 H
r meters in fli ndition Fig' 2'27 Pi st~~~~~~~~,"~~~t~~~achnumber inlet thsoat - ai total pressure recovery coefficient first flight + second flight
62 tr 1 I
0 0.25 0.5 0.75 i.0 1.25 1.50 1.75 2.0 Hydrogen combustlon efficiency coefficients at the experlmental scramjet Fig. 2.29 @ - fuel feedln for SU!ISO~~Ccombustion node through I1 and 111 injector rous: ob 0 - ice1 feedinc for supersonic. corchustim mode through I. I1 and I11 lnlecbor rous;#op- start chamber Bodes or near "shock out" modes 1. Recipient’s Reference 2. Originator’s Reference 3. Further Reference 4. Security Classification of Document AGARD-LS-194 ISBN 92-835-0732-0 UNCLASSIFIED/ UNLIMITED
6. Title RESEARCH AND DEVELOPMENT OF WSCRAMJETSAND TURBOWETSIN RUSSIA
8. Author(s)/Editor(s) 9. Date Various December 1993
LO. Author’s/Editor’s Address 1 1. Pages Various 232
I 12* DistributionStatement .This document is distributed in accordance with AGARD policies and regulations, which are outlined on the back covers of all AGARD publications. 13. Keywords/Descriptors Ramjet engines Computational fluid dynamics Solid rocket propellants Russia Liquid rocket propellants Flight tests
14. Abstract
Russia has a long tradition of achievement in ramjet research and development. This tradition, and aspirations toward new and effective products, have led to establish Russian priority in the ramjet field.
This Lecture Series will present and discuss the scientific problems of the development of ramjets/ scramjets and turboramjets.
Some specific aspects of liquid/solid ramjet development, the concepts of LH2 high efficiency RAM combustors, the results of full scale turboramjet testing, scramjet or CFD analyses and their ground flight tests will be studied.
This Lecture Series, endorsed by the Propulsion and Energetics Panel of AGARD, has been implemented by the Consultant and Exchange programme. t; E2 W ';J 00 :w 4 w p;' 0
t; E2 W ';J ca w
z4 w p;' 0 a
cn
e
cn
e CANADA LUXEMBOURO du tascimtifiqug voir&lglqoe Orman,Ontodo K1A OK2 NORVEGE DAhEMAUK N~DcfcnmRaaaFchEstlMishmsat DpnishDdcacc~~tab~t Acm: Bibliotetot Ryvmtg~AIUl P.O. Box 25 P.0.Bm 2715 N-2007 Kjdh DK-2100Cbpdap10 PAYS-BAS RSPAGNE NcmcrLndalM~@i~mAGARD JNTA AGARDPUblicptionr) L.borptory- phnorL34 """"'7P.O. Box 9050 28008 Madrid 1006BMAmstdmn ETATS-UNIS h~cF37,Diatributioac5lm 300 E Street, SW. Webinnton,D.C. 20546 BN92-835-07326