Sparse spectral methods for integral equations and equilibrium measures

In this thesis, we introduce new numerical approaches to two important types of integral equation problems using sparse spectral methods. First, linear as well as nonlinear Volterra integral and integro-differential equations and second, power-law integral equations on d-dimensional balls involved in the solution of equilibrium measure problems. These methods are based on ultraspherical spectral methods and share key properties and advantages as a result of their joint starting point: By working in appropriately weighted orthogonal Jacobi polynomial bases, we obtain recursively generated banded operators allowing us to obtain high precision solutions at low computational cost. This thesis consists of three chapters in which the background of the above-mentioned problems and methods are respectively introduced in the context of their mathematical theory and applications, the necessary results to construct the operators and obtain solutions are proved and the method's applicability and efficiency are showcased by comparing them with current state-of-the-art approaches and analytic results where available. The first chapter gives a general scope introduction to sparse spectral methods using Jacobi polynomials in one and higher dimensions. The second chapter concerns the numerical solution of Volterra integral equations. The introduced method achieves exponential convergence and works for general kernels, a major advantage over comparable methods which are limited to convolution kernels. The third chapter introduces an approximately banded method to solve power law kernel equilibrium measures in arbitrary dimensional balls. This choice of domain is suggested by the radial symmetry of the problem and analytic results on the supports of the resulting measures. For our method, we obtain the crucial property of computational cost independent of the dimension of the domain, a major contrast to particle simulations which are the current standard approach to these problems and scale extremely poorly with both the dimension and the number of particles.Open Acces

Fast Spectral Collocation Method for Solving Nonlinear Time-Delayed Burgers-Type Equations with Positive Power Terms

Since the collocation method approximates ordinary differential equations, partial differential equations, and integral equations in physical space, it is very easy to implement and adapt to various problems, including variable coefficient and nonlinear differential equations. In this paper, we derive a Jacobi-Gauss-Lobatto collocation method (J-GL-C) to solve numerically nonlinear time-delayed Burgers-type equations. The proposed technique is implemented in two successive steps. In the first one, we apply nodes of the Jacobi-Gauss-Lobatto quadrature which depend upon the two general parameters , and the resulting equations together with the two-point boundary conditions constitute a system of ordinary differential equations (ODEs) in time. In the second step, the implicit Runge-Kutta method of fourth order is applied to solve a system of ODEs of second order in time. We present numerical results which illustrate the accuracy and flexibility of these algorithms

On the Inversion of High Energy Proton

Inversion of the K-fold stochastic autoconvolution integral equation is an elementary nonlinear problem, yet there are no de facto methods to solve it with finite statistics. To fix this problem, we introduce a novel inverse algorithm based on a combination of minimization of relative entropy, the Fast Fourier Transform and a recursive version of Efron's bootstrap. This gives us power to obtain new perspectives on non-perturbative high energy QCD, such as probing the ab initio principles underlying the approximately negative binomial distributions of observed charged particle final state multiplicities, related to multiparton interactions, the fluctuating structure and profile of proton and diffraction. As a proof-of-concept, we apply the algorithm to ALICE proton-proton charged particle multiplicity measurements done at different center-of-mass energies and fiducial pseudorapidity intervals at the LHC, available on HEPData. A strong double peak structure emerges from the inversion, barely visible without it.Comment: 29 pages, 10 figures, v2: extended analysis (re-projection ratios, 2D

Numerical homogenization of time-dependent Maxwell\u27s equations with dispersion effects

This thesis studies the propagation of electromagnetic waves in heterogeneous structures such as metamaterials. The governing equations for this problem are Maxwell\u27s equations with highly oscillatory parameters. We use an analytic homogenization result which yields an effective Maxwell system that involves a convolution integral. This convolution represents dispersive effects that result from the interaction of the wave with the (locally) periodic microscopic structure. We discretize in space using the Finite Element Heterogeneous Multiscale Method (FE-HMM) and provide a semi-discrete error estimate. The rigorous error analysis in space is supplemented by a rather standard time discretization at the end of which an efficient, fully discrete method is proposed. This method uses a recursive approximation of the convolution that relies on the assumption that the convolution kernel is an exponential function. Eventually, we present numerical experiments both for the microscopic and the macroscopic scale

New Challenges Arising in Engineering Problems with Fractional and Integer Order

Mathematical models have been frequently studied in recent decades, in order to obtain the deeper properties of real-world problems. In particular, if these problems, such as finance, soliton theory and health problems, as well as problems arising in applied science and so on, affect humans from all over the world, studying such problems is inevitable. In this sense, the first step in understanding such problems is the mathematical forms. This comes from modeling events observed in various fields of science, such as physics, chemistry, mechanics, electricity, biology, economy, mathematical applications, and control theory. Moreover, research done involving fractional ordinary or partial differential equations and other relevant topics relating to integer order have attracted the attention of experts from all over the world. Various methods have been presented and developed to solve such models numerically and analytically. Extracted results are generally in the form of numerical solutions, analytical solutions, approximate solutions and periodic properties. With the help of newly developed computational systems, experts have investigated and modeled such problems. Moreover, their graphical simulations have also been presented in the literature. Their graphical simulations, such as 2D, 3D and contour figures, have also been investigated to obtain more and deeper properties of the real world problem

International Conference on Mathematical Analysis and Applications in Science and Engineering – Book of Extended Abstracts

The present volume on Mathematical Analysis and Applications in Science and Engineering - Book of Extended Abstracts of the ICMASC’2022 collects the extended abstracts of the talks presented at the International Conference on Mathematical Analysis and Applications in Science and Engineering – ICMA2SC'22 that took place at the beautiful city of Porto, Portugal, in June 27th-June 29th 2022 (3 days). Its aim was to bring together researchers in every discipline of applied mathematics, science, engineering, industry, and technology, to discuss the development of new mathematical models, theories, and applications that contribute to the advancement of scientific knowledge and practice. Authors proposed research in topics including partial and ordinary differential equations, integer and fractional order equations, linear algebra, numerical analysis, operations research, discrete mathematics, optimization, control, probability, computational mathematics, amongst others. The conference was designed to maximize the involvement of all participants and will present the state-of- the-art research and the latest achievements.info:eu-repo/semantics/publishedVersio

Trends, directions for further research, and some open problems of fractional calculus

The area of fractional calculus (FC) has been fast developing and is presently being applied in all scientific fields. Therefore, it is of key relevance to assess the present state of development and to foresee, if possible, the future evolution, or, at least, the challenges identified in the scope of advanced research works. This paper gives a vision about the directions for further research as well as some open problems of FC. A number of topics in mathematics, numerical algorithms and physics are analyzed, giving a systematic perspective for future research.info:eu-repo/semantics/publishedVersio