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

    Investigation of Propoxylation and Ethoxylation of 2-Ethylhexanol in the Presence of an Alkaline and DMC Type Catalyst at Initial Stages of the Syntheses

    2016-07-05 10:15:34WieslawHreczuchKarolinaDabrowskaArkadiuszChruciel
    China Detergent & Cosmetics 2016年1期

    Wieslaw Hreczuch, Karolina Dabrowska, Arkadiusz Chru?ciel

    MEXEO, Institute of Technology, Poland

    Katarzyna Materna, Agata Sznajdrowska

    Institute of Chemical Technology and Engineering, Poznan University of Technology, Poland

    Krystyna Czaja

    Faculty of Chemistry, Opole University, Poland

    Wanxu Wang, Qiuxiao Li, Yongqiang Sun

    China Research Institute of Daily Chemical Industry, China

    Introduction

    Surfactants have been applied in many products of the chemical industry, such as detergents, paints, dyestuffs,paper coatings, inks, plastics and fibers, personal care and cosmetics, agrochemicals, pharmaceuticals, food processing,etc. In addition, they play an important role in the oil industry. Hence, it seems to be extremely important to search for the appropriate sources of raw materials for their production. Such proposals include 2-ethylhexyl alcohol.[1-3]The compound is available on the market at competitive prices, in comparison with the detergent range alcohols.The discussed 2-ethylhexanol (EH) derived surfactants are synthesized according to the general formula (1):

    Where n is average ethoxylation grade, typically up to 15.

    Generally, homogeneous and heterogeneous type catalysts are used in the alkoxylation process. The former is represented by KOH which is commonly applied in the synthesis and leads to an undesirable, wide range homolog distribution of the products generated. The latter can be represented by the dimetalcyanide (DMC) type catalyst, a recently developed coordinative type propoxylation catalyst,which also attracts attention in ethoxylation.[4,5]

    The criteria for selecting an optimal catalyst include its kinetic performance in the studied synthesis, as well as composition and derived properties of the products obtained.

    Kinetic models of polyaddition

    The obtained polyaddition products are usually characterized by the average number of moles of the applied alkyloxirane (alkylene oxide) converted in the reaction per mole of the starter material. In fact, the synthesis follows a pattern of parallel, consecutive, second order reactions,leading to a polydispersed mixture of homologs within a certain range of polyaddition degrees around the average alkoxylation degree and some amount of the substrate usually remains unconverted. Generally, this is represented by the following set of kinetic Equations (2).

    Where k0, k1, ..., kndenote rate constants for the successive steps of the process, and X stands for oxygen, in the ethoxylation of an alcohol or another intermediate group,such as amine, amide or carboxyl.

    When the concentrations of the following homologs:RXH, RX(AO)H,..., RX(AO)nH are denoted as x0,x1, ..., xn, the alkyloxirane concentrations are described by xAOand the kinetic equations of the second order for the subsequent reactions are considered, then the reaction rates at the consecutive propagation steps can be expressed as follows (3) :

    Dividing these equations by the reaction rate at the first addition step (dx0/dt) and defining the distribution coefficients cn as the ratio of the successive rate constants:c1= k1/k0, c2= k2/k0,..., cn= kn/k0provides the following set of Equations (4) :

    Where the concentration xn is expressed as functions of the content of the x0(ROH) reagent. The above set of equations can be solved numerically, e.g., using the numerical integration technique of the Runge-Kutta method,optimization algorithms of the Monte Carlo method and the simplex Nelder & Mead method.[6,7]

    The Natta-Mantica Equation[8,9]can also be used for determination of the distribution coefficients(ci) (5):

    On the other hand, knowing an experimentally determined fractional content of an alkoxylation product(xi) and using the above Natta Mantica equation, one can compute the distribution constants (ci) and get results which are comparable to those of the analytical method. To that end, a dedicated computer software is required — and has also been developed and used by the present authors.[10]

    An intermediate simplified model was also reported. The amount of unconverted substrate and average polyaddition degree enable determination of the so-called Weibull-Nycander-Gold coefficient (WNG). It shows the ratio of the addition rate constant at the first step of the reaction and the following propagation steps, which are presumed equivalent(6):

    Where: kidenotes the reaction rate constant at the i-th addition step, whereas k1=k2=....kn, k0is the reaction rate constant at the first addition step, means average polyaddition degree and x0stands for the mole fraction of the unconverted substrate.

    With the cicoefficient, the corresponding fractional content model can be generated according to the following WNG Equation (7):

    The most general kinetic model reflecting the Poisson type distribution was described by Flory[11,12], with the assumption that an equivalent reaction rate occurs at the first and each of the following propagation steps (k0=k1=k2......kn)(8):

    The WNG distribution is much closer to that determined experimentally, in comparison with the Flory model.Moreover, the higher the value of CWNG=k/k0or ci=ki/k0or ci(NM),the broader the resulting homolog distribution, which is also shown in Figure 1.

    Figure 1. Illustration of the influence of the CWNG coefficient on the homolog distribution in polyaddition products

    Assuming accurate quantitative analysis results, the cicoefficients of the NM model reflect exactly the kinetic parameters (kik0); this has to take place in order to fix the results determined experimentally. The other theoretical model of WNG is also useful in practice for quantifying any existing difference in homolog distributions among the products compared, by contrasting the calculated CWNGvalues.

    Oxiranes readily react with diverse compounds by ring opening. Their typical reactions are those with nucleophiles,proceeding via the SN1 or SN2 mechanism both in acidic (weak nucleophiles: water, alcohols) and alkaline media (strong nucleophiles: OH—, RO—, NH3, RNH2, RR'NH, etc.). The general reaction pattern is as follows (9):

    The alkaline mechanism leads to the broadening of the obtained homolog distribution in the product mixture,as described by the Weibull-Tornquist effect,[13-16]through the semi-crown complexation of the alkaline cation by the growing oligo-ether chains Figure 2.

    Figure 2. Semi-crown complex of an alkaline metal active center

    Loosening the transient alcoholate-metal cation bond facilitates the connection with the next alkoxy-mere by decreasing the energetic barrier. It was further proved by calculating the distribution coefficients ci=ki/k0, that this implies increased reaction rate constants at the following addition steps until i=6, resulting in a broader distribution of the products formed.[10]The effect is blocked by bivalent or trivalent cations of the catalysts, yielding a narrower homolog mixture, that can also be attributed to the bivalent zinc active center of the complex DMC type catalyst.

    Despite the homolog distribution of the target polyaddition products, the applied catalyst is During alkoxylation, diols can be formed as by-products by the catalytic hydration of alkylene oxide (10):

    Where, R denotes Hor CH3.

    Further propagation generates dialkylene glycol,trialkylene glycol and higher polymers.

    When alkylene oxide is heated to about 150-300 ℃ in the presence of a catalyst, it can isomerize to form acetaldehyde(11):

    Moreover, it is a well known fact that, by eliminating the proton from the methyl group, methyloxirane can be converted into an allyl alcohol (12):

    Which is also capable of polyether chain propagation.

    The reactions described above affect the selectivity of the discussed polyaddition process as well as the quality and properties of the obtained products.

    A substantial difference appears in the structures of the symmetric oxirane and non-symmetric methyloxirane,which leads to isomerisation and steric hindrance effects during propagation. Sometimes, the methyl-branched structure and the formed secondary hydroxyl groups of poly-propoxy chains are responsible for their limited reactivity with comparison with the primary and linear polyethers formed by the oxirane addition.Therefore, ethoxylation and propoxylation should be considered separately, especially when they are applied in a production process one by one.

    The aim of this work is to compare the activity of a conventional alkaline (KOH) and a DMC type catalyst in the propoxylation and ethoxylation of 2-ethylhexanol.Moreover, the kinetics and selectivity at the initial stage of the syntheses were investigated. These parameters are important for the following propagation steps and properties of the final products.

    Experimental section

    Materials

    2-ethylhexanol was commercially produced by the K?dzierzyn Nitrogen Works, Azoty Group, K?dzierzyn-Ko?le, Poland. Ethylene oxide was obtained from PKN Orlen, P?ock, Poland. Propylene oxide was manufactured by PCC Rokita S.A., Brzeg Dolny, Poland. KOH, 85% pure,was purchased from Brenntag, K?dzierzyn-Ko?le, Poland.The applied DMC type catalyst[4]was manufactured by MEXEO, K?dzierzyn-Ko?le, Poland. NM-CAD, the computer program used for the calculation of kinetic parameters in polyaddition reactions, was developed by MEXEO.

    Methods

    Preparation of the DMC catalyst. K3[Co(CN)6] (7.5 g)was dissolved in a round-bottom flask, filled with a mixture of distilled water (40 mL) and tBuOH (50 mL) (Solution A). Solution B was prepared by dissolving ZnCl2(75 g) in a mixture of water (75 mL) and tBuOH (2 mL). Solution B was added to Solution A over 30 min with stirring using a mechanical stirrer. After combining them, stirring was continued for another 15min. 200 mL of a 1% aqueous solution of acrylic acid polyamide, a flocculent which enables a precipitate to be separated from aqueous solutions, was added to the prepared reaction mixture and stirred for 10 min more. Precipitation of a virtually amorphous precipitate of zinc hexacyanocobaltate(III) followed. The mixture was centrifuged and the resulting catalyst cake was dissolved in a mixture of tBuOH (130 mL) and 1% aqueous solution of the flocculent (55 mL). The reaction mixture was stirred vigorously for 15 min. The described operation of rinsing and filtration was performed once more and the obtained catalyst solid was dried to a constant weight at 50 ℃.[4,5]

    Polyaddition of methyloxirane/oxirane was carried out by means of a 1-L laboratory autoclave equipped with a programmable logic control system (PLCS)controlling its operation, as well as mechanical stirrer,heating jacket and cooling coil. The reactor was charged with an initiator and catalyst and dehydration followed through nitrogen purge over 30 min, at 130 ℃. The syntheses were performed through a one-step drop of the monomer into the starter material at a 1:1 molar ratio. In this case, the desired amount of the monomer was introduced at a time into a mixture of the starter with the catalyst and the reaction was continued at the desired temperature until overpressure in the reactor dropped to a constant level. Changes in overpressure reflected the kinetics of the monomer conversion.

    GPC analyses were performed using a GPC apparatus from Shimadzu, equipped with two LC-20AD pumps with a gradient, degassing unit Model DGU-20A3 and a fraction collector Model FRC-10A. Separation was carried out on the Shodex GPC column KF-805L. The flow rate of the eluent(tetrahydrofurane) was 1 cm3·min-1at a temperature of 40 ℃.The polymers’ molecular weights were found by calibration performed for polymethacrylate standards with narrowrange distributions of molecular weights (Sigma-Aldrich).The results were processed using the Shimadzu LCSolution software for determination of the molecular weights of the analyzed systems, their dispersity and purity.

    Gas chromatography (GC) was applied for analyzing the products’ compositions, using a gas chromatograph from Hewlett-Packard with flame ionization detector,equipped with a CHROMA computer integrator and a DB-5HT capillary column with the parameters 0.53 mm×15 m×0.1 μm. The analysis conditions: the column temperature programmed by 1min isotherm with an initial rate of 8 ℃ · min-1in the range of 50-340 ℃, the injector temperature was 340 ℃, and detector temperature was 350 ℃. The carrier gas (helium) flow rate was 2 mL · min-1.Interpretation of the results was performed using the internal standardization methodology.

    Results and discussion

    The experimental products examined in this work were manufactured by monopropoxylation of 2-ethylhexanol followed by its monoethoxylation. PO conversion was investigated in the presence of KOH and DMC catalysts.One mole of PO per one mole of 2-EH and catalyst was added to the autoclave feed at an ambient temperature. Then the reaction temperature was gradually increased. In both cases, initiations of the syntheses were observed at ca. 120℃ and 100 ℃ with KOH and DMC catalyst, respectively.From that moment on, overpressure began to decrease. The comparison of the initiation temperature and the kinetics of conversion during monopropoxylation of EH in the presence of KOH and DMC as the catalysts, where the catalyst concentration was 0.3% and 0.05% per product weight,respectively, is shown in Figure 3.

    The time of reduction of overpressure reflected the average rate of PO conversion. The amount of PO introduced into the reactor was consumed within 120min and 10min in the presence of KOH and DMC catalyst, respectively. Hence,the DMC catalyst’s activity was much higher, despite its significantly lower concentration in the reaction medium.

    The obtained products were further subjected to ethoxylation by 1 mole of EO per mole of EH-PO using the catalyst residue remaining in the starter. Dynamic trends from the synthesis are compared in Figure 4.

    The same initiation temperature was determined at c.-a.130℃ for both catalysts and the time of EO conversion was 80 and 10min, with KOH and DMC catalyst, respectively.

    Figure 3. The initiation temperature and the kinetics of conversion during monopropoxylation of EH in the presence of 0.3% KOH and 0.05% DMC as the catalysts (the red line depicts the reactor temperature, the green line – the overpressure, the blue line – the weight of PO in the feeding vessel)

    Figure 4. The initiation temperature for monoethoxylation of EH-PO in the presence of 0.3 % KOH and 0.05% DMC as the catalysts (the red line depicts the reactor temperature, green line - the overpressure, blue line - weight of PO in the feeding vessel)

    The differences in the kinetics of the compared PO conversion in the presence of KOH and DMC type catalyst can be quantified by calculating the average reaction rate factor Rp, expressed as g of PO consumed per minute, per g of catalyst and per g of starter, respectively Equation (13):

    A comparison of the calculated average reaction rate factors is shown in Table 1.

    Table 1. Average reaction rate factor Rp for the EH-PO and EH-PO-EO products

    A significant kinetic advantage of the DMC type catalyst over the compared alkaline catalyst was found in both of the cases shown in Table 1. Conversion was 60 to 70 times as fast in the presence of the DMC type catalyst, in comparison with KOH.

    The following systems investigated represent a model where one mole of oxirane or methyloxirane reacted directly with one mole of 2-ethylhexanol and the reaction products marked as EH-PO and EH-EO, respectively, were formed.Alkylene oxide was introduced into the reactor at the adjusted reaction temperature. This way, the influence of the catalyst’s concentration and the reaction temperature on the synthesis performance was examined. The obtained results are shown in Figure 5 and Figure 6.

    Figure 5. The reaction rate factor Rp depending on the catalyst type and its concentration, at the initial stage of the reaction(starter/monomer=1 mole/1 mole)

    The influence of the catalyst type investigated in this work and its concentration on the reaction rate Rp, is presented in Figure 5. It can be observed that the DMC catalyst’s activity decreases along the increase in its concentration while, with KOH, kinetics seems stable, although at a relatively low Rpvalue.

    Further research was performed to investigate the influence of the reaction temperature on the synthesis performance at selected concentrations for both catalysts.The obtained results are presented in Figure 6.

    Figure 6. The influence of reaction temperature on reaction rates of propoxylation and ethoxylation for both catalyst types at the initial stage of the reaction (starter/monomer = 1 mole/1 mole), at 0.3% of KOH and 0.05% DMC, respectively

    It was observed that, after initiating the reaction between 80℃ and 100℃, the kinetics of EO and PO addition seemed weakly influenced by a further increase in the reaction temperature with the DMC catalyst, even though the values of Rpwere relatively high. In the case of KOH, starting from 100℃, the conversion speeded-up more significantly in the case of ethoxylation as compared with propoxylation.

    Such significant differences between the activities of the compared catalysts, expressed by the values of Rpfactors, can be accounted for by different reaction mechanisms, while the DMC catalyst coordinates the ring opening attack at alkylene oxide of very high reactivity in the addition reaction, as compared with that of SN2 nucleophilic substitution, creating an active center at the less reactive alcoholate anion, further activated by the increase in temperature. Obviously, oxirane exhibits higher activity than does methyloxirane because of the lower energy barrier of the epoxy-ring opening in this case. It should be further stated that the apparent decrease in the activity of the DMC catalyst connected with its increased concentration and stable reactivity along the increased temperature may only prove that the reaction progress criterion was shifted from its kinetic phase to the diffusion rate of gaseous alkyloxirane into the liquid reaction medium.This means that the reaction rate is faster than dissolution of the gaseous reagents in the liquid reaction medium. As the result, the reaction rate with the DMC catalyst is fast enough to convert instantly all the available monomer, just dissolved.Consequently, the apparent decrease in the Rpfactor can only be attributed to the result of mathematical calculation, in fact.Therefore, a more efficient mass transfer reactor-absorber system should be used for investigating the reaction kinetics in this case.

    The influence of the catalyst concentration at a constant reaction temperature (130℃) and that of the reaction temperature at a constant catalyst concentration (KOH, 0.3%and DMC, 0.05%) on the obtained homolog distributions, as reflected by the CWNGdistribution coefficients, is presented in Figure 7.

    Figure 7. Comparison of the CWNG=k/k0 coefficients for EH-PO, marked by the symbols Δ? and EHEO denoted as ▲?, obtained in the presence of KOH (diamond) and DMC catalyst (triangle) in relationship to a) catalysts concentrations and b)reaction temperatures, respectively

    The obtained homolog distributions depend on the applied adduct (PO or EO) and catalyst (KOH or DMC).No influence of the reaction temperature or the applied catalyst concentration was noted during propoxylation. It appears, that the steric hindrance effect of the branched structure of methyloxirane and the generated second order hydroxyl groups as the synthesis active center prevail over the influence of the reaction mechanism, conditioned by the investigated catalyst type and reaction temperature, which play a dominant role in ethoxylation.

    For the EH-EO products, the observed differences in the calculated values of CWNGcoefficients are quite significant,taking into account the very initial stage of propagation.The applied DMC type catalyst gives narrower homolog distributions in comparison with KOH, as reflected by the lower values of CWNGwhere the ratio of reaction rate constants at the propagation steps to that at the first addition step is relatively lower. In consequence, the DMC based catalysis yields also a relatively higher conversion of the starter in the obtained products. Furthermore, the increased reaction temperature brings an additional narrowing effect when the DMC catalyst is used for ethoxylation, which is not observed with KOH. This seems an interesting issue, derived from the substantial differences in the reaction mechanisms.Its explanation will be the subject of further research.

    Initially, gas chromatography was selected for determining the obtained product content of a reasonably low molecular weight (weight balance average Mw<200).

    The issued chromatograms are shown in Figure 8.

    Figure 8. Gas chromatograms for the EH-PO products in the presence of a) KOH and b) DMC catalyst, respectively

    State-of-the-art knowledge indicates that addition of methyloxirane generates two isomeric forms, resulting in the formation of a first-order or a second-order hydroxyl group at the end of the polyether chain. Both isomer homologs can have the same number of mers so their peaks overlap or appear very close each other in the gas chromatograms.Additionally, polyglycols of the same number of mers are formed as by-products, which appear very close. For the sake of simplicity, the peak areas which are very close each other,presumably of the same number of mers, were additively unified as a single homolog in the calculations. The fractional contents, as quantified from the chromatograms in Figure 8,and the kinetic parameters, as calculated from the fractional contents, are shown in Table 2.

    Table 2. Peak areas for monopropoxylates of ethylhexanol (EH-PO,Figure 8) and the derived kinetic parameters

    Quite a similar situation occurs when monoethoxylation of 2-ethylhexanol is compared (Figure 9, Table 3).

    Figure 9. Gas chromatograms for the EH-EO products in the presence of a) KOH and b) DMC catalyst, respectively

    Table 3. Peak area for ethoxylation of EH (EH-EO, Fig.9) and the derived kinetic parameters

    The average polyaddition degrees (ν) for the products from Table 2 and 3 should be theoretically equal to Table 1, as calculated from the synthesis wt. balance. The values calculated based on the GC results at the level below ν=1,reflect the reaction’s lower selectivity, in addition to a possible analytical error. As the result, the GC results discussed above,describe only part of the reaction product. Consequently, the KOH catalyst seems to be less selective in comparison with DMC, while the calculated average oligoaddition degrees are lower in that case. The substantially lower average polyaddition degrees determined by GC, in comparison with that expected from the synthesis weight balance,prompted additional analyses by the GPC technique,allowing determination of the product contents over a much wider range of molecular masses. The obtained products’chromatograms are shown in Figure 10 and 11.

    Figure 10. Comparison of GPC chromatograms obtained for the EH-PO1 products in the presence of KOH and DMC catalyst, respectively (Mw - molecular weight)

    Figure 11. Comparison of GPC chromatograms obtained for the EH-EO1 products in the presence of KOH and DMC catalyst, respectively

    Figure 10 and 11 show that additional products or,rather, by-products of much higher molecular weights were found using the GPC technique within the range from up to 460 and 960 Da in the case of polypropoxylates and polyethoxylates, respectively.

    Although GPC is not a very quantitative determination in terms of molecular weight content, the products’ average molecular masses can reasonably be interpreted qualitatively.Presumably, they represent the polymers described by Equation (10) and, to some extent, the side products described in Equations (11,12). It is worth mentioning that their total contribution seems higher than that of the target polyadducts of 2-ethylhexanol. This can be explained by a much higher polymerization rate, with two active centers per molecule, in comparison with that of the polyaddition reaction. Such a low selectivity of polyaddition, especially in the presence of KOH, was not reported before for relatively low molecular weights. However, it has not been expected and not detectable by the GC technique, commonly applied in this case.

    Conclusions

    It was found that the investigated DMC type catalyst is significantly more active in propoxylation and ethoxylation of ethylhexanol, in comparison with KOH. Additionally,it generates a more narrow homolog distribution and more selective products. However, in spite of the target polyaddition, both of the compared catalysts promote parallel polymerization consuming a significant part of the converted alkylene oxides.

    Acknowledgment

    This work was supported by the National Centre for Research and Development (NCBR): Applied Research Program 2012-2015 (PBS 1), Project No. PBS1/A5/28/2012(EH-SURF).

    [1] ECOSURFTM EH. Specialty Surfactants. The DOW Chemical Company 2009, 119-02244-0309 AMS.

    [2] BERGSTROM K. An Alkoxylate Mixture and Its Use As A Cleaning Agent for Hard Surfaces. Patent EP 1 519 907 B1, 2003.

    [3] HRECZUCH W. et. al. The Method to Obtain Nonionic Surfactant of A Steric Amphiphilic Structure. Patent PL 215601,2011.

    [4] CHRU?CIEL A; HRECZUCH W; WOJDY?A H; CZAJA K; JANIK J; DOMARECKA D; DOMARECKI W. The Composition and Method to Obtain An Oxyalkylation Catalyst.Patent Application PL 398518, 2012.

    [5] CHRUSCIEL A; HRECZUCH W; JANIK J; CZAJA K;DZIUBEK K; FLISAK Z; SWINAREW A. Characterization of A Double Metal Cyanide (DMC)-Type Catalyst in the Polyoxypropylation Process: Effects of Catalyst Concentration.Ind. Eng. Chem. Res. 2014, 53, 6636-6646.

    [6] SZYMANOWSKI J; ATAMA?CZUK B; SZEWCZYK H.Molar Mass Distribution of the Polyoxyethylene Derivatives of Alcohols, Alkylphenols and Alkylamines. J. Chem. Tech.Biotechnol 1987, 40, 1-10.

    [7] KHORZOUGHI M; BEIRANVAND M; RASAEI M.Investigation of A New Multiphase ow Choke Correlation by Linear and Non-linear Optimization Methods and Monte Carlo Sampling. J. Petrol Explor. Prod. Technol 2013, 3, 279–285.

    [8] NATTA G; MANTICA E. The Distribution of Products in a Series of Consecutive Competitive Reactions., J. Am. Chem. Soc 1952, 74, 3152-3156.

    [9] POSKROBKO H; POSKROBKO J; HRECZUCH W;SOBCZY?SKA A; SZYMANOWSKI J. Oxyethylation and Oxypropylation of Alcohols of Low Relative Molecular Mass in the Presence of Amine-type Catalysts. J. Chem. Technol.Biotechnol 2000, 75, 547-552.

    [10] NM-CAD proprietary software, developed by MEXEO,K?dzierzyn-Ko?le, Poland.

    [11] FLORY P I. Molecular Size Distribution in Ethylene Oxide Polymers. J. Am. Oil Soc. 1940, 62, 1561-1565.

    [12] TATIANA S. ARTURI, NOEMI E. ZARITZKY, and EDGARDO M. CONTRERAS, Simple High-Performance Liquid Chromatography-Ultraviolet Method to Quantify the Molecular Size Distribution of Nonylphenol Ethoxylates., Ind.Eng. Chem. Res., Ind. Eng. Chem. Res 2014, 53, 1327-1333.

    [13] WEIBULL B; TORNQUIST J. Intern. Kongr. für Grenzflachenaktive Stoffe. Zurich 1972.

    [14] RUPP M; RUBACK W; KLEMM E. Alcohol Ethoxylation Kinetics: Proton Transfer Influence on Product Distribution in Microchannels. Chem. Eng. Process 2013, 74, 187-192.

    [15] HRECZUCH W; SZYMANOWSKI J. Synthesis of Surfactants with Narrow Range Distribution of Polyoxyethylene Chain J.Am. Oil. Chem. Soc. 1996, 73:73-78.

    [16] MATERNA K; TIC W J; SOKOLOWSKI A; HRECZUCH W.C12 Hydroxyester Ethoxylates as Nonionic Surfactants. Centr Eur. J. Chem. 2011, 9, 300-307.

    欧美zozozo另类| 免费黄色在线免费观看| 免费黄网站久久成人精品| 联通29元200g的流量卡| 亚洲av日韩在线播放| 成人三级黄色视频| 丰满乱子伦码专区| 久久人人爽人人片av| 国产精品永久免费网站| 精品不卡国产一区二区三区| 亚洲av免费在线观看| 成人无遮挡网站| 一级av片app| 中国美白少妇内射xxxbb| 久久精品91蜜桃| 色综合亚洲欧美另类图片| 亚洲av电影在线观看一区二区三区 | 观看免费一级毛片| 三级男女做爰猛烈吃奶摸视频| 精华霜和精华液先用哪个| 成人漫画全彩无遮挡| 免费黄色在线免费观看| 亚洲,欧美,日韩| 两个人的视频大全免费| 日本猛色少妇xxxxx猛交久久| 国产亚洲精品久久久com| av黄色大香蕉| 国产免费福利视频在线观看| 国产黄a三级三级三级人| 午夜精品在线福利| 99久久精品一区二区三区| 两个人的视频大全免费| 噜噜噜噜噜久久久久久91| 狂野欧美白嫩少妇大欣赏| 青春草视频在线免费观看| av.在线天堂| 国产成人91sexporn| 美女内射精品一级片tv| 久久久a久久爽久久v久久| 熟女电影av网| 久久精品影院6| 国产精品嫩草影院av在线观看| 国产又黄又爽又无遮挡在线| 精品一区二区三区人妻视频| 男女国产视频网站| 国产亚洲一区二区精品| 精品酒店卫生间| av国产免费在线观看| 一个人免费在线观看电影| 国产69精品久久久久777片| 亚洲精华国产精华液的使用体验| 亚洲乱码一区二区免费版| 我要看日韩黄色一级片| 免费看日本二区| 内射极品少妇av片p| 亚洲欧美日韩无卡精品| 男女视频在线观看网站免费| 中文在线观看免费www的网站| 最近中文字幕2019免费版| 99热网站在线观看| 国产精品久久久久久精品电影小说 | 中文欧美无线码| 成人漫画全彩无遮挡| 日本黄色片子视频| 亚洲色图av天堂| 久久精品久久久久久噜噜老黄 | 国产成人精品久久久久久| 国产成人免费观看mmmm| 午夜老司机福利剧场| 午夜久久久久精精品| 插逼视频在线观看| 久久草成人影院| 亚洲欧美日韩东京热| 五月伊人婷婷丁香| 狠狠狠狠99中文字幕| 亚洲av不卡在线观看| 九草在线视频观看| 精品国产三级普通话版| 成人无遮挡网站| 久久久a久久爽久久v久久| 日本黄大片高清| 又粗又硬又长又爽又黄的视频| 在线免费十八禁| 国产精品电影一区二区三区| 日韩欧美精品v在线| 久久精品久久精品一区二区三区| 国产在视频线在精品| 国产av码专区亚洲av| 精品国产三级普通话版| 99久国产av精品国产电影| 亚洲三级黄色毛片| 国产精品福利在线免费观看| 好男人视频免费观看在线| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩一区二区视频在线观看视频在线 | 夜夜爽夜夜爽视频| 神马国产精品三级电影在线观看| 成人午夜高清在线视频| 亚洲在线观看片| 午夜福利网站1000一区二区三区| 免费观看人在逋| 欧美成人精品欧美一级黄| 国产白丝娇喘喷水9色精品| 亚洲怡红院男人天堂| 国模一区二区三区四区视频| 欧美一级a爱片免费观看看| 亚洲精品aⅴ在线观看| 国国产精品蜜臀av免费| 欧美成人精品欧美一级黄| 亚洲成色77777| 99久久中文字幕三级久久日本| 国内精品宾馆在线| 高清视频免费观看一区二区 | 一二三四中文在线观看免费高清| 久久人妻av系列| 国产女主播在线喷水免费视频网站 | 国产黄片视频在线免费观看| 高清在线视频一区二区三区 | 性插视频无遮挡在线免费观看| 国产乱人偷精品视频| 午夜福利在线观看免费完整高清在| 中文精品一卡2卡3卡4更新| 在线免费十八禁| 中国美白少妇内射xxxbb| 少妇被粗大猛烈的视频| 一级毛片我不卡| 黄色欧美视频在线观看| eeuss影院久久| 美女内射精品一级片tv| 十八禁国产超污无遮挡网站| h日本视频在线播放| 寂寞人妻少妇视频99o| 看片在线看免费视频| 亚洲,欧美,日韩| 日韩国内少妇激情av| 午夜福利高清视频| 亚州av有码| 狠狠狠狠99中文字幕| 麻豆成人午夜福利视频| 高清在线视频一区二区三区 | 1000部很黄的大片| 3wmmmm亚洲av在线观看| 久久久久久久午夜电影| 国产精品久久久久久av不卡| 亚洲人与动物交配视频| 国产在视频线在精品| 伊人久久精品亚洲午夜| 又粗又爽又猛毛片免费看| 久久99热这里只有精品18| 免费观看人在逋| 99热全是精品| 国产精品电影一区二区三区| 国产午夜福利久久久久久| 国产亚洲午夜精品一区二区久久 | 国产精品久久久久久精品电影小说 | 免费不卡的大黄色大毛片视频在线观看 | 少妇熟女欧美另类| 九九在线视频观看精品| 我要看日韩黄色一级片| 最近最新中文字幕大全电影3| 国产乱人视频| 美女内射精品一级片tv| 亚洲精华国产精华液的使用体验| 国产精品国产三级国产av玫瑰| 水蜜桃什么品种好| 国内精品一区二区在线观看| 变态另类丝袜制服| 成年av动漫网址| 国产精品一区二区性色av| 青春草亚洲视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产av码专区亚洲av| 狂野欧美激情性xxxx在线观看| 中文字幕免费在线视频6| 日本免费在线观看一区| 国产高清有码在线观看视频| 99久久精品一区二区三区| www日本黄色视频网| 亚洲四区av| 国产伦一二天堂av在线观看| 久久人人爽人人爽人人片va| 九九爱精品视频在线观看| 精品一区二区三区人妻视频| 男女视频在线观看网站免费| 啦啦啦韩国在线观看视频| 亚洲精品日韩av片在线观看| 成人无遮挡网站| 国产成人午夜福利电影在线观看| 亚洲av.av天堂| 在线免费观看不下载黄p国产| 日韩,欧美,国产一区二区三区 | 国产伦理片在线播放av一区| 亚洲综合色惰| 一级av片app| 亚洲第一区二区三区不卡| 嫩草影院精品99| 高清在线视频一区二区三区 | 麻豆成人午夜福利视频| .国产精品久久| 久久久久久大精品| 国产av一区在线观看免费| 午夜福利在线在线| 国内精品宾馆在线| 亚洲经典国产精华液单| 亚洲av一区综合| 久久人妻av系列| 亚洲欧美精品专区久久| 国产精品美女特级片免费视频播放器| 国产高清不卡午夜福利| 校园人妻丝袜中文字幕| 中文字幕av在线有码专区| 久久久久精品久久久久真实原创| 青春草国产在线视频| 最近中文字幕2019免费版| 99九九线精品视频在线观看视频| 国产亚洲av片在线观看秒播厂 | 亚洲av成人精品一二三区| 精品人妻熟女av久视频| 亚洲成色77777| 亚洲精品国产成人久久av| 国产乱人偷精品视频| 2022亚洲国产成人精品| 久久久色成人| 不卡视频在线观看欧美| 九色成人免费人妻av| 久久久久久久国产电影| 久久99精品国语久久久| 国产真实乱freesex| 少妇熟女aⅴ在线视频| 天堂√8在线中文| 欧美极品一区二区三区四区| 亚洲成人精品中文字幕电影| 午夜爱爱视频在线播放| 成年女人看的毛片在线观看| 卡戴珊不雅视频在线播放| 天堂av国产一区二区熟女人妻| 国产不卡一卡二| 午夜爱爱视频在线播放| 欧美三级亚洲精品| 欧美xxxx黑人xx丫x性爽| 亚洲精华国产精华液的使用体验| 一边亲一边摸免费视频| 亚洲天堂国产精品一区在线| 在线观看一区二区三区| 亚洲精品456在线播放app| 日韩欧美 国产精品| 日韩av在线大香蕉| 国产伦在线观看视频一区| 国产成年人精品一区二区| 18禁动态无遮挡网站| 成人毛片a级毛片在线播放| 六月丁香七月| 少妇熟女aⅴ在线视频| 国国产精品蜜臀av免费| 国产探花在线观看一区二区| 爱豆传媒免费全集在线观看| 在线播放国产精品三级| 3wmmmm亚洲av在线观看| 国产精品伦人一区二区| 亚洲精品乱码久久久久久按摩| 能在线免费看毛片的网站| 国产伦一二天堂av在线观看| 久久久精品欧美日韩精品| 精品久久久噜噜| 26uuu在线亚洲综合色| 热99re8久久精品国产| 中文亚洲av片在线观看爽| 一夜夜www| 小蜜桃在线观看免费完整版高清| 国产av不卡久久| 国产在线一区二区三区精 | 51国产日韩欧美| 熟女电影av网| 性色avwww在线观看| 日韩,欧美,国产一区二区三区 | 一个人看的www免费观看视频| 美女xxoo啪啪120秒动态图| 啦啦啦观看免费观看视频高清| 久久这里只有精品中国| 亚洲人成网站在线播| av在线观看视频网站免费| 国产淫片久久久久久久久| 国产精品一区二区在线观看99 | 亚洲va在线va天堂va国产| 国内精品一区二区在线观看| 激情 狠狠 欧美| 成人综合一区亚洲| av又黄又爽大尺度在线免费看 | 日本色播在线视频| 干丝袜人妻中文字幕| 91久久精品电影网| 亚洲婷婷狠狠爱综合网| 看十八女毛片水多多多| 国产黄片视频在线免费观看| 久久精品国产亚洲av涩爱| or卡值多少钱| 久久久久久久久中文| 国产精品久久久久久精品电影| 18禁动态无遮挡网站| 亚洲精品,欧美精品| 乱系列少妇在线播放| 精品熟女少妇av免费看| 禁无遮挡网站| 国产一级毛片七仙女欲春2| 97热精品久久久久久| 日韩国内少妇激情av| 我要看日韩黄色一级片| 国产精品不卡视频一区二区| 丰满乱子伦码专区| 亚洲人成网站在线观看播放| 久久久久久国产a免费观看| 最近最新中文字幕大全电影3| 亚洲美女视频黄频| 日本三级黄在线观看| 少妇被粗大猛烈的视频| 成人亚洲精品av一区二区| 亚洲av中文av极速乱| 简卡轻食公司| 日韩人妻高清精品专区| 最近2019中文字幕mv第一页| 中文字幕av成人在线电影| 一本一本综合久久| 午夜福利网站1000一区二区三区| 在线免费十八禁| 国产极品精品免费视频能看的| 亚洲婷婷狠狠爱综合网| 一级黄色大片毛片| 能在线免费看毛片的网站| 人人妻人人澡人人爽人人夜夜 | 在现免费观看毛片| 国产成人精品久久久久久| 一本久久精品| 午夜福利在线观看免费完整高清在| 97热精品久久久久久| 亚洲经典国产精华液单| 91精品伊人久久大香线蕉| 国内精品一区二区在线观看| 亚洲国产精品成人综合色| 久久久久久久久久成人| 91在线精品国自产拍蜜月| 亚洲无线观看免费| 偷拍熟女少妇极品色| 91精品一卡2卡3卡4卡| 国国产精品蜜臀av免费| 97超碰精品成人国产| 亚洲欧洲日产国产| 联通29元200g的流量卡| 国产精品日韩av在线免费观看| 夜夜爽夜夜爽视频| 亚洲国产欧美人成| 国产在线男女| 听说在线观看完整版免费高清| 秋霞伦理黄片| or卡值多少钱| 国产免费男女视频| 日韩av在线大香蕉| 小说图片视频综合网站| 午夜福利在线观看吧| 免费av不卡在线播放| 18禁在线无遮挡免费观看视频| 大香蕉久久网| 国产免费视频播放在线视频 | 精品久久久久久电影网 | 亚洲欧洲日产国产| kizo精华| 日本av手机在线免费观看| 亚洲性久久影院| 超碰av人人做人人爽久久| 麻豆成人av视频| 少妇裸体淫交视频免费看高清| 亚洲国产日韩欧美精品在线观看| 国产精品精品国产色婷婷| 中国美白少妇内射xxxbb| 我要看日韩黄色一级片| 国产探花极品一区二区| 国产高清有码在线观看视频| 天堂影院成人在线观看| 日韩三级伦理在线观看| 久久这里有精品视频免费| 国产一区二区在线av高清观看| 联通29元200g的流量卡| 亚洲伊人久久精品综合 | 亚洲天堂国产精品一区在线| 国产精品电影一区二区三区| 最新中文字幕久久久久| 韩国av在线不卡| 国产老妇女一区| 国产黄色小视频在线观看| 欧美潮喷喷水| 国产精品,欧美在线| 国产精品伦人一区二区| 一区二区三区高清视频在线| 日本一二三区视频观看| 黄色配什么色好看| 国产一区二区在线av高清观看| 中文在线观看免费www的网站| 国产精品无大码| 免费在线观看成人毛片| 两个人的视频大全免费| 性色avwww在线观看| 久久鲁丝午夜福利片| 欧美精品国产亚洲| 欧美成人精品欧美一级黄| 国内精品美女久久久久久| 国产免费视频播放在线视频 | 亚洲精品一区蜜桃| 日本五十路高清| 精品欧美国产一区二区三| 九草在线视频观看| 国产精品爽爽va在线观看网站| 99热这里只有是精品50| 国产一区二区在线观看日韩| 麻豆精品久久久久久蜜桃| 国内精品美女久久久久久| 天天躁日日操中文字幕| 欧美日韩国产亚洲二区| 一级毛片电影观看 | 亚洲精品aⅴ在线观看| av国产久精品久网站免费入址| 国产伦理片在线播放av一区| 少妇的逼好多水| 十八禁国产超污无遮挡网站| 男女下面进入的视频免费午夜| 久久精品久久久久久久性| 男女啪啪激烈高潮av片| 日韩亚洲欧美综合| 欧美最新免费一区二区三区| 少妇的逼好多水| 亚洲国产精品久久男人天堂| 国产真实伦视频高清在线观看| 春色校园在线视频观看| 国产精品一区二区在线观看99 | 日韩制服骚丝袜av| 亚洲精品国产av成人精品| 黄色日韩在线| 亚洲欧美日韩卡通动漫| 搞女人的毛片| 中国国产av一级| 成人午夜精彩视频在线观看| 精品午夜福利在线看| 国产乱人偷精品视频| 国产精品av视频在线免费观看| 久久人人爽人人片av| 中文字幕av成人在线电影| 亚洲欧美中文字幕日韩二区| 国产一级毛片七仙女欲春2| 国产精品久久视频播放| 国产成人一区二区在线| 99久久精品热视频| 你懂的网址亚洲精品在线观看 | 国产免费男女视频| 91久久精品国产一区二区三区| 久久久久久久亚洲中文字幕| 超碰97精品在线观看| 亚洲一区高清亚洲精品| av在线蜜桃| 日韩强制内射视频| 一级爰片在线观看| 少妇人妻一区二区三区视频| 亚洲国产精品成人久久小说| 国语自产精品视频在线第100页| av在线播放精品| 看十八女毛片水多多多| 卡戴珊不雅视频在线播放| 亚洲国产精品合色在线| 超碰av人人做人人爽久久| 日韩欧美在线乱码| 国产亚洲91精品色在线| 中文字幕人妻熟人妻熟丝袜美| 99热网站在线观看| 国产精品麻豆人妻色哟哟久久 | 青青草视频在线视频观看| 看黄色毛片网站| 国产一区二区亚洲精品在线观看| 美女被艹到高潮喷水动态| 亚洲av福利一区| 99热这里只有是精品50| 亚洲av成人av| 亚洲三级黄色毛片| 色哟哟·www| 国产精品不卡视频一区二区| 国产亚洲最大av| 国产在线一区二区三区精 | 精品一区二区三区视频在线| 你懂的网址亚洲精品在线观看 | 久久精品人妻少妇| 亚洲成人久久爱视频| 亚洲精品日韩在线中文字幕| 欧美色视频一区免费| 午夜福利在线观看免费完整高清在| 国产亚洲精品久久久com| 永久免费av网站大全| 午夜爱爱视频在线播放| 波多野结衣高清无吗| 免费看av在线观看网站| 免费看a级黄色片| 亚洲人与动物交配视频| 日本五十路高清| 免费看av在线观看网站| 久久精品国产自在天天线| 超碰97精品在线观看| 免费观看a级毛片全部| 久久久久久久国产电影| 久久99热这里只有精品18| 又黄又爽又刺激的免费视频.| 亚洲熟妇中文字幕五十中出| 99热这里只有是精品50| 亚洲一级一片aⅴ在线观看| 99热精品在线国产| 国产一区二区在线av高清观看| 秋霞伦理黄片| 国产高清国产精品国产三级 | 97人妻精品一区二区三区麻豆| 韩国高清视频一区二区三区| 亚洲欧美日韩无卡精品| 久久精品国产亚洲av天美| 中文字幕av成人在线电影| 99热6这里只有精品| 中文欧美无线码| 免费黄网站久久成人精品| 亚洲av男天堂| 国产综合懂色| 国产精品,欧美在线| 日韩视频在线欧美| 欧美性猛交黑人性爽| 日韩av在线免费看完整版不卡| 精品欧美国产一区二区三| 麻豆久久精品国产亚洲av| 亚洲国产色片| 亚洲真实伦在线观看| 在线免费十八禁| 在线观看av片永久免费下载| 久久99热6这里只有精品| 男人舔女人下体高潮全视频| 国产综合懂色| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 不卡视频在线观看欧美| 性插视频无遮挡在线免费观看| 亚洲欧美精品自产自拍| 亚洲精品国产av成人精品| 高清午夜精品一区二区三区| 男人和女人高潮做爰伦理| 日本免费一区二区三区高清不卡| 国产精品国产三级国产专区5o | 久久精品国产99精品国产亚洲性色| 国产三级在线视频| 亚洲精品乱码久久久v下载方式| 色吧在线观看| 日本一本二区三区精品| АⅤ资源中文在线天堂| 国产精品电影一区二区三区| 亚洲av成人av| 国产探花在线观看一区二区| 又粗又硬又长又爽又黄的视频| 成人欧美大片| 国产精品精品国产色婷婷| 内射极品少妇av片p| 色播亚洲综合网| 欧美最新免费一区二区三区| 亚洲va在线va天堂va国产| 国产精品综合久久久久久久免费| 日韩 亚洲 欧美在线| 国产av在哪里看| 国产v大片淫在线免费观看| 国产精品av视频在线免费观看| 亚洲av免费高清在线观看| 级片在线观看| 亚洲av日韩在线播放| 免费观看的影片在线观看| 久久精品人妻少妇| 干丝袜人妻中文字幕| 日本wwww免费看| 在线播放国产精品三级| 国产精品,欧美在线| 看片在线看免费视频| 国产一区二区在线av高清观看| 看黄色毛片网站| 成人漫画全彩无遮挡| 国产精品99久久久久久久久| 国产成人a∨麻豆精品| 欧美日韩一区二区视频在线观看视频在线 | 毛片一级片免费看久久久久| 在线a可以看的网站| 国产淫片久久久久久久久| 2021天堂中文幕一二区在线观| 午夜福利成人在线免费观看| 国产精品国产三级国产av玫瑰| 国产精品野战在线观看| 偷拍熟女少妇极品色| 看十八女毛片水多多多| 亚洲最大成人手机在线| 久久欧美精品欧美久久欧美| 中文字幕亚洲精品专区| 亚洲欧美成人综合另类久久久 | 国产精品一区二区三区四区久久| 看黄色毛片网站| 女人十人毛片免费观看3o分钟| 中文亚洲av片在线观看爽| 国产美女午夜福利| 97人妻精品一区二区三区麻豆| 日韩欧美国产在线观看| 国国产精品蜜臀av免费| 看免费成人av毛片| 午夜免费激情av| 2021天堂中文幕一二区在线观| 国产av一区在线观看免费| 国产黄a三级三级三级人| 高清视频免费观看一区二区 | 亚洲三级黄色毛片| 欧美xxxx黑人xx丫x性爽| 亚洲国产精品sss在线观看| a级一级毛片免费在线观看| 春色校园在线视频观看| 久久精品久久久久久噜噜老黄 | 一级毛片我不卡|