VISUALIZATION RESULTS OF AXISYMMETRICAL WAKE FLOW WITH/WITHOUT ACOUSTIC EXCITATION

Visualization results demonstrate the evolution of Kelvin-Helmholtz unstable waves into vortex pairing in a separated shear layer of a blunf circular. The results with acoustic excitation are quite different from that without acoustic excitation, and the phenomenon with excitation in a separated shear layer follows the rule of Devil s staircase, which always occurs in a non-linear dynamical system of two coupling vibrators.

MILU-CG method and the numerical study on the flow around a rotating circular cylinder

A hybrid finite difference method and vortex method (HDV), which is based on domain decomposition and proposed by the authors (1992), is improved by using a modified incomplete LU decomposition conjugate gradient method (MILU-CG), and a high order implicit difference algorithm. The flow around a rotating circular cylinder at Reynolds number R-e = 1000, 200 and the angular to rectilinear speed ratio alpha is an element of (0.5, 3.25) is studied numerically. The long-time full developed features about the variations of the vortex patterns in the wake, and drag, lift forces on the cylinder are given. The calculated streamline contours agreed well with the experimental visualized flow pictures. The existence of critical states and the vortex patterns at the states are given for the first time. The maximum lift to drag force ratio can be obtained nearby the critical states.

Rheological effect on thermocapillary flow of a liquid film jet painted on a moving boundary

In the present paper, a liquid (or melt) film of relatively high temperature ejected from a vessel and painted on the-moving solid film is analyzed by using the second-order fluid model of the non-Newtonian fluid. The thermocapillary flow driven by the temperature gradient on the free surface of a Newtonian liquid film was discussed before. The effect of rheological fluid on thermocapillary flow is considered in the present paper. The analysis is based on the approximations of lubrication theory and perturbation theory. The equation of liquid height and the process of thermal hydrodynamics of the non-Newtonian liquid film are obtained, and the case of weak effect of the rheological fluid is solved in detail.

Streamline topology and dilute particle dynamics in a Karman vortex street flow

Three types of streamline topology in a Karman vortex street flow are shown under the variation of spatial parameters. For the motion of dilute particles in the Karman vortex street flow, there exist a route of bifurcation to a chaotic orbit and more attractors in a bifurcation diagram for the proportion of particle density to fluid density. Along with the increase of spatial parameters in the flow field, the bifurcation process is suspended, as well as more and more attractors emerge. In the motion of dilute particles, a drag term and gravity term dominate and result in the bifurcation phenomenon.

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.