Đếm lượng vật liệu mất trước khi dự báo chiều dày
Một kết quả xói mòn có thể được báo bằng kg/s, kg/(m²·s), mm/năm hoặc kg kim loại trên kg hạt va chạm. Các đại lượng này liên hệ với nhau nhưng không thể thay trực tiếp. Sai đơn vị hoặc nhầm lượng cát đi qua đường ống với lượng thực sự đập vào một vùng thành sẽ tạo sai số trước cả khi xét mô hình vật liệu.
Bài này dựng một vùng thành phẳng nhỏ với diện tích cố định và vật liệu đồng nhất. Chỉ xét mất vật liệu do hạt rắn va chạm, không cộng cơ chế ăn mòn điện hóa. Mọi hệ số dưới đây là giả định để so sánh lịch vận hành; chúng không phải hệ số thiết kế cho thép hoặc thông số của một vị trí thật.
Từ tỷ số xói mòn tới tốc độ mỏng thành
Gọi e là khối lượng vật liệu thành mất trên khối lượng hạt tới vùng xét, đơn vị kg/kg. COMSOL phân biệt các mô hình dùng tỷ số khối lượng này với mô hình dùng thể tích mất trên khối lượng hạt.[2] Vì vậy cần đọc đúng đại lượng đầu ra của mô hình trước khi nhân với mật độ.
Nếu m_dot,p là lưu lượng hạt thực sự tới vùng có diện tích A_w, khối lượng riêng thành ρ_w, thì:
m_dot,loss = e m_dot,p
h_dot = e m_dot,p/(ρ_w A_w)
h_dot có đơn vị m/s khi dùng kg, m và s. Ansys hướng dẫn đổi suất mất khối lượng trên diện tích sang tốc độ giảm chiều dày bằng cách chia cho khối lượng riêng thành.[1] Phép đổi này cho chiều dày tương đương trung bình trên vùng A_w; nó không tự xác định độ sâu lớn nhất của một hố mòn.
Ví dụ, suất mất 2,5×10⁻⁷ kg/(m²·s) trên vật liệu có ρ_w = 7.850 kg/m³ tương ứng 1,004331 mm trong một năm 365 ngày nếu giữ suất ấy liên tục. Con số “mm/năm” là cách chuẩn hóa đơn vị; không bảo đảm quá trình thực sẽ giữ nguyên trong cả năm.
Mô hình vận tốc cần một vận tốc tham chiếu
Đặt riêng mô hình khảo sát:
e = k(v/v_0)ⁿ
k là tỷ số khối lượng tại vận tốc hạt tham chiếu v_0, n không thứ nguyên. Chọn k = 3×10⁻⁷ kg/kg, v_0 = 10 m/s, n = 2,4. Tỷ số v/v_0 tránh việc hệ số thay đơn vị khi đổi m/s sang đơn vị khác. Góc va chạm, hình hạt và kích thước được giả định giữ nguyên trong phạm vi so sánh.
DNV mô tả phép thử xói mòn bằng việc cho lượng và kích thước hạt xác định tác động lên mẫu ở các vận tốc, góc khác nhau.[3] Điều này giải thích dữ liệu cần để hiệu chỉnh mô hình, nhưng không cung cấp giá trị k hoặc n của ví dụ. Đổi vật liệu hoặc dạng hạt phải được đánh giá bằng dữ liệu phù hợp.
Chọn A_w = 0,004 m², m_dot,p = 0,020 kg/s và ρ_w = 7.850 kg/m³. Tại v = v_0, tốc độ mỏng thành bằng 0,000687898 mm/h. Đây là tốc độ trung bình của đúng vùng đang nhận lưu lượng hạt đã cho, không phải tốc độ của toàn bộ bề mặt đường ống.
Tích lũy một lịch vận hành ba chế độ
Xét 1.000 giờ gồm ba đoạn có cùng hình học và cơ chế:
Trong mỗi đoạn, tính h_dot theo dữ liệu của đoạn rồi nhân thời gian tương ứng. Mất chiều dày lần lượt bằng 0,241597; 0,319654 và 0,563908 mm. Tổng bằng 1,125159 mm. Đoạn cuối chỉ chiếm 10% thời gian nhưng đóng góp 50,118% tổng mất chiều dày.
Kiểm cân bằng khối lượng bằng ρ_wA_wΔh cho tổng 0,03532999 kg, tức 35,32999 g vật liệu thành. Khối lượng hạt tới vùng trong cả lịch bằng 79.200 kg. Tỷ số mất vật liệu tích lũy là trung bình có trọng số theo khối lượng hạt, không phải trung bình đơn giản của ba giá trị e.
Cách cộng từng đoạn giả định hình dạng thành chưa đổi đủ để làm thay đổi quỹ đạo hạt, góc va chạm hoặc diện tích tập trung. Nếu hình học mòn làm dòng tái phân bố, phép cộng với k và A_w cố định cần được cập nhật; không cứ kéo dài bảng đến khi chiều dày trở thành âm.
Dùng vận tốc trung bình tạo ra kết quả khác
Trung bình theo thời gian của v bằng 10,2 m/s, còn trung bình m_dot,p bằng 0,022 kg/s. Đưa hai giá trị ấy vào mô hình trong 1.000 h chỉ cho 0,793519 mm, thấp hơn kết quả đúng theo lịch 29,475%.
Có hai nguyên nhân trong ví dụ. Quan hệ vⁿ là phi tuyến; đồng thời giai đoạn vận tốc cao cũng có lượng hạt cao. Phép nhân hai số trung bình đã xóa mối quan hệ đồng thời giữa các biến. Dữ liệu vận tốc và lượng cát cần được gắn cùng mốc thời gian.
Ngay cả khi lượng hạt không đổi, hai nửa thời gian ở 5 và 15 m/s có vận tốc trung bình 10 m/s nhưng mất vật liệu bằng 1,417821 lần trường hợp chạy đều 10 m/s. Phần tăng 41,782% đến từ độ cong của hàm lũy thừa, không phải một sai số đo.
Vận tốc tương đương phải có trọng số đúng
Có thể định nghĩa một v_eq tái tạo đúng tổng mất vật liệu trong chính mô hình này:
(v_eq/v_0)ⁿ = Σ[m_dot,p Δt (v/v_0)ⁿ]/Σ[m_dot,p Δt]
Với lịch ba đoạn, v_eq = 11,797512 m/s. Kết hợp v_eq với tổng khối lượng hạt sẽ khôi phục 1,125159 mm. Đây là một đại lượng rút gọn phục vụ mô hình đã chọn; nó không thay thế phổ vận tốc hạt, phân bố góc hoặc lịch dòng khi đổi sang mô hình khác.
Đặt trọng số theo thời gian chỉ đúng trong trường hợp lượng hạt tới vùng mỗi giây không thay đổi. Nếu hệ thống chỉ lưu khối lượng cát theo ngày nhưng vận tốc theo phút, việc ghép dữ liệu cần nêu rõ cách phân bổ lượng cát trong ngày; không thể khôi phục tương quan đã mất bằng một phép đổi đơn vị.
Giảm vận tốc: phải nói đại lượng nào giữ nguyên
Giữ m_dot,p tới vùng không đổi, giảm v 10% làm h_dot nhân 0,9²˙⁴ = 0,776573, tức giảm 22,343%. Muốn giảm một nửa tốc độ theo cùng giả thiết, v cần nhân 0,5¹⁄²˙⁴ = 0,749154.
Trong một kịch bản khác, giả định nồng độ hạt theo khối lượng trên thể tích dòng không đổi, vận tốc hạt tỷ lệ với vận tốc dòng và tỷ lệ hạt tới vùng không đổi. Khi ấy m_dot,p cũng tỷ lệ v, nên h_dot tỷ lệ v³˙⁴. Giảm v 10% cho hệ số 0,698915; giảm một nửa h_dot chỉ cần hệ số vận tốc 0,815571.
Hai kết quả khác nhau vì điều kiện giữ lượng hạt khác nhau. Không nhân thêm v nếu dữ liệu lưu lượng hạt đã phản ánh thay đổi đó; làm vậy sẽ đếm cùng tác động hai lần. Trong dòng thực, vận tốc hạt và vận tốc trung bình chất lưu cũng chưa chắc bằng nhau.
Độ không đảm bảo và kết luận có điều kiện
Giả định tại một chế độ, độ không đảm bảo tương đối độc lập của k, m_dot,p, v và A_w lần lượt là 15%, 10%, 5% và 8%; bỏ qua độ không đảm bảo ρ_w và n trong ví dụ. Truyền tuyến tính cho:
u_r² = u_k² + u_m² + n²u_v² + u_A²
u_r bằng 23,087%. Công thức chỉ dùng khi các đại lượng và mô hình đã xác định đúng. Nó không bao gồm sai số do dùng vận tốc chất lưu thay vận tốc hạt, do bỏ góc va chạm hoặc do mô hình chưa được hiệu chỉnh. Ansys cũng lưu ý các hệ số mặc định không tự cập nhật theo vật liệu thật.[1]
Nếu chỉ đặt một ngân sách mất chiều dày minh họa 0,40 mm ở chế độ 12 m/s và 0,020 kg/s, mô hình cho thời gian 375,406 h. Đây là phép chia có điều kiện, không phải khoảng kiểm tra hoặc thời gian vận hành được phê duyệt. Để chuyển sang quyết định kỹ thuật, cần dữ liệu chiều dày thực, vị trí sâu nhất, lịch hạt, tương tác cơ chế và tiêu chí đánh giá của thiết bị. Bảng tính nên giữ cả khối lượng mất, diện tích quy đổi và lịch vận hành để người đọc kiểm lại từng bước.
Nguồn tham khảo
[1] Fluent Theory Guide 12.10.6 — Accretion: particle erosion rate and conversion to thickness loss — Ansys — Fluent 2025 R1. Truy cập ngày 09/11/2026.
[2] Erosion Theory — Particle Tracing Module, COMSOL 6.4 — COMSOL. Truy cập ngày 09/11/2026.
[3] Erosion testing — Technology Centre Oslo, Norway — DNV. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
Một kết quả xói mòn có thể được báo bằng kg/s, kg/(m²·s), mm/năm hoặc kg kim loại trên kg hạt va chạm. Các đại lượng này liên hệ với nhau nhưng không thể thay trực tiếp. Sai đơn vị hoặc nhầm lượng cát đi qua đường ống với lượng thực sự đập vào một vùng thành sẽ tạo sai số trước cả khi xét mô hình vật liệu.
Bài này dựng một vùng thành phẳng nhỏ với diện tích cố định và vật liệu đồng nhất. Chỉ xét mất vật liệu do hạt rắn va chạm, không cộng cơ chế ăn mòn điện hóa. Mọi hệ số dưới đây là giả định để so sánh lịch vận hành; chúng không phải hệ số thiết kế cho thép hoặc thông số của một vị trí thật.
Từ tỷ số xói mòn tới tốc độ mỏng thành
Gọi e là khối lượng vật liệu thành mất trên khối lượng hạt tới vùng xét, đơn vị kg/kg. COMSOL phân biệt các mô hình dùng tỷ số khối lượng này với mô hình dùng thể tích mất trên khối lượng hạt.[2] Vì vậy cần đọc đúng đại lượng đầu ra của mô hình trước khi nhân với mật độ.
Nếu m_dot,p là lưu lượng hạt thực sự tới vùng có diện tích A_w, khối lượng riêng thành ρ_w, thì:
m_dot,loss = e m_dot,p
h_dot = e m_dot,p/(ρ_w A_w)
h_dot có đơn vị m/s khi dùng kg, m và s. Ansys hướng dẫn đổi suất mất khối lượng trên diện tích sang tốc độ giảm chiều dày bằng cách chia cho khối lượng riêng thành.[1] Phép đổi này cho chiều dày tương đương trung bình trên vùng A_w; nó không tự xác định độ sâu lớn nhất của một hố mòn.
Ví dụ, suất mất 2,5×10⁻⁷ kg/(m²·s) trên vật liệu có ρ_w = 7.850 kg/m³ tương ứng 1,004331 mm trong một năm 365 ngày nếu giữ suất ấy liên tục. Con số “mm/năm” là cách chuẩn hóa đơn vị; không bảo đảm quá trình thực sẽ giữ nguyên trong cả năm.
Mô hình vận tốc cần một vận tốc tham chiếu
Đặt riêng mô hình khảo sát:
e = k(v/v_0)ⁿ
k là tỷ số khối lượng tại vận tốc hạt tham chiếu v_0, n không thứ nguyên. Chọn k = 3×10⁻⁷ kg/kg, v_0 = 10 m/s, n = 2,4. Tỷ số v/v_0 tránh việc hệ số thay đơn vị khi đổi m/s sang đơn vị khác. Góc va chạm, hình hạt và kích thước được giả định giữ nguyên trong phạm vi so sánh.
DNV mô tả phép thử xói mòn bằng việc cho lượng và kích thước hạt xác định tác động lên mẫu ở các vận tốc, góc khác nhau.[3] Điều này giải thích dữ liệu cần để hiệu chỉnh mô hình, nhưng không cung cấp giá trị k hoặc n của ví dụ. Đổi vật liệu hoặc dạng hạt phải được đánh giá bằng dữ liệu phù hợp.
Chọn A_w = 0,004 m², m_dot,p = 0,020 kg/s và ρ_w = 7.850 kg/m³. Tại v = v_0, tốc độ mỏng thành bằng 0,000687898 mm/h. Đây là tốc độ trung bình của đúng vùng đang nhận lưu lượng hạt đã cho, không phải tốc độ của toàn bộ bề mặt đường ống.
Tích lũy một lịch vận hành ba chế độ
Xét 1.000 giờ gồm ba đoạn có cùng hình học và cơ chế:
- 600 h ở v = 8 m/s, m_dot,p = 0,020 kg/s.
- 300 h ở v = 12 m/s, m_dot,p = 0,020 kg/s.
- 100 h ở v = 18 m/s, m_dot,p = 0,040 kg/s.
Trong mỗi đoạn, tính h_dot theo dữ liệu của đoạn rồi nhân thời gian tương ứng. Mất chiều dày lần lượt bằng 0,241597; 0,319654 và 0,563908 mm. Tổng bằng 1,125159 mm. Đoạn cuối chỉ chiếm 10% thời gian nhưng đóng góp 50,118% tổng mất chiều dày.
Kiểm cân bằng khối lượng bằng ρ_wA_wΔh cho tổng 0,03532999 kg, tức 35,32999 g vật liệu thành. Khối lượng hạt tới vùng trong cả lịch bằng 79.200 kg. Tỷ số mất vật liệu tích lũy là trung bình có trọng số theo khối lượng hạt, không phải trung bình đơn giản của ba giá trị e.
Cách cộng từng đoạn giả định hình dạng thành chưa đổi đủ để làm thay đổi quỹ đạo hạt, góc va chạm hoặc diện tích tập trung. Nếu hình học mòn làm dòng tái phân bố, phép cộng với k và A_w cố định cần được cập nhật; không cứ kéo dài bảng đến khi chiều dày trở thành âm.
Dùng vận tốc trung bình tạo ra kết quả khác
Trung bình theo thời gian của v bằng 10,2 m/s, còn trung bình m_dot,p bằng 0,022 kg/s. Đưa hai giá trị ấy vào mô hình trong 1.000 h chỉ cho 0,793519 mm, thấp hơn kết quả đúng theo lịch 29,475%.
Có hai nguyên nhân trong ví dụ. Quan hệ vⁿ là phi tuyến; đồng thời giai đoạn vận tốc cao cũng có lượng hạt cao. Phép nhân hai số trung bình đã xóa mối quan hệ đồng thời giữa các biến. Dữ liệu vận tốc và lượng cát cần được gắn cùng mốc thời gian.
Ngay cả khi lượng hạt không đổi, hai nửa thời gian ở 5 và 15 m/s có vận tốc trung bình 10 m/s nhưng mất vật liệu bằng 1,417821 lần trường hợp chạy đều 10 m/s. Phần tăng 41,782% đến từ độ cong của hàm lũy thừa, không phải một sai số đo.
Vận tốc tương đương phải có trọng số đúng
Có thể định nghĩa một v_eq tái tạo đúng tổng mất vật liệu trong chính mô hình này:
(v_eq/v_0)ⁿ = Σ[m_dot,p Δt (v/v_0)ⁿ]/Σ[m_dot,p Δt]
Với lịch ba đoạn, v_eq = 11,797512 m/s. Kết hợp v_eq với tổng khối lượng hạt sẽ khôi phục 1,125159 mm. Đây là một đại lượng rút gọn phục vụ mô hình đã chọn; nó không thay thế phổ vận tốc hạt, phân bố góc hoặc lịch dòng khi đổi sang mô hình khác.
Đặt trọng số theo thời gian chỉ đúng trong trường hợp lượng hạt tới vùng mỗi giây không thay đổi. Nếu hệ thống chỉ lưu khối lượng cát theo ngày nhưng vận tốc theo phút, việc ghép dữ liệu cần nêu rõ cách phân bổ lượng cát trong ngày; không thể khôi phục tương quan đã mất bằng một phép đổi đơn vị.
Giảm vận tốc: phải nói đại lượng nào giữ nguyên
Giữ m_dot,p tới vùng không đổi, giảm v 10% làm h_dot nhân 0,9²˙⁴ = 0,776573, tức giảm 22,343%. Muốn giảm một nửa tốc độ theo cùng giả thiết, v cần nhân 0,5¹⁄²˙⁴ = 0,749154.
Trong một kịch bản khác, giả định nồng độ hạt theo khối lượng trên thể tích dòng không đổi, vận tốc hạt tỷ lệ với vận tốc dòng và tỷ lệ hạt tới vùng không đổi. Khi ấy m_dot,p cũng tỷ lệ v, nên h_dot tỷ lệ v³˙⁴. Giảm v 10% cho hệ số 0,698915; giảm một nửa h_dot chỉ cần hệ số vận tốc 0,815571.
Hai kết quả khác nhau vì điều kiện giữ lượng hạt khác nhau. Không nhân thêm v nếu dữ liệu lưu lượng hạt đã phản ánh thay đổi đó; làm vậy sẽ đếm cùng tác động hai lần. Trong dòng thực, vận tốc hạt và vận tốc trung bình chất lưu cũng chưa chắc bằng nhau.
Độ không đảm bảo và kết luận có điều kiện
Giả định tại một chế độ, độ không đảm bảo tương đối độc lập của k, m_dot,p, v và A_w lần lượt là 15%, 10%, 5% và 8%; bỏ qua độ không đảm bảo ρ_w và n trong ví dụ. Truyền tuyến tính cho:
u_r² = u_k² + u_m² + n²u_v² + u_A²
u_r bằng 23,087%. Công thức chỉ dùng khi các đại lượng và mô hình đã xác định đúng. Nó không bao gồm sai số do dùng vận tốc chất lưu thay vận tốc hạt, do bỏ góc va chạm hoặc do mô hình chưa được hiệu chỉnh. Ansys cũng lưu ý các hệ số mặc định không tự cập nhật theo vật liệu thật.[1]
Nếu chỉ đặt một ngân sách mất chiều dày minh họa 0,40 mm ở chế độ 12 m/s và 0,020 kg/s, mô hình cho thời gian 375,406 h. Đây là phép chia có điều kiện, không phải khoảng kiểm tra hoặc thời gian vận hành được phê duyệt. Để chuyển sang quyết định kỹ thuật, cần dữ liệu chiều dày thực, vị trí sâu nhất, lịch hạt, tương tác cơ chế và tiêu chí đánh giá của thiết bị. Bảng tính nên giữ cả khối lượng mất, diện tích quy đổi và lịch vận hành để người đọc kiểm lại từng bước.
Nguồn tham khảo
[1] Fluent Theory Guide 12.10.6 — Accretion: particle erosion rate and conversion to thickness loss — Ansys — Fluent 2025 R1. Truy cập ngày 09/11/2026.
[2] Erosion Theory — Particle Tracing Module, COMSOL 6.4 — COMSOL. Truy cập ngày 09/11/2026.
[3] Erosion testing — Technology Centre Oslo, Norway — DNV. Truy cập ngày 09/11/2026.
Rà soát nội dung: 09/11/2026.
