I. Giới thiệu chung về Tối Ưu Dạng Trong Cơ Học Chất Lỏng
Bài toán tối ưu dạng trong cơ học chất lỏng đã thu hút sự quan tâm lớn từ kỹ sư và nhà toán học từ những năm 1970. Ứng dụng thực tế của nó rất đa dạng, bao gồm thiết kế cánh máy bay tối ưu để giảm lực cản, thiết kế động mạch chủ trong phẫu thuật mạch vành nhân tạo, và nghiên cứu hình dạng tối ưu của vi sinh vật trong chất lỏng nhớt. Các bài toán tối ưu topology nhằm cực tiểu năng lượng tiêu tán của dòng chảy Stokes cũng được nghiên cứu. Bài toán tương tự cũng được xem xét cho dòng chảy Darcy-Stokes, Navier-Stokes với số Reynolds thấp đến trung bình. Mục đích chính của nghiên cứu này là xây dựng thuật toán giải số cho bài toán tối ưu dạng trong cơ học chất lỏng. Cho Ω là một miền mở, bị chặn trong không gian Rd ( d = 2 hoặc 3 ). Một dòng chảy chất lỏng đặt trên Ω được mô tả bởi hệ phương trình Stokes: − ν∆u + ∇ p = f trong Ω. Bài toán này đã được xem xét trong nhiều nghiên cứu trước đây [9], [15], [20]. Một cách hình thức, mô hình bài toán tối ưu dạng được phát biểu như sau: inf J ( Ω ). Các phương pháp giải bài toán tối ưu dạng trong cơ học tính toán đã được giới thiệu từ những năm 1960, trong đó, các kỹ thuật được phân biệt chủ yếu dựa trên cách biểu diễn dạng và các phân tích độ nhạy cảm (sensitivity analysis) của dạng.
1.1. Tầm quan trọng của tối ưu hình học chất lỏng trong kỹ thuật
Bài toán thiết kế tối ưu hóa chất lỏng đóng vai trò quan trọng trong nhiều lĩnh vực kỹ thuật, từ hàng không vũ trụ đến y sinh. Việc tối ưu hình dạng dòng chảy có thể dẫn đến cải thiện đáng kể hiệu suất, giảm tiêu thụ năng lượng và tăng tuổi thọ của thiết bị. Ví dụ, trong thiết kế cánh máy bay, việc tối ưu hình dạng có thể giúp giảm đáng kể lực cản và tăng hiệu suất khí động học. Trong lĩnh vực y sinh, thiết kế kênh dẫn chất lỏng tối ưu trong các thiết bị y tế có thể cải thiện hiệu quả điều trị và giảm thiểu tác dụng phụ. Pironneau đã có những nghiên cứu tiên phong trong lĩnh vực này từ thập niên 70, mở đường cho nhiều ứng dụng thực tế.
1.2. Tổng quan các phương pháp tối ưu ngược chất lỏng hiện nay
Các phương pháp tối ưu ngược chất lỏng hiện nay bao gồm phương pháp hiện, phương pháp đồng nhất và phương pháp ẩn. Phương pháp hiện biểu diễn dạng bởi tập các bậc tự do. Phương pháp đồng nhất dựa trên hàm mật độ. Phương pháp ẩn, với đại diện là phương pháp tập mức, xác định dạng dựa trên đường mức 0 của hàm mức đặc trưng. Các kỹ thuật biểu diễn dạng cần thoả mãn hai điều kiện: cho phép thực hiện các tính toán trên dạng được xét và đảm bảo đủ linh hoạt cho các phép biến dạng phức tạp của miền. Các nghiên cứu của [17], [22]–[25] đã trình bày chi tiết về các phương pháp đồng nhất.
II. Vấn đề và Thách thức trong Tối Ưu Hóa Hình Dạng Chất Lỏng
Bài toán tối ưu hóa hình dạng chất lỏng đặt ra nhiều thách thức về mặt toán học và tính toán. Thứ nhất, việc biểu diễn hình dạng và thực hiện các biến đổi hình học phức tạp đòi hỏi các phương pháp số hiệu quả và chính xác. Thứ hai, việc giải các phương trình Navier-Stokes hoặc Stokes để mô phỏng dòng chảy chất lỏng thường tốn kém về mặt tính toán, đặc biệt khi xét đến các dòng chảy rối và ba chiều. Thứ ba, việc tìm kiếm hình dạng tối ưu trong không gian tham số lớn có thể dẫn đến các bài toán tối ưu không lồi và đa cực trị, đòi hỏi các thuật toán tối ưu mạnh mẽ và khả năng tránh các cực trị cục bộ. Giới hạn trong nội dung luận văn, tác giả tìm hình dạng tối ưu nhằm cực tiểu hàm mục tiêu năng lượng tiêu tán của dòng chảy Stokes dưới ràng buộc về thể tích dạng. Bài toán được giải quyết nhờ phương pháp Lagrange tăng cường thông qua đạo hàm dạng bậc nhất của các hàm dạng tương ứng.
2.1. Khó khăn trong mô phỏng dòng chảy rối và ba chiều mô phỏng CFD
Mô phỏng dòng chảy rối và ba chiều bằng các phương pháp CFD tối ưu hóa hình dạng đòi hỏi lưới tính toán rất mịn và các mô hình rối phù hợp. Việc này làm tăng đáng kể chi phí tính toán và thời gian mô phỏng. Các phương pháp phần tử hữu hạn (FEM) và thể tích hữu hạn (FVM) cần được lựa chọn và triển khai cẩn thận để đảm bảo tính chính xác và ổn định của kết quả. Các nghiên cứu của [32], [33] đã thảo luận về các tính chất phức tạp của dòng chảy rối.
2.2. Vấn đề cực trị cục bộ trong bài toán tối ưu đa mục tiêu
Bài toán tối ưu đa mục tiêu trong cơ học chất lỏng thường có nhiều cực trị cục bộ, khiến việc tìm kiếm nghiệm tối ưu toàn cục trở nên khó khăn. Các thuật toán tối ưu như thuật toán di truyền hoặc thuật toán đàn kiến có thể được sử dụng để khám phá không gian tham số rộng hơn và tăng khả năng tìm thấy nghiệm tốt hơn. Tuy nhiên, việc lựa chọn tham số phù hợp cho các thuật toán này cũng là một thách thức. Cần phải có chiến lược thích hợp để giải quyết vấn đề này khi thực hiện tối ưu hóa hình học chất lỏng.
III. Phương Pháp Tối Ưu Dạng Dựa Trên Đạo Hàm Dạng
Nội dung luận văn được trình bày trong hai chương: i. Chương 1: trình bày một lược đồ mô phỏng số cho phương trình dòng chảy Navier- Stokes, dựa trên phương pháp đặc trưng và phương pháp phần tử hữu hạn. Phần này sẽ trình bày chi tiết các bước xây dựng lược đồ giải số cho một bài toán đạo hàm riêng cụ thể sử dụng phương pháp phần tử hữu hạn, là cơ sở phục vụ cho việc xây dựng các lược đồ giải số cho dòng chảy Stokes và bài toán đạo hàm riêng elliptic sẽ được trình bày trong chương 2. Chương 2: trình bày thuật toán giải số bài toán tối ưu dạng cho dòng chảy Stokes, dựa trên phương pháp biến phân Hadamard và phương pháp phần tử hữu hạn. Mã nguồn chương trình của các thuật toán và các thí nghiệm giải số được phát triển trên phần mềm Freefem++. Đây là một công cụ mã nguồn mở giải các bài toán đạo hàm riêng bằng phương pháp phần tử hữu hạn [31]. Các kết quả hình ảnh sử dụng trong luận văn được hỗ trợ bởi các phần mềm mã nguồn mở gnuplot (http://www. info/) và Medit (http://www.fr/frey/software.
3.1. Tính độ nhạy hình dạng và phân tích độ nhạy trong tối ưu topo chất lỏng
Việc tính toán độ nhạy hình dạng và thực hiện phân tích độ nhạy là bước quan trọng để xác định hướng biến đổi hình dạng nhằm cải thiện hàm mục tiêu. Các phương pháp như phương pháp biến phân Hadamard có thể được sử dụng để tính toán đạo hàm dạng của hàm mục tiêu. Các kết quả này sau đó được sử dụng để xây dựng thuật toán tối ưu. Trong tối ưu topo chất lỏng, việc xác định các vùng cần thêm hoặc loại bỏ vật liệu đòi hỏi phân tích độ nhạy chi tiết.
3.2. Sử dụng phương pháp gradient và thuật toán tối ưu hóa
Các thuật toán tối ưu hóa dựa trên phương pháp gradient có thể được sử dụng để tìm kiếm hình dạng tối ưu. Tuy nhiên, cần chú ý đến vấn đề cực trị cục bộ và lựa chọn bước di chuyển phù hợp để đảm bảo sự hội tụ của thuật toán. Các thuật toán tối ưu hiện đại như các phương pháp quasi-Newton hoặc các thuật toán dựa trên học máy cũng có thể được sử dụng để cải thiện hiệu quả của quá trình tối ưu.
IV. Ứng dụng của Tối Ưu Hình Dạng Cánh Trong Kỹ Thuật Hàng Không
Tối ưu dạng trong cơ học chất lỏng mang đến nhiều ứng dụng thực tế thú vị, có thể kể tới như thiết kế hình dạng tối ưu của cánh máy bay nhằm cực tiểu lực cản của dòng không khí xung quanh trong công nghiệp hàng không vũ trụ [3], [4], thiết kế hình dạng động mạch chủ trong các nghiên cứu y học về phẫu thuật mạch vành ghép tim nhân tạo, được trình bày trong [5]–[7]. Trong [8], các tác giả nghiên cứu hình dạng tối ưu của vi sinh vật chuyển động trong một chất lỏng nhớt dưới tác động của sức căng bề mặt. Ứng dụng của bài toán tối ưu topology nhằm cực tiểu năng lượng tiêu tán của dòng chảy Stokes trong không gian hai chiều được trình bày trong [9]. Các bài toán tương tự cũng được xem xét cho dòng chảy Darcy-Stokes [10], [11] và dòng chảy Navier-Stokes [12], với các số Reynolds từ thấp tới trung bình [13]–[15].
4.1. Tối ưu hóa thiết kế cánh máy bay để giảm lực cản
Một trong những ứng dụng quan trọng nhất của tối ưu hóa hình dạng cánh là trong kỹ thuật hàng không. Việc thiết kế cánh máy bay với hình dạng tối ưu có thể giúp giảm lực cản, tăng lực nâng và cải thiện hiệu suất bay. Các phương pháp CFD và các thuật toán tối ưu có thể được sử dụng để tìm kiếm hình dạng cánh tối ưu cho các điều kiện bay khác nhau. Các nghiên cứu của [3], [4] đã trình bày chi tiết về các phương pháp này.
4.2. Tăng hiệu suất năng lượng trong thiết kế máy bay
Việc tối ưu hình dạng không chỉ giúp giảm lực cản mà còn có thể tăng hiệu suất năng lượng của máy bay. Bằng cách giảm tiêu thụ nhiên liệu, thiết kế tối ưu hóa hình dạng chất lỏng có thể giúp giảm chi phí vận hành và giảm tác động đến môi trường. Các ràng buộc trong tối ưu hóa cần được xem xét để đảm bảo tính khả thi của thiết kế. Việc áp dụng trí tuệ nhân tạo trong cơ học chất lỏng đang ngày càng được quan tâm để giải quyết các bài toán phức tạp này.
V. Nghiên cứu Tối Ưu Dạng Kênh Trộn và Bộ Trao Đổi Nhiệt
Các nghiên cứu khác về tối ưu dạng, tối ưu topology và các ứng dụng của chúng có thể tham khảo thêm trong [3], [16]–[19] và các tài liệu tham khảo trong đó. Mục đích chính của luận văn này là xây dựng một lược đồ giải số cho bài toán tối ưu dạng trong cơ học chất lỏng. Cho Ω là một miền mở, bị chặn trong không gian Rd ( d = 2 hoặc 3 ). Một dòng chảy chất lỏng đặt trên Ω được mô tả bởi hệ phương trình Stokes: − ν∆u + ∇ p = f trong Ω. Bài toán này đã được xem xét trong nhiều nghiên cứu trước đây [9], [15], [20]. Một cách hình thức, mô hình bài toán tối ưu dạng được phát biểu như sau: inf J ( Ω ).
5.1. Tối ưu hóa kênh trộn để cải thiện hiệu quả trộn
Tối ưu hóa kênh trộn là một lĩnh vực quan trọng trong nhiều ngành công nghiệp, từ hóa chất đến thực phẩm. Việc thiết kế kênh dẫn chất lỏng với hình dạng tối ưu có thể giúp cải thiện hiệu quả trộn, giảm thời gian trộn và tiết kiệm năng lượng. Các phương pháp CFD và các thuật toán tối ưu có thể được sử dụng để tìm kiếm hình dạng kênh trộn tối ưu cho các ứng dụng cụ thể. Cần phải xem xét các ràng buộc trong tối ưu hóa.
5.2. Thiết kế bộ trao đổi nhiệt hiệu quả hơn thông qua tối ưu hóa hình học
Việc tối ưu hóa hình học của bộ trao đổi nhiệt có thể giúp tăng hiệu quả trao đổi nhiệt, giảm kích thước và trọng lượng của thiết bị và tiết kiệm năng lượng. Các phương pháp CFD và các thuật toán tối ưu có thể được sử dụng để tìm kiếm hình dạng bộ trao đổi nhiệt tối ưu cho các điều kiện vận hành khác nhau. Việc ứng dụng tối ưu hóa hình dạng trong lĩnh vực này mang lại nhiều lợi ích kinh tế và môi trường.
VI. Kết luận và hướng nghiên cứu tiếp theo về tối ưu dạng
Luận văn được hoàn thành trong chương trình Thạc sĩ Khoa học Toán ứng dụng tại Viện Toán ứng dụng và Tin học, Đại học Bách Khoa Hà Nội dưới sự hướng dẫn của TS. Tạ Thị Thanh Mai. Chương 2 của luận văn đã được tóm tắt và công bố trên tạp chí Applied Numerical Mathematics số 126 tháng 4 năm 2018, trang 160 - 179 với tiêu đề "Shape optimization for Stokes flows using sensitivity analysis and finite element method". Một số kết quả trong chương 1 được trình bày trong nghiên cứu "Characteristic finite element method for natural convection problems" đã gửi đăng. Tác giả vẫn đang tiếp tục nghiên cứu mở rộng thuật toán này cho bài toán tối ưu dạng trong các mô hình cơ học khác, chẳng hạn như dòng đối lưu tự nhiên, cấu trúc đàn hồi hay phương trình Navier-Stokes. Mặc dù được hoàn thành với nhiều cố gắng nhưng do những hạn chế về thời gian và kinh nghiệm, luận văn này không thể tránh khỏi những sai sót. Tác giả rất mong nhận được những ý kiến đóng góp quý báu từ thầy cô và các bạn học viên để luận văn được hoàn thiện hơn nữa.
6.1. Tổng kết các kết quả đạt được về thiết kế tối ưu hóa
Luận văn đã trình bày một thuật toán giải số cho bài toán tối ưu dạng trong cơ học chất lỏng, dựa trên phương pháp biến phân Hadamard và phương pháp phần tử hữu hạn. Mã nguồn chương trình của các thuật toán và các thí nghiệm giải số được phát triển trên phần mềm Freefem++. Các kết quả hình ảnh sử dụng trong luận văn được hỗ trợ bởi các phần mềm mã nguồn mở gnuplot và Medit. Các ví dụ mô phỏng số đã chứng minh tính hiệu quả của thuật toán trong việc tìm kiếm hình dạng tối ưu cho các ứng dụng khác nhau.
6.2. Các hướng nghiên cứu tiềm năng về tối ưu hóa tham số
Các hướng nghiên cứu tiềm năng bao gồm mở rộng thuật toán cho các mô hình cơ học khác, chẳng hạn như dòng đối lưu tự nhiên, cấu trúc đàn hồi hay phương trình Navier-Stokes. Việc tích hợp các kỹ thuật học máy và trí tuệ nhân tạo cũng có thể giúp cải thiện hiệu quả và độ chính xác của quá trình tối ưu. Hơn nữa, việc nghiên cứu các phương pháp tối ưu hóa tham số có thể được áp dụng đồng thời với tối ưu hóa hình dạng, tăng cường tiềm năng thiết kế tối ưu hóa chất lỏng.
