SIMULATION OF A SPILLED OIL SLICK WITH A SHALLOW WATER MODEL WITH FREE BOUNDARY

2007 ◽  
Vol 17 (03) ◽  
pp. 393-409 ◽  
Author(s):  
B. DI MARTINO ◽  
P. ORENGA ◽  
M. PEYBERNES
Keyword(s):  
Free Boundary ◽  
Shallow Water ◽  
Lower Layer ◽  
Water Model ◽  
Shallow Water Model ◽  
New Approach ◽  
Ale Formulation ◽  
Special Basis ◽  
Spilled Oil ◽  
Discretized Problem

In this paper we present a new approach to describe the behaviour of a pollutant slick at the sea surface. To this end, we consider that the pollutant and the water are immiscible and we propose a two-layer model where the lower layer corresponds to the water and the upper layer represents the pollutant. Since the dimension of the pollutant slick is generally much smaller than the domain occupied by the sea, we propose to compute the motion of the pollutant with a shallow water model with free boundary only in the domain occupied by the pollutant. To discretize in time the problem with free boundary, we use an ALE formulation coupled with the characteristic method. Then, to solve the space discretized problem, we approximate the pollutant velocity by using a Galerkin method with a special basis which verifies the boundary conditions and simplifies significantly the resolution. Finally we test this work in a real situation: the dam of Calacuccia (Corsica).

Instabilities of coupled density fronts and their nonlinear evolution in the two-layer rotating shallow-water model: influence of the lower layer and of the topography

2013 ◽  
Vol 716 ◽  
pp. 528-565 ◽  
Author(s):  
Bruno Ribstein ◽  
Vladimir Zeitlin
Keyword(s):  
Shallow Water ◽  
Nonlinear Evolution ◽  
Lower Layer ◽  
Phase Locking ◽  
Water Model ◽  
Shallow Water Model ◽  
Wave Instability ◽  
Rotating Shallow Water ◽  
New Generation ◽  
Rotating Shallow Water Equations

AbstractWe undertake a detailed analysis of linear stability of geostrophically balanced double density fronts in the framework of the two-layer rotating shallow-water model on the $f$-plane with topography, the latter being represented by an escarpment beneath the fronts. We use the pseudospectral collocation method to identify and quantify different kinds of instabilities resulting from phase locking and resonances of frontal, Rossby, Poincaré and topographic waves. A swap in the leading long-wave instability from the classical barotropic form, resulting from the resonance of two frontal waves, to a baroclinic form, resulting from the resonance of Rossby and frontal waves, takes place with decreasing depth of the lower layer. Nonlinear development and saturation of these instabilities, and of an instability of topographic origin, resulting from the resonance of frontal and topographic waves, are studied and compared with the help of a new-generation well-balanced finite-volume code for multilayer rotating shallow-water equations. The results of the saturation for different instabilities are shown to produce very different secondary coherent structures. The influence of the topography on these processes is highlighted.

scholarly journals Two Different Aperiodic Phases of Wind-Driven Ocean Circulation in a Double-Gyre, Two-Layer Shallow-Water Model

2006 ◽  
Vol 36 (7) ◽  
pp. 1265-1286 ◽  
Author(s):  
Tomonori Matsuura ◽  
Mitsutaka Fujita
Keyword(s):  
Power Spectrum ◽  
Shallow Water ◽  
Metastable State ◽  
Ocean Circulation ◽  
Lower Layer ◽  
Water Model ◽  
Barotropic Instability ◽  
Shallow Water Model ◽  
Gyre Circulation ◽  
Double Gyre

Abstract A two-layer shallow-water model is used to investigate the transition of wind-driven double-gyre circulation from laminar flow to turbulence as the Reynolds number (Re) is systematically increased. Two distinctly different phases of turbulent double-gyre patterns and energy trajectories are exhibited before and after at Re = 95: deterministic and fully developed turbulent circulations. In the former phase, the inertial subgyres vary between an asymmetric solution and an antisymmetric solution and the double-gyre circulations reach the aperiodic solution mainly due to their barotropic instability. An integrated kinetic energy in the lower layer is slight and the generated mesoscale eddies are confined in the upper layer. The power spectrum of energies integrated over the whole domain at Re = 70 has peaks at the interannual periods (4–7 yr) and the interdecadal period (10–20 yr). The loops of the attractors take on one cycle at those periods and display the blue-sky catastrophe. At Re = 95, the double-gyre circulation reaches a metastable state and the attracters obtained from the three energies form a topological manifold. In the latter, as Re increases, the double-gyre varies from a metastable state to a chaotic state because of the barotropic instability of the eastward jet and the baroclinic instability of recirculation retrograde flow, and the eastward jet meanders significantly with interdecadal variability. The generated eddies cascade to the red side of the power spectrum as expected in the geostrophic turbulence. The main results in the simulation may indicate essential mechanisms for the appearance of multiple states of the Kuroshio and for low-frequency variations in the midlatitude ocean.

scholarly journals Effects of Moisture in a Two-Layer Model of the Midlatitude Jet Stream

2019 ◽  
Vol 77 (1) ◽  
pp. 131-147
Author(s):  
Eric Bembenek ◽  
David N. Straub ◽  
Timothy M. Merlis
Keyword(s):  
Kinetic Energy ◽  
Shallow Water ◽  
Lower Layer ◽  
Eddy Kinetic Energy ◽  
Water Model ◽  
Energy Response ◽  
Shallow Water Model ◽  
Similar Reduction ◽  
Model Flow ◽  
Two Layer Model

Abstract The effects of moisture on the energetics of a statistically stationary, baroclinically unstable jet representing the midlatitude atmosphere are examined using a two-layer, β-plane shallow-water model. Flow is driven by a relaxation of the interface between the two layers to a baroclinically unstable profile. Moisture is input to the lower layer by evaporation. When supersaturation occurs, precipitation is triggered and the related latent heat release drives a mass transfer between the two layers. A comparison between dry and moist reference atmospheres shows that precipitation reduces eddy kinetic energy. This is related to the meridional distribution of precipitation, which occurs on the poleward side of the jet (where the interface field is raised). This latitudinal structure of precipitation is related to a correlation between poleward flow and ascent, which is analyzed using a shallow-water analog to the ω equation. The precipitation effect on the energy budget is predominately due to zonal- and time-averaged terms. Because of this, dry simulations in which the thermal forcing is modified to mimic the effect of zonally averaged precipitation are carried out and compared with their precipitating counterparts. These simulations show a similar reduction of baroclinic eddy kinetic energy; however, the barotropic eddy kinetic energy response shows a larger difference.

The instability of vertically curved fronts in a two-layer shallow water model

2020 ◽  
Vol 32 (12) ◽  
pp. 124117
Author(s):  
M. W. Harris ◽  
F. J. Poulin ◽  
K. G. Lamb
Keyword(s):  
Shallow Water ◽  
Water Model ◽  
Shallow Water Model

scholarly journals Modelling Weirs in Two-Dimensional Shallow Water Models

Water ◽  
2021 ◽  
Vol 13 (16) ◽  
pp. 2152
Author(s):  
Gonzalo García-Alén ◽  
Olalla García-Fonte ◽  
Luis Cea ◽  
Luís Pena ◽  
Jerónimo Puertas
Keyword(s):  
Experimental Data ◽  
Shallow Water ◽  
Shallow Water Equations ◽  
Water Model ◽  
Two Dimensional ◽  
Shallow Water Model ◽  
Experimental Dataset ◽  
River Hydraulics ◽  
Shallow Water Models ◽  
Water Models

2D models based on the shallow water equations are widely used in river hydraulics. However, these models can present deficiencies in those cases in which their intrinsic hypotheses are not fulfilled. One of these cases is in the presence of weirs. In this work we present an experimental dataset including 194 experiments in nine different weirs. The experimental data are compared to the numerical results obtained with a 2D shallow water model in order to quantify the discrepancies that exist due to the non-fulfillment of the hydrostatic pressure hypotheses. The experimental dataset presented can be used for the validation of other modelling approaches.

scholarly journals Development of a Practical Model for Predicting Soil Salinity in a Salt Marsh in the Arakawa River Estuary

Water ◽  
2021 ◽  
Vol 13 (15) ◽  
pp. 2054
Author(s):  
Naoki Kuroda ◽  
Katsuhide Yokoyama ◽  
Tadaharu Ishikawa
Keyword(s):  
Salt Marsh ◽  
Hydraulic Conductivity ◽  
Shallow Water ◽  
Soil Salinity ◽  
River Estuary ◽  
Flow Simulation ◽  
Design Tool ◽  
Water Model ◽  
Shallow Water Model ◽  
Practical Model

Our group has studied the spatiotemporal variation of soil and water salinity in an artificial salt marsh along the Arakawa River estuary and developed a practical model for predicting soil salinity. The salinity of the salt marsh and the water level of a nearby channel were measured once a month for 13 consecutive months. The vertical profile of the soil salinity in the salt marsh was measured once monthly over the same period. A numerical flow simulation adopting the shallow water model faithfully reproduced the salinity variation in the salt marsh. Further, we developed a soil salinity model to estimate the soil salinity in a salt marsh in Arakawa River. The vertical distribution of the soil salinity in the salt marsh was uniform and changed at almost the same time. The hydraulic conductivity of the soil, moreover, was high. The uniform distribution of salinity and high hydraulic conductivity could be explained by the vertical and horizontal transport of salinity through channels burrowed in the soil by organisms. By combining the shallow water model and the soil salinity model, the soil salinity of the salt marsh was well reproduced. The above results suggest that a stable brackish ecotone can be created in an artificial salt marsh using our numerical model as a design tool.

Diffusion Experiments with a Global Discontinuous Galerkin Shallow-Water Model

2009 ◽  
Vol 137 (10) ◽  
pp. 3339-3350 ◽  
Author(s):  
Ramachandran D. Nair
Keyword(s):  
Shallow Water ◽  
Discontinuous Galerkin ◽  
Water Model ◽  
Order System ◽  
Small Scale ◽  
Shallow Water Model ◽  
First Order ◽  
Advection Diffusion Equation ◽  
Diffusion Scheme ◽  
Cubed Sphere

Abstract A second-order diffusion scheme is developed for the discontinuous Galerkin (DG) global shallow-water model. The shallow-water equations are discretized on the cubed sphere tiled with quadrilateral elements relying on a nonorthogonal curvilinear coordinate system. In the viscous shallow-water model the diffusion terms (viscous fluxes) are approximated with two different approaches: 1) the element-wise localized discretization without considering the interelement contributions and 2) the discretization based on the local discontinuous Galerkin (LDG) method. In the LDG formulation the advection–diffusion equation is solved as a first-order system. All of the curvature terms resulting from the cubed-sphere geometry are incorporated into the first-order system. The effectiveness of each diffusion scheme is studied using the standard shallow-water test cases. The approach of element-wise localized discretization of the diffusion term is easy to implement but found to be less effective, and with relatively high diffusion coefficients, it can adversely affect the solution. The shallow-water tests show that the LDG scheme converges monotonically and that the rate of convergence is dependent on the coefficient of diffusion. Also the LDG scheme successfully eliminates small-scale noise, and the simulated results are smooth and comparable to the reference solution.

Sign in / Sign up

Export Citation Format

Share Document