Direct numerical simulations of flows with phase change

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76607/1/AIAA-1996-857-548.pd

An experimental and computational study of bouncing and deformation in droplet collision

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/76090/1/AIAA-1997-129-419.pd

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

A momentum-conserving, consistent, Volume-of-Fluid method for incompressible flow on staggered grids

The computation of flows with large density contrasts is notoriously difficult. To alleviate the difficulty we consider a consistent mass and momentum-conserving discretization of the Navier-Stokes equation. Incompressible flow with capillary forces is modelled and the discretization is performed on a staggered grid of Marker and Cell type. The Volume-of-Fluid method is used to track the interface and a Height-Function method is used to compute surface tension. The advection of the volume fraction is performed using either the Lagrangian-Explicit / CIAM (Calcul d'Interface Affine par Morceaux) method or the Weymouth and Yue (WY) Eulerian-Implicit method. The WY method conserves fluid mass to machine accuracy provided incompressiblity is satisfied which leads to a method that is both momentum and mass-conserving. To improve the stability of these methods momentum fluxes are advected in a manner "consistent" with the volume-fraction fluxes, that is a discontinuity of the momentum is advected at the same speed as a discontinuity of the density. To find the density on the staggered cells on which the velocity is centered, an auxiliary reconstruction of the density is performed. The method is tested for a droplet without surface tension in uniform flow, for a droplet suddenly accelerated in a carrying gas at rest at very large density ratio without viscosity or surface tension, for the Kelvin-Helmholtz instability, for a falling raindrop and for an atomizing flow in air-water conditions

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

A New Class of Nonsingular Exact Solutions for Laplacian Pattern Formation

We present a new class of exact solutions for the so-called {\it Laplacian Growth Equation} describing the zero-surface-tension limit of a variety of 2D pattern formation problems. Contrary to common belief, we prove that these solutions are free of finite-time singularities (cusps) for quite general initial conditions and may well describe real fingering instabilities. At long times the interface consists of N separated moving Saffman-Taylor fingers, with ``stagnation points'' in between, in agreement with numerous observations. This evolution resembles the N-soliton solution of classical integrable PDE's.Comment: LaTeX, uuencoded postscript file

Direct numerical simulations of fluid flow, heat transfer and phase changes

Direct numerical simulations of fluid flow, heat transfer, and phase changes are presented. The simulations are made possible by a recently developed finite difference/front tracking method based on the one-field formulation of the governing equations where a single set of conservation equations is written for all the phases involved. The conservation equations are solved on a fixed rectangular grid, but the phase boundaries are kept sharp by tracking them explicitly by a moving grid of lower dimension. The method is discussed and applications to boiling heat transfer and the solidification of drops colliding with a wall are shown

A phase-field model of Hele-Shaw flows in the high viscosity contrast regime

A one-sided phase-field model is proposed to study the dynamics of unstable interfaces of Hele-Shaw flows in the high viscosity contrast regime. The corresponding macroscopic equations are obtained by means of an asymptotic expansion from the phase-field model. Numerical integrations of the phase-field model in a rectangular Hele-Shaw cell reproduce finger competition with the final evolution to a steady state finger the width of which goes to one half of the channel width as the velocity increases