Medios Continuos Teor A
Medios Continuos Teor A
Medios Continuos Teor A
Aviso a navegantes:
Este documento tiene intención de redactar las preguntas de MMC desde el
principio de los tiempos hasta septiembre de 2022 para memorizar únicamente
(si, hay 34 paginas).
Suerte.
Tomaa
1
Contents
1 Pregunta 1: Trayectorias 3
3 Pregunta 3: Hidrostática. 8
2
1 Pregunta 1: Trayectorias
Desarrolle los siguientes conceptos: definición, obtención y posible
consecuencias:
1. Trayectoria de una partı́cula fluida
2. Lineas de corriente
3. Lineas de humo (solo definición)
Trayectoria: un punto material situado en ⃗x = ⃗xT (⃗x0 , t), con ⃗xT expresión
que nos muestra su posición en cada instante t, es decir, la trayectoria. Esta se
obtiene mediante el sistema de ec. diferenciales:
d⃗x
= ⃗v (⃗x, t)
dt
Con ⃗v el campo de velocidades.
h1 dx1 h2 dx2 h3 dx3
dt = = =
v1 (x1 , x2 , x3 , t) v2 (x1 , x2 , x3 , t) v3 (x1 , x2 , x3 , t)
Con lo que:
x1 = xT (x10 , x20 , x30 , t)
x2 = xT (x10 , x20 , x30 , t)
x3 = xT (x10 , x20 , x30 , t)
Eliminando t entre dos de estas expresiones, se obtiene una superficie que
representa la senda o curva descrita por el punto material en su movimiento.
Lineas de corriente: Las lineas de corriente en t son las lineas del campo de
velocidades en dicho instante. Estas serán las lineas tangentes en cada punto al
vector velocidad local.
3
• Si el movimiento fluido es estacionario, es decir que ⃗v = ⃗v (⃗x) por lo que
las ecuaciones diferenciales serás equivalentes, al no depender la velocidad
del tiempo.
• El segundo caso en el que coinciden se da cuando el campo de velocidades
es de la forma ⃗v = f (t)⃗v ′ (⃗x)
4
2 Pregunta 2: Fuerzas y tensor de esfuerzos.
Fuerzas de volumen y superficie sobre una partı́cula fluida. El tensor
de esfuerzos; deducción y criterio de signos (no se pide que demuestre
la simetrı́a de este tensor).
La fuerza másica que actúa sobre una partı́cula fluida δV será: ρf⃗m (⃗x, t)δV ,
donde:
⃗
dΩ
f⃗m (⃗x, t) = ⃗g − ⃗a0 − ⃗ ∧ (Ω ∧ ⃗x) − 2Ω
∧ ⃗x − Ω ⃗ ∧ ⃗v
dt
Donde f⃗m (⃗x, t) es la magnitud del campo de fuerzas por unidad de masa, a0
es la aceleración del origen del sistema de referencia no inercial y Ω ⃗ su velocidad
angular de giro.
Algunas de las fuerzas másicas derivan del gradiente de una función poten-
cial escalar f⃗m (⃗x, t) = −∇U .
5
Si el fluido es un lı́quido, aparecen contribuciones adicionales debidas a la
fuerza que se ejercen entre sı́ las moléculas situadas a un lado y otro de la su-
perficie.
Como las fuerzas de corto alcance decaen con d, su resultante sobre una
partı́cula fluida es proporcional a la superficie y no al volumen de dicha partı́cula
fluida δV , tal que (δV )1/3 >> d.
Dada la dependencia de f⃗n (⃗n, ⃗x, t) con la orientación, para calcular el campo
de esfuerzos f⃗n en función de una posición ⃗x y un instante t determinados, nece-
sitaremos una doble infinitud de fuerzas. Vamos a demostrar que en realidad
el número de grados de libertad es mucho menor, de modo que la dependencia
de f⃗n con se puede expresar en función de las seis componentes de un tensor
simétrico τ denominado tensor de esfuerzos.
6
Si llamamos f⃗i = (τ1 , τ2 , τ3 ) al esfuerzo que actúa sobre la cara perpen-
dicular a xi de área dAi , la fuerza que actúa sobre la cara perpendicular a
⃗n = (n1 , n2 , n3 ), cumplirá:
Donde:
τ11 τ12 τ13
τ = τ21 τ22 τ23
τ31 τ32 τ33
Es el denominado tensor de esfuerzos. (τii esfuerzos normales, τij esfuerzos tan-
genciales) La ecuación f⃗n = n1 f⃗1 + n2 f⃗2 + n3 f⃗3 nos dice que el esfuerzo (fuerza
por unidad de superficie) en la dirección ⃗n se puede expresar en función de las
nueve componentes del tensor de esfuerzos τ que son función de la posición y
del tiempo.
Este resultado puede ser más simplificado debido a la simetrı́a del tensor de
esfuerzos viscosos. Para ello aplicamos la ecuación de conservación del momento
cinético a un elemento fluido que forma inicialmente un cubo de arista dx.
Que en el limite de dx → 0 la ecuación se reduce al equilibrio del momento
de las fuerzas superficiales. Considerando por ejemplo la componente vertical,
la ecuación de equilibrio nos da dx3 (τ21 − τ12 ) = 0, indicando por tanto que
τ12 = τ21 , y de manera similar para los 4 restantes esfuerzos tangenciales, con lo
que se demuestra la simetrı́a del tensor de esfuerzos. Pudiéndose expresar por
tanto el esfuerzo superficial f⃗n (⃗n, ⃗x, t) como:
τ11 τ12 τ13
f⃗n = τ · ⃗n = τ12 τ22 τ23 · ⃗n
τ13 τ23 τ33
Donde los esfuerzos normales τ11 , τ22 y τ33 y los esfuerzos tangenciales τ12 ,
τ13 y τ23 son función de la posición y del tiempo, con una dependencia con el
estado termodinámico local del fluido y con la velocidad de deformación.
7
3 Pregunta 3: Hidrostática.
Nota: añado la 3º ecuación de Navier stokes, aunque en muchos años solo pre-
gunta por la primera y la segunda. Encuentra las ecuaciones de Navier
Stokes de balance de masa y de cantidad de movimiento para un
liquido en equilibrio mecánico (hidrostática) situado en un campo
gravitatorio en la superficie de la tierra.
−∇p + ρf⃗m = 0
∂e
ρ = ∇ · (k∇T ) + Qr + Qq
∂t
La primera nos indica que, para un fluido en reposo, la densidad no puede
variar con el tiempo, aún cuando si pueden existir variaciones espaciales. Por
otra parte, al ser nula tanto la aceleración del fluido como la velocidad de de-
formación, la ecuación de cantidad de movimiento (la segunda) se reduce al
equilibrio entre las fuerzas de presión y las fuerzas másicas.
dΩ
f⃗m (⃗x, t) = ⃗g − a⃗0 − ∧ ⃗x − Ω ∧ (Ω ∧ ⃗x)
dt
Donde la fuerza de Coriolis queda excluida de las fuerzas de inercia (de-
bido al movimiento no inercial de los sistemas de referencia) por ser v = 0. La
ecuación de la energı́a indica que el aumento local de energı́a interna se produce
a través de la adicción del calor de conducción, reacción quı́mica o radiación.
∇p/ρ = f⃗m
Con lo que las únicas fuerzas compatibles con el reposo son aquellas que derivan
de un potencial, tenemos entonces un potencial tal que:
⃗ ∧ ⃗x) · (Ω
U = −⃗g · ⃗x + ⃗a0 · ⃗x − (Ω ⃗ ∧ ⃗x)/2
Como:
dΩ
f⃗m = −∇U − ∧ ⃗x
dt
8
La condición anterior indica que dΩ
dt = 0 ha de cumplirse para el equilibrio
mecánico. Integrando la primera ecuación obtenemos:
p + ρU = C(t)
p + ρgz = pa
9
Z
⃗ =−
M ∇ ∧ [(pa − ρgz)(⃗x − x⃗0 )]dV
V
Estas dos ecuaciones constituyen el enunciado matemático del principio de
Arquı́medes: La fuerza de flotabilidad que ejerce un lı́quido sobre un cuerpo
sumergido en él es igual y contraria al peso del volumen del lı́quido desalojado,
y su punto de aplicación es el centro de gravedad del liquido desalojado.
10
4 Pregunta 4: Elemento cúbico fluido.
Deformación de un elemento cúbico de fluido (se pretende que en-
cuentre el significado fı́sico de cada una de las componentes de los
tensores velocidad de deformación unitaria y de rotación)+ variación
encuentre una expresión para la velocidad de dilatación cúbica a par-
tid del tensor de velocidades de deformación unitaria. Indica las
hipótesis y/o condiciones que utilice y las consecuencias del resultado
obtenido.
Significados:
11
• El angulo entre los ejes (inicialmente recto) variará una cantidad (∂1⃗v2 +
∂2⃗v1 )dt, por tanto ∂1⃗v2 +∂2⃗v1 es la velocidad de dilatación angular formado
por las lineas fluidas en las direcciones a y b donde a = 1, 2, 3 y b = 1, 2, 3 y
a ̸= b. Es decir, los términos no diagonales de T d es la unidad de velocidad
de dilatación angular.
• Si tomamos 2 ejes, por ejemplo x1 y x2 , han sufrido giros en su plano
de ángulos ∂1 v2 y ∂2 v1 , los cuales son la velocidad media de rotación
1/2(∂1 v2 − ∂2 v1 ) = ω3 /2, es decir, velocidad angular en torno al eje per-
pendicular al plano. Por tanto, T r es la rotación del elemento fluido.
La distorsión del elemento fluido cúbico ha dado también lugar a una posible
variación de su volumen. Si el volumen del cubo inicial es dl3 , el volumen del
paralelepipedo que se forma al pasar un tiempo dt pasa a ser:
1 + ∂1 v1 dt ∂1 v2 dt ∂1 v3 dt
dl3 ∂2 v1 dt 1 + ∂2 v2 dt ∂2 v3 dt = dl3 + dl3 dt(∂1 v1 + ∂2 v2 + ∂3 v3 )
∂3 v1 dt ∂3 v2 dt 1 + ∂3 v3 dt
12
Donde hemos despreciado términos pequeños de orden dt2 y dt3 . La ecuación
anterior revela que la divergencia de la velocidad ∇ · ⃗v = ∂1 v1 + ∂2 v2 + ∂3 v3
representa fı́sicamente la velocidad a la que varı́a el volumen de un elemento
fluido por unidad de volumen. A dicha cantidad, que coincide con la traza del
tensor gradiente de velocidades y del tensor de deformación, se la denomina
velocidad de dilatación cubica unitaria. En el caso de un lı́quido perfecto, para
el que el volumen de un elemento fluido no varia, es decir de densidad constante,
se debe por tanto cumplir que:
∇ · ⃗v = 0
13
5 Pregunta 5: Obtener el teorema de las fuerzas
vivas.
A partid de la segunda Ley de Navier Stokes,
D⃗v ∂⃗v ′
ρ =ρ + ∇(v /2) − ⃗v ∧ (∇ ∧ ⃗v = −∇p + ∇ · τ + ρf⃗m
2
Dt ∂t
Está ecuación, que deriva de la segunda ley de Newton, indica que la variación
de v 2 /2 es consecuencia directa del trabajo de las fuerzas másicas ρ⃗v · f⃗m y de
una parte del trabajo de las fuerzas de superficie, ⃗v · (∇ · τ ). Esta segunda
contribución resulta del producto escalar de la resultante de las fuerzas de su-
perficie por unidad de volumen, ∇ · τ , por la velocidad del centro de masas de
la partı́cula fluida, ⃗v , y corresponderı́a, por tanto, al trabajo total de las fuerzas
de superficie (por unidad de volumen y por unidad de tiempo) si la partı́cula
fluida no sufriera deformación.
14
6 Pregunta 6: Navier - Poisson demostración.
Enuncie la ley de Navier Poisson e indique las hipótesis que debe
cumplir un fluido para que se verifique dicha ley y las implicaciones
de cada una de esas hipótesis. A partid de dichas hipotesis, demuestre
la ley de Navier Poisson.
τi′ = (α − λ)γi + λ∇ · ⃗v
15
Con µ = (α −λ)/2 y λ primer y segundo coeficiente de viscosidad. Reescribimos
el tensor como:
′
τ = 2µT d + λ(∇ · ⃗v )I
Que en cualquier SR podemos escribir:
′ 2
τ = 2µT d + (µv − µ)(∇ · ⃗v )I
3
Que es la ley de Navier-Poisson.
Ya que la velocidad de dilatación cúbica unitaria ∇ · ⃗v es invariante y podemos
reemplazar los ejes principales por sistemas de referencia ortogonales.
16
7 Pregunta 7: Ley de Fourier.
Transporte de calor por conducción. Ley de Fourier.
De forma similar a como ocurre con los gradientes de velocidad que dan
lugar a fuerzas de superficie, vamos a ver como un gradiente de temperatura da
lugar a un flujo de calor por conducción.
qn = q1 n1 + q2 + n2 + q3 n3 (2)
qn = ⃗q · ⃗n (3)
Donde ⃗q es el vector flujo de calor ⃗q = (q1 , q2 , q3 ) y es función de ⃗x y t.
Ley de Fourier
17
Esta ley establece que el flujo de calor es proporcional al gradiente de tem-
peratura (hacia temperaturas decrecientes). Para los gases nos basamos en la
teorı́a cinética y para los liquidos en evidencia experimental.
⃗q = −k∇T (6)
Siendo k el coeficiente de conductividad térmica. k depende del estado ter-
modinámico del fluido, con el ETL podemos definirlo con dos variables, como p
y T . Podemos despreciar la presión y tendrı́amos que:
18
8 Pregunta 8: Tensor gradiente de velocidades.
Indique porqué el movimiento en el entorno de un punto fluido está
caracterizado por el tensor gradiente de velocidades. Diga en que dos
tensores se descompone este y lo que cada uno de ellos representan.
Por facilidad, utilice coordenadas cartesianas.
En ese instante, las velocidades de los extremos del elemento fluido diferirán
una cantidad d⃗v = d⃗x · ∇⃗v , con ∇⃗v ≡ tensor gradiente de velocidades.
Tras un tiempo dt los extremos pasan a valor ⃗x +⃗v dt y (⃗x + d⃗x) + (⃗v + d⃗v )dt.
Aquı́ tenemos una traslación ⃗v dt y una deformación del elemento fluido d⃗x +d⃗v dt
19
Como tenemos que d⃗v = d⃗x · ∇⃗v , tenemos:
Este tensor es el que contiene las propiedades de deformación propias del fluido,
pues si fuese igual a cero tendrı́amos un sólido rı́gido (rotación más traslación).
Esto nos permite escribir la velocidad de deformación de la siguiente forma:
d⃗vd = T d · ⃗ndS
con d⃗x = ⃗ndS. Como d⃗vd no suele tener la dirección de ⃗n, d⃗x sufre un giro
además de la dilatación lineal.
20
9 Pregunta 9: Descripciones eulerianas, lagrangianas
y derivada material
Descripción euleriana y lagrangiana.
Dando por admitida la HMC, para describir el movimiento del medio con-
tinuo fluido podemos usar 2 alternativas:
Es por ello por lo que usamos la descripción euleriana, la cual describe cada
magnitud fluidomecánica en el punto ⃗x y un instante t, es decir en cada instante
se da la distribución espacial en el interior del campo fluido f (⃗x, t).
Derivada material:
∂ρ
∂t X
Aquı́ solo tenemos la descripción espacial, no nos sirve. Tomemos una coorde-
nada x de la forma x(X, t) con la regla de la cadena queda:
∂ρ(⃗x, t) ∂ρ(x(X, t), t) ∂ρ(⃗x, t) ∂ρ(⃗x, t) ∂⃗x(X, t)
= = + · =
∂t x ∂t x ∂t x ∂⃗x ∂t X
21
∂ρ(⃗x, t) ∂ρ(⃗x, t) ∂ρ(⃗x, t)
= + · ⃗v (⃗x, t)
∂t x ∂t x ∂⃗x
Con lo que definimos ası́ el operador derivada material:
D() ∂()
= + ⃗v · ∇()
Dt ∂t
Donde el primer termino es la variación temporal y el segundo la variación por
el movimiento de la partı́cula fluida.
22
10 Pregunta 10: Balance de Ec y energı́a in-
terna.
Indique que representan cada una de las siguientes ecuaciones, co-
mente lo que representan cada uno de los términos que aparecen en
ellas y la influencia que estos términos tienen en la magnitud que la
ecuación caracteriza. Solución:
23
sobre la superficie fluida, dan lugar a un aumento de la energı́a cinética del
volumen fluido cuando no hay deformación (teorema de las fuerzas vivas).
24
11 Pregunta 11: Hipótesis del medio continuo.
Explique en qué consisten las hipótesis de medio continuo (HMC)
(partı́cula fluida) y de equilibrio termodinámico local (ETL), indi-
cando los parámetros que las caracterizan. ¿Por qué la segunda es
mucho más restrictiva que la primera? (Nota: no olvide incluir los
razonamientos relativos a la evolución temporal.)
25
En el caso de la mecánica de fluidos, el sistema no está en equilibrio, pueden
variar en el tiempo las variables y no tienen por qué ser uniformes.
Número de Kundsen.
De esta forma, cada molécula sufrirá muchas colisiones hasta alcanzar una
región de equilibrio estable apreciable macroscópicamente. Es como si en cada
punto del fluido estuviera casi en equilibrio correspondiente a las variables
termodinámicas locales. Para darse los choques, el volumen barrido por una
molécula en su movimiento ( d20 λ) debe ser igual al volumen que le corresponde
a cada molécula (d30 ). Por tanto, λ/d ≈ (d/d0 )2
Por último, indicar que el ETL, λ/L <<< 1 es más restrictivo que la HMC,
d/L <<< 1 ya que λ <<< d <<< L
26
12 Pregunta 12: Principio de la cantidad de
movimiento.
Sea la ecuación:
D⃗v ′
ρ = −∇p + ∇ · τ + ρf⃗m
Dt
a) Diga de qué ecuación se trata y lo que representa. Indique el nom-
bre y el significado fı́sico de cada uno de los términos que aparecen en
ella. ¿Por qué se debe usar la derivada material en esta ecuación? (en
algunas convocatorias añade lo siguiente) Escriba la ley de Navier Poisson
′
para el cálculo de τ , indicando el nombre y el significado fı́sico de
las magnitudes y parámetros que aparecen en ella. Indique las condi-
ciones que debe cumplir para poder aplicar dicha ley.
′ 2
τ = 2µT d + (µv − µ)(∇ · ⃗v )I (9)
3
′
• τ tensor de esfuerzos viscosos
• µ primer coeficiente de viscosidad; µ = 12 (α − λ)
27
• (µv − 23 µ)(∇ ·⃗v )I Termino que, junto con el primer sumando, nos deja con
′
la relación τ 1 = αγ1 + λγ2 + λγ3
La condición para poder aplicar esta ley es que el fluido es isótropo, newtoniano,
T d ̸= 0, considerar la dirección principal de T d y ∇ · ⃗v invariable, por lo que se
′
puede representar para direcciones ortogonales a la principal de τ
28
Por otro lado, al tener el gradiente de la presión, si añadimos un valor con-
stante a la presión nos generara variación, simplificando ası́ el resultado.
v2
c = p + ρU + ρ
2
Por lo que en fluidoestática, donde tomaremos sistemas de referencia tales que
la velocidad del fluido sea 0:
c = p + ρU
29
13 Pregunta 13: Cambios de entropı́a.
Los cambios de entropia en un volumen fluido Vf (t) se obtienen a
partir de la siguiente ecuación.
"Z # Z Z Z Z
d
ρsdV = (ϕv /T )dV − [(⃗q·∇T )/T 2 ]dV − (⃗q/T )·⃗ndσ+ (Qr /T +Qq /T )dV
dt Vf (t) Vf (t) Vf (t) Σf (t) Vf (t)
30
14 Pregunta 14: Ecuación de continuidad
a) Encuentre la ecuación de continuidad a partir del balance de masa
en un volumen fluido.
Vamos a buscar una ecuación que nos muestra que la masa no varı́a en un
volumen fluido con el tiempo.
Los subı́ndices c son de control la zona a estudiar y estos pueden variar con el
tiempo.
31
∇ · ⃗v = 0 Fluido incompresible
∇ · (ρ⃗v ) = 0 Movimiento estacionario de gases.
∇ · (ρ⃗v ) = ρ∇ · ⃗v = 0 = ∇ · ⃗v = 0 (16)
32
15 Pregunta 15: Ec de la energı́a interna.
Sea la ecuación
De ′
ρ = −p∇ · ⃗v + τ : ∇⃗v − ∇ · ⃗q + Qr + Qq
Dt
diga de qué ecuación se trata y lo que representa. Indique el nombre
y el significado fı́sico de cada uno de los términos que aparecen en
ella. ¿Por qué debe usarse la derivada material en esta ecuación?
33
16 Pregunta 16: Fluidos newtonianos.
¿Que son fluidos newtonianos?
′ 2
τ = 2µT d + (µv − µ)(∇ · ⃗v )
3
Qué son los coeficientes µ y µv
Ser isótropo. Para los fluidos newtonianos (la mayorı́a) se observa que los
esfuerzos viscosos son proporcionales a las velocidades de deformación.
′
τij = αijk γkl
34