scholarly journals Physics-Capturing Discretizations for Spectral Wind-Wave Models

Fluids ◽  
2021 ◽  
Vol 6 (2) ◽  
pp. 52
Author(s):  
Marcel Zijlema

This paper discusses the discretization methods that have been commonly employed to solve the wave action balance equation, and that have gained a renewed interest with the widespread use of unstructured grids for third-generation spectral wind-wave models. These methods are the first-order upwind finite difference and first-order vertex-centered upwind finite volume schemes for the transport of wave action in geographical space. The discussion addresses the derivation of these schemes from a different perspective. A mathematical framework for mimetic discretizations based on discrete calculus is utilized herein. A key feature of this algebraic approach is that the process of exact discretization is segregated from the process of interpolation, the latter typically involved in constitutive relations. This can help gain insight into the performance characteristics of the discretization method. On this basis, we conclude that the upwind finite difference scheme captures the wave action flux conservation exactly, which is a plus for wave shoaling. In addition, we provide a justification for the intrinsic low accuracy of the vertex-centred upwind finite volume scheme, due to the physically inaccurate but common flux constitutive relation, and we propose an improvement to overcome this drawback. Finally, by way of a comparative demonstration, a few test cases is introduced to establish the ability of the considered methods to capture the relevant physics on unstructured triangular meshes.

2017 ◽  
Vol 14 (03) ◽  
pp. 415-454 ◽  
Author(s):  
Ujjwal Koley ◽  
Nils Henrik Risebro ◽  
Christoph Schwab ◽  
Franziska Weber

This paper proposes a finite difference multilevel Monte Carlo algorithm for degenerate parabolic convection–diffusion equations where the convective and diffusive fluxes are allowed to be random. We establish a notion of stochastic entropy solutions to these equations. Our chief goal is to efficiently compute approximations to statistical moments of these stochastic entropy solutions. To this end, we design a multilevel Monte Carlo method based on a finite volume scheme for each sample. We present a novel convergence rate analysis of the combined multilevel Monte Carlo finite volume method, allowing in particular for low [Formula: see text]-integrability of the random solution with [Formula: see text], and low deterministic convergence rates (here, the theoretical rate is [Formula: see text]). We analyze the design and error versus work of the multilevel estimators. We obtain that the maximal rate (based on optimizing possibly the pessimistic upper bounds on the discretization error) is obtained for [Formula: see text], for finite volume convergence rate of [Formula: see text]. We conclude with numerical experiments.


Fluids ◽  
2020 ◽  
Vol 5 (4) ◽  
pp. 199
Author(s):  
Yong G. Lai ◽  
Han Sang Kim

The near-shore and estuary environment is characterized by complex natural processes. A prominent feature is the wind-generated waves, which transfer energy and lead to various phenomena not observed where the hydrodynamics is dictated only by currents. Over the past several decades, numerical models have been developed to predict the wave and current state and their interactions. Most models, however, have relied on the two-model approach in which the wave model is developed independently of the current model and the two are coupled together through a separate steering module. In this study, a new wave model is developed and embedded in an existing two-dimensional (2D) depth-integrated current model, SRH-2D. The work leads to a new wave–current model based on the one-model approach. The physical processes of the new wave model are based on the latest third-generation formulation in which the spectral wave action balance equation is solved so that the spectrum shape is not pre-imposed and the non-linear effects are not parameterized. New contributions of the present study lie primarily in the numerical method adopted, which include: (a) a new operator-splitting method that allows an implicit solution of the wave action equation in the geographical space; (b) mixed finite volume and finite difference method; (c) unstructured polygonal mesh in the geographical space; and (d) a single mesh for both the wave and current models that paves the way for the use of the one-model approach. An advantage of the present model is that the propagation of waves from deep water to shallow water in near-shore and the interaction between waves and river inflows may be carried out seamlessly. Tedious interpolations and the so-called multi-model steering operation adopted by many existing models are avoided. As a result, the underlying interpolation errors and information loss due to matching between two meshes are avoided, leading to an increased computational efficiency and accuracy. The new wave model is developed and verified using a number of cases. The verified near-shore wave processes include wave shoaling, refraction, wave breaking and diffraction. The predicted model results compare well with the analytical solution or measured data for all cases.


2010 ◽  
Vol 20 (02) ◽  
pp. 265-295 ◽  
Author(s):  
JÉRÔME DRONIOU ◽  
ROBERT EYMARD ◽  
THIERRY GALLOUËT ◽  
RAPHAÈLE HERBIN

We investigate the connections between several recent methods for the discretization of anisotropic heterogeneous diffusion operators on general grids. We prove that the Mimetic Finite Difference scheme, the Hybrid Finite Volume scheme and the Mixed Finite Volume scheme are in fact identical up to some slight generalizations. As a consequence, some of the mathematical results obtained for each of the methods (such as convergence properties or error estimates) may be extended to the unified common framework. We then focus on the relationships between this unified method and nonconforming Finite Element schemes or Mixed Finite Element schemes. We also show that for isotropic operators, on particular meshes such as triangular meshes with acute angles, the unified method boils down to the well-known efficient two-point flux Finite Volume scheme.


2019 ◽  
Vol 59 (1) ◽  
pp. 55-73
Author(s):  
Anurak Busaman ◽  
Khamron Mekchay ◽  
Suchada Siripant ◽  
Somporn Chuai-Aree ◽  
Rhysa McNeil

Mathematics ◽  
2019 ◽  
Vol 7 (10) ◽  
pp. 923 ◽  
Author(s):  
Abdul Ghafoor ◽  
Sirajul Haq ◽  
Manzoor Hussain ◽  
Poom Kumam ◽  
Muhammad Asif Jan

In this paper, a wavelet based collocation method is formulated for an approximate solution of (1 + 1)- and (1 + 2)-dimensional time fractional diffusion wave equations. The main objective of this study is to combine the finite difference method with Haar wavelets. One and two dimensional Haar wavelets are used for the discretization of a spatial operator while time fractional derivative is approximated using second order finite difference and quadrature rule. The scheme has an excellent feature that converts a time fractional partial differential equation to a system of algebraic equations which can be solved easily. The suggested technique is applied to solve some test problems. The obtained results have been compared with existing results in the literature. Also, the accuracy of the scheme has been checked by computing L 2 and L ∞ error norms. Computations validate that the proposed method produces good results, which are comparable with exact solutions and those presented before.


Water ◽  
2020 ◽  
Vol 12 (6) ◽  
pp. 1639
Author(s):  
Abdelkrim Aharmouch ◽  
Brahim Amaziane ◽  
Mustapha El Ossmani ◽  
Khadija Talali

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.


Sign in / Sign up

Export Citation Format

Share Document