Tổng quan về luận án

Trong bối cảnh ngành công nghiệp đóng tàu hiện đại đang chuyển dịch mạnh mẽ hướng tới tiêu chuẩn giảm phát thải khí nhà kính và tối ưu hóa hiệu suất năng lượng theo quy định của Tổ chức Hàng hải Quốc tế (IMO), việc hoàn thiện đường hình thân vỏ tàu thủy trở thành một nhiệm vụ then chốt. Tại Việt Nam, thực hiện chủ trương hiện đại hóa đội tàu đánh bắt xa bờ theo Nghị định 67/2014/NĐ-CP của Chính phủ, nhu cầu thiết kế các mẫu tàu cá vỏ thép hoạt động tại các ngư trường xa bờ với chiều dài trên 40 m ngày càng cấp thiết. Tuy nhiên, phần lớn các tàu cá trong nước trước đây đều thiết kế dạng mũi thẳng truyền thống do hạn chế về kinh phí thử nghiệm bể thử và thiếu dữ liệu thủy động lực học chuyên sâu, làm suy giảm đáng kể hiệu quả kinh tế và khả năng hành hải trong điều kiện sóng gió.

Luận án tiến sĩ kỹ thuật cơ khí động lực mang tên: "Ứ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 khoa học 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 TP. Hồ Chí Minh đã giải quyết trọn vẹn khoảng trống khoa học này. Nghiên cứu tập trung vào bài toán nâng cao hiệu suất thủy động học cho 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 72 và FAO 75) vốn đã được thử nghiệm tại các bể thử hàng đầu thế giới như NPL (National Physical Laboratory - Anh) và SSPA (Swedish State Shipbuilding Tank - Thụy Điển).

Khoảng trống nghiên cứu (Research Gap) cốt lõi được xác định gồm ba phương diện:

  1. Phương pháp đồ thị kinh điển của M. Kracht (1978) chỉ áp dụng cho các tàu hàng có hệ số béo thể tích $C_B$ nằm trong giới hạn hẹp $C_B \in [0.56, 0.82]$, hoàn toàn thiếu vắng công cụ thiết kế giải tích cho các dòng tàu cá thân béo thấp ($C_B = 0.524$) với tỷ số kích thước $L/B$ nhỏ và đoạn thân ống ngắn.
  2. Các nghiên cứu quốc tế đương đại thường tối ưu hóa cục bộ trên thân tàu hàng đã có sẵn quả lê thông qua hàm đơn mục tiêu tại một vận tốc thiết kế cố định, không phản ánh đúng chu kỳ khai thác thực tế của tàu cá vốn hoạt động ở nhiều chế độ vận tốc và mớn nước thay đổi liên tục.
  3. Thiếu một quy trình tích hợp khép kín từ mô hình hóa tham số 3D, mô phỏng số động lực học lưu chất (CFD) có kiểm soát độ nhạy tham số, đến giải thuật tối ưu hóa bề mặt đáp ứng (Surrogate Model) đảm bảo độ trơn nhẵn bề mặt vỏ tàu.

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

  • RQ1: Làm thế nào để mở rộng phạm vi áp dụng của phương pháp đồ thị Kracht cho các mẫu tàu cá có hệ số béo $C_B < 0.56$ nhằm thiết lập thông số hình học quả lê ban đầu?
  • RQ2: Quy chuẩn thiết lập miền tính toán và lựa chọn mô hình rối trong CFD như thế nào để sai số dự báo sức cản so với thực nghiệm bể thử kiểm soát dưới 5%?
  • 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 nhất hiệu quả giảm công suất có ích $P_e$ qua các chế độ vận tốc khai thác của tàu cá?
  • H1: Việc xây dựng các đường cong nội suy và ngoại suy toán học từ dữ liệu Kracht cho phép xác định tập thông số quả lê khả dĩ tiệm cận tối ưu cho tàu cá $C_B = 0.524$.
  • H2: Mô hình chảy rối SST $k-\omega$ kết hợp cấu trúc miền tính toán chuẩn hóa trong phần mềm CFD (XFlow) sẽ tái hiện chính xác trường áp suất và hiện tượng giao thoa sóng mũi.
  • H3: Tối ưu hóa đồng thời các kích thước hình học chính của quả lê thông qua mô hình thay thế (Surrogate Model) hai giai đoạn sẽ giảm trên 10% sức cản tổng của tàu.

Khung lý thuyết của luận án được xây dựng dựa trên sự tích hợp của lý thuyết triệt tiêu sóng cục bộ Inui (1962), định luật phân rã sức cản Froude, hệ phương trình trung bình Reynolds Navier-Stokes (RANSE) đóng kín bằng mô hình rối hai phương trình SST $k-\omega$, và phương pháp bề mặt đáp ứng (Response Surface Methodology - RSM). Kết quả đột phá của luận án đã chứng minh rằng: "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}$... và độ giảm sức cản tổng của tàu bằng khoảng 14%". Toàn bộ nghiên cứu được thực nghiệm số trên mẫu tàu FAO 72 ($L_{pp} = 22.0\text{ m}$) và FAO 75 ($L_{pp} = 27.5\text{ m}$, $C_B = 0.524$) tại các dải số Froude $Fn = 0.15 - 0.40$, đóng góp cơ sở dữ liệu học thuật chuẩn tắc cho ngành đóng tàu Việt Nam.

Literature Review và Positioning

Lịch sử khí thủy động lực học tàu thủy ghi nhận sự phát triển của mũi quả lê từ công trình tiên phong của D. Taylor (1910) trên tàu hải quân USS Delaware, tiếp nối qua các thế hệ tàu chở khách Bremen, Europa (1920) và Normandie (1935). Về mặt lý thuyết hàn lâm, bước ngoặt lớn diễn ra vào năm 1962 khi Giáo sư T. Inui công bố lý thuyết dạng tàu không sóng (waveless form theory). Theo Inui (1962), hệ thống sóng do mũi và đuôi tàu sinh ra có thể được triệt tiêu bằng cách bố trí các phân bố nguồn âm - dương (điểm kỳ dị toán học) tại vị trí quả lê nhằm 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 chính. Ứng dụng thực nghiệm của Inui trên các tàu Murasaki và Kurenai Maru (Nhật Bản) với đường kính quả lê gấp 3 lần bình thường đã minh chứng khả năng giảm đến 13.1% công suất máy chính.

Công trình thực nghiệm có quy mô toàn diện và phổ quát nhất thuộc về M. Kracht (1978). Tác giả đã phân tích dữ liệu kéo mô hình của hàng trăm dạng mũi tàu trong bể thử để thiết lập hệ thống đồ thị định lượng 6 hệ số hình học quả lê đặc trưng: hệ số chiều dài $C_{LPR}$, hệ số chiều rộng $C_{BB}$, hệ số chiều cao $C_{ZB}$, hệ số diện tích mặt cắt ngang $C_{ABT}$, hệ số diện tích cạnh bên $C_{ABL}$, và hệ số thể tích mũi quả lê $C_{\nabla PR}$. Mặc dù vậy, đồ thị Kracht tồn tại hạn chế nội tại: chỉ nội suy trong dải hệ số béo $C_B = 0.56 - 0.82$, không đưa ra thuật toán hình học để nối ghép trơn tru bề mặt quả lê vào tuyến hình thân tàu thực tế.

Trong kỷ nguyên số, các nghiên cứu quốc tế phân hóa thành hai trường phái chính:

  1. Trường phái tối ưu hóa theo phương án (Variant-based Optimization): Tiêu biểu là công trình 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 mẫu tàu container KCS (KRISO Container Ship). Nhóm tác giả thiết lập 9 phương án biến đổi kích thước chiều cao $H$ và chiều rộng $C$, sử dụng bộ giải mã nguồn mở OpenFOAM kết hợp mô hình thay thế Kriging để tìm điểm cực trị. Tuy nhiên, hạn chế của nghiên cứu này là việc lựa chọn 9 phương án ban đầu còn mang tính rời rạc, chưa có tiêu chí ràng buộc thủy tĩnh học chặt chẽ.
  2. Trường phái mô hình hóa tham số (Parametric Modeling Optimization): Đại diện là Weilin Luo và Linqiang Lan (2013, Đại học Phúc Châu) áp dụng môi trường tích hợp CAESES và ShipFlow trên tàu Ro-Ro ($L = 165.7\text{ m}, C_B = 0.65, Fn = 0.28$), sử dụng các đường cong B-Spline và mặt cong NURBS để tự động biến dạng mũi tàu. Hướng tiếp cận tương tự được triển khai bởi các nhà nghiên cứu của tổ chức Joint High Speed Sealift (JHSS, Hoa Kỳ) khi áp dụng hàm cơ sở xuyên tâm (RBF) để tối ưu hóa 4 biến thể hình học mũi: ST (Stem Bow), EB (Elliptical Bulb), GB (Gooseneck Bulb), và BB (Baseline Bulb). Ngoài ra, các nghiên cứu ứng dụng giải thuật di truyền (GA) trên tàu Series 60 (2012) cũng chứng minh hiệu quả giảm sức cản sóng $R_W$ nhưng đòi hỏi chi phí tính toán cực lớn.

Tại Việt Nam, các nghiên cứu CFD giai đoạn đầu của Đại học Bách Khoa Hà Nội (sử dụng phương pháp điểm kỳ dị cho tàu Wigley) và Đại học Nha Trang (luận án của Lê Văn Toàn, 2015 mô phỏng OpenFOAM cho tàu cá vỏ gỗ; ứng dụng biến đổi Lackenby) chủ yếu tập trung vào thân tàu trần. Nghiên cứu của PGS.TS Đỗ Quang Khải (Đại học Hàng hải Việt Nam / Đại học Liège) đã mô hình hóa đường nước tàu bằng hàm đa thức nhưng dừng lại ở mức độ biểu diễn hình học giải tích.

Luận án của NCS. Huỳnh Văn Chính định vị chính xác tại giao điểm của các dòng nghiên cứu trên: vừa khắc phục giới hạn của đồ thị thực nghiệm Kracht (1978) cho tàu cá thân béo thấp ($C_B = 0.524$), vừa chuẩn hóa công cụ CFD tiên tiến để thay thế thử nghiệm mô hình tốn kém, đồng thời giải quyết bài toán tối ưu hóa đa mục tiêu phù hợp với đặc thù khai thác thực tiễn của ngành thủy sản.

Đó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 đóng góp nền tảng cho lý thuyết thủy động lực học tàu thủy và thiết kế hình học thân tàu:

  • Mở rộng lý thuyết thiết kế Kracht: Thiết lập thành công các hàm toán học nội suy và ngoại suy, mở rộng đường cong thực nghiệm của Kracht sang miền giá trị hệ số béo thấp $C_B < 0.56$. Đây là đóng góp lý thuyết quan trọng giúp tính toán giải tích xác định bộ thông số quả lê ban đầu ($C_{LPR}, C_{BB}, C_{ZB}, C_{\nabla PR}$) cho các lớp tàu chuyên dụng như tàu cá xa bờ, tàu tuần tra cao tốc.
  • Mô hình toán học tối ưu hóa đa mục tiêu đặc thù: Xây dựng hàm mục tiêu tổng quát dựa trên độ giảm công suất có ích hiệu dụng $\Delta P_e$: $$\Delta P_e = \sum_{i=1}^{n} w_i \cdot \Delta P_{ei} = \sum_{i=1}^{n} w_i \frac{P_{ei} - P_{ebi}}{P_{ei}} \times 100%$$ trong đó $w_i$ là trọng số thời gian hoạt động tương ứng với 3 chế độ công tác chính của tàu cá ($U_1$: chạy hành trình, $U_2$: chạy dắt lưới, $U_3$: chạy kéo/thả lưới). Mô hình này thay thế hoàn toàn cách tiếp cận đơn mục tiêu truyền thống, đảm bảo hình dạng quả lê tối ưu đồng thời ở nhiều trạng thái tải trọng và mớn nước khai thác ($T = 4.57\text{ m}$).
  • Làm sáng tỏ bản chất vi mô của tương tác trường dòng: Bằng mô phỏng phân giải cao, luận án chứng minh cơ chế vật lý của việc triệt tiêu sóng: quả lê tối ưu tạo ra vùng hạ áp cục bộ ngay phía trước sống mũi, làm giảm gradient áp suất dương tại điểm dừng (stagnation point), triệt tiêu búp sóng mũi thứ nhất và tái phân bố trường vận tốc nhớt, hạn chế tối đa hiện tượng phân kỳ dòng chảy sinh xoáy.

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

Khung phân tích của luận án tích hợp một quy trình 4 giai đoạn logic và chặt chẽ:

  1. Giai đoạn 1 - Thiết kế sơ bộ: Mở rộng đồ thị Kracht qua giải tích toán học để xác định các kích thước quả lê ban đầu ($L_{PRo} = 1.50\text{ m}, B_{Bo} = 1.70\text{ m}, Z_{Bo} = 2.10\text{ m}$).
  2. Giai đoạn 2 - Mô hình hóa hình học NURBS: Sử dụng phần mềm chuyên dụng AutoShip để dựng bề mặt 3D, tái cấu trúc các đường sườn và đường nước, đảm bảo điều kiện biên tiếp giáp bậc cao ($G^1/G^2$ continuity) giữa quả lê và thân tàu vỏ dưa, bảo toàn tính trơn nhẵn thủy động học.
  3. Giai đoạn 3 - Động lực học lưu chất tính toán (CFD): Áp dụng phần mềm CFD tiên tiến XFlow giải hệ phương trình Navier-Stokes, thiết lập điều kiện biên không gian miền tính tối ưu, thu nhận dữ liệu sức cản tổng $R_T$, sức cản áp suất $R_P$ và sức cản ma sát $R_F$.
  4. Giai đoạn 4 - Tối ưu hóa bề mặt đáp ứng hai vòng lặp: Xây dựng ma trận thực nghiệm với các biến thiết kế kích thước quả lê ($\Delta L_{PR}, \Delta B_B, \Delta Z_B$), xấp xỉ hàm mục tiêu bằng mô hình bề mặt đáp ứng (RSM/Kriging) để quét toàn bộ không gian nghiệm và hội tụ tại cấu hình hình học tối ưu toàn cục.
+-------------------------------------------------------------------------+
|                  KHUNG PHÂN TÍCH TỐI ƯU HÓA QUẢ LÊ                      |
+-------------------------------------------------------------------------+
| [1. Thiết kế Kracht mở rộng] --> Tính L_PRo, B_Bo, Z_Bo (C_B = 0.524)   |
|                                     |                                   |
| [2. Tạo hình AutoShip NURBS]  --> Tích hợp trơn bề mặt vỏ dưa tàu cá     |
|                                     |                                   |
| [3. Mô phỏng số CFD - XFlow]  --> Giải RANSE, SST k-omega (Sai số < 5%) |
|                                     |                                   |
| [4. Mô hình thay thế RSM]     --> Quét ma trận & Tối ưu đa mục tiêu P_e |
+-------------------------------------------------------------------------+

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

Thiết kế nghiên cứu

Nghiên cứu tuân thủ chặt chẽ thế giới quan thực chứng (Positivism) kết hợp phương pháp luận diễn dịch - kiểm chứng số. Mô hình toán học mô phỏng dòng chảy quanh thân tàu dựa trên hệ phương trình trung bình Reynolds Navier-Stokes (RANSE) cho chất lưu nhớt không nén được: $$\frac{\partial \overline{u}_i}{\partial x_i} = 0$$ $$\rho \left( \frac{\partial \overline{u}_i}{\partial t} + \overline{u}_j \frac{\partial \overline{u}_i}{\partial x_j} \right) = -\frac{\partial \overline{p}}{\partial x_i} + \frac{\partial}{\partial x_j} \left[ \mu \left( \frac{\partial \overline{u}_i}{\partial x_j} + \frac{\partial \overline{u}_j}{\partial x_i} \right) - \rho \overline{u'_i u'_j} \right]$$ trong đó $-\rho \overline{u'_i u'_j}$ là tensor ứng suất Reynolds được đóng kín bằng mô hình rối hai phương trình SST $k-\omega$ (Shear Stress Transport) của Menter. Mô hình SST $k-\omega$ kết hợp ưu điểm của mô hình $k-\omega$ chuẩn ở vùng sát biên vỏ tàu và mô hình $k-\epsilon$ ở vùng dòng tự do xa vỏ tàu, cho phép dự báo chính xác hiện tượng bóc tách dòng chảy do gradient áp suất ngược tại khu vực mũi và đuôi tàu cá.

Quy trình nghiên cứu rigorous

Quy trình mô phỏng số và kiểm chứng được tiến hành chuẩn mực theo hướng dẫn của Hội nghị Quốc tế về Bể thử Tàu thủy (ITTC):

  1. Thiết lập miền tính toán (Domain Independence Test): Khảo sát độc lập kích thước miền tính toán với các phương án chiều dài, chiều rộng và chiều sâu khác nhau. Cấu hình miền tính tối ưu được lựa chọn đảm bảo loại bỏ hoàn toàn ảnh hưởng phản xạ sóng biên: khoảng cách từ mũi tàu đến biên vào (Inlet) là $1.5 L_{pp}$; từ đuôi tàu đến biên ra (Outlet) là $3.0 L_{pp}$; khoảng cách sang hai mạn (Sides) là $1.5 L_{pp}$; khoảng cách từ đáy tàu xuống đáy miền tính (Bottom) là $1.5 L_{pp}$.
  2. Xác lập điều kiện biên: Biên vào gán vận tốc dòng chảy $U_\infty$; biên ra gán áp suất thủy tĩnh $p = p_0$; bề mặt vỏ tàu áp dụng điều kiện biên không trượt (no-slip wall condition); mặt thoáng tự do được giải quyết chính xác bằng thuật toán nắm bắt mặt phân giới chất lỏng.
  3. Phân tích độ nhạy lưới và thông số mô hình rối: Đánh giá các hệ số suy giảm động năng rối $\beta^*, \sigma_k, \sigma_\omega$ nhằm tối ưu hóa độ phân giải lớp biên thủy động, đảm bảo bắt trọn các cấu trúc vi xoáy mà không gây phân kỳ nghiệm số.

Data và phân tích

Đối tượng dữ liệu thực nghiệm gồm 2 mẫu tàu cá tiêu chuẩn:

  • Tàu FAO 72: $L_{pp} = 22.00\text{ m}, B = 6.60\text{ m}, d = 2.70\text{ m}, C_B = 0.443$, góc vào nước $\alpha_E = 27.5^\circ$.
  • Tàu FAO 75: $L_{pp} = 27.50\text{ m}, B = 7.60\text{ m}, d = 2.90\text{ m}, C_B = 0.524, \nabla = 317.0\text{ m}^3$, góc vào nước $\alpha_E = 29.0^\circ$.

Dữ liệu tính toán từ phần mềm XFlow được đối chiếu trực tiếp với dữ liệu thực nghiệm kéo mô hình bể thử NPL và SSPA. Mức sai lệch tương đối $\delta R = \frac{|R_{XF} - R_{tn}|}{R_{tn}} \times 100%$ tại các dải vận tốc khảo sát đều đạt dưới 5%, khẳng định độ tin cậy và tính hợp lệ tuyệt đối của mô hình CFD.

Ma trận thực nghiệm tối ưu hóa được xây dựng với các biến hình học:

  • Chiều dài quả lê $L_{PR}$: thay đổi với bước nhảy $\Delta L_{PR} = 0.15\text{ m}$ (từ $1.35\text{ m}$ đến $1.80\text{ m}$).
  • Chiều rộng quả lê $B_B$: thay đổi với bước nhảy $\Delta B_B = 0.21\text{ m}$ (từ $1.49\text{ m}$ đến $2.12\text{ m}$).
  • Chiều cao tâm quả lê $Z_B$: thay đổi từ $1.89\text{ m}$ đến $2.31\text{ m}$.

Mô hình thay thế đa thức bậc 2 được hiệu chỉnh qua hai vòng lặp lồng ghép để xác định điểm hội tụ tối ưu cực đại của hàm $\Delta P_e$.

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

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

Nghiên cứu mang lại 5 phát hiện đột phá có giá trị thực chứng sâu sắc:

  1. Thiết kế hình học quả lê tối ưu toàn cục: Đối với tàu cá FAO 75, hình dạng quả lê tối ưu đạt được tại các kích thước: chiều dài $L_{PRop} = 1.65\text{ m}$ (tăng $0.15\text{ m}$ so với quả lê sơ bộ $L_{PRo} = 1.50\text{ m}$), chiều rộng lớn nhất $B_{Bop} = 1.91\text{ m}$ (tăng $0.21\text{ m}$ so với $B_{Bo} = 1.70\text{ m}$), chiều cao $Z_{Bop} = 2.10\text{ m}$.
  2. Hiệu quả giảm sức cản tổng vượt bậc: Quả lê tối ưu tạo ra hiện tượng giao thoa sóng ngược pha hoàn hảo tại vùng vận tốc khai thác danh định ($Fn = 0.35 - 0.38$), giúp giảm tới 14% sức cản tổng $R_T$ của tàu so với nguyên bản không có quả lê.
  3. Cơ chế cân bằng động giữa ma sát và sóng: Mặc dù việc lắp đặt quả lê làm tăng diện tích mặt ướt $S$ từ đó làm tăng nhẹ sức cản ma sát $R_F$, nhưng mức độ triệt tiêu của thành phần sức cản sinh sóng $R_W$ và sức cản bẻ gãy sóng mũi $R_{WB}$ chiếm ưu thế áp đảo, dẫn đến công suất kéo có ích $P_e$ giảm thuần trên toàn bộ các chế độ làm việc.
  4. Ảnh hưởng quyết định của chất lượng tiếp giáp bề mặt: Kết quả CFD chỉ ra rằng nếu quả lê có kích thước tối ưu nhưng không được xử lý đường biên chuyển tiếp trơn nhẵn với thân tàu (tạo góc gấp hoặc đột biến độ cong), dòng chảy sẽ bị tách lớp ngay sau quả lê, tạo ra vùng xoáy áp suất nhớt dư $R_{PV}$ làm triệt tiêu hoàn toàn lợi ích giảm lực cản sóng.
  5. Cải thiện tính năng lắc dọc và chúi mũi: Tàu trang bị mũi quả lê tối ưu có phân bố thể tích phần mũi đầy đặn hơn, làm tăng mômen quán tính diện tích đường nước mũi, giúp cải thiện rõ rệt độ ổn định lắc dọc (pitching) và độ chúi dọc (trim), tạo điều kiện an toàn vượt trội cho thao tác thả và thu lưới trên biển.

Implications đa chiều

  • Về mặt học thuật và phương pháp luận: Khẳng định tính khả thi của việc kết hợp CFD không lưới/phần tử hạt hoặc thể tích hữu hạn (FVM) với mô hình thay thế (Surrogate Model) trong thủy động lực học tàu thủy, mở ra phương pháp luận thiết kế tối ưu hình học tin cậy mà không phụ thuộc vào các bể thử mô hình đắt đỏ.
  • Về mặt kỹ thuật thực tiễn: Cung cấp bộ bản vẽ đường hình và quy chuẩn kích thước hình học chính xác để các viện thiết kế, nhà máy đóng tàu tại Việt Nam ứng dụng đóng mới hoặc cải hoán mũi quả lê cho các đội tàu cá xa bờ chiều dài $24 - 45\text{ m}$.
  • Về mặt chính sách và kinh tế biển: Việc giảm 14% lực cản đồng nghĩa với việc tiết kiệm từ 10 - 12% lượng dầu diesel tiêu thụ trên mỗi chuyến biển dài ngày. Điều này trực tiếp hạ giá thành khai thác, nâng cao năng lực cạnh tranh của ngư dân, đồng thời hiện thực hóa cam kết quốc gia về giảm phát thải CO₂ trong hoạt động thủy hải sản.

Limitations và Future Research

Nhằm duy trì tính khách quan và chuẩn mực học thuật, luận án thẳng thắn chỉ ra các giới hạn nghiên cứu:

  1. Giới hạn môi trường mô phỏng: Các phân tích CFD được thực hiện trong môi trường nước tĩnh, độ sâu không hạn chế, chưa mô phỏng tương tác động đa trường trong phổ sóng ngẫu nhiên (Seakeeping in irregular waves) và tác động của gió phức hợp.
  2. Giới hạn tương tác hệ động lực: Nghiên cứu tập trung vào sức cản thân tàu kéo theo định luật Archimede, chưa tích hợp mô phỏng tự hành (self-propulsion) có sự làm việc đồng thời của chân vịt và bánh lái phía đuôi để đánh giá sự biến thiên của hệ số dòng theo $w$ và hệ số hút $t$.
  3. Phạm vi hình học quả lê: Luận án khảo sát biến thể quả lê dạng giọt nước/ellip biến tính, chưa mở rộng phân tích các cấu trúc mũi phức tạp như mũi ngược (X-bow) hay mũi dao (axe bow).

Các hướng nghiên cứu tương lai được vạch ra rõ nét:

  • Ứng dụng mô phỏng CFD 6 bậc tự do (6-DoF) khảo sát tính năng hàng hải của tàu cá có mũi quả lê trong điều kiện sóng biển cấp 5 - cấp 7 theo thang Beaufort.
  • Xây dựng hệ thống tối ưu hóa đa ngành đồng thời (Multidisciplinary Design Optimization - MDO) tích hợp giữa thủy động lực học, độ bền kết cấu tấm vỏ mũi tàu khi chịu va đập sóng (slamming), và bài toán kinh tế năng lượng.
  • Tiến hành đo đạc thực nghiệm trường vận tốc bằng kỹ thuật PIV (Particle Image Velocimetry) trên mô hình tự hành để hoàn thiện bộ dữ liệu benchmark quốc tế.

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

Công trình nghiên cứu tạo ra những tác động lan tỏa sâu rộng:

  • Tác động học thuật: Đóng góp nguồn tài liệu tham khảo chuyên sâu và dữ liệu số hóa có giá trị trích dẫn cao trong các tạp chí chuyên ngành cơ khí thủy khí, kỹ thuật hàng hải và đóng tàu.
  • Chuyển giao công nghệ ngành đóng tàu: Cung cấp giải pháp công nghệ trực tiếp cho các nhà máy đóng tàu trong nước (như Nhà máy đóng tàu Ba Son, Công ty Đóng tàu Hồng Hà, Viện Khoa học Công nghệ Tàu thủy) để chuẩn hóa thiết kế mũi quả lê tiết kiệm năng lượng.
  • Hiệu quả kinh tế xã hội: Đối với một tàu cá công suất 800 - 1000 HP hoạt động 250 ngày biển/năm, mức tiết kiệm 14% sức cản giúp cắt giảm hàng chục nghìn lít dầu DO, tiết kiệm hàng trăm triệu đồng chi phí vận hành mỗi năm trên từng con tàu.

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

  • Nghiên cứu sinh và Giới học thuật: Tiếp cận phương pháp luận mở rộng lý thuyết Kracht và quy trình kiểm soát độ chính xác mô phỏng CFD chuẩn ITTC.
  • Kỹ sư và Chuyên gia R&D ngành đóng tàu: Sở hữu công cụ và bảng thông số tối ưu để ứng dụng trực tiếp vào các dự án thiết kế tàu cá vỏ thép, tàu dịch vụ hậu cần nghề cá.
  • Ngư dân và Doanh nghiệp khai thác hải sản: Vận hành những con tàu có tính năng hàng hải ưu việt, độ lắc êm hơn, tiết kiệm chi phí nhiên liệu và nâng cao độ an toàn sinh mạng trên biển.
  • Cơ quan quản lý nhà nước (Cục Hàng hải, Tổng cục Thủy sản): Có cơ sở khoa học xác đáng để ban hành các tiêu chuẩn kỹ thuật, định mức tiêu hao năng lượng và phê duyệt các mẫu tàu đóng mới theo chiến lược phát triển kinh tế biển bền vữ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ì? Đó là việc mở rộng biên độ ứng dụng của hệ thống đồ thị Kracht (1978) sang vùng hệ số béo thấp $C_B = 0.524$ bằng phương pháp ngoại suy giải tích toán học, phá vỡ giới hạn truyền thống $C_B \in [0.56, 0.82]$ vốn tồn tại hơn 4 thập kỷ trong ngành tàu thủy.

  2. Đột phá về phương pháp luận so với các nghiên cứu tiền nhiệm? So với phương pháp thử nghiệm bể thử đơn thuần tốn kém của Kracht (1978) hay tối ưu hóa 9 phương án rời rạc của Grzegorz Filip (2014), luận án đã thiết lập quy trình tự động hóa khép kín: Ngoại suy Kracht $\rightarrow$ Làm trơn NURBS trên AutoShip $\rightarrow$ Giải RANSE SST $k-\omega$ trên XFlow (sai số $<5%$) $\rightarrow$ Quét cực trị liên tục bằng mô hình thay thế (Surrogate Model) đa mục tiêu.

  3. Phát hiện bất ngờ và thú vị nhất về mặt số liệu? Mặc dù quả lê tối ưu có chiều dài và chiều rộng lớn hơn đáng kể so với kích thước ban đầu ($L_{PRop} = 1.65\text{ m}$ so với $1.50\text{ m}$; $B_{Bop} = 1.91\text{ m}$ so với $1.70\text{ m}$), làm tăng diện tích mặt ướt cản ma sát, nhưng hiệu ứng triệt tiêu sóng cục bộ và tái phân bố áp suất lại đạt đỉnh cực đại, giúp giảm đột biến tới 14% tổng lực cản toàn tàu.

  4. Luận án có cung cấp quy trình tái lập (Replication Protocol) không? Có. Quy trình tái lập được mô tả chi tiết từ tọa độ điểm đường hình tàu FAO 75, kích thước hình học không gian miền tính toán ($1.5L_{pp} \times 1.5L_{pp} \times 3.0L_{pp}$), các tham số đóng kín mô hình rối SST $k-\omega$, đến các bước thiết lập bề mặt đáp ứng trong mô hình thay thế.

  5. Chương trình nghiên cứu dài hạn trong 10 năm tới được vạch ra như thế nào? Định hướng phát triển tập trung vào: (1) Tích hợp CFD với trí tuệ nhân tạo (AI/Machine Learning) để tối ưu hóa hình dáng vỏ tàu theo thời gian thực; (2) Mở rộng bài toán Seakeeping đa trường tương tác sóng gió thực tế; (3) Chuyển đổi xanh cho đội tàu thủy nội địa và tàu vận tải ven biển Việt Nam.

Kết luận

Luận án tiến sĩ của Huỳnh Văn Chính là một công trình nghiên cứu công phu, nghiêm túc và có hàm lượng khoa học cao, tạo ra những dấu ấn đột phá:

  1. Hoàn thiện và mở rộng xuất sắc lý thuyết thiết kế mũi quả lê Kracht cho phân khúc tàu cá có hệ số béo thấp $C_B = 0.524$.
  2. Chuẩn hóa quy trình ứng dụng CFD giải hệ phương trình RANSE với mô hình rối SST $k-\omega$, đảm bảo độ tin cậy tuyệt đối với sai số dưới 5% so với thực nghiệm bể thử quốc tế NPL và SSPA.
  3. Thiết lập thành công mô hình toán học tối ưu hóa đa mục tiêu theo chu kỳ công tác thực tế của tàu cá, kết hợp hiệu quả với phương pháp bề mặt đáp ứng (Surrogate Model).
  4. Xác định cấu hình quả lê tối ưu toàn cục cho tàu FAO 75 ($L_{PRop} = 1.65\text{ m}, B_{Bop} = 1.91\text{ m}, Z_{Bop} = 2.10\text{ m}$), cắt giảm thành công khoảng 14% sức cản tổng.
  5. Giải quyết trọn vẹn bài toán tích hợp trơn nhẵn bề mặt quả lê vào thân tàu vỏ dưa trên nền tảng AutoShip, loại bỏ hiện tượng sinh xoáy cục bộ.
  6. Đóng góp trực tiếp vào mục tiêu hiện đại hóa đội tàu cá vỏ thép Việt Nam, mang lại giá trị kinh tế, môi trường và ý nghĩa chiến lược đối với chủ quyền biển đảo quốc gia.