Tổng quan nghiên cứu

Mô phỏng động lực học chất lưu trên mặt biển thời gian thực là một trong những bài toán phức tạp hàng đầu của khoa học tính toán hiện đại. Bài toán này đòi hỏi phải giải hệ phương trình Navier - Stokes phi tuyến với khối lượng tính toán khổng lồ lên tới hàng triệu phép tính ma trận trong mỗi bước thời gian. Trên các hệ thống máy tính tuần tự truyền thống (kiến trúc Von Neumann), việc giải các phương trình vi phân đạo hàm riêng rời rạc hóa thường gặp nút thắt cổ chai về băng thông dữ liệu, dẫn đến độ trễ xử lý kéo dài từ vài giây đến hàng chục phút. Sự chậm trễ này làm mất đi tính kịp thời trong các ứng dụng khẩn cấp như dự báo sóng thần, kiểm soát dòng chảy hải văn hay điều hướng tàu thuyền ven biển.

Nhằm khắc phục triệt để hạn chế trên, nghiên cứu tập trung vào giải pháp ứng dụng công nghệ Mạng nơron tế bào (Cellular Neural Network - CNN) kết hợp cấu hình phần cứng chuyên dụng FPGA để giải hệ phương trình Navier - Stokes cho dòng chảy không nén hai chiều trên mặt biển. Mục tiêu cụ thể của luận văn là xây dựng mô hình toán học tương đương, thiết lập thuật toán sai phân thích hợp, thiết kế các mẫu liên kết vô tính (cloning templates) cho mạng CNN và tiến hành mô phỏng kiểm chứng trên môi trường MATLAB. Phạm vi nghiên cứu được giới hạn trong không gian mặt biển 2 chiều với lưới sai phân chuẩn hóa tại Trường Đại học Công nghệ Thông tin và Truyền thông – Đại học Thái Nguyên. Kết quả nghiên cứu mang ý nghĩa đột phá khi khai thác kiến trúc tính toán song song vật lý tương tự (analogic), cho phép đạt tốc độ tính toán lý thuyết lên tới $10^{12}$ phép tính/giây và thời gian quá độ mạch điện ở mức $1{,}5\text{ }\mu\text{s}$, cắt giảm hơn 85% độ trễ tính toán so với các phương pháp xử lý số tuần tự trên máy tính cá nhân thông thường.

Cơ sở lý thuyết và phương pháp nghiên cứu

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên sự giao thoa giữa lý thuyết cơ học thủy khí hiện đại và lý thuyết mạch phi tuyến tính toán song song. Khung lý thuyết cốt lõi bao gồm hai trụ cột chính: Lý thuyết Mạng nơron tế bào (CNN) do Leon O. Chua và Lin Yang phát minh năm 1988 cùng Kiến trúc máy tính vạn năng CNN-UM (Universal Machine); và Lý thuyết phương trình vi phân đạo hàm riêng (PDE) phi tuyến mô tả dòng chảy thủy lực với trọng tâm là hệ phương trình Navier - Stokes hai chiều.

Hệ thống khái niệm chuyên ngành then chốt bao gồm:

  • Tế bào nơron chuẩn (Cell): Khối mạch điện tích hợp phi tuyến cục bộ gồm các phần tử điện trở, tụ điện tuyến tính và nguồn dòng phụ thuộc, hoạt động như một bộ xử lý vi phân tương tự độc lập.
  • Bán kính tương tác lân cận ($r=1$): Quy tắc ghép nối không gian giới hạn trong phạm vi $3 \times 3$ tế bào liền kề, giúp giảm thiểu độ phức tạp dây dẫn phần cứng nhưng vẫn bảo đảm khả năng lan truyền thông tin toàn cục.
  • Mẫu liên kết vô tính (Cloning Templates): Bộ tham số trọng số tương tác đồng nhất gồm mẫu hồi tiếp $A$, mẫu điều khiển $B$ và ngưỡng dịch $z$, được đặc tả bởi 19 tham số thực xác định đặc tính động học của mạng.
  • Phương trình sai phân vi phân tế bào (CPDDE): Cầu nối toán học chứng minh tính tương đương động học giữa phương trình đạo hàm riêng liên tục và mạng mạch tế bào rời rạc trong không gian.
  • Chuẩn phi tuyến hàm bão hòa: Hàm kích hoạt đầu ra giới hạn nghiêm ngặt trong khoảng biên độ $[-1, 1]$, bảo đảm tính ổn định tiệm cận của toàn bộ hệ thống động lực.

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

Luận văn kết hợp phương pháp phân tích giải tích, phương pháp sai phân hữu hạn (Finite Difference Method) và mô phỏng thực nghiệm trên máy tính. Dữ liệu đầu vào sử dụng các tham số thủy động lực học chuẩn hóa của môi trường nước biển bao gồm trường vận tốc dòng chảy hai chiều theo trục tọa độ và độ cao mực nước tương đối.

Miền không gian liên tục được rời rạc hóa thành lưới sai phân chữ nhật 2 chiều gồm $30 \times 30$ nút tính toán (tương đương 900 tế bào nơron xử lý đồng thời). Bước rời rạc không gian được thiết lập $\Delta x = \Delta y = 0{,}05\text{ m}$ và bước tích phân thời gian $\Delta t = 0{,}00001\text{ s}$ ($10^{-5}\text{ s}$). Phương pháp sai phân hữu hạn xấp xỉ khai triển Taylor cấp hai được lựa chọn nhờ khả năng chuyển đổi trực tiếp các toán tử vi phân không gian thành các trọng số liên kết cục bộ trong ma trận mẫu $A$ và $B$ của mạng CNN. Lý do lựa chọn phương pháp này là vì nó loại bỏ hoàn toàn các khâu biến đổi A/D và D/A trung gian, duy trì tính liên tục của biến trạng thái và tối ưu hóa tính toán song song phần cứng. Nghiên cứu được triển khai theo quy trình kéo dài 12 tháng: từ phân tích toán học, thiết kế mẫu CNN, mô phỏng thuật toán trên MATLAB, đến đặc tả phần cứng bằng ngôn ngữ VHDL trên chip Altera Cyclone II EP2C35.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

Quá trình nghiên cứu và thực nghiệm mô phỏng đã đạt được các phát hiện kỹ thuật nổi bật sau:

Thứ nhất, nghiên cứu đã thiết lập thành công mô hình toán học giải phương trình Navier - Stokes hai chiều bằng mạng CNN 3 lớp liên kết, trong đó mỗi lớp tế bào chịu trách nhiệm tính toán song song một biến trạng thái: độ cao mực nước $h(x,y,t)$, thành phần vận tốc $u(x,y,t)$ theo trục $x$, và thành phần vận tốc $v(x,y,t)$ theo trục $y$. Sai số tương đối giữa nghiệm xấp xỉ của mô hình CNN so với nghiệm giải tích chuẩn được kiểm soát chặt chẽ ở mức dưới $0{,}05%$.

Thứ hai, thời gian đáp ứng tính toán của mỗi bước lặp trên cấu trúc tế bào mạch điện đạt mức thời gian quá độ chỉ $1{,}5\text{ }\mu\text{s}$. Tốc độ này nhanh hơn khoảng 150 đến 200 lần so với việc thực thi các thuật toán lặp số học tuần tự trên bộ vi xử lý máy tính cá nhân cùng thời điểm.

Thứ ba, việc hiện thực hóa kiến trúc mạng CNN trên chip FPGA Altera Cyclone II (EP2C35) chứng minh tính khả thi cao về phần cứng. Thiết kế tận dụng tối ưu tài nguyên khối logic (Logic Elements) thông qua cơ chế đường ống (pipelines) và thanh ghi dịch, duy trì độ ổn định Lyapunov $100%$ trong suốt 1000 chu kỳ tính toán liên tục mà không xuất hiện hiện tượng rẽ nhánh hay phân kỳ nghiệm.

Thứ tư, mô phỏng trên môi trường MATLAB xác nhận rằng biên độ dao động sóng nước và trường vector vận tốc dòng chảy lan truyền hoàn toàn trơn tru, phản ánh chính xác các quy luật bảo toàn động lượng và bảo toàn khối lượng trong thủy lực học mặt biển.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự vượt trội về tốc độ xử lý của mạng CNN nằm ở cơ chế xử lý tương tự song song vật lý (asynchronous parallel analog computing). Khác với các máy tính số truyền thống phải nạp và lưu dữ liệu tuần tự qua bộ nhớ RAM, mỗi nút lưới trên mặt biển được gán cố định cho một tế bào nơron phần cứng. Các tế bào này tương tác tức thời với 8 tế bào láng giềng thông qua các dòng điện hồi tiếp được định hình bởi ma trận mẫu kích thước $3 \times 3$.

Khi so sánh với phương pháp phần tử hữu hạn (FEM) hoặc phương pháp sai phân truyền thống chạy trên CPU đa lõi, kiến trúc CNN loại bỏ được hơn $90%$ lưu lượng truyền thông toàn cục, giải quyết triệt để bài toán nghẽn băng thông dữ liệu.

Dữ liệu mô phỏng được trực quan hóa chi tiết thông qua các biểu đồ không gian 3 chiều (các hàm meshsurf trên MATLAB). Biểu đồ bề mặt biểu diễn phân bố độ cao mực nước $h(x,y)$ tại các thời điểm $t = 0{,}01\text{ s}$ và $t = 0{,}05\text{ s}$ cho thấy sóng dao động lan tỏa đồng đều từ tâm ra các vùng biên. Bảng dữ liệu ma trận so sánh các giá trị vận tốc $u_{i,j}$ và $v_{i,j}$ tại thời điểm ban đầu so với thời điểm $t+\Delta t$ thể hiện rõ sự suy giảm độ xoáy cục bộ do ma sát đáy biển, phù hợp với các kết quả thực nghiệm hải văn học được công bố trong các báo cáo chuyên ngành.

Đề xuất và khuyến nghị

Dựa trên kết quả nghiên cứu, luận văn đề xuất 4 nhóm giải pháp kỹ thuật và khuyến nghị ứng dụng thực tiễn mang tính khả thi cao:

  1. Phát triển vi mạch chuyên dụng ASIC/FPGA mật độ cao: Các đơn vị nghiên cứu phần cứng cần chuyển đổi thiết kế mạng CNN từ dòng chip Cyclone II sang các dòng vi mạch tiên tiến như Xilinx Virtex-6 hoặc Altera Stratix IV tiến trình 40nm. Mục tiêu cụ thể là tăng quy mô mảng tính toán lên $1024 \times 1024$ tế bào, nâng hiệu năng tính toán đạt mức trên $10^{13}$ phép tính/giây trong lộ trình 18 tháng do các kỹ sư vi mạch phụ trách.
  2. Mở rộng mô hình giải thuật sang bài toán dòng chảy 3 chiều: Nhóm nghiên cứu toán tin ứng dụng cần tiếp tục hoàn thiện thuật toán sai phân vi phân tế bào (CPDDE) cho hệ phương trình Navier - Stokes 3D đầy đủ bao gồm cả thành phần vận tốc thẳng đứng $w(x,y,z,t)$, nhằm tăng độ chính xác mô phỏng dòng đối lưu tầng sâu thêm $25%$ trong vòng 24 tháng.
  3. Tích hợp cảm biến môi trường thông minh thời gian thực: Các cơ quan quản lý khí tượng thủy văn cần tích hợp trực tiếp mảng cảm biến quang học (OPT) và cảm biến áp suất đáy biển vào các khối nhớ tương tự cục bộ (LAM) của chip CNN-UM, giúp giảm độ trễ thu thập và xử lý số liệu quan trắc xuống dưới $0{,}1\text{ ms}$ tại các trạm giám sát ven biển trước năm 2028.
  4. Chuẩn hóa bộ công cụ phần mềm tự động hóa thiết kế (EDA Tool): Các trường đại học kỹ thuật nên phối hợp xây dựng thư viện Toolbox chuyên dụng trên MATLAB/Simulink cho phép tự động chuyển đổi phương trình đạo hàm riêng bất kỳ thành cấu hình phần cứng VHDL cho mạng CNN trong thời gian 12 tháng, phục vụ công tác giảng dạy và nghiên cứu chuyên sâu.

Đối tượng nên tham khảo luận văn

Công trình nghiên cứu mang giá trị học thuật và ứng dụng thực tiễn cao, đặc biệt hữu ích cho 4 nhóm đối tượng trọng tâm sau:

  • Học viên cao học và nghiên cứu sinh ngành Khoa học máy tính, Toán tin: Tài liệu cung cấp phương pháp luận chặt chẽ về cách ánh xạ các phương trình đạo hàm riêng phi tuyến sang không gian trạng thái mạng nơron tế bào, làm tài liệu tham khảo cho các đề tài phát triển giải thuật tính toán hiệu năng cao.
  • Kỹ sư thiết kế vi mạch và kỹ sư lập trình phần cứng FPGA/ASIC: Cung cấp sơ đồ khối chi tiết, thiết kế phần cứng theo cơ chế pipeline và giải pháp tối ưu hóa thanh ghi dịch cho các mạng xử lý tín hiệu tương tự-số (Analogic).
  • Chuyên gia mô phỏng khí tượng thủy văn và hải dương học: Nắm bắt giải pháp công nghệ mới giúp giải quyết bài toán dự báo triều cường, tính toán dòng chảy cửa sông, xói lở bờ biển và cảnh báo sóng thần trong thời gian thực với độ chính xác cao.
  • Kỹ sư phát triển hệ thống thị giác máy tính và xử lý ảnh tốc độ cao: Tham khảo kiến trúc máy tính vạn năng CNN-UM để phát triển các dòng camera thông minh có khả năng xử lý và lọc nhiễu ảnh tốc độ trên 50.000 khung hình/giây trong công nghiệp tự động hóa.

Câu hỏi thường gặp

Mạng nơron tế bào (CNN) khác biệt như thế nào so với mạng nơron nhân tạo truyền thống (ANN)?

Khác với ANN sử dụng các liên kết toàn cục phức tạp giữa các lớp, CNN là mảng các tế bào xử lý tương tự chỉ kết nối cục bộ với các láng giềng liền kề trong bán kính $r=1$. Cấu trúc này giúp CNN đạt tốc độ tính toán song song cực cao lên tới $10^{12}$ phép tính/giây, hoạt động liên tục trong thời gian thực và cực kỳ tối ưu khi hiện thực hóa trên chip phần cứng VLSI.

Tại sao phương pháp sai phân lại thích hợp nhất để đưa phương trình Navier - Stokes lên mạng CNN?

Phương pháp sai phân chia miền tính toán thành lưới $30 \times 30$ nút rời rạc, trong đó mỗi điểm nút tương ứng trực tiếp với một tế bào nơron. Các đạo hàm riêng theo không gian được xấp xỉ thành các trọng số liên kết cục bộ kích thước $3 \times 3$, cho phép chuyển đổi hệ phương trình đạo hàm riêng phức tạp thành hệ phương trình trạng thái vi phân thường mà mạch điện CNN có thể giải đồng thời.

Ngôn ngữ MATLAB đóng vai trò gì trong công trình nghiên cứu này?

MATLAB đóng vai trò là môi trường mô phỏng thuật toán số và thẩm định độ chính xác của các mẫu liên kết $A, B$ trước khi chuyển giao sang phần cứng. Với các công cụ tính toán ma trận mạnh mẽ và chức năng vẽ đồ thị 2D/3D chuyên sâu, MATLAB đã mô phỏng thành công sự biến thiên động học của 900 điểm lưới sai phân với sai số dưới $0{,}05%$.

Nền tảng phần cứng FPGA Altera Cyclone II đáp ứng yêu cầu tính toán của mạng CNN ra sao?

Chip FPGA Altera Cyclone II EP2C35 sở hữu hơn 33.000 phần tử logic, cung cấp tài nguyên phần cứng dồi dào để cấu hình các khối tính toán tế bào song song. Kiến trúc này cho phép nạp lại cấu hình linh hoạt bằng ngôn ngữ VHDL và có thể dễ dàng mở rộng sang các dòng chip hiện đại như Stratix IV với 820.000 phần tử logic để xử lý lưới không gian lớn hơn gấp 10 lần.

Ứng dụng thực tiễn quan trọng nhất của luận văn là gì?

Đóng góp thực tiễn lớn nhất của đề tài là xây dựng nền tảng tính toán thời gian thực cho các hệ thống giám sát và cảnh báo thiên tai hải văn. Với thời gian đáp ứng ở mức $1{,}5\text{ }\mu\text{s}$, hệ thống cho phép dự báo biến động sóng biển và kiểm soát lưu lượng dòng chảy nhanh hơn từ 150 đến 200 lần so với các máy tính trạm xử lý số thông thường.

Kết luận

  • Hệ thống hóa toàn diện cơ sở lý thuyết về mạng nơron tế bào, kiến trúc máy tính vạn năng CNN-UM và mối quan hệ tương đương động học với phương trình vi phân đạo hàm riêng.
  • Xây dựng thành công thuật toán sai phân hữu hạn và xác định chính xác bộ mẫu liên kết $3 \times 3$ cho hệ phương trình Navier - Stokes hai chiều mô tả dòng chảy mặt biển.
  • Mô phỏng và kiểm chứng thuật toán trên lưới sai phân $30 \times 30$ điểm bằng MATLAB, đạt độ tin cậy và chính xác toán học trên $99{,}9%$.
  • Thiết kế hoàn chỉnh kiến trúc mạch tính toán song song trên phần cứng FPGA Cyclone II EP2C35, đạt thời gian trễ xử lý cực thấp chỉ $1{,}5\text{ }\mu\text{s}$.
  • Mở ra hướng tiếp cận đột phá trong việc ứng dụng công nghệ mạch tính toán tương tự phục vụ các bài toán cơ học thủy khí trong thời gian thực tại Việt Nam.

Luận văn thạc sĩ của tác giả Nguyễn Thị An Hải là công trình khoa học tiên phong, đặt nền móng vững chắc cho sự kết hợp giữa trí tuệ nhân tạo phần cứng và mô phỏng số hải dương học. Trong giai đoạn tiếp theo từ năm 2026 đến 2028, hướng nghiên cứu sẽ mở rộng sang mô hình 3D đa pha tích hợp chip bán dẫn chuyên dụng. Quý độc giả, các nhà nghiên cứu và kỹ sư vi mạch quan tâm hãy tham khảo toàn văn luận văn để tiếp cận trọn vẹn mô hình toán học và sơ đồ thiết kế chi tiết!