The positive floating potential of high-temperature particles in plasma

It has been predicted that the floating potential of particles in plasma may become positive when the particle surface temperature is high enough, but, to our knowledge, no positive floating potential has been obtained yet. In the present paper the floating potential theory of high-temperature particles in plasma is developed to cover the positive potential range for the first time, and a general approximate analytical formula for the positive floating potential with a thin plasma sheath and subsonic plasma flow is derived from the new model recently proposed by the authors. The results show that when the floating potential is positive, the net flux of charge incident on the particle approaches a constant similar to the 'electron saturation' phenomena in the case of the electric probes.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.

Numerical study of bluff body flow structures

Our recent progress in numerical studies of bluff body flow structures and a new method for the numerical analysis of near wake flow field for high Reynolds number flow are introduced. The paper consists of three parts. In part one, the evolution of wake vortex structure and variation of forces on a flat plate in harmonic oscillatory flows and in in-line steady-harmonic combined flows are presented by an improved discrete vortex method, as the Keulegan-Carpenter number (KC) varies from 2 to 40 and ratios of U-m to U-0 are of O(10(-1)), O(10) and O(10), respectively. In part 2, a domain decomposition hybrid method, combining the finite-difference and vortex methods for numerical simulation of unsteady viscous separated flow around a bluff body, is introduced. By the new method, some high resolution numerical visualization on near wake evolution behind a circular cylinder at Re = 10(2), 10(3) and 3 x 10(3) are shown. In part 3, the mechanism and the dynamic process for the three-dimensional evolution of the Karman vortex and vortex filaments in braid regions as well as the early features of turbulent structure in the wake behind a circular cylinder are presented numerically by the vortex dynamics method.

Vortex motion in the early stages of unsteady flow around a circular cylinder

In this paper, processes in the early stages of vortex motion and the development of flow structure behind an impulsively-started circular cylinder at high Reynolds number are investigated by combining the discrete vortex model with boundary layer theory, considering the separation of incoming flow boundary layer and rear shear layer in the recirculating flow region. The development of flow structure and vortex motion, particularly the formation and development of secondary vortex and a pair of secondary vortices and their effect on the flow field are calculated. The results clearly show that the flow structure and vortices motion went through a series of complicated processes before the symmetric main vortices change into asymmetric: development of main vortices induces secondary vortices; growth of the secondary vortices causes the main vortex sheets to break off and causes the symmetric main vortices to become “free” vortices, while a pair of secondary vortices is formed; then the vortex sheets, after breaking off, gradually extend downstream and the structure of a pair of secondary vortices becomes relaxed. These features of vortex motion look very much like the observed features in some available flow field visualizations. The action of the secondary vortices causes the main vortex sheets to break off and converts the main vortices into free vortices. This should be the immediate cause leading to the instability of the motion of the symmetric main vortices. The flow field structure such as the separation position of boundary layer and rear shear layer, the unsteady pressure distributions and the drag coefficient are calculated. Comparison with other results or experiments is also made.

Investigation of High-Speed Combustible Gas Ignited by a Hot Gas Jet Produced in the Shock Tube

The high-speed combustible gas ignited by a hot gas jet, which is induced by shock focusing, was experimentally investigated. By use of the separation mode of shock tube, the test section of a single shock tube is split into two parts, which provide the high-speed flow of combustible gas and pilot flame of hot gas jet, respectively. In the interface of two parts of test sections the flame of jet was formed and spread to the high-speed combustible gas. Two kinds of the ignitions, 3-D "line-flame ignition" and 2-D "plane-flame ignition", were investigated. In the condition of 3-D "line-flame ignition" of combustion, thicker hot gas jet than pure air jet, was observed in schlieren photos. In the condition of 2-D "plane-flame ignition" of combustion, the delay time of ignition and the angle of flame front in schlieren photos were measured, from which the velocity of flame propagation in the high-speed combustible gas is estimated in the range of 30-90m/s and the delay time of ignition is estimated in the range of 0.12-0.29ms.

Chemical Nonequilibrium Effects on Flow Field for Reusable Launch Vehicles

High temperature chemical non-equilibrium phenomena have a great effect on the flow field around a reentry vehicle. A set of three dimensional Navier-Stokes equations have been solved by implicit finite volume NND scheme. Both ideal gas viscous flow and chemical non-equilibrium flow are calculated for a spherical-cone at a small angle of attack. The results of the two flows have been compared and the effect of chemical non-equilibrium has been analyzed. The effect of wall material's properties, such as catalysis and radiation were studied. The results are in good agreement with the referenced paper.

Experimental Studies on Model Scramjet in a Mach 6 High Enthalpy Free-Jet Wind Tunnel

A side-wall compression scramjet model with different combustor geometries has been tested in a propulsion tunnel that typically provides the testing flow with Mach number of 5.8, total temperature of 1800K, total pressure of 4.5MPa and mass flow rate of 4kg/s. This kerosene-fueled scramjet model consists of a side-wall compression inlet, a combustor and a thrust nozzle. A strut was used to increase the contraction ratio and to inject fuels, as well as a mixing enhancement device. Several wall cavities were also employed for flame-holding. In order to shorten the ignition delay time of the kerosene fuel, a little amount of hydrogen was used as a pilot flame. The pressure along the combustor has an evident raise after ignition occurred. Consequently thrust was observed during the fuel-on period. However, the thrust was still less than the drag of the scramjet model. For this reason, the drag variation produced by different strut and cavities was tested. Typical results showed that the cavities do not influence the drag so much, but the length of the strut does.

A Multiscale-Domain Decomposition Method for High Reynolds Number Flows

The high Reynolds number flow contains a wide range of length and time scales, and the flow domain can be divided into several sub-domains with different characteristic scales. In some sub-domains, the viscosity dissipation scale can only be considered in a certain direction; in some sub-domains, the viscosity dissipation scales need to be considered in all directions; in some sub-domains, the viscosity dissipation scales are unnecessary to be considered at all. For laminar boundary layer region, the characteristic length scales in the streamwise and normal directions are L and L Re-1/ 2 , respectively. The characteristic length scale and the velocity scale in the outer region of the boundary layer are L and U, respectively. In the neighborhood region of the separated point, the length scale l<

Serrated Plastic Flow During Nanoindentation in Nd-Based Bulk Metallic Glasses

The microstructure of Nd_{60}Al_{10}Ni_{10}Cu_{20-x}Fex (x = 0, 5, 7, 10, 15, 20) alloys can change from homogeneous phase to a composite structure consisting of amorphous phase plus clusters or nanocrystals by adjusting the Fe content. The effect of microstructure on the plastic deformation behavior in this alloy system is studied by using nanoindentation. The alloys with homogeneous amorphous structure exhibit pronounced flow serrations during the loading process of nanoindentation. The addition of Fe changes the plastic deformation behavior remarkablely. No flow serration is observed in the alloys with high Fe content for the indentation depth of 500 nm. The mechanism for the change of plastic serrated flow behavior is discussed.

Numerical Study of Heat Transfer Mechanism in Turbulent Supercritical CO2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

Research progress on high-enthalpy and hypersonic flows

The research progress on high-enthalpy and hypersorlic flows having been achieved in the Institute of Mechanics, Chinese Academy of Sciences, is reported in this paper. The paper consists of three main parts: The first part is on the techniques to develop advanced hypersonic test facilities, in which the detonation-driven shock-reflected tunnel and the detonation-driven shock-expanded tube are introduced. The shock tunnel can be used for generating hypersonic flows of a Mach number ranging from 10 to 20, and the expansion tube is applicable to simulate the flows with a speed of 7 similar to 10km/s. The second part is dedicated to the shock tunnel nozzle flow diagnosis to examine properties of the hypersonic flows thus created. The third part is on experiments and numerical simulations. The experiments include measuring the aerodynamic pitching moment and heat transfer in hypersonic flows, and the numerical work reports nozzle flow simulations and flow non-equilibrium effects on the possible experiments that may be carried out on the above-mentioned hypersonic test facilities.

Numerical study on heavy rigid particle motion in a plane wake flow by spectral element method

Experimental particle dispersion patterns in a plane wake flow at a high Reynolds number have been predicted numerically by discrete vortex method (Phys. Fluids A 1992; 4:2244-2251; Int. J. Multiphase Flow 2000; 26:1583-1607). To address the particle motion at a moderate Reynolds number, spectral element method is employed to provide an instantaneous wake flow field for particle dynamics equations, which are solved to make a detail classification of the patterns in relation to the Stokes and Froude numbers. It is found that particle motion features only depend on the Stokes number at a high Froude number and depend on both numbers at a low Froude number. A ratio of the Stokes number to squared Froude number is introduced and threshold values of this parameter are evaluated that delineate the different regions of particle behavior. The parameter describes approximately the gravitational settling velocity divided by the characteristic velocity of wake flow. In order to present effects of particle density but preserve rigid sphere, hollow sphere particle dynamics in the plane wake flow is investigated. The evolution of hollow particle motion patterns for the increase of equivalent particle density corresponds to that of solid particle motion patterns for the decrease of particle size. Although the thresholds change a little, the parameter can still make a good qualitative classification of particle motion patterns as the inner diameter changes.

Acoustic Calculation for Supersonic Turbulent Boundary Layer Flow

An approach which combines direct numerical simulation (DNS) with the Lighthill acoustic analogy theory is used to study the potential noise sources during the transition process of a Mach 2.25 flat plate boundary layer. The quadrupole sound sources due to the flow fluctuations and the dipole sound sources due to the fluctuating surface stress are obtained. Numerical results suggest that formation of the high shear layers leads to a dramatic amplification of amplitude of the fluctuating quadrupole sound sources. Compared with the quadrupole sound source, the energy of dipole sound source is concentrated in the relatively low frequency range.

Fluid flow induced calcium response in bone cell network

In our previous work, bone cell networks with controlled spacing and functional intercellular gap junctions had been successfully established by using microcontact printing and self assembled monolayers technologies [Guo, X. E., E. Takai, X. Jiang, Q. Xu, G. M. Whitesides, J. T. Yardley, C. T. Hung, E. M. Chow, T. Hantschel, and K. D. Costa. Mol. Cell. Biomech. 3:95-107, 2006]. The present study investigated the calcium response and the underlying signaling pathways in patterned bone cell networks exposed to a steady fluid flow. The glass slides with cell networks were separated into eight groups for treatment with specific pharmacological agents that inhibit pathways significant in bone cell calcium signaling. The calcium transients of the network were recorded and quantitatively evaluated with a set of network parameters. The results showed that 18 alpha-GA (gap junction blocker), suramin (ATP inhibitor), and thapsigargin (depleting intracellular calcium stores) significantly reduced the occurrence of multiple calcium peaks, which were visually obvious in the untreated group. The number of responsive peaks also decreased slightly yet significantly when either the COX-2/PGE(2) or the NOS/nitric oxide pathway was disrupted. Different from all other groups, cells treated with 18 alpha-GA maintained a high concentration of intracellular calcium following the first peak. In the absence of calcium in the culture medium, the intracellular calcium concentration decreased slowly with fluid flow without any calcium transients observed. These findings have identified important factors in the flow mediated calcium signaling of bone cells within a patterned network.

Two-phase flow and pool boiling heat transfer in microgravity

Researches on two-phase flow and pool boiling heat transfer in microgravity, which included groundbased tests, flight experiments, and theoretical analyses, were conducted in the National Microgravity Laboratory/CAS. A semi-theoretical Weber number model was proposed to predict the slug-to-annular flow transition of two-phase gas–liquid flows in microgravity, while the influence of the initial bubble size on the bubble-to-slug flow transition was investigated numerically using the Monte Carlo method. Two-phase flow pattern maps in microgravity were obtained in the experiments both aboard the Russian space station Mir and aboard IL-76 reduced gravity airplane. Mini-scale modeling was also used to simulate the behavior of microgravity two-phase flow on the ground. Pressure drops of two-phase flow in microgravity were also measured experimentally and correlated successfully based on its characteristics. Two space experiments on pool boiling phenomena in microgravity were performed aboard the Chinese recoverable satellites. Steady pool boiling of R113 on a thin wire with a temperature-controlled heating method was studied aboard RS-22, while quasi-steady pool boiling of FC-72 on a plate was studied aboard SJ-8. Ground-based experiments were also performed both in normal gravity and in short-term microgravity in the drop tower Beijing. Only slight enhancement of heat transfer was observed in the wire case, while enhancement in low heat flux and deterioration in high heat flux were observed in the plate case. Lateral motions of vapor bubbles were observed before their departure in microgravity. The relationship between bubble behavior and heat transfer on plate was analyzed. A semi-theoretical model was also proposed for predicting the bubble departure diameter during pool boiling on wires. The results obtained here are intended to become a powerful aid for further investigation in the present discipline and development of two-phase systems for space applications.