Experimental Correlation for Splashing Condition of Droplets on Solid Substrates

Droplet splashing behaviors of water-ethanol binary mixture liquids on roughened solid surfaces were experimentally observed in order to investigate the effects of surface tension, viscosity, and wettability/surface roughness on the splashing occurrence. The range of the droplet volumes was from 1.7 μL to 32.6 μL. The ranges of the surface tension and the viscosity were from 21.1 mN/m to 71.9 mN/m, and from 1 mPas to 2.91 mPas, respectively. The surface roughness range was from 0.03 μm to 1.25 μm for Ra. The present experimental data were evaluated on the basis of the existing models. Resulting from these experiments, a simple model using the Ohnesorge number evaluated by the capillary length was proposed and the accuracy of the predicted critical values such as the critical Weber and Reynolds numbers were discussed. The result indicated that the liquid properties and the quantification of the surface condition such as surface roughness are important factors for the prediction of the splashing behavior.

Sonic Eddy Model of the Turbulent Boundary Layer

The effects of Mach number on the skin friction and velocity fluctuations of the turbulent boundary layer are considered through a sonic eddy model. Originally proposed for free shear flows, the model assumes that the eddies responsible for momentum transfer have a rotation Mach number of unity, with the entrainment rate limited by acoustic signaling. Under this assumption, the model predicts that the skin friction coefficient should go as the inverse Mach number in a regime where the Mach number is larger than unity but smaller than the square root of the Reynolds number. The velocity fluctuations normalized by the friction velocity should be the inverse square root of the Mach number in the same regime. Turbulent transport is controlled by acoustic signaling. The density field adjusts itself such that the Reynolds stresses correspond to the momentum transport. In contrast, the conventional van Driest–Morkovin view is that the Mach number effects are due to density variations directly. A new experiment or simulation is proposed to test this model using different gases in an incompressible boundary layer, following the example of Brown and Roshko in the free shear layer.

Instability and Transition of a Boundary Layer over a Backward-Facing Step

The development of secondary instabilities in a boundary layer over a backward-facing step is investigated numerically. Two step heights are considered, h/δo*=0.5 and 1.0 (where δo* is the displacement thickness at the step location), in addition to a reference flat-plate case. A case with a realistic freestream-velocity distribution is also examined. A controlled K-type transition is initiated using a narrow ribbon upstream of the step, which generates small and monochromatic perturbations by periodic blowing and suction. A well-resolved direct numerical simulation is performed. The step height and the imposed freestream-velocity distribution exert a significant influence on the transition process. The results for the h/δo*=1.0 case exhibit a rapid transition primarily due to the Kelvin–Helmholtz instability downstream of step; non-linear interactions already occur within the recirculation region, and the initial symmetry and periodicity of the flow are lost by the middle stage of transition. In contrast, case h/δo*=0.5 presents a transition road map in which transition occurs far downstream of the step, and the flow remains spatially symmetric and temporally periodic until the late stage of transition. A realistic freestream-velocity distribution (which induces an adverse pressure gradient) advances the onset of transition to turbulence, but does not fundamentally modify the flow features observed in the zero-pressure gradient case. Considering the budgets of the perturbation kinetic energy, both the step and the induced pressure-gradient increase, rather than modify, the energy transfer.

A Compressible Turbulence Model for Pressure—Strain

This work focuses on the performance and validation of compressible turbulence models for the pressure-strain correlation. Considering the Launder Reece and Rodi (LRR) incompressible model for the pressure-strain correlation, Adumitroaie et al., Huang et al., and Marzougui et al., used different modeling approaches to develop turbulence models, taking into account compressibility effects for this term. Two numerical coefficients are dependent on the turbulent Mach number, and all of the remaining coefficients conserve the same values as in the original LRR model. The models do not correctly predict the compressible turbulence at a high-speed shear flow. So, the revision of these models is the major aim of this study. In the present work, the compressible model for the pressure-strain correlation developed by Khlifi−Lili, involving the turbulent Mach number, the gradient, and the convective Mach numbers, is used to modify the linear mean shear strain and the slow terms of the previous models. The models are tested in two compressible turbulent flows: homogeneous shear flow and the newly developed plane mixing layers. The predicted results of the proposed modifications of the Adumitroaie et al., Huang et al., and Marzougui et al., models and of its universal versions are compared with direct numerical simulation (DNS) and experiment data. The results show that the important parameters of compressibility in homogeneous shear flow and in the mixing layers are well predicted by the proposal models.

Analysis of the Shear-Thinning Viscosity Behavior of the Johnson–Segalman Viscoelastic Fluids

This paper presents a numerical comparison of viscoelastic shear-thinning fluid flow using a generalized Oldroyd-B model and Johnson–Segalman model under various settings. Results for the standard shear-thinning generalization of Oldroyd-B model are used as a reference for comparison with those obtained for the same flow cases using Johnson–Segalman model that has specific adjustment of convected derivative to assure shear-thinning behavior. The modeling strategy is first briefly described, pointing out the main differences between the generalized Oldroyd-B model (using the Cross model for shear-thinning viscosity) and the Johnson–Segalman model operating in shear-thinning regime. Then, both models are used for blood flow simulation in an idealized stenosed axisymmetric vessel under different flow rates for various model parameters. The simulations are performed using an in-house numerical code based on finite-volume discretization. The obtained results are mutually compared and discussed in detail, focusing on the qualitative assessment of the most distinct flow field differences. It is shown that despite all models sharing the same asymptotic viscosities, the behavior of the Johnson–Segalman model can be (depending on flow regime) quite different from the predictions of the generalized Oldroyd-B model.

Characterization of Effective Diffusion within Viscoelastic Fluids with Elastic Instabilities

We considered effective diffusion, characterized by magnitudes of effective diffusion coefficients, in order to quantify mass transport due to the onset and development of elastic instabilities. Effective diffusion coefficient magnitudes were determined using different analytic approaches, as they were applied to tracked visualizations of fluorescein dye front variations, as circumferential advection was imposed upon a flow environment produced using a rotating Couette flow arrangement. Effective diffusion coefficient results were provided for a range of flow shear rates, which were produced using different Couette flow rotation speeds and two different flow environment fluid depths. To visualize the flow behavior within the rotating Couette flow environment, minute amounts of fluorescein dye were injected into the center of the flow container using a syringe pump. This dye was then redistributed within the flow by radial diffusion only when no disk rotation was used, and by radial diffusion and by circumferential advection when disk rotation was present. Associated effective diffusion coefficient values, for the latter arrangement, were compared to coefficients values with no disk rotation, which were due to molecular diffusion alone, in order to quantify enhancements due to elastic instabilities. Experiments were conducted using viscoelastic fluids, which were based on a 65% sucrose solution, with different polymer concentrations ranging from 0 ppm to 300 ppm. Associated Reynolds numbers based on the fluid depth and radially averaged maximum flow velocity ranged from 0.00 to 0.5. The resulting effective diffusion coefficient values for different flow shear rates and polymer concentrations quantified the onset of elastic instabilities, as well as significant and dramatic changes to local mass transport magnitudes, which are associated with the further development of elastic instabilities.

Design and Numerical Simulation of Biomimetic Structures to Capture Particles in a Microchannel

The study of separating different sizes of particles through a microchannel has been an interest in recent years and the primary attention of this study is to isolate the particles to the specific outlets. The present work highly focuses on the design and numerical analysis of a microchip and the microparticles capture using special structures like corrugated dragonfly wing structure and cilia walls. The special biomimetic structured corrugated wing is taken from the cross-sectional area of the dragonfly wing and cilia structure is obtained from the epithelium terminal bronchioles to the larynx from the human body. Parametric studies were conducted on different sizes of microchip scaled and tested up in the range between 2–6 mm and the thickness was assigned as 80 µm in both dragonfly wing structure and cilia walls. The microflow channel is a low Reynolds number regime and with the help of the special structures, the flow inside the microchannel is pinched and a sinusoidal waveform pattern is observed. The pinched flow with sinusoidal waveform carries the particles downstream and induces the particles trapped in desired outlets. Fluid particle interaction (FPI) with a time-dependent solver in COMSOL Multiphysics was used to carry out the numerical study. Two particle sizes of 5 µm and 20 µm were applied, the inlet velocity of 0.52 m/s with an inflow angle of 50° was used throughout the study and it suggested that: the microchannel length of 3 mm with corrugated dragonfly wing structure had the maximum particle capture rate of 20 µm at the mainstream outlet. 80% capture rate for the microchannel length of 3 mm with corrugated dragonfly wing structure and 98% capture rate for the microchannel length of 2 mm with cilia wall structure were observed. Numerical simulation results showed that the cilia walled microchip is superior to the corrugated wing structure as the mainstream outlet can conduct most of the 20 µm particles. At the same time, the secondary outlet can laterally capture most of the 5 µm particles. This biomimetic microchip design is expected to be implemented using the PDMS MEMS process in the future.

Characterization of the Ejector Pump Performance for the Assisted Bidirectional Glenn Procedure

This study introduces an algebraic model informed by computational fluid dynamics (CFD) simulations to investigate the performance of the assisted bidirectional Glenn (ABG) operation on a broad range of conditions. The performance of this operation, as measured by the superior vena cava (SVC) pressure, depends on the nozzle area in its ejector pump and the patient’s pulmonary vascular resistance (PVR). Using the developed algebraic model to explore this two-dimensional parameter space shows that the ejector pump can create a pressure difference between the pulmonary artery and the SVC as high as 5 mmHg. The lowest SVC pressure is produced at a nozzle area that decreases linearly with the PVR such that, at PVR =4.2 (Wood units-m2), there is no added benefit in utilizing the ejector pump effect (optimal nozzle area is zero, corresponding to the bidirectional Glenn circulation). At PVR =2 (Wood units-m2), the SVC pressure can be lowered to less than 4 mmHg by using an optimal nozzle area of ≈2.5 mm2. Regardless of the PVR, adding a 2 mm2 nozzle to the baseline bidirectional Glenn boosts the oxygen saturation and delivery by at least 15%. The SVC pressure for that 2 mm2 nozzle remains below 14 mmHg for all PVRs less than 7 Wood units-m2. The mechanical efficiency of the optimal designs consistently remains below 30%, indicating the potential for improvement in the future. A good agreement is observed between the algebraic model and high-fidelity CFD simulations.

Inlet and Outlet Boundary Conditions and Uncertainty Quantification in Volumetric Lattice Boltzmann Method for Image-Based Computational Hemodynamics

Inlet and outlet boundary conditions (BCs) play an important role in newly emerged image-based computational hemodynamics for blood flows in human arteries anatomically extracted from medical images. We developed physiological inlet and outlet BCs based on patients’ medical data and integrated them into the volumetric lattice Boltzmann method. The inlet BC is a pulsatile paraboloidal velocity profile, which fits the real arterial shape, constructed from the Doppler velocity waveform. The BC of each outlet is a pulsatile pressure calculated from the three-element Windkessel model, in which three physiological parameters are tuned by the corresponding Doppler velocity waveform. Both velocity and pressure BCs are introduced into the lattice Boltzmann equations through Guo’s non-equilibrium extrapolation scheme. Meanwhile, we performed uncertainty quantification for the impact of uncertainties on the computation results. An application study was conducted for six human aortorenal arterial systems. The computed pressure waveforms have good agreement with the medical measurement data. A systematic uncertainty quantification analysis demonstrates the reliability of the computed pressure with associated uncertainties in the Windkessel model. With the developed physiological BCs, the image-based computation hemodynamics is expected to provide a computation potential for the noninvasive evaluation of hemodynamic abnormalities in diseased human vessels.

Ultrasonic Atomization: New Spray Characterization Approaches

In particle engineering, spray drying is an essential technique that depends on producing sprays, ideally made of equal-sized droplets. Ultrasonic sprays appear to be the best option to achieve it, and Faraday waves are the background mechanism of ultrasonic atomization. The characterization of sprays in this atomization strategy is commonly related to the relation between characteristic drop sizes and the capillary length produced by the forcing frequency of wavy patterns on thin liquid films. However, although this atomization approach is practical when the intended outcome is to produce sprays with droplets of the same size, drop sizes are diverse in real applications. Therefore, adequate characterization of drop size is paramount to establishing the relations between empirical approaches proposed in the literature and the outcome of ultrasonic atomization in actual operating conditions. In this sense, this work explores new approaches to spray characterization applied to ultrasonic sprays produced with different solvents. The first two introduced are the role of redundancy in drop size measurements to avoid resolution limitation in the measurement technique and compare using regular versus variable bin widths when building the histograms of drop size. Another spray characterization tool is the Drop Size Diversity to understand the limitations of characterizing ultrasonic sprays solely based on representative diameters or moments of drop size distributions. The results of ultrasonic spray characterization obtained emphasize: the lack of universality in the relation between a characteristic diameter and the capillary length associated with Faraday waves; the variability on drop size induced by both liquid properties and flow rate on the atomization outcome, namely, lower capillary lengths produce smaller droplets but less efficiently; the higher sensibility of the polydispersion and heterogeneity degrees in Drop Size Diversity when using variable bin widths to build the histograms of drop size; the higher drop size diversity for lower flow rates expressed by the presence of multiple clusters of droplets with similar characteristics leading to multimodal drop size distributions; and the gamma and log-normal mathematical probability functions are the ones that best describe the organization of drop size data in ultrasonic sprays.