生物芯片微通道周期性电渗流特性

以双电层的Poisson-Boltzmann方程和黏性不可压缩流体运动的Navier-Stokes方程为 基础,提出二维均匀微通道周期电渗流的解析解. 分析结果表明,周期电渗流速度大 小不但与双电层特性和外电场有关, 而且与流动雷诺数(Re = \omega h^2/\nu )密切相关. 随雷诺数增加,双电层滑移速度下降. 当离开固壁距离增加时,双电层以外区域流动速度快 速衰减,速度滞后相位角明显增加. 研究发现在微通道有波浪状速度剖面. 给出在低雷 诺数时的周期电渗流渐近解,它的速度振幅与定常电渗流速度相同,并具有柱栓式速度分布 形态. 还得到在微通道宽对双电层厚的比值(\kappa h)很小时,Debye-H\"{u}ckel近似 的周期电渗流解, 并与解析解进行分析比较 微通道,双电层,周期电渗流,雷诺数     

Electroosmotic flow in slit microchannel containing salt-free solution

A theoretical study on the electroosmotic flow through a uniformly charged planar slit microchannel containing a salt-free solution is presented. Based on the exact analytical solutions for the electric potential distribution and the transient electroosmotic flow velocity, a systematic parametric study on the characteristics of the transient electroosmotic flow is then investigated. The results show that the general behavior of electroosmotic flow in a planar slit is similar to that in a cylindrical capillary; however, the characteristic time to reach the steady-state flow, the electroosmotic mobility and the correction factor to the Smoluchowski equation in a slit are larger than those in a cylinder with its diameter equal to the slit width. Furthermore, the osmotic pressure of counterion between two equally charged planar surfaces immersed in a salt-free solution is found to be always repulsive. In addition, the implications of these results on electroosmosis in a microchannel are discussed.     

微/纳流控系统电渗流研究进展

简要介绍了微/纳流控系统中双电层和电渗流的基本原理、当前研究热点以及最新进展.从胶体界面理论出发,基于流体连续性模型和微纳米尺度流动的多物理场数值分析法,介绍了电渗流特性, 如焦耳热效应、 反离子效应、 表面电场调控双电层、双流体输运和周期电渗输运等电渗流控制方法.最后对电渗流测量作简要介绍.      

Numerical simulation of non-Darcian flow through a porous medium

This work reports on fluid flow in a fluid-saturated porous medium, accounting for the boundary and inertial effects in the momentum equation. The flow is simulated by Brinkman-Forchheimer-extended Darcy formulation (DFB), using MAC (Marker And Cell) and Chorin pressure iteration method. The method is validated by comparison with analytic results. The effect of Reynolds number, Darcy number, porosity and viscosity ratio on velocity is investigated. As a result, it is found that Darcy number has a decisive influence on pressure as well as velocity, and the effect of viscosity ratio on velocity is very strong given the Darcy number. Additional key findings include unreasonable choice of effective viscosity can involve loss of important physical information.     

The evaluation of the influence of various parameters on the velocity of elastic wave propagation in a rock medium

A full waveform recording in a borehole during acoustic logging makes it possible to determine the elastic parameters of a medium under in-situ conditions.The velocity of elastic wave propagation in rocks and elastic moduli are influenced by factors connected with its macrostructure and microstructure, as well as with rock overburden and porous pressure and temperature.The results of the calculations of the relationships between the elastic and reservoir parameters of sedimentary rocks are presented in this paper. The theoretical Kuster and Toksöz model has been applied.The influence of the porosity, the pore space coefficient, and the saturation of different media of porous rocks on elastic moduli and on compressional and shear wave propagation have been considered in this model. The complex composition of the skeleton and the influence of clay material in the porous rock are taken into account.     

Filtration Behaviour of Cement-Based Grout in Porous Media

Filtration behaviour of cement particles, especially under the high grouting pressure with a rapid grout flow velocity, has a significant effect on the grout injection. However, there have been few studies on this field where the governing equation of this behaviour remains unclear. In the present study, a novel experimental procedure for grout injection was adopted to acquire the spatial and temporal variations in porosity and viscosity of high-speed grout flow in coarse sand. Experimental observations showed that there were dramatic variations in viscosity and porosity during the grout penetration within the first 50 s, suggesting that the high velocity had a significant influence on the distribution of the filtration coefficient. A model based on the Stokes–Brinkman (S–B) equation and advection–filtration equations was established to describe the filtration of grout flow in porous media. Meanwhile, numerical solutions from both the proposed model and traditional Darcy’s law were compared with experimental results. The comparative results showed that the proposed approach can match the laboratory tests well; the analysis indicated that Darcy’s law was unable to properly describe high-speed grout flow in porous media due to the lack of a shear force and the inertial term. Nonuniform filtration behaviour of cement grout flowing in porous media was revealed. Due to the nonuniform distribution of the pore velocity isoline caused by Poiseuille flow, it led to a heterogenous distribution of porosity as well. Parametric studies on the applicability of Darcy’s law and S–B equation for grout flow were discussed, in which an error of less than 10% was calculated when the Reynolds number was less than 2.5.     

Transient Pressure Behavior of Reservoirs with Discrete Conductive Faults and Fractures

Fractures and faults are common features of many well-known reservoirs. They create traps, serve as conduits to oil and gas migration, and can behave as barriers or baffles to fluid flow. Naturally fractured reservoirs consist of fractures in igneous, metamorphic, sedimentary rocks (matrix), and formations. In most sedimentary formations both fractures and matrix contribute to flow and storage, but in igneous and metamorphic rocks only fractures contribute to flow and storage, and the matrix has almost zero permeability and porosity. In this study, we present a mesh-free semianalytical solution for pressure transient behavior in a 2D infinite reservoir containing a network of discrete and/or connected finite- and infinite-conductivity fractures. The proposed solution methodology is based on an analytical-element method and thus can be easily extended to incorporate other reservoir features such as sealing or leaky faults, domains with altered petrophysical properties (for example, fluid permeability or reservoir porosity), and complicated reservoir boundaries. It is shown that the pressure behavior of discretely fractured reservoirs is considerably different from the well-known Warren and Root dual-porosity reservoir model behavior. The pressure behavior of discretely fractured reservoirs shows many different flow regimes depending on fracture distribution, its intensity and conductivity. In some cases, they also exhibit a dual-porosity reservoir model behavior.     

多孔介质——自由流界面应力跳跃条件下流动特性解析解

对非对称多孔介质--自由流复合通道内多孔介质内部及多孔介质与自由流体界面处复杂质量、动量输运特性进行研究. 在多孔介质区采用Brinkman-extended Darcy模型并结合速度连续,剪切应力跳跃的界面条件对此复合通道内流体的传递现象进行求解,提出了考虑界面应力跳跃时非对称复合通道各区域流体运动速度及摩擦系数的解析式,分析了界面应力跳跃系数,达西数及无量纲多孔层偏心厚度对流体速度及摩擦系数的影响. 结果表明:改变界面性质可在一定条件下明显控制各区域流体速度分布;在达西数、多孔层偏心厚度一定情况下,界面应力系数的增大会使界面流速减小,而使流体摩擦系数增大,特别是界面应力系数小于0的情况下变化更明显,此时若不考虑界面应力系数则会造成较大误差. 当界面应力系数及多孔层偏心厚度均为较小负数值时,改变多孔层偏心厚度对界面速度的影响要大于改变界面应力系数的情况;而当界面应力系数及多孔层偏心厚度为较大正数值时,情况则相反. 较大达西数下,界面应力系数及多孔层偏心厚度对流体摩擦系数的影响均较大,继续减小达西数至一定程度时,界面应力系数对流体摩擦系数的影响可忽略不计而认为只与多孔层偏心厚度相关,且对较大多孔层偏心厚度更敏感.      

Thickness-averaged model for numerical simulation of electroosmotic flow in three-dimensional microfluidic chips

The microfluidic system is a multi-physics interaction field that has attracted great attention. The electric double layers and electroosmosis are important flow-electricity interaction phenomena. This paper presents a thickness-averaged model to solve three-dimensional complex electroosmotic flows in a wide-shallow microchannel/chamber combined (MCC) chip based on the Navier-Stokes equations for the flow field and the Poisson equation to the electric field. Behaviors of the electroosmotic flow, the electric field, and the pressure are analyzed. The quantitative effects of the wall charge density (or the zeta potential) and the applied electric field on the electroosmotic flow rate are investigated. The two-dimensional thickness-averaged flow model greatly simplifies the three-dimensional computation of the complex electroosmotic flows, and correctly reflects the electrookinetic effects of the wall charge on the flow. The numerical results indicate that the electroosmotic flow rate of the thickness-averaged model agrees well with that of the three-dimensional slip-boundary flow model. The flow streamlines and pressure distribution of these two models are in qualitative agreement.     

Simulation of Pressure Transient Behavior for Asymmetrically Finite-Conductivity Fractured Wells in Coal Reservoirs

Based on Fick’s law in matrix and Darcy flow in cleats and hydraulic fractures, a new semi-analytical model considering the effects of boundary conditions was presented to investigate pressure transient behavior for asymmetrically fractured wells in coal reservoirs. The new model is more accurate than previous model proposed by Anbarci and Ertekin, SPE annual technical conference and exhibition, New Orleans, 27–30 Sept 1998 because new model is expressed in the form of integral expressions and is validated well through numerical simulation. (1) In this paper, the effects of parameters including fracture conductivity, coal reservoir porosity and permeability, fracture asymmetry factor, sorption time constant, fracture half-length, and coalbed methane (CBM) viscosity on bottomhole pressure behavior were discussed in detail. (2) Type curves were established to analyze both transient pressure behavior and flow characteristics in CBM reservoir. According to the characteristics of dimensionless pseudo pressure derivative curves, the process of the flow for fractured CBM wells was divided into six sub-stages. (3) This paper showed the comparison of transient steady state and pseudo steady state models. (4) The effects of parameters including transfer coefficient, wellbore storage coefficient, storage coefficient of cleat, fracture conductivity, fracture asymmetry factor, and rate coefficient on the shape of type curves were also discussed in detail, indicating that it is necessary to keep a bigger fracture conductivity and fracture symmetry for enhancing well production and reducing pressure depletion during the hydraulic fracturing design.     

Gas-gas displacement in a fractured porous reservoir

The results of an experimental investigation of the process of miscible displacement of a gas from a fractured porous reservoir by another gas are presented. It is established that the displacement process is influenced by convective and molecular diffusion. The displacement coefficient depends both on the reservoir properties of the formation and on the technological displacement parameters. The experiments revealed a decrease in the displacement coefficient with increase in the permeability ratio of the joints and the blocks and with decrease in formation pressure. The dependence of the displacement coefficient on the average gas velocity in the reservoir, i.e., the ratio of the gas flow rate to the product of the flow area and the total reservoir porosity coefficient (including the blocks and the joints), is established.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 74–79, May–June, 1991.     

煤系地层岩石渗透特性试验研究

基于MTS815.02型岩石伺服渗透试验系统,应用两种实验方法测定了岩石的渗透特性,给出煤系地层中多种岩石全应力-应变过程中渗透率变化范围,对几种岩石的渗透特征进行了分析和讨论。研究表明:1)煤系地层岩石全应力-应变过程的渗透率通常在10-9到10-4Darcy量级之间;2)不同岩石在应力峰值前后其渗透性呈现不同的变化趋势;3)岩石渗透率的离散性很大,即使在相同的试验条件下,渗透率变异系数大多在0.65以上,其离散性按粗砂岩、泥岩、砂质页岩、中砂岩、细砂岩、砾岩、石灰岩顺序逐渐增大;4)在MTS815.02系统上,当所测岩样渗透率小于10-4Darcy量级时,应采用瞬态法,当渗透率量级大于10-4Darcy时,可改用稳态法测定。岩石变形过程中应力场和渗流场的耦合作用十分复杂,渗透率受孔隙压力、围压、试件尺寸、饱和程度等因素影响。     

Thermal transport and performance analysis of pressure- and electroosmotically-driven liquid flow microchannel heat sink with wavy wall

The present study investigates the microchannel heat sinks (MCHSs) with smooth and wavy wall for pure electroosmotic flow (EOF), pressure-driven flow (PDF) and combined electroosmotic and pressure-driven flow (PDF + EOF). A three-dimensional numerical analysis was performed for EOF, PDF and combined flow (PDF + EOF) through finite volume analysis. The EOF was combined with the PDF to enhance the flow rate and to reduce the thermal resistance of the MCHS. The effect of wall waviness on electroosmosis and thermal performance of the MCHS was critically investigated for flow rate, friction factor, Nusselt number, thermal resistance and pumping power. The design variables related to the wavelength and amplitude and width of microchannel were investigated for their effect on the overall thermal performance and pumping power. The electroosmosis not only increases the flow rate but also suppresses the secondary flow developed due to the topology of the microchannel walls. The non-uniformity of the velocity and temperature is reduced due to the application of the EOF in a PDF and combined flow (PDF + EOF).     

页岩气在超低渗介质中的渗流行为

基于页岩气藏复杂孔隙结构和页岩气在纳米孔隙表面的滑脱和吸附-解吸附等现象,通过引入表观渗透率,修正Darcy渗流模型,得到了页岩气渗流本构方程. 将计算结果与Darcy模型计算结果进行了比照,结果表明,在产量定常情形下,基于非Darcy渗流模型得到的井底压力高于Darcy流结果; Darcy流模型得到的压力衰减速度较快,改进后的渗流模型更能准确描述与表征页岩气渗流过程;研究结果可为页岩气藏的经营与管理提供基础参数.      

Analytical study of mixed electroosmotic-pressure-driven flow in rectangular micro-channels

Operational state of many miniaturized devices deals with flow field in microchannels. Pressure-driven flow (PDF) and electroosmotic flow (EOF) can be recognized as the two most important types of the flow field in such channels. EOF has many advantages in comparison with PDF, such as being vibration free and not requiring any external mechanical pumps or moving parts. However, the disadvantages of this type of flow such as Joule heating, electrophoresis demixing, and not being suitable for mobile devices must be taken into consideration carefully. By using mixed electroosmotic/pressure-driven flow, the role of EOF in producing desired velocity profile will be reduced. In this way, the advantages of EOF can be exploited, and its disadvantages can be prevented. Induced pressure gradient can be utilized in order to control the separation in the system. Furthermore, in many complicated geometries such as T-shape microchannels, turns may induce pressure gradient to the electroosmotic velocity. While analytical formulas are completely essential for analysis and control of any industrial and laboratory microdevices, lack of such formulas in the literature for solving Poisson–Boltzmann equation and predicting electroosmotic velocity field in rectangular domains is evident. In the present study, first a novel method is proposed to solve Poisson–Boltzmann equation (PBE). Subsequently, this solution is utilized to find the electroosmotic and the mixed electroosmotic/pressure-driven velocity profile in a rectangular domain of the microchannels. To demonstrate the accuracy of the presented analytical method in solving PBE and finding electroosmotic velocity, a general nondimensional example is analyzed, and the results are compared with the solution of boundary element method. Additionally, the effects of different nondimensional parameters and also aspect ratio of channels on the electroosmotic part of the flow field will be investigated.     

基于REV尺度格子Boltzmann方法的页岩气流动数值模拟

结合页岩扫描电镜图像,提出页岩气藏物理模型,采用表征单元体积(representative elementary volume,REV)尺度格子Boltzmann方法,考虑滑脱效应,模拟页岩气在页岩气藏中的流动.模拟结果表明,页岩气主要沿着天然裂缝窜进,但在有机质和无机质中也存在缓慢的流动,且有机质中的流速要略大于无机质中的流速.通过改变地层压力,研究地层压力对页岩气渗流特性的影响.研究结果表明,整个流场的速度和渗透率均随着地层压力的下降而增加.     

Electro-magneto-hydrodynamic flow of couple stress nanofluids in micro-peristaltic channel with slip and convective conditions

This study explores the effects of electro-magneto-hydrodynamics,Hall currents,and convective and slip boundary conditions on the peristaltic propulsion of nanofluids(considered as couple stress nanofluids)through porous symmetric microchannels.The phenomena of energy and mass transfer are considered under thermal radiation and heat source/sink.The governing equations are modeled and non-dimensionalized under appropriate dimensionless quantities.The resulting system is solved numerically with MATHEMATICA(with an in-built function,namely the Runge-Kutta scheme).Graphical results are presented for various fluid flow quantities,such as the velocity,the nanoparticle temperature,the nanoparticle concentration,the skin friction,the nanoparticle heat transfer coefficient,the nanoparticle concentration coefficient,and the trapping phenomena.The results indicate that the nanoparticle heat transfer coefficient is enhanced for the larger values of thermophoresis parameters.Furthermore,an intriguing phenomenon is observed in trapping:the trapped bolus is expanded with an increase in the Hartmann number.However,the bolus size decreases with the increasing values of both the Darcy number and the electroosmotic parameter.     

Flow instability in a curved porous channel formed by two concentric cylindrical surfaces

Stability of laminar flow in a curved channel formed by two concentric cylindrical surfaces is investigated. The channel is occupied by a fluid saturated porous medium; the flow in the channel is driven by a constant azimuthal pressure gradient. The momentum equation takes into account two drag terms: the Darcy term that describes friction between the fluid and the porous matrix, and the Brinkman term, which allows imposing the no-slip boundary condition at the channel walls. An analytical solution for the basic flow velocity is obtained. Numerical analysis is carried out using the collocation method to investigate the onset of instability leading to the development of a secondary motion in the form of toroidal vortices. The dependence of the critical Dean number on porosity and the channel width is analyzed.     

Periodical streaming potential and electro-viscous effects in microchannel flow

This paper presents an analytical solution to periodical streaming potential, flow-induced electric field and velocity of periodical pressure-driven flows in twodimensional uniform microchannel based on the Poisson-Boltzmann equations for electric double layer and Navier-Stokes equation for liquid flow. Dimensional analysis indicates that electric-viscous force depends on three factors： （1） Electric-viscous number representing a ratio between maximum of electric-viscous force and pressure gradient in a steady state, （2） profile function describing the distribution profile of electro-viscous force in channel section, and （3） coupling coefficient reflecting behavior of arnplitude damping and phase offset of electro-viscous force. Analytical results indicate that flow-induced electric field and flow velocity depend on frequency Reynolds number （Re = wh^2/v）. Flow-induced electric field varies very slowly with Re when Re 〈 1, and rapidly decreases when Re 〉 1. Electro-viscous effect on flow-induced electric field and flow velocity are very significant when the rate of the channel width to the thickness of electric double layer is small.     

Periodical streaming potential and electro-viscous effects in microchannel flow

This paper presents an analytical solution to periodical streaming potential,flow-induced electric field and velocity of periodical pressure-driven flows in twodimensional uniform microchannel based on the Poisson.Boltzmann equations for electric double layer and Navier-Stokes equation for liquid flow.Dimensional analysis indicates that electric-Viscous force depends on three factors:(1)Electric-viscous number representing a ratio between maximum of electric-viscous force and pressure gradient in a steady state,(2)profile function describing the distribution profile of electro-viscous forcein channel section,and(3)coupling coefficient reflecting behavior of amplitude damping and phase Offset of electro-viscous force.Analytical results indicate that flow-induced electric field and flow velocity depend on frequency Reynolds number(Re=wh2/v).Flow-induced electric field varies very slowly with Re when Re＜1.and rapidly decreases when Re＞1.Electro-viscous effect on flow-induced electric field and flow velocity are very significant when the rate of the channel width to the thickness of electric double layer is small.