• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Adsorptive potential of Acacia nilotica based adsorbent for chromium(VI)from an aqueous phase

    2018-05-26 07:29:36RazaKhalidZaheerAslamAamirAbbasWaqarAhmadNaveedRamzanReyadShawabkeh

    Raza Khalid ,Zaheer Aslam ,*,Aamir Abbas ,Waqar Ahmad ,Naveed Ramzan ,Reyad Shawabkeh

    1 Department of Chemical Engineering,University of Engineering and Technology Lahore,54890,Pakistan

    2 Department of Chemical Engineering,King Fahd University of Petroleum and Minerals,Dhahran 31261,Saudi Arabia

    1.Introduction

    Excessive amounts of wastewater containing heavy metals and other pollutants are discharged into an aquatic media resulting from rapid industrialization and urbanization.Heavy metals in the waste aqueous streams pose a serious threat to the aquatic environment due to their non-degradable nature and easy absorption into the living cells.By this way these can be bio-accumulated and biomagnified in higher trophic levels[1–3].The presence of chromium heavy metal in the industrial effluents is also of great worry.Chromium salts are used in various applications like the manufacturing of corrosive paints,as wood preservatives and in the treatment of hides.Other sources of chromium include metal plating,mining,preparation of photographic materials and steel production.Chromium exists in different oxidation states but Cr(III)(trivalent chromium)and Cr(VI)(hexavalent chromium)are more important in relevance to their environmental and toxicological properties.Trivalent chromium is 500 times less toxic than hexavalent chromium and the allowable limits in the industrial discharge stream are 5 mg·L?1for Cr(III)and 0.05 mg·L?1for Cr(VI)[4–7].Different techniques which are available and may be applied for the removal of chromium from industrial waste water are comprised of chemical precipitation[8–10],reverse osmosis[11],Nanofiltration,ultra filtration[12–14],ion exchange[15–17],electrochemical process[18],biological processes[19]and adsorption[20–22].

    Among all,adsorption is considered as a cost effective,sustainable and ecofriendly technique for heavy metal removal from aqueous streams.This method also givespromising results even atlowerconcentrations of heavy metals.Commercial activated carbon has been widely used for the abatement of heavy metals because of its high surface area.But the increased price of activated carbons drives the research to look for low cost readily available adsorbents[23,24].Large numbers of potentially low cost adsorbents including agro-based material,microorganisms,seaweeds,industrial wastes,and natural residues have been reported in the literature[25–38].

    Feedstock like peelings of citrus,sugarcane waste,coffee husk,lemon peelings,gram husk,and rice husk have been reported in the literature[39–44].Adsorbent based on waste tamarind hull was used to remove hexavalent chromium from aqueous solution,different adsorption parameters were studied and 99%removal of chromium was reported at optimum conditions[45].Azadirachta indica(neem)leaf powder effectively removes chromiumup to 96.4%using 160μm adsorbent size at an optimum dosage[46].Vikrant and Pant used Eucalyptus bark for chromium removal from aqueous solution and found optimum conditions as pH 2,concentration and time were 200 μg·g?1and 2 h respectively,and the efficiency of adsorbent was more than 99%[47].Acacia nilotica bark had been used as an effective adsorbent in batch process to remove chromium metal.Activated carbon derived from A.nilotica provided 93.1%of chromium metal removal,while in the case of untreated A.nilotica bark 68.1%of chromium was removed from synthetic solution[48].A.nilotica sawdust,which contains less tannin concentration,is an abundantly available waste by-product of the timber industry.Exploring its potentiality for pollutant abatement will notonly help in waste management but will also provide a cheap source of adsorbent.In the past,different types of saw dust had been used to remove hexavalent chromium Cr(VI)from waste water like Hevea brasiliensis sawdust,treated robust sawdust,Indian rose wood sawdust,and Beech sawdust.As referenced earlier,there is great potential for sawdust to remove heavy metals particularly Cr(VI)from industrial effluents.Sawdust is comprised of electron rich functional groups which provide active sites for heavy metals and one step acid or alkali activation further enhances the adsorption efficiency of raw sawdust.The focus of this research is to functionalize further the acid activated sawdust with aldehyde functional agent to improve the pollutant's removal efficiency.The research results presented herein are part of investigations conducted to turn the inefficient raw A.nilotica sawdust to an economical and effective adsorbent to remove hazardous hexavalent chromium.The effect of different process parameters like initial solution pH,initial chromium concentration,contact time,temperature,and adsorbent dose on adsorption capacity was evaluated.Moreover,sorption kinetics,isotherms and thermodynamics were studied to analyze the various aspects of chromium removal from aqueous solution.

    2.Materials and Methods

    2.1.Materials

    Analytical grade Potassium dichromate(K2Cr2O7,294.2 g·mol?1)was obtained from Sigma Aldrich and it was used for Cr(VI)ion source in distilled water for experimentation.Hydrochloric acid(37%ACS Reagent grade,ACROS Organics)and formaldehyde(37%by mass,Fisher Scientific)were used for the modification of raw material.0.1 mol·L?1Sodium hydroxide(Merck Millipore,Germany)and 0.1 mol·L?1Nitric acid(65%,Merck Millipore,Germany)were used to adjust the pH of the adsorbate samples.1000 μg·g?1Cr(VI)atomic absorption spectrometry standard solution(Fluka Chemicals)was diluted to required concentrations with double distilled water(DDW)to run the calibration experiments.

    2.2.Adsorbent preparation

    The A.nilotica(keekar)sawdust(named as R-SD)was collected from a local saw mill situated in the vicinity.Raw sawdust was sun dried for 2 days and then sieved.The same sized particles were obtained by passing material through 20 mesh screens(US Taylor standard Sieve).Subsequently the sieved material was washed with distilled water several times to remove the dust from sample.Then washed sawdust was oven dried at 80°C for 12 h.The dried material was saved for subsequent chemical treatment.

    500 g dried sawdust was treated with 1 mol·L?1HCl in ratio of 1:15(sawdust:HCl,W/V)in a beaker with continuous stirring on a hot plate magnetic stirrer for 3 h at 40°C.After stirring for a required time the sample was removed and kept soaked for 24 h at room temperature.The acid treated sawdust was filtered and washed with distilled water to remove the acid contents until the pH of filtrate was 6 to ensure that all the acid had been washed.Filtered material was dried in an oven at 60°C for 8 h,and then transferred to sealed ceramic bottles for further use.The HCl washed sawdust was divided into two portions;the first one was used in experimental work without further treatment(named as H-SD)while the second one was soaked and impregnated with 1%formaldehyde solution with ratio of 1:3(sawdust:formaldehyde,W/V)for 24 h.After one day of soaking,the material was washed repeatedly and filtered using What man no.42 filter paper to remove formaldehyde.The filtered cake was dried in an oven at 80 °C and sample was named as “HF-SD”.Both samples(H-SD and HF-SD)were tested for their potential efficiency of removing chromium pollutant from water by doing batch experimentation.

    2.3.Adsorbent characterization

    Various instruments were used to characterize the prepared adsorbents.Surface texture and morphology and spot analysis was done using Scanning Electron Microscopy(SEM)on an EDX analyzer equipped coupled with Energy Dispersive X-ray Spectroscopy(EDX).SEMmicrogram adsorbent specimen was obtained under15 kV accelerating voltage.Spot analyses were conducted in triplicate runs to determine the elemental composition.The FT-IR spectra of samples were recorded by using a Nicolet Avatar 360 FT-IR spectrometer,and sample pellet was made by mixing 1 mg of adsorbent to 100 mg of KBr powder and hydraulically pressed to 10 t·m?2.The thin pellet was oven dried at 102°C to prevent the interference of moisture.The FTIR machine was set to transmission mode.Thermal stability studies were carried out by using a Thermo-Gravimetric Analyzer(TGA)over the temperature range 20–900 °C at a heating rate of 10 °C·min?1under a nitrogen flow rate of 20 ml·min?1.The BET(Brunauer–Emmett–Teller)characterization of raw and modified sawdust was done using Micromeritics ASAP2020.Accurately weighted quantity of an adsorbent was degassed at 333 K for 1 h under 666.61 Pa(5 mm Hg)pressure to desorb the entrapped gas molecules.Subsequently N2adsorption–desorption measurements were recorded and specific surface area was calculated by the BET method and pore size distribution was calculated by the Barrett–Joyner–Halenda(BJH)method.The test to determine the point of zero charge(pHpzc)was determined by taking 100 ml of 0.01 mol·L?1KNO3solution in different Erlenmeyer flasks.In each flask,initial value of pH was maintained varying from 1 to 10 and 1.02 g of an adsorbent was added in each flask and placed on an orbital shaker at 120 rpm for 24 h and then put idly for 15 min to settle the adsorbent.Final pH of the solution(pHf)was recorded and difference between initialand final pH(pHi?pHf)was calculated and it was plotted against initial pH to get the pHpzc.

    2.4.Preparation of adsorbate solution

    The stock solution of hexavalent chromium of concentration 500 μg·g?1was prepared by dissolving 1.415 g of analytical grade Potassium dichromate(K2Cr2O7)in 1000 ml of deionized water.The stock solution was further diluted to obtain desired concentrations for succeeding batch experiments.

    2.5.Batch adsorption studies

    Batch adsorption studies were carried out to remove hexavalent chromium from synthetic solution by using raw and activated sawdust as an adsorbent.Erlenmeyer flasks of capacity 250 ml were used in experimentation,in each flask 100 ml of chromium solution was taken ata required dosage,various parameters were studied which had an effect on the adsorption process such as dosage in grams(0.25,0.7,1.5,2.5,3.5,4.5,6)with 100 mg·L?1of chromium concentration at 30 °C temperature.Kinetic studies were carried out at different time intervals ranging from 1 to 540 min,with optimum dose and 100 mg·L?1of chromium.Concentration of chromium in synthetic solution varied in mg·L?1(300,200,100,80,60,30,15,5)with optimum dose and time.Effect of pH of the solution was examined in the range of 1–10 with optimum dose,time,and concentration.The effect of temperature was observed by running the experiments at four different temperatures i.e.(20–50 °C with a gap of 10 °C).All the runs were carried out at fixed 130 r·min?1.Afterwards the solutions were filtered and filtrate was collected into sample bottles to find the remaining concentration of hexavalent chromium in synthetic solution.Chromium concentrations were analyzed by an Atomic Absorption Spectrophotometer(AA-6800 SHIMADZU Japan).Percentage removal of hexavalent chromium(R,%)from synthetic solution with raw and activated sawdust was calculated by using the following equation,

    where Ciis the concentration of chromium(mg·L?1)before adsorption and Ctafter adsorption concentration when time t had been passed.To estimate the adsorption capacity of adsorbent the given expression was used,

    whereas,qe(mg·g?1)is the adsorption capacity,Ciand Cf(mg·L?1)are concentrations of Cr(VI)initially and after equilibrium time,V(L)is the volume and M(g)is the mass of adsorbent used.

    3.Results and Discussion

    3.1.Adsorbent characterizations

    3.1.1.Fourier transform infrared spectroscopy(FTIR)

    The presence of functional groups affects the adsorption capacity.Functional group distribution on the surface of adsorbent was analyzed by Fourier transform infrared spectroscopy(FTIR)in the range of 400 to 4000 cm?1.FTIR spectra were shown in Fig.1.Adsorbents had distinct functional groups which appeared on spectra.All the analyzed samples had wide bands at wavenumber from 3000 to 3650 cm?1(with maximum at 3413 cm?1),the absorption bands from 2800 to 3000 cm?1(with maximum at2920 cm?1).The bands centered at 3413 cm?1represented the hydroxyl group stretching vibrations at the surface of adsorbents while the band at 2920 cm?1was assigned to C--H stretching vibrations.The bands around 1735 cm?1were attributed to carbonyl groups(C=O)and peaks at 1630 cm?1occurred at spectra due to stretching of carbonyl groups[49].In the case of HF-SD these bands were strong due to the attachment of carbonyl groups on the surface of sawdust and a kink had also occurred in the range 2300 to 2900 cm?1as compared to H-SD which confirmed that carbonyl groups were present due to aldehyde group attached on the surface of adsorbent after treating it with formaldehyde as compared to H-SD.But in R-SD,these bands also have some comparable intensities due to the presence of lignin which have carbonyl groups in it.The bands at lower wavenumber range 1508 cm?1,1463 cm?1,and 1428 cm?1were associated with lignin present in sawdust[44].The absorption bands at 1508 cm?1were assigned for aromatic skeleton,at 1463 cm?1and 1428 cm?1for C--H deformation.The bands centered at wavenumber 1380 cm?1were attributed to bending of C--H in methane group,while the stretching vibrations of C─O─C linkages were represented with the bands centered at 1163 cm?1and 1115 cm?1and C--O stretching vibrations at 1030 cm?1(polysaccharide peaks in cellulose)[50–52].

    Fig.1.FTIR spectra of R-SD,H-SD and HF-SD.

    3.1.2.Point of zero charge

    In Fig.2 the point which intersects the horizontalline gave a value of pHpzc.The value ofpHpzcfor H-SDand HF-SDwas 4.5 and 5.5 respectively.These points indicated that before intersection with x-axis the adsorbent surface was positively charged while it was negatively charged after pHpzc.The value of pHpzcfor HF-SD was more than H-SD because the surface of adsorbent acquired additional positive charge due to carbonyl groups attached to the surface after the formaldehyde treatment.These additional carbonyl groups enhanced the adsorption capacity of adsorbent by enriching the positive charge over the surface of adsorbent,that developed more electrostatic interaction between active sites of adsorption and anions of chromium metal.The positively charged surface atlower pHwill favor the strong attraction ofadsorbate ions to the adsorbent surface.

    Fig.2.pHpzc for H-SD and HF-SD.

    3.1.3.SEM and EDX studies of an adsorbent

    Fig.3 shows the surface morphology of R-SD and modified materials.Fig.3(a)indicated that the cross-section of R-SD has some distorted pores.Porosity can be seen but it is non-homogeneous and super ficial and bulk of the material seems to have no porosity resulting in less surface area.Filament like pattern of R-SD filled with material imparted anisotropic character.When raw sawdust was treated with HCl,the isotropy of R-SD was eliminated.There was a clear difference in surface morphology after treatment with acid as shown in Fig.3(b).After HCL treatment the surface of sawdust became rough and ruptured due to the acid reaction with components of sawdust,resulting in large cavities produced as this had also been confirmed by BET results for area measurement[53].Fig.3(c)depicts that when H-SD was further function alized with formaldehyde the surface morphology changes again and the cracked and busted surface of sawdust got healed.A glassy film type layer had covered the surface which decreases the surface area of HFSD as confirmed from BET analysis.This transformation indicated that the increase in adsorption of metal for adsorbents was probably due to the surface functional groups instead of porosity of an adsorbent.

    Fig.3.Scanning electron micrographs of(a)R-SD,(b)H-SD and(c)HF-SD.

    The spot analysis of raw sawdust contained major mass percent of carbon and oxygen,with traces of Calcium(Ca),Iron(Fe),Sodium(Na),and Silica(Si).The high percentage of carbon(52.3 wt%)and oxygen(47.4 wt%)in R-SD corresponds to cellulose,hemicellulose and lignin with small quantities of Na(0.02 wt%),Si(0.08 wt%),Ca(0.11 wt%),and Fe(0.06 wt%)which came from the type of soil during growth of plant.When R-SD was treated with acid it ruptured the surface of sawdust and devolved the cavities by removing organic material asa result of which carbon contents were increased and oxygen was decreased and trace elements were completely leached out.Acid treated sawdust had elemental composition as C(59 wt%),O(40.5 wt%),Ca(0.35 wt%),and Fe(0.13 wt%).As acid treated sawdust was further functionalized with formaldehyde carbon and oxygen contents again readjusted with small variations of trace element percentages,due to carbon and oxygen source provided by formaldehyde chemical.The detailed analysis is shown in Table 1.

    Table 1 Chemical composition of raw and modified sawdust

    3.1.4.Thermo-gravimetric analysis

    Fig.4.Thermo-gravimetric(TG)and derivative thermo-gravimetric(DTG)analysis of R-SD,H-SD and HF-SD.

    Thermal stability of raw and prepared adsorbents was studied using thermo-gravimetric analysis.Fig.4 compares the TGA and DTGA curves for R-SD,H-SD and HF-SD.Sawdust is composed of three main components hemicellulose,cellulose and lignin[54].The TGA curves show the initial mass loss of 3%–4%for R-SD,H-SD and HF-SD samples is attributed to the loss of surface adsorbed moisture[55].Thermal oxidation of R-SD,H-SD and HF-SD started at 250°C in the second phase and corresponds to decarboxylation reactions related to degradation of hemicellulose and cellulose up to 350°C,while lignin was burned between 350 °C and 450 °C[56,57].Residual mass of around 3.5%for RSD at the end of experimentation provides the indication of some impurities which were removed by HClsoaking which lead to complete burning of carbonaceous material with no residual mass.This observation is in agreement with EDX results.Further modification of H-SD with formaldehyde provided a protective layer which induces some resistance to thermal heating and avoided complete decomposition of material.

    3.1.5.Surface porosity(BET)analysis

    The surface area of R-SD was 0.04 m2·g?1which rises to 0.54 m2·g?1after soaking the R-SD in an acid solution.The subsequent treatment with formaldehyde solution reduces the surface area to 0.44 m2·g?1.Fig.5 summarizes the pore size distribution of sawdust and its modified forms.The raw material has few pores of meso-and macro-range which were affected by HCl treatment and decreases in numbers in addition to formation of some pores of micro-range.Later,the modification with formaldehyde creates few micropores due to the reaction with carbon structure of sawdust.The formation of new pores in the micro-range is compensated by a decrease in mesopores and resultant decrease in overall surface area of an adsorbent.

    3.2.Adsorption studies

    3.2.1.Effect of adsorbent dose

    It is important to investigate the adsorption of Cr(VI)with different adsorbent doses,as it will illustrate the probability of adsorbent to adsorb chromium metal from solution having a particular initial concentration.As can be seen from Fig.6,with an increase of adsorbent dose the percentage removal of chromium increases for all three types of sawdust.At0.25 g of adsorbent dose,R-SD removes only 36%chromium while the other modified forms of sawdust give 53%and 75%for H-SD and HF-SD,respectively at the same dose as shown in Fig.6.Comparing the removal efficiency of modified sawdust shows that the HF-SD is superior to the H-SD.Probably,it is because of surface functional groups after the formaldehyde functionalization.The upturn of percentage removal was expected because with an increase of adsorbent dose the number of available active sites had been increased for adsorption of heavy metal ions[58].The analysis of experimental observations suggests 4.9 g can be considered as an optimum dose for the chromium removal for R-SD.On the other hand,after activation and functionalization of raw material the optimum dose goes much lower i.e.approximately 3.1 g for H-SD and 2.25 for HF-SD.The optimum dose required for HF-SD is 38%less than raw sawdust.As the adsorption potential of modified forms of sawdust is much higher than the raw sample,only activated and subsequently formaldehyde functionalized sawdust was taken for studying and optimizing the rest of the parameters.Comparing the increases in removal efficiency with small rise in BET surface area after modification of raws a wdust suggests the removal of metal ions is probably due to functional groups on the surface of modified sawdust instead of only pores on the surface.

    Fig.5.BJH(Barrett–Joyner–Halenda)pore size distribution of R-SD,H-SD and HF-SD.

    Fig.6.Effect of adsorbent dose on the percentage removal of chromium(temperature 30 °C,conc.100 mg·L?1,contact time 8 h.pH 6).

    3.2.2.Effect of pH

    pH of the aqueous solution of chromium is another most significant parameter which affects directly the adsorption of metal ions on the surface of an adsorbent.Fig.7 shows that the adsorption capacity decreases greatly with an increase in initial pH from 2 to 10.Chromium uptake was dropped down from 4.1 mg·g?1(at pH=2)to 2.8 mg·g?1for H-SD while it falls from 5.1 mg·g?1to 3.90 mg·g?1for the case of HF-SD,respectively.This may be due the factthatchromium anions behave differently at different pH[59].Hexavalent chromium exists in different forms like H2CrO4,CrO4?2,Cr2O7?2,and HCrO4?1in the aqueous solution.Available chromium exists in H2CrO4form below pH=2 and then chromium anions combined with hydrogen ions to form chromic acid instead of getting adsorbed on surface of adsorbent hence the percentage adsorption of chromium decreased.Above pH=2,chromiumtakes the formofHCrO4?1,Cr2O7?2and atneutral pH or higher it is present in the form of CrO4?2[60].At lower pH chromium adsorbed frequently in the form of HCrO4?1on surface of adsorbent,adsorbent was positively charged and electrostatic forces were developed between adsorbent surface and anions of chromium so,adsorption was high at lower pH.With the increase of pH,HCrO4?1ions of chromium start to deplete and surface charge also shifts from being positive to negative and as a result surface becomes repulsive towards chromium ions.At the same time the competition between OH?and CrO4?2would diminish the uptake to chromium ion onto the adsorbent surface.

    Fig.7.Effect of pH on uptake capacity of chromium at(temperature 30°C,conc.100 mg·L?1,optimum doses,time(300 min)).

    3.2.3.Kinetic studies of Cr(VI)adsorption

    According to Fig.8,it can be observed that at the start,the rate of metal uptake on both(H-SDand HF-SD)was high because of availability of large number of vacant sites.As the time continues available active sites were going to decrease which diminishes the concentration potential and as a result,rate of adsorption slows down.Fig.8 depicts that formaldehyde functionalized adsorbentneeds a long time to reach equilibrium as compared to its counterpart i.e.H-SD.This is probably due to the formation of a thin adhesive layer on the wood fiber after the posttreatment of H-SD with formalin solution resulting in deceleration of water penetration into the structure of HF-SD.It was concluded that 80 min would be sufficient to achieve equilibrium for H-SD while HFSD requires 160 min to reach equilibrium where the maximum uptake was 4.8 mg·g?1.

    The results of the effect of time were used to study the kinetics of adsorption for describing uptake rate of chromium on the modified sawdust.Different kinetic models were used to explore the mechanism involved and the rate controlling step of adsorption.Kinetic models like pseudo- first and second order,and intra-particle diffusion models were used.

    Fig.8.Effect of contact time on uptake capacity of chromium at(temperature 30°C,conc.100 mg·L?1,optimum doses,pH=6).

    Pseudo first order model was offered by Lagergren for a solid/liquid system[7].In most cases The pseudo- first order model does not satisfy the adsorption phenomenon for the entire range of time.The mathematical expression of the pseudo- first order model can be written as Eq.(3).where,qeand qtare the amount of chromiumions which were adsorbed on the unit mass of an adsorbent(mg·g?1)atequilibrium and atany instant of time t(minutes),respectively.k1is the rate constant(min?1)for the pseudo- first order adsorption model.The plot of ln(qe?qt)versus t in Fig.9(a),gives the values of rate constant(k1),and qefrom slope and intercept,respectively.The parameter values of pseudo first-order linear regression for each adsorbent are summarized in Table 2.The values of qe(experimental)were much higher than qe(calculated)for both H-SD and HF-SD and regression coefficient was around 90%for each adsorbent.

    Pseudo-second order model can be expressed as given below in Eq.(4)[59].

    where,k2is the rate constant(g·mg?1·min?1)for the pseudo second order model,qeis the adsorption capacity(mg·g?1)at equilibrium and qtis the instantaneous adsorption capacity(mg·g?1)at any time t(min).t/qtwas plotted against t and value of qewas determined from the slope of plot.The initial adsorption rate k2qe2(g·mg?1·min?1)had been used to determine the rate of adsorption,k2(g·mg?1·min?1)was determined from the intercept of the plot.The regression values are summarized in Table 2.Experimental values of qeare much closer to the calculated value and the regression coefficient is approximately equal to unity for both H-SD and HF-SD.It can be concluded that the best fitted pseudo second order model implies that the mechanism of the Cr(VI)sorption onto modified sawdust is the chemisorption process.

    Fig.9.Kinetic models of adsorption for hexavalent chromium:(a)pseudo first order,(b)pseudo second order,(c)intra-particle diffusion and(d)Boyd mo del.

    Table 2 Kinetic parameters for the adsorption of Cr(VI)onto H-SD and HF-SD

    To get a better insight of the mechanism involved in adsorption phenomenon,it is necessary to determine the rate limiting step.Film diffusion and/or intra-particle diffusion step may limit such type of adsorption processes.The Weber model and Boyd model were tested to infer whether the process is diffusional or film controlling.The Weber model for intra-particle diffusion can be expressed as in the following equation

    Here q is the amount of heavy metal adsorbed(mg·g?1)at time t(min),kidis intra-particle diffusion constant(g·mg?1·min1/2),and c is the intercept;t0.5is the square root of time(min0.5).The experimental data of q was plotted versus t0.5.It was observed that the plot was not a straight line.Instead,it presented a multi-linear behavior which directed about two steps involved in the adsorption phenomenon.Such type of results was obtained by R.Gottipati for adsorption of chromium on Marmelos fruit shell[60].First steep line in beginning of the plot(Fig.8)indicated sudden adsorption on external surfaces of adsorbent particles,while the line at later stage represents the intra-particle diffusion.Initially,abundant avail ability of free sites on the adsorbent surface with high concentration gradient of chromium resulted in the steep curve.At later stage,the slope of curve was small,because most of the external free sites on the adsorbent surface were occupied in the initial step,and metal ions now have to penetrate within the interstitial spaces to get adsorbed.For both,H-SD and HF-SD,initially the rate of adsorption was controlled by film diffusion and then in the second step rate was controlled by inter-and intra-particle diffusion after 10 min for H-SD with a value of q=2.45 mg·g?1and 5 min for HF-SD with a value of q=3.24 mg·g?1until equilibrium was achieved.After equilibrium,no further intra-particle diffusion controlled net adsorption occurred on adsorbent as shown in Fig.9(c).This describes the dual dependence of adsorption rate on film diffusion as well as intra-particle diffusion.Because of the double dependence of adsorption rate,the Boyd model was also applied for analysis of film diffusion kinetic data.The Boyd model for film diffusion can be expressed as,

    F is the fraction of metal amount adsorbed at any time t(min)with respect to amount adsorbed that could be adsorbed at equilibrium,and Bt is a mathematical function of F.Above Eq.(6)had a rearranged form as,

    While,F(xiàn)=q/qe

    where q is the amount of chromium adsorbed(mg·g?1)at any time t(min)and qerepresents the amount of chromium adsorbed(mg·g?1)at equilibrium.

    The Boyd model was plotted by using the values of[?0.4977 ?ln(1?F)]and t(min).It can be seen from Fig.9(d)that the plots were straight lines and did not passed through the origin having intercepts of?0.2446 and?0.6073 for H-SD and HF-SD.Straight lines with intercepts also described that initially the rate of adsorption was controlled by film diffusion;otherwise it was controlled by intra-particle diffusion.The value of diffusion coefficient Diwas calculated by using the given expression.

    where Diis the effective diffusion coefficient(m2·g?1)and r is the radius of adsorbent particles.The values calculated are summarized in Table 2.The value of diffusional coefficient for HF-SDis 2.3 times the diffusional coefficient of H-SD due to the adhesive layer of formalin on the adsorbent.In first 5–10 min approximately 75%of metal had adsorbed for both H-SD and HF-SD here,inter-and intra-particle diffusion was the rate limiting step.

    3.2.4.Adsorption isotherm

    Adsorption isotherms are the most significant to understand interaction of adsorbate with adsorbents,and for the optimum design ofan adsorption system[14].Experimental data were fitted by the Langmuir and Freundlich models.These models can be expressed as follows.

    Langmuir isotherm

    Freundlich isotherm

    Fig.10.Adsorption isotherm models fitting for removal of chromium.

    Here qe(mg·g?1)represents the adsorption capacity and ce(mg·L?1)indicates the concentration of adsorbate in solution at equilibrium,qmis the maximum adsorption capacity,and KLis the Langmuir adsorption equilibrium constant(L·mg?1),KFand n are Freundlich constants related to the adsorption capacity and adsorption intensity of the adsorbents.Fig.10 represents the adsorption isotherms for both adsorbents.The nonlinear fitting parameters from the Langmuir and Freundlich models were listed in Table 3.It could be seen from Fig.10 that the Langmuir model fitted the data better than the Freundlich model,which was confirmed with high value of regression factor(R2)as compared to Freundlich as provided in Table 3.This elaborated that adsorption of chromium was monolayer adsorption.Maximum adsorption capacity(qm)for HF-SD was greater than that for H-SD probably due to the attached carbonyl groups on the surface of adsorbent which leads to enhanced electrostatic force between anions of chromium and adsorbent surface.And due to functional ization with formalin HF-SD needs fewer amounts as compared to H-SD for the same adsorption removal of chromium ions.

    Table 3 Adsorption isotherms fitting parameters for H-SD and HF-SD

    3.2.5.Estimation of thermodynamic parameters

    The Gibbs free energy for the adsorption of Cr(VI)can be calculated by the Van't Hoff equation as provided here,

    The adsorption heat at constant temperature and change in entropy can also be related with the Gibbs free energy by the relationship as,

    By comparing and rearranging Eqs.(11)and(12),we get

    Fig.11.Van't Hoff plot of adsorption equilibrium constant K.

    In the above equations ΔG(J·mol?1)is the Gibbs free energy change,ΔS(J·mol?1·K?1)is the entropy change,ΔH(J·mol?1)is the change in enthalpy,R(8.314 J·mol?1·K?1)is universal gas constant,T(K)is the temperature,K is the constant of equilibrium,CAe(mg·L?1)is concentration of metal on the adsorbent at equilibrium,and Ce(mol·L?1)is the concentration of metal in solution at equilibrium.The plot(Fig.11)of ln K versus 1/T was used to calculate ΔH,ΔS.Negative values of ΔG describe that the adsorption process is spontaneous in nature and as a rule of thumb an increase in negative values of the Gibbs free energy with rise in temperature implies that the lower temperature favors the adsorption of chromium.A positive value of ΔH reflects that the adsorption process is endothermic and the adsorbate has to displace more than one water molecule for their adsorption and this results in the endothermicity of the adsorption process.The positive value ΔS may refer to the affinity of the adsorbent towards the adsorbate and suggests the increase in randomness at the solid–liquid interface with some structural changes in the adsorbate and the adsorbent.The adsorbed water molecules which were displaced by the chromium ions gain more translational entropy than is lost by the adsorbate chromium ions,thus gains more randomness in the system(Table 4).

    Table 4 Thermodynamic parameters for the adsorption of hexavalent chromium onto,H-SD and HF-SD at different temperatures

    In Table 5 comparison of adsorption capacity of A.nilotica sawdust with similar bio-sorbents illustrates that A.nilotica sawdust has almost similar or higher capability to adsorb the chromium metal from waste water.

    Table 5 Comparison of adsorption capacity of Acacia nilotica sawdust with other bio-sorbents

    4.Conclusions

    The research results indicate that the keekar sawdust could serve as an inexpensive source to treat chromium(VI)bearing wastewater.The characterization of synthesized adsorbents by SEM images reveals that the acid treatment of raw sawdust creates few pores on the surface which were later covered by an adhesive layer of formaldehyde and as a consequence surface area decreases.The porosity results through BET analysis support the findings of SEM images.The thermal analysis with TGA reflects that the HCl soaking re fined the raw sawdust and impurities were removed which lead to complete burning of carbonaceous material with no residual mass and elemental analysis by EDX is in agreement with thermal analysis.The treatment of raw sawdust with HCl solution and subsequent modification with formaldehyde improves the adsorption characteristics and as a result the optimum amount required for HF-SD adsorbent is 0.72 times the optimum dose required for H-SD.The kinetic study verifies that the pseudo-second order model is better fitted to experimental observations.The Langmuir isotherm model was followed by both modified forms of sawdust with fairly high regression coefficient(R2=0.995 for H-SD,and R2=0.998 for HF-SD)as compared to the Freundlich model.Thermodynamic analysis depicts the spontaneous and endothermic nature of chromium adsorption onto sawdust based adsorbent.Considering the low price of A.nilotica sawdust,it could be an excellent material for water purification application.

    Acknowledgments

    Authors are thankful to Chemical Engineering Department,University of Engineering and Technology Lahore for providing the facilities for conducting this research study.

    [1]N.K.Srivastava,C.B.Majumder,Novel bio filtration methods for the treatment of heavy metals from industrial wastewater,J.Hazard.Mater.151(2008)1–8.

    [2]N.Tewari,P.Vasudevan,B.K.Guha,Study on biosorption of Cr(VI)by Mucor hiemalis,Biochem.Eng.J.23(2005)185–192.

    [3]F.A.Al-Khaldi,B.Abu-Sharkh,A.M.Abulkibash,M.I.Qureshi,T.Laoui,M.A.Atieh,Effect of acid modification on adsorption of hexavalent chromium(Cr(VI))from aqueous solution by activated carbon and carbon nanotubes,Desalin.Water Treat.3994(2015)1–13.

    [4]S.A.Wanees,A.M.M.Ahmed,M.S.Adam,M.A.Mohamed,Adsorption studies on the removal of hexavalent chromium-contaminated wastewater using activated carbon and bentonite,Asian J.Chem.25(2013)8245–8252.

    [5]F.N.Acar,E.Malkoc,The removal of chromium(VI)from aqueous solutions by Fagus orientalis L,Bioresour.Technol.94(2004)13–15.

    [6]S.Kushwaha,B.Sreedhar,P.P.Sudhakar,A spectroscopic study for understanding the speciation of Cr on palm shell based adsorbents and their application for the remediation of chrome plating effluents,Bioresour.Technol.116(2012)15–23.

    [7]L.Khezami,R.Capart,Removal of chromium(VI)from aqueous solution by activated carbons:kinetic and equilibrium studies,J.Hazard.Mater.123(2005)223–231.

    [8]S.Golbaz,A.J.Jafari,M.Ra fiee,R.R.Kalantary,Separate and simultaneous removal of phenol,chromium,and cyanide from aqueous solution by coagulation/precipitation:mechanisms and theory,Chem.Eng.J.253(2014)251–257.

    [9]Y.Ku,I.L.Jung,Photocatalytic reduction of Cr(VI)in aqueous solutions by UV irradiation with the presence of titanium dioxide,Water Res.35(2001)135–142.

    [10]S.A.Mirbagheri,S.N.Hosseini,Pilot plant investigation on petrochemical wastewater treatment for the removal of copper and chromium with the objective of reuse,Desalination 171(2005)85–93.

    [11]J.S.George,A.Ramos,H.J.Shipley,Tanning facility wastewater treatment:analysis of physical–chemical and reverse osmosis methods,J.Environ.Chem.Eng.3(2015)969–976.

    [12]J.Bohdziewicz,Removal of chromium ions(VI)from underground water in the hybrid complexation–ultra filtration process,Desalination 129(2000)227–235.

    [13]J.Landaburu-Aguirre,E.Pongrácz,P.Per?m?ki,R.L.Keiski,Micellar-enhanced ultrafiltration for the removal of cadmium and zinc:use of response surface methodology to improve understanding of process performance and optimisation,J.Hazard.Mater.180(2010)524–534.

    [14]E.Abbasi-Garravand,C.N.Mulligan,Using micellar enhanced ultra filtration and reduction techniques for removal of Cr(VI)and Cr(III)from water,Sep.Purif.Technol.132(2014)505–512.

    [15]N.H.Mthombeni,M.S.Onyango,O.Aoyi,Adsorption of hexavalent chromium onto magnetic natural zeolite-polymer composite,J.Taiwan Inst.Chem.Eng.50(2015)242–251.

    [16]B.Alyüz,S.Veli,Kinetics and equilibrium studies for the removal of nickel and zinc from aqueous solutions by ion exchange resins,J.Hazard.Mater.167(2009)482–488.

    [17]S.Rengaraj,K.-H.Yeon,S.-H.Moon,Removal of chromium from water and wastewater by ion exchange resins,J.Hazard.Mater.87(2001)273–287.

    [18]D.Sbizzaro Filho,G.B.Bota,R.B.Borri,F.J.C.Teran,Electrocoagulation/ flotation followed by fluidized bed anaerobic reactor applied to tannery effluent treatment,Desalin.Water Treat.37(2012)359–363.

    [19]D.Sivakumar,Experimental and analytical model studies on leachate volume computation from solid waste,Int.J.Environ.Sci.Technol.10(2013)903–916.

    [20]K.Ranganathan,Chromium removal by activated carbons prepared from Casuarina equisetifolia leaves,Bioresour.Technol.73(2000)99–103.

    [21]J.-L.Gong,B.Wang,G.-M.Zeng,C.-P.Yang,C.-G.Niu,Q.-Y.Niu,W.-J.Zhou,Y.Liang,Removal of cationic dyes from aqueous solution using magnetic multi-wall carbon nanotube nanocomposite as adsorbent,Chem.Eng.J.164(2009)1517–1522.

    [22]P.Xu,G.M.Zeng,D.L.Huang,C.L.Feng,S.Hu,M.H.Zhao,C.Lai,Z.Wei,C.Huang,G.X.Xie,Z.F.Liu,Use of iron oxide nanomaterials in wastewater treatment:a review,Sci.Total Environ.424(2012)1–10.

    [23]A.Jusoh,L.Su Shiung,N.Ali,M.J.M.M.Noor,Asimulation study of the removal efficiency of granular activated carbon on cadmium and lead,Desalination 206(2007)9–16.

    [24]N.R.Bishnoi,M.Bajaj,N.Sharma,A.Gupta,Adsorption of Cr(VI)on activated rice husk carbon and activated alumina,Bioresour.Technol.91(2004)305–307.

    [25]V.S.Mane,P.V.V.Babu,Studies on the adsorption of Brilliant Green dye from aqueous solution onto low-cost NaOH treated saw dust,Desalination 273(2011)321–329.

    [26]K.Selvaraj,S.Manonmani,S.Pattabhi,Removal of hexavalent chromium using distillery sludge,Bioresour.Technol.89(2003)207–211.

    [27]V.Murphy,H.Hughes,P.McLoughlin,Comparative study of chromium biosorption by red,green and brown seaweed biomass,Chemosphere 70(2008)1128–1134.

    [28]Y.?ahin,A.?ztürk,Biosorption of chromium(VI)ions from aqueous solution by the bacterium Bacillus thuringiensis,Process Biochem.40(2005)1895–1901.

    [29]V.K.Gupta,A.Rastogi,A.Nayak,Adsorption studies on the removal of hexavalent chromium from aqueous solution using a low cost fertilizer industry waste material,J.Colloid Interface Sci.342(2010)135–141.

    [30]Y.Yun,D.Park,J.M.Park,B.Volesky,Biosorption of trivalent chromium on the brown seaweed biomass biosorption of trivalent chromium on the brown seaweed biomass,Environ.Sci.Technol.35(2001)4353–4358.

    [31]M.Bansal,D.Singh,V.K.Garg,A comparative study for the removal of hexavalent chromium from aqueous solution by agriculture wastes'carbons,J.Hazard.Mater.171(2009)83–92.

    [32]A.M.Zaheer,N.Ramzan,S.Naveed,N.Feroze,Ni(II)removal by biosorption using Ficus religiosa(Peepal)leaves,J.Chil.Chem.Soc.1(2010)81–84.

    [33]S.P.Dubey,K.Gopal,Adsorption of chromium(VI)on low cost adsorbents derived from agricultural waste material:a comparative study,J.Hazard.Mater.145(2007)465–470.

    [34]D.Mohan,S.Rajput,V.K.Singh,P.H.Steele,C.U.Pittman,Modeling and evaluation of chromium remediation from water using low cost bio-char,a green adsorbent,J.Hazard.Mater.188(2011)319–333.

    [35]S.Saraswat,J.P.N.Rai,Heavy metal adsorption from aqueous solution using Eichhornia crassipes dead biomass,Int.J.Miner.Process.94(2010)203–206.

    [36]H.Gao,Y.Liu,G.Zeng,W.Xu,T.Li,W.Xia,Characterization of Cr(VI)removal from aqueous solutions by a surplus agricultural waste—rice straw,J.Hazard.Mater.150(2008)446–452.

    [37]B.Tadesse,E.Teju,N.Megersa,The Teff straw:a novel low-cost adsorbent for quantitative removal of Cr(VI)from contaminated aqueous samples,Desalin.Water Treat.56(2014)1–12.

    [38]U.K.Garg,M.P.Kaur,V.K.Garg,D.Sud,Removal of hexavalent chromium from aqueous solution by agricultural waste biomass,J.Hazard.Mater.140(2007)60–68.

    [39]S.Schiewer,S.B.Patil,Modeling the effect of pH on biosorption of heavy metals by citrus peels,J.Hazard.Mater.157(2008)8–17.

    [40]M.Cobas,M.A.Sanromán,M.Pazos,Box–Behnken methodology for Cr(VI)and leather dyes removal by an eco-friendly biosorbent:F.vesiculosus,Bioresour.Technol.160(2014)166–174.

    [41]L.S.Oliveira,A.S.Franca,T.M.Alves,S.D.F.Rocha,Evaluation of untreated coffee husks as potential biosorbents for treatment of dye contaminated waters,J.Hazard.Mater.155(2008)507–512.

    [42]A.Bhatnagar,E.Kumar,A.K.Minocha,B.H.Jeon,H.Song,Y.C.Seo,Removal ofanionic dyes from water using Citrus limonum(lemon)peel:equilibrium studies and kinetic modeling,Sep.Sci.Technol.44(2009)316–334.

    [43]A.Saeed,M.Iqbal,M.W.Akhtar,Removal and recovery of lead(II)from single and multimetal(Cd,Cu,Ni,Zn)solutions by crop milling waste(black gram husk),J.Hazard.Mater.117(2005)65–73.

    [44]L.Agouborde,R.Navia,Heavy metals retention capacity of a non-conventional sorbent developed from a mixture of industrial and agricultural wastes,J.Hazard.Mater.167(2009)536–544.

    [45]A.Verma,S.Chakraborty,J.K.Basu,Adsorption study of hexavalent chromium using tamarind hull-based adsorbents,Sep.Purif.Technol.50(2006)336–341.

    [46]P.Venkateswarlu,M.V.Ratnam,D.S.Rao,M.V.Rao,Removal of chromium from an aqueous solution using Azadirachta indica(neem)leaf powder as an adsorbent,Int.J.Phys.Sci.2(2007)188–195.

    [47]V.Sarin,K.K.Pant,Removal of chromium from industrial waste by using eucalyptus bark,Bioresour.Technol.97(2006)15–20.

    [48]N.Rani,A.Gupta,A.K.Yadav,Removal of Cr(VI)from aqueous solutions by Acacia nilotica bark,Environ.Technol.27(2006)597–602.

    [49]S.Rangabhashiyam,N.Selvaraju,Efficacy of unmodified and chemically modified Swietenia mahagoni shells for the removal of hexavalent chromium from simulated wastewater,J.Mol.Liq.209(2015)487–497.

    [50]A.B.Wassie,V.C.Srivastava,Teff straw characterization and utilization for chromium removal from wastewater:kinetics,isotherm and thermodynamic modelling,J.Environ.Chem.Eng.4(2016)1117–1125.

    [51]J.B.Dima,C.Sequeiros,N.E.Zaritzky,Hexavalent chromium removal in contaminated water using reticulated chitosan micro/nanoparticles from seafood processing wastes,Chemosphere 141(2015)100–111.

    [52]M.H.Dehghani,D.Sanaei,I.Ali,A.Bhatnagar,Removal of chromium(VI)from aqueous solution using treated waste newspaper as a low-cost adsorbent:kinetic modeling and isotherm studies,J.Mol.Liq.215(2016)671–679.

    [53]M.H.Kalavathy,L.R.Miranda,Comparison of copper adsorption from aqueous solution using modified and unmodified Hevea brasiliensis saw dust,Desalination 255(2010)165–174.

    [54]S.I.Yang,M.S.Wu,C.Y.Wu,Application of biomass fast pyrolysis part I:pyrolysis characteristics and products,Energy 66(2014)162–171.

    [55]M.M.A.El-Latif,A.M.Ibrahim,Adsorption,kinetic and equilibrium studies on removal of basic dye from aqueous solutions using hydrolyzed Oak sawdust,Desalin.Water Treat.6(2009)252–268.

    [56]S.Ren,H.Lei,L.Wang,Q.Bu,S.Chen,J.Wu,Thermal behaviour and kinetic study for woody biomass torrefaction and torre fied biomass pyrolysis by TGA,Biosyst.Eng.116(2013)420–426.

    [57]H.Li,X.Liu,R.Legros,X.T.Bi,C.J.Lim,S.Sokhansanj,Torrefaction of sawdust in a fluidized bed reactor,Bioresour.Technol.103(2012)453–458.

    [58]V.K.Garg,R.Gupta,R.Kumar,R.K.Gupta,Adsorption of chromium from aqueous solution on treated sawdust,Bioresour.Technol.92(2004)79–81.

    [59]T.Karthikeyan,S.Rajgopal,L.R.Miranda,Chromium(VI)adsorption from aqueous solution by Hevea brasilinesis sawdust activated carbon,J.Hazard.Mater.124(2005)192–199.

    [60]R.Gottipati,S.Mishra,Preparation of microporous activated carbon from Aegle marmelos fruit shell by KOH activation,J.Ind.Eng.Chem.36(2016)355–363.

    国产成人欧美| 亚洲欧洲日产国产| 成在线人永久免费视频| 成年av动漫网址| 日韩免费高清中文字幕av| 黑人猛操日本美女一级片| 99久久国产精品久久久| 满18在线观看网站| 日韩视频一区二区在线观看| 下体分泌物呈黄色| 午夜福利在线观看吧| 国产淫语在线视频| 悠悠久久av| 成年人午夜在线观看视频| 在线观看免费高清a一片| 亚洲av男天堂| h视频一区二区三区| 男女高潮啪啪啪动态图| 国产不卡av网站在线观看| 热99国产精品久久久久久7| a级毛片在线看网站| 一区二区三区激情视频| 亚洲精品乱久久久久久| 美女国产高潮福利片在线看| 精品人妻熟女毛片av久久网站| 丝袜美足系列| 最新的欧美精品一区二区| 欧美日韩黄片免| 777米奇影视久久| 精品少妇黑人巨大在线播放| 精品第一国产精品| 男女之事视频高清在线观看| 老司机深夜福利视频在线观看 | 亚洲一码二码三码区别大吗| 午夜免费观看性视频| 悠悠久久av| 亚洲五月色婷婷综合| 一本大道久久a久久精品| 国产精品亚洲av一区麻豆| 亚洲天堂av无毛| 亚洲精品日韩在线中文字幕| 深夜精品福利| 欧美97在线视频| 99久久人妻综合| 中文字幕制服av| 亚洲av成人一区二区三| 日本av手机在线免费观看| 国产麻豆69| 欧美日韩视频精品一区| 肉色欧美久久久久久久蜜桃| 性高湖久久久久久久久免费观看| 999精品在线视频| 中文字幕精品免费在线观看视频| 国产精品一区二区免费欧美 | 老熟妇乱子伦视频在线观看 | 国产精品一二三区在线看| av一本久久久久| 伊人亚洲综合成人网| 18在线观看网站| 国产成人精品在线电影| 国产在线一区二区三区精| 中文字幕制服av| 亚洲欧美激情在线| 我要看黄色一级片免费的| 亚洲 欧美一区二区三区| 国产高清国产精品国产三级| 欧美日韩亚洲国产一区二区在线观看 | 久久精品成人免费网站| 狠狠婷婷综合久久久久久88av| 亚洲国产成人一精品久久久| 高清欧美精品videossex| 久久国产亚洲av麻豆专区| 男女高潮啪啪啪动态图| 日本wwww免费看| 大香蕉久久成人网| 亚洲精品久久久久久婷婷小说| 国产精品久久久久久人妻精品电影 | 欧美成狂野欧美在线观看| 在线永久观看黄色视频| 国产激情久久老熟女| 久久亚洲国产成人精品v| 久久国产亚洲av麻豆专区| 日韩人妻精品一区2区三区| videos熟女内射| 啦啦啦在线免费观看视频4| 亚洲精品日韩在线中文字幕| 日韩欧美一区视频在线观看| 日日摸夜夜添夜夜添小说| 后天国语完整版免费观看| 日韩人妻精品一区2区三区| 男女边摸边吃奶| 成年美女黄网站色视频大全免费| 国产福利在线免费观看视频| 啦啦啦免费观看视频1| 丝袜在线中文字幕| 国产精品免费视频内射| 精品少妇内射三级| 欧美国产精品一级二级三级| 老司机福利观看| 欧美一级毛片孕妇| 多毛熟女@视频| 欧美日韩一级在线毛片| 午夜老司机福利片| 另类亚洲欧美激情| 男人操女人黄网站| 亚洲天堂av无毛| 国产福利在线免费观看视频| 久久天躁狠狠躁夜夜2o2o| 成人手机av| 这个男人来自地球电影免费观看| 电影成人av| 丝袜人妻中文字幕| 亚洲美女黄色视频免费看| 狠狠婷婷综合久久久久久88av| 精品免费久久久久久久清纯 | 日韩中文字幕视频在线看片| 爱豆传媒免费全集在线观看| 婷婷丁香在线五月| 777米奇影视久久| 日韩大片免费观看网站| av片东京热男人的天堂| 日韩制服骚丝袜av| 各种免费的搞黄视频| 久久久久精品国产欧美久久久 | 法律面前人人平等表现在哪些方面 | 国产国语露脸激情在线看| 可以免费在线观看a视频的电影网站| 99热网站在线观看| 亚洲国产av影院在线观看| 成年女人毛片免费观看观看9 | 777米奇影视久久| 一区二区av电影网| 国产成人免费观看mmmm| 亚洲成国产人片在线观看| 亚洲成国产人片在线观看| √禁漫天堂资源中文www| 亚洲免费av在线视频| 国产极品粉嫩免费观看在线| 日本av免费视频播放| 两性夫妻黄色片| 午夜精品久久久久久毛片777| 亚洲人成电影免费在线| 99久久精品国产亚洲精品| 夜夜夜夜夜久久久久| 亚洲精品乱久久久久久| 高清av免费在线| 亚洲国产中文字幕在线视频| 狂野欧美激情性xxxx| 免费av中文字幕在线| 视频区图区小说| 国产精品偷伦视频观看了| 一区在线观看完整版| 亚洲免费av在线视频| 中文精品一卡2卡3卡4更新| 欧美人与性动交α欧美精品济南到| 一级片免费观看大全| a在线观看视频网站| 亚洲熟女精品中文字幕| 99香蕉大伊视频| av又黄又爽大尺度在线免费看| 国产片内射在线| 亚洲欧美日韩高清在线视频 | bbb黄色大片| 桃红色精品国产亚洲av| 日韩人妻精品一区2区三区| 久久性视频一级片| 国产熟女午夜一区二区三区| 欧美精品一区二区免费开放| 最新的欧美精品一区二区| 麻豆国产av国片精品| 亚洲国产av影院在线观看| 国产精品免费大片| 欧美中文综合在线视频| 国产高清videossex| 日日夜夜操网爽| 欧美日韩国产mv在线观看视频| 色婷婷久久久亚洲欧美| 日韩 亚洲 欧美在线| 亚洲国产成人一精品久久久| 欧美97在线视频| 深夜精品福利| 久久国产精品大桥未久av| 麻豆av在线久日| 久久久国产成人免费| 欧美一级毛片孕妇| 国产精品一区二区在线不卡| 亚洲精品国产色婷婷电影| 夫妻午夜视频| 国产亚洲午夜精品一区二区久久| 国产精品一区二区免费欧美 | videos熟女内射| 国产真人三级小视频在线观看| 宅男免费午夜| 最近中文字幕2019免费版| 国产精品自产拍在线观看55亚洲 | 久久久久久久久免费视频了| 国产精品欧美亚洲77777| 黑人欧美特级aaaaaa片| 看免费av毛片| 女性被躁到高潮视频| 欧美日韩av久久| 国产精品一二三区在线看| 18禁观看日本| 巨乳人妻的诱惑在线观看| 性高湖久久久久久久久免费观看| 日韩人妻精品一区2区三区| 国产区一区二久久| 久久国产精品人妻蜜桃| 精品久久久精品久久久| 精品少妇一区二区三区视频日本电影| 另类精品久久| 精品人妻一区二区三区麻豆| kizo精华| 大陆偷拍与自拍| 桃红色精品国产亚洲av| 自拍欧美九色日韩亚洲蝌蚪91| 国产成人啪精品午夜网站| 啦啦啦在线免费观看视频4| av不卡在线播放| 97人妻天天添夜夜摸| 国产极品粉嫩免费观看在线| 一本色道久久久久久精品综合| 9色porny在线观看| 下体分泌物呈黄色| 老熟女久久久| 亚洲男人天堂网一区| 啪啪无遮挡十八禁网站| 香蕉国产在线看| 菩萨蛮人人尽说江南好唐韦庄| 男女国产视频网站| 免费在线观看视频国产中文字幕亚洲 | 国产欧美日韩一区二区精品| 视频区图区小说| 色综合欧美亚洲国产小说| 国产成人免费无遮挡视频| 国产精品免费视频内射| 少妇粗大呻吟视频| 国产精品久久久av美女十八| 十八禁网站免费在线| av超薄肉色丝袜交足视频| 桃红色精品国产亚洲av| 久久久久视频综合| 啦啦啦中文免费视频观看日本| 亚洲午夜精品一区,二区,三区| 亚洲全国av大片| kizo精华| 99国产精品一区二区蜜桃av | 老司机在亚洲福利影院| 999精品在线视频| 9191精品国产免费久久| 啦啦啦在线免费观看视频4| 下体分泌物呈黄色| 国产亚洲av片在线观看秒播厂| 乱人伦中国视频| av欧美777| av在线app专区| 精品视频人人做人人爽| 亚洲av日韩精品久久久久久密| 电影成人av| 两人在一起打扑克的视频| 亚洲欧美日韩高清在线视频 | 熟女少妇亚洲综合色aaa.| 淫妇啪啪啪对白视频 | 亚洲五月色婷婷综合| 每晚都被弄得嗷嗷叫到高潮| 美女中出高潮动态图| 99精国产麻豆久久婷婷| 在线观看免费午夜福利视频| 99国产精品99久久久久| 国产精品99久久99久久久不卡| 国产深夜福利视频在线观看| 久久久欧美国产精品| 国内毛片毛片毛片毛片毛片| 波多野结衣一区麻豆| 亚洲精品国产色婷婷电影| 亚洲一区中文字幕在线| 国产亚洲av片在线观看秒播厂| 老汉色av国产亚洲站长工具| 一级片免费观看大全| 欧美另类亚洲清纯唯美| 不卡av一区二区三区| videos熟女内射| 国产精品.久久久| 精品熟女少妇八av免费久了| 大片电影免费在线观看免费| 日日摸夜夜添夜夜添小说| 亚洲熟女毛片儿| 日韩欧美一区视频在线观看| 亚洲精品日韩在线中文字幕| 欧美少妇被猛烈插入视频| 国产在视频线精品| 老熟女久久久| 性高湖久久久久久久久免费观看| 91国产中文字幕| 法律面前人人平等表现在哪些方面 | 久久久国产一区二区| 精品卡一卡二卡四卡免费| 每晚都被弄得嗷嗷叫到高潮| 久久av网站| 十八禁人妻一区二区| 黄片播放在线免费| 丝袜脚勾引网站| 丝袜喷水一区| 国产精品.久久久| 日本av免费视频播放| 免费av中文字幕在线| 精品福利观看| 天天影视国产精品| 天天操日日干夜夜撸| 国产区一区二久久| 精品久久久精品久久久| 97在线人人人人妻| 国产精品av久久久久免费| 一进一出抽搐动态| 久久国产精品人妻蜜桃| 精品少妇一区二区三区视频日本电影| 国产精品久久久人人做人人爽| 69av精品久久久久久 | 少妇人妻久久综合中文| av国产精品久久久久影院| 国产精品一区二区精品视频观看| 国产精品麻豆人妻色哟哟久久| 99精国产麻豆久久婷婷| 国产麻豆69| 91精品伊人久久大香线蕉| 老司机靠b影院| 国产在线视频一区二区| 免费人妻精品一区二区三区视频| 亚洲欧洲日产国产| 午夜激情av网站| 高清在线国产一区| 午夜激情av网站| 80岁老熟妇乱子伦牲交| 亚洲七黄色美女视频| 99热网站在线观看| 狂野欧美激情性bbbbbb| 亚洲久久久国产精品| 性少妇av在线| 咕卡用的链子| 精品人妻一区二区三区麻豆| 免费高清在线观看日韩| 一进一出抽搐动态| 少妇猛男粗大的猛烈进出视频| 午夜精品国产一区二区电影| 日韩熟女老妇一区二区性免费视频| 免费不卡黄色视频| 99九九在线精品视频| 91成人精品电影| 美女脱内裤让男人舔精品视频| 精品一区二区三区四区五区乱码| 啦啦啦在线免费观看视频4| 999精品在线视频| 国产三级黄色录像| 成人国产av品久久久| 在线av久久热| www日本在线高清视频| 亚洲精品久久成人aⅴ小说| 久久久久久久久久久久大奶| 多毛熟女@视频| 老司机靠b影院| 国产精品一区二区在线观看99| 国产成人一区二区三区免费视频网站| 咕卡用的链子| 国产在线观看jvid| 天堂中文最新版在线下载| 一本色道久久久久久精品综合| 麻豆乱淫一区二区| 欧美黄色片欧美黄色片| 久久性视频一级片| 色视频在线一区二区三区| 亚洲专区中文字幕在线| 国产一区二区激情短视频 | 麻豆乱淫一区二区| 亚洲精品粉嫩美女一区| 黄色片一级片一级黄色片| 又黄又粗又硬又大视频| 美女午夜性视频免费| 美女大奶头黄色视频| 9色porny在线观看| 久久久久久久大尺度免费视频| 日韩熟女老妇一区二区性免费视频| 老司机午夜福利在线观看视频 | 精品国产乱码久久久久久小说| 亚洲av片天天在线观看| 丝袜喷水一区| 欧美日韩成人在线一区二区| 无限看片的www在线观看| 国产欧美日韩一区二区三 | 一级毛片精品| 国产精品麻豆人妻色哟哟久久| 大香蕉久久成人网| 欧美日韩国产mv在线观看视频| 亚洲精品一二三| 菩萨蛮人人尽说江南好唐韦庄| 国产精品一二三区在线看| 精品国产乱码久久久久久男人| 黄片播放在线免费| 精品国产一区二区三区四区第35| 十八禁人妻一区二区| 91九色精品人成在线观看| 免费在线观看完整版高清| av片东京热男人的天堂| 在线十欧美十亚洲十日本专区| 亚洲成人免费av在线播放| 国产精品av久久久久免费| 人人妻人人澡人人看| 久久99热这里只频精品6学生| 亚洲午夜精品一区,二区,三区| 亚洲精品一区蜜桃| 国产精品二区激情视频| 国产成人影院久久av| 国产人伦9x9x在线观看| 侵犯人妻中文字幕一二三四区| 国产免费av片在线观看野外av| 亚洲av日韩在线播放| 91麻豆av在线| 老熟妇乱子伦视频在线观看 | 日韩欧美一区视频在线观看| 亚洲七黄色美女视频| 多毛熟女@视频| 亚洲专区国产一区二区| 精品人妻在线不人妻| 汤姆久久久久久久影院中文字幕| 高清在线国产一区| 狂野欧美激情性bbbbbb| 亚洲国产精品999| tocl精华| 99国产精品99久久久久| 伊人久久大香线蕉亚洲五| 午夜福利在线免费观看网站| 极品人妻少妇av视频| 丰满少妇做爰视频| 国产精品麻豆人妻色哟哟久久| 亚洲人成电影观看| 一本综合久久免费| kizo精华| 爱豆传媒免费全集在线观看| 亚洲三区欧美一区| 捣出白浆h1v1| 99久久精品国产亚洲精品| 美女国产高潮福利片在线看| 宅男免费午夜| 国产区一区二久久| a级毛片黄视频| 十八禁网站网址无遮挡| 日本撒尿小便嘘嘘汇集6| 日韩免费高清中文字幕av| 国产淫语在线视频| 黄片大片在线免费观看| 免费不卡黄色视频| 久久精品国产综合久久久| 欧美黄色片欧美黄色片| 免费不卡黄色视频| 国产精品99久久99久久久不卡| 久久久久精品人妻al黑| 极品人妻少妇av视频| 水蜜桃什么品种好| 国产亚洲欧美在线一区二区| 老司机福利观看| 丝袜在线中文字幕| 国产亚洲精品久久久久5区| 久久ye,这里只有精品| av在线app专区| 欧美日韩福利视频一区二区| 国产精品免费视频内射| 免费观看人在逋| 久久国产精品人妻蜜桃| videos熟女内射| 亚洲av国产av综合av卡| 国产1区2区3区精品| 午夜免费观看性视频| 日本wwww免费看| 久久国产精品影院| 黑丝袜美女国产一区| 制服人妻中文乱码| 国产日韩欧美亚洲二区| 汤姆久久久久久久影院中文字幕| 亚洲av片天天在线观看| 热99久久久久精品小说推荐| 色婷婷av一区二区三区视频| 国产视频一区二区在线看| 99re6热这里在线精品视频| 亚洲国产欧美在线一区| 黄色 视频免费看| 亚洲精品中文字幕一二三四区 | 黄色片一级片一级黄色片| 日韩 亚洲 欧美在线| 99国产综合亚洲精品| 久久热在线av| 亚洲国产看品久久| 日本wwww免费看| 国产精品一二三区在线看| 狠狠狠狠99中文字幕| 午夜福利视频精品| 亚洲第一青青草原| 亚洲精品国产一区二区精华液| 97在线人人人人妻| 免费观看a级毛片全部| 九色亚洲精品在线播放| 黄色视频,在线免费观看| xxxhd国产人妻xxx| 国产一卡二卡三卡精品| 国产男人的电影天堂91| 丝瓜视频免费看黄片| 亚洲色图 男人天堂 中文字幕| 国产精品99久久99久久久不卡| 国产亚洲精品久久久久5区| 久久国产亚洲av麻豆专区| 亚洲精品美女久久久久99蜜臀| 高清在线国产一区| 在线观看一区二区三区激情| 国产精品秋霞免费鲁丝片| 亚洲成人免费电影在线观看| 交换朋友夫妻互换小说| 99九九在线精品视频| 视频区欧美日本亚洲| av国产精品久久久久影院| 日韩欧美一区二区三区在线观看 | 久久精品国产综合久久久| 精品人妻熟女毛片av久久网站| www.熟女人妻精品国产| 在线观看一区二区三区激情| 国产成人av教育| 欧美亚洲日本最大视频资源| 中国美女看黄片| 欧美激情久久久久久爽电影 | 热99国产精品久久久久久7| 性高湖久久久久久久久免费观看| 色视频在线一区二区三区| 宅男免费午夜| 欧美国产精品一级二级三级| 精品视频人人做人人爽| 狂野欧美激情性xxxx| 涩涩av久久男人的天堂| 在线观看免费午夜福利视频| 性色av一级| 午夜福利在线观看吧| 欧美精品啪啪一区二区三区 | 欧美日韩一级在线毛片| 欧美日韩黄片免| 久久九九热精品免费| 亚洲欧美日韩另类电影网站| 亚洲av欧美aⅴ国产| 欧美日韩av久久| 国产成人精品无人区| 亚洲欧美一区二区三区久久| 欧美激情极品国产一区二区三区| 黄色a级毛片大全视频| 欧美激情 高清一区二区三区| 黑人猛操日本美女一级片| 青春草视频在线免费观看| 亚洲精品乱久久久久久| 女人高潮潮喷娇喘18禁视频| 精品熟女少妇八av免费久了| 精品国产超薄肉色丝袜足j| 国产黄色免费在线视频| 777久久人妻少妇嫩草av网站| 国产成人精品无人区| 女性被躁到高潮视频| 性色av乱码一区二区三区2| 亚洲精品国产av成人精品| 国产一区二区三区在线臀色熟女 | 91九色精品人成在线观看| 国产一区二区三区在线臀色熟女 | 9191精品国产免费久久| 91成年电影在线观看| 美女福利国产在线| av有码第一页| 99国产精品99久久久久| 欧美激情极品国产一区二区三区| 亚洲精品久久午夜乱码| 一本一本久久a久久精品综合妖精| 国产精品 欧美亚洲| 精品一区二区三卡| 成人国语在线视频| av线在线观看网站| 两人在一起打扑克的视频| 日本91视频免费播放| 各种免费的搞黄视频| 日本91视频免费播放| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕制服av| 人人妻人人澡人人看| 悠悠久久av| 精品国产乱码久久久久久小说| 在线亚洲精品国产二区图片欧美| 午夜老司机福利片| 中亚洲国语对白在线视频| 亚洲一区中文字幕在线| 美女大奶头黄色视频| av免费在线观看网站| 日韩人妻精品一区2区三区| 午夜激情av网站| 1024视频免费在线观看| 国产精品一二三区在线看| 后天国语完整版免费观看| 飞空精品影院首页| 三级毛片av免费| 老司机靠b影院| 日韩人妻精品一区2区三区| 免费在线观看视频国产中文字幕亚洲 | 9热在线视频观看99| 国产成人精品无人区| 亚洲精品成人av观看孕妇| 超色免费av| 69精品国产乱码久久久| 亚洲精品国产av蜜桃| 精品久久久久久电影网| 精品亚洲成国产av| 2018国产大陆天天弄谢| 久久天堂一区二区三区四区| 久久久久久久久免费视频了| 国产欧美日韩精品亚洲av| 丁香六月欧美| 亚洲精品国产精品久久久不卡|