Onset of motion of a three-dimensional droplet on a wall in shear flow at moderate Reynolds numbers

2008 ◽  
Vol 599 ◽  
pp. 341-362 ◽  
Author(s):  
HANG DING ◽  
PETER D. M. SPELT
Keyword(s):  
Contact Angle ◽  
Reynolds Number ◽  
Steady State ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Critical Value ◽  
Creeping Flow ◽  
Force Balance ◽  
Diffuse Interface ◽  
Critical Conditions

We investigate the critical conditions for the onset of motion of a three-dimensional droplet on a wall in shear flows at moderate Reynolds number. A diffuse-interface method is used for this purpose, which also circumvents the stress singularity at the moving contact line, and the method allows for a density and viscosity contrast between the fluids. Contact-angle hysteresis is represented by the prescription of a receding contact angle θRand an advancing contact angle value θA. Critical conditions are determined by tracking the motion and deformation of a droplet (initially a spherical cap with a uniform contact angle θ0). At sufficiently low values of a Weber number,We(based on the applied shear rate and the drop volume), the drop deforms and translates for some time, but subsequently reaches a stationary position and attains a steady-state shape. At sufficiently large values ofWeno such steady state is found. We present results for the critical value ofWeas a function of Reynolds numberRefor cases with the initial value of the contact angle θ0=θRas well as for θ0=θA. A scaling argument based on a force balance on the drop is shown to represent the results very accurately. Results are also presented for the static shape, transient motion and flow structure at criticality. It is shown that at lowReour results agree (with some qualifications) with those of Dimitrakopoulos & Higdon (1998,J. Fluid Mech. vol. 377, p. 189). Overall, the results indicate that the critical value ofWeis affected significantly by inertial effects at moderate Reynolds numbers, whereas the steady shape of droplets still shows some resemblance to that obtained previously for creeping flow conditions. The paper concludes with an investigation into the complex structure of a steady wake behind the droplet and the occurrence of a stagnation point at the upstream side of the droplet.

Transitions and instabilities of flow in a symmetric channel with a suddenly expanded and contracted part

2001 ◽  
Vol 434 ◽  
pp. 355-369 ◽  
Author(s):  
J. MIZUSHIMA ◽  
Y. SHIOTANI
Keyword(s):  
Reynolds Number ◽  
Steady State ◽  
Nonlinear Stability ◽  
Critical Reynolds Number ◽  
Reynolds Numbers ◽  
Critical Conditions ◽  
Stable Steady State ◽  
Weakly Nonlinear ◽  
Symmetric Channel ◽  
Steady State Solutions

Transitions and instabilities of two-dimensional flow in a symmetric channel with a suddenly expanded and contracted part are investigated numerically by three different methods, i.e. the time marching method for dynamical equations, the SOR iterative method and the finite-element method for steady-state equations. Linear and weakly nonlinear stability theories are applied to the flow. The transitions are confirmed experimentally by flow visualizations. It is known that the flow is steady and symmetric at low Reynolds numbers, becomes asymmetric at a critical Reynolds number, regains the symmetry at another critical Reynolds number and becomes oscillatory at very large Reynolds numbers. Multiple stable steady-state solutions are found in some cases, which lead to a hysteresis. The critical conditions for the existence of the multiple stable steady-state solutions are determined numerically and compared with the results of the linear and weakly nonlinear stability analyses. An exchange of modes for oscillatory instabilities is found to occur in the flow as the aspect ratio, the ratio of the length of the expanded part to its width, is varied, and its relation with the impinging free-shear-layer instability (IFLSI) is discussed.

Pulsatile Blood Flow in a Channel of Small Exponential Divergence—III. Unsteady Flow Separation

1977 ◽  
Vol 99 (2) ◽  
pp. 333-338 ◽  
Author(s):  
Daniel J. Schneck
Keyword(s):  
Reynolds Number ◽  
Blood Flow ◽  
Steady State ◽  
Flow Separation ◽  
Reynolds Numbers ◽  
Critical Value ◽  
Flow Oscillation ◽  
Pulsatile Blood Flow ◽  
Steady Streaming ◽  
Flow Through

Analysis of pulsatile flow through exponentially diverging channels reveals the existence of critical mean Reynolds numbers for which the flow separates at a downstream axial station. These Reynolds numbers vary directly with the frequency of flow oscillation and inversely with the rate of channel divergence. Increasing the Reynolds number above its critical value results in a rapid upstream displacement of the point of separation. For a tube of fixed geometry, periodic unsteadiness causes flow separation to occur at lower Reynolds numbers and upstream of a corresponding steady-state situation. The point of separation moves progressively downstream, however, towards its steady-state location, as the frequency of oscillation increases. These results are discussed as consequences of the nonlinear steady streaming phenomenon described in an earlier paper.

A numerical study of the motion of drops in Poiseuille flow. Part 1. Lateral migration of one drop

2000 ◽  
Vol 411 ◽  
pp. 325-350 ◽  
Author(s):  
SAEED MORTAZAVI ◽  
GRÉTAR TRYGGVASON
Keyword(s):  
Reynolds Number ◽  
Steady State ◽  
Poiseuille Flow ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Lateral Position ◽  
Physical Parameters ◽  
Computational Domain ◽  
Two Dimensional ◽  
Lateral Migration

The cross-stream migration of a deformable drop in two-dimensional Hagen–Poiseuille flow at finite Reynolds numbers is studied numerically. In the limit of a small Reynolds number (< 1), the motion of the drop depends strongly on the ratio of the viscosity of the drop fluid to the viscosity of the suspending fluid. For viscosity ratio 0.125 a drop moves toward the centre of the channel, while for ratio 1.0 it moves away from the centre until halted by wall repulsion. The rate of migration increases with the deformability of the drop. At higher Reynolds numbers (5–50), the drop either moves to an equilibrium lateral position about halfway between the centreline and the wall – according to the so-called Segre–Silberberg effect or it undergoes oscillatory motion. The steady-state position depends only weakly on the various physical parameters of the flow, but the length of the transient oscillations increases as the Reynolds number is raised, or the density of the drop is increased, or the viscosity of the drop is decreased. Once the Reynolds number is high enough, the oscillations appear to persist forever and no steady state is observed. The numerical results are in good agreement with experimental observations, especially for drops that reach a steady-state lateral position. Most of the simulations assume that the flow is two-dimensional. A few simulations of three-dimensional flows for a modest Reynolds number (Re = 10), and a small computational domain, confirm the behaviour seen in two dimensions. The equilibrium position of the three-dimensional drop is close to that predicted in the simulations of two-dimensional flow.

The dynamics and scaling law for particles suspended in shear flow with inertia

2000 ◽  
Vol 423 ◽  
pp. 317-344 ◽  
Author(s):  
E-JIANG DING ◽  
CYRUS K. AIDUN
Keyword(s):  
Reynolds Number ◽  
Steady State ◽  
Shear Flow ◽  
Scaling Law ◽  
Density Ratio ◽  
Reynolds Numbers ◽  
Creeping Flow ◽  
Elliptical Cylinder ◽  
Body Of Revolution ◽  
Time Periodic

The effect of inertia on the dynamics of a solid particle (a circular cylinder, an elliptical cylinder, and an ellipsoid) suspended in shear flow is studied by solving the discrete Boltzmann equation. At small Reynolds number, when inertia is negligible, the behaviour of the particle is in good agreement with the creeping flow solution showing periodic orbits. For an elliptical cylinder or an ellipsoid, the results show that by increasing the Reynolds number, the period of rotation increases, and eventually becomes infinitely large at a critical Reynolds number, Rec. At Reynolds numbers above Rec, the particle becomes stationary in a steady-state flow. It is found that the transition from a time-periodic to a steady state is through a saddle-node bifurcation, and, consequently, the period of oscillation near this transition is proportional to [mid ]p−pc[mid ]−1/2, where p is any parameter in the flow, such as the Reynolds number or the density ratio, which leads to this transition at p = pc. This universal scaling law is presented along with the physics of the transition and the effect of the inertia and the solid-to-fluid density ratio on the dynamics. It is conjectured that this transition and the scaling law are independent of the particle shape (excluding body of revolution) or the shear profile.

Flow instability of cylinder embedded within two-parallel moving walls of cuboidal cavity at different radii sizes

2020 ◽  
Vol 31 (05) ◽  
pp. 2050063
Author(s):  
Basma Souayeh
Keyword(s):  
Reynolds Number ◽  
Steady State ◽  
Flow Instability ◽  
Flow Behavior ◽  
Reynolds Numbers ◽  
Critical Value ◽  
The Other ◽  
Geometrical Parameters ◽  
Parallel Wall ◽  
Volume Method

A computational analysis has been performed to study the flow instability of two-parallel wall motions in a Cuboidal enclosure incorporated by a cylinder under different radii sizes. A numerical methodology based on the Finite Volume Method (FVM) and a full Multigrid acceleration is utilized in this paper. Left and right parallel walls of the cavity are maintained driven and all the other walls completing the domain are motionless. Different radii sizes ([Formula: see text], 0.1, 0.125, 0.15 and 0.175) are employed encompassing descriptive Reynolds numbers that range three orders of magnitude 100, 400 and 800 for the steady state. The obtained results show positions [Formula: see text] and [Formula: see text] of the inner cylinder promote cell distortion. Also, when the radius equates to [Formula: see text], it may lead to the birth of tertiary cells at [Formula: see text] which are more developed for [Formula: see text]. Thereafter, analysis of the flow evolution shows that with increasing Re beyond a certain critical value, the flow becomes unstable and undergoes a Hopf bifurcation. A nonuniform variation with the radius size of the inner cylinder is observed. Otherwise said, elongating the radius of the cylinder leads to decrease in the critical Reynolds number. Hence, the acceleration of the unsteadiness is realized. On the other hand, by further increasing Reynolds number more than the critical value from 1200 to 2100, we note that the kinetic energy is monotonously increasing with Reynolds number and a stronger motion in the velocity near the rear wall of the enclosure is observed. Furthermore, the symmetry of flow patterns observed in the steady state has been lost. Therefore, a systematic description of various effects illuminating the optimum geometrical parameters to achieve effective flow behavior in those systems has been successfully established through this paper.

Bifurcation Characteristics of Flows in Rectangular Sudden Expansion Channels

2005 ◽  
Author(s):  
Francine Battaglia ◽  
George Papadopoulos
Keyword(s):  
Reynolds Number ◽  
Laminar Flow ◽  
Critical Reynolds Number ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Expansion Ratio ◽  
Sudden Expansion ◽  
Two Dimensional ◽  
Effective Expansion ◽  
Three Dimensional Simulations

The effect of three-dimensionality on low Reynolds number flows past a symmetric sudden expansion in a channel was investigated. The geometric expansion ratio of in the current study was 2:1 and the aspect ratio was 6:1. Both experimental velocity measurements and two- and three-dimensional simulations for the flow along the centerplane of the rectangular duct are presented for Reynolds numbers in the range of 150 to 600. Comparison of the two-dimensional simulations with the experiments revealed that the simulations fail to capture completely the total expansion effect on the flow, which couples both geometric and hydrodynamic effects. To properly do so requires the definition of an effective expansion ratio, which is the ratio of the downstream and upstream hydraulic diameters and is therefore a function of both the expansion and aspect ratios. When the two-dimensional geometry was consistent with the effective expansion ratio, the new results agreed well with the three-dimensional simulations and the experiments. Furthermore, in the range of Reynolds numbers investigated, the laminar flow through the expansion underwent a symmetry-breaking bifurcation. The critical Reynolds number evaluated from the experiments and the simulations was compared to other values reported in the literature. Overall, side-wall proximity was found to enhance flow stability, helping to sustain laminar flow symmetry to higher Reynolds numbers in comparison to nominally two-dimensional double-expansion geometries. Lastly, and most importantly, when the logarithm of the critical Reynolds number from all these studies was plotted against the reciprocal of the effective expansion ratio, a linear trend emerged that uniquely captured the bifurcation dynamics of all symmetric double-sided planar expansions.

scholarly journals Nonlinear amplification in hydrodynamic turbulence

2021 ◽  
Vol 930 ◽  
Author(s):  
Kartik P. Iyer ◽  
Katepalli R. Sreenivasan ◽  
P.K. Yeung
Keyword(s):  
Reynolds Number ◽  
Isotropic Turbulence ◽  
Stokes Equations ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Navier Stokes ◽  
Navier Stokes Equations ◽  
Developed Turbulence ◽  
Advection Term ◽  
Nonlinear Amplification

Using direct numerical simulations performed on periodic cubes of various sizes, the largest being $8192^3$ , we examine the nonlinear advection term in the Navier–Stokes equations generating fully developed turbulence. We find significant dissipation even in flow regions where nonlinearity is locally absent. With increasing Reynolds number, the Navier–Stokes dynamics amplifies the nonlinearity in a global sense. This nonlinear amplification with increasing Reynolds number renders the vortex stretching mechanism more intermittent, with the global suppression of nonlinearity, reported previously, restricted to low Reynolds numbers. In regions where vortex stretching is absent, the angle and the ratio between the convective vorticity and solenoidal advection in three-dimensional isotropic turbulence are statistically similar to those in the two-dimensional case, despite the fundamental differences between them.

Time-Dependent Laminar Backward-Facing Step Flow in a Two-Dimensional Duct

1988 ◽  
Vol 110 (3) ◽  
pp. 289-296 ◽  
Author(s):  
F. Durst ◽  
J. C. F. Pereira
Keyword(s):  
Reynolds Number ◽  
Steady State ◽  
Separated Flow ◽  
Flow Region ◽  
Reynolds Numbers ◽  
Steady State Flow ◽  
Backward Facing Step ◽  
Step Flow ◽  
Recirculating Flow ◽  
State Flow

This paper presents results of numerical studies of the impulsively starting backward-facing step flow with the step being mounted in a plane, two-dimensional duct. Results are presented for Reynolds numbers of Re = 10; 368 and 648 and for the last two Reynolds numbers comparisons are given between experimental and numerical results obtained for the final steady state flow conditions. In the computational scheme, the convective terms in the momentum equations are approximated by a 13-point quadratic upstream weighted finite-difference scheme and a fully implicit first order forward differencing scheme is used to discretize the temporal derivatives. The computations show that for the higher Reynolds numbers, the flow starts to separate on the lower and upper corners of the step yielding two disconnected recirculating flow regions for some time after the flow has been impulsively started. As time progresses, these two separated flow regions connect up and a single recirculating flow region emerges. This separated flow region stays attached to the step, grows in size and approaches, for the time t → ∞, the dimensions measured and predicted for the separation region for steady laminar backward-facing flow. For the Reynolds number Re = 10 the separation starts at the bottom of the backward-facing step and the separation region enlarges with time until the steady state flow pattern is reached. At the channel wall opposite to the step and for Reynolds number Re = 368, a separated flow region is observed and it is shown to occur for some finite time period of the developing, impulsively started backward-facing step flow. Its dimensions change with time and reduce to zero before the steady state flow pattern is reached. For the higher Reynolds number Re = 648, the secondary separated flow region opposite to the wall is also present and it is shown to remain present for t → ∞. Two kinds of the inlet conditions were considered; the inlet mean flow was assumed to be constant in a first study and was assumed to increase with time in a second one. The predicted flow field for t → ∞ turned out to be identical for both cases. They were also identical to the flow field predicted for steady, backward-facing step flow using the same numerical grid as for the time-dependent predictions.

Experimental characterization of three-dimensional corner flows at low Reynolds numbers

2012 ◽  
Vol 707 ◽  
pp. 37-52 ◽  
Author(s):  
J. Sznitman ◽  
L. Guglielmini ◽  
D. Clifton ◽  
D. Scobee ◽  
H. A. Stone ◽  
...  
Keyword(s):  
Reynolds Number ◽  
Numerical Solutions ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Secondary Flows ◽  
Channel Cross Section ◽  
Experimental Characterization ◽  
Low Reynolds Numbers ◽  
Low Reynolds

AbstractWe investigate experimentally the characteristics of the flow field that develops at low Reynolds numbers ($\mathit{Re}\ll 1$) around a sharp $9{0}^{\ensuremath{\circ} } $ corner bounded by channel walls. Two-dimensional planar velocity fields are obtained using particle image velocimetry (PIV) conducted in a towing tank filled with a silicone oil of high viscosity. We find that, in the vicinity of the corner, the steady-state flow patterns bear the signature of a three-dimensional secondary flow, characterized by counter-rotating pairs of streamwise vortical structures and identified by the presence of non-vanishing transverse velocities (${u}_{z} $). These results are compared to numerical solutions of the incompressible flow as well as to predictions obtained, for a similar geometry, from an asymptotic expansion solution (Guglielmini et al., J. Fluid Mech., vol. 668, 2011, pp. 33–57). Furthermore, we discuss the influence of both Reynolds number and aspect ratio of the channel cross-section on the resulting secondary flows. This work represents, to the best of our knowledge, the first experimental characterization of the three-dimensional flow features arising in a pressure-driven flow near a corner at low Reynolds number.

Secondary instability and subcritical transition of the leading-edge boundary layer

2016 ◽  
Vol 792 ◽  
pp. 682-711 ◽  
Author(s):  
Michael O. John ◽  
Dominik Obrist ◽  
Leonhard Kleiser
Keyword(s):  
Boundary Layer ◽  
Reynolds Number ◽  
High Speed ◽  
Three Dimensional ◽  
Leading Edge ◽  
Reynolds Numbers ◽  
Secondary Instability ◽  
Bypass Transition ◽  
Wall Suction ◽  
Transition Mechanism

The leading-edge boundary layer (LEBL) in the front part of swept airplane wings is prone to three-dimensional subcritical instability, which may lead to bypass transition. The resulting increase of airplane drag and fuel consumption implies a negative environmental impact. In the present paper, we present a temporal biglobal secondary stability analysis (SSA) and direct numerical simulations (DNS) of this flow to investigate a subcritical transition mechanism. The LEBL is modelled by the swept Hiemenz boundary layer (SHBL), with and without wall suction. We introduce a pair of steady, counter-rotating, streamwise vortices next to the attachment line as a generic primary disturbance. This generates a high-speed streak, which evolves slowly in the streamwise direction. The SSA predicts that this flow is unstable to secondary, time-dependent perturbations. We report the upper branch of the secondary neutral curve and describe numerous eigenmodes located inside the shear layers surrounding the primary high-speed streak and the vortices. We find secondary flow instability at Reynolds numbers as low as$Re\approx 175$, i.e. far below the linear critical Reynolds number$Re_{crit}\approx 583$of the SHBL. This secondary modal instability is confirmed by our three-dimensional DNS. Furthermore, these simulations show that the modes may grow until nonlinear processes lead to breakdown to turbulent flow for Reynolds numbers above$Re_{tr}\approx 250$. The three-dimensional mode shapes, growth rates, and the frequency dependence of the secondary eigenmodes found by SSA and the DNS results are in close agreement with each other. The transition Reynolds number$Re_{tr}\approx 250$at zero suction and its increase with wall suction closely coincide with experimental and numerical results from the literature. We conclude that the secondary instability and the transition scenario presented in this paper may serve as a possible explanation for the well-known subcritical transition observed in the leading-edge boundary layer.

Sign in / Sign up

Export Citation Format

Share Document