FLUIDO

Con todos los supuestos dados anteriormente, el problema del fluido es el siguiente:
( ∂v
ρ ∂xii = 0
∂p ∂ 2 vi (1)
ρ ∂v ∂vi
∂t + ρvj ∂xj = − ∂xi + µ ∂xj ∂xj
i

Observe que para la primera ecuación de campo incompresible (o de divergencia nula) se puede reescribir de
la siguiente manera:
ρdiv(vi ) = 0
por otro lado, en la segunda ecuación tenemos:

∂vi ∂vi ∂p ∂ 2 vi
ρ + ρvj =− +µ
∂t ∂xj ∂xi ∂xj ∂xj

con v = (v1 , v2 ) y x = (x1 , x2 ). Comenzaremos a desarrollar para i = 1, 2 y j = 1, 2.

Para i = 1:
∂v1 ∂v1 ∂p ∂ 2 v1
ρ + ρvj =− +µ
∂t ∂xj ∂x1 ∂xj ∂xj
Para j = 1, 2:
∂v1 ∂v1 ∂v1 ∂p ∂ 2 v1 ∂ 2 v1
ρ + ρ[v1 + v2 ]=− + µ[ 2 + ]
∂t ∂x1 ∂x2 ∂x1 ∂x1 ∂x22
| {z }
laplacian

equivalentemente
∂v1 ∂v1 ∂v1 ∂p
ρ + ρ[(v1 , v2 ) · ( , )] = − + µ∆v1 (x1 , x2 )
∂t ∂x1 ∂x2 ∂x1
∂v1 ∂p
ρ + ρv · ∇v1 = − + µ · ∆v1
∂t ∂x1
Análogamente para i = 2, de manera que las ecuaciones Navier-Stokes bidimensional se expresan:
(
∂p
ρ ∂v
∂t + ρv · ∇v1 = − ∂x1 + µ · ∆v1
1

∂v2 ∂p
ρ ∂t + ρv · ∇v2 = − ∂x 2
+ µ · ∆v2

en forma matricial:  ∂v1     ∂p   

∂t ∇v1 ∂x1 ∆v1
ρ ∂v2
+ ρv =− ∂p
+µ
∂t
∇v2 ∂x2
∆v2
| {z }
gradp

Entonces la ecuación (1) se escribirá:


ρdiv(v) = 0
(2)
ρ ∂v
∂t + ρv · ∇v = −grad p + µ · ∆v

Observación. Considerando la siguiente ley constitutiva:

σf = −ρI + µf (∇v)

donde I es el tensor identidad, σf es el tensor de Cauchy, y µf = ρv es la viscosidad dinámica, mientras que

v es la viscosidad cinemática. Entonces, se cumple que

div(σf ) = div(−ρI + µf (∇v))

= div(−ρ · I) + div(µf (∇v))
= −[(grad p) · I + ρ ·   + µf · div(∇v)
div(I)]
= −grad p + µf · ∆v

1
luego, la ecuación (2) se escribirá:

ρ ∂v

∂t + ρ(v · ∇)v − div(σf ) = 0
ρdiv(v) = 0

donde
σf = −ρI + µf (∇v)
De manera particular, la formulación nos queda:
 ∂v
f
 ρ. + ρ(vO)vf − div(σt ) = 0 en Ωf (t)
∂t



ρ.div(vf ) = 0 en Ωf (t)


vf (x, 0) = vf0 (x) en Ωf (t) (3)

vf (x, t) = vfD (x, t) en ΓD
f (x, t)




ΓN

σf .~nf = gf (x, t) en f (t)

Donde el tensor de Cauchy σf satisface la ley constitutiva:

σ = −pI + µ(∇vf + ∇vfT )

2
 ESTRUCTURA
( s
∂σij
∂xj = 0 (4)
s E vE
σij = 1+v εij + (1+v)(1−2v) εkk δij
∂di ∂dj
donde εij = 21 ( ∂x j
+ ∂xi ) (5)

s
Aquı́, σij es el tensor de tensiones, E es el módulo de Young’s, v es el radio de Poisson, εij es el tensor
de deformación, δij es la delta de Kronecker y di es el desplazamiento en una dirección dada.
Observación. εkk = ε11 + ε22 = tr(ε)
En la primera ecuación de (4), tenemos que
s
∂σij
=0
∂xj

haciendo variar j = 1, 2.
s s
∂σi1 ∂σi2
+ =0
∂x1 ∂x2
ahora, variando i = 1, 2 ( s s
∂σ11 ∂σ12
∂xs1 + ∂xs2 =0
∂σ21 ∂σ22
∂x1 + ∂x2 =0
en forma matricial: s s
 ∂σ11 ∂σ12
+
  
∂x1 ∂x2 0
s
∂σ21 s
∂σ22 =
+ 0
∂x1 ∂x2

entonces que la primera ecuación de (4) se puede reescribir de la siguiente manera:

s
div(σij )=0 (5)

Por otro lado, en la segunda ecuación tenemos:

s E vE
σij = εij + εkk δij
1+v (1 + v)(1 − 2v)
teniendo en cuenta que estamos trabajando en el plano bidimensional, este tensor va a tener cuatro compo-
nentes, es decir:

s E vE
σ11 = ε11 + (ε11 + ε22 )δ11
1+v (1 + v)(1 − 2v)
s E vE ((((
σ12 = ε12 + ( ( (ε
( (+
11 ε22 )δ12
1+v (1(
( +(v)(1
(( − 2v)
s E vE ((((
σ21 = ε21 + (
( − 2v)( (ε
(11 (+ ε22 )δ21
1+v (1(
( +(v)(1(

s E vE
σ22 = ε22 + (ε11 + ε22 )δ22
1+v (1 + v)(1 − 2v)
s
*Nótese que δij es la delta de Kroenecker. Teniendo en cuenta la Observación, el tensor σij matricialmente
se expresará de la forma:
" #
E vE E
s 1+v ε11 + (1+v)(1−2v) (tr(ε)) 1+v ε12
σij = E E vE
1+v ε21 1+v ε22 + (1+v)(1−2v) (tr(ε))

equivalentemente
" #
E E vE
(1+v)(1−2v) (tr(ε)) 0
 
s
σij = 1+v ε11 1+v ε12 +
E E vE
1+v ε21 1+v ε22 0 (1+v)(1−2v) (tr(ε))

3
    
s E ε11 ε12 vE tr(ε) 0
σij = +
1 + v ε21 ε22 (1 + v)(1 − 2v) 0 tr(ε)
   
s E ε11 ε12 vE 1 0
σij = + tr(ε)
1 + v ε21 ε22 (1 + v)(1 − 2v) 0 1
de manera compacta el tensor para la estructura satisface la ecuación constitutiva
E vE
σs = ε+ tr(ε)δ (6)
1+v (1 + v)(1 − 2v)

Por lo tanto, considerando las ecuaciones(5) y (6), la ecuación (4) se escribirá como:
 s
div(σij )=0
s E vE (7)
σ = 1+v ε+ (1+v)(1−2v) tr(ε)δ

Acontinuación estudiaremos ahora el tensor de deformación εij

1 ∂di ∂dj
εij = ( + )
2 ∂xj ∂xi

en el plano de dos dimensiones, este tensor tendrá cuatro componentes:

 
ε11 ε12
ε=
ε21 ε22
" #
1 ∂d1 ∂d1 1 ∂d1 ∂d2
( + ) ( + )
= 21 ∂x 1
∂d2
∂x1
∂d1
2 ∂x2
1 ∂d2
∂x1
∂d2
2 ( ∂x1 + ∂x2 ) 2 ( ∂x2 + ∂x2 )
" # " #
∂d1 ∂d1 ∂d1 ∂d2
1 ∂x ∂x2 ∂x1 ∂x1
= ( ∂d2 ∂d2 + ∂d1 ∂d2 )
1
2 ∂x1 ∂x2 ∂x2 ∂x2
1
ε= (∇d + ∇dT )
2

Observación. Recuerde que se dijo que el tensor de tensiones es formulado en el marco Euleriano,
entonces en el marco lagrangiano, ∇d = ∇dˆ · F −1 y σ s se lee como P = Jˆσˆs F̂ −T , entonces, la última
ecuación se escribirá
1
ε̂ = (∇dˆ · F −1 + (∇dˆ · F −1 )T )
2
1
ε̂ = (∇dˆ · F −1 + F −T ∇dˆT )
2
el tensor para la estructura en el marco lagrangiano se escribirá de la siguiente manera: