Hybrid finite element/volume method for 3D incompressible flows with heat transfer

An implicit hybrid finite element (FE)/volume solver has been extended to incompressible flows coupled with the energy equation. The solver is based on the segregated pressure correction or projection method on staggered unstructured hybrid meshes. An intermediate velocity field is first obtained by solving the momentum equations with the matrix-free implicit cell-centred finite volume (FV) method. The pressure Poisson equation is solved by the node-based Galerkin FE method for an auxiliary variable. The auxiliary variable is used to update the velocity field and the pressure field. The pressure field is carefully updated by taking into account the velocity divergence field. Our current staggered-mesh scheme is distinct from other conventional ones in that we store the velocity components at cell centres and the auxiliary variable at vertices. The Generalized Minimal Residual (GMRES) matrix-free strategy is adapted to solve the governing equations in both FE and FV methods. The presented 2D and 3D numerical examples show the robustness and accuracy of the numerical method.     

An implementation of the Spalart–Allmaras DES model in an implicit unstructured hybrid finite volume/element solver for incompressible turbulent flow

In this paper, we describe an implicit hybrid finite volume (FV)/element (FE) incompressible Navier–Stokes solver for turbulent flows based on the Spalart–Allmaras detached eddy simulation (SA‐DES). The hybrid FV/FE solver is based on the segregated pressure correction or projection method. The intermediate velocity field is first obtained by solving the original momentum equations with the matrix‐free implicit cell‐centered FV method. The pressure Poisson equation is solved by the node‐based Galerkin FE method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centers and the auxiliary variable at vertices, making the current solver a staggered‐mesh scheme. The SA‐DES turbulence equation is solved after the velocity and the pressure fields have been updated at the end of each time step. The same matrix‐free FV method as the one used for momentum equations is used to solve the turbulence equation. The turbulence equation provides the eddy viscosity, which is added to the molecular viscosity when solving the momentum equation. In our implementation, we focus on the accuracy, efficiency and robustness of the SA‐DES model in a hybrid flow solver. This paper will address important implementation issues for high‐Reynolds number flows where highly stretched elements are typically used. In addition, some aspects of implementing the SA‐DES model will be described to ensure the robustness of the turbulence model. Several numerical examples including a turbulent flow past a flat plate and a high‐Reynolds number flow around a high angle‐of‐attack NACA0015 airfoil will be presented to demonstrate the accuracy and efficiency of our current implementation. Copyright © 2008 John Wiley & Sons, Ltd.     

Development of a hybrid finite volume/element solver for incompressible flows

In this paper, we report our development of an implicit hybrid flow solver for the incompressible Navier–Stokes equations. The methodology is based on the pressure correction or projection method. A fractional step approach is used to obtain an intermediate velocity field by solving the original momentum equations with the matrix‐free implicit cell‐centred finite volume method. The Poisson equation derived from the fractional step approach is solved by the node‐based Galerkin finite element method for an auxiliary variable. The auxiliary variable is closely related to the real pressure and is used to update the velocity field and the pressure field. We store the velocity components at cell centres and the auxiliary variable at cell vertices, making the current solver a staggered‐mesh scheme. Numerical examples demonstrate the performance of the resulting hybrid scheme, such as the correct temporal convergence rates for both velocity and pressure, absence of unphysical pressure boundary layer, good convergence in steady‐state simulations and capability in predicting accurate drag, lift and Strouhal number in the flow around a circular cylinder. Copyright © 2007 John Wiley & Sons, Ltd.     

Efficient solving method for unsteady incompressible interfacial flow problems

Unsteady interfacial problems, considered in an Eulerian form, are studied. The phenomena are modeled using the incompressible viscous Navier–Stokes equations to get the velocity field and an advection equation to predict interface evolutions. The momentum equation is solved by means of an implicit hybrid augmented Lagrangian–Projection method, whereas an explicit characteristic method coupled with a TVD SUPERBEE scheme is applied to the advection equation. The velocity components and the pressure are discretized on staggered grids with finite volumes. Emphasis is on the accuracy and robustness of the techniques described before. A precise explanation on the validation phase will be given, which uses such tests as the advection of a step function or Zalesak's problem to improve the calculation of the interface. The global approach is used on a physically hard interfacial test with strong disparities between viscosities and densities. Copyright © 1999 John Wiley & Sons, Ltd.     

A finite volume approach for unsteady viscoelastic fluid flows

A finite volume, time‐marching for solving time‐dependent viscoelastic flow in two space dimensions for Oldroyd‐B and Phan Thien–Tanner fluids, is presented. A non‐uniform staggered grid system is used. The conservation and constitutive equations are solved using the finite volume method with an upwind scheme for the viscoelastic stresses and an hybrid scheme for the velocities. To calculate the pressure field, the semi‐implicit method for the pressure linked equation revised method is used. The discretized equations are solved sequentially, using the tridiagonal matrix algorithm solver with under‐relaxation. In both, the full approximation storage multigrid algorithm is used to speed up the convergence rate. Simulations of viscoelastic flows in four‐to‐one abrupt plane contraction are carried out. We will study the behaviour at the entrance corner of the four‐to‐one planar abrupt contraction. Using this solver, we show convergence up to a Weissenberg number We of 20 for the Oldroyd‐B model. No limiting Weissenberg number is observed even though a Phan Thien–Tanner model is used. Several numerical results are presented. Smooth and stable solutions are obtained for high Weissenberg number. Copyright © 2002 John Wiley & Sons, Ltd.     

A volume-tracking method for incompressible multifluid flows with large density variations

A numerical technique (FGVT) for solving the time-dependent incompressible Navier–Stokes equations in fluid flows with large density variations is presented for staggered grids. Mass conservation is based on a volume tracking method and incorporates a piecewise-linear interface reconstruction on a grid twice as fine as the velocity–pressure grid. It also uses a special flux-corrected transport algorithm for momentum advection, a multigrid algorithm for solving a pressure-correction equation and a surface tension algorithm that is robust and stable. In principle, the method conserves both mass and momentum exactly, and maintains extremely sharp fluid interfaces. Applications of the numerical method to prediction of two-dimensional bubble rise in an inclined channel and a bubble bursting through an interface are presented. © 1998 John Wiley & Sons, Ltd.     

An implicit two‐dimensional non‐hydrostatic model for free‐surface flows

An implicit finite volume model in sigma coordinate system is developed to simulate two‐dimensional (2D) vertical free surface flows, deploying a non‐hydrostatic pressure distribution. The algorithm is based on a projection method which solves the complete 2D Navier–Stokes equations in two steps. First the pressure term in the momentum equations is excluded and the resultant advection–diffusion equations are solved. In the second step the continuity and the momentum equation with only the pressure terms are solved to give a block tri‐diagonal system of equation with pressure as the unknown. This system can be solved by a direct matrix solver without iteration. A new implicit treatment of non‐hydrostatic pressure, similar to the lower layers is applied to the top layer which makes the model free of any hydrostatic pressure assumption all through the water column. This treatment enables the model to evaluate both free surface elevation and wave celerity more accurately. A series of numerical tests including free‐surface flows with significant vertical accelerations and nonlinear behaviour in shoaling zone are performed. Comparison between numerical results, analytical solutions and experimental data demonstrates a satisfactory performance. Copyright © 2006 John Wiley & Sons, Ltd.     

A Godunov‐type fractional semi‐implicit method based on staggered grid for dam‐break modeling

A semi‐implicit finite volume model based upon staggered grid is presented for solving shallow water equation. The model employs a time‐splitting scheme that uses a predictor–corrector method for the advection term. The fluxes are calculated based on a Riemann solver in the prediction step and a downwind scheme in the correction step. A simple TVD scheme is employed for shock capturing purposes in which the Minmond limiter is used for flux functions. As a consequence of using staggered grid, an ADI method is adopted for solving the discretized equations for 2‐D problems. Several 1‐D and 2‐D flows have been modeled with satisfactory results when compared with analytical and experimental test cases. The model is also capable of simulating supercritical as well as subcritical flow. Copyright © 2010 John Wiley & Sons, Ltd.     

A pressure–velocity coupling approach for high void fraction free surface bubbly flows in overset curvilinear grids

A methodology for improved robustness in the simulation of high void fraction free surface polydisperse bubbly flows in curvilinear overset grids is presented. The method is fully two‐way coupled in the sense that the bubbly field affects the continuous fluid and vice versa. A hybrid projection approach is used in which staggered contravariant velocities at cell faces are computed for transport and pressure–velocity coupling while the momentum equation is solved on a collocated grid arrangement. Conservation of mass is formulated such that a strong coupling between void fraction, pressure, and velocity is achieved within a partitioned approach, solving each field separately. A pressure–velocity projection solver is iterated together with a predictor stage for the void fraction to achieve a robust coupling. The implementation is described for general curvilinear grids detailing particulars in the neighborhood to overset interfaces or a free surface. A balanced forced method to avoid the generation of spurious currents is extended for curvilinear grids. The overall methodology allows simulation of high void fraction flows and is stable even when strong packing forces accounting for bubble collisions are included. Convergence and stability in one‐dimensional (1D) and two‐dimensional (2D) configurations is evaluated. Finally, a full‐scale simulation of the bubbly flow around a flat‐bottom boat is performed demonstrating the applicability of the methodology to complex problems of engineering interest. Copyright © 2015 John Wiley & Sons, Ltd.     

A hybrid phase field multiple relaxation time lattice Boltzmann method for the incompressible multiphase flow with large density contrast

A hybrid phase field multiple relaxation time lattice Boltzmann method (LBM) is presented in this paper for simulation of multiphase flows with large density contrast. In the present method, the flow field is solved by a lattice Boltzmann equation. Concurrently, the interface of two fluids is captured by solving the macroscopic Cahn‐Hilliard equation using the upwind scheme. To be specific, for simulation of the flow field, an lattice Boltzmann equation (LBE) model developed in Shao et al. (Physical Review E, 89 (2014), 033309) for consideration of density contrast in the momentum equation is used. Moreover, in the present work, the multiple relaxation time collision operator is applied to this LBE to enable simulation of problems with large viscosity contrast or high Reynolds number. For the interface capturing, instead of solving another set of LBE as in many phase field LBMs, the macroscopic Cahn‐Hilliard equation is directly solved by using a weighted essentially non‐oscillatory scheme. In this way, the present hybrid phase field LBM shares full advantages of the phase field LBM while enhancing numerical stability. The ability of the present method to simulate multiphase flow problems with large density contrast is demonstrated by several numerical examples. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical analysis of density flows in adiabatic two-phase fluids through characteristic finite element method

The purpose of this study is to analyze the density flow in adiabatic two-phase fluids through the characteristic finite element method. The fluids are assumed to be liquids. The equations of conservations of mass and momentum for the adiabatic flows and the Birch–Murnaghan equation of state are employed as the governing equations. The employed finite element method is a combination of the characteristic method and the implicit method. The governing equations are divided into two parts: the advection part and the non-advection part. The characteristic method is applied to the advection part. The Hermite interpolation function, which is based on the complete third-order polynomial interpolation using triangular finite element is employed for the interpolation of both velocity and density. Using the discontinuity conditions, an interface translocation method can be derived. The interface of the two flow densities are interpolated through the third-order spline function, using which the curvature of the interface can be directly computed. For the numerical study, the development of density flow over the Tokyo bay is presented. It is detected out that high density area is abruptly diffused over the whole area. According to the differences in the two densities, various flow patterns are computed.     

An adaptive numerical scheme for solving incompressible 2‐phase and free‐surface flows

In this paper, we present a numerical scheme for solving 2‐phase or free‐surface flows. Here, the interface/free surface is modeled using the level‐set formulation, and the underlying mesh is adapted at each iteration of the flow solver. This adaptation allows us to obtain a precise approximation for the interface/free‐surface location. In addition, it enables us to solve the time‐discretized fluid equation only in the fluid domain in the case of free‐surface problems. Fluids here are considered incompressible. Therefore, their motion is described by the incompressible Navier‐Stokes equation, which is temporally discretized using the method of characteristics and is solved at each time iteration by a first‐order Lagrange‐Galerkin method. The level‐set function representing the interface/free surface satisfies an advection equation that is also solved using the method of characteristics. The algorithm is completed by some intermediate steps like the construction of a convenient initial level‐set function (redistancing) as well as the construction of a convenient flow for the level‐set advection equation. Numerical results are presented for both bifluid and free‐surface problems.     

A mass‐conservative staggered immersed boundary model for solving the shallow water equations on complex geometries

In this work, an approach is proposed for solving the 3D shallow water equations with embedded boundaries that are not aligned with the underlying horizontal Cartesian grid. A hybrid cut‐cell/ghost‐cell method is used together with a direction‐splitting implicit solver: Ghost cells are used for the momentum equations in order to prescribe the correct boundary condition at the immersed boundary, while cut cells are used in the continuity equation in order to conserve mass. The resulting scheme is robust, does not suffer any time step limitation for small cut cells, and conserves fluid mass up to machine precision. Moreover, the solver displays a second‐order spatial accuracy, both globally and locally. Comparisons with analytical solutions and reference numerical solutions on curvilinear grids confirm the quality of the method. Copyright © 2015 John Wiley & Sons, Ltd.     

Comparison of implicit and explicit AUSM‐family schemes for compressible multiphase flows

This paper makes the first attempt of extending implicit AUSM‐family schemes to multiphase flow simulations. Water faucet, air–water shock tube and oscillating manometer problems are used as benchmark tests with the generic four‐equation two‐fluid model. For solving the equations implicitly, Newton's method along with a sparse matrix solver (UMFPACK solver) is employed, and the numerical Jacobian matrix is calculated. Comparison between implicit and explicit AUSM‐family schemes is presented, indicating that similarly accurate results are obtained with both schemes. Furthermore, the water faucet problem is solved using both staggered and collocated grids. This investigation helps integrate high‐resolution schemes into staggered‐grid‐based computational algorithms. The influence of the interface pressure correction on the simulation results is also examined. Results show that the interfacial pressure correction introduces numerical dissipation. However, this dissipation cannot eliminate the overshoots because of the incompatibility of numerical discretization of the conservative and non‐conservative terms in the governing equations. The comparison of CPU time between implicit and explicit schemes is also studied, indicating that the implicit scheme is capable of improving the computational efficiency over its explicit counterpart. Copyright © 2014 John Wiley & Sons, Ltd.     

A divergence-free semi-implicit finite volume scheme for ideal,viscous, and resistive magnetohydrodynamics

In this paper, we present a novel pressure-based semi-implicit finite volume solver for the equations of compressible ideal, viscous, and resistive magnetohydrodynamics (MHD). The new method is conservative for mass, momentum, and total energy, and in multiple space dimensions, it is constructed in such a way as to respect the divergence-free condition of the magnetic field exactly, also in the presence of resistive effects. This is possible via the use of multidimensional Riemann solvers on an appropriately staggered grid for the time evolution of the magnetic field and a double curl formulation of the resistive terms. The new semi-implicit method for the MHD equations proposed here discretizes the nonlinear convective terms as well as the time evolution of the magnetic field explicitly, whereas all terms related to the pressure in the momentum equation and the total energy equation are discretized implicitly, making again use of a properly staggered grid for pressure and velocity. Inserting the discrete momentum equation into the discrete energy equation then yields a mildly nonlinear symmetric and positive definite algebraic system for the pressure as the only unknown, which can be efficiently solved with the (nested) Newton method of Casulli et al. The pressure system becomes linear when the specific internal energy is a linear function of the pressure. The time step of the scheme is restricted by a CFL condition based only on the fluid velocity and the Alfvén wave speed and is not based on the speed of the magnetosonic waves. Being a semi-implicit pressure-based scheme, our new method is therefore particularly well suited for low Mach number flows and for the incompressible limit of the MHD equations, for which it is well known that explicit density-based Godunov-type finite volume solvers become increasingly inefficient and inaccurate because of the more and more stringent CFL condition and the wrong scaling of the numerical viscosity in the incompressible limit. We show a relevant MHD test problem in the low Mach number regime where the new semi-implicit algorithm is a factor of 50 faster than a traditional explicit finite volume method, which is a very significant gain in terms of computational efficiency. However, our numerical results confirm that our new method performs well also for classical MHD test cases with strong shocks. In this sense, our new scheme is a true all Mach number flow solver.     

A hybrid vertex-centered finite volume/element method for viscous incompressible flows on non-staggered unstructured meshes

This paper proposes a hybrid vertex-centered finite volume/finite element method for solution of the two dimensional (2D) incompressible Navier-Stokes equations on unstructured grids.An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling.The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by joining the centroid of cells sharing the common vertex.For the temporal integration of the momentum equations,an implicit second-order scheme is utilized to enhance the computational stability and eliminate the time step limit due to the diffusion term.The momentum equations are discretized by the vertex-centered finite volume method (FVM) and the pressure Poisson equation is solved by the Galerkin finite element method (FEM).The momentum interpolation is used to damp out the spurious pressure wiggles.The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both velocity and pressure.The classic test cases,the lid-driven cavity flow,the skew cavity flow and the backward-facing step flow,show that numerical results are in good agreement with the published benchmark solutions.     

Convergence acceleration of segregated algorithms using dynamic tuning additive correction multigrid strategy

A convergence acceleration method based on an additive correction multigrid–SIMPLEC (ACM‐S) algorithm with dynamic tuning of the relaxation factors is presented. In the ACM‐S method, the coarse grid velocity correction components obtained from the mass conservation (velocity potential) correction equation are included into the fine grid momentum equations before the coarse grid momentum correction equations are formed using the additive correction methodology. Therefore, the coupling between the momentum and mass conservation equations is obtained on the coarse grid, while maintaining the segregated structure of the single grid algorithm. This allows the use of the same solver (smoother) on the coarse grid. For turbulent flows with heat transfer, additional scalar equations are solved outside of the momentum–mass conservation equations loop. The convergence of the additional scalar equations is accelerated using a dynamic tuning of the relaxation factors. Both a relative error (RE) scheme and a local Reynolds/Peclet (ER/P) relaxation scheme methods are used. These methodologies are tested for laminar isothermal flows and turbulent flows with heat transfer over geometrically complex two‐ and three‐dimensional configurations. Savings up to 57% in CPU time are obtained for complex geometric domains representative of practical engineering problems. Copyright © 1999 John Wiley & Sons, Ltd.     

Finite element solution to passive scalar transport behind line sources under neutral and unstable stratification

This study employed a direct numerical simulation (DNS) technique to contrast the plume behaviours and mixing of passive scalar emitted from line sources (aligned with the spanwise direction) in neutrally and unstably stratified open‐channel flows. The DNS model was developed using the Galerkin finite element method (FEM) employing trilinear brick elements with equal‐order interpolating polynomials that solved the momentum and continuity equations, together with conservation of energy and mass equations in incompressible flow. The second‐order accurate fractional‐step method was used to handle the implicit velocity–pressure coupling in incompressible flow. It also segregated the solution to the advection and diffusion terms, which were then integrated in time, respectively, by the explicit third‐order accurate Runge–Kutta method and the implicit second‐order accurate Crank–Nicolson method. The buoyancy term under unstable stratification was integrated in time explicitly by the first‐order accurate Euler method. The DNS FEM model calculated the scalar‐plume development and the mean plume path. In particular, it calculated the plume meandering in the wall‐normal direction under unstable stratification that agreed well with the laboratory and field measurements, as well as previous modelling results available in literature. Copyright © 2005 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.     

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.