You are currently browsing the monthly archive for febrero 2014.

La salida ahora para un tensor dos veces contravariante en la base ortonormal queda:

CovDerTen2BiSphCom1,

Para primera ecuación:

\boxed{\Delta \Theta_\beta = \frac{3}{2} \mathcal{D}_i \mathcal{D}_j (\alpha \psi^{-6} \hat{A}^{ij}) },

definimos como antes

V^i := \mathcal{D}_j \alpha \psi^{-6} \hat{A}^{ij},

de manera que la ecuación original la reescribimos como

\Delta \Theta_\beta = \frac{3}{2} \mathcal{D}_i V^i,

De esta manera, en nuestras coordenadas obtenemos:

\Delta \Theta_\beta = \frac{3}{2} \mathcal{D}_i V^i = \frac{3}{2} (\mathcal{D}_{\xi} V^{\xi} + \mathcal{D}_{\bar{\eta}} V^{\bar{\eta}} + \mathcal{D}_{\varphi} V^{\varphi}) =

div_biSphComNor1

donde

V^{\xi} = \mathcal{D}_{\xi} (\alpha \psi^{-6} \hat{A}^{\xi \xi}) + \mathcal{D}_{\bar{\eta}} ( \alpha \psi^{-6} \hat{A}^{\xi \bar{\eta}} ) + \mathcal{D}_{\varphi} ( \alpha \psi^{-6} \hat{A}^{\xi \varphi} ),

V^{\bar{\eta}} = \mathcal{D}_{\xi} (\alpha \psi^{-6} \hat{A}^{\bar{\eta} \xi}) + \mathcal{D}_{\bar{\eta}} ( \alpha \psi^{-6} \hat{A}^{\bar{\eta} \bar{\eta}} ) + \mathcal{D}_{\varphi} ( \alpha \psi^{-6} \hat{A}^{\bar{\eta} \varphi} ),

V^{\varphi} = \mathcal{D}_{\xi} (\alpha \psi^{-6} \hat{A}^{\varphi \xi}) + \mathcal{D}_{\bar{\eta}} ( \alpha \psi^{-6} \hat{A}^{\varphi \bar{\eta}} ) + \mathcal{D}_{\varphi} ( \alpha \psi^{-6} \hat{A}^{\varphi \varphi} ),

que desarrollando las covariantes quedan:

V^{\xi} = \frac{(1 - \bar{r})^2}{a} \partial_{\bar{r}} (\alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}}) +

+ \frac{1 - \bar{r}}{a \bar{r}} [ \partial_{\theta} ( \alpha \psi^{-6} \hat{A}^{\bar{r} \theta}) + \alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}} - \alpha \psi^{-6} \hat{A}^{\theta \theta} ] +

+ \frac{1 - \bar{r}}{a \bar{r}} [ \csc \theta \partial_{\varphi} ( \alpha \psi^{-6} \hat{A}^{\bar{r} \varphi}) + \alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}} + \cot \theta \alpha \psi^{-6} \hat{A}^{\bar{r} \theta} - \alpha \psi^{-6} \hat{A}^{\varphi \varphi}] ),

V^{\bar{\eta}} = \frac{(1 - \bar{r})^2}{a} \partial_{\bar{r}} (\alpha \psi^{-6} \hat{A}^{\theta \bar{r}}) +

+ \frac{1 - \bar{r}}{a \bar{r}} [ \partial_{\theta} ( \alpha \psi^{-6} \hat{A}^{\theta \theta}) + 2 \alpha \psi^{-6} \hat{A}^{\bar{r} \theta} ] +

+ \frac{1 - \bar{r}}{a \bar{r}} [ \csc \theta \partial_{\varphi} ( \alpha \psi^{-6} \hat{A}^{\theta \varphi}) + \alpha \psi^{-6} \hat{A}^{\bar{r} \theta} + \cot \theta \alpha \psi^{-6} \hat{A}^{\theta \theta} - \cot \theta \alpha \psi^{-6} \hat{A}^{\varphi \varphi} ] ),

V^{\varphi} = \frac{(1 - \bar{r})^2}{a} \partial_{\bar{r}} (\alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}}) +

+ \frac{1 - \bar{r}}{a \bar{r}} [ \partial_{\theta} ( \alpha \psi^{-6} \hat{A}^{\bar{r} \theta}) + \alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}} - \alpha \psi^{-6} \hat{A}^{\theta \theta} ] +

+ \frac{1 - \bar{r}}{a \bar{r}} [ \csc \theta \partial_{\varphi} ( \alpha \psi^{-6} \hat{A}^{\bar{r} \varphi}) + \alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}} + \cot \theta \alpha \psi^{-6} \hat{A}^{\bar{r} \theta} - \alpha \psi^{-6} \hat{A}^{\varphi \varphi}] ),

que combinandolo con la anterior, queda:

Finalmente, las ecuaciones:

\boxed{\Delta \beta^i = 2\mathcal{D}_j ( \alpha \psi^{-6} \hat{A}^{ij} ) - \frac{1}{3} \mathcal{D}^i \Theta_{\beta} },

con las que procederemos de manera similar a como hemos hecho con las X^i, es decir, calculando las fuentes en una base, haciendo un cambio de base que las desacople (cartesianas), resolviendolas de manera independiente y volviendo a la base original:

S^i_\beta (\bar{r},\theta,\varphi) := 2\mathcal{D}_j ( \alpha \psi^{-6} \hat{A}^{ij} ) - \frac{1}{3} \mathcal{D}_i \Theta_{\beta},

que quedan:

S^{\xi}_\beta= 2 \big [ \mathcal{D}_{\xi} (\alpha \psi^{-6} \hat{A}^{\xi \xi}) + \mathcal{D}_{\bar{\eta}} (\alpha \psi^{-6} \hat{A}^{\xi \bar{\eta}}) + \mathcal{D}_{\varphi} (\alpha \psi^{-6} \hat{A}^{\xi \varphi}) \big ] - \frac{1}{3} \mathcal{D}_{\xi} \Theta_\beta,

S^{\bar{\eta}}_\beta = 2 \big [ \mathcal{D}_{\xi} (\alpha \psi^{-6} \hat{A}^{\bar{\eta} \xi}) + \mathcal{D}_{\bar{\eta}} (\alpha \psi^{-6} \hat{A}^{\bar{\eta} \bar{\eta}}) + \mathcal{D}_{\varphi} (\alpha \psi^{-6} \hat{A}^{\bar{\eta} \varphi}) \big ] - \frac{1}{3} \mathcal{D}_{\bar{\eta}} \Theta_\beta,

S^{\varphi}_\beta = 2 \big [ \mathcal{D}_{\xi} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{x}}) + \mathcal{D}_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{y}}) + \mathcal{D}_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{z}}) \big ] - \frac{1}{3} \mathcal{D}_{\varphi} \Theta_\beta,

donde las derivadas covariantes del tensor dos veces contravariante:

T^{ij}:=\alpha \psi^{-6} \hat{A}^{ij}

son como acabamos de hacer en la ecuación anterior y las del escalar \Theta_\beta es como ya hicimos con las X^i:

S^{\xi} = 2 V^{\xi} - \frac{\mbox{\scriptsize cosh} \frac{b \bar{\eta}}{1 - \bar{\eta}} - \cos \xi}{3a} \partial_{\xi} \Theta_{\beta},

S^{\bar{\eta}} = 2 V^{\bar{\eta}} - \frac{\mbox{\scriptsize cosh} \frac{b \bar{\eta}}{1 - \bar{\eta}} - \cos \xi}{3a} \frac{(\bar{\eta} - 1)^2}{b} \partial_{\bar{\eta}} \Theta_{\beta},

S^{\varphi} = 2 V^{\varphi} - \frac{\mbox{\scriptsize cosh} \frac{b \bar{\eta}}{1 - \bar{\eta}} - \cos \xi}{3a} \csc \xi \partial_{\varphi} \Theta_{\beta}.

Hacemos a continuación el cambio:

[S^{\xi}(\xi,\bar{\eta},\varphi),S^{\bar{\eta}}(\xi,\bar{\eta},\varphi),S^{\varphi}(\xi,\bar{\eta},\varphi)] \rightarrow

\rightarrow [S^x(\xi,\bar{\eta},\varphi), S^y(\xi,\bar{\eta},\varphi), S^z(\xi,\bar{\eta},\varphi)],

y resolvemos:

\Delta \beta^{x} = S^{x}

\Delta \beta^{y} = S^{y}

\Delta \beta^{z} = S^{z},

deshaciendo el cambio:

[\beta^x(\xi,\bar{\eta},\varphi), \beta^y(\xi,\bar{\eta},\varphi), \beta^z(\xi,\bar{\eta},\varphi)] \rightarrow

\rightarrow [\beta^{\xi}(\xi,\bar{\eta},\varphi),\beta^{\bar{\eta}}(\xi,\bar{\eta},\varphi),\beta^{\varphi}(\xi,\bar{\eta},\varphi)]

para terminar.

Anuncios

Recordemos lo ya expuesto en este post: que en las coordenadas biesféricas (\xi, \eta, \varphi), las dos primeras (\xi, \eta) provienen de las coordenadas bipolares, donde la primera indica el ángulo entre las dos rectas que unen nuestro punto con los dos focos que necesitamos para determinar las bipolares y la segundo es el logartimo del ratio entre la longitud de estas dos rectas, mientras que la última proviene de rotarlas alrededor del eje que une los focos.

Compactificamos la segunda coordenada mediante \boxed{\eta = \frac{b \bar{\eta}}{1 - \bar{\eta}}}.

El Laplaciano, en estas coordenadas y con esta compactificación, queda:

\Delta = \frac{(\cos \xi - \mbox{\scriptsize cosh} \frac{b \bar{\eta}}{1-\bar{\eta}})^2}{a^2} \big [ \partial_{\xi \xi} + \csc \xi \frac{-1 + \cos \xi \, \mbox{\scriptsize cosh} \frac{b \bar{\eta}}{1-\bar{\eta}}}{\mbox{\scriptsize cosh} \frac{b \bar{\eta}}{1-\bar{\eta}} - \cos \xi} \partial_{\xi}

+\frac{(\bar{\eta} - 1)^4}{b^2} \partial_{\bar{\eta} \bar{\eta}} + \frac{(\bar{\eta} - 1)^2}{b^2} (2(\bar{\eta}-1) -\frac{b \, \mbox{\scriptsize sinh} \frac{b \bar{\eta}}{1 - \bar{\eta}}}{\mbox{\scriptsize cosh} \frac{b \bar{\eta}}{1 - \bar{\eta}} - \cos \xi}) \partial_{\bar{\eta}} + \csc^2 \xi \partial_{\varphi} \big ],

las derivadas covariantes de covectores (1-formas):

CovDer_BiSphComNor1

y las fuentes:

\boxed{\Delta \Theta_{X} = 6 \pi \mathcal{D}^j S^*_j}

\Delta \Theta_X = 6 \pi f^{ji} \mathcal{D}_i S^*_j = 6 \pi ( \mathcal{D}_{\xi} S^*_{\xi} + \mathcal{D}_{\bar{\eta}} S^*_{\bar{\eta}} + \mathcal{D}_{\varphi} S^*_{\varphi} ) =

s1_biSphComNor1

\boxed{\Delta X^{i} = 8 \pi f^{ij} S^*_j - \frac{1}{3} \mathcal{D}^i \Theta_X}

Pasando la derivada contravariante a covariante mediante la métrica, queda:

\Delta X^{i} = 8 \pi f^{ij} S^*_j - \frac{1}{3} f^{ik} \mathcal{D}_k \Theta_X.

Definimos ahora

S_X^i := 8 \pi f^{ij} S^*_j - \frac{1}{3} f^{ik} \mathcal{D}_k \Theta_X,

de manera que:

S_X^{\xi} = 8 \pi f^{\xi j} S^*_j - \frac{1}{3} f^{\xi k}\mathcal{D}_{k} \Theta_X = 8 \pi S^*_{\xi} - \frac{1}{3} \mathcal{D}_{\xi} \Theta_X =

= 8 \pi S^*_{\xi} - \frac{\mbox{\scriptsize cosh} \frac{b \bar{\eta}}{1 - \bar{\eta}} - \cos \xi}{a} \partial_{\xi} \Theta_X

S_X^{\bar{\eta}} = 8 \pi f^{\bar{\eta} j} S^*_j - \frac{1}{3} f^{\bar{\eta} k} \mathcal{D}_{k} \Theta_X = 8 \pi S^*_{\bar{\eta}} - \frac{1}{3} \mathcal{D}_{\bar{\eta}} \Theta_X =

= 8 \pi S^*_{\bar{\eta}} - \frac{\mbox{\scriptsize cosh} \frac{b \bar{\eta}}{1 - \bar{\eta}} - \cos \xi}{a} \frac{(\bar{\eta} - 1)^2}{b} \partial_{\bar{\eta}} \Theta_X

S_X^{\varphi} = 8 \pi f^{\varphi j} S^*_j - \frac{1}{3} f^{\varphi k} \mathcal{D}^{k} \Theta_X = 8 \pi S^*_{\varphi} - \frac{1}{3} \mathcal{D}_{\varphi} \Theta_X =

= 8 \pi S^*_{\varphi} - \frac{\mbox{\scriptsize cosh} \frac{b \bar{\eta}}{1 - \bar{\eta}} - \cos \xi}{a} \csc \xi \partial_{\varphi} \Theta_X

En este punto tenemos que el vector

(S_X^{\xi}(\xi,\bar{\eta},\varphi),S_X^{\bar{\eta}}(\xi,\bar{\eta},\varphi),S_X^{\varphi}(\xi,\bar{\eta},\varphi))

expresado en la base que resulta de normalizar la base coordenada \{ \partial_{\xi}, \partial_{\bar{\eta}}, \partial_{\varphi} \}. Lo que hacemos ahora es expresar este vector en la nueva base \{ \partial_x, \partial_y, \partial_z \}, de manera que obtenemos

(S_X^{x}(\xi,\bar{\eta},\varphi),S_X^{y}(\xi,\bar{\eta},\varphi),S_X^{z}(\xi,\bar{\eta},\varphi)).

y como es esta base las ecuaciones están desacopladas y \Theta_X es un campo escalar, resolvemos independientemente:

\Delta X^{x} = S_X^{x},

\Delta X^{y} = S_X^{y},

\Delta X^{z} = S_X^{z}.

Finalmente, con el cambio de base inverso, calculamos a partir de (X^{x},X^{y},X^{z}) el vector (X^{\xi},X^{\bar{\eta}},X^\varphi) .

\underline{\hat{A}^{ij} = \mathcal{D}^i X^j + \mathcal{D}^j X^i - \frac{2}{3} \mathcal{D}_k X^k f^{ij}}

Necesitamos ahora la derivada covariante de un vector (hasta ahora habían coincidido las derivadas covariantes de vectores y covectores, pero en este caso no):

CovDer_BiSphComNor1_vec

volvemos a pasar las derivadas contravariantes a covariantes:

\hat{A}^{ij} = f^{im} \mathcal{D}_m X^j + f^{jn} \mathcal{D}_n X^i - \frac{2}{3} f^{ij} \mathcal{D}_k X^{k}

y obtenemos:

\hat{A}^{\xi \xi} = f^{\xi m} \mathcal{D}_m X^{\xi} + f^{\xi n} \mathcal{D}_n X^{\xi} - \frac{2}{3} \mathcal{D}_k X^{k} = \frac{2}{3}( 2 \mathcal{D}_{\xi} X^{\xi} - \mathcal{D}_{\bar{\eta}} X^{\bar{\eta}} - \mathcal{D}_{\varphi} X^{\varphi}) =

A11_biSphComNor1

\hat{A}^{\xi \bar{\eta}} = f^{\xi m} \mathcal{D}_m X^{\bar{\eta}} + f^{\bar{\eta} n} \mathcal{D}_n X^{\xi} = \mathcal{D}_{\xi} X^{\bar{\eta}} + \mathcal{D}_{\bar{\eta}} X^{\xi} =

A12_biSphComNor1

\hat{A}^{\xi \varphi} = f^{\xi m} \mathcal{D}_m X^{\varphi} + f^{\varphi n} \mathcal{D}_n X^{\xi} = \mathcal{D}_{\xi} X^{\varphi} + \mathcal{D}_{\varphi} X^{\xi} =

A13_biSphComNor1

\hat{A}^{\bar{\eta} \bar{\eta}} = f^{\bar{\eta} m} \mathcal{D}_m X^{\bar{\eta}} + f^{\bar{\eta} n} \mathcal{D}_n X^{\bar{\eta}} - \frac{2}{3} \mathcal{D}_k X^{k} = \frac{2}{3}( - \mathcal{D}_{\xi} X^{\xi} + 2 \mathcal{D}_{\bar{\eta}} X^{\bar{\eta}} - \mathcal{D}_{\varphi} X^{\varphi}) =

A22_biSphComNor1

\hat{A}^{\bar{\eta} \varphi} = f^{\bar{\eta} m} \mathcal{D}_m X^{\varphi} + f^{\varphi n} \mathcal{D}_n X^{\bar{\eta}} = \mathcal{D}_{\bar{\eta}} X^{\varphi} + \mathcal{D}_{\varphi} X^{\bar{\eta}} =

A23_biSphComNor1

\hat{A}^{\varphi \varphi} = f^{\varphi m} \mathcal{D}_m X^{\varphi} + f^{\varphi n} \mathcal{D}_n X^{\varphi} - \frac{2}{3} \mathcal{D}_k X^{k} = \frac{2}{3}( - \mathcal{D}_{\bar{r}} X^{\bar{r}} - \mathcal{D}_{\theta} X^{\theta} +2 \mathcal{D}_{\varphi} X^{\varphi}) =

A33_biSphComNor1

Las dos ecuaciones no lineales correspondientes al factor conforme \psi y al lapse \alpha, como no contienen derivadas covariantes, quedan como las teniamos:

\boxed{\Delta \psi = -2 \pi E^* \psi^{-1} - \frac{1}{8}(f_{il} f_{jm} \hat{A}^{lm} \hat{A}^{ij}) \psi^{-7} }

\boxed{\Delta (\alpha \psi) = [ 2 \pi (E^* + 2 S^*) \psi^{-7} + \frac{1}{8}(f_{il} f_{jm} \hat{A}^{lm} \hat{A}^{ij}) \psi^{-8} ] (\alpha \psi) }

Finalmente, para el shift \beta y su ecuación auxiliar tenemos:

\boxed{\Delta \Theta_{\beta} = \frac{3}{4} \mathcal{D}_i \mathcal{D}_j (2 \alpha \psi^{-6} \hat{A}^{ij} )}

\boxed{\Delta \beta^i = \mathcal{D}_j ( 2 \alpha \psi^{-6} \hat{A}^{ij} ) - \frac{1}{3} \mathcal{D}^i \Theta_{\beta} }

que trataremos en el siguiente post.

El círculo de Apolonio de focos F_1 y F_2 de razón r es el lugar geométrico de los puntos de plano P tales que:

\overline{PF_1}/\overline{PF_2} = r.

Definimos las coordenadas bipolares (\xi,\eta) de la siguiente manera:

x = a \frac{\mbox{\scriptsize sinh} \eta}{\mbox{\scriptsize cosh} \eta - cos \xi},

y = a \frac{\sin \eta}{\mbox{\scriptsize cosh} \eta - cos \xi},

donde en los puntos F_1=(-a,0) y F_2=(a,0) tenemos definidos los dos focos.

La interpretación geométrica es la siguiente: dado un punto P=(x,y) del plano, la coordenada \xi es el ángulo que forman los vectores u:=\overline{PF_1} y v:=\overline{PF_2}:

\xi = \arccos \frac{u \cdot v}{|u||v|}, de manera que \xi \in [-\pi,\pi],

mientras que la coordenada \eta es el rátio entre los módulos de éstos:

\eta = \log \frac{|u|}{|v|}, por lo que \eta \in ]0,\infty[,

que si las expresamos de manera que aparezcan explicitamente las coordenadas cartesianas, quedan:

\xi = \arccos \frac{x^2 + y^2 - a^2}{\sqrt{(a-x)^2 + y^2}\sqrt{(a+x)^2+y^2}},

\eta = \frac{\log[(a+x)^2+y^2] - \log[(a-x)^2+y^2]}{2}.

A continuación, una imagen con lineas coordenadas:

BiSphCoo

Para pasar a coordenadas en tres dimensiones, lo único que hacemos es añadir una nueva coordenada \varphi que nos indica un ángulo de rotación respecto de un eje. Si el eje es la recta que une los focos, obtenemos las coordenadas biesféricas, que son las que nos interesarán y utilizaremos en posteriores entradas. Si lo que hacemos es rotar respecto a la recta perpendicular a la anterior, la que separa los dos focos, obtenemos las coordenadas toroidales.

La salida ahora para un tensor dos veces contravariante en la base ortonormal queda:

CovDerTen2SphCom1,

Para primera ecuación:

\boxed{\Delta \Theta_\beta = \frac{3}{2} \mathcal{D}_i \mathcal{D}_j (\alpha \psi^{-6} \hat{A}^{ij}) },

definimos como antes

V^i := \mathcal{D}_j \alpha \psi^{-6} \hat{A}^{ij},

de manera que la ecuación original la reescribimos como

\Delta \Theta_\beta = \frac{3}{2} \mathcal{D}_i V^i,

De esta manera, en nuestras coordenadas obtenemos:

\Delta \Theta_\beta = \frac{3}{2} \mathcal{D}_i V^i = \frac{3}{2} (\mathcal{D}_{\bar{r}} V^{\bar{r}} + \mathcal{D}_{\theta} V^{\theta} + \mathcal{D}_{\varphi} V^{\varphi}) =

= \frac{3 - 3\bar{r} }{2a \bar{r}} [ (\bar{r}-\bar{r}^2) \partial_{\bar{r}} V^{\bar{r}} + 2 V^{\bar{r}} + \partial_{\theta} V^{\theta} + \cot \theta V^{\theta} + \csc \theta \partial_{\varphi} V^{\varphi} ],

donde

V^{\bar{r}} = \mathcal{D}_{\bar{r}} (\alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}}) + \mathcal{D}_{\theta} ( \alpha \psi^{-6} \hat{A}^{\bar{r} \theta} ) + \mathcal{D}_{\varphi} ( \alpha \psi^{-6} \hat{A}^{\bar{r} \varphi} ),

V^{\theta} = \mathcal{D}_{\bar{r}} (\alpha \psi^{-6} \hat{A}^{\theta \bar{r}}) + \mathcal{D}_{\theta} ( \alpha \psi^{-6} \hat{A}^{\theta \theta} ) + \mathcal{D}_{\varphi} ( \alpha \psi^{-6} \hat{A}^{\theta \varphi} ),

V^{\varphi} = \mathcal{D}_{\bar{r}} (\alpha \psi^{-6} \hat{A}^{\varphi \bar{r}}) + \mathcal{D}_{\theta} ( \alpha \psi^{-6} \hat{A}^{\varphi \theta} ) + \mathcal{D}_{\varphi} ( \alpha \psi^{-6} \hat{A}^{\varphi \varphi} ),

que desarrollando las covariantes quedan:

V^{\bar{r}} = \frac{(1 - \bar{r})^2}{a} \partial_{\bar{r}} (\alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}}) +

+ \frac{1 - \bar{r}}{a \bar{r}} [ \partial_{\theta} ( \alpha \psi^{-6} \hat{A}^{\bar{r} \theta}) + \alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}} - \alpha \psi^{-6} \hat{A}^{\theta \theta} ] +

+ \frac{1 - \bar{r}}{a \bar{r}} [ \csc \theta \partial_{\varphi} ( \alpha \psi^{-6} \hat{A}^{\bar{r} \varphi}) + \alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}} + \cot \theta \alpha \psi^{-6} \hat{A}^{\bar{r} \theta} - \alpha \psi^{-6} \hat{A}^{\varphi \varphi}] ),

V^{\theta} = \frac{(1 - \bar{r})^2}{a} \partial_{\bar{r}} (\alpha \psi^{-6} \hat{A}^{\theta \bar{r}}) +

+ \frac{1 - \bar{r}}{a \bar{r}} [ \partial_{\theta} ( \alpha \psi^{-6} \hat{A}^{\theta \theta}) + 2 \alpha \psi^{-6} \hat{A}^{\bar{r} \theta} ] +

+ \frac{1 - \bar{r}}{a \bar{r}} [ \csc \theta \partial_{\varphi} ( \alpha \psi^{-6} \hat{A}^{\theta \varphi}) + \alpha \psi^{-6} \hat{A}^{\bar{r} \theta} + \cot \theta \alpha \psi^{-6} \hat{A}^{\theta \theta} - \cot \theta \alpha \psi^{-6} \hat{A}^{\varphi \varphi} ] ),

V^{\varphi} = \frac{(1 - \bar{r})^2}{a} \partial_{\bar{r}} (\alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}}) +

+ \frac{1 - \bar{r}}{a \bar{r}} [ \partial_{\theta} ( \alpha \psi^{-6} \hat{A}^{\bar{r} \theta}) + \alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}} - \alpha \psi^{-6} \hat{A}^{\theta \theta} ] +

+ \frac{1 - \bar{r}}{a \bar{r}} [ \csc \theta \partial_{\varphi} ( \alpha \psi^{-6} \hat{A}^{\bar{r} \varphi}) + \alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}} + \cot \theta \alpha \psi^{-6} \hat{A}^{\bar{r} \theta} - \alpha \psi^{-6} \hat{A}^{\varphi \varphi}] ),

que combinandolo con la anterior, queda (solo escribimos como empezaría debido a la longitud de la ecuación):

\Delta \Theta_\beta = \frac{3 - 3\bar{r} }{2a \bar{r}} \Big [ (\bar{r}-\bar{r}^2) \partial_{\bar{r}} [\frac{(1 - \bar{r})^2}{a} \partial_{\bar{r}} (\alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}}) + \ldots ] + \ldots \Big ]

Finalmente, las ecuaciones:

\boxed{\Delta \beta^i = 2\mathcal{D}_j ( \alpha \psi^{-6} \hat{A}^{ij} ) - \frac{1}{3} \mathcal{D}^i \Theta_{\beta} },

con las que procederemos de manera similar a como hemos hecho con las X^i, es decir, calculando las fuentes en una base, haciendo un cambio de base que las desacople (cartesianas), resolviendolas de manera independiente y volviendo a la base original:

S^i_\beta (\bar{r},\theta,\varphi) := 2\mathcal{D}_j ( \alpha \psi^{-6} \hat{A}^{ij} ) - \frac{1}{3} \mathcal{D}_i \Theta_{\beta},

que quedan:

S^{\bar{r}}_\beta= 2 \big [ \mathcal{D}_{\bar{r}} (\alpha \psi^{-6} \hat{A}^{\bar{r} \bar{r}}) + \mathcal{D}_{\theta} (\alpha \psi^{-6} \hat{A}^{\bar{r} \theta}) + \mathcal{D}_{\varphi} (\alpha \psi^{-6} \hat{A}^{\bar{r} \varphi}) \big ] - \frac{1}{3} \mathcal{D}_{\bar{r}} \Theta_\beta,

S^{\theta}_\beta = 2 \big [ \mathcal{D}_{\bar{r}} (\alpha \psi^{-6} \hat{A}^{\theta \bar{r}}) + \mathcal{D}_{\theta} (\alpha \psi^{-6} \hat{A}^{\theta \theta}) + \mathcal{D}_{\varphi} (\alpha \psi^{-6} \hat{A}^{\theta \varphi}) \big ] - \frac{1}{3} \mathcal{D}_{\theta} \Theta_\beta,

S^{\varphi}_\beta = 2 \big [ \mathcal{D}_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{x}}) + \mathcal{D}_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{y}}) + \mathcal{D}_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{z}}) \big ] - \frac{1}{3} \mathcal{D}_{\bar{y}} \Theta_\beta,

donde las derivadas covariantes del tensor dos veces contravariante:

T^{ij}:=\alpha \psi^{-6} \hat{A}^{ij}

son como acabamos de hacer en la ecuación anterior y las del escalar \Theta_\beta es como ya hicimos con las X^i:

S^{\bar{r}} = 2 V^{\bar{r}} - \frac{(1-\bar{r})^2}{3a} \partial_{\bar{r}} \Theta_{\beta},

S^{\theta} = 2 V^{\theta} - \frac{1-\bar{r}}{3a\bar{r}} \partial_{\theta} \Theta_{\beta},

S^{\varphi} = 2 V^{\varphi} - \frac{1-\bar{r}}{3a\bar{r}} \csc \theta \partial_{\varphi} \Theta_{\beta}.

Hacemos a continuación el cambio:

[S^{\bar{r}}(\bar{r},\theta,\varphi),S^{\theta}(\bar{r},\theta,\varphi),S^{\varphi}(\bar{r},\theta,\varphi)] \rightarrow

\rightarrow [S^x(\bar{r},\theta,\varphi), S^y(\bar{r},\theta,\varphi), S^z(\bar{r},\theta,\varphi)],

y resolvemos:

\Delta \beta^{x} = S^{x}

\Delta \beta^{y} = S^{y}

\Delta \beta^{z} = S^{z},

deshaciendo el cambio:

[\beta^x(\bar{r},\theta,\varphi), \beta^y(\bar{r},\theta,\varphi), \beta^z(\bar{r},\theta,\varphi)] \rightarrow

\rightarrow [\beta^{\bar{r}}(\bar{r},\theta,\varphi),\beta^{\theta}(\bar{r},\theta,\varphi),\beta^{\varphi}(\bar{r},\theta,\varphi)]

para terminar.

La salida ahora para un tensor dos veces contravariante en la base ortonormal queda:

CovDerTen2CarCom2,

Para primera ecuación:

\boxed{\Delta \Theta_\beta = \frac{3}{2} \mathcal{D}_i \mathcal{D}_j (\alpha \psi^{-6} \hat{A}^{ij}) },

definimos como antes

V^i := \mathcal{D}_j \alpha \psi^{-6} \hat{A}^{ij},

de manera que la ecuación original la reescribimos como

\Delta \Theta_\beta = \frac{3}{2} \mathcal{D}_i V^i,

De esta manera, en nuestras coordenadas obtenemos:

\Delta \Theta_\beta = \frac{3}{2} \mathcal{D}_i V^i = \frac{3}{2} (\mathcal{D}_{\bar{x}} V^{\bar{x}} + \mathcal{D}_{\bar{y}} V^{\bar{y}} + \mathcal{D}_{\bar{z}} V^{\bar{z}}) =

= \frac{3}{2} (\frac{1 + \cos(\pi \bar{x})}{a \pi} \partial_{\bar{x}} V^{\bar{x}}+ \frac{1 + \cos(\pi \bar{y})}{b \pi} \partial_{\bar{y}} V^{\bar{y}} + \frac{1 + \cos(\pi \bar{z})}{c \pi} \partial_{\bar{z}} V^{\bar{z}}),

donde

V^{\bar{x}} = \mathcal{D}_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{x}}) + \mathcal{D}_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{x} \bar{y}} ) + \mathcal{D}_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{x} \bar{z}} ),

V^{\bar{y}} = \mathcal{D}_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{x}}) + \mathcal{D}_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{y} \bar{y}} ) + \mathcal{D}_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{y} \bar{z}} ),

V^{\bar{z}} = \mathcal{D}_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{x}}) + \mathcal{D}_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{z} \bar{y}} ) + \mathcal{D}_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{z} \bar{z}} ),

que desarrollando las covariantes quedan:

V^{\bar{x}} = \frac{1 + \cos(\pi \bar{x})}{a \pi} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{x}}) +

+ \frac{1 + \cos(\pi \bar{y})}{b \pi} \partial_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{x} \bar{y}} ) + \frac{1 + \cos(\pi \bar{z})}{c \pi} \partial_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{x} \bar{z}} ),

V^{\bar{y}} = \frac{1 + \cos(\pi \bar{x})}{a \pi} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{x}}) +

+ \frac{1 + \cos(\pi \bar{y})}{b \pi} \partial_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{y} \bar{y}} ) + \frac{1 + \cos(\pi \bar{z})}{c \pi} \partial_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{y} \bar{z}} ),

V^{\bar{z}} = \frac{1 + \cos(\pi \bar{x})}{a \pi} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{x}}) +

+ \frac{1 + \cos(\pi \bar{y})}{b \pi} \partial_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{z} \bar{y}} ) + \frac{1 + \cos(\pi \bar{z})}{c \pi} \partial_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{z} \bar{z}} ).

Por tanto, combiando todo, tenemos:

\Delta \Theta_\beta =

= \frac{3 + 3cos(\pi \bar{x})}{2 a \pi} \partial_{\bar{x}} \big [ \frac{1 + \cos(\pi \bar{x})}{a \pi} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{x}}) +

+ \frac{1 + \cos(\pi \bar{y})}{b \pi} \partial_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{x} \bar{y}} ) +

+ \frac{1 + \cos(\pi \bar{z})}{c \pi} \partial_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{x} \bar{z}} ) \big ] +

+ \frac{3 + 3 \cos (\pi \bar{y})}{2b \pi} \partial_{\bar{y}} \big [ \frac{1 + \cos(\pi \bar{x})}{a \pi} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{x}}) +

+ \frac{1 + \cos(\pi \bar{y})}{b \pi} \partial_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{y} \bar{y}} ) +

+ \frac{1 + \cos(\pi \bar{z})}{c \pi} \partial_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{y} \bar{z}} ) \big ] +

+ \frac{3 + 3 \cos(\pi \bar{z})}{2c \pi} \partial_{\bar{z}} \big [ \frac{1 + \cos(\pi \bar{x})}{a \pi} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{x}}) +

\frac{1 + \cos(\pi \bar{y})}{b \pi} \partial_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{z} \bar{y}} ) +

+ \frac{1 + \cos(\pi \bar{z})}{c \pi} \partial_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{z} \bar{z}} ) \big ],

Finalmente, las ecuaciones:

\boxed{\Delta \beta^i = 2\mathcal{D}_j ( \alpha \psi^{-6} \hat{A}^{ij} ) - \frac{1}{3} \mathcal{D}^i \Theta_{\beta} }.

son:

\Delta \beta^{\bar{x}} = 2 \big [ \mathcal{D}_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{x}}) + \mathcal{D}_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{y}}) + \mathcal{D}_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{z}}) \big ] - \frac{1}{3} \mathcal{D}_{\bar{x}} \Theta_\beta,

\Delta \beta^{\bar{y}} = 2 \big [ \mathcal{D}_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{x}}) + \mathcal{D}_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{y}}) + \mathcal{D}_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{z}}) \big ] - \frac{1}{3} \mathcal{D}_{\bar{y}} \Theta_\beta,

\Delta \beta^{\bar{z}} = 2 \big [ \mathcal{D}_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{x}}) + \mathcal{D}_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{y}}) + \mathcal{D}_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{z}}) \big ] - \frac{1}{3} \mathcal{D}_{\bar{z}} \Theta_\beta,

y al sustituir las derivadas covariantes:

\Delta \beta^{\bar{x}} = \frac{2 + 2 \cos(\pi \bar{x})}{a \pi} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{x}}) +

+ \frac{2 + 2 \cos(\pi \bar{y})}{b \pi} \partial_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{y}}) +

+ \frac{2 + 2 \cos(\pi \bar{z})}{c \pi} \partial_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{z}}) - \frac{1 + \cos(\pi \bar{x})}{3a \pi} \partial_{\bar{x}} \Theta_\beta,

\Delta \beta^{\bar{y}} = \frac{2 + 2 \cos(\pi \bar{x})}{a \pi} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{x}}) +

+ \frac{2 + 2 \cos(\pi \bar{y})}{b \pi} \partial_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{y}}) +

+ \frac{2 + 2 \cos(\pi \bar{z})}{c} \partial_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{z}}) - \frac{1 + \cos(\pi \bar{y})}{3b \pi} \partial_{\bar{y}} \Theta_\beta,

\Delta \beta^{\bar{z}} = \frac{2 + 2 \cos(\pi \bar{x})}{a \pi} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{x}}) +

+ \frac{2 + 2 \cos(\pi \bar{y})}{b \pi} \partial_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{y}}) +

+ \frac{2 + 2 \cos(\pi \bar{z})}{c \pi} \partial_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{z}}) - \frac{1 + \cos(\pi \bar{z})}{3c \pi} \partial_{\bar{z}} \Theta_\beta.

Copio a continuación la salida generada por nuestra función en Mathematica que nos calcula todas las derivadas covariantes de tensores con dos índices (aunque en este caso particular no es excesivamente laborioso, si lo es para el resto de entradas, por lo que evitaremos morir en el intento de pasarlas a latex 😉 ) . En particular, aquí lo hacemos para un vector dos veces contravariante y para la base ortonormal:

CovDerTen2CarCom1,

donde primer término corresponde al factor que acompaña a la derivada parcial y la matriz contiene los factores que acopañan a cada par de valores de los índices.

Vamos a ver ahora, en este caso, como quedan las ecuaciones del shift. Para la primera:

\boxed{\Delta \Theta_\beta = \frac{3}{2} \mathcal{D}_i \mathcal{D}_j (\alpha \psi^{-6} \hat{A}^{ij}) },

como contraemos el índice j quedando libre únicamente el i, definimos

V^i := \mathcal{D}_j \alpha \psi^{-6} \hat{A}^{ij},

de manera que la ecuación original la reescribimos como

\Delta \Theta_\beta = \frac{3}{2} \mathcal{D}_i V^i,

que nos ayudará a no liarnos, ya que ésta última queda como una derivada covariante de un vector donde éste, a su vez, lo calcularemos a parte como la derivada covariante de un tensor dos veces contravariante.

De esta manera, en nuestras coordenadas tenemos:

\Delta \Theta_\beta = \frac{3}{2} \mathcal{D}_i V^i = \frac{3}{2} (\mathcal{D}_{\bar{x}} V^{\bar{x}} + \mathcal{D}_{\bar{y}} V^{\bar{y}} + \mathcal{D}_{\bar{z}} V^{\bar{z}}) =

= \frac{3}{2} (\frac{|\bar{x}^2-1|}{a} \partial_{\bar{x}} V^{\bar{x}}+ \frac{|\bar{y}^2-1|}{b} \partial_{\bar{y}} V^{\bar{y}} + \frac{|\bar{z}^2-1|}{c} \partial_{\bar{z}} V^{\bar{z}}),

donde

V^{\bar{x}} = \mathcal{D}_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{x}}) + \mathcal{D}_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{x} \bar{y}} ) + \mathcal{D}_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{x} \bar{z}} ),

V^{\bar{y}} = \mathcal{D}_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{x}}) + \mathcal{D}_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{y} \bar{y}} ) + \mathcal{D}_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{y} \bar{z}} ),

V^{\bar{z}} = \mathcal{D}_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{x}}) + \mathcal{D}_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{z} \bar{y}} ) + \mathcal{D}_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{z} \bar{z}} ),

que desarrollando las covariantes según lo encontrado al principio del post, quedan:

V^{\bar{x}} = \frac{|\bar{x}^2-1|}{a} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{x}}) + \frac{|\bar{y}^2-1|}{b} \partial_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{x} \bar{y}} ) + \frac{|\bar{z}^2-1|}{c} \partial_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{x} \bar{z}} ),

V^{\bar{y}} = \frac{|\bar{x}^2-1|}{a} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{x}}) + \frac{|\bar{y}^2-1|}{b} \partial_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{y} \bar{y}} ) + \frac{|\bar{z}^2-1|}{c} \partial_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{y} \bar{z}} ),

V^{\bar{z}} = \frac{|\bar{x}^2-1|}{a} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{x}}) + \frac{|\bar{y}^2-1|}{b} \partial_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{z} \bar{y}} ) + \frac{|\bar{z}^2-1|}{c} \partial_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{z} \bar{z}} ).

Por tanto, combiando todo, tenemos:

\Delta \Theta_\beta =

= \frac{3|\bar{x}^2-1|}{2a} \partial_{\bar{x}} \big [ \frac{|\bar{x}^2-1|}{a} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{x}}) +

+ \frac{|\bar{y}^2-1|}{b} \partial_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{x} \bar{y}} ) +

+ \frac{|\bar{z}^2-1|}{c} \partial_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{x} \bar{z}} ) \big ] +

+ \frac{3|\bar{y}^2-1|}{2b} \partial_{\bar{y}} \big [ \frac{|\bar{x}^2-1|}{a} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{x}}) +

+ \frac{|\bar{y}^2-1|}{b} \partial_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{y} \bar{y}} ) +

+ \frac{|\bar{z}^2-1|}{c} \partial_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{y} \bar{z}} ) \big ] +

+ \frac{3|\bar{z}^2-1|}{2c} \partial_{\bar{z}} \big [ \frac{|\bar{x}^2-1|}{a} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{x}}) +

\frac{|\bar{y}^2-1|}{b} \partial_{\bar{y}} ( \alpha \psi^{-6} \hat{A}^{\bar{z} \bar{y}} ) +

+ \frac{|\bar{z}^2-1|}{c} \partial_{\bar{z}} ( \alpha \psi^{-6} \hat{A}^{\bar{z} \bar{z}} ) \big ],

Para terminar, nos quedan la ecuaciónes:

\boxed{\Delta \beta^i = 2\mathcal{D}_j ( \alpha \psi^{-6} \hat{A}^{ij} ) - \frac{1}{3} \mathcal{D}^i \Theta_{\beta} }.

En primer lugar, bajamos el índice de la derivada contravariante:

\mathcal{D}^{\bar{x}} \Theta_\beta = f^{\bar{x} i} \mathcal{D}_i \Theta_\beta = f^{\bar{x} \bar{x}} \mathcal{D}_{\bar{x}} \Theta_\beta + f^{\bar{x} \bar{y}} \mathcal{D}_{\bar{y}} \Theta_\beta + f^{\bar{x} \bar{z}} \mathcal{D}_{\bar{z}} \Theta_\beta = \mathcal{D}_{\bar{x}} \Theta_\beta,

y de la misma manera:

\mathcal{D}^{\bar{y}} = \mathcal{D}_{\bar{y}} y \mathcal{D}^{\bar{z}} = \mathcal{D}_{\bar{z}}.

Así pues, lo que nos queda es:

\Delta \beta^{\bar{x}} = 2 \big [ \mathcal{D}_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{x}}) + \mathcal{D}_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{y}}) + \mathcal{D}_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{z}}) \big ] - \frac{1}{3} \mathcal{D}_{\bar{x}} \Theta_\beta,

\Delta \beta^{\bar{y}} = 2 \big [ \mathcal{D}_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{x}}) + \mathcal{D}_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{y}}) + \mathcal{D}_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{z}}) \big ] - \frac{1}{3} \mathcal{D}_{\bar{y}} \Theta_\beta,

\Delta \beta^{\bar{z}} = 2 \big [ \mathcal{D}_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{x}}) + \mathcal{D}_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{y}}) + \mathcal{D}_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{z}}) \big ] - \frac{1}{3} \mathcal{D}_{\bar{z}} \Theta_\beta,

que al sustituir las derivadas covariantes del tensor dos veces contravariante \hat{A}^{ab} y del escalar \Theta_\beta, por su valor calculado al principio del post, quedan:

\Delta \beta^{\bar{x}} = \frac{2|\bar{x}^2-1|}{a} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{x}}) +

+ \frac{2|\bar{y}^2-1|}{b} \partial_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{y}}) +

+ \frac{2|\bar{z}^2-1|}{c} \partial_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{x} \bar{z}}) - \frac{|\bar{x}^2-1|}{3a} \partial_{\bar{x}} \Theta_\beta,

\Delta \beta^{\bar{y}} = \frac{2|\bar{x}^2-1|}{a} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{x}}) +

+ \frac{2|\bar{y}^2-1|}{b} \partial_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{y}}) +

+ \frac{2|\bar{z}^2-1|}{c} \partial_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{y} \bar{z}}) - \frac{|\bar{y}^2-1|}{3b} \partial_{\bar{y}} \Theta_\beta,

\Delta \beta^{\bar{z}} = \frac{2|\bar{x}^2-1|}{a} \partial_{\bar{x}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{x}}) +

+ \frac{2|\bar{y}^2-1|}{b} \partial_{\bar{y}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{y}}) +

+ \frac{2|\bar{z}^2-1|}{c} \partial_{\bar{z}} (\alpha \psi^{-6} \hat{A}^{\bar{z} \bar{z}}) - \frac{|\bar{z}^2-1|}{3c} \partial_{\bar{z}} \Theta_\beta.

Hasta hace unos dias me creía que los únicos sistemas de coordenadas existentes eran los archiconocidos sistemas cartesiano y polar en dos dimensiones y  cartesiano, cilíndrico y esférico en tres dimensiones. Menuda sorpresa me lleve al descubrir que existen muchos mas. En esta entrada comentaré algunos sistemas de coordenadas alternativos para dos dimensiones y dedicaré otra para los tridimensionales.

Coordenadas parabólicas:

Es un sistema donde las lineas coordenadas son parábolas confocales. Lo que tenemos son dos conjuntos de parábolas:

2y = \frac{x^2}{\sigma^2} - \sigma^2, abiertas hacia el eje +y,

2y = -\frac{x^2}{\tau^2} + \tau^2, abiertas hacia el eje -y,

todas con foco en el origen. De esta manera, las coordenadas parabólicas (\sigma,\tau) vienen determinadas mediante las ecuaciones:

x = \sigma \tau

y = \frac{1}{2}(\tau^2 - \sigma^2)

En la siguiente imagen, vemos la apariencia de este sistema curvilineo:

ParCor2D

Iremos comentando en nuevos posts mas sistemas de coordenadas: coordenadas hiperbólicas, coordenadas elípticas, coordenadas bipolares, etc.

Solo referenciar este post de Francisco R. Villatoro sobre la demostración de Otelbaev del problema del milenio referente a las ecuaciones de Navier-Stokes.

Las ideas principales que comenta son:

  1. Se ha encontrado un contraejemplo (usuario “sup” de dxdy.ru + Stephen Montgomery-Smith + Terry Tao) al principal teorema del artículo, el teorema 6.1, lo que invalida completamente la demostración de Otelbaev.
  2. El matemático Terry Tao, en su artículo “Finite time blowup for an averaged three-dimensional Navier-Stokes equation“, demuestra que la técnica utilizada para atacar el problema no es suficiente, es decir, que no vamos a poder resolver nunca el problema con la técnica que utilizó Otelbaev.
  3. Parece que tenemos problema para rato, segun el propio Tao.

Para mas detalle, al post referido en el inicio de esta entrada 🙂

Compactificamos la primera coordenada mediante \boxed{r = \frac{a \bar{r}}{1 - \bar{r}}}.

El Laplaciano, con esta compactificación, queda:

\Delta = \frac{(a-\bar{r})^4}{a^4} \partial_{\bar{r} \bar{r}} + \frac{2}{\bar{r}} \frac{(a-\bar{r})^4}{a^4} \partial_{\bar{r}} + \frac{(a - \bar{r})^2}{a^2 \bar{r}^2} \partial_{\theta \theta} + \frac{(a-\bar{r})^2}{a^2 \bar{r}^2} \cot \theta \partial_{\theta} + \frac{(a - \bar{r})^2}{a^2 \bar{r}^2 } \csc \theta \partial_{\varphi \varphi}

y las fuentes:

\boxed{\Delta \Theta_{X} = 6 \pi \mathcal{D}^j S^*_j}

\Delta \Theta_X = 6 \pi f^{ji} \mathcal{D}_i S^*_j = 6 \pi ( \mathcal{D}_{\bar{r}} S^*_{\bar{r}} + \mathcal{D}_{\theta} S^*_{\theta} + \mathcal{D}_{\varphi} S^*_{\varphi} ) =

= 6 \pi \frac{1-\bar{r}}{a \bar{r}} ( (\bar{r}-\bar{r}^2) \partial_{\bar{r}} S^*_{\bar{r}} + 2 S^*_{\bar{r}} + \partial_{\theta} S^*_{\theta} + \cot \theta S^*_{\theta} + \csc \theta \partial_{\varphi} S^*_{\varphi} )

\boxed{\Delta X^{i} = 8 \pi f^{ij} S^*_j - \frac{1}{3} \mathcal{D}^i \Theta_X}

Pasando la derivada contravariante a covariante mediante la métrica, queda:

\Delta X^{i} = 8 \pi f^{ij} S^*_j - \frac{1}{3} f^{ik} \mathcal{D}_k \Theta_X.

Definimos ahora

S_X^i := 8 \pi f^{ij} S^*_j - \frac{1}{3} f^{ik} \mathcal{D}_k \Theta_X,

de manera que:

S_X^{\bar{r}} = 8 \pi f^{\bar{r} j} S^*_j - \frac{1}{3} f^{\bar{r} k}\mathcal{D}_{k} \Theta_X = 8 \pi S^*_{\bar{r}} - \frac{1}{3} \mathcal{D}_{\bar{r}} \Theta_X = 8 \pi S^*_{\bar{r}} - \frac{(1-\bar{r})^2}{a^2} \partial_{\bar{r}} \Theta_X,

S_X^{\theta} = 8 \pi f^{\theta j} S^*_j - \frac{1}{3} f^{\theta k} \mathcal{D}_{k} \Theta_X = 8 \pi S^*_{\theta} - \frac{1}{3} \mathcal{D}_{\theta} \Theta_X = 8 \pi S^*_{\theta} - \frac{1-\bar{r}}{a \bar{r}} \partial_{\theta} \Theta_X,

S_X^{\varphi} = 8 \pi f^{\varphi j} S^*_j - \frac{1}{3} f^{\varphi k} \mathcal{D}^{k} \Theta_X = 8 \pi S^*_{\varphi} - \frac{1}{3} \mathcal{D}_{\varphi} \Theta_X = 8 \pi S^*_{\varphi} - \frac{a-\bar{r}}{a \bar{r}} \csc \theta \partial_{\varphi} \Theta_X.

En este punto tenemos que el vector

(S_X^{\bar{r}}(\bar{r},\theta,\varphi),S_X^{\theta}(\bar{r},\theta,\varphi),S_X^{\varphi}(\bar{r},\partial_\theta,\partial_\varphi))

expresado en la base \{ \partial_{\bar{r}}, \theta, \varphi \}. Lo que hacemos ahora es expresar este vector en la nueva base \{ \partial_x, \partial_y, \partial_z \}, de manera que obtenemos

(S_X^{x}(\bar{r},\theta,\varphi),S_X^{y}(\bar{r},\theta,\varphi),S_X^{z}(\bar{r},\theta,\varphi)).

y como es esta base las ecuaciones están desacopladas y \Theta_X es un campo escalar, resolvemos independientemente:

\Delta X^{x} = S_X^{x},

\Delta X^{y} = S_X^{y},

\Delta X^{z} = S_X^{z}.

Finalmente, con el cambio de base inverso, calculamos a partir de (X^{x},X^{y},X^{z}) el vector (X^{\bar{r}},X^\theta,X^\varphi) .

\underline{\hat{A}^{ij} = \mathcal{D}^i X^j + \mathcal{D}^j X^i - \frac{2}{3} \mathcal{D}_k X^k f^{ij}}

volvemos a pasar las derivadas contravariantes a covariantes:

\hat{A}^{ij} = f^{im} \mathcal{D}_m X^j + f^{jn} \mathcal{D}_n X^i - \frac{2}{3} f^{ij} \mathcal{D}_k X^{k}

y obtenemos:

\hat{A}^{\bar{r} \bar{r}} = f^{\bar{r} m} \mathcal{D}_m X^{\bar{r}} + f^{\bar{r} n} \mathcal{D}_n X^{\bar{r}} - \frac{2}{3} \mathcal{D}_k X^{k} = \frac{2}{3}( 2 \mathcal{D}_{\bar{r}} X^{\bar{r}} - \mathcal{D}_{\theta} X^{\theta} - \mathcal{D}_{\varphi} X^{\varphi}) =

= \frac{2(1-\bar{r})}{3a\bar{r}} [ 2(\bar{r}-\bar{r}^2)\partial_{\bar{r}} X^{\bar{r}} - 2 X^{\bar{r}} - \partial_{\theta} X^{\theta} - \cot \theta X^{\theta} - \csc \theta \partial_{\varphi} X^{\varphi} ]

\hat{A}^{\bar{r} \theta} = f^{\bar{r} m} \mathcal{D}_m X^{\theta} + f^{\theta n} \mathcal{D}_n X^{\bar{r}} = \mathcal{D}_{\bar{r}} X^{\theta} + \mathcal{D}_{\theta} X^{\bar{r}} =

= \frac{1-\bar{r}}{a\bar{r}} [ (\bar{r}-\bar{r}^2) \partial_{\bar{r}} X^{\theta} - X^{\theta} + \partial_{\theta} X^{\bar{r}} ],

\hat{A}^{\bar{r} \varphi} = f^{\bar{r} m} \mathcal{D}_m X^{\varphi} + f^{\varphi n} \mathcal{D}_n X^{\bar{r}} = \mathcal{D}_{\bar{r}} X^{\varphi} + \mathcal{D}_{\varphi} X^{\bar{r}} =

= \frac{1-\bar{r}}{a\bar{r}} [ (\bar{r}-\bar{r}^2) \partial_{\bar{r}} X^{\varphi} - X^{\varphi} + \csc \theta \partial_{\varphi} X^{\bar{r}} ],

\hat{A}^{\theta \theta} = f^{\theta m} \mathcal{D}_m X^{\theta} + f^{\theta n} \mathcal{D}_n X^{\theta} - \frac{2}{3} \mathcal{D}_k X^{k} = \frac{2}{3}( - \mathcal{D}_{\bar{r}} X^{\bar{r}} + 2 \mathcal{D}_{\theta} X^{\theta} - \mathcal{D}_{\varphi} X^{\varphi}) =

= \frac{2(1-\bar{r})}{3 a \bar{r}} [ -(\bar{r}-\bar{r}^2) \partial_{\bar{r}} X^{\bar{r}} + X^{\bar{r}} + 2 \partial_{\theta} X^{\theta} - \cot \theta X^{\theta} - \csc \partial_{\varphi} X^{\varphi} ]

\hat{A}^{\theta \varphi} = f^{\theta m} \mathcal{D}_m X^{\varphi} + f^{\varphi n} \mathcal{D}_n X^{\theta} = \mathcal{D}_{\theta} X^{\varphi} + \mathcal{D}_{\varphi} X^{\theta} =

= \frac{1-\bar{r}}{a\bar{r}} [ \partial_{\theta} X^{\varphi} - \cot \theta X^{\varphi} + \csc \theta \partial_{\varphi} X^{\theta} ],

\hat{A}^{\varphi \varphi} = f^{\varphi m} \mathcal{D}_m X^{\varphi} + f^{\varphi n} \mathcal{D}_n X^{\varphi} - \frac{2}{3} \mathcal{D}_k X^{k} = \frac{2}{3}( - \mathcal{D}_{\bar{r}} X^{\bar{r}} - \mathcal{D}_{\theta} X^{\theta} +2 \mathcal{D}_{\varphi} X^{\varphi}) =

= \frac{2(1-\bar{r})}{3 a \bar{r}} [ -(\bar{r} - \bar{r}^2) \partial_{\bar{r}} X^{\bar{r}} + X^{\bar{r}} - \partial_{\theta} X^{\theta} + 2 \cot \theta X^{\theta} + 2 \csc \theta \partial_{\varphi} T^{\varphi} ]

Las dos ecuaciones no lineales correspondientes al factor conforme \psi y al lapse \alpha, como no contienen derivadas covariantes, quedan como las teniamos:

\boxed{\Delta \psi = -2 \pi E^* \psi^{-1} - \frac{1}{8}(f_{il} f_{jm} \hat{A}^{lm} \hat{A}^{ij}) \psi^{-7} }

\boxed{\Delta (\alpha \psi) = [ 2 \pi (E^* + 2 S^*) \psi^{-7} + \frac{1}{8}(f_{il} f_{jm} \hat{A}^{lm} \hat{A}^{ij}) \psi^{-8} ] (\alpha \psi) }

Finalmente, para el shift \beta y su ecuación auxiliar tenemos:

\boxed{\Delta \Theta_{\beta} = \frac{3}{4} \mathcal{D}_i \mathcal{D}_j (2 \alpha \psi^{-6} \hat{A}^{ij} )}

\boxed{\Delta \beta^i = \mathcal{D}_j ( 2 \alpha \psi^{-6} \hat{A}^{ij} ) - \frac{1}{3} \mathcal{D}^i \Theta_{\beta} }

que lo tratamos en este post.

Compactificamos ahora mediante \boxed{x_i =a_i \tan \frac{\pi \bar{x}_i}{2}}.

El Laplaciano, con esta nueva compactificación, queda:

\Delta = \frac{4}{a^2 \pi^2} \cos^4 \frac{\pi \bar{x}}{2} \partial_{\bar{x} \bar{x}} + \frac{\sin (\pi \bar{x})(1+\cos (\pi \bar{x}))}{a^2 \pi^2} \partial_{\bar{x}} +

+ \frac{4}{b^2 \pi^2} \cos^4 \frac{\pi \bar{y}}{2} \partial_{\bar{y} \bar{y}} + \frac{\sin (\pi \bar{y})(1+\cos (\pi \bar{y}))}{b^2 \pi^2} \partial_{\bar{y}} +

+ \frac{4}{c^2 \pi^2} \cos^4 \frac{\pi \bar{z}}{2} \partial_{\bar{z} \bar{z}} + \frac{\sin (\pi \bar{z})(1+\cos (\pi \bar{z}))}{c^2 \pi^2} \partial_{\bar{z}}

y las fuentes:

\boxed{\Delta \Theta_{X} = 6 \pi \mathcal{D}^j S^*_j}

\Delta \Theta_X = 6 \pi f^{ji} \mathcal{D}_i S^*_j = 6 \pi ( \mathcal{D}_{\bar{x}} S^*_{\bar{x}} + \mathcal{D}_{\bar{y}} S^*_{\bar{y}} + \mathcal{D}_{\bar{z}} S^*_{\bar{z}} ) =

= 6 \pi ( \frac{1 + \cos ( \pi \bar{x}) }{a \pi} \partial_{\bar{x}} S^*_{\bar{x}} + \frac{1 + \cos (\pi \bar{y})}{b \pi} \partial_{\bar{y}} S^*_{\bar{y}}+ \frac{1 + \cos (\pi \bar{z})}{c \pi} \partial_{\bar{z}} S^*_{\bar{z}})

\boxed{\Delta X^{i} = 8 \pi f^{ij} S^*_j - \frac{1}{3} \mathcal{D}^i \Theta_X}

Pasamos la derivada contravariante a covariante mediante la métrica:

\Delta X^{i} = 8 \pi f^{ij} S^*_j - \frac{1}{3} f^{ik} \mathcal{D}_k \Theta_X,

y como \Theta_X es un campo escalar, nos queda:

\Delta X^{\bar{x}} = 8 \pi f^{\bar{x} j} S^*_j - \frac{1}{3} f^{\bar{x} k}\mathcal{D}_{k} \Theta_X = 8 \pi S^*_{\bar{x}} - \frac{1}{3} \mathcal{D}_{\bar{x}} \Theta_X =

= 8 \pi S^*_{\bar{x}} - \frac{1+\cos(\pi \bar{x})}{3a \pi} \partial_{\bar{x}} \Theta_X,

\Delta X^{\bar{y}} = 8 \pi f^{\bar{y} j} S^*_j - \frac{1}{3} f^{\bar{y} k} \mathcal{D}_{k} \Theta_X = 8 \pi S^*_{\bar{y}} - \frac{1}{3} \mathcal{D}_{\bar{y}} \Theta_X =

= 8 \pi S^*_{\bar{y}} - \frac{1 +\cos(\pi \bar{y})}{3b \pi} \partial_{\bar{y}} \Theta_X,

\Delta X^{\bar{z}} = 8 \pi f^{\bar{z} j} S^*_j - \frac{1}{3} f^{\bar{z} k} \mathcal{D}^{k} \Theta_X = 8 \pi S^*_{\bar{z}} - \frac{1}{3} \mathcal{D}_{\bar{z}} \Theta_X =

= 8 \pi S^*_{\bar{z}} - \frac{1 + \cos(\pi \bar{z})}{3c \pi} \partial_{\bar{z}} \Theta_X.

\underline{\hat{A}^{ij} = \mathcal{D}^i X^j + \mathcal{D}^j X^i - \frac{2}{3} \mathcal{D}_k X^k f^{ij}}

volvemos a pasar las derivadas contravariantes a covariantes:

\hat{A}^{ij} = f^{im} \mathcal{D}_m X^j + f^{jn} \mathcal{D}_n X^i - \frac{2}{3} f^{ij} \mathcal{D}_k X^{k}

y obtenemos:

\hat{A}^{\bar{x} \bar{x}} = f^{\bar{x} m} \mathcal{D}_m X^{\bar{x}} + f^{\bar{x} n} \mathcal{D}_n X^{\bar{x}} - \frac{2}{3} \mathcal{D}_k X^{k} = \frac{2}{3}( 2 \mathcal{D}_{\bar{x}} X^{\bar{x}} - \mathcal{D}_{\bar{y}} X^{\bar{y}} - \mathcal{D}_{\bar{z}} X^{\bar{z}}) =

= \frac{4}{3 a \pi}(1+\cos(\pi \bar{x})) \partial_{\bar{x}} X^{\bar{x}} - \frac{2}{3 b \pi}(1 + \cos(\pi \bar{y})) \partial_{\bar{y}} X^{\bar{y}} - \frac{2}{3 c \pi}(1+\cos(\pi \bar{z})) \partial_{\bar{z}} X^{\bar{z}}

\hat{A}^{\bar{x} \bar{y}} = f^{\bar{x} m} \mathcal{D}_m X^{\bar{y}} + f^{\bar{y} n} \mathcal{D}_n X^{\bar{x}} = \mathcal{D}_{\bar{x}} X^{\bar{y}} + \mathcal{D}_{\bar{y}} X^{\bar{x}} =

= \frac{1+\cos (\pi \bar{x})}{a \pi} \partial_{\bar{x}} X^{\bar{y}} + \frac{1+ \cos(\pi \bar{y})}{b \pi} \partial_{\bar{y}} X^{\bar{x}},

\hat{A}^{\bar{x} \bar{z}} = f^{\bar{x} m} \mathcal{D}_m X^{\bar{z}} + f^{\bar{z} n} \mathcal{D}_n X^{\bar{x}} = \mathcal{D}_{\bar{x}} X^{\bar{z}} + \mathcal{D}_{\bar{z}} X^{\bar{x}} =

= \frac{1+\cos (\pi \bar{x})}{a \pi} \partial_{\bar{x}} X^{\bar{z}} + \frac{1+ \cos(\pi \bar{z})}{c \pi} \partial_{\bar{z}} X^{\bar{x}},

\hat{A}^{\bar{y} \bar{y}} = f^{\bar{y} m} \mathcal{D}_m X^{\bar{y}} + f^{\bar{y} n} \mathcal{D}_n X^{\bar{y}} - \frac{2}{3} \mathcal{D}_k X^{k} = \frac{2}{3}( - \mathcal{D}_{\bar{x}} X^{\bar{x}} + 2 \mathcal{D}_{\bar{y}} X^{\bar{y}} - \mathcal{D}_{\bar{z}} X^{\bar{z}}) =

= -\frac{2}{3a \pi}(1+\cos(\pi \bar{x})) \partial_{\bar{x}} X^{\bar{x}} + \frac{4}{3b \pi}(1 + \cos(\pi \bar{y})) \partial_{\bar{y}} X^{\bar{y}} - \frac{2}{3c \pi}(1+\cos(\pi \bar{z})) \partial_{\bar{z}} X^{\bar{z}}

\hat{A}^{\bar{y} \bar{z}} = f^{\bar{y} m} \mathcal{D}_m X^{\bar{z}} + f^{\bar{z} n} \mathcal{D}_n X^{\bar{y}} = \mathcal{D}_{\bar{y}} X^{\bar{z}} + \mathcal{D}_{\bar{z}} X^{\bar{y}} =

= \frac{1+\cos (\pi \bar{y})}{b \pi} \partial_{\bar{y}} X^{\bar{z}} + \frac{1+ \cos(\pi \bar{z})}{c \pi} \partial_{\bar{z}} X^{\bar{y}},

\hat{A}^{\bar{z} \bar{z}} = f^{\bar{z} m} \mathcal{D}_m X^{\bar{z}} + f^{\bar{z} n} \mathcal{D}_n X^{\bar{z}} - \frac{2}{3} \mathcal{D}_k X^{k} = \frac{2}{3}( - \mathcal{D}_{\bar{x}} X^{\bar{x}} - \mathcal{D}_{\bar{y}} X^{\bar{y}} +2 \mathcal{D}_{\bar{z}} X^{\bar{z}}) =

= -\frac{2}{3a \pi}(1+\cos(\pi \bar{x})) \partial_{\bar{x}} X^{\bar{x}} - \frac{2}{3b \pi}(1 + \cos(\pi \bar{y})) \partial_{\bar{y}} X^{\bar{y}} + \frac{4}{3c \pi}(1+\cos(\pi \bar{z})) \partial_{\bar{z}} X^{\bar{z}}

Las dos ecuaciones no lineales correspondientes al factor conforme \psi y al lapse \alpha, como no contienen derivadas covariantes, quedan como las teniamos:

\boxed{\Delta \psi = -2 \pi E^* \psi^{-1} - \frac{1}{8}(f_{il} f_{jm} \hat{A}^{lm} \hat{A}^{ij}) \psi^{-7} }

\boxed{\Delta (\alpha \psi) = [ 2 \pi (E^* + 2 S^*) \psi^{-7} + \frac{1}{8}(f_{il} f_{jm} \hat{A}^{lm} \hat{A}^{ij}) \psi^{-8} ] (\alpha \psi) }

Finalmente, para el shift \beta y su ecuación auxiliar tenemos:

\boxed{\Delta \Theta_{\beta} = \frac{3}{4} \mathcal{D}_i \mathcal{D}_j (2 \alpha \psi^{-6} \hat{A}^{ij} )}

\boxed{\Delta \beta^i = \mathcal{D}_j ( 2 \alpha \psi^{-6} \hat{A}^{ij} ) - \frac{1}{3} \mathcal{D}^i \Theta_{\beta} }

que comentamos en este post.

Compactificamos mediante \boxed{x_i = a_i \, \mbox{arctanh} \, \bar{x}_i }.

En este caso, el operador Laplaciano de todas las ecuaciones queda:

\Delta = \frac{(\bar{x}^2 - 1)^2}{a^2} \partial_{\bar{x} \bar{x}} + \frac{2 \bar{x}(\bar{x}^2 - 1)}{a^2} \partial_{\bar{x}} +

+ \frac{(\bar{y}^2 - 1)^2}{b^2} \partial_{\bar{y} \bar{y}} + \frac{2 \bar{y}(\bar{y}^2 - 1)}{b^2} \partial_{\bar{y}} +

+ \frac{(\bar{z}^2 - 1)^2}{c^2} \partial_{\bar{z} \bar{z}} + \frac{2 \bar{z}(\bar{z}^2 - 1)}{c^2} \partial_{\bar{z}}

y las fuentes:

\boxed{\Delta \Theta_{X} = 6 \pi \mathcal{D}^j S^*_j}

\Delta \Theta_X = 6 \pi f^{ji} \mathcal{D}_i S^*_j = 6 \pi ( \mathcal{D}_{\bar{x}} S^*_{\bar{x}} + \mathcal{D}_{\bar{y}} S^*_{\bar{y}} + \mathcal{D}_{\bar{z}} S^*_{\bar{z}} ) =

= 6 \pi ( \frac{|\bar{x}^2-1|}{a} \partial_{\bar{x}} S^*_{\bar{x}} + \frac{|\bar{y}^2-1|}{b} \partial_{\bar{y}} S^*_{\bar{y}}+ \frac{|\bar{z}^2-1|}{c} \partial_{\bar{z}} S^*_{\bar{z}})

\boxed{\Delta X^{i} = 8 \pi f^{ij} S^*_j - \frac{1}{3} \mathcal{D}^i \Theta_X}

Pasamos la derivada contravariante a covariante mediante la métrica:

\Delta X^{i} = 8 \pi f^{ij} S^*_j - \frac{1}{3} f^{ik} \mathcal{D}_k \Theta_X,

y como \Theta_X es un campo escalar, nos queda:

\Delta X^{\bar{x}} = 8 \pi f^{\bar{x} j} S^*_j - \frac{1}{3} f^{\bar{x} k}\mathcal{D}_{k} \Theta_X = 8 \pi S^*_{\bar{x}} - \frac{1}{3} \mathcal{D}_{\bar{x}} \Theta_X =

= 8 \pi S^*_{\bar{x}} - \frac{1}{3} \frac{|\bar{x}^2 - 1|}{a} \partial_{\bar{x}} \Theta_X,

\Delta X^{\bar{y}} = 8 \pi f^{\bar{y} j} S^*_j - \frac{1}{3} f^{\bar{y} k} \mathcal{D}_{k} \Theta_X = 8 \pi S^*_{\bar{y}} - \frac{1}{3} \mathcal{D}_{\bar{y}} \Theta_X =

= 8 \pi S^*_{\bar{y}} - \frac{1}{3} \frac{|\bar{y}^2 - 1|}{b} \partial_{\bar{y}} \Theta_X,

\Delta X^{\bar{z}} = 8 \pi f^{\bar{z} j} S^*_j - \frac{1}{3} f^{\bar{z} k} \mathcal{D}^{k} \Theta_X = 8 \pi S^*_{\bar{z}} - \frac{1}{3} \mathcal{D}_{\bar{z}} \Theta_X =

= 8 \pi S^*_{\bar{z}} - \frac{1}{3} \frac{|\bar{z}^2 - 1|}{c} \partial_{\bar{z}} \Theta_X.

\underline{\hat{A}^{ij} = \mathcal{D}^i X^j + \mathcal{D}^j X^i - \frac{2}{3} \mathcal{D}_k X^k f^{ij}}

volvemos a pasar las derivadas contravariantes a covariantes:

\hat{A}^{ij} = f^{im} \mathcal{D}_m X^j + f^{jn} \mathcal{D}_n X^i - \frac{2}{3} f^{ij} \mathcal{D}_k X^{k}

y obtenemos:

\hat{A}^{\bar{x} \bar{x}} = f^{\bar{x} m} \mathcal{D}_m X^{\bar{x}} + f^{\bar{x} n} \mathcal{D}_n X^{\bar{x}} - \frac{2}{3} \mathcal{D}_k X^{k} = \frac{2}{3}( 2 \mathcal{D}_{\bar{x}} X^{\bar{x}} - \mathcal{D}_{\bar{y}} X^{\bar{y}} - \mathcal{D}_{\bar{z}} X^{\bar{z}}) =

= \frac{2}{3}(\frac{2|\bar{x}^2-1|}{a} \partial_{\bar{x}} X^{\bar{x}} - \frac{|\bar{y}^2-1|}{b} \partial_{\bar{y}} X^{\bar{y}} - \frac{|\bar{z}^2-1|}{c} \partial_{\bar{z}} X^{\bar{z}}),

\hat{A}^{\bar{x} \bar{y}} = f^{\bar{x} m} \mathcal{D}_m X^{\bar{y}} + f^{\bar{y} n} \mathcal{D}_n X^{\bar{x}} = \mathcal{D}_{\bar{x}} X^{\bar{y}} + \mathcal{D}_{\bar{y}} X^{\bar{x}} =

= \frac{|\bar{x}^2-1|}{a} \partial_{\bar{x}} X^{\bar{y}} + \frac{|\bar{y}^2 - 1|}{b} \partial_{\bar{y}} X^{\bar{x}},

\hat{A}^{\bar{x} \bar{z}} = f^{\bar{x} m} \mathcal{D}_m X^{\bar{z}} + f^{\bar{z} n} \mathcal{D}_n X^{\bar{x}} = \mathcal{D}_{\bar{x}} X^{\bar{z}} + \mathcal{D}_{\bar{z}} X^{\bar{x}} =

= \frac{|\bar{x}^2-1|}{a} \partial_{\bar{x}} X^{\bar{z}} + \frac{|\bar{z}^2 - 1|}{c} \partial_{\bar{z}} X^{\bar{x}},

\hat{A}^{\bar{y} \bar{y}} = f^{\bar{y} m} \mathcal{D}_m X^{\bar{y}} + f^{\bar{y} n} \mathcal{D}_n X^{\bar{y}} - \frac{2}{3} \mathcal{D}_k X^{k} = \frac{2}{3}( - \mathcal{D}_{\bar{x}} X^{\bar{x}} + 2 \mathcal{D}_{\bar{y}} X^{\bar{y}} - \mathcal{D}_{\bar{z}} X^{\bar{z}}) =

= \frac{2}{3}(-\frac{|\bar{x}^2-1|}{a} \partial_{\bar{x}} X^{\bar{x}} + \frac{2|\bar{y}^2-1|}{b} \partial_{\bar{y}} X^{\bar{y}} - \frac{|\bar{z}^2-1|}{c} \partial_{\bar{z}} X^{\bar{z}}),

\hat{A}^{\bar{y} \bar{z}} = f^{\bar{y} m} \mathcal{D}_m X^{\bar{z}} + f^{\bar{z} n} \mathcal{D}_n X^{\bar{y}} = \mathcal{D}_{\bar{y}} X^{\bar{z}} + \mathcal{D}_{\bar{z}} X^{\bar{y}} =

= \frac{|\bar{y}^2-1|}{b} \partial_{\bar{y}} X^{\bar{z}} + \frac{|\bar{z}^2 - 1|}{c} \partial_{\bar{z}} X^{\bar{y}},

\hat{A}^{\bar{z} \bar{z}} = f^{\bar{z} m} \mathcal{D}_m X^{\bar{z}} + f^{\bar{z} n} \mathcal{D}_n X^{\bar{z}} - \frac{2}{3} \mathcal{D}_k X^{k} = \frac{2}{3}( - \mathcal{D}_{\bar{x}} X^{\bar{x}} - \mathcal{D}_{\bar{y}} X^{\bar{y}} +2 \mathcal{D}_{\bar{z}} X^{\bar{z}}) =

= \frac{2}{3}(-\frac{|\bar{x}^2-1|}{a} \partial_{\bar{x}} X^{\bar{x}} - \frac{|\bar{y}^2-1|}{b} \partial_{\bar{y}} X^{\bar{y}} + \frac{2|\bar{z}^2-1|}{c} \partial_{\bar{z}} X^{\bar{z}}).

Las dos ecuaciones no lineales correspondientes al factor conforme \psi y al lapse \alpha, como no contienen derivadas covariantes, quedan como las teniamos:

\boxed{\Delta \psi = -2 \pi E^* \psi^{-1} - \frac{1}{8}(f_{il} f_{jm} \hat{A}^{lm} \hat{A}^{ij}) \psi^{-7} }

\boxed{\Delta (\alpha \psi) = [ 2 \pi (E^* + 2 S^*) \psi^{-7} + \frac{1}{8}(f_{il} f_{jm} \hat{A}^{lm} \hat{A}^{ij}) \psi^{-8} ] (\alpha \psi) }

Finalmente, para el shift \beta y su ecuación auxiliar tenemos:

\boxed{\Delta \Theta_{\beta} = \frac{3}{4} \mathcal{D}_i \mathcal{D}_j (2 \alpha \psi^{-6} \hat{A}^{ij} )}

\boxed{\Delta \beta^i = \mathcal{D}_j ( 2 \alpha \psi^{-6} \hat{A}^{ij} ) - \frac{1}{3} \mathcal{D}^i \Theta_{\beta} }

que las tratamos es este nuevo post.

Los cálculos para cilíndricas, esféricas, esféricas compactificadas y cartesianas compactificadas ya están términados. Para hacer una pequeña comprobación de que son correctos, vamos a calcular, en cada caso, como quedaría la divergencia de un campo vectorial \mathcal{D}_i X^i utilizando las derivadas covariantes encontrados en los enlaces anteriores y compararla con el resultado que obtendriamos utilizando la fórmula para la divergencia en coordenadas curvilineas q^i:

\nabla \cdot \mathbf{X} = \frac{1}{\Pi_j h_j} \frac{\partial}{\partial q^i} (X^i \Pi_{j \neq i} h_j).

\{ \partial_r, \partial_{\theta}, \partial_z \}

\mathcal{D}_i X^i = \mathcal{D}_r X^r + \mathcal{D}_{\theta} X^{\theta} + \mathcal{D}_{z} X^{z} = \partial_r X^r + \partial_{\theta} X^{\theta} + \frac{1}{r} X^{\theta} + \partial_z X^z,

\{ \partial_r, \frac{1}{r} \partial_{\theta}, \partial_z \}

\mathcal{D}_i X^i = \partial_r X^r + \frac{1}{r} \partial_{\theta} X^{\theta} + \frac{1}{r} X^{\theta} + \partial_z X^z,

y con la fórmula para curvilineas obtenemos:

divCyl

\{ \partial_r, \partial_{\theta}, \partial_{\varphi} \}

\mathcal{D}_i X^i = \partial_r X^r + \frac{2}{r} X^r + \partial_{\theta} X^{\theta} + \cot \theta X^{\theta} + \partial_{\varphi} X^{\varphi},

\{ \partial_r, \frac{1}{r} \partial_{\theta}, \frac{\csc \theta}{r} \partial_{\varphi} \}

\mathcal{D}_i X^i = \frac{1}{r} \big [ r \partial_r X^r + 2 X^r + \partial_{\theta} X^{\theta} + \cot \theta X^{\theta} + \csc \theta \partial_{\varphi} X^{\varphi} \big ],

y con la fórmula para curvilineas obtenemos:

divSph

\{ \partial_{\bar{r}}, \partial_{\theta}, \partial_{\varphi} \}

\mathcal{D}_i X^i = \partial_{\bar{r}} X^r + \frac{1+\bar{r}}{1-\bar{r}} \frac{2}{\bar{r}} X^r + \partial_{\theta} X^{\theta} + \cot \theta X^{\theta} + \partial_{\varphi} X^{\varphi},

\{ \frac{(1-\bar{r})^2}{a} \partial_{\bar{r}}, \frac{1-\bar{r}}{a \bar{r}} \partial_{\theta}, \frac{1-\bar{r}}{a \bar{r}} \csc \theta \partial_{\varphi} \}

\mathcal{D}_i X^i = \frac{1 - \bar{r}}{a \bar{r}} \big [ (\bar{r} - \bar{r}^2) \partial_{\bar{r}} X^r + 2 X^r + \partial_{\theta} X^{\theta} + \cot \theta X^{\theta} + \csc \theta \partial_{\varphi} X^{\varphi} \big ],

y con la fórmula para curvilineas obtenemos:

divSphCom1

\{ \partial_{\bar{r}}, \partial_{\theta}, \partial_{\varphi} \}

\mathcal{D}_i X^i = \partial_{\bar{r}} X^r + \frac{4 \bar{r}}{(1-\bar{r})^2 \, \mbox{\scriptsize arctanh} \, \bar{r}} X^r + \partial_{\theta} X^{\theta} + \cot \theta X^{\theta} + \partial_{\varphi} X^{\varphi},

\{ \frac{1-\bar{r}^2}{a} \partial_{\bar{r}}, \frac{1}{a \, \mbox{\scriptsize arctanh} \, \bar{r} } \partial_{\theta}, \frac{\csc \theta}{a \, \mbox{\scriptsize arctanh} \, \bar{r}} \partial_{\varphi} \}

\mathcal{D}_i X^i = \frac{1}{a \, \mbox{\scriptsize arctanh} \, \bar{r}} \big [ (1-\bar{r}^2) \, \mbox{arctanh} \, \bar{r} \partial_{\bar{r}} X^r + 2 X^r +

+ \partial_{\theta} X^{\theta} + \cot \theta X^{\theta} + \csc \theta \partial_{\varphi} X^{\varphi} \big ],

y con la fórmula para curvilineas obtenemos:

divSphCom2

\{ \partial_{\bar{r}}, \partial_{\theta}, \partial_{\varphi} \}

\partial_{\bar{r}} X^r + ( \pi \tan \frac{\pi \bar{r}}{2} + 2 \pi \csc (\pi \bar{r} ) ) X^r + \partial_{\theta} X^{\theta} + \cot \theta X^{\theta} + \partial_{\varphi} X^{\varphi},

\{ \frac{1 + \cos ( \pi \bar{r})}{a \pi} \partial_{\bar{r}} , \frac{ \cot \frac{ \pi \bar{r} }{2} }{a} \partial_{\theta} , \frac{ \cot \frac{ \pi \bar{r} }{2} }{a} \csc \theta \partial_{\varphi} \}

\mathcal{D}_i X^i = \frac{ \cot \frac{ \pi \bar{r} }{2} }{a} \big [ \frac{1 + \cos ( \pi \bar{r} ) }{\bar{r}} \tan \frac{\pi \bar{r}}{2} \partial_{ \bar{r} } X^r + 2 X^r +

+ \partial_{\theta} X^{\theta} + \cot \theta X^{\theta} + \csc \theta \partial_{\varphi} X^{\varphi} \big ],

y con la fórmula para curvilineas obtenemos:

divSphCom3

\{ \partial_{\bar{x}} , \partial_{\bar{y}}, \partial_{\bar{z}} \}

\mathcal{D}_i X^i = \partial_{\bar{x}} X^{\bar{x}} + \frac{2 \bar{x}}{1-\bar{x}^2} X^{\bar{x}} + \partial_{\bar{y}} X^{\bar{y}} + \frac{2 \bar{y}}{1-\bar{y}^2} X^{\bar{y}} + \partial_{\bar{z}} X^{\bar{z}} + \frac{2 \bar{z}}{1-\bar{z}^2} X^{\bar{z}},

\{ \frac{|\bar{x}^2 -1|}{a} \partial_{\bar{x}} , \frac{|\bar{y}^2 -1|}{b} \partial_{\bar{y}}, \frac{|\bar{z}^2 -1|}{c} \partial_{\bar{z}} \}

\mathcal{D}_i X^i = \frac{|\bar{x}^2 -1|}{a} \partial_{\bar{x}} X^{\bar{x}} + \frac{|\bar{y}^2 - 1|}{b} \partial_{\bar{y}} X^{\bar{y}} + \frac{|\bar{z}^2 - 1|}{c} \partial_{\bar{z}} X^{\bar{z}},

y con la fórmula para curvilineas obtenemos:

divCarCom1

\{ \partial_{\bar{x}} , \partial_{\bar{y}}, \partial_{\bar{z}} \}

\mathcal{D}_i X^i = \partial_{\bar{x}} X^{\bar{x}} + \pi \tan \frac{\pi \bar{x}}{2} X^{\bar{x}} + \partial_{\bar{y}} X^{\bar{y}} + \pi \tan \frac{\pi \bar{y}}{2} X^{\bar{y}} + \partial_{\bar{z}} X^{\bar{z}} + \pi \tan \frac{\pi \bar{z}}{2} X^{\bar{z}}

\{ \frac{1 + \cos ( \pi \bar{x} ) }{a \pi} \partial_{\bar{x}} , \frac{1 + \cos ( \pi \bar{y} ) }{b \pi} \partial_{\bar{y}}, \frac{1 + \cos ( \pi \bar{z} ) }{c \pi} \partial_{\bar{z}} \}

\mathcal{D}_i X^i = \frac{1 + \cos ( \pi \bar{x} ) }{a \pi} \partial_{\bar{x}} X^{\bar{x}} + \frac{1 + \cos ( \pi \bar{y} ) }{b \pi} \partial_{\bar{y}} X^{\bar{y}} + \frac{1 + \cos ( \pi \bar{z} ) }{c \pi} \partial_{\bar{z}} X^{\bar{z}},

y con la fórmula para curvilineas obtenemos:

divCarCom2

Seguimos utilizando la misma función mencionada aquí.

Compactificaremos de dos manera diferentes:

\boxed{\boxed{x = a \, \mbox{arctanh} \, \bar{x}, y = b \, \mbox{arctanh} \, \bar{y}, z = c \, \mbox{arctanh} \, \bar{z} }}

Para la base \{ \partial_{\bar{x}}, \partial_{\bar{y}}, \partial_{\bar{z}}\}, los símbolos de Christoffel y las derivadas covariantes quedan (utilizamos X,Y,Z en Mathematica para representar \bar{x},\bar{y},\bar{z}):

ChrSym_CarCom1

CovDer_CarCom1

Para la base \{ \frac{|-1+\bar{x}^2|}{a} \partial_{\bar{x}}, \frac{|-1+\bar{y}^2|}{b} \partial_{\bar{y}}, \frac{|-1+\bar{z}^2|}{c} \partial_ {\bar{z}}\}, los coeficientes de rotación de Ricci y las derivadas covariantes quedan:

RotRic_CarComNor1

CovDer_CarComNor1

\boxed{\boxed{x = a \tan \frac{\pi \bar{x}}{2}, b \tan \frac{\pi \bar{y}}{2}, c \tan \frac{\pi \bar{z}}{2} }}

Para la base \{ \partial_{\bar{x}}, \partial_{\bar{y}}, \partial_{\bar{z}}\}, los símbolos de Christoffel y las derivadas covariantes quedan:

ChrSym_CarCom2

CovDer_CarCom2

Para la base \{ \frac{1+\cos (\pi \bar{x})}{a \pi} \partial_{\bar{x}}, \frac{1 + \cos(\pi \bar{y})}{b \pi} \partial_{\bar{y}}, \frac{1 + \cos(\pi \bar{z})}{c \pi} \partial_ {\bar{z}}\}, los coeficientes de rotación de Ricci y las derivadas covariantes quedan:

RotRic_CarComNor2

CovDer_CarComNor2

No hay que confundir los dos conceptos del título del post. Por un lado tenemos claramente definido el concepto de serie convergente: una serie es convergente si la sucesión de las sumas parciales tiene límite. Por otro, existen métodos de sumación, que pretenden extender, de manera consistente, la asignación de un valor de suma a una serie divergente en el sentido anterior.

Sin entrar en criterios de convergencia de series, existe una manera muy sencilla de saber si una serie diverge: si la sucesión de términos no tiende a cero, entonces la serie diverge:

\lim_{n \rightarrow \infty} a_n \neq 0 \Rightarrow \sum_{n=1}^{\infty} a_n = \infty.

Por ejemplo, las series:

\sum_{n=1}^{\infty} (-1)^{n-1}n = 1 - 2 + 3 - \ldots

\sum_{n=1}^{\infty} n = 1 + 2 + 3 + \ldots

son divergentes, ya que la sucesión de sus correspondientes sumas parciales (1,-1,2,-2\ldots; 1,3,6,10,\ldots ) no tiene límite o, por el criterio anterior, \lim_{n \rightarrow \infty} (-1)^{n-1} n \neq 0 y \lim_{n \rightarrow \infty} n \neq 0.

Sin embargo, podemos asignarle un valor a esta suma (jugando apropiadamente con los términos :-)). En el primer caso, ya Euler encontró que:

\sum_{n=1}^{\infty} n = 1 - 2 +3 -4 + \ldots = \frac{1}{4},

y posteriormente se establecieron métodos bien definidos para encontrar sumas generalizadas de series divergentes.

En el segundo caso, podemos utilizar la extensión analítica (de funciones analíticas y sus extensiones comenté algo en este post) de la función zeta de Riemann (una nota curiosa aquí):

\zeta(s):=\sum_{n=1}^{\infty} \frac{1}{n^s} = 1 + \frac{1}{2^s} + \frac{1}{3^s} + \ldots,

ya que nuestra serie no es mas que \zeta(-1) y, como se puede demostrar que:

\zeta(-s) = \frac{B_{s+1}}{s+1},

donde B_n son los números de Bernoulli, entonces tenemos que:

\sum_{n=1}^{\infty} n = 1 + 2 + 3 + 4 + \ldots = -\frac{1}{12},

En este post, Tao da una interpretación matemáticamente consistente de lo que son estos valores.

Seguimos utilizando la misma función mencionada aquí.

Compactificaremos de tres manera diferentes:

\boxed{\boxed{r = \frac{a \bar{r}}{1 - \bar{r}}}} (y no \frac{a \bar{r}}{a - \bar{r}} como escribimos en este post)

Para la base \{ \partial_{\bar{r}}, \partial_{\theta}, \partial_{\varphi}\}, los símbolos de Christoffel y las derivadas covariantes quedan (\bar{r} lo representamos mediante R en Mathematica):

ChrSym_SphCom

CovDer_SphCom

Para la base \{ \frac{(1-\bar{r})^2}{a} \partial_{\bar{r}}, \frac{1-\bar{r}}{a \bar{r}} \partial_{\theta}, \frac{1-\bar{r}}{a \bar{r}}\csc \theta \partial_ {\varphi}\}, los coeficientes de rotación de Ricci y las derivadas covariantes quedan:

RotRic_SphComNor

CovDer_SphComNor

\boxed{\boxed{r=a\, \mbox{arctanh} \bar{r}}}

Para la base \{ \partial_{\bar{r}}, \partial_{\theta}, \partial_{\varphi}\}, los símbolos de Christoffel y las derivadas covariantes quedan:

ChrSym_SphCom2

CovDer_SphCom2

Para la base \{ \frac{1-\bar{r}^2}{a}\partial_{\bar{r}}, \frac{1}{a\, \mbox{\scriptsize arctanh}\, \bar{r}} \partial_{\theta}, \frac{\csc \theta}{a\, \mbox{\scriptsize arctanh} \,\bar{r}} \partial_ {\varphi}\}, los coeficientes de rotación de Ricci y las derivadas covariantes quedan:

RotRic_SphComNor2

CovDer_SphComNor2

\boxed{\boxed{r = a \tan \frac{\pi \bar{r}}{2}}}

Para la base \{ \partial_{\bar{r}}, \partial_{\theta}, \partial_{\varphi}\}, los símbolos de Christoffel y las derivadas covariantes quedan:

ChrSym_SphCom3

CovDer_SphCom3

Para la base \{ \frac{1+\cos \pi \bar{r}}{a \pi} \partial_{\bar{r}}, \frac{\cot \frac{\pi \bar{r}}{2}}{a} \partial_{\theta}, \frac{\cot \frac{\pi \bar{r}}{2}}{a} \csc \theta \partial_ {\varphi}\}, los coeficientes de rotación de Ricci y las derivadas covariantes quedan:

RotRic_SphComNor3

CovDer_SphComNor3

Seguimos utilizando la misma función mencionada aquí.

Para la base \{ \partial_r, \partial_{\theta}, \partial_{\varphi}\}, los símbolos de Christoffel y las derivadas covariantes quedan:

   ChrSym_Sph

CovDer_Sph

Para la base \{ \partial_r, \frac{1}{r} \partial_{\theta}, \frac{\csc \theta}{r} \partial_ {\varphi}\}, los coeficientes de rotación de Ricci y las derivadas covariantes quedan:

RotRic_SphNor

CovDer_SphNor

Ya tenemos nuestra función lista para realizar todos estos cálculos de manera automática.

Para la base \{ \partial_r, \partial_{\theta}, \partial_z\}, los símbolos de Christoffel y las derivadas covariantes quedan:

ChrSym_Cyl

CovDer_Cyl

Para la base \{ \partial_r, \frac{1}{r} \partial_{\theta}, \partial_ z\}, los coeficientes de rotación de Ricci y las derivadas covariantes quedan:

RotRic_CylNor

CovDer_CylNor

febrero 2014
L M X J V S D
« Ene   Mar »
 12
3456789
10111213141516
17181920212223
2425262728