Teorema de la circulación de Kelvin

En la mecánica de fluidos, el teorema de circulación de Kelvin, llamado así por William Thomson, primer barón Kelvin que lo publicó en 1869, dice que: En un fluido barótropo ideal con fuerzas corporales conservadoras, la circulación alrededor de una curva cerrada (que encierra los mismos elementos del fluido) que se mueve con el fluido permanece constante con el tiempo. Dicho matemáticamente:

$\frac{\mathrm{D}\mathrm{\Gamma }}{\mathrm{D}t}=0$

donde $\mathrm{\Gamma }$ es la «circulación» alrededor de un contorno material $C(t)$. Dicho más simplemente este teorema dice que si se observa un contorno cerrado en un instante, y se sigue el contorno a lo largo del tiempo, siguiendo el movimiento de todos sus elementos fluidos, la circulación sobre los dos lugares de este contorno son iguales.

Este teorema no se sostiene en los casos de tensiones viscosas, fuerzas corporales no conservadoras, por ejemplo, una fuerza de coriolis, o relaciones de presión-densidad no barótropas.

Demostración matemática

Plantilla:Véase La circulación $\mathrm{\Gamma }$ alrededor de un conjunto material cerrado $C(t)$ viene definido por:

$\mathrm{\Gamma }(t)={\oint }_C𝒖\cdot \mathrm{d}𝒔$

donde u es el vector velocidad y ds es un elemento a lo largo del contorno cerrado

La ecuación gobernante para un fluido no viscoso con una fuerza corporal conservadora es:

$\frac{\mathrm{D}𝒖}{\mathrm{D}t}=-\frac{1}{\rho }\mathrm{\nabla }p+\mathrm{\nabla }\mathrm{\Phi }$

donde D/Dt es la derivada material, ρ es la densidad del fluido, p es la presión y Φ es el potencial de la fuerza del cuerpo. Estas son las ecuaciones de Euler con una fuerza corporal.

La condición de barotropicidad implica que la densidad es una función solo de la presión, es decir $\rho =\rho (p)$.

Tomando la derivada convectiva de la circulación se obtiene:

$\frac{\mathrm{D}\mathrm{\Gamma }}{\mathrm{D}t}={\oint }_C\frac{\mathrm{D}𝒖}{\mathrm{D}t}\cdot \mathrm{d}𝒔+{\oint }_C𝒖\cdot \frac{\mathrm{D}\mathrm{d}𝒔}{\mathrm{D}t}.$

Para el primer término, sustituimos de la ecuación gobernante, y luego aplicamos el teorema de Stokes, de la siguiente manera:

${\oint }_C\frac{\mathrm{D}𝒖}{\mathrm{D}t}\cdot \mathrm{d}𝒔=\underset{A}{\int }\mathrm{\nabla }\times (-\frac{1}{\rho }\mathrm{\nabla }p+\mathrm{\nabla }\mathrm{\Phi })\cdot 𝒏\mathrm{d}S=\underset{A}{\int }\frac{1}{{\rho }^2}(\mathrm{\nabla }\rho \times \mathrm{\nabla }p)\cdot 𝒏\mathrm{d}S=0.$

La igualdad final surge desde $\mathrm{\nabla }\rho \times \mathrm{\nabla }p=0$ debido a la barotropicidad. También hemos aprovechado el hecho de que el rotacional de cualquier gradiente es necesariamente 0 , o $\mathrm{\nabla }\times \mathrm{\nabla }f=0$ para cualquier función $f$.

Para el segundo término, observamos que la evolución del elemento de línea material viene dada por

$\frac{\mathrm{D}\mathrm{d}𝒔}{\mathrm{D}t}=(\mathrm{d}𝒔\cdot \mathrm{\nabla })𝒖.$

Por lo tanto

${\oint }_C𝒖\cdot \frac{\mathrm{D}\mathrm{d}𝒔}{\mathrm{D}t}={\oint }_C𝒖\cdot (\mathrm{d}𝒔\cdot \mathrm{\nabla })𝒖=\frac{1}{2}{\oint }_C\mathrm{\nabla }(|𝒖{|}^2)\cdot \mathrm{d}𝒔=0.$

La última igualdad se obtiene aplicando el teorema de Stokes. La última igualdad se obtiene aplicando el teorema de Stokes.

Como ambos términos son cero, obtenemos el resultado siguiente:

$\frac{\mathrm{D}\mathrm{\Gamma }}{\mathrm{D}t}=0.$

Teorema de la circulación de Poincaré-Bjerknes

Un principio similar que conserva una cantidad puede obtenerse también para el marco giratorio, conocido como el teorema de Poincaré-Bjerknes, llamado así por sus autores Henri Poincaré y Vilhelm Bjerknes, que derivaron la invariante en 1893. and 1898.

El teorema puede aplicarse a un marco giratorio que gira a una velocidad angular constante dada por el vector $\mathrm{\Omega }$, para la circulación modificada

$\mathrm{\Gamma }(t)={\oint }_C(𝒖+\mathrm{\Omega }\times 𝒓)\cdot \mathrm{d}𝒔$

Aquí, $𝒓$ es la posición del área del fluido. Según el teorema de Stokes, esto es:

$\mathrm{\Gamma }(t)=\underset{A}{\int }\mathrm{\nabla }\times (𝒖+\mathrm{\Omega }\times 𝒓)\cdot 𝒏\mathrm{d}S=\underset{A}{\int }(\mathrm{\nabla }\times 𝒖+2\mathrm{\Omega })\cdot 𝒏\mathrm{d}S$

Véase también

• Principio de Bernoulli
• Ecuaciones de Euler (fluidos)
• Teoremas de Helmholtz