Luận văn thạc sĩ về mô phỏng dòng chảy lưu chất nhớt không nén qua trụ tròn

Luận văn thạc sĩ nghiên cứu mô phỏng dòng chảy lưu chất nhớt không nén qua trụ tròn bằng phương pháp biên nhúng, cung cấp kiến thức chuyên sâu.

• Tác giả: Võ Trung Nam
• Người hướng dẫn: Ts. Phan Đức Huynh
• Trường học: Đại học Sư Phạm Kỹ Thuật TP.Hồ Chí Minh
• Chuyên ngành: Cơ Khí Chế Tạo Máy, Xây Dựng và Cơ Học Ứng Dụng
• Đề tài: Mô phỏng dòng chảy lưu chất nhớt qua trụ tròn bằng phương pháp biên nhúng
• Loại tài liệu: luận văn tốt nghiệp
• Năm xuất bản: 2013
• Địa điểm: Tp. Hồ Chí Minh

Mục lục chi tiết

• LỜI CAM ĐOAN
• CẢM TẠ
• TÓM TẮT
• ABSTRACT
• MỤC LỤC
• 1. CHƯƠNG 1: TỔNG QUAN
• 1.1. Tổng quan chung về lĩnh vực nghiên cứu
• 1.2. Mục đích của đề tài
• 1.3. Nhiệm vụ của đề tài và giới hạn của đề tài
• 1.4. Phương pháp nghiên cứu
• 1.5. Kế hoạch thực hiện
• 2. CHƯƠNG 2: PHƯƠNG PHÁP BIÊN NHÚNG CHO PHƯƠNG TRÌNH NAVIER-STOKES VỚI BIÊN CỨNG
• 2.1. Phương trình chuyển động
• 2.2. Phương pháp số
• 2.2.1. Sự rời rạc hóa theo không gian và thời gian
• 2.2.2. Phương pháp giải
• 2.2.3. Giải phương trình Navier-Stokes
• 2.2.3.1. Xử lý các thành phần phi tuyến, độ nhớt và thành phần lực khối
• 2.2.3.2. Sự hiệu chỉnh áp suất
• 2.2.3.3. Lưới xen kẽ
• 2.2.3.3.1. Đạo hàm xấp xỉ
• 2.2.3.3.2. Các điều kiện biên
• 2.2.3.3.3. Phương trình Poisson
• 3. CHƯƠNG 3: CẤU TRÚC HÀM DIRAC DELTA
• 4. CHƯƠNG 4: KẾT QUẢ MÔ PHỎNG SỐ
• 4.1. Dòng chảy qua trụ tròn cố định
• 4.2. Dòng chảy qua trụ tròn dao động
• 5. CHƯƠNG 5: KẾT LUẬN
• TÀI LIỆU THAM KHẢO

Tóm tắt

I. Tổng quan về mô phỏng dòng chảy lưu chất nhớt qua trụ tròn

Mô phỏng dòng chảy lưu chất nhớt qua trụ tròn là một lĩnh vực nghiên cứu quan trọng trong cơ học chất lỏng. Việc hiểu rõ về dòng chảy này không chỉ giúp tối ưu hóa thiết kế các cấu trúc mà còn hỗ trợ trong nhiều ứng dụng thực tiễn như trong ngành công nghiệp và y học. Phương pháp biên nhúng (IBM) đã được phát triển để giải quyết các bài toán phức tạp liên quan đến tương tác giữa lưu chất và kết cấu. Phương pháp này cho phép mô phỏng chính xác mà không cần phải chia lưới lại khi biên di chuyển.

1.1. Khái niệm về dòng chảy lưu chất nhớt

Dòng chảy lưu chất nhớt là hiện tượng mà trong đó lưu chất có độ nhớt cao, ảnh hưởng đến cách mà nó di chuyển qua các bề mặt. Đặc điểm này rất quan trọng trong việc mô phỏng và tính toán các ứng dụng thực tiễn.

1.2. Tầm quan trọng của mô phỏng trong nghiên cứu

Mô phỏng dòng chảy giúp các nhà nghiên cứu hiểu rõ hơn về hành vi của lưu chất trong các điều kiện khác nhau. Điều này có thể dẫn đến những cải tiến trong thiết kế và hiệu suất của các sản phẩm công nghiệp.

II. Vấn đề và thách thức trong mô phỏng dòng chảy lưu chất nhớt

Mô phỏng dòng chảy lưu chất nhớt qua trụ tròn gặp nhiều thách thức, đặc biệt là trong việc xử lý các hình dạng phức tạp và biên di chuyển. Các phương pháp truyền thống thường yêu cầu chia lưới lại mỗi khi biên di chuyển, điều này không chỉ tốn thời gian mà còn ảnh hưởng đến độ chính xác của kết quả. Phương pháp biên nhúng đã được đề xuất như một giải pháp hiệu quả để giải quyết những vấn đề này.

2.1. Các vấn đề trong việc chia lưới

Chia lưới là một bước quan trọng trong mô phỏng, nhưng việc này trở nên phức tạp khi biên dạng của vật thể thay đổi. Phương pháp biên nhúng giúp giảm thiểu sự phức tạp này.

2.2. Tác động của biên di chuyển đến mô phỏng

Biên di chuyển có thể gây ra sự thay đổi lớn trong dòng chảy, dẫn đến việc cần phải điều chỉnh lưới thường xuyên. Phương pháp biên nhúng cho phép xử lý vấn đề này một cách hiệu quả.

III. Phương pháp biên nhúng trong mô phỏng dòng chảy

Phương pháp biên nhúng (IBM) là một kỹ thuật mạnh mẽ cho phép mô phỏng dòng chảy lưu chất mà không cần phải chia lưới lại khi biên di chuyển. Phương pháp này sử dụng hàm Dirac delta để mô phỏng sự tương tác giữa lưu chất và biên, giúp giảm thiểu độ phức tạp trong tính toán. IBM đã chứng minh được tính hiệu quả và độ chính xác trong nhiều nghiên cứu trước đây.

3.1. Nguyên lý hoạt động của phương pháp biên nhúng

Phương pháp biên nhúng hoạt động bằng cách thay thế sự hiện diện của bề mặt bằng một lực khối, cho phép mô phỏng mà không cần thay đổi lưới tính toán.

3.2. Ưu điểm của phương pháp biên nhúng

Phương pháp này giúp giảm thiểu thời gian tính toán và tăng độ chính xác trong mô phỏng, đặc biệt là trong các bài toán có hình dạng phức tạp.

IV. Ứng dụng thực tiễn của mô phỏng dòng chảy lưu chất nhớt

Mô phỏng dòng chảy lưu chất nhớt qua trụ tròn có nhiều ứng dụng thực tiễn trong các lĩnh vực như thiết kế máy móc, y học và nghiên cứu môi trường. Các kết quả từ mô phỏng có thể được sử dụng để tối ưu hóa thiết kế và cải thiện hiệu suất của các sản phẩm. Nghiên cứu này không chỉ có giá trị lý thuyết mà còn mang lại lợi ích thực tiễn lớn.

4.1. Ứng dụng trong ngành công nghiệp

Trong ngành công nghiệp, mô phỏng dòng chảy giúp tối ưu hóa thiết kế của các thiết bị như tua bin, van và ống dẫn, từ đó nâng cao hiệu suất hoạt động.

4.2. Ứng dụng trong y học

Mô phỏng dòng chảy lưu chất nhớt cũng được áp dụng trong y học, đặc biệt là trong việc mô phỏng tuần hoàn máu và thiết kế các thiết bị y tế như tim nhân tạo.

V. Kết luận và tương lai của mô phỏng dòng chảy lưu chất nhớt

Mô phỏng dòng chảy lưu chất nhớt qua trụ tròn bằng phương pháp biên nhúng đã chứng minh được tính hiệu quả và độ chính xác trong nhiều nghiên cứu. Tương lai của lĩnh vực này hứa hẹn sẽ có nhiều tiến bộ hơn nữa, đặc biệt là trong việc phát triển các phương pháp mới và cải tiến các phương pháp hiện tại. Việc nghiên cứu sâu hơn về tương tác giữa lưu chất và kết cấu sẽ mở ra nhiều cơ hội mới trong các ứng dụng thực tiễn.

5.1. Triển vọng nghiên cứu trong tương lai

Nghiên cứu trong lĩnh vực mô phỏng dòng chảy sẽ tiếp tục phát triển, với nhiều ứng dụng mới và cải tiến trong các phương pháp hiện tại.

5.2. Tác động của công nghệ mới đến mô phỏng

Sự phát triển của công nghệ tính toán và phần mềm sẽ giúp nâng cao khả năng mô phỏng, cho phép giải quyết các bài toán phức tạp hơn trong tương lai.

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

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 Tp. Hồ Chí Minh, ngày 28 tháng 02 năm 2013 (Ký tên và ghi rõ họ tên) iii CẢM TẠ Tác giả xin chân thành cảm ơn quí Thầy Cô trong khoa Cơ Khí Chế Tạo Máy, khoa Xây Dựng và Cơ Học Ứng Dụng trường Đại học Sư Phạm Kỹ Thuật TP.Hồ Chí Minh đã tận tình giúp đỡ, hướng dẫn và tạo điều kiện thuận lợi để tác giả hoàn thành luận văn tốt nghiệp đúng tiến độ. Đặc biệt, tác giả xin chân thành cảm ơn TS. Phan Đức Huynh, dù rất bận rộn với công việc giảng dạy nhưng Thầy vẫn luôn dành thời gian quan tâm, hướng dẫn, chỉ bảo tận tình cho tôi trong suốt quá trình nghiên cứu, thực hiện luận văn. Tác giả cũng chân thành cám ơn NCS. Lê Quốc Cường, Ths. Nguyễn Hoàng Sơn đã nhiệt tình góp ý, giúp đỡ tác giả trong suốt quá trình nghiên cứu thực hiện luận văn. Hồ Chí Minh, ngày 28 tháng 02 năm 2013 Võ Trung Nam iv TÓM TẮT Sự tương tác giữa lưu chất – kết cấu hiện nay đang là một chủ đề nổi cộm và được nghiên cứu rất nhiều. Sư tương tác này làm cho biên di chuyển dẫn đến sự khó khăn trong việc chia lưới, tính toán, mô phỏng. Ngoài ra với những kết cấu có biên dạng phức tạp cũng xảy ra khó khăn như trên. Phương pháp phần tử hữu hạn được xem là giải pháp số toàn diện và thông dụng nhất tại thời điểm hiện nay để giải quyết những khó khăn đó. Đứng ở góc độ, khía cạnh khác thì phương pháp biên nhúng được đề nghị để giải quyết các bài toán này. Đối với phương pháp biên nhúng thì việc chia lưới được thực hiện dễ dàng và không phụ thuộc vào biên dạng hình học của vật thể cho dù vật thể đó có biên dạng hình học phức tạp. Ngoài ra vấn đề biên di chuyển cũng được giải quyết, các phương pháp số khác sẽ tái tạo lưới mỗi khi biên di chuyển ở mỗi bước thời gian và vấn đề này không cần thiết khi ta sử dụng phương pháp biên nhúng. Thông qua kết quả từ bài toán mô phỏng dòng chảy lưu chất nhớt không nén được qua trụ tròn đã chứng minh tính hợp lý, tính chính xác, tính hiệu quả, độ ổn định của lời giải so với các phương pháp số khác và so với phương pháp thực nghiệm. v ABSTRACT Fluid – Structure interaction is the hot topic at the moment and a lot of research. This interaction made the boundary move and lead to difficulties in meshing create, calculation and simulation. In addition, with structure has complex – shaped also occurs difficult above. Finite element method is considered a comprehensive solution and the most common at the present time to solve these difficulties.The other standpoint, immersed boundary method is recommended to solve this problem. For immersed boundary method, the mesh is made easy and does not depend on the boundary geometry of the body even if the body boundary complex geometry.In addition, moving boundary also well resolved, other methods of renewable mesh whenever boundary moves in each time step, and this is unnecessary when using the immersed boundary method. Through the results from simulations of incompressible viscous fluid flow past a circular cylinder has demonstrated the validity, accuracy, efficiency and stability of the solution compared with other numerical methods and experiment results. vi MỤC LỤC Trang tựa Trang Quyết định giao đề tài Lý lịch khoa học . i Lời cam đoan . ii Cảm tạ . iv Mục lục. vi Ký hiệu khoa học . viii Danh sách các hình. ix Danh sách các bảng . x Chƣơng 1 TỔNG QUAN . 1 Chƣơng 2 PHƢƠNG PHÁP BIÊN NHÚNG CHO PHƢƠNG TRÌNH NAVIER - STOKES VỚI BIÊN CỨNG .1 Phương trình chuyển động .2 Phương pháp số .1 Sự rời rạc hóa theo không gian và thời gian .2Phương pháp giải.3 Giải phương trình Navier-Stokes .1 Xử lý các thành phần phi tuyến, độ nhớt và thành phần lực khối .2 Sự hiệu chỉnh áp suất .3 Lưới xen kẽ .1 Đạo hàm xấp xỉ .2 Các điều kiện biên .3 Phương trình Poisson . 17 Chƣơng 3CẤU TRÚC HÀM DIRAC DELTA . 20 Chƣơng 4KẾT QUẢ MÔ PHỎNG SỐ .1 Dòng chảy qua trụ tròn cố định .2 Dòng chảy qua trụ tròn dao động . 39 Chƣơng 5 KẾT LUẬN . 44 TÀI LIỆU THAM KHẢO . 47 viii KÝ HIỆU KHOA HỌC  Xs, t    X s, t , Y s, t    X k , Yk  là hàm vector được cho bởi tọa độ của các điểm trên biên Γ(như là một hàm của độ dài cung s và thời gian t). k=0,1,2,…,m-1  F  Fx s, t , Fy s, t  là lực khối tác dụng trên biên  Us, t   U s, t ,V s, t   U k ,Vk  là vận tốc của điểm lưới Lagrangian  f   f x x, t , f y x, t là lực vật thể được lồng vào phương trình Navier-Stokes  x  x, y  là tọa độ theo lưới Eulerian  ux, t   ux, t , vx, t  là vận tốc của lưu chất ( theo 2 chiều x, y)  px, t  là áp suất của lưu chất.   là khối lượng riêng của lưu chất   là độ nhớt  u* là vận tốc trung gian  ∇p là gradient áp suất. 2 2    là toán tử Laplace x 2 y 2     grap     ,   x y      (.) xem như là một ký hiệu và được dùng : (  .     x y  (với grap, div là các hàm được sử dụng trong toán tử Laplace)  Lb là chiều dài của đường cong khép kín Γ ix  𝛿 𝐱 − 𝐗 𝑠, 𝑡 = 𝛿 𝑥 − 𝑋 𝛿 y − Y là hàm Dirac delta. DANH SÁCH CÁC HÌNH HÌNH TRANG Hình 2.1 a) Hệ lưu chất-kết cấu đơn giản. b) Lưới rời rạc Euler (đánh dấu sáng) và lưới Lagrange (đánh dấu tối) .1 c) Hình vẽ phác họa việc tính toán thành phần lực khối sử dụng hàm Dirac delta .2 Lưới xen kẽ với các ô biên .3 Sự dao động trong bài toán đối lưu-khuếch tán khi sử dụng sai phân trung tâm .1 Hàm Dirac delta .1 Đường dòng cho hệ số Re=20 và Re=40 .2 Trường áp suất cho hệ số Re=20 và Re=40 .3 Đường bao xoáy tại Re=100 và Re=200 .4 Hệ số cản tại Re=100 và Re=200 .5 Hệ số nâng tại Re=100 và Re=200 .6 Đường dòng tại Re=100 và Re=200 .7 Trường áp suất tại Re=100 và Re=200 .8 Đường bao xoáy và trường áp suất tại Re=300 .9 Hệ số cản và hệ số nâng tại Re=300 .9 Đường bao xoáy của trụ tròn dao dộng sau 40 chu kỳ dao động tại Re=200, (a: Zhang & Zheng; b: Hiện tại) .10 Hệ số cảnCDtại Re=200 đối với trụ tròn cố định và trụ tròn dao động .11 Hệ số nâng CLtại Re=200 đối với trụ tròn cố định và trụ tròn dao động 42 x Hình 4.12 Trường áp suất và đường dòng tại Re=200 cho trường hợp trụ tròn dao động . 43 DANH SÁCH CÁC BẢNG BẢNG TRANG Bảng 1.1 Kế hoạch thực hiện luận văn .1 Chiều dài vùng tuần hoàn (L/d), hệ số cản (CD) cho Re=20,Re=40 .2 Hệ số cản tại Re=100 và Re=200 .3 Hệ số nâng tại Re=100 và Re=200 .4 Hệ số Strouhal tại Re=80, 100, 200, 300 . 38 xi Chƣơng 1 TỔNG QUAN 1.1 Tổng quan chung về lĩnh vực nghiên cứu Trong tính toán động lực học lƣu chất, các vấn đề cần quan tâm nhất đó là sự chính xác, hiệu suất tính toán, độ ổn định của lời giải và đặc biệt là xử lý đƣợc các dạng hình học phức tạp. Có rất nhiều phƣơng pháp cho việc giải các bài toán dòng không nén đƣợc trong miền hình học phức tạp. Trong một số ứng dụng, phƣơng pháp phần tử hữu hạn (FEM) trên lƣới phi cấu trúc thì đƣợc sử dụng rộng rãi nhất hiện nay. Sự chính xác của FEM có thể đƣợc cải thiện bằng cách sử dụng các hàm nội suy bậc cao. Tuy nhiên, tạo ra lƣới phi cấu trúc theo yêu cầu bằng FEM đòi hỏi một kỹ thuật cao và khá tốn kém. Hơn nữa, để giải các bài toán với một biên chuyển động mà đặc biệt là trong lĩnh vực tƣơng tác giữa kết cấu và lƣu chất thật không dễ thực hiện, trƣờng hợp này ta phải chia lƣới lại sau mỗi lần di chuyển của biên. Đã có sự tiến bộ vuợt bậc đáng kể của các phƣơng pháp trong việc tính toán chính xác và hiệu quả đối với những vật thể có hình dạng phức tạp bất kỳ. Phƣơng pháp biên nhúng (Immersed Boundary Methods - IBM) đƣợc đề xuất bởi C. Peskin [1,2] gần đây đã đƣợc đƣa ra để giải quyết đối với các bài toán có những dạng hình học phức tạp trong khi yêu cầu tính toán ít hơn, tạo lƣới dễ hơn so với các phƣơng pháp khác mà vẫn đảm bảo đƣợc sự chính xác. Phƣơng pháp này quy định một lực khối để thay thế sự hiện diện của một bề mặt mà không làm thay đổi lƣới tính toán. Đặc biệt đối với trƣờng hợp biên di chuyển, cụ thể hơn là trong sự tƣơng tác giữa lƣu chất và kết cấu làm biên dịch chuyển. Lƣới bám theo vật thể phải đƣợc chia lại ở mỗi bƣớc thời gian khi tính toán, việc chia lƣới ảnh hƣởng trực tiếp đến chi phí tính toán, độ chính xác và sự ổn định của lời giải. Có thể nói đó là vấn đề cần xem xét một cách nghiêm túc. Lƣới Đềcác đƣợc xây dựng dễ dàng, nhanh chóng mà không 1 cần giải thuật phức tạp và đặc biệt là không bị ảnh hƣởng đáng kể đến sự phức tạp của hình dạng vật thể. Đây đƣợc xem là ƣu điểm nổi bật bậc nhất của phƣơng pháp biên nhúng. Vì vậy trong bài toán tƣơng tác giữa lƣu chất và kết cấu thì có thể nóiphƣơng pháp biên nhúng là phƣơng pháp số tốt nhất tại thời điểm hiện nay. Các kết quả nghiên cứu trong nƣớc và ngoài nƣớc đã công bố IBM lần đầu tiên đƣợc giới thiệu bởi C. Peskin (1972) với việc mô phỏng sự tƣơng tác giữa các cơ và dòng máu chảy trong lúc tim đang đập và nhìn chung kết quả phù hợp cho dòng lƣu chất với biên đàn hồi đƣợc nhúng. Trong phƣơng pháp này dòng lƣu chất đƣợc điều khiển bởi phƣơng trình Navier-Stokes không nén đƣợc và chúng đƣợc xử lý trên lƣới Đềcác cố định. Một số công trình nghiên cứu điển hình về IBM đã công bố trên các tạp chí khoa học uy tín trên thế giới: C. Peskin, The immersed boundary method, Acta Numerica.