Solución:: Demostrar La Ecuación de Navier Stokes: +

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 4

Demostrar la ecuación de Navier Stokes:

𝐷𝑣⃗
𝜌 = 𝜇∇2 𝑣⃗ − ∇𝑝⃗ + 𝜌𝑔⃗
𝐷𝑡
Solución:
Partiendo de la ecuación de movimiento para sus tres componentes en coordenadas rectangulares
𝜕 𝜕 𝜕 𝜕 𝜕 𝜕 𝜕 𝜕𝑝
𝜌𝑣𝑥 = − ( 𝜌𝑣𝑥 𝑣𝑥 + 𝜌𝑣𝑦 𝑣𝑥 + 𝜌𝑣𝑧 𝑣𝑥 ) − ( 𝜏𝑥𝑥 + 𝜏𝑦𝑥 + 𝜏𝑧𝑥 ) − + 𝜌𝑔𝑥 (3.2 − 5)
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥
𝜕 𝜕 𝜕 𝜕 𝜕 𝜕 𝜕 𝜕𝑝
𝜌𝑣𝑦 = − ( 𝜌𝑣𝑥 𝑣𝑦 + 𝜌𝑣𝑦 𝑣𝑦 + 𝜌𝑣𝑧 𝑣𝑦 ) − ( 𝜏𝑥𝑦 + 𝜏𝑦𝑦 + 𝜏𝑧𝑦 ) − + 𝜌𝑔𝑦 (3.2 − 6)
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑦
𝜕 𝜕 𝜕 𝜕 𝜕 𝜕 𝜕 𝜕𝑝
𝜌𝑣𝑧 = − ( 𝜌𝑣𝑥 𝑣𝑧 + 𝜌𝑣𝑦 𝑣𝑧 + 𝜌𝑣𝑧 𝑣𝑧 ) − ( 𝜏𝑥𝑧 + 𝜏𝑦𝑧 + 𝜏𝑧𝑧 ) − + 𝜌𝑔𝑧 (3.2 − 7)
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑧

Tenemos que:

𝜌𝑣⃗ = (𝜌𝑣𝑥 , 𝜌𝑣𝑦 , 𝜌𝑣𝑧 ) 𝑣𝑒𝑐𝑡𝑜𝑟 𝑣𝑒𝑙𝑜𝑐𝑖𝑑𝑎𝑑 𝑚á𝑠𝑖𝑐𝑎

𝜕𝑝 𝜕𝑝 𝜕𝑝
⃗∇⃗𝑝 = ( , , ) 𝑔𝑟𝑎𝑑𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑝𝑟𝑒𝑠𝑖ó𝑛
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜌𝑣𝑥 𝑣𝑥 𝜌𝑣𝑥 𝑣𝑦 𝜌𝑣𝑥 𝑣𝑧
[𝜌𝑣⃗𝑣⃗] = (𝜌𝑣𝑦 𝑣𝑥 𝜌𝑣𝑦 𝑣𝑦 𝜌𝑣𝑦 𝑣𝑧 ) 𝑡𝑒𝑛𝑠𝑜𝑟 𝑓𝑜𝑟𝑚𝑎𝑑𝑜 𝑝𝑜𝑟 𝑒𝑙 𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑜 𝑑𝑖á𝑑𝑖𝑐𝑜 𝑑𝑒 𝑣⃗𝑣⃗ (𝐴. 4 − 3)
𝜌𝑣𝑧 𝑣𝑥 𝜌𝑣𝑧 𝑣𝑦 𝜌𝑣𝑧 𝑣𝑧
𝜏𝑥𝑥 𝜏𝑦𝑥 𝜏𝑧𝑥
[𝜏⃗] = (𝜏𝑥𝑦 𝜏𝑦𝑦 𝜏𝑧𝑦 ) 𝑡𝑒𝑛𝑠𝑜𝑟 𝑒𝑠𝑓𝑢𝑒𝑟𝑧𝑜 𝑐𝑜𝑟𝑡𝑎𝑛𝑡𝑒
𝜏𝑥𝑧 𝜏𝑦𝑧 𝜏𝑧𝑧

Combinando las tres direcciones de la ecuación de movimiento, expresando mediante operadores diferenciales:
𝜕
⃗⃗ ∙ 𝜌𝑣⃗𝑣⃗] − ∇
𝜌𝑣⃗ = −[∇ ⃗⃗𝑝 − [∇
⃗⃗ ∙ 𝜏⃗] + 𝜌𝑔⃗
𝜕𝑡
Partiendo de la ecuación (3.2 − 5), considerando 𝜌 constante:
𝜕𝑣𝑥 𝜕 𝜕 𝜕 𝜕 𝜕 𝜕 𝜕𝑝
𝜌 + 𝜌 ( 𝑣𝑥 𝑣𝑥 + 𝑣𝑦 𝑣𝑥 + 𝑣𝑧 𝑣𝑥 ) = − ( 𝜏𝑥𝑥 + 𝜏𝑦𝑥 + 𝜏𝑧𝑥 ) − + 𝜌𝑔𝑥
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥
𝑑 𝑑 𝑑
𝑆𝑒𝑎𝑛 𝑢(𝑥) 𝑦 𝑣(𝑥): (𝑢𝑣) = 𝑢 𝑣 + 𝑣 𝑢 (𝑟𝑒𝑔𝑙𝑎 𝑑𝑒𝑙 𝑝𝑟𝑜𝑑𝑢𝑐𝑡𝑜)
𝑑𝑥 𝑑𝑥 𝑑𝑥
𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑥 𝜕𝑣𝑧 𝜕 𝜕 𝜕 𝜕𝑝
𝜌( + 𝑣𝑥 + 𝑣𝑥 + 𝑣𝑦 + 𝑣𝑥 + 𝑣𝑧 + 𝑣𝑥 ) = − ( 𝜏𝑥𝑥 + 𝜏𝑦𝑥 + 𝜏𝑧𝑥 ) − + 𝜌𝑔𝑥
𝜕𝑡 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥

𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕 𝜕 𝜕 𝜕𝑝


𝜌( + 𝑣𝑥 + 𝑣𝑦 + 𝑣𝑧 + 𝑣𝑥 ( + + )) = − ( 𝜏𝑥𝑥 + 𝜏𝑦𝑥 + 𝜏𝑧𝑥 ) − + 𝜌𝑔𝑥
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥

𝐷𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥


= + 𝑣𝑥 + 𝑣𝑦 + 𝑣𝑧 (3.0 − 2) 𝐷𝑒𝑟𝑖𝑣𝑎𝑑𝑎 𝑠𝑢𝑏𝑠𝑡𝑎𝑛𝑐𝑖𝑎𝑙
𝑑𝑡 𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
(∇ ∙ 𝑣⃗) = + + = 0 (𝜌 = 𝑐𝑡𝑡𝑒) (3.1 − 7) 𝐸𝑐. 𝑐𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑑𝑎𝑑
𝜕𝑥 𝜕𝑦 𝜕𝑧
𝐷𝑣𝑥 𝜕 𝜕 𝜕 𝜕𝑝
𝜌 = − ( 𝜏𝑥𝑥 + 𝜏𝑦𝑥 + 𝜏𝑧𝑥 ) − + 𝜌𝑔𝑥 (1)
𝑑𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥
Análogamente para las ecuaciones (3.2 − 6) y (3.2 − 7):
𝐷𝑣𝑦 𝜕 𝜕 𝜕 𝜕𝑝
𝜌 = − ( 𝜏𝑥𝑦 + 𝜏𝑦𝑦 + 𝜏𝑧𝑦 ) − + 𝜌𝑔𝑦 (2)
𝑑𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑦
𝐷𝑣𝑧 𝜕 𝜕 𝜕 𝜕𝑝
𝜌 = − ( 𝜏𝑥𝑧 + 𝜏𝑦𝑧 + 𝜏𝑧𝑧 ) − + 𝜌𝑔𝑧 (3)
𝑑𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑧

Combinando las ecuaciones (1), (2) (3):


𝐷𝑣⃗
𝜌 ⃗⃗𝑝 − [∇
= −∇ ⃗⃗ ∙ 𝜏⃗] + 𝜌𝑔⃗ (𝛼)
𝑑𝑡
Componentes del tensor esfuerzo cortante en coordenadas rectangulares tabla 3.4-5:
𝜕𝑣𝑥 2
𝜏𝑥𝑥 = −𝜇 [2 − (∇ ∙ 𝑣⃗)]
𝜕𝑥 3
𝜕𝑣𝑦 2
𝜏𝑦𝑦 = −𝜇 [2 − (∇ ∙ 𝑣⃗)]
𝜕𝑦 3
𝜕𝑣𝑧 2
𝜏𝑧𝑧 = −𝜇 [2 − (∇ ∙ 𝑣⃗)]
𝜕𝑧 3
𝜕𝑣𝑥 𝜕𝑣𝑦
𝜏𝑥𝑦 = 𝜏𝑦𝑥 = −𝜇 [ + ]
𝜕𝑦 𝜕𝑥

𝜕𝑣𝑦 𝜕𝑣𝑧
𝜏𝑦𝑧 = 𝜏𝑧𝑦 = −𝜇 [ + ]
𝜕𝑧 𝜕𝑦
𝜕𝑣𝑧 𝜕𝑣𝑥
𝜏𝑧𝑥 = 𝜏𝑥𝑧 = −𝜇 [ + ]
𝜕𝑥 𝜕𝑧
Ahora, como consideramos 𝜌 = 𝑐𝑡𝑡𝑒:
(∇ ∙ 𝑣⃗) = 0 (3.1 − 7)

Simplificando algunos componentes del tensor esfuerzo cortante:


𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥
𝜏𝑥𝑥 = −𝜇 [2 ] = −𝜇 [ + ]
𝜕𝑥 𝜕𝑥 𝜕𝑥
𝜕𝑣𝑦 𝜕𝑣𝑦 𝜕𝑣𝑦
𝜏𝑦𝑦 = −𝜇 [2 ] = −𝜇 [ + ]
𝜕𝑦 𝜕𝑦 𝜕𝑦
𝜕𝑣𝑧 𝜕𝑣𝑧 𝜕𝑣𝑧
𝜏𝑧𝑧 = −𝜇 [2 ] = −𝜇 [ + ]
𝜕𝑧 𝜕𝑧 𝜕𝑧
Ahora reemplazamos los componentes del tensor esfuerzo cortante en las ecuaciones (1), (2) y (3)
respectivamente:
𝐷𝑣𝑥 𝜕 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕 𝜕𝑣𝑧 𝜕𝑣𝑥 𝜕𝑝
𝜌 = − ( {−𝜇 [ + ]} + {−𝜇 [ + ]} + {−𝜇 [ + ]}) − + 𝜌𝑔𝑥
𝑑𝑡 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑥 𝜕𝑧 𝜕𝑥 𝜕𝑧 𝜕𝑥

𝐷𝑣𝑦 𝜕 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕 𝜕𝑣𝑦 𝜕𝑣𝑦 𝜕 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕𝑝


𝜌 = − ( {−𝜇 [ + ]} + {−𝜇 [ + ]} + {−𝜇 [ + ]}) − + 𝜌𝑔𝑦
𝑑𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑦 𝜕𝑦

𝐷𝑣𝑧 𝜕 𝜕𝑣𝑧 𝜕𝑣𝑥 𝜕 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕 𝜕𝑣𝑧 𝜕𝑣𝑧 𝜕𝑝


𝜌 = − ( {−𝜇 [ + ]} + {−𝜇 [ + ]} + {−𝜇 [ + ]}) − + 𝜌𝑔𝑧
𝑑𝑡 𝜕𝑥 𝜕𝑥 𝜕𝑧 𝜕𝑦 𝜕𝑧 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑧 𝜕𝑧

Ordenando
𝐷𝑣𝑥 𝜕 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕 𝜕𝑣𝑧 𝜕𝑣𝑥 𝜕𝑝
𝜌 = 𝜇( [ + ]+ [ + ]+ [ + ]) − + 𝜌𝑔𝑥
𝑑𝑡 𝜕𝑥 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑥 𝜕𝑧 𝜕𝑥 𝜕𝑧 𝜕𝑥

𝐷𝑣𝑦 𝜕 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕 𝜕𝑣𝑦 𝜕𝑣𝑦 𝜕 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕𝑝


𝜌 = 𝜇( [ + ]+ [ + ]+ [ + ]) − + 𝜌𝑔𝑦
𝑑𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑦 𝜕𝑦

𝐷𝑣𝑧 𝜕 𝜕𝑣𝑧 𝜕𝑣𝑥 𝜕 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕 𝜕𝑣𝑧 𝜕𝑣𝑧 𝜕𝑝


𝜌 = 𝜇( [ + ]+ [ + ]+ [ + ]) − + 𝜌𝑔𝑧
𝑑𝑡 𝜕𝑥 𝜕𝑥 𝜕𝑧 𝜕𝑦 𝜕𝑧 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑧 𝜕𝑧

Hallando las segundas derivadas y aplicando el teorema de Schwarz de las derivadas cruzadas:

𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑦
= ; = ; … (𝑇𝑒𝑜𝑟𝑒𝑚𝑎 𝑑𝑒 𝑆𝑐ℎ𝑤𝑎𝑟𝑧)
𝜕𝑥𝜕𝑦 𝜕𝑦𝜕𝑥 𝜕𝑧𝜕𝑦 𝜕𝑦𝜕𝑧

𝐷𝑣𝑥 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑧 𝜕 2 𝑣𝑥 𝜕𝑝
𝜌 = 𝜇( 2 + 2 + 2 + + + 2 )− + 𝜌𝑔𝑥 (4)
𝑑𝑡 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑥𝜕𝑦 𝜕𝑥𝜕𝑧 𝜕𝑧 𝜕𝑥

𝐷𝑣𝑦 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑧 𝜕𝑝
𝜌 = 𝜇( + 2 + 2 + 2 + 2 + )− + 𝜌𝑔𝑦 (5)
𝑑𝑡 𝜕𝑦𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑦𝜕𝑧 𝜕𝑦

𝐷𝑣𝑧 𝜕 2 𝑣𝑧 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑧 𝜕 2 𝑣𝑧 𝜕 2 𝑣𝑧 𝜕𝑝
𝜌 = 𝜇( 2 + + + 2 + 2 + 2 )− + 𝜌𝑔𝑧 (6)
𝑑𝑡 𝜕𝑥 𝜕𝑧𝜕𝑥 𝜕𝑧𝜕𝑦 𝜕𝑦 𝜕𝑧 𝜕𝑧 𝜕𝑧

𝜕 2 𝑠 𝜕 2𝑠 𝜕 2 𝑠
∇2 𝑠 = + + (𝐿𝑎𝑝𝑙𝑎𝑐𝑖𝑎𝑛𝑎 𝑑𝑒 𝑢𝑛 𝑐𝑎𝑚𝑝𝑜 𝑒𝑠𝑐𝑎𝑙𝑎𝑟)
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2

𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥
+ + = ∇2 𝑣𝑥
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2

𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑦
+ + = ∇2 𝑣𝑦
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2
𝜕 2 𝑣𝑧 𝜕 2 𝑣𝑧 𝜕 2 𝑣𝑧
+ + = ∇2 𝑣𝑧
𝜕𝑥 2 𝜕𝑦 2 𝜕𝑧 2
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧
(∇ ∙ 𝑣⃗) = + + = 0 (𝜌 = 𝑐𝑡𝑡𝑒) 𝐸𝑐. 𝑐𝑜𝑛𝑡𝑖𝑛𝑢𝑖𝑑𝑎𝑑
𝜕𝑥 𝜕𝑦 𝜕𝑧

𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑧 𝜕 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕


2
+ + = ( + + )= (0) = 0
𝜕𝑥 𝜕𝑥𝜕𝑦 𝜕𝑥𝜕𝑧 𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑥

𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑧 𝜕 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕


+ 2
+ = ( + + )= (0) = 0
𝜕𝑦𝜕𝑥 𝜕𝑦 𝜕𝑦𝜕𝑧 𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑦

𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑧 𝜕 𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧 𝜕


+ + = ( + + ) = (0) = 0
𝜕𝑧𝜕𝑥 𝜕𝑧𝜕𝑦 𝜕𝑧 2 𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑧 𝜕𝑦

Reemplazando en las ecuaciones (4), (5) y (6):


𝐷𝑣𝑥 𝜕𝑝
𝜌 = 𝜇∇2 𝑣𝑥 − + 𝜌𝑔𝑥
𝑑𝑡 𝜕𝑥
𝐷𝑣𝑦 𝜕𝑝
𝜌 = 𝜇∇2 𝑣𝑦 − + 𝜌𝑔𝑦
𝑑𝑡 𝜕𝑦
𝐷𝑣𝑧 𝜕𝑝
𝜌 = 𝜇∇2 𝑣𝑧 − + 𝜌𝑔𝑧
𝑑𝑡 𝜕𝑧
Combinando los tres componentes 𝑥, 𝑦, 𝑧 en una sola ecuación vectorial:
𝐷𝑣𝑥 𝐷𝑣𝑦 𝐷𝑣𝑧 𝜕𝑝 𝜕𝑝 𝜕𝑝
(𝜌 ,𝜌 ,𝜌 ) = (𝜇∇2 𝑣𝑥 − + 𝜌𝑔𝑥 , 𝜇∇2 𝑣𝑦 − + 𝜌𝑔𝑦 , 𝜇∇2 𝑣𝑧 − + 𝜌𝑔𝑧 )
𝑑𝑡 𝑑𝑡 𝑑𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
Por propiedades de vectores:
𝐷 𝜕𝑝 𝜕𝑝 𝜕𝑝
𝜌 (𝑣𝑥 , 𝑣𝑦 , 𝑣𝑧 ) = 𝜇(∇2 𝑣𝑥 , ∇2 𝑣𝑦 , 𝜌∇2 𝑣𝑧 ) − ( , , ) + 𝜌(𝑔𝑥 , 𝑔𝑦 , 𝑔𝑧 )
𝑑𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧

(𝑣𝑥 , 𝑣𝑦 , 𝑣𝑧 ) = 𝑣⃗ ; (𝑔𝑥 , 𝑔𝑦 , 𝑔𝑧 ) = 𝑔⃗

𝜕𝑠 𝜕𝑠 𝜕𝑠
⃗⃗𝑠 = ( , , ) (𝐺𝑟𝑎𝑑𝑖𝑒𝑛𝑡𝑒 𝑑𝑒 𝑢𝑛 𝑐𝑎𝑚𝑝𝑜 𝑒𝑠𝑐𝑎𝑙𝑎𝑟)

𝜕𝑥 𝜕𝑦 𝜕𝑧

∇2 𝑣⃗ = (∇2 𝑣𝑥 , ∇2 𝑣𝑦 , 𝜌∇2 𝑣𝑧 ) (𝐿𝑎𝑝𝑙𝑎𝑐𝑖𝑎𝑛𝑎 𝑑𝑒 𝑢𝑛 𝑐𝑎𝑚𝑝𝑜 𝑣𝑒𝑐𝑡𝑜𝑟𝑖𝑎𝑙)

𝐷𝑣⃗
𝜌 = 𝜇∇2 𝑣⃗ − ⃗∇⃗𝑝 + 𝜌𝑔⃗ 𝐸𝑐. 𝑁𝑎𝑣𝑖𝑒𝑟 𝑆𝑡𝑜𝑘𝑒𝑠
𝑑𝑡
Extra, demostrar la ecuación de Euler
Solución:
A partir de la ecuación (𝛼):
𝐷𝑣⃗
𝜌 ⃗⃗𝑝 − [∇
= −∇ ⃗⃗ ∙ 𝜏⃗] + 𝜌𝑔⃗ (𝛼)
𝑑𝑡
Consideramos que no existen fuerzas viscosas:
⃗⃗ ∙ 𝜏⃗] = 0
[∇

Reemplazando:
𝐷𝑣⃗
𝜌 ⃗⃗𝑝 + 𝜌𝑔⃗ 𝐸𝑐. 𝑑𝑒 𝐸𝑢𝑙𝑒𝑟
= −∇
𝑑𝑡

También podría gustarte