A film flow model for analysing gravity-driven,thin wavy fluid films

To analyse the physics underlying gravity-driven runoff of thin wavy films, a film flow model is developed, and is solved with computational fluid dynamics. This model is based on the lubrication theory, and takes into account the gravitational, wall shear and surface tension forces. A key characteristic of the model is that it assumes only one computational cell over the film height, which enables studying film flow on larger computational domains. A main aim of this study is to perform a detailed validation of the numerical model. The film flow model is validated against several experiments of gravity-driven, thin fluid films on smooth surfaces. The time-averaged film thickness and the fluid speed profiles predicted by the model show very good agreement with experimental results. Similarly, the film flow model is able to predict the wave speeds with sufficient accuracy. The energy spectra of the waves, where higher frequency waves are present in film flows at higher Reynolds numbers, show an exponentially decaying trend at these high frequencies. The model performs better than the Nusselt equation for film flows, which under-predicts the time-averaged film thickness and over-predicts the time-averaged fluid speeds, even for flows at low Reynolds numbers. The film flow model is compared qualitatively for fingering behaviour. This model also allows to investigate film flows on large surfaces, which can be rough, curved and of complex geometrical shape.     

Reynolds stress modelling of rectangular open‐channel flow

A Reynolds stress model for the numerical simulation of uniform 3D turbulent open‐channel flows is described. The finite volume method is used for the numerical solution of the flow equations and transport equations of the Reynolds stress components. The overall solution strategy is the SIMPLER algorithm, and the power‐law scheme is used to discretize the convection and diffusion terms in the governing equations. The developed model is applied to a flow at a Reynolds number of 77000 in a rectangular channel with a width to depth ratio of 2. The simulated mean flow and turbulence structures are compared with measured and computed data from the literature. The computed flow vectors in the plane normal to the streamwise direction show a small vortex, called inner secondary currents, located at the juncture of the sidewall and the free surface as well as the free surface and bottom vortices. This small vortex causes a significant increase in the wall shear stress in the vicinity of the free surface. A budget analysis of the streamwise vorticity is carried out. It is found that both production terms by anisotropy of Reynolds normal stress and by Reynolds shear stress contribute to the generation of secondary currents. Copyright © 2006 John Wiley & Sons, Ltd.     

On the stability of plane parallel viscoelastic shear flows in the limit of infinite Weissenberg and Reynolds numbers

Elastic effects on the hydrodynamic instability of inviscid parallel shear flows are investigated through a linear stability analysis. We focus on the upper convected Maxwell model in the limit of infinite Weissenberg and Reynolds numbers. We study the effects of elasticity on the instability of a few classes of simple parallel flows, specifically plane Poiseuille and Couette flows, the hyperbolic-tangent shear layer and the Bickley jet.The equation for stability is derived and solved numerically using the spectral Chebyshev collocation method. This algorithm is computationally efficient and accurate in reproducing the eigenvalues. We consider flows bounded by walls as well as flows bounded by free surfaces. In the inviscid, nonelastic case all the flows we study are unstable for free surfaces. In the case of wall bounded flow, there are instabilities in the shear layer and Bickley jet flows. In all cases, the effect of elasticity is to reduce and ultimately suppress the inviscid instability.     

Wall shear stress and void fraction in Poiseuille bubbly flows: Part I: simple analytic predictions

Upward, co-current bubbly flows in a vertical rectangular duct are investigated at low liquid Reynolds numbers. The conditions considered are such that the pseudo-turbulent stresses remain negligible compared to the viscous stresses. The void fraction transverse distribution is idealised as step-functions and is then inserted in the conservation equations supplemented by appropriate closure laws. Analytical expressions are then obtained for the axial velocity profiles, for the lineic gas fraction and for the wall friction. The sensitivity of these quantities to the void distribution, characterised by the void fraction and the width of the three layers introduced, is discussed. It is shown that differential buoyancy effects govern the modification of the liquid velocity profiles. Notably, void peaking near walls is able to induce a wall shear stress many times higher than its single-phase flow counterpart at the same liquid flow rate. Also, the presence of a near wall region free of gas favours the onset of downward directed secondary flows. All these features correspond to experimental observations, and a few quantitative comparisons are also presented which support the validity of the model even in case of void coring. A companion paper (part II) will be devoted to systematic comparisons between predictions and experiments in the case of axisymmetric Poiseuille bubbly flows.     

Stability and nonlinear wavy regimes in downward film flows on a corrugated surface

The linear and nonlinear stability of downward viscous film flows on a corrugated surface to freesurface perturbations is analyzed theoretically. The study is performed with the use of an integral approach in ranges of parameters where the calculated results and the corresponding solutions of Navier-Stokes equations (downward wavy flow on a smooth wall and waveless flow along a corrugated surface) are in good agreement. It is demonstrated that, for moderate Reynolds numbers, there is a range of corrugation parameters (amplitude and period) where all linear perturbations of the free surface decay. For high Reynolds numbers, the waveless downward flow is unstable. Various nonlinear wavy regimes induced by varying the corrugation amplitude are determined. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 1, pp. 110–120, January–February, 2007.     

Measurements of velocity,temperature and velocity fluctuation distributions in falling liquid films

The velocity, temperature and velocity fluctuation distributions within falling spindle oil films in an inclined rectangular channel were measured using hot-wire techniques and thin thermocouples. The interfacial shear was caused by cocurrent air flow.The results indicate that the liquid films are as a whole much more laminar-like than turbulent in a range of Reynolds numbers (4γ/μ) up to the experimental limit of 6000. Mixing motion occurs in the vicinity of the interface; however, the flow near the wall surface exhibits no sign of such eddy motions, as predicted by the wall law for single phase turbulent flow. Although velocity fluctuation is observed within films with interfacial shear, mean velocity profiles are approximately the same as those obtained by the laminar film prediction.     

Supersonic dusty-gas flows with Knudsen effect in interphase momentum exchange

On the basis of the two-continuum model of dilute gas-solid suspensions, the dynamic behavior of inertial particles in supersonic dusty-gas flows past a blunt body is studied for moderate Reynolds numbers, when the Knudsen effect in the interphase momentum exchange is significant. The limits of the inertial particle deposition regime in the space of governing parameters are found numerically under the assumption of the slip and free-molecule flow regimes around particles. As a model problem, the flow structure is obtained for a supersonic dusty-gas point-source flow colliding with a hypersonic flow of pure gas. The calculations performed using the full Lagrangian approach for the near-symmetry-axis region and the free-molecular flow regime around the particles reveal a multi-layer structure of the dispersed-phase density with a sharp accumulation of the particles in some thin regions between the bow and termination shock waves. The project supported by the National Natural Science Foundation of China (90205024), and the Russian Foundation for Basic Research (RFBR grant No. 02-01-00770 and joint RFBR-NSFC grant No. 03-01-39004)     

深槽表面气体静压润滑低Reynolds数层流节流效应实验研究

气体介质在润滑间隙流动过程中,沿气体流动方向开设的表面矩形深槽结构内部产生旋涡阻碍气体流动,形成节流效应.通过矩形深沟槽表面静压润滑实验,开展深槽表面和光滑表面的流量特性对比测试,研究了低Reynolds数层流状态下深槽表面间隙的气体节流效应.结果表明:深槽表面产生显著的节流效应,并且随Reynolds数增加而增强;在低Reynolds数层流状态下,深槽表面可以产生节流效应,但是节流效应强度不稳定,当处于完全湍流状态时,节流效应维持定值;间隙尺寸、槽深等结构参数对节流效应影响明显.     

Wall shear stress modified by bubbles in a horizontal channel flow of silicone oil in the transition region

We performed laboratory experiments on bubbly channel flows using silicone oil, which has a low surface tension and clean interface to bubbles, as a test fluid to evaluate the wall shear stress modification for different regimes of bubble migration status. The channel Reynolds numbers of the flow ranged from 1000 to 5000, covering laminar, transition and turbulent flow regimes. The bubble deformation and swarms were classified as packing, film, foam, dispersed, and stretched states based on visualization of bubbles as a bulk void fraction changed. In the dispersed and film states, the wall shear stress reduced by 9% from that in the single-phase condition; by contrast, the wall shear stress increased in the stretched, packing, and foam states. We carried out statistical analysis of the time-series of the wall shear stress in the transition and turbulent-flow regimes. Variations of the PDF of the shear stress and the higher order moments in the statistic indicated that the injection of bubbles generated pseudo-turbulence in the transition regime and suppressed drag-inducing events in the turbulent regime. Bubble images and measurements of shear stress revealed a correlated wave with a time lag, for which we discuss associated to the bubble dynamics and effective viscosity of the bubble mixture in wall proximity.     

Measurement of liquid film thickness in micro square channel

In micro channels, slug flow becomes one of the main flow regimes due to strong surface tension. In micro channel slug flow, elongated bubble flows with the thin liquid film confined between the bubble and the channel wall. Liquid film thickness is an important parameter in many applications, e.g., micro heat exchanger, micro reactor, coating process etc. In the present study, liquid film thickness in micro square channels is measured locally and instantaneously with the confocal method. Square channels with hydraulic diameter of D_h = 0.3, 0.5 and 1.0 mm are used. In order to investigate the effect of inertial force on the liquid film thickness, three working fluids, ethanol, water and FC-40 are used. At small capillary numbers, liquid film at the channel center becomes very thin and the bubble interface is not axisymmetric. However, as capillary number increases, bubble interface becomes axisymmetric. Transition from non-axisymmetric to axisymmetric flow pattern starts from lower capillary number as Reynolds number increases. An empirical correlation for predicting axisymmetric bubble radius based on capillary number and Weber number is proposed from the present experimental data.     

Active Control of Wave Instabilities in Three-Dimensional Compressible Flows

In many fluid flows of practical importance transition is caused by the linear growth of wave instabilities, such as Tollmien–Schlichting waves, which eventually grow to a finite size at which stage secondary instabilities come into play. If transition is to be delayed or even avoided in such flows, then the linear growth of the disturbances must be prevented since control in the nonlinear regime would be a considerably more difficult task. Here a strategy for active control of two-dimensional incompressible and compressible Tollmien–Schlichting waves and its use in controlling the more practically relevant problem of crossflow instability which arises in swept-wing flows is discussed. The control is through an active suction/blowing distribution at the wall though the same result could be achieved by variable wall heating. In order to control the instability it is assumed that the wall shear stress and pressure are known from measurements. It is shown that, certainly at finite Reynolds numbers, it is sufficient to know the flow properties at a finite number of points along the wall. The cases of high and finite Reynolds numbers are discussed using asymptotic and numerical methods respectively. It is shown that a control strategy can be developed to stop the growth of all two-dimensional Tollmien–Schlichting waves at finite and large Reynolds numbers. Some discussion of nonlinear effects in the presence of active control is given and the possible control of other instability mechanisms investigated. Received 1 May 1998 and accepted 24 September 1998     

饱和蒸发和热毛细力作用下低雷诺数液膜在波纹竖壁上的流动

考虑表面蒸发压力和热毛细力作用情况下,对饱和蒸发状态下低雷诺数自由降落液膜在小波幅正弦型波纹壁面上的流动进行理论分析。对控制微分方程及边界条件进行量纲一化并引入流函数,对微分方程及边界条件进行摄动展开,得到了这种情况下液膜流动的简化分析模型,求出了近似解析解。讨论了壁面波纹、表面张力、蒸发压力、热毛细力对液膜流动的影响。研究表明:液膜的波动幅度随蒸发强度和热毛细力的增大而增大;液膜波动与壁面波纹的相位差随蒸发强度增大而增大,随热毛细力增大而减小。     

Viscous liquid film flow down an inclined corrugated surface. Calculation of the flow stability to arbitrary perturbations using an integral method

Viscous liquid film flow along an inclined corrugated (sinusoidal) surface has been studied. Calculations were performed using an integral model. The stability of nonlinear steady-state flows to arbitrary perturbations was examined using the Floquet theory. It has been shown that for each type of corrugation there is a critical Reynolds number for which unstable perturbations occur. It has been found that this value greatly depends on the physical properties of the liquid and geometric parameters of the flow. In particular, in the case of film flow down a smooth wall, the critical waveformation parameter depends only on the angle of inclination of the flow surface. The values of the corrugation parameters (amplitude and period) were obtained for which the film flow down a wavy wall is stable to arbitrary perturbations up to moderate Reynolds numbers. Such parameter values exist for all investigated angles of inclination of the flow surface.     

Effect of Inertial Terms on Fluid–Porous Medium Flow Coupling

The study considers an effect of the nonlinear inertial terms in the Brinkman filtration equation on the characteristics of coupled flows in a pure fluid and porous medium in the frameworks of two independent problems. The first problem is the forced boundary-layer flow overlying the Darcy–Brinkman porous medium. The Prandtl theory is used, and the self-similar equations are built to describe it. It is shown that the inertial terms have a valuable effect on the boundary-layer structure because of the large velocity gradient in the transition zone. The boundary-layer thickness in a porous medium rapidly grows at large Reynolds numbers. The velocity magnitude and gradient at the interface also change. The second independent problem is an analysis of the inertial terms effect on the flow stability. The neutral curves of the full and linearized flow models are built using the shooting method. They have different short-wave asymptotic, but there are no significant changes in the critical Reynolds numbers and corresponding wave numbers.     

Hydrodynamic force and torque models for a particle moving near a wall at finite particle Reynolds numbers

This research work is aimed at proposing models for the hydrodynamic force and torque experienced by a spherical particle moving near a solid wall in a viscous fluid at finite particle Reynolds numbers. Conventional lubrication theory was developed based on the theory of Stokes flow around the particle at vanishing particle Reynolds number. In order to account for the effects of finite particle Reynolds number on the models for hydrodynamic force and torque near a wall, we use four types of simple motions at different particle Reynolds numbers. Using the lattice Boltzmann method and considering the moving boundary conditions, we fully resolve the flow field near the particle and obtain the models for hydrodynamic force and torque as functions of particle Reynolds number and the dimensionless gap between the particle and the wall. The resolution is up to 50 grids per particle diameter. After comparing numerical results of the coefficients with conventional results based on Stokes flow, we propose new models for hydrodynamic force and torque at different particle Reynolds numbers. It is shown that the particle Reynolds number has a significant impact on the models for hydrodynamic force and torque. Furthermore, the models are validated against general motions of a particle and available modeling results from literature. The proposed models could be used as sub-grid scale models where the flows between particle and wall can not be fully resolved, or be used in Lagrangian simulations of particle-laden flows when particles are close to a wall instead of the currently used models for an isolated particle.     

Effects of varying Reynolds number, blowing ratio, and internal geometry on trailing edge cutback film cooling

Three-dimensional mean velocity and concentration fields have been measured for a water flow in a pressure side cutback trailing edge film cooling geometry consisting of rectangular film cooling slots separated by tapered lands. Three-component mean velocities were measured with conventional magnetic resonance velocimetry, while time-averaged concentration distributions were measured with a magnetic resonance concentration technique for flow at two Reynolds numbers (Re) differing by a factor of 2, three blowing ratios, and with and without an internal pin fin array in the coolant feed channel. The results show that the flows are essentially independent of Re for the regime tested in terms of the film cooling surface effectiveness, normalized velocity profiles, and normalized mean streamwise vorticity. Blowing ratio changes had a larger effect, with higher blowing ratios resulting in surface effectiveness improvements at downstream locations. The addition of a pin fin array within the slot feed channel made the spanwise distribution of coolant at the surface more uniform. Results are compared with transonic experiments in air at realistic density ratios described by Holloway et al. (2002a).     

Large eddy simulation of boundary layer flow under cnoidal waves

Water waves in coastal areas are generally nonlinear, exhibiting asymmetric velocity profiles with different amplitudes of crest and trough. The behaviors of the boundary layer under asymmetric waves are of great significance for sediment transport in natural circumstances. While previous studies have mainly focused on linear or symmetric waves, asymmetric wave-induced flows remain unclear, particularly in the flow regime with high Reynolds numbers.Taking cnoidal wave as a typical example of asymmetric waves, we propose to use an infinite immersed plate oscillating cnoidally in its own plane in quiescent water to simulate asymmetric wave boundary layer. A large eddy simulation approach with Smagorinsky subgrid model is adopted to investigate the flow characteristics of the boundary layer. It is verified that the model well reproduces experimental and theoretical results. Then a series of numerical experiments are carried out to study the boundary layer beneath cnoidal waves from laminar to fully developed turbulent regimes at high Reynolds numbers, larger than ever studied before.Results of velocity profile, wall shear stress, friction coefficient, phase lead between velocity and wall shear stress, and the boundary layer thickness are obtained. The dependencies of these boundary layer properties on the asymmetric degree and Reynolds number are discussed in detail.     

Convective inertia effects in wall-bounded thin film flows

We apply the boundary layer equations to inertial flow in wall bounded films that might be characterized as ‘thin’, say ɛ ≤ 0.1 where ɛ is the ratio of the characteristic lengths, yet to which the lubrication approximation of Reynolds no longer applies. Two particular flow geometries are investigated, nominally parallel plates and nominally inclined plates, both with and without spatially periodic perturbation of the stationary plate. A Galerkin-B spline formulation of the governing equations is employed, and we rely on parametric continuation to obtain solutions at higher values of the Reynolds number. In particular, we are able to demonstrate that the boundary layer equations yield accurate results for a wide range of Reynolds numbers when the aspect ratio is less than 1/10. We also find that in both nominally parallel and nominally inclined geometries the sign of the inertial force correction is determined by the film contour in the neighborhood of the exit, this result might have implications in the design of MEMS devices.     

New interpretation of specific sign of Reynolds stress in the boundary layer on a flat plate

The specific sign of Reynolds stress in the boundary layer on a flat plate at zero incidence is newly interpreted in present paper based on the theory of vortex-induced vortex. It avoids some problems appeared in a traditional explanation, on the basis of relationship between mean and fluctuating flows due to the transport of momentum. Through the analysis of local flow field in the immediate neighborhood of wall, the characteristics of Reynolds stress are identified through introducing turbulence-induced small-scale streamwise eddies above the flat plate. The positive Reynolds stress is theoretically verified. And such new interpretation illustrates that the generation of Reynolds stress, as well as fluctuating velocity, is intrinsically independent of the mean flow. But its specific sign would be determined by the mean flow due to the inertial forces. Other features,such as the intensity relationship among three components of fluctuating velocity, are also presented.     

A Dissipation Rate Equation for Low-Reynolds-Number and Near-Wall Turbulence

Two-equation turbulence models are usually formulated for specific flow types and are seldom validated against a variety of flows to account for near-wall and low-Reynolds-number effects simultaneously. In addition to low-Reynolds-number effects, near-wall flows also experience wall blocking, which is absent in free flows. Consequently, near-wall modifications to two-equation models could be quite different from low-Reynolds-number corrections. Besides, it is known that existing two-equation models perform poorly when used to calculate plane wall jets and two-dimensional backstep flows. These problems could be traced to the modeling of the dissipation rate equation. In this paper an attempt is made to improve the modeling of the dissipation rate equation so that it could successfully predict both free and wall-bounded shear flows including plane wall jets and backstep flows. The predictions are compared with experimental and direct numerical simulation data whenever available. Most of the data used are obtained at low Reynolds numbers. Good correlation with data is obtained. Therefore, for the first time, a model capable of correctly predicting free and wall-bounded shear flows, backstep flows, and plane wall jets is available. Received: 12 December 1995 and accepted 12 November 1996