Numerical simulations of three-dimensional drop collisions

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/77346/1/AIAA-13136-185.pd

The Flow Induced by the Coalescence of Two Initially Stationary Drops

The coalescence of two initially stationary drops of different size is investigated by solving the unsteady, axisymmetric Navier-Stokes equations numerically, using a Front-Tracking/Finite Difference method. Initially, the drops are put next to each other and the film between them ruptured. Due to surface tension forces, the drops coalesce rapidly and the fluid from the small drop is injected into the larger one. For low nondimensional viscosity, or Ohnesorge number, little mixing takes place and the small drop fluid forms a blob near the point where the drops touched initially. For low Ohnesorge number, on the other hand, the small drop forms a jet that penetrates far into the large drop. The penetration depth also depends on the size of the drops and shows that for a given fluid of sufficiently low viscosity, there is a maximum penetration depth for intermediate size ratios

Numerical Simulations of Drop Collisions

Three-dimensional simulations of the off-axis collisions of two drops are presented. The full Navier-Stokes equations are solved by a Front-Tracking/Finite-Difference method that allows a fully deformable fluid interface and the inclusion of surface tension. The drops are accelerated towards each other by a body force that is turned off before the drops collide. Depending on whether the interface between the drops is ruptured or not, the drops either bounce or coalesce. For drops that coalesce, the impact parameter, which measures how far the drops are off the symmetry line, determines the eventual outcome of the collision. For low impact parameters, the drops coalesce permanently, but for higher impact parameters, a grazing collision, where the drops coalesce and then stretch apart again is observed. The results are in agreement with experimental observations

Numerical simulations of drop collisions

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76145/1/AIAA-1994-835-900.pd

Head‐on collision of drops—A numerical investigation

The head‐on collision of equal sized drops is studied by full numerical simulations. The Navier–Stokes equations are solved for the fluid motion both inside and outside the drops using a front tracking/finite difference technique. The drops are accelerated toward each other by a body force that is turned off before the drops collide. When the drops collide, the fluid between them is pushed outward leaving a thin layer bounded by the drop surface. This layer gets progressively thinner as the drops continue to deform, and in several of our calculations we artificially remove this double layer at prescribed times, thus modeling rupture. If no rupture takes place, the drops always rebound, but if the film is ruptured the drops may coalesce permanently or coalesce temporarily and then split again. Although the numerically predicted boundaries between permanent and temporary coalescence are found to be consistent with experimental observations, the exact location of these boundaries in parameter space is found to depend on the time of rupture. © 1996 American Institute of Physics.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/71337/2/PHFLE6-8-1-29-1.pd

Head-on collision of drops: A numerical investigation

The head-on collision of equal sized drops is studied by full numerical simulations. The Navier-Stokes equations are solved for fluid motion both inside and outside the drops using a front tracking/finite difference technique. The drops are accelerated toward each other by a body force that is turned off before the drops collide. When the drops collide, the fluid between them is pushed outward leaving a thin later bounded by the drop surface. This layer gets progressively thinner as the drops continue to deform and in several of the calculations this double layer is artificially removed once it is thin enough, thus modeling rupture. If no rupture takes place, the drops always rebound, but if the film is ruptured the drops may coalesce permanently or coalesce temporarily and then split again

Application of Electric Field to Developing Falling Films using Wire-Plate Electrode Configuration- An Experimental Study

Experimental investigation of Electrohydrodynamic developing falling film flow of transformer oil has been conducted within an inclined rectangular channel and hydrodynamic characteristics of the flow have been revealed. The electric field has been generated by five overhead thin wire electrodes connected to the positive high DC voltage on the air and the grounded plate electrodes which are placed upon the floor of the channel. It is the first time that the wavy behavior on a liquid falling film's interface has been created by this electrode configuration. A non-intrusive method has been used to measure the local flow structure by a high-speed camera, then statistical characteristics of the wavy falling film have been computed by image processing of the captured video frames. By applying 13-16 kV to the wire electrodes, the influence of EHD force on the wavy behavior of falling film has been conducted for Reynolds number 10-120 in the laminar-wavy regime at three different inclination angles 15◦, 30◦ and 45°. The vertical distance of the high-voltage wire electrodes to ground electrodes has been set to 14 mm. The liquid velocity, film thickness, and wave frequency have been measured for non-electrified and electrified falling film, and their results have been evaluated with other experimental studies and an acceptable agreement has been obtained. The results indicate that the proposed HV wire-grounded plate electrode configuration in this study does not disturb the original structure of the falling film and by intensifying the wavy behavior of laminar falling film can either suppress or enhance heat/mass transfer rate. The effects of the applied voltage on the frequency, velocity and film thickness of the falling liquid film have been also discussed in detail

Fluid-structure interaction simulation of prosthetic aortic valves : comparison between immersed boundary and arbitrary Lagrangian-Eulerian techniques for the mesh representation

In recent years the role of FSI (fluid-structure interaction) simulations in the analysis of the fluid-mechanics of heart valves is becoming more and more important, being able to capture the interaction between the blood and both the surrounding biological tissues and the valve itself. When setting up an FSI simulation, several choices have to be made to select the most suitable approach for the case of interest: in particular, to simulate flexible leaflet cardiac valves, the type of discretization of the fluid domain is crucial, which can be described with an ALE (Arbitrary Lagrangian-Eulerian) or an Eulerian formulation. The majority of the reported 3D heart valve FSI simulations are performed with the Eulerian formulation, allowing for large deformations of the domains without compromising the quality of the fluid grid. Nevertheless, it is known that the ALE-FSI approach guarantees more accurate results at the interface between the solid and the fluid. The goal of this paper is to describe the same aortic valve model in the two cases, comparing the performances of an ALE-based FSI solution and an Eulerian-based FSI approach. After a first simplified 2D case, the aortic geometry was considered in a full 3D set-up. The model was kept as similar as possible in the two settings, to better compare the simulations' outcomes. Although for the 2D case the differences were unsubstantial, in our experience the performance of a full 3D ALE-FSI simulation was significantly limited by the technical problems and requirements inherent to the ALE formulation, mainly related to the mesh motion and deformation of the fluid domain. As a secondary outcome of this work, it is important to point out that the choice of the solver also influenced the reliability of the final results