Gradient-augmented hybrid interface capturing method for incompressible two-phase flow

Motivated by inconveniences of present hybrid methods,a gradient-augmented hybrid interface capturing method(GAHM) is presented for incompressible two-phase flow.A front tracking method(FTM) is used as the skeleton of the GAHM for low mass loss and resources.Smooth eulerian level set values are calculated from the FTM interface,and are used for a local interface reconstruction.The reconstruction avoids marker particle redistribution and enables an automatic treatment of interfacial topology change.The cubic Hermit interpolation is employed in all steps of the GAHM to capture subgrid structures within a single spacial cell.The performance of the GAHM is carefully evaluated in a benchmark test.Results show significant improvements of mass loss,clear subgrid structures,highly accurate derivatives(normals and curvatures) and low cost.The GAHM is further coupled with an incompressible multiphase flow solver,Super CE/SE,for more complex and practical applications.The updated solver is evaluated through comparison with an early droplet research.     

A hybrid scheme for computing incompressible two-phase flows

We propose a hybrid scheme for computations of incompressible two-phase flows. The incompressible constraint has been replaced by a pressure Poisson-like equation and then the pressure is updated by the modified marker and cell method. Meanwhile, the moment equations in the incompressible Navier-Stokes equations are solved by our semidiscrete Hermite central-upwind scheme, and the interface between the two fluids is considered to be continuous and is described implicitly as the 0.5 level set of a smooth function being a smeared out Heaviside function. It is here named the hybrid scheme. Some numerical experiments are successfully carried out, which verify the desired efficiency and accuracy of our hybrid scheme.     

Anti-diffusion method for interface steepening in two-phase incompressible flow

In this paper, we present a method for obtaining sharp interfaces in two-phase incompressible flows by an anti-diffusion correction, that is applicable in a straight-forward fashion for the improvement of two-phase flow solution schemes typically employed in practical applications. The underlying discretization is based on the volume-of-fluid (VOF) interface-capturing method on unstructured meshes. The key idea is to steepen the interface, independently of the underlying volume-fraction transport equation, by solving a diffusion equation with reverse time, i.e. an anti-diffusion equation, after each advection time step of the volume fraction. As the solution of the anti-diffusion equation requires regularization, a limiter based on the directional derivative is developed for calculating the gradient of the volume fraction. This limiter ensures the boundedness of the volume fraction. In order to control the amount of anti-diffusion introduced by the correction algorithm we propose a suitable stopping criterion for interface steepening. The formulation of the limiter and the algorithm for solving the anti-diffusion equation are applicable to 3-dimensional unstructured meshes. Validation computations are performed for passive advection of an interface, for 2-dimensional and 3-dimensional rising-bubbles, and for a rising drop in a periodically constricted channel. The results demonstrate that sharp interfaces can be recovered reliably. They show that the accuracy is similar to or even better than that of level-set methods using comparable discretizations for the flow and the level-set evolution. Also, we observe a good agreement with experimental results for the rising drop where proper interface evolution requires accurate mass conservation.     

A conservative phase field method for solving incompressible two-phase flows

In this paper a conservative phase-field method based on the work of Sun and Beckermann [Y. Sun, C. Beckermann, Sharp interface tracking using the phase-field equation, J. Comput. Phys. 220 (2007) 626–653] for solving the two- and three-dimensional two-phase incompressible Navier–Stokes equations is proposed. The present method can preserve the total mass as the Cahn–Hilliard equation, but the calculation and implementation are much simpler than that. The dispersion-relation-preserving schemes are utilized for the advection terms while the Helmholtz smoother is applied to compute the surface-tension force term. To verify the proposed method, several benchmarks are examined and shown to have good agreements with previous results. It also shows that the satisfactions of mass conservations are guaranteed.     

A numerical method for incompressible viscous flow simulation

We describe a numerical scheme for computing time-dependent solutions of the incompressible Navier-Stokes equations in the primitive variable formulation. This scheme uses finite elements for the space discretization and operator splitting techniques for the time discretization. The resulting discrete equations are solved using specialized nonlinear optimization algorithms that are computationally efficient and have modest storage requirements. The basic numerical kernel is the preconditioned conjugate gradient method for symmetric, positive-definite, sparse matrix systems, which can be efficiently implemented on the architectures of vector and parallel processing supercomputers.     

A continuous Lagrangian sensitivity equation method for incompressible flow

A continuous Lagrangian sensitivity equation method (CLSEM) is presented as a cost effective alternative to the continuous (Eulerian) sensitivity equation method (CESEM) in the case of shape parameters. Boundary conditions for the CLSEM are simpler than those of the CESEM. However a mapping must be introduced to relate the undeformed and deformed configurations thus making the PDEs more complicated. We propose the use of pseudo-elasticity equations to provide a general framework to generate this mapping for unstructured meshes on complex geometries. The methodology is presented in details for the incompressible Navier–Stokes and sensitivity equations in variational form. The PDEs are solved with an adaptive FEM. Sensitivity data obtained with both approaches for a flow around a NACA 4512 are used to obtain estimates of flows around nearby geometries. Results indicate that the CLSEM produces significant improvements in terms of both accuracy and CPU time.     

A balanced force refined level set grid method for two-phase flows on unstructured flow solver grids

A balanced force refined level set grid method for two-phase flows on structured and unstructured flow solver grids is presented. To accurately track the phase interface location, an auxiliary, high-resolution equidistant Cartesian grid is introduced. In conjunction with a dual-layer narrow band approach, this refined level set grid method allows for parallel, efficient grid convergence and error estimation studies of the interface tracking method. The Navier–Stokes equations are solved on an unstructured flow solver grid with a novel balanced force algorithm for level set methods based on the recently proposed method by Francois et al. [M.M. Francois, S.J. Cummins, E.D. Dendy, D.B. Kothe, J.M. Sicilian, M.W. Williams, A balanced-force algorithm for continuous and sharp interfacial surface tension models within a volume tracking framework, J. Comput. Phys. 213 (2006) 141–173] for volume of fluid methods on structured grids. To minimize spurious currents, a second order converging curvature evaluation technique for level set methods is presented. The results of several different test cases demonstrate the effectiveness of the proposed method, showing good mass conservation properties and second order converging spurious current magnitudes.     

A stabilized MLPG method for steady state incompressible fluid flow simulation

In this paper, the meshless local Petrov–Galerkin (MLPG) method is extended to solve the incompressible fluid flow problems. The streamline upwind Petrov–Galerkin (SUPG) method is applied to overcome oscillations in convection-dominated problems, and the pressure-stabilizing Petrov–Galerkin (PSPG) method is applied to satisfy the so-called Babuška–Brezzi condition. The same stabilization parameter τ(τ_SUPG = τ_PSPG) is used in the present method. The circle domain of support, linear basis, and fourth-order spline weight function are applied to compute the shape function, and Bubnov–Galerkin method is applied to discretize the PDEs. The lid-driven cavity flow, backward facing step flow and natural convection in the square cavity are applied to validate the accuracy and feasibility of the present method. The results show that the stability of the present method is very good and convergent solutions can be obtained at high Reynolds number. The results of the present method are in good agreement with the classical results. It also seems that the present method (which is a truly meshless) is very promising in dealing with the convection- dominated problems.     

Force-coupling method for particulate two-phase flow: Stokes flow

In this paper we describe a force-coupling method for particle dynamics in fluid flows. The general principles of the model are described and it is tested on three different Stokes flow problems; a single isolated sphere, a pair of otherwise isolated spheres, and a single sphere in a channel. For sphere to sphere or sphere to wall distances larger than 1/4 of the sphere radius the force-coupling results compared well with analytical and accurate numerical values. For smaller distances the results agree qualitatively, but lubrication effects are not included and this leads to a quantitative discrepancy.     

A lattice Boltzmann method for incompressible two-phase flows on partial wetting surface with large density ratio

This paper reports a new numerical scheme of the lattice Boltzmann method for calculating liquid droplet behaviour on particle wetting surfaces typically for the system of liquid–gas of a large density ratio. The method combines the existing models of Inamuro et al. [T. Inamuro, T. Ogata, S. Tajima, N. Konishi, A lattice Boltzmann method for incompressible two-phase flows with large density differences, J. Comput. Phys. 198 (2004) 628–644] and Briant et al. [A.J. Briant, P. Papatzacos, J.M. Yeomans, Lattice Boltzmann simulations of contact line motion in a liquid–gas system, Philos. Trans. Roy. Soc. London A 360 (2002) 485–495; A.J. Briant, A.J. Wagner, J.M. Yeomans, Lattice Boltzmann simulations of contact line motion: I. Liquid–gas systems. Phys. Rev. E 69 (2004) 031602; A.J. Briant, J.M. Yeomans, Lattice Boltzmann simulations of contact line motion: II. Binary fluids, Phys. Rev. E 69 (2004) 031603] and has developed novel treatment for partial wetting boundaries which involve droplets spreading on a hydrophobic surface combined with the surface of relative low contact angles and strips of relative high contact angles. The interaction between the fluid–fluid interface and the partial wetting wall has been typically considered. Applying the current method, the dynamics of liquid drops on uniform and heterogeneous wetting walls are simulated numerically. The results of the simulation agree well with those of theoretical prediction and show that the present LBM can be used as a reliable way to study fluidic control on heterogeneous surfaces and other wetting related subjects.     

Three-dimensional multi-relaxation-time lattice Boltzmann front-tracking method for two-phase flow

We developed a three-dimensional multi-relaxation-time lattice Boltzmann method for incompressible and immiscible two-phase flow by coupling with a front-tracking technique. The flow field was simulated by using an Eulerian grid, an adaptive unstructured triangular Lagrangian grid was applied to track explicitly the motion of the two-fluid interface, and an indicator function was introduced to update accurately the fluid properties. The surface tension was computed directly on a triangular Lagrangian grid, and then the surface tension was distributed to the background Eulerian grid. Three benchmarks of two-phase flow, including the Laplace law for a stationary drop, the oscillation of a three-dimensional ellipsoidal drop,and the drop deformation in a shear flow, were simulated to validate the present model.     

A consistent and stabilized continuous/discontinuous Galerkin method for fourth-order incompressible flow problems

This paper presents a new consistent and stabilized finite-element formulation for fourth-order incompressible flow problems. The formulation is based on the C⁰-interior penalty method, the Galerkin least-square (GLS) scheme, which assures that the formulation is weakly coercive for spaces that fail to satisfy the inf-sup condition, and considers discontinuous pressure interpolations. A stability analysis through a lemma establishes that the proposed formulation satisfies the inf-sup condition, thus confirming the robustness of the method. This lemma indicates that, at the element level, there exists an optimal or quasi-optimal GLS stability parameter that depends on the polynomial degree used to interpolate the velocity and pressure fields, the geometry of the finite element, and the fluid viscosity term. Numerical experiments are carried out to illustrate the ability of the formulation to deal with arbitrary interpolations for velocity and pressure, and to stabilize large pressure gradients.     

HyPAM: A hybrid continuum-particle model for incompressible free-surface flows

Three major issues associated with numerical simulations of complex free-surface flows, viz. interface tracking, fragmentation and large physical jumps, are addressed by a new hybrid continuum-particle model (HyPAM). The new model consists of three parts: (1) the Polygonal Area Mapping method [Q. Zhang, P.L.-F. Liu, A new interface tracking method: the polygonal area mapping method, J. Comput. Phys. 227(8) (2008) 4063-4088]; (2) a new algorithm that decomposes the interested (water) phase into a continuum zone, a buffer zone and a particle zone, based on material topology and graph theory; (3) a ‘passive-response’ assumption, in which the air phase is assumed to respond passively to the continuum part of the water phase. The incompressible inviscid Euler equations and the equations describing the free fall of rigid bodies are used as the governing equations for the continuum-buffer zone and the particle zone, respectively, and separately. A number of examples, including water droplet impact, solitary wave propagation, and dam-break problems, are simulated for the illustration and validation of HyPAM. It is shown that HyPAM is more accurate and versatile than a continuum-based Volume-of-Fluid model. One major contribution of this work is the single-phase decomposition algorithm, useful for many other hybrid formulations. Neglecting surface tension, viscosity and particle interactions, HyPAM is currently limited to mildly-fragmented free-surface flows with high Reynolds and Weber numbers.     

A hybrid numerical method for interfacial fluid flow with soluble surfactant

We address a significant difficulty in the numerical computation of fluid interfaces with soluble surfactant that occurs in the physically representative limit of large bulk Peclet number Pe. At the high values of Pe in typical fluid-surfactant systems, there is a transition layer near the interface in which the surfactant concentration varies rapidly, and large gradients at the interface must be resolved accurately to evaluate the exchange of surfactant between the interface and bulk flow. We use the slenderness of the layer to develop a fast and accurate ‘hybrid’ numerical method that incorporates a separate, singular perturbation analysis of the dynamics in the transition layer into a full numerical solution of the interfacial free boundary problem. The accuracy and efficiency of the method is assessed by comparison with a more ‘traditional’ numerical approach that uses finite differences on a curvilinear coordinate system exterior to the bubble, without the separate transition layer reduction. The traditional method implemented here features a novel fast calculation of fluid velocity off the interface.     

A lattice Boltzmann method for a two-phase flow containing solid bodies with viscoelastic membranes

A lattice Boltzmann method (LBM) for two-phase flows containing solid bodies with viscoelastic membranes is proposed. The method is based on the two-phase LBM, in which one phase is regarded as the solid phase. In the present model, the membrane is assumed to be composed of identical particles that are connected to their neighboring particles by elastic springs to take account of stretching and compression effects. The method is applied to two representative problems, namely the behavior of a viscoelastic body under shear flow and the motion of a viscoelastic body in a Poiseuille flow. Tank-tread motion and axial migration, which are both characteristic of the motion of viscoelastic bodies, are simulated by using the method. These results indicate that the method is capable of simulating the complex behavior of viscoelastic bodies in capillaries, such as the motion of red blood cells in blood flows.     

Nested Cartesian grid method in incompressible viscous fluid flow

In this work, the local grid refinement procedure is focused by using a nested Cartesian grid formulation. The method is developed for simulating unsteady viscous incompressible flows with complex immersed boundaries. A finite-volume formulation based on globally second-order accurate central-difference schemes is adopted here in conjunction with a two-step fractional-step procedure. The key aspects that needed to be considered in developing such a nested grid solver are proper imposition of interface conditions on the nested-block boundaries, and accurate discretization of the governing equations in cells that are with block-interface as a control-surface. The interpolation procedure adopted in the study allows systematic development of a discretization scheme that preserves global second-order spatial accuracy of the underlying solver, and as a result high efficiency/accuracy nested grid discretization method is developed. Herein the proposed nested grid method has been widely tested through effective simulation of four different classes of unsteady incompressible viscous flows, thereby demonstrating its performance in the solution of various complex flow–structure interactions. The numerical examples include a lid-driven cavity flow and Pearson vortex problems, flow past a circular cylinder symmetrically installed in a channel, flow past an elliptic cylinder at an angle of attack, and flow past two tandem circular cylinders of unequal diameters. For the numerical simulations of flows past bluff bodies an immersed boundary (IB) method has been implemented in which the solid object is represented by a distributed body force in the Navier–Stokes equations. The main advantages of the implemented immersed boundary method are that the simulations could be performed on a regular Cartesian grid and applied to multiple nested-block (Cartesian) structured grids without any difficulty. Through the numerical experiments the strength of the solver in effectively/accurately simulating various complex flows past different forms of immersed boundaries is extensively demonstrated, in which the nested Cartesian grid method was suitably combined together with the fractional-step algorithm to speed up the solution procedure.     

Transient adaptivity applied to two-phase incompressible flows

An anisotropic adaptation process is applied to a three-dimensional incompressible two-phase flow solver. The solver uses a level set/finite element method on unstructured tetrahedral meshes. We show how the level set function can be used to build an anisotropic mesh with good properties. Some computations with a strong transient character and large densities ratios (1/1000) are presented. We show that the efficiency of the computations can be deeply enhanced by mesh adaptations.     

A high-order incompressible flow solver with WENO

A numerical method for solving the incompressible Navier–Stokes equations with a 5th-order weighted essentially non-oscillatory (WENO) scheme is presented. The method is not based on artificial compressibility and is free of tunable parameters such as the artificial compressibility parameter. The method makes use of the fractional-step method in conjunction with the low-dissipation and low-dispersion Runge–Kutta (LDDRK) scheme to improve temporal accuracy of the scheme. The use of a WENO scheme makes it possible to obtain stable solutions for discontinuous initial data and resolve difficult applications with strong shear such as Kelvin–Helmholtz instability or turbulence. Good convergence rate is established for the velocity variables and numerical solutions of the present method compare well with exact solutions and other numerical results.     

A lattice Boltzmann method for immiscible multiphase flow simulations using the level set method

We consider the lattice Boltzmann method for immiscible multiphase flow simulations. Classical lattice Boltzmann methods for this problem, e.g. the colour gradient method or the free energy approach, can only be applied when density and viscosity ratios are small. Moreover, they use additional fields defined on the whole domain to describe the different phases and model phase separation by special interactions at each node. In contrast, our approach simulates the flow using a single field and separates the fluid phases by a free moving interface. The scheme is based on the lattice Boltzmann method and uses the level set method to compute the evolution of the interface. To couple the fluid phases, we develop new boundary conditions which realise the macroscopic jump conditions at the interface and incorporate surface tension in the lattice Boltzmann framework. Various simulations are presented to validate the numerical scheme, e.g. two-phase channel flows, the Young–Laplace law for a bubble and viscous fingering in a Hele-Shaw cell. The results show that the method is feasible over a wide range of density and viscosity differences.     

耦合不可压流场输运方程的格子Boltzmann方法研究

基于格子Boltzmann方法,提出了求解耦合不可压缩流场输运方程的一种改进数值方法. 该方法使用格子Boltzmann方法求解流场方程,并根据流场格子模型的密度分布函数构建了输运方程的二阶离散格式. 通过二维平板通道流场输运系统验证了该方法的有效性.数值结果表明,该方法可以有效地减少计算过程中出现的非物理耗散, 并克服了传统模型所需巨大存储量的缺点.