Flujo Transitorio de Couette
Flujo Transitorio de Couette
Flujo Transitorio de Couette
Se tiene un fluido inicialmente en reposo, entre dos placas planas Las condiciones de frontera son:
paralelas de infinita extensión, separadas por una distancia δ . (1) en y 0 : v x 0
En t 0 , la placa superior se pone instantáneamente en
(2) en y δ : v x v0
movimiento con una velocidad v 0 , mientras que la placa inferior
se mantiene fija. Se desea determinar el perfil de velocidad del y se necesita también una condición inicial:
fluido en función de la posición y del tiempo, v x y , t .
(i) en t 0 : v x 0
v x v y v
z 0 donde se requiere que u y , t 0 cuando t . Entonces, se
x y z
sustituye v x s u en la ecuación diferencial:
donde se cancelan todos los términos según las suposiciones.
v x 2v x 2
ν s u ν 2 s u
La ecuación de conservación de momentum, para fluidos t y 2 t y
newtonianos de propiedades constantes, en coordenadas s u 2s 2u
ν 2 2
rectangulares, componente x es: t t y y
REVISIÓN 1
Página 1 de 4
86948.89
Solución para la parte estable s y Como Y no depende de t se puede sacar como constante de la
primera derivada. Lo mismo se puede decir de T en la segunda
La ecuación diferencial es:
derivada, porque no depende de y , con lo que se obtiene:
d 2s
0 T 2Y
dy 2 Y νT 2
t y
con condiciones de frontera:
Dividiendo entre νYT (la viscosidad cinemática ν se deja del lado
(1) en y 0 : s 0 más sencillo de la ecuación):
(2) en y δ : s v0
1 T 1 2Y
(nótese que la ecuación diferencial aparece ya en derivadas νT t Y y 2
ordinarias, puesto que s sólo depende de y , y que la viscosidad
cinemática ν se ha eliminado al pasarla dividiendo al cero).
El lado izquierdo de la ecuación tiene T y su derivada, que
dependen de t únicamente ( ν es constante), mientras que el
La ecuación diferencial se resuelve por simple integración, para
dar la solución general: lado derecho tiene Y y su segunda derivada, por lo que sólo
depende de y . El tiempo puede variar desde cero hasta infinito
s C1 y C2 ( t 0 ), la posición puede variar desde cero hasta la separación
entre las placas ( 0 y δ ), y de alguna manera los dos
Luego de aplicar las condiciones de frontera, se llega a la solución miembros de esta ecuación siempre deben ser iguales entre sí. La
particular: única manera en que puede mantenerse la igualdad, es que sean
iguales a una constante, llamada constante de separación.
y
s v0
δ
Estrictamente, deberían considerarse tres casos para dicha
constante de separación: cero, un valor positivo o un valor
Solución para la parte transitoria u y , t negativo, pero resulta que sólo una constante de separación
negativa da soluciones físicamente posibles, por lo que los otros
La ecuación diferencial es:
dos casos no se analizarán.
u 2u
ν 2 Tomando la constante de separación como λ2 , se iguala cada
t y
miembro de la ecuación con la constante de separación, y se
con condiciones de frontera (homogéneas): obtienen las dos ecuaciones diferenciales siguientes, cada una de
las cuales se debe resolver por separado:
(1) en y 0 : u 0
(2) en y δ : u 0 1 T dT
λ2 νλ2T 0
νT t dt
y se necesita también una condición inicial, obtenida a partir del 1 2Y d 2Y
λ2 λ2Y 0
planteamiento original v x s u , es decir, u v x s , donde s Y y 2 dy 2
es la solución estable obtenida anteriormente:
(nótese que ahora se emplean derivadas ordinarias puesto que
y y cada función depende únicamente de una variable).
en t 0 : vx 0 y s v0 u 0 v0
δ δ
Entonces, se llega a la condición inicial para u : Solución para la función del tiempo T t
y dT
(i) en t 0 : u v0 νλ2T 0
δ dt
Para resolver la ecuación diferencial, se aplica el método de Ésta es una ecuación diferencial ordinaria de primer orden, lineal,
separación de variables. Se asume que u y , t es el producto de homogénea, que se puede resolver por separación de variables:
una función de y y una función de t :
dT
νλ2T
u y, t Y y T t dt
dT
(en este método, se suele indicar estas funciones con las mismas νλ2dt
T
letras que las variables independientes, pero en mayúsculas). dT
REVISIÓN 1
Página 2 de 4
86948.89
Solución para la función de la posición Y y Aplicando la condición de frontera (1):
0 A cos 0 B sen 0 e νλ t
2
d 2Y
λ2Y 0
dy 2
cos 0 1 y sen 0 0 , con lo que se tiene:
Ésta es una ecuación diferencial ordinaria de segundo orden, 2
0 Ae νλ t
lineal, homogénea, de coeficientes constantes, que se puede
2
resolver por el método de la ecuación característica. Para ello, Al pasar dividiendo la función exponencial e νλ t se tiene A 0 ,
primero se escribe la ecuación empleando la notación D para con lo que la solución queda simplemente:
derivadas:
u y , t B sen λy e νλ t
2
2 2
D Y λY 0
Ahora se aplica la condición de frontera (2):
y se “factoriza” Y :
0 B sen λδ e νλ t
2
D 2 2
λ Y 0
2
Nuevamente, se pasa dividiendo e νλ t :
La ecuación se vuelve cero cuando cualquiera de los dos
“factores” es cero. Obviamente Y no puede ser cero (el perfil de 0 B sen λδ
temperatura no dependería de y ), por lo que debe ser D 2 λ2
lo que la haga cero. De aquí se obtiene la ecuación característica, pero ahora no se puede pasar también dividiendo sen λδ
remplazando el operador diferencial D por una variable m : porque entonces se tendría B 0 , y eso haría que u y , t 0
siempre, eliminando por completo la parte transitoria de la
m2 λ2 0 solución. Entonces, como B no puede ser cero, se llega a la
conclusión que sen λδ 0 . Sin embargo, más de un valor del
Resolviendo la ecuación característica:
argumento hace que la función seno sea cero:
m2 λ2 sen π 0
m λ2 sen 2π 0
m λi
sen 3π 0
Se obtienen raíces imaginarias ya que, por definición, λ2 es un
número negativo. Como son raíces imaginarias conjugadas, la
sen nπ 0
solución para Y está dada por funciones trigonométricas:
para cualquier valor de n entero. Igualando sen λδ sen nπ
Y C 4 cos λy C 5 sen λy
se deduce que λδ nπ , por lo que λ nπ / δ . Pero habrá un
número infinito de λ , una para cada valor de n , por lo que se
identifican con un subíndice:
Solución para la parte transitoria u y , t
(continuación) nπ
λn
Una vez que se obtuvieron las soluciones para T y Y , se puede δ
escribir la solución completa para u :
estos valores de λ n son los valores característicos de este
u y, t Y y T t sistema. Además, para cada valor de lambda se tiene una
solución u diferente, con su correspondiente constante B
u y , t C 4 cos λy C 5 sen λy C 3e νλ t
2
u y , t A cos λy B sen λy e νλ t
2
u y,t u y,t
n
Ésta es la solución general para u en la que se deben aplicar las n 1
condiciones de frontera e inicial, para encontrar las tres
u y,t B sen λ y e νλ n t
2
n 1
u y,t sen λ n y e νλ n t
2
δ λn
mπy n 1
y multiplicar toda la ecuación por sen :
δ
Solución final: perfil de velocidad v x y , t
B sen
y mπy nπy mπy
v0 sen n sen Recordando que v x y , t s y u y , t , lo único que falta es
δ δ n 1
δ δ
sumar las dos soluciones para obtener el perfil de velocidades:
e integrar con respecto a y de 0 a δ :
1 n
y 2v
vx y, t sen λ n y e νλ n t
2
δ v0 0
δ
y mπy nπy mπy
0
v 0 sen
δ δ
dy
0
n 1
Bn sen
δ
sen
δ
dy δ δ
n 1
λn
nπ
Aplicando la propiedad de linearidad de la integración respecto a donde los valores característicos son λ n .
δ
la sumatoria, y sacando v0 / δ y Bn como constantes:
Con esta solución analítica es posible calcular la velocidad del
v0 mπy
nπy mπy
δ δ
δ 0
y sen
δ
dy
n 1
Bn
0
sen
δ
sen
δ
dy
fluido en cualquier posición y y en cualquier tiempo t .
v0 δ
mπy δ
2 nπy
δ
0
y sen
δ
dy Bn
sen
0 δ
dy
Despejando Bn :
v0 δ
nπy
Bn
δ 0
y sen
δ
dy
δ
nπy
0
sen 2
δ
dy
2v 0 1
n
Bn
nπ
o bien, como nπ λ n δ :
2v 0 1
n
Bn
λ nδ
u y,t B sen λ e νλ
2
t
ny
n
n
n 1
REVISIÓN 1
Página 4 de 4
86948.89