Heating analysis of a droplet on stretchable hydrophilic surface

Heating of a droplet on a stretchable hydrophilic surface is investigated and fluid dynamics in the droplet under the heating load is assessed. Elastomer wafers are considered as the sample material and the fixture is designed and manufactured to assure uniform stretching of the droplet located elastomer surface. Droplet adhesion and possible slipping/sliding of the droplet are evaluated during stretching of the sample surface. Numerical simulations are carried out to predict thermal and flow response of the droplet fluid before and after stretching. The effect of droplet volume on heating enhancement is also included in the numerical simulations. Experiments are carried out using a high-speed recording system towards comparing the flow predictions. Findings reveal that predictions are in agreement with their counterparts of experiments. Stretching of sample surface increases wetting area and lowers height of the droplet while influencing thermal flow structures in the fluid. The Nusselt and the Bond numbers increase with enlarging stretching, which becomes more visible for large droplet volume (80 µl). Hence, stretching corresponding to 80% extension of elastomer surface gives rise to 60% improvement in the Nusselt number.     

Dynamic behaviors of droplet impact and spreading: A universal relationship study of dimensionless wetting diameter and droplet height

A notable universal relationship has been proposed in the literature for the evolution of dimensionless droplet height and wetting diameter during the initial spreading stage of droplet impingement. In this study, this universal relationship was investigated by employing three sets of measurements. Sequential images were recorded, and the whole droplet profile ensembles were plotted to facilitate this study. These sets of experiments were designed by changing impact velocity, surface hydrophobicity, or solution property. The experimental results illustrate that the importance of parameters causing the data variation is in the order of surface hydrophobicity > initial impact velocity > surfactant on wetting diameter, and surface hydrophobicity ≈ initial impact velocity > surfactant on droplet height. No universal relationship was observed for dimensionless droplet height and wetting diameter.     

Droplet size and velocity characteristics of water-air impinging jet atomizer

A water-air impinging jets atomizer is investigated in this study, which consists of flow visualization using high speed photography and mean droplet size and velocity distribution measurements of the spray using Phase Doppler Anemometry (PDA). Topological structures and break up details of the generated spray in the far and near fields are presented with and without air jet and for an impinging angle of 90°. Spray angle increases with the water jet velocity, air flow rate and impinging angle. PDA results indicate that droplet size is smallest in the spray center, with minimum value of Sauter mean diameter (SMD) of 50 µm at the air flow rate of Q_m = 13.50 g/min. SMD of droplets increases towards the spray outer region gradually to about 120 µm. The mean droplet velocity component W along the air-jet axis is highest in the spray center and decreases gradually with increasing distance from the spray center. SMD normalized by the air nozzle diameter is found firstly to decrease with gas-to-liquid mass ratio (GLR) and air-to-liquid momentum ratio (ALMR) and then remain almost constant. Its increasing with aerodynamic Weber number indicates an exponential variation. The study sheds light on the performance of water-air impinging jets atomizers providing useful information for future CFD simulation works.     

Mixed convection boundary layers in the stagnation-point flow toward a stretching vertical sheet

An analysis is made for the steady mixed convection boundary layer flow near the two-dimensional stagnation-point flow of an incompressible viscous fluid over a stretching vertical sheet in its own plane. The stretching velocity and the surface temperature are assumed to vary linearly with the distance from the stagnation-point. Two equal and opposite forces are impulsively applied along the x-axis so that the wall is stretched, keeping the origin fixed in a viscous fluid of constant ambient temperature. The transformed ordinary differential equations are solved numerically for some values of the parameters involved using a very efficient numerical scheme known as the Keller-box method. The features of the flow and heat transfer characteristics are analyzed and discussed in detail. Both cases of assisting and opposing flows are considered. It is observed that, for assisting flow, both the skin friction coefficient and the local Nusselt number increase as the buoyancy parameter increases, while only the local Nusselt number increases but the skin friction coefficient decreases as the Prandtl number increases. For opposing flow, both the skin friction coefficient and the local Nusselt number decrease as the buoyancy parameter increases, but both increase as Pr increases. Comparison with known results is excellent.     

Experimental investigation of inclined liquid water jet flow onto vertically located superhydrophobic surfaces

In this study, the behaviour of an inclined water jet, which is impinged onto hydrophobic and superhydrophobic surfaces, has been investigated experimentally. Water jet was impinged with different inclination angles (15°–45°) onto five different hydrophobic surfaces made of rough polymer, which were held vertically. The water contact angles on these surfaces were measured as 102°, 112°, 123°, 145° and 167° showing that the last surface was superhydrophobic. Two different nozzles with 1.75 and 4 mm in diameters were used to create the water jet. Water jet velocity was within the range of 0.5–5 m/s, thus the Weber number varied from 5 to 650 and Reynolds number from 500 to 8,000 during the experiments. Hydrophobic surfaces reflected the liquid jet depending on the surface contact angle, jet inclination angle and the Weber number. The variation of the reflection angle with the Weber number showed a maximum value for a constant jet angle. The maximum value of the reflection angle was nearly equal to half of the jet angle. It was determined that the viscous drag decreases as the contact angle of the hydrophobic surface increases. The drag force on the wall is reduced dramatically with superhydrophobic surfaces. The amount of reduction of the average shear stress on the wall was about 40%, when the contact angle of the surface was increased from 145° to 167°. The area of the spreading water layer decreased as the contact angle of the surface increased and as the jet inclination angle, Weber number and Reynolds number decreased.     

Thermal radiation effects on the flow and heat transfer in a liquid film on an unsteady stretching sheet

The influence of thermal radiation on the flow and heat transfer within Newtonian liquid film over an unsteady stretching sheet with and without thermocapillarity is examined. The governing non‐linear partial differential equations describing the problem are reduced to a system of nonlinear ordinary differential equations using similarity transformation, which is solved numerically for different values of the thermal radiation parameter and the thermocapillarity parameter. The results show that the dimensionless velocity, the film thickness and the local Nusselt number increase as the thermocapillarity parameter increases, while the free surface temperature decreases with increasing the thermocapillarity parameter. Also, both the dimensionless temperature and the free surface temperature increase and the local Nusselt number decreases as the thermal radiation parameter increases. Copyright © 2010 John Wiley & Sons, Ltd.     

Droplet heat transfer on micro-post arrays: Effect of droplet size on droplet thermal characteristics

Heat transfer towards a water droplet from hydrophobic micro-post array surface is considered while mimicking the environmental temperatures. Micro-post arrays are created on a silicon wafer surface via lithography technique. The textured surfaces are replicated by polydimethylsiloxane (PDMS) to achieve an optical transmittance. The droplet adhesion on micro-post array surface is presented and the influence of droplet size on the heat transfer and droplet internal flow characteristics is examined. The flow predictions are validated via the particle image velocimetry data. It is found that adhesion force between the water droplet and the micro-post arrays surface depends on the geometric size and the orientation of the micro-post arrays on the surface. Temperature and flow fields are influenced by the droplet size. The Nusselt and the Bond numbers increase with the droplet volume; however, the Bond number remains less than unity indicating that the Marangoni current dominates over the buoyancy current in the droplet. The Nusselt number attains larger values for micro-post array surface than that of the plain surface. This is because of temperature and velocity oscillations along the contact lines at the droplet bottom due to the pitches of the micro-post arrays.     

Particle velocimetry inside Newtonian and non-Newtonian droplets impacting a hydrophobic surface

A particle velocimetry technique is described which enables the measurement of the fluid velocity inside impacting drops. Using high speed photography of 2 μm fluorescent tracer particles suspended in the fluid, the velocity field was measured as a function of time and radial position. The potential of the technique is illustrated using velocimetry measurements of drops of pure water and aqueous solutions of 200 ppm poly-(ethylene oxide) (PEO). Dilute solutions of PEO have been known for some time to suppress the rebound of water from hydrophobic surfaces. The dissipation has traditionally been attributed to an increased extensional viscosity as the polymers stretch in the extensional flow of the droplet. Our results enable us to infer that the extensional viscosity of PEO drops, during both the spreading and retraction phase, is similar to that of pure water. The data suggest that the true source of dissipation lies at the droplet edge. We also show, by analysing the spreading of water drops, that the Roisman-Yarin theory for a droplet spreading on a surface is valid in the bulk of the droplet prior to the final stages of spreading.     

Atomization of liquids by disintegrating thin liquid films using gas jets

Liquid atomization is useful in many applications, such as engineering, science, pharmaceutics, medicine, forensics and others. In the present research, an innovative methodology and a new device for atomization of liquids into mists of micron and submicron droplets have been developed. The new liquid-atomization method exploits the physical phenomenon of fragmentation of thin liquid films into fine micron and submicron droplets by gas jets. For several tested prototypes, the direct observations using a high-speed visualization technique have demonstrated that bubbles were generated within a liquid and their shells have been subsequently destroyed by applying a mechanical impulse (pressure of a compressed air) once the bubbles came over the liquid surface. The main characteristics of the generated tap water mists have been experimentally measured by means of the laser diffraction technique under various conditions for each prototype. One of the prototype devices allowed obtaining mists containing 90–99% of droplets smaller than 1 µm, with the minimum arithmetic and Sauter mean droplet diameters of 1.48 µm and 2.66 µm, and the 2.64 ml/min of droplet flow rate for 3.5 bar manometer pressure of atomizing air. The gas to liquid mass ratios (GLR) in the new device are depending on the atomizing tube length and the number of perforated orifices in the tube: more the tube length, hence more the number of perforated orifices, and therefore more liquid droplets will form for the same gas flow rate. The measured GLR values related to 1 m length of the utilized atomizing tube were in the range of 0.65–1.06, and for the specifically utilized atomizing tube of 72 mm length were among 9.07–14.67. The results of this study demonstrate that the developed method of generation of very fine droplet mists has many advantages over the existing techniques and can be perspective for many practical applications.     

Experimental study of air natural convection on metallic foam-sintered plate

The natural convection on metallic foam-sintered plate at different inclination angles was experimentally studied. Seven copper foam samples with different pore densities (10–40 pore per inch), porosities (0.90–0.95), and aspect ratios (the ratio of foam thickness to sample length, 0.1–0.5) were measured at inclination angles of 0° (vertical orientation), 15°, 30°, 45°, 60°, 75°, 90° (horizontal orientation). The heat conduction and natural convection inside the foam both contributed to the total heat transfer. Although, the form and viscous drag, which are influenced by permeability and viscous friction in the thermal boundary layer respectively, tend to suppress the natural convection, the heat transfer was finally enhanced by the foam sintered surface due to large surface area extension. Optimum inclination range 60–75° corresponding to maximum average Nu number was found in the heat flux range of 600–1800 W/m². The sintered foam surface with lower porosity and pore density was recommended for heat transfer enhancement. Particularly, the sample with porosity 0.9, pore density of 10 PPI, aspect ratio of 0.5 offered the highest average Nu number among the studied samples. An empirical correlation for modified Nusselt number at isoflux boundary condition considering the foam morphology parameter and inclination angle was proposed within deviation ±15% between the correlation and the experimental data.     

Heat transfer and fluid flow analysis of artificially roughened ducts having rib and groove roughness

Artificially roughness is one of the well known methods of enhancing heat transfer from the heat transfer surface in the form of repeated ribs, grooves or combination of ribs and groove (compound turbulators). The artificial roughness produced on the heat transferring surface is used in cooling of gas turbine blades, nuclear reactor, solar air heating systems etc. Solar air heaters have wide applications in low to moderate temperature range, namely, drying of foods, agricultural crops, seasoning of wood and space heating etc. Solar air heaters have low value of convective heat transfer coefficient between the working fluid (air) and the heat transferring surface, due to the formation of thin laminar viscous sub-layer on its surface. The heat transfer from the surface can be increased by breaking this laminar viscous sub layer. Hence, in the present work compound turbulators in the form of integral wedge shaped ribs with grooves are used on the heat transfer surface, to study its effect on the heat transfer coefficient (Nusselt number) and friction factor in the range of Reynolds number 3,000–18,000. The roughness produced on the absorber plate forms the wetted side of upper broad wall of the rectangular duct of solar air heater. The relative groove position (g/p) was varied from 0.4 to 0.8 and the wedge angle (Φ) was varied from 10° to 25°, relative roughness pitch (p/e) and relative roughness height (e/D) was maintained as 8.0 and 0.033 respectively. The aspect ratio of the rectangular duct was maintained as 8. The Nusselt number and friction factor of the artificially roughened ducts were determined experimentally and the corresponding values were compared with that of smooth surface duct. It is observed that wedge-groove roughened surface shows more enhancement in heat transfer compared to only rib roughened surface arrangement. The investigation revealed that Nusselt number increases 1.5–3 times, while the friction factor increases two to three folds that of the smooth surface duct in the range of operating parameters. It is also observed that in rib–groove roughness arrangement with relative groove position of 0.65 shows the maximum enhancement in the heat transfer compared to the other rib-groove roughness arrangements. Statistical correlations for the Nusselt number and friction factor have been developed by the regression method in terms of the operating and roughness parameters. A program was also developed in MATLAB for the calculation of thermal efficiency and thermal effectiveness. It was observed that the thermal efficiency is more for wedge angle of 15° and relative groove position of 0.65 and its value ranges from 42 to 73 %. The uncertainties in the measurements due to various instruments for the Reynolds number, Nusselt number, and friction factor have been estimated as ±3.8, ±4.54 and ±7.6 % respectively in the range of investigation made.     

Computational study of bouncing and non-bouncing droplets impacting on superhydrophobic surfaces

We numerically investigate bouncing and non-bouncing of droplets during isothermal impact on superhydrophobic surfaces. An in-house, experimentally validated, finite element method-based computational model is employed to simulate the droplet impact dynamics and transient fluid flow within the droplet. The liquid–gas interface is tracked accurately in Lagrangian framework with dynamic wetting boundary condition at three-phase contact line. The interplay of kinetic, surface and gravitational energies is investigated via systematic variation of impact velocity and equilibrium contact angle. The numerical simulations demonstrate that the droplet bounces off the surface if the total droplet energy at the instance of maximum recoiling exceeds the initial surface and gravitational energy, otherwise not. The non-bouncing droplet is characterized by the oscillations on the free surface due to competition between the kinetic and surface energy. The droplet dimensions and shapes obtained at different times by the simulations are compared with the respective measurements available in the literature. Comparisons show good agreement of numerical data with measurements, and the computational model is able to reconstruct the bouncing and non-bouncing of the droplet as seen in the measurements. The simulated internal flow helps to understand the impact dynamics as well as the interplay of the associated energies during the bouncing and non-bouncing. A regime map is proposed to predict the bouncing and non-bouncing on a superhydrophobic surface with an equilibrium contact angle of 155°, using data of 86 simulations and the measurements available in the literature. We discuss the validity of the computational model for the wetting transition from Cassie to Wenzel state on micro- and nanostructured superhydrophobic surfaces. We demonstrate that the numerical simulation can serve as an important tool to quantify the internal flow, if the simulated droplet shapes match the respective measurements utilizing high-speed photography.     

Pool boiling enhancement with surface treatments

Pool boiling heat transfer characteristics on treated surfaces were investigated experimentally. Surface treatments were performed with sandpapers, building micro structures by etching, and micro-porous coating. Copper blocks of 20 mm × 20 mm were used as test sections and PF5060 was used as a working fluid. The effect of wall superheat (0–35 K), surface orientation (0°, 45°, 90°), and subcooling (0, 5, 10 K) on treated surfaces were also investigated. Heat transfer performance on inclined surface was better than that of horizontal surface in all test sections. Because bubble generation was suppressed by subcooling, higher wall superheat was needed to initiate boiling in all surfaces. Micro-porous coated surface showed the highest heat transfer enhancement among tested surfaces.     

Insight into the interaction between water and surfactant aerosol particles based on molecular dynamics simulations

The ubiquitous surfactant significantly influences the hygroscopic growth of atmospheric aerosol particles. However, knowledge on the morphology of surfactant particles after the adsorption of water is insufficient. In this study, the interaction between water and particles composed of surface active malonic acid (MA) or adipic acid (AA) are simulated based on the molecular dynamics method. The key point is the combined effect of temperature and water content on the structural properties of particles and the surface propensity of surfactants at the equilibrium state. Results show that demixed structure 1 with the adsorption of water clusters on acid grain, mixed structure and demixed structure 2 with acids coating on water droplet can be observed. With temperature increasing from 160 K to 330 K the surface propensity of MA and AA increases first and then weakens. Near the standard atmospheric temperature (280–330 K), the surface propensity of MA and AA increases with increasing water content and alkyl group, and its sensitivity to temperature and water content varies regularly. Moreover, all surfactants at the particle surface orient their hydrophobic groups toward the gas. These findings improve our insight into the surfactant partitioning and further assist in more accurate prediction of the particle hygroscopic growth.     

Unsteady MHD mixed convection flow over an impulsively stretched permeable vertical surface in a quiescent fluid

The unsteady mixed convection flow of an incompressible laminar electrically conducting fluid over an impulsively stretched permeable vertical surface in an unbounded quiescent fluid in the presence of a transverse magnetic field has been investigated. At the same time, the surface temperature is suddenly increased from the surrounding fluid temperature or a constant heat flux is suddenly imposed on the surface. The problem is formulated in such a way that for small time it is governed by Rayleigh type of equation and for large time by Crane type of equation. The non-linear coupled parabolic partial differential equations governing the unsteady mixed convection flow under boundary layer approximations have been solved analytically by using the homotopy analysis method as well as numerically by an implicit finite difference scheme. The local skin friction coefficient and the local Nusselt number are found to decrease rapidly with time in a small time interval and they tend to steady-state values for t*≥5. They also increase with the buoyancy force and suction, but decrease with injection rate. The local skin friction coefficient increases with the magnetic field, but the local Nusselt number decreases. There is a smooth transition from the unsteady state to the steady state.     

Measurement of extensional viscosity of polymer solutions

A filament stretching technique for measuring the extensional viscosity of polymer solutions at constant stretch rate is presented. The liquid sample is held between two coaxial discs and stretched by moving the bottom disc downwards with a speed that increases exponentially with time. This is illustrated using a constant viscosity, elastic fluid consisting of 0.185% polyisobutylene in a solvent of kerosene and polybutene. For the case of this particular fluid, two distinct stretch rate regions are found to arise. The stretch rate in the first region is much higher than in the second, which is, in most cases, close to the overall stretch rate imposed on the sample. Nonetheless, all the results of any given run can be represented using an average extensional rate. The extensional stress growth data, plotted as the Trouton ratio against time, show an initial linear viscoelastic region where T_R rises to a value of 3, independent of extensional rate. Beyond this region, T_R depends on the stretch rate and rises dramatically to values in excess of 10³; the higher the extensional rate, the faster is the increase in T_R. These data do not seem to reach a steady state and appear to be similar to polymer melt data obtained by others in the past. The reproducibility of the results is very good and all this suggests that it is now possible to obtain unambiguous constant-stretch-rate stress-growth data for polymer solutions stretched from a state of rest.     

Experimental study of slot jet impingement heat transfer on a wedge-shaped surface

An experimental investigation was conducted to study the convective heat transfer rate from a wedge-shaped surface to a rectangular subsonic air jet impinging onto the apex of the wedge. The jet Reynolds number, nozzle-to-surface distance and the wedge angle were considered as the main parameters. Jet Reynolds number was ranged from 5,000 to 20,000 and two dimensionless nozzle-to-surface distances h/w?=?4 and 10 were examined. The apex angle of the wedge ranged from 30° to 180° where the latter case corresponds with that of a flat surface. Velocity profile and turbulence intensity were provided for free jet flow using hot wire anemometer. Local and average Nusselt numbers on the impinged surface are presented for all the configurations. Based on the results presented, the local Nusselt number at the stagnation region increases as the wedge angle is decreased but, it then decreases over the remaining area of the impinged surface. Average Nusselt number over the whole surface is maximum when the wedge angle is 180° (i.e. plane surface) for any jet and nozzle-to-surface configuration.     

Boundary layer flow of a Jeffrey fluid with convective boundary conditions

The boundary layer stretched flow of a Jeffrey fluid subject to the convective boundary conditions was investigated. The governing dimensionless problems were computed by using the homotopy analysis approach. Convergence of the derived solutions was checked and the influence of embedded parameters was analyzed by plotting graphs. It was noticed that the velocity increases with an increase in the Deborah number. Furthermore, it was found that the temperature is also an increasing function of the Biot number. We further found that for fixed values of other parameters, the local Nusselt number increases by increasing the suction parameter and Deborah number. Numerical values of the skin friction coefficient and local Nusselt numbers were computed and examined. Copyright © 2011 John Wiley & Sons, Ltd.     

Atomization of glycerin with a twin-fluid swirl nozzle

Atomization of liquids with high viscosity is always a challenge, especially when small diameter droplets and high liquid flow rates are simultaneously required. In the present research, the performance of a Venturi–vortex twin-fluid swirl nozzle is examined, attending to its capabilities to generate droplets with diameters below 20 µm when atomizing pure glycerin at room temperature. In this nozzle, air is injected tangentially in a central convergent section, and discharges suctioning the liquid fed to a coaxial chamber, here using a gear pump. The resulting spray is visualized and analyzed. Droplet size distributions are measured with a laser diffractometer. As expected, droplet diameter increases with liquid flow rate, and quickly diminishes when air flow rate is increased. Sauter mean diameters (SMD) below 15 µm can be obtained even when atomizing pure glycerin. However, these values are obtained for relatively low glycerin flow rates (∼5 l/h), and with rather wide distributions. For 10 l/h and an air-to-liquid mass flow rate ratio (ALR) of 13.7 more than 26% of the glycerin volume is atomized in droplets smaller than 20 µm. Liquid ligaments are observed near the nozzle exit, but they tend to break up while moving downstream.     

Drag force on a free surface-piercing yawed circular cylinder in steady flow

The force distribution on a surface-piercing yawed cylinder surface differs significantly from that on a surface-piercing vertical cylinder. The established numerical model for flow past the surface-piercing yawed cylinder with yaw angles from −45° to 45° was solved by the standard large-eddy simulation (LES) methodology. Six cases at intervals of ±15° relative to the vertical were studied at the Reynolds number of 27 000 and the Froude number of 0.8 based on the cylinder diameter and free-stream velocity, among which the drag forces on four cylinders with yaw angles from −15° to 30° were tested for the validation of the LES approach. The results revealed that the time-averaged total drag coefficient for all cases increases with the increase of yaw angle compared to that of the surface-piercing vertical cylinder, even over 2.5 for the ±45°-yawed cylinders. The sectional drag coefficients for the negatively yawed cylinders are much greater than that for the vertical cylinder, and much less for the positively yawed cylinders. The unbalanced hydrostatic pressures on the inclined section are mainly responsible for those increment and decrement. Once the hydrostatic pressure was removed, the sectional drag coefficient on the mid-span of the positively yawed cylinder increases from the top section to the bottom, and decreases for the negatively yawed cylinder. The corresponding integrated total drag coefficient decreases with the increase of the yaw angle to ±15°, then increases with the further increase of the magnitude of yaw angle.