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…

April 28, 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.

September 30, 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; }

September 20, 2021