A Roe scheme for a compressible six‐equation two‐fluid model

We derive a partially analytical Roe scheme with wave limiters for the compressible six‐equation two‐fluid model. Specifically, we derive the Roe averages for the relevant variables. First, the fluxes are split into convective and pressure parts. Then, independent Roe conditions are stated for these two parts. These conditions are successively reduced while defining acceptable Roe averages. For the convective part, all the averages are analytical. For the pressure part, most of the averages are analytical, whereas the remaining averages are dependent on the thermodynamic equation of state. This gives a large flexibility to the scheme with respect to the choice of equation of state. Furthermore, this model contains nonconservative terms. They are a challenge to handle right, and it is not the object of this paper to discuss this issue. However, the Roe averages presented in this paper are fully independent from how those terms are handled, which makes this framework compatible with any treatment of nonconservative terms. Finally, we point out that the eigenspace of this model may collapse, making the Roe scheme inapplicable. This is called resonance. We propose a fix to handle this particular case. Numerical tests show that the scheme performs well. Copyright © 2012 John Wiley & Sons, Ltd.     

A pressure correction‐volume of fluid method for simulation of two‐phase flows

A pressure correction method coupled with the volume of fluid (VOF) method is developed to simulate two‐phase flows. A volume fraction function is introduced in the VOF method and is governed by an advection equation. A modified monotone upwind scheme for a conservation law (modified MUSCL) is used to solve the solution of the advection equation. To keep the initial sharpness of an interface, a slope modification scheme is introduced. The continuum surface tension (CST) model is used to calculate the surface tension force. Three schemes, central‐upwind, Parker–Youngs, and mixed schemes, are introduced to compute the interface normal vector and the gradient of the volume fraction function. Moreover, a height function technique is applied to compute the local curvature of the interface. Several basic test problems are performed to check the order of accuracy of the present numerical schemes for computing the interface normal vector and the gradient of the volume fraction function. Three physical problems, two‐dimensional broken dam problem, static drop, and spurious currents, and three‐dimensional rising bubble, are performed to demonstrate the efficiency and accuracy of the pressure correction method. Copyright © 2010 John Wiley & Sons, Ltd.     

The multi‐stage centred‐scheme approach applied to a drift‐flux two‐phase flow model

For two‐phase flow models, upwind schemes are most often difficult do derive, and expensive to use. Centred schemes, on the other hand, are simple, but more dissipative. The recently proposed multi‐stage (MUSTA ) method is aimed at coming close to the accuracy of upwind schemes while retaining the simplicity of centred schemes. So far, the MUSTA approach has been shown to work well for the Euler equations of inviscid, compressible single‐phase flow. In this work, we explore the MUSTA scheme for a more complex system of equations: the drift‐flux model, which describes one‐dimensional two‐phase flow where the motions of the phases are strongly coupled. As the number of stages is increased, the results of the MUSTA scheme approach those of the Roe method. The good results of the MUSTA scheme are dependent on the use of a large‐enough local grid. Hence, the main benefit of the MUSTA scheme is its simplicity, rather than CPU ‐time savings. Copyright © 2006 John Wiley & Sons, Ltd.     

Numerical study of wave propagation in compressible two‐phase flow

We propose a new model and a solution method for two‐phase two‐fluid compressible flows. The model involves six equations obtained from conservation principles applied to a one‐dimensional flow of gas and liquid mixture completed by additional closure governing equations. The model is valid for pure fluids as well as for fluid mixtures. The system of partial differential equations with source terms is hyperbolic and has conservative form. Hyperbolicity is obtained using the principles of extended thermodynamics. Features of the model include the existence of real eigenvalues and a complete set of independent eigenvectors. Its numerical solution poses several difficulties. The model possesses a large number of acoustic and convective waves and it is not easy to upwind all of these accurately and simply. In this paper we use relatively modern shock‐capturing methods of a centred‐type such as the total variation diminishing (TVD) slope limiter centre (SLIC) scheme which solve these problems in a simple way and with good accuracy. Several numerical test problems are displayed in order to highlight the efficiency of the study we propose. The scheme provides reliable results, is able to compute strong shock waves and deals with complex equations of state. Copyright © 2006 John Wiley & Sons, Ltd.     

The basic character of five two‐phase flow model equation sets

Three test problems were simulated using five different two‐phase flow model equation sets from the open literature. The test problems chosen were a fluidized bed, a batch settling, and a horizontal jet impingement on a vertical wall. These three problems demonstrate an important cross‐section of physical phenomena, such as fluidized bed voidage oscillations, phase separation, countercurrent flow, and jet formation. The dispersed flow regime is selected for all three problems. The study was performed to assess the basic character of the five‐field equation sets responding to the same initial and boundary conditions and using the same finite difference numerical scheme. The general performance of the five equation sets was found to be similar, even though one of them was ill posed as an initial‐value problem. Broad trends are the same and quantitative differences could be assessed by examining the fine structure of the results. None of the equation sets could be entirely rejected on the basis of producing physically impossible or unacceptable results. Copyright © 2000 John Wiley & Sons, Ltd.     

The numerical solution of the transient two‐phase flow in rigid pipelines

Consideration is given in this paper to the numerical solution of the transient two‐phase flow in rigid pipelines. The governing equations for such flows are two coupled, non‐linear, hyperbolic, partial differential equations with pressure dependent coefficients. The fluid pressure and velocity are considered as two principle dependent variables. The fluid is a homogeneous gas–liquid mixture for which the density is defined by an expression averaging the two‐component densities where a polytropic process of the gaseous phase is admitted. Instead of the void fraction, which varies with the pressure, the gas–fluid mass ratio (or the quality) is assumed to be constant, and is used in the mathematical formulation. The problem has been solved by the method of non‐linear characteristics and the finite difference conservative scheme. To verify their validity, the computed results of the two numerical techniques are compared for different values of the quality, in the case where the liquid compressibility and the pipe wall elasticity are neglected. Copyright © 1999 John Wiley & Sons, Ltd.     

Projection and partitioned solution for two‐phase flow problems

Multiphase flow through porous media is a highly nonlinear process that can be solved numerically with the aid of finite elements (FE) in space and finite differences (FD) in time. For an accurate solution much refined FE grids are generally required with the major computational effort consisting of the resolution to the nonlinearity frequently obtained with the classical Picard linearization approach. The efficiency of the repeated solution to the linear systems within each individual time step represents the key to improve the performance of a multiphase flow simulator. The present paper discusses the performance of the projection solvers (GMRES with restart, TFQMR, and BiCGSTAB) for two global schemes based on a different nodal ordering of the unknowns (ORD1 and ORD2) and a scheme (SPLIT) based on the straightforward inversion of the lumped mass matrix which allows for the preliminary elimination and substitution of the unknown saturations. It is shown that SPLIT is between two and three time faster than ORD1 and ORD2, irrespective of the solver used. Copyright © 2005 John Wiley & Sons, Ltd.     

Evaluation of some high‐order shock capturing schemes for direct numerical simulation of unsteady two‐dimensional free flows

The present study addresses the capability of a large set of shock‐capturing schemes to recover the basic interactions between acoustic, vorticity and entropy in a direct numerical simulation (DNS) framework. The basic dispersive and dissipative errors are first evaluated by considering the advection of a Taylor vortex in a uniform flow. Two transonic cases are also considered. The first one consists of the interaction between a temperature spot and a weak shock. This test emphasizes the capability of the schemes to recover the production of vorticity through the baroclinic process. The second one consists of the interaction of a Taylor vortex with a weak shock, corresponding to the framework of the linear theory of Ribner. The main process in play here is the production of an acoustic wave. The results obtained by using essentially non‐oscillatory (ENO), total variation diminishing (TVD), compact‐TVD and MUSCL schemes are compared with those obtained by means of a sixth‐order accurate Hermitian scheme, considered as reference. The results are as follows; the ENO schemes agree pretty well with the reference scheme. The second‐order accurate Upwind‐TVD scheme exhibits a strong numerical diffusion, while the MUSCL scheme behavior is very sensitive to the value on the parameter β in the limiter function minmod. The compact‐TVD schemes do not yield improvement over the standard TVD schemes. Copyright © 2000 John Wiley & Sons, Ltd.     

A further work on multi‐phase two‐fluid approach for compressible multi‐phase flows

This paper is to continue our previous work Niu (Int. J. Numer. Meth. Fluids 2001; 36 :351–371) on solving a two‐fluid model for compressible liquid–gas flows using the AUSMDV scheme. We first propose a pressure–velocity‐based diffusion term originally derived from AUSMDV scheme Wada and Liou (SIAM J. Sci. Comput. 1997; 18 (3):633—657) to enhance its robustness. The scheme can be applied to gas and liquid fluids universally. We then employ the stratified flow model Chang and Liou (J. Comput. Physics 2007; 225 :240–873) for spatial discretization. By defining the fluids in different regions and introducing inter‐phasic force on cell boundary, the stratified flow model allows the conservation laws to be applied on each phase, and therefore, it is able to capture fluid discontinuities, such as the fluid interfaces and shock waves, accurately. Several benchmark tests are studied, including the Ransom's Faucet problem, 1D air–water shock tube problems, 2D shock‐water column and 2D shock‐bubble interaction problems. The results indicate that the incorporation of the new dissipation into AUSM⁺‐up scheme and the stratified flow model is simple, accurate and robust enough for the compressible multi‐phase flows. Copyright © 2008 John Wiley & Sons, Ltd.     

Application of a thermodynamically compatible two‐phase flow model to the high‐resolution simulations of compressible gas–magma flow

This paper reports on the application and development of a fully hyperbolic and fully conservative two‐phase flow model for the simulation of gas and magma flow within volcanic processes. The model solves a set of mixture conservation equations for the gas and magma two‐phase flow with velocity non‐equilibrium. In this model, the effect of the relative velocity is introduced by a kinetic constitutive equation with other equations for volume and mass fractions of the gas phase. The model is examined numerically by the widely used finite volume Godunov methods of centered‐type. Using the Riemann problem, we numerically simulate wave propagation and the development of shocks and rarefactions in volcanic eruptions. These simulations are of magma fragmentation type where the relative velocity continues to dominate. A series of test cases whose solution contains features relevant to gas–magma mixtures are conducted. In particular, numerical results indicate that the model implementation predicts key features of the relative velocity within volcanic processes without any mathematical or physical simplifications. Simulation results are sharply and accurately provided without any spurious oscillations in all of the flow variables. The numerical methods and results are also compared with other numerical methods available in the literature. It is found that the provided resolutions are more accurate for the considered test cases. Copyright © 2014 John Wiley & Sons, Ltd.     

An optimizing implicit difference scheme based on proper orthogonal decomposition for the two‐dimensional unsaturated soil water flow equation

An optimizing reduced implicit difference scheme (IDS) based on singular value decomposition (SVD) and proper orthogonal decomposition (POD) for the two‐dimensional unsaturated soil water flow equation is presented. An ensemble of snapshots is compiled from the transient solutions derived from the usual IDS for a two‐dimensional unsaturated flow equation. Then, optimal orthogonal bases are reconstructed by implementing SVD and POD techniques for the ensemble of snapshots. Combining POD with a Galerkin projection approach, a new lower dimensional and highly accurate IDS for the two‐dimensional unsaturated flow equation is obtained. Error estimates between the true solution, the usual IDS solution, and the reduced IDS solution based on POD basis are derived. Finally, it is shown by means of a numerical example using the technology of local refined grids that the computational load is greatly diminished by using the reduced IDS. Also, the error between the POD approximate solution and the usual IDS solution is proved to be consistent with the derived theoretical results. Thus, both feasibility and efficiency of the POD method are validated. Copyright © 2011 John Wiley & Sons, Ltd.     

Mathematical modelling of two‐phase non‐Newtonian flow in a helical pipe

Governing equations for a two‐phase 3D helical pipe flow of a non‐Newtonian fluid with large particles are derived in an orthogonal helical coordinate system. The Lagrangian approach is utilized to model solid particle trajectories. The interaction between solid particles and the fluid that carries them is accounted for by a source term in the momentum equation for the fluid. The force‐coupling method (FCM), developed by M.R. Maxey and his group, is adopted; in this method the momentum source term is no longer a Dirac delta function but is spread on a numerical mesh by using a finite‐sized envelop with a spherical Gaussian distribution. The influence of inter‐particle and particle–wall collisions is also taken into account. Copyright © 2005 John Wiley & Sons, Ltd.     

A Riemann solver and upwind methods for a two‐phase flow model in non‐conservative form

We present a theoretical solution for the Riemann problem for the five‐equation two‐phase non‐conservative model of Saurel and Abgrall. This solution is then utilized in the construction of upwind non‐conservative methods to solve the general initial‐boundary value problem for the two‐phase flow model in non‐conservative form. The basic upwind scheme constructed is the non‐conservative analogue of the Godunov first‐order upwind method. Second‐order methods in space and time are then constructed via the MUSCL and ADER approaches. The methods are systematically assessed via a series of test problems with theoretical solutions. Copyright © 2005 John Wiley & Sons, Ltd.     

A surface remeshing technique for a Lagrangian description of 3D two‐fluid flow problems

Classical Lagrangian schemes applied to update the front position between two immiscible incompressible fluids have been long recognized to provide a sharp representation of the interface. However, the main drawback of these approaches is the progressive distortion in the distribution of the markers used to identify the material front. To avoid this problem, a 3D interface remeshing algorithm is proposed in this work. In addition, the remeshed front is enforced to preserve the global volume. These aspects are incorporated in an existing fluid dynamics formulation for the analysis of two‐fluid flows problems. The resulting formulation, called as the 3D‐moving Lagrangian interface remeshing technique, is applied in the numerical analysis of two‐fluid flow problems. Copyright © 2009 John Wiley & Sons, Ltd.     

A gas‐kinetic scheme for the modified Baer–Nunziato model of compressible two‐phase flow

Numerical methods for the Baer–Nunziato model of compressible two‐phase flow have attracted much attention in recent years. In this paper, a two‐phase Bhatnagar–Gross–Krook (BGK) model is constructed in which the non‐conservative terms in the Baer–Nunziato model are considered as the external forces and the collisions both with particles of their phases and other phases are taken into consideration. On the basis of this BGK model, the so‐called modified Baer–Nunziato model is derived and a gas‐kinetic scheme for this modified model is presented. The distribution functions are constructed at the cell interface based on the integral solutions of the BGK equations for both phases. Then, numerical fluxes can be obtained by taking moments of the distribution functions, and non‐conservative terms are explicitly introduced into the construction of numerical fluxes. In this method, not only the iterative processes in the exact Riemann solvers are eliminated but also the collisions with the particles of other phases are taken into account. Copyright © 2012 John Wiley & Sons, Ltd.     

Development of new finite volume schemes on unstructured triangular grid for simulating the gas–liquid two‐phase flow

This paper presents a coupled finite volume inner doubly iterative efficient algorithm for linked equations (IDEAL) with level set method to simulate the incompressible gas–liquid two‐phase flows with moving interfaces on unstructured triangular grid. The finite volume IDEAL method on a collocated grid is employed to solve the incompressible two‐phase Navier–Stokes equations, and the level set method is used to capture the moving interfaces. For the sake of mass conservation, an effective second‐order accurate finite volume scheme is developed to solve the level set equation on triangular grid, which can be implemented much easier than the classical high‐order level set solvers. In this scheme, the value of level set function on the boundary of control volume is approximated using a linear combination of a high‐order Larangian interpolation and a second‐order upwind interpolation. By the rotating slotted disk and stretching and shrinking of a circular fluid element benchmark cases, the mass conservation and accuracy of the new scheme is verified. Then the coupled method is applied to two‐phase flows, including a 2D bubble rising problem and a 2D dam breaking problem. The computational results agree well with those reported in literatures and experimental data. Copyright © 2015 John Wiley & Sons, Ltd.     

A two‐phase adaptive finite element method for solid–fluid coupling in complex geometries

In this paper we present a method to solve the Navier–Stokes equations in complex geometries, such as porous sands, using a finite‐element solver but without the complexity of meshing the porous space. The method is based on treating the solid boundaries as a second fluid and solving a set of equations similar to those used for multi‐fluid flow. When combined with anisotropic mesh adaptivity, it is possible to resolve complex geometries starting with an arbitrary coarse mesh. The approach is validated by comparing simulation results with available data in three test cases. In the first we simulate the flow past a cylinder. The second test case compares the pressure drop in flow through random packs of spheres with the Ergun equation. In the last case simulation results are compared with experimental data on the flow past a simplified vehicle model (Ahmed body) at high Reynolds number using large‐eddy simulation (LES). Results are in good agreement with all three reference models. Copyright © 2010 John Wiley & Sons, Ltd.     

Well‐balanced finite difference weighted essentially non‐oscillatory schemes for the blood flow model

The blood flow model maintains the steady‐state solutions, in which the flux gradients are non‐zero but exactly balanced by the source term. In this paper, we design high order finite difference weighted essentially non‐oscillatory (WENO) schemes to this model with such well‐balanced property and at the same time keeping genuine high order accuracy. Rigorous theoretical analysis as well as extensive numerical results all indicate that the resulting schemes verify high order accuracy, maintain the well‐balanced property, and keep good resolution for smooth and discontinuous solutions. Copyright © 2016 John Wiley & Sons, Ltd.