Research Article Volume 2 Issue 1
Department of Civil Engineering, Gregory University Uturu, Nigeria
Correspondence: Eluozo SN, Department of Civil Engineering, Gregory University Uturu (GUU)Abia State of Nigeria, Nigeria
Received: July 28, 2017 | Published: January 24, 2018
Citation: Eluozo SN, Afiibor BB, Amagbo L. Numerical modeling to monitor dispersion pressure on the deposition of arthrobacter in predominant phreatic bed. MOJ App Bio Biomech. 2018;2(1):23–27. DOI: 10.15406/mojabb.2018.02.00040
The deposition of Arthrobacter deposition were monitor base on the rate of dispersion in the study environment, the study express the behaviour of Arthrobacter deposition in phreatic bed through the influences of lithology observed in heterogeneous formation, such condition were observed to monitor the rate of pressure on the transport system, application of numerical simulation were applied, the results express fluctuation and exponential phase on the transport process to phreatic beds, the study developed series of concentration variation base of the refection from geochemistry depositions, thus reaction with the contaminant in terms of inhibition or supply of nutrient for its population. The study is imperative because it has express the behaviour of Arthrobacter in its depositions and transport process in phreatic beds.
Keywords: numerical modeling, dispersion, arthrobacter, phreatic bed
To determine if a given water supply is safe, the source needs to be protected and monitored regularly. There are two broad approaches to water quality monitoring for pathogen detection. The first approach is direct detection of the pathogen itself, for example, the protozoan Cryptosporidium parvum. While it will be more accurate and precise if specific disease-causing pathogens are detected directly for the determination of water quality, there are several problems with this approach. First, it would be practically impossible to test for each of the wide variety of pathogens that may be present in polluted water. Second, even though most of these pathogens can now be directly detected, the methods are often difficult, relatively expensive, and time-consuming (WHO, 1996). Instead, water monitoring for microbiological quality is primarily based on a second approach, which is to test for .indicator organisms. For a classification table created by the author of typical indicator organisms) The indicator microorganisms should fulfill the following criteria Stetler.1 The concept which explains the ultimate destination of rainwater is the sea either directly through run off or indirectly be infiltration and subsurface flow. A system of water movement in the atmosphere or rainfall, dews, hailstones or snowfalls over land as run off. Vertical and horizontal movement underground as infiltration or subsurface and continuous movement of all forms of water is the hydrogeology cycle. In the atmosphere, water vapours condense and may give rise to precipitation. However, not all this precipitation will reach the ground surface; some are intercepted by vegetation cover or surface of building and other structures and then evaporate back into the atmosphere. The precipitation that reaches the ground surface may flow in to stream, lake and ocean, where it will either be evaporated or form seepages intruding in to the ground likewise soil moisture and further percolate downward to underline aquifer where it may be held for several years longer. Groundwater in Nigeria is restricted by the fact that more than half of the country is underlain by crystalline basement rock of pre-cambian era. The main rock types in this geological terrain include igneous and metamorphic rock such as migmatites and granite gneisses. Generally in their unaltered form, they are characterized by low porosity and permeability. Porosity in basement rocks is by induction through weathering while secondary permeability induces by tectonic activities which manifest in form of that often act as conduct path facilitating water movement. In other words, aquiferous zones in the basement terrain include fractured/weathered rocks. The yielding capacity of well, drilled within such rock are always very enormous. (Shitta 2007) Groundwater is the main resource of drinking water in many parts of the world. Contamination resulting from industry, urbanization and agriculture poses a threat to the groundwater quality Amadi.2,3 The task of balancing groundwater protection and economic activities is challenging. Therefore, understanding the effects of different water management strategies and the role of climate change is essential for the sustainable use of coastal groundwater resources (Prasad and Narayana, 2004). According to Olobaniyi and Owoyemi (2006), the coastal regions of the world are the most densely populated areas in the world. More than one third of the world’s populations are living within 100 km of the coastline (Hughes, et al., 1998). At the same time, the coastal regions provide about one third of the world’s ecosystem services and natural capital Aris et al.4 Such growth is accompanied by increasing demand for water supply leading to the over-exploitation of the aquifer system and excessive drainage for land reclamation purposes. Contamination of the groundwater by natural means (seawater intrusion) and through anthropogenic means (human activities) cannot be ruled out in the area. The study is aimed at evaluating the quality of groundwater from the coastal plain-sand aquifer in Port-Harcourt area with the view of determining its suitability for domestic, irrigational and industrial purposes. The heavy industrial and human activities in the area lead to the present study. The aquifer system in the area is largely unconfined, highly porous and permeable and the possibility of anthropogenic interference cannot be completely ignored, hence the need for this study. Port-Harcourt, the ‘garden-city and treasure base of the nation’ is situated about 60km from the open sea lies between longitude 6o55’E to 7o10’E of the Greenwich meridian and latitude 4o38’N to 4o54’N of the Equator, covering a total distance of about 804km2 Akpokodje.5 In terms of drainage, the area is situated on the top of Bonny River and is entirely lowland with an average elevation of about 15m above sea level.6 The topography is under the influence of tides which results in flooding especially during rainy season.7
Climatically, the city is situated within the sub-equatorial region with the tropical monsoon climate characterized by high temperatures, low pressure and high relative humidity all the year round. The mean annual temperature, rainfall and relative humidity are 30oC, 2,300mm and 90% respectively (Ashton-Jones, 1998). The soil in the area is mainly silty-clay with interaction of sand and gravel while the vegetation is a combination of mangrove swamp forest and rainforest.8 Port-Harcourt falls within the Niger Delta Basin of Southern Nigeria which is defined geologically by three sub-surface sedimentary facies: Akata, Agbada and Benin formations.9,10 The Benin Formation (Oligocene to Recent) is the aquiferous formation in the study area with an average thickness of about 2100m at the centre of the basin and consists of coarse to medium grained sandstone, gravels and clay with an average thickness of about 2100m at the centre of the basin and consists of coarse to medium grained sandstone, gravels and clay Etu-Efeotor & Akpokodje.11 The Agbada Formation consists of alternating deltaic (fluvial, coastal, fluviomarine) and shale, while Akata Formation is the basal sedimentary unit of the entire Niger Delta, consisting of low density, high pressure shallow marine to deep water shale.12
∂C∂t=φV∂C∂x+D∂2C∂x2+qLINAC∂C∂t=φV∂C∂x+D∂2C∂x2+qLINAC(1) But
φV=φV=Velocity, v in meter per second (m/s), and porosity [-]. Thus equation (1) becomes:∂C∂t=v∂C∂x+D∂2C∂x2+qLINAC∂C∂t=v∂C∂x+D∂2C∂x2+qLINAC (2)
Converting the PDE to its algebraic equivalent equation by applying the finite different approximation technique for the implicit scheme, we obtain as follows.∂C∂t=Cj+1i−CjiΔt∂C∂t=Cj+1i−CjiΔt (3)
∂C∂x=Cj+1i+1−Cj+1i−12Δx∂C∂x=Cj+1i+1−Cj+1i−12Δx
(4)∂2C∂x2=Cj+1i+1−2Cj+1i+Cj+1i−1Δx2∂2C∂x2=Cj+1i+1−2Cj+1i+Cj+1i−1Δx2 (5)
Substituting equation (3) through (5) into (2) gives:
Cj=1i−CjiΔt=v[Cj+1i+1−Cj+1i−12Δx]+D[Cj+1i+1−2Cj+1i+Cj+1i−1Δx2]+qLINACj+1iCj=1i−CjiΔt=v[Cj+1i+1−Cj+1i−12Δx]+D[Cj+1i+1−2Cj+1i+Cj+1i−1Δx2]+qLINACj+1i
Cj+1i−Cji=Δt2Δxv[Cj+1i+1−Cj+1i−1]+ΔtDΔx2[Cj+1i+1−2Cj+1i+Cj+1i−1]+ΔtqLINACj+1lCj+1i−Cji=Δt2Δxv[Cj+1i+1−Cj+1i−1]+ΔtDΔx2[Cj+1i+1−2Cj+1i+Cj+1i−1]+ΔtqLINACj+1l
Cj+1i−Cji=λ(Cj+1i+1−Cj+1i−1)+K(Cj+1i+1−2Cj+1i+Cj+1i−1)+αCj+1iCj+1i−Cji=λ(Cj+1i+1−Cj+1i−1)+K(Cj+1i+1−2Cj+1i+Cj+1i−1)+αCj+1i
OrCji+(α−λ−2K−1)Cj+1i+(λ+K)Cj+1i+1+KCj+1i−1=0Cji+(α−λ−2K−1)Cj+1i+(λ+K)Cj+1i+1+KCj+1i−1=0 (6)
For cases where the initial and final conditions are given, boundary condition at the first node can be expressed as:
Cj+10=f0(tj+1)Cj+10=f0(tj+1)(7a)
Hence, first node equation is expressed as:Cji+(α−λ−2K−1)Cj+1i+(λ+K)Cj+1i+1=−Kf0(tj+1)Cji+(α−λ−2K−1)Cj+1i+(λ+K)Cj+1i+1=−Kf0(tj+1) (7b)
Cj+1l=fl+1(tj+1)Cj+1l=fl+1(tj+1)(8a)
Similarly, the last node boundary condition is:Cjl+(α−λ−2K−1)Cj+1l+KCj+1l−1=−(λ+K)fl+1(tj+1)Cjl+(α−λ−2K−1)Cj+1l+KCj+1l−1=−(λ+K)fl+1(tj+1) (8b)
For 1≤x≤9and0≤t≤4; and for the first instance, we obtain as follows:
Attime=0(i.ej=0):
i=1,
C01+KC10+(α−λ−2K−1)C11+(λ+K)C12=0 (9a)
i=2,
C02+KC11+(α−λ−2K−1)C12+(λ+K)C13=0(9b)
i=3 ,C03+KC12+(α−λ−2K−1)C13+(λ+K)C14=0(9c)
i=4 ,
C04+KC13+(α−λ−2K−1)C14+(λ+K)C15=0(9d)
i=5,
C05+KC14+(α−λ−2K−1)C15+(λ+K)C16=0 (9e)
i=6 ,
C06+KC15+(α−λ−2K−1)C16+(λ+K)C17=0(9f)
i=7,
C07+KC16+(α−λ−2K−1)C17+(λ+K)C18=0 (9g)
i=8 ,
C08+KC17+(α−λ−2K−1)C18+(λ+K)C19=0 (9h)
i=9
C09+KC18+(α−λ−2K−1)C19=−(λ+K)f10(t1) (9i)
Atime,t=0,C01=C02=C03=C04=C05=C06=C07=C08=C09=0
Arranging equations (6a) through (6i) in vector matrix gives:
[ωλ+K0000000Kωλ+K0000000Kωλ+K0000000Kωλ+K0000000Kωλ+K0000000Kωλ+K0000000Kωλ+K0000000Kωλ+K0000000Kω]{C11C12C13C14C15C16C17C18C19}={−Kf0(t1)0000000−(λ+K)f10(t1)}
Where:
ω=(α−λ−2K−1)
Hence, at any point with time, the general form of the above equation is presented as:
[ωλ+K0000000Kωλ+K0000000Kωλ+K0000000Kωλ+K0000000Kωλ+K0000000Kωλ+K0000000Kωλ+K0000000Kωλ+K0000000Kω]{Cj+11Cj+12Cj+13Cj+14Cj+15Cj+16Cj+17Cj+18Cj+19}={−Kf0(tj+1)0000000−(λ+K)fl+1(tj+1)}
Depth (m) |
Concentration (g/L) |
0 |
0.653 |
3 |
0.689 |
6 |
0.6526 |
9 |
0.512 |
12 |
0.316 |
15 |
0.125 |
18 |
0.064 |
21 |
0.027 |
24 |
0.008 |
27 |
0.001 |
30 |
0.00025 |
Table 1 Simulation Values from Arthrobacter Concentration at Different Depth
Depth (m) |
Concentration (g/L) |
0 |
0.5543 |
3 |
0.3675 |
6 |
0.2270 |
9 |
0.1263 |
12 |
0.0592 |
15 |
0.0193 |
18 |
0.0036 |
21 |
0.0024 |
24 |
0.0022 |
27 |
0.0015 |
30 |
0.0005 |
Table 2 Simulation Values from Arthrobacter Concentration at Different Depth
Depth (m) |
Concentration (g/L) |
0 |
0.653 |
3 |
0.659563 |
6 |
0.679649 |
9 |
0.714496 |
12 |
0.766303 |
15 |
0.838469 |
18 |
0.935964 |
21 |
1.065902 |
24 |
1.238402 |
27 |
1.467884 |
30 |
1.775038 |
Table 3 Simulation Values from Arthrobacter Concentration at Different Depth
Depth (m) |
Concentration (g/L) |
0 |
0.0956 |
3 |
0.0875 |
6 |
0.0821 |
9 |
0.0794 |
12 |
0.0794 |
15 |
0.0819 |
18 |
0.0869 |
21 |
0.0943 |
24 |
0.1041 |
27 |
0.1163 |
30 |
0.1308 |
Table 4 Simulation Values from Arthrobacter Concentration at Different Depth
Depth (m) |
Concentration (g/L) |
0 |
0.11882 |
3 |
0.1577 |
6 |
0.1691 |
9 |
0.1560 |
12 |
0.1371 |
15 |
0.1074 |
18 |
0.0777 |
21 |
0.0548 |
24 |
0.0457 |
27 |
0.0571 |
30 |
0.0960 |
Table 5Simulation Values from Arthrobacter Concentration at Different Depth
The rate of Arthrobacter transport has been thoroughly monitored applying numerical modeling and simulation. The study applied the concept to monitor the rate of Arthrobacter concentration in discretazed phase, these concept monitor the system thoroughly with change in concentration in discretazed condition. The study express various migration rate in the simulation, these was base on the lithology of the study environment, the system monitor the effect on heterogeneous lithology to phreatic bed, high to low concentration was experiences in few locations, while exponential phase of the contaminant migration were observed. Dispersion of the contaminant were observed to developed the rate of transport in different phase, the lithology were found to developed significant effect in the deposition, these are base on the degree of deposited formation characteristics.
Author declares that there is no conflict of interest.
©2018 Eluozo, 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.