Dynamics of simultaneously impinging drops on a dry surface: Role of impact velocity and air inertia.

Three dimensional simulations are performed to investigate the interaction dynamics between two drops impinging simultaneously on a dry surface. Of particular interest in this study is to understand the effects of impact velocity and surrounding gas density on droplet interactions. To simulate the droplet dynamics and morphologies, a computational framework based on the phase-field lattice Boltzmann formulation is employed for the two-phase flow computations involving high density ratio. Two different coalescence modes are identified when the impinging droplets have different impact speeds. When one of the droplet has a tangential impact velocity component, asymmetric ridge formation is observed. Influence of droplet impact angle on the interaction dynamics of the central ridge is further investigated. Traces of different fluid particles are seeded to analyse internal flow dynamics in oblique impact scenarios. Greater overlapping between the fluid particles is observed with increase in the impact angle. Finally, the present simulations indicate that the ambient gas density has a significant influence to determine the final outcome of the droplet interactions.

Rebound suppression of a droplet impacting on an oscillating horizontal surface.

The behavior of a droplet impinging onto a solid substrate can be influenced significantly by the horizontal motion of the substrate. The coupled interactions between the moving wall and the impacting droplet may result in various outcomes, which may be different from the usual normal droplet impact on a stationary wall. In this paper, we present a method to suppress drop rebound on hydrophobic surfaces via transverse wall oscillations, normal to the impact direction. The numerical investigation shows that the suppression of droplet rebound has a direct relationship with the oscillation phase, amplitude, and frequency. For a particular range of oscillation frequencies and amplitudes, a lateral shifting of the droplet position is observed along the oscillating direction. While large oscillation amplitude favors the process of droplet deposition, a high frequency promotes droplet rebound from the oscillating wall. A linear trend in the transition region between deposition and rebound is observed from a scaled phase diagram of the oscillation amplitude versus frequency. We provide a systematic investigation of drop deposition and elucidate the mechanism of rebound suppression through the temporal evolution of the nonaxial kinetic energy and the velocity flow field.

A review of numerical methods for red blood cell flow simulation.

In this review, we provide an overview of the simulation techniques employed for modelling the flow of red blood cells (RBCs) in blood plasma. The scope of this review omits the fluid modelling aspect while focusing on other key components in the RBC-plasma model such as (1) describing the RBC deformation with shell-based and spring-based RBC models, (2) constitutive models for RBC aggregation based on bridging theory and depletion theory and (3) additional strategies required for completing the RBC-plasma flow model. These include topics such as modelling fluid-structure interaction with the immersed boundary method and boundary integral method, and updating the variations in multiphase fluid property through the employment of index field methods. Lastly, we summarily discuss the current state and aims of RBC modelling and suggest some research directions for the further development of this field of modelling.

Effect of deformability difference between two erythrocytes on their aggregation.

In this study, we investigated the rheology of a doublet that is an aggregate of two red blood cells (RBCs). According to previous studies, most aggregates in blood flow consist of RBC doublet-pairs and thus the understanding of doublet dynamics has scientific importance in describing its hemodynamics. The RBC aggregation tendency can be significantly affected by the cell's deformability which can vary under both physiological and pathological conditions. Hence, we conducted a two-dimensional simulation of doublet dynamics under a simple shear flow condition with different deformability between RBCs. To study the dissociation process of the doublet, we employed the aggregation model described by the Morse-type potential function, which is based on the depletion theory. In addition, we developed a new method of updating fluid property to consider viscosity difference between RBC cytoplasm and plasma. Our results showed that deformability difference between the two RBCs could significantly reduce their aggregating tendency under a shear condition of 50 s(-1), resulting in disaggregation. Since even under physiological conditions, the cell deformability may be significantly different, consideration of the difference in deformability amongst RBCs in blood flow would be needed for the hemodynamic studies based on a numerical approach.

Numerical simulations of deformation and aggregation of red blood cells in shear flow.

This article reviews numerical simulations of red blood cells (RBCs) mainly using the lattice Boltzmann method (LBM), focusing on the 2-dimensional deformation and aggregation of the cells in simple shear flow. We outline the incorporation of the immersed boundary method into the LBM, in which the membrane forces are obtained from the membrane model. The RBCs are simulated as a single biconcave capsule and as a doublet of biconcave capsules. The transition from swinging to tumbling motions of the RBCs, as induced by reducing the shear rate or increasing the membrane bending stiffness, is discussed. Also discussed is the aggregation tendency of the doublet of RBCs, for which homogenous deformability maintained RBC aggregation, whereas an increased deformability difference resulted in RBC dissociation.

Mass transport in a microchannel bioreactor with a porous wall.

A two-dimensional flow model has been developed to simulate mass transport in a microchannel bioreactor with a porous wall. A two-domain approach, based on the finite volume method, was implemented. For the fluid part, the governing equation used was the Navier-Stokes equation; for the porous medium region, the generalized Darcy-Brinkman-Forchheimer extended model was used. For the porous-fluid interface, a stress jump condition was enforced with a continuity of normal stress, and the mass interfacial conditions were continuities of mass and mass flux. Two parameters were defined to characterize the mass transports in the fluid and porous regions. The porous Damkohler number is the ratio of consumption to diffusion of the substrates in the porous medium. The fluid Damkohler number is the ratio of the substrate consumption in the porous medium to the substrate convection in the fluid region. The concentration results were found to be well correlated by the use of a reaction-convection distance parameter, which incorporated the effects of axial distance, substrate consumption, and convection. The reactor efficiency reduced with reaction-convection distance parameter because of reduced reaction (or flux), and smaller local effectiveness factor due to the lower concentration in Michaelis-Menten type reactions. The reactor was more effective, and hence, more efficient with the smaller porous Damkohler number. The generalized results could find applications for the design of bioreactors with a porous wall.

Numerical simulation of mass transport in a microchannel bioreactor with cell micropatterning.

Micropatterning of two different cell types based on surface modification allows spatial control over two distinct cell subpopulations. This study considers a micropatterned coculture system, which has release and absorption parts alternately arranged at the base, and each part has a single cell type. A micropattern unit was defined and within each unit, there are one release part and one absorption part. The cells in the absorption parts consume species, which are secreted by the cells in the release parts. The species concentrations at the micropatterned cell base were computed from a three-dimensional numerical flow model incorporating mass transport. Different combined parameters were developed for the release and absorption parts to make the data collapse in each part. Combination of the collapse data in the release and absorption parts can be used to predict the concentration distribution through the whole channel. The correlated results were applied to predict the critical length ratio of the release and absorption parts for an actual micropatterned system (Bhatia et al., 1999, "Effect of Cell-Cell Interactions in Preservation of Cellular Phenotype: Co-Cultivation of Hepatocytes and Nonparenchymal Cell," FASEB J. 13, pp. 1883-1900) to avoid species insufficiency based on basic fibroblast growth factor (bFGF). The mass transfer effectiveness was found to be higher with more numbers of micropattern units. The optimal condition for micropatterned coculture bioreactors is achieved by having the product of the length ratio and the reaction ratio equal to 1. This condition was used to optimize the mass transfer in the micropatterned system (Bhatia et al., 1999, "Effect of Cell-Cell Interactions in Preservation of Cellular henotype: Co-Cultivation of Hepatocytes and Nonparenchymal Cell," FASEB J. 13, pp. 1883-1900) based on bFGF.

Numerical simulation on mass transport in a microchannel bioreactor for co-culture applications.

Microchannel bioreactors have applications for manipulating and investigating the fluid microenvironment on cell growth and functions in either single culture or co-culture. This study considers two different types of cells distributed randomly as a co-culture at the base of a microchannel bioreactor: absorption cells, which only consume species based on the Michaelis-Menten process, and release cells, which secrete species, assuming zeroth order reaction, to support the absorption cells. The species concentrations at the co-culture cell base are computed from a three-dimensional numerical flow-model incorporating mass transport. Combined dimensionless parameters are proposed for the co-culture system, developed from a simplified analysis under the condition of decreasing axial-concentration. The numerical results of species concentration at the co-culture cell-base are approximately correlated by the combined parameters under the condition of positive flux-parameter. Based on the correlated results, the critical value of the inlet concentration is determined, which depends on the effective microchannel length. For the flow to develop to the critical inlet concentration, an upstream length consisting only of release cells is needed; this upstream length is determined from an analytical solution. The generalized results may find applications in analyzing the mass transport requirements in a co-culture microchannel bioreactor.

Mass transport and shear stress in a microchannel bioreactor: numerical simulation and dynamic similarity.

Microchannel bioreactors have been used in many studies to manipulate and investigate the fluid microenvironment around cells. In this study, substrate concentrations and shear stresses at the base were computed from a three-dimensional numerical flow-model incorporating mass transport. Combined dimensionless parameters were developed from a simplified analysis. The numerical results of substrate concentration were well correlated by the combined parameters. The generalized results may find applications in design analysis of microchannel bioreactors. The mass transport and shear stress were related in a generalized result. Based on the generalized results and the condition of dynamic similarity, various means to isolate their respective effects on cells were considered.