-
Efficient numerical methods for the Maxey-Riley equations with Basset history term
Authors:
Julio Urizarna-Carasa,
Leon Schlegel,
Daniel Ruprecht
Abstract:
The Maxey-Riley equations (MRE) describe the motion of a finite-sized, spherical particle in a fluid. Because of wake effects, the force acting on a particle depends on its past trajectory. This is modelled by an integral term in the MRE, also called Basset force, that makes its numerical solution challenging and memory intensive. A recent approach proposed by Prasath, Vasan and Govindarajan explo…
▽ More
The Maxey-Riley equations (MRE) describe the motion of a finite-sized, spherical particle in a fluid. Because of wake effects, the force acting on a particle depends on its past trajectory. This is modelled by an integral term in the MRE, also called Basset force, that makes its numerical solution challenging and memory intensive. A recent approach proposed by Prasath, Vasan and Govindarajan exploits connections between the integral term and fractional derivatives to reformulate the MRE as a time-dependent partial differential equation on a semi-infinite pseudo-space. They also propose a numerical algorithm based on polynomial expansions. This paper develops a numerical approach based on finite difference instead, by adopting techniques by Koleva and Fazio and Janelli to cope with the issues of having an unbounded spatial domain. We compare convergence order and computational efficiency for particles of varying size and density of the polynomial expansion by Prasath et al., our finite difference schemes and a direct integrator for the MRE based on multi-step methods proposed by Daitche.
△ Less
Submitted 20 March, 2024;
originally announced March 2024.
-
CRPropa 3.2: a public framework for high-energy astroparticle simulations
Authors:
Sophie Aerdker,
Rafael Alves Batista,
Julia Becker Tjus,
Julien Dörner,
Andrej Dundovic,
Björn Eichmann,
Antonius Frie,
Christopher Heiter,
Mario Hoerbe,
Karl-Heinz Kampert,
Lukas Merten,
Gero Müller,
Patrick Reichherzer,
Simone Rossoni,
Andrey Saveliev,
Leander Schlegel,
Günter Sigl,
Arjen van Vliet,
Tobias Winchen
Abstract:
CRPropa is a Monte Carlo framework for simulating the propagation of (ultra-) high-energy particles in the Universe, including cosmic rays, gamma rays, electrons, and neutrinos. It covers energies from ZeV down to GeV for gamma rays and electrons, and TeV for cosmic rays and neutrinos, supporting various astrophysical environments such as the surroundings of astrophysical sources, galactic, and ex…
▽ More
CRPropa is a Monte Carlo framework for simulating the propagation of (ultra-) high-energy particles in the Universe, including cosmic rays, gamma rays, electrons, and neutrinos. It covers energies from ZeV down to GeV for gamma rays and electrons, and TeV for cosmic rays and neutrinos, supporting various astrophysical environments such as the surroundings of astrophysical sources, galactic, and extragalactic environments. The newest version, CRPropa 3.2, represents a significant leap forward towards a universal multi-messenger framework, opening up the possibility for many more astrophysical applications. This includes extensions to simulate cosmic-ray acceleration and particle interactions within astrophysical source environments, a full Monte Carlo treatment of electromagnetic cascades, improved ensemble-averaged Galactic propagation, significant performance improvements for cosmic-ray tracking through magnetic fields, and a user-friendly implementation of custom photon fields, among many more enhancements. This contribution will give an overview of the new features and present several applications to cosmic-ray and gamma-ray propagation.
△ Less
Submitted 18 August, 2023;
originally announced August 2023.
-
Multimessenger Picture of J1048+7143
Authors:
Emma Kun,
Ilja Jaroschewski,
Armin Ghorbanietemad,
Sándor Frey,
Julia Becker Tjus,
Silke Britzen,
Krisztina Éva Gabányi,
Vladimir Kiselev,
Leander Schlegel,
Marcel Schroller,
Patrick Reichherzer,
Lang Cui,
Xin Wang,
Yuling Shen
Abstract:
We draw the multimessenger picture of J1048+7143, a flat-spectrum radio quasar known to show quasi-periodic oscillations in the $γ$-ray regime. We generate the adaptively-binned Fermi Large Area Telescope light curve of this source above 168 MeV to find three major $γ$-ray flares of the source, such that all three flares consist of two-two sharp sub-flares. Based on radio interferometric imaging d…
▽ More
We draw the multimessenger picture of J1048+7143, a flat-spectrum radio quasar known to show quasi-periodic oscillations in the $γ$-ray regime. We generate the adaptively-binned Fermi Large Area Telescope light curve of this source above 168 MeV to find three major $γ$-ray flares of the source, such that all three flares consist of two-two sharp sub-flares. Based on radio interferometric imaging data taken with the Very Large Array, we find that the kpc-scale jet is directed towards west, while our analysis of $8.6$-GHz very long baseline interferometry data, mostly taken with the Very Long Baseline Array, revealed signatures of two pc-scale jets, one pointing towards east, one pointing towards south. We suggest that the misalignment of the kpc- and pc-scale jets is a revealing signature of jet precession. We also analyze the $5$-GHz total flux density curve of J1048+7143 taken with the Nanshan(Ur) and RATAN-600 single dish radio telescopes and find two complete radio flares, slightly lagging behind the $γ$-ray flares. We model the timing of $γ$-ray flares as signature of the spin-orbit precession in a supermassive black hole binary, and find that the binary could merge in the next $\sim 60-80$ years. We show that both the Pulsar Timing Arrays and the planned Laser Interferometer Space Antenna lack sensitivity and frequency coverage to detect the hypothetical supermassive black hole binary in J1048$+$7143. We argue that the identification of sources similar to J1048+7143 plays a key role to reveal periodic high-energy sources in the distant Universe.
△ Less
Submitted 22 October, 2022; v1 submitted 12 September, 2022;
originally announced September 2022.
-
CRPropa 3.2 -- an advanced framework for high-energy particle propagation in extragalactic and galactic spaces
Authors:
Rafael Alves Batista,
Julia Becker Tjus,
Julien Dörner,
Andrej Dundovic,
Björn Eichmann,
Antonius Frie,
Christopher Heiter,
Mario R. Hoerbe,
Karl-Heinz Kampert,
Lukas Merten,
Gero Müller,
Patrick Reichherzer,
Andrey Saveliev,
Leander Schlegel,
Günter Sigl,
Arjen van Vliet,
Tobias Winchen
Abstract:
The landscape of high- and ultra-high-energy astrophysics has changed in the last decade, largely due to the inflow of data collected by large-scale cosmic-ray, gamma-ray, and neutrino observatories. At the dawn of the multimessenger era, the interpretation of these observations within a consistent framework is important to elucidate the open questions in this field. CRPropa 3.2 is a Monte Carlo c…
▽ More
The landscape of high- and ultra-high-energy astrophysics has changed in the last decade, largely due to the inflow of data collected by large-scale cosmic-ray, gamma-ray, and neutrino observatories. At the dawn of the multimessenger era, the interpretation of these observations within a consistent framework is important to elucidate the open questions in this field. CRPropa 3.2 is a Monte Carlo code for simulating the propagation of high-energy particles in the Universe. This version represents a major leap forward, significantly expanding the simulation framework and opening up the possibility for many more astrophysical applications. This includes, among others: efficient simulation of high-energy particles in diffusion-dominated domains, self-consistent and fast modelling of electromagnetic cascades with an extended set of channels for photon production, and studies of cosmic-ray diffusion tensors based on updated coherent and turbulent magnetic-field models. Furthermore, several technical updates and improvements are introduced with the new version, such as: enhanced interpolation, targeted emission of sources, and a new propagation algorithm (Boris push). The detailed description of all novel features is accompanied by a discussion and a selected number of example applications.
△ Less
Submitted 29 July, 2022;
originally announced August 2022.
-
Shape Optimization for the Mitigation of Coastal Erosion via Smoothed Particle Hydrodynamics
Authors:
Luka Schlegel,
Volker Schulz
Abstract:
Adjoint-based shape optimization most often relies on Eulerian flow field formulations. However, since Lagrangian particle methods are the natural choice for solving sedimentation problems in oceanography, extensions to the Lagrangian framework are desirable. For the mitigation of coastal erosion, we perform shape optimization for fluid flows, that are described by Lagrangian shallow water equatio…
▽ More
Adjoint-based shape optimization most often relies on Eulerian flow field formulations. However, since Lagrangian particle methods are the natural choice for solving sedimentation problems in oceanography, extensions to the Lagrangian framework are desirable. For the mitigation of coastal erosion, we perform shape optimization for fluid flows, that are described by Lagrangian shallow water equations and discretized via smoothed particle hydrodynamics. The obstacle's shape is hereby optimized over an appropriate cost function to minimize the height of water waves along the shoreline based on shape calculus. Theoretical results will be numerically verified by exploring different scenarios.
△ Less
Submitted 6 September, 2022; v1 submitted 18 March, 2022;
originally announced March 2022.
-
Shape Optimization for the Mitigation of Coastal Erosion via Porous Shallow Water Equations
Authors:
Luka Schlegel,
Volker Schulz
Abstract:
Coastal erosion describes the displacement of land caused by destructive sea waves, currents or tides. Major efforts have been made to mitigate these effects using groynes, breakwaters and various other structures. We address this problem by applying shape optimization techniques on the obstacles. We model the propagation of waves towards the coastline using two-dimensional porous Shallow Water Eq…
▽ More
Coastal erosion describes the displacement of land caused by destructive sea waves, currents or tides. Major efforts have been made to mitigate these effects using groynes, breakwaters and various other structures. We address this problem by applying shape optimization techniques on the obstacles. We model the propagation of waves towards the coastline using two-dimensional porous Shallow Water Equations with artificial viscosity. The obstacle's shape, which is assumed to be permeable, is optimized over an appropriate cost function to minimize the height and velocities of water waves along the shore, without relying on a finite-dimensional design space, but based on shape calculus.
△ Less
Submitted 27 July, 2022; v1 submitted 29 November, 2021;
originally announced November 2021.
-
Shape Optimization for the Mitigation of Coastal Erosion via the Helmholtz Equation
Authors:
Luka Schlegel,
Volker Schulz
Abstract:
Coastal erosion describes the displacement of land caused by destructive sea waves, currents or tides. Major efforts have been made to mitigate these effects using groins, breakwaters and various other structures. We try to address this problem by applying shape optimization techniques on the obstacles. A first approach models the propagation of waves towards the coastline, using a 2D time-harmoni…
▽ More
Coastal erosion describes the displacement of land caused by destructive sea waves, currents or tides. Major efforts have been made to mitigate these effects using groins, breakwaters and various other structures. We try to address this problem by applying shape optimization techniques on the obstacles. A first approach models the propagation of waves towards the coastline, using a 2D time-harmonic system based on the famous Helmholtz equation in the form of a scattering problem. The obstacle's shape is optimized over an appropriate cost function to minimize the height of water waves along the shoreline, without relying on a finite-dimensional design space, but based on shape calculus.
△ Less
Submitted 17 January, 2022; v1 submitted 21 July, 2021;
originally announced July 2021.
-
Shape Optimization for the Mitigation of Coastal Erosion via Shallow Water Equations
Authors:
Luka Schlegel,
Volker Schulz
Abstract:
Coastal erosion describes the displacement of land caused by destructive sea waves, currents or tides. Major efforts have been made to mitigate these effects using groins, breakwaters and various other structures. We try to address this problem by applying shape optimization techniques to the obstacles. We model the propagation of waves towards the coastline, using two-dimensional shallow water eq…
▽ More
Coastal erosion describes the displacement of land caused by destructive sea waves, currents or tides. Major efforts have been made to mitigate these effects using groins, breakwaters and various other structures. We try to address this problem by applying shape optimization techniques to the obstacles. We model the propagation of waves towards the coastline, using two-dimensional shallow water equations. The obstacle's shape is optimized over an appropriate cost function to minimize the height and velocities of water waves along the shore, without relying on a finite-dimensional design space but based on shape calculus.
△ Less
Submitted 15 September, 2022; v1 submitted 20 July, 2021;
originally announced July 2021.
-
CRPropa 3.2: a framework for high-energy astroparticle propagation
Authors:
Rafael Alves Batista,
Julia Becker Tjus,
Julien Dörner,
Andrej Dundovic,
Björn Eichmann,
Antonius Frie,
Christopher Heiter,
Mario R. Hoerbe,
Karl-Heinz Kampert,
Lukas Merten,
Gero Müller,
Patrick Reichherzer,
Andrey Saveliev,
Leander Schlegel,
Günter Sigl,
Arjen van Vliet,
Tobias Winchen
Abstract:
The landscape of high- and ultra-high-energy astrophysics has changed in the last decade, in large part owing to the inflow of high-quality data collected by present cosmic-ray, gamma-ray, and neutrino observatories. At the dawn of the multimessenger era, the interpretation of these observations within a consistent framework is important to elucidate the open questions in this field. CRPropa 3.2 i…
▽ More
The landscape of high- and ultra-high-energy astrophysics has changed in the last decade, in large part owing to the inflow of high-quality data collected by present cosmic-ray, gamma-ray, and neutrino observatories. At the dawn of the multimessenger era, the interpretation of these observations within a consistent framework is important to elucidate the open questions in this field. CRPropa 3.2 is a Monte Carlo code for simulating the propagation of high-energy particles in the Universe. This new version represents a step further towards a more complete simulation framework for multimessenger studies. Some of the new developments include: cosmic-ray acceleration, support for particle interactions within astrophysical sources, full Monte Carlo treatment of electromagnetic cascades, improved ensemble-averaged Galactic propagation, and a number of technical enhancements. Here we present some of these novel features and some applications to gamma- and cosmic-ray propagation.
△ Less
Submitted 4 July, 2021;
originally announced July 2021.
-
Interpolation of turbulent magnetic fields and its consequences on diffusive cosmic ray propagation
Authors:
L. Schlegel,
A. Frie,
B. Eichmann,
P. Reichherzer,
J. Becker Tjus
Abstract:
Numerical simulations of the propagation of charged particles through magnetic fields solving the equation of motion often leads to the usage of an interpolation in case of discretely defined magnetic fields, typically given on a homogeneous grid structure. However, the interpolation method influences the magnetic field properties on the scales of the grid spacing and the choice of interpolation r…
▽ More
Numerical simulations of the propagation of charged particles through magnetic fields solving the equation of motion often leads to the usage of an interpolation in case of discretely defined magnetic fields, typically given on a homogeneous grid structure. However, the interpolation method influences the magnetic field properties on the scales of the grid spacing and the choice of interpolation routine can therefore change the result. At the same time, it provides an impact, i.e.\ error, on the spatial particle distribution. We compare three different interpolation routines -- trilinear, tricubic and nearest neighbor interpolation -- in the case of turbulent magnetic fields and show that there is no benefit in using trilinear interpolation. We show that in comparison, the nearest neighbor interpolation provides the best performance, i.e.\ requires least CPU time and results in the smallest error. In addition, we optimize the performance of an algorithm that generates a continuous grid-less turbulent magnetic field by more than an order of magnitude. This continuous method becomes practicable for the simulation of large particle numbers and its accuracy is only limited by the used number of wave-modes. We show that by using more than 100 wave-modes the diffusive behavior of the spatial particle distribution in form of the diffusion coefficient is determined with an error less than a few percentage.
△ Less
Submitted 12 July, 2021; v1 submitted 23 July, 2019;
originally announced July 2019.