A two‐dimensional finite volume morphodynamic model on unstructured triangular grids

We discuss the application of a finite volume method to morphodynamic models on unstructured triangular meshes. The model is based on coupling the shallow water equations for the hydrodynamics with a sediment transport equation for the morphodynamics. The finite volume method is formulated for the quasi‐steady approach and the coupled approach. In the first approach, the steady hydrodynamic state is calculated first and the corresponding water velocity is used in the sediment transport equation to be solved subsequently. The second approach solves the coupled hydrodynamics and sediment transport system within the same time step. The gradient fluxes are discretized using a modified Roe's scheme incorporating the sign of the Jacobian matrix in the morphodynamic system. A well‐balanced discretization is used for the treatment of source terms. We also describe an adaptive procedure in the finite volume method by monitoring the bed–load in the computational domain during its transport process. The method uses unstructured meshes, incorporates upwinded numerical fluxes and slope limiters to provide sharp resolution of steep bed gradients that may form in the approximate solution. Numerical results are shown for a test problem in the evolution of an initially hump‐shaped bed in a squared channel. For the considered morphodynamical regimes, the obtained results point out that the coupled approach performs better than the quasi‐steady approach only when the bed–load rapidly interacts with the hydrodynamics. Copyright © 2009 John Wiley & Sons, Ltd.     

A high‐order accurate method for two‐dimensional incompressible viscous flows

A high‐order accurate solution method for complex geometries is developed for two‐dimensional flows using the stream function–vorticity formulation. High‐order accurate spectrally optimized compact schemes along with appropriate boundary schemes are used for spatial discretization while a two‐level backward Euler implicit scheme is used for the time integration. The linear system of equations for stream function and vorticity are solved by an inner iteration while contravariant velocities constitute outer iterations. The effect of curvilinear grids on the solution accuracy is studied. The method is used to compute Cartesian and inclined driven cavity, flow in a triangular cavity and viscous flow in constricted channel. Benchmark‐like accuracy is obtained in all the problems with fewer grid points compared to reported studies. Copyright © 2006 John Wiley & Sons, Ltd.     

Pseudostress–velocity formulation for incompressible Navier–Stokes equations

This paper presents a numerical algorithm using the pseudostress–velocity formulation to solve incompressible Newtonian flows. The pseudostress–velocity formulation is a variation of the stress–velocity formulation, which does not require symmetric tensor spaces in the finite element discretization. Hence its discretization is greatly simplified. The discrete system is further decoupled into an H ( div ) problem for the pseudostress and a post‐process resolving the velocity. This can be done conveniently by using the penalty method for steady‐state flows or by using the time discretization for nonsteady‐state flows. We apply this formulation to the 2D lid‐driven cavity problem and study its grid convergence rate. Also, computational results of the time‐dependent‐driven cavity problem and the flow past rectangular problem are reported. Copyright © 2009 John Wiley & Sons, Ltd.     

A comparative study of high‐order variable‐property segregated algorithms for unsteady low Mach number flows

In this paper, five different algorithms are presented for the simulation of low Mach flows with large temperature variations, based on second‐order central‐difference or fourth‐order compact spatial discretization and a pressure projection‐type method. A semi‐implicit three‐step Runge–Kutta/Crank–Nicolson or second‐order iterative scheme is used for time integration. The different algorithms solve the coupled set of governing scalar equations in a decoupled segregate manner. In the first algorithm, a temperature equation is solved and density is calculated from the equation of state, while the second algorithm advances the density using the differential form of the equation of state. The third algorithm solves the continuity equation and the fourth algorithm solves both the continuity and enthalpy equation in conservative form. An iterative decoupled algorithm is also proposed, which allows the computation of the fully coupled solution. All five algorithms solve the momentum equation in conservative form and use a constant‐ or variable‐coefficient Poisson equation for the pressure. The efficiency of the fourth‐order compact scheme and the performances of the decoupling algorithms are demonstrated in three flow problems with large temperature variations: non‐Boussinesq natural convection, channel flow instability, flame–vortex interaction. 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.     

A comparison of time integration methods in an unsteady low‐Reynolds‐number flow

This paper describes three different time integration methods for unsteady incompressible Navier–Stokes equations. Explicit Euler and fractional‐step Adams–Bashford methods are compared with an implicit three‐level method based on a steady‐state SIMPLE method. The implicit solver employs a dual time stepping and an iteration within the time step. The spatial discretization is based on a co‐located finite‐volume technique. The influence of the convergence limits and the time‐step size on the accuracy of the predictions are studied. The efficiency of the different solvers is compared in a vortex‐shedding flow over a cylinder in the Reynolds number range of 100–1600. A high‐Reynolds‐number flow over a biconvex airfoil profile is also computed. The computations are performed in two dimensions. At the low‐Reynolds‐number range the explicit methods appear to be faster by a factor from 5 to 10. In the high‐Reynolds‐number case, the explicit Adams–Bashford method and the implicit method appear to be approximately equally fast while yielding similar results. Copyright © 2002 John Wiley & Sons, Ltd.     

An overlapping control volume method for the Navier–Stokes equations on non‐staggered grids

An algorithm, based on the overlapping control volume (OCV) method, for the solution of the steady and unsteady two‐dimensional incompressible Navier–Stokes equations in complex geometry is presented. The primitive variable formulation is solved on a non‐staggered grid arrangement. The problem of pressure–velocity decoupling is circumvented by using momentum interpolation. The accuracy and effectiveness of the method is established by solving five steady state and one unsteady test problems. The numerical solutions obtained using the technique are in good agreement with the analytical and benchmark solutions available in the literature. On uniform grids, the method gives second‐order accuracy for both diffusion‐ and convection‐dominated flows. There is little loss of accuracy on grids that are moderately non‐orthogonal. Copyright © 1999 John Wiley & Sons, Ltd.     

Application of a preconditioned high‐order accurate artificial compressibility‐based incompressible flow solver in wide range of Reynolds numbers

In the present study, the preconditioned incompressible Navier‐Stokes equations with the artificial compressibility method formulated in the generalized curvilinear coordinates are numerically solved by using a high‐order compact finite‐difference scheme for accurately and efficiently computing the incompressible flows in a wide range of Reynolds numbers. A fourth‐order compact finite‐difference scheme is utilized to accurately discretize the spatial derivative terms of the governing equations, and the time integration is carried out based on the dual time‐stepping method. The capability of the proposed solution methodology for the computations of the steady and unsteady incompressible viscous flows from very low to high Reynolds numbers is investigated through the simulation of different 2‐dimensional benchmark problems, and the results obtained are compared with the existing analytical, numerical, and experimental data. A sensitivity analysis is also performed to evaluate the effects of the size of the computational domain and other numerical parameters on the accuracy and performance of the solution algorithm. The present solution procedure is also extended to 3 dimensions and applied for computing the incompressible flow over a sphere. Indications are that the application of the preconditioning in the solution algorithm together with the high‐order discretization method in the generalized curvilinear coordinates provides an accurate and robust solution method for simulating the incompressible flows over practical geometries in a wide range of Reynolds numbers including the creeping flows.     

A non‐iterative implicit algorithm for the solution of advection–diffusion equation on a sphere

A numerical algorithm for the solution of advection–diffusion equation on the surface of a sphere is suggested. The velocity field on a sphere is assumed to be known and non‐divergent. The discretization of advection–diffusion equation in space is carried out with the help of the finite volume method, and the Gauss theorem is applied to each grid cell. For the discretization in time, the symmetrized double‐cycle componentwise splitting method and the Crank–Nicolson scheme are used. The numerical scheme is of second order approximation in space and time, correctly describes the balance of mass of substance in the forced and dissipative discrete system and is unconditionally stable. In the absence of external forcing and dissipation, the total mass and L₂‐norm of solution of discrete system is conserved in time. The one‐dimensional periodic problems arising at splitting in the longitudinal direction are solved with Sherman–Morrison's formula and Thomas's algorithm. The one‐dimensional problems arising at splitting in the latitudinal direction are solved by the bordering method that requires a prior determination of the solution at the poles. The resulting linear systems have tridiagonal matrices and are solved by Thomas's algorithm. The suggested method is direct (without iterations) and rapid in realization. It can also be applied to linear and nonlinear diffusion problems, some elliptic problems and adjoint advection–diffusion problems on a sphere. Copyright © 2015 John Wiley & Sons, Ltd.     

Numerical simulation of minimum B‐jumps at abrupt drops

A numerical simulation of minimum B‐jumps in horizontal rectangular channels having an abrupt drop is given. Gradually varied, steady, supercritical flow is assumed as the initial condition. An unsteady flow is created by increasing the downstream depth. One‐dimensional, unsteady Saint‐Venant equations are solved by using the MacCormack and the dissipative two–four explicit finite difference schemes. The steady flow solution is obtained by treating the time variable as an iteration parameter and letting the solution converge to the steady state. The abrupt drop is treated as an interior boundary and solved by the method of characteristics. The results are compared with experimental and analytical studies. Copyright © 2007 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.     

Preserving bounded and conservative solutions of transport in one‐dimensional shallow‐water flow with upwind numerical schemes: Application to fertigation and solute transport in rivers

This work intends to show that conservative upwind schemes based on a separate discretization of the scalar solute transport from the shallow‐water equations are unable to preserve uniform solute profiles in situations of one‐dimensional unsteady subcritical flow. However, the coupled discretization of the system is proved to lead to the correct solution in first‐order approximations. This work is also devoted to show that, when using a coupled discretization, a careful definition of the flux limiter function in second‐order TVD schemes is required in order to preserve uniform solute profiles. The work shows that, in cases of subcritical irregular flow, the coupled discretization is necessary but nevertheless not sufficient to ensure concentration distributions free from oscillations and a method to avoid these oscillations is proposed. Examples of steady and unsteady flows in test cases, river and irrigation are presented. Copyright © 2007 John Wiley & Sons, Ltd.     

A high‐order discontinuous Galerkin solver for low Mach number flows

In this work, we present a high‐order discontinuous Galerkin method (DGM) for simulating variable density flows at low Mach numbers. The corresponding low Mach number equations are an approximation of the compressible Navier–Stokes equations in the limit of zero Mach number. To the best of the authors'y knowledge, it is the first time that the DGM is applied to the low Mach number equations. The mixed‐order formulation is applied for spatial discretization. For steady cases, we apply the semi‐implicit method for pressure‐linked equation (SIMPLE) algorithm to solve the non‐linear system in a segregated manner. For unsteady cases, the solver is implicit in time using backward differentiation formulae, and the SIMPLE algorithm is applied to solve the non‐linear system in each time step. Numerical results for the following three test cases are shown: Couette flow with a vertical temperature gradient, natural convection in a square cavity, and unsteady natural convection in a tall cavity. Considering a fixed number of degrees of freedom, the results demonstrate the benefits of using higher approximation orders. Copyright © 2015 John Wiley & Sons, Ltd.     

On higher‐order mixed FEM for low Mach number flows: application to a natural convection benchmark problem

We consider higher‐order mixed finite elements with continuous pressures for the computation of stationary compressible flows at low Mach number. The proposed approach is based on a fully coupled treatment of the governing equations and therefore, for steady‐state calculations, does not rely on time‐stepping techniques. The non‐linear problem is solved by means of a quasi‐Newton iteration. The strongly coupled system resulting from higher‐order discretization of the linearized equations requires adequate solvers. We propose a new scheme based on multigrid methods with varying FEM ansatz orders on the grid hierarchy as well as multiplicative smoothers based on blocking techniques. Computational results are described for a benchmark configuration including a flow with heat transfer in the low Mach number regime. Furthermore, the issue of anisotropic grids is addressed in that context. Copyright © 2003 John Wiley & Sons, Ltd.     

High‐order numerical simulations of flow‐induced noise

In this paper, the flow/acoustics splitting method for predicting flow‐generated noise is further developed by introducing high‐order finite difference schemes. The splitting method consists of dividing the acoustic problem into a viscous incompressible flow part and an inviscid acoustic part. The incompressible flow equations are solved by a second‐order finite volume code EllipSys2D/3D. The acoustic field is obtained by solving a set of acoustic perturbation equations forced by flow quantities. The incompressible pressure and velocity form the input to the acoustic equations. The present work is an extension of our acoustics solver, with the introduction of high‐order schemes for spatial discretization and a Runge–Kutta scheme for time integration. To achieve low dissipation and dispersion errors, either Dispersion‐Relation‐Preserving (DRP) schemes or optimized compact finite difference schemes are used for the spatial discretizations. Applications and validations of the new acoustics solver are presented for benchmark aeroacoustic problems and for flow over an NACA 0012 airfoil. Copyright © 2010 John Wiley & Sons, Ltd.     

Implicit weighted essentially non‐oscillatory schemes for the incompressible Navier–Stokes equations

A class of lower–upper/approximate factorization (LUAF) implicit weighted essentially non‐oscillatory (ENO; WENO) schemes for solving the two‐dimensional incompressible Navier–Stokes equations in a generalized co‐ordinate system is presented. The algorithm is based on the artificial compressibility formulation, and symmetric Gauss–Seidel relaxation is used for computing steady state solutions while symmetric successive overrelaxation is used for treating time‐dependent flows. WENO spatial operators are employed for inviscid fluxes and central differencing for viscous fluxes. Internal and external viscous flow test problems are presented to verify the numerical schemes. The use of a WENO spatial operator not only enhances the accuracy of solutions but also improves the convergence rate for the steady state computation as compared with using the ENO counterpart. It is found that the present solutions compare well with exact solutions, experimental data and other numerical results. Copyright © 1999 John Wiley & Sons, Ltd.     

Three‐dimensional numerical modelling of free surface flows with non‐hydrostatic pressure

A three‐dimensional numerical model is developed for incompressible free surface flows. The model is based on the unsteady Reynolds‐averaged Navier–Stokes equations with a non‐hydrostatic pressure distribution being incorporated in the model. The governing equations are solved in the conventional sigma co‐ordinate system, with a semi‐implicit time discretization. A fractional step method is used to enable the pressure to be decomposed into its hydrostatic and hydrodynamic components. At every time step one five‐diagonal system of equations is solved to compute the water elevations and then the hydrodynamic pressure is determined from a pressure Poisson equation. The model is applied to three examples to simulate unsteady free surface flows where non‐hydrostatic pressures have a considerable effect on the velocity field. Emphasis is focused on applying the model to wave problems. Two of the examples are about modelling small amplitude waves where the hydrostatic approximation and long wave theory are not valid. The other example is the wind‐induced circulation in a closed basin. The numerical solutions are compared with the available analytical solutions for small amplitude wave theory and very good agreement is obtained. Copyright © 2002 John Wiley & Sons, Ltd.     

A Fourier–Chebyshev spectral collocation method for simulating flow past spheres and spheroids

An accurate Fourier–Chebyshev spectral collocation method has been developed for simulating flow past prolate spheroids. The incompressible Navier–Stokes equations are transformed to the prolate spheroidal co‐ordinate system and discretized on an orthogonal body fitted mesh. The infinite flow domain is truncated to a finite extent and a Chebyshev discretization is used in the wall‐normal direction. The azimuthal direction is periodic and a conventional Fourier expansion is used in this direction. The other wall‐tangential direction requires special treatment and a restricted Fourier expansion that satisfies the parity conditions across the poles is used. Issues including spatial and temporal discretization, efficient inversion of the pressure Poisson equation, outflow boundary condition and stability restriction at the pole are discussed. The solver has been validated primarily by simulating steady and unsteady flow past a sphere at various Reynolds numbers and comparing key quantities with corresponding data from experiments and other numerical simulations. Copyright © 1999 John Wiley & Sons, Ltd.     

Unsteady three‐dimensional boundary layer flow due to a stretching surface in a micropolar fluid

In this paper, the unsteady three‐dimensional boundary layer flow due to a stretching surface in a viscous and incompressible micropolar fluid is considered. The partial differential equations governing the unsteady laminar boundary layer flow are solved numerically using an implicit finite‐difference scheme. The numerical solutions are obtained which are uniformly valid for all dimensionless time from initial unsteady‐state flow to final steady‐state flow in the whole spatial region. The equations for the initial unsteady‐state flow are also solved analytically. It is found that there is a smooth transition from the small‐time solution to the large‐time solution. The features of the flow for different values of the governing parameters are analyzed and discussed. The solutions of interest for the skin friction coefficient with various values of the stretching parameter c and material parameter K are presented. Copyright © 2011 John Wiley & Sons, Ltd.