Viscous irrotational analysis of the deformation and break-up time of a bubble or drop in uniaxial straining flow

2011 ◽  
Vol 688 ◽  
pp. 390-421
Author(s):  
J. C. Padrino ◽  
D. D. Joseph

AbstractThe nonlinear deformation and break-up of a bubble or drop immersed in a uniaxial extensional flow of an incompressible viscous fluid is analysed by means of viscous potential flow. In this approximation, the flow field is irrotational and viscosity enters through the balance of normal stresses at the interface. The governing equations are solved numerically to track the motion of the interface by coupling a boundary-element method with a time-integration routine. When break-up occurs, the break-up time computed here is compared with results obtained elsewhere from numerical simulations of the Navier–Stokes equations (Revuelta, Rodríguez-Rodríguez & Martínez-Bazán J. Fluid Mech., vol. 551, 2006, p. 175), which thus keeps vorticity in the analysis, for several combinations of the relevant dimensionless parameters of the problem. For the bubble, for Weber numbers $3\leqslant \mathit{We}\leqslant 6$, predictions from viscous potential flow shows good agreement with the results from the Navier–Stokes equations for the bubble break-up time, whereas for larger $\mathit{We}$, the former underpredicts the results given by the latter. When viscosity is included, larger break-up times are predicted with respect to the inviscid case for the same $\mathit{We}$. For the drop, and considering moderate Reynolds numbers, $\mathit{Re}$, increasing the viscous effects of the irrotational motion produces large, elongated drops that take longer to break up in comparison with results for inviscid fluids. For larger $\mathit{Re}$, it comes as a surprise that break-up times smaller than the inviscid limit are obtained. Unfortunately, results from numerical analyses of the incompressible, unsteady Navier–Stokes equations for the case of a drop have not been presented in the literature, to the best of the authors’ knowledge; hence, comparison with the viscous irrotational analysis is not possible.

Author(s):  
Mohammad Taeibi-Rahni ◽  
Shervin Sharafatmand

The consistent behavior of non-dimensional parameters on the formation and break up of large cylindrical droplets has been studied by direct numerical simulations (DNS). A one-fluid model with a finite difference method and an advanced front tracking scheme was employed to solve unsteady, incompressible, viscous, immiscible, multi-fluid, two-dimensional Navier-Stokes equations. This time dependent study allows investigation of evolution of the droplets in different cases. For moderate values of Atwood number (AT), increasing Eotvos number (Eo) explicitly increases the deformation rate in both phenomena. Otherwise, raising the Ohnesorge number (Oh) basically amplifies the viscous effects.


2008 ◽  
Vol 605 ◽  
pp. 1-18 ◽  
Author(s):  
SADEGH DABIRI ◽  
WILLIAM A. SIRIGNANO ◽  
DANIEL D. JOSEPH

The flow in a plane liquid jet from an aperture is obtained by direct simulation of the Navier–Stokes equations. The gas–liquid interface is tracked using the level set method. Flows are calculated for different Reynolds and Weber numbers. When We = ∞, the maximum value of the discharge coefficient appears around Re = O(100). The regions that are vulnerable to cavitation owing to the total stress are identified from calculations based on Navier–Stokes equations and viscous potential flow; the two calculations yield similar results for high Reynolds numbers. We prove that the classical potential flow solution does not give rise to a normal component of the rate of strain at the free streamline. Therefore, the normal component of the irrotational viscous stresses also vanishes and cannot change the shape of the free surface. The results of calculations of flows governed by the Navier–Stokes equations are close to those for viscous potential flow outside the vorticity layers at solid boundaries. The Navier–Stokes solutions for the axisymmetric aperture are also given for two values of Reynolds numbers. The results for axisymmetric and planar apertures are qualitatively similar, but the axisymmetric apertures have a lower discharge coefficient and less contraction.


2009 ◽  
Vol 642 ◽  
pp. 421-443 ◽  
Author(s):  
R. M. J. KRAMER ◽  
D. I. PULLIN ◽  
D. I. MEIRON ◽  
C. PANTANO

The single-mode Richtmyer–Meshkov instability is investigated using a first-order perturbation of the two-dimensional Navier–Stokes equations about a one-dimensional unsteady shock-resolved base flow. A feature-tracking local refinement scheme is used to fully resolve the viscous internal structure of the shock. This method captures perturbations on the shocks and their influence on the interface growth throughout the simulation, to accurately examine the start-up and early linear growth phases of the instability. Results are compared to analytic models of the instability, showing some agreement with predicted asymptotic growth rates towards the inviscid limit, but significant discrepancies are noted in the transient growth phase. Viscous effects are found to be inadequately predicted by existing models.


Author(s):  
Pierre Ferrant ◽  
Lionel Gentaz ◽  
Bertrand Alessandrini ◽  
Romain Luquet ◽  
Charles Monroy ◽  
...  

This paper documents recent advances of the SWENSE (Spectral Wave Explicit Navier-Stokes Equations) approach, a method for simulating fully nonlinear wave-body interactions including viscous effects. The methods efficiently combines a fully nonlinear potential flow description of undisturbed wave systems with a modified set of RANS with free surface equations accounting for the interaction with a ship or marine structure. Arbitrary incident wave systems may be described, including regular, irregular waves, multidirectional waves, focused wave events, etc. The model may be fixed or moving with arbitrary speed and 6 degrees of freedom motion. The extension of the SWENSE method to 6 DOF simulations in irregular waves as well as to manoeuvring simulations in waves are discussed in this paper. Different illlustative simulations are presented and discussed. Results of the present approach compare favorably with available reference results.


2014 ◽  
Author(s):  
P. Bigay ◽  
A. Bardin ◽  
G. Oger ◽  
D. Le Touzé

In order to efficiently address complex problems in hydrodynamics, the advances in the development of a new method are presented here. This method aims at finding a good compromise between computational efficiency, accuracy, and easy handling of complex geometries. The chosen method is an Explicit Cartesian Finite Volume method for Hydrodynamics (ECFVH) based on a compressible (hyperbolic) solver, with a ghost-cell method for geometry handling and a Level-set method for the treatment of biphase-flows. The explicit nature of the solver is obtained through a weakly-compressible approach chosen to simulate nearly-incompressible flows. The explicit cell-centered resolution allows for an efficient solving of very large simulations together with a straightforward handling of multi-physics. A characteristic flux method for solving the hyperbolic part of the Navier-Stokes equations is used. The treatment of arbitrary geometries is addressed in the hyperbolic and viscous framework. Viscous effects are computed via a finite difference computation of viscous fluxes and turbulent effects are addressed via a Large-Eddy Simulation method (LES). The Level-Set solver used to handle biphase flows is also presented. The solver is validated on 2-D test cases (flow past a cylinder, 2-D dam break) and future improvements are discussed.


2019 ◽  
Vol 84 (4) ◽  
pp. 696-711 ◽  
Author(s):  
Qianxi Wang ◽  
WenKe Liu ◽  
David M Leppinen ◽  
A D Walmsley

Abstract This paper is concerned with microbubble dynamics in a viscous compressible liquid near a rigid boundary. The compressible effects are modelled using the weakly compressible theory of Wang & Blake (2010, Non-spherical bubble dynamics in a compressible liquid. Part 1. Travelling acoustic wave. J. Fluid Mech., 730, 245–272), since the Mach number associated is small. The viscous effects are approximated using the viscous potential flow theory of Joseph & Wang (2004, The dissipation approximation and viscous potential flow. J. Fluid Mech., 505, 365–377), because the flow field is characterized as being an irrotational flow in the bulk volume but with a thin viscous boundary layer at the bubble surface. Consequently, the phenomenon is modelled using the boundary integral method, in which the compressible and viscous effects are incorporated into the model through including corresponding additional terms in the far field condition and the dynamic boundary condition at the bubble surface, respectively. The numerical results are shown in good agreement with the Keller–Miksis equation, experiments and computations based on the Navier–Stokes equations. The bubble oscillation, topological transform, jet development and penetration through the bubble and the energy of the bubble system are simulated and analysed in terms of the compressible and viscous effects.


Sign in / Sign up

Export Citation Format

Share Document