SciLab toxic pool simulator

Arxleol — Sat, 19 Sep 2009 12:22:14 +0000

This SciLab program (function) is solving or rather simulating flow of acid matter trough drainage pools…It simulates entrance of toxic water in the number one pool and then flow trough it into another pool. After which pools fill until full capacity! After reaching full capacity spillover from second pool flows into river. This simulator calculates grams exiting trough spillover from second pool into river.

function [izhod] = pools(b1,b2,dt,tmax)
//--------------------------
// Setup parameters
//--------------------------
//vh in grams/m3
vh=30
Qin=0.1;
Qr=0.5*Qin;
K=((0.2/24)/60)/60;
t=0;
B1fill=0;
Kb1=0;
Kb2=0;
B2fill=0;
izhod=0;
//--------------------------
// The main time loop
//--------------------------
while tb1 then
tmp=B1fill-b1;
B1fill=b1;
Kb1=Kb1-tmp*vh;
end
B2fill=B2fill+Qr*dt;
tmp=Qr*dt*(Kb1/B1fill);
Kb1=Kb1-tmp;
Kb2=Kb2+tmp;
B1fill=B1fill-Qr*dt;
if B2fill-Qin*dt&gt;=b2 then
B2fill=B2fill-Qin*dt;
tmp=Qin*dt;
Kb2=Kb2-tmp*(Kb2/B2fill);
izhod=tmp*(Kb2/B2fill);
end
Kb1=Kb1-(Kb1*K*dt);
 Kb2=Kb2-(Kb2*K*dt);
end
endfunction

Function takes parameters capacity of pools change in time interval and maximum time for following flow! Always make sure to input larger maximum time interval, so that flow has chance to stabilize. Otherwise results will not be correct.

In case that you would like to have faster execution you should rewrite while time loop into for time loop.Similar Posts:

SciLab oxygen saturation

Arxleol — Fri, 18 Sep 2009 11:52:54 +0000

Here is small program I wrote about oxygen saturation of water in the pool. It is in fact really slow however it points what needs to be done in similar situations.

Though this program executes slowly if one would rewrite while loop in the for loop then program should execute faster however I haven’t tested it.

function [] = kisik(Qo,dt)
//--------------------------
// Setup parameters
//--------------------------
i=1;
V=10;
So(i)=0;
t(i)=0;
Ka=0.1;
Sosat=8;
dSo=0;
//--------------------------
// The main time loop
//--------------------------
while So<Sosat
i=i+1;
dSo=(Ka*Qo*(Sosat-So(i-1))*dt)/V;
So(i)=So(i-1)+dSo;
t(i)=t(i-1)+dt;
end
plot2d(t,So)
endfunction

As you can see function takes two parameters. First one represents volume of water passing trough water pool. While second one dt is change in time interval.

Function will run until saturation level in water is same or larger then maximum saturation of oxygen in water. Formula that is user to calculate change in saturation, is second line in while loop.

Almost to forget after saturation levels are calculated graph is drawn of saturation levels on y-axis while on x-axis is time.Similar Posts:

Axino.net » Scilab

SciLab toxic pool simulator

SciLab oxygen saturation