Modeling study on the flow, heat transfer and energy conversion characteristics of low-power arc-heated hydrogen/nitrogen thrusters

A modeling study is conducted to investigate the effect of hydrogen content in propellants on the plasma flow, heat transfer and energy conversion characteristics of low-power (kW class) arc-heated hydrogen/nitrogen thrusters (arcjets). 1:0 (pure hydrogen), 3:1 (to simulate decomposed ammonia), 2:1 (to simulate decomposed hydrazine) and 0:1 (pure nitrogen) hydrogen/nitrogen mixtures are chosen as the propellants. Both the gas flow region inside the thruster nozzle and the anode-nozzle wall are included in the computational domain in order to better treat the conjugate heat transfer between the gas flow region and the solid wall region. The axial variations of the enthalpy flux, kinetic energy flux, directed kinetic-energy flux, and momentum flux, all normalized to the mass flow rate of the propellant, are used to investigate the energy conversion process inside the thruster nozzle. The modeling results show that the values of the arc voltage, the gas axial-velocity at the thruster exit, and the specific impulse of the arcjet thruster all increase with increasing hydrogen content in the propellant, but the gas temperature at the nitrogen thruster exit is significantly higher than that for other three propellants. The flow, heat transfer, and energy conversion processes taking place in the thruster nozzle have some common features for all the four propellants. The propellant is heated mainly in the near-cathode and constrictor region, accompanied with a rapid increase of the enthalpy flux, and after achieving its maximum value, the enthalpy flux decreases appreciably due to the conversion of gas internal energy into its kinetic energy in the divergent segment of the thruster nozzle. The kinetic energy flux, directed kinetic energy flux and momentum flux also increase at first due to the arc heating and the thermodynamic expansion, assume their maximum inside the nozzle and then decrease gradually as the propellant flows toward the thruster exit. It is found that a large energy loss (31-52%) occurs in the thruster nozzle due to the heat transfer to the nozzle wall and too long nozzle is not necessary. Modeling results for the NASA 1-kW class arcjet thruster with hydrogen or decomposed hydrazine as the propellant are found to compare favorably with available experimental data.

Modeling study to compare the flow and heat transfer characteristics of low-power hydrogen, nitrogen and argon arc-heated thrusters

A modelling study is performed to compare the plasma °ow and heat transfer char- acteristics of low-power arc-heated thrusters (arcjets) for three di®erent propellants: hydrogen, nitrogen and argon. The all-speed SIMPLE algorithm is employed to solve the governing equa- tions, which take into account the e®ects of compressibility, Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. The temperature, veloc- ity and Mach number distributions calculated within the thruster nozzle obtained with di®erent propellant gases are compared for the same thruster structure, dimensions, inlet-gas stagnant pressure and arc currents. The temperature distributions in the solid region of the anode-nozzle wall are also given. It is found that the °ow and energy conversion processes in the thruster nozzle show many similar features for all three propellants. For example, the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appear- ing near the cathode tip; the °ow transition from the subsonic to supersonic regime occurs within the constrictor region; the highest axial velocity appears inside the nozzle; and most of the input propellant °ows towards the thruster exit through the cooler gas region near the anode-nozzle wall. However, since the properties of hydrogen, nitrogen and argon, especially their molecular weights, speci¯c enthalpies and thermal conductivities, are di®erent, there are appreciable di®er- ences in arcjet performance. For example, compared to the other two propellants, the hydrogen arcjet thruster shows a higher plasma temperature in the arc region, and higher axial velocity but lower temperature at the thruster exit. Correspondingly, the hydrogen arcjet thruster has the highest speci¯c impulse and arc voltage for the same inlet stagnant pressure and arc current. The predictions of the modelling are compared favourably with available experimental results.

Numerical modeling of heat transfer and plasma flow in a low power arcjet thruster

Modeling studies are preformed to investigate the plasma and heat transfer characteristics of a low power argon arcjet thruster. Computed temperature, velocity, static pressure, and Mach number distribution in arcjet thruster under typical operating condition are presented in this paper. It shows that the performance data from numerical modeling results are basically consistent with the experimental measured values.

Fully developed laminar flow and heat transfer in smooth trapezoidal microchannel

Numerical analysis of fully developed laminar slip flow and heat transfer in trapezoidal micro-channels has been studied with uniform wall heat flux boundary conditions. Through coordinate transformation, the governing equations are transformed from physical plane to computational domain, and the resulting equations are solved by a finite-difference scheme. The influences of velocity slip and temperature jump on friction coefficient and Nusselt number are investigated in detail. The calculation also shows that the aspect ratio and base angle have significant effect on flow and heat transfer in trapezoidal micro-channel. (c) 2005 Elsevier Ltd. All rights reserved.

Conjugate Model for Heat and Mass Transfer of Porous Wall in the High Temperature Gas Flow

Heat and mass transfer of a porous permeable wall in a high temperature gas dynamical flow is considered. Numerical simulation is conducted on the ground of the conjugate mathematical model which includes filtration and heat transfer equations in a porous body and boundary layer equations on its surface. Such an approach enables one to take into account complex interaction between heat and mass transfer in the gasdynamical flow and in the structure subjected to this flow. The main attention is given to the impact of the intraporous heat transfer intensity on the transpiration cooling efficiency.

Visualization of dust explosion under microgravity conditions

The objective of this work was to apply visualization methods to the experimental study of cornstarch dust-air mixture combustion in a closed vessel volume under microgravity conditions. A dispersion system with a small scale of turbulence was used in the experiments. A gas igniter initiated combustion of the dust-air mixture in the central or top part of the vessel. Flame propagation through the quiescent mixture was recorded by a high-speed video camera. Experiments showed a very irregular flame front and irregular distribution of the regions with local reactions of re-burning behind the flame front. at a later stage of combustion. Heat transfer from the hot combustion products to the walls is shown to have an important role in the combustion development. The maximum pressure and maximum rate of pressure rise were higher for flame propagation from the vessel center than for flame developed from the top pan of the vessel. The reason for smaller increase of the rate of pressure rise, for the flame developed from the top of the vessel. in comparison with that developed from the vessel center, was much faster increase of the contact surface of the combustion gases with the vessel walls. It was found that in dust flames only small part of hear was released at the flame front, the remaining part being released far behind it.

Coupling Of Evaporation And Thermocapillary Convection In A Liquid Layer With Mass And Heat Exchanging Interface

We propose and analyse a new model of thermocapillary convection with evaporation in a cavity subjected to horizontal temperature gradient, rather than the previously studied model without evaporation. The pure liquid layer with a top free surface in contact with its own vapour is considered in microgravity condition. The computing programme developed for simulating this model integrates the two-dimensional, time-dependent Navier-Stokes equations and energy equation by a second-order accurate projection method. We focus on the coupling of evaporation and thermocapillary convection by investigating the influence of evaporation Biot number and Marangoni number on the interfacial mass and heat transfer. Three different regimes of the coupling mechanisms are found and explained from our numerical results.

Basic equations for suspension flows with interphase mass-transfer

In the case of suspension flows, the rate of interphase momentum transfer M(k) and that of interphase energy transfer E(k), which were expressed as a sum of infinite discontinuities by Ishii, have been reduced to the sum of several terms which have concise physical significance. M(k) is composed of the following terms: (i) the momentum carried by the interphase mass transfer; (ii) the interphase drag force due to the relative motion between phases; (iii) the interphase force produced by the concentration gradient of the dispersed phase in a pressure field. And E(k) is composed of the following four terms, that is, the energy carried by the interphase mass transfer, the work produced by the interphase forces of the second and third parts above, and the heat transfer between phases. It is concluded from the results that (i) the term, (-alpha-k-nabla-p), which is related to the pressure gradient in the momentum equation, can be derived from the basic conservation laws without introducing the "shared-pressure presumption"; (ii) the mean velocity of the action point of the interphase drag is the mean velocity of the interface displacement, upsilonBAR-i. It is approximately equal to the mean velocity of the dispersed phase, upsilonBAR-d. Hence the work terms produced by the drag forces are f(dc) . upsilonBAR-d, and f(cd) . upsilonBAR-d, respectively, with upsilonBAR-i not being replaced by the mean velocity of the continuous phase, upsilonBAR-c; (iii) by analogy, the terms of the momentum transfer due to phase change are upsilonBAR-d-GAMMA-c, and upsilonBAR-d-GAMMA-d, respectively; (iv) since the transformation between explicit heat and latent heat occurs in the process of phase change, the algebraic sum of the heat transfer between phases is not equal to zero. Q(ic) and Q(id) are composed of the explicit heat and latent heat, so that the sum Q(ic) + Q(id)) is equal to zero.

Space experimental studies of microgravity fluid science in China

Microgravity fluid physics is an important part of microgravity sciences, which consists of simple fluids of many new systems, gas-liquid two-phase flow and heat transfer, and complex fluid mechanics. In addition to the importance of itself in sciences and applications, microgravity fluid physics closely relates to microgravity combustion, space biotechnology and space materials science, and promotes the developments of interdisciplinary fields. Many space microgravity experiments have been per- formed on board the recoverable satellites and space ships of China and pushed the rapid development of microgravity sciences in China. In the present paper, space experimental studies and the main re- sults of the microgravity fluid science in China in the last 10 years or so are introduced briefly.

Numerical simulations of interrupted and conventional microchannel heat sinks

We provide three-dimensional numerical simulations of conjugate heat transfer in conventional and the newly proposed interrupted microchannel heat sinks. The new microchannel heat sink consists of a set of separated zones adjoining shortened parallel microchannels and transverse microchambers. Multi-channel effect, physical property variations, and axial thermal conduction are considered. It is found that flow rate variations in different channels can be neglected, while heat received by different channels accounts for 2% deviations from the averaged value when the heat flux at the back surface of the silicon chip reaches 100 W/cm(2). The computed hydraulic and thermal boundary layers are redeveloping in each separated zone due to shortened flow length for the interrupted microchannel heat sink. The periodic thermal developing flow is responsible for the significant heat transfer enhancement. Two effects influence pressure drops across the newly proposed microchannel heat sink. The first one is the pressure recovery effect in the microchamber, while the second one is the head loss when liquid leaves the microchamber and enters the next zone. The first effect compensates or suppresses the second one, leading to similar or decreased pressure drop than that for the conventional microchannel heat sink, with the fluid Prandtl number larger than unity.

Combined Effects of Hot Curing Conditions and Reaction Heat on Rubber Vulcanization Efficiency and Vulcanizate Uniformity

A mathematical model of the chemical kinetics of silicone rubber Vulcanization is developed, with the thermal effects being computed using the increment method, and the hot Vulcanization process estimated with the finite element method. The results show that the reaction heat of rubber vulcanization is important for energy saving, and that a proper curing medium temperature is important when considering both vulcanization efficiency and vulcanizate uniformity. The results also indicate that increases in the forced convective heat transfer coefficient have no significant effect above a certain level. The validity of the numerical model is indirectly proven by comparison with existing data.

One century record of contamination by polycyclic aromatic hydrocarbons and polychlorinated biphenyls in core sediments from the southern Yellow Sea

Sixteen polycyclic aromatic hydrocarbons (PAHs) and 28 polychlorinated biphenyls (PCBs) were measured at a 2-cm interval in a core sample from the middle of the southern Yellow Sea for elucidating their historical variations in inflow and sources. The chronology was obtained using the Pb-210 method. PAHs concentrations decreased generally with depth and two climax values occurred in 14-16 cm and 20-22 cm layers, demonstrating that the production and usage of PAHs might reach peaks in the periods of 1956-1962 and 1938-1944. The booming economy and the navy battles of the Second World War might explain why the higher levels were detected in the two layers. The result of principal component analysis (PCA) revealed that PAHs were primarily owing to the combustion product. Down-cored variation of PCB concentrations was complex. Higher concentrations besides the two peaks being the same as PAHs were detected from 4 to 8 cm, depositing from 1980 to 1992, which probably resulted from the disposal of the out-dated PCB-containing equipment. The average Cl percentage of PCBs detected was similar to that of the mixture of Aroclor 1254 and 1242, suggesting they might origin from the dielectrical and heat-transfer fluid. The total organic carbon (TOC) content played a prevalent role in the adsorption of high molecular weight PAHs (>= 4-ring), while no obvious relationship among total PCBs, the concentration of congeners, and TOC was found.

Methanol steam reforming in a compact plate-fin reformer for fuel-cell systems

A compact plate-fin reformer (PFR) consisting of closely spaced plate-fins, in which endothermic and exothermic reactions take place in alternate chambers, has been studied. In the PFR, which was based on a plate-fin heat exchanger, catalytic combustion of the reforming gas, as a simulation of the fuel cell anode off gas (AOG), supplied the necessary heat for the reforming reaction. One reforming chamber, which was for hydrogen production, was integrated with two vaporization chambers and two combustion chambers to constitute a single unit of PFR. The PFR is very compact, easy to be placed and scaled up. The effect of the ratio of H2O/CH3OH on the performance of the PFR has been investigated, and temperature distributions in different chambers were studied. Besides, the stationary behavior of the PFR was also investigated. Heat transfer of the reformer was enhanced by internal plate-fins as well as by external catalytic combustion, which offer both high methanol conversion ratio and low CO concentration. In addition, the fully integrated reformer exhibited good test stability. Based on the PFR, a scale-up reformer was designed and operated continuously for 1000 h, with high methanol conversion ratio and low CO concentration. (c) 2004 International Association for Hydrogen Energy. Published by Elsevier Ltd. All rights reserved.

Effect of ambient pressure and radiation reabsorption of atmosphere on the flame spreading over thermally thin combustibles in microgravity

For the flame spread over thermally thin combustibles in an atmosphere, if the atmosphere cannot emit and absorb the thermal radiation (e.g. for atmosphere Of O-2-N-2), the conductive heat transfer from the flame to the fuel surface dominates the flame spread at lower ambient atmosphere. As the ambient pressure increases, the flame spread rate increases, and the radiant heat transfer from the flame to the fuel surface gradually becomes the dominant driving force for the flame spread. In contrast, if the atmosphere is able to emit and absorb the thermal radiation (e.g. for atmosphere Of O-2-CO2), at lower pressure, the heat transfer from flame to the fuel surface is enhanced by the radiation reabsorption of the atmosphere at the leading edge of the flame, and both conduction and thermal radiation play important roles in the mechanism of flame spread. With the increase in ambient pressure, the oxygen diffuses more quickly from ambient atmosphere into the flame, the chemical reaction in the flame is enhanced, and the flame spread rate increases. When the ambient pressure is greater than a critical value, the thermal radiation from the flame to the solid surface is hampered by the radiation reabsorption of ambient atmosphere with the further increase in ambient pressure. As a result, with the increase in ambient pressure, the flame spread rate decreases and the heat conduction gradually dominates the flame spread over the fuel surface.