Luận án tiến sĩ mô phỏng xói mòn ống dẫn - Zhihe Zhao

Dòng chảy hai pha quanh ống dẫn bao gồm pha lỏng (nước) và pha rắn (hạt cát). Tương tác giữa hai pha tạo nên độ phức tạp của quá trình xói mòn. Hạt rắn lơ lửng di chuyển theo nhiều chế độ khác nhau tùy thuộc vận tốc dòng chảy. Mô hình Euler-Euler xử lý cả hai pha như các môi trường liên tục xen kẽ nhau. Phương trình động lượng và liên tục được giải đồng thời cho mỗi pha. Lực tương tác giữa các hạt và giữa hạt với chất lỏng được tính toán chi tiết. Áp suất ống dẫn và trường vận tốc ảnh hưởng trực tiếp đến phân bố hạt rắn.

Mô hình rối k-ε được cải tiến để mô tả chính xác dòng chảy rối quanh ống. Năng lượng rối (k) và tốc độ tiêu tán (ε) được giải cho pha lỏng. Mô hình này phù hợp với dòng chảy có số Reynolds cao. Các hệ số trong phương trình được điều chỉnh cho môi trường hai pha. Độ nhớt rối được tính từ k và ε để đóng hệ phương trình. Mô hình cho phép dự đoán chính xác trường vận tốc và ứng suất Reynolds.

Ranh giới giữa đáy cát và nước được xác định bằng ngưỡng phân số thể tích. Giá trị ngưỡng được chọn dựa trên nồng độ hạt tối đa có thể. Phương pháp này cho phép theo dõi sự biến đổi địa hình theo thời gian. Giao diện di động được cập nhật liên tục trong quá trình mô phỏng. Kỹ thuật này hiệu quả hơn so với phương pháp Lagrangian theo dõi từng hạt.

Tải đáy xảy ra khi hạt cát lăn và trượt dọc theo đáy. Vận tốc dòng chảy gần đáy đủ lớn để khởi động hạt nhưng không nhấc chúng lên cao. Tải lơ lửng xuất hiện khi hạt rắn lơ lửng được nâng lên khỏi đáy bởi dòng rối. Nồng độ hạt trong lớp lơ lửng giảm dần theo độ cao. Tải dạng lớp là chế độ trung gian giữa tải đáy và tải lơ lửng. Hạt di chuyển theo các lớp mỏng sát đáy với nồng độ cao. Ba chế độ này đóng góp khác nhau vào quá trình xói mòn.

Hố xói bắt đầu hình thành ngay dưới ống do tăng vận tốc cục bộ. Giai đoạn đầu xói mòn diễn ra nhanh với tốc độ cao nhất. Sau đó quá trình chậm dần khi hố xói đạt độ sâu nhất định. Hình dạng hố xói thay đổi từ chữ V sang chữ U theo thời gian. Độ sâu xói cực đại phụ thuộc vào kích thước ống và vận tốc dòng chảy. Kết quả mô phỏng phù hợp tốt với dữ liệu thí nghiệm về cả hình dạng và độ sâu.

Các công thức dự đoán xói mòn cho vật thể trụ được kiểm chứng. Một số công thức cho kết quả phù hợp tốt với mô phỏng CFD. Mô hình chôn vùi mìn vận hành cho thấy sự khác biệt đáng kể. Nguyên nhân có thể do các giả định đơn giản hóa trong mô hình thực nghiệm. Mô phỏng hai pha cung cấp chi tiết hơn về cơ chế vật lý.

Bước một thiết lập xói mòn đáy sống quanh ống cố định. Dòng chảy mang trầm tích liên tục từ thượng nguồn tạo điều kiện đáy sống. Hố xói phát triển đến trạng thái cân bằng động. Bước hai bắt đầu khi hố xói đạt hình dạng ổn định. Ống được hạ xuống với vận tốc kiểm soát theo phương thẳng đứng. Quá trình võng tiếp tục cho đến khi ống chạm đáy hố xói. Phương pháp này mô phỏng chính xác điều kiện thực tế.

Ba vận tốc võng khác nhau được nghiên cứu: chậm, trung bình và nhanh. Vận tốc võng ảnh hưởng đến phân bố áp suất ống dẫn. Võng chậm cho phép dòng chảy thích ứng dần với vị trí ống mới. Võng nhanh tạo ra nhiễu loạn mạnh hơn trong trường dòng. Hình dạng hố xói cuối cùng phụ thuộc vào tốc độ võng. Nồng độ hạt lơ lửng cũng thay đổi theo vận tốc võng.

Kết quả được so sánh với nghiên cứu sử dụng mô phỏng xoáy lớn (LES). Cả hai phương pháp cho hình dạng xói tương tự nhau. Mô hình hai pha có ưu điểm về chi phí tính toán thấp hơn. LES cung cấp chi tiết hơn về cấu trúc rối nhưng tốn kém hơn. Sự phù hợp giữa hai phương pháp xác nhận độ tin cậy của mô hình hai pha.

Phương trình liên tục cho pha lỏng mô tả bảo toàn khối lượng nước. Phân số thể tích pha lỏng nhân với mật độ và vận tốc xuất hiện trong phương trình. Đạo hàm theo thời gian và không gian được tính cho trường dòng. Phương trình cho pha rắn tương tự nhưng với các tham số riêng. Không có số hạng nguồn khối lượng vì không có pha nào sinh ra hoặc mất đi. Tổng phân số thể tích hai pha luôn bằng một.

Phương trình động lượng cho pha lỏng bao gồm các số hạng đối lưu và khuếch tán. Gradient áp suất và ứng suất nhớt đóng vai trò quan trọng. Lực tương tác với pha rắn xuất hiện như số hạng nguồn. Phương trình cho pha rắn có thêm ứng suất va chạm giữa hạt. Áp suất pha rắn được tính từ lý thuyết động học hạt. Lực kéo giữa hai pha phụ thuộc vào vận tốc tương đối.

Lực kéo là thành phần chính của tương tác pha lỏng-pha rắn. Hệ số kéo phụ thuộc vào số Reynolds của hạt. Lực nâng Magnus và Saffman cũng được xem xét cho hạt nhỏ. Va chạm giữa các hạt được mô hình hóa qua hệ số khôi phục. Tương tác hạt-hạt tăng lên khi nồng độ hạt cao. Vật liệu ống không ảnh hưởng trực tiếp đến tương tác pha.

Lưới được làm mịn trong vùng gần ống và đáy cát. Tỷ lệ kích thước ô lưới được kiểm soát để tránh sai số số. Lớp biên được xây dựng với nhiều lớp ô mỏng sát tường. Giá trị y+ được duy trì trong khoảng phù hợp với mô hình rối. Lưới phi cấu trúc linh hoạt thích ứng với biên di động. Số lượng ô lưới được cân bằng giữa độ chính xác và chi phí tính toán.

Biên vào được thiết lập với phân bố vận tốc đều hoặc theo luật logarit. Nồng độ hạt tại biên vào được chỉ định cho trường hợp đáy sống. Biên ra áp dụng điều kiện áp suất hoặc đối lưu tự do. Tường ống và đáy được xử lý như biên không trượt. Hàm tường được sử dụng để kết nối vùng nhớt với vùng rối. Tham số mô hình như hệ số khuếch tán được hiệu chỉnh.

Sơ đồ upwind bậc hai được chọn cho độ chính xác cao. Thuật toán SIMPLE liên kết áp suất và vận tốc hiệu quả. Bước thời gian được điều chỉnh theo điều kiện Courant. Phần dư của phương trình được theo dõi đến khi đạt ngưỡng hội tụ. Cân bằng khối lượng toàn cục được kiểm tra định kỳ. Kết quả được lưu tại các bước thời gian để phân tích quá trình phát triển.

Kết quả mô phỏng hỗ trợ xác định độ cao đặt ống an toàn. Dự đoán độ sâu xói cực đại giúp tránh hư hỏng ống. Lựa chọn vật liệu ống và lớp bảo vệ dựa trên phân bố ứng suất. Vận tốc dòng chảy tới hạn được xác định cho từng cấu hình. Thiết kế có thể tối ưu hóa chi phí xây dựng và vận hành.

Mô hình có thể mở rộng cho đất kết dính với các tham số khác. Trường hợp nhiều ống song song cần lưới tính toán phức tạp hơn. Ống linh hoạt yêu cầu tích hợp mô hình cấu trúc. Sóng và dòng chảy kết hợp là điều kiện thực tế cần nghiên cứu. Địa hình đáy không đồng nhất tăng độ phức tạp mô phỏng.

Mô hình rối cho dòng hai pha cần cải tiến thêm. Tương tác hạt-hạt ở nồng độ cao cần nghiên cứu sâu hơn. Xác thực với dữ liệu thực địa nâng cao độ tin cậy. Tối ưu hóa thuật toán giảm chi phí tính toán. Tích hợp học máy có thể tăng tốc dự đoán.