Bối cảnh và vấn đề nghiên cứu

Trong kỹ thuật nhiệt - lạnh và nhiệt động lực học, các thông số trạng thái tại điểm tới hạn bao gồm nhiệt độ tới hạn ($T_c$), áp suất tới hạn ($p_c$) và thể tích phân tử tới hạn ($V_c$) hoặc khối lượng riêng tới hạn ($\rho_c$) đóng vai trò then chốt. Những thông số này là cơ sở trực tiếp để xác định áp suất bão hòa, khối lượng riêng lỏng bão hòa, thiết lập các phương trình trạng thái (như phương trình Van der Waals, phương trình Wagner, phương trình Rackett) và tính toán tính chất nhiệt động của môi chất theo phương pháp đồng dạng.

Về mặt thực nghiệm, việc đo đạc trực tiếp các thông số tới hạn gặp nhiều rào cản kỹ thuật phức tạp: nhiều hợp chất bị phân hủy nhiệt trước khi đạt tới trạng thái tới hạn, môi chất có tính ăn mòn thiết bị hoặc độc hại cao, thông số đo vượt quá thang đo của thiết bị chuyên dụng, và chi phí vận hành thí nghiệm lớn. Do đó, tính đến thời điểm thực hiện đề tài, chỉ có chưa đến 2.000 chất có đầy đủ số liệu đo thực nghiệm về điểm tới hạn, trong khi các thông số trạng thái khác như nhiệt độ sôi, áp suất bão hòa hay khối lượng riêng đã được công bố cho hàng trăm ngàn chất khác nhau. Tại Việt Nam, hướng nghiên cứu xác định tính chất và thông số nhiệt động ở điểm tới hạn còn rất mới, nguồn tài liệu chuyên khảo trong nước rất hạn chế. Đề tài được hình thành trực tiếp từ nhu cầu tìm kiếm, đánh giá các môi chất lạnh thân thiện với môi trường có hiệu suất biến đổi năng lượng cao trong khuôn khổ đề tài nghiên cứu do Quỹ Phát triển khoa học và công nghệ Quốc gia (NAFOSTED) tài trợ cho cán bộ hướng dẫn.

Mục tiêu và nhiệm vụ nghiên cứu của đồ án được xác định cụ thể gồm 5 nội dung:

  1. Tìm hiểu khái niệm nhiệt độ điểm tới hạn, vai trò của nhiệt độ điểm tới hạn trong nghiên cứu và ứng dụng kỹ thuật.
  2. Tìm hiểu, phân loại và hệ thống hóa các công thức, phương pháp xác định nhiệt độ điểm tới hạn đã được công bố trên thế giới.
  3. Tìm kiếm và tổng hợp cơ sở dữ liệu số liệu thực nghiệm chuẩn về điểm tới hạn của các chất.
  4. Phân tích, đánh giá độ chính xác, độ tin cậy của các phương pháp tính toán dự đoán dựa trên cấu tạo phân tử; từ đó nghiên cứu đề xuất phương pháp cải tiến nhằm nâng cao độ chính xác.
  5. Xây dựng phần mềm tính toán tự động xác định nhiệt độ tới hạn và các thông số nhiệt động liên quan.

Đối tượng nghiên cứu của đồ án là nhiệt độ tới hạn ($T_c$) của các hợp chất thuần và môi chất nhiệt - lạnh, cùng các phương pháp dự đoán theo cấu tạo phân tử (phương pháp đóng góp nhóm). Phạm vi nghiên cứu tập trung vào các hợp chất hữu cơ, đối chiếu và kiểm nghiệm trên các tập dữ liệu thực nghiệm chuẩn quốc tế (tập 24 chất khảo sát ban đầu, tập 173 chất có số liệu thực nghiệm chuẩn, cùng các tập dữ liệu từ 398 đến 572 chất theo tài liệu đối sánh). Thời gian thực hiện và bảo vệ đồ án tốt nghiệp vào năm 2015 tại Viện Khoa học và Công nghệ Nhiệt - Lạnh, Trường Đại học Bách khoa Hà Nội.

Cơ sở lý thuyết và phương pháp nghiên cứu

Về mặt nhiệt động lực học, điểm tới hạn là trạng thái mà tại đó ranh giới giữa hai pha lỏng và hơi hoàn toàn biến mất, khối lượng riêng của pha lỏng và pha hơi bằng nhau. Dưới nhiệt độ tới hạn ($T_c$), chất có thể tồn tại ở trạng thái cân bằng lỏng - hơi; trên nhiệt độ này, chất chỉ tồn tại ở pha khí và không thể hóa lỏng bằng việc tăng áp suất đơn thuần. Lý thuyết động học phân tử giải thích hiện tượng này thông qua sự cạnh tranh giữa lực hấp dẫn tiềm tàng (liên kết phân tử tạo pha lỏng) và năng lượng chuyển dịch động học (lực phân tán ngẫu nhiên tạo pha khí). Tại điểm tới hạn, năng lượng chuyển dịch động học cân bằng với năng lượng hấp dẫn tiềm tàng lớn nhất.

Các thông số điểm tới hạn được tác giả hệ thống hóa qua các mô hình nhiệt động kinh điển:

  • Hệ số nén tới hạn ($Z_c$): $$Z_c = \frac{p_c V_c}{R T_c}$$ Trong đó $p_c$ là áp suất tới hạn (bar), $V_c$ là thể tích tới hạn ($\text{cm}^3/\text{mol}$), $T_c$ là nhiệt độ tới hạn (K), $R = 8{,}31451 \text{ J}/(\text{mol}\cdot\text{K})$ là hằng số chất khí.
  • Phương trình trạng thái Van der Waals: $$p = \frac{RT}{v - b} - \frac{a}{v^2}$$ Với hai hằng số thực nghiệm được hiệu chỉnh trực tiếp từ thông số tới hạn: $$a = \frac{27 R^2 T_c^2}{64 p_c}, \quad b = \frac{R T_c}{8 p_c}$$
  • Phương trình áp suất bão hòa Wagner (1973, 1977): $$\ln P_r = \frac{A(1 - T_r) + B(1 - T_r)^{1{,}5} + C(1 - T_r)^3 + D(1 - T_r)^6}{T_r}$$ Với $T_r = T/T_c$, $P_r = P/p_c$, và $A, B, C, D$ là các hằng số Wagner.
  • Phương trình Rackett (1970) xác định thể tích mol lỏng bão hòa ($V_s$): $$V_s = \frac{R T_c}{p_c} Z_c^{\left[1 + (1 - T/T_c)^{2/7}\right]}$$

Về mặt phương pháp luận, đồ án tiếp cận theo hướng dự đoán tính chất nhiệt động dựa trên cấu tạo phân tử (phương pháp đóng góp nhóm - Group Contribution Method). Xuất phát từ quan sát của Guldberg (1890) rằng $T_c \approx 1{,}5 T_b$ (với $T_b$ là nhiệt độ sôi chuẩn ở áp suất khí quyển), Riedel (1949) và Lydersen (1955) đã chuẩn hóa công thức tính nhiệt độ tới hạn dưới dạng tỉ số phụ thuộc vào tổng đóng góp cấu trúc: $$T_c = \frac{T_b}{\theta}$$ Trong đó $\theta$ được tính từ tổng các số gia đóng góp của các nguyên tử hoặc nhóm nguyên tử cấu thành phân tử.

Tên phương pháp Năm công bố Số nhóm cấu trúc Công thức xác định $T_c$ Cỡ mẫu kiểm nghiệm (chất) Sai số tuyệt đối trung bình
Guldberg 1890 0 $T_c = 1{,}5 T_b$ — —
Riedel 1949, 1952 22 $T_c = \frac{T_b}{0{,}574 + \sum n_i \Delta T_i}$ — —
Lydersen 1955 32 $T_c = \frac{T_b}{0{,}567 + \sum n_i \Delta T_i - \left(\sum n_i \Delta T_i\right)^2}$ 557 $10{,}7 \text{ K } (1{,}71%)$
Ambrose 1978, 1979 30 $T_c = T_b \left[1 + \left(1{,}242 + \sum n_i \Delta T_i\right)^{-1}\right]$ (có hiệu chỉnh số Delta Platt) 528 $6{,}07%$
Daubert 1980 106 $T_c = 1{,}806 T_b - \sum n_i \Delta T_i$ 475 $23{,}9 \text{ K } (3{,}8%)$
Klincewicz & Reid 1984 35 $T_c = 45{,}40 + 0{,}77 M + 1{,}55 T_b - \sum n_i \Delta T_i$ 398 —
Joback 1984, 1987 41 $T_c = T_b \left[0{,}584 + 0{,}965 \sum \Delta T_c - \left(\sum \Delta T_c\right)^2\right]^{-1}$ 173 Đánh giá tại Chương 5
Somayajulu 1989 169 $T_c = T_b \left[1 - \left(1{,}242 + b_t \sum n_t \Delta T_i\right)^{-1}\right]$ (với $b_t = 0{,}138$) 517 $8{,}39 \text{ K } (1{,}46%)$
Constantinou & Gani 1994 Bậc 1 & 2 $\exp(T_c / 181{,}128) = \sum N_i C_i + W \sum M_j D_j$ 559 $17{,}2 \text{ K } (4{,}06%)$
Tu 1995 40 Dựa trên lý thuyết mạng lỗ trống Kurata-Isida 572 $23{,}3 \text{ K } (4{,}33%)$
Marrero & Pardillo 1999 39 $T_c = T_b \left[0{,}5851 + 0{,}9286 \sum N_i C_i - \left(\sum N_i C_i\right)^2\right]^{-1}$ (mô hình GIC) 458 $7{,}8 \text{ K } (1{,}21%)$
Wen & Quiang 2001 Nhóm liền kề + 13 cơ bản Dùng công thức hàm mũ khi biết hoặc chưa biết $T_b$ — —

Phương pháp thực hiện của tác giả bao gồm: thu thập dữ liệu nhiệt động từ các công bố quốc tế chuẩn; phân tích toán học các công thức đóng góp nhóm; thiết lập thuật toán hiệu chỉnh mô hình Joback dựa trên cấu trúc phân tử tương đồng của chất tham chiếu; kiểm chứng sai số độc lập; và tin học hóa thuật toán bằng ngôn ngữ lập trình. Nguồn dữ liệu thực nghiệm dùng trong đề tài gồm bảng giá trị chuẩn của 173 hợp chất có đầy đủ thông số điểm tới hạn, cơ sở dữ liệu nhiệt độ sôi của 18.000 chất, cùng số liệu tổng hợp kiểm tra từ nghiên cứu của Yash Nannoolal (2007).

Thiết kế và triển khai phần mềm tính toán

Phần mềm trong đồ án được thiết kế nhằm tin học hóa quy trình xác định nhiệt độ tới hạn và các thông số nhiệt động của các chất hữu cơ theo phương pháp Joback và phương pháp Joback cải tiến.

Cấu trúc chức năng của phần mềm được xây dựng gồm các khối xử lý:

  • Khối tính toán nhiệt độ tới hạn theo phương pháp Joback: Sử dụng công thức toán học thực nghiệm của Joback dựa trên tổng các số gia $\Delta T_c$ của 41 nhóm chức cơ bản. Trường hợp chưa có giá trị thực nghiệm của nhiệt độ sôi $T_b$, phần mềm tích hợp module tính $T_{b_\text{est}} = 198 + \sum \Delta T_b$.
  • Khối tính toán nhiệt độ tới hạn theo phương pháp Joback cải tiến: Cho phép người dùng nhập thông tin cấu trúc phân tử của chất cần tính và lựa chọn chất tham chiếu (Ref) có cấu trúc phân tử tương đồng để hiệu chỉnh sai số của mô hình gốc.
  • Khối tính toán 11 thông số nhiệt động và nhiệt hóa: Sử dụng bảng số liệu đóng góp nhóm mở rộng của Joback để tính toán tự động 11 thông số, bao gồm nhiệt dung riêng khí lý tưởng ($C_p$ dưới dạng hàm nhiệt độ), biến thiên enthalpy tạo thành, năng lượng tự do Gibbs, nhiệt hóa hơi, độ nhớt động học và các thông số trạng thái khác.
  • Khối tương tác và xuất nhập dữ liệu: Giao diện trực quan cho phép người dùng lựa chọn số lượng từng nhóm chức năng (mạch thẳng, mạch nhánh, mạch vòng, nhóm chứa oxy, nitơ, lưu huỳnh, halogen); hiển thị đồ thị đường đẳng nhiệt thực nghiệm; và tích hợp chức năng trích xuất toàn bộ bảng số liệu cùng kết quả tính toán ra định dạng bảng tính Microsoft Excel.
  • Triển khai và phân phối: Phần mềm được đóng gói để phục vụ miễn phí cho công tác đào tạo, nghiên cứu khoa học chuyên ngành Kỹ thuật Nhiệt - Lạnh và có khả năng truy cập qua môi trường mạng internet.

Nội dung chính theo từng chương

Chương 1: Tổng quan về nhiệt độ tới hạn

Chương này khái quát định nghĩa vật lý về điểm tới hạn, lịch sử phát hiện từ thí nghiệm nòng đại bác của Cagniard de La Tour (1822) đến các nghiên cứu hóa lỏng khí carbon dioxide ($CO_2$) ở 31,1°C của Thomas Andrews (1869), và các nghiên cứu thực nghiệm của Kobe và Lynn (1953). Tác giả trình bày vai trò cốt lõi của $T_c$ trong việc tính toán hệ số nén tới hạn $Z_c$, giải các phương trình trạng thái thực (Van der Waals, Wagner, Rackett) và ứng dụng trong nghiên cứu môi chất lạnh mới. Chương 1 phân loại các phương pháp xác định $T_c$ thành hai hướng chính: phương pháp thực nghiệm và phương pháp dự đoán (trong đó nhóm phương pháp dự đoán theo cấu tạo phân tử được lựa chọn làm trọng tâm nghiên cứu của đồ án).

Chương 2: Phương pháp đóng góp nhóm

Tác giả hệ thống hóa chi tiết 11 phương pháp đóng góp nhóm tiêu biểu theo trình tự thời gian: Riedel (1949), Lydersen (1955), Ambrose (1978, 1979), Daubert (1980), Klincewicz & Reid (1984), Joback (1984, 1987), Somayajulu (1989), Constantinou & Gani (1994), Tu (1995), Marrero-Morejon & Pardillo-Fontdevilla (1999), Marrero-Morejon & Gani (2001), và Wen & Quiang (2001). Chương 2 cung cấp toàn bộ các bảng giá trị số gia đóng góp cấu trúc của từng phương pháp. Tác giả tiến hành tính toán kiểm tra, so sánh độ chính xác của các phương pháp trên tập dữ liệu thực nghiệm gồm 24 chất điển hình, phân tích sai số tuyệt đối trung bình và chỉ ra ưu điểm đơn giản, trực quan của phương pháp Joback cũng như hạn chế sai số cần khắc phục.

Chương 3: Phương pháp Joback cải tiến

Chương 3 đi sâu phân tích ưu và nhược điểm của phương pháp Joback truyền thống. Phương pháp Joback có ưu điểm là số lượng nhóm cấu trúc vừa phải (41 nhóm), dễ áp dụng cho hầu hết các hợp chất hữu cơ, thuật toán tính toán đơn giản. Tuy nhiên, nhược điểm là phương pháp bỏ qua ảnh hưởng vị trí tương hỗ giữa các nhóm kề cận và sự phân nhánh phức tạp, dẫn đến sai số lớn ở một số nhóm chất đặc thù. Để giải quyết hạn chế này, tác giả đề xuất phương pháp Joback cải tiến: sử dụng nguyên lý cấu tạo phân tử kết hợp với chất tham chiếu (Ref) có cấu trúc gần tương đồng đã có sẵn số liệu thực nghiệm chuẩn. Chương này xây dựng chi tiết lưu đồ thuật toán, quy trình xác định và lập bảng kiểm tra độ chính xác của phương pháp Joback cải tiến.

Chương 4: Lập trình phần mềm

Chương này mô tả toàn bộ quá trình xây dựng phần mềm tính toán. Tác giả trình bày công cụ lập trình, thiết kế lưu đồ thuật toán chi tiết cho cả hai trường hợp (Joback gốc và Joback cải tiến), thiết kế giao diện đồ họa người dùng (GUI) trực quan. Chương 4 diễn giải cấu trúc bảng đóng góp nhóm cho 11 thông số nhiệt động được tích hợp trong mã nguồn, cung cấp ví dụ mẫu tính toán cụ thể cho từng chất, đồng thời trình bày các module chức năng: xác định $T_c$, xác định nhiệt dung riêng khí lý tưởng, xuất dữ liệu và kết quả ra Excel, và hiển thị đồ thị đường đẳng nhiệt thực nghiệm.

Chương 5: Kết quả và đánh giá

Chương này tập trung kiểm chứng và đánh giá sai số giữa kết quả dự đoán từ phần mềm với các giá trị đo thực nghiệm chuẩn. Nội dung đánh giá được chia làm hai trường hợp: trường hợp đã biết nhiệt độ sôi thực nghiệm ($T_b$) và trường hợp chưa biết nhiệt độ sôi (phải ước tính $T_b$ qua số gia nhóm). Tác giả thực hiện kiểm nghiệm chi tiết trên tập cơ sở dữ liệu gồm 173 chất có nguồn $T_c$ chuẩn. Đồng thời, chương 5 lập bảng phân tích chuyên sâu cho các tập hợp chất có độ lệch cao (gồm nhóm 64 chất và nhóm 70 chất có sai số lớn), đánh giá tương quan giữa hai đại lượng $T_{c_\text{exp_est}}$ và $T_{c_\text{est_exp}}$, qua đó khẳng định mức độ tin cậy và sự cải thiện độ chính xác của phương pháp cải tiến.

Chương 6: Kết luận và kiến nghị

Chương cuối tóm lược các kết quả đạt được: hoàn thành nghiên cứu lý thuyết về các phương pháp đóng góp nhóm; đề xuất phương pháp Joback cải tiến đạt độ chính xác cao hơn phương pháp gốc; xây dựng hoàn chỉnh công cụ phần mềm tính toán giao diện trực quan, miễn phí và hỗ trợ trích xuất Excel; đề xuất kiến nghị mở rộng cơ sở dữ liệu kiểm nghiệm và tích hợp thêm các phương pháp tính toán bậc cao trong tương lai.

Kết quả và đóng góp

Về mặt kết quả số liệu, đồ án đã thực hiện tính toán kiểm chứng trên các tập dữ liệu thực nghiệm độc lập:

  • Kiểm tra, so sánh độ chính xác của các phương pháp dự đoán $T_c$ trên tập 24 chất đại diện cho các nhóm hydrocarbon, dẫn xuất halogen, rượu, acid và este.
  • Đánh giá độ chính xác tổng thể của phương pháp Joback gốc và Joback cải tiến trên tập 173 chất hữu cơ có nguồn số liệu chuẩn về điểm tới hạn ($T_c^\text{ref}$).
  • Phân tích chi tiết độ lệch và nguyên nhân sai số trên hai tập hợp chất đặc thù: tập 64 chất và tập 70 chất có sai số tính toán lớn, chỉ ra rằng sai số chủ yếu xuất hiện ở các chất có cấu trúc mạch nhánh phức tạp hoặc chứa nhiều nhóm chức phân cực cao khi chưa biết nhiệt độ sôi thực nghiệm.
  • Kết quả đối sánh cho thấy việc sử dụng nhiệt độ sôi thực nghiệm ($T_b$) mang lại độ chính xác cao hơn rõ rệt so với trường hợp sử dụng giá trị nhiệt độ sôi ước tính ($T_{b_\text{est}} = 198 + \sum \Delta T_b$). Phương pháp Joback cải tiến giúp thu hẹp đáng kể độ lệch đối với các hợp chất có sẵn chất tham chiếu (Ref) tương đồng.

Về giải pháp và đóng góp khoa học, tác giả đã:

  • Đề xuất quy trình và thuật toán tính toán theo phương pháp Joback cải tiến, khai thác thông tin từ chất tham chiếu nhằm bù trừ sai số hệ thống của phương pháp đóng góp nhóm bậc một truyền thống.
  • Xây dựng hoàn chỉnh phần mềm tính toán tự động bằng giao diện đồ họa, tích hợp cơ sở dữ liệu đóng góp nhóm tính toán đồng thời 11 thông số nhiệt động - nhiệt hóa, cho phép kết nối xuất số liệu sang Microsoft Excel.
  • Đóng góp công cụ tính toán phục vụ trực tiếp cho đề tài NAFOSTED nghiên cứu tìm kiếm môi chất lạnh thân thiện môi trường tại Viện Khoa học và Công nghệ Nhiệt - Lạnh.

Hạn chế và hướng nghiên cứu tiếp

Hạn chế của nghiên cứu được xác định từ nội dung văn bản và phạm vi phân tích:

  • Phương pháp Joback và mô hình cải tiến vẫn dựa trên giả thiết cộng tính độc lập của các nhóm cấu trúc, chưa bao quát đầy đủ hiệu ứng tương tác không gian và hiệu ứng lập thể của các phân tử dị vòng hoặc hợp chất hữu cơ đa chức năng có cấu trúc phân tử lớn.
  • Khi chất nghiên cứu chưa có số liệu thực nghiệm về nhiệt độ sôi ($T_b$), sai số ước tính của $T_b$ sẽ lan truyền trực tiếp vào công thức tính $T_c$, làm tăng độ lệch của kết quả đầu ra.
  • Cơ sở dữ liệu kiểm chứng độc lập của đồ án mới dừng lại ở tập 173 chất có số liệu thực nghiệm chuẩn, chưa kiểm tra toàn bộ gần 2.000 chất đã có số liệu trên thế giới do hạn chế về nguồn tài liệu thu thập tại thời điểm thực hiện.

Hướng nghiên cứu tiếp theo do tác giả đề xuất:

  • Tiếp tục bổ sung và mở rộng tập dữ liệu thực nghiệm chuẩn để kiểm chứng phương pháp Joback cải tiến trên phạm vi môi chất rộng hơn.
  • Nghiên cứu tích hợp các phương pháp đóng góp nhóm bậc hai, bậc ba (như Constantinou & Gani, Marrero & Gani) vào phần mềm nhằm tăng độ chính xác cho các hợp chất dị vòng và phân tử chuỗi dài.
  • Tối ưu hóa thuật toán tự động nhận diện và đề xuất chất tham chiếu (Ref) tối ưu cho từng họ chất trong phương pháp Joback cải tiến.

Giá trị tham khảo

Đồ án là tài liệu tham khảo kỹ thuật hữu ích cho:

  • Sinh viên, học viên cao học và cán bộ nghiên cứu chuyên ngành Kỹ thuật Nhiệt - Lạnh, Máy và Thiết bị Nhiệt, Công nghệ Kỹ thuật Hóa học và Hóa dầu trong việc tra cứu tính chất nhiệt vật lý của môi chất.
  • Các nhóm nghiên cứu mô phỏng chu trình nhiệt động, tính toán thiết bị trao đổi nhiệt khi cần xác định thông số điểm tới hạn và phương trình trạng thái của các chất chưa có số liệu thực nghiệm.
  • Phần đáng tham khảo nhất là hệ thống bảng tổng hợp số gia đóng góp nhóm của 11 phương pháp quốc tế (Chương 2), thuật toán hiệu chỉnh chất tham chiếu của phương pháp Joback cải tiến (Chương 3), và bảng số liệu đóng góp nhóm của 11 thông số nhiệt động tích hợp trong phần mềm (Chương 4).

Câu hỏi thường gặp

Tại sao tác giả lựa chọn nghiên cứu phương pháp dự đoán theo cấu tạo phân tử thay vì tiến hành đo đạc thực nghiệm?

Đo đạc thực nghiệm điểm tới hạn đòi hỏi kỹ thuật phức tạp, thời gian dài, chi phí tốn kém, và nhiều chất dễ bị phân hủy nhiệt hoặc ăn mòn thiết bị đo ở vùng nhiệt độ cao. Do đó, trên thế giới chỉ có chưa đến 2.000 chất có số liệu đo thực nghiệm điểm tới hạn, trong khi các thông số khác đã được công bố cho hàng trăm ngàn chất. Tại Việt Nam, điều kiện thiết bị thực nghiệm còn hạn chế nên phương pháp dự đoán dựa trên các nhóm nguyên tử cấu thành phân tử là giải pháp tối ưu.

Phương pháp Joback gốc xác định nhiệt độ tới hạn và nhiệt độ sôi như thế nào?

Theo Joback (1984, 1987), nhiệt độ tới hạn được tính theo công thức: $$T_c = T_b \left[0{,}584 + 0{,}965 \sum \Delta T_c - \left(\sum \Delta T_c\right)^2\right]^{-1}$$ Trong đó $\sum \Delta T_c = \sum n_i \Delta T_{ci}$ là tổng đóng góp của các nhóm cấu trúc. Trong trường hợp chưa có nhiệt độ sôi thực nghiệm $T_b$, giá trị $T_b$ được ước lượng qua công thức: $$T_{b_\text{est}} = 198 + \sum \Delta T_b$$ với $\Delta T_b = \sum n_i \Delta T_{bi}$ lấy từ bảng số gia đóng góp của Joback.

Nguyên lý cơ bản của phương pháp Joback cải tiến được đề xuất trong đồ án là gì?

Phương pháp Joback cải tiến khắc phục nhược điểm bỏ qua tương tác cấu trúc của phương pháp Joback gốc bằng cách kết hợp thông tin cấu tạo phân tử của chất cần tính với một chất tham chiếu (Ref). Chất tham chiếu là chất có cấu trúc phân tử tương đồng và đã có sẵn số liệu thực nghiệm chuẩn về nhiệt độ tới hạn, từ đó dùng độ lệch của chất tham chiếu để hiệu chỉnh giá trị tính toán cho chất mục tiêu.

Phần mềm xây dựng trong đồ án có những tính năng kỹ thuật nổi bật nào?

Phần mềm cho phép: tính toán tự động $T_c$ theo phương pháp Joback gốc và Joback cải tiến; tính toán 11 thông số nhiệt động - nhiệt hóa dựa trên cấu tạo phân tử (bao gồm nhiệt dung riêng khí lý tưởng $C_p$); hiển thị đồ thị đường đẳng nhiệt thực nghiệm; và trích xuất trực tiếp toàn bộ dữ liệu cấu trúc cùng kết quả tính toán ra bảng tính Excel.

Tập dữ liệu nào đã được sử dụng để kiểm nghiệm và đánh giá độ chính xác của các mô hình trong đồ án?

Tác giả đã sử dụng tập dữ liệu 24 chất ban đầu để so sánh các phương pháp đóng góp nhóm; sau đó kiểm chứng toàn diện mô hình trên tập cơ sở dữ liệu 173 chất có số liệu $T_c$ chuẩn thực nghiệm; đồng thời lập danh mục phân tích chuyên sâu sai số trên hai nhóm hợp chất đặc thù gồm 64 chất và 70 chất có sai số lớn.

Kết luận

Đồ án tốt nghiệp của sinh viên Nguyễn Trường Phái đã hệ thống hóa toàn diện cơ sở lý thuyết và các mô hình đóng góp nhóm xác định nhiệt độ tới hạn, chỉ rõ vai trò của thông số tới hạn trong các phương trình trạng thái Van der Waals, Wagner và Rackett. Trên cơ sở phân tích ưu nhược điểm của phương pháp Joback, tác giả đã đề xuất phương pháp Joback cải tiến dựa trên chất tham chiếu tương đồng, giúp cải thiện độ chính xác dự đoán trên tập dữ liệu kiểm nghiệm 173 chất chuẩn. Công trình đã hiện thực hóa thành công thuật toán thành phần mềm tính toán chuyên dụng, hỗ trợ xác định đồng thời 11 thông số nhiệt động và trích xuất dữ liệu ra Excel, đóng góp thiết thực cho công tác đào tạo và nghiên cứu môi chất lạnh mới trong khuôn khổ đề tài NAFOSTED tại Đại học Bách khoa Hà Nội.