Bộ thuật ngữ mô phỏng phương pháp hạt (Particle Method)

Trong lĩnh vực mô phỏng động lực học chất lưu (CFD), Phương pháp Hạt (Particle Method) – đại diện tiêu biểu như MPS (Moving Particle Semi-implicit) hay SPH (Smoothed Particle Hydrodynamics) – đang ngày càng khẳng định thế mạnh vượt trội so với các phương pháp dựa trên lưới (Mesh-based) truyền thống khi xử lý các bài toán bề mặt tự do, biến dạng lớn, văng bắn chất lỏng hay tương tác đa vật lý.

Tuy nhiên, việc nắm vững hệ thống thuật ngữ chuyên ngành trong các phần mềm và tài liệu nghiên cứu về phương pháp hạt luôn là rào cản đối với nhiều kỹ sư. Bài viết này tổng hợp và hệ thống hóa toàn bộ thuật ngữ cốt lõi trong Particle Method theo các nhóm chuyên đề logic, giúp bạn tra cứu dễ dàng và áp dụng hiệu quả vào công việc.

1. Cơ sở lý thuyết & phương trình chi phối (Governing Equations)

Để mô phỏng chính xác hành vi của chất lưu, phương pháp hạt Particle Method dựa trên các định luật vật lý cơ bản và hệ quy chiếu động lực học.

• Governing equation (Phương trình chi phối / Phương trình chủ đạo): Hệ phương trình biểu diễn quy luật của các hiện tượng vật lý. Trong phương pháp hạt, hai phương trình chi phối chính là Phương trình liên tục và Phương trình Navier-Stokes.
• Conservation of mass (Bảo toàn khối lượng): Định luật bảo toàn khối lượng, được thể hiện toán học qua Phương trình liên tục (Equation of continuity).
• Equation of continuity (Phương trình liên tục): Phương trình mô tả mối quan hệ nhiệt động học giữa áp suất, khối lượng riêng và nhiệt độ. Trong MPS/SPH nén yếu thường sử dụng quan hệ giữa áp suất và khối lượng riêng.
• The law of conservation of momentum (Định luật bảo toàn động lượng): Định luật khẳng định tổng động lượng của hệ không đổi nếu không có ngoại lực tác dụng. Hệ phương trình Navier-Stokes chính là dạng biểu diễn của định luật này cho chất lưu.
• Navier-Stokes equations (Các phương trình Navier-Stokes): Phương trình chuyển động cơ bản cho chất lưu thực, mô tả mối quan hệ giữa gia tốc của chất lưu và các lực tác dụng lên nó (áp suất, độ nhớt, trọng lực). Navier-Stokes được xây dựng từ định luật bảo toàn động lượng kết hợp với mô hình ứng suất của chất lưu Newton.
• Equation of motion (Phương trình chuyển động): Phương trình mô tả sự chuyển động của vật chất dưới tác dụng của lực.
• Equation of state (Phương trình trạng thái): Phương trình thiết lập mối quan hệ giữa áp suất (P) và khối lượng riêng (p) của chất lưu.
• Lagrangian description (Mô tả Lagrangian): Hệ quy chiếu mô tả các đại lượng vật lý bằng cách di chuyển đi theo chính hạt/phần tử vật chất đó. Đây là nền tảng cốt lõi của Phương pháp Hạt.
• Euler description (Mô tả Euler): Hệ quy chiếu mô tả các đại lượng vật lý tại các vị trí cố định trong không gian (quan sát dòng chảy đi qua điểm cố định).
• Substantial derivative / Lagrangian derivative (Đạo hàm toàn phần / Đạo hàm Lagrangian): Đạo hàm theo thời gian của một đại lượng vật lý được đánh giá bằng cách di chuyển cùng với phần tử vật chất.

Các số hạng trong phương trình chuyển động:

• Advection term (Số hạng đối lưu): Xuất hiện khi chuyển đổi đạo hàm toàn phần sang hệ quy chiếu Euler, thể hiện sự biến đổi của đại lượng theo thời gian do sự dịch chuyển vị trí của dòng chảy.
• Pressure term (Số hạng áp suất): Thể hiện ảnh hưởng của độ dốc áp suất (pressure gradient) lên chuyển động của chất lưu.
• Viscous term (Số hạng độ nhớt): Thể hiện ảnh hưởng của độ nhớt, làm khuếch tán động lượng trong dòng chảy.
• Gravity term (Số hạng trọng lực): Thành phần thể hiện tác động của trường trọng lực lên hạt chất lưu.
• External force (Ngoại lực): Các lực tác động từ bên ngoài vào khối chất lưu (ví dụ: trọng lực, từ trường).

2. Toán tử vi phân & mô hình tương Tác Hạt (Differential Operators & Discretization)

Khác với phương pháp lưới sử dụng sai phân hay thể tích hữu hạn, phương pháp hạt xấp xỉ các toán tử vi phân dựa trên sự tương tác giữa các hạt lân cận.

3. Điều kiện biên & xử lý hạt tường (Boundary conditions & wall treatment)

Xử lý ranh giới tường và bề mặt tự do là yếu tố quyết định tính chính xác và độ ổn định của mô phỏng hạt.

• Boundary condition (Điều kiện biên): Các ràng buộc vật lý đặt tại ranh giới của miền tính toán.
• Dirichlet boundary condition (Điều kiện biên Dirichlet): Điều kiện gán trực tiếp giá trị của hàm số tại biên (ví dụ: gán áp suất P = 0 tại bề mặt tự do).
• Neumann boundary condition (Điều kiện biên Neumann): Điều kiện gán giá trị đạo hàm của hàm số tại biên (ví dụ: đạo hàm áp suất theo phương pháp tuyến bằng 0 tại tường cứng).
• Free surface (Bề mặt tự do): Bề mặt phân cách biến dạng tự do theo chuyển động chất lưu (ví dụ: mặt chất lỏng tiếp xúc không khí).
• Free surface particle (Hạt bề mặt tự do): Các hạt nằm ở lớp ngoài cùng của bề mặt chất lưu.
• Wall particle (Hạt tường): Hạt cấu tạo nên thành cứng. Hạt tường có tham số áp suất và cố định vị trí trong không gian.
• Dummy wall particle (Hạt tường ảo / Hạt tường đệm): Hạt đặt phía sau hạt tường chính để hỗ trợ tính toán mật độ số hạt nhưng không tham gia giải phương trình áp suất (độ dốc áp suất bằng 0).
• Polygon wall boundary (Biên tường đa giác): Kỹ thuật biểu diễn tường cứng bằng các mặt đa giác (CAD mesh) thay vì xếp hạt tường, giúp giảm số lượng hạt và mô tả bề mặt trơn mượt hơn.
• Weight function for wall (Hàm trọng số cho tường đa giác): Hàm trọng số đặc thù dùng để tính tương tác giữa hạt chất lưu và tường đa giác.
• Free-slip boundary (Biên trượt tự do): Điều kiện biên coi ma sát bằng 0 tại tường, chất lưu trượt tự do trên bề mặt.
• No-slip boundary (Biên không trượt): Điều kiện biên dính chặt, vận tốc chất lưu tại bề mặt tường bằng vận tốc của tường.
• Ghost particle (Hạt ma / Hạt ảo): Các hạt văng ra ngoài miền mô phỏng và không còn tham gia vào quá trình tính toán.

4. Phương pháp số, thuật toán & tối ưu tính toán

Để giải hệ phương trình vi phân phi tuyến theo thời gian, các thuật toán số và kỹ thuật lập trình song song được huy động tối đa.

Thuật toán thời gian & Phương pháp giải số:

• Explicit method (Phương pháp hiện): Phương pháp tìm trạng thái tương lai chỉ dựa trên các giá trị đã biết ở hiện tại. Đơn giản nhưng bị giới hạn nghiêm ngặt bởi bước thời gian Δt
• Implicit method (Phương pháp ẩn): Phương pháp tìm trạng thái tương lai bằng cách giải hệ phương trình đồng thời chứa các ẩn số ở tương lai. Cho phép dùng bước thời gian lớn hơn.
• Semi-implicit method (Phương pháp bán ẩn): Kết hợp cả hai quy trình (ví dụ: phương pháp MPS giải áp suất bằng phương pháp ẩn, các thành phần vận tốc khác giải bằng phương pháp hiện).
• Fractional step method (Phương pháp bước phân đoạn): Lược đồ tách quá trình tính toán thời gian thành các bước nhỏ nối tiếp nhau.
• Time marching method (Phương pháp tiến theo thời gian): Quy trình tích phân hệ phương trình theo các bước thời gian hữu hạn Δt
• Time increment (Δt – Bước thời gian): Khoảng chênh lệch thời gian giữa hai bước tính toán liên tiếp.
• Courant number (Chỉ số Courant – c = uΔt/Δx): Đại lượng không chiều khống chế bước thời gian Δt để đảm bảo độ ổn định trong phương pháp hiện.
• Diffusion number (Chỉ số khuếch tán – D = νΔt/(Δx)²): Đại lượng không chiều giới hạn bước thời gian dựa trên độ nhớt của chất lưu.
• SMAC method (Phương pháp SMAC): Viết tắt của Simplified Marker and Cell, một thuật toán tích phân thời gian phổ biến trong động lực học chất lưu.

Phương trình Poisson & Bộ giải hệ phương trình (Solvers):

• Poisson equation (Phương trình Poisson): Phương trình vi phân đạo hàm riêng dạng ∇²P = b, dùng để giải trường áp suất trong chất lưu không nén được.
• Source term (Số hạng nguồn): Vế phải b của phương trình Poisson.
• Conjugate gradient method (Phương pháp độ dốc liên hợp): Thuật toán giải lặp hiệu năng cao cho hệ phương trình đại số tuyến tính đối xứng.
• Gaussian elimination (Khử Gauss): Phương pháp giải trực tiếp hệ phương trình bằng cách biến đổi dòng.
• SOR method (Phương pháp lặp SOR – Successive Over-Relaxation): Phương pháp giải lặp cải tiến giúp tăng tốc độ hội tụ khi giải hệ phương trình.

Tối ưu hóa cấu trúc dữ liệu & Lập trình song song:

• Bucket (Ô phân vùng / Thuật toán Bucket): Chia miền tính toán thành các ô nhỏ để thu hẹp không gian tìm kiếm hạt lân cận, giúp giảm đáng kể chi phí tính toán.
• Neighboring particles (Các hạt lân cận): Danh sách các hạt nằm trong bán kính ảnh hưởng của hạt đang xét.
• Neighboring particle search (Tìm kiếm hạt lân cận): Thuật toán truy tìm các hạt lân cận xung quanh một hạt.
• MPI (Message Passing Interface): Thư viện lập trình song song dành cho hệ thống máy tính có bộ nhớ phân tán (distributed memory).
• OpenMP: Thư viện lập trình song song dành cho máy tính nhiều nhân sử dụng chung bộ nhớ (shared memory).
• Compiler (Trình biên dịch): Phần mềm chuyển mã nguồn (C++, Fortran) thành file thực thi (Visual Studio, gcc, Intel compiler…).

5. Hiện tượng vật lý & trực quan hóa kết quả (Visualization)

Phương pháp hạt đặc biệt xuất sắc trong việc tái hiện các hiện tượng dòng chảy phức tạp và tương tác đa vật lý.

Hiện tượng & Tính chất chất lưu:

• Collapse of a water column / Dam break (Sụp đổ cột nước / Vỡ đập): Bài toán kinh điển mô tả khối nước đứng yên bị đổ sụp do trọng lực, dùng để kiểm tra độ chính xác của thuật toán hạt.
• Sloshing (Hiện tượng sánh chất lưu): Dao động sóng sánh của chất lỏng trong bình chứa bị lắc mạn (xe bồn, tàu chở dầu).
• Surface tension (Sức căng bề mặt): Lực kéo bề mặt sinh ra do lực liên kết giữa các phân tử chất lỏng.
• CSF model (Mô hình CSF – Continuum Surface Force): Mô hình quy lực căng bề mặt về dạng lực thể tích tập trung tại biên.
• Potential model (Mô hình thế năng): Mô hình tính lực căng bề mặt dựa trên lực tương tác thế năng hút/đẩy trực tiếp giữa các hạt.
• Wetness (Độ dính ướt / Tính thấm ướt): Khả năng bám dính của chất lỏng lên bề mặt tường cứng, thể hiện qua góc tiếp xúc (contact angle).
• Highly viscous fluid (Chất lưu độ nhớt cao): Các chất lưu có độ nhớt lớn như keo, mật ong, nhựa đường, dung dịch polymer.
• Fluid friction (Ma sát chất lưu): Lực cản sinh ra do độ nhớt động năng.
• Turbulent flow (Dòng chảy rối): Dòng chảy xuất hiện các xoáy lốc bất quy tắc và làm tăng cường trao đổi động lượng.
• Multi-physics (Đa vật lý): Bài toán kết hợp nhiều hiện tượng vật lý cùng lúc (ví dụ: chất lưu tương tác với kết cấu biến dạng FSI, chất lưu kết hợp truyền nhiệt…).

Đại lượng đo lường & Trực quan hóa 3D:

• Particle number density (Mật độ số hạt): Tổng giá trị hàm trọng số của các hạt lân cận xung quanh một điểm, đại diện cho mật độ chất lưu cục bộ.
• Spatial resolution (Độ phân giải không gian): Kích thước hạt biểu thị cho độ chi tiết của miền tính toán.
• Initial distance between particles (Khoảng cách hạt ban đầu): Khoảng cách giữa các hạt khi khởi tạo lưới ban đầu.
• Fluid density (Khối lượng riêng chất lưu): Khối lượng trên một đơn vị thể tích (kg/m³).
• Mass of particle (Khối lượng của hạt): Khối lượng gán cho từng hạt đơn lẻ.
• Particle’s ID (ID của hạt): Mã số định danh duy nhất của từng hạt.
• Distance function (Hàm khoảng cách): Dữ liệu không gian lưu khoảng cách ngắn nhất từ một điểm đến bề mặt vật thể.
• Speed of sound (Vận tốc âm thanh): Vận tốc truyền sóng âm, tham số quan trọng trong phương pháp hạt nén được mờ (Weakly Compressible SPH).
• the MKS system of units (Hệ đơn vị MKS): Hệ đơn vị chuẩn dựa trên Mét, Kilôgam và Giây.
• Visualization (Trực quan hóa): Hiển thị dữ liệu số dưới dạng hình ảnh, biểu đồ hoặc video 3D.
• Marching cube (Thuật toán Marching Cubes): Thuật toán nối các điểm dữ liệu hạt thành bề mặt lưới đa giác (polygon mesh) trơn mượt.
• Metaball (Kỹ thuật Metaball): Kỹ thuật tạo dựng bề mặt đẳng mật độ giúp hiển thị bề mặt tự do của chất lỏng mượt mà.
• V&V (Verification & Validation – Kiểm chứng & Xác thực): Verification (Kiểm chứng): Đánh giá xem mô hình số giải đúng phương trình toán học hay chưa. Validation (Xác thực): Đánh giá xem kết quả mô phỏng có phản ánh đúng thực nghiệm thực tế hay không.

Hiểu rõ bản chất của các thuật ngữ trong Phương pháp Hạt (Particle Method) không chỉ giúp bạn đọc hiểu tài liệu chuyên ngành nhanh hơn mà còn là chìa khóa để thiết lập bài toán, lựa chọn tham số số học (solver parameters) và chẩn đoán lỗi mô phỏng CFD một cách chính xác. Nếu bạn đang ứng dụng các công cụ mô phỏng hạt như MPS hay SPH trong công việc, hãy lưu lại bảng tra cứu này để tham khảo khi cần thiết!

*SDE TECH biên soạn và xử lý nội dung từ Particleworks.