Tổng quan về luận án

Nghiên cứu hình dáng thân tàu tối ưu nhằm nâng cao hiệu quả năng lượng và giảm thiểu phát thải khí nhà kính đang là trọng tâm cốt lõi của ngành công trình cơ khí động lực thủy khí hiện đại. Trong bối cảnh vận tải thủy nội địa và ven biển Việt Nam đạt sản lượng trên 250 triệu tấn/năm (trong đó vận tải ven biển đóng góp hơn 60 triệu tấn/năm) và ghi nhận tốc độ tăng trưởng luân chuyển hàng hóa ven biển đạt 254%, phương thức vận tải bằng đội tàu container sông pha biển (cấp tàu SB) giữ vai trò huyết mạch trong chuỗi logistics đa phương thức. Tuy nhiên, một nghịch lý công nghệ tồn tại kéo dài là phần lớn đội tàu SB (với quy mô phát triển nhanh từ vài chục chiếc năm 2014 lên gần 2.000 tàu vào cuối năm 2020) chủ yếu được thiết kế dựa trên kinh nghiệm sao chép tàu mẫu hoặc hoán cải cơ học từ sà lan kéo ghép, chưa từng trải qua quy trình tối ưu hóa thủy động lực học có xét đến điều kiện biên luồng lạch hạn chế.

Khoảng trống nghiên cứu (research gap) then chốt nằm ở sự thiếu vắng một khung lý thuyết và công cụ toán học tích hợp có khả năng liên kết đồng thời ba giai đoạn: (1) Tối ưu hóa toàn cục các hệ số hình dáng thân tàu dựa trên hàm mục tiêu sức cản có xét ràng buộc độ sâu luồng cạn; (2) Tự động hóa quá trình tái cấu trúc hình học tuyến hình bằng hàm tham số giải quyết triệt để sự đứt gãy tiếp tuyến giữa các phân đoạn; (3) Kiểm chứng trường dòng thủy động lực học bằng mô phỏng số CFD độ chính xác cao. Để giải quyết khoảng trống này, luận án tập trung vào ba câu hỏi nghiên cứu cốt lõi:

  • RQ1: Làm thế nào để xây dựng hàm mục tiêu sức cản bán thực nghiệm tích hợp chính xác các hiệu ứng gia tăng lực cản trong vùng nước nông có độ sâu hạn chế ($H/T \le 2,0$) phục vụ giải thuật tối ưu hóa di truyền?
  • RQ2: Cơ chế toán học nào có thể xử lý triệt để hiện tượng mất liên tục hình học ($C^1, C^2$) phát sinh trên đường cong diện tích sườn (SAC) sau khi thực hiện các phép biến đổi tịnh tiến theo phương pháp Lackenby cổ điển?
  • RQ3: Mức độ suy giảm sức cản thực tế và sự phân bố lại trường áp suất, trường dòng bao quanh thân tàu container SB sau tối ưu hóa đạt được định lượng như thế nào dưới kiểm chứng của mô phỏng dòng nhớt RANS?

Hệ thống giả thuyết khoa học tương ứng được xác lập:

  • H1: Việc tích hợp chỉ số Froude độ sâu luồng ($F_{nH}$) và hệ số suy giảm vận tốc tới hạn vào hàm mục tiêu Holtrop–Mennen sẽ dẫn hướng giải thuật di truyền (GA) tìm được bộ thông số hình dáng toàn cục ($L/B, B/T, C_P, C_M, \text{LCB}$) làm giảm từ 3% đến 6% lực cản toàn tàu so với thiết kế truyền thống.
  • H2: Mô hình toán nghịch đảo Non-Uniform B-spline (NUBS) trên cơ sở ma trận không vuông cho phép tái lập đường cong diện tích sườn SAC và hệ thống đường sườn lý thuyết đạt độ trơn liên tục bậc hai ($C^2$) mà không làm sai lệch thể tích chiếm nước danh định ($\pm 1%$).
  • H3: Dưới kiểm chứng CFD RANS, thân tàu tối ưu hóa sẽ làm suy giảm cường độ gradient áp suất cục bộ ở vùng mũi và vùng vách lái, từ đó giảm đồng thời hệ số sức cản ma sát nhớt ($C_V$) và hệ số sức cản sóng ($C_W$) trong vùng luồng cạn.

Khung lý thuyết của công trình được định vị trên nền tảng tích hợp vòng xoắn thiết kế tàu Evans (1959), lý thuyết sức cản bán thực nghiệm Holtrop–Mennen (1982, 1984), lý thuyết biến đổi hình học Lackenby (1950), lý thuyết đường cong tham số B-spline của Rogers (1976), và lý thuyết động lực học dòng chảy nhớt trung bình Reynolds (RANS). Phạm vi nghiên cứu giới hạn trên nhóm tàu chở container sông biển có trọng tải đến 5.000 DWT (điển hình là mẫu tàu 128 TEU), hoạt động trên các tuyến luồng duyên hải và cửa sông phía Nam Việt Nam với độ sâu luồng $H = 2,5 - 4,5\text{ m}$, mớn nước thiết kế $T = 2,5 - 3,5\text{ m}$, tương ứng tỷ số độ sâu trên mớn nước $H/T$ từ 1,2 đến 2,0, vận tốc khai thác từ 10 đến 14 hải lý/giờ, thỏa mãn quy chuẩn kỹ thuật quốc gia QCVN 72:2013/BGTVT.

Literature Review và Positioning

Nghiên cứu thủy động lực học và hình dáng vỏ tàu trong lịch sử phát triển qua ba dòng tư tưởng chính:

Dòng nghiên cứu lý thuyết bán thực nghiệm và hồi quy thống kê: D. Taylor (1915), Guldhammer & Harvald (1974), Pophanken (1939) và đặc biệt là Holtrop & Mennen (1982, 1984) đã xây dựng các phương pháp hồi quy từ hàng ngàn chuỗi mô hình bể thử (như Series 60, BSRA, Wageningen). Phương pháp Holtrop–Mennen trở thành chuẩn mực công nghiệp cho giai đoạn thiết kế sơ bộ nhờ khả năng phân rã sức cản toàn phần $R_T$ thành sức cản ma sát $R_F$, sức cản nhớt hình dáng $R_V$, sức cản sóng $R_W$, sức cản vách lái chìm $R_{TR}$ và sức cản phần phụ $R_{APP}$. Tuy nhiên, hạn chế lớn nhất của trường phái này là giả định tàu hành hải trong môi trường nước sâu vô hạn ($H/T > 3,4$), dẫn đến sai số nghiêm trọng khi áp dụng cho tàu SB chạy trong vùng nước nông có sự xuất hiện của sóng đáy và hiện tượng nghẽn dòng (blockage effect).

Dòng nghiên cứu mô hình hóa toán học hình dáng vỏ tàu: Khởi đầu từ việc Ferguson (1960) và Bézier (1962) đưa ra các hàm cơ sở đa thức, Gordon & Riesenfeld (1974) phát triển B-spline, và Rogers (1976) hoàn thiện giải thuật Non-Uniform Rational B-Spline (NURBS). Trong kỹ thuật tàu thủy, Lackenby (1950) đề xuất phương pháp biến đổi tịnh tiến tọa độ sườn dựa trên đường cong SAC để điều chỉnh lượng chiếm nước ($\nabla$) và hoành độ tâm nổi (LCB). Tuy nhiên, các ứng dụng tiếp theo của Benson (1940), F. Taylor (1962), và Khai et al. (2009) đều bộc lộ khiếm khuyết: biến đổi Lackenby thuần túy tạo ra các điểm gãy đạo hàm bậc nhất ($C^1$) tại vị trí tiếp giáp giữa đoạn thân ống song song (parallel middle body) với phần mũi và phần lái, gây mất trơn cục bộ khi phóng dạng thực tế.

Dòng nghiên cứu tối ưu hóa hình dáng và mô phỏng số CFD: Ray et al. (1995) áp dụng thuật toán tôi luyện mô phỏng (Simulated Annealing) cho tàu container 1.336 TEU; Jun et al. (2004) và Kim et al. (2019) ứng dụng giải thuật di truyền (GA) tối ưu hóa sức cản sóng theo lý thuyết Michell; Tahara et al. (2006) kết hợp GA với phân tích CFD; Han et al. (2012) tối ưu hóa mũi quả lê tàu 14.000 TEU giúp giảm 5,7% sức cản; Cheng et al. (2012) điều khiển trực tiếp các điểm lưới NURBS cho tàu 1.300 TEU làm giảm 12% sức cản. Gần đây, Priftis et al. (2018) tại Đại học Strathclyde đã tích hợp GA vào mô hình Holtrop để tối ưu hóa tàu Post-Panamax đạt mức giảm 3% sức cản.

                    TIẾP CẬN LỊCH SỬ VỀ THIẾT KẾ VỎ TÀU
BÁN THỰC NGHIỆM             MÔ HÌNH TOÁN HỌC                TỐI ƯU HÓA & CFD
Holtrop-Mennen (1982)       Lackenby (1950)                 Ray et al. (1995)
Guldhammer-Harvald (1974)   Rogers NUBS (1976)              Priftis et al. (2018)
              KHOẢNG TRỐNG NGHIÊN CỨU HIỆN HỮU:
    - Thiếu tích hợp hiệu ứng nước nông (H/T <= 2.0) vào hàm mục tiêu GA.
    - Thiếu cơ chế xử lý mất liên tục C1/C2 sau biến đổi Lackenby.
    - Thiếu quy trình đóng kín từ GA -> NUBS nghịch đảo -> RANS CFD.
           ĐỊNH VỊ CỦA LUẬN ÁN (NGUYỄN THỊ NGỌC HOA, 2021)
    Khung tích hợp: GA (Holtrop + Nước cạn) -> Nghịch đảo NUBS ma trận
                    không vuông -> RANS CFD (SST k-w / RSM)

Tranh luận học thuật sâu sắc nhất tập trung vào xung đột phương pháp luận: Tối ưu hóa hình học trực tiếp bằng điểm điều khiển NURBS tự do (Free-form deformation) hay Tối ưu hóa gián tiếp dựa trên đường cong diện tích sườn SAC từ tàu mẫu. Trường phái biến dạng tự do (Cheng et al., 2012) linh hoạt nhưng dễ sinh ra các dạng hình học phi thực tế, vi phạm điều kiện xếp dỡ container tiêu chuẩn và mất ổn định thủy tĩnh sơ bộ. Ngược lại, tiếp cận dựa trên tàu mẫu (Khai et al., 2009) bảo tồn các đặc tính an toàn hàng hải nhưng bị bó hẹp trong không gian nghiệm cục bộ.

Luận án định vị tại điểm giao thoa đột phá: Duy trì tính ưu việt của thiết kế từ tàu mẫu thông qua đường cong SAC, nhưng phá vỡ giới hạn truyền thống bằng cách kết hợp thuật toán tối ưu hóa di truyền đa biến có ràng buộc độ sâu luồng, sau đó áp dụng giải thuật nghịch đảo NUBS ma trận không vuông để tự động xóa bỏ hoàn toàn các điểm gián đoạn $C^1, C^2$. So với hai nghiên cứu chuẩn quốc tế là mô hình tàu container KCS chuẩn mực của Hàn Quốc (KRISO Container Ship, $L_{PP} = 230\text{ m}$) và mẫu tàu container nội địa Rhine 86 của châu Âu (Hekkenberg et al., 2013), nghiên cứu này giải quyết đặc thù hình học có tỷ số $B/T$ rất lớn ($3,2 - 4,5$) và tỷ số $L/B$ nhỏ ($5,0 - 6,5$) điển hình của đội tàu SB Việt Nam trong điều kiện luồng cạn cực đoan.

Đóng góp lý thuyết và khung phân tích

Đóng góp cho lý thuyết

Luận án mở rộng hai lý thuyết nền tảng trong kỹ thuật đóng tàu:

  1. Mở rộng lý thuyết biến đổi hình dáng Lackenby (1950): Thiết lập cơ chế toán học liên tục hóa đường cong SAC bằng cách thay thế quá trình hiệu chỉnh thủ công bằng thuật toán giải tích nghịch đảo tham số, biến đổi hàm diện tích sườn gián đoạn thành một đường cong trơn $C^2$ duy nhất trên toàn bộ chiều dài tàu.
  2. Mở rộng mô hình tính toán sức cản Holtrop–Mennen (1982, 1984): Tích hợp thành công hiệu ứng thủy động lực học vùng nước nông thông qua việc hiệu chỉnh hệ số hình dáng $1+k_1$, bổ sung hàm số phụ thuộc tỷ số Froude độ sâu $F_{nH} = V / \sqrt{gH}$ và tỷ số mớn nước $H/T$, giải quyết sự thiếu hụt của mô hình Holtrop nguyên bản vốn chỉ áp dụng cho vùng nước sâu vô hạn.

Mô hình lý thuyết thiết lập mối quan hệ hàm số đa biến giữa sức cản toàn tàu $R_T$ với các biến số hình học: $$R_T = f\left(L_{PP}, B, T, C_B, C_P, C_M, \text{LCB}, L_{Pf}, L_{Pa}, i_E, H/T, F_{nH}\right)$$ Trong đó, các mệnh đề lý thuyết được chứng minh gồm:

  • Mệnh đề 1: Trong vùng nước nông ($H/T \le 1,5$), việc dịch chuyển hoành độ tâm nổi LCB về phía sau sườn giữa (khoảng $-1,0%$ đến $-1,5% L_{PP}$) sẽ làm giảm sự tích tụ áp suất dư ở khu vực vòm đuôi và hạn chế sự tách dòng cục bộ.
  • Mệnh đề 2: Tồn tại một tập hợp hàm cơ sở NUBS bậc 3 tối ưu với vector nút mở rộng cho phép biểu diễn chính xác $100%$ các điểm kiểm soát của đường sườn thực tế mà không gây ra hiện tượng dao động Runge.
                              KHUNG PHÂN TÍCH TỔNG THỂ

Khung phân tích độc đáo

Khung phân tích của luận án tích hợp ba lý thuyết công cụ:

  1. Lý thuyết tối ưu hóa tiến hóa (Evolutionary Computation): Sử dụng giải thuật di truyền (GA) để quét không gian nghiệm phi tuyến phức tạp, tránh rơi vào các cực tiểu địa phương.
  2. Lý thuyết hình học vi phân tính toán (Computational Differential Geometry): Sử dụng đường cong tham số B-spline phi đồng nhất (NUBS) để xây dựng hệ thống đường nước và đường sườn.
  3. Lý thuyết cơ học lưu chất tính toán (CFD): Giải hệ phương trình Navier–Stokes trung bình Reynolds (RANS) đóng kín bởi mô hình ứng suất dòng rối.

Điều kiện biên xác định tường minh: Thể tích chiếm nước $\nabla$ chỉ được dao động trong khoảng $\pm 1%$; hệ số béo lăng trụ $C_P \in [0,60; 0,75]$; hệ số sườn giữa $C_M \in [0,94; 0,985]$; hoành độ tâm nổi $\text{LCB} \in [-2,5%; +1,5%] L_{PP}$; tỷ số chiều sâu luồng trên mớn nước $H/T \in [1,2; 3,4]$.

Phương pháp nghiên cứu tiên tiến

Thiết kế nghiên cứu

Nghiên cứu được thiết kế theo trường phái thực chứng (Positivism) kết hợp mô hình hóa tính toán thực nghiệm (Computational Empiricism). Thiết kế nghiên cứu đa cấp độ gồm ba cấp bậc liên hoàn:

  • Cấp độ vĩ mô (Macro-level): Phân tích thống kê 18 mẫu tàu container đăng kiểm tại Việt Nam và cơ sở dữ liệu đường thủy nội địa để xác lập miền tham số hình học và điều kiện thủy hải văn biên.
  • Cấp độ trung mô (Meso-level): Xây dựng toán tử di truyền (GA) tích hợp hàm mục tiêu Holtrop cải biên trên phần mềm MATLAB để xác định bộ thông số tối ưu.
  • Cấp độ vi mô (Micro-level): Tái cấu trúc hình học bề mặt bằng giải thuật NUBS và giải bài toán thủy động lực học 3D bằng phương pháp thể tích hữu hạn (FVM) trên phần mềm ANSYS Fluent.

Đối tượng nghiên cứu thực nghiệm số là tàu container 128 TEU chạy tuyến sông biển với các thông số cơ bản ban đầu: Chiều dài thiết kế $L = 68,99\text{ m}$; chiều dài giữa hai đường vuông góc $L_{PP} = 65,00\text{ m}$; chiều rộng thiết kế $B = 12,80\text{ m}$; chiều cao mạn $D = 4,60\text{ m}$; mớn nước thiết kế $T = 3,00\text{ m}$; hệ số béo thể tích $C_B = 0,68$; hệ số béo lăng trụ $C_P = 0,705$; hoành độ tâm nổi $\text{LCB} = -0,38% L_{PP}$; vận tốc thiết kế $V = 10,0\text{ hải lý/giờ}$ ($F_n = 0,202$).

                    CẤU TRÚC NHIỄM SẮC THỂ TRONG GIẢI THUẬT GA

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu được thực hiện qua bốn bước chuẩn hóa:

Bước 1: Trích xuất và số hóa dữ liệu hình học. Sử dụng phương pháp tích phân số Simpson để rời rạc hóa đường cong diện tích sườn (SAC) của tàu mẫu thành 21 mặt cắt sườn lý thuyết (từ sườn số 0 đến sườn số 20), trích xuất chính xác diện tích mặt cắt ngang $A_i$ và diện tích mặt nước $A_{WPi}$.

Bước 2: Vận hành giải thuật tối ưu hóa di truyền (GA).

  • Cơ chế mã hóa: Mỗi cá thể là một vector thực chứa 3 gen biến số $[\text{LCB}, C_P, C_M]$.
  • Quần thể ban đầu: Khởi tạo ngẫu nhiên quy mô $N = 100$ cá thể qua 100 thế hệ tiến hóa.
  • Hàm thích nghi (Fitness Function): Cực tiểu hóa sức cản toàn tàu $R_T(V)$ ở vận tốc 10 hải lý/giờ: $$\text{Fitness} = \min \left[ R_T = R_F(1+k_1) + R_{APP} + R_W + R_B + R_{TR} + R_A \right]$$ Trong đó, hệ số sóng $R_W$ và hệ số hình dáng $k_1$ được điều chỉnh theo độ sâu luồng $H$ thông qua công thức hiệu chỉnh vận tốc tới hạn Schlichting: $$\frac{V_{\infty}}{V} = \left( \frac{A_{\text{mid}}}{b \cdot H} + 1 \right)^m$$
  • Toán tử di truyền: Lựa chọn theo phương pháp bánh xe Roulette (Roulette Wheel Selection); lai ghép số học đa điểm với xác suất $P_c = 0,85$; đột biến Gauss đồng nhất với xác suất $P_m = 0,05$.

Bước 3: Tái cấu trúc tuyến hình bằng giải thuật nghịch đảo NUBS. Đường cong tham số B-spline bậc $k$ được biểu diễn: $$P(u) = \sum_{i=0}^{n} N_{i,k}(u) B_i$$ Trong đó $N_{i,k}(u)$ là hàm cơ sở Cox-de Boor, $B_i$ là các điểm điều khiển (control points). Khi thực hiện hiệu chỉnh Lackenby, đường cong SAC mới xuất hiện $m$ điểm mẫu $D_j (j = 1, \dots, m)$ với $m > n+1$ (hệ phương trình siêu định). Luận án giải bài toán nghịch đảo bằng phương pháp bình phương tối thiểu thông qua ma trận không vuông: $$[N]^T [N] [B] = [N]^T [D] \implies [B] = \left( [N]^T [N] \right)^{-1} [N]^T [D]$$ Giải thuật này đảm bảo độ lệch hình học cực đại $\Delta_{\max} < 0,001\text{ mm}$, khôi phục hoàn toàn tính liên tục $C^2$ tại các điểm giao nối phân đoạn mũi và lái.

Bước 4: Thiết lập mô phỏng số CFD RANS.

  • Miền tính toán (Computational Domain): Theo chuẩn ITTC, biên trước đặt cách mũi tàu $1,5 L_{PP}$, biên sau cách lái $3,0 L_{PP}$, hai biên bên cách đường tim $1,5 L_{PP}$, biên dưới đáy điều chỉnh tương ứng với độ sâu luồng cạn $H/T = 1,2; 1,5; 2,0$ và $\infty$.
  • Mô hình dòng rối: Sử dụng mô hình $k-\omega$ SST (Shear Stress Transport) và mô hình ứng suất Reynolds (RSM) để nắm bắt chính xác hiện tượng tách dòng và xoáy thứ cấp ở vùng vòm đuôi.
  • Chia lưới: Lưới lăng trụ biên (Prism layers) được kiểm soát nghiêm ngặt với $y^+ < 1$ tại bề mặt vỏ tàu, tổng số lượng phần tử lưới dao động từ 3,8 đến 5,5 triệu phần tử sau khi đã vượt qua kiểm tra tính độc lập của lưới (Grid Convergence Index - GCI).

Data và phân tích

Độ tin cậy của mô hình tính toán Holtrop và giải thuật GA được đối chuẩn với dữ liệu thực nghiệm quốc tế của tàu mẫu container KCS ($L_{PP} = 230\text{ m}, B = 32,2\text{ m}, T = 10,8\text{ m}, C_B = 0,651$). Kết quả so sánh giữa phương pháp Holtrop cải biên trong luận án và số liệu thử bể kéo KCS cho thấy khoảng sai lệch chuẩn nằm trong biên độ cực kỳ hẹp $[-3%; +3%]$ xuyên suốt dải số Froude $F_n = 0,10 - 0,28$, xác nhận độ tin cậy tuyệt đối của công cụ toán học trước khi áp dụng vào tàu SB.

Thông số hình học Tàu mẫu ban đầu Tàu sau tối ưu hóa GA Mức độ thay đổi
Chiều dài giữa hai đường vuông góc $L_{PP}$ (m) 65,00 65,00 Không đổi
Chiều rộng thiết kế $B$ (m) 12,80 12,80 Không đổi
Mớn nước thiết kế $T$ (m) 3,00 3,00 Không đổi
Thể tích chiếm nước $\nabla$ ($\text{m}^3$) 1.701,42 1.698,85 $-0,15%$ (Bảo toàn)
Hệ số béo lăng trụ $C_P$ 0,705 0,688 $-2,41%$
Hệ số sườn giữa $C_M$ 0,965 0,982 $+1,76%$
Hoành độ tâm nổi LCB (% $L_{PP}$) $-0,38%$ $-1,15%$ Dịch về phía sau lái
Nửa góc vào nước đường nước thiết kế $i_E$ (độ) $16,5^\circ$ $13,2^\circ$ $-20,00%$

Phát hiện đột phá và implications

Những phát hiện then chốt

Phát hiện 1: Giảm thiểu sức cản toàn tàu vượt trội. Tại vận tốc thiết kế $V = 10\text{ hải lý/giờ}$ ($F_n = 0,202$), hình dáng thân tàu container 128 TEU tối ưu hóa bằng giải thuật GA tích hợp nghịch đảo NUBS đạt mức giảm lực cản toàn phần $R_T$ từ $46,85\text{ kN}$ xuống còn $44,32\text{ kN}$, tương ứng mức suy giảm $5,40%$. Khi vận tốc tăng lên dải khai thác cao $V = 14\text{ hải lý/giờ}$ ($F_n = 0,283$), mức độ giảm lực cản đạt tới $6,18%$, chứng minh hiệu quả giảm cản mở rộng trên toàn bộ dải vận tốc hoạt động.

Phát hiện 2: Cơ chế phân rã và tái cấu trúc thành phần sức cản. Phân tích CFD RANS chỉ ra rằng mức giảm sức cản tổng thể đạt được nhờ hai cơ chế vật lý đồng thời:

  • Thành phần sức cản sóng ($C_W$): Giảm mạnh từ $12,4%$ đến $16,8%$ do việc giảm hệ số béo lăng trụ $C_P$ từ 0,705 xuống 0,688 kết hợp với thu hẹp góc vào nước $i_E$ từ $16,5^\circ$ xuống $13,2^\circ$, triệt tiêu đáng kể hệ sóng mũi bậc một.
  • Thành phần sức cản nhớt ($C_V$): Giảm $3,15%$ nhờ việc tăng nhẹ hệ số sườn giữa $C_M$ lên 0,982 (tăng độ đầy ở khu vực giữa tàu) cho phép làm thon gọn đường dòng ở hai đầu mũi và lái, đồng thời việc dịch chuyển LCB về $-1,15% L_{PP}$ ngăn chặn hoàn toàn hiện tượng phân tách lớp biên sớm tại vùng hông lái.

Phát hiện 3: Tác động của điều kiện luồng hạn chế ($H/T$). Mô phỏng trường dòng trong điều kiện nước nông vạch rõ quy luật thủy động lực học phi tuyến: Khi tỷ số $H/T$ giảm từ điều kiện không hạn chế ($\infty$) xuống vùng nước cạn trung bình ($H/T = 1,5$) và nước rất nông ($H/T = 1,2$), lực cản toàn tàu tăng vọt lần lượt là $24,6%$ và $58,3%$. Nguyên nhân do hiệu ứng thắt hẹp đáy luồng làm gia tăng vận tốc dòng bao dưới đáy tàu, gây sụt giảm áp suất cục bộ (hiện tượng Squat). Tuy nhiên, thân tàu được tối ưu hóa theo mô hình toán NUBS duy trì được trường đẳng áp êm thuận hơn, làm giảm biên độ chúc mũi và giảm $4,8%$ độ gia tăng lực cản do đáy luồng so với tàu mẫu gốc.

Phát hiện 4: Loại bỏ 100% điểm kỳ dị hình học nhờ nghịch đảo NUBS. Khác với phương pháp Lackenby truyền thống luôn tạo ra sự nhảy bậc độ cong (curvature jump) tại mặt cắt tiếp giáp thân ống ($x/L = 0,35$ và $x/L = 0,65$), mô hình nghịch đảo ma trận không vuông NUBS mang lại đường cong SAC và các đường nước thiết kế hoàn toàn liên tục bậc hai ($C^2$). Phân tích độ cong vi phân xác nhận không còn bất kỳ điểm uốn hay điểm kỳ dị cục bộ nào trên toàn bộ vỏ tàu.

                      SO SÁNH CÁC THÀNH PHẦN SỨC CẢN
                (Tàu mẫu 128 TEU vs Tàu tối ưu GA-NUBS tại V = 10 kts)
  Lực cản (kN)

Implications đa chiều

Về mặt lý thuyết: Công trình cung cấp một phương pháp luận hoàn chỉnh kết hợp giải tích - tham số hóa - tối ưu hóa tiến hóa - thủy khí động lực học số trong thiết kế hình dáng phương tiện thủy, đóng góp phương pháp luận chuẩn xác cho các bài toán tối ưu hóa hình học phức tạp chịu ràng buộc vật lý đa trường.

Về mặt phương pháp: Giải thuật nghịch đảo NUBS ma trận không vuông được đóng gói thành công cụ tự động hóa trên MATLAB có khả năng chuyển giao và ứng dụng trực tiếp cho các loại tàu khác (tàu khách, tàu hàng rời, tàu chuyên dụng) mà không đòi hỏi chỉnh sửa cấu trúc giải thuật.

Về mặt thực tiễn kinh tế: Với mức suy giảm lực cản trung bình $5,4%$, công suất động cơ chính yêu cầu giảm tương ứng từ $750\text{ kW}$ xuống khoảng $710\text{ kW}$. Đối với một tàu container 128 TEU hoạt động trung bình 280 ngày/năm trên tuyến Hải Phòng – TP. Hồ Chí Minh hoặc TP. Hồ Chí Minh – Cần Thơ, lượng nhiên liệu dầu diesel (DO/FO) tiết kiệm được đạt xấp xỉ 22,5 đến 28,0 tấn/năm/tàu.

Về mặt chính sách và môi trường: Kết quả nghiên cứu là cơ sở khoa học trực tiếp hỗ trợ Cục Đăng kiểm Việt Nam và Bộ Giao thông Vận tải ban hành các tiêu chuẩn định mức năng lượng, thúc đẩy quá trình chuyển dịch xanh hóa đội tàu SB quốc gia, đóng góp thiết thực vào mục tiêu giảm $71,7\text{ triệu tấn } \text{CO}_2$ phát thải giao thông vận tải vào năm 2030 theo cam kết của Chính phủ.

Limitations và Future Research

Luận án thẳng thắn thừa nhận bốn giới hạn khoa học mang tính biên:

  1. Mô hình thủy động lực học tĩnh: Nghiên cứu mới khảo sát lực cản tàu ở trạng thái chạy thẳng trong nước tĩnh có xét độ sâu luồng, chưa tích hợp bài toán động lực học 6 bậc tự do (6-DOF Seakeeping) và sức cản gia tăng trong sóng bất quy tắc (added resistance in irregular waves) khi tàu hoạt động ở vùng biển hở ven bờ có chiều cao sóng $H_s \le 2,5\text{ m}$.
  2. Độ tin cậy của hàm Holtrop ở vùng nước cực nông: Khi tỷ số $H/T < 1,2$, tính phi tuyến của dòng chảy đáy luồng và hiệu ứng ba chiều trở nên cực kỳ phức tạp, vượt quá khả năng phản ánh chính xác của công thức bán thực nghiệm Holtrop cải biên, đòi hỏi phải dựa hoàn toàn vào mô phỏng CFD bậc cao.
  3. Quy mô mẫu thử: Đối tượng nghiên cứu đại diện tập trung vào phân khúc tàu container đơn chân vịt dưới 5.000 DWT hoạt động tại các luồng sông phía Nam, chưa khảo sát sâu cho nhóm tàu hai chân vịt hoặc tàu có tỷ số chiều dài trên chiều rộng đặc thù của luồng sông phía Bắc.
  4. Hiệu ứng chân vịt - vỏ tàu: Quá trình tối ưu hóa hình dáng thân tàu được tiến hành độc lập với chân vịt, chưa tích hợp mô hình tự đẩy (Self-propulsion) với chuyển động quay thực của chân vịt trong mô phỏng CFD.

Chương trình nghiên cứu tương lai định hướng năm trục mở rộng:

  • Tích hợp bài toán tối ưu hóa đa mục tiêu (Multi-Objective Optimization / NSGA-III) cân bằng giữa ba hàm mục tiêu xung đột: Giảm lực cản – Nâng cao tính ổn định lắc êm trong sóng – Tối đa hóa dung tích chứa container trên boong.
  • Mở rộng mô hình toán học từ đường cong tham số 2D (SAC) lên mô hình bề mặt tham số 3D trực tiếp (NURBS Surface Reverse Engineering).
  • Thiết lập mô phỏng số CFD tương tác vỏ - chân vịt - bánh lái ở chế độ tự hành trong luồng hẹp có bờ dốc không đối xứng.
  • Ứng dụng mạng nơ-ron sâu vật lý (Physics-Informed Neural Networks - PINNs) để xây dựng mô hình xấp xỉ bậc cao (Surrogate model), tăng tốc độ tính toán tối ưu hóa hình dáng theo thời gian thực.

Tác động và ảnh hưởng

Ảnh hưởng học thuật: Nghiên cứu mở ra một hướng đi tiên phong tại Việt Nam trong việc tích hợp sâu công nghệ trí tuệ nhân tạo (AI/GA) và đồ họa máy tính cao cấp vào kỹ thuật động lực học tàu thủy. Kết quả luận án tạo tiền đề cho hàng loạt đề tài nghiên cứu chuyên sâu về thủy khí động lực học tính toán và tự động hóa thiết kế công trình nổi.

Chuyển đổi công nghiệp đóng tàu: Cung cấp cho các viện nghiên cứu thiết kế tàu thủy (như Viện Khoa học Công nghệ Tàu thủy, các công ty tư vấn thiết kế hàng hải) một quy trình chuẩn hóa từ khâu thiết kế sơ bộ đến thiết kế kỹ thuật, rút ngắn thời gian thiết kế tuyến hình từ vài tuần xuống còn vài giờ, đồng thời nâng cao sức cạnh tranh của các nhà máy đóng tàu nội địa.

Lợi ích xã hội và môi trường: Việc tối ưu hóa hình dáng toàn bộ đội tàu SB ước tính có thể cắt giảm hàng chục ngàn tấn khí thải $\text{CO}2$ và bụi mịn $\text{PM}{10}$ mỗi năm. Đây là đòn bẩy kỹ thuật thúc đẩy sự chuyển đổi phương thức vận tải từ đường bộ sang đường thủy nội địa – một phương thức vận tải có suất tiêu hao năng lượng trên tấn-km chỉ bằng $27%$ so với đường bộ.

Đối tượng hưởng lợi

  • Nghiên cứu sinh và Giới học thuật Hàng hải: Tiếp cận phương pháp luận toán học chặt chẽ về lý thuyết nghịch đảo NUBS ma trận không vuông và quy trình kiểm chứng CFD chuẩn ITTC cho vùng nước nông.
  • Kỹ sư và Đơn vị Tư vấn Thiết kế Tàu thủy: Sở hữu bộ công cụ mã nguồn MATLAB và giải thuật tối ưu hóa có thể ứng dụng ngay vào việc thiết kế các mẫu tàu mới tiết kiệm năng lượng, thỏa mãn các tiêu chuẩn kỹ thuật nghiêm ngặt.
  • Doanh nghiệp Vận tải Thủy và Chủ tàu: Giảm thiểu trực tiếp từ 4,5% đến 6,0% chi phí nhiên liệu vận hành hàng năm, tăng tuổi thọ hệ động lực và nâng cao hiệu quả khai thác kinh tế trên các tuyến luồng phức tạp.
  • Cơ quan Quản lý Nhà nước (Bộ GTVT, Cục Đăng kiểm, Cục Đường thủy Nội địa): Có cơ sở khoa học định lượng vững chắc để xây dựng các quy chuẩn kỹ thuật thiết kế mẫu tàu SB chuẩn hóa quốc gia và hoạch định chính sách phát triển kết cấu hạ tầng luồng hàng hải xanh.

Câu hỏi chuyên sâu

1. Đóng góp lý thuyết độc đáo nhất của luận án là gì? Đó là việc thiết lập thành công mô hình toán học nghịch đảo tham số Non-Uniform B-spline (NUBS) dựa trên phương pháp bình phương tối thiểu ma trận không vuông $[N]^T [N] [B] = [N]^T [D]$. Mô hình này giải quyết triệt để vấn đề mất liên tục hình học ($C^1, C^2$) vốn tồn tại cố hữu trong phương pháp biến đổi hình dáng Lackenby (1950) suốt hơn 70 năm qua, tạo ra sự chuyển tiếp trơn tru tuyệt đối giữa đoạn thân ống và hai đầu mũi, lái tàu.

2. Điểm mới về phương pháp luận so với các nghiên cứu tiền nhiệm? Khác với Priftis et al. (2018) hay Han et al. (2012) chỉ tối ưu hóa cục bộ hoặc giả định môi trường nước sâu, luận án đã xây dựng một vòng lặp kín hoàn chỉnh: Tích hợp hàm mục tiêu sức cản Holtrop có ràng buộc vận tốc tới hạn nước cạn ($F_{nH}, H/T$) $\to$ Tối ưu hóa toàn cục bằng giải thuật di truyền (GA) $\to$ Tái cấu trúc hình học $C^2$ bằng nghịch đảo NUBS $\to$ Kiểm chứng dòng chảy 3D RANS CFD với mô hình rối $k-\omega$ SST.

3. Phát hiện bất ngờ nhất từ kết quả thực nghiệm số là gì? Phát hiện bất ngờ nhất là sự dịch chuyển của hoành độ tâm nổi LCB về phía sau lái (từ $-0,38%$ về $-1,15% L_{PP}$) kết hợp tăng nhẹ hệ số sườn giữa $C_M$ (từ 0,965 lên 0,982) không những không làm tăng sức cản ma sát nhớt do tăng diện tích mặt ướt mà ngược lại, làm giảm sức cản nhớt $3,15%$ và sức cản sóng $15,6%$. Điều này giải thích bởi việc tăng diện tích sườn giữa cho phép bóp thon góc vào nước ở mũi ($i_E$ giảm $20%$), tạo ra sự phân bố áp suất dọc thân tàu êm thuận hơn rất nhiều trong điều kiện luồng cạn.

4. Nghiên cứu có cung cấp quy trình tái lập (Replication Protocol) không? Có. Toàn bộ cấu trúc nhiễm sắc thể GA, hàm mục tiêu toán học, vector nút và ma trận hàm cơ sở NUBS bậc 3, cùng với các thông số điều kiện biên chia lưới CFD theo chuẩn ITTC và đề xuất của Hekkenberg cho vùng nước hạn chế đều được trình bày chi tiết và minh bạch trong các chương của luận án, cho phép tái lập $100%$ kết quả tính toán.

5. Chương trình nghị sự nghiên cứu 10 năm tới được phác thảo như thế nào? Chương trình nghiên cứu 10 năm định hướng phát triển hệ thống thiết kế tàu thông minh tự hành (Autonomous Intelligent Ship Design System) dựa trên ba cột mốc: (1) Tích hợp tối ưu hóa đa mục tiêu sức cản - lắc êm - kết cấu nhẹ (2022-2025); (2) Xây dựng mô hình số song sinh (Digital Twin) tương tác vỏ - chân vịt trong luồng cạn thời gian thực (2025-2028); (3) Tự động hóa hoàn toàn quy trình thiết kế vỏ tàu xanh không phát thải sử dụng AI và PINNs (2028-2032).

Kết luận

Luận án tiến sĩ của NCS. Nguyễn Thị Ngọc Hoa đã giải quyết trọn vẹn một bài toán kỹ thuật phức tạp mang tính thời sự cao bằng tư duy khoa học mẫu mực và phương pháp tiếp cận hiện đại:

  1. Tổng kết và lượng hóa thành công các đặc trưng ràng buộc thủy hải văn của hệ thống luồng sông pha biển Việt Nam tác động trực tiếp lên hình dáng tàu container.
  2. Xây dựng thành công hàm mục tiêu sức cản tích hợp hiệu ứng nước nông ($H/T, F_{nH}$) vào giải thuật tối ưu hóa di truyền (GA) trên MATLAB, tìm ra bộ thông số hình học tối ưu giúp giảm $5,40%$ lực cản toàn tàu ở vận tốc thiết kế.
  3. Đề xuất thành công mô hình toán nghịch đảo NUBS ma trận không vuông, xử lý triệt để hiện tượng mất liên tục hình học của phương pháp Lackenby, đảm bảo độ trơn liên tục $C^2$ cho toàn bộ tuyến hình tàu.
  4. Kiểm chứng và làm sáng tỏ bản chất thủy động lực học của thân tàu cải tiến thông qua mô phỏng số 3D RANS CFD, xác nhận sự suy giảm đồng thời của hệ số sức cản sóng $C_W$ và hệ số sức cản nhớt $C_V$.
  5. Đóng gói công cụ tính toán và thiết kế hoàn chỉnh, có khả năng ứng dụng thực tiễn cao, đóng góp nền tảng công nghệ quan trọng cho sự phát triển bền vững của ngành công nghiệp đóng tàu và vận tải thủy Việt Nam trong kỷ nguyên kinh tế biển xanh.