Mô hình hoá số trị quá trình lan truyền tạp chất vùng ven biển

ĐẠI HỌC QUỐC GIA HÀ NỘI
TRƯỜNG ĐẠI HỌC KHOA HỌC T ự N H IÊ N
TÊN ĐỂ TÀI
M Ô I I Ì M I H O Â S Ố T R Ị
T R Ì N H L A S T I R t r r Ể V T Ạ P C H A I
v r à ; Í E V ' lỉiỂ iv
MÃ s ố : QT-05-05
CHỦ ĐỂ TÀI: TS. TRẨN VÃN c ú c
’ OẠI HỌC QUỎC GlA HÀ NỘI
TRUNC ! ÀM riiỌ_r!ị?TlN THƯ ViỆN
HÀ NÔI - 2006

C5Ĩ- £ □ - BO

ĐẠI HỌC QUỐC GỈA HÀ NỘI
TRƯỜNG ĐẠI HỌC KHOA HỌC Tự NHIÊN
TẺN ĐẺ TẢI
MÔ HÌNH HOÁ SỐ TRỊ
QUÁ TRÌNH LAN TRUYÈN TẠP CHẤT
VÙNG VEN BIẾN
MÃ SỐ: QT-05- 05
CHỦ TRÌ ĐÈ TÀI: TS. TRẦN VĂN cúc
Những người tham gia chính:
TS. TRẦN VĂN TRẢN
NCS. LƯU QUANG HƯNG
CN. HOÀNG ĐÒNG PHƯỢNG
SVCNTN. NGUYẺN THỊ XUẲN
HÀ NỘI - 2006
2
TÓM TẮT ĐẺ TÀI
1- Tên đề tài: Mô hình hoá số trị quá trình lan truyền tạp chất vùng ven biển.

2- Mã số: QT05- 05
3- Chủ trì đè tài: TS.Trần Vãn Cúc
4- Những người tham gia chính: TS. Trần Văn Trản, NCS Lưu Quang Hưnơ, CN.
Hoàng Đồng Phượng, SV(CNTN)Nguyễn thị Xuân.
5- Mục tiêu và nội dung nghiên cứu:
- Tìm hiểu các quá trình lan truyền tạp chất trong nước
- Thiết lập mô hình toán học.
- Tìm lời giải giải tích một số trường hợp đặc biệt
- Xây dựng mô toán học cho bài toán lan truyền tạp chất vùng ven bờ.
- Thiết lập đúng đắn các điều kiện đầu và điều kiện biên cho bài toán.
- Giải sổ bài toán với một nhóm các điều kiện biên và điều kiện đầu .
6- Các kết quả đã đạt được:
- Tìm được lời giải giải tích một số bài toán cụ thể và mở rộng một số lời giải
tích đã có
- Đã xây dựng được mô hình toán học cho bài toán lan truyền tạp chất vùng
ven bờ và đã có lời giải sô bài toán đó.
- Các kết quả nghiên cứu được thể hiện qua một bài báo gửi đãng ở Tạp chí
Cơ học: “The diffusion problem o f compound in weak turbulent region” - “Bài
toán khuếch tán tạp chất trong miền rối yếu”
- M ột báo cáo K hoa học tại Hội nghị Cơ học chất lỏng toàn quốc ‘'Mô hình số
trị hai chiều bài toán khuếch tản tạp chất vừng ven biến” tháng 7- 2006 sap tới
- Hai tiểu luận tốt nghiệp của sinh viên.
3
7- Kinh phí đề tài:
a) Được cấp: 20.000.000đ
b) Chi tiêu:
- Thuê khoán chuyên môn : 12.000.000đ
- Hội nghị,hội thảo khoa học : 3.000.000đ
- Chi khác: 5.000.000đ
Tổng cộng: 20.000.000đ

KHOA QUÁN LÝ
(Ký và ghi rõ họ tên)
CHỦ TRÌ ĐÈ TÀI
(Ký và ghi rõ họ tên)
;
TRƯỜNG ĐẠI HỌC KHOA HỌC T ự NHIÊN
PHÓ HIỆU I KUỎNG
GS.T3.^4«m/ c t l U p / u
A
4
SUMMARY
1- Title of the project: Numerical Modelling the Pollutant Dispersion Process at
Coastal Zones.
2- Code of the project: QT-05- 05
3- Head of the research group: Tran Van Cue, PhD
4- P articipants: Tran Van Tran, PhD; Lưu Quans Hưnơ, PhD Student; Hoang Donơ
Phuong, BSc; Nguyen Thi Xuan, Honors Student;
5- Aims and contents of the project:
- Understanding the pollutant dispersion process in waters
- Reconstructring the mathematical model
- Finding the analytical solution in some special cases
- Establishing the mathematical model for the process of pollutant dispersion in
coastal zones.
- Numerical m odelling for some cases with selected boundary condition and
initial condition.
6- Obtained result:
- Finding analytical solutions for some well-known problems; extending the
result for more some generic case.
- Establishing the mathematical model for the process of pollutant dispersion in
coastal zones; and bluiding programs for these results.

- The result is to be appeared in the Vietnam Journal of Mechanics “The
problem o f pollutant dispersion in weak-turblulent region ”
- A paper in the proceedings of the national conferences on fluid mechanics
“Numerical modelling the pollutant dispersion in areas"
- Two student’s essays
5
7- Fiancials:
a) Receiving:
b) Spending:
- For research works :
- For conferences and seminars :
- For other works:
Total:
20.000.000d
12.000.000d
3.000.000d
5.000.000d
20.000.000d
6
MỤC LỤC
Trang
Mở đầu 9
Chương 1. Hiện tượng khuếch tán tạp chất: Cơ sở lý thuyết 10
1.1. Một số khái niệm 10
1.2. Hiện tượng lan truyền tạp chất 10
1.3. Thiết lập phương trình lan truyền tạp chât 11
1.3.1. Khuếch tán phân tử 11
1.3.2. Khuếch tán rối 12
1.4. Phương trình lan truyền tạp chất một chiều, hai chiều và ba chiều 13
1.4.1. Phương trình lan truyền tạp chất ba chiều 13

1.4.2. Phương trình lan truyền tạp chất hai chiều 13
1.4.3. Phương trình lan truyền tạp chất một chiều 13
Chương 2. M ột số nghiệm giải tích ỉ 4
2.1. Bài toán ba chiều khuếch tán tạp chất của một điểm nguồn 14
2.1.1. Đặt bài toán 14
2.1.2. Phương pháp giải 15
2.1.3. Áp dụng 16
2.1.4. Nhận xét 23
2.2. Bài toán khuếch tán tạp chất trong miền rối yếu 23
2.2.1. Xây dựng mô hình 23
2.2.2. Điều kiện đầu và các điều kiện biên của bài toán 24
2.2.3. M ột số phép biến đổi 25
2.2.4. Trường hợp riêng 28
2.2.5. Kết luận 31
Chương 3. Mô hình số trị hai chiều khuếch tán tạp chất vùng ven biển 32
3.1. Mô hình toán học 32
7
3.1.1. Hệ phương trình 32
3.1.2. Các tham số cùa mô hình 33
3.1.3. Các điều kiện biên và điều kiện ban đầu 33
3.2. Phương pháp số 34
3.2.1. Hệ phương trình dòng chảy 34
3.2.2. Phương trình khuếch tán tạp chất 35
3.3. Kết quả tính toán 36
3.3.1. Bài toán 1 36
3.3.2. Bài toán 2 38
3.4. Nhận xét 40
Kết luận 41
Tài liệu tham khảo 42
Phụ lục 44

8
M Ở ĐẦU
Vùng ven biển nước ta, đặc biệt là những nơi có cửa sông, vịnh, hàng ngày
có một số lượng lớn tạp chất đổ ra biển. Bài toán tìm quy luật phân bố tạp chất cho
những miền đỏ rất có ý nghĩa thực tiễn,vì thế đã có nhiều công trình nghiên cứu bài
toán này .Trong các hướng tiếp cận để giải bài toán, hướng nghiên cứu bàng phương
pháp mô hình hoá số trị cỏ nhiều ưu thế hơn, do dựa trên cơ sở toán học chặt chẽ và
kết quả nghiên cứu có khả năng triển khai ứng dụng. Cho đến nay,một số tác giả đã sử
dụng mô hình số trị sẵn có của nước ngoài [4], hoặc tự xây dựng nhưng sử dụns
phương pháp sai phân hiện đơn eiản [18] trong việc giải bài toán tìm phân bố tạp chất
trons chất lỏng nói chuns và trong vùng ven biển nói riêng.
Khi nghiên cứu bài toán trên, nhiều tác giả cho rằng các mô hình số trị hai
chiều(2D) cho kết đáng tin cậy và vì thế chủng vẫn đang được áp dụng rộng rãi không
nhữna ở nước ta mà cả trên thế giới.Khi sử dụng mô hình 2D,ta có thể tiếp cận bằng
sai phân lưới tam siác [12], bằng phương pháp đặc trưng cải tiến [17], hoặc bằng
phương pháp hệ thông tin địa lý GIS [12 ]
Trong đề tài này, ta sử dụng mô hình hai chiều (2D) đầy đủ, gồm một hệ
phưcmg trình động lực học có tính đến ảnh hưởng của gió, ma sát đáy và nhớt rối, kết
hợp đồng thời với việc giải phương trình khuếch tán tạp chất. Hệ phương trình thuỷ
động lực học trons trường hợp hai chiều được chứng minh là khá tốt cho các vùng
ven biển nước ta [18], Phương trình khuếch tán tạp chất trong mô hình có tính đến
yếu tố phân huỷ hoặc tạo thành tạp chất.
Đề tài đã tự xây dựng phần mềm trên cơ sở một sổ chương trình đã có.
9
Chương I.
HIỆN TƯỢNG KHUÉCH TÁN TẠP CHÁT
CỞ SỞ LÝ THUYÉT
1.1. Các khái niệm
Trước khi đi sâu vào tìm hiểu bản chất toán học của quá trình lan truyền tạp
chất, chúng ta hãy cùng làm quen với một số khái niêm sau đây:

Tải (advection): là quá trình vận chuyển vật chất do dòng chảy.
Đối lint (convection): là quá trình vận chuyển vật chất theo chiều thẳng đứng
sinh ra bởi sự mất ổn định thuỷ tĩnh.
Khuếch tán phân tử (molecular diffusion)', là quá trình sắp xếp vị trí phân tử do
chuyển độns ngẫu nhiên nội tại của các phân tử, được miêu tả bằng định luật Fick.
Khuếch tán rối (turbulent diffusionj: là quá trình sắp xếp ngẫu nhiên vị trí phân
tử vật chất do chuyển động rối.
Trượt (shear)-, là quá trình vận chuyển phần tử vật chất với vận tốc khác nhau
tại các ỉớp vị trí khác nhau.
Phân tán (dispersion)', là quá trình suy giảm mật độ phân tử hoặc đám mây vật
chất dưới tác dụng kết hợp cùa sự trượt và khuếch tán ngang.
Xáo trộn (mixing): là quá trình khuêeh tán hoặc phân tán vật chất; hoặc quá
trình trộn lẫn hay làm loãng các khối chất lỏne này với các khối chất lòng khác.
Bốc hơi (evapolation)'. là quá trình vận chuyển hơi nước từ bề mặt chất lỏng
hoặc chất ran vào không khí.
Lẳng (settling)-, là quá trình chìm xuống của các phần tử vật chất khác nhau về
khối lượng riêng.
Lan truyền (spread)-, là quá trình phân bố theo không gian.
Chất tan (solute)-, là chất phân rã được trong dune môi.
Chất ô nhiễm (pollutant)-, là chất gây hại cho môi trường khi nó được hoà trộn
với không khí, nước, và đất.
Chất bắn (contaminant)', là chất gây hại cho môi trường khi nó được hoà trộn
với không khí, nước, và đất; nó không tự nhiên xuất hiện trong môi trường.
Tạp chất (impurity)', là những chất có mặt trong chất khác, làm cho chất lượng
của chất đó giảm xuống.
1.2. Hiện tượng ian truyền tạp chất
Nếu trong chất lỏng có xuất hiện sự chênh lệch về mật độ vật chất, hoặc có
chuyển độns trona, chẩt lỏna. thì hiện tượng lan truyền xảy ra, Hiện tượng lan truyền
10
có thể do khuếch tán phân tử, khuếch tán rối, hay do tải. Đối với những chất tan, hiện

tượng khuếch tán phân tử liên quan mật thiết với định luật Fick. Trong nghiên cứu
này, chỉ xét đến những tạp chất tan được trong nuớc.
1.3. Thiết lập phưong trình lan truyền tạp chất
1.3.1. Khuếch tán phân tử
Phân tán phân tử của tạp chất (impurity dispersion) liên hệ với mật thiết với hai
quá trình: khuếch tán (diffusion) và tải (ađvection). Ảnh hưởng của khuếch tán và tải
được đặc trưng bởi số Peclet Pe = UL / D : nếu Pe > 1 thì quá trình tải chiếm ưu thế,và
ngược lại khi Pe < 1 quá trình khuếch tán chiếm ưu thế.
Xét một thể tích chất lỏng V giới hạn bởi một mặt kín s . Tạpchất chuyển qua
mặt kín s trong một khoảng thời gian A/ = f2 — được tính qua các quá trình sau:
1- Thay đổi nồng độ tạp chất do khuếch tán phân tử thoả mãn định luật Fick:
trong đó n là véc tơ pháp tuyến ngoài của mặt s, F = («,v,w )là vectơ vận tốc dòne
chảy,C=C(x,y,z,t) là nồng độ tạp chất trong chất lỏng, D là hệ số khuếch tán phân từ.
. 3- Thay đổi nồng độ tạp chất do các quá trình tạo thành tạp chất, hoặc là quá trình
biến đổi nồng độ tạp chất đó do các tác động sinh hoá, vật lý khác:
trong đó F(x,y,z,t) gọi là mật đô. của tạp chất được tạo thành trong một đơn vị thể
tích,trongmột đơn vị thời gian.
Do các quá trình trên,trong thể tích V sẽ diễn ra quá trình biến thiên nồng độ
tạp chất tại mỗi điểm của nó.Trong khoảng thời gian A t,trong V lượng tạp chất biến
thiên sẽ là:
2- Thay đổi nồng độ tạp chất do tài:
/2
ft = J ịịv„cdsdt (1.2)
/I s
a = jjfpw™ (1.3)
Theo định luật bảo toàn khối lượng :
(? 4 -
Q ì ~ ( Q ì + Q ĩ
)
Từ đó ta có:

11
] JJJf'dvá‘= I \\\Fdvdl - ( - ] +) \\V-Cdsd1
|| V U l 11 V V M s uri <1 s
Sử dụng công thức Ostrogradsky- Gauss và biến đổi toán tử g ra d :
j ịịị^ - d V d t = I ịị\F d V d t - - Ị \ị\div(D gradC)dVdt + j ị\ịdìv{vc)d v d t
(I
V dt
(| ị,- _ ,1 ( ,1
V
Do thể tích V và khoảng thời gian Aí được chọn một cách tuỳ ý nên ta thu
được phương trình mô tả quá trình thay đổi nồng độ tạp chất:
Ẽ£ Ẽỵ£. dvvC _ d fDdC'
õt õx õy õy õx V õx J õy
D
5C
õy
+ ■
õz
D õ- ê
õz
+ F (1.6)
1.3.2. Khuếch tản rối
Thực nghiệm cho thấy ràng ảnh hường của khuếch tán phân tử là không đáng
kể so với khuếch tán rối. Do đó trong các lý thuyết bán thực nghiệm, người ta thường
bỏ qua khuếch tán phân tử, chi xét đến khuếch tán rối.
Phương trình khuếch tán trong trường hợp bỏ qua hệ số khuếch tán phân tử và
không tính đến ảnh hường của nguồn có dang:
dC õuC õvC õwC

1

1

1

õt õx õy õy
0 (1.7)
Sử dụng biểu diễn các đặc trưng trong dòng rối theo Reynolds:
u - ũ + u’ ,v = V + v' ,w - W + w’, c = c + C' (1.8)
trong đó các đặc trưng có gạch đầu biểu thị các đặc trưng trung bình còn các đại
lương có dấu phẩy biểu thị các các đặc trưng biến đổi mạch độnơ.
Thay (1.8) vào (1.7) và lấy trung bình theo Reynolds,ta thu được:
õuC ôvC õwC
ÕC

1

Ôt õx
+
õy dy
d u 'C õv'C' õw'C'
õx
ôy
õz
(1-9)
Nếu chất lỏng là đẳng hưởng và không nén được, nhiều tác giả is „ iựng giả
thiết Bousinessq để đưa đến già thiết tỷ lệ grad nồng độ trung bình:
u'C' - -K
ÕC
XX

ôx
v'C' - - K
ÕC
yy
õy
ri
ộ
II
1
ÕC
GZ
12
trong đó K ^ K ^ tK được gọi là các hệ số khuếch tán rối theo các phương x,y,z.
Với các giả thiết đó phương trình (1.9) trở thành:
ÕC ^dC ^ÕC ^Õ C õ

h u —— + V — + w —— = —
dt õx õy õz õx
V
õx
+ -
õy
K.
ÕC
ày
ô_
õz
K.
ÕC
dz

( 1.10)
Trong trường hợp có xét đến nguồn và ảnh hưởng của quá trình láng, sau khi bò qua
ký hiệu trung bình,phương trình (1.10) có dang:
õ ( „ ÕC\ õ ( „ a c ì ÕC\ ôcủC
— K „ — +— K m— +— K = — + — — + F(1.11)
Sxl. d x ) cỊyl õ y ) ô z\ õ z ) õz
ÔC ÕC _ÕC ÕC
-

|_M

hv-

h w
-

dt dx dy õz
trong đó co là vận tốc lắng theo phương thẳng đứng z.
1.4. Phương trình lan truyền tạp chất một chiều, hai chiều và ba chiều
1.4.1. Phương trình lan truyền tạp chất ba chiểu
ÕC ÔC ÕC ÕC
-

V u ———h V ———hw -

ôt õx õy ôz OX
+ -
dy
K ễ £
” ẽy.

8 ( ạ ?
dz
õz
+ — + / ’(1.12)
ôz
1.4.2. Phương trình lan truyền tạp chất hai chiểu
Phương trình lan truyên tạp chất hai chiều có chứa phương thẳng đứna:
ÔC ÔC ÕC õ
(
acì
õ
(
ỔCÌ
+ u-

1- w

-
-
-

K
——
+ —
K :.

dt
õx C1
õx
XX

V
õx , õz
\
õz ,
õcữC . _
+ ^ - + F (1.13)
Õz
Phương trình lan truyền tạp chất hai chiều lấy trung bình theo phương ngang:
+ F (1.14)
ÕC dC ÕC
- — + U- — + V - —
õt dx õy
d_
õx
{ v ỔCÌ
õ ÔC^
Kxx —
+ —
K n.
V dx ,
õ y
>v
V
õ y J
ỉ.4.3. Phương trình lan truyền tạp chẩí một chiểu:
ÕC ÕC õ
õt dx õx
ÕC
õx
+ F (1.15)

Trong mô hình chính tới phương trình lan truyền tạp chất trung bình hai chiều theo
phương ngang (1.14).và phương trình lan íruyền tạp chất ba chiều
13
Chưcmg II.
MỘT SỐ NGHIỆM GIẢI TÍCH
Nghiệm giải tích đóng vai trò quan trọng trong việc bước đẩu đánh giá tính
chính xác của mô hình.
Dưới đây chúng tôi nghiên cứu hai trường hợp có nghiệm giải tích bàne hai
phương pháp khác nhau.
2.1. Bài toán ba chiều khuếch tán tạp chất của một điểm nguồn
Phương trình khuếch tán tạp chất có thể tìm được nghiệm eiải tích trong một số
trường hợp đơn giản.Nghiệm giải tích đóng vai trò quan trọng trons kiểm nghiệm mô
hình số, trong so sánh tính chính xác của thuật toán, cũn® như trong xây đựns các
thiết bị thí nghiệm phục vụ nghiên cứu lan truyền chất tan.
Có nhiều phương pháp đã được đưa ra để tìm nghiệm giải tích của phương
trình lan truyền tạp chất . Tuy nhiên, đa phần các phương pháp này hoặc tươna đối
phức tạp trong tỉnh toán, và đòi hỏi nhiều tài nguyên trong lập trình, hoặc chỉ áp dụng
được cho trường hợp 1- 2 chiều.
Trong phần này ta đề cập đến một phương pháp tìm nehiệm giải tích cho
phương trình khuếch tán 3 chiều. Phương pháp này có ưu điểm là nghiệm thu được
ngắn gọn, hiệu quả lập trình cao, có thể mở rộng cho trường hợp đa nguồn hoặc phân
bổ ban đầu liên tục bất kỳ.
2.1.1. Đặí bài toán
Một mô hình đơn giản cho địa hình được thiết lập dưới dạng nửa không gian,
chiều sâu của đáy là H . Dòng chảy chủ yếu theo hướng dọc theo trục X, vận tốc là
hằng số (u,v,w) = 0 ,0 ,0 ). Trong miền tính, phân bố nồng độ tạp chất là liên tục và giả
thiết các hệ số khuếch tán không đổi.Với các giả thiết đó,phươne trình 3 chiều mô tả
quá trình khuếhc tán tạp chất có dạng:
ÕC ÔC d7C d2C Ổ2C /T,
— + U— = ẢT —f f+/C—f (I)

õt õx “ õx2 w õy2 = ôz2
Tại các biên của nửa khôns gian gồm đáy, bờ và mặt thoáng, phân bố tạp
chất thoả mãn điều kiện không thấm (điều kiện Neumann):
— = 0 tại z = 0,
õy
— = 0 tại z = 0,H
dz
Nồno độ chất tan sẽ giảm bans không khi rất xa nguôn:
14
—— - » 0 k h i X - > ±00
õx
Điều kiện đầu được cho tại một điểm trong không gian qua hàm Dirac:
C(x,y,z,t0) = MS(xs,ys,zs)
2.1.2. Phương pháp giải
Fischer et al (1979) đã sử dụng phương pháp nguồn ào để tìm nghiệm giải
tích của phương trình lan truyền tạp chất một chiều. Phương pháp này dựa trên ý
tưởng: nếu tồn tại một nguồn ảo đối xứng gương qua biên, thì lượng chất tan chuyển
qua biên sẽ cân bàng. Lúc này, gradient nồng độ chất tan sẽ bằng không trên biên.
Ông đã áp dụng nó cho trường hợp với một biên, và trường hợp hai biên với nguồn
đối xứng. Ở đây, xét trường hợp hàm c có dạng bất kỳ:
Giả sử f(x) là nghiệm của bài toán đầu Cauchy. Khi đó hai nguồn ảo đối
xứng gương với f(x) qua biên là / ( 2X0-x) và f(2X0 - X + 2X) tác động đến phân
bố chất tan trong miền tính. Tuy nhiên, tại các biên của nguồn ảo, điều kiện biên
gradient cũng phải thoả mãn, nghĩa là ta thêm vào các nguồn ảo của chúng đại lượng
f(x- 2X) và f{x + 2X). Cứ như vậy chúng có dạng tuần hoàn với các nguồn ảo đan
xen nhau:
fm(x)= f(x ~2nX)
fĩn+Ảx)~ f { x 0 - x + 2(n + i)X)
Nghiệm thoả mãn các điều kiện biên trên là:
c(.v )- lim V {f{x-2nX)+ f(X ữ -x + 2(n + ỉ)x)}

v-*'rr:
Nếu chọn x0 =-L;X = 2Z,;/(x) = exp(-x2) thì ta thu được nghiệm của Fischer
H , \
c(x )= lim y 'e x p ( - ( x - 2 n l)2)
15
2.1.3. Áp dụng
Áp dụng cho bài toán 3 chiều trong kênh. Nếu coi f(x,y,z,t) là hàm mô tả
phân bố trong toàn không gian theo phương X. Nồng độ tạp chất trong miền được cho
bời hệ thức:
c (x ,y ,z,t) = lim lim V Y {/(x,y-2iW,z- 2 jH ,t)+ f{x,~y + 2(i + l)W,z- 2 jH,t)
I
-*5C
J
—>00
Á J
+ f(x,y - 2 iW-z + 20' +1 )H,t)+f(x,-y + 2(i +1 )W ,-z + 2( j +1 )H,t)}
Trở lại phương trình (I), chọn hàm 9[x,y,z,t) sao cho:
u(2x-ut)
c(x,y,z,t) = ớ(x,y,z,t)ex p
Từ đây thu được các biến đổi:
4 K.
(II)
ÔC 09 í \ a ( \
a a exp( )+0txp{ )
ÕC õỡ t \ a í \
^ = ^ exp( ) + ớexPl )
dx dx
( 2 A
u
V 4 KaJ

( \
u
+
6lexp( )
í \2
u
Õ2C _ô2e , , .09 ( >
i ^ =i ? exp( )+2i exp(>
d2c õ2e , ,
ỡ c a2ớ , S
i í =i í exp()
Nhân lần lượt 4 phương trình cuối các hệ số u,Ka,K ,K„ tương ứng rồi cân
bàng theo phương trình lan truyền chất tan, rồi chia cả hai phương trình cho expO, ta
co:
ôt
4 AT
+ u
XX /
õx 2Kr
= K.
d29 „ ôd u
■ + 2
ỏx 2K.
+ 0
í Ỹ
u
qv
/
Sau khi rút gọn :
dỡ

= K
Õ2Ỡ
+ K„
“ ôx2 ' “” ÔV2+K^ õzởi ờx' " ôy
Nghiệm của bài toán Cauchy này với điều kiện đau Dirac sẽlà:
16
0{x,y,z,t) =
M
rexp
K„K„KS
{X-Xs)2 {y-ys)2
4 K j AK t 4 Kjt
(III)
Thay (III) vào (II), ta thu được:
c(x,y,z,t) =
Hay:
c(x,y,z,t) =
M
rexp
K« K„K=
{x-Xsf_ u(2x-ut) ( y - ^ ) 2 {z-zsf
4 K j 4 Kr
4K t 4K j
Sử dụng tính toán ở trẽn ta dẫn đến:
c(x,y,zj) =
M
ẬAm)' K„K„K=
exp
(JC - xs -uíỴ - 2xsuí
AK.J

I J
lim lim 7 V
•i—1 i —J
/=-/
J=-J
exp
(y-2iỉV-ys)2 (z-2jH -zsf
4 K j
4 V
+
exp
exp
exp
(~y + 2W - ys)2 (z -2(j + Ỉ)H - zs)2
4 K j
4 K l
yy
(- ^ + 2(/ + l)^ - X s ) ' {z - 2jH ~zsf
4 K j
+
+
(-y + 2(/ + l)Fr - )2 (z-2 (ý +1)// - z5 y
4Kj
4 V
M
ix ~xs
-ui)2 - 2 xsut
{y-y5)2 (z~ zs)2
V (4«)J K„K„K= " r
4 KJ

4K j 4K J
Kết quả này đúng cho trường hợp hai chiều. Các hình học hai chiều được chúng tôi sử
dụng để minh hoạ cho kết quả tính toán
I OAI HOC QUỐC GIA HÀ MO, '
I ỈPƯNS íAM tHOnG fjNỊ ỉ V-‘FNj
HdmeUCM ỉ 23NƠV/2005 I LQH?0 Solute Dispersion Smitlatron
18
F ra n e o m I 23N Ơ V 2006 I LOH2D -S o lu te D isp e rsio n S m ijaio n
19
20
Fram e001 I 01 Dec 200 5 I LQH 2D - SoltrtBDispersionSfnmadon
f
rame 001 I 01 Dec 2005 I LQH2D Solute Dèsperstcn SfTiulalion
2 0 r
o 10 -
Frame 001 I 2 3 N o r 2 0 05 I LQ H2D - Solute Dispersion Srnulatĩon
z
2.1.4. Nhận xét
Ta nhận thấy ràng hàm phân bố nồng độ tạp chất là hàm e mũ. Phân bố tạp
chất trong miền phụ thuộc vào chiều rộng và độ sâu của miền; vị trí của nguồn trong
miền, cũng như tốc độ dòng chảy. Công thức này cho ta biết nồng độ cùa tạp chất một
cách chính xác tại bất kỳ điểm nào thuộc miền tính.
Neu cố định từng điểm x,y,z, có thể viết gọn công thức trên dưới dạng như
sau:
Ả,
c i( 0 = ^ r f exP
4 K
Rõ ràng giới hạn này tiến đến khôns khi thời gian tiến đến vô cùng. Điều đó
có nghĩa là nồng độ chất tan sẽ giảm đến không sau một khoảng thời gian đù lớn.
Tương tự, nếu như cố định thời gian t và vị trí y,z , công thức trên có thể được viết

lại như sau:
c2 exp
1
K
2 ■
V ^ 5
J
Đây là hàm phân phối chuẩn. Điều đó có nghĩa là dọc theo miền, phân bố tạp
chất vẫn giữ được hình dáng phân phối chuẩn. Bên cạnh đó, giới hạn C2 trona công
thức cũng giảm đến không tại các giá trị đủ lớn. Như vậy càng xa nguồn về phía hai
đầu chiều dài của miền, nồng độ chất tan càng nhỏ.
2.2. Bài toán khuếch tán tạp chất trong miền rối yếu
Trong phần này giới thiệu một mô hình xác định quy luật phân bố tạp chất trong miền
dòng chảy hai chiều có sự xáo trộn rôi yêu .Dòng chày sao cho thành phần vận tôc
theo chiều thẳng đứng biến đổi bé và gradien củathành phần vận tốc theo phưcmg
dòng chảy là rất bé được gọi ỉà miền rối yếu.
2.2.L X ây dựng m ô hình
Phương trình khuếch tán rối tạp chất írone không gian hai chiều,có chứa chiều íhẳne;
đứng theo (1.13) có dạng:
ÔC ÕC
——+U-—
õt Õx
ÕC_
ôz
- í
+
ị (
' k - £ \
Õx '
õx ) õz

, ôz J
ÕŨ)C n ,
+ ^ L + F ạ )
õz
Để đơn giản, ta xét bài toán với một sổ giả thiết sau:
1. Miền chất lỏng có phân bố tạp chất là một hình chữ nhật
£> = {0<x</,0<z</ỉ}
trong đó l,h hoàn toàn xác định.
23
2. Trong miền D phân bổ vận tốc dòng chảy thay đổi yếu theo hiềo sâu hay còn
gọi là miền có xáo trộn rối yếu [1],[3]
Nói cách khác trong miền D ta có thể giả thiết phân bố vận tốc của chất lỏng thoả
. 1 du du
mãn điêu kiện u » w, — » —
dz dx
Với già thiêt đó, ta có thê bỏ qua thành phân vân tôc w và giả thiêt — = £ trong đó £
õx
là hàng số bé khác O.Neu tích phân theo biến X hệ thức đó,ta có:
u(x,z,t) = €x + u0(z,t) (2)
với w0 (z,t) * 0 là phân bố vận tốc tại mặt cắt x=0.
3. Thành phần vận tốc lắng co không đổi.
4. Các thành phần cùa Tenxơ hệ số khuếch tán thoả mãn điều kiện
K:: =K = cosnt và theo [1],[2] có giả thiết
5. K„=au2 (3)
trong đó a là hệ số tỷ lệ, xác định từ những điều kiện cụ thể của bài toán còn
u (x,z,t) xác định theo (1.2)
ó.Trong miền đang xét không có phân bố nguồn.
Tóm lại ta cần xác định quy luật phân bố tạp chất thoả mãn phương trình:
ÔC ÕC ô
—— I- u -— = —

õt õx õx
( r ÕC^
(
õ2c) ÔC
+
K
õz2 ]
+ ú)
—
V õx )
V
õz
(4)
trong miền dòng chảy hai chiều D = {0 < x < /,0 < 2< /j} và u(x,/.,t),Kx thoả mãn các
giả thiết (A2), (A3).
2.2.2. Điều kiện đầu và các điều kiện biên của bài toán
1. Phân bố tạp chất ban đầu trona miền D, giả sử ta đã biết
4-0 =c°(*’z)
2.Tại mặt cắt đầu x=0 và mặt cắt cuối x=l, giả sử xác định được quy luật phán bố của
tạp chất:
cL = ci(z’0 (6)
c\x=l=c:(z,t) (7)
3.Tại đáy và mặt thoáng phân bổ tạp chất thoả mãn điều kiện
cL = c3(*»0
(8)
24
cL = c4(*,0
Hoặc
õc
Õx 2 = 0

= 0
(K^-+coc)
õx
z = h
= 0
(9)
(10)
(11)
trong đỏ các hàm c0, c,, c2, c3, c4 thoả mãn cá điều kiện đòi hỏi của bài toán
2.2.3. M ột số phép biến đổi
Đe giải bài toán đặt ra trên đây, ta thực hiện phép đổi biến về dạna không thứ nauyên
sau đây:
X. = ■
K 2
~ = l n [ l + g(z,t)x],z. ==-,t.
h'Ja n h
£
trong đó
■•=TĨ> (12)
g(z,0 =
u 0( z , t)
(13)
Phép đổi biến (A 12) là đơn trị vì Jacobien cùa phép biến đổi trên có
E ỉnỉn h2 t 0 (M)
D(x ,x,t )
Với phép đổi biến (3.1) phương trình (1.4) trở thành:
õc Ô7C d 2c õc õc , .
— = — ^ + -r^ - - 2 a — + 2b—- (15)
õt. õx. õz. õx. õz.
írong đó:

ị\Jrae)h Củh
2W ' b = ĩẺ
(16)
Miền D trở thành miền D, = {0<x. <l(z.,t.),0< z. < 1 }, với:
/0(z.,/.) = ^ = ln [l + lg(z.,r.)] (17)
shyja
25