Tổng quan nghiên cứu

Trong kỷ nguyên phát triển mạnh mẽ của khoa học máy tính và trí tuệ nhân tạo, việc giải các phương trình đạo hàm riêng (Partial Differential Equations - PDE) giữ vai trò then chốt trong mô phỏng các quy luật vật lý, cơ học chất lưu và các bài toán kỹ thuật phức tạp. Hầu hết các phương trình đạo hàm riêng trong thực tế không tồn tại nghiệm giải tích chính xác, buộc các nhà khoa học phải dựa vào các phương pháp giải số truyền thống. Tuy nhiên, các kỹ thuật số cổ điển thường đối mặt với rào cản chi phí tính toán khổng lồ khi số chiều không gian tăng cao hoặc khi phương trình chứa các thành phần phi tuyến mạnh. Nhằm khắc phục triệt để thách thức này, công trình nghiên cứu thạc sĩ chuyên ngành Toán Tin thực hiện vào tháng 06 năm 2021 tại Viện Toán ứng dụng và Tin học thuộc Trường Đại học Bách Khoa Hà Nội đã phát triển giải pháp giải xấp xỉ phương trình đạo hàm riêng bằng công nghệ Deep Learning.

Mục tiêu cụ thể của luận văn là xây dựng một khung thuật toán hoàn chỉnh, chuyển đổi bài toán giải hệ vi phân thành bài toán tìm giá trị nhỏ nhất của hàm mất mát không âm trên tập hữu hạn các điểm dữ liệu. Nghiên cứu thực hiện đánh giá toàn diện trên cả không gian 2 chiều và miền không gian - thời gian 3 chiều, bao quát từ phương trình Poisson tuyến tính đến hệ phương trình Navier-Stokes phi tuyến phức tạp. Về mặt ý nghĩa khoa học và thực tiễn, phương pháp này loại bỏ sự phụ thuộc vào việc phân chia lưới hình học phức tạp, tối ưu hóa sai số chuẩn L2 xuống ngưỡng xấp xỉ 0.0001 và xử lý mượt mà 4 bài toán kiểm định mẫu. Đây là bước tiến quan trọng đóng góp vào xu hướng phát triển của lĩnh vực Học máy Khoa học (Scientific Machine Learning) tại Việt Nam.

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

Khung lý thuyết áp dụng

Khung lý thuyết của đề tài được xây dựng trên 2 nền tảng toán học kinh điển về năng lực biểu diễn của mạng nơ-ron nhân tạo. Thứ nhất là Định lý xấp xỉ phổ quát do Kurt Hornik chứng minh năm 1991, khẳng định mạng nơ-ron một tầng ẩn có thể xấp xỉ mọi hàm số liên tục trên tập compact. Thứ hai là các định lý mở rộng của Zhou Lu và Boris Hanin năm 2017 chứng minh cho các mạng nơ-ron sâu (Deep Neural Networks) với số chiều rộng giới hạn sử dụng các hàm kích hoạt phi tuyến như Tanh, Sigmoid và Leaky ReLU. Về mặt mô hình hóa, nghiên cứu tiếp cận phương trình vi phân thông qua 3 khái niệm cốt lõi: quá trình lan truyền xuôi (Forward Propagation) để ước lượng hàm số, phương pháp đạo hàm tự động (Automatic Differentiation theo Evtushenko năm 1991) dựa trên quy tắc đạo hàm hàm hợp để tính toán đạo hàm chính xác mà không cần sai phân rời rạc, và quá trình lan truyền ngược (Backpropagation) để cập nhật tham số trọng số.

Hàm mất mát toàn cục được thiết lập dưới dạng tổng bình phương phần dư, bao gồm sai số phương trình bên trong miền tính toán, sai số điều kiện biên và sai số điều kiện ban đầu. Để giải quyết bài toán tối ưu hóa phi lồi này, luận văn đã khảo sát và phân tích toàn diện 8 thuật toán tối ưu hóa phổ biến nhất hiện nay: Gradient Descent cơ bản, Stochastic Gradient Descent, Momentum, Nesterov Accelerated Gradient (NAG), Adagrad, Adadelta, RMSprop, Adam và Nadam, làm cơ sở lựa chọn phương pháp huấn luyện thích hợp cho từng dạng bài toán vi phân.

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

Nghiên cứu sử dụng 2 kỹ thuật tạo dữ liệu chính: lấy mẫu ngẫu nhiên đồng đều Monte Carlo và lấy mẫu theo lưới vuông tọa độ Cartesian. Cỡ mẫu huấn luyện được phân bổ khoa học theo quy mô từng bài toán: đối với bài toán 2 chiều không phụ thuộc thời gian, tập dữ liệu gồm 400 điểm (lưới 20x20) và 2.500 điểm (lưới 50x50); đối với bài toán 3 chiều không gian - thời gian, quy mô mẫu mở rộng từ 40.000 điểm (lưới 20x20x100) đến 90.000 điểm (lưới 30x30x100). Tập dữ liệu kiểm định độc lập được trích xuất ngẫu nhiên từ 2.000 đến 10.000 điểm để tính toán sai số chuẩn L2 so với nghiệm đúng.

Lý do lựa chọn phương pháp phân tích thực nghiệm dựa trên ngôn ngữ Python kết hợp các framework TensorFlow và Keras là khả năng hỗ trợ tính toán song song trên GPU và tích hợp sẵn đạo hàm tự động mạnh mẽ. Thuật toán tối ưu Adam được lựa chọn làm công cụ điều chỉnh trọng số chủ đạo theo phương thức Mini-batch Gradient Descent với kích thước batch linh hoạt từ 16, 24 đến 48 điểm. Toàn bộ quá trình thử nghiệm mô hình, tinh chỉnh tham số và thu thập kết quả được thực hiện hoàn tất trong lộ trình nghiên cứu 6 tháng của năm 2021.

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

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

Kết quả thực nghiệm trên 4 bài toán phương trình đạo hàm riêng từ đơn giản đến phức tạp đã chứng minh tính hiệu quả vượt bậc của mô hình Deep Learning:

Thứ nhất, đối với phương trình Poisson 2 chiều không phụ thuộc thời gian, mạng nơ-ron gồm 3 lớp ẩn với 32 nơ-ron mỗi lớp đạt giá trị hàm mất mát tối ưu 0.0002041 chỉ sau 295 lần lặp trên lưới 2.500 điểm, rút ngắn số bước lặp khoảng 93% so với mạng 1 lớp ẩn 8 nơ-ron vốn cần tới 4.301 lần lặp để đạt mức mất mát 0.0002050.

Thứ hai, đối với phương trình nhiệt 3 chiều phụ thuộc thời gian, trên tập dữ liệu 90.000 điểm, mô hình 3 lớp ẩn với 32 nơ-ron đạt hàm mất mát 0.0002046 sau 194 lần lặp với thời gian huấn luyện 6.51 giây mỗi bước lặp, vượt trội hoàn toàn so với mô hình 2 lớp ẩn cần tới 3.141 lần lặp.

Thứ ba, ở bài toán Steady Navier-Stokes phi tuyến (bài toán Kovasznay với hệ số độ nhớt 0.025), mô hình đạt hàm mất mát 0.00010058 sau 328 bước lặp trên lưới 2.500 điểm mà không gặp bất kỳ trở ngại nào từ các số hạng phi tuyến.

Thứ tư, trong bài toán hốc chất lưu Navier-Stokes không dừng (Lid-driven Cavity với hệ số độ nhớt 0.01), mạng nơ-ron 4 lớp ẩn với 32 nơ-ron qua 648 bước lặp (khoảng 26 giây mỗi bước) đã xác định tọa độ tâm xoáy ổn định tại vị trí (0.5586, 0.7360), tương đồng ấn tượng với kết quả chuẩn của Ghia et al. là (0.5625, 0.7656), với độ sai lệch tọa độ trục hoành chỉ 0.69% và trục tung khoảng 3.86%.

Thảo luận kết quả

Nguyên nhân chính dẫn đến sự khác biệt về hiệu năng giữa các cấu trúc thử nghiệm là sự cân bằng giữa dung lượng mạng và khả năng tối ưu hóa. Mạng quá nông (1 lớp ẩn) không đủ khả năng nắm bắt độ cong phi tuyến của nghiệm, trong khi mạng quá sâu (từ 4 lớp ẩn trở lên ở bài toán nhiệt) khiến hàm mục tiêu trở nên gồ ghề, dễ làm thuật toán mắc kẹt tại các cực tiểu địa phương, dẫn đến sai số L2 tăng vượt ngưỡng 0.1. Dữ liệu thực nghiệm được trực quan hóa rõ nét qua 4 bảng tổng hợp thông số và các đồ thị đường dòng (streamline), cùng chuỗi hình ảnh mô phỏng 8 giai đoạn tiến triển của dòng xoáy cavity từ bước lặp 1 đến bước lặp 648.

So sánh với các nghiên cứu quốc tế của Raissi et al. và các kết quả thực nghiệm cơ học chất lưu tiêu chuẩn, giải pháp Deep Learning thể hiện ưu thế vượt trội ở khả năng giải quyết các phương trình phi tuyến nhiều chiều mà không cần thiết lập lưới phần tử phức tạp. Tuy nhiên, đối với các bài toán chưa biết trước trạng thái vật lý ổn định, việc dựa hoàn toàn vào việc giảm hàm mất mát đòi hỏi người nghiên cứu phải thận trọng kiểm chứng để tránh hội tụ vào các nghiệm giả không mang ý nghĩa vật lý.

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

Thứ nhất, tự động hóa tối ưu hóa siêu tham số (AutoML/NAS): Đề xuất các nhóm nghiên cứu Toán Tin và Khoa học Dữ liệu ứng dụng kỹ thuật tìm kiếm kiến trúc mạng nơ-ron tự động để xác định chính xác số lớp ẩn tối ưu (từ 2 đến 3 lớp) và số nơ-ron (16 đến 32 đơn vị) cho từng dạng phương trình, giúp cắt giảm khoảng 50% thời gian thử nghiệm thủ công trong giai đoạn 2026-2027.

Thứ hai, triển khai chiến lược tối ưu hóa lai hai giai đoạn (Hybrid Optimization): Khuyến nghị các kỹ sư mô phỏng kết hợp thuật toán Adam trong 500 bước lặp đầu để thoát nhanh khỏi các cực tiểu địa phương, sau đó chuyển sang thuật toán L-BFGS nhằm đẩy nhanh tốc độ hội tụ bậc hai và đưa hàm mất mát xuống dưới ngưỡng 0.000001 trong vòng 6 tháng tới.

Thứ ba, xây dựng cơ chế lấy mẫu thích nghi cục bộ (Residual-based Adaptive Refinement): Thay vì sử dụng lưới tĩnh phân bổ đều từ 400 đến 90.000 điểm, cần phát triển thuật toán tự động bổ sung điểm dữ liệu tại các khu vực có gradient biến thiên mạnh như lớp biên chất lưu và tâm xoáy xoay chiều, giúp tiết kiệm ít nhất 40% tài nguyên tính toán GPU.

Thứ tư, mở rộng ứng dụng vào các bài toán kỹ thuật thực tế quy mô lớn: Đề xuất các viện nghiên cứu chuyên ngành hàng không và khí tượng thủy văn ứng dụng khung mô hình này để giải các bài toán động lực học chất lưu Navier-Stokes 3D với số Reynolds lớn hơn 10.000 trong giai đoạn 2026-2030.

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

Nhóm 1: Học viên cao học và nghiên cứu sinh chuyên ngành Toán Tin, Toán ứng dụng: Nắm vững cơ sở toán học về xấp xỉ hàm liên tục, cơ chế hội tụ của 8 thuật toán tối ưu hóa Gradient Descent và phương pháp biến đổi phương trình vi phân thành bài toán tối ưu hàm mất mát.

Nhóm 2: Kỹ sư mô phỏng cơ học tính toán (CFD) và kỹ thuật thủy khí: Khai thác giải pháp mạng nơ-ron để mô phỏng dòng chảy Navier-Stokes và bài toán Kovasznay mà không cần phụ thuộc vào phần mềm thương mại chia lưới đắt đỏ, nâng cao độ chính xác dự báo động học chất lưu.

Nhóm 3: Chuyên gia phát triển phần mềm Trí tuệ nhân tạo và Machine Learning: Tham khảo kiến trúc mã nguồn trên TensorFlow và Keras, tận dụng cơ chế đạo hàm tự động để phát triển các ứng dụng Học máy Vật lý (Physics-Informed Neural Networks - PINNs) với độ chính xác đạt trên 99%.

Nhóm 4: Giảng viên và nhà nghiên cứu tại các trường đại học kỹ thuật: Sử dụng luận văn làm tài liệu giảng dạy chuyên đề sau đại học về Deep Learning ứng dụng trong khoa học tính toán, tạo tiền đề mở rộng các hướng nghiên cứu liên ngành Toán học và Trí tuệ nhân tạo.

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

Phương pháp giải phương trình đạo hàm riêng bằng Deep Learning khác gì so với phương pháp số truyền thống? Thay vì phải chia lưới không gian rời rạc và giải hệ đại số tuyến tính khổng lồ như phương pháp phần tử hữu hạn (FEM) hay sai phân hữu hạn (FDM), Deep Learning biểu diễn nghiệm dưới dạng hàm liên tục thông qua mạng nơ-ron và đạo hàm tự động, giúp giải quyết hiệu quả không gian 3 chiều mà không chịu ảnh hưởng của lời nguyền số chiều.

Tại sao thuật toán Adam lại được ưu tiên lựa chọn để huấn luyện mô hình trong luận văn? Thuật toán Adam là sự kết hợp tối ưu giữa cơ chế duy trì quán tính của Momentum và khả năng tự điều chỉnh tốc độ học của RMSprop. Sự phối hợp này giúp mô hình dễ dàng vượt qua các điểm cực tiểu địa phương và nhanh chóng hội tụ hàm mất mát về ngưỡng 0.00010058 sau khoảng 300 đến 600 bước lặp.

Hàm mất mát trong việc giải phương trình đạo hàm riêng được thiết lập như thế nào? Hàm mất mát toàn cục được xây dựng từ 3 thành phần chính: tổng bình phương phần dư của toán tử vi phân bên trong miền tính toán, sai số tại các điểm biên Dirichlet và sai số tại điều kiện ban đầu. Khi giá trị hàm mất mát tiệm cận về 0, đầu ra của mạng nơ-ron sẽ tiến sát về nghiệm chính xác của bài toán.

Cấu trúc số lớp ẩn và số nơ-ron ảnh hưởng như thế nào đến độ chính xác của nghiệm? Thực nghiệm chứng minh cấu trúc mạng có độ phức tạp vừa phải gồm 2 đến 3 lớp ẩn với 16 đến 32 nơ-ron đem lại kết quả tối ưu nhất. Mạng quá đơn giản (1 lớp ẩn) không đủ khả năng xấp xỉ, trong khi mạng quá phức tạp (trên 4 lớp ẩn) khiến sai số L2 tăng vượt mức 0.1 do quá trình tối ưu hóa tham số gặp khó khăn.

Nhược điểm lớn nhất của phương pháp giải phương trình đạo hàm riêng bằng mạng nơ-ron là gì? Hạn chế lớn nhất là khó kiểm soát sai số tuyệt đối khi giải các bài toán không có nghiệm giải tích để đối chiếu. Ngoài ra, việc lựa chọn kiến trúc mạng vẫn phụ thuộc nhiều vào thử nghiệm thực nghiệm, với thời gian huấn luyện có thể kéo dài tới 26 giây cho mỗi bước lặp ở các hệ phương trình phi tuyến phức tạp.

Kết luận

  • Đã thiết lập thành công quy trình 4 bước chuẩn mực để giải xấp xỉ các phương trình đạo hàm riêng bằng mạng nơ-ron sâu trên nền tảng Python, TensorFlow và Keras.
  • Hệ thống hóa và phân tích chuyên sâu lý thuyết cùng ưu nhược điểm của 8 thuật toán tối ưu hóa Gradient Descent, làm rõ cơ chế thích nghi tốc độ học trong bài toán vi phân.
  • Kiểm chứng thực nghiệm thành công trên 4 bài toán từ tuyến tính đến phi tuyến Navier-Stokes, đạt giá trị hàm mất mát tối ưu trong khoảng từ 0.0001 đến 0.0002.
  • Xác định cấu trúc mạng nơ-ron tối ưu gồm 2 đến 3 lớp ẩn với 16 đến 32 nơ-ron, giúp cân bằng hoàn hảo giữa năng lực biểu diễn và thời gian tính toán.
  • Mô phỏng chuẩn xác tọa độ tâm xoáy hiện tượng hốc chất lưu tại vị trí (0.5586, 0.7360) với độ sai lệch tọa độ dưới 4% so với các nghiên cứu thực nghiệm quốc tế.

Đóng góp then chốt của luận văn là đã đặt nền móng thực nghiệm vững chắc cho phương pháp giải phương trình vi phân bằng Deep Learning tại Đại học Bách Khoa Hà Nội, mở ra lộ trình phát triển các thuật toán tối ưu lai và lấy mẫu thích nghi trong giai đoạn 2026-2027. Bạn đọc hãy tham khảo toàn văn luận văn thạc sĩ này để tiếp cận trọn vẹn khung phương pháp luận hiện đại và ứng dụng hiệu quả Deep Learning vào các bài toán mô phỏng khoa học kỹ thuật chuyên sâu.