Free surface deformation due to a source and a sink of equal strength in Stokes flow

Two-dimensional Stokes flow due to a source and a sink of equal strength below the free surface is analyzed and free surface shape and cusp formation are discussed. The source-sink pair below the free surface are aligned vertical to the free surface. In the analysis, the Stokes' approximation is used and surface tension effects are included, but gravity is neglected. The solution is obtained by using conformal mapping and complex function theory. From the solution, typical free surface shapes are shown and formation of a cusp on the free surface is discussed. As the capillary number increases, the converging free surface shape becomes singular and tends to form a cusp for sufficiently large capillary number. Typically, streamline patterns for some capillary numbers are also shown. As the capillary number vanishes, the solution is reduced to the linearized potential flow solution.     

Flow above the free end of a surface-mounted finite-height circular cylinder: A review

The wake of a surface-mounted finite-height circular cylinder and the associated vortex patterns are strongly dependent on the cylinder aspect ratio and the thickness of the boundary layer on the ground plane relative to the dimensions of the cylinder. Above a critical aspect ratio, the mean wake is characterized by streamwise tip vortex structures and Kármán vortex shedding from the sides of the cylinder. Below a critical aspect ratio, a unique mean wake structure is observed. Recent experimental studies in the literature that used phase-averaged techniques, as well as recent numerical simulations, have led to an improved physical understanding of the near-wake vortex flow patterns. However, the flow above the free end of the finite circular cylinder, and its relationship to the near wake, has not been systematically studied. The effects of aspect ratio and boundary layer thickness on the free-end flow field are also not completely understood, nor has the influence of Reynolds number on the free-end flow field been fully explored. Common features associated with the free end include separation from the leading edge, a mean recirculation zone containing a prominent cross-stream arch (or mushroom) vortex, and reattachment onto the free-surface. Other flow features that remain to be clarified include a separation bubble near the leading edge, one or two cross-stream vortices within this separation bubble, the origins of the streamwise tip or trailing vortices, and various critical points in the near-surface flow topology. This paper reviews the current understanding of the flow above the free end of a surface-mounted finite-height circular cylinder, with a focus on models of the flow field, surface oil flow visualization studies, pressure and heat flux distributions on the free-end surface, measurements of the local velocity field, and numerical simulations, found in the literature.     

Instantaneous flow field above the free end of finite-height cylinders and prisms

The flow above the free ends of surface-mounted finite-height circular cylinders and square prisms was studied experimentally using particle image velocimetry (PIV). Cylinders and prisms with aspect ratios of AR = 9, 7, 5, and 3 were tested at a Reynolds number of Re = 4.2 × 10⁴. The bodies were mounted normal to a ground plane and were partially immersed in a turbulent zero-pressure-gradient boundary layer, where the boundary layer thickness relative to the body width was δ/D = 1.6. PIV measurements were made above the free ends of the bodies in a vertical plane aligned with the flow centreline. The present PIV results provide insight into the effects of aspect ratio and body shape on the instantaneous flow field. The recirculation zone under the separated shear layer is larger for the square prism of AR = 3 compared to the more slender prism of AR = 9. Also, for a square prism with low aspect ratio (AR = 3), the influence of the reverse flow over the free end surface becomes more significant compared to that for a higher aspect ratio (AR = 9). For the circular cylinder, a cross-stream vortex forms within the recirculation zone. As the aspect ratio of the cylinder decreases, the reattachment point of the separated flow on the free end surface moves closer to the trailing edge. For both the square prism and circular cylinder cases, the instantaneous velocity vector field and associated in-plane vorticity field revealed small-scale structures mostly generated by the separated shear layer.     

Free surface deformation due to a source or a sink in Stokes flow

Free surface shape and cusp formation are analyzed by considering two-dimensional viscous flow due to a line source or a line sink below the free surface where the strength of source/sink is given arbitrarily. In the analysis, the Stokes' approximation is used and surface tension effects are included, but gravity is neglected. The solution is obtained analytically by using conformal mapping and complex function theory. From the solution, shapes of the free surface are shown and the formation of a cusp on the free surface is discussed. As the capillary number decreases in negative, the free surface shape becomes singular and in a real fluid a cusp should form on the free surface below some negative critical capillary number. Typically, streamline patterns for some capillary numbers are also shown. As the small capillary number vanishes, the solution is reduced to a linearized potential flow solution.     

Residual mass and flow regimes for the immiscible liquid–liquid displacement in a plane channel

The motion of two immiscible liquids in a plane channel is analyzed for the case in which the flow conditions and the interactions between the liquids and the solid surface maintain the displaced fluid attached to the wall. The Galerkin Finite Element Method is used to compute the velocity field and the configuration of the interface between the two fluids. We compare the residual mass fraction left on the wall with its two counterparts in capillary tubes, namely residual mass fraction and dimensionless layer thickness of the displaced fluid. The main result of this comparison was that although there is a qualitative similarity concerning the layer thickness between the two cases, the residual fraction of mass presented an important difference, showing that when the aspect ratio of the capillary passage is large there is an increase in the displacement efficiency. The thickness of the displaced liquid film attached to the channel walls is a function of the capillary number (Ca) and the viscosity ratio (N_μ). A map of streamlines in the Cartesian space (Ca, N_μ) with the different flow regimes of the problem is presented. We also showed that we can adapt the available analytical results obtained for gas-displacement in capillary tubes to the plane channel case, for low values of Ca.     

Three‐dimensional transient free‐surface flow of viscous fluids inside cavities of arbitrary shape

The three‐dimensional transient free‐surface flow inside cavities of arbitrary shape is examined in this study. An adaptive (Lagrangian) boundary‐element approach is proposed for the general three‐dimensional simulation of confined free‐surface flow of viscous incompressible fluids. The method is stable as it includes remeshing capabilities of the deforming free‐surface, and thus can handle large deformations. A simple algorithm is developed for mesh refinement of the deforming free‐surface mesh. Smooth transition between large and small elements is achieved without significant degradation of the aspect ratio of the elements in the mesh. The method is used to determine the flow field and free‐surface evolution inside cubic, rectangular and cylindrical containers. These problems illustrate the transient nature of the flow during the mixing process. Surface tension effects are also explored. Copyright © 2003 John Wiley & Sons, Ltd.     

Mixed Convection in a Darcy–Brinkman Porous Medium with a Constant Convective Thermal Boundary Condition

The characteristics of the boundary layer flow past a plane surface adjacent to a saturated Darcy–Brinkman porous medium are investigated in this paper. The flow is driven by an external free stream moving with constant velocity. The surface is heated with a convective boundary condition with constant heat transfer coefficient. The problem is non-similar and is investigated numerically by a finite difference method. The problem is governed by four non-dimensional parameters, that is, the convective Darcy number, the convective Grashof number, the Prandtl number, and the axial distance along the plate. The influence of these parameters on the results is investigated, and the results are presented in tables and figures. The Darcy term and the Grashof term in the momentum equation contradict each other and this contradiction makes the problem complicated. However, the wall shear stress and the wall temperature increase continuously along the plate and the wall temperature always tends to 1.     

Double-diffusive parallel flow induced in a horizontal Brinkman porous layer subjected to constant heat and mass fluxes: analytical and numerical studies

Fluid flow and heat and mass transfer induced by double-diffusive natural convection in a horizontal porous layer subjected to vertical gradients of temperature and concentration are studied analytically and numerically using the Brinkman-extended Darcy model. Both cases of rigid and free horizontal boundaries are examined in the present study. The parameters governing the problem are the Rayleigh number R_T, the Lewis number Le, the buoyancy ratio N, the Darcy number Da and the aspect ratio Ar. The analytical solution is based on the parallel flow approximation. The critical Rayleigh number corresponding to the onset of the parallel flow in this system is determined analytically as a function of Le, N and Da. For sufficiently small Da, both free and rigid boundaries yield results which are identical to those predicted by the Darcy model. The present investigation shows that there exists a region in the plane (N, Le) where the convective flow is not possible in the layer regardless of the Rayleigh and Darcy numbers considered. Received on 21 December 1998     

Effect of span length and temperature on the 3-D confined flow around a vortex promoter

This article presents a numerical study on the influence of span length and wall temperature on the 3-D flow pattern around a square section vortex promoter located inside a micro-channel in the low Reynolds number regime. The first objective of the work is to quantify the critical Reynolds number that defines the onset of vortex shedding and to identify the different regimes that appear as a function of the channel aspect ratio (span to height ratio). We found that the critical Reynolds number for the onset of the Karman street regime increases as the aspect ratio decreases. In particular, for the aspect ratio of 1/2 the critical Reynolds number is nearly six times the critical Reynolds number of the 2-D problem. An intermediate oscillating regime between the steady and the Karman street solutions was also found to exist within a rather wide range of Reynolds numbers for small channel aspect ratios. The second objective was to investigate the influence of the vortex promoter wall temperature on both vortex shedding and flow pattern. This has practical engineering implications because the working fluid considered in the article is water that has a viscosity that depends significantly on temperature and promotes a strong coupling between the momentum and energy equations that influences the system behaviour. Results indicate that high surface temperature on the prism promotes the onset of the Karman street, suggesting design guidelines for micro-channel based heat sinks that make use of vortex promoters.     

Near-wake flow structure of elliptic cylinders close to a free surface: effect of cylinder aspect ratio

The flow fields behind elliptic cylinders adjacent to a free surface were investigated experimentally in a circulating water channel. A range of cylinder aspect ratios (AR=2, 3, 4) were considered, while the cross-sectional area of the elliptical cylinder was kept constant. The main objective of this study was to investigate the effect of cylinder aspect ratio and a free surface on the flow structure in the near-wake behind elliptic cylinders. For each elliptic cylinder, the flow structure was analyzed for various values of the submergence depth of the cylinder beneath the free surface. The flow fields were measured using a single-frame double-exposure PIV (Particle Image Velocimetry) system. For each experimental condition, 350 instantaneous velocity fields were obtained and ensemble-averaged to obtain the mean velocity field and spatial distribution of the mean vorticity statistics. The results show that near-wake can be classified into three typical flow patterns: formation of a Coanda flow, generation of substantial jet-like flow, and attachment of this jet flow to the free surface. The general flow structure observed behind the elliptic cylinders resembles the structure previously reported for a circular cylinder submerged near a free surface. However, the wake width and the angle of downward deflection of the shear layer developed from the lower surface of the elliptic cylinder differ from those observed for a circular cylinder. These trends are enhanced as cylinder aspect ratio is increased. In addition, the free surface distortion is also discussed in the paper.     

Momentum and heat transfer in the magnetohydrodynamic stagnation-point flow of a viscoelastic fluid toward a stretching surface

An analysis is made of the steady two-dimensional stagnation-point flow of an incompressible viscoelastic fluid over a flat deformable surface when the surface is stretched in its own plane with a velocity proportional to the distance from the stagnation-point. It is shown that for a viscoelastic conducting fluid of short memory (obeying Walters’ Bʹ model), a boundary layer is formed when the stretching velocity of the surface is less than the inviscid free-stream velocity and velocity at a point increases with increase in the Hartmann number. On the other hand an inverted boundary layer is formed when the surface stretching velocity exceeds the velocity of the free stream and the velocity decreases with increase in the Hartmann number. A novel result of the analysis is that the flow near the stretching surface is that corresponding to an inviscid stagnation-point flow when the surface stretching velocity is equal to the velocity of the free stream. Temperature distribution in the boundary layer is found when the surface is held at constant temperature and surface heat flux is determined. It is found that in the absence of viscous and Ohmic dissipation and strain energy in the flow, temperature at a point decreases with increase in the Hartmann number.     

Natural convection flow of a viscous incompressible fluid in a rectangular porous cavity heated from below with cold sidewalls

We consider the flow, which is induced by differential heating on the boundaries of a porous cavity heated from below. In particular we allow the sidewalls to have the same cold temperature as the upper surface, and thus the problem is a variant of the Darcy-Bénard convection problem, but one where there is flow at all non-zero Grashof numbers. Attention is focused on how the flow and heat transfer is affected by variations in the cavity aspect ratio, the Grashof number and the Darcy number. The flow becomes weaker as the Darcy number decreases from the pure fluid limit towards the Darcy-flow limit. In addition the number of cells which form in the cavity varies primarily with the aspect ratio and is always even due to the symmetry imposed by the cold sidewalls.     

Three-dimensional instability of two counter-rotating vortices in a rectangular cavity driven by parallel wall motion

The linear stability of two counter-rotating vortices driven by the parallel motion of two facing walls in a rectangular cavity is investigated by a finite volume method. Critical Reynolds and wave numbers are calculated for aspect ratios ranging from 0.1 to 5. This range is sufficient to find the asymptotic behavior of the critical parameters when the aspect ratio tends to zero and infinity, respectively. The critical curve is smooth for all aspect ratios and, hence, the character of the instability changes continuously. When the moving walls are far apart the mechanism is centrifugal, as in the classical lid-driven cavity. For aspect ratios near unity a combined mechanism, also involving strain near the vortex cores, leads to the instability which tends to asymmetrically displace the vortex cores, very similar to the cooperative short-wave instability of a free counter-rotating vortex pair. In the limit when plane Poiseuille flow is approached in the bulk, the three-dimensional perturbations are strongly localized near both downstream ends of the moving walls.     

Measurements of the flow over a low-aspect-ratio cylinder mounted on a ground plane

The flow over a finite-height cylinder of aspect ratio 1, with one end mounted on a ground plane and the other end free, has been studied by means of surface flow visualisation, particle image velocimetry (PIV) and surface pressure measurements. The diameter-based Reynolds number was 200,000. The mean flow topology has been identified in three areas: the horseshoe vortex system, the separated flow over the free-end and the wake region. Evidence is shown for the existence of a three-horseshoe vortex system, while the mean flow over the free-end consists of an arch vortex with its bases on the forward half of the free-end. There are two tip vortices coming off the free-end. The wake region is found to be highly unsteady, with considerable variation from the mean flow.     

A low‐dimensional spectral approach for transient axisymmetric free‐surface flow inside thin cavities of arbitrary shape

A low‐dimensional spectral method is used to solve the transient axisymmetric free surface flow inside thin cavities of arbitrary shape. The flow field is obtained on the basis of the lubrication equations, which are expanded in terms of orthonormal functions over the cavity gap. The formulation accounts for nonlinearities stemming from inertia and front location. The work is of close relevance to the filling stage during die casting, and injection molding, or the flow inside annular (extrusion) dies. Both flows under an imposed flow rate, and an imposed pressure at the cavity entrance are examined. The influence of inertia, aspect ratio, gravity, and wall geometry on the evolution of the front, flow rate, and pressure is assessed particularly in the early stage of flow, when a temporal behavior of the ‘boundary‐layer’ type develops. The multiple‐scale method is applied to obtain an approximate solution at small Reynolds number, Re. Comparison with the exact (numerical) solution indicates a wide range of validity for the multiple‐scale approach, including the moderately small Re range. Copyright © 2002 John Wiley & Sons, Ltd.     

Low-Reynolds-number motion of a deformable drop between two parallel plane walls

The motion of a three-dimensional deformable drop between two parallel plane walls in a low-Reynolds-number Poiseuille flow is examined using a boundary-integral algorithm that employs the Green’s function for the domain between two infinite plane walls, which incorporates the wall effects without discretization of the walls. We have developed an economical calculation scheme that allows long-time dynamical simulations, so that both transient and steady-state shapes and velocities are obtained. Results are presented for neutrally buoyant drops having various viscosity, size, deformability, and channel position. For nearly spherical drops, the decrease in translational velocity relative to the undisturbed fluid velocity at the drop center increases with drop size, proximity of the drop to one or both walls, and drop-to-medium viscosity ratio. When deformable drops are initially placed off the centerline of flow, lateral migration towards the channel center is observed, where the drops obtain steady shapes and translational velocities for subcritical capillary numbers. With increasing capillary number, the drops become more deformed and have larger steady velocities due to larger drop-to-wall clearances. Non-monotonic behavior for the lateral migration velocities with increasing viscosity ratio is observed. Simulation results for large drops with non-deformed spherical diameters exceeding the channel height are also presented.     

Viscous effects in liquid encapsulated liquid bridges

An analytical derivation of the surface deflections and the streamfunctions for the flow inside a liquid encapsulated liquid bridge has been derived using an asymptotic expansion about a small capillary number. The model assumes an initially flat and cylindrical interface under the assumption that the densities of both fluids are equal. To simplify the analysis, the top and bottom walls are assumed to be stress-free and the Reynolds number is assumed to be negligible. Flow is generated either by a moving outer wall (shear-driven flow) or by applying a temperature difference across the top and bottom walls (Marangoni-driven flow). The resulting equations show that for the shear-driven flow, as the viscosity ratio increases, the surface deflections increase monotonically. For the Marangoni-driven flow there exist values of the viscosity ratio where the surface deflections reach a minimum and then switch signs. This investigation shows that it may be possible in more realistic systems to use an outer encapsulating liquid of the proper viscosity ratio to stabilize the liquid–liquid interface during float zone crystal growth.     

Dynamics of viscoelastic liquid filaments: Low capillary number flows

Many applications of viscoelastic free surface flows requiring formation of drops from small nozzles, e.g., ink-jet printing, micro-arraying, and atomization, involve predominantly extensional deformations of liquid filaments. The capillary number, which represents the ratio of viscous to surface tension forces, is small in such processes when drops of water-like liquids are formed. The dynamics of extensional deformations of viscoelastic liquids that are weakly strain hardening, i.e., liquids for which the growth in the extensional viscosity is small and bounded, are here modeled by the Giesekus, FENE-P, and FENE-CR constitutive relations and studied at low capillary numbers using full 2D numerical computations. A new computational algorithm using the general conformation tensor based constitutive equation [M. Pasquali, L.E. Scriven, Theoretical modeling of microstructured liquids: a simple thermodynamic approach, J. Non-Newtonian Fluid Mech. 120 (2004) 101–135] to compute the time dependent viscoelastic free surface flows is presented. DEVSS-TG/SUPG mixed finite element method [M. Pasquali, L.E. Scriven, Free surface flows of polymer solutions with models based on conformation tensor, J. Non-Newtonian Fluid Mech. 108 (2002) 363–409] is used for the spatial discretization and a fully implicit second-order predictor–corrector scheme is used for the time integration. Inertia, capillarity, and viscoelasticity are incorporated in the computations and the free surface shapes are computed along with all the other field variables in a fully coupled way. Among the three models, Giesekus filaments show the most drastic thinning in the low capillary number regime. The dependence of the transient Trouton ratio on the capillary number in the Giesekus model is demonstrated. The elastic unloading near the end plates is investigated using both kinematic [M. Yao, G.H. McKinley, B. Debbaut, Extensional deformation, stress relaxation and necking failure of viscoelastic filaments, J. Non-Newtonian Fluid Mech. 79 (1998) 469–501] and energy analyses. The magnitude of elastic unloading, which increases with growing elasticity, is shown to be the largest for Giesekus filaments, thereby suggesting that necking and elastic unloading are related.     

Formation of the free surface of a viscous fluid volume inside a rotating horizontal cylinder

Two-dimensional viscous flow with a free surface in a horizontal cylinder rotating at a constant speed is investigated numerically using the boundary element method. It is shown that in the initial stage of rotation of the cylinder four different variants of the behavior of the free surface can be realized in the stage of transition from horizontal to steady-state form.     

Three-dimensional effects on regular reflection in steady supersonic flows

The reflection of shock waves between two symmetrical wedges is investigated for the case of three-dimensional flows. Oblique shadowgraphs at various optical angles of yaw and pitch were used to examine the nature of fully three-dimensional flows, with wedge aspect ratios as low as 0.25 being considered. These images were used to construct surface models of the overall flow field for various reflection patterns and aspect ratios, which provides a visual indication of the flow field shape. For a Mach number of 3.1, and suitable wedge angles, the flow field with regular reflection on the tunnel centreline and Mach reflection further out is examined. The point of transition from regular reflection to the peripheral Mach surfaces is identified for various wedge angles and aspect ratios. It is shown that the transition points move outwards from the central plane as the aspect ratio decreases. This shows that three-dimensional flows favor regular reflection, because of the increasing curvature of the incident shock as the wedge becomes narrower, causing a decrease in the local angle of incidence. The height of the Mach stem is shown to be highly dependent on the geometry of the test wedge models. The Mach stem height decreases with aspect ratio due to the three-dimensional relieving effect, where the increase in lateral flow relieves the pressure over the surfaces of the wedges. Experimental evidence of the existence of the strong oblique shock solution in steady flows is presented.Received: 7 July 2003, Revised: 20 October 2003, Accepted: 6 November 2003, Published online: 10 February 2004PACS: 47.40.Nm Correspondence to: B.W. Skews