Numerical Sensitivity Analysis and Partial Similarity of Porous Media Flows

A numerical optimisation approach to identify dominant dimensionless variables in porous media flows by sensitivity analysis is proposed. We have validated the approach at first by examining a simple oil reservoir theoretically and numerically as well. A more complex water-flooding reservoir is examined based on sensitivity analysis of oil recovery to the similarity parameters, thus demonstrating the feasibility of the proposed approach to identify dominant similarity parameters for water-oil two-phase flows.

Assessment of bursting detection methods with LES database of turbulent channel flows

Validated by comparison with DNS, numerical database of turbulent channel flows is yielded by Large Eddy Simulation (LES). Three conventional techniques: uv quadrant 2, VITA and mu-level techniques for detecting turbulent bursts are applied to the identification of turbulent bursts. With a grouping parameter introduced by Bogard & Tiedemann (1986) or Luchik & Tiederman (1987), multiple ejections detected by these techniques which originate from a single burst can be grouped into a single-burst event. The results are compared with experimental results, showing that all techniques yield reasonable average burst period. However, uv quadrant 2 and mu-level are found to be superior to VITA in having large threshold-independent range.

Effect of Pipe Diameter on Flow Pattern Transition in Vertical Oil Water Flows

This paper investiges the effect of pipe diameter on flow pattern transition boundary in oil water vertical flows, and proposes a model to determine the maximum inner diameter (D_{infty s}) of a pipe in which the slug flow would not occur When pipe inner diameter D>D_{infty s}, only bubble flow exists, while D

Simulation of Vortex-dominant and Turbulent Flows in Natural Environment

In the present talk, the simulation of vortex dominant and turbulent flows are primarily addressed. To cope with complicated circumstances in environmental flows we illustrate the strategy of combining simplified physical model and suitable algorithm by a few examples.

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.

Monte Carlo Modeling of Micro Scale Gas Flows

An information preservation (IP) method has been used to simulate many micro scale gas flows. It may efficiently reduce the statistical scatter inherent in conventional particle approaches such as the direct simulation Monte Carlo (DSMC) method. This paper reviews applications of IP to some benchmark problems. Comparison of the IP results with those given by experiment, DSMC, and the linearized Boltzmann equation, as well as the Navier-Stokes equations with a slip boundary condition, and the lattice Boltzmann equation, shows that the IP method is applicable to micro scale gas flows over the entire flow regime from continuum to free molecular.

IP Simulation of Micro Gas Flows under 3-D Head Sliders

Gas film lubrication of a three-dimensional flat read-write head slider is calculated using the information preservation (IP) method and the direct simulation Monte Carlo (DSMC) method, respectively. The pressure distributions on the head slider surface at different velocities and flying heights obtained by the two methods are in excellent agreement. IP method is also employed to deal with head slider with three-dimensional complex configuration. The pressure distribution on the head slider surface and the net lifting force obtained by the IP method also agree well with those of DSMC method. Much less (of the order about 10(2) less) computational time (the sum of the time used to reach a steady stage and the time used in sampling process) is needed by the IP method than the DSMC method and such an advantage is more remarkable as the gas velocity decreases.

A Breakdown Criterion of Free Molecular Flows and an Optimum Analysis of EBPVD

Two important issues in electron beam physical vapor deposition (EBPVD) are addressed. The first issue is a validity condition of the classical cosine law widely used in the engineering context. This requires a breakdown criterion of the free molecular assumption on which the cosine law is established. Using the analytical solution of free molecular effusion flow, the number of collisions (N-c) for a particle moving from an evaporative source to a substrate is estimated that is proven inversely proportional to the local Knudsen number at the evaporation surface. N-c = 1 is adopted as a breakdown criterion of the free molecular assumption, and it is verified by experimental data and DSMC results. The second issue is how to realize the uniform distributions of thickness and component over a large-area thin film. Our analysis shows that at relatively low evaporation rates the goal is easy achieved through arranging the evaporative source positions properly and rotating the substrate.

Analytical calculations of Eulerian and Lagrangian time correlations in turbulent shear flows

The Taylor series expansion method is used to analytically calculate the Eulerian and Lagrangian time correlations in turbulent shear flows. The short-time behaviors of those correlation functions can be obtained from the series expansions. Especially, the propagation velocity and sweeping velocity in the elliptic model of space-time correlation are analytically calculated and further simplified using the sweeping hypothesis and straining hypothesis. These two characteristic velocities mainly determine the space-time correlations.

An Atomistic-Continuum Hybrid Simulation Of Fluid Flows Over Superhydrophobic Surfaces

Recent experiments have found that slip length could be as large as on the order of 1 mu m for fluid flows over superhydrophobic surfaces. Superhydrophobic surfaces can be achieved by patterning roughness on hydrophobic surfaces. In the present paper, an atomistic-continuum hybrid approach is developed to simulate the Couette flows over superhydrophobic surfaces, in which a molecular dynamics simulation is used in a small region near the superhydrophobic surface where the continuum assumption is not valid and the Navier-Stokes equations are used in a large region for bulk flows where the continuum assumption does hold. These two descriptions are coupled using the dynamic coupling model in the overlap region to ensure momentum continuity. The hybrid simulation predicts a superhydrophobic state with large slip lengths, which cannot be obtained by molecular dynamics simulation alone.

Space-time correlations of fluctuating velocities in turbulent shear flows

Space-time correlations or Eulerian two-point two-time correlations of fluctuating velocities are analytically and numerically investigated in turbulent shear flows. An elliptic model for the space-time correlations in the inertial range is developed from the similarity assumptions on the isocorrelation contours: they share a uniform preference direction and a constant aspect ratio. The similarity assumptions are justified using the Kolmogorov similarity hypotheses and verified using the direct numerical simulation DNS of turbulent channel flows. The model relates the space-time correlations to the space correlations via the convection and sweeping characteristic velocities. The analytical expressions for the convection and sweeping velocities are derived from the Navier-Stokes equations for homogeneous turbulent shear flows, where the convection velocity is represented by the mean velocity and the sweeping velocity is the sum of the random sweeping velocity and the shearinduced velocity. This suggests that unlike Taylor’s model where the convection velocity is dominating and Kraichnan and Tennekes’ model where the random sweeping velocity is dominating, the decorrelation time scales of the space-time correlations in turbulent shear flows are determined by the convection velocity, the random sweeping velocity, and the shear-induced velocity. This model predicts a universal form of the spacetime correlations with the two characteristic velocities. The DNS of turbulent channel flows supports the prediction: the correlation functions exhibit a fair good collapse, when plotted against the normalized space and time separations defined by the elliptic model.

An atomistic-continuum hybrid simulation of fluid flows over superhydrophobic surfaces

Recent experiments have found that slip length could be as large as on the order of 1 mu m for fluid flows over superhydrophobic surfaces. Superhydrophobic surfaces can be achieved by patterning roughness on hydrophobic surfaces. In the present paper an atomistic-continuum hybrid approach is developed to simulate the Couette flows over superhydrophobic surfaces in which a molecular dynamics simulation is used in a small region near the superhydrophobic surface where the continuum assumption is not valid and the Navier-Stokes equations are used in a large region for bulk flows where the continuum assumption does hold. These two descriptions are coupled using the dynamic coupling model in the overlap region to ensure momentum continuity. The hybrid simulation predicts a superhydrophobic state with large slip lengths which cannot be obtained by molecular dynamics simulation alone.

Kinetic Study Of The Rayleigh-Benard Flows

The present paper employs the direct simulation Monte Carlo (DSMC) method to study the Rayleigh-Benard flows, where the temperature ratio of the upper to lower plate is fixed to 0.1. For a Knudsen number (Kn) of 0.01, as the Rayleigh number (Ra) increases, the flow changes from the thermal conductive state to the convective state at about Ra=1700, and the calculated relation of heat flux through the lower plate versus Ra is in good agreement with classical experimental and theoretical results. For Kn=0.05, the thermal conductive state remains stable, and the increase of Ra cannot trigger thermal instability.

Kinetic studies of gas flows

Our recent studies on kinetic behaviors of gas flows are reviewed in this paper. These flows have a wide range of background, but share a common feature that the flow Knudsen number is larger than 0.01. Thus kinetic approaches such as the direct simulation Monte Carlo method are required for their description. In the past few years, we studied several micro/nano-scale flows by developing novel particle simulation approach, and investigated the flows in low-pressure chambers and at high altitude. In addition, the microscopic behaviors of a couple of classical flow problems were analyzed, which shows the potential for kinetic approaches to reveal the microscopic mechanism of gas flows.

A Simple Model For Predicting The Void Fraction Of Gas/Non-Newtonian Fluid Intermittent Flows In Upward Inclined Pipes

In this work, a simple correlation, which incorporates the mixture velocity, drift velocity, and the correction factor of Farooqi and Richardson, was proposed to predict the void fraction of gas/non-Newtonian intermittent flow in upward inclined pipes. The correlation was based on 352 data points covering a wide range of flow rates for different CMC solutions at diverse angles. A good agreement was obtained between the predicted and experimental results. These results substantiated the general validity of the model presented for gas/non-Newtonian two-phase intermittent flows.