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.

Evolution of ideal micro-crack system

The ideal micro-cracks are treated with the number-density function; the characteristics of their evolution are investigated; a deterministic model is applied to the discussion of their extension. It is discovered that under certain conditions saturation may occur in the number-density. The main features of the statistical formulation are illustrated by several examples and compared with those observed in experiments.

Excitation mechanism of thermocapillary oscillatory convection

It is suggested that the oscillation of thermocapillary convection may be excited by the buoyancy instability. By means of numerical simulation of the finite-element method, the temperature distributions in the liquid bridge are qualitatively analyzed. The temperature gradient in a certain flow region of liquid bridge may turn to be parallel to the direction of gravity when the temperature difference △T between two boundary rods of liquid bridge is larger than the critical value. The buoyancy instability may be excited, and then the thermocapillary oscillatory convection appears, as the temperature difference increases further. The distribution of the critical Marangoni number in the micro-gravity environment is derived from the data on the ground experiments. The results show that the onset of thermocapillary oscillatory convection is delayed in the case of smaller typical scale of liquid bridge and lower gravity environment.

Influence of suspended particles and particle cover over the density interface

The effect of the particle cover over the density interface between two layers of fluids and of the suspended solid particles in the upper turbulcnt layer on the turbulent entrainment has been studied experimentally. The entrainment distance D is a function of the time of power: D=kt, where =0.200-0.130p. For suspended particles in the upper layer and pure 2-layer fluid is equal to 0.200, but the value of k for the suspended particles is smaller than that for the pure 2-layer fluid. The non-dimensional entrainment velocity is E=KRiln, where n=1.50+0.93 p. It is shown that the particle cover over the interface changes the power of Ril in the entrainment and hinders the turbulent entrainment. The variation rule of E for the suspended particles is the same as that for the pure 2-layer fluid, but the K value of the former is smaller than that of the latter. The turbulent mixing mechanism has been discussed.

Stratification by solid particle suspensions

Turbulence was generated by an oscillating grid above a bed of sediment of spherical glass beads. As expected, part of the sediment was lifted up by the grid action and a suspension layer of depth D formed above the grid. This depth was found remaining independent of grid action but varying with the sediment layer depth when the grid action was kept constant. Volume concentration measurements show the existence of only weak concentration gradients over the layer depth with a rapid fall off in concentration at the outer edge. The theoretical analysis based on a concentration flux model is in good qualitative agreement with observations.

Magnetostatic equilibrium of isolated magnetic-flux tube

In the present paper, an isolated axisymmetric flux tube is discussed for slender magnetic configuration. The magnetostatic model and the stratified atmospheric model are applied, respectively, to the regions inside and outside the flux tube. The problem is described mathematically by the nonlinear partial differential equations under the nonlinear boundary condition at the free boundary of flux tube. According to the approximation of a small expansive angle, the solutions of series expressions are obtained formally. The model of polytropic plasma is discussed in detail especially. The results show the distributions of thermodynamic quantities and magnetic field extending from the high β region to the low β region, and the flux tube may be either divergent or convergent according to the pressure difference outside and inside the flux tube.

Approximate theory of natural-convection with small grashof number

To gain some insight into the behaviour of low-gravity flows in the material processing in space, an approximate theory has been developed for the convective motion of fluids with a small Grashof number Gr. The expansion of the variables into a series of Gr reduces the Boussinesq equation to a system of weakly coupled linearly inhomogeneous equations. Moreover, the analogy concept is proposed and utilized in the study of the plate bending problems in solid mechanics. Two examples are investigated in detail, i. e. the 2-dimensional steady flows in either circular or square infinite closed cylinder, which is horizontally imposed at a specified temperature of linear distribution on the boundaries. The results for stream function ψ, velocity u and temperature T are provided. The analysis of the influences of some parameters such as the Grashof number Gr and the Prandtl number Pr, on motions will lead to several interesting conclusions. The theory seems to be useful for seeking for an analytical solutions. At least, it will greatly simplify the complicated problems originally governed by the Navier-Stokes equation including buoyancy. It is our hope that the theory might be applicable to unsteady or 3-dimensional cases in future.

Fracture analysis for 3-dimensional bodies with a surface crack

This paper deals with fracture analyses in 3-dimensional bodies containing a surface crack. A general solution of stress-strain fields at crack tip is proposed. Based on the stress-strain fields obtained, a high-order 3-dimensional special element is established to calculate the stress intensity factors in a plate with a surface crack. The variation of stress intensity factors with geometric parameters is investigated.

Hierarchal structure of the simplified Navier-Stokes equations (SNSE) and its mechanical connotation and application

The hierarchial structure and mathematical property of the simplified Navier-Stokesequations (SNSE) are studied for viscous flow over a sphere and a jet of compressible flu-id. All kinds of the hierarchial SNSE can be divided into three types according to theirmathematical property and also into five groups according to their physical content. Amultilayers structure model for viscous shear flow with a main stream direction is pre-sented. For the example of viscous incompressible flow over a flat plate there existthree layers for both the separated flow and the attached flow; the character of thetransition from the three layers of attached flow to those of separated flow is elucidated.A concept of transition layer being situated between the viscous layer and inviscidlayer is introduced. The transition layer features the interaction between viscous flow andinviscid flow. The inner-outer-layers-matched SNSE proposed by the present author inthe past is developed into the layers matched (LsM)-SNSE.

The elasto-damage theory of the components assembling model

The potential energy in materials is well approximated by pair functional which is composed of pair potentials and embedding energy. During calculating material potential energy, the orientational component and the volumetric component are derived respectively from pair potentials and embedding energy. The sum of energy of all these two kinds of components is the material potential. No matter how microstructures change, damage or fracture, at the most level, they are all the changing and breaking atomic bonds. As an abstract of atomic bonds, these components change their stiffness during damaging. Material constitutive equations have been formulated by means of assembling all components' response functions. This material model is called the component assembling model. Theoretical analysis and numerical computing indicate that the proposed model has the capacity of reproducing some results satisfactorily, with the advantages of great conceptual simplicity, physical explicitness, and intrinsic induced anisotropy, etc.

Large-eddy simulation of flows past a flapping airfoil using immersed boundary method

The numerical simulation of flows past flapping foils at moderate Reynolds numbers presents two challenges to computational fluid dynamics: turbulent flows and moving boundaries. The direct forcing immersed boundary (IB) method has been devel- oped to simulate laminar flows. However, its performance in simulating turbulent flows and transitional flows with moving boundaries has not been fully evaluated. In the present work, we use the IB method to simulate fully developed turbulent channel flows and transitional flows past a stationary/plunging SD7003 airfoil. To suppress the non-physical force oscillations in the plunging case, we use the smoothed discrete delta function for interpolation in the IB method. The results of the present work demonstrate that the IB method can be used to simulate turbulent flows and transitional flows with moving boundaries.

New Numerical Algorithms in SUPER CE/SE and Their Applications in Explosion Mechanics

The new numerical algorithms in SUPER/CESE and their applications in explosion mechanics are studied. The researched algorithms and models include an improved CE/SE (space-time Conservation Element and Solution Element) method, a local hybrid particle level set method, three chemical reaction models and a two-fluid model. Problems of shock wave reflection over wedges, explosive welding, cellular structure of gaseous detonations and two-phase detonations in the gas-droplet system are simulated by using the above-mentioned algorithms and models. The numerical results reveal that the adopted algorithms have many advantages such as high numerical accuracy, wide application field and good compatibility. The numerical algorithms presented in this paper may be applied to the numerical research of explosion mechanics.

Numerical Study On Transient Flow In The Deep Naturally Fractured Reservoir With High Pressure

According to the experimental results and the characteristics of the pressure-sensitive fractured formation, a transient flow model is developed for the deep naturally-fractured reservoirs with different outer boundary conditions. The finite element equations for the model are derived. After generating the unstructured grids in the solution regions, the finite element method is used to calculate the pressure type curves for the pressure-sensitive fractured reservoir with different outer boundaries, such as the infinite boundary, circle boundary and combined linear boundaries, and the characteristics of the type curves are comparatively analyzed. The effects on the pressure curves caused by pressure sensitivity module and the effective radius combined parameter are determined, and the method for calculating the pressure-sensitive reservoir parameters is introduced. By analyzing the real field case in the high temperature and pressure reservoir, the perfect results show that the transient flow model for the pressure-sensitive fractured reservoir in this paper is correct.

Control of lunar and martian dust - Experimental insights from artificial and natural cyanobacterial and algal crusts in the desert of Inner Mongolia, China

Studies on the colonization of environmentally extreme ground surfaces were conducted in a Mars-like desert area of Inner Mongolia, People's Republic of China, with microalgae and cyanobacteria. We collected and mass-cultured cyanobacterial strains from these regions and investigated their ability to form desert crusts artificially. These crusts had the capacity to resist sand wind erosion after just 15 days of growth. Similar to the surface of some Chinese deserts, the surface of Mars is characterized by a layer of fine dust, which will challenge future human exploration activities, particularly in confined spaces that will include greenhouses and habitats. We discuss the use of such crusts for the local control of desert sands in enclosed spaces on Mars. These experiments suggest innovative new directions in the applied use of microbe-mineral interactions to advance the human exploration and settlement of space.

A simple closed aquatic ecosystem (CAES) for space

A closed aquatic ecosystem (CAES) was developed to stud), the effects of microgravity on the function of closed ecosystems aboard the Chinese retrieved satellite and on the spacecraft SHENZHOU-II. These systems housed a small freshwater snail (Bulinus australianus) and an autotrophic green algae (Chlorella pyrenoidosa). The results of the test on the satellite were that the concentration of algae changed little, but that the snails died during the experiments. We then sought to optimize the function of the control system, the cultural conditions and the data acquisition system and carried out an experiment on the spacecraft SHENZHOU-II. Using various sensors to monitor the CAES, real-time data regarding the operation of the CAES in microgravity was acquired. In addition, all on-board Ig centrifuge was included to identify gravity-related factors. It was found that microgravity is the major factor affecting the operation of the CAES in space. The change in biomass of the primary producer during each day in microgravity was larger than that of the control groups. The mean biomass concentration per day in the microgravity group decreased, but that of the control groups increased for several days and then leveled off. Space effects on the biomass of a primary producer may be a result of microgravity effects leading to increasing metabolic rates of the consumer combined with decreases in photosynthesis. (c) 2007 COSPAR. Published by Elsevier Ltd. All rights reserved.