Đang chuẩn bị nút TẢI XUỐNG, xin hãy chờ
Tải xuống
Trong công trình này chúng tôi đưa ra các kết quả mô phỏng cân bằng lỏng hơi của flo lỏng bằng kỹ thuật mô phỏng Monte Carlo toàn cục sử dụng thế cặp tương tác phân tử 5 vị trí mới ab initio đã được chúng tôi đưa ra từ các tính toán cơ học lượng tử đối với dime F2-F2. Các kết quả mô phỏng bao gồm tỷ trọng, áp suất hơi và entanpi hóa hơi đã được so sánh với các số liệu thực nghiệm và với những số liệu lấy từ các tài liệu. | Mô phỏng Monte Carlo cân bằng lỏng hơi của Flo sử dụng các thế tương tác phân tử AB Initio mới Tạp chí KHOA HỌC ĐHSP TP. HCM Phạm Văn Tất MÔ PHỎNG MONTE CARLO CÂN BẰNG LỎNG HƠI CỦA FLO SỬ DỤNG CÁC THẾ TƯƠNG TÁC PHÂN TỬ AB INITIO MỚI Phạm Văn Tất* 1. Giới thiệu Mô phỏng máy tính trở thành công cụ bắt buộc cần thiết cho nghiên cứu chất lỏng và hỗn hợp chất lỏng [1]. Một trong các cố gắng đầu tiên của Nasrabad và Deiters dự đoán cân bằng pha lỏng hơi áp suất cao của hỗn hợp khí hiếm [2,3] từ các mô phỏng toàn cục sử dụng thế tương tác phân tử. Leonhard và Deiters đã sử dụng thế Morse 5 vị trí để mô tả thế cặp của nitơ [4] và đã dự đoán áp suất hơi và tỷ trọng cân bằng pha bằng mô phỏng Monte Carlo toàn cục[5]. Trong công trình này chúng tôi đưa ra các kết quả mô phỏng cân bằng lỏng hơi của flo lỏng bằng kỹ thuật mô phỏng Monte Carlo toàn cục sử dụng thế cặp tương tác phân tử 5 vị trí mới ab initio đã được chúng tôi đưa ra từ các tính toán cơ học lượng tử đối với dime F2-F2. Các kết quả mô phỏng bao gồm tỷ trọng, áp suất hơi và entanpi hóa hơi đã được so sánh với các số liệu thực nghiệm và với những số liệu lấy từ các tài liệu. 2. Phương pháp tính toán 2.1. Phương pháp mô phỏng Mô phỏng GEMC-NPT được sử dụng để kiểm tra sự chính xác của các thế cặp. Mô phỏng này được khảo sát ở điều kiện đẳng áp 1,0 MPa và 10,0 MPa trong khoảng nhiệt độ từ 90,0 K đến 270,0 K. Mô phỏng GEMC-NVT được tiến hành để nhận được tỷ trọng hơi và lỏng, áp suất hơi trong khoảng 60,0 K đến 140,0 K với thay đổi 10,0 K. Các thế cặp 5 vị trí mới của flo dime được dùng cho cả hai trường hợp mô phỏng: 5 5 ij ( rij ij ) 2 C ij qi q j u1 Deij ((1 e ) 1) f1a (rij ) nn f 2 ( rij ) (1) i 1 i 1 n 6,8,10 r 4 0 rij ij * TS. – Trường ĐH Đà Lạt 63 Tạp chí KHOA HỌC ĐHSP TP. HCM Số 16 năm 2009 5 5 ( rij ij ) 2 Cnij qi q j u 2 Deij ((1 e ij ) 1) f1b (rij ) f