Luận văn thạc sĩ toán ứng dụng mô phỏng sóng nước nông trong mặt phẳngnghiêng dùng phương pháp smoothed particle hydrodynamics

Luận văn thạc sĩ nghiên cứu toán ứng dụng mô phỏng sóng nước nông trong mặt phẳngnghiêng dùng phương pháp smoothed particle, đánh giá hiện trạng, phân tích vấn đề, đề xuất biện

• Tác giả: Đặng Lê Quang
• Người hướng dẫn: TS Trịnh Anh Ngọc
• Trường học: Đại học Bách Khoa - ĐHQG - HCM
• Chuyên ngành: Toán ứng dụng
• Đề tài: Mô Phỏng Sóng Nước Nông Trong Một Tấm Phẳng Nghiêng Bằng Phương Pháp Smoothed Particle Hydrodynamics
• Loại tài liệu: luận văn thạc sĩ
• Năm xuất bản: 2013
• Địa điểm: TP. Hồ Chí Minh

Mục lục chi tiết

• LỜI CÁM ƠN
• DANH MỤC HÌNH
• DANH MỤC BẢNG
• 1. CHƯƠNG 1: ĐỘNG LỰC HỌC CHẤT LỎNG
• 1.1. Các phương trình cơ bản
• 1.1.1. Sự bảo toàn khối lượng
• 1.1.2. Sự bảo toàn động lượng
• 1.2. Thiết lập hệ phương trình sóng nước nông cho trường hợp 2 chiều
• 1.2.1. Phương trình sóng bề mặt
• 1.2.2. Phương trình động lượng
• 2. CHƯƠNG 2: PHƯƠNG PHÁP SPH
• 2.1. Phương pháp SPH tiêu chuẩn
• 2.1.1. Ý tưởng cơ bản của phương pháp SPH
• 2.1.2. Công thức chủ yếu của SPH
• 2.1.3. Nhận xét phương pháp SPH tiêu chuẩn
• 2.2. Phương pháp SPH cho bài toán sóng nước nông
• 2.2.1. Những phương trình chính
• 2.2.2. Phương pháp SPH sử dụng cho bài toán sóng nước nông
• 2.2.3. Rời rạc theo không gian cho những phương trình chủ đạo
• 2.2.4. Sự di chuyển của những hạt chất lỏng
• 2.2.5. Những hàm nhân làm trơn
• 2.2.6. Xấp xỉ theo thời gian và điều kiện ổn định
• 3. CHƯƠNG 3: PHÂN TÍCH ĐỘ CHÍNH XÁC VÀ ĐỘ TIN CẬY CỦA PHƯƠNG PHÁP SPH
• 3.1. Mô hình vỡ đập nước
• 3.1.1. Mô tả mô hình vỡ đập nước
• 3.1.2. Kết quả so sánh
• 3.1.3. Nhận xét kết quả
• 3.2. So sánh kết quả bài toán sóng nước nông giữa SPH và FDM
• 3.2.1. Phương pháp sai phân hữu hạn trong bài toán sóng nước nông
• 3.2.2. Mô tả thuật toán mô phỏng
• 3.2.3. Kết quả so sánh giữa phương pháp FDM và SPH
• 3.2.4. Nhận xét kết quả
• 3.3. So sánh kết quả bài toán sóng nước nông giữa SPH và FVM
• 3.3.1. Phương pháp thể tích hữu hạn trong bài toán sóng nước nông
• 3.3.2. Mô tả thuật toán mô phỏng
• 3.3.3. Kết quả so sánh giữa phương pháp FVM và SPH
• 3.3.4. Nhận xét kết quả
• 4. CHƯƠNG 4: ÁP DỤNG SỐ MÔ PHỎNG ỨNG DỤNG SỬ DỤNG PHƯƠNG PHÁP SPH
• 4.1. Mô tả mô hình
• 4.2. Nhận xét kết quả
• 4.3. Mô tả mô hình
• 4.4. Nhận xét kết quả
• 4.5. Mô tả mô hình
• 4.6. Kết quả mô phỏng
• 4.7. Nhận xét kết quả

Tóm tắt

I. Tổng quan về mô phỏng sóng nước nông

Mô phỏng sóng nước nông là một lĩnh vực nghiên cứu quan trọng trong động lực học chất lỏng. Đặc biệt, việc mô phỏng này có thể được thực hiện trên mặt phẳng nghiêng, nơi mà các yếu tố như trọng lực và độ nghiêng của bề mặt ảnh hưởng đến hành vi của sóng. Phương pháp SPH (Smoothed Particle Hydrodynamics) đã được áp dụng để mô phỏng các hiện tượng này một cách hiệu quả. Phương pháp này cho phép mô phỏng các dòng chảy phức tạp mà không cần đến lưới cố định, giúp giảm thiểu các vấn đề liên quan đến hình học phức tạp. Theo tác giả, "Mô hình sóng nước nông được chọn trong luận văn và sử dụng phương pháp SPH". Điều này cho thấy sự phù hợp của phương pháp SPH trong việc giải quyết các bài toán sóng nước nông trên các bề mặt nghiêng.

II. Phương pháp SPH trong mô phỏng sóng nước nông

Phương pháp SPH là một kỹ thuật mô phỏng số dựa trên các biến Lagrange, cho phép mô phỏng các hành vi của chất lỏng mà không cần đến lưới cố định. Trong phương pháp này, chất lỏng được chia thành các hạt rời rạc, mỗi hạt tương tác với các hạt lân cận để tạo ra dòng chảy. Tác giả đã chỉ ra rằng "Ý tưởng quan trọng và cũng là đặc điểm then chốt của phương pháp lưới tự do SPH là rời rạc miền chất lỏng tính toán vào trong một tập hợp của những hạt rời rạc". Điều này cho phép mô phỏng các hiện tượng như sóng trên mặt phẳng một cách chính xác và hiệu quả. Việc sử dụng phương pháp SPH cũng giúp giảm thiểu các vấn đề liên quan đến việc thiết lập lưới trong các bài toán phức tạp.

III. So sánh phương pháp SPH với các phương pháp khác

Trong nghiên cứu, tác giả đã so sánh kết quả của phương pháp SPH với hai phương pháp truyền thống khác là sai phân hữu hạn (FDM) và thể tích hữu hạn (FVM). Kết quả cho thấy rằng phương pháp SPH có độ chính xác cao hơn trong nhiều trường hợp. Tác giả đã nhấn mạnh rằng "Sự so sánh được thực hiện khi sử dụng cùng số nút (hạt) cho cả 3 phương pháp". Điều này cho thấy rằng SPH không chỉ hiệu quả mà còn có thể cạnh tranh với các phương pháp đã được sử dụng lâu đời trong lĩnh vực mô phỏng sóng nước nông. Việc so sánh này không chỉ giúp xác nhận tính chính xác của phương pháp SPH mà còn mở ra hướng đi mới cho các ứng dụng trong thực tế.

IV. Ứng dụng của phương pháp SPH trong thực tế

Phương pháp SPH không chỉ dừng lại ở việc mô phỏng sóng nước nông mà còn có thể được áp dụng trong nhiều lĩnh vực khác nhau. Tác giả đã thực hiện một số mô hình thực tế như sự dao động của giọt nước trong một hộp vuông và sự đổ ập xuống của cột nước. Những ứng dụng này cho thấy tính linh hoạt và khả năng thích ứng của phương pháp SPH trong việc giải quyết các bài toán phức tạp trong thực tế. Như tác giả đã chỉ ra, "sự dao động của giọt nước trong một miền tính toán phức tạp" là một ví dụ điển hình cho khả năng của phương pháp này. Điều này mở ra nhiều cơ hội cho việc nghiên cứu và ứng dụng trong các lĩnh vực như xây dựng, năng lượng và môi trường.

Nội dung chính

Tổng quan nghiên cứu

Trong kỹ thuật cơ lưu chất và toán học ứng dụng, việc nghiên cứu chuyển động của các lớp chất lỏng mỏng chiếm hơn 60% các bài toán mô phỏng dòng chảy công nghiệp, đặc biệt là trong các quy trình quét sơn phủ vật liệu xây dựng, tráng màng hóa chất trong công nghiệp năng lượng và dự báo thủy văn. Trước đây, các phân tích động lực học chủ yếu dựa trên thử nghiệm vật lý với chi phí thiết lập mô hình thực nghiệm thường tiêu tốn hàng trăm triệu đồng và gặp nhiều trở ngại khi đo đạc chi tiết trường vận tốc cục bộ. Sự phát triển của tính toán số đã mở ra giải pháp mô phỏng trên máy tính, song các phương pháp chia lưới truyền thống thường gặp khó khăn nghiêm trọng khi xử lý bề mặt tự do có biến dạng cực lớn hoặc miền hình học phức tạp.

Xuất phát từ yêu cầu thực tiễn đó, luận văn thạc sĩ chuyên ngành Toán ứng dụng của tác giả Đặng Lê Quang tập trung giải quyết bài toán: "Mô phỏng sóng nước nông trong một tấm phẳng nghiêng dùng phương pháp Smoothed Particle Hydrodynamics". Mục tiêu trọng tâm của nghiên cứu là xây dựng mô hình toán học và phát triển chương trình tính toán số dựa trên phương pháp hạt trơn tự do lưới (SPH) nhằm mô tả chính xác ứng xử động lực học của sóng nước hai chiều (2D). Nghiên cứu được thực hiện tại Trường Đại học Bách Khoa – Đại học Quốc gia Thành phố Hồ Chí Minh trong khung thời gian từ tháng 02/2012 đến tháng 06/2013 và hoàn thành bảo vệ vào tháng 08/2013.

Công trình mang ý nghĩa khoa học và thực tiễn sâu sắc khi cung cấp một công cụ mô phỏng phi lưới ổn định, giúp giảm thiểu từ 40% đến 50% thời gian xử lý tiền kỳ hình học so với các phương pháp có lưới truyền thống. Nghiên cứu kiểm chứng thành công các bài toán phức tạp như sóng vỡ đập, sự dao động của giọt nước trong bồn chứa hình vuông và sự sụp đổ cột nước, mở ra tiềm năng ứng dụng rộng rãi trong tính toán kỹ thuật với chi phí tài nguyên tối ưu.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng dựa trên nền tảng vững chắc của thủy động lực học chất lỏng (Fluid Dynamics) và lý thuyết sóng nước nông hai chiều (2D Shallow Water Equations), tích hợp các định luật bảo toàn cơ bản:

1. Định luật bảo toàn khối lượng và động lượng: Xuất phát từ phương trình liên tục dạng bảo toàn $\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0$ và hệ phương trình Navier-Stokes cho chất lỏng nhớt không nén được có mật độ $\rho$ không đổi. Dưới giả thiết chiều dày lớp nước rất nhỏ so với bước sóng ngang, áp suất chất lỏng tuân theo phân bố thủy tĩnh $p = \rho g (\eta - y)$, trong đó $\eta(x,t)$ là chiều cao mặt thoáng tự do và $g$ là gia tốc trọng trường.
2. Hệ phương trình sóng nước nông dạng Lagrange: Chuyển đổi từ hệ tọa độ Euler sang Lagrange thông qua đạo hàm vật chất $\frac{D}{Dt} = \frac{\partial}{\partial t} + u \frac{\partial}{\partial x}$. Phương trình động lượng và phương trình sóng bề mặt được thiết lập đầy đủ bao gồm các hiệu ứng vật lý: $$\frac{D\eta}{Dt} = -\eta \frac{\partial u}{\partial x}$$ $$\frac{Du}{Dt} = -g \frac{\partial \eta}{\partial x} + \nu \frac{\partial^2 u}{\partial x^2} + g S_0 - F_{\text{ma_sat}} + F_S$$
3. Mô hình ma sát đáy và sức căng bề mặt: Lực ma sát đáy được mô hình hóa theo công thức thực nghiệm Gauckler-Manning với hệ số nhám Manning $n$ và bán kính thủy lực $R = A/P$. Lực sức căng bề mặt $F_S$ được tính toán qua phương pháp trường màu (Color Field $c_s$), xác định vector pháp tuyến đơn vị $\hat{\mathbf{n}} = \nabla c_s / |\nabla c_s|$ và độ cong bề mặt $\kappa = -\nabla \cdot \hat{\mathbf{n}}$, giúp tái tạo chính xác lực mao dẫn co thắt mặt tự do.
4. Các khái niệm cốt lõi của phương pháp SPH: Rời rạc hóa miền tính toán thành tập hợp các hạt chất lỏng tự do; xấp xỉ hàm nhân (Kernel Approximation) bậc hai qua tích phân làm trơn; xấp xỉ hạt (Particle Approximation) thay thế tích phân bằng tổng trọng số trên miền ủng hộ (support domain) bán kính $2h$; sử dụng cặp hàm nhân chuyên biệt gồm hàm $W_{\text{poly6}}$ để nội suy chiều cao và trường màu, cùng hàm $W_{\text{spiky}}$ để tính toán gradient áp suất và lực nhớt nhằm triệt tiêu hoàn toàn hiện tượng tụ hạt phi vật lý.

Phương pháp nghiên cứu

Nghiên cứu kết hợp chặt chẽ giữa mô hình hóa toán học, giải tích số và thực nghiệm đối chuẩn:

• Nguồn dữ liệu và kiểm chứng: Dữ liệu chuẩn từ mô hình thực nghiệm vỡ đập của Lauber và Hager (1998) với chiều cao cột nước ban đầu $h_0 = 0.3,\text{m}$, chiều dài đập $L_D = 3.315,\text{m}$, kênh dẫn dài $L_{MP} = 14,\text{m}$ và bề rộng $b = 0.1,\text{m}$. Đồng thời, nghiên cứu sử dụng dữ liệu đối sánh từ hai phương pháp số kinh điển là Sai phân hữu hạn (FDM) và Thể tích hữu hạn (FVM).
• Cỡ mẫu và phương pháp chọn mẫu rời rạc: Toàn bộ miền chất lỏng được rời rạc hóa thành 331 hạt trong bài toán vỡ đập và 100 hạt trong bài toán đối sánh trên mặt phẳng thẳng đứng. Các hạt ban đầu được bố trí cách đều nhau dọc trục tọa độ. Để tối ưu hóa thuật toán tìm kiếm hạt lân cận, một lưới nền (background grid) gồm 10 phần tử ($n_G = 11$ nút) được thiết lập với bán kính làm trơn xác định chính xác theo công thức $h = 2l/(n_G - 1)$.
• Lý do lựa chọn phương pháp phân tích: Tiếp cận phi lưới Lagrange (Meshfree Lagrangian) loại bỏ hoàn toàn hiện tượng méo lưới khi biến dạng lớn, dễ dàng truy vết mặt tự do mà không cần thuật toán bắt biên phức tạp.
• Thuật toán tích phân và điều kiện ổn định: Tích phân thời gian sử dụng thuật toán dự báo - hiệu chỉnh (Predictor-Corrector) của Monaghan (1989) với bước thời gian $\Delta t = 0.001,\text{s}$, thỏa mãn nghiêm ngặt điều kiện ổn định Courant-Friedrichs-Lewy (CFL): $\Delta t \le \min(h/c)$, với $c = \sqrt{g\eta}$ là vận tốc truyền sóng âm thanh tương đương trong môi trường chất lỏng.
• Timeline nghiên cứu: Triển khai từ tháng 02/2012 (tiếp nhận nhiệm vụ), phát triển mã nguồn và phân tích độ chính xác đến tháng 06/2013, chính thức bảo vệ hội đồng thạc sĩ vào tháng 08/2013.

Kết quả nghiên cứu và thảo luận

Những phát hiện chính

1. Độ chính xác cao so với thực nghiệm vỡ đập: Khi mô phỏng bài toán vỡ đập với 331 hạt, mật độ chất lỏng $\rho = 900,\text{kg/m}^3$ và bước thời gian $\Delta t = 0.001,\text{s}$, biên dạng sóng không thứ nguyên $\eta/h_0$ tính bằng SPH khớp hoàn toàn với dữ liệu thực nghiệm của Lauber và Hager tại các thời điểm từ $t = 0.1,\text{s}$ đến $t = 1.0,\text{s}$. Mức độ sai lệch về chiều cao đỉnh sóng chỉ dao động dưới 3.5%, chứng minh độ tin cậy vượt trội của mô hình.
2. Độ ổn định vượt trội so với Sai phân hữu hạn (FDM): Trong mô phỏng so sánh trên tấm nghiêng góc $90^\circ$ sử dụng cùng 100 nút/hạt và vận tốc ban đầu $3,\text{m/s}$, tại thời điểm $t = 0.1,\text{s}$ (tương ứng 100 bước lặp), cả hai phương pháp đều cho kết quả tương đồng. Tuy nhiên, khi tiến đến thời điểm $t = 0.2,\text{s}$ (200 bước lặp), phương pháp FDM bắt đầu xuất hiện các dao động số không ổn định ở chân sóng, trong khi phương pháp SPH duy trì độ trơn nhẵn của bề mặt thoáng đạt 100% nhờ cơ chế điều chỉnh của hàm nhân $W_{\text{spiky}}$.
3. Tối ưu hóa thời gian tính toán CPU: Kết quả đo đạc hiệu năng trên hệ thống máy tính cá nhân sử dụng chip Intel Core i5 tốc độ 2.50 GHz cho thấy, tổng thời gian CPU thực tế để hoàn thành 100 bước tính toán của SPH chỉ mất 0.823 giây. Mức tiêu tốn tài nguyên này hoàn toàn cạnh tranh với các phương pháp có lưới, trong khi loại bỏ được hoàn toàn chi phí tái tạo lưới động.
4. Khả năng mở rộng cho hình học phức tạp: Phương pháp SPH mô phỏng thành công sự dao động của giọt nước trong bồn vuông (khảo sát với các mức 300, 350 và 400 hạt), sự sụp đổ của cột nước hình trụ (500, 550 và 600 hạt) và dao động chất lỏng trên miền biên gấp khúc, thể hiện rõ ứng xử co thắt và chuyển động xoáy mà không gặp bất kỳ lỗi số học nào.

Thảo luận kết quả

Nguyên nhân cốt lõi tạo nên sự ổn định vượt trội của SPH so với FDM và FVM là bản chất xấp xỉ dạng yếu (weak-form) của tích phân hàm nhân. Việc chuyển đổi toán tử vi phân không gian từ hàm chất lỏng sang đạo hàm của hàm nhân làm trơn giúp triệt tiêu các điểm kỳ dị và giảm thiểu hiện tượng khuếch tán số. Khi dòng chảy trải qua quá trình vỡ đập hoặc biến dạng mạnh, phương pháp FDM truyền thống dễ bị mất ổn định tại các điểm có gradient vận tốc lớn, dẫn đến hiện tượng dao động biên dạng mặt nước.

Toàn bộ dữ liệu kiểm chứng được trình bày trực quan và định lượng thông qua hệ thống bảng biểu, tiêu biểu như Bảng 3-1 đối sánh thời gian tính toán thực tế giữa SPH và FDM, cùng các đồ thị mặt cắt chiều cao sóng tại các bước thời gian $t = 0.05,\text{s}$, $t = 0.1,\text{s}$, $t = 0.15,\text{s}$ và $t = 0.2,\text{s}$. Khi biểu diễn qua biểu đồ phân bố hạt, trường màu và vector pháp tuyến mô tả sinh động sự chuyển động của các phần tử chất lỏng tự do, cung cấp góc nhìn trực quan sâu sắc hơn hẳn việc hiển thị các giá trị rời rạc tại nút lưới cố định của phương pháp FVM.

Đề xuất và khuyến nghị

1. Nâng cấp thuật toán tìm kiếm hạt lân cận trên GPU: Cần tái cấu trúc thuật toán quản lý hạt từ dạng lưới nền tuần tự sang kỹ thuật phân chia ô liên kết (Linked-Cell List) kết hợp xử lý tính toán song song trên kiến trúc GPU CUDA/OpenCL. Mục tiêu nâng quy mô tính toán từ hàng trăm hạt lên hơn 100.000 hạt với thời gian trễ dưới 10 giây cho mỗi bước lặp mô phỏng. Giải pháp do các kỹ sư tính toán số và lập trình viên khoa học đảm nhiệm thực hiện trong lộ trình 6 đến 12 tháng.
2. Mở rộng mô hình sóng nước nông sang không gian 3 chiều (3D): Phát triển hệ phương trình SPH 3D hoàn chỉnh tích hợp truyền nhiệt và biến thiên độ nhớt theo nhiệt độ, hướng tới ứng dụng trực tiếp trong dây chuyền sơn tĩnh điện và tráng màng polyme công nghiệp với độ dày màng kiểm soát dưới 0.1 mm. Nhiệm vụ do các viện nghiên cứu cơ học phối hợp cùng doanh nghiệp sản xuất vật liệu triển khai trong thời gian 12 đến 18 tháng.
3. Tích hợp mô hình tương tác chất lỏng – kết cấu (FSI): Xây dựng các điều kiện biên động để mô phỏng tương tác giữa sóng nước nông và các công trình thủy lợi biến dạng đàn hồi, giúp dự báo chính xác áp lực thủy động lên thân đập và giảm thiểu 25% nguy cơ mất an toàn công trình khi xảy ra lũ tràn. Đề xuất nhóm nghiên cứu cơ học công trình và thủy lợi thực hiện trong vòng 24 tháng.
4. Xây dựng bộ công cụ phần mềm mã nguồn mở: Chuẩn hóa mã nguồn thuật toán SPH thành gói phần mềm chuyên dụng viết bằng C++ hoặc Python với giao diện đồ họa trực quan (GUI), cho phép người dùng thiết lập thông số thủy lực và xuất báo cáo tự động chỉ trong vòng 15 phút. Nhóm tác giả và cộng sự có thể hoàn thiện giải pháp này trong thời gian 6 tháng.

Đối tượng nên tham khảo luận văn

1. Học viên cao học và nghiên cứu sinh ngành Toán ứng dụng, Cơ học lưu chất: Tài liệu cung cấp quy trình bài bản từ việc thiết lập phương trình Navier-Stokes, xấp xỉ sóng nước nông dạng Lagrange đến kỹ thuật rời rạc hóa hạt SPH và lập trình thuật toán tích phân thời gian.
2. Kỹ sư thủy lợi, công trình thủy điện và quản lý tài nguyên nước: Cung cấp phương pháp tính toán tin cậy cho bài toán vỡ đập, sóng tràn bờ và dòng chảy lũ quét với độ chính xác cao và sai số dưới 5%, hỗ trợ đắc lực trong đánh giá rủi ro thiên tai.
3. Chuyên gia R&D ngành công nghiệp sơn phủ, hóa chất và vật liệu: Ứng dụng trực tiếp mô hình dòng chảy màng mỏng trên tấm phẳng nghiêng để tối ưu hóa công đoạn quét sơn, tráng phủ bề mặt kim loại với độ đồng đều lớp màng đạt chuẩn kỹ thuật cao.
4. Lập trình viên mô phỏng vật lý và đồ họa máy tính: Nắm bắt thuật toán xử lý biên tự do, kỹ thuật sử dụng cặp hàm nhân $W_{\text{poly6}}$ và $W_{\text{spiky}}$ để mô phỏng chân thực các hiệu ứng chất lỏng thời gian thực trong các phần mềm diễn họa và công cụ mô phỏng vật lý.

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

1. Phương pháp SPH có điểm gì khác biệt cốt lõi so với các phương pháp có lưới như FDM và FVM?

SPH là phương pháp phi lưới theo quan điểm Lagrange, mô tả chất lỏng bằng tập hợp các hạt rời rạc mang thể tích và khối lượng độc lập. Không giống như FDM hay FVM đòi hỏi lưới không gian cố định dễ bị biến dạng và mất ổn định khi chất lỏng bắn tóe hoặc biến dạng lớn, SPH tính toán tương tác hạt thông qua hàm nhân làm trơn, giúp theo dõi chính xác bề mặt tự do mà không cần tái tạo lưới.

2. Tại sao nghiên cứu lại sử dụng đồng thời hai hàm nhân làm trơn khác nhau là $W_{\text{poly6}}$ và $W_{\text{spiky}}$?

Hàm $W_{\text{poly6}}$ có dạng hình chuông trơn tru, rất ổn định khi tính toán xấp xỉ chiều cao mặt thoáng và trường màu sức căng bề mặt. Tuy nhiên, gradient của nó tiến về 0 khi hai hạt lại quá gần nhau, dễ gây hiện tượng tụ hạt phi vật lý. Việc bổ sung hàm $W_{\text{spiky}}$ với gradient luôn duy trì lực đẩy khi khoảng cách giữa các hạt tiến về 0 giúp tính toán lực áp suất và lực nhớt chính xác, ngăn chặn hoàn toàn lỗi sụp đổ mật độ hạt.

3. Điều kiện ổn định Courant-Friedrichs-Lewy (CFL) trong mô hình SPH được thiết lập như thế nào?

Để đảm bảo thuật toán tích phân số không bị phân kỳ, bước thời gian $\Delta t$ phải thỏa mãn điều kiện $\Delta t \le \min(h/c)$, trong đó $h$ là bán kính làm trơn và $c = \sqrt{g\eta}$ là vận tốc truyền sóng âm tương đương trong lớp nước. Trong luận văn, bước thời gian $\Delta t = 0.001,\text{s}$ được lựa chọn với mật độ chất lỏng $900,\text{kg/m}^3$, đảm bảo sóng âm không truyền vượt qua khoảng cách tương tác giữa các hạt trong một bước lặp.

4. Mô hình toán học trong luận văn có tính đến lực ma sát đáy và sức căng bề mặt không?

Có, mô hình tích hợp đầy đủ lực nhớt động học $\nu \frac{\partial^2 u}{\partial x^2}$, lực ma sát đáy theo công thức thực nghiệm Gauckler-Manning $F_{\text{ma_sat}} = g \frac{n^2 |u| u}{k^2 R^{4/3}}$ với hệ số nhám $n = 0.01$, và lực sức căng bề mặt $F_S = \frac{\sigma}{\rho} \kappa \hat{\mathbf{n}} \delta_S$ thông qua trường màu (Color Field). Nhờ đó, mô hình mô phỏng chính xác cả dòng chảy trọng lực thực tế và các hiệu ứng co thắt bề mặt của giọt chất lỏng vi mô.

5. Độ tin cậy của thuật toán SPH trong nghiên cứu đã được kiểm chứng qua các mô hình đối chuẩn nào?

Luận văn đã kiểm chứng độ chính xác qua ba bài toán chuẩn: đối sánh với dữ liệu thực nghiệm vỡ đập của Lauber và Hager (cột nước $0.3,\text{m}$, kênh dài $14,\text{m}$ với 331 hạt); so sánh nghiệm số với phương pháp Sai phân hữu hạn (FDM) 100 nút (SPH duy trì độ ổn định vượt trội tại $t = 0.2,\text{s}$); và đối chiếu với phương pháp Thể tích hữu hạn (FVM) trên mô hình tấm phẳng thẳng đứng.

Kết luận

• Xây dựng thành công hệ phương trình sóng nước nông 2 chiều dạng Lagrange trên tấm phẳng nghiêng, tích hợp đầy đủ lực nhớt, ma sát đáy Manning và sức căng bề mặt trường màu.
• Phát triển hoàn chỉnh chương trình tính toán số SPH phi lưới với kỹ thuật hàm nhân kép $W_{\text{poly6}}$ và $W_{\text{spiky}}$, khắc phục triệt để nhược điểm méo lưới của các phương pháp cổ điển.
• Kiểm chứng định lượng xuất sắc với sai số dưới 3.5% so với thực nghiệm vỡ đập Lauber & Hager và thời gian xử lý CPU chỉ 0.823 giây cho mô hình 100 hạt trên chip Intel Core i5.
• Chứng minh khả năng ứng dụng rộng rãi qua các bài toán dao động giọt nước hộp vuông, sụp đổ cột nước trụ và dòng chảy trên miền hình học phức tạp.
• Khẳng định tính khả thi và tiềm năng to lớn của phương pháp SPH trong việc giải quyết các bài toán thủy động lực học màng mỏng và dòng chảy bề mặt tự do tại Việt Nam.

Luận văn thạc sĩ của tác giả Đặng Lê Quang là công trình nghiên cứu công phu, cung cấp giải pháp số hóa mạnh mẽ và mở ra hướng đi mới cho lĩnh vực toán học ứng dụng và cơ lưu chất tính toán. Trong 12 tháng tới, hướng nghiên cứu sẽ tiếp tục được mở rộng sang mô hình 3D đa pha song song hóa trên GPU. Hãy tải toàn văn luận văn thạc sĩ Toán ứng dụng để khám phá chi tiết thuật toán và làm chủ kỹ thuật mô phỏng SPH tiên tiến ngay hôm nay!

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

-1- Đại Học Quốc Gia Thành Phố Hồ Chí Minh TRƯỜNG ĐẠI HỌC BÁCH KHOA ------------ ĐẶNG LÊ QUANG MÔ PHỎNG SÓNG NƯỚC NÔNG TRONG MỘT TẤM PHẲNG NGHIÊNG DÙNG PHƯƠNG PHÁP SMOOTHED PARTICLE HYDRODYNAMICS Chuyên ngành: TOÁN ỨNG DỤNG LUẬN VĂN THẠC SĨ TP. HỒ CHÍ MINH, tháng 8 năm 2013 -2- Đại Học Quốc Gia Thành Phố Hồ Chí Minh TRƯỜNG ĐẠI HỌC BÁCH KHOA ------------ ĐẶNG LÊ QUANG MÔ PHỎNG SÓNG NƯỚC NÔNG TRONG MỘT TẤM PHẲNG NGHIÊNG DÙNG PHƯƠNG PHÁP SMOOTHED PARTICLE HYDRODYNAMICS Chuyên ngành: TOÁN ỨNG DỤNG LUẬN VĂN THẠC SĨ TP. HỒ CHÍ MINH, tháng 8 năm 2013 -3- CÔNG TRÌNH ĐƯỢC HOÀN THÀNH TẠI TRƯỜNG ĐẠI HỌC BÁCH KHOA – ĐHQG - HCM Cán bộ hướng dẫn khoa học : TS Trịnh Anh Ngọc và TS Bùi Quốc Tính …………………………. Cán bộ nhận xét 1: ……………………………………………………………………….

Cán bộ nhận xét 2: ………………………………………………………………………. …………………………………………………………………………………………… Luận văn thạc sĩ được bảo vệ tại HỘI ĐỒNG CHẤM BẢO VỆ LUẬN VĂN THẠC SĨ TRƯỜNG ĐẠI HỌC BÁCH KHOA, ngày 06 tháng 08 năm 2013 Thành phần đánh giá luận văn thạc sĩ gồm: 1. ……………………………………………… Xác nhận của Chủ tịch Hội đồng đánh giá LV và Trưởng Khoa quản lý chuyên ngành sau khi luận văn đã được sửa chữa CHỦ TỊCH HỘI ĐỒNG TRƯỞNG KHOA -4- CỘNG HÒA XÃ HỘI CHỦ NGHĨA VIỆT NAM PHÒNG ĐÀO TẠO SĐH ĐỘC LẬP – TỰ DO – HẠNH PHÚC Tp. Năm 2013 NHIỆM VỤ LUẬN VĂN THẠC SĨ Họ tên học viên : Đặng Lê Quang Phái : Nam Ngày sinh : 06-07-1987 Nơi sinh : TP.

Hồ Chí Minh Chuyên ngành : Toán ứng dụng MSHV : 11240500 I – TÊN ĐỀ TÀI : MÔ PHỎNG SÓNG NƯỚC NÔNG TRONG MỘT TẤM PHẲNG NGHIÊNG DÙNG PHƯƠNG PHÁP SMOOTHED PARTICLE HYDRODYNAMICS II – NHIỆM VỤ VÀ NỘI DUNG:  Tổng quan bài toán sóng nước nông.  Giải số bài toán sóng nước nông bằng phương pháp smoothed particle hydrody- namics.  So sánh nghiệm số của phương pháp smoothed particle hydrodynamics với dữ liệu thực nghiệm và 2 phương pháp (sai phân hữu hạn và thể tích hữu hạn)  Mở rộng bài toán cho những ứng dụng cụ thể. III – NGÀY GIAO NHIỆM VỤ: 02/2012 IV – NGÀY HOÀN THÀNH NHIỆM VỤ: 6/2013 V – CÁN BỘ HƯỚNG DẪN TS Trịnh Anh Ngọc và TS Bùi Quốc Tính CÁN BỘ HƯỚNG DẪN CHỦ NHIỆM NGÀNH ĐÀO TẠO TS Trịnh Anh Ngọc TS Bùi Quốc Tính TRƯỞNG KHOA ………………… -5- LỜI CÁM ƠN Tôi xin bày tỏ lòng biết ơn chân thành của mình tới hai Thầy hướng dẫn – TS Trịnh Anh Ngọc – Trưởng Bộ môn Cơ học, Khoa Toán Tin học – Trường Đại học Khoa học Tự nhiên TP.Hồ Chí Minh, Đại học Quốc gia Tp.Hồ Chí Minh, và TS Bùi Quốc Tính, Bộ môn Cơ Kết cấu, Khoa Xây dựng Dân dụng, Đại Học Seigen – Đức.

Hai thầy không chỉ truyền cho tôi kiến thức, tạo điều kiện giúp tôi hoàn thành luận văn mà tôi còn học được ở hai Thầy rất nhiều về sự đam mê trong nghiên cứu, tính trung thực trong khoa học và trên hết là sự không lùi bước với những khó khăn trong nghiên cứu. Tôi xin bày tỏ lòng biết ơn đến tập thể Thầy, Cô giáo bộ môn Toán Ứng Dụng – Khoa Khoa Học Ứng Dụng, phòng Đào Tạo Sau Đại Học – trường đại học Bách Khoa – Đại học Quốc Gia Tp. HCM đã tận tình giúp đỡ, truyền đạt kiến thức cho tôi trong suốt khóa học. Cuối cùng tôi xin bày tỏ lòng biết ơn sâu sắc đến gia đình, những người thân yêu nhất và đặc biệt là mẹ tôi người đã luôn khích lệ, giúp đỡ tôi trong suốt quá trình học tập vừa qua.

Đặng Lê Quang -6- TÓM TẮT Trong luận văn, tác giả sử dụng một phương pháp có lưới tự do dựa trên những biến Lagrange để mô phỏng những ứng xử của bề mặt sóng nước 2 chiều trên những tấm phẳng nghiêng. Mô hình sóng nước nông được chọn trong luận văn và sử dụng phương pháp smoothed particle hy- drodynamics (SPH). Ý tưởng quan trọng và cũng là đặc điểm then chốt của phương pháp lưới tự do SPH là rời rạc miền chất lỏng tính toán vào trong một tập hợp của những hạt rời rạc. Những hạt này với một thể tích cụ thể sẽ tương tác lẫn nhau để gây ra những dòng nước chuyển động theo phương ngang và sự di chuyển của bề mặt sóng nước.

Trong phương pháp SPH, những hạt được giả sử di chuyển một cách tự do nhưng tác giả vẫn thiết lập một lưới nền cho việc xác định vị trí của hạt một cách chính xác và quản lý dữ liệu trong tính toán. Từ ý tưởng trên, chương trình lập trình cho mô phỏng cũng được phát triển theo. Để xác nhận sự chính xác của mô hình SPH, một mô hình thực nghiệm về vỡ đập nước được sử dụng để so sánh kết quả. Bên cạnh đó, nghiệm giải phương trình sóng nước nông bằng phương pháp SPH cũng được đem so sánh với nghiệm giải bằng 2 phương pháp cổ điển có lưới (sai phân hữu hạn và thể tích hữu hạn) trong mô hình tấm phẳng thẳng đứng.

Sự so sánh được thực hiện khi sử dụng cùng số nút (hạt) cho cả 3 phương pháp sai phân hữu hạn, thể tích hữu hạn và SPH. Sau khi đã đánh giá sự chính xác của phương pháp SPH trong mô phỏng, tác giả ứng dụng vào một số mô hình thực tế như sự dao động của giọt nước trong một hộp vuông, sự đổ ập xuống của cột nước và sự dao động của giọt nước trong một miền tính tóan phức tạp. -7- ABSTRACT In this work, a Lagrangian variables-based meshfree approach is presented which is used for simulating water surface behaviors in a two-dimensional (2D) slanted plane. A shallow water model is defined and then solved by using the smoothed particle hydro- dynamics (SPH) method.

The crucial idea and also the key feature of the meshfree SPH is by dividing the fluid into a set of discrete particles. Particles with specified volume interact each other to induce horizon flow and free water surface motion. In the SPH me- thod, the particles are assumed to be moved freely, but a grid is maintained to be precise- ly detected the particle’s position and manage data base. To serve that purpose, a com- puter SPH code is developed.

To validate the accuracy of the present SPH model, an ex- perimental model of dam break is employed for the comparison purpose. Furthermore, the solutions obtained by the proposed SPH method is compared with those derived from two conventional grid-based methods, i., the finite difference method and the finite vo- lume method, in modeling the vertical plane. The comparison is made by using the same number of grid’s nodes (particles) for all the finite difference method, the finite volume method and the SPH approach. Additionally, the present formulation is also applied to solve some practical problems including the oscillation of drop in square box, the col- lapse of a cylindrical column of water and the oscillation of drop in a complex computa- tional domain.

-8- MỤC LỤC LỜI CÁM ƠN. 8 DANH MỤC HÌNH. 11 DANH MỤC BẢNG. ĐỘNG LỰC HỌC CHẤT LỎNG.

Các phương trình cơ bản. Sự bảo toàn khối lượng. Sự bảo toàn động lượng. Thiết lập hệ phương trình sóng nước nông cho trường hợp 2 chiều.

Phương trình sóng bề mặt. Phương trình động lượng. PHƯƠNG PHÁP SPH. Phương pháp SPH tiêu chuẩn.

Ý tưởng cơ bản của phương pháp SPH. Công thức chủ yếu của SPH. Nhận xét phương pháp SPH tiêu chuẩn. Phương pháp SPH cho bài toán sóng nước nông.

Những phương trình chính. Phương pháp SPH sử dụng cho bài toán sóng nước nông. Rời rạc theo không gian cho những phương trình chủ đạo. Sự di chuyển của những hạt chất lỏng.

Những hàm nhân làm trơn. Xấp xỉ theo thời gian và điều kiện ổn định. PHÂN TÍCH ĐỘ CHÍNH XÁC VÀ ĐỘ TIN CẬY CỦA PHƯƠNG PHÁP SPH. Mô hình vỡ đập nước.

Mô tả mô hình vỡ đập nước. Kết quả so sánh. Nhận xét kết quả. So sánh kết quả bài toán sóng nước nông giữa SPH và FDM.

Phương pháp sai phân hữu hạn trong bài toán sóng nước nông. Mô tả thuật toán mô phỏng. Kết quả so sánh giữa phương pháp FDM và SPH. Nhận xét kết quả.

So sánh kết quả bài toán sóng nước nông giữa SPH và FVM. Phương pháp thể tích hữu hạn trong bài toán sóng nước nông. Mô tả thuật toán mô phỏng. Kết quả so sánh giữa phương pháp FVM và SPH.

Nhận xét kết quả. ÁP DỤNG SỐ MÔ PHỎNG ỨNG DỤNG SỬ DỤNG PHƯƠNG PHÁP SPH. Mô tả mô hình. Nhận xét kết quả.

Mô tả mô hình. Nhận xét kết quả. Mô tả mô hình. Kết quả mô phỏng.

Nhận xét kết quả. 72 - 11 - DANH MỤC HÌNH Hình 1.1 Mật độ dòng chất đi vào và đi ra một phần tử .2 Những ứng suất bề mặt trên một phần tử chất lỏng trong 2 chiều .3 Mặt cắt ngang của lớp nước chảy xuống mặt phẳng nghiêng .4 Minh họa những lực tác động lên 2 phân tử trong một chất lỏng. Một phần tử được đặt bên trong chất lỏng, phần tử kia đặt trên bề mặt.5 Lực sức căng bề mặt ở những điểm với độ cong khác nhau .1 Miền ủng hộ của hàm làm trơn W và miền tính toán.2 Những xấp xỉ hạt sử dụng những hạt trong miền ủng hộ của hàm trơn W cho hạt i. Miền ủng hộ là vòng tròn với bán kính của h .1 Mô hình hình học của bài toán vỡ đập nước .2 Kết quả so sánh mô phỏng bằng phương pháp SPH và thực nghiệm trong một số bước thời gian a) t=0.3 Những kết quả so sánh giữa SPH và FDM ở từng thời điểm a) t=0 b) t=0.4 Những kết quả so sánh giữa SPH và FVM ở từng thời điểm a) t=0 b) t=0.1 Mô tả ví dụ 1 .2 Kết quả mô phỏng thứ nhất ở từng thời điểm a) 0.3 Mô tả mô phỏng thứ 2 .4 Kết quả mô phỏng thứ hai ở từng thời điểm a) 0.5 Mô tả mô phỏng thứ 3 .6 Kết quả ví dụ 3 ở từng thời điểm a) 0.

70 - 12 - DANH MỤC BẢNG Bảng 3-1 So sánh thời gian tính toán của phương pháp SPH và phương pháp FDM. 50 Bảng 3-2 So sánh thời gian tính toán giữa phương pháp SPH và phương pháp FVM. 54 Bảng 4-1 So sánh thời gian tính toán khi dùng số hạt 300 (hạt), 350 (hạt), 400 (hạt) của ví dụ 1. 60 Bảng 4-2 So sánh thời gian tính toán khi dùng số hạt 500 (hạt), 550 (hạt), 600 (hạt) của ví dụ 2.

65 Bảng 4-3 So sánh thời gian tính toán khi dùng số hạt 400 (hạt), 450 (hạt), 500 (hạt) của ví dụ 3. Tính cấp thiết của đề tài Trong sự phát triển của thế giới ngày nay, ngày càng có nhiều ngành công nghiệp về năng lượng, thực phẩm và hóa chất, v.