A Conservative Scheme with Optimal Error Estimates for a Multidimensional Space-Fractional Gross-Pitaevskii Equation

The present work departs from an extended form of the classical multi-dimensional Gross-Pitaevskii equation, which considers fractional derivatives of the Riesz type in space, a generalized potential function and angular momentum rotation. It is well known that the classical system possesses functionals which are preserved throughout time. It is easy to check that the generalized fractional model considered in this work also possesses conserved quantities, whence the development of conservative and efficient numerical schemes is pragmatically justified. Motivated by these facts, we propose a finite-difference method based on weighted-shifted Grünwald differences to approximate the solutions of the generalized Gross-Pitaevskii system. We provide here a discrete extension of the uniform Sobolev inequality to multiple dimensions, and show that the proposed method is capable of preserving discrete forms of the mass and the energy of the model. Moreover, we establish thoroughly the stability and the convergence of the technique, and provide some illustrative simulations to show that the method is capable of preserving the total mass and the total energy of the generalized system. © 2019 Ahmed S. Hendy et al., published by Sciendo 2019

Exponential integrators: tensor structured problems and applications

The solution of stiff systems of Ordinary Differential Equations (ODEs), that typically arise after spatial discretization of many important evolutionary Partial Differential Equations (PDEs), constitutes a topic of wide interest in numerical analysis. A prominent way to numerically integrate such systems involves using exponential integrators. In general, these kinds of schemes do not require the solution of (non)linear systems but rather the action of the matrix exponential and of some specific exponential-like functions (known in the literature as phi-functions). In this PhD thesis we aim at presenting efficient tensor-based tools to approximate such actions, both from a theoretical and from a practical point of view, when the problem has an underlying Kronecker sum structure. Moreover, we investigate the application of exponential integrators to compute numerical solutions of important equations in various fields, such as plasma physics, mean-field optimal control and computational chemistry. In any case, we provide several numerical examples and we perform extensive simulations, eventually exploiting modern hardware architectures such as multi-core Central Processing Units (CPUs) and Graphic Processing Units (GPUs). The results globally show the effectiveness and the superiority of the different approaches proposed

Spectral and High Order Methods for Partial Differential Equations ICOSAHOM 2018

This open access book features a selection of high-quality papers from the presentations at the International Conference on Spectral and High-Order Methods 2018, offering an overview of the depth and breadth of the activities within this important research area. The carefully reviewed papers provide a snapshot of the state of the art, while the extensive bibliography helps initiate new research directions

Generalized averaged Gaussian quadrature and applications

A simple numerical method for constructing the optimal generalized averaged Gaussian quadrature formulas will be presented. These formulas exist in many cases in which real positive GaussKronrod formulas do not exist, and can be used as an adequate alternative in order to estimate the error of a Gaussian rule. We also investigate the conditions under which the optimal averaged Gaussian quadrature formulas and their truncated variants are internal

MS FT-2-2 7 Orthogonal polynomials and quadrature: Theory, computation, and applications

Quadrature rules find many applications in science and engineering. Their analysis is a classical area of applied mathematics and continues to attract considerable attention. This seminar brings together speakers with expertise in a large variety of quadrature rules. It is the aim of the seminar to provide an overview of recent developments in the analysis of quadrature rules. The computation of error estimates and novel applications also are described