scholarly journals Well-balanced finite volume multi-resolution schemes for solving the Ripa models

2021 ◽  
Vol 13 (3) ◽  
pp. 168781402110034
Author(s):  
Asad Rehman ◽  
Ishtiaq Ali ◽  
Saqib Zia ◽  
Shamsul Qamar
Keyword(s):  
Steady State ◽  
Finite Volume ◽  
Numerical Test ◽  
Temperature Gradients ◽  
Water Model ◽  
Test Problems ◽  
Numerical Schemes ◽  
Weno Schemes ◽  
Horizontal Temperature Gradients ◽  
Fifth Order

In this article, fifth order well-balanced finite volume multi-resolution weighted essentially non-oscillatory (FV MR-WENO) schemes are constructed for solving one-dimensional and two-dimensional Ripa models. The Ripa system generalizes the shallow water model by incorporating horizontal temperature gradients. The presence of temperature gradients and source terms in the Ripa models introduce difficulties in developing high order accurate numerical schemes which can preserve exactly the steady-state conditions. The proposed numerical methods are capable to exactly preserve the steady-state solutions and maintain non-oscillatory property near the shock transitions. Moreover, in the procedure of derivation of the FV MR-WENO schemes unequal central spatial stencils are used and linear weights can be chosen any positive numbers with only restriction that their total sum is one. Various numerical test problems are considered to check the validity and accuracy of the derived numerical schemes. Further, the results obtained from considered numerical schemes are compared with those of a high resolution central upwind scheme and available exact solutions of the Ripa model.

High Order Fixed-Point Sweeping WENO Methods for Steady State of Hyperbolic Conservation Laws and Its Convergence Study

2016 ◽  
Vol 20 (4) ◽  
pp. 835-869 ◽  
Author(s):  
Liang Wu ◽  
Yong-Tao Zhang ◽  
Shuhai Zhang ◽  
Chi-Wang Shu
Keyword(s):  
Fixed Point ◽  
Steady State ◽  
Conservation Laws ◽  
Hyperbolic Conservation Laws ◽  
High Order ◽  
Efficient Computation ◽  
Weno Schemes ◽  
Hyperbolic Conservation ◽  
Time Marching ◽  
Fifth Order

AbstractFixed-point iterative sweeping methods were developed in the literature to efficiently solve static Hamilton-Jacobi equations. This class of methods utilizes the Gauss-Seidel iterations and alternating sweeping strategy to achieve fast convergence rate. They take advantage of the properties of hyperbolic partial differential equations (PDEs) and try to cover a family of characteristics of the corresponding Hamilton-Jacobi equation in a certain direction simultaneously in each sweeping order. Different from other fast sweeping methods, fixed-point iterative sweeping methods have the advantages such as that they have explicit forms and do not involve inverse operation of nonlinear local systems. In principle, it can be applied in solving very general equations using any monotone numerical fluxes and high order approximations easily. In this paper, based on the recently developed fifth order WENO schemes which improve the convergence of the classical WENO schemes by removing slight post-shock oscillations, we design fifth order fixed-point sweeping WENO methods for efficient computation of steady state solution of hyperbolic conservation laws. Especially, we show that although the methods do not have linear computational complexity, they converge to steady state solutions much faster than regular time-marching approach by stability improvement for high order schemes with a forward Euler time-marching.

scholarly journals A finite volume method for two-moment cosmic-ray hydrodynamics on a moving mesh

2021 ◽  
Author(s):  
T Thomas ◽  
C Pfrommer ◽  
R Pakmor
Keyword(s):  
Finite Volume Method ◽  
Finite Volume ◽  
Source Term ◽  
Cosmic Ray ◽  
Moving Mesh ◽  
Test Problems ◽  
Integration Step ◽  
Custom Made ◽  
Anisotropic Transport ◽  
Volume Method

Abstract We present a new numerical algorithm to solve the recently derived equations of two-moment cosmic ray hydrodynamics (CRHD). The algorithm is implemented as a module in the moving mesh Arepo code. Therein, the anisotropic transport of cosmic rays (CRs) along magnetic field lines is discretised using a path-conservative finite volume method on the unstructured time-dependent Voronoi mesh of Arepo. The interaction of CRs and gyroresonant Alfvén waves is described by short-timescale source terms in the CRHD equations. We employ a custom-made semi-implicit adaptive time stepping source term integrator to accurately integrate this interaction on the small light-crossing time of the anisotropic transport step. Both the transport and the source term integration step are separated from the evolution of the magneto-hydrodynamical equations using an operator split approach. The new algorithm is tested with a variety of test problems, including shock tubes, a perpendicular magnetised discontinuity, the hydrodynamic response to a CR overpressure, CR acceleration of a warm cloud, and a CR blast wave, which demonstrate that the coupling between CR and magneto-hydrodynamics is robust and accurate. We demonstrate the numerical convergence of the presented scheme using new linear and non-linear analytic solutions.

scholarly journals A Finite-Volume Icosahedral Shallow-Water Model on a Local Coordinate

2009 ◽  
Vol 137 (4) ◽  
pp. 1422-1437 ◽  
Author(s):  
Jin-Luen Lee ◽  
Alexander E. MacDonald
Keyword(s):  
Shallow Water ◽  
Finite Volume ◽  
Spherical Surface ◽  
Second Order ◽  
Water Model ◽  
Numerical Dissipation ◽  
Shallow Water Model ◽  
Finite Volume Model ◽  
Volume Model ◽  
Cartesian Coordinate

Abstract An icosahedral-hexagonal shallow-water model (SWM) on the sphere is formulated on a local Cartesian coordinate based on the general stereographic projection plane. It is discretized with the third-order Adam–Bashforth time-differencing scheme and the second-order finite-volume operators for spatial derivative terms. The finite-volume operators are applied to the model variables defined on the nonstaggered grid with the edge variables interpolated using polynomial interpolation. The projected local coordinate reduces the solution space from the three-dimensional, curved, spherical surface to the two-dimensional plane and thus reduces the number of complete sets of basis functions in the Vandermonde matrix, which is the essential component of the interpolation. The use of a local Cartesian coordinate also greatly simplifies the mathematic formulation of the finite-volume operators and leads to the finite-volume integration along straight lines on the plane, rather than along curved lines on the spherical surface. The SWM is evaluated with the standard test cases of Williamson et al. Numerical results show that the icosahedral SWM is free from Pole problems. The SWM is a second-order finite-volume model as shown by the truncation error convergence test. The lee-wave numerical solutions are compared and found to be very similar to the solutions shown in other SWMs. The SWM is stably integrated for several weeks without numerical dissipation using the wavenumber 4 Rossby–Haurwitz solution as an initial condition. It is also shown that the icosahedral SWM achieves mass conservation within round-off errors as one would expect from a finite-volume model.

Chemical equilibrium systems as numerical test problems

1990 ◽  
Vol 16 (2) ◽  
pp. 143-151 ◽  
Author(s):  
Keith Meintjes ◽  
Alexander P. Morgan
Keyword(s):  
Chemical Equilibrium ◽  
Numerical Test ◽  
Test Problems

A Quasi-three-dimensional Finite-volume Shallow Water Model for Green Water on Deck

2013 ◽  
Vol 57 (03) ◽  
pp. 125-140
Author(s):  
Daniel A. Liut ◽  
Kenneth M. Weems ◽  
Tin-Guen Yen
Keyword(s):  
Shallow Water ◽  
Finite Volume ◽  
Time Domain ◽  
Degrees Of Freedom ◽  
Initial Conditions ◽  
Three Dimensional ◽  
Green Water ◽  
Water Model ◽  
Ship Motions ◽  
Shallow Water Model

A quasi-three-dimensional hydrodynamic model is presented to simulate shallow water phenomena. The method is based on a finite-volume approach designed to solve shallow water equations in the time domain. The nonlinearities of the governing equations are considered. The methodology can be used to compute green water effects on a variety of platforms with six-degrees-of-freedom motions. Different boundary and initial conditions can be applied for multiple types of moving platforms, like a ship's deck, tanks, etc. Comparisons with experimental data are discussed. The shallow water model has been integrated with the Large Amplitude Motions Program to compute the effects of green water flow over decks within a time-domain simulation of ship motions in waves. Results associated to this implementation are presented.

scholarly journals A Well-balanced Finite Volume Scheme for Shallow Water Equations with Porosity: Application to Modelling Flow through Rigid Vegetation

2018 ◽  
Vol 40 ◽  
pp. 05032
Author(s):  
Minh H. Le ◽  
Virgile Dubos ◽  
Marina Oukacine ◽  
Nicole Goutal
Keyword(s):  
Shallow Water ◽  
Finite Volume ◽  
Shallow Water Equations ◽  
Simple Wave ◽  
Water Model ◽  
Strong Interactions ◽  
Finite Volume Scheme ◽  
Volume Scheme ◽  
Rigid Vegetation ◽  
Vertical Cylinders

Strong interactions exist between flow dynamics and vegetation in open-channel. Depth-averaged shallow water equations can be used for such a study. However, explicit representation of vegetation can lead to very high resolution of the mesh since the vegetation is often modelled as vertical cylinders. Our work aims to study the ability of a single porosity-based shallow water model for these applications. More attention on flux and source terms discretizations are required in order to archive the well-balancing and shock capturing properties. We present a new Godunov-type finite volume scheme based on a simple-wave approximation and compare it with some other methods in the literature. A first application with experimental data was performed.

Sign in / Sign up

Export Citation Format

Share Document