Mô hình hydist và phương pháp xử lý nồng độ phù sa tại biên lỏng đối với bài toán hai chiều

Bài viết trình bày một phương pháp xử lý nồng độ phù sa trên biên lỏng khi dòng chảy từ miền tính chảy ra biên trong mô hình dòng chảy hai chiều kết hợp đồng nhất, xen kẽ giữa tính toán thủy lực và chuyển tải phù sa - xử lý bằng phương pháp đường đặc trưng

BÀI BÁO KHOA HỌC MƠ HÌNH HYDIST VÀ PHƯƠNG PHÁP XỬ LÝ NỒNG ĐỘ PHÙ SA TẠI BIÊN LỎNG ĐỐI VỚI BÀI TOÁN HAI CHIỀU Nguyễn Thị Bảy1, Trần Thị Kim2, Huỳnh Cơng Hồi1, Phạm Anh Tài1, Nguyễn Đàm Quốc Huy2, Nguyễn Kỳ Phùng3 Tóm tắt: Bài báo trình bày phương pháp xử lý nồng độ phù sa biên lỏng dòng chảy từ miền tính chảy biên mơ hình dòng chảy hai chiều kết hợp đồng nhất, xen kẽ tính tốn thủy lực chuyển tải phù sa - xử lý phương pháp đường đặc trưng Tồn mơ hình tích hợp phần mềm HYDIST, tính áp dụng cho đoạn sông Tiền, thị trấn Tân Châu, tỉnh An Giang Kết phản ánh tính hiệu phương pháp, kết thu lan truyền phù sa từ thượng lưu sông Tiền thị trấn Tân Châu truyền khỏi miền tính hạ lưu, mà không tồn đọng lại miền, miền tính giới hạn, với phương pháp xử lý này, phù sa miền tính trao đổi tốt với miền ngồi chảy ngài miền tính, mà khơng thiết phải kéo dài miền tính Các kết cho thấy tính khả thi việc xử lý biên phương pháp đường đặc trưng áp dụng tính tốn lan truyền chất sơng, đặc biệt đoạn tính tương đối giới hạn Từ khóa: Mơ hình chuyển tải phù sa 2D, Xử lý biên, Phương pháp đường đặc trưng, Thủy lực Ban Biên tập nhận bài: 12/6/2019 Ngày phản biện xong: 24/7/2019 Giới thiệu Hiện nay, hầu hết nghiên cứu dòng chảy kết hợp với chuyển tải chất nói chung, phù sa, diễn biến lòng dẫn nói riêng mơ hình tốn phát triển theo ba hướng: (i) Áp dụng phần mềm xây dựng sẵn thương mại hóa Có thể kể đến phần mềm MIKE, xây dựng viện DHI, Đan Mạch Bộ phần mềm bao gồm nhiều module, sử dụng phổ biến Việt Nam nghiên cứu tính tốn động lực học dòng chảy Bên cạnh phần mềm CCHE 2, 3D Đại học Công Nghệ Mississippi biên soạn ứng dụng mơ q trình truyền thủy lực, chuyển động bùn cát lơ lửng, bùn cát đáy diễn biến lòng dẫn Ngồi nhiều phần mềm khác thương mại hóa, nhiên giá thành đắt đỏ Đại học Bách Khoa TP Hồ Chí Minh Đại học Tài Ngun Mơi trường TP HCM Viện Khoa học Cơng nghệ tính tốn HCM Email: ntbay@hcmut.edu.vn, kyphungng@gmail.com Ngày đăng bài: 25/08/2019 (ii) Hướng thứ hai hướng sử dụng phần mềm mã nguồn mở, tải miễn phí từ mạng Telemac (là chương trình viết ngơn ngữ Fortran biên soạn tập đoàn Điện lực Pháp), chương trình khơng cần quyền chưa có phần giao diện, việc khai thác sử dụng phức tạp, nên chưa phổ biến rộng rãi, mà thông qua số chuyên gia, tổ chức mở khóa tập huấn phần mềm Bộ phần mềm DELFT 2, 3D sản xuất Deltares - Hà Lan mơ 2, 3D cho dòng chảy, vận chuyển trầm tích, hình thái học,…được du nhập vào Việt Nam năm gần đây, nhiên chưa phổ biến rộng rãi (iii) Hướng tiếp cận thứ ba mô hình tính thủy lực hình thái sơng tự xây dựng nước, năm gần kể đến mơ hình HydroGIS PGS.TS Nguyễn Hữu Nhân có tích hợp cơng cụ GIS, demo kết giao diện tốt MK4 PGS TS Lê Song Giang với phần giao diện tốt phát triển nên chưa phổ biến rộng rãi Mơ TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 03 - 2017 57 BÀI BÁO KHOA HỌC hình HYDIST tác giả tự xây dựng từ năm 2004, đến dần hồn thiện Mơ hình tính dòng chảy 2D tích hợp sóng gió, kết hợp nhất, xen kẽ với chuyển tải phù sa diễn biến đáy Ban đầu, mơ hình biên nước chảy áp C dụng S  (C nồng độ phù sa, s phương x y) Tuy nhiên kết áp dụng điều kiện phù sa khơng trao đổi với bên ngồi vùng tính, thích hợp để giải tốn với miền tính lớn, để điều kiện biên khơng ảnh hưởng nhiều đến kết vùng tính Còn vùng tính giới hạn, điều kiện bộc lộ nhược điểm Với mục tiêu phát triển tiếp mô hình HYDIST, tác giả phát triển phương pháp xử lý biên phương pháp đường đặc trưng Phương pháp nghiên cứu thu thập tài liệu 2.1 Giới thiệu khu vực nghiên cứu đoạn sông đo đạc năm 1999 (Hình 1b) Trên Hình 3b biễu diễn mặt cắt qua trạm thủy văn Quốc gia Tân Châu Chế độ thủy lực khu vực Tân Châu thuộc chế độ chịu ảnh hưởng lũ từ thượng nguồn sơng Mekong đổ về, dòng chảy quanh năm hầu hết theo hướng từ thượng lưu đổ hướng biển (số thời gian dòng đổ ngược thượng lưu ít, đặc biệt mùa lũ, dòng chảy theo hướng chảy hạ lưu Trên Hình 2, đường đứt nét biểu diễn lưu lượng đo đạc theo trạm Tân Châu vào mùa lũ năm 1999 [7] Đây sở để tác giả chọn đoạn sông để tính tốn áp dụng mơ hình HYDIST sử lý biên phù sa hạ lưu phương pháp đường đặc trưng 2.2 Cơ sở lý thuyết mơ hình Mơ hình thiết lập dựa lời giải hệ phương trình chuyển động Reynolds [10], kết hợp với phương trình chuyển tải hai chiều sau: Hệ phương trình Reynolds: u u u  u  v2  u  v  g  Ku  A u t x y x h v  u  v2 v v  u  v  g  Kv  A v x y y h t    t x Hình Vị trí đoạn sơng Tiền qua thị trấn Tân Châu (a) địa hình đoạn sơng tính tốn (b) với đương đồng mức độ sâu ề 58 ấ Đoạn sông Tiền chảy qua thị trấn Tân châu nằm thượng nguồn sông Tiền Đồng sông Cửu Long Trên Hình giới thiệu vị trí đoạn sơng tính tốn (Hình 1a) độ sâu TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 08 - 2019 h   u    y h   v   (1) (2) (3) Trong u, v thành phần vận tốc trung bình theo độ sâu hai phương x y (m/s); cao độ mực nước so với chuẩn “0” (m); h độ sâu (m) Chuẩn “0” chọn mặt thoáng yên lặng, trục hướng lên; K hệ số ma sát đáy; A hệ số nhớt rối phương ngang (m2/s); g gia tốc trọng trường (m/s2) Phương trình chuyển tải: C    C  C C    C  S (4)  HK     u  v    HK x  y t  x y  H x  x  H y  y  H Trong C nồng độ phù sa trung bình theo chiều sâu (kg/m3); Kx, Ky hệ số phân tán phù sa theo phương phương x, y tính trung bình theo chiều sâu (m2/s); H độ sâu từ mặt thoáng xuống đáy (H=h+  ) (m); y hệ số phân BÀI BÁO KHOA HỌC bố vận tốc theo chiều sâu; S hàm nguồn, mô tả bốc lên hay lắng xuống hạt (kg/(m2.s)) Hàm số nguồn tính tự động chương trình dựa vào cơng thức thực nghiệm Van Rijn (1993) [8-9] S hàm số xói: M  b  e (5) S Khi    s b e S hàm số bồi: S  s C b d  b d e Khi  b   d (6) Và hàm số nguồn 0: S  Khi  e   b   d (7) Với s tốc độ lắng đọng hạt (m/s); M khả bốc hạt lên vùng tính (kg/m2/s);  b ứng suất tiếp đáy (N/m2) (8)   f U ( b fw  w b 8gn H1/ (9) Trong Ub vận tốc đáy (m/s);  khối lượng riêng nước, (kg/m3); Cb nồng độ phù sa đáy (kg/m3); fw hệ số ma sát đát đáy, tính theo Chezy (cơng thức (9)); n hệ số nhám; ,e , d ứng suất tới hạn xói bồi (N/m2) Các phương trình giải phương pháp sai phân hữu hạn, sơ đồ sai phân ẩn luân hướng ADI (Alternating Direction Implicit Method đề xuất Peaceman, Rachfor) [1], lưới sơ đồ trình bày Hình Hình Lưới sai phân theo sơ đồ ADI Trong thành phần u, v,  bố trí đặc biệt, cụ thể mực nước nồng độ phù sa C đặt tâm ô lưới (i, j), vận tốc u đặt vị trí (i+1/2, j) vận tốc v đặt vị trí (i, j+1/2) (với i, j = 1, 2, …) Chiều rộng ô lưới x  y Lưới đánh số số i (đối với phương x) từ đến N j (đối với phương y) từ đến M [6] 2.3 Điều kiện ban đầu điều kiện biên mơ hình 2.3.1 Điều kiện ban đầu Trong mơ hình, bắt đầu tính từ t0 = 0, tốn thủy lực gắn trạng thái tĩnh tồn miền, tốn chuyển tải phù sa gắn nồng độ số ban đầu Trong trường hợp tốn tính tiếp từ thời điểm t = t1 đó, điều kiện ban đầu trường vận tốc u,v(x,y) nồng độ C(x,y) thời điểm t1 toàn miền tính tốn 2.3.2 Điều kiện biên Biên lỏng: - Đối với tốn thủy lực: Hình Quỹ đạo hạt vật chất (di chuyển từ lớp thời điểm t đến t+ t Biên thượng nguồn dòng chảy chuỗi liệu lưu lượng đổ miền tính theo thời gian Q(t), biên hạ lưu dòng chảy cho dạng dao động mực nước  (t) - Đối với tốn chuyển tải phù sa: Khi dòng chảy hướng từ ngồi vào miền tính (thượng lưu), nồng độ phù sa biên gán nồng độ cho trước C(t), thường nội suy từ dãy giá trị đo đạc tính tương quan theo thời gian Riêng biên lỏng hạ lưu dòng chảy, mà dòng chảy từ miền tính hướng sử dụng  2C điều kiện S2  0, với S phương quỹ đạo mà hạt vật chất truyền Khi đó, nồng độ biên tính thơng qua q trình tải, q trình khuếch tán bỏ qua bước tính Nồng độ phù sa biên nghiệm phương trình truyền tải, giải theo phương pháp đường đặc trưng, trình bày kỹ mục 2.4 Biên cứng đường bờ: TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 08 - 2019 59 BÀI BÁO KHOA HỌC Áp dụng điều kiện biên khơng thấm tốn thủy lực với thành phần vận tốc vng góc bờ un = áp dụng điều kiện phản xạ toàn C phần toán chuyển tải phù sa: S  2.4 Xử lý nồng độ phù sa biên lỏng phương pháp đường đặc trưng Như trình bày trên, dòng chảy từ miền tính chảy biên lỏng, phù sa tính trước miền ảnh hưởng tồn lên biên, đặc biệt mùa lũ, mà vận tốc dòng chảy qua biên lớn, lúc chuỗi liệu biên C(t) cho trước biên trở nên khơng hiệu sử dụng làm điều kiện biên tính tốn cho tốn chuyển tải phù sa Lúc việc bỏ qua trình khuếch tán lại trở nên hữu hiệu  2C Theo đó, người ta chấp nhận biên:  S Bài toán chuyển tải giải luân hướng theo sơ đồ ADI (đã trình bày trên) theo phương x, sau phương y Trên phương, nồng độ biên giải từ phương trình truyền tải [2]: C C S U  t  (x, y) H (10) C( B , t + t ) = f(C( A , t)) (11) Hình mơ tả vết hạt vật chất từ A (lớp thời điểm trước) đến vị trí biên lỏng B (lớp thời điểm sau) Trong trình tải túy, xem hạt vật chất chất lỏng mang với vận tốc không đổi khoảng thời gian  t, nồng độ biên B hàm số tính theo nồng độ A: Để xác định nồng độ điểm B biên, cần ngược thời gian theo quỹ đạo điểm A, ta biết nồng độ (nồng độ A biết trước nội suy từ nồng độ nút lân cận): Q trình vừa mơ tả sở lý thuyết phương pháp đường đặc trưng: Trong phương trình (10), thay: U  dx dt Ta có: C C dx S (12)   hay: 60 t x dt dC S  dt H H TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 08 - 2019 (13) Lời giải phương trình (13) với điều kiện nồng độ A B sau: S (14) C B  C A  t H Trường hợp xem khoản thời gian dt, bốc lên hay lắng đọng hạt vật chất lân cận biên không đáng kể (S=0), ta có : (15) CB = CA Trong mơ hình HYDIST, điều kiện (15) áp dụng để xác định nồng độ phù sa biên lỏng (CB) dựa vào nồng độ chân đường đặc trưng CA dòng chảy hướng từ miền tính Trình tự tính: + Bước 1: Trong khoản thời gian  t, chấp nhận vận tốc số, cụ thể mơ hình tính trung bình cộng vận tốc hai lớp thời gian n n+1, từ xác định ngược lại vị trí điểm A chân đường đặc trưng + Bước 2: Nội suy giá trị nồng độ chân đường đặc trưng (tại A) theo giá trị biết nút + Bước 3: Tính giá trị B lớp thời gian n+1: CB = CA Các tốn kiểm định mơ hình HYDIST tính toán thủy lực truyền tải phù sa thực trình bày báo đăng [4-5] Sau xin giới thiệu kết ứng dụng tính tốn dòng chảy chuyển tải phù sa cho đoạn sông Tiền chảy qua thị trấn Tân Châu Các thông số áp dụng mơ hình hiệu chỉnh kiểm định mơ hình 3.1 Các thơng số đầu vào mơ hình Như trình bày mục 2.1, đoạn sơng Tiền qua thị trấn Tân Châu chọn để áp dụng tính mơ hình Lưới tính đoạn sơng lưới vng góc có kích thước: 464x324 phần tử, bao gồm phần đất liền lẫn lòng sơng, x =  y = 10m; bước thời gian tính  t = 2s Biên lỏng miền tính: dạng chuỗi lưu lượng nồng độ phù sa Q(t), C(t) mùa lũ 1999, từ đầu mùa lũ 10/5/1999 đến hết tháng 12/1999 Biên lỏng miền tính: dao động BÀI BÁO KHOA HỌC mực nước theo thời gian tương ứng, cho từ số liệu đo đạc trạm Tân Châu dời xuống trễ Chương trình tính tốn viết ngơn ngữ VB6, để tiết kiệm thời gian tính, module tính tốn, qt theo phương x y, có chương trình nhận dạng tự động biên bờ, để giới hạn biên đầu biên cuối tính tốn, tính tốn đến lòng sơng Trong chương trình, để tổng quát cho trường hợp nước chảy ngược lên, nhập chuỗi phù sa C(t) cho sẵn biên Tuy nhiên khí tính tốn trực tiếp cho đoạn sơng Tiền (Tân Châu), dòng chảy theo chiều từ thượng lưu tới hạ lưu, nên chương trình tự động xử lý lại biên phương pháp đường đặc trưng trình bày (mục 2.4), mà khơng sử dụng chuỗi liệu phù sa biên cho sẵn 3.2 Hiệu chỉnh kiểm định mơ hình Thơng số hiệu chỉnh mơ hình hệ số nhám Hệ số nhám vùng tính hiệu chỉnh thay đổi tỷ lệ nghịch với độ sâu lòng dẫn, dao động khoảng từ 0,005 đến 0,06 tương ứng với độ sâu lòng sơng thay đổi từ 41m đến 0,1m Khi tính tốn, lưu lượng trích xuất phục vụ kiểm định hiệu chỉnh vị trí trạm Tân Châu để đảm bảo chế độ thủy lực Kết tính lưu lượng cho mùa lũ 1999 từ mơ hình trạm Tân Châu so với kết thực đo xác (trên Hình 4, đường liền nét kết tính tốn, đường đứt nét kết đo đạc từ ngày 5/7/1999 đến cuối tháng 12/1999) Hình Lưu lượng tính từ mơ hình tốn thực đo vào sườn lũ lên từ 10/5/1999 đến 24/6/1999 trạm Tân Châu Hình Kết phân bố hệ số nhám n sông Tiền-Tân Châu sau hiệu chỉnh TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 08 - 2019 61 BÀI BÁO KHOA HỌC Hình Kết tính tốn trường vận tốc đoạn sơng Tiền (Tân Châu) lúc đỉnh lũ lúc 0h 6/10/1999 (a) (b) (d) (e) (f) (g) (h) (i) (j) 62 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 08 - 2019 (c) Hình Lan truyền phù sa đoạn sơng Tiền (Tân Châu) kết tính tốn bắt đầu tính từ lúc 10 ngày 10/5/1999 a)10h30’; b)11h30’; c)12h30’; d)13h30’; e)14h30’; f)15h30’; g)16h30’; h)17h30’; i)18h30’; j)19h30’ BÀI BÁO KHOA HỌC Trên Hình trình bày phân bố hệ số nhám miền tính sau hiệu chỉnh, với màu thay đổi từ nhạt (n = 0,005) tâm hố sâu, đến đậm (n = 0,06) vùng sát bờ Thảo luận trường vận tốc kết lan truyền phù sa miền tính Trường vận tốc đoạn sơng áp dụng tính lúc đỉnh lũ (ngày 0h ngày 6/10/1999) trình bày Hình Đối chiếu với độ sâu miền tính (hình 3), vị trí mặt cắt qua trạm thủy văn Tân Châu), vận tốc lớn gần trùng với vách hố sâu, đạt giá trị lên đến 3,18 m/s (thể vùng véc tơ sẫm màu nhất) Hình ảnh trường vận tốc thể rõ dòng chảy qua vùng khuỷu tay gấp khúc Tân Châu bám sát vách hố Để phân tích đánh giá lan truyền phù sa từ thượng lưu vào miền tính khỏi miền tính phía hạ lưu, kết lan truyền phù sa theo thời gian xuất tương ứng với thời điểm sau lúc bắt đầu tính lan truyền, vẽ thành hình ảnh đường đồng mức có tơ màu theo độ lớn giá trị phù sa (Hình 7a-7j) cách 1g từ ngày tính (10h ngày 10/5/1999) Trên hình 7a thể phân bố phù sa sau 30 phút tính tốn lan truyền (lúc 10h30’), lúc thấy phù sa bắt đầu truyền vơ vùng tính đạt giá trị lớn đầu vào lạch sâu (đạt tới lớn 0,008 kg/m3) Sau phù sa tiếp tục truyền sâu vơ miền tính, đến 11h30’ 12h30’ (Hình 7b-7c), ta thấy phù sa lan truyền vơ đến phân nửa miền tính Các kết dạng đường đồng mức phù sa hình 7d7f rõ phù sa truyền đến tận biên di chuyển khỏi vùng tính sau (hình 7g-7j) với nồng độ phù sa đạt lên đến 0,007 kg/m3 khu vực biên hạ lưu Nhìn vào dạng đường đồng mức phân bố phù sa miền tính, nhận thấy phù sa theo dòng chảy khỏi miền tính Các kết nói lên tính ưu việt việc xác định biên phù sa phương pháp đường đặc trưng trình bày mục 2.2, miền tính giới hạn phạm vi nhỏ (so với tồn hệ thống sơng Tiền sông Hậu) phù sa truyền hợp lý từ thượng lưu tới hạ lưu khỏi miền tính, mà khơng tồn đọng miền Nó cho thấy có trao đổi tốt với bên ngồi miền tính nước chảy Kết luận Việc áp dụng phương pháp đường đặc trưng xác định nồng độ phù sa biên lỏng khí dòng chảy khỏi miền tính tác giả áp dụng để tính tốn phù sa biên cách hiệu quả, giúp chương trình tính tốn thực dễ dàng nhanh gọn Mơ hình HYDIST kiểm tra tốt với tóan giải lý thuyết phòng thí nghiệm, áp dụng tính tốn thực tế đoạn sơng, kết tính khả thi mơ hình áp dụng tính tốn thủy lực, chuyển tải hai chiều đoạn sông khác Lời cảm ơn: Nghiên cứu tài trợ đề tài mã số: NĐT.28.KR/17 khuôn khổ thực nhiệm vụ khoa học cơng nghệ theo Nghi Định thư Nhóm tác giả chân thành cảm ơn hỗ trợ Tài liệu tham khảo Nguyễn Kỳ Phùng, Nguyễn Thị Bảy (2007), Mơ hình hóa chất lượng nước mặt, NXB Đại học Quốc gia Tp HCM Nguyễn Tất Đắc (2005), Mơ hình tốn cho dòng chảy chất lượng nước hệ thống kênh sông, Nhà xuất Nông nghiệp Tp HCM Nguyễn Thị Bảy (2004), Báo cáo tổng kết đề tài Khoa học - Công nghệ cấp Bộ: “Nghiên cứu dòng bùn cát, phù sa bồi lắng nó, Đề tài MS: B2000- 20- 82, Tp.HCM Nguyễn Thị Bảy, Mạch Quỳnh Trang (2006), Mơ hình tốn tính chuyển tải bùn cát kết dính TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 08 - 2019 63 BÀI BÁO KHOA HỌC vùng ven biển Phần 1: Mô hình tốn, Tạp chí Phát triển Khoa học Cơng nghệ Đại học QG Tp HCM, 9-2/2006, 53-60 Nguyễn Thị Bảy, Mạch Quỳnh Trang (2006), Mơ hình tốn tính chuyển tải bùn cát kết dính vùng ven biển Phần 2: Áp dung tính tốn mơ dòng bùn cát vùng ven biển Cần Giờ Tạp chí Phát triển Khoa học Công nghệ, 9-4/2006, 31-40 Nguyen, T.B., Tran, T.T., et al (2011), Numerical investigation on the sediment transport trend of Can Gio coastal area (Southern Vietnam) Journal of Marine Environmental Engineering, 9, 191-210 Trung tâm dự báo Khí tượng Thủy văn An Giang (2000), Báo cáo nghiên cứu thủy văn phục vụ xây dựng kè sông Tiền chảy qua thi trấn Tân Châu-An Giang Van Rijn, L.C., (1989), Handbook of Sediment transport by currents and waves Vols and 2, Delft Hydraulics, Delft Van Rijn, L.C., (1993), Principles of Sediment Transport in rivers, estuaries and coastal seas, Delft Hydraulic June 10 Vorobiev, V.N., Smimov, N.P., (2006), Hải dương học đại cương: Phần - Các trình động lực học (Phạm Văn Huấn dịch), NXB Đại học Quốc gia Hà Nội HYDIST MODEL AND THE APPROACH OF SOLVING SEDIMENT CONCENTRATION AT OPEN BOUNDARIES Nguyen Thi Bay1*, Tran Thi Kim2, Pham Anh Tai1, Nguyen Dam Quoc Huy2, Huynh Cong Hoai1, Nguyen Ky Phung3* Ho Chi Minh City University of Technology Ho Chi Minh University of Natural Resources and Environment Institute for Computational Science and Technology Abstract: This paper presents a method that solves the sediment concentration at the open boundaries when the flow direction from the computational domain to open boundaries in a two-dimensional area - characteristic curves method It combines alternatingly the hydraulic and sediment transport modules All of these features are built and available in HYDIST model and applied for calculation in a branch of Tien River passing Tan Chau Town, An Giang Province The results showed the effectiveness of the method, the results of spreading sediment from upstream of Tien river, Tan Chau town were transferred out of the cpmputational area Although the study is limited, with this treatment the sediments in the area are exchanged well with the outside area and flow out into the domain, without necessarily extending the domain The results showed the feasibility of boundary processing by the characteristic curves method in the application of sediment transport in the river, especially when the segment is relatively limited Keywords: 2-D Sediment transport module, processing the open boundaries, characteristic curves method, hydraulics 64 TẠP CHÍ KHÍ TƯỢNG THỦY VĂN Số tháng 08 - 2019 ... điều kiện biên không thấm toán thủy lực với thành phần vận tốc vng góc bờ un = áp dụng điều kiện phản xạ toàn C phần toán chuyển tải phù sa: S  2.4 Xử lý nồng độ phù sa biên lỏng phương pháp đường... mức phù sa hình 7d7f rõ phù sa truyền đến tận biên di chuyển khỏi vùng tính sau (hình 7g-7j) với nồng độ phù sa đạt lên đến 0,007 kg/m3 khu vực biên hạ lưu Nhìn vào dạng đường đồng mức phân bố phù. .. lưu dòng chảy cho dạng dao động mực nước  (t) - Đối với toán chuyển tải phù sa: Khi dòng chảy hướng từ ngồi vào miền tính (thượng lưu), nồng độ phù sa biên gán nồng độ cho trước C(t), thường