7.1 Giới thiệu
Giao điểm giữa sông ngòi và đường giao thông đặt ra nhiều thách thức thiết kế khiến mô hình thủy lực 2D trở nên ưu việt. FHWA đã công bố tài liệu Two-Dimensional Hydraulic Modeling for Highways in the River Environment Reference Document (FHWA, 2019) nhằm hỗ trợ các kỹ sư thủy lực nắm vững việc ứng dụng mô hình 2D vào các dự án giao thông. Cùng với các tài liệu hướng dẫn riêng cho từng phần mềm mô hình, tài liệu tham khảo này là nguồn được khuyến nghị dành cho các kỹ sư áp dụng mô hình 2D cho các dự án liên quan đến giao thông trong môi trường sông ngòi.
Chương 4 trình bày các nguyên lý về những khái niệm cơ bản trong hầu hết các phép tính mô hình thủy lực thường gặp trong dòng hở và phân tích thủy lực cầu. Chương 5 cung cấp thông tin về cách chọn phương pháp mô hình thủy lực phù hợp nhất và tổng quan về cách xây dựng dữ liệu điều kiện biên đầu vào cùng các yếu tố cần cân nhắc khác trong các bài toán thủy lực cầu. Chương này cung cấp cái nhìn tổng quan về các chủ đề liên quan đến việc xây dựng hình học mô hình 2D chuyên biệt cho phân tích thủy lực cầu.
7.2 Các phương pháp số trong mô hình 2D
Ba phương pháp số phổ biến nhất được sử dụng để giải các phương trình vi phân riêng phần 2D là sai phân hữu hạn (finite difference), phần tử hữu hạn (finite element), và thể tích hữu hạn (finite volume). Mỗi phương pháp chia nhỏ (rời rạc hóa) miền mô hình thành các hình nhỏ, gọi là phần tử (elements) hoặc ô lưới (cells). Tập hợp các phần tử này tạo thành một lưới (mesh) hoặc mạng lưới biểu diễn miền mô hình. Nếu các phần tử có hình dạng và kích thước đồng nhất (ví dụ, hình vuông hoặc hình chữ nhật), thì được gọi là lưới có cấu trúc (structured grid). Lưới không có cấu trúc (unstructured grid) được sử dụng khi các phần tử có hình dạng và kích thước khác nhau (ví dụ, tam giác, tứ giác và các đa giác khác). Thuật ngữ “mesh” là tên gọi khác của lưới không cấu trúc.
7.2.1 Sai phân hữu hạn (Finite Difference)
Phương pháp sai phân hữu hạn là cách tiếp cận trực tiếp nhất để ánh xạ các phương trình vi phân thành các ô riêng lẻ. Thay vì biểu diễn phương trình như một hàm liên tục, các giá trị tại điểm đơn lẻ của các biến chính được dùng để ước lượng các thành phần trong các phương trình điều khiển. Việc tính toán thường được thực hiện tại tâm của mỗi ô trong một lưới đều đặn (lưới có cấu trúc). Kiểu lưới có cấu trúc này giúp việc giải phương trình hiệu quả và ổn định hơn, nhưng khó đại diện cho các đặc điểm địa hình có hình dạng bất quy tắc trong miền mô hình, như lòng sông hoặc các bờ đắp đường.
7.2.2 Phần tử hữu hạn (Finite Element)
Phương pháp phần tử hữu hạn sử dụng lưới không cấu trúc gồm các phần tử có hình dạng và kích thước đa dạng. Phương pháp này tích hợp một dạng của các phương trình điều khiển trên từng phần tử trong khi xét đến giá trị tại các phần tử lân cận để phát triển nghiệm. Cách tiếp cận này tạo ra nghiệm liên tục hoặc trơn trên mỗi phần tử. Việc sử dụng lưới không cấu trúc phù hợp hơn với các bài toán thủy lực liên quan đến cơ sở hạ tầng giao thông so với lưới có cấu trúc. Tuy nhiên, phương pháp phần tử hữu hạn có thể bị giới hạn về tính ổn định và có thể gặp các vấn đề về bảo toàn khối lượng so với phương pháp thể tích hữu hạn hoặc sai phân hữu hạn.
7.2.3 Thể tích hữu hạn (Finite Volume)
Phương pháp thể tích hữu hạn cũng sử dụng lưới không cấu trúc. Tuy nhiên, cơ sở số học của phương pháp thể tích hữu hạn thừa nhận rằng nhiều quá trình trong thực tế có thể được biểu diễn bằng định luật bảo toàn khối lượng và động lượng trong vật lý. Mỗi phần tử được xem là một thể tích điều khiển riêng biệt, và tốc độ thay đổi của các biến được tính tại mỗi cạnh phần tử. Phương pháp thể tích hữu hạn chỉ cần thực hiện các tính toán thông lượng tại các ranh giới phần tử. Phương pháp này đảm bảo bảo toàn khối lượng trong tất cả các ô lưới, giúp duy trì khối lượng cục bộ và toàn cục trong lưới. Phương pháp thể tích hữu hạn thường cho nghiệm ổn định hơn so với phương pháp phần tử hữu hạn, và thể hiện tốt hơn trong các điều kiện dòng chảy hỗn hợp thường gặp trong các tình huống lũ bị giới hạn như tại vị trí cầu.
7.2.4 Ẩn (Implicit) so với Tường minh (Explicit)
Các sơ đồ nghiệm số có thể là tường minh hoặc ẩn. Một sơ đồ tường minh giải trực tiếp bước thời gian kế tiếp bằng cách sử dụng giá trị biến đã biết tại bước thời gian hiện tại. Một sơ đồ ẩn sử dụng giá trị biến tại cả bước thời gian hiện tại và bước kế tiếp, và áp dụng ma trận hoặc kỹ thuật lặp để tìm nghiệm. Tài liệu của FHWA Two Dimensional Hydraulic Modeling for Highways in the River Environment (FHWA 2019) cung cấp giải thích chi tiết và ví dụ về các sơ đồ nghiệm ẩn và tường minh. Các mô hình sử dụng sơ đồ tường minh cần xét đến điều kiện Courant. Mục 8.4.4 của tài liệu HDS-7 này giải thích chi tiết về điều kiện đó.
7.3 Yêu cầu hình học đối với mô hình lưới không cấu trúc
7.3.1 Lưới (Meshes)
Mỗi mô hình số đều có giới hạn về loại lưới được sử dụng và loại phần tử cấu thành lưới. Lưới có cấu trúc thường hiệu quả hơn về tính toán nhưng có giới hạn trong cách chúng thể hiện các đặc điểm địa hình và mức độ chi tiết thủy lực do giới hạn về độ phân giải lưới. Dạng lưới có cấu trúc nhất là lưới Cartesian thuần túy với các phần tử có kích thước tương đương nhau (Hình 7.1). Dù vẫn giữ được sự đơn giản của một lưới, một số mô hình cho phép độ phân giải thay đổi bằng cách cho phép phần tử hình chữ nhật hoặc thay đổi kích thước. Lưới kiểu Quad-tree cung cấp sự linh hoạt lớn nhất trong nhóm lưới có cấu trúc. Dạng lưới này cho phép chia nhỏ thành các phần tử nhỏ hơn để tăng độ phân giải ở các khu vực có hình học phức tạp. Tuy nhiên, hướng sắp xếp của phần tử vẫn bị ràng buộc theo hướng của lưới.
Loại lưới linh hoạt nhất bao gồm một mạng biến đổi rộng rãi của các tam giác và tứ giác, còn gọi là lưới không cấu trúc. Các hướng dẫn chung về hình dạng phần tử được đưa ra nhằm tối đa hóa tính ổn định và hiệu suất số; tuy nhiên, không có ràng buộc tuyệt đối nào về hình dạng và kích thước phần tử. Dạng lưới này là một bề mặt phân mảnh đại diện cho mặt đất. Khi các điểm nhỏ lại (gần nhau hơn), khả năng biểu diễn được cải thiện. Một số mô hình hỗ trợ các tùy chọn phần tử bổ sung. Một ví dụ là phần tử thích nghi (adaptive elements), không yêu cầu phần tử chuyển tiếp như tam giác. Một ví dụ khác là phần tử đa giác Voronoi, cho phép có nhiều cạnh hơn.
7.3.2 Độ phân giải lưới (Mesh Resolution)
Một lưới thường bao gồm hàng nghìn đến hàng chục nghìn phần tử. Do số lượng phần tử có thể rất lớn, việc tạo lưới từ trước đến nay chiếm một phần đáng kể thời gian cần thiết để thực hiện mô phỏng mô hình 2D. Ngày nay, các công cụ và tiện ích hỗ trợ tạo lưới đã được cải tiến, giúp quá trình phát triển lưới bớt tốn công sức hơn.
Mức độ chi tiết mà lưới thể hiện các đặc điểm hình học và địa lý quan trọng trong miền mô hình ảnh hưởng đáng kể đến tính hợp lệ của kết quả tính toán từ mô hình số. Việc bỏ sót một yếu tố quan trọng như một kênh dẫn dòng có thể khiến mực nước trong mô hình cao hơn thực tế (xem Hình 7.2). Ngược lại, nếu không mô phỏng đúng một vật cản như bờ đắp chắn ngang một phần vùng ngập lũ, kết quả có thể là mực nước thấp hơn hoặc sự phân bố dòng không chính xác. Việc thể hiện chính xác các đặc điểm hình học là yếu tố then chốt của một lưới được xây dựng tốt. Quá trình tinh chỉnh thêm, hoặc tăng độ phân giải lưới, có thể cần thiết để mô phỏng chính xác điều kiện thủy lực.
Khi có sẵn dữ liệu địa hình đầy đủ, kỹ sư sẽ xem xét một số câu hỏi trong quá trình phát triển lưới, bao gồm: cần bao nhiêu phần tử và nên phân bố các điểm ra sao để thể hiện bề mặt? Câu trả lời cho các câu hỏi này phụ thuộc vào các công trình và cấu trúc có trong miền mô hình, biến đổi dòng chảy trong miền, và mục tiêu của mô hình. Trường dòng có thể thay đổi đáng kể trong một khu vực nhỏ, chẳng hạn như tại vị trí co hẹp qua cầu hoặc xung quanh trụ cầu. Một lưới phù hợp bao gồm đủ độ phân giải để thể hiện được các đặc điểm hình học quan trọng ảnh hưởng đến khả năng dẫn dòng và trường dòng/mực nước tính toán được.
Một lưới với phần tử rất nhỏ có thể thể hiện được tất cả các đặc điểm quan trọng nhưng lại tạo ra số lượng phần tử lớn, có thể dẫn đến thời gian tính toán và xử lý hậu kỳ dài hơn so với một lưới hiệu quả hơn. Cực đoan ngược lại là tạo ra một lưới có độ phân giải đủ để thể hiện các đặc điểm quan trọng nhưng không hơn. Bắt đầu từ một lưới tối thiểu như vậy, thực hành tốt là tinh chỉnh thêm ở những khu vực có biến đổi hình học hoặc thủy lực lớn, cho đến khi việc tinh chỉnh không còn làm thay đổi kết quả mong muốn, chẳng hạn như trường vận tốc qua một khẩu độ cầu. Cách tiếp cận này cân bằng giữa hiệu suất tính toán và chất lượng kết quả.
Các mô hình sử dụng độ phân giải dưới lưới (sub-grid), được mô tả tại Mục 7.4, đôi khi cho phép dùng lưới thô hơn. Tuy nhiên, ở các khu vực có biến đổi thủy lực lớn, ngay cả mô hình có sub-grid cũng cần các phần tử nhỏ hơn để thể hiện chính xác sự thay đổi nhanh của điều kiện thủy lực. Với mọi loại lưới và phần tử, việc căn chỉnh cạnh phần tử theo các đặc điểm thủy lực quan trọng — như chân và đỉnh của bờ đắp đường, bờ kênh, hoặc đê — là rất cần thiết.
7.3.3 Chất lượng lưới (Mesh Quality)
Trong lịch sử, các mô hình 2D thể hiện sự nhạy cảm số học đối với hình dạng phần tử. Nếu các phần tử không gần giống tam giác đều hoặc tứ giác đều, mô hình có thể phân kỳ và không tìm ra nghiệm. Các mô hình 2D thế hệ mới hiện nay không còn nhạy cảm với hình dạng phần tử như các mô hình cũ. Tuy vậy, việc kiểm tra chất lượng lưới vẫn là một thực hành tốt để nâng cao hiệu suất và độ tin cậy của mô hình.
Các mô hình sử dụng phần tử đều, phần tử thích nghi (adaptive), hoặc phần tử Voronoi thường có hình dạng và chất lượng lưới tốt. Lưới không cấu trúc yêu cầu cần được thiết lập cẩn thận để tránh các phần tử có thể gây bất ổn cho mô hình. Các hình dạng phần tử thường gây sự cố bao gồm các góc trong quá nhỏ hoặc quá lớn.
Vì mọi phương pháp số đều sử dụng lý thuyết ma trận ở một mức độ nào đó, nên chỉ cần một bất ổn nhỏ cũng có thể khiến mô hình không tìm được nghiệm. Do đó, việc kiểm tra chất lượng lưới là quan trọng. Một số phần mềm mô hình cung cấp công cụ để hỗ trợ đánh giá các vấn đề về chất lượng lưới. Có thể xem vấn đề chất lượng lưới tương tự như lỗi, cảnh báo và ghi chú trong mô hình HEC-RAS 1D. Thông lệ chuẩn là loại bỏ các lỗi và cảnh báo nghiêm trọng. Tuy nhiên, không phải tất cả các cảnh báo và ghi chú đều mang tính nghiêm trọng, và một số có thể không cần điều chỉnh mô hình.
7.4 Biểu diễn địa hình (Representing Terrain)
Dữ liệu địa hình và hình học được phát triển cho mô hình có ảnh hưởng lớn nhất đến độ chính xác của kết quả mô phỏng. Các yếu tố điều khiển thủy lực là những phần của địa hình có ảnh hưởng đến thủy lực. Các yếu tố điều khiển thủy lực bao gồm các đặc điểm dẫn dòng như kênh dẫn và các vùng cao như bờ đắp đường. Việc phát triển mô hình 2D bao gồm bước cần thiết là kiểm tra xem bề mặt địa hình số có thể hiện đúng các yếu tố điều khiển thủy lực phù hợp với mục tiêu của mô hình hay không. Bắt đầu với một bề mặt địa hình chính xác, kỹ sư có thể xây dựng được hình học lưới mô phỏng thể hiện đúng các yếu tố điều khiển thủy lực trong miền mô hình.
Các đặc điểm địa hình kiểm soát thủy lực ở lưu lượng thấp có thể không còn ảnh hưởng gì ở lưu lượng cao. Tương tự, một số đặc điểm địa hình kiểm soát thủy lực ở lưu lượng lớn có thể không ảnh hưởng gì ở lưu lượng nhỏ — ví dụ như các đặc điểm địa hình vùng bãi sông nằm ngoài phạm vi mặt nước khi lưu lượng thấp.
Ở bất kỳ kích thước phần tử nào, việc biểu diễn bề mặt địa hình đều có thể bao gồm một số sai số. Một số mô hình (ví dụ: HEC-RAS 2D) lưu trữ thông tin bổ sung dọc theo cạnh phần tử để bù cho sự xấp xỉ trong việc biểu diễn địa hình bằng lưới. Điều này thường được gọi là biểu diễn đặc điểm ở cấp độ dưới lưới (sub-grid scale feature representation). Thay vì chỉ lưu một cao độ duy nhất, phần tử sub lưới lưu trữ một cao độ tối thiểu và một đường cong biểu diễn mối quan hệ giữa cao độ và diện tích theo cạnh hoặc giữa cao độ và thể tích theo phần tử.
Tuy nhiên, cách tiếp cận cấp độ sub lưới này không thể hiện được các yếu tố điều khiển thủy lực bên trong phần tử. Do đó, mô hình không tính đến tác động mà các đặc điểm sub lưới có thể gây ra đối với hướng dòng chảy bên trong phần tử. Đây là lý do tại sao việc sử dụng các đường breakline — buộc các cạnh phần tử đi theo các yếu tố điều khiển thủy lực — là điều tối quan trọng khi xây dựng hình học mô hình.
Hình 7.3 minh họa hai phương pháp biểu diễn một bờ đắp trong lưới. Ở phía bên trái, một hàng cạnh phần tử duy nhất đại diện cho đỉnh bờ đắp. Ở phía bên phải, bờ đắp có hai hàng phần tử dọc theo đỉnh. Các hình phía dưới minh họa một mặt cắt ngang cắt qua bờ đắp. Đường nét đứt biểu diễn hình dạng thực tế của bờ đắp, đường liền biểu diễn mặt cắt lấy ra từ lưới, và các khối được tô bóng biểu diễn cao độ phần tử mà mô hình sử dụng.
Hầu hết các mô hình đánh giá dòng chảy qua các cạnh phần tử. Một số mô hình (ví dụ, SRH-2D) sử dụng cao độ tại tâm phần tử đồng thời cũng đánh giá các cao độ được xác định tại các cạnh để định nghĩa các rào cản thủy lực. Các đường đen giữa các khối đại diện cho phần rào cản của bờ đắp. Việc tăng độ phân giải lưới là hợp lý nếu cần có kết quả thủy lực chi tiết hơn để thể hiện chính xác mặt cắt mực nước và vận tốc, độ sâu nước tràn qua.
Lưới cần bao gồm các đặc điểm điều khiển thủy lực quan trọng như kênh dẫn và bờ đắp. Điều này không đồng nghĩa với việc phải tăng độ phân giải ở mọi nơi — hoặc thậm chí không cần tăng ở những nơi không cần thiết. Điều quan trọng là phải bố trí các đường breakline gồm một hàng các nút nơi tồn tại các yếu tố điều khiển. Việc lựa chọn độ phân giải một cách có hiểu biết để thể hiện các yếu tố này và đạt được độ phân giải mong muốn của kết quả thủy lực cho phép tạo ra lưới hiệu quả mà vẫn thể hiện đúng địa hình.
7.5 Biểu diễn sức cản dòng chảy
Giá trị Manning’s n, đã được đề cập trong Chương 4, là phương pháp phổ biến nhất dùng để tính đến sức cản dòng chảy do thảm thực vật và các đặc điểm bề mặt khác gây ra. Để gán các thông số độ nhám cho từng phần tử, các kỹ sư xây dựng một tập dữ liệu đa giác biểu diễn các loại hình sử dụng đất khác nhau trong toàn bộ khu vực dự án (Hình 7.4), sau đó gán các giá trị Manning’s n phù hợp cho các vùng đa giác được chỉ định.
Các tài liệu tham khảo về Manning’s n hiện có được phát triển để áp dụng cho các mô hình phân tích và mô hình 1D. Mô hình hai chiều (2D) vốn dĩ đã tính đến một số giả định được tích hợp trong hướng dẫn giá trị n hiện tại. Do đó, có thể suy luận hợp lý rằng các giá trị n trong mô hình 2D nên thấp hơn so với mô hình 1D. Các mô hình 2D đã được hiệu chỉnh cho thấy giá trị n có thể thấp hơn khoảng 10% so với giá trị dùng trong mô hình 1D. Tuy nhiên, hiện chưa có đủ nghiên cứu được công bố về chủ đề này để đưa ra kết luận dứt khoát. Ngoài ra, tiêu chuẩn hiện tại trong việc xác định giá trị n không đủ chính xác để đảm bảo rằng không cần hiệu chỉnh giữa mô hình 2D và 1D trừ khi có dữ liệu hiệu chỉnh chứng minh điều đó là cần thiết. Do đó, các giá trị tham khảo Manning’s n hiện có được trình bày trong Chương 4 là phù hợp để sử dụng cho cả mô hình 1D và 2D.
Hầu hết các mô hình 2D đều cho phép sử dụng Manning’s n thay đổi theo độ sâu, tương tự như các mô hình 1D. Nhiều loại hình sử dụng đất có thể dẫn đến sự biến đổi độ nhám tùy theo độ sâu dòng chảy. Tuy nhiên, hiện tại có rất ít tài liệu tham khảo sẵn có về chủ đề này. Việc áp dụng giá trị n thay đổi theo độ sâu đòi hỏi phải có đánh giá kỹ thuật.
Giá trị Manning’s n là đủ để đại diện cho sức cản dòng chảy đối với phần lớn miền mô hình. Tuy nhiên, đôi khi cần phải biểu diễn một số đặc điểm chi tiết hơn. Mô hình hai chiều bao gồm các phương pháp bổ sung để biểu diễn các công trình khác nhau. Một phương pháp là tạo một khoảng trống trong lưới như một bức tường không cho phép dòng chảy đi qua khu vực đó trong mô hình. Cách làm này tương tự như các vật cản bị chặn trong mô hình 1D. Một số mô hình cũng cho phép gán loại sử dụng đất “không xác định” giống như một khoảng trống. Cả hai phương pháp này đều phù hợp để mô phỏng các công trình xây dựng trong trường hợp cần thông tin chi tiết hơn về dòng chảy xung quanh công trình.
Áp dụng các vật cản là một tùy chọn khác mà một số mô hình cung cấp để biểu diễn sức cản dòng chảy. Vật cản bao gồm các đặc điểm điểm và đường nằm trong một hoặc nhiều phần tử. Các tham số đầu vào bao gồm chiều dài, chiều rộng và chiều cao. Mô hình sau đó tính toán lực cản liên quan đến vật cản và giảm khả năng dẫn dòng qua các phần tử chứa vật cản đó. Điều này biểu thị ảnh hưởng đến khả năng dẫn dòng trong phần tử do vật cản gây ra, nhưng không mô phỏng chi tiết các mô hình dòng chảy xung quanh vật cản cụ thể đó.
7.6 Mô hình hóa giao cắt dòng chảy với đường (Modeling Highway Crossings)
Một trong những lợi ích chính mà mô hình 2D cung cấp so với mô hình 1D hoặc các phương pháp phân tích đơn giản là khả năng mô tả chính xác hơn các mô hình dòng chảy qua cầu, gần cống, và trên đập tràn. Hầu hết các công trình đều có hình dạng phức tạp với các đặc điểm thay đổi đột ngột. Chúng cũng ảnh hưởng đến các dạng dòng chảy và có thể tạo ra các điều kiện dòng phức tạp và thay đổi nhanh chóng. Các điều kiện dòng phức tạp bao gồm chuyển tiếp giữa các dạng dòng, thay đổi không gian nhanh chóng về hướng dòng chảy, hình thành các mô hình xoáy, phân dòng, dòng chảy áp lực, dòng qua lỗ, dòng qua đập tràn, và các dạng khác.
Trong khi việc biểu diễn các đặc điểm tự nhiên đòi hỏi sự cẩn trọng khi xây dựng lưới, thì các công trình thủy lực lại mang đến những thách thức bổ sung. Các công trình như vậy bao gồm các tường đứng, các kết cấu kín như cống và cầu, và trụ cầu. Lưới mô hình hai chiều không thể trực tiếp biểu diễn tất cả sự phức tạp này trong hình học lưới, nhưng một số phương pháp trong các mô hình cho phép tích hợp các đặc điểm công trình.
Thông thường, các mô hình sử dụng ba phương pháp để biểu diễn công trình thủy lực, gồm:
- Xây dựng công trình trực tiếp trong lưới (biểu diễn tường minh – explicit representation).
- Đặc trưng hóa ảnh hưởng của công trình như một thuộc tính của phần tử chứa công trình, ví dụ như: thêm lực cản bởi vật cản trong phần tử, thuộc tính đặc điểm dưới lưới, hoặc các vùng áp lực buộc giới hạn mực nước tối đa dựa trên cao độ đáy dầm được người dùng nhập vào.
- Biểu diễn công trình bằng một mô hình riêng biệt (thường là 1D), hoặc phép tính riêng, sau đó liên kết mô hình đó, hoặc kết quả của nó, với mô hình tính toán 2D.
7.6.1 Biểu diễn nền đường và mố cầu
Biểu diễn các công trình một cách tường minh bằng các phần tử lưới 2D (phương pháp 1) là lựa chọn ưu tiên khi có thể thực hiện được. Một ví dụ phổ biến là dòng tràn tại bờ đắp đường dẫn đầu cầu. Mặc dù hình dạng dọc theo đỉnh của bờ đắp có thể được biểu diễn dưới dạng đập tràn 1D (phương pháp 3) trong miền mô hình 2D, việc mô hình hóa hình học của bờ đắp bằng nhiều phần tử lưới, theo cả chiều ngang và chiều dọc, sẽ cho kết quả mô phỏng chính xác hơn về đặc điểm hình học bên trong mô hình và tạo ra kết quả thực tế hơn.
Tương tự, việc mô hình hóa mố cầu trong lưới là cách tiếp cận phổ biến, như minh họa trong Hình 7.5. Cũng như với khu vực dưới cầu, tính toán dòng chảy 2D bao gồm cả dòng chảy xung quanh mố cầu. Việc đưa vào toàn bộ khu vực bờ đắp đường và mố cầu là hợp lý nếu mực nước trong mô phỏng vượt quá cao độ phần đường đầu cầu.
7.6.2 Biểu diễn trụ cầu
Mô hình hóa trụ cầu bằng cách tạo lỗ trong lưới (phương pháp 1), như minh họa trong Hình 7.5, hoặc bằng cách tính lực cản trong phần tử chứa trụ (phương pháp 2), là hai phương pháp phổ biến để mô hình hóa trụ cầu. Tăng độ phân giải lưới để đưa trụ vào trong định nghĩa hình học là hợp lý nếu cần mô tả trường dòng quanh trụ hoặc phân bố dòng qua khẩu độ cầu. Việc coi trụ như một lực cản bổ sung lên phần tử chứa nó có thể được áp dụng nếu kích thước phần tử lớn hơn kích thước trụ, hoặc nếu trụ có hình học phức tạp và việc mô hình hóa bằng một phần tử duy nhất là không đủ.
Vì các giả định thủy lực 2D và đặc tính ba chiều của dòng chảy thực quanh trụ, nên việc mô phỏng chi tiết dòng chảy vượt quá mức hai đến ba phần tử qua trụ là không hợp lý.
Một số mô hình không cho phép tạo lỗ trong lưới. Trong trường hợp đó, có thể tăng độ phân giải lưới quanh trụ và nâng cao độ các nút bên trong trụ để tạo hình học đại diện cho vật cản dòng chảy.
Phương pháp 2 phù hợp để phân tích vật cản trôi (debris) mắc vào trụ nếu cần. Bằng cách ước lượng kích thước và hình dạng khối vật cản, kỹ sư (hoặc một phần mềm mô hình) có thể xác định gần đúng lực cản và áp dụng giá trị n tương đương tăng lên tại các phần tử bị ảnh hưởng.
7.6.3 Biểu diễn mặt cầu, dòng chảy áp lực và dòng tràn qua mặt cầu
Mặt cầu ảnh hưởng đến thủy lực khi mực nước trong mô hình đạt đến phần kết cấu thấp nhất của công trình. Mặc dù lưới trong khu vực này của mô hình đã mô phỏng dòng chảy dưới mặt cầu, vẫn cần một phương pháp khác để biểu diễn chính phần mặt cầu. Việc này thường được thực hiện theo một trong hai cách:
- Gán thêm thuộc tính cho các nút hoặc phần tử trong khu vực cầu.
- Tách riêng vùng cầu ra khỏi mô hình và mô phỏng nó như một mô hình dòng riêng biệt.
Phương pháp đầu tiên tách ảnh hưởng của mặt cầu thành hai thành phần: dòng chảy áp lực dưới cầu và dòng chảy tràn qua cầu. Một số mô hình cho phép xác định cao độ trần cho từng nút hoặc phần tử dưới cầu. Hình 7.6 minh họa cấu trúc cầu phủ lên một lưới. Kỹ sư xác định cao độ đáy của mặt cầu cho từng nút dưới kết cấu đó, và mô hình sẽ tích hợp các thành phần áp lực vào các phương trình vi phân riêng phần cho dòng chảy nông. Một số mô hình còn có chức năng mô phỏng hiệu ứng qua lỗ (orifice) tại cửa vào vùng dòng chảy áp lực.
Nếu dòng chảy trong mô phỏng đạt đến đáy mặt cầu nhưng không vượt quá cao độ cần thiết để tràn qua, thì phương pháp này đủ để mô phỏng ảnh hưởng của mặt cầu. Thường thì, do độ nghiêng siêu cao của cầu hoặc chiều cao của phần kết cấu phía trên mặt cầu (lan can và rào chắn), tất cả dòng tràn xảy ra trên mặt đường chứ không phải trên mặt cầu. Trong các trường hợp như vậy, mô hình không cần thay đổi gì để mô phỏng dòng chảy qua mặt cầu.
Để mô phỏng dòng chảy kỳ vọng qua mặt cầu, người dùng cần xác định loại đập tràn và cao độ đỉnh để cho phép tính toán dòng chảy theo phương pháp 1D (sử dụng phương trình đập tràn), trong đó dòng chảy rời khỏi miền mô hình 2D tại mặt thượng lưu của cầu và quay trở lại miền mô hình 2D tại mặt hạ lưu.
Một tùy chọn khác trong một số mô hình là biểu diễn cầu như một vật cản. Cách tiếp cận này có thể tạo ra hiệu ứng backwater gắn với cầu một cách hiệu quả, nhưng nhìn chung cần hiệu chỉnh và không cung cấp thêm thông tin nào về dòng chảy vượt qua hoặc dòng chảy xuyên qua khẩu độ cầu. Phương pháp này dẫn đến vận tốc giảm tại vị trí vật cản vì đã bao gồm hiệu ứng của vật cản lên mực nước dâng, nhưng không xét đến ảnh hưởng của việc giảm diện tích khẩu độ đến vận tốc dòng. Điều này có nghĩa là kỹ sư sẽ phải tự xác định vận tốc thực tế qua khẩu độ cầu bằng các tính toán riêng biệt.
Có những lý do quan trọng để đưa cầu vào mô hình thủy lực. Lý do đầu tiên là để mô phỏng chi tiết dòng chảy và phân bố vận tốc gần cầu. Lý do thứ hai là để hiểu rõ đường đi và mô hình dòng chảy. Điều này rất cần thiết khi thiết kế cầu mới hoặc phân tích một cây cầu hiện hữu, phục vụ cho cả thiết kế thủy lực và đánh giá xói lở. Trong những tình huống này, lưới mô hình cần thể hiện các đặc tính vật lý của cầu. Không phải tất cả các phần mềm mô hình đều hỗ trợ các chức năng này. Đối với hầu hết các phân tích thủy lực liên quan đến giao thông, việc thể hiện hình học của các thành phần cầu và ảnh hưởng của chúng đến dòng chảy càng sát với thực tế càng tốt là điều lý tưởng. Việc lựa chọn một phần mềm mô hình có các tính năng và khả năng phù hợp với từng dự án là rất quan trọng. Đồng thời, cũng quan trọng không kém là phải hiểu rõ cách biểu diễn các thành phần của cầu đúng cách trong mô hình đã chọn.
7.6.4 Cống (Culverts)
Ba phương pháp được mô tả trước đó để mô hình hóa các công trình thủy lực cũng áp dụng cho cống. Cống có thể:
(1) là trọng tâm chính của mô hình 2D,
(2) là một yếu tố thủy lực được đưa vào mô hình, hoặc
(3) quá nhỏ để ảnh hưởng đáng kể nên không cần mô hình hóa.
Việc biểu diễn cống dưới dạng các yếu tố 1D trong mô hình 2D đã trở nên phổ biến. Với cách tiếp cận này, mô hình 2D tính toán mực nước (WSEs) tại hai đầu thượng lưu và hạ lưu của cống. Kết hợp mực nước hoặc đường năng lượng này với các đặc điểm của cống (số lượng ống, kích thước và hình dạng ống, kiểu cửa vào, loại vật liệu, v.v.), phần mềm sẽ tính toán dòng chảy qua cống. Kết quả trở thành điều kiện biên cục bộ trong mô hình 2D: dòng chảy rời khỏi miền mô hình tại đầu vào cống và quay trở lại tại đầu ra. Do cống có thể ảnh hưởng đến mực nước thượng và hạ lưu, việc tính toán thường là một quá trình lặp. Cách tiếp cận này về bản chất là mô hình 1D hoặc phân tích, giao tiếp với hoặc tích hợp vào mô hình 2D. Độ chính xác bị giới hạn bởi phương pháp mà mô hình sử dụng để tính toán dòng chảy qua cống. Ví dụ, phần mềm HY-8 của FHWA giả định điều kiện ao nước ổn định (pooling) phía thượng lưu cống và do đó bỏ qua vận tốc đầu vào.
Một lựa chọn khác trong mô hình hóa cống là biểu diễn chúng tương tự như cầu, như mô tả trong phần trên. Trong trường hợp này, người dùng đưa phần dưới hoặc đáy cống vào hình học lưới và mô phỏng phần trên của cống như mặt cầu (xem Hình 7.7). Các mặt đứng của cống có thể được biểu diễn như bờ dốc lớn trong lưới hoặc dưới dạng các lỗ, tùy theo phần mềm. Phương pháp này đòi hỏi nhiều công sức hơn khi xây dựng lưới, nhưng có ưu điểm là biểu diễn được dòng chảy 2D qua cống và tính đến vận tốc dòng ở thượng lưu.
Hỗ trợ duy trì trang:
Tôi xây dựng trang này để chia sẻ các tài liệu kỹ thuật cốt lõi trong thiết kế hạ tầng giao thông.
Nếu bạn thấy nội dung hữu ích và muốn góp phần duy trì trang hoạt động bền vững, tôi rất trân trọng mọi sự ủng hộ.