I The workers and entrepreneurs of the recycling chain
2. Mô hình toán học của bể lọc nước
Hầu hết các mô tả toán học của quá trình lọc qua lớp vật liệu lọc hạt được tạo cho trường hợp tốc độ không đổi [5-8]. Điều này do, thực tế rõ ràng là các phương trình cho tốc độ biến đổi phức tạp hơn nhiều so với trường hợp của V = const. Việc thiếu các mô tả toán học đáng tin cậy về lọc với tốc độ thay đổi khiến việc chọn các tham số công nghệ chính của bể lọc - kích thước hạt, chiều cao lớp vật liệu lọc, tốc độ lọc, thời gian chu trình lọc.
Phương trình cơ bản cho quá trình lọc qua một lớp vật liệu lọc hạt được hiển thị bên dưới (một số trong số chúng ở dạng hàm). Phương trình cân bằng cho các hạt lơ lửng trong nước và trong cặn được giữ lại là:
F1(C(x,t), ρ(x,t),m(x,t),V(t)) = 0. (1)
Ở đây C = C(x, t) - nồng độ các chất lơ lửng trong nước; ρ = ρ (x, t) - Nồng độ cặn trong vật liệu lọc; x, t - không gian và thời gian tọa độ; V(t) - tốc độ lọc; m(x, t) - độ xốp vật liệu lọc.
Phương trình động học:
F2(C(x,t), ρ(x,t),m(x,t),V(t),a,b) = 0. (2)
Trong đó b, a là hệ số động học, xác định tương ứng cường độ bám dính và tách rời các hạt lơ lửng khỏi bề mặt tải (hoặc từ các hạt bám dính trước đó).
Độ rỗng vật liệu lọc thay đổi theo độ sâu lớp vật liệu lọc và theo thời gian:
m(х,t) = mo - ρ( )x t,
γ (3)
Trong đó γ là nồng độ khối lượng của các hạt rắn trên một đơn vị thể tích cặn.
Phương trình (1) - (3) xác định động lực của làm trong nước với một lớp vật liệu lọc hạt. Khối thủy động lực của bài
toán được xác định bởi các phương trình: Tổn thất áp lực trong một lớp: hc(t) = ( ) 0 , L I x t dx ∫ (4) Độ dốc thủy lực: I(x,t)=F3(V(t),m(x,t), ν, dx), (5)
Trong đó L - chiều cao của lớp vật liệu lọc, cm; V - tốc độ lọc, cm/s;
ν - hệ số độ nhớt động của nước, cm2/s; đương lượng đường kính hạt,
Khi tăng tổn thất áp lực trong vật liệu lọc thì mực nước H tăng lên trong bể lọc, giảm cung cấp nước vào bể lọc V1 và giảm tốc độ lọc V (xem hình 2).
Phương trình cân bằng nước đến và đi từ bể lọc: V1 - V = dH t( )
dt (6)
Giá trị của V1 được xác định từ điều kiện phụ thuộc phi tuyến tổn thất áp lực từ tốc độ, tức là h = SV2, từ đó
V1=[(Z1-H)/S1]0,5. (7)
Mực nước trong bể lọc bằng – H = Z2 + hс(t) + S2V2. (8)
Các điều kiện ban đầu và điều kiện biên như sau: 0 0 0 0 0 0 ( ), , , , 0, , , ρ ρ α ρ ρ = = = = = = − > = = np x C C t x V V H H m m dc t C C dx (9)
Trong đó ρnp là độ bão hòa tối đa của không gian lỗ rỗng được giữ lại bởi hệ thống lơ lửng.
a - với bộ điều chỉnh tốc độ bộ lọc; b - với dòng chảy gián đoạn; 1 - điều chỉnh công suất; 2 – đồng hồ đo lưu lượng; 3 - bộ điều khiển điều chỉnh; 4 - máng nước tràn
Hình 1. Các sơ đồ bể lọc có tốc độ không đổi:
a. b.
Z1 và Z2 – cốt áp kế trong bộ thu gom nguồn nước đầu vào và thu nước lọc;
H – cốt mực nước trong bể lọc (H0 là cốt mức ở ban đầu lọc);
S1 – sức kháng của ống cung cấp nước cho bể lọc; S2 – sức kháng của ống dẫn nước đã lọc (bao gồm cả hệ thống thu nước ở đáy);
V1 – lưu lượng của nước vào bể lọc trên mỗi đơn vị diện tích của nó,
V - tốc độ lọc
Hình 2. Sơ đồ tính toán bể lọc:
Điều kiện biên cuối cùng có nghĩa là với thời gian lọc dài, lớp này sẽ “kiếm được” và ngừng làm trong nước. Đương nhiên, lúc đầu điều này xảy ra trong các lớp đầu tiên của vật liệu lọc theo đường đi của nước.
Hệ phương trình vi phân, tích phân và đại số (1) - (8) với điều kiện ban đầu và biên (9) là mô hình toán học của bể lọc nhanh với tốc lọc độ biến đổi. Có 8 ẩn số trong hệ phương trình này - C(x,t), ρ(x,t), m(x,t), I(x,t), hc(t), V1(t), V(t), H(t), tương ứng với số phương trình này. Do đó, hệ phương trình được đóng lại và về nguyên tắc, giải được. Tuy nhiên, sự phức tạp của hệ phương trình này không có lý do để hy vọng để có được các giải pháp phân tích. Cách tự nhiên là sử dụng phương pháp số.
Phương trình (1) - (2) được trình bày dưới dạng sai phân hữu hạn, với kết quả sau khi một số đơn giản hóa rõ ràng, sử dụng phương trình động học theo [5], chúng ta có: C V t x ρ ∆ = − ∆ ∆ ∆ (1’) C bC a x V ρ ∆ = − + ∆ (2’)
Từ (1’) và (2’) chúng ta có được các quan hệ lặp lại để tính toán ρ và C: ρij+1 = ρij – Vj+1 t(Ci 1– ij), x − C ∆ ∆ (10) Ci+1,j = Cij + ∆x(-bCij+ ij) j a V ρ , (11)
Trong đó Δx và Δt - bước đếm ở độ sâu của lớp và theo thời gian; i - số lớp theo chiều sâu;
j - số của lớp “tạm thời”.
Tổn thất áp lực trong lớp, dựa trên phương trình (5), bằng hc = Δx
0 ( , )
L
I x t
∑ . (12)
Sự thay đổi mực nước trong bể lọc được tính theo tỷ số - V1 – V = H t( ) t ∆ ∆ Từ đó Hj+1 = Hj + Δt (V1j –Vj). (13)
Cấp nước cho bể lọc và tốc độ lọc được xác định bởi công thức, tiếp theo từ (7) đến (8):
V1j = [(Z1 – Hj)/S1]0,5; (14)
Vj = (Hj –Z2–hcj)/S2 . (15)
Tổn thất áp suất trong lớp cơ bản của độ dày Δx bằng:
hij = ΔxF3(Vj,mij). (16)
Lập thuật toán để tính toán mô hình kết quả, thử nghiệm thành công về vấn đề lọc với tốc độ không đổi [9]:
1. Cho dữ liệu ban đầu – C0, Cф, Lo, Δx, Δt, ρo(x), dэ, tb, a, b, mo, γ, Z1, Z2, H0.
2. Tính V10 và V0.
3. Lấy t = 0 và tạo ra một phép tính từng lớp của ρ(x,0), C(x,0), m(x,0), h(x,0) đến giá trị x=L.
4. Tìm tổng tổn thất áp suất trong vật liệu lọc của bể lọc hс(0).
5. Chuyển đến lần tiếp theo tj + 1 = tj + Δt, trong đó các giá trị mới của H, V1, V và Vп được xác định tuần tự.
6. Khi giá trị tìm thấy của tốc độ lọc V thực hiện tính toán lớp ρ, C, m, h.
7. Tính toán được lặp lại cho đến khi đạt được giá trị được chỉ định khả năng lọc (thường không quá 48 giờ).
3. Kết luận
Thuật toán được mô tả đã được triển khai bằng Microsoft Excel sử dụng macro. Một số tính toán đã được thực hiện cho thấy sự tương ứng đầy đủ của kết quả của chúng với các khái niệm vật lý của các quá trình xảy ra trong quá trình lọc với tốc độ giảm. Ví dụ, hình 3 cho thấy động lực của những thay đổi trong tốc độ lọc, tổn thất áp suất trong vật liệu lọc và mực nước trong bể lọc. Ở đây có thể thấy rằng khi bắt đầu lọc mực nước trong bể lọc và tổn thất áp suất trong vật liệu lọc nhanh chóng tăng, và tốc độ tăng trưởng của chúng sẽ chậm lại. Tốc độ lọc cũng tăng nhanh, và sau đó, đạt mức tối đa, giảm dần. Điều này được giải thích bởi thực tế là ở cốt đầu mức trong bể lọc thấp (nó tương ứng với cốt của cạnh của máng thu nước rửa lọc). Do đó, tỷ lệ cấp nước cho bể lọc lớn (sự khác biệt Z1 - Ho là lớn), tốc độ lọc thấp, kết quả là bể lọc được nạp đầy nhanh chóng. Vì vậy, chúng tôi đã có được một mô hình toán học của một bể lọc nhanh, với tốc độ giảm, và một thuật toán số được phát triển thực hiện mô hình. Tính toán đã được thực hiện cho thấy sự tuân thủ đầy đủ kết quả của chúng với khái niệm vật lý về các quy trình lọc tốc độ lọc giảm.
Các nhiệm vụ chính của nghiên cứu sâu hơn là:
- Nghiên cứu ảnh hưởng của các tham số công nghệ (tốc độ ban đầu chiều cao lớp, kích thước hạt, v.v.), cũng như các đặc tính động học ảnh hưởng đến biến động lọc, tốc độ biến thiên;
- Phê duyệt thử nghiệm mô hình toán học thu được./.
T¿i lièu tham khÀo
1. Клячко В.А., Апельцин И.Э. Очистка природных вод. – М.: Стройиздат, 1971. – 579 с.
2. Hudson H.E. Declining rate filtration // JAWWA. Vol.51, №11,1959. – p. 42-50.
3. Cleasby J.L. Water filtration through deep granular media // Public Works, №6, 1970.– p.36-5. 4. Сысоев М.Н., Казакова Л.П., Богданова С.И., Круглов Л.С. Работа фильтрующих сооружений с переменной скоростью // Водоснабжение и санитарная техника. –1968.– №2. – С.15-19. 5. Минц Д. М.Теоретические основы технологии очистки воды.– М.:Стройиз-дат,1964.–56с. 6. Венецианов Е. В., Рубинштейн Р. Н. Динамика сорбции из жидких сред.–М.:Наука, 1983. – 237 с. 7. Олейник А. Я., Тугай А. М. Моделирование процессов кольматажа и суффозии вприфильтровой зоне скважины // Докл. НАН Украины. – 2001. – №9. – С.190 – 194. 8. Поляков В. Л. О фильтровании суспензий при заданном напоре // Докл. НАН Украины. – 2005. – №4. – С.48-54. 9. Грабовський П. О., Гурінчик Н.О. Чисельна реалізація математичної моделі фі- льтрування // Проблеми водопостачання, водовідведення та гідравліки: Наук.-техн. зб.Вип.6. – К., 2005. – С.4-13.