Tổng quan về luận án

Bài toán tương tác rắn – lỏng (Fluid-Structure Interaction – FSI) là một trong những thách thức tính toán phức tạp nhất trong cơ học kỹ thuật hiện đại, chi phối các hiện tượng vật lý đa trường từ khí động lực học công trình cầu đường, tuabin gió, dao động khí đàn hồi của cánh máy bay đến huyết động học tim mạch. Luận án tiến sĩ kỹ thuật của tác giả Lê Quốc Cường (2019) tại Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh mang tính tiên phong khi giải quyết triệt để rào cản tính toán trong mô phỏng FSI: hiện tượng bùng nổ bậc tự do (curse of dimensionality) khi chia lưới mịn và chi phí tái tạo lưới (remeshing) cực lớn khi vật thể chuyển động hoặc biến dạng phi tuyến.

Khoảng trống nghiên cứu (research gap) trọng tâm xuất phát từ việc các phương pháp phần tử hữu hạn (FEM) hoặc sai phân hữu hạn (FDM) truyền thống khi giải phương trình đạo hàm riêng đa chiều trên miền hình học phức tạp đòi hỏi số bậc tự do tăng theo hàm mũ $M^N$ (với $N$ là số chiều không gian, $M$ là số nút lưới mỗi chiều; ví dụ $100^3 = 10^6$ biến cho bài toán 3D). Khi tích hợp phương pháp biên nhúng (IBM), việc giải hệ phương trình Navier–Stokes và phương trình Poisson áp suất trên toàn miền 2D/3D vẫn tạo ra các ma trận đại số khổng lồ, đòi hỏi bộ nhớ và thời gian tính toán rất lớn. Luận án đặt ra hệ thống câu hỏi nghiên cứu và giả thuyết khoa học cốt lõi:

  1. Câu hỏi nghiên cứu 1 (RQ1): Làm thế nào để giải các phương trình vi phân đạo hàm riêng bậc cao (Poisson, Biharmonic) và hệ Navier–Stokes phi tuyến mà không chịu chi phí tính toán cấp số nhân theo chiều không gian?
  2. Câu hỏi nghiên cứu 2 (RQ2): Cơ chế kết hợp nào giữa kỹ thuật lực cưỡng bức trực tiếp của IBM và thuật toán tách biến PGD đảm bảo điều kiện biên không trượt trên bề mặt vật cản mà vẫn bảo toàn tính phân ly một chiều của toán tử vi phân?
  3. Giả thuyết khoa học (H1 & H2): Áp dụng mô hình giảm bậc PGD cho phép quy đổi bài toán $N$ chiều về chuỗi các bài toán 1D giải tích cực nhanh với số bậc tự do chỉ còn $Q \cdot N \cdot M$ ($Q$ là số mode làm giàu nghiệm, $Q \ll M$), từ đó giảm thời gian tính toán tới hàng chục lần so với FDM truyền thống trong khi vẫn giữ vững độ chính xác bậc hai và bậc bốn.

Khung lý thuyết xây dựng dựa trên việc tích hợp Lý thuyết cơ học môi trường liên tục, Thủy khí động lực học chất lưu nhớt không nén (Navier–Stokes), Phương pháp biên nhúng của Peskin (1977) và Mô hình giảm bậc tiên nghiệm Proper Generalized Decomposition (Ammar et al., 2006, 2007). Phạm vi nghiên cứu bao quát các bài toán 2D dòng chảy tầng nhớt không nén ở dải số Reynolds thấp đến trung bình ($Re = 20 - 5000$), kiểm chứng từ bài toán cơ sở (hốc trượt Lid-driven cavity, dòng chảy qua bậc ngược Backward-facing step) đến tương tác phức tạp với vật thể biên cứng đứng yên/dao động và vật thể biên đàn hồi (sợi đàn hồi đơn, đôi và màng kín).

Literature Review và Positioning

Nghiên cứu về FSI trong lịch sử chia thành hai dòng tiếp cận chính: dòng phương pháp lưới tương thích (Body-fitted/Arbitrary Lagrangian-Eulerian - ALE) và dòng phương pháp lưới cố định (Fixed-grid/Immersed Boundary). Tiếp cận ALE truyền thống (Donea et al., 1982; Hughes et al., 1981) giải quyết chính xác điều kiện biên tại bề mặt phân cách nhưng chịu gánh nặng tái tạo lưới liên tục khi vật thể dao động lớn, dẫn đến suy thoái chất lượng phần tử và chi phí xử lý hình học chiếm tới 70-80% thời gian tính toán.

Để khắc phục điều này, Peskin (1977, 2002) khởi xướng phương pháp biên nhúng (IBM) trong mô phỏng van tim, sử dụng hàm phân bố Dirac delta rời rạc để truyền lực giữa lưới cấu trúc Lagrange và lưới lưu chất Euler. Từ nền tảng này, các trường phái IBM phát triển mạnh mẽ:

  • Trường phái Cưỡng bức liên tục (Continuous forcing) của Beyer & LeVeque (1992), Lai & Peskin (2000) dùng lò xo ảo, hay hồi tiếp phản hồi của Goldstein et al. (1993), Saiki & Biringen (1996) thường gặp hiện tượng dao động lực phi vật lý và giới hạn bước thời gian khắt khe.
  • Trường phái Cưỡng bức rời rạc (Discrete/Direct forcing) do Mohd-Yosuf (1997), Fadlun et al. (2000) và Uhlmann (2005) phát triển, tính toán trực tiếp lực biên nhằm thỏa mãn điều kiện biên không trượt tại lưới Euler.
  • Các biến thể cải tiến như Phương pháp mặt phân cách nhúng (Immersed Interface Method - IIM) của LeVeque & Li (1994, 2003), Phương pháp ô ảo (Ghost-Cell IBM) của Tseng & Ferziger (2003), Mittal et al. (2008), và Phương pháp cắt ô (Cut-Cell IBM) của Ye et al. (1999), Udaykumar et al. (2001) duy trì độ chính xác bậc hai tại biên nhưng gặp khó khăn nghiêm trọng khi mở rộng lên không gian 3D do hình học đa giác phức tạp.

Về mặt giải thuật số, các kỹ thuật giảm bậc mô hình (Reduced-Order Models - ROM) như Proper Orthogonal Decomposition (POD) đòi hỏi phải có dữ liệu mô phỏng (snapshots) từ trước thông qua các phép tính toàn phần tốn kém. Ngược lại, Proper Generalized Decomposition (PGD) do Ammar, Chinesta và cộng sự (2006, 2007) đề xuất là một phương pháp giảm bậc tiên nghiệm (a priori), tự động xây dựng các cơ sở không gian tách biệt mà không cần dữ liệu quá khứ. Các nghiên cứu quốc tế của Dumon et al. (2011, 2013) từng kết hợp PGD với thể tích hữu hạn và phương pháp phổ cho bài toán dòng chảy trong hốc vuông đơn giản không có vật cản.

Luận án định vị chính xác khoảng trống học thuật chưa từng được giải quyết: Chưa có công trình nào tích hợp thành công PGD tách biến không gian với kỹ thuật biên nhúng IBM trên lưới so le để mô phỏng tương tác FSI phi tuyến của vật thể biên cứng và màng đàn hồi. Nghiên cứu tạo ra bước tiến vượt bậc khi so sánh với công trình của Uhlmann (2005) về chi phí giải hệ phương trình động lượng và vượt trội hơn nghiên cứu của Dumon et al. (2011) nhờ khả năng xử lý hình học biên phức tạp tùy ý.

Đó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 của trường phái giải tích số và cơ học tính toán qua 3 đóng góp bản chất:

  1. Mở rộng lý thuyết tách biến PGD cho phương trình vi phân bậc cao: Phát triển công thức biến phân dạng yếu (weak form) và thuật toán lặp điểm bất động luân phiên (Fixed-point alternating direction strategy) cho phương trình Poisson 2D/3D ($\nabla^2 u = f$) và phương trình Biharmonic 3D ($\nabla^4 u = f$). Nghiệm $u(x,y,z)$ được phân rã chính xác thành: $$u(x, y, z) \approx \sum_{i=1}^{N} X_i(x) Y_i(y) Z_i(z)$$ Cho phép chứng minh hội tụ tiệm cận của sai số thặng dư trực giao trên không gian Hilbert $L^2(X) \otimes L^2(Y) \otimes L^2(Z)$.
  2. Khung lý thuyết ghép nối IBM-PGD: Đóng góp một nguyên lý mới trong cơ chế truyền tương tác đa trường. Lực biên nhúng $\mathbf{F}(\mathbf{X},t)$ từ miền kết cấu Lagrange được ánh xạ thành lực khối $\mathbf{f}(\mathbf{x},t)$ trên nền lưới Euler thông qua hàm Dirac delta rời rạc 4 điểm trơn của Peskin (2002), chuyển đổi bài toán tương tác hai pha phức tạp thành bài toán dòng lưu chất đồng nhất có nguồn lực ngoài, hoàn toàn tương thích với toán tử tách biến không gian của PGD.
  3. Mô hình toán học cho tương tác biên đàn hồi: Thiết lập hệ phương trình liên hợp mô tả lực căng kéo (tensile force) $T = K_s \left( \left| \frac{\partial \mathbf{X}}{\partial s} \right| - 1 \right)$ và lực kháng uốn (bending force) $E_b \frac{\partial^4 \mathbf{X}}{\partial s^4}$ dọc theo tọa độ cung Lagrange $s$, đóng góp công thức biến phân ổn định giải quyết bài toán biến dạng lớn của kết cấu mỏng ngâm trong dòng chất lưu.

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

Khung phân tích của luận án kết hợp chặt chẽ 3 lý thuyết nền tảng: Lý thuyết thủy động lực học Navier–Stokes không nén, Phương pháp phân tách phân đoạn thời gian (Fractional Step/Projection Method của Chorin 1968, Temam 1969, Perot 1993) và Phân rã tensor dạng tích suy rộng PGD.

Điểm độc đáo nằm ở chỗ: thay vì giải phương trình Poisson áp suất bằng các thuật toán lặp ma trận thưa tốn kém như BiCGStab hay GMRES trên lưới 2D/3D, luận án áp dụng PGD để phân rã trường áp suất thành chuỗi $p(x,y) \approx \sum_{i=1}^N R_i(x)S_i(y)$. Tại mỗi bước làm giàu nghiệm (enrichment step), các hàm số 1D $R(x)$ và $S(y)$ được giải luân phiên bằng sơ đồ sai phân hữu hạn 5 điểm đối xứng bậc 4 ($O(h^4)$). Điều này biến một hệ phương trình đại số kích thước $(N_x \times N_y) \times (N_x \times N_y)$ thành hai hệ phương trình dạng 5 đường chéo kích thước $N_x \times N_x$ và $N_y \times N_y$, triệt tiêu hoàn toàn sự tắc nghẽn bộ nhớ khi số điểm lưới lên đến hàng nghìn điểm mỗi chiều.

Điều kiện biên hạn định (boundary conditions) được xác lập nghiêm ngặt: miền tính toán Euler sử dụng biên vận tốc Dirichlet/Neumann, trong khi bề mặt vật thể ngâm tuân thủ điều kiện không trượt $\mathbf{u}(\mathbf{X}) = \mathbf{U}_b$ thông qua bước hiệu chỉnh lực cưỡng bức.

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 thực nghiệm tính toán (Computational Positivism & Deterministic Continuum Mechanics). Mô hình thiết kế nghiên cứu là sự kết hợp đa tầng (Multi-level computational design) giữa phương pháp toán tử liên tục và rời rạc hóa số bậc cao:

  • Cấp độ toàn cục (Eulerian Field): Hệ tọa độ Descartes trực giao cố định đại diện cho miền lưu chất.
  • Cấp độ cục bộ (Lagrangian Interface): Hệ tọa độ cong $s$ một chiều di động đại diện cho biên vật thể.
  • Cấp độ khử chiều (PGD Decomposition Level): Chiếu biến phân Galerkin trên các không gian tích vô hướng 1D $L^2$.

Quy trình nghiên cứu rigorous

Quy trình mô phỏng được chuẩn hóa theo sơ đồ thuật toán nghiêm ngặt:

  1. Rời rạc hóa không gian: Áp dụng lưới so le (Staggered Grid / MAC grid của Harlow & Welch 1965), trong đó biến áp suất $p$ và lực khối $f$ được định vị tại tâm ô lưới $(i,j)$, còn các thành phần vận tốc $u$ và $v$ được đặt tương ứng tại trung điểm các cạnh thẳng đứng $(i+1/2, j)$ và nằm ngang $(i, j+1/2)$. Cách bố trí này ngăn chặn hoàn toàn hiện tượng dao động áp suất răng cưa (checkerboard oscillations).
  2. Xấp xỉ đạo hàm bậc cao: Sử dụng sơ đồ sai phân hữu hạn 5 điểm chính xác bậc bốn ($O(h^4)$) cho các đạo hàm bậc nhất, bậc hai và bậc bốn:
    • Điểm nút bên trong ($i = 3, \dots, N-2$): $$\left. \frac{d^2 f}{dx^2} \right|i = \frac{-\frac{1}{12} f{i-2} + \frac{4}{3} f_{i-1} - \frac{5}{2} f_i + \frac{4}{3} f_{i+1} - \frac{1}{12} f_{i+2}}{h^2} + O(h^4)$$
    • Điểm nút cận biên ($i = 1, 2, N-1, N$): Áp dụng sơ đồ sai phân lệch tâm tiến/lùi bậc 4 tối ưu để duy trì tính nhất quán toàn cục.
  3. Tiêu chuẩn hội tụ kép: Luận án áp dụng đồng thời hai tiêu chuẩn dừng:
    • Hội tụ lặp nội tại cho mode PGD: $\epsilon = | R^{(q)} S^{(q)} - R^{(q-1)} S^{(q-1)} |_2 \le 10^{-8}$.
    • Hội tụ toàn cục nghiệm PDE: $E = | \nabla^2 u - f |_2 \le 10^{-6}$.

Data và phân tích

Toàn bộ mã nguồn thuật toán và mô hình số được tự lập trình phát triển độc lập trên nền tảng MATLAB (phiên bản 2016a) chạy trên phần cứng CPU Intel Core i7-3632QM @ 2.20GHz. Quá trình kiểm định sai số và bậc hội tụ số học được thực hiện dựa trên chuẩn sai số không gian: $$e_u = \max_{i,j} |u_{i,j} - u_{i,j}^{ref}|, \quad e_{u,L2} = \sqrt{\frac{1}{N_x N_y} \sum_{i=1}^{N_x} \sum_{j=1}^{N_y} (u_{i,j} - u_{i,j}^{ref})^2}$$ Bậc hội tụ không gian được tính toán tự động: $$\text{Order} = \frac{\log(e_{u1} / e_{u2})}{\log(h_1 / h_2)}$$ Kết quả kiểm chứng nghiệm giải tích cho thấy thuật toán đạt tiệm cận bậc hội tụ lý thuyết bậc 4 cho phương trình đạo hàm riêng và bậc 2 cho hệ Navier–Stokes kết hợp IBM.

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

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

Luận án mang lại 5 kết quả đột phá mang tính định lượng rõ ràng:

  1. Kiểm chứng phương trình vi phân bậc cao: Giải phương trình Poisson giải tích $\nabla^2 u = -8\pi^2 \sin(2\pi x)\cos(2\pi y)$ trên lưới $100 \times 100$, phương pháp PGD đạt nghiệm cực kỳ chính xác với sai số $e_{u,\infty} = 1.05 \times 10^{-6}$ và $e_{u,L2} = 7.42 \times 10^{-7}$, thời gian giải chỉ mất 0.082 giây. Với phương trình Biharmonic trên tấm mỏng chịu uốn với 4 cạnh ngàm/gối tựa, PGD giải quyết chính xác trong chưa đầy 0.15 giây.
  2. Dòng chảy trong hốc vuông (Lid-driven Cavity): Khảo sát dải Reynolds từ $Re = 100$ đến $Re = 5000$. Tại $Re = 5000$, phương pháp PGD bắt trọn cấu trúc xoáy chính tại tâm $(x_c = 0.5117, y_c = 0.5352)$ cùng 3 cấu trúc xoáy thứ cấp ở các góc đáy và góc trên, khớp hoàn toàn với dữ liệu chuẩn kinh điển của Ghia, Ghia & Shin (1982) với độ lệch vận tốc cực đại dưới 1.2%.
  3. Dòng chảy qua bậc ngược (Backward-facing step): Mô phỏng từ $Re = 100$ đến $Re = 800$. Luận án chứng minh quy luật phi tuyến của chiều dài vùng xoáy tái bám $x_r/S$: tại $Re = 800$, chiều dài xoáy đạt $x_r/S \approx 12.1$, hoàn toàn tương đồng với thực nghiệm của Armaly et al. (1983) và mô phỏng của Gartling (1990).
  4. Dòng chảy qua trụ tròn cố định và dao động:
    • Tại $Re = 20$ và $Re = 40$: Vùng xoáy ổn định đối xứng sau trụ đạt chiều dài tương ứng $L_w/D = 0.93$ và $L_w/D = 2.24$, hệ số cản $C_d = 2.06$ ($Re=20$) và $C_d = 1.54$ ($Re=40$), trùng khớp dữ liệu của Tritton (1959) và Dennis & Chang (1970).
    • Tại $Re = 100$ và $Re = 200$: Nắm bắt chính xác hiện tượng phát tán xoáy Von Kármán bất ổn định; số Strouhal tính được $St = 0.165$ ($Re=100$) và $St = 0.198$ ($Re=200$), biên độ lực nâng $C_l = \pm 0.335$ và $C_l = \pm 0.69$, hoàn toàn phù hợp với công bố của Williamson (1996) và Liu et al. (1998).
    • Với trụ dao động trực tuyến trong lưu chất tĩnh ($KC = 5$), đồ thị phân bố áp suất và lực cản tái hiện chính xác kết quả thực nghiệm của Dütsch et al. (1998).
  5. Tương tác FSI với màng và sợi đàn hồi:
    • Bài toán sợi đàn hồi đơn gắn sau trụ trong dòng chảy: Tần số vẫy và quỹ đạo đầu tự do sợi đàn hồi có khối lượng đạt trạng thái dao động tuần hoàn phi tuyến bền vững, tái hiện chính xác hiện tượng "vẫy cờ" (flapping flag) của Zhang et al. (2000).
    • Bài toán hai sợi đàn hồi song song: Khám phá sự chuyển pha động lực học từ trạng thái dao động cùng pha (in-phase) ở khoảng cách hẹp $d = 0.5$ sang dao động ngược pha (out-of-phase) khi khoảng cách mở rộng $d = 1.5$, do cơ chế giao thoa xoáy thủy động học khe hẹp chi phối.
    • Bài toán màng đàn hồi khép kín: Độ mất mát diện tích (area loss) tại $t = 0.020$s kiểm soát ở mức dưới 0.8%, chứng minh tính bảo toàn khối lượng vượt trội so với các thuật toán IBM truyền thống (vốn mất mát từ 3% đến 5%).

So sánh định lượng hiệu năng và thời gian tính toán

Hiệu quả vượt bậc của mô hình kết hợp IBM-PGD được chứng minh qua bảng tổng hợp so sánh trực tiếp với phương pháp sai phân hữu hạn truyền thống (FDM):

Bài toán kiểm chuẩn Cấu hình lưới Euler / Bậc tự do Thời gian FDM chuẩn Thời gian đề xuất IBM-PGD Tốc độ gia tốc (Speedup) Độ chính xác chuẩn $L_2$
Poisson 2D miền chữ nhật $100 \times 100$ (10,000 DOFs) 1.84 s 0.082 s 22.4 lần $7.42 \times 10^{-7}$
Biharmonic 3D (Tấm mỏng) $100^3$ ($10^6$ DOFs) ~ 320 s 4.15 s 77.1 lần $1.15 \times 10^{-5}$
Hốc trượt Lid Cavity ($Re=1000$) $128 \times 128$ (16,384 DOFs) 85.6 s/chu kỳ 9.2 s/chu kỳ 9.3 lần Khớp Ghia et al. 99%
Trụ tròn cố định ($Re=100$) $512 \times 256$ (131,072 DOFs) 412.0 s/bước 28.4 s/bước 14.5 lần $C_d=1.35, St=0.165$
FSI Sợi đàn hồi trong dòng lưu chất $256 \times 128$ + 50 nút Lagrange 640.5 s/chu kỳ 52.1 s/chu kỳ 12.3 lần Area loss < 0.8%

Implications đa chiều

  • Ý nghĩa phương pháp luận: Khẳng định tính khả thi tuyệt đối của việc tích hợp các phương pháp giảm bậc mô hình tiên nghiệm (ROM/PGD) vào các bài toán biên chuyển động phức tạp, mở ra hướng đi mới cho mô phỏng CFD thời gian thực (Real-time CFD).
  • Ý nghĩa thực tiễn công nghệ: Cho phép các kỹ sư mô phỏng kết cấu cầu nhịp lớn, cánh turbine điện gió và thiết bị y sinh học (van tim nhân tạo) với độ phân giải lưới cực cao mà không cần đến các siêu máy tính chuyên dụng đắt tiền.

Limitations và Future Research

Mặc dù đạt được những bước tiến mang tính đột phá, luận án thẳng thắn chỉ rõ 4 giới hạn nội tại:

  1. Giới hạn hình học tensor của PGD: Phương pháp PGD tách biến không gian phát huy hiệu quả tối đa trên các miền tính toán quy chuẩn (hình hộp chữ nhật, hình vuông Đề-các). Khi mở rộng sang các miền biên cong tổng quát ngoài Euler, việc phân rã tensor đòi hỏi các phép ánh xạ tọa độ hoặc kỹ thuật PGD không gian tham số phức tạp hơn.
  2. Phạm vi thủy động lực học: Nghiên cứu hiện tại mới chỉ tập trung vào dòng chảy tầng nhớt không nén ở hệ số Reynolds thấp và trung bình ($Re \le 5000$). Chưa tích hợp các mô hình rối (RANS, LES hoặc DNS) cho các dòng chảy rối FSI ở số Reynolds cao ($Re > 10^6$).
  3. Mô hình kết cấu 1D/2D: Các bài toán tương tác đàn hồi mới khảo sát trên các sợi và màng đàn hồi 1D/2D ngâm trong dòng chảy phẳng 2D, chưa mở rộng sang kết cấu phi tuyến vật liệu 3D siêu đàn hồi (Hyperelastic) hoặc tương tác dị hướng phức tạp.

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

  • Tích hợp mô hình dòng chảy rối LES (Large Eddy Simulation) vào thuật toán IBM-PGD để mô phỏng tương tác khí đàn hồi cầu treo và nhà cao tầng trong bão.
  • Phát triển PGD tham số hóa (Parametric PGD) nhằm giải quyết bài toán tối ưu hóa hình học khí động học và phân tích độ tin cậy kết cấu FSI theo thời gian thực.
  • Mở rộng thuật toán lên không gian 3 chiều hoàn chỉnh cho bài toán van tim cơ học 3D tương tác với dòng máu phi Newton.

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

  • Tác động học thuật (Academic Impact): Tạo ra cầu nối đột phá giữa hai cộng đồng nghiên cứu: Cơ học tính toán ứng dụng phương pháp biên nhúng (IBM) và Toán học ứng dụng mô hình giảm bậc (PGD). Dự báo tạo tiềm năng trích dẫn lớn trong các tạp chí ISI đầu ngành như Journal of Computational Physics, Computer Methods in Applied Mechanics and Engineering, Computers & Structures.
  • Chuyển đổi công nghiệp (Industry Transformation): Đóng góp trực tiếp giải pháp phần mềm tối ưu cho các ngành công nghiệp mũi nhọn:
    • Năng lượng tái tạo: Phân tích dao động khí đàn hồi mỏi của cánh quạt turbine gió ngoài khơi chịu tải trọng gió và sóng ngẫu nhiên.
    • Kỹ thuật hàng không & Xây dựng công trình: Dự báo sớm hiện tượng flutter khí đàn hồi của mặt cắt dầm cầu dây văng và cánh máy bay không người lái (UAV).
    • Y sinh học (Biomedical Engineering): Mô phỏng chuyển động của hồng cầu trong vi mạch máu và huyết động học qua van tim nhân tạo.

Đố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 thuật toán rõ ràng, chi tiết từ dạng yếu toán học đến các đoạn mã nguồn MATLAB mẫu được đính kèm ở phụ lục luận án, rút ngắn thời gian phát triển code mô phỏng FSI từ hàng năm xuống vài tuần.
  • Kỹ sư R&D công nghiệp: Sở hữu một công cụ tính toán số học tốc độ cao, tiết kiệm tài nguyên phần cứng, có thể tích hợp vào các chu trình thiết kế tối ưu hóa hình học tự động (Design Optimization Loops).
  • Cơ quan quản lý và Tiêu chuẩn xây dựng: Cung cấp cơ sở khoa học định lượng chuẩn xác để hoàn thiện các tiêu chuẩn thiết kế kháng chấn khí động học cho các công trình cao tầng và cầu nhịp lớn vượt biển.

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ì?
Đóng góp độc đáo nhất là việc thiết lập thành công mô hình toán học tích hợp giữa phương pháp biên nhúng trực tiếp (Direct-forcing IBM) và phương pháp phân tách suy rộng chính tắc PGD. Luận án đã mở rộng lý thuyết PGD vốn chỉ dùng cho các phương trình đạo hàm riêng trên miền trơn đơn giản sang giải quyết thành công bài toán FSI phi tuyến có lực kỳ dị biên nhúng, đưa hệ phương trình Poisson áp suất và Navier–Stokes 2D/3D về chuỗi các phương trình vi phân thường 1D giải tích, phá vỡ rào cản bùng nổ số chiều tính toán.

2. Điểm mới về phương pháp luận khi so sánh với các nghiên cứu quốc tế tiền nhiệm?
So với các nghiên cứu PGD cho dòng chảy của Dumon et al. (2011, 2013) vốn bị giới hạn trong miền hốc vuông kín không có vật cản, luận án lần đầu tiên tích hợp lực khối Dirac delta rời rạc của Peskin (2002) và Uhlmann (2005) vào cấu trúc tensor tách biến của PGD. So với các phần mềm FSI truyền thống dùng lưới ALE (ANSYS Fluent, COMSOL), giải pháp đề xuất không cần tái tạo lưới khi biên di chuyển và giảm thời gian giải phương trình áp suất tới 10-20 lần.

3. Phát hiện số học nào gây bất ngờ và có ý nghĩa cơ học sâu sắc nhất?
Phát hiện bất ngờ nhất là tính ổn định và bảo toàn khối lượng vượt trội của màng đàn hồi khép kín trong chất lưu tĩnh: độ mất mát diện tích tính toán sau thời gian $t = 0.020$s chỉ là 0.78% mà không cần bổ sung các thuật toán sửa sai thể tích nhân tạo phức tạp. Đồng thời, luận án phát hiện sự đảo pha dao động (từ cùng pha sang ngược pha) của hai sợi đàn hồi ngâm song song khi thay đổi khoảng cách không gian $d$, phản ánh chính xác hiệu ứng khóa tần số thủy động lực học phi tuyến.

4. Luận án có cung cấp quy trình tái lập (Replication Protocol) minh bạch không?
Hoàn toàn minh bạch. Luận án trình bày chi tiết từng bước toán tử rời rạc sai phân hữu hạn bậc 4, sơ đồ lưới so le MAC, thuật toán lặp luân phương tích phân không gian, và công bố công khai toàn bộ mã nguồn chương trình MATLAB tại Phụ lục 1 (Giải Poisson 2D bằng PGD) và Phụ lục 2 (Mô phỏng dòng chảy nhớt qua trụ tròn bằng IBM-PGD), cho phép cộng đồng khoa học tái lập kết quả một cách tuyệt đối.

5. Lộ trình phát triển 10 năm của hướng nghiên cứu này được vạch ra như thế nào?
Lộ trình 10 năm bao gồm 3 giai đoạn: (1) Mở rộng giải thuật lên không gian 3D song song hóa trên kiến trúc GPU (2020-2023); (2) Tích hợp mô hình dòng chảy rối LES và vật liệu đàn hồi phi tuyến đa lớp (2024-2027); (3) Xây dựng nền tảng phần mềm "Digital Twin FSI" mô phỏng thời gian thực phục vụ y học cá thể hóa tim mạch và thiết kế khí động học hàng không (2028-2030).

Kết luận

  1. Phát triển thành công hệ phương pháp lai đột phá IBM-PGD: Tích hợp xuất sắc phương pháp biên nhúng và mô hình giảm bậc PGD, giải quyết triệt để bài toán FSI dòng chảy nhớt không nén qua vật cản biên cứng và đàn hồi.
  2. Quy đổi hiệu quả số bậc tự do: Chuyển đổi toàn diện việc giải các PDE đa chiều phức tạp về chuỗi các ODE một chiều giải bằng FDM bậc 4, giảm bậc tự do từ $M^N$ xuống $Q \cdot N \cdot M$, giúp tiết kiệm bộ nhớ máy tính và gia tốc tốc độ tính toán gấp từ 10 đến 70 lần so với các phương pháp số truyền thống.
  3. Giải quyết bài toán vật thể di động và biến dạng không cần tái tạo lưới: Sử dụng kỹ thuật lực cưỡng bức trực tiếp và hàm Dirac delta 4 điểm trơn, loại bỏ 100% chi phí remeshing trên nền lưới Descartes so le cố định.
  4. Kiểm chứng toàn diện và chính xác cao: Kiểm định thành công trên 8 bài toán chuẩn kinh điển (Poisson, Biharmonic, Lid-driven cavity, Backward-facing step, Trụ tròn cố định/dao động, Sợi và Màng đàn hồi), đạt độ chính xác tương đương thực nghiệm quốc tế và sai số thể tích dưới 0.8%.
  5. Mở ra 3 nhánh nghiên cứu học thuật mới: Khởi xướng các hướng nghiên cứu sâu về CFD thời gian thực (Real-time CFD), PGD tham số hóa cho tối ưu hóa khí đàn hồi, và mô phỏng tương tác vi lưu - huyết động học đa trường trong y sinh.
  6. Di sản khoa học đo lường được: Công trình của tác giả Lê Quốc Cường cùng nhóm nghiên cứu GACES tại Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh đã cung cấp một bộ công cụ tính toán hoàn chỉnh, mã nguồn mở, đóng góp một giải pháp học thuật vững chắc cho nền cơ học tính toán Việt Nam hội nhập chuẩn mực quốc tế.