Absorbing boundary technique for open channel flows

An absorbing boundary condition is formulated and applied to the one‐dimensional open channel flow equations in conjunction with an explicit MacCormack scheme. The physical flow domain has been truncated by introducing an artificial pseudo‐boundary. By using an appropriate boundary condition on a truncated domain, it is shown that, for flow containing shocks, the solution can be accelerated to its stationary profile with no loss of accuracy. Copyright © 2000 John Wiley & Sons, Ltd.     

A confrontation of 1D and 2D RKDG numerical simulation of transitional flow at open‐channel junction

In this study, a comparison between the 1D and 2D numerical simulation of transitional flow in open‐channel networks is presented and completely described allowing for a full comprehension of the modeling water flow. For flow in an open‐channel network, mutual effects exist among the channel branches joining at a junction. Therefore, for the 1D study, the whole system (branches and junction) cannot be treated individually. The 1D Saint Venant equations calculating the flow in the branches are then supplemented by various equations used at the junction: a discharge flow conservation equation between the branches arriving and leaving the junction, and a momentum or energy conservation equation. The disadvantages of the 1D study are that the equations used at the junction are of empirical nature due to certain parameters given by experimental results and moreover they often present a reduced field of validity. On the contrary, for the 2D study, the entire network is considered as a single unit and the flow in all the branches and junctions is solved simultaneously. Therefore, we simply apply the 2D Saint Venant equations, which are solved by a second‐order Runge–Kutta discontinuous Galerkin method. Finally, the experimental results obtained by Hager are used to validate and to compare the two approaches 1D and 2D. Copyright © 2009 John Wiley & Sons, Ltd.     

Flux difference splitting for 1D open channel flow equations

An upwind finite difference scheme based on flux difference splitting is presented for the solution of the equations governing unsteady open channel hydraulics. An approximate Jacobian needed for splitting the flux differences is defined that satisfies the conditions required to construct a first-order upwind conservative discretization of the equations. Added limited second-order corrections make the resulting scheme robust and accurate for the computation of all regimes of open channel flow. Some numerical results and comparisons with other classical schemes under exacting conditions are presented.     

Central schemes for open‐channel flow

The resolution of the Saint‐Venant equations for modelling shock phenomena in open‐channel flow by using the second‐order central schemes of Nessyahu and Tadmor (NT) and Kurganov and Tadmor (KT) is presented. The performances of the two schemes that we have extended to the non‐homogeneous case and that of the classical first‐order Lax–Friedrichs (LF) scheme in predicting dam‐break and hydraulic jumps in rectangular open channels are investigated on the basis of different numerical and physical conditions. The efficiency and robustness of the schemes are tested by comparing model results with analytical or experimental solutions. Copyright © 2003 John Wiley & Sons, Ltd.     

Explicit formulations for evaluation of velocity gradients using boundary‐domain integral equations in 2D and 3D viscous flows

In this paper, explicit boundary‐domain integral equations for evaluating velocity gradients are derived from the basic velocity integral equations. A free term is produced in the new strongly singular integral equation, which is not included in recent formulations using the complex variable differentiation method (CVDM) to compute velocity gradients (Int. J. Numer. Meth. Fluids 2004; 45 :463–484; Int. J. Numer. Meth. Fluids 2005; 47 :19–43). The strongly singular domain integrals involved in the new integral equations are accurately evaluated using the radial integration method (RIM). Considerable computational time for evaluating integrals of velocity gradients can be saved by using present formulation than using CVDM. The formulation derived in this paper together with those presented in reference (Int. J. Numer. Meth. Fluids 2004; 45 :463–484) for 2D and in (Int. J. Numer. Meth. Fluids 2005; 47 :19–43) for 3D problems constitutes a complete boundary‐domain integral equation system for solving full Navier–Stokes equations using primitive variables. Three numerical examples for steady incompressible viscous flow are given to validate the derived formulations. Copyright © 2007 John Wiley & Sons, Ltd.     

Analytical solutions for transverse distributions of stream-wise velocity in turbulent flow in rectangular channel with partial vegetation

The theory of poroelasticity is introduced to study the hydraulic properties of the steady uniform turbulent flow in a partially vegetated rectangular channel. Plants are assumed as immovable media. The resistance caused by vegetation is expressed by the theory of poroelasticity. Considering the influence of a secondary flow, the momentum equation can be simplified. The momentum equation is nondimensionalized to obtain a smooth solution for the lateral distribution of the longitudinal velocity. To verify the model, an acoustic Doppler velocimeter (ADV) is used to measure the velocity field in a rectangular open channel partially with emergent artificial rigid vegetation. Comparisons between the measured data and the computed results show that the method can predict the transverse distributions of stream-wise velocities in turbulent flows in a rectangular channel with partial vegetation.     

SIMPLE‐type preconditioners for cell‐centered,colocated finite volume discretization of incompressible Reynolds‐averaged Navier–Stokes equations

This paper contains a comparison of four SIMPLE‐type methods used as solver and as preconditioner for the iterative solution of the (Reynolds‐averaged) Navier–Stokes equations, discretized with a finite volume method for cell‐centered, colocated variables on unstructured grids. A matrix‐free implementation is presented, and special attention is given to the treatment of the stabilization matrix to maintain a compact stencil suitable for unstructured grids. We find SIMPLER preconditioning to be robust and efficient for academic test cases and industrial test cases. Compared with the classical SIMPLE solver, SIMPLER preconditioning reduces the number of nonlinear iterations by a factor 5–20 and the CPU time by a factor 2–5 depending on the case. The flow around a ship hull at Reynolds number 2E9, for example, on a grid with cell aspect ratio up to 1:1E6, can be computed in 3 instead of 15 h.Copyright © 2012 John Wiley & Sons, Ltd.     

Development of a generalized multi‐layer model for 3‐D simulation of free surface flows

A mathematical model was developed for three‐dimensional (3‐D) simulation of free surface flows. In this model, the flow depth is divided into a number of layers and shallow water equations are integrated in each layer to derive the hydrodynamic equations. To give a general form to this model, each layer is assumed to be non‐horizontal with varying thickness in the flow domain. A non‐orthogonal curvilinear coordinate system is employed in the model, to allow for flexibility in dealing with the irregular geometry of natural watercourses. Due to the similarity in governing equations, two‐dimensional (2‐D) depth averaged programs can be developed into a multi‐layer model. The development for a depth averaged program and its numerical scheme is described in this paper. Experimental data and semi‐analytical solutions are used to evaluate the performance of the model. Three different cases of open channel flow are tested: 1‐flow in a straight open channel, 2‐the flow development region in a channel, and 3‐flow in a meandering channel. It is shown that the model has the capability to predict velocity distribution and secondary flows in complex 3‐D flow conditions. Copyright © 2004 John Wiley & Sons, Ltd.