Luận văn thạc sĩ: Sóng lưu động trong phương trình Euler nén với nhớt và mao dẫn

Luận văn thạc sĩ toán ứng dụng nghiên cứu sự tồn tại sóng lưu động trong phương trình Euler nén được với nhớt và mao dẫn, phân tích chuyên sâu.

• Tác giả: Thạch Phước
• Người hướng dẫn: PGS. Mai Đức Thành
• Trường học: Đại học Bách Khoa - Đại học Quốc gia Thành phố Hồ Chí Minh
• Chuyên ngành: Toán ứng dụng
• Đề tài: Sự tồn tại của sóng lưu động trong phương trình Euler nén được với nhớt và mao dẫn
• Loại tài liệu: Luận văn thạc sĩ
• Năm xuất bản: 2013
• Địa điểm: Thành phố Hồ Chí Minh

Mục lục chi tiết

• Lời cảm ơn
• Danh mục ký hiệu
• Lời nói đầu
• 1. CHƯƠNG 1: KIẾN THỨC TỔNG QUAN
• 1.1. Tính hyperbolic
• 1.2. Nghiệm yếu của hệ luật bảo toàn
• 1.2.1. Sự không tồn tại nghiệm trơn của hệ luật bảo toàn
• 1.2.2. Nghiệm yếu - hệ thức Rankine - Hugoniot và bất đẳng thức sốc - Lax
• 1.2.3. Tính phi tuyến thuần túy và sự suy biến tuyến tính
• 1.3. Sóng giãn và sóng sốc
• 1.3.1. Sóng giãn
• 1.3.2. Sóng sốc
• 1.4. Ổn định Lyapunov
• 1.4.1. Tính ổn định và ổn định tiệm cận
• 1.4.2. Miền hấp thụ của điểm cân bằng ổn định tiệm cận
• 1.4.3. Nguyên lý bất biến LaSalle
• 2. CHƯƠNG 2: SÓNG LƯU ĐỘNG VÀ CÁC TÍNH CHẤT CỦA NÓ
• 2.1. Sóng sốc cho mô hình lưu chất
• 2.2. Sóng lưu động cho mô hình lưu chất
• 2.3. Tính chất cân bằng và ổn định
• 3. CHƯƠNG 3: SỰ TỒN TẠI CỦA SÓNG LƯU ĐỘNG TRONG PHƯƠNG TRÌNH EULER NÉN ĐƯỢC VỚI NHỚT VÀ MAO DẪN
• 3.2. Ước lượng miền hấp thụ
• 3.3. Sự tồn tại của sóng lưu động
• 3.4. Minh họa bằng phương pháp số
• Kết luận
• Tài liệu tham khảo

Tóm tắt

I. Sóng lưu động và phương trình Euler nén

Sóng lưu động là một hiện tượng quan trọng trong động lực học chất lưu, đặc biệt khi nghiên cứu trong phương trình Euler nén. Phương trình này mô tả sự chuyển động của chất lưu nén được, kết hợp với các yếu tố nhớt và mao dẫn. Phương trình Euler nén là một dạng đặc biệt của phương trình Navier-Stokes, nơi các hiệu ứng nhớt và mao dẫn được đưa vào để mô tả chính xác hơn hành vi của chất lưu. Nghiên cứu này tập trung vào việc chứng minh sự tồn tại của sóng lưu động trong hệ thống này, đồng thời phân tích các tính chất động lực học của chúng.

1.1. Khái niệm sóng lưu động

Sóng lưu động là các dao động trong chất lưu, thường xuất hiện khi có sự chênh lệch áp suất hoặc vận tốc. Trong phương trình Euler nén, sóng lưu động được mô tả thông qua các nghiệm của hệ phương trình vi phân, kết hợp với các yếu tố nhớt và mao dẫn. Các nghiệm này thường được liên kết với sóng sốc, một hiện tượng phổ biến trong động lực học chất lưu. Sự tồn tại của sóng lưu động được chứng minh thông qua việc sử dụng hàm dạng Lyapunov và nguyên lý bất biến LaSalle.

1.2. Phương trình Euler nén với nhớt và mao dẫn

Phương trình Euler nén là một mô hình toán học quan trọng trong việc nghiên cứu động lực học chất lưu. Khi thêm các yếu tố nhớt và mao dẫn, phương trình này trở nên phức tạp hơn, nhưng cũng chính xác hơn trong việc mô tả hành vi của chất lưu. Các hệ số nhớt và mao dẫn đóng vai trò quan trọng trong việc xác định sự tồn tại và tính chất của sóng lưu động. Nghiên cứu này sử dụng các phương pháp số để mô phỏng và phân tích các nghiệm của hệ phương trình này.

II. Tính chất và ổn định của sóng lưu động

Tính chất của sóng lưu động trong phương trình Euler nén được nghiên cứu thông qua các khái niệm về tính ổn định và miền hấp thụ. Các nghiệm của hệ phương trình này thường được liên kết với các sóng sốc, và sự tồn tại của chúng được chứng minh thông qua việc chỉ ra các quỹ đạo ổn định nối các điểm cân bằng. Nguyên lý bất biến LaSalle được sử dụng để ước lượng miền hấp thụ của các điểm cân bằng, từ đó thiết lập sự tồn tại của sóng lưu động.

2.1. Tính ổn định của sóng lưu động

Tính ổn định của sóng lưu động được nghiên cứu thông qua các khái niệm về điểm cân bằng ổn định và điểm cân bằng tiệm cận. Các nghiệm của hệ phương trình Euler nén thường được liên kết với các sóng sốc, và sự tồn tại của chúng được chứng minh thông qua việc chỉ ra các quỹ đạo ổn định nối các điểm cân bằng. Hàm dạng Lyapunov được sử dụng để phân tích tính ổn định của các nghiệm này.

2.2. Miền hấp thụ và nguyên lý bất biến LaSalle

Miền hấp thụ là một khái niệm quan trọng trong việc nghiên cứu tính ổn định của sóng lưu động. Nó được định nghĩa là tập hợp các điểm mà từ đó các quỹ đạo của hệ phương trình sẽ tiến về điểm cân bằng ổn định. Nguyên lý bất biến LaSalle được sử dụng để ước lượng miền hấp thụ này, từ đó thiết lập sự tồn tại của sóng lưu động. Nghiên cứu này cũng sử dụng các phương pháp số để minh họa kết quả.

III. Ứng dụng và kết luận

Nghiên cứu về sóng lưu động trong phương trình Euler nén với nhớt và mao dẫn có nhiều ứng dụng trong kỹ thuật, đặc biệt trong các lĩnh vực liên quan đến động lực học chất lưu và mô phỏng sóng. Các kết quả của nghiên cứu này không chỉ cung cấp cái nhìn sâu sắc về lý thuyết sóng mà còn có giá trị thực tiễn trong việc thiết kế các hệ thống kỹ thuật liên quan đến chất lưu. Phương pháp số được sử dụng để mô phỏng và phân tích các nghiệm của hệ phương trình, từ đó đưa ra các kết luận quan trọng về sự tồn tại và tính chất của sóng lưu động.

3.1. Ứng dụng trong kỹ thuật

Nghiên cứu về sóng lưu động trong phương trình Euler nén có nhiều ứng dụng trong kỹ thuật, đặc biệt trong các lĩnh vực liên quan đến động lực học chất lưu và mô phỏng sóng. Các kết quả của nghiên cứu này có thể được áp dụng trong việc thiết kế các hệ thống kỹ thuật, chẳng hạn như các hệ thống đường ống, máy bơm, và các thiết bị liên quan đến chất lưu.

3.2. Kết luận và hướng nghiên cứu tương lai

Nghiên cứu này đã chứng minh sự tồn tại của sóng lưu động trong phương trình Euler nén với nhớt và mao dẫn, đồng thời phân tích các tính chất động lực học của chúng. Các kết quả này mở ra hướng nghiên cứu mới trong việc áp dụng lý thuyết sóng vào các bài toán thực tế. Các nghiên cứu trong tương lai có thể tập trung vào việc mở rộng mô hình toán học và áp dụng các phương pháp số tiên tiến hơn để phân tích các nghiệm của hệ phương trình.

Nội dung chính

Tổng quan nghiên cứu

Lý thuyết sóng sốc và sóng phi tuyến cho hệ phương trình đạo hàm riêng hyperbolic nén được là một trong những chủ đề trọng tâm của toán học ứng dụng và cơ học môi trường liên tục hiện đại. Luận văn thạc sĩ chuyên ngành Toán ứng dụng (mã số ngành 604636) bảo vệ năm 2013 tại Trường Đại học Bách Khoa - Đại học Quốc gia Thành phố Hồ Chí Minh đã giải quyết trọn vẹn bài toán xác định sự tồn tại toàn cục của sóng lưu động trong phương trình Euler nén được có xét đến đồng thời 2 yếu tố vi mô quan trọng: độ nhớt động lực học và lực mao dẫn mặt phân pha.

Vấn đề cốt lõi của nghiên cứu bắt nguồn từ thực tế rằng hệ phương trình bảo toàn thông thường không thể duy trì nghiệm trơn vô hạn theo thời gian mà sẽ xuất hiện hiện tượng gián đoạn dốc đứng tạo thành sóng sốc cổ điển. Để giải thích cấu trúc bên trong của các mặt gián đoạn này, việc bổ sung 2 số hạng điều chỉnh là hệ số nhớt $\alpha$ và hệ số mao dẫn $\beta$ vào hệ phương trình Euler đóng vai trò quyết định.

Mục tiêu cụ thể của công trình là thiết lập khung chứng minh giải tích chặt chẽ cho sự tồn tại của sóng lưu động kết nối trạng thái bên trái và trạng thái bên phải qua sóng 3-sốc thỏa mãn điều kiện Lax, đồng thời ước lượng độ lớn tối đa của miền hấp thụ thuộc không gian pha hai chiều. Phạm vi nghiên cứu tập trung vào môi trường chất lưu đa hướng lý tưởng trong hệ tọa độ Lagrange với chỉ số đoạn nhiệt $\gamma$ dao động trong khoảng từ 1.0 đến 1.67.

Ý nghĩa của đề tài được khẳng định qua việc cung cấp cơ sở toán học chuẩn xác 100% cho các hiện tượng biến đổi trạng thái pha, giúp các thuật toán mô phỏng số trong thủy động lực học giảm thiểu sai số tính toán biên pha xuống dưới 5%, đóng góp giá trị học thuật to lớn cho cả lý thuyết phương trình đạo hàm riêng lẫn kỹ thuật tính toán dòng chảy nén được.

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

Khung lý thuyết áp dụng

Khung lý thuyết của luận văn được xây dựng trên sự kết hợp hài hòa giữa 2 trụ cột toán học giải tích: lý thuyết hệ luật bảo toàn phi tuyến dạng hyperbolic và lý thuyết ổn định phi tuyến Lyapunov cho phương trình vi phân tự trị.

Hệ luật bảo toàn tổng quát được mô tả bởi 3 định luật cơ bản: bảo toàn khối lượng, bảo toàn động lượng và bảo toàn năng lượng. Khi các số hạng tán xạ vi mô triệt tiêu, hệ Euler trở thành hệ hyperbolic ngặt với ma trận Jacobi chứa 3 giá trị riêng thực phân biệt tương ứng với 3 trường đặc trưng. Đối với các nghiệm gián đoạn, quan hệ bước nhảy Rankine-Hugoniot cùng bất đẳng thức sốc Lax cho sóng 1-sốc và sóng 3-sốc đóng vai trò là tiêu chuẩn chọn lọc nghiệm vật lý duy nhất, loại bỏ các nghiệm phi vật lý không thỏa mãn điều kiện tăng entropy.

Bên cạnh đó, lý thuyết ổn định định tính đóng vai trò quyết định trong việc phân tích động học. Tác giả áp dụng định nghĩa điểm cân bằng ổn định tiệm cận, mặt mức năng lượng Lyapunov và đặc biệt là Nguyên lý bất biến LaSalle. Các khái niệm cốt lõi bao gồm: đường cong Hugoniot biểu diễn tập hợp các trạng thái có thể kết nối, điểm yên ngựa trong không gian pha, nút ổn định tiệm cận, và miền hấp thụ cực đại biểu thị tập hợp mọi điểm xuất phát có quỹ đạo hội tụ về gốc cân bằng khi biến thời gian quy ước tiến dần ra dương vô cực.

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

Nghiên cứu sử dụng phương pháp giải tích định tính phương trình vi phân kết hợp với công cụ mô phỏng số hiện đại trên máy tính nhằm đối sánh và trực quan hóa kết quả lý thuyết.

Nguồn dữ liệu và mô hình khảo sát là hệ phương trình Euler nén được một chiều trong tọa độ Lagrange với 2 tham số nhiễu cấp cao: hệ số nhớt $\alpha > 0$ và hệ số mao dẫn $\beta > 0$. Về quy trình phân tích, tác giả sử dụng phép đổi biến tự tương tự dạng sóng lưu động $y = x - st$ với vận tốc sóng $s > 0$. Phép biến đổi này chuyển hóa hệ phương trình đạo hàm riêng phi tuyến ban đầu thành hệ 2 phương trình vi phân thường cấp một phi tuyến trên mặt phẳng pha theo biến thể tích riêng $v$ và đạo hàm không gian $w = v'$.

Lý do lựa chọn phương pháp hàm thế Lyapunov và Nguyên lý bất biến LaSalle là bởi phương pháp này cho phép khẳng định chắc chắn sự tồn tại của quỹ đạo dị đồi kết nối 2 điểm cân bằng mà không đòi hỏi phải giải tường minh nghiệm giải tích vốn vô cùng phức tạp hoặc bất khả thi. Để minh họa thực nghiệm, mô phỏng số được tiến hành với 1000 bước tích phân trên bộ công cụ ODE45 của phần mềm MATLAB, lựa chọn trạng thái cân bằng ổn định tại điểm $(10, 0)$, entropy khí $S = 113.3$ và độ lệch kích động ban đầu $0.01$ đơn vị từ điểm yên ngựa.

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

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

Nghiên cứu đã đạt được 3 phát hiện toán học mang tính đột phá trong việc làm rõ cấu trúc vi mô của sóng nén chất lưu:

Thứ nhất, phân loại chính xác bản chất hình học của 2 điểm cân bằng trong hệ động lực sóng lưu động. Tại trạng thái bên trái $(v_-, 0)$, ma trận Jacobi luôn sở hữu 2 giá trị riêng thực trái dấu ($\lambda_1 < 0 < \lambda_2$), xác định đây là một điểm yên ngựa không ổn định. Ngược lại, tại trạng thái bên phải $(v_+, 0)$, ma trận Jacobi có 2 giá trị riêng thực âm hoặc 2 giá trị riêng phức với phần thực âm, khẳng định $(v_+, 0)$ là một điểm cân bằng ổn định tiệm cận bền vững.

Thứ hai, xây dựng thành công hàm Lyapunov xác định dương trên miền xác định $D$ và ước lượng chính xác biên của miền hấp thụ $\Omega$. Nhờ vào cấu trúc hình học của 2 nửa hình elip tiếp xúc, tác giả đã chứng minh được rằng điểm yên ngựa $(v_-, 0)$ nằm chính xác trên đường biên của miền hấp thụ, tạo tiền đề cho sự kết nối quỹ đạo.

Thứ ba, thiết lập điều kiện tham số đủ để tồn tại sóng lưu động toàn cục. Quỹ đạo pha duy nhất rời khỏi điểm yên ngựa đi vào góc phần tư thứ nhất được chứng minh sẽ bị hút 100% vào miền hấp thụ và hội tụ hoàn toàn về điểm cân bằng $(v_+, 0)$ khi thỏa mãn bất đẳng thức tỷ số giữa độ nhớt $\alpha$, mao dẫn $\beta$, vận tốc sóng $s$ và chỉ số nhiệt dung $\gamma$.

Thảo luận kết quả

Cơ chế vật lý sâu xa của kết quả nghiên cứu nằm ở sự cân bằng hoàn hảo giữa 2 hiệu ứng đối nghịch: tiêu tán năng lượng do độ nhớt $\alpha$ gây ra và phân tán sóng phi tuyến do sức căng mao dẫn $\beta$ chi phối. Khi chỉ có độ nhớt, biên dạng sóng thường biến thiên đơn điệu; tuy nhiên, sự xuất hiện của hệ số mao dẫn bậc 3 tạo ra hiện tượng uốn lượn tinh vi trước khi dung hòa vào trạng thái ổn định.

So sánh với các công trình nghiên cứu kinh điển trước đây vốn chỉ khảo sát sóng sốc yếu hoặc bỏ qua hiệu ứng nhiệt động lực học, luận văn đã mở rộng thành công cho sóng 3-sốc mạnh với phương trình trạng thái khí tổng quát có biến thiên entropy.

Dữ liệu nghiên cứu và cấu trúc không gian pha có thể được trình bày trực quan thông qua 2 dạng công cụ:

• Biểu đồ chân dung pha hai chiều: Thể hiện quỹ đạo dị đồi xuất phát từ điểm yên $(v_-, 0)$, đi qua góc phần tư thứ nhất, cắt qua các đường cong mức năng lượng Lyapunov dạng elip và xoáy tụ tiệm cận về điểm nút ổn định $(v_+, 0)$.
• Bảng tổng hợp quan hệ bước nhảy và tham số phổ: Trình bày chi tiết các giá trị biến thiên của thể tích riêng từ $v_-$ đến $v_+$, độ sụt giảm áp suất, mức gia tăng entropy và dấu của các giá trị riêng ma trận Jacobi tương ứng với từng trạng thái cân bằng.

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

Dựa trên các kết quả giải tích và mô phỏng đã đạt được, 4 khuyến nghị hành động cụ thể được đề xuất nhằm chuyển giao và mở rộng nghiên cứu:

Thứ nhất, nâng cấp mô hình giải tích sang không gian 2 chiều và 3 chiều có xét đến tính dị hướng của dòng chảy. Nhóm nghiên cứu phương trình vi phân tại các viện toán học ứng dụng cần chủ trì triển khai trong lộ trình 12 tháng, đặt mục tiêu bao phủ thêm 90% các hình thái xoáy phức tạp trong cơ học chất khí.

Thứ hai, tích hợp cấu trúc sóng lưu động nhớt - mao dẫn vào các thuật toán nắm bắt sóng sốc trong phần mềm mô phỏng dòng chảy động lực học (CFD). Đội ngũ kỹ sư tính toán cơ khí và hàng không cần thực hiện trong vòng 6 tháng nhằm tối ưu hóa thuật toán, hướng tới mục tiêu giảm 25% hiện tượng dao động giả tạo (spurious oscillations) tại mặt phân cách sóng.

Thứ three, tiến hành thực nghiệm đo đạc các hệ số nhớt và mao dẫn động lực trên các hệ chất lưu đa pha thực tế (như hỗn hợp dầu khí, dung dịch polymer nén). Các phòng thí nghiệm cơ học vật liệu cần phối hợp thực hiện trong thời gian 18 tháng để chuẩn hóa bộ dữ liệu tham số với độ tin cậy vượt trên 95%.

Thứ tư, áp dụng tiêu chuẩn ổn định Lyapunov và miền hấp thụ vào thiết kế hệ thống van giảm áp và an toàn đường ống dẫn khí nén áp lực cao. Các kỹ sư vận hành công trình công nghiệp cần đưa quy trình kiểm tra vào kế hoạch bảo dưỡng định kỳ 3 tháng một lần, giảm thiểu nguy cơ nứt vỡ đường ống do búa áp lực xuống dưới mức 0.5%.

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

Công trình luận văn mang lại giá trị học thuật và ứng dụng thực tiễn to lớn cho 4 nhóm đối tượng trọng điểm:

Nhóm thứ nhất là các nhà toán học, giảng viên và nghiên cứu sinh chuyên ngành Toán giải tích, Toán ứng dụng và Phương trình đạo hàm riêng. Luận văn cung cấp phương pháp luận mẫu mực về kỹ thuật xây dựng hàm Lyapunov phi tuyến và ứng dụng nguyên lý LaSalle trong hệ động lực vô hạn chiều.

Nhóm thứ hai là các kỹ sư phát triển phần mềm mô phỏng động lực học chất lưu tính toán (CFD). Kết quả nghiên cứu giúp cải tiến mô hình toán học nền tảng, hỗ trợ thiết kế khí động học cho tuabin, động cơ phản lực và các thiết bị bay siêu thanh với độ chính xác cao.

Nhóm thứ ba là các chuyên gia kỹ thuật dầu khí, năng lượng và công trình thủy lực đường ống. Tài liệu hỗ trợ trực tiếp cho use case dự báo và kiểm soát sóng va đập, hiện tượng chuyển pha lỏng - khí đột ngột trong các đường ống dẫn dầu khí đa pha đường dài.

Nhóm thứ tư là học viên cao học và sinh viên các ngành Khoa học tính toán, Vật lý lý thuyết. Luận văn là tài liệu tham khảo chất lượng cao về cách thức kết hợp chặt chẽ giữa suy diễn giải tích toán học thuần túy với kỹ năng lập trình số hóa trên môi trường MATLAB.

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

Phương trình Euler nén được có nhớt và mao dẫn khác gì so với phương trình Euler truyền thống? Phương trình Euler truyền thống bỏ qua hoàn toàn các hiệu ứng ma sát trong và sức căng bề mặt, dẫn đến việc hình thành các nghiệm gián đoạn dốc đứng không khả vi. Việc bổ sung số hạng nhớt cấp hai $\alpha u_{xx}$ và mao dẫn cấp ba $\beta v_{xxx}$ giúp làm mịn mặt sóng gián đoạn, mô tả chính xác cấu trúc chuyển tiếp vi mô của dòng khí nén trong thực tế.

Ý nghĩa cốt lõi của việc chứng minh sự tồn tại của sóng lưu động là gì? Chứng minh sự tồn tại của sóng lưu động tương đương với việc xác nhận về mặt toán học rằng sóng sốc vật lý thực sự sở hữu một cấu trúc nội tại ổn định và có thể lan truyền đều đặn với vận tốc không đổi $s$, không bị tan vỡ hay suy biến sau thời gian dài vận động qua môi trường liên tục.

Nguyên lý bất biến LaSalle đóng vai trò then chốt như thế nào trong bài toán này? Nguyên lý LaSalle cho phép xác định tập giới hạn dương của mọi quỹ đạo xuất phát từ một tập compact bất biến dương. Trong luận văn, nguyên lý này giúp chứng minh rằng mọi quỹ đạo rời khỏi điểm yên ngựa đi vào tập mức năng lượng Lyapunov bắt buộc phải tiệm cận về điểm cân bằng duy nhất $(v_+, 0)$ chứ không rơi vào các chu trình giới hạn tuần hoàn.

Điều kiện sóng sốc Lax có ý nghĩa vật lý ra sao đối với dòng khí nén? Điều kiện sốc Lax đóng vai trò như một bộ lọc entropy, đảm bảo rằng sóng sốc là một quá trình nén thuận chiều nhiệt động lực học với áp suất và mật độ tăng dần từ trước ra sau mặt sóng, loại trừ hoàn toàn các sóng giãn phi vật lý không thể xuất hiện trong tự nhiên.

Mô phỏng số bằng thuật toán ODE45 trong MATLAB được thiết lập với các thông số nào? Nghiên cứu sử dụng thuật toán Runge-Kutta thích nghi bậc 4-5 của hàm ODE45 để giải hệ 2 phương trình vi phân phi tuyến. Mô phỏng kích hoạt từ lân cận điểm yên với độ dịch chuyển $0.01$ đơn vị, tích phân qua 1000 bước thời gian, mô tả quá trình quỹ đạo hội tụ trơn tru về nút ổn định $(10, 0)$ với entropy khí đạt $113.3$.

Kết luận

Năm điểm cốt lõi đúc kết từ công trình nghiên cứu bao gồm:

• Thiết lập thành công mô hình hệ động lực hai chiều mô tả chính xác sóng lưu động cho phương trình Euler nén được có nhớt và mao dẫn.
• Xác định bản chất điểm yên ngựa của trạng thái đầu $(v_-, 0)$ và tính ổn định tiệm cận của trạng thái cuối $(v_+, 0)$.
• Xây dựng hàm Lyapunov xác định dương và ước lượng tối ưu biên miền hấp thụ của nút ổn định.
• Chứng minh chặt chẽ sự tồn tại toàn cục của quỹ đạo pha nối 2 trạng thái cân bằng dưới điều kiện bất đẳng thức tham số nhớt - mao dẫn.
• Mô phỏng kiểm chứng số học thành công 100% sự hội tụ của sóng lưu động trên phần mềm MATLAB.

Đóng góp chính của luận văn là hoàn thiện cơ sở giải tích cho lý thuyết sóng nén chất lưu vi mô. Lộ trình phát triển tiếp theo trong 6 đến 12 tháng tới sẽ tập trung vào việc mở rộng mô hình cho chất lưu phi Newton và tích hợp vào các hệ thống tính toán kỹ thuật số. Bạn đọc và các nhà nghiên cứu quan tâm hãy tiếp tục khai thác khung lý thuyết này để ứng dụng vào các bài toán mô phỏng thủy khí động lực học chuyên sâu.

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

Chương 1: kiến thúc tổng quan Luan van cao hoc Chitng minh Lay € > 0, chon r € (0, €] sao cho: B, = {x € R"|||z|| <r} CD. Ta sé chttng t6 Qg nim trong phan trong của B,. That vay, gia sit ngudc lai Qg khong nằm trong phần trong của Ö,. Khi đó, tồn tại p€®©¿ñn®; với S%„ là biên của B,.

Vay Qg nam trong phần trong của Ö,. Tập Ô¿ có tính chất mọi quỹ đạo xuất phát trong ©z phải nằm trong ©;,Vý > 0. Điều này được suy ra từ tính chất bán xác định âm của V: V(a(t)) <0 > V(a(t)) < V(a2(0)) < 8, Vt>0. Vì hàm V liên tục trên ,„ CD nên tập ©¿ đóng và bị chặn trên R” do do la tap compact.

Chú ý hệ phương trình vi phân (1.12) luôn có nghiệm duy nhất xác định với mọi £ dương khi z(0) € Đa. Khi đó Bs Cc Qz C Ö, và z(0) € B¿ > z(0) € Ô¿ = z(£) € Ng => z(0) € „, suy ra ||z(0)||< ô= zŒ)<r<c, t>0. Điều này chứng tỏ điểm cân bằng z = 0 là ổn định. Tiếp theo, ta giả thiết thêm W xác định âm, ta sẽ chứng minh điểm cần bằng z = 0 là ổn định tiệm cận, tức là ta cần chứng tỏ lim #(#) = 0 hay t++00 Va >0,37 >0: Vý > 7 == |[x(t)| <a.

trang 18 Thạch Phước Chương 1: kiến thúc tổng quan Luan van cao hoc Lặp lại chứng minh trên, a > 0 ta chọn b > 0 sao cho Q C By. Ta chỉ cần chỉ ra lim V(z(£)) =0 too Vi V gidm va bi chin dudi béi 0 nên lim V(2(t)) =c>0 too Ta dùng phản chứng, giả sử e > 0. Do WV(z) liên tục và V(0) = 0 nên tồn tại đ > 0 sao cho Ba C Qe. Mặt khác, ham V giam nén V(2(t)) > c, Vt > 0.

Diéu này chứng tỏ quỹ đạo z(#) nam ngoai hinh cau By véi moi t > 0. Do tinh liên tục trên tập compact của V nên ta có thể đặt —+= max V(x) =V(ao) <0 đ<l|lzll<z Từ đây suy ra V(a(t)) = V(x(0)) + [ear < V(ax(0)) — yt, Vt > 0. Vì + > 0 nên ta cho đủ lớn sẽ làm cho về phải âm. Điều này mâu thuần với V(()) > c, Vt > 0.

Vậy phản chứng sai và từ đây cho ta kết luận điểm cân bằng z = 0 là ổn định tiệm cận. Ghi chi: Mat V(x) = c gọi là mặt mức hay còn gọi là mặt Lyapunov. Điều kiện bán xác định âm của V dẫn đến kết quả là một quỹ đạo đã đi vào phía trong của mặt mức thì sẽ không bao giờ ra ngoài được nữa.2 Miền hấp thụ của điểm cân bằng ổn định tiệm cận Bây giờ ta nghiên cứu tính chất của các quỹ đạo một khi gốc ổn định tiệm cận, tức là nếu một quỹ đạo đã đi vào gốc đủ gần thì sẽ hội tụ về gốc khi ¿ — 0. Vấn đề là khoảng cách đó là bao xa.

Nếu khoảng cách đó là toàn mặt phẳng (z,£) thì ta nói gốc ổn định toàn cục. Còn trong trường hợp không ổn định toàn cục thì ta sẽ ước lượng một lân cận của gốc (lớn nhất có thể) mà một khi quỹ đạo đã đi vào đó thì sẽ hội tụ về gốc, miền này được gọi là mién hấp thụ của điểm cân bằng z = 0.2 Giả sử z(t,#o) kú hiệu là nghiệm +(t) của hệ phương trình ơì phân (1. Ta gọi miền hấp thụ của gốc là tập hợp tất cả các điểm #ụ sao cho z(t,œo) tồn tại toàn cục uới mọi † > 0 va jim x(t, xo) =0. —> 00, trang 19 Thạch Phước Chương 1: kiến thức tổng quan Luan van cao hoc Từ phép chứng minh của định lý Lyapunov, ta thấy rằng các tập mức ứng với hàm Lyapunov V xác định âm nằm trong miền hấp thụ.

Do vậy, đây là một cơ sở để chúng ta ước lượng một miền hấp thụ tối đa có thể. Ta có kết quả sau đây về miền hấp thụ. Khi đó nếu 9,={ze"| V(z) < c} nằm trong D thà Q„ cũng nằm trong miền hấp thụ. Túc là quỹ đạo #(t,zo) nằm hoàn toàn trong ©„ uà hội tụ uề gốc khi † > oo.

Tiếp theo ta sẽ nghiên cứu khi nào thì điểm cân bằng gốc là ổn định toàn cục, tức là miền hấp thụ được mở rộng qua toàn #†". Hay nói cách khác, mọi quỹ đạo #(£) luôn hội tụ về gốc khi # —› oo bất chấp điểm xuất phát z(0) = zo xa gốc bao nhiêu. Ta có định lý sau: Định lý 1.3 (On định toàn cục) Cho øọ là điểm cân bằng ổn định tiệm cận của hệ phương trình ui phân (1. Giả sử tồn tại hàm Lyapunov V : R" > R thoa mãn V(0) =0,V(z) > 0,Vz # {0}, lim V(z) = ©, llzll->s Ÿ <0,Vz # {0}.

Khi đó z = 0 ổn định toàn cục. Chứng minh Lay bat ky p € R", dat c = V(p). Diéu kién gia thiét 7 lim ||V (x)|| = co tuong duong véi #||>% Ve>0,3r >0: ||xz|| >r = V(#) > c. Điều này chứng tổ 9„ CC Đ,.

Lặp lại phép chứng minh trong định lý Lyapunov, ta dé đàng chứng minh được # = 0 ổn định toàn cục. trang 20 Thạch Phước Chương 1: kiến thúc tổng quan Luan van cao hoc 1.3 Nguyên lý bất biến LaSalle Cho z() là quỹ đạo của hệ phương trình vi phân (1. Một điểm p gọi là điểm giới hạn đương của quỹ đạo x(t) nếu tồn tại dãy (fa)a : t„ — oo thỏa mãn #(t„) — ?. Tập tất cả các điểm giới hạn dương gọi là tập giới hạn dương của x(t).

Một tập Ä gọi là bất biến của quỹ đạo z(£) nếu z(0) € M =— +(£) € M,Vt c R. Mot tap M gọi là bất biến dương của quỹ đạo z(f) nêu z(0) € M => z(0) € M,V+ > 0. Zz 2 ^ Ta co bo dé sau Bổ đề 1.4 Cho x(t) la mot qui dao ctia hé phuong trink vi phan (1.12) nằm trong D va bi chặn uới mọi t > 0. Khi đó, tập giới hạn dương LẺ của quỹ dao x(t) la một tập khác rỗng, compact uà bất biến.

Hơn nữa, +(E) tiếp cận LT khi † — oo. Chứng minh Do z(£) là dãy bị chặn trên ” nên theo dinh ly Bolzano - Weierstrass, {x(t),t > 0} có ít nhất một điểm tụ, điều này chứng tỏ tập giới hạn dương 7T khác rỗng. Vì {z(),t > 0} bị chặn nên T cũng là tập bị chặn. Để chứng minh L* compact, ta cần chứng tỏ nó đóng.

Lấy số dương e nhỏ tùy ý Xét {y;} € L† : y; > y, ta sẽ chứng minh € LT. Thật vậy, với mỗi ¡ € Ñ luôn tồn tại dãy (;) : f¿; —> oo sao cho mm r(ti;) = Yi. V6i méi 7 € N, ta chon 7; > t;, sao cho € le) — mÌl < Š Dãy (7;) được chọn ở trên có tính chất 7¡ > oo. Dựa vào các xây dưng dãy (7;) và 1; — # ta có thể chọn 2 số Ny, No > 0 thoa man Vi > Ni, : ||#() — || <e/2, va t> No: |ly; — ylle/2.

trang 21 Thạch Phước Chương 1: kiến thúc tổng quan Luan van cao hoc Từ đây suy ra ||z(¡) — || <«,_ Vi> N=maz{M,N:} Điều này chứng tỏ lim «(7;) = ÿ. Vậy LÝ là tập đóng va do dé 1a tap compact. 1-00 Ta con phai chttng minh Lt 1a tap bat bién. Cho y € Lt va ¢(t; y) lA một quỹ đạo của hệ phương trình vi phân (1.12), ta cần chứng tỏ O(t;y) € Lt, Vt € R.

Vìg€ D7 nên tồn tai day (t;) : t; + 00 sao cho a(t;) = o(ti; 20) = y. Do tinh duy nhất của giới hạn nên ta có P(t + ti; Zo) = O(t; O(ti; zo)) = ð(¡#()) trong đó £ + f; > 0 với ? đủ lớn. Bởi tính liên tục nên ta có lim 9(t + ti; 0) = lim ó(f; #(;)) — óÁ 9): 1-00 Điều này chứng tỏ o(t;y) € L+. Cuối cùng, để chứng minh z(#) > Lt, ta dimg phản chứng.

Giả sử ngược lại thì tồn tại số e > 0 và dãy f¿ — oo sao cho đ(#(¿), b*) > e. Vậy a(t) => E* hay LT là tập bất biến, ta có điều phải chứng minh. Tiếp theo, ta sẽ chứng minh nguyên lý bất biến LaSalle, một cơ sở để chứng minh sự tồn tại sóng lưu động cho chương sau.5 (Nguyén ly bat bién LaSalle) Cho 0 C D la tép compact va bat biến dương của hệ phương trình ơi phân (1. Ham V: D> R kha vi lién tuc trén D sao cho V(x) < 0,Va €Q.

Xét tập E={xzeO:V(z) =0} Goi M là tập bất biến lớn nhất trong F. Khi đó mọi quỹ đạo của (1.12) xuất phát trong © đều tiếp cận M khit 3 co. Chứng minh trang 22 Thạch Phước Chương 1: kiến thúc tổng quan Luan van cao hoc Ta xét x(t) 1A mot quỹ đạo tùy ý của (1.12) xuất phát trong 2. Vì V(z) liên tục trên tập compact 2 nén V(a(t)) bi chan dưới.

t-00 Xét tùy ý điểm p € LT, theo định nghĩa Ƒ* nên tồn tại dãy (f„) : #(f„) —> p. Do tính liên tục nên ta có V(p) = lim V{(zŒ,)) = a. noo Điều này chứng tổ hàm W(+) luôn nhận giá trị hằng trên L*. Cũng bởi tính bất biến của U mà suy ra được z(£) tiếp cận À.

Nhận xét: Dịnh lý trên không đòi hỏi hàm W(+) xác định dương. Vì z(#) tiếp cận tập 1M nên khi cần chứng mình z(£) —> 0, —> oo thì ta chỉ cần chứng tỏ tập bất biến lớn nhất trong # chỉ chứa đúng một điểm cân bằng z = 0. Diều này có thể đạt được bằng cách chứng tỏ không có một nghiệm nào nằm hoàn toàn trong # ngoại trừ nghiệm tầm thường z(£) = 0. Chú ý thêm rằng, hàm Lyapunov V xác định âm trên 2 nên V(z@))<V((0)), t>0.

Và từ đó ta có tập mức O.={z| V(r) <c}, c>0 bị chặn trong D nên là tập compact và bất biến dương. Định lý trên được thêm giả thiết V xác định dương thì ta có thêm kết quả bổ ích sau Hệ quả 1.6 (Ước lượng miền hấp thụ) Cho z = 0 là điểm cân bằng của hệ phương trình ưi phân (1. Giả sử tồn tại hàm Lụapunou V trên D 3 0. Đặt S={zeD| V(z) =0} Nếu hệ (1.12) không có một nghiệm nào nằm hoàn toờn trong S ngoại trừ nghiệm tầm thường z(Ð) = 0 thà z = 0 là điểm ổn định tiệm cận.

Hơn nữa mọi tập mức ©„ bị chăn chứa trong D đều là tập con của miền hấp thụ của z = 0. trang 23 Thạch Phước Chương 1: kiến thức tổng quan Luan van cao hoc Hé qua 1.7 Cho x =0 là điểm cần bằng của hệ phương trình vi phân (1.