Tổng quan luận án

Dòng chảy một chiều trong lòng dẫn hở đóng vai trò thiết yếu trong công tác phát triển nguồn nước, tính toán thủy lực cho sông, hồ, biển và bảo vệ môi trường, đặc biệt trong các bài toán dự báo dòng chảy kiệt và dòng chảy lũ khi chưa tràn bờ. Hệ phương trình vi phân Saint-Venant (1871) là công cụ nền tảng được ứng dụng phổ biến để mô phỏng dòng chảy hở một chiều không ổn định. Tuy nhiên, hệ phương trình cổ điển này được xây dựng trên giả thuyết vận tốc phân bố đều trên mặt cắt ngang, phân bố áp suất tuân theo quy luật thủy tĩnh và gia tốc theo phương đứng không đáng kể.

Khoảng trống nghiên cứu được tác giả xác định là sự thiếu vắng mô hình toán một chiều có khả năng mô tả dòng chảy chịu ảnh hưởng của trường trọng lực khi xuất hiện thành phần vận tốc tương đối lớn theo phương thẳng đứng ở đáy lòng dẫn. Trong thực tế kỹ thuật thủy lợi và thủy văn, hiện tượng này xuất hiện tại các khu vực lòng dẫn có nước trồi, nước ngầm có áp phun lên từ đáy (như dòng chảy mùa kiệt ở kênh vùng trung du, kênh đi qua lòng hồ Phú Ninh), kênh có bố trí đập ngầm để nâng cao trình mặt nước phục vụ tưới tự chảy, hoặc dòng chảy tràn qua đê biển khi đoạn bờ biển thẳng. Các phần mềm một chiều thương mại hiện hành (như HEC-RAS, MIKE11, KOD1, SAL) chỉ cho phép khai báo nguồn bổ sung nhập lưu dạng ngang (lateral flow) với vận tốc nhỏ, vốn không phản ánh đúng bản chất vật lý của vận tốc đứng tại đáy – yếu tố gây cản trở dòng chảy chính và làm dâng mực nước thượng lưu rõ rệt. Mặt khác, việc giải quyết bằng các mô hình hai chiều đứng (2DV) hoặc ba chiều (3D) bằng phần mềm như ANSYS lại đòi hỏi khối lượng tính toán lớn, cấu trúc phức tạp và phương pháp thể tích hữu hạn (FVM) có bậc chính xác thấp hơn.

Mục tiêu nghiên cứu của luận án được xác định gồm ba nội dung:

  1. Xây dựng phương trình một chiều (1D) suy rộng, tổng quát hơn phương trình 1D cổ điển, cho phép mô tả dòng chảy có tốc độ theo phương thẳng đứng ở đáy lòng dẫn; đáp ứng yêu cầu tính toán cho các bài toán thực tế như lòng dẫn có nước trồi hoặc có vật nhô cao ở đáy.
  2. Áp dụng phương pháp phần tử hữu hạn Taylor-Galerkin và lập trình bằng ngôn ngữ Fortran 90 để thiết lập chương trình tính giải hệ phương trình một chiều suy rộng.
  3. Kiểm nghiệm thuật toán và chương trình tính đã thiết lập thông qua đối chiếu với số liệu thực nghiệm trên mô hình vật lý máng thủy lực.

Đối tượng nghiên cứu của luận án là dòng chảy hở một chiều.

Phạm vi nghiên cứu:

  • Không gian nghiên cứu: Tập trung vào dòng chảy hở một chiều trong lòng dẫn kênh, sông có điều kiện biên vận tốc thẳng đứng tại đáy; thực nghiệm kiểm chứng được tiến hành trên mô hình vật lý máng thủy lực.
  • Nội dung và phương pháp: Thiết lập hệ phương trình suy rộng thông qua tích phân hệ phương trình vi phân chuyển động hai chiều đứng (2DV); giải phương trình vi phân bằng phương pháp số phần tử hữu hạn Taylor-Galerkin đạt bậc chính xác bậc 3 theo thời gian; xây dựng phần mềm tính toán chuyên dụng và kiểm chứng qua thực nghiệm mô hình vật lý.

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

Văn bản luận án đã tổng hợp hệ thống các công trình nghiên cứu trong nước và quốc tế liên quan đến lý thuyết dòng chảy kênh hở và các phương pháp giải số, bao gồm:

  1. Lý thuyết cơ bản và các dạng xấp xỉ dòng chảy một chiều:
  • Saint-Venant (1871): Đưa ra hệ hai phương trình vi phân đầy đủ mô tả chuyển động không ổn định thay đổi chậm một chiều trên cơ sở các định luật bảo toàn khối lượng, động lượng và năng lượng.
  • Bakhmeteff (1932): Đề xuất hệ thống phân loại đường mặt nước cho dòng chảy ổn định không đều dựa trên độ dốc đáy, độ dốc phân giới, năng lượng đơn vị và số Froude.
  • Hayami (1951): Đề xuất phương pháp xấp xỉ sóng khuếch tán một chiều để diễn toán dòng chảy lũ không ổn định.
  • Chow (1971), Strelkoff (1969), Yen (1973, 1984): Thực hiện tích phân phương trình Navier-Stokes trên mặt cắt/chiều sâu để rút ngắn khoảng cách giữa cơ học chất lỏng và thủy lực kênh hở, làm rõ sự khác biệt giữa phương trình động lực và phương trình năng lượng.
  • Daluz Vieira (1983): So sánh lời giải số phương trình Saint-Venant với các phép xấp xỉ sóng động học, sóng khuếch tán và trọng lực dưới hai điều kiện biên hạ lưu.
  • J. Cunge và cộng sự (1981), Croley (1977): Hệ thống hóa thủ tục tính toán dòng chảy kênh hở và phát triển các mô hình thủy văn, thủy lực.
  1. Các nghiên cứu chuyên sâu về phân bố lưu tốc, lực cản và cấu trúc dòng chảy:
  • Xingwei Chen và Yee Meng Chiew (2004): Nghiên cứu lý thuyết và thực nghiệm phân bố vận tốc của dòng chảy rối trong kênh hở có dòng thấm vào đáy, đề xuất phân bố theo luật logarit sửa đổi.
  • Xavier Litrico và cộng sự (2005, 2006): Ứng dụng phương trình Saint-Venant kiểm soát tự động hồ kênh thủy lợi và điều tra chế độ dao động do sóng phản xạ trên kênh thí nghiệm tại Bồ Đào Nha.
  • Hitoshi Sugiyama và cộng sự (2006): Phân tích số dòng chảy rối trong kênh hỗn hợp uốn cong bằng mô hình đại số ứng suất Reynolds.
  • Jing Yan và cộng sự (2011): Khảo sát hiện tượng nhịp vận tốc trong dòng chảy hở qua thí nghiệm máng và máy đo Doppler Laser.
  • Zhihua Xie, Binliang Lin và nhóm cộng sự (2012): Phân tích biên độ xung vận tốc tức thời, số Reynolds và Froude liên quan đến quá trình xói mòn bùn cát và sạt lở đê kè.
  • Rashwan (2013): Thiết lập phương trình nước nhảy trong kênh hở hình bán nguyệt dựa trên phương trình động lượng.
  • Volpe Plantamura (2014): Áp dụng mô hình Entropy để xác định phân bố vận tốc mặt cắt ngang sông.
  • Xiao-guang Liu và Yu-hong Zeng (2016): Xác định hệ số cản của thực vật cứng trong lòng dẫn sông ngòi.
  • Agua và cộng sự (2017): Cải tiến thuật toán tính tốc độ dòng chảy dựa trên cơ cấu thủy lực và phương trình Saint-Venant 1D, nâng cao độ chính xác đo lường.
  • Vito Ferro (2018): Đề xuất phương trình kháng dòng chảy cho kênh đáy sỏi dựa trên biểu đồ vận tốc - năng lượng qua 416 thử nghiệm và 83 phép đo kiểm chứng.
  • Preti (2018): Thí nghiệm quy mô thực trên đoạn kênh thoát nước dài 300 m để ước tính hệ số nhám của thảm thực vật linh hoạt.
  • Lai và cộng sự (2014): Giải hệ phương trình Saint-Venant có xét đến trường trọng lực và lực Coriolis.
  1. Phương pháp sai phân hữu hạn (FDM) và mô hình hóa trong nước:
  • Crank-Nicholson: Sơ đồ sai phân hiện trung tâm trên 3 điểm không gian ở bước thời gian $j-1$ và 1 điểm ở bước $j$; mô hình DAMBRK áp dụng với trọng số $\theta = 0{,}6$.
  • Mac Cormack (1969): Phương pháp sai phân hiện khai triển chuỗi Taylor theo hai bước dự báo và hiệu chỉnh, đạt độ chính xác bậc 2 theo thời gian.
  • Preissmann (1961): Sơ đồ sai phân ẩn 4 điểm kinh điển trên hai lớp thời gian, được phân tích và ứng dụng bởi Abbott, Amein, Fang, Cunge.
  • Abbott và Ionescu (1967): Sơ đồ sai phân xen kẽ 8 điểm giải đồng thời độ sâu và lưu lượng.
  • Liggett và Woolhiser (1967): So sánh sơ đồ hiện (Leap-frog, Lax-Wendroff) với sơ đồ ẩn và phương pháp đường đặc trưng.
  • Vreugdenhil (1973): Sơ đồ ẩn tính xen kẽ mực nước và lưu lượng trên mắt lưới chữ nhật.
  • Francesco Greco và Lorenzo Panattoni (1974): Phương pháp sai phân ẩn giải phương trình Saint-Venant ứng dụng cho sông Arno.
  • Federici (1979), Ireneusz Stepien (1983 - ứng dụng sông Vistula), Chalfen và Niemiec (1986 - ứng dụng sông Nysa), Glaister (2005): Đánh giá tính ổn định và bảo toàn của các sơ đồ sai phân.
  • Nguyễn Như Khuê (cuối thập niên 1970): Xây dựng mô hình toán thủy lực VRSAP mô phỏng hệ thống kênh sông và ô đồng bằng.
  • Nguyễn Ân Niên: Phát triển phần mềm KOD1 dựa trên sơ đồ sai phân hiện.
  • Nguyễn Tất Đắc (2005): Nghiên cứu các sơ đồ sai phân, thuật toán truy đuổi, phát triển hệ thống phần mềm SAL (các phiên bản FWQ87, SAL790, TLUC, SAL1193, SAL99, SALBOD, VRSAP_SAL) giải bài toán dòng chảy và xâm nhập mặn.
  • Trần Đình Hợi, Lê Văn Nghị (2002): Giải sai phân phương trình dòng không ổn định để tính toán khả năng lấy nước hệ thống sông.
  • Vũ Đức Thái (2011): Ứng dụng mạng nơ-ron tế bào CNN trên chip FPGA EP2C35 để giải hệ phương trình đạo hàm riêng của dòng chảy một chiều.
  1. Phương pháp thể tích hữu hạn (FVM) và phương pháp đường đặc trưng:
  • Bojan Crnkovic và cộng sự (2009): Phát triển sơ đồ FVM nửa ẩn cân bằng có xét đến điều kiện CFL cục bộ cho phương trình Saint-Venant.
  • Dolz và cộng sự (2012): Nghiên cứu tích hợp mô hình FVM 1D và 2D trong tính toán lũ kênh tự nhiên.
  • Phương pháp đặc trưng: J. Massau khởi xướng; Amein, Fletcher, Liggett, Lai phát triển lưới đặc trưng cố định; Yao Hsin Hwang (2013) nghiên cứu phương pháp phân mảnh đặc trưng cho phương trình Saint-Venant.
  1. Phương pháp phần tử hữu hạn (FEM) trong cơ học chất lưu:
  • Alexander Hrennikoff (1941), Richard Courant (1942), Strang (1973): Thiết lập nền tảng toán học cho FEM.
  • Thomas J. Hughes và cộng sự (1984): Đề xuất phương pháp Taylor-Galerkin giải bài toán bình lưu - khuếch tán.
  • Quartapelle (1985, 1997, 1998): Phát triển sơ đồ Taylor-Galerkin bậc 2 và bậc 3 giải phương trình hyperbolic phi tuyến, phương trình Burgers nhớt và sóng nước phân tán trên lưới phi cấu trúc.
  • Nwaogazie (1985, 1986, 1987): Xây dựng chương trình Fortran WICFEM kết hợp Galerkin và lặp Newton-Raphson cho hệ Saint-Venant 1D và 2D.
  • Selmin (1987), Litsa Anastasiadou Partheniou và Terzidis (1988), Szymkiewicz (1991), Blandford và Ormsbee (1993), Tabarrok và Su (1993), Yi Zhang (2005), Chun Bin Liu (2005), Phongthanapanich và cộng sự (2006), Zhou Yi-Lin và cộng sự (2006), Li (2007), Zarmehi và cộng sự (2011), Xu Sun và cộng sự (2012): Mở rộng ứng dụng FEM cho các bài toán thủy lực và dòng chảy mặt tự do.
  • Nguyễn Thế Hùng (1989, 2001), Lê Văn Nghị, Trần Đình Hợi, Nguyễn Thế Hùng (2004): Ứng dụng Galerkin và sai phân giải hệ phương trình dòng chảy hai chiều đứng có thành phần giả nén và hệ phương trình Reynolds hai chiều đứng.

Khoảng trống luận án lựa chọn giải quyết: Mặc dù các phương pháp số và lý thuyết dòng chảy 1D đã phát triển phong phú, các nghiên cứu trước đây chưa thiết lập hệ phương trình một chiều tổng quát có tích hợp tường minh thành phần vận tốc thẳng đứng tại đáy lòng dẫn ($w^*$). Luận án chọn lấp đầy khoảng trống này bằng cách xây dựng hệ phương trình 1D suy rộng từ phương trình 2DV, giải bằng phương pháp Taylor-Galerkin bậc 3 trên ngôn ngữ Fortran 90 và kiểm nghiệm bằng thực nghiệm mô hình vật lý.

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

Luận án vận dụng các lý thuyết, khái niệm và công cụ phương pháp luận sau:

  1. Cơ sở lý thuyết và khung phân tích:
  • Mô hình rối chiều dài xáo trộn của Prandtl: Được lựa chọn nhằm đơn giản hóa quá trình diễn toán nhưng vẫn đảm bảo khả năng mô tả dòng chảy thực tế. Ứng suất tiếp rối được xác định theo hệ thức: $$\tau_t = \rho l^2 \left|\frac{du}{dz}\right| \frac{du}{dz}$$ Hệ số nhớt rối động lực $\eta_t = \rho l^2 \left|\frac{du}{dz}\right|$ và hệ số nhớt rối động học $\nu_t = l^2 \left|\frac{du}{dz}\right|$. Chiều dài xáo trộn $l$ phụ thuộc vào độ sâu tính từ đáy theo quan hệ $l = kz$, trong đó hằng số Von Karman $k = 0{,}4$ cho dòng chảy trong ống và $k = 0{,}54$ đối với sông thiên nhiên.
  • Ứng suất tiếp tổng cộng trong dòng chảy rối bao gồm ma sát nhớt phân tử và ma sát rối: $$\tau = \mu \frac{du}{dz} + \eta_t \frac{du}{dz}$$
  • Hệ phương trình xuất phát hai chiều đứng (2DV): Gồm phương trình liên tục và phương trình động lượng Navier-Stokes theo phương trục dọc $x$ và phương thẳng đứng $z$: Hệ trục tọa độ đặt tại đáy ở thượng lưu, trục $x$ dương hướng về hạ lưu, trục $z$ dương hướng lên trên.
  • Điều kiện biên vật lý:
    • Tại mặt thoáng ($z = h$): Áp suất dư bằng 0; điều kiện động học mặt thoáng $\frac{dh}{dt} = w_m$.
    • Tại đáy lòng dẫn ($z = 0$): Thành phần vận tốc dọc đáy $u_0 = 0$ do tính nhớt dính bám của chất lỏng thực; thành phần vận tốc đứng tại vị trí có dòng chảy phun lên là $w = w^*$.
    • Điều kiện biến đổi chậm: Giả thiết $w^*$ và chiều sâu $h$ biến đổi chậm theo thời gian, đạo hàm $\frac{\partial w}{\partial z} > \frac{\partial w}{\partial x}$.
  • Tích phân theo phương đứng: Sử dụng quy tắc tích phân Leibniz và phương pháp tích phân từng phần để chuyển đổi trường áp suất phi thủy tĩnh $p$ và vận tốc $u, w$ thành các đại lượng trung bình mặt cắt ngang trong hệ phương trình 1D.
  1. Phương pháp thực tế:
  • Phương pháp lý thuyết giải tích và biến đổi toán học: Thiết lập hệ phương trình 1D suy rộng thông qua tích phân giải tích dọc theo chiều sâu; phân tích bậc (order of magnitude analysis) của từng số hạng trong phương trình vi phân để loại bỏ các đại lượng vi phân vô cùng bé bậc cao.
  • Phương pháp số: Lựa chọn phương pháp phần tử hữu hạn Taylor-Galerkin. Phương pháp này thực hiện rời rạc hóa thời gian bằng khai triển chuỗi Taylor đạt độ chính xác bậc 3 trước khi rời rạc hóa không gian bằng phương pháp Galerkin tiêu chuẩn.
  • Phương pháp lập trình: Lập trình thuật toán Taylor-Galerkin bằng ngôn ngữ Fortran 90, xây dựng phần mềm chuyên dụng mang tên TG1D.
  • Phương pháp thực nghiệm: Thiết lập mô hình vật lý trong máng thủy lực mô phỏng dòng chảy hở một chiều có vận tốc thẳng đứng tại đáy lòng dẫn để thu thập chuỗi số liệu đối chứng.
  1. Nguồn tư liệu và dữ liệu:
  • Luận án sử dụng dữ liệu thực nghiệm đo đạc phân bố vận tốc của Xingwei Chen và Yee Meng Chiew (2004) tại các mặt cắt $x = -10\text{ cm}$, $25\text{ cm}$, $50\text{ cm}$ và $75\text{ cm}$ (được công bố trong bài báo về kênh hở có thấm đáy) để làm căn cứ tính toán phân tích bậc, tính tích phân hình thang và kiểm định mức độ ảnh hưởng của các số hạng vi phân.
  • Dữ liệu thực nghiệm đo đạc mực nước và vận tốc trên mô hình vật lý máng thủy lực do tác giả trực tiếp thực hiện.

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

Chương 1: Tổng quan về dòng chảy hở một chiều và phương pháp giải số

Chương 1 hệ thống hóa cơ sở lý luận của dòng chảy hở một chiều và đánh giá so sánh toàn diện các phương pháp giải số phương trình Saint-Venant cổ điển.

  1. Luận điểm về hệ phương trình Saint-Venant cổ điển và các giả thiết: Hệ hai phương trình vi phân Saint-Venant mô tả chuyển động không ổn định thay đổi chậm: $$\frac{\partial A}{\partial t} + \frac{\partial Q}{\partial x} = q$$ $$\frac{\partial Q}{\partial t} + \frac{\partial}{\partial x}\left(\beta \frac{Q^2}{A}\right) + gA\left(\frac{\partial h}{\partial x} + S_f - S_0\right) = 0$$ Các giả thiết cơ bản gồm: dòng chảy một chiều theo vận tốc trung bình mặt cắt; đường dòng có độ cong nhỏ; gia tốc thẳng đứng không đáng kể nên áp suất phân bố theo quy luật thủy tĩnh; chiều sâu biến đổi chậm; độ dốc đáy nhỏ ($\cos\theta \approx 1$); ma sát đáy tuân theo công thức dòng ổn định.

  2. Đánh giá ưu nhược điểm của các phương pháp số:

Phương pháp giải số Bậc chính xác Ưu điểm chính Nhược điểm / Hạn chế
Sai phân hữu hạn (FDM) Bậc 1 đến bậc 2 Thuật toán đơn giản, dễ lập trình, thời gian tính nhanh. Độ chính xác giới hạn ở bậc 2; nghiệm số không đảm bảo bảo toàn nghiêm ngặt; lưới không linh hoạt cho biên phức tạp.
Thể tích hữu hạn (FVM) Bậc 1 đến bậc 2 Bảo toàn khối lượng và thông lượng tốt; thuật toán tương đối đơn giản; lưới tương đối mềm dẻo. Chỉ áp dụng cho chất lưu; rất khó giải các phương trình vi phân bậc cao (bậc $\ge 3$); miền tính khó xử lý đồng thời nhiều vùng có tính chất vật lý khác nhau.
Đường đặc trưng (MOC) Tùy thuộc lược đồ Làm rõ bản chất vật lý truyền sóng; trực quan trong phân biệt vùng chảy êm và chảy xiết. Bước thời gian bị khống chế ngặt nghèo bởi điều kiện Courant-Friedrichs-Lewy; việc nội suy trên lưới biến đổi làm giảm độ chính xác.
Phần tử hữu hạn (FEM / Taylor-Galerkin) Bậc cao (đạt bậc $\ge 3$) Sai số nhỏ, độ chính xác cao; lưới mềm dẻo bám sát biên hình học; giải được phương trình vi phân bậc cao; đưa được điều kiện biên tự nhiên vào bài toán. Thuật toán phức tạp; khối lượng tính toán lớn; đòi hỏi kỹ thuật lập trình nâng cao.
  1. Tồn tại của các công cụ phần mềm hiện nay:
  • Các phần mềm 1D chuyên dụng (HEC-RAS, MIKE11, KOD1, SAL) chỉ tiếp nhận nguồn bổ sung ngang $q$ dạng dòng chảy bên với lưu tốc nhỏ, không thể mô phỏng thành phần vận tốc thẳng đứng $w^*$ ở đáy.
  • Các phần mềm như ANSYS 2D/3D giải quyết được nhưng phải chuyển sang bài toán nhiều chiều và áp dụng FVM có bậc chính xác thấp hơn FEM.
  • Do đó, chương 1 đi đến kết luận: Cần thiết phải phát triển hệ phương trình 1D suy rộng có kể đến vận tốc đứng tại đáy và sử dụng phương pháp phần tử hữu hạn Taylor-Galerkin bậc 3 trên nền tảng ngôn ngữ Fortran 90.

Chương 2: Mô hình toán dòng chảy hở một chiều suy rộng khi có kể đến vận tốc chiều đứng ở đáy lòng dẫn

Chương 2 trình bày chi tiết quá trình thiết lập toán học hệ phương trình vi phân dòng chảy hở một chiều suy rộng khi có kể đến vận tốc chiều đứng ở đáy lòng dẫn ($w^*$).

  1. Thiết lập biểu thức vận tốc đứng và áp suất phi thủy tĩnh:
  • Từ phương trình liên tục 2DV, tích phân từ đáy ($z = 0$) đến mặt thoáng ($z = h$) kết hợp quy tắc Leibniz, tác giả xác định biểu thức vận tốc đứng tại mặt thoáng $w_m$: $$w_m = -\int_{0}^{h} \frac{\partial u}{\partial x} dz + w^* = -\frac{\partial}{\partial x}\int_{0}^{h} u dz + u_m \frac{\partial h}{\partial x} + w^*$$
  • Tích phân phương trình động lượng theo phương đứng từ cao độ $z$ bất kỳ đến mặt thoáng $h$, có xét đến điều kiện biên áp suất dư tại mặt thoáng bằng 0, nhận được biểu thức phân bố áp suất phi thủy tĩnh $p(z)$. Áp suất tại mỗi điểm không chỉ phụ thuộc vào trọng lượng cột nước mà còn chứa các số hạng đạo hàm của vận tốc đứng và vận tốc dọc.

  1. Thiết lập phương trình chuyển động suy rộng: Tác giả thế biểu thức áp suất phi thủy tĩnh vào phương trình động lượng theo phương $x$, tiến hành tích phân từng phần trên toàn bộ chiều sâu dòng chảy từ 0 đến $h$. Quá trình biến đổi sử dụng các điều kiện biên dính bám tại đáy ($u_0 = 0$), điều kiện động học mặt thoáng và đạo hàm hàm hợp.

  2. Phân tích bậc các số hạng bằng số liệu thực nghiệm: Để loại bỏ các số hạng vi phân bậc cao có giá trị vô cùng bé mà không làm mất tính chính xác của phương trình, tác giả sử dụng số liệu phân bố lưu tốc thực nghiệm từ bài báo của Xingwei Chen và Yee Meng Chiew (2004) tại 4 mặt cắt không gian:

Tọa độ mặt cắt $x$ (cm) Mục đích phân tích Phương pháp tính phân tích bậc
$x = -10\text{ cm}$ Xác định phân bố lưu tốc thực nghiệm $u(z)$ thượng lưu vùng đáy thấm Tích phân số trị theo công thức hình thang
$x = 25\text{ cm}$ Khảo sát biến thiên lưu tốc tại vùng bắt đầu chịu ảnh hưởng đáy Nội suy và tích phân diện tích hình thang
$x = 50\text{ cm}$ Đánh giá gradient vận tốc $\frac{\partial u}{\partial x}$ và $\frac{\partial u}{\partial z}$ tại tâm vùng ảnh hưởng So sánh độ lớn định lượng của các số hạng đạo hàm riêng
$x = 75\text{ cm}$ Kiểm chứng trạng thái phân bố vận tốc ở hạ lưu vùng ảnh hưởng Đánh giá độ tắt dần của các số hạng bậc cao

Dựa trên kết quả tích phân hình thang từ bảng số liệu thực nghiệm trên, tác giả chứng minh được rằng tích của các đạo hàm riêng bậc cao và các số hạng chứa $(\frac{\partial w}{\partial x})$ có giá trị rất nhỏ so với các thành phần chính và hoàn toàn có thể triệt tiêu. Kết quả dẫn đến hệ phương trình vi phân một chiều suy rộng hoàn chỉnh có xét đến thành phần vận tốc thẳng đứng tại đáy lòng dẫn $w^*$.

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

Các đóng góp mới của luận án được văn bản tự xác định bao gồm ba điểm cốt lõi:

  1. Đóng góp về lý luận toán học thủy lực: Luận án đã xây dựng thành công hệ phương trình một chiều suy rộng cho dòng chảy hở khi có xét đến thành phần vận tốc tương đối lớn theo phương thẳng đứng ở đáy lòng dẫn ($w^*$). Hệ phương trình này tổng quát hơn hệ phương trình vi phân Saint-Venant cổ điển năm 1871, khắc phục giả thiết phân bố áp suất thủy tĩnh và cho phép mô hình hóa các hiện tượng thủy lực phức tạp tại đáy lòng dẫn trong khuôn khổ mô hình 1D.

  2. Đóng góp về phương pháp giải số và công cụ tính toán: Luận án đã thiết lập thuật toán và xây dựng chương trình tính số trị bằng ngôn ngữ Fortran 90, đặt tên là TG1D. Chương trình ứng dụng phương pháp phần tử hữu hạn Taylor-Galerkin đạt bậc chính xác bậc 3 theo thời gian, cho phép giải hệ phương trình 1D suy rộng với độ ổn định cao và khả năng mô tả chính xác các biến đổi cục bộ của dòng chảy.

  3. Đóng góp về thực nghiệm và cấu trúc dòng chảy: Luận án đã thực hiện thí nghiệm trên mô hình vật lý máng thủy lực dưới điều kiện dòng chảy một chiều có vận tốc theo phương thẳng đứng ở đáy lòng dẫn. Chuỗi số liệu thực nghiệm này đã cung cấp cơ sở kiểm chứng độ tin cậy của thuật toán và phần mềm TG1D, đồng thời đóng góp dữ liệu thực nghiệm phục vụ việc phân tích cấu trúc thủy lực của dòng chảy hở khi có dòng phun hoặc nước trồi từ đáy.

Kiến nghị và giải pháp ứng dụng thực tiễn: Luận án kiến nghị sử dụng mô hình và phần mềm TG1D để tính toán, thiết kế cho các bài toán thực tế thuộc lĩnh vực tài nguyên nước và công trình thủy lợi, bao gồm:

  • Tính toán dòng chảy mùa kiệt trên các hệ thống kênh, sông đi qua khu vực địa chất có hiện tượng nước trồi hoặc dòng ngầm có áp phun lên từ đáy (như kênh vùng trung du, kênh chạy qua thung lũng, triền đồi, kênh dẫn qua lòng hồ Phú Ninh).
  • Tính toán nâng cao trình mặt nước trên kênh dẫn có bố trí đập ngầm nhân tạo nhằm phục vụ mục đích cấp nước tưới tự chảy.
  • Mô phỏng dòng chảy tràn qua thân đê biển tại các đoạn bờ biển thẳng.

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

Hạn chế của nghiên cứu:

  • Giả thiết mô hình toán bỏ qua ảnh hưởng của trường gió và lực Coriolis trên bề mặt dòng chảy.
  • Giả thiết độ cong của đường dòng là nhỏ và độ dốc đáy trung bình đủ nhỏ sao cho $\cos\theta \approx 1$.
  • Giả thiết bỏ qua dung tích phần trữ của mặt cắt ngang lòng dẫn ($S = 0$).
  • Giả thiết thành phần vận tốc đáy $w^*$ và chiều sâu $h$ biến đổi chậm theo thời gian, đạo hàm vận tốc đứng thỏa mãn điều kiện $\frac{\partial w}{\partial z} > \frac{\partial w}{\partial x}$.
  • Mô hình rối áp dụng dừng lại ở mức mô hình giải tích chiều dài xáo trộn của Prandtl với hằng số Von Karman kinh nghiệm ($k = 0{,}54$ cho sông tự nhiên và $k = 0{,}4$ cho ống), chưa áp dụng các mô hình rối hai phương trình vi phân bậc cao.

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

  • Mở rộng hệ phương trình suy rộng cho trường hợp dòng chảy xét đến ảnh hưởng của lực Coriolis và trường gió mặt thoáng.
  • Nghiên cứu tích hợp mô hình rối bậc cao hơn vào việc xác định phân bố ứng suất tiếp dọc theo chiều sâu lòng dẫn.
  • Bổ sung dung tích phần trữ mặt cắt ngang và áp dụng mô hình cho các đoạn sông có hình học mặt cắt tự nhiên phức tạp.

Giá trị tham khảo

Luận án mang lại giá trị tham khảo cụ thể cho các nhóm đối tượng học thuật và kỹ thuật:

  1. Đối tượng thụ hưởng:
  • Nghiên cứu sinh, học viên cao học và giảng viên chuyên ngành Kỹ thuật tài nguyên nước, Thủy lực và Thủy văn, Xây dựng công trình thủy lợi, Cơ học chất lưu.
  • Các kỹ sư thiết kế công trình thủy lợi, chuyên gia tư vấn quy hoạch nguồn nước và cán bộ kỹ thuật tại các cơ quan quản lý đê điều, khai thác công trình thủy nông.
  1. Các nội dung đáng tham khảo:
  • Quy trình tích phân giải tích chuyển đổi từ hệ phương trình vi phân chuyển động hai chiều đứng (2DV) sang hệ phương trình một chiều suy rộng có áp suất phi thủy tĩnh và vận tốc đáy.
  • Kỹ thuật phân tích bậc dựa trên số liệu thực nghiệm phân bố lưu tốc để tinh giản hệ phương trình vi phân.
  • Thuật toán số và cấu trúc chương trình nguồn Fortran 90 của mô hình phần tử hữu hạn Taylor-Galerkin bậc 3 (chương trình TG1D).
  • Chuỗi số liệu đo đạc thực nghiệm trên mô hình vật lý máng thủy lực phục vụ kiểm chứng các mô hình số trị dòng chảy hở.

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

  1. Điểm khác biệt cơ bản giữa hệ phương trình một chiều của luận án và hệ phương trình Saint-Venant (1871) cổ điển là gì? Hệ phương trình Saint-Venant cổ điển giả thiết gia tốc thẳng đứng không đáng kể, áp suất phân bố theo quy luật thủy tĩnh và bỏ qua vận tốc đứng tại đáy. Ngược lại, hệ phương trình của luận án là mô hình 1D suy rộng, tích hợp tường minh thành phần vận tốc tương đối lớn theo phương thẳng đứng ở đáy lòng dẫn ($w^*$) và xét đến sự phân bố áp suất phi thủy tĩnh dọc theo chiều sâu cột nước.

  2. Tại sao các phần mềm thương mại 1D hiện hành như HEC-RAS hay MIKE11 không giải quyết được bài toán luận án đặt ra? Các phần mềm thương mại 1D hiện nay chỉ cho phép đưa vào điều kiện biên nguồn bổ sung ngang ($q$ - lateral flows) với vận tốc nhỏ. Nguồn bổ sung ngang có bản chất cơ học khác với vận tốc đứng ở đáy ($w^*$), vì vận tốc thẳng đứng tại đáy trực tiếp cản trở dòng chuyển động chính dọc trục và gây hiện tượng dâng mực nước ở thượng lưu lớn hơn nhiều so với nhập lưu ngang.

  3. Thuật toán số nào được lựa chọn để giải hệ phương trình suy rộng và ưu điểm của thuật toán này là gì? Thuật toán được lựa chọn là phương pháp phần tử hữu hạn Taylor-Galerkin, lập trình bằng ngôn ngữ Fortran 90. Phương pháp này đạt độ chính xác bậc 3 theo thời gian, có sai số nhỏ, khả năng thích ứng lưới rất mềm dẻo bám sát biên bài toán và xử lý tốt các bài toán có hệ phương trình vi phân bậc cao mà phương pháp sai phân hoặc thể tích hữu hạn gặp hạn chế.

  4. Tác giả đã sử dụng dữ liệu nào để thực hiện phân tích bậc nhằm loại bỏ các số hạng bé trong phương trình vi phân? Tác giả sử dụng số liệu thực nghiệm đo đạc phân bố vận tốc dòng chảy rối theo chiều sâu của Xingwei Chen và Yee Meng Chiew (2004) tại các mặt cắt $x = -10\text{ cm}$, $25\text{ cm}$, $50\text{ cm}$ và $75\text{ cm}$. Phương pháp tích phân hình thang được áp dụng trên các dữ liệu này để tính toán độ lớn định lượng của các số hạng đạo hàm riêng, làm căn cứ loại bỏ các số hạng vô cùng bé.

  5. Hằng số Von Karman được thiết lập với giá trị bao nhiêu trong mô hình rối của luận án? Trong mô hình rối chiều dài xáo trộn $l = kz$ của Prandtl được tác giả áp dụng, hằng số Von Karman $k$ được lấy giá trị bằng $0{,}4$ đối với dòng chảy trong ống và bằng $0{,}54$ đối với dòng chảy trong sông thiên nhiên.

Kết luận

Luận án đã giải quyết bài toán dòng chảy hở một chiều khi có kể đến vận tốc tương đối lớn theo phương thẳng đứng ở đáy lòng dẫn thông qua việc thiết lập hệ phương trình vi phân 1D suy rộng từ hệ phương trình 2DV. Nghiên cứu đã ứng dụng phương pháp phần tử hữu hạn Taylor-Galerkin bậc 3 theo thời gian và lập trình hoàn chỉnh chương trình tính TG1D bằng ngôn ngữ Fortran 90. Độ tin cậy của mô hình toán và thuật toán số đã được kiểm chứng bằng thực nghiệm vật lý trong máng thủy lực, đồng thời cung cấp công cụ tính toán thủy lực cho các bài toán kỹ thuật thực tế như dòng chảy qua vùng nước trồi, đập ngầm và đê biển.