4 - Balance Integrales Ultimo
4 - Balance Integrales Ultimo
4 - Balance Integrales Ultimo
Estos enfoques tienen que ver con el observador respecto al fenómeno cuyo movimiento se
estudia.
DESCRIPCION LAGRANGIANA
En la descripción Lagrangiana estudia cómo evolucionan en el tiempo todas y cada una de
las principales variables asociadas a las partículas de un fluido (como lo es la posición, la
cantidad de movimiento y la energía) que intervienen en un determinado problema.
En este caso las variables independientes para describir el movimiento de un fluido son su
vector posición 𝑋⃗ elegido con respecto a alguna ( en algún sistema de referencia elegido ) y el
tiempo t . Se identifica a la particula X con su vector posición 𝑋⃗.
Por ejemplo la posición de una partícula se escribe como 𝑥̅ = (𝑥̅0 , t) y donde 𝑥̅0 es la posición
de la partícula en un tiempo de referencia 𝑡0 que por lo general 𝑡0 = 0.
FORMULACIÓN LAGRANGIANA
Las cantidades intensivas pueden ser expresadas como una función de las coordenadas de la
partícula 𝑋𝑖 y del tiempo presente t. Las cuatro coordenadas 𝑋1 , 𝑋2 , 𝑋3 , 𝑡 también se las
llama “coordenadas materiales” porque identifican una porción de materia o “coordenadas
Lagrangianas”.
Si para un tiempo 𝑡0 una partícula se encuentra en la posición 𝑟̅ = 𝑟̅0 entonces para un
tiempo 𝑡 ≠ 𝑡0 la posición de la partícula estará dada por:
𝑟̅ = 𝑟̅ (𝑟̅0 , t)
donde 𝑟̅0 es la posición de la partícula en 𝑡 = 𝑡0 , donde 𝑟̅0 la etiqueta que identifica a la
partícula.
DESCRIPCION EULERIANA
En la descripción Euleriana, el observador se concentra en lo que sucede en una región del
espacio (un punto espacial) definido por el vector 𝑥̅ por lo tanto las variables independientes
son 𝑥̅ y t.
O sea, me paro en un punto y veo pasar un paquete de fluidos.
La posición de la partícula se fija en un lugar determinado del espacio, simbolizado por 𝑥⃗(en
minúscula).
En esta descripción se pone atención a un lugar determinado del espacio (identificado por y
no a la partícula). Esta es la descripción utilizada en Mecánica de Fluidos, porque en la
caracterización de un flujo lo que interesa es la región del espacio donde las diferentes
partículas ingresan y salen o sea pasan.
Posteriormente a esta región del espacio se la denominará “volumen control” (VC) lo que se
definirá oportunamente.
𝑄 = ∫ (𝑣⃗. 𝑛⃗⃗)𝑑𝐴
𝑠
La ecuación (6.3) determina el “flujo volumétrico” o caudal. Para calcular el flujo de masa
(“flujo másico”) solo es suficiente multiplicar (6.3) por la densidad ρ entonces:
𝑑𝑚
= 𝑚̇ = ∫ 𝜌 (𝑣⃗. 𝑛⃗⃗)𝑑𝐴
𝑑𝑡
𝑠
ECUACIONES DE BALANCE
Para estudiar un determinado sistema debemos “identificar los procesos involucrados”
Magnitudes intensivas: Son aquéllas cuyo valor no depende de la cantidad de materia que
interviene en el fenómeno en cuestión. Están asociadas a formulaciones diferenciales.
Ejemplos: velocidad, temperatura, etc.
Si llamamos a B como una propiedad extensiva cualquiera del fluido que se conserve (masa,
cantidad de movimiento, etc.).
Y si definimos β como la variación de B respecto de la masa
𝑑𝐵
𝛽=
𝑑𝑚
La cantidad total de B en el volumen de control es:
𝑑
Variación de β en el interior del VC: (∫𝑉𝐶 𝛽 𝜌 𝑑𝑉 ) (integral triple)
𝑑𝑡
Por lo tanto, podemos decir que : “El cambio instantáneo de B en el sistema es la suma
de la variación interior más el flujo que sale menos el que entra”.
𝑑𝐵𝑠𝑖𝑠𝑡 𝑑
= (∫ 𝛽 𝜌 𝑑𝑉 ) + ∫ 𝛽 𝜌 (𝑣̅ . 𝑛̅) 𝑑𝐴𝑠𝑎𝑙𝑖𝑑𝑎 − ∫ 𝛽 𝜌 (𝑣̅ . 𝑛̅) 𝑑𝐴𝑒𝑛𝑡𝑟𝑎𝑑𝑎
𝑑𝑡 𝑑𝑡 𝑉𝐶 𝑆𝐶 𝑆𝐶
Esta es la expresión del teorema del transporte de Reynolds para un volumen de control fijo
arbitrario.
Cuando la propiedad B es la masa, cantidad de movimiento, momento cinético o energía,
tenemos las leyes básicas en forma integral.
La forma compacta del transporte de Reynolds es pues:
𝑑𝐵𝑠𝑖𝑠𝑡 𝑑
= (∫ 𝛽 𝜌 𝑑𝑉 ) + ∫ 𝛽 𝜌 (𝑣 ̅. 𝑛 ̅ ) 𝑑𝐴
𝑑𝑡 𝑑𝑡 𝑉𝐶 𝑆𝐶
VELOCIDADES RELATIVAS
Hasta el momento se ha supuesto un volumen de control fijo y que no se mueve. En el caso de
que el volumen de control esté en movimiento, con velocidad uniforme 𝝎̅ , un observador fijo
̅𝒓 definida por:
al VC verá pasar el flujo a una velocidad relativa 𝒗
Velocidad relativa= velocidad del fluido- velocidad volumen control
̅𝒓 = 𝒗
𝒗 ̅−𝝎
̅
Por lo que la ecuación nos queda:
𝒅 𝒅
𝑩𝒔𝒊𝒔𝒕 = ̅𝒓 . 𝒏 ̅ ) 𝒅𝑨
(∫ 𝜷 𝝆 𝒅𝑽) + ∫ 𝜷 𝝆 (𝒗
𝒅𝒕 𝒅𝒕 𝑽𝑪 𝑺𝑪
CONSERVACIÓN DE LA MASA
Considerando en el teorema de transporte de Reynolds que la propiedad B es la masa y que
β = 1 = dm/dm se tiene:
𝑑𝑚 𝑑
= (∫ 𝜌 𝑑𝑉 ) + ∫ 𝜌 (𝑣̅𝑟 . 𝑛 )𝑑𝐴 = 0
𝑑𝑡 𝑑𝑡 𝑉𝐶 𝑆𝐶
Para volumen de control fijo, estado estacionario y para flujo es incompresible, la ecuación
nos queda.
∫ (𝑣̅𝑟 . 𝑛 )𝑑𝐴 = 0
𝑆𝐶
𝑑𝑚 𝑑
= (∫ 𝜌 𝑑𝑉 ) + ∫ 𝜌 (𝑣̅ . 𝑛̅)𝑑𝐴 = 0
𝑑𝑡 𝑑𝑡 𝑉𝐶 𝑆𝐶
𝜕𝜌
∫ 𝑑𝑉 + ∫ 𝜌 (𝑣̅ . 𝑛̅)𝑑𝐴 = 0
𝑉𝐶 𝜕𝑡 𝑆𝐶
𝜕𝜌
∫ 𝑑𝑉 + ∫ ∇. (𝜌𝑣̅ ) 𝑑𝑉 = 0
𝑉𝐶 𝜕𝑡 𝑉𝐶
𝜕𝜌
∫ [ + ∇. (𝜌𝑣̅ )] 𝑑𝑉 = 0
𝑉𝐶 𝜕𝑡
𝜕𝜌
+ ∇. (𝜌𝑣̅ ) = 0
𝜕𝑡
Trabajando la última expresión:
𝜕𝜌
+ 𝑣̅ . ∇𝜌 + 𝜌 ∇. 𝑣̅ = 0
𝜕𝑡
Expresando en función de la derivada material:
𝐷𝜌
+ 𝜌 ∇. 𝑣̅ = 0
𝐷𝑡
Para flujos incompresibles densidad constante:
∇. 𝑣̅ = 0
𝑑(𝑚𝑣̅ ) 𝑑
= (∫ 𝑣̅ 𝜌 𝑑𝑉 ) + ∫ 𝑣̅ 𝜌 (𝑣̅𝑟 . 𝑛 ̅ ) 𝑑𝐴
𝑑𝑡 𝑑𝑡 𝑉𝐶 𝑆𝐶
𝑑
∑ 𝐹̅𝑒𝑥𝑡𝑒𝑟𝑛𝑎𝑠 = (∫ 𝑣̅ 𝜌 𝑑𝑉 ) + ∫ 𝑣̅ 𝜌 (𝑣̅𝑟 . 𝑛 ̅ ) 𝑑𝐴
𝑑𝑡 𝑉𝐶 𝑆𝐶
Fuerzas superficiales, provienen del exterior del sistema, aquí estarían las fuerzas de presión
y las fuerzas viscosas.
Fuerzas volumétricas, provienen del interior del sistema, aquí estaría la fuerza peso.
Como 𝑡𝑛̅ tiene unidades de presión, debe estar relacionado con la presión. En hidrostática se
vio que el vector esfuerzo se expresa como:
𝐷𝑣̅
∫ 𝜌 𝑑𝑉 = ∫ 𝜌𝑔̅ 𝑑𝑉 + ∫ ∇. 𝑇𝑑𝑉
𝑉𝑚(𝑡) 𝐷𝑡 𝑉𝑚(𝑡) 𝑉𝑚(𝑡)
𝐷𝑣̅
∫ [𝜌 − 𝜌𝑔̅ − ∇. 𝑇] 𝑑𝑉 = 0
𝑉𝑚(𝑡) 𝐷𝑡
Por lo tanto:
𝐷𝑣̅
𝜌 − 𝜌𝑔̅ − ∇. 𝑇 = 0
𝐷𝑡
Donde
Si analizamos las fuerzas por unidad de volumen que actúan sobre un elemento de fluido en
movimiento, resulta:
̅ = −𝛁𝒑 + 𝝆𝒈
𝝆𝒂 ̅ 𝒗𝒊𝒔𝒄𝒐𝒔𝒂𝒔
̅+𝒇
Si el fluido es ideal (viscosidad cero) la ecuación nos queda:
̅ = −𝛁𝒑 + 𝝆𝒈
𝝆𝒂 ̅
Que es la ecuación de Euler.
Ahora si consideramos un fluido newtoniano, la fuerza viscosa por unidad de volumen será
igual a:
̅ 𝒗𝒊𝒔𝒄𝒐𝒔𝒂𝒔 = 𝝁 𝛁𝟐 𝒗
𝒇 ̅
Entonces la ecuación toma la forma:
𝝆𝒂 ̅ + 𝝁 𝛁𝟐 𝒗
̅ = −𝛁𝒑 + 𝝆𝒈 ̅
Y esta ecuación se conoce como ecuación de Navier Stokes, que se utiliza generalmente para
encontrar el perfil de velocidades de un campo de flujo.
A continuación, se dan las ecuaciones de la continuidad y de Navier Stokes en coordenadas
cartesianas y cilíndricas:
BALANCE DE ENERGIA
La ecuación de balance de energía se obtiene aplicando el teorema de transporte de
Reynolds a la primera ley dela termodinámica.
En este caso la variable B es ahora la energía y la energía por unidad de masa es
β = dE/dm = e
𝒅𝑬𝒔𝒊𝒔𝒕 𝒅𝒒 𝒅𝑾 𝒅
= + = ̅𝒓 . 𝒏 ̅ ) 𝒅𝑨
(∫ 𝒆 𝝆 𝒅𝑽) + ∫ 𝒆 𝝆 (𝒗
𝒅𝒕 𝒅𝒕 𝒅𝒕 𝒅𝒕 𝑽𝑪 𝑺𝑪
Recordar que q es positivo significa calor comunicado al sistema y que W positivo significa
trabajo realizado por el sistema.
La energía del sistema e por unidad de masa puede ser de varios tipos:
𝑒 = 𝑒𝑖𝑛𝑡𝑒𝑟𝑛𝑎 + 𝑒𝑐𝑖𝑛𝑒𝑡𝑖𝑐𝑎 + 𝑒𝑝𝑜𝑡𝑒𝑛𝑐𝑖𝑎𝑙 + 𝑒𝑜𝑡𝑟𝑎𝑠
En otras energías se refiere a cambios de composición química, reacciones nucleares y
efectos electrostáticos y magnéticos. Aquí despreciamos este término y solo consideramos
los tres primeros.
1
𝑒 = 𝑢 + 𝑣 2 + 𝑔𝑧
2
𝑑𝐸𝑠𝑖𝑠𝑡 𝑑𝑞 𝑑𝑊 𝑑
= + = (∫ 𝑒 𝜌 𝑑𝑉 ) + ∫ 𝑒 𝜌 (𝑣̅𝑟 . 𝑛 ̅ ) 𝑑𝐴
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑉𝐶 𝑆𝐶
𝑑𝑞 𝑑𝑊
+ = ∫ 𝑒 𝜌 (𝑣̅ . 𝑛̅) 𝑑𝐴
𝑑𝑡 𝑑𝑡 𝑆𝐶
Reemplazando:
𝑑𝑄 𝑑𝑊 1
+ = ∫ ( 𝑣 2 + 𝑔𝑧 ) 𝜌 (𝑣̅ . 𝑛̅) 𝑑𝐴
𝑑𝑡 𝑑𝑡 𝑆𝐶 2
𝑑𝑄 𝑑𝑊 1
+ = 𝜌 ∫ ( 𝑣 2 + 𝑔𝑧 ) (𝑣̅ . 𝑛̅) 𝑑𝐴
𝑑𝑡 𝑑𝑡 𝑆𝐶 2
𝑑𝑞 𝑑𝑊 1
+ = 𝜌 ∫ ( 𝑣 2 + 𝑔𝑧 ) 𝑣 𝑑𝐴
𝑑𝑡 𝑑𝑡 𝑆𝐶 2
𝑑𝑞 𝑑𝑊 1
+ = 𝜌 ( 𝑣 2 + 𝑔𝑧 ) 𝑣 ∫ 𝑑𝐴
𝑑𝑡 𝑑𝑡 2 𝑆𝐶
𝑑𝑞 𝑑𝑊 1 1
+ = 𝜌 ( 𝑣 2 + 𝑔𝑧 ) 𝑣 𝐴 (𝑠𝑎𝑙𝑖𝑑𝑎) − 𝜌 ( 𝑣 2 + 𝑔𝑧 ) 𝑣 𝐴 (𝑒𝑛𝑡𝑟𝑎𝑑𝑎)
𝑑𝑡 𝑑𝑡 2 2
Si el sistema no hay flujo de transferencia de calor con los alrededores, sistema adiabático,
𝑑𝑞
= 0.
𝑑𝑡
𝑑𝑊 1 1
= 𝜌 ( 𝑣 2 + 𝑔𝑧 ) 𝑣 𝐴 (𝑠𝑎𝑙𝑖𝑑𝑎) − 𝜌 ( 𝑣 2 + 𝑔𝑧 ) 𝑣 𝐴 (𝑒𝑛𝑡𝑟𝑎𝑑𝑎)
𝑑𝑡 2 2
El trabajo como dijimos, se divide mecánico, trabajo de las fuerzas de presión y trabajo de
esfuerzos viscosos.
Trabajo de las fuerzas de presión:
Considere el flujo de un segmento de un fluido ideal a través de una tubería no uniforme en
un intervalo de tiempo Δt. Durante el intervalo de tiempo, el extremo izquierdo del
segmento se mueve hacia la derecha una distancia Δx1, que es la longitud de la porción
sombreada azul a la izquierda.
Mientras tanto, el extremo derecho del
segmento se mueve hacia la derecha una
distancia Δx2.
La fuerza que ejerce el fluido sobre el extremo
izquierdo tiene una magnitud:
𝐹1 = 𝑃1 𝐴1
El trabajo invertido por esta fuerza sobre el
segmento en un intervalo de tiempo Δt es:
𝑊1 = 𝐹1 𝛥𝑥1 𝑐𝑜𝑠0 = 𝐹1 𝛥𝑥1 = 𝑃1 𝐴1 𝛥𝑥1 = 𝑃1 𝑉1
De forma similar, el trabajo invertido por el fluido a la derecha del segmento en el mismo
intervalo de tiempo Δt es:
𝐹2 = 𝑃2 𝐴2
𝑊2 = 𝐹2 𝛥𝑥2 𝑐𝑜𝑠180° = −𝐹2 𝛥𝑥2 = −𝑃2 𝐴2 𝛥𝑥2 = − 𝑃2 𝑉2
El trabajo de la fuerza debido a la presión de fluido invertido en el segmento por dichas
fuerzas en el mismo intervalo Δt es:
𝑊𝐹 𝑑𝑒 𝑝𝑟𝑒𝑠𝑖𝑜𝑛 = 𝑃1 𝑉1 − 𝑃2 𝑉2
Y como el volumen del elemento de fluido es el mismo:
𝑊𝐹 𝑑𝑒 𝑝𝑟𝑒𝑠𝑖𝑜𝑛 = (𝑃1 − 𝑃2 )𝑉
𝑾𝑭 𝒅𝒆 𝒑𝒓𝒆𝒔𝒊𝒐𝒏 = −∆𝒑𝑽
Derivando en función del tiempo
𝑑𝑊𝑝𝑟𝑒𝑠𝑖𝑜𝑛 𝑑𝑉
= −∆𝑝 = −∆𝑝𝑄
𝑑𝑡 𝑑𝑡
Donde Q es el caudal y ∆𝑝 = 𝑃𝑠𝑎𝑙𝑖𝑑𝑎 − 𝑃𝑒𝑛𝑡𝑟𝑎𝑑𝑎
Ahora veamos el trabajo mecánico que se divide en el trabajo que aporta energía al flujo de
fluido por ejemplo una bomba y en el trabajo que extraemos del flujo por ejemplo el trabajo
de una turbina.
Entonces:
𝑑𝑊𝑝𝑟𝑒𝑠𝑖𝑜𝑛 𝑑𝑊𝑏𝑜𝑚𝑏𝑎 𝑑𝑊𝑡𝑢𝑟𝑏𝑖𝑛𝑎 𝑑𝑊𝑓.𝑣𝑖𝑠𝑐𝑜𝑠𝑎𝑠
+ + +
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝑑𝑡
1 1
= (𝜌 ( 𝑣 2 + 𝑔𝑧 ) 𝑄) − (𝜌 ( 𝑣 2 + 𝑔𝑧 ) 𝑄)
2 𝑠𝑎𝑙𝑖𝑑𝑎 2 𝑒𝑛𝑡𝑟𝑎𝑑𝑎
𝑑𝑊𝑏𝑜𝑚𝑏𝑎 𝑑𝑊𝑡𝑢𝑟𝑏𝑖𝑛𝑎 𝑑𝑊𝑓.𝑣𝑖𝑠𝑐𝑜𝑠𝑎𝑠
−∆𝑝𝑄 + + +
𝑑𝑡 𝑑𝑡 𝑑𝑡
1 2 1
= 𝜌 ( 𝑣 + 𝑔𝑧 ) 𝑄 (𝑠𝑎𝑙𝑖𝑑𝑎) − 𝜌 ( 𝑣 2 + 𝑔𝑧 ) 𝑄 (𝑒𝑛𝑡𝑟𝑎𝑑𝑎)
2 2
1 1 𝑑𝑊𝑏𝑜𝑚𝑏𝑎 𝑑𝑊𝑡𝑢𝑟𝑏𝑖𝑛𝑎
𝜌 ( 𝑣 2 + 𝑔𝑧 ) 𝑄 (𝑒𝑛𝑡𝑟𝑎𝑑𝑎) − 𝜌 ( 𝑣 2 + 𝑔𝑧 ) 𝑄 (𝑠𝑎𝑙𝑖𝑑𝑎) − ∆𝑝𝑄 + +
2 2 𝑑𝑡 𝑑𝑡
𝑑𝑊𝑓.𝑣𝑖𝑠𝑐𝑜𝑠𝑎𝑠
+ =0
𝑑𝑡
Dividiendo todo por Q
1 1 𝑑𝑊𝑏𝑜𝑚𝑏𝑎
𝜌 ( 𝑣 2 + 𝑔𝑧 ) − 𝜌 ( 𝑣 2 + 𝑔𝑧 ) − (𝑃𝑠𝑎𝑙𝑖𝑑𝑎 − 𝑃𝑒𝑛𝑡𝑟𝑎𝑑𝑎 ) +
2 (𝑒𝑛𝑡𝑟𝑎𝑑𝑎) 2 (𝑠𝑎𝑙𝑖𝑑𝑎) 𝑄𝑑𝑡
𝑑𝑊𝑡𝑢𝑟𝑏𝑖𝑛𝑎 𝑑𝑊𝑓.𝑣𝑖𝑠𝑐𝑜𝑠𝑎𝑠
+ + =0
𝑄𝑑𝑡 𝑄𝑑𝑡
Dividiendo por 𝜌𝑔
𝑣2 𝑝 𝑣2 𝑝 𝑑𝑊𝑏𝑜𝑚𝑏𝑎 𝑑𝑊𝑡𝑢𝑟𝑏𝑖𝑛𝑎 𝑑𝑊𝑓.𝑣𝑖𝑠𝑐𝑜𝑠𝑎𝑠
( +𝑧+ ) −( +𝑧+ ) + + + =0
2𝑔 𝜌𝑔 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 2𝑔 𝜌𝑔 𝑒𝑛𝑡𝑟𝑎𝑑 𝜌𝑔𝑄𝑑𝑡 𝜌𝑔𝑄𝑑𝑡 𝜌𝑔𝑄𝑑𝑡
𝑑𝑊𝑏𝑜𝑚𝑏𝑎
= +ℎ𝑏𝑜𝑚𝑏𝑎 (𝑒𝑛𝑒𝑟𝑔𝑖𝑎 𝑔𝑎𝑛𝑎𝑑𝑎 𝑝𝑜𝑟 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑒𝑠 𝑠𝑜𝑙𝑖𝑑𝑎𝑠 𝑚𝑜𝑣𝑖𝑙𝑒𝑠)
𝜌𝑔𝑄𝑑𝑡
𝑑𝑊𝑡𝑢𝑟𝑏𝑖𝑛𝑎
= −ℎ𝑡𝑢𝑟𝑏𝑖𝑛𝑎 (𝑒𝑛𝑒𝑟𝑔𝑖𝑎 𝑝𝑒𝑟𝑑𝑖𝑑𝑎 𝑝𝑜𝑟 𝑠𝑢𝑝𝑒𝑟𝑓𝑖𝑐𝑖𝑒𝑠 𝑠𝑜𝑙𝑖𝑑𝑎𝑠 𝑚𝑜𝑣𝑖𝑙𝑒𝑠)
𝜌𝑔𝑄𝑑𝑡
𝑑𝑊𝑓 𝑣𝑖𝑠𝑐𝑜𝑠𝑎𝑠
= −ℎ𝑓𝑟𝑖𝑐𝑐𝑖𝑜𝑛 (𝑒𝑛𝑒𝑟𝑔𝑖𝑎 𝑝𝑒𝑟𝑑𝑖𝑑𝑎 𝑝𝑜𝑟 𝑓𝑟𝑖𝑐𝑐𝑖𝑜𝑛)
𝜌𝑔𝑄𝑑𝑡
𝑣2 𝑝 𝑣2 𝑝
( +𝑧+ ) −( +𝑧+ ) + ℎ𝑏𝑜𝑚𝑏𝑎 − ℎ𝑡𝑢𝑟𝑏𝑖𝑛𝑎 − ℎ𝑓𝑟𝑖𝑐𝑐𝑖𝑜𝑛 = 0
2𝑔 𝜌𝑔 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 2𝑔 𝜌𝑔 𝑠𝑎𝑙𝑖𝑑𝑎
𝑣2 𝑝 𝑣2 𝑝
( +𝑧+ ) =( +𝑧+ ) − ℎ𝑏𝑜𝑚𝑏𝑎 + ℎ𝑡𝑢𝑟𝑏𝑖𝑛𝑎 + ℎ𝑓𝑟𝑖𝑐𝑐𝑖𝑜𝑛
2𝑔 𝜌𝑔 𝑒𝑛𝑡𝑟𝑎𝑑𝑎 2𝑔 𝜌𝑔 𝑠𝑎𝑙𝑖𝑑𝑎
𝒑𝟏 𝒗𝟐𝟏 𝒑𝟐 𝒗𝟐𝟐
+ + 𝒛𝟏 = + + 𝒛𝟐 + 𝒉𝒇𝒓𝒊𝒄𝒄𝒊𝒐𝒏 + 𝒉𝒕𝒖𝒓𝒃𝒊𝒏𝒂 − 𝒉𝒃𝒐𝒎𝒃𝒂
𝝆𝒈 𝟐𝒈 𝝆𝒈 𝟐𝒈
ECUACION DE MOVIMIENTO
Coordenadas Cilíndricas
𝜕𝑣𝑟 𝜕𝑣𝑟 𝑣𝜃 𝜕𝑣𝑟 𝑣𝜃2 𝜕𝑣𝑟 𝜕𝑝 1 𝜕 𝜕𝑣𝑟 𝑣𝑟 1 𝜕 2 𝑣𝑟 2 𝜕𝑣𝑟 𝜕 2 𝑣𝑟
𝜌( + 𝑣𝑟 + − + 𝑣𝑧 ) = 𝜌𝑔𝑟 − +𝜇( (𝑟 )− 2+ 2 − + )
𝜕𝑡 𝜕𝑟 𝑟 𝜕𝜃 𝑟 𝜕𝑧 𝜕𝑟 𝑟 𝜕𝑟 𝜕𝑟 𝑟 𝑟 𝜕𝜃 2 𝑟 2 𝜕𝜃 𝜕𝑧 2
Martin Rodriguez