CoCoNuT es un código que permite realizar simulaciones de colapso estelar. Reescribimos las ecuaciones CFC, que son un caso particular de la aproximación FCF haciendo que las $latex h^{ij}$ sean cero, en terminos de las variables que éste utiliza. Empezamos con una auxilar: $latex Delta X^i = 8 pi f^{ij}S_j^* – frac{1}{3}mathcal{D}^i mathcal{D}_j X^j$ donde: …
Archivos mensuales:febrero 2013
Diagnostic Tools (from A Multigrid Tutorial)
Textualmente del libro: «As with any numerical code, debugging can be the most difficult part of creating a successful program. For multigrid, this situation is exacerbated in two ways. First, the interaction between the various multigrid components are very subtle, and it can be difficult to determine which part of a code is defective. Even …
Seguir leyendo «Diagnostic Tools (from A Multigrid Tutorial)»
El lado bueno de las cosas
No tiene nada que ver con los temas habituales del blog, pero bueno, la vida es algo mas que ciencia y tecnología ¿no?. Simplemente unas frases que me gustaron de «El lado bueno de las cosas»: «Estoy convencido de ello. Tienes que hacer todo lo que puedes y esforzarte al máximo y, si mantienes el …
La derivación y la integración desde un punto de vista algebraico. Utilización en la teoría cuántica de campos.
Simplemente un reblog de nuestro crack Terry Tao sobre como capturar los conceptos esenciales de la derivación y la integración de manera algebraica para permitir su utilización sobre otros sistemas numéricos distintos de aquellos que soportan el concepto de límite, o sea, los reales y los complejos. Posteriormente comenta como puede utilizar éstas cuando trabaja …
PDE de tipo Poisson con solución analítica
Una manera sencilla de tener una ecuación de Poisson en $latex 3D$ de la que conocer su solución analítica es la siguiente. Para empezar, consideramos una función: $latex u(x,y,z)$ a la que le aplicamos el operador $latex Delta$ y obtendremos otra función: $latex s(x,y,z)$. Ya tenemos $latex Delta u = s$, es decir, $latex frac{partial^2}{partial …
Diferencias finitas en sistemas no rectangulares para ecuaciones elípticas.
Aunque siempre podemos hacer cambios de coordenadas, vamos a ver como quedan los esquemas de diferencias finitas en sistemas no rectangulares: coordenadas cilíndricas, $latex (rho,phi, z)$, y coordenadas esféricas, $latex (r,theta,phi)$. Nos centraremos en la ecuación de Poisson aunque la técnica se puede extender de manera inmediata a cualquier tipo de PDE. En coordenadas cilíndricas …
Seguir leyendo «Diferencias finitas en sistemas no rectangulares para ecuaciones elípticas.»
Operador Laplaciano n-dimensional. Discretización y fronteras mediante tensores.
En $latex n$ dimensiones, el operador Laplaciano queda como: $latex Delta u= sum_{i=1}^n frac{partial^2}{partial x_i^2}u$ en coordenadas cartesianas, y como: $latex Delta u = frac{partial}{partial r^2}u + frac{n-1}{r}frac{partial}{partial r}u + frac{1}{r^2}Delta_{S^{n-1}}u$ en esféricas, donde $latex Delta_{S^{n-1}}$ es el operador de Laplace-Beltrami, una generalización del Laplaciano para funciones definidas sobre variedades, en la $latex (n-1)$-esfera ($latex …
Seguir leyendo «Operador Laplaciano n-dimensional. Discretización y fronteras mediante tensores.»
¿Qué pasa con tres o todas las fronteras Neumann en 2D?
Vamos a suponer $latex n=3$ para reducir el tamaño de las matrices. Empezamos suponiendo que conocemos: $latex frac{partial}{partial x}|_{0,0,}u, frac{partial}{partial x}|_{0,1}u, frac{partial}{partial x}|_{0,2}u$ $latex frac{partial}{partial y}|_{0,0}u, frac{partial}{partial y}|_{1,0}u$ $latex frac{partial}{partial y}|_{0,2}u, frac{partial}{partial y}|_{1,2}u$ $latex u|_{2,0}, u|_{2,1}, u|_{2,2}$ Discretizamos: $latex frac{u_{-1,0}-2u_{0,0}+u_{1,0}}{h^2} + frac{u_{0,-1}-2u_{0,0}+u_{0,1}}{h^2} = f_{0,0}$ $latex frac{u_{-1,1}-2u_{0,1}+u_{1,1}}{h^2} + frac{u_{0,0}-2u_{0,1}+u_{0,2}}{h^2} = f_{0,1}$ $latex frac{u_{-1,2}-2u_{0,2}+u_{1,2}}{h^2} + frac{u_{0,1}-2u_{0,2}+u_{0,3}}{h^2} = …
Seguir leyendo «¿Qué pasa con tres o todas las fronteras Neumann en 2D?»
Discretización de PDEs, matrices por bloques, simetrizaciones y rangos.
Suponemos $latex Delta u = f$ en $latex 2D$, es decir, $latex frac{partial^2}{partial x^2}u(x,y) + frac{partial^2}{partial y^2}u(x,y) = f(x,y)$. Miraremos como queda la matriz del sistema al discretizar, como simetrizarla y su rango en tres casos: condición Neuman respecto $latex x$ en una frontera, condición Neumann respecto $latex y$ en una frontera y condición Neumann …
Seguir leyendo «Discretización de PDEs, matrices por bloques, simetrizaciones y rangos.»
Ejemplo sencillo de frontera Neumann en 2D
Suponemos $latex n=5$. En el caso de tener todas las fronteras con condiciones Dirichlet: $latex frac{u_{0,1} -2u_{1,1} + u_{2,1}}{h^2} + frac{u_{1,0} -2u_{1,1} + u_{1,2}}{h^2} = f_{1,1}$ para $latex i,j=1,1$, $latex frac{u_{0,2} -2u_{1,2} + u_{2,2}}{h^2} + frac{u_{1,1} -2u_{1,2} + u_{1,3}}{h^2} = f_{1,2}$ para $latex i,j=1,2$, $latex frac{u_{0,3} -2u_{1,3} + u_{2,3}}{h^2} + frac{u_{1,2} -2u_{1,3} + u_{1,4}}{h^2} = …
Condiciones de contorno tipo Neumann
En el post anterior hablamos sobre condiciones de frontera y su transferencia entre mallas pero no comentamos en el caso de que las condición haga referencia al valor de la derivada y no al de la función: condición de Neumann. En $latex 1D$ supongamos que ahora tenemos $latex frac{partial^2}{partial x^2}u = f$ en $latex [a,b]$ …
Vaya con las fronteras…
A ver si nos aclaramos sobre como van las condiciones frontera en las transiciones entre mallas… Empezamos en $latex 1D$. Tenemos $latex Delta u = f$, o lo que es lo mismo, $latex u_{xx} = f$ (en este caso es una ODE pero bueno…) definida en $latex [a,b]$ con $latex u(a) = u_a$ y $latex …