Boundary conditions at the solid-liquid surface over the multiscale channel size from nanometer to micron

The boundary condition at the solid surface is one of the important problems for the microfluidics. In this paper we study the effects of the channel sizes on the boundary conditions (BC), using the hybrid computation scheme adjoining the molecular dynamics (MD) simulations and the continuum fluid mechanics. We could reproduce the three types of boundary conditions (slip, no-slip and locking) over the multiscale channel sizes. The slip lengths are found to be mainly dependent on the interfacial parameters with the fixed apparent shear rate. The channel size has little effects on the slip lengths if the size is above a critical value within a couple of tens of molecular diameters. We explore the liquid particle distributions nearest the solid walls and found that the slip boundary condition always corresponds to the uniform liquid particle distributions parallel to the solid walls, while the no-slip or locking boundary conditions correspond to the ordered liquid structures close to the solid walls. The slip, no-slip and locking interfacial parameters yield the positive, zero and negative slip lengths respectively. The three types of boundary conditions existing in "microscale" still occur in "macroscale". However, the slip lengths weakly dependent on the channel sizes yield the real shear rates and the slip velocity relative to the solid wall traveling speed approaching those with the no-slip boundary condition when the channel size is larger than thousands of liquid molecular diameters for all of the three types of interfacial parameters, leading to the quasi-no-slip boundary conditions.

Study on methane hydration process in a semi-continuous stirred tank reactor

The methane hydration process is investigated in a semi-continuous stirred tank reactor. Liquid temperatures and reaction rates without stirrer are compared with those occurring with stirrer, while at the same time better stirring conditions of the methane hydration process are given by the experiments. Some basic data of fluid mechanics, for example, stirring Reynolds number, Froucle number and stirrer power, are calculated during the methane hydration process, which can be applied to evaluate stirrer capacity and provide some basic data for a scaled up reactor. Based on experiment and calculations in this work, some conclusions are drawn. First, the stirrer has great influence on the methane hydration process. Batch stirring is helpful to improve the mass transfer and heat transfer performances of the methane hydration process. Second, induction time can be shortened effectively by use of the stirrer. Third, in this paper, the appropriate stirring velocity and stirring time were 320 rpm and 30 min, respectively, at 5.0 MPa, for which the storage capacity and reaction time were 159.1 V/V and 370 min, respectively. Under the condition of the on-flow state, the initial stirring Reynolds number of the fluid and the stirring power were 12,150 and 0.54 W, respectively. Fourth, some suggestions, for example, the use of another type of stirrer or some baffles, are proposed to accelerate the methane hydration process. Comparing with literature data, higher storage capacity and hydration rate are achieved in this work. Moreover, some fluid mechanics parameters are calculated, which can provide some references to engineering application.

A New Criterion Number for the Boundary Conditions at the Solid/Liquid Interface in Nanoscale

A simple, but important three-atom model was proposed at the solid/liquid interface, leading to a new criterion number, lambda, governing the boundary conditions (BCs) in nanoscale. The solid wall is considered as the face-centered-cubic (fcc) structure. The fluid is the liquid argon with the well-known LJ potential. Based on the concept, the two micro-systems have the same BCs if they have The same criterion number. The degree of the locking BCs is enhanced when lambda equals to 0.757. Such critical criterion number results in the substantial epitaxial ordering and one, two, or even three liquid layers are locked by the solid wall, depending on the coupling energy scale ratio of the solid and liquid atoms. With deviation from the critical criterion number, the flow approaches the slip BCs and there are little ordering structures within the liquid. Always at the same criterion number, the degree of the slip is decreased or the locking is enhanced with increasing the coupling energy scale ratio of the solid and liquid atoms. The above analysis is well confirmed by the molecular dynamics (MD) simulation. The slip length is well correlated in terms of the new criterion number. The future work is suggested to extend the present theory for other microstructures of the solid wall atoms and quasi-LJ potentials.

Space grown semi-insulating gallium arsenide single crystal and its application

Low noise field effect transistors and analogue switch integrated circuits (ICs) have been fabricated in semi-insulating gallium arsenide (SI-GaAs) wafers grown in space by direct ion-implantation. The electrical behaviors of the devices and the ICs have surpassed those fabricated in the terrestrially grown SI-GaAs wafers. The highest gain and the lowest noise of the transistors made from space-grown SI-GaAs wafers are 22.8 dB and 0.78 dB, respectively. The threshold back-gating voltage of the ICs made from space-grown SI-GaAs wafers is better than 8.5 V The con-elation between the characterizations of materials and devices is studied systematically. (C) 2002 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

Numerical simulation of a bubble rising in shear-thinning fluids

The motion of a single bubble rising freely in quiescent non-Newtonian viscous fluids was investigated experimentally and computationally. The non-Newtonian effects in the flow of viscous inelastic fluids are modeled by the Carreau theological model. An improved level set approach for computing the incompressible two-phase flow with deformable free interface is used. The control volume formulation with the SIMPLEC algorithm incorporated is used to solve the governing equations on a staggered Eulerian grid. The simulation results demonstrate that the algorithm is robust for shear-thinning liquids with large density (rho(1)/rho(g) up to 10(3)) and high viscosity (eta(1)/eta(g) up to 10(4)). The comparison of the experimental measurements of terminal bubble shape and velocity with the computational results is satisfactory. It is shown that the local change in viscosity around a bubble greatly depends on the bubble shape and the zero-shear viscosity of non-Newtonian shear-thinning liquids. The shear-rate distribution and velocity fields are used to elucidate the formation of a region of large viscosity at the rear of a bubble as a result of the rather stagnant flow behind the bubble. The numerical results provide the basis for further investigations, such as the numerical simulation of viscoelastic fluids. (C) 2010 Elsevier B.V. All rights reserved.

Linear spatio-temporal instability analysis of ice growth under a falling water film

A linear spatio-temporal stability analysis is conducted for the ice growth under a falling water film along an inclined ice plane. The full system of linear stability equations is solved by using the Chebyshev collocation method. By plotting the boundary curve between the linear absolute and convective instabilities (AI/CI) of the ice mode in the parameter plane of the Reynolds number and incline angle, it is found that the linear absolute instability exists and occurs above a minimum Reynolds number and below a maximum inclined angle. Furthermore, by plotting the critical Reynolds number curves with respect to the inclined angle for the downstream and upstream branches, the convectively unstable region is determined and divided into three parts, one of which has both downstream and upstream convectively unstable wavepackets and the other two have only downstream or upstream convectively unstable wavepacket. Finally, the effect of the Stefan number and the thickness of the ice layer on the AI/CI boundary curve is investigated.

Hydroelastic dynamic characteristics of a slender axis-symmetric body

The slender axis-symmetric submarine body moving in the vertical plane is the object of our investigation. A coupling model is developed where displacements of a solid body as a Euler beam (consisting of rigid motions and elastic deformations) and fluid pressures are employed as basic independent variables, including the interaction between hydrodynamic forces and structure dynamic forces. Firstly the hydrodynamic forces, depending on and conversely influencing body motions, are taken into account as the governing equations. The expressions of fluid pressure are derived based on the potential theory. The characteristics of fluid pressure, including its components, distribution and effect on structure dynamics, are analyzed. Then the coupling model is solved numerically by means of a finite element method (FEM). This avoids the complicacy, combining CFD (fluid) and FEM (structure), of direct numerical simulation, and allows the body with a non-strict ideal shape so as to be more suitable for practical engineering. An illustrative example is given in which the hydroelastic dynamic characteristics, natural frequencies and modes of a submarine body are analyzed and compared with experimental results. Satisfactory agreement is observed and the model presented in this paper is shown to be valid.

Sedimentary processes and development of the Zenisu deep-sea channel, Philippine Sea

Zenisu deep-sea channel originated from a volcanic arc region, Izu-Ogasawara Island Arc, and vanished in the Shikoku Basin of the Philippine Sea. According to the swath bathymetry, the deep-sea channel can be divided into three,segments. They are Zenisu canyon, E-W fan channel and trough-axis channel. A lot of volcanic detritus were deposited in the Zenisu Trough via the deep-sea channel because it originated from volcanic arc settings. On the basis of the swath bathymetry, submersible and seismic reflection data, the deposits are characterized by turbidite and debrite deposits as those in the other major deep-sea channels. Erosion or few sediments were observed in the Zenisu canyon, whereas a lot of turbidites and debrites occurred in the E-W channel and trough axis channel. Cold seep communities, active fault and fluid flow were discovered along the lower slope of the Zenisu Ridge. Vertical sedimentary sequences in the Zenisu Trough consist of the four post-rift sequence units of the Shikoku Basin, among which Units A and B are two turbidite units. The development of Zenisu canyon is controlled by the N-S shear fault, the E-W fan channel is related to the E-W shear fault, and the trough-axis channel is related to the subsidence of central basin.

Dynamics of the buoyant plume off the Pearl River Estuary in summer

Field measurements of salinity, wind and river discharge and numerical simulations of hydrodynamics from 1978 to 1984 are used to investigate the dynamics of the buoyant plume off the Pearl River Estuary (PRE), China during summer. The studies have shown that there are four major horizontal buoyant plume types in summer: Offshore Bulge Spreading (Type I), West Alongshore Spreading (Type II), East Offshore Spreading (Type III), and Symmetrical Alongshore Spreading (Type IV). River mouth conditions, winds and ambient coastal currents have inter-influences to the transport processes of the buoyant plume. It is found that all of the four types are surface-advected plumes by analysing the vertical characteristic of the plumes, and the monthly variations of the river discharge affect the plume size dominantly. The correlation coefficient between the PRE plume size and the river discharge reaches 0.85 during the high river discharge season. A wind strength index has been introduced to examine the wind effect. It is confirmed that winds play a significant role in forming the plume morphology. The alongshore wind stress and the coastal currents determine the alongshore plume spreading. The impact of the ambient currents such as Dongsha Current and South China Sea (SCS) Warm Current on the plume off the shelf has also assessed. The present study has demonstrated that both the river discharge and wind conditions affect the plume evolution.

东营凹陷岩性油气藏成藏动力学机制和预测

This paper is belonging to Chinese Petrochemical Industry Corporation's key project. Although it is very difficult, it has important theoretical and practical value. Its targets is to make lithological petroleum pool exploration great breakthrough in Dongying sag, by applying advanced theories, the last-minute methods and technology in highly explored zones. By using synthetically multi- discipline theories, methods and technology such as petroleum geology, sedimentology, structure geology, rock mechanics, dynamics of petroleum pool formation, geochemistry, geophysics and so on, and by making full use of computer , the process of petroleum pool forming and distribution rules of lithological petroleum pools have been thoroughly investigated and analyzed in sharp-slope, gentle-slope as well as low-lying region of Dongying sag including dynamic and static. With the study of tectonic stress field, fluid potential field and pressure field, we revealed dynamics condition, distribution rule, control factors and petroleum forming mechanism of lithological pool, and established the forming mode of lithological pool of Dongying sag. The main conclusion as follow: Strata framework, structure framework and sedimentary system of Dongying sag have been established which were the basis of petroleum prediction. There are three kinds of oil source which were from Es4,Es3 and mixed type, also three petroleum forming phases which were the telophase of Dongying stage, Guantao stage and Minghuazhen group, which occur in different geological environment. By using of most advanced numerical modeling software, the space distribution and time evolve of stress field and fluid potential field have been revealed from Esl up to the present. The region with low earth stress and low fluid potential were enrichment region of lithological petroleum pool and fault-block pool. The dynamics mechanism of Lithological petroleum pool in Dongying sag was collocating seal box, abnormity pressure, index number of petroleum forming and static factors on time and space, which was the most important factor of controlling petroleum pool forming, distribution and enrichment. The multi phase active and evolve of seal and unseal about different order fault were main factors of controlling petroleum pool forming of Dongying sag, which have important value for predicting lithological petroleum pool. It is revealed the lithological petroleum pool forming mode that included respective character, forming mechanism and distribution rule in four structural belt, which was a base for lithological petroleum pool prediction. The theories, technology and methods of studying, description, characterize and prediction lithological petroleum pool were established, which have important popularization value. Several lithological pool have been predicted in stress transform, zone, abrupt slope zone, fractured surface changed zone, tosional stress growth zone and abnormity pressure zone with noticeable economic benefit after exploration.

胜坨油田沙二段油藏流场剩余油形成机理和分布模式

By applying multi-discipline theory and methods comprehensively and with full use of computer, the paper deeps into studying changing rule and control factor of fluid field of ES2 Shengtuo oil field during waterflood development, physical and chemical function, and stress. Matrix field, network field, fluid field, stress field and physical chemistry field and fluid model for dynamic function were established. Macroscopic and microscopic genesis mechanics, distribution rule and control factor of remaining oil were revealed. Remaining oil and emulate model were established. Macroscopic and microscopic distribution rule of mover remaining oil were predicted, several results were achieved as following: The distribution of remaining oil was controlled by micro-structure. At the same development stage, remaining oil saturation of the wells located in higher position of micro-structure is higher than the average saturation in the same layer. The water content ratio has same law. It is the enrichment district that the high position of micro-structure controlled by seal faults. The remaining oil distribution was affected by sedimentary micro-facies, micro-structure, fault sealing, reservoir heterogeneity and affusion-oil extraction. On the plane, the zone owning higher saturation of remaining oil is the area that at the edge miacro-facies and sand-body distribution discontinuously; on the section, the content of waterflood of the upper or middle-upper oil layer of positive rhythm and positive comprehensive rhythm is lower, middle and weak waterflood is main, remaining oil is in enrichment relatively. The remaining oil is relative enrichment at the zones of well network of affusion and oil extraction not affected. 4D dynamic model of reservoir of Es2in Shengtuo oil field was established. Macroscopic and microscopic forming mechanics, distribution rule and control factor were revealed. The emulate model of dynamic function of Shengtuo oil field was established, the space distribution of remaining oil were predicted. Reservoir flow field, matrix field, network field, seep field, physical and chemical field, stress field and fluid field models were established. Reservoir flow field character and distribution were revealed. An improvement of the development geology theory in continental fault depression continental basin was brought on.

文33块沙二下段油藏动态模型及剩余油预测

By applying synthetically multi-subject theories, methods and technology, such as petroleum geology, sedimentology, seep mechanics, geochemistry, geophysics and so on; and by making full use of computer; combining quantity and quality, macroscopic and microscopic, intensive static and active description, comprehensive studying and physical modeling, 3 dimension and 4 dimension description; the paper took Wen-33 block of Zhongyuan oil field as an example; and studied reservoir macroscopic and microscopic parameter changing rule and evolve mechanics in different water-blood stage. The reservoir dynamic model and remaining-oil distribution mode was established, and several results were achieved as follows: (1) Three types of parameter gaining, optimizing and whole data body of Wen33th reservoir were established. Strata framework, structure framework, reservoir types and distribution of Wen33th reservoir were discussed. Reservoir genesis types, space distribution law and evolve rule of Wen33th reservoir were explained. 4D dynamic model of macroscopic parameter of reservoir flow dynamic geologic function of Wen33th reservoir was established. The macroscopic remaining-oil distribution and control factor was revealed. The models of the microscopic matrix field, pore-throat network field, fluid field, clay mineral field of Wen-33 block were established. The characters, changing rules and controlled factors in different water stage were revealed. The evolve rule and mechanics of petroleum fluid field in Wen-33 block reservoir were revealed. Macroscopic and microscopic remaining oil distribution mode of Wen-33 block were established. Seven types, namely 12 shapes of dynamic model of microscopic remaining oil were discussed, and the distribution of mover remaining oil was predicted. Emulation model: mathematical model and prediction model of Wen-33 block were established. The changing mechanics of reservoir parameter and distribution of remaining-oil were predicted. Firstly, the paper putting forward that the dynamic geologic function of petroleum development is the factor of controlling remaining-oil, which is the main factor leading to matrix field, network field, clay mineral field, fluid field, physic and chemical field, stress field and fluid field forming and evolving. (10) A set of theories, methods and technologies of investigating, describing, characterizing and predicting complex fault-block petroleum were developed.