Patron de lineas de corinteinstantineas para fujo alrededor de un clindo en rotaciéa con rnimero de Reynolds 1000, tomado de Nair et al. [19]. Esta configuracién, que oscila'debido al desprendimiento de los vértices, fue obtenida por Dindmica de Fluidos Computacional (CFD) y ‘coincide con los resultados de los experimentos de visualizacién de fujos (20). En este capitulo se obtendrin y se discutirin las ecuaciones de movimiento resueltas en este ejemplo por CFD. (Figura cortesta del profesor Tapan K. Sengupta.)Capitulo Capt) : ones diferenciales | para una particula fluida Motivacién, Cuando analizamos el movimiento de los fluidos podemos seguir dos caminos distintos: (1) buscar una estimacién de los efectos globales (flujo mésico, fuerza aplicada, intercambio de energia) sobre una regién finita o volumen de control, o (2) analizar unto a punto los detalles del campo fluido analizando una regién infinitesimal del flujo. El primer enfoque, de tipo global, se tratd en el Capitulo 3 Este capitulo trata de la segunda de las técnicas para analizar el movimiento de los fluidos:-el andlisis a pequefia escala, 0 diferencial. Esto es, aplicamos las cuatro leyes de conservacién bésicas a un volumen de control infinitesimal o, alternativamente, a un sistema fluido infinitesimal. En ambos casos se obtienen las ecuaciones diferenciales bbasicas del movimiento de un fluido. También se desarrollan las condiciones de contorno apropiadas. En su forma mas bésica, estas ecuaciones diferenciales del movimiento son bastante dificiles de resolver, y se conoce muy poco sobre sus propiedades mateméticas generales. Sin embargo, se pueden mostrar ciertos aspectos que tienen un gran valor educativo. En primer lugar, como se muestra en el Capitulo 5, las ecuaciones (aunque no se resuel- van) revelan los pardmetros adimensionales bdsicos que gobiernan el movimiento de los fluidos. En segundo lugar, como se muestra en el Capitulo 6, se puede encontrar un gran-niimero de soluciones titiles si se hacen dos hipstesis simplificadoras: (1) flujo estacionario y (2) flujo incompresible. Una tercera simplificacién bastante més dristica, la de flujo no viscoso, hace que sea vélida la ecuacién de Bemoulli y proporeiona una ‘gran variedad de soluciones ideales, o de fluido perfecto, posibles. Estos flujos idealizados se tratan en el Capitulo 8; se debe ser cuidadoso ¢ indagar si estas soluciones son de hhecho realistas cuando se comparan con el movimiento real del fluido, Finalmente, a pesar de su gran dificultad, las ecuaciones diferenciales generales se pueden resolver hhoy en dia mediante ia técnica aproximada del andlisis numérico, donde las derivadas se sustituyen por relaciones algebraicas entre un nimero finito de puntos del campo fluido que pueden resolverse:posteriormente- mediante un ordenador. La Referencia | es un ejemplo de un libro de texto dedicado integramente al andlisis numérico del mo- vimiento de los fluidos.lo 4. Relaciones diferenciales para una partcula fuida campo deacleraciones Ena Secon 15 exablecimas la forma vectorial cartesian de un campo de velocdades funcién de la posicién y del tiempo: Vor, 2) = in(x, y, 2, ) + jul, y, gt) + kee, y, 2 1) (14) 5 la Variable mids importante de-la Mecénica de Fluidos. Conocer el campo de velo- Cidades es a menudo equivalente a resolver el problema. Nuestras coordenadas estén fijas enel espacio y observamos cémo pasa el fluido: como si hubiéramos tallado un conjunto de Iineas de coordenadas sobre la ventana de cristal de un tinel de viento, Este es el mé- todo descriptivo euleriano, que es distinto al método lagrangiano, en el cual se sigue el movimiento de las particulas individuales. La aceleracién a también es fundamental en Mecénica de Fluidos, ya que aparece al aplicar la segunda ley de Newton a un sistema fluido infinitesimal. Por tanto, necesitamos calcular la derivada total del vector velocidad con respecto al tiempo: Como cada componente escalar(u, v, w) es una funcién de las cuatro variables (x,y, z, ‘), utiligamos la regla de la cadena para obtener la derivada temporal de cada escalar. Por ejemplo, duly, 8) _ 3 4 ud, dude dt a ay dt oc dt Pero, por definici6n, dx/dt es la componente u de la velocidad local, y dy/dt = v y deldt = ww. Ast pues, la derivada total de u se puede escribir en la siguiente forma compacta: du au uy uy au ae ae BE yen Bey. om Puy tis 9 "a a 0 ™ an Y sustituyendo u por vo w, se obtienen expresiones similares para duldt o dwidt. Suman- 4o estas expresiones para formar un vector, obtenemos la aceleracién total: (4.2) El término 2V/at se denomina aceleracién local y se anula cuando el flujo es estacionario, esto es, independiente del tiempo. Los tres términos entre paréntesis forman la acelera- cién convectiva, que aparece cuando la particula se mueve a través de regiones donde la ‘velocidad varia, como en una toberao en un difusor. En flujos nominalmente “estacionarios”, el fluido puede sufrir grandes aceleraciones a consecuencia de los téminos convectivos. Obsérvese el uso que hacemos del producto escalar entre V y el operador gradiente V: Fe ae dea cade ax a El concepto de la derivada total temporal —a veces llamada derivada sustancial 0 mate- riai— puede aplicarse a cualquier variable, como por ejemplo la presién: ap _ op vy oa 2 ° = 24k Ete +(V- Vp (43) Siempre que aparecen efectos convectivos en oe a bésicas de conservacién de la masa, Cantidad de movimiento 0 energia, las ecuaciones diferenciales basicas se vuelven4.2. La ecuacién diferencial de conservacién de la masa 227 10 lineales, lo que origina dificultades mateméticas que las hacen més complicadas que en los flujos que no suften cambios convectivos. Recalcamos que esta derivada temporal total sigue a una particula con una identidad fj, lo cual es conveniente para expresar las leyes de la mecénica en la descripcién eule- ele denominar derivada sustancial o material y a menudo se jo especial D/Dt como recordatorio de que tiene cuatro términos y sigue ‘una particula determinada ‘EJEMPLO 4.1 Dado el campo vectorial de velocidades euleriano V = 3a + xj + 7k determine [a aceleracién total de.una particula Solucién * Consideraciones: Conocemos las. componentes no estacionarias de la velocidad, ee + Procedimienio: Caleulamos todas las derivadas necesarias con respecto a (xy, 2 0 las susttuimos en el vector de aceleracibn total, Ecuacién (4.2), y agrupamos términos. +» Paso’ I> Primero obtenemos la aceleracién local W/O: e wv ar '* Paso-2: De forma similar, os términos de la aceleracién convectiva de: la-Ecuacién (4.2) on: a a Oya he qn. td a alt ko) = +i tek = e208 + xj + 9°) = GOI + J+ Ok) = 3h oa a 5G + x + 07k) = (az}(Oi + Oj + 2syke) = Dory k w a ate 2) wae Og Oa td + Fk) Ot a + Ob) = By ‘© Paso 3: Agrupando todos los términos, obtenemos la derivada “total” o.“sustancial”: v aes Few = Gi A) ed Dock + Guy 23+ Ga + ay t+ OF + Dk Rep. © Comentarios: Suponiendo que la expresiém dada para V es vélida-en todas partes, e! vector aceleracién total dV/dt es aplicable a todos los puntos instantes del campo ffuido: SS 4.2. La ecuacién diferencial Todas las ecuaciones diferenciales bésicas pueden deducirse considerando un volumen de conservacién de la masa de control elemental o un sistema elemental. Elegiremos aqui un volumen de controi infinitesimal fijo (dx, dy, dz) como el de la Figura 4.1, y utilizaremos las relaciones228 Capitulo 4, Relaciones diferenciales para una partcula fuida Figura 4.1. Volumen de control infinitesimal fo en coordenadas cartesianas mostrando los fljos mésicos de entrada y de salida en las caras perpendiculares alee x. bésicas para volmenes de control del Capftulo 3. El flujo a través de cada cara del clemento es aproximadamente unidimensional y la relacién de conservacién de la masa apropinda es aqui 4 aay + LOAM a ~ DOA ea= 0 6.22) Bl elemento es tan pequefo que la integral de volumen se reduce al témino diferencial: ] ap ap 88 yma 8 a, i. at OF poe Los términos del flujo mésico aparecen en las seis cara, tres de entrada y tres de salida. Hacemos uso del concepto de continuo del Capitulo 1, donde todas las propiedades flui- das se consideran descritas por funciones que varfan uniformemente con el tiempo y la posiciGn, tal como p = p (x,y, 2,4). Por tanto, si Tes la temperatura en la cara izquierda del elemento de la Figura 4.1, la cara derecha tendré una temperatura ligeramente dife- rente, T + (@7/@x) dr. Para la conservaciGn de la masa, si pu es dato en la cara izquierda, el valor de este producto en la cara derecha es pu + (Gpuldx) dx. [La Figura 4.1 muestra tinicamente los flujos en las caras izquierda y derecha. Los flu jos en las caras perpendicular a los ejes y (inferior y superior) y (anterior y posterior) se han omitido para mis claridad en el dibujo, Haremos un listado de estos seis flujos como sigue: Caras Flujo misieo de entrada Flujo méisicn de sada x pudy de [av Zana] dye y poaréc [ro Zona are = t pw dedy [ov + Lew nay Introduciendo estos términos en la Ecuacién (3.22), tenemos: i a # & e sy ded de + =) ded de += ip cd dy d+ + (ow) de dy de = 0