This paper proposes a new version of the high-resolution entropy-consistent(EC-Limited)flux for hyperbolic conservation laws based on a new minmod-type slope limiter.Firstly,we identify the numerical entropy productio...This paper proposes a new version of the high-resolution entropy-consistent(EC-Limited)flux for hyperbolic conservation laws based on a new minmod-type slope limiter.Firstly,we identify the numerical entropy production,a third-order differential term deduced from the previous work of Ismail and Roe[11].The corresponding dissipation term is added to the original Roe flux to achieve entropy consistency.The new,resultant entropy-consistent(EC)flux has a general and explicit analytical form without any corrective factor,making it easy to compute and a less-expensive method.The inequality constraints are imposed on the standard piece-wise quadratic reconstruction to enforce the pointwise values of bounded-type numerical solutions.We design the new minmod slope limiter as combining two separate limiters for left and right states.We propose the EC-Limited flux by adding this reconstruction data method to the primitive variables rather than to the conservative variables of the EC flux to preserve the equilibrium of the primitive variables.These resulting fluxes are easily applied to general hyperbolic conservation laws while having attractive features:entropy-stable,robust,and non-oscillatory.To illustrate the potential of these proposed fluxes,we show the applications to the Burgers equation and the Euler equations.展开更多
1 Introduction China has a vast area of continental shelf and is very rich in marine resources,but because of the complex geological environment and frequent geological disasters,the utilization of marine resources an...1 Introduction China has a vast area of continental shelf and is very rich in marine resources,but because of the complex geological environment and frequent geological disasters,the utilization of marine resources and the construction of marine engineering are limited(Zhu et al.,2016).As the展开更多
A two-dimensional (2D) dam-break flow numerical model was developed based on the finite-volume total variation diminishing (TVD) and monotone upstream-centered scheme for conservation laws (MUSCL)-Hancock scheme...A two-dimensional (2D) dam-break flow numerical model was developed based on the finite-volume total variation diminishing (TVD) and monotone upstream-centered scheme for conservation laws (MUSCL)-Hancock scheme, which has second-order accuracy in both time and space. A Harten-Lax-van Leer-contact (HLLC) approximate Riemann solver was used to evaluate fluxes. The TVD MUSCL-Hancock numerical scheme utilizes slope limiters, such as the minmod, double minmod, superbee, van Albada, and van Leer limiters, to prevent spurious oscillations and maintain monotonicity near discontinuities. A comparative study of the impact of various slope limiters on the accuracy of the numerical flow model was conducted with several dam-break examples including wet and dry bed cases. The numerical results of the superbee and double minmod limiters agree better with the theoretical solution and have higher accuracy than other limiters in one-dimensional (1D) space. The ratio of the downstream water depth to the upstream water depth was used to select the proper slope limiter. For the 2D numerical model, the superbee limiter should not be used, owing to significant numerical dispersion.展开更多
One-dimensional open channel flows are simulated using the discontinuous Galerkin finite element method. Three different explicit time marching schemes, including multistep/multistage schemes, are evaluated for differ...One-dimensional open channel flows are simulated using the discontinuous Galerkin finite element method. Three different explicit time marching schemes, including multistep/multistage schemes, are evaluated for different channel shapes for accuracy and efficiency. The Forward Euler, second-order Adam-Bashforth (multistep), and second-order total variation diminishing (TVD) Runge-Kutta (multistage) time marching schemes are utilized. The role of monotonized central, minmod, and zero TVD slope limiters for each of the time marching scheme is investigated. The numerical flux is approximated using HLL function. The accuracy and robustness of different time marching schemes are evaluated for steady and unsteady flows using analytical and measured data. The unsteady flows include dam break tests with wet and dry beds downstream of the dam in prismatic (rectangular, trapezoidal, triangular, and parabolic cross-sections) and non-prismatic (natural river) channels. The steady flow test involves simulation of hydraulic jump in a diverging rectangular channel. The various schemes are evaluated by comparing accuracy using statistical measures and efficiency using maximum possible time step size as well as CPU runtime. The second-order Adam-Bashforth time marching scheme is found to have the best accuracy and efficiency among the time stepping schemes tested.展开更多
Groundwater flows play a key role in the recharge of aquifers, the transport of solutes through subsurface systems or the control of surface runoff. Predicting these processes requires the use of groundwater models wi...Groundwater flows play a key role in the recharge of aquifers, the transport of solutes through subsurface systems or the control of surface runoff. Predicting these processes requires the use of groundwater models with their applicability directly linked to their accuracy and computational efficiency. In this paper, we present a new method to model water dynamics in variably- saturated porous media. Our model is based on a fully-explicit discontinuous-Galerkin formulation of the 3D Richards equation, which shows a perfect scaling on parallel architectures. We make use of an adapted jump penalty term for the discontinuous-Galerkin scheme and of a slope limiter algorithm to produce oscillation-free exactly conservative solutions. We show that such an approach is particularly well suited to infiltration fronts. The model results are in good agreement with the reference model Hydrus-lD and seem promising for large scale applications involving a coarse representation of saturated soil.展开更多
We discuss the development,verification,and performance of a GPU accelerated discontinuous Galerkin method for the solutions of two dimensional nonlinear shallow water equations.The shallow water equations are hyperbo...We discuss the development,verification,and performance of a GPU accelerated discontinuous Galerkin method for the solutions of two dimensional nonlinear shallow water equations.The shallow water equations are hyperbolic partial differential equations and are widely used in the simulation of tsunami wave propagations.Our algorithms are tailored to take advantage of the single instruction multiple data(SIMD)architecture of graphic processing units.The time integration is accelerated by local time stepping based on a multi-rate Adams-Bashforth scheme.A total variational bounded limiter is adopted for nonlinear stability of the numerical scheme.This limiter is coupled with a mass and momentum conserving positivity preserving limiter for the special treatment of a dry or partially wet element in the triangulation.Accuracy,robustness and performance are demonstrated with the aid of test cases.Furthermore,we developed a unified multi-threading model OCCA.The kernels expressed in OCCA model can be cross-compiled with multi-threading models OpenCL,CUDA,and OpenMP.We compare the performance of the OCCA kernels when cross-compiled with these models.展开更多
基金the National Natural Science Found Project of China through project number 11971075.
文摘This paper proposes a new version of the high-resolution entropy-consistent(EC-Limited)flux for hyperbolic conservation laws based on a new minmod-type slope limiter.Firstly,we identify the numerical entropy production,a third-order differential term deduced from the previous work of Ismail and Roe[11].The corresponding dissipation term is added to the original Roe flux to achieve entropy consistency.The new,resultant entropy-consistent(EC)flux has a general and explicit analytical form without any corrective factor,making it easy to compute and a less-expensive method.The inequality constraints are imposed on the standard piece-wise quadratic reconstruction to enforce the pointwise values of bounded-type numerical solutions.We design the new minmod slope limiter as combining two separate limiters for left and right states.We propose the EC-Limited flux by adding this reconstruction data method to the primitive variables rather than to the conservative variables of the EC flux to preserve the equilibrium of the primitive variables.These resulting fluxes are easily applied to general hyperbolic conservation laws while having attractive features:entropy-stable,robust,and non-oscillatory.To illustrate the potential of these proposed fluxes,we show the applications to the Burgers equation and the Euler equations.
基金supported by NSFC Open Research Cruise (Cruise No. NORC2015-05 and Cruise No. NORC2015-06)funded by Shiptime Sharing Project of NSFC
文摘1 Introduction China has a vast area of continental shelf and is very rich in marine resources,but because of the complex geological environment and frequent geological disasters,the utilization of marine resources and the construction of marine engineering are limited(Zhu et al.,2016).As the
基金supported by the National Natural Science Foundation of China(Grants No.51679170,51379157,and 51439007)
文摘A two-dimensional (2D) dam-break flow numerical model was developed based on the finite-volume total variation diminishing (TVD) and monotone upstream-centered scheme for conservation laws (MUSCL)-Hancock scheme, which has second-order accuracy in both time and space. A Harten-Lax-van Leer-contact (HLLC) approximate Riemann solver was used to evaluate fluxes. The TVD MUSCL-Hancock numerical scheme utilizes slope limiters, such as the minmod, double minmod, superbee, van Albada, and van Leer limiters, to prevent spurious oscillations and maintain monotonicity near discontinuities. A comparative study of the impact of various slope limiters on the accuracy of the numerical flow model was conducted with several dam-break examples including wet and dry bed cases. The numerical results of the superbee and double minmod limiters agree better with the theoretical solution and have higher accuracy than other limiters in one-dimensional (1D) space. The ratio of the downstream water depth to the upstream water depth was used to select the proper slope limiter. For the 2D numerical model, the superbee limiter should not be used, owing to significant numerical dispersion.
文摘One-dimensional open channel flows are simulated using the discontinuous Galerkin finite element method. Three different explicit time marching schemes, including multistep/multistage schemes, are evaluated for different channel shapes for accuracy and efficiency. The Forward Euler, second-order Adam-Bashforth (multistep), and second-order total variation diminishing (TVD) Runge-Kutta (multistage) time marching schemes are utilized. The role of monotonized central, minmod, and zero TVD slope limiters for each of the time marching scheme is investigated. The numerical flux is approximated using HLL function. The accuracy and robustness of different time marching schemes are evaluated for steady and unsteady flows using analytical and measured data. The unsteady flows include dam break tests with wet and dry beds downstream of the dam in prismatic (rectangular, trapezoidal, triangular, and parabolic cross-sections) and non-prismatic (natural river) channels. The steady flow test involves simulation of hydraulic jump in a diverging rectangular channel. The various schemes are evaluated by comparing accuracy using statistical measures and efficiency using maximum possible time step size as well as CPU runtime. The second-order Adam-Bashforth time marching scheme is found to have the best accuracy and efficiency among the time stepping schemes tested.
基金funded by the Fond de la Recherche Scientifique de Belgique (FRSFNRS)
文摘Groundwater flows play a key role in the recharge of aquifers, the transport of solutes through subsurface systems or the control of surface runoff. Predicting these processes requires the use of groundwater models with their applicability directly linked to their accuracy and computational efficiency. In this paper, we present a new method to model water dynamics in variably- saturated porous media. Our model is based on a fully-explicit discontinuous-Galerkin formulation of the 3D Richards equation, which shows a perfect scaling on parallel architectures. We make use of an adapted jump penalty term for the discontinuous-Galerkin scheme and of a slope limiter algorithm to produce oscillation-free exactly conservative solutions. We show that such an approach is particularly well suited to infiltration fronts. The model results are in good agreement with the reference model Hydrus-lD and seem promising for large scale applications involving a coarse representation of saturated soil.
基金The authors gratefully acknowledge travel grants from Pan-American Advanced Studies Institute,grant from DOE and ANL(ANL Subcontract No.1F-32301 on DOE grant No.DE-AC02-06CH11357)grant from ONR(Award No.N00014-13-1-0873)fellowships from Ken Kennedy Institute of technology at Rice University and support from Shell(Shell Agreement No.PT22584),NVIDIA,and AMD.The authors also acknowledge Dr.Jesse Chan for fruitful discussions during the preparation of this manuscript.
文摘We discuss the development,verification,and performance of a GPU accelerated discontinuous Galerkin method for the solutions of two dimensional nonlinear shallow water equations.The shallow water equations are hyperbolic partial differential equations and are widely used in the simulation of tsunami wave propagations.Our algorithms are tailored to take advantage of the single instruction multiple data(SIMD)architecture of graphic processing units.The time integration is accelerated by local time stepping based on a multi-rate Adams-Bashforth scheme.A total variational bounded limiter is adopted for nonlinear stability of the numerical scheme.This limiter is coupled with a mass and momentum conserving positivity preserving limiter for the special treatment of a dry or partially wet element in the triangulation.Accuracy,robustness and performance are demonstrated with the aid of test cases.Furthermore,we developed a unified multi-threading model OCCA.The kernels expressed in OCCA model can be cross-compiled with multi-threading models OpenCL,CUDA,and OpenMP.We compare the performance of the OCCA kernels when cross-compiled with these models.