Bối cảnh và vấn đề nghiên cứu

Trong vật lý lý thuyết và ứng dụng, việc giải các phương trình đạo hàm riêng đóng vai trò then chốt trong việc mô tả định lượng các quá trình tự nhiên. Các phương trình đạo hàm riêng kinh điển bao gồm: phương trình truyền sóng mô tả dao động của dây đàn hồi, sóng âm, sóng nước hay sóng điện từ; phương trình truyền nhiệt mô tả quy luật biến thiên nhiệt độ theo thời gian và không gian trong một môi trường xác định; phương trình Laplace mô tả trạng thái dừng của các trường thế như thế năng trọng trường, điện thế trong tĩnh điện học hay dòng chảy dừng trong cơ học chất lỏng.

Để giải quyết các phương trình này, toán lý thường sử dụng hai nhánh tiếp cận chính:

  • Phương pháp số: Có khả năng giải quyết các bài toán có điều kiện biên và hình học phức tạp, tuy nhiên kết quả thu được chỉ là nghiệm gần đúng rời rạc.
  • Phương pháp giải tích: Cho phép xác định nghiệm chính xác dưới dạng hàm số giải tích tường minh hoặc chuỗi hàm, phản ánh rõ bản chất vật lý và sự phụ thuộc của nghiệm vào các tham số của bài toán. Mặc dù gặp giới hạn khi cấu trúc biên quá phức tạp, phương pháp giải tích là công cụ nền tảng trong đào tạo đại học chuyên ngành vật lý, nơi các mô hình nghiên cứu thường được chuẩn hóa.

Tại các cơ sở giáo dục đại học, sinh viên ngành Vật lý được trang bị nhiều phương pháp giải tích khác nhau. Các phương pháp phổ biến như phương pháp tách biến (phương pháp Fourier) và phương pháp đa thức d'Alembert thường được giảng dạy rộng rãi nhờ yêu cầu kiến thức toán học ở mức độ vừa phải. Bên cạnh đó, các phương pháp nâng cao hơn như phương pháp biến đổi tích phân (Fourier, Laplace) và phương pháp hàm Green mang lại công cụ tìm nghiệm trực tiếp và tổng quát hơn nhưng đòi hỏi khối lượng kiến thức giải tích toán học phức tạp. Sự đa dạng về dạng phương trình (thuần nhất, không thuần nhất, có nguồn, biên cố định, biên tự do, biên phụ thuộc thời gian) và sự phân tán của các kỹ thuật giải khiến sinh viên gặp khó khăn trong việc xâu chuỗi logic. Đề tài khóa luận được thực hiện nhằm hệ thống hóa một cách đồng bộ các phương pháp giải tích giải phương trình vật lý - toán.

Mục tiêu nghiên cứu:

  1. Đưa ra được hệ thống các phương pháp giải tích để giải phương trình đạo hàm riêng ứng dụng rộng rãi trong vật lý: phương trình truyền sóng, phương trình truyền nhiệt, phương trình Laplace.
  2. Đưa ra hệ thống bài giải các bài tập phương trình đạo hàm riêng đã nêu bằng các phương pháp giải tích: phương pháp tách biến, phương pháp đa thức d'Alembert, phương pháp biến đổi tích phân và phương pháp hàm Green.

Đối tượng nghiên cứu:

  • Các bài toán phương trình đạo hàm riêng ứng dụng trong vật lý.
  • Các phương pháp giải tích áp dụng giải các bài toán đạo hàm riêng trong vật lý.

Nhiệm vụ nghiên cứu:

  1. Tìm hiểu các bài toán đạo hàm riêng thường gặp trong vật lý thông qua các giáo trình, sách và tài liệu chuyên ngành liên quan.
  2. Phân tích các ưu điểm, nhược điểm của các phương pháp giải tích áp dụng giải các bài toán vật lý – toán.

Phạm vi nghiên cứu:

  • Tập trung vào các phương pháp giải tích giải các phương trình vật lý - toán thường gặp: phương trình truyền sóng trên dây hữu hạn, vô hạn và nửa vô hạn; phương trình truyền nhiệt trên thanh hữu hạn (có nguồn và không chứa nguồn); phương trình Laplace hai chiều trong hệ tọa độ Descartes và các hệ tọa độ cong (tọa độ cực, tọa độ cầu).

Cơ sở lý thuyết và phương pháp

Khóa luận thiết lập hệ thống cơ sở toán học giải tích gồm hai nhóm công cụ chính: các hàm đặc biệt và các phép biến đổi tích phân.

1. Khung lý thuyết và các công cụ toán học bổ trợ

  • Hàm delta Dirac ($\delta(x)$): Sử dụng để mô tả các phân bố điểm trong vật lý (mật độ khối lượng điểm, mật độ điện tích điểm). Được định nghĩa bởi $\delta(x) = 0$ khi $x \neq 0$, $\delta(0) = \infty$, thỏa mãn điều kiện chuẩn hóa $\int_{-\infty}^{+\infty} \delta(x)dx = 1$ và tính chất lọc: $$\int_{-\infty}^{+\infty} \delta(x - x_0)f(x)dx = f(x_0)$$
  • Hàm Heaviside ($H(t)$): Hàm bậc thang đơn vị không liên tục, nhận giá trị $H(t) = 0$ khi $t < 0$ và $H(t) = 1$ khi $t > 0$, ứng dụng trong mô hình hóa tín hiệu và mạch điện.
  • Hàm Bessel: Xuất hiện khi giải phương trình vi phân Bessel bậc $m$: $$x^2 \frac{d^2 y}{dx^2} + x \frac{dy}{dx} + (x^2 - m^2)y = 0$$ Nghiệm tổng quát được xác định qua tổ hợp tuyến tính của hàm Bessel loại 1 ($J_m(x)$) và loại 2 ($Y_m(x)$ hoặc $J_{-m}(x)$), áp dụng cho các bài toán vật lý trong miền đối xứng trụ.
  • Đa thức Legendre ($P_n(x)$): Nghiệm của phương trình vi phân Legendre bậc $n$ trên miền $x \in (-1, 1)$: $$(1 - x^2) \frac{d^2 y}{dx^2} - 2x \frac{dy}{dx} + n(n+1)y = 0$$ Biểu diễn qua công thức Rodrigues: $P_n(x) = \frac{1}{2^n n!} \frac{d^n}{dx^n}(x^2 - 1)^n$, đóng vai trò cơ sở trực giao trong các bài toán đối xứng cầu.
  • Các phép biến đổi tích phân:
    • Biến đổi Fourier thực, Fourier Sin ($\mathcal{F}_s$), Fourier Cos ($\mathcal{F}_c$) và các công thức biến đổi ngược tương ứng.
    • Biến đổi Fourier phức: $\mathcal{F}(f(t)) = F(\omega) = \frac{1}{\sqrt{2\pi}} \int_{-\infty}^{+\infty} f(t)e^{-i\omega t}dt$, kèm theo định lý tích chập $\mathcal{F}(f * g) = \mathcal{F}(f)\mathcal{F}(g)$.
    • Biến đổi Laplace: $\mathcal{L}(f(x)) = F(p) = \int_0^{+\infty} f(x)e^{-px}dx$, biến đổi ngược qua tích phân Bromwich $\mathcal{L}^{-1}(F(p)) = \frac{1}{2\pi i} \int_{\gamma - i\infty}^{\gamma + i\infty} F(p)e^{px}dp$ được tính bằng định lý thặng dư Cauchy trên chu tuyến đóng $l + L$ trong mặt phẳng phức: $\int e^{px}F(p)dp = 2\pi i \sum \text{Res}(e^{px}F(p))$.
Hàm số gốc $f(x)$ Ảnh Laplace $F(p)$ Hàm số gốc mở rộng $f(x)$ Ảnh Laplace mở rộng $F(p)$
$a$ $\frac{a}{p}$ $e^{ax}x^n$ $\frac{n!}{(p-a)^{n+1}}$
$e^{ax}$ $\frac{1}{p-a}$ $e^{ax}\sin\omega x$ $\frac{\omega}{(p-a)^2 + \omega^2}$
$x^n$ $\frac{n!}{p^{n+1}}$ $e^{ax}\cos\omega x$ $\frac{p-a}{(p-a)^2 + \omega^2}$
$\sin\omega x$ $\frac{\omega}{p^2 + \omega^2}$ $\sinh\beta x$ $\frac{\beta}{p^2 - \beta^2}$
$\cos\omega x$ $\frac{p}{p^2 + \omega^2}$ $\cosh\beta x$ $\frac{p}{p^2 - \beta^2}$
$\delta(x - T)$ $e^{-pT}$ — —

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

  • Phương pháp phân tích và tổng hợp lý thuyết: Thu thập, đối chiếu và hệ thống hóa các tài liệu toán lý chuyên khảo để xây dựng cấu trúc các phương pháp.
  • Phương pháp mô hình hóa giải tích: Chuyển đổi các bài toán vật lý biên về dạng phương trình vi phân đạo hàm riêng tương ứng, phân loại theo điều kiện biên (Dirichlet, Neumann, Robin hoặc hỗn hợp) và điều kiện ban đầu (Cauchy).
  • Phương pháp so sánh: Đánh giá phạm vi áp dụng, độ phức tạp thuật toán và yêu cầu toán học của từng phương pháp giải.
  • Nguồn dữ liệu: Các giáo trình toán cho vật lý, phương trình đạo hàm riêng và tài liệu chuyên khảo vật lý toán trong nước và quốc tế.

Nội dung chính theo từng chương

Mở đầu

Phần này trình bày bối cảnh học thuật, lý do lựa chọn đề tài xuất phát từ nhu cầu hệ thống hóa kiến thức cho sinh viên sư phạm vật lý; xác định rõ mục đích nghiên cứu, đối tượng, 2 nhiệm vụ nghiên cứu, phạm vi khảo sát và bố cục tổng thể của khóa luận gồm 3 chương chính và phần kết luận.

Chương 1: Cơ sở lý thuyết của đề tài nghiên cứu

Chương 1 tập trung xây dựng hệ thống công cụ toán học hỗ trợ. Tác giả trình bày định nghĩa, tính chất giải tích và ý nghĩa vật lý của các hàm đặc biệt: hàm delta Dirac (mật độ điểm), hàm Heaviside (hàm bước nhảy), hàm Bessel (nghiệm phương trình vi phân trụ), và đa thức Legendre (nghiệm phương trình vi phân cầu). Tiếp theo, chương trình bày chi tiết lý thuyết chuỗi Fourier, tích phân Fourier, biến đổi Fourier Sin/Cos, biến đổi Fourier phức cùng định lý tích chập. Phép biến đổi Laplace được phân tích kèm điều kiện khả tích, công thức tính biến đổi ngược thông qua thặng dư Cauchy tại các cực điểm trên mặt phẳng phức và bảng tra cứu các cặp biến đổi Laplace cơ bản và mở rộng.

Chương 2: Các phương pháp giải tích giải các phương trình vật lý – toán

Chương 2 là trọng tâm lý thuyết của khóa luận, phân tích chi tiết quy trình giải của 4 phương pháp giải tích:

  1. Phương pháp tách biến (Separation of Variables):

    • Nguyên lý: Giả thiết nghiệm có dạng tích các hàm độc lập theo từng biến số (ví dụ $u(x,t) = X(x)T(t)$), đưa phương trình đạo hàm riêng về hệ phương trình vi phân thường thông qua hằng số tách biến $\lambda$. Phân tích chi tiết 3 trường hợp của $\lambda$ ($\lambda < 0, \lambda = 0, \lambda > 0$) kết hợp điều kiện biên để loại bỏ nghiệm tầm thường và xác định các giá trị riêng cùng hàm riêng.
    • Phương trình truyền sóng một chiều: $\frac{\partial^2 u}{\partial x^2} - \frac{1}{c^2}\frac{\partial^2 u}{\partial t^2} = 0$. Khảo sát chi tiết 4 dạng điều kiện biên trên đoạn $[0, l]$: hai đầu cố định ($u(0,t)=u(l,t)=0$), hai đầu tự do ($u'(0,t)=u'(l,t)=0$), đầu $x=0$ cố định và $x=l$ tự do, đầu $x=0$ tự do và $x=l$ cố định. Nghiệm tổng quát được thiết lập dưới dạng chuỗi Fourier lượng giác với các hệ số $A_k, B_k$ tính từ điều kiện đầu $u(x,0)=f(x)$ và $\frac{\partial u}{\partial t}(x,0)=F(x)$.
    • Bài toán truyền sóng cưỡng bức: Phương trình $\frac{\partial^2 u}{\partial t^2} = a^2 \frac{\partial^2 u}{\partial x^2} + g(x,t)$. Giải quyết trường hợp hai đầu buộc chặt bằng cách khai triển hàm nguồn $g(x,t) = \sum G_k(t)\sin\frac{k\pi x}{l}$ để tìm nghiệm riêng của phương trình vi phân thường cho $T_k(t)$. Đối với trường hợp hai đầu dây chuyển động theo quy luật cho trước $u(0,t)=\varphi_1(t)$ và $u(l,t)=\varphi_2(t)$, tác giả áp dụng phép đổi biến hàm số $u(x,t) = v(x,t) + \varphi_1(t) + \frac{x}{l}[\varphi_2(t) - \varphi_1(t)]$ để chuyển bài toán về dạng có điều kiện biên thuần nhất cho hàm $v(x,t)$.
    • Phương trình truyền nhiệt một chiều: $\frac{\partial u}{\partial t} = c^2 \frac{\partial^2 u}{\partial x^2}$. Khảo sát thanh hữu hạn không chứa nguồn với 4 cấu hình biên (hai đầu giữ $0^\circ\text{C}$, hai đầu cách nhiệt, một đầu cách nhiệt một đầu $0^\circ\text{C}$). Khảo sát thanh có chứa nguồn $g(x,t)$ và bài toán có điều kiện biên phụ thuộc thời gian bằng phương pháp đổi biến khử biên.
    • Phương trình Laplace hai chiều: $\nabla^2 u = \frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = 0$ trên miền chữ nhật $[0, l] \times [0, m]$ với các điều kiện biên Dirichlet thuần nhất và không thuần nhất.
    • Phương pháp tách biến trong các hệ tọa độ khác: Giải phương trình Laplace trong hệ tọa độ cực trên miền tròn bán kính $r_0$ (kết hợp điều kiện biên tuần hoàn $\Phi(\varphi) = \Phi(\varphi + 2\pi)$ và điều kiện nghiệm hữu hạn tại tâm $r=0$ để loại bỏ số hạng logarit $C_1 \ln r$); giải phương trình Laplace trong hệ tọa độ cầu có đối xứng trục trên miền cầu bán kính $a$, đưa về phương trình vi phân Legendre và nghiệm chuỗi đa thức Legendre $P_n(\cos\theta)$.

  2. Phương pháp đa thức d'Alembert:

    • Trình bày kỹ thuật giải phương trình truyền sóng trên dây dài vô hạn $(-\infty, +\infty)$ bằng cách chuyển đổi hệ tọa độ đặc trưng $(\xi = x - at, \eta = x + at)$, đưa ra công thức nghiệm d'Alembert: $$u(x,t) = \frac{f(x-at) + f(x+at)}{2} + \frac{1}{2a} \int_{x-at}^{x+at} F(s)ds$$
    • Khảo sát bài toán sóng trên dây nửa vô hạn $[0, +\infty)$ bằng phương pháp phản xạ sóng (kéo dài hàm điều kiện đầu chẵn hoặc lẻ tùy thuộc điều kiện biên cố định hay tự do tại gốc $x=0$).

  3. Phương pháp biến đổi tích phân:

    • Ứng dụng biến đổi Fourier và Laplace để triệt tiêu các đạo hàm riêng theo biến không gian hoặc thời gian, chuyển phương trình đạo hàm riêng thành phương trình đại số hoặc phương trình vi phân thường trên miền ảnh. Sau đó giải tìm nghiệm ảnh và thực hiện phép biến đổi tích phân ngược để thu về nghiệm thực.

  4. Phương pháp hàm Green:

    • Trình bày định nghĩa hàm Green như đáp ứng của hệ thống đối với nguồn điểm dạng hàm delta Dirac. Xây dựng nghiệm hàm Green cho phương trình sóng độc lập với thời gian (phương trình Helmholtz), phương trình sóng không thuần nhất trong không gian 3 chiều và ứng dụng cho hệ phương trình Maxwell phụ thuộc thời gian (thế trễ điện từ).

Chương 3: Áp dụng các phương pháp giải tích trong việc giải các phương trình vật lý - toán

Chương 3 trình bày hệ thống các bài tập cụ thể minh họa cho từng phương pháp đã xây dựng ở Chương 2:

  • Hệ thống bài tập áp dụng phương pháp tách biến: Giải chi tiết bài toán truyền sóng trên dây hữu hạn với các điều kiện đầu xác định; bài toán phân bố nhiệt độ trên thanh đồng tính khi biết phân bố nhiệt ban đầu; bài toán phân bố điện thế trong miền chữ nhật, miền tròn và miền cầu.
  • Hệ thống bài tập áp dụng phương pháp đa thức d'Alembert: Xác định dạng sóng và đồ thị lan truyền sóng $y = u(x,t)$ tại các thời điểm khác nhau.
  • Hệ thống bài tập áp dụng phương pháp biến đổi tích phân (Fourier, Laplace) cho các bài toán truyền sóng và truyền nhiệt trên miền vô hạn hoặc nửa vô hạn.
  • Hệ thống bài tập áp dụng phương pháp hàm Green: Tính toán điện thế và trường sóng tạo bởi các nguồn điểm và phân bố nguồn liên tục.

Kết quả và đóng góp

Khóa luận đã tổng hợp và chuẩn hóa quy trình giải toán cho các bài toán biên vật lý. Các kết quả giải tích chính được so sánh tổng hợp trong bảng sau:

Phương pháp giải tích Loại phương trình và miền khảo sát tối ưu Kỹ thuật xử lý toán học cốt lõi Ưu điểm chính Hạn chế kỹ thuật
Phương pháp tách biến Phương trình sóng, nhiệt, Laplace trên miền hữu hạn có biên đối xứng (Descartes, cực, cầu). Tách nghiệm $u = X \cdot Y$ hoặc $X \cdot T$, giải bài toán Sturm-Liouville tìm giá trị riêng và hàm riêng. Nghiệm dạng chuỗi tường minh, dễ phân tích chế độ dao động riêng/mode nhiệt. Chỉ áp dụng được khi miền biên có thể tách biến tọa độ; hội tụ chậm ở một số điều kiện biên gián đoạn.
Phương pháp đa thức d'Alembert Phương trình truyền sóng 1D trên miền vô hạn $(-\infty, +\infty)$ hoặc nửa vô hạn $[0, +\infty)$. Sử dụng tọa độ đặc trưng $\xi = x-at, \eta = x+at$ và phương pháp thác triển điều kiện đầu. Công thức nghiệm dạng tích phân đóng trực tiếp, thể hiện rõ bản chất sóng chạy. Bị giới hạn chủ yếu cho phương trình truyền sóng một chiều, khó áp dụng cho không gian nhiều chiều có biên phức tạp.
Phương pháp biến đổi tích phân Phương trình sóng, nhiệt trên miền vô hạn hoặc nửa vô hạn ($t \ge 0$, $x \in [0, \infty)$). Sử dụng toán tử tích phân Fourier (Sin/Cos/phức) hoặc Laplace, kết hợp định lý tích chập và thặng dư Cauchy. Chuyển PDE trực tiếp về ODE hoặc phương trình đại số, giải quyết tốt điều kiện đầu không tuần hoàn. Việc tính tích phân ngược đòi hỏi kỹ thuật giải tích phức nâng cao (tính thặng dư trên chu tuyến phức).
Phương pháp hàm Green Phương trình vi phân không thuần nhất, phương trình Helmholtz, bài toán điện động lực học (Maxwell). Tìm nghiệm cơ bản ứng với nguồn xung $\delta(\vec{r}-\vec{r}')$, sau đó lấy tích phân chập với hàm nguồn $g(\vec{r})$. Cung cấp nghiệm tổng quát cho mọi dạng hàm nguồn; bản chất vật lý tường minh (nguyên lý chồng chất). Việc xây dựng hàm Green thỏa mãn đồng thời các điều kiện biên là bài toán toán học rất khó.

Các kết quả đạt được:

  • Xây dựng hoàn chỉnh thuật toán giải cho 4 cấu hình điều kiện biên của phương trình truyền sóng và 4 cấu hình biên của phương trình truyền nhiệt trên miền 1D.
  • Trình bày giải pháp kỹ thuật khử điều kiện biên không thuần nhất phụ thuộc thời gian ($\varphi_1(t), \varphi_2(t)$) thông qua phép tịnh tiến nghiệm tuyến tính theo biến không gian.
  • Thiết lập công thức nghiệm tường minh cho phương trình Laplace trong hệ tọa độ cực (chuỗi Fourier theo góc $\varphi$) và hệ tọa độ cầu đối xứng trục (chuỗi đa thức Legendre $P_n(\cos\theta)$).

Đóng góp của đề tài:

  • Cung cấp một tài liệu học tập có tính hệ thống cao, kết nối chặt chẽ giữa công cụ giải tích toán học trừu tượng và các bài toán mô hình hóa vật lý thực tế.
  • Phân loại rõ ràng phạm vi ứng dụng và cấu trúc giải của từng phương pháp, giúp người học định hướng chính xác phương pháp tiếp cận khi đối mặt với từng dạng bài toán biên cụ thể.

Hạn chế và hướng nghiên cứu tiếp

  • Hạn chế:
    • Khóa luận tập trung hoàn toàn vào các bài toán tuyến tính với các miền hình học chuẩn (đoạn thẳng 1D, miền chữ nhật 2D, miền tròn và miền cầu đối xứng trục); chưa mở rộng sang các hệ phương trình phi tuyến hoặc các bài toán trong môi trường dị hướng, không đồng nhất.
    • Các ví dụ minh họa chủ yếu tập trung vào không gian một chiều và hai chiều đối xứng; chưa đi sâu vào các bài toán không gian ba chiều tổng quát không có tính đối xứng.
  • Hướng nghiên cứu tiếp:
    • Tiếp tục mở rộng nghiên cứu phương pháp hàm Green và biến đổi tích phân cho các bài toán phương trình sóng và nhiệt trong không gian 3 chiều không đối xứng.
    • Khảo sát sự kết hợp giữa phương pháp giải tích và các thuật toán tính toán số hoặc sử dụng các hệ thống đại số máy tính (Computer Algebra Systems) để giải các bài toán biên có dạng hình học phức tạp hơn.

Giá trị tham khảo

  • Đối tượng sử dụng:
    • Sinh viên ngành Sư phạm Vật lý, Vật lý học, Kỹ thuật và Toán ứng dụng cần tài liệu học tập, ôn tập cho học phần Phương trình Vật lý - Toán hoặc Phương pháp Toán cho Vật lý.
    • Giảng viên phụ trách các học phần liên quan có thể sử dụng cấu trúc phân loại bài toán và hệ thống bài giải mẫu làm tài liệu tham khảo giảng dạy.
  • Các nội dung có giá trị tham khảo nổi bật:
    • Quy trình chi tiết biện luận hằng số tách biến $\lambda$ và cách loại bỏ nghiệm tầm thường trong phương pháp tách biến.
    • Kỹ thuật chuyển đổi bài toán biên không thuần nhất thành thuần nhất trong bài toán dao động cưỡng bức và truyền nhiệt có nguồn.
    • Bảng tra cứu biến đổi Laplace và quy trình tính tích phân Laplace ngược bằng định lý thặng dư Cauchy.

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

1. Trong phương pháp tách biến cho phương trình truyền sóng hoặc truyền nhiệt, vì sao trường hợp hằng số tách biến $\lambda < 0$ và $\lambda = 0$ luôn dẫn đến nghiệm tầm thường?
Khi $\lambda < 0$ (đặt $\lambda = -\alpha^2$), nghiệm của phương trình không gian $X(x)$ có dạng hàm mũ $X(x) = Ae^{-\alpha x} + Be^{\alpha x}$. Khi áp dụng các điều kiện biên thuần nhất (như hai đầu cố định $X(0)=X(l)=0$), ta thu được hệ phương trình đại số tuyến tính có định thức khác không, dẫn đến $A = B = 0$. Tương tự, khi $\lambda = 0$, nghiệm có dạng bậc nhất $X(x) = Ax + B$, điều kiện biên $X(0)=X(l)=0$ cũng buộc $A = B = 0$. Cả hai trường hợp đều dẫn đến $X(x) \equiv 0$, tức $u(x,t) \equiv 0$ (nghiệm tầm thường), do đó chỉ có $\lambda > 0$ mới sinh ra các nghiệm dao động riêng phi tầm thường.

2. Khóa luận giải quyết bài toán biên không thuần nhất phụ thuộc thời gian $u(0,t)=\varphi_1(t)$ và $u(l,t)=\varphi_2(t)$ như thế nào?
Tác giả sử dụng phép đổi biến hàm số: $$u(x,t) = v(x,t) + \varphi_1(t) + \frac{x}{l}[\varphi_2(t) - \varphi_1(t)]$$ Phép biến đổi tuyến tính theo biến $x$ này giúp triệt tiêu các điều kiện biên không thuần nhất tại hai đầu $x=0$ và $x=l$, đưa bài toán về tìm hàm $v(x,t)$ thỏa mãn điều kiện biên thuần nhất $v(0,t) = v(l,t) = 0$. Khi đó, bài toán cho $v(x,t)$ trở thành bài toán truyền sóng hoặc truyền nhiệt có nguồn điều chỉnh $g_1(x,t)$, có thể giải tiếp bằng phương pháp tách biến hoặc khai triển Fourier.

3. Phương trình Laplace trong hệ tọa độ cầu đối xứng trục dẫn đến phương trình vi phân nào và có nghiệm tổng quát ra sao?
Khi thế $u(r,\theta) = R(r)\Phi(\theta)$ vào phương trình Laplace đối xứng trục và đặt $\omega = \cos\theta$, phương trình cho phần góc trở thành phương trình vi phân Legendre: $$(1-\omega^2)\frac{d^2\Phi}{d\omega^2} - 2\omega\frac{d\Phi}{d\omega} + n(n+1)\Phi = 0$$ Nghiệm của phương trình này là các đa thức Legendre $P_n(\cos\theta)$. Kết hợp với nghiệm của phương trình phần bán kính $R_n(r) = A_n r^n + B_n r^{-(n+1)}$ (với $B_n = 0$ để đảm bảo nghiệm hữu hạn tại tâm cầu), nghiệm tổng quát trong miền cầu bán kính $a$ có dạng: $$u(r,\theta) = \sum_{n=0}^{+\infty} A_n r^n P_n(\cos\theta)$$

4. Kỹ thuật biến đổi Laplace ngược được thực hiện như thế nào trong khóa luận?
Biến đổi Laplace ngược được xác định qua tích phân Bromwich trên mặt phẳng phức: $$\mathcal{L}^{-1}(F(p)) = f(x) = \frac{1}{2\pi i} \int_{\gamma - i\infty}^{\gamma + i\infty} F(p)e^{px}dp$$ Tác giả áp dụng định lý thặng dư Cauchy bằng cách dựng chu tuyến đóng $l + L$ bao quanh toàn bộ các điểm bất thường cô lập (cực điểm $a_j$) của hàm dưới dấu tích phân. Tích phân được tính thông qua tổng thặng dư: $$f(x) = \sum_{j=1}^k \text{Res}_{z=a_j} \left[ e^{zx} F(z) \right]$$

Kết luận

Khóa luận tốt nghiệp của tác giả Huỳnh Trúc Phương đã hệ thống hóa một cách toàn diện và chặt chẽ 4 phương pháp giải tích chính (tách biến, đa thức d'Alembert, biến đổi tích phân, hàm Green) dùng để giải các bài toán phương trình đạo hàm riêng kinh điển trong vật lý. Đề tài đã xây dựng thành công quy trình giải chuẩn tắc cho các bài toán truyền sóng, truyền nhiệt và phương trình Laplace trên các miền không gian một chiều, hai chiều Descartes, hệ tọa độ cực và tọa độ cầu. Công trình là tài liệu học tập và tham khảo có giá trị sư phạm cao, phục vụ hiệu quả cho việc giảng dạy và học tập học phần Phương pháp toán cho vật lý tại các trường đại học sư phạm và đào tạo khoa học cơ bản.