Một số phương pháp symplectic cho hệ Hamilton

Luận văn thạc sĩ: Một số phương pháp symplectic cho hệ Hamilton - ĐH Khoa học Thái Nguyên 2021

• Tác giả: Trần Văn Cẩn
• Người hướng dẫn: TS. Nguyễn Thanh Sơn
• Trường học: Trường Đại học Khoa học - Đại học Thái Nguyên
• Chuyên ngành: Toán ứng dụng
• Đề tài: Một số phương pháp symplectic cho hệ Hamilton
• Loại tài liệu: Luận văn thạc sĩ
• Năm xuất bản: 2021
• Địa điểm: Thái Nguyên

Tóm tắt

I. Tổng quan về phương pháp symplectic cho hệ Hamilton

Hệ Hamilton là mô hình toán học quan trọng trong cơ học và vật lý. Hệ mô tả nhiều hiện tượng tự nhiên, từ dao động con lắc đến truyền sóng phi tuyến. Điểm cốt lõi của hệ là hàm Hamilton, biểu diễn năng lượng toàn phần của hệ. Năng lượng này được bảo toàn theo thời gian. Đây là tính chất đặc trưng quan trọng nhất.

Trong thực tế, nghiệm tường minh của hệ Hamilton hiếm khi tồn tại. Các phương pháp số thông thường, chẳng hạn Runge-Kutta cổ điển, không đảm bảo bảo toàn cấu trúc hình học của nghiệm. Chúng có thể dẫn đến sai lệch tích lũy theo thời gian, làm mất đi ý nghĩa vật lý của kết quả mô phỏng.

Phương pháp symplectic ra đời để giải quyết vấn đề này. Các phương pháp này được thiết kế chuyên biệt, có khả năng bảo toàn tính chất symplectic của dòng Hamilton. Luận văn thạc sĩ của tác giả Trần Văn Cẩn, thực hiện tại Trường Đại học Khoa học Thái Nguyên năm 2021, trình bày hệ thống ba phương pháp symplectic chính: Euler symplectic, Điểm giữa và Störmer-Verlet. Các phương pháp được kiểm chứng trên ba mô hình cụ thể: phương trình sóng, sine-Gordon và Schrödinger.

1.1. Hàm Hamilton và tính bảo toàn năng lượng

Hàm Hamilton H(q, p) đặc trưng cho năng lượng toàn phần của hệ, trong đó q là tọa độ suy rộng và p là động lượng suy rộng. Tính bảo toàn năng lượng có nghĩa là H dọc theo nghiệm luôn là hằng số. Đây là ràng buộc vật lý quan trọng. Mọi phương pháp số giải hệ Hamilton đều cần tôn trọng ràng buộc này, nếu không kết quả mô phỏng dài hạn sẽ mất ý nghĩa. Phương pháp symplectic bảo toàn hoặc bảo toàn xấp xỉ tính chất này một cách có kiểm soát.

1.2. Tính symplectic của dòng Hamilton là gì

Tính symplectic của dòng Hamilton mô tả khả năng bảo toàn tổng diện tích có hướng của các hình chiếu lên các mặt phẳng tọa độ trong không gian pha. Về mặt toán học, dòng Hamilton φt thỏa mãn điều kiện (Dφt)^T J (Dφt) = J, với J là ma trận Poisson. Đây là điều kiện symplectic. Tính chất hình học này phân biệt hệ Hamilton với các hệ vi phân thông thường và là cơ sở lý thuyết để xây dựng các phương pháp tích phân số chuyên dụng.

II. Phân tích hạn chế của phương pháp số thông thường

Các phương pháp số cổ điển, như Euler tường minh hay Runge-Kutta bậc cao, được xây dựng để đảm bảo độ chính xác cục bộ và tính ổn định. Tuy nhiên, chúng không được thiết kế để bảo toàn cấu trúc hình học của hệ Hamilton. Khi áp dụng các phương pháp này trong thời gian dài, năng lượng tính toán thường bị trôi dạt, tăng hoặc giảm liên tục.

Sự trôi dạt năng lượng này không phải sai số ngẫu nhiên. Đây là hệ quả tất yếu khi lược đồ số không bảo toàn tính symplectic. Trong các bài toán vật lý dài hạn, chẳng hạn mô phỏng quỹ đạo thiên văn hay dao động phân tử, sai lệch này tích lũy và làm nghiệm số hoàn toàn sai lệch so với nghiệm thực.

Các tiêu chí đánh giá lược đồ số gồm: sự hội tụ, tính chính xác và tính ổn định tuyệt đối. Với hệ Hamilton, cần thêm tiêu chí thứ tư: khả năng bảo toàn bất biến. Phương pháp Runge-Kutta thông thường đáp ứng ba tiêu chí đầu nhưng thất bại ở tiêu chí thứ tư. Đây là lý do cần phát triển lớp phương pháp symplectic riêng biệt.

2.1. Sự hội tụ và ổn định của lược đồ sai phân

Một lược đồ sai phân được gọi là hội tụ nếu sai số toàn cục tiến về không khi bước thời gian h tiến về không. Tính ổn định tuyệt đối yêu cầu nghiệm số không phát sinh dao động tăng trưởng không kiểm soát. Các phương pháp Runge-Kutta đảm bảo hai tính chất này ở mức độ cao. Tuy nhiên, đối với hệ Hamilton, hội tụ và ổn định là chưa đủ. Cần thêm điều kiện bảo toàn cấu trúc symplectic để nghiệm số phản ánh đúng động lực học dài hạn của hệ.

2.2. Sai lệch năng lượng trong mô phỏng dài hạn

Khi dùng phương pháp Euler tường minh để giải hệ Hamilton, năng lượng tính toán thường tăng đơn điệu theo thời gian mô phỏng. Ngược lại, Euler ẩn khiến năng lượng giảm. Cả hai đều phản ánh sai thực tế vật lý. Hiện tượng này gọi là trôi dạt năng lượng. Trong bài toán con lắc đơn hay phương trình sóng phi tuyến, trôi dạt năng lượng dẫn đến biên độ dao động sai lệch sau vài nghìn bước thời gian. Phương pháp symplectic khắc phục vấn đề này bằng cách giữ năng lượng dao động gần giá trị ban đầu trong suốt quá trình tính toán.

III. Các phương pháp symplectic chính cho hệ Hamilton

Luận văn trình bày ba phương pháp symplectic cụ thể. Mỗi phương pháp có cơ sở lý thuyết riêng và đặc điểm tính toán khác nhau. Tất cả đều bảo toàn tính chất symplectic của dòng Hamilton theo nghĩa toán học chặt chẽ.

Phương pháp Euler symplectic là đơn giản nhất. Nó cập nhật tọa độ và động lượng theo thứ tự xen kẽ, không đồng thời. Phương pháp này bậc một, đơn giản nhưng hiệu quả vì chi phí tính toán thấp.

Phương pháp Điểm giữa (Midpoint Rule) sử dụng giá trị trung bình tại điểm giữa của bước thời gian. Phương pháp này bậc hai và bảo toàn năng lượng chính xác hơn. Tuy nhiên, nó đòi hỏi giải hệ phi tuyến ẩn ở mỗi bước, tốn thời gian tính toán hơn.

Phương pháp Störmer-Verlet được đánh giá là tốt nhất trong ba phương pháp. Nó cũng bậc hai, tường minh với phần tọa độ, bảo toàn năng lượng rất tốt và có thời gian tính toán ngắn hơn phương pháp Điểm giữa. Kết quả số trên các mô hình sóng, sine-Gordon và Schrödinger xác nhận kết luận này.

3.1. Phương pháp Euler symplectic và quy tắc điểm giữa

Phương pháp Euler symplectic tách hệ Hamilton thành hai bước: cập nhật động lượng p trước, sau đó cập nhật tọa độ q sử dụng giá trị p mới. Sự tách biệt này tạo ra tính symplectic của lược đồ. Phương pháp có bậc chính xác bằng một. Quy tắc điểm giữa đánh giá vế phải của hệ tại điểm giữa của bước thời gian. Lược đồ này bậc hai và bảo toàn năng lượng tốt hơn. Tuy nhiên, cần giải hệ phi tuyến ẩn ở mỗi bước bằng lặp Newton, làm tăng đáng kể thời gian tính toán, đặc biệt với bài toán lớn.

3.2. Phương pháp Störmer Verlet và ưu điểm thực tiễn

Phương pháp Störmer-Verlet cập nhật theo ba bước tuần tự: bán bước cho p, bước đầy đủ cho q dùng p mới, rồi bán bước tiếp theo cho p. Lược đồ này tường minh hoàn toàn khi hàm Hamilton có dạng tách biệt H = T(p) + V(q). Bậc chính xác là hai. Kết quả thực nghiệm trên ba mô hình trong luận văn cho thấy Störmer-Verlet đạt cân bằng tốt nhất giữa độ chính xác bảo toàn năng lượng và thời gian tính toán. Với bài toán sóng, thời gian chạy chỉ 0.6 giây, thấp hơn nhiều so với phương pháp Điểm giữa mất 1248 giây.

IV. Kết quả ứng dụng và kết luận nghiên cứu

Luận văn kiểm chứng ba phương pháp symplectic trên ba mô hình phương trình vi phân đạo hàm riêng quan trọng. Mỗi mô hình đại diện cho một lớp bài toán vật lý khác nhau.

Phương trình sóng mô tả dao động cơ học. Phương trình sine-Gordon xuất hiện trong vật lý trường lượng tử và liên kết siêu dẫn. Phương trình Schrödinger là nền tảng của cơ học lượng tử. Cả ba đều có thể viết dưới dạng hệ Hamilton sau khi rời rạc hóa không gian.

Kết quả số khẳng định: phương pháp Điểm giữa và Störmer-Verlet bảo toàn giá trị hàm Hamilton tốt hơn rõ rệt so với Euler symplectic. Độ lệch tương đối của năng lượng với hai phương pháp này rất nhỏ và ổn định theo thời gian. Trong khi đó, Euler symplectic có độ lệch lớn hơn, mặc dù vẫn tốt hơn các phương pháp cổ điển không symplectic.

Xét tổng thể về độ chính xác lẫn hiệu quả tính toán, phương pháp Störmer-Verlet là lựa chọn tốt nhất. Đây là kết luận thực tiễn có giá trị cho các bài toán mô phỏng vật lý dài hạn.

4.1. So sánh hiệu quả trên các mô hình vật lý

Trên mô hình phương trình sóng với n=500 nút không gian và 5000 bước thời gian, phương pháp Störmer-Verlet chạy trong 0.6 giây trong khi Điểm giữa cần tới 1248 giây. Cả hai đều bảo toàn hàm Hamilton tốt hơn Euler symplectic. Với mô hình sine-Gordon và Schrödinger, kết quả tương tự. Phương pháp Điểm giữa đôi khi chính xác hơn một chút nhưng chi phí tính toán quá cao khiến nó kém khả thi trong ứng dụng thực tế với bài toán quy mô lớn.

4.2. Định hướng phát triển và ứng dụng thực tiễn

Các phương pháp symplectic trong luận văn mở ra hướng nghiên cứu quan trọng. Có thể phát triển các phương pháp bậc cao hơn, chẳng hạn phương pháp tổ hợp symplectic bậc bốn. Ngoài ra, có thể mở rộng sang các hệ Hamilton suy biến hoặc hệ có ràng buộc. Ứng dụng thực tiễn rất đa dạng: mô phỏng động lực học phân tử, tính toán quỹ đạo vệ tinh, mô phỏng plasma và thiết bị lượng tử. Trong tất cả các lĩnh vực này, bảo toàn năng lượng dài hạn là yêu cầu bắt buộc, làm cho phương pháp symplectic trở thành công cụ không thể thiếu.

Nội dung chính

Tổng quan nghiên cứu

Hệ Hamilton là mô hình toán học tổng quát giữ vai trò hạt nhân trong cơ học giải tích, vật lý plasma, thiên văn học và lý thuyết trường lượng tử. Trong không gian pha 2n chiều xác định bởi tọa độ suy rộng và xung lượng, hơn 90% các hệ động lực Hamilton phi tuyến trong thực tế không thể giải bằng phương pháp giải tích để thu được nghiệm tường minh. Khi đó, việc sử dụng các phương pháp số rời rạc hóa trên máy tính trở thành lựa chọn tất yếu. Tuy nhiên, các phương pháp số giải phương trình vi phân cổ điển thường bỏ qua bản chất hình học của hệ, dẫn đến sai số trôi dạt năng lượng nghiêm trọng. Cụ thể, khi thử nghiệm phương pháp Euler hiện trên mô hình con lắc đơn qua 1000 bước tính với khoảng thời gian T = 40, mức năng lượng toàn phần liên tục tăng vọt và độ lệch tương đối khỏi trạng thái ban đầu mất kiểm soát hoàn toàn.

Nhằm giải quyết triệt để thách thức trên, luận văn thạc sĩ Toán học chuyên ngành Toán ứng dụng (mã số 8 46 01 12) của tác giả Trần Văn Cẩn, thực hiện tại Trường Đại học Khoa học - Đại học Thái Nguyên năm 2021 dưới sự hướng dẫn của TS. Nguyễn Thanh Sơn, đã tập trung nghiên cứu đề tài "Một số phương pháp symplectic cho hệ Hamilton". Mục tiêu cụ thể của công trình là xây dựng cơ sở lý thuyết, phân tích tính ổn định và kiểm nghiệm hiệu năng thực tế của ba lược đồ tích phân bảo toàn hình học gồm Euler symplectic, Điểm giữa và Störmer-Verlet. Phạm vi nghiên cứu khảo sát từ bài toán con lắc đơn trong không gian 2 chiều đến 3 hệ phương trình đạo hàm riêng phức tạp gồm phương trình sóng, sine-Gordon và Schrödinger phi tuyến. Kết quả nghiên cứu mang ý nghĩa học thuật và ứng dụng to lớn khi chứng minh các phương pháp symplectic có khả năng kiểm soát độ lệch năng lượng trong ngưỡng bậc 2, đồng thời rút ngắn thời gian tính toán thực tế tới hơn 99% so với các lược đồ số truyền thống.

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

Khung lý thuyết áp dụng

Nghiên cứu được xây dựng vững chắc trên hai trụ cột lý thuyết lớn: lý thuyết hệ động lực Hamilton trên đa tạp đối xứng và lý thuyết tích phân số bảo toàn hình học (Geometric Numerical Integration). Hệ Hamilton chính tắc được xác định thông qua cặp hàm vectơ tọa độ q(t) và xung lượng p(t) trong không gian Euclid n chiều, liên kết trực tiếp với hàm năng lượng Hamilton H(q, p) theo hệ phương trình vi phân đạo hàm riêng: đạo hàm cấp một theo thời gian của q bằng đạo hàm riêng của H theo p, và đạo hàm cấp một theo thời gian của p bằng trừ đạo hàm riêng của H theo q. Khi đặt vectơ trạng thái y bao gồm q và p, hệ phương trình được biểu diễn cô đọng qua tích của ma trận Poisson chính tắc J cỡ 2n nhân với toán tử nabla của hàm Hamilton.

Khung lý thuyết của luận văn tập trung vào 4 khái niệm nền tảng:

• Thứ nhất, dòng Hamilton, tức ánh xạ biến đổi tọa độ trạng thái từ thời điểm ban đầu đến thời điểm t.
• Thứ hai, tính bảo toàn năng lượng toàn phần dọc theo đường cong nghiệm, tức đạo hàm toàn phần của H theo thời gian luôn triệt tiêu bằng 0.
• Thứ ba, dạng symplectic chính tắc omega xác định bởi tích vô hướng của vectơ trạng thái với ma trận Poisson, thể hiện tính bảo toàn tổng diện tích có hướng của các hình chiếu lên các mặt phẳng pha tọa độ.
• Thứ tư, khái niệm ánh xạ bước bảo toàn symplectic, đảm bảo Jacobian của phép chuyển dịch pha luôn duy trì cấu trúc ma trận trực giao symplectic.

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

Luận văn kết hợp phương pháp nghiên cứu giải tích toán học với mô phỏng số thực nghiệm trên máy tính. Để chuyển đổi các phương trình đạo hàm riêng vô hạn chiều thành hệ Hamilton hữu hạn chiều, tác giả áp dụng phương pháp sai phân hữu hạn trung tâm đối xứng theo biến không gian. Kích thước mẫu thực nghiệm số được thiết lập với cỡ mạng lưới từ 256 nút đối với phương trình Schrödinger phi tuyến đến 500 nút lưới đều đối với phương trình sóng và phương trình sine-Gordon. Tiến trình mô phỏng thời gian được thực hiện qua các chuỗi bước lặp từ 1000 đến 5000 bước, với cỡ bước thời gian h cố định trong khoảng từ 0.01 đến 0.02.

Phương pháp phân tích được lựa chọn là phân tích sai số làm gọn địa phương bậc p và tính 0-ổn định theo điều kiện liên tục Lipschitz. Lý do lựa chọn cách tiếp cận này là nhằm bảo đảm tính hội tụ toàn cục của lược đồ sai phân trên toàn bộ miền thời gian khảo sát, đồng thời cho phép liên kết chặt chẽ sai số địa phương với phân tích sai số ngược (backward error analysis). Việc chọn mẫu thực nghiệm dựa trên 4 mô hình tiêu biểu đại diện cho các lớp bài toán cơ bản: hệ dao động phi tuyến 1 bậc tự do (con lắc đơn), phương trình hyperbolic tuyến tính (sóng), phương trình sóng phi tuyến mạnh (sine-Gordon) và phương trình tán xạ lượng tử (Schrödinger). Toàn bộ quá trình nghiên cứu lý thuyết và tính toán thực nghiệm được triển khai và hoàn tất trong năm 2021.

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

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

Quá trình kiểm nghiệm lý thuyết và tính toán mô phỏng trên các mô hình đã mang lại 4 phát hiện quan trọng có giá trị khoa học cao:

Thứ nhất, phương pháp Euler symplectic phân hoạch biến đạt độ chính xác bậc 1 với sai số năng lượng O(h). Lược đồ này bảo tồn chính xác tính symplectic của dòng nhưng mức dao động biên độ năng lượng còn tương đối rõ nét ở các bước thời gian lớn.

Thứ hai, phương pháp Điểm giữa và phương pháp Störmer-Verlet đều đạt bậc chính xác số 2 với sai số năng lượng xấp xỉ O(h mũ 2). Năng lượng toàn phần của hệ thống không bị trôi dạt vô hạn mà được giữ ổn định trong một dải biên độ hẹp được chặn chặt chẽ theo thời gian.

Thứ ba, sự vượt trội áp đảo về hiệu năng thời gian tính toán của phương pháp Störmer-Verlet so với phương pháp Điểm giữa. Trong bài toán phương trình sóng với kích thước mạng 500 nút và 5000 bước tích phân thời gian (bước h = 0.01), phương pháp Điểm giữa mất tới 1248 giây (tương đương 20.8 phút), trong khi Störmer-Verlet chỉ tiêu tốn 0.31 giây và Euler symplectic tốn 0.27 giây. Như vậy, Störmer-Verlet rút ngắn thời gian tính toán tới 4025 lần, tương đương tiết kiệm 99.97% chi phí CPU.

Thứ tư, tính ổn định vượt bậc ở các mô hình phi tuyến phức tạp. Tại bài toán sine-Gordon (500 nút, 1000 bước), Störmer-Verlet hoàn thành trong 0.42 giây so với 874 giây của Điểm giữa (nhanh hơn xấp xỉ 2080 lần). Đối với phương trình Schrödinger (256 nút, 3000 bước), Störmer-Verlet chỉ mất 0.52 giây so với 306 giây của Điểm giữa.

Thảo luận kết quả

Nguyên nhân cốt lõi dẫn đến sự chênh lệch thời gian khổng lồ giữa hai phương pháp bậc hai bắt nguồn từ cấu trúc giải thuật. Phương pháp Điểm giữa là một lược đồ ẩn hoàn toàn trên cả hai biến tọa độ và xung lượng. Tại mỗi bước lặp, thuật toán bắt buộc phải giải một hệ phương trình đại số phi tuyến có số chiều gấp đôi số biến không gian rời rạc bằng các phép lặp giải tích ma trận tốn kém. Ngược lại, phương pháp Störmer-Verlet kết hợp khéo léo hai nửa bước cập nhật theo nguyên lý Euler symplectic đối xứng. Khi áp dụng cho các hệ Hamilton có hàm năng lượng tách biệt dạng động năng cộng thế năng, Störmer-Verlet trở thành một lược đồ hiện hoàn toàn, loại bỏ 100% nhu cầu giải phương trình phi tuyến ở mỗi bước thời gian.

So sánh với các nghiên cứu kinh điển của Ernst Hairer, Christian Lubich và Gerhard Wanner (2006) hay Benedict Leimkuhler và Sebastian Reich (2004), kết quả của luận văn đã kiểm chứng hoàn hảo định lý về sự bảo tồn năng lượng của các phương pháp symplectic: sai số năng lượng dao động điều hòa quanh nghiệm thực chứ không tích lũy tăng dần theo thời gian như các họ Runge-Kutta cổ điển.

Toàn bộ dữ liệu nghiên cứu có thể được thể hiện một cách trực quan thông qua hai dạng đồ thị khoa học:

• Dạng thứ nhất là biểu đồ đường biểu diễn độ lệch tương đối của hàm Hamilton theo trục thời gian, trong đó đường cong của Euler hiện dốc đứng tăng theo hàm mũ, trong khi các đường cong của Điểm giữa và Störmer-Verlet duy trì dạng dải sóng nằm ngang phẳng tuyệt đối.
• Dạng thứ hai là đồ thị không gian pha minh họa tính bảo toàn diện tích có hướng. Khi khởi tạo một hình tròn bán kính 0.4 và một hình elip bán trục 0.8 trên mặt phẳng pha trong khoảng thời gian từ 0 đến 3pi/2, dòng Hamilton làm các hình này bị kéo dãn và uốn cong phức tạp qua 4 mốc thời gian, nhưng tổng diện tích hình học của chúng được bảo tồn trọn vẹn không hề thay đổi.

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

Dựa trên các kết quả nghiên cứu và dữ liệu thực nghiệm, tác giả đề xuất 4 nhóm khuyến nghị mang tính ứng dụng cao:

Thứ nhất, ưu tiên tích hợp lược đồ Störmer-Verlet làm thuật toán lõi trong các phần mềm tính toán động lực học và mô phỏng vật lý.

• Hành động cụ thể: Thay thế các bộ giải số Euler và Runge-Kutta bậc thấp bằng thuật toán Störmer-Verlet trong các module mô phỏng cơ học.
• Chỉ số mục tiêu: Cắt giảm từ 90% đến 99% thời gian CPU khi chạy các bài toán dài hạn, đồng thời duy trì độ lệch năng lượng dưới 0.01%.
• Khung thời gian: Thực hiện chuyển giao và cài đặt trong vòng 6 tháng.
• Chủ thể thực hiện: Các kỹ sư lập trình công cụ mô phỏng vật lý, trung tâm tính toán hiệu năng cao và các nhà phát triển phần mềm mô hình hóa kỹ thuật.

Thứ hai, chuẩn hóa phương pháp rời rạc hóa sai phân trung tâm bảo toàn cấu trúc cho các phương trình đạo hàm riêng.

• Hành động cụ thể: Thiết lập thư viện ma trận sai phân hữu hạn đối xứng D theo biến không gian để bảo đảm tính chất Hamilton vô hạn chiều được kế thừa chính xác trên lưới hữu hạn.
• Chỉ số mục tiêu: Đạt 100% tính bảo toàn dạng symplectic trên không gian pha sau rời rạc hóa.
• Khung thời gian: Nghiên cứu và chuẩn hóa trong 9 tháng.
• Chủ thể thực hiện: Các nhóm nghiên cứu giải tích số và mô hình hóa toán học tại các trường đại học chuyên ngành.

Thứ ba, phát triển các thuật toán symplectic thích nghi cỡ bước thời gian.

• Hành động cụ thể: Nghiên cứu kết hợp biến đổi thời gian giải tích kiểu Poincaré để thay đổi bước lưới h linh hoạt tại các vùng có độ cong lớn mà không phá hủy tính bảo toàn cấu trúc pha.
• Chỉ số mục tiêu: Tăng tốc độ tính toán thêm 25% đến 35% trên các quỹ đạo có độ lệch tâm cao.
• Khung thời gian: Triển khai đề tài nghiên cứu giai đoạn 12 tháng.
• Chủ thể thực hiện: Giảng viên, nghiên cứu sinh chuyên ngành Toán ứng dụng.

Thứ tư, mở rộng ứng dụng cho các mô hình quang học phi tuyến và cơ học lượng tử phức tạp.

• Hành động cụ thể: Áp dụng phương pháp symplectic bậc hai cho phương trình Schrödinger phi tuyến 2 chiều và 3 chiều trong tính toán chất ngưng tụ Bose-Einstein.
• Chỉ số mục tiêu: Kiểm soát sai số chuẩn sóng bảo toàn dưới 0.001% qua hàng chục nghìn chu kỳ dao động.
• Khung thời gian: Lộ trình nghiên cứu ứng dụng từ năm 2022 đến năm 2025.
• Chủ thể thực hiện: Các phòng thí nghiệm vật lý lý thuyết và viện nghiên cứu vật liệu tiên tiến.

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

Luận văn thạc sĩ này là tài liệu chuyên khảo đặc biệt hữu ích cho 4 nhóm đối tượng chính:

Nhóm thứ nhất: Học viên cao học và nghiên cứu sinh chuyên ngành Toán học, Toán ứng dụng và Khoa học tính toán. Luận văn cung cấp nền tảng bài bản từ cơ sở vi phân của dòng Hamilton đến kỹ thuật phân tích sai số ngược (backward error analysis), giúp người học có cơ sở lý thuyết chuẩn xác để phát triển các đề tài luận án chuyên sâu.

Nhóm thứ hai: Kỹ sư mô phỏng cơ học thiên văn, hàng không vũ trụ và động lực học phân tử. Những người trực tiếp giải quyết bài toán tính toán quỹ đạo hành tinh, vệ tinh hay tương tác hạt có thể ứng dụng ngay công thức sai phân Störmer-Verlet để giải quyết triệt để vấn đề mất ổn định năng lượng sau hàng triệu chu kỳ lặp.

Nhóm thứ ba: Lập trình viên phát triển công cụ vật lý (Physics Engine) trong đồ họa máy tính và trò chơi điện tử. Luận văn cung cấp giải pháp thực tế giúp tái hiện chân thực các chuyển động dây treo, con lắc và sóng âm với tốc độ thời gian thực (real-time) mà không lo ngại hiện tượng nổ mô phỏng (simulation blow-up).

Nhóm thứ tư: Giảng viên đại học và các nhà nghiên cứu trong lĩnh vực Phương trình Vật lý Toán. Luận văn là nguồn tư liệu giảng dạy mẫu mực với đầy đủ phương pháp luận, bảng so sánh tham số thời gian và các ví dụ minh họa trực quan từ sóng đến phương trình Schrödinger trên mạng lưới 256 nút.

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

Thứ nhất: Phương pháp symplectic khác biệt bản chất như thế nào so với các phương pháp giải số truyền thống?

• Trả lời: Điểm khác biệt mấu chốt nằm ở tính bảo toàn cấu trúc hình học. Các phương pháp cổ điển như Euler hay Runge-Kutta chỉ tập trung xấp xỉ quỹ đạo nghiệm địa phương nên làm tích lũy sai số năng lượng khiến hệ thống bị trôi dạt dần theo thời gian. Ngược lại, phương pháp symplectic bảo toàn chính xác diện tích pha và dạng vi phân ngoại, giúp năng lượng của hệ chỉ dao động trong một ngưỡng cố định không đổi kể cả khi mô phỏng tới 5000 bước.

Thứ hai: Tại sao phương pháp Störmer-Verlet lại có tốc độ tính toán nhanh hơn phương pháp Điểm giữa hàng nghìn lần?

• Trả lời: Mặc dù cùng đạt độ chính xác bậc 2, Störmer-Verlet có cấu trúc giải thuật tường minh khi động năng và thế năng tách biệt. Phương pháp Điểm giữa là lược đồ ẩn, đòi hỏi giải hệ phương trình phi tuyến ở mỗi bước thời gian khiến bài toán sóng mất tới 1248 giây. Trong khi đó, Störmer-Verlet chỉ cập nhật tuần tự các biểu thức hiện nên hoàn thành bài toán tương tự chỉ sau 0.31 giây.

Thứ ba: Sai số năng lượng trong các phương pháp symplectic có bị tích tụ vô hạn khi thời gian kéo dài không?

• Trả lời: Hoàn toàn không. Theo định lý của Hairer và cộng sự, năng lượng tính toán bởi phương pháp symplectic tương ứng chính xác với một hàm Hamilton cải biên nằm sát hàm năng lượng gốc ở mức O(h mũ 2). Do đó, sai số năng lượng chỉ dao động tuần hoàn quanh giá trị ban đầu và bị chặn trong suốt khoảng thời gian mô phỏng theo hàm mũ của nghịch đảo bước thời gian.

Thứ tư: Khi nào thì phương pháp Störmer-Verlet có thể áp dụng dưới dạng một lược đồ hiện hoàn toàn?

• Trả lời: Lược đồ Störmer-Verlet trở thành phương pháp hiện khi hàm năng lượng Hamilton có dạng phân hoạch tách biến, nghĩa là đạo hàm riêng của H theo p chỉ phụ thuộc vào p, và đạo hàm riêng của H theo q chỉ phụ thuộc vào q. Đây là đặc điểm cố hữu trong hầu hết các hệ vật lý thực tế như phương trình sóng hay con lắc dao động.

Thứ năm: Phương trình Schrödinger phi tuyến phức tạp được đưa về hệ Hamilton hữu hạn chiều bằng cách nào?

• Trả lời: Tác giả tách hàm sóng phức u thành hai thành phần thực và ảo q và p, biến phương trình Schrödinger thành một hệ hai phương trình vi phân thực tương hỗ. Tiếp đó, phương pháp sai phân trung tâm trên mạng 256 nút lưới không gian chuyển hóa hệ vi phân thành hệ Hamilton kích thước 512 biến số hoàn chỉnh, bảo tồn cấu trúc Poisson chính tắc.

Kết luận

Luận văn thạc sĩ "Một số phương pháp symplectic cho hệ Hamilton" của tác giả Trần Văn Cẩn đã hoàn thành xuất sắc các mục tiêu nghiên cứu với 5 đóng góp học thuật và thực tiễn nổi bật:

• Hệ thống hóa toàn diện cơ sở toán học của hệ động lực Hamilton, làm sáng tỏ tính chất bất biến của năng lượng toàn phần và đặc trưng bảo toàn diện tích của dòng symplectic.
• Xây dựng và chứng minh toán học chặt chẽ tính symplectic cùng bậc chính xác cho ba phương pháp sai phân: Euler symplectic bậc 1, Điểm giữa bậc 2 và Störmer-Verlet bậc 2.
• Rời rạc hóa thành công ba bài toán phương trình đạo hàm riêng phi tuyến kinh điển trong vật lý toán gồm phương trình sóng, sine-Gordon và Schrödinger về dạng hệ Hamilton hữu hạn chiều.
• Thực nghiệm số chứng minh sự vượt trội toàn diện của phương pháp Störmer-Verlet, đạt hiệu năng tính toán nhanh hơn phương pháp Điểm giữa tới 4025 lần trên cùng một mức độ chính xác bậc 2.
• Cung cấp các kiểm chứng hình ảnh trực quan khẳng định tính bảo toàn hình học của pha và độ ổn định lâu dài của năng lượng qua hàng nghìn bước lặp.

Về lộ trình phát triển tiếp theo, nghiên cứu mở ra hướng đi đầy tiềm năng trong việc xây dựng các lược đồ symplectic bậc cao (bậc 4 và bậc 6) kết hợp kỹ thuật thích nghi bước thời gian cho các bài toán cơ học thiên văn đa vật thể trong giai đoạn 2022 - 2025. Để tiếp cận toàn văn công trình, các chứng minh định lý chi tiết và dữ liệu mô phỏng, quý độc giả và các nhà khoa học có thể tham khảo trực tiếp luận văn tại thư viện Trường Đại học Khoa học - Đại học Thái Nguyên.

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

Chương 1 Phương pháp số cho phương trình vi phân tổng quát Chương này trình bày một số phương pháp thông dụng cho phương trình vi phân thường tổng quát. Các phương pháp được nhắc đến chủ yếu là Euler hiện, Euler ẩn, Điểm giữa, phương pháp hình thang và họ phương pháp Runge-Kutta. Thêm vào đó, chúng tôi cũng trình bày chi tiết khái niệm về bậc chính xác, tính ổn định và sự hội tụ của một lược đồ sai phân. Tài liệu tham khảo chính cho chương này là cuốn sách chuyên khảo [3].1 Khái niệm và phương pháp cơ bản Ta xét bài toán giá trị ban đầu cho phương trình vi phân thường dạng tổng quát ẏ = f (t, y), t ∈ [0, T ], y(0) = y0 , (1.1) trong đó ẩn hàm y và hàm vế phải f (t, y) là các đại lượng véc tơ trong Rn.

Ta giả sử rằng f (t, y) là hàm đủ trơn và bị chặn để đảm bảo rằng bài toán (1.1) luôn tồn tại nghiệm, và nghiệm đó khả vi tới mức cần thiết để phục vụ cho các lập luận lý thuyết. Để giải bài toán trên máy tính, người ta phải rời rạc hóa miền xác định của ẩn hàm. Xét lưới chia 0 = t0 < t1 < .2) 4 và đặt các cỡ bước thứ k là hk := tk − tk−1. Ta sẽ lần lượt xây dựng các xấp xỉ y0 , y1 ,.

, yK−1 , yK của các giá trị ẩn hàm cần tìm y(tk ), k = 0,. Do đã biết ẩn hàm tại t0 , dựa vào ràng buộc phương trình, ta có thể tìm xấp xỉ tại t1. Điều này cứ tiếp tục cho đến khi toàn bộ xấp xỉ được tìm. Người ta còn gọi các bước tìm này là các bước cập nhật do nó cho ta biết trạng mới dựa trên thông tin của trạng thái cũ.

Những phương pháp này được xếp vào loại phương pháp một bước. Một phương pháp đa bước là phương pháp mà khi tính yk , nó đòi hỏi thông tin về giá trị của hàm tại ít nhất hai bước trước đó yk−1 , yk−2 ,. Luận văn này chỉ xét các phương pháp một bước. Người ta cũng gọi tương ứng từ trạng thái cũ với trạng thái mới là ánh xạ bước và ký hiệu là φhk (y).

Để thực hiện điều này, ta xét khai triển Taylor của ẩn hàm tại điểm lưới tk bất kỳ với giả thiết rằng ta đã biết giá trị của ẩn hàm tại điểm trước đó tk−1 1 y(tk ) = y(tk−1 ) + hk ẏ(tk−1 ) + h2k ÿ(tk−1 ) + · · · .3) 2 Ta sẽ thường xuyên sử dụng khái niệm vô cùng bé cùng bậc, ký hiệu là O(hm ). Ta nói đại lượng y là một vô cùng bé bậc m của h, hay y tương đương với hm , nếu tồn tại hằng số C > 0 sao cho |y| ≤ Chm , với h > 0 đủ nhỏ. Sử dụng ký hiệu này, ta có thể viết lại phương trình trên ở dạng y(tk ) = y(tk−1 ) + hk ẏ(tk−1 ) + O(h2k ).4) Phương pháp Euler hiện (explicit Euler) hay phương pháp Euler tiến (for- ward Euler) chính là hệ quả của việc bỏ đi đại lượng vô cùng bé bậc hai trong 5 (1.4) và thay ẏ bởi f để thu được yk = yk−1 + hk f (tk−1 , yk−1 ).5) Đây là phương pháp số đơn giản nhất để giải bài toán giá trị ban đầu cho phương trình vi phân. Nó là một phương pháp hiện vì để tìm nghiệm xấp xỉ tại k , ta chỉ cần sử dụng một công thức hiện và thông tin nghiệm tại bước trước đó k − 1.

Phương pháp Euler ẩn (implicit Euler) hay phương pháp Euler lùi (back- ward Euler) được xây dựng tương tự như phương pháp Euler hiện nhưng thay vì dựa vào thông tin tại tk−1 để tính toán tại tk , ta lại thực hiện ngược lại, tức là y(tk−1 ) = y(tk ) − hk ẏ(tk ) + O(h2k ).6) Từ đây, ta có lược đồ ẩn yk = yk−1 + hk f (tk , yk ).7) Có thể quan sát thấy, (1.7) không đơn giản như (1.5) vì giá trị cần tìm lại nằm trong hàm f , Để tìm yk từ yk−1 , ta phải giải phương trình đại số (1. Nếu f (t, y) là một hàm tuyến tính thì (1.7) là một hệ phương trình tuyến tính. Ngược lại, ta đối mặt với một hệ phi tuyến mà muốn giải được, ta phải dùng một phương pháp riêng, chẳng hạn Newton. Điều này lý giải thuật ngữ “ẩn" trong tên gọi của nó.

Câu hỏi đặt ra ở đây là tại sao Euler ẩn khó và phức tạp hơn Euler hiện nhưng người ta vẫn sử dụng. Lí do của việc này là trong một số lớp bài toán cương, việc giải bằng phương pháp hiện thường không chính xác do tính dao động rất mạnh quanh nghiệm giải tích trừ khi bước thời gian cực kỳ nhỏ. Nhưng việc giải với bước thời gian nhỏ sẽ gặp bất lợi nếu ta cần xét phương trình vi phân trên một khoảng thời gian dài. Khi đó, người ta thường sử dụng các phương pháp ẩn.

Cả hai phương pháp Euler đều sử dụng các xấp xỉ bậc một nên đôi khi gặp bất lợi cho việc áp dụng cho những mô hình tính toán đòi hỏi tính chính xác cao. 6 Điều này có thể được cải thiện khi sử dụng phương pháp hình thang. Ý tưởng của phương pháp này là “mượn" thông tin của ẩn hàm ở điểm giữa tk−1/2 , trung điểm của tk−1 và tk , để xây dựng một lược đồ xấp xỉ bậc hai. Cụ thể, ta viết hk h2 h3 .8) Trừ hai phương trình trên vế theo vế rồi chia cho hk , ta thu được y(tk ) − y(tk−1 ) h2k .9) hk 24 Tiếp đến, áp dụng các khai triển giống như ở (1.8) nhưng cho ẏ rồi cộng hai phương trình lại với nhau, ta thu được h2k .10), ta thu được y(tk ) − y(tk−1 ) ẏ(tk ) + ẏ(tk−1 ) h2k .11) hk 2 12 Mối quan hệ (1.11) gợi ý cho ta công thức xấp xỉ hk yk = yk−1 + (f (tk , yk ) + f (tk−1 , yk−1 )).12) còn được gọi là phương pháp hình thang.

Có thể thấy đây cũng là một lược đồ ẩn.2 Sự hội tụ, tính chính xác, và tính ổn định của lược đồ sai phân Ta có thể nhận thấy có nhiều sai số trong việc giải xấp xỉ một phương trình vi phân. Trước tiên, tại mỗi bước, người ta xấp xỉ đạo hàm bởi sai phân nên sinh 7 ra sai số. Tiếp đến, người ta lại sử dụng kết quả tính được, vốn chứa sai số, để tính ở bước tiếp theo. Như vậy, nhìn chung, các sai số sẽ tích lũy trong quá trình giải và có xu hướng càng về sau càng lớn.

Hầu hết các phương pháp số đều phải quan tâm đến hai vấn đề này. Để thuận tiện cho phân tích, ta sẽ trình bày những phân tích này cho phương pháp Euler hiện. Ý tưởng và phương pháp luận của nó được áp dụng tương tự cho các lược đồ khác. Ký hiệu toán tử sai phân tương ứng với lược đồ Euler hiện và cỡ bước h áp dụng cho phương trình (1.1) với hàm u bất kỳ định nghĩa trên lưới chia (1.13) hk Giả sử yh là một hàm lưới có giá trị yk tại mỗi điểm lưới tk , k = 0,.

Khi đó, lược đồ cập nhật tương ứng có thể được viết ngắn gọn ở dạng Nh yh (tk ) = 0, y0 cho trước. Khi đó, ta định nghĩa sai số làm gọn địa phương, do việc cắt bỏ các đại lượng vô cùng bé, là phần dư khi áp dụng toán tử sai phân vào nghiệm chính xác, tức là dk = Nh y(tk ).14) Có thể thấy, sai số làm gọn địa phương cho ta biết sai khác giữa toán tử sai phân và toán tử vi phân. Một phương pháp sai phân được gọi là chính xác bậc p nếu dk = O(hpk ), (1.15) với p là một số dương cho trước. Chẳng hạn, từ khai triển Taylor (1.3), đối với phương pháp Euler hiện, ta có hk dk = ÿ(tk ) + O(h2k ).16) 2 8 Do đó, Euler hiện có độ chính xác bậc một.

Phương pháp Euler ẩn cũng có độ chính xác bậc một. Bằng phân tích tương tự, ta có thể chỉ ra được rằng phương pháp hình thang có độ chính xác bậc hai. Tuy nhiên, phân tích sai số sẽ không có ý nghĩa nếu ta chỉ ước lượng được sai số làm gọn địa phương mà không nói được sai số của nghiệm thu được trên toàn lưới sai phân. Để đạt được điều này, trước tiên ta ký hiệu cỡ bước lớn nhất h = max{hk , k = 1,.

, K} và giả sử đại lượng Kh là bị chặn mà không phụ thuộc vào số điểm của lưới sai phân K. Ta ký hiệu ek = yk − y(tk ), k = 0,. , K và gọi nó là sai số toàn cục. Phương pháp sai phân được gọi là hội tụ bậc p nếu sai số toàn cục thỏa mãn ek = O(hp ), k = 0,.

Nhìn chung, bậc của hội tụ thường thấp hơn bậc của sai số làm gọn địa phương. Điều kiện sau đây đảm bảo bậc hội tụ bằng bậc của sai số làm gọn địa phương. Phương pháp sai phân được nói là 0-ổn định (đọc là “ổn định bậc không") nếu tồn tại các hằng số h0 và C sao cho với mọi hàm lưới xh và zh , h ≤ h0 , bất đẳng thức sau được thỏa mãn |xk −zk | ≤ C(|x0 −z0 |+ max |Nh xh (tj )−Nh zh (tj ))|, 1 ≤ k ≤ K.17) 1≤j≤K Thảo luận ở trên được chính xác hóa trong khẳng định sau đây. Nếu một phương pháp sai phân có độ chính xác bậc p và 0-ổn định thì nó hội tụ bậc p.

Cụ thể, ta có |e| ≤ C max |dj | = O(hp ).18) j Chứng minh chi tiết định lý này có thể được tìm thấy ở [3, Định lý 3. 9 Ta đã biết phương pháp Euler hiện có độ chính xác bậc 1, áp dụng Định lý 1.1, ta sẽ chỉ ra rằng phương pháp này hội tụ bậc 1 bằng cách chứng minh nó 0-ổn định. Thật vậy, Ký hiệu sk = xk − zk , θ = max |Nh xh (tj ) − Nh zh (tj ))|. 1≤j≤K Với mỗi k , ta có sk − sk−1 θ≥ − (f (tk−1 , xk−1 ) − f (tk−1 , zk−1 )) hk |sk | sk−1 ≥ − + (f (tk−1 , xk−1 ) − f (tk−1 , zk−1 )).

hk hk Sử dụng điều kiện liên tục Lipschitz của f , ta thu được sk−1 |sk−1 | + (f (tk−1 , xk−1 ) − f (tk−1 , zk−1 )) ≤ + L|sk−1 | hk hk   1 ≤ + L |sk−1 |.