DEPARTAMENTO DE FÍSICA MATEMÁTICA Y DE FLUIDOS

Examen de Métodos Matemáticos IV

3er curso de Grado en Fı́sica
(códigos: carrera 6104, asignatura 61043012)
CURSO 2014-2015 Prueba de Evaluación Continua

PRUEBA DE EVALUACIÓN CONTINUA (10 puntos)

Estudio de superficies generales de revolución en el espacio tridimensional Euclı́deo habitual:
Consideremos una curva genérica (γ(s)) contenida en el plano xz definida por las ecuaciones paramétricas

x = f (s), z = g(s)

siendo s el parámetro longitud de arco (es decir, por hipótesis se cumple f˙2 + ġ 2 = 1 ∀s [f˙ ≡ df /ds]) donde s
toma valores en un cierto intervalo I. Suponiendo que esta curva no corta en ningún punto al eje z (es decir,
suponiendo que f (s) > 0 ∀s ∈ I) consideremos la superficie de revolución que se obtiene al hacer girar dicha
curva paramétrica alrededor del eje z. Dicha superficie está descrita en coordenadas paramétricas por medio
de las relaciones
x = f (s) cos θ, y = f (s) sen θ, z = g(s)
con s ∈ I, θ ∈ [0, 2π].
En este ejercicio vamos realizar un estudio de las propiedades generales de las superficies de revolución
regulares parametrizadas en términos de los parámetros

y 1 = s ∈ I, y 2 = θ ∈ [0, 2π]

suponemos por tanto que f (s) y g(s) son funciones suficientemente suaves y que f (s) > 0 ∀s ∈ I.

γ(s)

Figura 1: Superficie de revolución.

 Problema 1 (1 punto).
Calcule la base natural (formada por los vectores t1 y t2 ) del espacio tangente a esta superficie genérica
correspondiente al sistema de coordenadas generalizadas dado por s y θ.
Calcule la primera forma fundamental de esta superficie y explique de manera justificada si es ortogonal
el sistema de coordenadas propuesto.
Problema 2 (1 punto).
Calcule la curvatura y torsión de las lı́neas coordenadas.

Problema 3 (1 punto).
Indique el número de sı́mbolos de Christoffel Γkij (sin tener en cuenta relaciones de simetrı́a como Γkij =
Γkji ) que empleamos para calcular derivadas sobre esta variedad.
Calcule los sı́mbolos de Christoffel Γkij de esta variedad en el sistema de coordenadas indicado.
Problema 4 (1 punto).
Dado un móvil que se desplaza sobre esta superficie con velocidad v = t1 + t2 ¿es constante la velocidad
de este móvil? razone su respuesta. Calcule su aceleración en caso contrario.

Problema 5 (1 punto).
Calcule el vector normal unitario y la segunda forma fundamental de esta superficie.
Problema 6 (1 punto).
Calcule las curvaturas gaussiana y media de esta superficie. Explique de manera razonada las condi-
ciones que deben cumplirse para que un punto de esta superficie sea parabólico, elı́ptico o hiperbólico
(respectivamente).

Problema 7 (1 punto).
Explique de manera razonada las condiciones que deben cumplirse para que sean geodésicas las lı́neas
coordenadas s = cte y θ = cte.

Problema 8 (1.5 puntos).

Escriba la ecuación de continuidad de un fluido incompresible (∇ · v = 0) para un fluido con velocidad
v cuyo movimiento está confinado a esta superficie (dado por v = v 1 t1 + v 2 t2 ).

Problema 9 (1.5 puntos).

Dado el campo de velocidades v del ejercicio anterior determine la forma que toma en el presente sistema
de coordenadas la condición de flujo irrotacional, dada por

∇ × v = 0.

R
Para ello deberá extender la definición del anterior campo de velocidades al espacio 3 , lo cual puede
hacerse por medio v ≡ v 1 t1 + v 2 t2 + v 3 N , donde N es el vector unitario normal a la superficie y donde
hemos definido la tercera componente de v como v 3 ≡ 0.

Solución:

Problema 1
Base natural:
{ }
t1 = f˙(s) cos θ, f˙(s) sen θ, ġ(s) , t2 = {−f (s) sen θ, f (s) cos θ, 0} ,

Primera forma fundamental:

E = t1 · t1 = 1, F = t1 · t2 = 0, G = t2 · t2 = f 2 (s)

El sistema de coordenadas propuesto es ortogonal ya que F = t1 · t2 = 0 ∀s y ∀θ.

Métodos Matemáticos IV CURSO 2014-2015 Prueba de Evaluación Continua 2/6

 Problema 2
Curvatura de las lı́neas coordenadas θ = cte: { }
El vector tangente unitario es t1 . El vector de curvatura es k = ṫ1 = f¨(s) cos θ, f¨(s) sen θ, g̈(s) y la
( )1/2
curvatura es k = f¨2 (s) + g̈ 2 (s) .
Curvatura de las lı́neas coordenadas s = cte:
En este caso debemos tener en cuenta que las curvas s = cte no están descritas en términos del parámetro
longitud de arco. El vector tangente unitario es
t2
= {− sen θ, cos θ, 0}
∥t2 ∥

el vector de curvatura es
1 d t2 {cos θ, sen θ, 0}
k= =−
∥t2 ∥ dθ ∥t2 ∥ f (s)
y la curvatura es k = 1/f (s) tal y como cabı́a esperar, ya que por definición el radio de curvatura de
estas lı́neas es f (s).
Ambos conjuntos de lı́neas coordenadas son curvas planas, por tanto la torsión es en ambos casos nula.
Problema 3
Al ser una variedad bidimensional existen 8 sı́mbolos de Christoffel, dados por:

f˙(s)
Γ122 = −f (s)f˙(s), Γ212 = Γ221 =
f (s)

todos los restantes son nulos.

Problema 4
Tenemos un móvil con velocidad v = t1 + t2 , esto implica que para este móvil ds/dt = dθ/dt = 1 y
d2 s/dt2 = d2 θ/dt2 = 0 (siendo t el tiempo). Teniendo esto en cuenta la aceleración de este móvil es
( ) ( )
a = Γ111 + Γ112 + Γ121 + Γ122 t1 + Γ211 + Γ212 + Γ221 + Γ222 t2

A partir de la tabla de valores de Γkij calculada en el ejercicio anterior encontramos

f˙(s)
a = −f (s)f˙(s)t1 + 2 t2
f (s)

Por tanto la aceleración de este móvil sólo es nula si f˙(s) = 0, es decir, si la función f (s) es una constante,
en caso contrario la velocidad del móvil no es constante.
Es importante darse cuenta de que tanto el vector velocidad como el vector aceleración considerados
aquı́ son los medidos desde la variedad de manera intrı́nseca, es decir, sin considerar nada externo a
la propia variedad. Desde un punto de vista extrı́nseco, sin embargo, si consideramos que esta variedad
R
está inmersa en un espacio vectorial de dimensión superior ( 3 en este caso), tendremos que el vector
aceleración tiene una tercera componente perpendicular a la variedad considerada, dada por la aceleración
centrı́peta, cuyo valor es distinto de cero.
Problema 5
El vector normal unitario es { }
N = −ġ(s) cos θ, −ġ(s) sen θ, f˙(s)

y la segunda forma fundamental

L = f˙(s)g̈(s) − ġ(s)f¨(s), M = 0, N = f (s)ġ(s)

Métodos Matemáticos IV CURSO 2014-2015 Prueba de Evaluación Continua 3/6

 Problema 6
Calculamos en primer lugar los determinantes de la primera
( forma fundamental
) EG − F 2 = f 2 (s) y de
la segunda forma fundamental LN − M = f (s)ġ(s) f˙(s)g̈(s) − ġ(s)f¨(s) . Recordando la hipótesis de
2

partida f˙2 (s) + ġ 2 (s) = 1 deducimos

f˙(s)f¨(s) + ġ(s)g̈(s) = 0

lo que permite simplificar el determinante de la segunda forma fundamental a

LN − M 2 = −f (s)f¨(s)

una vez hecho esto calculamos:

Curvatura media ( )
EN + GL − 2F M 1 ġ(s) ˙ ¨
H= = + f (s)g̈(s) − ġ(s)f (s)
2 (EG − F )
2 2 f (s)
Curvatura de Gauss
LN − M 2 f¨(s)
K= = −
EG − F 2 f (s)
La discusión de casos (parabólico, elı́ptico, hiperbólico) depende del signo de la curvatura de Gauss, igual
al signo del determinante de la segunda forma fundamental dado que la primera forma fundamental tiene
siempre determinante positivo.
De los resultados anteriores se deduce inmediatamente que todos los puntos donde f¨(s) = 0 son puntos
parabólicos (al menos una de las curvaturas principales se anula y la superficie es localmente [en un
entorno pequeño alrededor del valor de s para el que f¨ se anula] equivalente a un “tronco de cono”). Los
puntos donde f¨(s) > 0 son hiperbólicos (las curvaturas principales tienen signos opuestos) y los puntos
donde f¨(s) < 0 son elı́pticos (ambas curvaturas principales tienen el mismo signo).
Es interesante observar que estos resultados son independientes de la coordenada θ, tal y como cabı́a
esperar al tratarse de una superficie de revolución.
Problema 7
A partir de los resultados anteriores es inmediato calcular que la curvatura geodésica de las curvas θ =
cte es nula, por tanto las lı́neas coordenadas de s (los “meridianos” de la superficie) siempre son curvas
geodésicas (en todos los puntos de la superficie, independientemente de las funciones f (s) y g(s)).
En cuanto a las lı́neas coordenadas correspondientes a s = cte (los “paralelos” de la superficie) su
curvatura geodésica es −f˙(s)/f (s), de modo que estas curvas solo son geodésicas en los máximos y
mı́nimos locales del radio (distancia al eje z) de la superficie, es decir, para aquellos valores de s para los
que f˙(s) = 0, independientemente del valor de θ.
Problema 8
La divergencia de un campo de velocidades v definido sobre esta variedad es
1 ∂ (√ i )
∇ · v = v;ii = √ gv
g ∂y i
que en este caso se simplifica a
[ ]
1 ∂ ( ) ∂v 2 1 ∂ ( ) ∂v̂ 2
∇·v = f (s)v 1 + = f (s)v̂ 1 +
f (s) ∂s ∂θ f (s) ∂s ∂θ

donde la última igualdad está escrita en componentes fı́sicas (dadas por: v̂ 1 = v 1 , v̂ 2 = f (s)v 2 , v̂ 3 = v 3 ).
La ecuación de continuidad ∇ · v = 0 queda finalmente como
∂ ( ) ∂v̂ 2
f (s)v̂ 1 + =0
∂s ∂θ
Por ejemplo, en el caso de un flujo axisimétrico (∂θ = 0) la condición de incompresibilidad se reduce a
f (s)v̂ 1 = cte.

Métodos Matemáticos IV CURSO 2014-2015 Prueba de Evaluación Continua 4/6

 Problema 9
Para calcular el rotacional del campo de velocidades tenemos que considerar la variedad desde el punto
R
de vista extrı́nseco, como un sub-conjunto de 3 . Sin embargo no es necesario extender las definiciones
R
del campo vectorial v y el sistema de coordenadas a todo 3 , sino que nos basta con hacer esto en
un entorno infinitesimal alrededor de la superficie, ya que solo vamos a evaluar ∇ × v sobre la propia
superficie.
Teniendo esto en cuenta definimos el tercer vector de la base como N (tal y como indica el enunciado) y la
tercera coordenada generalizada, y 3 , como la distancia recorrida en la dirección normal a la superficie (es
decir, en la dirección de N ). Como indicábamos antes solo es necesario definir esta tercera componente en
un entorno diferencial alrededor de la superficie, es decir, los valores posibles de y 3 se reducen al intervalo
y 3 ∈ [−dy, +dy]. Esto
{ } supone una simplificación importante, ya que extender el sistema de coordenadas
R
y 1 = s, y 2 = θ, y 3 a todo 3 podrı́a ser complicado dependiendo de la forma de la superficie, la cual
de momento no ha sido especificada (solo sabemos que tiene simetrı́a de revolución).
Parametrizamos entonces el espacio en un entorno infinitesimal alrededor de la superficie de acuerdo a

r = f (s) cos θ i + f (s) sen θ j + g(s) k + y 3 N

siendo N el vector normal unitario a la superficie calculado previamente. Con esto las coordenadas
y 1 = s e y 2 = θ especifican el punto de la superficie donde nos encontramos y la tercera coordenada
(y 3 ∈ [−dy, +dy]) representa la “altura” sobre dicho punto de la superficie.
De acuerdo a esta parametrización la base natural está dada por los vectores t1 , t2 y N . Es evidente
que esta base es ortogonal en todos los puntos, siendo los factores de escala
[( )2 ( )2 ]1/2
h1 = f˙(s) − y 3 g̈(s) + ġ(s) + y 3 f¨(s) , h2 = f (s) − y 3 ġ(s), h3 = 1

que evaluados sobre la superficie (y 3 = 0) se reducen a

h1 = 1, h2 = f (s), h3 = 1

Para calcular el rotacional de v aplicamos

ϵijk 1
∇ × v = εijk vk;j ti = √ vk;j ti = [(v3;2 − v2;3 ) t1 − (v3;1 − v1;3 ) t2 + (v2;1 − v1;2 ) t3 ]
g f (s)
√
donde hemos evaluado g en y 3 = 0.
En esta ecuación aparecen derivadas de la velocidad respecto de y 3 . En principio v solo está definido
sobre la superficie (es decir, en y 3 = 0), de modo que para calcular estas derivadas debemos extender
la definición de v a un entorno infinitesimal alrededor de la superficie. La manera menos arbitraria de
realizar esta extensión es por medio del transporte paralelo del vector v en la dirección de la normal, es
decir imponemos que vi;3 = 0 con i = 1, 2, 3:

∂v1
v1;3 = − vk Γk13 = 0 con y 3 ∈ [−dy, +dy]
∂y 3
∂v2
v2;3 = − vk Γk23 = 0 con y 3 ∈ [−dy, +dy]
∂y 3
∂v3
v3;3 = − vk Γk33 = 0 con y 3 ∈ [−dy, +dy]
∂y 3
Este conjunto de ecuaciones (junto con las condiciones iniciales dadas por los valores conocidos de vi
sobre la superficie, incluyendo v3 = 0 en y 3 = 0) permite calcular la extensión del campo v a un entorno
infinitesimal alrededor de la superficie. En nuestro caso este sistema nos permite despejar las derivadas
respecto a y 3 que aparecen en la expresión del rotacional de v.
Sustituyendo las condiciones anteriores y recordando que v3 es nulo en la superficie encontramos
[ ( ) ]
1 ∂v2 ∂v1
∇×v = −vk Γk32 t1 + vk Γk31 t2 + − N
f (s) ∂s ∂θ

Métodos Matemáticos IV CURSO 2014-2015 Prueba de Evaluación Continua 5/6

 Para completar el cálculo necesitamos los sı́mbolos de Christoffel Γk3j con j, k = 1, 2, evaluados sobre la
superficie. Derivando los factores de escala hi y sustituyendo posteriormente y 3 = 0 encontramos

Γ113 = Γ131 = −f˙(s)g̈(s) + ġ(s)f¨(s), Γ122 = −f (s)f˙(s)

f˙(s) ġ(s)
Γ212 = Γ221 = , Γ223 = Γ232 = −
f (s) f (s)
Γ3 = f˙(s)g̈(s) − ġ(s)f¨(s),
11 Γ322 = f (s)ġ(s)

siendo nulos los restantes.

A partir de este resultado llegamos finalmente a
[ ( ) ( ) ]
1 ġ(s) ∂v2 ∂v1
∇×v = v2 t1 − v1 f˙(s)g̈(s) − ġ(s)f¨(s) t2 + − N
f (s) f (s) ∂s ∂θ

que en componentes fı́sicas puede escribirse como

( ) ( )
ġ(s) 1 ∂f (s)v̂ 2 ∂v̂ 1
∇ × v = v̂ 2 t̂1 − v̂ 1 f˙(s)g̈(s) − ġ(s)f¨(s) t̂2 + − N
f (s) f (s) ∂s ∂θ

donde se ha tenido en cuenta que las componentes fı́sicas covariantes de v coinciden con las contrava-
riantes (v̂i = v̂ i ), dado que el sistema de coordenadas empleado es ortogonal en todos los puntos de la
superficie.
Como puede verse el rotacional de v contiene una parte en la dirección de la normal y otra en el plano
tangente. La parte de ∇ × v contenida en el plano tangente describe giros impuestos en el campo de
velocidades por la curvatura de la propia superficie, de modo que desde el punto de vista de un observador
que mida el campo de velocidades v desde la superficie solo es relevante la componente de ∇ × v según
la normal.
Para un observador situado sobre la superficie, por tanto, la condición de flujo irrotacional se reduce a
(∇ × v) · N = 0, que en nuestro caso queda como

∂f (s)v̂ 2 ∂v̂ 1
− =0
∂s ∂θ

Por último es interesante observar que las componentes de ∇ × v según t1 y t2 dependen de la forma
especı́fica que decidamos emplear para extender la definición de v a un entorno infinitesimal alrededor
de la superficie. En el cálculo anterior hemos decidido hacer esto por medio de vi;3 = 0, que parece la
elección menos arbitraria. Sin embargo, si realizamos esta extensión por medio de vi,3 = 0 (algo más
arbitrario pero considerablemente más sencillo) llegamos directamente a
( )
1 ∂f (s)v̂ 2 ∂v̂ 1
∇×v = − N
f (s) ∂s ∂θ

A efectos de calificar esta PEC damos por válida cualquiera de las 2 posibilidades anteriores.

Métodos Matemáticos IV CURSO 2014-2015 Prueba de Evaluación Continua 6/6