Bojan Crnković

High order accurate semi-implicit WENO schemes for hyperbolic balance laws

Abstract In this paper we propose a family of well-balanced semi-implicit numerical schemes for hyperbolic conservation and balance laws. The basic idea of the proposed schemes lies in the combination of the finite volume WENO discretization with Roe’s solver and the strong stability preserving (SSP) time integration methods, which ensure the stability properties of the considered schemes [S. Gottlieb, C.-W. Shu, E. Tadmor, Strong stability-preserving high-order time discretization methods, SIAM Rev. 43 (2001) 89–112]. While standard WENO schemes typically use explicit time integration methods, in this paper we are combining WENO spatial discretization with optimal SSP singly diagonally implicit (SDIRK) methods developed in [L. Ferracina, M.N. Spijker, Strong stability of singly diagonally implicit Runge–Kutta methods, Appl. Numer. Math. 58 (2008) 1675–1686]. In this way the implicit WENO numerical schemes are obtained. In order to reduce the computational effort, the implicit part of the numerical scheme is linearized in time by taking into account the complete WENO reconstruction procedure. With the proposed linearization the new semi-implicit finite volume WENO schemes are designed. A detailed numerical investigation of the proposed numerical schemes is presented in the paper. More precisely, schemes are tested on one-dimensional linear scalar equation and on non-linear conservation law systems. Furthermore, well-balanced semi-implicit WENO schemes for balance laws with geometrical source terms are defined. Such schemes are then applied to the open channel flow equations. We prove that the defined numerical schemes maintain steady state solution of still water. The application of the new schemes to different open channel flow examples is shown.

Generalized distributed order diffusion equations with composite time fractional derivative

Abstract In this paper we investigate the solution of generalized distributed order diffusion equations with composite time fractional derivative by using the Fourier–Laplace transform method. We represent solutions in terms of infinite series in Fox H -functions. The fractional and second moments are derived by using Mittag-Leffler functions. We observe decelerating anomalous subdiffusion in case of two composite time fractional derivatives. Generalized uniformly distributed order diffusion equation, as a model for strong anomalous behavior, is analyzed by using Tauberian theorem. Some previously obtained results are special cases of those presented in this paper.

Ergodicity-Based Cooperative Multiagent Area Coverage via a Potential Field

This paper considers a problem of area coverage where the objective is to achieve given coverage density by use of multiple mobile agents. We present an ergodicity-based coverage algorithm which enables a centralized feedback control for multiagent system based on radial basis function (RBF) representation of the ergodicity problem and a solution of an appropriately designed stationary heat equation for the potential field. The heat equation uses a source term that depends on the difference between the given goal density distribution and the current coverage density (time average of RBFs along trajectories). The agent movement is directed using the gradient of that potential field. The heat equation driven area coverage has a built-in cooperative behavior of agents which includes collision avoidance and coverage coordination. The algorithm is robust, scalable, and computationally inexpensive.

Binary weighted essentially non-oscillatory (BWENO) approximation

Weighted essentially non-oscillatory (WENO) schemes have been mainly used for solving hyperbolic partial differential equations (PDEs). Such schemes are capable of high order approximation in smooth regions and non-oscillatory sharp resolution of discontinuities. The base of the WENO schemes is a non-oscillatory WENO approximation procedure, which is not necessarily related to PDEs. The typical WENO procedures are WENO interpolation and WENO reconstruction. The WENO algorithm has gained much popularity but the basic idea of approximation did not change much over the years. In this paper, we first briefly review the idea of WENO interpolation and propose a modification of the basic algorithm. New approximation should improve basic characteristics of the approximation and provide a more flexible framework for future applications. New WENO procedure involves a binary tree weighted construction that is based on key ideas of WENO algorithm and we refer to it as the binary weighted essentially non-oscillatory (BWENO) approximation. New algorithm comes in a rational and a polynomial version. Furthermore, we describe the WENO reconstruction procedure, which is usually involved in the numerical schemes for hyperbolic PDEs, and propose the new reconstruction procedure based on the described BWENO interpolation. The obtained numerical results show that the newly proposed procedures perform very well on the considered test examples.