1. Áp suất chưa cho biết đầy đủ năng lượng
Hai hệ có cùng áp suất nhưng thể tích khí khác nhau có thể chứa mức năng lượng có khả năng sinh công rất khác nhau. Ngay cả khi áp suất và thể tích giống nhau, cách gọi “năng lượng tích trữ” vẫn cần một định nghĩa: công giãn nở theo đường nào, có nhận nhiệt từ môi trường không và phần công đẩy khí quyển có được tính là hữu ích không?
Tài liệu nhiệt động học của MIT nhấn mạnh công phụ thuộc đường quá trình; tích phân áp suất trong hệ theo thể tích chỉ biểu diễn công biên khi điều kiện gần tĩnh phù hợp.[1] Bài này lập các mô hình lý tưởng để so sánh định nghĩa. Không dùng kết quả làm năng lượng nổ, khoảng cách an toàn hoặc quy trình xả áp.
2. Bộ đầu vào và quy ước dấu
Giả định một lượng không khí lý tưởng ở trạng thái đầu có thể tích V_1 = 0,300 m³, áp suất dư p_g = 9,00 bar, áp suất môi trường p_0 = 1,00 bar tuyệt đối và nhiệt độ T_1 = T_0 = 300 K. Chọn R = 287 J/(kg·K), γ = 1,40 không đổi cho bài tập.
p_1 = p_g + p_0 = 10,00 bar tuyệt đối = 1.000.000 Pa
m = p_1 V_1/(R T_1) = 3,484321 kg
1 bar = 100.000 Pa; 1 Pa·m³ = 1 J
Công do khí thực hiện ra ngoài được lấy dương. Các tỷ số áp suất và phương trình trạng thái phải dùng áp suất tuyệt đối. T_1 dùng K, không đưa nhiệt độ theo °C trực tiếp vào pV = mRT. Thể tích 0,300 m³ là thể tích khí thực trong trạng thái đầu, không phải thể tích khí đã quy đổi về điều kiện tham chiếu.
3. Giãn nở đẳng nhiệt thuận nghịch
Hình dung lượng khí được đặt trong một hệ piston lý tưởng, giữ nguyên khối lượng và nhận nhiệt để luôn ở 300 K trong khi giãn nở từ p_1 xuống p_0. Từ pV = p_1V_1:
V_2,iso = V_1 p_1/p_0 = 3,000 m³
W_iso = ∫ p dV = p_1 V_1 ln(V_2,iso/V_1)
W_iso = p_1 V_1 ln(p_1/p_0) = 690,775528 kJ
Khí lý tưởng có nội năng phụ thuộc nhiệt độ, nên ΔU = 0 trong quá trình này. Cân bằng năng lượng cho nhiệt nhận vào bằng công phát ra. Vì vậy, 690,776 kJ không thể được diễn giải là toàn bộ năng lượng ban đầu trong khí tự chuyển thành công mà không có trao đổi nhiệt.
MIT trình bày dạng logarit này cho giãn nở đẳng nhiệt gần tĩnh.[1] Các giá trị số và trường hợp đang xét được tự xây dựng ở đây. Một phép xả nhanh từ bình cứng không có cùng ranh giới và đường quá trình với piston đẳng nhiệt.
4. Tách công đẩy môi trường
Một phần công biên dùng để chiếm thêm thể tích chống lại áp suất môi trường p_0. Khi p_0 không đổi, phần đó bằng p_0(V_2−V_1). Bài giảng của IIT Kanpur phân biệt công cực đại với phần công hữu ích còn lại sau khi trừ công đẩy môi trường.[2]
W_khí_quyển,iso = p_0(3,000 − 0,300) = 270,000 kJ
W_hữu_ích,iso = W_iso − W_khí_quyển,iso
W_hữu_ích,iso = 420,775528 kJ
Với trạng thái đầu có T_1 = T_0, mô hình khí lý tưởng và quá trình thuận nghịch đưa hệ về cân bằng với môi trường, đây là phần khả năng sinh công hữu ích của hệ kín theo các giả thiết đã chọn. Dạng gọn là p_1V_1 ln(p_1/p_0) − (p_1−p_0)V_1.
Tích áp suất dư–thể tích p_gV_1 cũng bằng 270 kJ trong ví dụ, nhưng không bằng công giãn nở và không bằng công hữu ích vừa tính. Cùng đơn vị J không làm hai đại lượng trở thành cùng một ý nghĩa vật lý. Cần ghi tên đại lượng ngay cạnh kết quả trong bảng tính.
5. Giãn nở đoạn nhiệt thuận nghịch
Xét một đường khác: hệ kín không trao đổi nhiệt, quá trình thuận nghịch và γ không đổi. NASA suy ra quan hệ áp suất–thể tích và nhiệt độ cho quá trình đẳng entropy của khí lý tưởng.[3] Khi kết thúc ở cùng p_0:
pV^γ = hằng số
V_2,ad = V_1(p_1/p_0)^(1/γ) = 1,553842 m³
T_2,ad = T_1(p_0/p_1)^((γ−1)/γ) = 155,384240 K
W_ad = (p_1V_1 − p_0V_2,ad)/(γ−1) = 361,539399 kJ
Công đẩy môi trường ở đường này bằng 125,384240 kJ; phần còn lại bằng 236,155159 kJ. Hai kết quả công hữu ích khác nhau vì đường nhiệt động và trạng thái cuối về nhiệt độ khác nhau. Khí sau đường đoạn nhiệt vẫn lạnh hơn môi trường, chưa ở trạng thái cân bằng nhiệt cuối cùng dùng trong định nghĩa khả năng sinh công đầy đủ.
Nhiệt độ 155,38 K là kết quả của mô hình γ hằng và không trao đổi nhiệt. Nó không phải dự báo chắc chắn nhiệt độ khí hoặc thành bình của một thiết bị thật. Nếu tính chất biến đổi, có ngưng tụ hoặc trao đổi nhiệt đáng kể, cần dùng mô hình tính chất và quá trình phù hợp.
6. Kiểm độc lập bằng nội năng
Từ các giả thiết khí lý tưởng, c_v = R/(γ−1) = 717,5 J/(kg·K). Với đường đoạn nhiệt Q = 0, định luật thứ nhất cho W_ad = m c_v(T_1−T_2). Thay số thu lại 361,539399 kJ, giống kết quả tích phân p(V).
Hai phép tính dùng các biểu thức khác nhau nhưng cùng giả thiết. Sự trùng khớp là kiểm tra đại số và đơn vị, chưa xác nhận mô hình mô tả được sự cố thực tế. Nếu phần mềm dùng nhiệt dung biến đổi còn bảng tính dùng γ hằng, sai khác kết quả có thể bắt nguồn từ giả thiết vật tính chứ không phải lỗi lập trình.
7. Độ nhạy với áp suất và thể tích bị bỏ sót
Nếu nhầm 9 bar dư thành 9 bar tuyệt đối nhưng vẫn dùng p_0 = 1 bar, công hữu ích đẳng nhiệt tính sai thành 353,250636 kJ, thấp hơn 16,048%. Lỗi chỉ một bar ở đầu vào đã đi vào cả hệ số trước logarit lẫn tỷ số trong logarit.
Giữ p_0 và V_1, đạo hàm của công hữu ích đẳng nhiệt theo p_1 là V_1 ln(p_1/p_0). Tại bộ số đang xét, sai lệch nhỏ ±0,20 bar tương ứng thay đổi tuyến tính khoảng ±13,816 kJ. Đây là độ nhạy cục bộ; với thay đổi lớn nên tính lại toàn công thức.
Nếu hệ còn 50 m ống có đường kính trong 25 mm chứa cùng trạng thái khí, thể tích ống bằng π × 0,025² × 50/4 = 0,024544 m³. Phần này làm tổng thể tích tăng 8,181% so với bình 0,300 m³. Trong mô hình cùng p và T, cả các công đang xét tăng cùng tỷ lệ. Nhưng nếu van cô lập chia các khoang ở áp suất khác nhau, cần lập riêng từng miền thay vì cộng mọi thể tích rồi dùng một áp suất chung.
8. Bình cứng xả khí không phải piston kín
Một bình cứng có thể tích hình học gần như không đổi nhưng khối lượng bên trong giảm khi xả. Với thể tích kiểm soát bao quanh bình, năng lượng đi ra theo dòng khí, có thể trao đổi nhiệt qua thành; công biên của thành cứng gần bằng không. Kết luận “dV = 0 nên hệ không còn năng lượng nguy hiểm” là sai vì đã bỏ mất dòng năng lượng qua cửa ra.
Ngược lại, gán toàn bộ lượng khí cho một piston giãn nở thuận nghịch rồi gọi kết quả là năng lượng sóng áp cũng không đúng. Sự cố còn liên quan quá trình phá hủy, dòng không ổn định và phân bố năng lượng vào nhiều dạng. Bài tính hiện tại không mô hình hóa các cơ chế đó.
9. Đọc một con số năng lượng trong hồ sơ áp lực
Hướng dẫn nội bộ về an toàn áp lực của Berkeley Lab xem áp suất, thể tích, độ phức tạp của hệ và hậu quả giải phóng ngoài kiểm soát là những yếu tố cần xét; tài liệu cũng có một cách tính năng lượng giãn nở đẳng entropy.[4] Các ngưỡng phân loại của phòng thí nghiệm không được áp dụng mặc nhiên cho nhà máy hoặc công trường khác.
Một kết quả có thể kiểm toán nên đi kèm sơ đồ miền khí, áp suất tuyệt đối, nhiệt độ, thể tích thực, vật tính, đường quá trình và định nghĩa công. Đặc biệt, phải ghi rõ đã trừ công đẩy môi trường hay chưa. Nhờ đó người đọc có thể so sánh các phép tính cùng cơ sở, thay vì chọn công thức cho con số nhỏ nhất hoặc xem mọi giá trị có đơn vị kJ là tương đương.
Nguồn tham khảo
[1] Thermodynamics and Propulsion, 1.3 — Changing the State of a System with Heat and Work — Massachusetts Institute of Technology. Truy cập ngày 09/11/2026.
[2] Lecture 26 — Availability Function for a non-flow Process — Indian Institute of Technology Kanpur. Truy cập ngày 09/11/2026.
[3] Isentropic Compression — NASA Glenn Research Center. Truy cập ngày 09/11/2026.
[4] ES&H Manual, Chapter 07 — Pressure Safety, appendices E–F — Lawrence Berkeley National Laboratory, revised 07/2024. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Hai hệ có cùng áp suất nhưng thể tích khí khác nhau có thể chứa mức năng lượng có khả năng sinh công rất khác nhau. Ngay cả khi áp suất và thể tích giống nhau, cách gọi “năng lượng tích trữ” vẫn cần một định nghĩa: công giãn nở theo đường nào, có nhận nhiệt từ môi trường không và phần công đẩy khí quyển có được tính là hữu ích không?
Tài liệu nhiệt động học của MIT nhấn mạnh công phụ thuộc đường quá trình; tích phân áp suất trong hệ theo thể tích chỉ biểu diễn công biên khi điều kiện gần tĩnh phù hợp.[1] Bài này lập các mô hình lý tưởng để so sánh định nghĩa. Không dùng kết quả làm năng lượng nổ, khoảng cách an toàn hoặc quy trình xả áp.
2. Bộ đầu vào và quy ước dấu
Giả định một lượng không khí lý tưởng ở trạng thái đầu có thể tích V_1 = 0,300 m³, áp suất dư p_g = 9,00 bar, áp suất môi trường p_0 = 1,00 bar tuyệt đối và nhiệt độ T_1 = T_0 = 300 K. Chọn R = 287 J/(kg·K), γ = 1,40 không đổi cho bài tập.
p_1 = p_g + p_0 = 10,00 bar tuyệt đối = 1.000.000 Pa
m = p_1 V_1/(R T_1) = 3,484321 kg
1 bar = 100.000 Pa; 1 Pa·m³ = 1 J
Công do khí thực hiện ra ngoài được lấy dương. Các tỷ số áp suất và phương trình trạng thái phải dùng áp suất tuyệt đối. T_1 dùng K, không đưa nhiệt độ theo °C trực tiếp vào pV = mRT. Thể tích 0,300 m³ là thể tích khí thực trong trạng thái đầu, không phải thể tích khí đã quy đổi về điều kiện tham chiếu.
3. Giãn nở đẳng nhiệt thuận nghịch
Hình dung lượng khí được đặt trong một hệ piston lý tưởng, giữ nguyên khối lượng và nhận nhiệt để luôn ở 300 K trong khi giãn nở từ p_1 xuống p_0. Từ pV = p_1V_1:
V_2,iso = V_1 p_1/p_0 = 3,000 m³
W_iso = ∫ p dV = p_1 V_1 ln(V_2,iso/V_1)
W_iso = p_1 V_1 ln(p_1/p_0) = 690,775528 kJ
Khí lý tưởng có nội năng phụ thuộc nhiệt độ, nên ΔU = 0 trong quá trình này. Cân bằng năng lượng cho nhiệt nhận vào bằng công phát ra. Vì vậy, 690,776 kJ không thể được diễn giải là toàn bộ năng lượng ban đầu trong khí tự chuyển thành công mà không có trao đổi nhiệt.
MIT trình bày dạng logarit này cho giãn nở đẳng nhiệt gần tĩnh.[1] Các giá trị số và trường hợp đang xét được tự xây dựng ở đây. Một phép xả nhanh từ bình cứng không có cùng ranh giới và đường quá trình với piston đẳng nhiệt.
4. Tách công đẩy môi trường
Một phần công biên dùng để chiếm thêm thể tích chống lại áp suất môi trường p_0. Khi p_0 không đổi, phần đó bằng p_0(V_2−V_1). Bài giảng của IIT Kanpur phân biệt công cực đại với phần công hữu ích còn lại sau khi trừ công đẩy môi trường.[2]
W_khí_quyển,iso = p_0(3,000 − 0,300) = 270,000 kJ
W_hữu_ích,iso = W_iso − W_khí_quyển,iso
W_hữu_ích,iso = 420,775528 kJ
Với trạng thái đầu có T_1 = T_0, mô hình khí lý tưởng và quá trình thuận nghịch đưa hệ về cân bằng với môi trường, đây là phần khả năng sinh công hữu ích của hệ kín theo các giả thiết đã chọn. Dạng gọn là p_1V_1 ln(p_1/p_0) − (p_1−p_0)V_1.
Tích áp suất dư–thể tích p_gV_1 cũng bằng 270 kJ trong ví dụ, nhưng không bằng công giãn nở và không bằng công hữu ích vừa tính. Cùng đơn vị J không làm hai đại lượng trở thành cùng một ý nghĩa vật lý. Cần ghi tên đại lượng ngay cạnh kết quả trong bảng tính.
5. Giãn nở đoạn nhiệt thuận nghịch
Xét một đường khác: hệ kín không trao đổi nhiệt, quá trình thuận nghịch và γ không đổi. NASA suy ra quan hệ áp suất–thể tích và nhiệt độ cho quá trình đẳng entropy của khí lý tưởng.[3] Khi kết thúc ở cùng p_0:
pV^γ = hằng số
V_2,ad = V_1(p_1/p_0)^(1/γ) = 1,553842 m³
T_2,ad = T_1(p_0/p_1)^((γ−1)/γ) = 155,384240 K
W_ad = (p_1V_1 − p_0V_2,ad)/(γ−1) = 361,539399 kJ
Công đẩy môi trường ở đường này bằng 125,384240 kJ; phần còn lại bằng 236,155159 kJ. Hai kết quả công hữu ích khác nhau vì đường nhiệt động và trạng thái cuối về nhiệt độ khác nhau. Khí sau đường đoạn nhiệt vẫn lạnh hơn môi trường, chưa ở trạng thái cân bằng nhiệt cuối cùng dùng trong định nghĩa khả năng sinh công đầy đủ.
Nhiệt độ 155,38 K là kết quả của mô hình γ hằng và không trao đổi nhiệt. Nó không phải dự báo chắc chắn nhiệt độ khí hoặc thành bình của một thiết bị thật. Nếu tính chất biến đổi, có ngưng tụ hoặc trao đổi nhiệt đáng kể, cần dùng mô hình tính chất và quá trình phù hợp.
6. Kiểm độc lập bằng nội năng
Từ các giả thiết khí lý tưởng, c_v = R/(γ−1) = 717,5 J/(kg·K). Với đường đoạn nhiệt Q = 0, định luật thứ nhất cho W_ad = m c_v(T_1−T_2). Thay số thu lại 361,539399 kJ, giống kết quả tích phân p(V).
Hai phép tính dùng các biểu thức khác nhau nhưng cùng giả thiết. Sự trùng khớp là kiểm tra đại số và đơn vị, chưa xác nhận mô hình mô tả được sự cố thực tế. Nếu phần mềm dùng nhiệt dung biến đổi còn bảng tính dùng γ hằng, sai khác kết quả có thể bắt nguồn từ giả thiết vật tính chứ không phải lỗi lập trình.
7. Độ nhạy với áp suất và thể tích bị bỏ sót
Nếu nhầm 9 bar dư thành 9 bar tuyệt đối nhưng vẫn dùng p_0 = 1 bar, công hữu ích đẳng nhiệt tính sai thành 353,250636 kJ, thấp hơn 16,048%. Lỗi chỉ một bar ở đầu vào đã đi vào cả hệ số trước logarit lẫn tỷ số trong logarit.
Giữ p_0 và V_1, đạo hàm của công hữu ích đẳng nhiệt theo p_1 là V_1 ln(p_1/p_0). Tại bộ số đang xét, sai lệch nhỏ ±0,20 bar tương ứng thay đổi tuyến tính khoảng ±13,816 kJ. Đây là độ nhạy cục bộ; với thay đổi lớn nên tính lại toàn công thức.
Nếu hệ còn 50 m ống có đường kính trong 25 mm chứa cùng trạng thái khí, thể tích ống bằng π × 0,025² × 50/4 = 0,024544 m³. Phần này làm tổng thể tích tăng 8,181% so với bình 0,300 m³. Trong mô hình cùng p và T, cả các công đang xét tăng cùng tỷ lệ. Nhưng nếu van cô lập chia các khoang ở áp suất khác nhau, cần lập riêng từng miền thay vì cộng mọi thể tích rồi dùng một áp suất chung.
8. Bình cứng xả khí không phải piston kín
Một bình cứng có thể tích hình học gần như không đổi nhưng khối lượng bên trong giảm khi xả. Với thể tích kiểm soát bao quanh bình, năng lượng đi ra theo dòng khí, có thể trao đổi nhiệt qua thành; công biên của thành cứng gần bằng không. Kết luận “dV = 0 nên hệ không còn năng lượng nguy hiểm” là sai vì đã bỏ mất dòng năng lượng qua cửa ra.
Ngược lại, gán toàn bộ lượng khí cho một piston giãn nở thuận nghịch rồi gọi kết quả là năng lượng sóng áp cũng không đúng. Sự cố còn liên quan quá trình phá hủy, dòng không ổn định và phân bố năng lượng vào nhiều dạng. Bài tính hiện tại không mô hình hóa các cơ chế đó.
9. Đọc một con số năng lượng trong hồ sơ áp lực
Hướng dẫn nội bộ về an toàn áp lực của Berkeley Lab xem áp suất, thể tích, độ phức tạp của hệ và hậu quả giải phóng ngoài kiểm soát là những yếu tố cần xét; tài liệu cũng có một cách tính năng lượng giãn nở đẳng entropy.[4] Các ngưỡng phân loại của phòng thí nghiệm không được áp dụng mặc nhiên cho nhà máy hoặc công trường khác.
Một kết quả có thể kiểm toán nên đi kèm sơ đồ miền khí, áp suất tuyệt đối, nhiệt độ, thể tích thực, vật tính, đường quá trình và định nghĩa công. Đặc biệt, phải ghi rõ đã trừ công đẩy môi trường hay chưa. Nhờ đó người đọc có thể so sánh các phép tính cùng cơ sở, thay vì chọn công thức cho con số nhỏ nhất hoặc xem mọi giá trị có đơn vị kJ là tương đương.
Nguồn tham khảo
[1] Thermodynamics and Propulsion, 1.3 — Changing the State of a System with Heat and Work — Massachusetts Institute of Technology. Truy cập ngày 09/11/2026.
[2] Lecture 26 — Availability Function for a non-flow Process — Indian Institute of Technology Kanpur. Truy cập ngày 09/11/2026.
[3] Isentropic Compression — NASA Glenn Research Center. Truy cập ngày 09/11/2026.
[4] ES&H Manual, Chapter 07 — Pressure Safety, appendices E–F — Lawrence Berkeley National Laboratory, revised 07/2024. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
