Home / Luận văn thạc sĩ / Đề tài: Phương pháp RBIEM với miền địa phương tròn giải hệ, HOT

Đề tài: Phương pháp RBIEM với miền địa phương tròn giải hệ, HOT

Chương 1 Giới thiệu tổng quan Phương pháp phần tử biên (BEM) để giải phương trình Navier-Stokes là một trong những bài toán được các nhà khoa học quan tâm. Khi dùng phương trình tích phân biên, số hạng phi tuyến xuất hiện trong tích phân miền. Có nhiều phương pháp khác nhau để giải số hạng phi tuyến đó như Zheng et al. [11] dùng phương pháp nghiệm riêng, Power và Partridge [7] sử dụng phương pháp đối ngẫu tương hỗ (DRM). Nhưng kết hợp giữa BEM và DRM chỉ giải được các bài toán dòng chảy phức tạp với số Reynolds nhỏ bằng 40 hay 100. Bằng phương pháp phân chia miền con [4, 8] Power và Mingo đã giải bài toán cho số Reynolds cao hơn với độ chính xác cao hơn. Tuy nhiên phương pháp BEM-DRM đã xấp xỉ đạo hàm của vận tốc trong số hạng phi tuyến thông qua hàm bán kính cơ sở và tạo ra phương trình đại số tuyến tính với số phương trình lơn hơn số ẩn làm tăng độ phức tạp của bài toán. Bên cạnh đó, phương pháp không lưới kết hợp với phương trình tích phân biên đang được quan tâm rộng rãi bởi tính chính xác mà phương trình tích phân biên mang lại. Trong đó phương pháp không lưới tích phân miền địa phương (LBIE) đưa ra bởi Zhu et al. [12, 13] giải bài toán Poison và bài toán phi tuyến dựa trên xấp xỉ dịch chuyển bình phương tối thiểu với ý tưởng tạo ra biên địa phương trên mỗi nút. Sau đó Sellountos và Sequeira [10] dùng LBIE để giải phương trình Navier-Stokes với cách tiếp cận dùng phương pháp nghiệm đi kèm để xấp xỉ số hạng phi tuyến. Gần đây, Popov và Bui [5] đưa ra phương pháp không lưới dựa trên phương trình tích phân biên và hàm bán kính cơ sở (RBIEM) để giải bài toán khuếch tán nhiễu, trong đó phương trình tích phân biên được áp dụng trên mỗi miền con địa 3
6. phương tương ứng với mỗi nút. Khi đó RBIEM tạo ra hệ phương trình đại số tuyến tính với số phương trình bằng số ẩn để giải, ma trận hệ số là ma trận thưa. RBIEM được áp dụng để giải hệ phương trình Navier-Stokes, trong đó với mỗi nút trong miền tính toán, có bảy ẩn số tương ứng với bảy phương trình tích phân biên. Thay vì phải xấp xỉ biến đạo hàm riêng của vận tốc ∂ui ∂xh hàm bán kính cơ sở. Ý tưởng của phương pháp RBIEM là xây dựng một miền con địa phương ứng với mỗi nút bên trong và trên biên miền tính toán. Về lý thuyết, những miền con địa phương này có thể có hình dạng bất kỳ. Khi đó để tích phân trên biên của miền bất kỳ, RBIEM phân rã biên thành những phần tử, tích phân trên biên địa phương sẽ được tính trên từng phần tử và sau đó được ghép lại. Trên thực tế, để thuận tiện trong quá trình tính toán, miền con được RBIEM tạo ra là những miền tròn. Nhưng khi đó, để tính tích phân biên có thể dùng phương pháp khác đơn giản hiệu quả hơn việc phân rã biên. Trong luận văn này, phương pháp không lưới RBIEM cải tiến được đề xuất. Để thuận tiện, ta gọi phương pháp RBIEM cải tiến là m-RBIEM (modified RBIEM). Để tính tích phân trên biên của miền con, thay việc rời rạc biên thành các phần tử bằng cách thêm vào các nút trên biên, phương pháp không lưới m-RBIEM sẽ sử dụng hệ tọa độ cực để tính trực tiếp các tích phân khi miền con có dạng hình tròn. Phương pháp m-RBIEM đưa ra lời giải số chính xác hơn, tiết kiệm thời gian tính toán hơn và dễ dàng hơn trong việc lập trình giải các bài toán thực tế. Cấu trúc luận văn được trình bày như sau: – Chương 1: Giới thiệu tổng quan về phương pháp không lưới dùng phương trình tích phân biên. – Chương 2: Đề cập phương pháp không lưới RBIEM giải phương trình Navier-Stokes. – Chương 3: Phương pháp RBIEM với miền địa phương tròn giải hệ phương trình Navier- Stokes. – Chương 4: Kết quả số. 4
7. Chương 2 Phương pháp không lưới RBIEM giải phương trình Navier-Stokes 2.1 Phương trình tích phân biên và phương pháp đối ngẫu tương hỗ Phương pháp đối ngẫu tương hỗ DRM (Dual Reciprocity Method) được kết hợp với phương pháp phương trình tích phân biên BEM (Boundary Element Method) dùng để chuyển số hạng tích phân miền thành tích phân trên biên khi giải phương trình Navier-Stokes. Xét phương trình Navier-Stokes cho chất lỏng không nén được: ρ ∂ui ∂t +ρuj ∂ui ∂xj = ∂σij ∂xj +ρFi; ∂ui ∂xi = 0, (2.1) trong đó: ui: là thành phần vectơ vận tốc theo hướng i; ρ: là mật độ; Fi: là lực tác động theo hướng i; σij: là tensơ ứng suất tương ứng trường vận tốc và áp suất (ui, p). 5
8. 2.1. PHƯƠNG TRÌNH TÍCH PHÂN BIÊN VÀ PHƯƠNG PHÁP ĐỐI NGẪU TƯƠNG HỖ Với chất lỏng Newton ta có: σij = −pδij + µ ( ∂ui ∂xj + ∂uj ∂xi ) , (2.2) trong đó: p: là áp suất chất lỏng; δij: là ký hiệu Kronecker; µ: là hệ số nhớt. Phương trình Navier-Stokes cho một điểm x trong miền Ω đóng bởi biên S dưới dạng tích phân được đưa ra bởi Ladyzhenskaya (1963): uk (x) = ∫ S t∗ ki (x,y)ui (y)dSy − ∫ S u∗ ki (x,y)ti (y)dSy + ∫ Ω u∗ ki (x,y)gidΩ, (2.3) trong đó: gi = ρujui,j: là số hạng phi tuyến; ti = σijnj, nj: là vectơ pháp tuyến hướng ra ngoại miền S; uk i : là trường nghiệm vectơ vận tốc của phương trình Stokes. Trong trường hợp hai chiều nghiệm u∗ ki và qk có dạng: u∗ ki (x,y) = − 1 4πµ [ ln ( 1 r ) δik + (xi −yi)(xk −yk) r2 ] ; qk (x,y) = − 1 2π (xk −yk) r2 , (2.4) trong đó r = |x−y|. Nghiệm cơ bản t∗ ki có dạng: t∗ ki = − 1 πr (xi −yi)(xk −yk) ( xj −yj ) r3 nj. (2.5) Khai triển số hạng gi (x) để xấp xỉ tích phân miền trong phương trình (2.3) thành tích phân biên dạng: gi (x) = ND ∑ m=1 fm (x)αm l δil, (2.6) 6
9. 2.1. PHƯƠNG TRÌNH TÍCH PHÂN BIÊN VÀ PHƯƠNG PHÁP ĐỐI NGẪU TƯƠNG HỖ trong đó fm (x) là hàm bán kính cơ sở phụ thuộc vào bán kính điểm cần xấp xỉ x và điểm lân cận ym, m = 1,…,N. Hàm fm (x) chỉ phụ thuộc vào giá trị R = |x−ym| là khoảng cách từ điểm x đến điểm

MÃ TÀI LIỆU:3589

 

  • PHÍ TÀI LIỆU: 25.000
  • ĐỊNH DẠNG: WORD+PDF
  • THANH TOÁN QUA MOMO, CHUYỂN KHOẢN, THẺ CÀO ĐIỆN THOẠI (X2)
  • NỘI DUNG: MÃ TÀI LIỆU – EMAIL NHẬN ( VÍ DỤ: 0324 – trinhnam34gmailcom) có thể bỏ chữ @ mới gửi được)
  • CHECK EMAIL (1-15 PHÚT)

  • Đăng nhập MOMO
  • Quét mã QR
  • Nhập số tiền
  • Nội dung: Mã Tài liệu – Email 
  • Check mail (1-15p)

  • Mua thẻ cào chỉ Viettel,  Vinaphone
  • Mệnh giá gấp 2 phí tài liệu (vì phí nhà mạng 50%) 
  • Add Zalo 0932091562
  • Nhận file qua zalo, email

  • Đăng nhập Internet Mobile
  • Chuyển tiền
  • Nhập số tiền
  • Nội dung: Mã Tài liệu – Email
  • Check mail (1-15p)

NẾU CHỜ QUÁ 15 PHÚT CHƯA THẤY MAIL VUI LÒNG NHẮN ZALO: 0932091562

 

 

NHẬP TÀI LIỆU BẠN CẦN TÌM VÀO ĐÂY


Notice: Undefined index: hide_title in /home/cdmhewql/tailieumau.vn/wp-content/plugins/wp-google-search/wgs-widget.php on line 26

Notice: Undefined index: title in /home/cdmhewql/tailieumau.vn/wp-content/plugins/wp-google-search/wgs-widget.php on line 28

 

 

About Trang web mặc định

Check Also

Đề tài: Biện pháp nâng cao lợi nhuận tại Công ty cổ phần Đông Đô

LỜI NÓI ĐẦU 1. Tính cấpthiết của đề tài Việc chuyển mình từ nền kinh …

Trả lời

Email của bạn sẽ không được hiển thị công khai. Các trường bắt buộc được đánh dấu *