Experiments were performed in order to examine the stability of hydrocarbon-fueled flames in cavity flameholders in supersonic airflows. Methane and ethylene were burned in two different cavity configurations having aft walls ramped at 22.5 and 90$^{\circ}$. Air stagnation temperatures were 590 K at Mach 2 and 640 K at Mach 3. Lean blowout limits showed dependence on the air mass flowrates. Visual observations, planar laser induced fluorescence (PLIF) of nitric oxide (NO), and Schlieren imaging were used to investigate these phenomena. Large differences were noted between cavity floor and cavity ramp injection schemes. Cavity ramp injection provided better performance in most cases. Ethylene pilots have a wider range of stable operation than methane. Fuel flowrates at ignition showed similar trends as lean blowout limits, but higher flowrates were required.

This paper outlines a new convergence acceleration de-signed to solve scramjet flowfields with zones of re-circulation. Named the “marching-window”, the algorithm consists of performing pseudo-time iterations on a minimal width subdomain composed of a sequence of cross-stream planes of nodes. The upstream boundary of the subdomain is positioned such that all nodes upstream exhibit a residual smaller than the user-specified convergence threshold. The advancement of the downstream boundary follows the advancement of the upstream boundary, except in zones of significant streamwise ellipticity where a streamwise ellipticity sensor ensures its continuous progress. Compared to the standard pseudo-time marching approach, the march-ing-window is here seen to decrease the work required for convergence by up to 24 times for supersonic flows with little streamwise ellipticity and by up to 8 times for supersonic flows with large streamwise separated regions. The memory requirements are observed to be reduced sixfold by not allocating memory to the nodes not included in the computational subdomain. The marching-window satisfies the same convergence criterion as the standard pseudo-time stepping methods, hence resulting in the same converged solution within the tolerance of the user-specified convergence threshold. The extension of the marching-window to the weakly-ionized Navier-Stokes equations is also discussed.

In this study, a boundary-layer bleeding and a two-staged fuel injection were applied to a scramjet engine for suppressing unstart transition and improving the thrust performance under Mach 6 flight conditions. With the boundary-layer bleeding, the engine could operate without unstart transition around at the fuel equivalence ratio of unity ($\Phi$ = 1). The thrust increment from the no fuel condition (dF) increased to 2460 N, which was about 1.4 times as large as that of the case without the bleeding and maximum in our Mach 6 tests. It was confirmed that the boundary-layer bleeding suppressed the separation during the engine operation. The two-staged fuel injection was less effective for improving the thrust performance com-pared with the single-staged one with the bleeding at Mach 6.

Two different type scramjet models with side-wall compression and top-wall compression inlets have been tested in HPTF (Hypersonic Propulsion Test Facility) under the experimental conditions of Mach number 5.8, total temperature 1700K, total pressure 4.5㎫ and mass flow rate 3.5kg/s. The liquid kerosene was used as main fuel for the scramjets. In order to get fast ignition in the combustor, a small amount of hydrogen was used as a pilot. A strut with alternative tail was employed for increasing the compression ratio and for mixing enhancement in the side-wall compression case. Recessed cavities were used as a flameholder for combustion stability. The combustion efficiency was estimated by one dimensional theory. The uniformity of the facility nozzle flow was verified by a scanning pitot rake. The experimental results showed that the kerosene fuel was successfully ignited and stable combustion was achieved for both scramjet models. However the thrusts were still less than the model drags due to the low combustion efficiencies.

At the University of Queensland (UQ), research into the performance of high-speed (in excess of 7000km/hr) air-breathing engines in the form of supersonic combustion ramjets (or scramjets) has been made for almost 20 years. This has been possible because the T3 and T4 shock tunnels, located at the ANU and UQ, respectively can simulate these conditions. However, like all facilities, there are differences between the flow generated in these facilities and that, which occurs in flight. The correlation between the two has not been determined for these facilities, or indeed for any shock tunnel performing supersonic combustion experiments. The aim of the HyShot flight program is to obtain this correlation by undertaking a sounding rocket program based at Woomera in South Australia. The seminar will discuss new approach taken by the UQ researches in developing this cost effective flight program, as well as the triumphs and disappointments which have been encountered so far in completing this somewhat ambitious program.

Air-fuel mixing and flame-holding are two important factors that have to be considered in the design of an injection system. Different injection strategies have been proposed with particular concern for rapid air-fuel mixing and flame-holding. Two representative injection techniques can be applied in a supersonic combustor. One of the simplest approaches is a transverse(normal) injection. The cavity flame holder, an integrated fuel injection/flame-holding approach, has been proposed as a new concept for flame holding and air-fuel mixing in a supersonic combustor. This paper describes numerical efforts to characterize the flame-holding and air-fuel mixing process of a model scramjet engine combustor, where hydrogen is injected into a supersonic cross flow and a cavity. The combustion phenomena in a model scramjet engine, which has been experimentally studied at University of Queensland and Australian National University using a free-piston shock tunnel, were observed around the separation region of the transverse injector upstream and the inside cavity. The results show that this flow separation generates recirculation regions which increase air-fuel mixing. Self-ignition occurs in the separation-freestream and cavity-fteestream interfaces.

The effects of injector spacing s and injector diameter d on mixing are numerically investigated in supersonic combustor with perpendicular injection behind a backward-facing step. Simulations are reported for airstream Mach number of 2.4. Parameters are changed on following 4 cases to investigate the effects of injector configuration on mixing efficiency $\eta_m$. In the case of varying d or s, dynamic pressure ratio $Rq(=(pu^2)_j/(pu^2)_a)$ is also varied to keep bulk equivalence ratio $\Phi({\oe})Rq.d^2/s)$ constant. (l) Injector spacing s is varied at constant $\Phi$=0.5, 1, 2 for injector diameter d=6mm. In the case of $\Phi$=1, $\eta_m$ has its maximum value at s=24mm. The reason is that increase of $\eta_m$. , by widening spacing at Rq=constant competes with decrease of $\eta_m$ by increasing Rq at s=constant. When spacing is narrow, the flow field of vicinity of injector becomes two-dimensional because adjacent jets interferes each other. By widening spacing, air is easily entrained by three-dimensional effect. This mechanism also appears in the case of $\Phi$=0.5, 2 for d=6mm, and $\eta_m$. reaches its maximum value at s=24mm for $\Phi$=0.5 and at s=42mm for $\Phi$=2. (2) In the case of injector diameter d varied at $\Phi$=1 for s=30mm, $\eta_m$. has its maximum value at d=3mm. The reason is that decrease of $\eta_m$ by increasing injector diameter competes with increase of $\eta_m$ by decreasing Rq at d=constant.(3) In the case of s varied at $\Phi$=0.5, 1,2 for d=3mm, the injector spacing at which mixing efficiency has its maximum value is s= 18mm for $\Phi$=0.5, s=24mm for $\Phi$=1, s=24mm for $\Phi$=2. Therefore it is found that d=3mm and s=24mm can be optimum configuration over a range of $\Phi$=0.5~2.(4) The effect of h on the optimum spacing is investigated. s is varied for d=6mm at step height h=4, 6, 8mm. The simulation results do not show significant change on the step height.

Numerical simulations were conducted in a rectangular scramjet combustor with a cavity and/or a step in order to investigate their performances for flame-holding. Flow structures and OH radical profiles in the cavity and the step were calculated. The calculated results showed that the cavity generated a larger recirculation zone than the step that had the same depth. Additionally, the combustor with a cavity could make a large low-velocity area than the combustor with a step. The cavity performance was determined by its depth and length. The cavities with too large or too short length did not work effectively, and a certain aspect ratio showed high performance for flame-holding. There was a minimal depth under which the cavity did not work as flame-holder. The fuel injections upstream the cavity and inside the cavity were also tested to investigate the effects on the cavity performance. The result showed that the fuel injection inside the cavity reduced reaction areas and residence time. Therefore, the upstream injection was preferable to the inside injection.

A comprehensive numerical analysis has been carried out for both non-reacting and reacting flows in a scramjet engine combustor with and without a cavity. The theoretical formulation treats the complete conservation equations of chemically reacting flows with finite-rate chemistry of hydrogen-air. Turbulence closure is achieved by means of a k-$\omega$ two-equation model. The governing equations are discretized using a MUSCL-type TVD scheme, and temporally integrated by a second-order accurate implicit scheme. Transverse injection of hydrogen is considered over a broad range of injection pressure. The corresponding equivalence ratio of the overall fuel/air mixture ranges from 0.167 to 0.50. The work features detailed resolution of the flow and flame dynamics in the combustor, which was not typically available in most of the previous studies. In particular, the oscillatory flow characteristics are captured at a scale sufficient to identify the .underlying physical mechanisms. Much of the flow unsteadiness is related not only to the cavity, but also to the intrinsic unsteadiness in the flow-field. The interactions between the unsteady flow and flame evolution may cause a large excursion of flow oscillation. The roles of the cavity, injection pressure, and heat release in determining the flow dynamics are examined systematically.

In scramjet engines with sidewall compression inlet, it is well known that a non-uniform flow appears since a separated region is generated near the flow centerline on the body side. The separated region is caused by shock-boundary layer interaction and likely to cause un-start phenomena since the flow in the separated region is subsonic and acts as a communication path between the isolator and the combustor. In the present study, the non-uniform flow characteristics in the scramjet inlet-isolator region are numerically studied in detail. Effect of flow suction from body sidewall surface on the non-uniform flow field numerically examined to clarify the flow mechanism to suppress the un-start transition.

Various jet engines (Turbine engine family and RAM Jet engine) have been developed for high speed aircrafts. but their application to hypersonic flight is restricted by principle problems such as increase of total pressure loss and thermal stress. Therefore, the development of next generation propulsion system for hypersonic aircraft is a very important subject in the aerospace engineering field, SCRAM Jet engine based on a key technology, Supersonic Combustion. is supposed as the best choice for the hypersonic flight. Since Supersonic Combustion requires both rapid ignition and stable flame holding within supersonic air stream, much attention have to be given on the mixing state between air stream and fuel flow. However. the wider diffusion of fuel is expected with less total pressure loss in the supersonic air stream. So. in this study the direction of fuel injection is inclined 30 degree to downstream and the total pressure of jet is controlled for lower penetration height than thickness of boundary layer. Under these flow configuration both streams, fuel and supersonic air stream, would not mix enough. To spread fuel wider into supersonic air an aerodynamic force, baroclinic torque, is adopted. Baroclinic torque is generated by a spatial misalignment between pressure gradient (shock wave plane) and density gradient (mixing layer). A wedge is installed in downstream of injector orifice to induce an oblique shock. The schlieren optical visualization from side transparent wall and the total pressure measurement at exit cross section of combustor estimate how mixing is enhanced by the incidence of shock wave into supersonic boundary layer composed by fuel and air. In this study non-combustionable helium gas is injected with total pressure 0.66㎫ instead of flammable fuel to clarify mixing process. Mach number 1.8. total pressure O.5㎫, total temperature 288K are set up for supersonic air stream.

This paper report our preliminary results of characterizing the jet structures of kerosene injection into quiescent atmosphere and a Mach 2.5 crossflow at various preheat temperature. A heating system has been designed and tested that can prepare heated kerosene of 0.8 kg up to 670 K at a pressure of 5.5 ㎫. Temperature measurement near the injector shows that the temperature of pressurized kerosene can be kept constant during the experimental duration. Comparison of kerosene jet structures in the preheat temperature range of 290-550 K demonstrates that with injection pressure of 4 ㎫ the jet plume turns into vapor phase completely at injection temperature of 550 K, while keeping the penetration depth essentially unchanged. The results suggest that the injection of vaporized fuel would improve the performance of a liquid hydrocarbon-fueled supersonic combustor because the evaporation process is now omitted.

As a research to develop a SCRAM jet engine is actively conducted, a necessity to produce a high-enthalpy flow in a laboratory is increasing. In order to develop the SCRAM-jet engine, stabilized combustion in a supersonic flow-field should be attained, in which a duration time of flow is extremely short. Therefore, a mixing process of breathed air and fuel, which is injected into supersonic flow-fields is one of the most important problem. Since, the flow inside SCRAM jet engine has high-enthalpy, an experimental facility is required to produce such high-enthalpy flow-field. In this study, a detonation-driven shock tunnel was built and was used to produce high-enthalpy flow. Further-more, SCRAM jet engine model equipped backward-facing step was installed at test section and flow-fields were visualized using color-schlieren technique and high speed video camera. The fuel was injected perpendicular to the flow of Mach number three behind backward-facing step. The height of the step, distance of injection and injection pressure were changed to investigate the effects of step on a mixing characteristic between air and fuel. The schlieren photograph and pressure histories show that the fuel was ignited behind the step.

An experimental research on supersonic combustion of kerosene in a model scramjet combustor has been conducted. Kerosene was injected normally into a Mach 2 vitiated airstream either at an atomized liquid state or at a gaseous state. The atomization of kerosene was achieved by the “effervescent atomization” method, and the gaseous kerosene was supplied by passing kerosene inside a heated pipe. The results are discussed and are also compared to those in our previous experiment, in which no atomization nor vaporization methods has been conducted to the kerosene.

The effects of pressure on combustion characteristics of the conical flameholder were investigated experimentally in the pressure range from 0.11 ㎫ to 0.40 ㎫. The result shows that the total equivalence ratio of lean limit becomes lower as the combustor inlet pressure rises and NOx emission is proportional to the pressure to the 0.5th power.

Generally, flow around a bluff body such as a circular cylinder is complicated compared with that around a streamlined body because of the existence of separated shear layers. Long bluff body such as a flat blunt plate is more complicated than short bluff body, because of separated-and-reattaching flow on the after bodies.(omitted)

An experiment was carried out to confirm the validity of time series evaluation of supersonic mixing condition by using catalytic reaction on a platinum wire. Geseous hydrogen was injected parallel to supersonic freestream (M$\infty$ $\approx$ 1.81) from a slit injector, which was located at backward facing step. Time series condition of supersonic mixing was evaluated by using W-type probe which has a platinum wire and reference wire (nickel wire). The evaluation was by simultaneously measuring each electric circuit which kept the temperature of wire constant. Investigations were also conducted for helium, air and no secondary injectant cases to compare with the hydrogen injectant case. The results indicated that it was possible to measure the time series behavior of air and hydrogen supersonic mixing layer or coherent motion of turbulence by using this evaluation.

The mixing enhancement is one of the most important problems for the development of scramjet engines. The influence of the streamwise vortices produced by a ramp in a unheated supersonic flow on the mixing of twin jets injected from its base was experimentally investigated. Nominal Mach number of the main airstream and of the twin jets at the nozzle exits were 2.35 and 2.0, respectively. Three dimensional velocity distributions near the ramp with and without injection were measured by Particle Image Velocimetry (PIV). A pair of counter rotating streamwise vortices could be seen behind the injector without injection. On the other hand, two pairs of streamwise vortices could be seen with injection. The outer one had the same direction as the vortex pair produced by the ramp, but they were stronger than those produced by the ramp. The inner ones had the opposite directions to the outer ones. It is considered that these vortices enhance the mixing near the injector.

Experimental study on combustion characteristics of double swirl coaxial injectors has been conducted for the assessment of critical injector design parameters. A reusable, unielement thrust chamber has been fabricated with a water-cooled copper nozzle. Two principle design parameters, a swirl angle and a recess length, have been investigated through hot firing tests for the understanding of their effects on high pressure combustion. Clearly, both parameters considerably affect the combustion efficiency, dynamics and hydraulic characteristics of an injector. Internal mixing of propellants in a recess region increases combustion efficiency along with the increase of a pressure drop required for flowing the same amount of mass flow rates. It is concluded that pressure buildup due to flame can be released by the increase of LOx flow axial momentum or the reduction of a recess length. Dynamic pressure measurements of the thrust chamber show varied dynamic behaviors depending on injector configurations.

The ongoing developmental studies on the application of hydrogen peroxide for propulsion are briefly reviewed. A detailed design-study of a laboratory scale facility of a hydrogen peroxide mono-propellant engine of 100-N thrust is presented. For the preparation of concentrated hydrogen peroxide, a distillation facility has been realized. Results of water analogy tests are presented. Initial firings using the concentrated hydrogen peroxide were not successful. Low environmental temperature, low contact area of the catalyst pack, and contamination in the hydrogen peroxide were considered to be the reasons. Addressing the first two points resulted in successful firing of the rocket engine.

Many monopropellant thrusters use a catalyst for decompose the propellant, hydrazine. The catalyst directly affects the thruster performances and lifetime. Therefore, it is important to confirm that the catalyst is suitable for our thrusters. Until 2002, we used She1l405 catalyst, for satellite RCS thrusters, and H-IIA and M-V launch vehicle upper-stage RCS thrusters. In 2002, however, Shell Chemical Inc. ceased manufacturing She1l405 catalyst and transferred the product to AEROJET, where it was renamed S405. We subsequently investigated the characteristics of AEROJET's S405 catalyst and SOLVAY's KC12GA catalyst, (SOLVAY is a Belgian chemical company, and KC12GA is used for ASTRIUM's thruster) and conducted thruster firing tests using the new catalysts. After conducting, we confirm that the KC12GA catalyst was suitable for our thrusters, and decided to use KC12GA for two satellite programs.

The flow in the LOX manifold of liquid rocket has been investigated using a CAE technique with an objective of economical modeling of injection holes in order to reduce the overall computational cost of flow analysis during the optimal rocket design procedure. The computational geometry is very close to that of the actual rocket design and the flow condition through the injection holes resembles that in the actual manifold of the liquid rocket. The result shows that the flow in the plane just above the injection holes is not uniformly distributed in terms of pressure and mass flow rate and this is attributed to the large-scale flow patterns present the LOX manifold. Thus, the flow physics should be understood correctly before making any attempt to model the injection holes. In the present study, several boundary conditions which were designed to effectively replace the presence of injection holes have been tested and it was found that a simple modeling can be possible by mimicking the actual geometry of the injection holes. By using this simple injection hole modeling, it was able to obtain about 30% reduction in computational cost but it was still able to reproduce the flow patterns correctly. Also the flow has been analyzed after incorporating a couple of different types of pre-distributors in LOX manifold and the effect of those will be discussed.

As a rocket propellent of hydrocarbon fuels, the characteristics of liquefied natural gas was evaluated with the viewpoint of the constituents and content, the cooling performance as a coolant, and characteristic velocity and specific impulse as parameters of the engine performance. Content of methane was a principal factor to determine the characteristics as a rocket propellant and more than 90 % of it was needed as a fuel and coolant in the regenerative cooled liquid rocket engine. Some constituents of the liquefied natural gas can be frozen by the pre-cooling of the pipe lines, therefore they can be a factor disturbing the normal working of engine. In case the content of methane is around 90% in the liquefied natural gas, a normalized stoichiometric O/F mixture ratio of 0.75 is suggested for a nominal operation condition to get the maximum specific impulse and characteristic velocity.

In order for studying pressure-coupled dynamic responses of droplet vaporization, open-loop experiment of an isolated droplet vaporization exposed to pressure perturbations in stagnant gaseous environment is numerically conducted, Governing equations are solved for flow parameters at gas and liquid phases separately and thermodynamic parameters at the interfacial boundary are matched for problem closure. For high-pressure effects, vapor-liquid interfacial thermodynamics is rigorously treated. A series of parametric calculations in terms of mean pressure level and wave frequencies are carried out employing a n-pentane droplet in stagnant gaseous nitrogen. Results show that wave instability in view of pressure-coupled vaporization response seems more susceptible at higher pressures and higher wave frequencies. Mass evaporation rate responding to pressure waves is amplified with increase in pressure due to substantial reduction in latent heat of vaporization. Augmentation of perturbation frequency also enhances amplification due to the reduction of phase differences between pressure perturbation and surface temperature fluctuation.

LE-7A is the main engine of the H-IIA launch vehicle. Under its development, the nozzle suffered from two troubles during startup and shutdown transients of the engine. One is a large side load, which damages the actuator of the nozzle, and the other is damage on regenerative cooling tubes due to high heat load. It has been considered that these problems are caused by a peculiar separation pattern called Restricted Shock Separation (RSS). RSS is observed in several rocket nozzles, for example, LE-7A nozzle, Vulcain nozzle and so on. Their contours are not conventional truncated perfect (TP) nozzle - LE-7A nozzle is a compressed truncated perfect (CTP) nozzle and Vulcain nozzle is a thrust optimized (TO) nozzle. Although it is believed that the occurrence of RSS is affected by the nozzle contour, the mechanisms are not clarified sufficiently yet. In the present paper, a parametric numerical study is carried out to investigate the mechanisms of the occurrence of RSS in CTP nozzles during startup transient. The results show that RSS is caused by the adverse pressure gradient downstream of the Mach disk. The adverse pressure gradient is caused by the interaction of the pressure wave and Mach disk. The method to avoid the occurrence of RSS is also examined. A small step inside the nozzle affects the position of the separation point and prevents RSS. The result shows that the possibility that RSS can be suppressed by controlling the position of the separation point.

Breakup characteristics of liquid sheets formed by the impingement of two water jets, such as a breakup length and a breakup wavelength of sheet, were investigated as increasing the injection velocity up to 30m/s and the ambient gas pressure up to 4.0㎫. While round edged orifices formed a laminar sheet which has no waves on the sheet when the injection velocity is low, sharp edged orifices formed a turbulent sheet which has impact waves irrespective of the injection velocity. Thus we compared the differences of breakup characteristics between them. The results showed that the aerodynamic force significantly affects the breakup of laminar sheet when the gas based Weber number is higher than unity. It was also found that the turbulent sheets have three breakup regimes, i.e. expansion regime, wave breakup regime and catastrophic breakup regime according to the gas based Weber number.

Our previous study showed that although the hybrid rocket engine with swirling gaseous oxygen had high performance, a direct injection of LOX with swirl into the combustion chamber of the hybrid rocket engine lowered the performance of the engine, compared to that with gaseous oxygen. In order to clarify this reason, combustion tests of a small PMMA combustor with an inner port diameter of 2 mm were conducted in liquid oxygen flow by comparison with gaseous oxygen flow. Although the oxygen mass fluxes of LOX were about two orders of magnitude larger than those of gaseous oxygen, the fuel regression rate of LOX were remarkably smaller than those of gaseous oxygen. For both liquid and gaseous oxygen, diffusion flames in the port of the grain controlled the combustion process of PMMA in oxygen flow. These results may be explained by the fact that only small amount of LOX vaporized and consumed in the combustion with PMMA while flowing through the port due to relatively larger latent heat of injected liquid oxygen compared to the heat of release by combustion which depended on the burning surface area of PMMA.

The reactive flowfield of the transverse injecting combustor has been studied using Euler-Lagrange method in order to develop an efficient solution procedure for the understanding of liquid spray combustion in the transverse injecting combustor which has been widely used in ramjets and turbojet afterburners. The unsteady two-dimensional gas-phase equations have been represented in Eulerian coordinates and the liquid-phase equations have been formulated in Lagrangian coordinates. The gas-phase equations based on the conservation of mass, momentum, and energy have been supplemented by combustion. The vaporization model takes into account the transient effects associated with the droplet heating and the liquid-phase internal circulation. The droplet trajectories have been determined by the integration of the Lagrangian equation in the flow field obtained from the separate calculation without considering the iterative effect between liquid and gas phases. The reported droplet trajectories had been found to deviate from the initial conical path toward the flow direction in the very end of its lifetime when the droplet size had become small due to evaporation. The integration scheme has been based on the TEACH algorithm for gas-phase equation, the second order Runge-Kutta method for liquid-phase equations and the linear interpolation between the two coordinate systems. The calculation results has shown that the characteristics of the droplet penetration and recirculation have been strongly influenced by the interaction between gas and liquid phases in such a way that most of the vaporization process has been confined to the wake region of the injector, thereby improving the flame stabilization properties of the flowfield.

Recently, fluidic thrust vectoring methods have been preferably employed to control the movement of propulsive systems due to relatively simpler design and lower cost than mechanical thrust vectoring methods. For An application of the thrust vectoring to flight bodies, it is necessary to understand very complicated exhaust flows which are often subject to shock waves and boundary layer separation. But researches for the thrust vector control using counterflow have been few. In the present study, experiments have been performed to investigate the characteristics of supersonic jets controlled by a thrust vectoring method using counterflow. The primary jet is expanded through a two-dimensional primary nozzle shrouded by collars, and is deflected by the suction of the air near nozzle into an upper slot placed between the primary nozzle and the upper collar. A shadowgraph method is used to visualize the supersonic jet flowfields. Primary nozzle pressure ratios and suction nozzle pressure ratios are varied from 3.0 to 5.0, and from 0.2 to 1.0 respectively. The present experimental results showed that, for a given primary nozzle pressure ratio, a decrease in the suction nozzle pressure ratio produced an increased thrust vector angle. As the suction nozzle pressure ratios were increased and decreased, the hysteresis of the thrust vectoring was observed through the wall pressure distributions

액화천연가스를 연료로 사용하여 물 냉각 및 천연가스와 액화천연가스 재생냉각 연소시험을 수행하였다. 연소시험과 CEC86을 이용한 연소해석 결과를 액체로켓엔진 성능인자로서, 특성속도와 비추력 관점에서, 추진제 혼합비와 연료의 연소실 유입온도의 영향을 분석하였으며, 엔진성능이 추진제 혼합비와 연료의 연소실 유입온도의 영향을 크게 받고 있음을 알 수 있었다. 엔진 성능으로서 특성속도는 추진제 혼합비가 0.72∼0.75일 때, 이론적 비추력은 추진제 혼합비가 0.75일 때 최대 값을 보여주었으며, 연료의 연소실 유입온도의 증가에 비례하여 엔진 성능이 향상되는 경향에서 재생냉각이 엔진 성능을 증대시키는 경향을 확인하였다.

The pulse detonation engine (PDE) has recently expected as a new aerospace propulsion system. The PDE system has high thermal efficiency because of its constant-volume combustion and its simple tube structure. We measured thrust of single-tube pulse detonation rocket (PDR) by two methods using the PDR-Engineering Model (full scale model) for ground testing. The first involved measuring the displacement of the PDR-EM by laser displacement meter, and the second involved measuring the time-averaged thrust by combining a load cell and a spring-damper system. From these two measurements, we obtained 130.1 N of time-averaged thrust, which corresponds to 321.2 sec of effective specific impulse (ISP). As well, we measured the heat flux in the wall of PDE tubes. The heat flux was approximately 400 ㎾/$m^2$. We constructed the PDR-Flight Mode] (PDR-FM). In the vertical flight test in a laboratory, the PDR-FM was flying and keeping its altitude almost constant during 0.3 sec.

For design of cryogenic propellant feeding system, one of the main requirements is to meet temperature requirement for satisfying turbo-pump NPSH requirement. In this paper improved method of estimating the thermal stratification in liquid oxygen tank is presented to help design. In the case of liquid rocket using turbo-pump, the inner pressure of liquid oxygen tank is maintained low, so vaporization of liquid oxygen is generally occurred. In this paper, inner process of LOX tank is analyzed by two phase flow modeling. The vaporization rate and required helium mass is investigated.

A new cavitating model by using bubble size distribution based on bubbles-mass has been proposed. Both liquid and vapor phases are treated with Eulerian framework as a mixture containing minute cavitating bubbles. In addition vapor phase consists of various sizes of vapor bubbles, which are distributed to classes based on their mass. The bubble number-density for each class was solved by considering the change of the bubble-mass due to phase change as well as generation of new bubbles due to heterogeneous nucleation. In this method, the bubble-mass is treated as an independent variable, and the other dependent variables are solved in spatial coordinates and bubble-mass coordinate. Firstly, we employed this method to calculate bubble nucleation and growth in stationary super-heated liquid nitrogen, and bubble collapse in stationary sub-cooled one. In the case of bubble growth in super-heated liquid, bubble number-density of the smallest class based on its mass is increased due to the nucleation. These new bubbles grow with time, and the bubbles shift to larger class. Therefore void fraction of each class is increased due to the growth in the whole class. On the other hand, in the case of bubble collapse in sub-cooled liquid, the existing bubbles are contracted, and then they shift to smaller class. It finally becomes extinct at the smallest one. Secondly, the present method is applied to a cavitating flow around NACA00l5 foil. Liquid nitrogen and liquid oxygen are employed as working fluids. Cavitation number, $\sigma$, is fixed at 0.15, inlet velocities are changed at 5, 10, 20 and 50m/s. Inlet temperatures are 90K in case of liquid nitrogen, and 90K and 1l0K in case of liquid oxygen. 110K of oxygen is corresponding to the 90K of nitrogen because of the same relative temperature to the critical one, $T_{r}$=$T/T_c^{+}$. Cavitating flow around the NACA0015 foils was properly analyzed by using bubble size distribution. Finally, the method is applied to a cavitating flow in an inducer of the LE-7A hydrogen turbo-pump. This inducer has 3 spiral foils. However, for simplicity, 2D calculation was carried out in an unrolled channel at 0.9R cross-section. The channel moves against the fluid at a peripheral velocity corresponding to the inducer revolutions. Total inlet pressure, $Pt_{in}$, is set at l00KPa, because cavitation is not generated at a design point, $Pt_{in}$=260KPa. The bubbles occur upstream of the foils and collapse between them. Cavitating flow in the inducer was successfully predicted by using the bubble size distribution.

The spray and acoustic characteristics by the self-oscillation of a swirl coaxial injector were experimentally studied. The self-oscillation of a swirl coaxial injector is defined as pressure and flowrate oscillations by a time-delayed feedback between liquid and gas phase and has strong influences on atomization and mixing processes. Hence the occurrence and effect of the self-oscillation are measured using shadow photography technique, acoustic test and PDPA. The occurrence of self-oscillation largely depends on the injection conditions, such as pressure drop of liquid phase and relative momentum ratio. From the experimental results, self-oscillation occurs when the momentum of gas phase is enough large and the smaller the pressure drop of liquid phase is, the better self-oscillation occurs at the same momentum ratio. The self-oscillation is also affected by injector geometries, increasing the recess length results in the expansion of self-oscillation region and the increase of sound pressure level. The self-oscillation of a swirl coaxial injector accompanies a high intensity scream and this scream may provide harmful disturbances to combustion processes. Self-oscillation leads to strong changes in the drop size distribution and smoothly varies the slope of radial SMD distribution.

In the development of rocket turbo-pump, cavitation at the inducer is one of the major problems. Cryogenic fluids are commonly used for rocket propellant, therefore, thermodynamic effect becomes noticeable compared to conventional water cavitation. In the present study, a numerical simulation method for cryogenic cavitation is proposed, which reveals the difference between cryogenic and water cavitation.

Characteristics of steady and unsteady cavitation in a turbopump inducer were investigated in this paper. To see the effect of tip clearance on the inducer performance, three cases of tip clearance were tested. The helical inducer, which has two blades with inlet tip blade angle of 7.8 degree and tip solidity of 2.7, was tested in the water. In the non-cavitating condition, the inducer head decreased with increase in the tip clearance. Rotating cavitation and cavitation surge were observed through unsteady pressure measurements at the inducer inlet. The cell number and propagation speed of the rotating cavitation were determined through cross-correlation analysis. During the rotating cavitation one cell rotated at the same rotational speed as that of the inducer rotation and the cavitation surge did not rotate. The critical cavitation number increased with increase in the tip clearance at the same flow rate, but the change of critical cavitation number was small at the nominal flow rate.

The discharge current oscillation has been measured for various hollow anode widths and its axial positions using a 1㎾-class anode layer hall thruster. As a result, there were thresholds of magnetic flux density for stable discharge. The plasma structure inside the hollow anode was numerically analyzed using the fully kinetic 2D3V Particle-in-Cell (PIC) and Direct Simulation Monte Carlo (DSMC) methods. The results reproduced both stable and unstable operation modes. In the stable operation case, which corresponds to the case with low magnetic flux, the plasma penetrated into the hollow anode deeper than the case with higher magnetic flux density case. This suggests that comparably large substantial anode area should contribute to stable operation.

A concept in which laser-sustained plasmas (LSPs) are combined with inductively coupled plasmas (ICPs) is proposed. The concept is aiming at extensions of operative conditions of a CW laser thruster due to the fact that the ICP has some characteristics which are in contrast to those of LSPs. An estimation confirmed that the concept would effectively work. And a fundamental experiment was conducted. The results showed that the radio frequency magnetic field induced by a alternate current of 13.56 MHz coupled inductively with LSPs, resulting in the enlargement of the plasma region and the attainment of the enthalpy. It is expected that some improvements will enable to transfer the RF power to the work gas more effectively and to demonstrate the synergy effect between the LSPs and the ICPs.

Computational fluid dynamics analysis was carried out for thermo-chemical flow field in Arcjet thruster with mono-propellant Hydrazine ($N_2$H$_4$) as a working fluid. The theoretical formulation is based on the Reynolds Averaged Navier-Stokes equations for compressible flows with thermal radiation. The electric potential field governed by Maxwell equation is loosely coupled with the fluid dynamics equations through the Ohm heating and Lorentz force. Chemical reactions were assumed being infinitely fast due to the high temperature field inside the arcjet thruster. An equilibrium chemistry module for nitrogen-hydrogen mixture and a thermal radiation module for optically thin media were incorporated with the fluid dynamics code. Thermo-physical process inside the arcjet thruster was understood from the flow field results and the performance prediction shows that the thrust force is increased by amount of 3 times with 0.6KW arc heating.

Basic experiments were carried out using the THT-IV low-power Hall thruster to examine the influences of magnetic field shape and strength, and acceleration channel length on thruster performance and to establish guidelines for design of high-performance Hall thrusters. Thrusts were measured with varying magnetic field and channel structure. Exhaust plasma diagnostic measurement was also made to evaluate plume divergent angles and voltage utilization efficiencies. Ion current spatial profiles were measured with a Faraday cup, and ion energy distribution functions were estimated from data with a retarding potential analyzer. The thruster was stably operated with a highest performance under an optimum acceleration channel length of 20 mm and an optimum magnetic field with a maximum strength of about 150 Gauss near the channel exit and with some shape considering ion acceleration directions. Accordingly, an optimum magnetic field and channel structure is considered to exist under an operational condition, related to inner physical phenomena of plasma production, ion acceleration and exhaust plasma feature. A new Hall thruster was designed with basic research data of the THT-IV thruster. With the thruster with many considerations, long stable operations were achieved. In all experiments at 200-400 V with 1.5-3 mg/s, the thrust and the specific impulse ranged from 15 to 70 mN and from 1100 to 2300 see, respectively, in a low electric power range of 300~1300 W. The thrust efficiency reached 55 %. Hence, a large map of the thruster performance was successfully made. The thermal characteristics were also examined with data of both measured and calculated temperatures in the thruster body. Thermally safe conditions were achieved with all input powers.

As one of the concepts of the laser/electric hybrid propulsion system, a feasibility study on possibilities of electrostatic acceleration of a laser ablation plasma induced from a solid target was conducted. Energy distributions of accelerated ions were measured by a Faraday cup. A time-of-flight measurement was also conducted for ion velocity measurement. It was found that an average speed of ions from a pure laser ablation in this case was about 20 km/sec for pulse energy of 40 $\mu$J/pulse with pulse width of 250 psec. On the other hand, through an electrostatic field with a + I ,000 V electrode, the speed could be accelerated up to 40 km/sec. It was shown that the electrode with positive potential was more effective than that with negative potential for positive-ion acceleration in laser induced plasma, or pulsed plasma, in which ions were induced with the Coulomb explosion following electrons. In addition, the ion-acceleration or deceleration strongly depended on conditions of pairs of inner diameter and electrodes gap.

To realize magnetic sails, momentum of the solar wind should be efficiently transferred to a spacecraft via magnetic field, which is produced around a spacecraft. In this paper, two important physical processes are addressed: 1) diffusive processes caused by plasma turbulence at the magnetospheric boundary around the spacecraft; and 2) field aligned current loops that will electrically connect the magnetospheric boundary and the spacecraft. The idea of the magnetic sails will be demonstrated by an experimental simulator, in which a fast plasma beam will penetrate into a dipole magnetic field. For that purpose, the two important physical processes should be scaled down to a small laboratory experiment in a space chamber. From the scaling considerations, the interaction can be scaled down if high-speed and high-density $(10^{19}m^{-3})$ plasma jet is used with 1-T-class magnetic field.

Discharge current oscillation in a 20KHz range is a serious problem for Hall thruster performance, In our previous work, a possibility of controlling the oscillation amplitude by increasing the speed of neutral particles incoming to the ionization zone was predicted in our previous work. In this paper, the effects of diameter of anode orifice on the oscillation phenomena and the thrust performance were evaluated experimentally. The experimental results show that the measured amplitude of oscillation becomes smaller as the diameter of anode orifice. However, the larger orifice makes thrust performance lower. The results of numerical analysis of neutral particles show that these tendencies have much to do with neutral properties.

A numerical simulation is presented for investigating the effects of pressure ratio of $D_2$ injector to supersonic nozzle on the population inversion in the DF chemical laser cavity, while a lasing concurrently takes place. The laser beam is generated between the mirrors in the cavity and it is important to obtain stronger population inversion and more uniform distribution of the excited molecules in the laser cavity in order to produce high power laser beam with good quality. In this study, these phenomena are investigated by means of analyzing the distributions of the DF excited molecules and the F atom used as an oxidant, while simultaneously estimating the maximum small signal and saturated gains and power in the DF chemical laser cavity. For the numerical solution, an 11-species (including DF molecules in various excited states of energies), 32-step chemistry model is adopted for the chemical reaction of the DF chemical laser system. The results are discussed by comparison with two $D_2$injector pressure cases; 192 torr and 388.64 torr. Major results reveal that in the resonator, stronger population inversions occur in the all transitions except DF(1)-DF(0), when the $D_2$injection pressure is lower. But, the higher $D_2$injection pressure provides a favorable condition for DF(1)-DF(0) transition to generate the higher power laser beam. In other words, as the pressure of $D_2$injector increases, the maximum small signal gain in the $V_{1-0}$ transition, which is in charge of generating most of laser power, becomes higher. Therefore, the total laser beam power becomes higher.r.

An assessment of a novel laser-electric hybrid propulsion system was conducted, in which a laser-induced plasma was induced through laser beam irradiation onto a solid target and accelerated by electrical means instead of direct acceleration only by using a laser beam. A fundamental study of newly developed rectangular laser-assisted pulsed-plasma thruster (PPT) was conducted. On discharge characteristics and thrust performances with increased peak current compared to our previous study to increase effects of electromagnetic forces on plasma acceleration. Maximum peak current increased for our early study by increasing electromagnetic effects in a laser assisted PPT. At 8.65 J discharge energy, the maximum current reached about 8000 A. Plasma behaviors emitted from a thruster in various cases were observed with an ICCD camera. It was shown that the plasma behaviors were almost identical between low and high voltage cases in initial several hundred nanoseconds, however, plasma emission with longer duration was observed in higher voltage cases. Canted current sheet structures were also observed in the higher voltage cases using a larger capacitor. With a newly developed torsion-balance type thrust stand, thrust performances of laser assisted PPT could be estimated. The impulse bit and specific impulse linearly increased. On the other hand, coupling coefficient and the thrust efficiency did not increase linearly. The coupling coefficient decreased with energy showing maximum value (20.8 ？Nsec/J) at 0 J, or in a pure laser ablation cases. Thrust efficiency first decreased with energy from 0 to 1.4 J and then increased linearly with energy from 1.4 J to 8.6 J. At 8.65 J operation, impulse bit of 38.1 ？Nsec, specific impulse of 3791 sec, thrust efficiency of 8 %, and coupling coefficient of 4.3 ？Nsec/J were obtained.

A preliminary study on. field emitter array cathodes was conducted aiming at applying for electrodymanic tether (EDT) propulsion systems. The EDT propulsion systems are assumed to use for active removal systems of post-mission spacecraft, which would otherwise become space debris. A survey on field emit-ter array cathode technology was conducted, and it showed that carbon nanotube (CNT) emitters are suit-able to EDT application. Trial fabrications and evaluation tests of CNT emitters were conducted, which demonstrated a target emission current density of 10 ㎃/$\textrm{cm}^2$. It was found out that the most important technical issue for developing CNT emitters is to improve the performance against voltage breakdown between the emitter and the opposite electrode.

The microwave discharge ion engine generates plasmas of both the main ion source and the neutralizer using 4㎓ microwave without discharge electrodes and hollow cathodes, so that long life and durability against oxygen and air are expected. The MUSES-C “HAYABUSA” asteroid explorer installing four microwave discharge ion engines “$\mu$10s” was launched into deep space by M-V rocket No.5 on May 9, 2003. After vacuum exposure and several runs of baking for reduction of residual gas the ion engine system established the continuous acceleration of the spacecraft toward the asteroid “ITOKAWA”. The Doppler shift measurement of the communication microwave revealed the performance of ion engines, which is 8mN thrust force for a single unit with 3,200sec specific impulse at 23mN/㎾ thrust power ratio. At the end of 2003 the accumulated operational time exceeded 8,000 hour and unit. HAYABUSA will execute the Earth swing-by on June 2004 and arrive at the asteroid in 2005 and return to Earth in 2007.

On 9/5/2003, the planet probe “HAYABUSA” as MUSES-C project was launched by The Institute of Space and Astronautical Science. “HAYABUSA” has microwave discharge ion engines and these engines are characterized by their high efficiency and specific impulse in comparison with chemical engine. A large ion engine can be used as a planet explorer, while a small ion engine can be used as attitude control of small satellite. We have been developing a high thrust density microwave discharge ion engine using “Multi-Monopole Antenna”. The performance of this engine are: ion cost of 344W/A, propellant utilization efficiency of 52% and thrust density of 0.055mN/$\textrm{cm}^2$ for Kr gas flow rate of 2.5sccm, microwave(2.45㎓) power of 32W and acceleration voltage of l.4㎸.

This paper shows an application of single crystal metals and Single Shell Polymer Concentrator (SSPC) to Solar Thermal Propulsion (STP). Based on it, we fabricated a breadboard model of STP system (STP-BBM) for microsatellites. We also proposed Eco-Friendly End-of-Life De-Orbiting (EFELDO) by using such a high performance STP system.

Laser Propulsion is a device that generates thrust using laser energy. Laser-driven In-Tube Accelerator (LITA) has been developed at Tohoku University. LITA is a laser propulsion system that accelerates an object not in an open air but in a tube. Experiments of vertical launching and pressure measurement on the tube wall were carried out and in order to observe the initial state of plasma and blast wave, the visualization experiment was carried out using the shadowgraph method. In this study, the time variation of pressure on the tube wall is numerically simulated solving Euler equation. In order to model the laser energy, heat source function added to the frozen flow Euler equation. Plasma size from the shadowgraph images was used for the initial condition of laser energy input. For verification of the modeling, these results were compared with the previous experimental and numerical results. From these verifications, an analysis of LITA performance will be investigated.

Experiment of laser propulsion in free flight has never been conducted. At Institute of Fluid Science (IFS), Tohoku University, propulsive impulse generation by focusing on a rest projectile was demonstrated. Based on the ideas obtained from this experiment, experiment of laser propulsion of a projectile in flight by focusing $CO_2$ laser beam is being prepared for. The objective velocity increment in experiment is about 50 m/s.

Spectroscopic and electrostatic probe measurements were made to examine plasma characteristics with or without a metal plate for a 10-㎾-class direct-current arcjet Heat fluxes into the plate from the plasma were also evaluated with a Nickel slug and thermocouple arrangement. Ammonia and mixtures of nitrogen and hydrogen were used. The NH$_3$ and $N_2$+3H$_2$ plasmas in the nozzle and in the downstream plume without a plate were in thermodynamical nonequilibrium states. As a result, the H-atom electronic excitation temperature and the $N_2$ molecule-rotational excitation temperature intensively decreased downstream in the nozzle although the NH molecule-rotational excitation temperature did not show an axial decrease. Each temperature was kept in a small range in the plume without a plate except for the NH rotational temperature for NH$_3$ gas. On the other hand, as approaching the plate, the thermodynamical nonequilibrium plasma came to be a temperature-equilibrium one because the plasma flow tended to stagnate in front of the plate. The electron temperature had a small radial variation near the plate. Both the electron number density and the heat flux decreased radially outward, and an increase in H$_2$ mole fraction raised them at a constant radial position. In cases with NH$_3$ and $N_2$+3H$_2$ a large number of NH radical with a radially wide distribution was considered to cause a large amount of energy loss, i.e., frozen flow loss, for arcjet thrusters.

On a liquid propellant PPT, the discharge processes that discharge was initiated and plasma was accelerated was observed by using a ultra high speed camera. Liquid propellant PPT is a pulsed plasma thruster using liquid as propellant. Our past study showed the successful operation of liquid propellant PPT and the thruster showed high specific impulse. However, its acceleration mechanism has not been clarified. In this study we observed the plasma acceleration processes in order to deepen our understanding of the acceleration mechanism.

Motion of a bubble inside narrow tube is numerically studied. The numerical code assumes axi-symmetric incompressible flow field. The surface of the bubble is captured by VOF (Volume Of Fluid) method, and it is advected by MARS (Multiphase Advection and Reconstruction Scheme). Air bubble inside water is first studied, and it was found that a strong vortex, which is induced by the pressure difference caused by the surface tension, is formed at the rear part of the bubble. Then flow parameters are parametrically varied to understand the correlation between the bubble shape, the bubble velocity, and the flow parameters. The parametric study revealed that the aspect ratio of the bubble mainly depends on We number, and the oscillation of the bubble speeds is dependent on Re number.

The unsteady aerodynamic characteristics of an aerofoil gradually accelerating or decelerating at subsonic speeds are investigated through two-dimensional, unsteady, compressible Navier-Stokes simulations. An acceleration factor is defined to provide various acceleration or deceleration characters of the time-dependent flow over the aerofoil. The results show that an increase in the absolute value of the non-dimensional acceleration factor leads to a lesser change in the location and range of flow featues such as shockwave and boundary layer separation in a specific time range. Generally, the gradual speed-up and speed-down of the subsonic aerofoil results in different aerodynamic characteristics whose changes are more significant at angles of attack.

When a body falls in fluid, the body often experiences autorotations, namely, various kind of rotating motions, such as tumbling, flat spin and coming. Tumbling is a rotating motion with an axis perpendicular to a falling direction. Tumbling is a very important phenomenon in aeronautical and space engineering, ballistics and meteorology. For example, when an satellite re-en-tries into the atomosphere, its body collapses into many fragments which are disperse in the wide range of field. Some fragments fall in tumbling motion. Then tumbling is useful to predict fragment's motion.(omitted)

A computational study on the supersonic flow around the lateral jet controlled missile has been performed. Case studies have been performed by comparing the normal force coefficient and the moment coefficient of a missile body for several different jet flow conditions, angle of attacks, circumferential jet locations, and spouting jet angles. For the several different jet flow conditions, which include the jet pressure, the jet Mach number, and the corresponding jet mass flow rate, the results show that the normal force coefficient is almost proportional to the jet thrust but the moment coefficient is not. Distinctly different flow phenomena can be noticed as the pressure ratio and the jet Mach number increase. By investigating the angle of attack effect to the normal force and the pitching moment, it has been identified that the normal force and the pitching moment show nonlinearity with respect to the angle of attack. From the detailed flow field analyses with respect to the jet flow conditions and the angle of attacks, it is verified that most of the normal force loss and the pitching moment generation are taken place at the low-pressure region behind the jet nozzle. Furthermore, the normal force and the pitching moment characteristics of the missile have been identified by comparing different circumferential jet locations and spouting jet angles.

The multi-phase flow analysis in a solid rocket motor is very important when performing the performance of a motor, and prediction of nozzle ablation. However, only in consideration of regular power, it has analyzed as power which a metal particle receives from a flow until now. We conduct analysis and an experiment about the virtual mass clause which will influence at the place where acceleration is big. We aim at the improvement in accuracy of multi-phase flow analysis from the result.

Gas turbine engine simulation in terms of transient, steady state performance and operational characteristics is complex work at the various engineering functions of aero engine manufacturers. Especially, efficiency of control system design and development in terms of cost, development period and technical relevance implies controlling diverse simulation and identification activities. The previous engine simulation has been accomplished within a limited analysis area such as fan, compressor, combustor, turbine, controller, etc. and this has resulted in improper engine performance and control characteristics because of limited interaction between analysis areas. In this paper, we propose a new simulation methodology for gas turbine engine performance analysis as well as its digital controller to solve difficulties as mentioned above. The novel method has particularities of (ⅰ) resulting in the integrated control simulation using almost every component/module analysis, (ⅱ) providing automated math model generation process of engine itself, various engine subsystems and control compensators/regulators, (ⅲ) presenting total sophisticated output results and easy understandable graphic display for a final user. We call this simulation system GT3GS (Gas Turbine 3D Graphic Simulator). GT3GS was built on both software and hardware technology for total simulation capable of high calculation flexibility as well as interface with real engine controller. All components in the simulator were implemented using COTS (Commercial Off the Shelf) modules. In addition, described here includes GT3GS main features and future works for better gas turbine engine simulation.

The purpose of this work is to develop an effective active control system for combustion instabilities of premixed combustors. For the first step, the natural modes of combustion oscillation were investigated for a methane-air premixed combustor and the controls by secondary fuel injection were examined. The main premixed flame is stabilized by a swirler with orifices for secondary injection installed on the central hub. For sensing purposes, a pressure transducer and a chemiluminescence sensor were placed on the appropriate positions. The acoustic characteristics and the source of the oscillation were analyzed by those signals. To test the controllability, two methods of actuations by secondary fuel injection were examined. One is the open loop control and the other is the closed loop control. The comparison of the reduction levels of p $_{rms}$ shows that the closed loop control with a phase-shift injection performs best in this condition.ition.n.

A fuel spiking test was performed to measure the surge margin of the compressor in a gas turbine engine. During the test, fuel spiking signal was superimposed on the engine controller demand and the mixed signals were used to control a fuel line servo-valve. For the superimposition, a subsystem composed of a fuel controller and a function generator was used. During the fuel spiking test, the original scheduled fuel signals and the modified signals were compared to guarantee the consistency excluding the spiking signals. The spiking signals were carefully selected to maintain the engine speed constant. The fuel spiking effects were checked by three dynamic pressure sensors. Sensors were placed before the servo-valve, after the servo-valve, and after the compressor location, respectively. The modulations of the spiking signal duration and fuel flow rate were examined to make the- operating point approach the surge region. The real engine test was performed at the Altitude Engine Test Facility (AETF) in Korea Aerospace Research Institute (KARI). In the real engine test, fuel spiking signals with 25~50 ㎳ of spiking signal time and 17~46 % of base fuel flow rate condition were used. The dithering signal was 5~6 ㎃ at 490 Hz. The test results showed good agreement between the fuel spiking signals and the fuel line pressure signals. Also, the compressor discharge pressure signals showed fuel spiking effects and the changes of the operating point on the compressor characteristic map could be traced.

An unique ceramic material produced through unidirectional solidification with eutectic composition of two-phase oxides was introduced recently. This composite material has the microstructure of coupled networks of two single crystals interpenetrate each other without grain boundaries. Depending on this microstructure this material, called Melt Growth Composite (MGC), can sustain its room temperature strength up to 1$700^{\circ}C$ (near its melting point) and offer strong oxidization-resistant ability, making its characteristics quite ideal for the gas turbine application. The research project on MGC started in 2001 with the objective of establishing component technologies for MGC application to the high temperature components of the gas turbine engine. MGC turbine nozzles are expected to improve efficiency of gas turbine. However, reduction of the thermal stress is required since high thermal stress is easily generated in MGC turbine nozzles due to temperature distribution. Firstly, the hollow nozzle shape was optimized to reduce thermal stress using numerical analysis. From the results of the first hot gas flow tests, the thermal stress due to span-wise temperature distribution was required to be reduced, and separated nozzle to three pieces was designed. This was tested in hot gas flow at 140$0^{\circ}C$ level, and temperature distributions on the nozzle surface were obtained and stress field was evaluated.

The performance characteristics of partial admission supersonic turbines are analyzed by using the commercial CFD program FLUENT6.0. The governing equations were discretized with Euler implicit method in time and 2nd-order upwind scheme of FVM in space. The k-$\varepsilon$ turbulence model was utilized to describe the turbulent flow field. In order to investigate the nozzle--rotor interactions and the effect of partial admission, the flows in supersonic turbine rotor cascades with a nozzle are computed. Extensive computations of partial admission supersonic turbines provide the shock structures and flow patterns in the nozzle and rotor. It is clearly shown that the nozzle flow is highly affected by the shocks or expansion waves propagated from the rotor leading edge. And the rotor flow is also affected by the shocks or wakes originated from the nozzle.

The three-dimensional woven fabric C/SiC composites blisk turbine rotor model was evaluated. The spin tests of the blisk model were performed to measure strain distributions at the room temperature. The rotational strength of the blisk model could be improved by the fiber addition. But, there are still more researches to be done.

This study proposes a interim development result for the l-㎾ class small wind turbine system, which is applicable to relatively low wind speed regions like Korea and has the variable pitch control mechanism. In the aerodynamic design of the wind turbine blade, parametric studies were carried out to determine an optimum aerodynamic configuration which is not only more efficient at low wind speed but whose diameter is not much larger than similar class other blades. A light composite structure, which can endure effectively various loads, was newly designed. In order to evaluate the structural design of the composite blade, the structural analysis was performed by the finite element method. Moreover both structural safety and stability were verified through the full-scale structural test.

Numerical and experimental investigation were car-ried out to clarify the flow structure of underexpanded jet from a square nozzle. The square nozzle rep-resents one of the clustered combustors of a linear aerospike engine. From the numerical results, the three-dimensional shock wave of the underexpanded square jet was found to be composed of two shocks. One is the intercepting shock which corresponds to the shock observed in two-dimensional planar jet. The other is the recompression shock divided into two types. The expansion fans coming from the nozzle edges interact with each other at the comers of the nozzle exit, and overexpanded regions are generated. Therefore one of the two recompression shocks is formed at the comers of the nozzle exit behind the overexpanded regions. As the jet goes downstream, the overexpanded regions grow larger to coalesce at the symmetry planes. Then, the other type of the recompression shock is generated. The three-dimensional shock structure formed by the intercepting shock and the recompression shocks dominates the expansion of the jet boundary. The shock detection algorithm us-ing CFD results was developed to reveal the relation between the shock waves and the jet boundary, and it was found that the cross-sectional jet shape becomes cross-shape. The key features observed in the numerical investigation were verified by the experimental results. The shock structure at the diagonal plane was in good agreement with the experimental schlieren images. Moreover, the cross-sections visualized by the Mie scattering method confirmed that the cross-section of the jet becomes cross-shape.

A cone cylinder is used to obtain variable operation conditions for the sonic ejector-diffuser system. The cone cylinder is designed to be shifted upstream and downstream to change the ejector throat area ratio, thus obtaining variable mass flow rates. The present study investigates the effects of ejector throat area ratio and operating pressure ratio on the entrainment of secondary stream for the variable sonic ejector system. The study is carried out experimentally. The ejector throat area is varied at the range from Ψ= 11.88 to 66.69, and the operating pressure ratio is changed from $P_{op}$ / $P_{a}$=1.25 to 9.0. The results show that the variable sonic ejector system can be operated to obtain specific entrainment ratio of secondary stream by altering the ejector throat area ratio and operating pressure ratio.o.

Characteristics of a liquid-vapor interface where a nonequilibrium condensation flow exists are considered based on molecular dynamics simulations, The condensation coefficient, the velocity distributions of the reflected and evaporated molecules and the number flux of the evaporated molecules are compared with those under the liquid-vapor equilibrium. The comparison shows that the condensation coefficient under the nonequilibrium condensation is slightly larger and the number flux of the evaporated molecules is considerably smaller than those under the liquid-vapor equilibrium. The net condensation flux under the nonequilibrium condensation is underestimated if it is evaluated from the condensation coefficient and the number flux of the evaporated molecules under the liquid-vapor equilibrium. However the underestimation is relatively small.

Large-Eddy Simulation (LES) is applied for the simulation of compressible flat plate boundary with Reynolds number up to 5 X 10$^{5}$ . Numerical examples include shock/boundary layer interaction and boundary layer transition, aiming future application to the analysis of transonic fan/compressor cascades. The present LES code uses hybrid com-pact/WENO scheme for the spatial discretization and compact diagonalized implicit scheme for the time integration. The present code successfully predicted the bypass transition of subsonic boundary layer. As for supersonic turbulent boundary layer, mean and fluctuation velocity of the attached boundary, as well as the evolution of the friction coefficient and the displacement thickness both upstream and downstream of the separation region are all in good agreement with experiment. The separation point also agreed with the experiment. In the simulation of the shock/laminar boundary layer interaction, the dependence of the transition upon the shock strength is reproduced qualitatively, but the extent of the separation region is overpredicted. These numerical examples show that LES can predict the behavior of boundary layer including transition and shock interaction, which are hardly managed by the conventional Reynolds-averaged Navier-Stokes approach, although there needs to be more effort before achieving quantitative agreement.

This paper reports the grid resolution issues on the supersonic mixing simulation inside the engine for future aerospace vehicles. Unstructured finite volume method is used for the simulations. Three types of grids are used, namely, hybrid unstructured grids composed of prism and tetrahedron cells, locally refined grids, and hexahedral grids. Hexahedral grids are used to take advantage of fine distribution naturally behind the edge of the ramp where the vortex is generated. These latter two grids show much improved evaluations of the vortex motion and the mixing of the injected and the main flows.

Toasting the 100 year anniversary of controlled, powered flight, the propulsion system used on today's aircraft represents the evolution of jet propulsion based on the gas turbine, first conceived by Whittle and Von Ohain about 70 years ago. In that period, propulsion system concepts have evolved through turbo-props, turbo-jets, low by-pass ratio(BPR) turbofans to today's high BPR 2-shaft and 3-shaft turbofans. Also, this period has seen remarkable progress in the performance, reliability environmental compatibility of these propulsion systems.(omitted)

“Flameless combustion” of lean to ultra lean mixtures, supported by high-temperature burned gas, can resolve the dilemma between complete combustion versus ultra-low NOx emissions in gas turbine combustors. The characteristics of NOx emissions and combustion in “lean-lean” two-stage combustion were investigated for fuel vapor and droplets / air mixture jets injected from the main injection tube that was placed perpendicular to the combustor wall into the primary hot burned gas prepared by combustion of lean mixtures on a perforated flame holder. The present results clearly show that the ultra-low NOx combustion supported by the reaction of lean mixtures well mixed with the hot burned gas from the primary stage is much more advantageous in achieving ultra-low NOx emissions while maintaining high combustion efficiency.

An axially staged combustor equipped with an LPP combustion system and CMC liner walls has been investigated for stable combustion and low NOx emissions for the ESPR project. Several fuel injectors were designed and manufactured for the LPP burner, and single sector combustor tests were conducted to evaluate fundamental combustion characteristics such as emissions, instabilities, auto-ignition, and flash back at typical operating conditions from idle to Mn 2.2 cruise. The latest test results showed that the LPP burner had a good potential for the low NOx target. It was also found that the NOx emission level was greatly affected by a distortion in the air flow velocity field upstream of the LPP burner due to the diffuser and fuel feed arm. The CMC material was investigated to apply for the high temperature and low NOx combustor. Annular combustor liner walls were manufactured with the CMC material, and they have been tested at low pressure conditions to evaluate the soundness of the material and the mounting and seal system. This paper reports the latest research activities on the LPP combustion system and CMC liner walls for the ESPR project.

ESPR project started in 1999 with METI and NEDO support as five-years program in order to develop necessary technologies for the next-generation SST engine. In ESPR program, jet noise reduction technologies are focused as environmentally compatible technologies, which are critical to realize next-generation SST. In designing a lobed mixer nozzle which is a jet noise suppression system, there are many difficulties to understand the detailed flow phenomena occurred in the system because of its complexity. Large Eddy Simulation (LES) was applied to the lobed mixer nozzle flow analysis in ESPR project. The results demonstrated that LES approach was capable of predicting mixing characteristics of a complicated flow.

In gas turbines, excess air for combustion is available and therefore lean premixed combustion is the most promising approach to the significant reduction of thermal NOx emissions. At lean conditions, however, flame stability is inherently worse and hence combustion tends to be incomplete. Efforts have been devoted toward extending the operating range of complete combustion at leaner conditions. One of them is the lean-lean two-stage combustion where lean to ultra-lean secondary mixtures are mixed with the hot burned gas from the primary stage. Conventional flame combustion or flameless reaction are initiated depending on the conditions of the secondary zone. In the first part of the present study, the effects of fuel injection on the emissions and flame stability were investigated for a single tubular flame, In the second part, the emissions and flame stability were studied for a two-stage combustor with secondary mixture injected through the tangential slots on a cylindrical combustor wall. The effects of the ratio of air flow rates to the primary and secondary zones on the emissions and combustion characteristics were investigate.

This paper describes the characteristics of the unsteady pressure on the stator surface induced by rotor viscous wakes. The primary object of this study is to investigate the effects of axial spacing between the rotor and the stator and three-dimensional vane geometries such as stator sweep and stator lean on the unsteady pressure fluctuations on the stator vane. To predict these fluctuations, unsteady three-dimensional Navier-Stokes analyses are performed. Furthermore, a three-dimensional analytical method using unsteady lifting-surface theory is also used to elucidate the mechanism of interaction of passing rotor wakes with downstream stator. Five different fan configurations with three sets of stator geometries, which are three radial stator configurations with different axial spacing, the swept stator and the swept and leaned stator, are used for this study. It is found that, in axial spacing between rotor and stator, the effect of radial phase skew of incoming rotor wake is important for the reduction of the induced unsteady pressure in addition to the rotor wake decay. It is also shown that incorporation of stator sweep and lean is effective to reduce this unsteady pressure.

Reactions of NO$_{x}$, HO$_{x}$ and $O_3$ chemistry in a diffusion process of the exhaust plume under a stratospheric condition were investigated numerically. Expanding Box method was used to assess the effects of exhaust gases from a stratospheric flight system on $O_3$ depletions. Sensitivity analysis was also performed to identify prime reactions of $O_3$ depletions in an exhaust plume right after the nozzle. In addition, a calculation of reactive flows in stratospheric condition was performed to investigate the characteristics of reactions in a plume. As a result of this study, prime reactions of NO$_{x}$, HO$_{x}$ and $O_3$chemistry in an exhaust plume were identified, and fundamental behavior of chemical species were examined in a exhaust plume.t plume.

Sawtooth mixing device used in a non-premixed burner is evaluated for flame stabilization and NO$_{x}$ reduction. Three mixers with different blade angles are tested. Methane is delivered through the fuel jet and air passes through the co-flow annulus. The flame mode changes (attached flame, lifted flame and extinction) against the fuel flow speed are measured, and the stability diagram is drawn. Moreover, by traversing thermocouple and sampling probe in the flame, the distribution of temperature and NO$_{x}$ mole fraction are measured. With the change in blade angle, flame shape, flame stabilization, the distribution of temperature and NO$_{x}$ mole fraction are changed considerably.rably.

Lift Fan Engines of JAXA's conceptual Jet VTOL aircraft have a very small bellmouse shape air intake, which make some differences in aerodynamic design of the blades. To obtain a better rotor or stator blade design, this paper performs a numerical analysis of the throughflow on a lift fan as a two-dimensional axisymmetrical flow. Based on the last report focusing on the air intake's influence on the throughflow, a more realistic bellmouse air intake case is treated to reconsider the influence on the throughflow by the small bellmouse air intake. Three work input patterns are tested to reduce some problematic influences on the throughflow or blade designs. The obtained result shows one of acceptable blade designs for the lift fan engine.

A Propulsion System of the CRW(Canard Rotor Wing) type UAV(Unmanned Aerial Vehicle) was composed of the turbojet engine to generate the propulsive exhaust gas, and the duct system including straight bent ducts, tip-jet nozzles, a master valve and a variable main nozzle for three flight modes such as lift/landing mode, low speed transition flight mode and high speed forward flight mode. In this study, in order to operate safely the propulsion system, the dynamic Performance behavior of the system was modeled and simulated using the SIMULIN $K^{ }$, which is the user-friendly GUI type dynamic analysis tool provided by MATLA $B^{ }$. In the transient performance model, the inter-component volume model was used. The performance analysis using the developed models was performed at various flight condition, valve angle positions and fuel flow schedules, and these results could set the safe flight mode transition region to satisfy the inlet temperature overshoot limitation as well as the compressor surge margin. Performance analysis results using the SIMULIN $K^{ }$ performance program were compared with them using the commercial program GSP.m GSP.

We analytically estimated the propulsive performance of pulse detonation engines (PDEs) in three cases, which were (1) a fully-fueled simplified PDE, (2) a partially-fueled simplified PDE, and (3) a PDE optimized as a system. The results of the model analyses in the cases of (1) and (2) were in good agreement with published experimental data which were obtained by using simplified PDEs. The comparison between the results of the analyses of simplified PDEs and those of optimized PDE systems showed that specific impulse would become higher by about 10-20% due to PDE-system optimization.

A steady-state/transient performance simulation model was newly developed for the propulsion system of the CRW (Canard Rotor Wing) type UAV (Unmanned Aerial Vehicle) during flight mode transition. The CRW type UAV has a new concept RPV (Remotely Piloted Vehicle) which can fly at two flight modes such as the take-off/landing and low speed forward flight mode using the rotary wing driven by engine bypass exhaust gas and the high speed forward flight mode using the stopped wing and main engine thrust. The propulsion system of the CRW type UAV consists of the main engine system and the duct system. The flight vehicle may generally select a proper type and specific engine with acceptable thrust level to meet the flight mission in the propulsion system design phase. In this study, a turbojet engine with one spool was selected by decision of the vehicle system designer, and the duct system is composed of main duct, rotor duct, master valve, rotor tip-jet nozzles, and variable area main nozzle. In order to establish the safe flight mode transition region of the propulsion system, steady-state and transient performance simulation should be needed. Using this simulation model, the optimal fuel flow schedules were obtained to keep the proper surge margin and the turbine inlet temperature limitation through steady-state and transient performance estimation. Furthermore, these analysis results will be used to the control optimization of the propulsion system, later. In the transient performance model, ICV (Inter-Component Volume) model was used. The performance analysis using the developed models was performed at various flight conditions and fuel flow schedules, and these results could set the safe flight mode transition region to satisfy the turbine inlet temperature overshoot limitation as well as the compressor surge margin. Because the engine performance simulation results without the duct system were well agreed with the engine manufacturer's data and the analysis results using a commercial program, it was confirmed that the validity of the proposed performance model was verified. However, the propulsion system performance model including the duct system will be compared with experimental measuring data, later.

The Japan Aerospace Exploration Agency (JAXA) has proposed new vertical take-off and landing (VTOL) aircraft known as the Jet-VTOL aircraft shown in Fig.1. The Jet-VTOL aircraft is based on a canard wing configuration. The aircraft has the clustered lift-fans mounted near the center of gravity for vertical flight, and has the clustered fans mounted beside the vertical tail for cruise flight. Both fans are driven by the core engine mounted inside the aft end of fuselage. The propulsion system is innovative and attractive not to be seen even in the world.

Hybrid rocket had many advantage with compared to solid and liquid rockets. However, the engines have not yet been used in practical rocket systems, due mainly to the disadvantage of hybrid combustion, such as low fuel regression rate. In this study, lab-scale hybrid motor was designed and manufactured. And the methods of regression rate improvement were considered. Test firings with thrusts up to 300 N were conducted with GOX and transparent PMMA. Thrust was calculated with the pressure of the combustion chamber and the regression rate was measured in with variation of oxidizer flow rate. The regression rates showed a strong dependency on GOX mass flux. The frequency analysis technique of the bulk-mode oscillation of motor was applied to a hybrid rocket motor and was based on the principle that this frequency was inversely proportional to the square root of the chamber volume. Several problems and solutions of operating hybrid rocket were presented.

For a simple one-dimensional ignition problem a mathematical model is described to investigate the difficulties in numerical simulations. Some computation results are obtained and comparison is made with analytical solution. Discussions are made on topics such as 1) coordinate transformation, 2) gas-phase and solid-phase analysis; (divergence form of the governing system, a finite-volume discretization, implicit time integration, upwind split flux, spatial accuracy improvement are described. Mass, reagent mass, and energy conservations are solved.), and 3) method to determine quantities on the burning surface (matching). Results obtained for small values of the non-dimensional pressure show a steady-combustion and good agreement with the analytical solution. Numerical instability appeared for larger values of the pressure, discussion on the cause of the problem is made. This effort is a part of a study of flame spread phenomena on solid propellant surface.

The present work have been developed the interpretation processor including the behavior of material failure and the separation phenomena under transient dynamic loading (the operation of explosive bolt) using AUTODYN V4.3, SoildWork 2003 and TrueGrid V2.1 programs. It has been demonstrated that the interpretation in ridge-cut explosive bolt under dynamic loading condition should be necessary to the appropriate failure model and the basic stress of bolt failure is the principal stress. The use of this interpretation processor developing the present work could be extensively helped to design the shape and the amount of explosives in the explosive bolt having a complex geometry. It is also proved that the interpretation processor approach is an accurate and effective analysis technique to evaluate the separation mechanism in explosive bolts.

하이브리드 로켓 엔진은 일반적으로 저비용ㆍ저공해ㆍ고안전성 등의 이점이 있지만 연료 후퇴 속도가 느려 연비를 제어할 수 없는 문제점이 일어 아직 실용화되지 않고 있다.(중략)

Theoretical studies have been carried out to examine the influence of the grain geometry-dependent driving forces, which control the internal flow pattern of solid rockets. Numerical studies have been executed with the help of a two-dimensional code. This code solves standard k-omega turbulence equations using the coupled second order implicit unsteady formulation. It has been concluded that the grain port divergence angles have significant leverage on the formation of recirculation bubbles leading for pressure oscillations, flow separation and reattachment. In solid rockets flow reattachment will favour secondary ignition and that will add to the complexity of the starting transient prediction.

The optimal design and combustion analysis of the gas generator for Liquid Rocket Engine (LRE) were performed. A fuel-rich gas generator in open cycle turbopump system was designed for 10ton$_{f}$ in thrust with RP-1/Lox propellant. The optimal design was done for maximizing specific impulse of main combustion chamber with constraints of combustion temperature and power matching required by turbopump system. Design variables were selected as total mass flow rate to gas generator, O/F ratio in gas generator, turbine injection angle, partial admission ratio, and turbine rotational speed. Results of optimal design show the dimension of length, diameter, and contraction ratio of gas generator. Also, the combustion test was conducted to evaluate the performance of injector and combustion chamber. And the effect of the turbulence ring was investigated on the mixing enhancement in the chamber.r.

Numerical stability analysis of one-dimensional axial flow in solid rocket motors is performed based on the Euler equation coupled with an unsteady combustion equation of solid propellant. In order to check the numerical scheme, behavior of a standing wave in a closed tube is examined. A standing wave in solid rocket motor decays or grows depending on the total effect of propellant combustion, nozzle flow, and so on. The stability boundary of the fundamental mode standing wave is determined by changing one of the combustion parameters. In addition growth rates of the wave are calculated numerically in relatively low Mach number flow region for the motors with different port and nozzle throat diameters. The results obtained here agree well with the approximate solution. The same scheme is applied to a motor with shorter length and L＊-instability is observed.

KARI is achieving the KSLV program according to National Space Technology Development Program. In this paper, the authors are intended to introduce the Integrated Power Plant (abb. IPP) test facility which will be constructed for the variety of tests on KSLV program. IPP test facility refers to comprehensive testing equipment for liquid rocket launch vehicle. Using this facility, KARl can verify the adaptiveness of parts and subsystems for launch vehicle and finally can qualify the system characteristics of launch vehicle doing kinds of test including hot firing test. Using this facility, KARI can simulate the vehicle launching circumstances and it make to predict the performance of launch vehicle when its flight test.

A coaxial pulsed plasma thruster (PPT) with a Teflon cavity was designed, and its performance characteristics were examined varying stored energy, cavity length and capacitance. The PPT was tested as the entire system including the discharge circuit, and the results were explained with both the transfer efficiency and the acceleration efficiency. The transfer efficiency is defined as the fraction of energy in capacitors supplied into plasma, and the acceleration efficiency as the fraction of energy supplied into plasma converted to thrust energy. To estimate these efficiencies, the equivalent plasma resistance was defined and calculated using energy conservation during discharge. The equivalent plasma resistance proportionally increased with cavity length, and therefore the current peak increased with decreasing cavity length. The energy density calculated by the transfer efficiency was increased with decreasing cavity length. As a result, higher acceleration efficiency and lower transfer efficiency were obtained with shorter cavity length. Accordingly, there was an optimal cavity length for the thrust efficiency. The specific impulse and the impulse bit per unit stored energy ranged from 390 s and 50 $\mu$ Ns/J for a cavity length of 34 mm to 825 s and 11 $\mu$ Ns/J for a cavity length of 4 mm when the stored energy was fixed to 21.4J. Thus, it was showed that the performance of this PPT approached that of electromagnetic-acceleration-type PPT with decreasing cavity length. The PPT achieved thrust efficiencies of 10-12% at 21.4 J and 6-7% at 5.35 J at cavity lengths between 14 mm and 29 mm.

This paper presents the preliminary fabrication results of colloid thrusters which can provide thrust of the order of micro to milli-Newtons. MEMS technology has been used for fabrication, and four essential fabrication techniques - deep etching with nested masks, isotropic plasma etching, anisotropic reactive ion etching, and direct fusion wafer bonding - have been newly developed. Among diverse models which have been designed and fabricated, the fabrication results of 4-inch wafer-based colloid thrusters are presented.

We are developing a micro-solid propellant rocket array thruster for simple attitude control of a 10 kg class micro-spacecraft. The prototype has ø 0.8 mm solid propellant micro-rockets arrayed at a pitch of 1.2 mm on a 22 x 22 mm substrate. In previous studies, an impulse thrust of 4.6 x 10$^{-4}$ Ns was obtained in vacuum, but we found the problems of unacceptably low ignition success rate and incomplete combustion. This paper describes experiments to improve the ignition rate. In order to achieve this goal, we tried to solidify paste-like ignition aid (RK) on the ignition heaters with strong adhesion. To make the paste-like RK, isoamyl acetate was added to RK powder. We tested 9 rockets, but only 2 rockets were ignited with huge ignition energy. This is because the heat con-duction between the ignition heater and the RK was too low to ignite the RK, since dried RK had a lot of pores. Also, a large cavity was sometimes found just above the ignition heater.

Fabrication technique and performance test of catalytic micro propulsion are treated based on MEMS technology. This propulsion is designed to use hydrogen peroxide as liquid mono-propellant for attitude control of pica-satellite. The propellant is fed into the micro reactor channel and decomposed into hot gas yielding controllable thrust by catalyst. In order to increase the efficiency of the reaction that depends on the contact area of propellant and catalyst, porous surface formation on the channel accompanied by platinum particle deposition has been performed using H$_2$PtCl$_{6}$ solution as a precursor. Several thrusters were fabricated in different concentration of H$_2$PtCl$_{6}$ solution to determine the best quantity of Pt particles. For the comparison of the performance of each thruster, the volume of oxygen generated by the decomposition of hydrogen peroxide and the thrust were measured.red.

Precise position and attitude control of pico-satellite requires huge number of impulses of the order of 10$^{-6}$ Ns. MEMS solid rocket array is a promising propulsion system but the higher degree of miniaturization causes unreliable operation mainly due to quenching. In order to breakthrough this situation, a novel design of solid micro-rocket is proposed, which generates tiny impulses repetitively from a single rocket not from array. This unique micro-rocket is based on the utilization of quenching, which causes propellant reaction to sustain only in a small area. A test chip of a micro solid propellant tank and micro heater array is fabricated and ignition test is conducted. Obtained results show the feasibility of this concept and future direction of this quenching-based propulsion is discussed.

Solid propellants allow thrusters to be light-weight, com-pact and robust because they require neither tank nor valve, Moreover, the solid propellant will not leak, spill or slosh. Consequently, the solid propellant thruster is one of the potential candidates for the microthruster. On the other hand, the control of the solid propellant combustion is difficult, since the conventional solid propellant continues to bum until all the stored propellant is consumed. Although particular devices like thrust reverser were designed to control the combustion, these devices were rarely used in the practical rocket motors. These devices rise thruster weight as well as complicate the thruster operation. In this study, a solid propellant microthruster using laser sustained combustion was designed in order to develop a high-efficiency microthruster overcoming the previously-mentioned difficulty. This designed thruster has semiconductor lasers and non-self-combustible solid propellants in addition to the conventional solid propellant thruster. In this designed thruster, the semiconductor laser controls the combustion of the non-self-combustible solid propellant. In order to demonstrate that the solid propellant combustion is controllable with laser, some non-self-combustible solid propellants were irradiated with the laser at a back-pressure of about 1㎪. A 40-W class Neodymium Yttrium Aluminum Garnet (ND:YAG) laser was used as a tentative alternate to the semiconductor laser. This experiment has shown that the solid propellant combustion was controllable with 10- W class laser irradiation.

In this study, microfabrication of a micro-arcjet nozzle with Fifth-harmonic generation Nd：YAG pulses (wavelength 213 nm) and its thrust performance tests were conducted. A micro-arcjet nozzle was machined in a 1.2 mm thick quartz plate. Sizes of the nozzle were 0.44 mm in width of the nozzle exit and constrictor diameter of 0.1 mm. For an anode, a thin film of Au (~100 nm thick) was deposited by DC discharge PVD in vacuum on divergent part of the nozzle. As for a cathode, an Au film was also coated on inner wall surface. In operational tests, a stable discharge was observed for mass flow of 1.0mg/sec, discharge current of 6 ㎃, discharge voltage of 600 V, or 3.6 W input power (specific power of 3.6 MW/kg). In this case, plenum pressure of the discharge chamber was 80 ㎪. With 3.6 W input power, thrust obtained was 1.4 mN giving specific impulse of 138 sec with thrust efficiency of 24 %.

Miniaturization of subsystems including propulsion systems is recent trends in spacecraft technology. Small space vehicle propulsion is not only a technological challenge of a scaling system down, but also a combination of fundamental flow/combustion constraints. In this paper, physical constraints of micronozzle for cold gas micro-thruster are reviewed and discussed. Method to measure small thrust are also described.

There are many difficulties in realizing Ultra-micro gas turbine system. Among them, the effects of tip clearance upon the micro turbine flowfield are discussed in this paper. The flowfield was investigated numerically with the Reynolds-averaged three-dimensional thin-layer Navier-Stokes equations. Calculations were conducted with clearance height from 0％ to 10% of the passage height. Leakage mass flow and deterioration of efficiency are proportional to the clearance height for the clearance height larger than 4%. However, in the case of 2% clearance, leakage flow is significantly reduced due to relative motion of the casing and as a result deterioration of efficiency is very small. It is difficult to control tip clearance in micro turbines, but the results of this study indicate that if the clearance height is controlled within a few per-cent of passage height, deterioration of stage performance will be small.

초소형위성은 개발기간의 단축에 따른 발사 비용의 저감과 업무와 기능의 분산화에 따른 발사 실패의 리스크 저감 등 여러 장점이 있기 때문에 그 필요성이 높아지고 있다.(중략)

In this paper are presented a concept of a new supersonic air inlet, which is designated a Multi-Row Disk (MRD) inlet, aiming at performance improvement under off-design conditions, and results of wind tunnel tests examined performance characteristics of the MRD inlet. The MRD inlet is frequently called ‘a skeleton inlet’ because of its appearance. The performance of a conventional axisymmetric inlet with a solid center body (spike) deteriorates under off-design Mach number conditions. It is due to the fact that total pressure recovery (TPR) governed by the throat area of inlet and mass capture ratio (MCR) governed by an incidence position of an oblique shock from the spike tip into the cowl can not be controlled independently in such air inlet. The MRD inlet has the spike that is composed of a tip cone and several disks arranged downstream of it, based on the experimental fact that several deep cavities on a conical surface have little negative effect on the boundary layer growth. The overall spike length of the MRD inlet is adjustable to the given flight speed by changing space between disks so that a spillage flow can be controlled independently from controlling the throat area. It could be made clear from the result of wind tunnel tests that the MRD inlet improves TPR by 10% compared with a conventional inlet with a solid spike under off-design conditions.

In the present study were examined numerically and experimentally the off-design performance characteristics on an axisymmetric plug nozzle with variable throat area. In this nozzle concept, its throat area can be changed by translating the plug into the axial direction. First, a mixed-expansion plug nozzle, in which two expansion parts are arranged both inside and outside, was designed by means of the method of characteristics. Second, the CFD analysis was verified by the cold-flow wind tunnel test. Third, its performance characteristics were evaluated over a wide range of pressure ratio from half to double throat area through the design point, using the CFD code verified by the wind tunnel tests. It was made clear from the study that not so critical thrust efficiency losses were found and the maximum thrust efficiency loss was at most approximately 5 % under off-design conditions without external flow. This result shows that a plug nozzle can give the altitude compensation even under off-design geometry operations. However, shock waves were observed in the inner expansion part under the doubled throat area operation and thus some thermal problems may be caused on the plug surface. Furthermore, collapse of cell structure on the plug surface was observed with external flow (around Mach number 2.0) as it became lower pressure ratio below the design point and the fact may result in big efficiency loss regardless of geometrical configuration.

Hypersonic wind tunnel test of the rectangular variable geometry intake is performed. For realization of a Precooled turbojet engine, development of a hypersonic ramjet engine is planned. To investigate performance of the intake of the hypersonic ramjet engine, wind tunnel test is done with freestream Mach number of 5.1. The total pressure recovery was 18 % with 12.9 % of ramp bleed. Several reasons for low total pressure recovery are shown. Supersonic internal compression is not enough. Then, the throat Mach number is high (M2.61) and total pressure losses at the terminal shock is large. Supersonic flow at the throat and position of the terminal shock is sensitive to a difference of the second ramp's throat height and the third ramp's throat height. Flow separations at the second ramp's trailing edge and the third ramp's leading edge are seen those could result in the trigger of unstart. The seal mechanism between the ramps and the sidewalls is important.

In this study, a fundamental experiment was carried out to investigate the frost formation on a cryogenic flat plate with/without temperature distribution from 230K to 160K under the convective flow. The effects of mixing ethanol as a condensable substance were also researched. From the test results, when surface temperature at the upstream is 230K, mass flux is high. On the other hand, when surface temperature at the downstream is 160K, mass flux is low. The degree of improvement to restrain frost formation by ethanol mixing is relatively larger at the upstream than at the downstream.

This paper reviews the latest studies of the expander cycle Air Turbo Ramjet engine (ATREX) conducted in JAXA. First, a system analysis including the vehicle and trajectory was conducted to optimize the engine cycle and turbo-machine configuration. We selected the precooled turbo-jet cycle for a prototype engine using the near term technologies. Second, a system ground-firing test was conducted to verify a defrosting system for the precooler. Methanol injection with its particles atomization could compensate 80 % of pressure loss caused by the frost. Thirdly, a feasibility of carbon/carbon composites for the engine components was investigated by making complex shapes such as a heat exchanger and a plug nozzle. Basic technologies on the gas leakage, the junction and bonding were also studied. The end of the paper, some basic studies such as wind tunnel tests of a new type air inlet and a plug nozzle are described.

The shock-induced combustion ramjet (shcramjet) is a hypersonic airbreathing propulsion concept which over-comes the drawbacks of the long, massive combustors present in the scramjet by using a standing oblique detonation wave (a coupled shock-combustion front) as a means of nearly instantaneous heat addition. A novel shcramjet combustor design that makes use of wedge-shaped flameholders to avoid detonation wave-wall interactions is proposed and analyzed with computational fluid dynamics (CFD) simulations in this study. The laminar, two-dimensional Navier-Stokes equations coupled with a non-equilibrium hydrogen-air combustion model based on chemical kinetics are used to represent the physical system. The equations are solved with the WARP (window-allocatable resolver for propulsion) CFD code (see: Parent, B. and Sislian, J. P., “The Use of Domain Decomposition in Accelerating the Convergence of Quasihyperbolic Systems”, J. of Comp. Physics, Vol. 179, No. 1,2002, pages 140-169). The solver was validated with experimental results found in the literature. A series of steady-state numerical simulations was conducted using WARP and it was deter-mined by means of thrust potential calculations that this combustor design is a viable one for shcramjet propulsion: assuming a shcramjet flight Mach number of twelve at an altitude of 36,000 m, the geometrical dimensions used for the combustor give rise to an operational range for combustor inlet Mach numbers between six and eight. Different shcramjet flight Mach numbers would require different combustor dimensions and hence a variable geometry system in or-der to be viable.

For the prediction of sloshing in the propellant tank of rocket vehicle utilized in RVT (reusable rocket vehicle testing) conducted by ISAS/JAXA, the flow field in the propellant tank during the ballistic flight was experimentally reproduced with the sub-scale model of it. The lateral acceleration as large as about 0.8 G was provided with a mechanical exciter and the deformation of liquid surface in the vessel was visualized with a high-speed camera. The several con-figurations of damping devices were installed and tested in the vessel, which should keep the ullage gas away from the outlet port. It was consequently suggested that the combination of a baffle plate and a perforated cylinder could be effective against the gas suction before the re-ignition of the engine. The sloshing phenomena were also simulated with the CFD code, called CIP-LSM. The numerical results showed good agreement with the corresponding data obtained in the experiment.

This paper describes development status and program of ATREX engine as a propulsion system of future spaceplane. Development activities using ATREX-500 engine from 1990 were finished in 2003 with large number of outcomes. We made system-level validation of the hydrogen fuel turbojet engine with air precooling device under sea level static condition. As a next step, we started design of the flight-type ATREX engine with large thrust and lightweight.

In this study, the affect of mounting axisymmetrical supersonic inlet to airfoil, which has 65 degree swept angle was numerically investigated. The parameter for this calculation are tree stream Mach number M=2.0 and 2.5, the distance between inlet spike and airfoil lower surface $L_{sw}$/$R_{cowl}$ = 1.21-1.54 and angle of attack to the airfoil 0-4. The mass capture ratio improved 3points in M=2.0 condition and 1points in M=2.5 while the mass capture ratio without airfoil surface was 57% and 71 % for each case. These are the result from increase of density and change of velocity deflection by the shock wave structure formed between inlet and airfoil surface. On the other hand, the distortion of Mach number at cowl lip plane increased by 13% in M=2.0, 3% in M=2.5 condition. The effects of the angle attack on the mass capture ratio is greater than that of the shock wave interaction between inlet and cowl, but the effects to the distortion is smaller in the range of this calculation condition. In the condition of M=2.0 with 4 degrees of angle of attack, inlet distortion of Mach number is mainly caused by the affects of the shock wave interaction between inlet and airfoil surface, while the largest angle of the velocity vector in the radial direction at cowl lip plane is caused by the affect of angle of attack. This large velocity vector made the flow inside the cowl subsonic and caused spillage, which interfere with the boundary layer of airfoil surface.

A variant of the magnetoplasma jet engine (magjet) is here proposed for airbreathing flight in the hypersonic regime. As shown in Figure 1, the engine consists of two distinct ducts: the high-speed duct, in which power is added electromagnetically to the incoming air by a momentum addition device, and the fuel cell duct in which the flow stagnation temperature is reduced by extracting energy through the use of a magnetoplas-madynamic (MPD) generator. The power generated is then used to accelerate the flow exiting the fuel cells with a fraction bypassed to the high-speed duct. The analysis is performed using a quasi one-dimensional model neglecting the Hall and ion slip effects, and fix-ing the fuel cell efficiency to 0.6. Results obtained show that the specific impulse of the magjet is at least equal to and up to 3 times the one of a turbojet, ram-jet, or scramjet in their respective flight Mach number range. Should the air stagnation temperature in the fuel cell compartment not exceed 5 times the incoming air static temperature, the maximal flight Mach number possible would vary between 6.5 and 15 for a magnitude of the ratio between the Joule heating and the work interaction in the MPD generator varied between 0.25 and 0.01, respectively. Increasing the mass flow rate ratio between the high speed and fuel cell ducts from 0.2 to 20 increases the engine efficiency by as much as 3 times in the lower supersonic range, while resulting in a less than 10% increase for a flight Mach number exceeding 8.

Variable air intake and variable exhaust nozzle of hypersonic engines are designed and tested in this study. Dimensions for variable geometry air intake, ram combustor and variable geometry exhaust nozzle are defined based on the requirements of a pre-cooled turbojet engine. Hypersonic Ramjet Engine is designed as a scaled test bed for each component. Actuation forces of moving parts for variable intake and variable nozzle are reduced by balancing the other force in the opposite direction. A demonstrator engine which includes variable intake and variable nozzle is designed and the components are fabricated. Composite material with silicone carbide is applied for high temperature parts under oxidation environment such as leading edge of the variable intake and combustor liner. Internal cooling structure is adopted for both moving and static parts of the variable nozzle. Pressure recovery and mass capture ratio of the variable intake at Mach 5 is obtained by a hypersonic wind tunnel test. Flow characteristics of the variable nozzle are obtained by a low temperature flow test. Wall temperature and heat flux of the nozzle at Mach 3 is obtained by a firing test. As results, the intake and the nozzle are proved to be used at designed pressure and temperature environment.

An experimental investigation was conducted to investigate the effects of the equivalence ratio and air mass flux on the combustion efficiency in a solid fuel ramjet used fuel grains which were highly loaded with boron carbide. Combustion efficiency increased with increasing equivalence ratio (grain length), and decreasing air mass flux. Higher inlet air temperature produced higher combustion efficiencies, apparently the result of enhanced combustion of the larger boron particles those bum in a diffusion controlled regime. Short grains which considered primarily of the recirculation region produced larger particles and lower combustion efficiencies. The result of the normalized combustion efficiency increased with inlet air temperature, is coincident with the result of the Brayton cycle thermal and the total efficiency relating to the heat input.

Experiments and numerical calculations were conducted on the flow field of a model ejector ramjet configuration to investigate fundamental fluid dynamic aspects of its pumping and mixing effects. Also a one-dimensional flight performance analysis program was constructed with a simple ejector modeling. After comparing the model with some of the previous experimental and numerical results, a flight performance analysis was conducted with the program. The present states of the program and some features to be improved are presented.

Boundary layer ingestion in airframe-integrated scramjet engines causes engine stall (“engine un start” hereafter) and restricts engine performance. To improve the unstart characteristics in engines, boundary layer bleed and a two-staged injection of fuel were examined in Mach 4 and Mach 6 engine tests. A boundary layer bleed system consisting of a porous plate, an air coolers, a metering orifice and an ON/OFF valve, was designed for each of the engines. First, a method to determine bleed rate requirements was developed. Porous plates were designed to suck air out of the Mach 4 engine at a rate of 200 g/s and out of the Mach 6 engine at a rate of 30 g/s. Air coolers were then optimized based on the bleed airflow rates. The exhaust air temperature could be cooled below 600 K in the porous plates and the compact air coolers. The Mach 4 engine tests showed that a small bleed rate of 3% doubled the engine operating range and thrust. With the assistance of two-staged fuel injection of H2, the engine operating range was extended to Ф0.95 and the maximum thrust was tripled to 2560 N. The Mach 6 tests showed that a bleed of 30 g/s (0.6% of captured air in the engine) extended the start limit from Ф0.48 to Ф1 to deliver a maximum thrust of 2460 N.

In this paper are presented TSTO system analysis including some controlled variables on the engine operation such as a fuel flow rate and a pressure ratio of compressor, as well as variables on the trajectory. TSTO studied here is accelerated up to Mach 6 by a fly-back booster powered by air breathing engines. Three different types of engine cycle were treated for propulsion system of the booster, such as a turbo ramjet, a precooled turbojet and an EXpander cycle Air Turbo Ramjet (ATREX). The history of the controlled variables on the engine operation was optimized by Sequential Quadratic Programming (SQP) to accomplish the minimum fuel consumption. The trajectory was also optimized simultaneously. The results showed that the turbo ramjet gave the best fuel consumption. The optimal trajectory was almost the same except in the transonic range and just before reaching to Mach 6. The history of the pressure ratio of compressor considerably depended on the engine type. It is concluded that simultaneous optimization for engine control and trajectory is effective especially for a high-speed airplane propelled by turbojets like the TSTO booster.

In a liquid-fueled ramjet engine, the insufficient mixing and evaporation result in the low combustion efficiency and combustion instability. Improving its spray characteristics and devising a means of mixing fuel droplets with air may compensate these disadvantages of liquid fuel ramjet engine. The jet penetrations of various fuel injectors were measured to investigate the spray characteristics of a liquid-fueled ramjet engine under high pressure air-stream conditions. The penetrations in high pressure conditions are smaller than the values calculated from Inamura's or Lee's equations, and, in the high pressure conditions, the jet penetrations are similar each other. In the dual hole injectors, the jet penetrations of rear orifice is rapidly increased due to the reduction of the drag, which is created by the jet column of front orifice. The jet penetration of rear orifice is increased because of the drag reduction created by the jet column of the front orifice. And, because of the drag reduction formed by the column of jet, the jet penetration in the rear orifice of dual hole injector is much larger than the jet penetration of single hole injector. As the distances of the orifice are increased, the jet penetrations of the rear orifice decrease.

Conceptual studies of a combined-cycle engine have been conducted. Herein, the results are presented. The engine is composed of ejector-jet, ramjet, scramjet and rocket modes, and will be mounted on the Single-Stage-to-Orbit aerospace plane. Propellants are hydrogen and oxygen. Calculated engine thrust performances and cooling requirement of the engine are presented. Pitching moment of the plane with the engine will be balanced even in the vacuum condition. The experimental results of the inlet and the ejector-jet, ramjet and scramjet modes are presented. The effect of the airframe configuration on the engine performance and the thermal environment in the in-side of the plane are also presented. Through the investigations, possibilities of the combined-cycle engine and the aerospace plane are being made clear now.

A possible and practical engine system research method is proposed. Varieties of objectives of the engine component and system technology developments are fulfilled by the small scale rig and engine demonstration. Some research applications of small jet engines in National Aerospace Laboratory of Japan (NAL) are presented together with historical overview.

When a gas expands through a convergent nozzle in which the ratio of the ambient to the stagnation pressures is higher than that of the critical one, the issuing jet from the nozzle is underexpanded. If a flat plate is placed normal to the jet at a certain distance from the nozzle, a detached shock wave is formed at a region between the nozzle exit and the plate. In general, supersonic moist air jet technologies with nonequilibrium condensation are very often applied to industrial manufacturing processes. In spite of the importance in major characteristics of the supersonic moist air jets impinging to a solid body, its qualitative characteristics can not even know. In the present study, the effect of the nonequilibrium condensation on the underexpanded moist air jet impinging on a vertical flat plate is investigated experimentally. Flow visualization and impact pressure measurement are performed for various relative humidities and flat plate positions. The obtained results show the plate shock and Mach disk are dependent on the nozzle pressure ratio and the relative humidity, but for a given nozzle pressure ratio, the diameters of the plate shock and Mach disk depend on the stagnation relative humidity. The impact pressure deviation from the flow of without condensation is large, as the relative stagnation humidity increases.

The atomization process of a circular $SF_{6}$ liquid jet issued into an otherwise quiescent, high-pressure $N_2$ gas was observed to explore the breakup mechanism of liquid ligaments involved in turbulent atomization. Both liquid and gas temperatures were fixed at a room temperature but the gas pressure was elevated to more than twice the critical pressure of $SF_{6}$. Therefore, the liquid surface was in a thermodynamic state close to a critical mixing condition with suppressed vaporization. Since the surface tension and the surface gas density approach zero and the surface liquid density, respectively, phenomena equivalent to those which would appear when a very high speed laminar flow of water were injected into the atmospheric-pressure air can be observed by issuing $SF_{6}$ liquid at low speeds in micro-gravity environment which avoid disturbances due to gravity forces. The instability ob near-critical mixing surface jet was quantitatively characterized using a newly developed device, which could issue a very small amount of $SF_{6}$ liquid at small constant velocity into a very high-pressure $N_2$ gas.

First this paper introduces an advanced FADEC (Full Authority Digital Electric Control) for current and future jet engines.It is designed to realize not only stable thrust control, but also performance improvement, reliability enhancement, service life extension, etc. It can be built by using current micro-processor with high computational power and there exists no difficulties but reliability problem of the micro- processor. Next, the simulation results of SFC minimization control are shown. The target engine is a supersonic, low-bypass ratio, 2-spool, combined cycle turbofan, designated as HYPR90T, which consists of a turbo engine for under Mach 3 flight and a ram engine for over Mach 3 flight. he results can then be used for performance optimization of the engine, which plays important role in the advanced FADEC.

A visualization study of shock formation of the supersonic jet nozzle using a Shadowgraph Method (SM) was carried out to investigate the effect of the longitudinal variation of coaxial pipe end tip position inside the supersonic nozzle. The experiment was performed for the Mach number range from 1.1 to 1.2 at nozzle exit. The well known shock cell structure was shown with the pipe end located deep inside the nozzle for the studied Mach number. With the pipe end approaches nozzle exit, it was found that the shock cell structure disappeared and turned into complex formation. In order to understand the mechanism of the shock structural change, computational simulation was carried out using the Navier-Stokes solver, FLUENT. Topological sketch was added with an aid of the visualization and the numerical simulation.

The pressure sensitive paint (PSP) technique has been well established in external flow field. However, there are still unresolved issues in internal flow field. This work was focused on the application to unsteady pressure measurement of fan flow field. The PSP measurement system was established and the image processing software was developed. First, the performance of PSP was investigated at the static cell. Then the unsteady PSP measurement was carried out at fan test facility. PSP data images were acquired from the suction and pressure surface of stator vanes. Pressure distributions on the surface of the stator vane were detected non-intrusively. The issues of image acquisition and image processing were clarified through the practical PSP application to fan flow field.

An experimental study has been carried out to examine heat-transfer characteristics of an axisymmetric, under-expanded, sonic jet impinging on a flat plate and the local measurement of surface pressures and heat transfer coefficients on a plate have been achieved together with a visualization test of shock structure in a jet. As a result, it has been found that the Nusselt number distribution has different aspects depending on the under-expansion ratios and the nozzle-to-plate distances.

For the protection of the local air quality and the global atmosphere, the emissions of trace species including nitric oxides (NO and NO$_2$) from gas turbines are regulated by local governments and by the International Civil Aviation Organization. In-situ measurements of such species are needed not only for the development of advanced low-emission combustion concepts but also for providing emissions data required for the sound assessment of the effects of the emissions on environment. We have been developing a laser absorption system that has a capability of simultaneous determination of NO and NO$_2$concentrations in the exhaust jets from aero gas turbines. A diode laser operating near 1.8 micrometer is used for the detection of NO while a separated visible tunable diode laser operating near 676 nanometers is used for NO$_2$. The sensitivities at elevated temperature conditions were determined for simulated gas mixtures heated up to 500K in a heated cell of a straight 0.5 m optical path. Sensitivity limits estimated as were 30 ppmv-m and 3.7 ppmv-m for NO and NO$_2$, respectively, at a typical exhaust gas temperature of 800K. Experiments using the simulated exhaust flows have proven that $CO_2$ and $H_2O$ vapor - both major combustion products - do not show any interference in the NO or NO$_2$ measurements. The measurement system has been applied to the NO/NO$_2$ measurements in NO and NO$_2$ doped real combustion gas jets issuing from a rectangular nozzle having 0.4 m optical path. The lower detection limits of the system were considerably decreased by using a multipass optical cell. A pair of off-axis parabola mirrors successfully suppressed the beam steering in the combustion gas jets by centralizing the fluctuating beam in sensor area of the detectors.

In this paper, a small supersonic wind tunnel was designed and built to study the flow characteristics of a supersonic impulse turbine cascade by experiment. The flow was visualized by means of a single pass Schlieren system. The supersonic cascade with 3-dimensional supersonic nozzle was tested over a wide range of pressure ratio. Highly complicated flow patterns including shocks, nozzle-cascade interaction and shock boundary layer interactions were observed.

The paper studies the effect of neighboring blade rows on flutter characteristics of cascading blades. For this purpose the computation program to calculate the unsteady blade loading based on the un-steady lifting surface theory for contra-rotating annular cascades was formulated and coded. Then a computation program to solve the coupled bending-torsion flutter equation for the contra-rotating annular cascades was also developed. Some results of the flutter analysis are presented. The presence of the neighboring blade row gives rise to significant change in the critical flutter condition when the main acoustic duct mode is of cut-on state.

A numerical analysis based on two-dimensional and three-dimensional incompressible Navier-Stokes equations has been carried out for double-circular-arc (DCA) compressor cascades. Two types of double-circular-arc cascades were used in this analysis. The appropriate turbulence model for compressor analysis was selected among the conventional turbulence models such as Baldwin-Lomax, k-$\varepsilon$ and k-$\varepsilon$ models. The results of current study were compared with available experimental data at various incidence angles. The 2-D and 3-D computational codes based on SIMPLE/PWIM algorithm for collocated grid and hybrid scheme for the convective terms were the main features of numerical tools. As commonly known, turbulence modeling is very important for the prediction of cascade flows, which are extremely complex with separation and reattachment by adverse pressure gradient. For selection of turbulence model, 2-D analysis was performed. And then, k-$\varepsilon$ turbulence model with wall function chosen as the reasonable turbulence model for 3-D calculation was used to increase the efficiency of computation times. A reasonable result of 3-D flow pattern passing through the double-circular-arc cascade was obtained.

A fluid-solid conjugate solver has been newly developed and applied to an actual engine disk system. Most of the currently available conjugate solvers lack the special thermal modeling for turbomachinery disk system applications. In the present new code, these special models are implemented to expand the applicability of the conjugate method and to reduce the required computational resources. Most of the conjugate analysis work so far are limited to the axisymmetric framework. However, the actual disk system includes several non-axisymmetric components which inevitably affect the local heat transfer phenomena. Also the previous work devoted to this area usually concentrate their efforts on the steady-state thermal field, although the one in the transient condition is more critical to the engine components. This paper presents full 3D conjugate analysis of a single stage high pressure turbine rotor-stator disk system to assess the three-dimensional effects (Fig. 1). The analysis is carried out not only in the steady-state but also in the engine accelerating transient condition. The predicted temperatures shows good agreement with measured data.

In this study, the flow characteristics within supersonic cascades are numerically investigated by using Fine Turbo, a commercial CFD code. Cascade flows are computed for three different inlet conditions. : a uniform supersonic inlet condition, a linear nozzle and a converging-diverging nozzle located in front of cascades. The effect of inlet conditions is compared and flow characteristics including shock patterns and shock-boundary layer interaction are analyzed. Also the effect of design parameters such as pitch-chord ratio, blade angle and blade surface curvature on the flow within supersonic cascades are studied.

Regarding to the project SUAV (Smart Unmanned Aerial Vehicle) in KARI (Korea Aerospace Research Institute), several engine configurations has been evaluated. However it's not an easy task to collect all the necessary data of each engine for the analysis. Usually, some kind of modeling technique is required in order to determine the unknown data. In the present paper a qualitative method for reverse engineering is proposed, in order to identify some design patterns and relationships between parameters. The method can be used to estimate several parameters that usually are not provided by the manufacturer. The method consists of modeling an existing engine and through a simulation, compare its transient behavior with its operating envelope. In the simulation several parameters such as thermodynamics, performance, safety and mechanics concerning to the definition of operation-envelope, have been discussed qualitatively. With the model, all engine parameters can be estimated with acceptable accuracy, making possible the study of dependencies among different parameters such as power-turbine total inertia, TIT, take-off time and part load, in order to check if the engine transient performance is within the design criteria. For more realistic approach and more detailed design requirements, it will be necessary to enhance the compressor map first, and more realistic estimated values must be taken into account for intake-loss, bleed-air and auxiliary power extraction. The relative importance of these “unknown” parameters must be evaluated using sensitivity analysis in the future evaluation. Moreover, fluid dynamics, thermal analysis and stress analysis necessary for the resulting life assessment of en engine, will not be addressed here but in a future paper. With the methodology presented in the paper was possible to infer the relationships between operation-envelope and engine parameters.

"Research and Development of Melt-Growth Composite (MGC) Ultra High Efficiency Gas Turbine System Technology" program has been started in JFY2001. The main objective of the program is to establish basic component technologies to apply MGC material to an efficient gas turbine system successfully. It is known that MGC material maintains its mechanical strength at room temperature up to about 2000 K, which is ideal for the high temperature gas turbine. The purposes of the present study are to develop the cooling structure of the gas turbine combustor liner where MGC material is applied as the heat shield panel, also to develop the low NOx combustion system for a 1970 K (1700 deg.C) class gas turbine combustor. To start with, basic heat transfer characteristics were investigated by one-dimensional calculation and heat transfer experiment for the cooling structure. Axially staged configuration and fuel preparation were investigated by CFD calculation and experiments for the low NOx combustor.

A numerical flow analysis has been performed on the partial admission turbine of KARI turbopump to support the aerodynamic and structural dynamic assessments. The flow-field in a partial admission turbine is essentially three dimensional and unsteady because of a tip clearance and a finite number of nozzles. Therefore the mixing plane method is generally not appropriate. To avoid heavy computational load due to an unsteady three dimensional calculation, a frozen rotor method was implemented in steady calculation. It adopted a rotating frame in the grid block of a rotor blade by adding some source terms in governing equations. Its results were compared with a mixing plane method. The frozen rotor method can detect the variation of flow-field dependent upon the blade's circumferential position relative to the nozzle. It gives a idea of wake loss mechanism starting from the lip of a nozzle. This wake loss was assumed to be one of the most difficult issues in turbine designers. Thus, the frozen rotor approach has proven to be an efficient and robust tool in design of a partial admission turbine.

A smart blade conception has been proposed by the authors. With stretching-twisting coupling effect, the blade is twisted by centrifugal load or ambient temperature change. In this paper, the blades, made by three kinds of anti-symmetric laminates, are investigated by rotational tests. The results show the angle of smart blade tips increases in proportion to the 2nd power of a rotating speed and is well in agreement with the numerical results by FEM.

In Japan, a long-waited civil aero engine development project has been recently started by the New Energy and Industrial Technology Development Organization, NEDO. “High efficiency,” “environmentally friendliness” and “low-cost” are the key words of the target engine. The target engine is of l0000-lb thrust with project consists of three phases: Feasibility studies and market research in the first phase, FY 2003, engine component development in the second phase, FY 2004-2006, and core and full engine demonstrators in the third phase, FY 2007-2009. In league with this government/industry joint funded project, Institute of Space Technology and Aeronautics, JAXA, has initiated “Clean Engine Technology” project.