scholarly journals Zonostrophic Beta-plumes, Breaking Waves and Zonal Jets in Locally-forced, Large-Scale Shallow Water Experiments

2021 ◽  
Author(s):  
Peter Read ◽  
Roland Young ◽  
Hélène Scolan ◽  
Enrico Ferrero ◽  
Boris Galperin ◽  
...  
Keyword(s):  
Shallow Water ◽  
Large Scale ◽  
Breaking Waves ◽  
Zonal Jets

Waves, jets and vortices: Regime transitions in shallow water turbulence

2021 ◽  
Author(s):  
Nikos Bakas
Keyword(s):  
Shallow Water ◽  
Large Scale ◽  
Rossby Waves ◽  
Phase Speed ◽  
Critical Value ◽  
Turbulent Regime ◽  
Wave Regime ◽  
Frequency Spectra ◽  
Spectra Analysis ◽  
Zonal Jets

<p>Forced-dissipative beta-plane turbulence in a single-layer shallow-water fluid has been widely considered as a simplified model of planetary turbulence as it exhibits turbulence self-organization into large-scale structures such as robust zonal jets and strong vortices. In this study we perform a series of numerical simulations to analyze the characteristics of the emerging structures as a function of the planetary vorticity gradient and the deformation radius. We report four regimes that appear as the energy input rate ε of the random stirring that supports turbulence in the flow increases. A homogeneous turbulent regime for low values of ε, a regime in which large scale Rossby waves form abruptly when ε passes a critical value, a regime in which robust zonal jets coexist with weaker Rossby waves when ε passes a second critical value and a regime of strong materially coherent propagating vortices for large values of ε. The wave regime which is not predicted by standard cascade theories of turbulence anisotropization and the vortex regime are studied thoroughly. Wavenumber-frequency spectra analysis shows that the Rossby waves in the second regime remain phase coherent over long times. The coherent vortices are identified using the Lagrangian Averaged Deviation (LAVD) method. The statistics of the vortices (lifetime, radius, strength and speed) are reported as a function of the large scale parameters. We find that the strong vortices propagate zonally with a phase speed that is equal or larger than the long Rossby wave speed and advect the background turbulence leading to a non-dispersive line in the wavenumber-frequency spectra.</p>

Numerical Simulations of Forced Shallow-Water Turbulence: Effects of Moist Convection on the Large-Scale Circulation of Jupiter and Saturn

2007 ◽  
Vol 64 (9) ◽  
pp. 3132-3157 ◽  
Author(s):  
Adam P. Showman
Keyword(s):  
Numerical Simulations ◽  
Shallow Water ◽  
Shallow Water Equations ◽  
Large Scale ◽  
Spherical Geometry ◽  
Small Deformation ◽  
Moist Convection ◽  
Zonal Jets ◽  
Strong Control ◽  
Water Turbulence

Abstract To test the hypothesis that the zonal jets on Jupiter and Saturn result from energy injected by thunderstorms into the cloud layer, forced-dissipative numerical simulations of the shallow-water equations in spherical geometry are presented. The forcing consists of sporadic, isolated circular mass pulses intended to represent thunderstorms; the damping, representing radiation, removes mass evenly from the layer. These results show that the deformation radius provides strong control over the behavior. At deformation radii <2000 km (0.03 Jupiter radii), the simulations produce broad jets near the equator, but regions poleward of 15°–30° latitude instead become dominated by vortices. However, simulations at deformation radii >4000 km (0.06 Jupiter radii) become dominated by barotropically stable zonal jets with only weak vortices. The lack of midlatitude jets at a small deformation radii results from the suppression of the beta effect by column stretching; this effect has been previously documented in the quasigeostrophic system but never before in the full shallow-water system. In agreement with decaying shallow-water turbulence simulations, but in disagreement with Jupiter and Saturn, the equatorial flows in these forced simulations are always westward. In analogy with purely two-dimensional turbulence, the size of the coherent structures (jets and vortices) depends on the relative strengths of forcing and damping; stronger damping removes energy faster as it cascades upscale, leading to smaller vortices and more closely spaced jets in the equilibrated state. Forcing and damping parameters relevant to Jupiter produce flows with speeds up to 50–200 m s−1 and a predominance of anticyclones over cyclones, both in agreement with observations. However, the dominance of vortices over jets at deformation radii thought to be relevant to Jupiter (1000–3000 km) suggests that either the actual deformation radius is larger than previously believed or that three-dimensional effects, not included in the shallow-water equations, alter the dynamics in a fundamental manner.

On the vorticity transport due to dissipating or breaking waves in shallow-water flow

2000 ◽  
Vol 407 ◽  
pp. 235-263 ◽  
Author(s):  
OLIVER BÜHLER
Keyword(s):  
Numerical Simulations ◽  
Shallow Water ◽  
Gravity Waves ◽  
Large Scale ◽  
Bottom Topography ◽  
Breaking Waves ◽  
Small Scale ◽  
Practical Consideration ◽  
Mean Flow ◽  
Vorticity Transport

Theoretical and numerical results are presented on the transport of vorticity (or potential vorticity) due to dissipating gravity waves in a shallow-water system with background rotation and bottom topography. The results are obtained under the assumption that the flow can be decomposed into small-scale gravity waves and a large-scale mean flow. The particle-following formalism of ‘generalized Lagrangian-mean’ theory is then used to derive an ‘effective mean force’ that captures the vorticity transport due to the dissipating waves. This can be achieved without neglecting other, non-dissipative, effects which is an important practical consideration. It is then shown that the effective mean force obeys the so-called ‘pseudomomentum rule’, i.e. the force is approximately equal to minus the local dissipation rate of the wave's pseudomomentum. However, it is also shown that this holds only if the underlying dissipation mechanism is momentum-conserving. This requirement has important implications for numerical simulations, and these are discussed.The novelty of the results presented here is that they have been derived within a uniform theoretical framework, that they are not restricted to small wave amplitude, ray-tracing or JWKB-type approximations, and that they also include wave dissipation by breaking, or shock formation. The theory is tested carefully against shock-capturing nonlinear numerical simulations, which includes the detailed study of a wavetrain subject to slowly varying bottom topography. The theory is also cross-checked in the appropriate asymptotic limit against recently formulated weakly nonlinear theories. In addition to the general finite-amplitude theory, detailed small-amplitude expressions for the main results are provided in which the explicit appearance of Lagrangian fields can be avoided. The motivation for this work stems partly from an on-going study of high-altitude breaking of internal gravity waves in the atmosphere, and some preliminary remarks on atmospheric applications and on three-dimensional stratified versions of these results are given.

Towards Control of an Estuary

1987 ◽  
Vol 19 (9) ◽  
pp. 155-174
Author(s):  
Henk L. F. Saeijs
Keyword(s):  
Shallow Water ◽  
Storm Surge ◽  
Salt Marshes ◽  
Large Scale ◽  
Western Europe ◽  
Oxygen Depletion ◽  
Negative Effects ◽  
Intertidal Flats ◽  
Environmental Consequences ◽  
Short Time

The Delta Project is in its final stage. In 1974 it was subjected to political reconsideration, but it is scheduled now for completion in 1987. The final touches are being put to the storm-surge barrier and two compartment dams that divide the Oosterschelde into three areas: one tidal, one with reduced tide, and one a freshwater lake. Compartmentalization will result in 13% of channels, 45% of intertidal flats and 59% of salt marshes being lost. There is a net gain of 7% of shallow-water areas. Human interventions with large scale impacts are not new in the Oosterschelde but the large scale and short time in which these interventions are taking place are, as is the creation of a controlled tidal system. This article focusses on the area with reduced tide and compares resent day and expected characteristics. In this reduced tidal part salt marshes will extend by 30–70%; intertidal flats will erode to a lower level and at their edges, and the area of shallow water will increase by 47%. Biomass production on the intertidal flats will decrease, with consequences for crustaceans, fishes and birds. The maximum number of waders counted on one day and the number of ‘bird-days' will decrease drastically, with negative effects for the wader populations of western Europe. The net area with a hard substratum in the reduced tidal part has more than doubled. Channels will become shallower. Detritus import will not change significantly. Stratification and oxygen depletion will be rare and local. The operation of the storm-surge barrier and the closure strategy chosen are very important for the ecosystem. Two optional closure strategies can be followed without any additional environmental consequences. It was essential to determine a clearly defined plan of action for the whole area, and to make land-use choices from the outset. How this was done is briefly described.

Assessing the Numerical Accuracy of Complex Spherical Shallow-Water Flows

2007 ◽  
Vol 135 (11) ◽  
pp. 3876-3894 ◽  
Author(s):  
Ali R. Mohebalhojeh ◽  
David G. Dritschel
Keyword(s):  
Shallow Water ◽  
Fourth Order ◽  
Large Scale ◽  
Numerical Accuracy ◽  
Conservation Of Energy ◽  
Water Flows ◽  
Shallow Water Flows ◽  
Low Froude Number ◽  
Contour Advection ◽  
High Reynolds

Abstract The representation of nonlinear shallow-water flows poses severe challenges for numerical modeling. The use of contour advection with contour surgery for potential vorticity (PV) within the contour-advective semi-Lagrangian (CASL) algorithm makes it possible to handle near-discontinuous distributions of PV with an accuracy beyond what is accessible to conventional algorithms used in numerical weather and climate prediction. The emergence of complex distributions of the materially conserved quantity PV, in the absence of forcing and dissipation, results from large-scale shearing and deformation and is a common feature of high Reynolds number flows in the atmosphere and oceans away from boundary layers. The near-discontinuous PV in CASL sets a limit on the actual numerical accuracy of the Eulerian, grid-based part of CASL. For the spherical shallow-water equations, the limit is studied by comparing the accuracy of CASL algorithms with second-order-centered, fourth-order-compact, and sixth-order-supercompact finite differencing in latitude in conjunction with a spectral treatment in longitude. The comparison is carried out on an unstable midlatitude jet at order one Rossby number and low Froude number that evolves into complex vortical structures with sharp gradients of PV. Quantitative measures of global conservation of energy and angular momentum, and of imbalance as diagnosed using PV inversion by means of Bolin–Charney balance, indicate that fourth-order differencing attains the highest numerical accuracy achievable for such nonlinear, advectively dominated flows.

scholarly journals Turbulence of Capillary Waves on Shallow Water

Fluids ◽  
2021 ◽  
Vol 6 (5) ◽  
pp. 185
Author(s):  
Natalia Vladimirova ◽  
Ivan Vointsev ◽  
Alena Skoba ◽  
Gregory Falkovich
Keyword(s):  
Shallow Water ◽  
Large Scale ◽  
Capillary Waves ◽  
Total Momentum ◽  
Reflection Symmetry ◽  
Spontaneous Breakdown ◽  
Developed Turbulence ◽  
Degree Of Anisotropy ◽  
Zero Momentum ◽  
Square Box

We consider the developed turbulence of capillary waves on shallow water. Analytic theory shows that an isotropic cascade spectrum is unstable with respect to small angular perturbations, in particular, to spontaneous breakdown of the reflection symmetry and generation of nonzero momentum. By computer modeling we show that indeed a random pumping, generating on average zero momentum, produces turbulence with a nonzero total momentum. A strongly anisotropic large-scale pumping produces turbulence whose degree of anisotropy decreases along a cascade. It tends to saturation in the inertial interval and then further decreases in the dissipation interval. Surprisingly, neither the direction of the total momentum nor the direction of the compensated spectrum anisotropy is locked by our square box preferred directions (side or diagonal) but fluctuate.

Advanced CFD Simulations of free-surface flows around modern sailing yachts using a newly developed openFOAM solver

2016 ◽  
Author(s):  
Janek Meyer ◽  
Hannes Renzsch ◽  
Kai Graf ◽  
Thomas Slawig
Keyword(s):  
Free Surface ◽  
Large Scale ◽  
Breaking Waves ◽  
Free Surface Flows ◽  
Body Motion ◽  
Computational Time ◽  
Efficient Computation ◽  
Surface Flows ◽  
Scale Free ◽  
Wide Range

While plain vanilla OpenFOAM has strong capabilities with regards to quite a few typical CFD-tasks, some problems actually require additional bespoke solvers and numerics for efficient computation of high-quality results. One of the fields requiring these additions is the computation of large-scale free-surface flows as found e.g. in naval architecture. This holds especially for the flow around typical modern yacht hulls, often planing, sometimes with surface-piercing appendages. Particular challenges include, but are not limited to, breaking waves, sharpness of interface, numerical ventilation (aka streaking) and a wide range of flow phenomenon scales. A new OF-based application including newly implemented discretization schemes, gradient computation and rigid body motion computation is described. In the following the new code will be validated against published experimental data; the effect on accuracy, computational time and solver stability will be shown by comparison to standard OF-solvers (interFoam / interDyMFoam) and Star CCM+. The code’s capabilities to simulate complex “real-world” flows are shown on a well-known racing yacht design.

A Coupled-Mode Technique for the Run-Up of Non-Breaking Dispersive Waves on Plane Beaches

2006 ◽  
Author(s):  
K. A. Belibassakis ◽  
G. A. Athanassoulis
Keyword(s):  
Shallow Water ◽  
Water Waves ◽  
Wave Equations ◽  
Neumann Boundary ◽  
Breaking Waves ◽  
Numerical Results ◽  
Coupled Mode Theory ◽  
Mode Theory ◽  
Coupled Mode ◽  
Run Up

A coupled-mode model is developed and applied to the transformation and run-up of dispersive water waves on plane beaches. The present work is based on the consistent coupled-mode theory for the propagation of water waves in variable bathymetry regions, developed by Athanassoulis & Belibassakis (1999) and extended to 3D by Belibassakis et al (2001), which is suitably modified to apply to a uniform plane beach. The key feature of the coupled-mode theory is a complete modal-type expansion of the wave potential, containing both propagating and evanescent modes, being able to consistently satisfy the Neumann boundary condition on the sloping bottom. Thus, the present approach extends previous works based on the modified mild-slope equation in conjunction with analytical solution of the linearised shallow water equations, see, e.g., Massel & Pelinovsky (2001). Numerical results concerning non-breaking waves on plane beaches are presented and compared with exact analytical solutions; see, e.g., Wehausen & Laitone (1960, Sec. 18). Also, numerical results are presented concerning the run-up of non-breaking solitary waves on plane beaches and compared with the ones obtained by the solution of the shallow-water wave equations, Synolakis (1987), Li & Raichlen (2002), and experimental data, Synolakis (1987).

scholarly journals A shallow water analogue of asymmetric core-collapse, and neutron star kick/spin

2011 ◽  
Vol 7 (S279) ◽  
pp. 134-137
Author(s):  
Thierry Foglizzo ◽  
Frédéric Masset ◽  
Jérôme Guilet ◽  
Gilles Durand
Keyword(s):  
Neutron Star ◽  
Shallow Water ◽  
Acoustic Waves ◽  
Large Scale ◽  
Massive Stars ◽  
Core Collapse ◽  
Hydraulic Jumps ◽  
Core Collapse Supernova ◽  
Accretion Shock

AbstractMassive stars end their life with the gravitational collapse of their core and the formation of a neutron star. Their explosion as a supernova depends on the revival of a spherical accretion shock, located in the inner 200km and stalled during a few hundred milliseconds. Numerical simulations suggest that the large scale asymmetry of the neutrino-driven explosion is induced by a hydrodynamical instability named SASI. Its non radial character is able to influence the kick and the spin of the resulting neutron star. The SWASI experiment is a simple shallow water analog of SASI, where the role of acoustic waves and shocks is played by surface waves and hydraulic jumps. Distances in the experiment are scaled down by a factor one million, and time is slower by a factor one hundred. This experiment is designed to illustrate the asymmetric nature of core-collapse supernova.

Predicting mesh density for adaptive modelling of the global atmosphere

2009 ◽  
Vol 367 (1907) ◽  
pp. 4523-4542 ◽  
Author(s):  
Hilary Weller
Keyword(s):  
Shallow Water ◽  
Shallow Water Equations ◽  
Large Scale ◽  
Mesh Adaptation ◽  
Similar Proportion ◽  
Test Case ◽  
Uniform Mesh ◽  
Coarse Mesh ◽  
Mesh Density ◽  
The Cost

The shallow water equations are solved using a mesh of polygons on the sphere, which adapts infrequently to the predicted future solution. Infrequent mesh adaptation reduces the cost of adaptation and load-balancing and will thus allow for more accurate mapping on adaptation. We simulate the growth of a barotropically unstable jet adapting the mesh every 12 h. Using an adaptation criterion based largely on the gradient of the vorticity leads to a mesh with around 20 per cent of the cells of a uniform mesh that gives equivalent results. This is a similar proportion to previous studies of the same test case with mesh adaptation every 1–20 min. The prediction of the mesh density involves solving the shallow water equations on a coarse mesh in advance of the locally refined mesh in order to estimate where features requiring higher resolution will grow, decay or move to. The adaptation criterion consists of two parts: that resolved on the coarse mesh, and that which is not resolved and so is passively advected on the coarse mesh. This combination leads to a balance between resolving features controlled by the large-scale dynamics and maintaining fine-scale features.

Sign in / Sign up

Export Citation Format

Share Document