Tổng quan về luận án

Mô phỏng tương tác rắn - lỏng (Fluid-Structure Interaction - FSI) là bài toán vật lý đa trường phi tuyến có ý nghĩa thực tiễn cốt lõi trong các ngành hàng không, công trình ven biển, năng lượng tái tạo (turbin gió) và cơ sinh học dòng chảy (van tim, huyết học). Thách thức toán - cơ học lớn nhất của mô phỏng FSI truyền thống nằm ở sự xung đột giữa miền chất lưu Euler và kết cấu di động biến dạng Lagrange. Khi giải quyết bài toán FSI bằng các phương pháp mạng lưới tương thích (Arbitrary Lagrangian-Eulerian - ALE) kết hợp với các kỹ thuật phần tử hữu hạn (FEM) hay thể tích hữu hạn (FVM), chi phí tái tạo lưới (re-meshing) và giải hệ đại số ma trận toàn cục tăng theo quy luật lũy thừa số chiều $O(M^N)$, gây ra hiện tượng bế tắc tính toán (curse of dimensionality).

Luận án tiến sĩ kỹ thuật "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" của NCS. Lê Quốc Cường (ngành Cơ kỹ thuật, Trường Đại học Sư phạm Kỹ thuật TP. Hồ Chí Minh) đã giải quyết triệt để rào cản này. Nghiên cứu xác định khoảng trống học thuật (research gap) trọng yếu: các biến thể của Phương pháp Biên nhúng (Immersed Boundary Method - IBM) từ Peskin (1977) đến Uhlmann (2005) đã loại bỏ được việc chia lại lưới nhờ hàm Dirac delta rời rạc hóa lực cưỡng bức, nhưng vẫn phụ thuộc vào việc giải phương trình vi phân đạo hàm riêng (PDEs) trên lưới Euler kích thước lớn; ngược lại, phương pháp giảm bậc mô hình tách biến PGD (Ammar et al., 2006; Dumon et al., 2011) giải quyết được bài toán $N$ chiều về chuỗi bài toán 1 chiều $O(Q \cdot N \cdot M)$ nhưng mới chỉ giới hạn ở các miền hình học phẳng đơn giản, chưa xử lý được bài toán có biên di chuyển phức tạp.

Nghiên cứu đặt ra 4 câu hỏi và giả thuyết khoa học tương ứng:

  • RQ1: Làm thế nào để xây dựng giải thuật PGD kết hợp sai phân hữu hạn (FDM) giải chính xác các phương trình vi phân đạo hàm riêng bậc cao (Poisson, Biharmonic) đa chiều? (Hypothesis H1: Giải thuật lặp luân phiên điểm cố định PGD hội tụ nghiệm cấp phân rã hữu hạn $N$ hàm 1D với sai số chuẩn $L_2 < 10^{-5}$).
  • RQ2: Làm thế nào để tích hợp PGD vào thuật toán giải hệ phương trình Navier-Stokes nhớt không nén qua phương pháp chiếu (Projection Method)? (Hypothesis H2: Áp dụng PGD giải phương trình Poisson áp suất giảm đáng kể bộ nhớ RAM và CPU so với FDM truyền thống).
  • RQ3: Giải thuật liên hợp IBM-PGD có mô phỏng chính xác đáp ứng thủy động lực của vật thể biên cứng cố định và dao động ở các dải số Reynolds khác nhau không? (Hypothesis H3: Lực cưỡng bức trực tiếp kết hợp PGD tái tạo hoàn hảo xoáy Von Kármán, lực nâng $C_l$, lực cản $C_d$ và tần số Strouhal $St$).
  • RQ4: Cơ chế ghép nối lực Lagrange-Euler trong IBM-PGD có duy trì tính ổn định số đối với biên đàn hồi phi tuyến, sợi uốn và màng khép kín không? (Hypothesis H4: Mô hình biến dạng liên tục đảm bảo bảo toàn diện tích màng và bắt đúng tần số dao động đồng pha/ngược pha của đa sợi đàn hồi).

Phạm vi nghiên cứu tập trung vào dòng chảy nhớt hai chiều (2D) không nén được ở chế độ chảy tầng ($Re \le 5000$). Kết quả đạt được mức giảm thời gian tính toán lên tới $60% - 75%$ so với phương pháp FDM truyền thống trên cùng mật độ lưới mịn, đồng thời triệt tiêu hoàn toàn chi phí tái cấu trúc lưới trong các bài toán biên động.

Literature Review và Positioning

Nghiên cứu tương tác rắn - lỏng đã hình thành ba dòng tư tưởng học thuật chính:

                          CÁC DÒNG TIẾP CẬN TRONG MÔ PHỎNG FSI
                                           │
         ┌─────────────────────────────────┼─────────────────────────────────┐
         ▼                                 ▼                                 ▼
   MẠNG LƯỚI BÁM BIÊN             PHƯƠNG PHÁP BIÊN NHÚNG          GIẢM BẬC MÔ HÌNH (ROM)
      (Conforming ALE)                      (IBM)                        (POD / PGD)
   • Hughes et al. (1981)         • Peskin (1977): Liên tục       • Ammar et al. (2006): PGD
   • Donea et al. (2004)          • Fadlun et al. (2000): Ép buộc • Dumon et al. (2011): Navier-Stokes
   • Hạn chế: Re-meshing lớn,        trực tiếp (Direct Forcing)   • Hạn chế: Chỉ áp dụng trên
     biến dạng lớn dễ suy biến    • Uhlmann (2005): Lực rời rạc     miền chữ nhật không vật cản
         │                                 │                                 │
         └─────────────────────────────────┼─────────────────────────────────┘
                                           ▼
                                LUẬN ÁN ĐỀ XUẤT: IBM - PGD
                    (Khắc phục re-meshing bằng IBM + Tối ưu hóa giải
                       hệ Navier-Stokes đa chiều thành chuỗi 1D bằng PGD)
  1. Dòng nghiên cứu mạng lưới bám biên (Conforming Mesh / ALE): Bắt nguồn từ các công trình kinh điển của Hughes et al. (1981) và Donea et al. (2004), phương pháp ALE đồng bộ lưới tính toán với bề mặt biên vật thể. Tuy nhiên, khi kết cấu dao động biên độ lớn hoặc biến dạng phức tạp, lưới biên bị xoắn vặn nghiêm trọng, đòi hỏi thuật toán tái tạo lưới cục bộ liên tục, làm tăng độ phức tạp thuật toán và chi phí tính toán lên bậc $O(N^3)$ ma trận toàn cục.
  2. Dòng nghiên cứu Phương pháp Biên nhúng (Immersed Boundary Method - IBM): Khởi xướng bởi Peskin (1977) trong mô phỏng van tim cơ học, giải quyết sự hiện diện của vật thể bằng việc đưa lực cưỡng bức $f$ vào phương trình Navier-Stokes trên lưới Euler cố định thông qua hàm Dirac delta rời rạc. Tiếp nối dòng nghiên cứu này, xuất hiện hai trường phái đối nghịch:
    • Cưỡng bức liên tục (Continuous Forcing): Peskin (1977), Lai & Peskin (2000) sử dụng lò xo ảo để tạo liên kết cứng, song đòi hỏi hệ số cứng rất lớn, dẫn đến hệ phương trình vi phân cứng (stiff system) và bước thời gian tích phân cực nhỏ.
    • Cưỡng bức rời rạc (Discrete Forcing / Direct Forcing): Mohd-Yusof (1997), Fadlun et al. (2000), Uhlmann (2005) xác định lực cưỡng bức trực tiếp từ việc hiệu chỉnh vận tốc tại các điểm lưới biên nhúng. Dù khắc phục tính bất ổn định của lò xo ảo, phương pháp này đòi hỏi giải phương trình áp suất Poisson trên toàn bộ lưới Euler mịn với khối lượng ẩn số ma trận khổng lồ.
  3. Dòng nghiên cứu Giảm bậc mô hình (Reduced-Order Modeling - ROM): Phương pháp Proper Orthogonal Decomposition (POD) (Sirovich, 1987; Holmes et al., 1996) yêu cầu phải có dữ liệu mô phỏng từ trước (snapshots). Để vượt qua hạn chế này, Ammar et al. (2006, 2007) và Chinesta et al. (2010) phát triển Proper Generalized Decomposition (PGD), xây dựng nghiệm tách biến dạng tích $u(x,y,z) \approx \sum_{i=1}^N X_i(x)Y_i(y)Z_i(z)$ a priori mà không cần dữ liệu huấn luyện trước. Dumon et al. (2011, 2013) mở rộng PGD cho hệ phương trình Navier-Stokes kết hợp phương pháp thể tích hữu hạn và phương pháp phổ. Tuy nhiên, các công trình quốc tế của Dumon et al. chỉ giải quyết được miền hình học chữ nhật đơn giản thuần túy, bế tắc khi có sự xuất hiện của vật thể chuyển động hoặc biên đàn hồi phức tạp.

So sánh với hai nghiên cứu quốc tế tiêu biểu:

  • Nghiên cứu của Taira & Colonius (2007) (Immersed boundary projection method): Giải quyết liên kết vận tốc - áp suất bằng cách ghép nhân tử Lagrange của lực biên vào hệ phương trình Poisson mở rộng. Tuy đạt độ chính xác cao nhưng kích thước ma trận mở rộng tăng vọt, tiêu tốn dung lượng bộ nhớ lớn và ma trận toán tử bị suy biến khi biên di chuyển.
  • Nghiên cứu của Dumon et al. (2011) (Proper Generalized Decomposition for incompressible Navier-Stokes equations): Ứng dụng PGD giải dòng chảy trong hốc vuông (lid-driven cavity) thành công nhưng không thể thiết lập điều kiện biên cho các vật thể nội tại phức tạp nếu không có cơ chế tích hợp lực thể.

Vị trí học thuật của luận án: Luận án định vị tại điểm giao thoa đột phá: sử dụng IBM để ảo hóa hình học vật thể phức tạp thành trường lực thể $f$ tác động lên nền chất lưu Descartes trực giao, sau đó tận dụng PGD để phân rã hệ toán tử vi phân nhiều chiều thành chuỗi các phương trình vi phân 1 chiều, triệt tiêu bài toán ma trận kích thước lớn.

Đó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ơ học tính toán qua 3 trụ cột:

  1. Mở rộng lý thuyết xấp xỉ biến PGD cho toán tử vi phân bậc cao: Mở rộng lý thuyết của Chinesta và Ammar cho phương trình Poisson ($\nabla^2 u = f$) và phương trình Biharmonic ($\nabla^4 u = f$) trong không gian 2D và 3D. Bằng cách thiết lập dạng biến phân yếu (weak form) kết hợp phép chiếu Galerkin đa chiều, luận án dẫn xuất tường minh hệ phương trình vi phân thường 1D chứa các tích phân hình học liên kết: $$a_y = \int_Y S^2 dy, \quad b_y = \int_Y S \frac{d^2 S}{dy^2} dy, \quad c_y = \int_Y S \frac{d^4 S}{dy^4} dy$$
  2. Hình thành lý thuyết tách biến cho trường áp suất dòng nhớt: Đóng góp mô hình toán học giải phương trình Poisson áp suất trong phương pháp chiếu Chorin-Perot trên lưới so le (staggered grid). Hệ PDEs 2D/3D liên kết phi tuyến được chuyển vị thành chuỗi phương trình 1 chiều có ma trận hệ số tam đường chéo (tridiagonal) hoặc ngũ đường chéo (pentadiagonal), giải trực tiếp qua thuật toán Thomas với độ phức tạp $O(M)$ thay vì $O(M^2)$ hay $O(M^3)$.
  3. Mô hình hóa tương tác phi tuyến biên đàn hồi dạng Euler-Lagrange liên tục: Thiết lập tường minh trường lực kéo $F_\sigma$ và lực uốn $F_b$ trên hệ tọa độ cong cục bộ $s$ của sợi đàn hồi: $$\mathbf{F}(s, t) = \frac{\partial}{\partial s} \left( \sigma(s,t) \boldsymbol{\tau}(s,t) \right) - \frac{\partial^2}{\partial s^2} \left( E_b \frac{\partial^2 \mathbf{X}(s,t)}{\partial s^2} \right)$$ với $\sigma(s,t) = K_s \left( \left| \frac{\partial \mathbf{X}}{\partial s} \right| - 1 \right)$.
   NGUYÊN LÝ KHUNG PHÂN TÍCH LIÊN HỢP IBM - PGD
   
   ┌────────────────────────────────────────────────────────┐
   │ BƯỚC 1: LƯỚI LAGRANGE KẾT CẤU (X_k)                   │
   │ Tính lực biên: Cứng (Direct Forcing) / Đàn hồi (F_s,b)  │
   └─────────────────────────┬──────────────────────────────┘
                             │ Phân phối lực qua Dirac Delta:
                             │ f(x,t) = ∫ F(s,t) δ_h(x - X(s,t)) ds
                             ▼
   ┌────────────────────────────────────────────────────────┐
   │ BƯỚC 2: HỆ PHƯƠNG TRÌNH NAVIER - STOKES (LƯỚI EULER)   │
   │ Vận tốc trung gian: u* = u^n + Δt [ - (u·∇)u + ν∇²u + f]│
   └─────────────────────────┬──────────────────────────────┘
                             │ Vận tốc không thỏa điều kiện ∇·u = 0
                             ▼
   ┌────────────────────────────────────────────────────────┐
   │ BƯỚC 3: GIẢI PHƯƠNG TRÌNH POISSON ÁP SUẤT BẰNG PGD     │
   │ ∇² p^{n+1} = (ρ/Δt) ∇·u*                               │
   │ Tách biến: p(x,y) ≈ ∑ R_i(x) S_i(y)                    │
   │ Giải luân phiên các ODE 1D bằng thuật toán Thomas      │
   └─────────────────────────┬──────────────────────────────┘
                             │ Hiệu chỉnh trường vận tốc và nội suy:
                             │ u^{n+1} = u* - (Δt/ρ) ∇p^{n+1}
                             │ U(X_k) = ∫ u^{n+1} δ_h(x - X_k) dx
                             ▼
   ┌────────────────────────────────────────────────────────┐
   │ BƯỚC 4: CẬP NHẬT TỌA ĐỘ VẬT THỂ BIÊN ĐÀN HỒI           │
   │ X^{n+1} = X^n + Δt · U_b^{n+1}                         │
   └────────────────────────────────────────────────────────┘

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

Khung phân tích của luận án tích hợp 4 mô hình lý thuyết:

  • Lý thuyết Động lực học chất lưu Navier-Stokes: Mô tả chất lưu thực không nén được.
  • Lý thuyết Phân rã trực giao PGD: Thuật toán lặp luân phiên điểm cố định (Fixed-point alternating direction strategy) xác định tuần tự các mode hàm độc lập.
  • Lý thuyết Biên nhúng lực cưỡng bức rời rạc: Truyền tương tác đa trường không cần bám biên.
  • Lý thuyết Cơ học tấm vỏ phi tuyến và sợi đàn hồi Euler-Bernoulli: Xác định lực đàn hồi cục bộ.

Điều kiện biên giới hạn (boundary conditions) được xác lập nghiêm ngặt: chất lưu Newton, không xét quá trình chuyển pha, dòng chảy đẳng nhiệt và tỷ số giãn dài kết cấu tuân thủ định luật Hooke biến dạng uốn tuyến tính cục bộ.

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

Thiết kế nghiên cứu

Nghiên cứu được xây dựng trên lập trường triết học Thực nghiệm tính toán (Computational Positivism / Computational Empiricism). Bản chất của phương pháp là xây dựng mô hình toán - cơ học tiền định, rời rạc hóa bằng giải thuật số và kiểm chứng tính đúng đắn (verification and validation) thông qua đối chuẩn trực tiếp với các nghiệm giải tích chuẩn (analytical solutions) và các tập dữ liệu thực nghiệm chuẩn quốc tế.

Thành phần thiết kế Đặc tả thông số kỹ thuật
Không gian miền tính toán Hình chữ nhật / Hình vuông 2D Cartesian ($[0, L_x] \times [0, L_y]$)
Cơ chế lưới Euler chất lưu Lưới so le (Staggered Grid / MAC mesh)
Cơ chế lưới Lagrange vật thể Chuỗi hạt chất điểm rời rạc cách đều $\Delta s \approx 0.5 \Delta x \sim 1.0 \Delta x$
Kích thước mẫu kiểm định $64 \times 64$, $100 \times 100$, $128 \times 128$, $256 \times 256$, $512 \times 512$ điểm nút lưới
Dải thông số dòng chất lưu Hệ số Reynolds $Re \in [20, 5000]$; Hệ số Keulegan-Carpenter $KC \in [5, 10]$
Ngôn ngữ & Môi trường MATLAB R2018a / Linux Workstation Environment

Quy trình nghiên cứu rigorous

Quy trình nghiên cứu triển khai qua 4 tầng kiểm định đa mức (multi-level validation protocol):

       QUY TRÌNH NGHIÊN CỨU KIỂM ĐỊNH ĐA TẦNG (RIGOROUS VALIDATION)
       
       TẦNG 1: KIỂM CHỨNG TOÁN HỌC THUẦN TÚY (BENCHMARK PDEs)
       ├─ Phương trình Poisson 2D/3D (Nghiệm hàm lượng giác, đa thức)
       └─ Phương trình Biharmonic 2D/3D (Uốn tấm mỏng biên ngàm/gối tựa)
                                    │
                                    ▼
       TẦNG 2: KIỂM CHỨNG DÒNG CHẢY KHÔNG CÓ VẬT CẢN (CFD THUẦN)
       ├─ Dòng chảy hốc vuông Lid-driven cavity (Re = 100 -> 5000)
       └─ Dòng chảy qua bậc Backward-facing step (Re = 100 -> 800)
                                    │
                                    ▼
       TẦNG 3: KIỂM CHỨNG TƯƠNG TÁC BIÊN CỨNG (IBM - PGD RIGID)
       ├─ Trụ tròn cố định (Re = 20, 40, 100, 200)
       ├─ Trụ tròn dao động trực tuyến trong chất lưu tĩnh (KC = 5)
       └─ Trụ tròn dao động cắt ngang dòng tự do (Lock-in, fe/f0 = 0.8 -> 1.2)
                                    │
                                    ▼
       TẦNG 4: KIỂM CHỨNG TƯƠNG TÁC BIÊN ĐÀN HỒI (IBM - PGD ELASTIC)
       ├─ Một sợi đàn hồi không khối lượng / có khối lượng
       ├─ Hai sợi đàn hồi song song (Khảo sát khoảng cách khe hở d)
       └─ Màng đàn hồi khép kín (Đánh giá mức độ bảo toàn diện tích)
  1. Giao thức rời rạc hóa thời gian - không gian:
    • Bước thời gian $\Delta t$ được khống chế bởi điều kiện ổn định Courant-Friedrichs-Lewy: $CFL = \max \left( \frac{u \Delta t}{\Delta x}, \frac{v \Delta t}{\Delta y} \right) \le 0.2$.
    • Số hạng phi tuyến đối lưu $(\mathbf{u}\cdot\nabla)\mathbf{u}$ được rời rạc bằng sơ đồ sai phân trung tâm bậc hai (Central Difference Scheme) kết hợp sơ đồ sai phân ngược (Upwind Scheme) để đảm bảo tính ổn định ở $Re$ cao.
    • Hàm Dirac delta rời rạc 4 điểm của Peskin (2002) được ứng dụng: $$\delta_h(x) = \frac{1}{h} \phi\left(\frac{x}{h}\right)$$ $$\phi(r) = \begin{cases} \frac{1}{8}\left(3 - 2|r| + \sqrt{1 + 4|r| - 4r^2}\right), & |r| \le 1 \ \frac{1}{8}\left(5 - 2|r| - \sqrt{-7 + 12|r| - 4r^2}\right), & 1 < |r| \le 2 \ 0, & |r| > 2 \end{cases}$$
  2. Độ tin cậy và chuẩn sai số:
    • Tiêu chuẩn dừng PGD cục bộ: $\varepsilon_{RST} = |R^{(q)}S^{(q)}T^{(q)} - R^{(q-1)}S^{(q-1)}T^{(q-1)}|_{L_2} \le 10^{-6}$.
    • Tiêu chuẩn dừng PGD toàn cục: $E = |\nabla^2 u - f|_{L_2} \le 10^{-5}$.
    • Sai số chuẩn hóa đánh giá theo chuẩn $L_2$ và $L_\infty$: $E_{L_2} = \sqrt{\frac{\sum_{i,j} (u_{i,j}^{PGD} - u_{i,j}^{ex})^2}{\sum_{i,j} (u_{i,j}^{ex})^2}}$.

Data và phân tích

Dữ liệu mô phỏng thu thập từ các trường biến số tức thời (vận tốc $u, v$, áp suất $p$, độ xoáy $\omega = \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y}$) và các đại lượng khí động lực tích phân:

  • Lực cản $F_x$ và hệ số cản: $C_d = \frac{2 F_x}{\rho U_\infty^2 D}$
  • Lực nâng $F_y$ và hệ số nâng: $C_l = \frac{2 F_y}{\rho U_\infty^2 D}$
  • Số Strouhal không thứ nguyên: $St = \frac{f_s D}{U_\infty}$ (với $f_s$ là tần số nhả xoáy xác định qua biến đổi Fourier nhanh FFT).

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

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

  1. Khả năng triệt tiêu sự bùng nổ ẩn số của PGD trên phương trình vi phân bậc cao:
    • Trong phương trình Poisson 3D với lưới $100 \times 100 \times 100$ ($10^6$ bậc tự do), PGD chỉ mất 6 đến 8 cặp mode phân rã ($N=8$) để đạt độ chính xác sai số $L_2 = 3.24 \times 10^{-5}$.
    • Đối với phương trình Biharmonic uốn tấm mỏng 2D lưới $100 \times 100$, PGD tái tạo độ võng với thời gian chỉ mất 0.84 giây, đạt độ khớp $99.98%$ so với nghiệm giải tích Navier và Timoshenko.
┌────────────────────────────────────────────────────────────────────────────┐
│ KẾT QUẢ ĐỐI CHUẨN ĐỘ DÀI VÙNG XOÁY (L_r/D) VÀ HỆ SỐ CẢN C_d Ở Re THẤP     │
├──────────────┬───────────────────┬───────────────────┬─────────────────────┤
│ Thông số     │ Re = 20 (Nghiên   │ Re = 40 (Nghiên   │ Dữ liệu đối sánh    │
│              │ cứu Luận án)      │ cứu Luận án)      │ Quốc tế             │
├──────────────┼───────────────────┼───────────────────┼─────────────────────┤
│ L_r / D      │ 0.92              │ 2.24              │ 0.91 & 2.25 (Tritton│
│              │                   │                   │ 1959; Dennis 1970)  │
│ C_d          │ 2.06              │ 1.54              │ 2.05 & 1.52 (He &   │
│              │                   │                   │ Doolen, 1997)       │
└──────────────┴───────────────────┴───────────────────┴─────────────────────┘
  1. Tái tạo chính xác cấu trúc xoáy nhả bất ổn định và hiện tượng Lock-in:
    • Dòng qua trụ tròn tại $Re = 100$: Luận án thu được $C_d = 1.34 \pm 0.012$, $C_l = \pm 0.315$, và số Strouhal $St = 0.165$. Tại $Re = 200$: $C_d = 1.36 \pm 0.045$, $C_l = \pm 0.685$, $St = 0.196$, hoàn toàn trùng khớp với kết quả thực nghiệm của Williamson (1996) và mô phỏng của Liu et al. (1998).
    • Trụ tròn dao động cắt ngang ($Re=185$): Mô phỏng thành công hiện tượng khóa tần số (lock-in) khi tỷ số tần số cưỡng bức tiệm cận tần số nhả xoáy tự nhiên $f_e/f_0 \approx 1.0$, biên độ lực nâng $C_l$ tăng vọt gấp 3.2 lần so với trạng thái trụ đứng yên.
  CÁC TRẠNG THÁI DAO ĐỘNG CỦA CẶP SỢI ĐÀN HỒI (HÀNG MÃ MÔ PHỎNG CHƯƠNG 5)
  
  Trường hợp 1: Khoảng cách hẹp d = 0.3 (Hiện tượng Đập lệch pha/Đồng pha bất đối xứng)
  Sợi 1: ~~~~~~~~~~~~~~~~~~ (Dao động lệch pha biên độ cực đại)
  Sợi 2: ~~~~~~~~~~~~~~~~~~ (Trường xoáy tương tác mạnh, áp suất giao thoa)
  
  Trường hợp 2: Khoảng cách rộng d = 0.8 (Hiện tượng Tách biệt động lực học độc lập)
  Sợi 1: ------------------ (Dao động tuần hoàn điều hòa đơn lẻ)
  Sợi 2: ------------------ (Xoáy Von Kármán độc lập phía sau mỗi đuôi sợi)
  1. Hiện tượng ghép cặp phi tuyến của hai sợi đàn hồi đặt song song trong dòng nhớt:
    • Khi khoảng cách giữa hai sợi $d = 0.3$, lực thủy động tạo ra sự tương tác giao thoa luồng xoáy, hai sợi dao động lệch pha (out-of-phase flapping) với tần số biến thiên phức tạp.
    • Khi mở rộng khoảng cách $d \ge 0.8$, hiệu ứng tương tác thủy động suy giảm, hai sợi chuyển sang trạng thái dao động độc lập và dần chuyển sang đồng pha (in-phase flapping).
  2. Bảo toàn diện tích màng đàn hồi khép kín:
    • Mô phỏng màng đàn hồi hình elip co giãn về trạng thái cân bằng tròn trong lưu chất tĩnh đạt mức mất mát diện tích cực thấp: tại thời điểm $t = 0.020\text{ s}$, diện tích ban đầu $S_0 = 0.1963$, diện tích sau biến dạng $S = 0.1948$, tỷ lệ mất mát chỉ $\Delta S/S_0 = 0.76%$, vượt trội so với mức suy hao $2.1%$ của các phương pháp IBM cổ điển dùng hàm liên kết lò xo tuyến tính.

Implications đa chiều

  • Về mặt lý thuyết: Khẳng định tính đúng đắn của việc kết hợp biểu diễn trường Euler không liên tục với các nghiệm suy rộng PGD; mở ra hướng chứng minh toán học giải tích cho tính hội tụ của PGD trên hệ PDEs ghép nối không thuần nhất.
  • Về mặt phương pháp học: Tạo ra bộ khung thuật toán mở (open framework) cho phép tích hợp linh hoạt các module vật lý khác như truyền nhiệt đối lưu hoặc điện từ thủy động lực học (MHD) mà không làm tăng bậc độ phức tạp tính toán.
  • Về mặt ứng dụng thực tiễn: Cung cấp công cụ tối ưu hóa hình dáng cánh quạt turbine gió cỡ nhỏ, phân tích hiện tượng rung do nhả xoáy (Vortex-Induced Vibrations - VIV) trên các đường ống dẫn dầu ngầm biển sâu.
  • Về mặt chính sách và R&D công nghiệp: Giúp các trung tâm nghiên cứu cơ khí - hàng không trong nước làm chủ công nghệ phần mềm mô phỏng độc lập (in-house solver code), giảm phụ thuộc vào các gói phần mềm thương mại đắt đỏ như ANSYS Fluent hay COMSOL Multiphysics.

Limitations và Future Research

  1. Giới hạn về số chiều và chế độ dòng: Mô hình hiện tại mới triển khai thuật toán ghép nối IBM-PGD trên không gian 2D, chưa mở rộng cho bài toán FSI 3D hoàn chỉnh do độ phức tạp trong việc xây dựng toán tử uốn bề mặt màng đàn hồi 3D.
  2. Giới hạn số Reynolds: Nghiên cứu giới hạn ở chế độ dòng chảy tầng ($Re \le 5000$). Chưa tích hợp mô hình mô phỏng xoáy lớn (Large Eddy Simulation - LES) hoặc Reynolds-Averaged Navier-Stokes (RANS) để nắm bắt hoàn toàn hiện tượng rối ở $Re > 10^5$.
  3. Biến dạng kết cấu phi tuyến lớn: Giả định sợi đàn hồi tuân theo mô hình Euler-Bernoulli; chưa xét đến bài toán kết cấu siêu đàn hồi phi tuyến hình học cấp cao (hyperelastic material) hoặc hiện tượng nứt gãy mỏi do dòng va đập.

Định hướng phát triển 5 - 10 năm tới:

  • Nâng cấp bộ giải lên không gian ba chiều (3D) kết hợp kiến trúc tính toán song song đa luồng trên GPU (CUDA C++/OpenCL).
  • Tích hợp mô hình rối LES-PGD để mô phỏng dòng chảy qua turbine gió thực tế.
  • Ứng dụng mô phỏng huyết học chính xác: Tương tác giữa hồng cầu, huyết tương và van tim sinh học hai lá nhân tạo.

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

  • Học thuật: Tạo tiền đề công bố trên các tạp chí quốc tế Q1 chuyên ngành cơ học tính toán (Journal of Computational Physics, Computer Methods in Applied Mechanics and Engineering). Ước tính công trình có tiềm năng thu hút từ 80 đến 120 trích dẫn học thuật quốc tế trong vòng 5 năm đầu công bố.
  • Công nghiệp: Định hình giải pháp mô phỏng nhanh (rapid simulation tool) trong thiết kế thiết bị trao đổi nhiệt dạng tấm uốn và tối ưu khí động lực học cho thiết bị bay không người lái (UAV) hoạt động tại dải Reynolds thấp.
  • Xã hội và Giáo dục: Đóng góp mã nguồn mở giáo khoa (Phụ lục MATLAB) phục vụ đào tạo sau đại học các chuyên ngành Cơ học tính toán, Toán ứng dụng và Khoa học máy tính kỹ thuật.

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

                              CÁC NHÓM THỤ HƯỞNG NGHIÊN CỨU
                                             │
      ┌──────────────────────┬───────────────┴──────────────┬──────────────────────┐
      ▼                      ▼                              ▼                      ▼
  NCS / HỌC GIẢ         GIẢNG VIÊN ĐÀO TẠO            KỸ SƯ R&D CÔNG NGHIỆP    CƠ QUAN QUẢN LÝ /
  KHOA HỌC DỮ LIỆU       ĐẠI HỌC/SAU ĐẠI HỌC          THIẾT KẾ KẾT CẤU          CHÍNH SÁCH
  Khai thác khung        Sử dụng bộ mã nguồn          Ứng dụng thiết kế hệ      Tham chiếu xây dựng
  toán PGD mở rộng       mẫu (MATLAB Phụ lục)         thống tiêu tán năng       tiêu chuẩn thiết kế
  để phát triển mô       làm học liệu giảng dạy       lượng VIV, turbine        khí động học công trình
  hình giảm bậc cho      chuyên đề Cơ học chất        gió và cánh robot         và an toàn thiết bị y
  các bài toán FSI 3D    lưu tính toán (CFD)          mềm dưới nước             tế dòng chảy
  • Nghiên cứu sinh & Nhà nghiên cứu trẻ: Tiếp cận khung phương pháp luận kết hợp ROM và IBM chi tiết để giải các bài toán đa vật lý đa trường (Multiphysics).
  • Kỹ sư R&D Hàng không & Năng lượng: Sử dụng giải thuật PGD để rút ngắn $70%$ chu kỳ mô phỏng thử nghiệm tối ưu kết cấu cánh khí động.
  • Nhà khoa học Y sinh: Áp dụng mô hình màng đàn hồi kín để tính toán áp lực thành mạch và mô phỏng tương tác huyết động học vi mô.

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 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 đổi giảm bậc PGD với phương pháp lực thể biên nhúng Direct Forcing IBM trên lưới so le. Luận án đã mở rộng lý thuyết PGD của Chinesta và Dumon, chứng minh rằng toán tử Poisson áp suất suy biến do lực cưỡng bức biên nhúng hoàn toàn có thể phân rã thành tích các hàm độc lập 1D mà không làm suy giảm độ chính xác nghiệm bảo toàn khối lượng ($\nabla \cdot \mathbf{u} = 0$).
  2. Sự đổi mới về phương pháp luận so với các nghiên cứu trước đây thể hiện như thế nào? So với phương pháp ALE truyền thống (Hughes et al., 1981), nghiên cứu triệt tiêu chi phí $O(N^3)$ của việc chia lại lưới. So với phương pháp IB Chiếu của Taira & Colonius (2007), phương pháp của luận án không tạo ra ma trận hệ số áp suất mở rộng khổng lồ. So với PGD thuần túy của Dumon et al. (2011), luận án lần đầu tiên cho phép PGD xử lý được biên hình học tùy ý, biên dao động cưỡng bức và biên đàn hồi tự do.
  3. Phát hiện bất ngờ và thú vị nhất thu được từ dữ liệu số là gì? Sự chuyển đổi đột ngột trạng thái pha dao động của hai sợi đàn hồi đặt song song trong dòng chảy: tồn tại một ngưỡng khoảng cách tới hạn $d_c \approx 0.5 - 0.6$. Khi $d < d_c$, áp suất dòng khe hẹp tạo lực hút đẩy phi tuyến khiến hai sợi dao động hoàn toàn lệch pha; khi $d \ge d_c$, hệ thống tự tổ chức đồng bộ hóa trường xoáy chuyển sang dao động đồng pha hoàn toàn.
  4. Luận án có cung cấp giao thức tái lập (Replication Protocol) minh bạch không? Hoàn toàn minh bạch và chi tiết. Luận án cung cấp toàn văn các sơ đồ thuật toán tổng quát (Hình 3.3, 4.4, 5.3) và công khai toàn bộ mã nguồn chương trình tính toán nguyên bản viết trên môi trường MATLAB tại Phụ lục 1 (Giải phương trình Poisson 2D bằng PGD) và Phụ lục 2 (Mô phỏng dòng chảy qua trụ tròn bằng IBM-PGD), cho phép cộng đồng học thuật kiểm chứng và phát triển tiếp nối độc lập.
  5. Chương trình nghiên cứu 10 năm tiếp theo được định hình ra sao? Lộ trình 10 năm bao gồm: Giai đoạn 1 (1-3 năm): Hoàn thiện thư viện PGD-FSI 3D hướng đối tượng mã nguồn mở; Giai đoạn 2 (3-6 năm): Mở rộng giải quyết bài toán dòng siêu âm / nén được kết hợp truyền nhiệt; Giai đoạn 3 (6-10 năm): Ứng dụng PGD xây dựng bản sao số (Digital Twin) thời gian thực phục vụ phẫu thuật cấy ghép van tim nhân tạo.

Kết luận

  1. Phát triển thành công khung tính toán liên hợp IBM - PGD: Thiết lập giải pháp đột phá kết hợp sức mạnh xử lý biên phức tạp của IBM với hiệu năng giảm bậc mô hình siêu việt của PGD.
  2. Giải quyết triệt để bài toán phương trình vi phân đạo hàm riêng bậc cao: Dẫn xuất hoàn chỉnh hệ phương trình phân rã 1D cho phương trình Poisson và Biharmonic 2D/3D với độ hội tụ nghiệm cao ($L_2 < 10^{-5}$).
  3. Tối ưu hóa hiệu năng tính toán hệ Navier - Stokes: Giảm $60% - 75%$ thời gian tính toán và tiết kiệm bộ nhớ RAM vượt trội so với các phương pháp sai phân/thể tích hữu hạn truyền thống trên cùng độ mịn lưới.
  4. Mô phỏng chính xác đáp ứng động lực học dòng qua vật cản biên cứng: Tái hiện chuẩn xác quy luật nhả xoáy Von Kármán, lực nâng, lực cản và hiện tượng khóa tần số dao động (lock-in) tại các dải $Re = 20 \div 200$.
  5. Làm chủ bài toán tương tác chất lưu - kết cấu biên đàn hồi phi tuyến: Nắm bắt trọn vẹn động lực học dao động phức tạp của sợi đơn, hiện tượng tương tác ghép cặp của đa sợi song song và duy trì bảo toàn diện tích màng khép kín với sai số dưới $1%$.
  6. Mở ra 3 dòng nghiên cứu khoa học mới: Mô phỏng FSI 3D thời gian thực dựa trên PGD; Cơ học dòng chảy sinh học vi mô; và Tối ưu hóa khí động học công trình giảm tiêu hao năng lượng. Nghiên cứu tạo dấu ấn khoa học vững chắc, khẳng định năng lực phát triển các công cụ tính toán cơ học tiên tiến ngang tầm chuẩn mực quốc tế.