Prashant Valluri

Convective rolls and hydrothermal waves in evaporating sessile drops.

Recent experiments on the evaporation of sessile droplets have revealed the spontaneous formation of various patterns including the presence of hydrothermal waves. These waves had previously been observed, in the absence of evaporation, in thin liquid layers subjected to an imposed, uniform temperature gradient. This is in contrast to the evaporating droplet case wherein these gradients arise naturally due to evaporation and are spatially and temporally varying. In the present paper, we present a theory of evaporating sessile droplets deposited on a heated surface and propose a candidate mechanism for the observed pattern formation using a linear stability analysis in the quasi-steady-state approximation. A qualitative agreement with experimental trends is observed.

Linear and nonlinear stability of hydrothermal waves in planar liquid layers driven by thermocapillarity

A shallow planar layer of liquid bounded from above by gas is set into motion via the thermocapillary effect resulting from a thermal gradient applied along its interface. Depending on the physical properties of the liquid and the strength of the gradient, the system is prone to departure from its equilibrium state and to the consequent development of an oscillatory regime. This problem is numerically investigated for the first time by means of two-phase direct numerical simulations fully taking into account the presence of a deformable interface. Obliquely travelling hydrothermal waves (HTWs), similar to those first described by Smith and Davis [J. Fluid Mech. 132, 119–144 (1983)]10.1017/S0022112083001512, are reported presenting good agreement with linear stability theory and experiments. The nonlinear spatiotemporal growth of the instabilities is discussed extensively along with the final bulk flow for both the liquid and gas phases. Our study reveals the presence of interface deformations which accomp...

On phase change in Marangoni-driven flows and its effects on the hydrothermal-wave instabilities

This paper investigates the effects of phase change on the stability of a laterally heated liquid layer for the first time. The interface is open to the atmosphere and vapor diffusion is the rate-limiting mechanism for evaporation. In this configuration, the planar layer is naturally vulnerable to the formation of travelling thermal instabilities, i.e., hydrothermal waves (HTWs), due to the presence of temperature gradients along the gas-liquid interface. Recent work carried out for deformable interfaces and negligible evaporation indicates that the HTWs additionally give rise to interface deformations of similar features, i.e., physical waves. The study presented here reveals that phase change plays a dual role through its effect on these instabilities: the latent energy required during the evaporation process tends to inhibit the HTWs while the accompanying level reduction enhances the physical waves by minimizing the role of gravity. The dynamics of the gas phase are also discussed. The HTW-induced convective patterns in the gas along with the travelling nature of the instabilities have a significant impact on the local evaporation flux and the vapor distribution above the interface. Interestingly, high (low) concentrations of vapor are found above cold (hot) spots. The phase-change mechanism for stable layers is also investigated. The Marangoni effect plays a major role in the vapor distribution generating a vacuum effect in the warm region and vapor accumulations at the cold boundary capable of inverting the phase change, i.e., the capillary flow can lead to local condensation. This work also demonstrates the inefficiencies of the traditional phase change models based on pure vapor diffusion to capture the dynamics of thermocapillary flows.

Linear and nonlinear instability in vertical counter-current laminar gas-liquid flows

We consider the genesis and dynamics of interfacial instability in vertical gas-liquid flows, using as a model the two-dimensional channel flow of a thin falling film sheared by counter-current gas. The methodology is linear stability theory (Orr-Sommerfeld analysis) together with direct numerical simulation of the two-phase flow in the case of nonlinear disturbances. We investigate the influence of two main flow parameters on the interfacial dynamics, namely the film thickness and pressure drop applied to drive the gas stream. To make contact with existing studies in the literature, the effect of various density contrasts is also examined. Energy budget analyses based on the Orr-Sommerfeld theory reveal various coexisting unstable modes (interfacial, shear, internal) in the case of high density contrasts, which results in mode coalescence and mode competition, but only one dynamically relevant unstable interfacial mode for low density contrast. A study of absolute and convective instability for low densi...

High-performance computational fluid dynamics: a custom-code approach

We introduce a modified and simplified version of the pre-existing fully parallelized three-dimensional Navier–Stokes flow solver known as TPLS. We demonstrate how the simplified version can be used as a pedagogical tool for the study of computational fluid dynamics (CFDs) and parallel computing. TPLS is at its heart a two-phase flow solver, and uses calls to a range of external libraries to accelerate its performance. However, in the present context we narrow the focus of the study to basic hydrodynamics and parallel computing techniques, and the code is therefore simplified and modified to simulate pressure-driven single-phase flow in a channel, using only relatively simple Fortran 90 code with MPI parallelization, but no calls to any other external libraries. The modified code is analysed in order to both validate its accuracy and investigate its scalability up to 1000 CPU cores. Simulations are performed for several benchmark cases in pressure-driven channel flow, including a turbulent simulation, wherein the turbulence is incorporated via the large-eddy simulation technique. The work may be of use to advanced undergraduate and graduate students as an introductory study in CFDs, while also providing insight for those interested in more general aspects of high-performance computing.

Manufacturing of microcirculation phantoms using rapid prototyping technologies

In this paper, we describe a method for the manufacturing of a microcirculation phantom that may be used to investigate hemodynamics using optics based methods. We made an Acrylonitrile Butadiene Styrene (ABS) negative mold, manufactured in a Fused Deposition Modelling (FDM) printer, embedded it in Polydimethysilioxane (PDMS) and dissolved it from within using acetone. We successfully made an enlarged three-dimensional (3D) network of microcirculation, and tested it using red blood cell (RBC) analogues. This phantom may be used for testing medical imaging technology.

How does blood regulate cerebral temperatures during hypothermia

Macro-modeling of cerebral blood flow can help determine the impact of thermal intervention during instances of head trauma to mitigate tissue damage. This work presents a bioheat model using a 3D fluid-porous domain coupled with intersecting 1D arterial and venous vessel trees. This combined vascular porous (VaPor) model resolves both cerebral blood flow and energy equations, including heat generated by metabolism, using vasculature extracted from MRI data and is extended using a tree generation algorithm. Counter-current flows are expected to increase thermal transfer within the brain and are enforced using either the vascular structure or flow reversal, represented by a flow reversal constant, CR. These methods exhibit larger average brain cooling (from 0.56 °C ± <0.01 °C to 0.58 °C ± <0.01 °C) compared with previous models (0.39 °C) when scalp temperature is reduced. An greater reduction in core brain temperature is observed (from 0.29 °C ± <0.01 °C to 0.45 °C ± <0.01 °C) compared to previous models (0.11 °C) due to the inclusion of counter-current cooling effects. The VaPor model also predicts that a hypothermic average temperature (<36 °C) can be reached in core regions of neonatal models using scalp cooling alone.

Parallel Computing: On the Road to Exascale

We introduce TPLS (Two-Phase Level Set), an MPI-parallel Direct Numerical Simulation code for two-phase flows in channel geometries. Recent developments to the code are discussed which improve the performance of the solvers and I/O by using the PETSc and NetCDF libraries respectively. Usability and functionality improvements enabled by code refactoring and merging of a separate OpenMP-parallelized version are also outlined. The overall scaling behaviour of the code is measured, and good strong scaling up to 1152 cores is observed for a 5.6 million element grid. A comparison is made between the legacy serial textformatted I/O and new NetCDF implementations, showing speedups of up to 17x. Finally, we explore the effects of output file striping on the Lustre parallel file system on ARCHER, a Cray XC30 supercomputer, finding performance gains of up to 12% over the default striping settings.