đánh giá đặc tính khí động học của máy bay bình thuốc trừ sâu tích hợp

Bạn đang xem bản rút gọn của tài liệu. Xem và tải ngay bản đầy đủ của tài liệu tại đây (3.58 MB, 117 trang )

<span class="text_page_counter">Trang 1</span><div class="page_container" data-page="1">

---

NGUYỄN THANH HẢO

<b>ĐÁNH GIÁ ĐẶC TÍNH KHÍ ĐỘNG HỌC CỦA MÁY BAY BÌNH THUỐC TRỪ SÂU TÍCH HỢP </b>

Chun ngành: Kỹ Thuật Hàng KhơngMã số : 8520120

LUẬN VĂN THẠC SĨ

TP. HỒ CHÍ MINH, tháng 01 năm 2024

</div><span class="text_page_counter">Trang 2</span><div class="page_container" data-page="2">

Cán bộ hướng dẫn khoa học :

- Cán bộ hướng dẫn 1 : TS. Lưu Văn Thuần

- Cán bộ hướng dẫn 2 : TS. Nguyễn Song Thanh Thảo

Cán bộ chấm nhận xét 1 : PGS.TS Vũ Ngọc Ánh

Cán bộ chấm nhận xét 2 : PGS. TS Lê Tuấn Phương Nam

Luận văn thạc sĩ được bảo vệ tại Trường Đại học Bách Khoa, ĐHQG Tp. HCM ngày 27 tháng 01 năm 2024

Thành phần Hội đồng đánh giá luận văn thạc sĩ gồm:

(Ghi rõ họ, tên, học hàm, học vị của Hội đồng chấm bảo vệ luận văn thạc sĩ) 1. Chủ tịch : PGS. TS Ngô Khánh Hiếu

2. Thư ký : TS. Lê Thị Hồng Hiếu3. Phản biện 1 : PGS. TS Vũ Ngọc Ánh

4. Phản biện 2 : PGS. TS Lê Tuấn Phương Nam5. Ủy viên : TS. Phạm Minh Vương

Xác nhận của Chủ tịch Hội đồng đánh giá LV và Trưởng Khoa quản lý chuyên ngành sau khi luận văn đã được sửa chữa (nếu có).

<b>CHỦ TỊCH HỘI ĐỒNG TRƯỞNG KHOA KỸ THUẬT GIAO THÔNG </b>

</div><span class="text_page_counter">Trang 3</span><div class="page_container" data-page="3">

ĐẠI HỌC QUỐC GIA TP.HCM

<b>TRƯỜNG ĐẠI HỌC BÁCH KHOA </b>

<b>CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM Độc lập - Tự do - Hạnh phúc</b>

<b>NHIỆM VỤ LUẬN VĂN THẠC SĨ </b>

<b>Họ tên học viên: NGUYỄN THANH HẢO </b> MSHV: 2170161 Ngày, tháng, năm sinh: 06/09/1998 Nơi sinh: TP.HCM Chuyên ngành: Kỹ thuật Hàng không Mã số: 8520120

Tên tiếng Việt: Đánh giá đặc tính khí động học của máy bay bình thuốc trừ sâu tích hợp Tên tiếng Anh: Evaluation of aerodynamic characteristics of the integrated-pesticide tank drone

<b>II. NHIỆM VỤ VÀ NỘI DUNG:</b>

Thực hiện luận văn với mục tiêu nghiên cứu là: Mô phỏng và đánh giá đặc tính khí động của drone với bình thuốc trừ sâu tích hợp trong 3 chế độ bay (bay treo, bay lấy độ cao và bay tiến), thơng qua việc phân tích lực và phân bố của các trường áp suất, vận tốc với các tốc độ vịng quay của chong chóng khác nhau.

<b>III. NGÀY GIAO NHIỆM VỤ: 06/02/2023</b>

<b>IV. NGÀY HOÀN THÀNH NHIỆM VỤ: 10/12/2023</b>

<i><b>V. CÁN BỘ HƯỚNG DẪN </b></i>

Hướng dẫn 1: TS. Lưu Văn Thuần

Hướng dẫn 2: TS. Nguyễn Song Thanh Thảo

<b>TS. Lưu Văn Thuần TS. Nguyễn Song Thanh Thảo </b>

<b>TRƯỞNG KHOA KỸ THUẬT GIAO THÔNG </b>

(Họ tên và chữ ký)

</div><span class="text_page_counter">Trang 4</span><div class="page_container" data-page="4">

<b> LỜI CẢM ƠN </b>

<i>Luận văn thạc sĩ với đề tài “Đánh giá đặc tính khí động học của máy bay bình thuốc </i>

<i>trừ sâu tích hợp” nằm trong khối kiến thức Luận văn thạc sĩ của ngành Kỹ thuật Hàng </i>

không nhằm trang bị cho học viên sau đại học khả năng tự nghiên cứu, tìm hiểu và đề xuất giải pháp giải quyết các vấn đề liên quan đến ngành học.

Để hoàn thành luận văn này, trước hết tôi xin chân thành cảm ơn TS. Lưu Văn Thuần – Phụ trách Khoa Kỹ thuật Hàng không của Học viện Hàng không Việt Nam, là giảng viên hướng dẫn chính của tơi – người đã hướng dẫn, giải đáp và cung cấp thêm cho tôi các nền tảng kiến thức lý thuyết cơ sở để thực hiện luận văn cũng như chỉ ra các mặt hạn chế trong quá trình nghiên cứu để tơi hồn thiện kết quả của mình. Tiếp theo, tơi xin gửi lời cảm ơn đến TS. Nguyễn Song Thanh Thảo - Giảng viên Bộ môn Kỹ thuật Hàng khơng, Phó Khoa Kỹ thuật Giao thông, Trường Đại học Bách Khoa, Đại học Quốc gia Thành phố Hồ Chí Minh đã hướng dẫn và hỗ trợ tôi trong việc lựa chọn tên đề tài, hướng phát triển của đề tài, cách thức trình bày và hoàn thiện bài báo cáo luận văn.

Bên cạnh đó, tơi xin cảm ơn các nhóm nghiên cứu ngành Kỹ thuật Hàng không tại Học viện Hàng không VN và Trường Đại học Bách Khoa, Đại học Quốc gia TPHCM đã hỗ trợ tôi trong việc thảo luận, tìm hiểu lỗi để xử lý bài tốn mơ phỏng một cách tốt nhất. Nhờ đó, tơi có thể hiểu và hoàn thành luận văn như hiện nay. Cuối cùng, tôi xin chân thành cảm ơn các thầy cô trong Bộ môn Kỹ thuật Hàng không - Khoa Kỹ thuật Giao thông, đặc biệt là TS. Lê Thị Hồng Hiếu – Trưởng Bộ môn Kỹ thuật Hàng không và TS. Vương Thị Hồng Nhi, đã truyền đạt cho tôi những kiến thức và tài liệu quý báu về khí động lực học và CFD trong suốt quá trình học tập cao học 2 năm qua.

Do kiến thức và kinh nghiệm của bản thân còn hạn chế, trong quá trình thực hiện luận văn có thể khơng tránh khỏi những sai sót. Vì vậy, rất mong nhận được những ý kiến đóng góp của thầy cơ Hội đồng để tôi bổ sung và trau đổi thêm kiến thức về lĩnh vực này.

Tôi xin chân thành cảm ơn!

</div><span class="text_page_counter">Trang 5</span><div class="page_container" data-page="5">

<b> TÓM TẮT </b>

Với sự phát triển của công nghệ kỹ thuật trong lĩnh vực nông nghiệp, nhu cầu sử dụng máy bay không người lái ngày cảng trở nên phổ biến hơn. Hiện nay, trên thị trường, máy bay không người lái đã được sử dụng để phục vụ công tác phun thuốc trừ sâu cho các đồng ruộng với bình thuốc được gắn dưới thân. Nhằm mang lại những cải tiến đáng kể cho nền nông nghiệp với công nghệ máy bay khơng người lái, nhóm sinh viên bộ mơn Kỹ thuật hàng không của Trường Đại học Bách Khoa, Đại học Quốc gia TPHCM đã đưa ra mơ hình máy bay khơng người lái có tích hợp bình thuốc vào bên trong thân máy bay và đã đạt được một số thành tựu: đưa ra được thiết kế, phân tích tĩnh kết cấu và từ đó đánh giá tính khả thi của việc tích hợp bình thuốc vào máy bay, mơ phỏng số động lực học (CFD) cho mơ hình máy bay này. Trong bày luận văn này, tác giả tiếp tục thực hiện việc đánh giá đặc tính khí động của máy bay thông qua mô phỏng số động lực học với các chế độ bay (bay treo, bay lấy độ cao và bay tiến) bằng phần mềm Ansys Fluent. Thông qua quá trình phân tích lực và phân bố của các trường áp suất, vận tốc với các tốc độ vòng quay của chong chóng khác nhau, các kết quả nghiên cứu sẽ được thu thập, và có thể sử dụng để chỉ ra ưu điểm về mặt khí động trong thiết kế mơ hình máy bay cùng các bài tốn liên quan đến hệ thống phun ở các nghiên cứu sau này.

</div><span class="text_page_counter">Trang 6</span><div class="page_container" data-page="6">

<b> ABSTRACT </b>

With the development of technology and techniques in the agricultural field, the demand for unmanned aerial vehicles (UAV) has become more popular. Nowadays, the quadcopters with the removable-pesticide tanks have been usually used for spraying pesticides on fields. To bring significant improvements to agriculture with unmanned aircraft technology, a group of students from the Aerospace Engineering Department of Ho Chi Minh City University of Technology, Vietnam National University of Ho Chi Minh City, has launched a model of an unmanned aerial vehicle with an integrated tank inside the fuselage and has been achieved several achievements: providing the design and the static analysis of the structure; from there evaluating the feasibility of integrated-pesticides tank into the drone, simulating by computational fluids dynamics (CFD) method for this aircraft model. In this thesis, the author continues to evaluate the aerodynamic characteristics of the integrated-pesticides tank drone through numerical simulation in different flight modes (including hovering, climbing, and forward flight) by using Ansys Fluent software. Research results will be collected by analyzing force and distribution of pressure and velocity fields at different rotational speeds. They can be used to point out the advantages of aerodynamics in this model design and the problems related to pesticide spraying systems in future research.

</div><span class="text_page_counter">Trang 7</span><div class="page_container" data-page="7">

<b> LỜI CAM ĐOAN </b>

Tôi xin cam đoan đây là bài làm luận văn độc lập của riêng tơi, do chính tơi thực hiện dưới sự hướng dẫn của GVHD – TS. Lưu Văn Thuần và TS. Nguyễn Song Thanh Thảo. Các nội dung, tài liệu sử dụng trong đề tài này có nguồn gốc xác thực và được trích dẫn đầy đủ. Các kết quả phân tích trong luận văn sẽ được trình bày trung thực, không sao chép từ bất kỳ đề tài nghiên cứu nào.

Tôi sẽ chịu trách nhiệm về việc làm của mình!

<b>Ngày 10 tháng 12 năm 2023 Học viên thực hiện </b>

Nguyễn Thanh Hảo

</div><span class="text_page_counter">Trang 8</span><div class="page_container" data-page="8">

<b>MỤC LỤC </b>

NHIỆM VỤ LUẬN VĂN THẠC SĨ ... i

LỜI CẢM ƠN ...ii

1.3. Các nghiên cứu liên quan ... 2

1.4. Mục tiêu nghiên cứu ... 6

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

1.6. Đối tượng và phạm vi nghiên cứu ... 7

1.7. Cấu trúc luận văn ... 7

Chương 2. CƠ SỞ LÝ THUYẾT ... 9

2.1. Giới thiệu máy bay UAV (Drone) ... 9

2.2. Phương pháp tính tốn mơ phỏng số động lực học lưu chất (Computational Fluid Dynamics - CFD) ... 10

2.3. Lưới trong CFD ... 15

2.4. Mơ phỏng dịng rối ... 22

2.5. Xử lý gần tường ... 25

2.6. Các giải thuật tính tốn trong Ansys Fluent ... 28

Chương 3. THIẾT LẬP BÀI TỐN MƠ PHỎNG 1 CHONG CHĨNG ĐƠN ... 31

3.1. Xây dựng miền hình học ... 31

</div><span class="text_page_counter">Trang 9</span><div class="page_container" data-page="9">

3.2. Xây dựng lưới ... 34

3.3. Thiết lập bài tốn mơ phỏng ... 38

3.4. Phân tích kết quả mơ phỏng ... 42

Chương 4. THIẾT LẬP BÀI TỐN MƠ PHỎNG TỒN BỘ MƠ HÌNH DRONE 52 4.1. Xây dựng miền hình học ... 52

4.2. Xây dựng lưới ... 55

4.3. Thiết lập bài tốn mơ phỏng ... 57

Chương 5. PHÂN TÍCH KHÍ ĐỘNG HỌC CỦA DRONE KHI HOẠT ĐỘNG ... 70

5.1. Kết quả mô phỏng trường hợp bay treo ... 70

5.2. Kết quả mô phỏng trường hợp bay lấy độ cao ... 77

5.3. Kết quả mô phỏng trường hợp bay tiến ... 81

5.4. Nhận xét ... 88

Chương 6. KẾT LUẬN VÀ HƯỚNG PHÁT TRIỂN ... 90

6.1. Kết luận ... 90

6.2. Hạn chế của đề tài ... 90

6.3. Hướng phát triển của đề tài ... 91

TÀI LIỆU THAM KHẢO ... 92

LÝ LỊCH TRÍCH NGANG ... 100

</div><span class="text_page_counter">Trang 10</span><div class="page_container" data-page="10">

<b> DANH MỤC HÌNH ẢNH </b>

Hình 1.1 Ví dụ máy bay phun thuốc trừ sâu trên thị trường [1] ... 1

Hình 1.2 Dịng chuyển động của lưu chất đối với hệ quy chiếu quay (bên trái) và phân bố dịng rối (bên phải) [3] ... 3

Hình 1.3 Phân bố trường vận tốc theo 3 hướng của mặt cắt YZ [6] ... 4

Hình 1.4 Kết quả tính trường dòng của UAV với các độ cao treo khác nhau: (a) 2 m, (b) 3 m, (c) 5 m, và (d) độ cao không giới hạn tương đối [7] ... 5

Hình 2.1 Cơ chế hoạt động của máy bay UAV 4 chong chóng: ... 9

Hình 2.2 Các lực tạo ra trên drone: (a) Bay treo; (b) Bay tiến [11] ... 10

Hình 2.3 Một vài ứng dụng thực tiễn của CFD [15] ... 11

Hình 2.4 Ví dụ về quy trình các bước mơ phỏng trong Ansys Fluent ... 15

Hình 2.5 Thành phần cơ bản của lưới [17] ... 16

Hình 2.6 Phần tử lưới phổ biến (a) dạng 2D; (b) dạng 3D [17] ... 16

Hình 2.7 (a) Lưới cấu trúc; (b) Lưới phi cấu trúc [18] ... 17

Hình 2.8 Lưới hỗn hợp (lưới lai) [19] ... 18

Hình 2.9 Chất lượng phần tử lưới dựa trên tiêu chuẩn AR ... 20

Hình 2.10 Ví dụ về các vectơ trong công thức đánh giá độ khơng trực giao [20] ... 21

Hình 2.11 Ví dụ về 2 phần tử lưới liền kề [17] ... 21

Hình 2.12 Các phương pháp mơ phỏng dịng rối phổ biến [23] ... 23

Hình 2.13 Các khu vực gần tường theo y<sup>+</sup> [24] ... 26

Hình 2.14 Độ cao của phần tử lưới đầu tiên so với tường [25] ... 27

Hình 2.15 Sơ đồ khối bộ giải dựa trên áp suất [26] ... 29

Hình 3.1 Thơng số chong chóng T-Motor G28-9.2 từ nhà sản xuất ... 31

Hình 3.2 Thơng số mơ hình một chong chóng ... 32

Hình 3.3 Miền xoay chong chóng ... 33

Hình 3.4 Miền tĩnh chong chóng... 34

Hình 3.5 Biểu đồ sự độc lập giữa nghiệm và lưới ... 36

</div><span class="text_page_counter">Trang 11</span><div class="page_container" data-page="11">

Hình 3.6 Lưới miền tĩnh và miền quay chong chóng ... 37

Hình 3.7 Lưới lớp biên gần bề mặt chong chóng ... 37

Hình 3.8 Thơng số thiết lập điều kiện biên “pressure oulet” ... 39

Hình 3.9 Điều kiện biên của các mặt miền chính trường hợp mơ phỏng bay treo ... 40

Hình 3.10 Biểu đồ sai số Residuals bài tốn mơ phỏng chong chóng đơn với tốc độ 3000 vịng/phút ... 42

Hình 3.11 Lưu lượng đầu ra và vào của bài toán mơ phỏng chong chóng đơn với tốc độ 3000 vịng/phút ... 43

Hình 3.12 Lực đẩy do chong chóng tạo ra khi quay với tốc độ 3000 vòng/phút ... 43

Hình 3.13 Đường đặc tính vận hành của chong chóng ... 44

Hình 3.14 Phân bố vận tốc ở mặt phẳng YZ qua tâm chong chóng đơn với tốc độ vịng quay 3000 vịng/phút: (a) tồn bộ miền tính tốn; (b) một phần miền tính tốn ... 45

Hình 3.15 Lý thuyết động lượng qua một chong chóng [32] ... 45

Hình 3.16 Phân bố vận tốc ở mặt phẳng YZ qua tâm chong chóng đơn với tốc độ vịng quay 3000 vòng/phút khi chia lưới mịn hơn ở bề mặt chong chóng ... 47

Hình 3.17 Phân bố vận tốc ở mặt phẳng YZ qua tâm chong chóng đơn với tốc độ vòng quay 3000 vòng/phút đối với miền tính tốn mở rộng xuống dưới ... 48

Hình 3.18 Kết quả mơ phỏng trường vận tốc: (a) Chong chóng đơn với chiều cao miền tính tốn 5D [33]; (b) Chong chóng đơn với chiều cao miền tính tốn 7,4 D [30] ... 49

Hình 3.19 Kết quả thực nghiệm và mơ phỏng CFD đối với 1 mơ hình drone ở cùng 1 khoảng cách: (a) thực nghiệm; (b) mô phỏng ... 49

Hình 3.20 Phân bố áp suất tại mặt phẳng chứa chong chóng khi bay treo với tốc độ 3000 vịng/phút ... 50

Hình 3.21 Sự thay đổi áp suất qua chong chóng theo phương Y ... 50

Hình 4.1 Mơ hình Drone sử dụng để mơ phỏng ... 52

Hình 4.2 Kết cấu chi tiết của thân drone ... 53

Hình 4.3 (a) Kích thước miền xoay; (b) Kích thước miền tĩnh tham khảo nghiên cứu của Mauro Ghirardelli ... 54

</div><span class="text_page_counter">Trang 12</span><div class="page_container" data-page="12">

Hình 4.4 Kích thước miền tĩnh thu hẹp ... 54

Hình 4.5 Phân bố vận tốc tại mặt phẳng YZ đi qua tâm chong chóng với tốc độ quay 3000 vịng/phút: (a) Vùng hình học rộng; (b) Vùng hình học thu hẹp ... 55

Hình 4.6 Lưới trên bề mặt thân máy bay và vùng quay chong chóng ... 56

Hình 4.7 Lưới xung quang chong chóng... 57

Hình 4.8 Lưới xung quang phần thân máy bay ... 57

Hình 4.9 Điều kiện biên của các mặt trong trường hợp bay treo ... 59

Hình 4.10 Điều kiện biên của các mặt trong trường hợp bay lấy độ cao ... 62

Hình 4.11 Góc nghiêng của drone khi bay tới ... 65

Hình 4.12 Điều kiện biên của các mặt trong trường hợp bay tiến ... 65

Hình 5.1 Biểu đồ sai số Residuals bài tốn mơ phỏng bay treo của drone ở tốc độ vịng quay 3000 vịng/phút ... 70

Hình 5.2 Tổng lực đẩy do 4 chong chóng tạo ra ở tốc độ quay 3000 vịng/phút ... 71

Hình 5.3 Lưu lượng đầu ra và vào của bài tốn mơ phỏng bay treo ... 71

Hình 5.4 Lực đẩy và lực cản tác dụng drone theo số vòng quay ... 72

Hình 5.5 Phân bố áp suất trên bề mặt 4 chong chóng khi bay treo tốc độ 3000 rpm .... 73

Hình 5.6 Phân bố áp suất trên bề mặt chong chóng khi bay treo tốc độ 3000 vịng/phút: (a) mặt trên; (b) mặt dưới ... 74

Hình 5.7 Phân bố áp suất tại mặt cắt ngang chong chóng khi bay treo tốc độ 3000 vịng/phút ... 74

Hình 5.8 Phân bố áp suất tại thân máy bay khi bay treo tốc độ 3000 vịng/phút: ... 75

Hình 5.9 Phân bố vận tốc tại mặt phẳng YZ đi qua tâm chong chóng: ... 76

Hình 5.10 Phân bố vận tốc của khơng khí qua drone nhìn từ trên xuống: ... 77

Hình 5.11 Biểu đồ sai số Residuals bài toán bay lấy độ cao với vận tốc 1 m/s ... 78

Hình 5.12 Tổng lực đẩy do 4 chong chóng tạo ra ở khi bay lấy độ cao với vận tốc 1 m/s ... 78

Hình 5.13 Lưu lượng ra vào của bài tốn mơ phỏng bay lấy độ cao với vận tốc 1 m/s 78 Hình 5.14 Lực cản trên thân máy bay theo tốc độ lấy độ cao ... 79

</div><span class="text_page_counter">Trang 13</span><div class="page_container" data-page="13">

Hình 5.15 Phân bố áp suất trên bề mặt chong chóng khi bay lấy độ cao với tốc độ ... 80 Hình 5.16 Phân bố vận tốc tại mặt phẳng chứa chong chóng với vận tốc lấy độ cao 3 m/s: (a) mặt cắt YZ, nhìn trực diện; (b) mặt cắt ZX, nhìn từ trên xuống ... 81 Hình 5.17 Phân bố trường vận tốc quanh drone khi bay lấy độ cao với vận tốc 3 m/s . 81 Hình 5.18 Biểu đồ sai số Residuals bài tốn mơ phỏng bay tiến với vận tốc 5 m/s, góc nghiêng 5 độ ... 82 Hình 5.19 Lưu lượng đầu ra và vào của bài tốn mơ phỏng bay tiến với vận tốc 5 m/s, góc nghiêng 5 độ ... 82 Hình 5.20 Tổng lực đẩy chong chóng với vận tốc 5 m/s, góc nghiêng 5 độ: ... 83 Hình 5.21 Phân bố trường vận tốc quanh drone với vận tốc bay tiến 5 m/s và góc nghiêng 5 độ: (a) nhìn nghiêng; (b) nhìn từ trên xuống; (c) nhìn mặt sau ... 86 Hình 5.22 Phân bố áp suất trên thân khi bay tiến với góc nghiêng 5 độ: (a) vận tốc bay 5 m/s; (b) vận tốc bay 8 m/s ... 87 Hình 5.23 Phân bố vận tốc tại mặt XY đi qua tâm máy bay với vận tốc 5 m/s, ở các góc nghiêng: (a) 5 độ, (b) 15 độ ... 87 Hình 5.24 Phân bố vận tốc tại mặt XY đi qua tâm máy bay với vận tốc 8 m/s và góc nghiêng 5 độ ... 88

</div><span class="text_page_counter">Trang 14</span><div class="page_container" data-page="14">

<b> DANH MỤC BẢNG BIỂU </b>

Bảng 2.1 Chất lượng phần tử lưới dựa trên các khoảng giá trị độ cân xứng lưới [20] .. 19

Bảng 2.2 Chất lượng phần tử lưới dựa trên các khoảng giá trị độ không trực giao [20] ... 21

Bảng 2.3 Bảng phân loại mơ hình rối theo số phương trình vận chuyển ... 23

Bảng 3.1 Kích thước miền tĩnh chong chóng ... 33

Bảng 3.2 Thông số chi tiết chia lưới lớp biên trên chong chóng ... 35

Bảng 3.3 Lực đẩy chong chóng theo mật độ lưới ... 36

Bảng 3.4 Thông số chất lượng lưới bài tốn mơ phỏng một chong chóng ... 37

Bảng 3.5 Bảng thiết lập điều kiện giải ... 40

Bảng 3.6 Các giá trị khởi tạo trong trường hợp bay treo ... 41

Bảng 3.7 Lực đẩy chong chóng theo tốc độ quay ... 43

Bảng 3.8 Kích thước miền tĩnh chong chóng mở rộng xuống dưới ... 47

Bảng 4.1 Thông số và chất lượng lưới bài tốn mơ phỏng drone ... 56

Bảng 4.2 Thiết lập bài tốn mơ phỏng ... 58

Bảng 4.3 Phương pháp giải bài tốn mơ phỏng ... 58

Bảng 4.4 Thiết lập điều kiện biên trường hợp mô phỏng bay treo ... 60

Bảng 4.5 Các giá trị khởi tạo trường hợp bay treo ... 61

Bảng 4.6 Thiết lập điều kiện biên trường hợp mô phỏng bay lấy độ cao ... 62

Bảng 4.7 Các giá trị khởi tạo trường hợp bay lấy độ cao với các vận tốc khác nhau .... 63

Bảng 4.8 Thiết lập điều kiện biên trường hợp mơ phỏng bay tiến với các góc nghiêng khác nhau, vận tốc 5 m/s ... 66

Bảng 4.9 Thiết lập điều kiện biên trường hợp mô phỏng bay tiến với các góc nghiêng khác nhau, vận tốc 8 m/s ... 68

Bảng 5.1 Tổng lực đẩy tạo ra bởi chong chóng và lực cản trên thân theo tốc độ quay . 71 Bảng 5.2 Tốc độ quay của chong chóng và lực cản trên thân khi bay lấy độ cao ... 78

</div><span class="text_page_counter">Trang 15</span><div class="page_container" data-page="15">

Bảng 5.3 Lực đẩy phương thẳng đứng và phương ngang theo góc nghiêng khi bay tiến với vận tốc 5 m/s ... 84 Bảng 5.4 Lực đẩy phương thẳng đứng và phương ngang theo góc nghiêng khi bay tiến với vận tốc 8 m/s ... 84 Bảng A.1 Bảng số liệu chong chóng T-motor G28-9.2 của RCcopter ... 96 Bảng A.2 Bảng số liệu chong chóng T-motor G28-9.2 của Artcopter ... 97 Bảng B.1 Bảng tính tốn sai số kết quả mô phỏng lực đẩy bằng Ansys Fluent so với kết quả đo thử nghiệm của Artcopter và Rccopter ... 98

</div><span class="text_page_counter">Trang 16</span><div class="page_container" data-page="16">

<b> BẢNG CHỮ VIẾT TẮT </b>

<b>STT CHỮ VIẾT TẮT GIẢI NGHĨA </b>

(Phương pháp sai phân hữu hạn)

(Phương pháp thể tích hữu hạn)

(Phương pháp Lattice Boltzmann)

(Mơ phỏng xốy lớn)

(Nhiều khung tham chiếu)

</div><span class="text_page_counter">Trang 17</span><div class="page_container" data-page="17">

14 UAV <sup>Unmanned Aerial Vehicle </sup>

(Phương tiện bay/Máy bay không người lái)

(Hệ thống máy bay không người lái)

16 WALE <sup>Wall-adapting Local-eddy Viscosity </sup>(Độ nhớt xoáy cục bộ thích ứng với tường)

</div><span class="text_page_counter">Trang 18</span><div class="page_container" data-page="18">

<b>Chương 1. TỔNG QUAN ĐỀ TÀI 1.1. Giới thiệu </b>

Với sự bùng nổ dân số trên thế giới hiện nay, việc đầu tư phát triển nền nông nghiệp năng suất cao, hiệu suất cao và bền vững ngày càng trở nên quan trọng hơn. Trong môi trường hiện đại, nông nghiệp là điều cần thiết cho sự tồn tại của hơn 60 phần trăm tổng dân số thế giới. Như vậy, hiện đại hóa nơng nghiệp được coi là điều tất yếu khi mà nhu cầu và nguồn cung lương thực tăng lên. Và máy bay không người lái là một trong những thiết bị có lợi nhất cho nông nghiệp hiện đại, đặc biệt là ứng dụng của nó trong việc phun thuốc trừ sâu cho đồng ruộng bởi lẽ phun thuốc trừ sâu bằng tay tồn tại rất nhiều hạn chế như tốn nhiều thời gian, công sức đồng thời có nguy cơ gây ra khối u, mẫn cảm, dị ứng và các bệnh khác ở người [1].

<i>Hình 1.1 Ví dụ máy bay phun thuốc trừ sâu trên thị trường [1] </i>

Điều đáng nói là khi sử dụng máy bay không người lái, người dùng chỉ cần đứng yên một chỗ và điều khiển máy từ xa dựa theo bản đồ đã được lập trình sẵn. Cụ thể hơn, ưu điểm của việc phun thuốc trừ sâu bằng máy bay có thể kể đến như:

− Máy có khả năng tự động cao với hệ thống phun cực kỳ chính xác và đồng đều; − Thời gian làm việc nhanh và hiệu quả;

− Tiết kiệm được chi phí nhân cơng và nhiên liệu; − An toàn cho người dùng và môi trường;

− Sử dụng được trên mọi địa hình khác nhau.

</div><span class="text_page_counter">Trang 19</span><div class="page_container" data-page="19">

<b>1.2. Đặt vấn đề </b>

Một trong các máy bay không người lái được sử dụng phổ biến trên thị trường hiện nay để phun thuốc trừ sâu là quadcopter. Quadcopters hoạt động bằng cách sử dụng chong chóng, là thiết bị có thể biến chuyển động quay thành lực đẩy tuyến tính, theo thứ tự, tạo ra lực nâng theo hướng của trục quay. Chong chóng máy bay không người lái cung cấp lực nâng cho máy bay bằng cách quay và tạo ra luồng không khí thơng qua sự chênh lệch áp suất giữa bề mặt trên và dưới của chong chóng [2]. Thiết kế hiện tại của chong chóng được sử dụng trong máy bay không người lái (UAV) sẽ được mô phỏng bằng công cụ Ansys Fluent và dữ liệu kết quả sẽ được sử dụng cho mục đích đánh giá đặc tính khí động.

Khi xét đến phần thân của máy bay, thông thường đối với máy bay UAV phun thuốc trừ sâu trên thị trường, bình thuốc sẽ được gắn dưới phần thân máy bay. Nhận thấy việc tích hợp bình thuốc vào thân của drone có thể mang lại ưu điểm về thiết kế hình học và lợi ích về đặc tính khí động, thiết kế drone với mơ hình có bình thuốc trừ sâu tích hợp vào thân đã được đề xuất, từ đó xem xét, đánh giá đặc tính khí động của mơ hình này.

Cuối cùng, tác động của mơ phỏng hình dạng của thiết kế mới này sẽ được kiểm tra kỹ lưỡng về mặt khí động và những thay đổi, thay thế (nếu có) sẽ được đề xuất để tối ưu hóa hiệu suất khí động học của máy bay, nhằm mang lại lợi ích cao cho việc ứng dụng vào công việc phun thuốc trừ sâu của nền nông nghiệp trên thực tế.

<b>1.3. Các nghiên cứu liên quan </b>

Đề tài nghiên cứu của Konstantinos Christodoulou cùng 5 cộng sự của ông tập trung vào việc phân tích khí động học của 1 chong chóng máy bay khơng người lái loại 4 cánh (Quadcopter) bằng phương pháp CFD. Cụ thể là, một thiết kế chong chóng được mô phỏng cho máy bay quadcopter để nghiên cứu hiệu suất của nó theo các phương pháp mơ hình hóa và các mơ hình dịng rối khác nhau. Mơ hình độ nhớt xốy tuyến tính SST k-ω, mơ hình BSL-RSM và mơ hình dịng rối đại số tường minh đã được sử dụng cho thí nghiệm này. Nghiên cứu cho thấy khơng có sự khác biệt đáng kể nào giữa kết quả mà

</div><span class="text_page_counter">Trang 20</span><div class="page_container" data-page="20">

các mơ hình dịng rối cung cấp liên quan đến đặc tính khí động học của chong chóng và sự phân bố áp suất xung quanh chong chóng [3].

<i>Hình 1.2 Dịng chuyển động của lưu chất đối với hệ quy chiếu quay (bên trái) và phân bố dòng rối (bên phải) [3] </i>

Một bài nghiên cứu khác về Mơ hình khí động học của các phương tiện quadcopter loại nhỏ, Seokkwanyoon cùng 4 nhà nghiên cứu đã thực hiện một mô phỏng tính tốn có độ chính xác cao để quan sát các tương tác khí động học giữa các chong chóng và thân của quadcopter đó. Phương pháp mơ phỏng được sử dụng thơng qua các phương trình Navier-Stoke khơng ổn định bằng các sơ đồ bậc cao, số Mach thấp và mơ hình dịng rối hỗn hợp. Nghiên cứu chỉ ra rằng các kết quả tính tốn cho thấy sự gắn kết cao giữa mô phỏng với dữ liệu thực nghiệm. Ngồi ra, bài nghiên cứu cịn so sánh giữa các chong chóng được gắn ở mặt trên so với thân và mặt dưới so với thân, cho thấy kết quả là chong chóng được gắn ở mặt dưới thân khơng tốt hơn chóng chóng được lắp ở mặt trên như thường lệ [4]. Tại Trường Đại học Vigo, Tây Ban Nha, C. Paz và các đồng nghiệp của ông đã thực hiện một nghiên cứu về đánh giá phương pháp mô phỏng CFD đối với máy bay UAV 4 chong chóng DJI Phantom 3. Mục đích của nghiên cứu này là đánh giá hiệu suất mơ phỏng dịng chảy 3D thông qua 2 phương pháp: phương pháp Nhiều khung tham chiếu (MRF) và Lưới trượt. Bên cạnh đó, tác giả cũng đánh giá ảnh hưởng của khoảng cásch drone đến mặt đất để xem xét độ cao bay của drone. Kết quả thu được cho thấy vị trí tương đối giữa cánh quạt và khung được chứng minh là yếu tố then chốt, bởi lẽ khi so sánh giữa 2 phương pháp trong trường hợp mô phỏng nguyên drone, tốc độ thay đổi lực

</div><span class="text_page_counter">Trang 21</span><div class="page_container" data-page="21">

đẩy tương tự đã đạt được khi giảm thiểu sự chồng chất của chong chóng lên các cánh tay drone trong trường hợp mô phỏng bằng phương pháp MRF. Tuy nhiên, cường độ lực đẩy khác nhau ít nhất 11% giữa các vị trí mơ phỏng bất kỳ. Mặc dù có độ lệch này, chi phí tính tốn bằng MRF thấp hơn đáng kể, do đó phương pháp này trở thành một lựa chọn được ưa chuộng [5].

Quan tâm đến việc ứng dụng drone trong lĩnh vực nông nghiệp, một nghiên cứu nổi bật trong lĩnh vực máy bay UAV phun thuốc trừ sâu là nghiên cứu về động học và phân bố trường dịng chuyển động xi của máy bay không người lái sáu cánh quạt SLK-5, được thực hiện bởi Yang Fengbo và nhóm nghiên cứu. Nghiên cứu sử dụng phương pháp mơ phỏng trung bình RANS (Reynolds-averaged Navier-Stokes) ở chế độ bay treo, với mô hình rối RNG <i>k - ε và kỹ thuật lưới động. Một số kết luận thu được là: vận tốc hướng </i>

z chủ yếu nằm trong quy luật phân bố vận tốc của dịng khí xi dịng; giá trị vận tốc tối đa của dịng chuyển động xi gần 10 m/s; sai số tương đối nhỏ hơn 12% giữa kết quả vận tốc hướng z thu được từ thí nghiệm vận tốc gió và từ mơ phỏng tại các điểm được đánh dấu; và chiều cao khuyến nghị cho UAV khi hoạt động là 0,6 m để tránh ảnh hưởng của cánh quạt khi phun thuốc [6].

<i>Hình 1.3 Phân bố trường vận tốc theo 3 hướng của mặt cắt YZ [6] </i>

</div><span class="text_page_counter">Trang 22</span><div class="page_container" data-page="22">

Để nghiên cứu ảnh hưởng của trường dịng chảy xi của UAV đến các hoạt động phun thực tế, Yongjun Zheng và các cộng sự đã xây dựng mơ hình 3D của máy bay UAV nông nghiệp sáu chong chóng JF01-10. Qua q trình mơ phỏng và phân tích trường dịng xi ở các độ cao bay treo khác nhau, có thể kết luận rằng hiệu quả vận hành phụ thuộc rất nhiều vào độ cao bay treo và độ cao bay treo phù hợp cho máy bay không người lái JF01-10 được coi là 3 m trên thực tế [7].

<i>Hình 1.4 Kết quả tính trường dịng của UAV với các độ cao treo khác nhau: (a) 2 m, (b) 3 m, (c) 5 m, và (d) độ cao không giới hạn tương đối [7] </i>

Nghiên cứu của tác giả Haiyan Zhang và Yubin Lan cùng các nhà khoa học khác, mơ phỏng trường dịng chuyển động xi của UAV bốn chong chóng P20 ở các tốc độ bay khác nhau, sử dụng Phương pháp Lattice Boltzmann (LBM) và mơ hình độ nhớt xốy cục bộ thích ứng với tường (WALE). Kết quả thu được cho thấy vận tốc trung bình của dịng khơng khí dưới 1,6 m đủ cao để làm xáo trộn hầu hết tán cây và tăng cường tốc độ xâm nhập của thuốc trừ sâu vào cây trồng. Ngoài ra, trường dịng chuyển động xi bị ảnh hưởng đáng kể bởi tốc độ bay vì tốc độ máy bay càng tăng thì tốc độ xốy càng tăng. Một mơ phỏng khác về máy bay nông nghiệp cũng sử dụng phương pháp Lattice Boltzmann, liên quan đến việc điều khiển tốc độ quay của rotor nhằm điều chỉnh chế độ vận hành phun chính xác [8]. Nghiên cứu được thực hiện bởi nhóm tác giả của Kun Chang, sử dụng mơ hình máy bay 4 chong chóng 410S và phần mềm Xflow để mơ

</div><span class="text_page_counter">Trang 23</span><div class="page_container" data-page="23">

phỏng trường dịng vận tốc rotor ở các tốc độ quay khác nhau; đồng thời sử dụng máy đo vận tốc hình ảnh hạt tốc độ cao (PIV) để thu thập và xác minh trạng thái chuyển động của giọt phun. Kết quả mô phỏng số cho thấy vận tốc và độ xoáy cực đại của trường dịng dưới chong chóng tăng khi tốc độ quay chong chóng tăng. Tuy nhiên, khi thời gian tăng lên, nhiễu loạn được tạo ra và các giá trị tối đa của vận tốc cùng độ xoáy giảm xuống. Trong khi đó, kết quả thử nghiệm PIV cho thấy khi tốc độ chong chóng tăng lên, tổng góc phun và diện tích vùng tốc độ cao của khu vực phun bị co lại và giảm dần dưới tác động của vùng gió xi dịng, cịn vùng tốc độ thấp ở diện tích phun dần dần mở rộng. Bên cạnh đó, thí nghiệm cũng chỉ ra rằng kết quả thử nghiệm PIV xác minh một cách hiệu quả độ tin cậy của kết quả mô phỏng số [9].

<b>1.4. Mục tiêu nghiên cứu </b>

Luận văn sử dụng 1 mô hình máy bay khơng người lái 4 chong chóng với bình thuốc trừ sâu tích hợp vào thân máy bay, đã được thiết kế và hiệu chỉnh bởi 1 nhóm nghiên cứu [10]. Mục tiêu của luận văn là sử dụng phương pháp mơ phỏng số CFD để phân tích các đặc tính khí động của mơ hình drone này trong các chế độ bay: chế độ bay treo (hovering), chế độ bay lấy độ cao (climbing), và chế độ bay tiến (forward flight) thơng qua việc phân tích lực và phân bố áp suất, vận tốc với các tốc độ quay của chong chóng khác nhau.

Tuy nhiên, đề tài chỉ giới hạn ở 1 số trường hợp mô phỏng cơ bản (không xem xét được hết các trường hợp bay ở điều kiện thực) như: mô phỏng bay treo, bay lấy độ cao không gia tốc, bay tiến với sự thay đổi góc nghiêng, nhằm mục đích bước đầu đánh giá về sự phù hợp của thiết kế khí động trong các chế độ bay.

<b>1.5. Phương pháp nghiên cứu </b>

<b>1.5.1. Phương pháp nghiên cứu lý thuyết </b>

- Nghiên cứu lý thuyết về khí động lực học, lý thuyết về động lực học lưu chất (CFD) và nguyên lý hoạt đông của máy bay không người lái;

</div><span class="text_page_counter">Trang 24</span><div class="page_container" data-page="24">

- Tham khảo các báo cáo, bài báo nghiên cứu khoa học về máy bay không người lái nơng nghiệp được sử dụng với mục đích phun thuốc trừ sâu đã được công bố trong các hội nghị khoa học, các tạp chí khoa học trong và ngồi nước.

<b>1.5.2. Phương pháp mơ phỏng </b>

Luận văn sử dụng phương pháp mô phỏng số Động lực học lưu chất (Computational Fluid Dynamics – CFD) với phương pháp Nhiều khung tham chiếu (Multiple Reference Frames – MRF), là một trong những phương pháp CFD phổ biến để mơ hình hóa dịng lưu chất trong các hệ thống chuyển động quay. Phần mềm được sử dụng là Ansys Fluent phiên bản 2022 R1, một gói phần mềm thương mại dựa trên Phương pháp Phần tử Hữu hạn (Finite Element Method - FEM) để phân tích các bài tốn vật lý - cơ học, chuyển các phương trình vi phân và phương trình đạo hàm riêng từ dạng giải tích về dạng số.

<b>1.6. Đối tượng và phạm vi nghiên cứu </b>

- Đối tượng nghiên cứu chính của luận văn là máy bay không người lái 4 chong chóng (quadcopter) có tích hợp bình thuốc trừ sâu vào thân.

- Phạm vi nghiên cứu: Luận văn tập trung nghiên cứu, đưa ra các số liệu để đánh giá về đặc tính khí động của mơ hình máy bay đã thiết kế (drone có tích hợp bình thuốc trừ sâu vào thân), từ đó chỉ ra sự phù hợp của mơ hình này (về mặt khí động) trong các chế độ bay.

<b>1.7. Cấu trúc luận văn </b>

- Chương 2: Cơ sở lý thuyết

Chương 2 trình bày khái quát về phương pháp mô phỏng số CFD sử dụng trong luận văn, các loại lưới sử dụng trong mô phỏng, các mơ hình rối, lý thuyết xử lý gần tường và các giải thuật sử dụng trong phần mềm.

- Chương 3: Thiết lập bài tốn mơ phỏng 1 chong chóng đơn

Chương 3 trình bày các bước mơ phỏng đối với 1 chong chóng đơn, đồng thời so sánh kết quả mô phỏng với thực nghiệm.

- Chương 4: Thiết lập bài tốn mơ phỏng tồn bộ mơ hình drone

</div><span class="text_page_counter">Trang 25</span><div class="page_container" data-page="25">

Từ các thơng số đầu vào lựa chọn đối với bài tốn 1 chong chóng đơn, Chương 4 tiến hành các bước thiết lập đối với bài tốn mơ phỏng tồn bộ mơ hình drone, từ bước xây dựng miền hình học, chia lưới cho đến thiết lập các thông số mô phỏng ở 3 chế độ bay: bay treo, bay lấy độ cao và bay tiến.

- Chương 5: Phân tích khí động học của drone khi hoạt động

Từ kết quả mô phỏng thu được, nội dung Chương 5 sẽ phân tích cụ thể ở 3 khía cạnh: Phân tích lực, Trường áp suất, và Trường vận tốc ở 3 chế độ bay đã đề cập.

- Chương 6: Kết luận và hướng phát triển

Chương 6 cũng là chương cuối cùng của bài luận văn, trình bày tóm tắt kết quả đạt được, hạn chế và hướng phát triển của đề tài.

</div><span class="text_page_counter">Trang 26</span><div class="page_container" data-page="26">

<b>Chương 2. CƠ SỞ LÝ THUYẾT 2.1. Giới thiệu máy bay UAV (Drone) </b>

Dựa theo số chong chóng, có thể phân ra nhiều loại máy bay không người lái (từ 2 đến 8 chong chóng). Nhiều chong chóng hơn giúp cải thiện độ ổn định của máy bay và khả năng chịu tải, nhưng những máy bay không người lái như vậy cần nhiều năng lượng pin hơn để điều khiển nhiều động cơ nhằm đạt được công suất cao. Trên thực tế, người ta thường sử dụng phổ biến là máy bay không người lái 4 chong chóng (gọi là quadcopter). Quadcopter có bốn chong chóng ở bốn góc của khung và được đặt cách nhau một khoảng bằng nhau. Đối với mỗi chong chóng, tốc độ và hướng quay được điều khiển độc lập để cân bằng và tạo ra chuyển động của máy bay. Tuy nhiên, để duy trì sự cân bằng của hệ thống máy bay, một cặp chong chóng quay theo chiều kim đồng hồ và cặp còn lại quay theo hướng ngược chiều kim đồng hồ (hình 2.1a) [11].

<i>Hình 2.1 Cơ chế hoạt động của máy bay UAV 4 chong chóng: </i>

<i>(a) Chiều quay chong chóng; (b) Di chuyển của drone từ bay treo sang bay tiến [12] </i>

Khi drone bay ở chế độ bay treo, lực đẩy sẽ cùng hướng với lực nâng; tuy nhiên, khi bay tiến (hoặc lùi), máy bay sẽ nghiêng 1 góc so với phương ngang, khi đó lực đẩy tạo ra trên drone sẽ thay đổi theo hướng chuyển động của máy bay (hình 2.2).

</div><span class="text_page_counter">Trang 27</span><div class="page_container" data-page="27">

<i>Hình 2.2 Các lực tạo ra trên drone: (a) Bay treo; (b) Bay tiến [11] </i>

<b>2.2. Phương pháp tính tốn mơ phỏng số động lực học lưu chất (Computational Fluid Dynamics - CFD) </b>

Phương pháp tính tốn mô phỏng số động lực học lưu chất (CFD) là một ngành khoa học và là một nhánh của cơ lưu chất, sử dụng các máy tính kỹ thuật số để đưa ra các dự đốn định tính và định lượng về các hiện tượng dòng lưu chất dựa trên các định luật bảo toàn (bảo toàn khối lượng, động lượng và năng lượng) chi phối chuyển động của lưu chất [13]. Máy tính được dùng để thực hiện các phép tính cần thiết nhằm mơ phỏng dịng chuyển động tự do của lưu chất và sự tương tác của lưu chất (chất lỏng và chất khí) với các bề mặt được xác định bởi các điều kiện biên. Với các siêu máy tính tốc độ cao, việc ứng dụng CFD có thể giúp mang lại các giải pháp tốt hơn và thường được yêu cầu để giải các bài tốn khó và phức tạp [14].

CFD được áp dụng cho nhiều vấn đề nghiên cứu và kỹ thuật trong nhiều lĩnh vực nghiên cứu và công nghiệp, bao gồm phân tích khí động học và hàng khơng vũ trụ, siêu âm, mô phỏng thời tiết, khoa học tự nhiên và kỹ thuật mơi trường, thiết kế và phân tích hệ thống cơng nghiệp, kỹ thuật sinh học, dịng chất lỏng và nhiệt chuyển giao, phân tích động cơ và đốt cháy, và hiệu ứng hình ảnh cho phim và trị chơi [14].

</div><span class="text_page_counter">Trang 28</span><div class="page_container" data-page="28">

<b>2.2.2. Các phương trình điều khiển của CFD </b>

Tất cả các phương pháp CFD đều dựa trên các phương trình điều khiển cơ bản của động lực học lưu chất, đó là phương trình liên tục, phương trình động lượng và phương trình năng lượng – được phát biểu bởi 3 định luật vật lý lần lượt là: Định luật bảo toàn khối lượng, Định luật II Newton và Định luật bảo tồn năng lượng.

Các phương trình điều khiển có thể thu được ở các định dạng khác nhau. Có rất ít sự khác biệt giữa các phương trình này đối với lý thuyết khí động học, nhưng đối với một số thuật toán nhất định trong CFD, việc lựa chọn các phương trình điều khiển là rất quan trọng do việc sử dụng các phương trình khác nhau dẫn đến các kết quả khác nhau. Trong số này có thể là kết quả chính xác, có thể là kết quả khơng chính xác và thậm chí là khơng ổn định hoặc khác nhau. Do đó, cần cẩn thận để đảm bảo sử dụng phương trình điều khiển phù hợp để giải quyết vấn đề thực tế trong quá trình CFD. Các đại lượng vật

<i>lý được sử dụng trong các phương trình này bao gồm áp suất p; khối lượng riêng </i><small></small> ;

<i>nhiệt độ T; vận tốc theo 3 phương u , v , w; … </i>

- Phương trình liên tục (Phương trình bảo toàn khối lượng) [16]:

</div><span class="text_page_counter">Trang 30</span><div class="page_container" data-page="30">

<b>2.2.3. Các phương pháp rời rạc hóa [14] </b>

<i>2.2.3.1. Phương pháp thể tích hữu hạn (Finite Volume Method – FVM) </i>

Phương pháp thể tích hữu hạn (FVM) là một kỹ thuật rời rạc hóa cho các phương trình đạo hàm từng phần, đặc biệt là những phương trình phát sinh từ các định luật bảo toàn vật lý. FVM sử dụng cơng thức tích phân khối của bài tốn với một tập hợp phân vùng hữu hạn của khối lượng để rời rạc hóa các phương trình, đặc biệt là các phương trình động lực học lưu chất.

Phương pháp thể tích hữu hạn có lợi trong việc sử dụng bộ nhớ và tốc độ giải, đặc biệt đối với các vấn đề lớn, dòng hỗn loạn với số Reynolds cao và dòng bị chi phối bởi giới hạn nguồn. Nó được trình bày dưới dạng phương trình điều khiển như sau:







<i><b>trong đó Q là vectơ của các biến được bảo toàn, F là vectơ của dịng lưu chất (xem </b></i>

phương trình Euler hoặc phương trình Navier–Stokes), V là thể tích của phần tử thể tích kiểm sốt và A là diện tích bề mặt của phần tử thể tích kiểm sốt.

<i>2.2.3.2. Phương pháp phần tử hữu hạn (Finite Element Method - FEM) </i>

Phương pháp phần tử hữu hạn (FEM) là một phương pháp số để giải phương trình vi phân hoặc tích phân, được áp dụng cho một số bài toán vật lý, trong đó có sẵn các phương trình vi phân điều khiển. Về cơ bản, phương pháp này bao gồm việc giả sử hàm liên tục từng phần để giải và thu được các tham số của các hàm theo cách làm giảm lỗi trong khi giải. Phương pháp này ổn định hơn nhiều so với cách tiếp cận thể tích hữu hạn. Tuy nhiên, FEM có thể u cầu nhiều bộ nhớ hơn và có thời gian giải quyết chậm hơn so với FVM.

Trong phương pháp này, một phương trình phần dư có trọng số (weighted residual equation) được hình thành:

</div><span class="text_page_counter">Trang 31</span><div class="page_container" data-page="31">

trong đó <i>R<sub>i</sub> là phương trình phần dư tại một đỉnh i của phần tử, </i>

<i>Q</i>

<i><sub>i</sub></i>là phương trình bảo tồn được biểu thị trên cơ sở phần tử,

<i>W</i>

<i><sub>i</sub></i> là hệ số trọng lượng và <i>V<sup>e</sup></i>là thể tích của phần tử.

<i>2.2.3.3. Phương pháp sai phân hữu hạn (Finite Difference Method - FDM) </i>

Một trong các phương pháp được biết đến nhiều nhất mà tương đối đơn giản và ổn định là phương pháp sai phân hữu hạn, bao gồm việc thay thế từng đạo hàm bằng một thương số khác nhau tương ứng trong cơng thức cổ điển. Nói cách khác, công thức sai phân hữu hạn đưa ra cách tiếp cận trực tiếp hơn đối với nghiệm số của phương trình đạo hàm riêng so với phương pháp dựa trên các công thức khác. Hạn chế của phương pháp sai phân hữu hạn là độ chính xác và tính linh hoạt. Các phương pháp sai phân hữu hạn tiêu chuẩn địi hỏi tính quy tắc hơn của nghiệm

(

( ) )

<i>u</i><i>C</i>  và tam giác (lưới đồng nhất). Khó khăn cũng nảy sinh trong việc áp đặt các điều kiện biên.

0 <sub>+</sub> <sub>+</sub> <sub>+</sub> <sub>=</sub>

<i>t</i> <i>x</i> <i>y</i> <i>z</i>

<i><b>trong đó Q là vectơ của các biến được bảo toàn và F, G và H lần lượt là vector của các </b></i>

dòng lưu chất theo các hướng x, y và z.

<b>2.2.4. Quy trình mơ phỏng CFD </b>

Q trình mơ phỏng CFD là q trình giải quyết các bài tốn về dịng lưu chất và q trình này được cấu trúc quanh các thuật tốn số để phục vụ mục đích trên. Q trình tính toán động lực học lưu chất đều bao gồm ba giai đoạn là: tiền xử lý (pre – processing), xử lý (processing) và hậu xử lý (post – processing).

<i><b>- Tiền xử lý (Pre – processing): sẽ bao gồm các bước chính như sau: </b></i>

+ Xác định vấn đề: xác định cụ thể mục tiêu phân tích và kết quả dự báo đạt được, các yếu tố tác động đến bài toán như khả năng trao đổi năng lượng, loại dòng chảy của lưu chất (tầng hay rối), điều kiện hoạt động, các điều kiện biên…

+ Xây dựng và xử lý hình học: là bước xây dựng mơ hình vật thể và mơ hình vùng lưu chất bao quanh vật thể cần mơ phỏng. Hai mơ hình này có vai trị quan trọng trong

</div><span class="text_page_counter">Trang 32</span><div class="page_container" data-page="32">

việc diễn tả đúng bản chất bài toán, tối ưu hóa q trình mơ phỏng và là nền tảng của q trình thiết lập mơ hình lưới tính tốn.

+ Xác định các điều kiện biên: là việc thiết lập các điều kiện biên vật lý trên các bề mặt của vùng lưu chất hữu hạn, nhằm xác định các yếu tố đầu vào của mơ hình mơ phỏng và kết nối mơ hình với mơi trường xung quanh.

+ Xây dựng mơ hình lưới (chia lưới): là q trình chia nhỏ mơ hình vật thể và vùng lưu chất thành các ô lưới nhỏ kết nối với nhau, khơng chồng chéo lên nhau. Đây là q trình tiêu tốn nhiều thời gian nhất và có vai trị cực kỳ quan trọng trong mô phỏng CFD. Đồng thời cũng góp phần quyết định độ chính xác của kết quả mơ phỏng, thời gian tính tốn, tài ngun máy tính.

<i><b>- Xử lý (Processing): là việc thiết lập quá trình mơ phỏng bao gồm thiết lập tính chất </b></i>

dịng chuyển động lưu chất về mặt không gian (2D, 3D), thời gian (ổn định, khơng ổn định), mơ hình rối, giải thuật tính tốn, …

<i><b>- Hậu xử lý (Post – processing): sau khi hồn thành q trình mơ phỏng, ta sẽ nhận </b></i>

được các kết quả dưới dạng hình ảnh, đồ thị hoặc số liệu. Các kết quả sẽ được phân tích, so sánh nhằm đánh giá ảnh hưởng, tính chất và mức độ thành cơng của bài tốn.

<i>Hình 2.4 Ví dụ về quy trình các bước mơ phỏng trong Ansys Fluent </i>

<b>2.3. Lưới trong CFD </b>

<b>2.3.1. Định nghĩa lưới </b>

Trong một bài tốn CFD, để tính tốn kết quả chính xác thì mơ hình sẽ được chia thành nhiều ô lưới với mỗi ô là các phần tử dòng chuyển động được giải quyết. Lưới

</div><span class="text_page_counter">Trang 33</span><div class="page_container" data-page="33">

chính là đại diện rời rạc về hình học của vấn đề cần giải quyết. Ơ lưới có tác động đáng kể đến tỷ lệ hội tụ, độ chính xác của giải thuật và thời gian tính tốn của CPU máy tính. Nếu chia lưới q thưa thì kết quả tính được sẽ khơng chính xác, cịn nếu ô lưới nhỏ quá thì sẽ tiêu hao nhiều tài nguyên, mất nhiều thời gian.Vì vậy, ta cần phải chia lưới một cách hiệu quả và chính xác. Các thành phần cơ bản của lưới (2D và 3D) bao gồm [17]:

• Ơ lưới hoặc phần tử lưới (Cell/Element); • Nút lưới (Node);

• Trung tâm một ơ lưới (Cell center); • Mặt (Face);

• Cạnh (Edge).

<i>Hình 2.5 Thành phần cơ bản của lưới [17] </i>

Các loại hình dạng của phần tử lưới 2D và 3D phổ biến được minh họa như hình 2.6:

</div><span class="text_page_counter">Trang 34</span><div class="page_container" data-page="34">

- Lưới phi cấu trúc (unstructured mesh): là các mắt lưới có kết nối chung và có cấu trúc là tự do và do đó kết nối của các phần tử phải được xác định và lưu trữ. Vì tính tự do này mà lưới có thể dễ dàng được tinh chỉnh mức độ phân giải của hệ lưới sao cho phù hợp với hình học mà cần tính tốn. Các loại phần tử kết nối chung là khơng trực giao, chẳng hạn như tam giác (2D) và tứ diện (3D).

- Lưới lai hay còn gọi là lưới hỗn hợp (hybrid mesh): là sự kết hợp hai loại lưới cấu trúc và không cấu trúc vô một hình học để tận dụng ưu điểm của hai loại và làm cho việc tính tốn mơ hình được tối ưu hơn.

</div><span class="text_page_counter">Trang 35</span><div class="page_container" data-page="35">

Để đảm bảo cho một bài tốn mơ phỏng đạt hội tụ thì việc chia lưới và đạt chuẩn các tiêu chuẩn về lưới là điều khơng thể thiếu. Việc chia lưới địi hỏi sự tỉ mỉ về từng chất lượng, kích thước gữa các ơ lưới phải phù hợp bài tốn mình đang xét. Sau đây là những phương pháp đánh giá chất lượng lưới được được sử dụng trong mơ phỏng bài tốn này:

- Tiêu chuẩn đánh giá chất lượng từng phần tử lưới: Skewness (độ cân xứng các phần tử lưới) và Aspect ratio (tỉ lệ kích thước phần tử lưới).

- Tiêu chuẩn đánh giá tổng thể chất lượng lưới: Orthogonal quality (độ không trực giao) và Smoothness (độ mịn lưới).

<b>a) Độ cân xứng của lưới (skewness) [17] </b>

Công thức của độ cân xứng lưới dựa vào độ lệch góc đều:

</div><span class="text_page_counter">Trang 36</span><div class="page_container" data-page="36">



<i><sub>e</sub></i>: là góc mà phần tử ơ lưới đang xét (đặc biệt nếu chia lưới tam giác là 60 độ, chia lưới tứ giác là 90<small>0</small>)

Phạm vi xét độ cân xứng của lưới là từ 0 đến 1. Độ cân xứng tốt nhất là bằng 0 và không đạt khi bằng 1. Tiêu chuẩn này được xem là tiêu chuẩn quan trọng nhất trong ba tiêu chuẩn, nếu nó tốt thì các tiêu chuẩn khác cũng tốt.

<i>Bảng 2.1 Chất lượng phần tử lưới dựa trên các khoảng giá trị độ cân xứng lưới [20] </i>

<b>Xuất sắc Rất tốt Tốt Trung bình Xấu Rất xấu </b>

0 – 0,25 0,25 – 0,50 0,50 – 0,80 0,80 – 0,94 0,95 – 0,97 0,98 – 1,00

<b>b) Tỉ lệ kích thước các cạnh của phần tử (AR) [21] </b>

Định nghĩa: là tỉ lệ độ dài giữa cạnh lớn nhất và cạnh nhỏ nhất của một ô lưới. Công thức dùng đánh giá cho phương pháp này:

- Đối với phần tử lưới dạng tam giác và tứ diện:

<i>f =</i> đánh giá 2D, <sup>1</sup>3

<i>f =</i> đánh giá 3D

- Đối với phần tử lưới dạng hình chóp và lục giác:

(

)

max , ,...,min , ,...,

<i>• e<small>i</small></i> là trung bình cộng 2 cạnh đối của 1 phần tử lưới

• n là tổng số hướng (n = 2 với dạng hình chóp; n = 3 với dạng lục giác) Phạm vi xét tiêu chuẩn này với giá trị lý tưởng là 1, tuy nhiên khó có thể đạt được.

</div><span class="text_page_counter">Trang 37</span><div class="page_container" data-page="37">

<i>Hình 2.9 Chất lượng phần tử lưới dựa trên tiêu chuẩn AR </i>

<b>c) Độ không trực giao (orthogonal quality) [20] </b>

Định nghĩa: tiêu chí này đo góc giữa 2 tâm của ơ lưới và pháp tuyến của đường giao tuyến của hai ô lưới đó.

- Đối với phần tử lưới, độ không trực giao được xét theo công thức sau:

• <i><small>A</small><sub>i</sub></i> là vectơ pháp tuyến của mặt

• <i><small>f</small><sub>i</sub></i> là vectơ nối từ tâm phần tử đến tâm của mặt đó • <i>c là vectơ nối tâm 2 phần tử liên kề <sub>i</sub></i>

• <i>e là vectơ nối từ trọng tâm của mặt đến trọng tâm của cạnh <sub>i</sub></i>

</div><span class="text_page_counter">Trang 38</span><div class="page_container" data-page="38">

<i>Hình 2.10 Ví dụ về các vectơ trong công thức đánh giá độ không trực giao [20] </i>

Phạm vi xét tiêu chuẩn này cũng chạy từ 0 đến 1. Tiêu chuẩn tốt nhất khi đạt ở gần 1, và tệ nhất là gần 0.

<i>Bảng 2.2 Chất lượng phần tử lưới dựa trên các khoảng giá trị độ không trực giao [20] </i>

<b>Rất xấu Xấu Trung bình Tốt Rất tốt Xuất sắc </b>

0 – 0,001 0,001 – 0,14 0,15 – 0,20 0,20 – 0,69 0,70 – 0,95 0,95 – 1,00

<b>d) Độ mịn lưới (smoothness) </b>

Tiêu chuẩn này xét đến sự thay đổi kích thước của các phần tử liền kề nhau. Lưới được xem là đảm bảo tiêu chuẩn về độ mịn khi sự thay đổi kích thước của 2 phần tử lưới liên tiếp không được vượt quá 20%.

</div><span class="text_page_counter">Trang 39</span><div class="page_container" data-page="39">

<b>2.4. Mô phỏng dịng rối </b>

<b>2.4.1. Phương pháp mơ phỏng dịng rối </b>

Dòng chuyển động được xem là rối khi các phần tử chất lỏng hay chất khí chuyển động khơng có trật tự và hỗn loạn. Ngược lại, dòng chuyển động được xem là tầng khi các phần tử di chuyển thành từng lớp và trật tự. Người ta thường dùng số Reynolds của dịng để dự đốn một dịng là tầng hay rối. Số Reynolds của dòng là tỉ số giữa lực quán tính (inertia force) của một phần tử chất lỏng hay chất khí và lực ma sát nhớt (viscous force) tác dụng lên phần tử đó [22]. Cơng thức tính số Reynolds như sau:

<i>Số Reynolds (Re) = Lực qn tính/Lực nhớt(2-15)</i>

Từ số Reynolds của dịng, ta đem so sánh với số Reynolds tới hạn – Critical Reynolds

<i>number (Re<small>crit</small>) để xác định trạng loại dòng. Nếu Re < Re<small>crit</small> thì dịng này được gọi là dịng tầng. Nếu Re > Re<small>crit</small> thì dịng này được gọi là dịng rối. Nếu Re ≈ Re<small>crit</small> thì được </i>

<i>gọi là dòng chuyển tiếp, đây là dòng hỗn hợp của dòng tầng và dòng rối, với dòng rối ở </i>

trung tâm và dòng tầng ở gần biên. Phần lớn các vấn đề về lưu chất gặp phải khi mô phỏng đều là dòng rồi, và với mỗi loại dòng rối thì ta cần phải có các cách tiếp cận khác nhau. Hiện nay, có ba phương pháp mơ phỏng dịng rối phố biến: Mô phỏng trực tiếp – Direct Numerical Simulation (DNS); Mô phỏng Navier – Stokes sử dụng số Reynolds trung bình – Reynolds Averaged Navier – Stokes (RANS); Mơ phỏng xốy lớn – Large Eddy Simulation (LES).

<b>- Phương pháp DNS: là phương pháp mô phỏng dòng 3D và phụ thuộc vào thời gian. </b>

Đặc biệt, với phương pháp này có thể mơ phỏng được tất cả xốy lớn, xốy nhỏ trong dịng được tính tốn. Do mơ phỏng được đầy đủ dịng rối nên phương pháp này cung cấp được những giá trị khó lấy được từ thực nghiệm, ngồi ra có thể dùng để kiểm tra với các mơ hình rối khác. Tuy nhiên, phương pháp này có một số hạn chế như: chi phí tính tốn lớn, u cầu bộ nhớ lớn, cung cấp số lượng thông tin quá lớn, không khả thi trong công nghiệp, …

</div><span class="text_page_counter">Trang 40</span><div class="page_container" data-page="40">

<b>- Phương pháp RANS: là phương pháp sử dụng phương pháp lấy trung bình, do đó, </b>

khơng tính tốn cụ thể mọi xốy lớn, nhỏ trong dịng mà chỉ mơ phỏng và đưa ra giá trị trung bình, vì thế tiết kiệm được thời gian và khối lượng tính tốn. Hiện nay, đây là phương pháp được sử dụng nhiều nhất trong công nghiệp và trong các nghiên cứu.

<b>- Phương pháp LES: là phương pháp lai giữa DNS và RANS. Cách tiếp cận của phương </b>

pháp này là giải trực tiếp các xốy lớn và mơ phỏng các xốy nhỏ. Nhờ đó, độ chính xác và đầy đủ của kết quả tính tốn là vượt trội hơn hẳn so với phương pháp RANS và gần tương đương với phương pháp DNS, trong khi vẫn đảm bảo được khối lượng và chi phí tính tốn khơng q cao.

<i>Hình 2.12 Các phương pháp mơ phỏng dịng rối phổ biến [23] </i>

Trong phương pháp mơ phỏng RANS, ta phân loại các mơ hình rối dựa vào số lượng phương trình vận chuyển cộng thêm vào mơ hình:

<i>Bảng 2.3 Bảng phân loại mơ hình rối theo số phương trình vận chuyển </i>

<b>Loại mơ hình rối <sup>Số lượng phương trình </sup>vận chuyển thêm </b>

3) k-epsilon (standard, RNG, realizable) 2

</div>