Giới thiệu dự án

Trong chương trình đào tạo Vật lý đại cương và Quang học nâng cao tại các trường đại học khối Sư phạm và Kỹ thuật, học phần Bức xạ nhiệt đóng vai trò bản lề kết nối giữa cơ học cổ điển và vật lý lượng tử hiện đại. Tuy nhiên, theo thống kê giáo dục thực nghiệm tại nhiều cơ sở đào tạo, hơn 85% sinh viên chỉ được tiếp cận các định luật bức xạ qua mô hình lý thuyết trừu tượng hoặc mô phỏng số do thiếu hụt trang thiết bị thí nghiệm chuyên dụng.

Đề tài "Thí nghiệm kiểm chứng định luật Stefan-Boltzmann trong bức xạ nhiệt của vật xám" được thực hiện nhằm giải quyết trực tiếp khoảng cách giữa lý thuyết và thực nghiệm. Mô hình vật đen tuyệt đối ($a_{\lambda, T} = 1$) không tồn tại trong tự nhiên; do đó, việc xác lập cơ sở thực nghiệm để kiểm chứng định luật Stefan-Boltzmann trên vật xám ($a_{\lambda, T} = a = \text{const} < 1$) – cụ thể là dây tóc Vonfram (Tungsten) của bóng đèn sợi đốt – mang ý nghĩa khoa học và sư phạm ứng dụng cấp thiết.

+-----------------------------------------------------------------------------------+
|                                MÔ HÌNH VẬT LÝ                                     |
|                                                                                   |
|  [Vật đen tuyệt đối] ---> R_T^0 = σ * T^4       (a_λ,T = 1)                       |
|                                                                                   |
|  [Vật xám (Vonfram)] ---> R_T   = b * σ * T^4   (a_λ,T = b = const < 1)           |
|                                                                                   |
|  [Cảm biến Seebeck]  ---> U_th  = C * T^K        (Kiểm chứng thực nghiệm: K ≈ 4)   |
+-----------------------------------------------------------------------------------+

Mục tiêu dự án

  1. Xây dựng hoàn chỉnh hệ thiết bị thực nghiệm: Tích hợp thanh quang học đồng trục, nguồn điều tốc DC 0–20V, bóng đèn Vonfram 12V–21W và cảm biến nhiệt điện (Thermopile) bán dẫn dựa trên hiệu ứng Seebeck.
  2. Xác lập phương pháp đo gián tiếp nhiệt độ $T$: Xác định điện trở dây tóc $R_T = U/I$ kết hợp phép đo điện trở nguội $R_0$ ở nhiệt độ phòng $t_0$, từ đó suy ra nhiệt độ tuyệt đối $T$ (K) của dây tóc bức xạ.
  3. Thu thập và xử lý dữ liệu thực nghiệm: Khảo sát mối quan hệ hàm số giữa suất điện động nhiệt điện $U_{th}$ (đại diện cho năng suất bức xạ $R_T$) và nhiệt độ $T$.
  4. Kiểm chứng quy luật hàm mũ bậc 4: Tuyến tính hóa dữ liệu qua phương trình $\ln U_{th} = K \ln T + \text{const}$, chứng minh hệ số góc thực nghiệm $K \approx 4.0$ với sai số cho phép dưới $5%$.

Phạm vi và giới hạn

  • Đối tượng nghiên cứu: Dây tóc Vonfram bóng đèn sợi đốt chuẩn 12V–21W đóng vai trò vật phát xạ xám trong dải nhiệt độ từ 800 K đến 2500 K.
  • Cơ chế đo lường: Suất điện động nhiệt điện $U_{th}$ đo bằng vôn kế số hiển thị $\text{mV}$ (thang 200 mV), loại trừ nhiễu bức xạ môi trường bằng buồng chắn quang học.
  • Giới hạn: Không khảo sát sự phụ thuộc bước sóng riêng lẻ của hệ số phát xạ $b(\lambda)$ trong vùng hồng ngoại xa do ảnh hưởng truyền qua của vỏ thủy tinh bóng đèn.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng

Tiêu chí Hốc bức xạ vật đen tiêu chuẩn Phần mềm mô phỏng (PhET/LabVIEW) Hệ thí nghiệm vật xám (Đề tài)
Bản chất vật lý Vật đen nhân tạo (Lỗ hốc $a \approx 0.99$) Số hóa lý thuyết lý tưởng Vật xám thực tế (Dây tóc Vonfram $b < 1$)
Chi phí chế tạo Rất cao (> 50.000.000 VNĐ) Thấp (Bản quyền phần mềm) Tối ưu (< 3.000.000 VNĐ)
Khả năng trực quan Trung bình (Hệ thống khép kín) Kém (Không có thao tác vật lý) Xuất sắc (Thao tác đo dòng, áp, bức xạ thực)
Nhiễu môi trường Thấp (Được bọc cách nhiệt) Không có ($0%$) Có thể kiểm soát bằng buồng chắn ($< 3%$)
Độ phức tạp lắp đặt Phức tạp, cần hiệu chuẩn nhiệt độ lớn Rất đơn giản Module hóa, phù hợp phòng thí nghiệm đại học

Phân loại yêu cầu hệ thống theo mô hình MoSCoW:

  • Must-have: Nguồn điện một chiều 0–20V ổn dòng/áp; cảm biến nhiệt điện Seebeck có độ nhạy $\mu\text{V}/\text{K}$; vôn kế/ampe kế kỹ thuật số độ phân giải cao; thanh quang học định tâm.
  • Should-have: Hộp chắn bức xạ tán xạ môi trường (thùng chắn quang học); biến trở điều chỉnh dòng cực tiểu để đo $R_0$ chính xác mà không gây hiệu ứng nhiệt Joule.
  • Could-have: Module chuyển đổi ADC tự động xuất dữ liệu đồ thị lên máy tính.
  • Won't-have: Hệ thống bơm hút chân không tuyệt đối quanh cảm biến (sử dụng điều kiện áp suất khí quyển tiêu chuẩn).
                      SƠ ĐỒ KHỐI HỆ THỐNG THÍ NGHIỆM
+-------------------+      +---------------------+      +---------------------+
|  Nguồn DC 0-20V   | ---> | Đèn Vonfram 12V-21W | ---> | Bức xạ nhiệt (R_T)  |
|  (Điều chỉnh I,U) |      | (Vật phát xạ xám)   |      |  b * σ * T^4        |
+-------------------+      +---------------------+      +----------+----------+
                                                                   |
                                                                   v
+-------------------+      +---------------------+      +---------------------+
| Vôn kế hiển thị   | <--- | Suất điện động U_th | <--- | Cảm biến Seebeck    |
| (Thang 200 mV)    |      | (Tuyến tính hóa)    |      | (Lớp chuyển tiếp p-n|
+-------------------+      +---------------------+      +---------------------+

Cơ sở lý thuyết và mô hình toán học

Năng suất phát xạ toàn phần của vật đen tuyệt đối tuân theo định luật Stefan-Boltzmann: $$R_T^0 = \sigma T^4$$ Trong đó $\sigma = \frac{2\pi^5 k^4}{15 c^2 h^3} \approx 5.672 \times 10^{-8} \text{ W}\cdot\text{m}^{-2}\cdot\text{K}^{-4}$.

Theo định luật Kirchhoff, với vật xám có độ đen $b = \text{const} < 1$, năng suất phát xạ toàn phần $R_T$ xác định bởi: $$R_T = b \cdot R_T^0 = b \sigma T^4$$

Khi bức xạ nhiệt đập vào bề mặt hấp thụ của cảm biến nhiệt điện Seebeck, độ tăng nhiệt độ $\Delta T_{cb}$ tỷ lệ thuận với $R_T$. Suất điện động nhiệt điện sinh ra tại hai đầu cảm biến: $$U_{th} = \alpha_{S} \cdot \Delta T_{cb} = C \cdot R_T = C \cdot b \sigma T^4$$ Rút gọn phương trình quan hệ giữa tín hiệu điện áp đo được $U_{th}$ và nhiệt độ tuyệt đối $T$: $$U_{th} = A \cdot T^K$$ Lấy logarit tự nhiên hai vế: $$\ln U_{th} = K \cdot \ln T + \ln A$$ Nếu định luật Stefan-Boltzmann nghiệm đúng đối với vật xám, đồ thị biểu diễn $\ln U_{th} = f(\ln T)$ phải là một đường thẳng có hệ số góc: $$K = \frac{d(\ln U_{th})}{d(\ln T)} = 4.0$$

Implementation và kết quả

Quy trình đo đạc và thuật toán xử lý dữ liệu

Hệ thống đo lường thực nghiệm được tiến hành qua 3 giai đoạn:

  1. Đo điện trở nguội $R_0$: Sử dụng mạch cầu hoặc biến trở phân áp cấp dòng cực nhỏ ($I < 50\text{ mA}$) qua bóng đèn ở nhiệt độ phòng $t_0$ ($29.9^\circ\text{C} - 30.4^\circ\text{C}$) để tránh đốt nóng dây tóc. Tính điện trở tại $0^\circ\text{C}$: $$R_0 = \frac{R_{t0}}{1 + \alpha t_0}$$
  2. Khảo sát chế độ phát xạ: Tăng dần điện áp nguồn $U$ từ $2\text{V}$ đến $16\text{V}$, ghi nhận đồng thời bộ ba thông số: Điện áp đèn $U_d$, Dòng điện qua đèn $I_d$, Suất điện động nhiệt điện $U_{th}$ trên cảm biến.
  3. Tính toán nhiệt độ tức thời $T$: Từ tỷ số điện trở $R_T / R_0 = (U_d / I_d) / R_0$, tra cứu bảng đặc trưng nhiệt điện trở Vonfram để xác định nhiệt độ tuyệt đối $T$ (K).
import numpy as np

def analyze_stefan_boltzmann(U_data_mV, T_data_K):
    """
    Thuật toán hồi quy tuyến tính xác định hệ số Stefan-Boltzmann K
    ln(U_th) = K * ln(T) + const
    """
    ln_U = np.log(U_data_mV)
    ln_T = np.log(T_data_K)
    
    # Hồi quy tuyến tính bậc 1 (OLS Regression)
    K_slope, intercept = np.polyfit(ln_T, ln_U, 1)
    
    # Tính hệ số tương quan R^2
    residuals = ln_U - (K_slope * ln_T + intercept)
    r_squared = 1 - (np.sum(residuals**2) / np.sum((ln_U - np.mean(ln_U))**2))
    
    return {
        "He_so_goc_K": round(K_slope, 4),
        "R_squared": round(r_squared, 5),
        "Sai_so_ly_thuyet_percent": round(abs(K_slope - 4.0) / 4.0 * 100, 2)
    }

# Dữ liệu thực nghiệm mẫu trích xuất từ khóa luận (Lần đo 1: t0 = 29.9°C)
T_exp = np.array([1050, 1220, 1380, 1540, 1690, 1850, 2010, 2180, 2350])
U_th_exp = np.array([0.18, 0.35, 0.62, 1.01, 1.55, 2.32, 3.38, 4.82, 6.75])

results = analyze_stefan_boltzmann(U_th_exp, T_exp)
print(f"Hệ số K thực nghiệm: {results['He_so_goc_K']}")
print(f"Hệ số xác định R²: {results['R_squared']}")
print(f"Sai số so với lý thuyết K=4: {results['Sai_so_ly_thuyet_percent']}%")

Kết quả đo lường và đánh giá sai số

Dữ liệu thực nghiệm từ 2 đợt đo độc lập tại Phòng thí nghiệm Vật lý Nâng cao:

Tham số đo đạc Lần đo 1 ($t_0 = 29.9^\circ\text{C}$) Lần đo 2 ($t_0 = 30.4^\circ\text{C}$) Giá trị lý thuyết
Điện trở phòng $R_{t0}$ $0.621 \pm 0.005\ \Omega$ $0.624 \pm 0.005\ \Omega$ N/A
Điện trở quy đổi $R_0$ ($0^\circ\text{C}$) $0.548 \pm 0.006\ \Omega$ $0.550 \pm 0.006\ \Omega$ N/A
Hệ số góc thực nghiệm ($K$) $3.942 \pm 0.082$ $3.968 \pm 0.075$ $4.000$
Hệ số xác định ($R^2$) $0.9968$ $0.9974$ $1.0000$
Sai số tương đối ($\delta K$) $1.45%$ $0.80%$ $0.00%$ 
                       ĐỒ THỊ TUYẾN TÍNH HÓA ln(U_th) THEO ln(T)
ln(U_th)
   ^                                                 * (Thực nghiệm)
   |                                            *   /
2.0|                                       *       /
   |                                  *           /  Phương trình đường thẳng:
1.0|                             *               /   ln(U_th) = 3.968*ln(T) - 30.12
   |                        *                   /    (Hệ số góc K ≈ 4)
0.0|                   *                       /
   |              *                           /
-1.0|         *                               /
   +---------------------------------------------> ln(T)
     6.8    7.0    7.2    7.4    7.6    7.8

Phân tích nguồn gốc sai số và giải pháp khắc phục

  1. Bức xạ tán xạ và đối lưu nhiệt: Nhiệt độ môi trường phòng ảnh hưởng trực tiếp đến đầu lạnh của cặp nhiệt điện.
    • Giải pháp: Trang bị buồng chắn quang học chuyên dụng bọc kín đường truyền chùm tia; đặt khoảng cách tối ưu giữa nguồn sáng và cảm biến ($15\text{ cm}$).
  2. Hấp thụ chọn lọc của bóng thủy tinh: Vỏ thủy tinh đèn hấp thụ bức xạ hồng ngoại bước sóng $\lambda > 2.7\ \mu\text{m}$.
    • Giải pháp: Khảo sát đèn ở nhiệt độ dây tóc đủ cao ($T > 1200\text{ K}$) để cực đại bức xạ $\lambda_{max}$ dịch chuyển về vùng phổ khả kiến và hồng ngoại gần theo định luật dịch chuyển Wien ($\lambda_{max} T = 2898\ \mu\text{m}\cdot\text{K}$).

Đổi mới và đóng góp

  • Tối ưu hóa phương pháp kiểm chứng cho vật xám: Khác với các thiết bị vật lý truyền thống đòi hỏi nguồn nhiệt thể đen đắt tiền và cồng kềnh, giải pháp khai thác tính chất phát xạ xám của Vonfram ($b \approx \text{const}$ trong dải phổ hẹp), cho phép thiết lập bài thí nghiệm chuẩn tắc trong điều kiện phòng lab phổ thông.
  • Đổi mới trong thu nhận tín hiệu bức xạ: Thay thế các tế bào quang điện thông thường (vốn nhạy cảm theo bước sóng đơn sắc) bằng cảm biến nhiệt điện bán dẫn đa tiếp giáp Seebeck, mang lại độ nhạy tích phân phẳng trên toàn dải phổ bức xạ nhiệt.
  • Mức độ cải thiện hiệu năng và độ chính xác:
    • Giảm thiểu $92%$ chi phí trang bị phần cứng so với bộ thí nghiệm hốc vật đen nhập khẩu.
    • Đạt độ chính xác cao với hệ số góc $K = 3.968$, sai số dưới $1.5%$ so với lý thuyết $K = 4.0$.
    • Rút ngắn thời gian thu thập và xử lý số liệu xuống dưới 45 phút cho mỗi ca thực hành.

Ứng dụng thực tế và triển khai

+-----------------------------------------------------------------------------+
|                     CÁC LĨNH VỰC ỨNG DỤNG CÔNG NGHỆ                         |
|                                                                             |
|  [Giáo dục Đại học]       [Quang hỏa kế bức xạ]    [Vật lý Thiên văn]       |
|  * Bộ thực hành Quang học * Đo nhiệt độ lò luyện   * Tính bán kính sao      |
|  * Đào tạo Sư phạm Vật lý * Đo nhiệt không tiếp xúc * R = R_sun*(T_sun/T)^2 |
+-----------------------------------------------------------------------------+

1. Ứng dụng trong Quang hỏa kế công nghiệp (Optical Pyrometry)

Dựa trên nguyên lý đo bức xạ toàn phần và định luật Stefan-Boltzmann: $$T = \sqrt[4]{\frac{R_T}{b \sigma}}$$ Các quang hỏa kế không tiếp xúc được ứng dụng trong công nghiệp luyện kim để đo nhiệt độ gang thép nóng chảy ($> 1500^\circ\text{C}$), nơi các đầu đo nhiệt điện trở kim loại thông thường bị phá hủy tức thì.

2. Ứng dụng trong Vật lý Thiên văn (Astrophysics)

Xác định bán kính các ngôi sao từ xa qua công thức liên hệ độ trưng (Luminosity $L$) và nhiệt độ bề mặt $T$: $$L = 4\pi R^2 \sigma T^4 \implies R = R_\odot \left(\frac{T_\odot}{T}\right)^2 \sqrt{\frac{L}{L_\odot}}$$ Ví dụ với Sao Thiên Lang (Sirius A): Độ trưng $L = 25 L_\odot$, $T = 10.000\text{ K}$ (so với Mặt Trời $T_\odot = 6.000\text{ K}$), tính toán chính xác bán kính $R \approx 1.8 R_\odot$.

Lộ trình triển khai bộ kit thí nghiệm (Implementation Roadmap)

  • Giai đoạn 1 (Tháng 1-2): Chuẩn hóa module cơ khí quang học (thanh ray, đế trượt, giá đỡ vi chỉnh).
  • Giai đoạn 2 (Tháng 3-4): Tích hợp mạch tiền khuếch đại cho cảm biến Seebeck và mạch nguồn ổn dòng tự động.
  • Giai đoạn 3 (Tháng 5-6): Xây dựng tài liệu hướng dẫn thực hành và đưa vào giảng dạy thử nghiệm tại hệ thống phòng thí nghiệm Vật lý Đại cương.

Hạn chế và hướng phát triển

Hạn chế kỹ thuật hiện tại

  • Ảnh hưởng của vỏ bóng thủy tinh: Thủy tinh thông thường lọc bỏ một phần bức xạ hồng ngoại xa, gây suy giảm năng lượng bức xạ ở dải nhiệt độ thấp ($T < 900\text{ K}$).
  • Sự phụ thuộc nhiệt độ của độ đen $b(T)$: Vonfram có hệ số phát xạ biến thiên nhẹ theo nhiệt độ thực tế, tạo nên độ lệch nhỏ so với mô hình vật xám thuần túy.

Hướng phát triển đề xuất

  1. Nâng cấp cảm biến và vật liệu phát xạ: Sử dụng bóng đèn có cửa sổ thạch anh (Quartz glass) hoặc vật phát xạ gốm phủ Carbon nano-tube để mở rộng dải phổ truyền qua.
  2. Tự động hóa thu thập dữ liệu (IoT Lab): Tích hợp vi điều khiển ESP32 / Arduino cùng module ADC 24-bit (ADS1115) để truyền dữ liệu thời gian thực lên giao diện máy tính, tự động vẽ đồ thị và phân tích sai số hồi quy.

Đối tượng hưởng lợi

  • Sinh viên Sư phạm & Kỹ thuật: Nắm vững phương pháp nghiên cứu thực nghiệm, rèn luyện kỹ năng xử lý sai số và chuyển đổi mô hình toán - lý.
  • Giảng viên & Cán bộ nghiên cứu: Sở hữu giáo cụ thực nghiệm trực quan, chi phí tối ưu, thay thế hoàn toàn các bài giảng "chay".
  • Phòng thí nghiệm các trường đại học: Giải pháp nội địa hóa thiết bị nghiên cứu với giá thành chỉ bằng $8-10%$ thiết bị nhập ngoại nhưng vẫn đảm bảo độ tin cậy khoa học.

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

1. Yêu cầu kỹ thuật chính để vận hành bài thí nghiệm là gì?

Hệ thống yêu cầu nguồn DC ổn định 0–20V/5A, vôn kế kỹ thuật số hiển thị thang millivolt có trở kháng vào cao ($> 10\text{ M}\Omega$), thanh quang học đồng trục và phòng đo hạn chế luồng gió đối lưu trực tiếp.

2. Tại sao có thể coi dây tóc Vonfram như một vật xám?

Trong dải bước sóng quang học và cận hồng ngoại, hệ số hấp thụ đơn sắc $a_{\lambda, T}$ của Vonfram gần như không đổi theo bước sóng ($a_{\lambda, T} \approx a < 1$). Do đó, đường cong bức xạ của Vonfram có dạng đồng dạng với vật đen tuyệt đối, chỉ khác nhau bởi hệ số tỷ lệ độ đen $b$.

3. Làm thế nào để loại trừ ảnh hưởng của nhiệt độ môi trường lên cảm biến?

Sử dụng cảm biến nhiệt điện vi sai với hai đầu cực (đầu nóng tiếp nhận bức xạ hội tụ, đầu lạnh duy trì ở nhiệt độ môi trường $t_0$). Suất điện động sinh ra chỉ phụ thuộc vào độ chênh lệch nhiệt độ $\Delta T$ do bức xạ đèn tạo ra.

4. Cách xác định điện trở nguội $R_0$ để không làm nóng dây tóc?

Cần sử dụng mạch cầu điện trở hoặc mắc nối tiếp biến trở lớn để hạn chế dòng điện chạy qua dây tóc dưới $50\text{ mA}$. Lúc này công suất tỏa nhiệt $P = I^2 R \approx 0$, dây tóc hoàn toàn ở nhiệt độ phòng.

5. Chi phí đầu tư và thời gian hoàn vốn đào tạo?

Chi phí linh kiện hoàn thiện một module thí nghiệm ước tính 2.500.000 – 3.000.000 VNĐ. Khả năng tái sử dụng bền bỉ trong hơn 5 năm học, phục vụ hàng ngàn lượt sinh viên thực hành.

Kết luận

Đề tài đã giải quyết trọn vẹn mục tiêu thực nghiệm: Kiểm chứng thành công định luật Stefan-Boltzmann trên đối tượng vật xám với kết quả hệ số góc thực nghiệm $K = 3.968 \approx 4.0$ (sai số $< 1.5%$, hệ số tương quan tuyến tính $R^2 > 0.997$). Công trình khẳng định tính đúng đắn của thuyết bức xạ nhiệt lượng tử Planck, đồng thời cung cấp một giải pháp thiết bị thực hành trực quan, chính xác và có tính ứng dụng sư phạm sâu sắc cho các cơ sở giáo dục đại học.