Tổng quan luận án
Luận án tiến sĩ với đề tài "Ứ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 nghiên cứu sinh Huỳnh Văn Chính thực hiện tại Trường Đại học Giao thông vận tải TP. Hồ Chí Minh, dưới sự hướng dẫn khoa học của PGS.TS Trần Gia Thái và TS. Bùi Hồng Dương, thuộc ngành Kỹ thuật cơ khí động lực (mã số: 9520116), hoàn thành năm 2022.
Về tính cấp thiết và khoảng trống nghiên cứu, việc trang bị mũi quả lê phía trước thân tàu là giải pháp khí động - thủy động học nhằm tạo ra hệ thống sóng phụ giao thoa ngược pha với hệ thống sóng mũi tàu, từ đó làm giảm thành phần sức cản sinh sóng và giảm từ 10% đến 15% sức cản tổng của tàu biển. Đối với tàu cá, mũi quả lê còn có tác dụng cải thiện độ chúi dọc và giảm dao động lắc dọc khi hoạt động trên biển. Tuy nhiên, nếu hình dạng và vị trí quả lê không phù hợp, hiện tượng giao thoa tiêu cực sẽ làm gia tăng sóng tổng hợp và làm tăng sức cản chuyển động. Tại Việt Nam, triển khai Nghị định 67/2014/NĐ-CP đã thúc đẩy việc đóng mới hàng loạt tàu cá vỏ thép dưới 30 m và định hướng phát triển tàu cá trên 40 m. Phần lớn các mẫu tàu cá hiện hành được thiết kế dạng mũi thẳng truyền thống, chưa trang bị mũi quả lê và chưa qua thử nghiệm mô hình trong bể thử do chi phí tốn kém. Mặt khác, các phương pháp thiết kế kinh điển như đồ thị Kracht (1978) bị giới hạn ở dải hệ số béo lớn ($C_B = 0{,}55 - 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{,}524$), đồng thời các nghiên cứu tối ưu hóa hiện đại thường chỉ áp dụng hàm đơn mục tiêu sức cản trên tàu hàng mà chưa xét đến đặc thù đa chế độ khai thác của tàu cá.
Mục tiêu chung của luận án là ứng dụng phương pháp CFD để thiết kế và xác định hình dạng hình học tối ưu của mũi quả lê cho tàu cá vỏ thép nhằm đạt độ giảm sức cản tổng lớn nhất, đồng thời thỏa mãn các điều kiện ràng buộc về hình học và tính năng hàng hải. Luận án đặt ra ba mục tiêu cụ thể:
- Phân tích tổng quan các công trình nghiên cứu liên quan đến phương pháp CFD, tính toán sức cản và tối ưu hóa mũi quả lê nhằm xác lập cơ sở phương pháp luận phù hợp.
- Ứng dụng CFD ước tính sức cản của tàu nghiên cứu đạt độ chính xác mong đợi (sai lệch dưới 5% so với thực nghiệm) làm cơ sở cho bài toán tối ưu.
- Xây dựng mô hình toán và phương pháp tối ưu hóa hình dạng mũi quả lê phù hợp với điều kiện khai thác thực tế của tàu cá.
Đối tượng nghiên cứu của luận án là các mẫu tàu 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ể gồm mẫu tàu FAO 72 và FAO 75 đã được thử nghiệm sức cản tại các bể thử quốc tế NPL (Anh) và SSPA (Thụy Điển). Phạm vi nghiên cứu được giới hạn ở chế độ tàu bơi theo định luật Archimede, chuyển động thẳng đều trong môi trường nước tĩnh không giới hạn độ sâu, không xét tác động của sóng gió khí quyển; lưu chất được giả định là đồng chất, có tính nhớt và không nén được; các dữ liệu thực nghiệm bể thử có sẵn được dùng làm chuẩn kiểm chứng và hiệu chỉnh mô hình tính.
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ử phát triển và các công trình nghiên cứu mũi quả lê qua hai giai đoạn:
Giai đoạn nghiên cứu truyền thống dựa trên thử nghiệm bể thử: Bắt đầu từ ứng dụng đầu tiên của D. Taylor (1910) trên tàu quân sự USS Delaware, tàu khách Bremen và Europa (Đức, 1920), Hoover và Coolidge (1931), Normandie (1935), tàu tuần dương Yamato (1940) và Yamashiro Maru (1963). Các dạng mũi quả lê truyền thống gồm dạng ống trụ tròn kiểu C (Đại học Michigan) và kiểu bóng đèn tròn 3A, 3B, 3C. Năm 1962, T. Inui phát triển lý thuyết tàu không sóng (waveless form) bằng cách bố trí các điểm kỳ dị toán học (singularities) để triệt tiêu sóng mũi và lái, áp dụng thử nghiệm trên tàu Murasaki và Kurenai Maru (giảm 13,1% công suất máy chính). Năm 1978, M. Kracht công bố công trình thực nghiệm tổng hợp trên nhiều mô hình quả lê, thiết lập hệ thống đồ thị thiết kế Kracht ứng với hệ số béo $C_B = 0{,}55 - 0{,}82$.
Giai đoạn nghiên cứu hiện đại dựa trên công cụ CFD và kỹ thuật tối ưu hóa:
- Phương pháp phương án: Grzegorz Filip và cộng sự tại Đại học Michigan phối hợp Đăng kiểm Mỹ (ABS) (2014) tối ưu hóa mũi quả lê cho tàu container KCS của KRISO (Hàn Quốc) bằng phần mềm OpenFOAM, mô hình rối SST $k-\omega$, khảo sát 9 phương án kích thước chiều cao và chiều rộng, sau đó sử dụng mô hình thay thế nội suy Kriging với đa thức bậc hai để xây dựng bề mặt đáp ứng tìm điểm giảm công suất tối đa.
- Phương pháp tham số: Weilin Luo và Linqiang Lan tại Đại học Phúc Châu (2016) tối ưu hóa mũi quả lê cho tàu Ro-Ro ($L = 165{,}7\text{ m}$, $C_B = 0{,}65$, $Fn = 0{,}28$) bằng cách kết hợp CAESES và ShipFlow, mô tả hình học bằng đường cong B-Spline và mặt NURBS với 8 biến thiết kế ($D_{cf}, H_f, F_{tl}, C_{ftl}, F_{ll}, C_{fll}, F_{hb}, B_h$), làm giảm hệ số sức cản sóng $C_W$ đi 4,320%. Nghiên cứu trên tạp chí Journal of Marine Science and Application (2012) áp dụng thuật toán di truyền GA trên nền tảng MATLAB kết hợp RANSE SST $k-\omega$ để tối ưu tàu Seri 60 ($L_{pp} = 7{,}37\text{ m}$, $C_B = 0{,}60$). Phương pháp hàm cơ sở xuyên tâm RBF (Radial Basis Function) được áp dụng để tối ưu 4 mẫu tàu của tổ chức JHSS (ST, EB, GB, BB).
- Các nghiên cứu tại Việt Nam: Viện Cơ khí Động lực (Trường ĐH 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 tàu 6 bậc tự do; Trường ĐH Hàng Hải Việt Nam và ĐH Giao thông vận tải TP.HCM nghiên cứu tương tác dòng chảy chân vịt - bánh lái và lực bẻ lái; Trường ĐH Nha Trang có luận án của Lê Văn Toàn ứng dụng OpenFOAM tính sức cản tàu cá vỏ gỗ và các nghiên cứu tối ưu LCB bằng phương pháp chuyển đổi hình học Lackenby kết hợp CFD; PGS.TS Đỗ Quang Khải xây dựng hàm đa thức phân đoạn biểu diễn đường nước và mũi quả lê trong đồ án thiết kế tàu thủy.
Khoảng trống nghiên cứu được luận án xác định: Phần lớn các nghiên cứu hiện đại tập trung vào tàu hàng có hệ số béo lớn, tỷ số $L/B > 6{,}0$, thân ống dài, dải số Froude hẹp và tải trọng ít thay đổi. Tàu cá vỏ thép có tỷ số $L/B$ nhỏ, thân ống ngắn, chạy chậm và hoạt động ở nhiều chế độ vận tốc, mớn nước khác nhau. Đồ thị Kracht không bao phủ hệ số béo của tàu cá FAO 75 ($C_B = 0{,}524 < 0{,}55$) và không hướng dẫn kỹ thuật ghép nối quả lê vào thân tàu. Đồng thời, việc tối ưu hóa quả lê tàu cá đòi hỏi hàm đa mục tiêu công suất thay vì hàm đơn mục tiêu sức cản đơn thuần.
Cơ sở lý thuyết và phương pháp nghiên cứu
Luận án sử dụng khung lý thuyết động lực học lưu chất tính toán (CFD) dựa trên hệ phương trình Navier-Stokes trung bình Reynolds (RANSE) cho dòng chảy nhớt, không nén được:
$$\frac{\partial U_i}{\partial x_i} = 0$$
$$\rho \left( \frac{\partial U_i}{\partial t} + U_j \frac{\partial U_i}{\partial x_j} \right) = -\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]$$
Mô hình rối được lựa chọn là mô hình Shear Stress Transport ($SST\ k-\omega$), kết hợp ưu điểm của mô hình $k-\omega$ ở vùng sát thành vỏ tàu và mô hình $k-\epsilon$ ở vùng dòng chảy tự do xa vỏ tàu. Phương pháp số sử dụng kỹ thuật thể tích hữu hạn (FVM) tích hợp trong phần mềm thương mại XFlow để giải trường vận tốc và trường áp suất.
Hình học mũi quả lê được xác định qua 6 tham số kích thước cơ bản: chiều dài nhô ra $L_{PR}$, chiều rộng lớn nhất $B_B$, chiều cao $Z_B$, diện tích mặt cắt ngang tại trụ mũi $A_{BT}$, diện tích hình chiếu cạnh $A_{BL}$, thể tích phần nhô $V_{PR}$; và 6 hệ số hình học tương ứng:
$$C_{LPR} = \frac{L_{PR}}{L_{pp}}, \quad C_{BB} = \frac{B_B}{B}, \quad C_{ZB} = \frac{Z_B}{T_{FP}}, \quad C_{ABT} = \frac{A_{BT}}{A_{MS}}, \quad C_{ABL} = \frac{A_{BL}}{A_{MS}}, \quad C_{\nabla PR} = \frac{V_{PR}}{\nabla}$$
Quy trình phương pháp nghiên cứu gồm:
- Phương pháp mô phỏng số: Xây dựng mô hình hình học 3D bề mặt thân tàu trên phần mềm AutoShip bằng kỹ thuật NURBS, xuất sang định dạng STL để nạp vào phần mềm XFlow. Khảo sát độc lập kích thước miền tính toán và tham số mô hình rối để đảm bảo sai số tính toán sức cản dưới 5% so với dữ liệu thực nghiệm bể thử.
- Phương pháp giải tích và đồ thị: Mở rộng phạm vi đồ thị Kracht thông qua việc thiết lập các đường cong nội suy và ngoại suy giải tích để tính toán bộ thông số quả lê ban đầu ($L_{PRo}, B_{Bo}, Z_{Bo}$) cho tàu có $C_B = 0{,}524$.
- Phương pháp tối ưu hóa: Xây dựng bài toán tối ưu đa mục tiêu với hàm mục tiêu là độ giảm công suất có ích trung bình có trọng số $\Delta P_e$ qua 3 chế độ khai thác của tàu cá (chạy hành trình $U_1$, chạy dắt lưới $U_2$, chạy kéo/thả lưới $U_3$):
$$\max \Delta P_e = \sum_{i=1}^{3} w_i \cdot \Delta P_{ei}$$
trong đó $\Delta P_{ei} = \frac{P_{ei} - P_{ebi}}{P_{ei}} \cdot 100%$ và $\sum w_i = 1$. Thiết lập ma trận các phương án thay đổi kích thước quả lê kết hợp mô hình thay thế (Surrogate Model / Response Surface Method - RSM) để xác định điểm tối ưu toàn cục.
Nguồn dữ liệu sử dụng là hồ sơ thiết kế đường hình và số liệu đo sức cản mô hình trong bể thử của hai mẫu tàu cá FAO 72 và FAO 75 do Tổ chức FAO công bố.
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 trình bày tổng quan lịch sử, cơ chế thủy động của mũi quả lê và phân tích chi tiết các công trình nghiên cứu trong và ngoài nước theo hai nhóm truyền thống và hiện đại. Tác giả phân tích cơ chế giao thoa sóng: khi tàu chuyển động, áp lực nước ở mũi tạo ra hệ thống sóng mũi; quả lê gắn ở chân sống mũi tạo ra hệ thống sóng phụ lệch pha, làm triệt tiêu đỉnh sóng mũi và hạ thấp biên độ sóng tổng hợp, giảm sức cản sinh sóng $R_W$ và sức cản bẻ gãy sóng $R_{WB}$, dù làm tăng diện tích ướt dẫn đến tăng nhẹ sức cản ma sát $R_F$.
| Tác giả / Tổ chức |
Đối tượng tàu |
Phương pháp nghiên cứu |
Kết quả chính |
| M. Kracht (1978) |
Nhiều dãy mô hình tàu biển |
Thực nghiệm bể kéo mô hình |
Thiết lập hệ thống đồ thị thiết kế quả lê với $C_B = 0{,}55 - 0{,}82$ |
| G. Filip và CS (2014) |
Tàu container KCS |
CFD (OpenFOAM, $SST\ k-\omega$) + Kriging |
Tối ưu 9 phương án $(H, C)$, xác định bề mặt đáp ứng giảm công suất |
| W. Luo, L. Lan (2016) |
Tàu Ro-Ro ($L = 165{,}7\text{ m}$) |
CAESES + ShipFlow + B-Spline |
Giảm hệ số sức cản sóng $C_W$ đi 4,320%, lượng chiếm nước đổi 0,112% |
| JMSA (2012) |
Tàu Seri 60 ($L_{pp} = 7{,}37\text{ m}$) |
GA (MATLAB) + CFD (RANSE) |
Giảm sức cản sóng $R_W$ và độ dâng mặt sóng $\sum \eta_i^2$ |
| JHSS (2007) |
4 mẫu tàu JHSS ($L_{pp} = 289{,}79\text{ m}$) |
RBF biến dạng hình học + CFD |
Đánh giá 4 dạng mũi: ST, EB, GB, BB trên lượng chiếm nước 37.074 tấn |
Trên cơ sở phân tích khoảng trống học thuật và thực tiễn đội tàu cá Việt Nam, tác giả xác định hướng nghiên cứu: sử dụng mẫu tàu FAO 72 và FAO 75, ứng dụng CFD kiểm chứng độ chính xác tính sức cản, ngoại suy đồ thị Kracht và tối ưu hóa đa mục tiêu hình dạng quả lê bằng mô hình thay thế.
Chương 2. Tính sức cản tàu bằng CFD
Chương 2 tập trung giải quyết bài toán mô phỏng sức cản vỏ tàu bằng công cụ CFD với độ chính xác kiểm chứng qua thực nghiệm. Tác giả khảo sát các đặc điểm hình học của đội tàu cá Việt Nam (Bảng 2.1) và lựa chọn hai mẫu tàu chuẩn: FAO 72 (Bảng 2.2) và FAO 75 (Bảng 2.3).
Quy trình mô phỏng số được thực hiện qua các bước:
- Xây dựng mô hình 3D vỏ tàu trong AutoShip, hiệu chỉnh lưới mặt cong và kiểm tra các thông số hình học (diện tích ướt $S$, hoành độ tâm nổi $LCB$, thể tích chiếm nước $\nabla$) so với số liệu thiết kế gốc (sai số hình học dưới 0,5%).
- Xuất bề mặt vỏ tàu sang định dạng STL và nạp vào phần mềm XFlow.
- Khảo sát kích thước không gian miền tính toán: thiết lập 4 phương án kích thước miền tính (phía trước mũi, phía sau lái, hai bên mạn, phía dưới đáy và phía trên mặt thoáng) để tìm miền tính tối ưu, tại đó giá trị sức cản tổng $R_T$ đạt trạng thái hội tụ và không phụ thuộc vào kích thước biên.
- Xác định điều kiện biên: vận tốc dòng vào (Inlet Velocity), áp suất dòng ra (Outlet Pressure), điều kiện trượt/không trượt tại thành vỏ tàu (Moving Wall / No-slip), mặt thoáng tự do.
- Khảo sát độ nhạy tham số mô hình rối $SST\ k-\omega$ để lựa chọn bộ hệ số tối ưu.
Kết quả tính toán sức cản từ XFlow được đối chiếu trực tiếp với kết quả đo thực nghiệm trong bể thử:
| Mẫu tàu |
Trường hợp thử nghiệm |
Vận tốc tính toán / Số Froude |
Sức cản XFlow $R_{XF}$ (KG) |
Sức cản thực nghiệm $R_{tn}$ (KG) |
Độ sai lệch $\Delta R$ (%) |
| FAO 75 |
Trường hợp I ($T = 4{,}57\text{ m}$) |
Khảo sát theo dải tốc độ $Fn$ |
Số liệu trích xuất từ mô phỏng |
Số liệu bể thử NPL/SSPA |
Sai lệch $< 5{,}0%$ trên toàn dải tốc độ |
| FAO 72 |
Trường hợp II |
Khảo sát theo dải tốc độ $Fn$ |
Số liệu trích xuất từ mô phỏng |
Số liệu bể thử NPL/SSPA |
Sai lệch $< 5{,}0%$ trên toàn dải tốc độ |
Kết quả kiểm chứng tại Chương 2 khẳng định mô hình số thiết lập trên XFlow đạt độ tin cậy cao, sai số nằm trong giới hạn cho phép ($< 5%$), đủ điều kiện làm công cụ đánh giá sức cản trong bài toán tối ưu hóa.
Chương 3. Thiết kế tối ưu mũi quả lê tàu cá
Chương 3 là nội dung cốt lõi của luận án, trình bày quy trình thiết kế quả lê ban đầu và phương pháp tối ưu hóa hình dạng quả lê cho tàu cá FAO 75.
1. Thiết kế mũi quả lê ban đầu bằng đồ thị Kracht mở rộng:
Tàu FAO 75 có hệ số béo thể tích $C_B = 0{,}524$ và số Froude thiết kế $Fn = 0{,}377$. Do $C_B = 0{,}524$ nằm ngoài dải tra cứu của đồ thị Kracht ($C_B = 0{,}55 - 0{,}82$), tác giả đã xây dựng các hàm giải tích nội suy và ngoại suy dựa trên mối quan hệ giữa hệ số giảm công suất dư $C_{\Delta P R}$ cực đại với các hệ số hình học quả lê ($C_{LPR}, C_{BB}, C_{ZB}$). Từ đó tính được các thông số hình học quả lê ban đầu cho tàu FAO 75:
- Chiều dài nhô ra: $L_{PRo} = 1{,}50\text{ m}$
- Chiều rộng lớn nhất: $B_{Bo} = 1{,}70\text{ m}$
- Chiều cao quả lê: $Z_{Bo} = 2{,}10\text{ m}$
Tác giả sử dụng phần mềm AutoShip để dựng đường cong biên dạng dọc, hiệu chỉnh các đường sườn ngang và ghép nối trơn tru bề mặt quả lê vào thân tàu ban đầu, đảm bảo tính liên tục của mặt vỏ và giữ nguyên các thông số hình học đã xác định.
2. Xây dựng mô hình toán tối ưu hóa đa mục tiêu:
- Biến thiết kế: Kích thước hình học của quả lê gồm chiều dài $L_{PR}$, chiều rộng $B_B$, chiều cao $Z_B$.
- Hàm mục tiêu: Cực đại hóa độ giảm công suất có ích tổng hợp $\Delta P_e$ qua 3 chế độ hoạt động chính của tàu cá: chế độ hành trình (vận tốc $U_1$, trọng số $w_1$), chế độ dắt lưới (vận tốc $U_2$, trọng số $w_2$) và chế độ kéo/thả lưới (vận tốc $U_3$, trọng số $w_3$).
- Điều kiện ràng buộc: Biến thiên hệ số béo thể tích $|\Delta C_B| \le \epsilon_1$, độ dịch chuyển hoành độ tâm nổi $|\Delta LCB| \le \epsilon_2$, chiều cao tâm ổn định ban đầu $MG \ge MG_{min}$, và các giới hạn kích thước hình học quả lê nhằm đảm bảo tầm nhìn, độ bền kết cấu và thao tác thả lưới.
3. Phương pháp giải và kết quả tối ưu:
Tác giả thiết lập ma trận các phương án kích thước quả lê (Bảng 3.8) bằng cách thay đổi đồng thời từng cặp kích thước với các gia số xác định:
- Cặp chiều dài - chiều rộng ($L_{PR} - B_B$): Bảng 3.10 đến 3.15
- Cặp chiều dài - chiều cao ($L_{PR} - Z_B$): Bảng 3.16 đến 3.21
- Cặp chiều rộng - chiều cao ($B_B - Z_B$): Bảng 3.22 đến 3.27
Mô phỏng CFD được thực hiện cho các phương án trong ma trận. Dữ liệu công suất thu được được dùng để huấn luyện mô hình thay thế dạng bề mặt đáp ứng (RSM) qua 2 bước lặp tính toán để tìm cực trị chính xác.
| Thông số hình học |
Quả lê ban đầu (Kracht ngoại suy) |
Quả lê tối ưu (CFD + RSM) |
Mức độ thay đổi kích thước |
| Chiều dài quả lê $L_{PR}$ |
$1{,}50\text{ m}$ |
$1{,}65\text{ m}$ |
Tăng $+0{,}15\text{ m}$ |
| Chiều rộng lớn nhất $B_B$ |
$1{,}70\text{ m}$ |
$1{,}91\text{ m}$ |
Tăng $+0{,}21\text{ m}$ |
| Chiều cao quả lê $Z_B$ |
$2{,}10\text{ m}$ |
$2{,}10\text{ m}$ |
Giữ nguyên |
| Độ giảm sức cản tổng $\Delta R_T$ |
- |
Khoảng $14{,}0%$ |
Giảm đáng kể so với tàu gốc |
Phân tích trường dòng chảy cho thấy quả lê tối ưu giúp phân bố lại trường áp suất ở vùng mũi, hạ thấp áp suất đỉnh sóng mũi và làm suy giảm đáng kể hệ thống sóng tổng hợp dọc thân tàu.
Chương 4. Kết luận và khuyến nghị
Chương 4 tổng hợp các kết luận khoa học rút ra từ quá trình nghiên cứu, khẳng định việc hoàn thành đầy đủ các mục tiêu đề ra và đưa ra các khuyến nghị ứng dụng thực tế cho ngành đóng tàu cá vỏ thép tại Việt Nam.
Kết quả và những đóng góp mới
Luận án đã đạt được các đóng góp mới về mặt lý luận khoa học và giá trị thực tiễn:
Về mặt khoa học:
- Xây dựng hoàn chỉnh quy trình ứng dụng lý thuyết CFD để tính toán sức cản cho tàu cá vỏ thép cỡ vừa và nhỏ (chạy tốc độ chậm, thân ống ngắn) có và không có mũi quả lê, kiểm soát các thông số đầu vào (mô hình 3D NURBS, kích thước miền tính, mô hình rối $SST\ k-\omega$) đảm bảo sai số tính toán dưới 5% so với dữ liệu thực nghiệm bể thử.
- Bổ sung và hoàn thiện phương pháp thiết kế quả lê bằng đồ thị Kracht thông qua việc thiết lập các phương trình đường cong nội suy và ngoại suy giải tích, cho phép xác định các thông số quả lê cho tàu có hệ số béo nhỏ ($C_B = 0{,}524$) nằm ngoài dải dữ liệu gốc của Kracht, kết hợp quy trình tích hợp trơn tru bề mặt quả lê vào thân tàu trên phần mềm AutoShip.
- Xây dựng mô hình toán tối ưu hóa đa mục tiêu mũi quả lê dựa trên hàm độ giảm công suất có ích tổng hợp thích ứng với các chế độ khai thác đặc thù của tàu cá, kết hợp ma trận phương án kích thước có ràng buộc kỹ thuật hàng hải và mô hình thay thế (RSM) để giải bài toán tối ưu với chi phí tính toán hợp lý.
Về mặt thực tiễn:
- Xác định được bộ thông số hình học mũi quả lê tối ưu cho mẫu tàu cá vỏ thép FAO 75 ($L_{PRop} = 1{,}65\text{ m}$, $B_{Bop} = 1{,}91\text{ m}$, $Z_{Bop} = 2{,}10\text{ m}$), giúp làm giảm khoảng 14% sức cản tổng của tàu.
- Cung cấp cơ sở khoa học, dữ liệu và công cụ tính toán tin cậy hỗ trợ các đơn vị tư vấn thiết kế, nhà máy đóng tàu trong nước thiết kế mũi quả lê cho các mẫu tàu cá vỏ thép cỡ lớn (chiều dài trên 40 m) phục vụ khai thác xa bờ theo chủ trương hiện đại hóa đội tàu quốc gia.
- Đóng góp tài liệu tham khảo chuyên ngành phục vụ công tác giảng dạy, nghiên cứu sau đại học trong lĩnh vực kỹ thuật tàu thủy và thủy động lực học tàu thuyền.
Các khuyến nghị và giải pháp đề xuất: Áp dụng quy trình tính toán và tối ưu hóa của luận án vào thực tế thiết kế các mẫu tàu cá vỏ thép đóng mới tại Việt Nam; mở rộng áp dụng cho các loại tàu chuyên dụng khác.
Hạn chế và hướng nghiên cứu tiếp
Về mặt hạn chế, các kết quả của luận án được xác lập trong phạm vi giả định: tàu chuyển động thẳng đều trong môi trường nước tĩnh đồng chất, không nén được, độ sâu vô hạn, chưa xét đến tương tác thực tế của sóng biển ngẫu nhiên, gió khí quyển và sức cản không khí phần thượng tầng.
Hướng nghiên cứu tiếp theo được tác giả đề xuất:
- Nghiên cứu mô phỏng số tính toán sức cản và chuyển động lắc của tàu cá có trang bị mũi quả lê trong điều kiện sóng biển thực tế (sóng đều và sóng ngẫu nhiên).
- Tiến hành chế tạo mô hình vật lý và thử nghiệm kéo trong bể thử đối với mẫu tàu cá mang mũi quả lê tối ưu để kiểm chứng thực nghiệm bổ sung trước khi triển khai đóng thực tế.
- Mở rộng phương pháp tối ưu hóa đa mục tiêu cho các dạng mũi tàu và phần đuôi tàu (tối ưu hóa toàn diện đường hình thân tàu cá).
Giá trị tham khảo
Luận án mang lại giá trị tham khảo cụ thể cho các nhóm đố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 cơ khí động lực, Kỹ thuật tàu thủy, Thủy động lực học: Tham khảo quy trình thiết lập bài toán CFD trên phần mềm XFlow, kỹ thuật xử lý miền tính, lựa chọn tham số mô hình rối $SST\ k-\omega$, và phương pháp ứng dụng mô hình thay thế (RSM) trong tối ưu hóa hình học.
- Kỹ sư thiết kế vỏ tàu tại các viện nghiên cứu và nhà máy đóng tàu: Sử dụng phương pháp ngoại suy đồ thị Kracht kết hợp công cụ AutoShip để thiết kế nhanh kết cấu mũi quả lê cho các tàu có hệ số béo nhỏ, đặc biệt là các dòng tàu cá vỏ thép khai thác xa bờ.
- Cơ quan quản lý và hoạch định chính sách ngành thủy sản, đăng kiểm tàu thủy: Tham khảo các căn cứ kỹ thuật về hiệu quả tiết kiệm năng lượng và tính năng hàng hải của mũi quả lê để xây dựng tiêu chuẩn, quy chuẩn thiết kế mẫu tàu cá hiện đại hóa tại Việt Nam.
Câu hỏi thường gặp
1. Tại sao đồ thị thiết kế mũi quả lê Kracht truyền thống không thể áp dụng trực tiếp cho tàu cá FAO 75?
Đồ thị thiết kế của M. Kracht (1978) được xây dựng từ dữ liệu thực nghiệm trên các mẫu tàu hàng có hệ số béo thể tích nằm trong dải $C_B = 0{,}55 - 0{,}82$. Tàu cá FAO 75 có hệ số béo $C_B = 0{,}524$, nằm ngoài phạm vi bao phủ của đồ thị này. Do đó, luận án phải xây dựng các đường cong giải tích nội suy và ngoại suy để xác định các hệ số hình học ban đầu của quả lê.
2. Luận án đã kiểm chứng độ chính xác của phương pháp mô phỏng CFD như thế nào?
Tác giả đã mô phỏng sức cản của hai mẫu tàu FAO 72 và FAO 75 trên phần mềm XFlow (sử dụng mô hình rối $SST\ k-\omega$) và so sánh trực tiếp với kết quả đo thực nghiệm trong bể thử NPL (Anh) và SSPA (Thụy Điển). Kết quả cho thấy độ sai lệch sức cản giữa CFD và thực nghiệm đều nằm dưới ngưỡng 5% trên toàn dải tốc độ khảo sát.
3. Vì sao luận án xây dựng hàm đa mục tiêu thay vì dùng hàm đơn mục tiêu sức cản để tối ưu hóa quả lê tàu cá?
Tàu cá có đặc thù vận hành qua nhiều chế độ khác nhau với vận tốc và mớn nước thay đổi liên tục: chế độ hành trình ($U_1$), dắt lưới ($U_2$) và kéo/thả lưới ($U_3$). Việc tối ưu hóa theo hàm đơn mục tiêu tại một vận tốc cố định có thể làm tăng sức cản ở các chế độ làm việc khác. Hàm đa mục tiêu cực đại hóa độ giảm công suất có ích trung bình có trọng số $\Delta P_e$ phản ánh đúng hiệu quả khai thác thực tế của tàu cá.
4. Kích thước hình học của mũi quả lê tàu FAO 75 thay đổi như thế nào sau khi tối ưu và hiệu quả đạt được là bao nhiêu?
Từ quả lê ban đầu ($L_{PRo} = 1{,}50\text{ m}$, $B_{Bo} = 1{,}70\text{ m}$, $Z_{Bo} = 2{,}10\text{ m}$), phương án tối ưu xác định được kích thước: chiều dài $L_{PRop} = 1{,}65\text{ m}$ (tăng $0{,}15\text{ m}$), chiều rộng $B_{Bop} = 1{,}91\text{ m}$ (tăng $0{,}21\text{ m}$), chiều cao giữ nguyên $Z_{Bop} = 2{,}10\text{ m}$. Kết cấu quả lê tối ưu giúp giảm khoảng 14% sức cản tổng của tàu.
Kết luận
Luận án tiến sĩ của tác giả 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 tính toán động lực học lưu chất (CFD) kết hợp mô hình thay thế để tối ưu hóa hình dạng mũi quả lê cho tàu cá vỏ thép. Công trình đã hoàn thiện phương pháp thiết kế đồ thị Kracht cho dải hệ số béo nhỏ ($C_B = 0{,}524$) và đề xuất mô hình toán tối ưu đa mục tiêu phù hợp với các chế độ khai thác thực tế của tàu cá. Kết quả tối ưu hóa trên mẫu tàu FAO 75 đạt độ giảm 14% sức cản tổng, cung cấp cơ sở khoa học và dữ liệu thiết kế có giá trị ứng dụng cao cho chương trình hiện đại hóa đội tàu cá vỏ thép tại Việt Nam.
