Coordenadas ortogonales

De testwiki
Ir a la navegación Ir a la búsqueda
Una transformación conforme que actúa sobre una cuadrícula rectangular. Téngase en cuenta que se conserva la ortogonalidad de la cuadrícula curva

Un sistema de coordenadas ortogonales es un sistema de coordenadas tal que en cada punto los vectores tangentes a las curvas coordenadas son ortogonales entre sí. Este tipo de coordenadas pueden definirse sobre un espacio euclídeo o más generalmente sobre una variedad riemanniana o pseudoriemanniana.

Motivación

Si bien las operaciones vectoriales y las leyes físicas son normalmente más fáciles de interpretar en coordenadas cartesianas, las coordenadas ortogonales no cartesianas se utilizan a menudo en su lugar para la solución de varios problemas, especialmente en el caso de condiciones de contorno específicas, como los que surgen en mecánica cuántica, fluidodinámica, electromagnetismo, el estudio del comportamiento del plasma, la difusión de sustancias químicas o el calor.

La principal ventaja de las coordenadas no cartesianas es que se pueden elegir para que coincidan con la simetría del problema. Por ejemplo, la onda de presión debido a una explosión lejos del suelo (u otras barreras) depende del espacio 3D en coordenadas cartesianas, sin embargo, la presión se aleja predominantemente del centro, de modo que en coordenadas esféricas el problema se vuelve casi unidimensional (ya que la onda de presión depende predominantemente solo del tiempo y de la distancia desde el centro). Otro ejemplo es el flujo (no turbulento) de un líquido en una tubería circular recta: en coordenadas cartesianas se tiene que resolver un difícil problema de valor límite bidimensional que involucra una ecuación diferencial parcial, pero en coordenadas cilíndricas el problema se vuelve unidimensional y se puede resolver mediante una ecuación diferencial ordinaria en lugar de una ecuación en derivadas parciales.

La razón para preferir coordenadas ortogonales en lugar de coordenadas curvilíneas en general es su simplicidad: surgen muchas complicaciones cuando las coordenadas no son ortogonales. Por ejemplo, en coordenadas ortogonales muchos problemas pueden resolverse utilizando separación de variables. La separación de variables es una técnica matemática que convierte un problema complejo de d dimensiones en d problemas unidimensionales que pueden resolverse en términos de funciones conocidas. Muchas ecuaciones pueden reducirse a la ecuación de Laplace o a la ecuación de Helmholtz. La ecuación de Laplace es separable en 13 de los 14Plantilla:Esdsistemas de coordenadas ortogonales mencionados más adelante (con la única excepción de las coordenadas toroidales), y la ecuación de Helmholtz es separable en 11 de los mencionados sistemas de coordenadas ortogonales.[1][2]

Las coordenadas ortogonales nunca tienen términos fuera de la diagonal en sus tensores métricos. En otras palabras, la distancia infinitesimal al cuadrado ds2 siempre se puede escribir como una suma escalada de los desplazamientos infinitesimales al cuadrado de las coordenadas

ds2=k=1d(hkdqk)2

donde d es la dimensión y las funciones de escala (o factores de escala)

hk(𝐪) =def gkk(𝐪)=|𝐞k|

son iguales a las raíces cuadradas de las componentes diagonales del tensor métrico, o las longitudes de los vectores base locales 𝐞k descritos a continuación. Estas funciones de escala hi se utilizan para calcular operadores diferenciales en las nuevas coordenadas, como por ejemplo el gradiente, el laplaciano, la divergencia y el rotacional.

Un método simple para generar sistemas de coordenadas ortogonales en dos dimensiones es mediante una transformación conforme de una cuadrícula bidimensional estándar de coordenadas cartesianas Plantilla:Nowrap. Se puede formar un número complejo z = x + iy a partir de las coordenadas reales x e y, donde i representa la unidad imaginaria. Cualquier función holomorfa w = f(z) con derivada compleja distinta de cero producirá una transformación conforme. Si el número complejo resultante se escribe Plantilla:Nowrap, entonces las curvas de u y v constantes se intersecan en ángulos rectos, tal como lo hacen las líneas originales de x e y constantes.

Se pueden generar coordenadas ortogonales en tres o más dimensiones a partir de un sistema de coordenadas ortogonal bidimensional, ya sea proyectándolo a una nueva dimensión (coordenadas cilíndricas) o rotando el sistema bidimensional sobre uno de sus ejes de simetría. Sin embargo, existen otros sistemas de coordenadas ortogonales en tres dimensiones que no se pueden obtener proyectando o rotando un sistema bidimensional, como las coordenadas elipsoidales. Se pueden obtener coordenadas ortogonales más generales comenzando con algunas superficies de coordenadas necesarias y considerando sus trayectorias ortogonales.

Definición

Dada una variedad de (pseudo)riemanniana , un conjunto abierto O del mismo y un punto dentro de dicho conjunto abierto mO, una carta local o "sistema de coordenadas" local puede representarse por una función: Plantilla:Ecuación Donde d es la dimensión del espacio donde se define el sistema de coordenadas local. Las d curvas coordenadas Ci(t) y sus vectores tangentes vienen definidas por las ecuaciones: Plantilla:Ecuación El sistema de coordenadas será ortogonal si los vectores tangentes a las curvas coordenadas xi son ortogonales, es decir, si: Plantilla:Ecuación Donde g(, ) es el tensor métrico del espacio donde se definen las coordenadas.

Vectores base

Base covariante

En coordenadas cartesianas, las bases son fijas (constantes). En el contexto más general de las coordenadas curvilíneas, un punto en el espacio se especifica mediante las coordenadas, y en cada uno de esos puntos hay un conjunto de vectores base, que generalmente no son constantes: esta es la esencia de las coordenadas curvilíneas en general y es un concepto muy importante. Lo que distingue a las coordenadas ortogonales es que, aunque los vectores base varían, siempre son ortogonales entre sí. En otras palabras,

𝐞i𝐞j=0siij

Estos vectores base son, por definición, los vectores tangentes de las curvas obtenidas al variar una coordenada, manteniendo las otras fijas:

Visualización de coordenadas ortogonales en 2D. Se muestran las curvas obtenidas manteniendo constantes todas las coordenadas menos una, junto con los vectores base. Debe tenerse en cuenta que los vectores base no tienen la misma longitud: no es necesario que la tengan, solo que sean ortogonales entre sí
𝐞i=𝐫qi

donde r es un punto y qi es la coordenada para la que se extrae el vector base. En otras palabras, se obtiene una curva al fijar todas las coordenadas menos una. La coordenada no fija se varía como en una ecuación paramétrica, y la derivada de la curva con respecto al parámetro (la coordenada variable) es el vector base para esa coordenada.

Debe tenerse en cuenta que los vectores no tienen necesariamente la misma longitud. Las funciones útiles conocidas como factores de escala de las coordenadas son simplemente las longitudes hi de los vectores base 𝐞i (consúltese la tabla siguiente). Los factores de escala a veces se denominan coeficientes de Lamé, que no deben confundirse con los parámetros de Lamé.

Los vectores de la base normalizados se anotan con un sombrero y se obtienen dividiendo por la longitud:

𝐞^i=𝐞ihi=𝐞i|𝐞i|

Un campo vectorial puede especificarse por sus componentes con respecto a los vectores base o los vectores base normalizados, y se debe estar seguro de qué caso se trata. Las componentes en la base normalizada son más comunes en aplicaciones para facilitar la claridad de las cantidades (por ejemplo, se puede querer tratar con la velocidad tangencial en lugar de con la velocidad tangencial multiplicada por un factor de escala); en las derivaciones, la base normalizada es menos común, ya que las expresiones son más complicadas.

Base contravariante

Los vectores de la base que se muestran arriba son vectores de la base covariante (porque covarían con los vectores). En el caso de coordenadas ortogonales, los vectores de la base contravariantes son fáciles de encontrar, ya que estarán en la misma dirección que los vectores covariantes pero tendrán la longitud recíproca (por esta razón, se dice que los dos conjuntos de vectores base son recíprocos entre sí):

𝐞i=𝐞^ihi=𝐞ihi2

Esto se deduce del hecho de que, por definición, 𝐞i𝐞j=δij, utilizando la delta de Kronecker. Nótese que:

𝐞^i=𝐞ihi=hi𝐞i𝐞^i

Por lo tanto, se cuenta con tres conjuntos de bases diferentes que se utilizan comúnmente para describir vectores en coordenadas ortogonales: la base covariante ei, la base contravariante ei y la base normalizada êi. Si bien un vector es una cantidad objetiva, es decir, su identidad es independiente de cualquier sistema de coordenadas, las componentes de un vector dependen de la base en la que se represente el vector.

Para evitar confusiones, las componentes del vector x con respecto a la base ei se representan como xi, mientras que las componentes con respecto a la base ei se representan como xi:

𝐱=ixi𝐞i=ixi𝐞i

La posición de los índices representa cómo se calculan las componentes (los índices superiores no deben confundirse con los exponentes de una potencia). Téngase en cuenta que los símbolos del sumatorio Σ (sigma mayúscula) y el rango de la suma, que indica la suma de todos los vectores base (i = 1, 2, ..., d), a menudo son omitidos. Las componentes están relacionadas simplemente por:

hi2xi=xi

No existe una notación general que distinga y se utilice para las componentes vectoriales con respecto a la base normalizada. En este artículo se utilizarán subíndices para las componentes vectoriales y se indicará que las componentes se calculan en la base normalizada.

Álgebra vectorial

La suma y el cambio de signo de vectores se realizan componente por componente, tal como en coordenadas cartesianas, sin complicaciones. Es posible que se requieran consideraciones adicionales para otras operaciones con vectores.

Sin embargo, debe tenerse en cuenta que todas estas operaciones suponen que dos vectores en un campo vectorial están ligados al mismo punto (en otras palabras, las colas de los vectores coinciden). Dado que los vectores base generalmente varían en coordenadas ortogonales, si se suman dos vectores cuyas componentes se calculan en diferentes puntos del espacio, es necesario considerar los diferentes vectores base.

Producto escalar

El producto escalar en coordenadas cartesianas (en el espacio euclídeo con una base ortonormal) es simplemente la suma de los productos de las componentes. En coordenadas ortogonales, el producto escalar de dos vectores x e y adopta esta forma familiar cuando las componentes de los vectores se calculan en la base normalizada:

𝐱𝐲=ixi𝐞^ijyj𝐞^j=ixiyi

Esta es una consecuencia inmediata del hecho de que la base normalizada en algún punto puede formar un sistema de coordenadas cartesianas: el conjunto de vectores de la base es ortonormal.

Para las componentes en las bases covariantes o contravariantes,

𝐱𝐲=ihi2xiyi=ixiyihi2=ixiyi=ixiyi

Esta expresión se puede deducir fácilmente escribiendo los vectores en forma de componentes, normalizando los vectores base y tomando el producto escalar. Por ejemplo, en 2D:

𝐱𝐲=(x1𝐞1+x2𝐞2)(y1𝐞1+y2𝐞2)=(x1h1𝐞^1+x2h2𝐞^2)(y1𝐞^1h1+y2𝐞^2h2)=x1y1+x2y2

donde se ha utilizado el hecho de que las bases covariantes y contravariantes normalizadas son iguales.

Producto vectorial

El producto vectorial en coordenadas cartesianas 3D es:

𝐱×𝐲=(x2y3x3y2)𝐞^1+(x3y1x1y3)𝐞^2+(x1y2x2y1)𝐞^3

La fórmula anterior sigue siendo válida en coordenadas ortogonales si las componentes se calculan en la base normalizada.

Para construir el producto vectorial en coordenadas ortogonales con bases covariantes o contravariantes, simplemente se deben normalizar los vectores de la base, por ejemplo:

𝐱×𝐲=ixi𝐞i×jyj𝐞j=ixihi𝐞^i×jyjhj𝐞^j

que, escrito expandido, toma la forma

𝐱×𝐲=(x2y3x3y2)h2h3h1𝐞1+(x3y1x1y3)h1h3h2𝐞2+(x1y2x2y1)h1h2h3𝐞3

La notación concisa para el producto vectorial, que simplifica la generalización a coordenadas no ortogonales y dimensiones superiores, es posible con el símbolo de Levi-Civita, que tendrá componentes distintas de ceros y unos si los factores de escala no son todos iguales a uno.

Cálculo vectorial

Diferenciación

Si se observa un desplazamiento infinitesimal desde algún punto, es evidente que

d𝐫=i𝐫qidqi=i𝐞idqi

Según su definición, el gradiente de una función debe satisfacer (esta definición sigue siendo cierta si ƒ es cualquier tensor)

df=fd𝐫df=fi𝐞idqi

De ello se deduce que el operador gradiente debe ser:

=i𝐞iqi

y esto sigue siendo cierto en coordenadas curvilíneas generales. Cantidades como el gradiente y el laplaciano se obtienen mediante la aplicación adecuada de este operador.

Fórmulas de vectores de una base

A partir de dr y los vectores de una base normalizados êi, se puede construir lo siguiente.[3][4]

Elemento diferencial Vectores Escalares
Elemento de línea Vector tangente a la curva coordenada qi:

d=hidqi𝐞^i=𝐫qidqi

Longitud infinitesimal

d=d𝐫d𝐫=(h1dq1)2+(h2dq2)2+(h3dq3)2

Elemento de superficie Normal a la superficie coordenada qk= constante:

d𝐒=(hidqi𝐞^i)×(hjdqj𝐞^j)=dqidqj(𝐫qi×𝐫qj)=hihjdqidqj𝐞^k

Superficie infinitesimal

dSk=hihjdqidqj

Elemento de volumen N/A Volumen infinitesimal

dV=|(h1dq1𝐞^1)(h2dq2𝐞^2)×(h3dq3𝐞^3)|=|𝐞^1𝐞^2×𝐞^3|h1h2h3dq1dq2dq3=h1h2h3dq1dq2dq3=Jdq1dq2dq3

donde

J=|𝐫q1(𝐫q2×𝐫q3)|=|(x,y,z)(q1,q2,q3)|=h1h2h3

es el determinante jacobiano, que tiene la interpretación geométrica de la deformación en volumen desde el cubo infinitesimal dxdydz hasta el volumen curvo infinitesimal en las coordenadas ortogonales.

Integración

Usando el elemento de línea que se muestra arriba, la integral de arco en una trayectoria 𝒫 de un vector F es:

𝒫𝐅d𝐫=𝒫iFi𝐞ij𝐞jdqj=i𝒫Fidqi

Un elemento infinitesimal de área para una superficie descrita manteniendo constante una coordenada qk es:

dAk=ikdsi=ikhidqi

De manera similar, el elemento de volumen es:

dV=idsi=ihidqi

donde el símbolo Π grande (pi mayúscula) indica un productorio, de la misma manera que una Σ grande indica un sumatorio. Nótese que el producto de todos los factores de escala es el determinante jacobiano.

A modo de ejemplo, la integral de superficie de una función vectorial F sobre una superficie q1 = constante 𝒮 en 3D es:

𝒮𝐅d𝐀=𝒮𝐅𝐧^ dA=𝒮𝐅𝐞^1 dA=𝒮F1h2h3h1dq2dq3

Debe tenerse en cuenta que F1/h1 es la componente de F normal a la superficie.

Propiedades

La elección de uno u otro sistema depende de las simetrías del problema geométrico o físico planteado. Al ser todos estos sistemas de coordenas ortogonales, el tensor métrico tiene en todos ellos la forma: Plantilla:Ecuación Donde las tres componentes no nulas son los llamados factores de escala son funciones de las tres coordenadas.

Operadores vectoriales en coordenadas ortogonales

Los operadores vectoriales pueden expresarse fácilmente en términos de estas componentes del tensor métrico.

Plantilla:Ecuación

Plantilla:Ecuación

  • El rotacional viene dado por el desarrollo del siguiente determinante:

Plantilla:Ecuación

  • El laplaciano de una magnitud escalar viene dado por:

Plantilla:Ecuación

Las expresiones anteriores se pueden escribir de forma más compacta utilizando la simbología de Levi-Civita ϵijk y el determinante jacobiano J=h1h2h3, suponiendo la suma sobre índices repetidos:

Plantilla:Ecuación

Plantilla:Ecuación

Plantilla:Ecuación

Plantilla:Ecuación

También debe observarse que el gradiente de un campo escalar se puede expresar en términos de la matriz jacobiana J que contiene las derivadas parciales canónicas:

𝐉=[ϕq1,ϕq2,ϕq3]

para un cambio de base:

ϕ=𝐒𝐑𝐉T

donde las matrices de rotación y escala son:

𝐑=[𝐞1,𝐞2,𝐞3]
𝐒=diag([h11,h21,h31]).

Ejemplos en el espacio euclídeo

En el espacio euclídeo tridimensional se emplean diferentes sistemas de coordenadas, a veces, combinando tipos de coordenadas ortogonales y angulares: Plantilla:Lista de columnas

Ejemplos en variedades diferenciales

La coordenadas usadas en la teoría de la relatividad general son el ejemplo físico más conocido de sistemas de coordenadas sobre un espacio globalmente no euclídeo.

En un espacio-tiempo estático siempre es posible escoger alrededor de cualquier punto del espacio-tiempo un sistema de coordenadas ortogonal.Plantilla:Cita requerida

Tabla de sistemas de coordenadas ortogonales bidimensionales

Sistema Transformación compleja

x+iy=f(u+iv)

Forma de las isolíneas u y v Comentario
Cartesiano u+iv Recta, recta
Log-polar exp(u+iv) Circunferencia, recta Para u=lnr equivalen a las coordenadas polares
Parabólico 12(u+iv)2 Parábola, parábola
Dipolo puntual (u+iv)1 Circunferencia, circunferencia
Elíptico cosh(u+iv) Elipse, hipérbola Similar al log-polar para distancias grandes
Bipolar coth(u+iv) Circunferencia, circunferencia Similar al dipolo puntual para distancias grandes
u+iv Hipérbola, hipérbola Campo en el interior de un límite
u=x2+2y2, y=vx2 Elipse, parábola

Plantilla:Multiple image

Tabla de coordenadas ortogonales tridimensionales

Además de las coordenadas cartesianas habituales, a continuación se tabulan varias otras.[4] Se utiliza la notación de intervalos para que la columna de coordenadas sea compacta: los paréntesis indican intervalos abiertos (es dedir, el valor al que acompañan es un límite del intervalo), y los corchetes indican intervalos cerrados (por el contrario, el valor al que acompañan forma parte del intervalo). Por ejemplo, en el caso de las coordenadas esféricas, los tres intervalos de los parámetros que las definen son CA-CC-CA (siendo C cerrado y A abierto).

Los tres factores de escala (h1h2h3) son los valores de la diagonal del tensor métrico que definepermite el cálculo de distintas operaciones dentro del sistema de coordenadas (véase la sección Propiedades).

Coordenadas curvilíneas (q1, q2, q3) Transformación a cartesianas (x, y, z) Factores de escala
Coordenadas esféricas

(r,θ,ϕ)[0,)×[0,π]×[0,2π)

x=rsinθcosϕy=rsinθsinϕz=rcosθ h1=1h2=rh3=rsinθ
Coordenadas parabólicas

(u,v,ϕ)[0,)×[0,)×[0,2π)

x=uvcosϕy=uvsinϕz=12(u2v2) h1=h2=u2+v2h3=uv
Coordenadas bipolares cilíndricas

(u,v,z)[0,2π)×(,)×(,)

x=asinhvcoshvcosuy=asinucoshvcosuz=z h1=h2=acoshvcosuh3=1
Coordenadas elipsoidales

(λ,μ,ν)[0,c2)×(c2,b2)×(b2,a2)λ<c2<b2<a2,c2<μ<b2<a2,c2<b2<ν<a2,

x2a2qi+y2b2qi+z2c2qi=1

donde (q1,q2,q3)=(λ,μ,ν)

hi=12(qjqi)(qkqi)(a2qi)(b2qi)(c2qi)
Coordenadas paraboloidales

(λ,μ,ν)[0,b2)×(b2,a2)×(a2,)b2<a2

x2qia2+y2qib2=2z+qi

donde (q1,q2,q3)=(λ,μ,ν)

hi=12(qjqi)(qkqi)(a2qi)(b2qi)
Coordenadas cilíndricas

(r,ϕ,z)[0,)×[0,2π)×(,)

x=rcosϕy=rsinϕz=z h1=h3=1h2=r
Coordenadas elípticas cilíndricas

(u,v,z)[0,)×[0,2π)×(,)

x=acoshucosvy=asinhusinvz=z h1=h2=asinh2u+sin2vh3=1
Coordenadas esferoidales oblatas

(ξ,η,ϕ)[0,)×[π2,π2]×[0,2π)

x=acoshξcosηcosϕy=acoshξcosηsinϕz=asinhξsinη h1=h2=asinh2ξ+sin2ηh3=acoshξcosη
Coordenadas esferoidales prolatas

(ξ,η,ϕ)[0,)×[0,π]×[0,2π)

x=asinhξsinηcosϕy=asinhξsinηsinϕz=acoshξcosη h1=h2=asinh2ξ+sin2ηh3=asinhξsinη
Coordenadas biesféricas

(u,v,ϕ)(π,π]×[0,)×[0,2π)

x=asinucosϕcoshvcosuy=asinusinϕcoshvcosuz=asinhvcoshvcosu h1=h2=acoshvcosuh3=asinucoshvcosu
Coordenadas toroidales

(u,v,ϕ)(π,π]×[0,)×[0,2π)

x=asinhvcosϕcoshvcosuy=asinhvsinϕcoshvcosuz=asinucoshvcosu h1=h2=acoshvcosuh3=asinhvcoshvcosu
Coordenadas parabólicas cilíndricas

(u,v,z)(,)×[0,)×(,)

x=12(u2v2)y=uvz=z h1=h2=u2+v2h3=1
Coordenadas cónicas

(λ,μ,ν)ν2<b2<μ2<a2λ[0,)

x=λμνaby=λa(μ2a2)(ν2a2)a2b2z=λb(μ2b2)(ν2b2)b2a2 h1=1h22=λ2(μ2ν2)(μ2a2)(b2μ2)h32=λ2(μ2ν2)(ν2a2)(ν2b2)

Referencias

Plantilla:Listaref

Plantilla:Control de autoridades

  1. Plantilla:Cite web
  2. Plantilla:Harvnb
  3. Mathematical Handbook of Formulas and Tables (3rd edition), S. Lipschutz, M.R. Spiegel, J. Liu, Schuam's Outline Series, 2009, Plantilla:Isbn.
  4. 4,0 4,1 Vector Analysis (2nd Edition), M.R. Spiegel, S. Lipschutz, D. Spellman, Schaum’s Outlines, McGraw Hill (USA), 2009, Plantilla:Isbn