ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, VOL 17, NO 1.2, 2019 13 NGHIÊN CỨU BƠM NHIỆT VỚI HỆ THỐNG LẠNH SỬ DỤNG CHU TRÌNH RANKINE HỮU CƠ STUDY ON HEAT PUMP WITH REFRIGERANT SYSTEM USING ORGANIC RANKINE CYCLE PROCESS Hoàng Thành Đạt , Hồ Trần Anh Ngọc Trường Đại học Sư phạm Kỹ thuật – Đại học Đà Nẵng; hoangthanhdat1976@gmail.com, anhngoctr@yahoo.com Tóm tắt - Nghiên cứu sử dụng nguồn nhiệt thừa cấp nhiệt cho hệ thống bơm nhiệt với hệ thống lạnh sử dụng chu trình Rankine hữu ORC (Organic Rankine Cycle) Trên tính tốn theo mơ hình để chọn mơi chất thích hợp dùng cho hệ thống, đưa kết tính tốn hệ số làm nóng làm lạnh COP môi chất R22, R600, R601, R123, R1234ze, R134a, R152a, R227ea, R245fa, R717, Rượu, R718 Nghiên cứu tính năng, đặc tính thay đổi hệ thống bơm nhiệt kết hợp làm lạnh tùy theo thay đổi nhiệt độ thiết bị sinh hơi, nhiệt độ bay nhiệt độ ngưng tụ, kết tính tốn hiệu suất vịng tuần hồn ηorc, hệ số làm việc hệ thống COPS Các kết luận rút từ kết tính tốn hệ thống bơm nhiệt với hệ thống lạnh chu trình Rankine hữu Abstract - The research is on the use of extra heat sources for heat pump systems combined with refrigeration through the ORC (Organic Rankine Cycle) system based on the new model to choose the most suitable for the system The results of calculating the heating and cooling coefficients for COPs for R22, R600, R601, R123, R152a, R227ea, R245fa, R717, Alcohol, R718, and R718 are given By studying the features, characteristics of the change of the heat pump combined cooling according to the temperature change in the steam generator, evaporative temperature and condensing temperature, we see that the results of the calculation of cycle efficiency complete ηorc, the working coefficient of the COPS system The conclusions are drawn from the results of calculating the heat pump system combined cooling through the ORC system Từ khóa - ORC; bơm nhiệt; môi chất lạnh; COP; thu hồi nhiệt Key words - ORC; heat pump; refrigerant; COP; heat recovery Đặt vấn đề Hiện nay, nguồn lượng ngày cạn kiệt, việc nghiên cứu sử dụng nguồn lượng thừa, tiết kiệm nhiều nhà khoa học quan tâm Việc sử dụng nguồn nhiệt thải cao để sản xuất điện nghiên cứu thực thục, nguồn nhiệt thải trung bình thấp nghiên cứu ứng dụng gặp phải khó khăn định số lượng lớn, tiềm sử dụng nhiều nên việc nghiên cứu sử dụng nguồn lượng thải có ý nghĩa to lớn Đã có nhiều kỹ thuật để thu hồi sử dụng nhiệt thải dùng hệ thống tuần hồn mơi chất hữu cơ, sử dụng lượng mặt trời, phương tiện vận tải xe ô tô tàu thuyền, ứng dụng kỹ thuật ống nhiệt, kỹ thuật dùng nhiệt thải để phát điện, kỹ thuật ứng dụng nhiệt thải cách thay đổi thiết bị nhiệt hệ thống bơm nhiệt hệ thống làm lạnh Hệ thống tuần hoàn môi chất hữu dùng loại chất hữu làm mơi chất tuần hồn cho hệ thống, thích hợp cho việc ứng dụng thu hồi nguồn nhiệt thải lượng thấp, ví dụ lượng mặt trời, lượng địa nhiệt, lượng sinh khối Từ thập niên 60 nhà nghiên cứu khoa học quan tâm đến môi chất hữu có nhiệt độ sơi thấp Trung Quốc bắt đầu nghiên cứu vào thập niên 80 kỷ 20 Đại học Wenshou Đài Loan, Hung TC Nghiên cứu phân tích dùng mơi chất benzen (C6H6), C6H5CH3, R123 R113 cho hệ thống tuần hồn mơi chất hữu [1-2] Cơng ty đóng tàu Mitsui Nhật Bản nghiên cứu sử dụng nguồn lượng khói thải lị nhiệt độ 340℃ dùng chu trình ORC để phát điện với cơng suất máy phát 14MW, hiệu suất hoạt động thực tế hiệu suất lượng hữu ích đạt 16,1% 42,8% [3] Công ty chế tạo máy Mỹ dùng môi chất R123 làm mơi chất cho hệ thống tuần hồn hữu cơ, sử dụng lượng thải nhiệt độ 120~220℃ nhà máy lọc dầu, nhà máy hóa học để phát điện công suất 1500~3500kW [4] Demierre nhiều tác giả [5] người dùng môi chất R134a cho hệ thống ORC-VCR với công suất ngưng tụ 20kW, Wang [6-7] nghiên cứu đưa hệ thống ORC-VCR với công suất lạnh 5kW, Bu xianbiao, Li huashan, Wang lingbao [8] năm 2013 nghiên cứu, phân tích sử dụng nguồn nhiệt thải tàu thuyền cấp nhiệt cho chu trình Rankine hữu dùng chạy hệ thống điều hịa khơng khí Để thu hồi nhiệt lượng khí thải động tàu nhiệt độ nước làm mát động cơ, họ dùng chu trình Rankine hữu hệ thống lạnh để làm điều hòa nhiệt độ, xây dựng hệ thống nhiệt động học, phân tích số loại mơi chất thích hợp với tính hệ thống Nhiệt độ sinh nhiệt độ ngưng tụ có ảnh hưởng lớn đến tính hệ thống Thơng qua thay đổi lưu lượng nước nóng điều khiển điều tiết nhiệt độ nước nóng từ ưu việt hóa tính hệ thống Hệ thống kết hợp tuần hồn mơi chất hữu hệ thống lạnh sử dụng chung loại môi chất lạnh, dùng chung hệ thống ngưng tụ, mặt kết cấu đơn giản, mặt khác tượng rị rỉ mơi chất ảnh hưởng đến hệ thống khơng nhiều Ở nội dung xây dựng mơ hình thệ thống bơm nhiệt dùng chung môi chất cho chu trình Rankine hữu vịng tuần hồn làm lạnh Tiến hành nghiên cứu tính ảnh hưởng hệ thống, thay đổi nhiệt độ sinh hơi, nhiệt độ bay hơi, nhiệt độ ngưng tụ… Kết tính tốn làm tư liệu để chế tạo hệ thống, bước nghiên cứu Nguyên lý làm việc hệ thống Hình biểu thị hệ thống bơm nhiệt chu trình Rankine hữu với hệ thống lạnh, hệ thống gồm có thiết bị Hồng Thành Đạt , Hồ Trần Anh Ngọc 14 sau: I - Thiết bị sinh hơi, II – Máy dãn nở, III – Thiết bị ngưng tụ, IV – Bơm môi chất, V – Van tiết lưu, VI – Thiết bị bay hơi, VII – Máy nén lạnh Hệ thống gồm hai phận hợp thành gồm chu trình Rankine hữu vịng tuần hồn máy lạnh, hai phận thông qua máy dãn nở Máy nén lạnh hợp thành Chu trình Rankine hữu bao gồm: Thiết bị sinh hơi, máy dãn nở, thiết bị ngưng tụ, bơm mơi chất; Vịng tuần hồn máy lạnh bao gồm: Thiết bị bay hơi, máy nén lạnh, thiết bị ngưng tụ, van tiết lưu, hai vịng tuần hồn dùng chung thiết bị ngưng tụ Quá trình hoạt động sau: Môi chất thiết bị sinh nhận nhiệt bay trở thành có áp suất nhiệt độ cao đưa đến máy dãn nở, dãn nở sinh công áp suất nhiệt độ giảm xuống, khỏi máy dãn nở môi chất đưa vào thiết bị ngưng tụ nhả nhiệt cho môi trường làm mát ngưng tụ thành lỏng bơm môi chất bơm lên áp suất cao đưa lại thiết bị sinh nhận nhiệt sinh tiếp tục chu trình (hồn thành vịng tuần hồn thứ nhất) Vịng tuần hồn hệ thống lạnh (vịng tuần hồn thứ hai): Tại thiết bị ngưng tụ, lỏng môi chất khỏi thiết bị chia thành hai phần, phần qua bơm môi chất phần qua van tiết lưu giảm áp suất nhiệt độ sau vào thiết bị bay hấp thụ nhiệt môi trường làm lạnh bay có áp suất nhiệt độ thấp đưa máy nén, nén đoạn nhiệt máy nén lên áp suất nhiệt độ cao sau vào thiết bị ngưng tụ nhả nhiệt cho môi trường làm mát ngưng tụ thành lỏng, lỏng ngưng tụ khỏi thiết bị ngưng tụ tiếp tục chu trình II III Hình Sơ đồ nguyên lý hệ thống bơm nhiệt kết hợp máy lạnh dùng chu trình Rankine hữu Hình biểu diễn đồ thị T-s lgp-h hệ thống bơm nhiệt ORC kết hợp làm lạnh T 5' mvcr h1 h2 morc h7 h6 Hệ số làm lạnh COP: COP 4s 6' 2s h6 h5 mvcr h1 h4 morc (h6 h5 )(h1 h2 ) (h1 h4 )(h7 h6 ) 7s (3) Hệ số làm nóng COP: COP (h2 h3 )morc (h7 h8 )mvcr (h1 h4 )morc (h h ) (h h )(h h ) (h1 h4 ) (h1 h4 )(h7 h6 ) (4) khỏi bơm, Wpump : công suất bơm Pb morc h1 h4 5' 6' (7) 2s 7s orc S h Hình Đồ T-s lgp-h hệ thống bơm nhiệt kết hợp làm lạnh chu trình Rankine hữu 1→2→3→4→5’→6’→1 biểu thị vịng tuần hồn hữu cơ, 1→2s biểu thị trình dãn nở đẳng entropy, q trình 3→4s bơm mơi chất đẳng enthalpy, 7→2s→3→5→6→7 biểu thị chu trình máy lạnh, 6→7s trình nén đoạn nhiệt (8) Với: Wnet: công thực sinh chu trình Rankine hữu (2) Wnet Wexp Wpump Với morc: lưu lượng môi chất chu trình Rankine hữu cơ, mvcr: lưu lượng mơi chất vịng tuần hồn làm lạnh, h1, h2: enthalpy vào máy dãn nở, h6, h7: enthalpy vào khỏi máy nén Tỉ số lưu lượng: Với: Pb: công suất cấp nhiệt cho thiết bị sinh lgp 4s (1) h h Wpump morc 4s (6) pump Với: pump : công suất bơm h4s: đẳng enthalpy V IV h1 h2 morc h7 h6 mvcr Công thức trên: Wexp: công suất máy dãn nở sinh h2s đẳng enthalpy khỏi máy dãn nở VI I Tính tốn mơ hình Cơng suất máy dãn nở sinh cung cấp hoàn toàn cho máy nén lạnh hoạt động ta có cơng thức sau: Trong cơng thức h3: enthalpy vào bơm môi chất, h5: enthalpy vào dàn lạnh, h4: enthalpy khỏi bơm môi chất - Tính cho chu trình Rankine hữu (5) Wexp morc h1 h2s VII máy nén lạnh, 7→8 trình ngưng tụ đẳng áp thiết bị ngưng tụ, 8→5 trình tiết lưu, 5→6 trình hấp thu nhiệt thiết bị bay Wnet Pb Với: ηorc: hiệu suất chu trình Rankine hữu W h h Wm net h1 h2s exp 4s pump morc (9) (10) Với: Wm: với tỉ số Wnet: ηexp, morc.: hiệu suất, đẳng enthalpy máy dãn nở - Tính cho vịng tuần hồn lạnh Wcom Wexp (11) ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, VOL 17, NO 1.2, 2019 Công thức Wcom: công suất máy lạnh: Wcom mvcr h7s h6 (12) Công thức h7s: đẳng enthalpy khỏi máy nén Hệ số làm lạnh máy lạnh: Qeva h h Wcom h7s h6 (13) Công thức trên,Qeva: công suất làm lạnh Tỉ số áp suất vào máy nén R P7s P6 (14) Công thức trên,P7s: áp suất khỏi máy nén;P6: áp suất vào máy nén, CR COPvcr R (15) Công thức trên,CR tỉ số COPvcr R Hệ số COPs toàn hệ thống: COPs orc COPvcr Hình biểu thị mối quan hệ loại môi chất với hệ số làm lạnh COP tùy theo nhiệt độ thay đổi Từ Hình ta thấy, nhiệt độ sinh tăng lên hệ số làm lạnh COP tất môi chất tăng theo Tại nhiệt độ sinh môi chất có hệ số làm lạnh COP khác Mơi chất R717, Rượu R718 đạt hệ số làm lạnh COP cao Môi chất R134a,R227ea R245fa đạt hệ số làm lạnh thấp cụ thể nhiệt độ sinh 91℃, COP R718 đạt trung bình 0,8 R134a đạt trung bình 0,66 Hình thể tùy theo nhiệt độ sinh khác mà hệ số làm nóng COP mơi chất khác Từ Hình ta nhận thấy, nhiệt độ sinh tăng cao hệ số làm nóng COP môi chất tăng cao Tại nhiệt độ sinh khác hệ số làm nóng COP môi chất khác Các môi chất R717, Rượu, R718 có hệ số làm nóng COP cao tương đương Các môi chất R227ea, R1234ze R134a có hệ số làm nóng COP thấp Cụ thể nhiệt độ thiết bị sinh 95℃ hệ số làm nóng mơi chất R718 1,8 R227ea nhỏ 1,58 Chọn thông số làm việc:Nhiệt độ thiết bị sinh 55℃~95℃ Nhiệt độ ngưng tụ 35℃, 40℃, 45℃ 50℃ nhiệt độ bay 1℃~7℃ Máy dãn nở máy nén lạnh, bơm môi chất với hiệu suất đẳng enthalpy 0,85, 0,80 0,90 Nhiệt lượng cung cấp cho thiết bị sinh Pb = 2kW Chọn mơi chất tính toán: R22, R600, R601, R123 R1234ze, R134a, R152a, R227ea, R245fa, R717, Rượu, R718 Dùng phần mềm NIST để xác định tính vật lý như: Nhiệt độ, áp suất, enthalpy, entropy, thể tích thơng số khác Giả định: (1) Hệ thống hoạt động ổn định; (2) Không tính đến tổn thất áp suất, tổn thất nhiệt đường ống thiết bị trao đổi nhiệt; (3) Q trình tiết lưu đẳng enthalpy Tính tốn phân tích kết 4.1 Hệ số làm nóng lạnh COP môi chất 0.9 t1=55℃ t1=61℃ t1=67℃ 0.8 t1=79℃ t1=85℃ t1=91℃ 1.5 1.4 1.3 1.2 1.1 Hình Hệ số làm nóng COP với mơi chất khác 4.2 Các thông số ảnh hưởng đến hệ thống Hình biểu thị ảnh hưởng nhiệt độ bay hiệu suất chu trình Rankine hữu Từ Hình ta biết hiệu suất chu trình Rankine hữu ηorc tăng lên nhiệt độ thiết bị bay tăng lên Nhiệt độ sinh tăng từ 55℃ đến 95℃ hiệu suất chu trình Rankine hữu tăng từ 0,01 đến 0,11 0.12 R1234ze 0.08 ηorc Hệ số làm lạnh COP R123 0.10 R134a R152a 0.06 0.4 0.04 0.3 0.02 Hình Hệ số làm lạnh COP với mơi chất khác t1=73℃ 1.6 0.6 0.1 t1=67℃ t1=91℃ 1.7 0.7 0.2 t1=61℃ t1=85℃ 1.8 t1=73℃ 0.5 t1=55℃ t1=79℃ 1.9 (16) Hệ số làm nóng COP COPvcr 15 R227ea R22 0.00 55 65 75 tb/℃ 85 95 Hình Ảnh hưởng nhiệt độ sinh hiệu suất chu trình Rankine hữu Hồng Thành Đạt , Hồ Trần Anh Ngọc 16 Hình biểu thị Wm với tỉ số Wnet morc tùy thuộc vào nhiệt độ sinh Từ hình ta nhận thấy Wm tăng lên nhiệt độ bay tăng lên Cụ thể môi chất R152a nhiệt độ sinh 55℃ Wm=2,5 nhiệt độ sinh 95℃ Wm= 23, tăng trung bình gấp 10 lần Nguyên nhân sau: từ công thức (10) nhận thấy nhiệt độ sinh tăng lên cao dường nhiệt độ ngưng tụ không đổi, h1 tăng với h2s khoảng cách lớn h4s h3 thay đổi không lớn dẫn đến Wm tăng nhiều nhiệt độ ngưng tụ tăng cao Cụ thể sau: Tại nhiệt độ ngưng tụ 35℃ môi chất R123 R152a R22.R134a R1234ze R227ea có COPvcr phân biệt 9,07, 8,84 8,60, 8,60, 8,58 8,17 Nhưng nhiệt độ ngưng tụ 50℃ môi chất có COPvcr giảm 5,52, 5,30, 5,08, 5,02, 4,99 4,51 0.14 R123 R134a R227ea 0.13 R1234ze R152a R22 0.12 25 Wm/kW.s.kg-1 R123 15 10 ηorc 0.11 20 0.10 R1234ze 0.09 R134a 0.08 R152a 0.07 R227ea 0.06 35 R22 55 65 75 tb/℃ 85 Hình Ảnh hưởng nhiệt độ ngưng tụ hệ số ηorc 95 10 Hình biểu thị hệ số làm việc hệ thống COPs tùy theo nhiệt độ sinh Đối với tất mơi chất hệ số COPs tùy theo nhiệt độ tb tăng cao tăng cao, R227ea có hệ số COPs thấp nhất, nhiệt độ thiết bị sinh tăng lên 1℃ hệ số COPs tăng bình qn 2% Thứ tự mơi chất có COPs từ cao đến thấp sau: R123>R152a>R22>R134a>R1234ze>R227ea COPvcr Hình Ảnh hưởng nhiệt độ sinh Wm R123 R134a R227ea 35 0.70 R1234ze R152a R22 0.80 40 45 50 Nhiệt độ ngưng tụ (℃) Hình Ảnh hưởng nhiệt độ ngưng tụ hệ số COPvcr 0.60 R123 0.50 R1234ze 0.40 R134a 0.30 R152a 0.20 R227ea 0.10 R22 0.00 55 65 75 85 95 tb/℃ Hình Ảnh hưởng nhiệt độ sinh hệ số COPs nhiệt độ ngưng tụ hệ số làm việc hệ thống Hình biểu thị mối quan hệ hiệu suất ηorc chu trình Rankine hữu nhiệt độ ngưng tụ Từ Hình ta nhận thấy nhiệt độ ngưng tụ tăng lên ηorc giảm xuống Cụ thể nhiệt độ 35℃ môi chất R123, R152a, R22, R134a, R1234ze R227ea ηorc có giá trị tương ứng: 0,1285, 0,1238, 0,1239, 0,1142, 0,1112 0,0973 nhiệt độ ngưng tụ 50℃ mơi chất có ηorc giảm sau: 0,0958, 0,0895m, 0,0879, 0,0816, 0,0798 0,0687 Hình biểu thị mối quan hệ hệ số làm lạnh COPvcr thay đổi nhiệt độ ngưng tụ Từ Hình ta nhận thấy, tất mơi chất lạnh có COPvcr giảm xuống COPs COPs 50 45 40 Nhiệt độ ngưng tụ (℃) 1.2 1.1 1.0 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 R123 R134a R227ea R1234ze R152a R22 35 40 45 Nhiệt độ ngưng tụ (℃) 50 Hình 10 Ảnh hưởng nhiệt độ ngưng tụ hệ số COPs Hình 10 biểu thị mối quan hệ hệ số làm lạnh toàn hệ thống COPs thay đổi nhiệt độ ngưng tụ Từ Hình 10 ta thấy, tất mơi chất có hệ số làm việc COPs tùy theo nhiệt độ ngưng tụ tăng lên giảm xuống Cụ thể nhiệt độ ngưng tụ từ 35℃ đến 55℃ mơi chất R123 có COPs có giá trị 1,18 giảm xuống 0,55 tỉ lệ trung bình giảm 53,3% 4.3 Ảnh hưởng nhiệt độ bay hệ số làm lạnh chu trình COPvcr Hình 11 biểu thị mối mối quan hệ COPvcr với thay đổi nhiệt độ bay Từ Hình 11 ta thấy, tất ISSN 1859-1531 - TẠP CHÍ KHOA HỌC VÀ CƠNG NGHỆ ĐẠI HỌC ĐÀ NẴNG, VOL 17, NO 1.2, 2019 môi chất hệ số COPvcr tăng lên tùy theo nhiệt độ bay tăng cao Cụ thể sau: Tại nhiệt độ bay 1℃ môi chất R123, R152a, R22, R134a, R1234ze R227ea COPvcr phân biệt 5,28, 5,10, 4,91, 4,85, 4,81 4,37 Nhưng nhiệt độ bay 7℃ mơi chất có hệ số COPvcr có giá trị trình tự tăng sau: 6,40, 6,18, 5,96, 5,91, 5,89 5,43 6.5 COPvcr 6.0 R123 R1234ze R134a R152a R227ea R22 5.5 5.0 4.5 4.0 Nhiệt độ bay (℃) Hình 11 Ảnh hưởng nhiệt độ bay hệ số COPvcr 0.75 0.70 0.65 R123 R134a R227ea R1234ze R152a R22 17 chất COPs có giá trị tăng lên sau: 0,6843, 0,6251, 0,5972, 0,5484, 0,5327 0,4270 Kết luận Ở nội dung đề suất nguyên lý hệ thống bơm nhiệt tuần hoàn hữu kết hợp làm lạnh đồng thời xây dựng tính tốn tính chất mơ hình Căn việc tính tốn mơ hình hệ thống nghiên cứu đặc tính biến đổi loại môi chất khác tùy theo nhiệt độ sinh hơi, nhiệt độ bay hơi, nhiệt độ ngưng tụ, rút kết sau đây: 1) Đối với hệ số làm lạnh làm nóng COP R717, Rượu R718 đạt hệ số cao hơn, cịn mơi chất R227ea, R245fa với R134a đạt hệ số thấp 2) Nhiệt độ thiết bị sinh ảnh hưởng lớn tổng hiệu hệ thống, nhiệt độ sinh tăng 1℃ tổng hiệu hệ thống tăng lên trung bình 2% 3) Mức độ hiệu trung bình mơi chất tăng lên giảm xuống tùy thuộc vào nhiệt độ bay tăng lên hay giảm xuống Do tính tốn thiết kế hệ thống dùng chung mơi chất tính tốn hệ thống bơm nhiệt kết hợp hai vịng tuần hồn nên nâng cao nhiệt độ sinh hơi, giảm nhiệt độ ngưng tụ, nâng cao nhiệt độ bay COPs 0.60 TÀI LIỆU THAM KHẢO 0.55 0.50 0.45 0.40 0.35 0.30 Nhiệt độ bay (℃) Hình 12 Ảnh hưởng nhiệt độ bay hệ số COPs Hình 12 biểu thị mối quan hệ nhiệt độ bay với hệ số làm việc tồn hệ thống COPs Từ Hình 12 ta nhận thấy, tất mơi chất có nhiệt độ bay tăng lên hệ số COPs tăng lên Tại nhiệt độ bay 1℃ môi chất R123, R152a, R22, R134a, R1234ze R227ea COPs giá trị 0,5647, 0,5158, 0,4927, 0,4498, 0,4355 0,3435 Nhưng nhiệt độ bay 7℃ với môi [1] Hung T.C,Shai T.Y,Wang S.K A review of organic rankine cycles (ORCs) for the recovery of low-grade waste heat [J] . Energy, 1997, 22(7): 66-667. [2] Hung T.C Waste heat recovery oforganic rankine cycle using dryfluids Energy Conversion and Management 2001, 42(5): 539-553. [3] Hirakawa Y 14MW ORC plant installed at Nippon steel [J] 1981 [4] Legmann H Recovery of industrial heat in the cement industry by means of the ORC process [C] Cement Industry Technical Confernece IEEE-IAS/PCA 44th IEEE 2002: 29-35 [5] LIM S M Economies of ship size: A new evaluation [J] Maritime Policy anh Managenment 1994 21 (2):149-166 [6] RYDER S C CHAPPELL D Optimal speed anh ship size for liner trades[J] Maritime Policy anh Management 1980, (1): 55-57 [7] TALLEY W K Optimal Container ship size [J] Maritime Policy anh Management 1980, 17 (3): 165-175 [8] 卜宪标,李华山,王令宝等 船舶余热驱动的有机朗肯-蒸汽压 缩空调性能分析与工质选择 大连海事大学学报 2013, 4:101-103 (BBT nhận bài: 19/9/2018, hoàn tất thủ tục phản biện: 30/10/2018)