Computation of unsteady viscous incompressible flows in generalized non‐inertial co‐ordinate system using Godunov‐projection method and overlapping meshes

Time‐dependent incompressible Navier–Stokes equations are formulated in generalized non‐inertial co‐ordinate system and numerically solved by using a modified second‐order Godunov‐projection method on a system of overlapped body‐fitted structured grids. The projection method uses a second‐order fractional step scheme in which the momentum equation is solved to obtain the intermediate velocity field which is then projected on to the space of divergence‐free vector fields. The second‐order Godunov method is applied for numerically approximating the non‐linear convection terms in order to provide a robust discretization for simulating flows at high Reynolds number. In order to obtain the pressure field, the pressure Poisson equation is solved. Overlapping grids are used to discretize the flow domain so that the moving‐boundary problem can be solved economically. Numerical results are then presented to demonstrate the performance of this projection method for a variety of unsteady two‐ and three‐dimensional flow problems formulated in the non‐inertial co‐ordinate systems. Copyright © 2002 John Wiley & Sons, Ltd.     

Parallel computation of viscous incompressible flows using Godunov‐projection method on overlapping grids

The Godunov‐projection method is implemented on a system of overlapping structured grids for solving the time‐dependent incompressible Navier–Stokes equations. This projection method uses a second‐order fractional step scheme in which the momentum equation is solved to obtain the intermediate velocity field which is then projected on to the space of divergence‐free vector fields. The Godunov procedure is applied to estimate the non‐linear convective term in order to provide a robust discretization of this terms at high Reynolds number. In order to obtain the pressure field, a separate procedure is applied in this modified Godunov‐projection method, where the pressure Poisson equation is solved. Overlapping grids are used to discretize the flow domain, as they offer the flexibility of simplifying the grid generation around complex geometrical domains. This combination of projection method and overlapping grid is also parallelized and reasonable parallel efficiency is achieved. Numerical results are presented to demonstrate the performance of this combination of the Godunov‐projection method and the overlapping grid. Copyright © 2002 John Wiley & Sons, Ltd.     

Numerical analysis of axisymmetric and planar sudden expansion flows for laminar regime

In this study, the Navier–Stokes equations are solved numerically for axisymmetric and planar sudden expansion flows. The flow is considered laminar and steady state, and the fluid is incompressible. Finite difference equations are obtained using a control volume method in a non‐staggered grid arrangement, and solved by line‐by‐line TDMA technique using the SIMPLEM (SIMPLEM‐Modified) algorithm. Calculations are performed for higher expansion ratios, β, ranging from 1.5 to 10 and Reynolds numbers from 0.1 to 500. Results are presented in terms of streamlines, relative eddy intensity, location of the eddy center, and the eddy reattachment length depending on Re number and β values for both axisymmetric and planar sudden expansions. It is aimed to provide a picture of the effects of high values of expansion ratio and Reynolds number on the sudden expansion flow. As a result, it is found that the flow characteristics keep their structure for both higher expansion ratios and higher Reynolds numbers. Further, correlations were developed for the nondimensional eddy reattachment length, location of the eddy center and the relative eddy intensity, which have agreeable results to the computed results available from the literature. Copyright © 2010 John Wiley & Sons, Ltd.     

Numerical implementation of Aristov–Pukhnachev's formulation for axisymmetric viscous incompressible flows

In the present work a finite‐difference technique is developed for the implementation of a new method proposed by Aristov and Pukhnachev (Doklady Phys. 2004; 49 (2):112–115) for modeling of the axisymmetric viscous incompressible fluid flows. A new function is introduced that is related to the pressure and a system similar to the vorticity/stream function formulation is derived for the cross‐flow. This system is coupled to an equation for the azimuthal velocity component. The scheme and the algorithm treat the equations for the cross‐flow as an inextricably coupled system, which allows one to satisfy two conditions for the stream function with no condition on the auxiliary function. The issue of singularity of the matrix is tackled by adding a small parameter in the boundary conditions. The scheme is thoroughly validated on grids with different resolutions. The new numerical tool is applied to the Taylor flow between concentric rotating cylinders when the upper and lower lids are allowed to rotate independently from the inner cylinder, while the outer cylinder is held at rest. The phenomenology of this flow is adequately represented by the numerical model, including the hysteresis that takes place near certain specific values of the Reynolds number. Thus, the present results can be construed to demonstrate the viability of the new model. The success can be attributed to the adequate physical nature of the auxiliary function. The proposed technique can be used in the future for in‐depth investigations of the bifurcation phenomena in rotating flows. Copyright © 2009 John Wiley & Sons, Ltd.     

A conservative semi‐implicit method for coupled surface–subsurface flows in regional scale

A semi‐implicit method for coupled surface–subsurface flows in regional scale is proposed and analyzed. The flow domain is assumed to have a small vertical scale as compared with the horizontal extents. Thus, after hydrostatic approximation, the simplified governing equations are derived from the Reynolds averaged Navier–Stokes equations for the surface flow and from the Darcy's law for the subsurface flow. A conservative free‐surface equation is derived from a vertical integral of the incompressibility condition and extends to the whole water column including both, the surface and the subsurface, wet domains. Numerically, the horizontal domain is covered by an unstructured orthogonal grid that may include subgrid specifications. Along the vertical direction a simple z‐layer discretization is adopted. Semi‐implicit finite difference equations for velocities and a finite volume approximation for the free‐surface equation are derived in such a fashion that, after simple manipulation, the resulting discrete free‐surface equation yields a single, well‐posed, mildly nonlinear system. This system is efficiently solved by a nested Newton‐type iterative method that yields simultaneously the pressure and a non‐negative fluid volume throughout the computational grid. The time‐step size is not restricted by stability conditions dictated by friction or surface wave speed. The resulting algorithm is simple, extremely efficient, and very accurate. Exact mass conservation is assured also in presence of wetting and drying dynamics, in pressurized flow conditions, and during free‐surface transition through the interface. A few examples illustrate the model applicability and demonstrate the effectiveness of the proposed algorithm. Copyright © 2015 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.     

Numerical simulation of the unsteady behaviour of cavitating flows

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

Assessment of a high-order accurate Discontinuous Galerkin method for turbomachinery flows

In this work the capabilities of a high-order Discontinuous Galerkin (DG) method applied to the computation of turbomachinery flows are investigated. The Reynolds averaged Navier–Stokes equations coupled with the two equations k-ω turbulence model are solved to predict the flow features, either in a fixed or rotating reference frame, to simulate the fluid flow around bodies that operate under an imposed steady rotation. To ensure, by design, the positivity of all thermodynamic variables at a discrete level, a set of primitive variables based on pressure and temperature logarithms is used. The flow fields through the MTU T106A low-pressure turbine cascade and the NASA Rotor 37 axial compressor have been computed up to fourth-order of accuracy and compared to the experimental and numerical data available in the literature.     

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.     

Finite element solution of the streamfunction-vorticity equations for incompressible two-dimensional flows

The streamfunction-vorticity equations for incompressible two-dimensional flows are uncoupled and solved in sequence by the finite element method. The vorticity at no-slip boundaries is evaluated in the framework of the streamfunction equation. The resulting scheme achieves convergence, even for very high values of the Reynolds number, without the traditional need for upwinding. The stability and accuracy of the approach are demonstrated by the solution of two well-known benchmark problems: flow in a lid-driven cavity at Re ? 10,000 and flow over a backward-facing step at Re = 800.     

Numerical study of flows with multiple free surfaces

This paper demonstrates that a numerical method based on the generalized simplified marker and cell (GENSMAC) flow solver and Youngs' volume of fluid (Y‐VOF) surface‐tracking technique is an effective tool for studying the basic mechanics of hydraulic engineering problems with multiple free surfaces and non‐hydrostatic pressure distributions. Two‐dimensional flow equations in a vertical plane are solved numerically for this purpose. The numerical results are compared with experimental data and earlier numerical results based on a higher‐order depth‐averaged flow model available in the literature. Two classical problems, (i) flow in a free overfall and (ii) flow past a floor slot, are considered. The numerical results correspond very well with the experimental data for both sub‐critical and supercritical flows. Copyright © 2001 John Wiley & Sons, Ltd.     

Effects of Reynolds number on physiological‐type pulsatile flows in a pipe with ring‐type constrictions

The effects of Reynolds number on the physiological‐type of laminar pulsatile flow fields within the vicinity of mechanical ring‐type constriction in small pipes were studied numerically. The parameters considered are: the Reynolds number (Re) in the range of 50–1500; Strouhal number (St) in the range of 0.00156–3.98; Womersley number (Nw) from 0.0 to 50.0. The pulsatile flows considered were physiological‐type of simulated flows. Within a pulsating cycle, detailed flow characteristics were studied through the pulsating contours of streamline (ψ), vorticity (Ω), shear stress (τ) and isobar. The relations between the instantaneous flow rate (Q) and instantaneous pressure gradients (dp/dz) are observed to be elliptic. The relations between the instantaneous flow rate (Q) and pressure loss (P_loss) are quadratic. Linear relations were observed between the instantaneous flow rate (Q) and the maximum velocity, maximum vorticity and maximum shear stress. The Reynolds number of the flow in a pulsating cycle was found to have significant effects on the recirculation length and the pressure gradient within the pulsatile flow regime. Copyright © 1999 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.     

Three‐dimensional numerical simulation of compound meandering open channel flow by the Reynolds stress model

Turbulent flow in a compound meandering open channel with seminatural cross sections is one of the most complicated turbulent flows as the flow pattern is influenced by the combined action of various forces, such as centrifugal force, pressure, and shear stresses. In this paper, a three‐dimensional (3D) Reynolds stress model (RSM) is adopted to simulate the compound meandering channel flows. Governing equations of the flow are solved numerically with finite‐volume method. The velocity fields, wall shear stresses, and Reynolds stresses are calculated for a range of input conditions. Good agreement between the simulated results and measurements indicates that RSM can successfully predict the complicated flow phenomenon. Copyright © 2008 John Wiley & Sons, Ltd.     

An implicit numerical algorithm for solving non‐hydrostatic free‐surface flow problems

Details are given of the development of a two‐dimensional vertical numerical model for simulating unsteady free‐surface flows, using a non‐hydrostatic pressure distribution. In this model, the Reynolds equations and the kinematic free‐surface boundary condition are solved simultaneously, so that the water surface elevation can be integrated into the solution and solved for, together with the velocity and pressure fields. An efficient numerical algorithm has been developed, deploying implicit parameters similar to those used in the Crank–Nicholson method, and generating a block tri‐diagonal algebraic system of equations. The model has been applied to simulate a range of unsteady flow problems involving relatively strong vertical accelerations. The results show that the numerical algorithm described is able to produce accurate predictions and is also easy to apply. Copyright © 2001 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.     

The DRBEM solution of incompressible MHD flow equations

This paper presents a dual reciprocity boundary element method (DRBEM) formulation coupled with an implicit backward difference time integration scheme for the solution of the incompressible magnetohydrodynamic (MHD) flow equations. The governing equations are the coupled system of Navier‐Stokes equations and Maxwell's equations of electromagnetics through Ohm's law. We are concerned with a stream function‐vorticity‐magnetic induction‐current density formulation of the full MHD equations in 2D. The stream function and magnetic induction equations which are poisson‐type, are solved by using DRBEM with the fundamental solution of Laplace equation. In the DRBEM solution of the time‐dependent vorticity and current density equations all the terms apart from the Laplace term are treated as nonhomogeneities. The time derivatives are approximated by an implicit backward difference whereas the convective terms are approximated by radial basis functions. The applications are given for the MHD flow, in a square cavity and in a backward‐facing step. The numerical results for the square cavity problem in the presence of a magnetic field are visualized for several values of Reynolds, Hartmann and magnetic Reynolds numbers. The effect of each parameter is analyzed with the graphs presented in terms of stream function, vorticity, current density and magnetic induction contours. Then, we provide the solution of the step flow problem in terms of velocity field, vorticity, current density and magnetic field for increasing values of Hartmann number. Copyright © 2010 John Wiley & Sons, Ltd.     

A comparative study of central and upwind difference schemes using the primitive variables

The use of the velocity-pressure formulation of the Navier-Stokes equations for the numerical solution of fluid flow problems is favoured for free-surface problems, more involved flow configurations, and three-dimensional flows. Many engineering problems involve such features in addition to strong inertial effects. The computational instabilities arising from central-difference schemes for the convective terms of the governing equations impose serious limitations on the range of Reynolds numbers that can be investigated by the numerical method. Solutions for higher Reynolds numbers Re > 1000 could be reached using upwind-difference schemes. A comparative study of both schemes using a method based on the primitive variables is presented. The comparison is made for the model problem of the driven flow in a square cavity. Using a central scheme stable solutions of the pressure and velocity fields were obtained for Reynolds numbers up to 5000. The streamfunction and vorticity fields were calculated from the resulting velocity field and compared with previous solutions. It is concluded that total upwind differencing results in a considerable change in the flow pattern due to the false diffusion. For practical calculations, by a proper choice of a small amount of partial upwind differencing the vorticity diffusion near the walls and the global features of the solutions are not sigificantly altered.     

A VORTICITY–STREAMFUNCTION FORMULATION FOR STEADY INCOMPRESSIBLE TWO-DIMENSIONAL FLOWS

A vorticity–streamfunction formulation for incompressible planar viscous flows is presented. The standard kinematic field equations are discretized using centred finite difference schemes and solved in a coupled way via a Newton-like linearization scheme. The linearized system of partial differential equations is handled through the restarting linear GMRES algorithm, preconditioned by means of an incomplete LU approximate factorization. The proposed solution technique constitutes a fast and robust algorithm for treating laminar flows at high Reynolds numbers. The pressure field is obtained at a subsequent step by solving a convection– diffusion equation in terms of the stagnation pressure, which presents certain advantages compared with the widely used static pressure Poisson equation. Results are shown for a wide variety of applications including internal and external flows.     

A local‐analytic‐based discretization procedure for the numerical solution of incompressible flows

An Erratum has been published for this article in International Journal for Numerical Methods in Fluids 2005, 49(8): 933. We present a local‐analytic‐based discretization procedure for the numerical solution of viscous fluid flows governed by the incompressible Navier–Stokes equations. The general procedure consists of building local interpolants obtained from local analytic solutions of the linear multi‐dimensional advection–diffusion equation, prototypical of the linearized momentum equations. In view of the local analytic behaviour, the resulting computational stencil and coefficient values are functions of the local flow conditions. The velocity–pressure coupling is achieved by a discrete projection method. Numerical examples in the form of well‐established verification and validation benchmarks are presented to demonstrate the capabilities of the formulation. The discretization procedure is implemented alongside the ability to treat embedded and non‐matching grids with relative motion. Of interest are flows at high Reynolds number, ??(10⁵)–??(10⁷), for which the formulation is found to be robust. Applications include flow past a circular cylinder undergoing vortex‐induced vibrations (VIV) at high Reynolds number. Copyright © 2005 John Wiley & Sons, Ltd.