Mô phỏng xoáy lớn (LES) khác gì so với RANS và DNS?

RANS mô hình hóa toàn bộ các quy mô xoáy rối qua giả thuyết thống kê, DNS giải trực tiếp toàn bộ phổ xoáy không cần mô hình hóa nhưng tốn chi phí cực lớn, còn LES phân giải trực tiếp các xoáy lớn và chỉ mô hình hóa các xoáy nhỏ dưới mắt lưới.

Tensor ứng suất dưới lưới (SGS) trong LES phát sinh từ đâu?

Tensor ứng suất dưới lưới phát sinh từ tích phân chập phi tuyến của các số hạng đối lưu khi áp dụng toán tử lọc không gian lên phương trình Navier-Stokes, đại diện cho tương tác động lượng giữa trường dòng chảy đã lọc và các xoáy chưa phân giải.

Tại sao việc phân giải lớp biên sát thành trong LES lại tốn kém chi phí tính toán?

Gần bề mặt thành rắn, các cấu trúc xoáy rối trở nên cực kỳ nhỏ và phụ thuộc vào độ nhớt phân tử, khiến số lượng điểm lưới cần thiết cho LES giải quyết sát thành tăng theo quy luật lũy thừa 13/7 (xấp xỉ 1.86) của số Reynolds.

Mô phỏng xoáy lớn là gì? Cơ sở toán học và mô hình hóa

Mô phỏng xoáy lớn là phương pháp mô phỏng số trong cơ học chất lưu tính toán nhằm phân tích chính xác các dòng chảy rối thông qua việc phân giải trực tiếp các cấu trúc xoáy quy mô lớn và mô hình hóa các cấu trúc xoáy quy mô nhỏ dưới kích thước lưới. Phương pháp này đóng vai trò cầu nối quan trọng giữa các mô hình trung bình thống kê truyền thống và phương pháp mô phỏng trực tiếp toàn phần, mang lại độ tin cậy vật lý vượt trội trong nghiên cứu khí động học, âm học chất lưu và các quá trình truyền nhiệt - cháy phức tạp. Bài viết này trình bày toàn diện về cơ sở toán học của phép lọc không gian, các mô hình xoáy dưới lưới tiêu biểu, yêu cầu tính toán và phạm vi ứng dụng của phương pháp.

Vị trí của mô phỏng xoáy lớn trong cơ học dòng rối

Dòng chảy rối trong tự nhiên và kỹ thuật được đặc trưng bởi sự tồn tại đồng thời của một phổ xoáy trải rộng trên nhiều bậc quy mô không gian và thời gian khác nhau. Theo lý thuyết dòng rối đồng nhất của Kolmogorov, năng lượng động học được truyền từ các cấu trúc xoáy quy mô lớn nhất (phụ thuộc chặt chẽ vào hình học biên của dòng chảy) xuống các cấu trúc xoáy trung gian và cuối cùng bị tiêu tán thành nhiệt năng bởi lực nhớt phân tử tại thang vi mô Kolmogorov.

Trong cơ học chất lưu tính toán, ba cách tiếp cận chính đối với dòng chảy rối bao gồm:

Mô phỏng phương trình Navier-Stokes trung bình Reynolds (RANS): Thực hiện lấy trung bình thời gian hoặc trung bình tập hợp toàn bộ trường dòng chảy, mô hình hóa toàn bộ dao động rối thông qua các phương trình độ nhớt xoáy bổ sung. Phương pháp này có chi phí tính toán thấp nhưng thường đánh mất các đặc tính tức thời phi dừng quan trọng của dòng chảy.
Mô phỏng số trực tiếp (DNS): Phân giải toàn bộ phổ xoáy rối từ kích thước hình học vĩ mô đến thang đo Kolmogorov nhỏ nhất mà không sử dụng bất kỳ giả định mô hình hóa dòng rối nào. Tuy nhiên, khối lượng tính toán của DNS tăng lũy thừa theo số Reynolds, khiến phương pháp này hầu như chỉ giới hạn trong nghiên cứu cơ bản ở số Reynolds thấp và trung bình.
Mô phỏng xoáy lớn (LES): Nằm ở vị trí trung gian giữa RANS và DNS. Phương pháp này phân giải trực tiếp các xoáy lớn mang phần lớn năng lượng động học và chịu trách nhiệm cho quá trình vận chuyển chất lượng - nhiệt lượng, trong khi các xoáy nhỏ hơn mắt lưới (mang tính đẳng hướng và có tính chất phổ quát cao hơn) được đại diện thông qua các mô hình xoáy dưới lưới (subgrid-scale models).

Cơ sở toán học của phép lọc không gian

Nguyên lý cốt lõi của mô phỏng xoáy lớn dựa trên việc áp dụng một toán tử lọc không gian lên phương trình Navier-Stokes cho chất lưu không nén được. Vận tốc lọc $\bar{u}_i(\boldsymbol{x}, t)$ tại tọa độ không gian $\boldsymbol{x}$ được xác định qua tích phân chập với hàm nhân lọc $G(\boldsymbol{x}, \boldsymbol{\xi}; \Delta)$ :

\bar{u}_i(\boldsymbol{x}, t) = \int_{\Omega} u_i(\boldsymbol{\xi}, t) G(\boldsymbol{x}, \boldsymbol{\xi}; \Delta) \, d\boldsymbol{\xi}

Trong đó $\Delta$ biểu diễn độ rộng của bộ lọc, thường được chọn tỷ lệ thuận với kích thước mắt lưới tính toán. Khi áp dụng toán tử lọc này lên phương trình bảo toàn động lượng Navier-Stokes cho dòng không nén được, ta thu được hệ phương trình bảo toàn vận tốc lọc:

\frac{\partial \bar{u}_i}{\partial t} + \frac{\partial (\bar{u}_i \bar{u}_j)}{\partial x_j} = -\frac{1}{\rho} \frac{\partial \bar{p}}{\partial x_i} + \nu \frac{\partial^2 \bar{u}_i}{\partial x_j \partial x_j} - \frac{\partial \tau_{ij}}{\partial x_j}

Sự xuất hiện của số hạng $\tau_{ij}$ trong phương trình trên đại diện cho tensor ứng suất dưới lưới (subgrid-scale stress tensor), được định nghĩa bởi:

\tau_{ij} = \overline{u_i u_j} - \bar{u}_i \bar{u}_j

Tensor $\tau_{ij}$ biểu diễn tương tác động lượng giữa các cấu trúc xoáy nhỏ chưa được phân giải và trường dòng chảy quy mô lớn đã được lọc. Để khép kín hệ phương trình, các nhà nghiên cứu cần xây dựng các mô hình toán học nhằm biểu diễn $\tau_{ij}$ theo các biến số đã biết của trường vận tốc lọc.

Lịch sử hình thành và công trình đột phá của Deardorff

Ý tưởng nền tảng về việc mô hình hóa các quy mô nhỏ trong mô phỏng khí quyển được Joseph Smagorinsky đề xuất đầu tiên trong các nghiên cứu khí quyển ban đầu. Tuy nhiên, việc áp dụng thành công phương pháp này vào một bài toán dòng chảy giới hạn thành trong cơ học chất lưu tính toán được ghi nhận bởi nhà khoa học James Deardorff.

Trong công trình mang tính bước ngoặt công bố năm 1970 trên tạp chí Journal of Fluid Mechanics, Deardorff đã thực hiện mô phỏng số thành công cho dòng chảy kênh ba chiều ở số Reynolds cao. Công trình của Deardorff năm 1970 đã chứng minh tính khả thi của việc kết hợp việc phân giải trực tiếp các xoáy cấu trúc lớn với mô hình hóa các xoáy dưới lưới, mở ra một kỷ nguyên mới cho việc ứng dụng mô phỏng xoáy lớn trong cơ học chất lưu hiện đại.

Mô hình xoáy dưới lưới động và đóng góp của Lilly

Mô hình đại số kinh điển nhất để khép kín tensor ứng suất dưới lưới là mô hình Smagorinsky, dựa trên giả thuyết độ nhớt xoáy của Boussinesq:

\tau_{ij} - \frac{1}{3} \tau_{kk} \delta_{ij} = -2 \nu_{sgs} \bar{S}_{ij}

Trong đó $\bar{S}_{ij} = \frac{1}{2} \left( \frac{\partial \bar{u}_i}{\partial x_j} + \frac{\partial \bar{u}_j}{\partial x_i} \right)$ là tensor tốc độ biến dạng của trường vận tốc lọc, và độ nhớt xoáy dưới lưới được tính theo công thức:

\nu_{sgs} = (C_s \Delta)^2 |\bar{S}|

với $|\bar{S}| = \sqrt{2 \bar{S}_{ij} \bar{S}_{ij}}$ . Nhược điểm lớn của mô hình Smagorinsky tiêu chuẩn là hệ số $C_s$ là một hằng số cố định, dẫn đến hiện tượng tiêu tán năng lượng quá mức trong các vùng dòng chảy tầng hoặc vùng chuyển tiếp, đồng thời không tự động triệt tiêu tại bề mặt thành rắn.

Để khắc phục hạn chế này, Germano và các cộng sự đã đề xuất mô hình xoáy dưới lưới động. Tiếp nối ý tưởng này, trong bài báo công bố năm 1992 trên tạp chí Physics of Fluids A, Douglas Lilly đã đề xuất phương pháp bình phương tối thiểu để tối ưu hóa việc xác định hệ số mô hình:

C(\boldsymbol{x}, t) = \frac{1}{2} \frac{\langle L_{ij} M_{ij} \rangle}{\langle M_{kl} M_{kl} \rangle}

Trong đó $L_{ij} = \widehat{\bar{u}_i \bar{u}_j} - \hat{\bar{u}}_i \hat{\bar{u}}_j$ là tensor Germano đại diện cho ứng suất rối giải quyết được giữa cấp độ lọc lưới và cấp độ lọc kiểm tra, $M_{ij}$ biểu diễn sự khác biệt giữa hai toán tử biến dạng tương ứng, và dấu ngoặc nhọn biểu thị phép lấy trung bình không gian dọc theo các hướng dòng chảy đồng nhất. Phương pháp của Lilly năm 1992 cho phép xác định hệ số mô hình biến thiên linh hoạt theo không gian và thời gian, triệt tiêu độ nhớt xoáy một cách tự nhiên tại các vùng dòng không có nhiễu loạn rối mà không cần điều chỉnh hệ số thủ công.

Mô hình WALE và kiểm soát ứng suất sát thành

Mặc dù mô hình động Germano-Lilly mang lại độ chính xác cao, việc tính toán lấy trung bình dọc theo các hướng đồng nhất thường gặp khó khăn trong các hình học dòng chảy phức tạp phi đối xứng.

Để giải quyết bài toán này, Nicoud và Ducros trong công trình công bố năm 1999 trên tạp chí Flow, Turbulence and Combustion đã giới thiệu mô hình WALE (Wall-Adapting Local Eddy-viscosity). Mô hình WALE được xây dựng dựa trên bình phương của tensor gradient vận tốc, cho phép phân biệt rõ ràng giữa chuyển động quay của dòng xoáy và chuyển động biến dạng trượt thuần túy. Biểu thức độ nhớt xoáy dưới lưới của mô hình WALE được xác định bởi:

\nu_{sgs} = (C_w \Delta)^2 \frac{(S_{ij}^d S_{ij}^d)^{3/2}}{(\bar{S}_{ij} \bar{S}_{ij})^{5/2} + (S_{ij}^d S_{ij}^d)^{5/4}}

Trong đó $S_{ij}^d = \frac{1}{2} (\bar{g}_{ik}\bar{g}_{kj} + \bar{g}_{jk}\bar{g}_{ki}) - \frac{1}{3} \delta_{ij} \bar{g}_{kk}^2$ là tensor đối xứng không vết của bình phương gradient vận tốc và $\bar{g}_{ij} = \frac{\partial \bar{u}_i}{\partial x_j}$ . Đặc tính nổi bật của mô hình Nicoud và Ducros năm 1999 là độ nhớt xoáy dưới lưới tự động suy giảm về không theo bậc lũy thừa ba khi tiến sát bề mặt thành rắn, tái hiện chính xác hành vi vật lý của lớp đệm dính mà không cần áp dụng các hàm suy giảm khoảng cách nhân tạo.

Yêu cầu tính toán và thách thức phân giải lưới

Một trong những yếu tố quyết định sự thành công và chi phí tính toán của mô phỏng xoáy lớn là việc thiết kế lưới tính toán tại các vùng lớp biên sát thành.

Trong bài tổng quan chuyên sâu công bố năm 1999 trên tạp chí Progress in Aerospace Sciences, Ugo Piomelli đã phân tích toàn diện các thành tựu và thách thức tính toán của mô phỏng xoáy lớn. Piomelli năm 1999 chỉ ra rằng thách thức lớn nhất của LES nằm ở lớp biên sát thành, nơi quy mô của các xoáy rối trở nên vô cùng nhỏ bé và phụ thuộc trực tiếp vào độ dày nhớt của chất lưu.

Vấn đề này được định lượng hóa rõ nét trong công trình nghiên cứu của Choi và Moin công bố năm 2012 trên tạp chí Physics of Fluids. Bằng các phân tích giải tích và thống kê số liệu dòng chảy lớp biên, Choi và Moin năm 2012 đã chứng minh rằng đối với phương pháp LES giải quyết đầy đủ lớp biên sát thành (wall-resolved LES), số lượng điểm lưới cần thiết tăng theo quy luật tỷ lệ chính xác là $Re_{L_x}^{13/7}$ (xấp xỉ $Re_{L_x}^{1.86}$ ). Trong khi đó, nếu áp dụng kỹ thuật mô hình hóa lớp biên sát thành (wall-modeled LES), số điểm lưới chỉ tăng xấp xỉ tuyến tính theo số Reynolds. Kết quả này lý giải vì sao các mô phỏng công nghiệp ở số Reynolds lớn thường kết hợp LES với các mô hình thành hoặc chuyển sang các phương pháp lai ghép.

So sánh toàn diện các phương pháp mô phỏng dòng rối

Bảng dưới đây tóm tắt các đặc tính kỹ thuật chính để phân biệt ba trường phái mô phỏng số dòng chảy rối trong cơ học chất lưu:

Tiêu chí đánh giá	RANS	LES	DNS
Bản chất toán học	Lấy trung bình thống kê toàn bộ	Lọc không gian theo quy mô xoáy	Giải trực tiếp toàn bộ phổ không lọc
Quy mô phân giải	Chỉ phân giải chuyển động trung bình	Phân giải trực tiếp các xoáy lớn	Phân giải toàn bộ đến thang Kolmogorov
Mô hình hóa dòng rối	Mô hình hóa toàn bộ dao động rối	Chỉ mô hình hóa xoáy dưới lưới	Không sử dụng bất kỳ mô hình rối nào
Tính chất trường dòng	Trường dòng trung bình dừng/phi dừng	Trường dòng tức thời 3 chiều biến thiên	Trường dòng vật lý tức thời không chứa sai số mô hình hóa dòng rối (chỉ tồn tại sai số rời rạc hóa số trị)
Chi phí tính toán	Thấp, phù hợp thiết kế công nghiệp	Trung bình - cao, đòi hỏi máy tính hiệu năng cao	Cực lớn, tăng lũy thừa theo số Reynolds

Phạm vi ứng dụng và định hướng phát triển

Nhờ khả năng mô tả chính xác sự phát triển và tương tác của các cấu trúc xoáy kết hợp trong không gian và thời gian, mô phỏng xoáy lớn được ứng dụng sâu rộng trong các lĩnh vực kỹ thuật then chốt:

Khí động lực học hàng không và phương tiện giao thông: Phân tích hiện tượng dòng chảy tách rời khỏi cánh máy bay ở góc tấn lớn, hiện tượng hình thành vệt xoáy đuôi tàu cao tốc và xe hơi nhằm tối ưu hóa lực cản và tăng tính ổn định khí động.
Âm học chất lưu (aeroacoustics): Dự đoán chính xác tiếng ồn sinh ra do luồng phản lực máy bay, cánh quạt tuabin gió hoặc tiếng ồn gió quanh gương chiếu hậu xe hơi, nơi các nguồn âm phụ thuộc trực tiếp vào sự dao động áp suất tức thời của xoáy lớn.
Quá trình hòa trộn và buồng đốt phản ứng: Mô phỏng sự tương tác giữa trường xoáy rối và ngọn lửa trong động cơ tuabin khí và động cơ tên lửa, nơi hiện tượng tắt lửa cục bộ hoặc bất ổn định nhiệt âm phụ thuộc vào cấu trúc xoáy tức thời.
Khí tượng đô thị và phát tán ô nhiễm: Dự báo sự phân tán của các chất gây ô nhiễm không khí xung quanh các tòa nhà cao tầng và địa hình phức tạp dưới tác động của gió giật khí quyển.