METODO DE RUNGE KUTTA para los primeros 5 min

Datos iniciales:to=0 sDo=0.12mh=5msQe=140 l/s

k1k2k3k4kt(s)h(m)

-----05

0.048491930.047920740.047927450.047366160.04792575205.04792575

0.047366200.046814820.046821220.046279240.04681959405.09474534

0.046279280.045746710.045752820.045229190.04575126605.14049659

0.045229220.044714540.044720380.044214190.04471887805.18521547

0.044214220.043716550.043722140.043232550.043720691005.22893616

0.043232590.042751120.042756460.042282700.042755071205.27169123

0.042282730.041816710.041821830.041363160.041820491405.31351173

0.041363190.040911890.040916800.040472520.040915521605.35442724

0.040472550.040035320.040040030.039609500.040038791805.39446603

0.039609530.039185730.039190250.038772860.039189062005.43365509

0.038772890.038361930.038366280.037961460.038365132205.47202022

0.037961480.037562820.037566990.037174200.037565882405.50958610

0.037174220.036787320.036791340.036410060.036790272605.54637637

0.036410080.036034450.036038320.035668080.036037292805.58241365

0.035668100.035303280.035307000.034947350.035306003005.61771966

TABLA N1

Se observa que pasado los 5 primeros minutos la altura h=6.61771966m, el cual ser la condicin inicial para la segunda parte del anlisis donde el caudal aumentar 5 l/s y se tomar desde el tiempo de 5min hasta los 10 min.

METODO DE RUNGE KUTTA para el intervalo de 5 min a 10 min

Datos iniciales:to=300 sDo=0.12mh=5msQe=145 l/s

TABLA N2

k1k2k3k4kt(s)h(m)

-----3005.61771966

0.043164810.042712220.042716960.042271250.042715743205.66043539

0.042271270.041832460.0418370.041404760.041835823405.70227122

0.041404790.040979140.040983510.040564150.040982373605.74325359

0.040564180.040151140.040155330.039748310.040154243805.78340783

0.039748340.039347360.039351390.038956190.039350344005.82275816

0.038956210.038566790.038570680.038186790.038569664205.86132782

0.038186810.037808480.037812220.037439190.037811234405.89913905

0.037439210.037071510.037075110.03671250.037074164605.93621321

0.036712520.036355030.03635850.036005890.036357584805.97257079

0.036005910.035658220.035661570.035318570.035660685006.00823147

Se observa que hay presencia de desborde entre los tiempo t=480 seg y t=500 seg, por lo cual para poder determinar el tiempo en el cual se desborda se procede a interpolar los valores en dichos tiempos.

El sistema se desbordara en t= 495.3834 seg