BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ NGUYỄN THỊ QUỲNH TRÂM ỨNG DỤNG PHƯƠNG PHÁP PROPER GENERALIZED DECOMPOSITION VÀ PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN CHO BÀI TOÁN LƯU CHẤT NGÀNH: CÔNG NGHỆ CHẾ TẠO MÁY - 605204 S KC 0 0 4 0 4 2 Tp. Hồ Chí Minh, tháng 09 năm 2013 Luan van BỘ GIÁO DỤC VÀ ĐÀO TẠO TRƯỜNG ĐẠI HỌC SƯ PHẠM KỸ THUẬT THÀNH PHỐ HỒ CHÍ MINH LUẬN VĂN THẠC SĨ NGUYỄN THỊ QUỲNH TRÂM ỨNG DỤNG PHƯƠNG PHÁP PROPER GENERALIZED DECOMPOSITION VÀ PHƯƠNG PHÁP PHẦN TỬ HỮU HẠN CHO BÀI TOÁN LƯU CHẤT NGÀNH: CÔNG NGHỆ CHẾ TẠO MÁY – 605204 Hướng dẫn khoa học: TS. PHAN ĐỨC HUYNH Tp. Hồ Chí Minh, tháng 09 / 2013 Luan van LÝ LỊCH KHOA HỌC I.
LÝ LỊCH SƠ LƢỢC Họ & tên: NGUYỄN THỊ QUỲNH TRÂM Giới tính: Nữ Ngày, tháng, năm sinh: 08/10/1984 Nơi sinh: T-T-Huế Quê quán: Xuân Thủy, Lệ Thủy, Quảng Bình Dân tộc: Kinh Chỗ ở riêng hoặc địa chỉ liên lạc: 41A Chu Văn An, Hiệp Phú, Q.HCM Điện thoại: 016 587 787 08 Email: lovelytram84@gmail.QUÁ TRÌNH ĐẠO TẠO 1.Đại học: Hệ đào tạo: Chính quy Thời gian đào tạo: từ 09/2003 đến 02/2008 Nơi học ( trường, thành phố): trường Đại học Sư phạm Kỹ thuật tp.HCM Ngành học: Kỹ thuật công nghiệp Tên đồ án, luận án hoặc môn thi tốt nghiệp:Công nghệ gia công gỗ trên máy cưa Ngày & nơi bảo vệ đồ án, luận án hay thi tốt nghiệp: ĐH Sư phạm Kỹ thuật tp.HCM Người hướng dẫn: Th.Thạc sĩ: Hệ đào tạo: Chính quy Thời gian đạo tạo: từ 02/2011 đến 02/2013 Nơi học ( trường, thành phố): trường Đại học Sư phạm Kỹ thuật tp.HCM Ngành học: Công nghệ chế tạo máy i Luan van Tên đồ án, luận án hoặc môn thi tốt nghiệp:Ứng dụng phương pháp Proper Generalized Decomposition và phương pháp phần tử hữu hạn cho bài toán lưu chất. Ngày & nơi bảo vệ đồ án, luận án hay thi tốt nghiệp: trường ĐH Sư phạm Kỹ thuật tp.HCM Người hướng dẫn: TS. Phan Đức Huynh Ngày 18 tháng 09 năm 2013 Người khai ký tên ii Luan van LỜI CAM ĐOAN Tôi cam đoan đây là công trình nghiên cứu của tôi. Các số liệu, kết quả nêu trong luận văn là trung thực và chưa từng được ai công bố trong bất kỳ công trình nào khác.
Hồ Chí Minh, ngày 18 tháng 09 năm 2013 NGUYỄN THỊ QUỲNH TRÂM iii Luan van CẢM TẠ Trong suốt quá trình nghiên cứu đề tài dưới sự hướng dẫn cua thấy Phan Đức Huynh, tôi đã nhận được sự hướng dẫn chu đáo từ phía thầy và đặc biệt nhất là sự quan tâm tận tình vô cùng chân quý từ phía anh Lê Quốc Cường thông qua sự giới thiệu của thầy Phan Đức Huynh, hiện đang theo nghiên cứu bậc Tiến sĩ tại trường Đai học Sư phạm Kỹ thuật. Tôi xin được chuyển đến những dòng biết ơn chân thành và lòng kính trọng sâu sắc đến thầy Phan Đức Huynh và anh Lê Quốc Cường. Đồng thời, tôi cũng xin gửi lời cảm ơn đến ngôi trường ĐH Sư phạm Kỹ thuật mà tôi đã gắn bó trong suốt quãng đời sinh viên và học viên, cùng với đội ngũ các giảng viên-giáo viên-nhân viên của Trường mà tôi đã được theo học và hỗ trợ nhiệt tình,cuối cùng gửi lời thân thương đến những người anh chị,người bạn mà tôi đã được quen biết, trao đổi và giúp đỡ. iv Luan van ABSTRACT Nowadays, numerical techniques become the effective tools to solve the problems in science and engineering.
Eventhough the impressive recent progresses attained in computer technologies and computational simulation techniques, numerous models intractable when the usual and well-experienced discretization techniques are applied for their numerical simulation due to their high complexity and requirements. One of the typical difficulties is highly multi-dimensional models arising from quantum mechanics or kinetic theory descriptions of solids and complex fluids,… When one applies standard mesh based discretization techniques the number of degrees of freedom involved scales exponentially with the dimension of the space concerned. In order to overcome the drawbacks above, one lastest technique in recent years proposed to support, activate in using the mesh-based discretization techniques -FEM -is called Proper Generalized Decomposition (PGD). This is a powerful model reduction technique by means of successive enrichment a separated representation of the unknown field, so the computational complexity of the PGD scales linearly with the dimension of the space.
And a coupling Proper Generalized Decomposition and Finite Element Method – PGD-FEM briefly – will open a new approach in searching a powerful kind of simulation techinique in both terms of computing time and accuracy. Therefore, the topic “ Coupling Proper Generalized Decomposition and Finite Element Method for fluid problem” was born here. Eventhough the topic just started to invest PGD-FEM for fluid problem in a small term of the pressure Poisson equation from 2D unsteady imcompressive Navier-Stokes flow, the comparative results speaked out the outstanding innovative property of PGD-FEM in both computing time and accuracy from the traditional v Luan van discretization technique (FEM). Moreover, in order to overcome its remaining drawbacks and enlarge, develop further research trends, I also provided to solve unsteady imcompressive Navier-Stokes equations by FEM based on the Chorin- Temam projection method.
vi Luan van TÓM TẮT Ngày nay phương pháp số là một công cụ đắc lực giúp giải quyết hầu hết các bài toán trong khoa học và kỹ thuật. Mặc dù với những tiến bộ, phát triển vượt bậc đạt được trong công nghệ máy tính, kỹ thuật t nh toán nhưng những khó khăn để giải quyết nhiều bài toán vẫn còn bị thách thức khi mà phương pháp rời rạc truyền thống đang bị hạn chế do tính phức tạp và mức độ yêu cầu đòi hỏi ngày càng cao của bài toán. Có thể nêu một trong những khó khăn điển hình, nổi cộm là bài toán có số chiều không gian lớn thường gặp trong cơ lượng tử, thuyết động học của cơ lưu chất phức tạp,…Khi sử dụng phương pháp rời rạc thông thường thì độ phức tạp của bài toán tăng theo tỉ lệ hàm mũ với số chiều không gian của bài toán. Để nhằm khắc phục tính hạn chế trên, một phương pháp rất mới đã ra đời trong vài năm gần đây đã góp phần bổ trợ, thúc đẩy trong quá trình phối hợp với phương pháp rời rạc, cụ thể là phương pháp phần tử hữu hạn (FEM) mà được nghiên cứu ở đây, với tên gọi là phương pháp Proper Generalized Decomposition (PGD).
Đây là một công cụ giảm bậc mô hình bài toán dựa trên cơ sở tách biến giúp độ phức tạp của bài toán giảm xuống với tỉ lệ tuyến tính theo số chiều của bài toán. Vì thế sự kết hợp giữa phương pháp PGD và FEM (gọi tắt PGD-FEM) sẽ bước đầu mở ra một hướng tiếp cận mới trong việc tìm kiếm một loại hình phương pháp số mới với tính năng ưu việt hơn về mặt thời gian xử lí mà vẫn đảm bảo độ chính xác so với phương pháp rời rạc truyền thống đã có. Và “ứng dụng phương pháp PGD-FEM cho bài toán lưu chất” đã ra đời trong đề tài nghiên cứu ở đây. Mặc dù đề tài chỉ mới bước đầu khái thác phương pháp PGD-FEM cho lĩnh vực bài toán lưu chất ở một khía cạnh hẹp là giải quyết phương trình Poisson áp suất 2D cho bài toán Navier-Stokes của dòng chảy nhớt không nén phụ thuộc vào thời gian trong hai trường hợp điều kiên biên Dirichlet đồng nhất và điều kiện biên vii Luan van hỗn hợp ( Dirchlet-Neumann) nhưng những kết quả đạt được đã cho thấy sự ưu việt khi giải quyết bằng phương pháp PGD-FEM về mặt thời gian t nh toán và độ chính xác so với phương pháp rời rạc truyền thống (FEM).
Đồng thời với mong muốn tạo một sự thuận lợi trong việc hoàn thiện cũng như mở rộng, phát triển hơn cho đề tài trong tương lại, tác giả cũng đã đề cập đến việc giải phương trình Navier-Stokes cho dòng chảy nhớt không nén phụ thuộc thời gian với điều kiện biên lid-driven cavity bằng phương pháp FEM dựa trên kỹ thuật tham chiếu Chorin-Temam. viii Luan van MỤC LỤC TRANG Trang tựa Quyết định giao đề tài Lý lịch cá nhân. i Lời cam đoan. iii Cảm tạ.
ix Danh mục kí hiệu-từ viết tắt. xi Danh mục hình vẽ. xii Chƣơng 1.1 Tổng quan về hướng nghiên cứu.2 Mục đ ch nghiên cứu, khách thể và đối tượng nghiên cứu .3 Xác định nhiệm vụ và phạm vi nghiên cứu của đề tài .4 Phương pháp nghiên cứu.1 Phương pháp Proper Generalized Decomposition (PGD) .2 Phương pháp phần tử hữu hạn (FEM). 11 ix Luan van Chƣơng 3.
ỨNG DỤNG PHƢƠNG PHÁP PGD và FEM CHO BÀI TOÁN LƢU CHẤT.1 Giới thiệu phương trình Navier-Stokes .2 Giải phương trình Poisson bằng phương pháp PGD-FEM.1 Trường hợp điều kiên biên Dirichlet đồng nhất.1 Tiến trình giải bài toán bằng phương pháp PGD-FEM .2 Sơ đồ giải thuật tổng quát .3 Kết quả - nhận xét .2 Trường hợp điều kiên biên hỗn hợp.1 Tiến trình giải bài toán bằng phương pháp PGD-FEM .2 Kết quả - nhận xét .3 Phƣơng trình Navier-Stokes không nén phụ thuộc vào thời gian .1 Mô hình bài toán .2 Điều kiên biên của bài toán .3 Tiến trình các bước giải bằng phương pháp FEM.4 Sơ đồ giải thuật tổng quát .5 Kết quả - nhận xét. KẾT LUẬN VÀ HƢỚNG PHÁT TRIỂN. 48 TÀI LIỆU THAM KHẢO. 51 x Luan van DANH MỤC KÍ HIỆU, TỪ VIẾT TẮT <.> L2(Ω) Tích trong của hai hàm trên L2 trong miền Ω ||.
||2 Chuẩn vec-tơ trong L2 hay chuẩn Euclide Resn Phần sai số giữa nghiệm chính xác và nghiệm xấp xỉ Ω Miền khảo sát của bài toán Ωx Miền khảo sát theo phương x Ωy Miền khảo sát theo phương y ∂Ω Biên của miền Ω ∂Ω D Biên Dirichlet của ∂Ω ∂Ω N Biên Neumann của ∂Ω H1 Không gian hàm Sobolev mà có giá trị triệt tiêu trên ∂Ω D n Vec-tơ pháp tuyến (hướng ra ngoài) trên biên ∂Ω N V Vec-tơ hàm dạng tại các nút trên miền Ω M Vec-tơ hàm dạng tại các nút trên miền Ωx N Vec-tơ hàm dạng tại các nút trên miền Ωy Nnod Tổng số nút trên miền Ω Nnod_x Tổng số nút trên miền Ωx Nnod_y_ Tổng số nút trên miền Ωy X,R,F Hàm phụ thuộc trên miền Ωx Y,S,G Hàm phụ thuộc trên miền Ωy X,R,F Vec-tơ giá trị hàm X,R,F tại các nút trên miền Ωx xi Luan van Y,S,G Vec-tơ giá trị hàm Y,S,G tại các nút trên miền Ωy p p p Toán tử gradient , x y u u .u Toán tử divergence x y 2 u 2u 2u Toán tử Laplace 2 x y 2 FEM Finite Element Method PGD Proper Generalized Decoposition PGD-FEM coupling Proper Generalized Decomposition and Finite Element Method MBS Multi-Bead Spring ROM Reduced-Order Model LATIN LArger Time INcremential POD Proper Orthogonal Decomposition SVD Singular Value Decomposition PDE Partial Differential Equations xii Luan van DANH MỤC HÌNH VẼ Hình 2.1 Miền khảo sát và điều kiện biên của phương trình Poisson Hình 3.