The double-mass model of drop deformation and secondary breakup

In the paper, a new analytical model of drop deformation and secondary breakup is presented. The model is a direct extension of the TAB (Taylor Analogy Breakup) model of O’Rourke and Amsden [P. O’Rourke, A.A. Amsden, The TAB method for numerical calculation of spray droplet breakup, SAE Paper No 872089, 1987] [9]. The drop is represented by the system of two masses connected by a spring, allowed to oscillate and move along a specified axis. Two versions of the model are analyzed: linear – offering analytical solution for drop oscillations, and nonlinear – defined in terms of parameters with clear, physical interpretation, and more interesting from the point of view of applications. Conditions of stability of a drop subjected to impulsive acceleration by ambient flow are discussed and a new criterion is introduced including droplet Weber number, Ohnesorge number and density ratio. The role of the density ratio proves to be prominent for large Ohnesorge numbers or when the drop density approaches the ambient density.     

固壁近旁Stokes流中粘性液滴的运动和变形

发展了一种模拟固壁近旁轴对称Stokes流中粘性液滴的运动和变形及直接计算固壁上应力的边界积分方法．用此方法对不同的液滴-固壁初始相对间距、粘度比、表面张力和浮力联合参数以及环境流动参数情况进行了数值实验．数值结果显示,由于环境流动和浮力的作用,随着时间的推进,液滴在轴向压缩,在径向拉伸．当环境流动的作用弱于浮力作用时,随着时间的推移,液滴上升并向上弯,固壁上由液滴运动所引起的应力不断减小．当环境流动的作用强于浮力作用时,随着时间的推移,液滴变得越来越扁．在这种情形,当大初始间距时,壁面上的应力随液滴的演变而增大;当小初始间距时,由环境流动、浮力及壁面对流动的较强作用的联合影响,此应力随液滴的演变而减小．由于液滴运动所引起的壁面应力的有效作用仅限于对称轴附近的一个小范围内,且此范围随液滴与固壁的初始间距增大而增大．应力的大小随初始间距增大而大为减小．表面张力对液滴变形有阻止作用．液滴粘性会减小液滴的变形和位置迁移．     

Drop formation from a wettable nozzle

The process of drop formation from a nozzle can be seen in many natural systems and engineering applications. However, previous research focuses on the pinch-off mechanism of drops from a non-wettable nozzle. Here we investigate the formation of a liquid droplet from a wettable nozzle. In the experiments, drops forming from a wettable nozzle initially climb the outer walls of the nozzle due to surface tension. Then, when the weight of the drops gradually increases, they eventually fall due to gravity. By changing the parameters such as the nozzle size and fluid flow rate, we have observed different behaviors of the droplets. Such oscillatory behavior is characterized by an equation that consists of capillary force, viscous drag, and gravity. Two asymptotic solutions in the initial and later stages of drop formation are obtained and show good agreement with the experimental observations.     

Liquid drops in a viscous fluid under the influence of gravity and surface tension

We consider the steady fall of a drop of fluid under its own weight in an infinite reservoir of another viscous fluid; the shape of the drop is determined by surface tension. For small data we prove existence and uniqueness of a classical solution to this problem.     

Deformation of a drop in a viscous stream and conditions of existence of its equilibrium form

The difference method is used for obtaining a solution of the problem of unsteady motion of a drop in a stream, taking into account its deformation under conditions of axial symmetry. The fluid inside and outside the drop is assumed viscous and incompressible. The stable forms of drop are represented for various Reynolds and Weber numbers of external stream. By analyzing the conditions for normal stresses at the drop boundary, the critical Weber number was obtained, which establishes the conditions of existence of equilibrium form of the drop.     

The stationary motion of a weakly deformed drop in the flow of a viscous fluid

The flow of a drop of viscous, incompressible fluid in another viscous, incompressible fluid is studied with values of the viscosities, densities, surface tension, and drop size which guarantee small deviation of the shape of the drop from spherical. The effect of the deviation on the character of the flow, the drag coefficient, etc. is considered. The cell model for the motion of a weakly deformed drop is investigated. An approximate formula is given for finding the drag coefficient in constrained flow. Some cases of nonisothermic motion of a spherical drop are studied.Translated from Zapiski Nauchnykh Seminarov Leningradskogo Otdeleniya Matematicheskogo Instituta im. V. A. Steklova AN SSSR, Vol. 69, pp. 157–170, 1977.     

On the Deformation of a Viscous Droplet Caused by Variable Surface Tension

The surface tension of a small viscous droplet is made variable by the deposition of a surface active material. The resultant shear forces produce a deformation of and a fluid motion within the droplet. The process, including the nonlinear convection or spread of surfactant, is examined theoretically. Some experimental results are reported, and application to cell division is discussed.     

Numerical simulation of a flexible fiber deformation in a viscous flow by the immersed boundary-lattice Boltzmann method

In this paper, deformation of a mass-less elastic fiber with a fixed end, immersed in a two-dimensional viscous channel flow, is simulated numerically. The lattice-Boltzmann method (LBM) is used to solve the Newtonian flow field and the immersed-boundary method (IBM) is employed to simulate the deformation of the flexible fiber interacting with the flow. The results of this unsteady simulation including fiber deformation, fluid velocity field, and variations of the fiber length are depicted in different time-steps through the simulation time. Similar trends are observed in plots representing length change of fibers with different values of stretching constant. Also, the numerical solution reaches a steady state equivalent to the fluid channel flow over a flat plate.     

A numerical procedure for viscous free surface flows

We present a numerical procedure for two-dimensional unsteady viscous free surface flow problems with surface tension. The procedure is based on a finite difference approach to a primitive variable formulation; a coordinate transformation is used to transform the irregularly shaped flow domain onto a fixed rectangular domain. The procedure is tested on a standing wave problem and a cavity flow problem with a free surface. Satisfactory numerical solutions are obtained for both problems for Reynolds numbers up to 200.     

Axisymmetric spreading of a thin drop under gravity and time-dependent non-uniform surface tension

A second-order non-linear partial differential equation modelling the gravity driven spreading of a thin viscous liquid film with time-dependent non-uniform surface tension Σ(t,r) is considered. The problem is specified in cylindrical polar coordinates where we assume the flow is axisymmetric. Similarity solutions describing the spreading of a thin drop and the flattening of a thin bubble are investigated.     

Effect of viscosity of a thermoplastic prepreg and matrix upon winding of rings

The problem of compression of a unidirectional layer and shear of a polymer interlayer during winding of rings is considered. The equations determining the dependence of the layer thickness and stresses on the parameters entering into the power flow law for a prepreg and polymer matrix and on the basic parameters of the winding process—the initial tension of the prepreg, its placement rate, and the radius of a mandrel—are derived. The ring thickness measurements obtained at various temperatures and initial tension forces of plies confirm the adequacy of the model offered. It is found that the viscous properties of the prepreg and matrix upon winding affect the relative change in the layer thickness to a greater extent than the stresses in these layers. With increase in temperature and tension force upon winding, the effect of viscous deformations of the prepreg and matrix increases. A decrease in viscosity and an increase in the tension force of the tape lead to a higher strength of the ring in tension and interlaminar shear; however, the growing percolation of the polymer melt leads to a greater inhomogeneity of the structure of the composite in the ring and to a lower reinforcing effect of the factors mentioned. Presented at the 11th International Conference on Mechanics of Composite Materials (Riga, June 11–15, 2000). Translated from Mekhanika Kompozitnykh Materialov, Vol. 36, No. 3, pp. 419–428, 2000.     

Theoretical analysis of Rayleigh–Taylor instability on a spherical droplet in a gas stream

A linear analysis of the Rayleigh–Taylor (R–T) instability on a spherical viscous liquid droplet in a gas stream is presented. Different from the most previous studies in which the external acceleration is usually assumed to be radial, the present study considers a unidirectional acceleration acting on a spherical droplet with arbitrary initial disturbances and therefore can provide insights into the influence of R–T instability on the atomization of spherical droplets. A general recursion relation coupling different spherical modes is derived and two physically prevalent limiting cases are discussed. In the limiting case of inviscid droplet, the critical Bond numbers to excite the instability and the growth rates for a given Bond number are obtained by solving two eigenvalue problems. In the limiting case of large droplet acceleration, different spherical modes are asymptotically decoupled and an explicit dispersion relation is derived. For given Bond number and Ohnesorge numbers, the critical size of stable droplet, the most-unstable mode and its corresponding growth rate are determined theoretically.     

液滴的脉动雾化

应用稳定性分析,本文显示一球形液滴可以由参数共振而促使形成比原液滴直径小得多的小微滴。表面张力和周围气体的密度将延缓这一雾化过程。与达朗伯表观体力有关的两重迭平面流体层的界面不稳定性与参数共振引起的界面不稳定性有着密切的关系。     

Free boundary fluid flow in a trough: A numerical approach

This paper examines the slow flow of a viscous liquid in an open rectangular container, one side (the base) of which moves steadily along its own plane, thereby providing the driving force the liquid needs. Unlike the two vertical sides that are rigid and stationary, the top side is left open so that the upper part of the liquid is in contact with air and is being controlled by surface tension and gravity. A numerical procedure for obtaining solutions for the cases when the capillary numbers are small is provided and the curves of the free boundaries obtained here are presented for some flow parameters. The deviation of the shape of the free boundary is observed to be strongly dependent on the aspect ratio of the boundary (i.e., the ratio of the vertical to horizontal spread of the liquid) with its curvature changing sign in the interval [1, 1.5].     

An immersed boundary-lattice Boltzmann method combined with a robust lattice spring model for solving flow–structure interaction problems

An immersed boundary (IB)-lattice Boltzmann method (LBM) combined with a robust lattice spring model (LSM) was developed for modeling fluid–elastic body interactions. To include the effects of viscous flow forces on the deformation of a flexible body, rotational invariant springs were connected regularly inside the deformable body with square lattices. Fluid–solid interactions were due to an additional force density in the lattice Boltzmann equation enhanced by the split-forcing approach. To check the validity and accuracy of the numerical method, the flow over a rigid plate and the deformation of a cantilever beam were investigated. To demonstrate the capability of the new method, different test cases were examined. The deformation of a two-dimensional flexible vertical plate in a laminar cross-flow stream at different conditions was analyzed. The simulations were performed for different boundary conditions imposed on the elastic plate, namely, fixed-end corners and fixed middle point. Different flow conditions such as “steady flow regime”, “vortex shedding flow regime”, and the limit of “rigid body motion” were examined using the new IB-LBM-LSM approach. A general formulation for evaluating the deformation of the elastic body was also introduced, in which the position of the LSM nodes (inside the body) was updated implicitly at each time step. Two dimensionless groups, namely capillary number (Ca) and Reynolds number (Re), were used for parametric study of the behavior of the flow around the deformable plate. It was found that for low Reynolds numbers (Re < 50) and when the middle of the plate was fixed, decreasing the capillary number led to a decrease in the drag coefficient. The fluctuation of the plate during the vortex shedding flow regime was also explored. It was found that when the middle of the plate was fixed, the critical Reynolds number for the initiation of vortex shedding increased. For Re > 100, the Strouhal number was observed to increase with the decrease in capillary number.     

Dynamics of surfactant-covered drops in linear flows

Surfactants modify interfacial properties and significantly affect drop behavior in flow. We study the dynamics of a drop, which is covered with a monolayer of insoluble surfactant, in linear viscous flows, both unbounded and in the presence of a wall. The effect of viscosity contrast is included. Analytical results in a form of expansions for weak flows and high viscosity drops are developed. Numerical results with three-dimensional boundary integral simulations are used to explore large deformations. The results show that surfactant generally enhances drop deformation, certainly under small-deformation conditions. The steady-state drop shape and surfactant distribution are independent of viscosity contrast in straining flows (e.g. hyperbolic, axisymmetric strain). The drop shape and surfactant distribution are insensitive to viscosity contrast under small deformation conditions for any linear flow; the effect of the rotational component at higher-order. The theory quantifies the effect of surfactant on drop migration away from a bounding wall as well as the rheology of dilute emulsions. It predicts non-Newtonian features such as shear thinning viscosity and normal stress. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Spatial Velocity and Temperature Measurements in Non‐Isothermal Liquid Flows

A surface tension driven flow in the liquid vicinity of an air bubble on a heated wall is studied experimentally. The liquid flow caused by the temperature gradient along the surface of the bubble is termed thermocapillary convection. The surface tension force and the buoyancy force oppose one another. The measurement technique is the 3D particle tracking velocity and thermometry, 3D PTV/T, using thermochromic liquid crystals and digital image processing. The paper describes the method in some detail and presents quantitative results for different Marangoni numbers. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The effect of surface elasticity and residual surface stress in an elastic body with an elliptic nanohole

The deformation of an elastic plane with an elliptic hole in a uniform stress field is considered, taking into account the surface elasticity and the residual surface tension. The solution of the problem, based on the use of the linearized Gurtin–Murdoch surface elasticity relations and the complex Goursat–Kolosov potentials, is reduced to a singular integrodifferential equation. Using the example of a circular hole, for which an exact solution of the equation is obtained in closed form, the effect of the residual surface tension and the surface elasticity on the stress state close to and on the boundary of a nanohole is analysed for uniaxial tension. It is shown that the effect of the residual surface stress and the surface tension, due to deformation of the body, depends on the elastic properties of the surface, the value of the stretching load and the dimensions of the hole.     

Thermal and surface effects on the pull-in characteristics of circular nanoplate NEMS actuator based on nonlocal elasticity theory

This paper aims to investigate the coupling influences of thermal loading and surface effects on pull-in instability of electrically actuated circular nanoplate based on Eringen's nonlocal elasticity theory, where the electrostatic force and thermally corrected Casimir force are considered. By utilizing the Kirchhoff plate theory, the nonlinear equilibrium equation of axisymmetric circular nanoplate with variable coefficients and clamped boundary conditions is derived and analytically solved. The results describe the influences of surface effect and thermal loading on pull-in displacements and pull-in voltages of nanoplate under thermal corrected Casimir force. It is seen that the surface effect becomes significant at the pull-in state with the decrease of nanoplate thicknesses, and the residual surface tension exerts a greater influence on the pull-in behavior compared to the surface elastic modulus. In addition, it is found that temperature change plays a great role in the pull-in phenomenon; when the temperature change grows, the circular nanoplate without applied voltage is also led to collapse.     

Effect of magnetic field on the liquid-gas interfacial deformation in thermocapillary convection problem

Liquid contained in an open cavity flows owing to a temperature gradient applied along its free surface. The thermal variation of surface tension induces a steady viscous flow directed from hot end to cold end. For small aspect ratios, giving flow in thin, two dimensional slots, an asymptotic theory valid for (A ? 0)(A \rightarrow 0) is used to show that the interface undergoes an O(1) deformation from its flat position. Further, it is observed that for some values of the Hartmann number, the deformed interface can be made almost flat.