We consider theoretically, computationally and experimentally spontaneous evaporation of water and isopropanol drops on smooth silicon wafers. In contrast to a number of previous works, the solid ...surface we consider is smooth and therefore the droplets' evolution proceeds without contact line pinning. We develop a theoretical model for evaporation of pure liquid drops that includes Marangoni forces due to the thermal gradients produced by non-uniform evaporation, and heat conduction effects in both liquid and solid phases. The key ingredient in this model is the evaporative flux. We consider two commonly used models: one based on the assumption that the evaporation is limited by the processes originating in the gas (vapour diffusion-limited evaporation), and the other one which assumes that the processes in the liquid are limiting. Our theoretical model allows for implementing evaporative fluxes resulting from both approaches. The required parameters are obtained from physical experiments. We then carry out fully nonlinear time-dependent simulations and compare the results with the experimental ones. Finally, we discuss how the simulation results can be used to predict which of the two theoretical models is appropriate for a particular physical experiment.
Metallic nanoparticles, liquified by fast laser irradiation, go through a rapid change of shape attempting to minimize their surface energy. The resulting nanodrops may be ejected from the substrate ...when the mechanisms leading to dewetting are sufficiently strong, as in the experiments involving gold nanoparticles Habenicht et al., Science 309, 2043 (2005). We use a direct continuum-level approach to accurately model the process of liquid nanodrop formation and the subsequent ejection from the substrate. Our computations show a significant role of inertial effects and an elaborate interplay of initial geometry and wetting properties: e.g., we can control the direction of ejection by prescribing appropriate initial shape and/or wetting properties. The basic insight regarding ejection itself can be reached by considering a simple effective model based on an energy balance. We validate our computations by comparing directly with the experiments specified above involving the length scales measured in hundreds of nanometers and with molecular dynamics simulations on much shorter scales measured in tens of atomic diameters, as by M. Fuentes-Cabrera et al. Phys. Rev. E 83, 041603 (2011). The quantitative agreement, in addition to illustrating how to control particle ejection, shows utility of continuum-based simulation in describing dynamics on nanoscale quantitatively, even in a complex setting as considered here.
We consider two (2D) and three (3D) dimensional granular systems exposed to compression, and ask what is the influence of the number of physical dimensions on the properties of the interaction ...networks that spontaneously form as these systems evolve. The study is carried out based on discrete element simulations of frictional disks in 2D and spheres in 3D. Within the constraints of the considered simulation protocols, the main finding is that both the number of physical dimensions and the type of particle-particle interaction significantly influence the properties of interaction networks. These networks play an important role in bridging the microscale (particle size) and macroscale (system size), thus both aspects (the interaction model and dimensionality) are carefully considered. Our work uses a combination of tools and techniques, including percolation study, statistical analysis, as well as algebraic topology-based techniques. In many instances, different techniques and measures provide complementary information that, when combined, allow for gaining better insight into the properties of interaction networks in compressed particulate systems.
Graphic Abstract
In this work we consider a new class of oscillatory instabilities that pertain to thermocapillary destabilization of a liquid film heated by a solid substrate. We assume the substrate thickness and ...substrate–film thermal conductivity ratio are large so that the effect of substrate thermal diffusion is retained at leading order in the long-wave approximation. As a result, the system dynamics is described by a nonlinear partial differential equation for the film thickness that is non-locally coupled to the full substrate heat equation. Perturbing about a steady quiescent state, we find that its stability is described by a non-self-adjoint eigenvalue problem. We show that, under appropriate model parameters, the linearized eigenvalue problem admits complex eigenvalues that physically correspond to oscillatory (in time) instabilities of the thin-film height. As the principal results of our work, we provide a complete picture of the susceptibility to oscillatory instabilities for different model parameters. Using this description, we conclude that oscillatory instabilities are more relevant experimentally for films heated by insulating substrates. Furthermore, we show that oscillatory instability where the fastest-growing (most unstable) wavenumber is complex, arises only for systems with sufficiently large substrate thicknesses. Finally, we discuss adaptation of our model to a practical setting and make predictions of conditions at which the reported instabilities can be observed.
We study the influence of a membrane filter's internal pore structure on its flow and adsorptive fouling behaviour. Membrane performance is measured via (1) comparison between volumetric flow rate ...and throughput during filtration and (2) control of concentration of foulants at membrane pore outlets. Taking both measures into account, we address the merits and drawbacks of selected membrane pore structures. We first model layered planar membrane structures with intra-layer pore connections, and present comparisons between non-connected and connected structures. Our model predicts that membrane filters with connected pore structures lead to higher total volumetric throughput than those with non-connected structures, over the filter lifetime. We also provide a sufficient criterion for the concentration of particles escaping the filter to achieve a maximum in time (indicative of a membrane filter whose particle retention capability can deteriorate). Additionally, we find that the influence of intra-layer heterogeneity in pore-size distribution on filter performance depends on the connectivity properties of the pores.
Thin liquid films in a funnel Lin, T.-S.; Dijksman, J.A.; Kondic, L.
Journal of fluid mechanics,
10/2021, Letnik:
924
Journal Article
Recenzirano
Odprti dostop
We explore flow of a completely wetting fluid in a funnel, with particular focus on contact line instabilities at the fluid front. While the flow in a funnel may be related to a number of other flow ...configurations as limiting cases, understanding its stability is complicated due to the presence of additional azimuthal curvature, as well as due to convergent flow effects imposed by the geometry. The convergent nature of the flow leads to thickening of the film, therefore influencing its stability properties. In this work, we analyse these stability properties by combining physical experiments, asymptotic modelling, self-similar type of analysis and numerical simulations. We show that an appropriate long-wave-based model, supported by the input from experiments, simulations and linear stability analysis that originates from the flow down an incline plane, provides a basic insight allowing an understanding of the development of contact line instability and emerging length scales.
We consider a sheared granular system experiencing intermittent dynamics of stick-slip type via discrete element simulations. The considered setup consists of a two-dimensional system of soft ...frictional particles sandwiched between solid walls, one of which is exposed to a shearing force. The slip events are detected using stochastic state space models applied to various measures describing the system. The amplitudes of the events spread over more than four decades and present two distinctive peaks, one for the microslips and the other for the slips. We show that the measures describing the forces between the particles provide earlier detection of an upcoming slip event than the measures based solely on the wall movement. By comparing the detection times obtained from the considered measures, we observe that a typical slip event starts with a local change in the force network. However, some local changes do not spread globally over the force network. For the changes that become global, we find that their size strongly influences the further behavior of the system. If the size of a global change is large enough, then it triggers a slip event; if it is not, then a much weaker microslip follows. Quantification of the changes in the force network is made possible by formulating clear and precise measures describing their static and dynamic properties.
Recent experiments by Sengupta et al. (Phys. Rev. Lett. 2013) 9 revealed interesting transitions that can occur in flow of nematic liquid crystal under carefully controlled conditions within a long ...microfluidic channel of width much larger than height, and homeotropic anchoring at the walls. At low flow rates the director field of the nematic adopts a configuration that is dominated by the surface anchoring, being nearly parallel to the channel height direction over most of the cross-section; but at high flow rates there is a transition to a flow-dominated state, where the director configuration at the channel centerline is aligned with the flow (perpendicular to the channel height direction). We analyze simple channel-flow solutions to the Leslie–Ericksen model for nematics. We demonstrate that two solutions exist, at all flow rates, but that there is a transition between the elastic free energies of these solutions: the anchoring-dominated solution has the lowest energy at low flow rates, and the flow-dominated solution has lowest energy at high flow rates.
•We present a simple model, and exact solutions, for Poiseuille flow of nematic liquid crystal in a microfluidic.•We demonstrate the existence of two distinct solutions, and argue that these are analogous to two distinct states observed in experiments.•We propose an argument, based on energetics considerations, that explains the observed transitions between these two states.