AAA Abstract and Applied Analysis 1687-0409 1085-3375 Hindawi Publishing Corporation 384640 10.1155/2013/384640 384640 Research Article Modeling and Simulation of Flow and Formation Damage of Asphalt-Paved Roads Alawi M. H. 1 El-Qadi M. M. 2 El-Ameen M. A. 3 El-Amin Mohamed Fathy 1 Civil Engineering Department Collage of Engineering Umm Al-Qura University Saudi Arabia uqu.edu.sa 2 Mathematics Department Faculty of Science Helwan University Egypt helwan.edu.eg 3 Preparing Year Department Collage of Engineering Umm Al-Qura University Saudi Arabia uqu.edu.sa 2013 20 11 2013 2013 29 09 2013 05 10 2013 2013 Copyright © 2013 M. H. Alawi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Porous asphalt is a standard asphalt built on aggregate storage bed which allows water to drain through it and reduces stormwater runoff. However, porosity of the porous asphalt and the storage bed may be effectively reduced due to trapping suspended solids from the water or from the asphalt damage. In this paper, we present mathematical modeling and numerical simulation of flow and damage of porous asphalt-paved roads. A mathematical model to describe the fine-particles transport carried by a two-phase flow in a porous medium is presented. The buoyancy, capillarity, and mixed relative permeabilities correlations to fit with the mixed-wet system are considered. Throughout this investigation, we monitor the changing of the fluids properties such as water saturation and solid properties such as porosity and permeability due to trapping the fine-particles.

1. Introduction

The impervious asphalt-paved roads may store significant amounts of thermal energy during summer. The stored thermal energy may be transferred to stream waters during runoff events. Moreover, as a result of fluctuations in ambient air temperatures—diurnal and seasonal, intensity of solar radiation of the asphalt-paved, a significant deformation of asphalt-paved may occur due to the heat and mass transfer between asphalt and water. On the other hand porous (pervious, permeable, or open-graded) asphalt is standard asphalt with reduced sand or fines that allows water to drain through it. Pervious asphalt built on aggregate storage bed reduces stormwater runoff. In addition to reducing runoff, this effectively traps suspended solids and filters pollutants from the water. The stormwater flows through the asphalt to the layer of crushed stone aggregate bedding and base that supports the asphalt while providing storage and runoff treatment. The use of porous asphalt can potentially reduce additional expenditures and land consumption for conventional collection, conveyance, and detention stormwater infrastructure. Compared to the lifetime of dense graded asphalt concrete roads, the lifetime of porous asphalt concrete is less . The loss of stones from the road surface, called raveling, is mostly reported as the dominant defect in porous asphalt wearing course [2, 3]. Raveling is the start of major defects like potholes, because once a stone is gone, the surrounding stones will follow for lacking support in at least one direction . Raveling also has negative influence on the noise reduction function and skid resistance of porous asphalt pavement. During service life, the pores tend to be clogged by dirt, dust, or other clogging agents. The formulation of fine-particles transport in two-phase flow in porous media has been studied experimentally and numerically in . Formation damage is a common problem in reservoir development. The particles migrate through the porous media, deposit on the pore surfaces, and become trapped at pore constrictions to reduce the rock porosity and permeability.

Researchers have done many attempts to predict these phenomena experimentally and numerically; however, there is no study that considers the possible porosity and permeability reduction. So, in order to keep porous asphalt and its storage bed efficient, considering these kinds of effects may lead to interesting results. In the current work, we introduce modeling and numerical simulation of fine-particle transport in two-phase flow in porous asphalt-paved roads.

2. Mathematical Modeling

The basic equations that govern the flow of the two-phase flow in porous media are mass conservation equation and constitutive equation (Darcy’s law). The two-dimensional governing equations may be written as(1a)(φρwSw)t=-·(ρwvw),(1b)(φρaSa)t=-·(ρava),(2a)vw=-Kkrwμw(pw-ρwg),(2b)va=-Kkraμa(pa-ρag),where S is the saturation and v [m/s] is the velocity. w stands for the wetting phase (water), and a stands for the nonwetting phase (air). φ is the porosity of the medium, and =(/x,/z) is the divergence operator. K [m2] is the absolute permeability, krα [dimensionless] is the relative permeability, ρ [kg m3] is the density, p [Pa] is the pressure, and g=(0,-g)T is the gravitational acceleration. μ is the viscosity. The fluid saturations for the two-phase flow are interrelated by (3)Sw+Sa=1. The total velocity is (4)vt=vw+va. Summing the saturation equation for water phase and the oil phase, one obtains (5)·vt=0. Moreover, adding the constitutive equations for each phase, (1a)–(2b), and substituting into (5), we end up with (6)vt=-K[λt(S)pw+λa(S)pc]+Kχ(S)g, where λα(S)=krα(S)/μα is the mobility, λt(S)=λa(S)+λw(S) is the total mobility, and χ=λa(S)ρa+λw(S)ρw. In order to derive the pressure equation, substitute (6) into (5). We obtain (7)·[λt(S)pw+λa(S)pc(S)-χ(S)g]=0. Substituting the constitutive equation of the water phase, (2a) and (2b), into (1a) and (1b) gives (8)(φSw)t+·vw=0. Therefore, water velocity may be written as (9)vw=favt+λwK[fapc]+λwfaKΔρg, where Δρ=ρw-ρa.

On the other hand, a mathematical model is developed to describe the fine-particles transport carried by two-phase flow in porous media. Assuming that we have a number m of size intervals of particles in the water phase, the transport equation for each size interval i of the particles in the water phase can be written as (10)(ϕSwCi)t+vw·Ci=·(φSwDiCi)+Ri, where i=1,2,,m. Di is the dispersion coefficients of particles in size interval i in the water phase. Ri is the net rate of loss of particles in size interval i in the water phase. The net rate of loss of particles may be written as  follows: (11)Ri=(δϕ)it, where (δϕ)i=vi+vi* is the porosity change due to release or retention of particles of interval i in the water phase. vi is the volume of the particles of interval size i in contact with the water phase available on the pore surfaces per unit bulk volume of sandstone. vi* is the volume of the particles of interval size i entrapped in pore throats from the water phase per unit bulk volume of sandstone due to plugging and bridging. At the critical velocity of the surface deposition only particle retention occurs while above it retention and entrainment of the particles take place simultaneously . A modified Gruesbeck and Collins’s model for the surface deposition is expressed by  the following: (12)vit={αd,i|vw|Ciwhen  vwvc,αd,i|vw|Ci-αe,ivi|vw-vc|when  vw>vc,  where αd,i is the rate coefficients for surface retention of the particles in interval i. αe,i is the rate coefficients for entrainment of the particles in interval i. vc is the critical velocity. Similarly, the rate of entrapment of the particles in interval i is (13)vi*t=αpt,i|vw|Ci, where αpt,i is the pore throat blocking constants. Porosity may be changed because of particles deposition on the pore surfaces or blocking of pore throats. The porosity variation may be expressed by  the following: (14)ϕ=ϕ0-(δϕ)i, where ϕ0 is the initial porosity. Also, the permeability variation due to particles deposition on the pore surfaces or blocking of pore throats may be expressed as  follows: (15)K=K0((1-f)kf+fϕϕ0)l, where K0 is the initial permeability, kf is a constant for fluid seepage allowed by the plugged pores, f=1-iαf,ivi* is the fraction of the original cross-sectional area open to flow, and the value of the exponent l has range from 2.5 to 3.5. For the particles transport carried by fluid stream in the porous media, deposition on pore surfaces and blockage in pore throats may occur. The retained particles on pore surfaces may desorb for hydrodynamic forces and then possibly adsorb on other sites of the pore bodies or get entrapped at other pore throats.

3. Results and Discussion

In order to get physical insights for the problem under consideration, we consider that typical cross-section for pervious pavement system consists of six porous layers (Figure 1) that can be simplified to one-dimensional system for calculations. The porous asphalt layer taken is 15 cm, the choker course layer is 20 cm, the filter course is 30 cm, the filter blanket is 8 cm, the reservoir course is 10 cm, and finally the native material is optional that we consider it 20 cm. We consider the top boundary as inlet of water-particles suspension. So water imbibes into the porous asphalt and replaces the air that leaves the medium from the same boundary. This type of flow in porous media is called countercurrent imbibition . In countercurrent imbibition both wetting and nonwetting phases flow through one inflow-outflow boundary. Therefore, the total velocity becomes zero, vt=0. Therefore, the flow equations may be written as (16)Pwz=-fw(S)Pc(S)z-χ(S)g,vw=Kλwfa(Pcz+Δρg). Therefore, (17)(φSw)t+Kz[λwfa(Pcz-Δρg)]=0. The relative permeabilities are given in terms of the normalized wetting phase saturation S, given as (18)krw=krw0Sa,kra=kra0(1-S)b, where a and b are positive numbers, krw0=krw(S=1) is the endpoint relative permeability to wetting phase, and kra0=kra(S=0) is the endpoint relative permeability to nonwetting phase. The normalized wetting phase saturation is given by (19)S=Sw-Siw1-Sra-Siw,0S1, where Siw is the irreducible water saturation and Sra is the residual air saturation.

Typical cross-section for pervious pavement system.

For the capillary pressure we used the general correlation (20)pc=cw(Sw-Swr1-Swr)-aw+ca(Sa-Sar1-Sar)-aa, where cw and ca are constants that represent the entry pressure for imbibition and drainage, respectively. The constants 1/aw and 1/aa are the pore size distribution indexex for imbibition and drainage, respectively.

The transport equation for the interval i of the particles can be written as (21)(φSwCi)t+vwCiz=z(φSwDiCiz)+Ri,Ri={(αd,i+αpt,i)vwCi,uwuc(αd,i+αpt,i)vwCi-αe,ivivw-vcCi,uw>uc. The surface deposition of the particles in the interval i is (22)vit={αd,ivwCi,vwucαd,ivwCi-αe,ivivw-vcCi,vw>uc. The rate of entrapment of the particles in interval i is (23)vi*t=αpt,ivwCi. The initial conditions are (24)Sw=S0w,Ci=vi=vi*=0,t=0,0zH, where H is the depth and S0w is the initial water saturation. The boundary conditions are (25)Sw=1-Sor,Ci=Ci,0,vi=vi*=0,t0,z=0,dSwdz=dCidz=dvidz=dvi*dz=0,t0,z=H, where Ci,0 is the concentration of particles in the particles suspension at the inlet boundary. The governing equations (17)–(23) are solved numerically along with the initial and boundary conditions, (24)–(25). An efficient algorithm is used to solve the above high-nonlinear parabolic partial differential equation in one space variable z and time t. The Galerkin method is used for spatial discretization , while the time integration for the resulting ordinary differential equation is done with an adaptive time step. 100 points of the spatial grid were used during calculations and were enough to provide an acceptable accuracy. Now, we consider one-size particles suspension in the water phase at the inlet, with the following parameter [6, 7], αd,i=8 m−1, αpt,i=17 m−1, αe,i=400 m−1, vc=4.6×10-3 m/s, and Di=5.6×10-3 m2/s. The inlet particles concentration values are Ci,0=0.0 (without particles), 0.0005, and 0.01. The remaining model parameters are Swr=Sor=0.001, aw=aa=0.5, cw=500, ca=-500, φ0=0.2, l=3, kf=0.015, γf=0.01, krw0=kra0=1, and a=b=4.

Figure 2 shows the water saturation against the dimensionless distance with various imbibition times and inlet particles concentrations. It is notable that the concentration of particles increases the water saturation in particular after significant time of imbibition. Moreover, the saturation of water increases with the time imbibition.

Water saturation against the height with various imbibition times and inlet concentrations.

The normalized particles concentration is plotted in Figure 3 against the dimensionless distance with various inlet concentrations and imbibition times. From this figure it can be seen that the concentration of particles increases as the concentration at the inlet increases for positive Ri. An interesting phenomenon can be observed from Figure 3 that the concentration increases with time but after a certain time of imbibition concentration of particles start to decrease with time. This may be interpreted by increasing the rate of particles precipitation on the pore surface of the porous medium. The opposite is true for minus Ri. It is noteworthy that, in Figure 3, the dimensionless concentration goes larger than 1 because the concentration of particles increases by erosion of more particles from the medium (-Ri). On the other hand, the reference concentration was the inlet concentration of the particles coming with water from outside.

Normalized particles concentration against height with positive/minus Ri and various imbibition times.

Figures 4 and 5 illustrate the ratios of the permeability and porosity against the dimensionless distance with various imbibition times, respectively. Both the permeability and the porosity are reduced due to the precipitation of the particles on the pores walls. It is interesting to note that the reduction rates of the permeability and porosity depend on the layer permeability of the medium.

Permeability ratio against the height with various imbibition times.

Porosity ratio against the height with various imbibition times.

4. Conclusions

In this pape, we presented numerical modeling and simulation of particle transport in two-phase flow in porous asphalt-paved roads. Numerical experiments have been performed to explore these phenomena and to study the possible porosity and permeability variations. We found that both the permeability and the porosity are reduced due to the precipitation of the particles on the pores walls. These results may help engineers to keep porous asphalt and its storage bed efficient by looking for cleaning methodologies to avoid media blocking.

Acknowledgment

The authors would like to thank the Institute of Scientific Research and Revival of Islamic Heritage, Umm Al-Qura University, for supporting the Project no. 43208015.

Voskuilen J. L. M. Verhoef P. N. W. Cause of premature ravelling failure of porous asphalt Proceedings of the 6th International RILEM Symposium on Performance Testing and Evaluation of Bituminous Materials 2003 191 197 Padmos C. Over ten years experience with porous road surfaces Proceedings of the ISAP 9th International Conference on Asphalt Pavements 2002 Copenhagen, Denmark Voskuilen J. L. M. Huurman M. Conversations, Centre For Transport and Navigation of the Dutch Ministry of Transport 2009 Delft, The Netherlands Public Wprks and Water Management Kneepkens A. van Hoof T. Schaefer H. van Keulen W. VIA-RAL for Porous Asphalt: A Result of research and Development, But Most of all of implementation 2004 Delft, The Netherlands Wegbouwkundige Werkdagen Liu X. H. Civian F. Characterization and prediction of formation damage in two-phase flow systems Proceedings of the SPE Production Operations Symposium 1993 Oklahoma City, Okla, USA 10.2118/25429-MS Liu X. Civan F. Multiphase mud fluid infiltration and filter cake formation model Proceedings of the SPE International Symposium on Oilfield Chemistry March 1989 607 621 2-s2.0-0027206302 Liu X. Farouk C. Formation damage and skin factor due to filter cake formation and fines migration in the near-wellbore region Proceedings of the International Symposium on Formation Damage Control February 1994 Lafayette, Ind, USA 259 273 2-s2.0-0028255501 Gruesbeck C. Collins R. E. Entrainment and deposition of fines particles in porous media Society of Petroleum Engineers Journal 1982 22 6 847 856 2-s2.0-0020293353 El-Amin M. F. Salama A. Sun S. Numerical and dimensional investigation of two-phase countercurrent imbibition in porous media Journal of Computational and Applied Mathematics 2013 242 285 296 10.1016/j.cam.2012.09.035 MR2997443 ZBL06127790 Sun S. Salama A. El Amin M. F. Matrix-oriented implementation for the numerical solution of the partial differential equations governing flows and transport in porous media Computers & Fluids 2012 68 38 46 10.1016/j.compfluid.2012.07.027 MR2979705 El-Amin M. Salama A. Sun S. Effects of gravity and inlet location on a two-phase countercurrent imbibition in porous media International Journal of Chemical Engineering 2012 2012 7 210128 10.1155/2012/210128 El-Amin M. F. Sun S. Effects of gravity and inlet/outlet location on a two-phase cocurrent imbibition in porous media Journal of Applied Mathematics 2011 2011 18 673523 10.1155/2011/673523 MR2771482 ZBL1213.76191 Skeel R. D. Berzins M. A method for the spatial discretization of parabolic equations in one space variable SIAM Journal on Scientific and Statistical Computing 1990 11 1 1 32 10.1137/0911001 MR1032224 ZBL0701.65065