scholarly journals The Finite Volume WENO with Lax–Wendroff Scheme for Nonlinear System of Euler Equations

Mathematics ◽  
2018 ◽  
Vol 6 (10) ◽  
pp. 211 ◽  
Author(s):  
Haoyu Dong ◽  
Changna Lu ◽  
Hongwei Yang
Keyword(s):  
Finite Volume ◽  
Time Discretization ◽  
Euler System ◽  
Weno Scheme ◽  
Local Characteristic ◽  
Finite Volume Scheme ◽  
Mesh Structure ◽  
Hyperbolic Conservation Law ◽  
Higher Derivative ◽  
Accuracy Order

We develop a Lax–Wendroff scheme on time discretization procedure for finite volume weighted essentially non-oscillatory schemes, which is used to simulate hyperbolic conservation law. We put more focus on the implementation of one-dimensional and two-dimensional nonlinear systems of Euler functions. The scheme can keep avoiding the local characteristic decompositions for higher derivative terms in Taylor expansion, even omit partly procedure of the nonlinear weights. Extensive simulations are performed, which show that the fifth order finite volume WENO (Weighted Essentially Non-oscillatory) schemes based on Lax–Wendroff-type time discretization provide a higher accuracy order, non-oscillatory properties and more cost efficiency than WENO scheme based on Runge–Kutta time discretization for certain problems. Those conclusions almost agree with that of finite difference WENO schemes based on Lax–Wendroff time discretization for Euler system, while finite volume scheme has more flexible mesh structure, especially for unstructured meshes.

scholarly journals The Simple Finite Volume Lax-Wendroff Weighted Essentially Nonoscillatory Schemes for Shallow Water Equations with Bottom Topography

2018 ◽  
Vol 2018 ◽  
pp. 1-15 ◽  
Author(s):  
Changna Lu ◽  
Luoyan Xie ◽  
Hongwei Yang
Keyword(s):  
Finite Difference ◽  
Shallow Water ◽  
Finite Volume ◽  
Shallow Water Equations ◽  
Bottom Topography ◽  
Time Discretization ◽  
Weno Scheme ◽  
Weighted Essentially Nonoscillatory ◽  
Accuracy Order ◽  
Essentially Nonoscillatory

A Lax-Wendroff-type procedure with the high order finite volume simple weighted essentially nonoscillatory (SWENO) scheme is proposed to simulate the one-dimensional (1D) and two-dimensional (2D) shallow water equations with topography influence in source terms. The system of shallow water equations is discretized using the simple WENO scheme in space and Lax-Wendroff scheme in time. The idea of Lax-Wendroff time discretization can avoid part of characteristic decomposition and calculation of nonlinear weights. The type of simple WENO was first developed by Zhu and Qiu in 2016, which is more simple than classical WENO fashion. In order to maintain good, high resolution and nonoscillation for both continuous and discontinuous flow and suit problems with discontinuous bottom topography, we use the same idea of SWENO reconstruction for flux to treat the source term in prebalanced shallow water equations. A range of numerical examples are performed; as a result, comparing with classical WENO reconstruction and Runge-Kutta time discretization, the simple Lax-Wendroff WENO schemes can obtain the same accuracy order and escape nonphysical oscillation adjacent strong shock, while bringing less absolute truncation error and costing less CPU time for most problems. These conclusions agree with that of finite difference Lax-Wendroff WENO scheme for shallow water equations, while finite volume method has more flexible mesh structure compared to finite difference method.

scholarly journals Construction of a low Mach finite volume scheme for the isentropic Euler system with porosity

2021 ◽  
Author(s):  
Pascal Omnes ◽  
Stéphane Dellacherie ◽  
Jonathan Jung
Keyword(s):  
Mach Number ◽  
Finite Volume ◽  
Euler System ◽  
Finite Volume Scheme ◽  
Linear Wave ◽  
Finite Volume Schemes ◽  
Linear Wave Equation ◽  
Continuous Space ◽  
Low Mach Number ◽  
Non Linear System

Classical finite volume schemes for the Euler system  are not accurate at low Mach number and some fixes have to be used and were developed in a vast literature over the last two decades. The question we are interested in in this article is: What about if the porosity is no longer uniform? We first show that this problem may be understood on the linear wave equation taking into account porosity. We explain the influence of the cell geometry on the accuracy property at low Mach number. In the triangular case, the stationary space of the Godunov scheme approaches well enough the continuous space of constant pressure and divergence-free velocity, while this is not the case in the Cartesian case. On Cartesian meshes, a fix is proposed and accuracy at low Mach number is proved to be recovered. Based on the linear study, a numerical scheme and a low Mach fix for the non-linear system, with a non-conservative source term due to the porosity variations, is proposed and tested.

scholarly journals A finite volume scheme for the Euler system inspired by the two velocities approach

2019 ◽  
Vol 144 (1) ◽  
pp. 89-132 ◽  
Author(s):  
Eduard Feireisl ◽  
Mária Lukáčová-Medvid’ová ◽  
Hana Mizerová
Keyword(s):  
Finite Volume ◽  
Euler System ◽  
Finite Volume Scheme ◽  
Volume Scheme

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.

scholarly journals A Hydrodynamic-Based Robust Numerical Model for Debris Hazard and Risk Assessment

2021 ◽  
Vol 13 (14) ◽  
pp. 7955
Author(s):  
Yongde Kang ◽  
Jingming Hou ◽  
Yu Tong ◽  
Baoshan Shi
Keyword(s):  
Numerical Model ◽  
Debris Flow ◽  
Finite Volume ◽  
Graphics Processing Unit ◽  
Morphological Changes ◽  
Finite Volume Scheme ◽  
Processing Unit ◽  
Algebraic Equations ◽  
Flow Process ◽  
Computing Technique

Debris flow simulations are important in practical engineering. In this study, a graphics processing unit (GPU)-based numerical model that couples hydrodynamic and morphological processes was developed to simulate debris flow, transport, and morphological changes. To accurately predict the debris flow sediment transport and sediment scouring processes, a GPU-based parallel computing technique was used to accelerate the calculation. This model was created in the framework of a Godunov-type finite volume scheme and discretized into algebraic equations by the finite volume method. The mass and momentum fluxes were computed using the Harten, Lax, and van Leer Contact (HLLC) approximate Riemann solver, and the friction source terms were calculated using the proposed splitting point-implicit method. These values were evaluated using a novel 2D edge-based MUSCL scheme. The code was programmed using C++ and CUDA, which can run on GPUs to substantially accelerate the computation. After verification, the model was applied to the simulation of the debris flow process of an idealized example. The results of the new scheme better reflect the characteristics of the discontinuity of its movement and the actual law of the evolution of erosion and deposition over time. The research results provide guidance and a reference for the in-depth study of debris flow processes and disaster prevention and mitigation.

Sign in / Sign up

Export Citation Format

Share Document