Seepage analysis of a long-polluted river in the steady-state using finite difference method

Abstract

It is very important to analyze the polluted water that streams transfer around them by diffusion. Dirty water leaking from streams like this causes soil and water pollution. The protection of surface and underground water resources is very important given the decreasing water resources. How deep a pollutant has infiltrated can be found by solving the diffusion equation. In this study, the concentration amount of the pollutant passing from a long-polluted stream to the soil was investigated by the finite difference method. Instead of using expensive packet programs, a widespread program Matlab used by engineers and scientists to analyze data, develop algorithms, solve problems analytically and numerically, and create models is preferred in this study. The cross-section of the stream is rectangular and its length is considered to be quite long. Therefore, the problem has been examined in two dimensions. First, the necessary finite difference equations are derived and this problem is analyzed using appropriate boundary conditions. A numerical solution has been obtained for two cases in which the bottom of the stream is permeable and it is not. The mass flow rate of both of the cases was calculated and it is found that the mass flow rate in the first case in which the bottom of the stream is permeable is almost double of the second case. The flow velocities on the boundaries have been found to be less in the case in which the bottom of the stream is impermeable.

Akarsuların etraflarına difüzyon yolu ile aktardıkları kirli suyun analizi oldukça önemlidir. Bu akarsulardan sızan kirli su toprak ve su kirliliğine neden olmaktadır. Yüzey ve yeraltı su kaynaklarının korunması giderek azalan su kaynakları dikkate alındığında çok önemlidir. Bir kirleticinin ne kadar derine sızdığı difüzyon denklemi çözülerek bulunabilir. Bu çalışmada atıksularla kirlenmiş uzun bir dereden toprağa geçen kirleticinin konsantrasyon miktarı sonlu farklar yöntemi ile incelenmiştir. Bu çalışmada pahalı paket programlar kullanmak yerine mühendisler ve bilim insanlarının verileri analiz etmek, algoritma geliştirmek, problemleri analitik ve sayısal olarak çözmek ve model oluşturmak için kullandıkları yaygın bir program olan Matlab tercih edilmiştir. Derenin kesiti diktörtgen alınmış ve boyu oldukça uzun kabul edilmiştir. Bundan dolayı problem iki boyutta incelenmiştir. Önce gerekli sonlu farklar denklemleri türetilmiş ve uygun sınır şartları kullanılarak problem analiz edilmiştir. Akarsuyun tabanının geçirimli olduğu ve olmadığı iki durum için sayısal çözüm elde edilmiştir. Her iki durumun kütlesel debisi hesaplanmış ve akarsuyun tabanının geçirimli olduğu birinci durumdaki kütlesel debinin, ikinci durumun neredeyse iki katı olduğu bulunmuştur. Akarsu tabanının geçirimsiz olması durumunda sınırlardaki akış hızlarının daha az olduğu tespit edilmiştir.