blob: a209b74e02bd1d2657e899ba5840e910174a26c3 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
//Mass flow//
pathname=get_absolute_file_path('4.02.sce')
filename=pathname+filesep()+'4.02-data.sci'
exec(filename)
//If I=integral of(pV.dA):
//For system: ICS=Iab+Ibc+Icd+Ida
//But ICS=0
//For Aab:
function p=f(y),p=-d*U*w*y^0,endfunction
IAab=intg(0,t,f)
//For Acd:
function q=g(y),q=d*U*w*(2*y/t-(y/t)^2),endfunction
IAcd=intg(0,t,g)
//Mass flow rate across surface bc(in kg/sec):
mbc=(-IAab-IAcd)/1000
printf("\n\nRESULTS\n\n")
printf("\n\nMass flow rate across surface bc: %.4f kg/sec\n\n",mbc)
|
