Tổng quan nghiên cứu

Thống kê từ Dược điển Hoa Kỳ chỉ ra rằng có khoảng 40% các hoạt chất dược phẩm mới phát hiện thuộc nhóm kém tan hoặc thực tế không tan trong nước, trong khi hơn 80% hoạt chất điều trị hiện nay tồn tại ở dạng rắn. Vấn đề nan giải này trực tiếp làm suy giảm sinh khả dụng và hạn chế hiệu quả điều trị của thuốc trong cơ thể người. Kích thước hạt đóng vai trò quyết định đến tốc độ hòa tan và con đường dẫn truyền dược chất: các hạt thuốc kích thước dưới 5 µm mới có khả năng đi sâu vào phế nang phổi, trong khi kích thước khoảng 0.3 µm là tiêu chuẩn lý tưởng cho hệ dẫn thuốc dạng huyền phù tiêm tĩnh mạch.

Mục tiêu trọng tâm của nghiên cứu là xây dựng mô hình mô phỏng số thủy động lực học nhằm làm sáng tỏ cơ chế hình thành và dự đoán kích thước hạt mầm acetaminophen (paracetamol) thông qua công nghệ giãn nở nhanh dung dịch siêu tới hạn (RESS). Nghiên cứu được triển khai thực hiện từ tháng 06 năm 2013 đến tháng 05 năm 2014 tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh. Luận văn mang ý nghĩa khoa học và thực tiễn sâu sắc khi cung cấp công cụ tính toán mô phỏng CFD chuẩn xác, giúp giảm thiểu hơn 60% chi phí thử nghiệm thực nghiệm tốn kém, đồng thời mở ra hướng phát triển công nghệ sản xuất dược phẩm nano sạch, không sử dụng dung môi hữu cơ độc hại và bảo toàn dược tính sinh học tối ưu.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng trên nền tảng kết hợp giữa các nguyên lý truyền khối, nhiệt động lực học dung dịch siêu tới hạn và động học kết tinh hạt:

  • Phương trình Noyes-Whitney: Thiết lập mối quan hệ định lượng giữa tốc độ hòa tan của dược chất với diện tích bề mặt tiếp xúc và nồng độ bão hòa, khẳng định việc thu nhỏ kích thước hạt xuống cấp độ nano giúp tăng diện tích bề mặt riêng theo cấp số nhân, từ đó cải thiện đáng kể tốc độ hòa tan.
  • Phương trình trạng thái Peng-Robinson: Ứng dụng mô tả chính xác các đặc tính nhiệt động, tỷ trọng và độ chịu nén của lưu chất carbon dioxide (CO2) ở trạng thái siêu tới hạn khi biến thiên qua các miền áp suất và nhiệt độ khác nhau.
  • Mô hình độ hòa tan cân bằng Octavio Elizalde-Solis và Luis A. Galicia-Luna: Xác định mối tương quan giữa độ hòa tan của acetaminophen trong pha CO2 siêu tới hạn thông qua hệ số thực nghiệm chuyên biệt của hệ hai cấu tử.
  • Thuyết tạo mầm đồng thể cổ điển (Classical Nucleation Theory - CNT): Đóng vai trò giải thích cơ chế hình thành các mầm tinh thể tới hạn khi dung dịch đạt trạng thái siêu bão hòa tức thời trong dòng phun siêu âm.

Phương pháp nghiên cứu

Nghiên cứu ứng dụng phương pháp mô phỏng số động lực học dòng chảy (CFD) sử dụng phần mềm ANSYS FLUENT v15 kết hợp đối chứng thực nghiệm:

  • Phương pháp xây dựng lưới tính toán: Mô hình hình học 2D đối xứng trục được thiết kế chính xác theo kích thước thực tế của vòi phun và buồng giãn nở. Cỡ mẫu lưới tính toán bao gồm 101.277 phần tử (elements) dạng tứ giác (Quadrilaterals) với 102.069 nút (nodes), kích thước phần tử cục bộ đạt mức siêu mịn 10^-4 m nhằm thu bắt trọn vẹn hiện tượng sóng xung kích và sự biến thiên gradient áp suất cực lớn.
  • Phương pháp giải số trị: Áp dụng hệ phương trình Navier-Stokes cho dòng lưu chất chịu nén, kết hợp giải đồng thời các phương trình cân bằng liên tục, động lượng và năng lượng trong chế độ dòng chảy dừng.
  • Nguồn dữ liệu và thực nghiệm kiểm chứng: Thí nghiệm tạo hạt được tiến hành trên hệ thống thiết bị RESS của hãng Thar Technologies đặt tại Phòng thí nghiệm trọng điểm Công nghệ Hóa học và Dầu khí, Đại học Bách khoa TP.HCM. Cấu trúc hình thái và kích thước hạt được phân tích thông qua kính hiển vi điện tử quét (SEM) và phân tích phổ hồng ngoại biến đổi Fourier (FT-IR) trong mốc thời gian hoàn thành từ cuối năm 2013 đến giữa năm 2014.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Mô phỏng thủy động lực học và tính toán động học tạo hạt đã cung cấp những kết quả định lượng chi tiết về hành vi của hệ lưu chất trong quá trình RESS:

  • Sự sụt giảm áp suất và nhiệt độ cục bộ: Khi dòng lưu chất CO2 siêu tới hạn đi qua vòi phun và thoát vào buồng giãn nở, áp suất giảm đột ngột từ 120 bar xuống xấp xỉ 1 bar chỉ trong khoảng cách vài milimét, kéo theo sự hạ nhiệt độ mạnh mẽ do hiệu ứng giãn nở đoạn nhiệt Joule-Thomson.
  • Tạo độ quá bão hòa cực đại: Độ quá bão hòa của acetaminophen trong pha khí tăng vọt đạt cấp số cực đại xấp xỉ 10^15 tại vị trí ngay sau khi dòng ra khỏi vòi phun và tiến vào vùng phun tự do siêu âm.
  • Hình thành mầm tinh thể kích thước siêu nhỏ: Tốc độ tạo mầm đồng thể đạt giá trị cực lớn, cho phép tạo ra các hạt mầm tới hạn có bán kính cực mịn chỉ khoảng 0.05 nm với nồng độ hạt mầm phân bố đồng nhất.
  • Thu nhỏ kích thước hạt thực nghiệm: Kết quả phân tích SEM cho thấy sản phẩm hạt acetaminophen thu được sau quá trình giãn nở chuyển dịch từ cấu trúc hình kim thô ban đầu có kích thước 20.8 µm sang dạng hình cầu đồng đều với kích thước trung bình chỉ 0.46 µm, tương đương mức giảm kích thước ấn tượng hơn 97.7%.

Thảo luận kết quả

Độ quá bão hòa khổng lồ 10^15 chính là động lực then chốt thúc đẩy tốc độ tạo mầm đồng đều diễn ra áp đảo so với tốc độ phát triển tinh thể. Hiện tượng này xuất hiện do dung môi CO2 chuyển từ trạng thái lỏng siêu tới hạn có tỷ trọng cao sang trạng thái khí có tỷ trọng thấp trong thời gian phần triệu giây, làm triệt tiêu hoàn toàn khả năng hòa tan của chất tan. Kết quả phân bố số Mach và trường vận tốc dọc đường tâm cho thấy dòng lưu chất đạt tốc độ siêu âm ngay tại miệng vòi phun, tạo nên đĩa Mach và các vùng sóng xung kích tương tự như quan sát của B. Helfgen (2001) và P. Hirunsit (2005).

Dữ liệu mô phỏng được biểu diễn trực quan qua biểu đồ đường cong biến thiên áp suất - nhiệt độ dọc theo trục tọa độ không gian và đồ thị phân bố kích thước hạt SEM. Sự tương thích cao giữa dữ liệu tính toán CFD và ảnh chụp SEM khẳng định tính chính xác của mô hình lý thuyết. Phổ FT-IR cũng xác nhận các liên kết hóa học đặc trưng của acetaminophen không bị biến tính hay phân hủy nhiệt sau quá trình RESS.

Đề xuất và khuyến nghị

Nhằm nâng cao hiệu quả công nghệ và đưa quy trình mô phỏng RESS vào ứng dụng thực tiễn quy mô công nghiệp, 4 giải pháp trọng tâm được khuyến nghị:

  • Tối ưu hóa thiết kế hình học vòi phun: Tiến hành hiệu chỉnh tỷ số chiều dài trên đường kính (L/D) của ống mao dẫn vòi phun ở mức dưới 200 nhằm hạn chế sự kết tinh sớm trong ống và kiểm soát độ phân tán kích thước hạt dưới 5%, thực hiện bởi nhóm kỹ sư cơ khí chế tạo trong vòng 6 tháng.
  • Ứng dụng chất đồng dung môi thân thiện môi trường: Bổ sung từ 0.5% đến 2.0% phần mol đồng dung môi ethanol để gia tăng độ hòa tan của paracetamol trong CO2 siêu tới hạn lên gấp 10 đến 15 lần, do bộ phận nghiên cứu và phát triển (R&D) công nghệ hóa dược triển khai trong thời hạn 9 tháng.
  • Nâng cấp mô hình mô phỏng 3D đa pha: Phát triển mô hình CFD 3D toàn diện tích hợp phương trình cân bằng tổng quát (GDE) với quy mô lưới trên 500.000 phần tử nhằm tính toán chính xác hiện tượng đông tụ do va chạm cơ học trong buồng giãn nở, do các chuyên gia mô phỏng tính toán thực hiện trong 12 tháng.
  • Đánh giá sinh khả dụng và thử nghiệm lâm sàng: Triển khai các thử nghiệm hòa tan in-vitro và hấp thu sinh học in-vivo trên mô hình động vật đối với các lô hạt paracetamol 0.46 µm nhằm chứng minh độ tăng hấp thu trên 85%, thực hiện bởi các viện kiểm nghiệm dược phẩm trong thời gian 18 tháng.

Đối tượng nên tham khảo luận văn

Nội dung và kết quả nghiên cứu của luận văn mang lại giá trị thiết thực cho 4 nhóm đối tượng chuyên ngành:

  • Nhà nghiên cứu công nghệ nano dược phẩm: Khai thác cơ sở dữ liệu về động học tạo mầm và công nghệ RESS để phát triển các dạng thuốc tiêm hoặc thuốc hít phế nang thế hệ mới cho các hoạt chất khó tan.
  • Kỹ sư thiết kế thiết bị quá trình hóa học: Sử dụng các phương trình nhiệt động, thông số thủy động lực học vòi phun và điều kiện biên CFD để tính toán thiết kế chi tiết buồng giãn nở và vòi phun áp suất cao.
  • Doanh nghiệp sản xuất dược phẩm và thực phẩm chức năng: Ứng dụng công nghệ xanh không dung môi độc hại để nâng cao giá trị thương mại, tối ưu hóa độ đồng đều viên nén và cải thiện sinh khả dụng của sản phẩm xuất khẩu.
  • Học viên cao học và sinh viên kỹ thuật hóa học: Sử dụng công trình như một tài liệu tham khảo mẫu mực về phương pháp kết hợp phần mềm mô phỏng ANSYS FLUENT với thực nghiệm kiểm chứng trong các đề tài nghiên cứu lưu chất siêu tới hạn.

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

Tại sao CO2 siêu tới hạn lại là dung môi lý tưởng trong phương pháp RESS? CO2 có nhiệt độ tới hạn tương đối thấp là 31.1°C và áp suất tới hạn 73.8 bar, rất an toàn đối với các hoạt chất sinh học kém bền nhiệt. Khí CO2 hoàn toàn không độc hại, không cháy nổ, chi phí thấp và tự động bay hơi hoàn toàn sau khi giảm áp, giúp thu được hạt thuốc khô tinh khiết mà không để lại bất kỳ dư lượng dung môi nào.

Phương pháp RESS khác biệt gì so với phương pháp nghiền cơ học truyền thống? Phương pháp nghiền truyền thống (top-down) sinh nhiệt ma sát lớn dễ gây biến tính hoạt chất, tạo phân bố kích thước hạt rộng và tiêu hao nhiều năng lượng. Ngược lại, RESS là phương pháp từ dưới lên (bottom-up), tận dụng độ quá bão hòa cực lớn 10^15 để tạo hạt đồng nhất cấp độ nano từ 0.46 µm mà không làm thay đổi cấu trúc hóa học của dược chất.

Mô hình Peng-Robinson có vai trò gì trong việc mô phỏng CFD dòng RESS? Dòng lưu chất trong vòi phun RESS trải qua sự biến thiên áp suất khổng lồ từ trên 120 bar xuống 1 bar trong điều kiện chịu nén mạnh. Phương trình trạng thái Peng-Robinson cung cấp các thông số chính xác về khối lượng riêng, entanpy và hệ số nén của CO2 siêu tới hạn, tạo điều kiện biên chuẩn xác cho hệ phương trình Navier-Stokes trong ANSYS FLUENT.

Làm thế nào để kiểm soát kích thước hạt tạo thành không bị đông tụ lớn hơn? Kích thước hạt có thể được kiểm soát chặt chẽ bằng cách điều chỉnh áp suất hòa tan ở mức 12 đến 16 MPa, gia nhiệt vòi phun lên khoảng 403 K để tránh tắc nghẽn và duy trì khoảng cách phun từ vòi đến đĩa thu hồi trong khoảng 5 đến 7 cm nhằm tối ưu hóa thời gian lưu và giảm thiểu va chạm đông tụ.

Mô hình toán học trong luận văn có thể áp dụng cho các loại dược phẩm khác không? Mô hình hoàn toàn có thể mở rộng áp dụng cho nhiều hoạt chất hữu cơ khác như ibuprofen, aspirin, salicylic acid hoặc cholesterol. Người nghiên cứu chỉ cần thay thế các hệ số thực nghiệm về độ hòa tan cân bằng tương ứng của chất tan trong pha CO2 siêu tới hạn và thiết lập lại các thông số tạo mầm nhiệt động học.

Kết luận

  • Xây dựng thành công mô hình mô phỏng số 2D CFD trên ANSYS FLUENT v15 mô tả chính xác trường thủy động lực học của CO2 siêu tới hạn qua vòi phun và buồng giãn nở.
  • Xác định được độ quá bão hòa cực đại đạt xấp xỉ 10^15 và bán kính hạt mầm tới hạn ở mức 0.05 nm trong vùng phun tự do siêu âm.
  • Kiểm chứng thực nghiệm thành công việc chế tạo hạt acetaminophen hình cầu đồng đều kích thước 0.46 µm từ nguyên liệu thô 20.8 µm trên thiết bị Thar RESS.
  • Khẳng định tính ưu việt của phương pháp RESS trong việc bảo toàn 100% liên kết hóa học đặc trưng của dược chất thông qua kiểm tra phổ FT-IR.
  • Lộ trình 12 đến 18 tháng tới cần tập trung mở rộng mô hình 3D đa pha và thử nghiệm sinh khả dụng trên động vật nhằm thương mại hóa công nghệ sản xuất thuốc nano.

Quý độc giả, các nhà nghiên cứu và doanh nghiệp quan tâm đến công nghệ tạo hạt nano siêu tới hạn hãy kết nối ngay với bộ môn Quá trình và Thiết bị Công nghệ Hóa học để tiếp cận trọn vẹn dữ liệu mô phỏng và thúc đẩy hợp tác chuyển giao công nghệ.