In this paper, the transient settling dynamics of a spherical particle sedimenting in a linearly stratified fluid is investigated by performing fully resolved direct numerical simulations. The settling behaviour is quantified for different values of Reynolds, Froude and Prandtl numbers. It is demonstrated that the transient settling dynamics is correlated to the induced Lagrangian drift of flow around the settling particle. A simplified model is provided to predict the maximum velocity of the settling particle in linearly stratified fluids. The peak velocity can be followed by the oscillation of the settling velocity and the particle can even reverse its direction of motion before reaching to its neutrally buoyant level. The frequency of oscillation of settling velocity scales with the Brunt–Väisälä frequency and the motion of the particle can lead to the formation of secondary and tertiary vortices following the primary vortex.

Purpose -This paper presents a first step toward developing a comprehensive methodology for fully resolved numerical simulations of fusion deposition modeling. Design/methodology/approach -A front-tracking/finite volume method previously developed for simulations of multiphase flows is extended to model the injection of hot polymer and its cooling down. Findings -The accuracy and convergence properties of the new method are tested by grid refinement and the method is shown to produce convergent solutions for the shape of the filament, the temperature distribution, contact area and reheat region when new filaments are deposited on top of previously laid down filaments. Research limitations/implications -The present paper focuses on modeling the fluid flow and the cooling. The modeling of solidification, volume changes and residual stresses will be described in Part II. Practical implications -The ability to carry our fully resolved numerical simulations of the fusion deposition process is expected to help explore new deposition strategies and to provide the "ground truth" for the development of reduced order models. Originality/value -The present paper is the first fully resolved simulation of the deposition in fusion filament modeling.

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, focusing on the mass transfer from bubbles in liquids with low mass diffusivity, suggests one approach. We conclude by discussing immediate needs for progress on the theoretical framework for describing the large-scale motion of multiphase flows and the need for multiscale methods to capture physical processes taking place at diverse length and time scales.

No abstract

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 Eötvös 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.

No abstract

In this paper, we study the dynamics of deformable cells in a channel flow of Newtonian and polymeric fluids and unravel the effects of deformability, elasticity, inertia, and size on the cell motion. We investigate the role of polymeric fluids on the cell migration behavior and the performance of inertial microfluidic devices. Our results show that the equilibrium position of the cell is on the channel diagonal, in contrast to that of rigid particles, which is on the center of the channel faces for the same range of Reynolds number. A constant-viscosity polymeric fluid, modeled using an Oldroyd-B constitutive equation, drives the cells toward the channel centerline, while a shear-thinning polymeric fluid, modeled using a Giesekus constitutive equation, pushes the cells toward the channel wall. The findings of this paper suggest that the addition of polymers in microfluidic devices can be used to enhance the throughput of cell focusing and separation devices at a low cost. This study provides an insight on the role of rheological properties of the fluid and the ways that they can be tuned to control the focal position of the cells.

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 number between 100 and 2000. The orifice has a rounded upstream corner and a sharp downstream corner with length-to-diameter ratio between 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 is studied. The Navier-Stokes equations are solved numerically using a finitevolume method and a boundary-fitted orthogonal grid that comes from the streamlines and potential lines of an axisymmetric equipotential flow in the same geometry. A level-set formulation is used to track the interface and model the surface tension.