Se deducen expresiones para las componentes viscosa e inercial de la potencia en flujo laminar impermanente en tubos rectos, circulares, para fluidos newtonianos incompresibles. Se evalúan dichas componentes para flujo unidireccionalmente acelerado y flujo oscilatorio, como ilustración del método. Se discute la importancia relativa de las componentes de la potencia para los ejemplos analizados

Durán B., P.

Letelier S., M.

Russo A., A.

Contribuciones Científicas y Tecnológicas

Contribuciones Científicas y Tecnológicas Nº 64 - Area Mecánica VI - Enero 1984 potencia en flujo laminar impermanente en tuberías* M. Letelier S.••, A. Russo A.•• y P. Durán B.•• RESUMEN: Se deducen expresiones para las componentes vis-cosa e inercial de la potencia en flujo laminar impermanente en tubos rectos, circulares, para fluidos newtonianos incompresi-bles. Se evalúan dichas componentes para flujo unidireccional-mente acelerado y flujo oscilatorio, como ilustración del méto-do. Se discute la importancia relativa de las componentes de la potencia para los ejemplos analizados. SUMMARY: Analytical expressions are presented for the vis-cous and inertial components of pipe flow power for laminar unsteady flow in straight, circular tubes when the fluid is Newtonian incompressible. These components are evaluated for unidirectionally accelerated flow and for oscillatory flow . The relative importance of the power components are discussed for the examples herein analysed. Manuscrito revisado y aprobado en forma definitiva en noviembre de 1983. * ** Departamento de Ingeniería Mecánica, Facultad de Ingeniería, Universidad de Santiago de Chile (USACH). 5 l. INTRODUCCION El conocimiento de la potencia requerida por un sistema fluido-tubería es importante para el diseño de duetos y sistemas de impulsión, y también para el entendimiento de escurrimientos complejos como los que, por ejemplo, se encuentran en forma natural en arterias y otros conductos fisiológicos. Para los fines de este trabajo, se definirá la potencia específica consumida por una cañería como la diferencia de energía en el fluido por unidad de tiempo y peso entre los extremos de la cañería, para una longitud de ésta dada. La potencia así definida es, en general, sólo una parte de la potencia que debe aportar el sistema de impulsión . Si el flujo es isotérmico, incompresible, y la sección de la tubería es constante, la variación de energía específica H entre los extremos de entrada ( l) y salida (2) de la cañería está dada por la variación de la suma de energías específicas de presión y potencial, a saber H = _E!_ + h 1 - ( ~ + h2 ) , (1) 'Y 'Y donde p es la presión, hes la altura y 'Y es el peso específico del fluido. La potencia consumida W, para un caudal volumétrico Q, es W=-yQH (2) Cuando el flujo es temporalmente permanente, la energía y potencia gastadas son utilizadas exclusivamente en vencer el roce viscoso en el interior de la cañería. Un balance de cantidad de movimiento indica que, en este caso, para una cañería recta, circular de radio a y longitud L, existe la relación 2L H = -- T0 , (3) a-y donde T0 es el esfuerzo cortante en la pared del dueto. La relación (3) permite expresar alternativamente la potencia como W = 2L Qrº . a (4) La energía específica H puede ser también descrita por la ecuación de Darcy-Weisbach, a saber L _U2 H =f -2a 2g (5) en la cual fes el coeficiente de fricción, U es la velocidad media (U = Q/7r a2) y g es la aceleración de gravedad. De (5) y (2) resulta 1rapLfU2 w =----4g (6) expresión que es válida tanto para flujo laminar como turbulento y en la cual pes la densidad. Para flujo laminar, con f = 64v/2aU (v = viscosidad cinemática) , se tiene W = 81rLµU 2, siendo µ la viscosidad dinámica. 6 (7) Si el flujo es impermanente, es decir, si la velocidad varía temporalmente, expresiones universales similares a (4) y (7), son más difíciles de obtener aun en el caso más simple de flujo , esto es, laminar y paralelo. Es el objetivo de este trabajo presentar algunas expresiones originales válidas para flujo laminar impermanente que expresan la potencia consumida en términos de diversas variables del flujo. Los resultados que se plantean en lo que sigue se han derivado del modelo general de análisis para flujo laminar impermanente en duetos circulares publicado en referencia 1• 2. EXPRESIONES PARA LA POTENCIA CONSUMIDA EN FLUJO IMPERMANENTE LAMINAR EN TUBOS CIRCULARES El movimiento impermanente de un fluido incompresible newtoniano, en reg1men laminar, a través de un tubo recto y largo, de sección constante, es gobernado por la siguiente forma de las ecuaciones de Navier-Stokes, en coordenadas cilíndricas, V ( ..-ª=!!_ + _1 ..2.!!_) , ar2 r ar (8) en la cual r, x, t son las coordenadas radial, axial y temporal, respectivamente, u es la velocidad axial y p = p + -y h = presión piezométrica. Dado que el flujo es paralelo, la presión piezométrica no varía espacialmente en la sección transversal del dueto y sólo es función del tiempo. La potencia total es ahora también función del tiempo y obedece a la relación (9) combinando (8) y (9) se obtiene a a W = 21r p L t u :~ rdr - 21r µ L t (10) La ecuación (10) indica que en flujo impermanente la potencia tiene una componente inercial y otra viscosa, en contraste al movimiento estacionario, donde sólo existe la segunda2; dichas componentes son, respectivamente, W 1 = 21rp L tª u :~ r dr, ( 11) I · a au (12) Wv = - 21rµ L O u élr (r élr ) dr. Se definirán, por conveniencia, variables adimensionales en lo que sigue. Sean U0 y T 0 valores· constantes de referencia para la velocidad y el tiempo; con esto se tiene t* = tff0 ; r* = r/a; u* = U/U0 W* = W/1r µ L U~ ; <f>* = a2 (P1 - P2)/µ U0 L; O= a2/vT0 • (13) En lo anterior, la misma adimensionalización vale para W1 y W.,, <f>* es la llamada función forzante y O es el número de impermanencia. De acuerdo a1, y respetando la condición de no-deslizamiento en la pared del dueto, la velocidad se puede expresar como 7 u = dA dt (l - r6) + ... , (14) A es una función del tiempo relacionada con <I> a través de la ecuación n dA n2 d2A ~ A + - -- + + ... = (15) 22 dt 2242 dt2 4 A partir de la ecuación ( 14) sólo se utilizarán variables adimensionales, razón por la cual se omitirán los asteriscos en lo que resta de este trabajo. Las formas adimensionales del esfuerzo cortante en la pared y de la velocidad media son, respectivamente, _ au ( l) ar ' U= 2 Lª urdr, las cuales, cuando deducidas de ( 14), se convierten en O dA T =2A+- -+ 0 4 dt U = ~ + __Q_ dA + .... 2 24 dt (16) (17) (18) (19) Series del tipo ( 15), ( 18) y ( 19) pueden ser invertidas para expresar A en función de <I>, T0 y U. Si estas nuevas series se introducen en las formas adimensionales de (11) y ( 12), se obtienen series para W1 y Wv en términos de <I> , T y U. Las fórmulas que se plantean a continuación son el resultado de esta operación, donde también se ha incluido como variable dependiente la velocidad en el centro del tubo Uc , la cual se obtiene de ( 14) poniendo r =O.Así se encuentra 4 dU n2 d2U ( dU )2 J W1=- nu --- [U + + 3 dt 144 dt2 dt + !}3 [U d3U + 2 dU ~ ] + ... , 4320 dt3 dt dt2 (20) Wv = 8 U2 + n2 ( dU )2 _ ___!!:_ dU ~ + --144 dt 2160 dt dt2 (21) w = 8 u 2 + 4 nu dU n2 u d2U + . "' - ---3 dt 144 dt2 (22) (23) 8 w = <!>2 - o <l> d<j> + 1102 [( d<j> )2 V 8 24 dt 3072 dt + 2"' d2<j> l + 'I' dt2 . . . ' W, = O T dTo 12 ° dt 02 1[( dTo )2 + To dd2t~o] + ... , 256 dt W = T~ + 0 T dTo 2 24 ° dt 02 [2 (dTo )2 + To d2T0 ] + 768 dt ~ ... , O duc + 0 2 W1 = -u 3 e dt 192 + Uc d2uc ] + dt2 .. • , w = 2 u~+ o duc 02 d2uc _ 3 ( duc )2] + -u -- [2u V 12 e dt 1152 e dt2 dt ... ' w 2u~ + 50 duc + 02 [ 9 ( duc )2 + 4 Uc d2uc J + ... = 12 u., dt 1152 dt dt2 Las fórmulas (25), (28) y (31) se han obtenido haciendo W = W 1 + Wv. (24) (25) (26) (27) (28) (29) (30) (31) La integración espacial de la forma adimensional de (8) conduce a la ecuación unidimen-sional de movimiento, a saber dU 0-+2T =A,. dt o 'I' Si (32) se multiplica por U resulta W = <j>U = ou ~~ + 2 To U (32) (33) donde <j>U es la potencia total exacta pero OU dU/dt y 2 T0 U representan solamente aproxima-ciones para las potencias inicial y viscosa. En particular, si O -o, 2 T0 U es exactamente la forma adimensional de (4). Combinando (18) y (19) resulta T0 = 4 U + O dU 6 dt 0 2 d2U 288 dt2 + ... (34) expresión que, al ser sustituida en (33) , genera la fórmula para la potencia viscosa aproximada según (33). Así se tiene 9 w, = nu dU dt w~ = 2 ToU = 8 U2 + n u dU 3 dt (35) 02 U d2U + 144 dt2 (36) Las ecuaciones (35) y (36) representan, respectivamente, las potencias inercial y viscosa asociadas directamente a la velocidad media, en tanto (20) y (21) entregan estas magnitudes referidas a la velocidad local u. 3. EV ALUACION DE LA POTENCIA CONSUMIDA EN FLUJO UNIDIRECCIONALMENTE ACELERADO Y EN FLUJO OSCILATORIO Se considerarán dos casos especiales de flujo impermanente con el objeto de ilustrar algunas de las fórmulas deducidas en el párrafo precedente. Flujo unidireccionalmente acelerado: Se definirá por la expresión U = 1 + t (37) Las formas apropiadas para determinar la potencia están dadas por (20), (21) , (22), (35) y (38). Dado que en este ejemplo U tiene un número finito de derivadas no nulas , las expresio-nes para la potencia son válidas para todo valor de n. En la figura 1 se describe la evolución temporal de las potencias inercial (20), viscosa (2 1) y total (22) exactas para O= 3 y O= 20. W1 w Wt w Wv Wv .n.: 3.0 \ .O.: 20 180 \ u : 1 • t 4 Sx10' 180 \ u : 1 • t 1,8 \ ~\ W=W'-/ \ 160 Wv ""2.\ 1,6 160 \ '!!J. 'º \ \f Wv \ 140 35 140 ',,¡ 1,·, \ 120 \ 30 12 O 1,2 \ \ ¡ · ', 100 ' 2S 100 1,0 ' ' 80 ' 20 80 0,8 ' / ' &O 15 60 0,6 'º 10 40 ~ º·' 20 W¡ 5 20 0,2 ~---- ---=------ Wi o o o o o 1,0 2,0 3,0 t o 1,0 2,0 3,0 Figura 1. Componentes exactas viscosa (W v) e inercial (W 1), y aproximadas viscosa cw:) e inercial (W;) de la potencia W en flujo unidireccionalmente acelerado. 10 Se aprecia que para flujo lentamente acelerado (0 = 3), la potencia viscosa predomina, en tanto la contribución de la potencia inercial es pequeña. Para flujo fuertemente acelerado (0 = 20) la mayor parte de la potencia se gasta inicialmente en vencer la inercia del fluido. Las formas aproximadas de las potencias inercial y viscosa entregan valores cercanos a los exactos para O = 3, pero existe una discrepancia apreciable entre valores aproximados y exactos para O = 20. Debe hacerse resaltar que en este modelo simple de flujo acelerado, la potencia, en todas sus formas, aumenta con el tiempo, lo cual no es necesariamente válido para otros sistemas mecánicos, tales como aquellos compuestos por masas sólidas. En estos últimos, a veces el torque inicial requerido para poner en movimiento el sistema en un tiempo dado, exige que el mayor consumo de potencia se origine al comienzo del movimiento, siendo necesario posteriormente solamente suministrar potencia para vencer las fuerzas de roce, cuando éste es el caso, lo cual exige una magnitud de potencia menor que la inicial. En relación al movimiento de un fluido, el roce viscoso aumenta con la velocidad, por lo cual, al menos la componente viscosa de la potencia, siempre crece con la velocidad . Flujo oscilatorio: Se considerará una variación temporal de la velocidad media dada por U = cost (38) La expresión para U, en este ejemplo, admite infinitas derivadas diferentes de cero, lo que hace a su vez infinitas las series para la potencia. Dichas series, en particular las ecuaciones (20), (21), (22), (35) y (36), fueron evaluadas para O = 2, valor del número de impermanencia para el cual existe buena convergencia. En la figura 2 se han graficado las diversas componentes de la potencia, a partir de la ecuación (38). En el caso cuando O = 2, la contribución de la potencia viscosa es fundamental existiendo desfase entre ésta, la potencia inercial y la potencia total. Se ilustra además la discrepancia entre las formas aproximadas y exactas de las componentes de la potencia. w .o.= 2 u=cost .n. = 2 u= e os t 8,0 w 7,0 j ·~ wj Wv 6,0 5,0 4,0 3,0 2,0 1.0 o \ / \ '-·· '/ ,,.. 1 \ \ / , ••• W¡ •'/ -1,0 ' / ' " '~-- ___ .. :., ,::.:/ --,,. --__ ,,, __ ,. W¡ W¡ O 1.0 2,0 3,0 4,0 1,0 2,0 3,0 4, 0 Figura 2. Componentes exactas viscosa (W .) e inercial (W 1), y aproximadas viscosa (W:) e inercial (W;) de la potencia W en flujo oscilatorio. 11 La evaluación de la potencia para valores más altos den se hace bastante más compleja y no es fácil obtener expresiones explícitas como las aquí presentadas. No obstante, cálculos numéricos indican que para valores mayores de O la importancia relativa de la potencia inercial aumenta notablemente. En arterias de seres vivos3 con diámetro superior a 4 mm, por ejemplo, se tienen > 10. En estas condiciones, deberá esperarse que la potencia gastada por el organismo en generar el flujo sanguíneo, en gran medida esté destinada a vencer la inercia de la sangre, al menos en partes del ciclo cardíaco. En arterias de pequeño diámetro, en cambio, la fórmula cuasiperma-nente (7) de la potencia será una buena aproximación. Si bien el ciclo cardíaco se refleja en una dinámica circulatoria más compleja que la descrita por la ecuación (38), se estima que los resultados predichos por este modelo simple de flujo oscilatorio, son extensibles cualitativa-mente en cierta medida a aquélla. 4. CONCLUSIONES Se han obtenido expresiones explícitas para la evaluación de las componentes inercial y viscosa de la potencia en flujo impermanente laminar en tubos circulares rectos con las restricciones adicionales indicadas en el texto. Las fórmulas se han aplicado a flujo unidireccionalmente acelerado y oscilatorio, lo cual ha permitido poner de manifiesto la importancia relativa de las componentes de la potencia para diferentes valores del número de impermanencia. Se estima que estos resultados son de interés tanto para el diseño de líneas fluidas como para el conocimiento de la dinámica de flujos fisiológicos, en particular del sistema circula-torio. REFERENCIAS l. LETELIER S .. M.F. and H.J . LEUTHEUSSER, " Unified Approach to the Solution of Problems of Unsteady Laminar Flow in Long Pipes", Joumal of Applied Mechanics. Vol. 50. Nº 1 (1983). 2. MORAGA B., N. , F. ESTEFFAN M. y M. LETELIER S .. " Pérdida de energía en flujo laminar transiente". Contribuciones, Nº 20 (1976). 3. LETELIER S., M.F. and H .J. LEUTHEUSSER, "Analytical Deduction of the lnstantaneous Velocity Distribu-tion, Wall Shear Stress and Pressure Gradient from Transcutaneous Measurements of the Time-varying Rate of Blood Flow", Medica! & Biological Engineering & Computing. Vol. 19 (1981) . 12 

Potencia en flujo laminar impermanente en tuberías

http://www.revistas.usach.cl/ojs/index.php/contribuciones/article/download/4215/26003283/

Potencia en flujo laminar impermanente en tuberías

Abstract

Similar works

Full text

Available Versions

Contribuciones Científicas y Tecnológicas