Bối cảnh và vấn đề nghiên cứu
Trong vật lý lý thuyết và vật lý toán, các mô hình toán học giải tích mô tả hiện tượng tự nhiên thường bắt đầu từ việc khảo sát các bài toán đơn giản hóa nhằm tìm kiếm nghiệm chính xác dưới dạng giải tích. Các nghiệm giải tích này đóng vai trò làm hệ hàm cơ sở và phép thử nghiệm chuẩn xác để đánh giá, lựa chọn các thuật toán số học khi giải quyết các bài toán vật lý phức tạp hơn. Nhiều bài toán quan trọng trong vật lý liên quan đến quá trình truyền nhiệt, tương tác bức xạ với vật chất, sự lan truyền của sóng điện từ, sóng âm, cũng như cấu trúc bên trong của các sao và phản ứng hạt nhân đều dẫn về bài toán Sturm-Liouville chứa phương trình vi phân đạo hàm riêng Laplace hoặc Helmholtz.
Đối với các miền không gian hình học đơn giản như đoạn thẳng, hình chữ nhật hoặc hình bình hành, nghiệm của phương trình vi phân có thể biểu diễn qua các hàm sơ cấp cơ bản. Tuy nhiên, khi khảo sát các miền có tính đối xứng cong như hình tròn, hình trụ, hình cầu hoặc các mặt bậc hai phức tạp, hệ hàm riêng bắt buộc phải được biểu diễn thông qua các hàm toán đặc biệt. Việc nghiên cứu có hệ thống và ứng dụng các hàm này là nhiệm vụ thiết yếu để giải quyết bài toán biên trong các ngành khoa học tự nhiên.
Khóa luận tốt nghiệp của tác giả Vũ Hoàng Thanh Trang xác định các mục tiêu và nhiệm vụ nghiên cứu cụ thể như sau:
- Hệ thống hóa định nghĩa, tính chất giải tích, hệ thức truy hồi, tính trực giao và dạng tiệm cận của một số hàm toán đặc biệt thường dùng: hàm Bessel (cùng họ hàm trụ tổng quát), đa thức trực giao Legendre cổ điển, đa thức Legendre liên hiệp và hệ hàm cầu.
- Ứng dụng các hàm toán đặc biệt đã khảo sát vào việc thiết lập và giải chính xác nghiệm của 3 bài toán biên vật lý điển hình:
- Bài toán truyền nhiệt làm nguội của hình trụ tròn dài vô hạn.
- Bài toán dao động và sóng dừng trên bề mặt màng trống tròn đàn hồi.
- Bài toán tán xạ sóng vô hướng có tính đối xứng trục trên vật thể phỏng cầu dài (prolate spheroid).
- Lập trình thuật toán giải tích và mô phỏng số trên phần mềm máy tính để trực quan hóa sự phân bố nhiệt độ, các mode dao động màng và trường sóng tán xạ.
Đối tượng và phạm vi nghiên cứu của đề tài:
- Đối tượng nghiên cứu: Các hàm toán đặc biệt (Bessel loại một $J_\nu$, Bessel loại hai $Y_\nu$, Hankel $H_\nu^{(1)}, H_\nu^{(2)}$, Neumann $N_\nu$, Infeld $I_\nu$, Macdonald $K_\nu$, đa thức Legendre $P_n$, đa thức Legendre liên hiệp $P_n^m$, hàm cầu $Y_n^m$) và các phương trình vi phân đạo hàm riêng của vật lý toán (phương trình truyền nhiệt, phương trình sóng hai chiều, phương trình Laplace và phương trình Helmholtz).
- Phạm vi nghiên cứu: Khóa luận tốt nghiệp bậc đại học ngành Sư phạm Vật lý (thuộc tổ bộ môn Toán lý, Khoa Vật lý, Trường Đại học Sư phạm TP. Hồ Chí Minh, năm 2020), dung lượng toàn văn gồm 49 trang, 15 hình vẽ, 3 bảng số liệu tính toán và 3 phụ lục chương trình.
Cơ sở lý thuyết và phương pháp
Khung lý thuyết và mô hình toán học
Đề tài vận dụng lý thuyết bài toán biên Sturm-Liouville và phương trình vi phân tuyến tính cấp hai có điểm kỳ dị:
- Phương trình vi phân Bessel: $x^2 u'' + xu' + (x^2 - \nu^2)u = 0$, xuất phát từ việc giải phương trình Laplace hoặc Helmholtz trong hệ tọa độ cực và tọa độ trụ.
- Phương trình vi phân Legendre: $\frac{d}{dx}\left[(1-x^2)\frac{dy}{dx}\right] + \lambda y = 0$ trên đoạn $[-1, 1]$, với trị riêng $\lambda = n(n+1)$.
- Phương trình vi phân Legendre liên hiệp: $\frac{d}{dx}\left[(1-x^2)\frac{dP_n^m}{dx}\right] + \left[n(n+1) - \frac{m^2}{1-x^2}\right]P_n^m(x) = 0$, là cơ sở trực tiếp để xây dựng phần hàm góc của toán tử Laplace trên mặt cầu.
- Lý thuyết hàm cầu và Định lý Steklov: Khai triển chuỗi Laplace của một hàm khả tích tuyệt đối và biến thiên hữu hạn trên mặt cầu đơn vị theo hệ hàm cầu trực giao đầy đủ $Y_n^m(\theta, \phi)$.
- Phương trình phỏng cầu dài: Mô hình hình học ellipsoid quay quanh trục lớn $Oz$ trong tọa độ cầu $r(\theta) = \sqrt{1 + \varepsilon^2 \cos^2\theta}$, với $\varepsilon$ là tham số biến dạng hình học biểu thị tỉ số giữa trục lớn và trục nhỏ.
Phương pháp nghiên cứu
- Phương pháp phân tích giải tích: Sử dụng phương pháp tách biến Fourier để phân ly phương trình vi phân đạo hàm riêng nhiều chiều thành hệ các phương trình vi phân thường đối với từng tọa độ độc lập (bán kính, góc cực, góc kinh độ, thời gian).
- Phương pháp chuỗi lũy thừa và công thức Frobenius: Tìm nghiệm giải tích của phương trình vi phân quanh điểm kỳ dị, thiết lập các hệ thức truy hồi cho hệ số khai triển, từ đó xây dựng biểu thức giải tích của hàm Bessel và đa thức Legendre.
- Phương pháp giải tích phức và khai triển tiệm cận: Thiết lập các công thức xấp xỉ cho hàm Bessel, hàm Hankel khi đối số $x \to \infty$ và $x \to 0$.
- Phương pháp tính toán số và mô phỏng đồ họa: Tác giả sử dụng phần mềm toán học Maple để lập trình tính toán các nghiệm số (trị riêng, điểm không của hàm Bessel), tính ma trận hệ số tán xạ và vẽ đồ thị 2D, 3D của các trường vật lý.
- Nguồn dữ liệu: Nghiên cứu thuần túy lý thuyết và giải tích số; toàn bộ số liệu bắt nguồn từ các hằng số toán học, công thức vật lý chuẩn và kết quả giải số trực tiếp trên phần mềm Maple.
Nội dung chính theo từng chương
Chương 1: Hệ thống một số hàm toán đặc biệt
Chương này trình bày cơ sở giải tích của các lớp hàm toán đặc biệt cốt lõi:
-
Hàm Bessel và các hàm trụ:
- Phương trình Bessel tổng quát được giải bằng phương pháp chuỗi lũy thừa $u(x) = \sum_{k=0}^\infty a_k x^{\rho+k}$. Khảo sát hệ thức truy hồi dẫn đến định nghĩa hàm Bessel loại một cấp $\nu$: $$J_\nu(x) = \sum_{k=0}^\infty \frac{(-1)^k}{k!\Gamma(\nu+k+1)} \left(\frac{x}{2}\right)^{2k+\nu}$$
- Hàm Bessel loại hai $Y_\nu(x)$ được định nghĩa để tạo thành hệ nghiệm độc lập tuyến tính với $J_\nu(x)$ khi $\nu$ không nguyên; khi $\nu=n$ nguyên thì $Y_n(x) = (-1)^n J_{-n}(x)$.
- Hệ thức truy hồi liên hệ các cấp lân cận: $J_{\nu+1}(x) = \frac{2\nu}{x}J_\nu(x) - J_{\nu-1}(x)$.
- Các hàm Bessel cấp bán nguyên được biểu diễn qua hàm sơ cấp: $J_{1/2}(x) = \sqrt{\frac{2}{\pi x}}\sin x$ và $Y_{1/2}(x) = -\sqrt{\frac{2}{\pi x}}\cos x$.
- Tác giả tính toán chính xác 3 điểm không đầu tiên của $J_0(x)$ đến 4 chữ số thập phân: $\mu_1^{(0)} = 2.4048$, $\mu_2^{(0)} = 5.5201$, $\mu_3^{(0)} = 8.6537$.
- Khảo sát các hàm trụ mở rộng: hàm Hankel loại một $H_\nu^{(1)}(x) = J_\nu(x) + iN_\nu(x)$, Hankel loại hai $H_\nu^{(2)}(x) = J_\nu(x) - iN_\nu(x)$, hàm Neumann $N_\nu(x)$, hàm Bessel biến tính Infeld $I_\nu(x)$ và Macdonald $K_\nu(x)$.
-
Đa thức trực giao cổ điển Legendre:
- Nghiệm của phương trình vi phân Legendre trên đoạn $[-1, 1]$ trở thành đa thức khi tham số $\lambda = n(n+1)$ ($n \in \mathbb{N}$).
- Đa thức Legendre $P_n(x)$ thỏa mãn tính trực giao $\int_{-1}^1 P_n(x)P_m(x)dx = 0$ khi $n \ne m$ với trọng số $p(x)=1$.
- Biểu diễn thông qua công thức vi phân Rodrigues: $P_n(x) = \frac{1}{2^n n!} \frac{d^n}{dx^n}(x^2-1)^n$.
- Tác giả lập biểu thức giải tích tường minh cho 11 đa thức đầu tiên từ $P_0(x)$ đến $P_{10}(x)$.
-
Đa thức Legendre liên hiệp:
- Định nghĩa thông qua đạo hàm cấp $m$ của đa thức Legendre: $P_n^m(x) = (-1)^m (1-x^2)^{m/2} \frac{d^m}{dx^m}(P_n(x))$.
- Tính trực giao khi cố định bậc $m$: $\int_{-1}^1 P_k^m(x) P_n^m(x) dx = \frac{2(n+m)!}{(2n+1)(n-m)!}\delta_{k,n}$.
- Liệt kê biểu thức đại số và đồ thị của các hàm $P_n^m(x)$ với $m = 1, 2, 3$ và $n = 0, 1, 2, 3$.
-
Định nghĩa và tính chất cơ bản của hàm cầu:
- Tách biến phương trình Laplace trong tọa độ cầu dẫn đến phương trình góc $\Delta_{\theta,\phi}Y + n(n+1)Y = 0$.
- Hàm cầu được xác định bởi: $Y_n^m(\theta, \phi) = P_n^{|m|}(\cos\theta) e^{im\phi}$ hoặc viết dưới dạng hàm lượng giác thực $\cos(m\phi), \sin(m\phi)$.
- Thiết lập bảng biểu thức giải tích tường minh cho các hàm cầu từ cấp 0 đến cấp 2: $Y_0^{(0)}, Y_1^{(-1)}, Y_1^{(0)}, Y_1^{(1)}, Y_2^{(-2)}, Y_2^{(-1)}, Y_2^{(0)}, Y_2^{(1)}, Y_2^{(2)}$ cùng đồ thị 2D và 3D.
- Trình bày Định lý Steklov về khai triển chuỗi Laplace của hàm $f(\theta, \phi)$ trên mặt cầu.
-
Hàm riêng của quả cầu:
- Phân tích phương trình Euler cho hàm bán kính $r^2 R'' + 2rR' - n(n+1)R = 0$, xác định hệ nghiệm riêng tổng quát trong quả cầu $u_{nm}(r, \theta, \phi) = r^n Y_n^m(\theta, \phi)$ và ngoài quả cầu $u_{nm}(r, \theta, \phi) = r^{-(n+1)} Y_n^m(\theta, \phi)$.
Chương 2: Ứng dụng của các hàm toán đặc biệt trong việc giải các bài toán biên
Chương này triển khai giải 3 bài toán biên vật lý:
-
Bài toán về sự làm nguội của hình trụ tròn dài vô hạn:
- Khảo sát phương trình truyền nhiệt đối xứng trục $\frac{\partial u}{\partial t} = a^2 \Delta u$ trên hình trụ bán kính $r_0$, điều kiện biên Dirichlet $u(r_0, t) = 0$, nhiệt độ ban đầu $u(r, 0) = \phi(r)$.
- Bằng phương pháp tách biến $u(r,t) = R(r)T(t)$, phương trình không gian chuyển về phương trình Bessel cấp 0: $x^2 R'' + x R' + x^2 R = 0$ (với $x = \lambda r$).
- Điều kiện hữu hạn tại tâm $|R(0)| < \infty$ loại bỏ nghiệm Neumann $N_0$, giữ lại nghiệm Bessel loại một $J_0(\lambda r)$.
- Điều kiện biên dẫn đến phương trình trị riêng $J_0(\lambda r_0) = 0 \implies \lambda_k = \frac{\mu_k^{(0)}}{r_0}$.
- Nghiệm tổng quát có dạng chuỗi Fourier-Bessel: $$u(r,t) = \sum_{k=1}^\infty A_k J_0\left(\frac{\mu_k^{(0)}}{r_0}r\right) \exp\left[-\left(\frac{\mu_k^{(0)}}{r_0}a\right)^2 t\right]$$
- Tác giả thực hiện tính toán số với bộ thông số $r_0 = 1$, $a = 1$, $\phi(r) = 100$, mô phỏng profile nhiệt độ tại các thời điểm $t = 0; 0.1; 0.2; 0.3; 0.4; 0.5$ (Hình 2.2; mã nguồn Maple tại Phụ lục A, trang 44).
-
Bài toán khảo sát sự rung động của bề mặt trống:
- Khảo sát màng tròn đàn hồi bán kính $r_0$, cố định tại biên $u(r_0, \phi, t) = 0$, tuân theo phương trình sóng hai chiều $\frac{\partial^2 u}{\partial t^2} = c^2 \Delta u$.
- Xét trường hợp đối xứng trục: nghiệm biểu diễn qua $J_0(\lambda_k r)\cos(c\lambda_k t)$ với $\lambda_k = \mu_k^{(0)}/r_0$.
- Xét trường hợp tổng quát phụ thuộc góc $\phi$: nghiệm hàm bán kính là $J_m(\lambda_k^{(m)} r)$ với $J_m(\lambda_k^{(m)} r_0) = 0$. Nghiệm tổng quát có dạng: $$u_k^{(m)}(r, \phi, t) = \left[A\cos(c\lambda_k^{(m)}t) + B\sin(c\lambda_k^{(m)}t)\right] J_m(\lambda_k^{(m)} r) \left[C\cos(m\phi) + D\sin(m\phi)\right]$$
- Đồ thị hóa các mode dao động tương ứng với các bộ chỉ số $(m,n)$ (Hình 2.5; mã nguồn Maple tại Phụ lục B, trang 45).
-
Bài toán tán xạ vô hướng trên phỏng cầu dài:
- Khảo sát phương trình Helmholtz $(\Delta + k^2)\psi = 0$ ngoài miền phỏng cầu $r > \sqrt{1 + \varepsilon^2}$.
- Trường sóng tới $\psi_0 = \exp(ikr\cos\theta)$ và sóng tán xạ được khai triển qua chuỗi hàm cầu và hàm Hankel: $$\psi = \sum_{n=0}^\infty \left[c_n \psi_n^-(r) + c_n^+ \psi_n^+(r)\right] Y_n(\theta)$$
- Tác giả tính toán giá trị biên độ sóng tới $c_n = i^n \sqrt{\frac{2}{2n+1}}$ và bình phương module $|c_n|^2$ (Bảng 2.1).
- Thiết lập hệ phương trình biên để giải các hệ số tán xạ $c_n^+$, tính đại lượng tiêu hao dòng sóng $S_N = \sum_{n=0}^N |c_n^+|^2$ tương ứng các độ lệch trục $\varepsilon = 0, 1, 2, 3$ cho trường hợp cắt cụt $N=4$ (Bảng 2.2) và $N=5$ (Bảng 2.3).
- Tính toán tiết diện tán xạ toàn phần $\sigma = \int_0^{2\pi}\int_0^\pi I(\theta)\sin\theta d\theta d\phi$, hàm tán xạ chuẩn hóa $T_\omega$, mật độ xác suất $|\psi(r,\theta)|^2$, phần thực $\text{Re}(\psi)$ và phần ảo $\text{Im}(\psi)$ (Hình 2.8, 2.9; mã nguồn Maple tại Phụ lục C, trang 46).
Kết quả và đóng góp
Số liệu và kết quả tính toán định lượng
Tác giả đã tính toán chính xác các hệ số sóng tới $c_n$ và năng lượng tương ứng $|c_n|^2$ trong bài toán tán xạ vô hướng:
| Bậc $n$ | Hệ số sóng tới $c_n$ | Bình phương module $|c_n|^2$ |
|---|---|---|
| 0 | 0.7071 | 0.5 |
| 1 | 1.2247i | 1.5 |
| 2 | 1.5811 | 2.5 |
| 3 | 1.8708i | 3.5 |
| 4 | 2.1213 | 4.5 |
| 5 | 2.3452i | 5.5 |
Bảng 2.1: Bảng giá trị của $c_n$ và $|c_n|^2$ ứng với $n = 0, 1, 2, 3, 4, 5$ (trích xuất từ trang 37 của khóa luận).
Tổng hợp các bài toán biên đã giải quyết trong khóa luận:
| Tên bài toán biên | Phương trình chi phối | Hệ tọa độ sử dụng | Lớp hàm đặc biệt xuất hiện trong nghiệm |
|---|---|---|---|
| Làm nguội hình trụ dài vô hạn | Truyền nhiệt $\frac{\partial u}{\partial t} = a^2 \Delta u$ | Tọa độ trụ $(r, t)$ | Hàm Bessel loại một cấp không $J_0$ |
| Rung động của màng trống tròn | Sóng 2D $\frac{\partial^2 u}{\partial t^2} = c^2 \Delta u$ | Tọa độ cực $(r, \phi, t)$ | Hàm Bessel loại một cấp nguyên $J_m$ |
| Tán xạ vô hướng trên phỏng cầu dài | Helmholtz $(\Delta + k^2)\psi = 0$ | Tọa độ cầu $(r, \theta)$ | Hàm Hankel cấp bán nguyên $H_{n+1/2}^{(1),(2)}$, Đa thức Legendre $P_n$ |
Đóng góp học thuật của đề tài
- Hệ thống hóa hoàn chỉnh lý thuyết giải tích từ các phương trình vi phân nền tảng đến các hệ thức biểu diễn tường minh của họ hàm trụ, hàm Legendre và hàm cầu.
- Thiết lập quy trình mẫu mực giải bài toán biên Sturm-Liouville bằng phương pháp tách biến trên các hình học đối xứng cong.
- Đưa ra giải pháp tính toán giải tích kết hợp số trị hoàn chỉnh cho bài toán tán xạ trên vật thể biến dạng hình học phỏng cầu dài (có xét đến tham số độ lệch trục $\varepsilon$).
- Xây dựng 3 bộ mã nguồn chương trình Maple (tại Phụ lục A, B, C) phục vụ tính toán số và trực quan hóa kết quả vật lý.
Hạn chế và hướng nghiên cứu tiếp
- Hạn chế: Khóa luận chỉ tập trung khảo sát 3 lớp hàm toán đặc biệt cơ bản (Bessel, Legendre, hàm cầu) áp dụng cho môi trường đồng nhất và hình học giải tích chuẩn (trụ tròn, màng phẳng tròn, phỏng cầu dài). Đề tài chưa xét đến môi trường dị hướng, điều kiện biên hỗn hợp phi tuyến hoặc bài toán tán xạ trường vector (như hệ phương trình Maxwell trong điện từ trường).
- Hướng phát triển: Tác giả kiến nghị mở rộng nghiên cứu sang các lớp hàm toán đặc biệt khác thường gặp trong vật lý lượng tử và vật lý hạt nhân nguyên tử, bao gồm: hàm Airy, đa thức Laguerre, đa thức Chebyshev, hàm Kelvin và đa thức Hermite.
Giá trị tham khảo
- Đối tượng thụ hưởng: Sinh viên, học viên cao học và giảng viên các ngành Sư phạm Vật lý, Vật lý học, Toán ứng dụng, Cơ học kỹ thuật, Kỹ thuật nhiệt và Kỹ thuật viễn thông.
- Phần nội dung đáng tham khảo:
- Quy trình khai triển nghiệm chuỗi Fourier-Bessel và chuỗi Laplace trong bài toán biên có điều kiện Dirichlet.
- Phương pháp xử lý bài toán tán xạ sóng phức tạp trên vật thể phỏng cầu dài.
- Các đoạn mã lập trình Maple mẫu tại Phụ lục A (trang 44), Phụ lục B (trang 45) và Phụ lục C (trang 46) giúp ích trực tiếp cho việc tính toán số trị và vẽ đồ thị hàm đặc biệt.
Câu hỏi thường gặp
1. Khóa luận tính toán các điểm không của hàm Bessel loại một $J_0(x)$ đạt độ chính xác bao nhiêu và giá trị cụ thể là gì?
Tác giả đã tính toán 3 điểm không thực dương đầu tiên của hàm $J_0(x)$ với độ chính xác đến 4 chữ số thập phân, lần lượt là: $\mu_1^{(0)} = 2.4048$, $\mu_2^{(0)} = 5.5201$, và $\mu_3^{(0)} = 8.6537$.
2. Trong mô phỏng bài toán làm nguội hình trụ dài vô hạn, các tham số đầu vào được thiết lập cụ thể như thế nào?
Các thông số được chọn để tính toán và vẽ đồ thị phân bố nhiệt độ gồm: bán kính hình trụ $r_0 = 1$, hệ số khuếch tán nhiệt $a = 1$, phân bố nhiệt ban đầu đồng đều $\phi(r) = 100$, và thời điểm khảo sát tại $t = 0; 0.1; 0.2; 0.3; 0.4; 0.5$.
3. Phương trình toán học mô tả hình học phỏng cầu dài (prolate spheroid) trong hệ tọa độ cầu được xác định ra sao?
Phương trình của phỏng cầu dài được biểu diễn dưới dạng: $r(\theta) = \sqrt{1 + \varepsilon^2 \cos^2\theta}$, trong đó trục lớn trùng với trục $Oz$, trục nhỏ được chuẩn hóa bằng 1 đơn vị, và $\varepsilon$ là tham số phụ thuộc vào tỉ số giữa trục lớn và trục nhỏ.
4. Khóa luận sử dụng phần mềm nào để lập trình tính toán và các đoạn mã nguồn được bố trí ở đâu?
Tác giả sử dụng phần mềm Maple để lập trình các thuật toán giải số và xuất đồ thị. Toàn bộ mã nguồn được đặt tại phần phụ lục cuối khóa luận: Phụ lục A (bài toán truyền nhiệt hình trụ, trang 44), Phụ lục B (bài toán dao động màng trống, trang 45) và Phụ lục C (bài toán tán xạ trên phỏng cầu dài, trang 46).
Kết luận
Khóa luận tốt nghiệp của tác giả Vũ Hoàng Thanh Trang đã hệ thống hóa lý thuyết giải tích về các hàm toán đặc biệt quan trọng gồm hàm Bessel, đa thức Legendre và hàm cầu. Nghiên cứu đã ứng dụng các lớp hàm này để giải quyết chính xác 3 bài toán biên vật lý thực tế về truyền nhiệt, dao động cơ học và tán xạ sóng vô hướng trên vật thể hình học phức tạp. Kết quả giải tích kết hợp cùng chương trình mô phỏng số trên Maple tạo nên một tài liệu tham khảo có giá trị học thuật cao cho việc nghiên cứu và giảng dạy môn Toán lý và Vật lý lý thuyết.