Cooling Effect of Water Injection on a High-Temperature Supersonic Jet

The high temperature and high pressure supersonic jet is one of the key problems in the design of solid rocket motors. To reduce the jet temperature and noise, cooling water is typically injected into the exhaust plume. Numerical simulations for the gas-liquid multiphase flow field with mixture multiphase model were developed and a series of experiments were carried out. By introducing the energy source terms caused by the vaporization of liquid water into the energy equation, a coupling solution was developed to calculate the multiphase flow field. The temperature data predictions agreed well with the experimental results. When water was injected into the plume, the high temperature core region area was reduced, and the temperature on the head face was much lower than that without water. The relationship between the reduction of temperature on the bottom plate and the momentum ratio is developed, which can be used to predict the cooling effect of water injection in many cases.


Introduction
Supersonic high-temperature exhaust plumes cause rocket engine damage on the launcher and surrounding equipment.Thermal erosion due to high temperature is one of the main factors damaging the launcher.Generally, the launcher, the diversion trench and even the missile are directly affected by jet exhaust.Usually water is employed as the coolant for blast deflectors [1] owing to its low-cost and ready availability.The injection of water into the jet exhaust as a solution to prolong the launcher life has been suggested.This solution can not only lower the temperature, but also suppress the jet noise during launch [2] and nitrous oxide emissions [3].
On account of its significant latent heat of vaporization and higher heat capacity than air, water injection represents an effective approach to reduce the plume temperature.In 1989, Mathias [4] and co-workers developed a one-dimensional thermal model to confirm the cooling effect of water spray on the bottom surface of the steel case in the redesigned system.In 1990s, Miller and Koo [5] studied the plume flow field by injecting water.The results showed that the loss of heat flux caused by the vaporization and the momentum of water flow were the main factors that affected the exhaust jet.Through a series of subscale experiments the cooling effect of water had been quantified in the gas and an equation to validate the relationship between reduction in heat flux and momentum ratio was proposed.The research laid a foundation for further study of the effect of water to the rocket exhaust environment.Santangelo and Kennedy [6] found that the jet structure was retained with and without the spray from the Mie scattering images and no gross changes happened in gas phase structure from Energies 2015, 8, 13194-13210; doi:10.3390/en81112363www.mdpi.com/journal/energiesplanar laser induced fluorescence (PLIF) measurements.Golden [3] conducted a series of experiments with water fog injection in a gas turbine and measured the temperature of the combustor to study the reduction of nitrous oxide emissions.Engblom [7] validated a new conjugate heat-transfer capability against thermal experimental data for three water-cooled supersonic flow experiments.Kandula [8] developed a one-dimensional control volume formulation to predict turbulent jet mixing noise reduction of the exhaust plume by water injection and the mechanism of interaction of water injection and gas.Cho [9] investigated the effect of coolant injection to small-scale diffuser, and the results indicated that the ratio between gas and water mass flow rates is an important parameter in determining the cooling characteristics.Krothapalli and his co-workers [10][11][12] studied the effect of different size droplets.A series of experiments were performed to confirm the favorable effect on the atomization by using water sprays with micro-nozzles.The results showed that the effect of water injection was to reduce the magnitude of the source terms.They also measured the effect of water injection on the turbulence intensity of jet and thickness of shear layer.
In the past 50 years, liquid injection to hypersonic flow was frequently put forward in studying supersonic combustion ramjets.The gas-liquid multiphase field of rocket exhaust plumes with water injection was studied by subscale testing and numerical simulation.For the studies on jet gas, accuracy of numerical simulation and reliability of testing at high Mach numbers is required.
In this work, the mixture multiphase model was used to simulate gas-liquid multiphase field in a coupling way by inducing the energy source items that were caused by the vaporization of liquid water into the energy equation.The vaporization energy increments were added to the energy equation, and the vaporization mechanism was expressed in vaporization equation to investigate the liquid water vaporization process.From photographs using high speed cameras and infrared thermography, the characteristics of the flow field can be seen clearly.At the same time, experimental data served to validate the accuracy of the numerical simulation.Lastly, the results of numerical simulation were compared with the experimental results.The cooling effect of water injection was quantified.

Multiphase Flow Model
There are two approaches for the numerical calculation of multiphase flow: the Euler-Lagrange approach and the Euler-Euler approach.A fundamental assumption made in Euler-Lagrange mode is that the dispersed second phase occupies a low volume fraction, even though high mass loading is acceptable.The fluid phase is treated as a continuous phase while the motion of the discontinuous particulate phase is obtained by integrating the equation of motion of individual particles along their trajectories.In the Euler-Euler approach, the different phases are modeled as interpenetrating continua, with the probability of existence of each phase at a point in the computational domain specified by its respective volume fraction.Conservation equations for each phase are derived to obtain a set of equations, which have similar structure for all phases.These equations are closed by providing constitutive relations that are obtained from empirical information or kinetic theory.
In the present problem, the volume fraction of the second phase (water) is modeled as a continuous phase and cannot be neglected.Thus, Euler-Euler approach has been chosen for this case.Three different Euler-Euler multiphase models are provided by Fluent commercial code [13]: the Euler model, the mixture model and volume of fluid (VOF) model.In some cases, it is possible to use a simplified formulation of the Euler model.
The mixture model (or "algebraic slip model") represents a simplified formulation of the Euler model, especially when the interphase drag laws are unknown or their applicability to the system is questionable.The mixture model may be applicable, e.g., for a relatively homogeneous suspension of one or more species of dispersed phase that closely follow the motion of the continuous carrier fluid.In 1966, Wen and Yu [14] created the Wen-Yu model to describe momentum exchange coefficient and Energies 2015, 8, 13194-13210 interphase drag force equation in mixture model.Based on this model, a complete mixture model was built by Bowen [15].In 1980, Lee [16] derived detailed and complete vaporization equation and condensation equation, and combined them with the mixture model to simulate gas-liquid flow.Yamazaki [17] analyzed the effect factors about the diameter of liquid drop in gas-liquid by applying mixture model.In 2010, Ma and Jiang [18] used mixture model to analyze the cooling effect of a wet concentric canister launcher.
In this paper, the gas flow field with water injection is calculated with the usage of the MIXTURE multiphase model.In the calculations, the gas phase flow and liquid phase flow are full fused in every cell.The physical parameters of the mixture are analyzed by using the weighted average method.With the mixture model we solve the continuity, momentum and the energy equations for the mixture, and the volume fraction equation for the secondary phases, as well as algebraic expressions for the relative velocities [13].

Mass conservation equation
The continuity equation for the mixture is: where Ñ ν m is the mass-averaged velocity of the mixture: and ρ m is the mixture density: α k is the volume fraction of phase k.

Momentum Equation
The momentum equation for the mixture can be obtained by summing the individual momentum equations for all phases.It can be expressed as: where n is the number of phases; Ñ F is the body force; µ m is the viscosity of the mixture: Ñ ν dr,k is the drift velocity of secondary phase k:

Energy Equation
The energy equation for the mixture is defined as: where the effective conductivity k e f f is defined as Equation (8), and the first term on the right-hand side of Equation ( 7) represents energy transfer due to heat conduction.S E includes any other volumetric heat sources: where k t is the turbulent thermal conductivity.
In Equation (7), for a compressible phase, E k is defined as: where h k is the sensible enthalpy for phase k.

Volume Fraction Equation for the Secondary Phases
From the continuity equation for secondary phase, the volume fraction equation for the secondary phase can be obtained: where .m 12 is the mass transfer from main phase to secondary phase.
. m 21 is the mass transfer from secondary phase to main phase.Ñ ν dr,2 is the drift velocity of secondary phase.Ñ ν 2 is the velocity of secondary phase.Ñ ν m is the mass-averaged velocity of the mixture and is shown in Equation ( 2): Thus the volume fraction of main phase α 1 is calculated from Equations ( 10) and (12).

Vaporization and Condensation Equation
The vaporization and condensation of water has fundamental effects on the flow field by the interaction of liquid water and the high temperature combustion gas.The vapor-liquid equilibrium compositions may exhibit a nonlinear behavior.To include the influence of the vaporization process, energy and mass source terms, caused by the vaporization of liquid water, were introduced to the conservation equations.The vaporization and condensation processes are calculated in every finite cell.The vaporization ratio of water is calculated according to its local temperature and pressure.When the temperature of the mixture is higher than the saturation temperature of water, the water transforms into vapor by absorbing heat; otherwise, the vapor condenses to water by releasing heat (each local pressure has a corresponding saturation temperature).The transformation mechanism of vapor phase and liquid phase in each cell is expressed by the following relations: Water vaporization: .
where .m l and .m ν are vaporization rate of liquid and condensation rate of vapor, respectively, and they are reported in units of kilogram at one cubic meter per second.λ is the time relaxation factor taken as 0.1 in this study, T sat is saturated temperature of water, T l is the transient temperature of liquid, T ν is the transient temperature of gas, and α l and α ν are the volume fraction of liquid and gas.
The transient vaporization rate of water at time t for one cell is calculated by the above formulas: . m " .
It should be pointed out that the change in energy is caused by the phase transition: where S k is energy source term, which is one of S E , and ∆H is latent heat of vaporization of saturated water.

Computational Model
Figure 1 shows the computational domain and boundary conditions.Four temperature monitoring points on the bottom plate are marked.The distances of those points to the symmetry axis are 0.2, 0.3, 0.4 and 0.5 m, respectively.The angle between the water pipe and the symmetry axis is 60 ˝.Because of the symmetry of jet flow, a quarter of three-dimensional computational domain is used to solve the flow field.This model is suitable for both two-pipe case and four-pipe case.

Computational Model
Figure 1 shows the computational domain and boundary conditions.Four temperature monitoring points on the bottom plate are marked.The distances of those points to the symmetry axis are 0.2, 0.3, 0.4 and 0.5 m, respectively.The angle between the water pipe and the symmetry axis is 60°.Because of the symmetry of jet flow, a quarter of three-dimensional computational domain is used to solve the flow field.This model is suitable for both two-pipe case and four-pipe case.The nozzle wall, water pipes' walls, and bottom plate are set as solid walls.The inlet boundary type of the nozzle is pressure inlet, and chamber pressure is 7 MPa, and temperature is 3000 K.The outlet boundary type of the plume is pressure outlet boundary, pressure is 1 bar, and temperature is 300 K.The mass flow rate of the rocket motor Qg is 1.5 kg/s, which is invariable.The velocity of the combustion gas at the exit plate of the nozzle Vg is invariable, which is 2200 m/s.The outlet water velocity, number of pipes, and total mass flow rate of water pipes are variable in different cases.

Grid Independence
In transient numerical simulation, grid size and density directly have an effect on the accuracy of numerical calculations.Verifying the independence of grid is necessary before calculation.In this study, eight kinds of models with different grid sizes and densities are simulated.These models are divided into 100 thousand, 350 thousand, 600 thousand, 850 thousand, 1100 thousand, 1600 thousand, 2100 thousand and 2600 thousand grids, respectively.In Figure 2, there are 600 thousand grids in Model A, 1100 thousand grids in model B, 2600 thousand grids in model C. By refining the grids in gas flow zone and zone nearby water pipe, the simulation results are close to the experimental results.
The two-phase flow coupled interaction is closely related to temperature in this study.By comparing the temperature, the appropriate grid size is determined.In Figure 3, the dimensionless temperature difference ratio is : where is the temperature at the center of the bottom plate in model C and is the temperature at the center of the bottom plate in other models.The difference ratio decreases with increasing grid number and then tends to be slow.The ratio is within 3% in the 1100 thousand, 1600 thousand and 2100 thousand grids models.The nozzle wall, water pipes' walls, and bottom plate are set as solid walls.The inlet boundary type of the nozzle is pressure inlet, and chamber pressure is 7 MPa, and temperature is 3000 K.The outlet boundary type of the plume is pressure outlet boundary, pressure is 1 bar, and temperature is 300 K.The mass flow rate of the rocket motor Q g is 1.5 kg/s, which is invariable.The velocity of the combustion gas at the exit plate of the nozzle V g is invariable, which is 2200 m/s.The outlet water velocity, number of pipes, and total mass flow rate of water pipes are variable in different cases.

Grid Independence
In transient numerical simulation, grid size and density directly have an effect on the accuracy of numerical calculations.Verifying the independence of grid is necessary before calculation.In this study, eight kinds of models with different grid sizes and densities are simulated.These models are divided into 100 thousand, 350 thousand, 600 thousand, 850 thousand, 1100 thousand, 1600 thousand, 2100 thousand and 2600 thousand grids, respectively.In Figure 2, there are 600 thousand grids in Model A, 1100 thousand grids in model B, 2600 thousand grids in model C. By refining the grids in gas flow zone and zone nearby water pipe, the simulation results are close to the experimental results.
The two-phase flow coupled interaction is closely related to temperature in this study.By comparing the temperature, the appropriate grid size is determined.In Figure 3, the dimensionless temperature difference ratio is R: R " |pT C ´Tother q {T C | (17 where T C is the temperature at the center of the bottom plate in model C and T other is the temperature at the center of the bottom plate in other models.The difference ratio decreases with increasing grid number and then tends to be slow.The ratio is within 3% in the 1100 thousand, 1600 thousand and 2100 thousand grids models.On account of the different grid sizes in these models, the temperature differences are compared in Figure 3.The Peclet number can be taken to explain the grid partitioning.It is a dimensionless number.It is defined to be the ratio of the rate of convection of a physical quantity by the flow to the rate of diffusion of the same quantity driven by an appropriate gradient.In 1976, Roache put forward the non-dimensional cell Peclet number as a measure of the relative strengths of convection and diffusion (Equation ( 18)) [19]:  On account of the different grid sizes in these models, the temperature differences are compared in Figure 3.The Peclet number can be taken to explain the grid partitioning.It is a dimensionless number.It is defined to be the ratio of the rate of convection of a physical quantity by the flow to the rate of diffusion of the same quantity driven by an appropriate gradient.In 1976, Roache put forward the non-dimensional cell Peclet number as a measure of the relative strengths of On account of the different grid sizes in these models, the temperature differences are compared in Figure 3.The Peclet number can be taken to explain the grid partitioning.It is a dimensionless number.It is defined to be the ratio of the rate of convection of a physical quantity by the flow to the rate of diffusion of the same quantity driven by an appropriate gradient.In 1976, Roache put forward the non-dimensional cell Peclet number as a measure of the relative strengths of convection and diffusion (Equation ( 18)) [19]: Pe " ρu{ pΓ{δxq (18) where u is the local flow velocity and δx is the cell width.ρ is the density.Γ is a constant relevant to flow field physical property on the compute node.Pe is proportional to δx.When Pe is greater than 1, convection is superior to diffusion and the calculation result is greatly influenced by the upstream node.Otherwise, diffusion is superior to convection and the calculation result is greatly influenced by the downstream node.
In this numerical calculation, u nearby the bottom plate is very small.With small δx the diffusion-dominated flow has a more accurate numerical result.The errors by numerical calculation in models with small grid size and refined grids are less than that in other models.
Also, the computing efficiency in the 1600 thousand, 2100 thousand and 2600 thousand grids models is about two-thirds of that in model B, so the 1100 thousand grid model was chosen for numerical calculation.

Numerical Method
The finite volume method was used to discretize the governing equations.The RNG k ´ε model was used, as it is more accurate and reliable for a wider class of flows than the standard k ´ε model.The adiabatic standard wall function was chosen for the working time of the engine.With coupling calculation, the energy source items were applied to the mixture model by UDF (User Defined Function) in gas-liquid multiphase flow filed.

Experiment System
The schematic of the experimental system is shown in Figure 4. Water was supplied by a water pump to ensure constant water pressure.The mass flow rate was measured by a flow meter.
Energies 2015, 8, page-page where is the local flow velocity and δ is the cell width.ρ is the density.Γ is a constant relevant to flow field physical property on the compute node.
is proportional to δ .When is greater than 1, convection is superior to diffusion and the calculation result is greatly influenced by the upstream node.Otherwise, diffusion is superior to convection and the calculation result is greatly influenced by the downstream node.
In this numerical calculation, nearby the bottom plate is very small.With small δ the diffusion-dominated flow has a more accurate numerical result.The errors by numerical calculation in models with small grid size and refined grids are less than that in other models.
Also, the computing efficiency in the 1600 thousand, 2100 thousand and 2600 thousand grids models is about two-thirds of that in model B, so the 1100 thousand grid model was chosen for numerical calculation.

Numerical Method
The finite volume method was used to discretize the governing equations.The RNG k − ε model was used, as it is more accurate and reliable for a wider class of flows than the standard k − ε model.The adiabatic standard wall function was chosen for the working time of the engine.With coupling calculation, the energy source items were applied to the mixture model by UDF (User Defined Function) in gas-liquid multiphase flow filed.

Experiment System
The schematic of the experimental system is shown in Figure 4. Water was supplied by a water pump to ensure constant water pressure.The mass flow rate was measured by a flow meter.The angle between the water pipes (two or four pipes were chosen in different cases) and the centerline axis is 60°.The distance between the point of intersection of the centers of two water-injection pipes and that of nozzle exit is 520 mm.The temperature of the water injection is 300 K and water flow velocity is 16 m/s.The diameter of the round pipe is constant and equal to 12 mm.Four temperature sensors were fixed to the bottom plate of the bench to monitor the Energies 2015, 8,[13194][13195][13196][13197][13198][13199][13200][13201][13202][13203][13204][13205][13206][13207][13208][13209][13210] The angle between the water pipes (two or four pipes were chosen in different cases) and the centerline axis is 60 ˝.The distance between the point of intersection of the centers of two water-injection pipes and that of nozzle exit is 520 mm.The temperature of the water injection is 300 K and water flow velocity is 16 m/s.The diameter of the round pipe is constant and equal to 12 mm.Four temperature sensors were fixed to the bottom plate of the bench to monitor the temperature.The locations of these sections are shown in Figure 4.An XSJ flow totalizer (Shanghai MeiKai YouDi Instrument Co., Shanghai, China) was employed for flow measurements.In addition to this contact measurement, non-contact measurement was chosen to detect the characters of the multiphase flow field by high-speed camera and infrared thermograph.To capture clearer images from the highly unsteady interaction between the water and combustion gas, a Phantom v10 high-speed camera and TVS2000 (NEC-AVIO Co., Tokyo, Japan) infrared thermograph were used to observe the flow field structure and the radiation intensity emitted from the exhaust plume.The sampling frequency of the Phantom v10 high-speed camera is 480 images/s and the resolution of the images is 2400 ˆ1800 pixels.The resolution and the sampling frequency of the TVS2000 infrared thermograph is 0.1 ˝C and 30 images/s.The E12-1-C-U temperature sensors (Nanmac Co., Washington, DC, USA) are used to monitor the temperature of the four points.The measuring limit of E12-1-C-U temperature sensor is 3300 K, which could satisfy the demands of exhaust plumes temperature test.
A horizontal motor test bench is generally selected for the engine test, as it is easy to design and assures testing security.However, in consideration of the large mass of water, which may result in much more gravitational influence, a vertical test bench is chosen to simulate the rocket launching more truly.Furthermore, a vertical motor test bench makes the measurements more accurate and the test equipment such as high-speed camera and infrared thermograph is more easily located.
Too high an experiment mounting plate tends to render the entire device unsteady.Too low an experiment mounting plate makes it difficult to perform a complete observation.In this study, the height of the mounting plate depends on the length of the core region of the jet: x r e " 6.9 p1 `0.38Maq 2 (19) where x is the length of the core region; r e is the radius of nozzle exit; Ma represents Mach number and equals 3.5 in this study.Calculated according to this formula [20], the length of the core region is 1.33 m.The simulation results show that the length of the core region is between 1.25 and 1.4 m.The simulation confirms the effectiveness of Equation (19).Then, considering the length of the motor nozzle, the height of the mounting plate is 1.9 m.
To simulate typical rocket exhaust conditions, the designed Mach number of the free jet in the test is about 3.5 and the outlet temperature is approximately 2000 K. To realize this condition, a solid rocket motor was chosen.The total temperature and total pressure of the combustion chamber are 3000 K and 7 MPa, respectively.The jet flow of this rocket motor is an axisymmetric over-expanded flow, so the rocket motor is designed based on the basic formula of one-dimensional isentropic flow and the combustion characteristic of solid motor grain.The rocket motor is shown in Figure 5.
Energies 2015, 8, page-page multiphase flow field by high-speed camera and infrared thermograph.To capture clearer images from the highly unsteady interaction between the water and combustion gas, a Phantom v10 high-speed camera and TVS2000 (NEC-AVIO Co., Tokyo, Japan) infrared thermograph were used to observe the flow field structure and the radiation intensity emitted from the exhaust plume.The sampling frequency of the Phantom v10 high-speed camera is 480 images/s and the resolution of the images is 2400 × 1800 pixels.The resolution and the sampling frequency of the TVS2000 infrared thermograph is 0.1 °C and 30 images/s.The E12-1-C-U temperature sensors (Nanmac Co., Washington, DC, USA) are used to monitor the temperature of the four points.The measuring limit of E12-1-C-U temperature sensor is 3300 K, which could satisfy the demands of exhaust plumes temperature test.
A horizontal motor test bench is generally selected for the engine test, as it is easy to design and assures testing security.However, in consideration of the large mass of water, which may result in much more gravitational influence, a vertical test bench is chosen to simulate the rocket launching more truly.Furthermore, a vertical motor test bench makes the measurements more accurate and the test equipment such as high-speed camera and infrared thermograph is more easily located.
Too high an experiment mounting plate tends to render the entire device unsteady.Too low an experiment mounting plate makes it difficult to perform a complete observation.In this study, the height of the mounting plate depends on the length of the core region of the jet: 6.9 1 0.38 (19) where is the length of the core region; is the radius of nozzle exit; represents Mach number and equals 3.5 in this study.
Calculated according to this formula [20], the length of the core region is 1.33 m.The simulation results show that the length of the core region is between 1.25 and 1.4 m.The simulation confirms the effectiveness of Equation (19).Then, considering the length of the motor nozzle, the height of the mounting plate is 1.9 m.
To simulate typical rocket exhaust conditions, the designed Mach number of the free jet in the test is about 3.5 and the outlet temperature is approximately 2000 K. To realize this condition, a solid rocket motor was chosen.The total temperature and total pressure of the combustion chamber are 3000 K and 7 MPa, respectively.The jet flow of this rocket motor is an axisymmetric over-expanded flow, so the rocket motor is designed based on the basic formula of one-dimensional isentropic flow and the combustion characteristic of solid motor grain.The rocket motor is shown in Figure 5.With composite modified double base propellant, the screw-compressed motor grain is available.To ensure a steady thrust, the combustion of internal propellant occurs simultaneously with the external propellant, and both sides of the motor grain are covered by flame-resistant material under constant temperature and pressure.With composite modified double base propellant, the screw-compressed motor grain is available.To ensure a steady thrust, the combustion of internal propellant occurs simultaneously with the external propellant, and both sides of the motor grain are covered by flame-resistant material under constant temperature and pressure.

Free Jet without Water Injection
As shown in Figure 6, two photographs were taken under the combustion steady process conditions.With the over-expanded nature of the jet, five shock cells are clearly seen in the jet core region of free jet flow in.The location of first shock cell is x 1 = 2.1 d e , and the location of the second one is x 2 = 4.9 d e .The temperature in IR image ranges from 0 to 1000 K (Figure 7a).The flow pattern of free jet in high-speed camera and infrared thermograph resemble each other, especially the high-temperature core region, which extends from the exit to the bottom plate.

Free jet without Water Injection
As shown in Figure 6, two photographs were taken under the combustion steady process conditions.With the over-expanded nature of the jet, five shock cells are clearly seen in the jet core region of free jet flow in.The location of first shock cell is 2.1 , and the location of the second one is 4.9 .The temperature in IR image ranges from 0 to 1000 K (Figure 7a).The flow pattern of free jet in high-speed camera and infrared thermograph resemble each other, especially the high-temperature core region, which extends from the exit to the bottom plate.As shown in Figure 6 the high-speed photograph is in accordance with the predicted temperature contour.Moreover, the same number of clear Mach disks is noted.Temperature distribution without water injection is clearly displayed and the temperature of central line is the highest.The temperatures of the testing bench from experiment and that in simulation are consistent with each other in Figure 8.The peak temperature at the center point of the testing bench reaches 2400 K.The monitoring point that is 0.2 m from the center of the testing reaches up to about 1350 K. Obviously, there is extreme shock and ablation on the testing bench owing to the free jet (Figure 7b).As shown in Figure 6 the high-speed photograph is in accordance with the predicted temperature contour.Moreover, the same number of clear Mach disks is noted.Temperature distribution without water injection is clearly displayed and the temperature of central line is the highest.The temperatures of the testing bench from experiment and that in simulation are consistent with each other in Figure 8.The peak temperature at the center point of the testing bench reaches 2400 K.The monitoring point that is 0.2 m from the center of the testing reaches up to about 1350 K. Obviously, there is extreme shock and ablation on the testing bench owing to the free jet (Figure 7b).As shown in Figure 6 the high-speed photograph is in accordance with the predicted temperature contour.Moreover, the same number of clear Mach disks is noted.Temperature distribution without water injection is clearly displayed and the temperature of central line is the highest.The temperatures of the testing bench from experiment and that in simulation are consistent with each other in Figure 8.The peak temperature at the center point of the testing bench reaches 2400 K.The monitoring point that is 0.2 m from the center of the testing reaches up to about 1350 K. Obviously, there is extreme shock and ablation on the testing bench owing to the free jet (Figure 7b).

Jet with Water Injection
The conditions for simulation are the same as those for the experiments.The photos given are digital video images and high-speed photographs.The flat deflector plate was shocked by jet flow without water to cause the plate to be seriously ablated in Figure 9a.As shown in Figure 9b, jets flowing symmetrically from two pipes intersected directly under the engine nozzle, spraying small drops of water all around.A large amount of water vapor began to be produced after injecting water into the jet flow in Figure 9c.There are only two shock cells, which are clearly visible in Figure 10a for the water injection.It is easy to find that the first shock cell is at the same location as the case without water injection 2.1 .However, the interaction between the water spray and the main supersonic combustion gas stream is highly unsteady, so the second one is not fixed as the impact of the water injection, varied between 4.7 and 5.1 .The jet core region is compressed to an isosceles triangle, the peak of which is not stable and the apex angle varied rapidly between 0.5β and 1β.β = 15° is expansion cone half angle.The shape of the jet core region shown in IR (Infrared Radiation) image is similar to that shown in Figure 10b.

Jet with Water Injection
The conditions for simulation are the same as those for the experiments.The photos given are digital video images and high-speed photographs.The flat deflector plate was shocked by jet flow without water to cause the plate to be seriously ablated in Figure 9a.As shown in Figure 9b, jets flowing symmetrically from two pipes intersected directly under the engine nozzle, spraying small drops of water all around.A large amount of water vapor began to be produced after injecting water into the jet flow in Figure 9c.

Jet with Water Injection
The conditions for simulation are the same as those for the experiments.The photos given are digital video images and high-speed photographs.The flat deflector plate was shocked by jet flow without water to cause the plate to be seriously ablated in Figure 9a.As shown in Figure 9b, jets flowing symmetrically from two pipes intersected directly under the engine nozzle, spraying small drops of water all around.A large amount of water vapor began to be produced after injecting water into the jet flow in Figure 9c.There are only two shock cells, which are clearly visible in Figure 10a for the water injection.It is easy to find that the first shock cell is at the same location as the case without water injection 2.1 .However, the interaction between the water spray and the main supersonic combustion gas stream is highly unsteady, so the second one is not fixed as the impact of the water injection, varied between 4.7 and 5.1 .The jet core region is compressed to an isosceles triangle, the peak of which is not stable and the apex angle varied rapidly between 0.5β and 1β.β = 15° is expansion cone half angle.The shape of the jet core region shown in IR (Infrared Radiation) image is similar to that shown in Figure 10b.There are only two shock cells, which are clearly visible in Figure 10a for the water injection.It is easy to find that the first shock cell is at the same location as the case without water injection .
x 1 " 2.1d e .However, the interaction between the water spray and the main supersonic combustion gas stream is highly unsteady, so the second one is not fixed as the impact of the water injection, .
x 2 varied between 4.7d e and 5.1d e .The jet core region is compressed to an isosceles triangle, the peak of which is not stable and the apex angle varied rapidly between 0.5β and 1β.β = 15 ˝is expansion cone half angle.The shape of the jet core region shown in IR (Infrared Radiation) image is similar to that shown in Figure 10b.Moreover, the numerical results agree well with the experimental results at monitoring points (Table 1).The maximum error between the two results is within 10%.Through comparison of data and temperature-distribution pictures between the experiments and simulations, the simulation model is suitable for verifying the experiment.

Temperature Comparison
As shown in Figure 12a, the total temperature of free jet is compared with that of jet with water injection in the plane-XOY.The temperature comparison between free jet and jet with water injection In experiment and simulation, mass flow rate of each water pipe is Q w = 1.8 kg/s and mass flow rate of jet is Q g = 1.5 kg/s.Outlet water velocity is 16 m/s and there are two pipes.Because of the water injection, the length of the jet core region in the state of water injection is shorter than that of free jet (Figure 11).The jet flow shape and temperature distribution is consistent in the two pictures.Moreover, the numerical results agree well with the experimental results at monitoring points (Table 1).The maximum error between the two results is within 10%.Through comparison of data and temperature-distribution pictures between the experiments and simulations, the simulation model is suitable for verifying the experiment.

Temperature Comparison
As shown in Figure 12a, the total temperature of free jet is compared with that of jet with water injection in the plane-XOY.The temperature comparison between free jet and jet with water injection Moreover, the numerical results agree well with the experimental results at monitoring points (Table 1).The maximum error between the two results is within 10%.Through comparison of data and temperature-distribution pictures between the experiments and simulations, the simulation model is suitable for verifying the experiment.

Temperature Comparison
As shown in Figure 12a, the total temperature of free jet is compared with that of jet with water injection in the plane-XOY.The temperature comparison between free jet and jet with water injection in the plane-XOZ is shown in Figure 12b.The X-axis is vertical axis of the nozzle.The plane of two water pipes is named plane-XOY and its perpendicular plane is named plane-XOZ.
Energies 2015, 8, page-page 12 in the plane-XOZ is shown in Figure 12b.The X-axis is vertical axis of the nozzle.The plane of two water pipes is named plane-XOY and its perpendicular plane is named plane-XOZ.The external zone of the jet with injected water is a subsonic one, but the internal zone is a supersonic one, so the influence of water on external zone is greater than that on internal zone.The length of flow core region is shortened by injecting water.The energy of the core region hardly transmits outward through heat transfer means that the temperature distribution of the free jet is similar to that of the jet with water injection near the X-axis.The diameter of the pipe is constant and two pipes are only installed in the XOY plane.The location of pipes causes jet asymmetry, so the interphase interaction between liquid and gas in the XOZ plane decreases and the temperature distribution of the flow is wider when compared with that of the XOY plane.
The axial temperature of jet flow shows a serious decrease in the second half of the flow in Figure 13a.The point A is the water injection location on the X-axis.This point and the nozzle exit are 0.592 m apart.The maximum temperature of the axial temperature exhibits two peaks.The first one is located at the same place as in the free jet case; the second one lies on the place of 0.75 m of distance from the exit plate, which is the apex of an isosceles triangle jet core.The temperature decreased rapidly after the second peak value due to the endothermic water vaporization process.The maximum reduction of temperature is nearly 70% of the free jet case and the mean reduction of temperature is about 50% in the center line.The predicted average temperature at the central line of the jet core region illustrates that the cooling effect is obvious.This is clearly reported by infrared thermography and the temperature contour.The temperature of the bench without water is shown as a line-free in Figure 13b.In the case with water injection, line Y is parallel to the XOY plane and line Z is parallel to XOZ plane.There is The external zone of the jet with injected water is a subsonic one, but the internal zone is a supersonic one, so the influence of water on external zone is greater than that on internal zone.The length of flow core region is shortened by injecting water.The energy of the core region hardly transmits outward through heat transfer means that the temperature distribution of the free jet is similar to that of the jet with water injection near the X-axis.The diameter of the pipe is constant and two pipes are only installed in the XOY plane.The location of pipes causes jet asymmetry, so the interphase interaction between liquid and gas in the XOZ plane decreases and the temperature distribution of the flow is wider when compared with that of the XOY plane.
The axial temperature of jet flow shows a serious decrease in the second half of the flow in Figure 13a.The point A is the water injection location on the X-axis.This point and the nozzle exit are 0.592 m apart.The maximum temperature of the axial temperature exhibits two peaks.The first one is located at the same place as in the free jet case; the second one lies on the place of 0.75 m of distance from the exit plate, which is the apex of an isosceles triangle jet core.The temperature decreased rapidly after the second peak value due to the endothermic water vaporization process.The maximum reduction of temperature is nearly 70% of the free jet case and the mean reduction of temperature is about 50% in the center line.The predicted average temperature at the central line of the jet core region illustrates that the cooling effect is obvious.This is clearly reported by infrared thermography and the temperature contour.The external zone of the jet with injected water is a subsonic one, but the internal zone is a supersonic one, so the influence of water on external zone is greater than that on internal zone.The length of flow core region is shortened by injecting water.The energy of the core region hardly transmits outward through heat transfer means that the temperature distribution of the free jet is similar to that of the jet with water injection near the X-axis.The diameter of the pipe is constant and two pipes are only installed in the XOY plane.The location of pipes causes jet asymmetry, so the interphase interaction between liquid and gas in the XOZ plane decreases and the temperature distribution of the flow is wider when compared with that of the XOY plane.
The axial temperature of jet flow shows a serious decrease in the second half of the flow in Figure 13a.The point A is the water injection location on the X-axis.This point and the nozzle exit are 0.592 m apart.The maximum temperature of the axial temperature exhibits two peaks.The first one is located at the same place as in the free jet case; the second one lies on the place of 0.75 m of distance from the exit plate, which is the apex of an isosceles triangle jet core.The temperature decreased rapidly after the second peak value due to the endothermic water vaporization process.The maximum reduction of temperature is nearly 70% of the free jet case and the mean reduction of temperature is about 50% in the center line.The predicted average temperature at the central line of the jet core region illustrates that the cooling effect is obvious.This is clearly reported by infrared thermography and the temperature contour.The temperature of the bench without water is shown as a line-free in Figure 13b.In the case with water injection, line Y is parallel to the XOY plane and line Z is parallel to XOZ plane.There is The temperature of the bench without water is shown as a line-free in Figure 13b.In the case with water injection, line Y is parallel to the XOY plane and line Z is parallel to XOZ plane.There is an interesting phenomenon that the monitoring line temperature of line Z is slightly below that of line Y (Figure 13b).The reason is that the contact area between water and the expanded jet flow without water shock is great.The energy loss increases with more vaporization than that of the XOZ plane.The monitoring point temperature of the flow center is about 480 K and is down to 20%.The monitoring point temperature at 0.2 m drops to 34.6%.The average temperature of the bench is 100-200 K.It is concluded that the effect of cooling in the surrounding places as well as that in the center and the temperature distribution of the bottom plate does produce a drastic change like the free jet case.

Results and Analysis
Table 2 shows different cases of water injection considered.The cooling effects caused in the various conditions are shown below.

Effect of Water Injection Velocity
Cases A, B, C and D have the same mass flow rate and pipe number but different velocities.Figure 14a shows the vaporization rate in the four cases.The vaporization rate gradually rises and becomes stable, and in the stable stage, there is pulsating change periodically.
Energies 2015, 8, page-page 13 an interesting phenomenon that the monitoring line temperature of line Z is slightly below that of line Y (Figure 13b).The reason is that the contact area between water and the expanded jet flow without water shock is great.The energy loss increases with more vaporization than that of the XOZ plane.The monitoring point temperature of the flow center is about 480 K and is down to 20%.The monitoring point temperature at 0.2 m drops to 34.6%.The average temperature of the bench is 100-200 K.It is concluded that the effect of cooling in the surrounding places as well as that in the center and the temperature distribution of the bottom plate does produce a drastic change like the free jet case.

Results and Analysis
Table 2 shows different cases of water injection considered.The cooling effects caused in the various conditions are shown below.

Effect of Water Injection Velocity
Cases A, B, C and D have the same mass flow rate and pipe number but different velocities.Figure 14a shows the vaporization rate in the four cases.The vaporization rate gradually rises and becomes stable, and in the stable stage, there is pulsating change periodically.13206 The vast amount of water vaporizes instantly impacting the main stream.The ratio of transient impulse changes and the contact area also changes between the gas phase and water phase.Then, the momentum ratio of the two phases is back to the initial one because of constant water resupply, so the vaporization ratio changes periodically.In cases A, B and C, the lower the water velocity is, the lower the vaporization rate is.However, the vaporization rate of case D is close to that of case B. The phenomenon is caused by the higher water velocity in case D. There are two branching flows in Figure 15, so the contact area decreases and vaporization rate is lower than that of case C. As shown in Figure 14b, the mass of remaining water does not change when the vaporization rate tends to be stable.The temperature along the jet axis decreases as the temperature decreases.It indicates that it is available to improve the efficiency of the vaporization rate by appropriately choosing the water velocity.
Energies 2015, 8, page-page 14 The vast amount of water vaporizes instantly impacting the main stream.The ratio of transient impulse changes and the contact area also changes between the gas phase and water phase.Then, the momentum ratio of the two phases is back to the initial one because of constant water resupply, so the vaporization ratio changes periodically.In cases A, B and C, the lower the water velocity is, the lower the vaporization rate is.However, the vaporization rate of case D is close to that of case B. The phenomenon is caused by the higher water velocity in case D. There are two branching flows in Figure 15, so the contact area decreases and vaporization rate is lower than that of case C. As shown in Figure 14b, the mass of remaining water does not change when the vaporization rate tends to be stable.The temperature along the jet axis decreases as the temperature decreases.It indicates that it is available to improve the efficiency of the vaporization rate by appropriately choosing the water velocity.

Effect of Water Mass Flow Rate
In Table 2, the water mass flow rate in case D is about 2.5 times more than that in case C, and mass flow rate of water in case E is lower than that in case F. It is obvious that cooling effect relates to this condition in Figures 15 and 16

Effect of Exhaust Jet Momentum
Table 3 shows that the vaporization rate in case H is the highest when compared with those in other cases.Before the formation of the two branching flows, the vaporization rate rises with increasing outlet water velocity, but after passing the tipping point, the vaporization rate gradually drops with increasing outlet water velocity.

Effect of Water Mass Flow Rate
In Table 2, the water mass flow rate in case D is about 2.5 times more than that in case C, and mass flow rate of water in case E is lower than that in case F. It is obvious that cooling effect relates to this condition in Figures 15 and 16.The temperature decreases with increasing mass flow rate.
Energies 2015, 8, page-page 14 The vast amount of water vaporizes instantly impacting the main stream.The ratio of transient impulse changes and the contact area also changes between the gas phase and water phase.Then, the momentum ratio of the two phases is back to the initial one because of constant water resupply, so the vaporization ratio changes periodically.In cases A, B and C, the lower the water velocity is, the lower the vaporization rate is.However, the vaporization rate of case D is close to that of case B. The phenomenon is caused by the higher water velocity in case D. There are two branching flows in Figure 15, so the contact area decreases and vaporization rate is lower than that of case C. As shown in Figure 14b, the mass of remaining water does not change when the vaporization rate tends to be stable.The temperature along the jet axis decreases as the temperature decreases.It indicates that it is available to improve the efficiency of the vaporization rate by appropriately choosing the water velocity.

Effect of Water Mass Flow Rate
In Table 2, the water mass flow rate in case D is about 2.5 times more than that in case C, and mass flow rate of water in case E is lower than that in case F. It is obvious that cooling effect relates to this condition in Figures 15 and 16

Effect of Exhaust Jet Momentum
Table 3 shows that the vaporization rate in case H is the highest when compared with those in other cases.Before the formation of the two branching flows, the vaporization rate rises with increasing outlet water velocity, but after passing the tipping point, the vaporization rate gradually drops with increasing outlet water velocity.

Effect of Exhaust Jet Momentum
Table 3 shows that the vaporization rate in case H is the highest when compared with those in other cases.Before the formation of the two branching flows, the vaporization rate rises with increasing outlet water velocity, but after passing the tipping point, the vaporization rate gradually drops with increasing outlet water velocity.
The dimensional polylines between the momentum ratio and vaporization rate are based on Table 3.In Figure 17, the jet flow showa compressive deformation and two branching flows were formed in all simulations, so the polylines are in accordance with the relevant physical laws.With the same conditions except outlet water velocity, the best cooling effect appears on around momentum ratio 50.
The dimensional polylines between the momentum ratio and vaporization rate are based on Table 3.In Figure 17, the jet flow showa compressive deformation and two branching flows were formed in all simulations, so the polylines are in accordance with the relevant physical laws.With the same conditions except outlet water velocity, the best cooling effect appears on around momentum ratio 50.

Conclusions
This investigation addresses the concept of injecting water into the rocket exhaust plume to reduce the jet temperature.A series of experiments of injecting water into the exhaust plume, which were performed in a controlled laboratory environment, were run.A mixture multiphase model, which introduces the energy source items caused by the vaporization of liquid water into the energy equation, was developed to simulate the flow field of the complex gas-liquid multiphase flow.This multiphase model was validated by the experimental results.The flow field, distribution of temperature in both simulation and experiment, and cooling effect were discussed.
It is seen that the high temperature core region of the plume is compressed to an isosceles triangle by water injection in both the numerical simulations and experiments.The length and the area of the core are smaller than those of the free jet, respectively, especially the region whose temperature is higher than 1000 K.The main energy of the plume jet is concentrated up the high temperature core region, where it is feasible to cool the rocket exhaust plume by evaporation of water injection.
The stagnation temperature of the plume sometimes leads to damage of launcher structures due to high temperature erosion.As a result, the cooling effect on the bottom plate warrants further study.The stagnation temperature disappears in the center of the bottom plate for the water injection, and Figure 14 shows that the effect of cooling in the center is better, and the temperature distribution of the bottom plate departs considerably from the free jet case.
In the different cases considered, the varying parameters like mass flow rate, pipe number, and outlet water velocity can produce temperature distribution variations on the bottom plate.The

Conclusions
This investigation addresses the concept of injecting water into the rocket exhaust plume to reduce the jet temperature.A series of experiments of injecting water into the exhaust plume, which were performed in a controlled laboratory environment, were run.A mixture multiphase model, which introduces the energy source items caused by the vaporization of liquid water into the energy equation, was developed to simulate the flow field of the complex gas-liquid multiphase flow.This multiphase model was validated by the experimental results.The flow field, distribution of temperature in both simulation and experiment, and cooling effect were discussed.
It is seen that the high temperature core region of the plume is compressed to an isosceles triangle by water injection in both the numerical simulations and experiments.The length and the area of the core are smaller than those of the free jet, respectively, especially the region whose temperature is higher than 1000 K.The main energy of the plume jet is concentrated up the high temperature core region, where it is feasible to cool the rocket exhaust plume by evaporation of water injection.
The stagnation temperature of the plume sometimes leads to damage of launcher structures due to high temperature erosion.As a result, the cooling effect on the bottom plate warrants further study.The stagnation temperature disappears in the center of the bottom plate for the water injection, and Figure 14 shows that the effect of cooling in the center is better, and the temperature distribution of the bottom plate departs considerably from the free jet case.
In the different cases considered, the varying parameters like mass flow rate, pipe number, and outlet water velocity can produce temperature distribution variations on the bottom plate.The momentum of water was chosen as an important parameter, which influences the cooling effect by comparative analysis.A scaling equation of a dimensionless relationship between reduction in mean temperature and momentum ratio is calculated to predict the cooling effect of subscale or full-scale launching test.
It should be noted that this study has not examined the cooling effect of the injection angle and the water injection orifice diameter.Further study is needed to investigate the effect of the variable multiphase flow field.Notwithstanding its limitations, this study does confirm that water injection is an effective coolant for rocket motor exhaust.At the same time, a preliminary study of the mechanism of the interaction of water and combustion gas and the relationship between cooling effect and the momentum ratio were developed, which could be used to predict the results in a full-scale launch test.

Figure 1 .
Figure 1.Schematic of the boundary condition.

Figure 1 .
Figure 1.Schematic of the boundary condition.

Figure 2 .Figure 3 .
Figure 2. Mesh structure at the symmetry plane of the missile model.(a) Model A; (b) Model B; (c) Model C.

Figure 2 .
Figure 2. Mesh structure at the symmetry plane of the missile model.(a) Model A; (b) Model B; (c) Model C.

Figure 2 .Figure 3 .
Figure 2. Mesh structure at the symmetry plane of the missile model.(a) Model A; (b) Model B; (c) Model C.

Figure 8 .
Figure 8.Comparison chart of monitor values and computational values.

Figure 9 .
Figure 9. (a) Jet flow model without water; (b) Water injection model without jet flow; (c) Water injection model with jet flow.

Figure 8 .
Figure 8.Comparison chart of monitor values and computational values.

Figure 9 .
Figure 9. (a) Jet flow model without water; (b) Water injection model without jet flow; (c) Water injection model with jet flow.

Figure 9 .
Figure 9. (a) Jet flow model without water; (b) Water injection model without jet flow; (c) Water injection model with jet flow.

Figure 10 .Figure 11 .
Figure 10.(a) High-speed photograph; (b) IR image.In experiment and simulation, mass flow rate of each water pipe is 1.8 kg/s and mass flow rate of jet is 1.5 kg/s.Outlet water velocity is 16 m/s and there are two pipes.Because of the water injection, the length of the jet core region in the state of water injection is shorter than that of free jet (Figure11).The jet flow shape and temperature distribution is consistent in the two pictures.

Figure 10 .Figure 11 .
Figure 10.(a) High-speed photograph; (b) IR image.In experiment and simulation, mass flow rate of each water pipe is 1.8 kg/s and mass flow rate of jet is 1.5 kg/s.Outlet water velocity is 16 m/s and there are two pipes.Because of the water injection, the length of the jet core region in the state of water injection is shorter than that of free jet (Figure11).The jet flow shape and temperature distribution is consistent in the two pictures.

Figure 11 .
Figure 11.(a) Experimental jet with water injection; (b) Temperature contour with water injection.

Figure 12 .
Figure 12.(a) Temperature comparison between free jet and jet with water injection in the XOY plane; (b) Temperature comparison between free jet and jet with water injection in the XOZ plane.

Figure 12 .
Figure 12.(a) Temperature comparison between free jet and jet with water injection in the XOY plane; (b) Temperature comparison between free jet and jet with water injection in the XOZ plane.

Energies 2015, 8 ,Figure 12 .
Figure 12.(a) Temperature comparison between free jet and jet with water injection in the XOY plane; (b) Temperature comparison between free jet and jet with water injection in the XOZ plane.

Figure 14 .
Figure 14.(a) Vaporization rate; (b) Mass of the remaining water; (c) Temperature of axial gas stream.

Figure 14 .
Figure 14.(a) Vaporization rate; (b) Mass of the remaining water; (c) Temperature of axial gas stream.

Figure 15 .
Figure 15.Temperature contour of the XOZ plane in case D.

Figure 16 .
Figure 16.(a) Temperature of axial gas stream; (b) Temperature of monitoring line.

Figure 15 .
Figure 15.Temperature contour of the XOZ plane in case D.

Figure 15 .
Figure 15.Temperature contour of the XOZ plane in case D.

Figure 16 .
Figure 16.(a) Temperature of axial gas stream; (b) Temperature of monitoring line.

Figure 16 .
Figure 16.(a) Temperature of gas stream; (b) Temperature of monitoring line.

Figure 17 .
Figure 17.Dimensional polylines between momentum ratio and vaporization rate.

Table 1 .
Comparison between experimental data and numerical data.

Table 1 .
Comparison between experimental data and numerical data.

Table 1 .
Comparison between experimental data and numerical data.
Figure 17.Dimensional polylines between momentum ratio and vaporization rate.