Aproximación del campo gravitatorio creado por un sistema de masas puntuales con Chombo

Para empezar, vamos a calcular numéricamente, y en mecánica clásica,  el campo creado por una partícula con masa en el espacio. Si colocamos las fronteras de nuestro dominio $latex Omega in mathbb{R}^3$ lo suficientemente lejos, el cálculo del potencial gravitatorio se reduce a resolver la ecuación de Poisson $latex Delta u(x,y,z) = 4 pi G …

Experiencias elípticas con Chombo (II)

En el libro de Análisis numérico de Burden y Faires aparecen una serie de ejemplos y ejercicios resueltos de PDEs elípticas en 2D. Vamos a intentar resolverlos utilizando Chombo… La primera ecuación corresponde con un ejemplo y es $latex frac{partial^2}{partial x^2}u(x,y) + frac{partial^2}{partial y^2}u(x,y) = 0$ en $latex Omega = { (x,y): 0<x<0.5, 0<y<0.5}$ y …

Experiencias elípticas con Chombo

El software Chombo, del Berkeley Lab, combina los métodos en diferencias finitas con las mallas adaptativas (AMR) para resolver, entre otras, PDEs elípticas. Las siguientes imágenes, en 2D y 3D, se han obtenido a partir de su AMRPoisson: