The splash/non-splash boundary upon a dry surface and thin fluid film

On the basis of empirical data, power-law boundary relations are formulated to delineate the splash and non-splash regions on dry surfaces or thin films under isothermal conditions, using the Ohnesorge and Reynolds numbers. Approximation of the relations permits cancellation of fundamental fluid physical constants to give simplified formulas which provide insight into the governing parameters describing splashing and non-splashing behaviors. Thus, for a droplet impinging upon a dry solid surface, the splash/non-splash border is well described by √Ca = 0.35. For a drop impinging upon a thin fluid film, the analytical simplification yields a boundary described by √We = 20. For both expressions, values greater than the numerical value result in splashing.     

Steady and transient liquid sphere heating in a convective gas stream

 A finite-difference scheme has been developed to solve the equations governing the laminar forced convection heat transfer around and inside a spherical fluid droplet moving steadily in another immiscible fluid for both steady and transient thermal conditions. For large values of the external flow Reynolds number (Re), results not available in the literature have been obtained for circulating droplets at intermediate and high interior-to-exterior viscosity ratios (μ^*). Detailed results over a wide range of viscosity ratio (μ^*) and for 200≤Re≤1000 are presented for the temperature profiles outside and inside the sphere, Nusselt number, the time required to attain a uniform surface temperature and the time required to reach the steady-state temperature. Results show that convective heating is dependent on the external flow Reynolds number (Re) and the interior-to-exterior viscosity ratio (μ^*) where increasing Re or decreasing μ^* result in increasing heat transfer rate convected to the liquid sphere. Received on 1 March 1999     

The impact of a single drop on a wetted solid surface

 The impact of single drops on a thin liquid film was studied to understand the mechanism of secondary atomisation of sprays colliding on a wetted, cold, solid surface. To span a wide range of conditions various mixtures of water and glycerol were used. The use of Weber number, Ohnesorge number and non-dimensional film thickness to describe the peculiarities of the phenomenon allowed to carry out the experiments under appropriate similarity conditions. The impact of millimetric drops was analysed in detail by photographic means, using both still photography to study impact morphology, and laser sheet visualisation to investigate secondary droplet formation. Two mechanisms of splash were identified, depending essentially on the liquid viscosity (Ohnesorge number), a parameter which appears to play an important role also in defining the splash morphology. A photographic documentation is annexed. The characteristic times of the crown formation, the non-linear evolution of cusps (jet formation) and the surface roughness influence are further discussed. The experimental results allow to propose an empirical correlation for the splashing/deposition limit, for a wide range of conditions, and a comparison to available previous works is presented. The influence of the film thickness and liquid viscosity on the splash is confirmed and quantified. Received: 1 March 1996/Accepted: 12 November 1996     

Collision dynamics of high-speed droplets upon layers of variable thickness

The collision dynamics between a droplet and a film has been studied with high-impact energy that can be grouped in a dimensionless Weber number, We, as normalized by surface energy. To accomplish this, we have developed a technique based on cutting of a high-speed jet, which can generate a single droplet with speed up to 23 m/s and We on the order of thousands. It was found that the boundaries indicating the occurrence of a central jet and that of a secondary droplet disintegrated from the jet decreased monotonically with increased dimensionless film thickness, H, and remained constant when the film thickness was larger than the crater depth. However, the transition designating multiple droplets that are originated from a central jet shows a non-monotonic trend with the variation of H, with a minimum We being at H ≈ 3, which is about the maximum crater depth, owing to a tuning behavior. The critical We for splashing that occurs at an early phase immediately after the impact is relatively sensitive to the film thickness only when H is between 1 and 2, which increases with reduced H. At large We (≳2,570 for high H), the ejected crown is closed to form a bubble and the transition boundary reveals a similar dependence on H as that for creation of a central jet.     

Heat transfer in a power-law fluid film over a unsteady stretching sheet

In this paper the problem of momentum and heat transfer in a thin liquid film of power-law fluid on an unsteady stretching surface has been studied. Numerical solutions are obtained for some representative values of the unsteadiness parameter S and the power-law index n for a wide range of the generalized Prandtl number, 0.001 ≤ Pr ≤ 1000. Typical temperature and velocity profiles, the dimensionless film thickness, free-surface temperature, and the surface heat fluxes are presented at selected controlling parameters. The results show that increasing the value of n tends to increase the boundary-layer thickness and broadens the temperature distributions. The free-surface temperature of a shear thinning fluid is larger than that of a Newtonian fluid, but the opposite trend is true for a shear thickening fluid. For small generalized Prandtl numbers, the surface heat flux increases with a decrease in n, but the impacts of n on the heat transfer diminish for Pr greater than a moderate value (approximately 1 ≤ Pr ≤ 10, depending on the magnitude of S).     

Breakup of a droplet at high velocity impacting a solid surface

We have studied the collision between a droplet of different liquids with high impact energy and a solid plate with varied surface roughness, which is characterized by a dimensionless Weber number (We, defined as the impact inertia of the droplet normalized by its surface force) extending up to 12,000 for water. To make such collision, a technique was developed to generate a single droplet with speed up to 42 m/s, which was initially driven by upstream air flow through a nozzle and accelerated to nearly the same velocity of the high-speed flow downstream. Via a high-speed photographing system, the various splashing mechanisms were investigated and a specific prompt splash on a smooth plate was found at sufficiently high We, which was different somehow from the conventionally defined one that was generally believed to occur only on a rough surface. The radius when multiple secondary droplets were shed out of the rim of the expanding lamella was found to scale almost invariantly with We at large values, whereas the coupled effect of liquid viscosity might affect the ultimate value.     

Non-Darcian and Anisotropic Effects on Free Convection in a Porous Enclosure

The free convective flow and heat transfer, within the framework of Boussinesq approximation, in an anisotropic fluid filled porous rectangular enclosure subjected to end-to-end temperature difference have been investigated using Brinkman extended non-Darcy flow model. The studies involve simultaneous consideration of hydrodynamic and thermal anisotropy. The flow and temperature fields in general are governed by, Ra, the Rayleigh number, AR, the aspect ratio of the slab, K*, the permeability ratio and k*, the thermal conductivity ratio, and Da, Darcy number. Numerical solutions employing the successive accelerated replacement (SAR) scheme have been obtained for 100 ≤ Ra ≤ 1000, 0.5 ≤ AR ≤ 5, 0.5 ≤ K* ≤ 5, 0.5 ≤ k* ≤ 5, and 0 ≤ Da ≤ 0.1. It has been found that [`(Nu)]{\overline {Nu}}, average Nusselt number increases with increase in K* and decreases as k* increases. However, the magnitude of the change in [`(Nu)]{\overline {Nu}} depends on the parameter Da, characterizing the Brinkman extended non-Darcy flow.     

Vortex ring generation due to the coalescence of a water drop at a free surface

 Results are presented of an experimental investigation of vortex ring formation by a fluid drop contacting a free surface with negligible velocity. The pool fluid is mixed with fluorescein dye, and a laser sheet is used to illuminate a plane of the flow. A series of representative images is recorded by a CCD camera and speculation is made regarding specific sources of vorticity flux through the free surface. Two scaling analyses previously presented by other investigators are demonstrated to be equivalent under the assumptions of this experiment, and they provide the motivation for a series of test runs in which the duration of the coalescence process, τ_*, is related to variations in drop diameter L and fluid surface tension σ. Experimental results are in agreement with the analyses, showing τ_*∼σ^-1/2 and τ_*∼L ^3/2. Received: 22 December 1995 / Accepted: 15 October 1996     

The combined influence of a rough surface and thin fluid film upon the splashing threshold and splash dynamics of a droplet impacting onto them

Surface roughness can have a critical effect upon the splashing threshold and dynamics of a drop impacting on either a dry or rough solid surface or one coated by a thin fluid film. As most coating applications and spray systems quickly evolve to a state where the droplets impinge upon fluid deposited by preceding droplets, the combined contributions of surface roughness and a pre-deposited thin liquid film of comparable thickness upon droplet impingement dynamics are examined. For comparison, we include results for droplets impacting on a smooth, dry surface and a smooth surface wetted by a thin fluid film. The inclusion of surface roughness considerably lowers the splashing threshold and alters the splashing dynamics such that differences in fluid surface tensions between 20.1 and 72.8 dynes/cm or viscosities between 0.4 and 3.3 cP have little effect.     

Transient combined natural convection-conduction in open-ended vertical concentric annuli

Transient conjugate natural convection heat transfer in open-ended vertical concentric annuli is investigated numerically. The governing equations of an induced laminar flow for a fluid of Pr=0.7 are solved using a finite-difference technique. The heating is achieved by a step change in the temperature of the outer surface of the outer tube while the inner surface of the inner tube is kept adiabatic. The range of Grashof number considered is 500≤Gr^*≤10⁵. The effects of solid-fluid conductivity ratio and diffusivity ratio on the transient induced flow characteristics are presented. Received on 25 March 1998     

Droplet deformation of a strongly shear thinning dense suspension of polymeric micelles

We investigated the deformation of a strong shear thinning droplet undergoing simple shear flow in a Newtonian liquid. The droplet was an aqueous solution of poly(ethylene oxide) end capped with an alkyl group that forms spherical micelles in aqueous solution. At high concentrations and below a critical temperature, the jammed micelles showed strong shear thinning behaviour, and neither a yield stress nor a Newtonian viscosity was observed. At small shear rates, the droplet rotated and aligned in the flow, but did not deform or only very weakly. At high shear rates, the droplet deformation increased with increasing shear rate. The deformed droplet did not relax after the shear was stopped except for a modest rounding of the edges. For each shear rate, an apparent viscosity, η _ad, of the equivalent Newtonian droplet was calculated assuming affine deformation. η _ad showed a power law dependence on the capillary number Ca with an exponent of − 1.8 and was larger than the shear viscosity of the micelle suspension at the same shear rates. The results were explained by the existence of a strong gradient of the viscosity inside the droplet leading to a very low viscosity fluid layer near the droplet/matrix interface.     

Study of airborne particles produced by normal impact of millimetric droplets onto a liquid film

The aim of this work is to carry out an experimental investigation into the generation of airborne microparticles when millimetric droplets of aqueous solutions impact onto a liquid film. Impact experiments using 3.9 mm diameter droplets were carried out for Weber numbers between 159 and 808, with a fixed Ohnesorge number of 2 × 10⁻³ and film parameters S _f (the ratio between the thickness of the liquid film h _film and the diameter of the impacting droplet d _i) between 0.3 and 1. Observed results show that the deposition/splashing threshold is independent of the parameter S _f in agreement with the data in the literature. The aerosol measurement results demonstrate the production of solid particles from the evaporation of secondary microdroplets with diameters less than 30 μm formed when splash occurs. The median diameter of these microdroplets is around 20 μm, corresponding to a value of d ₅₀/d _i = 5 × 10⁻³. Taken together, the results show that the mass and the number of particles emitted increase as the Weber number increases. Moreover, at a Weber number of 808, the results show that the mass and number of particles emitted increases as the parameter S _f decreases. In this case, the mean number of microdroplets emitted per impact is equal to 14 for S _f = 1 and equal to 76 for S _f = 0.3.     

Non-Darcy buoyancy driven flows in a fluid saturated porous medium: the use of asymptotic computational fluid dynamics (ACFD) approach

Natural convection in a fluid saturated porous medium has been numerically investigated using a generalized non-Darcy approach. The governing equations are solved by using Finite Volume approach. First order upwind scheme is employed for convective formulation and SIMPLE algorithm for pressure velocity coupling. Numerical results are presented to study the influence of parameters such as Rayleigh number (10⁶ ≤Ra ≤10⁸), Darcy number (10⁻⁵ ≤ Da ≤ 10⁻²), porosity (0.4 ≤ ɛ ≤ 0.9) and Prandtl number (0.01 ≤ Pr ≤ 10) on the flow behavior and heat transfer. By combining the method of matched asymptotic expansions with computational fluid dynamics (CFD), so called asymptotic computational fluid dynamics (ACFD) technique has been employed to generate correlation for average Nusselt number. The technique is found to be an attractive option for generating correlation and also in the analysis of natural convection in porous medium over a fairly wide range of parameters with fewer simulations for numerical solutions.     

A new algorithm for surface tension model in moving particle methods

A new algorithm for the surface tension model was developed for moving particle methods. The algorithm is based on the link‐list search algorithm and the continuum surface tension (CST) model. The developed algorithm with the CST model was implemented to a kind of moving particle approach, the finite volume particle (FVP) method. The FVP method with the new algorithm was tested by oscillatory behaviour of a two‐dimensional droplet. The oscillatory period agrees well with analytical one, and the transient shape of the droplet is also in good agreement with that obtained by other numerical methods. The droplet impact on a liquid surface was also studied using the new algorithm. The deposition and splashing phenomena were clearly reproduced. Simulated spread radius of the splashing phenomena was consistent with a power law. Copyright © 2007 John Wiley & Sons, Ltd.     

Droplet formation in the rupture of a liquid film under the action of a thermal load

The deformation and rupture of a liquid film, suspended between two solid walls, under the action of a localized thermal load is considered. For studying this process, a two-dimensional model is used, which describes the motion of a thin layer of a viscous non-isothermal liquid under microgravity conditions. For modeling the dynamics of the liquid, the Navier-Stokes equations in the “vorticity—stream function” variables are used. A numerical analysis of the influence of thermal loads on the deformation and the mechanism of the rupture of freely suspended films is performed. Is is shown that for a certain width of the thermal beam acting on the film free surface the rupture of the film may occur with the formation of a droplet. The results of the solution of model problems are presented.     

Well-posedness and Stability of a Multi-dimensional Tumor Growth Model

We study a moving boundary problem modeling the growth of in vitro tumors. This problem consists of two elliptic equations describing the distribution of the nutrient and the internal pressure, respectively, and a first-order partial differential equation describing the evolution of the moving boundary. An important feature is that the effect of surface tension on the moving boundary is taken into account. We show that this problem is locally well-posed for a large class of initial data by using analytic semi-group theory. We also prove that if the surface tension coefficient γ is larger than a threshold value γ _* then the unique flat equilibrium is asymptotically stable, whereas in the case γ  < γ _* this flat equilibrium is unstable.     

Effects of solid and slotted baffles on the convection characteristics of backward-facing step flow in a channel

This study investigates the enhancement of the laminar forced convection characteristics of backward-facing step flow in a two-dimensional channel through the installation of solid and slotted baffles onto the channel wall. The effects of the height of baffle H _b, inclination of baffle installation ϕ_b, height of slot in baffle H _t, inclination of slot in baffle ϕ_t, and distance between the backward-facing step and baffle D on the flow structure, temperature distribution and Nusselt number variation for the system at various Re are numerically explored. Results show that a slotted baffle can enhance the average Nusselt number for the heating section of channel plate by the maximum 190% when Pr=0.7, H _s=0.5, L=5, H _b ≤ 0.3, W _b ≤ 0.2, 0.1 ≤ D ≤ 0.5, 0° ≤ ϕ_b ≤ 45°, H _t ≤ 0.1, 0° ≤ ϕ_t ≤ 45° and 50 ≤ Re ≤ 400. As for the solid baffle, the enhancement may be up by 230%. The solid baffle might cause the re-separation of main stream, and consequently result in poor local heat transfer coefficient in the end region of heating section. This disadvantage can be obviously improved as the baffle is slotted. Besides the penalty of increase in pressure drop due to the baffle installation is much higher for the situation with solid baffle.     

An experimental investigation of negative wakes behind spheres settling in a shear-thinning viscoelastic fluid

We present detailed experimental results examining “negative wakes” behind spheres settling along the centerline of a tube containing a viscoelastic aqueous polyacrylamide solution. Negative wakes are found for all Deborah numbers (2.43≤De(˙γ)≤8.75) and sphere-to-tube aspect ratios (0.060≤a/R≤0.396) examined. The wake structures are investigated using laser-Doppler velocimetry (LDV) to examine the centerline fluid velocity around the sphere and digital particle image velocimetry (DPIV) for full-field velocity profiles. For a fixed aspect ratio, the magnitude of the most negative velocity, U _min , in the wake is seen to increase with increasing De. Additionally, as the Deborah number becomes larger, the location of this minimum velocity shifts farther downstream. When normalized with the sphere radius and the steady state velocity of the sphere, the axial velocity profiles become self-similar to the point of the minimum velocity. Beyond this point, the wake structure varies weakly with aspect ratio and De, and it extends more than 20 radii downstream. Inertial effects at high Reynolds numbers are observed to shift the entire negative wake farther downstream. Using DPIV to investigate the transient kinematic response of the fluid to the initial acceleration of the sphere from rest, it is seen that the wake develops from the nonlinear fluid response at large strains. Measurements of the transient uniaxial extensional viscosity of this weakly strain-hardening fluid using a filament stretching rheometer show that the existence of a negative wake is consistent with theoretical arguments based on the opposing roles of extensional stresses and shearing stresses in the wake of the sphere. Received: 10 November 1997 Accepted: 1 May 1998     

Numerical investigation of natural convection of Al<Subscript>2</Subscript>O<Subscript>3</Subscript> nanofluid in vertical annuli

This paper presents the numerical study of internal free convection of Al₂O₃ water nanofluid in vertical annuli. Vertical walls are maintained at constant temperatures and horizontal walls are adiabatic. Results are validated by experimental data. Effect of nanofluids on natural convection is investigated as a function of geometrical and physical parameters and particle fractions for aspect ratio of 1 ≤ H/L ≤ 5, Grashof number of 10³ ≤ Gr ≤ 10⁵ and concentration of 0 ≤ ϕ ≤ 0.06. More than 330 different numerical cases are investigated to develop a new correlation for the Nusselt number. This correlation is presented as a function of Nusselt number of base fluid and particle fraction which is a linear decreasing function of particle fraction. The developed correlation for annuli is also valid for the natural convection of Al₂O₃ water nanofluid in a square cavity. Furthermore, the effect of the viscosity and conductivity models on the Nusselt number of nanofluids in cylindrical cavities are discussed.     

Shear and normal stress relaxation in polymer blends: stress predictions from interface velocity and Laplace pressure terms

We derived for the first time the relationships among shear stress and normal stress differences for ellipsoidal interfaces under large step shear strains considering interface velocity term and Laplace pressure term in the expression of the stress tensor for mixtures of two Newtonian fluids. In the derivation, orientation angle of the interface is assumed to be given by the affine deformation assumption and is independent of time based on experimental results for blends with 0.048 ≤ K ≤ 0.54 where K is the ratio of droplet viscosity to matrix viscosity. For ellipsoidal droplets, the shear stress is only proportional to the first normal stress difference. On the other hand, for spheroidal droplets, proportionality among the shear stress, the first and the second normal stress differences was derived, and the ratio of the second normal stress difference to the first normal stress difference was given as a function of step strain. The shear stress and the first normal stress difference obtained experimentally satisfy the derived relationship, indicating applicability of the stress expression for polymer blends.