A control volume finite element method for three-dimensional three-phase flows

A novel control volume finite element method with adaptive anisotropic unstructured meshes is presented for three-dimensional three-phase flows with interfacial tension. The numerical framework consists of a mixed control volume and finite element formulation with a new P₁DG-P₂ elements (linear discontinuous velocity between elements and quadratic continuous pressure between elements). A “volume of fluid” type method is used for the interface capturing, which is based on compressive control volume advection and second-order finite element methods. A force-balanced continuum surface force model is employed for the interfacial tension on unstructured meshes. The interfacial tension coefficient decomposition method is also used to deal with interfacial tension pairings between different phases. Numerical examples of benchmark tests and the dynamics of three-dimensional three-phase rising bubble, and droplet impact are presented. The results are compared with the analytical solutions and previously published experimental data, demonstrating the capability of the present method.     

Invariant discrete flows

In this paper, we investigate the evolution of joint invariants under invariant geometric flows using the theory of equivariant moving frames and the induced invariant discrete variational complex. For certain arc length preserving planar curve flows invariant under the special Euclidean group , the special linear group , and the semidirect group , we find that the induced evolution of the discrete curvature satisfies the differential‐difference mKdV, KdV, and Burgers' equations, respectively. These three equations are completely integrable, and we show that a recursion operator can be constructed by precomposing the characteristic operator of the curvature by a certain invariant difference operator. Finally, we derive the constraint for the integrability of the discrete curvature evolution to lift to the evolution of the discrete curve itself.     

Precise measurement and control of the pressure-driven flows for microfluidic systems

The pressure-driven device is designed and the flow rates of the microfluidic systems can be supplied by the pressure-driven flows, which can significantly reduce the flow-rate fluctuations coming from the pump source. For pressure-driven flows, the flow rates of the fluids can be predicted by measuring the pressure drop along a polytetrafluoroethylene (PTFE) tubing. Especially, by varying the geometrical parameters of the PTFE tubing, the predicted flow rates of the fluids are compared with the experimental measurements, and the testing precision of the pressure-driven flows can be obtained. Meanwhile, the dynamic characteristics of the open-loop and closed-loop control pressure-driven device are comparatively studied. Particularly, a proportional and integral (PI) controller is integrated with the closed-loop control pressure-driven device, and the effects of the parameters of the PI controller on the dynamic characteristics of the pressure-driven devices are mainly discussed. Most importantly, by improving the dynamic characteristics of the pressure-driven devices, precise measurement and control of the pressure-driven flows can be achieved for microfluidic systems.     

Simulating bubble dynamics in a buoyant system

Multiphase flows are critical components of many physical systems; however, numerical models of multiphase flows with large parameter gradients can be challenging. Here, two different numerical methods, volume of fluid (VOF) and smoothed particle hydrodynamics (SPH), are used to model the buoyant rise of isolated gas bubbles through quiescent fluids for a range of Bond and Reynolds numbers. The VOF is an Eulerian grid–based method, whereas the SPH is Lagrangian and mesh free. Each method has unique strengths and weaknesses, and a comparison of the two approaches as applied to multiphase phenomena has not previously been performed. The VOF and SPH simulations are compared, verified, and validated. Results using two-dimensional VOF and SPH simulations are similar to each other and are able to reproduce numerical benchmarks and experimental results for sufficiently large Morton and Reynolds numbers. It is also shown that at low Reynolds numbers, the two methods, SPH and VOF, diverge in the transient regime of the bubble rise. Regimes that require simulations capable of representing three-dimensional drag are identified as well as regimes in which results from VOF and SPH diverge.     

Existence and asymptotic properties of solutions for a nonlinear Schrödinger elliptic equation from geophysical fluid flows

In this paper, we study the existence and asymptotic behavior of radial solutions for a class of nonlinear Schrödinger elliptic equations on infinite domains describing the gyre of geophysical fluid flows. The existence theorem and asymptotic properties of radial positive solutions are established by using a new renormalization technique.