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

Trong chương trình đào tạo đại học các ngành Sư phạm Vật lý và Vật lý học, học phần Phương pháp toán - lý giữ vai trò cầu nối quan trọng giữa các công cụ toán học cao cấp với các mô hình vật lý lý thuyết. Một trong những nội dung trọng tâm của học phần này là giải các bài toán về phương trình đạo hàm riêng mô tả hiện tượng vật lý, điển hình là quá trình truyền dẫn nhiệt trong vật thể.

Theo tác giả khóa luận, đối với các bài toán truyền nhiệt nhiều chiều, các phương pháp giải giải tích truyền thống như phương pháp tách biến Fourier hay phương pháp biến đổi Laplace thường bộc lộ nhiều hạn chế, dẫn đến việc thiết lập và tìm nghiệm giải tích gặp nhiều trở ngại và biến đổi phức tạp. Phương pháp hàm Bessel là một công cụ giải tích hữu hiệu để giải các bài toán biên nhiều chiều, đặc biệt đối với các hệ vật lý có tính đối xứng trục hoặc cấu trúc hình trụ. Tuy nhiên, các giáo trình lý thuyết hiện hành thường ít đi sâu vào việc hướng dẫn phương pháp hàm Bessel, đồng thời thiếu hệ thống bài tập mẫu và bài tập áp dụng cụ thể, gây khó khăn cho sinh viên trong quá trình tiếp cận và vận dụng. Vì vậy, việc nghiên cứu có hệ thống và xây dựng tài liệu hoàn chỉnh về phương pháp này cho học phần Phương pháp toán - lý là yêu cầu thực tiễn cần thiết.

Mục tiêu nghiên cứu của khóa luận được xác định gồm ba nội dung chính:

  1. Tìm hiểu các bài toán truyền nhiệt.
  2. Nghiên cứu cơ sở toán học cho hàm Bessel.
  3. Sử dụng hàm Bessel để tìm nghiệm của bài toán phương trình truyền nhiệt.

Nhiệm vụ nghiên cứu được phân định thành hai phần:

  1. Xây dựng phương pháp hàm Bessel để tìm nghiệm của bài toán truyền nhiệt.
  2. Giải một số bài toán truyền nhiệt bằng phương pháp hàm Bessel.

Đối tượng nghiên cứu của đề tài bao gồm ba thành tố:

  1. Các bài toán truyền nhiệt.
  2. Cơ sở lý luận về bài tập vật lý trong dạy học bộ môn.
  3. Cơ sở toán học về phương trình Bessel và hàm Bessel.

Phạm vi nghiên cứu của khóa luận giới hạn trong các bài toán truyền nhiệt giải tích trong vật dẫn. Giả thuyết khoa học được đặt ra là: việc sử dụng phương pháp hàm Bessel có thể tìm được nghiệm giải tích của các bài toán truyền nhiệt tương ứng.

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

Khóa luận thuộc loại hình nghiên cứu lý thuyết giải tích và phương pháp toán lý. Để giải quyết bài toán đặt ra, tác giả kết hợp chặt chẽ giữa các quy luật nhiệt động lực học và các công cụ toán học giải tích hiện đại.

Về cơ sở vật lý, tác giả xuất phát từ các định luật cơ bản:

  • Định luật Fourier về quá trình truyền nhiệt: Vectơ dòng nhiệt được xác định qua hệ thức $\vec{q} = -k \nabla u$, trong đó $k$ là hệ số dẫn nhiệt và dấu trừ chỉ chiều dòng nhiệt truyền từ nơi có nhiệt độ cao sang nơi có nhiệt độ thấp.
  • Định lý Gauss: Chuyển đổi tích phân mặt kín của thông lượng nhiệt sang tích phân thể tích $\iiint_V \text{div}\vec{q}, d\tau = \iint_S \vec{q}\cdot\vec{n}, d\sigma$.
  • Định luật bảo toàn năng lượng: Thiết lập cân bằng nhiệt lượng biến thiên trong thể tích kiểm soát với lượng nhiệt đi qua mặt biên và nguồn nhiệt nội sinh, dẫn đến phương trình vi phân truyền nhiệt tổng quát $\text{div}(k \text{grad} u) + H = \frac{\partial}{\partial t}(c\rho u)$.
  • Định luật Fick về khuếch tán: Mở rộng mô hình truyền nhiệt sang quá trình khuếch tán vật chất với phương trình dòng khuếch tán $\vec{J} = -D \nabla C$.

Về cơ sở toán học, đề tài xây dựng trên nền tảng:

  • Lý thuyết phương trình vi phân đạo hàm riêng tuyến tính và phương pháp tách biến (Fourier separation of variables).
  • Lý thuyết bài toán biên và các điều kiện biên kinh điển.
  • Chuỗi Fourier, định lý Dirichlet về sự hội tụ, chuỗi lũy thừa và chuỗi Taylor/Maclaurin.
  • Tích phân suy rộng loại 1 (cận vô tận) và loại 2 (hàm không bị chặn), hàm Gamma Euler $\Gamma(k)$.
  • Lý thuyết toán tử tuyến tính, hàm riêng, trị riêng và nguyên lý chồng chất trạng thái.
  • Hệ phương trình vi phân Bessel, phương trình vi phân Legendre và phương trình hàm cầu trong hệ tọa độ cầu.
Loại điều kiện biên Bản chất toán lý Biểu thức điều kiện biên
Điều kiện loại I (Dirichlet) Nhiệt độ được cố định trên biên bề mặt $S$ $u(x,y,z,t)\vert_S = f_1(x,y,z,t)$
Điều kiện loại II (Neumann) Dòng nhiệt qua biên được cho trước; biên cách nhiệt khi $f_2 = 0$ $\frac{\partial u}{\partial n}\vert_S = \vec{n}\cdot\text{grad},u\vert_S = f_2(x,y,z,t)$
Điều kiện loại III (Robin) Dòng nhiệt qua biên tỷ lệ với trao đổi nhiệt môi trường ngoài $\left[\frac{\partial u}{\partial n} + h u\right]_S = f_3(x,y,z,t)\ (h > 0)$
Điều kiện hỗn hợp Kết hợp các điều kiện biên loại I và II trên từng phần bề mặt Kết hợp các biểu thức trên biên

Phương pháp nghiên cứu thực tế được tác giả sử dụng bao gồm:

  • Phương pháp đọc sách và tham khảo tài liệu: Khai thác các giáo trình phương pháp toán - lý, sách giải tích chuyên sâu và tài liệu vật lý lý thuyết liên quan.
  • Phương pháp toán học: Vận dụng các thuật toán giải tích, vi phân, tích phân và khai triển chuỗi để thiết lập công thức và tìm nghiệm phương trình.
  • Phương pháp phân tích: Phân tách hiện tượng và phương trình nhiều chiều thành các bài toán một chiều đơn giản hơn theo từng biến không gian và thời gian.
  • Phương pháp đàm thoại trao đổi ý kiến với giáo viên: Trao đổi học thuật định kỳ để định hướng cấu trúc nội dung và kiểm tra tính chính xác của các dẫn giải giải tích.

Nguồn dữ liệu của khóa luận là các định lý, công thức giải tích chuẩn mực và hệ thống mô hình bài toán biên được trích xuất từ các tài liệu tham khảo chuyên ngành phương pháp toán lý và cơ học nhiệt.

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

CHƯƠNG I: CƠ SỞ LÝ LUẬN CỦA ĐỀ TÀI

Chương I gồm 4 nội dung lý luận và công cụ toán học:

  1. Lý luận về bài tập vật lý: Xác định bản chất bài tập vật lý là phương tiện nghiên cứu hiện tượng, hình thành khái niệm, phát triển tư duy logic, rèn luyện kỹ năng vận dụng vào thực tiễn, giáo dục ý chí vượt khó và kiểm tra mức độ nắm vững kiến thức. Tác giả phân tích hai thao tác tư duy chủ đạo là phương pháp phân tích (đi từ cái cần tìm đến cái đã biết) và phương pháp tổng hợp (đi từ cái đã biết đến cái cần tìm), đồng thời khái quát hóa tiến trình giải bài toán qua 4 bước: (1) Tìm hiểu đề bài; (2) Xác lập mối liên hệ cơ bản; (3) Rút ra kết quả cần tìm; (4) Kiểm tra xác nhận kết quả (tính toán, thứ nguyên, so sánh cách giải khác).
  2. Bài toán biên: Khảo sát cấu trúc nghiệm của phương trình vi phân tuyến tính cấp $n$, nghiệm tổng quát $y = y_c + y_p$. Trình bày bài toán giá trị ban đầu và phân tích sự phụ thuộc của nghiệm vào các dạng điều kiện biên hai điểm và điều kiện biên hỗn hợp. Giới thiệu phương pháp biến thiên hằng số để giải phương trình vi phân cấp 1.
  3. Các chuỗi và hệ trực giao: Trình bày định nghĩa chuỗi lượng giác, chuỗi Fourier, định lý Dirichlet về điều kiện đủ để khai triển hàm, chuỗi hàm chẵn/lẻ và hàm chu kỳ $2L$. Trình bày lý thuyết chuỗi lũy thừa, tính chất giải tích của tổng chuỗi, chuỗi Taylor và Maclaurin. Định nghĩa hệ hàm trực giao trên đoạn $[a, b]$, khảo sát các hệ hàm lượng giác trực giao cơ bản. Giới thiệu hàm Gamma Euler $\Gamma(k)$ với tính chất $\Gamma(k+1) = k\Gamma(k)$ và lý thuyết tích phân suy rộng (cận vô tận, hàm không bị chặn cùng các dấu hiệu so sánh, tiêu chuẩn Dirichlet).
  4. Khái niệm toán tử, hàm riêng, trị riêng: Định nghĩa toán tử tuyến tính gắn liền với nguyên lý chồng chất trạng thái; nêu ví dụ về toán tử tọa độ, toán tử vi phân và toán tử Laplace. Xây dựng phương trình trị riêng $\hat{A}\Psi = A\Psi$ và giải bài toán trị riêng mẫu cho toán tử vi phân $\hat{A} = -i\frac{\partial}{\partial x}$ với điều kiện biên tuần hoàn trên khoảng $(0, L)$, xác định phổ trị riêng gián đoạn $A_n = \frac{2n\pi}{L}$.

CHƯƠNG II: XÂY DỰNG PHƯƠNG TRÌNH BESSEL VÀ HÀM BESSEL

Chương II tập trung hoàn chỉnh cơ sở giải tích về phương trình Bessel và các hàm đặc biệt liên quan:

  1. Khái niệm hàm Bessel: Xuất phát từ phương trình Laplace $\Delta u = 0$ khi xét bài toán dừng của phương trình sóng và phương trình truyền nhiệt. Chuyển toán tử Laplace sang tọa độ trụ $(r, \phi, z)$, áp dụng phương pháp tách biến $u(r, \phi, z) = V(r, z)\Phi(\phi)$ và $V(r, z) = R(r)Z(z)$, từ đó thu được phương trình vi phân Bessel trong tọa độ trụ đối với biến bán kính: $\frac{d^2 R}{dx^2} + \frac{1}{x}\frac{dR}{dx} + \left(1 - \frac{\nu^2}{x^2}\right)R = 0$ với $x = \lambda r$.
  2. Cơ sở cho việc xây dựng hàm Bessel, phương trình hàm Bessel: Tìm nghiệm của phương trình Bessel quanh điểm kỳ dị $x = 0$ dưới dạng chuỗi Frobenius $y = x^\rho \sum_{k=0}^\infty a_k x^k$. Thiết lập phương trình chỉ số $\rho = \pm\nu$. Biểu diễn hệ số qua hàm Gamma và thu được hàm Bessel loại I 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}$$ Khảo sát tính phụ thuộc tuyến tính khi $\nu$ là số nguyên $n$: $J_{-n}(x) = (-1)^n J_n(x)$. Xây dựng nghiệm độc lập thứ hai là hàm Bessel loại II cấp $\nu$ ($Y_\nu(x)$), sử dụng quy tắc L'Hôpital khi $\nu \to n$ để xác định tường minh $Y_n(x)$ và trường hợp riêng $Y_0(x)$. Ngoài ra, chương này mở rộng khảo sát phương trình Legendre, xây dựng đa thức Legendre $P_n(x)$ bằng công thức Rodrigues, chứng minh tính trực giao và giải phương trình hàm cầu $Y_n(\theta, \phi)$ trong tọa độ cầu.
  3. Tính trực giao của hàm Bessel: Thiết lập công thức trực giao thứ nhất đối với các nghiệm dương $\mu_i, \mu_j$ của phương trình $J_\nu(x) = 0$: $$\int_0^L x J_\nu\left(\mu_i \frac{x}{L}\right) J_\nu\left(\mu_j \frac{x}{L}\right) dx = \begin{cases} 0 & \text{khi } i \neq j \ \frac{L^2}{2} J_{\nu+1}^2(\mu_i) & \text{khi } i = j \end{cases}$$ Thiết lập công thức trực giao thứ hai khi $\mu_i, \mu_j$ là nghiệm của phương trình $\alpha J_\nu(x) + \beta x J'_\nu(x) = 0$.
  4. Các hệ thức liên quan đến hàm Bessel: Xây dựng phương pháp khai triển hàm tùy ý thành chuỗi Fourier - Bessel và chuỗi Dyni - Bessel, đưa ra công thức giải tích xác định hệ số khai triển $a_i$.

CHƯƠNG III: SỬ DỤNG HÀM BESSEL ĐỂ GIẢI MỘT SỐ BÀI TOÁN TRUYỀN NHIỆT

Chương III áp dụng trực tiếp hàm Bessel để giải quyết các bài toán vật lý cụ thể:

  1. Thiết lập phương trình truyền nhiệt: Dẫn giải phương trình truyền nhiệt từ định luật bảo toàn năng lượng và định luật Fourier thông qua định lý Gauss. Xét các trường hợp đặc biệt cho vật liệu đồng chất, phương trình trạng thái dừng và tổng hợp các dạng phương trình trong hệ tọa độ Descartes. Giới thiệu tương đương hình thức giữa phương trình truyền nhiệt và phương trình khuếch tán dựa trên định luật Fick.
  2. Các bài toán cho các tọa độ:
    • Tọa độ trụ xuyên tâm: Giải bài toán truyền nhiệt trong thanh trụ dài có tiết diện tròn bán kính $r_0$, nhiệt độ chỉ phụ thuộc bán kính $r$ và thời gian $t$: $\frac{\partial u}{\partial t} = a^2 \left(\frac{\partial^2 u}{\partial r^2} + \frac{1}{r}\frac{\partial u}{\partial r}\right)$ với điều kiện biên loại I trên bề mặt $u(r_0, t) = 0$ và điều kiện ban đầu $u(r, 0) = f(r)$. Bằng phương pháp tách biến $u(r, t) = R(r)T(t)$, bài toán theo biến $r$ quy về phương trình Bessel cấp không. Do nghiệm phải hữu hạn tại tâm $r = 0$, loại bỏ thành phần $Y_0(\omega r)$, thu được hàm riêng $R_i(r) = J_0(\omega_i r)$ với $\omega_i = \mu_i^{(0)}/r_0$, trong đó $\mu_i^{(0)}$ là các không điểm của $J_0$. Nghiệm tổng quát có dạng chuỗi: $$u(r, t) = \sum_{i=1}^\infty A_i J_0\left(\frac{\mu_i^{(0)}}{r_0} r\right) e^{-\left(\frac{\mu_i^{(0)} a}{r_0}\right)^2 t}$$ Nhờ tính trực giao của hàm Bessel, hệ số $A_i$ được xác định tường minh: $$A_i = \frac{2}{r_0^2 J_1^2(\mu_i^{(0)})} \int_0^{r_0} r f(r) J_0\left(\frac{\mu_i^{(0)}}{r_0} r\right) dr$$
    • Tọa độ trụ không xuyên tâm: Xét ống trụ tròn dài vô hạn bán kính $r_0$, nhiệt độ ban đầu phụ thuộc cả bán kính và góc cực $u(r, \phi, 0) = f(r, \phi)$, bề mặt duy trì nhiệt độ bằng không $u(r_0, \phi, t) = 0$. Sử dụng nghiệm tách biến ba thành phần $u(r, \phi, t) = R(r)\Phi(\phi)T(t)$. Tính tuần hoàn theo góc cực $\Phi(\phi + 2\pi) = \Phi(\phi)$ dẫn đến bậc $n$ nguyên. Phương trình bán kính trở thành phương trình Bessel cấp $n$. Nghiệm tổng quát được thiết lập dưới dạng chuỗi hai chỉ số Bessel - lượng giác: $$u(r, \phi, t) = \sum_{n=0}^\infty \sum_{k=0}^\infty J_n\left(\frac{\mu_k^{(n)}}{r_0} r\right) e^{-\left(\frac{\mu_k^{(n)} a}{r_0}\right)^2 t} (A_{nk} \cos n\phi + B_{nk} \sin n\phi)$$ Các hệ số $A_{nk}, B_{nk}$ được tính toán chính xác dựa trên tích phân trực giao kết hợp giữa hàm Bessel và các hàm lượng giác.
  3. Các bài toán cho biên và Một số bài toán dừng: Phân loại và khảo sát các phương án giải phương trình truyền nhiệt khi điều kiện biên chuyển sang loại II, loại III hoặc bài toán ở trạng thái ổn định (dừng) theo thời gian.

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

Khóa luận đạt được các kết quả giải tích cụ thể:

  • Thiết lập hệ thống bảng dạng toán tử và dạng một chiều của phương trình truyền nhiệt trong hệ tọa độ Descartes:
Trường hợp Dạng toán tử Dạng một chiều
Tổng quát $\nabla(k\nabla u) + H = \rho c \frac{\partial u}{\partial t}$ $\frac{\partial}{\partial x}\left(k\frac{\partial u}{\partial x}\right) + H = \rho c \frac{\partial u}{\partial t}$
Vật liệu đồng chất $\nabla^2 u + \frac{H}{k} = \frac{\rho c}{k} \frac{\partial u}{\partial t}$ $\frac{\partial^2 u}{\partial x^2} + \frac{H}{k} = \frac{\rho c}{k} \frac{\partial u}{\partial t}$
Trạng thái dừng $\nabla(k\nabla u) + H = 0$ $\frac{\partial}{\partial x}\left(k\frac{\partial u}{\partial x}\right) + H = 0$
Trạng thái dừng với vật liệu đồng chất $\nabla^2 u + \frac{H}{k} = 0$ $\frac{d^2 u}{dx^2} + \frac{H}{k} = 0$
  • Xây dựng nghiệm giải tích tường minh cho hai mô hình truyền nhiệt trong hệ tọa độ trụ:
    1. Nghiệm bài toán truyền nhiệt xuyên tâm đối xứng trục với hệ số $A_i$ được tính bằng tích phân trực giao chứa hàm $J_1^2(\mu_i^{(0)})$.
    2. Nghiệm bài toán truyền nhiệt không xuyên tâm phụ thuộc góc cực $\phi$ kết hợp giữa chuỗi hàm Bessel $J_n$ và hệ hàm trực giao lượng giác $\cos n\phi, \sin n\phi$.

Đóng góp của đề tài do chính tác giả xác định:

  • Cung cấp tài liệu tham khảo có cấu trúc hoàn chỉnh, gắn kết từ lý thuyết chuỗi, toán tử đến các bài toán vật lý thực tế.
  • Đưa ra phương pháp giải hiệu quả bài toán truyền nhiệt nhiều chiều thay thế cho các biến đổi phức tạp của biến đổi Laplace hay phương pháp Fourier thông thường.
  • Góp phần nâng cao kết quả học tập và rèn luyện kỹ năng giải bài tập phương pháp toán lý cho sinh viên ngành Vật lý.

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

Khóa luận phản ánh một số điểm giới hạn dựa trên phạm vi nghiên cứu:

  • Đề tài tập trung giải tích vào vật dẫn có dạng hình học đối xứng trụ với các điều kiện biên thuần nhất; chưa đi sâu giải các bài toán có nguồn nhiệt phụ thuộc phi tuyến vào nhiệt độ hoặc bài toán trong môi trường vật liệu dị hướng phức tạp.
  • Khóa luận dừng lại ở các bước suy diễn giải tích và thiết lập công thức nghiệm giải tích dạng chuỗi; chưa xây dựng các thuật toán tính toán số để mô phỏng hình ảnh phân bố trường nhiệt độ theo thời gian.

Hướng nghiên cứu tiếp theo được tác giả gợi mở:

  • Mở rộng ứng dụng hàm Bessel để tìm nghiệm giải tích cho các bài toán biên nhiều chiều có cấu trúc phức tạp hơn.
  • Khảo sát sâu hơn các bài toán dừng và bài toán không dừng với các điều kiện biên loại II (Neumann), loại III (Robin) và điều kiện biên hỗn hợp trong các hệ tọa độ cong khác.

Giá trị tham khảo

Văn bản khóa luận mang lại giá trị học thuật và thực hành giảng dạy đối với:

  • Sinh viên chuyên ngành Sư phạm Vật lý, Vật lý học và Toán ứng dụng: Khóa luận là tài liệu học tập trực tiếp cho học phần Phương pháp toán - lý, giúp sinh viên nắm vững phương pháp giải phương trình đạo hàm riêng bằng hàm đặc biệt.
  • Giảng viên bộ môn Phương pháp toán - lý: Tài liệu cung cấp hệ thống dẫn giải giải tích chi tiết, có thể trích xuất làm bài giảng lý thuyết và bài tập mẫu về hàm trụ Bessel.
  • Phần đáng tham khảo nhất: Quy trình tách biến chuyển sang phương trình Bessel và kỹ thuật ứng dụng tính chất trực giao của hàm Bessel để tính toán hệ số chuỗi Fourier - Bessel trong Chương II và Chương III.

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

  1. Tại sao phương pháp hàm Bessel được đánh giá là ưu việt hơn phương pháp biến đổi Laplace hay tách biến Fourier thuần túy trong bài toán truyền nhiệt nhiều chiều? Theo tác giả ở phần mở đầu, khi giải các bài toán biên nhiều chiều (nhất là trong các miền có tính đối xứng trục), việc dùng phương pháp tách biến Fourier thông thường hoặc biến đổi Laplace làm cho việc tìm nghiệm gặp khó khăn và quy trình biến đổi toán học rất phức tạp. Khi chuyển bài toán sang tọa độ trụ và sử dụng hàm Bessel, cấu trúc trực giao của hàm này giúp việc xác định nghiệm giải tích và tính toán các hệ số trở nên đơn giản hơn nhiều.

  2. Vì sao trong nghiệm phương trình truyền nhiệt đối với thanh trụ, thành phần hàm Bessel loại II $Y_0(\omega r)$ lại bị triệt tiêu? Tại phương trình (3.26) của Chương III, nghiệm tổng quát của phương trình Bessel cấp không là $R(r) = C_1 J_0(\omega r) + C_2 Y_0(\omega r)$. Do bài toán xét trên toàn bộ miền thanh trụ $0 \le r \le r_0$, bao gồm cả tâm $r = 0$, mà hàm Bessel loại II $Y_0(x)$ có điểm kỳ dị tiến tới $-\infty$ khi $x \to 0$, nên để nghiệm nhiệt độ có ý nghĩa vật lý (hàm nhiệt độ phải hữu hạn trong toàn miền), tác giả phải chọn hằng số $C_2 = 0$.

  3. Hệ số $A_i$ trong nghiệm bài toán truyền nhiệt xuyên tâm được tính toán dựa trên tính chất toán học nào? Hệ số $A_i$ được xác định dựa trên tính trực giao thứ nhất của hàm Bessel loại I trên đoạn $[0, r_0]$ với hàm trọng là $r$. Khi khai triển hàm phân bố nhiệt ban đầu $f(r) = \sum_{i=1}^\infty A_i J_0(\omega_i r)$, tác giả nhân hai vế với $r J_0(\omega_j r)$ và lấy tích phân từ $0$ đến $r_0$. Nhờ tính chất trực giao, các số hạng với $i \neq j$ đều triệt tiêu, chỉ còn lại số hạng $i = j$, từ đó suy ra biểu thức tính $A_i$ qua tích phân xác định chứa $J_1^2(\mu_i^{(0)})$.

  4. Khóa luận kết hợp những định luật vật lý nào để thiết lập phương trình vi phân truyền nhiệt tổng quát? Trong Chương III, phương trình truyền nhiệt tổng quát $\text{div}(k \text{grad} u) + H = \frac{\partial}{\partial t}(c\rho u)$ được thiết lập dựa trên ba yếu tố: (1) Định luật Fourier về dòng nhiệt dẫn $\vec{q} = -k\text{grad},u$; (2) Định luật bảo toàn năng lượng viết cho thể tích $V$ bất kỳ bao quanh bởi mặt kín $S$; (3) Định lý tích phân Gauss dùng để chuyển tích phân mặt của thông lượng nhiệt sang tích phân thể tích.

  5. Phương trình khuếch tán vật chất được liên hệ như thế nào với phương trình truyền nhiệt trong khóa luận? Tác giả chỉ ra rằng phương trình khuếch tán có cấu trúc toán học tương đương với phương trình truyền nhiệt. Trong đó, nhiệt độ $u(x, y, z, t)$ được thay bằng nồng độ chất khuếch tán $C(x, y, z, t)$, và định luật Fourier được thay thế tương ứng bằng định luật Fick $\vec{J} = -D\nabla C$, với $D$ là hệ số khuếch tán. Các dạng điều kiện biên Dirichlet, Neumann và Robin của phương trình khuếch tán cũng tương thích hoàn toàn với bài toán nhiệt.

Kết luận

Khóa luận tốt nghiệp đã giải quyết trọn vẹn mục tiêu nghiên cứu thông qua việc hệ thống hóa chặt chẽ cơ sở toán học của hàm Bessel và ứng dụng hàm này để giải các bài toán truyền nhiệt trong hệ tọa độ trụ. Bằng phương pháp tách biến kết hợp với tính trực giao của hàm Bessel loại I, tác giả đã thiết lập thành công nghiệm giải tích dạng chuỗi tường minh cho cả hai trường hợp truyền nhiệt đối xứng xuyên tâm và không xuyên tâm. Công trình không chỉ cung cấp tài liệu tham khảo có giá trị học thuật cao cho sinh viên ngành Vật lý trong học phần Phương pháp toán - lý mà còn minh chứng cho tính hiệu quả của công cụ hàm đặc biệt trong việc giải quyết các bài toán biên phức tạp của vật lý toán.