A direct reinitialization approach of level‐set/splitting finite element method for simulating incompressible two‐phase flows

Computation of a moving interface by the level‐set (LS) method typically requires reinitialization of LS function. An inaccurate execution of reinitialization results in incorrect free surface capturing and thus errors such as mass gain/loss so that an accurate and robust reinitialization process in the LS method is essential for the simulation of free surface flows. In the present study, we pursue further development of the reinitialization process, which directly corrects the LS function after advection is carried out by using the normal vector to the interface instead of solving the reinitialization equation of hyperbolic type. The Taylor–Galerkin method is adopted to discretize the advection equation of the LS function and the P1P1 splitting finite element method is applied to solve the Navier–Stokes equation. The proposed algorithm is validated with the well‐known benchmark problems, i.e. stretching of a circular fluid element, time‐reversed single‐vortex, solitary wave propagation, broken dam flow and filling of a container. The simulation results of these flows are in good agreement with previously existing experimental and numerical results. Copyright © 2010 John Wiley & Sons, Ltd.     

A Q2Q1 finite element/level‐set method for simulating two‐phase flows with surface tension

A Q2Q1 (quadratic velocity/linear pressure) finite element/level‐set method was proposed for simulating incompressible two‐phase flows with surface tension. The Navier–Stokes equations were solved using the Q2Q1 integrated FEM, and the level‐set variable was linearly interpolated using a ‘pseudo’ Q2Q1 finite element when calculating the density and viscosity of a fluid to avoid an unbounded density/viscosity. The advection of the level‐set function was calculated through the Taylor–Galerkin method, and the direct approach method is employed for reinitialization. The proposed method was tested by solving several benchmark problems including rising bubbles exhibiting a large density difference and the surface tension effect. The numerical results of the rising bubbles were compared with the existing results to validate the benchmark quantities such as the centroid, circularity, and rising velocity. Furthermore, we focused our attention mainly on mass conservation and time‐step. We observed that the present method represented a convergence rate between 1.0 and 1.5 orders in terms of mass conservation and provided more stable solutions even when using a larger time‐step than the critical time‐step that was imposed because of the explicit treatment of surface tension. Copyright © 2011 John Wiley & Sons, Ltd.     

Two‐level finite element method with a stabilizing subgrid for the incompressible Navier–Stokes equations

We consider the Galerkin finite element method for the incompressible Navier–Stokes equations in two dimensions. The domain is discretized into a set of regular triangular elements and the finite‐dimensional spaces employed consist of piecewise continuous linear interpolants enriched with the residual‐free bubble functions. To find the bubble part of the solution, a two‐level finite element method with a stabilizing subgrid of a single node is described, and its application to the Navier–Stokes equation is displayed. Numerical approximations employing the proposed algorithm are presented for three benchmark problems. The results show that the proper choice of the subgrid node is crucial in obtaining stable and accurate numerical approximations consistent with the physical configuration of the problem at a cheap computational cost. Copyright © 2008 John Wiley & Sons, Ltd.     

An operator splitting scheme for the stream‐function formulation of unsteady Navier–Stokes equations

A fictitious time is introduced into the unsteady equation of the stream function rendering it into a higher‐order ultra‐parabolic equation. The convergence with respect to the fictitious time (we call the latter ‘internal iterations’) allows one to obtain fully implicit nonlinear scheme in full time steps for the physical‐time variable. For particular choice of the artificial time increment, the scheme in full time steps is of second‐order of approximation. For the solution of the internal iteration, a fractional‐step scheme is proposed based on the splitting of the combination of the Laplace, bi‐harmonic and advection operators. A judicious choice for the time staggering of the different parts of the nonlinear advective terms allows us to prove that the internal iterations are unconditionally stable and convergent. We assess the number of operations needed per time step and show computational effectiveness of the proposed scheme. We prove that when the internal iterations converge, the scheme is second‐order in physical time and space, nonlinear, implicit and absolutely stable. The performance of the scheme is demonstrated for the flow created by oscillatory motion of the lid of a square cavity. All theoretical findings are demonstrated practically. Copyright © 2006 John Wiley & Sons, Ltd.     

Simulation of single‐ and two‐phase flows on sliding unstructured meshes using finite volume method

The paper presents an efficient finite volume method for unstructured grids with rotating sliding parts composed of arbitrary polyhedral elements for both single‐ and two‐phase flows. Mathematical model used in computations is based on the ensemble averaged conservation equations. These equations are solved for each phase and in case of single‐phase flow reduce to the transient Reynolds‐averaged Navier–Stokes (TRANS) equations. Transient flow induced by rotating impellers is thus resolved in time. The use of unstructured grids allows an easy and flexible meshing for the entire flow domain. Polyhedral cell volumes are created on the arbitrary mesh interface placed between rotating and static parts. Cells within the rotating parts move each time step and the new faces are created on the arbitrary interfaces only, while the rest of the domain remain ‘topologically’ unchanged. Implicit discretization scheme allows a wide range of time‐step sizes, which further reduce the computational effort. Special attention is given to the interpolation practices used for the reconstruction of the face quantities. Mass fluxes are recalculated at the beginning of each time step by using an interpolation scheme, which enhances the coupling between the pressure and velocity fields. The model has been implemented into the commercially available CFD code AVL SWIFT (AVL AST, SWIFT Manual 3.1, AVL List GmbH, Graz, Austria, 2002). Single‐phase flow in a mixing vessel stirred by a six‐bladed Rushton‐type turbine and two‐phase flow in aerated stirred vessel with the four‐blade Rushton impeller are simulated. The results are compared with the available experimental data, and good agreement is observed. The proposed algorithm is proved to be both stable and accurate for single‐phase as well as for the two‐phase flows calculations. Copyright 2004 John Wiley & Sons, Ltd.     

Fitted finite element discretization of two‐phase Stokes flow

We propose a novel fitted finite element method for two‐phase Stokes flow problems that uses piecewise linear finite elements to approximate the moving interface. The method can be shown to be unconditionally stable. Moreover, spherical stationary solutions are captured exactly by the numerical approximation. In addition, the meshes describing the discrete interface in general do not deteriorate in time, which means that in numerical simulations, a smoothing or a remeshing of the interface mesh is not necessary. We present several numerical experiments for our numerical method, which demonstrate the accuracy and robustness of the proposed algorithm. Copyright © 2016 John Wiley & Sons, Ltd.     

Spurious currents in finite element based level set methods for two‐phase flow

A study of spurious currents in continuous finite element based simulations of the incompressible Navier–Stokes equations for two‐phase flows is presented on the basis of computations on a circular drop in equilibrium. The conservative and the standard level set methods are used. It is shown that a sharp surface tension force, expressed as a line integral along the interface, can give rise to large spurious currents and oscillations in the pressure that do not decrease with mesh refinement. If instead a regularized surface tension representation is used, exact force balance at the interface is possible, both for a fully coupled discretization approach and for a fractional step projection method. However, the numerical curvature calculation introduces errors that cause spurious currents. Different ways to extend the curvature from the interface to the whole domain are discussed and investigated. The impact of using different finite element spaces and stabilization methods is also considered. Copyright © 2011 John Wiley & Sons, Ltd.     

A method for computing compressible,highly stratified flows in astrophysics based on operator splitting

A new numerical approach based on consistent operator splitting is presented for computing compressible, highly stratified flows in astrophysics. The algorithm is particularly designed to search for steady or almost steady solutions for the time-dependent Navier–Stokes equations, describing viscous flow under the influence of a strong gravitational field. The algorithm proposed is multidimensional and works in Cartesian, cylindrical or spherical co-ordinates. It uses a second-order finite volume scheme with third-order upwinding and a second-order time discretization. An adaptive time step control and monotonic multilevel grid distribution has been incorporated to speed up convergence. This method has been incorporated into a hydrodynamical code by which, for the first time, for two-dimensional models the dynamics of the boundary layer in the accretion disk around a compact star could be computed over the whole viscous time scale. © 1998 John Wiley & Sons, Ltd.     

A semi‐Lagrangian level set method for incompressible Navier–Stokes equations with free surface

In this paper, we formulate a level set method in the framework of finite elements‐semi‐Lagrangian methods to compute the solution of the incompressible Navier–Stokes equations with free surface. In our formulation, we use a quasi‐monotone semi‐Lagrangian scheme, which is both unconditionally stable and essentially non oscillatory, to compute the advective terms in the Navier–Stokes equations, the transport equation and the equation of the reinitialization stage for the level set function. The method we propose is quite robust and flexible with regard to the mesh and the geometry of the domain, as well as the magnitude of the Reynolds number. We illustrate the performance of the method in several examples, which range from a benchmark problem to test the volume conservation property of the method to the flow past a NACA0012 foil at high Reynolds number. Copyright © 2005 John Wiley & Sons, Ltd.     

GENSMAC3D: a numerical method for solving unsteady three‐dimensional free surface flows

A numerical method for solving three‐dimensional free surface flows is presented. The technique is an extension of the GENSMAC code for calculating free surface flows in two dimensions. As in GENSMAC, the full Navier–Stokes equations are solved by a finite difference method; the fluid surface is represented by a piecewise linear surface composed of quadrilaterals and triangles containing marker particles on their vertices; the stress conditions on the free surface are accurately imposed; the conjugate gradient method is employed for solving the discrete Poisson equation arising from a velocity update; and an automatic time step routine is used for calculating the time step at every cycle. A program implementing these features has been interfaced with a solid modelling routine defining the flow domain. A user‐friendly input data file is employed to allow almost any arbitrary three‐dimensional shape to be described. The visualization of the results is performed using computer graphic structures such as phong shade, flat and parallel surfaces. Results demonstrating the applicability of this new technique for solving complex free surface flows, such as cavity filling and jet buckling, are presented. Copyright © 2001 John Wiley & Sons, Ltd.     

Finite‐element simulation of incompressible fluid flow in an elastic vessel

Finite‐element simulation was performed to predict the incompressible Navier–Stokes flow in a domain, partly bounded by an elastic vessel, which is allowed to vary with time. Besides satisfying the physical conservation laws, both surface and the volume conservation laws are satisfied at the discrete level for ensuring the balance between physical and geometrical variables. Several problems which are amenable to analytical solutions were tested for validating the method. The simulated results are observed to agree favourably with analytical solutions. Having verified the applicability of the finite‐element code to problems involving moving grids, we consider an incompressible fluid flow bounded by rigid and elastic vessel walls. Our emphasis was placed on the validation of the formulation developed within the moving‐grid framework. Copyright © 2003 John Wiley & Sons, Ltd.     

Developing a unified FVE‐ALE approach to solve unsteady fluid flow with moving boundaries

In this study, an arbitrary Lagrangian–Eulerian (ALE) approach is incorporated with a mixed finite‐volume–element (FVE) method to establish a novel moving boundary method for simulating unsteady incompressible flow on non‐stationary meshes. The method collects the advantages of both finite‐volume and finite‐element (FE) methods as well as the ALE approach in a unified algorithm. In this regard, the convection terms are treated at the cell faces using a physical‐influence upwinding scheme, while the diffusion terms are treated using bilinear FE shape functions. On the other hand, the performance of ALE approach is improved by using the Laplace method to improve the hybrid grids, involving triangular and quadrilateral elements, either partially or entirely. The use of hybrid FE grids facilitates this achievement. To show the robustness of the unified algorithm, we examine both the first‐ and the second‐order temporal stencils. The accuracy and performance of the extended method are evaluated via simulating the unsteady flow fields around a fixed cylinder, a transversely oscillating cylinder, and in a channel with an indented wall. The numerical results presented demonstrate significant accuracy benefits for the new hybrid method on coarse meshes and where large time steps are taken. Of importance, the current method yields the second‐order temporal accuracy when the second‐order stencil is used to discretize the unsteady terms. Copyright © 2009 John Wiley & Sons, Ltd.     

Numerical method for unsteady viscous hydrodynamical problem with free boundaries

A finite element method for the transient incompressible Navier–Stokes equations with the ability to handle multiple free boundaries is presented. Problems of liquid–liquid type are treated by solving two coupled Navier–Stokes problems for two separate phases. The possibility to solve problems of liquid–gas, liquid–liquid–gas or liquid–liquid–liquid type is demonstrated too. Surface tension effects are included at deformable interfaces. The method is of Lagrangian type with mesh redefinition. A predictor-corrector scheme is used to compute the position of the deformable interface with automatic control of its accuracy and smoothness. The method is provided with an automatic choice of the time integration step and an optional spline filtration of the truncation error at the free surface. In order to show the accuracy of the method, tests and comparisons are presented. Numerical examples include motion of bubbles and multiple drops.     

A consistent splitting scheme for unsteady incompressible viscous flows I. Dirichlet boundary condition and applications

A well‐recognized approach for handling the incompressibility constraint by operating directly on the discretized Navier–Stokes equations is used to obtain the decoupling of the pressure from the velocity field. By following the current developments by Guermond and Shen, the possibilities of obtaining accurate pressure and reducing boundary‐layer effect for the pressure are analysed. The present study mainly reports the numerical solutions of an unsteady Navier–Stokes problem based on the so‐called consistent splitting scheme (J. Comput. Phys. 2003; 192 :262–276). At the same time the Dirichlet boundary value conditions are considered. The accuracy of the method is carefully examined against the exact solution for an unsteady flow physics problem in a simply connected domain. The effectiveness is illustrated viz. several computations of 2D double lid‐driven cavity problems. Copyright © 2005 John Wiley & Sons, Ltd.     

A hybrid Lagrangian–Eulerian particle‐level set method for numerical simulations of two‐fluid turbulent flows

A coupled Lagrangian interface‐tracking and Eulerian level set (LS) method is developed and implemented for numerical simulations of two‐fluid flows. In this method, the interface is identified based on the locations of notional particles and the geometrical information concerning the interface and fluid properties, such as density and viscosity, are obtained from the LS function. The LS function maintains a signed distance function without an auxiliary equation via the particle‐based Lagrangian re‐initialization technique. To assess the new hybrid method, numerical simulations of several ‘standard interface‐moving’ problems and two‐fluid laminar and turbulent flows are conducted. The numerical results are evaluated by monitoring the mass conservation, the turbulence energy spectral density function and the consistency between Eulerian and Lagrangian components. The results of our analysis indicate that the hybrid particle‐level set method can handle interfaces with complex shape change, and can accurately predict the interface values without any significant (unphysical) mass loss or gain, even in a turbulent flow. The results obtained for isotropic turbulence by the new particle‐level set method are validated by comparison with those obtained by the ‘zero Mach number’, variable‐density method. For the cases with small thermal/mass diffusivity, both methods are found to generate similar results. Analysis of the vorticity and energy equations indicates that the destabilization effect of turbulence and the stability effect of surface tension on the interface motion are strongly dependent on the density and viscosity ratios of the fluids. Copyright © 2007 John Wiley & Sons, Ltd.     

Three‐dimensional finite element modelling of coupled free/porous flows: applications to industrial and environmental flows

Conjunctive modelling of free/porous flows provides a powerful and cost‐effective tool for designing industrial filters used in the process industry and also for quantifying surface–subsurface flow interactions, which play a significant role in urban flooding mechanisms resulting from sea‐level rise and climate changes. A number of well‐established schemes are available in the literature for simulation of such regimes; however, three‐dimensional (3D) modelling of such flow systems still presents numerical and practical challenges. This paper presents the development of a fully 3D, transient finite element model for the prediction and quantitative analyses of the hydrodynamic behaviour encountered in industrial filtrations and environmental flows represented by coupled flows. The weak‐variational formulation in this model is based on the use of C⁰ continuous equal‐order Lagrange polynomial functions for velocity and pressure fields represented by 3D hexahedral finite elements. A mixed UVWP finite element scheme based on the standard Galerkin technique satisfying the Ladyzhenskaya–Babuska–Brezzi stability criterion through incorporation of an artificial compressibility term in the continuity equation has been employed for the solution of coupled partial differential equations. We prove that the discretization generates unified stabilization for both the Navier–Stokes and Darcy equations and preserves the geometrical flexibility of the computational grids. A direct node‐linking procedure involving the rearrangement of the global stiffness matrix for the interface elements has been developed by the authors, which is utilized to couple the governing equations in a single model. A variety of numerical tests are conducted, indicating that the model is capable of yielding theoretically expected and accurate results for free, porous and coupled free/porous problems encountered in industrial and environmental engineering problems representing complex filtration (dead‐end and cross‐flow) and interacting surface–subsurface flows. The model is computationally cost‐effective, robust, reliable and easily implementable for practical design of filtration equipments, investigation of land use for water resource availability and assessment of the impacts of climatic variations on environmental catastrophes (i.e. coastal and urban floods). The model developed in this work results from the extension of a multi‐disciplinary project (AEROFIL) primarily sponsored by the European aerospace industries for development of a computer simulation package (Aircraft Cartridge Filter Analysis Modelling Program), which was successfully utilized and deployed for designing hydraulic dead‐end filters used in Airbus A380.Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical simulation of free‐surface flow using the level‐set method with global mass correction

A new numerical method that couples the incompressible Navier–Stokes equations with the global mass correction level‐set method for simulating fluid problems with free surfaces and interfaces is presented in this paper. The finite volume method is used to discretize Navier–Stokes equations with the two‐step projection method on a staggered Cartesian grid. The free‐surface flow problem is solved on a fixed grid in which the free surface is captured by the zero level set. Mass conservation is improved significantly by applying a global mass correction scheme, in a novel combination with third‐order essentially non‐oscillatory schemes and a five stage Runge–Kutta method, to accomplish advection and re‐distancing of the level‐set function. The coupled solver is applied to simulate interface change and flow field in four benchmark test cases: (1) shear flow; (2) dam break; (3) travelling and reflection of solitary wave and (4) solitary wave over a submerged object. The computational results are in excellent agreement with theoretical predictions, experimental data and previous numerical simulations using a RANS‐VOF method. The simulations reveal some interesting free‐surface phenomena such as the free‐surface vortices, air entrapment and wave deformation over a submerged object. Copyright © 2009 John Wiley & Sons, Ltd.     

Global volume conservation in unsteady free surface flows with energy absorbing far‐end boundaries

A wave absorption filter for the far‐end boundary of semi‐infinite large reservoirs is developed for numerical simulation of unsteady free surface flows. Mathematical model is based on finite volume solution of the Navier–Stokes equations and depth‐integrated continuity equation to track the free surface. The Sommerfeld boundary condition is applied at the far‐end of the truncated computational domain. A dissipation zone is formed by applying artificial pressure on water surface to dissipate the kinetic energy of the outgoing waves. The computational scheme is tested to verify the conservation of total fluid volume in the domain for long simulation durations. Combination of the Sommerfeld boundary and dissipation zone can effectively minimize reflections and prevent cumulative changes in total fluid volume in the domain. Solitary wave, nonlinear periodic waves and irregular waves are simulated to illustrate the numerical developments. Earthquake excited surface waves and nonlinear hydrodynamic pressures in a dam–reservoir are computed. Copyright © 2009 John Wiley & Sons, Ltd.     

Two‐level consistent splitting methods based on three corrections for the time‐dependent Navier–Stokes equations

Three kinds of two‐level consistent splitting algorithms for the time‐dependent Navier–Stokes equations are discussed. The basic technique of two‐level type methods for solving the nonlinear problem is first to solve a nonlinear problem in a coarse‐level subspace, then to solve a linear equation in a fine‐level subspace. Hence, the two‐level methods can save a lot of work compared with the one‐level methods. The approaches to linearization are considered based on Stokes, Newton, and Oseen corrections. The stability and convergence demonstrate that the two‐level methods can acquire the optimal accuracy with the proper choice of the coarse and fine mesh scales. Numerical examples show that Stokes correction is the simplest, Newton correction has the best accuracy, while Oseen correction is preferable for the large Reynolds number problems and the long‐time simulations among the three methods. Copyright © 2015 John Wiley & Sons, Ltd.     

A coupling strategy based on anisotropic mesh adaptation for solving two‐fluid flows

We propose a coupling strategy for solving efficiently bifluid flows based on the Stokes equations. Our approach relies on a level set formulation of the interface‐capturing problem, and involves a finite element discretization for the fluid resolution, the method of characteristics for solving the advection of the interface and the anisotropic mesh adaptation of the computational domain in the vicinity of the interface for better accuracy. Copyright © 2010 John Wiley & Sons, Ltd.