Shear cell rupture of nematic liquid crystal droplets in viscous fluids

We model the hydrodynamics of a shear cell experiment with an immiscible nematic liquid crystal droplet in a viscous fluid using an energetic variational approach and phase-field methods [86]. The model includes the coupled system for the flow field for each phase, a phase-field function for the diffuse interface and the orientational director field of the liquid crystal phase. An efficient numerical scheme is implemented for the two-dimensional evolution of the shear cell experiment for this initial data. The same model reduces to an immiscible viscous droplet in a viscous fluid, which we simulate first to compare with other numerical and experimental behavior. Then we simulate drop deformation by varying capillary number (independent of liquid crystal physics), liquid crystal interfacial anchoring energy and Oseen–Frank distortional elastic energy. We show the number of eventual droplets (one to several) and “beads on a string” behavior are tunable with these three physical parameters. All stable droplets possess signature quadrupolar shear and normal stress distributions. The liquid crystal droplets always possess a global surface defect structure, called a boojum, when tangential surface anchoring is imposed. Boojums [79], [32] consist of degree +1/2 and ?1/2 surface defects within a bipolar global orientational structure.     

Droplet relaxation in blends with one viscoelastic component: bulk and confined conditions

Using a counter rotating parallel plate shear flow cell, the shape relaxation of deformed droplets in a quiescent matrix is studied microscopically. Both the effects of geometrical confinement and component viscoelasticity are systematically explored at viscosity ratios of 0.45 and 1.5. The flow conditions are varied from a rather low to a nearly critical Ca number. Under all conditions investigated, viscoelasticity of the droplet phase has no influence on shape relaxation, whereas matrix viscoelasticity and geometrical confinement result in a slower droplet retraction. Up to high confinement ratios, the relaxation curves for ellipsoidal droplets can be superposed onto a master curve. Confined droplets with a sigmoidal shape relax in two stages: the first consists of a shape change to an ellipsoid with a limited amount of retraction, and the second is the retraction of this ellipsoid. The latter stage can be described by means of one single relaxation time that can be obtained from the relaxation of initially ellipsoidal droplets. The experimental results are compared to the predictions of a recently published phenomenological model for droplet dynamics in confined systems with viscoelastic components (Minale et al., Langmuir 26:126–132, 2010). However, whereas the model predicts additive effects of geometrical confinement and component viscoelasticity, the experimental data reveal more complex interactions.     

Droplet retraction in the presence of nanoparticles with different surface modifications

We studied the influence of nanoparticles with different surface modifications on the interfacial tension and relaxation of model polymer blend after cessation of different strains. The droplet retraction experiments were carried out on a model system composed of polydimethylsiloxane (PDMS) as the suspending fluid and polyisobutylene (PIB) as droplet at room temperature in the presence of hydrophobic and hydrophilic nanosilica. Different weight fractions of particles were dispersed in the PIB droplet before forming a dispersed droplet by using a microsyringe in shear cell. We found that applied strain, nanoparticle concentration and their thermodynamically preferred localization affect both nominal interfacial tension and droplet retraction process. By addition of nanoparticles at a concentration as low as 0.2%wt, the nominal interfacial tension decreases from 3.12?±?0.15 mN/m for neat PIB-PDMS interface depending on the surface characteristics of nanosilica. Hydrophilic nanosilica has the most effect on nominal interfacial tension and decreases it as low as 0.2?±?0.21 mN/m at 1 wt.% loading under a strain of 7. The results show that the retraction process in this system is mainly controlled by interfacial phenomena rather than bulk rheological properties. Additionally, the shape evolution of droplets changes and the retraction rate slows down in the presence of nanoparticles.     

Shape and structure recovery of an LCP droplet under a large step strain: observation and stress calculation

Shape recovery of a droplet of liquid crystalline polymer (LCP) hydroxypropylcellulose in a matrix of poly(dimethyl siloxane) subjected to a step shear strain has been studied via optical microscopy. Just after application of a large strain, the LCP droplet shape is flat ellipsoid, and then the droplet takes cylindrical shape and band texture perpendicular to the flow direction appears. The band texture fades away before emergence of poly-domain structure. In the final process with the shape of spheroid, poly-domain structure recovers very slowly. Except for the final process, the shape change is identical with that of isotropic droplet at strains smaller than 3, when the LCP viscosity in Region II is taken as an equivalent viscosity for normalization. For a 20:80 blend, the excess relaxation modulus is calculated based on the Doi-Ohta theory, taking account of the distribution of droplet size and compared with experimental modulus data.     

A parametric study of droplet deformation through a microfluidic contraction: Shear thinning liquids

Numerical simulations of a droplet passing through an axisymmetric microfluidic contraction are presented, focusing on systems where one of the two liquids present is shear thinning. The simulations are performed using a transient Volume of Fluid (VOF) algorithm. When the droplet is shear thinning and the surrounding phase Newtonian, droplets deform in a similar way to Newtonian droplets that have a viscosity equal to the average viscosity of the shear thinning fluid while it is within the contraction. When the surrounding phase is shear thinning and the droplet Newtonian, droplets deform in a similar way to droplets contained within a Newtonian liquid that has a viscosity that is lower than that of the droplet. In both cases the behaviour of the shear thinning fluid can be broadly described in terms of a ‘characteristic’ Newtonian viscosity: However, determining the exact value of this viscosity without performing a full shear thinning simulation is not possible.     

基于SPH方法的湿润性固壁边界条件模拟研究

固体壁面由于表面特殊结构和材料属性,时常表现出对交界面上水体的吸附作用,而这一特征对微小水体作用尤为明显。本文提出了一种湿润性固壁边界条件的计算方法,即假设壁面粒子的亲水性以及毛细吸附作用统一表现为对支持域内流体粒子的吸附力。基于光滑粒子流体动力学（SPH）方法,模拟了静态液滴在不同湿润性壁面上的变形至稳定过程。模拟了液滴撞击疏水壁面的过程,将液滴的运动过程分为碰撞、铺展、回缩和回弹四个阶段,分析各阶段壁面受力分布情况。研究表明：根据模拟液滴静态接触角的变化特点,本文湿润性固壁边界条件可以较好的反映出壁面湿润性;液滴撞击输水表面的模拟数据与试验结果趋势上吻合良好;壁面压力波伴随着液滴的铺展和回缩传播并衰减;只有在回弹后期液滴即将脱离壁面时壁面拉力起主导作用,其余各时刻壁面均以压力为主。     

基于LS-DYNA的热成型钢断裂失效预测研究

固体壁面由于表面特殊结构和材料属性,时常表现出对交界面上水体的吸附作用,而这一特征对微小水体作用尤为明显。本文提出了一种湿润性固壁边界条件的计算方法,即假设壁面粒子的亲水性以及毛细吸附作用统一表现为对支持域内流体粒子的吸附力。基于光滑粒子流体动力学（SPH）方法,模拟了静态液滴在不同湿润性壁面上的变形至稳定过程。模拟了液滴撞击疏水壁面的过程,将液滴的运动过程分为碰撞、铺展、回缩和回弹四个阶段,分析各阶段壁面受力分布情况。研究表明：根据模拟液滴静态接触角的变化特点,本文湿润性固壁边界条件可以较好的反映出壁面湿润性;液滴撞击输水表面的模拟数据与试验结果趋势上吻合良好;壁面压力波伴随着液滴的铺展和回缩传播并衰减;只有在回弹后期液滴即将脱离壁面时壁面拉力起主导作用,其余各时刻壁面均以压力为主。     

Modelling impingement of hollow metal droplets onto a flat surface

Despite many theoretical and experimental works dealing with the impact of dense melt droplets on the substrate during the process of thermal spray coating, the dynamics of the impingement of hollow melt droplet and the subsequent splat formation are not well addressed. In this paper a model study for the dynamic impingement of hollow droplet is presented. The hollow droplet is modelled such that it consists of a liquid shell enclosing a gas cavity. The impingement model considers the transient flow dynamics during impact, spreading and solidification of the droplet using the volume of fluid surface tracking method (VOF) coupled with a solidification model within a one-domain continuum formulation. The results for spreading, solidification and formation of splats clearly show that the impingement process of hollow droplet is distinctly different from the dense droplet. Study with different droplet void fractions and void distribution indicates that void fraction and void distribution have a significant influence on the flow dynamics during impact and on the final splat shape. The results are likely to provide insights for the less-explored behaviour of hollow melt droplets in thermal spray coating processes.     

镓铟锡液滴撞击泡沫金属表面的运动学特性研究

常温下为液态的镓铟锡合金以其优异的导热性能在具有特殊要求的传热领域有着重要的应用价值,与传统流动介质相比较大的表面张力使得其产生的流动现象必有所区别.本文研究镓铟锡所形成的液滴撞击泡沫金属表面后所产生的铺展、回缩及回弹现象.采用高速相机拍摄液滴投影轮廓随液滴运动的变化过程,并通过图像处理获得不同撞击速度、底板表面孔径下的液滴铺展系数、中心位置轮廓高度以及液滴回弹后在空中的振动特性.研究结果表明:具有较高表面张力的镓铟锡液滴的铺展系数随无量纲时间的变化在铺展初始阶段仍满足常规流体的1/2次幂关系,只在铺展后期与底板的无量纲孔径有关系;液滴的最大铺展系数在较小无量纲孔径底板大于在光滑镍板,且随底板无量纲孔径增大而逐渐减小;在回弹过程,由于底板孔隙结构的存在使得液滴回弹后在空中的振动呈现3种形态:规则的横向和纵向振动、带旋转的横向和纵向振动以及旋转振动;最后,通过对振动频率的拟合和分析,进一步拓展了传统振动频率理论公式在非规则振动过程预测中的应用.      

Role of droplet bridging on the stability of particle-containing immiscible polymer blends

The effect of micron-sized hydrophobic calcium carbonate particles on the stabilization of polydimethylsiloxane (PDMS)/polyisobutylene (PIB) immiscible model blends is investigated in this study. The analytical splitting of bulk and liquid–liquid interface contributions from the droplet bridging one is successfully performed due to the negligible contribution of hydrophobic microparticles to the bulk rheology of phases. The presence of particles at the fluid–fluid interface is supported by wetting parameter calculation and verified by optical microscopy observations. Moreover, direct visualizations shows that particles are able to form clusters of droplets by simultaneously adsorbing at two fluid–fluid interfaces and glue-dispersed droplets together, probably due to the patchy interactions induced by heterogeneous distribution of particles along the interface. Rheological studies show that the flow-induced coalescence is slowed down upon addition of particles and almost suppressed with the addition of 4 wt% particles. The linear viscoelastic response is modeled to estimate interfacial tension by considering the contribution of particle-induced droplet aggregation in addition to bulk and droplet deformation ones. From linear and nonlinear viscoelastic responses, the improved stability of filled polymer blends is attributed to the interfacial rheology and/or the bridged structure of droplets, even though the interfacial area is not fully covered by particles. Furthermore, Doi–Ohta scaling relations are investigated by employing stress growth response upon step-up of shear flow.     

Drop shape under slow steady shear flow and during relaxation. Experimental results and comparison with theory

The three-dimensional small deformation of a single Newtonian drop immersed in an immiscible Newtonian liquid was investigated, both in slow steady shear and during retraction after cessation of shear. The experiments were performed in a parallel plate apparatus equipped with video- enhanced microscopy. The drop was observed from two perpendicular directions, and accurate measurements were obtained in each view by image analysis. The results were compared to existing theoretical predictions from the exact fluid dynamic problem, obtained perturbatively for the case of small deformations of the drop. Excellent agreement between data and theory was found, thus providing the first assessment of drop shape predictions in steady flow and relaxation for Newtonian fluids. Received: 6 February 2000 Accepted: 24 May 2000     

Viscosity of a Newtonian fluid calculated from the deformation of droplets covered with a surfactant under a linear shear flow

The viscosity of small fluid droplets covered with a surfactant is determined using drop deformation techniques. This method, proposed by Hu and Lips, is here extended to the case of the presence of a surface-active adsorpted at the liquid–liquid interface, to consider more general scenarios. In these experiments, a droplet is sheared by another immiscible fluid of known viscosity, both Newtonian liquids. From the steady-state deformation and retraction mechanisms, the droplet viscosity is calculated using an equation derived from the theories of Taylor and Rallison. Although these theories were expressed for surfactant-free interfaces, they can be applied when a surfactant is present in the system if the sheared droplet reaches reliable steady-state deformations and the surfactant attains its equilibrium adsorption concentration. These determinations are compared to bulk viscosities measured in a rheometer for systems with different viscosity ratios and surfactant concentrations. Very good agreement between both determinations is found for drops more viscous than the continuous phase.     

Electrokinetic Mixing and Displacement of Charged Droplets in Hydrogels

Mixing in droplets is an essential task in a variety of microfluidic systems. Inspired by electrokinetic mixing, electric field-induced hydrodynamic flow inside a charged droplet embedded in an unbounded polyelectrolyte hydrogel is investigated theoretically. In this study, the polyelectrolyte hydrogel is modeled as a soft, and electrically charged porous solid saturated with a salted Newtonian fluid, and the droplet is considered an incompressible Newtonian fluid. The droplet-hydrogel interface is modeled as a surface, which is located at the plane of shear, with the electrostatic potential \(\zeta \) . The fluid inside the droplet attains a finite velocity owing to hydrodynamic coupling with the electroosmotic flow arising from the droplet and polymer charge. The fluid velocity inside the droplet is linearly proportional to the electroosmotic flow velocity in the charged gel and the electroosmotic flow velocity beyond the electrical double layer of a charged interface. It is found that the polymer boundary condition at the droplet surface and the viscosities of the fluids inside and outside the droplet significantly modulate the interior fluid flow. The ionic strength and the permeability of the polymer network impact the flow differently depending on whether the flow arises from the droplet or polymer charge. Finally, the displacement of a charged droplet embedded in a gel under the influence of an external electric field is undertaken. This work is motivated by experimental attempts, which can register sub-nanometer-scale inclusion displacements in hydrogels, to advance electrical microrheology as a diagnostic tool for probing inclusion-hydrogel interfaces. In the absence of polymer charge, a close connection is found between the electrical response of a charged droplet when it is immobilized in an uncharged incompressible gel and when it is dispersed in a Newtonian electrolyte.     

Investigation the effects of shear rate on stationary droplets coalescence by lattice Boltzmann

Coalescence is the process by which two or more droplets merge to form a single droplet. It can take place in many processes, ranging from meteorology to astrophysics. When two stationary droplets are suspended in a bulk vapor, if the gap between the droplets has been smaller than a critical value, the two droplets will coalesce. In this paper, single component, two phase flow is modeled under shear flow using a free energy lattice Boltzmann approach and the coalescence of stationary droplets are investigated for different cases of radius and shear rate. The results show that there is a critical gap between droplets and for the values larger than that they will not coalesce. Also in the case of constant thermophysical properties, this critical gap is a function of droplet radius and shear rate.     

The translation of a Newtonian droplet in a 4-constant Oldroyd fluid

The viscoelastic correction to the Hadamard—Rybczynski terminal velocity was calculated for the retarded motion translation of a finite viscosity Newtonian droplets in a 4-constant Oldroyd fluid. The droplet translational velocity can be enhanced or hindered relative to the Hadamard—Rybczinski value depending upon the degree of shear thinning, and elongational and memory effects in the viscoelastic fluid. Specifically, it is found that for a large droplet viscosity, for which surface mobility is negligible, shear thinning and fluid memory enhance translational motion. At the other extreme, when the droplet surface is very mobile (i.e. for gas bubbles) the translational velocity can be enhanced or hindered relative to the Hadamard—Rybczynski value depending on the relative influence of the memory and the elongational properties of the viscoelastic phase.     

An icing physics study by using lifetime-based molecular tagging thermometry technique

An experimental investigation was conducted to quantify the unsteady heat transfer and phase changing process within small icing water droplets in order to elucidate underlying physics to improve our understanding of the important micro-physical process of icing phenomena. A novel, lifetime-based molecular tagging thermometry (MTT) technique was developed and implemented to achieve temporally-and-spatially resolved temperature distribution measurements to reveal the time evolution of the unsteady heat transfer and dynamic phase changing process within micro-sized water droplets in the course of icing process. It was found that, after a water droplet impinged onto a frozen cold surface, the liquid water at the bottom of the droplet would be frozen and turned to solid ice rapidly, while the upper portion of the droplet was still in liquid state. As the time goes by, the interface between the liquid phase water and solid phase ice was found to move upward continuously with more and more liquid water within the droplet turned to solid ice. Interestingly, the averaged temperature of the remaining liquid water within the small icing droplet was found to increase, rather than decrease, continuously in the course of icing process. The temperature increase of the remaining liquid water is believed to be due to the heat release of the latent heat during solidification process. The volume expansion of the water droplet during the icing process was found to be mainly upward to cause droplet height growth rather than radial to enlarge the contact area of the droplet on the test plate. As a result, the spherical-cap-shaped water droplet was found to turn to a prolate-spheroid-shaped ice crystal with cusp-like top at the end of the icing process. The required freezing time for the water droplets to turn to ice crystals completely was found to depend on the surface temperature of the test plate strongly, which would decrease exponentially as the surface temperature of the frozen cold test plate decreases.     

A modified SPH method to model the coalescence of colliding non-Newtonian liquid droplets

This paper develops a modified smoothed particle hydrodynamics (SPH) method to model the coalescence of colliding non-Newtonian liquid droplets. In the present SPH, a van der Waals (vdW) equation of state is particularly used to represent the gas-to-liquid phase transition similar to that of a real fluid. To remove the unphysical behavior of the particle clustering, also known as tensile instability, an optimized particle shifting technique is implemented in the simulations. To validate the numerical method, the formation of a Newtonian vdW droplet is first tested, and it clearly demonstrates that the tensile instability can be effectively removed. The method is then extended to simulate the head-on binary collision of vdW liquid droplets. Both Newtonian and non-Newtonian fluid flows are considered. The effect of Reynolds number on the coalescence process of droplets is analyzed. It is observed that the time up to the completion of the first oscillation period does not always increase as the Reynolds number increases. Results for the off-center binary collision of non-Newtonian vdW liquid droplets are lastly presented. All the results enrich the simulations of the droplet dynamics and deepen understandings of flow physics. Also, the present SPH is able to model the coalescence of colliding non-Newtonian liquid droplets without tensile instability.     

Convective heat transfer from molten salt droplets in a direct contact heat exchanger

This paper presents a new predictive model of droplet flow and heat transfer from molten salt droplets in a direct contact heat exchanger. The process is designed to recover heat from molten CuCl in a thermochemical copper–chlorine (Cu–Cl) cycle of hydrogen production. This heat recovery occurs through the physical interaction between high temperature CuCl droplets and air. Convective heat transfer between droplets and air is analyzed in a counter-current spray flow heat exchanger. Numerical results for the variations of temperature, velocity and heat transfer rate are presented for two cases of CuCl flow. The optimal dimensions of the heat exchanger are found to be a diameter of 0.13 m, with a height of 0.6 and 0.8 m, for 1 and 0.5 mm droplet diameters, respectively. Additional results are presented and discussed for the heat transfer effectiveness and droplet solidification during heat recovery from the molten CuCl droplets.     

Droplet-based flows in serpentine microchannels: Chemical reactions and secondary flows

Mixing is an essential operation in many microfluidic devices. Droplet-based micromixers utilize droplets for mixing enhancement. In the present study, a novel three-dimensional simulation is conducted which has the ability to capture not only the mixing process, but also the chemical reactions inside liquid droplets. This two-phase model is used for simulating the reacting flow inside a serpentine microchannel and explores the effects of droplet size and reaction rate on the production and consumption of species in droplets. It is observed that the chemical reaction in each droplet, begins from its front area. Furthermore, it is shown that the production of species does not depend on water fraction (ratio of water flow rate to total flow rate) and for all droplet sizes, only depends on the reaction rate. Moreover, different transient generated vortices and secondary flows are studied in the presence and absence of droplets and explained in details based on the position of droplets. It is found that during the passage of droplets through a microchannel bend, the number of cross-sectional vortices changes and the velocity magnitude in these secondary flows increases dramatically.