Với những kết quả thu được đề tài có thể được ứng dụng vào thực tiễn nghiên cứu ảnh hưởng của hình dạng tới lực cản, từ đó tối ưu hóa trong việc thiết kế, chế tạo và sản xuất.
Có thể mở rộng và phát triển theo hướng sâu và rộng hơn như:
Không chỉ áp dụng cho đối tượng là các mô hình đơn giản dạng 2D như đã tính toán ở phần trước, mà có thể tính toán cho các đối tượng khác có liên quan tới vấn đề động học lưu chất như: xe gắn máy, tàu thủy, máy bay, ôtô...
Tính toán trong điều kiện phức tạp hơn, xây dựng mô hình và tính toán cho bài toán không gian ba chiều (bài toán 3D) và kiểm tra, so sánh trong trường hợp 3D và 2D.
Qua những hướng phát triển vừa nêu trên, nếu bài toán được thực hiện trong điều kiện và môi trường tốt hơn thì chắc hẳn nó sẽ có vị trí ở tầm cao hơn.
PHỤ LỤC
CÁC NHÓM CÔNG CỤ TRONG PHẦN MỀM GAMBIT VÀ FLUENT 1. Phần mềm Gambit
1.1. Các nhóm phần tử trong phần mềm
Khởi động GAMBIT: StartprogramFluent Inc Products Gambit 2.2 Giao diện làm việc chính của Gambit 2.2
Mô hình trong GAMBIT 2.4 bao gồm thanh menu, thanh Operation, thanh Global control và vùng mô tả lệnh.
Thanh Operation
Biểu tượng Mô tả
Geometry
Biểu tượng tạo điểm
Biểu tượng tạo đường thẳng Thanh Operation
Thanh Global control Vùng mô tả lệnh
Biểu tượng tạo mặt phẳng
Biểu tượng tạo hình khối
Nhóm hình học
Mesh
Biểu tượng chia lưới cho biên
Biểu tượng chia lưới trên đường
Biểu tượng chia lưới trên mặt
Biểu tượng chia lưới trên hình khối
Biểu tượng chia lưới nhóm hình học
Zones
Biểu tượng tạo điều kiện biên
Biểu tượng khai báo vật liệu
Tool
Tạo ra, sửa đổi, trình bày, những hệ thống tạo độ, lưới, thước đo
Tạo ra những hàm điều khiển các phần tử mắt lưới.
Tạo ra và đan lưới của những thể tích cùng nhau đại diện vùng luồng bao vây trong turbomachine
Biểu tượng xóa các cấu trúc hình học
Các đầu nối
Thanh Global control gồm các lệnh hiệu chỉnh màn hình làm việc của Gambit
Biểu tượng Mô tả
Tỷ lệ những hiển thị của màn hình phù hợp với những ranh giới.
Chỉ định vị trí của trục quay cho chuyển động của mô hình bằng phương pháp con chuột.
Thể hiện các cửa sổ đồ họa.
Tạo mủi tên, đường thẳng, text trong màn hình hiển thị
Chỉnh độ sáng hoặc tối của màn hình làm việc
1.2. Các lệnh sử dụng để xây dựng mô hình trong Gambit
Tạo điểm
Hiển thị mắt lưới
Các đặc trưng của màn hình hiển thị
Chọn màu của mô hình
4. Vùng đặt tên
1. Biểu tượng cấu trúc hình học
2. Biểu tượng tạo điểm 3. Vùng nhập tọa độ điểm
5. Nhấn Apply để xác lập
Tạo đường
Tạo mặt
1. Biểu tượng tạo mặt
2. Menu các lệnh tạo mặt 1. Biểu tượng tạo cạnh
Chia lưới trên đường
Chia lưới cho mặt
3. Kiểu lưới chia
1. Biểu tượng chia lưới cho mặt
2. Vùng khai báo các thông số chia
1. Biểu tượng chia lưới
2. Biểu tượng chia lưới trên đường
3. Vùng khai báo các thông số của đườngđược chia
Áp đặt điều kiện biên cho bài toán
Khai báo môi trường làm việc
1. Biểu áp đặt điều kiện biên 2. Vùng khai báo thông sốđiều kiện biên 1. Biểu tượng khai báo vật liệu 2. Vùng khai báo các thuộc tính
Xuất sang File lưới (.msh)
Xuất sang file lưới (Mesh)
Click chọn Export 2-D(x-y)
Click Accept lưu File
2. Các nhóm phần tử trong Fluent 6.2
Khởi động FLUENT: StartprogramFluent Inc Products Fluent 6.2
Lựa chọn giạng bài toán mô phỏng
Trong đó:
- 2d: Bộ giải chính xác 2 chiều đơn
- 2ddp: Bộ giải chính xác 2 chiều, chính xác kép - 3d: Bộ giải chính xác 3 chiều đơn
- 3ddp: Giải chính xác 3 chiều, chính xác kép
2.1. Menu File
- Read: Đọc những file .msh được lưu từ GAMBIT, và những file .cas, hoặc .dat là những kết quả được lưu trong FLUENT.
- Write: Lưu flie Fluent dạng đuôi “.cas”, “.dat”
- Import: Nhập dữ liệu vào từ các file : ABQUS, ANSYS, GAMBIT, CFX .... - Export: Xuất dữ liệu ra các file của : ABQUS, ANSYS imput, ACCII, CFX .... - Hardcopy: Xuất kết quả mô phỏng
2.2. Menu Gid
Hiển thị, kiểm tra lưới trước khi giải. - Check: Kiểm tra lưới
- Info: Cho biết tất cả những thông tin của lưới - Scale, translate, rotate: Tỷ lệ, di chuyển, quay lưới. - Smooth/swap: Có tác dụng làm nhắn bề mặt lưới.
2.3. Menu Define
Định nghĩa tất cả các điều kiện trước khi giải.
- Model: Định nghĩa mô hình tính toán trong Fluent. - Model: Định nghĩa mô hình tính toán trong Fluent. - Material: Định nghĩa thuộc tính vật chất.
- Oprating conditions: Chỉ định điều kiện hoạt động ban đầu.
- Boundary Conditions: Chỉ định những điều kiện ranh giới hay điều kiện biên của bài toán, các điều kiện cụ thể cho một phương trình vi phân, giúp giới hạn các kết quả ra của phương trình.
- Priodic Conditions: Định nghĩa điều kiện hoạt động có tính chu kỳ. - Grid interface : Định nghĩa bề mặt trong của lưới.
- Units: Định nghĩa đơn vị trong tính toán mô phỏng.
2.4. Menu Solve
- Control: Điều khiển quá trình giải. - Initialize: Khởi đầu quá trình giải.
- Monitor: Hiển thị kết quả các giá trị vận tốc , các hệ số,.. cho phép ta có thể in hoặc plot đồ thị - một cách hiển thị kết quả trong xử lý kết quả mô phỏng bằng FLUENT.
- Iterate: Bắt đầu quá trình giải.
2.5. Menu Adapt
Kỹ thuật được dùng để hoàn thiện chất lượng lưới chia, cho phép làm mịn hoặc làm thô lưới tính toán dựa trên mô hình hình học và dữ liệu số của phương pháp giải.
Ngoài ra, Fluent còn cung cấp các công cụ để tạo và quan sát miền lưới được thích ứng tuỳ theo các ứng dụng đặc biệt.
- Boundary Adapt: Biên khảo sát. - Region: Miền khảo sát.
2.6. Menu Suface
Trong menu này ta có thể tạo các khu vực, điểm, đường... trên mô hình, từ đó ta có thể lấy các giá trị áp suất, nhiệt độ, vận tốc..tại các điểm đó và hiển thị nó ra.
2.7. Menu Display
2.8. Menu plot
Xuất kết quả mô phỏng ra dạng đồ thị, biểu đồ, từ đó có thể so sánh với thực nghiệm bằng việc chồng đồ thị.
Ta có thể xuất ra những kết quả các giá trị của áp suất, vận tốc, nhiệt độ, rối loạn ...trong mô hình, từ đó ta có thể nhìn thấy trực quan.
2.9. Menu parallel
Xử lý song song trên máy tính, quá trình xử lý đồng thời sử dụng hơn 1 bộ xử lý hoặc một nhân xử lý để thực hiện chạy một chương trình hoặc tính toán đa luồng. Ý tưởng xử lý song song làm các chương trình chạy nhanh hơn bởi vì có nhiều phương tiện (Bộ xử lý hoặc nhân) để chạy các chương trình đó.
2.10. Menu Help
DANH MỤC CÁC BẢNG
Bảng 5.1: Điều kiện biên vận tốc vào (Velocity-inlet) ... 38
Bảng 5.2: Điều kiện biên áp suất ra (Pressure-outlet) ... 38
Bảng 5.3: Ảnh hưởng của mật độ lưới tới hệ số CD ... 39
Bảng 6.1: Hệ số lực cản CD của các vật thể có hình dạng khác nhau, kết quả thực nghiệm [5],[7] ... 46
Bảng 6.2: Điều kiện biên vận tốc vào (Velocity-inlet) ... 48
Bảng 6.3: Điều kiện biên áp suất ra (Pressure-outlet) ... 48
Bảng 6.4: Điều kiện biên vận tốc vào (Velocity-inlet) ... 53
Bảng 6.5: Điều kiện biên áp suất ra (Pressure-outlet) ... 53
Bảng 6.6: Điều kiện biên vận tốc vào (Velocity-inlet) ... 58
Bảng 6.7: Điều kiện biên áp suất ra (Pressure-outlet) ... 58
Bảng 6.8: Điều kiện biên vận tốc vào (Velocity-inlet) ... 62
Bảng 6.9: Điều kiện biên áp suất ra (Pressure-outlet) ... 62
Bảng 6.10: Điều kiện biên vận tốc vào (Velocity-inlet) ... 67
Bảng 6.11: Điều kiện biên áp suất ra (Pressure-outlet) ... 67
Bảng 6.12: Điều kiện biên vận tốc vào (Velocity-inlet) ... 72
Bảng 6.13: Điều kiện biên áp suất ra (Pressure-outlet) ... 72
Bảng 6.14: Điều kiện biên vận tốc vào (Velocity-inlet) ... 76
Bảng 6.15: Điều kiện biên áp suất ra (Pressure-outlet) ... 76
Bảng 6.16: Bảng so sánh kết quả tính toàn và hệ số thực nghiệm của các mô hình khác nhau trong cùng điều kiện với số 4 10 Re ... 81
DANH MỤC CÁC HÌNH ẢNH
Hình 2.1: Ô tô chuyển động trong không khí với vận tốc V ... 3
Hình 2.2: Nghiên cứu dòng lưu chất chuyển động qua khối trụ tròn bằng thực nghiệm ... 4
Hình 2.3: Mô hình nghiên cứu thực nghiệm dòng qua khối trụ tròn ... 4
Hình 2.4: Đường hầm khí động thí nghiệm ... 5
Hình 2.5: Đồ thị hệ số CD phụ thuộc vào số Re của một số vật thể có hình dạng khác nhau [7] ... 5
Hình 2.6: Hệ số lực cản của một số vật thể 2D [5] ... 6
Hình 2.7: Hệ số lực cản CD của một số vật thể dạng 3D [5] ... 7
Hình 2.8: Hệ số lực cản CD của ô tô [5] ... 7
Hình 2.9: Trường áp suất phân bố trên cánh và thân máy bay trực thăng ... 8
Hình 2.10: Mô phỏng hệ thống tuần hoàn máu trong y học ... 8
Hình 2.11: Mô phỏng dòng lưu chất qua mũ bảo hiểm của vận động viên. Công trình nghiên cứu của Đại Học Shefield, Anh ... 9
Hình 2.12: Mô phỏng dòng chảy qua cơ thể vận động viên bơi lội. Công trình nghiên cứu của công ty sản xuất đồ dùng thể thao hàng đầu thế giới ... 9
Hình 2.13: Mô phỏng ứng suất tiếp sinh ra trên cơ thể vận động viên bơi lôi. Công trình nghiên cứu của công ty sản xuất đồ dùng thể thao hàng đầu thế giới Speedo, Mỹ ... 10
Hình 2.14: Thử nghiệm sản phẩm xe đua thể thức 1 trong đường hầm khí động, tại phòng thí nghiệm mô phỏng khí động... 11
Hình 2.15: Xem xét đặc tính dòng không khí qua xe đua thể thức 1, mô phỏng bằng phần mềm Fluent ... 11
Hình 16: Xem xét dòng khí qua mô tô thể thao của hai nhà thiết kế Werner Seibert
người Đức và Robert Lewis người Anh, 2004 ... 12
Hình 2.17: Mô hình máy bay AEW&C trong ống khí động, Trung Quốc ... 12
Hình 2.18: Mô hình cánh máy bay trong phòng thí nghiệm ... 13
Hình 2.19: Phân bố áp suất trên cánh máy bay mô phỏng bằng phần mềm Fluent... 13
Hình 2.20: Phân bố áp suất tĩnh trên một nửa bề mặt máy bay ARJ21 (a), đồ thị so sánh kết quả thực nghiệm với kết quả tính toán bằng phần mềm Fluent, nghiên cứu của Zhu Jie, công ty hàng không thương mại, Trung Quốc [9] ... 13
Hình 2.21: Kết quả thu được từ mô phỏng buồng đốt sinh học bằng phần mềm Fluent ... 14
Hình 2.22: Mô phỏng bằng phần mềm Fluent sự phân bố áp suất trên bề mặt cánh (a); phân bố túi hơi trên bề mặt cánh (b) của thiết bị đẩy tàu thủy (chân vịt) ... 14
Hình 3.1: Hình minh họa lớp biên từ vùng chảy tầng sang vùng chảy rối ... 16
Hình 3.2: Lực cản và lực nâng trên cánh ... 17
Hình 3.3: Thành phần lực cản ma sát và lực cản áp suất trên bề mặt vật thể ... 17
Hình 3.4: Mốii quan hệ giữa lực cản áp suất và lực cản ma sát của một số vật thể, kết quả nghiên cứu thực nghiệm ... 18
Hình 4.1: Cấu trúc của phần mềm Fluent ... 25
Hình 4.2: Sơ đồ thuật toán trong phần mềm Fluent ... 27
Hình 4.3: Sơ đồ phương pháp giải lặp trong phần mềm Fluent ... 31
Hình 4.4: Rời rạc hóa (chia lưới) trong Gambit với khoảng 1000 phần tử lưới ... 33
Hình 5.1: Mô hình và điều kiện biên cho bài toán dòng qua khối trụ tròn ... 37
Hình 5.2: Rời rạc hóa miền không gian tính cho bài toán dòng qua khối trụ tròn ... 38
Hình 5.3: Đồ thị biểu diễn hệ số lực cản CD theo các mật độ lưới khác nhau ... 40
Hình 5.4: Phân bố trường vận tốc quanh khối trụ dưới dạng contours ... 41
Hình 5.5: Phân bố trường vận tốc quanh khối trụ dưới dạng vector ... 41
Hình 5.6: Phương của vector vận tốc trên bề mặt khối trụ tròn ... 42
Hình 5.7: Vị trí tách lớp biên của dòng khí ... 42
Hình 5.9: Đồ thị phân bố hệ số áp suất Cptại nửa dưới của khối trụ tròn ... 44
Hình 6.1: Mô hình cho bài toán dòng qua tấm phẳng ... 47
Hình 6.2: Rời rạc hóa miền không gian tính ... 48
Hình 6.3: Phân bố vận tốc quanh tấm phẳng dưới dạng contours ... 49
Hình 6.4: Phân bố vận tốc quanh tấm phẳng biểu diễn dưới dạng vector ... 49
Hình 6.5: Sự tạo xoáy của dòng khí phía sau tấm phẳng đứng dưới dạng vector .... 50
Hình 6.6: Hiện tượng dòng chuyển động ngược và tách rời lớp biên trên bề mặt tấm phẳng ... 50
Hình 6.7: Trường phân bố áp suất quanh tấm phẳng ... 51
Hình 6.8: Đồ thị phân bố hệ số áp suất Cp tại nửa trên của tấm phẳng ... 52
Hình 6.9: Đồ thị phân bố hệ số áp suất Cp tại nửa dưới của tấm phẳng ... 52
Hình 6.10: Mô hình bài toán nửa trụ tròn thuận dòng ... 53
Hình 6.11: Rời rạc hóa miền không gian tính ... 54
Hình 6.12: Phân bố trường vận tốc quanh khối trụ dưới dạng contours ... 54
Hình 6.13: Phân bố trường vận tốc quanh khối trụ dưới dạng vector... 55
Hình 6.14: Vị trí tách lớp biên trên bề mặt khối trụ ... 55
Hình 6.15: Phân bố trường áp suất quanh khối trụ dưới dạng contours ... 56
Hình 6.16: Đồ thị phân bố hệ số áp suất Cptại nửa trên của khối trụ ... 57
Hình 6.17: Đồ thị phân bố hệ số áp suất Cptại nửa dưới của khối trụ ... 57
Hình 6.18: Mô hình bài toán nửa khối trụ tròn ngược dòng ... 58
Hình 6.19: Rời rạc hóa miền không gian tính... 59
Hình 6.20: Phân bố trường vận tốc quanh khối trụ dưới dạng contours ... 59
Hình 6.21: Phân bố trường vận tốc quanh khối trụ dưới dạng vector... 60
Hình 6.22: Dòng chuyển động ngược và vị trí tách lớp biên trên bề mặt khối trụ .... ... 60
Hình 6.23: Phân bố trường áp suất quanh khối trụ dưới dạng contours ... 61
Hình 6.24: Đồ thị phân bố hệ số áp suất Cp tại nửa trên của khối trụ ... 61
Hình 6.26: Rời rạc hóa miền không gian tính cho bài toán ... 63
Hình 6.27: Phân bố trường vận tốc quanh khối trụ dưới dạng contours ... 63
Hình 6.28: Tách dòng khí tại vị trí đỉnh của tam giác... 64
Hình 6.29: Phân bố trường vận tốc quanh khối trụ dưới dạng vector... 64
Hình 6.30: Phân bố trường áp suất tĩnh quanh khối trụ tam giác thuận dòng ... 65
Hình 6.31: Đồ thị phân bố áp suất Cptại nửa dưới của khối trụ ... 66
Hình 6.32: Đồ thị phân bố áp suất Cptại nửa trên của khối trụ ... 66
Hình 6.33: Mô hình bài toán dòng qua khối trụ tam giác ngược dòng ... 67
Hình 6.34: Rời rạc hóa miền không gian tính... 68
Hình 6.35: Phân bố trường vận tốc quanh khối trụ dưới dạng contours ... 68
Hình 6.36: Phân bố trường vận tốc quanh khối trụ dưới dạng vector... 69
Hình 6.37: Dòng chuyển động ngược và tách lớp biên trên bề mặt khối trụ ... 69
Hình 6.38: Phân bố trường áp suất quanh khối trụ dưới dạng contours ... 70
Hình 6.39: Đồ thị phân bố hệ số áp suất Cptại nửa dưới của khối trụ ... 71
Hình 6.40: Mô hình bài toán dòng qua khối trụ vuông ... 71
Hình 6.41: Rời rạc hóa miền không gian tính... 72
Hình 6.42: Phân bố trường vận tốc quanh khối trụ vuông ... 73
Hình 6.43: Phân bố trường vận tốc quanh khối trụ dưới dạng vector... 73
Hình 6.44: Dòng chảy ngược và tách lớp biên trên bề mặt của khối trụ ... 74