analysis of water injection into high-temperature mixture of ...

More documents

Recommendations

Info

92-0.69353550E-14 -0.14245228E+05 0.63479156E+010.31001901E+01 0.51119464E-03 0.52644210E-07 -0.34909973E-10 ...0.36945345E-14 -0.87738042E+03 -0.19629421E+010.25000000E+01 0.00000000E+00 0.00000000E+00 0.00000000E+00 ...0.00000000E+00 0.25471627E+05 -0.46011763E+000.25420596E+01 -0.27550619E-04 -0.31028033E-08 0.45510674E-11 ...-0.43680515E-15 0.29230803E+05 0.49203080E+010.29106427E+01 0.95931650E-03 -0.19441702E-06 0.13756646E-10 ...0.14224542E-15 0.39353815E+04 0.54423445E+010.31890000E+01 0.13382281E-02 -0.52899318E-06 0.95919332E-10 ...-0.64847932E-14 0.98283290E+04 0.67458126E+012.41594293E+00 1.74890600E-04 -1.19023667E-07 3.02262387E-11 ...-2.03609790E-15 5.61337748E+04 4.64960986E+00];EndD-3. M-file for fueldata.mfunction [alpha,beta,gamma,delta,Afuel]=fueldata(fuel);switch fuelcase 'gasoline' % Fergusonalpha=7; beta=17; gamma=0; delta=0;Afuel=[4.0652 6.0977E-02 -1.8801E-05 0 0 -3.5880E+04 15.45];case 'diesel' % Fergusonalpha=14.4; beta=24.9; gamma=0; delta=0;Afuel=[7.9710 1.1954E-01 -3.6858E-05 0 0 -1.9385E+04 -1.7879];case 'methane' % Fergusonalpha=1; beta=4; gamma=0; delta=0;Afuel=[1.971324 7.871586E-03 -1.048592E-06 0 0 -9.930422E+03 8.873728];case 'methanol' % Fergusonalpha=1; beta=4; gamma=1; delta=0;Afuel=[1.779819 1.262503E-02 -3.624890E-06 0 0 -2.525420E+04 1.50884E+01];case 'nitromethane' % Fergusonalpha=1; beta=3; gamma=2; delta=1;Afuel=[1.412633 2.087101E-02 -8.142134E-06 0 0 -1.026351E+041.917126E+01];case 'benzene' % Fergusonalpha=6; beta=6; gamma=0; delta=0;Afuel=[-2.545087 4.79554E-02 -2.030765E-05 0 0 8.782234E+03 3.348825E+01];EndD-4. M-file for enginedata.m% ***** engine geometry **********************************************b=0.1; % engine bore (m)stroke=0.08; % engine stroke (m)eps=0.25; % half stroke to rod ratio, s/2lr=10; % compression ratioVtdc=pi/4*b^2*stroke/(r-1); % volume at TDC
93Vbdc=pi/4*b^2*stroke+Vtdc; % volume at BDC% ***** engine thermofluids parameters *******************************Xv=1.0; %Ratio of mole water injection-fuel 0.1-1.0Tv=400; %Temperature of vapor injection > 373 (K)Pr=1.1;% %Pressure Vaporr injection per Pressure in cyclinderCblowby=0.8; % piston blowby constant (s^-1)f=0.1; % residual fractionfueltype='gasoline';airscheme='GMcB';phi=1.2; % equivalence ratiothetas=-35*pi/180; % start of burningthetab=60*pi/180; % burn duration anglethetaWa=10*pi/180; % Vapor injection duration anglethetaperiod=5*pi/180; %Period angle Before Vapor injectionRPM=2000;omega=RPM*pi/30; % engine speed in rad/sheattransferlaw='constant'; % 'constant', or 'Woschni'hcu=500; % unburned zone heat transfer coefficient/weightinghcb=500; % burned zone heat transfer coefficient/weightingTw=420; % engine surface temperature% ***** initial conditions *******************************************p1=100e3;T1=350;theta1=-pi;V1=Vbdc;[h1,u1,v1,s1,Y1,cp1,dlvlT1,dlvlp1]=farg(p1,T1,phi,f,fueltype,airscheme);mass1=Vbdc/v1;U1=u1*mass1;D-5. M-file for farg.mfunction [h,u,v,s,Y,cp,dlvlT,dlvlp]=fargw(Xv,p,T,phi,f,fueltype,airscheme);[alpha,beta,gamma,delta,Afuel]=fueldata(fueltype);switch airschemecase 'GMcB'A=airdata('GMcB_low');case 'Chemkin'A=airdata('Chemkin_low');end%X=0.125; % Mole fractions of fuel-waterRu=8314.34; % J/kmolKtable=[-1 1 0 0 1 -1]';M=[44.01 18.02 28.008 32.000 28.01 2.018]'; % kg/kmolMinMol=1e-25;dlvlT=1; dlvlp=-1;eps=0.210/(alpha+0.25*beta-0.5*gamma);if phi
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 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 122 and 123: 108p=pTarray(:,1);T=pTarray(:,2);hc
Page 124 and 125: 110gamma_der_tautau = 0;for i = 1 :
show all

analysis of water injection into high-temperature mixture of ...

Create successful ePaper yourself

Delete template?

Save as template?