Một gradient phải gắn với đúng pha vật liệu
Khi tính sự xuyên qua của một cấu tử qua màng hoặc lớp phủ, biết chiều dày vẫn chưa đủ. Cần phân biệt khả năng di chuyển bên trong vật liệu với mức cấu tử có thể hòa tan vào vật liệu ấy. Hai lớp có cùng hệ số khuếch tán nhưng khác mức hấp thụ ở mặt ngoài có thể truyền lượng chất rất khác nhau. Ghép trực tiếp các giá trị nồng độ đo trong các pha khác nhau còn có thể tạo một gradient không đúng ý nghĩa.
Bài này xét khuếch tán một chiều, đẳng nhiệt, không đối lưu và không phản ứng trong hai lớp phẳng. Mọi hệ số là giả định phục vụ tính toán, không đại diện độ chống thấm của một sản phẩm. Ta bắt đầu ở trạng thái ổn định, sau đó chỉ rõ điều kiện nào khiến cách lấy thông lượng nhân thời gian không còn phù hợp.
Thông lượng khác tổng lưu lượng
Với một cấu tử loãng và D không đổi trong từng lớp, định luật Fick viết:
J = -D dc/dx
J là thông lượng mol, mol/(m²·s); D là hệ số khuếch tán, m²/s; c là nồng độ trong chính pha đang xét, mol/m³; x tính bằng m. Dấu âm làm dòng đi theo chiều giảm c trong mô hình này. COMSOL giải thích D thông qua tỷ số giữa thông lượng và gradient, đồng thời lưu ý tính phụ thuộc môi trường và khả năng cần mô hình nhiều cấu tử.[1]
Tổng lưu lượng qua diện tích A là n_dot = AJ, mol/s. Khi A gấp đôi nhưng các nồng độ biên và cấu tạo lớp giữ nguyên, J không đổi còn n_dot gấp đôi. Không đặt diện tích vào cả công thức J và bước nhân cuối, vì như vậy sẽ đếm diện tích hai lần.
Cân bằng trong một lát vật liệu không phản ứng cho ∂c/∂t = -∂J/∂x. Ở trạng thái ổn định, ∂J/∂x = 0, nên thông lượng xuyên qua hai lớp bằng nhau nếu không tích tụ tại mặt nối. Tài liệu Penn State trình bày mối liên hệ giữa bảo toàn lượng chất và hai định luật Fick.[3] Điều kiện liên tục này không yêu cầu hai nồng độ thuộc hai vật liệu phải bằng nhau.
Đưa các nồng độ về cùng một thang tham chiếu
Đặt một đại lượng tham chiếu φ có đơn vị mol/m³ và giả định cân bằng phân bố tuyến tính c_i = K_iφ. K_i là hệ số không thứ nguyên của lớp i theo cùng pha tham chiếu. Trong bài, φ trùng với nồng độ ở hai bình chứa ngoài khi tiếp giáp đã cân bằng; bên trong mỗi lớp, c_i vẫn là nồng độ vật lý của lớp đó.
Ansys dùng nồng độ chuẩn hóa để xử lý trường hợp các vật liệu có mức bão hòa khác nhau và nồng độ gián đoạn ở mặt tiếp giáp.[2] Mô hình K_i dưới đây là cách đặt giả thiết đơn giản riêng của ví dụ; không đồng nhất trị số K_i với một trường nhập phần mềm nếu định nghĩa chuẩn hóa khác nhau.
Với K_i hằng và không có trở lực riêng của mặt nối:
J = -D_iK_i dφ/dx
R_i = L_i/(D_iK_i)
J = (φ_L-φ_R)/(R_1+R_2)
R_i có đơn vị s/m. Ta cộng được các R_i vì chúng dùng cùng thang φ, cùng diện tích truyền và cùng thông lượng. Nếu lớp bị co diện tích hoặc dòng đi vòng ở mép, phải viết lại tổng lưu lượng theo hình học thực.
Ví dụ hai lớp có độ hòa tan khác nhau
Giả định lớp 1 dày 0,40 mm, D_1 = 2×10⁻¹¹ m²/s, K_1 = 2,0. Lớp 2 dày 0,60 mm, D_2 = 5×10⁻¹² m²/s, K_2 = 0,50. Hai phía ngoài được duy trì φ_L = 10 và φ_R = 0,50 mol/m³. Đổi chiều dày sang m trước khi tính:
R_1 = 1,00×10⁷ s/m
R_2 = 2,40×10⁸ s/m
J = 9,50/(2,50×10⁸)
J = 3,80×10⁻⁸ mol/(m²·s)
Với A = 12 m², lượng truyền ổn định trong một ngày bằng 0,0393984 mol. Nếu chọn khối lượng mol giả định M = 0,032 kg/mol, kết quả tương ứng 1,2607488 g/ngày. Đây là khối lượng của cấu tử đang theo dõi, không phải khối lượng toàn bộ dung dịch hoặc vật liệu.
Nồng độ tham chiếu ở mặt nối là φ_int = 10-JR_1 = 9,62 mol/m³. Ở phía lớp 1, c_1,int = 19,24 mol/m³; ở phía lớp 2, c_2,int = 4,81 mol/m³. Nồng độ hai phía khác nhau bốn lần nhưng φ liên tục theo giả thiết cân bằng phân bố.
Kiểm riêng lớp 1: hai mặt có c bằng 20 và 19,24; D_1×0,76/0,0004 cho đúng J. Lớp 2 có hai mặt 4,81 và 0,25; D_2×4,56/0,0006 cũng cho đúng J. Phép kiểm này đồng thời xác nhận chiều dòng, hệ số phân bố và việc đổi mm sang m.
Bỏ hệ số phân bố tạo sai số có hệ thống
Nếu tự đặt K_1 = K_2 = 1, tổng trở thành 1,40×10⁸ s/m và J thành 6,785714×10⁻⁸ mol/(m²·s), cao hơn 78,57%. Đây không phải sai số làm tròn. Bảng tính đã thay vật liệu bằng hai lớp có tính hấp thụ khác.
Nếu phần mềm hoặc tài liệu cung cấp trực tiếp tích D_iK_i như một hệ số thấm theo đúng thang tham chiếu, không nhân thêm K_i lần nữa. Ngược lại, một hệ số thấm được định nghĩa theo chênh áp suất có đơn vị khác tích đang dùng theo nồng độ. Muốn chuyển đổi phải có quan hệ cân bằng phù hợp, không chỉ đổi tên cột.
Lớp chi phối quyết định độ nhạy
Trong ví dụ, lớp 2 chiếm 96% tổng trở. Giữ các đại lượng khác, tăng gấp đôi chiều dày lớp 1 làm J giảm còn 3,653846×10⁻⁸, tức giảm 3,846%. Tăng gấp đôi riêng lớp 2 làm J giảm còn 1,938776×10⁻⁸, tức giảm 48,98%. Hai phương án cùng gọi là “tăng chiều dày” nhưng tác động khác nhau.
Đạo hàm logarit cho quan hệ cục bộ:
∂lnJ/∂lnD_i = R_i/ΣR
∂lnJ/∂lnL_i = -R_i/ΣR
Vì vậy sai lệch nhỏ của D_2 có ảnh hưởng khoảng 24 lần sai lệch tương đối tương đương của D_1 ở trạng thái gốc. Các đạo hàm dùng cho bước nhỏ; với thay đổi 100%, phải tính lại như hai phép so vừa làm. Không lấy độ nhạy tuyến tính để thay kết quả hữu hạn.
Trở lực ngoài và trạng thái bình chứa
Giả định thêm hai lớp phim ngoài có hệ số truyền theo φ là k_L = 2×10⁻⁶ và k_R = 1×10⁻⁶ m/s. Tổng trở ngoài bằng 1/k_L+1/k_R = 1,50×10⁶ s/m. Khi ấy J giảm còn 3,777336×10⁻⁸, khoảng 0,596% dưới kết quả trước. Trong bộ số này, cải thiện khuấy phía ngoài có tác động nhỏ vì phần lớn trở nằm trong lớp 2.
Một kịch bản khác dùng hai bình kín trộn đều V_L = 0,020 và V_R = 0,010 m³, màng A = 0,50 m². Giả định màng đáp ứng gần ổn định và bỏ lượng tích trữ trong màng. Bảo toàn mol ở hai bình cho chênh nồng độ Δφ giảm theo hàm mũ:
Δφ(t) = Δφ(0)exp(-t/τ)
τ = [A(1/V_L+1/V_R)/ΣR]⁻¹
τ bằng 3.333.333 s, tức 38,5802 ngày. Sau 10 ngày, hai bình có φ_L = 9,276951 và φ_R = 1,946098 mol/m³; J chỉ còn 2,932341×10⁻⁸. Tổng mol vẫn bằng 0,205 mol theo giả thiết bỏ tích trữ màng. Nhân thông lượng ban đầu với 10 ngày sẽ đánh giá cao lượng truyền.
Phạm vi của mô hình cần đi cùng kết quả
Xấp xỉ hai bình vừa dùng cần thời gian cân bằng trong lớp nhỏ so với thời gian biến đổi bình và lượng giữ trong màng đủ nhỏ. Nó không giải giai đoạn nạp chất vào lớp phủ ban đầu. Nếu D phụ thuộc mạnh nồng độ, có phản ứng, nứt xuyên lớp hoặc điện trường, mô hình hệ số hằng cần được thay tương ứng.
Vật liệu xốp còn đòi hỏi phân biệt hệ số của pha lỗ rỗng với hệ số hiệu dụng tính trên toàn tiết diện; Penn State nêu rõ ảnh hưởng đường đi ngoằn ngoèo và cấu trúc lỗ rỗng.[3] Khi bàn giao một kết quả truyền chất, cần giữ định nghĩa c, K, D, diện tích, nồng độ biên và trạng thái thời gian. Chỉ một con số độ thấm không đủ để tái tạo phép tính.
Nguồn tham khảo
[1] Diffusion Coefficient — COMSOL — Multiphysics Cyclopedia. Truy cập ngày 09/11/2026.
[2] 9.2. Normalized Concentration Approach — Ansys — Theory Reference, Release 2025 R2. Truy cập ngày 09/11/2026.
[3] 6.3 Diffusion — CE 574: Reactive Transport in Porous Media — Pennsylvania State University. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Khi tính sự xuyên qua của một cấu tử qua màng hoặc lớp phủ, biết chiều dày vẫn chưa đủ. Cần phân biệt khả năng di chuyển bên trong vật liệu với mức cấu tử có thể hòa tan vào vật liệu ấy. Hai lớp có cùng hệ số khuếch tán nhưng khác mức hấp thụ ở mặt ngoài có thể truyền lượng chất rất khác nhau. Ghép trực tiếp các giá trị nồng độ đo trong các pha khác nhau còn có thể tạo một gradient không đúng ý nghĩa.
Bài này xét khuếch tán một chiều, đẳng nhiệt, không đối lưu và không phản ứng trong hai lớp phẳng. Mọi hệ số là giả định phục vụ tính toán, không đại diện độ chống thấm của một sản phẩm. Ta bắt đầu ở trạng thái ổn định, sau đó chỉ rõ điều kiện nào khiến cách lấy thông lượng nhân thời gian không còn phù hợp.
Thông lượng khác tổng lưu lượng
Với một cấu tử loãng và D không đổi trong từng lớp, định luật Fick viết:
J = -D dc/dx
J là thông lượng mol, mol/(m²·s); D là hệ số khuếch tán, m²/s; c là nồng độ trong chính pha đang xét, mol/m³; x tính bằng m. Dấu âm làm dòng đi theo chiều giảm c trong mô hình này. COMSOL giải thích D thông qua tỷ số giữa thông lượng và gradient, đồng thời lưu ý tính phụ thuộc môi trường và khả năng cần mô hình nhiều cấu tử.[1]
Tổng lưu lượng qua diện tích A là n_dot = AJ, mol/s. Khi A gấp đôi nhưng các nồng độ biên và cấu tạo lớp giữ nguyên, J không đổi còn n_dot gấp đôi. Không đặt diện tích vào cả công thức J và bước nhân cuối, vì như vậy sẽ đếm diện tích hai lần.
Cân bằng trong một lát vật liệu không phản ứng cho ∂c/∂t = -∂J/∂x. Ở trạng thái ổn định, ∂J/∂x = 0, nên thông lượng xuyên qua hai lớp bằng nhau nếu không tích tụ tại mặt nối. Tài liệu Penn State trình bày mối liên hệ giữa bảo toàn lượng chất và hai định luật Fick.[3] Điều kiện liên tục này không yêu cầu hai nồng độ thuộc hai vật liệu phải bằng nhau.
Đưa các nồng độ về cùng một thang tham chiếu
Đặt một đại lượng tham chiếu φ có đơn vị mol/m³ và giả định cân bằng phân bố tuyến tính c_i = K_iφ. K_i là hệ số không thứ nguyên của lớp i theo cùng pha tham chiếu. Trong bài, φ trùng với nồng độ ở hai bình chứa ngoài khi tiếp giáp đã cân bằng; bên trong mỗi lớp, c_i vẫn là nồng độ vật lý của lớp đó.
Ansys dùng nồng độ chuẩn hóa để xử lý trường hợp các vật liệu có mức bão hòa khác nhau và nồng độ gián đoạn ở mặt tiếp giáp.[2] Mô hình K_i dưới đây là cách đặt giả thiết đơn giản riêng của ví dụ; không đồng nhất trị số K_i với một trường nhập phần mềm nếu định nghĩa chuẩn hóa khác nhau.
Với K_i hằng và không có trở lực riêng của mặt nối:
J = -D_iK_i dφ/dx
R_i = L_i/(D_iK_i)
J = (φ_L-φ_R)/(R_1+R_2)
R_i có đơn vị s/m. Ta cộng được các R_i vì chúng dùng cùng thang φ, cùng diện tích truyền và cùng thông lượng. Nếu lớp bị co diện tích hoặc dòng đi vòng ở mép, phải viết lại tổng lưu lượng theo hình học thực.
Ví dụ hai lớp có độ hòa tan khác nhau
Giả định lớp 1 dày 0,40 mm, D_1 = 2×10⁻¹¹ m²/s, K_1 = 2,0. Lớp 2 dày 0,60 mm, D_2 = 5×10⁻¹² m²/s, K_2 = 0,50. Hai phía ngoài được duy trì φ_L = 10 và φ_R = 0,50 mol/m³. Đổi chiều dày sang m trước khi tính:
R_1 = 1,00×10⁷ s/m
R_2 = 2,40×10⁸ s/m
J = 9,50/(2,50×10⁸)
J = 3,80×10⁻⁸ mol/(m²·s)
Với A = 12 m², lượng truyền ổn định trong một ngày bằng 0,0393984 mol. Nếu chọn khối lượng mol giả định M = 0,032 kg/mol, kết quả tương ứng 1,2607488 g/ngày. Đây là khối lượng của cấu tử đang theo dõi, không phải khối lượng toàn bộ dung dịch hoặc vật liệu.
Nồng độ tham chiếu ở mặt nối là φ_int = 10-JR_1 = 9,62 mol/m³. Ở phía lớp 1, c_1,int = 19,24 mol/m³; ở phía lớp 2, c_2,int = 4,81 mol/m³. Nồng độ hai phía khác nhau bốn lần nhưng φ liên tục theo giả thiết cân bằng phân bố.
Kiểm riêng lớp 1: hai mặt có c bằng 20 và 19,24; D_1×0,76/0,0004 cho đúng J. Lớp 2 có hai mặt 4,81 và 0,25; D_2×4,56/0,0006 cũng cho đúng J. Phép kiểm này đồng thời xác nhận chiều dòng, hệ số phân bố và việc đổi mm sang m.
Bỏ hệ số phân bố tạo sai số có hệ thống
Nếu tự đặt K_1 = K_2 = 1, tổng trở thành 1,40×10⁸ s/m và J thành 6,785714×10⁻⁸ mol/(m²·s), cao hơn 78,57%. Đây không phải sai số làm tròn. Bảng tính đã thay vật liệu bằng hai lớp có tính hấp thụ khác.
Nếu phần mềm hoặc tài liệu cung cấp trực tiếp tích D_iK_i như một hệ số thấm theo đúng thang tham chiếu, không nhân thêm K_i lần nữa. Ngược lại, một hệ số thấm được định nghĩa theo chênh áp suất có đơn vị khác tích đang dùng theo nồng độ. Muốn chuyển đổi phải có quan hệ cân bằng phù hợp, không chỉ đổi tên cột.
Lớp chi phối quyết định độ nhạy
Trong ví dụ, lớp 2 chiếm 96% tổng trở. Giữ các đại lượng khác, tăng gấp đôi chiều dày lớp 1 làm J giảm còn 3,653846×10⁻⁸, tức giảm 3,846%. Tăng gấp đôi riêng lớp 2 làm J giảm còn 1,938776×10⁻⁸, tức giảm 48,98%. Hai phương án cùng gọi là “tăng chiều dày” nhưng tác động khác nhau.
Đạo hàm logarit cho quan hệ cục bộ:
∂lnJ/∂lnD_i = R_i/ΣR
∂lnJ/∂lnL_i = -R_i/ΣR
Vì vậy sai lệch nhỏ của D_2 có ảnh hưởng khoảng 24 lần sai lệch tương đối tương đương của D_1 ở trạng thái gốc. Các đạo hàm dùng cho bước nhỏ; với thay đổi 100%, phải tính lại như hai phép so vừa làm. Không lấy độ nhạy tuyến tính để thay kết quả hữu hạn.
Trở lực ngoài và trạng thái bình chứa
Giả định thêm hai lớp phim ngoài có hệ số truyền theo φ là k_L = 2×10⁻⁶ và k_R = 1×10⁻⁶ m/s. Tổng trở ngoài bằng 1/k_L+1/k_R = 1,50×10⁶ s/m. Khi ấy J giảm còn 3,777336×10⁻⁸, khoảng 0,596% dưới kết quả trước. Trong bộ số này, cải thiện khuấy phía ngoài có tác động nhỏ vì phần lớn trở nằm trong lớp 2.
Một kịch bản khác dùng hai bình kín trộn đều V_L = 0,020 và V_R = 0,010 m³, màng A = 0,50 m². Giả định màng đáp ứng gần ổn định và bỏ lượng tích trữ trong màng. Bảo toàn mol ở hai bình cho chênh nồng độ Δφ giảm theo hàm mũ:
Δφ(t) = Δφ(0)exp(-t/τ)
τ = [A(1/V_L+1/V_R)/ΣR]⁻¹
τ bằng 3.333.333 s, tức 38,5802 ngày. Sau 10 ngày, hai bình có φ_L = 9,276951 và φ_R = 1,946098 mol/m³; J chỉ còn 2,932341×10⁻⁸. Tổng mol vẫn bằng 0,205 mol theo giả thiết bỏ tích trữ màng. Nhân thông lượng ban đầu với 10 ngày sẽ đánh giá cao lượng truyền.
Phạm vi của mô hình cần đi cùng kết quả
Xấp xỉ hai bình vừa dùng cần thời gian cân bằng trong lớp nhỏ so với thời gian biến đổi bình và lượng giữ trong màng đủ nhỏ. Nó không giải giai đoạn nạp chất vào lớp phủ ban đầu. Nếu D phụ thuộc mạnh nồng độ, có phản ứng, nứt xuyên lớp hoặc điện trường, mô hình hệ số hằng cần được thay tương ứng.
Vật liệu xốp còn đòi hỏi phân biệt hệ số của pha lỗ rỗng với hệ số hiệu dụng tính trên toàn tiết diện; Penn State nêu rõ ảnh hưởng đường đi ngoằn ngoèo và cấu trúc lỗ rỗng.[3] Khi bàn giao một kết quả truyền chất, cần giữ định nghĩa c, K, D, diện tích, nồng độ biên và trạng thái thời gian. Chỉ một con số độ thấm không đủ để tái tạo phép tính.
Nguồn tham khảo
[1] Diffusion Coefficient — COMSOL — Multiphysics Cyclopedia. Truy cập ngày 09/11/2026.
[2] 9.2. Normalized Concentration Approach — Ansys — Theory Reference, Release 2025 R2. Truy cập ngày 09/11/2026.
[3] 6.3 Diffusion — CE 574: Reactive Transport in Porous Media — Pennsylvania State University. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
