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

    Experimental and quantum chemical studies of carboxylates as corrosion inhibitors for AM50 alloy in pH neutral NaCl solution

    2022-07-13 08:25:24JunjieYngPingliJingYouminQiuChihYuJoCrstenBlwertSvitlnLmkAnissBouliXiopengLuMikhilZhelukevichWeiLi
    Journal of Magnesium and Alloys 2022年2期

    Junjie Yng ,Pingli Jing ,Youmin Qiu ,Chih-Yu Jo ,Crsten Blwert ,Svitln Lmk ,Aniss Bouli,Xiopeng Lu,Mikhil L.Zhelukevich,e,Wei Li

    a Institute of Advanced Wear &Corrosion Resistant and Functional Materials,Jinan University,Guangzhou 510632,PR China

    b Helmholtz-Zentrum Hereon,Institute of Surface Science,Geesthacht 21502,Germany

    c School of Materials Science and Engineering,Zhengzhou University,Zhengzhou 450001,PR China

    d Shenyang National Laboratory for Materials Science,Northeastern University,3-11 Wenhua Road,Shenyang 110819,PR China

    e Faculty of Engineering,University of Kiel,Kaiserstra?e 2,Kiel 24143,Germany

    Abstract Sodium salts of mono-and di-carboxylic acids (glycolic,fumaric and benzoic acid) were studied as corrosion inhibitors for AM50 alloy in pH neutral aqueous NaCl environment.Hydrogen evolution,electrochemical and surface characterization techniques were employed to reveal their corrosion inhibition mechanism,whilst the molecular features of inhibitors were investigated by quantum chemical calculation.All inhibitors reduced anodic dissolution of AM50 and their efficiency generally increased with time and concentration from 5 mM to 100 mM.The inhibition mechanism can be described as physisorption of inhibitive molecules on the surface of the intrinsic oxide layer followed by chemisorption with Mg2+ and Al3+,and the difference in inhibition action among these inhibitors was explained on the molecular scale.

    Keywords: Mg-Al magnesium alloy;Neutral NaCl solution;Corrosion inhibitors;Electrochemistry;Density functional theory.

    1.Introduction

    Magnesium-aluminum (Mg-Al) alloys,containing 2-10 wt.% Al and minor manganese (Mn),are the most commercially used Mg alloys,especially in automotive and aerospace industries where a weight reduction can significantly minimize the fuel consumption and improve the performance of the vehicles and aircrafts [1].Apart from the enhancement in mechanical properties (e.g.tensile strength and hardness),castability and weldability,the corrosion resistance of Mg is significantly improved by adding proper amount of Al and Mn[2].Explanation for such corrosion performance of AM series alloys can be ascribed to two factors.Firstly,Mn reduces the harmful effect imposed by transition metallic impurities(e.g.iron,copper and nickel) through forming intermetallic compound that settles to the melt bottom or encapsulating impurities that remains within the metal [3].Secondly,the dissolution of Al-richα-Mg matrix favors Al enrichment on the metallic surface which improves corrosion resistance of the primary passivation layer,whereasβ-phase (Mg-Al) network precipitated along the grain boundary can act as a barrier against the progression of corrosion attack [4-7].In spite of these advantages,anticorrosion actions are still indispensable for Mg-Al alloys for industrial application.

    Corrosion inhibitors are either organic or inorganic compounds that can effectively reduce the degradation rate or even retard the corrosion process of metallic-based materials at a relatively low concentration [8].Compared with other corrosion protection techniques,reducing the degradation rate of metallic materials by using corrosion inhibitors is a more easy-going and promising strategy in practice.Corrosion inhibitors can be classified into anodic,cathodic and mixed type depending on their impact is mainly in retarding the anodic,cathodic of corrosion process or both of them [8].In principle,anodic inhibitors that hinder the anodic dissolution of metals are expected to be the best candidates for effective corrosion inhibition,since they can straightforwardly passivate,precipitate or adsorb on the surface by forming protective layer(s) on the surface.The inhibition efficiency of precipitating inorganic anodic inhibitors is closely related with the solubility product constant (Ksp) of compounds formed between the inhibitive anions and Mg2+,which determines the stability of the corrosion products layer.For instance,the superior corrosion inhibition performance of phosphate can be ascribed to the lowKspof Mg3(PO4)2(1.0×10-25,where[Mg2+]=4.9×10-6M),which is about two orders lower than that of Mg(OH)2(1.8×10-11,where[Mg2+]=1.6×10-4M)[9,10].Besides,dissolvable carbonates (6.8×10-6) [10] and fluoride (5.2×10-11) [11] were also reported as corrosion inhibitors for AZ31 Mg alloy.However,their inhibition efficiencies are not high enough especially at relative low concentrations,since the stability of their Mg containing compounds is close to that of Mg(OH)2,which can hardly change the conditions of corrosion products layer.On another hand,organic inhibitors prevent corrosion mainly through their adsorption on the metal or hydroxide surface via electrostatic interaction (physisorption) or forming adsorptive products with metallic cations released during corrosion (chemisorption).However,the out most molecular orbital (3 s) of Mg is fully occupied by a pair of spinning electrons,which implies that the adsorption of organic inhibitors on Mg surface through donating or accepting their electrons to a vacant orbital is more difficult compared to many other metallic materials.Albeit,a few exceptions e.g.2-hydroxy-4-methoxy-acetophenone(paeonol)[12],tetraphenylporphyrin(TPP) [13] and 8-hydroxyquinoline (8-HQ) [14] can complex with Mg2+via multi substitutes forming insoluble precipitates within the defective areas of the porous Mg(OH)2film.

    As a logical and promising alternative strategy in corrosion inhibition of Mg alloys,cathodic inhibitors decrease the H2O and O2reduction rate (cathodic reaction:2H2O+O2+4e-→4OH-),which in turn limit the anodic dissolution of Mg [15].The cathodic sites within Mg alloys are those electrochemically nobler components including intermetallics and impurities.Anionic surfactants (e.g.sodium dodecyl sulfate [16],sodium salt of N-lauroyl-Nmethyltaurine and dodecylbenzensulphonic acid[17])are well studied inhibitors,since they have negatively charged hydrophilic substitutes selectively adsorb on the surface of micro cathodes with leaving their long hydrophobic tails suspending towards the solution.Thus,it is understandable that a satisfactory corrosion inhibition efficiency of anionic surfactants can only be achieved at relatively high concentration when the surface of cathodic sites is fully covered by inhibitive molecules.For instance,selective adsorption of sodium dodecyl sulfate (SDS) on AlMn intermetallic of AZ91 is confirmed by Lu et al.employing electron probe micro-analyzer(EPMA) [16].They suggested that adsorption of SDS was realized via interaction between the unshared paired electrons of S-O or S=O bonds and the vacant orbital of Mn,and the inhibition efficiencyIEpreached 88.8% when 0.05 M of the inhibitor is added in 3.5 wt.% NaCl solution.Except for the intermetallics,metallic impurities (e.g.nickel (Ni),copper (Cu) and iron (Fe)) introduced during the alloy preparation process may also serve as micro cathodes and promote the corrosion of Mg alloys [18,19].A reasonable explanation about the corrosion promotion mechanism of trace Fe particles was proposed by H?che et al.[20].They claimed that Fe inclusions undergo “self-corrosion” after their loss of direct contact with bulk Mg,and the resultant Fe2+/3+cations or Fe oxide(s) would re-plated back on the Mg surface,which promote the corrosion of the Mg substrate.In view of the role of Fe impurities during Mg corrosion,Lamaka et al.elucidated that Fe2+/3+complexing agents (e.g.cyanide,salicylate,oxalate,methylsalicylate and thiocyanate) can also serve as inhibitors for Mg alloys.The inhibition efficiency of the latter complexing agents was in line with their capability of forming complexes with Fe3+[21].Unlike anodic and cathodic inhibitors,mixed type inhibitors (e.g.2,5-pyridine dicarboxylate and fumarate [22,23]) were characterized by their adsorption or formation of coordination precipitation on the surface of MgO/Mg(OH)2via carboxyl groups,which enhances the compactness of the corrosion products layer and reduces both anodic and cathodic reaction rates.

    Among the known inhibitors,carboxylates are widely accepted as effective inhibitors for Mg alloys and their inhibition mechanism was generally ascribed to the interaction between the negatively charged carboxyl groups and positively charged metallic ions on the surface [24].Nevertheless,the inhibitions efficiency of carboxylates significantly varies with respect to the composition of alloy as well as the type and position of functional groups within the inhibitor molecule[25-27].As reported by Lamaka et al.[27],the inhibition efficiency of 2,5-pyridine dicarboxylate (2,5-PDCA) on pure Mg,regardless of the content of Fe,was more than 80%,which was almost twice that of ZE41.Moreover,the inhibition efficiency for AM50 increased from minus 97% (for salicylic acid) to 29%,41% and 66% when -CH3substituted at para-(5-methylsalicylic acid),meta-(4-methylsalicylic acid)and ortho-(3-methylsalicylic acid) positions relative to -OH of salicylic acid.In light of this fact,it is of great importance to start the study from inhibitors with simple molecular structure,which can provide (1) the most basic understanding of the influence of functional groups and their position on the corrosion inhibition performance,and (2) fundamental information for machine learning approaches to design new inhibitors or predict the efficiency of untested compounds[28,29].In the present work,three characteristic compounds containing mono-or bi-carboxyl group were used as corrosion inhibitors for Mg-Al-Mn (AM50) alloy.The influence of inhibitor concentration and presence of secondary phases on the corrosion inhibition actions of each inhibitor werestudied through hydrogen evolution,electrochemical measurement and subsequent surface analysis.As follows,quantum chemical calculation based upon density functional theory(DFT) was performed to disclose the relationship between the inhibitive actions of the inhibitors and their molecular configuration which contributes to the setup of corresponding inhibition mechanism.

    Fig.1.Molecular structure of (a) glycolic acid,(b) fumaric acid and (c)benzoic acid.

    2.Experimental

    2.1.Materials and reagents

    Cast AM50 alloy was used as substrate material in the present study.Spark optical emission spectroscopy (Spark OES,Spark analyzer M9,Spectro Ametek,Germany) was employed to determine the elemental composition of the AM50 alloy,and the results are shown in Table 1.For hydrogen evolution test,the AM50 ingot was machined into chips with a larger surface area (area-to-mass ratio of 56.86±16 cm2/g) using a turning lathe equipped with a diamond cutting tool.Details about the preparation process can be found in previous work [24].Rectangular coupons with dimension of 15×15×4 mm3sliced from the ingot were used for electrochemical corrosion testing and surface characterization.Prior to the tests and characterizations,all specimens were ground using silicon carbide (SiC) papers successively up to 1200 grit.For optical microscopy observation,the specimens were further polished with a mixture of alumina (Al2O3) slurry and 1 μm non-aqueous diamond suspension.Eventually,the specimens were rinsed using deionized water and ethanol,and then dried in a pressured air flow at room temperature.

    Table 1Elemental composition of AM50 alloy (in wt.%).

    All chemicals used in present study were in analytical reagent (AR) grade and no further purification was conducted before application.Three corrosion inhibitors(Fig.1),including glycolic acid (HOCH2COOH,CAS:79-14-1,Sigma-Aldrich Chemie GmbH,Germany),fumaric acid (HOOCCH=CHCOOH,CAS:110-17-8,Alfa Aesar,Germany) and benzoic acid (C6H5COOH,CAS:65-85-0,Sigma-Aldrich Chemie GmbH,Germany),were dissolved in 0.5 wt.% NaCl solution to obtain different concentrations (i.e.5 mM,10 mM,50 mM and 100 mM).3 M NaOH was used to adjust the pH of the inhibitor-added solutions until the neutral range (pH:6.8-7.2).The pH value of the prepared solutions was monitored by a pH meter (Metrohm 691) during preparation.

    2.2.Hydrogen evolution

    Hydrogen evolution was recorded during immersion of 0.5 g AM50 chips (with specific surface area to mass ratio of 56.86 cm2/g) in 0.5 wt.% NaCl solution with and without addition of inhibitors using eudiometers (NS45/27,Neubert-Glas,Germany).A series of different concentrations of inhibitor solutions were prepared to investigate the influence of the concentration on inhibition efficiency.Magnetic stirrers were used to avoid aggregation of chips on the bottom of the bottle,and the speed was 200 rpm.All measurements were conducted for 24 h at room temperature,and the air resided in the eudiometer (less than 10 cm3) was not purged before running the experiments.The inhibition efficiency (IEH) derived from hydrogen evolution results was calculated as follows:

    2.3.Electrochemical measurements

    Interface 1010B potentiostat (Gamry,USA) was employed for potentiodynamic polarization and electrochemical impedance spectroscopy (EIS) measurements.A conventional three-electrode cell was used.In this cell,bulk AM50 specimens with 1 cm2exposed area served as the working electrode,a platinum (Pt) wire as the counter electrode and a saturated Ag/AgCl electrode as a reference electrode.All the electrochemical measurements were carried out in 0.5 wt.%NaCl solution with and without 50 mM each inhibitor.Each electrochemical measurement was repeated at least three times and all measurements were performed at room temperature.

    Prior to the polarization measurement,the working electrode (AM50 alloy) was immersed in the solution for 10 min to obtain a relatively stable open circuit potential (OCP).After that,the potential was scanned from OCP to -300 mV and 300 mV (vs.Ag/AgCl) respectively with a scan rate of 0.5 mV/s.The EIS tests were performed in the frequency range from 105Hz to 10-2Hz with applying an amplitude of 10 mV rms sinusoidal perturbation with respect to OCP.Commercial software Zview was employed for fittin the experimental impedance spectra.The inhibition efficiency derived from the EIS results can be calculated through the equation as follow:

    Fig.2.Metallographic microstructure (a),phase composition (b) and SEM microstructure (c and d) of as-cast AM50.

    2.4.AFM/SKPFM analysis

    The Volta potential measurement of AM50 alloy was carried out using atomic force microscopy(AFM,JPK NanoWizard) equipped with scanning Kelvin probe force microscopy(SKPFM) in AC mode.Cr/Pt coated silicon (Si) probe was employed for AFM/SKPFM scanning,and the resonance frequency and force constant of the probe are 75 kHz and 3 N/m,respectively.Polished samples were immediately used for SKPFM measurement to avoid oxidation.Measurement was conducted in open air at room temperature.A conductive paste was placed between the sample and the AFM stage to achieve electrical connection.The lift height of the probe was set as 100 nm above the sample surface.Scanning was performed on an area of 40×40 μm with a resolution of 512×512 pixels.The data was analyzed using JPKSPM Data Processing software.

    2.5.Morphology and composition characterization

    An optical microscope (OM,DM2500M,Leica,Germany)was used to observe the phases on the well-polished sample surface,whereas a scanning electron microscope (SEM,Tescan Vega3) equipped with an energy dispersive X-ray spectrometer (EDS) was used for further analyzing the surface morphology and elemental distribution before and after corrosion.XRD measurements were carried out using a D8 Advance AXS diffractometer (Bruker,Germany).X-ray diffraction patterns (2θ) were recorded from 10° to 80° using Cu Kα(λ=1.54,056 °A)radiation at room temperature.The tube voltage and tube current were 40 kV and 40 mA respectively,whereas the scanning speed was 0.1°/s.

    2.6.Quantum chemical calculation

    Geometrical optimization of deprotonated glycolic acid,fumaric acid and benzoic acid was conducted using Gaussian 16 program package at the B3LYP/6-31+G??level of density functional theory (DFT).Effects of water solvation are included by using a self-consistent reaction field (SCRF) with polarizable continuum model (PCM).

    3.Results

    3.1.Microstructure,phase composition and AFM/SKPFM analysis

    The investigation of metallographic (Fig.2a) and SEM(Fig.2c and d) microstructure and phase composition(Fig.2b) of AM50 was conducted on the as-polished AM50 alloy.As shown in the optical microscopy (OM) image,the investigated AM50 alloy exhibits a typical metallographic microstructure featured by an average grain size of 1.44±0.41 μm determined according to ASTM E112-88.In addition to the primaryα-Mg phase,two different types of intermetallic precipitates (β-Mg17Al12and Al8Mn5) are identified in the XRD patterns,which can be clearly observed by different contrast in backscattered electron mode since the electron scattering ability of these phases is different.According to the EDS analysis (Fig.2c inset),the darker phases contain Al and Mg,which can be correlated withβ-Mg17Al12peaks in the XRD pattern.Meanwhile,the light color phases contain Al and Mn,which are preferentially referred to Al8Mn5phase considering the characteristic XRD pattern.In addition,the large partially divorced eutecticβ-Mg17Al12exhibits an island-like morphology containing lamellar eutectic Mg sphere residues,whilst the Al8Mn5precipitate is commonly in polygonal shape.

    Fig.3.Topography (a) and VPD maps (b) acquired upon polished AM50 substrate employing AFM/SKPFM.

    Fig.3 illustrates topography and Volta potential difference(VPD) measured on a specific area (40×40 μm) that contains bothβ-Mg17Al12and Al8Mn5phases.As shown in the topography map (Fig.3a),obvious contrast can be observed betweenβ-Mg17Al12and Al8Mn5phases implying significant height difference between the two phases.The peak heights of Al8Mn5are 209.5 nm and 305.6 nm respectively,whereas the maximum height value ofβ-Mg17Al12is only 61.5 nm.Such a height difference can be mainly resulted from the hardness differences of these phases.β-Mg17Al12is much softer compared to Al8Mn5intermetallic but still harder than theα-Mg matrix.Hence,the hard Al8Mn5intermetallic is less removed by grinding during the preparation process leaving such a comet-like trajectory along the AFM/SKPFM scanning direction [30].Following topography measurement,the VPD measurement is performed over the same area.As shown in Fig.3b,the VPD illustrates a trend similar to the topography results.The VPD values of P1 and P2 (Al8Mn5intermetallic)are 165.6 mV and 131.2 mV respectively,whilst theβ-phase(P3) shows an average of 52.5 mV in the VPD map.Also,it is important to notice that the largeβ-Mg17Al12particles are not uniform in VPD containing regions featuring higher values than the average.This is attributed to the casting process that solidifies close to equilibrium and generates largerβ-Mg17Al12grains [31].As a result,Al partially dissolves into theβ-Mg17Al12intermetallic.Based on the VPD results,microgalvanic corrosion is expected,especially between Al-Mn intermetallic andα-Mg,which is the dominant corrosion mechanism for AM50 alloy [32].

    3.2.Hydrogen evolution analysis

    Fig.4 shows the variation of H2volume as function of immersion time of 0.5 g AM50 chips(the overall exposed area is estimated to be 28.43 cm2)in 0.5 wt.%NaCl solution without and with addition of different concentrations (5 mM,10 mM,50 mM and 100 mM) of the three inhibitors.Hydrogen evolution was used as a basic technique for screening inhibitors and metallic chips with increased the surface area were employed with the aim of obtaining results in a relatively short period,as shown in Fig.4a [33,34].Fig.4b-d show the H2evolution curves gained from immersing 0.5 g AM50 chips in pH neutral NaCl solution (0.5 wt.%) containing various concentrations of inhibitors.Similar H2release trend can be seen for glycolate (Fig.4b) and benzoate (Fig.4d),where the hydrogen evolution rate (HER) generally decreases with increasing concentration of the inhibitor except for 100 mM glycolate addition,which promotes H2evolution.This unexpected H2volume increase witnessed upon the 100 mM glycolate addition specimen may be attributed to the further deprotonation of hydroxyl induced by the pH increase due to alloy corrosion.In contrast,the corrosion of AM50 chips is instantly inhibited by the addition of fumarate into the solution,and increasing the concentration from 5 mM to 100 mM does not affect the amount of evolved H2to a large extent.The calculated corrosion inhibition efficiencies at 4 h and 15 h immersion with the different inhibitors’ addition is given in Fig.5.As expected,the efficiency values obtained at 15 h are higher than those at 4 h for all inhibitors,and fumarate reveals the highest inhibition efficiencies among three tested inhibitors.It has been shown by multiple examples before that the absolute value of inhibiting efficiency is often dependent on immersion time [35].The inhibition efficiency of glycolate,fumarate and benzoate after stabilization is 57.5%(50 mM),81.5% (100 mM) and 59% (100 mM),respectively.

    3.3.Potentiodynamic polarization

    DC potentiodynamic polarization was conducted in two separate sweeps from OCP to -300 mV and from OCP to 300 mV (vs.Ag/AgCl) for the bulk AM50 specimens immersed in NaCl solutions (0.5 wt.%) in absence and presence of 50 mM inhibitors (Fig.6).Cathodic Tafel extrapolation,as described in [36],is adopted to estimate the corrosion current density (icorr) from the polarization curves,and the critical parameters derived from the polarization curves are listed in Table 2.Since only one representative curve was selected from three or more repetitive measurements,slight difference can be identified between the plots in Fig.6 and averaged values (with errors) in Table 2.As shown in Fig.6a,theEcorrof reference group measured in blank NaCl solution without inhibitor addition shifts gradually from -1.54 V after 10 min immersion to -1.49 V after 24 h,whilst the corrosion current density stabilizes around 10 μA/cm2after 3 h immersion.Moreover,both anodic and cathodic slopes of the polarization curves are slightly suppressed during immersion,which suggest that the intrinsic hydroxide layer on the surface of Mg alloy would impede both anodic Mg dissolution and cathodic hydrogen evolution during the corrosion of AM50 alloy.With addition of carboxylates in the blank NaCl solution,obvious corrosion inhibition mainly occurs on the anodic branches.As for glycolate and benzoate,the initial(10 min)polarization curves shift towards nobler potential and lower current density direction comparing to that of the control specimen,which indicate suppression of corrosion by both inhibitors.More importantly,pseudo-passive behaviors,characterized by moderate increase of current density with potential shift,appear on the anodic branches of glycolate and benzoate specimens after 6 h and 3 h,respectively,which indicate the barrier property of hydr(oxide) is strengthened by inhibitors.However,breakdown of the inhibitor incorporated layers occurs at specific potentials (Ebd),which suggests such an enhancement is limited.Up to 24 h,benzoate registers lowericorr(2.5±3.0 μA/cm2) than that (4.6±3.5 μA/cm2)of glycolate,which suggests better corrosion inhibition ability.Similar potentiodynamic behaviors are also observed upon the specimen with 50 mM fumarate addition,whose polarization curves shift towards positive potential and lower current density with achieving anicorrof 2.5±3.0 μA/cm2by 24 h.Besides,more significant passivation behavior can be define upon the anodic branches of fumarate specimen since 3 h,and even theEbdis not shown since 6 h immersion.Such a remarkable suppression of anodic behavior indicates improved barrier property of the fumarate modified hydr(oxide) layer,which would result in better corrosion performance.

    Table 2Critical polarization parameters for AM50 immersed in NaCl solution (0.5 wt.%) without and with addition of 50 mM inhibitors.

    Fig.4.Released hydrogen volume in function of time of AM50 in NaCl solution (pH=6.8-7.2) without and with addition of inhibitors.(a) Comparison of hydrogen volume between bulk and chipped AM50,(b) glycolate,(c) fumarate and (d) benzoate.

    Fig.5.Corrosion inhibition efficiency (%) of three tested inhibitors (glycolate,fumarate and benzoate) in the concentration range from 5 mM to 100 mM after 4 h and 15 h immersion.

    Fig.6.Potentiodynamic polarization curves of AM50 in NaCl solution (0.5 wt.%) without and with addition of 50 mM inhibitors.(a) Blank NaCl solution,(b) glycolate,(c) fumarate and (d) benzoate.

    3.4.Electrochemical impedance spectroscopy (EIS)

    Electrochemical impedance spectroscopy (EIS) was employed to study the corrosion process and interface conditions in NaCl solution without and with addition of 50 mM inhibitors.Nyquist and Bode plots recorded after consecutive immersion periods (10 min,1 h,3 h,6 h,12 h and 24 h)are shown in Fig.7 at the same scale.It can be observed from the Nyquist that all inhibitor protected specimens reveal larger semi-circles than the control one,which indicate superior corrosion performance.In the blank NaCl solution(Fig.7a and b),a double semi-circle pattern is shown in the Nyquist spectrum shortly after 10 min immersion,which then merge into a larger one after 1 h immersion.In general,the low-frequency (LF:10°-10-2Hz) and high-frequency (HF:103-10° Hz) loops can be associated to the electrochemical processes occurring at the metal/electrolyte interface and hydr(oxide)layer,respectively.Hence,the combination of two loops might be attributed to the loss of the barrier property of the hydr(oxide) layer.Based on this thought,the contribution of inhibitors on the corrosion performance of AM50 alloy is mainly resulted from their adsorptive modifi cation of the hydr(oxide) layer.Similar evolution trend of EIS curves can be found on the glycolate and fumarate specimens(Fig.7c-f).The diameters of the double semi-circles in the Nyquist plots of both specimens demonstrate a sharp increase by 1 h immersion with achieving their highest impedance of 5 kohm·cm2and 4.9 kohm·cm2at 10-2Hz,respectively.Since then,the impedance decreases but fluctuate indicating the glycolate and fumarate modified hydr(oxide) cannot completely retard Cl-from penetrating this layer and the substrate is still vulnerable to corrosion attack.In contrast,50 mM benzoate added specimen demonstrates more effective corrosion inhibition action as manifested by the consecutive shift of the middle frequency (MF:103-10° Hz) phase towards lower frequency and growth of overall impedance from 2.8 kohm·cm2(10 min) to 6.9 kohm·cm2(24 h)(Fig.7h).These results suggest that the barrier property is enhanced with time,which may result from the increase of compactness of the hydr(oxide) layer due to the adsorption of inhibitive molecules.

    Fig.7.EIS curves of AM50 alloy immersed in blank NaCl solution without and with 50 mM inhibitors within 24 h.(a)(b) Blank NaCl solution,(c)(d)Glycolate,(e)(f) Fumarate and (g)(h) Benzoate.

    Fig.8.Equivalent circuit employed to fi the experimental EIS data in blank NaCl solution without and with addition of 50 mM inhibitors.

    To quantify the physical-chemical processes involved within different corrosion cells,a typical equivalent circuit(Fig.8) comprising of two time-constants is used to fi the experimental EIS results.In this circuit,Rs stands for the resistance of NaCl solution with and without 50 mM inhibitors,CPEfandRfare associated with the capacitance and resistance of the protective layer on the top of AM50,andCPEdlandRctrepresent the double layer capacitance and charge transfer resistance occurred at the metal/electrolyte interface.Considering the imperfect surface or interface condition caused by the non-uniform formation of corrosion products,distribution of charge and adsorption of inhibitive molecules [26],constant phase element,rather than pure capacitor,was used for the simulation.In addition,effective capacitance (Cf) was calculated depending on the following equation:

    whereωmaxis the frequency at which the imaginary impedance achieves the maximum value for the respective constant [37].The evolution of film resistance (Rf),film effective capacitance (Cf),charge transfer resistance (Rct) and corrosion inhibition efficiency (IEE) in function of time is depicted in Fig.9.As shown in Fig.9a,the blank and inhibitor protected specimens register quite similarRfvalues (approx.1.3 kohm·cm2) at the beginning of immersion,which implies that the electrochemical reactions and property of hydroxide layer are not significantly affected by the inhibitors.But after 1 h,theRfof the control one remains almost unchanged,whereas the inhibitor incorporated groups reveal significant increase (to 2.5 kohm·cm2(glycolate),2.8 kohm·cm2(fumarate) and 3.7 kohm·cm2(benzoate) respectively followed by slight fluctuations In the meantime,all inhibitor protected specimens demonstrate constant increase ofCfwith time achieving 16.0±2.2 μF/cm2,23.2±2.6 μF/cm2and 15.0±0.7 μF/cm2by 24 h (Fig.9b).These consistent trend inCfmay result from hydrolytically degradation of the alloy and oxide [38] or growth and compactness of the hydr(oxide)layer.Thus,the faster adsorption of fumarate on the corroding surface would result in thinner and/or more compact hybrid hydr(oxide) layer.Furthermore,influence by the resistance of the hydr(oxide)layers,a similar variation tendency is found upon theRct(Fig.9c).TheRctof the control specimen meets a consecutive reduction from 0.8 kohm·cm2(10 min)to 0.2 kohm·cm2(24 h).Simultaneously,glycolate and fumarate,having lower initialRctvalues (0.3 kohm·cm2and 0.6 kohm·cm2),exhibit a fast increase to 1.4 kohm·cm2and 1.3 kohm·cm2at 1 h,and then stabilize until 24 h.Moreover,the addition of 50 mM benzoate results in a higherRct(1.0 kohm·cm2) since the beginning,which is followed by a moderate increase to 3.4 kohm·cm2by the end of immersion.Fig.9d illustrates the corrosion inhibition efficiency(IEE)calculated from comparing the overall impedance difference between the specimens with and without inhibitorto that of the former partAs shown,the efficiency values surge significantly from-43%(glycolate),6%(fumarate)and 28% (benzoate) at 10 min to 51%,60% and 72% (1 h),and then remain slight increase tendency 66%,65% and 88%till 24 h.The initial negative value in efficiency for glycolate indicates corrosion acceleration,which may be induced by the reduced rate of hydroxide thickness and/or formation of soluble complex with Mg2+[22,23].

    3.5.Morphology of corroded surface

    In order to investigate the inhibition mechanism,the surface and cross-section morphology of corroded specimens obtained after 24 h immersion in NaCl electrolyte without and with addition of 50 mM inhibitors are shown in Fig.10 As can be seen in Fig.10a,AM50 alloy illustrates a conventional corrosion morphology,which is characterized by a relatively thick layer and inhomogeneous skeletal surface morphology.Someβ-Mg17Al12phases can still be vaguely observed along the grain boundaries,whilst the smaller AlMn intermetallics are covered by corrosion products and are hard to be seen at current resolution.Taking a closer look at the surface,one can realize the corroded surface is not compact.Instead,islands of large loose areas decorated with scattered domes can be distinguished on the surface,which might be a result of dehydration of Mg(OH)2in atmosphere and/or H2evolution during immersion [39].Much more information about the mentioned features can be obtained from the cross-section view in Fig.10b and c.It is obvious that theα-Mg phase around the AlMn intermetallics(Fig.10b) is more preferentially corroded as a result of microgalvanic coupling caused by the high potential difference.As a result,the corrosion product of Mg(OH)2predominately forms in these areas even completely burying the secondary phase particles.However,as clearly indicated in the left inset of Fig.10b,the layer of Mg(OH)2precipitate reveals a porous and laminated microstructure that can inevitably be penetrated by the electrolyte containing Cl-species.Hence,it can be expected that severer localized damage is about to occur (right inset of Fig.10b) with increasing immersion period.It is interesting to notice that Mg17Al12phases may drop out of Mg substrate being surrounded by corrosion products on the surface.In contrast,corrosion products are rarely observable on the top of Mg17Al12phases,which indicates the galvanic coupling betweenβ-Mg17Al12andα-Mg may be not fully activated in the relatively diluted NaCl solution.

    Fig.9.Evolution of film resistance Rf (a),film capacitance (b),charge transfer resistance Rct (c) and corrosion inhibition efficiency IEI (d) during immersion in 0.5 wt.% NaCl solution without and with addition of 50 mM inhibitors.

    The corrosion inhibition of AM50 alloy due to addition of 50 mM inhibitors in NaCl solution are also manifested by a corroded surface morphology (Fig.10d,g and j).Primarily,paralleled scratches induced during the polishing process can still be identified on the surfaces,which indicates the anodic dissolution of Mg is retarded by the inhibitors and less Mg(OH)2is formed.Cracks appear on the dark areas,which highlights the places where severe galvanic corrosion occurs.Both intermetallicsβ-Mg17Al12and AlMn are observable on the corroded surface in the case of fumarate and benzoate addition,but they are hard to see on the glycolate protected specimen.These surface observations are further confirmed by their cross-section views.The local areas containing AlMn inclusion illustrate absence of galvanic corrosion attack for fumarate and benzoate specimens (Fig.10e and h),whereas glycolate still reveals slightα-Mg dissolution caused by its galvanic coupling with AlMn (Fig.10k).Moreover,it can be expected that the corrosion behaviors ofβ-Mg17Al12phase is not much influence by the inhibitors (Fig.10f,i and l),since almost no corrosion can be distinguished in theβ-Mg17Al12included regions and this situation is the same as in the control group.Considering the corroded surface morphology analysis,the corrosion inhibition functions of present three inhibitors can be ascribed to their adsorption ability on the oxide surface and specifically on the AlMn intermetallics also impeding their galvanic coupling withα-Mg.

    4.Discussion

    4.1.Effect of carboxylates on the corrosion of AM50 in NaCl solution

    Combining the metallographic observation (Fig.2),phase VPD determination (Fig.3) and corroded surface analysis(Fig.10a-c),it is not hard to figure out that the microgalvanic coupling betweenβ-Mg17Al12phase/Al8Mn5intermetallic andα-Mg is the dominant factor that triggers the corrosion of AM alloys in pH neutral Cl-containing solution.Within these galvanic couples,theα-Mg is continuously corroding during the aqueous exposure with leaving cathodicβ-Mg17Al12phase and Al8Mn5intermetallic unattacked.As corrosion proceeds,Mg(OH)2precipitates on the surface with a less compact and inhomogeneous surface microstructure.Both H2release and pH increase caused by the cathodic reaction (2H2O+2e-→H2+2OH-)may contribute to these morphology features.Al enriches in the skeletal network areas left after the dissolution ofα-Mg and in the primary corrosion product layer,which acts not only as boundaries that impedes the propagation of corrosion among the surface but also improves the corrosion resistance of the corrosion product(s) layer [40].

    Fig.10.Corroded surface (a,d,g and j) and cross-section (b,c,e,f,h,I,k and l) observation of AM50 after immersion in NaCl solution without (a,b and c) and with addition of 50 mM of glycolate (d,e and f),fumarate (g,h and i) and benzoate (j,k and l).

    Despite the fact that all three inhibitors can effectively strengthen the resistive capability of AM50 against corrosion attack in NaCl solution,their adsorption behaviors vary significantly .As shown in Fig.4b-d,the initial (within 0 h-5 h) H2evolution tendency of glycolate and benzoate reveal slight difference comparing to that of the control group,whereas fumarate reveals significant suppression even at the lowest concentration.Such a difference can be ascribed to the presence of double carboxyl moieties within fumarate molecule that benefit its physisorption at the metal (negative charge)/electrolyte (positive charged species) interface.Further,the corrosion inhibition efficiency generally increases consistently with increasing concentration and immersion duration (Figs.4 and 5).Given an exception that the inhibition efficiency of glycolate dramatically decreases from 50.8% to 1.2% along with the concentration increase from 50 mM to 100 mM at 4 h may result from a further deprotonation of-OH group,which is stimulated by the increasing pH as a result of the cathodic reaction (2H2O+2e-→H2+2OHand/or 2H2O+O2+4e-→4OH-).Another plausible explanation might be formation of soluble Mg-glycolate complex at higher concentration of glycolate that results in corrosion acceleration.However,the corrosion rate gradually reduces when chemisorption of glycolate on the hydroxide occurs,which modified the microstructure of the hydroxide layer and improves its compactness [41].

    The difference in corrosion inhibition behavior of these studied carboxylates can be further understood from the potentiodynamic polarization results in Fig.6.All three investigated compounds behave as anodic inhibitors judging from the consecutive shift ofEcorr,which implies that the anodic dissolution ofα-Mg is restricted.Besides,both fumarate and benzoate incorporated specimens reveal a “pseudo-passive”behavior in the anodic branch of polarization curve shortly after 10 min and 24 h exposure respectively to the inhibitor containing NaCl solution,which can hardly be define upon the glycolate one.These difference in anodic branch of polarization curves indicate that the barrier property of hydroxide layer is gradually reinforced by fumarate and benzoate whilst such a strengthening effect on hydroxide is much weak for glycolate.

    Fig.11.Mullikin charge (a,d and g) and the frontier molecule orbital density distributions (HOMO (b,e and h) and LUMO (c,f and i)) of deprotonated glycolic,fumaric and benzoic acid.

    This hydroxide enhancement effect can be also elucidated by the EIS results,in which the diameters of HF semi-circle in Nyquist spectra of benzoate group are not only much larger than those of glycolate and fumarate but also keep constant increase with immersion time.However,discrepancy in the inhibition efficiency evaluation can be distinguished between the results obtained from electrochemical measurements and hydrogen evolution.In details,the efficiency of glycolate,fumarate and benzoate (all in 50 mM,@15 h) from hydrogen evolution is 57.7%,79.1% and 50.6% respectively,whereas the efficiency values interpreted from EIS are 65.4%,65.0%and 80.1% (@12 h).These inconsistency in efficiency estimation can be mainly ascribed to the state of material used in measurement.For hydrogen evolution test,Mg chips characterized with much higher surface area (56.86 cm2/g) are immersed in the testing solutions and kept stirring,which benefit the physisorption of organic molecules on the surface and reduction of the degradation rate of the alloy.However,due to the constant mechanical stirring evolved during hydrogen evolution,formation of hydroxide on the surface of Mg chips and subsequent chemisorption of inhibitor on the hydr(oxide) may be impaired,which in turn induces relatively lower corrosion inhibition rate than that obtained from EIS (i.e.glycolate and benzoate).From this point of view,inhibition efficiency calculated from EIS is much closer to the real atmospheric corrosion condition since only a sinusoidal 10 mV amplitude potential perturbation is exerted during the measurement,which is believed to cause much weak impact on the adsorption behaviors of inhibitors.As a result,evaluating the performance of inhibitor employing electrochemical approaches conducted on specific surface area is more reliable especially when the investigated parts serve under free corrosion conditions.

    4.2.Corrosion inhibition mechanism of studied carboxylate inhibitors

    To further explain the molecular difference on the corrosion inhibition performance of carboxylates,molecular analysis was conducted employing quantum chemical calculation method.Fig.11 shows the optimized molecular structure and distribution of HOMO and LUMO of glycolate,fumarate and benzoate,respectively.It is worth to mention that all inhibitors are calculated in their deprotonated state and the solvent influence (i.e.,water) is considered.It is not surprising that the most negative charge can be found upon theOatoms within the carboxyl and hydroxyl moieties of the inhibitors,which denotes the binding sites that most likely adsorb on the metal surface through electrostatic attraction.TheEHOMOandELUMOelectron clouds distribution within each inhibitor define the occurrence probability of electrons,which determines if electrons may be donated or accepted.As shown,the electron donating sites within present inhibitive molecules are conformably define around their carboxyl groups,whereas their electron accepting locations can be found in the vicinity of hydroxyl (glycolate),C=C bond (fumarate) and benzene ring (benzoate).

    Fig.12.Corrosion inhibition mechanisms of carboxylates on AM50 in NaCl solution.(a) Adsorption on passivation layer and (b) Formation of complex during hydroxide precipitation.

    Considering the experimental and DFT study on corrosion inhibition behaviors and molecular features of glycolate,fumarate and benzoate,the feasible inhibition mechanism can be described as shown in Fig.12.As supported by the electrochemical results(Figs.6 and 7),the firs scenario of corrosion inhibition action of these inhibitors can be descripted as their simultaneous adsorption on intrinsic passivation film with a few nanometers thickness,which increases the stability of this barrier film However,such an adsorption behavior varies with respect to the charge state of molecule.As illustrated in the hydrogen evolution results in Fig.4,fumarate with double equal deprotonated carboxyl group can more easily absorb on the oxide surface in comparison with glycolate and benzoate employing only one carboxyl group.Nevertheless,it is reasonable to suggest that fumarate incorporated film registers higher thickness or compactness than that of the other candidates judging from the evolution of film capacitance(Fig.9b).However,the highest thickness cannot be equivalent to the best inhibition performance in current situation since the resistance of adsorptive inhibitor layer against corrosion attack can be also affected by molecular structure and properties of inhibitor.Based on this thought,benzoate employing lower thickness exhibits higher and stable impedance since it tends to be deposited in plane with the surface due to the uniform distribution ofπbonds in benzene ring may also form bonds with hydroxides,which demonstrates higher surface coverage than the other two inhibitors in the same concentration.Despite of the adsorption of inhibitors on the oxide surface,hydrolytical degradation of the oxide would also occur with time.Since then,the inhibitors tend to be involved during the formation of hydroxide products.Chemisorption of fumarate was found as a layer in the middle of corrosion product(s) via in situ Raman by Maltseva et al.[22] (but was not confirmed by in situ ATR-FTIR characterization [23]),who suggested fumarate improved the corrosion performance by forming bridging bidentate carboxylate coordination complex that enhances the barrier property of the corrosion product(s)layer.Moreover,function of inhibitors may be assigned during the hydroxide formation through considering their ion complexing capability Mg2+and Al3+(from the dissolution of Al enrichedα-Mg),as shown in Fig.12b.The stability constant logKfor Mg2+-glycolate complex is 0.92 (ML/M.L@ 25 °C,0.1 ionic strength),whilst that for Mg2+-benzoate complex is 0.1 (ML/M.L @ 30 °C,0.4) and Al3+-benzoate complex is 12.09(MOHL/M.OH.L@25°C,0.5)[42],which indicate the formation of Men+-Inhibitor complex(s) may also contribute the enhancement of corrosion.Unfortunately,this assumption cannot experimentally be envisaged in the present study due to the limit in skill to characterize the trace amount of complex in aqueous solution,which deserves a further investigation.

    5.Conclusions

    In the present work three characteristic carboxylates were studied as corrosion inhibitors for reducing the corrosion rate of AM50 Mg alloy in neutral NaCl solution.The optimum concentrations of inhibitors were preselected via hydrogen evolution measurements.Experimental (including electrochemical methods and surface analysis techniques) and theoretical calculations (DFT) were employed to understand their inhibition mechanism.Based on the obtained results,three critical conclusions can be drawn:

    1.The sodium salts of glycolate,fumarate and benzoate can effectively reduce the corrosion rate of AM50 in 0.5 wt.%NaCl solution and their inhibition efficiency generally increases with concentration and immersion time.

    2.All investigated carboxylates reveal anodic corrosion inhibition action,and their inhibition behaviors are ascribed to initial physisorption on the passivation and hydr(oxide)surface,followed by chemisorption with Mg2+and Al3+.

    3.Adsorption of carboxylate inhibitors on the surface of corroding AM50 alloy mainly occurs at the carboxyl (electron donating) and hydroxyl group,C=C bond and benzene ring (electron accepting).Increasing the number of equally deprotonated carboxyl would facilitate the adsorption of carboxylate inhibitors,whilst the inhibition efficiency is mainly determined by the stability of complex formed between inhibitor and ions formed during corrosion.

    Declaration of Competing Interest

    The authors declare that they have no conflict of interest to this work.We declare that we do not have any commercial or associative interest that represents a conflict of interest in connection with the work submitted.

    Acknowledgment

    Junjie Yang appreciates the financial support by Guangdong Basic and Applied Basic Research Founding (Grand No.2020A1515110754),MOE Key Lab of Disaster Forest and Control in Engineering,Jinan University (Grand No.20200904008),Educational Commission of Guangdong Province (Grand No.2020KTSCX012),the Fundamental Research Funds for Central Universities (Grand No.21620342)and the High Performance Public Computing Service Platform of Jinan University.Xiaopeng Lu acknowledges the financial support from National Natural Science Foundation of China(Grand No.52071067) and the Fundamental Research Funds for the Central Universities (Grand No.N2002009).

    Data availability

    The raw/processed data required to reproduce these findings cannot be shared at this time as the data also forms part of an ongoing study.

    淫妇啪啪啪对白视频| 成人亚洲精品av一区二区 | 嫩草影视91久久| 日本欧美视频一区| 久久久久久久精品吃奶| 狂野欧美激情性xxxx| 高潮久久久久久久久久久不卡| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲中文字幕日韩| 免费av毛片视频| 亚洲黑人精品在线| 男女高潮啪啪啪动态图| 亚洲专区国产一区二区| 精品人妻在线不人妻| 天堂中文最新版在线下载| 一级作爱视频免费观看| 色播在线永久视频| 久久欧美精品欧美久久欧美| 亚洲性夜色夜夜综合| 视频在线观看一区二区三区| 91国产中文字幕| 亚洲精品中文字幕一二三四区| 1024视频免费在线观看| 久久久国产成人免费| xxx96com| 国产欧美日韩一区二区三区在线| 久久久久久大精品| 人人妻人人澡人人看| 国产精品一区二区在线不卡| 97碰自拍视频| 久久久精品欧美日韩精品| 在线观看免费午夜福利视频| 中文字幕高清在线视频| 999久久久国产精品视频| 亚洲成av片中文字幕在线观看| 免费看十八禁软件| 黄色a级毛片大全视频| 亚洲熟女毛片儿| 在线永久观看黄色视频| 久久精品亚洲av国产电影网| 国产激情欧美一区二区| 婷婷精品国产亚洲av在线| 欧美日韩中文字幕国产精品一区二区三区 | 国产亚洲精品久久久久5区| 免费在线观看亚洲国产| 国产伦一二天堂av在线观看| 人人妻,人人澡人人爽秒播| 久久久久久久午夜电影 | 亚洲国产精品999在线| 国产精品自产拍在线观看55亚洲| 极品人妻少妇av视频| 美女福利国产在线| 自线自在国产av| 久久香蕉国产精品| 免费在线观看黄色视频的| 欧美最黄视频在线播放免费 | 亚洲精华国产精华精| 免费av中文字幕在线| 国产成+人综合+亚洲专区| 国产av在哪里看| 日韩视频一区二区在线观看| av网站在线播放免费| 一本综合久久免费| 久久亚洲精品不卡| 亚洲国产欧美一区二区综合| 久久人妻av系列| 成人特级黄色片久久久久久久| 国产一卡二卡三卡精品| 精品熟女少妇八av免费久了| netflix在线观看网站| 亚洲伊人色综图| 亚洲狠狠婷婷综合久久图片| 91麻豆av在线| 又黄又爽又免费观看的视频| 亚洲激情在线av| a在线观看视频网站| 国产精品一区二区在线不卡| 老汉色∧v一级毛片| 亚洲精品美女久久av网站| 别揉我奶头~嗯~啊~动态视频| av国产精品久久久久影院| 亚洲国产欧美日韩在线播放| 国产av精品麻豆| 欧美黄色淫秽网站| 欧美性长视频在线观看| 婷婷六月久久综合丁香| 欧美精品一区二区免费开放| 怎么达到女性高潮| 村上凉子中文字幕在线| 电影成人av| 中文字幕最新亚洲高清| 久久久久久人人人人人| 欧美日韩亚洲综合一区二区三区_| 免费不卡黄色视频| 亚洲色图av天堂| 国产av一区二区精品久久| 18美女黄网站色大片免费观看| av天堂久久9| 精品国产超薄肉色丝袜足j| 他把我摸到了高潮在线观看| 国产欧美日韩综合在线一区二区| 欧美性长视频在线观看| 午夜福利,免费看| 久久午夜亚洲精品久久| 美女福利国产在线| 亚洲精品在线观看二区| 久久亚洲真实| 丰满迷人的少妇在线观看| 国产亚洲精品久久久久久毛片| 欧美老熟妇乱子伦牲交| 久久久精品国产亚洲av高清涩受| 欧美日韩精品网址| 成年人免费黄色播放视频| 国产欧美日韩综合在线一区二区| 国产精品乱码一区二三区的特点 | 人成视频在线观看免费观看| 最近最新中文字幕大全电影3 | 国产99久久九九免费精品| 新久久久久国产一级毛片| 亚洲性夜色夜夜综合| 精品久久蜜臀av无| 亚洲一区中文字幕在线| 欧美黑人欧美精品刺激| 人人妻人人澡人人看| 国产成+人综合+亚洲专区| 一区福利在线观看| 在线观看免费午夜福利视频| a级毛片黄视频| 国产精华一区二区三区| 国产熟女xx| 色在线成人网| 中国美女看黄片| 黄色怎么调成土黄色| 亚洲 欧美 日韩 在线 免费| 一进一出抽搐gif免费好疼 | 欧美中文日本在线观看视频| 一级毛片高清免费大全| 亚洲国产欧美日韩在线播放| 日本黄色日本黄色录像| 十八禁人妻一区二区| 可以免费在线观看a视频的电影网站| 丁香六月欧美| 老司机在亚洲福利影院| 国产在线观看jvid| 他把我摸到了高潮在线观看| 亚洲精品在线美女| 99久久人妻综合| 国产精品国产av在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 日韩欧美在线二视频| 人人妻人人添人人爽欧美一区卜| a级毛片黄视频| 伊人久久大香线蕉亚洲五| 深夜精品福利| 啪啪无遮挡十八禁网站| 日韩一卡2卡3卡4卡2021年| 黄色a级毛片大全视频| 国产精品免费视频内射| 国产三级在线视频| 亚洲精品久久成人aⅴ小说| 香蕉丝袜av| 热99re8久久精品国产| 免费在线观看日本一区| 精品国产超薄肉色丝袜足j| 涩涩av久久男人的天堂| 久久香蕉国产精品| 国产精品偷伦视频观看了| 欧美中文综合在线视频| 精品国内亚洲2022精品成人| 亚洲精品一卡2卡三卡4卡5卡| 99国产极品粉嫩在线观看| 亚洲成av片中文字幕在线观看| av网站免费在线观看视频| 中出人妻视频一区二区| 一级毛片精品| 亚洲 欧美 日韩 在线 免费| 午夜a级毛片| 日本免费a在线| 性少妇av在线| 国产伦一二天堂av在线观看| 亚洲精品在线观看二区| 国产欧美日韩精品亚洲av| 青草久久国产| 亚洲欧美日韩高清在线视频| 侵犯人妻中文字幕一二三四区| 日韩 欧美 亚洲 中文字幕| 精品免费久久久久久久清纯| 人人妻人人澡人人看| 一级毛片精品| 18禁观看日本| 宅男免费午夜| 亚洲精华国产精华精| 一级片'在线观看视频| 香蕉久久夜色| 日韩精品青青久久久久久| 成年人免费黄色播放视频| 成人手机av| 高清在线国产一区| 国产精品98久久久久久宅男小说| 精品一区二区三卡| 日本免费a在线| 亚洲一区高清亚洲精品| 久久国产亚洲av麻豆专区| 悠悠久久av| 亚洲第一青青草原| 丝袜人妻中文字幕| 91大片在线观看| 国产精品98久久久久久宅男小说| 村上凉子中文字幕在线| 国产一区二区三区在线臀色熟女 | 午夜两性在线视频| 男女做爰动态图高潮gif福利片 | 午夜福利一区二区在线看| 国产免费现黄频在线看| 亚洲人成伊人成综合网2020| 人妻丰满熟妇av一区二区三区| 国产精品99久久99久久久不卡| 国产又色又爽无遮挡免费看| 免费女性裸体啪啪无遮挡网站| 三级毛片av免费| 99久久99久久久精品蜜桃| 日本免费一区二区三区高清不卡 | 免费观看人在逋| e午夜精品久久久久久久| 午夜91福利影院| 久久国产乱子伦精品免费另类| 国产一区二区三区在线臀色熟女 | av天堂在线播放| 十八禁网站免费在线| 国产区一区二久久| 纯流量卡能插随身wifi吗| 欧美国产精品va在线观看不卡| 黄片小视频在线播放| 极品人妻少妇av视频| 黄网站色视频无遮挡免费观看| 国产一区二区三区在线臀色熟女 | 亚洲狠狠婷婷综合久久图片| xxxhd国产人妻xxx| 一a级毛片在线观看| 国产国语露脸激情在线看| 国产精品免费视频内射| 亚洲一码二码三码区别大吗| 韩国av一区二区三区四区| 99精品欧美一区二区三区四区| 精品久久蜜臀av无| 最好的美女福利视频网| 超色免费av| 超碰成人久久| 亚洲人成77777在线视频| 搡老乐熟女国产| 亚洲国产精品999在线| 久久久国产成人免费| 色婷婷久久久亚洲欧美| 一进一出好大好爽视频| 如日韩欧美国产精品一区二区三区| 国产亚洲欧美精品永久| www日本在线高清视频| 国产不卡一卡二| a在线观看视频网站| 免费在线观看日本一区| 国产高清videossex| 亚洲一区二区三区不卡视频| 久久欧美精品欧美久久欧美| 久久精品成人免费网站| 热re99久久国产66热| av在线播放免费不卡| 制服诱惑二区| 黑人操中国人逼视频| 精品国产美女av久久久久小说| 在线观看午夜福利视频| 免费一级毛片在线播放高清视频 | 夜夜看夜夜爽夜夜摸 | 无人区码免费观看不卡| 亚洲一区二区三区不卡视频| 久久国产精品男人的天堂亚洲| 精品福利永久在线观看| 国内久久婷婷六月综合欲色啪| 亚洲精品成人av观看孕妇| 在线免费观看的www视频| 午夜精品久久久久久毛片777| 9191精品国产免费久久| 久久久国产欧美日韩av| 51午夜福利影视在线观看| 人妻久久中文字幕网| 不卡av一区二区三区| 男女下面进入的视频免费午夜 | 欧美国产精品va在线观看不卡| 国产又爽黄色视频| 亚洲一区二区三区欧美精品| 久久久久久久久免费视频了| 999精品在线视频| 嫩草影院精品99| 免费久久久久久久精品成人欧美视频| 国产精品98久久久久久宅男小说| 亚洲男人天堂网一区| 亚洲人成电影免费在线| 午夜日韩欧美国产| 搡老岳熟女国产| 别揉我奶头~嗯~啊~动态视频| 精品少妇一区二区三区视频日本电影| 日本免费a在线| 在线国产一区二区在线| 久久久久久久午夜电影 | 香蕉丝袜av| 成人精品一区二区免费| 美女高潮到喷水免费观看| 亚洲在线自拍视频| 国产成年人精品一区二区 | 国产一区二区激情短视频| 最近最新中文字幕大全电影3 | 国产亚洲av高清不卡| 满18在线观看网站| 国产精品国产高清国产av| 日韩欧美免费精品| 中文字幕精品免费在线观看视频| 少妇裸体淫交视频免费看高清 | 中文字幕人妻丝袜一区二区| 亚洲在线自拍视频| 国产欧美日韩精品亚洲av| 久久影院123| 在线永久观看黄色视频| 美女扒开内裤让男人捅视频| av欧美777| 亚洲第一欧美日韩一区二区三区| 国产精品久久久av美女十八| 国产97色在线日韩免费| 在线视频色国产色| 久久久久国产一级毛片高清牌| 精品人妻1区二区| 99久久国产精品久久久| 自线自在国产av| 精品熟女少妇八av免费久了| 男人舔女人的私密视频| 99久久国产精品久久久| 国产精品影院久久| 国产区一区二久久| 亚洲国产精品sss在线观看 | 亚洲av成人不卡在线观看播放网| 99久久国产精品久久久| 亚洲伊人色综图| 99久久国产精品久久久| 久久久久亚洲av毛片大全| 国产欧美日韩综合在线一区二区| ponron亚洲| 别揉我奶头~嗯~啊~动态视频| 亚洲精品中文字幕在线视频| 别揉我奶头~嗯~啊~动态视频| 丝袜在线中文字幕| 国产精品美女特级片免费视频播放器 | 国产欧美日韩一区二区三区在线| 女警被强在线播放| 99国产极品粉嫩在线观看| 国产精品秋霞免费鲁丝片| 午夜视频精品福利| 久久这里只有精品19| 三上悠亚av全集在线观看| 后天国语完整版免费观看| 妹子高潮喷水视频| 后天国语完整版免费观看| 国产一区二区三区视频了| 国产色视频综合| 妹子高潮喷水视频| 欧美日韩福利视频一区二区| 乱人伦中国视频| 黄色片一级片一级黄色片| 男人舔女人的私密视频| 日本撒尿小便嘘嘘汇集6| 亚洲伊人色综图| 色婷婷av一区二区三区视频| 免费看a级黄色片| 别揉我奶头~嗯~啊~动态视频| 欧美国产精品va在线观看不卡| 成人免费观看视频高清| 一级毛片精品| 成在线人永久免费视频| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美精品一区二区免费开放| 欧美乱色亚洲激情| 国产xxxxx性猛交| 国产色视频综合| 国产亚洲精品一区二区www| 国产高清激情床上av| 18禁黄网站禁片午夜丰满| 男人的好看免费观看在线视频 | 色老头精品视频在线观看| 在线免费观看的www视频| 性欧美人与动物交配| 亚洲精品成人av观看孕妇| 亚洲狠狠婷婷综合久久图片| 国产亚洲欧美精品永久| 日本免费一区二区三区高清不卡 | 香蕉国产在线看| 国产免费男女视频| 88av欧美| 可以免费在线观看a视频的电影网站| 亚洲中文字幕日韩| 嫩草影院精品99| 丰满迷人的少妇在线观看| 天天躁狠狠躁夜夜躁狠狠躁| av超薄肉色丝袜交足视频| 亚洲中文日韩欧美视频| 黄片小视频在线播放| 亚洲男人的天堂狠狠| 亚洲少妇的诱惑av| 欧美色视频一区免费| 国产精品1区2区在线观看.| 99精国产麻豆久久婷婷| 国产99久久九九免费精品| 久久亚洲精品不卡| 国产亚洲欧美98| 波多野结衣高清无吗| 亚洲五月婷婷丁香| 99国产精品一区二区三区| 日韩三级视频一区二区三区| 国产精品永久免费网站| 99国产极品粉嫩在线观看| 在线观看午夜福利视频| 亚洲自偷自拍图片 自拍| 亚洲精品美女久久久久99蜜臀| 大码成人一级视频| 欧美日韩亚洲高清精品| 免费人成视频x8x8入口观看| 久久 成人 亚洲| 一进一出好大好爽视频| 亚洲人成伊人成综合网2020| 一边摸一边做爽爽视频免费| 中文字幕另类日韩欧美亚洲嫩草| 日韩大尺度精品在线看网址 | 狠狠狠狠99中文字幕| 欧美日韩精品网址| 国产色视频综合| 亚洲国产欧美日韩在线播放| bbb黄色大片| 涩涩av久久男人的天堂| 高清av免费在线| 在线观看免费视频网站a站| svipshipincom国产片| 欧美老熟妇乱子伦牲交| 国产精品99久久99久久久不卡| av在线天堂中文字幕 | 一区二区日韩欧美中文字幕| 日本撒尿小便嘘嘘汇集6| 免费在线观看影片大全网站| 男人的好看免费观看在线视频 | 精品久久久久久,| 国产极品粉嫩免费观看在线| 女同久久另类99精品国产91| 性色av乱码一区二区三区2| 精品高清国产在线一区| 老汉色av国产亚洲站长工具| 国产精品一区二区在线不卡| 午夜a级毛片| 可以免费在线观看a视频的电影网站| 神马国产精品三级电影在线观看 | 久久热在线av| 午夜a级毛片| 狂野欧美激情性xxxx| 亚洲精品一区av在线观看| 99久久国产精品久久久| 亚洲 国产 在线| 欧美丝袜亚洲另类 | 国产伦人伦偷精品视频| 精品少妇一区二区三区视频日本电影| av福利片在线| 一级a爱片免费观看的视频| 国产av一区在线观看免费| 91精品三级在线观看| 亚洲午夜理论影院| 母亲3免费完整高清在线观看| 精品一区二区三区视频在线观看免费 | 午夜免费鲁丝| 久久人妻福利社区极品人妻图片| 亚洲精品一卡2卡三卡4卡5卡| 午夜精品在线福利| 国产男靠女视频免费网站| 亚洲精品国产一区二区精华液| 免费在线观看完整版高清| 91精品国产国语对白视频| 中文字幕另类日韩欧美亚洲嫩草| 黑人巨大精品欧美一区二区mp4| 久久香蕉激情| 日韩成人在线观看一区二区三区| 两人在一起打扑克的视频| 99久久人妻综合| 人人妻人人澡人人看| 欧美一区二区精品小视频在线| 欧美日韩视频精品一区| 亚洲精品国产区一区二| 欧美色视频一区免费| 欧美日韩黄片免| 中文欧美无线码| 国产亚洲精品久久久久久毛片| 国产精品 国内视频| 亚洲色图av天堂| 一个人观看的视频www高清免费观看 | 成年人免费黄色播放视频| 一级毛片女人18水好多| 啦啦啦 在线观看视频| 脱女人内裤的视频| 女生性感内裤真人,穿戴方法视频| 国产激情欧美一区二区| 久久久国产精品麻豆| 热re99久久国产66热| 露出奶头的视频| 精品国产国语对白av| 99精国产麻豆久久婷婷| 亚洲自偷自拍图片 自拍| 天堂中文最新版在线下载| 亚洲va日本ⅴa欧美va伊人久久| 精品免费久久久久久久清纯| 日韩精品中文字幕看吧| 两性夫妻黄色片| 亚洲成人精品中文字幕电影 | 午夜a级毛片| av国产精品久久久久影院| 久久精品国产亚洲av高清一级| 妹子高潮喷水视频| 老司机午夜福利在线观看视频| 亚洲av美国av| 超碰成人久久| 麻豆国产av国片精品| 成人三级黄色视频| 亚洲精品久久午夜乱码| 色婷婷久久久亚洲欧美| 国产欧美日韩一区二区精品| 免费少妇av软件| 免费观看精品视频网站| 男女下面插进去视频免费观看| 老司机靠b影院| 国产亚洲欧美精品永久| av欧美777| 国内久久婷婷六月综合欲色啪| 久久精品国产清高在天天线| 久久午夜综合久久蜜桃| 99精品欧美一区二区三区四区| 精品国产乱子伦一区二区三区| 国产精品自产拍在线观看55亚洲| 亚洲人成网站在线播放欧美日韩| 两人在一起打扑克的视频| a级片在线免费高清观看视频| 欧美日韩国产mv在线观看视频| 午夜福利在线观看吧| 人人妻,人人澡人人爽秒播| 男女高潮啪啪啪动态图| 成在线人永久免费视频| 成人特级黄色片久久久久久久| 成人国语在线视频| 亚洲男人天堂网一区| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美日韩av久久| 黄色片一级片一级黄色片| 少妇被粗大的猛进出69影院| 精品卡一卡二卡四卡免费| 嫁个100分男人电影在线观看| 亚洲美女黄片视频| 亚洲欧美日韩高清在线视频| 在线看a的网站| 欧美一区二区精品小视频在线| 欧美中文综合在线视频| 免费高清在线观看日韩| 亚洲 欧美 日韩 在线 免费| 两人在一起打扑克的视频| 国产不卡一卡二| 精品第一国产精品| 午夜福利一区二区在线看| 叶爱在线成人免费视频播放| 一本综合久久免费| 男女下面插进去视频免费观看| 一级a爱视频在线免费观看| 丝袜人妻中文字幕| 国产一区二区三区视频了| 在线天堂中文资源库| 成熟少妇高潮喷水视频| 日日摸夜夜添夜夜添小说| 日韩欧美一区视频在线观看| 狂野欧美激情性xxxx| 久久亚洲精品不卡| 久久人人爽av亚洲精品天堂| 欧美激情极品国产一区二区三区| 国产99久久九九免费精品| 亚洲人成77777在线视频| 女性被躁到高潮视频| 日韩高清综合在线| 一进一出抽搐gif免费好疼 | 亚洲人成网站在线播放欧美日韩| 久久久国产成人免费| 女人高潮潮喷娇喘18禁视频| 99久久综合精品五月天人人| 亚洲片人在线观看| xxx96com| 久久久久久久午夜电影 | 午夜91福利影院| 久久久久久久久久久久大奶| 看黄色毛片网站| 最近最新中文字幕大全免费视频| 亚洲av五月六月丁香网| 午夜精品在线福利| 久久人人97超碰香蕉20202| 久久精品国产99精品国产亚洲性色 | 丰满迷人的少妇在线观看| 日韩三级视频一区二区三区| 韩国av一区二区三区四区| av网站免费在线观看视频| 不卡一级毛片| 色婷婷av一区二区三区视频| 亚洲精品久久成人aⅴ小说| 精品国产乱子伦一区二区三区| 亚洲性夜色夜夜综合| 九色亚洲精品在线播放| 成人黄色视频免费在线看| 亚洲欧美日韩无卡精品| 国产aⅴ精品一区二区三区波| 欧美日韩中文字幕国产精品一区二区三区 | 在线观看免费日韩欧美大片| 亚洲av成人不卡在线观看播放网|