Source author record

Rajneesh Bhardwaj

Rajneesh Bhardwaj appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2020arXiv

A finite-element model for computing fluid flow inside a sessile evaporating droplet on a solid surface

A finite element model was developed to compute the fluid flow inside a sessile evaporating droplet on hydrophilic substrate in ambient conditions. The evaporation is assumed as quasi-steady and the flow is considered as axisymmetric with a pinned contact line. The Navier-Stokes equations in cylindrical coordinates were solved inside the droplet. Galerkin weight residual approach and velocity pressure formulation was used to discretise the governing equations. Six node triangular mesh and quadratic shape functions were used to obtain higher accuracy solutions. Radial velocity profiles in axial directions calculated by the FEM solver were compared with a existing analytical model and were found in excellent agreement. The contours of velocity magnitude and streamlines show the characteristic flow i.e. radially outward inside the evaporating droplet.

preprint2020arXiv

Computational Modeling and Analysis of Flow-induced Vibration of an Elastic Splitter Plate Using a Sharp-interface Immersed Boundary Method

We present the development and benchmarking of an in-house fluid-structure interaction (FSI) solver. An implicit partitioned approach is utilized to couple a sharp-interface immersed boundary (IB) method based flow solver and a finite-element method based structural solver. In the present work, the coupling is accelerated using a dynamic under-relaxation scheme. The revised coupling is around two to three times faster and numerically stable, as compared to the one that uses a constant under-relaxation parameter. The solver is validated against two FSI benchmarks in which a thin, finite thickness, elastic splitter plate is attached to the lee side of a circular or square rigid cylinder, subjected to laminar flow. In these two-dimensional benchmarks, the flow induces a wave-like deformation in the plate, and it attains a periodic self-sustained oscillation. We employ the FSI solver to analyze the flow-induced vibration (FIV) of the plate in a uniform laminar free-stream flow for a wide range of mass ratio and bending stiffness at Reynolds number ($Re$) of 100, based on the diameter of the cylinder. At the given $Re$, two-dimensional numerical simulations show that the FIV of the plate effectively depends only on the mass ratio and bending stiffness. The largest displacement of the plate vibration is found to occur in the lock-in region, where the vortex shedding frequency of the coupled fluid-structure system is close to the natural frequency of the splitter plate. We briefly discuss wake structures and phase plots for different cases of mass ratio and bending stiffness.

preprint2020arXiv

Evaporation of initially heated sessile droplets and the resultant dried colloidal deposits on substrates held at ambient temperature

The present study experimentally and numerically investigates the evaporation and resultant patterns of dried deposits of aqueous colloidal sessile droplets, when the droplets are initially elevated to a high temperature before being placed on a substrate held at ambient temperature. The system is then released for natural evaporation without applying any external perturbation. Infrared thermography and optical profilometry were used as essential tools for interfacial temperature measurements and quantification of the coffee-ring dimensions, respectively. Initially, a significant temperature gradient exists along the liquid-gas interface as soon as the droplet is deposited on the substrate which triggers a Marangoni stress-induced recirculation flow directed from the top of the droplet towards the contact line along the liquid-gas interface. Thus, the flow is in the reverse direction to that seen in the conventional substrate heating case. Interestingly, this temperature gradient decays rapidly -- within the first 10% of the total evaporation time and the droplet-substrate system reaches thermal equilibrium with ambient thereafter. Despite fast decay of the temperature gradient, the coffee-ring dimensions significantly diminish, leading to an inner deposit. This suppression of the coffee-ring effect is attributed to the fact that the initial Marangoni stress-induced recirculation flow continues until the last stage of the evaporation, even after the interfacial temperature gradient vanishes. This is essentially a consequence of liquid inertia. Overall, together with a new experimental condition, the present investigation discloses a distinct nature of Marangoni stress-induced flow in the drying droplet and its role in influencing the associated colloidal deposits, which was not explored previously.

preprint2020arXiv

Likelihood of survival of coronavirus in a respiratory droplet deposited on a solid surface

We predict and analyze the drying time of respiratory droplets from a COVID-19 infected subject, which is a crucial time to infect another subject. The drying of the droplet is predicted by diffusion-limited evaporation model for a sessile droplet placed on a partially-wetted surface with a pinned contact line. The variation of droplet volume, contact angle, ambient temperature, and humidity are considered. We analyze the chances of the survival of the viruses present in the droplet, based on the lifetime of the droplets in several conditions, and find that the chances of survival of the virus are strongly affected by each of these parameters. The magnitude of shear stress inside the droplet computed using the model is not large enough to obliterate the virus. We also explore the relationship between the drying time of a droplet and the growth rate of the spread of COVID-19 for five different cities, and find that they are weakly correlated.

preprint2016arXiv

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.

preprint2016arXiv

Droplet Impact Dynamics on Micropillared Hydrophobic Surfaces

The effect of pitch of the pillars and impact velocity are studied for the impact dynamics of a microliter water droplet on a micropillared hydrophobic surface. The results are presented qualitatively by the high-speed photography and quantitatively by the temporal variation of wetted diameter and droplet height. A characterization of the transient quantitative results is a novel aspect of our work. Three distinct regimes, namely, non-bouncing, complete bouncing and partial bouncing are presented. A critical pitch as well as impact velocity exists for the transition from one regime to another. This is explained with a demonstration of Cassie to Wenzel wetting transition in which the liquid penetrates in the grooves between the pillars at larger pitch or impact velocity. The regimes are demarcated on a map of pitch and impact velocity. A good agreement is reported between the present measurements and published analytical models.

preprint2016arXiv

Effects of Substrate Heating and Wettability on Evaporation Dynamics and Deposition Patterns for a Sessile Water Droplet Containing Colloidal Particles

Effects of substrate temperature, substrate wettability and particles concentration are experimentally investigated for evaporation of a sessile water droplet containing colloidal particles. Time-varying droplet shapes and temperature of the liquid-gas interface are measured using high-speed visualization and infrared thermography, respectively. The motion of the particles inside the evaporating droplet is qualitatively visualized by an optical microscope and profile of final particle deposit is measured by an optical profilometer. On a non-heated hydrophilic substrate, a ring-like deposit forms after the evaporation, as reported extensively in the literature; while on a heated hydrophilic substrate, a thinner ring with an inner deposit is reported in the present work. The latter is attributed to Marangoni convection and recorded motion of the particles as well as measured temperature gradient across the liquid-gas interface confirms this hypothesis. The thinning of the ring scales with the substrate temperature and is reasoned to stronger Marangoni convection at larger substrate temperature. In case of a non-heated hydrophobic substrate, an inner deposit forms due to very early depinning of the contact line. On the other hand, in case of a heated hydrophobic substrate, the substrate heating as well as larger particle concentration helps in the pinning of the contact line, which results in a thin ring with an inner deposit. We propose a regime map for predicting three types of deposits namely, ring, thin ring with inner deposit and inner deposit - for varying substrate temperature, substrate wettability and particles concentration. A first-order model corroborates the liquid-gas interface temperature measurements and variation in the measured ring profile with the substrate temperature.

preprint2016arXiv

Numerical study of Heat Transfer Enhancement by Deformable Twin Plates in Laminar Heated Channel flow

Fluid-structure interaction (FSI) of thin flexible fins coupled with convective heat transfer has applications in energy harvesting and in understanding functioning of several biological systems. We numerically investigate FSI of the thin flexible fins involving large-scale flow-induced deformation as a potential heat transfer enhancement technique. An in-house, strongly-coupled fluid-structure interaction (FSI) solver is employed in which flow and structure solvers are based on sharp-interface immersed boundary and finite element method, respectively. We consider twin flexible fins in a heated channel with laminar pulsating cross flow. The vortex ring past the fin sweep higher sources of vorticity generated on the channel walls out into the downstream - promoting the mixing of the fluid. The moving fin assists in convective mixing, augmenting convection in bulk and at the walls; and thereby reducing thermal boundary layer thickness and improving heat transfer at the channel walls. The thermal augmentation is quantified in term of instantaneous Nusselt number at the wall. Results are presented for two limiting cases of thermal conductivity of the fin - an insulated fin and a highly conductive fin. We discuss feasibility of flexible fins and effect of flow cycle-to-cycle variation on the Nusselt number. Finally, we investigate the effect of important problem parameters - Young's Modulus, flow frequency and Prandtl number - on the thermal augmentation.

preprint2016arXiv

The Response of an Elastic Splitter Plate Attached to a Cylinder to Laminar Pulsatile Flow

The flow-induced deformation of a thin, elastic splitter plate attached to the rear of a circular cylinder and subjected to laminar pulsatile inflow is investigated. The cylinder and elastic splitter plate are contained within a narrow channel and the Reynolds number is mostly restricted to Re = 100, primarily covering the two-dimensional flow regime. An in-house fluid-structure interaction code is employed for simulations, which couples a sharp-interface immersed boundary method for the fluid dynamics with a finite-element method to treat the structural dynamics. The structural solver is implicitly (two-way) coupled with the flow solver using a partitioned approach. This implicit coupling ensures numerical stability at low structure-fluid density ratios. A power spectrum analysis of the time-varying plate displacement shows that the plate oscillates at more than a single frequency for pulsatile inflow, compared to a single frequency observed for steady inflow. The multiple frequencies obtained for the former case can be explained by beating between the applied and plate oscillatory signals. The plate attains a self-sustained time-periodic oscillation with a plateau amplitude in the case of steady flow, while the superimposition of pulsatile inflow with induced plate oscillation affects the plateau amplitude. Lock-in of the plate oscillation with the pulsatile inflow occurs at a forcing frequency that is twice of the plate natural frequency in a particular mode and this mode depends on the plate length. The plate displacement as well as pressure drag increases at the lock-in condition. The percentage change in the maximum plate displacement, and skin-friction and pressure drag coefficients on the plate, due to pulsatile inflow is quantified. The non-linear dynamics of the plate and its coupling with the pulsatile flow are briefly discussed.

preprint2015arXiv

Flow-Induced Deformation of a Flexible Thin Structure as Manifestation of Heat Transfer Enhancement

Flow-induced deformation of thin structures coupled with convective heat transfer has potential applications in energy harvesting and is important for understanding functioning of several biological systems. We numerically demonstrate large-scale flow-induced deformation as an effective passive heat transfer enhancement technique. An in-house, strongly-coupled fluid-structure interaction (FSI) solver is employed in which flow and structure solvers are based on sharp-interface immersed boundary and finite element method, respectively. In the present work, we validate convective heat transfer module of the in-house FSI solver against several benchmark examples of conduction and convective heat transfer including moving structure boundaries. The thermal augmentation is investigated as well as quantified for the flow-induced deformation of an elastic thin plate attached to lee side of a rigid cylinder in a heated channel laminar flow. We show that the wake vortices past the plate sweep higher sources of vorticity generated on the channel walls out into the high velocity regions, promoting the mixing of the fluid. The self-sustained motion of the plate assists in convective mixing, augmenting convection in bulk and near the walls, and thereby reducing thermal boundary layer thickness as well as improving Nusselt number at the channel walls. We quantify the thermal improvement with respect to channel flow without any bluff body and analyze the role of Reynolds number, Prandtl number and material properties of the plate in the thermal augmentation.

preprint2010arXiv

A numerical investigation on the influence of liquid properties and interfacial heat transfer during microdroplet deposition onto a glass substrate

This work investigates the impingement of a liquid microdroplet onto a glass substrate at different temperatures. A finite-element model is applied to simulate the transient fluid dynamics and heat transfer during the process. Results for impingement under both isothermal and non-isothermal conditions are presented for four liquids: isopropanol, water, dielectric fluid (FC-72) and eutectic tin-lead solder (63Sn-37Pb). The objective of the work is to select liquids for a combined numerical and experimental study involving a high resolution, laser-based interfacial temperature measurement to measure interfacial heat transfer during microdroplet deposition. Applications include spray cooling, micro-manufacturing and coating processes, and electronics packaging. The initial droplet diameter and impact velocity are 80 μm and 5 m/s, respectively. For isothermal impact, our simulations with water and isopropanol show very good agreement with experiments. The magnitude and rates of spreading for all four liquids are shown and compared. For non-isothermal impacts, the transient drop and substrate temperatures are expressed in a non-dimensional way. The influence of imperfect thermal contact at the interface between the drop and the substrate is assessed for a realistic range of interfacial Biot numbers. We discuss the coupled influence of interfacial Biot numbers and hydrodynamics on the initiation of phase change.

preprint2010arXiv

Interfacial temperature measurements, high-speed visualization and finite-element simulations of droplet impact and evaporation on a solid surface

The objective of this work is to investigate the coupling of fluid dynamics, heat transfer and mass transfer during the impact and evaporation of droplets on a heated solid substrate. A laser-based thermoreflectance method is used to measure the temperature at the solid-liquid interface, with a time and space resolution of 100 μs and 20 μm, respectively. Isopropanol droplets with micro- and nanoliter volumes are considered. A finite-element model is used to simulate the transient fluid dynamics and heat transfer during the droplet deposition process, considering the dynamics of wetting as well as Laplace and Marangoni stresses on the liquid-gas boundary. For cases involving evaporation, the diffusion of vapor in the atmosphere is solved numerically, providing an exact boundary condition for the evaporative flux at the droplet-air interface. High-speed visualizations are performed to provide matching parameters for the wetting model used in the simulations. Numerical and experimental results are compared for the transient heat transfer and the fluid dynamics involved during the droplet deposition. Our results describe and explain temperature oscillations at the drop-substrate interface during the early stages of impact. For the first time, a full simulation of the impact and subsequent evaporation of a drop on a heated surface is performed, and excellent agreement is found with the experimental results. Our results also shed light on the influence of wetting on the heat transfer during evaporation.

preprint2010arXiv

Non-isothermal wetting during impact of millimeter size water drop on a flat substrate: numerical investigation and comparison with high-speed visualization experiments

The objective of this work is to develop and validate a numerical model to study wetting during the impact of millimeter-size drops on a flat, smooth, solid substrate under isothermal or non-isothermal conditions. A finite-element modeling is used to simulate the transient fluid dynamics and heat transfer, considering Laplace forces on the liquid-gas boundary. The Lagrangian scheme allows a very precise tracking of the free surface deformation. In this work, the numerical model is extended to account for a temperature-dependent viscosity and for dynamic wetting at the contact line. Numerical results are presented to study the influence of the kinetic wetting parameter on the wetting incipience and behavior. Our results show the influence of wetting on the spreading and the transient drop shape. Also, numerical results are compared with high-speed visualization, for cases of isothermal and non-isothermal impact. Matching between simulations and high-speed visualization allows the determination of the value of the kinetic wetting parameter. Our main finding is that warm drops spread more than cold drops because of a reduction of viscous forces, and not because of an increase of wetting.

preprint2010arXiv

Pattern formation during the evaporation of a colloidal nanoliter drop: a numerical and experimental study

An efficient way to precisely pattern particles on solid surfaces is to dispense and evaporate colloidal drops, as for bioassays. The dried deposits often exhibit complex structures exemplified by the coffee ring pattern, where most particles have accumulated at the periphery of the deposit. In this work, the formation of deposits during the drying of nanoliter colloidal drops on a flat substrate is investigated numerically and experimentally. A finite-element numerical model is developed that solves the Navier-Stokes, heat and mass transport equations in a Lagrangian framework. The diffusion of vapor in the atmosphere is solved numerically, providing an exact boundary condition for the evaporative flux at the droplet-air interface. Laplace stresses and thermal Marangoni stresses are accounted for. The particle concentration is tracked by solving a continuum advection-diffusion equation. Wetting line motion and the interaction of the free surface of the drop with the growing deposit are modeled based on criteria on wetting angles. Numerical results for evaporation times and flow field are in very good agreement with published experimental and theoretical results. We also performed transient visualization experiments of water and isopropanol drops loaded with polystyrene microsphere evaporating on respectively glass and polydimethylsiloxane substrates. Measured evaporation times, deposit shape and sizes, and flow fields are in very good agreement with the numerical results. Different flow patterns caused by the competition of Marangoni loops and radial flow are shown to determine the deposit shape to be either a ring-like pattern or a homogeneous bump.

preprint2010arXiv

Self-assembly of colloidal particles from evaporating droplets: role of DLVO interactions and proposition of a phase diagram

The shape of deposits obtained from drying drops containing colloidal particles matters for technologies such as inkjet printing, microelectronics and bioassay manufacturing. In this work, the formation of deposits during the drying of nanoliter drops containing colloidal particles is investigated experimentally with microscopy and profilometry, and theoretically with an inhouse finite-element code. The system studied involves aqueous drops containing titania nanoparticles evaporating on a glass substrate. Deposit shapes from spotted drops at different pH values are measured using a laser profilometer. Our results show that the pH of the solution influences the dried deposit pattern, which can be ring-like or more uniform. The transition between these patterns is explained by considering how DLVO interactions such as the electrostatic and van der Waals forces modify the particle deposition process. Also a phase diagram is proposed to describe how the shape of a colloidal deposit results from the competition between three flow patterns: a radial flow driven by evaporation at the wetting line, a Marangoni recirculating flow driven by surface tension gradients, and the transport of particles towards the substrate driven by DLVO interactions. This phase diagram explains three types of deposits commonly observed experimentally, such as a peripheral ring, a small central bump, or a uniform layer. Simulations and experiments are found in very good agreement.