Related papers: A C++ Code to Solve the DGLAP Equations Applied to…
Ultra high energy cosmic rays (UHECR) pose a problem either for particle physics or for astrophysics (or for both) by the unexpectedly high number of cosmic ray showers observed with energy above about 5x10^{19}eV, the…
Solving the question of the origin of ultra-high energy cosmic rays (UHECRs) requires the development of detailed simulation tools in order to interpret the experimental data and draw conclusions on the UHECR universe. CRPropa is a public…
*The gradient discretisation method (GDM) is a generic framework, covering many classical methods (Finite Elements, Finite Volumes, Discontinuous Galerkin, etc.), for designing and analysing numerical schemes for diffusion models. In this…
In this paper, we introduce novel discontinuous Galerkin (DG) schemes for the Cahn-Hilliard equation, which arises in many applications. The method is designed by integrating the mixed DG method for the spatial discretization with the…
In this work numerical methods for solving Einstein's equations are developed and applied to the study of inhomogeneous cosmological models. A two-dimensional computer code is described which implements two advanced numerical methods:…
We use a kinetic-equation approach to describe the propagation of ultra high energy cosmic ray protons and nuclei comparing theoretical results with the observations of the Pierre Auger Observatory.
Dense, discrete Graphical Models with pairwise potentials are a powerful class of models which are employed in state-of-the-art computer vision and bio-imaging applications. This work introduces a new MAP-solver, based on the popular Dual…
State of the art domain decomposition algorithms for large-scale boundary value problems (with $M\gg 1$ degrees of freedom) suffer from bounded strong scalability because they involve the synchronisation and communication of workers…
A new computational algorithm, the discrete singular convolution (DSC), is introduced for computational electromagnetics. The basic philosophy behind the DSC algorithm for the approximation of functions and their derivatives is studied.…
Understanding the anisotropy of ultra high-energy cosmic rays (UHECRs) is crucial for unraveling the origins and propagation mechanisms of these enigmatic particles. In this work, we studied the dipolar anisotropy of UHECRs in the diffusive…
Superheavy particles of masses $\simeq 10^{13}-10^{14} GeV$ with lifetimes $\simeq 10^{10}-10^{22} years$ are very interesting, since their decays may account for the ultra-high energy (UHE) cosmic rays discovered beyond the…
A method of numerically solving the Maxwell equations is considered for modeling harmonic electromagnetic fields. The vector finite element method makes it possible to obtain a physically consistent discretization of the differential…
Recently found all the fundamental solutions of a multidimensional singular elliptic equation are expressed in terms of the well-known Lauricella hypergeometric function in many variables. In this paper, we find a unique solution of the…
In the first chapter we give an introduction to hard diffractive scattering in QCD to introduce basic concepts and terminology. In the second chapter we make predictions for the evolution of skewed parton distributions in a proton in the…
The modified evolution equation for parton distributions of Dokshitzer, Marchesini and Salam is extended to non-singlet Deep Inelastic Scattering coefficient functions and the physical evolution kernels which govern their scaling violation.…
In this paper, we present convergence analysis of high-order finite element based methods, in particular, we focus on a discontinuous Galerkin scheme using summation-by-parts operators. To this end, it is crucial that structure preserving…
In Generalized Linear Estimation (GLE) problems, we seek to estimate a signal that is observed through a linear transform followed by a component-wise, possibly nonlinear and noisy, channel. In the Bayesian optimal setting, Generalized…
We propose two new alternative numerical schemes to solve the coupled Einstein-Euler equations in the Generalized Harmonic formulation. The first one is a finite difference (FD) Central Weighted Essentially Non-Oscillatory (CWENO) scheme on…
A compiler approach for generating low-level computer code from high-level input for discontinuous Galerkin finite element forms is presented. The input language mirrors conventional mathematical notation, and the compiler generates…
In this paper, we study superconvergence properties of the discontinuous Galerkin (DG) method for one-dimensional linear hyperbolic equation when upwind fluxes are used. We prove, for any polynomial degree $k$, the $2k+1$th (or $2k+1/2$th)…