M ụcălục
4.3. Mơ hình sự cháy tia phun Diesel
Sựăcháyăđ ợcăđiềuăkhiểnăb iăcácăyếuătốăsau, theo thứătự:
• Sựăbayăhơiăcủaăgiọt nhiên liệu. • Hỗnăhợpăkhơng khí - nhiên liệu • Sựăcháy
• Cung c p khơng khí • Dậpătắt
81
4.3.1. Sự bay h i c a giọt nhiên li u.
Sự bayăhơiăcủa các giọt nhiên liệu đặt một h n chế b sung về tỷ lệ hỗn hợp
khơng khí - nhiên liệu. Mốiă t ơngă quan thực nghiệm đ ợc sử dụngă để xácă định đ ng kính trung bình của một kỳ cho các giọt nhiên liệu trong một gĩi nhiên liệu - khơng khí, cùng với tốcăđộ bayăhơi. Khơng cĩ dữ liệu là cần thiết cho một phần của tính tốn. Các tính ch t nhiên liệu đ ợc bao g m trong ch ơngătrình.
4.3.2. H n h p khơng khí - nhiên li u.
Hỗnăhợpăkhơng khí - nhiên liệuăđ ợcămơăt ăbằngăcácăcơăchếăhình thành phun.
Sau mộtăth iăgianănghỉăngắnăngủi nhiên liệuăt oăthành mộtăphun,ăsựăxâmănhậpăcũng đ ợcăxácăđịnhăb iămộtăsựăt ơngăquanăthựcănghiệm.ăTr ớcăkhiăthâmănhậpăphá vỡălà tuyếnătínhătheoăth iăgian,ăsauăkhiăphá vỡănĩ thayăđ iăvớiăcĕnăbậcăhaiăcủaăth iăgian.ă Cơng thứcătiêu chuẩnăchoăcácăcơăchếăphun,ădựaătrên độngăl ợngădịng ch y liên tục,ă
đ ợcăsửădụng. Tuy nhiên, sựăbiếnăđ iăđểăđ aăvàoămộtăl ợngăkhơng khí chuyểnăđộngă
trong bu ngăcháy,ănh ăxốy.ăHaiăyếuătốăng iăsửădụngăquyăđịnhăkhơng khí hút vào
đ ợcăsửădụngănh ănhânătốăcuốnătheoăkhíătrongăch ơngătrình. Yếuătố "sựăcháyătr ớcă
khi hút khơng khí vào"ăđ ợcăsửădụngătrongăth iăgiană khơngă khíăhoặcăbayăhơiăsựă cháy cĩ kiểmăsốtătrongăkhiăcác yếuătố "sựăcháyătrongăquáătrình hút khơng khí vào"
đ ợcăsửădụngătrongăth iăgianăkhơngăcháy, ngọnălửaăth iăraăngoài, hoặcăch mădứtă
quá trình đốtăcháy.ăGiáătrịămặcăđịnhăchoăhaiăthơngăsốănày là 1,8. Mộtăgiáătrịăcủaă1ăcĩă
thểăđ ợcăsửădụngăđểămơăt ămộtăphunătrongăkhơngăkhíătĩnhălặng.
Việc sử dụng hai yếu tố này đ ợc tiếp tục sửaăđ i b i yếu tố "sự cháy khơng khí hút vào", trongăđĩ điều chỉnh tỷ lệ hỗn hợp khơng khí s ch để khí cháy bị cuốn theo vào phun. σĩăt ơngăđ ơngăvới tỷ lệ của khí cháy để khơng khí s ch ch a bị cuốn theo trong quá trình đốt cháy khơng gian bên ngồi phun. Nếu phun đập vào t ng nĩ đ ợc gi định di chuyển dọc theo nĩ, sự hút khơng khí vào và xâm nhập t i tỷ lệ của sự phun tĩnh, nhân với yếu tố "phun thâm nhập vàoăt ng". σg i dùng nhập vào hoặc dự đốn "Tỷ lệ xốy lốc" đ ợc sử dụngăđể làm gi m sự phun th m và tốc độ.
82
4.3.3. Sự cháy.
Sự chậm trễ cháy, , đ ợc l y từ một mốiăt ơngăquan thực nghiệm. Sự cháy x y ra khi:
∫ = 1 ( 4.11) với:ăăăăăăăτ = AP ∅ e /
Trongăđĩ: P áp su t xy lanh.
Gĩi tỷ lệt ơngăđ ơng.
T Nhiệtăđộ xy lanh.
A, B, C, D hằng số thực nghiệm.
Sự chậm trễ cháy tối đa cho phép là:
= 1
( 540∗ ) [ ] ( 4.12)
Trongăđĩ:ă FIGN ng i dùng nhập vào yếu tố chậm trễ cháy.
4.3.4. Sự hình thành NOx.
Tr ng thái của khí trong các gĩi đ ợc mơ t b i cácăđiều kiện cân bằng của các d ng tiêu chuẩn của WAVE áp su t và nhiệtăđộ hiện t i. Thơng tin này là cần thiết cho việc tính tốn nhiệtăđộng lực của phát hành nhiệt, phân tử phân ly là một yếu tố quan trọng.
Nĩ cũng đ ợc sử dụng nh ămột phần của cơă chế Zeldovich m rộng cho sự hình thành NOx. Các thành phần của khơng khí s ch đ a vào độngăcơ là cần thiết nh ălàădữ liệu, về phân số mol của oxy, nitrơ, carbon dioxide vàăn ớc. Sựthayăđ i này của WAVE về khí th i tái l uăthơng (nếu dữ liệu WAVE chứa thơng tin về một hệ thống tái tuần hồn) và khí cịn l i bị mắc kẹt trong xy lanh khi van hútăđĩng. Hình thành oxit nitrơđ ợc tính cho mỗi gĩi nhiên liệu - khơng khí dùng riêng nhiệt độ cục bộ. Cơăchế Zeldovich m rộng đ ợc sử dụng:
N + NO → N + O ( 4.13) N + O → NO + O ( 4.14) N + OH → NO + H ( 4.15)
Việc tính tốn số cĩ thể đ ợc tìm th y là khơng năđịnh. Một triệu chứng điển hình của số b t n là một giá trị oxit nitrơ r t th p (ví dụ, ít hơnă10 ppm). năđịnh
83
cĩ thểđ ợc phục h i bằng cách thiết lập ng i dùng nhập vào "yếu tố cháy trễ" lớn hơnăgiáătrị mặc định của nĩ.
4.3.5. Sự hình thành b hĩng.
Tốcăđộ hình thành l ới b hĩng bao g m hai phần: hình thành và tốcăđộ quá trình oxy hĩa:
= − ( 4.16)
Trongăđĩ:ă Ms hình thành b ăhĩngăl ới.
Msf sựăhình thành b ăhĩng.
Mso quá trình oxy hĩa b ăhĩng.
Sự hình thành b hĩng đ ợc dựa trên cơng thức Hiroyasubanăđầu, tức là: = S ∗k M vớiăk = A e ( 4.17)
Trongăđĩ:ă Sf yếuătốăquyămơăhình thành b ăhĩngăng iădùng nhậpăvào.
Mfv khốiăl ợngănhiên liệuăbốcăhơi.
Af hằngăsốămũătr ớcăArrheniusă(150ăin/sec).
Tf nhiệtăđộăkíchăho tăchoănhiên liệuăDieselă6290ă0K.
Mơ hình quá trình oxy hĩa b hĩng đ ợc dựa trên cơăchế Nagle và Strickland- Constable. Cho quá trình oxy hĩa carbon cĩ hai lo i vị trí trên bề mặt carbon: lo i ph n ứng A hơn một vị trí và lo i ít ph n ứng B. Với sựgiaătĕng của nhiệtăđộ, vị trí lo i A đ ợc sắp xếp l i để thành vị trí lo i B. Tốcă độ ph n ứng bề mặt (in g carbon/s/cm2) là:
=
1 + ∗ + ( 1− ) ( 4.18)đây,ăphầnăcủaăbềămặtăx yăraăb iăvịătríălo iăA,ăxăđ ợcăchoăb i: đây,ăphầnăcủaăbềămặtăx yăraăb iăvịătríălo iăA,ăxăđ ợcăchoăb i:
= 1 + ( 4.19)
Trongăđĩ:ă Po2 phầnăápăsu tăcủaăoxyă[atm].
kA hằngăsốătốcăđộă=ă20exp(-15100/T) [g/cm2/s/atm].
84
kT hằngăsốătốcăđộă=ă1,51*105exp(-48830/T) [g/cm2/s].
kZ hằngăsốătốcăđộă=ă21,3exp2064/T) [1/atm].
Tốcăđộăoxyăhĩaăb ăhĩngăđ ợcătínhăb i:
= ∗6 ( 4.20)
Trongăđĩ:ă Ms khốiăl ợngăl ớiăb ăhĩng.
MWC trọngăl ợngăphânătửăcủaăcarbon.
khốiăl ợngăriêng củaăb ăhĩng.
Dsđ ngăkínhătrungăbình củaăb ăhĩng.
Sc yếuătốăquyămơătốcăđộăoxyăhĩaăb ăhĩngăng iădùng nhậpăvào.
Nh n xét chung.
Số l ợng chỉ định báo cáo trong các tập tin đầu ra và tĩm tắt (ví dụ IMEP, ISFC, Chỉ định cơng su t, NOx) sử dụng cơng việc đ ợc chỉ định từ cácăđộngăcơă diesel phun mơ hình xy lanh, cĩ thể là khác nhau hơnăsoăvới WAVE, kết qu xy lanh báo cáo cho mơ phỏng t ơngătự.
Mơ hình phun cháy Diesel dựa trên ng iădùngăxácăđịnh c u hình khốiăl ợng nhiên liệu l u cho các vịi phun đínhăkèm và đ ng kính vịi phun, đ ợcăđịnh nghĩaă trong phun đầu vào mơ hình độngăcơăDiesel, để cĩăđ ợc cácăđặc tính phun nh áp su t vịi phun nhỏ giọt, vận tốc phun, và phun nhiên liệu kíchăth ớc giọt . Nĩ khơng sử dụng ng i dùng nhập vào c u hình áp lực phun của vịi phun.
Quy mơ C u hình WAVE tự động phun l uăl ợng theo chiều dọc dựa trên l uă l ợng nhiên liệu độngă cơ, nh ng khơng cĩ kh nĕng m rộng quy mơ th i gian phun. Điều quan trọng là c u hình phun (th iăgian)ăđ ợc sử dụng đĩălà phù hợp với tiếp nhiên liệu (t i) khi mơ phỏng. Nếu một th i gian phun nhiên liệu dài đ ợc sử dụng trongăđiều kiện t i nhẹ mà độngăcơ thực tế sẽ cĩ một th i gian phun ngắn, mơ hình phun sẽ tính tốn giá trị sai của đặcăđiểm phun với kết qu m tăđộ chính xác và cĩ thể khơng năđịnh [37].
85
4.4. Mơ hình phân tích sự gi i phĩng nhi t c a quá trình cháy. 4.4.1. Sự Nén.
Từ khi bắtăđầu nén và ng i dùng nhập vào "Bắtăđầu từ gĩc quay trục khuỷu cho đến việc tính tốn gi i thốt nhiệt", áp su t và nhiệtăđộ đ ợc tính tốn sử dụng mơ hình một khu vực nh ăsau:
1. B ớc về phíaătr ớc trong th i gian sử dụng ng i dùng nhập vào " Kích th ớc B ớc cho b ng tra cứu áp su t xy lanh".
2. Xácă định nĕngă l ợng bên trong bằng cách tính tốn cơng việc, truyền nhiệt, và enthalpy trong dịng ch y do phun nhiên liệu và sựbayăhơi nhiên liệu lỏng. 3. Tính tốn áp su t và nhiệtăđộ từ nĕngăl ợng nội bộ trên.
4. Lặp l iăcácăb ớc từ1ăđến 3ăchoăđến khi mơ hình đ tăđến ng i dùng nhập
vào "Bắtăđầu từ gĩc quay trục khuỷu cho đến việc tính tốn sự thốt nhiệt". 5. T i ng i dùng nhập vào " Bắtăđầu từ gĩc quay trục khuỷu cho đến việc tính tốn sự thốt nhiệt", so sánh áp su tăđoăvà tính tốn. Nếu lỗi áp lực lớnăhơn 0,01%, quay tr l iăb ớc 1 vàăđiều chỉnh áp lực bị mắc kẹt (và nhiệtăđộ) khi bắtăđầu nén để phù hợp với áp lực đoăt i " Bắtăđầu từ gĩc quay trục khuỷu cho đến việc tính tốn sự thốt nhiệt ".
6. Lặpăđiălặp l i quá trình này choăđến khi áp lực tính tốn phù hợp với áp
su tăđo t i "Bắtăđầu từ gĩc quay trục khuỷu cho đến việc tính tốn sự thốt nhiệt" trong vịng 0,01 % [37].
4.4.2. Sự Cháy.
Giữa ng i dùng nhập vào "Bắtăđầu từ gĩc quay trục khuỷu cho đến việc tính tốn sự thốt nhiệt" và EVO, áp su tăđo đ ợc sử dụng trong khi nhiệtăđộ và tốcăđộ gi i phĩng nhiệt đ ợc xácăđịnhănh ăsau:
1. B ớc về phíaătr ớc trong th i gian sử dụng ng i dùng nhập vào
"B ớc Kíchăth ớc b ng tra cứu áp su t xy lanh".
2. Tính tốn nhiệtăđộ từ ph ơngătrình tr ng thái.
3. Xácăđịnh nĕngăl ợng nội bộ bằng cách tính tốn cơng việc, truyền nhiệt, và enthalpy trong dịng ch y do phun nhiên liệu và sự bayăhơi nhiên liệu lỏng.
86
4. Tính tốn việc thốt nhiệt dựa trên quy luậtăđầu tiên của nhiệtăđộng lực học.
5. Xác định tỷ lệ nhiên liệu cháy từ tỷ lệ gi i thốt nhiệt dựa trên việc đốt cháy tỷ lệ khơng khí - nhiên liệu tức th i và các thuộc tính nhiệtăđộng lực học.
Đối với phần kín của chu kỳ độngăcơ, sự gi i phĩng nhiệt thuăđ ợc bằng cách
áp dụng định luật đầu tiên của nhiệtă động lực học cho hỗn hợp trong xy lanh, th ng nĩi rằng sựthayăđ i trong nội nĕngă(dE) bằng sự truyền nhiệt vào hệ thống ( Ql) trừđiăcác cơng việc đ ợc thực hiện b i hệ thống ( W).
= − − ( 4.21) Sự thayăđ i này trong nội nĕngăcĩăthể đ ợc chia làm hai phần: phần hình thành, dEf, và một phần hợp lý, dEs. Sự thayăđ i nĕngăl ợng trong phần hình thành là do ph n ứng đốt cháy, trong khi sựthayăđ i nĕngăl ợng trong các phần hợp lý là b i sự thayăđ i về nhiệtăđộ. Định luật đầu tiên sauăđĩăcĩăthểđ ợc viết l i nh ăsau:
+ = − − ( 4.22) Sự thayăđ i nĕngăl ợng trong các phần hợp lý do sự thayăđ i về nhiệtăđộ cũngă đ ợc biếtăđến là:
= ∆ ( 4.23)
Trongăđĩ:ă m khốiăl ợngătrongăxy lanh.
Cv nhiệtădungăriêng t iăthểătíchăkhơngăđ i.
T sựăthayăđ iănhiệtăđộ.
Từph ơngătrình tr ng thái, sựthayăđ i của nhiệtăđộ () chúng tơi xácăđịnh b i sựthayăđ i trong áp lực ( P), thể tích ( V), và hằng số khí ( R):
∆ = ∆
+ ∆
−∆ ( 4.24) T i mỗi b ớc th i gian, cơng việc thực hiện đ ợc tính bằng cách sử dụng các giá trị áp lực cung c p trong các ng i dùng nhập vào d u hiệu áp lực:
∆W = 1
87
Từđịnh nghĩa của hệ số truyền nhiệt, m t nhiệt trong quá trình từng b ớcătĕng gĩc quay để ống lĩt xy lanh, piston, và đầu xy lanh cĩ thể đ ợc tính bằng:
∆Q = Ah T −T ∆t ( 4.26)
Trongăđĩ:ă A diệnătíchăcủaăpiston,ăốngălĩt,ăhoặcăđầuăxyălanh.
h hệăsốătruyềnănhiệt.
Tw nhiệtăđộăbềămặtăt ơngăứng.
Tg nhiệtăđộăkhíăcháy.
t kíchăth ớcăb ớcăth iăgian.
Hệ số truyền nhiệt, h, đ ợc tính bằng cách chọn sử dụng mơ hình truyền nhiệt WAVE. Sauăđĩ, sựthayăđ i nĕngăl ợng của phần hình thành (sự thốt nhiệt rõ ràng) cĩ thểđ ợc tìm th y từ các ph ơngătrình sauăđây:
= − − − ( 4.27) Từ việc gi i phĩng nhiệt rõ ràng chúng ta cĩ thểxácăđịnh khốiăl ợng nhiên liệu đ ợc tiêu thụ. WAVE sử dụng ngay lập tức tỷ lệ A / F theo tính tốn của mơ hình hĩa họcăl ợng pháp và các d ng tính ch t nhiệtăđộng lực họcăđể tính tốn bao nhiêu khốiăl ợng nhiên liệu và khơng khí s ch ph iăđ ợcăđốtăcháyăđể đ tăđ ợc những thay đ i cần thiết trong sự hình thành nĕngăl ợng .
Cho quá trình đốt cháy khơng hịa trộn tr ớc, dữ liệu phun nhiên liệuăđ ợc sử dụngăđể đánhăgiáăviệcăđốt cháy tỷ lệ khơng khí - nhiên liệu tức th i, và chuyểnăđ i tốcăđộ nhiên liệu cháy vào tốcăđộ gi i phĩng nhiệt. Nếuăng i dùng chỉ định dữ liệu phun nhiên liệu khơng phù hợp với tỷ lệ gi i phĩng nhiệtăthuăđ ợc từ những áp su t xy lanhăđo đ ợc, nhiên liệu cĩ thể t m th i ch y ra ngồi trong quá trình đốt cháy, trongăđĩăsẽ (trong Wave 7.1 hay tr lên) ch m dứt việc tính tốn tỷ lệ cháy sớm. Khiăđĩătốcăđộ cháy ch m dứt sớmăđ ợc áp dụngănh ăđịnh luật sự cháy, nĩ lầnăl ợt, cho áp su t xy lanh đầuăvàoăcaoăhơnă áp su t banăđầu nhập vào.ăĐể tránh xu t hiện
x yăraăđiềuăđĩ, WAVE gi phát hiện số l ợng nhiên liệu cĩ sẵn trong phân tích sự
gi i phĩng nhiệt. Nếu nhiên liệu t m th i ch y ra x y ra, lập l i mơ hình tỷ lệ phun bắtăđầu sớmăhơnăvàăth i gian phun ngắnăhơn.ăCácăquá trình trên chỉ áp dụng cho các
88
phân tích thốt nhiệt và lặpăđiălặp l i mỗi chu kỳ độngăcơăchoăđến khi khơng cĩ nhiên liệu t m th i ch y ra ngồi. Cách tiếp cận này giúp c i thiện mơ hình m nh mẽ, bỏ qua các hiệu ứng phụ của sựbayăhơiănhiênăliệu trên sự cân bằng nhiệt.
Cho các mơ hình phun trực tiếp, dữ liệu phun nhiên liệuăđ ợc sử dụngăđể đánhă giá việcăđốt cháy tỷ lệ khơng khí - nhiên liệu tức th iăvàăsauăđĩăchuyểnăđ i tỷ lệ thốt nhiệt thành tỷ lệ nhiên liệuăđốt cháy. Nếu dữ liệu phun nhiên liệuăng i dùng quyăđịnh khơng phù hợp với tỷ lệ gi i phĩng nhiệtăthuăđ ợc từ những áp su t xy lanh ng i dùng chỉ định, cĩ thể cĩ một thiếu hụt t m th i của nhiên liệu trong quá trình đốt cháy, mà sẽ ch m dứt việc tính tốn tỷ lệ cháy sớm. Khi một biên d ng đ ợc t o ra từ một phân tích gi i phĩng nhiệt với một ch m dứt sớmăđ ợc thơng qua vào xy lanhăđộngăcơănh ăbiên d ng cháy, nĩ t o ra một áp su t xy lanh caoăhơnă so với biên d ng áp lực xy lanhăbanăđầu.
Để tránh v năđề này, WAVE theo dõi l ợng nhiên liệu cĩ sẵn của nhiên liệu trong phân tích sự gi i phĩng nhiệt. Nếu nhiên liệu t m th i ch y ra ngồi, mơ hình điều chỉnh biên d ng phun nhiên liệu bằng cách sử dụng thuậtătốnăsauăđâyă(ăvới n đ i diện cho số chu kỳ tĕngălênăkể từ khi thiếu nhiên liệu bắtăđầu ):
1. Thay thế ng i dùng nhập vào biên d ng phun nhiên liệu với một biên d ng hình chữ nhật.
2. Thiết lập sự kh iăđầu của th i gian phun bằng sự bắtăđầu của th i gian cháy.
3. Thiết lập th i gian biên d ng phun bằng 25 CA °(Crank Angle).
4. Điều chỉnh sự bắtăđầu của th i gian phun bằng cách sử dụng ph ơngătrình sauăđây (với một giới h n của SOC-20):
SOI = SOC - n * 2 deg (4.28)
5. Điều chỉnh th i gian biên d ng phun bằng cách sử dụng ph ơngătrình (với một th i gian khơngăítăhơnă10 ° CA):
Th iăl ợng phun = 25 * 0,9 * (n - 1) (4.29)
Các quá trình trên chỉ áp dụng cho các phân tích gi i phĩng nhiệt và lặpăđiălặp l i mỗi chu kỳ độngăcơ choăđến khi khơng cĩ tình tr ng thiếu nhiên liệu trong phân
89
tích gi i phĩng nhiệt. Cách tiếp cận này giúp c i thiện mơ hình m nh mẽ, bỏ qua các hiệu ứng phụ của sự bayăhơi nhiên liệu trên sự cân bằng nhiệt [37].
4.4.3. Xác đ nh đúngđiểm ch t trên (TDC).
Nếu một biên d ng áp su t ơ tơ đ ợc nhập và lựa chọn ơ tơ đã đ ợc thiết lập cho các mơ hình, các mơ hình gi i phĩng nhiệt cĩ thể đ ợc sử dụng để tính tốn điểm chết trên (TDC) bùăđắp, đĩ là sự khác biệt giữa vị trí điểm chết trên (TDC) đúngăsự thật và cĩ hệ thống dữ liệu đĕngănhập chỉ đ ợc điểm chết trên. Điều quan trọng là biếtăđ ợc vị trí TDC đúng nh chính xác càng tốt vì hiệu su t chỉ định của độngăcơ là r t nh y c m với TDC lỗi - một lỗi 1 độ vị trí TDC cĩ thể gây ra một lỗi 10 % trong tính tốn IMEP.
Khi mộtăđộngăcơ đ ợc kh iăđộng, áp su t xy lanh tốiăđa sẽ x yăraătr ớc khi đúngăTDC. Điều này là b i vì khi piston tiến gầnăđỉnh của kỳ đĩăthì sự m t mát nhiệt từ nén (vì nĩng) khí vào các bứcăt ng bu ng cháy sẽ làm gi m áp lực trong đĩ, áp lựcăđĩălà caoăhơnăsoăvới sựgiaătĕng áp lực do sự chuyểnăđộng nhỏ của piston. Sự khác biệt trong gĩc giữa đoăáp lực tốiăđacĩăđộngăcơ và thể tích xy lanh tối thiểu đ ợc gọi là gĩc t n th t nhiệtăđộng lực học. Điểm chết trên bùăđắp do truyền nhiệt cĩ thểđ ợcă ớc tính trực tiếp bằng cách gi i ph ơngătrình nĕngăl ợng (4.30) [37]:
(γ −1)dQ
d + γPdV
d = −VdP
d ( 4.30) Trongăđĩ:ă gĩc quay trụcăkhuỷu, tỷălệănhiệtădungăriêng.
áp su t tốiăđa, dP / d = 0,ădoăđĩph ơngătrình tr thành:
(γ −1)dQ
d = −γPdV
d ( 4.31) Các điều kho n của truyền nhiệt và sựthayăđ i trong thể tích cĩ thểđ ợcăđánhă giá nh ăsau:
= ∑
= ∑ ℎ −
( 4.32) Và
90
=
= + 1− + 1 −
≈ 1 + ∆ ( 4.33)
Trongăđĩ:ă vậnătốcăgĩcăquayătrụcăkhuỷu.
Acyl diệnătíchălịng xy lanh.
S hành trình piston. A bán kính trụcăkhuỷu. l chiềuădài thanh truyền.
Doăđĩ, điểm chết trên (TDC) bùăđắp do truyền nhiệt cĩ thểthuăđ ợc từph ơngă