A 5‐equation,transient, hyperbolic, 1‐dimensional model for slug capturing in pipes

A novel numerical scheme for slug capturing in pipes using a 1‐dimensional transient hyperbolic 5‐equation 2‐fluid model is presented. Previous work has shown that 1‐dimensional 2‐fluid models are able to capture slug flow automatically. In this work, a similar approach is further developed using a new numerical scheme, applied to a hyperbolic 5‐equation 2‐fluid model. Starting from a finite volume discretisation of a 5‐equation 2‐fluid hyperbolic model and adding appropriate closure relations, a second‐order code is implemented and applied to air‐water flows in horizontal pipes, simulating the 2‐phase to 1‐phase flow process. The code is evaluated in some common standard test cases. A slug capturing application is also discussed. We show, in an air/water horizontal pipe, slug initiation, growth, and development. Moreover, a grid refinement analysis is performed showing that the method is grid independent and we show the code capability to take into account eventual surface tension effects, through the instantaneous pressure relaxation process. Finally, a prediction of flow regime transitions is shown and compared with a well‐known theoretical flow pattern map in addition to a preliminary comparison of computed slug characteristics against well‐known empirical correlations.     

Two‐dimensional non‐Newtonian injection molding with a new control volume FEM/volume of fluid method

A new method based on volume of fluid for interface tracking in the simulation of injection molding is presented. The proposed method is comprised of two main stages: accumulation and distribution of the volume fraction. In the first stage the equation for the volume fraction with a noninterfacial flux condition is solved. In the second stage the accumulated volume of fluid that arises as a consequence of the application of the first one is dispersed. This procedure guarantees that the fluid fills the available space without dispersion of the interface. The mathematical model is based on two‐phase transport equations that are numerically integrated through the control volume finite element method. The numerical results for the interface position are successfully verified with analytical results and numerical data available in the literature for one‐dimensional and two‐dimensional domains. The transient position of the advance fronts showed an effective and consistent simulation of an injection molding process. The nondispersive volume of fluid method here proposed is implemented for the simulation of nonisothermal injection molding in two‐dimensional cavities. The obtained results are represented as transient interface positions, isotherms and pressure distributions during the injection molding of low density polyethylene. Copyright © 2012 John Wiley & Sons, Ltd.     

A compressible single‐temperature conservative two‐phase model with phase transitions

A model for multidimensional compressible two‐phase flow with pressure and velocity relaxations based on the theory of thermodynamically compatible system is extended to study liquid–gas flows with cavitation. The model assumes for each phase its own pressure and velocity, while a common temperature is considered. The governing equations form a hyperbolic system in conservative form and are derived through the theory of a thermodynamically compatible system. The phase pressure‐equalizing process and the interfacial friction are taken into account in the balance laws for the volume fractions of one phase and for the relative velocity by adding two relaxation source terms, while the phase transition is introduced into the model considering in the balance equation for the mass of one phase the relaxation of the Gibbs free energies of the two phases. A modification of the central finite‐volume Kurganov–Noelle–Petrova method is adopted in this work to solve the homogeneous hyperbolic part, while the relaxation source terms are treated implicitly. In order to investigate the effect of the mass transfer in the solution, a 1D cavitation tube problem is presented. In addition, two 2D numerical simulations regarding cavitation problem are also studied: a cavitating Richtmyer–Meshkov instability and a laser‐induced cavitation problem. Copyright © 2014 John Wiley & Sons, Ltd.     

An immersed‐boundary finite‐volume method for direct simulation of flows with suspended paramagnetic particles

In the present study, we have proposed an immersed‐boundary finite‐volume method for the direct numerical simulation of flows with inertialess paramagnetic particles suspended in a nonmagnetic fluid under an external magnetic field without the need for any model such as the dipole–dipole interaction. In the proposed method, the magnetic field (or force) is described by the numerical solution of the Maxwell equation without current, where the smoothed representation technique is employed to tackle the discontinuity of magnetic permeability across the particle–fluid interface. The flow field, on the other hand, is described by the solution of the continuity and momentum equations, where the discrete‐forcing‐based immersed‐boundary method is employed to satisfy the no‐slip condition at the interface. To validate the method, we performed numerical simulations on the two‐dimensional motion of two and three paramagnetic particles in a nonmagnetic fluid subjected to an external uniform magnetic field and then compared the results with the existing finite‐element and semi‐analytical solutions. Comparison shows that the proposed method is robust in the direct simulation of such magnetic particulate flows. This method can be extended to more general flows without difficulty: three‐dimensional particulate flows, flows with a great number of particles, or flows under an arbitrary external magnetic field. Copyright © 2010 John Wiley & Sons, Ltd.     

Implementation of phase change thermodynamic probability for unsteady simulation of cavitating flows

The aim of this work is to investigate the non‐equilibrium effects of phase change in cavitating flows. For this purpose, the concept of phase change thermodynamic probability is used along with homogeneous model to simulate two‐phase cavitating flows. For simulation of unsteady behaviors of cavitation, which have practical applications, unsteady PISO algorithm based on the non‐conservative approach is utilized. For multi‐phase simulation, single‐fluid Navier–Stokes equations, along with the volume fraction transport equation, are employed. In this paper, phase change thermodynamics probabilities and cavitation model is briefly summarized. Thus, derivation of the cavitation model, starting from the basic thermodynamic equations to the mass and momentum conservation equations at a liquid–vapor two‐phase flow, is presented to explain the numerical model. Unsteady simulations of cavitation around a flat plate normal to flow direction are presented to clarify the accuracy of the model. The accuracy of the numerical results is good, and it is possible to apply this method to more complex geometries. Copyright © 2010 John Wiley & Sons, Ltd.     

Generalized Fourier analyses of the advection–diffusion equation—Part I: one‐dimensional domains

This paper presents a detailed multi‐methods comparison of the spatial errors associated with finite difference, finite element and finite volume semi‐discretizations of the scalar advection–diffusion equation. The errors are reported in terms of non‐dimensional phase and group speed, discrete diffusivity, artificial diffusivity, and grid‐induced anisotropy. It is demonstrated that Fourier analysis provides an automatic process for separating the discrete advective operator into its symmetric and skew‐symmetric components and characterizing the spectral behaviour of each operator. For each of the numerical methods considered, asymptotic truncation error and resolution estimates are presented for the limiting cases of pure advection and pure diffusion. It is demonstrated that streamline upwind Petrov–Galerkin and its control‐volume finite element analogue, the streamline upwind control‐volume method, produce both an artificial diffusivity and a concomitant phase speed adjustment in addition to the usual semi‐discrete artifacts observed in the phase speed, group speed and diffusivity. The Galerkin finite element method and its streamline upwind derivatives are shown to exhibit super‐convergent behaviour in terms of phase and group speed when a consistent mass matrix is used in the formulation. In contrast, the CVFEM method and its streamline upwind derivatives yield strictly second‐order behaviour. In Part II of this paper, we consider two‐dimensional semi‐discretizations of the advection–diffusion equation and also assess the affects of grid‐induced anisotropy observed in the non‐dimensional phase speed, and the discrete and artificial diffusivities. Although this work can only be considered a first step in a comprehensive multi‐methods analysis and comparison, it serves to identify some of the relative strengths and weaknesses of multiple numerical methods in a common analysis framework. Published in 2004 by John Wiley & Sons, Ltd.     

流体饱和两相多孔介质拟静态问题的有限元解法

给出基于混合物理论的流体饱和两相多孔介质模型,该模型由一可变形固体 一流体相组成。采用Ｇａｌｅｒｋｉｎ加权残值法导出求解拟静态问题的有限元公式,并编制了二维有限元程序。用程序分析了一维和二维问题,得到合理的结果。     

Simulation of sharp gas–liquid interface using VOF method and adaptive grid local refinement around the interface

The volume of fluid (VOF) method is used to perform two‐phase simulations (gas–liquid). The governing Navier–Stokes conservation equations of the flow field are numerically solved on two‐dimensional axisymmetric or three‐dimensional unstructured grids, using Cartesian velocity components, following the finite volume approximation and a pressure correction method. A new method of adaptive grid local refinement is developed in order to enhance the accuracy of the predictions, to capture the sharp gas–liquid interface and to speed up the calculations. Results are compared with experimental measurements in order to assess the efficiency of the method. Copyright © 2004 John Wiley & Sons, Ltd.     

Front‐tracking by the level‐set and the volume penalization methods in a two‐phase microfluidic network

Two‐phase immiscible fluids in a two‐dimensional micro‐channels network are considered. The incompressible Stokes equations are used to describe the Newtonian fluid flow, while the Oldroyd‐B rheological model is used to capture the viscoelastic behavior. In order to perform numerical simulations in a complex geometry like a micro‐channels network, the volume penalization method is implemented. To follow the interface between the two fluids, the level‐set method is used, and the dynamics of the contact line is modeled by Cox law. Numerical results show the ability of the method to simulate two‐phase flows and to follow properly the contact line between the two immiscible fluids. Finally, simulations with realistic parameters are performed to show the difference when a Newtonian fluid is pushed by a viscoelastic fluid instead of a Newtonian one. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical study of breaking waves by a two‐phase flow model

A two‐phase flow model, which solves the flow in the air and water simultaneously, is presented for modelling breaking waves in deep and shallow water, including wave pre‐breaking, overturning and post‐breaking processes. The model is based on the Reynolds‐averaged Navier–Stokes equations with the k ?ε turbulence model. The governing equations are solved by the finite volume method in a Cartesian staggered grid and the partial cell treatment is implemented to deal with complex geometries. The SIMPLE algorithm is utilised for the pressure‐velocity coupling and the air‐water interface is modelled by the interface capturing method via a high resolution volume of fluid scheme. The numerical model is validated by simulating overturning waves on a sloping beach and over a reef, and deep‐water breaking waves in a periodic domain, in which good agreement between numerical results and available experimental measurements for the water surface profiles during wave overturning is obtained. The overturning jet, air entrainment and splash‐up during wave breaking have been captured by the two‐phase flow model, which demonstrates the capability of the model to simulate free surface flow and wave breaking problems.Copyright © 2011 John Wiley & Sons, Ltd.     

Finite‐volume‐type VOF method on dynamically adaptive quadtree grids

This paper describes a finite‐volume volume‐of‐fluid (VOF) method for simulating viscous free surface flows on dynamically adaptive quadtree grids. The scheme is computationally efficient in that it provides relatively fine grid resolution at the gas–liquid interface and coarse grid density in regions where flow variable gradients are small. Special interpolations are used to ensure volume flux conservation where differently sized neighbour cells occur. The numerical model is validated for advection of dyed fluid in unidirectional and rotating flows, and for two‐dimensional viscous sloshing in a rectangular tank. Copyright © 2004 John Wiley & Sons, Ltd.     

A high‐resolution method for compressible two‐fluid flows and simulation of three‐dimensional shock–bubble interactions

A high‐resolution method is developed to capture the material interfaces of compressible two‐fluid flows in multiple dimensions. A fluid mixture model system with single velocity and pressure is used, and viscous effect can also be taken into account. A consistent thermodynamic law based on the assumption of pressure equilibrium is employed to describe the thermodynamic behaviors of the pure fluids and mixture of two components. The splitting and unsplit Eulerian formulations of piecewise parabolic method are extended to numerically integrate the hyperbolic part of the model system, whereas the system of diffusion equations is solved using an explicit, central difference scheme. The block‐structured adaptive mesh refinement (AMR) capability is built in the hydrodynamic code to locally improve grid resolution. The resulting method is verified to be at least second‐order accurate in space. Numerical results show that the discontinuities, particularly contact discontinuities, can be resolved sharply. The use of AMR allows flow features at disparate scales to be resolved sufficiently. In addition, three‐dimensional shock–bubble interactions are simulated to investigate effects of Mach number on bubble evolution. The flow structures including those peculiar to three‐dimensional bubble are resolved correctly, and some physical phenomena with increasing Mach number are reported. Copyright © 2012 John Wiley & Sons, Ltd.     

The intersection marker method for 3D interface tracking of deformable surfaces in finite volumes

Currently, the majority of computational fluid dynamics (CFD) codes use the finite volume method to spatially discretise the computational domain, sometimes as an array of cubic control volumes. The Finite volume method works well with single‐phase flow simulations, but two‐phase flow simulations are more challenging because of the need to track the surface interface traversing and deforming within the 3D grid. Surface area and volume fraction details of each interface cell must be accurately accounted for, in order to calculate for the momentum exchange and rates of heat and mass transfer across the interface. To attain a higher accuracy in two‐phase flow CFD calculations, the intersection marker (ISM) method is developed. The ISM method is a hybrid Lagrangian–Eulerian front‐tracking algorithm that can model an arbitrary 3D surface within an array of cubic control volumes. The ISM method has a cell‐by‐cell remeshing capability that is volume conservative and is suitable for the tracking of complex interface deformation in transient two‐phase CFD simulations. Copyright © 2015 John Wiley & Sons, Ltd.     

Extension of a high‐resolution scheme to 1D liquid–gas flow

This paper presents the extension of a high‐resolution conservative scheme to the one‐dimensional one‐pressure six‐equation two‐fluid flow model. Only mixtures of water and air have been considered in this study, both fluids have been characterized using simple equations of state, namely stiffened gas for the liquid phase and perfect gas for the gas phase. The resulting scheme is explicit and first‐order accurate in space and time. A second‐order version of the scheme has also been derived using the MUSCL strategy and slope limiters. Some numerical results show the good capabilities of this type of schemes in the solution of discontinuities in two‐fluid flow problems, all of them are based on water/air numerical benchmarks widely used in the two‐phase flow literature. Copyright © 2005 John Wiley & Sons, Ltd.     

A Lagrangian–Eulerian model of particle dispersion in a turbulent plane mixing layer

A Lagrangian–Eulerian model for the dispersion of solid particles in a two‐dimensional, incompressible, turbulent flow is reported and validated. Prediction of the continuous phase is done by solving an Eulerian model using a control‐volume finite element method (CVFEM). A Lagrangian model is also applied, using a Runge–Kutta method to obtain the particle trajectories. The effect of fluid turbulence upon particle dispersion is taken into consideration through a simple stochastic approach. Validation tests are performed by comparing predictions for both phases in a particle‐laden, plane mixing layer airflow with corresponding measurements formerly reported by other authors. Even though some limitations are detected in the calculation of particle dispersion, on the whole the validation results are rather successful. Copyright © 2002 John Wiley & Sons, Ltd.     

Coupling of finite volume method and thermal lattice Boltzmann method and its application to natural convection

On the basis of the existing density distribution function reconstruction operator, the temperature distribution operator was derived to calculate heat transfer by coupling the lattice Boltzmann method (LBM) with the finite volume method. The present coupling model was validated by two‐dimensional natural convection flows with and without an isolated internal vertical plate. The results from the coupling model agree well with those from the pure‐finite volume method, pure‐LBM and references, and all the physical quantities cross the coupled interface smoothly. On the basis of residual history curves, it is likely that the convergence property and the numerical stability of the present model are better than those of the pure‐LBM at fine grid numbers and high Rayleigh numbers. Copyright © 2011 John Wiley & Sons, Ltd.     

Numerical simulation of the unsteady behaviour of cavitating flows

A 2D numerical model is proposed to simulate unsteady cavitating flows. The Reynolds‐averaged Navier–Stokes equations are solved for the mixture of liquid and vapour, which is considered as a single fluid with variable density. The vapourization and condensation processes are controlled by a barotropic state law that relates the fluid density to the pressure variations. The numerical resolution is a pressure‐correction method derived from the SIMPLE algorithm, with a finite volume discretization. The standard scheme is slightly modified to take into account the cavitation phenomenon. That numerical model is used to calculate unsteady cavitating flows in two Venturi type sections. The choice of the turbulence model is discussed, and the standard RNG k–εmodel is found to lead to non‐physical stable cavities. A modified k–εmodel is proposed to improve the simulation. The influence of numerical and physical parameters is presented, and the numerical results are compared to previous experimental observations and measurements. The proposed model seems to describe the unsteady cavitation behaviour in 2D geometries well. Copyright © 2003 John Wiley & Sons, Ltd.     

A meshless method for numerical simulation of depth‐averaged turbulence flows using a k‐ϵ model

A three‐dimensional numerical model is developed to analyze free surface flows and water impact problems. The flow of an incompressible viscous fluid is solved using the unsteady Navier–Stokes equations. Pseudo‐time derivatives are introduced into the equations to improve computational efficiency. The interface between the two phases is tracked using a volume‐of‐fluid interface tracking algorithm developed in a generalized curvilinear coordinate system. The accuracy of the volume‐of‐fluid method is first evaluated by the multiple numerical benchmark tests, including two‐dimensional and three‐dimensional deformation cases on curvilinear grids. The performance and capability of the numerical model for water impact problems are demonstrated by simulations of water entries of the free‐falling hemisphere and cone, based on comparisons of water impact loadings, velocities, and penetrations of the body with experimental data. For further validation, computations of the dam‐break flows are presented, based on an analysis of the wave front propagation, water level, and the dynamic pressure impact of the waves on the downstream walls, on a specific container, and on a tall structure. Extensive comparisons between the obtained solutions, the experimental data, and the results of other numerical simulations in the literature are presented and show a good agreement. Copyright © 2015 John Wiley & Sons, Ltd.     

Shallow flow simulation on dynamically adaptive cut cell quadtree grids

A computationally efficient, high‐resolution numerical model of shallow flow hydrodynamics is described, based on dynamically adaptive quadtree grids. The numerical model solves the two‐dimensional non‐linear shallow water equations by means of an explicit second‐order MUSCL‐Hancock Godunov‐type finite volume scheme. Interface fluxes are evaluated using an HLLC approximate Riemann solver. Cartesian cut cells are used to improve the fit to curved boundaries. A ghost‐cell immersed boundary method is used to update flow information in the smallest cut cells and overcome the time step restriction that would otherwise apply. The numerical model is validated through simulations of reflection of a surge wave at a wall, a low Froude number potential flow past a circular cylinder, and the shock‐like interaction between a bore and a circular cylinder. The computational efficiency is shown to be greatly improved compared with solutions on a uniform structured grid implemented with cut cells. Copyright © 2006 John Wiley & Sons, Ltd.