Chiều sâu cần một định nghĩa đo được
Câu “khuếch tán sâu 0,5 mm” chưa đủ để tái tạo kết quả. Đó có thể là vị trí đạt một nồng độ xác định, vị trí vượt giới hạn phát hiện của phép đo hoặc chỉ là một thang chiều dài từ phân tích thứ nguyên. Những định nghĩa khác nhau cho những con số khác nhau ngay cả khi cùng dùng một phương trình và cùng bộ dữ liệu.
Bài này dựng một vật liệu đồng nhất bán vô hạn, ban đầu có nồng độ đều, rồi thay đổi nồng độ tại mặt phẳng x = 0. Các hệ số và nồng độ đều giả định; ví dụ không mô tả thời gian xử lý nhiệt hoặc thấm hóa học cho một mác thép cụ thể. Mục tiêu là nối điều kiện biên với kết quả định lượng có thể kiểm tra.
Nghiệm nào đi cùng điều kiện nào
Khi D hằng, không có phản ứng hoặc vận chuyển đối lưu, phương trình một chiều là:
∂c/∂t = D ∂²c/∂x²
Cho c(x,0) = c_0, c(0,t) = c_s với t > 0, và c tiến về c_0 khi x rất lớn. Nghiệm theo hàm sai số bù được trình bày trong bài giảng MIT.[1] Đặt f là phần tăng nồng độ so với bước tăng ở mặt:
f = (c-c_0)/(c_s-c_0)
f(x,t) = erfc[x/(2√Dt)]
D có đơn vị m²/s, x bằng m, t bằng s và đối số trong erfc không có đơn vị. erfc(0) = 1, còn erfc tiến về 0 khi đối số tăng lớn. Không thay t bằng giờ trong khi giữ D theo giây; sai số khi ấy đi vào căn bậc hai và có thể làm kết quả vẫn trông có vẻ hợp lý.
Bhadeshia chỉ ra biến tương tự x/√t làm gộp các đường nồng độ theo không gian–thời gian dưới các điều kiện phù hợp.[2] Đó là tính chất của bài toán không có một chiều dài hữu hạn chi phối, không phải quy luật áp dụng nguyên dạng cho mọi lớp phủ hoặc mẫu nhiều lớp.
Ba mức nồng độ cho ba chiều sâu
Giả định D = 1×10⁻¹¹ m²/s, c_0 = 2 và c_s = 20 mol/m³. Sau t = 3.600 s, ta có √Dt = 0,189737 mm và 2√Dt = 0,379473 mm. Với một mức f* được chọn, chiều sâu là:
x_f = 2√Dt × erfc⁻¹(f*)
erfc⁻¹ biểu thị hàm ngược, không phải nghịch đảo 1/erfc. Dùng các giá trị đối số 0,4769363; 1,1630872; 1,8213864 tương ứng f* = 0,50; 0,10; 0,01, thu được:
Tại x = 0,30 mm, f = 0,2635525 và c = 6,743945 mol/m³. Điểm này nằm sâu hơn x_50 nhưng nông hơn x_10, nhất quán với nồng độ giảm theo chiều sâu. Đây là một kiểm tra thứ tự đơn giản trước khi tin bảng tính nhiều chữ số.
Trong mô hình lý tưởng, không có một mặt sắc mà phía ngoài hoàn toàn chưa đổi nồng độ. Giá trị tăng giảm rất nhanh nhưng không bằng đúng không ở mọi x hữu hạn. Một giới hạn đo sẽ tạo “mặt phát hiện” phụ thuộc thiết bị; không nên diễn giải nó là ranh giới vật lý tuyệt đối.
Thang chiều dài không phải cùng một đại lượng
Phân tích thứ nguyên cho chiều dài tỷ lệ √Dt và thời gian tỷ lệ L²/D; COMSOL giải thích các thang này cùng vai trò của số Fourier.[3] Hệ số đứng trước căn phụ thuộc định nghĩa và điều kiện biên. Trong ví dụ, 2√Dt bằng 0,379473 mm, khác cả ba chiều sâu vừa tính.
Nếu xét riêng sự lan rộng của một phân bố hạt ban đầu tập trung trong không gian một chiều vô hạn, độ lệch chuẩn của nghiệm Gaussian là √(2Dt), ở đây bằng 0,268328 mm. Đây là một bài toán điều kiện đầu khác. Dùng trị số đó làm chiều sâu đạt 10% trong bài bước nồng độ sẽ tạo sai khác dù cả hai biểu thức đều có đúng đơn vị.
Giữ D và f* không đổi, tăng chiều sâu mục tiêu từ 0,50 lên 1,00 mm đòi thời gian gấp bốn. Với f* = 0,10, thời gian đến 0,50 mm là 4.620,144 s, khoảng 1,28337 giờ. D tăng gấp đôi làm thời gian ấy giảm một nửa; thời gian tăng gấp đôi chỉ làm chiều sâu tăng √2 lần.
Kiểm lượng tích lũy bằng thông lượng tại mặt
Lấy đạo hàm nghiệm theo x tại x = 0 rồi dùng Fick cho:
J_s(t) = (c_s-c_0)√[D/(πt)]
Tích phân phần tăng nồng độ trên nửa không gian cho lượng tăng trên đơn vị diện tích:
Q_A(t) = 2(c_s-c_0)√(Dt/π)
Sau một giờ, J_s = 5,352372×10⁻⁷ mol/(m²·s), còn Q_A = 0,003853708 mol/m². Đạo hàm Q_A theo thời gian đúng bằng J_s, nên cân bằng lượng chất được khép kín. Nhân thông lượng tại phút cuối với toàn bộ 3.600 s chỉ cho một nửa lượng đúng, vì thông lượng trước đó lớn hơn.
J_s tiến về vô hạn khi t tiến về không là hệ quả của bước nồng độ tức thời tại mặt và mô hình liên tục lý tưởng. Tổng lượng trong một khoảng thời gian hữu hạn vẫn hữu hạn. Thiết bị thật có thời gian tăng nồng độ và trở lực mặt; không lấy giá trị kỳ dị này làm thông lượng vận hành có thể đạt được.
Một mẫu hữu hạn cần kiểm mặt phía sau
Giả định bề dày thực H = 1 mm. Nghiệm bán vô hạn dự báo tại x = H sau một giờ phần tăng f chỉ khoảng 0,00019394. Nếu người phân tích tự chọn mức 1% làm tiêu chí khảo sát ảnh hưởng tới mặt sau, nghiệm chạm mức đó vào khoảng 7.535,912 s, tức 2,09331 giờ.
Mức 1% là lựa chọn minh họa, không phải tiêu chuẩn chấp nhận sai số. Nó chỉ giúp quyết định lúc nào cần giải bài toán hữu hạn với điều kiện mặt sau thật. Nếu mặt sau không thấm, được giữ ở nồng độ khác hoặc tiếp xúc với một bình nhỏ, các nghiệm dài hạn sẽ khác nhau. Một kết quả tốt ở thời gian sớm không xác nhận việc kéo dài mô hình bán vô hạn vô hạn.
Ngừng cấp ở mặt không làm bên trong dừng ngay
Xét mặt được giữ ở 20 mol/m³ trong t_p = 3.600 s rồi trở lại 2 mol/m³. Do phương trình đang tuyến tính, sau khi ngừng cấp có thể chồng một bước tăng và một bước giảm:
f_p = erfc[x/(2√Dt)] - erfc[x/(2√D(t-t_p))]
Công thức này dùng cho t > t_p; trước đó chỉ giữ số hạng thứ nhất. Tại x = 0,30 mm và t = 7.200 s, f_p = 0,1656428, nên c = 4,981571 mol/m³. Vật liệu còn chứa cấu tử dù mặt đã trở về mức ban đầu.
Đặc biệt, ở cùng chiều sâu, nồng độ vẫn tăng trong một thời gian sau khi tắt bước mặt. Giải điều kiện đạo hàm theo t bằng không cho cực đại tại t = 4.288,339 s, tức khoảng 688 s sau khi ngừng cấp; c_max = 7,311667 mol/m³. Gradient bên trong tiếp tục phân phối chất. Không suy thời điểm cực đại ở chiều sâu từ thời điểm cực đại của nồng độ mặt.
Khi D thay đổi đều theo thời gian
Nếu D chỉ phụ thuộc t, đồng nhất theo x và không phụ thuộc c, đặt Θ(t) = ∫D(t)dt. Quy tắc đạo hàm đổi phương trình thành ∂c/∂Θ = ∂²c/∂x²; nghiệm bước khi ấy dùng √Θ thay √Dt. Đây là phép đổi biến, không phải quy tắc lấy trung bình nhiệt độ.
Giả định 1.800 s đầu D = 1×10⁻¹¹ và 1.800 s sau D = 4×10⁻¹¹ m²/s. Θ bằng 9×10⁻⁸ m², nên x_10 sau một giờ bằng 0,697852 mm. Nếu dùng D ban đầu cho toàn giờ, chỉ được 0,441361 mm. Đổi nhiệt độ còn có thể đổi nồng độ mặt, hệ số phân bố, tổ chức hoặc cơ chế khuếch tán; lúc ấy phép thay một đồng hồ Θ chưa đủ. Kết quả chiều sâu chỉ có thể hiểu đúng khi giữ cùng định nghĩa nồng độ và toàn bộ lịch sử biên.
Nguồn tham khảo
[1] Archived Lecture Notes No. 9 — Diffusion, 3.091 Introduction to Solid State Chemistry — Massachusetts Institute of Technology — OpenCourseWare. Truy cập ngày 09/11/2026.
[2] The Exponential and Error Function Solutions in Diffusion — H. K. D. H. Bhadeshia — tài liệu trên máy chủ University of Cambridge. Truy cập ngày 09/11/2026.
[3] Diffusion Length and Timescales — COMSOL — Multiphysics Cyclopedia. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Câu “khuếch tán sâu 0,5 mm” chưa đủ để tái tạo kết quả. Đó có thể là vị trí đạt một nồng độ xác định, vị trí vượt giới hạn phát hiện của phép đo hoặc chỉ là một thang chiều dài từ phân tích thứ nguyên. Những định nghĩa khác nhau cho những con số khác nhau ngay cả khi cùng dùng một phương trình và cùng bộ dữ liệu.
Bài này dựng một vật liệu đồng nhất bán vô hạn, ban đầu có nồng độ đều, rồi thay đổi nồng độ tại mặt phẳng x = 0. Các hệ số và nồng độ đều giả định; ví dụ không mô tả thời gian xử lý nhiệt hoặc thấm hóa học cho một mác thép cụ thể. Mục tiêu là nối điều kiện biên với kết quả định lượng có thể kiểm tra.
Nghiệm nào đi cùng điều kiện nào
Khi D hằng, không có phản ứng hoặc vận chuyển đối lưu, phương trình một chiều là:
∂c/∂t = D ∂²c/∂x²
Cho c(x,0) = c_0, c(0,t) = c_s với t > 0, và c tiến về c_0 khi x rất lớn. Nghiệm theo hàm sai số bù được trình bày trong bài giảng MIT.[1] Đặt f là phần tăng nồng độ so với bước tăng ở mặt:
f = (c-c_0)/(c_s-c_0)
f(x,t) = erfc[x/(2√Dt)]
D có đơn vị m²/s, x bằng m, t bằng s và đối số trong erfc không có đơn vị. erfc(0) = 1, còn erfc tiến về 0 khi đối số tăng lớn. Không thay t bằng giờ trong khi giữ D theo giây; sai số khi ấy đi vào căn bậc hai và có thể làm kết quả vẫn trông có vẻ hợp lý.
Bhadeshia chỉ ra biến tương tự x/√t làm gộp các đường nồng độ theo không gian–thời gian dưới các điều kiện phù hợp.[2] Đó là tính chất của bài toán không có một chiều dài hữu hạn chi phối, không phải quy luật áp dụng nguyên dạng cho mọi lớp phủ hoặc mẫu nhiều lớp.
Ba mức nồng độ cho ba chiều sâu
Giả định D = 1×10⁻¹¹ m²/s, c_0 = 2 và c_s = 20 mol/m³. Sau t = 3.600 s, ta có √Dt = 0,189737 mm và 2√Dt = 0,379473 mm. Với một mức f* được chọn, chiều sâu là:
x_f = 2√Dt × erfc⁻¹(f*)
erfc⁻¹ biểu thị hàm ngược, không phải nghịch đảo 1/erfc. Dùng các giá trị đối số 0,4769363; 1,1630872; 1,8213864 tương ứng f* = 0,50; 0,10; 0,01, thu được:
- Tăng bằng 50% bước mặt: x_50 = 0,180985 mm, c = 11 mol/m³.
- Tăng bằng 10% bước mặt: x_10 = 0,441361 mm, c = 3,8 mol/m³.
- Tăng bằng 1% bước mặt: x_1 = 0,691168 mm, c = 2,18 mol/m³.
Tại x = 0,30 mm, f = 0,2635525 và c = 6,743945 mol/m³. Điểm này nằm sâu hơn x_50 nhưng nông hơn x_10, nhất quán với nồng độ giảm theo chiều sâu. Đây là một kiểm tra thứ tự đơn giản trước khi tin bảng tính nhiều chữ số.
Trong mô hình lý tưởng, không có một mặt sắc mà phía ngoài hoàn toàn chưa đổi nồng độ. Giá trị tăng giảm rất nhanh nhưng không bằng đúng không ở mọi x hữu hạn. Một giới hạn đo sẽ tạo “mặt phát hiện” phụ thuộc thiết bị; không nên diễn giải nó là ranh giới vật lý tuyệt đối.
Thang chiều dài không phải cùng một đại lượng
Phân tích thứ nguyên cho chiều dài tỷ lệ √Dt và thời gian tỷ lệ L²/D; COMSOL giải thích các thang này cùng vai trò của số Fourier.[3] Hệ số đứng trước căn phụ thuộc định nghĩa và điều kiện biên. Trong ví dụ, 2√Dt bằng 0,379473 mm, khác cả ba chiều sâu vừa tính.
Nếu xét riêng sự lan rộng của một phân bố hạt ban đầu tập trung trong không gian một chiều vô hạn, độ lệch chuẩn của nghiệm Gaussian là √(2Dt), ở đây bằng 0,268328 mm. Đây là một bài toán điều kiện đầu khác. Dùng trị số đó làm chiều sâu đạt 10% trong bài bước nồng độ sẽ tạo sai khác dù cả hai biểu thức đều có đúng đơn vị.
Giữ D và f* không đổi, tăng chiều sâu mục tiêu từ 0,50 lên 1,00 mm đòi thời gian gấp bốn. Với f* = 0,10, thời gian đến 0,50 mm là 4.620,144 s, khoảng 1,28337 giờ. D tăng gấp đôi làm thời gian ấy giảm một nửa; thời gian tăng gấp đôi chỉ làm chiều sâu tăng √2 lần.
Kiểm lượng tích lũy bằng thông lượng tại mặt
Lấy đạo hàm nghiệm theo x tại x = 0 rồi dùng Fick cho:
J_s(t) = (c_s-c_0)√[D/(πt)]
Tích phân phần tăng nồng độ trên nửa không gian cho lượng tăng trên đơn vị diện tích:
Q_A(t) = 2(c_s-c_0)√(Dt/π)
Sau một giờ, J_s = 5,352372×10⁻⁷ mol/(m²·s), còn Q_A = 0,003853708 mol/m². Đạo hàm Q_A theo thời gian đúng bằng J_s, nên cân bằng lượng chất được khép kín. Nhân thông lượng tại phút cuối với toàn bộ 3.600 s chỉ cho một nửa lượng đúng, vì thông lượng trước đó lớn hơn.
J_s tiến về vô hạn khi t tiến về không là hệ quả của bước nồng độ tức thời tại mặt và mô hình liên tục lý tưởng. Tổng lượng trong một khoảng thời gian hữu hạn vẫn hữu hạn. Thiết bị thật có thời gian tăng nồng độ và trở lực mặt; không lấy giá trị kỳ dị này làm thông lượng vận hành có thể đạt được.
Một mẫu hữu hạn cần kiểm mặt phía sau
Giả định bề dày thực H = 1 mm. Nghiệm bán vô hạn dự báo tại x = H sau một giờ phần tăng f chỉ khoảng 0,00019394. Nếu người phân tích tự chọn mức 1% làm tiêu chí khảo sát ảnh hưởng tới mặt sau, nghiệm chạm mức đó vào khoảng 7.535,912 s, tức 2,09331 giờ.
Mức 1% là lựa chọn minh họa, không phải tiêu chuẩn chấp nhận sai số. Nó chỉ giúp quyết định lúc nào cần giải bài toán hữu hạn với điều kiện mặt sau thật. Nếu mặt sau không thấm, được giữ ở nồng độ khác hoặc tiếp xúc với một bình nhỏ, các nghiệm dài hạn sẽ khác nhau. Một kết quả tốt ở thời gian sớm không xác nhận việc kéo dài mô hình bán vô hạn vô hạn.
Ngừng cấp ở mặt không làm bên trong dừng ngay
Xét mặt được giữ ở 20 mol/m³ trong t_p = 3.600 s rồi trở lại 2 mol/m³. Do phương trình đang tuyến tính, sau khi ngừng cấp có thể chồng một bước tăng và một bước giảm:
f_p = erfc[x/(2√Dt)] - erfc[x/(2√D(t-t_p))]
Công thức này dùng cho t > t_p; trước đó chỉ giữ số hạng thứ nhất. Tại x = 0,30 mm và t = 7.200 s, f_p = 0,1656428, nên c = 4,981571 mol/m³. Vật liệu còn chứa cấu tử dù mặt đã trở về mức ban đầu.
Đặc biệt, ở cùng chiều sâu, nồng độ vẫn tăng trong một thời gian sau khi tắt bước mặt. Giải điều kiện đạo hàm theo t bằng không cho cực đại tại t = 4.288,339 s, tức khoảng 688 s sau khi ngừng cấp; c_max = 7,311667 mol/m³. Gradient bên trong tiếp tục phân phối chất. Không suy thời điểm cực đại ở chiều sâu từ thời điểm cực đại của nồng độ mặt.
Khi D thay đổi đều theo thời gian
Nếu D chỉ phụ thuộc t, đồng nhất theo x và không phụ thuộc c, đặt Θ(t) = ∫D(t)dt. Quy tắc đạo hàm đổi phương trình thành ∂c/∂Θ = ∂²c/∂x²; nghiệm bước khi ấy dùng √Θ thay √Dt. Đây là phép đổi biến, không phải quy tắc lấy trung bình nhiệt độ.
Giả định 1.800 s đầu D = 1×10⁻¹¹ và 1.800 s sau D = 4×10⁻¹¹ m²/s. Θ bằng 9×10⁻⁸ m², nên x_10 sau một giờ bằng 0,697852 mm. Nếu dùng D ban đầu cho toàn giờ, chỉ được 0,441361 mm. Đổi nhiệt độ còn có thể đổi nồng độ mặt, hệ số phân bố, tổ chức hoặc cơ chế khuếch tán; lúc ấy phép thay một đồng hồ Θ chưa đủ. Kết quả chiều sâu chỉ có thể hiểu đúng khi giữ cùng định nghĩa nồng độ và toàn bộ lịch sử biên.
Nguồn tham khảo
[1] Archived Lecture Notes No. 9 — Diffusion, 3.091 Introduction to Solid State Chemistry — Massachusetts Institute of Technology — OpenCourseWare. Truy cập ngày 09/11/2026.
[2] The Exponential and Error Function Solutions in Diffusion — H. K. D. H. Bhadeshia — tài liệu trên máy chủ University of Cambridge. Truy cập ngày 09/11/2026.
[3] Diffusion Length and Timescales — COMSOL — Multiphysics Cyclopedia. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
