Multi‐material incompressible flow simulation using the moment‐of‐fluid method

This paper compares the numerical performance of the moment‐of‐fluid (MOF) interface reconstruction technique with Youngs, LVIRA, power diagram (PD), and Swartz interface reconstruction techniques in the context of a volume‐of‐fluid (VOF) based finite element projection method for the numerical simulation of variable‐density incompressible viscous flows. In pure advection tests with multiple materials MOF shows dramatic improvements in accuracy compared with the other methods. In incompressible flows where density differences determine the flow evolution, all the methods perform similarly for two material flows on structured grids. On unstructured grids, the second‐order MOF, LVIRA, and Swartz methods perform similarly and show improvement over the first‐order Youngs' and PD methods. For flow simulations with more than two materials, MOF shows increased accuracy in interface positions on coarse meshes. In most cases, the convergence and accuracy of the computed flow solution was not strongly affected by interface reconstruction method. Published in 2009 by John Wiley & Sons, Ltd.     

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.     

Numerical errors of the volume‐of‐fluid interface tracking algorithm

One of the important limitations of the interface tracking algorithms is that they can be used only as long as the local computational grid density allows surface tracking. In a dispersed flow, where the dimensions of the particular fluid parts are comparable or smaller than the grid spacing, several numerical and reconstruction errors become considerable. In this paper the analysis of the interface tracking errors is performed for the volume‐of‐fluid method with the least squares volume of fluid interface reconstruction algorithm. A few simple two‐fluid benchmarks are proposed for the investigation of the interface tracking grid dependence. The expression based on the gradient of the volume fraction variable is introduced for the estimation of the reconstruction correctness and can be used for the activation of an adaptive mesh refinement algorithm. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical approach for generic three-phase flow based on cut-cell and ghost fluid methods

In this paper, we introduce numerical methods that can simulate complex multiphase flows. The finite volume method, applying Cartesian cut-cell is used in the computational domain, containing fluid and solid, to conserve mass and momentum. With this method, flows in and around any geometry can be simulated without complex and time consuming meshing. For the fluid region, which involves liquid and gas, the ghost fluid method is employed to handle the stiffness of the interface discontinuity problem. The interaction between each phase is treated simply by wall function models or jump conditions of pressure, velocity and shear stress at the interface. The sharp interface method “coupled level set (LS) and volume of fluid (VOF)” is used to represent the interface between the two fluid phases. This approach will combine some advantages of both interface tracking/capturing methods, such as the excellent mass conservation from the VOF method and good accuracy of interface normal computation from the LS function. The first coupled LS and VOF will be generated to reconstruct the interface between solid and the other materials. The second will represent the interface between liquid and gas.     

A stencil-like volume of fluid （VOF） method for tracking free interface

A stencil-like volume of fluid （VOF） method is proposed for tracking free interface. A stencil on a grid cell is worked out according to the normal direction of the interface, in which only three interface positions are possible in 2D cases, and the interface can be reconstructed by only requiring the known local volume fraction information. On the other hand, the fluid-occupying-length is defined on each side of the stencil, through which a unified fluid-occupying volume model and a unified algorithm can be obtained to solve the interface advection equation. The method is suitable for the arbitrary geometry of the grid cell, and is extendible to 3D cases. Typical numerical examples show that the current method can give ＂sharp＂ results for tracking free interface.     

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.     

An efficient method for capturing free boundaries in multi‐fluid simulations

An easy‐to‐use front capturing method is devised by directly solving the transport equation for a volume of fluid (VOF) function. The key to this method is a semi‐Lagrangian conservative scheme, namely CIP_CSL3, recently proposed by the author. In the CIP_CSL3 scheme, the first‐order derivative of the interpolation polynomial at each cell centre is used to control the shape of the reconstructed profile. We show in the present paper that the first‐order derivative, which plays a crucial role in reconstructing the interpolation profile, can also be used to eliminate numerical diffusion. The resulting algorithm can be directly used to compute the VOF‐like function and retain the compact thickness of the moving interface in multi‐fluid simulations. No surface reconstruction based on the value of VOF function is required in the method, which makes it quite economical and easy to use. The presented method has been tested with various interfacial flows including pure rotation, vortex shearing, multi‐vortex deformation and the moving boundaries in real fluid as well. The method gives promising results to all computed problems. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical prediction of turbulent thermocapillary convection in superposed fluid layers with a free interface

The present paper introduces a new numerical method for predicting the characteristics of thermocapillary turbulent convection in a differentially-heated rectangular cavity with two superposed and immiscible fluid layers. The unsteady Reynolds form of the Navier–Stokes equations and energy equation are solved by using the control volume approach on a staggered grid system using SIMPLE algorithm. The turbulence quantities are predicted by applying the standard k–ε turbulence model. The level set formulation is applied for predicting the topological changes of the interface separating the two fluid layers and to provide an accurate and robust modeling of the interfacial normal and tangential stresses. The computational results obtained showed good agreement when compared with the previous experimental, numerical and analytical benchmark data for different validation cases in both laminar and turbulent regimes. The present numerical method is then applied to predict the velocity and temperature distribution in two immiscible liquid layers with undeformable interface for a wide range of Marangoni numbers. The laminar-turbulent transition is demonstrated by obtaining the turbulence features at high interfacial temperature gradient which is characterized by high Marangoni number. The effect of increasing Marangoni number on the interface dynamics in turbulent regime is also investigated.     

Numerical prediction of turbulent thermocapillary convection in superposed fluid layers with a free interface

The present paper introduces a new numerical method for predicting the characteristics of thermocapillary turbulent convection in a differentially-heated rectangular cavity with two superposed and immiscible fluid layers. The unsteady Reynolds form of the Navier–Stokes equations and energy equation are solved by using the control volume approach on a staggered grid system using SIMPLE algorithm. The turbulence quantities are predicted by applying the standard k–ε turbulence model. The level set formulation is applied for predicting the topological changes of the interface separating the two fluid layers and to provide an accurate and robust modeling of the interfacial normal and tangential stresses. The computational results obtained showed good agreement when compared with the previous experimental, numerical and analytical benchmark data for different validation cases in both laminar and turbulent regimes. The present numerical method is then applied to predict the velocity and temperature distribution in two immiscible liquid layers with undeformable interface for a wide range of Marangoni numbers. The laminar-turbulent transition is demonstrated by obtaining the turbulence features at high interfacial temperature gradient which is characterized by high Marangoni number. The effect of increasing Marangoni number on the interface dynamics in turbulent regime is also investigated.     

Explicit/implicit fluid/structure staggered procedures with a structural predictor and fluid subcycling for 2D inviscid aeroelastic simulations

Field time integrators with second-order-accurate numerical schemes for both the fluid and the structure are considered for unsteady Euler aeroelastic computations. We show that if these schemes are simply coupled and used straightforwardly with subcycling, then accuracy and stability properties may be lost. We present new coupling staggered procedures where momentum conservation is enforced at the interface. This is done by using a structural predictor. Continuity of structural and fluid grid displacements is not satisfied at the fluid/structure interface. However, we show on a two-degree-of-freedom aerofoil that this new type of method has many advantages, e.g. accuracy of conservation at the interface and extended stability. The supersonic flutter of a flat panel is simulated in order to numerically prove that the algorithm gives accurate results with arbitrary subcycling for the fluid in the satisfying limit of 30 time steps per period of coupled oscillation. © 1997 John Wiley & Sons, Ltd.     

A numerical method for flows in porous and homogenous fluid domains coupled at the interface by stress jump

A numerical method was developed for flows involving an interface between a homogenous fluid and a porous medium. The numerical method is based on the finite volume method with body‐fitted and multi‐block grids. A generalized model, which includes Brinkman term, Forcheimmer term and non‐linear convective term, was used to govern the flow in the porous medium region. At its interface, a shear stress jump that includes the inertial effect was imposed, together with a continuity of normal stress. Furthermore, the effect of the jump condition on the diffusive flux was considered, additional to that on the convective part which has been usually considered. Numerical results of three flow configurations are presented. The method is suitable for coupled problems with regions of homogeneous fluid and porous medium, which have complex geometries. Copyright © 2006 John Wiley & Sons, Ltd.     

Development of a compressive surface capturing formulation for modelling free‐surface flow by using the volume‐of‐fluid approach

With the aim of accurately modelling free‐surface flow of two immiscible fluids, this study presents the development of a new volume‐of‐fluid free‐surface capturing formulation. By building on existing volume‐of‐fluid approaches, the new formulation combines a blended higher resolution scheme with the addition of an artificial compressive term to the volume‐of‐fluid equation. This reduces the numerical smearing of the interface associated with explicit higher resolution schemes while limiting the contribution of the artificial compressive term to ensure the integrity of the interface shape is maintained. Furthermore, the computational efficiency of the the higher resolution scheme is improved through the reformulation of the normalised variable approach and the implementation of a new higher resolution blending function. The volume‐of‐fluid equation is discretised via an unstructured vertex‐centred finite volume method and solved via a Jacobian‐type dual time‐stepping approach. Copyright © 2012 John Wiley & Sons, Ltd.     

On the composite Riemann problem for multi‐material fluid flows

We develop an Eulerian fixed grid numerical method for calculating multi‐material fluid flows. This approach relates to the class of interface capturing methods. The fluid is treated as a heterogeneous mixture of constituent materials, and the material interface is implicitly captured by a region of mixed cells that have arisen owing to numerical diffusion. To suppress this numerical diffusion, we propose a composite Riemann problem (CRP), which describes the decay of an initial discontinuity in the presence of a contact point between two different fluids, which is located off the initial discontinuity point. The solution to the CRP serves to calculate multi‐material no mixed numerical flux without introducing any material diffusion. We discuss the CRP solution and its implementation in the multi‐material fluid Godunov method. Numerical results show that a simple framework of the CRP greatly improves capturing material interfaces in the Godunov method and reproduces many of the advantages of more complicated interface tracking multi‐material treatments. Copyright © 2014 John Wiley & Sons, Ltd.     

A stable unstructured finite volume method for parallel large-scale viscoelastic fluid flow calculations

A new stable unstructured finite volume method is presented for parallel large-scale simulation of viscoelastic fluid flows. The numerical method is based on the side-centered finite volume method where the velocity vector components are defined at the mid-point of each cell face, while the pressure term and the extra stress tensor are defined at element centroids. The present arrangement of the primitive variables leads to a stable numerical scheme and it does not require any ad-hoc modifications in order to enhance the pressure–velocity–stress coupling. The log-conformation representation proposed in [R. Fattal, R. Kupferman, Constitutive laws for the matrix–logarithm of the conformation tensor, J. Non-Newtonian Fluid Mech. 123 (2004) 281–285] has been implemented in order improve the limiting Weissenberg numbers in the proposed finite volume method. The time stepping algorithm used decouples the calculation of the polymeric stress by solution of a hyperbolic constitutive equation from the evolution of the velocity and pressure fields by solution of a generalized Stokes problem. The resulting algebraic linear systems are solved using the FGMRES(m) Krylov iterative method with the restricted additive Schwarz preconditioner for the extra stress tensor and the geometric non-nested multilevel preconditioner for the Stokes system. The implementation of the preconditioned iterative solvers is based on the PETSc library for improving the efficiency of the parallel code. The present numerical algorithm is validated for the Kovasznay flow, the flow of an Oldroyd-B fluid past a confined circular cylinder in a channel and the three-dimensional flow of an Oldroyd-B fluid around a rigid sphere falling in a cylindrical tube. Parallel large-scale calculations are presented up to 523,094 quadrilateral elements in two-dimension and 1,190,376 hexahedral elements in three-dimension.     

A finite volume method for fluid flow in polar cylindrical grids

In the numerical simulation of fluid flows using a polar cylindrical grid, grid lines meet at a single point on the axis of the polar cylindrical grid system; this makes the grids around the axis degenerate from being general quadrilaterals into triangles. Therefore, a special treatment must be performed when the axis has to be included in the computational domain in order to solve a non-axisymmetrical fluid flow problem. In this paper a new numerical method has been developed to deal with the difficulty of the axis when the control volume technique is used with a non-staggered grid arrangement. Two illustrative examples of the proposed method are presented for simulating the fluid flows on the axis and all the numerical results obtained for the two examples are shown to be in good agreement with the available analytical solutions. © 1998 John Wiley & Sons, Ltd.     

Convected level set method for the numerical simulation of fluid buckling

‘Fluid buckling’ is a phenomenon characterized mainly by the existence of fluid toroidal oscillations during flow. It appears when a high viscosity fluid flows vertically against a flat surface and may occur in industrial applications, as in injection molding of a propergol in complex‐shaped cavities. These coiling or folding oscillations appear during the mold filling stage, leading to air entrapment. To understand and to model this free surface flow problem, a convected level set method is proposed. First, a sinus filter is applied to the distance function to get a smooth truncation far from the interface. Second, the reinitialization is embedded in the transport equation model, avoiding it as a separate step during calculation. In order to validate the method, numerical results are presented on classical interface capturing benchmarks. Finally, results are shown on two‐dimensional and three‐dimensional viscous jet buckling problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Volume preserving immersed boundary methods for two‐phase fluid flows

In this article, we propose a simple area‐preserving correction scheme for two‐phase immiscible incompressible flows with an immersed boundary method (IBM). The IBM was originally developed to model blood flow in the heart and has been widely applied to biofluid dynamics problems with complex geometries and immersed elastic membranes. The main idea of the IBM is to use a regular Eulerian computational grid for the fluid mechanics along with a Lagrangian representation of the immersed boundary. Using the discrete Dirac delta function and the indicator function, we can include the surface tension force, variable viscosity and mass density, and gravitational force effects. The principal advantage of the IBM for two‐phase fluid flows is its inherent accuracy due in part to its ability to use a large number of interfacial marker points on the interface. However, because the interface between two fluids is moved in a discrete manner, this can result in a lack of volume conservation. The idea of an area preserving correction scheme is to correct the interface location normally to the interface so that the area remains constant. Various numerical experiments are presented to illustrate the efficiency and accuracy of the proposed conservative IBM for two‐phase fluid flows. Copyright © 2011 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 numerical continuous model for the hydrodynamics of fluid particle systems

In order to understand the hydrodynamic interactions that can appear in a fluid particle motion, an original method based on the equations governing the motion of two immiscible fluids has been developed. These momentum equations are solved for both the fluid and solid phases. The solid phase is assumed to be a fluid phase with physical properties, such as its behaviour can be assimilated to that of pseudo‐rigid particles. The only unknowns are the velocity and the pressure defined in both phases. The unsteady two‐dimensional momentum equations are solved by using a staggered finite volume formulation and a projection method. The transport of each particle is solved by using a second‐order explicit scheme. The physical model and the numerical method are presented, and the method is validated through experimental measurements and numerical results concerning the flow around a circular cylinder. Good agreement is observed in most cases. The method is then applied to study the trajectory of one settling particle initially off‐centred between two parallel walls and the corresponding wake effects. Different particle trajectories related to particulate Reynolds numbers are presented and commented. A two‐body interaction problem is investigated too. This method allows the simulation of the transport of particles in a dilute suspension in reasonable time. One of the important features of this method is the computational cost that scales linearly with the number of particles. Copyright © 1999 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.