Coordenadas curvilíneas ortogonales en 3 dimensiones [ editar ]
Coordenadas, base, y los vectores [ editar ]
Fig. 1 - Coordinar superficies, líneas de coordenadas y ejes de coordenadas de coordenadas curvilíneas generales.
Fig. 2 - Coordinar superficies, coordinar líneas y coordinar ejes de coordenadas esféricas. Superficies: r - esferas, θ - conos, - semiplanos; Líneas: r - vigas rectas, θ - semicírculos verticales, φ - círculos horizontales; Ejes: r - vigas rectas, - tangentes a semicírculos verticales, φ - tangentes a círculos horizontales
Por ahora, consideremos el espacio 3-D . Un punto P en el espacio 3d (o su vector de posición r ) se puede definir usando coordenadas cartesianas ( x , y , z ) [escrito de manera equivalente ( x 1 , x 2 , x 3)], por, donde e x , e y , e zson los vectores de 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 solo punto de manera inequívoca. La relación entre las coordenadas viene dada por las funciones de transformación invertible:
Las superficies q 1 = constante, q 2 = constante, q 3 = constante se denominan superficies de coordenadas ; y las curvas de espacio formadas por su intersección en pares se denominan curvas de coordenadas . Los ejes de coordenadas están determinados por las tangentes a las curvas de coordenadas en la intersección de tres superficies. No son, en general, direcciones fijas en el espacio, lo que sucede en el caso de las coordenadas cartesianas simples, por lo que generalmente no existe una base global natural para las coordenadas curvilíneas.
En el sistema cartesiano, los vectores de base estándar pueden derivarse de la derivada de la ubicación del punto P con respecto a la coordenada local
La aplicación de los mismos derivados al sistema curvilíneo localmente en el punto P define los vectores de base natural:
Dicha base, cuyos vectores cambian su dirección y / o magnitud de un punto a otro se denomina base local . Todas las bases asociadas con coordenadas curvilíneas son necesariamente locales. Los vectores de base que son iguales en todos los puntos son bases globales , y pueden asociarse solo con sistemas de coordenadaslineales o afines .
Nota: para este artículo, e está reservado para la base estándar (cartesiana) y h o b es para la base curvilínea.
Estos pueden no tener una unidad de longitud, y también pueden no ser ortogonales. En el caso de que seanortogonales en todos los puntos donde los derivados están bien definidos, definimos los coeficientes de Lamé(después de Gabriel Lamé ) por
y los vectores de base ortonormal curvilínea por
Es importante tener en cuenta que estos vectores de base pueden depender de la posición de P ; por lo tanto, es necesario que no se asuma que son constantes en una región. (Técnicamente forman una base para el paquete tangente deen P , y por lo tanto son locales a P. )
En general, coordenadas curvilíneas permiten que los vectores de la base natural de h i no todos mutuamente perpendiculares entre sí, y no se requiere para ser de unidad de longitud: pueden ser de magnitud y dirección arbitraria. El uso de una base ortogonal hace que las manipulaciones de vectores sean más simples que las no ortogonales. Sin embargo, algunas áreas de la física y la ingeniería , en particular la mecánica de fluidos y la mecánica continua , requieren bases no ortogonales para describir las deformaciones y el transporte de fluidos para tener en cuenta las dependencias direccionales complicadas de las cantidades físicas. Una discusión del caso general aparece más adelante en esta página.
Cálculo vectorial [ editar ]
Elementos diferenciales [ editar ]
entonces los factores de escala son
En coordenadas no ortogonales la longitud de es la raíz cuadrada positiva de (Con la convención sumatoria de Einstein ). Los seis productos escalares independientes g ij = h i . h j de los vectores de base natural 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 h3 .
Covariantes y contravariantes bases [ editar ]
Un vector v ( rojo ) representado por • una base vectorial ( amarillo, izquierda: e 1 , e 2 , e 3 ), vectores tangentes para coordinar curvas ( negro ) y • una base covector o cobasis ( azul , derecha: e 1 , e 2 , e 3), vectores normales para coordinar superficies ( gris ) en coordenadas curvilíneas generales (no necesariamente ortogonales ) ( q 1 , q 2 , q3 ). Tenga en cuenta la base y las cobasis no coinciden a menos que el sistema de coordenadas sea ortogonal. [1]
Los gradientes espaciales, las distancias, las derivadas de tiempo y los factores de escala están interrelacionados dentro de un sistema de coordenadas por dos grupos de vectores básicos:
- Vectores de base que son localmente tangentes a su línea de ruta de coordenadas asociada:
que se transforma como vectores covariantes (denotados por índices reducidos), o
- Vectores de base que son localmente normales a la isosuperficie creada por las otras coordenadas:
que se transforma como vectores contravariantes (denotados por índices elevados), ∇ es el operador del .
En consecuencia, un sistema de coordenadas curvilíneas general tiene dos conjuntos de vectores de base para cada punto: { b 1 , b 2 , b 3 } es la base covariante, y { b 1 , b 2 , b 3 } es la contravariante (también conocida como recíproca) base. Los tipos de vectores de base covariante y contravariante 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:
-
Un vector v puede especificarse en términos de cualquier base, es decir,
Usando la convención de suma de Einstein, los vectores de base se relacionan con los componentes por [2] ( pp30–32 )
y
donde g es el tensor métrico (ver más abajo).
Un vector se puede especificar con coordenadas covariantes (índices bajos, escrito v k ) o coordenadas contravariantes (índices elevados, escrito v k ). De las sumas de vectores anteriores, se puede ver que las coordenadas contravariantes están asociadas con vectores de base covariantes, y las coordenadas covariantes están asociadas con vectores de base contravariantes.
Una característica clave de la representación de vectores y tensores en términos de componentes indexados y vectores básicos es la invariabilidad en el sentido de que los componentes vectoriales que se transforman de manera covariante (o contravariante) están emparejados con vectores de base que se transforman en contravariante (o de manera covariante).
Base covariante [ editar ]
Construyendo una base covariante en una dimensión [ editar ]
Fig. 3 - Transformación de bases covariantes locales en el caso de coordenadas curvilíneas generales
Considere la curva unidimensional que se muestra en la Fig. 3. En el punto P , tomada como origen , xes una de las coordenadas cartesianas, y q 1 es una de las coordenadas curvilíneas. El (no unidad) base vector local es b 1 (anotado h 1 anterior, con breservado para vectores unitarios) y que está construido sobre la q 1 eje que es una tangente a la recta de coordenadas en el punto P . El eje q 1 y por tanto el vector b 1 forman un ángulo.con el eje xcartesiano y el vector de base cartesiana e 1 .
Se puede ver desde el triángulo PAB que
donde | e 1 |, | b 1 | son las magnitudes de los dos vectores de base, es decir, las intercepciones escalares PB y PA . Tenga en cuenta que PA también es la proyección de b 1 en el eje x .
Sin embargo, este método para transformaciones de vectores de base que utilizan cosenos direccionales no es aplicable a las coordenadas curvilíneas por las siguientes razones:
- Al aumentar la distancia desde P , el ángulo entre la línea curva q 1 y el eje cartesiano x se desvía cada vez más de.
- A la distancia PB, el ángulo verdadero es el que forma la tangente en el punto C con el eje x, y este último ángulo es claramente diferente de.
Los ángulos que el q 1 línea y que forma eje con los x eje se vuelven más cerca en valor cuanto más cerca se mueve hacia el punto P y se convierten en exactamente igual a P .
Deje que el punto E se ubique muy cerca de P , tan cerca que la distancia PE sea infinitamente pequeña. Entonces PE mide en los q 1 eje casi coincide con PE medida en el q 1 línea. Al mismo tiempo, la relación PD / PE ( PD es la proyección de PE en el eje x ) se vuelve casi exactamente igual a.
Deje que las intercepciones infinitesimalmente pequeñas, PD y PE, se etiqueten, respectivamente, como dx y d q 1 . Entonces
- .
Por lo tanto, los cosenos direccionales se pueden sustituir en transformaciones con las relaciones más exactas entre las intercepciones de coordenadas infinitesimalmente pequeñas. De ello se deduce que la componente (proyección) de b 1 en 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 relaciones de transformación se pueden escribir como y . Es decir, esas relaciones son derivadas parciales de coordenadas que pertenecen a un sistema con respecto a las coordenadas que pertenecen a otro sistema.
Construyendo una base covariante en tres dimensiones [ editar ]
Haciendo lo mismo para las coordenadas en las otras 2 dimensiones, b 1 se puede expresar como:
Ecuaciones similares son válidas 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 } por el siguiente sistema de ecuaciones
Por un razonamiento análogo, se puede obtener la transformación inversa de una base local a otra estándar:
Jacobiano de la transformación [ editar ]
- .
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 usar 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 (sistema de segunda ecuación), las incógnitas son los vectores de base curvilínea. Para cualquier ubicación específica, solo puede existir uno y solo un 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 solución única, del álgebra lineal , un sistema de ecuaciones lineales tiene una solución única (no trivial) solo si el determinante de su matriz de sistema no es cero:
que muestra la razón detrás del requisito anterior con respecto al determinante jacobiano inverso.
Generalización a n dimensiones [ editar ]
El formalismo se extiende a cualquier dimensión finita de la siguiente manera.
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 la norma base conjunto de vectores para el espacio R n , y i = 1, 2, ... n es un componentes de etiquetado índice. Cada vector tiene exactamente un componente en cada dimensión (o "eje") y son mutuamente ortogonales ( perpendiculares ) y normalizados (tienen una unidad de magnitud ).
Más generalmente, podemos definir los vectores de base b i para que dependan de q = ( q 1 , q 2 , ..., q n ), es decir, cambian de punto a punto: 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 también, es decir v i = v i ( x ). Entonces un vectorv en este espacio, con respecto a estas coordenadas alternativas y los vectores de base, se puede expandir como una combinación lineal en esta base (lo que simplemente significa multiplicar cada vector de base e i por un número v i - multiplicación escalar ):
La suma de vectores que describe v en la nueva base se compone de diferentes vectores, aunque la suma en sí permanece igual.
Transformación de coordenadas [ editar ]
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] Tenga en cuenta que dos parches de coordenadas difeomórficas en un colector diferencial no necesitan solaparse de manera diferente. 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 de uno a uno entre los puntos en las coordenadas "antiguo" y "nuevo", es decir, esas funciones son bijections y cumplen los siguientes requisitos dentro de sus dominios :
- Son funciones suaves : q i = q i ( x )
- El determinante jacobiano inverso.
no es cero; Es decir, la transformación es invertible : x i ( q ).
De acuerdo con el teorema de la función inversa . La condición de que el determinante jacobiano no sea cero refleja el hecho de que tres superficies de diferentes familias se intersecan en un solo punto y, por lo tanto, determinan la posición de este punto de una manera única. [4]
Álgebra vectorial y tensorial en coordenadas curvilíneas tridimensionales [ editar ]
- Nota: la convención de suma de Einstein de sumar en índices repetidos se utiliza a continuación.
El vector elemental y el álgebra tensorial en coordenadas curvilíneas se utilizan en la literatura científica más antigua en mecánica y física y pueden ser indispensables para comprender el trabajo desde 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 los contenidos son principalmente de Ogden, [6] Naghdi, [7] Simmonds, [2] Green y Zerna, [5] Basar y Weichert, [8] y Ciarlet. [9]
Tensores en coordenadas curvilíneas [ editar ]
Un tensor de segundo orden se puede expresar como
dónde Denota el producto tensorial . Los componentes S ij se denominan componentes contravariantes , S i jlos componentes mixtos covariantes derechos , S i j los componentes mixtos covariantes izquierdos y S ij los componentes covariantes del tensor de segundo orden. Los componentes del tensor de segundo orden están relacionados por
El tensor métrico en coordenadas curvilíneas ortogonales [ editar ]
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 la métrica del espacio , dado por:
- .
La siguiente porción de la ecuación anterior.
Relación con los coeficientes de Lamé [ editar ]
Definiendo los factores de escala h i por
Da una relación entre el tensor métrico y los coeficientes de Lamé. Tenga en cuenta también que
donde h ij son los coeficientes de lamé. Para una base ortogonal también tenemos:
Ejemplo: coordenadas polares [ editar ]
Si consideramos coordenadas polares para R 2 , ten en cuenta que
(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 ortogonal son b r = (cos θ, sin θ), b θ = (−sin θ, 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.
El tensor alterno [ editar ]
En una base ortonormal con la mano derecha, 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
Símbolos de Christoffel [ editar ]
- Símbolos de Christoffel del primer tipo.
Ya que
el uso de estos para reorganizar las relaciones anteriores da
- Símbolos de Christoffel del segundo tipo.
Esto implica que
Otras relaciones que siguen son
Operaciones vectoriales [ editar ]
- Producto de punto :
El producto escalar de dos vectores en coordenadas curvilíneas es [2] ( p32 )
- Producto vectorial :
dónde es el símbolo de permutación yEs un vector de base cartesiana. En coordenadas curvilíneas, la expresión equivalente es
dónde Es el tensor alterno de tercer orden .
Cálculo vectorial y tensorial en coordenadas curvilíneas tridimensionales [ editar ]
- Nota: la convención de suma de Einstein de sumar en índices repetidos se utiliza a continuación.
Los ajustes deben realizarse en el cálculo de las integrales de línea , superficie y volumen . Para simplificar, lo siguiente se restringe a tres dimensiones y coordenadas curvilíneas ortogonales. Sin embargo, los mismos argumentos se aplican a los espacios en n dimensiones. Cuando el sistema de coordenadas no es ortogonal, hay algunos términos adicionales en las expresiones.
La magia de esta teoría no dejará de imponerse a quien la haya entendido realmente; representa un auténtico triunfo del método de cálculo diferencial absoluto, fundado por Gauss, Riemann, Ricci y Levi-Civita.
En esta sección se dan algunas relaciones útiles en el cálculo de vectores y tensores de segundo orden en coordenadas curvilíneas. La notación y los contenidos son 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 ... sean parámetros de las coordenadas
Elementos geométricos [ editar ]
- Vector tangente : si x ( λ ) parametriza una curva C en coordenadas cartesianas, entonces
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 eso para la métrica g ij = 0 cuando i ≠ j , la magnitud es:
- Elemento plano tangente : si x ( λ 1 , λ 2 ) parametriza una superficie S en coordenadas cartesianas, el siguiente producto cruzado de vectores tangentes es un vector normal a S con la magnitud del elemento plano infinitesimal, en coordenadas curvilíneas. Utilizando el resultado anterior,
Integración [ editar ]
Operador | Campo escalar | Campo vectorial |
Integral de linea | | |
Integral de superficie | | |
Volumen integral | | |
La diferenciación [ editar ]
Las expresiones para gradiente, divergencia y laplaciano pueden extenderse 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 discutió al principio de este artículo, también se denomina base curvilínea covariante. También podemos definir una base recíproca , o una base curvilínea contravariante , b i . Todas las relaciones algebraicas entre los vectores de base, como se discutió en la sección sobre álgebra tensorial, se aplican para la base natural y su recíproco en cada punto x .
-
Fuerzas ficticias en coordenadas curvilíneas generales [ editar ]
Por definición, si una partícula sin fuerzas que actúan sobre ella tiene su posición expresada en un sistema de coordenadas de inercia ( x 1 , x 2 , x 3 , t ), entonces 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" debido a un eje de tiempo no recto o a ejes espaciales no rectos (o ambos). En otras palabras, los vectores de base de las coordenadas pueden variar en el tiempo en posiciones fijas, o pueden variar con la posición en tiempos fijos, o en 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 elegir continuar 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 llaman fuerzas ficticias.[16] La componente de cualquier fuerza ficticia normal a la trayectoria de la partícula y en el plano de la curvatura de la trayectoria se llama fuerza centrífuga . [17]
Este contexto más general aclara la correspondencia entre los conceptos de fuerza centrífuga en sistemas de coordenadas de rotación y en sistemas de coordenadas curvilíneas estacionarias. (Ambos de estos conceptos aparecen con frecuencia en la literatura. [18] [19] [20] ) Para un ejemplo simple, considere una partícula de masa m que se mueve en un círculo de radio r con velocidad angular w relativa a un sistema de coordenadas polares que gira con velocidad angular W . La ecuación radial del movimiento es mr ”= F r + mr ( w + W ) 2. Así, la fuerza centrífuga es mr veces el cuadrado de la velocidad de rotación absoluta A = w + W de la partícula. Si elegimos un sistema de coordenadas que gira 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 seleccionamos un sistema de coordenadas estacionario tenemos W = 0 y w = A , en cuyo caso la fuerza centrífuga vuelve a ser mrA 2.. El motivo de esta igualdad de resultados es que en ambos casos los vectores base en la ubicación de la partícula están cambiando en el tiempo exactamente de la misma manera. Por lo tanto, estas son realmente dos formas diferentes de describir exactamente lo mismo, una descripción en términos de coordenadas rotativas y la otra en términos de coordenadas curvilíneas estacionarias, las cuales no son inerciales según el significado más abstracto de ese término .
Cuando se describe el movimiento general, las fuerzas reales que actúan sobre una partícula a menudo se refieren al círculo oscilante 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 centrífuga y Coriolis Los componentes están cambiando constantemente. Esto es cierto independientemente de si el movimiento se describe en términos de coordenadas estacionarias o giratorias.
No hay comentarios:
Publicar un comentario