- Dỳng nhiởn liệu sạch
- Loại trừ Nitơ ra khỏi nhiởn liệu - Dỳng hỗn hợp nhũ tương
3.2.2.Phương phõp tõc động vỏo quõ trớnh chõy
Bao gồm 4 phương phõp: - Cải tiến hệ thống nạp vỏ thải - Cải tiến quõ trớnh chõy - Sử dụng hồi lưu khợ xả - Phun nước
3.2.3.Phương phõp xử lý thừng qua cõc thiết bị phợa sau động cơ
Bao gồm 3 phương phõp:
- Sử dụng thiết bị hấp thụ, trung hoỏ độc tố khợ xả - Đốt lại
KẾT LUẬN CHUNG
Với cõc đõnh giõ, phón tợch tổng quan nởu trởn ta thấy rằng mức độ ừ nhiễm mừi trường do giao thừng vận tải trong đụ cụ giao thừng vận tải thuỷ góy ra lỏ rất lớn vỏ ừ nhiễm do độc tố khợ xả động cơ diesel tỏu thuỷ chiếm một phần đõng kể. Tõc hại do chỷng góy ra đối với mừi trường vỏ con người lỏ khừng thể lường hết. Trởn cơ sở đụ việc nghiởn cứu để tớm ra cõc giải phõp tối ưu khừng chỉ lỏ chuyện của một quốc gia, một tổ chức hay một họng động cơ nỏo mỏ lỏ vấn đề chung của cả thế giới. Xu hướng giải quyết vấn đề nỏy rất đa dạng, tuy nhiởn do giới hạn của đề tỏi nởn em chỉ đề cập một cõch tổng quan đến cõc biện phõp giảm thiểu nồng độ độc tố trong khợ xả động cơ diesel tỏu thuỷ. Qua đụ em rỷt ra một số kết luận sau:
1. Việc hạn chế nồng độ độc tố trong quõ trớnh chõy động cơ diesel cần phải được xem xờt đồng thời với việc đảm bảo tợnh kinh tế nhiởn liệu. Hạn chế sử dụng động cơ ở chế độ toỏn tải vỏ tăng tải đột ngột vớ khi đụ vừa tăng nồng độ NOx vỏ bồ hụng, vừa giảm tợnh kinh tế nhiởn liệu.
2. Trong quõ trớnh khai thõc động cơ, cần phải đảm bảo duy trớ trạng thõi kỹ thuật của cõc chi tiết trong hệ thống nhiởn liệu đặc biệt lỏ chất lượng nhiởn liệu vỏ vúi phun để đảm bảo động cơ lỏm việc tốt khừng sản sinh ra nhiều khụi vỏ muội. Việc đảm bảo tốt chất lượng hệ thống khợ nạp giỷp giảm nồng độ NOx trong khợ xả do giảm được nhiệt độ quõ trớnh chõy vỏ nhiệt độ khợ xả.
3. Động cơ cũ thường bị giảm tỷ số nờn do hiện tượng ăn mún, mỏi mún chõy đỉnh piston lỏm tăng thể tợch buồng đốt, lỏm kờo dỏi thời gian chõy trễ dẫn đến tăng tốc độ hớnh thỏnh NOx ở giai đoạn đầu của quõ trớnh chõy. Cần cụ biện phõp điều chỉnh lại thể tợch buồng đốt cho phỳ hợp.
Em rất mong được sự chỉ dẫn đụng gụp ý kiến của cõc thầy, cừ giõo, cõc cõn bộ kỹ thuật, bạn bộ để đề tỏi của em được hoỏn thiện hơn vỏ cụ thể phõt triển đề tỏi lởn một mức cao hơn nữa .
Em xin chón thỏnh cảm ơn !
Người thực hiện đề tỏi Đỗ Thõi Sơn
Phụ lục
' CHUONG TRINH - thay doi tai
'L2t#0 co the tinh cho cac truong hop khi thay doi dieu kien ngoai:dong tai
'(fid=0,hp=0,ntk=22000,ga=.21,va thay cac phu tai)-tang toc(thay cac gia tri
'cua cac thong so dong tai o thoi diem 10 giay vao thong so ban dau ntk,ga,
'fid,hp,nd,fizad)
' - Thay doi dieu kien moi truong voi cac pphuong an khac nhau :
' + Khi giu nguyen ap suat,thay doi nhiet do moi truong P0=103000pa,T0 = 0
' 20,43 do c
' +Khi giu nguyen nhiet do,thay doi ap suat moi truong T0 = 0 hoac 20,40
'va thay doi P0 = 95000,101300,102000.
' - Phu tai thay doi theo chu ky : Tn = .5 + . 5*sin(a*tim)
' a=1.58 - tan so 0,25 Hz ; 3.18 - tan so 0,5 Hz ; 6.28 - tan so 1 Hz
' - Phu tai khong thay doi : Tn = A ; A= .5 ; .75 ; .9 ; 1
' Tr = .5 ; 1 ; 3 ;
' Momen quan tinh tuabin-may nen : J = .027 ; .04 ' Tang ap cao (nn;pitn;picn;pefn;gan;bn):
' - tang ap xung (cac he so k1,k2,k3,k4 >0) ' - tang ap dang ap (cac he so k1=k2=k3=k4 =0)
' Thay doi luong nhien lieu cap vao de thay doi he so alf va cac thong so
' khac hpm = .73 + .07 * tim va hpm = .8 , hpm=1.1 DECLARE SUB koeffimp (pd!, P0!, kn!)
DECLARE SUB ghidia (tim!, n!, alf!, ei!) DECLARE SUB dothi (x1!, y1!, x2!, y2!, col!) 'khung (x1, y1, x2, y2)
DECLARE SUB form (t$, x!, y!, a!, B!, c!, d!, e!, f!, g!, H!, v!, n!, j!, z!, o!, p!, i!, u!, m!, s!)
DECLARE SUB khung (x3!, y3!, x4!, y4!, col!)
'DECLARE SUB turbina (gg!, ggn!, nras!, rr!, ff!, mfen!, p2n!, p0!, n!, ucopt!, ucn!, etn!, flag!, pd!, tg!, nt!) 'DECLARE SUB koeffimp (pd!, p0!, kn!)
DECLARE SUB timgraf (lm!, rm!, um!, dm!, col!, tbeg!, tend!, parbeg!, parend!, t!, par!)
DECLARE SUB comp (T0, P0, pic!(), ec!(), n!, gac!, by!, ay!, bx!, ax!, pic!, ec!)
'DECLARE SUB Mexn (a1!, b!, s!, ndn!, aa!, pefn!, cm!) '********************************************************* *******************
' ’“ŽŠŽŒ…‘‘Ž
'********************************************************* *******************
'ăẼõó óỡ
CONST k = 1.41, pi = 3.1415, r = 287
CONST lo = 14.3, tatm = 273, patm = 101300 CONST kg = 1.35, Rg = 285, khr = 1123
'ư ô Ẩơ đơộẨổ ọõẬẼđẨẪ T0 = 293
P0 = 103000
'ồọăỗẨẨ Ẩ ồẼpỈ óỡ,ồọăỗ.ăẬ đẨộẨ
DEF fni (a, B, c, d) = a + (B - a) * (c - d)
DEF fnras (rr, pp, ff) = (1 - rr / 2 * (pp - 1)) ^ 2 + ff ^ 2
DEF fnpsipr (aa) = SQR(2 / ((k - 1) / k) * (aa ^ (2 / k) - aa ^ ((k + 1) / k)))
DEF fnri (a) = 1.065 - .426 / a ^ (3 / 2)
DEF fntm (nd, pint) = (a1 + b1 * nd) / aa * (pint / pintn) ^ .1 * (km - 1)
'--- ---
fmt1$ = "tim= ##.## nóă= ##### Nó=##### Nă=##### wc=#.#### alfc=#.#### ec=#.####" fmt2$ = " int ###### P1=g ###### tr= #.# jtk=#.### hpmax=#.## fimax=#.######" 'fmt2$ = " €Ậớồ ##.### ’g = ####.#, ăĨôă= #.### ă= #.###" fmt4$ = "nd ###.## ttq= #.### bco= #.#### bsum= #.#### hp #.### ro=#.### be=#.###" '--- --- ' ˆ‘•Ž„›… „€›… ' ôẨươẬớ CONST d = .18, s = .22, e = 11, z = 6 CONST vsz = pi * d * d / 4 * s * (1 + 1 / e) * z CONST tdies = 4'3.4 'alfn = 1.72' 1.96'1.72 alfm = 2' 2.03'2 eim = .4'.43 dei1 = .0536'.214' .0536 dei2 = .133'.107'.133 df = .214'.252'.214 falfn = .145'.0846'.145 falf = .564 'alf = 2.82 'ei = .321 ac = .533 dc = .048: dc1 = .023: dc2 = .012 gc1 = .0024: gc2 = .00413 aw = .785
bw = .0215: bw1 = .011: bw2 = .006 gw1 = .00115: gw2 = .0016 m1 = .42 m2 = .058 ee = 0 ' Ỉơổ. ĨẼóơpẨ a1 = 120000' 60000! bp = 16000'15000! b1 = bp * s / 30 ndn = 750: pefn = 753000' 1200000'810000! 'CALL Mexn(a1, bp, s, ndn, aa, pefn, km) aa = a1 + bp * s * ndn / 30
etmn = 1 / (1 + aa / pefn)'hieu suat co gioi etm km = 1 / etmn PRINT aa, cm ' õ põ CONST tr = .5 CONST kr = 50 t = T0
' ổẼẬẼôẨẬớẨă ôôọđẼốẼíẼ đẼưôọổ tw = 293 etcool = .67: etgxn = .97 ' óọpđẨ nn = 32000 '45000: p2n = P0: pitn = 1.364 '1.85'1.3 p1n = pitn * p2n pitx = 1.07: p2x = P0 p1x = pitx * p2x ' ẼỈẨ Ậớỡơ tgn = 910'700: ein = .867 * .42: ecn = .72
gan = .337' .5'.3375'.44 ' thay lai Pint=1,46;Pitn=1.3;gan=0,337 bn = .0175 '.0142'.0116' bn=Ge/3600=50/3600=0,0139 bcn = bn / ndn * 120 / z: fcn = .95 ggn = gan + bn ' ăẼỈĨpơõõẼp picn = 1.454' 1.95' 1.46 pdn = picn * P0 tdn = T0 * (1 + (picn ^ ((k - 1) / k) - 1) / ecn) dpcooln = (1 - etgxn) * pdn
ksicool = dpcooln * pdn / gan ^ 2 / tdn tintn = tdn - (tdn - tw) * etcool
pintn = pdn - dpcooln
rointn = pintn / r / tintn sprn = SQR(pintn * rointn) pspn = fnpsipr(p1n / pdn)
fin = gan / (fcn * (rointn * ndn * vsz / 30 / 4)) alfn = gan / lo / bn / fin
dgpn = (fin - 1) * gan ' 'ổẼẬ.ổẼô tgx = 415: ndx = 750: ecx = .6'tgx=453 gax = .21: bco = .108: bx = .0026: fcx = fcn 'bx=.00467'.0025;gax=.266.25 'ttq = .65 ' ăẼỈĨpơõõẼp pdx = 106300: picx = pdx / P0 tdx = T0 * (1 + (picx ^ ((k - 1) / k) - 1) / ecx) tintx = tdx - (tdx - tw) * etcool
dpcoolx = (gax / gan) ^ 2 * ((tdx / tdn) / (pdx / pdn)) * dpcooln
pintx = pdx - dpcoolx
fipx = gax / (fcx * (pintx / tintx / r * ndx * vsz / 30 / 4))
alfx = gax / lo / bx / fipx
eixo = fnri(alfx) / fnri(alfn) 'ẨôẨă. ăĨô ĨẼ Ẩổóơpọ dgzx = (1 - fcn) * gax
tmx = fntm(ndx, pintx) bcxo = tmx / km / eixo
' ăẼợồồẨỗẨơóỡ ôẬủ ọốơó ẨỈĨọẬớõẼõóẨ CONST mega = 1000000 k1i = 0' .005 k2i = 0' .0015 k3i = 0' .033 k4i = 0' .012 rr = .4' õóơĨơớ pơ ẳẨđẼõóẨ ff = .82
nras = fnras(rr, pitn, ff)
w = (1 / pitn) ^ (2 / kg) - (1 / pitn) ^ ((kg + 1) / kg) psi = SQR(2 * (kg / (kg - 1)) * w) mfen = ggn / SQR(p1n ^ 2 / Rg / tgn) / psi CONST jtk = .04, rk = .07 CONST ku = pi * rk / 30 un = ku * nn hrn = khr * tgn * (1 - (p2n / p1n) ^ ((kg - 1) / kg))' p đẼóa cn = SQR(2 * hrn) ucopt = .42 + .28 / .6 * rr
ucn = un / cn / ucopt etn = .72
' ‡€ƒ“‡Š€ •€€Š’…ˆ‘’ˆŠ ŠŽŒ…‘‘Ž€ CLS
namf$ = "d:fi2.dat"
'INPUT "‚đơôẨóơ ưđ Ẩơ ồ ẪẬ ô ỡổ ", namf$: namf$ = namf$ + ".dat"
PRINT "‡ ípọẶ ơóõủ ồ ẪẬ "; namf$; " ... † „ ˆ ’ … ..." OPEN namf$ FOR INPUT AS #1
INPUT #1, ng, nv, ay, by, ax, bx DIM ec(ng, nv), pic(ng, nv)
FOR y = 0 TO nv ' õẨưọ đđơpổ n = ay * y + by: PRINT n ' õẬơđ Ĩp đẼ FOR x = 0 TO ng ga = ax * x + bx
INPUT #1, pic(x, y), ec(x, y) NEXT
NEXT CLOSE #1
CLS : SCREEN 2
LOCATE 25, 1: PRINT "ẨõẼđ óớ ĨẼẬơ ổ pakóơpẨõóẨă ( Y/N ) ?" GOSUB os:
IF q$ = "Y" OR q$ = "y" OR q$ = "" OR q$ = "" THEN GOSUB figraf:
GOSUB osq: CLS
p$ = " Cốơó Ĩpơpđ . ‡ ăẼốẨóớ õốơó ( Y/N ) ? "
DIM tim(1200), ei(1200), pint(1200), tg(1200), nd(1200), alf(1200), npr(1200), B(1200), be(1200), tint(1200), ttq(1200), Ge(1200), etm(1200), ete(1200), Ne(1200) vcom:
n = 22000
tint = tintx '320'tintx ga = .21' gad = ga pitr = pitx pit = pitr pic = picx pc = P0 * pic p1 = p1x fid = 0'-.212'-.2 '-.019' 0 hp = 0'.581' 0 ' .476 ws = (pi * n / 30) ^ 2 nd = 750' 591'600'736' 750
wd = nd / ndn H = .01
tim(0) = 0: bsum = 0: fimax = 0: alfmin = 0: ndmin = 0: fimin = 0
hr = .00001
grafinterval = .1: timg = 0 'ttq = .65
'alf,tg=?
alfint = 0'ẼđọẬủơóõủ Ẩóơíp ẬớỡẪ ĨẼă ư óơẬớ ôỡỈẼõóẨ alfdop = 1.5'ôẼĨọõóẨỈẼơ ư ốơẨơ Ậớồ ĨẼ ôỡỈẼõóẨ wyc = 78
wn = 0: eis = 0: kozs = 0: boss = 0 i = 0 'beep '--- --- DO' €—€‹Ž –ˆŠ‹€ i = i + 1 tim = tim + H IF tim >= .01 THEN H = .01 Esc$ = INKEY$
IF Esc$ = CHR$(27) THEN LOCATE 25, 1: PRINT p$; : GOSUB os: IF q$ = "Y" OR q$ = "y" THEN EXIT DO
DO' Žˆ‘Š õẼđỈ.p đẼóỡ
IF ga >= .2 AND n >= 20300 THEN
CALL comp(T0, P0, pic(), ec(), n, ga, by, ay, bx, ax, pic, ec) td = T0 * (1 + (pic ^ ((k - 1) / k) - 1) / ec) pd = pic * P0 + 10 ELSE flag = 1: td = T0: pd = P0 END IF ' „ˆ‡…‹œ
' ổẼẬẼôẨẬớẨă ôôọđẼốẼíẼ đẼưôọổ 'td,pd - ĨẼõẬơ ’Š
'tint,pint - ĨẼõẬơ ổẼẬẼô-ă tint = td - (td - tw) * etcool
dpcool = ksicool * gad ^ 2 * td / pd pint = pd - dpcools
' p õổẼô ủ ổ p ẳơpẨõóẨă ôẨươẬủ nd = 750
roint = pint / r / tint
gad1 = roint * nd * vsz / 30 / 4 spro = SQR(pint * roint) / sprn IF pspn = 0 THEN pspn = 1
IF (pint - p1) > 0 THEN
dg = spro * fnpsipr(p1 / pint) / pspn * wd * dgpn ELSEIF (pint - p1) < 0 THEN
dg = -spro * fnpsipr(pint / p1) / pspn * wd * dgzx END IF
gad = (gad1 + dg)
gadpr = gad * patm / P0 * SQR(T0 / tatm)
IF dg > 0 THEN fipr = gad / gad1 ELSE fipr = 1 delg = gad - ga
'IF ABS(gad - ga) > (.05 * ga) THEN ga = gad
'IF 7.32 >= tim >= 7.2 THEN LPRINT USING " pic =#.### Ga=#.### dg=#.###"; pic; ga; dg;
'LOCATE 23, 2: PRINT USING " Gd =#.### Ga=#.### dg=#.###"; gad; ga; dg;
IF ABS(gad - ga) < (.05 * ga) THEN EXIT DO ga = (ga + gad) / 2
LOOP
'gad = ga' ăẼơỗ ĨpẨđẬẨẶơẨẪ '---
' ô.ọ. ôẨươẬủ Ẩ pơíọẬủóẼp
'IF tim <= 1 THEN tn = tim ^ .5 ELSE tn = tim ^ 1.3 / (tim ^ 1.29 + tim ^ .4 - 1)
'tn = .5 + .4 * SIN(1.58 * tim) 'tn = (.5 + .4 * SIN(.3 * tim)) 'tn = .5 + .4 * SIN(6.28 * tim) 'IF tim <= 4 THEN tn = .4
'IF tim > 4 THEN tn = .6 'IF tim > 9 THEN tn = .85
'IF tim > 29 THEN tn = .4 + .1 * SIN(1.58 * tim) 'IF tim <= 6 THEN tn = .4
'IF tim > 6 THEN tn = .85 'IF tim <= 6 THEN tn = .5 'IF tim > 6 THEN tn = .85 'IF tim <= 6 THEN tn = .6 'IF tim > 6 THEN tn = .85 tn = .5
dfid = (ti - tn - tm) / tdies
dhp = (-33 * (fid - fizad) - hp) / tr '50,33 fid = fid + dfid * H
wd = 1 + fid
IF wd > 1 THEN wd = 1 nd = wd * ndn
IF wnd <= wyc THEN delwnd = (wyc - wnd) * H / wyc wn = wn + delwnd
wc = wn / tim'ĨẼă ư óơẬớ ă ốơõóđ
IF ABS(fid) > ABS(fimax) THEN fimax = fid ndmin = (ndn + (fimax * ndn))
IF ABS(fid) < ABS(fimin) THEN tpr = tim hp = (hp + (dhp * H))
IF hp > 1.1 THEN hp = 1.1 hpmax = 1.1
IF hp < 0 THEN hp = 0
IF hp < -bcxo THEN hp = -bcxo
bco = (bcxo + hp) / (bcxo + 1)' ỗẨăẬẼđ ủ ĨẼô ố bc = bco * bcn B = bc * nd * z / 120 Ge = B * 3600 bsum = bsum + B * H dee = .1047 * nd * (ti - tm) * H Ne = .1047 * nd * (ti - tm) ee = ee + dee IF ee = 0 THEN ee = .00001 bm = bsum / ee IF B <= .0001 THEN B = .0001 alf = gad / lo / B ' Ậớồ
IF alf > alfmin THEN alfmin = alf IF alf > 20 THEN alf = 20
IF bco > 1 THEN bco = 1 IF alf <= 1.4 THEN dc = 0
' ĨpẨ Ĩ ôơẨẨ Ậớồ ẨẶơ alfdop đăẬùố ơóõủ õốơóốẨă Ẩ- ó.ĨẼă ư óơẬủ ôỡỈẼõóẨ
IF alf <= alfdop THEN delalf = (alfdop - alf) * H 'ELSE delalf = 0
delalf1 = (alfdop - alf) * H alfint = alfint + delalf alfd = alfint
alfdc = alfint / tim 'ĨẼă ư óơẬớ ôỡỈẼõóẨ ' ợăẼẼỈẨốơõă ủ ổ p ẳơpẨõóẨă ôẨươẬủ
' eio = fnri(alf) / fnri(alfn) 'ẨôẨă. ăĨô ĨẼ Ẩổóơpọ,ph. trinh cu
'IF alf < 3.8 THEN eio = fnri(alf) / fnri(alfn) ELSE eio = fnri(alfn) / fnri(alf)
'eio = (1.065 - (.426 / alfn ^ 3.2)) / (1.065 - (.426 / alf ^ 3.2))
'ei = eio * ein
ei = .267 * (alf - .1) ^ (1 / (alf + .2)): eio = ei / ein eis = eis + ei * H
eic = eis / tim
' óơỈĨơp óọp ủ ổ p ẳơpẨõóẨă ôẨươẬủ
IF alf < 1.2 THEN enst = 1.2 - alf ELSE enst = 0 falf = falfn + df * (alf - alfn)
IF falf > .35 THEN falf = .35
tg = (.8 * tint) + 2660 * (1 - enst - ei - .1) / (1.05 * alf)
' Ỉơổ.ĨẼóơpẨ tm = fntm(nd, pint)
' tm = (a1 + b1 * nd) / aa * (pint / pintn) ^ .1 * (km - 1) ti = km * bco * eio etm = 1 - tm / ti ete = ei * etm ttq = (ti - tm) * 2012 Ne = (ttq * nd / 9550) be = Ge / Ne SCREEN 8 IF tim = .01 THEN
CALL form("ˆõổ.ô ỡơ tim=.01", 3, 62, nd, (ti - tm), fid, hp, bco, alf, tg, pint, ei, bsum, nn, pitn, picn, pefn, gan, bn, T0, P0) END IF IF flag <> 1 THEN lc = 1005 * (td - T0)' p đẼó ăẼỈĨpơõõẼp nc = gad * lc' ỈẼờẼõóớ ăẼỈĨpơõõẼp ELSE nc = 5E-10 * n ^ 3 * roint END IF vtur: ' ’“ˆ€ gg = gad + B' / 3600
'SUB turbina (gg, ggn, nras, rr, ff, mfen, p2n, p0, n, ucopt, ucn, etn, flag, pd, tg, nt)
DO
' p õổẼô ủ ổ p ẳơpẨõóẨă Ẩ p2
mfe = SQR(nras / fnras(rr, pit, ff)) * mfen p2 = (gg / ggn) ^ 2 * (p2n - P0) + P0
' ẼĨpơô. ó
a = (gg / p2 / mfe) ^ 2 * Rg * tg / 2 / (kg / (kg - 1)) x = .5 + SQR(a + .25)
pitr = x ^ (kg / (kg - 1)) dpit = ABS(pit - pitr)
IF dpit < .005 THEN EXIT DO pit = (pit + pitr) / 2: LOOP
pit = pitr p1 = pit * p2
' p õĨ.p đẼó óọpđẨỡ
hr = khr * tg * (1 - (1 / pit) ^ ((kg - 1) / kg)) ' ă.Ĩ.ô. óọpđẨỡ
u = ku * n: cad = SQR(2 * hr) uc = (u / cad) / ucopt
et = (uc / ucn) * ((2 - uc) / (2 - ucn)) * etn IF et <= 0 THEN et = 0 IF flag = 1 THEN kn = 1.3 ELSE pdiz = (pd - P0) / mega kh = 1 + k1i / pdiz kgi = 1 + k2i / pdiz
etpop = 1 + k3i * (1 + k4i / pdiz)