An Implementation of MAC Grid-Based IIM-Stokes Solver for Incompressible Two-Phase Flows

2011 ◽  
Vol 10 (5) ◽  
pp. 1333-1362 ◽  
Author(s):  
Zhijun Tan ◽  
K. M. Lim ◽  
B. C. Khoo

AbstractIn this paper, a novel implementation of immersed interface method combined with Stokes solver on a MAC staggered grid for solving the steady two-fluid Stokes equations with interfaces. The velocity components along the interface are introduced as two augmented variables and the resulting augmented equation is then solved by the GMRES method. The augmented variables and /or the forces are related to the jumps in pressure and the jumps in the derivatives of both pressure and velocity, and are interpolated using cubic splines and are then applied to the fluid through the jump conditions. The Stokes equations are discretized on a staggered Cartesian grid via a second order finite difference method and solved by the conjugate gradient Uzawa-type method. The numerical results show that the overall scheme is second order accurate. The major advantages of the present IIM-Stokes solver are the efficiency and flexibility in terms of types of fluid flow and different boundary conditions. The proposed method avoids solution of the pressure Poisson equation, and comparisons are made to show the advantages of time savings by the present method. The generalized two-phase Stokes solver with correction terms has also been applied to incompressible two-phase Navier-Stokes flow.

2014 ◽  
Vol 580-583 ◽  
pp. 2890-2893
Author(s):  
Ru Quan Liang ◽  
Zhi Hui Zhang ◽  
Tai Yin Gao ◽  
Fu Sheng Yan

In this paper, the vertical vibration influence on the surface velocity of a 5cSt silicone oil liquid bridge has been investigated numerically. The Navier-Stokes equations coupled with the energy conservation equation are solved on a staggered grid, and the two-phase surface is captured by using the mass conserving level set method. The present results indicate that the axial and radial surface velocities of the liquid bridge are suppressed by the external vertical vibration.


Author(s):  
Hessam Babaee ◽  
Sumanta Acharya

An accurate and efficient finite difference method for solving the three dimensional incompressible Navier-Stokes equations on curvilinear grids is developed. The semi-staggered grid layout has been used in which all three components of velocity are stored on the corner vertices of the cell facilitating a consistent discretization of the momentum equations as the boundaries are approached. Pressure is stored at the cell-center, resulting in the exact satisfaction the discrete continuity. The diffusive terms are discretized using a second-order central finite difference. A third-order biased upwind scheme is used to discretize the convective terms. The momentum equations are integrated in time using a semi-implicit fractional step methodology. The convective and diffusive terms are advanced in time using the second-order Adams-Bashforth method and Crank-Nicolson method respectively. The Pressure-Poisson is discretized in a similar approach to the staggered gird layout and thus leading to the elimination of the spurious pressure eigen-modes. The validity of the method is demonstrated by two standard benchmark problems. The flow in driven cavity is used to show the second-order spatial convergence on an intentionally distorted grid. Finally, the results for flow past a cylinder for several Reynolds numbers in the range of 50–150 are compared with the existing experimental data in the literature.


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.


Analysis ◽  
2020 ◽  
Vol 40 (3) ◽  
pp. 127-150
Author(s):  
Tania Biswas ◽  
Sheetal Dharmatti ◽  
Manil T. Mohan

AbstractIn this paper, we formulate a distributed optimal control problem related to the evolution of two isothermal, incompressible, immiscible fluids in a two-dimensional bounded domain. The distributed optimal control problem is framed as the minimization of a suitable cost functional subject to the controlled nonlocal Cahn–Hilliard–Navier–Stokes equations. We describe the first order necessary conditions of optimality via the Pontryagin minimum principle and prove second order necessary and sufficient conditions of optimality for the problem.


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.


1993 ◽  
Vol 115 (4) ◽  
pp. 678-685 ◽  
Author(s):  
M. J. Braun ◽  
F. K. Choy ◽  
Y. M. Zhou

The flow in a hydrostatic pocket is described by a mathematical model that uses the Navier-Stokes equations written in terms of the primary variables, u, v, and p. Using the conservative formulation, a finite difference method is applied through a staggered grid. The power law scheme is applied in the treatment of the convective terms for this highly recirculating flow. The discussion pertaining to the convergence of the numerical scheme and the computational error, shows that the strict convergence criteria applied to both velocities and pressure were successfully statisfied. The numerical model is applied in a parametric mode to the study of the velocities, the pressure patterns, and shear forces that characterize the flow in a square (α = 1), deep (α>1), and shallow (α≪1) hydrostatic pocket. The effects of the variation of the location and angle of the hydrostatic jet are also investigated.


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.


Sign in / Sign up

Export Citation Format

Share Document