Examination of the LBM in simulation of microchannel flow in transitional regime

The gas flows in micro-electro-mechanical systems possess relatively large Knudsen number and usually belong to the slip flow and transitional flow regimes. Recently the lattice Boltzmann method (LBM) was proposed by Nie et al. in Journal of Statistical Physics, vol. 107, pp. 279-289, in 2002 to simulate the microchannel and microcavity flows in the transitional flow regime. The present article intends to test the feasibility of doing so. The results of using the lattice Boltzmann method and the direct simulation Monte Carlo method show good agreement between them for small Kn (Kn = 0.0194), poor agreement for Kn = 0.194, and large deviation for Kn = 0.388 in simulating microchannel flows. This suggests that the present version of the lattice Boltzmann method is not feasible to simulate the transitional channel flow.

A geometry model for tortuosity of flow path in porous media

A simple geometry model for tortuosity of flow path in porous media is proposed based on the assumption that some particles in a porous medium are unrestrictedly overlapped and the others are not. The proposed model is expressed as a function of porosity and there is no empirical constant in this model. The model predictions are compared with those from available correlations obtained numerically and experimentally, both of which are in agreement with each other. The present model can also give the tortuosity with a good approximation near the percolation threshold. The validity of the present tortuosity model is thus verified.

Two-phase flow patterns in a 90 degrees bend at microgravity

Bends are widely used in pipelines carrying single- and two-phase fluids in both ground and space applications. In particular, they play more important role in space applications due to the extreme spatial constraints. In the present study, a set of experimental data of two-phase flow patterns and their transitions in a 90degrees bend with inner diameter of 12.7 mm. and curvature radius of 76.5 mm at microgravity conditions are reported. Gas and liquid superficial velocities are found to range from (1.0 similar to 23.6) m/s for gas and (0.09 similar to 0.5) m/s for liquid, respectively. Three major flow patterns, namely slug, slug-annular transitional, and annular flows, are observed in this study. Focusing on the differences between flow patterns in bends and their counterparts in straight pipes, detailed analyses of their characteristics are made. The transitions between adjoining flow patterns are found to be more or less the same as those in straight pipes, and can be predicted using Weber number models satisfactorily. The reasons for such agreement are carefully examined.

Modeling, numerical methods, and simulation for particle-fluid two-phase flow problems

In this paper, we study the issues of modeling, numerical methods, and simulation with comparison to experimental data for the particle-fluid two-phase flow problem involving a solid-liquid mixed medium. The physical situation being considered is a pulsed liquid fluidized bed. The mathematical model is based on the assumption of one-dimensional flows, incompressible in both particle and fluid phases, equal particle diameters, and the wall friction force on both phases being ignored. The model consists of a set of coupled differential equations describing the conservation of mass and momentum in both phases with coupling and interaction between the two phases. We demonstrate conditions under which the system is either mathematically well posed or ill posed. We consider the general model with additional physical viscosities and/or additional virtual mass forces, both of which stabilize the system. Two numerical methods, one of them is first-order accurate and the other fifth-order accurate, are used to solve the models. A change of variable technique effectively handles the changing domain and boundary conditions. The numerical methods are demonstrated to be stable and convergent through careful numerical experiments. Simulation results for realistic pulsed liquid fluidized bed are provided and compared with experimental data. (C) 2004 Elsevier Ltd. All rights reserved.

Numerical estimate of the stability curve for the flow past a rotating cylinder

It is demonstrated that the primary instability of the wake of a two-dimensional circular cylinder rotating with constant angular velocity can be qualitatively well described by the Landau equation. The coefficients of the Landau equation are determined by means of numerical simulations for the Navier-Stokes equations. The critical Reynolds numbers, which depend on the angular velocity of the cylinder, are evaluated correctly by linear regression. (C) 2004 American Institute of Physics.

Rheological flow from a die and painting on a moving solid wall

Die swell is an important, phenomenon. in polymer processing, and is explained usually by rheological properties of the fluid. Because of the nonuniform of temperature distribution on the free surface of the liquid jet, the thermo capillary convection driven by surface tension gradient exists. The rheological fluid flowing out of a die and painting on a moving solid wall is studied by the numerical finite element method of a two-dimensional and unsteady model in the present paper, and both the rheological effect of a non-Newtonian fluid and the thermocapillary effect are considered. The results show that both,effects; will enlarge the cross-section of the fluid jet, and the rheological effect of non-Newtonian fluid dominates the process in general.

Numerical Study on Flow Field and Temperature Distribution in Growth Process of 200 mm Czochralski Silicon Crystals

The melt flow and temperature distribution in a 200 mm silicon Czochralski furnace with a cusp magnetic field was modeled and simulated by using a finite-volume based FLUTRAPP ( Fluid Flow and Transport Phenomena Program) code. The melt flow in the crucible was focused, which is a result of the competition of buoyancy, the centrifugal forces caused by the rotations of the crucible and crystal, the thermocapillary force on the free surfaces and the Lorentz force induced by the cusp magnetic field. The zonal method for radiative heat transfer was used in the growth chamber, which was confined by the crystal surface, melt surface, crucible, heat shield, and pull chamber. It was found that the cusp magnetic field could strength the dominant counter-rotating swirling flow cell in the crucible and reduce the flow oscillation and the pulling-rate fluctuation. The fluctuation of dopant and oxygen concentration in the growing crystal could thus be smoothed.

Critical serrated flow features during nanoindentation in La-based bulk metallic glasses

La57.6Al17.5(Cu,Ni)(24.9) and La64Al14(Cu,Ni)(22) bulk metallic glasses (BMGs) were prepared by copper-mould casting method. Plastic deformation behavior of the two BMGs at various loading rates was studied by nanoindentation. The results showed that the La57.6Al17.5(Cu,Ni)(24.9) BMG with a glass transition temperature of 423 K exhibited prominent serrated flow at low loading rates, whereas less pronounced serrated flow at high rates during nanoindentation. In contrast, the La64Al14(Cu,Ni)(22) BMG with a glass transition temperature of 401 K exhibited prominent serrated flow at high loading rates. The different rate dependency of serrated flow in the two La-based BMGs is related to the different glass transition temperature, and consequently the degree of viscous flow during indentation at room temperature. A smoother flow occurs in the alloy with relatively lower glass transition temperature, due to the relaxation of stress concentration.

Proper orthogonal decomposition of oscillatory Marangoni flow in half-zone liquid bridges of low-Pr fluids

Proper orthogonal decomposition (POD) using method of snapshots was performed on three different types of oscillatory Marangoni flows in half-zone liquid bridges of low-Pr fluid (Pr = 0.01). For each oscillation type, a series of characteristic modes (eigenfunctions) have been extracted from the velocity and temperature disturbances, and the POD provided spatial structures of the eigenfunctions, their oscillation frequencies, amplitudes, and phase shifts between them. The present analyses revealed the common features of the characteristic modes for different oscillation modes: four major velocity eigenfunctions captured more than 99% of the velocity fluctuation energy form two pairs, one of which is the most energetic. Different from the velocity disturbance, one of the major temperature eigenfunctions makes the dominant contribution to the temperature fluctuation energy. On the other hand, within the most energetic velocity eigenfuction pair, the two eigenfunctions have similar spatial structures and were tightly coupled to oscillate with the same frequency, and it was determined that the spatial structures and phase shifts of the eigenfunctions produced the different oscillatory disturbances. The interaction of other major modes only enriches the secondary spatio-temporal structures of the oscillatory disturbances. Moreover, the present analyses imply that the oscillatory disturbance, which is hydrodynamic in nature, primarily originates from the interior of the liquid bridge. (C) 2007 Elsevier B.V. All rights reserved.

Study on flow characteristics of solid/liquid system in lysozyme crystal growth

During the process of lysozyme protein crystallization with batch method, the macroscopic flow field of solid/liquid system was observed by particle image velocimetry (PIV). Furthermore, a normal growth rate of (110) face and local flow field around a single protein crystal were obtained by a long work distance microscope. The experimental results showed that the average velocity, the maximal velocity of macroscopic solid/liquid system and the velocity of local flow field around single protein crystal were fluctuant. The effective boundary layer thickness delta(eff), the concentration at the interface Q and the characteristic velocity V were calculated using a convection-diffusion model. The results showed that the growth of lysozyme crystal in this experiment was dominated by interfacial kinetics rather than bulk transport, and the function of buoyancy-driven flow in bulk transport was small, however, the effect of bulk transport in crystal growth had a tendency to increase with the increase of lysozyme concentration. The calculated results, also showed that the order of magnitude of shear force was about 10(-21) N, which was much less than the bond force between the lysozyme molecules. Therefore the shear force induced by buoyancy-driven flows cannot remove the protein molecules from the interface of crystal.

The characterization of plastic flow in three different bulk metallic glass systems

Plastic deformation behaviors of Zr52.5Al10Ni10Cu15Be12.5, Mg65Cu25Gd10 and Pd43Ni10Cu27P20 bulk metallic glasses (BMGs) are studied by using the depth-sensing nanoindentation, macroindentation and uniaxial compression. The significant difference in plastic deformation behavior cannot be correlated to the Poisson's ratio or the ratio of shear modulus to bulk modulus of the three BMGs, but can be explained by the free volume model. It is shown that the nucleation of local shear band is easy and multiple shear bands can be activated in the Zr52.5Al10Ni10Cu15Be12.5 alloy, which exhibits a distinct plastic strain during uniaxial compression and less serrated flow during nanoindentation. (c) 2006 Elsevier B.V. All rights reserved.

Attractor crisis and bursting in a fluid flow with two no-slip directions

Characteristic burtsing behavior is observed in a driven, two-dimensional viscous flow, confined to a square domain and subject to no-slip boundaries. Passing a critical parameter value, an existing chaotic attractor undergoes a crisis, after which the flow initially enters a transient bursting regime. Bursting is caused by ejections from and return to a limited subdomain of the phase space, whereas the precrisis chaotic set forms the asymptotic attractor of the flow. For increasing values of the control parameter the length of the bursting regime increases progressively. Passing another critical parameter value, a second crisis leads to the appearance of a secondary type of bursting, of very large dynamical range. Within the bursting regime the flow then switches in irregular intervals from the primary to the secondary type of bursting. Peak enstrophy levels for both types of bursting are associated to the collapse of a primary vortex into a quadrupolar state.