108p=pTarray(:,1);T=pTarray(:,2);hc=hcu;C2=0;case 'combu'p=pTarray(:,1);T=pTarray(:,3);hc=hcu;C2=3.24e-3;case 'combb'p=pTarray(:,1);T=pTarray(:,2);hc=hcb;C2=3.24e-3;case 'peri'p=pTarray(:,1);T=pTarray(:,2);hc=hcb;C2=3.24e-3;case 'in<strong>water</strong>'p=pTarray(:,1);T=pTarray(:,2);hc=hcb;C2=3.24e-3;case 'exp'p=pTarray(:,1);T=pTarray(:,2);hc=hcb;C2=3.24e-3;endswitch heattransferlawcase 'constant'hcoeff=hc;case 'Woschni'V=Vtdc*(1+(r-1)/2*(1-cos(theta)+ ...1/eps*(1-(1-eps^2*sin(theta).^2).^0.5)));upmean=omega*stroke/pi; % mean piston velocityVs=Vbdc-Vtdc;k=1.3;C1=2.28;pm=p1*(V1./V).^k; % motoring pressurehcoeff=hc*130*b^(-0.2)*T.^(-0.53).*(p/100e3).^(0.8).* ...(C1*upmean+C2*Vs*T1/p1/V1*(p-pm)).^(0.8);endq=hcoeff.*(T-Tw);D-14. M-file for Xstream.m
109function Out=XSteam(fun,In1,In2)%Cp_pT Specific isobaric heat capacity as a function <strong>of</strong> pressure and <strong>temperature</strong>.case 'cp_pt'p = toSIunit_p(In1);T = toSIunit_T(In2);Region = region_pT(p, T);switch Regioncase 1Out = fromSIunit_Cp(Cp1_pT(p, T));case 2Out = fromSIunit_Cp(Cp2_pT(p, T));case 3hs = h3_pT(p, T);rhos = 1 / v3_ph(p, hs);Out = fromSIunit_Cp(Cp3_rhoT(rhos, T));case 4Out = NaN;case 5Out = fromSIunit_Cp(Cp5_pT(p, T));otherwiseOut = NaN;endfunction Cv1_pT = Cv1_pT(p, T)%Release on the IAPWS Industrial formulation 1997 for the ThermodynamicProperties <strong>of</strong> Water and Steam, September 1997%5 Equations for Region 1, Section. 5.1 Basic Equation%Eqution 7, Table 3, Page 6I1 = [0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 4, 4, 4, 5, 8, 8, 21, 23, 29,30, 31, 32];J1 = [-2, -1, 0, 1, 2, 3, 4, 5, -9, -7, -1, 0, 1, 3, -3, 0, 1, 3, 17, -4, 0, 6, -5, -2, 10, -8, -11,-6, -29, -31, -38, -39, -40, -41];n1 = [0.14632971213167, -0.84548187169114, -3.756360367204, 3.3855169168385,-0.95791963387872, 0.15772038513228, -0.016616417199501, 8.1214629983568E-04, 2.8319080123804E-04, -6.0706301565874E-04, -0.018990068218419, -0.032529748770505, -0.021841717175414, -5.283835796993E-05, -4.7184321073267E-04, -3.0001780793026E-04, 4.7661393906987E-05, -4.4141845330846E-06, -7.2694996297594E-16, -3.1679644845054E-05, -2.8270797985312E-06, -8.5205128120103E-10, -2.2425281908E-06, -6.5171222895601E-07, -1.4341729937924E-13, -4.0516996860117E-07, -1.2734301741641E-09, -1.7424871230634E-10, -6.8762131295531E-19,1.4478307828521E-20, 2.6335781662795E-23, -1.1947622640071E-23,1.8228094581404E-24, -9.3537087292458E-26];R = 0.461526; Pi = p / 16.53;tau = 1386 / T;gamma_der_pi = 0;gamma_der_pipi = 0;gamma_der_pitau = 0;
- Page 1 and 2:
ANALYSIS OF WATER INJECTION INTO HI
- Page 3 and 4:
Thesis CertificateThe Graduate Coll
- Page 5 and 6:
ชื่อ : นายปรม
- Page 7 and 8:
TABLE OF CONTENTSPageAbstract (in E
- Page 9 and 10:
LIST OF TABLESTablePage3-1 Solution
- Page 11 and 12:
LIST OF FIGURES (CONTINUED)FigurePa
- Page 13 and 14:
LIST OF ABBREVIATIONS, SYMBOLS AND
- Page 15 and 16:
CHAPTER 1INTRODUCTION1.1 Background
- Page 17 and 18:
31.6.6 Water injected is assumed to
- Page 19 and 20:
6for these working fluid models can
- Page 22 and 23:
CHAPTER 3METHODOLOGY FOR ANALYSIS O
- Page 24 and 25:
11perform the necessary calculation
- Page 26 and 27:
13h b P T w−0.2 0.8 −0.55 0.8=
- Page 28 and 29:
15the procedure required Nitrogen(
- Page 30 and 31:
171N22 NEq.3-461 1O+2N2 NOEq.3-472
- Page 32 and 33:
19∂y1 cy ∂y1 cy ∂y 1 c ∂y 1
- Page 34 and 35:
211/2( cy )∂y ∂∂c ∂y∂T
- Page 36 and 37:
23[ A][ ∂y/ ∂ P] + [ ∂f / ∂
- Page 38 and 39:
25⎛ • • •ln ln ⎞⎛⎞( )
- Page 40 and 41:
27( θ= −π)θ>θ bθ>θ W( θ=π
- Page 42 and 43:
CHAPTER 4RESULTS AND DISCUSSIONThis
- Page 44 and 45:
31FIGURE 4-1 Comparison of an actua
- Page 46 and 47:
33FIGURE 4-4 Schematic of the port
- Page 48 and 49:
35Temperature (K)250023002100190017
- Page 50 and 51:
37cylinder temperature. So, a more
- Page 52 and 53:
39Thermal efficiency (%)43424140393
- Page 54 and 55:
41Theoretically, the high useful co
- Page 56 and 57:
43Thermal efficiency (%)46454443424
- Page 58 and 59:
45injection-fuel ratio increases gr
- Page 60 and 61:
47Thermal efficiency (%)44434241403
- Page 62 and 63:
49When considering relative tempera
- Page 64 and 65:
51Temperature (K)240021001800150012
- Page 66 and 67:
53When considering relative tempera
- Page 68 and 69:
55Temperature (K)220019001600130010
- Page 70 and 71:
CHAPTER 5CONCLUSIONS AND SUGGESTION
- Page 72 and 73: REFERENCES1. Ricardo, H.R. The High
- Page 74: APPENDIX ADerivative equations of i
- Page 78 and 79: 64TABLE A-3 Curve fit coefficients
- Page 80 and 81: 66TABLE A-5 Curve fit coefficients
- Page 82 and 83: 68The following is the derivative v
- Page 84 and 85: 70From the definition of entropyh =
- Page 86 and 87: 72• ⎛ ⎞ • ⎛Tb∂u •b∂
- Page 88 and 89: 74• ⎛1 ⎞ •∂u ⎛∂ ∂
- Page 90 and 91: 76⎛ • • •ln ln ⎞⎛⎞( )
- Page 92 and 93: 78The following is the derivative v
- Page 94 and 95: 80From the definition of entropy h
- Page 96 and 97: 82• ⎛ Tb u ⎞ • • ⎛bmTb
- Page 98 and 99: 84• ⎛1 ⎞ • •∂u ⎛ ⎞
- Page 100 and 101: 86• • • ⎛ u ⎞bmb( hW ub)
- Page 102 and 103: 88Appendix D MATLAB program scripts
- Page 104 and 105: 90[thetawater,pTbWQlHl2]=ode45('Rat
- Page 106 and 107: 92-0.69353550E-14 -0.14245228E+05 0
- Page 108 and 109: 940.21*(1-phi) 0 0]';dcdT=0;else %
- Page 110 and 111: 96dfdp=zeros(4,1);dYdT=zeros(11,1);
- Page 112 and 113: 98dcdT(1)=-dKdT(1)*sqrt(patm)/K(1)^
- Page 114 and 115: 100Iter=Iter+1;[hb,u,v,s,Y,cp,dlvlT
- Page 116 and 117: 102yprime(2)=-Const1/cpb/x*Qconvb+v
- Page 118 and 119: 104masswaterin=massfuel*mwpmf;% Tot
- Page 120 and 121: 106savefile = 'Volume.mat';RTWV=RTW
- Page 124 and 125: 110gamma_der_tautau = 0;for i = 1 :