scholarly journals Euler–Lagrange Simulations of Bubble Cloud Dynamics Near a Wall

2015 ◽  
Vol 137 (4) ◽  
Author(s):  
Jingsen Ma ◽  
Chao-Tsung Hsiao ◽  
Georges L. Chahine

We present in this paper a two-way coupled Eulerian–Lagrangian model to study the dynamics of clouds of microbubbles subjected to pressure variations and the resulting pressures on a nearby rigid wall. The model simulates the two-phase medium as a continuum and solves the Navier–Stokes equations using Eulerian grids with a time and space varying density. The microbubbles are modeled as interacting singularities representing moving and oscillating spherical bubbles, following a modified Rayleigh–Plesset–Keller–Herring equation and are tracked in a Lagrangian fashion. A two-way coupling between the Euler and Lagrange components is realized through the local mixture density determined by the bubbles' volume change and motion. Using this numerical framework, simulations involving a large number of bubbles were conducted under driving pressures at different frequencies. The results show that the frequency of the driving pressure is critical in determining the overall dynamics: either a collective strongly coupled cluster behavior or nonsynchronized weaker multiple bubble oscillations. The former creates extremely high pressures with peak values orders of magnitudes higher than that of the excitation pressure. This occurs when the driving frequency matches the natural frequency of the bubble cloud. The initial distance between the bubble cloud and the wall also affects significantly the resulting pressures. A bubble cloud collapsing very close to the wall exhibits a cascading collapse, with the bubbles farthest from the wall collapsing first and the nearest ones collapsing last, thus the energy accumulates and this results in very high pressure peaks at the wall. At farther cloud distances from the wall, the bubble cloud collapses quasi-spherically with the cloud center collapsing last.

Author(s):  
Jingsen Ma ◽  
Chao-Tsung Hsiao ◽  
Georges L. Chahine

We present in this paper a two-way coupled Eulerian-Lagrangian model to study the dynamics of microbubble clouds exposed to incoming pressure waves and the resulting pressure loads on a nearby rigid wall. The model simulates the two-phase medium as a continuum and solves the N-S equations using Eulerian grids with a time and space varying density. The microbubbles are modeled as interacting spherical bubbles, which follow a modified Rayleigh-Plesset-Keller-Herring equation and are tracked in a Lagrangian fashion. A two-way coupling between the Euler and Lagrange components is realized through the local mixture density associated with the bubbles volume change and motion. Using this numerical framework, simulations involving a large number of bubbles were conducted under driving pressures of different frequencies. The results show that the frequency of the driving pressure is critical in determining the overall dynamics: either a collective strongly coupled cluster behavior or non-synchronized weaker multiple bubble oscillations. The former creates extremely high pressures with peak values orders of magnitudes higher than that of the excitation pressures. This occurs when the driving frequency matches the natural frequency of the bubble cloud. The initial distance between the bubble cloud and the wall is also critical on the resulting pressure loads. A bubble cloud collapsing very close to the wall exhibits a cascading collapse with the bubbles farthest from the wall collapsing first and the nearest ones collapsing last, thus the energy accumulates and then results in very violent pressure peaks at the wall. Farther from the wall, the bubble cloud collapses quasi spherically with the cloud center collapsing last.


2006 ◽  
Vol 4 ◽  
pp. 224-236
Author(s):  
A.S. Topolnikov

The paper is devoted to numerical modeling of Navier–Stokes equations for incompressible media in the case, when there exist gas and liquid inside the rectangular calculation region, which are separated by interphase boundary. The set of equations for incompressible liquid accounting for viscous, gravitational and surface (capillary) forces is solved by finite-difference scheme on the spaced grid, for description of interphase boundary the ideology of Level Set Method is used. By developed numerical code the set of hydrodynamic problems is solved, which describe the motion of two-phase incompressible media with interphase boundary. As a result of numerical simulation the solutions are obtained, which are in good agreement with existing analytical and experimental solutions.


Water ◽  
2019 ◽  
Vol 11 (12) ◽  
pp. 2511
Author(s):  
Jintao Liu ◽  
Di Xu ◽  
Shaohui Zhang ◽  
Meijian Bai

This paper investigates the physical processes involved in the water filling and air expelling process of a pipe with multiple air valves under water slow filling condition, and develops a fully coupledwater–air two-phase stratified numerical model for simulating the process. In this model, the Saint-Venant equations and the Vertical Average Navier–Stokes equations (VANS) are respectively applied to describe the water and air in pipe, and the air valve model is introduced into the VANS equations of air as the source term. The finite-volume method and implicit dual time-stepping method (IDTS) with two-order accuracy are simultaneously used to solve this numerical model to realize the full coupling between water and air movement. Then, the model is validated by using the experimental data of the pressure evolution in pipe and the air velocity evolution of air valves, which respectively characterize the water filling and air expelling process. The results show that the model performs well in capturing the physical processes, and a reasonable agreement is obtained between numerical and experimental results. This agreement demonstrates that the proposed model in this paper offers a practical method for simulating water filling and air expelling process in a pipe with multiple air valves under water slow filling condition.


Author(s):  
Marion W. Vance ◽  
Kyle D. Squires

An approach to parallel solution of an Eulerian-Lagrangian model of dilute gas-solid flows is presented. Using Lagrangian treatments for the dispersed phase, one of the principal computational challenges arises in models in which inter-particle interactions are taken into account. Deterministic treatment of particle-particle collisions in the present work pose the most computationally intensive aspect of the simulation. Simple searches lead to algorithms whose cost is O(N2p) where Np is the particle population. The approach developed in the current effort is based on localizing collision detection neighborhoods using a cell-index method and spatially distributing those neighborhoods for parallel solution. The method is evaluated using simulations of the gas-solid turbulent flow in a vertical channel. The instantaneous position and the velocity of any particle is obtained by solving the equation of motion for a small rigid sphere assuming that the resulting force induced by the fluid reduces to the drag contribution. Binary particle collisions without energy dissipation or inter-particle friction are considered. The carrier flow is computed using Large Eddy Simulation of the incompressible Navier-Stokes equations. The entire dispersed-phase population is partitioned via static spatial decomposition of the domain to maximize parallel efficiency. Simulations on small numbers of distributed memory processors show linear speedup in processing of the collision detection step and nearly optimal reductions in simulation time for the entire solution.


2018 ◽  
Vol 8 (9) ◽  
pp. 1413 ◽  
Author(s):  
Dan Yao ◽  
Kwongi Lee ◽  
Minho Ha ◽  
Cheolung Cheong ◽  
Inhiug Lee

A new pump, called the hybrid airlift-jet pump, is developed by reinforcing the advantages and minimizing the demerits of airlift and jet pumps. First, a basic design of the hybrid airlift-jet pump is schematically presented. Subsequently, its performance characteristics are numerically investigated by varying the operating conditions of the airlift and jet parts in the hybrid pump. The compressible unsteady Reynolds-averaged Navier-Stokes equations, combined with the homogeneous mixture model for multiphase flow, are used as the governing equations for the two-phase flow in the hybrid pump. The pressure-based methods combined with the Pressure-Implicit with Splitting of Operators (PISO) algorithm are used as the computational fluid dynamics techniques. The validity of the present numerical methods is confirmed by comparing the predicted mass flow rate with the measured ones. In total, 18 simulation cases that are designed to represent the various operating conditions of the hybrid pump are investigated: eight of these cases belong to the operating conditions of only the jet part with different air and water inlet boundary conditions, and the remaining ten cases belong to the operating conditions of both the airlift and jet parts with different air and water inlet boundary conditions. The mass flow rate and the efficiency are compared for each case. For further investigation into the detailed flow characteristics, the pressure and velocity distributions of the mixture in a primary pipe are compared. Furthermore, a periodic fluctuation of the water flow in the mass flow rate is found and analyzed. Our results show that the performance of the jet or airlift pump can be enhanced by combining the operating principles of two pumps into the hybrid airlift-jet pump, newly proposed in the present study.


2007 ◽  
Vol 18 (04) ◽  
pp. 536-545 ◽  
Author(s):  
NAOKI TAKADA ◽  
AKIO TOMIYAMA

For interface-tracking simulation of two-phase flows in various micro-fluidics devices, we examined the applicability of two versions of computational fluid dynamics method, NS-PFM, combining Navier-Stokes equations with phase-field modeling for interface based on the van der Waals-Cahn-Hilliard free-energy theory. Through the numerical simulations, the following major findings were obtained: (1) The first version of NS-PFM gives good predictions of interfacial shapes and motions in an incompressible, isothermal two-phase fluid with high density ratio on solid surface with heterogeneous wettability. (2) The second version successfully captures liquid-vapor motions with heat and mass transfer across interfaces in phase change of a non-ideal fluid around the critical point.


Author(s):  
Christian Rohde ◽  
Lars von Wolff

We consider the incompressible flow of two immiscible fluids in the presence of a solid phase that undergoes changes in time due to precipitation and dissolution effects. Based on a seminal sharp interface model a phase-field approach is suggested that couples the Navier–Stokes equations and the solid’s ion concentration transport equation with the Cahn–Hilliard evolution for the phase fields. The model is shown to preserve the fundamental conservation constraints and to obey the second law of thermodynamics for a novel free energy formulation. An extended analysis for vanishing interfacial width reveals that in this limit the sharp interface model is recovered, including all relevant transmission conditions. Notably, the new phase-field model is able to realize Navier-slip conditions for solid–fluid interfaces in the limit.


Water ◽  
2021 ◽  
Vol 13 (21) ◽  
pp. 3092
Author(s):  
Lourenço Sassetti Mendes ◽  
Javier L. Lara ◽  
Maria Teresa Viseu

Spillway design is key to the effective and safe operation of dams. Typically, the flow is characterized by high velocity, high levels of turbulence, and aeration. In the last two decades, advances in computational fluid dynamics (CFD) made available several numerical tools to aid hydraulic structures engineers. The most frequent approach is to solve the Reynolds-averaged Navier–Stokes equations using an Euler type model combined with the volume-of-fluid (VoF) method. Regardless of a few applications, the complete two-phase Euler is still considered to demand exorbitant computational resources. An assessment is performed in a spillway offset aerator, comparing the two-phase volume-of-fluid (TPVoF) with the complete two-phase Euler (CTPE). Both models are included in the OpenFOAM® toolbox. As expected, the TPVoF results depend highly on the mesh, not showing convergence in the maximum chute bottom pressure and the lower-nappe aeration, tending to null aeration as resolution increases. The CTPE combined with the k–ω SST Sato turbulence model exhibits the most accurate results and mesh convergence in the lower-nappe aeration. Surprisingly, intermediate mesh resolutions are sufficient to surpass the TPVoF performance with reasonable calculation efforts. Moreover, compressibility, flow bulking, and several entrained air effects in the flow are comprehended. Despite not reproducing all aspects of the flow with acceptable accuracy, the complete two-phase Euler demonstrated an efficient cost-benefit performance and high value in spillway aerated flows. Nonetheless, further developments are expected to enhance the efficiency and stability of this model.


2009 ◽  
Vol 74 ◽  
pp. 139-142
Author(s):  
Ting Ye ◽  
Hua Li

A modeling of two-phase system is presented for investigation of the cell motion and deformation in the microchannel subject to the mechanical and electrical coupled forces. In order to evaluate the mechanical force developed by cell membrane, it is treated as an incompressible and elastic shell with uniform thickness capable of shearing and bending deformation. Due to the irregular and complex cell configuration after deformation, the Maxwell stress tensor (MST) method is successfully employed to analyze the dielectrophoretic force. The modified particle binary level set (MPBLS) method is presented to accurately track the moving interface between the two phases, which is vital for a modeling of two-phase system. Afterwards the modified SIMPLER coupled with SIMPLEC is used to numerically solve the incompressible Navier-Stokes equations governing the entire flow field. On basis of the series of methods, the motion and deformation of red blood cell (RBC) in the microchannel under the mechanical and electrical forces are simulated to demonstrate the deformation process and the moving trajectory of RBC. The present study is not only of great value for deeper understanding of some diseases caused by cell abnormality, but also of practical significance for cell manipulation and separation.


Sign in / Sign up

Export Citation Format

Share Document