On barotropic and baroclinic tides over an arbitrary sloping topography

Some theoretical concepts about the frictionless dynamics of propagation of the barotropic tide over two-dimensional continental slopes of arbitrary shape are developed. A numerical procedure which generalizes the exact solution obtained over a rectilinear sloping topography is given. This technique can be applied to compute the harmonic components of the barotropic tide everywhere over sloping bottom contours of any shape. It permits in particular the avoidance of discontinuities at the boundaries of rectilinear-continental-slope profiles. The barotropic tidal results are used afterwards to calculate the barotropic forcing for the generation of internal tides. Numerical experiments are performed to study the interaction between the tide and some typical sloping topographies. A three-layered model is used for this purpose. Results are compared with those previously obtained over a rectilinear continental slope.     

SAR浅海水下地形遥感探测技术综述

SAR已成为浅海水下地形探测的重要技术手段之一.与传统浅海水下地形探测技术相比,SAR浅海水下地形遥感探测技术具有明显的经济效益.该水深探测技术通过对浅海水下地形SAR图像仿真模型的反演求解,从SAR图像中提取水下地形信息.本文回顾了SAR浅海水下地形遥感探测技术的不同数值模型和应用实例,并针对目前SAR浅海水下地形遥感探测技术存在的问题和今后研究方向进行了探讨和总结.     

Laboratory experiment on the 3D tide-induced Lagrangian residual current using the PIV technique

The 3D structure of the tide-induced Lagrangian residual current was studied using the particle image velocimetry (PIV) technique in a long shallow narrow tank in the laboratory. At the mouth of the tank, a wave generator was used to make periodic wave which represents the tide movement, and at the head of the tank, a laterally sloping topography with the length of one fifth of the water tank was installed, above which the tide-induced Lagrangian residual current was studied. Under the weakly nonlinear condition in the present experiment setup, the results show that the Lagrangian residual velocity (LRV) field has a three-layer structure. The residual current flows inwards (towards the head) in the bottom layer and flows outwards in the middle layer, while in the surface layer, it flows inwards along the shallow side of the sloping topography and outwards along the deep side. The depth-averaged and breadth-averaged LRV are also analyzed based on the 3D LRV observations. Our results are in good agreement with the previous experiment studies, the analytical solutions with similar conditions and the observational results in real bays. Moreover, the volume flux comparison between the Lagrangian and Eulerian residual currents shows that the Eulerian residual velocity violates the mass conservation law while the LRV truly represents the inter-tidal water transport. This work enriches the laboratory studies of the LRV and offers valuable references for the LRV studies in real bays.     

On the exact shape of the horizontal profile of a topographically rectified tidal flow

Abstract

Starting from the nonlinear shallow water equations of a homogeneous rotating fluid we derive the equation describing the evolution of vorticity by a fluctuating bottom topography of small amplitude, using a multiple scale expansion in a small parameter, which is the topographic length scale relative to the tidal wave length. The exact response functions of residual vorticity for a sinusoidal bottom topography are compared with those obtained by a primitive perturbation series and by harmonic truncation, showing the former to be invalid for small topographic length scales and the latter to be only a fair approximation for vorticity produced by planetary vortex stretching. In deriving the exact shape of the horizontal residual velocity profile at a step-like break in the bottom topography, it is shown that the Lagrangian profile only exists in a strip having the width of the amplitude of the tidal excursion at both sides of the break, and that it vanishes outside that interval. Moreover, in the limit of small amplitude topography at least, it vanishes altogether for the generation mechanism by means of planetary vortex stretching. The Eulerian profile is shown to extend over twice the interval of the Lagrangian profile both for production by vortex stretching and by differential bottom friction. These finite intervals over which the residual velocity profiles exist for a step-like topography are not reproduced by harmonic truncation of the basic equation. This method gives exponentially decaying profiles, indicating spurious horizontal diffusion of vorticity. In terms of orders of magnitude, the method of harmonic truncation is reliable for residual velocity produced by vortex stretching but it overestimates the residual velocity produced by differential bottom friction by a factor 2.     

Well-balanced high-order centred schemes for non-conservative hyperbolic systems. Applications to shallow water equations with fixed and mobile bed

This paper concerns the development of high-order accurate centred schemes for the numerical solution of one-dimensional hyperbolic systems containing non-conservative products and source terms. Combining the PRICE-T method developed in [Toro E, Siviglia A. PRICE: primitive centred schemes for hyperbolic system of equations. Int J Numer Methods Fluids 2003;42:1263–91] with the theoretical insights gained by the recently developed path-conservative schemes [Castro M, Gallardo J, Parés C. High-order finite volume schemes based on reconstruction of states for solving hyperbolic systems with nonconservative products applications to shallow-water systems. Math Comput 2006;75:1103–34; Parés C. Numerical methods for nonconservative hyperbolic systems: a theoretical framework. SIAM J Numer Anal 2006;44:300–21], we propose the new PRICE-C scheme that automatically reduces to a modified conservative FORCE scheme if the underlying PDE system is a conservation law. The resulting first-order accurate centred method is then extended to high order of accuracy in space and time via the ADER approach together with a WENO reconstruction technique. The well-balanced properties of the PRICE-C method are investigated for the shallow water equations. Finally, we apply the new scheme to the shallow water equations with fix bottom topography and with variable bottom solving an additional sediment transport equation.     

Improved bounds on linear instability of barotropic zonal flow within the shallow water equations

Here we develop mathematical results to describe the location of linear instability of a parallel mean flow within the framework of the shallow water equations; growth estimates of near neutral modes (for disturbances subcritical with respect to gravity wave speed) in the cases of non-rotating and rotating shallow water. The bottom topography is taken to be one-dimensional and the isobaths are parallel to the mean flow. In the case of a rotating fluid, the isobaths and the mean flow are assumed to be zonal. The flow is front-like: there is a monotonic increase of mean flow velocity. Our results show that for barotropic flows the location of instabilities will be a semi-ellipse region in the complex wave velocity plane, that is based on the wave-number, Froude number, and depth of the fluid layer. We also explore the instability region for the case of spatially unbounded mean velocity profiles for non-rotating shallow water.     

Three-dimensional numerical modelling of wind-driven circulation in a homogeneous lake

A three-dimensional numerical model has been developed to study wind-induced circulation patterns in a shallow homogeneous lake with a complex bathymetry. The governing equations are the unsteady Reynolds-Averaged Navier–Stokes equations in which the non-hydrostatic pressure distribution has been included. The model was tested against analytical solutions and laboratory data for wind-induced currents and then applied to Esthwaite Water, a small lake in Cumbria, UK. The model was used to study the main model parameters and to generate typical circulation patterns for a variety of conditions in the lake. Simulations showed that a non-hydrostatic pressure distribution did not have any noticeable influence on the overall circulation pattern in the lake. However, comparisons with field data at some measurement stations in the near-shore region with sharply varying bottom topography showed that the hydrodynamic pressure component had some influence on the vertical velocity profile.     

Two-dimensional numerical modeling of dam-break flows over natural terrain using a central explicit scheme

A two-dimensional (2D) numerical model has been developed to solve shallow water equations for simulation of dam-break flows. The spatial derivatives are discretized using a well-balanced explicit central upwind conservative scheme. The scheme is Riemann solver free and guarantees the positivity of the flow depth over complex topography if the Courant number is kept less than 0.25. The time integration is performed by Euler’s scheme. The model is verified against analytical results for water surface elevation and discharge for three benchmark test cases. A good agreement between analytical solutions and computed results is observed. The property of well-balancing in still water over an uneven bottom is also confirmed. The model is then validated by simulating a laboratory experiment in which a dam break flow propagates over a triangular obstacle. The model performance was found to be satisfactory. A dam break laboratory experimental test case on a frictionless horizontal bottom is also simulated for 2D validation of the model, and good agreement between simulation and the experimental data is observed. The suitability of the proposed model for real life applications is demonstrated by simulating the Malpasset dam-break event, which occurred in 1959 in France. The computed arrival time of the flood wave front and the maximum flow depths at various observation points matched well with the measurements on a 1/400 scale physical model. The overall performance indicates that this model can be applied for simulation of dam-break waves in real life cases.     

Positivity-preserving high order well-balanced discontinuous Galerkin methods for the shallow water equations

Shallow water equations with a non-flat bottom topography have been widely used to model flows in rivers and coastal areas. An important difficulty arising in these simulations is the appearance of dry areas where no water is present, as standard numerical methods may fail in the presence of these areas. These equations also have still water steady state solutions in which the flux gradients are nonzero but exactly balanced by the source term. In this paper we propose a high order discontinuous Galerkin method which can maintain the still water steady state exactly, and at the same time can preserve the non-negativity of the water height without loss of mass conservation. A simple positivity-preserving limiter, valid under suitable CFL condition, will be introduced in one dimension and then extended to two dimensions with rectangular meshes. Numerical tests are performed to verify the positivity-preserving property, well-balanced property, high order accuracy, and good resolution for smooth and discontinuous solutions.     

Balancing the source terms in a SPH model for solving the shallow water equations

A shallow flow generally features complex hydrodynamics induced by complicated domain topography and geometry. A numerical scheme with well-balanced flux and source term gradients is therefore essential before a shallow flow model can be applied to simulate real-world problems. The issue of source term balancing has been exhaustively investigated in grid-based numerical approaches, e.g. discontinuous Galerkin finite element methods and finite volume Godunov-type methods. In recent years, a relatively new computational method, smooth particle hydrodynamics (SPH), has started to gain popularity in solving the shallow water equations (SWEs). However, the well-balanced problem has not been fully investigated and resolved in the context of SPH. This work aims to discuss the well-balanced problem caused by a standard SPH discretization to the SWEs with slope source terms and derive a corrected SPH algorithm that is able to preserve the solution of lake at rest. In order to enhance the shock capturing capability of the resulting SPH model, the Monotone Upwind-centered Scheme for Conservation Laws (MUSCL) is also explored and applied to enable Riemann solver based artificial viscosity. The new SPH model is validated against several idealized benchmark tests and a real-world dam-break case and promising results are obtained.     

Difference between the Lagrangian trajectories and Eulerian residual velocity fields in the southwestern Yellow Sea

The responses to tidal and/or wind forces of Lagrangian trajectories and Eulerian residual velocity in the southwestern Yellow Sea are investigated using a high-resolution circulation model. The simulated tidal harmonic constants agree well with observations and existing studies. The numerical experiment reproduces the long-range southeastward Eulerian residual current over the sloping bottom around the Yangtze Bank also shown in previous studies. However, the modeled drifters deployed at the northeastern flank of the Yangtze Bank in the simulation move northeastward, crossing over this strong southeastward Eulerian residual current rather than following it. Additional sensitivity experiments reveal that the influence of the Eulerian tidal residual currents on Lagrangian trajectories is relatively weaker than that of the wind driven currents. This result is consistent with the northeastward movement of ARGOS surface drifters actually released in the southwestern Yellow Sea. Further experiments suggest that the quadratic nature of the bottom friction is the crucial factor, in the southwestern Yellow Sea, for the weaker influence of the Eulerian tidal residual currents on the Lagrangian trajectories. This study demonstrates that the Lagrangian trajectories do not follow the Eulerian residual velocity fields in the shallow coastal regions of the southwestern Yellow Sea.     

Numerical resolution of well-balanced shallow water equations with complex source terms

This paper presents a well-balanced numerical scheme for simulating frictional shallow flows over complex domains involving wetting and drying. The proposed scheme solves, in a finite volume Godunov-type framework, a set of pre-balanced shallow water equations derived by considering pressure balancing. Non-negative reconstruction of Riemann states and compatible discretization of slope source term produce stable and well-balanced solutions to shallow flow hydrodynamics over complex topography. The friction source term is discretized using a splitting implicit scheme. Limiting value of the friction force is derived to ensure stability. This new numerical scheme is validated against four theoretical benchmark tests and then applied to reproduce a laboratory dam break over a domain with irregular bed profile.     

Tides and mixing in the northwestern East China Sea Part I: Rotating and reversing tidal flows

Bottom-mounted ADV and ADCP instruments in combination with CTD profiling measurements taken along the Chinese coast of the East China Sea were used to study the vertical structure of temperature, salinity, and velocity in reversing tidal currents on a shallow inner shelf and in rotating tidal flows over a deeper sloping bottom of the outer shelf. These two regimes of barotropic tide affect small-scale dynamics in the lower part of the water column differently. The reversing flow was superimposed by seiches of ∼2.3 h period generated in semienclosed Jiaozhou Bay located nearby. As the tidal vector rotates over the sloping bottom, the height of the near-bottom logarithmic layer is subjected to tidal-induced variations. A maximum of horizontal velocity U_max appears at the upper boundary of the log layer during the first half of the current vector rotation from the minor to the major axis of tidal ellipse. In rotating tidal flow, vertical shear generated at the seafloor, propagated slowly to the water interior up to the height of U_max, with a phase speed of ∼5 m/h. The time-shifted shear inside the water column, relative to the shear at the bottom, was associated with periodically changing increases and decreases of the tidal velocity above the log layer toward the sea surface. In reversing flows, the shear generated near the bottom and the shear at the upper levels were almost in phase.     

High-order finite volume WENO schemes for the shallow water equations with dry states

The shallow water equations are used to model flows in rivers and coastal areas, and have wide applications in ocean, hydraulic engineering, and atmospheric modeling. These equations have still water steady state solutions in which the flux gradients are balanced by the source term. It is desirable to develop numerical methods which preserve exactly these steady state solutions. Another main difficulty usually arising from the simulation of dam breaks and flood waves flows is the appearance of dry areas where no water is present. If no special attention is paid, standard numerical methods may fail near dry/wet front and produce non-physical negative water height. A high-order accurate finite volume weighted essentially non-oscillatory (WENO) scheme is proposed in this paper to address these difficulties and to provide an efficient and robust method for solving the shallow water equations. A simple, easy-to-implement positivity-preserving limiter is introduced. One- and two-dimensional numerical examples are provided to verify the positivity-preserving property, well-balanced property, high-order accuracy, and good resolution for smooth and discontinuous solutions.     

Two-dimensional numerical model of two-layer shallow water equations for confluence simulation

This study presents a finite-volume explicit method to solve 2D two-layer shallow water equations. This numerical model is intended to describe two-layer shallow flows in which the superposed layers differ in velocity, density and rheology in a two-dimensional domain. The rheological behavior of mudflow or debris flow is called the Bingham fluid. Thus, the shear stress on rigid bed can be derived from the constitutive equation. The computational approach adopts the HLL scheme, a novel approach for the purpose of computing a Godunov flux and solving the Riemann problem approximately proposed by Harten, Lax and van Leer, as a basic building block, treats the bottom slope by lateralizing the momentum flux, and refines the scheme using the Strang splitting to manage the frictional source term. This study successfully performed 2D two-layer shallow water computations on a rigid bed. The proposed numerical model can describe the variety of depths and velocities of substances including water and mud, when the hyperconcentrated tributary flows into the main river. The analytical results in this study will be valuable for further advanced research and for designing or planning hydraulic engineering structures.     

Space discontinuous Galerkin method for shallow water flows—kinetic and HLLC flux,and potential vorticity generation

In this paper, a second order space discontinuous Galerkin (DG) method is presented for the numerical solution of inviscid shallow water flows over varying bottom topography. Novel in the implementation is the use of HLLC and kinetic numerical fluxes¹ in combination with a dissipation operator, applied only locally around discontinuities to limit spurious numerical oscillations. Numerical solutions over (non-)uniform meshes are verified against exact solutions; the numerical error in the L₂-norm and the convergence of the solution are computed. Bore–vortex interactions are studied analytically and numerically to validate the model; these include bores as “breaking waves” in a channel and a bore traveling over a conical and Gaussian hump. In these complex numerical test cases, we correctly predict the generation of potential vorticity by non-uniform bores. Finally, we successfully validate the numerical model against measurements of steady oblique hydraulic jumps in a channel with a contraction. In the latter case, the kinetic flux is shown to be more robust.     

Vertical structure of tidal current in a typically coastal raft-culture area

In this paper vertical structure of tidal current in a typically coastal raft-culture area is discussed by field measurement and a numerical model. The observations show that the vertical structure changed dramatically. A tidal surface boundary layer (SBL) is well formed due to the frictional effects induced by extensive, high-density suspended culture as surface obstruction. Both the aquaculture drag and the bottom friction are much higher than those in non-raft-culture areas, and show an obvious variation with tidal flow. The significant earlier ebbing and earlier flooding appear in the upper water column instead of the seabed. And the maximal phase lag is about 1 h within one tide cycle. A 1D hydrodynamic model was modified to include the SBL and parameterized with the field data. It replicated the observed velocity profile and was then used to investigate the impacts of varying culture density and bottom friction on the vertical tidal-current structure. Modeling results indicate that the surface current velocity was largely damped because culture activities enhanced the frictional effects on flow intensively. The magnitude and vertical structure of tidal current are determined together with aquaculture drag and bottom friction. In addition, the vertical velocity structure has a nonlinear trend along with culture density and bottom friction. This study is a theoretical foundation for optimizing aquaculture configuration through regulating culture density and species distribution.     

Numerical calculation for bed variation in compound-meandering channel using depth integrated model without assumption of shallow water flow

In a compound meandering channel, patterns of flow structures and bed variations change with increasing water depth owing to complex momentum exchange between high-velocity flow in a main channel and low-velocity flows in flood plains. We have developed a new quasi-three-dimensional model without the shallow water assumption, i.e., hydrostatic pressure distribution; our method is known as the general bottom velocity computation (BVC) method. In this method, a set of depth-integrated equations, including depth-integrated momentum and vorticity equations, are prepared for evaluating bottom velocity and vertical velocity distributions. The objective of this study is to develop a bed variation calculation method for both single and compound meandering channels by using the BVC method coupled with a sediment transport model. This paper shows that the BVC method can reproduce the pattern change of bed variation in a compound meandering channel flow with increasing relative depth. The variation in sediment transport rate due to overbank flow is explained by experimental and computational results.     

Numerical simulation of a dam break for an actual river terrain environment

A two‐dimensional (2D) finite‐difference shallow water model based on a second‐order hybrid type of total variation diminishing (TVD) approximate solver with a MUSCL limiter function was developed to model flooding and inundation problems where the evolution of the drying and wetting interface is numerically challenging. Both a minimum positive depth (MPD) scheme and a non‐MPD scheme were employed to handle the advancement of drying and wetting fronts. We used several model problems to verify the model, including a dam break in a slope channel, a dam break flooding over a triangular obstacle, an idealized circular dam‐break, and a tide flow over a mound. Computed results agreed well with the experiment data and other numerical results available. The model was then applied to simulate the dam breaking and flooding of Hsindien Creek, Taiwan, with the detailed river basin topography. Computed flooding scenarios show reasonable flow characteristics. Though the average speed of flooding is 6–7 m s^?1, which corresponds to the subcritical flow condition (Fr < 1), the local maximum speed of flooding is 14·12 m s^?1, which corresponds to the supercritical flow condition (Fr ≈ 1·31). It is necessary to conduct some kind of comparison of the numerical results with measurements/experiments in further studies. Nevertheless, the model exhibits its capability to capture the essential features of dam‐break flows with drying and wetting fronts. It also exhibits the potential to provide the basis for computationally efficient flood routing and warning information. Copyright © 2006 John Wiley & Sons, Ltd.     

Effect of stratification on the frequency of bounded Rossby modes over a non-flat bottom

This work attempts to express and analyze the challenges, induced by stratification, affecting the Rossby-topographic eigenmodes of a closed domain with a general uneven bottom of arbitrary shape filled with a uniform fluid in the unperturbed configuration. The modified eigenmodes have been computed analytically: stratification is introduced in the mathematical form of a perturbation of a homogeneous fluid over a non-flat bottom. The eigenmodes lose their barotropic character and differences appear in the dynamical fields (velocity and pressure) from upper to lower layer, as expected. Expressions for the baroclinic and ageostrophic velocity components due to the perturbation are given. The analysis is carried out in the frame of linear shallow water approximation. All terms have been retained apart from nonlinear advection in the governing equations. We find that the frequencies of the eigenmodes change; an analytical expression of frequency correction as a function of layer density difference and interface depth is found. Initial results for some elementary geometrical settings with a waveguide bottom are determined and expressed in a concise, easily readable closed form. The results obtained in the shallow water approximation are expanded in series with respect to the Rossby number. Next, they are compared with the frequency correction obtained in an alternative framework in which the quasi-geostrophic approximation is used, and a purely baroclinic perturbation is imposed from the outset as the result of the introduction of stratification in the otherwise homogeneous fluid. In this scenario, reduced gravity and the ratio of upper to lower layer depth are, in turn, used as the expansion parameters in lieu of the Rossby number.