Numerical study of bifurcation solutions of spherical Taylor-Couette flow

The steady bifurcation flows in a spherical gap (gap ratio sigma=0.18) with rotating inner and stationary outer spheres are simulated numerically for Re(c1)less than or equal to Re less than or equal to 1 500 by solving steady axisymmetric incompressible Navier-Stokes equations using a finite difference method. The simulation shows that there exist two steady stable flows with 1 or 2 vortices per hemisphere for 775 less than or equal to Re less than or equal to 1 220 and three steady stable flows with 0, 1, or 2 vortices for 1 220conditions which on be generated by different accelerations of the inner sphere. Generation of zero-or two-vortex flow depends mainly on the acceleratio n, but that of one-vortex flow also depends on the perturbation breaking the equatorial symmetry. The mechanism of development of a saddle point in the meridional plane at higher Re number and its role in the formation of two-vortex flow are analyzed.

Water and salt movements in simultaneous flood-irrigation and well-drainage operations

This paper describes a new technology for solonchak soil reclamation in which surface flood irrigation of fresh water and pumped wells drainage of salty groundwater are combined. The comprehensive investigation of water and salt movement has been conducted through field test, laboratory simulation and numerical calculation. The dependence of desalination on irrigation water quantity, drainage quantity, leaching time and other parameters is obtained based on the field tests. The entire desalination process under the flood-irrigation and well-drainage operations was experimentally simulated in a vertical soil column. The water and salt movement has been numerically analysed for both the field and laboratory conditions. The present work indicates that this new technology can greatly improve the effects of desalination.

Preliminary-results of gaas single-crystal growth under high gravity conditions

GaAs single crystals have been grown under high gravity conditions, up to 9g0, by a recrystallization method with decreasing temperature. The impurity striations in GaAs grown under high gravity become weak and indistinct with smaller striation spacings. The dislocation density of surcharge-grown GaAs increases with increase of centrifugal force. The cathodoluminescence results also show worse perfection in the GaAs grown at high gravity than at normal earth gravity.

Complex stress functions and plastic zone sizes for notch cracks subjected to various loading conditions

In this paper, the conformal mapping method is used to solve the plane problem of an infinite plate containing a central lip-shaped notch subjected to biaxial loading at a remote boundary or a surface uniform pressure on the notch. The stress intensity factors KI and KII are obtained by the derived complex stress functions. The simple analytical expressions can be applied to the situation of cracks originating from a circular or an elliptical notch. The plastic zone sizes for such notch cracks are subsequently evaluated in light of the Dugdale strip yield concept. The results are consistent with available numerical data.

Information Preservation (IP) Method in Simulation

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Modified Simplified Rate Equation Model of Flowing Chemical Oxygen-Iodine Laser and its Application

A modified simplified rate equation (RE) model of flowing chemical oxygen-iodine laser (COIL), which is adapted to both the condition of homogeneous broadening and inhomogeneous broadening being of importance and the condition of inhomogeneous broadening being predominant, is presented for performance analyses of a COIL. By using the Voigt profile function and the gain-equal-loss approximation, a gain expression has been deduced from the rate equations of upper and lower level laser species. This gain expression is adapted to the conditions of very low gas pressure up to quite high pressure and can deal with the condition of lasing frequency being not equal to the central one of spectral profile. The expressions of output power and extraction efficiency in a flowing COIL can be obtained by solving the coupling equations of the deduced gain expression and the energy equation which expresses the complete transformation of the energy stored in singlet delta state oxygen into laser energy. By using these expressions, the RotoCOIL experiment is simulated, and obtained results agree well with experiment data. Effects of various adjustable parameters on the performances of COIL are also presented.

Information preservation (IP) method in simulation of internal rarefied gas flows in MEMS

This paper reviews firstly methods for treating low speed rarefied gas flows: the linearised Boltzmann equation, the Lattice Boltzmann method (LBM), the Navier-Stokes equation plus slip boundary conditions and the DSMC method, and discusses the difficulties in simulating low speed transitional MEMS flows, especially the internal flows. In particular, the present version of the LBM is shown unfeasible for simulation of MEMS flow in transitional regime. The information preservation (IP) method overcomes the difficulty of the statistical simulation caused by the small information to noise ratio for low speed flows by preserving the average information of the enormous number of molecules a simulated molecule represents. A kind of validation of the method is given in this paper. The specificities of the internal flows in MEMS, i.e. the low speed and the large length to width ratio, result in the problem of elliptic nature of the necessity to regulate the inlet and outlet boundary conditions that influence each other. Through the example of the IP calculation of the microchannel (thousands m ? long) flow it is shown that the adoption of the conservative scheme of the mass conservation equation and the super relaxation method resolves this problem successfully. With employment of the same measures the IP method solves the thin film air bearing problem in transitional regime for authentic hard disc write/read head length ( 1000 L m ? = ) and provides pressure distribution in full agreement with the generalized Reynolds equation, while before this the DSMC check of the validity of the Reynolds equation was done only for short ( 5 L m ? = ) drive head. The author suggests degenerate the Reynolds equation to solve the microchannel flow problem in transitional regime, thus provides a means with merit of strict kinetic theory for testing various methods intending to treat the internal MEMS flows.

Discontinuous growth mechanisms of mode II crack under high-speed impact conditions

A recoverable plate impact testing technology has been used for studying the growth mechanisms of mode II crack. The results show that interactions of microcracks ahead of a crack tip cause the crack growth unsteadily. Failure mode transitions of materials were observed. Based on the observations, a discontinuous crack growth model was established. Analysis shows that the shear crack grows unsteady as the growth speed is between the Rayleigh wave speed c(R) and the shear wave speed c(s); however, when the growth speed approaches root 2c(s), the crack grows steadily. The transient microcrack growth makes the main crack speed to jump from subsonic to intersonic and the steady growth of all the sub-cracks leads the main crack to grow stably at an intersonic speed.

Simulated pathogenic conformational switch regions matched well with the biochemical findings

Pathogenic conformational conversion is a general causation of many disease, such as transmissible spon- giform encephalopathy (TSE) caused by misfolding of prion, sickle cell anemia, and etc. In such structural changes, misfolding occurs in regions important for the stability of native structure firstly. This destabi- lizes the normal conformation and leads to subsequent errors in folding pathway. Sites involved in the first stage can be deemed switch regions of the protein, and are vital for conformational conversion. Namely it could be a switch of disease at residue level. Here we report an algorithm that can identify such sites computationally with an accuracy of 93%, by calculating the probability of the native structure of a short segment jumping to a mistake one. Knowledge of such switch sites could be used to target clinical therapy, study physiological and pathologic mechanism of protein, and etc.

An Improved Ce/Se Scheme For Multi-Material Elastic-Plastic Flows And Its Applications

In this paper, a new definition of SE and CE, which is based on the hexahedron mesh and simpler than Chang's original CE/SE method (the space-time Conservation Element and Solution Element method), is proposed and an improved CE/SE scheme is constructed. Furthermore, the improved CE/SE scheme is extended in order to solve the elastic-plastic flow problems. The hybrid particle level set method is used for tracing the interfaces of materials. Proper boundary conditions are presented in interface tracking. Two high-velocity impact problems are simulated numerically and the computational results are carefully compared with the experimental data, as well as the results from other literature and LS-DYNA software. The comparisons show that the computational scheme developed currently is clear in physical concept, easy to be implemented and high accurate and efficient for the problems considered. (C) 2008 Elsevier Ltd. All rights reserved.

Catalytic Cracking and Heat Sink Capacity of Aviation Kerosene Under Supercritical Conditions

Catalytic cracking of China no. 3 aviation kerosene using a zeolite catalyst was investigated under supercritical conditions. A three-stage heating/cracking system was specially designed to be capable of heating 0.8 kg kerosene to a temperature of 1050 K and pressure of 7.0 MPa with maximum mass flow rate of 80 g/s. Sonic nozzles of different diameters were used to calibrate and monitor the mass flow rate of the cracked fuel mixture. With proper experiment arrangements, the mass flow rate per unit throat area of the cracked fuel mixture was found to well correlate with the extent of fuel conversion. The gaseous products obtained from fuel cracking under different conditions were also analyzed using gas chromatography. Composition analysis showed that the average molecular weight of the resulting gaseous products and the fuel mass conversion percentage were a strong function of the fuel temperature and were only slightly affected by the fuel pressure. The fuel conversion was also shown to depend on the fuel residence time in the reactor, as expected. Furthermore, the heat sink levels due to sensible heating and endothermic cracking were determined and compared at varying test conditions. It was found that at a fuel temperature of similar to 1050 K, the total heat sink reached similar to 3.4 MJ/kg, in which chemical heat sink accounted for similar to 1.5 MJ/kg.

Heat Transfer of Aviation Kerosene at Supercritical Conditions

The heat transfer characteristics of China no. 3 kerosene were investigated experimentally and analytically under conditions relevant to a regenerative cooling system for scramjet applications. A test facility developed for the present study can handle kerosene in a temperature range of 300-1000 K, a pressure range of 2.6-5 MPa, and a mass How rate range of 10-100 g/s. In addition, the test section was uniquely designed such that both the wall temperature and the bulk fuel temperature were measured at the same location along the flowpath. The measured temperature distributions were then used to analytically deduce the local heat transfer characteristics. A 10-component kerosene surrogate was proposed and employed to calculate the fuel thermodynamic and transport properties that were required in the heat transfer analysis. Results revealed drastic changes in the fuel flow properties and heat transfer characteristics when kerosene approached its critical state. Convective heat transfer enhancement was also found as kerosene became supercritical. The heat transfer correlation in the relatively low-fuel-temperature region yielded a similar result to other commonly used jet fuels, such as JP-7 and JP-8, at compressed liquid states. In the high-fuel-temperature region, near and beyond the critical temperature, heat transfer enhancement was observed; hence, the associated correlation showed a more significant Reynolds number dependency.

Configuration optimization of hypersonic vehicles under transitional flow conditions

The rarefied gas effects on several configurations are investigated under hypersonic flow conditions using the direct simulation Mont Carlo method. It is found that the Knudsen number, the Mach number, and the angle of attack all play a mixed role in the aerodynamics of a flat plate. The ratio of lift to drag decreases as the Knudsen number increases. Studies on 3D delta wings show that the ratio of lift to drag could be increased by decreasing the wing thickness and/or by increasing the wing span. It is also found that the waveriders could produce larger ratio of lift to drag as compared with the delta wing having the same length, wing span, and cross section area.