Lattice Boltzmann modeling of shallow water flows over discontinuous beds

Numerical modeling of shallow water flows over discontinuous beds is presented. The flows are described with the shallow water equations and the equations are solved using the lattice Boltzmann method (LBM) with single relaxation time (Bhatnagar–Gross–Krook‐LBM (BGK‐LBM)) and the multiple relaxation time (MRT‐LBM). The weighted centered scheme for force term together with the bed height for a bed slope is described to improve simulation of flows over discontinuous bed. Furthermore, the resistance stress is added to include the local head loss caused by flow over a step. Four test cases, one‐dimensional tidal over regular bed and steps, dam‐break flows, and two‐dimensional shallow water flow over a square block, are considered to verify the present method. Agreements between predictions and analytical solutions are satisfactory. Furthermore, the performance and CPU cost time of BGK‐LBM and MRT‐LBM are compared and studied. The results have shown that the lattice Boltzmann method is simple and accurate for simulating shallow water flows over discontinuous beds. This demonstrates the capability and applicability of the lattice Boltzmann method in modeling shallow water flows on bed topography with a discontinuity in practical hydraulic engineering. Copyright © 2014 John Wiley & Sons, Ltd.     

A spectral approach for the stability analysis of turbulent open-channel flows over granular beds

A novel Orr–Sommerfeld-like equation for gravity-driven turbulent open-channel flows over a granular erodible bed is here derived, and the linear stability analysis is developed. The whole spectrum of eigenvalues and eigenvectors of the complete generalized eigenvalue problem is computed and analyzed. The fourth-order eigenvalue problem presents singular non-polynomial coefficients with non-homogenous Robin-type boundary conditions that involve first and second derivatives. Furthermore, the Exner condition is imposed at an internal point. We propose a numerical discretization of spectral type based on a single-domain Galerkin scheme. In order to manage the presence of singular coefficients, some properties of Jacobi polynomials have been carefully blended with numerical integration of Gauss–Legendre type. The results show a positive agreement with the classical experimental data and allow one to relate the different types of instability to such parameters as the Froude number, wavenumber, and the roughness scale. The eigenfunctions allow two types of boundary layers to be distinguished, scaling, respectively, with the roughness height and the saltation layer for the bedload sediment transport.     

Boundary layer interaction in three-dimensional flows

An approach known from the theory of matched asymptotic expansions involving the isolation of subregions with different scales is used to study flows which are assumed to be described by the boundary layer equations almost everywhere near the surface except for a fairly narrow zone in which the inflowing boundary layers interact. Two characteristic types of interaction are identified. An approximate theory describing the flow in the interaction zone, which makes it possible to locate the position of the interaction zone on the surface, is proposed. The interaction flow on the end wall of a vane channel is calculated subject to certain simplifications. The results of an experimental investigation of this flow are presented and it is shown that the theoretical model proposed describes the three-dimensional corner separation which occurs in the neighborhood of the line of intersection of the end wall and the convex edge of the vane.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 116–123, May–June, 1988.     

Horizontal mixing layer in shallow water flows

Horizontal-shear thin-layer homogeneous fluid flow in the open channel is considered. A one-dimensional mathematical model of the development and evolution of the horizontal mixing layer is derived within the framework of the three-layer scheme. The steady-state solutions of the equations of motion are constructed and investigated. In particular, supercritical (subcritical)-in-average flow concepts are introduced and the problem of the mixing layer structure is solved. The proposed model is verified on the basis of comparison with a numerical solution of two-dimensional equations of shallow water theory.     

Adverse pressure gradient turbulent flows over rough walls

An experimental study of a fully developed turbulent channel flow and an adverse pressure gradient (APG) turbulent channel flow over smooth and rough walls has been performed using a particle image velocimetry (PIV) technique. The rough walls comprised two-dimensional square ribs of nominal height, k = 3 mm and pitch, p = 2k, 4k and 8k. It was observed that rib roughness enhanced the drag characteristics, and the degree of enhancement increased with increasing pitch. Similarly, rib roughness significantly increased the level of turbulence production, Reynolds stresses and wall-normal transport of turbulence kinetic energy and Reynolds shear stress well beyond the roughness sublayer. On the contrary, the distributions of the eddy viscosity, mixing length and streamwise transport of turbulence kinetic energy and Reynolds shear stress were reduced by wall roughness, especially in the outer layer. Adverse pressure gradient produced a further reduction in the mean velocity (in comparison to the results obtained in the parallel section) but increased the wall-normal extent across which the mean flow above the ribs is spatially inhomogeneous in the streamwise direction. APG also reinforced wall roughness in augmenting the equivalent sand grain roughness height. The combination of wall roughness and APG significantly increased turbulence production and Reynolds stresses except in the immediate vicinity of the rough walls. The transport velocities of the turbulence kinetic energy and Reynolds shear stress were also augmented by APG across most part of the rough-wall boundary layer. Further, APG enhanced the distributions of the eddy viscosity across most of the boundary layer but reduced the mixing length outside the roughness sublayer.     

Boundary of monotonic and oscillatory convective stability of a horizontal fluid layer

Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, No. 6, pp. 44–50, November–December, 1991.     

Stability of jet flows in a rotating shallow water layer

The problem of the stability of an isolated jet flow and two counter-streaming jet flows in a rotating shallow-water layer is considered. These flows are described by exact solutions of the Charny–Obukhov equation with one or two discontinuities of the potential vorticity, respectively. The isolated jet flow is shown to be stable. For the system consisting of two jet flows the dependence of the characteristics of the unstable wave modes on a geometric parameter, namely, the ratio of the spacing between the jet axes to the deformation radius, is determined. On the basis of the contour dynamics method a weakly-nonlinear model of the longwave instability is developed.     

Simulation of shallow flows over variable topographies using unstructured grids

Simulation of shallow flows over variable topographies is a challenging case for most available shock‐capturing schemes. This problem arises because the source terms and flux gradients are not balanced in the numerical computations. Treatments for this problem generally work well on structured grids, but they are usually too expensive, and most of them are not directly applicable to unstructured grids. In this paper we propose two efficient methods to treat the source terms without upwinding and to satisfy the compatibility condition on unstructured grids. In the first method, the calculation of the bed slope source term is performed by employing a compatible approximation of water depth at the cell interfaces. In the second one, different components of the bed slope term are considered separately and a compatible discretization of the components is proposed. The present treatments are applicable for most schemes including the Roe's method without changing the performance of the original scheme for smooth topographies. Copyright © 2006 John Wiley & Sons, Ltd.     

A turbulent boundary layer over a two-dimensional rough wall

Particle image velocimetry (PIV) measurements and planar laser induced fluorescence (PLIF) visualizations have been made in a turbulent boundary layer over a rough wall. The wall roughness consisted of square bars placed transversely to the flow at a pitch to height ratio of λ/k = 11 for the PLIF experiments and λ/k = 8 and 16 for the PIV measurements. The ratio between the boundary layer thickness and the roughness height k/δ was about 20 for the PLIF and 38 for the PIV. Both the PLIF and PIV data showed that the near-wall region of the flow was populated by unstable quasi-coherent structures which could be associated to shear layers originating at the trailing edge of the roughness elements. The streamwise mean velocity profile presented a downward shift which varied marginally between the two cases of λ/k, in agreement with previous measurements and DNS results. The data indicated that the Reynolds stresses normalized by the wall units are higher for the case λ/k = 16 than those for λ/k = 8 in the outer region of the flow, suggesting that the roughness density effects could be felt well beyond the near-wall region of the flow. As expected the roughness disturbed dramatically the sublayer which in turn altered the turbulence production mechanism. The turbulence production is maximum at a distance of about 0.5k above the roughness elements. When normalized by the wall units, the turbulence production is found to be smaller than that of a smooth wall. It is argued that the production of turbulence is correlated with the form drag.     

Simple efficient algorithm (SEA) for shallow flows with shock wave on dry and irregular beds

An explicit Godunov‐type solution algorithm called SEA (simple efficient algorithm) has been introduced for the shallow water equations. The algorithm is based on finite volume conservative discretisation method. It can deal with wet/dry and irregular beds. Second‐order accuracy, in both time and space, is achieved using prediction and correction steps. A very simple and efficient flux limiting technique is used to equip the algorithm with total variation dimensioning property for shock capturing purposes. In order to make sure about the balance between the flux gradient and the bed slope, treatment of the source term has been done using a new procedure inspired mainly by the physical rather than mathematical consideration. SEA has been applied to one‐dimensional problems, although it can equally be applied to multi‐dimensional problems. In order to assess the capability of proposed algorithm in dealing with practical applications, several test cases have been examined. Copyright © 2007 John Wiley & Sons, Ltd.     

Experimental and numerical study of the turbulent boundary layer over shallow dimples

A very convincing experimental paper from TU Delft, of 2016, motivated this joint study. They directly measured a reduction in drag of up to 3% for the developed turbulent boundary layer flowing over an array of shallow dimples, of diameter comparable with the boundary-layer thickness. The Reynolds numbers were not high, but nor marginal either. Our experimental set-up has been well validated by previous riblets studies. The simulation approach contains more new features, which are described, and passed the essential tests. It is a Direct Numerical Simulation, at the correct Reynolds number. The key finding is that both of our approaches predict a drag increase of at least 1%, with an uncertainty of the order of 1%. Not having identified any flaw in the 2016 work even through discussions with the authors, we are currently unable to explain the contradiction. We do note however that these initial results had a Reynolds-number dependence which we consider puzzling, and also that a 2018 paper from TU Delft by almost the same set of authors presents some less favourable findings.     

Solving the fully nonlinear weakly dispersive Serre equations for flows over dry beds

We describe a numerical method for solving the Serre equations that can simulate flows over dry bathymetry. The method solves the Serre equations in conservation law form with a finite volume method. A finite element method is used to solve the auxiliary elliptic equation for the depth‐averaged horizontal velocity. The numerical method is validated against the lake at rest analytic solution, demonstrating that it is well‐balanced. Since there are currently no known nonstationary analytical solutions to the Serre equation that involve bathymetry, a nonstationary forced solution, involving bathymetry was developed. The method was further validated and its convergence rate established using the developed nonstationary forced solution containing the wetting and drying of bathymetry. Finally, the method is also validated against experimental results for the run‐up of a solitary wave on a sloped beach. The finite‐volume finite‐element approach to solving the Serre equation was found to be accurate and robust.     

Turbulent boundary layer velocity profiles over a smooth surface

Mean velocity measurements made for two-dimensional incompressible turbulent boundary layers, for 3 pressure gradients on a smooth wall have been used to evaluate the pressure gradient parameter, the skin-friction and the roughness function by the method of “curve-fit” to the entire profile. Local skin-friction coefficients so obtained do not agree with the experimental results but the integral parameters are well represented.     

Turbulent boundary layer flow with a step change from smooth to rough surface

A direct numerical simulation (DNS) dataset of a turbulent boundary layer (TBL) with a step change from a smooth to a rough surface is analyzed to examine the characteristics of a spatially developing flow. The roughness elements are periodically arranged two-dimensional (2-D) spanwise rods, with the first rod placed 80θ_in downstream from the inlet, where θ_in denotes the inlet momentum thickness. Based on an accurate estimation of relevant parameters, clear evidence for mean flow universality is provided when scaled properly, even for the present roughness configuration, which is believed to have one of the strongest impacts on the flow. Compared to previous studies, it is shown that overshooting behavior is present in the first- and second-order statistics and is locally created either within the cavity or at the leading edge of the roughness depending on the type of statistics and the wall-normal measurement location. Inspection of spatial two-point correlations of the streamwise velocity fluctuations shows a continuous increase of spanwise length scales of structures over the rough wall after the step change at a greater growth rate than that over smooth wall TBL flow. This is expected because spanwise energy spectrum shows presence of much energetic wider structures over the rough wall. Full images of the DNS data are presented to describe not only predominance of hairpin vortices but also a possible spanwise scale growth mechanism via merging over the rough wall.     

A modified lattice Boltzmann model for shallow water flows over complex topography

A modified lattice Boltzmann model is proposed to describe shallow water flows over complex topography. In the proposed model, the quadratic depth term is excluded from the equilibrium distribution functions (EDFs), and the hydrostatic pressure term is combined with the bed slope term to be treated as a part of the sourcing term in the lattice Boltzmann equation (LBE). Therefore, it is unnecessary to match the coefficients of the quadratic depth term in the EDFs with those of the bed slope term in the sourcing terms in the LBE. This would bring more flexibility to the treatment of the sourcing terms in the LBE. In order to recover the shallow water equations (SWEs), the basic constraints are redefined, and under these constraints, the coefficients of the EDFs are derived afterwards. Several benchmark problems are used to validate the proposed model, including stationary case, steady flows over a two‐dimensional bump and tidal wave flows over irregular bed elevation. The computed results are in excellent agreement with the results of the other numerical methods and the analytical solutions, indicating that the proposed model is capable of simulating shallow water flows over complex bathymetry. It also proves that the proposed model has potential to produce competitive solutions to shallow water flows over complex bed topography. Copyright © 2014 John Wiley & Sons, Ltd.     

Fluidized and vibrofluidized shallow beds of fresh leaves

The fluid dynamics behavior of shallow fluidized and vibrofluidized beds operating with fresh leaves was investigated with the aim of exploring drying applications in a modified conveyor belt (MCB) system,which may be operated in a fixedor fluidized-bed mode.Leaves of the specimens Duranta repens,Schinus molle,Coleus barbatus,Buxus sempervirens,and Bougainvillea spectabilis were tested with a range of sphericities from 0.063 to 0.213,bulk densities from 0.038 to 0.251 g/cm 3,apparent densities from 0.52 to ...     

Development of a turbulent boundary layer after a step from smooth to rough surface

The flow developing downstream of a step change from smooth to rough surface condition is studied in the light of Townsend’s wall similarity hypothesis. Previous studies seem to support the hypothesis for channel and pipe flows, but there are considerable controversies about its application to boundary layers and in particular to surface roughness formed by spanwise bars. It has been suggested that this controversy arises from insufficient separation of scales between the boundary layer thickness and the roughness length scale. An experimental investigation has therefore been undertaken where the flow evolves from a fully developed smooth wall boundary layer at high Reynolds numbers over a step in surface roughness (Re _θ = 13,400 at the step). The flow is mapped through the development of the internal layer until the flow is fully developed over the rough wall. The internal layer is found to grow as δ ∼ X ^0.73, and after about 15 boundary layer thicknesses at the step, the internal layer has reached the outer edge of the incoming layer. At the last rough wall measurement station, the Reynolds number has grown to Re _θ ≈ 32,600 and the ratio of boundary layer to roughness length scales is δ/k ≈ 140. The outer layer differences between the smooth and the rough wall data were found to be sufficiently small to conclude that for this setup the Townsend’s wall similarity hypothesis appears to hold.