Exact solution for entrance region flow between parallel plates

A method of obtaining an exact solution of the complete Navier-Stokes equations for incompressible flow has been developed to study the development of fluid flow between parallel plates. An analytical solution has been developed in the form of a convergent infinite series and a numerical solution obtained for Reynolds numbers of 10 and 500. These results have been compared with data available elsewhere. The method developed here is expected to have general application     

Heat transfer in tube assemblies under conditions of laminar axial,transverse and inclined flow

This paper considers laminar flow heat transfer in tube assemblies. The main interest is focused on the virtually unexplored cases of heat transfer under conditions of fully-developed flow inclined to the axes of the tubes and of purely transverse developing flow. The limiting cases of purely axial or purely transverse fully-developed flow are also examined. In all cases, the thermal boundary condition on the tubes is constant heat flux. Governing differential equations are expressed in terms of curvilinear-orthogonal coordinates and solved using finite-differences. Results are compared with available theoretical and experimental data. The effect of the transverse component of the flow on the temperature distribution is found to remain very strong even in nearly-axial flows and therefore considerably higher heat transfer coefficients are exhibited by a nearly-axial flow than a purely axial one.     

Visualisation and modelling studies of churn flow in a vertical pipe

Churn flow is an important flow regime intermediate between slug flow and annular flow. A feature of this regime is the occurrence of very large waves travelling upwards over a liquid film substrate which may intermittently travel downwards. These waves are often formed close to the liquid inlet where their behaviour is usually difficult to observe. This paper describes a series of experiments using a test section with a specially constructed transparent liquid inlet. High-speed video recordings show clearly the process of wave formation and analysis of the recordings gave data on wave frequencies and typical velocities. Also, predictions of velocity and distance travelled by the waves were obtained via the application of a simple mathematical model that takes into account the forces acting on a circumferentially coherent wave.     

Improving the robustness of the control volume finite element method with application to multiphase porous media flow

Control volume finite element methods (CVFEMs) have been proposed to simulate flow in heterogeneous porous media because they are better able to capture complex geometries using unstructured meshes. However, producing good quality meshes in such models is nontrivial and may sometimes be impossible, especially when all or parts of the domains have very large aspect ratio. A novel CVFEM is proposed here that uses a control volume representation for pressure and yields significant improvements in the quality of the pressure matrix. The method is initially evaluated and then applied to a series of test cases using unstructured (triangular/tetrahedral) meshes, and numerical results are in good agreement with semianalytically obtained solutions. The convergence of the pressure matrix is then studied using complex, heterogeneous example problems. The results demonstrate that the new formulation yields a pressure matrix than can be solved efficiently even on highly distorted, tetrahedral meshes in models of heterogeneous porous media with large permeability contrasts. The new approach allows effective application of CVFEM in such models.     

Simulations for gas flows in microgeometries using the direct simulation Monte Carlo method

Micro gas flows are often encountered in MEMS devices and classical CFD could not accurately predict the flow and thermal behavior due to the high Knudsen number. Therefore, the gas flow in microgeometries was investigated using the direct simulation Monte Carlo (DSMC) method. New treatments for boundary conditions are verified by simulations of micro-Poiseuille flow, compared with the previous boundary treatments, and slip analytical solutions of the continuum theory. The orifice flow and the corner flow in microchannels are simulated using the modified DSMC codes. The predictions were compared with existing experimental phenomena as well as predictions using continuum theory. The results showed that the channel geometry significantly affects the microgas flow. In the orifice flow, the flow separation occurred at very small Reynolds numbers. In the corner flow, no flow separation occurred even with a high driving pressure. The DSMC results agreed well with existing experimental information.     

On the comparison of two different solutions in the form of series of the governing equation of an unsteady flow of a second grade fluid

In this paper, two different solutions in the form of series of the governing equation of unsteady flow of a second grade fluid are considered. These are series expansions with respect to inverse power of time and a perturbation expansion. Two illustrative examples are given. One of them is the unsteady flow of a second grade fluid over a plane wall suddenly set in motion and the other is the diffusion of a line vortex in a fluid of second grade. It is a remarkable fact that the expression of the series expansion with respect to inverse power of time is exactly in the same form as that of the perturbation expansion. Thus, it is possible to replace a series expansion with respect to inverse power of time with a perturbation expansion.     

气液混输悬链线立管系统两相流特性实验

严重段塞流是海洋工程气液混输管线--立管系统中常见的一种特殊有害流动现象, 采用水平--下倾--悬链线立管气液混输组合管道系统, 通过系列实验在悬链线立管中获得了严重段塞流、间歇流和震荡流等流型, 阐述了这些流动现象的形成机理, 提出了能够产生严重段塞流的判定准则. 结果表明, 悬链线立管严重段塞流具有明显周期性, 在一个周期内的流动特征可分为液塞形成、液体出流、液气喷发及液体回流等4个阶段, 进而给出了各阶段中相关流动参数的变化规律. 在实验中同时还对悬链线与垂直立管中严重段塞流形成机理进行了比较分析, 发现两者在液塞形成阶段有显著差别. 其中, 在悬链线立管中液塞形成之前首先需要经历一个气液混合液塞形成过程, 而垂直立管则没有这个过程.      

ON INITIATION OF GULLY DEBRIS FLOW BY LIMIT EQUILIBRIUM ANALYSIS1)

沟道泥石流启动的力学本质为沟道中的堆积体所受剪切力克服摩擦力的过程。本文通过堆积体临界状态的受力平衡分析确立颗粒堆积体的临界平衡条件,发现堆积结构稳定性受坡度、颗粒内摩擦系数、堆积高度和水流速度共同作用。在此基础上通过文献中的系列实验数据分析拟合得到了一个泥石流启动的临界流速条件,相关结果和文献中的两个系列实验结果吻合良好,有望为沟道泥石流启动提供一定的支撑。     

On the flows produced by sudden application of a constant pressure gradient or by impulsive motion of a boundary

Some properties of the time-dependent Navier-Stokes equations are discussed for flows impulsively started from rest by sudden application of a constant pressure gradient or by the impulsive motion of a boundary. Five illustrative examples are given. They are: unsteady flow in a circular cylinder moving parallel to its length, starting flow in a circular pipe, unsteady flow in a rotating cylinder, starting flow in a rectangular channel moving parallel to its length and unsteady flow in a channel of rectangular cross-section. It is found that the expressions of the quantities such as velocity, flux and skin friction are in series forms which may be rapidly convergent for large values of the time but slowly convergent for small values of the time or vice versa. It is shown that if their expressions can be found for one of large values of the time or small values of the time, these expressions can be used for the other.     

高速三维边界层的横流不稳定性

用两点四阶差分格式研究旋转圆锥超音速三维边界层的横流不稳定性和壁面冷却对稳定性的影响数值结果表明，与二维边界层相比横流使三维边界层第一模式增长率增大，对第二模式影响很小；Ｍｅ＜４３第一模式最不稳定，Ｍｅ＞４３第二模式最不稳定；三维边界层最不稳定第二模式是三维波，二维边界层则为二维波；壁面冷却对第一模式起稳定作用，对第二模式起不稳定作用     

A THREE-FLUID MODEL OF THE SAND-DRIVEN FLOW

ＡＴＨＲＥＥ－ＦＬＵＩＤＭＯＤＥＬＯＦＴＨＥＳＡＮＤ－ＤＲＩＶＥＮＦＬＯＷ￥(刘大有，董飞)ＬｉｕＤａｙｏｕ；ＤｏｎｇＦｅｉ（ＩｎｓｔｉｔｕｔｅｏｆＭｅｃｈａｎｉｃｓ，ＡｃａｄｅｍｉａＳｉｎｉｃａ，Ｂｅｉｊｉｎｇ１０００８０，Ｐ．Ｒ．Ｃｈｉｎａ）Ａｂｓ...     

An investigation of adiabatic spiral vortex flow in wide annular gaps by visualisation and digital analysis

This paper, the third in a series describing experimental investigations into spiral vortex flow, presents visual evidence illustrating adiabatic transition modes in a wide gap of radius ratio 0.848. Also, power spectra, relating to velocity fluctuations, are found to be comparable for two working fluids, oil and air. Good agreement has been found between these results and those published by the authors in which diabatic transition modes were related to the heat transfer characteristics of the flow at various axial Reynolds numbers     

Analytic solution for axisymmetric flow over a nonlinearly stretching sheet

An analysis is performed for the boundary-layer flow of a viscous fluid over a nonlinear axisymmetric stretching sheet. By introducing new nonlinear similarity transformations, the partial differential equations governing the flow are reduced to an ordinary differential equation. The resulting ordinary differential equation is solved using the homotopy analysis method (HAM). Analytic solution is given in the form of an infinite series. Convergence of the obtained series solution is explicitly established. The solution for an axisymmetric linear stretching sheet is obtained as a special case.     

An efficient method for solving the mixed direct-inverse problem of the transonic rotational flow in plane casades

The method in [1] has been extended to the case of rotational flow in this paper. A new method for dealing with the shock wave is presented. This method has the advantages of both the shock-fitting and the shock capturing methods. The direct problem and the mixed direct-inverse problem of the rotational flow in a transonic plane cascade at both design and off design conditions are solved, and the results show that the present method has rapid convergence rate and high accuracy even for the flow with moderately strong shocks. The calculations have been carried out on the DPS-8 computer, and for the direct problem, only 50–80 iterations are needed, and 50–80 seconds of CPU time are required.     

Simulation of particle/fluid flows in vertical circular pipes

A two fluid continuum model is applied to the simulation of steady fully developed particle/fluid flow in a vertical circular pipe. Both closed form and numerical solutions are obtained to the associated governing equations. These solutions are compared to the predictions of another two fluid continuum model. These comparisons are used to illustrate the differences in the predictions of two widely used models, even in the fundamental problem of steady flow in a circular pipe.     

Transient effects of orthogonal pipe oscillations on laminar developing incompressible flow

This paper presents a numerical study of the transient developing laminar flow of a Newtonian incompressible fluid in a straight horizontal pipe oscillating around the vertical diameter at its entrance. The flow field is influenced by the tangential and Coriolis forces, which depend on the through‐flow Reynolds number, the oscillation Reynolds number and the angular amplitude of the pipe oscillation. The impulsive start of the latter generates a transient pulsating flow, whose duration increases with axial distance. In any cross‐section, this flow consists of a pair of symmetrical counter‐rotating vortices, which are alternatively clockwise and anti‐clockwise. The circumferentially averaged friction factor and the axial pressure gradient fluctuate with time and are always larger than the corresponding values for a stationary pipe. On the other hand, local axial velocities and local wall shear stress can be smaller than the corresponding stationary pipe values during some part of the pipe oscillation. The fluctuation amplitude of these local variables increases with axial distance and can be as high as 50% of the corresponding stationary pipe value, even at short distances from the pipe entrance. Eventually, the flow field reaches a periodic regime that depends only on the axial position. The results show that the transient flow field depends on the pipe oscillation pattern (initial position and/or direction of initial movement). Copyright © 2000 John Wiley & Sons, Ltd.     

大气绕山体的运动及其输移特性的研究

对分层大气山体绕流的流动模式及扩散输移特性进行了数值模拟。采用隐式时间离散方法在贴体网格系统下求解雷诺平均的N-S方程,计算结果描述出大气流动的特征,证实了分层(以Froude数为特征参数,定义为F=U/NH,U为来流风速,N为Brunt-Vaisala频率,H为山体高度)变化对山体绕流流态的影响。数值结果表明:当Froude数大于4.0时,山体绕流的流动不再依赖于大气分层的变化。当Froude数介于4.0和1.0之间时,流场中出现了Lee波,并随着Froude数的进一步减小,流动分离发生及Lee波破碎现象。同时模型也预测了在各种流动模式下大气中夹杂着的污染物绕山体的传输特性,表明大气的分层现象对污染物的分布有着非常重要的影响。     

Numerical solution of the turbulent flow of an incompressible fluid in curved planar and axially symmetrical channels

The governing equations for axially symmetric flow, where the Reynolds stresses are expressed by scalar turbulent viscosity, are the Reynolds equations. The turbulence model k, ? is used in the well-known form for fully developed turbulent flow.The numerical method, a continuation of the MAC system¹, is adapted so that even for high Reynolds cell numbers precision (δx²) can be achieved for the steady flow. Irregular cells join the rectangular network on the curved surface. Von Neumann's stability condition of the linearised numerical system is investigated. Special problems concerning the numerical solution of the turbulence model equations are stated and a special procedure is worked out to ensure that the fields k, ? do not converge to physically meaningless values. The program for the computer is universal in that the boundary problems can be assigned by input data.As an example, an axially symmetrical diffuser with an area ratio of widening 1.40 is computed. Fields of velocity and pressure at the wall as well as fields v_T and k are assessed. The results are compared with an experiment. The conclusion is that this method is suitable for the problems mentioned in this study as well as for unsteady flow.     

Unsteady theory of the hydrofoil of finite span

A three-dimensional unsteady flow around the hydrofoil has been dealt with in this paper. Under the assumption of infinitesimal wave and thin hydrofoil, analytic solution of perturbation velocity potential induced by the unsteady flow of a hydrofoil which is simplified as a series of rectangular vortex rings is achieved and then the unsteady forces can be evaluated. The paper had been accepted by the X VI th International Congress of IUTAM, Lyngby, Denmark, Aug. 19–25, 1985.     

流动聚焦及射流不稳定性

流动聚焦是一种有效的微细射流产生方法,其原理可以描述为从毛细管流出的流体由另一种高速运动的流体驱动,经小孔聚焦后形成稳定的锥–射流结构,射流因不稳定性破碎成单分散的液滴.自从1998年流动聚焦被提出以来,陆续发展了单轴流动聚焦、电流动聚焦、复合流动聚焦和微流控流动聚焦等毛细流动技术.这些技术稳定、易操作、没有苛刻的环境条件的要求,能够制备单分散性较好的微纳米量级的液滴、颗粒和胶囊,在科学研究和实际应用中具有重要价值.流动聚焦涉及了多尺度、多界面和多场耦合的复杂流体力学问题,其中稳定的锥形是形成稳定射流的先决条件,过程参数是影响射流界面扰动发展的关键因素,而射流不稳定性分析是揭示射流破碎的最主要理论工具.该文回顾了近二十年来不同结构流动聚焦的研究进展,概述这些技术涉及的过程控制、流动模式、尺度律和不稳定性分析等关键力学问题,总结射流不稳定性的研究方法和已取得的成果,最后展望流动聚焦的研究方向和应用前景.