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

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