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.

Numerical Simulation of a Low-power Hydrazine Arcjet Thruster

A modeling study is conducted to investigate the plasma flow and heat transfer characteristics of low-power (kW class) arc-heated thrusters (arcjets) with 2:1 hydrogen/nitrogen to simulate decomposed hydrazine as the propellant. The all-speed SIMPLE algorithm is employed to solve the governing equations, which take into account the effects of compressibility, the Lorentz force and Joule heating, as well as the temperature- and pressure-dependence of the gas properties. Typical computed results about the temperature, velocity and Mach number distributions within arcjet thruster are presented for the case with arc current of 9 A and inlet stagnant pressure of 3.3×105 Pa to show the flow and heat transfer characteristics. It is found that the propellant is heated mainly in the near-cathode and constrictor region, with the highest plasma temperature appearing near the cathode tip, and the flow transition from the subsonic to supersonic regime occurs within the constrictor region. The effect of gas viscosity on the plasma flow within arcjet thruster is examined by an additional numerical test using artificially reduced values of gas viscosity. The test results show that the gas viscosity appreciably affects the plasma flow and the performance of the arcjet thruster for the cases with the hydrazine or hydrogen as the propellant. The integrated axial Lorentz force in the thruster nozzle is also calculated and compared with the thrust force of the arcjet thruster. It is found that the integrated axial Lorentz force is much smaller than the thrust force for the low-power arcjet thruster. Modeling results for the NASA 1-kW class arcjet thruster with simulated hydrazine as the propellant are found to be reasonably consistent with available experimental data.

A 2-kW COIL with a square pipe-array JSOG and nitrogen buffer gas

A 2-kW-class chemical oxygen-iodine laser (COIL) using nitrogen buffer gas has been developed and tested since industrial applications of COIL devices will require the use of nitrogen as the buffer gas. The laser, with a gain length of 11.7 cm, is energized by a square pipe-array jet-type singlet oxygen generator (SPJSOG) and employs a nozzle bank with a designed Mach number of 2.5. The SPJSOG has advantages over the traditional plate-type JSOG in that it has less requirements on basic hydrogen peroxide (BHP) pump, and more important, it has much better operational stability. The SPJSOG without a cold trap and a gas-liquid separator could provide reliable operations for a total gas flow rate up to 450 mmol/s and with a low liquid driving pressure of around 0.7 atm or even lower. The nozzle bank was specially designed for a COIL using nitrogen as the buffer gas. The cavity was designed for a Mach number of 2.5, in order to provide a gas speed and static temperature in the cavity similar to that for a traditional COIL with helium buffer gas and a Mach 2 nozzle. An output power of 2.6 kW was obtained for a chlorine flow rate of 140 mmol/s, corresponding to a chemical efficiency of 20.4%. When the chlorine flow rate was reduced to 115 mmol/s, a higher chemical efficiency of 22.7% was attained. Measurements showed that the SPJSOG during normal operation could provide a singlet oxygen yield Y greater than or equal to 55%, a chlorine utilization U greater than or equal to 85%, and a relative water vapor concentration w = [H2O]/([O-2] + [Cl-2]) less than or equal to 0.1.

An analysis on overall crack-number-density of short-fatigue-cracks

The evolution of dispersed short-fatigue-cracks is analysed based on the equilibrium of crack-number-density (CND). By separating the mean value and the stochastic fluctuation of local CND, the equilibrium equation of overall CND is derived. Comparing with the mean-field equilibrium equation, the equilibrium equation of overall CND has different forms in the expression of crack-nucleation-rate or crack-growth-rate. The simulation results are compared with experimental measurements showing the stochastic analyses provide consistent tendency with experiments. The discrepancy in simulation results between overall CND and mean-field CND is discussed.

Generalization of Response Number for Dynamic Plastic Response of Shells Subjected to impulsive Loading

A dimensionless number, termed as response number in Zhao [Archive of Applied Mechanics 68 (1998) 524], has been suggested for the dynamic plastic response of beams and plates made up of rigidly perfect plastic materials subjected to dynamic loading. Many theoretical and experimental results can be reformulated into new concise forms with the response number. The concept of a new dimensionless number, response number, termed as Rn(n), is generalized in Zhao [Forschung im Ingenieurwesen 65 (1999) 107] to study the elastic, plastic, dynamic elastic as well as dynamic plastic buckling problems of columns, plates as well as shells. The response number Rn(n) is generalized to the dynamic behaviour of shells of various shapes in the present paper.

Application of response number for dynamic plastic response of plates subjected to impulsive loading

A dimensionless number, termed response number, is applied to the dynamic plastic response of plates subjected to dynamic loading. Many theoretical and experimental results presented by different researchers are reformulated into new concise forms with the response number. The advantage of the new forms is twofold: (1) they are more physically meaningful, and (2) they are independent of the choice of units, thus, they have wider range of applications.

Effect of liquid bridge volume on the instability in small-Prandtl-number half zones

A perturbation method is used to examine the linear instability of thermocapillary convection in a liquid bridge of floating half-zone filled with a small Prandtl number fluid. The influence of liquid bridge volume on critical Marangoni number and flow features is analyzed. The neutral modes show that the instability is mainly caused by the bulk flow that is driven by the nonuniform thermocapillary forces acting on the free surface. The hydrodynamic instability is dominant in the case of small Prandtl number fluid and the first instability mode is a stationary bifurcation. The azimuthal wave number for the most dangerous mode depends on the liquid bridge volume, and is not always two as in the case of a cylindrical liquid bridge with aspect ratio near 0.6. Its value may be equal to unity when the liquid bridge is relatively slender.

气相爆轰波在障碍物上Mach反射后流场的分析

详细分析了气相爆轰波在斜劈上Mach反射后的流场。根据胞格结构的烟迹记录，得到了爆轰波Mach反射时三波点迹线与楔面的夹角X与楔角θ的定量关系，依此推算出爆轰波绕射楔块 时从规则反射向Mach反射转变的临界角范围。论文也分析了影响X角和Mach杆后压力的主要 因素。该结果有助于更深入了解爆轰波的本质，也为数值模拟气相爆轰波在障碍物上Mach反射现象提供了可靠的对照依据。

The Mach reflection of a detonation based on soot track measurements

This paper presents a series of soot tracks formed by gaseous detonation waves diffracting around wedges with different wedge angles. These cellular structure patterns describe the Mach-reflection processes of a detonation and reveal some unique characteristics. They can be used to analyze the relationship between the trajectory angle of the triple point, wedge angle, and initial pressure in Mach reflection. Compared to the Mach-reflected one-dimensional shock wave in nonreactive air, all these unique characteristics for a Mach-reflected detonation should be attributed to the transverse-wave structure of the detonation front; meanwhile, the precursor shock wave and transverse wave influence the Mach-reflected detonation, respectively. The experimental results support the recently published numerical simulation of this complex phenomenon.

The Effect of Conical Flowfields on the Performance of Waveriders at Mach 6

The performance of 23 kinds of waveriders, derived from different conical flowfields, is analyzed by the numerical computation under the conditions of fight speed of Mach 6, attack angle of 0° and flight altitude of 30 km. These results indicate that the performance is influenced by the shapes and the width to height ratios (W/H ) of generating cones. The geometrical parameter and the lift coefficient are proportional to W/H, while the drag coefficient and the lift to drag ratio (L/D ) have extreme values. Considering the base drag and the computation errors, the waverider with the highest L/D is cut from the elliptical cone’s flowfield (W/H = 1.5―1.618), and the configuration with the lowest drag can also be obtained at W/H = 1:1.5. Accordingly, good suggestions are proposed for practical design based on these computational results.

气相爆轰波与空气冲击波在斜劈上Mach反射的特征参数比较

以传统的烟迹技术记录胞格结构，证实了气相爆轰波在管道中与斜劈相互作用发生Mach反射，得到了系统的实验结果。在此基础上详细地分析了爆轰波反射后流场的特征参数(三波点迹线与楔面的夹角X与斜劈楔角Θ及气体初压P_0的关系等),并与空气冲击波在斜劈上Mach反射后流场的特征参数进行了比较，初步探讨了参数偏差的原因。这一研究有助于了解爆轰波的前导冲击波和横波在Mach反射中所起的作用。

Suggestion of a new dimensionless number for dynamic plastic response of beams and plates

A dimensionless number, termed response number in the present paper, is suggested for the dynamic plastic response of beams and plates made of rigid-perfectly plastic materials subjected to dynamic loading. It is obtained at dimensional reduction of the basic governing equations of beams and plates. The number is defined as the product of the Johnson's damage number and the square of the half of the slenderness ratio for a beam; the product of the damage number and the square of the half of the aspect ratio for a plate or membrane loaded dynamically. Response number can also be considered as the ratio of the inertia force at the impulsive loading to the plastic limit load of the structure. Three aspects are reflected in this dimensionless number: the inertia of the applied dynamic loading, the resistance ability of the material to the deformation caused by the loading and the geometrical influence of the structure on the dynamic response. For an impulsively loaded beam or plate, the final dimensionless deflection is solely dependent upon the response number. When the secondary effects of finite deflections, strain-rate sensitivity or transverse shear are taken into account, the response number is as useful as in the case of simple bending theory. Finally, the number is not only suitable to idealized dynamic loads but also applicable to dynamic loads of general shape.

Towards an Understanding of the Influence of Sedimentation on Colloidal Aggregation by Peclet Number

The Peclet number is a useful index to estimate the importance of sedimentation as compared to the Brownian motion. However, how to choose the characteristic length scale for the Peclet number evaluation is rather critical because the diffusion length increases as the square root of the time whereas the drifting length is linearly related to time. Our Brownian dynamics simulation shows that the degree of sedimentation influence on the coagulation decreases when the dispersion volume fraction increases. Therefore using a fixed length, such as the diameter of particle, as the characteristic length scale for Peclet number evaluation is not a good choice when dealing with the influence of sedimentation on coagulation. The simulations demonstrated that environmental factors in the coagulation process, such as dispersion volume fraction and size distribution, should be taken into account for more reasonable evaluation of the sedimentation influence.

Instability from steady and axisymmetric to steady and asymmetric floating half zone convection in a fat liquid bridge of larger Prandtl number

The linear instability analysis of the present paper shows that the thermocapillary convection in a half floating zone of larger Prandtl number has a steady instability mode w(i) = 0 and m = 1 for a fat liquid bridge V = 1.2 with small geometrical aspect ratio A = 0.6. This conclusion is different from the usual idea of hydrothermal instability, and implies that the instability of the system may excite a steady and axial asymmetric state before the onset of oscillation in the ease of large Prandtl number.