Submit manuscript...
MOJ
eISSN: 2573-2919

Ecology & Environmental Sciences

Research Article Volume 4 Issue 2

Contribution to the wastewater by use of a numerical study of a secondary clarifier comsol multiphysics under

Fares Redouane,1 Aissa Abdelrahmane,2 Bouadi Abed,3 Lounis Mourad1

1Laar, University of Science and Technology, Algeria
2(LPQ3M), University Mustapha Stambouli Mascara, Algeria
3Universitére Ahmed Zabana RELIZANE Center, Algeria

Correspondence: Dr. Fares Redouane, Ahmed ZABANA University Center of Relizane, Laar, University of Science and Technology, Mohamed Boudiaf, BP 1505, Oran, Algeria, Tel 00213771728504

Received: December 10, 2018 | Published: March 13, 2019

Citation: Redouane F, Abdelrahmane A, Abed B, et al. Contribution to the wastewater by use of a numerical study of a secondary clarifier comsol multiphysics under. MOJ Eco Environ Sci. 2019;4(2):44-47. DOI: 10.15406/mojes.2019.04.00131

Download PDF

Abstract

The environmental assessment of industrial effluents (containing organic pollutants and minerals without proper treatment), generated by the mineral fertilizer factory and emptied into the seacoast, revealed the concentrations intense of suspended solids (SS).The wastewater treatment is a multi-step process to remove contaminants. At first, large solid particles were removed by sedimentation, flotation and filtration. In a second step, the biological treatment is aggregated by forming the smaller particles, which is called flocs. These flocs can be easily removed by sedimentation. The present investigation studying the separation of water flocs in a circular secondary clarifier. A mixture model is used to model the multiphase turbulent flow in the tank.

Keywords: comsol, simulation, pollution, clarifier

Introduction

Our modern life imposes increasingly important requirements on what to eat, dress and how to manage time and space by finding daily means and methods that are increasingly sophisticated to meet the maximum comfort.1,2 Any economic activity caused by the human needs in the entire planet, generates more technological development and human living conditions, pollution problems accompany most of the industrial and agricultural processes including environmental and ecological consequences are very serious and can in some situations simply turn against the objectives initially set by these companies.3,4 The degradation of the Earth's environment caused by the intense industrial activities, especially in some developed countries thus enter into the concerns of the nations. Algeria is also affected by pollution in general, and its coast is relatively polluted as it was reported in several research5,6 than untreated industrial effluents discharged directly into the Mediterranean sea are the primary causes.7 They also showed a growing pollution in the entire Algerian coast, especially in the coastal areas of large cities. Primary causes.7,8 They also showed that a growing pollution in the entire Algerian coast, especially in the coastal areas of large cities.

We have tried our best to provide solutions and respond to pollution issues raised during the experimental measurements of physico-chemical parameters features exploiting the numerical simulation tool.5 These answers are powerful and efficient, and enable the Arzew Fertial-run business in clean environmental standards. To contributed to reducing pollution after Fertial-Arzew complex, proposed a mathematical model of turbulence kε. To follow the mixing speed, pressure and the phase volume fractions, we propose the application of the mixture model. The latter is a multi phasic flow pattern, particularly suitable for suspensions or mixtures of solid or liquid particles.

Model definition

By analyzing all released liquids of Fertial-Arzew company directly at its channels evacuations,9 it appears to be the the overruns legislative standards. For reasons of efficiency, it is attempted to provide a solution for the parameter measuring solids MY suspension, whose values ​​measured at the three liquid effluent discharge channels exceed a very exaggerated manner standards.

The values ​​from the MES parameter of treated wastewater were significantly lower than those discharged directly into the sea. The wastewater treatment is a process that involves several steps to remove contaminants. The most important of these steps is one in which large solid particles are removed by sedimentation, flotation and filtration obtaining flocs or sludge on the one hand and lightly loaded waste water solids suspension on the other hand. The present investigation studying the separation of water in a circular flocs settling tank used in sewage treatment systems. Figure 1 shows the representative diagram in longitudinal section of the structure of a circular decanter.

Figure 1 Longitudinal section diagram showing 2/3 a circular clarifier.

Liquid wastes from Fertial-Arzew company arrive via a channel flowing at the entrance of the clarifier at its center. The treated water leaves at its ends leaving the clarifier decant the bottom solids suspensions. In recent exit through the bottom of clarifier under the effect of gravity to evacuate to dedicated sites or also for a suitable adjustment. To accentuate the wall settling phenomenon ("S") is provided in the center of the clarifier to prevent the waste water jet arriving through inlet re-mix sludge already settled to the bottom of clarifier.

The clarifier in question has an axial symmetry which allows us to reduce the number calculations operations by merely half the system in 2D, as shown in Figure 2. The first figure shows the axisymmetric diagram of half of the clarifier 2D and the second mesh in its triangular cells.

Figure 2 Representation of the diagram of axisymmetric geometry and its triangular mesh.

Our model is simulated on a circular secondary clarifier,10 the flocs are removed from water by sedimentation and this through gravity. The purpose of this study is to investigate the complex turbulent poly phasic flow in a circular secondary clarifier.

FIG 1montre geometry of the clarifier

The sludge used, consisting of a mixture of solid flakes and water enters through the inlet of the reservoir. The fuel tank has two outputs. One is located in the center of the tank bottom. The purpose of the outlet port is to remove the flocs sedimented from the reservoir. There is also a peripheral outlet purified water as shown in Figure 2, which also shows the corresponding model in 2D axisymmetric. We modeled the floc by solid circular of equal size. The model consists of a set of equations, boundary conditions and initial conditions. The set is defined on a geometry decomposed into sub-areas and limited by borders. This model is simulated with COMSOL Multiphysics, or we first made the physical description of the problems (variables, equations, initial limits and conditions), followed by the implementation of our model, and finally made discussions on the results.

Mathematical model

Based on the conservation of mass and momentum of each phase, the mixture model uses the following equations:

ρ u t +ρ( u. )u=p.( ρ c d ( 1 c d ) u slip 2 )+. τ Gm +ρg MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaeqyWdiNaamyDa8aadaWgaaWcbaWdbiaadshaa8aabeaak8qacqGH RaWkcqaHbpGCdaqadaWdaeaapeGaamyDaiaac6cacqGHhis0aiaawI cacaGLPaaacaWG1bGaeyypa0JaeyOeI0Iaey4bIeTaamiCaiabgkHi TiabgEGirlaac6cadaqadaWdaeaapeGaeqyWdiNaam4ya8aadaWgaa WcbaWdbiaadsgaa8aabeaak8qadaqadaWdaeaapeGaaGymaiabgkHi TiaadogapaWaaSbaaSqaa8qacaWGKbaapaqabaaak8qacaGLOaGaay zkaaGaamyDa8aadaqhaaWcbaWdbiaadohacaWGSbGaamyAaiaadcha a8aabaWdbiaaikdaaaaakiaawIcacaGLPaaacqGHRaWkcqGHhis0ca GGUaGaeqiXdq3damaaBaaaleaapeGaam4raiaad2gaa8aabeaak8qa cqGHRaWkcqaHbpGCcaWGNbaaaa@6853@     (1)

( ρ c ρ d )[ .( ϕ d ( 1 c d ) u slip D md ϕ d ) ]+ ρ c ( .u )=0 MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape WaaeWaa8aabaWdbiabeg8aY9aadaWgaaWcbaWdbiaadogaa8aabeaa k8qacqGHsislcqaHbpGCpaWaaSbaaSqaa8qacaWGKbaapaqabaaak8 qacaGLOaGaayzkaaWaamWaa8aabaWdbiabgEGirlaac6cadaqadaWd aeaapeGaeqy1dy2damaaBaaaleaapeGaamizaaWdaeqaaOWdbmaabm aapaqaa8qacaaIXaGaeyOeI0Iaam4ya8aadaWgaaWcbaWdbiaadsga a8aabeaaaOWdbiaawIcacaGLPaaacaWG1bWdamaaBaaaleaapeGaam 4CaiaadYgacaWGPbGaamiCaaWdaeqaaOWdbiabgkHiTiaadseapaWa aSbaaSqaa8qacaWGTbGaamizaaWdaeqaaOWdbiabgEGirlabew9aM9 aadaWgaaWcbaWdbiaadsgaa8aabeaaaOWdbiaawIcacaGLPaaaaiaa wUfacaGLDbaacqGHRaWkcqaHbpGCpaWaaSbaaSqaa8qacaWGJbaapa qabaGcpeWaaeWaa8aabaWdbiabgEGirlaac6cacaWG1baacaGLOaGa ayzkaaGaeyypa0JaaGimaaaa@6845@     (2)

t ( ϕ d ρ d )+. ϕ d ρ d u d =0 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape WaaSaaa8aabaWdbiabgkGi2cWdaeaapeGaeyOaIyRaamiDaaaadaqa daWdaeaapeGaeqy1dy2damaaBaaaleaapeGaamizaaWdaeqaaOWdbi abeg8aY9aadaWgaaWcbaWdbiaadsgaa8aabeaaaOWdbiaawIcacaGL PaaacqGHRaWkcqGHhis0caGGUaGaeqy1dy2damaaBaaaleaapeGaam izaaWdaeqaaOWdbiabeg8aY9aadaWgaaWcbaWdbiaadsgaa8aabeaa k8qacaWG1bWdamaaBaaaleaapeGaamizaaWdaeqaaOWdbiabg2da9i aaicdaaaa@51BB@     (3)

Where, u represents the velocity of the mixture (m/s), ρ is the mixture density (kg/m3), p is the pressure (Pa), and Cd is the mass fraction of the solid phase (kg/kg). In addition, uslip is the relative velocity between the two phases (m/s), is the sum of the viscous and turbulent stress (kg/(m·s2)), and g is the gravity (m/s2). The mixing speed (m/s) is defined as τ Gm MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaeqiXdq3damaaBaaaleaapeGaam4raiaad2gaa8aabeaaaaa@3C28@

u= ϕ c ρ c u c + ϕ d ρ d u d ρ MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaamyDaiabg2da9maalaaapaqaa8qacqaHvpGzpaWaaSbaaSqaa8qa caWGJbaapaqabaGcpeGaeqyWdi3damaaBaaaleaapeGaam4yaaWdae qaaOWdbiaadwhapaWaaSbaaSqaa8qacaWGJbaapaqabaGcpeGaey4k aSIaeqy1dy2damaaBaaaleaapeGaamizaaWdaeqaaOWdbiabeg8aY9 aadaWgaaWcbaWdbiaadsgaa8aabeaak8qacaWG1bWdamaaBaaaleaa peGaamizaaWdaeqaaaGcbaWdbiabeg8aYbaaaaa@4E4B@     (4)

Ouet represent respectively the volume fractions of the liquid phase (continuous phase) and the solid phase (discontinuous phase) (m3/m3), the velocity of the liquid phase (m/s) the speed of the solid phase (m/s), the density ρc liquid phase (kg/m3), the density ρd solid phase (kg/m3), ρ is the density and mixture of (kg/m3). The relationship between the two phases of speeds is defined by ϕ c   ϕ d u c u d MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape Gaeqy1dy2damaaBaaaleaapeGaam4yaaWdaeqaaOWdbiaabckacqaH vpGzpaWaaSbaaSqaa8qacaWGKbaapaqabaGcpeGaamyDa8aadaWgaa WcbaWdbiaadogaa8aabeaak8qacaWG1bWdamaaBaaaleaapeGaamiz aaWdaeqaaaaa@444A@

u cd = u d u c = u slip D md ( 1 c d ) ϕ d ϕ d MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaamyDa8aadaWgaaWcbaWdbiaadogacaWGKbaapaqabaGcpeGaeyyp a0JaamyDa8aadaWgaaWcbaWdbiaadsgaa8aabeaak8qacqGHsislca WG1bWdamaaBaaaleaapeGaam4yaaWdaeqaaOWdbiabg2da9iaadwha paWaaSbaaSqaa8qacaWGZbGaamiBaiaadMgacaWGWbaapaqabaGcpe GaeyOeI0YaaSaaa8aabaWdbiaadseapaWaaSbaaSqaa8qacaWGTbGa amizaaWdaeqaaaGcbaWdbmaabmaapaqaa8qacaaIXaGaeyOeI0Iaam 4ya8aadaWgaaWcbaWdbiaadsgaa8aabeaaaOWdbiaawIcacaGLPaaa cqaHvpGzpaWaaSbaaSqaa8qacaWGKbaapaqabaaaaOWdbiabgEGirl abew9aM9aadaWgaaWcbaWdbiaadsgaa8aabeaaaaa@59E9@     (5)

For the solid particles sliding speed, we can use the drag act Hadamard-Rybczynski. In the mode of application of the mixture model of the Hadamard-Rybczynski slip speed is provided as a predefined template. The mode of application then calculates the slip speed as a result:

u slip = ( ρ ρ d ) d d 2 18ρη p MathType@MTEF@5@5@+= feaagKart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaamyDa8aadaWgaaWcbaWdbiaadohacaWGSbGaamyAaiaadchaa8aa beaak8qacqGH9aqpcqGHsisldaWcaaWdaeaapeWaaeWaa8aabaWdbi abeg8aYjabgkHiTiabeg8aY9aadaWgaaWcbaWdbiaadsgaa8aabeaa aOWdbiaawIcacaGLPaaacaWGKbWdamaaDaaaleaapeGaamizaaWdae aapeGaaGOmaaaaaOWdaeaapeGaaGymaiaaiIdacqaHbpGCcqaH3oaA aaGaey4bIeTaamiCaaaa@519E@     (6)

Where denotes the diameter of the solid particles (m). To the mixture density and viscosity, we offer the following expressions dd

η= η c ( 1 ϕ d ϕ max ) 2.5 ϕ max MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape Gaeq4TdGMaeyypa0Jaeq4TdG2damaaBaaaleaapeGaam4yaaWdaeqa aOWdbmaabmaapaqaa8qacaaIXaGaeyOeI0YaaSaaa8aabaWdbiabew 9aM9aadaWgaaWcbaWdbiaadsgaa8aabeaaaOqaa8qacqaHvpGzpaWa aSbaaSqaa8qacaWGTbGaamyyaiaadIhaa8aabeaaaaaak8qacaGLOa GaayzkaaWdamaaCaaaleqabaWdbiabgkHiTiaaikdacaGGUaGaaGyn aiabew9aM9aadaWgaaadbaWdbiaad2gacaWGHbGaamiEaaWdaeqaaa aaaaa@5218@     (7)

ρ= ϕ c ρ c + ϕ d ρ d MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaeqyWdiNaeyypa0Jaeqy1dy2damaaBaaaleaapeGaam4yaaWdaeqa aOWdbiabeg8aY9aadaWgaaWcbaWdbiaadogaa8aabeaak8qacqGHRa WkcqaHvpGzpaWaaSbaaSqaa8qacaWGKbaapaqabaGcpeGaeqyWdi3d amaaBaaaleaapeGaamizaaWdaeqaaaaa@485B@     (8)

The turbulence model k-ε is used to determine the eddy viscosity. The mode of application calculates the dispersion of particles of Dmd coefficient (m2/s) from

D md = η T ρ σ T MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape Gaamira8aadaWgaaWcbaWdbiaad2gacaWGKbaapaqabaGcpeGaeyyp a0ZaaSaaa8aabaWdbiabeE7aO9aadaWgaaWcbaWdbiaadsfaa8aabe aaaOqaa8qacqaHbpGCcqaHdpWCpaWaaSbaaSqaa8qacaWGubaapaqa baaaaaaa@4447@     (9)

Where ηT is the eddy viscosity (Pa.s) and σT is the (dimensionless) turbulent Schmidt number,

The mass flow of the solid phase through the inlet is given by

M in = inlet | v in | ϕ din ρ d dS MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape Gaamyta8aadaWgaaWcbaWdbiaadMgacaWGUbaapaqabaGcpeGaeyyp a0ZaaybCaeqal8aabaWdbiaadMgacaWGUbGaamiBaiaadwgacaWG0b aapaqaaaqdbaWdbiabgUIiYdaakmaaemaapaqaa8qacaWG2bWdamaa BaaaleaapeGaamyAaiaad6gaa8aabeaaaOWdbiaawEa7caGLiWoacq aHvpGzpaWaaSbaaSqaa8qacaWGKbGaamyAaiaad6gaa8aabeaak8qa cqaHbpGCpaWaaSbaaSqaa8qacaWGKbaapaqabaGcpeGaamizaiaado faaaa@545B@     (10)

Correspondingly to the flow through the peripheral outlet

M out = outlet | v out | ϕ dout ρ d dS MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape Gaamyta8aadaWgaaWcbaWdbiaad+gacaWG1bGaamiDaaWdaeqaaOWd biabg2da9maawahabeWcpaqaa8qacaWGVbGaamyDaiaadshacaWGSb Gaamyzaiaadshaa8aabaaaneaapeGaey4kIipaaOWaaqWaa8aabaWd biaadAhapaWaaSbaaSqaa8qacaWGVbGaamyDaiaadshaa8aabeaaaO WdbiaawEa7caGLiWoacqaHvpGzpaWaaSbaaSqaa8qacaWGKbGaam4B aiaadwhacaWG0baapaqabaGcpeGaeqyWdi3damaaBaaaleaapeGaam izaaWdaeqaaOWdbiaadsgacaWGtbaaaa@5873@     (11)

The difference between the latter two entities provides the remaining mass at the bottom clarifier. Note that the mode of application of the model mix automatically configures all the modeling equations.

Convergence of the solution

The numerical method is to interpolate the solution on a discretization of the geometry (mesh). A digital solution is acceptable when the convergence criterion is established, that is to say when the numerical result is very close to the exact solution. This is precisely known, it is only available in some cases for simple geometries. Therefore an error estimate is constructed from a Taylor expansion of the second order differential of the equation system operator. Convergence is reached when the estimated error value is less than a threshold value.

The error for the non-linear case is where N is the number of nodes, E is the estimation error, and W is the weight of each node (equal to 1 by default). This value must be less than a factor defined by the user (the default value is 10-6).

err= ( 1 N i ( | E i | W i ) 2 ) MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaamyzaiaadkhacaWGYbGaeyypa0ZaaOaaa8aabaWdbmaabmaapaqa a8qadaWcaaWdaeaapeGaaGymaaWdaeaapeGaamOtaaaadaGfqbqabS WdaeaapeGaamyAaaqab0WdaeaapeGaeyyeIuoaaOWaaeWaa8aabaWd bmaalaaapaqaa8qadaabdaWdaeaapeGaamyra8aadaWgaaWcbaWdbi aadMgaa8aabeaaaOWdbiaawEa7caGLiWoaa8aabaWdbiaadEfapaWa aSbaaSqaa8qacaWGPbaapaqabaaaaaGcpeGaayjkaiaawMcaa8aada ahaaWcbeqaa8qacaaIYaaaaaGccaGLOaGaayzkaaaaleqaaaaa@4DFB@

For a time calculation, it is necessary that the solution satisfies the criterion E at each time step, where Ai is the absolute tolerance (default 10-3), and Ui is the relative tolerance (default 10-2).

rr= ( 1 N i ( | E i | A i +R| U i | ) 2 ) <1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaamOCaiaadkhacqGH9aqpdaGcaaWdaeaapeWaaeWaa8aabaWdbmaa laaapaqaa8qacaaIXaaapaqaa8qacaWGobaaamaawafabeWcpaqaa8 qacaWGPbaabeqdpaqaa8qacqGHris5aaGcdaqadaWdaeaapeWaaSaa a8aabaWdbmaaemaapaqaa8qacaWGfbWdamaaBaaaleaapeGaamyAaa WdaeqaaaGcpeGaay5bSlaawIa7aaWdaeaapeGaamyqa8aadaWgaaWc baWdbiaadMgaa8aabeaak8qacqGHRaWkcaWGsbWaaqWaa8aabaWdbi aadwfapaWaaSbaaSqaa8qacaWGPbaapaqabaaak8qacaGLhWUaayjc SdaaaaGaayjkaiaawMcaa8aadaahaaWcbeqaa8qacaaIYaaaaaGcca GLOaGaayzkaaaaleqaaOGaeyipaWJaaGymaaaa@55FA@

Convergence reflects the adequacy of the numerical solution with the approximate solution of the model. The quality (digital) of the obtained solution and its stability depend on the mesh refinement and time step to satisfy the CFL condition (Courant - Friedrichs -Lewy) presented in the following equation:

C= Δt ΔL max( | v m | )1 MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape Gaam4qaiabg2da9maalaaapaqaa8qacqGHuoarcaWG0baapaqaa8qa cqGHuoarcaWGmbaaaiaad2gacaWGHbGaamiEamaabmaapaqaa8qada abdaWdaeaapeGaaeODa8aadaWgaaWcbaWdbiaad2gaa8aabeaaaOWd biaawEa7caGLiWoaaiaawIcacaGLPaaacqGHKjYOcaaIXaaaaa@4B8C@

Where, Δt MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaeyiLdqKaamiDaaaa@3AAB@ is the time interval ΔL MathType@MTEF@5@5@+= feaagKart1ev2aqatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn hiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr 4rNCHbGeaGqkY=wjYJH8sqFD0xXdHaVhbbf9v8qqaqFr0xc9pk0xbb a9q8WqFfeaY=biLkVcLq=JHqpepeea0=as0Fb9pgeaYRXxe9vr0=vr 0=vqpWqaaeaabiGaciaacaqabeaadaqaaqaaaOqaaabaaaaaaaaape GaeyiLdqKaamitaaaa@3A83@ is the size range, and vm is the slow component of the velocity.

Conditions to the limits

Once the geometry of our clarifier system is defined and meshed, it fixed the boundary conditions of the system. After setting the initial mass concentration of the sewage sludge to 0003Kg/m3, is fed the input speed of the waste water to 1.25m/s and pressure is fixed at the outlet of treated water at atmospheric pressure. Then, we choose the most representative type of flow of what is happening inside the circular clarifier and was chosen in this case to turbulent flow k type ε multiphasic. Comsol Multiphysics after execution of the computer code on our system that consists of solving the set of equations governing the digital decanter and after verification of convergence curves that provide information on the reliability of the results obtained solutions were recovered in different forms. To demonstrate the appearance phase separation in the decanter under the effect of gravity, the mass concentrations showing the waste water in Kg/m3 dans the assumed two-phase basin as shown in Figure 3 during the settling period, which lasts for about 24 hours. Can be obtained from solutions of other information on other physical parameters such as speed lines, the current flow lines.

Figure 3 Streamlines and concentration of the solid phase during the settling.

Results and discussion

Figure 3 shows stream lines of the mixing rate and the solid phase concentration, measured in kg/m3, after 24 hours. At this point, the flow has reached a stable state. We can interpret the speed of the average mixing time. As expected, the solid particles fall down. However, the turbulent fluctuations results in the diffusion of the solid phase have a negative effect on the separation (Figure 4).

Figure 4 Mass flow of the dispersed phase at the inlet (blue), output device (green) and central outlet (red).

The clarifier and removes solid particles Min rates-Mut. Calculation of the rate of withdrawal from the results shows that the clarifier removes 0518-.151= .367kg of solid particles per second.

The results of numerical simulations clearly show the effectiveness of the clarifier, which forces the driven solid particles of gravity to get to the pond bottom, leaving the surface relatively light-weight solid effluent water. Also, it confirms what was reported previously on the importance of the wall introduced into the clarifier, where the impact is observed on device performance by preventing the clearest water from coming into contact with the sludge piled deep. Depending on requirements, we can get more details on a particular region or regions of the decanter to improve and optimize there after the digital system.

Conclusion

Marine pollution by wastewater is a serious problem in the Mediterranean Sea. The country's defining reject most of their liquid effluents that this is from the sanitation water towns or discharges from industry without any treatment and Algeria is one. Failure prediction and lack of stations and sewage systems waste water just before the discharge of liquid waste into the sea, causing the deterioration of the overall health and air quality of coastal waters. This causes significant disruption of ecosystems.

Clearances of Fertial-Arzew society as industrial waste water discharged directly into the sea. Although it has the ISO 14001 certificate, the release of this company pollution results of measurements of parameters indicative reveal that study sites are more or less polluted by wastewater discharges. However, prolonged contamination of sea water with these high values ​​of MES up to 260mg/L degrade the marine ecosystem and harm wildlife and existing marine flora.

Measurements of MY show the need for pretreatment of raw wastewater before sending them into the sea and improve their quality according to the required standards as well as to protect the environment and human health. To remedy this, we proposed a wastewater treatment system called "clarifier" or "settler", which allows a considerable reduction in the values ​​of MES. The values ​​found by numerical simulation under the calculation code Comsol Multiphysics, clearly showed the improvement of the quality of treated water from those who were not treated at all. The clarifier separates the suspended impurities of clean water. The latter is sent to the sea via the output device and recovered at the bottom outlet of the purification system sediments as sludge. Thus, the sludges are prevented from going contaminate the sea water and can be upgraded later. Our well designed system reduced the amount of suspended solids up to 60%. The fact that contaminated water were retained in the clarifier tank for 24 hours, before sending them into the sea for settling.

Acknowledgments

None.

Conflicts of interest

Authors declare that there is no conflict of interest.

References

  1. Teh CY, Wu TY, Juan JC. Potential use of rice starch in coagulation-flocculation process of agro-industrial wastewater: Treatment performance and floc characterization. Ecological Engineering. 2014 ;71:509–519.
  2. Mompoint M. Evaluation of environmental hazards generated by urban liquid effluents on the ecosystem of the bay of Port-au-Prince: First methodological approach. Quisqueya University-Civil Engineer ; 2004.
  3. Renuka N, Sood A, Prasanna R, et al. Phycoremediation of wastewaters: a synergistic approach using microalgae for bioremediation and biomass generation. International Journal of Environmental Science and Technology. 2015;12(4):1443–1460.
  4. Espinosa F, Garcia-Garcia JM, García-Gómez JC. Sewage pollution and exinction risk: an endangered limpet as bioindicator? Biodiversity and Conservation. 2007;16(2) :377–397.
  5. Fares R, Lounis M. Determination of the quality of sea waters Arzew, Algeria Gulf. Phytochem & BioSub Journal. 2017;11(2):110–117.
  6. Remili Kerfouf SA. Environmental monitoring and coastal sustainable development (case of Oran coast). Proceedings of the 5th International Symposium on the theme "Energy, Climate Change and Sustainable Development", Tunisia : Hammamet ; 2009.
  7. Taghezout Fatima. Environmental impact on the water releases the long littoral Western Algerian, Magister in: environmental sciences, faculty of natural sciences and life biology department laboratory networking environmental monitoring (ERA will) USTO-MB. 2015.
  8. Ridha Hachicha. Determination of priority measures for pollution of the coastal area south of Grand Sfax except SIAPE, Scenarios & remediation Environmental Regulations. Tunisia SMAP III project (2006-2008), Strategy Management Integrated Coastal Zone South Grand Sfax; 2008.
  9. Fares R, Louns M. Pollution characterization of liquid wastes of the factory complex Fertial (Arzew, Algeria). J Air Wast Manage Assoc. 2016;66(3):260–266.
  10. Comsol Multiphysics User's Guide © Copyright 1994-2008. Comsol AB.
Creative Commons Attribution License

©2019 Redouane, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.