Development of level set method with good area preservation to predict interface in two‐phase flows

A two‐step conservative level set method is proposed in this study to simulate the gas/water two‐phase flow. For the sake of accuracy, the spatial derivative terms in the equations of motion for an incompressible fluid flow are approximated by the coupled compact scheme. For accurately predicting the modified level set function, the dispersion‐relation‐preserving advection scheme is developed to preserve the theoretical dispersion relation for the first‐order derivative terms shown in the pure advection equation cast in conservative form. For the purpose of retaining its long‐time accurate Casimir functionals and Hamiltonian in the transport equation for the level set function, the time derivative term is discretized by the sixth‐order accurate symplectic Runge–Kutta scheme. To resolve contact discontinuity oscillations near interface, nonlinear compression flux term and artificial damping term are properly added to the second‐step equation of the modified level set method. For the verification of the proposed dispersion‐relation‐preserving scheme applied in non‐staggered grids for solving the incompressible flow equations, three benchmark problems have been chosen in this study. The conservative level set method with area‐preserving property proposed for capturing the interface in incompressible fluid flows is also verified by solving the dam‐break, Rayleigh–Taylor instability, bubble rising in water, and droplet falling in water problems. Good agreements with the referenced solutions are demonstrated in all the investigated problems. Copyright © 2010 John Wiley & Sons, Ltd.     

Exactly well‐balanced positivity preserving nonstaggered central scheme for open‐channel flows

In this paper, we construct and study an exactly well‐balanced positivity‐preserving nonstaggered central scheme for shallow water flows in open channels with irregular geometry and nonflat bottom topography. We introduce a novel discretization of the source term based on hydrostatic reconstruction to obtain the exactly well‐balanced property for the still water steady‐state solution even in the presence of wetting and drying transitions. The positivity‐preserving property of the cross‐sectional wet area is obtained by using a modified “draining" time‐step technique. The current scheme is also Riemann‐solver‐free. Several classical problems of open‐channel flows are used to test these properties. Numerical results confirm that the current scheme is robust, exactly well‐balanced and positivity‐preserving.     

An unstructured finite volume model for dam‐break floods with wet/dry fronts over complex topography

A robust, well‐balanced, unstructured, Godunov‐type finite volume model has been developed in order to simulate two‐dimensional dam‐break floods over complex topography with wetting and drying. The model is based on the nonlinear shallow water equations in hyperbolic conservation form. The inviscid fluxes are calculated using the HLLC approximate Riemann solver and a second‐order spatial accuracy is achieved by implementing the MUSCL reconstruction technique. To prevent numerical oscillations near shocks, slope‐limiting techniques are used for controlling the total variation of the reconstructed field. The model utilizes an explicit two‐stage Runge–Kutta method for time stepping, whereas implicit treatments for friction source terms. The novelties of the model include the flux correction terms and the water depth reconstruction method both for partially and fully submerged cells, and the wet/dry front treatments. The proposed flux correction terms combined with the water depth reconstruction method are necessary to balance the bed slope terms and flux gradient in the hydrostatical steady flow condition. Especially, this well‐balanced property is also preserved in partially submerged cells. It is found that the developed wet/dry front treatments and implicit scheme for friction source terms are stable. The model is tested against benchmark problems, laboratory experimental data, and realistic application related to dam‐break flood wave propagation over arbitrary topography. Numerical results show that the model performs satisfactorily with respect to its effectiveness and robustness and thus has bright application prospects. Copyright © 2010 John Wiley & Sons, Ltd.     

Genuinely multidimensional characteristic‐based scheme for incompressible flows

This paper presents a novel multidimensional characteristic‐based (MCB) upwind method for the solution of incompressible Navier–Stokes equations. As opposed to the conventional characteristic‐based (CB) schemes, it is genuinely multidimensional in that the local characteristic paths, along which information is propagated, are used. For the first time, the multidimensional characteristic structure of incompressible flows modified by artificial compressibility is extracted and used to construct an inherent multidimensional upwind scheme. The new proposed MCB scheme in conjunction with the finite‐volume discretization is employed to model the convective fluxes. Using this formulation, the steady two‐dimensional incompressible flow in a lid‐driven cavity is solved for a wide range of Reynolds numbers. It was found that the new proposed scheme presents more accurate results than the conventional CB scheme in both their first‐ and second‐order counterparts in the case of cavity flow. Also, results obtained with second‐order MCB scheme in some cases are more accurate than the central scheme that in turn provides exact second‐order discretization in this grid. With this inherent upwinding technique for evaluating convective fluxes at cell interfaces, no artificial viscosity is required even at high Reynolds numbers. Another remarkable advantage of MCB scheme lies in its faster convergence rate with respect to the CB scheme that is found to exhibit substantial delays in convergence reported in the literature. The results obtained using new proposed scheme are in good agreement with the standard benchmark solutions in the literature. Copyright © 2007 John Wiley & Sons, Ltd.     

An eigenvector‐based linear reconstruction scheme for the shallow‐water equations on two‐dimensional unstructured meshes

This paper presents a new approach to MUSCL reconstruction for solving the shallow‐water equations on two‐dimensional unstructured meshes. The approach takes advantage of the particular structure of the shallow‐water equations. Indeed, their hyperbolic nature allows the flow variables to be expressed as a linear combination of the eigenvectors of the system. The particularity of the shallow‐water equations is that the coefficients of this combination only depend upon the water depth. Reconstructing only the water depth with second‐order accuracy and using only a first‐order reconstruction for the flow velocity proves to be as accurate as the classical MUSCL approach. The method also appears to be more robust in cases with very strong depth gradients such as the propagation of a wave on a dry bed. Since only one reconstruction is needed (against three reconstructions in the MUSCL approach) the EVR method is shown to be 1.4–5 times as fast as the classical MUSCL scheme, depending on the computational application. Copyright © 2006 John Wiley & Sons, Ltd.     

The composite finite volume method on unstructured meshes for the two‐dimensional shallow water equations

A composite finite volume method (FVM) is developed on unstructured triangular meshes and tested for the two‐dimensional free‐surface flow equations. The methodology is based on the theory of the remainder effect of finite difference schemes and the property that the numerical dissipation and dispersion of the schemes are compensated by each other in a composite scheme. The composite FVM is formed by global composition of several Lax–Wendroff‐type steps followed by a diffusive Lax–Friedrich‐type step, which filters out the oscillations around shocks typical for the Lax–Wendroff scheme. To test the efficiency and reliability of the present method, five typical problems of discontinuous solutions of two‐dimensional shallow water are solved. The numerical results show that the proposed method, which needs no use of a limiter function, is easy to implement, is accurate, robust and is highly stable. Copyright © 2001 John Wiley & Sons, Ltd.     

Finite‐volume multi‐stage schemes for shallow‐water flow simulations

A finite‐volume multi‐stage (FMUSTA) scheme is proposed for simulating the free‐surface shallow‐water flows with the hydraulic shocks. On the basis of the multi‐stage (MUSTA) method, the original Riemann problem is transformed to an independent MUSTA mesh. The local Lax–Friedrichs scheme is then adopted for solving the solution of the Riemann problem at the cell interface on the MUSTA mesh. The resulting first‐order monotonic FMUSTA scheme, which does not require the use of the eigenstructure and the special treatment of entropy fixes, has the generality as well as simplicity. In order to achieve the high‐resolution property, the monotonic upstream schemes for conservation laws (MUSCL) method are used. For modeling shallow‐water flows with source terms, the surface gradient method (SGM) is adopted. The proposed schemes are verified using the simulations of six shallow‐water problems, including the 1D idealized dam breaking, the steady transcritical flow over a hump, the 2D oblique hydraulic jump, the circular dam breaking and two dam‐break experiments. The simulated results by the proposed schemes are in satisfactory agreement with the exact solutions and experimental data. It is demonstrated that the proposed FMUSTA schemes have superior overall numerical accuracy among the schemes tested such as the commonly adopted Roe and HLL schemes. Copyright © 2007 John Wiley & Sons, Ltd.     

Solution of the 2‐D shallow water equations with source terms in surface elevation splitting form

A vertex‐centred finite‐volume/finite‐element method (FV/FEM) is developed for solving 2‐D shallow water equations (SWEs) with source terms written in a surface elevation splitting form, which balances the flux gradients and source terms. The method is implemented on unstructured grids and the numerical scheme is based on a second‐order MUSCL‐like upwind Godunov FV discretization for inviscid fluxes and a classical Galerkin FE discretization for the viscous gradients and source terms. The main advantages are: (1) the discretization of SWE written in surface elevation splitting form satisfies the exact conservation property (??‐Property) naturally; (2) the simple centred‐type discretization can be used for the source terms; (3) the method is suitable for both steady and unsteady shallow water problems; and (4) complex topography can be handled based on unstructured grids. The accuracy of the method was verified for both steady and unsteady problems, including discontinuous cases. The results indicate that the new method is accurate, simple, and robust. Copyright © 2007 John Wiley & Sons, Ltd.     

High‐accuracy upwind method using improved characteristics speeds for incompressible flows

In this study, the Nervier–Stokes equations for incompressible flows, modified by the artificial compressibility method, are investigated numerically. To calculate the convective fluxes, a new high‐accuracy characteristics‐based (HACB) scheme is presented in this paper. Comparing the HACB scheme with the original characteristic‐based method, it is found that the new proposed scheme is more accurate and has faster convergence rate than the older one. The second order averaging scheme is used for estimating the viscose fluxes, and spatially discretized equations are integrated in time by an explicit fourth‐order Runge–Kutta scheme. The lid driven cavity flow and flow in channel with a backward facing step have been used as benchmark problems. It is shown that the obtained results using HACB scheme are in good agreement with the standard solutions. Copyright © 2015 John Wiley & Sons, Ltd.     

Accuracy of high‐order density‐based compressible methods in low Mach vortical flows

At low Mach numbers, Godunov‐type approaches, based on the method of lines, suffer from an accuracy problem. This paper shows the importance of using the low Mach number correction in Godunov‐type methods for simulations involving low Mach numbers by utilising a new, well‐posed, two‐dimensional, two‐mode Kelvin–Helmholtz test case. Four independent codes have been used, enabling the examination of several numerical schemes. The second‐order and fifth‐order accurate Godunov‐type methods show that the vortex‐pairing process can be captured on a low resolution with the low Mach number correction applied down to 0.002. The results are compared without the low Mach number correction and also three other methods, a Lagrange‐remap method, a fifth‐order accurate in space and time finite difference type method based on the wave propagation algorithm, and fifth‐order spatial and third‐order temporal accurate finite volume Monotone Upwind Scheme for Conservation Laws (MUSCL) approach based on the Godunov method and Simple Low Dissipation Advection Upstream Splitting Method (SLAU) numerical flux with low Mach capture property. The ability of the compressible flow solver of the commercial software, ANSYS FLUENT , in solving low Mach flows is also demonstrated for the two time‐stepping methods provided in the compressible flow solver, implicit and explicit. Results demonstrate clearly that a low Mach correction is required for all algorithms except the Lagrange‐remap approach, where dissipation is independent of Mach number. © 2013 Crown copyright. International Journal for Numerical Methods in Fluids. © 2013 John Wiley & Sons, Ltd.     

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

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

Numerical solution of three‐dimensional velocity–vorticity Navier–Stokes equations by finite difference method

This paper describes the finite difference numerical procedure for solving velocity–vorticity form of the Navier–Stokes equations in three dimensions. The velocity Poisson equations are made parabolic using the false‐transient technique and are solved along with the vorticity transport equations. The parabolic velocity Poisson equations are advanced in time using the alternating direction implicit (ADI) procedure and are solved along with the continuity equation for velocities, thus ensuring a divergence‐free velocity field. The vorticity transport equations in conservative form are solved using the second‐order accurate Adams–Bashforth central difference scheme in order to assure divergence‐free vorticity field in three dimensions. The velocity and vorticity Cartesian components are discretized using a central difference scheme on a staggered grid for accuracy reasons. The application of the ADI procedure for the parabolic velocity Poisson equations along with the continuity equation results in diagonally dominant tri‐diagonal matrix equations. Thus the explicit method for the vorticity equations and the tri‐diagonal matrix algorithm for the Poisson equations combine to give a simplified numerical scheme for solving three‐dimensional problems, which otherwise requires enormous computational effort. For three‐dimensional‐driven cavity flow predictions, the present method is found to be efficient and accurate for the Reynolds number range 100?Re?2000. Copyright © 2004 John Wiley & Sons, Ltd.     

High‐order implicit Runge–Kutta time integrators for fluid‐structure interactions

This paper presents an approach to develop high‐order, temporally accurate, finite element approximations of fluid‐structure interaction (FSI) problems. The proposed numerical method uses an implicit monolithic formulation in which the same implicit Runge–Kutta (IRK) temporal integrator is used for the incompressible flow, the structural equations undergoing large displacements, and the coupling terms at the fluid‐solid interface. In this context of stiff interaction problems, the fully implicit one‐step approach presented is an original alternative to traditional multistep or explicit one‐step finite element approaches. The numerical scheme takes advantage of an arbitrary Lagrangian–Eulerian formulation of the equations designed to satisfy the geometric conservation law and to guarantee that the high‐order temporal accuracy of the IRK time integrators observed on fixed meshes is preserved on arbitrary Lagrangian–Eulerian deforming meshes. A thorough review of the literature reveals that in most previous works, high‐order time accuracy (higher than second order) is seldom achieved for FSI problems. We present thorough time‐step refinement studies for a rigid oscillating‐airfoil on deforming meshes to confirm the time accuracy on the extracted aerodynamics reactions of IRK time integrators up to fifth order. Efficiency of the proposed approach is then tested on a stiff FSI problem of flow‐induced vibrations of a flexible strip. The time‐step refinement studies indicate the following: stability of the proposed approach is always observed even with large time step and spurious oscillations on the structure are avoided without added damping. While higher order IRK schemes require more memory than classical schemes (implicit Euler), they are faster for a given level of temporal accuracy in two dimensions. Copyright © 2015 John Wiley & Sons, Ltd.     

Multigrid third‐order least‐squares solution of Cauchy–Riemann equations on unstructured triangular grids

In this paper, a multigrid algorithm is developed for the third‐order accurate solution of Cauchy–Riemann equations discretized in the cell‐vertex finite‐volume fashion: the solution values stored at vertices and the residuals defined on triangular elements. On triangular grids, this results in a highly overdetermined problem, and therefore we consider its solution that minimizes the residuals in the least‐squares norm. The standard second‐order least‐squares scheme is extended to third‐order by adding a high‐order correction term in the residual. The resulting high‐order method is shown to give sufficiently accurate solutions on relatively coarse grids. Combined with a multigrid technique, the method then becomes a highly accurate and efficient solver. We present some results to demonstrate its accuracy and efficiency, including both structured and unstructured triangular grids. Copyright © 2006 John Wiley & Sons, Ltd.     

High‐order compact numerical schemes for non‐hydrostatic free surface flows

The development of a numerical scheme for non‐hydrostatic free surface flows is described with the objective of improving the resolution characteristics of existing solution methods. The model uses a high‐order compact finite difference method for spatial discretization on a collocated grid and the standard, explicit, single step, four‐stage, fourth‐order Runge–Kutta method for temporal discretization. The Cartesian coordinate system was used. The model requires the solution of two Poisson equations at each time‐step and tridiagonal matrices for each derivative at each of the four stages in a time‐step. Third‐ and fourth‐order accurate boundaries for the flow variables have been developed including the top non‐hydrostatic pressure boundary. The results demonstrate that numerical dissipation which has been a problem with many similar models that are second‐order accurate is practically eliminated. A high accuracy is obtained for the flow variables including the non‐hydrostatic pressure. The accuracy of the model has been tested in numerical experiments. In all cases where analytical solutions are available, both phase errors and amplitude errors are very small. Copyright © 2006 John Wiley & Sons, Ltd.     

Vorticity–velocity formulation of the 3D Navier–Stokes equations in cylindrical co‐ordinates

A finite difference method is presented for solving the 3D Navier–Stokes equations in vorticity–velocity form. The method involves solving the vorticity transport equations in ‘curl‐form’ along with a set of Cauchy–Riemann type equations for the velocity. The equations are formulated in cylindrical co‐ordinates and discretized using a staggered grid arrangement. The discretized Cauchy–Riemann type equations are overdetermined and their solution is accomplished by employing a conjugate gradient method on the normal equations. The vorticity transport equations are solved in time using a semi‐implicit Crank–Nicolson/Adams–Bashforth scheme combined with a second‐order accurate spatial discretization scheme. Special emphasis is put on the treatment of the polar singularity. Numerical results of axisymmetric as well as non‐axisymmetric flows in a pipe and in a closed cylinder are presented. Comparison with measurements are carried out for the axisymmetric flow cases. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical implementation of the Crank–Nicolson/Adams–Bashforth scheme for the time‐dependent Navier–Stokes equations

This article considers numerical implementation of the Crank–Nicolson/Adams–Bashforth scheme for the two‐dimensional non‐stationary Navier–Stokes equations. A finite element method is applied for the spatial approximation of the velocity and pressure. The time discretization is based on the Crank–Nicolson scheme for the linear term and the explicit Adams–Bashforth scheme for the nonlinear term. Comparison with other methods, through a series of numerical experiments, shows that this method is almost unconditionally stable and convergent, i.e. stable and convergent when the time step is smaller than a given constant. Copyright © 2009 John Wiley & Sons, Ltd.     

The Lagrange–Galerkin method for the two‐dimensional shallow water equations on adaptive grids

The weak Lagrange–Galerkin finite element method for the two‐dimensional shallow water equations on adaptive unstructured grids is presented. The equations are written in conservation form and the domains are discretized using triangular elements. Lagrangian methods integrate the governing equations along the characteristic curves, thus being well suited for resolving the non‐linearities introduced by the advection operator of the fluid dynamics equations. An additional fortuitous consequence of using Lagrangian methods is that the resulting spatial operator is self‐adjoint, thereby justifying the use of a Galerkin formulation; this formulation has been proven to be optimal for such differential operators. The weak Lagrange–Galerkin method automatically takes into account the dilation of the control volume, thereby resulting in a conservative scheme. The use of linear triangular elements permits the construction of accurate (by virtue of the second‐order spatial and temporal accuracies of the scheme) and efficient (by virtue of the less stringent Courant–Friedrich–Lewy (CFL) condition of Lagrangian methods) schemes on adaptive unstructured triangular grids. Lagrangian methods are natural candidates for use with adaptive unstructured grids because the resolution of the grid can be increased without having to decrease the time step in order to satisfy stability. An advancing front adaptive unstructured triangular mesh generator is presented. The highlight of this algorithm is that the weak Lagrange–Galerkin method is used to project the conservation variables from the old mesh onto the newly adapted mesh. In addition, two new schemes for computing the characteristic curves are presented: a composite mid‐point rule and a general family of Runge–Kutta schemes. Results for the two‐dimensional advection equation with and without time‐dependent velocity fields are illustrated to confirm the accuracy of the particle trajectories. Results for the two‐dimensional shallow water equations on a non‐linear soliton wave are presented to illustrate the power and flexibility of this strategy. Copyright © 2000 John Wiley & Sons, Ltd.     

A well‐balanced explicit/semi‐implicit finite element scheme for shallow water equations in drying–wetting areas

Numerical solutions of the shallow water equations can be used to reproduce flow hydrodynamics occurring in a wide range of regions. In hydraulic engineering, the objectives include the prediction of dam break wave propagation, fluvial floods and other catastrophic flooding phenomena, the modeling of estuarine and coastal circulations, and the design and optimization of hydraulic structures. In this paper, a well‐balanced explicit and semi‐implicit finite element scheme for shallow water equations over complex domains involving wetting and drying is proposed. The governing equations are discretized by a fractional finite element method using a two‐step Taylor–Galerkin scheme. First, the intermediate increment of conserved variable is obtained explicitly neglecting the pressure gradient term. This is then corrected for the effects of pressure once the pressure increment has been obtained from the Poisson equation. In order to maintain the ‘well‐balanced’ property, the pressure gradient term and bed slope terms are incorporated into the Poisson equation. Moreover, a local bed slope modification technique is employed in drying–wetting interface treatments. The proposed model is well validated against several theoretical benchmark tests. Copyright © 2014 John Wiley & Sons, Ltd.     

Time‐accurate stabilized finite‐element model for weakly nonlinear and weakly dispersive water waves

Introduction of a time‐accurate stabilized finite‐element approximation for the numerical investigation of weakly nonlinear and weakly dispersive water waves is presented in this paper. To make the time approximation match the order of accuracy of the spatial representation of the linear triangular elements by the Galerkin finite‐element method, the fourth‐order time integration of implicit multistage Padé method is used for the development of the numerical scheme. The streamline‐upwind Petrov–Galerkin (SUPG) method with crosswind diffusion is employed to stabilize the scheme and suppress the spurious oscillations, usually common in the numerical computation of convection‐dominated flow problems. The performance of numerical stabilization and accuracy is addressed. Treatments of various boundary conditions, including the open boundary conditions, the perfect reflecting boundary conditions along boundaries with irregular geometry, are also described. Numerical results showing the comparisons with analytical solutions, experimental measurements, and other published numerical results are presented and discussed. Copyright © 2007 John Wiley & Sons, Ltd.