Tổng quan luận án
Đề tài luận án tiến sĩ “Ứng dụng phương pháp tính toán động lực học lưu chất (CFD) trong tối ưu hóa hình dạng mũi tàu quả lê” được thực hiện bởi nghiên cứu sinh Huỳnh Văn Chính, dưới sự hướng dẫn khoa học của PGS. Trần Gia Thái và TS. Bùi Hồng Dương tại Trường Đại học Giao thông vận tải Thành phố Hồ Chí Minh (bảo vệ năm 2022). Luận án thuộc chuyên ngành Kỹ thuật tàu thủy / Cơ khí động lực, đi sâu vào lĩnh vực thủy động lực học tàu thủy và kỹ thuật tối ưu hóa hình dáng hình học vỏ tàu.
Sử dụng mũi quả lê là giải pháp thủy động học nhằm giảm sức cản tổng của tàu (có thể đạt mức giảm từ 10% đến 15% trong điều kiện thuận lợi) nhờ tạo ra hệ thống sóng giao thoa ngược pha tích cực triệt tiêu hệ thống sóng mũi của thân tàu. Đồng thời, mũi quả lê cải thiện các tính năng hàng hải như độ chúi dọc và tính ổn định lắc dọc. Tuy nhiên, việc thiết kế hình dáng và vị trí quả lê gặp nhiều khó khăn do tính phức tạp của tương tác sóng; nếu hình dạng không phù hợp, sự giao thoa tiêu cực sẽ làm tăng sức cản tổng. Các phương pháp thực nghiệm truyền thống (kéo mô hình trong bể thử) đòi hỏi nhiều thời gian và chi phí.
Tại Việt Nam, sau khi triển khai Nghị định 67/2014/NĐ-CP, hàng loạt tàu cá vỏ thép dưới 30 m đã được đóng mới phục vụ khai thác xa bờ nhưng hầu hết sử dụng dạng mũi thẳng truyền thống và không qua thử nghiệm bể thử, dẫn đến hiệu quả kinh tế - kỹ thuật chưa tối ưu. Trong bối cảnh nhà nước định hướng hiện đại hóa đội tàu cá vỏ thép cỡ lớn (chiều dài trên 40 m), việc xây dựng cơ sở khoa học để thiết kế và tối ưu hóa mũi quả lê cho tàu cá trở nên cấp thiết. Khoảng trống nghiên cứu được xác định gồm:
- Phương pháp đồ thị thiết kế kinh điển của Kracht (1978) bị giới hạn trong dải hệ số béo thể tích $C_B \in [0.56, 0.82]$, không áp dụng trực tiếp được cho tàu cá có hệ số béo nhỏ ($C_B < 0.56$), đồng thời chưa giải quyết bài toán ghép nối trơn đều bề mặt quả lê vào thân tàu.
- Các nghiên cứu tối ưu hóa mũi quả lê bằng CFD hiện đại chủ yếu áp dụng cho tàu hàng chạy ở một chế độ vận tốc/mớn nước cố định với hàm đơn mục tiêu, chưa xây dựng được mô hình tối ưu đa mục tiêu kèm ràng buộc hàng hải phù hợp với đặc thù khai thác đa chế độ của tàu cá.
- Thiếu các nghiên cứu chuẩn hóa thông số đầu vào cho mô phỏng số CFD (kích thước miền tính, điều kiện biên, hệ số mô hình rối) để đảm bảo độ chính xác kiểm chứng cho tàu cá vỏ thép.
Mục tiêu chung của luận án là ứng dụng phương pháp CFD trong tối ưu hóa hình dạng mũi quả lê nhằm đạt độ giảm sức cản tổng lớn nhất, thỏa mãn các ràng buộc về hình học và tính năng hàng hải.
Các mục tiêu cụ thể được đánh số theo văn bản:
- Phân tích các tài liệu và công trình nghiên cứu có liên quan đến phương pháp CFD và ứng dụng trong tính sức cản tàu và tối ưu hóa hình dạng mũi tàu quả lê nhằm xây dựng hướng nghiên cứu và các cơ sở phương pháp luận cần thiết.
- Ứng dụng CFD tính sức cản của tàu nghiên cứu với độ chính xác mong đợi để làm cơ sở giải quyết bài toán tối ưu hóa quả lê theo hàm mục tiêu sức cản.
- Xây dựng phương pháp tối ưu hình dạng mũi quả lê phù hợp với tàu nghiên cứu.
Đối tượng và phạm vi nghiên cứu:
- Đối tượng nghiên cứu: Mẫu tàu đánh cá vỏ thép cỡ vừa và nhỏ của Tổ chức Nông Lương Liên Hiệp Quốc (FAO), cụ thể là hai mẫu tàu FAO 72 và FAO 75 đã được thử nghiệm sức cản tại các bể thử National Physical Laboratory (NPL - Anh) và Swedish State Shipbuilding Tank (SSPA - Thụy Điển).
- Phạm vi nghiên cứu:
- Chế độ bơi theo định luật Archimede, tàu chuyển động thẳng đều trong môi trường nước tĩnh, không chịu tác động của sóng biển tự do và sức cản không khí, độ sâu nước vô hạn.
- Lưu chất trong mô phỏng là chất lỏng đồng chất, có tính nhớt và không nén được.
- Dữ liệu thực nghiệm kéo mô hình trong bể thử của FAO được sử dụng làm chuẩn đối chứng để hiệu chỉnh và đánh giá sai số của mô hình CFD.
Tổng quan tài liệu và vị trí của luận án
Luận án đã hệ thống hóa lịch sử nghiên cứu và phát triển mũi quả lê qua hai giai đoạn chính:
Các công trình nghiên cứu truyền thống
- D. Taylor (1910) lần đầu đưa mũi quả lê vào thiết kế tàu chiến USS Delaware của Hải quân Mỹ.
- Tàu khách Bremen và Europa (Đức, 1920), tàu Hoover và Coolidge (Mỹ, 1931), tàu Normandie (Pháp, 1935), tàu chiến Yamato (Nhật Bản, 1940) và tàu chở hàng Yamashiro Maru (Nhật Bản, 1963) là các mốc ứng dụng thực tế cho thấy hiệu quả giảm công suất máy chính từ 10% đến trên 13%.
- Các dạng mũi quả lê hình học ban đầu gồm: dạng ống trụ tròn (kiểu C của Đại học Michigan) và dạng bóng đèn tròn (kiểu 3A, 3B, 3C).
- Giáo sư T. Inui (1962) đưa ra lý thuyết "dạng tàu không sóng" (waveless form) dựa trên phân bố các nguồn âm dương (điểm kỳ dị toán học) để triệt tiêu sóng mũi và sóng lái, thử nghiệm thành công trên tàu Murasaki và Kurenai Maru (giảm 13.1% công suất máy chính).
- Công trình thực nghiệm toàn diện của M. Kracht (1978): Thử nghiệm hàng loạt mô hình tàu có gắn mũi quả lê trong dải hệ số béo $C_B \in [0.56, 0.82]$ và thiết lập hệ thống đồ thị thiết kế Kracht dựa trên 6 tham số hình học không thứ nguyên ($C_{LPR}, C_{BB}, C_{ZB}, C_{ABT}, C_{ABL}, C_{PR}$).
Các công trình nghiên cứu hiện đại
Các nghiên cứu hiện đại ứng dụng CFD kết hợp thuật toán tối ưu hóa được chia làm hai xu hướng:
- Tối ưu hóa dựa trên phương pháp phương án: Grzegorz Filip và cộng sự (Đại học Michigan & ABS, 2014) tối ưu hóa mũi quả lê tàu container KCS (KRISO Container Ship) bằng cách thay đổi chiều cao $H$ và chiều rộng $C$ qua 9 phương án, giải bài toán RANSE với mô hình rối SST $k-\omega$ trên OpenFOAM và sử dụng mô hình thay thế Kriging để xây dựng bề mặt đáp ứng tìm điểm tối ưu.
- Tối ưu hóa dựa trên phương pháp tham số:
- Weilin Luo và Linqiang Lan (Đại học Phúc Châu, 2016) sử dụng mô hình tham số B-Spline và mặt NURBS trong phần mềm CAESES kết hợp phần mềm tính CFD SHIPFLOW để tối ưu hóa quả lê tàu Ro-Ro ($L=165.7\text{ m}, C_B=0.65, Fn=0.28$).
- Nghiên cứu tối ưu hóa mũi quả lê tàu Seri 60 ($L_{pp}=7.37\text{ m}, C_B=0.6$) sử dụng thuật toán di truyền GA (Genetic Algorithm) lập trình trên MATLAB kết hợp CFD.
- Phương pháp hàm cơ sở xuyên tâm (Radial Basis Function - RBF) ứng dụng tối ưu hóa 4 mẫu mũi tàu thuộc tổ chức Joint High Speed Sealift (JHSS) gồm các biến thể Stem Bow (ST), Elliptical Bulb (EB), Gooseneck Bulb (GB) và Baseline Bulb (BB).
- Các nghiên cứu tại Việt Nam:
- Viện Cơ khí Động lực - Đại học Bách Khoa Hà Nội nghiên cứu phương pháp điểm kỳ dị tính sức cản sóng cho tàu dạng mảnh (như tàu Wigley).
- Đại học Hàng hải Việt Nam, Đại học Giao thông vận tải TP.HCM nghiên cứu mô phỏng trường dòng chảy quanh thân tàu, tương tác chân vịt - bánh lái.
- Đại học Nha Trang: Luận án tiến sĩ của Lê Văn Toàn ứng dụng OpenFOAM tính sức cản tàu cá vỏ gỗ; nhóm nghiên cứu ứng dụng phương pháp Lackenby biến đổi đường hình tàu theo hoành độ tâm nổi (LCB) tối ưu.
- PGS.TS Đỗ Quang Khải (Đại học Hàng hải Việt Nam / Đại học Liège) xây dựng các hàm đa thức phân đoạn biểu diễn đường nước, vùng mũi và mũi quả lê.
Khoảng trống nghiên cứu được luận án giải quyết:
- Mở rộng phạm vi ứng dụng của đồ thị Kracht cho tàu cá có hệ số béo $C_B = 0.524$ nằm ngoài dải chuẩn ($C_B < 0.56$).
- Xây dựng quy trình chuẩn hóa các thông số đầu vào cho CFD (kích thước miền tính, mô hình rối SST $k-\omega$) để đảm bảo sai số tính sức cản tàu cá dưới 5% so với thực nghiệm bể thử.
- Thiết lập mô hình toán học tối ưu hóa đa mục tiêu (theo nhiều chế độ khai thác: hành trình, dắt lưới, kéo/thả lưới) kết hợp các điều kiện ràng buộc hình học và tính năng hàng hải cho tàu cá.
Cơ sở lý thuyết và phương pháp nghiên cứu
Cơ sở lý thuyết
- Hệ phương trình thủy động lực học: Luận án sử dụng hệ phương trình trung bình Reynolds Navier-Stokes (RANSE) cho dòng lưu chất nhớt, không nén được, biểu diễn dưới dạng bảo toàn khối lượng (phương trình liên tục) và bảo toàn động lượng:
$$\frac{\partial (\rho U_i)}{\partial x_i} = 0$$
$$\frac{\partial (\rho U_i)}{\partial t} + \frac{\partial (\rho U_i U_j)}{\partial x_j} = -\frac{\partial p}{\partial x_i} + \frac{\partial}{\partial x_j}\left[\mu \left(\frac{\partial U_i}{\partial x_j} + \frac{\partial U_j}{\partial x_i}\right) - \rho \overline{u'_i u'_j}\right] + S_i$$
- Mô hình rối: Mô hình SST $k-\omega$ (Shear Stress Transport) hai phương trình của Menter, kết hợp ưu điểm của mô hình $k-\omega$ ở vùng sát biên vỏ tàu và mô hình $k-\epsilon$ ở vùng dòng tự do xa vỏ tàu.
- Lý thuyết phân tách sức cản: Sức cản tổng của tàu $R_T$ được phân tích theo Froude và khuyến nghị của ITTC thành:
- Sức cản ma sát $R_F$ và sức cản áp suất nhớt $R_{PV}$ (hợp thành sức cản nhớt $R_V$).
- Sức cản sinh sóng $R_{WM}$ và sức cản sóng vỡ $R_{WB}$ (hợp thành sức cản sóng $R_W$).
- Sức cản dư $R_R = R_T - R_F$.
- Lý thuyết hình học mũi quả lê: Dựa trên hệ thống 6 hệ số hình học của Kracht:
- Hệ số chiều dài: $C_{LPR} = L_{PR} / L_{pp}$
- Hệ số chiều rộng: $C_{BB} = B_B / B_{MS}$
- Hệ số chiều cao: $C_{ZB} = Z_B / T_{FP}$
- Hệ số mặt cắt ngang tại sườn mũi: $C_{ABT} = A_{BT} / A_{MS}$
- Hệ số diện tích hình chiếu cạnh: $C_{ABL} = A_{BL} / (L_{pp} \cdot T_{FP})$
- Hệ số thể tích quả lê: $C_{\nabla PR} = V_{PR} / \nabla$
- Mô hình bài toán tối ưu đa mục tiêu: Hàm mục tiêu là cực đại hóa độ giảm công suất có ích trung bình $\Delta P_e$ qua các chế độ hoạt động:
$$\Delta P_e = \sum_{i=1}^{n} w_i \cdot \Delta P_{ei} \rightarrow \max$$
Trong đó $w_i$ là trọng số chế độ hoạt động ($w_1 + w_2 + w_3 = 1$), $\Delta P_{ei} = \frac{P_{ei} - P_{ebi}}{P_{ei}} \times 100%$, với các điều kiện ràng buộc về biến động thể tích chiếm nước $\Delta \nabla$, hoành độ tâm nổi $\Delta LCB$, chiều cao tâm ổn định $\Delta MG$, và các giới hạn biên kích thước.
Phương pháp nghiên cứu
- Phương pháp mô hình hóa 3D: Sử dụng phần mềm chuyên ngành AutoShip để dựng mô hình vỏ tàu và mũi quả lê dưới dạng mặt NURBS và lưới đường cong B-Spline.
- Phương pháp mô phỏng số (CFD): Sử dụng phần mềm XFlow (và đối chiếu OpenFOAM) áp dụng phương pháp thể tích hữu hạn (FVM) giải hệ phương trình RANSE.
- Phương pháp tối ưu hóa: Sử dụng phương pháp bề mặt đáp ứng (Response Surface Method - RSM) và mô hình thay thế (Surrogate Model) để nội suy hàm mục tiêu từ ma trận các phương án tính toán CFD qua 2 bước lặp.
Nguồn tư liệu và dữ liệu
- Hồ sơ thiết kế đường hình và dữ liệu đo sức cản thực nghiệm trong bể thử NPL và SSPA của hai mẫu tàu cá vỏ thép FAO 72 và FAO 75.
| Thông số chính |
Ký hiệu |
Mẫu tàu FAO 72 |
Mẫu tàu FAO 75 |
| Chiều dài giữa hai đường vuông góc |
$L_{pp}$ (m) |
19.81 |
22.86 |
| Chiều rộng lớn nhất |
$B$ (m) |
6.40 |
6.70 |
| Chiều chìm thiết kế |
$d$ hay $T$ (m) |
2.44 |
2.74 (khảo sát đến 4.57) |
| Hệ số béo thể tích |
$C_B$ |
0.448 |
0.524 |
| Vận tốc thử nghiệm / Số Froude |
$U$ (m/s) / $Fn$ |
Dải vận tốc khai thác |
$Fn = 0.377$ (tốc độ thiết kế) |
Nội dung chính theo từng chương
Chương 1. Tổng quan các công trình nghiên cứu liên quan
Chương 1 tổng thuật chi tiết các nghiên cứu về mũi quả lê từ các thử nghiệm mô hình kinh điển (D. Taylor, T. Inui, M. Kracht) đến các nghiên cứu hiện đại sử dụng công cụ CFD (OpenFOAM, SHIPFLOW, CAESES) và các thuật toán tối ưu hóa (Kriging, RBF, GA). Chương này cũng điểm lại tình hình nghiên cứu CFD và thiết kế tàu cá tại Việt Nam (các công trình tại ĐHBK Hà Nội, ĐH Hàng hải, ĐH Giao thông vận tải TP.HCM, ĐH Nha Trang).
Từ việc phân tích tổng quan, tác giả xác định 4 vấn đề kỹ thuật trọng tâm:
- Thiếu các nghiên cứu tối ưu hóa mũi quả lê chuyên biệt cho tàu cá vốn có tỷ số $L/B$ nhỏ, thân ống ngắn và hoạt động đa chế độ.
- Các công trình CFD thường bỏ qua bước chuẩn hóa thông số mô phỏng đầu vào dẫn đến khó kiểm soát sai số.
- Đồ thị Kracht cần được mở rộng để tính toán cho tàu có $C_B < 0.56$.
- Bài toán tối ưu hóa quả lê tàu cá bắt buộc phải là bài toán đa mục tiêu có kèm các điều kiện ràng buộc hàng hải.
Chương 2. Tính sức cản tàu bằng CFD
Chương 2 giải quyết bài toán xác định sức cản tàu bằng mô phỏng số CFD trên phần mềm XFlow với mô hình rối SST $k-\omega$, thiết lập quy trình đảm bảo độ tin cậy của lời giải số.
Quy trình thực hiện:
- Xây dựng mô hình 3D: Tác giả sử dụng AutoShip phục dựng đường hình 3D của tàu FAO 72 và FAO 75 từ bản vẽ mặt cắt ngang, kiểm tra sai lệch hình học bề mặt, sau đó xuất file định dạng STL sang phần mềm XFlow.
| Hạng mục kiểm tra hình học |
Tàu thực tế FAO 75 |
Mô hình 3D AutoShip |
Sai lệch (%) |
| Thể tích chiếm nước $\nabla$ ($m^3$) |
225.00 |
225.12 |
+0.05% |
| Diện tích mặt ướt $S$ ($m^2$) |
205.40 |
205.18 |
-0.11% |
| Hoành độ tâm nổi $LCB$ (m) |
-0.48 |
-0.48 |
0.00% |
- Xác định kích thước miền tính toán: Tác giả khảo sát các phương án kích thước không gian miền tính để triệt tiêu ảnh hưởng dội sóng từ các biên. Miền tính tối ưu được xác lập: phía trước mũi cách $1.5 L_{pp}$, phía sau đuôi cách $3.0 L_{pp}$, hai bên mạn cách $1.5 L_{pp}$, đáy cách $1.5 L_{pp}$, và chiều cao không khí phía trên là $1.0 L_{pp}$.
- Hiệu chỉnh tham số mô hình rối: Xác lập các hệ số thực nghiệm cho mô hình SST $k-\omega$ ($k_1, \omega_1, \omega_2, \sigma_{k1}, \sigma_{\omega1}, P_k$) phù hợp với dòng chảy bao quanh tàu cá vỏ thép.
- Kiểm chứng độ chính xác: So sánh giá trị sức cản tổng $R_{XF}$ (tính từ XFlow) với $R_{tn}$ (đo từ bể thử NPL/SSPA) ở các chế độ vận tốc:
| Mẫu tàu thử nghiệm |
Vận tốc / Trường hợp |
Sức cản bể thử $R_{tn}$ (KG) |
Sức cản XFlow $R_{XF}$ (KG) |
Sai lệch $\Delta R$ (%) |
| FAO 75 |
Trường hợp I (vận tốc thấp) |
1420.5 |
1455.2 |
+2.44% |
| FAO 75 |
Trường hợp I (vận tốc thiết kế) |
3850.0 |
3740.0 |
-2.85% |
| FAO 72 |
Trường hợp II (dải vận tốc khai thác) |
2150.0 |
2230.0 |
+3.72% |
Tất cả các sai số đều nằm trong giới hạn cho phép ($\Delta R < 5%$), khẳng định độ tin cậy của mô hình CFD đã thiết lập.
Chương 3. Thiết kế tối ưu mũi quả lê tàu cá
Chương 3 là nội dung trọng tâm của luận án, thực hiện tính toán thiết kế quả lê ban đầu, tích hợp hình học và giải bài toán tối ưu hóa.
-
Mở rộng đồ thị thiết kế Kracht:
Tàu FAO 75 có hệ số béo $C_B = 0.524$ nằm ngoài dải tra cứu của đồ thị Kracht ($0.56 - 0.82$). Tác giả đã thiết lập các đường cong nội suy và ngoại suy thể hiện mối quan hệ giữa các hệ số hình học quả lê với hệ số béo ở số Froude thiết kế $Fn = 0.377$. Từ đó, xác định được kích thước quả lê ban đầu cho tàu FAO 75:
- Chiều dài quả lê ban đầu: $L_{PR0} = 1.50\text{ m}$
- Chiều rộng lớn nhất quả lê ban đầu: $B_{B0} = 1.70\text{ m}$
- Chiều cao đỉnh quả lê ban đầu: $Z_{B0} = 2.10\text{ m}$
-
Tích hợp hình học quả lê vào thân tàu:
Sử dụng AutoShip hiệu chỉnh các đường sườn và đường sống mũi, tạo vùng chuyển tiếp trơn tru giữa quả lê và vỏ tàu, đảm bảo các thông số hình học của quả lê không bị biến dạng và diện tích mặt ướt tăng ở mức tối thiểu.
-
Thiết lập ma trận phương án và giải bài toán tối ưu:
- Các biến thiết kế: Chiều dài ($L_{PR}$), chiều rộng ($B_B$) và chiều cao ($Z_B$) của quả lê.
- Các trường hợp khảo sát thay đổi kích thước đồng thời: $(L_{PR}, B_B)$, $(L_{PR}, Z_B)$, và $(B_B, Z_B)$ với các bước nhảy gia số $\Delta L_{PR}, \Delta B_B, \Delta Z_B$.
- Hàm đa mục tiêu $\Delta P_e$ được đánh giá qua 3 chế độ: $U_1$ (chạy hành trình), $U_2$ (chạy dắt lưới), và $U_3$ (kéo/thả lưới) ở các mớn nước khảo sát (bao gồm $T = 4.57\text{ m}$).
- Sử dụng mô hình thay thế (Surrogate Model) dạng bề mặt đáp ứng qua 2 bước lặp để tìm điểm cực trị của hàm mục tiêu.
| Thông số hình học |
Quả lê ban đầu (Kracht mở rộng) |
Quả lê tối ưu (CFD + Mô hình thay thế) |
Mức thay đổi |
| Chiều dài quả lê $L_{PR}$ (m) |
1.50 |
1.65 |
+0.15 m |
| Chiều rộng quả lê $B_B$ (m) |
1.70 |
1.91 |
+0.21 m |
| Chiều cao quả lê $Z_B$ (m) |
2.10 |
2.10 |
0.00 m |
| Độ giảm sức cản tổng $\Delta R_T$ (%) |
~8.5% |
~14.0% |
Tăng hiệu quả giảm cản thêm 5.5% |
Phân tích trường dòng chảy quanh thân tàu cho thấy mũi quả lê tối ưu làm giảm biên độ sóng mũi, hạn chế hiện tượng sóng vỡ cục bộ và hạ thấp gradient áp suất tại khu vực sống mũi.
Chương 4. Kết luận và khuyến nghị
Chương 4 tổng hợp các kết quả nghiên cứu lý thuyết và ứng dụng, khẳng định tính đúng đắn của phương pháp đề xuất và đưa ra các khuyến nghị ứng dụng thực tiễn.
Kết quả và những đóng góp mới
Đóng góp mới về mặt khoa học (Lý luận)
- Quy trình chuẩn hóa mô phỏng CFD cho tàu cá: Xây dựng cơ sở phương pháp luận ứng dụng mô hình RANSE và mô hình rối SST $k-\omega$ để tính toán sức cản cho tàu cá vỏ thép chạy chậm, thân béo, thân ống ngắn. Xác định bộ thông số đầu vào chuẩn (kích thước miền tính, điều kiện biên, hệ số mô hình rối) đảm bảo độ sai lệch so với thực nghiệm bể thử dưới 5%.
- Bổ sung và hoàn thiện phương pháp đồ thị Kracht: Xây dựng các đường cong nội suy và ngoại suy cho phép áp dụng phương pháp Kracht để thiết kế mũi quả lê cho tàu có hệ số béo nhỏ ($C_B = 0.524$) nằm ngoài giới hạn áp dụng truyền thống ($C_B < 0.56$), kết hợp kỹ thuật tích hợp trơn đều quả lê vào thân tàu 3D trên AutoShip.
- Mô hình toán tối ưu hóa mũi quả lê đa mục tiêu: Đề xuất mô hình toán học tối ưu hóa hình dạng quả lê với hàm đa mục tiêu về độ giảm công suất có ích ($\Delta P_e$) tương ứng với các chế độ khai thác thực tế của tàu cá, kết hợp các điều kiện ràng buộc thủy tĩnh ($\Delta C_B, \Delta LCB, \Delta MG$) và giải bài toán bằng sự kết hợp giữa CFD và mô hình thay thế (Surrogate Model).
Đóng góp mới về mặt thực tiễn
- Số liệu thiết kế tối ưu cho mẫu tàu cá FAO 75: Xác định bộ kích thước hình học tối ưu cho mũi quả lê tàu FAO 75 ($L_{PRop} = 1.65\text{ m}, B_{Bop} = 1.91\text{ m}, Z_{Bop} = 2.10\text{ m}$), giúp giảm khoảng 14% sức cản tổng của tàu so với khi không trang bị quả lê.
- Ứng dụng trong công nghiệp đóng tàu cá: Cung cấp tài liệu thiết kế và công cụ tính toán phục vụ các cơ sở thiết kế đóng mới tàu cá vỏ thép cỡ lớn tại Việt Nam, nâng cao hiệu quả tiết kiệm nhiên liệu và cải thiện tính năng lắc dọc.
- Phục vụ đào tạo và nghiên cứu: Kết quả luận án cung cấp tài liệu tham khảo chuyên sâu cho công tác giảng dạy sau đại học và nghiên cứu khoa học chuyên ngành Kỹ thuật Tàu thủy.
Kiến nghị và giải pháp đề xuất
- Áp dụng quy trình tính toán CFD và phương pháp tối ưu hóa mũi quả lê của luận án vào việc thiết kế các mẫu tàu cá vỏ thép xa bờ có chiều dài trên 40 m tại Việt Nam.
- Kết hợp công cụ AutoShip và phần mềm CFD thương mại/mã nguồn mở trong các viện nghiên cứu và nhà máy đóng tàu để giảm thiểu chi phí thử nghiệm mô hình trong bể thử.
Hạn chế và hướng nghiên cứu tiếp
Hạn chế
- Mô hình CFD được giới hạn trong điều kiện tàu chuyển động thẳng đều trong môi trường nước tĩnh, chưa xét đến tác động đồng thời của phổ sóng ngẫu nhiên trên biển, gió và sức cản không khí.
- Nghiên cứu tập trung vào chế độ bơi hoàn toàn theo định luật Archimede ở vùng nước sâu vô hạn, chưa khảo sát ảnh hưởng của hiệu ứng nước nông.
- Quá trình tối ưu hóa chủ yếu biến thiên 3 kích thước hình học chính ($L_{PR}, B_B, Z_B$) của mũi quả lê, chưa kết hợp biến dạng tự do đồng thời toàn bộ đường hình phần mũi và phần đuôi tàu.
Hướng nghiên cứu tiếp
- Mở rộng bài toán mô phỏng CFD và tối ưu hóa hình dạng quả lê trong điều kiện tàu chuyển động trên sóng biển thực tế (bài toán Seakeeping, đánh giá chuyển động lắc dọc và va đập sóng slamming).
- Nghiên cứu ảnh hưởng của mũi quả lê đến dòng chảy vào đĩa chân vịt, xác định sự thay đổi của hệ số hút ($t$) và hệ số dòng theo ($w$) để đánh giá toàn diện hiệu suất đẩy của tàu.
- Ứng dụng quy trình tối ưu hóa cho các loại tàu chuyên dụng khác của Việt Nam (tàu kiểm ngư, tàu vận tải ven biển).
Giá trị tham khảo
Luận án mang lại giá trị tham khảo chuyên môn cho nhiều đối tượng:
- Nghiên cứu sinh, học viên cao học và giảng viên ngành Kỹ thuật Tàu thủy / Thủy khí kỹ thuật: Tham khảo phương pháp thiết lập miền tính, xử lý lưới và hiệu chỉnh tham số mô hình rối SST $k-\omega$ trong CFD; phương pháp xây dựng mô hình thay thế (Surrogate Model) trong tối ưu hóa thủy động lực học.
- Kỹ sư thiết kế vỏ tàu tại các viện nghiên cứu, nhà máy đóng tàu: Sử dụng phương pháp ngoại suy đồ thị Kracht và bảng kích thước hình học tối ưu để thiết kế mũi quả lê cho tàu cá vỏ thép cỡ vừa và lớn.
- Cơ quan quản lý nhà nước (Tổng cục Thủy sản, Cục Đăng kiểm Việt Nam): Cơ sở dữ liệu khoa học hỗ trợ xây dựng quy chuẩn, tiêu chuẩn kỹ thuật thiết kế tàu cá vỏ thép xa bờ tiết kiệm năng lượng.
Câu hỏi thường gặp
1. Tại sao đồ thị Kracht (1978) kinh điển không thể áp dụng trực tiếp cho tàu cá FAO 75?
Đồ thị Kracht được xây dựng trên cơ sở dữ liệu thực nghiệm của các mẫu tàu hàng có hệ số béo thể tích nằm trong khoảng $C_B \in [0.56, 0.82]$. Mẫu tàu cá FAO 75 có hệ số béo $C_B = 0.524$, nằm ngoài biên dưới của đồ thị. Luận án đã giải quyết hạn chế này bằng cách thiết lập các đường cong nội suy và ngoại suy toán học tại số Froude thiết kế $Fn = 0.377$ để xác định các hệ số hình học quả lê ban đầu.
2. Mô hình CFD của luận án đã được kiểm chứng độ chính xác như thế nào?
Mô hình CFD (giải hệ phương trình RANSE với mô hình rối SST $k-\omega$ trên XFlow) được kiểm chứng bằng cách so sánh kết quả tính sức cản tổng với dữ liệu đo thực nghiệm trong bể thử NPL (Anh) và SSPA (Thụy Điển) của hai mẫu tàu FAO 72 và FAO 75. Sai lệch giữa kết quả mô phỏng số và thực nghiệm ở các dải vận tốc đều đạt dưới 5% (cụ thể ở FAO 75 là $+2.44%$ và $-2.85%$; FAO 72 là $+3.72%$).
3. Sự khác biệt giữa hàm mục tiêu tối ưu hóa của luận án với các nghiên cứu trước đây là gì?
Các nghiên cứu trước đây thường sử dụng hàm đơn mục tiêu (giảm sức cản tổng tại một vận tốc cố định), chỉ phù hợp với tàu hàng chạy đơn chế độ. Luận án xây dựng hàm đa mục tiêu cực đại hóa độ giảm công suất có ích trung bình $\Delta P_e$ kết hợp trọng số qua 3 chế độ khai thác đặc thù của tàu cá: chạy hành trình ($U_1$), chạy dắt lưới ($U_2$), và kéo/thả lưới ($U_3$), đồng thời tích hợp các điều kiện ràng buộc về ổn định hàng hải ($\Delta C_B, \Delta LCB, \Delta MG$).
4. Kích thước và hiệu quả của mũi quả lê tối ưu trên tàu FAO 75 đạt được là bao nhiêu?
Quả lê tối ưu đạt được có chiều dài $L_{PRop} = 1.65\text{ m}$ (tăng 0.15 m so với ban đầu), chiều rộng $B_{Bop} = 1.91\text{ m}$ (tăng 0.21 m), và chiều cao $Z_{Bop} = 2.10\text{ m}$ (giữ nguyên). Mũi quả lê tối ưu giúp giảm sức cản tổng của tàu khoảng 14% so với dạng mũi thẳng ban đầu.
Kết luận
Luận án của nghiên cứu sinh Huỳnh Văn Chính đã giải quyết trọn vẹn bài toán ứng dụng phương pháp CFD kết hợp mô hình thay thế để thiết kế tối ưu hóa mũi quả lê cho tàu cá vỏ thép. Công trình đã mở rộng phạm vi ứng dụng của phương pháp đồ thị Kracht cho các tàu có hệ số béo thấp ($C_B = 0.524$) và chuẩn hóa quy trình mô phỏng số RANSE với mô hình rối SST $k-\omega$ đạt sai số dưới 5% so với thực nghiệm bể thử. Kết quả tối ưu hóa trên tàu FAO 75 đạt mức giảm 14% sức cản tổng, cung cấp cơ sở khoa học và thực tiễn phục vụ công tác thiết kế hiện đại hóa đội tàu cá vỏ thép tại Việt Nam.
