Một trị số hydro cần có cơ sở đo
Khi báo cáo hydro trong kim loại hàn, cần biết giá trị nói về phần có thể khuếch tán ra trong điều kiện thử nào, khối lượng kim loại nào và cơ sở quy đổi thể tích khí. TWI phân biệt hydro khuếch tán với phần hydro còn bị giữ tại các vị trí trong kim loại, đồng thời mô tả cách báo thể tích khí quy về điều kiện chuẩn trên100 g kim loại đắp. [1]
Bài này không suy chế độ xử lý nhiệt từ một trị số hydro. Nó xét bài toán truyền khối lý tưởng để hiểu thời gian, chiều dày và điều kiện bề mặt ảnh hưởng lượng còn lại ra sao. Hệ số khuếch tán và mọi kích thước ví dụ đều tự đặt, không đại diện một mác thép hoặc quy trình hàn đã xác nhận.
Mô hình tấm thoát qua hai mặt
Xét tấm phẳng rất rộng, chiều dày2L, tọa độ x tính từ mặt giữa, −L≤x≤L. Ban đầu nồng độ c bằng c₀ đồng đều. Từ t=0, hai mặt được giữ c=0; bỏ trở lực phản ứng bề mặt, bẫy hydro, ứng suất, cạnh tấm và biến đổi thể tích.
Chọn c có đơn vị mol/m³ và hệ số khuếch tán D không đổi, m²/s. Phương trình một chiều là:
∂c/∂t=D∂²c/∂x².
Điều kiện c=0 là một bể thu lý tưởng, không có nghĩa mọi bề mặt thép ngoài không khí luôn thỏa điều kiện ấy. Nếu lớp phủ, phản ứng tái hợp hoặc môi trường cản trở thoát khí, cần điều kiện biên truyền khối khác.
Bài giảng MIT dùng tách biến và chuỗi Fourier cho vật hữu hạn, đồng thời nhấn mạnh điều kiện không thông lượng sẽ bảo toàn tổng lượng chất. [2] Ở đây ta tự áp dụng cùng phương pháp với hai mặt nồng độ bằng0, nên tổng lượng trong tấm giảm theo thời gian.
Tự dựng các mode khuếch tán
Do đối xứng, dùng các hàm cos(λ_nx). Điều kiện tại x=±L yêu cầu λ_n=(2n+1)π/(2L), với n=0,1,2,… Hệ số từ khai triển nồng độ đầu đồng đều là A_n=4(−1)^n/[π(2n+1)].
c/c₀=ΣA_n cos(λ_nx)exp(−Dλ_n²t).
Mỗi mode có tốc độ giảm khác nhau; mode biến thiên nhanh theo không gian mất nhanh hơn. Tích phân theo chiều dày rồi chia2L cho tỷ phần lượng còn lại R=c_trung bình/c₀:
R=(8/π²)Σ exp[−(2n+1)²π²Fo/4]/(2n+1)².
Fo=Dt/L² là thời gian không thứ nguyên. Chuỗi cộng từ n=0 tới vô hạn. Tại t=0, tổng hệ số trung bình bằng1; tại t lớn, mọi số hạng về0. Tại mặt giữa, cos bằng1 nhưng hệ số còn dấu luân phiên, nên nồng độ tâm khác trung bình.
Ví dụ một giờ và hai giờ
Chọn2L=10,0 mm, tức L=0,00500 m, D=1,00×10⁻⁹ m²/s. Sau3.600 s, Fo=0,144000. Tính chuỗi cho R=0,571860, còn c_tâm/c₀=0,875185. Tấm đã mất khoảng42,81% lượng ban đầu nhưng tâm mới giảm khoảng12,48%.
Nếu c₀=1,00 mol/m³, nồng độ trung bình và tâm lần lượt0,571860 và0,875185 mol/m³. Không dùng kết quả trung bình từ lượng khí thoát để khẳng định mọi điểm trong tấm đã giảm cùng tỷ lệ.
Sau7.200 s, Fo=0,288000, R=0,398418 và tỷ phần tại tâm0,624890. Sau10.800 s, chúng là0,279176 và0,438490. Các cặp số đều thuộc mô hình hai mặt thoát hoàn hảo, không là dữ liệu chứng minh mức hydro an toàn của một liên kết.
Kiểm xấp xỉ một số hạng
Giữ riêng n=0 cho R₁=(8/π²)exp(−π²Fo/4). Tại một giờ, R₁=0,568176, thấp hơn chuỗi đầy đủ khoảng0,6442% tương đối. Tại hai giờ, R₁=0,398268, gần hơn vì các mode bậc cao đã giảm mạnh.
Không dùng xấp xỉ ấy tại t=0: nó chỉ cho8/π²≈0,810569 thay vì1. Các mode bị bỏ đang mang gần19% lượng trung bình ban đầu. Một công thức đúng ở thời gian đủ lớn có thể sai ngay khi khởi đầu.
Thời gian đặc trưng của mode chậm nhất là τ₀=4L²/(π²D)=10.132,118364 s. Đây là thời gian để biên độ mode ấy giảm theo hệ số exp(−1); không phải thời gian để toàn tấm mất hết hydro hoặc giảm còn một tỷ lệ cố định nào nếu chưa nêu trạng thái ban đầu.
Giải R=0,100 trong chuỗi cho Fo≈0,848085408 và t≈21.202,135201 s, tức5,889482 giờ. Đây chỉ là thời gian đạt trung bình10% trong bài tập. Một mục tiêu về nồng độ tâm sẽ cho thời gian khác.
Chiều dày và mặt thoát tạo hệ số lớn
Vì Fo=Dt/L², giữ cùng tỷ phần R thì thời gian tỷ lệ L²/D. Tăng toàn bộ chiều dày từ10 lên20 mm, giữ D và hai mặt thoát, làm thời gian tăng bốn lần. Tăng D gấp đôi làm thời gian giảm một nửa trong mô hình này.
Nếu tấm vẫn dày10 mm nhưng một mặt bịt hoàn toàn, một mặt giữ c=0, mặt bịt đóng vai trò mặt đối xứng. Chiều dài đặc trưng lúc đó là10 mm thay vì5 mm. Sau một giờ, Fo tương ứng0,036 và R≈0,785905, cao hơn rõ0,571860 của hai mặt thoát.
Đây cũng là lỗi dễ gặp khi người tính nhập chiều dày toàn tấm vào công thức dùng nửa chiều dày. Kết quả thời gian lệch bốn lần dù số liệu D không đổi. Cần vẽ hoặc mô tả hai mặt biên trước khi chọn L.
Khi hệ số thay đổi theo thời gian
Nếu D chỉ phụ thuộc t và đồng đều trong toàn tấm tại mỗi thời điểm, đặt Θ(t)=∫D(t)dt. Quy tắc dây chuyền chuyển phương trình về dạng khuếch tán với thời gian Θ; trong chuỗi thay Dt bằng Θ. Điều này không đúng nguyên dạng nếu D phụ thuộc vị trí, nồng độ hoặc trạng thái bẫy phức tạp.
Giả sử1.800 s đầu D=1×10⁻⁹ m²/s và1.800 s sau D=4×10⁻⁹ m²/s. Khi ấy Θ=9×10⁻⁶ m², Fo=0,360000 và R=0,333474. D trung bình theo thời gian bằng2,5×10⁻⁹ m²/s cho cùng kết quả trong đúng giả thiết này.
Dùng trung bình nhân2×10⁻⁹ m²/s lại cho Fo=0,288 và R=0,398418. Cũng không lấy nhiệt độ trung bình rồi mặc nhiên cho rằng D ở nhiệt độ ấy bằng D trung bình. Quan hệ giữa D và nhiệt độ thường cần mô hình riêng và dữ liệu phù hợp.
Liên hệ với trị số thể tích khí
Nếu riêng một mẫu120 g có trị số ban đầu5,00 ml/100 g trên cùng cơ sở khí chuẩn, tổng lượng biểu diễn bằng thể tích tương đương là6,00 ml. Giả sử toàn bộ phần đang xét tuân mô hình tuyến tính trên, sau một giờ lượng còn tương đương6×0,571860=3,431160 ml.
Đây là phép tỷ lệ lượng chất; trong kim loại không có nghĩa tồn tại một túi khí thể tích3,431160 ml. Cần giữ điều kiện quy đổi khí và mẫu số kim loại. Phần hydro bị giữ không tuân cùng động học không được tự cộng vào c₀ rồi xử lý bằng D này.
Vật liệu thật cần thêm dữ liệu
Nghiên cứu của Rhode và cộng sự trên kim loại hàn P91 cho thấy trạng thái sau hàn và sau nhiệt luyện có hành vi khuếch tán khác nhau; bẫy hydro và hướng vi cấu trúc ảnh hưởng kết quả. [3] Điều này là bằng chứng cần gắn D với vật liệu và trạng thái, không phải cơ sở để dùng số D giả định của bài cho P91.
Để dùng mô hình trong kỹ thuật thực, còn cần phân bố ban đầu, điều kiện bề mặt, lịch sử nhiệt, ứng suất và dữ liệu bẫy. Chuỗi lý tưởng giúp kiểm đơn vị, mức thời gian và sự khác nhau giữa tâm với trung bình. Việc chọn xử lý nhiệt hoặc đánh giá nguy cơ nứt phải dựa trên hồ sơ vật liệu, quy trình và bằng chứng kiểm chứng riêng.
Nguồn tham khảo
[1] How do I measure the diffusible hydrogen level in my ferritic steel weld? — TWI. Truy cập ngày 09/11/2026.
[2] 3.21 Kinetic Processes in Materials — Examination 1, Instructor’s Solution, câu3 — Massachusetts Institute of Technology, 2006. Truy cập ngày 09/11/2026.
[3] Hydrogen diffusion in creep-resistant 9% Cr P91 multi-layer weld metal — Michael Rhode và cộng sự; Welding in the World64,267–281(2020), hồ sơ tác giả tại TUM. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Khi báo cáo hydro trong kim loại hàn, cần biết giá trị nói về phần có thể khuếch tán ra trong điều kiện thử nào, khối lượng kim loại nào và cơ sở quy đổi thể tích khí. TWI phân biệt hydro khuếch tán với phần hydro còn bị giữ tại các vị trí trong kim loại, đồng thời mô tả cách báo thể tích khí quy về điều kiện chuẩn trên100 g kim loại đắp. [1]
Bài này không suy chế độ xử lý nhiệt từ một trị số hydro. Nó xét bài toán truyền khối lý tưởng để hiểu thời gian, chiều dày và điều kiện bề mặt ảnh hưởng lượng còn lại ra sao. Hệ số khuếch tán và mọi kích thước ví dụ đều tự đặt, không đại diện một mác thép hoặc quy trình hàn đã xác nhận.
Mô hình tấm thoát qua hai mặt
Xét tấm phẳng rất rộng, chiều dày2L, tọa độ x tính từ mặt giữa, −L≤x≤L. Ban đầu nồng độ c bằng c₀ đồng đều. Từ t=0, hai mặt được giữ c=0; bỏ trở lực phản ứng bề mặt, bẫy hydro, ứng suất, cạnh tấm và biến đổi thể tích.
Chọn c có đơn vị mol/m³ và hệ số khuếch tán D không đổi, m²/s. Phương trình một chiều là:
∂c/∂t=D∂²c/∂x².
Điều kiện c=0 là một bể thu lý tưởng, không có nghĩa mọi bề mặt thép ngoài không khí luôn thỏa điều kiện ấy. Nếu lớp phủ, phản ứng tái hợp hoặc môi trường cản trở thoát khí, cần điều kiện biên truyền khối khác.
Bài giảng MIT dùng tách biến và chuỗi Fourier cho vật hữu hạn, đồng thời nhấn mạnh điều kiện không thông lượng sẽ bảo toàn tổng lượng chất. [2] Ở đây ta tự áp dụng cùng phương pháp với hai mặt nồng độ bằng0, nên tổng lượng trong tấm giảm theo thời gian.
Tự dựng các mode khuếch tán
Do đối xứng, dùng các hàm cos(λ_nx). Điều kiện tại x=±L yêu cầu λ_n=(2n+1)π/(2L), với n=0,1,2,… Hệ số từ khai triển nồng độ đầu đồng đều là A_n=4(−1)^n/[π(2n+1)].
c/c₀=ΣA_n cos(λ_nx)exp(−Dλ_n²t).
Mỗi mode có tốc độ giảm khác nhau; mode biến thiên nhanh theo không gian mất nhanh hơn. Tích phân theo chiều dày rồi chia2L cho tỷ phần lượng còn lại R=c_trung bình/c₀:
R=(8/π²)Σ exp[−(2n+1)²π²Fo/4]/(2n+1)².
Fo=Dt/L² là thời gian không thứ nguyên. Chuỗi cộng từ n=0 tới vô hạn. Tại t=0, tổng hệ số trung bình bằng1; tại t lớn, mọi số hạng về0. Tại mặt giữa, cos bằng1 nhưng hệ số còn dấu luân phiên, nên nồng độ tâm khác trung bình.
Ví dụ một giờ và hai giờ
Chọn2L=10,0 mm, tức L=0,00500 m, D=1,00×10⁻⁹ m²/s. Sau3.600 s, Fo=0,144000. Tính chuỗi cho R=0,571860, còn c_tâm/c₀=0,875185. Tấm đã mất khoảng42,81% lượng ban đầu nhưng tâm mới giảm khoảng12,48%.
Nếu c₀=1,00 mol/m³, nồng độ trung bình và tâm lần lượt0,571860 và0,875185 mol/m³. Không dùng kết quả trung bình từ lượng khí thoát để khẳng định mọi điểm trong tấm đã giảm cùng tỷ lệ.
Sau7.200 s, Fo=0,288000, R=0,398418 và tỷ phần tại tâm0,624890. Sau10.800 s, chúng là0,279176 và0,438490. Các cặp số đều thuộc mô hình hai mặt thoát hoàn hảo, không là dữ liệu chứng minh mức hydro an toàn của một liên kết.
Kiểm xấp xỉ một số hạng
Giữ riêng n=0 cho R₁=(8/π²)exp(−π²Fo/4). Tại một giờ, R₁=0,568176, thấp hơn chuỗi đầy đủ khoảng0,6442% tương đối. Tại hai giờ, R₁=0,398268, gần hơn vì các mode bậc cao đã giảm mạnh.
Không dùng xấp xỉ ấy tại t=0: nó chỉ cho8/π²≈0,810569 thay vì1. Các mode bị bỏ đang mang gần19% lượng trung bình ban đầu. Một công thức đúng ở thời gian đủ lớn có thể sai ngay khi khởi đầu.
Thời gian đặc trưng của mode chậm nhất là τ₀=4L²/(π²D)=10.132,118364 s. Đây là thời gian để biên độ mode ấy giảm theo hệ số exp(−1); không phải thời gian để toàn tấm mất hết hydro hoặc giảm còn một tỷ lệ cố định nào nếu chưa nêu trạng thái ban đầu.
Giải R=0,100 trong chuỗi cho Fo≈0,848085408 và t≈21.202,135201 s, tức5,889482 giờ. Đây chỉ là thời gian đạt trung bình10% trong bài tập. Một mục tiêu về nồng độ tâm sẽ cho thời gian khác.
Chiều dày và mặt thoát tạo hệ số lớn
Vì Fo=Dt/L², giữ cùng tỷ phần R thì thời gian tỷ lệ L²/D. Tăng toàn bộ chiều dày từ10 lên20 mm, giữ D và hai mặt thoát, làm thời gian tăng bốn lần. Tăng D gấp đôi làm thời gian giảm một nửa trong mô hình này.
Nếu tấm vẫn dày10 mm nhưng một mặt bịt hoàn toàn, một mặt giữ c=0, mặt bịt đóng vai trò mặt đối xứng. Chiều dài đặc trưng lúc đó là10 mm thay vì5 mm. Sau một giờ, Fo tương ứng0,036 và R≈0,785905, cao hơn rõ0,571860 của hai mặt thoát.
Đây cũng là lỗi dễ gặp khi người tính nhập chiều dày toàn tấm vào công thức dùng nửa chiều dày. Kết quả thời gian lệch bốn lần dù số liệu D không đổi. Cần vẽ hoặc mô tả hai mặt biên trước khi chọn L.
Khi hệ số thay đổi theo thời gian
Nếu D chỉ phụ thuộc t và đồng đều trong toàn tấm tại mỗi thời điểm, đặt Θ(t)=∫D(t)dt. Quy tắc dây chuyền chuyển phương trình về dạng khuếch tán với thời gian Θ; trong chuỗi thay Dt bằng Θ. Điều này không đúng nguyên dạng nếu D phụ thuộc vị trí, nồng độ hoặc trạng thái bẫy phức tạp.
Giả sử1.800 s đầu D=1×10⁻⁹ m²/s và1.800 s sau D=4×10⁻⁹ m²/s. Khi ấy Θ=9×10⁻⁶ m², Fo=0,360000 và R=0,333474. D trung bình theo thời gian bằng2,5×10⁻⁹ m²/s cho cùng kết quả trong đúng giả thiết này.
Dùng trung bình nhân2×10⁻⁹ m²/s lại cho Fo=0,288 và R=0,398418. Cũng không lấy nhiệt độ trung bình rồi mặc nhiên cho rằng D ở nhiệt độ ấy bằng D trung bình. Quan hệ giữa D và nhiệt độ thường cần mô hình riêng và dữ liệu phù hợp.
Liên hệ với trị số thể tích khí
Nếu riêng một mẫu120 g có trị số ban đầu5,00 ml/100 g trên cùng cơ sở khí chuẩn, tổng lượng biểu diễn bằng thể tích tương đương là6,00 ml. Giả sử toàn bộ phần đang xét tuân mô hình tuyến tính trên, sau một giờ lượng còn tương đương6×0,571860=3,431160 ml.
Đây là phép tỷ lệ lượng chất; trong kim loại không có nghĩa tồn tại một túi khí thể tích3,431160 ml. Cần giữ điều kiện quy đổi khí và mẫu số kim loại. Phần hydro bị giữ không tuân cùng động học không được tự cộng vào c₀ rồi xử lý bằng D này.
Vật liệu thật cần thêm dữ liệu
Nghiên cứu của Rhode và cộng sự trên kim loại hàn P91 cho thấy trạng thái sau hàn và sau nhiệt luyện có hành vi khuếch tán khác nhau; bẫy hydro và hướng vi cấu trúc ảnh hưởng kết quả. [3] Điều này là bằng chứng cần gắn D với vật liệu và trạng thái, không phải cơ sở để dùng số D giả định của bài cho P91.
Để dùng mô hình trong kỹ thuật thực, còn cần phân bố ban đầu, điều kiện bề mặt, lịch sử nhiệt, ứng suất và dữ liệu bẫy. Chuỗi lý tưởng giúp kiểm đơn vị, mức thời gian và sự khác nhau giữa tâm với trung bình. Việc chọn xử lý nhiệt hoặc đánh giá nguy cơ nứt phải dựa trên hồ sơ vật liệu, quy trình và bằng chứng kiểm chứng riêng.
Nguồn tham khảo
[1] How do I measure the diffusible hydrogen level in my ferritic steel weld? — TWI. Truy cập ngày 09/11/2026.
[2] 3.21 Kinetic Processes in Materials — Examination 1, Instructor’s Solution, câu3 — Massachusetts Institute of Technology, 2006. Truy cập ngày 09/11/2026.
[3] Hydrogen diffusion in creep-resistant 9% Cr P91 multi-layer weld metal — Michael Rhode và cộng sự; Welding in the World64,267–281(2020), hồ sơ tác giả tại TUM. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
