METODOS NUME

ejemplo 1

Un paracaidista con una masa de 68.1 kg salta de

un globo aerostático fijo. Aplique la ecuación g t
(1.10) para calcular la velocidad antes de que se 9.8 m/s^2
abra el paracaídas. Considere que el coeficiente m 68.1 Kg 0
de resistencia es igual a 12.5 kg/s. cr 12.5 kg/s 2
4
6
8
10
60 12
14
50

40

30

20

10

0
0 2 4 6 8 10 12 14 16
NUMERICOS
v(t)=(g*m/cr)*(1-e(-cr/m)*t)
0
16.4049808028706
27.7692914638702
35.6417515631291
41.0952832258206
44.8731375713484
47.4901909486417
49.3031164225182

12 14 16
 x f(x)=(x*100)/60
peso total peso del agua en el
cuerpo 100% cuerpo 60%
plasma 4.50% 7.50% 60.00%

tejidos conectivos densos

4.50% 7.50% 50.00%
y los cartilagos
linfa intersticial 12% 20%
40.00%
agua inaccesible en los
4.50% 7.50%
huesos
30.00%
agua intracelular 33% 55%
agua transcelular 1.50% 2.50% 20.00%

10.00%

0.00%
0.00% 5.00% 1
60.00%

50.00%

40.00%

30.00%

20.00%

10.00%

0.00%
0.00% 5.00% 10.00% 15.00% 20.00% 25.00% 30.00% 35.00%
estudiantes 30
volumen de cada estudiante 0.075 m^3
volumen del salon 240 m^3
calor de cada estudiante 80 W
Cv 0.718 Kj/(kg K)
temperatura inicial 20 °C
presion 101.325 kPa
peso molecular del aire 28.97 Kg/kmol
constante del gas ideal 8.314 kPa m^3/(kmol K)

V=Vs-(e*Ve) mA=P*V*Mwt/Re*Ti
237.75 4197.06424818078 kg

q=Ce*e 120000000
2400
100000000

80000000

60000000

40000000

20000000

0
0 2 4 6 8 10 12

Column I Column J
 t ΔT(t)=(mA*Cv*q*t)-(q*t)
1 7229981.11246512
2 14459962.2249302
3 21689943.3373954
4 28919924.4498605
5 36149905.5623256
6 43379886.6747907
7 50609867.7872558
8 57839848.899721
9 65069830.0121861
10 72299811.1246512
11 79529792.2371163
12 86759773.3495814
13 93989754.4620466
14 101219735.574512
15 108449716.686977

8 10 12 14 16

I Column J
 Este es un cálculo transitorio

balance = balance previo + deposito - retiros fecha depositos

1-May
balance jun = 1512.33 + 220.13 - 327.26 = 1405.2 220.13
1-Jun
balance jul = 1405.2 + 216.8 - 378.61 = 1243.39 216.8
1-Jul
balance ago= 1243.39 + 450.25 - 106.8 = 1586.84 450.25
1-Aug
balance set= 1586.84 + 127.31 - 350.61 = 1363.54 127.31
1-Sep

1800
1600
1400
1200
1000
800
600
400
200
0
19-Apr 9-May 29-May 18-Jun 8-Jul 28-Jul 17-Aug 6-Sep 26-Sep
retiros balance
1512.33
327.26
1405.2
378.61
1243.39
106.8
1586.84
350.61
1363.54
La tasa de flujo volumetrico a traves de un tubo esta dado por la
ecuacion , donde V es la velocidad promedio y A.

RECUERDA Q=V*A V=VELOCIDAD

Q1 = Q2 + Q3 A=AREA DE LA SECCION
TRANSVERSAL
∑▒ 〖𝐹𝐿𝑈𝐽𝑂 𝐸𝑁𝑇𝑅𝐴𝐷𝐴 = 〗 ∑▒〖𝐹𝐿𝑈𝐽𝑂 𝑆𝐴𝐿𝐼𝐷𝐴〗
40m3/s = 20m3/s + 6m/s*A

A=10/3m2
A=3.33m2
 ENTRADA SALIDA
40 20
0 20
EQUILIBRIO 40 40
 ENTRADAS
1.2
1
0.05
0.3
0
EQULIBRIO.SS 2.55

∑▒ 〖𝐹𝐿𝑈𝐽𝑂 𝐸𝑁𝑇𝑅𝐴𝐷𝐴 = 〗 ∑▒〖𝐹𝐿𝑈𝐽𝑂

𝐴𝐿𝐼𝐷𝐴〗
X+1+0,05+0,3+0=1,4+0,35+0,2+0,4+0,2
X=1,2
 SALIDAS
1.4
0.35
0.2
0.4
0.2
2.55

𝑅𝐴𝐷𝐴 = 〗 ∑▒〖𝐹𝐿𝑈𝐽𝑂 𝑆

+0=1,4+0,35+0,2+0,4+0,2
X=1,2
 PRIMER SUJETO
MASA= 70KG
COEF ARRASTRE= 12
tiempo velocidad
0 0
1 9.00601142
2 16.5932198
3 22.9851427
4 28.3700849
5 32.9066857
6 36.7285924
7 39.9483976
8 42.6609561
9 44.9461793
10 46.8713895
 SEGUNDO SUJETO
COEF ARRAST15
MASA= 75kg GRAFICA DEL ALCANCE
50
tiempo velocidad 45
0 0 40 velocidad tiempo velocidad
1 8.8821931 35
VELOCIDAD

2 16.1543177 30
25
3 22.1082298 20
4 26.9828808 15
5 30.9739074 10
6 34.2414836 5
0
7 36.9167488 0 2 4 6 8 10 12 14 16 18
8 39.1070706 TIEMPO
9 40.9003545
10 42.3685711
11 43.5706452
12 44.5548203
13 45.3605947
14 46.0203069
15 46.5604336
16 47.002652
ad

14 16 18
Utilice el cálculo para resolver la ecuación (1.9) para el caso
en que la velocidad inicial, v(0) es diferente de cero.

Conociendo los datos reemplazamos : v(t+1) = v(ti)+[g-(c/m)v(ti)](ti+1 - ti)

t,s v, m/s
0 1 veloc
1 10.61644640235 50
2 18.467759470934 45
3 24.8779357795 40
4 30.111501165054 35
5 34.384426795551 30
6 37.873041554077

v, m/s
25
7 40.721308522859
20
8 43.046765842451
15
9 44.94537710485
10
10 46.495491439496
5
0
0 2 4
1) = v(ti)+[g-(c/m)v(ti)](ti+1 - ti)

velocidad terminal
50
45
40
35
30
v, m/s

25
20
15
10
5
0
0 2 4 6 8 10 12
t,s
 Conociendo los datos reemplazamos : v(t+1) = v(ti)+[g-(c/m)v(ti)](ti+1 - ti)

El tiempo de paso de 1

t,s v, m/s velocidad terminal

0 0
60
1 9.8
2 17.8011747
50
3 24.3337051
4 29.667166
40
5 34.0216509
6 37.5768545

v, m/s
30
7 40.4794877
8 42.8493321
20
9 44.784183
10 46.3638851
10
11 47.6536272
12 48.7066325
0
0 2 4 6 8 10
t,s

En el tiempo de 10 segundos la velocidad es 46.3638851

/m)v(ti)](ti+1 - ti)

velocidad terminal

4 6 8 10 12 14
t,s
 Conociendo los datos reemplazamos : v(t+1) = v(ti)+[g-(c/m)v(ti)](ti+1 - ti)

El tiempo pasa de 0.5

t,s v, m/s
0 0 velovidad terminal
0.5 4.9 50
1 9.35291971
1.5 13.3995511 45
2 17.0769643
40
2.5 20.4188472
3 23.4558137 35
3.5 26.2156847
4 28.7237426 30

4.5 31.0029632
v, m/s
25
5 33.0742257
5.5 34.9565044 20
6 36.6670423
15
6.5 38.2215092
7 39.6341452 10
7.5 40.9178911
8 42.0845069 5
8.5 43.1446796
0
9 44.1081213 0 2 4 6
9.5 44.9836576 t,s
10 45.7793093
10.5 46.502365

En el tiempo de 10 segundos la velocidad es 45.7793093

El error :

a) Velocidad b) Velocidad

46.3638851 45.7793093
Error = 46.3638851 - 45.7793093 = 0.5845758 m/s
/m)v(ti)](ti+1 - ti)

velovidad terminal

4 6 8 10 12

t,s

793093
a) Con el empleo del cálculo, obtenga la solución de forma
cerrada para el caso en que al inicio el saltador se encuentra
en reposo (v = 0 en t = 0).

FU = –c′v^2

v(ti+1)=v(ti)+(g-(c/m)(v(ti))^2)*(ti+1+ti)

Reemplazando en la funcion:

v(1) = 0 + (9.8 - (12.5/68.1)*0^2)*0 = 0

b) Repita el cálculo numérico en el ejemplo 1.2 con los mismos
valores de condición inicial y de parámetros. Utilice un valor
de 0.225 kg/m para c′.

FU = –c′v^2

v(ti+1)=v(ti)+(g-(c/m)(v(ti))^2)*(ti+1+ti)

t,s v, m/s
0 0 velocidad
2 19.6 60
4 36.661497797357
6 47.380008232959 50
8 52.146097350876
10 53.77771319497 40
12 54.267300619108
14 54.407345223835 v, m/s
30
16 54.446821781287
18 54.45790282479
20
20 54.461009564387
22 54.461880294967
10
24 54.462124313051

0
0 5 10
 velocidad terminal

5 10 15 20 25 30
t,s
1.10. Calcule la velocidad de un paracaidista en caída libre con el empleo del método de Euler para el caso en que m = 80 kg
Lleve a cabo el cálculo desde t = 0 hasta t = 20 s con un tamaño de paso de 1 s. Use una condición inicial en que el paracaidis
velocidad hacia arriba de 20 m/s en t = 0. Suponga que el paracaídas se abre instantáneamente en t = 10 s, de modo que el c
arrastre sube a 50 kg/s.

i V(i) (m/s) m (kg) g (m/s) c (kg/s) e t(i+1)

0 20 80 9.8 10 2.71828182846 1
1 27.3 80 9.8 10 2.71828182846 2
2 33.6875 80 9.8 10 2.71828182846 3
3 39.2765625 80 9.8 10 2.71828182846 4
4 44.16699219 80 9.8 10 2.71828182846 5
5 48.44611816 80 9.8 10 2.71828182846 6
6 52.19035339 80 9.8 10 2.71828182846 7
7 55.46655922 80 9.8 10 2.71828182846 8
8 58.33323932 80 9.8 10 2.71828182846 9
9 60.8415844 80 9.8 10 2.71828182846 10
10 63.03638635 80 9.8 50 2.71828182846 11
11 33.43864488 80 9.8 50 2.71828182846 12
12 22.33949183 80 9.8 50 2.71828182846 13
13 18.17730944 80 9.8 50 2.71828182846 14
14 16.61649104 80 9.8 50 2.71828182846 15
15 16.03118414 80 9.8 50 2.71828182846 16
16 15.81169405 80 9.8 50 2.71828182846 17
17 15.72938527 80 9.8 50 2.71828182846 18
18 15.69851948 80 9.8 50 2.71828182846 19
19 15.6869448 80 9.8 50 2.71828182846 20
20 15.6826043 80 9.8 50 2.71828182846 21

Grafico del movimiento del paracaidista co

70

60

50

40

30
50

40

30

20

10

0
0 5 10 15
l caso en que m = 80 kg y c = 10 kg/s.
cial en que el paracaidista tiene una
= 10 s, de modo que el coeficiente de

Metodo Analitico
t (i) V(t+1) = V(i) + [g - (c/m)V(i)].(t(i+1)-t(i))
0 27.3
1 33.6875
2 39.2765625
3 44.1669921875
4 48.4461181640625
5 52.1903533935547
6 55.4665592193603
7 58.3332393169403
8 60.8415844023228
9 63.0363863520324
10 33.4386448820122
11 22.3394918307546
12 18.177309436533
13 16.6164910386999
14 16.0311841395124
15 15.8116940523172
16 15.7293852696189
17 15.6985194761071
18 15.6869448035402
19 15.6826043013276
20 15.6809766129978

del paracaidista con el Metodo Numerico

0 15 20 25
 1.11 En el ejemplo del paracaidista en caida libre se supuso que la b)Para el caso en que e
aceleracion debida a la gravedad era un valor constante de 9.8 m/s. Aunque regla de la cadaena par
esta es una buena aproxmacion cuando se estudiarn objetos e caida cerca de funcion de la altitud en
la superficie de la tierra la fuerza gravitacional disminuye conforme se es:
acerca al nivel del mar.𝑔_((𝑥)=) 𝑔_((0) ) (𝑅^2/(𝑅+𝑥)^2 ) 𝑑𝑣/𝑑𝑥=(𝑑𝑣.𝑑𝑥 )/(𝑑𝑥.𝑑
donde: SOLUCION: 𝑑𝑣/𝑑𝑥=−𝑔

g(x):aceleracion gravitacional
aceleracion gravitacional en la superficie :
r: radio de la superfifice 6.37x10^6m

c) Use el calculo para obtener la

a)En forma similar en que se obtuvo la ecuacion
𝑑𝑣/𝑑𝑡=𝑔−𝑐/𝑚 𝑣 use un balance de fuerzas para obtener una ∫128▒ 〖 1/(𝑎+𝑥)^2
ecuacion difeerencial para la velocidad como funcion de tiempo ⅆ𝑥 〗 =∫128▒ 〖 (𝑎+𝑥)^(−2
que itilice esta representacion mas completa de la gravitacion
.Sin embargo para esta obtencion suponga como positiva la ∫128▒ 〖𝑎𝑥 ^2 ⅆ𝑥 〗 =
velocidad hacia arriba SOLUCION:
𝑣=√((29101𝑅^2)/(𝑅+𝑥)+𝑣
 b)Para el caso en que el arrastre es despreciable ,utilice la
regla de la cadaena para expresar la ecuacion diferencial como
funcion de la altitud en luda de tiempo.recuerde que la regla
es:
𝑑𝑣/𝑑𝑥=(𝑑𝑣.𝑑𝑥 )/(𝑑𝑥.𝑑𝑡 )
SOLUCION: 𝑑𝑣/𝑑𝑥=−𝑔(0) 𝑅^2/(𝑣(𝑅+𝑥)^2 )

c) Use el calculo para obtener la forma cerra donde 𝑣=𝑣_0,en x=0

∫128▒ 〖 1/(𝑎+𝑥)^2
ⅆ𝑥 〗 =∫128▒ 〖 (𝑎+𝑥)^(−2)=(𝑎+𝑥)^(−1)+𝐶=1/(𝑎+𝑥)+𝐶 〗

∫128▒ 〖𝑎𝑥 ^2 ⅆ𝑥 〗 =𝑎𝑥/2+𝐶

SOLUCION:
𝑣=√((29101𝑅^2)/(𝑅+𝑥)+𝑣_0^2− 〖 2g_0 〗⁡𝑅 )
 1.12 La cantidad de un contaminante radiactivo distribuido uniformemente que se encuentra contenido en un reactor cerra
c (becquerel/litro, o Bq/L). El contaminante disminuye con una tasa de decaimiento proporcional a su concentración, es de
tasa de decaimiento = –kc
donde k es una constante con unidades de día–1. Entonces, de acuerdo con la ecuación (1.13), puede escribirse un balance
dc/dt = – kc
(cambio de la masa )=( disminución por decaimiento )
a) Use el método de Euler para resolver esta ecuación desde t = 0 hasta 1 d, con k = 0.2 d–1. Emplee un tamaño de paso de
Δt = 0.1. La concentración en t = 0 es de 10 Bq/L.
b) Grafique la solución en papel semilogarítmico (p.ej., ln c versus t) y determine la pendiente. Interprete sus resultados.

A) SOLUCIÓN B) SOLUCIÓN
c (i+1)=c (i)+Δt(-kc(i)) -1.82
El resultado expresa que por

t c Apox
0 10
MÉTODO
0.1 9.8 12
0.2 9.604
10
0.3 9.41192
0.4 9.2236816
8
0.5 9.039207968
0.6 8.8584238086 6
0.7 8.6812553325
0.8 8.5076302258 4
0.9 8.3374776213
1 8.1707280689 2

0
PENDIENTE Aprox -1.82837687 0 0.2 0.4 0.
entra contenido en un reactor cerrado, se mide por su concentración
porcional a su concentración, es decir:
1.13), puede escribirse un balance de masa para el reactor, así:

–1. Emplee un tamaño de paso de

iente. Interprete sus resultados.

SOLUCIÓN
.82
resultado expresa que por cada 0.1 de tiempo varia -1.82 su concentración

MÉTODOS EULER

0.2 0.4 0.6 0.8 1 1.2

PROBLEMA 1.13

Del problema se obtiene que: DATOS

(𝑑(𝐴𝑦))/𝑑𝑡=3𝑄 〖𝑠𝑒𝑛〗 ^2 (𝑡)−𝑄

Como el area es constante, entonces queda: A 1200 m^2

Q 500 m^3/d
𝑑𝑦/𝑑𝑡=(3𝑄 〖𝑠𝑒𝑛〗 ^2 (𝑡))/𝐴−𝑄" " /𝐴 Δt 0.5 d

Usando el metodo de Euler:

NUEVO VALOR = VALOR ANTIGUO + PENDIENTE*TAM. DE PASO

Obtenemos:

(𝑦_𝑓−𝑦_0)/Δ𝑡=(3𝑄 〖𝑠𝑒𝑛〗 ^2 (𝑡))/𝐴−𝑄" " /𝐴

𝑦_𝑓−𝑦_0=(3𝑄 〖𝑠𝑒𝑛〗 ^2 (𝑡)Δ𝑡)/𝐴−(𝑄" " Δ𝑡)/𝐴

𝑦_𝑓=𝑦_0+0.625∗ 〖𝑠𝑒𝑛〗 ^2 (𝑡)−0.208

TIEMPO(d) PROFUNDIDAD(m)
0 0
0.5 -0.06434447058379
1 0.170201415837188
1.5 0.584074071024828 Gráfica de Profundidad x Tiempo
2 0.892837702544706 2.5
2.5 0.908693269587448
3 0.713140055009209 2
3.5 0.582045600526926
4 0.732014361092118 1.5

4.5 1.12124256793108
1
5 1.48795242076747
5.5 1.5910693901462
0.5
6 1.4318650280423
6.5 1.25278790883911 0
7 1.31455752814917 0 2 4 6 8

7.5 1.65646000091755 -0.5

8 2.06022858851861
8.5 2.25071713165973
9 2.14886816033345
9.5 1.94439796715012
10 1.92137232283344
de Profundidad x Tiempo

4 6 8 10 12
PROBLEMA 1.14

Del problema se obtiene que: DATOS

𝑑𝑦/𝑑𝑡=(3𝑄 〖𝑠𝑒𝑛〗 ^2 (𝑡))/𝐴−(α 〖 (1+𝑦) 〗 ^1.5 " " )/𝐴

Usando el metodo de Euler: A 1200 m^2

NUEVO VALOR = VALOR ANTIGUO + PENDIENTE*TAM. DE PASO Q 500 m^3/d
Δt 0.5 d
Obtenemos:
α 300
(𝑦_𝑓−𝑦_0)/Δ𝑡=(3𝑄 〖𝑠𝑒𝑛〗 ^2 (𝑡))/𝐴−(α 〖 (1+𝑦_0) 〗 ^1.5 " " )/𝐴

𝑦_𝑓−𝑦_0=(3𝑄 〖𝑠𝑒𝑛〗 ^2 (𝑡)Δ𝑡)/𝐴−(α 〖 (1+𝑦_0) 〗 ^1.5 " " Δ𝑡)/𝐴

𝑦_𝑓=𝑦_0+0.625∗ 〖𝑠𝑒𝑛〗 ^2 (𝑡)−0.125 〖 (1+𝑦_0) 〗 ^1.5

TIEMPO(d) PROFUNDIDAD(m)
0 0
0.5 0.018655529416206 Gráfica de Profundidad x Tiempo
1 0.332687240595024 1.4
1.5 0.762249672286781
1.2
2 0.986590544554688
2.5 0.860442473124662 1
3 0.555688379946327
0.8
3.5 0.390048032716575
4 0.543157946234893 0.6
4.5 0.900764763531092
0.4
5 1.14790578478372
5.5 1.06553352008742 0.2
6 0.743256988082482
0
6.5 0.484471738250459 0 2 4 6 8 10
7 0.528158354273516
7.5 0.841924636332188
8 1.14121640977463
8.5 1.1480524939402
9 0.860673972902294
9.5 0.5469436940607
10 0.491414343884073
dad x Tiempo

8 10 12
PROBLEMA 1.15

Considerando la función: DATOS

𝑑𝑉/𝑑𝑡=−𝑘𝐴
k 0.1
Volumen de la gota: 𝑉=4/3 𝜋𝑟^3 Δt 0.25
r(0) 3
Area de la gota: 𝐴=4𝜋𝑟^2
v(0) 113.097
Al resolver 𝑟 y sustituir se obtiene: t(0) 0
𝐴= 〖 (4𝜋) 〗 ^(1/3) 〖 (3) 〗 ^(2/3) 𝑉^(2/3)
Reemplazando en la funcion:

𝑑𝑉/𝑑𝑡=−𝑘 〖 (4𝜋) 〗 ^(1/3) 〖 (3) 〗 ^(2/3) 𝑉^(2/3)

𝑑𝑉/𝑑𝑡=−0.4835𝑉^(2/3)

-Determinamos el volumen incial de la gota sabiedo que al inicio

tiene un 𝑟_0=3mm:

𝑉_0=4/3 𝜋 〖 (3) 〗 ^3
𝑉_0=113.097
El metodo de Euler sería:

𝑡_𝑛=𝑡_0+𝑛∗𝛥𝑡

𝑡_𝑛=0.25𝑛

𝑉_(𝑛+1)=𝑉_𝑛+𝛥𝑡∗"ƒ(" 𝑡_𝑛,𝑉_𝑛 ") n=0,1,2,3,..."

𝑉_(𝑛+1)=𝑉_𝑛−0.120875∗𝑉_𝑛^(2/3)
 TIEMPO(min) VOLUMEN(mm^3)
0 113.097 Grafico VOLUMEN
mm/min 0.25 110.270142756708 120
min 0.5 107.490588783019
mm 0.75 104.757942320528
mm^3 1 102.071807582663 100
min 1.25 99.4317887541975
1.5 96.83748999075
80
1.75 94.28851541827
2 91.7844691325112
2.25 89.3249551984898 60
2.5 86.9095776499288
2.75 84.5379404886859
40
3 82.2096476841667
3.25 79.9243031727205
3.5 77.6815108570198 20
3.75 75.4808746054216
4 73.3219982513107
0
4.25 71.2044855924234 0 2 4 6
4.5 69.1279403901519
4.75 67.0919663688275
5 65.096167214983
5.25 63.140146576592
5.5 61.2235080622856
5.75 59.3458552405448
6 57.5067916388668
6.25 55.7059207429068
6.5 53.9428459955904
6.75 52.2171707961989
7 50.528498499424
7.25 48.8764324143919
7.5 47.2605758036542
7.75 45.6805318821454
8 44.1359038161047
8.25 42.6262947219604
8.5 41.1513076651757
8.75 39.7105456590536
9 38.3036116634992
9.25 36.930108583737
9.5 35.5896392689817
9.75 34.2818065110592
10 33.0062130429765
 Grafico VOLUMEN X TIEMPO

2 4 6 8 10 12
problema 1.16
𝑑𝑇/𝑑𝑡=−𝑘(𝑇−𝑇_𝑎)
t T(t)
𝑑𝑇/((𝑇−21) )=−0.017𝑑𝑡 0 22
integramos 22.05
1 21.98314368
∫1▒ 〖𝑑𝑇 /((𝑇−21) )
2 21.9665715
=∫1▒ 〖− 0.017𝑑𝑡" " 〗〗 22
3 21.95027867
ln⁡〖 (𝑇−21) 〗 =−0.017𝑡 4 21.93426047 21.95
5 21.91851228
𝑒^(−0.017𝑡)=𝑇−21 21.9
6 21.90302955
𝑒^(−0.017𝑡)+21=𝑇 7 21.8878078
21.85
8 21.87284263
9 21.85812972
21.8
10 21.84366482
21.75
0 2 4
 T(t)
05

22

95

1.9

85

1.8

75
0 2 4 6 8 10 12
problema 1.17
a)
𝑑𝑁/𝑑𝑡=𝑢𝑁
𝑑𝑁/𝑁=𝑢𝑑𝑡
integramos
∫1▒ 〖𝑑𝑁 /𝑁=∫1▒𝑢𝑑𝑡" " 〗
ln⁡〖 (𝑁) 〗 =𝑢(𝑡−𝑡_𝑜)
𝑒^(𝑢(𝑡−𝑡_𝑜))=𝑁
t=20horas
a)
2𝑁=𝑒^(𝑢(20)" " ) taza de incremento (ln(2N))/20
ln⁡(2𝑁)=20𝑢 volumen de crecimiento
𝑢=(ln(2𝑁))/20 tiempo
b)
vol=4/3×𝜋×𝑟^3

𝑑𝑉/𝑑𝑡=𝑢
∫_(𝑟_𝑜)^(𝑟_ )▒𝑑𝑉=∫_(𝑡_𝑜)^(𝑡_ )▒𝑢𝑑𝑡

4/3 𝜋 〖 (𝑟−𝑟_𝑜 ) 〗 ^3=𝑢(𝑡−𝑡_𝑜 )+c

〖 (𝑟−𝑟_𝑜 ) 〗 ^3=3𝑢(𝑡−𝑡_𝑜 )/4𝜋

〖 (𝑟−𝑟_𝑜 ) 〗 ^1=∛(3𝑢(𝑡−𝑡_𝑜 )/4𝜋)

𝑟^1=∛(3𝑢(𝑡−𝑡_𝑜 )/4𝜋)+𝑟_𝑜
𝑟^1=∛(3𝑢(𝑡−𝑡_𝑜 )/4𝜋)+20
c)
500=∛(3𝑢/4𝜋)+20
480=∛(3𝑢/4𝜋)

480^3=3𝑢/4𝜋
𝑡=(480^3 4𝜋)/3𝑢
t=463246686.3/𝑢
b) c)

∛(3𝑢(𝑡−𝑡_𝑜 )/4𝜋) +21

463246686.3/u
problema 1.18
Se bombea un fluido por la red que se ilustra en la figura P1.18. Si Q2 = 0.6, Q
0.4, Q7 = 0.2 y Q8 = 0.3 m3/s, determine los otros flujos.
en la figura P1.18. Si Q2 = 0.6, Q3 =
tros flujos.
 En este caso utilizamos que :

El fluido que entra sera igual al fluido que

Desarrollando

Q1 1
Q2 0.6
Q3 0.4
Q4 0.1
Q5 0.3
Q6 0.1
Q7 0.2
Q8 0.3
Q9 0.4
Q10 1
igual al fluido que sale.

esarrollando

Tenemos
Q1 = Q2 + Q3
Q3 = Q1 - Q2
Q5 = Q3 - Q4
Q7 = Q5 - Q6
Q8 = Q6 + Q7
Q9 = Q4 + Q8
Q10 = Q2 + Q9