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

    Development of low power non-thermal plasma jet and optimization of operational parameters for treating dyes and emerging contaminants

    2022-09-06 13:04:38DeepchandraJOSHIVedaPRAKASHShaikhZiauddinAHAMMADSatyanandaKARandSREEKRISHNAN
    Plasma Science and Technology 2022年10期

    Deepchandra JOSHI,G Veda PRAKASH,Shaikh Ziauddin AHAMMAD,?,Satyananda KAR,?and T R SREEKRISHNAN

    1 Department of Biochemical Engineering and Biotechnology,Indian Institute of Technology Delhi,New Delhi 110016,India

    2 Department of Energy Science and Engineering,Indian Institute of Technology Delhi,New Delhi 110016,India

    Abstract Non-thermal plasma has emerged as an effective treatment system against the latest class of highly recalcitrant and toxic environmental pollutants termed emerging contaminants(ECs).In the present work,a detailed experimental study is carried out to evaluate the efficacy of a non-thermal plasma jet with two dyes,Rd.B and Met.Blue,as model contaminants.The plasma jet provided a complete dye decoloration in 30 min with an applied voltage of 6.5 kVp?p.·OH,having the highest oxidation potential,acts as the main reactive species,which with direct action on contaminants also acts indirectly by getting converted into H2O2 and O3.Further,the effect of critical operational parameters viz,sample pH,applied voltage(4.5-6.5 kVp?p),conductivity(5-20 mS cm?1),and sample distance on plasma treatment efficacy was also examined.Out of all the assessed parameters,the applied voltage and sample conductivity was found to be the most significant operating parameters.A high voltage and low conductivity favored the dye decoloration,while the pH effect was not that significant.To understand the influence of plasma discharge gas on treatment efficacy,all the experiments are conducted with argon and helium gases under the fixed geometrical configuration.Both the gases provided a similar dye decoloration efficiency.The DBD plasma system with complete dye removal also rendered maximum mineralization of 73%for Rd.B,and 60%for Met.Blue.Finally,the system’s efficiency against the actual ECs (four pharmaceutical compounds,viz,metformin,atenolol,acetaminophen,and ranitidine) and microbial contaminant (E.coli) was also tested.The system showed effectivity in the complete removal of targeted pharmaceuticals and a log2.5 E.coli reduction.The present systematic characterization of dye degradation could be of interest to large communities working towards commercializing plasma treatment systems.

    Keywords: advanced oxidation processes (AOPs),non-thermal plasma,dye-degradation,ECs treatment,reactive species quantification,mineralization

    1.Introduction

    Emerging contaminants (ECs) are the latest class of environmental contaminants that have been recently detected in the environment.The reason for their increased environmental concentration is mainly due to the increased use of chemicals and unresponsible anthropogenic activities [1].The ECs mainly found in urban wastewater include pharmaceuticals,personal care compounds,endocrine disrupting compounds,pesticides,and antibiotic-resistant bacteria (ARB) [2].These contaminants are generally found in a very narrow concentration range of a few ng l?1.However,in such low concentrations also,they are extremely harmful to human and aquatic health [3].The conventional wastewater treatment systems(WWTs) were not designed considering ECs removal;therefore,they show limited efficiency against them [2].This necessitates advanced treatment systems,which are highly effective in taking care of ECs.The membrane filtration,adsorption process,nanotechnology,and chemical oxidation process are effective advanced processes.However,they have a major limitation that either they are a phase transfer method or require toxic chemicals [4].These issues limit their commercial use against the ECs.The advanced oxidation process(AOPs),which are based on the high oxidation potential of·OH radical,has emerged as a highly effective treatment process [5-8].Moreover,like other advanced treatment technologies,AOPs do not have major limitations when used on a large scale [9].The·OH is highly reactive and has the highest oxidation potential among the common disinfectants: ·OH(E°= 2.80 V) > O3(E° = 2.07 V) > H2O2(E°= 1.8 V) > Cl2(E°= 1.36 V).Therefore the ·OH can oxidize and degrade organic contaminants most easily and in a non-specific manner.The conventional AOPs mainly include Fenton and Fenton-like reactions,UV/H2O2/O3based processes,and photocatalysis[9].Recently cold plasma has been found to be an effective AOPs and is already under extensive laboratory trial for finding its large-scale applications in medicine,industrial,and wastewater treatment [10,11].Plasma is the mixture of charged electronic species that could be used to produce various reactive oxygen and nitrogen species(RONS)such as·OH,NO,·O2?,O3,N2+,NO3?,NO2?etc.In addition,it also provides other action mechanisms such as UV radiation,heat energy,sonication,etc [12-14].All these multiple action mechanisms act synergistically and enhance plasma treatment efficiency.In addition,the plasma operation is also energy efficient,easy to operate,and eco-friendly.Multiple types of plasma systems are reported in the literature,such as pulsed discharge plasma reactor,gliding arc discharge reactor,contact glow discharge electrolysis plasma reactor,atmospheric pressure glow discharge (APGD) plasma reactor,and dielectric barrier discharge(DBD)reactor,etc[15-19],each of these has their advantages and setbacks [20].Although several studies have been reported on dye treatment with non-thermal plasma,the systematic characterization in view of treating ECs is scarce.Thus,more systematic studies are needed to assess the plasma efficacy used for the treatment of wastewater contaminants.Finding the effect of different gases under the same experimental configuration is undoubtedly interesting,considering the future commercialization of the method.

    In this work,cold atmospheric pressure plasma jet system is developed,and its feasibility study is done using the Rhodamine B (Rd.B) and Methylene Blue(Met.Blue) dyes.The major emphasis of the present study is to understand the efficacy of developed plasma jet in the treatment of ECs.For this,a systematic study is performed to obtain optimized parameters using dye as a model contaminant.The selected dyes are common textile water pollutants with reported carcinogenic and mutagenic effects and have a strong chemical structure matching the ECs [21,22].The concentration of dyes was selected according to the ECs concentration range present in sewage water.Another reason for choosing the dyes is due to their simple spectrophotometric analysis.All the experiments were performed using two different plasma gases (argon and helium) under the same reactor configuration to evaluate the effect of plasma discharge gas on treatment efficiency.The effects of other parameters such as sample pH,conductivity,applied voltage,plasma-sample distance were also verified and optimized.Physicochemical characterization of the plasma was performed,and the main plasma reactive species were quantified in the liquid phase.The reactive species analysis helped in understanding the plasma action mechanism.

    After understanding the reaction mechanism and attaining the system’s optimized conditions,treatment of actual ECs was performed.For this,the removal of four pharmaceuticals viz,metformin (C4H11N5),atenolol (C14H22N2O3),acetaminophen(C8H9NO2),and ranitidine(C19H30BiN4O10S)from secondary treated sewage water was performed.The removal of targeted pharmaceuticals in sequencing batch reactor(SBR)and SBR plasma-based hybrid reactor was also evaluated for comparative analysis.On the other hand,the antibiotic resistance bacteria (ARB) are an important class of microbial ECs whose degradation is also very difficult.ARB has a vital role in the development and proliferation of antibiotic resistance,which has been declared the biggest problem of the 21st century by WHO [23].Therefore,the system’s efficiency is also verified against the microbial contaminant removal,and for this,Escherichia coli (E.coli) was selected as model ARB,and its treatment was also done in the same optimized condition used for the pharmaceuticals.

    2.Experimental section

    Figure 1.Schematic of the atmospheric pressure cold plasma jet treatment system.

    2.1.Experimental setup

    The schematic of the experimental setup is shown in figure 1.The plasma jet system is a combination of a high voltage sinusoidal power supply and a plasma reactor.The indigenously developed high voltage sinusoidal AC power supply of 0-6.5 kVp?p,25 kHz was used to generate the plasma.The plasma reactor was made up of a glass tube with an inner and outer diameter of 4 mm and 6 mm,respectively.A 1.5 mm diameter copper rod was used as a high voltage electrode,and a copper wire of 1 mm diameter was wrapped around the glass tube and used as a ground electrode.Electrodes are arranged concentrically,which make a cross-field configuration of the plasma reactor.Argon and helium gases were used independently for the generation of plasma.The position of the electrodes and other plasma reactor configurations was fixed for the entire experiment.Under fixed geometrical conditions,with the applied voltage of 6.5 kVp?p,25 kHz frequency,and gas flow of 3 l min?1for argon and 4 l min?1for helium,a plasma jet length of~3 cm was obtained in the ambient air with both the gases.Typical voltage and current profiles for argon plasma jet are shown in figure S1(available online at stacks.iop.org/PST/24/105501/mmedia).The single jet DBD plasma system used in this study can be scaled up following linear scale-up criteria.Multiple arrays of this plasma jet and appropriate power supply can be used to develop a larger system with increased capacity.

    2.2.Materials and instruments

    The dye and all other chemicals used in this study are of analytical grade and were procured from Sigma-Aldrich,India.The media for microbial analysis was from HiMedia.The colored dye solutions were prepared in RO water(conductivity~0.1 mS cm?1),and for fluorescent dye,Milli Q(conductivity~0.05 μS cm?1),water was used.The commercially available compressed argon and helium gas (99.99%) was used as the plasma discharge gas.The spectrophotometer(JENWAY 7305),pH meter(EUTECH instruments),conductivity meter (Hanna Instruments),Fluorescence meter (Spectra plus 384),Ion chromatography(Metrohm),and TOC analyzer(SHIMADZU)were used for performing the appropriate analysis.The applied voltage was measured using a Tektronix high voltage probe P6015A.The emission profile from the plasma discharge was acquired by the Ocean Optics spectrometer HR-4000 model connected with the optical fiber cable(200 μm)at a 5 mm distance radially near the glass nozzle.

    2.3.Experimental procedure

    The colored dye samples were prepared in RO water,and initially,a stock concentration of 100 ppm was developed.A 1 ppm sample was used for plasma treatment from the stock solution.Experiments were performed with an operating voltage of 6.5 kVp?p(~5 W average power) and 25 kHz sinusoidal voltage.The dye feasibility study was performed with the DBD plasma jet system,and then the other parameters such as pH,conductivity,input voltage,and plasma sample distance were optimized to obtain a complete dye decoloration.All the experiments were performed with fixed applied voltage until specifically mentioned.For the pH experiments,the sample’s initial pH was adjusted using the HCl and NaOH,and for the conductivity experiment,NaCl was used to adjust the sample’s initial conductivity.All the experiments were done in triplicate and then averaged.The decoloration efficiency was calculated using equation (1) after measuring the absorbance of pre-and post-treated samples.The Rd.B and Met.Blue shows maximum UV-visible absorption at a wavelength of 554 and 663 nm,respectively

    For the kinetic study,the results obtained were fitted in the first-order rate equation(equation(2)),from which the apparent rate constant was also calculated

    To determine the extent of dye mineralization under the specific treatment conditions,total organic carbon (TOC) analysis was performed.Finally,to verify the treatment efficiency of the plasma system against the actual ECs,4 pharmaceuticals(metformin,atenolol,acetaminophen,and ranitidine) and a microbial contaminant (E.coli) were treated in the optimized condition obtained from the dye study.

    2.4.Reactive species determination

    A variety of short-lived (few nanoseconds) and long-lived(more than 1 s)reactive oxygen and nitrogen species(RONS)get produced during plasma action at the gas-liquid interface,and the chemical reactions involved in their formation are shown in equations (5)-(17) [12,21,24].The high voltage plasma generates the free electrons (e?) and other excited species,which then ionize the interacting water molecule.This starts a chain reaction,in which the ionized water molecule may break down to form the hydronium ion,free H atom,and ·OH radical.The preference of these interfacial species in air or liquid can be known from their dimensionless Henry constant (KH) value.The molecule with a KHvalue greater than one has more preference in the air compared to liquid.The KHvalue of important reactive species formed in plasma are ·OH 6.92 × 102,H2O21.92 × 106and O31 × 10?1[25].The short-lived·OH initiates the formation of long-lived reactive species.Therefore,identifying the longlived species further affirms the existence of the short-lived species.The bulk diffusion of these reactive species is also crucial as it acts as a rate-limiting step.Those plasma species that interact with the liquid compounds have a faster diffusion than non-reactive species.While,the non-reactive species diffusion occurs simply by the absorption process [26].To understand the mechanism of plasma interaction with contaminants,the analysis of these reactive species is essential.In plasma,·OH,is one of the main reactive species,and its quantification was done using the coumarin-based chemical probe method.It is an indirect method in which the coumarin,a non-fluorescent dye,reacts with ·OH,to form a fluorescent compound 7-hydroxy coumarin as per the mechanism shown in (supplementary figure S2) [27].The fluorescence intensity of 7 hydroxy coumarins is directly proportional to the ·OH concentration [28].Coumarin gives a maximum UV-visible absorption at 277 nm,and 7 hydroxy coumarin,gives the maximum fluorescence at 455 nm after initially exciting at 332 nm [28].1 mM aqueous coumarin sample was used to determine the aqueous ·OH concentration.In the plasma alone,·OH presence was analyzed using optical emission spectroscopy(OES).The H2O2is another essential and longlived reactive species,and its quantification was done using the potassium iodide(KI)based colorimetric method given by Ovenston[29].In this method,the H2O2reacts with KI in the presence of ammonium molybdate ((NH4)6Mo7O24),and iodine (I2) is liberated,which makes the solution yellowcolored (equation (3))

    The intensity of the developed color is directly proportional to the H2O2concentration,and by measuring the absorbance of color at 353 nm and comparing it with a standard H2O2curve,the H2O2concentration in unknown samples can be determined.Ozone(O3)analysis was done using the commercially available O3test kit procured from the Hanna instrument,which is based on the Diethyl-p-phenylene diamine (DPD) method.The nitrogen species(NO3?and NO2?)and other ions/nutrient formations were quantified using ion chromatography.

    2.5.ECs analysis

    Liquid chromatography-mass spectroscopy(LCMS)was used for the pharmaceutical analysis.The difference in pre-and post-treated sample concentration on a percentage basis was used to evaluate the removal efficiency.

    2.5.1.Microbial sample.A overnight grown E.coli cell culture in nutrient broth was centrifuged at 10 000 r min?1for 10 min.The obtained cell pellet was then resuspended in 10 ml PBS buffer.The resuspended solution in triplicate was used as the samples for plasma treatment.The treated sample was then serially diluted using PBS buffer,and suitable dilution was plated on the nutrient agar plate to observe the effect of plasma action on E.coli deactivation.For comparison,the non-treated and control samples were also plated similarly.The E.coli inoculated in PBS buffer acted as the positive control,and PBS buffer without inoculation was used as a negative control.The plates were incubated overnight at 35 ±2 °C.The difference in log CFU ml?1calculated using equation(4)provided the system’s microbial inactivation efficiency.

    where Ci= initial microbial colony count in log CFU ml?1and Cf= microbial colony count in log CFU ml?1after plasma treatment.

    3.Results and discussion

    3.1.RONS quantification

    The ·OH formation during plasma action has been shown in figures 2(a) and (b).The figure shows an increase in the sample’s fluorescence with increased plasma treatment time due to the conversion of coumarin into 7 hydroxy coumarin.This conversion decreases the coumarin concentration in the sample,which can be seen from the decrease in the sample’s absorption spectrum shown in supplementary figures S3(a)and (b).In figures 2(a) and (b),an unusual behavior (instead of a continuous increase,the fluorescence intensity starts decreasing after a particular time) can be seen after 5 min(with argon plasma) and 10 min (with helium plasma).This behavior was observed because the plasma action becomes so prominent that it also starts degrading the 7 hydroxy coumarin molecules after this specific time.Therefore,the sample’s fluorescence intensity decreases.Using a higher concentration dye sample can be a solution,but with coumarin,it is not possible due to its solubility limitation of around 1 mM in the aqueous solution.The conversion of ·OH into other reactive species also occurs,and this also limits the available number of ·OH molecules getting trapped by coumarin to form 7 hydroxy coumarin.This limits the detection of an actual number of ·OH molecules,resulting in a difference in the observed and actual ·OH concentration after some time.This problem is also present with other chemical probes,and previously,Tampieri et al used Terephthalate (TPA) as a probe for ·OH analysis and reported that the TPA formation also deviated from linearity after a certain time [30].Therefore,for maintaining the uniformity in data,the·OH and other plasma species analysis was performed for up to 20 min only.The·OH formation in the plasma jet was detected using OES,and the obtained spectra have been shown in supplementary figure S4 for argon plasma and figure S5 for helium plasma.The characteristic peak at 309 nm(prominent in the·OH band of 306-312 nm) in the OES spectrum is of the ·OH,and the nearby other peaks are of other RONS.The ·OH concentration in the liquid sample (RO water and neutral pH) in a treatment time of 1-20 min is shown in table 1.The maximum estimated ·OH concentration with argon and helium plasma is 14.04 μM and 15.04 μM,respectively.The concentration obtained is in correlation with other similar studies[13,31].The study performed in the reference[32],has used 20 W input power and in 10 min the trend of ·OH concentration obtained with different plasma gases was O2(1.19 mM) > air (1.03 mM) > N2(0.64 mM) > Ar (0.59 mM).Further,Kova?evi? et al,compared five different gases (Ar,He,air,O2,and N2) for ·OH formation in a DBD reactor in direct contact with water.The highest·OH concentration was obtained with oxygen as the plasma feed gas,while a comparable concentration gets produced with the Ar and He gases[33].Due to difficulty in ·OH detection,many studies have also used the fluorescence value or the extent of contaminant degradation to express the ·OH concentration [34,35].

    Figure 2.The fluorescence spectrum of 1 mM coumarin sample treated with (a) argon and (b) helium DBD plasma.(Voltage 6.5 kVp?p).

    In the present study,both the argon and helium gases provided almost a similar concentration of ·OH.In such a case,argon gas is more advantageous as it is a heavier gas,which helps it penetrate more into the liquid and promote the reactive species diffusion in the bulk liquid.Still,we used both gases in all our dye studies to evaluate the difference in treatment efficiency.The ·OH with a half-life of 1 μs and diffusion coefficient,Dgof 2.0 × 10?9m s?2can only travel a maximum distance of 10 μm[32].The diffusivity of·OH in liquid gets increased by its conversion into H2O2(equation (8)),which has more stability and can diffuse deeply inside the liquid [36].Then as per equations (12) and(13),again ·OH gets formed by the breakdown of H2O2.It is reported that 99.6% of total H2O2formed in the plasma system is by the recombination of·OH[25].Therefore,H2O2concentration can also be used as a measure to find the plasma system’s efficiency in producing·OH[37,38].Table 1 shows the H2O2concentration obtained in this study,and it shows an increase in H2O2concentration with increasing treatment time.In 20 min treatment time,we have obtained 562 μM and 448 μM of H2O2with the argon and helium plasma,respectively.With 5 W of average applied power,the H2O2concentration obtained in the present experiments is in a similar range to that obtained by other researchers,such as 680 μM in 10 min with an input power of 20 W [32].Whereas higher concentrations of H2O2i.e.1765 μM and 1176 μM with 20 min treatment time for 150 W,and 100 W applied power,respectively were obtained in [38].In the present work,the H2O2concentration obtained was very high compared to ·OH because before trapping the ·OH molecule by coumarin,it gets converted into H2O2or any other reactive species.Due to the interconversion of reactive species,a major portion of the generated H2O2also gets lost by decomposition to form the H2O and O2(equation (11)) [32].

    Table 1.Formation of RONS in DBD plasma.(Sample DI water,voltage 6.5 kVp?p.)

    Ozone (O3) is another important stable reactive species formed during the plasma action (equation (15)).The ozone formation was also evaluated in the RO water,and similar to H2O2,a continuous increase in its concentration was obtained.A maximum O3concentration of 320 μM and 111 μM in 20 min of plasma treatment was obtained with the argon and helium plasma,respectively (table 1).Multiple studies have reported that the formation of O3in plasma reactors takes place mainly when O2or air has been used as the plasma feed gas[33,39].In a study by Kinandana et al,it was also found that O3formation has also been getting significantly affected by the system geometry.In their studies,the researchers have compared the series and parallel DBD configurations and reported that the series configuration was better as it provided 2000 μM more O3concentration than the parallel configuration [40].Due to the operation of the reactor in an open atmosphere,a few O2molecules may entrain into the ionization region and convert into O3in the gaseous phase.Due to the limited solubility of O3in neutral pH,a smaller number of O3molecules get dissolved into the liquid sample.Whereas,higher dissolved O3concentration has been reported in those plasma systems which changes the sample pH to acidic [41].On the other hand,during plasma action,O3is also produced from the conversion of ·OH.The series of reactions involving this conversion process has been represented in the equations(5)-(17).In a study by Teng et al,the influence of O3on the formation of NO2?and NO3?was studied.In this study,the researchers have used different ratios of O2:N2in plasma feed gas [42]and found a positive effect of O3on nitrogen species,especially NO3?.Whereas in the present experiments,we have obtained a low O3concentration,and also,there was no direct source of nitrogen.Therefore,very less NO3?and NO2?formation occurred,and it was beyond the detectable range of ion chromatography.In literature also,NO3?and NO2?formation was mainly reported when nitrogen gas or air was used as the plasma feed gas[33,43,44].The low nitrogen formation by our system provides an additional advantage,as these nitrogen species themselves are an important water contaminants [45].In the study by Meropoulis et al,it has been concluded that these nitrogen species do not have much role in pollutant degradation and cause a decrease in sample pH during plasma action [46].Similar findings have also been reported in another study [33].Due to the formation of these nitrogen species in undetectable concentrations,no significant pH change in a post-treated sample was observed in our case.After 30 min of plasma treatment time,the sample pH changed from neutral to slightly acidic (pH 6) which was due to the formation of organic acids formed from the dye degradation by-products [47].The concentration of these organic acids is also very less,due to the system’s high mineralization efficiency,about which a detailed description has been provided in section 3.9.Finally,after analyzing the results of reactive species formation and understating their chemistry with the contaminants[9,34,44],it is apparent that the ·OH and H2O2are the main reactive species responsible for the dye degradation during plasma action.

    3.2.Effect of pH

    In any of the AOPs,the sample pH is considered to influence the treatment efficiency.To evaluate the effect of pH on dye decoloration,the samples with acidic (pH 4),neutral (pH 7),and basic (pH 10) pH were treated.The results have been shown in figures 3(a)-(d).The data represents a minor or almost negligible effect of pH on the degradation of the model dyes.Specifically,samples treated with helium have shown marginally higher decolorization in a shorter treatment time(in 5 min > 85%),whereas it took a little longer for argon.

    Overall,the acidic pH was found to slightly favor the dye decoloration process.The acidic pH samples of Rd.B and Met.Blue showed 100% removal,while complete color removal was not obtained in the basic and neutral pH samples.A higher contaminant degradation in an acidic pH was also reported by other researchers[48,49].While the basic pH favors those reactions in which the important RONS gets decomposed.Such as the reaction between the ·OH and H2O2in which hydroperoxyl (HO2·) radicals get formed(equation (9)),which has less oxidation potential than ·OH[50].The ·OH radical has 100 times faster reactivity(k = 7.1 × 109M?1s?1) with HO·2than H2O,and in this reaction,few ·OH molecules get consumed (equation (12)).The self-decomposition of H2O2is also favored in basic pH according to equations(9)and(11).In basic pH,the OH?,CO3?,HCO3?,etc,ions are present,and they also act as the·OH scavenger.With this in basic pH,there is an easy conversion of ·OH into its conjugate base O·?.The ·OH acts as an electrophile while reacting with organic molecules,and the O·?is a nucleophile.Therefore,these two radicals can form different intermediate products,leading to different degradation pathways [51].Further,it has been reported that plasma treatment shifts the treated sample pH towards the acidic side due to nitric/nitrous acid formation.These acids are formed from the nitrogen species(NO2?and NO3?),mainly generated when air or nitrogen is used as the plasma feed gas [52,53].A very less pH variation has been observed when nitrogen was not present in the plasma feed gas [43,54].Bolouki et al compared the pH variation in the plasma-treated sample when argon,oxygen,nitrogen,and air were used as the plasma feed gas [54].They reported that no significant pH changed with argon and oxygen,while pH decreased sharply with air and nitrogen plasma.In the same study,it was also found that the nitrogen species also acts as the ·OH quencher,which was evidenced by the absence of·OH peak in air and nitrogen plasma OES spectra.As we have also used argon and helium as the plasma feed gas,not much pH drop in the treated sample was observed.

    3.3.Effect of conductivity

    The sample conductivity strongly affects the plasma discharge characteristics,which further influences the reactive species generation.The conductivity becomes more critical when the discharge is applied directly over the liquid solution.In this case,the plasma behavior such as discharge channel,current,length,and propagation depends on the sample conductivity [55].

    Figure 3.Dye decoloration at different pH values with increasing treatment time.(a) Rd.B-argon,(b) Rd.B-helium,(c) Met.Blue-argon,and (d) Met.Blue-helium DBD plasma.

    Our study target is to check the DBD cold plasma feasibility for wastewater treatment.The wastewater generally has concentrated salts and inorganics,increasing its conductivity to the range of tens of mS cm?1.Therefore,the conductivity experiment was done at a neutral pH and in the conductivity range of non-altered (n.a),5,10,and 20 mS cm?1.The n.a sample has a conductivity value of~0.16 mS cm?1for Rd.B and~0.28 mS cm?1for Met.Blue.The results have been shown in figures 4(a)-(d),which depict that for the first 5 min,a high dye removal was observed with both the dyes and the gases,and from 5 to 15 min the dye decoloration increased gradually.After this,almost a saturation is obtained,and for non-altered (n.a) and 5 mS cm?1conductivity samples,complete dye decoloration was obtained in 30 min.In the same duration with 10 and 20 mS cm?1conductivity samples,dye decoloration was obtained only up to 90%.This shows a decrease in treatment efficiency with increasing conductivity.Such behavior was observed because the plasma discharge becomes more adverse at higher conductivity,resulting in plasma conversion to spark discharge with an increment in plasma intensity and discharge current for both argon and helium gases.Similar observations were reported in other studies where the decrease in reactive species formation was noticed with an increase in intensity and discharge current as the sample conductivity increased[44,56].This reduction in reactive species could be the possible reason for reduced treatment efficiency in the present case.An opposite behavior is also reported in a study where a much higher conductivity sample (up to 100 mS cm?1) was used [49].The high conductivity was favorable in their case because the other plasma effect(ozone and UV)became more prominent,which ultimately improved the plasma efficiency,but with a high energy requirement.

    3.4.Effect of applied voltage

    The applied voltage directly influences the formation of highenergy electrons,which results in more reactive species generation.Therefore,the decoloration efficiency should increase with an increase in applied voltage.Therefore,the present system’s potential in different applied voltage values was evaluated.The chosen voltage values were 4.5,5.5,and 6.5 kVp?p,and their corresponding average discharge power values were 3.1,3.9,and 4.9 W respectively.The results have been shown in figures 5(a)-(d),which shows the complete removal of Rd.B and Met.Blue in 30 min with 6.5 kVp?papplied voltage,while it was only 91%for Rd.B with argon,and 93% for Met.Blue with argon at 4.5 kVp?papplied voltage.With lower power to achieve the same level of removal,a longer time is needed.This shows that the plasma efficiency depends on the applied voltage,which is an indirect measure of applied power.In higher power,more dense charged electron species are generated with increased physical effects resulting in increasing the efficiency of the treatment process [55].

    Figure 4.Effect of sample’s initial conductivity on dye decoloration with increasing treatment time.(a)Rd.B-argon,(b)Rd.B-helium.(c)Met.Blue-argon and (d) Met.Blue-helium DBD plasma.

    Figure 5.Applied voltage effect on dye decoloration with increasing treatment time.(a) Rd.B-argon,(b) Rd.B-helium,(c) Met.Blue-argon and (d) Met.Blue-helium DBD plasma.

    3.5.Effect of sample distance from the plasma

    With the increase in plasma sample distance,the gaseous phase short-lived RONS gets decayed before they get diffused inside the liquid[30].Therefore,finding an optimum distance helps to obtain a point at which the maximum plasma penetration in the liquid could be obtained.The highest removal was obtained when the sample was placed at 2 cm,i.e.at the middle of the plasma plume.At 1 cm distance(top of plasma plume),the dye decoloration gets decreases,while at 3 cm distance (tip of plasma),the least removal was observed.Few reactive species are present at the plasma tip,and plasma penetration efficiency is also not strong enough if the sample plasma distance increases.While in a much closer distance,the coming discharge gas gets bounced back after hitting the water layer,as of which suitable plasma plume formation did not take place.In addition to this,narrower plasma-sample distance fluctuation in the liquid surface causes improper Tylor cone formation,resulting in instability in the top layers of liquid,which also affects the formation of plasma discharge [52,57].Therefore,keeping the sample at an optimum distance,which in our case was around the middle of the plasma plume,would provide maximum efficiency.

    3.6.Kinetic study

    The dye samples were treated for a time interval of 0-30 min to find the effect of increasing treatment time on the dye decoloration process.For both the dyes,during the initial 5 min,the dye decoloration was exponential,which becomes gradually up to 15 min.After this,a saturating behavior was observed until complete removal,as shown in figures 6(a) and (b).A similar decoloration pattern of both the dye was obtained with both the argon and helium gas.The 90% removal of both the dyes was obtained in 20 min of plasma treatment,but it took an extra 10 min for complete dye decoloration.The extra time was needed to degrade chromophores and resistive intermediates released from the breakdown of dye.In an extended treatment time,the direct point treatment of liquid results in mild heating of the samples.As a result,a slight decrease in liquid volume due to evaporation was observed.In this study,during 30 min of treatment,the maximum temperature of the sample reached was around 50 °C.

    The ln(C0/Ct)versus time graph was plotted,which has been shown in figure 6(c),and the apparent or observed rate constant k was calculated,using which the half-life of the dye molecule during plasma treatment was calculated.The dye decoloration results were consistent with the first-order rate equation,and therefore the rate constant depends on the dye concentration at any given time.Table 2 lists the rate constant for the initial 5 min and the overall dye decoloration reaction.In all cases,the rate constant for the initial 5 min is more than the overall reaction,which is also reflected in the t1/2value of the dye molecule.Initially,breakdown of the complete dye molecule occurs,but as the process continues,resistive intermediates get formed,due to which the dye decoloration rate decreases.In the extended time,the dye decoloration also loses its linearity,and the reaction gets shifted from first order to pseudo-first-order rate equation.

    Table 2.Kinetic parameters of the dye degradation process for the initial 5 min and complete dye decoloration process.

    Table 3.The removal of the pharmaceutical compounds in a conventional method and DBD plasma treatment system.(SBR treatment time-24 h,and volume-1 l).

    Figure 6.Kinetic study of dye decoloration (a) for Rd.B,(b) for Met.Blue and (c) logarithmic representation of (C0/Ct) versus time graph for determining kinetic constants of argon plasma-treated Rd.B sample.

    3.7.Effect of dye concentration

    To evaluate the efficiency of plasma on a different level of contamination,samples with an increasing dye concentration were treated.The result of increasing dye concentration has been depicted in figure S6,from which it can be concluded that increasing dye concentration limits the plasma treatment efficiency.As in all cases,the maximum dye decoloration occurred during the initial 5 min of treatment.Therefore,for this experiment,sample treatment was done only for 5 min.Around 80% removal was observed during the selected treatment time in the 1,5,and 10 ppm samples with both argon and helium gas.In the 50 ppm sample,the removal gets lowered by 50%,and an average removal of 40% was observed.The majority of the ECs in wastewater were found to be in the range of a few ppm;therefore,the plasma treatment will be effective against them.If,in some specific type of wastewater,a higher concentration of ECs will be present,then by increasing the treatment time or input voltage of the DBD plasma system,a high treatment efficiency can be achieved.

    3.8.Individual effects of RONS

    It is reported that certain molecules with higher reactivity with RONS act as the scavenger molecule and consume the RONS before they could oxidize the contaminants [58,59].Methanol is such a scavenger molecule for ·OH,and its presence would decrease the plasma effectivity.Here we used this property of methanol for evaluating the individual effect of·OH on the reactor’s treatment efficiency [60].For this,10%methanol acting as the·OH radical quencher was added to the dye sample,and dye decoloration was compared with the sample without any added methanol.Like this,the individual effect of other main reactive species was also determined.To find out the individual effect of H2O2,the samples were treated with externally added H2O2.A major limitation in determining the individual effect of H2O2during plasma treatment is maintaining the H2O2concentration at a constant level.Due to the continuous formation of H2O2,its concentration gets changed.Still,to get a rough idea of the effect of H2O2,the highest reported concentration of H2O2obtained during plasma treatment in the current study was used as the sole treating reagent in a batch reactor.In the batch reactor,only H2O2in a concentration of 600 μM was added to the dye samples,and treatment was performed for 5 min.The results obtained from this experiment were compared with the experiment conducted using plasma for the same treatment time.The results obtained from these two studies have been compared in figure S7.In the sample where the ·OH was quenched,the least dye decoloration was observed,42.45%(Rd.B-Ar) and 44.21% (Met.Blue-He).While the H2O2treated samples have shown a decoloration of 68.33%(Rd.B)and 66.95% (Met.Blue).In figure S8 the results of dye decoloration with increasing H2O2concentration (50-600 μM) could also be found.As,it can be seen that in ·OH quenched samples least dye decoloration was obtained,which affirms that ·OH has the highest contribution in plasma treatment.The·OH with a direct attack on dye molecule also has a role in the formation of other reactive species.

    3.9.Dye mineralization and presumption of degradation mechanism

    Mineralization is the efficiency of a treatment process to convert the pollutant into mineral acids,H2O,and CO2.In the absence of proper mineralization,the degradation by-products having similar or higher toxicity can get formed.Therefore,the decoloration experiment was followed by the mineralization experiment to assess the system’s mineralization efficiency,and it was expressed in terms of total organic carbon (TOC).The decrease in the TOC value shows a decrease in the effluent toxicity level and the formation of mineralized products.The results of the mineralization experiment have been shown in table 2,which shows that argon as plasma discharge gas provided more mineralization 73% (Rd.B-Ar),60% (Met.Blue-Ar),compared to helium with 45% (Rd.B-He),and 35% (Met.Blue-He).Compared to the argon plasma,low mineralization was achieved with the helium plasma.This shows that although the helium has high dye decoloration efficiency,its mineralization efficiency is low compared to argon.The mineralization efficiency obtained in this study was higher than other similar studies performed on dye treatment [61].The results obtained show that the present treatment system with high decoloration can also provide good dye mineralization.The degradation pathway by which the mineralization of Rd.B [62]and Met.Blue [4,32]occurs has been well explained in the respective literature.

    3.10.Pharmaceuticals and microbial contamination removal

    The detailed dye decoloration study showed that both the plasma discharge gases provided similar efficiency under the same reactor configuration.Compared to helium,the argon plasma provided a higher dye mineralization efficacy,and it is also a relatively cheaper gas.Therefore,argon gas was selected to perform further studies.The system’s efficiency against the actual ECs present in wastewater was determined by treating the sequencing batch reactor (SBR) treated secondary wastewater and targeting the removal of 4 pharmaceuticals.Efficient removal of all the selected pharmaceuticals was observed,and the results are shown in table 3.The table also includes the comparative results of the SBR,plasma,and SBR-hybrid plasma system.The SBR has a working volume of 1 l,a cycle time of 24 h,and 10% sludge,and it shows a high removal because of the contaminants bioadsorption by the sludge.The 5 min argon plasma treatment also showed a high removal of pharmaceuticals except for the metformin.Almost a complete removal is observed when the plasmahybrid treatment system was used.This shows that a combination of plasma and conventional wastewater treatment systems is an effective ECs treatment method.

    Figure 7.Microbial treatment with plasma.(a) The change in microbial concentration with time,(b) media plates showing the removal of microbial colonies with increasing treatment time,and pre and post-treated (c) E.coli and,(d) B.subtilis.(C0-log8 CFU ml?1).

    Finally,the plasma system efficiency was also evaluated with the microbial contaminants considering E.coli as model antibiotic resistance bacteria (ARB).The E.coli was treated using both the argon and helium plasma,and argon plasma was found to be more effective with log2.5 E.coli deactivation in 6 min of plasma treatment (figure 7(a)).The typical image of argon plasma-treated sample with increasing treatment time has also been shown in figure 7(b),and the absence of E.coli colony can be seen in 6 min treated sample plate.Cell wall degradation is an important mechanism of microbial deactivation by plasma action [25,63].The E.coli (Gramnegative)has a thin cell wall(<10 nm),and to check plasma action on other cell wall types Bacillus subtilis,a Grampositive bacteria with a thick cell wall (20-80 nm),was also treated under the same treatment condition.The plasmatreated samples show a log removal of 2.45 for E.coli and 2.42 for B.subtilis(figures 7(c)and(d)).Therefore,a similar removal was observed with both the bacterial types,from which it can be concluded that plasma acts non-specifically and is equally effective against different bacterial types.

    4.Conclusion

    In this study,the DBD plasma jet system assessment against the ECs treatment was performed using Rd.B and Met.Blue as model contaminants.The system’s operational parameters were optimized using the dye as the model contaminant,the fastest dye decoloration up to 80% in all cases occurred during the first 5 min,and a complete dye decoloration was achieved within 30 min of treatment.By the analysis of RONS species,the formation of ·OH,H2O2,and O3was confirmed.With further detailed analysis,it was also found that the ·OH and H2O2act as the main action mechanism during plasma action.The system was also efficient in providing high dye mineralization which is extremely useful in ensuring the absence of any toxicity in the effluent.A comparative study of argon and helium as plasma discharge gas was performed,and it was observed that both gases provide a similar treatment efficiency.Considering the treatment efficiency and cost-effectivity,the utilization of argon gas is suggested.Finally,under the optimized conditions,the plasma system was highly effective against the plasma treatment of sewage water (consisting of four pharmaceutical chemicals and microbial) contaminants.The pharmaceutical removal was also checked in SBR (conventional treatment system) and SBR plasma hybrid reactor for the comparative study.The complete removal of the pharmaceuticals was observed in a plasma SBR hybrid reactor.These results show the efficacy of a developed non-thermal plasma system for the treatment of ECs.Moreover,the present systematic study of plasma system operational parameters and optimization for dye decoloration will aid in establishing the technology for large scale ECs treatment.

    Acknowledgments

    This work was supported by grants from the IIT Delhi FIRP program grant(No.MI02081).Authors thank the referees for their useful suggestions.DJ and GVP thank Mahreen for technical help on spectroscopy analysis.

    ORCID iDs

    舔av片在线| 亚洲国产欧美在线一区| 22中文网久久字幕| 欧美xxⅹ黑人| 国产高清不卡午夜福利| 在线看a的网站| 91午夜精品亚洲一区二区三区| 2022亚洲国产成人精品| 街头女战士在线观看网站| av免费观看日本| 一级毛片 在线播放| xxx大片免费视频| 久久久精品欧美日韩精品| 少妇猛男粗大的猛烈进出视频 | 欧美日本视频| 亚洲av国产av综合av卡| 舔av片在线| 国产美女午夜福利| 亚洲精品一二三| 91午夜精品亚洲一区二区三区| 日韩av免费高清视频| 婷婷色综合www| 内地一区二区视频在线| 欧美日韩综合久久久久久| 少妇丰满av| 观看美女的网站| 亚洲人与动物交配视频| 中文资源天堂在线| 精品人妻视频免费看| 青春草亚洲视频在线观看| 最近2019中文字幕mv第一页| 亚洲欧美成人精品一区二区| 欧美日韩视频精品一区| 超碰av人人做人人爽久久| 下体分泌物呈黄色| 欧美日本视频| 精品国产乱码久久久久久小说| 一级毛片久久久久久久久女| 香蕉精品网在线| 特级一级黄色大片| 视频区图区小说| 亚洲av一区综合| av卡一久久| 国产美女午夜福利| 久久精品久久久久久久性| 天堂网av新在线| 日本免费在线观看一区| 自拍欧美九色日韩亚洲蝌蚪91 | 日本爱情动作片www.在线观看| 麻豆精品久久久久久蜜桃| 精品久久久久久久人妻蜜臀av| av在线蜜桃| 欧美精品国产亚洲| 51国产日韩欧美| 成人午夜精彩视频在线观看| 久久99精品国语久久久| 国产精品女同一区二区软件| 色视频www国产| av在线亚洲专区| 九草在线视频观看| 看免费成人av毛片| 亚洲国产精品999| 搞女人的毛片| 国产毛片a区久久久久| 国产毛片a区久久久久| 黄色一级大片看看| 免费大片黄手机在线观看| 97超碰精品成人国产| av黄色大香蕉| 国产久久久一区二区三区| 少妇被粗大猛烈的视频| 赤兔流量卡办理| 亚洲天堂国产精品一区在线| 高清毛片免费看| 亚洲图色成人| 丝袜喷水一区| 80岁老熟妇乱子伦牲交| 久久久久国产网址| 大话2 男鬼变身卡| av播播在线观看一区| 成人国产麻豆网| 大片电影免费在线观看免费| 午夜亚洲福利在线播放| 好男人视频免费观看在线| 乱系列少妇在线播放| 亚洲美女搞黄在线观看| 久久久久久久亚洲中文字幕| 精品国产一区二区三区久久久樱花 | 亚洲四区av| 99久久人妻综合| 夫妻性生交免费视频一级片| 久久久国产一区二区| 日本欧美国产在线视频| 成人午夜精彩视频在线观看| 少妇被粗大猛烈的视频| 国产探花在线观看一区二区| 国产高清不卡午夜福利| 精品久久国产蜜桃| 中文字幕亚洲精品专区| 91精品一卡2卡3卡4卡| 最近中文字幕2019免费版| 国产乱人视频| 我的老师免费观看完整版| 91久久精品电影网| 亚洲精品,欧美精品| 一级爰片在线观看| 老女人水多毛片| 人妻一区二区av| 国产精品一区二区性色av| 国内揄拍国产精品人妻在线| 男女国产视频网站| 一级毛片我不卡| 国产一级毛片在线| 精品久久久久久久人妻蜜臀av| 99re6热这里在线精品视频| 免费观看av网站的网址| 综合色av麻豆| 欧美高清性xxxxhd video| 最近中文字幕2019免费版| 男女下面进入的视频免费午夜| 五月伊人婷婷丁香| eeuss影院久久| 亚洲综合色惰| 成人特级av手机在线观看| 69人妻影院| 热99国产精品久久久久久7| 久久人人爽人人片av| 秋霞在线观看毛片| 欧美 日韩 精品 国产| 国产v大片淫在线免费观看| 国产亚洲一区二区精品| 色视频www国产| 亚洲不卡免费看| 久热这里只有精品99| 中文资源天堂在线| 久久久久久久久久久免费av| 中文天堂在线官网| 国产伦精品一区二区三区视频9| 色哟哟·www| 一区二区三区乱码不卡18| 男人狂女人下面高潮的视频| 成人午夜精彩视频在线观看| 寂寞人妻少妇视频99o| 最新中文字幕久久久久| 午夜免费鲁丝| 少妇 在线观看| 男人添女人高潮全过程视频| 亚洲av二区三区四区| 久久国产乱子免费精品| 国产亚洲最大av| 五月天丁香电影| 国产色爽女视频免费观看| 久久女婷五月综合色啪小说 | 丝袜脚勾引网站| 日韩一区二区视频免费看| 国产精品福利在线免费观看| 午夜福利视频1000在线观看| 欧美高清成人免费视频www| 欧美日韩视频精品一区| 日本一二三区视频观看| 国产成人一区二区在线| 成人国产麻豆网| 亚洲欧美清纯卡通| 国产成人一区二区在线| 国产成年人精品一区二区| 国产美女午夜福利| 亚洲成色77777| 中文精品一卡2卡3卡4更新| 少妇熟女欧美另类| 欧美区成人在线视频| 大片电影免费在线观看免费| 亚洲精品国产av蜜桃| 人体艺术视频欧美日本| 欧美一区二区亚洲| 最近最新中文字幕大全电影3| 成人毛片60女人毛片免费| 舔av片在线| 亚洲人成网站在线播| 国产av码专区亚洲av| 夫妻午夜视频| 亚洲天堂国产精品一区在线| 国产 一区精品| 能在线免费看毛片的网站| 少妇的逼水好多| 欧美日韩在线观看h| 亚洲精品影视一区二区三区av| 成人综合一区亚洲| 黄片wwwwww| 日本一二三区视频观看| 不卡视频在线观看欧美| 国产男人的电影天堂91| 亚洲人成网站在线播| 亚洲av日韩在线播放| 99久久人妻综合| 免费观看av网站的网址| av在线播放精品| 国产高清三级在线| 97人妻精品一区二区三区麻豆| 亚洲综合精品二区| 国产伦在线观看视频一区| 亚洲av.av天堂| 少妇熟女欧美另类| av在线蜜桃| 69人妻影院| 丝袜喷水一区| 97人妻精品一区二区三区麻豆| av国产久精品久网站免费入址| 久久久精品免费免费高清| 黄片wwwwww| 午夜爱爱视频在线播放| 人人妻人人看人人澡| 日本-黄色视频高清免费观看| 亚洲av欧美aⅴ国产| 久久久久久久久大av| 久久人人爽av亚洲精品天堂 | 亚洲精品国产成人久久av| 欧美xxxx黑人xx丫x性爽| 五月天丁香电影| 久久久久久久国产电影| 真实男女啪啪啪动态图| 尾随美女入室| 欧美日韩一区二区视频在线观看视频在线 | 岛国毛片在线播放| 免费观看在线日韩| 国内揄拍国产精品人妻在线| 久久影院123| 在线免费十八禁| 成人综合一区亚洲| 麻豆久久精品国产亚洲av| 日韩av免费高清视频| 熟女人妻精品中文字幕| 午夜视频国产福利| 亚洲一区二区三区欧美精品 | 美女xxoo啪啪120秒动态图| 午夜免费鲁丝| 搡老乐熟女国产| 国产精品国产三级国产专区5o| 国产成人午夜福利电影在线观看| 亚洲欧美日韩无卡精品| 亚洲精品影视一区二区三区av| 国产一区有黄有色的免费视频| 身体一侧抽搐| 国产高清国产精品国产三级 | 少妇人妻精品综合一区二区| 国产精品人妻久久久久久| 国产精品无大码| 国产免费一区二区三区四区乱码| 少妇的逼水好多| 少妇人妻 视频| 久久久色成人| 免费人成在线观看视频色| 欧美亚洲 丝袜 人妻 在线| 女人十人毛片免费观看3o分钟| 直男gayav资源| 五月开心婷婷网| 成人国产麻豆网| 最后的刺客免费高清国语| 久久精品国产亚洲av天美| av福利片在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精品中文字幕在线视频 | 久久影院123| 国产视频内射| 欧美xxⅹ黑人| 亚洲婷婷狠狠爱综合网| 日韩,欧美,国产一区二区三区| 国产男女内射视频| 国产色婷婷99| 最近最新中文字幕免费大全7| 特大巨黑吊av在线直播| 精品人妻视频免费看| 乱系列少妇在线播放| 美女脱内裤让男人舔精品视频| 少妇裸体淫交视频免费看高清| 亚洲丝袜综合中文字幕| 国产欧美另类精品又又久久亚洲欧美| 中文在线观看免费www的网站| av免费观看日本| 国产中年淑女户外野战色| 日本免费在线观看一区| 一级毛片我不卡| 黄色配什么色好看| 免费大片18禁| 人妻一区二区av| 97在线视频观看| 99热这里只有是精品50| 国内揄拍国产精品人妻在线| 欧美极品一区二区三区四区| 少妇人妻一区二区三区视频| 日本一二三区视频观看| 亚洲无线观看免费| 欧美少妇被猛烈插入视频| 在线亚洲精品国产二区图片欧美 | 丝袜脚勾引网站| 一级毛片黄色毛片免费观看视频| 一级毛片我不卡| 九草在线视频观看| 亚洲,一卡二卡三卡| 国产伦理片在线播放av一区| 男女国产视频网站| 人人妻人人看人人澡| 久久人人爽人人爽人人片va| 国产高清三级在线| 日韩三级伦理在线观看| 99热这里只有是精品50| 国产一区亚洲一区在线观看| 国产乱人偷精品视频| 日韩伦理黄色片| 国产精品久久久久久精品古装| 亚洲激情五月婷婷啪啪| 最近的中文字幕免费完整| 三级国产精品欧美在线观看| 高清欧美精品videossex| 网址你懂的国产日韩在线| 在线观看国产h片| 亚洲丝袜综合中文字幕| 国产一区二区亚洲精品在线观看| a级一级毛片免费在线观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲精品,欧美精品| 国产黄色视频一区二区在线观看| 男女啪啪激烈高潮av片| 狂野欧美白嫩少妇大欣赏| 一级毛片黄色毛片免费观看视频| 国语对白做爰xxxⅹ性视频网站| 国产精品.久久久| 91在线精品国自产拍蜜月| 国产精品不卡视频一区二区| 听说在线观看完整版免费高清| 一级毛片aaaaaa免费看小| 亚洲国产欧美在线一区| 日韩电影二区| 国产成人精品婷婷| 亚洲三级黄色毛片| 日本猛色少妇xxxxx猛交久久| 免费黄频网站在线观看国产| 亚洲怡红院男人天堂| 亚洲人成网站在线播| 日韩中字成人| 久久久久久久国产电影| 草草在线视频免费看| 如何舔出高潮| 精品国产乱码久久久久久小说| 欧美3d第一页| 国产 一区精品| 成人鲁丝片一二三区免费| 大话2 男鬼变身卡| 一区二区三区四区激情视频| 国产色婷婷99| 搡女人真爽免费视频火全软件| 国产精品久久久久久精品电影小说 | 免费黄网站久久成人精品| 国产乱人视频| 午夜老司机福利剧场| 色5月婷婷丁香| 2018国产大陆天天弄谢| videos熟女内射| 欧美bdsm另类| 午夜激情福利司机影院| 在线看a的网站| 好男人视频免费观看在线| 精品视频人人做人人爽| 亚洲熟女精品中文字幕| av国产久精品久网站免费入址| 人人妻人人爽人人添夜夜欢视频 | 涩涩av久久男人的天堂| 美女主播在线视频| 亚洲国产精品成人久久小说| 欧美最新免费一区二区三区| 在线观看美女被高潮喷水网站| 老女人水多毛片| 又大又黄又爽视频免费| 久久精品国产亚洲网站| 两个人的视频大全免费| 成人无遮挡网站| 久久久久国产网址| av在线天堂中文字幕| 国产色爽女视频免费观看| 日韩精品有码人妻一区| 国产v大片淫在线免费观看| freevideosex欧美| 婷婷色综合www| 夫妻性生交免费视频一级片| 黄色怎么调成土黄色| 波多野结衣巨乳人妻| 国产男人的电影天堂91| 最近最新中文字幕大全电影3| 精品国产三级普通话版| 国产精品伦人一区二区| 日日啪夜夜撸| 国产老妇女一区| 搡老乐熟女国产| 久久精品熟女亚洲av麻豆精品| av免费在线看不卡| 女人被狂操c到高潮| 日韩欧美 国产精品| 欧美日韩国产mv在线观看视频 | 国产免费一级a男人的天堂| 少妇猛男粗大的猛烈进出视频 | 99九九线精品视频在线观看视频| 成人综合一区亚洲| 高清日韩中文字幕在线| 久久久久久久国产电影| 亚洲自拍偷在线| 亚洲精品成人av观看孕妇| 午夜福利网站1000一区二区三区| 69人妻影院| 亚洲,欧美,日韩| 嫩草影院精品99| 国产色婷婷99| 哪个播放器可以免费观看大片| 秋霞在线观看毛片| 亚洲av欧美aⅴ国产| 日本一二三区视频观看| 久久久久久久久久久免费av| 欧美性猛交╳xxx乱大交人| 国产色婷婷99| 一级毛片 在线播放| 久久精品久久久久久久性| 亚洲精品aⅴ在线观看| 色播亚洲综合网| 精品久久久久久久人妻蜜臀av| 亚洲国产av新网站| 亚洲成人一二三区av| 永久免费av网站大全| 能在线免费看毛片的网站| 免费看光身美女| 丰满乱子伦码专区| 丝袜喷水一区| 搡女人真爽免费视频火全软件| 免费看av在线观看网站| 中文天堂在线官网| 日韩免费高清中文字幕av| 国产成人精品福利久久| 青春草亚洲视频在线观看| 一级毛片aaaaaa免费看小| 99久久中文字幕三级久久日本| 99热这里只有精品一区| 亚洲精品乱码久久久久久按摩| 在线免费十八禁| 国产亚洲5aaaaa淫片| 香蕉精品网在线| 国产精品久久久久久精品电影| 丝袜喷水一区| 国产亚洲最大av| 不卡视频在线观看欧美| 性色avwww在线观看| 边亲边吃奶的免费视频| 成年免费大片在线观看| 少妇人妻 视频| 亚洲人成网站高清观看| 激情 狠狠 欧美| 亚洲国产日韩一区二区| 高清欧美精品videossex| 国产永久视频网站| 极品教师在线视频| 爱豆传媒免费全集在线观看| 欧美高清成人免费视频www| 啦啦啦中文免费视频观看日本| av天堂中文字幕网| 亚洲欧美中文字幕日韩二区| 只有这里有精品99| 国产在线一区二区三区精| 国产亚洲91精品色在线| 国产成人一区二区在线| 亚洲国产最新在线播放| av网站免费在线观看视频| 亚洲精品第二区| 黄片wwwwww| av在线老鸭窝| 精品一区在线观看国产| 女人久久www免费人成看片| 2021少妇久久久久久久久久久| 国内少妇人妻偷人精品xxx网站| 国产精品一二三区在线看| 国内揄拍国产精品人妻在线| 男人添女人高潮全过程视频| 午夜免费鲁丝| 亚洲国产高清在线一区二区三| 久久久久国产精品人妻一区二区| 男女下面进入的视频免费午夜| 久久久久久久亚洲中文字幕| av在线天堂中文字幕| 人人妻人人澡人人爽人人夜夜| 国产精品嫩草影院av在线观看| 亚洲av一区综合| 国产白丝娇喘喷水9色精品| 在现免费观看毛片| 自拍偷自拍亚洲精品老妇| 99久久人妻综合| 亚洲内射少妇av| 一个人看视频在线观看www免费| 欧美成人一区二区免费高清观看| 中文字幕av成人在线电影| 日韩一区二区三区影片| 亚洲伊人久久精品综合| 日韩一区二区视频免费看| 69人妻影院| 免费观看在线日韩| 热re99久久精品国产66热6| 精品久久久久久电影网| 少妇的逼水好多| 久久韩国三级中文字幕| 国内少妇人妻偷人精品xxx网站| 久热久热在线精品观看| 在线观看三级黄色| 伊人久久精品亚洲午夜| 午夜福利视频1000在线观看| 国产精品成人在线| 真实男女啪啪啪动态图| 99九九线精品视频在线观看视频| 最近中文字幕2019免费版| 在线免费观看不下载黄p国产| 黄片wwwwww| 亚洲三级黄色毛片| 亚洲欧美日韩卡通动漫| 亚洲av欧美aⅴ国产| 中文字幕制服av| 午夜激情久久久久久久| 九草在线视频观看| 国产人妻一区二区三区在| 亚洲精品日韩av片在线观看| 国产精品麻豆人妻色哟哟久久| 欧美一区二区亚洲| 天堂中文最新版在线下载 | 人妻夜夜爽99麻豆av| 又爽又黄无遮挡网站| 亚洲av.av天堂| 在线播放无遮挡| 一区二区三区乱码不卡18| 免费高清在线观看视频在线观看| 中文在线观看免费www的网站| 久久影院123| 精品国产露脸久久av麻豆| 国产中年淑女户外野战色| www.色视频.com| 亚洲人成网站在线观看播放| 国模一区二区三区四区视频| 亚洲精品久久久久久婷婷小说| 亚洲久久久久久中文字幕| 久久精品综合一区二区三区| 国产精品.久久久| 国产成人aa在线观看| 国产精品秋霞免费鲁丝片| 黄色欧美视频在线观看| 激情五月婷婷亚洲| 色5月婷婷丁香| 日本猛色少妇xxxxx猛交久久| 18禁裸乳无遮挡免费网站照片| 在线看a的网站| 最近2019中文字幕mv第一页| 日日啪夜夜撸| 国产精品久久久久久精品古装| 亚洲熟女精品中文字幕| 日韩欧美精品免费久久| 小蜜桃在线观看免费完整版高清| 我要看日韩黄色一级片| 成年版毛片免费区| 2021少妇久久久久久久久久久| 亚洲精品自拍成人| 欧美日韩亚洲高清精品| 国产精品国产三级国产av玫瑰| 男人添女人高潮全过程视频| 18禁在线播放成人免费| 成年版毛片免费区| 免费观看的影片在线观看| 日韩人妻高清精品专区| 国产精品久久久久久精品古装| 亚洲怡红院男人天堂| 黄色欧美视频在线观看| 我要看日韩黄色一级片| 五月开心婷婷网| 国产亚洲精品久久久com| 成人毛片a级毛片在线播放| 国产精品一二三区在线看| 精品国产三级普通话版| 在线观看一区二区三区激情| 精品一区在线观看国产| 久久久久国产网址| 亚洲精品国产av成人精品| 精品人妻偷拍中文字幕| 国产视频内射| 黄片wwwwww| tube8黄色片| 亚洲精品国产成人久久av| 五月开心婷婷网| 国产老妇伦熟女老妇高清| 国产真实伦视频高清在线观看| 丝瓜视频免费看黄片| 一边亲一边摸免费视频| 国产又色又爽无遮挡免| 欧美变态另类bdsm刘玥| 精品久久国产蜜桃| 在线播放无遮挡| 毛片一级片免费看久久久久| 人体艺术视频欧美日本| 精品人妻熟女av久视频| 亚洲真实伦在线观看| 欧美成人午夜免费资源| 欧美精品国产亚洲| 久久久久久久久久成人| 精品久久久精品久久久| 日韩人妻高清精品专区| 国产精品蜜桃在线观看| 国产一区亚洲一区在线观看| 国产成人一区二区在线| 丰满少妇做爰视频| 男人舔奶头视频| 亚洲精品乱码久久久v下载方式| 大片电影免费在线观看免费| 欧美激情国产日韩精品一区| 性色avwww在线观看| 亚洲激情五月婷婷啪啪| 蜜桃久久精品国产亚洲av| 黑人高潮一二区|