Thông tin tài liệu
bộ giáo dục đào tạo trờng đại học bách khoa hà nội - luận văn thạc sĩ khoa học nghiên cứu mô hình hoá trình tạo nhiên liệu giàu hydrô để giảm phát thải độc hại động đốt ngành : kỹ thuật động nhiệt m số:23.04.3898 phạm ngọc anh Ngời hớng dẫn khoa học: PTS hoàng đình long Hà Nội 2009 LI CAM OAN Tụi xin cam oan õy l ti nghiờn cu ca riờng tụi Cỏc s liu kt qu nờu lun l trung thc v cha tng c cụng b cỏc cụng trỡnh no khỏc! H ni , thỏng 11 nm 2009 Tỏc gi Phm Ngc Anh DANH MC CC Kí HIU V CH VIT TT CO : Mụnụxớt cỏcbon CO2 : Cỏcboniụxớt CmHn v H-C : Hyụcỏcbon NOX : Cỏc loi ụxớtnit SO2 : Sunfua dioxit P-M : Cht thi dng ht N2 : Nit O2 : ễxy CR : Common Rail CNG : Nhiờn liu khớ thiờn nhiờn nộn LNG : Khớ thiờn nhiờn hoỏ lng LPG : Khớ du m hoỏ lng H2 : Nhiờn liu khớ hydrụ T0 : Nhit u vo Tc : Nhit cui quỏ trỡnh nộn k : H s on nhit Vo/Vc : T s nộn A/F : T s khụng khớ trrờn nhiờn liu : Hiu sut nhit lý thuyt pa : p sut cui quỏ trỡnh np po : p sut du quỏ trỡnh np pr : p sut khớ thi Tr : Nhit khớ sút Q : Lng nhit cn truyn cho kột G : Lu lng hn hp khớ giu hydrụ DANH MC BNG Cỏc bng Trang E j Bng 3-1: Cỏc hng s ng hc phn ng kj = k j e RT 59 Bng 3-2: Cỏc hng s cõn bng 60 Bng 3-3: Cỏc hng s hp th 60 DANH MC HèNH V V TH Cỏc hỡnh v v th Trang Hỡnh 1-1: B x lý xỳc tỏc ba ng 21 Hỡnh 1-2: Lc mui than khớ thi ng c Diesel 26 Hỡnh 2-1: Khớ thi ca ng c hydrụ 40 Hỡnh 2-2: Khớ thi ca ng c xng Hỡnh 2-3: So sỏnh th tớch bung chỏy v hm lng nng lng 41 gia ng c xng v ng c s dng nhiờn liu hydro Hỡnh 3-1: S 2D ca lũ phn ng xỳc tỏc nghiờn cu 68 Hỡnh 3-2: ng c diesel D243 phc v th nghim 74 Hỡnh 3-3: S nguyờn lý b xỳc tỏc 76 Hỡnh 3-4: Bn v thit k b xỳc tỏc 80 Hỡnh 3-5: Kớch thc mt ngang ca ng lm mỏt khớ giu hydrụ 81 Hỡnh 3-6: Kớch thc cỏnh tn nhit 81 Hỡnh 3-7: Kớch thc chung ca kột tn nhit 82 th 1: Phn trm sn phm F/W = 2/0 V=5000l/h 68 th 2: Phn trm sn phm F/W = 2/2 V=5000l/h th 3: Phn trm sn phm F/W = 2/4 V=5000l/h 69 th 4: Phn trm sn phm F/W = 2/0 V=10000l/h th 5: Phn trm sn phm F/W = 2/2 V=10000l/h 70 th 6: Phn trm sn phm F/W = 2/4 V=10000l/h th 7: Phn trm sn phm F/W = 2/0 V=20000l/h th 8: Phn trm sn phm F/W = 2/2 V=20000l/h th 9: Phn trm sn phm F/W = 2/4 V=20000l/h th 10: Phn trm sn phm Vf = 0,5l/h th 11: Phn trm sn phm Vf = 1l/h th 12: Phn trm sn phm Vf = 2l/h th 13: Phn trm hydrụ sinh thay i chiu di b xỳc tỏc 70 71 72 73 74 MC LC LI CAM OAN DANH MC CC Kí HIU V CH VIT TT DANH MC BNG DANH MC HèNH V V TH LI M U CHNG 1: GII THIU CHUNG 1.1 S phỏt thi ca ng c t v phng phỏp gim phỏt thi 1.1.1 S phỏt thi c hi ca ng c t 1.1.1.1 Cỏc thnh phn phỏt thi c hi ca ng c t 1.1.1.2 Tỡnh hỡnh ụ nhim mụi trng khớ thi ng c gõy Vit Nam 1.1.1.3 Cỏc ngun phỏt thi ca ng c t 1.1.2 Cỏc bin phỏp gim phỏt thi cỏc thnh phn c hi 10 1.1.2.1 Gim thnh phn c hi khớ thi nh ti u hoỏ thit k ng c 10 1.1.2.2 Gim thnh phn c hi bng bin phỏp x lý khớ thi 14 1.1.2.3 Thay i cỏch s dng nhiờn liu v dựng nhiờn liu thay th 21 1.2.3 Mc ớch v phm vi nghiờn cu ca ti 24 1.2.3.1 Mc ớch ca ti 24 1.2.3.2 Phm vi nghiờn cu ca ti 24 1.2.3.3 í ngha khoa hc ca ti 24 CHNG 2: LI CH CA NHIấN LIU GIU HYDRO V CC PHNG PHP TO RA T NHIấN LIU TRUYN THNG 26 2.1 Vn chung v s dng hydro v nhiờn liu giu H2 trờn th gii 26 2.1.1 Mt s c im liờn quan n nhiờn liu hydro s dng trờn ng c 31 2.1.1.1 Phm vi chỏy rng 31 2.1.1.2 Khong dp tt nh 31 2.1.1.3 Nhit t chỏy cao 32 2.1.1.4 Tc chỏy cao 32 2.1.1.5 khuờch tỏn cao 32 2.1.1.6 Mt thp 33 2.2 Phng phỏp sn xut hydro v nhiờn liu giu hydro 38 2.2.1 Phn ng nhit húa nhiờn liu hydrocacbon vi hi nc (phn ng thu nhit ) 39 2.2.2 Phn ng oxy húa nhiờn liu khụng hon ton (phn ng ta nhit): 40 2.2.3 Phn ng nhit húa hydrocacbon vi khớ cacbonic (phn ng thu nhit) 41 2.2.4 Bin i nhit húa nhiờn liu hydrocacbon vi ng thi hi nc, oxy v cacbonic (quỏ trỡnh phn ng trung tớnh v nhit ) 41 CHNG 44 NGHIấN CU TN DNG NHIT KH THI TO NHIấN LIU GIU HYDRễ 44 3.1 Nhit khớ thi v kh nng tn dng 44 3.1.1 Nhit khớ thi ngoi 44 3.1.2 Tớnh toỏn cỏc quỏ trỡnh cụng tỏc: 45 3.2 Mụ hỡnh toỏn quỏ trỡnh to nhiờn liu giu hydro t nhiờn liu truyn thng 49 3.2.1 Gii thiu chung 49 3.2.2 Cỏc phn ng xỳc tỏc bin i nhit húa ca diesel 50 3.2.3 Tc cỏc phn ng bin i nhit húa ca diesel 53 3.2.4 Cỏc phng trỡnh c bn 55 3.2.4.1 Tớnh cht ca khớ 59 3.2.4.2 H s trao i nhit v trao i cht 59 3.3 Kt qu nghiờn cu 61 3.4 Thit k ch to h thng xỳc tỏc 68 3.4.1 Cỏc thụng s ca ng c D243 68 3.4.2 Thit k h thng xỳc tỏc cho ng c diesel D243 69 3.4.2.1 Yờu cu i vi b xỳc tỏc 69 3.4.2.2 Thit k h thng xỳc tỏc 69 3.4.2.3 Tớnh toỏn, thit k kột lm mỏt hn hp khớ giu hydrụ 75 CHNG IV KT LUN V HNG PHT TRIN CA TI 80 4.1 Kt lun 80 4.2 Hng phỏt trin ca ti 81 TI LIU THAM KHO 1 LI M U Trong thi k cụng nghip húa v hin i húa t nc cựng vi s i lờn ca cỏc ngnh kinh t, ng c t ngy cng gi mt vai trũ quan trng v c ng dng nhiu trờn cỏc lnh vc khỏc nh nụng nghip, cụng nghip, giao thụng ti cng nh nhiu ngnh cú kinh t khỏc S gia tng nhanh v s lng ng c t s dng ó v ang gõy ụ nhim mụi trng trm trng khớ thi c hi Chớnh vỡ vy t cho cỏc nh nghiờn cu l lm th no gim c ụ nhim khớ thi Mt nhng hng nghiờn cu hin l dựng nhiờn liu thay th gõy phỏt thi c hi ớt, vớ d nhiờn liu hydrụ v nhiờn liu giu hydụ Nhiờn liu hyrụ l loi nhiờn liu t lõu ó c cỏc nh nghiờn cu cõn nhc s dng cho ng c t Nhiờn liu ny cú u im l chỏy nhanh, tr s c tan cao nờn cho phộp ng c cú th lm vic tc rt cao, nh ú m d dng tng cụng sut ng c Thờm na, khớ thi ca ng c hydro rt sch, gii hn thnh phn hn hp m bo kh nng chỏy tt rt rng nờn ng c cú th lm vic vi hn hp rt loóng, gúp phn lm tng tớnh kinh t s dng ng c Tuy nhiờn, nhiờn liu hydro cú nhc im so vi nhiờn liu xng v diesel l nhit tr mole rt thp nờn nu khụng thay i kt cu ng c chuyn t ng c chy xng hoc diesel sang ng c chy hon ton bng hydro thỡ cụng sut ng c s b gim nhiu Thờm na, vic sn xut, chuyn v tớch tr bo qun nhiờn liu hydro thay th hon ton xng hoc diesel khỏ khú khn v tn kộm nhiờn liu cú t trng rt thp Chớnh vỡ vy, nhiu nh nghiờn cu quan tõm n vic s dng hydro nh mt thnh phn ph gia cho nhiờn liu xng hoc diesel Vi phng phỏp ny, ch cn cp mt t l nh hydrụ vo ng c hũa trn vi nhiờn liu chớnh l xng hoc diesel Nh c tớnh chỏy nhanh, hydro s giỳp t chỏy hn hp nhiờn liu chớnh tt hn nờn gim c thnh phn c hi khớ thi Vn t l cn cú ngun cung cp hydro n nh v tin li trờn xe ng c hot ng liờn tc vi hn hp nhiờn liu núi trờn Trong ti tt nghip ny Em a hng nghiờn cu gim thnh phn c hi khớ thi ng c diesel bng cỏch to v cung cp nhiờn liu giu hydro cho ng c bng phng phỏp phn ng xỳc tỏc mt phn nhiờn liu diesel nh chớnh khớ thi ng c i tng nghiờn cu l ng c D243 (Nga) Brarut sn xut vi mc ớch l nghiờn cu mụ hỡnh hoỏ quỏ trỡnh phn ng xỳc tỏc nhiờn liu vi khớ thi to nhiờn liu giu hydrụ v t ú thit k ti u b xỳc tỏc Ni dung tt nghip gm hai phn, th nht l nghiờn cu li ớch ca nhiờn liu giu hydro v cỏc phng phỏp to nhiờn liu ny t nhiờn liu truyn thng Th hai l nghiờn cu thit k b xỳc tỏc tn dng nhit khớ thi to nhiờn liu giu hydro Em xin chõn thnh cỏm n TS Hong ỡnh Long ó giỳp v ch bo cho em rt nhiu Cỏm n cỏc thy cụ b mụn ng c t trong, phũng thớ nghim ng c t ó to mi iu kin v thi gian cng nh trang thit b thớ nghim em cú th hon thnh bn lun ny Hc viờn cao hc KS Phm Ngc Anh & (cpctd(ig)*100.,ig=1,igmax-1),(cpct(ig)*100.,ig=1,igmax-1) end if if (it eq 18000 and ir eq and mod(iz,4) eq 0) then ctotal=c(1,ir,iz)+c(2,ir,iz)+c(3,ir,iz)+ & c(4,ir,iz)+c(5,ir,iz)+c(6,ir,iz)+c(7,ir,iz) ig=1,igmax cpct(ig)=c(ig,ir,iz)/ctotal end do ig=1,igmax cpctd(ig)=cpct(ig)/(1.-cpct(4)) end HCcvt=(cin(1)-cpct(1))/cin(1)*100 write(7,62)iz*delz,Tg(iz),T(ir,iz),HCcvt, & (cpctd(ig)*100.,ig=1,igmax-1),(cpct(ig)*100.,ig=1,igmax-1) end if if (it eq 30000 and ir eq and mod(iz,4) eq 0) then ctotal=c(1,ir,iz)+c(2,ir,iz)+c(3,ir,iz)+ & c(4,ir,iz)+c(5,ir,iz)+c(6,ir,iz)+c(7,ir,iz) ig=1,igmax cpct(ig)=c(ig,ir,iz)/ctotal end do ig=1,igmax cpctd(ig)=cpct(ig)/(1.-cpct(4)) end HCcvt=(cin(1)-cpct(1))/cin(1)*100 write(8,62)iz*delz,Tg(iz),T(ir,iz),HCcvt, & (cpctd(ig)*100.,ig=1,igmax-1),(cpct(ig)*100.,ig=1,igmax-1) end if C**** Tinh he so truyen nhiet qua if (ir eq irmax-1) then hiner=0.27*kga/2./(irmax*delr)*(2*(irmax*delr)*u*rog/ !0.027 & void/mvis)**0.8*(mvis*cpg/kga)**0.33*4.121*(0.001*Theat)**5.5 Raa=dout**3*9.81/Theat*100.*0.7/(21.e-6)**2 houter=0.031/dout*(0.6+0.387*Raa**(1/6.)/ !tinh cho doi luu tu nhien & (1.+(0.559/0.7)**(9/16.))**(8/27.))**2 htt(iz)=1./(1./hiner+bw/50.7)*heat end if u=uin*Rgmc*T(0,iz)/Rgmin/Tgin thay doi theo z 300 continue !end loop z ********************************************** 200 continue !end loop r ********************************************** iz=-1,izmax+1 ig=1,igmax c(ig,0,iz)=c(ig,1,iz) end end C**** Print data C -Data trung binh doc theo z Do iz=0,izmax if (it eq itmax and mod(iz,4) eq 0) then Tgtb=0 ! toc do ir=1,irmax Tgtb=Tgtb+s(ir)*T(ir,iz)/(irmax**2*1.) end do ig=1,igmax cg(ig)=0 ir=1,irmax cg(ig)=cg(ig)+s(ir)*c(ig,ir,iz)/(irmax**2*1.) end end ctotal=cg(1)+cg(2)+cg(3)+cg(4)+cg(5)+cg(6)+cg(7) ig=1,igmax cpct(ig)=cg(ig)/ctotal end do ig=1,igmax cpctd(ig)=cpct(ig)/(1.-cpct(4)) end HCcvt=(cin(1)-cpct(1))/cin(1)*100 write(41,62)iz*delz,Tgtb,T(0,iz),HCcvt, & (cpctd(ig)*100.,ig=1,igmax-1),(cpct(ig)*100.,ig=1,igmax-1) end if if (it eq 3000 and mod(iz,4) eq 0) then Tgtb=0 ir=1,irmax Tgtb=Tgtb+s(ir)*T(ir,iz)/(irmax**2*1.) end do ig=1,igmax cg(ig)=0 ir=1,irmax cg(ig)=cg(ig)+s(ir)*c(ig,ir,iz)/(irmax**2*1.) end end ctotal=cg(1)+cg(2)+cg(3)+cg(4)+cg(5)+cg(6)+cg(7) ig=1,igmax cpct(ig)=cg(ig)/ctotal end do ig=1,igmax cpctd(ig)=cpct(ig)/(1.-cpct(4)) end HCcvt=(cin(1)-cpct(1))/cin(1)*100 write(51,62)iz*delz,Tgtb,T(0,iz),HCcvt, & (cpctd(ig)*100.,ig=1,igmax-1),(cpct(ig)*100.,ig=1,igmax-1) end if if (it eq 6000 and mod(iz,4) eq 0) then Tgtb=0 ir=1,irmax Tgtb=Tgtb+s(ir)*T(ir,iz)/(irmax**2*1.) end do ig=1,igmax cg(ig)=0 ir=1,irmax cg(ig)=cg(ig)+s(ir)*c(ig,ir,iz)/(irmax**2*1.) end end ctotal=cg(1)+cg(2)+cg(3)+cg(4)+cg(5)+cg(6)+cg(7) ig=1,igmax cpct(ig)=cg(ig)/ctotal end do ig=1,igmax cpctd(ig)=cpct(ig)/(1.-cpct(4)) end HCcvt=(cin(1)-cpct(1))/cin(1)*100 write(61,62)iz*delz,Tgtb,T(0,iz),HCcvt, & (cpctd(ig)*100.,ig=1,igmax-1),(cpct(ig)*100.,ig=1,igmax-1) end if if (it eq 18000 and mod(iz,4) eq 0) then Tgtb=0 ir=1,irmax Tgtb=Tgtb+s(ir)*T(ir,iz)/(irmax**2*1.) end do ig=1,igmax cg(ig)=0 ir=1,irmax cg(ig)=cg(ig)+s(ir)*c(ig,ir,iz)/(irmax**2*1.) end end ctotal=cg(1)+cg(2)+cg(3)+cg(4)+cg(5)+cg(6)+cg(7) ig=1,igmax cpct(ig)=cg(ig)/ctotal end do ig=1,igmax cpctd(ig)=cpct(ig)/(1.-cpct(4)) end HCcvt=(cin(1)-cpct(1))/cin(1)*100 write(71,62)iz*delz,Tgtb,T(0,iz),HCcvt, & (cpctd(ig)*100.,ig=1,igmax-1),(cpct(ig)*100.,ig=1,igmax-1) end if if (it eq 30000 and mod(iz,4) eq 0) then Tgtb=0 ir=1,irmax Tgtb=Tgtb+s(ir)*T(ir,iz)/(irmax**2*1.) end do ig=1,igmax cg(ig)=0 ir=1,irmax cg(ig)=cg(ig)+s(ir)*c(ig,ir,iz)/(irmax**2*1.) end end ctotal=cg(1)+cg(2)+cg(3)+cg(4)+cg(5)+cg(6)+cg(7) ig=1,igmax cpct(ig)=cg(ig)/ctotal end do ig=1,igmax cpctd(ig)=cpct(ig)/(1.-cpct(4)) end HCcvt=(cin(1)-cpct(1))/cin(1)*100 write(81,62)iz*delz,Tgtb,T(0,iz),HCcvt, & (cpctd(ig)*100.,ig=1,igmax-1),(cpct(ig)*100.,ig=1,igmax-1) end if end C -Data theo z tai cac r khac va thoi gian khac iz=0,izmax if (it eq itmax and mod(iz,4) eq 0) then ir=0,irmax ctr(ir)=c(1,ir,iz)+c(2,ir,iz)+c(3,ir,iz)+c(4,ir,iz)+ & c(5,ir,iz)+c(6,ir,iz)+c(7,ir,iz) cir1(ir)=c(1,ir,iz)/ctr(ir) cir2(ir)=c(2,ir,iz)/ctr(ir) cir3(ir)=c(3,ir,iz)/ctr(ir) cir4(ir)=c(4,ir,iz)/ctr(ir) cir5(ir)=c(5,ir,iz)/ctr(ir) cir6(ir)=c(6,ir,iz)/ctr(ir) cird1(ir)=c(1,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird2(ir)=c(2,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird3(ir)=c(3,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird4(ir)=c(4,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird5(ir)=c(5,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird6(ir)=c(6,ir,iz)/(ctr(ir)-c(4,ir,iz)) end write(42,72)iz*delz,(T(ir,iz),ir=0,irmax,jr), & (cird1(ir)*100.,ir=0,irmax,jr), & (cird2(ir)*100.,ir=0,irmax,jr), & (cird3(ir)*100.,ir=0,irmax,jr), & (cird4(ir)*100.,ir=0,irmax,jr), & (cird5(ir)*100.,ir=0,irmax,jr), & (cird6(ir)*100.,ir=0,irmax,jr), & (cir1(ir)*100.,ir=0,irmax,jr), & (cir2(ir)*100.,ir=0,irmax,jr), & (cir3(ir)*100.,ir=0,irmax,jr), & (cir4(ir)*100.,ir=0,irmax,jr), & (cir5(ir)*100.,ir=0,irmax,jr), & (cir6(ir)*100.,ir=0,irmax,jr) end if if (it eq 3000 and mod(iz,4) eq 0) then ir=0,irmax ctr(ir)=c(1,ir,iz)+c(2,ir,iz)+c(3,ir,iz)+c(4,ir,iz)+ & c(5,ir,iz)+c(6,ir,iz)+c(7,ir,iz) cir1(ir)=c(1,ir,iz)/ctr(ir) cir2(ir)=c(2,ir,iz)/ctr(ir) cir3(ir)=c(3,ir,iz)/ctr(ir) cir4(ir)=c(4,ir,iz)/ctr(ir) cir5(ir)=c(5,ir,iz)/ctr(ir) cir6(ir)=c(6,ir,iz)/ctr(ir) cird1(ir)=c(1,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird2(ir)=c(2,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird3(ir)=c(3,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird4(ir)=c(4,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird5(ir)=c(5,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird6(ir)=c(6,ir,iz)/(ctr(ir)-c(4,ir,iz)) end write(52,72)iz*delz,(T(ir,iz),ir=0,irmax,jr), & (cird1(ir)*100.,ir=0,irmax,jr), & (cird2(ir)*100.,ir=0,irmax,jr), & (cird3(ir)*100.,ir=0,irmax,jr), & (cird4(ir)*100.,ir=0,irmax,jr), & (cird5(ir)*100.,ir=0,irmax,jr), & (cird6(ir)*100.,ir=0,irmax,jr), & (cir1(ir)*100.,ir=0,irmax,jr), & (cir2(ir)*100.,ir=0,irmax,jr), & (cir3(ir)*100.,ir=0,irmax,jr), & (cir4(ir)*100.,ir=0,irmax,jr), & (cir5(ir)*100.,ir=0,irmax,jr), & (cir6(ir)*100.,ir=0,irmax,jr) end if if (it eq 6000 and mod(iz,4) eq 0) then ir=0,irmax ctr(ir)=c(1,ir,iz)+c(2,ir,iz)+c(3,ir,iz)+c(4,ir,iz)+ & c(5,ir,iz)+c(6,ir,iz)+c(7,ir,iz) cir1(ir)=c(1,ir,iz)/ctr(ir) cir2(ir)=c(2,ir,iz)/ctr(ir) cir3(ir)=c(3,ir,iz)/ctr(ir) cir4(ir)=c(4,ir,iz)/ctr(ir) cir5(ir)=c(5,ir,iz)/ctr(ir) cir6(ir)=c(6,ir,iz)/ctr(ir) cird1(ir)=c(1,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird2(ir)=c(2,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird3(ir)=c(3,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird4(ir)=c(4,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird5(ir)=c(5,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird6(ir)=c(6,ir,iz)/(ctr(ir)-c(4,ir,iz)) end write(62,72)iz*delz,(T(ir,iz),ir=0,irmax,jr), & (cird1(ir)*100.,ir=0,irmax,jr), & (cird2(ir)*100.,ir=0,irmax,jr), & (cird3(ir)*100.,ir=0,irmax,jr), & (cird4(ir)*100.,ir=0,irmax,jr), & (cird5(ir)*100.,ir=0,irmax,jr), & (cird6(ir)*100.,ir=0,irmax,jr), & (cir1(ir)*100.,ir=0,irmax,jr), & (cir2(ir)*100.,ir=0,irmax,jr), & (cir3(ir)*100.,ir=0,irmax,jr), & (cir4(ir)*100.,ir=0,irmax,jr), & (cir5(ir)*100.,ir=0,irmax,jr), & (cir6(ir)*100.,ir=0,irmax,jr) end if if (it eq 18000 and mod(iz,4) eq 0) then ir=0,irmax ctr(ir)=c(1,ir,iz)+c(2,ir,iz)+c(3,ir,iz)+c(4,ir,iz)+ & c(5,ir,iz)+c(6,ir,iz)+c(7,ir,iz) cir1(ir)=c(1,ir,iz)/ctr(ir) cir2(ir)=c(2,ir,iz)/ctr(ir) cir3(ir)=c(3,ir,iz)/ctr(ir) cir4(ir)=c(4,ir,iz)/ctr(ir) cir5(ir)=c(5,ir,iz)/ctr(ir) cir6(ir)=c(6,ir,iz)/ctr(ir) cird1(ir)=c(1,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird2(ir)=c(2,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird3(ir)=c(3,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird4(ir)=c(4,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird5(ir)=c(5,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird6(ir)=c(6,ir,iz)/(ctr(ir)-c(4,ir,iz)) end write(72,72)iz*delz,(T(ir,iz),ir=0,irmax,jr), & (cird1(ir)*100.,ir=0,irmax,jr), & (cird2(ir)*100.,ir=0,irmax,jr), & (cird3(ir)*100.,ir=0,irmax,jr), & (cird4(ir)*100.,ir=0,irmax,jr), & (cird5(ir)*100.,ir=0,irmax,jr), & (cird6(ir)*100.,ir=0,irmax,jr), & (cir1(ir)*100.,ir=0,irmax,jr), & (cir2(ir)*100.,ir=0,irmax,jr), & (cir3(ir)*100.,ir=0,irmax,jr), & (cir4(ir)*100.,ir=0,irmax,jr), & (cir5(ir)*100.,ir=0,irmax,jr), & (cir6(ir)*100.,ir=0,irmax,jr) end if if (it eq 30000 and mod(iz,4) eq 0) then ir=0,irmax ctr(ir)=c(1,ir,iz)+c(2,ir,iz)+c(3,ir,iz)+c(4,ir,iz)+ & c(5,ir,iz)+c(6,ir,iz)+c(7,ir,iz) cir1(ir)=c(1,ir,iz)/ctr(ir) cir2(ir)=c(2,ir,iz)/ctr(ir) cir3(ir)=c(3,ir,iz)/ctr(ir) cir4(ir)=c(4,ir,iz)/ctr(ir) cir5(ir)=c(5,ir,iz)/ctr(ir) cir6(ir)=c(6,ir,iz)/ctr(ir) cird1(ir)=c(1,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird2(ir)=c(2,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird3(ir)=c(3,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird4(ir)=c(4,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird5(ir)=c(5,ir,iz)/(ctr(ir)-c(4,ir,iz)) cird6(ir)=c(6,ir,iz)/(ctr(ir)-c(4,ir,iz)) end write(82,72)iz*delz,(T(ir,iz),ir=0,irmax,jr), & (cird1(ir)*100.,ir=0,irmax,jr), & (cird2(ir)*100.,ir=0,irmax,jr), & (cird3(ir)*100.,ir=0,irmax,jr), & (cird4(ir)*100.,ir=0,irmax,jr), & (cird5(ir)*100.,ir=0,irmax,jr), & (cird6(ir)*100.,ir=0,irmax,jr), & (cir1(ir)*100.,ir=0,irmax,jr), & (cir2(ir)*100.,ir=0,irmax,jr), & (cir3(ir)*100.,ir=0,irmax,jr), & (cir4(ir)*100.,ir=0,irmax,jr), & (cir5(ir)*100.,ir=0,irmax,jr), & end if end (cir6(ir)*100.,ir=0,irmax,jr) C -Data theo thoi gian if (mod(it,500) eq 0) then ig=1,igmax cg(ig)=0 ir=1,irmax cg(ig)=cg(ig)+s(ir)*c(ig,ir,izmax)/(irmax**2*1.) end end ctotal=cg(1)+cg(2)+cg(3)+cg(4)+cg(5)+cg(6)+cg(7) ig=1,igmax cpct(ig)=cg(ig)/ctotal end if (cpct(1) gt cin(1)-1.e-3) cpct1=cin(1)-1.e-3 !de cho an toan c1cvt=(cin(1)-cpct(1))/cin(1)*100 c6pc1in=cpct(6)/cin(1) c6pc1cvt=cpct(6)/(cin(1)-cpct(1)) !HC converted (%) !mol H2 per mol HC inlet !mol H2 per mol HC converted c5pc1in=cpct(5)/cin !mol CO per mol HC c5pc1cvt=cpct(5)/(cin(1)-cpct(1)) !mol CO per mol HC inlet converted ig=1,igmax cpctd(ig)=cpct(ig)/(1.-cpct(4)) end write(1,56)time,T(irmax,2),T(irmax,izmax), & & (cpctd(ig)*100.,ig=1,igmax-1),(cpct(ig)*100.,ig=1,igmax-1), c1cvt,c6pc1in,c6pc1cvt,c5pc1in,c5pc1cvt write(2,57)time,(T(ir,izmax),ir=0,irmax,jr) ir=0,irmax ctr(ir)=c(1,ir,izmax)+c(2,ir,izmax)+c(3,ir,izmax)+ & c(5,ir,izmax)+c(6,ir,izmax)+c(7,ir,izmax) cir6(ir)=c(6,ir,izmax)/ctr(ir) end write(3,58)time,(cir6(ir)*100.,ir=0,irmax,jr) ir=0,irmax cir1(ir)=c(1,ir,izmax)/ctr(ir) end write(31,58)time,(cir1(ir)*100.,ir=0,irmax,jr) write(*,56)time,T(irmax,izmax),Tg(izmax), & (cpct(ig)*100.,ig=1,igmax-1),cpctd(6)*100 end if C**** Tinh Tt(irmax,iz) for current time step iz=0,izmax Tt(0,iz)=1./(1.+htt(iz)*delr/K)* & (Tt(1,iz)+htt(iz)*delr/K*Theat) end do ir=0,irmax iz=0,izmax T(ir,iz)=Tt(ir,iz) end T(ir,-1)=T(ir,0) T(ir,izmax+1)=T(ir,izmax) end 100 continue !end loop !************************************************** 50 format(1x,A6,2(1x,A6),13(1x,A6),5(1x,A10)) 51 format(1x,A7,7(1x,A8)) 52 format(1x,A7,7(1x,f8.6)) 53 format(1x,A7,3(1x,f15.5)) 54 format(1x,f7.0,2(1x,f8.0),6(1x,f8.2)) 55 format(1x,f6.2,2(1x,f6.0),6(1x,f6.2)) 56 format(1x,f6.0,2(1x,f6.0),13(1x,f6.2),5(1x,f10.3),1x,f5.2) 57 format(1x,f7.2,7(1x,f8.2)) 58 format(1x,f7.3,7(1x,f8.2)) 59 format(1x,f7.0,7(1x,f8.2)) 60 format(1x,I3,1x,f10.0,1x,f15.6,6(1x,f8.2)) 61 format(1x,A5,2(1x,A5),14(1x,A7),1x,A10) 62 format(1x,f5.3,2(1x,f5.0),14(1x,f7.3),1x,f10.8) 63 format(1x,I2,1x,f6.0,6(1x,f6.3),2(1x,f6.0),2(1x,f6.3)) 71 format(92(1x,A6)) 72 format(1x,f6.3,7(1x,f6.0),7(1x,f6.2),7(1x,f6.3), & end 35(1x,f6.2),7(1x,f6.3),28(1x,f6.2)) t ... cỏc loi hydrocarbon thm thng rt c, vớ d nh hydrocarbon cú nhõn benzene cú th gõy ung th n gin a cỏc tiờu chun v mụi trng, ngi ta ch a thnh phn hydrocarbon tng cng khớ thi ( vit tt l TH ) Hydrocarbon... cao 32 2.1.1.6 Mt thp 33 2.2 Phng phỏp sn xut hydro v nhiờn liu giu hydro 38 2.2.1 Phn ng nhit húa nhiờn liu hydrocacbon vi hi nc (phn ng thu nhit ) 39 2.2.2 Phn... th gõy nguy him tớnh mng cho ngi - HC: (hydrocacbon, ụi cũn c ký hiu l CmHn) l cỏc loi hydrocarbon cú nhiờn liu hoc du bụi trn khụng chỏy ht cha khớ thi Hydrocarbon cú rt nhiu loi Mi loi cú mc
Ngày đăng: 24/07/2017, 22:39
Xem thêm: Nghiên cứu mô hình hóa quá trình tạo nhiên liệu giàu hydro để giảm phát thải độc hại của động cơ đốt trong , Nghiên cứu mô hình hóa quá trình tạo nhiên liệu giàu hydro để giảm phát thải độc hại của động cơ đốt trong