Bối cảnh và vấn đề nghiên cứu

Trong cơ học chất lưu và kỹ thuật khí động học, việc mô tả và dự báo các đặc tính chuyển động của dòng chất lưu (chất lỏng và chất khí) đóng vai trò then chốt trong nhiều ứng dụng thực tiễn, từ tính toán lực nâng, mô-men tác động lên thân máy bay, xác định lưu lượng vận chuyển dầu mỏ trong đường ống dẫn đến các mô hình dự báo thời tiết khí tượng. Về mặt toán học, các quá trình này được mô hình hóa thông qua các đặc trưng vận tốc, mật độ và nhiệt độ phụ thuộc vào biến không gian và thời gian.

Hệ phương trình Euler một chiều dạng đẳng entropy (isentropic Euler equations) là mô hình thu gọn của dòng chảy chất lưu trong ống dẫn, được thiết lập từ các định luật bảo toàn khối lượng, động lượng và năng lượng. Khi giả thiết đại lượng entropy $S$ (đo lường mức độ rối loạn trong hệ) là một hằng số đồng nhất trên toàn miền khảo sát, phương trình bảo toàn năng lượng được rút gọn, đưa hệ về dạng hai phương trình đạo hàm riêng phi tuyến:

$$\begin{cases} \rho_t + (\rho v)_x = 0 \ (\rho v)_t + (\rho v^2 + p)_x = 0 \end{cases}, \quad x \in \mathbb{R}, ; t > 0$$

trong đó $\rho(x,t)$ là mật độ chất lưu, $v(x,t)$ là vận tốc dòng chảy, và $p(x,t)$ là áp suất thủy tĩnh tuân theo phương trình trạng thái khí lý tưởng $p = \kappa \rho^\gamma$.

Hệ phương trình Euler đẳng entropy có thể biểu diễn dưới dạng hệ định luật bảo toàn hyperbolic phi tuyến $u_t + f(u)_x = 0$. Một đặc trưng cơ bản của hệ này là sự hình thành các gián đoạn (như sóng sốc) sau một khoảng thời gian hữu hạn ngay cả khi điều kiện ban đầu hoàn toàn trơn, dẫn đến sự không tồn tại của nghiệm cổ điển toàn cục. Nhằm giải quyết bài toán này, khóa luận tập trung xây dựng phương pháp sai phân hữu hạn dạng bảo toàn để tìm nghiệm xấp xỉ cho hệ phương trình.

Mục đích nghiên cứu của đề tài bao gồm 3 nhiệm vụ cụ thể được đánh số trong văn bản:

  1. Thiết lập mô hình hệ phương trình Euler đẳng entropy từ các định luật bảo toàn khối lượng, bảo toàn động lượng và bảo toàn năng lượng.
  2. Xây dựng lược đồ sai phân hữu hạn cân bằng Lax-Friedrichs cho hệ phương trình Euler đẳng entropy.
  3. Kiểm định tính hiệu quả của phương pháp thông qua tính toán độ chênh lệch của nghiệm xấp xỉ và nghiệm chính xác bằng phần mềm Matlab.

Đối tượng nghiên cứu là hệ phương trình Euler đẳng entropy một chiều, bài toán Cauchy và bài toán Riemann tương ứng. Phạm vi nghiên cứu tập trung vào việc khảo sát mô hình lý thuyết hệ hyperbolic một chiều và tìm nghiệm xấp xỉ số cho các cấu hình sóng đặc trưng (sóng tĩnh, sóng giãn, sóng sốc) thông qua lược đồ sai phân Lax-Friedrichs.

Cơ sở lý thuyết và phương pháp

Khóa luận dựa trên nền tảng lý thuyết phương trình đạo hàm riêng hyperbolic, lý thuyết hệ các định luật bảo toàn và giải tích số cho phương trình vi phân:

  • Khung lý thuyết hệ hyperbolic các định luật bảo toàn: Xét hệ $u_t + f(u)_x = 0$ với $u \in \Omega \subset \mathbb{R}^m$ và ma trận Jacobi $A(u) = \frac{\partial f}{\partial u}$. Hệ được gọi là hyperbolic ngặt khi $A(u)$ có $m$ giá trị riêng thực phân biệt $\lambda_1(u) < \lambda_2(u) < \dots < \lambda_m(u)$ cùng hệ véctơ riêng tương ứng độc lập tuyến tính. Đối với hệ Euler đẳng entropy với véctơ trạng thái $u = (\rho, v)^T$, hai giá trị riêng là $\lambda_1 = v - c$ và $\lambda_2 = v + c$, trong đó $c = \sqrt{p'(\rho)} = \sqrt{\kappa \gamma \rho^{\gamma-1}}$ là vận tốc âm địa phương.
  • Lý thuyết nghiệm yếu và điều kiện entropy: Khi nghiệm cổ điển mất tính khả vi, khái niệm nghiệm yếu được định nghĩa qua không gian hàm thử khả vi vô hạn có giá mang compact $C_0^\infty$. Tại các mặt gián đoạn, nghiệm yếu phải thỏa mãn điều kiện bước nhảy Rankine-Hugoniot: $-s[u] + [f(u)] = 0$, với $s$ là tốc độ lan truyền gián đoạn. Để khắc phục tính không duy nhất của nghiệm yếu, tiêu chuẩn nghiệm entropy được áp dụng thông qua cặp hàm lồi entropy - thông lượng entropy $(U, F)$ thỏa mãn bất đẳng thức $U(u)_t + F(u)_x \le 0$ theo nghĩa phân bố.
  • Cấu trúc nghiệm bài toán Riemann: Nghiệm tự đồng dạng $u(x,t) = v(x/t)$ được phân tích thành các dạng sóng cơ bản: sóng $k$-giãn (rarefaction wave) tương ứng với trường đặc trưng phi tuyến thực sự $\nabla \lambda_k \cdot r_k \neq 0$, sóng $k$-sốc (shock wave) thỏa mãn bất đẳng thức Lax $\lambda_k(u_R) < s < \lambda_k(u_L)$, và sóng $k$-tiếp xúc gián đoạn tương ứng với trường suy biến tuyến tính.
  • Lý thuyết hội tụ Lax-Wendroff: Định lý khẳng định nếu một lược đồ sai phân dạng bảo toàn và tương thích hội tụ bị chặn theo biến phân toàn phần ($TV \le R$) khi bước lưới tiến về 0, thì hàm giới hạn chính là một nghiệm yếu của hệ định luật bảo toàn.

Phương pháp nghiên cứu kết hợp giữa nghiên cứu lý thuyết giải tích và thực nghiệm tính toán số:

  • Phương pháp phân tích - tổng hợp: Kế thừa và hệ thống hóa các kết quả từ tài liệu chuyên khảo về giải tích số cho hệ hyperbolic.
  • Phương pháp giải tích chính xác: Thiết lập phương trình đại số phi tuyến xác định đường cong trạng thái sóng tĩnh $g(\rho) = 0$ và sử dụng phương pháp chia đôi (bisection method) để tìm nghiệm phân ly với sai số nhỏ hơn $\varepsilon$.
  • Phương pháp sai phân hữu hạn dạng bảo toàn: Rời rạc hóa miền tính toán bằng lưới sai phân $(x_j, t_n)$ với bước không gian $h = \Delta x$ và bước thời gian $k = \Delta t$, sử dụng hàm thông lượng số Lax-Friedrichs: $$F(u, v) = \frac{1}{2}[f(u) + f(v)] - \frac{\Delta x}{2\Delta t}(v - u)$$
  • Mô phỏng số: Sử dụng ngôn ngữ lập trình Matlab để cài đặt thuật toán kiểm soát bước thời gian tự thích nghi theo điều kiện Courant-Friedrichs-Lewy (CFL).

Nguồn dữ liệu và thông số mô hình sử dụng phương trình trạng thái khí lý tưởng $p(\rho) = \kappa \rho^\gamma$ với các hằng số tiêu chuẩn $\kappa = 1$ và chỉ số đoạn nhiệt $\gamma = 1.4$.

Nội dung chính theo từng chương

Chương 1: Những kiến thức chuẩn bị

Chương 1 trình bày hệ thống cơ sở toán học về hệ phương trình hyperbolic một chiều và xây dựng nghiệm giải tích cho bài toán Riemann của hệ Euler đẳng entropy:

  • Phân tích tính chất hyperbolic và phương pháp đường đặc trưng $dx/dt = a(u)$, minh họa hiện tượng giao cắt của các đường đặc trưng dẫn đến sự hình thành sóng sốc trong phương trình Burgers $u_t + u u_x = 0$ tại thời điểm $t = 1$.
  • Định nghĩa không gian nghiệm yếu, thiết lập điều kiện bước nhảy Rankine-Hugoniot và chứng minh tính không duy nhất của nghiệm yếu qua bài toán Riemann của phương trình Burgers; từ đó đưa ra điều kiện cặp entropy $(U, F)$ nhằm xác định nghiệm vật lý duy nhất.
  • Thiết lập mô hình hệ Euler từ các định luật bảo toàn liên tục, phân tích phương trình trạng thái khí lý tưởng $p = R\rho T$, nội năng riêng $e = c_v T$, entanpi $h = c_p T$, và rút gọn hệ khi entropy $S = c_v \ln(p/\rho^\gamma) + \text{const}$ không đổi.
  • Khảo sát chi tiết ba dạng cấu hình sóng cho hệ Euler đẳng entropy:
    1. Sóng tĩnh: Nghiệm dừng thỏa mãn $\rho v = \rho_0 v_0$ và $\rho v^2 + p = \rho_0 v_0^2 + p_0$. Biến đổi đưa về khảo sát hàm số phi tuyến: $$g(\rho) = \kappa \rho^{\gamma+1} - (\rho_0 v_0^2 + \kappa \rho_0^\gamma)\rho + \rho_0^2 v_0^2 = 0$$ Tác giả chứng minh điều kiện có nghiệm phân biệt phụ thuộc vào giá trị cực tiểu tại $\rho' = \left(\frac{\rho_0 v_0^2 + \kappa \rho_0^\gamma}{\kappa(\gamma+1)}\right)^{1/\gamma}$, từ đó xác lập khoảng phân ly nghiệm $[\rho', \rho_{\max}]$ hoặc $[0, \rho']$ để giải bằng thuật toán chia đôi.
    2. Sóng giãn: Xây dựng đường cong 1-giãn $R_1(u_0)$ và 2-giãn $R_2(u_0)$ liên kết trạng thái bên trái $u_L$ và bên phải $u_R$: $$v = v_0 \mp \frac{2\sqrt{\kappa \gamma}}{\gamma - 1} \left(\rho^{\frac{\gamma-1}{2}} - \rho_0^{\frac{\gamma-1}{2}}\right)$$
    3. Sóng sốc: Áp dụng hệ thức Rankine-Hugoniot kết hợp điều kiện bước nhảy Lax, thiết lập đường cong 1-sốc $S_1(u_0)$ và 2-sốc $S_2(u_0)$: $$v = v_0 \mp \sqrt{\frac{(p - p_0)(\rho - \rho_0)}{\rho \rho_0}}$$

Chương 2: Xây dựng lược đồ sai phân hữu hạn cho hệ hyperbolic các định luật bảo toàn

Chương 2 tập trung vào việc thiết kế và phân tích tính chất toán học của các phương pháp số:

  • Khảo sát bài toán đối lưu tuyến tính $u_t + A u_x = 0$, phân tích sai số chặt cụt địa phương $L_k(x,t)$ bằng khai triển chuỗi Taylor, chứng minh lược đồ Lax-Friedrichs đạt độ chính xác cấp 1 theo thời gian $O(\Delta t)$. Phân tích tính ổn định Lax-Richtmyer và điều kiện CFL với số Courant $\nu = \max |\lambda_p| \frac{\Delta t}{\Delta x} \le 1$.
  • Chỉ ra nguy cơ của việc sử dụng các lược đồ sai phân dạng không bảo toàn (non-conservative schemes) đối với bài toán phi tuyến: phương pháp có thể hội tụ rất tốt về một hàm số nhưng hàm đó lại không phải là nghiệm yếu của bài toán (minh họa qua phương trình Burgers với nghiệm gián đoạn lan truyền sai vận tốc sóng).
  • Định nghĩa dạng sai phân bảo toàn thông qua hàm thông lượng số $F(u, v)$ và chứng minh lược đồ Lax-Friedrichs phi tuyến: $$U_j^{n+1} = \frac{1}{2}(U_{j-1}^n + U_{j+1}^n) - \frac{\Delta t}{2\Delta x}\left[f(U_{j+1}^n) - f(U_{j-1}^n)\right]$$ hoàn toàn viết được dưới dạng bảo toàn tổng quát.
  • Chứng minh tính tương thích của lược đồ thông qua điều kiện liên tục Lipschitz của hàm thông lượng số và trình bày chi tiết chứng minh Định lý Lax-Wendroff: khẳng định nghiệm số $U_h(x,t)$ khi tinh mịn lưới ($\Delta x, \Delta t \to 0$) sẽ hội tụ về dạng tích phân yếu của hệ định luật bảo toàn.
  • Áp dụng lược đồ Lax-Friedrichs cho hệ Euler đẳng entropy, thiết lập công thức xác định bước thời gian tự thích nghi ở mỗi bước lặp nhằm bảo đảm điều kiện CFL: $$\Delta t = \text{CFL} \cdot \frac{\Delta x}{\max_{j} \left(|\lambda_1(U_j^n)|, |\lambda_2(U_j^n)|\right)}$$

Chương 3: Các thử nghiệm số

Chương 3 tiến hành cài đặt chương trình tính toán trên Matlab để kiểm định tính hiệu quả của lược đồ Lax-Friedrichs qua 4 kịch bản thử nghiệm:

  • Thử nghiệm số 1 (Sóng tĩnh):
    • Trường hợp 1 (điều kiện tầm thường): Chọn $u_L = u_R = (1, 2)^T$, nghiệm chính xác là hàm hằng $\rho(x,t) = 1, v(x,t) = 2$. Lược đồ bắt hoàn toàn sóng tĩnh ở mọi độ chia lưới $N = 500, 1000, 4000$.
    • Trường hợp 2 (bước nhảy gián đoạn): Chọn trạng thái $u_L = (1, 2)^T$, tính toán trạng thái bên phải qua phương pháp chia đôi thu được $u_R \approx (2.34298, 0.85387)^T$. Nghiệm đúng có dạng hàm từng khúc gián đoạn tại $x = 0$.
  • Thử nghiệm số 2 (Sóng 1-giãn): Thiết lập tham số $\kappa = 1, \gamma = 1.4$, điều kiện đầu $u_L = (1, 2)^T, \rho_R = 0.5$, tính toán vận tốc theo đường cong 1-giãn thu được $v_R \approx 2.76583$. Nghiệm chính xác liên tục nối hai trạng thái qua vùng sóng dãn.
  • Thử nghiệm số 3 (Sóng 1-sốc): Thiết lập điều kiện đầu $u_L = (1, 2)^T, \rho_R = 1.5$, tính vận tốc theo đường cong 1-sốc thu được $v_R \approx 1.49532$. Vận tốc truyền sóng sốc giải tích xác định bởi $s = \frac{\rho_R v_R - \rho_L v_L}{\rho_R - \rho_L}$.
  • Thử nghiệm số 4 (Khảo sát tham số CFL): Cố định số điểm chia lưới $N = 1000$, thực hiện biến thiên số CFL từ 0.1 đến các giá trị tiệm cận 1 để đánh giá tương quan giữa thời gian tính toán và độ chính xác của nghiệm số.

Kết quả và đóng góp

Kết quả tính toán số

Khóa luận ghi nhận kết quả mô phỏng số chi tiết trên các mức phân giải lưới $N \in {500, 1000, 4000}$:

Kịch bản thử nghiệm Số điểm chia ($N$) Đánh giá sai số và hiệu năng tính toán
Thử nghiệm 1
(Sóng tĩnh gián đoạn)		 $N = 500$
$N = 1000$
$N = 4000$		 - Tại $N = 500$: Sai số vượt phạm vi $10^{-2}$, sóng tĩnh chưa được bắt tốt.
- Tại $N = 1000$ và $N = 4000$: Sai số giảm xuống phạm vi $10^{-3}$, sóng tĩnh được nắm bắt ổn định (thời gian tính toán tại $N=1000$ là $\approx 10.21\text{s}$).
Thử nghiệm 2
(Sóng 1-giãn)		 $N = 500$
$N = 1000$
$N = 4000$		 - Bắt sóng hiệu quả cao nhất tại $N = 4000$ với sai số tương đối đạt $\approx 0.0024$, thời gian tính toán $\approx 155.35\text{s}$.
- Độ hiệu quả bắt sóng 1-giãn tương đồng với trường hợp sóng tĩnh (sai số so lệch giữa hai sóng ở $N=500$ là $\approx 0.0027$).
Thử nghiệm 3
(Sóng 1-sốc)		 $N = 500$
$N = 1000$
$N = 4000$		 - Sóng 1-sốc được bắt tốt, sai số tương đối nằm trong phạm vi $10^{-2}$ (sai số tuyệt đối ghi nhận ở mức $\approx 0.0056$).
Thử nghiệm 4
(Ảnh hưởng số CFL)		 $N = 1000$
$\text{CFL} \in [0.1, 0.9]$		 - Giảm CFL xuống $0.1$ khiến thời gian tính toán tăng gấp 10 lần, đồng thời làm sai số tăng mạnh từ $10^{-2}$ lên $10^{-1}$ do khuếch tán số.
- Tăng số CFL tiệm cận 1 giúp cải thiện đồng thời cả độ chính xác lẫn tốc độ tính toán.

Đóng góp học thuật của công trình

  1. Hệ thống hóa phương pháp luận giải tích và số: Xây dựng quy trình khép kín từ việc chứng minh tính hyperbolic ngặt, thiết lập giải tích các đường cong sóng Riemann đến việc rời rạc hóa bảo toàn bằng lược đồ sai phân.
  2. Chứng minh tính chất toán học của lược đồ: Chứng minh tường minh tính tương thích Lipschitz của hàm thông lượng số Lax-Friedrichs đối với hệ luật bảo toàn phi tuyến và liên hệ trực tiếp với Định lý Lax-Wendroff để đảm bảo tính hội tụ về nghiệm yếu.
  3. Xác lập quy tắc tối ưu hóa tham số tính toán: Khẳng định bằng dữ liệu thực nghiệm rằng trong giới hạn ổn định $\text{CFL} < 1$, việc lựa chọn số CFL càng lớn sẽ giảm thiểu sai số khuếch tán số và tối ưu hóa thời gian thực thi thuật toán.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế: Các thử nghiệm số trong khóa luận chỉ mới tập trung kiểm định trên từng dạng sóng đơn lẻ (sóng tĩnh, sóng 1-giãn, sóng 1-sốc) của bài toán Riemann. Tác giả chưa thực hiện mô phỏng cho bài toán Riemann tổng quát với cấu hình nghiệm phức hợp đầy đủ gồm 4 trường hợp kết hợp:
    1. Hai sóng giãn (2-rarefaction).
    2. Hai sóng sốc (2-shock).
    3. Sóng 1-giãn kết hợp sóng 2-sốc.
    4. Sóng 2-sốc kết hợp sóng 1-giãn.
  • Hướng nghiên cứu tiếp:
    • Mở rộng thuật toán để giải bài toán Riemann tổng quát cho hệ Euler đẳng entropy với các cấu hình đa sóng phức tạp.
    • Nghiên cứu các kỹ thuật nâng cao bậc chính xác của lược đồ sai phân nhằm giảm thiểu thời gian tính toán khi tinh chỉnh lưới không gian $N$ lớn.

Giá trị tham khảo

Khóa luận là tài liệu tham khảo chuyên môn phù hợp cho sinh viên, học viên cao học và nghiên cứu viên thuộc các chuyên ngành Toán học, Toán ứng dụng, Toán - Tin học, Cơ học tính toán và Kỹ thuật thủy khí động lực học:

  • Phần xây dựng nghiệm giải tích bài toán Riemann (Mục 1.3 - 1.5): Cung cấp phương pháp chi tiết để xác định phương trình sóng tĩnh, sóng sốc, sóng giãn và thuật toán chia đôi phân ly nghiệm.
  • Phần phân tích tính tương thích và Định lý Lax-Wendroff (Mục 2.2 - 2.3): Cung cấp tài liệu tham khảo mẫu mực về kỹ thuật chứng minh sự hội tụ của lược đồ sai phân bảo toàn cho phương trình đạo hàm riêng phi tuyến.
  • Phần mã nguồn mô phỏng Matlab (Phụ lục): Cung cấp thuật toán tính toán bước thời gian tự thích nghi theo điều kiện CFL và quy trình cập nhật thông lượng số cho hệ phương trình hai biến trạng thái.

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

1. Tại sao hệ phương trình Euler trong khóa luận lại được gọi là "đẳng entropy" (isentropic)? Hệ được gọi là đẳng entropy vì tác giả giả thiết đại lượng entropy $S$ đo lường sự rối loạn của hệ là một hằng số đồng nhất ($S = \text{const}$). Khi đó, phương trình bảo toàn năng lượng toàn phần được rút gọn và áp suất $p$ trở thành hàm duy nhất phụ thuộc vào mật độ $\rho$ theo phương trình trạng thái $p = \kappa \rho^\gamma$, biến hệ Euler tổng quát ba phương trình thành hệ hai phương trình gồm bảo toàn khối lượng và bảo toàn động lượng.

2. Ma trận Jacobi của hệ Euler đẳng entropy có những giá trị riêng nào và có ý nghĩa gì? Khi chọn véctơ trạng thái $u = (\rho, v)^T$, ma trận Jacobi của hệ có hai giá trị riêng thực phân biệt là $\lambda_1 = v - c$ và $\lambda_2 = v + c$, với $c = \sqrt{p'(\rho)} = \sqrt{\kappa \gamma \rho^{\gamma-1}}$ là vận tốc âm địa phương. Vì $\lambda_1 \neq \lambda_2$ khi $c > 0$, hệ phương trình luôn đạt tính chất hyperbolic ngặt.

3. Tại sao không thể áp dụng trực tiếp các lược đồ sai phân dạng phi bảo toàn cho bài toán này? Như tác giả đã chứng minh qua phương trình Burgers ở Chương 2, các lược đồ sai phân dạng phi bảo toàn (tựa tuyến tính) dù thỏa mãn khai triển Taylor trơn nhưng khi gặp gián đoạn sẽ hội tụ về một nghiệm số sai lệch, lan truyền sóng với vận tốc không đúng theo điều kiện Rankine-Hugoniot. Do đó, bắt buộc phải sử dụng dạng sai phân bảo toàn để đảm bảo hội tụ đúng về nghiệm yếu vật lý theo Định lý Lax-Wendroff.

4. Trạng thái sóng tĩnh nối với trạng thái ban đầu được xác định như thế nào trong khóa luận? Để tồn tại sóng tĩnh nối với trạng thái $(\rho_0, v_0)$, mật độ $\rho$ phải là nghiệm của phương trình phi tuyến $g(\rho) = \kappa \rho^{\gamma+1} - (\rho_0 v_0^2 + \kappa \rho_0^\gamma)\rho + \rho_0^2 v_0^2 = 0$. Tác giả chứng minh phương trình có nghiệm $\rho_1 = \rho_0$ và nghiệm thứ hai $\rho_2$ nằm trong khoảng phân ly $[\rho', \rho_{\max}]$ (hoặc $[0, \rho']$), sau đó dùng phương pháp chia đôi để tính xấp xỉ $\rho_2$ trước khi đưa vào bài toán Riemann.

5. Tham số CFL ảnh hưởng như thế nào đến độ chính xác và thời gian mô phỏng trong Thử nghiệm số 4? Khi giảm số CFL xuống thấp (ví dụ $\text{CFL} = 0.1$), bước thời gian $\Delta t$ bị chia rất nhỏ dẫn đến số bước lặp tăng vọt, thời gian tính toán tăng gấp 10 lần. Đồng thời, việc lặp quá nhiều bước làm tích lũy sai số khuếch tán số, khiến sai số tăng từ $10^{-2}$ lên $10^{-1}$. Do đó, việc duy trì CFL tiệm cận 1 (nhưng vẫn thỏa mãn $\text{CFL} < 1$) mang lại độ chính xác cao nhất và thời gian thực thi ngắn nhất.

Kết luận

Khóa luận của tác giả Nguyễn Thị Ngân Trúc đã thiết lập thành công mô hình toán học và xây dựng lược đồ sai phân hữu hạn dạng bảo toàn Lax-Friedrichs cho hệ phương trình Euler đẳng entropy một chiều. Công trình đã hoàn thành việc chứng minh tính tương thích, phân tích điều kiện ổn định CFL và kiểm chứng định lượng khả năng bắt các cấu hình sóng tĩnh, sóng giãn, sóng sốc thông qua các thử nghiệm số trên Matlab. Đây là một nghiên cứu có kết cấu chặt chẽ, đóng vai trò nền tảng quan trọng cho việc tiếp tục mở rộng giải quyết bài toán Riemann tổng quát trong động lực học chất khí.