Tổng quan về luận án

Nghiên cứu của tác giả 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 (2022) tập trung giải quyết bài toán thủy động lực học mang tính thời 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ê”. Trong bối cảnh chuyển dịch năng lượng hàng hải quốc tế theo các tiêu chuẩn giảm phát thải khí nhà kính của Tổ chức Hàng hải Quốc tế (IMO) và chủ trương hiện đại hóa đội tàu khai thác xa bờ của Việt Nam (theo định hướng Nghị định 67/2014/NĐ-CP), việc tối ưu hóa đường hình thân tàu để giảm thiểu sức cản và tiết kiệm nhiên liệu đóng vai trò quyết định đến hiệu quả kinh tế - kỹ thuật.

Khoảng trống nghiên cứu (research gap) cốt lõi xuất phát từ thực tiễn thiết kế tàu thủy: các phương pháp giải tích và đồ thị kinh điển, tiêu biểu là đồ thị thiết kế của M. Kracht (1978), chỉ áp dụng hiệu quả cho các mẫu tàu hàng có hệ số béo thể tích trong phạm vi hẹp $C_B \in [0.56, 0.82]$ và tỷ số chiều dài trên chiều rộng $L/B > 6.0$. Ngược lại, các loại tàu chuyên dụng như tàu đánh cá vỏ thép cỡ vừa và nhỏ có tỷ số $L/B$ nhỏ, thân ống ngắn, hệ số béo thấp ($C_B = 0.524$) và vận hành trong dải số Froude biến thiên rộng ($F_n = 0.15 - 0.40$) với nhiều chế độ khai thác (chạy tự do, dắt lưới, kéo thả lưới) hoàn toàn vắng bóng trong các chỉ dẫn thiết kế quả lê chuẩn mực. Bên cạnh đó, các nghiên cứu CFD đương đại chủ yếu tối ưu hóa cục bộ trên thân tàu hàng đã có sẵn quả lê với hàm đơn mục tiêu ở một mớn nước cố định, thiếu vắng quy trình tích hợp trơn tru hình học 3D (CAD-NURBS), kiểm soát sai số mô phỏng CFD và giải bài toán tối ưu đa chế độ vận hành.

Hệ thống câu hỏi và giả thuyết nghiên cứu được xác lập chặt chẽ:

  • RQ1: Làm thế nào để thiết lập quy chuẩn mô phỏng số CFD (kích thước miền tính, điều kiện biên, tham số mô hình rối) nhằm dự báo chính xác sức cản tàu cá vỏ thép với sai số kiểm chứng dưới 5% so với thực nghiệm bể thử?
  • RQ2: Phương pháp mở rộng và ngoại suy đồ thị Kracht có thể áp dụng để xác định các thông số hình học ban đầu cho thân tàu cá có hệ số béo $C_B = 0.524$ nằm ngoài dải chuẩn hay không?
  • RQ3: Thuật toán kết hợp giữa mô phỏng CFD độ chính xác cao và mô hình thay thế (Surrogate Model) đa mục tiêu có khả năng tìm ra cấu hình mũi quả lê tối ưu giúp giảm cực đại công suất có ích trên toàn bộ chu kỳ khai thác của tàu cá hay không?
  • H1: Việc hiệu chỉnh chính xác mô hình rối Shear Stress Transport ($k-\omega$ SST) và miền tính toán trong CFD sẽ đảm bảo độ lệch sức cản tổng giữa mô phỏng và thực nghiệm bể thử (NPL, SSPA) đạt dưới 5%.
  • H2: Hình dạng quả lê ngoại suy từ đồ thị Kracht kết hợp làm trơn biên dạng 3D bằng Non-Uniform Rational B-Splines (NURBS) trong AutoShip cho phép tạo ra sự giao thoa sóng triệt tiêu khả thi mà không làm biến dạng các đặc tính nổi và ổn định ban đầu.
  • H3: Việc tối ưu hóa đồng thời các biến kích thước hình học chính ($L_{PR}, B_B, Z_B$) qua hàm mục tiêu giảm công suất có ích đa chế độ sẽ mang lại mức giảm sức cản tổng trên 10% so với tàu nguyên bản không trang bị quả lê.

Khung lý thuyết của công trình tích hợp lý thuyết triệt tiêu sóng của Inui (1962), lý thuyết phân rã sức cản của Froude, hệ phương trình Navier-Stokes trung bình Reynolds (RANS) kết hợp mô hình rối hai phương trình SST $k-\omega$ của Menter, cùng lý thuyết xấp xỉ bề mặt đáp ứng Kriging. Đóng góp đột phá được định lượng rõ nét: giảm khoảng 14% sức cản tổng của mẫu tàu cá điển hình FAO 75, thiết lập quy chuẩn mô phỏng CFD kiểm chứng thành công trên tập dữ liệu thử nghiệm mô hình bể thử quốc tế (NPL - Anh và SSPA - Thụy Điển) cho hai mẫu tàu FAO 72 và FAO 75. Phạm vi nghiên cứu bao hàm các chế độ chuyển động thẳng đều của tàu trong môi trường nước tĩnh không nén được, ở ba trạng thái mớn nước và dải vận tốc khai thác thực tế.

Literature Review và Positioning

Lịch sử nghiên cứu khí động và thủy động lực học mũi quả lê ghi nhận bước ngoặt từ công trình thử nghiệm của D. Taylor (1910) trên tàu chiến USS Delaware, tiếp nối bằng các ứng dụng trên tàu khách Bremen, Europa (1920) và Normandie (1935). Về mặt lý thuyết giải tích, T. Inui (1962) đã tạo nên cuộc cách mạng với lý thuyết "dạng tàu không sóng" (waveless form), chứng minh rằng việc bố trí mũi quả lê tương đương với việc đưa vào các điểm kỳ dị toán học (singularities - nguồn âm và nguồn dương) tạo ra hệ thống sóng phụ trợ có cùng biên độ nhưng ngược pha $180^\circ$ với hệ thống sóng mũi của thân tàu, từ đó triệt tiêu sóng mũi và làm giảm đáng kể sức cản sinh sóng ($R_W$). Trên bình diện thực nghiệm kinh điển, công trình đồ sộ của M. Kracht (1978) đã phân tích dữ liệu kéo mô hình của hàng trăm biến thể quả lê để xây dựng hệ thống đồ thị thiết kế chuẩn hóa định lượng 6 hệ số không thứ nguyên: hệ số chiều dài ($C_{LPR}$), chiều rộng ($C_{BB}$), chiều cao ($C_{ZB}$), diện tích mặt cắt ngang ($C_{ABT}$), diện tích cạnh bên ($C_{ABL}$) và thể tích ($C_{VPR}$).

Trong y văn tồn tại hai trường phái đối thoại và tranh luận học thuật sâu sắc:

  1. Trường phái Thế lưu (Potential Flow) cổ điển vs. Trường phái Động lực học lưu chất nhớt (Viscous CFD - RANS): Các nghiên cứu truyền thống ứng dụng phương pháp phần tử biên (BEM), phương pháp Panel giải phương trình Laplace giả định chất lỏng lý tưởng phi nhớt (như các bài toán giải tích của Michell hay Wigley). Nhược điểm chí mạng của trường phái này là hoàn toàn bỏ qua sức cản nhớt ($R_V$), lớp biên động năng và hiện tượng bẻ gãy sóng cục bộ ($R_{WB}$) kèm tách dòng ở khu vực mũi tàu. Ngược lại, trường phái RANS hiện đại mô hình hóa đầy đủ tương tác nhớt - sóng, đặc biệt quan trọng đối với các tàu béo chạy chậm ($F_n < 0.2$), nơi sức cản sinh sóng chỉ chiếm phần nhỏ nhưng thành phần sức cản áp suất nhớt ($R_{PV}$) do xoáy mũi lại chiếm tỷ trọng áp đảo.
  2. Trường phái Tối ưu hóa đơn điểm (Single-point/Single-objective) vs. Tối ưu hóa đa chế độ khai thác (Multi-operating condition): Hầu hết các công trình quốc tế tối ưu hóa hình học quả lê tại một vận tốc thiết kế và mớn nước định mức duy nhất. Tuy nhiên, các nhà nghiên cứu thủy động học phản biện rằng một quả lê tối ưu ở chế độ chạy tự do (free running) có thể gây ra giao thoa sóng tiêu cực (tăng sức cản) khi tàu thay đổi mớn nước hoặc khi chuyển sang chế độ kéo lưới công suất cao ở dải tốc độ thấp.

Vị trí của luận án được định vị tiên phong so với các nghiên cứu quốc tế tiêu biểu:

  • So với nghiên cứu của 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) trên tàu container KCS: Nhóm tác giả Filip chỉ khảo sát 9 phương án biến đổi hình học rời rạc trên tàu container có sẵn quả lê, sử dụng phần mềm OpenFOAM với mô hình SST $k-\omega$. Luận án của NCS Huỳnh Văn Chính vượt trội hơn khi giải quyết bài toán thiết kế mới từ đầu cho tàu chưa từng có quả lê, xây dựng cơ sở toán học cho việc mở rộng đồ thị Kracht ngoài biên $C_B = 0.524$, đồng thời xác lập ma trận ràng buộc hình học và ổn định hàng hải chặt chẽ.
  • So với công trình của Weilin Luo và Linqiang Lan (Đại học Phúc Châu, 2017) trên tàu Ro-Ro: Nghiên cứu của Luo & Lan sử dụng phần mềm thương mại CAESES tích hợp ShipFlow để tối ưu hóa tham số hóa đường cong B-Spline và mặt NURBS, nhưng chỉ tập trung vào tối thiểu hóa hệ số sức cản sóng $C_W$ trong điều kiện thế lưu. Luận án của Huỳnh Văn Chính tích hợp mô phỏng nhớt 3D trong XFlow, đánh giá trực tiếp sức cản tổng $R_T$ và chuyển đổi thành hàm công suất có ích đa mục tiêu $P_e$, bảo đảm tính ứng dụng thực tế vượt trội cho đội tàu khai thác thủy hải sản.
       +-------------------------------------------------------------+
       |             TIẾN TRÌNH PHÁT TRIỂN LÝ THUYẾT                 |
       +-------------------------------------------------------------+
                                     |
       [Inui (1962)]                 |           [Kracht (1978)]
       Lý thuyết triệt tiêu sóng     |           Đồ thị thiết kế thực nghiệm
       (Singularities ngược pha)     |           (Giới hạn: C_B = 0.56 - 0.82)
                                     v
       +-------------------------------------------------------------+
       |             TRANH BIỆN PHƯƠNG PHÁP HỌC THUẬT                |
       |  Thế lưu phi nhớt (BEM)   vs.   Dòng nhớt thực tế (RANS CFD)|
       |  Tối ưu đơn điểm (1 Fn)   vs.   Đa chế độ khai thác (Tàu cá)|
       +-------------------------------------------------------------+
                                     |
                                     v
       +-------------------------------------------------------------+
       |                ĐÓNG GÓP CỦA LUẬN ÁN (2022)                  |
       |  1. Ngoại suy phi tuyến Kracht cho C_B = 0.524              |
       |  2. Tích hợp CAD (AutoShip) - RANS CFD (SST k-w) - Kriging  |
       |  3. Tối ưu hóa đa mục tiêu theo chu kỳ khai thác thực tế     |
       +-------------------------------------------------------------+

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

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

Luận án mang lại những bước tiến lý thuyết thực chất trong chuyên ngành Kỹ thuật Tàu thủy và Thủy động lực học tính toán:

  1. Mở rộng lý thuyết thiết kế bán thực nghiệm Kracht (1978): Công trình đã chứng minh tính khả thi và thiết lập các hàm giải tích nội suy/ngoại suy đường cong phi tuyến cho các hệ số hình học quả lê ($C_{LPR}, C_{BB}, C_{ZB}, C_{VPR}, C_{ABT}, C_{ABL}$) tại điểm kỳ dị hình học $C_B = 0.524$ nằm ngoài miền xác định $C_B \in [0.56, 0.82]$. Điều này phá vỡ rào cản tồn tại hơn 4 thập kỷ của phương pháp Kracht đối với các tàu có hệ số béo nhỏ.
  2. Lý thuyết tương tác sóng - nhớt trên thân tàu có tỷ số $L/B$ nhỏ: Luận án bổ sung luận cứ khoa học vào mô hình phân rã sức cản tổng $R_T = R_F(1+k) + R_W + R_{WB}$, chỉ ra rằng ở dải số Froude trung bình thấp của tàu cá, mũi quả lê được tối ưu hóa không chỉ làm giảm biên độ sóng tự do $R_W$ theo nguyên lý giao thoa Inui, mà còn có tác dụng dẫn dòng, giảm độ dốc mặt sóng mũi, triệt tiêu thành phần năng lượng tiêu tán do bẻ gãy sóng $R_{WB}$ và tái phân bố trường áp suất nhớt $R_{PV}$ dọc phần mũi tàu.
  3. Mô hình toán học tối ưu hóa đa mục tiêu có trọng số thích ứng: Luận án đề xuất mô hình giải tích chuyển đổi từ tối ưu hóa sức cản tĩnh sang tối ưu hóa mức giảm công suất có ích tổng hợp $\Delta P_e$: $$\Delta P_e = \sum_{i=1}^{n} w_i \cdot \Delta P_{ei} = \sum_{i=1}^{n} w_i \left( \frac{R_{Ti} - R_{Tbi}}{R_{Ti}} \right) \times 100%$$ trong đó $w_i$ là trọng số thời gian vận hành thực tế của tàu tại từng chế độ công tác (hành trình $w_1 = 0.5$, dắt lưới $w_2 = 0.3$, kéo thả lưới $w_3 = 0.2$), tạo nên bước chuyển dịch mô thức từ tối ưu thuần túy hình học sang tối ưu hóa theo chu kỳ kinh tế - kỹ thuật.

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

Khung phân tích của luận án kết hợp một cách hệ thống ba trụ cột phương pháp luận:

  • Cơ sở hình học 3D tham số hóa (Parametric Geometric Modeling): Sử dụng AutoShip xây dựng bề mặt trơn nhẵn bằng lưới NURBS bậc cao, giải quyết triệt để bài toán nối ghép chuyển tiếp trơn đều (smooth transition surface) giữa phần vỏ thân tàu ban đầu và quả lê thiết kế mới, bảo đảm không tạo ra các điểm gãy góc gây xoáy cục bộ và giữ nguyên thể tích chiếm nước ban đầu.
  • Mô hình động học chất lưu nhớt RANS giải bằng Thể tích hữu hạn (FVM): Thiết lập hệ thống phương trình bảo toàn khối lượng và động lượng dòng chảy hai pha nước - không khí có xét đến mặt thoáng tự do: $$\frac{\partial \rho}{\partial t} + \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] + \rho g_i$$
  • Mô hình bề mặt đáp ứng (Response Surface Methodology - RSM) và Nội suy Kriging: Đóng vai trò là mô hình thay thế (surrogate model), ánh xạ phi tuyến giữa véc-tơ biến thiết kế $\mathbf{X} = [L_{PR}, B_B, Z_B]^T$ và hàm mục tiêu $\Delta P_e$, giúp giảm thiểu số lượng ca tính CFD tốn kém tài nguyên tính toán nhưng vẫn bảo đảm hội tụ chính xác tới cực trị toàn cục.
       +-------------------------------------------------------------+
       |                 KHUNG PHÂN TÍCH TOÀN DIỆN                   |
       +-------------------------------------------------------------+
                                     |
       [BƯỚC 1: HÌNH HỌC]            |           [BƯỚC 2: MÔ PHỎNG CFD]
       - Ngoại suy đồ thị Kracht     |           - Lưới biên dạng FVM/LBM
       - Mô hình hóa NURBS AutoShip  |           - Lời giải RANS SST k-w
       - Nối ghép tiếp xúc trơn G1   |           - Xác lập miền tính & biên
                                     |
                                     +-------------------------------+
                                                                     |
                                                                     v
       +-------------------------------------------------------------+
       |         [BƯỚC 3: MÔ HÌNH THAY THẾ & TỐI ƯU HÓA]             |
       |  - Ma trận biến thiết kế: L_PR, B_B, Z_B                    |
       |  - Xây dựng bề mặt đáp ứng Kriging đa thức bậc 2            |
       |  - Ràng buộc: Delta V ~ 0; Delta LCB <= eps; MG >= MG_min   |
       |  - Hàm mục tiêu: Maximize Delta P_e (Đa chế độ)             |
       +-------------------------------------------------------------+

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

Thiết kế nghiên cứu

Nghiên cứu đứng trên lập trường nhận thức luận Thực chứng khoa học (Positivism) kết hợp với Hiện thực phản biện (Critical Realism), xem các quy luật cơ học thủy khí là khách quan và có thể lượng hóa chính xác thông qua mô phỏng số nếu các điều kiện vật lý biên và mô hình toán được chuẩn hóa. Thiết kế nghiên cứu phân tầng (multi-level design) bao gồm:

  • Tầng 1 - Chuẩn hóa hình học: Dựng lại chính xác mô hình 3D của tàu cá mẫu FAO 72 và FAO 75 từ tập dữ liệu tọa độ sườn thực nghiệm của FAO.
  • Tầng 2 - Hiệu chuẩn và Kiểm chứng số (Verification and Validation - V&V): So khớp kết quả tính toán thủy động lực học trên phần mềm CFD chuyên dụng (XFlow/OpenFOAM) với dữ liệu thử nghiệm mô hình thực tế kéo trong bể thử NPL (National Physical Laboratory - Anh) và SSPA (Swedish State Shipbuilding Tank - Thụy Điển).
  • Tầng 3 - Tối ưu hóa thích ứng: Thiết lập không gian tìm kiếm biến hình học quả lê, xây dựng kế hoạch thực nghiệm số (Design of Experiments - DOE) và tìm điểm tối ưu thông qua mô hình thay thế.
Mẫu khảo sát định lượng:
- Mẫu tàu FAO 75: Chiều dài L_pp = 22.86 m, Chiều rộng B = 6.70 m, Mớn nước d = 2.44 m, Lượng chiếm nước V = 188.4 m3, Hệ số béo C_B = 0.524.
- Mẫu tàu FAO 72: Chiều dài L_pp = 20.00 m, Chiều rộng B = 6.40 m, Mớn nước d = 2.40 m, Hệ số béo C_B = 0.584.

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ nghiêm ngặt các hướng dẫn quốc tế của Hội nghị các Bể thử Tàu thủy Quốc tế (ITTC):

  1. Chiến lược lấy mẫu không gian và thiết lập miền tính: Miền tính toán (computational domain) được tối ưu hóa để loại trừ ảnh hưởng của hiệu ứng thành bể và phản xạ sóng. Khoảng cách từ mũi tàu đến biên cửa vào (Inlet) là $1.5 L_{pp}$; từ đuôi tàu đến biên cửa ra (Outlet) là $3.5 L_{pp}$; độ sâu miền tính dưới đáy tàu là $2.5 L_{pp}$; khoảng cách sang hai bên mạn là $1.5 L_{pp}$; và chiều cao không gian không khí phía trên mặt nước là $1.0 L_{pp}$.
  2. Điều kiện biên vật lý:
    • Biên cửa vào (Inlet): Vận tốc dòng vào xác định $U = U_{ship}$, mức độ rối $I = 1%$, tỷ lệ độ nhớt rối $\mu_t/\mu = 10$.
    • Biên cửa ra (Outlet): Áp suất thủy tĩnh tương ứng với mực nước mặt thoáng tự do ($p = 0\text{ Pa}$ tại mặt thoáng).
    • Bề mặt vỏ tàu (Hull wall): Điều kiện không trượt (No-slip condition), tích hợp hàm thành (Standard Wall Functions) kiểm soát $y^+$ phù hợp.
    • Biên đáy và biên bên: Điều kiện trượt tự do (Free-slip condition) mô phỏng môi trường nước vô tận.
  3. Tam giác đạc phương pháp (Methodological Triangulation): Kết quả CFD được đối chiếu chéo giữa phương pháp lưới thể tích hữu hạn truyền thống (RANS) và phương pháp Boltzmann mạng (LBM) trên phần mềm XFlow, sau đó thẩm định trực tiếp với đường cong thực nghiệm bể thử SSPA/NPL.
       +-------------------------------------------------------------+
       |               KÍCH THƯỚC MIỀN TÍNH TOÁN CFD                 |
       +-------------------------------------------------------------+
                               Top (+1.0 Lpp)
       +-------------------------------------------------------------+
       |                                                             |
       |  Inlet             [MÔ HÌNH VỎ TÀU 3D]               Outlet |
       | (1.5 Lpp) <============ Lpp ============>          (3.5 Lpp)|
       |                                                             |
       +-------------------------------------------------------------+
                              Bottom (-2.5 Lpp)
                      <==== Rộng mạn: 1.5 Lpp mỗi bên ====>

Data và phân tích

Quá trình mô phỏng số áp dụng mô hình rối hai phương trình $k-\omega$ SST với các hằng số thực nghiệm được hiệu chỉnh tối ưu cho bài toán vỏ tàu cá:

  • Phương trình động năng rối $k$: $$\frac{\partial (\rho k)}{\partial t} + \frac{\partial (\rho k u_i)}{\partial x_i} = \frac{\partial}{\partial x_j} \left[ \left( \mu + \frac{\mu_t}{\sigma_k} \right) \frac{\partial k}{\partial x_j} \right] + G_k - Y_k + S_k$$
  • Phương trình tốc độ tiêu tán động năng rối riêng $\omega$: $$\frac{\partial (\rho \omega)}{\partial t} + \frac{\partial (\rho \omega u_i)}{\partial x_i} = \frac{\partial}{\partial x_j} \left[ \left( \mu + \frac{\mu_t}{\sigma_\omega} \right) \frac{\partial \omega}{\partial x_j} \right] + G_\omega - Y_\omega + D_\omega + S_\omega$$

Kết quả kiểm tra tính độc lập của lưới và tham số mô hình rối trên tàu FAO 75 ở trường hợp thử nghiệm I cho thấy sự hội tụ vượt trội. Sai số tương đối giữa sức cản tính toán từ mô phỏng CFD ($R_{XF}$) và sức cản đo trong bể thử ($R_{tn}$) tại dải vận tốc khảo sát đạt mức: $$\Delta R = \left| \frac{R_{XF} - R_{tn}}{R_{tn}} \right| \times 100% < 4.65%$$ Thỏa mãn hoàn hảo tiêu chuẩn kiểm chứng của ITTC ($\Delta R \le 5%$), khẳng định độ tin cậy tuyệt đối của công cụ tính toán trước khi đưa vào bài toán tối ưu hóa.

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

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

  1. Xác lập cấu hình mũi quả lê tối ưu toàn cục cho tàu cá FAO 75: Từ quả lê thiết kế ban đầu dựa trên ngoại suy Kracht có kích thước chiều dài $L_{PRo} = 1.50\text{ m}$, chiều rộng lớn nhất $B_{Bo} = 1.70\text{ m}$, chiều cao $Z_{Bo} = 2.10\text{ m}$, thuật toán tối ưu hóa kết hợp mô hình thay thế hai bước lặp đã xác định được cấu hình quả lê tối ưu:

"Kết quả đã tính được quả lê tối ưu của tàu FAO 75 có chiều dài $L_{PRop} = 1.65\text{ m}$, chiều rộng $B_{Bop} = 1.91\text{ m}$, chiều cao $Z_{Bop} = 2.10\text{ m}$ dựa trên việc thay đổi đồng thời chiều dài và chiều rộng quả lê ban đầu với các gia số $\Delta L_{PR} = 0.21\text{ m}$, và độ giảm sức cản tổng của tàu bằng khoảng 14%."

  1. Quy luật triệt tiêu sóng và biến đổi trường áp suất: Phân tích trường vận tốc và trường áp suất quanh thân tàu trước và sau khi gắn quả lê chứng minh rằng quả lê tối ưu tạo ra một vùng áp suất dương cục bộ ngay dưới đường ngấn nước mũi, triệt tiêu đỉnh sóng mũi đầu tiên ($Wave\ Crest$). Sự giao thoa sóng ngược pha làm giảm độ cao sóng tổng hợp dọc thân tàu từ 18% đến 25%, qua đó giảm sức cản sinh sóng $R_W$.

  2. Hiện tượng phi trực quan về tương tác kích thước quả lê: Dữ liệu thực nghiệm số chỉ ra rằng việc gia tăng thể tích quả lê (tăng đồng thời $L_{PR}$ thêm 10% và $B_B$ thêm 12.35%) tuy làm tăng diện tích mặt ướt $S$ (làm tăng nhẹ sức cản ma sát $R_F$), nhưng mức độ sụt giảm của sức cản áp suất $R_P$ và sức cản bẻ gãy sóng $R_{WB}$ vượt trội hơn rất nhiều. Hiện tượng này giải thích lý do tại sao các quả lê dạng giọt nước đáy đầy (bulb with lower center of area) phát huy hiệu quả giảm cản tốt nhất trên các tàu có mớn nước lớn.

  3. Độ chính xác và tính hội tụ của Mô hình thay thế Kriging: So sánh giá trị độ giảm công suất cực đại giữa mô hình thay thế $\Delta P_{e,\text{max}}^{\text{surrogate}}$ và mô phỏng CFD trực tiếp $\Delta P_{e,\text{max}}^{\text{CFD}}$ trong lần tính thứ hai cho thấy độ lệch tuyệt đối dưới $0.35%$, chứng minh mô hình Kriging bậc hai hoàn toàn đủ năng lực thay thế các vòng lặp CFD phức tạp.

+------------------------------------------------------------------------------------+
|               BẢNG SO SÁNH THÔNG SỐ VÀ HIỆU QUẢ THỦY ĐỘNG LỰC TÀU FAO 75           |
+------------------------------------+-----------------------+-----------------------+
| Thông số khảo sát                  | Trạng thái ban đầu    | Cấu hình tối ưu (CFD) |
+------------------------------------+-----------------------+-----------------------+
| Chiều dài quả lê L_PR (m)          | 1.50                  | 1.65 (+10.0%)         |
| Chiều rộng quả lê B_B (m)          | 1.70                  | 1.91 (+12.35%)        |
| Chiều cao quả lê Z_B (m)           | 2.10                  | 2.10 (Giữ nguyên)     |
| Sức cản tổng R_T tại Fn = 0.377    | 100% (Baseline)       | Giảm ~14.0%           |
| Sai số mô phỏng so với bể thử      | < 5.0%                | Đạt chuẩn ITTC        |
| Tương tác mặt sóng tự do           | Sóng mũi cao, dốc     | Giao thoa triệt tiêu  |
+------------------------------------+-----------------------+-----------------------+

Implications đa chiều

  • Về mặt học thuật và lý thuyết: Luận án cung cấp dữ liệu định lượng và cơ sở khoa học lấp đầy khoảng trống nghiên cứu về thủy động lực học tàu cá vỏ nhỏ có $C_B$ thấp. Mở rộng biên độ áp dụng của đồ thị Kracht từ $C_B \ge 0.56$ xuống $C_B = 0.524$.
  • Về mặt phương pháp luận: Xây dựng hoàn chỉnh chuỗi quy trình khép kín: Số hóa đường hình 3D (AutoShip) $\rightarrow$ Mô phỏng CFD tự động hóa (XFlow/OpenFOAM) $\rightarrow$ Tối ưu hóa bề mặt đáp ứng (Kriging). Quy trình này có thể chuyển giao và ứng dụng trực tiếp cho các bài toán tối ưu hóa hình học phức tạp khác như vây giảm lắc, bánh lái và đạo lưu chân vịt.
  • Về mặt thực tiễn và kinh tế - kỹ thuật: Giảm 14% sức cản tổng đồng nghĩa với việc tiết kiệm trực tiếp từ 10% đến 12% chi phí nhiên liệu dầu DO trong mỗi chuyến biển xa bờ kéo dài 20-30 ngày. Điều này giúp nâng cao biên lợi nhuận cho ngư dân, cải thiện độ chúi dọc và tính năng lắc dọc ($Pitching$) khi tàu chạy ngược sóng.
  • Về chính sách và môi trường: Cung cấp cơ sở khoa học để Tổng cục Thủy sản, Cục Đăng kiểm Việt Nam và Bộ Giao thông Vận tải ban hành các quy chuẩn kỹ thuật thiết kế mẫu tàu cá vỏ thép tiêu chuẩn, hiện thực hóa các cam kết quốc tế về cắt giảm phát thải khí nhà kính $CO_2, NO_x$.

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 ghi nhận các giới hạn nghiên cứu:

  1. Môi trường mô phỏng nước tĩnh: Nghiên cứu mới khảo sát sức cản tàu ở chế độ bơi theo định luật Archimede trong môi trường nước tĩnh, chưa xét đến tác động ngẫu nhiên của phổ sóng biển thực tế ($Irregular\ Waves$), gió và dòng chảy tầng mặt.
  2. Mô hình tàu vỏ cứng không biến dạng và bỏ qua xâm thực: Chưa xét đến tương tác kết cấu - lưu chất (Fluid-Structure Interaction - FSI) khi quả lê chịu va đập sóng lớn ($Slamming$) và nguy cơ xâm thực ($Cavitation$) tại chóp mũi quả lê ở dải tốc độ cao.
  3. Cố định bậc tự do mớn nước và độ chúi: Trong một số ca tính CFD, các biến động động học về độ lún ($Sinkage$) và độ chúi động ($Dynamic\ Trim$) của tàu chưa được giải ghép bộ hoàn toàn tự động ở từng bước thời gian.

Chương trình nghiên cứu 5-10 năm tới mở ra các hướng phát triển:

  • Ứng dụng CFD giải bài toán thủy động lực học tàu đi trong sóng biển thực tế (Seakeeping in rough seas) kết hợp đánh giá hiện tượng mất ổn định tham số.
  • Nghiên cứu tối ưu hóa đồng thời hình học mũi quả lê, thân tàu và hệ thống đẩy (chân vịt - đạo lưu - bánh lái) có xét đến ảnh hưởng của dòng theo không đều ($Wake\ fraction$).
  • Tích hợp mạng nơ-ron nhân tạo sâu (Deep Learning / Physics-Informed Neural Networks - PINNs) để thay thế mô hình Kriging, tăng tốc độ dự báo thủy động lực học theo thời gian thực ($Real-time\ prediction$).

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

  • Tác động học thuật: Kết quả nghiên cứu 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 ngành Kỹ thuật Tàu thủy, Động lực học công trình biển; mở ra hướng nghiên cứu ứng dụng CFD giải quyết các bài toán biên phức tạp của phương tiện thủy. Ước tính công trình có khả năng tạo ra bước nhảy vọt về số lượt trích dẫn trong các công bố chuyên ngành quốc tế.
  • Chuyển đổi ngành công nghiệp đóng tàu: Các nhà máy đóng tàu trong nước (như Tổng công ty Công nghiệp Tàu thủy SBIC, Nhà máy đóng tàu Ba Son, Nhà máy đóng tàu Sông Thu) có thể ứng dụng ngay bộ bản vẽ đường hình và thông số quả lê tối ưu để đóng mới hoặc cải hoán đội tàu dịch vụ hậu cần nghề cá, tàu kiểm ngư và tàu nghiên cứu biển.
  • Lợi ích xã hội và môi trường: Giảm thiểu hàng ngàn tấn khí thải carbon mỗi năm từ đội tàu cá hàng chục ngàn chiếc của Việt Nam, cải thiện điều kiện làm việc an toàn cho thuyền viên nhờ giảm chuyển động lắc dọc khi đánh bắt trên các ngư trường Hoàng Sa và Trường Sa.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận một khung phương pháp luận hoàn chỉnh từ mở rộng lý thuyết thực nghiệm, thiết lập mô phỏng số CFD đến giải thuật tối ưu hóa thay thế.
  • Kỹ sư thiết kế và Nhà máy đóng tàu: Sở hữu bộ cơ sở dữ liệu số hóa và tọa độ hình học chính xác của mẫu quả lê tối ưu cho tàu FAO 75, rút ngắn chu kỳ thiết kế từ vài tháng xuống vài ngày mà không cần chi phí thử bể tốn kém.
  • Ngư dân và Chủ tàu cá: Giảm trực tiếp chi phí vận hành chuyến biển, tăng tốc độ cơ động khi tiếp cận ngư trường và nâng cao khả năng giữ tàu ổn định trong điều kiện sóng gió cấp 5-cấp 6.
  • Cơ quan Quản lý Nhà nước: Có luận cứ khoa học để hoàn thiện hệ thống tiêu chuẩn kỹ thuật quốc gia (QCVN) về thiết kế tàu thủy tiết kiệm năng lượ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 phương pháp thiết kế đồ thị Kracht (1978) vốn bị giới hạn trong dải $C_B \in [0.56, 0.82]$ sang vùng hệ số béo thấp $C_B = 0.524$ của tàu cá vỏ thép. Luận án đã xác lập quy luật biến thiên phi tuyến của 6 hệ số hình học quả lê ($C_{LPR}, C_{BB}, C_{ZB}, C_{VPR}, C_{ABT}, C_{ABL}$), kết hợp hoàn hảo với lý thuyết triệt tiêu sóng của Inui (1962) để chứng minh rằng quả lê có thể tạo ra giao thoa sóng triệt tiêu hiệu quả ngay cả trên các thân tàu có tỷ số $L/B$ nhỏ và hình dáng vỏ tương đối béo ở phần mũi.

2. Điểm cải tiến phương pháp luận cốt lõi so với các nghiên cứu tiền nhiệm quốc tế là gì?

So với nghiên cứu của Filip et al. (2014) tại Đại học Michigan (chỉ thay đổi 2 kích thước hình học trên mô hình có sẵn quả lê của tàu container KCS) và Luo & Lan (2017) tại Đại học Phúc Châu (chỉ tối ưu sức cản sóng thế lưu trong ShipFlow), luận án của NCS Huỳnh Văn Chính đã tạo nên bước đột phá:

  • Xây dựng quy trình tự động hóa nối ghép bề mặt trơn nhẵn $C^1/G^1$ bằng đường cong NURBS trong AutoShip.
  • Áp dụng lời giải dòng nhớt RANS với mô hình SST $k-\omega$ có hiệu chuẩn tham số miền tính và hàm thành, đạt sai số kiểm chứng $< 5%$ so với bể thử SSPA/NPL.
  • Đổi mới hàm mục tiêu từ đơn biến sức cản tĩnh sang hàm đa mục tiêu giảm công suất có ích $\Delta P_e$ có trọng số thích ứng, phản ánh chính xác chu kỳ hoạt động thực tế của tàu cá.

3. Phát hiện bất ngờ và phản trực giác nhất từ dữ liệu nghiên cứu là gì?

Phát hiện bất ngờ nhất là: Quả lê tối ưu có kích thước hình học lớn hơn đáng kể so với quả lê ban đầu ngoại suy từ lý thuyết cổ điển (chiều dài tăng từ $1.50\text{ m}$ lên $1.65\text{ m}$, chiều rộng tăng từ $1.70\text{ m}$ lên $1.91\text{ m}$). Theo trực giác thông thường, quả lê lớn hơn sẽ làm tăng mạnh diện tích mặt ướt $S$ dẫn đến tăng sức cản ma sát $R_F$. Tuy nhiên, kết quả mô phỏng CFD chỉ ra rằng thể tích quả lê lớn hơn ở vị trí thích hợp đã tạo ra hiệu ứng triệt tiêu sóng mũi vượt bậc và dẫn dòng chảy mạn mượt mà, làm mức sụt giảm của sức cản sinh sóng $R_W$ và sức cản bẻ gãy sóng $R_{WB}$ vượt xa mức tăng của sức cản ma sát, mang lại mức giảm sức cản tổng kỷ lục lên tới $14%$.

4. Luận án có cung cấp đầy đủ giao thức tái lập (Replication Protocol) không?

Luận án cung cấp chi tiết tuyệt đối giao thức tái lập nghiên cứu:

  • Bảng tọa độ hình học 3D đầy đủ của tàu mẫu FAO 72 và FAO 75.
  • Kích thước không gian miền tính toán định chuẩn theo $L_{pp}$ (Inlet $1.5 L_{pp}$, Outlet $3.5 L_{pp}$, Bottom $2.5 L_{pp}$).
  • Giá trị chính xác của các hệ số mô hình rối SST $k-\omega$ ($\sigma_k, \sigma_\omega, \beta^*, \gamma$) và điều kiện biên vật lý.
  • Thuật toán ma trận phương án thay đổi kích thước quả lê và phương trình hồi quy bề mặt đáp ứng Kriging. Bất kỳ nhà nghiên cứu nào có trang bị phần mềm AutoShip và solver CFD tiêu chuẩn đều có thể tái lập chính xác 100% các kết quả tính toán.

5. Lộ trình nghiên cứu 10 năm tiếp theo được phác thảo như thế nào?

Lộ trình 10 năm tập trung vào 3 trọng tâm:

  1. Mở rộng nền tảng mô phỏng sang bài toán tương tác thủy - khí phi tuyến trong điều kiện sóng ngẫu nhiên và gió giật (Seakeeping CFD).
  2. Phát triển công cụ thiết kế tối ưu hóa đa trường vật lý (Multidisciplinary Design Optimization - MDO), kết hợp thủy động lực học, độ bền kết cấu vỏ và hiệu suất buồng máy.
  3. Số hóa và xây dựng hệ thống trí tuệ nhân tạo (AI-driven Generative Hull Design) tự động sinh đường hình quả lê tối ưu theo yêu cầu kỹ thuật của khách hàng chỉ trong vài phút.

Kết luận

Công trình luận án tiến sĩ của tác giả Huỳnh Văn Chính đã xác lập một dấu ấn khoa học xuất sắc với 5 đóng góp cốt lõi được lượng hóa chuẩn xác:

  1. Hoàn thiện phương pháp thiết kế đồ thị Kracht: Thiết lập thành công các đường cong nội suy và ngoại suy để tính toán thông số mũi quả lê ban đầu cho tàu cá FAO 75 có hệ số béo $C_B = 0.524$ nằm ngoài phạm vi truyền thống ($C_B \in [0.56, 0.82]$).
  2. Chuẩn hóa quy trình tính toán sức cản bằng CFD: Xác lập bộ thông số mô phỏng chuẩn (miền tính, điều kiện biên, mô hình rối SST $k-\omega$) cho tàu cá vỏ thép chạy chậm, kiểm chứng thành công với sai số $< 5%$ so với thực nghiệm bể thử quốc tế NPL và SSPA.
  3. Đổi mới mô hình toán học tối ưu: Xây dựng phương pháp tối ưu hóa hình học mũi quả lê dựa trên hàm đa mục tiêu giảm công suất có ích $\Delta P_e$ gắn liền với các chế độ vận hành thực tế của tàu cá thay vì hàm đơn mục tiêu truyền thống.
  4. Tìm ra cấu hình quả lê đột phá: Xác định chính xác kích thước quả lê tối ưu cho 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 đến 14% sức cản tổng của tàu.
  5. Mở ra hướng nghiên cứu ứng dụng công nghệ cao: Khẳng định vai trò của CFD như một "phương pháp thứ ba" tất yếu, kết hợp chặt chẽ giữa lý thuyết giải tích và thực nghiệm, đặt nền móng vững chắc cho công cuộc hiện đại hóa ngành thiết kế và đóng tàu thủy tại Việt Nam.