P. K. Nag. Engineering Thermodynamics, страница 7

1 b a rhf = 191.83; hfg = 2392.8; sf = 0.6493;sg = 8.1502;// At 100 d e g r e ehf100 = 419.04; hfg100 = 2257.0; sf100 = 1.3069;sg100 = 7.3549;// At 150 d e g r e ehf150 = 632.20; hfg150 = 2114.3; sf150 = 1.8418;sg150 = 6.8379;x2 = ( s1 - sf150 ) /4.9961;h2 = hf150 +( x2 * hfg150 ) ;x3 = ( s1 - sf100 ) /6.0480;h3 = hf100 +( x3 * hfg100 ) ;x4 = ( s1 - sf ) /7.5010;h4 = hf +( x4 * hfg ) ;h5 = hf ; h6 = h5 ;h7 = hf100 ; h8 = h7 ;h9 = 632.2; h10 = h9 ;m1 = ( h9 - h7 ) /( h2 - h7 ) ;m2 = ((1 - m1 ) *( h7 - h6 ) ) /( h3 - h6 ) ;Wt = 1*( h1 - h2 ) +(1 - m1 ) *( h2 - h3 ) +(1 - m1 - m2 ) *( h3 - h4 ) ;Q1 = h1 - h9 ;Wp = 0 ; // Pump work i s n e g l e c t e dn_cycle = 100*( Wt - Wp ) / Q1 ;sr = 3600/( Wt - Wp ) ;disp ( ” kJ / kg ” ,Wt , ” Net work p e r kg o s stem i s ” )disp ( ”%” , n_cycle , ” C y c l e e f f i c i e n c y i s ” )disp ( ” kg /kW h ” ,sr , ” Stream r a t e i s ” )Scilab code Exa 12.7 Calculations on expansion of steam in a turbine871234567891011121314151617181920212223Ti = 2000;Te = 450;T0 = 300;Q1_dot = 100 e03 ; // i n kWcpg = 1.1;wg = Q1_dot /( cpg *( Ti - Te ) ) ;af1 = wg * cpg * T0 *(( Ti / T0 ) -1 - log ( Ti / T0 ) ) ;af2 = wg * cpg * T0 *(( Te / T0 ) -1 - log ( Te / T0 ) ) ;afi = af1 - af2 ;h1 = 2801; h3 = 169; h4 = 172.8; h2 = 1890.2;s1 = 6.068; s2 = s1 ; s3 = 0.576; s4 = s3 ;Wt = h1 - h2 ;Wp = h4 - h3 ;Q1 = h1 - h4 ;Q2 = h2 - h3 ;Wnet = Wt - Wp ;ws = Q1_dot /2628;afu = 38*( h1 - h4 - T0 *( s1 - s3 ) ) ;I_dot = afi - afu ;Wnet_dot = ws * Wnet ;afc = ws *( h2 - h3 - T0 *( s2 - s3 ) ) ;n2 = 100* Wnet_dot / af1 ;disp ( ”%” ,n2 , ” The s e c o n d law e f f i c i e n c y i s ” )Scilab code Exa 12.8 Calculations on steam power plant123456789// P a r t ( a )h1 = 2758; h2 = 1817; h3 = 192; h4 = 200;Wt = h1 - h2 ; Wp = h4 - h3 ;Q1 = h1 - h4 ; Wnet = Wt - Wp ;n1 = Wnet / Wt ;WR = Wnet / Wp ;Q1_ = 100;PO = n1 * Q1_ ;cp = 1000;8810111213141516171819202122232425262728293031323334353637383940414243444546wg = ( Q1_ /(833 -450) ) ;EIR = wg * cpg *((833 -300) -300*( log (833/300) ) ) ;n2 = PO / EIR ;disp ( ” P a r t ( a ) ” )disp ( ”%” , n1 *100 , ” n1 i s ” )disp ( ”%” , n2 *100 , ” n2 i s ” )disp ( WR , ”Work r a t i o i s ” )// P a r t ( b )h1b = 3398; h2b = 2130; h3b = 192; h4b = 200;Wtb = 1268; Wpb = 8; Q1b = 3198;n1b = ( Wt - Wp ) / Q1 ;WRb = ( Wt - Wp ) / Wt ;EIRb = 59.3;Wnetb = Q1b * n1b ;n2b = Wnetb / EIRb ;disp ( ” P a r t ( b ) ” )disp ( ”%” , n1b *100 , ” n1 i s ” )disp ( ”%” , n2b *100 , ” n2 i s ” )disp ( WRb , ”Work r a t i o i s ” )// P a r t ( c )h1c = 3398; h2c = 2761; h3c = 3482; h4c = 2522; h5c= 192; h6c = 200;Wt1 = 637; Wt2 = 960; Wtc = Wt1 + Wt2 ; Wpc = 8;Wnetc = 1589; Q1c = 3198+721;n1c = Wnetc / Q1c ;WRc = Wnetc / Wtc ;POc = Q1_ * n1c ;EIRc = 59.3;n2c = POc / EIRc ;disp ( ” P a r t ( c ) ” )disp ( ”%” , n1c *100 , ” n1 i s ” )disp ( ”%” , n2c *100 , ” n2 i s ” )disp ( WRc , ”Work r a t i o i s ” )// P a r t ( d )T3 = 318.8; T1 = 568;n1d = 1 -( T3 / T1 ) ;Q1d = 2758 -1316;Wnet = Q1d * n1d ;8947 Wpd = 8; Wtd = 641;48 WRd = ( Wt - Wp ) / Wt ;49 POd = Q1_ *0.439;50 EIRd = ( Q1_ /(833 -593) ) * cpg *((833 -300) -300*( log5152535455(833/300) ) ) ;n2d = POd / EIRd ;disp ( ” P a r t ( d ) ” )disp ( ”%” , n1d *100 , ” n1 i s ” )disp ( ”%” , n2d *100 , ” n2 i s ” )disp ( WRd , ”Work r a t i o i s ” )Scilab code Exa 12.9 Calculations on steam in a chemical plant12345678910111213141516171819202122hfg = 2202.6;Qh = 5.83;ws = Qh / hfg ;eg = 0.9; // e f f i c i e n c y o f g e n e r a t o rP = 1000;Wnet = 1000/0.9;nbrake = 0.8;h1_2s = Wnet /( ws * nbrake ) ; // h1−h 2 sn_internal = 0.85;h12 = n_internal * h1_2s ;hg = 2706.3; h2 = hg ;h1 = h12 + h2 ;h2s = h1 - h1_2s ;hf = 503.71;x2s = ( h2s - hf ) / hfg ;sf = 1.5276; sfg = 5.6020;s2s = sf +( x2s * sfg ) ;s1 = s2s ;P1 = 22.5; // i n b a r from M o i l l e r c h a r tt1 = 360;disp ( ” d e g r e e ” ,t1 , ” T e m p e r a t u r e o f t h e steam i s ” )disp ( ” b a r ” ,P1 , ” P r e s s u r e o f t h e steam i s ” )90Scilab code Exa 12.10 Calculation of oil consumption per day in a factory12345678910111213141516171819h1 = 3037.3; h2 = 561+(0.96*2163.8) ;s2 = 1.6718+(0.96*5.3201) ;s3s = s2 ;x3s = ( s3s -0.6493) /7.5009;h3s = 191.83+( x3s *2392.8) ;h23 = 0.8*( h2 - h3s ) ; // h2−h3h3 = h2 - h23 ;h5 = 561.47; h4 = 191.83;Qh = 3500; // i n kJ / sw = Qh /( h2 - h5 ) ;Wt = 1500;ws = ( Wt + w *( h2 - h3 ) ) /( h1 - h3 ) ;ws_ = 3600* ws ; // i n kg / hh6 = (( ws - w ) * h4 + w * h5 ) / ws ;h7 = h6 ;n_boiler = 0.85;CV = 44000; // i n kJ / kgwf = (1.1* ws_ *( h1 - h7 ) ) /( n_boiler * CV ) ;disp ( ” kg / h ” ,wf , ” F u e l b u r i n g r a t e i s ” )Scilab code Exa 12.11 Calculations on a steam turbine1234567h1 = 3285; h2s = 3010; h3 = 3280; h4s = 3030;h4 = h3 -0.83*( h3 - h4s ) ;h5s = 2225;h5 = h4 -0.83*( h4 - h5s ) ;h6 = 162.7; h7 = h6 ;h8 = 762.81;h2 = h1 -0.785*( h1 - h2s ) ;918 m = ( h8 - h7 ) /( h4 - h7 ) ;9 n_cycle = (( h1 - h2 ) +( h3 - h4 ) +(1 - m ) *( h4 - h5 ) ) /(( h1 - h8 ) +(h3 - h2 ) )10 disp ( ” kg / s ” ,m , ” Steam f l o w a t t u r b i n e i n l e t i s ” )11 disp ( ”%” , n_cycle *100 , ” c y c l e e f f i c i e n c y i s ” )Scilab code Exa 12.12 Calculations on a binary vapour cycle123456789101112131415161718192021222324// From t a b l e and g r a p hh1 = 2792.2;h4 = 122.96;hb = 254.88;hc = 29.98;ha = 355.98;hd = hc ;h2 = 1949.27;//m = ( h1 - h4 ) /( hb - hc ) ; // Amount o f m e r c u r ycirculatingQ1t = m *( ha - hd ) ;W1t = m *( ha - hb ) + ( h1 - h2 ) ;Nov = W1t / Q1t ;disp ( ”%” , Nov *100 , ” O v e r a l l e f f i c i e n c y o f t h e c y c l e ” )S = 50000; // Stem f l o w r a t e t h r o u g h t u r b i n e i n kg /hwm = S * m ;disp ( ” kg /h ” ,wm , ” Flow t h r o u g h t h e m e r c u r y t u r b i n e i s ”)Wt = W1t * S /3600;disp ( ”kW” ,Wt , ” U s e f u l work done i n b i n a r y v a p o u rc y c l e i s ”)nm = 0.85; // I n t e r n a l e f f i c i e n c y o f m e r c u r y t u r b i n ens = 0.87; // I n t e r n a l e f f i c i e n c y o f steam t u r b i n eWTm = nm *( ha - hb ) ;hb_ = ha - WTm ; // hb ’m_ = ( h1 - h4 ) /( hb_ - hc ) ; // m’922526272829303132h1_ = 3037.3; // h ’Q1t = m_ *( ha - hd ) +( h1_ - h1 ) ;x2_ = (6.9160 -0.4226) /(8.47 -0.4226) ;h2_ = 121+(0.806*2432.9) ;WTst = ns *( h1_ - h2_ ) ;WTt = m_ *( ha - hb_ ) + WTst ;Nov = WTt / Q1t ;disp ( ”%” , Nov *100 , ” O v e r a l l e f f i c i e n c y i s ” )93Chapter 13Gas power cycleScilab code Exa 13.1 Calculations on otto cycle123456789101112131415161718192021T1 = 273+35;P1 = 100 e03 ; // i n kN/m2Q1 = 2100;R = 0.287;v1 = 0.884; v2 = 0.11; v3 =rk = 8; g = 1.4; // gamman_cycle = 1 -(1/ rk ^(1.4 -1) ) ;v12 = 8; // v1 / v2v1 = ( R * T1 ) / P1 ;v2 = v1 /8;T2 = T1 *( v1 / v2 ) ^( g -1) ;cv = 0.718;T3 = Q1 / cv + T2P21 = ( v1 / v2 ) ^ g ;P2 = P21 * P1 ;P3 = P2 *( T3 / T2 ) ;Wnet = Q1 * n_cycle ;Pm = Wnet /( v1 - v2 ) ;disp ( ”MPa” , P3 /1 e06 , ”Maximumdisp ( ”K” ,T3 , ” T e m p e r a t u r e o fdisp ( ”%” , n_cycle *100 , ” C y c l e94v2 ;p r e s s u r e i s ”)the c y c l e i s ”)e f f i c i e n c y i s ”)22disp ( ”MPa” , Pm /1 e06 , ”Mean e f f e c t i v e p r e s s u r e i s ” )Scilab code Exa 13.2 Calculations on a diesel engine12345rk = 14;k = 0.06rc = k *(14 -1) +1;g = 1.4;n_diesel = 1 -((1/ g ) ) *(1/ rk ^( g -1) ) *(( rc ^( g -1) ) /( rc -1));6 disp ( ”%” , n_diesel *100 , ” A i r s t a n d a r d e f f i c i e n c y i s ” )Scilab code Exa 13.3 Calculations on air standard diesel cycle123456789101112131415161718rk = 16;T1 = 273+15;P1 = 100; // i n kN/m2T3 = 1480+273;g = 1.4; // gammaR = 0.287;T2 = 288*( rk ^( g -1) ) ;rc = T3 / T2 ;cp = 1.005; cv = 0.718;Q1 = cp *( T3 - T2 ) ;T4 = T3 *(( rc / rk ) ^( g -1) ) ;Q2 = cv *( T4 - T1 ) ;n = 1 -( Q2 / Q1 ) ; // c y c l e e f f i c i e n c yn_ = 1 -((1/ g ) ) *(1/ rk ^( g -1) ) *(( rc ^( g -1) ) /( rc -1) ) ; //c y c l e e f f i c i e n c y from a n o t h e r f o r m u l aWnet = Q1 * n ;v1 = ( R * T1 ) / P1 ;v2 = v1 / rk ;Pm = Wnet /( v1 - v2 ) ;9519202122disp ( rc , ” cut − o f f r a t i o i s ” )disp ( ” kJ / kg ” ,Q1 , ” Heat s u p p l i e d p e r kg o f a i r i s ” )disp ( ”%” ,n *100 , ” C y c l e e f f i c i e n c y i s ” )disp ( ”KPa” ,Pm , ”Mean e f f e c t i v e p r e s s u r e i s ” )Scilab code Exa 13.4 Calculations on air standard dual cycle1234567891011121314151617181920212223T1 = 273+50;rk = 16;g = 1.4; // gammaP3 = 70; cv = 0.718; cp = 1.005; R = 0.287;T2 = T1 *(( rk ^( g -1) ) ) ;P1 = 1; // i n b a rP2 = P1 *( rk ) ^ g ;T3 = T2 *( P3 / P2 ) ;Q23 = cv *( T3 - T2 ) ;T4 = ( Q23 / cp ) + T3 ;v43 = T4 / T3 ; // v4 / v3v54 = rk / v43 ; // v5 / v4 = ( v1 / v2 ) ∗ ( v3 / v4 )T5 = T4 *(1/ v54 ) ^( g -1) ;P5 = P1 *( T5 / T1 ) ;Q1 = cv *( T3 - T2 ) + cp *( T4 - T3 ) ;Q2 = cv *( T5 - T1 ) ;n_cycle = 1 -( Q2 / Q1 ) ;v1 = ( R * T1 ) / P1 ;v12 = (15/16) * v1 ; // v1−v2Wnet = Q1 * n1 ;Pm = Wnet /( v12 ) ;disp ( ”%” ,n *100 , ” E f f i c i e n c y o f t h e c y c l e i s ” )disp ( ” b a r ” ,Pm , ”Mean e f f e c t i v e p r e s s u r e i s ” )Scilab code Exa 13.5 finding the increase in cycle efficiency of gas turbineplant9612345678910111213141516171819202122P1 = 0.1 e06 ;T1 = 303;T3 = 1173;PR = 6; // P r e s s u r e r a t i orp = 6; nt = 0.8; nc = 0.8;g = 1.4; cv = 0.718; cp = 1.005; R = 0.287;j = ( PR ) ^(( g -1) / g ) ;T2s = j * T1 ;T4s = T3 / j ;T21 = ( T2s - T1 ) / nc ; // T2−T1T34 = nt *( T3 - T4s ) ; // T3−T4Wt = cp * T34 ;Wc = cp * T21 ;T2 = T21 + T1 ;Q1 = cp *( T3 - T2 ) ;n = ( Wt - Wc ) / Q1 ;T4 = T3 -375;T6 = 0.75*( T4 - T2 ) + T2 ;Q1_ = cp *( T3 - T6 ) ;n_ = ( Wt - Wc ) / Q1_ ;I = ( n_ - n ) / n ;disp ( ”%” ,I *100 , ” The p e r c e n t a g e e f f i c i e n c y i n c y c l ee f f i c i e n c y due t o r e g e n e r a t i o n i s ” )Scilab code Exa 13.6 Calculations on gas turbine plant operating on bryton cycle1 cp = 1.005;2 Tmax = 1073; Tmin = 300;3 Wnet_max = cp *( sqrt ( Tmax ) - sqrt ( Tmin ) ) ^2;4 n_cycle = 1 - sqrt ( Tmin / Tmax ) ;5 n_carnot = 1 -( Tmin / Tmax ) ;6 r = n_cycle / n_carnot ;7 disp ( ” kJ / kg ” , Wnet_max , ”Maximum work done p e r kg o fa i r i s ”)9789disp ( ”%” , n_cycle *100 , ” c y c l e e f f i c i e n c y i s ” )disp (r , ” r a t i o o f b r a y t o n and c a r n o t e f f i c i e n c yi s ”)Scilab code Exa 13.7 Calculations on an ideal bryton cycle123456789101112131415rp = 6;g = 1.4; cv = 0.718; cp = 1.005; R = 0.287;T1 = 300; T3 = 1100; T0 = 300;n_cycle = 1 -(1/ rp ^(( g -1) / g ) ) ;j = rp ^(( g -1) / g ) ;T2 = T1 * j ;T4 = T3 / j ;Wc = cp *( T2 - T1 ) ;Wt = cp *( T3 - T4 ) ;WR = ( Wt - Wc ) / Wt ;Q1 = 100; // i n MWPO = n_cycle * Q1 ;m_dot = ( Q1 *1 e06 ) /( cp *( T3 - T2 ) ) ;R = m_dot * cp * T0 *(( T4 / T0 ) -1 - log ( T4 / T0 ) ) ;disp ( ”%” , n_cycle *100 , ” The t h e r m a l e f f i c i e n c y o f t h ec y c l e i s ”)16 disp ( WR , ”Work r a t i o i s ” )17 disp ( ”MW” ,PO , ” Power o u t p u t i s ” )18 disp ( ”MW” ,R /1 e06 , ” Energy f l o w r a t e o f t h e e x h a u s tgas stream i s ”)Scilab code Exa 13.8 Calculations on stationary gas turbine12345nc = 0.87; nt = 0.9; T1 = 311;rp = 8; // P2/P1P1 = 1; P2 = 8; P3 = 0.95* P2 ; P4 = 1;g = 1.4; cv = 0.718; cp = 1.005; R = 0.287;// With no c o o l i n g9867891011121314T2s = T1 *(( P2 / P1 ) ^(( g -1) / g ) ) ;T2 = T1 + ( T2s - T1 ) /0.87;T4s = T3 *( P4 / P3 ) ^(( g -1) / g ) ;n = ((( T3 - T4s ) * nt ) -(( T2s - T1 ) / nc ) ) /( T3 - T2 ) ;// With c o o l i n gn_cycle = n -0.05;x = 0.13;r = 0.13/1.13;disp ( ”%” ,r *100 , ” P e r c e n t a g e o f a i r t h a t may be t a k e nfrom t h e c o m p r e s s o r i s ” )Scilab code Exa 13.10 Calculations on air flying through the engine of aturbojet aircraft1234567891011121314151617T1 = 233; V1 = 300; cp = 1.005; g = 1.4;T2 = T1 +(( V1 ^2) /(2* cp ) ) *1 e -03 ;P1 = 35;P2 = P1 *( T2 / T1 ) ^( g /( g -1) ) ;rp = 10; // P r e s s u r e r a t i oP3 = rp * P2 ;T3 = T2 *( P3 / P2 ) ^(( g -1) / g ) ;T4 = 1373;T5 = T4 - T3 + T2 ;P4 = P3 ;P5 = P4 *( T5 / T4 ) ^( g /( g -1) ) ;disp ( ”K” ,T5 , ” T e m p e r a t u r e a t t h e t u r b i n e e x i t i s ” )disp ( ” kPa ” ,P5 , ” P r e s s u r e a t t h e t u r b i n e e x i t i s ” )P6 = P1 ;T6 = T5 *( P6 / P5 ) ^(( g -1) / g ) ;V6 = (2* cp *1000*( T5 - T6 ) ) ^0.5 ;disp ( ”m/ s ” ,V6 , ” V e l o c i t y o f t h e g a s a t t h e n o z z l ee x i t i s ”)18 w = 50;19 Ve = V6 ; Vi = 300;20 Wp_dot = w * Vi *( Ve - Vi ) ;9921 h4 = 1373; h3 = 536.66;22 Q1 = w * cp *( h4 - h3 ) ; // i n kJ / kg23 np = Wp_dot /( Q1 *1000) ;24 disp ( ”%” , np *100 , ” The p r o p u l s i v ee f f i c i e n c y of thec y c l e i s ”)Scilab code Exa 13.11 Calculations on a combined GT ST plant1 Ta = 288;2 rp = 8; // Pb/Pa3 g = 1.33; g1 = 1.44; cv = 0.718; cpa = 1.005; cpg =45678910111213141516171819202122231.11; R = 0.287;Tb = Ta *( rp ) ^(( g1 -1) / g1 ) ;Tc = 1073; Tm = 800+273; Tmin = 100+273;Td = Tc /( rp ^(( g -1) / g ) ) ;Wgt = cpg *( Tc - Td ) - cpa *( Tb - Ta ) ;Q1 = cpg *( Tc - Tb ) ;Q1_ = cpg *( Tc - Td ) ;h1 = 3775; h2 = 2183; h3 = 138; h4 = h3 ;Q1_st = h1 - h3 ; // Q1 ’Q_fe = cpg *( Tm - Tmin ) ;was = Q1_st / Q_fe ; // wa/ wsWst = h1 - h2 ;PO = 190 e03 ; // i n kWws = PO /( was * Wgt + Wst ) ;wa = was * ws ;CV = 43300; // i n kJ / kgwaf = CV /( Q1 + Q1_ ) ;FEI = ( wa / waf ) * CV ;noA = PO / FEI ;disp ( waf , ” A i r f u e l r a t i o i s ” )disp ( ”%” , noA *100 , ” O v e r a l l e f f i c i e n c y o f combinedplant i s ”)100Chapter 14Refrigeration cycleScilab code Exa 14.1 Finding the power required to drive a cold storageplant123456T2 = 268; T1 = 308;COP = T2 /( T1 - T2 ) ;ACOP = COP /3; // A c t u a l COPQ2 = 29; // i n kWW = Q2 / ACOP ;disp ( ”kW” ,W , ” Power r e q u i r e d t o d e r i v e t h e p l a n e i s ” )Scilab code Exa 14.2 Heat calculations on a refrigerator12345678h1 = 236.04; s1 = 0.9322; s2 = s1 ;P2 = 0.8; // i n MPah2 = 272.05; h3 = 93.42; h4 = h3 ;m = 0.06; // mass f l o w r a t eQ2 = m *( h1 - h4 ) ;Wc = m *( h2 - h1 ) ;Q1 = m *( h2 - h4 ) ;COP = Q2 / Wc ;101disp ( ”kW” ,Q2 , ” The r a t e o f h e a t r e m o v a l i s ” )disp ( ”kW” ,Wc , ” Power i n p u t t o t h e c o m p r e s s o r i s ” )disp ( ”kW” ,Q1 , ” The h e a t r e j e c t i o n r a t e i n t h econdenser i s ”)12 disp ( COP , ”COP i s ” )91011Scilab code Exa 14.3 Calculations on refrigeration by a simple R 12 plant12345678910111213141516171819202122232425h1 = 183.19; h2 = 209.41; h3 = 74.59; h4 = h3 ;T1 = 313; T2 = 263;W = 70000/3600; // P l a n t c a p a c i t y i n kWw = W /( h1 - h4 ) ; // R e f r i g e r a n t f l o w r a t ev1 = 0.077;VFR = w * v1 ;T = 48; // i n d e g r e eP2 = 9.6066; P1 = 2.1912;rp = P2 / P1 ; // P r e s s u r e r a t i oQ1 = w *( h2 - h3 ) ;hf = 26.87; hfg = 156.31;x4 = ( h4 - hf ) / hfg ;COP = ( h1 - h4 ) /( h2 - h1 ) ;PI = w *( h2 - h1 ) ;COP = T2 /( T1 - T2 ) ;COP_v = 4.14;r = COP_v / COP ;disp ( ” kg / s ” ,w , ” R e f r i g e r a n t f l o w r a t e i s ” )disp ( ”m3/ s ” ,VFR , ” Volume f l o w r a t e i s ” )disp ( ” d e g r e e ” ,T , ” C o m p r e s s o r d i s c h a r g e t e m p e r a t u r e i s”)disp ( rp , ” P r e s s u r e r a t i o i s ” )disp ( ”kW” ,Q1 , ” Heat r e j e c t e d t o t h e c o n d e n s e r i s ” )disp ( ”%” , x4 *100 , ” F l a s h g a s p e r c e n t a g e i s ” )disp ( COP , ”COP i s ” )disp ( ”kW” ,PI , ” Power r e q u i r e d t o d r i v e t h e c o m p r e s s o ri s ”)10226disp (r , ” R a t i o o f COP o f c a r n o t r e f r i g e r a t o ri s ”)Scilab code Exa 14.4 Calculations on R 12 vapour compression plant1234567891011121314151617181920212223h3 = 882; h2 = 1034;h6 = 998; h1 = 1008;v1 = 0.084;h4 = h3 - h1 + h6 ; h5 = h4 ;t4 = 25+273;disp ( ” kJ / kg ” ,h6 - h5 , ” R e f r i g e r a t i o n e f f e c t i s ” )m = 10;w = ( m *14000) /(( h6 - h5 ) *3600) ; // i n kg / sdisp ( ” kg / s ” ,w , ” R e f r i g e r a n t f l o w r a t e i s ” )v1 = 0.084;VFR = w *3600* v1 ; // i n kg / hve = 0.8; // v o l u m e t r i c e f f i c i e n c yCD = VFR /( ve *60) ; // i n m3/ minN = 900;n = 2;D = (( CD *4) /( %pi *1.1* N * n ) ) ^(1/3) ; // L = 1 .