Abstract: Urban flooding resulting from a sudden release of
water due to dam-break or excessive rainfall is a serious threatening
environment hazard, which causes loss of human life and large
economic losses. Anticipating floods before they occur could
minimize human and economic losses through the implementation
of appropriate protection, provision, and rescue plans. This work
reports on the numerical modelling of flash flood propagation
in urban areas after an excessive rainfall event or dam-break.
A two-dimensional (2D) depth-averaged shallow water model is
used with a refined unstructured grid of triangles for representing
the urban area topography. The 2D shallow water equations are
solved using a second-order well-balanced discontinuous Galerkin
scheme. Theoretical test case and three flood events are described
to demonstrate the potential benefits of the scheme: (i) wetting and
drying in a parabolic basin (ii) flash flood over a physical model of
the urbanized Toce River valley in Italy; (iii) wave propagation on
the Reyran river valley in consequence of the Malpasset dam-break
in 1959 (France); and (iv) dam-break flood in October 1982 at the
town of Sumacarcel (Spain). The capability of the scheme is also
verified against alternative models. Computational results compare
well with recorded data and show that the scheme is at least as
efficient as comparable second-order finite volume schemes, with
notable efficiency speedup due to parallelization.
Abstract: We present and analyze reliable numerical techniques
for simulating complex flow and transport phenomena related to
natural gas transportation in pipelines. Such kind of problems
are of high interest in the field of petroleum and environmental
engineering. Modeling and understanding natural gas flow and
transformation processes during transportation is important for the
sake of physical realism and the design and operation of pipeline
systems. In our approach a two fluid flow model based on a system
of coupled hyperbolic conservation laws is considered for describing
natural gas flow undergoing hydratization. The accurate numerical
approximation of two-phase gas flow remains subject of strong
interest in the scientific community. Such hyperbolic problems are
characterized by solutions with steep gradients or discontinuities, and
their approximation by standard finite element techniques typically
gives rise to spurious oscillations and numerical artefacts. Recently,
stabilized and discontinuous Galerkin finite element techniques
have attracted researchers’ interest. They are highly adapted to the
hyperbolic nature of our two-phase flow model. In the presentation
a streamline upwind Petrov-Galerkin approach and a discontinuous
Galerkin finite element method for the numerical approximation of
our flow model of two coupled systems of Euler equations are
presented. Then the efficiency and reliability of stabilized continuous
and discontinous finite element methods for the approximation is
carefully analyzed and the potential of the either classes of numerical
schemes is investigated. In particular, standard benchmark problems
of two-phase flow like the shock tube problem are used for the
comparative numerical study.
Abstract: This paper deals with a high-order accurate Runge
Kutta Discontinuous Galerkin (RKDG) method for the numerical
solution of the wave equation, which is one of the simple case of a
linear hyperbolic partial differential equation. Nodal DG method is
used for a finite element space discretization in 'x' by discontinuous
approximations. This method combines mainly two key ideas which
are based on the finite volume and finite element methods. The
physics of wave propagation being accounted for by means of
Riemann problems and accuracy is obtained by means of high-order
polynomial approximations within the elements. High order accurate
Low Storage Explicit Runge Kutta (LSERK) method is used for
temporal discretization in 't' that allows the method to be nonlinearly
stable regardless of its accuracy. The resulting RKDG
methods are stable and high-order accurate. The L1 ,L2 and L∞ error
norm analysis shows that the scheme is highly accurate and effective.
Hence, the method is well suited to achieve high order accurate
solution for the scalar wave equation and other hyperbolic equations.
Abstract: This paper is concerned with the numerical minimization
of energy functionals in BV (
) (the space of bounded variation
functions) involving total variation for gray-scale 1-dimensional inpainting
problem. Applications are shown by finite element method
and discontinuous Galerkin method for total variation minimization.
We include the numerical examples which show the different recovery
image by these two methods.
Abstract: A water surface slope limiting scheme is tested and
compared with the water depth slope limiter for the solution of one
dimensional shallow water equations with bottom slope source term.
Numerical schemes based on the total variation diminishing Runge-
Kutta discontinuous Galerkin finite element method with slope
limiter schemes based on water surface slope and water depth are
used to solve one-dimensional shallow water equations. For each
slope limiter, three different Riemann solvers based on HLL, LF, and
Roe flux functions are used. The proposed water surface based slope
limiter scheme is easy to implement and shows better conservation
property compared to the slope limiter based on water depth. Of the
three flux functions, the Roe approximation provides the best results
while the LF function proves to be least suitable when used with
either slope limiter scheme.
Abstract: We deal with the numerical solution of time-dependent convection-diffusion-reaction equations. We combine the local projection stabilization method for the space discretization with two different time discretization schemes: the continuous Galerkin-Petrov (cGP) method and the discontinuous Galerkin (dG) method of polynomial of degree k. We establish the optimal error estimates and present numerical results which shows that the cGP(k) and dG(k)- methods are accurate of order k +1, respectively, in the whole time interval. Moreover, the cGP(k)-method is superconvergent of order 2k and dG(k)-method is of order 2k +1 at the discrete time points. Furthermore, the dependence of the results on the choice of the stabilization parameter are discussed and compared.