• <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

    宅男免费午夜| 在线播放国产精品三级| 激情在线观看视频在线高清| 男人舔女人下体高潮全视频| 欧美性感艳星| 在线免费观看不下载黄p国产 | 网址你懂的国产日韩在线| 日本与韩国留学比较| 五月玫瑰六月丁香| www.www免费av| 免费黄网站久久成人精品 | 免费搜索国产男女视频| 亚洲美女搞黄在线观看 | 国产精品电影一区二区三区| 成人av一区二区三区在线看| 国产一区二区在线av高清观看| 亚洲欧美清纯卡通| 亚洲在线自拍视频| 国产一区二区三区视频了| 欧美午夜高清在线| 18美女黄网站色大片免费观看| 精品午夜福利视频在线观看一区| 亚洲国产精品合色在线| 午夜福利欧美成人| 日日干狠狠操夜夜爽| 色视频www国产| 亚洲中文日韩欧美视频| 在线观看av片永久免费下载| 欧美性猛交黑人性爽| 成人av在线播放网站| 精品人妻一区二区三区麻豆 | 99热这里只有精品一区| 国产午夜精品论理片| 97碰自拍视频| 少妇裸体淫交视频免费看高清| 国产中年淑女户外野战色| 欧美最黄视频在线播放免费| 日韩免费av在线播放| 日本精品一区二区三区蜜桃| 99riav亚洲国产免费| 色尼玛亚洲综合影院| 国产精品国产高清国产av| 国产真实伦视频高清在线观看 | 日本撒尿小便嘘嘘汇集6| 亚洲av中文字字幕乱码综合| 欧美又色又爽又黄视频| 少妇裸体淫交视频免费看高清| 夜夜夜夜夜久久久久| 看黄色毛片网站| 日本免费一区二区三区高清不卡| 变态另类丝袜制服| 国产中年淑女户外野战色| 麻豆一二三区av精品| ponron亚洲| 床上黄色一级片| 最近在线观看免费完整版| 亚洲av成人精品一区久久| 在线天堂最新版资源| 亚洲激情在线av| 国产精品久久视频播放| 日本在线视频免费播放| 变态另类成人亚洲欧美熟女| 丰满人妻一区二区三区视频av| 国模一区二区三区四区视频| 国产大屁股一区二区在线视频| 国产精品亚洲一级av第二区| 国产国拍精品亚洲av在线观看| 久久九九热精品免费| 中文在线观看免费www的网站| 欧美最新免费一区二区三区 | 欧美黄色淫秽网站| 成人高潮视频无遮挡免费网站| 亚洲国产精品sss在线观看| 欧美一区二区亚洲| 桃色一区二区三区在线观看| 黄色配什么色好看| 宅男免费午夜| 能在线免费观看的黄片| 国产伦在线观看视频一区| av在线观看视频网站免费| 舔av片在线| 亚洲欧美激情综合另类| 免费av观看视频| 午夜福利在线在线| 欧美日韩乱码在线| 国内久久婷婷六月综合欲色啪| 亚洲经典国产精华液单 | 免费观看精品视频网站| 成人国产一区最新在线观看| 国产伦人伦偷精品视频| 又爽又黄a免费视频| 又紧又爽又黄一区二区| 国产三级黄色录像| 精品无人区乱码1区二区| 欧美日本视频| 午夜福利18| 国产高清有码在线观看视频| 国产 一区 欧美 日韩| 一区二区三区激情视频| 国产成人福利小说| 99热这里只有是精品在线观看 | 性色avwww在线观看| 九色国产91popny在线| 亚洲精品在线美女| 国产成人福利小说| 国产精品久久久久久久久免 | or卡值多少钱| 校园春色视频在线观看| 精品一区二区三区视频在线观看免费| 国产精品爽爽va在线观看网站| 美女高潮的动态| 亚洲va日本ⅴa欧美va伊人久久| 国产一区二区在线观看日韩| 国产亚洲av嫩草精品影院| 窝窝影院91人妻| 中文字幕av成人在线电影| 亚洲性夜色夜夜综合| 99久国产av精品| 国产精品一区二区三区四区免费观看 | 国产成人a区在线观看| 悠悠久久av| 一个人免费在线观看的高清视频| 欧美日韩亚洲国产一区二区在线观看| 国产乱人伦免费视频| 免费无遮挡裸体视频| 99久久精品热视频| 免费观看的影片在线观看| 怎么达到女性高潮| 无遮挡黄片免费观看| 欧美在线一区亚洲| 人妻丰满熟妇av一区二区三区| 亚洲 国产 在线| 在线观看美女被高潮喷水网站 | 99国产极品粉嫩在线观看| 日本免费a在线| 别揉我奶头 嗯啊视频| 日本黄色视频三级网站网址| 免费看美女性在线毛片视频| 乱码一卡2卡4卡精品| 亚洲中文字幕一区二区三区有码在线看| 精品久久国产蜜桃| 午夜福利成人在线免费观看| 老熟妇仑乱视频hdxx| 亚洲最大成人av| 国产成人欧美在线观看| 999久久久精品免费观看国产| 成人特级av手机在线观看| 欧美丝袜亚洲另类 | 午夜a级毛片| 十八禁人妻一区二区| 丁香六月欧美| 亚洲欧美日韩无卡精品| 国产亚洲av嫩草精品影院| 成人精品一区二区免费| 在线观看一区二区三区| 麻豆国产av国片精品| 亚洲中文日韩欧美视频| 国产欧美日韩一区二区三| 99精品久久久久人妻精品| 国产亚洲精品av在线| 一a级毛片在线观看| 俺也久久电影网| 国产精品电影一区二区三区| 国产午夜精品久久久久久一区二区三区 | 变态另类成人亚洲欧美熟女| 变态另类丝袜制服| 久久伊人香网站| 夜夜看夜夜爽夜夜摸| 噜噜噜噜噜久久久久久91| 特大巨黑吊av在线直播| 99久久久亚洲精品蜜臀av| 听说在线观看完整版免费高清| 国产精品久久电影中文字幕| 婷婷精品国产亚洲av在线| 欧美激情在线99| 老女人水多毛片| 成人国产一区最新在线观看| 国模一区二区三区四区视频| 久久久久久久久久成人| 少妇的逼好多水| 成年免费大片在线观看| 全区人妻精品视频| 日本a在线网址| 真实男女啪啪啪动态图| 亚洲成人中文字幕在线播放| 成人av在线播放网站| 日韩国内少妇激情av| 一进一出好大好爽视频| 一进一出抽搐动态| 一级黄色大片毛片| 免费大片18禁| 午夜久久久久精精品| 亚洲专区中文字幕在线| 在线观看午夜福利视频| 成人永久免费在线观看视频| 在线免费观看的www视频| 欧美xxxx性猛交bbbb| 欧美日韩瑟瑟在线播放| 久久久久久久午夜电影| 蜜桃久久精品国产亚洲av| 亚洲成人久久性| 天美传媒精品一区二区| 国产色爽女视频免费观看| 日本五十路高清| 岛国在线免费视频观看| av专区在线播放| 91av网一区二区| 禁无遮挡网站| 国产精品不卡视频一区二区 | 99久久精品一区二区三区| 91av网一区二区| 在现免费观看毛片| 国产精品亚洲美女久久久| 精品一区二区免费观看| 国产69精品久久久久777片| 亚洲成人免费电影在线观看| av视频在线观看入口| 亚洲av成人精品一区久久| 午夜福利18| 亚洲最大成人手机在线| 亚洲精品一卡2卡三卡4卡5卡| 国产一区二区在线av高清观看| 欧美另类亚洲清纯唯美| 亚洲精品在线观看二区| 欧美色视频一区免费| 黄色配什么色好看| 综合色av麻豆| 亚洲在线自拍视频| 婷婷色综合大香蕉| 久久久久性生活片| 免费在线观看影片大全网站| 国产成人影院久久av| 亚洲,欧美,日韩| 亚洲国产欧洲综合997久久,| 国产久久久一区二区三区| 欧美+亚洲+日韩+国产| 久久午夜福利片| 在线观看舔阴道视频| 99热只有精品国产| 亚洲国产精品久久男人天堂| 无遮挡黄片免费观看| 久久久色成人| 免费无遮挡裸体视频| 午夜福利免费观看在线| 精品国内亚洲2022精品成人| 国产精品,欧美在线| 久久久国产成人精品二区| 亚洲中文字幕一区二区三区有码在线看| 亚洲av第一区精品v没综合| 欧美最黄视频在线播放免费| 久久中文看片网| 成年女人毛片免费观看观看9| 国产午夜福利久久久久久| 国产黄色小视频在线观看| 脱女人内裤的视频| 精品一区二区免费观看| 人妻久久中文字幕网| 久久热精品热| 麻豆国产av国片精品| 一个人看的www免费观看视频| 国产精品爽爽va在线观看网站| 最近最新中文字幕大全电影3| 美女大奶头视频| 国产人妻一区二区三区在| 三级男女做爰猛烈吃奶摸视频| 舔av片在线| 成人美女网站在线观看视频| 国产私拍福利视频在线观看| 国产久久久一区二区三区| 日本与韩国留学比较| 麻豆国产av国片精品| 亚洲精品一卡2卡三卡4卡5卡| 极品教师在线免费播放| 国产在线男女| 国产av不卡久久| 国产精品永久免费网站| eeuss影院久久| 国产麻豆成人av免费视频| 麻豆成人av在线观看| 深爱激情五月婷婷| 看黄色毛片网站| 一夜夜www| 亚洲av中文字字幕乱码综合| 天堂√8在线中文| www.熟女人妻精品国产| avwww免费| 老司机午夜十八禁免费视频| 一进一出抽搐gif免费好疼| www.www免费av| 欧美激情久久久久久爽电影| 欧美日韩乱码在线| 色综合婷婷激情| 国产高清视频在线观看网站| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 99精品久久久久人妻精品| h日本视频在线播放| 午夜福利免费观看在线| 精品久久久久久久人妻蜜臀av| 精品乱码久久久久久99久播| 亚洲久久久久久中文字幕| 日本五十路高清| 简卡轻食公司| 国产精品综合久久久久久久免费| 禁无遮挡网站| 亚洲电影在线观看av| 日本 欧美在线| 两性午夜刺激爽爽歪歪视频在线观看| 偷拍熟女少妇极品色| 国产伦精品一区二区三区四那| 国产私拍福利视频在线观看| 亚洲美女搞黄在线观看 | 午夜福利高清视频| 久9热在线精品视频| 日本成人三级电影网站| 午夜福利在线观看吧| 亚洲va日本ⅴa欧美va伊人久久| 淫秽高清视频在线观看| 亚洲av中文字字幕乱码综合| 欧美精品国产亚洲| 99热这里只有是精品50| 国产成年人精品一区二区| 国产成人啪精品午夜网站| 成人三级黄色视频| 一级av片app| 男人舔奶头视频| 成人性生交大片免费视频hd| 我的老师免费观看完整版| 国模一区二区三区四区视频| 欧美日韩综合久久久久久 | 麻豆国产av国片精品| 全区人妻精品视频| 久久性视频一级片| 亚洲 欧美 日韩 在线 免费| 高清毛片免费观看视频网站| 国产精品亚洲一级av第二区| 夜夜爽天天搞| av天堂中文字幕网| 最近最新免费中文字幕在线| .国产精品久久| 免费人成视频x8x8入口观看| 两性午夜刺激爽爽歪歪视频在线观看| 日本a在线网址| 国产精品99久久久久久久久| 免费观看人在逋| 国产精品人妻久久久久久| 床上黄色一级片| 人妻丰满熟妇av一区二区三区| 欧美黄色片欧美黄色片| 午夜精品一区二区三区免费看| 久久久久久久久大av| 色尼玛亚洲综合影院| 日本 av在线| 极品教师在线视频| 国产一区二区在线观看日韩| 国产精品久久电影中文字幕| 久久久成人免费电影| 中文字幕免费在线视频6| 深夜a级毛片| 欧美中文日本在线观看视频| 婷婷色综合大香蕉| 日韩欧美在线二视频| 中文字幕人成人乱码亚洲影| 亚洲精品乱码久久久v下载方式| 一a级毛片在线观看| 精品日产1卡2卡| 亚洲精品成人久久久久久| 免费看日本二区| 欧美日本亚洲视频在线播放| 首页视频小说图片口味搜索| 2021天堂中文幕一二区在线观| 国产中年淑女户外野战色| 国产三级中文精品| 欧美色欧美亚洲另类二区| 99国产极品粉嫩在线观看| 免费在线观看成人毛片| 2021天堂中文幕一二区在线观| 99久久精品国产亚洲精品| 精品人妻1区二区| 一级黄片播放器| 男人舔奶头视频| 国内精品久久久久久久电影| 精品久久久久久,| 成人鲁丝片一二三区免费| 亚洲成人精品中文字幕电影| 欧美成人免费av一区二区三区| 亚洲欧美日韩无卡精品| 一卡2卡三卡四卡精品乱码亚洲| 精华霜和精华液先用哪个| 亚洲精品456在线播放app | 狠狠狠狠99中文字幕| 国产探花极品一区二区| 一个人观看的视频www高清免费观看| 99久久九九国产精品国产免费| 免费搜索国产男女视频| 亚洲人成网站在线播放欧美日韩| 嫁个100分男人电影在线观看| 日本黄色视频三级网站网址| 国产真实乱freesex| 91久久精品电影网| 日本撒尿小便嘘嘘汇集6| 欧美日本亚洲视频在线播放| 成熟少妇高潮喷水视频| 亚洲综合色惰| 国产色爽女视频免费观看| 国产成年人精品一区二区| 中文字幕久久专区| 琪琪午夜伦伦电影理论片6080| 少妇的逼水好多| 亚洲av日韩精品久久久久久密| 免费av观看视频| 免费大片18禁| 亚洲色图av天堂| 久久天躁狠狠躁夜夜2o2o| 一级黄色大片毛片| 国产欧美日韩一区二区精品| 国产精品久久久久久亚洲av鲁大| 免费看日本二区| 久久精品夜夜夜夜夜久久蜜豆| 精品久久久久久久久久久久久| 亚洲综合色惰| 99精品在免费线老司机午夜| 综合色av麻豆| 少妇熟女aⅴ在线视频| 女人十人毛片免费观看3o分钟| 白带黄色成豆腐渣| 18禁裸乳无遮挡免费网站照片| 精品一区二区三区视频在线观看免费| bbb黄色大片| 午夜福利在线在线| 色av中文字幕| 色吧在线观看| 精品久久久久久久人妻蜜臀av| 欧美性猛交╳xxx乱大交人| 搞女人的毛片| 国产精品久久视频播放| 人妻夜夜爽99麻豆av| 可以在线观看的亚洲视频| 男人舔奶头视频| 成人鲁丝片一二三区免费| 午夜福利欧美成人| 九色国产91popny在线| 亚洲内射少妇av| 日韩欧美在线乱码| 欧美zozozo另类| 亚洲专区国产一区二区| 久久午夜亚洲精品久久| 亚洲精品在线观看二区| 亚洲久久久久久中文字幕| 国产私拍福利视频在线观看| a级一级毛片免费在线观看| 中文字幕熟女人妻在线| 人妻制服诱惑在线中文字幕| 99国产精品一区二区三区| 欧美成人性av电影在线观看| 在线天堂最新版资源| 国产真实乱freesex| av中文乱码字幕在线| 欧美黄色淫秽网站| 色综合亚洲欧美另类图片| 日本一本二区三区精品| or卡值多少钱| 白带黄色成豆腐渣| 1024手机看黄色片| 亚洲精品一区av在线观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲av免费在线观看| av专区在线播放| 超碰av人人做人人爽久久| 国产白丝娇喘喷水9色精品| 国产亚洲精品久久久com| ponron亚洲| 国产成人欧美在线观看| 我的女老师完整版在线观看| 亚洲av中文字字幕乱码综合| 日日干狠狠操夜夜爽| 国产不卡一卡二| 一本精品99久久精品77| 亚洲国产欧美人成| 国内精品一区二区在线观看| 亚洲成av人片免费观看| 亚洲精品亚洲一区二区| 欧美日韩亚洲国产一区二区在线观看| 欧美黑人欧美精品刺激| 亚洲午夜理论影院| 色播亚洲综合网| 国产午夜福利久久久久久| 亚洲欧美日韩高清专用| 村上凉子中文字幕在线| 国产精品永久免费网站| 国产精品日韩av在线免费观看| 久久精品国产亚洲av涩爱 | 日韩精品青青久久久久久| 国产v大片淫在线免费观看| 免费av毛片视频| 99久久精品一区二区三区| 嫩草影院新地址| 亚洲专区国产一区二区| 狂野欧美白嫩少妇大欣赏| eeuss影院久久| 国产成人a区在线观看| 乱人视频在线观看| 小说图片视频综合网站| 中文字幕精品亚洲无线码一区| 亚洲成人中文字幕在线播放| 亚洲在线观看片| 中文字幕熟女人妻在线| 一级黄色大片毛片| 国产成人av教育| 亚洲精品成人久久久久久| 国产在线精品亚洲第一网站| 欧美一区二区亚洲| 中文亚洲av片在线观看爽| 一级黄片播放器| 草草在线视频免费看| 成人欧美大片| 丝袜美腿在线中文| ponron亚洲| 欧美极品一区二区三区四区| 国产av不卡久久| 亚洲va日本ⅴa欧美va伊人久久| 日本一二三区视频观看| av在线观看视频网站免费| 黄色丝袜av网址大全| 精品久久国产蜜桃| 亚洲色图av天堂| 久99久视频精品免费| 久久精品国产自在天天线| 免费观看精品视频网站| 一区二区三区四区激情视频 | xxxwww97欧美| 免费电影在线观看免费观看| 国产熟女xx| 日韩欧美精品免费久久 | 精品久久久久久成人av| 国产精品三级大全| 国产白丝娇喘喷水9色精品| 能在线免费观看的黄片| 久久亚洲精品不卡| 赤兔流量卡办理| 久久亚洲真实| 此物有八面人人有两片| 直男gayav资源| a级毛片a级免费在线| 精品福利观看| 51国产日韩欧美| 九九久久精品国产亚洲av麻豆| 午夜视频国产福利| 精品久久久久久久末码| 99国产极品粉嫩在线观看| 天堂影院成人在线观看| 舔av片在线| 亚洲成人免费电影在线观看| 亚洲一区高清亚洲精品| www.999成人在线观看| 少妇熟女aⅴ在线视频| 亚洲国产色片| 国产主播在线观看一区二区| 人人妻人人看人人澡| 成人无遮挡网站| av中文乱码字幕在线| 最近最新中文字幕大全电影3| 深爱激情五月婷婷| 国产一区二区三区视频了| 亚洲欧美精品综合久久99| 少妇的逼好多水| 看十八女毛片水多多多| 美女 人体艺术 gogo| 高清日韩中文字幕在线| 搞女人的毛片| 国产真实乱freesex| 国内揄拍国产精品人妻在线| 看片在线看免费视频| 中国美女看黄片| 午夜影院日韩av| 欧美一级a爱片免费观看看| 日日干狠狠操夜夜爽| 国内少妇人妻偷人精品xxx网站| 久久伊人香网站| 国产精品精品国产色婷婷| 亚洲国产高清在线一区二区三| 中出人妻视频一区二区| 国产一区二区在线av高清观看| 国产精品1区2区在线观看.| 美女 人体艺术 gogo| 国产精品免费一区二区三区在线| 老司机深夜福利视频在线观看| 精品欧美国产一区二区三| avwww免费| 99久久99久久久精品蜜桃| 久久性视频一级片| 久久精品国产自在天天线| 少妇的逼水好多| 午夜免费激情av| 精品人妻视频免费看| 可以在线观看毛片的网站| 亚洲精品影视一区二区三区av| av天堂在线播放| 男女那种视频在线观看| 亚洲aⅴ乱码一区二区在线播放| 两个人的视频大全免费| 欧美国产日韩亚洲一区| 在线观看美女被高潮喷水网站 | 久久久久九九精品影院| 亚洲国产精品sss在线观看| av专区在线播放| 性色avwww在线观看| 精华霜和精华液先用哪个| 赤兔流量卡办理| 禁无遮挡网站| 一区福利在线观看| 色5月婷婷丁香| 别揉我奶头~嗯~啊~动态视频| 国产伦精品一区二区三区视频9| 国产伦在线观看视频一区| 国产人妻一区二区三区在|