1.

problema autosimilar (Problema de conducción transitoria unidimensional, en forma

adimensional)

∂2 T 1 ∂ T
2
=
∂x α ∂t

Condiciones iniciales:
T ( 0 , t )=T S

T ( x →∞ , t ) =T o

T ( x , 0 )=T o

∆ T =T −T o

∆ T S =T S−T o

Desarrollo de la adimensionalización de la ecuación:

∆ T o ∂2 T ¿ 1 ∆ To ∂ T
¿

2 ¿2
= ¿ ¿
lo ∂ x αo α to ∂ t

α o t o ∂2 T ¿ 1 ∂ T ¿
2 ¿2
= ¿ ¿
lo ∂ x α ∂t

α o to
Fo= 2 (Número de Fourier)
lo

∂ 2 T ¿ 1 ∂T ¿
Fo ¿2
= ¿ ¿
∂x α ∂t

A continuación de hallan los números adimensionales (𝜋), mediante el método PI

∆ T (x ,t , α , ∆ T S )
 [ ]
∆ T x t α ∆T S
L
0 10 2 0
t
0 0 1 −1 0
T
1 00 0 1

Para la cantidad de números adimensionales (𝜋):

n=i−k =5−3=2

Se considera la matriz A para hallar si el Det(A) es L.I.

[ ]
t α ∆TS
L
0 2 0
A= t
1 −1 0
T
0 0 1

Det ( A )=0 ( 0 )−2 ( 1−0 ) + 0(0)

Det ( A )=−2≠ 0 → L . I .

π=t n α n ∆ T Sn La t b T c
1 2 3

¿
¿

{
n1 +b−n2=0
2 n2 +a=0
n3 + c=0

−a
n1 =n2−b → n1= −b
2
−a
n2 =
2
n 3=−c

1 n1 n2 n3 0 0 −1
π 1=∆T t α ∆ T S =∆ T t α ∆ T S

∆T
π 1= =1
∆TS
−1 −1
1 n1 n2 n3 2 2 0
π 2=x t α ∆ T S =x t α ∆TS
 x
π 2=
√t √α
x
π 2=
√t α

De la ecuación adimensional hallada anteriormente:

∂ 2 T ¿ 1 ∂T ¿
Fo = ¿ ¿
∂ x ¿2 α ∂ t
Se hallan las derivadas parciales de esta ecuación con respecto al número adimensional (
π 2).

( )
2 ¿
∂ T
Para la doble derivada parcial ¿2 se tiene:
∂x

¿
x =π 2

∂T ∂ T ∂ π 2
¿ ¿

¿ = ¿
∂x ∂ π2 ∂ x

∂ ∂T ∂π2
[ ]
2 ¿ ¿
∂T
=
∂x
¿2
∂ x¿ ∂ π2 ∂ x¿

( )
2
∂2 T ¿ ∂ ∂T ¿ ∂ π 2 ∂T ¿ ∂ π 2
= ¿ +
∂x
¿2
∂ x ∂ π 2 ∂ x ¿ ∂ π 2 ∂ x ¿2

∂ π (∂ π ) ∂ x ∂ π ∂x
2
∂ ∂T ∂π ∂T ∂ π
2 ¿ ¿ ¿
∂T 2 2
=¿2
+ ¿ ¿2
∂x 2 2 2

( )
2 2
∂ T ∂ T ∂π
2 ¿
∂T ∂ π
2 ¿
2
¿
2
= ¿2
+ 2 ¿ ¿2
∂x ∂π ∂x ∂π ∂ x
2 2
 x¿
π 2=
√ t¿ α ¿
∂ π2 1
¿=
∂ x √ t¿ α¿
2

( √ t1α )=0
∂ π2 ∂
¿2
= ¿ ¿ ¿
∂x ∂x

∂2 T ¿ ∂2 T ¿
=
∂ x ¿2 ∂ π 22 ( +
)
1 2 ∂T¿
√ t ¿ α¿ ∂ π 2
( 0)

2 ¿ 2 ¿
∂T ∂ T 1
= ¿ ¿
∂ x ¿2 ∂ π 22 t α

Para la derivada parcial

∂ T¿
∂t
¿ se tiene: ( )
∂T ¿ ∂ T ¿ ∂ π 2
=
∂ t ¿ ∂ π2 ∂ t ¿

∂ π 2 ∂ x ¿ t ¿−0,5
∂t
¿=
∂t
¿
√α¿ [ ]
∂ π 2 −x ¿ t ¿−0,5−1
¿=
∂t 2 √α ¿
∂ π2 −x
¿
=
∂ t ¿ 2t ¿ √ t ¿ α ¿

¿ ¿ ¿
∂T ∂ T −x
¿=
∂t ∂ π 2 2 t¿ √ t ¿ α¿

De esta forma se tiene que:

∂2 T ¿ 1 1 ∂ T ¿ −x¿
Fo 2 ¿ ¿
= ¿
∂ π 2 t α α ∂ π 2 2t √ t α
¿ ¿ ¿
 2 ¿ ¿
∂ T −1 ∂T
Fo = π2
∂ π 22 2 ∂ π 2

2 ¿ ¿
d T −1 d T
→ Fo = π
d π2
2
2 2 d π2
¿
dT
→ u=
d π2
2 ¿
d T
→ du= 2
d π2
d π2

du d 2 T ¿
→ =
d π 2 d π 22

du −1
Fo = π u
d π2 2 2

du −1
→ Fo = π dπ
u 2 2 2
du −1
→ Fo ∫ = ∫ π2 d π2
u 2

[ ]
2
−1 π 2
→ FoLn(u)= +c
2 2
−1 2
→ FoLn(u)= π +c
4 2

→ u=e
[
1 −1 2
Fo 4
π2 + c
]

→
dT
=e Fo
[¿ 1 −1 2
4
π 2 +c
]
d π2

→∫ d T =∫ A e
¿ Fo 4[
1 −1 2
π2
]d π
2

→ T ¿ =∫ A e Fo 4[
1 −1 2
π2
]d π
2

Donde (A) es una constante.

De las condiciones iniciales:
T ( 0 , t )=T S (Superficie) para t > 0

T ( x →∞ , t ) =T o

T ( x , 0 )=T o (Inicial)

( )
T ( x , t ) −T s x
→ =erf
T o−T s 2 √t α

La función error está dada por:

x
2 √tα

( 2√xtα )= √2π ∫ e
2
−η
erf dη
❑

La distribución de temperatura se da:

x
2 √ tα
T ( x , t )−T s 2
∫
2

= e−η d η
T o−T s √π ❑

Para el flujo de calor en una posición (x) se puede dar así:

∂T
q x =−K A Á rea
∂x

Adoptando esta solución a nuestro problema, nos queda:

T ( x , t )−T s 2
= A
[ ]
2 √ tα

∫ e dπ
−π 22
4 Fo
T o−T s √π ❑
2
Aplicando la derivada parcial:

[ ]
2
−π 2
∂T
∂x
=( T o−T s )
2
√π
Ae
4 Fo ∂
(
x
∂ x 2 √ tα )
[ ]
2
−π 2

→
∂T
∂x
=( T o−T s )
2A
√π
e
4 Fo
( 2 √1tα )
[ ]
2
−π 2
4 Fo
∂T Ae
→ =( T o−T s )
∂x √tα π
Donde (A) es una constante.

2. Adimensionalización de las ecuaciones de continuidad y de Navier-Stokes en

coordenadas cartesianas.

Componente X de la ecuación de Navier-Stokes de flujo incompresible:

( ) ( )
2 2 2
∂u ∂u ∂u ∂ u −∂ P ∂ u ∂ u ∂ u
⍴ +u +v +w = +⍴ g x + µ + 2+ 2
∂t ∂x ∂y ∂z ∂x 2
∂x ∂ y ∂z
Las variables (u , v y w ) son las componentes de la velocidad en las coordenadas ( x , y , z)
respectivamente.
Desarrollo adimensional:

∂u
∂t
+u
∂u
∂x
+v
∂u
∂y
+w = + +
(
∂ u −1 ∂ P µ ∂ 2 u ∂2 u ∂2 u
+
∂ z ⍴ ∂ x ⍴ ∂ x 2 ∂ y 2 ∂ z2
+ gx
)
Los parámetros (⍴ , µ y g x ) no se adimencionalizan ya que son parámetros propios del
fluido.

+ u + v
t o ∂t ¿ l o ∂ x ¿ l o ∂ y ¿ l o
+ w = + + (
v o ∂ u¿ v o2 ¿ ∂u ¿ v o2 ¿ ∂ u¿ v o2 ¿ ∂u ¿ −∆ Po ∂ P¿ µ v o ∂2 u ¿ ∂2 u¿ ∂2 u ¿
+
∂ z¿ ⍴ l o ∂ x ¿ ⍴ l o2 ∂ x ¿2 ∂ y ¿2 ∂ z ¿2
+ gx )
Similitud:
lo
t o= , ∆ P o=⍴ v o2
vo

[ ] ( )
2 2
v o ∂ u¿ ¿ ∂ u¿ ¿ ∂u ¿ ¿ ∂u
¿
−⍴ v o ∂ P¿ µ v o ∂2 u ¿ ∂2 u¿ ∂2 u¿
→ +u ¿+v ¿+w ¿ = + + + + gx
l o ∂t ¿ ∂x ∂y ∂z ⍴ l o ∂ x ¿ ⍴ l o2 ∂ x ¿2 ∂ y ¿2 ∂ z ¿2

( )
g x lo
¿ ¿ ¿ ¿ ¿ 2 ¿ 2 ¿ 2 ¿
∂u ¿ ∂u ¿ ∂u ¿ ∂u −∂ P µ ∂u ∂ u ∂u
→ ¿ +u ¿ +v ¿ +w ¿= ¿ + + + +
∂t ∂x ∂y ∂z ∂x ⍴ v o l o ∂ x ¿2 ∂ y ¿2 ∂ z ¿2 vo
2

⍴ v o lo 1 µ
Re = → = (Número de Reynolds)
µ R e ⍴ v o lo

vo vo2
2
1 g x lo
F r= → Fr = → = 2 (Número de Froude)
√ g x lo g x lo Fr
2
vo

Solución adimensional:

∂u ¿ ¿ ∂ u¿ ¿ ∂ u¿
( )
¿
¿ ∂u −∂ P¿ 1 ∂2 u ¿ ∂ 2 u¿ ∂2 u ¿ 1
¿ + u ¿ + v ¿ +w ¿ = ¿ + ¿2
+ ¿2 + ¿2 + 2
∂t ∂x ∂y ∂z ∂x Re ∂ x ∂y ∂z Fr

Componente Y de la ecuación de Navier-Stokes de flujo incompresible:

 ⍴ ( ∂v
∂t
+u
∂v
∂x
+v
∂v
∂y
+w
∂ v −∂ P
∂z
=
∂y )
+ ⍴ gy + µ + (
∂2 v ∂ 2 v ∂ 2 v
+
∂ x2 ∂ y 2 ∂ z2 )
Las variables (u , v y w ) son las componentes de la velocidad en las coordenadas ( x , y , z)
respectivamente.
Desarrollo adimensional:

∂v
∂t
+u
∂v
∂x
+v
∂v
∂y
+w
∂ v −1 ∂ P µ ∂2 v ∂2 v ∂2 v
= + + +
∂ z ⍴ ∂ y ⍴ ∂ x2 ∂ y2 ∂ z2 (
+gy
)
Los parámetros (⍴ , µ y g y ) no se adimencionalizan ya que son parámetros propios del
fluido.

( )
2 2 2
v o ∂ v ¿ v o ¿ ∂ v ¿ v o ¿ ∂ v ¿ v o ¿ ∂ v ¿ −∆ P o ∂ P¿ µ v o ∂2 v ¿ ∂2 v ¿ ∂ 2 v ¿
+ u ¿+ v + w ¿= + + + + gy
t o ∂t ¿ l o ∂x lo ∂ y ¿ lo ∂z ⍴ l o ∂ y ¿ ⍴ l o2 ∂ x ¿2 ∂ y ¿2 ∂ z ¿2

Similitud:
lo
t o= , ∆ P o=⍴ v o2
vo

[ ] ( )
2 2
v o ∂ v¿ ¿ ∂ v ¿ ¿ ∂ v ¿ ¿ ∂v
¿
−⍴ v o ∂ P¿ µ v o ∂2 v ¿ ∂2 v ¿ ∂2 v ¿
→ +u +v +w = + + + + gy
l o ∂t ¿ ∂x
¿
∂y
¿
∂z
¿
⍴l o ∂ y ¿ ⍴ l o2 ∂ x ¿2 ∂ y ¿2 ∂ z ¿2

( g y lo
)
¿ ¿ ¿ ¿ ¿ 2 ¿ 2 ¿ 2 ¿
∂v ¿ ∂v ¿∂v ¿ ∂v −∂ P µ ∂ v ∂ v ∂ v
→ ¿ +u ¿ +v ¿ +w ¿= ¿ + ¿2
+ ¿2 + ¿2 + 2
∂t ∂x ∂y ∂z ∂y ⍴ v o lo ∂ x ∂y ∂z vo

⍴ v o lo 1 µ
Re = → = (Número de Reynolds)
µ R e ⍴ v o lo

vo vo2 1 g y lo
F r= 2
→ Fr = → = (Número de Froude)
√ g y lo gylo F r2 v o2

Solución adimensional:

( )
¿ ¿ ¿ ¿ ¿ 2 ¿ 2 ¿ 2 ¿
∂v ¿ ∂v ¿∂v ¿ ∂v −∂ P 1 ∂ v ∂ v ∂ v 1
¿ +u ¿ +v ¿ +w ¿= ¿ + ¿2
+ ¿2 + ¿2 + 2
∂t ∂x ∂y ∂z ∂y Re ∂ x ∂y ∂z Fr
Componente Z de la ecuación de Navier-Stokes de flujo incompresible:

( ) ( )
2 2 2
∂w ∂w ∂w ∂ w −∂ P ∂ w ∂w ∂ w
⍴ +u +v +w = +⍴ g z + µ + +
∂t ∂x ∂y ∂z ∂z ∂ x 2 ∂ y 2 ∂ z2

Las variables (u , v y w ) son las componentes de la velocidad en las coordenadas ( x , y , z)

respectivamente.
Desarrollo adimensional:

( )
2 2 2
∂w ∂w ∂w ∂ w −1 ∂ P µ ∂ w ∂ w ∂ w
+u +v +w = + + + + gz
∂t ∂x ∂y ∂ z ⍴ ∂ z ⍴ ∂ x2 ∂ y2 ∂ z2

Los parámetros (⍴ , µ y g z ) no se adimencionalizan ya que son parámetros propios del

fluido.

( )
2 2 2
v o ∂ w ¿ v o ¿ ∂ w ¿ v o ¿ ∂ w¿ v o ¿ ∂ w ¿ −∆ Po ∂ P ¿ µ v o ∂2 w ¿ ∂2 w ¿ ∂2 w ¿
+ u ¿+ v ¿+ w ¿ = + + + + gz
t o ∂ t ¿ lo ∂x lo ∂y lo ∂z ⍴l o ∂ z ¿ ⍴ l o2 ∂ x¿2 ∂ y ¿2 ∂ z ¿2

Similitud:
lo
t o= , ∆ P o=⍴ v o2
vo

[ ] ( )
2 2
v o ∂ w ¿ ¿ ∂ w¿ ¿ ∂ w ¿ ¿ ∂ w ¿ −⍴ v o ∂ P ¿ µ v o ∂ 2 w ¿ ∂ 2 w ¿ ∂ 2 w ¿
→ +u ¿ +v ¿ +w ¿ = + + + + gz
lo ∂ t ¿ ∂x ∂y ∂z ⍴ l o ∂ z ¿ ⍴ l o2 ∂ x¿2 ∂ y ¿2 ∂ z ¿2

( g z lo
)
¿ ¿ ¿ ¿ ¿ 2 ¿ 2 ¿ 2 ¿
∂w ¿ ∂w ¿ ∂w ¿ ∂w −∂ P µ ∂ w ∂ w ∂ w
→ ¿ +u ¿ + v ¿ + w ¿ = ¿ + ¿2
+ ¿2
+ ¿2
+ 2
∂t ∂x ∂y ∂z ∂z ⍴ vo l o ∂ x ∂y ∂z vo

⍴ v o lo 1 µ
Re = → = (Número de Reynolds)
µ R e ⍴ v o lo

vo v o2 1 gl
F r= → Fr =
2
→ = z 2o (Número de Froude)
√ g z lo g z lo Fr
2
vo
Solución adimensional:

( )
¿ ¿ ¿ ¿ ¿ 2 ¿ 2 ¿ 2 ¿
∂w ¿ ∂w ¿ ∂w ¿ ∂w −∂ P 1 ∂ w ∂ w ∂ w 1
¿ +u ¿ +v ¿ +w ¿= ¿ + + + +
∂t ∂x ∂y ∂z ∂z Re ∂ x ¿2 ∂ y ¿2 ∂ z ¿2 F r2