Numerical solutions for 2D depth-averaged shallow water equations

2018 ◽  
Vol 13 (2) ◽  
pp. 79-90 ◽  
Author(s):  
Huda Altaie ◽  
Pierre Dreyfuss
2011 ◽  
Vol 675 ◽  
pp. 1-26 ◽  
Author(s):  
R. ROTUNNO ◽  
J. B. KLEMP ◽  
G. H. BRYAN ◽  
D. J. MURAKI

Nearly all analytical models of lock-exchange flow are based on the shallow-water approximation. Since the latter approximation fails at the leading edges of the mutually intruding fluids of lock-exchange flow, solutions to the shallow-water equations can be obtained only through the specification of front conditions. In the present paper, analytic solutions to the shallow-water equations for non-Boussinesq lock-exchange flow are given for front conditions deriving from free-boundary arguments. Analytic solutions are also derived for other proposed front conditions – conditions which appear to the shallow-water system as forced boundary conditions. Both solutions to the shallow-water equations are compared with the numerical solutions of the Navier–Stokes equations and a mixture of successes and failures is recorded. The apparent success of some aspects of the forced solutions of the shallow-water equations, together with the fact that in a real fluid the density interface is a free boundary, shows the need for an improved theory of lock-exchange flow taking into account non-hydrostatic effects for density interfaces intersecting rigid boundaries.


Symmetry ◽  
2020 ◽  
Vol 12 (3) ◽  
pp. 345
Author(s):  
Sudi Mungkasi ◽  
Stephen Gwyn Roberts

This paper proposes some formulations of weak local residuals of shallow-water-type equations, namely, one-, one-and-a-half-, and two-dimensional shallow water equations. Smooth parts of numerical solutions have small absolute values of weak local residuals. Rougher parts of numerical solutions have larger absolute values of weak local residuals. This behaviour enables the weak local residuals to detect parts of numerical solutions which are smooth and rough (non-smooth). Weak local residuals that we formulate are implemented successfully as refinement or coarsening indicators for adaptive mesh finite volume methods used to solve shallow water equations.


2013 ◽  
Vol 714 ◽  
pp. 634-643 ◽  
Author(s):  
Chris Garrett ◽  
Patrick Cummins

AbstractThe maximum power that can be obtained from a confined array of turbines in steady or tidal flows is considered using the two-dimensional shallow-water equations and representing the turbine farm by a uniform local increase in friction within a circle. Analytical results supported by dimensional reasoning and numerical solutions show that the maximum power depends on the dominant term in the momentum equation for flows perturbed on the scale of the farm. If friction dominates in the basic flow, the maximum power is a fraction (half for linear friction and 0.75 for quadratic friction) of the dissipation within the circle in the undisturbed state; if the advective terms dominate, the maximum power is a fraction of the undisturbed kinetic energy flux into the front of the turbine farm; if the acceleration dominates, the maximum power is similar to that for the linear frictional case, but with the friction coefficient replaced by twice the tidal frequency.


Author(s):  
Setyo Nugroho ◽  
Mohamad Riyadi

ABSTRACT In this paper , the propagation of waves with an initial condition on the foundation are varied( varying bottom ) simulated by the model Linear Shallow Water Equations ( LSWE ) 1D . 1D LSWE solution with an initial condition was approached with numerical solutions , which use spectral methods . To reduce waves so as not to repeat back to the spatial domain need to be defined damping zone. The simulation results showed that the more superficial level it will increase the amplitude of the wave . Keywords : Wave Propagation , LSWE ID , Varying Bottom , MetodeSspectral , Damping Zone


Author(s):  
Kamel Al-Khaled ◽  
M. K. Al-Safeen

In this paper, the homotopy perturbation method is adopted to find explicit and numerical solutions for systems of non-linear fractional shallow water equations. The fractional derivatives are described in the Caputo sense. We apply both the homotopy perturbation method and the homotopy analysis method, to solve  certain shallow water equations with time-fractional derivatives, and explicitly construct convergent power series solutions. The  results obtained reveal that these  methods are  both very effective and simple for finding approximate solutions. Some numerical examples and plots are presented to illustrate the efficiency and reliability of these methods.  


2011 ◽  
Vol 2011 ◽  
pp. 1-17 ◽  
Author(s):  
Montri Maleewong

A modified predictor-corrector scheme combining with the depth gradient method (DGM) and the weighted average flux (WAF) method has been presented to solve the one-dimensional shallow water equations with source terms. Approximate solutions in the predictor step are obtained by the DGM with piecewise-linear reconstructions in each cell volume. The source terms can then be calculated directly by these predicted values at the corresponding half-time step. In the corrector step, the TVD version of the WAF method is applied to calculate the numerical fluxes at the same half-time step for each cell face. The accuracy of numerical solutions is shown by applying the method to solve various test cases in both steady and unsteady problems with and without source terms. It shows that the numerical results are in good agreement with the existing analytical solutions as well as experimental data in some test cases.


2017 ◽  
Vol 817 ◽  
pp. 406-438 ◽  
Author(s):  
Jim Thomas ◽  
K. Shafer Smith ◽  
Oliver Bühler

We investigate theoretically and numerically the modulation of near-inertial waves by a larger-amplitude geostrophically balanced mean flow. Because the excited wave is initially trapped in the mixed layer, it projects onto a broad spectrum of vertical modes, each mode $n$ being characterized by a Burger number, $Bu_{n}$, proportional to the square of the vertical scale of the mode. Using numerical simulations of the hydrostatic Boussinesq equations linearized about a prescribed balanced background flow, we show that the evolution of the wave field depends strongly on the spectrum of $Bu_{n}$ relative to the Rossby number of the balanced flow, $\unicode[STIX]{x1D716}$, with smaller relative $Bu_{n}$ leading to smaller horizontal scales in the wave field, faster accumulation of wave amplitude in anticyclones and faster propagation of wave energy into the deep ocean. This varied behaviour of the wave may be understood by considering the dynamics in each mode separately; projecting the linearized hydrostatic Boussinesq equations onto modes yields a set of linear shallow water equations, with $Bu_{n}$ playing the role of the reduced gravity. The wave modes fall into two asymptotic regimes, defined by the scalings $Bu_{n}\sim O(1)$ for low modes and $Bu_{n}\sim O(\unicode[STIX]{x1D716})$ for high modes. An amplitude equation derived for the former regime shows that vertical propagation is weak for low modes. The high-mode regime is the basis of the Young & Ben Jelloul (J. Mar. Res., vol. 55, 1997, pp. 735–766) theory. This theory is here extended to $O(\unicode[STIX]{x1D716}^{2})$, from which amplitude equations for the subregimes $Bu_{n}\sim O(\unicode[STIX]{x1D716}^{1/2})$ and $Bu_{n}\sim O(\unicode[STIX]{x1D716}^{2})$ are derived. The accuracy of each approximation is demonstrated by comparing numerical solutions of the respective amplitude equation to simulations of the linearized shallow water equations in the same regime. We emphasize that since inertial wave energy and shear are distributed across vertical modes, their overall modulation is due to the collective behaviour of the wave field in each regime. A unified treatment of these regimes is a novel feature of this work.


Sign in / Sign up

Export Citation Format

Share Document