Particle manipulation via integration of electroosmotic flow of power-law fluids with standing surface acoustic waves (SSAW)

Standing surface acoustic wave (SSAW) based microfluidic devices have shown great promise toward fluid and particle manipulation applications in medicine, chemistry, and biotechnology. In this article, we present an analytical model for investigating continuous manipulation of particles (both synthetic and biological) within electroosmotic flow of non-Newtonian bio-fluids in a microfluidic channel under the influence of standing surface acoustic waves (SSAW). The particles are injected along the center of channel into the electroosmotically driven flow of power-law fluids, wherein their transport through the SSAW region is dictated by the hydrodynamic, electrophoretic, and acoustic forces. We first present a mathematical model to analyze the characteristics of electroosmotic flow of non-Newtonian power-law fluids in a hydrophobic slit microchannel. Next, we investigate the trajectories of particles in the flow field due to the combined effect of electroosmotic, electrophoretic, and acoustophoretic forcing mechanisms. The effect of key parameters such as particle size, their physical properties, input power, flow rate, and flow behavior index on the particle trajectories is examined while including the effect of the channel walls. The presented model delineates the methodologies of improving SSAW-based particle separation technology by considering the fluid rheology as well as the surface properties of the channel walls. Therefore, we believe that this model can serve as an efficient tool for device design and quick optimizations to explore novel applications concerning the integration of electroosmotic flows with acoustofluidic technologies.     

Effects of particle size and concentration on bubble coalescence and froth formation in a slurry bubble column

A new approach for simulating the formation of a froth layer in a slurry bubble column is proposed. Froth is considered a separate phase, comprised of a mixture of gas, liquid, and solid. The simulation was carried out using commercial flow simulation software (FIRE v2014) for particle sizes of 60–150 μm at solid concentrations of 0–40 vol%, and superficial gas velocities of 0.02–0.034 m/s in a slurry bubble column with a hydraulic diameter of 0.2 m and height of 1.2 m. Modelling calculations were conducted using a Eulerian–Eulerian multiphase approach with k–ε turbulence. The population balance equations for bubble breakup, bubble coalescence rate, and the interfacial exchange of mass and momentum were included in the computational fluid dynamics code by writing subroutines in Fortran to track the number density of different bubble sizes. Flow structure, radial gas holdup, and Sauter mean bubble diameter distributions at different column heights were predicted in the pulp zone, while froth volume fraction and density were predicted in the froth zone. The model was validated using available experimental data, and the predicted and experimental results showed reasonable agreement. To demonstrate the effect of increasing solid concentration on the coalescence rate, a solid-effect multiplier in the coalescence efficiency equation was used. The solid-effect multiplier decreased with increasing slurry concentration, causing an increase in bubble coalescence efficiency. A slight decrease in the coalescence efficiency was also observed owing to increasing particle size, which led to a decrease in Sauter mean bubble diameter. The froth volume fraction increased with solid concentration. These results provide an improved understanding of the dynamics of slurry bubble reactors in the presence of hydrophilic particles.     

A cost‐effective curvature calculation approach for interfacial flows on unstructured meshes

We present a simple and cost‐effective curvature calculation approach for simulations of interfacial flows on structured and unstructured grids. The interface is defined using volume fractions, and the interface curvature is obtained as a function of the gradients of volume fractions. The gradient computation is based on a recently proposed gradient recovery method that mimicks the least squares approach without the need to solve a system of equations and is quite easy to implement on arbitrary polygonal meshes. The resulting interface curvature is used in a continuum surface force formulation within the framework of a well‐balanced finite‐volume algorithm to simulate multiphase flows dominated by surface tension. We show that the proposed curvature calculation is at least as accurate as some of the existing approaches on unstructured meshes while being straightforward to implement on any mesh topology. Numerical investigations also show that spurious currents in stationary problems that are dependent on the curvature calculation methodology are also acceptably low using the proposed approach. Studies on capillary waves and rising bubbles in viscous flows lend credence to the ability of the proposed method as an inexpensive, robust, and reasonably accurate approach for curvature calculation and numerical simulation of multiphase flows.     

Hybrid Isothermal Model for the Ferrohydrodynamic Chemically Reactive Species

A hybrid isothermal model for the homogeneous-heterogeneous reactions in ferrohydrodynamic boundary layer flow is established. The characteristics of Newtonian heating and magnetic dipole in a ferrofluid due to a stretchable surface is analyzed for three chemical species. It is presumed that the isothermal cubic autocatalator kinetic gives the homogeneous reaction and the first order kinetics gives the heterogeneous (surface) reaction. The analysis is carried out for equal diffusion coefficients of all autocatalyst and reactions. Heat flux is examined by incorporating Fourier's law of heat conduction. Characteristics of materialized parameters on the magneto-thermomechanical coupling in the flow of a chemically reactive species are investigated. Further, the heat transfer rate and friction drag are depicted for the ferrohydrodynamic chemically reactive species. It is evident that the Schmidt number has increasing behavior on the rate of heat transfer in the boundary layer. Comparison with available results for specific cases is found an excellent agreement.     

Gibbsian dynamics and ergodicity of magnetohydrodynamics and related systems forced by random noise

The magnetohydrodynamics system consists of the Navier-Stokes equations from fluid mechanics, coupled with the Maxwell’s equations from electromagnetism through multiples of non-linear terms involving derivatives. Following the approach of [1 Weinan, E., Mattingly, J. C., Sinai, Y. (2001). Gibbsian dynamics and ergodicity for the stochastically forced Navier-Stokes equation. Commun. Math. Phys. 224:83–106. DOI:10.1007/s002201224083.[Crossref], [Web of Science ®] , [Google Scholar]], we prove the existence of a unique invariant measure in case the forcing terms consist of the cylindrical Wiener processes with only low modes. Its proof requires taking advantage of the structure of the non-linear terms carefully and is extended to various other related models such as the magnetohydrodynamics-Boussinesq system from fluid mechanics in atmosphere and oceans, as well as the magneto-micropolar fluid system from the theory of microfluids.     

The Effects of Wettability on Primary Vortex and Secondary Flow in Three-Dimensional Rotating Fluid

The secondary flow driven by the primary vortex in a cylinder,generating the so called "tea leaf paradox",is fundamental for understanding many natural phenomena,industrial applications and scientific researches.In this work,the effect of wettability on the primary vortex and secondary flow is investigated by the three-dimensional multiphase lattice Boltzmann method based on a chemical potential.We find that the surface wettability strongly affects the shape of the primary vortex.With the increase of the contact angle of the cylinder,the sectional plane of the primary vortex gradually changes from a steep valley into a saddle with two raised parts.Because the surface friction is reduced correspondingly,the core of the secondary vortex moves to the centerline of the cylinder and the vortex intensity also increases.The stirring force has stronger effects to enhance the secondary flow and push the vortex up than the surface wettability.Interestingly,a small secondary vortex is discovered near the three-phase contact line when the surface has a moderate wettability,owing to the interaction between the secondary flow and the curved gas/liquid interface.     

基于有限体积虚拟区域方法的两相流动直接模拟

为提高计算效率，提出有限体积法离散下的虚拟区域颗粒两相流动直接模拟方法.在控制方程中加入相应的虚拟区域源项，保证了颗粒内部的刚体运动特性.该源项中含有颗粒信息部分及流体信息部分.在每次迭代后，对源项中的流体信息部分进行更新，从而更好地保证颗粒内速度的刚体分布.计算静止颗粒圆柱绕流及单个颗粒的沉降过程，验证了算法的准确性.     

Helicity and other conservation laws in perfect fluid motion

In this review paper, we discuss helicity from a geometrical point of view and see how it applies to the motion of a perfect fluid. We discuss its relation with the Hamiltonian structure, and then its extension to arbitrary space dimensions. We also comment about the existence of additional conservation laws for the Euler equation, and its unlikely integrability in Liouville's sense.     

The application of the Lattice Boltzmann method to the one-dimensional modeling of pulse waves in elastic vessels

The one-dimensional nonlinear equations for the blood flow motion in distensible vessels are considered using the kinetic approach. It is shown that the Lattice Boltzmann (LB) model for non-ideal gas is asymptotically equivalent to the blood flow equations for compliant vessels at the limit of low Knudsen numbers. The equations of state for non-ideal gas are transformed to the pressure-luminal area response. This property allows to model arbitrary pressure-luminal area relations. Several test problems are considered: the propagation of a sole nonlinear wave in an elastic vessel, the propagation of a pulse wave in a vessel with varying mechanical properties (artery stiffening) and in an artery bifurcation, in the last problem Resistor–Capacitor–Resistor (RCR) boundary conditions are considered. The comparison with the previous results shows a good precision.     

超声搅拌磁流变抛光液的声场仿真分析

为了探究超声搅拌磁流变抛光液的制备及优化工艺,利用多物理场数值计算方法,建立了超声搅拌磁流变抛光液的声场仿真模型,并进行了频域分析。研究了不同液位深度、超声变幅杆探入深度,不同功率下磁流变抛光液的声场分布。通过测量磁流变抛光液的声场强度对声场仿真进行了验证。结果表明：随着距变幅杆距离的增加,声强逐渐减弱,高声强区域主要分布在换能器轴线附近。声强在距变幅杆20mm范围内急剧衰减,变幅杆最佳探入深度为10mm,增大功率有助于空化区域的扩大。声场仿真结果与实验测量结果基本一致,对磁流变抛光液的制备提供了数值计算基础。