Numerical solutions for thin film flow down the outside and inside of a vertical cylinder

We consider a model for thin film flow down the outside and inside of a vertical cylinder. Our focus is to study the effect that the curvature of the cylinder has on the gravity-driven instability of the advancing contact line and to simulate the resulting fingering patterns that form due to this instability. The governing partial differential equation is fourth order with a nonlinear degenerate diffusion term that represents the stabilising effect of surface tension. We present numerical solutions obtained by implementing an efficient alternating direction implicit scheme. When compared to the problem of flow down a vertical plane, we find that increasing substrate curvature tends to increase the fingering instability for flow down the outside of the cylinder, whereas flow down the inside of the cylinder substrate curvature has the opposite effect. Further, we demonstrate the existence of nontrivial travelling wave solutions which describe fingering patterns that propagate down the inside of a cylinder at constant speed without changing form. These solutions are perfectly analogous to those found previously for thin film flow down an inclined plane.

The hard-particle model for a dense granular flow down an inclined plane

Perfectly hard particles are those which experience an infinite repulsive force when they overlap, and no force when they do not overlap. In the hard-particle model, the only static state is the isostatic state where the forces between particles are statically determinate. In the flowing state, the interactions between particles are instantaneous because the time of contact approaches zero in the limit of infinite particle stiffness. Here, we discuss the development of a hard particle model for a realistic granular flow down an inclined plane, and examine its utility for predicting the salient features both qualitatively and quantitatively. We first discuss Discrete Element simulations, that even very dense flows of sand or glass beads with volume fraction between 0.5 and 0.58 are in the rapid flow regime, due to the very high particle stiffness. An important length scale in the shear flow of inelastic particles is the `conduction length' delta = (d/(1 - e(2))(1/2)), where d is the particle diameter and e is the coefficient of restitution. When the macroscopic scale h (height of the flowing layer) is larger than the conduction length, the rates of shear production and inelastic dissipation are nearly equal in the bulk of the flow, while the rate of conduction of energy is O((delta/h)(2)) smaller than the rate of dissipation of energy. Energy conduction is important in boundary layers of thickness delta at the top and bottom. The flow in the boundary layer at the top and bottom is examined using asymptotic analysis. We derive an exact relationship showing that the a boundary layer solution exists only if the volume fraction in the bulk decreases as the angle of inclination is increased. In the opposite case, where the volume fraction increases as the angle of inclination is increased, there is no boundary layer solution. The boundary layer theory also provides us with a way of understanding the cessation of flow when at a given angle of inclination when the height of the layer is decreased below a value h(stop), which is a function of the angle of inclination. There is dissipation of energy due to particle collisions in the flow as well as due to particle collisions with the base, and the fraction of energy dissipation in the base increases as the thickness decreases. When the shear production in the flow cannot compensate for the additional energy drawn out of the flow due to the wall collisions, the temperature decreases to zero and the flow stops. Scaling relations can be derived for h(stop) as a function of angle of inclination.

Kinetic theory for the density plateau in the granular flow down an inclined plane

The striking lack of observable variation of the volume fraction with height in the center of a granular flow down an inclined plane is analysed using constitutive relations obtained from kinetic theory. It is shown that the rate of conduction in the granular energy balance equation is O(delta(2)) smaller than the rate of production of energy due to mean shear and the rate of dissipation due to inelastic collisions, where the small parameter delta = (d/(1 - e(n))H-1/2), d is the particle diameter, en is the normal coefficient of restitution and H is the thickness of the flowing layer. This implies that the volume fraction is a constant in the leading approximation in an asymptotic analysis in small delta. Numerical estimates of both the parameter delta and its pre-factor are obtained to show that the lack of observable variation of the volume fraction with height can be explained by constitutive relations obtained from kinetic theory.

Dense granular flow down an inclined plane A comparison between the hard particle model and soft particle simulations

The granular flow down an inclined plane is simulated using the discrete element (DE) technique to examine the extent to which the dynamics of an unconfined dense granular flow can be well described by a hard particle model First, we examine the average coordination number for the particles in the flow down an inclined plane using the DE technique using the linear contact model with and without friction, and the Hertzian contact model with friction The simulations show that the average coordination number decreases below 1 for values of the spring stiffness corresponding to real materials, such as sand and glass, even when the angle of inclination is only 10 larger than the angle of repose Additional measures of correlations in the system, such as the fraction of particles with multibody contact, the force ratio (average ratio of the magnitudes of the largest and the second largest force on a particle), and the angle between the two largest forces on the particle, show no evidence of force chains or other correlated motions in the system An analysis of the bond-orientational order parameter indicates that the flow is in the random state, as in event-driven (ED) simulations V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] The results of the two simulation techniques for the Bagnold coefficients (ratio of stress and square of the strain rate) and the granular temperature (mean square of the fluctuating velocity) are compared with the theory V Kumaran, J Fluid Mech 632, 107 (2009), J Fluid Mech 632, 145 (2009)] and are found to be in quantitative agreement In addition, we also conduct a comparison of the collision frequency and the distribution of the precollisional relative velocities of particles in contact The strong correlation effects exhibited by these two quantities in event-driven simulations V Kumaran, J Fluid Mech 632, 145 (2009)] are also found in the DE simulations (C) 2010 American Institute of Physics doi 10 1063/1 3504660]

Effect of base dissipation on the granular flow down an inclined plane

The effect of base dissipation on the granular flow down an inclined plane is examined by altering the coefficient of restitution between the moving and base particles in discrete element (DE) simulations. The interaction laws between two moving particles are kept fixed, and the coefficient of restitution (damping constant in the DE simulations) between the base and moving particles are altered to reduce dissipation, and inject energy from the base. The energy injection does result in an increase in the strain rate by up to an order of magnitude, and the temperature by up to two orders of magnitude at the base. However, the volume fraction, strain rate and temperature profiles in the bulk (above about 15 particle diameters from the base) are altered very little by the energy injection at the base. We also examine the variation of h(stop), the minimum height at the cessation of flow, with energy injection from the base. It is found that at a fixed angle of inclination, h(stop) decreases as the energy dissipation at the base decreases.

Transition due to base roughness in a dense granular flow down an inclined plane

Particle simulations based on the discrete element method are used to examine the effect of base roughness on the granular flow down an inclined plane. The base is composed of a random configuration of fixed particles, and the base roughness is decreased by decreasing the ratio of diameters of the base and moving particles. A discontinuous transition from a disordered to an ordered flow state is observed when the ratio of diameters of base and moving particles is decreased below a critical value. The ordered flowing state consists of hexagonally close packed layers of particles sliding over each other. The ordered state is denser (higher volume fraction) and has a lower coordination number than the disordered state, and there are discontinuous changes in both the volume fraction and the coordination number at transition. The Bagnold law, which states that the stress is proportional to the square of the strain rate, is valid in both states. However, the Bagnold coefficients in the ordered flowing state are lower, by more than two orders of magnitude, in comparison to those of the disordered state. The critical ratio of base and moving particle diameters is independent of the angle of inclination, and varies very little when the height of the flowing layer is doubled from about 35 to about 70 particle diameters. While flow in the disordered state ceases when the angle of inclination decreases below 20 degrees, there is flow in the ordered state at lower angles of inclination upto 14 degrees. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4710543]

Stability of developing film flow down an inclined surface

Film flows on inclined surfaces are often assumed to be of constant thickness, which ensures that the velocity profile is half-Poiseuille. It is shown here that by shallow water theory, only flows in a portion of Reynolds number-Froude number (Re-Fr) plane can asymptotically attain constant film thickness. In another portion on the plane, the constant thickness solution appears as an unstable fixed point, while in other regions the film thickness seems to asymptote to a positive slope. Our simulations of the Navier-Stokes equations confirm the predictions of shallow water theory at higher Froude numbers, but disagree with them at lower Froude numbers. We show that different regimes of film flow show completely different stability behaviour from that predicted earlier. Supercritical decelerating flows are shown to be always unstable, whereas accelerating flows become unstable below a certain Reynolds number for a given Froude number. Subcritical flows on the other hand are shown to be unstable above a certain Reynolds number. In some range of parameters, two solutions for the base flowexist, and the attached profile is found to be more stable. All flows except those with separation become more stable as they proceed downstream. (C) 2012 American Institute of Physics. http://dx.doi.org/10.1063/1.4758299]

The effect of base roughness on the development of a dense granular flow down an inclined plane

The development of the flow of a granular material down an inclined plane starting from rest is studied as a function of the base roughness. In the simulations, the particles are rough frictional spheres interacting via the Hertz contact law. The rough base is made of a random configuration of fixed spheres with diameter different from the flowing particles, and the base roughness is decreased by decreasing the diameter of the base particles. The transition from an ordered to a disordered flowing state at a critical value of the base particle diameter, first reported by Kumaran and Maheshwari Phys. Fluids 24, 053302 (2012)] for particles with the linear contact model, is observed for the Hertzian contact model as well. The flow development for the ordered and disordered flows is very different. During the development of the disordered flow for the rougher base, there is shearing throughout the height. During the development of the ordered flow for the smoother base, there is a shear layer at the bottom and a plug region with no internal shearing above. In the shear layer, the particles are layered and hexagonally ordered in the plane parallel to the base, and the velocity profile is well approximated by Bagnold law. The flow develops in two phases. In the first phase, the thickness of the shear layer and the maximum velocity increase linearly in time till the shear front reaches the top. In the second phase, after the shear layer encompasses the entire flow, there is a much slower increase in the maximum velocity until the steady state is reached. (C) 2013 AIP Publishing LLC.

Effect of surfactants on the instability of a two-layer film flow down an inclined plane

The effect of insoluble surfactants on the instability of a two-layer film flow down an inclined plane is investigated based on the Orr-Sommerfeld boundary value problem. The study, focusing on Stokes flow P. Gao and X.-Y. Lu, ``Effect of surfactants on the inertialess instability of a two-layer film flow,'' J. Fluid Mech. 591, 495-507 (2007)], is further extended by including the inertial effect. The surface mode is recognized along with the interface mode. The initial growth rate corresponding to the interface mode accelerates at sufficiently long-wave regime in the presence of surface surfactant. However, the maximum growth rate corresponding to both interface and surface modes decelerates in the presence of surface surfactant when the upper layer is more viscous than the lower layer. On the other hand, when the upper layer is less viscous than the lower layer, a new interfacial instability develops due to the inertial effect and becomes weaker in the presence of interfacial surfactant. In the limit of negligible surface and interfacial tensions, respectively, two successive peaks of temporal growth rate appear in the long-wave and short-wave regimes when the interface mode is analyzed. However, in the case of the surface mode, only the long-wave peak appears. (C) 2014 AIP Publishing LLC.

Linear stability of miscible two-fluid flow down an incline

We show that miscible two-layer free-surface flows of varying viscosity down an inclined substrate are different in their stability characteristics from both immiscible two-layer flows, and flows with viscosity gradients spanning the entire flow. New instability modes arise when the critical layer of the viscosity transport equation overlaps the viscosity gradient. A lubricating configuration with a less viscous wall layer is identified to be the most stabilizing at moderate miscibility (moderate Peclet numbers). This also is in contrast with the immiscible case, where the lubrication configuration is always destabilizing. The co-existence that we find under certain circumstances, of several growing overlap modes, the usual surface mode, and a Tollmien-Schlichting mode, presents interesting new possibilities for nonlinear breakdown. © 2013 AIP Publishing LLC.

Gravity-driven fingering simulations for a thin liquid film flowing down the outside of a vertical cylinder

A numerical study is presented to examine the fingering instability of a gravity-driven thin liquid film flowing down the outer wall of a vertical cylinder. The lubrication approximation is employed to derive an evolution equation for the height of the film, which is dependent on a single parameter, the dimensionless cylinder radius. This equation is identified as a special case of that which describes thin film flow down an inclined plane. Fully three-dimensional simulations of the film depict a fingering pattern at the advancing contact line. We find the number of fingers observed in our simulations to be in excellent agreement with experimental observations and a linear stability analysis reported recently by Smolka & SeGall (Phys Fluids 23, 092103 (2011)). As the radius of the cylinder decreases, the modes of perturbation have an increased growth rate, thus increasing cylinder curvature partially acts to encourage the contact line instability. In direct competition with this behaviour, a decrease in cylinder radius means that fewer fingers are able to form around the circumference of the cylinder. Indeed, for a sufficiently small radius, a transition is observed, at which point the contact line is stable to transverse perturbations of all wavenumbers. In this regime, free surface instabilities lead to the development of wave patterns in the axial direction, and the flow features become perfectly analogous to the two-dimensional flow of a thin film down an inverted plane as studied by Lin & Kondic (Phys Fluids 22, 052105 (2010)). Finally, we simulate the flow of a single drop down the outside of the cylinder. Our results show that for drops with low volume, the cylinder curvature has the effect of increasing drop speed and hence promoting the phenomenon of pearling. On the other hand, drops with much larger volume evolve to form single long rivulets with a similar shape to a finger formed in the aforementioned simulations.

Vortex motion in the early stages of unsteady flow around a circular cylinder

In this paper, processes in the early stages of vortex motion and the development of flow structure behind an impulsively-started circular cylinder at high Reynolds number are investigated by combining the discrete vortex model with boundary layer theory, considering the separation of incoming flow boundary layer and rear shear layer in the recirculating flow region. The development of flow structure and vortex motion, particularly the formation and development of secondary vortex and a pair of secondary vortices and their effect on the flow field are calculated. The results clearly show that the flow structure and vortices motion went through a series of complicated processes before the symmetric main vortices change into asymmetric: development of main vortices induces secondary vortices; growth of the secondary vortices causes the main vortex sheets to break off and causes the symmetric main vortices to become “free” vortices, while a pair of secondary vortices is formed; then the vortex sheets, after breaking off, gradually extend downstream and the structure of a pair of secondary vortices becomes relaxed. These features of vortex motion look very much like the observed features in some available flow field visualizations. The action of the secondary vortices causes the main vortex sheets to break off and converts the main vortices into free vortices. This should be the immediate cause leading to the instability of the motion of the symmetric main vortices. The flow field structure such as the separation position of boundary layer and rear shear layer, the unsteady pressure distributions and the drag coefficient are calculated. Comparison with other results or experiments is also made.

Simulation of gas flow over a micro cylinder

Gas flow over a micro cylinder is simulated using both a compressible Navier-Stokes solver and a hybrid continuum /particle approach. The micro cylinder flow has low Reynolds number because of the small length scale and the low speed, which also indicates that the rarefied gas effect exists in the flow. A cylinder having a diameter of 20 microns is simulated under several flow conditions where the Reynolds number ranges from 2 to 50 and the Mach number varies from 0.1 to 0.8. It is found that the low Reynolds number flow can be compressible even when the Mach number is less than 0.3, and the drag coefficient of the cylinder increases when the Reynolds number decreases. The compressible effect will increase the pressure drag coefficient although the friction coefficient remains nearly unchanged. The rarefied gas effect will reduce both the friction and pressure drag coefficients, and the vortex in the flow may be shrunk or even disappear.

An efficient shock capturing scheme for the pseudo-unsteady equations representing steady inviscid flow

An efficient algorithm is presented for the solution of the steady Euler equations of gas dynamics. The scheme is based on solving linearised Riemann problems approximately and in more than one dimension incorporates operator splitting. The scheme is applied to a standard test problem of flow down a channel containing a circular arc bump for three different mesh sizes.

Dynamics of dense sheared granular flows. Part II. The relative velocity distributions

The distribution of relative velocities between colliding particles in shear flows of inelastic spheres is analysed in the Volume fraction range 0.4-0.64. Particle interactions are considered to be due to instantaneous binary collisions, and the collision model has a normal coefficient of restitution e(n) (negative of the ratio of the post- and pre-collisional relative velocities of the particles along the line joining the centres) and a tangential coefficient of restitution e(t) (negative of the ratio of post- and pre-collisional velocities perpendicular to line joining the centres). The distribution or pre-collisional normal relative velocities (along the line Joining the centres of the particles) is Found to be an exponential distribution for particles with low normal coefficient of restitution in the range 0.6-0.7. This is in contrast to the Gaussian distribution for the normal relative velocity in all elastic fluid in the absence of shear. A composite distribution function, which consists of an exponential and a Gaussian component, is proposed to span the range of inelasticities considered here. In the case of roughd particles, the relative velocity tangential to the surfaces at contact is also evaluated, and it is found to be close to a Gaussian distribution even for highly inelastic particles.Empirical relations are formulated for the relative velocity distribution. These are used to calculate the collisional contributions to the pressure, shear stress and the energy dissipation rate in a shear flow. The results of the calculation were round to be in quantitative agreement with simulation results, even for low coefficients of restitution for which the predictions obtained using the Enskog approximation are in error by an order of magnitude. The results are also applied to the flow down an inclined plane, to predict the angle of repose and the variation of the volume fraction with angle of inclination. These results are also found to be in quantitative agreement with previous simulations.