scholarly journals A mass conservative TR-BDF2 semi-implicit semi-Lagrangian DG discretization of the shallow water equations on general structured meshes of quadrilaterals

2016 ◽  
Vol 7 (3) ◽  
pp. 165-190 ◽  
Author(s):  
Giovanni Tumolo
Keyword(s):  
Shallow Water ◽  
Shallow Water Equations ◽  
Time Integration ◽  
Mass Conservation ◽  
Momentum Equation ◽  
Implicit Time ◽  
Time Step ◽  
Mesh Distortion ◽  
Structured Meshes ◽  
Conservative Discretization

Abstract As an extension of a previous work considering a fully advective formulation on Cartesian meshes, a mass conservative discretization approach is presented here for the shallow water equations, based on discontinuous finite elements on general structured meshes of quadrilaterals. A semi-implicit time integration is performed by employing the TR-BDF2 scheme and is combined with the semi-Lagrangian technique for the momentum equation only. Indeed, in order to simplify the derivation of the discrete linear Helmoltz equation to be solved at each time-step, a non-conservative formulation of the momentum equation is employed. The Eulerian flux form is considered instead for the continuity equation in order to ensure mass conservation. Numerical results show that on distorted meshes and for relatively high polynomial degrees, the proposed numerical method fully conserves mass and presents a higher level of accuracy than a standard off-centered Crank Nicolson approach. This is achieved without any significant imprinting of the mesh distortion on the solution.

MAPPING TIDAL CURRENTS AND RESIDUAL CURRENTS BY USE OF COASTAL ACOUSTIC TOMOGRAPHY

2017 ◽  
Author(s):  
Xiao-Hua Zhu ◽  
Xiao-Hua Zhu ◽  
Ze-Nan Zhu ◽  
Ze-Nan Zhu ◽  
Xinyu Guo ◽  
...  
Keyword(s):  
Shallow Water ◽  
Shallow Water Equations ◽  
Mean Amplitude ◽  
Momentum Equation ◽  
Tidal Currents ◽  
Acoustic Tomography ◽  
Mean Values ◽  
Residual Currents ◽  
Water Depths ◽  
Deep Area

A coastal acoustic tomography (CAT) experiment for mapping the tidal currents in the Zhitouyang Bay was successfully carried out with seven acoustic stations during July 12 to 13, 2009. The horizontal distributions of tidal current in the tomography domain are calculated by the inverse analysis in which the travel time differences for sound traveling reciprocally are used as data. Spatial mean amplitude ratios M2 : M4 : M6 are 1.00 : 0.15 : 0.11. The shallow-water equations are used to analyze the generation mechanisms of M4 and M6. In the deep area, velocity amplitudes of M4 measured by CAT agree well with those of M4 predicted by the advection terms in the shallow water equations, indicating that M4 in the deep area where water depths are larger than 60 m is predominantly generated by the advection terms. M6 measured by CAT and M6 predicted by the nonlinear quadratic bottom friction terms agree well in the area where water depths are less than 20 m, indicating that friction mechanisms are predominant for generating M6 in the shallow area. Dynamic analysis of the residual currents using the tidally averaged momentum equation shows that spatial mean values of the horizontal pressure gradient due to residual sea level and of the advection of residual currents together contribute about 75% of the spatial mean values of the advection by the tidal currents, indicating that residual currents in this bay are induced mainly by the nonlinear effects of tidal currents.

scholarly journals A Fully Implicit Finite Volume Scheme for a Seawater Intrusion Problem in Coastal Aquifers

Water ◽  
2020 ◽  
Vol 12 (6) ◽  
pp. 1639
Author(s):  
Abdelkrim Aharmouch ◽  
Brahim Amaziane ◽  
Mustapha El Ossmani ◽  
Khadija Talali
Keyword(s):  
Finite Volume ◽  
Coastal Aquifers ◽  
Nonlinear Partial Differential Equations ◽  
Time Integration ◽  
Mass Conservation ◽  
Coupled System ◽  
Finite Volume Scheme ◽  
Time Step ◽  
Volume Scheme ◽  
Fully Implicit

We present a numerical framework for efficiently simulating seawater flow in coastal aquifers using a finite volume method. The mathematical model consists of coupled and nonlinear partial differential equations. Difficulties arise from the nonlinear structure of the system and the complexity of natural fields, which results in complex aquifer geometries and heterogeneity in the hydraulic parameters. When numerically solving such a model, due to the mentioned feature, attempts to explicitly perform the time integration result in an excessively restricted stability condition on time step. An implicit method, which calculates the flow dynamics at each time step, is needed to overcome the stability problem of the time integration and mass conservation. A fully implicit finite volume scheme is developed to discretize the coupled system that allows the use of much longer time steps than explicit schemes. We have developed and implemented this scheme in a new module in the context of the open source platform DuMu X . The accuracy and effectiveness of this new module are demonstrated through numerical investigation for simulating the displacement of the sharp interface between saltwater and freshwater in groundwater flow. Lastly, numerical results of a realistic test case are presented to prove the efficiency and the performance of the method.

Semi-Lagrangian exponential time-integration method for the shallow water equations on the cubed sphere grid

2020 ◽  
Vol 35 (6) ◽  
pp. 355-366
Author(s):  
Vladimir V. Shashkin ◽  
Gordey S. Goyman
Keyword(s):  
Shallow Water ◽  
Gravity Waves ◽  
Shallow Water Equations ◽  
Time Integration ◽  
Standard Test ◽  
Integration Scheme ◽  
Test Problems ◽  
Cubed Sphere ◽  
Lagrangian Scheme ◽  
Inertia Gravity Waves

AbstractThis paper proposes the combination of matrix exponential method with the semi-Lagrangian approach for the time integration of shallow water equations on the sphere. The second order accuracy of the developed scheme is shown. Exponential semi-Lagrangian scheme in the combination with spatial approximation on the cubed-sphere grid is verified using the standard test problems for shallow water models. The developed scheme is as good as the conventional semi-implicit semi-Lagrangian scheme in accuracy of slowly varying flow component reproduction and significantly better in the reproduction of the fast inertia-gravity waves. The accuracy of inertia-gravity waves reproduction is close to that of the explicit time-integration scheme. The computational efficiency of the proposed exponential semi-Lagrangian scheme is somewhat lower than the efficiency of semi-implicit semi-Lagrangian scheme, but significantly higher than the efficiency of explicit, semi-implicit, and exponential Eulerian schemes.

scholarly journals Investigation of Air Pocket Behavior in Pipelines Using Rigid Column Model and Contributions of Time Integration Schemes

Water ◽  
2021 ◽  
Vol 13 (6) ◽  
pp. 785
Author(s):  
Arman Rokhzadi ◽  
Musandji Fuamba
Keyword(s):  
Transient Flow ◽  
Time Integration ◽  
Driving Pressure ◽  
Implicit Time ◽  
Numerical Dissipation ◽  
Integration Scheme ◽  
Time Step ◽  
Column Model ◽  
Integration Schemes ◽  
Backward Euler

This paper studies the air pressurization problem caused by a partially pressurized transient flow in a reservoir-pipe system. The purpose of this study is to analyze the performance of the rigid column model in predicting the attenuation of the air pressure distribution. In this regard, an analytic formula for the amplitude and frequency will be derived, in which the influential parameters, particularly, the driving pressure and the air and water lengths, on the damping can be seen. The direct effect of the driving pressure and inverse effect of the product of the air and water lengths on the damping will be numerically examined. In addition, these numerical observations will be examined by solving different test cases and by comparing to available experimental data to show that the rigid column model is able to predict the damping. However, due to simplified assumptions associated with the rigid column model, the energy dissipation, as well as the damping, is underestimated. In this regard, using the backward Euler implicit time integration scheme, instead of the classical fourth order explicit Runge–Kutta scheme, will be proposed so that the numerical dissipation of the backward Euler implicit scheme represents the physical dissipation. In addition, a formula will be derived to calculate the appropriate time step size, by which the dissipation of the heat transfer can be compensated.

scholarly journals A Semi-Implicit Version of the MPAS-Atmosphere Dynamical Core

2015 ◽  
Vol 143 (9) ◽  
pp. 3838-3855 ◽  
Author(s):  
Steven Sandbach ◽  
John Thuburn ◽  
Danail Vassilev ◽  
Michael G. Duda
Keyword(s):  
Time Integration ◽  
Standard Test ◽  
Newton Iteration ◽  
Atmospheric Modeling ◽  
Implicit Time ◽  
Integration Scheme ◽  
Baroclinic Wave ◽  
Time Step ◽  
Integration Schemes ◽  
Helmholtz Problem

Abstract An important question for atmospheric modeling is the viability of semi-implicit time integration schemes on massively parallel computing architectures. Semi-implicit schemes can provide increased stability and accuracy. However, they require the solution of an elliptic problem at each time step, creating concerns about their parallel efficiency and scalability. Here, a semi-implicit (SI) version of the Model for Prediction Across Scales (MPAS) is developed and compared with the original model version, which uses a split Runge–Kutta (SRK3) time integration scheme. The SI scheme is based on a quasi-Newton iteration toward a Crank–Nicolson scheme. Each Newton iteration requires the solution of a Helmholtz problem; here, the Helmholtz problem is derived, and its solution using a geometric multigrid method is described. On two standard test cases, a midlatitude baroclinic wave and a small-planet nonhydrostatic gravity wave, the SI and SRK3 versions produce almost identical results. On the baroclinic wave test, the SI version can use somewhat larger time steps (about 60%) than the SRK3 version before losing stability. The SI version costs 10%–20% more per step than the SRK3 version, and the weak and strong scalability characteristics of the two versions are very similar for the processor configurations the authors have been able to test (up to 1920 processors). Because of the spatial discretization of the pressure gradient in the lowest model layer, the SI version becomes unstable in the presence of realistic orography. Some further work will be needed to demonstrate the viability of the SI scheme in this case.

scholarly journals WAVETRISK-2.1: an adaptive dynamical core for ocean modelling

2021 ◽  
Author(s):  
Nicholas Keville-Reynolds Kevlahan ◽  
Florian Lemarié
Keyword(s):  
Free Surface ◽  
Shallow Water ◽  
Time Integration ◽  
Ocean Model ◽  
Mesh Adaptivity ◽  
Time Step ◽  
Vertical Coordinate ◽  
Slip Boundary ◽  
Modelling Framework ◽  
Atmosphere Model

Abstract. This paper introduces WAVETRISK-2.1 (i.e. WAVETRISK-OCEAN), an incompressible version of the atmosphere model wavetrisk-1.x with free-surface. This new model is built on the same wavelet-based dynamically adaptive core as wavetrisk, which itself uses DYNANICO's mimetic vector-invariant multilayer rotating shallow water formulation. Both codes use a Lagrangian vertical coordinate with conservative remapping. The ocean variant solves the incompressible multilayer shallow water equations with inhomogeneous density layers. Time integration uses barotropic--baroclinic mode splitting via an semi-implicit free surface formulation, which is about 34–44 times faster than an unsplit explicit time-stepping. The barotropic and baroclinic estimates of the free surface are reconciled at each time step using layer dilation. No slip boundary conditions at coastlines are approximated using volume penalization. The vertical eddy viscosity and diffusivity coefficients are computed from a closure model based on turbulent kinetic energy (TKE). Results are presented for a standard set of ocean model test cases adapted to the sphere (seamount, upwelling and baroclinic turbulence). An innovative feature of wavetrisk-ocean is that it could be coupled easily to the wavetrisk atmosphere model, thus providing a first building block toward an integrated Earth-system model using a consistent modelling framework with dynamic mesh adaptivity and mimetic properties.

Sign in / Sign up

Export Citation Format

Share Document