7. Bố cục của luận án
2.2.2 Các bước trong tính toán lý thuyết và các giả thiết khi thực hiện tính toán
Phân phối năng lượng electron đặc trưng theo thực nghiệm “Steady-State
Townsend” là: F (ε) +αF (ε) +α2F (ε) +α3F (ε) +... F (ε) = 0 1 s 1+ α 2 A +α2A 3 +α3A +... 1 2 3 (2.26) Khi đó: Tần số ion hóa là: R ∞ i = V1N∫ 0 ε1/2Q F dε. (2.26a) Tần số kết hợp là: ∞ Ra = V1N∫ 0 ε1/2Q F dε. (2.26b)
2.2.2 Các bước trong tính toán lý thuyết và các giả thiết khi thực hiện tính toán
Phương pháp hữu hạn đã được sử dụng để xác định các tham số [6, 74]: -Các hệ số Ak (k = 0, 1, 2, 3…) được xác định từ giải Fk(ε);
- Các điều kiện: F0(0) và F0(∞) = 0 và fn(0) = 0 cho n ≥ 1[79]với fn(ε) = Fn(ε)ε-1/2.
Hàm phân phối năng lượng electron fn (ε) được tính theo dải năng lượng của electron (ε) có giá trị 0÷ (εmax) trong đó, f0(εmax)/f0(ε)max<10-10, và các khoảng phân chia năng lượng Δε được xác định xấp xỉ sao cho thỏa mãn mứccân bằng năng lượng trong khoảng 0,1%. Các giá trị của εmax và Δε thay đổi tùy thuộc vào các giá trị E/N khác nhau.
- Các điều kiện ban đầu và dữ liệu đầu vào để tính toán là: i s
+ Điều kiện ban đầu: mật độ khối khí N = 3,5353 × 1016cm-3 (ở 00C và áp suất là 1Torr) với tỷ lệ phân chia năng lượng còn lại sau khi va chạm ion hóa là 0,5;
+ Dữ liệu đầu vào: Các bộ tiết diện va chạm của electron các chất khí xem xét, nhiệt độ khí T, giá trị của E/N, năng lượng ε, và tỷ lệ chia vùng năng lượng electron từ 0 ÷ εmax.
Lưu đồ thuật toán phục vụ việc xác định các hệ số chuyển động của electron sử dụng phương trình xấp xỉ bậc hai Boltzmann được thể hiện trên hình 2.4. Bảng hiện thị kết quả tính toán sử dụng thuật toán phương trình xấp xỉ bậc hai Boltzmann được lấy ra từ kết quả chạy chương trình trình bày trong phụ lục 1.
Hình 2.4 Lưu đồ thuật toán tính toán hàm phân bố năng lượng của phương trình xấp xỉ bậc hai Boltzmann
Với thuật toán hàm phân bố năng lượng của phương trình xấp xỉ bậc hai Boltzmann trong luận án, nghiên cứu sinh có sử dụng và tính toán được hàm phân phối năng lượng electron fn (ε) được tính trong dải năng lượng của electron với giá trị từ 0 ÷ (εmax); trong đó: f0 (εmax )/f0 (ε)max<10-10. Kết quả thể hiện tính toán được thể hiện như phụ lục 1.