Instability Of Plane Poiseuille Flow In A Fluid-Porous System

The instability of Poiseuille flow in a fluid-porous system is investigated. The system consists of a fluid layer overlying porous media and is subjected to a horizontal plane Poiseuille flow. We use Brinkman's model instead of Darcy's law to describe the porous layer. The eigenvalue problem is solved by means of a Chebyshev collocation method. We study the influence of the depth ratio (d) over cap and the Darcy number delta on the instability of the system. We compare systematically the instability of Brinkman's model with the results of Darcy's model. Our results show that no satisfactory agreement between Brinkman's model and Darcy's model is obtained for the instability of a fluid-porous system. We also examine the instability of Darcy's model. A particular comparison with early work is made. We find that a multivalued region may present in the (k, Re) plane, which was neglected in previous work. Here k is the dimensionless wavenumber and Re is the Reynolds number. (C) 2008 American Institute of Physics. [DOI: 10.1063/1.3000643]

Non-modal instabilities of two-dimensional disturbances in plane Couette flow of a power-law fluid

Instabilities of fluid flows have traditionally been investigated by normal mode analysis, i.e. by linearizing the equations of flow and testing for unstable eigenvalues of the linearized problem. However, the results of eigenvalue analysis agree poorly in many cases with experiments, especially for shear flows. In this paper we study the instabilities of two-dimensional Couette flow of a polymeric fluid in the framework of non-modal stability theory rather than normal mode analysis. A power-law model is used to describe the polymeric liquid. We focus on the response to external excitations and initial conditions by examining the pseudospectra structures and the transient energy growths. For both Newtonian and non-Newtonian flows, the results show that there can be a rather large transient growth even though the linear operator of Couette flow has no unstable eigenvalue. The effects of non-Newtonian viscosity on the transient behaviors are examined in this study. The results show that the "shear-thinning/shear-thickening" effect increases/decreases the amplitude of responses to external excitations and initial conditions. (C) 2010 Elsevier B.V. All rights reserved.

Generalized rayleigh principle and its applications

In this paper, we mainly deal with cigenvalue problems of non-self-adjoint operator. To begin with, the generalized Rayleigh variational principle, the idea of which was due to Morse and Feshbach, is examined in detail and proved more strictly in mathematics. Then, other three equivalent formulations of it are presented. While applying them to approximate calculation we find the condition under which the above variational method can be identified as the same with Galerkin's one. After that we illustrate the generalized variational principle by considering the hydrodynamic stability of plane Poiseuille flow and Bénard convection. Finally, the Rayleigh quotient method is extended to the cases of non-self-adjoint matrix in order to determine its strong eigenvalne in linear algebra.

Instabilities of a liquid film flowing down an inclined porous plane

The problem of a film flowing down an inclined porous layer is considered. The fully developed basic flow is driven by gravitation. A careful linear instability analysis is carried out. We use Darcy's law to describe the porous layer and solve the coupling equations of the fluid and the porous medium rather than the decoupled equations of the one-sided model used in previous works. The eigenvalue problem is solved by means of a Chebyshev collocation method. We compare the instability of the two-sided model with the results of the one-sided model. The result reveals a porous mode instability which is completely neglected in previous works. For a falling film on an inclined porous plane there are three instability modes, i.e., the surface mode, the shear mode, and the porous mode. We also study the influences of the depth ratio d, the Darcy number delta, and the Beavers-Joseph coefficient alpha(BJ) on the instability of the system.

光滑粒子动力学方法中粒子分布与数值稳定性分析

光滑粒子动力学(SPH)作为一种拉格朗日型无网格粒子方法,已经成功地应用于包括含多相流动界面以及移动边界的可压缩和不可压缩流体运动的研究中.通过对Poiseuille流动的深入研究,探索了SPH方法中粒子分布对计算精度的影响,揭示了一种因为粒子不规则分布而导致的数值不稳定现象.研究显示,这种数值不稳定性起源于SPH方法粒子近似过程中的不连续性.使用了一种新的粒子近似格式以确保SPH方法中粒子近似的连续性.计算结果表明,这种新的粒子近似格式对于规则和不规则的粒子分布都能得到稳定精度的结果.

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.

耗散粒子动力学处理复杂固体壁面的一种有效方法

耗散粒子动力学(dissipative particle dynamics,DPD)作为一种介观尺度拉格朗日型粒子方法,已经成功地应用于微纳米流动和生化科技的研究中. 复杂固体壁面的处理和壁面边界条件的实施一直是DPD方法发展及应用的一个障碍. 提出了处理复杂固体壁面的一种新的方法. 复杂固体区域通过冻结随机分布并且达到平衡状态的DPD粒子代表;所冻结的DPD粒子位于临近流动区域的一个截距内;在靠近固体壁面的流动区域中设置流动反弹层,当流动DPD粒子进入此流动层后反弹回流动区域. 应用这种固体壁面处理方法.

Flow Characteristics of Liquid with Pressure-Dependent Viscosities in Microtubes

It is obvious that the pressure gradient alone, the axial direction in a pipe flow keeps constant according to the Haoen-Poiseuille equation. However, recent experiments indicated that the distribution of the pressure seemed no longer linear for liquid flows in microtubes driven by high pressure (1-30MPa). Based on H-P equation with slip boundary condition and Bridgman's relation of viscosity vs. static pressure, the nonlinear distribution of pressure along the axial direction is analyzed in this paper. The revised standard Poiseuille number with the effect of pressure-dependent viscosity taken into account agrees well with the experimental results. Therefore, the dependence of the viscosity on the pressure is one of the dominating, factors under high driven pressure, and is represented by an important property coefficient et of the liquid.

Flow Characteristics of Non-Polar Organic Liquids with Small Molecules in a Microchannel

用去离子水及有机液体在内径约为25μm的石英圆管内进行了流量特性实验.液体分子量范围为18～160,动力黏性系数的范围为0.5～1 mPa.s.实验雷诺数范围为Re<8.所用有机液体为:四氯化碳、乙基苯及环己烷都是非极性液体,其分子结构尺度小于1 nm.实验结果表明,在定常层流条件下,圆管内的液体流量与两端压力差成正比,其压力-流量关系仍符合经典的Hagen-Poiseuille流动.这说明非极性小分子有机液体在本实验所用微米尺度管道中其流动规律仍符合连续介质假设.鉴于微尺度流动实验的特殊性,文中还介绍了微流动实验装置,分析了微尺度流动测量误差来源及提高测量精度的措施.

Flow Characteristics of Liquids in Microtubes Driven by a High Pressure

The flow characteristics of liquids in microtubes driven by a high pressure ranging from 1 MPa to 30 MPa are studied in this paper. The diameter of the microtube is from 3 μm to 10 μm and liquids composed of simple small molecules are chosen as the working fluids. The Reynolds number ranges from 0. 1 to 24. The behavior of isopropanol and carbon tetrachloride under high pressure is found different from the prediction from conventional Hagen-Poiseuille (HP) equation. The normalized friction coefficient C* increases significantly with the pressure. From an analysis of the microtube deformation, liquid compressibility, viscous heating and wall slip, it may be seen that the viscosity at high pressure plays an important role here. An exponential function of viscosity vs pressure is introduced into the HP equation to counteract the difference between experimental and theoretical values. However, this difference is not so marked for di-water.

Pressure Drop of Bubbly Two-Phase Flow in a Square Channel at Reduced Gravity

A new set of experimental pressure drop data, collected aboard the Russian IL-76MDK, is reported for bubbly airwater two-phase flow in a square channel with a cross-sectional area of 12x 12mm(2). The present data are compared to several frequently used empirical models, e.g. homogeneous model, Lockhart-Martinelli-Chisholm correlation and Friedel's model. It is shown that the predictions of the models mentioned above are generally not satisfied. A new homogeneous model is developed based on the analysis of the characteristics of bubbly two-phase flow at reduced gravity. It is tested with the present data and other data collected by other researchers in circular pipes. Some questions related to the present model are also discussed. (C) 2002 COSPAR. Published by Elsevier Science Ltd. All rights reserved.

Electro-Osmotic Flow and Mixing in Heterogeneous Microchannels

Analytical and numerical studies of secondary electro-osmotic flow EOF and its mixing in microchannels with heterogeneous zeta potentials are carried out in the present work. The secondary EOFs are analyzed by solving the Stokes equation with heterogeneous slip velocity boundary conditions. The analytical results obtained are compared with the direct numerical simulation of the Navier-Stokes equations. The secondary EOFs could transport scalar in larger areas and increase the scalar gradients, which significantly improve the mixing rate of scalars. It is shown that the heterogeneous zeta potentials could generate complex flow patterns and be used to enhance scalar mixing.

Flow Stress of Biomorphous Metal-Matrix Composites

For metal-matrix composites (MMCs), interfacial debonding between the ductile matrix and the reinforcing hard inclusions is an important failure mode. A fundamental approach to improving the properties of MMCs is to optimize their microstructure to achieve maximum strength and toughness. Here, we investigate the flow stress of a MMC with a nanoscale microstructure similar to that of bone. Such a 'biomorphous' MMC would be made of staggered hard and slender nanoparticles embedded in a ductile matrix. We show that the large aspect ratio and the nanometer size of inclusions in the biomorphous MMC lead to significantly improved properties with increased tolerance of interfacial damage. In this case, the partially debonded inclusions continue to carry mechanical load transferred via longitudinal shearing of the matrix material between neighboring inclusions. The larger the inclusion aspect ratio, the larger is the flow stress and work hardening rate for the composite. Increasing the volume concentration of inclusion also makes the biomorphous MMC more tolerant of interfacial damage.

Experimental Study on Vortex-Induced Vibrations of a Pipeline Near a Plane Boundary in Steady Flow

对单向水流作用下近壁管道横向涡激振动进行了实验模拟,重点探讨了管道与壁面间隙比(e/D)对管道涡激振动幅值和涡激振动频率响应特性的影响规律.实验结果表明,管道与壁面间隙宽度对管道涡激振动特性有较明显影响.在较大间隙比(e/D＞0.66)下,管道振幅随着Vr数的增大先快速增长到最大值,然后平缓下降;在振动初期(即Vr数较小时),管道振动频率变化基本符合Strouhal规律;在振动中后期(即Vr数较大时),管道振动频率变化不符合Strouhal规律,而在管道固有频率附近缓慢增长.在较小间隙比(e/D＜0.30)下,管道振幅随Vr数的增大先平缓上升到最大值,随后较快速下降;在振动初期,管道振动频率变化不遵循Strouhal规律;在整个振动范围内,与较大间隙比情况相比,随着Vr数增加,管道振动频率增长幅度明显较大.

Underwater Acoustics and Cavitating Flow of Water Entry

The fluid mechanics of water entry is studied through investigating the underwater acoustics and the supercavitation. Underwater acoustic signals in water entry are extensively measured at about 30 different positions by using a PVDF needle hydrophone. From the measurements we obtain (1) the primary shock wave caused by the impact of the blunt body on free surface; (2) the vapor pressure inside the cavity; (3) the secondary shock wave caused by pulling away of the cavity from free surface; and so on. The supercavitation induced by the blunt body is observed by using a digital high-speed video camera as well as the single shot photography. The periodic and 3 dimensional motion of the supercavitation is revealed. The experiment is carried out at room temperature.