Orbits in the Field of a Gravitating Magnetic Monopole

Orbits of test particles and light rays are an important tool to study the properties of space-time metrics. Here we systematically study the properties of the gravitational field of a globally regular magnetic monopole in terms of the geodesics of test particles and light. The gravitational field depends on two dimensionless parameters, defined as ratios of the characteristic mass scales present. For critical values of these parameters the resulting metric coefficients develop a singular behavior, which has profound influence on the properties of the resulting space-time and which is clearly reflected in the orbits of the test particles and light rays.Comment: 24 pages, 15 figures. Accepted for publication in GR

Meshfree finite differences for vector Poisson and pressure Poisson equations with electric boundary conditions

We demonstrate how meshfree finite difference methods can be applied to solve vector Poisson problems with electric boundary conditions. In these, the tangential velocity and the incompressibility of the vector field are prescribed at the boundary. Even on irregular domains with only convex corners, canonical nodal-based finite elements may converge to the wrong solution due to a version of the Babuska paradox. In turn, straightforward meshfree finite differences converge to the true solution, and even high-order accuracy can be achieved in a simple fashion. The methodology is then extended to a specific pressure Poisson equation reformulation of the Navier-Stokes equations that possesses the same type of boundary conditions. The resulting numerical approach is second order accurate and allows for a simple switching between an explicit and implicit treatment of the viscosity terms.Comment: 19 pages, 7 figure

Adaptive step ODE algorithms for the 3D simulation of electric heart activity with graphics processing units

In this paper we studied the implementation and performance of adaptive step methods for large systems of ordinary differential equations systems in graphics processing units, focusing on the simulation of three-dimensional electric cardiac activity. The Rush-Larsen method was applied in all the implemented solvers to improve efficiency. We compared the adaptive methods with the fixed step methods, and we found that the fixed step methods can be faster while the adaptive step methods are better in terms of accuracy and robustness. (c) 2013 Elsevier Ltd. All rights reserved.This work has been partially funded by Universitat Politecnica de Valencia through Programa de Apoyo a la InvestigaciOn y Desarrollo (PAID-06-11) and (PAID-05-12), by Generalitat Valenciana through projects PROMETEO/2009/013 and Ayudas para la realizacion de proyectos de I+D para grupos de investigacion emergentes GV/2012/039, and by Ministerio Espafiol de Economia y Competitividad through project TEC2012-38142-004.García Mollá, VM.; Liberos Mascarell, A.; Vidal Maciá, AM.; Guillem Sánchez, MS.; Millet Roig, J.; González Salvador, A.; Martínez Zaldívar, FJ.... (2014). Adaptive step ODE algorithms for the 3D simulation of electric heart activity with graphics processing units. Computers in Biology and Medicine. 44:15-26. https://doi.org/10.1016/j.compbiomed.2013.10.023S15264

Indirect measurements of neutron-induced reaction cross sections at storage rings

Neutron-induced reaction cross sections of unstable nuclei are essential for understanding the synthesis of heavy elements in stars. However, their measurement is very difficult due to the radioactivity of the targets involved. We propose to circumvent this problem by using for the first time the surrogate reaction method in inverse kinematics at heavy-ion storage rings. In this contribution, we describe the developments we have done to perform surrogate-reaction studies at the storage rings of GSI/FAIR. In particular, we present the first results of the proof of principle experiment, which we conducted recently at the Experimental Storage Ring (ESR)