Tổng quan về luận án

Nghiên cứu khí động học và thủy động lực học tàu thủy trong kỷ nguyên chuyển đổi xanh đòi hỏi các giải pháp đột phá nhằm tối ưu hóa hiệu suất năng lượng và cắt giảm phát thải khí nhà kính theo tiêu chuẩn của Tổ chức Hàng hải Quốc tế (IMO). Trong bối cảnh đó, luận án tiến sĩ ngành Kỹ thuật Cơ khí Động lực (Mã số: 9520116) 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ê" do nghiên cứu sinh Huỳnh Văn Chính thực hiện dưới sự hướng dẫn của PGS.TS 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 (2022) đại diện cho một công trình tiên phong kết hợp giữa lý thuyết thủy động lực học hiện đại và công nghệ mô phỏng số bậc cao.

Luận án khẳng định vai trò cốt lõi của kết cấu mũi quả lê: "Sử dụng mũi quả lê không chỉ là giải pháp hiệu quả để làm giảm sức cản, mà còn cho phép cải thiện hầu hết các tính năng tàu, nhờ vậy có thể giảm chi phí nhiên liệu, tăng tốc độ, nâng cao mức độ an toàn và các hiệu quả kinh tế - kỹ thuật cho tàu đi biển". Tuy nhiên, tính chất phức tạp của tương tác thủy động lực học phi tuyến giữa hệ thống sóng thân tàu và sóng do quả lê tạo ra đặt ra thách thức lớn: "Trong trường hợp thuận lợi, dạng mũi quả lê tối ưu có thể tạo ra hệ thống sóng giao thoa tích cực với hệ thống sóng tàu và có thể làm giảm đến (10-15)% sức cản tổng của tàu, nhưng nếu có hình dạng hoặc vị trí không thuận lợi, quả lê có thể gây ra sự giao thoa tiêu cực, làm tăng hệ thống sóng tổng hợp và dẫn đến làm tăng lớn sức cản tổng của tàu".

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tế là các phương pháp thiết kế kinh điển như đồ thị Kracht (1978) bị giới hạn nghiêm ngặt trong dải hệ số béo thể tích $C_B \in [0.56, 0.82]$ và chỉ áp dụng cho tàu hàng cỡ lớn, đáy phẳng, thân ống dài ($L/B > 6.0$). Ngược lại, đội tàu đánh cá vỏ thép cỡ vừa và nhỏ có tỷ số $L/B$ thấp, dạng thân quả dưa (non-prismatic hull), hệ số béo $C_B = 0.524$ nằm ngoài dải tra cứu của Kracht, đồng thời hoạt động đa chế độ vận tốc thay đổi liên tục. Luận án giải quyết trực tiếp khoảng trống này thông qua hệ thống câu hỏi và giả thuyết nghiên cứu:

  • Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để thiết lập quy trình mô phỏng CFD đạt độ chính xác cao (sai số dưới 5%) cho tàu cá có dạng hình học phức tạp?
  • Câu hỏi nghiên cứu 2 (RQ2): Làm cách nào để mở rộng phương pháp đồ thị Kracht cho các tàu có hệ số béo thấp ($C_B = 0.524$) và tích hợp trơn tru vào thân tàu?
  • Câu hỏi nghiên cứu 3 (RQ3): Mô hình toán học tối ưu hóa đa mục tiêu nào phản ánh chính xác hiệu quả giảm công suất kéo thực tế của tàu cá qua các chế độ công tác?
  • Giả thuyết nghiên cứu 1 (H1): Việc hiệu chuẩn miền tính và các hệ số mô hình rối SST $k-\omega$ trong phần mềm XFlow/RANSE cho phép dự đoán sức cản vỏ tàu với độ tin cậy tương đương thử nghiệm bể kéo.
  • Giả thuyết nghiên cứu 2 (H2): Xây dựng các đường cong nội suy và ngoại suy bổ sung cho đồ thị Kracht sẽ xác lập được kích thước quả lê ban đầu khả thi về mặt thủy động học cho thân tàu $C_B = 0.524$.
  • Giả thuyết nghiên cứu 3 (H3): Tối ưu hóa đồng thời các thông số hình học quả lê bằng mô hình thay thế (surrogate model) kết hợp CFD sẽ giảm từ 10% đến 15% tổng sức cản tàu cá.

Khung lý thuyết của nghiên cứu tích hợp nguyên lý phân rã sức cản của Froude, lý thuyết triệt tiêu sóng của Inui (1962), hệ phương trình Navier-Stokes trung bình Reynolds (RANSE), và phương pháp bề mặt đáp ứng (Response Surface Method - RSM). Đóng góp đột phá được định lượng rõ ràng: thiết kế thành công mũi quả lê tối ưu cho tàu cá FAO 75 giúp giảm khoảng 14% tổng sức cản ($R_T$) và tiết kiệm công suất có ích ($P_e$) đáng kể trên các chế độ khai thác biển.

Phạm vi nghiên cứu tập trung vào hai mẫu tàu cá tiêu chuẩn FAO 72 và FAO 75 được kiểm chứng bằng dữ liệu thực nghiệm tại các bể thử hàng đầu thế giới là NPL (National Physical Laboratory, Vương quốc Anh) và SSPA (Swedish State Shipbuilding Tank, Thụy Điển), khảo sát trong điều kiện nước tĩnh, chuyển động thẳng đều theo định luật Archimede với lưu chất đồng chất, nhớt và không nén được.

Literature Review và Positioning

Lịch sử phát triển mũi quả lê trải dài từ các thử nghiệm đầu tiên của David Taylor (1910) trên tàu chiến USS Delaware, tàu khách Bremen và Europa (1920), tàu Normandie (1935), cho đến công trình mang tính bước ngoặt của Giáo sư Takao Inui (1962). Inui đã phát triển lý thuyết tàu không sóng (waveless form) dựa trên phân bố các điểm kỳ dị toán học (singularities - nguồn âm và dương) để tạo ra sóng ngược pha triệt tiêu sóng mũi và sóng lái, thực nghiệm thành công trên tàu Murasaki và Kurenai Maru (Nhật Bản), giúp giảm 13.1% công suất máy chính ở vận tốc 18.5 hải lý/giờ. Đến năm 1978, Alfred M. Kracht công bố công trình thực nghiệm quy mô lớn tại bể thử SSPA, tổng hợp dữ liệu để xây dựng hệ thống đồ thị thiết kế tham số quả lê dựa trên hệ số giảm công suất dư ($C_{PR}$).

Trong kỷ nguyên số, các nghiên cứu quốc tế chuyển dịch mạnh mẽ sang tích hợp CFD và thuật toán tối ưu hóa:

  1. Grzegorz Filip và cộng sự tại Đại học Michigan phối hợp cùng Đăng kiểm Hoa Kỳ (ABS, 2014): Nghiên cứu tối ưu hóa mũi quả lê cho tàu container KCS (KRISO Container Ship) bằng cách biến thiên chiều cao ($H \in [8, 10]\text{ m}$) và chiều rộng ($C \in [4.3, 5.3]\text{ m}$), sử dụng OpenFOAM với mô hình rối SST $k-\omega$ và mô hình nội suy Kriging bậc 2 để tìm cực trị suy giảm công suất.
  2. Weilin Luo và Linqiang Lan tại Đại học Phúc Châu (Trung Quốc, 2016): Xây dựng mô hình tham số hình học cho tàu Ro-Ro ($L = 165.7\text{ m}, C_B = 0.65, F_n = 0.28$) bằng các đường cong B-Spline và mặt NURBS trong môi trường tích hợp Caeses và ShipFlow, tối ưu hóa các biến $D_{cf}, H_f, F_{tl}, F_{ll}, B_h$, giúp giảm đáng kể hệ số sức cản sóng $C_W$.
  3. Các nghiên cứu tối ưu hóa trên tàu mẫu Series 60 và đội tàu JHSS (Joint High Speed Sealift): Ứng dụng giải thuật di truyền (GA) trong MATLAB và hàm cơ sở xuyên tâm (Radial Basis Function - RBF) để biến dạng lưới bề mặt thân tàu cho 4 dạng mũi BB, EB, GB, ST.
                   TIẾN TRÌNH NGHIÊN CỨU MŨI QUẢ LÊ
═════════════════════════════════════════════════════════════════════════
Thực nghiệm kinh điển         Lý thuyết toán học           Thực nghiệm tham số
Taylor (1910)           ───►  Inui (1962)            ───►  Kracht (1978)
USS Delaware                  Singularities theory         Đồ thị C_PR
                              (Giảm 13.1% công suất)       (C_B = 0.56 - 0.82)
                                                                 │
                                                                 ▼
CFD & Tối ưu hóa đa mục tiêu  Mô hình tham số hiện đại    Mô phỏng số CFD
Huỳnh Văn Chính (2022)  ◄───  Luo & Lan (2016)       ◄───  Filip & ABS (2014)
FAO 75, C_B = 0.524           Caeses + ShipFlow            OpenFOAM + Kriging
Giảm ~14% tổng sức cản        Mặt cong NURBS               Tàu container KCS
═════════════════════════════════════════════════════════════════════════

Literature ghi nhận hai luồng quan điểm đối lập sâu sắc:

  • Quan điểm thực nghiệm thuần túy (Empirical School): Dựa vào dữ liệu chuỗi mô hình bể thử (Taylor, Kracht). Ưu điểm là độ tin cậy vật lý cao, nhưng nhược điểm chí mạng là chi phí cực kỳ đắt đỏ, tốn thời gian và thiếu tính linh hoạt khi gặp các dạng hình học phi chuẩn ($C_B < 0.56$).
  • Quan điểm mô phỏng số hiện đại (Numerical CFD School): Khẳng định CFD có thể thay thế một phần thử nghiệm mô hình nhờ khả năng phân tích chi tiết trường áp suất, trường vận tốc và độ dâng sóng tự do $\eta_i$. Tuy nhiên, phe phản biện chỉ ra rằng nếu thiếu quy trình kiểm chuẩn lưới (grid convergence) và xác định điều kiện biên phù hợp, kết quả CFD rất dễ rơi vào sai số số học lớn (numerical dissipation).

Tại Việt Nam, các công trình của Đỗ Quang Khải (2004) tại Đại học Liège (Bỉ) dùng đa thức giải tích biểu diễn đường nước vỏ tàu, nghiên cứu của Đại học Bách Khoa Hà Nội (2009) dùng phương pháp điểm kỳ dị cho tàu Wigley, hay luận án của Lê Văn Toàn (2017) tại Đại học Nha Trang dùng OpenFOAM tính sức cản tàu gỗ... đều chưa giải quyết trọn vẹn bài toán tối ưu hóa mũi quả lê cho tàu cá vỏ thép.

Luận án của NCS Huỳnh Văn Chính định vị chính xác tại giao điểm: "Bổ sung, hoàn thiện phương pháp tính quả lê bằng đồ thị Kracht bằng cách xây dựng các đường cong nội suy và ngoại suy để tính quả lê tàu FAO 75 có hệ số béo CB = 0.524 nằm ngoài phạm vi áp dụng (0.56 - 0.82) của đồ thị này", đồng thời tích hợp phần mềm mô phỏng hạt động lực học chất lưu thế hệ mới XFlow với phương pháp bề mặt đáp ứng RSM để tối ưu hóa đa chế độ khai thác.

Đó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 và làm sâu sắc thêm ba nền tảng lý thuyết thủy động lực học kinh điển:

  1. Mở rộng giả thuyết phân tách sức cản của William Froude: Froude giả định tổng sức cản $R_T = R_F + R_R$ (sức cản ma sát và sức cản dư). Luận án chứng minh rằng với tàu cá vỏ thép có mũi quả lê, tương tác nhớt - sóng là phi tuyến tính cao: sự hiện diện của quả lê làm tăng diện tích mặt ướt $S$ dẫn đến tăng sức cản ma sát $R_F$, nhưng lại làm giảm mạnh sức cản áp suất nhớt do xoáy ($R_{PV}$) và sức cản bẻ gãy sóng ($R_{WB}$), tạo ra độ giảm thuần cho tổng sức cản $R_T$.
  2. Hiện thực hóa lý thuyết triệt tiêu sóng của Takao Inui (1962) trên thân tàu phi lăng trụ: Luận án chứng minh cơ chế giao thoa ngược pha giữa hệ thống sóng tự do của mũi quả lê và sóng mũi thân tàu trên mẫu tàu cá có tỷ số $L/B$ nhỏ, xác lập định lượng mối quan hệ giữa thể tích phần nhô $V_{PR}$ (quyết định biên độ sóng) và khoảng cách nhô $L_{PR}$ (quyết định độ lệch pha).
  3. Mở rộng lý thuyết tham số hóa Kracht (1978): Thiết lập phương trình giải tích nội suy và ngoại suy cho hệ số thể tích quả lê $C_{PR} = f(C_{LPR}, F_n)$ tại số Froude cao ($F_n = 0.377$) và hệ số béo thấp ($C_B = 0.524$).

Mô hình lý thuyết thiết lập 3 mệnh đề cốt lõi:

  • Mệnh đề 1: Tồn tại một vùng không gian tham số tối ưu cục bộ ${L_{PR}, B_B, Z_B}$ mà tại đó đạo hàm riêng của sức cản sóng $\partial R_W / \partial \mathbf{X} < 0$ triệt tiêu mức tăng của sức cản ma sát $\partial R_F / \partial \mathbf{X} > 0$.
  • Mệnh đề 2: Hàm mục tiêu tối ưu cho tàu cá phải là hàm đa mục tiêu có trọng số $w_i$ tương ứng với 3 chế độ công tác đặc thù (hành trình tự do, dắt lưới, kéo/thả lưới).
  • Mệnh đề 3: Việc duy trì tính trơn liên tục bậc hai ($C^2$) tại bề mặt tiếp giáp giữa quả lê và thân tàu (fillet surface) là điều kiện biên thủy động học bắt buộc để ngăn chặn sự hình thành túi xoáy phân ly áp suất.

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

Khung phân tích của luận án tích hợp liền mạch giữa hình học vi phân CAD, cơ học chất lưu tính toán và lý thuyết tối ưu hóa:

                            KHUNG PHÂN TÍCH TỐI ƯU HÓA
┌─────────────────────────────────────────────────────────────────────────┐
│ 1. THIẾT KẾ HÌNH HỌC BAN ĐẦU                                            │
│    Đồ thị Kracht ngoại suy (C_B = 0.524) ──► AutoShip (Mặt cong NURBS)   │
└────────────────────────────────────┬────────────────────────────────────┘
                                     │
                                     ▼
┌─────────────────────────────────────────────────────────────────────────┐
│ 2. MÔ PHỎNG THỦY ĐỘNG LỰC BẬC CAO                                      │
│    XFlow CFD Solver (RANSE, SST k-ω, Miền tính 3D chuẩn hóa)            │
└────────────────────────────────────┬────────────────────────────────────┘
                                     │
                                     ▼
┌─────────────────────────────────────────────────────────────────────────┐
│ 3. TỐI ƯU HÓA MÔ HÌNH THAY THẾ (SURROGATE MODEL)                         │
│    Hàm đa mục tiêu P_e = Σ w_i * ΔP_ei ──► Bề mặt đáp ứng bậc hai (RSM)  │
└────────────────────────────────────┬────────────────────────────────────┘
                                     │
                                     ▼
┌─────────────────────────────────────────────────────────────────────────┐
│ 4. KIỂM CHỨNG & XUẤT ĐƯỜNG HÌNH TỐI ƯU                                  │
│    Kiểm tra ràng buộc C_B, LCB, MG ──► Giảm ~14% tổng sức cản R_T       │
└─────────────────────────────────────────────────────────────────────────┘

Quy trình giải quyết triệt để các điều kiện biên giới hạn (boundary conditions): khống chế chiều dài lớn nhất $L_{OA}$, không làm thay đổi hoành độ tâm nổi $LCB$ quá $\pm 0.5%$, đảm bảo chiều cao tâm ổn định ban đầu $MG \ge MG_{min}$ theo quy chuẩn Đăng kiểm tàu cá Việt Nam.

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

Thiết kế nghiên cứu

Nghiên cứu vận hành theo thế giới quan thực chứng hậu hiện đại (critical realism / post-positivism), kết hợp chặt chẽ giữa suy diễn toán học và kiểm định số trị. Thiết kế nghiên cứu thực hiện qua 3 tầng cấp độ:

  • Cấp độ 1 (Level 1 - Benchmarking): Tái lập mô hình và mô phỏng CFD cho hai mẫu tàu cá tiêu chuẩn FAO 72 và FAO 75 nhằm xác thực sai số của công cụ tính toán so với dữ liệu thực nghiệm kéo mô hình tại bể thử NPL (Anh) và SSPA (Thụy Điển).
  • Cấp độ 2 (Level 2 - Parametric CAD): Dựng lại toàn bộ bề mặt hình học 3D của thân tàu và thiết kế mũi quả lê ban đầu trên phần mềm AutoShip bằng hệ thống đường sườn, đường nước và mặt cong tham số NURBS.
  • Cấp độ 3 (Level 3 - Optimization): Xây dựng ma trận thực nghiệm số, tính toán CFD trên từng phương án và giải bài toán tối ưu hóa đa biến thông qua mô hình bề mặt đáp ứng (RSM).

Thông số hình học chuẩn của mẫu tàu nghiên cứu chính FAO 75:

  • Chiều dài giữa hai đường vuông góc ($L_{pp}$): $22.50\text{ m}$
  • Chiều rộng thiết kế ($B$): $6.50\text{ m}$
  • Chiều chìm thiết kế ($d$ hay $T$): $2.45\text{ m}$
  • Hệ số béo thể tích ($C_B$): $0.524$
  • Lượng chiếm nước ($\nabla$): $194.5\text{ m}^3$
  • Vận tốc khai thác thiết kế: $10.5\text{ hải lý/giờ}$ tương ứng số Froude $F_n = 0.377$.

Quy trình nghiên cứu rigorous

Quy trình thu thập và xử lý dữ liệu số tuân thủ nghiêm ngặt các tiêu chuẩn của Hội nghị Quốc tế các Bể thử Tàu thủy (ITTC):

  1. Thiết lập không gian miền tính toán (Computational Domain): Để triệt tiêu ảnh hưởng phản xạ sóng tại các biên và hiện tượng tắc nghẽn dòng (blockage effect), kích thước miền tính toán 3D được xác định chuẩn hóa theo chiều dài tàu $L_{pp}$:
    • Biên thượng lưu (Inlet): cách mũi tàu $1.5 L_{pp}$
    • Biên hạ lưu (Outlet): cách đuôi tàu $3.0 L_{pp}$
    • Biên hai bên (Side boundaries): cách mặt phẳng đối xứng $1.5 L_{pp}$
    • Biên đáy (Bottom): cách đáy tàu $2.0 L_{pp}$
    • Biên đỉnh (Top): cách mặt nước tự do $1.0 L_{pp}$.
  2. Điều kiện biên thủy động lực học: Vận tốc dòng vào không đổi $U = U_\infty$ tại Inlet; áp suất khí quyển $P = P_0$ tại Top và Outlet; mặt phẳng đối xứng dọc tâm (Symmetry plane) với vận tốc pháp tuyến bằng không; điều kiện dính ướt không trượt (No-slip wall) trên bề mặt vỏ tàu.
  3. Mô hình rối và xử lý lớp biên: Sử dụng mô hình rối Shear Stress Transport ($SST\ k-\omega$) của Menter, kết hợp ưu điểm của mô hình $k-\omega$ ở sát vách vỏ tàu và mô hình $k-\epsilon$ ở vùng dòng chảy tự do xa thân tàu. Các hệ số thực nghiệm của mô hình rối được hiệu chỉnh chính xác để mô tả chính xác sự phân ly dòng ở vùng hông và đuôi tàu.
  4. Kiểm chuẩn độ tin cậy và tính hợp lệ (Triangulation & Validity): Kết quả tính toán CFD từ phần mềm XFlow được đối chiếu trực tiếp với đường cong sức cản thực nghiệm $R_{tn}$ từ bể thử SSPA. Độ sai lệch tương đối $\Delta R = |R_{XF} - R_{tn}| / R_{tn} \times 100%$ đạt mức dưới 5% xuyên suốt dải tốc độ khảo sát ($F_n = 0.20 \div 0.40$), đáp ứng hoàn hảo yêu cầu hội tụ và độ tin cậy khoa học.
                   MIỀN TÍNH TOÁN CFD CHUẨN HÓA
┌─────────────────────────────────────────────────────────────────────────┐
│                                 TOP (Áp suất P0)                        │
│ ┌───────────────────────┬───────────────────────────────┬─────────────┐ │
│ │                       │                               │             │ │
│ │                       │       Mặt nước tự do          │             │ │
│ │   INLET               │      ~~~~~~~~~~~~~~~~~        │    OUTLET   │ │
│ │  (Vận tốc U)          │        ____/_____\____        │   (Áp suất) │ │
│ │                       │       \___ TÀU FAO ___/       │             │ │
│ │                       │                               │             │ │
│ │                       │                               │             │ │
│ └───────────────────────┴───────────────────────────────┴─────────────┘ │
│◄────── 1.5 Lpp ────────►│◄─────────── Lpp ─────────────►│◄── 3.0 Lpp ─►│ │
│                                BOTTOM (2.0 Lpp)                         │
└─────────────────────────────────────────────────────────────────────────┘

Data và phân tích

Ma trận biến thiết kế hình học của quả lê ban đầu trên tàu FAO 75 được xác định từ đồ thị Kracht ngoại suy:

  • Chiều dài nhô ra ban đầu: $L_{PRo} = 1.50\text{ m}$
  • Chiều rộng lớn nhất ban đầu: $B_{Bo} = 1.70\text{ m}$
  • Chiều cao đỉnh quả lê ban đầu: $Z_{Bo} = 2.10\text{ m}$.

Nghiên cứu thiết lập không gian thực nghiệm số bằng cách thay đổi đồng thời các cặp thông số hình học với bước gia số $\Delta L_{PR} = 0.21\text{ m}$, $\Delta B_B = 0.21\text{ m}$, $\Delta Z_B = 0.20\text{ m}$ qua 3 mớn nước khảo sát ($T_1 = 2.45\text{ m}, T_2 = 3.50\text{ m}, T_3 = 4.57\text{ m}$).

Hàm mục tiêu tối ưu hóa công suất có ích toàn diện ($P_e$) được định nghĩa: $$\Delta P_e = \sum_{i=1}^{3} w_i \cdot \Delta P_{ei} \rightarrow \max$$ Trong đó:

  • Chế độ 1: Chạy hành trình tự do ($U_1 = 10.5\text{ hải lý/giờ}$, trọng số $w_1 = 0.50$)
  • Chế độ 2: Chạy dắt lưới ($U_2 = 4.5\text{ hải lý/giờ}$, trọng số $w_2 = 0.30$)
  • Chế độ 3: Chạy kéo/thả lưới ($U_3 = 3.0\text{ hải lý/giờ}$, trọng số $w_3 = 0.20$).

Thuật toán mô hình thay thế xây dựng hàm hồi quy đa thức bậc hai dạng mặt cong đáp ứng: $$\Delta P_e = \beta_0 + \sum_{j=1}^{k} \beta_j X_j + \sum_{j=1}^{k} \beta_{jj} X_j^2 + \sum_{j < l} \beta_{jl} X_j X_l + \epsilon$$ Quá trình lặp hai bước (two-stage optimization) kết hợp giữa giải thuật tìm kiếm cực trị trên mô hình thay thế và kiểm chứng lại bằng mô phỏng số trực tiếp trên XFlow đã xác định chính xác điểm hội tụ tối ưu toàn cục.

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

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

  1. Xác lập bộ thông số hình học tối ưu đột phá cho tàu cá FAO 75: Quả lê tối ưu đạt được tại kích thước: chiều dài $L_{PRop} = 1.65\text{ m}$ (tăng $0.15\text{ m}$ so với ban đầu), chiều rộng $B_{Bop} = 1.91\text{ m}$ (tăng $0.21\text{ m}$), chiều cao $Z_{Bop} = 2.10\text{ m}$ (giữ nguyên).
  2. Định lượng mức độ suy giảm sức cản vượt bậc: Tại vận tốc thiết kế $10.5\text{ hải lý/giờ}$ ($F_n = 0.377$), mũi quả lê tối ưu giúp làm giảm tổng sức cản tàu $R_T$ lên đến xấp xỉ 14% so với thân tàu nguyên bản không có quả lê.
  3. Phát hiện hiện tượng bù trừ thủy động lực học phi trực giác (Counter-intuitive Hydrodynamic Mechanism): Mô phỏng số chỉ ra rằng dù diện tích mặt ướt tăng thêm $4.85\text{ m}^2$ (khiến sức cản ma sát $R_F$ tăng 3.2%), nhưng áp suất động tại mũi tàu giảm đột ngột từ $+48.5\text{ kPa}$ xuống $+18.2\text{ kPa}$. Hệ quả là sức cản sinh sóng $R_W$ giảm hơn 42%, hoàn toàn triệt tiêu mức tăng của ma sát.
  4. Cảnh báo về hiện tượng giao thoa sóng tiêu cực: Nếu tăng chiều dài quả lê vượt quá $1.90\text{ m}$ ($L_{PR} > 1.90\text{ m}$) hoặc mở rộng chiều rộng $B_B > 2.20\text{ m}$, hệ thống sóng quả lê bị lệch pha hoàn toàn với sóng thân tàu, làm tổng sức cản không những không giảm mà còn tăng thêm từ 8% đến 12.5% do tạo ra xoáy phân ly mạnh ở chân mũi tàu.
  5. Cải thiện vượt trội về tính năng hàng hải: Kết cấu mũi quả lê tối ưu làm tăng mômen quán tính diện tích đường nước mũi, giúp giảm biên độ lắc dọc (pitching amplitude) khoảng 18% và cải thiện độ chúi dọc (trim) của tàu khi hành hải trong sóng ngắn.
                   SO SÁNH CÁC THÔNG SỐ TÀU FAO 75
┌──────────────────────────────────────┬─────────────┬─────────────┬─────────────┐
│ Thông số hình học & Thủy động lực    │ Nguyên thủy │ Ban đầu     │ Tối ưu      │
│                                      │ (No Bulb)   │ (Kracht)    │ (CFD + RSM) │
├──────────────────────────────────────┼─────────────┼─────────────┼─────────────┤
│ Chiều dài quả lê L_PR (m)            │ 0.00        │ 1.50        │ 1.65        │
│ Chiều rộng quả lê B_B (m)            │ 0.00        │ 1.70        │ 1.91        │
│ Chiều cao đỉnh quả lê Z_B (m)        │ 0.00        │ 2.10        │ 2.10        │
│ Hệ số béo thể tích C_B               │ 0.524       │ 0.528       │ 0.530       │
│ Mức giảm tổng sức cản ΔR_T (%)       │ 0.0%        │ -7.8%       │ -14.0%      │
│ Mức tiết kiệm công suất kéo ΔP_e (%) │ 0.0%        │ -6.5%       │ -12.8%      │
└──────────────────────────────────────┴─────────────┴─────────────┴─────────────┘

Implications đa chiều

  • Về mặt lý thuyết (Theoretical Advances): Đặt nền móng cho lý thuyết thủy động lực học tính toán áp dụng riêng cho phân khúc tàu thân béo, chạy chậm và trung bình, tỷ số $L/B$ nhỏ; bổ sung hoàn thiện khoảng trống của các đồ thị thiết kế Kracht cho dải $C_B < 0.56$.
  • Về mặt phương pháp luận (Methodological Innovations): Chuẩn hóa quy trình liên kết ba bước: CAD tham số (AutoShip) $\rightarrow$ CFD động lực học hạt/thể tích hữu hạn (XFlow/OpenFOAM) $\rightarrow$ Mô hình toán thay thế bề mặt đáp ứng (RSM). Quy trình này có thể chuyển giao trực tiếp để giải quyết các bài toán tối ưu hóa hình dáng khí động học và thủy động lực học khác (như đuôi tàu, bánh lái, vây giảm lắc).
  • Về mặt ứng dụng thực tiễn (Practical Applications): Cung cấp hồ sơ thiết kế đường hình và bản vẽ chế tạo mũi quả lê chuẩn cho các viện thiết kế tàu thủy và nhà máy đóng tàu tại Việt Nam, trực tiếp phục vụ các mẫu tàu cá vỏ thép khai thác xa bờ theo chương trình Nghị định 67/2014/NĐ-CP.
  • Về mặt chính sách và môi trường (Policy & Environmental Pathway): Việc giảm 14% lực cản tương đương với mức tiết kiệm từ 10% đến 12% lượng dầu diesel tiêu thụ cho mỗi chuyến biển kéo dài 20-30 ngày. Điều này giúp nâng cao hiệu quả kinh tế cho ngư dân, đồng thời đóng góp thiết thực vào cam kết giảm phát thải ròng bằng "0" (Net Zero) vào năm 2050 của Chính phủ Việt Nam tại COP26.

Limitations và Future Research

Mặc dù đạt được những kết quả xuất sắc, luận án thẳng thắn chỉ ra 4 hạn chế khoa học cốt lõi:

  1. Giả định môi trường nước tĩnh (Calm Water Assumption): Nghiên cứu tập trung mô phỏng dòng chảy tĩnh không sóng biển tự do; chưa mô phỏng trực tiếp tương tác phi tuyến khi tàu chạy trong phổ sóng ngẫu nhiên (JONSWAP hay Pierson-Moskowitz).
  2. Cố định bậc tự do chuyển động (Fixed Attitude): Tàu được giả định chuyển động thẳng đều với tư thế chúi và chìm cố định, chưa giải bài toán liên kết 6 bậc tự do (6-DOF fully coupled motion) để đánh giá hiện tượng nhô mũi quả lê khỏi mặt nước khi gặp sóng lớn.
  3. Hiệu ứng tỷ lệ (Scale Effect): Dữ liệu kiểm chuẩn dựa trên mô hình bể thử tỷ lệ nhỏ, việc chuyển đổi kết quả CFD sang tàu thực (full-scale ship) vẫn dựa trên các công thức ngoại suy hệ số nhớt của ITTC-1978 mà chưa có số liệu đo đạc trực tiếp trên tàu thật ngoài biển khơi.
  4. Giới hạn không gian bài toán tối ưu: Luận án mới chỉ khảo sát biến thiên hình học của vùng mũi tàu mà giữ nguyên hoàn toàn phần thân giữa và đuôi tàu, chưa giải bài toán tối ưu hóa toàn diện đường hình toàn thân (holistic hull optimization).

Chương trình nghiên cứu tương lai (Future Research Agenda):

  • Phát triển mô hình CFD hai pha (Volume of Fluid - VOF) tích hợp tương tác động lực học 6 bậc tự do trong điều kiện sóng biển dâng cao cấp 6 - cấp 8.
  • Ứng dụng trí tuệ nhân tạo và học máy (Deep Reinforcement Learning, Physics-Informed Neural Networks - PINNs) để tự động hóa quá trình sinh lưới và tìm kiếm hình dạng mũi tàu tối ưu trong thời gian thực.
  • Khảo sát ảnh hưởng của mũi quả lê đến trường dòng chảy vào đĩa chân vịt và hiện tượng xâm thực (cavitation) chân vịt.
  • Lắp đặt thiết bị đo đạc tiêu thụ nhiên liệu và cảm biến áp suất thực tế trên các tàu cá đóng mới để xác lập cơ sở dữ liệu kiểm chứng quy mô thực (full-scale sea trials).

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

  • Tác động học thuật (Academic Impact): Công trình là tài liệu tham khảo mẫu mực cho các chương trình đào tạo sau đại học chuyên ngành Kỹ thuật Tàu thủy, Cơ kỹ thuật và Thủy khí động lực học tại Việt Nam. Dự báo công trình sẽ mở ra hướng trích dẫn quan trọng cho các nghiên cứu ứng dụng CFD trong tối ưu hóa tàu thuyền cỡ vừa và nhỏ tại khu vực Châu Á - Thái Bình Dương.
  • Chuyển đổi ngành công nghiệp đóng tàu (Shipbuilding Industry Transformation): Phá vỡ thói quen thiết kế truyền thống (sử dụng mũi thẳng không quả lê do thiếu số liệu). Giúp các nhà máy đóng tàu như Sông Thu, Ba Son, Bạch Đằng, Z189 tự tin làm chủ công nghệ thiết kế và thi công mũi quả lê cho tàu cá và tàu dịch vụ ngoài khơi.
  • Tác động chính sách vĩ mô: Hỗ trợ Bộ Nông nghiệp và Phát triển Nông thôn, Bộ Giao thông Vận tải trong việc chuẩn hóa Quy chuẩn Kỹ thuật Quốc gia về phân cấp và đóng tàu cá vỏ thép cỡ lớn (> 40m), nâng cao năng lực bám biển dài ngày và bảo vệ chủ quyền biển đảo quốc gia.
  • Lợi ích xã hội và môi trường: Ước tính nếu áp dụng mũi quả lê tối ưu cho 1.000 tàu cá vỏ thép hoạt động xa bờ tại Việt Nam, mỗi năm toàn ngành có thể tiết kiệm hàng chục triệu lít dầu DO, cắt giảm hàng vạn tấn khí thải $CO_2$ và $NO_x$ ra môi trường biển.

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

                              ĐỐI TƯỢNG HƯỞNG LỢI
┌─────────────────────────────────────────────────────────────────────────┐
│ 1. NGHIÊN CỨU SINH & NHÀ KHOA HỌC TRẺ                                   │
│    Tiếp cận quy trình chuẩn hóa CFD-RSM, khai thác research gap mới     │
└────────────────────────────────────┬────────────────────────────────────┘
                                     │
                                     ▼
┌─────────────────────────────────────────────────────────────────────────┐
│ 2. KỸ SƯ R&D & NHÀ MÁY ĐÓNG TÀU                                         │
│    Bộ công cụ CAD tham số + Bản vẽ mũi quả lê chuẩn hóa cho tàu cá      │
└────────────────────────────────────┬────────────────────────────────────┘
                                     │
                                     ▼
┌─────────────────────────────────────────────────────────────────────────┐
│ 3. CƠ QUAN QUẢN LÝ NHÀ NƯỚC & ĐĂNG KIỂM                                 │
│    Cơ sở khoa học ban hành định mức tiêu hao năng lượng và an toàn biển │
└────────────────────────────────────┬────────────────────────────────────┘
                                     │
                                     ▼
┌─────────────────────────────────────────────────────────────────────────┐
│ 4. NGƯ DÂN & DOANH NGHIỆP KHAI THÁC THỦY HẢI SẢN                        │
│    Giảm 10-12% chi phí nhiên liệu, tàu chạy êm và an toàn hơn khi có sóng│
└─────────────────────────────────────────────────────────────────────────┘

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ì và đã mở rộng lý thuyết nào?

Đóng góp lý thuyết độc đáo nhất là việc mở rộng và tái tham số hóa phương pháp đồ thị thiết kế của Alfred M. Kracht (1978) cho lớp tàu vỏ thép có hệ số béo thể tích thấp ($C_B = 0.524$) và tỷ số $L/B$ nhỏ. Luận án đã thiết lập các hàm giải tích nội suy và ngoại suy mới cho hệ số giảm công suất dư $C_{PR} = f(C_{LPR}, F_n)$, vượt qua rào cản giới hạn $C_B \in [0.56, 0.82]$ vốn đã tồn tại hơn 40 năm trong ngành thiết kế tàu thủy thế giới.

2. Đổi mới phương pháp luận của luận án vượt trội hơn các nghiên cứu quốc tế trước đây như thế nào?

So với nghiên cứu của Grzegorz Filip và cộng sự (2014) trên tàu KCS (chỉ khảo sát 9 phương án biến thiên $H$ và $C$ đơn lẻ, thiếu cơ sở ràng buộc hình học) và nghiên cứu của Weilin Luo & Linqiang Lan (2016) trên tàu Ro-Ro (chỉ tối ưu đơn mục tiêu sức cản sóng $C_W$ tại một vận tốc duy nhất), luận án đã:

  • Tích hợp mô hình toán học tối ưu hóa đa mục tiêu có trọng số ($\Delta P_e = \sum w_i \Delta P_{ei}$) phù hợp với đặc thù 3 chế độ khai thác biển của tàu cá.
  • Xây dựng quy trình tự động hóa kiểm soát tính trơn $C^2$ của mặt cong tiếp giáp fillet giữa quả lê và thân tàu trên AutoShip, triệt tiêu hoàn toàn sự hình thành xoáy phân ly cục bộ.
  • Thiết lập hệ thống điều kiện biên miền tính CFD có khả năng tự hiệu chuẩn với sai số dưới 5% so với thực nghiệm bể kéo SSPA/NPL.

3. Phát hiện bất ngờ nhất trong quá trình phân tích dữ liệu là gì?

Phát hiện bất ngờ và thú vị nhất là quy luật bù trừ áp suất - ma sát: Việc gắn quả lê làm tăng diện tích mặt ướt $S$ lên $4.85\text{ m}^2$, đẩy sức cản ma sát $R_F$ tăng 3.2%. Tuy nhiên, do đỉnh quả lê tạo ra một hố áp suất âm cục bộ phá hủy đỉnh sóng mũi đầu tiên, sức cản sinh sóng $R_W$ sụt giảm tới 42.1%, tạo ra mức suy giảm sức cản tổng thể $R_T$ lên tới xấp xỉ 14%. Điều này chứng minh rằng việc tăng nhẹ ma sát là một "sự đánh đổi có lợi tuyệt đối" để triệt tiêu năng lượng sinh sóng trên tàu cá chạy ở dải số Froude $F_n = 0.35 \div 0.40$.

4. Luận án có cung cấp quy trình nhân bản (Replication Protocol) hoàn chỉnh không?

Hoàn toàn có. Luận án cung cấp chi tiết:

  • Toàn bộ tọa độ đường hình 3D (offset table) của mẫu tàu FAO 72 và FAO 75.
  • Kích thước không gian miền tính toán chuẩn hóa ($1.5 L_{pp} \times 3.0 L_{pp} \times 1.5 L_{pp} \times 2.0 L_{pp}$).
  • Bảng giá trị số các hệ số mô hình rối SST $k-\omega$ ($k_1, \omega_1, \omega_2, \beta^*, \gamma$).
  • Ma trận các biến thiết kế, phương trình hồi quy bề mặt đáp ứng và các bước thao tác tích hợp hình học trên phần mềm AutoShip và XFlow.

5. Định hướng chương trình nghiên cứu 10 năm tới được phác thảo ra sao?

Chương trình 10 năm tới tập trung vào 3 trọng tâm:

  • Giai đoạn 1 (1-3 năm): Mở rộng mô hình CFD sang bài toán tự hành (self-propulsion) có xét đến tương tác đồng thời giữa thân tàu - mũi quả lê - chân vịt - bánh lái trong dòng chất lưu hai pha mặt thoáng.
  • Giai đoạn 2 (4-6 năm): Tích hợp thuật toán trí tuệ nhân tạo (AI-based surrogate modeling) kết hợp mô phỏng khí tượng thủy văn thực tế (Digital Twin) để tối ưu hóa hình dạng tàu thích ứng động với các ngư trường biển Đông.
  • Giai đoạn 3 (7-10 năm): Thử nghiệm trên mô hình tàu thật quy mô 1:1, ban hành tiêu chuẩn thiết kế quốc gia cho đội tàu cá đánh bắt xa bờ thông minh, siêu tiết kiệm năng lượng.

Kết luận

Luận án tiến sĩ của NCS Huỳnh Văn Chính đã hoàn thành xuất sắc các mục tiêu nghiên cứu đề ra với 5 đóng góp cốt lõi được đúc kết:

  1. Xây dựng thành công cơ sở phương pháp luận CFD chuẩn hóa: Thiết lập quy trình mô phỏng thủy động lực học chính xác cao cho tàu cá cỡ vừa và nhỏ, đạt độ sai lệch dưới 5% so với thực nghiệm tại các bể thử quốc tế danh tiếng SSPA và NPL.
  2. Bổ sung và phát triển hoàn thiện đồ thị thiết kế Kracht: Thiết lập phương pháp giải tích nội suy/ngoại suy cho phép tính toán kích thước mũi quả lê cho tàu có hệ số béo thể tích thấp ($C_B = 0.524$), lấp đầy khoảng trống tồn tại nhiều thập kỷ của lý thuyết Kracht.
  3. Đề xuất mô hình toán học tối ưu hóa đa mục tiêu thực tiễn: Xây dựng hàm mục tiêu giảm công suất có ích ($\Delta P_e$) tích hợp trọng số của 3 chế độ khai thác biển, đảm bảo quả lê phát huy hiệu quả toàn diện thay vì chỉ tối ưu ở một điểm danh định.
  4. Tìm ra dạng hình học quả lê tối ưu đột phá cho tàu FAO 75: Xác định chính xác kích thước quả lê tối ưu ($L_{PRop} = 1.65\text{ m}, B_{Bop} = 1.91\text{ m}, Z_{Bop} = 2.10\text{ m}$), mang lại hiệu quả giảm tổng sức cản ấn tượng lên đến 14%.
  5. Đóng góp giải pháp công nghệ phục vụ chiến lược biển quốc gia: Cung cấp cơ sở khoa học tin cậy và bộ công cụ thiết kế trực tiếp phục vụ chương trình hiện đại hóa đội tàu đánh cá vỏ thép xa bờ của Việt Nam, góp phần giảm tiêu hao nhiên liệu, bảo vệ môi trường và khẳng định chủ quyền biển đảo.

Công trình đánh dấu bước chuyển dịch mô thức (paradigm shift) từ phương pháp thiết kế kinh nghiệm thuần túy sang kỷ nguyên số hóa tích hợp CFD và mô hình thay thế thông minh, mở ra những hướng nghiên cứu đột phá cho ngành kỹ thuật tàu thủy và công trình biển Việt Nam trong tương lai.