Mô Phỏng Số Xói Mòn Quanh Đường Ống Cố Định và Võng Sử Dụng Mô Hình Hai Pha

Luận án tiến sĩ: Mô phỏng số xói mòn quanh đường ống cố định và võng bằng mô hình hai pha. Nghiên cứu kỹ thuật và ứng dụng mô hình số trong xây dựng.

• Tác giả: Zhihe Zhao
• Người hướng dẫn: Professor Harindra J. Fernando
• Trường học: Arizona State University
• Đề tài: Numerical Simulation Of Scour Around Fixed And Sagging Pipelines Using A Two-Phase Model
• Loại tài liệu: Dissertation
• Năm xuất bản: 2006

Mục lục chi tiết

• ACKNOWLEDGEMENT
• ABSTRACT
• INFORMATION TO USERS
• LIST OF FIGURES
• 1. CHƯƠNG 1: SIMULATION OF SCOUR AROUND A FIXED PIPELINE
• 1.1. Overview of Scour Models
• 1.2. MATHEMATICAL DESCRIPTION OF THE TWO-PHASE MODEL
• 1.3. Turbulence Closure for Fluid Phase
• 1.4. Transport Equation for Granular Temperature
• 1.5. NUMERICAL SIMULATION AND VALIDATION: FIXED PIPELINES
• 1.5.1. Mao’s Experimental Set-up
• 1.5.2. Simulation with Fluent
• 1.6. RESULTS AND DISCUSSION: SIMULATION WITH A FIXED PIPELINE
• 1.6.1. Clear-water Scour Simulation
• 1.6.2. Sediment Transport Modes
• 1.6.2.1. Bed-load, suspended-load and laminated-load
• 1.6.2.2. Calculation of bed-load and suspended-load
• 1.6.2.3. Comparison with NBURY and DRAMBUIE Models
• 1.6.3. Conclusions: Scour Under Fixed Pipelines
• 2. SIMULATION OF SCOUR BELOW A SAGGING PIPELINE
• 2.1. Overview of Sagging Pipeline Studies
• 2.2. TWO-PHASE MODEL, SIMULATION AND VALIDATION: SAGGING PIPELINES
• 2.2.1. Fredsge et al.’s Experimental Set-up
• 2.3. RESULTS AND DISCUSSIONS: SAGGING PIPELINES
• 2.3.1. Live-bed scour around a fixed pipeline
• 2.3.2. Scour around a sagging pipeline
• 2.3.3. Comparisons with Cheng and Li’s Simulation
• 2.3.4. Vortex shedding in scouring process
• 2.3.5. Conclusions: Scour under Sagging Pipelines
• REFERENCES

Tóm tắt

I. Tổng Quan Mô Phỏng Số Xói Mòn Đường Ống Tầm Quan Trọng

Xói mòn liên tục quanh đường ống dưới tác động của sóng và dòng chảy ảnh hưởng lớn đến sự ổn định cấu trúc của đường ống và môi trường xung quanh. Hiểu rõ quá trình xói mòn và khả năng dự đoán xói mòn là rất quan trọng trong thiết kế đường ống ngoài khơi. Xói mòn quanh các vật thể đặt trên đáy cát rất phức tạp, thường liên quan đến dòng chảy rối hai pha và nhiều phương thức vận chuyển trầm tích. Giao diện giữa nước và đáy cát cũng phức tạp. Dòng chảy thay đổi hình dạng đáy, từ đó ảnh hưởng đến dòng chảy. Thêm vào sự phức tạp này là cơ học tương tác dòng chảy-hạt và hạt-hạt. Các phương trình Navier-Stokes và các lược đồ đóng gói độ rối thích hợp cần được sửa đổi để tính đến các hiện tượng phức tạp này. Mục tiêu của nghiên cứu này là mô phỏng xói mòn quanh các vật thể hình trụ dài bằng mô hình hai pha Eulerian với hy vọng hiểu rõ xói mòn quanh và sự chôn vùi của các loại mìn thường có hình dạng tương tự.

1.1. Mục Tiêu Nghiên Cứu Mô Phỏng Xói Mòn Đường Ống

Nghiên cứu này là một phần của nỗ lực của chương trình Dự đoán Chôn vùi Mìn (MBP) do Văn phòng Nghiên cứu Hải quân Hoa Kỳ tài trợ, với mục tiêu dài hạn là phát triển các mô hình chôn vùi mìn kết hợp các quá trình môi trường động, các đặc tính vật liệu đáy biển và các loại mìn khác nhau. Các mô hình dự đoán chôn vùi mìn hiện có hoạt động kém, do đó cần thiết phải điều tra hiệu quả của các công thức xói mòn được sử dụng trong các mô hình này [41]. Trong hai thập kỷ qua, một số lượng lớn các mô hình số đã được phát triển để dự đoán xói mòn, nhưng nỗ lực hiện tại được xây dựng trên lý thuyết dòng chảy hai pha được phát triển gần đây.

1.2. Tổng Quan Các Mô Hình Mô Phỏng Xói Mòn Hiện Có

Mao [28] đã áp dụng lý thuyết dòng thế đã sửa đổi và phương trình liên tục trầm tích để mô phỏng xói mòn bên dưới các hình trụ dài đặt trên đáy cát chịu tác động của dòng chảy trung bình. Các dự đoán của mô hình được so sánh với các thí nghiệm trong phòng thí nghiệm bổ sung, và mô hình dòng thế có thể mô phỏng phần thượng nguồn của hố xói mòn một cách thỏa đáng. Li & Cheng [23] cũng đã phát triển một mô hình xói mòn dựa trên lý thuyết dòng thế. Thay vì sử dụng công thức vận chuyển trầm tích thực nghiệm, họ tính toán kích thước hố xói mòn cân bằng bằng cách giả định rằng ứng suất cắt đáy ở mọi nơi trên đáy biển bằng hoặc nhỏ hơn ứng suất cắt trường xa khi trạng thái cân bằng đạt được.

II. Thách Thức Độ Chính Xác Tin Cậy Mô Hình Xói Mòn Số

Mặc dù có nhiều tiến bộ, việc xây dựng mô hình dòng chảy hai pha đáng tin cậy cho pha dòng chảy vẫn là một thách thức do sự hiểu biết hạn chế về tương tác hạt-độ rối. Các công thức vận chuyển trầm tích hiện có có thể không nắm bắt đầy đủ các quá trình vật lý phức tạp liên quan đến xói mòn đường ống. Ngoài ra, việc xác thực các mô hình mô phỏng bằng dữ liệu thực nghiệm là rất quan trọng nhưng đòi hỏi nhiều thời gian và chi phí. Độ chính xác của các mô hình mô phỏng phụ thuộc vào độ phân giải lưới, độ chính xác của các thông số đầu vào và các giả định được đưa ra trong mô hình. Việc sử dụng CFD trong phân tích xói mòn là một hướng đi đầy hứa hẹn, nhưng cần phải giải quyết những thách thức liên quan đến độ chính xác và hiệu quả tính toán.

2.1. Vấn Đề Với Các Mô Hình Dòng Thế Trong Mô Phỏng Xói Mòn

Một kỹ thuật điều chỉnh biên dựa trên phương pháp Newton-Raphson đã được sử dụng cho các mô phỏng. Mô hình của họ cũng dự đoán độ sâu xói mòn thượng nguồn gần đúng một cách hợp lý, nhưng không đưa ra dự đoán chính xác cho phần hạ nguồn của hố xói mòn do những hạn chế của lý thuyết dòng thế. Sau đó, các tác giả tương tự [24] đã giải các phương trình dòng chảy bằng cách sử dụng công thức đóng lưới con (SGS) Smagorinsky. Hố xói mòn cân bằng được dự đoán phù hợp với kết quả thực nghiệm. Tuy nhiên, phương pháp này chỉ dựa vào giả định nói trên về ứng suất cắt đáy và nó chỉ tạo ra các cấu hình xói mòn cân bằng, không phải sự phát triển của xói mòn.

2.2. Hạn Chế Của Mô Hình k epsilon và Yêu Cầu Về Làm Mịn

Leeuwestein et al. [22] đã sử dụng mô hình độ rối k—£ kết hợp với phương trình vận chuyển trầm tích để mô phỏng xói mòn xung quanh các đường ống. Trong nghiên cứu này, chỉ có vận chuyển tải đáy được xem xét đầu tiên, và trong trường hợp này, vận chuyển trầm tích chính xảy ra bởi © gợn sóng. Dự đoán này không thể được chứng minh bằng các thí nghiệm trong phòng thí nghiệm. Các trầm tích lơ lửng sau đó đã được đưa vào mô hình, khi đó các gợn sóng biến mất, cho thấy rằng chỉ vận chuyển tải đáy là không đủ để đại diện cho các quá trình vận chuyển trầm tích phức tạp xung quanh các vật thể rắn. Liang et al. [26] đã thực hiện mô phỏng xói mòn xung quanh các đường ống bằng một phương pháp tương tự. Hai mô hình độ rối, mô hình k-¢ tiêu chuẩn và mô hình Smagorinsky SGS, đã được áp dụng. Để tránh sự xuất hiện của các cấu hình xói mòn không đều sắc nét một cách phi thực tế và sự bất ổn số trong quá trình tính toán của họ, một kỹ thuật làm mịn đặc biệt, được gọi là mô hình trượt cát, đã được sử dụng trong mô phỏng.

2.3. Đánh Giá Khả Năng Ứng Dụng Thực Tế Của Các Mô Hình Hiện Tại

Simulations with the SGS model showed intense vortex shedding during scour, but scour profiles obtained with the k-é€ model were more realistic when compared with Mao's 4 measurements. Therefore, Liang et al. recommended the use of a standard &—e£ model for scour predictions. Dupuis & Chopard [10] proposed a Lattice Boltzman Method to simulate scour around pipelines. In this method, fictitious fluid and sediment particles moved on a regular lattice synchronously at discrete time steps, and time-dependent erosion processes involved were simulated. However, only a portion of the equilibrium scour hole could be quantitatively compared well with laboratory measurements. The model needed to be “tuned” to match various other cases, thus limiting its utility as a robust predictive tool. Ali & Karim [1] used CFD software FLUENT to predict the three-dimensional flow and bed shear stress over a rigid bed. By employing experimental data and the one- dimensional sediment continuity equation, they derived the variation of maximum scour depth with time as a function of the dimensionless bed-shear stress and streamwise distance. Field measurements of scour around bridges were also applied to further verify the latter result. Since the numerical simulation was only limited to a rigid bed, the scouring was not simulated explicitly.

III. Phương Pháp Mô Hình Hai Pha Eulerian Cho Xói Mòn Đường Ống

Nghiên cứu sử dụng mô hình hai pha Eulerian để mô phỏng xói mòn quanh đường ống. Mô hình này xem xét cả pha lỏng và pha rắn, cũng như tương tác giữa chúng. Mô hình này giải các phương trình bảo toàn khối lượng và động lượng cho cả hai pha. Tương tác giữa các hạt được tính đến bằng cách sử dụng lý thuyết động học của dòng hạt. Mô hình độ rối k-epsilon sửa đổi được sử dụng để mô tả độ rối của pha lỏng. Giao diện giữa cát và nước được xác định bằng cách sử dụng một phân số thể tích ngưỡng của cát. Sự phát triển của hình thái đáy được nghiên cứu chi tiết. Cách tiếp cận này cho phép dự đoán chính xác hơn về xói mòn đường ống so với các mô hình chỉ xem xét một pha.

3.1. Ưu Điểm Của Mô Hình Hai Pha Eulerian So Với Mô Hình Kinh Nghiệm

Trong những năm gần đây, các mô hình hai pha, xem xét động lực học của các pha hạt và chất lỏng cũng như các tương tác của chúng, đã được sử dụng để tính toán vận chuyển trầm tích trong khuôn khổ các phương trình Navier-Stokes. Các mô hình như vậy dự đoán vận chuyển trầm tích từ các phương trình động lực học cơ bản hơn (mặc dù được mô hình hóa), do đó tránh được việc sử dụng các công thức vận chuyển trầm tích thuần túy thực nghiệm. Các công thức như vậy tràn lan trong tài liệu và đã được chứng minh là phụ thuộc vào trường hợp, do đó hạn chế việc sử dụng chung của chúng để bao gồm một loạt các cấu hình dòng chảy. Các công thức hai pha, mặt khác, được phát triển dựa trên các khái niệm cơ bản hơn, mặc dù tự nhiên một số tham số hóa là cần thiết để đóng gói. Do đó, các mô hình như vậy được kỳ vọng sẽ có khả năng áp dụng chung hơn cho một loạt các vấn đề.

3.2. Ứng Dụng Của Mô Hình Hai Pha Trong Tính Toán Vận Chuyển Trầm Tích

Liên quan đến các mô hình hai pha, Yeganeh et al. [48] đã sử dụng mô hình hai pha kết hợp Euler-Lagrange để mô phỏng vận chuyển tải đáy dưới ứng suất cắt đáy cao. Mặc dù các kết quả thực nghiệm đã chỉ ra sự tồn tại của cấu hình vận tốc loại ba lớp, mô hình chỉ tạo ra cấu hình vận tốc hai lớp. Các tác giả quy sự khác biệt này cho việc bỏ qua các va chạm giữa các hạt trong mô hình. Hsu et al. [18] đã sử dụng mô hình hai pha để mô phỏng vận chuyển trầm tích lơ lửng và chứng minh khả năng của các mô hình như vậy để dự đoán nồng độ trung bình theo thời gian trong một loạt các điều kiện.

IV. Ứng Dụng Thực Tiễn Mô Phỏng Xói Mòn Quanh Đường Ống Cố Định

Nghiên cứu mô phỏng xói mòn quanh đường ống cố định đặt trên đáy cát không kết dính. Kết quả cho thấy sự phát triển của hình thái đáy phù hợp với các phép đo trong phòng thí nghiệm trước đó. Ba phương thức vận chuyển trầm tích (tải đáy, tải lơ lửng và tải lớp) liên quan đến sự phát triển của xói mòn. Một số công thức phát triển xói mòn được đề xuất trước đó cho các vật thể hình trụ phù hợp với các mô phỏng. Tuy nhiên, các dự đoán xói mòn dựa trên một số mô hình chôn mìn hoạt động cho thấy sự khác biệt với các mô phỏng hiện tại.

4.1. So Sánh Kết Quả Mô Phỏng Với Thí Nghiệm Trong Phòng Thí Nghiệm

Để điều tra quá trình xói mòn, các mô hình và phương pháp tính toán đã được so sánh với các phép đo và thí nghiệm trong phòng thí nghiệm trước đó. Mục đích của việc so sánh này là để đánh giá độ chính xác của mô hình và liệu nó có thể tái tạo trung thực các hiện tượng xói mòn trong môi trường thực tế hay không. Kết quả mô phỏng cho thấy sự phát triển của hình thái đáy phù hợp với các phép đo trong phòng thí nghiệm trước đó. Điều này cung cấp sự tin tưởng vào khả năng của mô hình trong việc dự đoán xói mòn trong điều kiện tương tự.

4.2. Các Phương Thức Vận Chuyển Trầm Tích Trong Quá Trình Xói Mòn

Trong quá trình phát triển xói mòn, ba phương thức vận chuyển trầm tích đã được xác định: tải đáy, tải lơ lửng và tải lớp. Tải đáy là lớp trầm tích di chuyển gần đáy. Tải lơ lửng là trầm tích di chuyển lơ lửng trong chất lỏng do tác động của độ rối. Tải lớp là một lớp mỏng trầm tích di chuyển trên bề mặt đáy. Sự hiểu biết về các phương thức vận chuyển trầm tích này là quan trọng để mô hình hóa chính xác quá trình xói mòn.

V. Nghiên Cứu Xói Mòn Quanh Đường Ống Võng Ảnh Hưởng Vận Tốc

Nghiên cứu cũng điều tra xói mòn dưới đường ống võng. Dòng chảy và đường ống phát triển qua hai bước: (1) xói mòn live-bed cục bộ xung quanh đường ống phát triển xung quanh đường ống cố định; và (2) đường ống được hạ xuống hố xói mòn theo cách có kiểm soát cho đến khi nó chạm đáy hố xói mòn. Ba vận tốc võng được mô phỏng, và các cấu hình xói mòn được dự đoán phù hợp với dữ liệu trong phòng thí nghiệm. Các đặc điểm chung của trường dòng chảy, bao gồm độ rối, sự lơ lửng của các hạt và vận chuyển trầm tích, được mô tả chú ý đến sự phụ thuộc của chúng vào độ võng của đường ống.

5.1. Các Giai Đoạn Phát Triển Xói Mòn Quanh Đường Ống Võng

Quá trình xói mòn quanh đường ống võng diễn ra theo hai giai đoạn chính. Đầu tiên, xuất hiện xói mòn live-bed cục bộ xung quanh đường ống khi nó vẫn ở trạng thái cố định. Sau đó, đường ống từ từ hạ xuống hố xói mòn đã hình thành, tiếp tục thay đổi cấu trúc dòng chảy và quá trình vận chuyển trầm tích. Việc phân tích cả hai giai đoạn này là rất quan trọng để hiểu đầy đủ về hiện tượng xói mòn phức tạp này.

5.2. So Sánh Kết Quả Mô Phỏng Với Nghiên Cứu Dựa Trên LES

Các cấu hình xói mòn được mô phỏng cũng phù hợp với một nghiên cứu số dựa trên LES được báo cáo trước đó. Cấu trúc liên kết dòng chảy và sự phụ thuộc của nó vào độ võng của đường ống được khảo sát. Nghiên cứu này cung cấp thông tin chi tiết về vai trò của sự hình thành xoáy trong các cơ chế xói mòn. Phương pháp LES (Large Eddy Simulation) được sử dụng trong nghiên cứu trước đó cung cấp độ phân giải cao hơn của các cấu trúc độ rối so với các mô hình độ rối RANS (Reynolds-Averaged Navier-Stokes) được sử dụng trong các nghiên cứu khác. Sự phù hợp giữa các kết quả mô phỏng từ các phương pháp khác nhau này tăng cường sự tin tưởng vào độ chính xác của các mô hình.

VI. Kết Luận Ưu Điểm và Hướng Phát Triển Mô Phỏng Xói Mòn

Nghiên cứu này trình bày một mô hình số để mô phỏng xói mòn quanh đường ống cố định và võng. Mô hình này dựa trên phương pháp dòng chảy hai pha Eulerian, cho phép tính đến tương tác giữa pha lỏng và pha rắn. Kết quả của mô phỏng phù hợp với dữ liệu thực nghiệm, cho thấy rằng mô hình có thể được sử dụng để dự đoán xói mòn trong các điều kiện khác nhau. Tuy nhiên, cần nghiên cứu thêm để cải thiện độ chính xác của mô hình và để hiểu rõ hơn về các quá trình vật lý liên quan đến xói mòn.

6.1. Hướng Nghiên Cứu Tương Lai Về Mô Phỏng Xói Mòn Đường Ống

Để nâng cao hơn nữa khả năng mô phỏng xói mòn đường ống, các nghiên cứu trong tương lai nên tập trung vào việc phát triển các mô hình độ rối chính xác hơn, đặc biệt là trong vùng gần tường. Việc cải thiện độ chính xác của việc mô hình hóa tương tác hạt-độ rối cũng rất quan trọng. Nghiên cứu thêm là cần thiết để phát triển các điều kiện biên chính xác hơn cho tải trầm tích. Cuối cùng, cần thực hiện nhiều nghiên cứu hơn để xác thực các mô hình bằng dữ liệu thực địa.

6.2. Tầm Quan Trọng Của Nghiên Cứu Tiếp Theo Về Xói Mòn Đường Ống

Nghiên cứu xói mòn là điều cần thiết để đảm bảo an toàn và tuổi thọ của các đường ống ngầm. Việc phát triển các mô hình chính xác và đáng tin cậy cho mô phỏng xói mòn sẽ cho phép các kỹ sư thiết kế các đường ống có khả năng chống xói mòn và thực hiện các biện pháp khắc phục để giảm thiểu các tác động của xói mòn. Nghiên cứu này sẽ có ý nghĩa quan trọng đối với ngành công nghiệp dầu khí ngoài khơi và cho việc bảo vệ cơ sở hạ tầng dưới nước.

Nội dung chính

Tổng quan luận án

Hiện tượng xói lở cục bộ xung quanh các đường ống đặt dưới đáy biển dưới tác động của sóng và dòng chảy là yếu tố trực tiếp chi phối độ ổn định kết cấu công trình biển cũng như quá trình tự chôn lấp của các vật thể hình trụ (như thủy lôi chống hạm). Về mặt vật lý, xói lở đáy cát là bài toán phức tạp gồm dòng chảy rối hai pha (nước - hạt cát), tương tác động lực giữa chất lưu và hạt, va chạm giữa các hạt, cùng nhiều chế độ vận chuyển bùn cát diễn ra đồng thời. Hầu hết các mô hình số trị kinh điển trước đây đều dựa vào các phương trình thủy động lực học một pha kết hợp với các công thức vận chuyển bùn cát thực nghiệm. Các công thức thực nghiệm này phụ thuộc nhiều vào điều kiện thí nghiệm cụ thể, thiếu tính tổng quát và thường bỏ qua tương tác hai chiều giữa trường dòng chảy và sự biến dạng của đáy. Đồng thời, các mô hình dự báo chôn lấp thủy lôi hiện hành (như NBURY và DRAMBUIE) bộc lộ nhiều sai số khi so sánh với thực tế.

Mục tiêu tổng quát của luận án là ứng dụng và hoàn thiện mô hình hai pha Euler-Euler để mô phỏng số trị quá trình xói lở xung quanh các vật thể hình trụ dài đặt trên đáy cát không kết dính. Cụ thể, nghiên cứu giải quyết hai mục tiêu chính được đánh số như sau:

1. Mô phỏng quá trình xói lở nước trong (clear-water scour) xung quanh đường ống cố định nằm ngang sát đáy cát dưới tác dụng của dòng chảy ổn định, làm sáng tỏ các cơ chế vận chuyển bùn cát và đánh giá độ chính xác của các mô hình dự báo vận hành.
2. Mô phỏng quá trình xói lở đáy động (live-bed scour) xung quanh đường ống cố định và quá trình phát triển hố xói khi đường ống bị võng (sagging pipeline) hạ dần xuống đáy hố xói dưới các vận tốc võng khác nhau.

Đối tượng nghiên cứu là trường dòng chảy rối hai pha, sự biến đổi hình học của đáy cát, và chuyển động của đường ống ngầm hình trụ tròn. Phạm vi nghiên cứu về không gian là miền dòng chảy hai chiều (2D) xung quanh ống trụ nằm ngang trong kênh dẫn hẹp và máng thí nghiệm. Phạm vi thời gian nghiên cứu tập trung vào toàn bộ quá trình phát triển xói lở từ trạng thái ban đầu ($t = 0$) qua các giai đoạn chuyển tiếp ($t = 10, 30, 60, 100, 200\text{ phút}$) cho đến khi đạt trạng thái cân bằng động ($t \ge 300\text{ phút}$).

---

Tổng quan tài liệu và vị trí của luận án

Luận án đã phân tích và hệ thống hóa các hướng tiếp cận mô hình hóa xói lở đáy biển và vận chuyển bùn cát trong hai thập kỷ trước thời điểm công bố:

Hướng tiếp cận	Tác giả, công trình tiêu biểu	Đặc điểm và kết quả chính	Hạn chế ghi nhận trong văn bản
Lý thuyết thế năng (Potential Flow)	Mao (1988) [28]; Li & Cheng (1999) [23]	Kết hợp dòng thế với phương trình liên tục bùn cát hoặc giả thiết ứng suất cắt đáy cân bằng với trường xa; dùng thuật toán Newton-Raphson điều chỉnh biên.	Dự báo tốt phần thượng lưu hố xói nhưng thất bại ở phần hạ lưu do không mô phỏng được dòng tách rời và xoáy sau ống.
Mô hình RANS kết hợp công thức thực nghiệm	Leeuwestein et al. (1985) [22]; Brørs (1999) [3]; Liang et al. (2005) [26]	Sử dụng mô hình rối $k-\epsilon$ giải phương trình RANS; ghép nối tải đáy và tải lơ lửng bằng phương sai hữu hạn hoặc phần tử hữu hạn; dùng mô hình trượt cát (sand-slide).	Mô hình của Brørs dừng phát triển xói sau 100 phút (thực tế thí nghiệm kéo dài sau 300 phút); bất ổn định số học do công thức tải đáy phi tuyến; tải đáy đơn lẻ gây gợn sóng nhân tạo.
Mô hình xoáy lớn (LES / SGS)	Li & Cheng (2000) [24]; Liang et al. (2005) [26]; Cheng & Li (2003)	Ứng dụng mô hình độ nhớt xoáy dưới lưới Smagorinsky; ghi nhận hiện tượng phát tán xoáy mạnh (vortex shedding) ở vùng xói lở.	Chỉ tính toán được hình dạng cân bằng tĩnh hoặc đòi hỏi chi phí tính toán rất lớn; biên dạng xói từ $k-\epsilon$ đôi khi thực tế hơn khi so với số liệu đo đạc của Mao.
Phương pháp mạng Boltzmann (LBM)	Dupuis & Chopard (2002) [10]	Hạt chất lưu và hạt cát giả định chuyển động đồng bộ trên mạng lưới rời rạc; mô phỏng xói phụ thuộc thời gian.	Chỉ so khớp định lượng được một phần hố xói cân bằng; mô hình phải hiệu chỉnh (tune) tham số cho từng trường hợp, tính tổng quát hạn chế.
Mô hình dòng hai pha cho bùn cát	Yeganeh et al. (2000) [48]; Hsu et al. (2004) [18]; Greimann (2004) [16]; Wanker et al. (2001) [43]	Tiếp cận Euler-Lagrange hoặc Euler-Euler; mô phỏng tải đáy dưới ứng suất cao, tải lơ lửng và sự dịch chuyển của đụn cát dựa trên tương tác hạt - chất lưu.	Bỏ qua va chạm giữa các hạt (Yeganeh et al.); chỉ áp dụng phương trình hai pha cho pha rắn mà không giải trọn vẹn cho pha lỏng do thiếu tham số hóa tương tác hạt - rối (Greimann).

Khoảng trống nghiên cứu mà luận án lựa chọn giải quyết là: xây dựng một khuôn khổ mô phỏng số trị dựa trên lý thuyết dòng hai pha Euler-Euler hoàn chỉnh, tích hợp tương tác hai chiều giữa chất lưu và hạt cát, va chạm phi đàn hồi giữa các hạt và ma sát hạt ở nồng độ cao, từ đó mô phỏng trực tiếp quá trình tiến triển xói lở theo thời gian mà không cần sử dụng các công thức vận chuyển bùn cát kinh nghiệm.

---

Cơ sở lý thuyết và phương pháp nghiên cứu

Mô hình toán học của luận án dựa trên hệ phương trình vi phân chi phối dòng chảy hai pha Euler-Euler, coi pha lỏng ($f$) và pha rắn ($s$) là hai môi trường liên tục thâm nhập lẫn nhau. Tỷ lệ thể tích của mỗi pha ký hiệu là $\alpha_t$ ($t = s, f$), thỏa mãn điều kiện ràng buộc $\alpha_s + \alpha_f = 1$.

1. Hệ phương trình chi phối

Phương trình bảo toàn khối lượng (liên tục) cho từng pha: $$\frac{\partial}{\partial t}(\alpha_t \rho_t) + \nabla \cdot (\alpha_t \rho_t \vec{v}_t) = 0$$

Phương trình bảo toàn động lượng cho pha lỏng: $$\frac{\partial}{\partial t}(\alpha_f \rho_f \vec{v}_f) + \nabla \cdot (\alpha_f \rho_f \vec{v}_f \vec{v}_f) = -\alpha_f \nabla P + \nabla \cdot \bar{\bar{\tau}}f + \alpha_f \rho_f \vec{g} + K{sf}(\vec{v}_s - \vec{v}_f)$$

Phương trình bảo toàn động lượng cho pha rắn: $$\frac{\partial}{\partial t}(\alpha_s \rho_s \vec{v}_s) + \nabla \cdot (\alpha_s \rho_s \vec{v}_s \vec{v}_s) = -\alpha_s \nabla P - \nabla P_s + \nabla \cdot \bar{\bar{\tau}}s + \alpha_s \rho_s \vec{g} + K{fs}(\vec{v}_f - \vec{v}_s)$$

Trong đó, $P$ là áp suất chung; $P_s$ là áp suất pha rắn sinh ra do va chạm và tiếp xúc hạt; $K_{sf} = K_{fs}$ là hệ số trao đổi động lượng liên pha (sử dụng mô hình lực cản đối xứng); $\bar{\bar{\tau}}_f$ và $\bar{\bar{\tau}}_s$ lần lượt là tensor ứng suất của pha lỏng và pha rắn.

Độ nhớt cắt của pha rắn $\mu_s$ được tổng hợp từ ba thành phần: độ nhớt va chạm ($\mu_{s,col}$), độ nhớt động học ($\mu_{s,kin}$) và độ nhớt ma sát ($\mu_{s,fric}$): $$\mu_s = \mu_{s,col} + \mu_{s,kin} + \mu_{s,fric}$$ Thành phần ma sát $\mu_{s,fric}$ được xác định từ lý thuyết thế dẻo (plastic potential theory) của Schaeffer (1987) nhằm phản ánh lực ma sát giữa các hạt cát ở vùng nồng độ cao sát đáy. Hàm phân bố xuyên tâm biểu diễn xác suất va chạm giữa các hạt khi thể tích đạt giới hạn $\alpha_{s,max} = 0.63$: $$g_{0,ss} = \left[1 - \left(\frac{\alpha_s}{\alpha_{s,max}}\right)^{1/3}\right]^{-1}$$

2. Mô hình rối và nhiệt độ hạt

• Pha lỏng: Sử dụng mô hình rối $k-\epsilon$ tiêu chuẩn của Launder và Spalding (1974), bổ sung các số hạng nguồn $\Pi_{k_f}$ và $\Pi_{\epsilon_f}$ phản ánh tác động trao đổi động lượng rối liên pha phụ thuộc vào vận tốc trôi $v_{dr}$ và hệ số khuếch tán rối $D_{t,fs}$.
• Pha rắn: Ứng dụng lý thuyết của Tchen (Hinze, 1975) về sự phân tán hạt rời rạc trong dòng rối đồng nhất, xác định động năng rối $k_s$ và độ nhớt xoáy $\nu_{t,s}$ thông qua tương quan giữa thang thời gian quán tính hạt $\tau_{F,fs}$ và thang thời gian tương tác xoáy - hạt $\tau_{t,fs}$.
• Nhiệt độ hạt ($\Theta_s$): Biểu thị động năng dao động ngẫu nhiên của các hạt cát, được giải bằng phương trình truyền tải suy dẫn từ lý thuyết động học chất khí đậm đặc của Ding và Gidaspow (1990), bao gồm các số hạng sinh năng do ứng suất rắn, khuếch tán $k_{\Theta_s} \nabla \Theta_s$, tiêu tán do va chạm không đàn hồi $\gamma_{\Theta_s}$ (với hệ số phục hồi $e_{ss}$) và trao đổi năng lượng vi mô với pha lỏng $\phi_{fs} = -3 K_{fs} \Theta_s$.

3. Phương pháp số và kỹ thuật tính toán

Mô hình được triển khai trong phần mềm CFD thương mại FLUENT kết hợp bộ tạo lưới GAMBIT. Để khắc phục hiện tượng suy giảm dòng cục bộ phi thực tế do tham số hóa tương tác hạt - rối quá mức và sự chênh lệch thang thời gian giữa dòng chảy và biến dạng đáy, tác giả sử dụng quy trình tính toán tách rời:

1. Giải trường vận tốc đơn pha ổn định của pha lỏng bằng phương trình Navier-Stokes và mô hình $k-\epsilon$.
2. Sử dụng trường vận tốc này làm đầu vào để chạy mô hình hai pha Euler-Euler nhằm tính toán sự dịch chuyển bùn cát và cập nhật biên dạng đáy.
3. Mặt phân giới nước - cát được xác định tại đường đồng mức phân số thể tích $\alpha_s = 0.5$ (dựa trên quan sát của Wang & Chien (1985) về hành vi chảy tầng của pha rắn khi $\alpha_s \ge 0.5$).
4. Tái tạo lưới tự động đáp ứng biên dạng đáy mới bằng cách thực thi tệp nhật trình (journal file) trong GAMBIT, sau đó lặp lại quy trình cho đến khi đạt trạng thái xói lở cân bằng.

---

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

Chương 1. Mô phỏng xói lở xung quanh đường ống cố định (Simulation of Scour Around a Fixed Pipeline)

Chương 1 tập trung mô phỏng quá trình xói lở nước trong (clear-water scour) xung quanh đường ống cố định nằm ngang sát đáy cát, đối chiếu trực tiếp với dữ liệu thực nghiệm chuẩn của Mao (1988).

          Dòng chảy vào (U_∞ = 0.31 m/s)
      
                     [ Đường ống (D = 0.1 m) ]
-----------------------·-·-·-·-·-·-·-·-·-·-·-·-·-·-·-·-------- Mặt phân giới cát (α_s = 0.5)
░░░░░░░░░░░░░░░░░░░░░ ( Hố xói sâu S_e ≈ 0.6D )   /\ Đụn cát

• Thông số thiết lập: Kênh dẫn dài $23\text{ m}$, rộng $2\text{ m}$, chiều sâu nước $H_0 = 0.4\text{ m}$. Đường ống có đường kính $D = 0.1\text{ m}$ đặt sát lớp cát dày $\delta_s \approx 1D$ ($0.1\text{ m}$) và $\delta_s = 1.5D$ ($0.15\text{ m}$), đường kính hạt cát $d_s = 0.36\text{ mm}$, tham số Shields ban đầu $\theta = 0.048$, vận tốc dòng tới $U_\infty = 0.31\text{ m/s}$. Lưới tính toán gồm $9.803$ nút và $9.575$ ô lưới ($105 \times 60$ ở pha nước, $105 \times 31$ ở tầng cát). Biên dòng vào cách tâm ống $5D$, biên dòng ra cách $15D$. Một nhiễu loạn dạng sóng sin biên độ $0.1D$ được đưa vào tại $t = 0$.
• Động lực học hố xói: Dòng chảy qua khe hẹp dưới đáy ống bị co thắt và gia tốc mạnh, cuốn trôi các hạt cát phía dưới và đẩy về phía hạ lưu hình thành đụn cát (mound). Cường độ rối chuẩn hóa $u'/w_s$ tại vị trí $1\text{ cm}$ trên đáy đạt giá trị cực đại ban đầu rồi giảm dần theo thời gian khi hố xói mở rộng, tiệm cận giá trị $u'/w_s \approx 1$ (bằng vận tốc lắng hạt $w_s$), duy trì trạng thái lơ lửng của cát.
• Độ sâu xói cân bằng: Độ sâu hố xói cực đại không thứ nguyên đạt $S_e/D \approx 0.6$, phù hợp với công thức thực nghiệm của Sumer & Fredsøe (1990) ($S_e/D = 0.6 \pm 0.1$) và công thức tiến triển thời gian $S = S_e [1 - \exp(-t/T)]$.

Thời gian $t$ (phút)	Trạng thái phát triển hố xói	Đặc điểm đụn cát hạ lưu	Cường độ rối $u'/w_s$ ($1\text{ cm}$ trên đáy)
$t = 10$	Xói khởi phát nhanh dưới đáy ống	Đụn cát hình thành sát sau ống	Đạt đỉnh do vận tốc dòng co thắt lớn
$t = 30$	Hố xói mở rộng về độ sâu và hạ lưu	Đụn cát cao dần, dịch chuyển ra xa	Bắt đầu suy giảm tại vùng dưới tâm ống
$t = 100$	Tốc độ xói giảm dần	Đụn cát tiếp tục di chuyển về hạ lưu	Phân bố rối mở rộng theo vết trượt dòng
$t = 200$	Tiệm cận trạng thái cân bằng	Đỉnh đụn cát hạ thấp và trải dài	$u'/w_s \approx 1$ duy trì cát lơ lửng
$t = 300$	Đạt cân bằng động ($S/D \approx 0.6$)	Đụn cát ổn định ở vị trí xa ống	Ổn định quanh giá trị vận tốc lắng $w_s$

• Nhận diện ba chế độ vận chuyển bùn cát:
  1. Tải đáy (Bed-load, $q_b$): Nằm trong lớp sát mặt phân giới có chiều dày $\Delta_b \approx 3d_s$, chuyển động do va chạm giữa các hạt, tính bằng $q_b = \alpha_s \Delta_b U_s$.
  2. Tải lơ lửng (Suspended-load, $q_s$): Nằm phía trên lớp tải đáy, các hạt được nâng đỡ bởi chuyển động rối của chất lưu, tính bằng tích phân $q_s = \int_{Y_0 + \Delta_b}^H U_s \alpha_s dY$.
  3. Tải phân lớp (Laminated-load): Nằm phía dưới mặt phân giới ($\alpha_s \ge 0.5$). Do ứng suất cắt truyền sâu và gradient áp suất âm hướng về đụn cát, các lớp hạt sát dưới bề mặt trượt theo dạng lớp tựa chảy tầng và hình thành một vùng tuần hoàn ngầm (recirculation zone) với vận tốc hạt tính toán đạt khoảng $0.2 U_\infty$.
• So sánh với mô hình NBURY và DRAMBUIE: Mô hình NBURY dự báo tốc độ xói ban đầu quá nhanh và hố xói sâu hơn số liệu mô phỏng hai pha. Mô hình DRAMBUIE dự báo thời gian đạt cân bằng kéo dài quá mức ($1.400\text{ phút}$ mới đạt cân bằng), cho thấy các mô hình bán kinh nghiệm này chưa mô tả chính xác động lực học xói cục bộ.

---

Chương 2. Mô phỏng xói lở dưới đường ống võng (Simulation of Scour Below a Sagging Pipeline)

Chương 2 mở rộng nghiên cứu sang điều kiện xói đáy động (live-bed scour) và hiện tượng đường ống tự hạ chìm (sagging) vào hố xói theo cấu hình thí nghiệm của Fredsøe et al. (1988).

• Thiết lập hai giai đoạn:
  1. Giai đoạn 1: Đường ống được giữ cố định sát đáy cát cho xói lở đáy động phát triển tự do trong $60\text{ phút}$ để tạo biên dạng hố xói ban đầu.
  2. Giai đoạn 2: Đường ống được điều khiển hạ dần theo phương thẳng đứng với các vận tốc võng không đổi $V_s = 1.1\text{ mm/phút}$, $3.3\text{ mm/phút}$ và $12.5\text{ mm/phút}$ cho đến khi chạm đáy hố xói.
• Tương tác thủy động và vận chuyển cát khi ống võng:
  • Khi ống bắt đầu hạ xuống, khoảng cách khe hở (gap) giữa đáy ống và hố xói thu hẹp, làm vận tốc dòng chảy trung bình qua khe hở gia tăng trở lại.
  • Vận tốc dòng tăng kích thích xói thứ cấp dưới đáy ống, duy trì lưu lượng bùn cát tải qua khe hở ngay cả khi ống đang đi xuống.
  • Vận tốc võng $V_s$ càng nhỏ ($1.1\text{ mm/phút}$), dòng chảy càng có đủ thời gian để quét sạch bùn cát đáy hố, tạo ra hố xói rộng và sâu hơn trước khi ống tiếp đáy. Ngược lại, ở vận tốc võng cao ($12.5\text{ mm/phút}$), ống chìm nhanh hơn tốc độ xói cục bộ, khiến khe hở bị bít kín sớm hơn.
• Hiện tượng phát tán xoáy (Vortex Shedding): Trong quá trình xói cố định và giai đoạn đầu của ống võng, sự phát tán xoáy luân phiên phía sau ống tạo ra các dao động áp suất và duy trì tải lơ lửng ở hạ lưu. Khi ống hạ sâu vào hố xói, thành hố xói hai bên ức chế cơ chế cuộn xoáy tự do, làm suy giảm cường độ phát tán xoáy.
• So sánh với mô phỏng LES của Cheng & Li (2003): Biên dạng hố xói và xu hướng phát triển độ sâu từ mô hình hai pha Euler-Euler đạt sự tương đồng cao với kết quả mô phỏng bằng phương pháp xoáy lớn (LES) của Cheng & Li, đồng thời khớp tốt với các mặt cắt đo đạc thực nghiệm của Fredsøe et al. (1988).

---

Kết quả và những đóng góp mới

1. Đóng góp về lý luận và phương pháp học thuật

• Ứng dụng lý thuyết hai pha hoàn chỉnh vào bài toán xói lở: Khác với các nghiên cứu truyền thống ghép nối phương trình dòng chảy một pha với công thức tải cát thực nghiệm rời rạc, luận án đã ứng dụng hệ phương trình vi phân hai pha Euler-Euler tích hợp lý thuyết động học chất khí đậm đặc và độ nhớt ma sát dẻo để mô tả đồng thời cả pha lỏng và pha rắn.
• Phát hiện và lý giải cơ chế tải phân lớp (Laminated-load): Luận án chứng minh rằng lớp cát nằm dưới mặt phân giới ($\alpha_s \ge 0.5$) không hoàn toàn bất động mà vẫn tham gia chuyển động dạng tầng dưới tác dụng của ứng suất cắt bề mặt và gradient áp suất, hình thành vùng tuần hoàn ngầm trong lớp trầm tích. Đây là cơ sở lý giải vì sao các mô hình kinh điển (như của Brørs) bị dừng phát triển xói lở sớm hơn thực tế.
• Quy trình giải số tách rời hiệu quả: Đề xuất quy trình giải xen kẽ trường vận tốc đơn pha ổn định làm đầu vào cho mô hình hai pha, kết hợp tái tạo lưới thích ứng tự động qua GAMBIT journal file, giúp triệt tiêu hiện tượng lắng đọng nhân tạo và khắc phục sự bất tương thích thang thời gian.

2. Đóng góp về thực tiễn và giải pháp kỹ thuật

• Đánh giá và chỉ rõ hạn chế của các mô hình vận hành quân sự: Cung cấp bằng chứng số trị chi tiết về sai số của hai mô hình dự báo chôn lấp thủy lôi NBURY (dự báo xói quá nhanh) và DRAMBUIE (dự báo thời gian cân bằng quá chậm), phục vụ định hướng cải tiến cho chương trình Mine Burial Prediction (MBP) của Hải quân Hoa Kỳ.
• Dữ liệu định lượng cho bài toán ổn định đường ống biển: Cung cấp bức tranh toàn diện về trường vận tốc, cường độ rối, lưu lượng bùn cát và quy luật biến thiên hố xói theo các tốc độ võng ống ($V_s = 1.1 - 12.5\text{ mm/phút}$), hỗ trợ kỹ sư công trình biển đánh giá rủi ro võng gãy và tối ưu hóa giải pháp chôn lấp bảo vệ đường ống.

---

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

1. Hạn chế được ghi nhận từ văn bản

• Độ lớn vận tốc trong vùng tuần hoàn ngầm: Vận tốc hạt trong vùng tuần hoàn ngầm tính toán được ($0.2 U_\infty$) tương đối cao so với trực giác vật lý, xuất phát từ việc thiếu dữ liệu đo đạc trực tiếp bên trong tầng cát để hiệu chỉnh tham số độ nhớt ma sát dẻo $\mu_{s,fric}$ cũng như tương tác hạt - rối ở nồng độ cao.
• Giả định hình học hai chiều (2D): Mô hình chưa xét đến các hiệu ứng ba chiều dọc theo trục ống như hiện tượng xói truyền dọc (3D scour propagation) hay sự uốn cong không gian của đường ống thực tế.
• Giả định dòng chảy điều chỉnh: Quy trình coi biên dạng đáy không đổi trong thời gian chất lưu điều chỉnh vận tốc (flow adjustment) là một phép đơn giản hóa có thể làm biến dạng nhẹ tốc độ xói tức thời.

2. Hướng nghiên cứu tiếp

• Tiến hành các thí nghiệm chuyên sâu đo đạc vận tốc hạt và ứng suất bên trong lớp cát dưới mặt phân giới để kiểm chứng định lượng cơ chế tải phân lớp.
• Hoàn thiện tham số hóa mô hình độ nhớt ma sát và tương tác rối liên pha trong dòng hai pha nồng độ cao.
• Mở rộng mô hình hai pha sang không gian 3D và bổ sung tác động kết hợp của trường sóng và dòng chảy dao động.

---

Giá trị tham khảo

• Đối với nghiên cứu sinh và nhà nghiên cứu thủy khí / công trình thủy: Cung cấp phương pháp luận chi tiết về việc ứng dụng lý thuyết dòng hai pha Euler-Euler, kỹ thuật đóng kín độ nhớt ma sát bằng lý thuyết thế dẻo và kỹ thuật tái tạo lưới thích ứng động trong phần mềm CFD thương mại (FLUENT/GAMBIT).
• Đối với kỹ sư thiết kế đường ống dầu khí và cáp ngầm biển: Cung cấp tập số liệu tham chiếu về trường dòng chảy qua khe hẹp, áp lực xoáy và sự phụ thuộc của hình học hố xói vào tốc độ võng ống, phục vụ tính toán dao động do xoáy (VIV) và khả năng tự chôn lấp của công trình.
• Đối với các chuyên gia hải quân và mô hình hóa địa kỹ thuật biển: Cung cấp các phân tích so sánh định lượng giúp nhận diện điểm yếu trong thuật toán của các mô hình chôn lấp thủy lôi thực chiến như NBURY và DRAMBUIE.

---

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

1. Mặt phân giới giữa nước và cát được xác định trong mô hình số trị bằng tiêu chí nào?

Mặt phân giới nước - cát được xác định tại đường đồng mức có phân số thể tích hạt cát $\alpha_s = 0.5$. Tiêu chí này dựa trên các quan sát thực nghiệm của Wang & Chien (1985), chỉ ra rằng khi $\alpha_s \ge 0.5$, các hạt cát liên kết chặt chẽ làm triệt tiêu chuyển động rối của chất lưu và pha rắn chuyển động theo dạng chảy tầng.

2. Tải phân lớp (laminated-load) hình thành do những cơ chế động lực học nào?

Tải phân lớp hình thành do ứng suất cắt bề mặt truyền sâu xuống lớp hạt qua va chạm liên hạt (thể hiện qua độ nhớt va chạm $\mu_{s,col}$) kết hợp với gradient áp suất âm sinh ra bởi dòng chảy co thắt hướng về phía đụn cát hạ lưu. Sự kết hợp này kích thích một lớp cát sát dưới mặt phân giới chuyển động và tạo thành vùng tuần hoàn ngầm trong trầm tích.

3. Tại sao tác giả phải sử dụng trường vận tốc đơn pha để hỗ trợ tính toán mô hình hai pha trong FLUENT?

Nếu giải trực tiếp mô hình hai pha hoàn chỉnh ngay từ đầu, các số hạng tương tác hạt - rối hiện hành trong phần mềm quá mạnh làm suy giảm vận tốc cục bộ phi thực tế, khiến hạt lắng đọng quá nhanh quanh ống. Đồng thời, thang thời gian điều chỉnh của dòng chảy qua miền tính toán không tương thích với thang thời gian tương tác hạt - rối. Việc giải trường vận tốc đơn pha ổn định trước giúp duy trì trường vận tốc thực tế và đảm bảo tính ổn định số học.

4. Kết quả mô phỏng độ sâu xói cân bằng dưới đường ống cố định đạt giá trị bao nhiêu và so với thực nghiệm như thế nào?

Độ sâu hố xói cân bằng không thứ nguyên đạt $S_e/D \approx 0.6$ đối với lớp cát có chiều dày $\delta_s = 1D$. Kết quả này hoàn toàn phù hợp với công thức tổng hợp thực nghiệm của Sumer & Fredsøe (1990) ($S_e/D = 0.6 \pm 0.1$) và dữ liệu đo đạc của Mao (1988).

5. Tốc độ võng của đường ống ($V_s$) ảnh hưởng như thế nào đến quá trình xói lở hố xói?

Vận tốc võng quyết định thời gian lưu của dòng co thắt qua khe hở đáy:

• Ở vận tốc nhỏ ($V_s = 1.1\text{ mm/phút}$), dòng chảy duy trì vận tốc quét đáy đủ lâu để tạo ra hố xói rộng và sâu trước khi ống tiếp đáy.
• Ở vận tốc lớn ($V_s = 12.5\text{ mm/phút}$), ống hạ xuống đáy hố nhanh hơn tốc độ bồi xói tự nhiên, làm khe hở nhanh chóng bị bít kín và hạn chế sự phát triển tiếp theo của hố xói.

---

Kết luận

Luận án của Zhihe Zhao đã ứng dụng mô hình dòng hai pha Euler-Euler kết hợp lý thuyết động học chất khí đậm đặc và mô hình rối $k-\epsilon$ cải tiến để mô phỏng quá trình xói lở cục bộ xung quanh đường ống ngầm cố định và đường ống võng. Nghiên cứu giải thích cơ chế vận chuyển bùn cát thông qua ba thành phần (tải đáy, tải lơ lửng và tải phân lớp) cùng sự xuất hiện của vùng tuần hoàn ngầm trong lớp cát dưới đáy hố xói. Kết quả mô phỏng độ sâu xói cân bằng ($S_e/D \approx 0.6$) và biên dạng hố xói đạt độ tương thích cao với các dữ liệu thực nghiệm của Mao (1988), Fredsøe et al. (1988) cũng như mô hình LES của Cheng & Li (2003). Đồng thời, luận án cung cấp các đối chiếu định lượng quan trọng chỉ ra sai số của các mô hình dự báo chôn lấp thủy lôi hiện hành NBURY và DRAMBUIE.

Trích đoạn nội dung tài liệu

NUMERICAL SIMULATION OF SCOUR AROUND FIXED AND SAGGING PIPELINES USING A TWO-PHASE MODEL by - Zhihe Zhao A Dissertation Presented in Partial Fulfillment _ of the Requirements for the Degree Doctor of Philosophy ARIZONA STATE UNIVERSITY December 2006 UMI Number: 3241370 INFORMATION TO USERS The quality of this reproduction is dependent upon the quality of the copy submitted. Broken or indistinct print, colored or poor quality illustrations and photographs, print bleed-through, substandard margins, and improper alignment can adversely affect reproduction. In the unlikely event that the author did not send a complete manuscript and there are missing pages, these will be noted. Also, if unauthorized copyright material had to be removed, a note will indicate the deletion.

® UMI UMI Microform 3241370 Copyright 2007 by ProQuest Information and Learning Company. All rights reserved. This microform edition is protected against unauthorized copying under Title 17, United States Code. ProQuest Information and Learning Company 300 North Zeeb Road P.

Box 1346 Ann Arbor, MI 48106-1346 NUMERICAL SIMULATION OF SCOUR AROUND FIXED AND SAGGING PIPELINES USING A TWO-PHASE MODEL Zhihe Zhao has been approved September 2006 Par Cue | ee— — OL aa 875772 `ae YX Supervisory Committee ACCEPTED: xui Mb Dean, Division of Graduate Studies ABSTRACT A key aspect of design and the maintenance of underwater pipelines is the assessment of local scour and its propagation. Scouring around objects placed on a sandy bottom is very complex because it involves two-phase turbulent flows and a myriad of sediment transport modes. This dissertation addresses two principal configurations of scour around pipelines in two parts. First, clear-water scour around a long fixed pipeline placed just above a non-cohesive sandy bed is numerically simulated.

Second, live-bed scour around a fixed pipeline and scour below a sagging pipeline are investigated. These two simulations are conducted by using an Eulerian two-phase model that implements Euler- Euler coupled governing equations for fluid and solid phases and a modified k-¢ turbulence closure for the fluid phase, the modeling system being a part of software FLUENT. Both flow-particle and particle-particle interactions are considered in the model. During the simulations, the interface between sand and water is specified using a threshold volume fraction of sand, and the evolution of the bedforms is studied in detail.

For clear-water scour around a fixed pipeline, the predictions of bedform evolution are in agreement with previous laboratory measurements. Investigations into the mechanisms of scour reveal that three sediment transport modes (bed-load, suspended-load and laminated-load) are associated with the scour development. While some previously proposed scour development formulae for cylindrical objects are in good agreement with the simulations, scour predictions based on some operational mine-burial models show disparities with present simulations. iii For investigations of live-bed scour and scour under a sagging pipeline, the flow and pipeline evolve in two steps: (1) the local live-bed scour around the pipeline developed around a fixed pipeline; and (2) the pipeline is lowered to the scour hole in controlled fashion until it reaches the bottom of the scour hole.

Three sagging velocities are simulated, and predicted scour profiles agree well with the laboratory data. General characteristics of flow fields, including turbulence, suspension of particles and sediment transport, are described paying attention to their dependence on pipeline sagging. Scour profiles simulated are also in agreement with a LES-based numerical study reported earlier. iv ACKNOWLEDGEMENT I would like to extend my thanks to my advisor Professor Harindra J.

Fernando for his guidance, patience and support on my research work. I specially thank him for the many hours he spent going through my drafts and giving me resourceful thoughts to improve my work. I also thank him for showing me what it takes to be successful in the graduate studies and beyond. I am thankful for Professor Don L.

Boyer for his insightful suggestions on my research work. In addition, it is also a privilege for me to have Professor Ronald Calhoun, Professor Kangping Chen and Professor Mark Schmeeckle as committee members. Finally, I appreciate the administrative supports from Ms. Gabrielle Stidham, Ms.

Jennifer McCulley and Mr. This work was funded by the U. office of Naval Research through its Coastal Geosciences and Mine Burial Programs.