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

    Electrochemical Removal of Chlorophenol Pollutants by Reactive Electrode Membranes: Scale-Up Strategy for Engineered Applications

    2022-04-22 11:42:38ShuzhoPeiYiWngShijieYouZhnguoLiNnqiRen
    Engineering 2022年2期

    Shuzho Pei, Yi Wng, Shijie You,*, Zhnguo Li, Nnqi Ren

    a State Key Laboratory of Urban Water Resource and Environment, School of Environment, Harbin Institute of Technology, Harbin 150090, China

    b State Key Laboratory of NBC Protection for Civilian, Beijing 102205, China

    Keywords:Chlorophenols Reactive electrode membrane Tubular concentric electrode Scale-up

    ABSTRACT Chlorophenols(CPs)are significant refractory pollutants that are highly toxic to humans and other organisms.Reactive electrode membranes(REMs)show considerable potential in the electrochemical removal of refractory pollutants by allowing flow-through operations with convection-enhanced mass transfer.However, relevant studies are commonly performed on the laboratory scale, and there is no straightforward method that guarantees success in scaling up engineered REM reactors. In this study, we demonstrated that a tubular concentric electrode (TCE) configuration with a titanium suboxide ceramic anode and a stainless-steel cathode is suitable for large-scale CPs removal. Both theoretical and experimental results showed that the TCE configuration not only allows the electrode surface to be orthogonal to electric field lines everywhere, but also has an ohmic resistance that is inversely proportional to the length of the electrode. In addition, the TCE configuration can be operated in either the anode-tocathode(AC)or the cathode-to-anode (CA)mode based on the flow direction,creating adjustable conditions for selective degradation of CPs. This was confirmed by 98% removal of 2,4-dichlorophenol (2,4-DCP)and 72.5%removal of chemical oxygen demand(COD)in the CA mode,in which the kinetic constant was one order of magnitude higher than that for the AC mode under flow-through single-pass operations.This can be explained by the lower activation energy and free energy in the CA mode,as revealed by theoretical calculations and experimental measurements. The TCE configuration is also suitable for a numbering-up strategy to scale up the electrochemical reactor without increasing the ohmic resistance or decreasing the specific electrode area, achieving 99.4% removal of 2,4-DCP with an energy consumption of 1.5 kW·h·m-3 when three TCE modules were employed. This study presents a suitable electrode design configuration for the REM reactor, offering effective strategies to bridge the ‘‘Valley of Death”encountered when scaling up the electrochemical removal of CP pollutants.

    1. Introduction

    Chlorophenols (CPs) are quantitatively significant pollutants released into the environment by chemical and pharmaceutical industry activities. They are also by-products of industrial processes in the manufacturing of polymers, dyes, textiles, resins,and explosives. CPs and their derivatives are mutagenic, carcinogenic, and immunogenic for humans and other organisms, and their treatment, disposal, and management have become serious challenges to stakeholders and health sectors. Most types of CPs are refractory to microbial degradation, making it imperative to develop more efficient methods for minimizing their pollution in water environments and ecological systems. Electrochemical oxidation has drawn increasing research interest for decentralized water decontamination because the application of electrons as traceless reactants offers several unique advantages, such as high efficiency,cleanliness,process integration,and ease of automation[1-3].Electrochemical oxidation leads to the destruction of CP pollutants via direct electron transfer (DET) and indirect oxidation mediated by·OH produced by water electrolysis.

    Because low-concentration CPs can be oxidized faster than they can be transported to the anode surface,electrochemical oxidation of CPs is commonly controlled by diffusion. Reactive electrode membranes (REMs) offer a promising alternative for improving mass-transfer efficiency by flow-through electrolysis [4,5].Compared with mechanical agitation and flow-by operations,REMs can conduct electrolysis in the flow-through mode with convection-enhanced mass transfer [6]. Thus, membrane electrodes show great promise for potential applications in water purification, such as pollutant removal, self-cleaning membrane,resource recovery, and seawater desalination [7-9].

    Most studies on REMs are performed at the laboratory scale under well-controlled conditions. Consequently, a gap remains between experimental studies and practical applications. The process of scaling up a reactor from a lab prototype to an industrial plant is often referred to as the ‘‘Valley of Death,” and there is no straightforward method that guarantees success in the scale-up of membrane electrodes. Unlike conventional chemical and biochemical reactors, where the geometric, kinematic, and thermal features are similar between the model and practical units, the electrode configuration of REMs is of prime importance [2,10,11].Most flow-through electrochemical reactors are configured for plug-flow filtration [12,13] based on two parallel plate electrodes(PPEs) made of porous conductive materials, such as carbon nanotubes(CNTs)or Ti4O7[14,15].Although this mode provides a convenient and easy-to-handle platform to test the performance of electrode materials, challenges remain to scale up the electrochemical setup because increasing the electrode spacing or reactor volume increases the ohmic resistance and decreases the area-tovolume ratio. Hence, there is a compelling need to develop more efficient, reliable, and scalable electrode configurations to enable well-engineered industrial applications.

    In principle, the appropriate electrode configuration should①have a reasonable electric field distribution and low ohmic resistance, ②be effective for the removal of CP pollutants, and③allow modularization for easy-to-handle scale-up and engineering applications. Tubular electrodes with a cylindrical geometry are more favorable for increasing current density and power generation in energy conversion systems, such as solid oxide fuel cells(SOFCs), polymer electrolyte membrane fuel cells (PEMFCs), and vanadium redox flow batteries (VRFBs) [13,16-18]. Inspired by these successful applications, we attempted to bridge the ‘‘Valley of Death” for scaling up REMs by developing a tubular concentric electrode (TCE) configuration for the engineered application of electrochemical CPs removal.

    In this study, we developed a TCE configuration containing an inner titanium suboxide (TiSO) tube (ceramic membrane anode)and outer stainless-steel (SS) tube (cathode), which allows the wastewater stream to flow through the two tubes in either the anode-to-cathode (AC) or cathode-to-anode (CA) direction. The tubular TiSO was selected as the anode material because its conductivity approaches that of metals, and its corrosion resistance is close to that of ceramics.The unique advantages of high conductivity,good stability,high oxygen evolution potential,and low cost make it a promising anode material for water purification [19].First, we investigated the electric field distribution for different electrode configurations using software simulation and electrochemical validation. Second, we tested the electrochemical removal of typical biorefractory CP pollutants (i.e., 2,4-dichlorophenol (2,4-DCP)) under AC and CA flow modes, and the relevant mechanisms were elucidated using theoretical calculations and experimental data. Finally, we proposed a numberingup strategy for scaling up the TCE module and flow-through reactor and tested its suitability for engineered water decontamination in single-pass mode. The primary objective of this study was to determine a reasonable electrode configuration for the scale-up of the REM reactor for CP pollutants removal.

    2. Materials and methods

    2.1. Materials

    All purchased reagents were of analytical grade if there is no special statement and used without further purification. Sodium perchlorate(NaClO4),2,4-DCP,and chromatographic-grade methanol were purchased from Sinopharm Chemical Reagent Co., Ltd.(China), and the solution was prepared using ultrapure water obtained from a Milli-Q system (18.2 MΩ·cm at (22 ± 1) °C). The TiSO electrode was designed and manufactured according to the procedures described in our previous study [19].

    2.2. Experimental setup and operation

    The electrolysis experiments were conducted in a Plexiglas cylindrical electrolytic cell with a total volume of 3.8 L(220 mm×φ 150 mm)at room temperature((25±1)°C).The cell contained a tubular TiSO anode(outer diameter,30 mm;thickness,3 mm; length, 143 mm; projected area, 132.9 cm2) and a SS cathode (140 cm2). The electrodes were connected to a direct current(DC) power supply by copper wires. Prior to the tests, the TiSO electrodes were cleaned for 10 min in a 0.5 mol·L-1NaClO4solution at a current density of 30 mA·cm-2. The initial influent contained 0.1 mol·L-1NaClO4and 10 mg·L-12,4-DCP at a neutral pH(7.2 ± 0.2), and the electrolytic cell was operated in single-pass mode without recirculation. The experiments were performed with 10 min of stabilization without the application of current,and 5 mL aliquots were withdrawn from the effluent at specified intervals, followed by immediate filtration with cellulose acetate filters (0.22 μm) before chemical analysis.

    The hydraulic retention time (HRT, s) of the system was calculated by

    where J (m·s-1) is the permeation flux, d (m) is the thickness of a single TiSO anode, A (m2) is the anode area, and Q (m3·s-1) is the flow rate.

    2.3. Electrochemical characterization

    Electrochemical characterization was performed with a SS counter electrode and Ag/AgCl (saturated KCl aqueous solution(KClaq), 0.198 V vs standard hydrogen electrode (SHE)) reference electrode.Cyclic voltammetry(CV)was performed at a sweep rate of 100 mV·s-1in a three-electrode cell in the 2,4-DCP degradation phase. Electrochemical impedance spectroscopy (EIS) was conducted between 10 mHz and 100 kHz using a sinusoidal potential perturbation with an amplitude of 5 mV.

    2.4. Analyses

    The 2,4-DCP concentration was determined using a ultrahighperformance liquid chromatography (UPLC) instrument (Agilent 1200, USA) equipped with a photodiode array detector (PDA) and a reversed-phase C18 column (50 mm × φ 2.1 mm, Waters,USA). The mobile phase consisted of water/methanol (30:70, v/v)at a flow rate of 0.1 mL·min-1. The degraded intermediates of 2,4-DCP were analyzed qualitatively based on the mass-to-charge ratio using UPLC with electrospray ionization mass spectrometry(UPLC-ESI/MS)(Agilent 6460,USA)with mobile phase A(ultrapure water) and B (chromatographic-grade acetonitrile). Highresolution mass spectra were collected for range of the mass-tocharge ratio (m/z) 80-220 using negative continuum mode with a scanning time of 0.15 s, and the suspected intermediates were fragmented in the collision cell. The ionization source was set up with a capillary voltage of 3 kV, temperature of 200 °C, and desolvation temperature of 550 °C. The samples (20 μL) were analyzed in duplicate to ensure the accuracy and repeatability of the mass spectroscopy measurements. Chemical oxygen demand (COD)was determined using an LH-3BA spectrophotometer (Lianhua Technology Co., Ltd. (China), chemical reagent type number: LHD/E).

    Energy consumption normalized per order(EEO)of contaminant was calculated by Ref. [20]:

    where CODf(mg·L-1) and CODp(mg·L-1) are the COD concentrations of the feed and permeate, respectively, b (4) is the stoichiometric coefficient for number of electron transfer, F is the Faraday constant(96 485 C·mol-1),V(L)is the volume treated over duration t, and MO2is the equivalent molar mass of oxygen (32 g·mol-1).

    2.5. Electric field simulation

    Simulations were conducted on the PPE, tubular/planar electrode (TPE), and TCE configurations using COMSOL Multiphysics software (version 5.5; Palo Alto, USA). The models were used to simulate the electric field distribution for different electrode configurations using the definition of electric field strength:

    2.6. Theoretical calculations

    Density functional theory (DFT) calculations were conducted using Gaussian software. The activation energy (Ea, kJ·mol-1) for DET oxidation as a function of electrode potential was calculated according to Marcus theory. The calculation for the DET reduction reaction is similar and not shown.

    3. Results and discussion

    3.1. Electric field distributions of different electrode configurations

    The spatial distribution of the electric field significantly affects electrochemical reactions because it is closely linked to charged ion transport[13,22-25].In other words,electrochemical reactions only occur in the region of the electrode surface where the electric field lines reach,and vice versa.Because the electric field distribution is closely related to the electrode configuration,we used COMSOL modeling to visualize and quantitatively analyze the electrostatic field of three typical electrode configurations for wastewater treatment [13]: PPE, TPE, and TCE. As shown in Fig. 1, the symmetric PPE and TCE configurations exhibited the most uniform electric field distributions, whereas an asymmetric electric field was observed for the TPE configuration,which is similar to the results found by Hankin et al. [23] in a photoelectrochemical reactor and Sun et al. [26] in an electrochemical reactor[26]. This clearly suggests that complete utilization of the anode area could be achieved with the PPE and TPE configurations.However, less than 1/3 of the area of the tubular anode was active for the TPE configuration because the equipotential lines are always perpendicular to the electric field at each of their crossings. The electric field points away from the anode plane and toward the cathode plane. Because the charges are equal and opposite, the electric fields sum in the interpolar region, whereas they cancel each other out to zero in the region outside of the two electrodes[26]. In other words, only the electrode surfaces normal to each other are electrochemically active, and no transformation occurs in other regions. This can be easily visualized experimentally by observing the electrodeposition of a Cu film on only the front side of an SS disc(Fig.S1 in Appendix A).Note that the small amount of Cu film deposited at the edges of the back side was caused by the molecular diffusion of Cu2+under a concentration gradient.

    For the same anode, the cathode geometry can have a significant impact on the electric field distribution. Clearly, a uniform interelectrode electric field distribution is essential to effectively utilize the electrochemically active area because the ions can diffuse along the field lines, which has been proven to be efficient for water disinfection [24]. These results suggest that when designing a flow-through electrolytic cell, the surface of the electrode should be orthogonal to the electric field lines everywhere to maximize the effective electrochemically active area of the electrode. In accordance with this principle, the PPE design is more suitable for plug-flow electrolysis,whereas the TCE design is more rational for encompassed-flow electrolysis.

    Fig. 1. COMSOL simulations of the electrostatic field for the (a) PPE, (b) TPE, and(c) TCE configurations with front and sectional views; (d) optical image of the TCE module.

    3.2. Dependence of ohmic resistance on electrode configuration

    When the anode is polarized by passing current through the electrolytic cell, the voltage drop U (V) is determined by the equilibrium potential difference and overpotential [27]:

    where Umin(V) is the theoretical minimum voltage required for electrolysis, ηA(V) is the activation overpotential dependent on a suitable catalyst,ηD(V)is the diffusion overpotential due to concentration polarization,and ηohm(V)is the ohmic overpotential caused by the electrode material, configuration, and electrolyte [27]. For a given electrolytic system under flow-through operation,ηohmis the dominant contributor to the cell voltage and thus energy consumption under constant-current electrolysis.Herein,the ohmic overpotential is calculated according to Ohm’s law:

    where ΣRohm(Ω) is the sum of the ohmic resistances.

    It is worth noting that the electrode configuration may have a significant impact on the magnitude of the ohmic resistance. For example, based on an identical anode diameter (29 mm), axial length (142 mm), and electrolyte (0.1 mol·L-1NaClO4, pH 6.5),the ΣRohmvalue of the TCE configuration (2.1 Ω) was 89.3% and 83.3% lower than those of the PPE (19.8 Ω) and TPE (12.6 Ω) configurations, respectively, according to the x-intercept in the Nyquist plot of the EIS data (Fig. 2(a)). This difference could result from the difference in ion transport in the electrolyte in the interpolar space for the two systems. For the PPE configuration, the ohmic resistance is positively proportional to the axial length of the reactor, as follows:

    where r1is the radius of the tubular anode, and r2is the radius of the tubular cathode.

    Fig.2. EIS spectra for(a)different electrode configurations and(b)PPE and(c)TCE tubular reactors with various axial lengths; (d) total internal resistance (ΣRohm) as function of TCE length.

    Despite minor differences in ΣRohmvalues between the PPE and TCE configurations in small lab-scale prototype cells, scaling up these two configurations would lead to opposite behaviors of their ΣRohmvalues with increasing reactor volume (Figs. 2(b) and (c)).The decreased ohmic resistance of the TCE system becomes more significant when the conductivity of the solution is low (e.g.,domestic sewage and groundwater) [28]. As shown in Fig. 2(d),the calculated ΣRohmvalue of the PPE configuration was 16 times that of the TCE configuration when the axial length of the tubular reactor was increased from 1 to 14 cm(Eqs.(9)and(10)).The TCE configuration was compact with a small electrode spacing and large electroactive area, considerably enhancing electron transfer and lowering the internal resistance [29].

    Together,the combination of experimental results and theoretical verification demonstrate that the TCE configuration is more favorable for scaling up REMs owing to its uniform electric field and low ohmic resistance. Similarly, Lei et al. [30] reported efficient phosphate recovery (>50 %) using a scalable prototype with a tubular SS cathode in a columnar electrochemical reactor over 173 days. The unique advantages of the amplified TCE configuration enabled the application of a scaled-up membrane electrode system[31].Next,we shifted our focus to the engineered application potential of the TCE system for wastewater treatment.

    3.3. Electrochemical removal of 2,4-DCP by the REM module

    The electric field between the anode and cathode in the TCE configuration exhibited a uniform distribution, suggesting that efficient faradaic reactions occurred. The stream flow sequence,upstream from the anode to the cathode(AC mode)or downstream from the cathode to the anode (CA mode), has a crucial impact on the electrochemical reactions. Thus, we examined a tubular concentric TiSO anode/SS cathode REM reactor for the electrochemical removal of 2,4-DCP(10 mg·L-1)with a NaClO4(0.1 mol·L-1,pH 6.5)electrolyte under flow-through operation with different flow directions, current densities, and flow rates. Electrocatalytic reduction of 2,4-DCP occurs by accepting electrons at the cathode, whereas oxidation reactions occur at the TiSO anode. The anodic oxidation and cathodic reduction are coupled by the hydraulic connection of the flow stream between the electrodes.

    A series of experiments were conducted to investigate the removal of 2,4-DCP at a current density of 5 to 30 mA·cm-2and a flow rate of 15 to 120 mL·min-1.As shown in Fig.3(a),increasing the current density from 5 to 30 mA·cm-2increased 2,4-DCP removal from 5%-28% to 18%-71% in AC mode. Because the limiting current density (0.6 mA·cm-2) for the anodic oxidation of 10 mg·L-12,4-DCP was considerably lower than the applied current density, 2,4-DCP was degraded under diffusion-controlled conditions. Nevertheless, when the flow velocity was increased in the range 15-120 mL·min-1, the mass transfer coefficient increased from 2.5×10-6to 2×10-5m·s-1.The decrease in retention time from 160 to 20 s under continuous operation resulted in a 70.7%-78.5% decrease in 2,4-DCP removal under AC mode. Similar results were also reported by Wang et al. [32], who found enhanced electrochemical oxidation efficiency at lower flow rates and higher applied current densities. Similarly, for CA mode, 2,4-DCP removal increased from 30%-85% (5 mA·cm-2) to 81%-99%(30 mA·cm-2)with increasing current density and decreased from 85%-99%(15 mL·min-1)to 30%-81%(120 mL·min-1)with increasing flow rate. The rate constant obtained under CA mode was one order of magnitude greater than that under AC mode(Figs.3(c)and(d)),indicating the inherent dependence of synergistic cathode/anode interactions on the flow direction. Fig. S2 in Appendix A compares 2,4-DCP removal under the two operational modes at 20 mA·cm-2and 60 mL·min-1. Under CA mode, almost complete removal of 2,4-DCP was achieved after 5 min of electrolysis, compared with only 16.8% removal under AC mode. Based on a flow rate of 60 mL·min-1, raising the current density from 5 to 30 mA·cm-2increased the rate constant from 0.19 to 0.21 min-1in CA mode, which is approximately 11.2-78.1 times that in AC mode (Fig. 3(c)). When the current density was fixed at 20 mA·cm-2, increasing the flow rate from 15 to 120 mL·min-1enhanced the oxidation rate because of improved flow-through convective diffusion in the porous structure of the anode. This accounts for the increase in the rate constant from 0.13 to 0.28 min-1in CA mode, which is approximately 6.7-11.2 times that in AC mode(Fig.3(d)).The kinetic behavior of 2,4-DCP oxidation in CA mode appeared to be more sensitive to the flow rate, as indicated by the nearly one order of magnitude higher areanormalized apparent rate constant (0.036-0.079 m·s-1) compared to that for AC mode (Table S1 in Appendix A). This implies that·OH-mediated oxidation was predominantly limited by mass transfer following cathodic activation (discussed below) when the wastewater stream flowed from the cathode to the anode.

    Fig. 3. 2,4-DCP removal performance of the TCE reactor as a function of flow direction, current density, and flow rate. Removal efficiency in (a) AC and (b) CA modes as a function of current density and flow rate.Removal kinetics as a function of (c) current density at a flow rate of 60 mL·min-1 and (d) flow rate at a current density of 20 mA·cm-2 for both flow modes(Kobs is the observed rate constant).(e)COD removal and (f)current efficiency in AC mode(hollow symbols) and CA mode(solid symbols) as a function of current density and flow rate (AC: waterstream from anode to cathode; CA: waterstream from cathode to anode).Influent solution([2,4-DCP] = 10 mg·L-1, [NaClO4] = 100 mmol·L-1) was electrolyzed.

    In addition,water electrolysis resulted in the release of H+near the anode and OH-near the cathode, and the difference between CA and AC modes was also associated with differences in the migration of H+or OH-, particularly at high current densities. In CA mode, OH-produced from the cathodic hydrogen evolution reaction (HER) diffused to the anode, resulting in a slight increase in the pH value near the anode due to neutralization of protons produced from water oxidation(Fig.S3 in Appendix A).In this case,deprotonated 2,4-DCP (pKa= 7.8) was more likely to be adsorbed on the anode via electrostatic interactions, which facilitated the subsequent anodic oxidation. As shown in Fig. S4 in Appendix A,the CV curve exhibited a minimum potential of 1.66 V vs SHE for 2,4-DCP oxidation, and the intensity of the oxidative peaks was sensitive and positively related to the pH value.This indicates that 2,4-DCP oxidation is more favorable under higher pH conditions.Similarly, the typical oxidative intermediates, such as 2-CP and 4-CP, had pKavalues of 8 to 9 [33], making them more reactive at higher pH values. Together, these factors allowed the TCE reactor to degrade 2,4-DCP with a maximum overall COD removal of 72.5% (Fig. 3(e)) and a current efficiency of 35.5% (Fig. 3(f)) under CA operation,which are considerably higher than the values under AC operation.

    3.4. Mechanistic insight into the impact of hydraulic flow

    To further elucidate the 2,4-DCP degradation mechanism of the TiSO-based TCE reactor, the intermediate products were analyzed using UPLC-quadrupole-time-of-flight (UPLC-Q-TOF) mass spectrometry based on the mass-to-charge ratio (m/z). For AC mode,the results (Fig. S5 in Appendix A) matched those reported in previous studies, where six types of typical intermediates were detected during the electrochemical oxidation of 2,4-DCP [34,35].In contrast, different intermediate products were detected for CA mode, indicating the reaction pathway and selectivity depend on the reaction sequence determined by the flow direction. Fig. S6 in Appendix A shows the accumulation of chlorinated dihydroxybenzenes(2 or 3)(m/z=142.9342)and chlorinated benzoquinones(4 or 5) (m/z = 140.9376), which were characterized as the main intermediates in AC mode (Fig. 4). The reactions involve direct and indirect anodic oxidation of 2,4-DCP and direct cathodic reduction of chlorinated benzoquinone.In comparison,the characteristic intermediates observed for CA mode were phenol (m/z = 93.1581)and CPs (8 or 9) (m/z = 126.7408), whose relative abundance increased within the first 2 min,followed by a rapid decline without noticeable accumulation. Therefore, the cathodic reduction of 2,4-DCP and the anodic oxidation of phenol (10), hydroquinone(11), or catechol (12) are the most probable relevant reactions.

    The mechanistic insight into the impact of the flow mode can also be explained in the context of the difference in free energy for the two operational modes predicted by DFT calculations(Fig. 4). As shown in Fig. S7 in Appendix A, the activation energy for electro-reduction was lower than that for electro-oxidation.Specifically,the minimum potentials required for anodic oxidation of 2,4-DCP,phenol(10),hydroquinone(11),and catechol(12)were 1.90,1.85,1.40,and 1.50 V vs SHE,respectively,whereas the cathodic reduction of CPs (8 or 9), chlorinated benzoquinones (4 or 5)occurred at maximum potentials of -1.55, -1.55, -0.55, and-0.75 V vs SHE,respectively.Considering the similar redox nature of the DET of dihydroxybenzenes(11 or 12)and benzoquinones(6 or 7), the calculation of the overall free energy was simplified based on 1,4-benzoquinone (7). As revealed in Table S2 in Appendix A, the free energy was calculated to be -488.4 and -684.7 kJ·mol-1for AC and CA modes, respectively, at anode and cathode potentials of 2.38 and -1.45 V vs SHE, respectively, indicating a more thermodynamically favorable pathway for 2,4-DCP degradation under CA mode. As shown in Fig. 4(c), the conversion of ochloro-dihydroxybenzene (3) to o-chloro-p-benzoquinone (5) was the limiting step for AC mode,which could be alleviated by implementing cathodic reduction before anodic oxidation.Owing to the considerably lower energy required for the dechlorination of CPs(8 or 9) (-250.4 kJ·mol-1), cathodic reduction allows the more thermodynamically feasible dechlorination process to occur under a low applied cathode potential(E0= -1.15 V vs SHE).The formation of phenol as a key intermediate increases the activity for subsequent anodic DET oxidation (Fig. S7) and·OH attacks, which provides the most likely explanation for the high activity and high selectivity.Compared with that in AC mode,a high electric energy input is not required because cathodic reduction greatly decreases the minimum potential and activation energy (Uminin Eq. (7)) of 2,4-DCP,resulting in enhanced current efficiency and electrochemical mineralization (Fig. 3). Another advantage of CA operation is the lower biotoxicity of the dechlorinated intermediates (Fig. S8 in Appendix A). In this way, we were able to effectively destroy more CP pollutants with the removal efficiency more than 85%(Fig. S9 in Appendix A; 20 mA·cm-2and 60 mL·min-1).

    3.5. Numbering-up strategy for scaling up the REM module

    Based on the above results, we attempted to scale up the flowthrough electrolytic reactor based on the TCE configuration. In an electrochemical system, the area-to-volume ratio (A/V), that is,the specific electrode area(AS), is the key parameter that determines the inherent conversion capacity of the electrochemical system[36]. The importance of ASmust be emphasized when the REM is scaled up to a larger size.

    If the electrochemical system is scaled up along the axial length of the tubular reactor(PPE),the size and volume of the reactor will increase at the expense of a decline in AS:

    where L is the axial length of the electrolytic cell,and r is the radius of the planar electrode.Although it is theoretically possible to scale up a PPE reactor by increasing the electrode area,this not only adds complexity to electrode fabrication and process integration but also causes nonuniformity in the primary current distribution on the electrode [37].

    Because the TCE configuration allows water to flow through the side wall of the anode, the cylindrical geometry provides a higher surface-to-volume ratio than that of the planar geometry. Scaling up the electrolytic cell can be achieved simply by numbering up the TCE arrays, in which ASis independent of the axial length and volume of the reactor:

    where r1is the radius of the inner tube of a single TCE module,R is the radius of the reactor, and n is the number of TCE modules. This design enables modularization of TCE arrays with arbitrary axial lengths and numbers to create integrated electrolytic reactors that depend on the amount of wastewater to be treated in engineered applications.

    Fig.4. Mechanism under(a)AC and(b)CA mode,and(c)DFT calculation for 2,4-DCP degradation at anode and cathode potentials of 2.38 and-1.45 V vs SHE,respectively.

    To evaluate the effectiveness of the numbering-up strategy,we designed an electrochemical reactor containing TiSO TCE arrays with n = 1, 2, and 3 (Fig. 5(a)) and tested the overall performance of 2,4-DCP degradation(10 mg·L-1).Increasing the n value resulted in a higher removal efficiency,higher current efficiency,and lower energy consumption due to the increased area-to-volume ratio,particularly in AC mode,which offered slower kinetics.At the same current intensity(4 A),increasing n from 1 to 3 caused the removal efficiency to increase from 27.3% to 51.6% in AC mode and from 92.5%to 99.4%in CA mode(Fig.5(b))and the energy consumption per order(EEO)to decrease from 24.6 to 10.8 kW·h·m-3in AC mode and from 3.0 to 1.5 kW·h·m-3in CA mode(Fig.5(c)).Moreover,the feasibility of the numbering-up strategy for scale-up was confirmed by the stable 2,4-DCP removal efficiency(Fig.S10 in Appendix A).Despite the similar cell voltages(4.7,3.8,and 3.5 V for n=1,2,and 3,respectively)for the two modes,the energy consumption in CA mode was nearly one order of magnitude greater than that in AC mode because the decisive factor was effluent residual concentration(Cp;Eq.(2))rather than the removal efficiency.The considerably higher energy consumption could be interpreted in the context of a higher Cpdue to less efficient destruction of organic pollutants in AC mode. This corresponds to the lower current efficiency and COD removal in AC mode (68.7%-87.0% and 64.2%-87.9% lower, respectively) than those in CA mode (Fig. S11 in Appendix A). The lowest EEOobtained at 1.5 kW·h·m-3(99.4%removal)is in the range of values for the advanced electrochemical oxidation process of the membrane electrode on a laboratory scale(<2.0 kW·h·m-3)[20].Compared with conventional electrochemical treatment (summarized in Table S3 in Appendix A), a higher energy consumption (several orders of kilowatt-hours per cubic meter) (Table S3 [37,38]) was the consequence of the higher thermodynamic potential necessary for·OH production(2.38 V vs SHE)and side parasitic reactions, such as oxygen evolution by water electrolysis. Despite this, the TCE-based modularized process is still competitive for larger-scale decentralized treatment of highsalinity industrial wastewater.

    Fig. 5. (a) Optical images of the modularized TCE electrolytic reactor (n = 1, 2, 3);(b) removal efficiency of 2,4-DCP and (c) energy consumption per order under AC and CA modes.Tests were performed at[2,4-DCP]of 10 mg·L-1,total current of 4 A,and flow rate of 60 mL·min-1.

    4. Conclusions

    Based on the above results, the following conclusions can be drawn. The TCE configuration with a titanium suboxide ceramic anode and stainless-steel cathode was shown to be suitable for scale-up of the reactor.In the TCE configuration,the electrode surface is orthogonal to electric field lines everywhere,and the ohmic resistance is inversely proportional to the length of the electrode.The TCE configuration can be operated in either AC or CA mode based on the flow direction, creating adjustable conditions for selective degradation of pollutants. In CA mode, 98% removal of 2,4-DCP and 72.5% removal of COD were achieved because the kinetic constant was one order of magnitude higher than that in AC mode under flow-through single-pass operation. The TCE configuration was also shown to be suitable for a numbering-up strategy to scale up the electrochemical reactor without increasing the ohmic resistance or reducing the specific electrode area, which achieved 99.4% removal of 2,4-DCP with an energy consumption of 1.5 kW·h·m-3when three TCE modules were employed. This study provides a reasonable, easy-to-handle, and reliable strategy for scaling up REM reactors,which can also be extended to realize the removal of other recalcitrant organic pollutants and decentralized electrochemical wastewater treatment.

    Acknowledgments

    This work was supported by the National Natural Science Foundation of China(U21A20161 and 51822806), State Key Laboratory of Urban Water Resource and Environment,China(Harbin Institute of Technology) (2020DX07), and Heilongjiang Touyan Innovation Team Program, China (HIT-SE-01).

    Compliance with ethics guidelines

    Shuzhao Pei, Yi Wang, Shijie You, Zhanguo Li, and Nanqi Ren declare that they have no conflict of interest or financial conflicts to disclose.

    Appendix A. Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.eng.2021.11.017.

    国产女主播在线喷水免费视频网站 | 自拍偷自拍亚洲精品老妇| 99热这里只有是精品50| 国产美女午夜福利| 亚洲最大成人av| 免费av毛片视频| 97人妻精品一区二区三区麻豆| 日韩欧美 国产精品| 亚洲欧美成人精品一区二区| 亚洲三级黄色毛片| 国产v大片淫在线免费观看| 亚洲国产色片| 国内精品宾馆在线| 国产成人a∨麻豆精品| 免费大片18禁| 国产精品日韩av在线免费观看| 九九久久精品国产亚洲av麻豆| 亚洲精品成人久久久久久| 国产极品天堂在线| 久久精品91蜜桃| 中文在线观看免费www的网站| 国产单亲对白刺激| 午夜激情欧美在线| 亚洲最大成人手机在线| 春色校园在线视频观看| 亚洲av一区综合| 国产精品日韩av在线免费观看| 精品久久国产蜜桃| 男女边吃奶边做爰视频| 亚洲18禁久久av| 亚洲av成人av| 日日撸夜夜添| 日韩欧美一区二区三区在线观看| 黄片无遮挡物在线观看| 黄片无遮挡物在线观看| 久久99蜜桃精品久久| 亚洲激情五月婷婷啪啪| 日韩精品青青久久久久久| 神马国产精品三级电影在线观看| 欧美+亚洲+日韩+国产| 日本黄色视频三级网站网址| 黄色日韩在线| 精品不卡国产一区二区三区| 免费无遮挡裸体视频| 在线观看66精品国产| 九九在线视频观看精品| 九九在线视频观看精品| 1000部很黄的大片| 男的添女的下面高潮视频| 久久精品91蜜桃| 免费无遮挡裸体视频| 日本熟妇午夜| 九草在线视频观看| 美女脱内裤让男人舔精品视频 | 寂寞人妻少妇视频99o| 国产爱豆传媒在线观看| 一进一出抽搐gif免费好疼| 一级黄色大片毛片| 黄色日韩在线| 老熟妇乱子伦视频在线观看| 色综合色国产| 亚洲性久久影院| 欧美激情国产日韩精品一区| 伦理电影大哥的女人| 夜夜看夜夜爽夜夜摸| 亚洲精品自拍成人| 欧美精品一区二区大全| 久久久午夜欧美精品| 国产极品天堂在线| 亚洲精品久久久久久婷婷小说 | 高清毛片免费看| 一本精品99久久精品77| 老女人水多毛片| 在线播放国产精品三级| 真实男女啪啪啪动态图| 国产高清三级在线| 天堂av国产一区二区熟女人妻| av专区在线播放| 国产精品人妻久久久久久| 少妇熟女aⅴ在线视频| 欧美不卡视频在线免费观看| 在线观看免费视频日本深夜| 两个人的视频大全免费| 久久欧美精品欧美久久欧美| 久久欧美精品欧美久久欧美| av免费观看日本| 欧美bdsm另类| 免费不卡的大黄色大毛片视频在线观看 | or卡值多少钱| 一级毛片久久久久久久久女| 一个人观看的视频www高清免费观看| 男插女下体视频免费在线播放| av女优亚洲男人天堂| 国产午夜精品一二区理论片| 18+在线观看网站| 国产在线男女| 白带黄色成豆腐渣| www.av在线官网国产| 黄色配什么色好看| 99久久九九国产精品国产免费| 白带黄色成豆腐渣| 免费观看在线日韩| 日本免费一区二区三区高清不卡| 中文字幕av成人在线电影| 人人妻人人澡欧美一区二区| 亚洲精品久久久久久婷婷小说 | 亚洲婷婷狠狠爱综合网| 成人av在线播放网站| 在线免费观看不下载黄p国产| 亚洲精品久久国产高清桃花| 国内久久婷婷六月综合欲色啪| 天堂√8在线中文| 三级经典国产精品| 99热网站在线观看| 亚洲精品亚洲一区二区| 欧美又色又爽又黄视频| 精品一区二区三区视频在线| 熟女人妻精品中文字幕| 日本撒尿小便嘘嘘汇集6| 99在线视频只有这里精品首页| 日韩一区二区三区影片| 亚洲欧美精品综合久久99| 亚洲精品乱码久久久v下载方式| 国产综合懂色| 精品久久久久久久久久久久久| 欧美一级a爱片免费观看看| 中文资源天堂在线| 国语自产精品视频在线第100页| 看免费成人av毛片| 国产极品精品免费视频能看的| 男女下面进入的视频免费午夜| 看免费成人av毛片| 欧美日韩综合久久久久久| 寂寞人妻少妇视频99o| 91aial.com中文字幕在线观看| 久久精品夜夜夜夜夜久久蜜豆| 亚州av有码| 深爱激情五月婷婷| 插阴视频在线观看视频| 最近的中文字幕免费完整| 如何舔出高潮| 精品久久久久久久久久免费视频| 久久久欧美国产精品| 91精品一卡2卡3卡4卡| 在线观看66精品国产| 成人毛片60女人毛片免费| 2021天堂中文幕一二区在线观| 色综合亚洲欧美另类图片| 久久久久久伊人网av| 亚洲最大成人av| 天美传媒精品一区二区| 午夜福利在线观看吧| 简卡轻食公司| 啦啦啦啦在线视频资源| 精品国内亚洲2022精品成人| 亚洲国产精品合色在线| 国产精品一区二区性色av| 亚洲欧美日韩无卡精品| 国产精品一区二区三区四区久久| 乱人视频在线观看| 男人和女人高潮做爰伦理| 高清毛片免费看| 久久久欧美国产精品| 全区人妻精品视频| 久久人人精品亚洲av| 国产精品爽爽va在线观看网站| 国产精品久久视频播放| 熟女人妻精品中文字幕| 内地一区二区视频在线| 看片在线看免费视频| 亚洲国产精品合色在线| 日韩av不卡免费在线播放| 日韩欧美一区二区三区在线观看| 亚洲无线在线观看| 禁无遮挡网站| 亚洲精品456在线播放app| 天堂av国产一区二区熟女人妻| 人妻夜夜爽99麻豆av| 久久精品国产清高在天天线| 狠狠狠狠99中文字幕| 一进一出抽搐动态| 久久精品国产鲁丝片午夜精品| 久久久久免费精品人妻一区二区| 身体一侧抽搐| 精华霜和精华液先用哪个| 精品久久国产蜜桃| 国模一区二区三区四区视频| 婷婷色av中文字幕| 日韩欧美 国产精品| 精品人妻偷拍中文字幕| 神马国产精品三级电影在线观看| 草草在线视频免费看| 欧美一区二区亚洲| 春色校园在线视频观看| 欧美成人免费av一区二区三区| 青春草视频在线免费观看| .国产精品久久| 日本与韩国留学比较| av又黄又爽大尺度在线免费看 | 两个人的视频大全免费| 高清日韩中文字幕在线| 搡老妇女老女人老熟妇| 一本一本综合久久| 精品久久久久久久久久久久久| 99久久九九国产精品国产免费| 国产乱人视频| 免费不卡的大黄色大毛片视频在线观看 | 国产亚洲欧美98| 成人一区二区视频在线观看| 国产精品国产高清国产av| 国产高潮美女av| 国产精品,欧美在线| 欧美日韩乱码在线| 简卡轻食公司| 国产亚洲av片在线观看秒播厂 | 成人国产麻豆网| 久久综合国产亚洲精品| 亚洲国产精品成人久久小说 | 婷婷六月久久综合丁香| 搡老妇女老女人老熟妇| 简卡轻食公司| 又爽又黄a免费视频| 黄色一级大片看看| 国产av在哪里看| 日韩制服骚丝袜av| 亚洲七黄色美女视频| 成年免费大片在线观看| 午夜精品一区二区三区免费看| 婷婷精品国产亚洲av| 村上凉子中文字幕在线| 哪里可以看免费的av片| 久久九九热精品免费| 99热这里只有是精品50| 国产高清不卡午夜福利| 看非洲黑人一级黄片| 国产精品女同一区二区软件| 国产精品一区www在线观看| 又黄又爽又刺激的免费视频.| 成年版毛片免费区| 村上凉子中文字幕在线| 欧美一区二区亚洲| 非洲黑人性xxxx精品又粗又长| 久久久精品大字幕| 成年版毛片免费区| 18禁黄网站禁片免费观看直播| 国产精品一区二区三区四区久久| 久久韩国三级中文字幕| 国产精品久久久久久久久免| 久久鲁丝午夜福利片| 在线观看av片永久免费下载| 国产精品一及| av卡一久久| 午夜福利成人在线免费观看| 久久亚洲国产成人精品v| 国产乱人视频| 中文精品一卡2卡3卡4更新| 久久午夜亚洲精品久久| 亚洲精品乱码久久久v下载方式| 国产免费男女视频| 国内少妇人妻偷人精品xxx网站| 久久人人爽人人爽人人片va| 最近手机中文字幕大全| 99热这里只有精品一区| 国产精品乱码一区二三区的特点| 久久99蜜桃精品久久| 国产精品女同一区二区软件| 非洲黑人性xxxx精品又粗又长| 久久人妻av系列| 国产精品乱码一区二三区的特点| 欧美日韩在线观看h| 乱人视频在线观看| 日韩一区二区三区影片| 欧美变态另类bdsm刘玥| 亚洲人成网站在线播放欧美日韩| 久久久成人免费电影| 午夜福利在线观看吧| 国产精品久久电影中文字幕| 看十八女毛片水多多多| 亚洲成人中文字幕在线播放| 亚洲aⅴ乱码一区二区在线播放| 日本黄色片子视频| 特级一级黄色大片| 最好的美女福利视频网| 噜噜噜噜噜久久久久久91| 亚洲成a人片在线一区二区| 国产伦一二天堂av在线观看| 蜜臀久久99精品久久宅男| 国产一区二区三区在线臀色熟女| 亚洲欧美中文字幕日韩二区| 日韩视频在线欧美| 麻豆久久精品国产亚洲av| 国产精品国产三级国产av玫瑰| 性色avwww在线观看| 亚洲欧美日韩高清在线视频| 久久精品91蜜桃| 亚洲内射少妇av| 最近的中文字幕免费完整| 在线播放无遮挡| 久久久午夜欧美精品| 亚洲三级黄色毛片| 在线a可以看的网站| 亚洲综合色惰| 午夜福利视频1000在线观看| 国产又黄又爽又无遮挡在线| 中文字幕av成人在线电影| 欧美精品国产亚洲| 99在线视频只有这里精品首页| 精品国内亚洲2022精品成人| 日韩av不卡免费在线播放| 波多野结衣高清作品| 一本一本综合久久| 99热6这里只有精品| 啦啦啦观看免费观看视频高清| 亚洲国产日韩欧美精品在线观看| 在线观看66精品国产| 国产精品av视频在线免费观看| 亚洲五月天丁香| 观看免费一级毛片| 人体艺术视频欧美日本| 亚洲国产欧美在线一区| 成人性生交大片免费视频hd| 赤兔流量卡办理| 人人妻人人澡欧美一区二区| 黄色视频,在线免费观看| 蜜桃亚洲精品一区二区三区| 久久久久久久久久成人| 亚洲国产精品sss在线观看| 一个人观看的视频www高清免费观看| 日韩欧美国产在线观看| av.在线天堂| 国产黄a三级三级三级人| 狂野欧美白嫩少妇大欣赏| 97超碰精品成人国产| 美女大奶头视频| 免费人成视频x8x8入口观看| 日本在线视频免费播放| 成人特级黄色片久久久久久久| АⅤ资源中文在线天堂| 人妻少妇偷人精品九色| 成人永久免费在线观看视频| 午夜福利在线在线| 综合色丁香网| 亚洲不卡免费看| 身体一侧抽搐| 乱人视频在线观看| 亚洲av不卡在线观看| 少妇高潮的动态图| 亚洲av男天堂| 午夜激情欧美在线| 日日摸夜夜添夜夜添av毛片| 亚洲精品日韩av片在线观看| 最后的刺客免费高清国语| 免费人成在线观看视频色| 波多野结衣高清无吗| 黄色视频,在线免费观看| 老司机福利观看| 欧美最新免费一区二区三区| 精品午夜福利在线看| 国产精品美女特级片免费视频播放器| 国产亚洲精品av在线| 日韩中字成人| 99riav亚洲国产免费| 久久这里有精品视频免费| 精品人妻视频免费看| 丝袜美腿在线中文| 国产精品麻豆人妻色哟哟久久 | 国产中年淑女户外野战色| 内射极品少妇av片p| 亚洲精品成人久久久久久| 一本久久中文字幕| 欧美最黄视频在线播放免费| 日本与韩国留学比较| 欧美激情在线99| 日本一本二区三区精品| 22中文网久久字幕| 一区二区三区四区激情视频 | 亚洲成人精品中文字幕电影| 少妇猛男粗大的猛烈进出视频 | 可以在线观看毛片的网站| 亚洲成人中文字幕在线播放| 国产亚洲av嫩草精品影院| 看片在线看免费视频| 成人三级黄色视频| 插阴视频在线观看视频| 麻豆一二三区av精品| eeuss影院久久| 国产精品嫩草影院av在线观看| 色5月婷婷丁香| 亚洲一级一片aⅴ在线观看| 日本撒尿小便嘘嘘汇集6| 国产精品1区2区在线观看.| 亚洲熟妇中文字幕五十中出| 美女高潮的动态| 欧美最新免费一区二区三区| 九色成人免费人妻av| 女的被弄到高潮叫床怎么办| 国产精品一区二区性色av| 亚洲av第一区精品v没综合| 久久久a久久爽久久v久久| 久久精品人妻少妇| 免费人成视频x8x8入口观看| 日韩大尺度精品在线看网址| 日日啪夜夜撸| 亚洲五月天丁香| 免费人成在线观看视频色| 97人妻精品一区二区三区麻豆| 18禁裸乳无遮挡免费网站照片| 赤兔流量卡办理| 18禁黄网站禁片免费观看直播| 午夜精品在线福利| 一个人看视频在线观看www免费| 色吧在线观看| 亚洲成a人片在线一区二区| 小说图片视频综合网站| 九九热线精品视视频播放| 精品久久久久久久久久久久久| 亚洲激情五月婷婷啪啪| 国产成人a∨麻豆精品| 亚洲丝袜综合中文字幕| 少妇被粗大猛烈的视频| 日本欧美国产在线视频| 亚洲,欧美,日韩| 国产亚洲91精品色在线| 美女高潮的动态| 三级国产精品欧美在线观看| 99久久无色码亚洲精品果冻| 啦啦啦观看免费观看视频高清| 精品人妻视频免费看| 99riav亚洲国产免费| 尾随美女入室| av又黄又爽大尺度在线免费看 | 亚洲欧美日韩卡通动漫| 亚洲av成人精品一区久久| 色噜噜av男人的天堂激情| 国产美女午夜福利| 白带黄色成豆腐渣| 变态另类丝袜制服| 欧美高清性xxxxhd video| 成人三级黄色视频| 给我免费播放毛片高清在线观看| 91狼人影院| 国产黄a三级三级三级人| 亚洲精品乱码久久久v下载方式| 久久久午夜欧美精品| 毛片女人毛片| 在线播放无遮挡| 精品午夜福利在线看| 如何舔出高潮| 99在线视频只有这里精品首页| 日韩av在线大香蕉| 久久精品久久久久久久性| 欧美区成人在线视频| 97超视频在线观看视频| 能在线免费观看的黄片| 日本成人三级电影网站| 嫩草影院精品99| 又爽又黄a免费视频| 91aial.com中文字幕在线观看| 欧美色视频一区免费| 亚洲经典国产精华液单| 亚洲最大成人av| 免费人成在线观看视频色| 亚洲av男天堂| 精品日产1卡2卡| 色哟哟哟哟哟哟| 男女视频在线观看网站免费| 亚洲在久久综合| 男女边吃奶边做爰视频| 国产 一区 欧美 日韩| 3wmmmm亚洲av在线观看| 美女 人体艺术 gogo| 欧美日韩一区二区视频在线观看视频在线 | 级片在线观看| 欧美日韩在线观看h| 狂野欧美白嫩少妇大欣赏| 国产老妇女一区| 女人十人毛片免费观看3o分钟| 亚洲丝袜综合中文字幕| 黑人高潮一二区| 午夜激情福利司机影院| 亚洲第一区二区三区不卡| 免费搜索国产男女视频| 欧美激情国产日韩精品一区| 有码 亚洲区| 十八禁国产超污无遮挡网站| 少妇丰满av| 青春草国产在线视频 | 99热网站在线观看| 美女高潮的动态| 日韩精品青青久久久久久| 国产老妇伦熟女老妇高清| 波野结衣二区三区在线| 一夜夜www| 日韩在线高清观看一区二区三区| 国产成人精品婷婷| 真实男女啪啪啪动态图| 热99在线观看视频| 久久精品国产99精品国产亚洲性色| 国产乱人视频| 国产片特级美女逼逼视频| av福利片在线观看| 18禁在线播放成人免费| 麻豆一二三区av精品| 91麻豆精品激情在线观看国产| 免费黄网站久久成人精品| 自拍偷自拍亚洲精品老妇| 国产麻豆成人av免费视频| 欧美三级亚洲精品| 夜夜爽天天搞| 国产精品乱码一区二三区的特点| 亚洲最大成人av| a级一级毛片免费在线观看| 欧美精品国产亚洲| 国产乱人偷精品视频| 在线观看一区二区三区| 女人被狂操c到高潮| 免费看a级黄色片| 99热6这里只有精品| 欧美一区二区精品小视频在线| 给我免费播放毛片高清在线观看| 欧美日韩综合久久久久久| 国产成人a∨麻豆精品| 九草在线视频观看| 综合色丁香网| 国产伦理片在线播放av一区 | 麻豆成人午夜福利视频| 小说图片视频综合网站| 夜夜爽天天搞| 色综合色国产| 精品免费久久久久久久清纯| 亚洲性久久影院| 欧美性感艳星| 看十八女毛片水多多多| 亚洲精品粉嫩美女一区| av在线老鸭窝| 免费电影在线观看免费观看| 别揉我奶头 嗯啊视频| 久久国产乱子免费精品| 又爽又黄无遮挡网站| 国产三级中文精品| 日韩 亚洲 欧美在线| 日韩av不卡免费在线播放| 可以在线观看的亚洲视频| 我的老师免费观看完整版| 一本一本综合久久| 寂寞人妻少妇视频99o| 亚洲四区av| 久久精品91蜜桃| 一本久久精品| 全区人妻精品视频| 亚洲最大成人手机在线| 婷婷精品国产亚洲av| 免费搜索国产男女视频| 日韩国内少妇激情av| 狠狠狠狠99中文字幕| 特级一级黄色大片| av视频在线观看入口| 欧美+日韩+精品| 黄色视频,在线免费观看| 亚洲,欧美,日韩| av视频在线观看入口| 欧美不卡视频在线免费观看| 久久6这里有精品| 欧美人与善性xxx| 一夜夜www| 亚洲最大成人中文| 女人十人毛片免费观看3o分钟| 婷婷亚洲欧美| 好男人视频免费观看在线| .国产精品久久| 能在线免费看毛片的网站| 亚洲中文字幕日韩| 精品欧美国产一区二区三| 12—13女人毛片做爰片一| 日韩 亚洲 欧美在线| 最好的美女福利视频网| 女人被狂操c到高潮| 3wmmmm亚洲av在线观看| 成人亚洲精品av一区二区| 国产精品av视频在线免费观看| 久久精品国产鲁丝片午夜精品| 国产精品久久久久久精品电影| 国产av一区在线观看免费| 国产精品人妻久久久影院| 欧美最黄视频在线播放免费| 欧美性猛交黑人性爽| 尾随美女入室| 听说在线观看完整版免费高清| 波多野结衣巨乳人妻| 国产精品女同一区二区软件| 搞女人的毛片| 国产精品精品国产色婷婷| 国产在线男女| 日韩欧美一区二区三区在线观看| 天堂影院成人在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 国内揄拍国产精品人妻在线| 久久久国产成人免费| 看免费成人av毛片| 国产国拍精品亚洲av在线观看| 久久久国产成人精品二区| 成人亚洲欧美一区二区av| 欧美性猛交╳xxx乱大交人| 国内精品久久久久精免费| 天堂√8在线中文| 国产伦精品一区二区三区四那| 中国国产av一级| 人妻制服诱惑在线中文字幕| 美女内射精品一级片tv| 色哟哟哟哟哟哟| 麻豆乱淫一区二区| 看非洲黑人一级黄片| 老师上课跳d突然被开到最大视频| 波多野结衣高清作品| 赤兔流量卡办理|