En geometría , las coordenadas curvilíneas son un sistema de coordenadas para el espacio euclidiano en el que las líneas de coordenadas pueden ser curvas. Estas coordenadas pueden derivarse de un conjunto de coordenadas cartesianas mediante una transformación que es localmente invertible (una función biunívoca) en cada punto. Esto significa que se puede convertir un punto dado en un sistema de coordenadas cartesianas a sus coordenadas curvilíneas y viceversa. El nombre de coordenadas curvilíneas , acuñado por el matemático francés Lamé , deriva del hecho de que las superficies de coordenadas de los sistemas curvilíneos son curvas.
Ejemplos bien conocidos de sistemas de coordenadas curvilíneas en el espacio euclidiano tridimensional ( R 3 ) son las coordenadas cilíndricas y esféricas . Una superficie de coordenadas cartesianas en este espacio es un plano de coordenadas ; por ejemplo, z = 0 define el plano x - y . En el mismo espacio, la superficie de coordenadas r = 1 en coordenadas esféricas es la superficie de una esfera unitaria , que es curva. El formalismo de las coordenadas curvilíneas proporciona una descripción unificada y general de los sistemas de coordenadas estándar.
Las coordenadas curvilíneas se utilizan a menudo para definir la ubicación o distribución de magnitudes físicas que pueden ser, por ejemplo, escalares , vectores o tensores . Las expresiones matemáticas que involucran estas magnitudes en el cálculo vectorial y el análisis tensorial (como el gradiente , la divergencia , el rotacional y el laplaciano ) se pueden transformar de un sistema de coordenadas a otro, de acuerdo con las reglas de transformación para escalares, vectores y tensores. Dichas expresiones se vuelven entonces válidas para cualquier sistema de coordenadas curvilíneas.
Un sistema de coordenadas curvilíneas puede ser más sencillo de utilizar que el sistema de coordenadas cartesianas para algunas aplicaciones. El movimiento de partículas bajo la influencia de fuerzas centrales suele ser más fácil de resolver en coordenadas esféricas que en coordenadas cartesianas; esto es cierto para muchos problemas físicos con simetría esférica definida en R 3 . Las ecuaciones con condiciones de contorno que siguen superficies de coordenadas para un sistema de coordenadas curvilíneas particular pueden ser más fáciles de resolver en ese sistema. Si bien uno podría describir el movimiento de una partícula en una caja rectangular utilizando coordenadas cartesianas, es más fácil describir el movimiento en una esfera con coordenadas esféricas. Las coordenadas esféricas son los sistemas de coordenadas curvilíneas más comunes y se utilizan en ciencias de la Tierra , cartografía , mecánica cuántica , relatividad e ingeniería .
Por ahora, consideremos el espacio 3-D . Un punto P en el espacio 3-D (o su vector de posición r ) se puede definir usando coordenadas cartesianas ( x , y , z ) [escritas de manera equivalente ( x 1 , x 2 , x 3 )], por , donde e x , e y , e z son los vectores base estándar .
También se puede definir por sus coordenadas curvilíneas ( q 1 , q 2 , q 3 ) si este triplete de números define un único punto de forma inequívoca. La relación entre las coordenadas viene dada entonces por las funciones de transformación invertibles:
Las superficies q 1 = constante, q 2 = constante, q 3 = constante se denominan superficies coordenadas ; y las curvas espaciales formadas por su intersección en pares se denominan curvas coordenadas . Los ejes de coordenadas están determinados por las tangentes a las curvas coordenadas en la intersección de tres superficies. En general, no son direcciones fijas en el espacio, lo que sucede en el caso de las coordenadas cartesianas simples, y por lo tanto, generalmente no existe una base global natural para las coordenadas curvilíneas.
En el sistema cartesiano, los vectores base estándar se pueden derivar de la derivada de la ubicación del punto P con respecto a la coordenada local.
Aplicando las mismas derivadas al sistema curvilíneo localmente en el punto P se definen los vectores base naturales:
Una base de este tipo, cuyos vectores cambian su dirección y/o magnitud de un punto a otro, se denomina base local . Todas las bases asociadas a coordenadas curvilíneas son necesariamente locales. Los vectores de base que son iguales en todos los puntos son bases globales y pueden asociarse únicamente a sistemas de coordenadas lineales o afines .
Para este artículo, e está reservado para la base estándar (cartesiana) y h o b son para la base curvilínea.
Estos pueden no tener longitud unitaria y también pueden no ser ortogonales. En el caso de que sean ortogonales en todos los puntos donde las derivadas estén bien definidas, definimos los coeficientes de Lamé.(según Gabriel Lamé ) por
y los vectores base ortonormales curvilíneos por
Estos vectores base pueden depender de la posición de P ; por lo tanto, es necesario que no se suponga que son constantes en una región. (Técnicamente, forman una base para el fibrado tangente de en P , y por lo tanto son locales a P ).
En general, las coordenadas curvilíneas permiten que los vectores de base naturales h i no sean todos mutuamente perpendiculares entre sí, y no se requiere que tengan una longitud unitaria: pueden tener una magnitud y dirección arbitrarias. El uso de una base ortogonal hace que las manipulaciones de vectores sean más simples que para las no ortogonales. Sin embargo, algunas áreas de la física y la ingeniería , particularmente la mecánica de fluidos y la mecánica de medios continuos , requieren bases no ortogonales para describir deformaciones y transporte de fluidos para dar cuenta de dependencias direccionales complicadas de cantidades físicas. Más adelante en esta página aparece una discusión del caso general.
En coordenadas curvilíneas ortogonales, dado que el cambio diferencial total en r es
Entonces los factores de escala son
En coordenadas no ortogonales, la longitud de es la raíz cuadrada positiva de (con la convención de suma de Einstein ). Los seis productos escalares independientes g ij = h i . h j de los vectores base naturales generalizan los tres factores de escala definidos anteriormente para coordenadas ortogonales. Los nueve g ij son los componentes del tensor métrico , que tiene solo tres componentes distintos de cero en coordenadas ortogonales: g 11 = h 1 h 1 , g 22 = h 2 h 2 , g 33 = h 3 h 3 .
Los gradientes espaciales, las distancias, las derivadas temporales y los factores de escala están interrelacionados dentro de un sistema de coordenadas mediante dos grupos de vectores base:
Nótese que, debido a la convención de suma de Einstein, la posición de los índices de los vectores es opuesta a la de las coordenadas.
En consecuencia, un sistema de coordenadas curvilíneas general tiene dos conjuntos de vectores base para cada punto: { b 1 , b 2 , b 3 } es la base contravariante, y { b 1 , b 2 , b 3 } es la base covariante (también conocida como recíproca). Los tipos de vectores base covariantes y contravariantes tienen una dirección idéntica para los sistemas de coordenadas curvilíneas ortogonales, pero como es habitual, tienen unidades invertidas entre sí.
Tenga en cuenta la siguiente igualdad importante: donde denota el delta de Kronecker generalizado .
En el sistema de coordenadas cartesianas , podemos escribir el producto escalar como:
Consideremos un desplazamiento infinitesimal . Sean dq 1 , dq 2 y dq 3 los cambios infinitesimales correspondientes en las coordenadas curvilíneas q 1 , q 2 y q 3 respectivamente.
Por la regla de la cadena, dq 1 se puede expresar como:
Si el desplazamiento d r es tal que dq 2 = dq 3 = 0, es decir, el vector de posición r se mueve una cantidad infinitesimal a lo largo del eje de coordenadas q 2 =const y q 3 =const, entonces:
Dividiendo por dq 1 y tomando el límite dq 1 → 0:
o equivalentemente:
Ahora bien, si el desplazamiento d r es tal que dq 1 = dq 3 = 0, es decir, el vector de posición r se mueve una cantidad infinitesimal a lo largo del eje de coordenadas q 1 = constante y q 3 = constante, entonces:
Dividiendo por dq 2 , y tomando el límite dq 2 → 0:
o equivalentemente:
Y así sucesivamente para los demás productos escalares.
Prueba alternativa:
y la convención de suma de Einstein está implícita.
Un vector v se puede especificar en términos de cualquiera de las bases, es decir,
Utilizando la convención de suma de Einstein, los vectores base se relacionan con los componentes mediante [2] : 30–32
y
donde g es el tensor métrico (ver más abajo).
Un vector puede especificarse con coordenadas covariantes (índices reducidos, escritos v k ) o coordenadas contravariantes (índices elevados, escritos v k ). A partir de las sumas de vectores anteriores, se puede ver que las coordenadas contravariantes están asociadas con vectores base covariantes, y las coordenadas covariantes están asociadas con vectores base contravariantes.
Una característica clave de la representación de vectores y tensores en términos de componentes indexados y vectores base es la invariancia en el sentido de que los componentes vectoriales que se transforman de manera covariante (o contravariante) se emparejan con vectores base que se transforman de manera contravariante (o covariante).
Considere la curva unidimensional que se muestra en la Fig. 3. En el punto P , tomado como origen , x es una de las coordenadas cartesianas, y q 1 es una de las coordenadas curvilíneas. El vector base local (no unitario) es b 1 (anotado h 1 arriba, con b reservado para vectores unitarios) y está construido sobre el eje q 1 que es una tangente a esa línea de coordenadas en el punto P . El eje q 1 y por lo tanto el vector b 1 forman un ángulo con el eje cartesiano x y el vector base cartesiano e 1 .
Se puede ver desde el triángulo PAB que
donde | e 1 |, | b 1 | son las magnitudes de los dos vectores base, es decir, las intersecciones escalares PB y PA . PA es también la proyección de b 1 sobre el eje x .
Sin embargo, este método para transformaciones de vectores base utilizando cosenos direccionales no es aplicable a coordenadas curvilíneas por las siguientes razones:
Los ángulos que la línea q 1 y ese eje forman con el eje x se acercan en valor cuanto más nos acercamos al punto P y se vuelven exactamente iguales en P .
Sea el punto E situado muy cerca de P , tan cerca que la distancia PE sea infinitesimalmente pequeña. Entonces, PE medido en el eje q 1 casi coincide con PE medido en la línea q 1 . Al mismo tiempo, la relación PD/PE ( siendo PD la proyección de PE en el eje x ) se vuelve casi exactamente igual a .
Sean los puntos de corte infinitesimalmente pequeños PD y PE etiquetados, respectivamente, como dx y d q 1 . Entonces
De este modo, los cosenos direccionales pueden sustituirse en las transformaciones con las razones más exactas entre los puntos de corte de coordenadas infinitesimalmente pequeños. De ello se deduce que el componente (proyección) de b 1 sobre el eje x es
Si q i = q i ( x 1 , x 2 , x 3 ) y x i = x i ( q 1 , q 2 , q 3 ) son funciones suaves (continuamente diferenciables), las razones de transformación se pueden escribir como y . Es decir, esas razones son derivadas parciales de coordenadas pertenecientes a un sistema con respecto a coordenadas pertenecientes al otro sistema.
Haciendo lo mismo para las coordenadas en las otras 2 dimensiones, b 1 se puede expresar como:
Se cumplen ecuaciones similares para b 2 y b 3 de modo que la base estándar { e 1 , e 2 , e 3 } se transforma en una base local (ordenada y normalizada ) { b 1 , b 2 , b 3 } mediante el siguiente sistema de ecuaciones:
Por razonamiento análogo, se puede obtener la transformación inversa de base local a base estándar:
Los sistemas de ecuaciones lineales anteriores se pueden escribir en forma matricial utilizando la convención de suma de Einstein como
Esta matriz de coeficientes del sistema lineal es la matriz jacobiana (y su inversa) de la transformación. Estas son las ecuaciones que se pueden utilizar para transformar una base cartesiana en una base curvilínea y viceversa.
En tres dimensiones, las formas expandidas de estas matrices son
En la transformación inversa (segundo sistema de ecuaciones), las incógnitas son los vectores base curvilíneos. Para cualquier posición específica solo puede existir un único conjunto de vectores base (de lo contrario, la base no está bien definida en ese punto). Esta condición se cumple si y solo si el sistema de ecuaciones tiene una única solución. En álgebra lineal , un sistema de ecuaciones lineales tiene una única solución (no trivial) solo si el determinante de su matriz de sistema no es cero:
lo que muestra la razón detrás del requisito anterior relativo al determinante jacobiano inverso.
El formalismo se extiende a cualquier dimensión finita de la siguiente manera.
Consideremos el espacio euclidiano real n -dimensional, es decir R n = R × R × ... × R ( n veces) donde R es el conjunto de números reales y × denota el producto cartesiano , que es un espacio vectorial .
Las coordenadas de este espacio se pueden denotar por: x = ( x 1 , x 2 ,..., x n ). Como se trata de un vector (un elemento del espacio vectorial), se puede escribir como:
donde e 1 = (1,0,0...,0), e 2 = (0,1,0...,0), e 3 = (0,0,1...,0),..., e n = (0,0,0...,1) es el conjunto base estándar de vectores para el espacio R n , e i = 1, 2,... n es un índice que etiqueta los componentes. Cada vector tiene exactamente un componente en cada dimensión (o "eje") y son mutuamente ortogonales ( perpendiculares ) y normalizados (tienen magnitud unitaria ).
En términos más generales, podemos definir vectores base b i de modo que dependan de q = ( q 1 , q 2 ,..., q n ), es decir, que cambien de un punto a otro: b i = b i ( q ). En cuyo caso, para definir el mismo punto x en términos de esta base alternativa: las coordenadas con respecto a esta base v i también dependen necesariamente de x , es decir, v i = v i ( x ). Entonces, un vector v en este espacio, con respecto a estas coordenadas alternativas y vectores base, se puede desarrollar como una combinación lineal en esta base (lo que simplemente significa multiplicar cada vector base e i por un número v i – multiplicación escalar ):
La suma vectorial que describe v en la nueva base está compuesta de diferentes vectores, aunque la suma en sí sigue siendo la misma.
Desde una perspectiva más general y abstracta, un sistema de coordenadas curvilíneas es simplemente un parche de coordenadas en la variedad diferenciable E n ( espacio euclidiano n-dimensional ) que es difeomorfo al parche de coordenadas cartesianas en la variedad. [3] Dos parches de coordenadas difeomorfos en una variedad diferencial no necesitan superponerse de manera diferenciable. Con esta definición simple de un sistema de coordenadas curvilíneas, todos los resultados que siguen a continuación son simplemente aplicaciones de teoremas estándar en topología diferencial .
Las funciones de transformación son tales que existe una relación uno a uno entre los puntos en las coordenadas "antiguas" y "nuevas", es decir, esas funciones son biyecciones y cumplen los siguientes requisitos dentro de sus dominios :
El álgebra elemental de vectores y tensores en coordenadas curvilíneas se utiliza en algunas de las publicaciones científicas más antiguas sobre mecánica y física y puede ser indispensable para comprender trabajos de principios y mediados del siglo XX, por ejemplo, el texto de Green y Zerna. [5] En esta sección se dan algunas relaciones útiles en el álgebra de vectores y tensores de segundo orden en coordenadas curvilíneas. La notación y el contenido provienen principalmente de Ogden, [6] Naghdi, [7] Simmonds, [2] Green y Zerna, [5] Basar y Weichert, [8] y Ciarlet. [9]
Un tensor de segundo orden se puede expresar como
donde denota el producto tensorial . Los componentes S ij se denominan componentes contravariantes , S i j componentes mixtos covariantes por la derecha , S i j componentes mixtos covariantes por la izquierda y S ij componentes covariantes del tensor de segundo orden. Los componentes del tensor de segundo orden están relacionados por
En cada punto, se puede construir un pequeño elemento de línea d x , por lo que el cuadrado de la longitud del elemento de línea es el producto escalar d x • d x y se llama métrica del espacio , dado por:
La siguiente parte de la ecuación anterior
es un tensor simétrico llamado tensor fundamental (o métrico) del espacio euclidiano en coordenadas curvilíneas.
Los índices se pueden aumentar o disminuir mediante la métrica:
Definiendo los factores de escala h i mediante
da una relación entre el tensor métrico y los coeficientes de Lamé, y
donde h ij son los coeficientes de Lamé. Para una base ortogonal también tenemos:
Si consideramos coordenadas polares para R 2 ,
(r, θ) son las coordenadas curvilíneas, y el determinante jacobiano de la transformación ( r ,θ) → ( r cos θ, r sin θ) es r .
Los vectores de base ortogonales son b r = (cos θ, sin θ), b θ = (−r sin θ, r cos θ). Los factores de escala son h r = 1 y h θ = r . El tensor fundamental es g 11 =1, g 22 = r 2 , g 12 = g 21 =0.
En una base dextrógira ortonormal, el tensor alterno de tercer orden se define como
En una base curvilínea general el mismo tensor puede expresarse como
También se puede demostrar que
donde la coma denota una derivada parcial (ver Cálculo de Ricci ). Para expresar Γ kij en términos de g ij ,
Desde
Al utilizarlos para reorganizar las relaciones anteriores se obtiene
Esto implica que
Otras relaciones que siguen son
El producto escalar de dos vectores en coordenadas curvilíneas es [2] : 32
El producto vectorial de dos vectores se expresa mediante [2] : 32–34
donde es el símbolo de permutación y es un vector de base cartesiana. En coordenadas curvilíneas, la expresión equivalente es
Es necesario realizar ajustes en el cálculo de las integrales de línea , superficie y volumen . Para simplificar, lo siguiente se limita a tres dimensiones y coordenadas curvilíneas ortogonales. Sin embargo, los mismos argumentos se aplican para espacios n -dimensionales. Cuando el sistema de coordenadas no es ortogonal, hay algunos términos adicionales en las expresiones.
Simmonds, [2] en su libro sobre análisis tensorial , cita a Albert Einstein diciendo [10]
La magia de esta teoría difícilmente dejará de imponerse a cualquiera que la haya comprendido verdaderamente; representa un auténtico triunfo del método del cálculo diferencial absoluto, fundado por Gauss, Riemann, Ricci y Lévi-Civita.
El cálculo vectorial y tensorial en coordenadas curvilíneas generales se utiliza en el análisis tensorial en variedades curvilíneas de cuatro dimensiones en la relatividad general , [11] en la mecánica de capas curvas , [9] en el examen de las propiedades de invariancia de las ecuaciones de Maxwell , lo que ha sido de interés en metamateriales [12] [13] y en muchos otros campos.
En esta sección se presentan algunas relaciones útiles en el cálculo de vectores y tensores de segundo orden en coordenadas curvilíneas. La notación y el contenido provienen principalmente de Ogden, [14] Simmonds, [2] Green y Zerna, [5] Basar y Weichert, [8] y Ciarlet. [9]
Sea φ = φ( x ) un campo escalar bien definido y v = v ( x ) un campo vectorial bien definido, y λ 1 , λ 2 ... parámetros de las coordenadas
es un vector tangente a C en coordenadas curvilíneas (usando la regla de la cadena ). Usando la definición de los coeficientes de Lamé, y la de la métrica g ij = 0 cuando i ≠ j , la magnitud es:
donde es el símbolo de permutación . En forma determinante:
Las expresiones para el gradiente, la divergencia y el laplaciano se pueden extender directamente a n dimensiones, sin embargo el rizo solo se define en 3D.
El campo vectorial b i es tangente a la curva de coordenadas q i y forma una base natural en cada punto de la curva. Esta base, como se explicó al principio de este artículo, también se denomina base curvilínea covariante . También podemos definir una base recíproca o base curvilínea contravariante , b i . Todas las relaciones algebraicas entre los vectores base, como se explicó en la sección sobre álgebra tensorial, se aplican a la base natural y su recíproco en cada punto x .
Por definición, si una partícula sin fuerzas actuando sobre ella tiene su posición expresada en un sistema de coordenadas inercial, ( x 1 , x 2 , x 3 , t ), entonces allí no tendrá aceleración (d 2 x j /d t 2 = 0). [15] En este contexto, un sistema de coordenadas puede dejar de ser "inercial" ya sea debido a que el eje del tiempo no es recto o a que el eje del espacio no es recto (o ambos). En otras palabras, los vectores base de las coordenadas pueden variar en el tiempo en posiciones fijas, o pueden variar con la posición en tiempos fijos, o ambos. Cuando las ecuaciones de movimiento se expresan en términos de cualquier sistema de coordenadas no inercial (en este sentido), aparecen términos adicionales, llamados símbolos de Christoffel. Estrictamente hablando, estos términos representan componentes de la aceleración absoluta (en mecánica clásica), pero también podemos optar por seguir considerando d 2 x j /d t 2 como la aceleración (como si las coordenadas fueran inerciales) y tratar los términos adicionales como si fueran fuerzas, en cuyo caso se denominan fuerzas ficticias. [16] El componente de cualquier fuerza ficticia normal a la trayectoria de la partícula y en el plano de curvatura de la trayectoria se denomina entonces fuerza centrífuga . [17]
Este contexto más general deja en claro la correspondencia entre los conceptos de fuerza centrífuga en sistemas de coordenadas rotatorios y en sistemas de coordenadas curvilíneas estacionarias. (Ambos conceptos aparecen frecuentemente en la literatura. [18] [19] [20] ) Para un ejemplo simple, considere una partícula de masa m moviéndose en un círculo de radio r con velocidad angular w relativa a un sistema de coordenadas polares que rota con velocidad angular W . La ecuación radial de movimiento es mr ” = F r + mr ( w + W ) 2 . Por lo tanto, la fuerza centrífuga es mr por el cuadrado de la velocidad rotacional absoluta A = w + W de la partícula. Si elegimos un sistema de coordenadas que rota a la velocidad de la partícula, entonces W = A y w = 0, en cuyo caso la fuerza centrífuga es mrA 2 , mientras que si elegimos un sistema de coordenadas estacionario tenemos W = 0 y w = A , en cuyo caso la fuerza centrífuga es nuevamente mrA 2 . La razón de esta igualdad de resultados es que en ambos casos los vectores base en la posición de la partícula cambian en el tiempo exactamente de la misma manera. Por lo tanto, en realidad son sólo dos formas diferentes de describir exactamente lo mismo, una descripción en términos de coordenadas rotatorias y la otra en términos de coordenadas curvilíneas estacionarias, ambas no inerciales según el significado más abstracto de ese término.
Al describir el movimiento general, las fuerzas reales que actúan sobre una partícula se refieren a menudo al círculo osculador instantáneo tangente a la trayectoria del movimiento, y este círculo en el caso general no está centrado en una ubicación fija, por lo que la descomposición en componentes centrífugos y de Coriolis cambia constantemente. Esto es así independientemente de si el movimiento se describe en términos de coordenadas estacionarias o rotatorias.