Sadegh Dabiri

Cavitation in an orifice flow

The purpose of this study is to identify the potential locations for cavitation induced by total stress on the flow of a liquid through an orifice of an atomizer. A numerical simulation of two-phase incompressible flow is conducted in an axisymmetric geometry of the orifice for Reynolds numbers between 100 and 2000. The orifice has a rounded upstream corner and a sharp downstream corner with length-to-diameter ratio between 0.1 and 5. The total stress including viscous stress and pressure has been calculated in the flow field and, from there, the maximum principal stress is found. The total-stress criterion for cavitation is applied to find the regions where cavitation is likely to occur and compared with those of the traditional pressure criterion. Results show that the viscous stress has significant effects on cavitation. The effect of geometry and occurrence of hydraulic flip in the orifice on the total stress are studied. The Navier-Stokes equations are solved numerically using a finite-volume method ...

Transition between regimes of a vertical channel bubbly upflow due to bubble deformability

The effect of the deformability of viscous bubbles on the flow rate of bubbly upflow in a vertical channel is examined using direct numerical simulations. A sharp transition between two different flow regimes has been observed. At large bubble deformability, characterized by large Eotvos number (Eo), the flow rate is close to the single phase flow rate, with adjusted pressure gradient, and the bubbles are almost uniformly distributed in the middle of the channel. On the other hand, at low Eo the bubbles are concentrated near channel walls and flow rates are much smaller than the single phase flow. The transition from high flow rate to low flow rate occurs rather abruptly. It is found that the transition occurs when the less deformable bubbles enter the viscous sublayer due to the lateral lift force on the bubbles. This leads to an increase in the viscous dissipation near the wall which leads to a decrease in the flow rate.

Multiscale considerations in direct numerical simulations of multiphase flowsa)

Direct Numerical Simulations of multiphase flows have progressed rapidly over the last decade and it is now possible to simulate, for example, the motion of hundreds of deformable bubbles in turbulent flows. The availability of results from such simulations should help advance the development of new and improved closure relations and models of the average or large-scale flows. We review recent results for bubbly flow in vertical channels, discuss the difference between upflow and downflow and the effect of the bubble deformability and how the resulting insight allowed us to produce a simple description of the large scale flow, for certain flow conditions. We then discuss the need for the development of numerical methods for more complex situations, such as where the flow creates spontaneous thin films and threads, or where additional physical processes take place at a rate that is very different from the fluid flow. Recent work on capturing localized small-scale processes using embedded analytical models,...

Interaction between a cavitation bubble and shear flow

The deformation of a cavitation bubble in shear and extensional flows is studied numerically. The Navier–Stokes equations are solved to observe the three-dimensional behaviour of the bubble as it grows and collapses. During the collapse phase of the bubble, two re-entrant jets are observed on two sides of the bubble. The re-entrant jets are not the result of interaction with a solid wall or free surface; rather, they are formed due to interaction of the bubble with the background flow. Effects of the viscosity, surface tension and shear rate on the formation and strength of re-entrant jets are investigated. Re-entrant jets with enough strength break up the bubble into smaller bubbles. Post-processing and analysis of the results are done to cast the disturbance by the bubble on the liquid velocity field in terms of spherical harmonics. It is found that quadrupole moments are created in addition to the monopole source.

Two-dimensional and axisymmetric viscous flow in apertures

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.

On the rising motion of a drop in stratified fluids

The rising dynamics of a deformable drop in a linearly stratified fluid is numerically obtained using a finite-volume/front-tracking method. Our results show that the drag coefficient of a spherical drop in a stratified fluid enhances as Cd,s/Cd,h−1∼Frd−2.86 for drop Froude numbers in the range of 4 < Frd < 16. The role of the deformability of the drop on the temporal evolution of the motion is investigated along with stratification and inertial effects. We also present the important role of stratification on the transient rising motion of the drop. It is shown that a drop can levitate in the presence of a vertical density gradient. The drop undergoes a fading oscillatory motion around its neutrally buoyant position except for high viscosity ratio drops where the oscillation occurs around a density level lighter than the neutral buoyancy level. In addition, a detailed characterization of the flow signature of a rising drop in a linearly stratified fluid including the buoyancy induced vortices and the resu...

Deformation of a droplet in a particulate shear flow

A numerical simulation using a distributed-Lagrange-multiplier-based computational method is conducted in order to investigate the deformation and breakup of a droplet in a particulate shear flow. A level-set method is used to track the interface and model the surface tension. The results show that the presence of particles leads to larger droplet deformation and a perforation in the center of the droplet. It is found that the critical Stokes number above which a perforation occurs increases linearly with the inverse of the capillary number and viscosity ratio.

Scaling law for bubbles rising near vertical walls

This paper examines the rising motion of a layer of gas bubbles next to a vertical wall in a liquid in the presence of an upward flow parallel to the wall to help with the understanding of the fluid dynamics in a bubbly upflow in vertical channels. Only the region near the wall is simulated with an average pressure gradient applied to the domain that balances the weight of the liquid phase. The upward flow is created by the rising motion of the bubbles. The bubbles are kept near the wall by the lateral lift force acting on them as a result of rising in the shear layer near the wall. The rise velocity of the bubbles sliding on the wall and the average rise velocity of the liquid depend on three dimensionless parameters, Archimedes number, Ar, Eotvos number, Eo, and the average volume fraction of bubbles on the wall. In the limit of small Eo, bubbles are nearly spherical and the dependency on Eo becomes negligible. In this limit, the scaling of the liquid Reynolds number with Archimedes number and the void fraction is presented. A scaling argument is presented based on viscous dissipation analysis that matches the numerical findings. Viscous dissipation rates are found to be high in a thin film region between the bubble and the wall. A scaling of the viscous dissipation and steady state film thickness between the bubble and the wall with Archimedes number is presented.

Fully Resolved Numerical Simulations of Fused Deposition Modeling

Fully resolved numerical simulations of the deposition of a filament of hot, viscous liquid, and its solidification, are presented. The computations are done using a front tracking/finite volume code designed for direct numerical simulations of multiphase flows, but here we add a volume source that follows a prescribed trajectory to inject the melt. The cooling of the melt is simulated and once it is cool enough the melt is solidified. The result suggest that the deposition is fully captured and that remelting of already solidified material results in a compaction of the deposited filaments.Copyright

A numerical study on the effects of cavitation on orifice flow

Previous experimental studies have shown better atomization of sprays generated by high-pressure liquid injectors when cavitation occurs inside the nozzle. It has been proposed that the collapse of traveling cavitation bubbles increases the disturbances inside the liquid flow. These disturbances will later trigger the instabilities in the emerged jet and cause a shorter breakup distance. In this paper, effects of cavitation on increasing the disturbances in the flow through the orifice of an atomizer are studied. In previous cavitation models, spherical cavitation bubbles are considered. Here, the cavitation bubbles are allowed to deform as they travel through the orifice. Dynamics of the cavitation bubble, traveling in the separated shear layer in the orifice, is analyzed through a one-way coupling between the orifice flow and bubble dynamics. Effects of shear strain, normal strain, and pressure variation are examined. Three mechanisms are suggested that could be responsible for the increase in disturban...