Gruppe B224 D. Kildekode <strong>til</strong> BEV-modellen h[3]=enthalpy(air;v=v[3];p=p[3]) Specifik entalpi efter varme<strong>til</strong>førsel, [kJ/kg] Q ind=m luft*(h[3]-h[2]) Tilført varme fra forbrændingsproces, [kJ] ”!Ekspansionsproces ved konstant volumen” v[4]=v[1] Isochor ekspansionsproces s[3]=entropy(air;v=v[3];p=p[3]) Entropi efter varme<strong>til</strong>førsel, [kJ/kg*K] s[4]=s[3] Antagelse om isentropisk ekspansion p[4]=pressure(air;v=v[4];s=s[4]) Tryk efter ekspansion, [kPa] u[4]=intenergy(air;v=v[4];s=s[4]) Specifik indre energi efter ekspansion, [kJ/kg] T[4]=temperature(air;v=v[4];s=s[4]) Temperatur efter ekspansion, [K] Q ud=m luft*(u[4]-u[1]) Afgivet arbejde, [kJ] W net arbejde=Q ind-Q ud Net arbejdsoutput fra proces, [kJ] ”!Resultater vedrørende virkningsgrad og motoreffekt” Q ind=m diesel*bv*convert(MJ/kg;kJ/kg) Relater Q ind <strong>til</strong> faktuelt <strong>brændstof</strong>inputenergi, [kJ] m diesel=m luft/bra Brændstofandel i cylinderen - ikke korrigeret for volumetrisk virkningsgrad, [kg] m dot diesel=(rpm/60)*m diesel Tilført massestrøm af <strong>brændstof</strong> per sekund [kg/s] eta motor=convert(’-’;’%’)*W net arbejde/Q ind Motorens teoretiske virkningsgrad (fraregnet varmetab og mekaniske tab), [-] W dot motor=(W net arbejde/m luft)*m dot diesel W dot motor hk=W dot motor*convert(kW;hp) ”!Udregner Emissioner:” Lmin=1,86*((c*pro)+5,55*(h*pro)+0,7*(s*pro)-0,7*(o*pro))/0,21Luftforbruget[mˆ3 pr. kg Brændsel] vCO2=((c*pro)/12,0112)*22,26Volumen af CO2 i røggassen [mˆ3 pr. kg Brændsel] mCO2=vCO2*1,9768massen af CO2 i røggassen [kg/kg brændsel] vH2O=((h*pro)/2,0158)*22,4+1,24*(w*pro)Volumen af H2O i røggassen [mˆ3 pr. kg Brændsel] mH2O=vH2O*0,804massen af HO2 i røggassen [kg/kg brændsel] vSO2=((s*pro)/32,064)*21,89Volumen af SO2 i røggassen [mˆ3 pr. kg Brændsel] mSO2=vSO2*2,9263massen af SO2 i røggassen [kg/kg brændsel] vN2=0,79*lambda*LminVolumen af N2 i røggassen mN2=vN2*1,2505massen af N2 i røggassen [kg/kg brændsel] vO2=0,21*(lambda-1)*LminVolumen af O2 i røggassen [mˆ3 pr. kg Brændsel] mO2=vO2*1,42895massen af O2 i røggassen [kg/kg brændsel] ”!Uregning af tør og fugtig røggas” vt=vCO2+vSO2+vN2+vO2Volumen af tør røggas [mˆ3 pr. kg Brændsel] vf=vt+vH2OVolumen af fugtig røggas [mˆ3 pr. kg Brændsel] ”!Udregn kuldioxid CO2 i den tørre røggas udregnes i procent” CO2t=vCO2/vt 86
”!Brændværdi” HHV=(337*c-126*o+1420*h+93*s)/1000 Øvre brændværdi [MJ/kg] ”!Procedure” CALL brande(HHV;Br:bv) <strong>Aalborg</strong> <strong>Universitet</strong> 87