Mates

Schwarzschild - Absorción de ondas escalares

Agujero negro de SchwarzschildAbsorción de ondas escalares - Primeros pasos Parte 1 - Coordenadas tortuga Voy a hacer una serie de posts hablando de cómo desarrollé mi TFM, y éste es el primero. El objetivo del TFM era poder distinguir entre diferentes ECOs (Extreme / Exotic Compact Objects), como un agujero negro o un agujero de gusano. Los papers clave para el desarrollo fueron Absorption by black hole remnants in metric-affine gravity y Reissner-Nordström black holes in extended Palatini theories, aunque hay muchos más…

Leer más →

28 de abril de 2022

Coordenadas tortuga y GSL

Coordenadas tortuga y GSLPor cuestiones que aún no puedo desvelar (misterio misterioso) estoy necesitando obtener las coordenadas tortuga, r(r *), de forma numérica. Por lo visto para el estudio de Bounce Black holes y Wormholes es relativamente fácil obtener r *( r) de forma analítica pero muy difícil obtener r(r *). Para entrenar un poco antes de ir al ring podemos estudiar el caso de Schwarzschild. SchwarzschildEn el caso de Schwarzschild sabemos que r *( r) = (1 - rS/r)^-1 y por tanto r(r *) = (1 - rS/r), o sea, que tenemos ambas expresiones analíticas lo que es perfecto para probar nuestros métodos.

Leer más →

30 de septiembre de 2021

GSL

GSL (GNU Scientific Library) Hace ya un tiempo que estoy utilizando la librería GSL y la verdad que me ha parecido una librería fantástica con una buena documentación. Aún así, sobre las ODE voy a dar un par de apuntes. Código algo más legibleVamos a tomar el mismo ejemplo que en la documentación, un oscilador de Van der Pol con la ecuación Ecuación del oscilador de Van der Pol y el sistema Sistema Van del Pol Para hacer el código más legible lo haremos más cercano a nuestras ecuaciones int func (double t, const double y[], double yp[], void* params) { (void)(t); /* Evitar warning no usado */ double mu = *(double*)params; double u = y[0]; double v = y[1]; double ut = v; double vt = -u + mu*v*(1 - u*u); yp[0] = ut; yp[1] = vt; return GSL_SUCCESS; }

Leer más →

20 de septiembre de 2021