“Año del diálogo y la Reconciliación”

UNIVERSIDAD NACIONAL MAYOR DE

SAN MARCOS
(Universidad del Perú, Decana de América)

FACULTAD DE CIENCIAS FISICAS

E.A.P. DE INGENIERIA MECANICA DE FLUIDOS

“TAREAS 1, 2 y 3 DE FLUJO NO PERMANENTE”

CURSO: Flujo No Permanente

PROFESOR: William Chauca Nolasco

Alumno:
-Pariona Tenorio, Gabriel Masiño
Código:
-14130168
CICLO: Ciclo 8ctavo (Verano).

Ciudad Universitaria, 29 de enero del 2018

Tarea 1
Integrar la ecuación

De la ecuación de diferencias parcial para el caso de resistencia turbulenta

1 ∂P ∂z ∂ v ∂ v f v2
+g +v + + =0
ρ ∂s ∂s ∂s ∂t 2 D

∂V
Integrando la ecuación entre dos puntos para P=0, =0
∂s

Por lo tanto:

∂v fL v 2
g ( z 2−z 1) + L+ =0
∂t 2D

Donde: z 2−z 1=2 z

∂v fL v 2
2 gz + L+ =0
∂t 2D

Podemos expresar reordenando (÷ L)

dv fv |v| 2 gz
+ + =0
dt 2 D L

dz
Donde v=
dt

d dz
( ) | |( )
+
f dz dz 2 gz
dt dt 2 D dt dt
+
L
=0

2
dz
2
+
d t 2 D dt dt | |( )
f dz dz 2 gz
+
L
=0

Para la ecuación, se puede Integrar una vez Haciendo

dz
P=
dt

De donde
2
d z dP dP
= =
d t 2 dt dz ( )( dzdt )=P ( dPdz )
Reemplazando en la ecuación restringiendo el movimiento a la dirección (-z).
2
dP P 2 gz
P ( )
dz
−f +
2D L
=0

páá g. 2
Integrando la ecuación
−fz
Multiplicando ambos extremos por el factor integrante e D

−fz −fz
P2 2 gz
e D
× P (( )
dP
dz
−f +
2D L
=0 ×e ) D

−fz −fz −fz

dP P2 2 gz
e D
.P ( ) dz
−f
2D
.e D
+
L
.e D
=0

−fz −fz −fz

D P2 D 2 gz D
e . PdP−f .e dz + .e dz=0
2D L
−fz −fz −fz
P2 2g
e D
. ∫ PdP−f .∫ e D
dz + .∫ z e D
dz=0
2D L
−fz −fz −fz
P2 2g
D
e⏟ . ∫ PdP Integral
− 1f . ∫ e D dz + . ∫ z e D dz =0
⏟ ⏟
2 D L
Integral 2 Integral 3

Integrando la Integral 1
−fz −fz
P2
e D
. ∫ PdP=e D +K 1
2

Integrando la Integral 2
−fz −fz −fz
P2 P 2 −D P2
f
2D ∫
. e D
dz=f
2D
.
f
e ( D
) + K 2 =−e D
2
+K2

Integrando la Integral 3, Integrando por el método por partes

−fz −fz −fz
D2 D
2g
L ∫
. ze D
dz=
−2 g zD
L
.
f
e ( D
+ 2 e +K3
f )
Entonces reemplazando en la ecuación anterior tenemos:
−fz −fz −fz −fz
P2 P2 D2
e D
2
+ K 1− −e ( D
2
+ K 2 +( )
−2 g zD
L
.
f
e ( D
+ 2e
f
D
) + K 3 )=0

Entonces factorizando y ordenando tenemos:

−fz 2 2 2
e D
(
P P 2 g zD D
+ −
2 2 L
.
f (
+ 2 )+ K 1−K 2 + K 3=0
f )
−fz
2 g D2 fz
e D
( P 2−
Lf 2 ( )
. 1+ )+ K 1−K 2 + K 3=0
D

Donde K 1−K 2 + K 3=C

páá g. 3
 −fz
2 g D2
e D
( P 2−
Lf 2
. 1+
fz
D (=C )) , Entonces

fz
2 g D2 fz
2
P−
Lf 2 (D )
. 1+ =C . e D

dz
Como P= , Reemplazando en la ecuación anterior
dt
fz
dz 2 2 g D2 fz
( ) dt
=
Lf 2
. 1+ +C . e D
D ( )
La Constante de Integración se calcula bajo las condiciones en la cual z=z m

dz
Para =0
dt
fz
2 g D2 f zm
( )
m

0= . 1+ +C . e D
Lf 2 D
−f zm
−2 g D 2 f zm
C=
Lf
2
. 1+
D ( +C . e ) D

Finalmente:

m
z−z ¿
¿
f¿
¿
fz
fz
1+ − 1+ m . e¿
D
2
D ( )
dz 2 g D2
dt( ) =
Lf2
.¿

Tarea 2

Implementar un programa para un tiempo final de t=10” (Método de RK-4)

Un tubo en U tiene un diámetro de 500mm, f=0.03 y tiene una oscilación máxima,

Calcular como varían las posiciones y las velocidades hasta el tiempo final de t=10”.

Zm=6m, D=0.5m, F=0.03, L=20m

g=9.81m/s2, ∆ T =0.1 seg, c1=2*9.81/20=0.9810

c2=0.03/2*0.5=0.03

Condiciones iniciales: Zm=6m, Vm=0, t=0

De la ecuación diferencial
páá g. 4
 dv fv |v| 2 gz
+ + =0
dt 2 D L

dz
=v
dt

Entonces:

dv −fv |v| 2 gz
= −
dt 2D L

Haciendo

2 gz f
C1 = ,C 2=
L 2D

Finalmente

dv
=−C 1 z−C 2 v|v|
dt

Donde se implementará el algoritmo de RK4 para el sistema de ecuaciones (2) y (4)

dz
=v
dt

dv
=−C 1 z−C 2 v|v|
dt

1
Z n +1=Zn + [k 1 +2 k 2 +2 k 3 +k 4 ]
6

1
V n+1 =V n+ [ L1 +2 L2 +2 L3+ L4 ]
6

Haciendo unos cambios de variable, donde Z=x e V=y

Donde z es posición e y es la velocidad

Código en Matlab
format short
t0=input('Ingrese el valor de t0 inicial : ');
tf=input('Ingrese el valor de tf final : ');
h=input('Ingrese el tamaño de paso : ');
syms x y xi yi t z
f=input('Ingrese la función f : ');
g=input('Ingrese la función g : ');
d=(tf-t0)/h;
disp(' i t x y')
disp(' ----------------------')
r=zeros(d,4);
p=zeros(d,3);
t=t0;
x=input('Ingrese el valor inicial de x: ');
páá g. 5
y=input('Ingrese el valor inicial de y: ');
for i=1:1:d+1;
x0=x;
y0=y;
t0=t;

r(i,1)=i;
r(i,2)=t;
r(i,3)=x;
r(i,4)=y;
p(i,1)=t;
p(i,2)=x;
p(i,3)=y;
%PARA K1 L1
k1=eval(f);
L1=eval(g);
%PARA K2 L2
k2=[t0+(h/2),x0+(h/2)*k1,y0+(h/2)*L1];
t=k2(1);
x=k2(2);
y=k2(3);
k2=eval(f);
L2=[t0+(h/2),x0+(h/2)*k1,y0+(h/2)*L1];
L2=eval(g);
%Para K3 L3
k3=[t0+(h/2),x0+(h/2)*k2,y0+(h/2)*L2];
t=k3(1);
x=k3(2);
y=k3(3);
k3=eval(f);
L3=[t0+(h/2),x0+(h/2)*k2,y0+(h/2)*L2];
L3=eval(g);
%Para K4 L4
k4=[t0+(h),x0+(h)*k3,y0+(h)*L3];
t=k4(1);
x=k4(2);
y=k4(3);
k4=eval(f);
L4=[t0+(h),x0+(h)*k3,y0+(h)*L3];
L4=eval(g);
%PARA t x e y
x=x0+(h/6)*(k1+2*k2+2*k3+k4);
y=y0+(h/6)*(L1+2*L2+2*L3+L4);
t=t0+h;
end
%IMPRESION DE DATOS
r
plot(p)
xlabel('t')
ylabel('x ,y ')
title('Rk4 de 2 funciones')
Resultados

>> Untitled

Ingrese el valor de t0 inicial: 0

Ingrese el valor de tf final: 10

Ingrese el tamaño de paso: 0.1

Ingrese la función f: y

páá g. 6
Ingrese la función g: -(0.981*x+0.03*y*abs (y))

Ingrese el valor inicial de x: 6

Ingrese el valor inicial de y: 0

r=

i t x y

------------------------------------------

1.0000 0 6.0000 0

2.0000 0.1000 5.9706 -0.5873

3.0000 0.2000 5.8828 -1.1668

4.0000 0.3000 5.7378 -1.7308

5.0000 0.4000 5.5374 -2.2722

%Para ahorrar espacio quitamos algunos tramos de los resultados de 6 a 10

11.0000 1.0000 3.3602 -4.7170

12.0000 1.1000 2.8763 -4.9528

13.0000 1.2000 2.3715 -5.1339

14.0000 1.3000 1.8513 -5.2600

15.0000 1.4000 1.3213 -5.3314

%Para ahorrar espacio quitamos algunos tramos de los resultados de 16 a 20

21.0000 2.0000 -1.7737 -4.6998

22.0000 2.1000 -2.2311 -4.4404

23.0000 2.2000 -2.6606 -4.1448

24.0000 2.3000 -3.0589 -3.8164

25.0000 2.4000 -3.4229 -3.4584

%Para ahorrar espacio quitamos algunos tramos de los resultados de 26 a 30

31.0000 3.0000 -4.7527 -0.8775

32.0000 3.1000 -4.8169 -0.4064

33.0000 3.2000 -4.8339 0.0675

páá g. 7
 34.0000 3.3000 -4.8035 0.5403

35.0000 3.4000 -4.7261 1.0062

%Para ahorrar espacio quitamos algunos tramos de los resultados de 36 a 40

41.0000 4.0000 -3.3673 3.3570

42.0000 4.1000 -3.0174 3.6337

43.0000 4.2000 -2.6420 3.8691

44.0000 4.3000 -2.2450 4.0617

45.0000 4.4000 -1.8311 4.2104

%Para ahorrar espacio quitamos algunos tramos de los resultados de 46 a 50

51.0000 5.0000 0.7646 4.1865

52.0000 5.1000 1.1762 4.0403

53.0000 5.2000 1.5715 3.8585

54.0000 5.3000 1.9468 3.6435

55.0000 5.4000 2.2991 3.3978

%Para ahorrar espacio quitamos algunos tramos de los resultados de 56 a 60

61.0000 6.0000 3.7837 1.4334

62.0000 6.1000 3.9080 1.0512

63.0000 6.2000 3.9937 0.6610

64.0000 6.3000 4.0400 0.2660

65.0000 6.4000 4.0468 -0.1311

%Para ahorrar espacio quitamos algunos tramos de los resultados de 66 a 70

71.0000 7.0000 3.2848 -2.3172

72.0000 7.1000 3.0382 -2.6093

73.0000 7.2000 2.7639 -2.8715

74.0000 7.3000 2.4650 -3.1014

75.0000 7.4000 2.1448 -3.2969

%Para ahorrar espacio quitamos algunos tramos de los resultados de 76 a 80

páá g. 8
 81.0000 8.0000 -0.0207 -3.6991

82.0000 8.1000 -0.3879 -3.6385

83.0000 8.2000 -0.7473 -3.5440

84.0000 8.3000 -1.0956 -3.4171

85.0000 8.4000 -1.4297 -3.2596

%Para ahorrar espacio quitamos algunos tramos de los resultados de 86 a 90

91.0000 9.0000 -2.9846 -1.7929

92.0000 9.1000 -3.1486 -1.4837

93.0000 9.2000 -3.2810 -1.1628

94.0000 9.3000 -3.3808 -0.8328

95.0000 9.4000 -3.4473 -0.4962

96.0000 9.5000 -3.4799 -0.1558

97.0000 9.6000 -3.4784 0.1858

98.0000 9.7000 -3.4429 0.5251

99.0000 9.8000 -3.3736 0.8583

100.0000 9.9000 -3.2715 1.1813

101.0000 10.0000 -3.1378 1.4906 %Para t=10” (tiempo final)

Gráfica

páá g. 9
Tarea 3

Implementar un programa utilizando Euler modificado para el caso de la ecuación

de movimiento Resistencia Laminar

Código en Matlab
format short
t0=input('Ingrese el valor de t0 inicial: ');
tf=input('Ingrese el valor de tf final: ');
h=input('Ingrese el tamaño de paso: ');
syms x y xi yi t z
f=input('Ingrese la función f : ');
g=input('Ingrese la función g : ');
d=(tf-t0)/h;
disp(' i t x y ')
disp(' ----------------------')
r=zeros(d,4);
p=zeros(d,3);
t=t0;
x=input('Ingrese el valor inicial de x: ');
y=input('Ingrese el valor inicial de y: ');
for i=1:1:d+1;
x0=x;
y0=y;
t0=t;
r(i,1)=i;
r(i,2)=t;
r(i,3)=x;
r(i,4)=y;
p(i,1)=t;
p(i,2)=x;
p(i,3)=y;
%PARA K1 L1
páá g. 10
 k1=eval(f);
L1=eval(g);
%PARA K2 L2
k2=[t0+(h/2),x0+(h/2)*k1,y0+(h/2)*L1];
t=k2(1);
x=k2(2);
y=k2(3);
k2=eval(f);
L2=[t0+(h/2),x0+(h/2)*k1,y0+(h/2)*L1];
L2=eval(g);
x=x0+(h/2)*(k1+k2);
y=y0+(h/2)*(L1+L2);
t=t0+h;
end
%IMPRESION DE DATOS
r
plot(p)
xlabel('t')
ylabel('x , y ')
title('euler modificado de 2 funciones')

Aplicando a un problema de Resistencia laminar dado en el cuaderno

ft 2
Un tubo en U de 1 in de diámetro contiene aceite de υ=10 con una longitud
−4
s
total de 120 in se crea una diferencia manométrica de 16 in al aplicar aire
comprimido a uno de los meniscos, Al retirar súbitamente la presión del Aire hace
que la columna de Aceite Oscile libremente, Calcular como varía las posiciones
hasta encontrar la Velocidad máxima y en qué tiempo Donde V 0=1.9371ft/s y

16∈ ¿
12 ¿
ft
¿/2=0.667 ft
z=¿

De la ecuación diferencial de resistencia laminar

dv 32υ dz 2 gz
+ + =0
dt D2 dt L

Donde

dz
=v
dt

dv −32υ 2 gz
= v−
dt D
2
L

De Donde
páá g. 11
 dz dv
=v , =−C 1 v−C 2 z
dt dt

32 υ 2g
Donde: C1 = =0.4608 , C2 = =6.44
D2 L

Haciendo un cambio de variable Z=x y V=y, Haremos para un t 0=0 hasta t=10” con
tamaño de paso ∆ T =0.1

Resultados del Matlab

>> Untitled2

Ingrese el valor de t0 inicial: 0

Ingrese el valor de tf final: 10

Ingrese el tamaño de paso: 0.1

Ingrese la función f: y

Ingrese la función g: -(0.4608*y+6.44*x)

Ingrese el valor inicial de x: 0.667

Ingrese el valor inicial de y: 1.9371

r=

i t x y

-------------------------------------------

1.0000 0 0.6670 1.9371

2.0000 0.1000 0.8477 1.3931

3.0000 0.2000 0.9718 0.7675

4.0000 0.3000 1.0320 0.1016

5.0000 0.4000 1.0254 -0.5616

%Para ahorrar espacio quitamos algunos tramos de los resultados de 6 a 10

11.0000 1.0000 -0.0791 -2.4788

12.0000 1.1000 -0.3228 -2.2757

13.0000 1.2000 -0.5426 -1.9299

14.0000 1.3000 -0.7246 -1.4655

páá g. 12
 15.0000 1.4000 -0.8578 -0.9139

%Para ahorrar espacio quitamos algunos tramos de los resultados de 16 a 20

21.0000 2.0000 -0.4585 2.1346

22.0000 2.1000 -0.2401 2.2948

23.0000 2.2000 -0.0094 2.3062

24.0000 2.3000 0.2187 2.1700

25.0000 2.4000 0.4297 1.8970

%Para ahorrar espacio quitamos algunos tramos de los resultados de 26 a 30

31.0000 3.0000 0.7753 -1.1329

32.0000 3.1000 0.6508 -1.5566

33.0000 3.2000 0.4865 -1.8749

34.0000 3.3000 0.2933 -2.0690

35.0000 3.4000 0.0841 -2.1281

%Para ahorrar espacio quitamos algunos tramos de los resultados de 36 a 40

41.0000 4.0000 -0.7958 -0.1113

42.0000 4.1000 -0.7940 0.4022

43.0000 4.2000 -0.7415 0.8829

44.0000 4.3000 -0.6423 1.3004

45.0000 4.4000 -0.5034 1.6291

%Para ahorrar espacio quitamos algunos tramos de los resultados de 46 a 50

51.0000 5.0000 0.5519 1.1571

52.0000 5.1000 0.6574 0.7344

53.0000 5.2000 0.7194 0.2706

54.0000 5.3000 0.7346 -0.2040

55.0000 5.4000 0.7026 -0.6590

%Para ahorrar espacio quitamos algunos tramos de los resultados de 56 a 60

páá g. 13
 61.0000 6.0000 -0.1581 -1.6862

62.0000 6.1000 -0.3223 -1.4816

63.0000 6.2000 -0.4635 -1.1851

64.0000 6.3000 -0.5732 -0.8170

65.0000 6.4000 -0.6447 -0.4017

%Para ahorrar espacio quitamos algunos tramos de los resultados de 66 a 70

71.0000 7.0000 -0.2365 1.5919

72.0000 7.1000 -0.0753 1.6442

73.0000 7.2000 0.0884 1.5908

74.0000 7.3000 0.2443 1.4365

75.0000 7.4000 0.3823 1.1924

%Para ahorrar espacio quitamos algunos tramos de los resultados de 76 a 80

81.0000 8.0000 0.5020 -0.9861

82.0000 8.1000 0.3964 -1.2448

83.0000 8.2000 0.2670 -1.4204

84.0000 8.3000 0.1223 -1.5028

85.0000 8.4000 -0.0282 -1.4880

%Para ahorrar espacio quitamos algunos tramos de los resultados de 86 a 90

91.0000 9.0000 -0.5678 0.1346

92.0000 9.1000 -0.5454 0.4878

93.0000 9.2000 -0.4884 0.8049

94.0000 9.3000 -0.4009 1.0661

95.0000 9.4000 -0.2891 1.2556

96.0000 9.5000 -0.1603 1.3623

97.0000 9.6000 -0.0231 1.3803

98.0000 9.7000 0.1137 1.3099

99.0000 9.8000 0.2414 1.1568

páá g. 14
 100.0000 9.9000 0.3518 0.9318

101.0000 10.0000 0.4383 0.6504

Gráfica de los resultados de Matlab

páá g. 15