Tổng quan về luận án

Bài toán tương tác giữa dòng lưu chất và kết cấu rắn (Fluid-Structure Interaction - FSI) là một trong những thách thức trung tâm của cơ học tính toán đương đại, hiện diện phức tạp trong khí động lực học công trình, dao động cánh turbine gió, phản ứng khí phiếm của máy bay và động lực học dòng chảy sinh học như van tim nhân tạo hay sinh thể bơi lội. Nghiên cứu sinh Lê Quốc Cường (chuyên ngành Cơ kỹ thuật, mã số 9520101, Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh) đã thực hiện luận án tiến sĩ với đề tài: "Phát triển phương pháp biên nhúng kết hợp phương pháp Proper Generalized Decomposition (PGD) cho bài toán dòng chảy nhớt không nén được qua vật thể biên cứng và biên đàn hồi" dưới sự hướng dẫn khoa học của TS. Phan Đức Huynh. Luận án giải quyết một điểm nghẽn học thuật cốt lõi: sự bùng nổ chi phí tính toán (Curse of Dimensionality) và sự phụ thuộc nặng nề vào quá trình tái tạo lưới (remeshing) khi giải hệ phương trình Navier-Stokes liên kết đa trường trên các miền hình học phức tạp có biên biến dạng hoặc dao động phi tuyến.

Khoảng trống nghiên cứu then chốt được xác định rõ ràng: các phương pháp biên nhúng truyền thống (Immersed Boundary Method - IBM) dù loại bỏ được nhu cầu tạo lưới động tương thích theo mặt phân cách, nhưng khi tích hợp với các kỹ thuật rời rạc hóa cổ điển (sai phân hữu hạn FDM, thể tích hữu hạn FVM, phần tử hữu hạn FEM) thì số bậc tự do vẫn tăng theo hàm số mũ $\mathcal{O}(M^D)$ trong không gian $D$ chiều (với $M$ là số nút trên một trục tọa độ). Trong khi đó, các mô hình giảm bậc (Reduced-Order Models - ROM) như Proper Orthogonal Decomposition (POD) lại phụ thuộc hoàn toàn vào dữ liệu nghiệm quá khứ (snapshots) thu được từ các mô phỏng tốn kém trước đó. Luận án đặt ra hệ thống câu hỏi nghiên cứu và giả thuyết khoa học:

  • RQ1: Làm thế nào để xây dựng một thuật toán xấp xỉ tiên nghiệm (a priori) tách biến không gian cho các phương trình vi phân đạo hàm riêng bậc cao (Poisson, Biharmonic) mà không cần dữ liệu snapshot ban đầu?
  • RQ2: Cơ chế ghép nối giữa trường lực tương tác Lagrange-Euler trong IBM và sơ đồ giải PGD cho phương trình Poisson áp suất trong hệ Navier-Stokes hoạt động như thế nào để đảm bảo tính ổn định và điều kiện không trượt tại biên vật thể?
  • RQ3: Mức độ tối ưu hóa về thời gian xử lý CPU và dung lượng bộ nhớ RAM của thuật toán kết hợp IBM-PGD đạt được bao nhiêu phần trăm so với phương pháp sai phân hữu hạn truyền thống trên cùng một độ mịn của lưới?
  • H1: Việc biểu diễn nghiệm dưới dạng chuỗi tách biến $u(x,y,z) \approx \sum_{i=1}^N X_i(x)Y_i(y)Z_i(z)$ thông qua giải thuật lặp luân phiên điểm cố định (Fixed-point alternating direction strategy) sẽ chuyển đổi hệ phương trình đạo hàm riêng 2D/3D thành một chuỗi các phương trình vi phân thường 1D, giúp giảm cấp độ phức tạp tính toán từ $\mathcal{O}(M^3)$ xuống $\mathcal{O}(Q \cdot N \cdot M)$.
  • H2: Sự kết hợp đồng bộ giữa lực cưỡng bức biên nhúng $f$ và giải thuật chiếu PGD cho phép mô phỏng chính xác các hiện tượng thủy động lực học phi tuyến (tách xoáy von Kármán, dao động khóa pha lock-in của vật thể đàn hồi) với thời gian thực thi giảm từ 60% đến 85%.

Phạm vi nghiên cứu bao gồm mô hình hóa dòng chảy tầng của chất lưu nhớt không nén được ở số Reynolds thấp và trung bình ($Re \le 5000$), khảo sát tương tác với vật thể biên cứng cố định, vật thể dao động tịnh tiến và cắt ngang, cùng các kết cấu sợi và màng đàn hồi chịu lực uốn và lực kéo phi tuyến trong không gian hai chiều.

Literature Review và Positioning

Lịch sử phát triển của các phương pháp tính toán FSI chứng kiến sự phân hóa sâu sắc giữa hai trường phái: lưới tương thích biên (Body-conformal grid methods như ALE-FEM) và lưới Descartes phi cấu trúc (Non-conformal Cartesian grid methods). Trường phái IBM được khởi xướng bởi Peskin (1977) để nghiên cứu cơ học huyết động qua van tim sinh học, sử dụng hàm Dirac delta xấp xỉ trơn để truyền lực từ lưới Lagrange của biên đàn hồi sang lưới Euler của lưu chất. Sau đó, Mohd-Yosuf (1997) cùng Fadlun và cộng sự (2000) đã phát triển phương pháp biên nhúng cưỡng bức trực tiếp (Direct Forcing IBM), tái cấu trúc trường vận tốc cục bộ để giải quyết bài toán biên cứng. Uhlmann (2005) thiết lập sơ đồ cưỡng bức trực tiếp hiện trên lưới Euler-Lagrange, trong khi Taira và Colonius (2007) đề xuất phương pháp chiếu IB (Immersed Boundary Projection Method) đưa lực biên thành nhân tử Lagrange giải đồng thời với áp suất qua phương trình Poisson mở rộng.

Tuy nhiên, các tranh luận học thuật lớn tồn tại giữa hai cách tiếp cận:

  • Quan điểm 1 (Mittal & Iaccarino, 2005; Tseng & Ferziger, 2003): Ủng hộ phương pháp ô ảo (Ghost-cell IBM) và cắt ô (Cut-cell IBM) nhằm bảo toàn độ chính xác không gian bậc hai sắc nét tại mặt phân cách, nhưng đối mặt với sự bất ổn định số học do trọng số âm khi điểm biên quá gần nút lưới Euler và bế tắc hình học khi mở rộng sang không gian 3D phức tạp.
  • Quan điểm 2 (Kempe & Fröhlich, 2012; Su et al., 2007): Ủng hộ sơ đồ cưỡng bức trực tiếp ẩn để thỏa mãn chính xác điều kiện không trượt tại biên nhúng, song lại đòi hỏi tài nguyên tính toán khổng lồ để nghịch đảo các ma trận đại số lớn liên kết toàn cục.

Song song đó, nhóm phương pháp giảm bậc mô hình truyền thống như POD (Sirovich, 1987; Holmes et al., 1996) đòi hỏi phải có cơ sở dữ liệu nghiệm phong phú (a posteriori), dẫn đến chi phí tiền xử lý cực lớn. Nhằm vượt qua rào cản này, Ammar, Mokdad, Chinesta và Keunings (2006, 2007) đã đề xuất phương pháp Proper Generalized Decomposition (PGD), xây dựng lời giải tách biến tiên nghiệm không gian-thời gian hoặc đa thông số. Dumon và cộng sự (2011, 2013) đã bước đầu thử nghiệm PGD kết hợp thể tích hữu hạn và phương pháp phổ cho hệ Navier-Stokes nhưng chỉ giới hạn trong miền tính toán hình chữ nhật thuần túy không có vật cản.

Vị trí học thuật của luận án được xác lập tại điểm giao thoa đột phá: lần đầu tiên thiết lập khung tính toán tích hợp IBM-PGD toàn diện. Khác với các nghiên cứu quốc tế của Dumon et al. (2011) vốn bất lực trước các miền hình học phức tạp, và khác với phương pháp IBM truyền thống của Kim & Choi (2000) chịu gánh nặng chi phí lưới mịn, mô hình của luận án cho phép xử lý mọi hình học biên cứng và biên đàn hồi phức tạp trên lưới Descartes cố định, đồng thời PGD đảm nhận vai trò giải phương trình Poisson áp suất với chi phí thuật toán tối thiểu.

Đó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 nền tảng lý thuyết phân rã không gian của Ammar et al. (2006) và lý thuyết lực cưỡng bức biên nhúng của Peskin (2002) sang cơ học chất lưu tính toán đa trường:

  1. Mở rộng lý thuyết PGD cho các toán tử vi phân bậc cao: Thiết lập thành công dạng yếu biến phân và sơ đồ lặp luân phiên tìm nghiệm cho phương trình Poisson $\nabla^2 u = f(x,y,z)$ và phương trình Biharmonic $\nabla^4 u = f(x,y,z)$ trong không gian 2D và 3D. Nghiệm tách biến được chuẩn hóa: $$u(x, y, z) \approx \sum_{i=1}^N X_i(x) Y_i(y) Z_i(z)$$ Chuyển hóa toán tử đạo hàm riêng bậc 4 thành hệ các phương trình vi phân thường bậc 4 một chiều giải tích - số học chuẩn xác.
  2. Khung lý thuyết tương tác màng đàn hồi Lagrange phi tuyến: Xây dựng hệ thức liên kết giữa mật độ lực biến dạng trên biên nhúng thông qua năng lượng biến dạng đàn hồi phi tuyến, bao gồm lực kéo căng $\mathbf{F}_s$ và lực kháng uốn $\mathbf{F}_b$: $$\mathbf{F}(s, t) = \frac{\partial}{\partial s} \left( T(s) \boldsymbol{\tau}(s) \right) - \frac{\partial^2}{\partial s^2} \left( E_b \frac{\partial^2 \mathbf{X}(s, t)}{\partial s^2} \right)$$ trong đó $T(s) = K_s \left( \left| \frac{\partial \mathbf{X}}{\partial s} \right| - 1 \right)$ là lực căng với độ cứng kéo $K_s$, và $E_b$ là độ cứng chống uốn của sợi đàn hồi.

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

Khung phân tích của luận án tích hợp ba lý thuyết nền tảng:

  • Lý thuyết môi trường liên tục Navier-Stokes cho pha lưu chất nhớt không nén.
  • Lý thuyết cơ học kết cấu phi tuyến 1D/2D cho biên đàn hồi theo tọa độ vật chất Lagrange.
  • Không gian hàm Hilbert $L^2(\Omega)$ và giải thuật phân rã ten-sơ PGD.

Quy trình giải tích số học được thiết lập chặt chẽ:

[Bước thời gian t^n]
[1. Nội suy vận tốc Euler sang biên Lagrange qua Dirac Delta: U_b = I(u^n)]
[2. Tính lực phản hồi biên đàn hồi F(s) hoặc lực cưỡng bức biên cứng f(X)]
[3. Phân bố lực lên lưới Euler: f(x) = S(F)]
[4. Tính trường vận tốc trung gian u* qua phương pháp phân bước Chorin]
[5. Áp dụng PGD giải phương trình Poisson áp suất: ∇²p^(n+1) = (ρ/Δt) ∇·u*]
[6. Cập nhật vận tốc solenoidal u^(n+1) và vị trí kết cấu X^(n+1)]

Điều kiện biên hạn định (Boundary conditions) được kiểm soát nghiêm ngặt: điều kiện Dirichlet/Neumann trên biên ngoài của miền lưu chất $\Omega$, và điều kiện không trượt $\mathbf{u}(\mathbf{X}(s,t), t) = \mathbf{U}_b(s,t)$ được thỏa mãn tự nhiên tại mặt phân cách nhúng $\Gamma$.

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 hoàn toàn trên hệ hình nhận thức luận thực chứng (Positivism) kết hợp phương pháp luận mô phỏng số tất định (Deterministic Numerical Simulation). Thiết kế nghiên cứu là sự tích hợp đa tầng giữa giải tích toán học, phương pháp rời rạc sai phân hữu hạn và kỹ thuật lập trình khoa học. Toàn bộ thuật toán được hiện thực hóa độc lập trên nền tảng ngôn ngữ ma trận cấp cao MATLAB, không phụ thuộc vào bất kỳ phần mềm thương mại hộp đen nào (như ANSYS Fluent hay COMSOL).

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu tuân thủ các chuẩn mực nghiêm ngặt của cơ học tính toán:

  1. Rời rạc hóa không gian trên lưới so le (Staggered MAC Grid): Biến áp suất $p$ được định vị tại tâm ô lưới $(i, j)$, trong khi các thành phần vận tốc $u$ và $v$ lần lượt được định vị tại các mặt bên $(i \pm 1/2, j)$ và $(i, j \pm 1/2)$, loại bỏ hoàn toàn hiện tượng phân kỳ áp suất răng cưa (checkerboard pressure oscillations).
  2. Xấp xỉ hàm Dirac Delta: Luận án áp dụng hàm delta rời rạc trơn 4 điểm của Peskin: $$d_h(r) = \begin{cases} \frac{1}{4h} \left( 1 + \cos\left(\frac{\pi r}{2h}\right) \right), & |r| \le 2h \ 0, & |r| > 2h \end{cases}$$ đảm bảo bảo toàn tuyệt đối tổng lực và moment chuyển đổi giữa hai hệ tọa độ Euler và Lagrange.
  3. Tiêu chuẩn hội tụ kép: Luận án áp dụng hai vòng kiểm soát sai số:
    • Vòng lặp PGD cục bộ: Hội tụ khi $|R^{(q)}S^{(q)} - R^{(q-1)}S^{(q-1)}|{L_2} \le \varepsilon{RST} = 10^{-6}$.
    • Vòng lặp PGD toàn cục: Dừng khi sai số thặng dư của chuỗi đạt $|E_N|_{L_2} \le \varepsilon_u = 10^{-6}$.

Data và phân tích

Ma trận kiểm chứng dữ liệu được triển khai chi tiết qua các bài toán benchmark kinh điển:

  • Phương trình vi phân đạo hàm riêng bậc cao: Lưới tính toán biến thiên từ $32 \times 32$, $64 \times 64$ đến $100 \times 100$ điểm trên mỗi chiều. Trong bài toán phương trình Poisson 2D $\nabla^2 u = -2\pi^2 \sin(\pi x)\sin(\pi y)$, lời giải PGD đạt độ chính xác sai số tuyệt đối cực đại $|u_{PGD} - u_{ex}| \le 4.12 \times 10^{-5}$ chỉ sau 1 mode PGD duy nhất ($N=1$).
  • Bài toán hốc vuông Lid-Driven Cavity: Kích thước miền $\Omega = [0, 1] \times [0, 1]$, khảo sát các số Reynolds $Re = 100, 400, 1000, 3200, 5000$. Dữ liệu vận tốc $u(0.5, y)$ và $v(x, 0.5)$ được đối sánh trực tiếp với bộ dữ liệu chuẩn quốc tế của Ghia, Ghia và Shin (1982).
  • Bài toán dòng chảy qua bước thụt (Backward-Facing Step): Khảo sát dòng chảy tầng trong dải $Re = 100 - 800$, kiểm chứng sự phát triển của chiều dài vùng xoáy hoàn lưu hạ lưu $x_1/S$ và $x_2/S$ so với số liệu thực nghiệm của Armaly và cộng sự (1983).

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

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

  1. Sự triệt tiêu hiện tượng thắt cổ chai tính toán trong phương trình Poisson áp suất: Việc áp dụng PGD cho phương trình Poisson áp suất trong dòng chảy nhớt không nén cho phép rút ngắn thời gian giải ma trận đại số từ bậc $\mathcal{O}(M^4)$ (đối với phương pháp trực tiếp) xuống $\mathcal{O}(N \cdot M)$, giúp tổng thời gian tính toán của toàn bộ chu trình Navier-Stokes giảm trung bình từ 65% đến 82% so với sơ đồ sai phân hữu hạn giải lặp truyền thống (Gauss-Seidel hoặc SOR).
  2. Độ chính xác cao trong bài toán dòng chảy qua trụ tròn cố định: Tại $Re = 20$ và $Re = 40$, chiều dài vùng xoáy tuần hoàn sau trụ đạt $L_w/D = 0.92$ (tại $Re=20$) và $L_w/D = 2.24$ (tại $Re=40$), hệ số lực cản đạt $C_d = 2.05$ và $C_d = 1.54$, hoàn toàn trùng khớp với kết quả giải tích và mô phỏng của Tritton (1959) và Kim & Choi (2000). Tại $Re = 100$ và $Re = 200$, hiện tượng tách dòng không ổn định sinh ra đường xoáy von Kármán tuần hoàn với tần số không thứ nguyên hóa qua số Strouhal đạt $St = 0.165$ ($Re=100$) và $St = 0.198$ ($Re=200$), biên độ dao động hệ số nâng $C_l = \pm 0.33$ và $C_l = \pm 0.68$.
  3. Hiện tượng khóa pha (Lock-in) và chuyển đổi trạng thái dòng của hai sợi đàn hồi: Trong khảo sát dòng chảy qua hai sợi đàn hồi đặt song song trong kênh dẫn, luận án phát hiện khi khoảng cách tương đối giữa hai sợi $d/L \le 0.5$, hai sợi dao động cùng pha (in-phase) do sự tương tác thủy động lực học ghép cặp mạnh mẽ. Khi nới rộng khoảng cách $d/L \ge 1.0$, trạng thái chuyển dịch sang dao động ngược pha (out-of-phase) với tần số dao động bám sát tần số tách xoáy tự nhiên.
  4. Bảo toàn diện tích màng đàn hồi khép kín: Khảo sát bài toán màng đàn hồi hình quả bóng biến dạng trong chất lưu tĩnh chứng minh rằng thuật toán IBM-PGD duy trì mức độ mất mát diện tích (area loss) cực thấp, chỉ $0.081%$ tại thời điểm ổn định $t = 0.020\text{ s}$, khắc phục nhược điểm rò rỉ thể tích cố hữu của phương pháp biên nhúng cổ điển Peskin.

Implications đa chiều

  • Về mặt lý thuyết: Chứng minh tính khả thi của việc áp dụng đại số tuyến tính giảm bậc PGD vào các hệ phương trình đạo hàm riêng phi tuyến ghép nối biên động, mở đường cho việc xây dựng các công thức biến phân tổng quát cho bài toán tương tác đa trường.
  • Về mặt phương pháp luận: Thiết lập một nền tảng thuật toán mã nguồn mở linh hoạt, cho phép tích hợp dễ dàng các mô hình cấu lượng đàn hồi phi tuyến tính cao (Hyperelasticity, Viscoelasticity) mà không cần can thiệp tái tạo cấu trúc lưới không gian Euler.
  • Về mặt ứng dụng thực tiễn: Tạo tiền đề cho việc xây dựng các công cụ tính toán thủy động lực học thời gian thực (Real-time CFD) phục vụ thiết kế y sinh (mô phỏng huyết học trong vi mạch, van tim cơ học) và tối ưu hóa kết cấu cánh turbine trong kỹ thuật năng lượng tái tạo.

Limitations và Future Research

Luận án thẳng thắn thừa nhận các giới hạn kỹ thuật nội tại:

  1. Giới hạn số Reynolds: Mô hình hiện tại mới chỉ giải quyết dòng chảy tầng đẳng nhiệt ở số Reynolds thấp và trung bình ($Re \le 5000$). Khi chuyển tiếp sang dòng chảy rối ($Re > 10^5$), các cấu trúc xoáy vi mô đòi hỏi các mô hình rối (LES hoặc RANS) kết hợp hàm vách (wall functions) chưa được tích hợp vào sơ đồ PGD hiện hành.
  2. Không gian hình học 2D: Mặc dù cơ sở lý thuyết PGD cho phương trình Poisson và Biharmonic đã được phát triển cho không gian 3D, phần ghép nối tương tác IBM-PGD cho bài toán FSI trong luận án chủ yếu mới được thử nghiệm toàn diện trên các mô hình hai chiều.
  3. Xử lý hội tụ PGD với tỷ số co giãn hình học lớn: Khi miền tính toán có tỷ số co giãn (aspect ratio) cực lớn hoặc dị hướng mạnh, số lượng mode PGD ($N$) yêu cầu để đạt hội tụ tăng lên, làm giảm phần nào ưu thế tốc độ tính toán của phương pháp.

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

  • Phát triển mô hình IBM-PGD 3D song song hóa trên nền tảng tính toán hiệu năng cao GPU/CUDA.
  • Mở rộng thuật toán cho dòng chảy rối đa pha có tính nén nhẹ kết hợp truyền nhiệt đối lưu.
  • Ứng dụng mô phỏng chi tiết tương tác dòng máu và van tim nhân tạo 3D với cơ tính dị hướng phi tuyến.

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

Công trình của nghiên cứu sinh Lê Quốc Cường tạo ra tác động học thuật và kỹ thuật đáng kể:

  • Tác động học thuật: Cung cấp phương pháp luận mới kết hợp ROM-PGD và IBM, mở ra hướng trích dẫn tiềm năng lớn trong các tạp chí quốc tế chuyên ngành cơ học tính toán (như Journal of Computational Physics, Computer Methods in Applied Mechanics and Engineering).
  • Chuyển giao công nghệ và R&D công nghiệp: Cung cấp thuật toán cốt lõi cho các trung tâm nghiên cứu thiết kế thiết bị trao đổi nhiệt, turbine gió trục đứng và kỹ thuật hàng hải mô phỏng động lực học chân vịt tàu thủy.
  • Giáo dục đào tạo sau đại học: Đóng góp bộ mã nguồn mẫu MATLAB chuẩn mực cho giảng dạy giải tích số bậc cao và cơ học chất lưu tính toán tại các trường đại học kỹ thuật.

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

  • Nghiên cứu sinh và học viên cao học ngành Cơ kỹ thuật/Toán ứng dụng: Tiếp cận một cấu trúc thuật toán PGD minh bạch, có thể tái lập (reproducible) để phát triển các đề tài liên quan đến phương trình vi phân phi tuyến.
  • Kỹ sư R&D Khí - Thủy động lực học: Sở hữu công cụ mô phỏng sơ bộ nhanh chóng, cho phép khảo sát hàng loạt thông số thiết kế (parametric study) trong thời gian ngắn mà không tốn chi phí bản quyền phần mềm đắt đỏ.
  • Nhà hoạch định và quản lý dự án y sinh học tính toán: Sử dụng dữ liệu mô phỏng tương tác sợi và màng đàn hồi để hỗ trợ thẩm định các thiết bị y tế cấy ghép tim mạch.

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 thành công lý thuyết xấp xỉ phân rã chính tắc PGD của Ammar et al. (2006) để giải phương trình Biharmonic 3D dạng yếu $\nabla^4 u = f(x,y,z)$ và tích hợp giải thuật này vào bước giải phương trình Poisson áp suất của hệ phương trình Navier-Stokes kết hợp biên nhúng, biến một bài toán tương tác phức tạp đa chiều thành chuỗi các bài toán 1D giải tích-số học hiệu quả cao.

  2. Phương pháp luận của luận án có gì đổi mới khi so sánh với hai nghiên cứu quốc tế điển hình?

    • So với Dumon và cộng sự (2011): Dumon áp dụng PGD cho Navier-Stokes trên miền chữ nhật đơn giản và không thể xử lý vật cản; luận án đã khắc phục triệt để bằng cách tích hợp lực biên nhúng $f$, cho phép giải bài toán qua mọi hình học phức tạp mà vẫn giữ nguyên ưu thế tách biến của PGD.
    • So với Taira & Colonius (2007) (Immersed Boundary Projection Method): Phương pháp của Taira tạo ra ma trận Poisson áp suất hiệu chỉnh khổng lồ và bị suy biến điều kiện số khi biên di chuyển; phương pháp của luận án tách biệt bước tính áp suất thuần túy qua PGD, giữ ma trận hệ số 1D cố định theo thời gian, giúp tiết kiệm bộ nhớ vượt trội.

  3. Phát hiện số học nào gây bất ngờ nhất trong quá trình thực nghiệm mô phỏng? Đó là việc chỉ cần duy nhất $N=1$ đến $N=2$ mode PGD là thuật toán đã tái tạo chính xác trên 99.9% nghiệm của phương trình Poisson áp suất trên lưới so le, đồng thời giúp giảm thời gian chạy chương trình tới 82% trong bài toán dòng chảy qua trụ tròn dao động ngang mà không hề gây mất ổn định số học.

  4. Luận án có cung cấp giao thức tái lập nghiên cứu (Replication Protocol) không? Có. Luận án cung cấp đầy đủ mã nguồn chương trình MATLAB hoàn chỉnh trong phần Phụ lục (Phụ lục 1: Chương trình giải phương trình Poisson 2D bằng PGD; Phụ lục 2: Chương trình mô phỏng dòng chảy qua trụ tròn bằng IBM-PGD), cùng toàn bộ các thông số vật lý, kích thước miền lưới và bước thời gian $\Delta t$.

  5. Chương trình nghiên cứu 10 năm tiếp theo từ nền tảng luận án được định hình ra sao? Lộ trình 10 năm hướng tới: (1) Xây dựng thư viện mã nguồn mở song song IBM-PGD 3D; (2) Tích hợp học máy vật lý (Physics-Informed Neural Networks - PINNs) với PGD để tối ưu hóa việc tìm các hàm cơ sở; (3) Ứng dụng mô phỏng động lực học dòng máu vi tuần hoàn tương tác với tế bào hồng cầu biến dạng đàn hồi trong các vi kênh sinh học.

Kết luận

  1. Luận án đã phát triển thành công phương pháp biên nhúng (IBM) kết hợp phương pháp Proper Generalized Decomposition (PGD) giải quyết trọn vẹn bài toán dòng chảy nhớt không nén qua vật thể biên cứng và biên đàn hồi.
  2. Thuật toán PGD giải phương trình Poisson và Biharmonic 2D/3D chứng minh tính chính xác vượt trội với sai số tuyệt đối nhỏ hơn $10^{-4}$ và giảm bậc phức tạp tính toán từ đa chiều về một chiều.
  3. Việc tích hợp IBM-PGD vào hệ Navier-Stokes giải quyết triệt để vấn đề chia lại lưới động, đồng thời giảm thời gian tính toán từ 65% đến 82% so với phương pháp sai phân hữu hạn thông thường.
  4. Kết quả mô phỏng các bài toán chuẩn (Lid-driven cavity, Backward-facing step, Dòng chảy qua trụ tròn cố định và dao động) đạt sự tương thích hoàn hảo với các dữ liệu thực nghiệm và giải tích quốc tế ($Re \le 5000$).
  5. Khảo sát tương tác với một sợi, hai sợi đàn hồi và màng khép kín làm sáng tỏ các cơ chế thủy động lực học phi tuyến phức tạp (khóa pha, bảo toàn diện tích màng với độ mất mát chỉ 0.081%).
  6. Công trình tạo nền tảng vững chắc cho sự phát triển của các công cụ mô phỏng FSI hiệu năng cao trong kỹ thuật cơ khí, hàng không và y sinh học tính toán hiện đại.