(LUẬN văn THẠC sĩ) các phương pháp tạo lưới tự động và ứng dụng trong tính toán cơ học

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN Nguyễn Thị Thủy CÁC PHƯƠNG PHÁP TẠO LƯỚI TỰ ĐỘNG VÀ ỨNG DỤNG TRONG TÍNH TỐN CƠ HỌC LUẬN VĂN THẠC SỸ KHOA HỌC Hà Nội – 2011 TIEU LUAN MOI download : skknchat@gmail.com ĐẠI HỌC QUỐC GIA HÀ NỘI TRƯỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN Nguyễn Thị Thủy CÁC PHƯƠNG PHÁP TẠO LƯỚI TỰ ĐỘNG VÀ ỨNG DỤNG TRONG TÍNH TỐN CƠ HỌC LUẬN VĂN THẠC SỸ KHOA HỌC Chuyên ngành: Cơ học chất lỏng Mã số: 60 44 22 Người hướng dẫn khoa học: PGS.TS Trần Văn Trản Hà Nội – 2011 TIEU LUAN MOI download : skknchat@gmail.com Mục lục Lời cảm ơn Lời mở đầu Một số phương pháp chia lưới tự động không cấu trúc 1.1 Phương pháp chia lưới Delaunay Triangulation 1.1.1 Giới thiệu 1.1.2 Cơ sở hình học 1.1.3 Thiết lập hệ tam giác ban đầu 1.1.4 Thuật toán Bowyer - Watson 1.1.5 Các phương pháp chèn điểm 1.1.6 Hạn chế hình dạng hệ tam giác Delaunay 1.1.7 Phép chia lưới Delaunay triangulation không gian ba chiều 1.2 Phương pháp tịnh tiến biên (Advancing Front) 1.2.1 Giới thiệu 1.2.2 Điều khiển lưới 1.2.3 Thuật toán AFT 1.2.4 Sự thích nghi không gian tham số 1.2.5 Cải thiện chất lượng lưới 1.3 Phương pháp tạo lưới khơng cấu trúc sử dụng thuật tốn chèn điểm tự động tái kết nối địa phương 1.3.1 Giới thiệu 1.3.2 Phương pháp chia lưới AFLR 1.4 Các phương pháp chồng tạo lưới tứ giác lục giác 1.4.1 Giới thiệu 1.4.2 Các phương pháp chồng 1.5 Một số phương pháp phát triển 6 10 10 13 18 20 21 21 22 25 35 36 37 37 38 45 45 46 52 Áp dụng số toán học 2.1 Bài toán 2.2 Bài toán 2.3 Bài toán 54 54 58 59 TIEU LUAN MOI download2: skknchat@gmail.com 2.4 Bài toán 63 2.5 Bài toán 66 Tài liệu tham khảo 71 TIEU LUAN MOI download3: skknchat@gmail.com Lời mở đầu Ngày phương pháp chia lưới tự động trở thành công cụ phổ biến việc sử dụng nghiệm số phương trình vi phân phần miền có hình dạng Trong năm gần đây, phương pháp chia lưới tự động phát triển rộng rãi, học kết cấu học vật rắn, sau tính tốn động lực học chất lỏng Phương pháp chia lưới tự động cung cấp chìa khóa để giải vấn đề hình dạng biên từ phương pháp sai phân hữu hạn sử dụng việc xác định điểm nút phương pháp phần tử hữu hạn Với mạng lưới tất thuật toán số, sai phân hữu hạn hay phần tử hữu hạn thực miền tính tốn có hình dạng đường biên Các phương pháp chia lưới tự động bao gồm hai loại: phương pháp chia lưới tự động có cấu trúc phương pháp chia lưới tự động không cấu trúc Một ví dụ đơn giản lưới có cấu trúc lưới Đề-Các Các nút lưới lưới (hình chữ nhật) giao đường lưới Lưới tạo cách dễ dàng nhanh chóng Tuy nhiên miền tính tốn có biên biên ngồi đường cong để giải lưới Đề-Các yêu cầu phải xác định sơ đồ số gần biên Điều thường khó Hình 1: (a) Lưới Đề-Các, (b) body-fitted grid Một loại khác lưới có cấu trúc body-fitted (trùng khít với miền tính tốn) Về lưới giống lưới Đề-Các lưới tứ giác có hình dạng phù hợp với biên Các phương pháp tính tốn phát triển cho lưới Đề-Các chỉnh sửa để áp dụng cho lưới body-fitted thông TIEU LUAN MOI download4: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc qua ánh xạ lưới mà khơng làm thay đổi chất chương trình số Khó khăn việc thực phương pháp nằm việc tạo nút lưới có nhiều vật cản bên miền tính tốn Khơng giống lưới có cấu trúc, nút lưới ô lưới lưới không cấu trúc giao đường thẳng song song Một ưu điểm lưới không cấu trúc điểm nút đặt biên miền tính tốn áp dụng phương pháp số lưới khơng cấu trúc cho độ xác cao Hơn nữa, lưới không cấu trúc cho phép phần tử có kích thước khác biểu diễn xác biên miền tính tốn mà không cần số lượng lớn nút lưới lưới Khi miền tính tốn có hình dạng phức tạp so với lưới có cấu trúc, lưới khơng cấu trúc giúp cải thiện độ xác nghiệm tổng thể, lưới khơng cấu trúc chứng minh có khả thích nghi cao toán biên chuyển dịch dịng thời Về ngun tắc, lưới khơng cấu trúc bao gồm có hình dạng bất kỳ, xây dựng cách kết nối điểm cho trước đến số tùy ý điểm khác, nói chung hợp thành từ tam giác tứ giác không gian hai chiều, tứ diện lục giác không gian ba chiều Ưu điểm lưới chúng có khả thích nghi cách cho phép chèn thêm điểm vào, chúng làm việc với miền tính tốn có hình dạng phức tạp Ở thời điểm phương pháp chia lưới không cấu trúc áp dụng không gian ba chiều Tuy nhiên phương pháp chia lưới không cấu trúc dựa tam giác Delaunay thường xuyên sử dụng chứng minh có khả thích nghi cao, phù hợp với miền tính toán phức tạp Các lưới cho phép bổ sung thêm nút lưới vào hệ tam giác tồn ảnh hưởng đến cấu trúc lưới địa phương mà không ảnh hưởng đến cấu trúc lưới tổng thể Trong luận văn giới thiệu hai phương pháp chia lưới không cấu trúc hay sử dụng: phương pháp chia lưới Delaunay Triangulation phương pháp tịnh tiến biên (Advancing Front method) Ngoài hai phương pháp phát triển: Phương pháp tạo lưới không cấu trúc sử dụng thuật toán chèn điểm tự động tái kết nối địa phương phương pháp chồng tạo lưới tứ giác lục giác trình bày luận văn * Nội dung luận văn gồm chương: Chương 1: Một số phương pháp chia lưới tự động không cấu trúc Chương 2: Áp dụng số toán học TIEU LUAN MOI download5: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc Chương Một số phương pháp chia lưới tự động không cấu trúc 1.1 1.1.1 Phương pháp chia lưới Delaunay Triangulation Giới thiệu Phương pháp Delaunay Triangulation phương pháp chia lưới không cấu trúc phát triển từ sớm [8] Phương pháp dựa tiêu chuẩn Delaunay (hay gọi tiêu chuẩn vòng tròn ngoại tiếp trống): Siêu cầu đơn hình khơng gian n - chiều xác định n + điểm khơng có điểm nút khác lưới Ví dụ không gian ba chiều, bốn đỉnh tứ diện xác định mặt cầu mặt cầu không chứa điểm nút khác lưới tứ diện Trong không gian hai chiều, hệ tam giác thu dựa tiêu chuẩn Delaunay gọi hệ tam giác Delaunay Hệ tam giác áp dụng phổ biến thực hành chúng có đặc điểm tối ưu sau: • Các tam giác Delaunay tam giác xấp xỉ đều; • Góc lớn tam giác cực tiểu hóa; • Góc nhỏ tam giác cực đại hóa Với đặc điểm hệ tam giác Delaunay không bị biến dạng méo mó Tiêu chuẩn Delaunay không đưa dẫn điểm lưới nên định nghĩa liên kết với nào? Một hạn chế tiêu chuẩn Delaunay có khả khơng thể áp dụng tiêu chuẩn tồn miền tính tốn với tam giác biên xác định trước Nhược điểm đưa hai cách tiếp cận chia lưới tam giác có bảo tồn liên kết biên áp dụng tiêu chuẩn Delaunay Trong cách tiếp cận thứ tiêu chuẩn Delaunay bỏ qua điểm gần biên hệ TIEU LUAN MOI download6: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc biên lưới trước nguyên vẹn Để kết hợp với kỹ thuật này, điểm thêm vào dạng sơ đồ để đảm bảo không xảy phá hủy biên Cách tiếp cận thứ hai, áp dụng tiêu chuẩn Delaunay tồn miền tính tốn, sau khơi phục lại biên ban đầu cách bỏ đơn hình nằm bên ngồi miền tính tốn [1] Có nhiều thuật tốn tạo lưới khơng cấu trúc dựa tiêu chuẩn Delaunay, chẳng hạn có số thuật tốn sử dụng phương pháp chia lưới có cấu trúc tạo phân bố điểm nút lưới trước sau điểm nút lưới kết nối để thu tam giác thỏa mãn tiêu chuẩn hình học (tương đương với tiêu chuẩn Delaunay) Tuy nhiên thuật toán hay sử dụng thuật tốn Bowyer - Watson Thuật tốn áp dụng với không gian n - chiều Thuật toán hệ tam giác vài điểm, sau tiếp tục bước ta thêm điểm vào hệ tam giác tái thiết lập hệ tam giác cách địa phương Quá trình cho phép cải thiện chất lượng lưới khuôn khổ tiêu chuẩn Delaunay Điểm khác biệt thuật toán vị trí điểm liên kết tính tốn cách đồng thời 1.1.2 Cơ sở hình học • Định nghĩa ô lồi Một ô lồi n - chiều S bao lồi n + k điểm P1 , , Pn+k (k>1) mà điểm không nằm mặt phẳng (n − 1) - chiều Như S bao gồm điểm x ∈ Rn thỏa mãn n+k x= ∑ αi Pi , i =1 n+k ∑ αi = 1, i =1 ≥ αi ≥ Gọi tất điểm Pl tập Pi , i = 1, ., n + k nằm biên S đỉnh ô lồi S Một mặt m - chiều ô lồi n - chiều S (n > m) gọi bao lồi m + đỉnh Pl , bao lồi không chứa đỉnh khác S Ta nói S lồi mạnh khơng có hai mặt nằm mặt phẳng m - chiều với m < n Nếu P điểm nằm ô lồi mạnh S với đỉnh P1 , , Pn+k n+k P= i ∑ α Pi , i =1 n+k ∑ αi = 1, i =1 αi ≥ 0, i = 1, , n + k, TIEU LUAN MOI download7: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc Hình 1.1: Ơ lồi (bên trái) tứ diện lồi mạnh (bên phải) • Đơn hình đơn hình Một phần tử đơn giản n - chiều sử dụng để rời rạc hóa miền tính tốn gọi n - chiều Ô bao n + điểm x1 , ., xn+1 mà không nằm mặt phẳng (n − 1) - chiều Những gọi đơn hình Như đơn hình tạo nên điểm x ∈ Rn thỏa mãn n +1 x= ∑ αi xi , i = 1, , n + 1, i =1 n +1 ∑ αi = 1, i =1 αi ≥ 0, Đơn hình lồi mạnh có đỉnh x1 , ., xn+1 Chẳng hạn đơn hình ba chiều tứ diện có đỉnh x1 , x2 , x3 , x4 , đơn hình hai chiều tam giác, đơn hình chiều đoạn thẳng Mỗi mặt m chiều đơn hình đơn hình m - chiều định nghĩa qua m + đỉnh Điểm x điểm đơn hình α > với i = 1, , n + Trong thực hành, để rời rạc hóa miền tính tốn, ta hay sử dụng lồi có mặt biên đơn hình Những ô gọi ô đơn hình Gọi Ni , i = 1, , n, số mặt đơn hình i - chiều S, N0 số đỉnh S Ta có   n −1 i +  Ni = (−1) n−1 Nk , ∑ (−1)i  k+1 i =k k = −1, , n − 2, N−1 = 1,     l l  = l (l − 1) ( l − m + 1) , m ≥ 1,   = 1,  m! m TIEU LUAN MOI download8: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc • Tính qn lưới Bằng phép rời rạc phù hợp thu tập hợp điểm V ∈ Rn tập ô lồi mạnh T thỏa mãn điều kiện sau: Tập hợp đỉnh ô T trùng với V; Nếu hai ô khác S1 S2 giao nhau, miền giao mặt chung hai Hình 1.2: Các ô giao chấp nhận (a) không chấp nhận (b, c, d) Tập hợp ô phép rời rạc phù hợp tạo thành miền kết nối đơn giản n - chiều Gọi Ni , i > số lượng mặt biên i - chiều miền rời rạc, N0 số đỉnh biên, theo định lý Euler ta có: n −1 ∑ (−1)i Ni = + (−1)n−1 i =0 biểu thức sử dụng để xác định tính qn lưới • Lưới tổ ong Dirichlet Xét tập tùy ý gồm điểm Pi , i = 1, 2, , N miền xác định n - chiều Với điểm Pi xác định miền V ( Pi ) Rn bao gồm điểm có khoảng cách tới Pi nhỏ tới điểm Pj khác 
Vi = x ∈ Rn | d ( x, Pi ) ≤ d x, Pj , i 6= j, j = 1, ., N , d( a, b) khoảng cách hai điểm a, b Các miền Vi gọi khối đa diện Voronoi Do khối đa diện giao bán không gian nên chúng đa diện lồi, không cần thiết bị chặn Mặt biên chung hai Voronoi hai điểm V ( Pi ) V ( Pj ) đa giác (n − 1) - chiều Cặp điểm Pi Pj gọi cặp cấu hình khối đa diện Voronoi chúng có mặt chung Bằng cách kết nối điểm kề ta thu lưới Trong lưới này, tập hợp n + điểm kề với điểm khác tạo thành đơn hình n - chiều Tâm đơn hình TIEU LUAN MOI download9: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc • Trường hợp 3: chia lưới dày nguồn nhiệt, cửa hút cửa xả File input: 10 0.0 0.0 4.0 0.0 10 0.0 0.0 0 4.0 0.0 6.0 0.0 40 0.0 0.0 0 6.0 0.0 10 0.0 10 0.0 0.0 0 10 0.0 10 4.0 10 0.0 0.0 0 10 4.0 10 6.0 40 0.0 0.0 0 10 6.0 10 10 15 0.0 0.0 0 10 10 0.0 10 30 0.0 0.0 0 0.0 10 0.0 6.0 10 0.0 0.0 0 0.0 6.0 0.0 4.0 40 0.0 0.0 0 0.0 4.0 0.0 0.0 15 0.0 0.0 0 TIEU LUAN MOI download57: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc 2.2 Bài tốn 2: Miền tính tốn hình chữ nhật có hình chữ nhật nhỏ File input: 11 0.0 0.0 6.0 0.0 15 0.0 0.0 0 6.0 0.0 6.0 10 20 0.0 0.0 0 6.0 10 4.0 10 10 0.0 0.0 0 4.0 10 4.0 8.0 0.0 0.0 0 4.0 8.0 4.0 2.0 60 0.0 0.0 0 4.0 2.0 2.0 2.0 25 0.0 0.0 0 2.0 2.0 2.0 8.0 60 0.0 0.0 0 2.0 8.0 4.0 8.0 25 0.0 0.0 0 4.0 8.0 4.0 10 0.0 0.0 4.0 10 0.0 10 15 0.0 0.0 0.0 10 0.0 0.0 20 0.0 0.0 0 TIEU LUAN MOI download58: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc 2.3 Bài tốn 3: Miền tính tốn hình chữ nhật có trụ trịn TIEU LUAN MOI download59: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc • Trường hợp 1: chia lưới dày nửa vòng tròn bên File input: 0.0 0.0 14 0.0 20 0.0 0.0 0 14 0.0 14 6.0 15 0.0 0.0 0 14 6.0 8.0 6.0 0 0.0 0.0 8.0 6.0 6.0 6.0 10 7.0 6.0 −180 6.0 6.0 8.0 6.0 40 7.0 6.0 −180 8.0 6.0 14 6.0 0.0 0.0 14 6.0 14 12 15 0.0 0.0 14 12 0.0 12 20 0.0 0.0 0 0.0 12 0.0 0.0 30 0.0 0.0 0 TIEU LUAN MOI download60: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc • Trường hợp 2: chia lưới dày nửa vòng tròn bên File input: 0.0 0.0 14 0.0 20 0.0 0.0 0 14 0.0 14 6.0 15 0.0 0.0 0 14 6.0 8.0 6.0 0 0.0 0.0 8.0 6.0 6.0 6.0 40 7.0 6.0 −180 6.0 6.0 8.0 6.0 10 7.0 6.0 −180 8.0 6.0 14 6.0 0.0 0.0 14 6.0 14 12 15 0.0 0.0 14 12 0.0 12 20 0.0 0.0 0 0.0 12 0.0 0.0 30 0.0 0.0 0 TIEU LUAN MOI download61: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc • Trường hợp 3: chia lưới dày xung quanh vòng tròn File input: 0.0 0.0 14 0.0 20 0.0 0.0 0 14 0.0 14 6.0 15 0.0 0.0 0 14 6.0 8.0 6.0 0 0.0 0.0 8.0 6.0 6.0 6.0 40 7.0 6.0 −180 6.0 6.0 8.0 6.0 40 7.0 6.0 −180 8.0 6.0 14 6.0 0.0 0.0 14 6.0 14 12 15 0.0 0.0 14 12 0.0 12 20 0.0 0.0 0 0.0 12 0.0 0.0 30 0.0 0.0 0 TIEU LUAN MOI download62: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc 2.4 Bài toán 4: Miền tính tốn hình chữ nhật có hai trụ trịn nằm ngang • Trường hợp 1: chia lưới dày xung quanh hình trụ thứ File input: 13 0.0 0.0 19 0.0 25 0.0 0.0 0 19 0.0 19 6.0 15 0.0 0.0 0 19 6.0 14 6.0 0 0.0 0.0 14 6.0 12 6.0 30 13 6.0 −180 12 6.0 7.0 6.0 0.0 0.0 0 7.0 6.0 5.0 6.0 10 6.0 6.0 −180 5.0 6.0 7.0 6.0 10 6.0 6.0 −180 7.0 6.0 12 6.0 0.0 0.0 12 6.0 14 6.0 30 13 6.0 −180 14 6.0 19 6.0 0.0 0.0 19 6.0 19 12 15 0.0 0.0 19 12 0.0 12 25 0.0 0.0 0 0.0 12 0.0 0.0 20 0.0 0.0 0 TIEU LUAN MOI download63: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc • Trường hợp 2: chia lưới dày xung quanh hình trụ thứ hai File input: 13 0.0 0.0 19 0.0 25 0.0 0.0 0 19 0.0 19 6.0 15 0.0 0.0 0 19 6.0 14 6.0 0 0.0 0.0 14 6.0 12 6.0 10 13 6.0 −180 12 6.0 7.0 6.0 0.0 0.0 0 7.0 6.0 5.0 6.0 30 6.0 6.0 −180 5.0 6.0 7.0 6.0 30 6.0 6.0 −180 7.0 6.0 12 6.0 0.0 0.0 12 6.0 14 6.0 10 13 6.0 −180 14 6.0 19 6.0 0.0 0.0 19 6.0 19 12 15 0.0 0.0 19 12 0.0 12 25 0.0 0.0 0 0.0 12 0.0 0.0 20 0.0 0.0 0 TIEU LUAN MOI download64: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc • Trường hợp 3: chia lưới dày xung quanh hai hình trụ File input: 13 0.0 0.0 19 0.0 25 0.0 0.0 0 19 0.0 19 6.0 15 0.0 0.0 0 19 6.0 14 6.0 0 0.0 0.0 14 6.0 12 6.0 30 13 6.0 −180 12 6.0 7.0 6.0 0.0 0.0 0 7.0 6.0 5.0 6.0 30 6.0 6.0 −180 5.0 6.0 7.0 6.0 30 6.0 6.0 −180 7.0 6.0 12 6.0 0.0 0.0 12 6.0 14 6.0 30 13 6.0 −180 14 6.0 19 6.0 0.0 0.0 19 6.0 19 12 15 0.0 0.0 19 12 0.0 12 25 0.0 0.0 0 0.0 12 0.0 0.0 20 0.0 0.0 0 TIEU LUAN MOI download65: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc 2.5 Bài tốn 5: Miền tính tốn hình chữ nhật có hai trụ trịn thẳng đứng • Trường hợp 1: chia lưới dày xung quanh hình trụ thứ TIEU LUAN MOI download66: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc File input: 13 0.0 0.0 6.0 0.0 15 0.0 0.0 0 6.0 0.0 6.0 5.0 0 0.0 0.0 6.0 5.0 6.0 7.0 30 6.0 6.0 −180 6.0 7.0 6.0 12 0.0 0.0 0 6.0 12 6.0 14 10 6.0 13 −180 6.0 14 6.0 12 10 6.0 13 −180 6.0 12 6.0 7.0 0.0 0.0 6.0 7.0 6.0 5.0 30 6.0 6.0 −180 6.0 5.0 6.0 0.0 0.0 0.0 6.0 0.0 12 0.0 15 0.0 0.0 12 0.0 12 19 30 0.0 0.0 0 12 19 0.0 19 30 0.0 0.0 0 0.0 19 0.0 0.0 30 0.0 0.0 0 TIEU LUAN MOI download67: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc • Trường hợp 2: chia lưới dày xung quanh hình trụ thứ hai File input: 13 0.0 0.0 6.0 0.0 15 0.0 0.0 0 6.0 0.0 6.0 5.0 0 0.0 0.0 6.0 5.0 6.0 7.0 10 6.0 6.0 −180 6.0 7.0 6.0 12 0.0 0.0 0 6.0 12 6.0 14 30 6.0 13 −180 6.0 14 6.0 12 30 6.0 13 −180 6.0 12 6.0 7.0 0.0 0.0 6.0 7.0 6.0 5.0 10 6.0 6.0 −180 6.0 5.0 6.0 0.0 0.0 0.0 6.0 0.0 12 0.0 15 0.0 0.0 12 0.0 12 19 30 0.0 0.0 0 12 19 0.0 19 30 0.0 0.0 0 0.0 19 0.0 0.0 30 0.0 0.0 0 TIEU LUAN MOI download68: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc • Trường hợp 3: chia lưới dày xung quanh hai hình trụ File input: 13 0.0 0.0 6.0 0.0 15 0.0 0.0 0 6.0 0.0 6.0 5.0 0 0.0 0.0 6.0 5.0 6.0 7.0 30 6.0 6.0 −180 6.0 7.0 6.0 12 0.0 0.0 0 6.0 12 6.0 14 30 6.0 13 −180 6.0 14 6.0 12 30 6.0 13 −180 6.0 12 6.0 7.0 0.0 0.0 6.0 7.0 6.0 5.0 30 6.0 6.0 −180 6.0 5.0 6.0 0.0 0.0 0.0 6.0 0.0 12 0.0 15 0.0 0.0 12 0.0 12 19 30 0.0 0.0 0 12 19 0.0 19 30 0.0 0.0 0 0.0 19 0.0 0.0 30 0.0 0.0 0 TIEU LUAN MOI download69: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc Kết Luận Luận văn trình bày số phương pháp chia lưới khơng cấu trúc tự động, hai phương pháp hay sử dụng phương pháp chia lưới Delaunay Triangulation phương pháp tịnh tiến biên (Advancing Front method) Áp dụng vào số toán học cho kết tốt Tuy dừng trường hợp hai chiều thời gian có hạn kết luận văn sử dụng thực tế giải toán học phương pháp phần tử hữu hạn thể tích hữu hạn Hướng phát triển luận văn: Tiếp tục mở rộng sang toán ba chiều tìm hiểu thêm số phương pháp chia lưới tự động phát triển TIEU LUAN MOI download70: skknchat@gmail.com (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc (LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc(LUAN.van.THAC.si).cac.phuong.phap.tao.luoi.tu.dong.va.ung.dung.trong.tinh.toan.co.hoc