An object-oriented numerical integration template implementation is presented on the basis of the CCC programming language. Aiming its straightforward application in finite and boundary element methods, the design supports integrand objects of scalar, vector or matrix types, so that a single programming statement is able to integrate element matrices and vectors. The integrand can contain singu...

Historically the finite volume methods have been developed for the numerical integration of conservation laws. In this study we present some recent results on the application of such schemes to dispersive PDEs. Namely, we solve numerically a representative of Boussinesq type equations in view of important applications to the coastal hydrodynamics. Numerical results of the runup of a moderate wa...

A numerical method to solve the fractional diffusion equation, which could also be easily extended to many other fractional dynamics equations, is considered. These fractional equations have been proposed in order to describe anomalous transport characterized by non-Markovian kinetics and the breakdown of Fick’s law. In this paper we combine the forward time centered space (FTCS) method, well k...

In this paper Gauss-Radau and Gauss-Lobatto quadrature rules are presented to evaluate the rational integrals of the element matrix for a general quadrilateral. These integrals arise in finite element formulation for second order Partial Differential Equation via Galerkin weighted residual method in closed form. Convergence to the analytical solutions and efficiency are depicted by numerical ex...

In [16, 17], we constructed uniformly high order accurate discontinuous Galerkin (DG) schemes which preserve positivity of density and pressure for the Euler equations of compressible gas dynamics with the ideal gas equation of state. The technique also applies to high order accurate finite volume schemes. For the Euler equations with various source terms (e.g., gravity and chemical reactions),...

We consider a standard model for mass transport across an evolving interface. The solution has to satisfy a jump condition across an evolving interface. In the recent paper [C. Lehrenfeld, A. Reusken, SIAM J. Num. Anal., 51 (2013)] a new finite element discretization method for this mass transport problem is presented and analyzed. This method is based on a space-time approach in which a discon...

This paper describes a simple but effective technique for reducing dispersion errors in finite element solutions of timeharmonic wave propagation problems. The method involves a simple shift of the integration points to locations away from conventional Gauss or Gauss–Lobatto integration points. For bilinear rectangular elements, such a shift results in fourth-order accuracy with respect to disp...

We derive a posteriori error estimates for two classes of explicit finite difference schemes for ordinary differential equations. To facilitate the analysis, we derive a systematic reformulation of the finite difference schemes as finite element methods. The a posteriori error estimates quantify various sources of discretization errors, including effects arising from explicit discretization. Th...

In this paper, the extended finite element method (XFEM) is employed to investigate the statics and vibration problems of cracked isotropic bars and beams. Three kinds of elements namely the standard, the blended and the enriched elements are utilized to discretize the structure and model cracks. Two techniques referred as the increase of the number of Gauss integration points and the rectangle...

in this paper, an enhanced computational code was developed using finite-volume method for solving the incompressible natural convection flow within the cylindrical cavities. grids were generated by an easy method with a view to computer program providing. an explicit integration algorithm was applied to find the steady state condition. also instead of the conventional algorithms of simple, sim...