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.     

Development of quasi two-dimensional structures in a shallow free-surface mixing layer

 The large quasi two-dimensional turbulence structures that emerge in a shallow mixing layer are studied experimentally using Laser Doppler Anemometry. Velocity profiles and turbulence intensities are measured in the first two meters downstream of the splitter plate. In contradistinction with previous experiments, it is shown that the initial growth rate and the turbulence intensities of the shallow mixing layer compare well with deep-water plane mixing layers. Two-point measurements allowed for the determination of spatial correlations of the fluctuating velocity components. The large eddies were found to extend from one tenth of the water depth up to the free surface while the streamwise size of the eddies was found to be three times the mixing layer width. The two-dimensional character of the large structures and the associated reversed energy cascade is inferred from the power spectra of the lateral velocity component. Received: 2 April 1997/Accepted: 25 August 1997     

A stabilized SPH method for inviscid shallow water flows

In this paper, the smoothed particle hydrodynamics (SPH) method is applied to the solution of shallow water equations. A brief review of the method in its standard form is first described then a variational formulation using SPH interpolation is discussed. A new technique based on the Riemann solver is introduced to improve the stability of the method. This technique leads to better results. The treatment of solid boundary conditions is discussed but remains an open problem for general geometries. The dam‐break problem with a flat bed is used as a benchmark test. Copyright © 2004 John Wiley & Sons, Ltd.     

A POD-based reduced-order model for uncertainty analyses in shallow water flows

ABSTRACT

This work presents an intrusive reduced-order model (IROM) for uncertainty propagation analyses for flood flows. The 2D shallow water equations are reduced using Galerkin’s projection onto bases obtained from the snapshot-based proper orthogonal decomposition technique. To speed up the computations, the non-polynomial and nonlinear momentum and friction terms are judiciously approximated and the time accuracy issues are addressed using the principal interval decomposition technique. The performance of the IROM is investigated in some test cases. Also, this model is applied to the study of uncertainty propagation for a hypothetical flood in a real river, to derive a probabilistic flood map. The upstream discharge and the Manning roughness coefficient are considered as the uncertain parameters. For relatively small variations around the mean of the inputs, the comparisons of the statistical moments (mean and standard deviation) of the water depth show errors, between the reduced and full models, less than 0.72%. These simulations were completed at up to 50 times faster using the proposed reduced model.     

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.     

High-resolution numerical model for shallow water flows and pollutant diffusions

ＩｎｔｒｏｄｕｃｔｉｏｎＳｈａｌｌｏｗｗａｔｅｒｆｌｏｗｅｘｉｓｔｓｖｅｒｙｗｉｄｅｌｙｉｎｈｙｄｒａｕｌｉｃｅｎｇｉｎｅｅｒｉｎｇ .Ｐｏｌｌｕｔａｎｔｄｉｆｆｕｓｉｏｎｔｈａｔｉｓｍａｉｎｌｙｃｏｎｔｒｏｌｌｅｄｂｙｔｈｅｓｈａｌｌｏｗｗａｔｅｒｆｌｏｗｈａｓｂｒｏｕｇｈｔｉｎｃｒｅａｓｉｎｇｌｙａｔｔｅｎｔｉｏｎｉｎｅｎｖｉｒｏｎｍｅｎｔａｌｅｎｇｉｎｅｅｒｉｎｇ.Ｔｈｅｔｗｏ＿ｄｉｍｅｎｓｉｏｎａｌｍｏｄｅｌｆｏｒｄｅｓｃｒｉｂｉｎｇｔｈｅｆｌｏｗｃｈａｒａｃｔｅｒｉｓｔｉｃｓｈａｓｂｅｅ…     

Over-reflection and instabilities in shear flows of shallow water with bottom friction

In a two-dimensional shear flow of shallow water, the bottom friction relates uniquely the spanwise profile of the depth-averaged velocity to the bottom topography. If the basic flow varies weakly in the spanwise direction, the local analysis of stability at every spanwise position gives the region of the flow parameters for which the classic hydraulic instability due to the bottom friction cannot occur. In this region, the linear analyses of the waves scattering and instability due to the lateral shear can be performed effectively by means of the frictionless linearized equations if both the bottom slope and friction are equally small.The energy of the total perturbed flow can be split into three main parts that correspond to the basic flow, small amplitude wave motion and induced mean flow. The waves can be either amplified or damped near the critical layers, where their streamwise phase velocity equals the velocity of the basic flow. Two physical mechanisms of this amplification exist. The first one is similar to that suggested by Takehiro and Hayashi for a linear frictionless shallow water flow. The incident and transmitted waves carry energy of opposite signs, which results in an increase in the amplitude of the reflected wave compared to that of the incident one. This mechanism of over-reflection operates for any combination of the flow parameters. The other mechanism is similar to Landau damping in plasma flows; it is related to the energy exchange between the waves and fluid particles at the critical layers due to the velocity synchronism. It may lead to either additional amplification or damping of the waves for different flow conditions. In particular, its significance can be reduced by stronger bottom friction. If the basic flow has uniform potential vorticity, Landau damping is negligible, and over-reflection always occurs. If the feed-back is provided by another critical layer, the net over-reflection results in the formation of trapped modes.     

Boundary layer structure of oscillatory open-channel shallow flows over smooth and rough beds

The boundary layer structure of oscillatory shallow open channel flows has been studied in a wide flume. Fluorescence solution was released at a porous rough bed through a diffuser covered by gravel of 0.5 cm grain size. A planar laser-induced fluorescence (PLIF) system was used to visualise the dye plumes in both vertical and horizontal planes for a qualitative understanding of the roles of large-scale flow structures in mass transport. A variety of tests were conducted for a range of oscillatory periods (30–240 s), water depths (3–16 cm) and velocity amplitudes (0.027–0.325 m/s), which cover a wide range of oscillatory flows with Reynolds numbers Re _a varied from 0.3 × 10⁴ (laminar) to 2.1 × 10⁶ (fully turbulent). For quantitative investigation, a novel technique, namely combined laser-induced fluorescence (LIF) and 2D laser Doppler velocimetry (LDV) (LIF/LDV), was developed and used to measure the velocity and solute concentration simultaneously in a vertical plane over 50 cycles. From the dye plumes revealed by the PLIF in transitional flows, there are different patterns of flow structure and solute transport with three representative stages of acceleration, deceleration and flow reversal. In the acceleration stage, turbulence was suppressed with dye layers adhering to the surface with little vertical mass transport. In the deceleration stage, flame-like turbulent structures occurred when turbulence generation was prominent. This was investigated quantitatively by recording the percentage occurrence of the adhered smooth layers per cycle. For those smooth bed cases with Re _a < 1.8 × 10⁵, the adhered smooth dye layers type of boundary layer occupied 100% of the oscillation period. Over a sufficiently high Re _a , a rough bed can generate fully turbulent oscillatory flows without the appearance of adhering dye layers. Between these two extremes, a transitional flow regime occurs in a wide range of flow conditions: Re _a > 2.7 × 10⁴ over the rough bed and Re _a > 8.3 × 10⁶ over a smooth bed.     

An hp‐adaptive discontinuous Galerkin method for shallow water flows

An adaptive spectral/hp discontinuous Galerkin method for the two‐dimensional shallow water equations is presented. The model uses an orthogonal modal basis of arbitrary polynomial order p defined on unstructured, possibly non‐conforming, triangular elements for the spatial discretization. Based on a simple error indicator constructed by the solutions of approximation order p and p?1, we allow both for the mesh size, h, and polynomial approximation order to dynamically change during the simulation. For the h‐type refinement, the parent element is subdivided into four similar sibling elements. The time‐stepping is performed using a third‐order Runge–Kutta scheme. The performance of the hp‐adaptivity is illustrated for several test cases. It is found that for the case of smooth flows, p‐adaptivity is more efficient than h‐adaptivity with respect to degrees of freedom and computational time. Copyright © 2010 John Wiley & Sons, Ltd.     

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.     

A Non Homogeneous Riemann Solver for shallow water and two phase flows

In this work we consider a two steps finite volume scheme, recently developed to solve nonhomogeneous systems. The first step of the scheme depends on a diffusion control parameter which we modulate, using the limiters theory. Results on Shallow water equations and two phase flows are presented.     

A staggered conservative scheme for every Froude number in rapidly varied shallow water flows

This paper proposes a numerical technique that in essence is based upon the classical staggered grids and implicit numerical integration schemes, but that can be applied to problems that include rapidly varied flows as well. Rapidly varied flows occur, for instance, in hydraulic jumps and bores. Inundation of dry land implies sudden flow transitions due to obstacles such as road banks. Near such transitions the grid resolution is often low compared to the gradients of the bathymetry. In combination with the local invalidity of the hydrostatic pressure assumption, conservation properties become crucial. The scheme described here, combines the efficiency of staggered grids with conservation properties so as to ensure accurate results for rapidly varied flows, as well as in expansions as in contractions. In flow expansions, a numerical approximation is applied that is consistent with the momentum principle. In flow contractions, a numerical approximation is applied that is consistent with the Bernoulli equation. Both approximations are consistent with the shallow water equations, so under sufficiently smooth conditions they converge to the same solution. The resulting method is very efficient for the simulation of large‐scale inundations. Copyright © 2003 John Wiley & Sons, Ltd.     

Investigation of turbulent flows via pseudo flow visualization part I: Axisymmetric jet mixing layer

Hot-wire anemometer measurements obtained in the near-field axisymmetric jet mixing layer by Glauser and George [1] are examined using a pseudo flow visualization (PFV) technique. Pseudo flow visualization is a visualization procedure used to manipulate data obtained from an array of probes to create a graphical representation of the instantaneous and fluctuating velocity components of a flow field. An indicator function was employed to identify the frequency content of each velocity-time trace, giving insight into the analysis of the visualizations. From this application, the natural shedding frequency, or preferred mode, of the large-scale structures was determined and compared with the conventional streamwise and radial spectral measurements acquired by Glauser and George [1]. Furthermore, the wavelength of the preferred mode, nondimensionalized by the jet exit diameter, was determined to be approximately 2.4, a result consistent with the work of Crowe and Champagne [2]. In Part 1 the technique is developed and discussed for the fundamental and fairly well-researched mixing layer of the axisymmetric jet. Our aim is to verify the effectiveness of PFV in the context of a well-documented flow. In Part 2, this technique is then applied to an industrial flow field, namely, the mixing region of a lobed mixer     

Numerical study on the stability and mixing of vertical round buoyant jet in shallow water

IntroductionAsimplewaytodisposeofthelargequantitiesofwasteheatresultingfromsteamelectricpowergeneration ,istodischargetheheatedcondenserwaterthroughasubmergedroundoutfalllocatedatthebottomofthereceivingwater.Inrecentyears,suchkindsofhighvelocityinjectionsystemshaveincreasinglybeenusedaseffectivemixingdevices,frequently ,thesedevicesarelocatedinwatersofonlyafewportdiametersdeep ,anddischargingareaisrelativelysmall.Thusanunderstandingoftheinducedexcesstemperaturedistribution ,inrelationshiptothe…     

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.