反渗透膜及水溶液内扩散过程的分子模拟研究
|
摘要
本论文针对目前传统试验手段用在研究分子扩散过程以及扩散系数方面的局限性,通过分子动力学模拟(Molecular Dynamic Simulation)的手段,深入详细地研究了小分子在反渗透复合膜以及水溶液内的扩散过程。主要进行了四方面的研究工作:分子在体系内扩散系数的分子模拟方法研究;甘氨酸等12种氨基酸分子在水溶液内扩散系数的模拟;氧气等9种气体分子在水溶液内扩散系数的模拟;水和盐分子(NaCl, MgCl2,CaSO4,K2SO4)在8种反渗透复合膜中扩散状态及扩散系数的研究,膜的结构单体与水和盐分子在膜内扩散系数关系的分析与讨论。
1)论文首次从模型的分子数样本以及模拟步长出发,对影响分子在体系内扩散系数的主要因素进行了研究。并针对动力学模拟数据提出了新的分析方法。
①对浓度分别为0.7mol/L和0.9mol/L的甘氨酸水体系进行了动力学模拟。建立模型时分别加入了不同数目的分子:对0.7mol/L的甘氨酸模型的构建,分别在立方晶格中加入了3、5、7、10、12个甘氨酸分子;对于0.9mol/L的甘氨酸则分别在立方晶格中加入了4、7、9、13、15个甘氨酸分子。扩散系数模拟值与实验值的对比表明分子数样本对于扩散系数的模拟结果有较大的影响。分子数越多,则模拟结果与实验值越接近,当分子数达到一定值时,扩散系数趋于一个稳定的值。
论文还针对不同温度下氧气在水中的扩散过程作了关于样本数的模拟。分别考察了温度为303K和313K条件下的气体扩散。模拟体系中氧气分子数分别采用1,2,3个。结果表明随着模型中分子数的增多,得到的扩散系数的值也逐渐增大,并逐渐恒定于正确值。
②通过对比传统分析方法以及微分-区限变分法用在均方位移(MSD)曲线上的优缺点,本论文提出了新的数据分析方法-聚类分析方法。该方法克服了传统方法用在MSD线性不好时给分析结果带来误差的缺点,对微分-区限变分法进行了改进。将聚类方法和微分-区限变分法应用在氧气在水中的扩散系数的分析中,实验值,微分-区限变分法值和聚类方法值分别为2.10×10-9m2s-1,1.80×10-9m2s-1,2.00×10-9m2s-1,聚类分析方法的误差要小于微分-区限变分法,充分证明聚类方法是一个有效的分析方法。
③模拟甲烷和氨气在水中的扩散过程时,对0.1fs的步长值延长了2-3倍的模拟时间,结果表明模拟时间的延长对最后的模拟结果没有明显的影响,从而排除了模拟时间对扩散系数模拟结果的影响。另外,论文对甲烷,氨气和二氧化碳气体在水中的扩散过程采用了0.1fs-1.5fs范围内几个不同的步长值,模拟结果表明,步长值对扩散系数的模拟结果有较大影响。随着步长的增加,扩散系数有一个极值,与该极值相邻的两个步长对应的扩散系数与实验值最为接近。这两个步长值为极值对应的步长值的70%和150%。结合步长在模拟计算中的作用,本论文认为对于一个特定体系,应在0.1fs到1.5fs之间取多个值进行模拟,取模拟值最大的步长的70%即为模拟的合理值。
2)本论文首次通过分子模拟的手段研究了氨基酸在水溶液内的扩散过程,并利用微分-区限变分法计算了氨基酸在水中的扩散系数。目前氨基酸在水中扩散系数主要通过实验手段来获得,对于氨基酸在水中的扩散系数的模拟,目前尚未见报道。本文通过分子动力学模拟的手段模拟计算了甘氨酸、丙氨酸等12种氨基酸在水中的扩散系数。模拟值与实验值的误差小于7%,证明了分子模拟方法应用在扩散系数计算方面的有效性。
3)本论文通过分子模拟的手段模拟了氧气等9种气体在水溶液内扩散过程,并利用聚类分析方法计算了气体在水中的扩散系数。对于O2,CO2,CH4,NH3四种气体在水中的扩散系数的模拟,已有文献的模拟值与实验值的误差最大超过50%,而本论文的模拟数据与实验值的误差小于30%。充分证明了本文所采用的分子模拟途径的可靠性。
4)本文以8种反渗透复合膜为主要研究对象,模拟了海水中主要成分(H2O,NaCl, MgCl2,CaSO4和K2SO4)在膜内的扩散过程并计算了各粒子在膜内的扩散系数。
①模拟了TMC/MPD膜在体系中的扩散系数,模拟结果为0.0115A2/ps,与已有的对于该膜的模拟结果0.015 A2/ps基本吻合。证明了本文模拟过程中对于模型建立,模拟方法以及数据分析方法的有效性。
②水分子在TMC/MPD膜内的扩散状态表明水分子与高分子链上的功能基团有明显的分子间相互作用。水分子在靠近苯环的区域内与苯环有排斥运动;在靠近酰胺基的区域内有吸引运动。
③在所有膜中,水分子在膜HT/MMPD中的扩散系数最大,远远高于在其他膜内的扩散系数。对比水分子在膜TMC/MMPD和TMC-HT/MMPD的扩散系数发现水分子在加入了HT的膜中扩散系数要高一些了;对比水分子在膜TMC/MPD和膜TMC/MMPD内的扩散系数,水分子在膜TMC/MMPD的扩散系数要小很多。由于MMPD单体在苯环上含有一个甲基,在同样数量和比例生成聚合物膜的情况下,MMPD比MPD的刚性要强一些,从而不利于水分子的通通过。对比水分子在膜CFIC/MMPD和膜CFIC/MPD的扩散系数也可得出同样的结论。
④NaCl分子在膜HT/MMPD中的扩散系数最大,而在膜ICIC-HT/MMPD中的扩散系数则最小,其次分别是膜ICIC-IPC/MPD ,膜TMC-HT/MMPD,膜CFIC/MPD以及膜TMC/MMPD。模拟结果中只有ICIC-IPC/MPD与实验结果有较大差距。
⑤在所模拟的8种膜内,Na+和Cl-在膜内的扩散系数值相差较大,NaCl分子的运动受限于运动较慢的离子。膜的种类不同,制约NaCl在膜内扩散的离子便不同:对于膜TMC/MPD以及膜CFIC/MMPD,NaCl分子的扩散受Cl-的制约;而其他6种类型的膜则受Na+的制约。
⑥对于膜ICIC-HT/MPD,水分子在MgCl2,CaSO4和K2SO4体系中的扩散系数与在NaCl体系中的结果基本一致,该结果表明在同种膜中,水分子的扩散过程不受体系中盐分子类型的影响。
⑦在ICIC-HT/MPD反渗透膜中,阴离子(Cl-,SO42-)的扩散系数均远大于阳离子(Mg2+,Ca2+,K+)的扩散系数。这与NaCl水体系中,Cl-的扩散系数远大于Na+的扩散系数的运动状况一致,该结果说明当盐分子单独存在时,其在膜内的扩散过程受哪个离子制约只与膜种类有关,与盐分子本身无关。
This paper used molecular Dynamic Simulation method to detailedly study the diffusion process and diffusion coefficient of small molecules in reverse osmosis membrane and water solution, and got some inaccessible results with traditional experimental method. There are mainly four aspects discussed in the paper: (1) Research of molecular Dynamic Simulation method for molecular diffusion coefficients. (2) The simulation of diffusion coefficients of Glycine and other 11 kinds of amino acids in water solution. (3) The simulation of diffusion coefficients of oxygen and other 8 kinds of gas molecules in water solution. (4) The diffusion process and diffusion coefficients of H2O, NaCl, MgCl2, CaSO4, K2SO4 in 8 kinds of reverse osmosis membrane, and relationships of the membranous structure monomer and the diffusion coefficients.
1) This paper firstly researched the factors which affect diffusion coefficients in the system from the molecular number sample and simulation time step of the model, and advanced a new analytic method with the Dynamic Simulation data.
①This paper simulated diffusion process of Glycine solutions of 0.7mol/L and 0.9mol/L. In the Dynamic Simulation, the models of 0.7mol/L Glycine solution were constructed by adding 3, 5, 7, 10, 12 Glycine molecules into the cell, and 4,7,9,13,15 Glycine molecules were added into the cell of 0.9mol/L Glycine solution models. Comparing the simulation values with the experimental results, the sample number of molecules had some effect on the simulation of diffusion coefficients, which was closer to the experimental results with more molecules added into the cell, and approached a steady value when the molecule number exceeds a line.
This paper also simulated the diffusion process of oxygen in water of different temperature, such as 303K and 313K. The models were separately constructed by adding 1,2,3 oxygen molecules into the cell . Simulation values were also found approaching experimental results with more molecules in the cell.
②Through comparing the use of traditional analysis method in MSD curves with differentiation interval variation method, this paper advanced a new data analysis method-clustering method, which overcome shortcomings of traditional method used in some MSD curves with bad linear direction and modifies differentiation interval variation method. Using clustering method and differentiation interval variation method separately in the analysis of diffusion coefficient of oxygen in the water, the results of them and the experimental value were separately 2.00×10-9m2s-1,1.80×10-9m2s-1,2.10×10-9m2s-1, which obviously proved clustering method effective.
③This paper increased simulation time of 0.1fs time step to 2~3 times longer than before in the simulation of diffusion process in water of methane and ammonia, the results showed the increasing time of simulation did not have obvious effect. Then several different time step between 0.1fs to 0.5 fs were used in the simulation of diffusion process in water of methane, ammonia and carbon dioxide, results showed time step was the effective factor to the diffusion coefficient, which rose to a maximum with the increasing of time step. The diffusion coefficients corresponding to the two time steps close to the maximum were most close to the experimental value. This two time steps were 70 and 150 percent of the time step corresponding to the maximum diffusion coefficient value. Considering the effect of time step in the simulation, this paper advised several time steps between 0.1fs to 0.5 fs should be taken in the simulation of a special system, and took 70 percent of the time step corresponding to maximum diffusion coefficient value as the reasonable value.
2) The paper dealt with for the first time the diffusion process of amino acid in the water solution by the way of molecular simulation, at the same time, the diffusion coefficient of amino acid in water was calculated by the method of differentiation interval variation. At present, the method of gaining the diffusion coefficient of amino acid in the water mainly based on the experiments, on the contrary, the way of simulation for the diffusion coefficient of amino acid in the water wasn’t published preciously. In this thesis, the diffusion coefficient of amino acid, such as Glycine and other 11 kinds of amino acid, was calculated using the method of molecular dynamic simulation. The error between the simulated results and the experimental results was less than seven percent, and therefore the way of calculating diffusion coefficient by molecular dynamic simulation was proved to be effective.
3) In this paper, the diffusion process of oxygen and other 8 kinds of gases in the water solution was simulated by the way of molecular simulation, meanwhile, the diffusion coefficient of gas in the water was calculated by the method of clustering method. The work present in this thesis about the simulation of getting diffusion coefficient of the oxygen, carbon dioxide, methanol and ammonia in the water solution was enough to prove the reliability of the method of molecular simulation, in which the error between the simulated data in this paper and the experimental results was less than 30 percent, however, the biggest error between simulated data in the literature and experimental result was more than 50 percent.
4) In this paper, eight kinds of reverse osmosis compound membrane were used as the main base case, in which the diffusion process of the main components (water, sodium chloride, magnesium chloride, calcium sulphate and kalium sulphate ) in the seawater within the membrane were simulated and the diffusion coefficient of every particle in the membrane is calculated.
①The diffusion coefficient of TMC/MPD membrane in the system was simulated with the result of 0.0115A2/ps, which was almost inosculated with the existing simulated result 0.015 A2/ps for this membrane. In a word, the foundation of the model, the method of simulation and the analytic method for the data were effective in the simulation process for this paper.
②The diffusion state of water molecules in the TMC/MPD indicates that there were obvious interaction between water molecules and Functional groups of the polymer. When the Water molecules were in the field nearing benzene fragments there was a significant repulsion action. On the other way, when water molecules near the amido-group there was an attracting tend movement.
③The diffusion coefficient of water molecule in the membranes of HT/MMPD was the biggest in all of the membranes we studied. It will be founded that the diffusion coefficient of water molecule in the membrane of TMC-HT/MMPD was higher somewhat than TMC/MMPD which was added TH; Analysis of the water molecule diffusion coefficient in the membrane of TMC/MPD and TMC/MMPD displayed that the water molecule diffusion coefficient in the membrane of TMC/MPD was much bigger than in TMC/MMPD. For the monomer of the membrane MMPD had a methyl substitution in the benzene ring, the rigidity of the membrane MMPD was stronger than MPD when the membranes consist of the same account and scale of the monomer. As a result, this did not benefit the molecules of water transitting through the membrane MMPD. The same conclusion was obtained by comparing the diffusion coefficient of the molecules of water in the membrane CFIC/MMPD and CFIC/MPD.
④NaCl molecule had the biggest diffusion coefficient in the membrane of HT/MMPD. Then it took turns to the membrane of ICIC-IPC/MPD, TMC-HT/MMPD, CFIC/MPD , TMC/MMPD, and the least in the membrane of ICIC-HT/MMPD. Only ICIC-IPC/MPD had the bigger difference between simulation result and the experiment result.
⑤In eight kinds of membranes we simulated, the diffusion coefficient value of Na+ and Cl- had much difference. The movement of NaCl molecule was restricted by the ions which moves slowly. The kind of the ions which constrained the molecule NaCl was variety in different membranes: in the membranes of TMC/MPD and CFIC/MMPD, the constrain ion was Cl-, while was Na+ in other six membranes.
⑥For membrane ICIC-HT/MPD ,the diffusion coefficients of water in the systems containing MgCl2,CaSO4 and K2SO4 were consistent with that in the NaCl system, this result demonstrated that movement of molecule of water was not affected by the kinds of salt molecules in the same membrane.
⑦The diffusion coefficients of anions (Cl-,SO42-) in the reverse osmosis membrane ICIC-HT/MPD were much bigger than cations (Mg2+,Ca2+,K+). This result was consistent with the diffusion coefficient of the anions Cl- and cations Na+ in the NaCl-H2O system, which indicated that diffusion process of salt molecule in the membrane was only restricted by the membrane species,without relationship with salt molecules itself.
1)论文首次从模型的分子数样本以及模拟步长出发,对影响分子在体系内扩散系数的主要因素进行了研究。并针对动力学模拟数据提出了新的分析方法。
①对浓度分别为0.7mol/L和0.9mol/L的甘氨酸水体系进行了动力学模拟。建立模型时分别加入了不同数目的分子:对0.7mol/L的甘氨酸模型的构建,分别在立方晶格中加入了3、5、7、10、12个甘氨酸分子;对于0.9mol/L的甘氨酸则分别在立方晶格中加入了4、7、9、13、15个甘氨酸分子。扩散系数模拟值与实验值的对比表明分子数样本对于扩散系数的模拟结果有较大的影响。分子数越多,则模拟结果与实验值越接近,当分子数达到一定值时,扩散系数趋于一个稳定的值。
论文还针对不同温度下氧气在水中的扩散过程作了关于样本数的模拟。分别考察了温度为303K和313K条件下的气体扩散。模拟体系中氧气分子数分别采用1,2,3个。结果表明随着模型中分子数的增多,得到的扩散系数的值也逐渐增大,并逐渐恒定于正确值。
②通过对比传统分析方法以及微分-区限变分法用在均方位移(MSD)曲线上的优缺点,本论文提出了新的数据分析方法-聚类分析方法。该方法克服了传统方法用在MSD线性不好时给分析结果带来误差的缺点,对微分-区限变分法进行了改进。将聚类方法和微分-区限变分法应用在氧气在水中的扩散系数的分析中,实验值,微分-区限变分法值和聚类方法值分别为2.10×10-9m2s-1,1.80×10-9m2s-1,2.00×10-9m2s-1,聚类分析方法的误差要小于微分-区限变分法,充分证明聚类方法是一个有效的分析方法。
③模拟甲烷和氨气在水中的扩散过程时,对0.1fs的步长值延长了2-3倍的模拟时间,结果表明模拟时间的延长对最后的模拟结果没有明显的影响,从而排除了模拟时间对扩散系数模拟结果的影响。另外,论文对甲烷,氨气和二氧化碳气体在水中的扩散过程采用了0.1fs-1.5fs范围内几个不同的步长值,模拟结果表明,步长值对扩散系数的模拟结果有较大影响。随着步长的增加,扩散系数有一个极值,与该极值相邻的两个步长对应的扩散系数与实验值最为接近。这两个步长值为极值对应的步长值的70%和150%。结合步长在模拟计算中的作用,本论文认为对于一个特定体系,应在0.1fs到1.5fs之间取多个值进行模拟,取模拟值最大的步长的70%即为模拟的合理值。
2)本论文首次通过分子模拟的手段研究了氨基酸在水溶液内的扩散过程,并利用微分-区限变分法计算了氨基酸在水中的扩散系数。目前氨基酸在水中扩散系数主要通过实验手段来获得,对于氨基酸在水中的扩散系数的模拟,目前尚未见报道。本文通过分子动力学模拟的手段模拟计算了甘氨酸、丙氨酸等12种氨基酸在水中的扩散系数。模拟值与实验值的误差小于7%,证明了分子模拟方法应用在扩散系数计算方面的有效性。
3)本论文通过分子模拟的手段模拟了氧气等9种气体在水溶液内扩散过程,并利用聚类分析方法计算了气体在水中的扩散系数。对于O2,CO2,CH4,NH3四种气体在水中的扩散系数的模拟,已有文献的模拟值与实验值的误差最大超过50%,而本论文的模拟数据与实验值的误差小于30%。充分证明了本文所采用的分子模拟途径的可靠性。
4)本文以8种反渗透复合膜为主要研究对象,模拟了海水中主要成分(H2O,NaCl, MgCl2,CaSO4和K2SO4)在膜内的扩散过程并计算了各粒子在膜内的扩散系数。
①模拟了TMC/MPD膜在体系中的扩散系数,模拟结果为0.0115A2/ps,与已有的对于该膜的模拟结果0.015 A2/ps基本吻合。证明了本文模拟过程中对于模型建立,模拟方法以及数据分析方法的有效性。
②水分子在TMC/MPD膜内的扩散状态表明水分子与高分子链上的功能基团有明显的分子间相互作用。水分子在靠近苯环的区域内与苯环有排斥运动;在靠近酰胺基的区域内有吸引运动。
③在所有膜中,水分子在膜HT/MMPD中的扩散系数最大,远远高于在其他膜内的扩散系数。对比水分子在膜TMC/MMPD和TMC-HT/MMPD的扩散系数发现水分子在加入了HT的膜中扩散系数要高一些了;对比水分子在膜TMC/MPD和膜TMC/MMPD内的扩散系数,水分子在膜TMC/MMPD的扩散系数要小很多。由于MMPD单体在苯环上含有一个甲基,在同样数量和比例生成聚合物膜的情况下,MMPD比MPD的刚性要强一些,从而不利于水分子的通通过。对比水分子在膜CFIC/MMPD和膜CFIC/MPD的扩散系数也可得出同样的结论。
④NaCl分子在膜HT/MMPD中的扩散系数最大,而在膜ICIC-HT/MMPD中的扩散系数则最小,其次分别是膜ICIC-IPC/MPD ,膜TMC-HT/MMPD,膜CFIC/MPD以及膜TMC/MMPD。模拟结果中只有ICIC-IPC/MPD与实验结果有较大差距。
⑤在所模拟的8种膜内,Na+和Cl-在膜内的扩散系数值相差较大,NaCl分子的运动受限于运动较慢的离子。膜的种类不同,制约NaCl在膜内扩散的离子便不同:对于膜TMC/MPD以及膜CFIC/MMPD,NaCl分子的扩散受Cl-的制约;而其他6种类型的膜则受Na+的制约。
⑥对于膜ICIC-HT/MPD,水分子在MgCl2,CaSO4和K2SO4体系中的扩散系数与在NaCl体系中的结果基本一致,该结果表明在同种膜中,水分子的扩散过程不受体系中盐分子类型的影响。
⑦在ICIC-HT/MPD反渗透膜中,阴离子(Cl-,SO42-)的扩散系数均远大于阳离子(Mg2+,Ca2+,K+)的扩散系数。这与NaCl水体系中,Cl-的扩散系数远大于Na+的扩散系数的运动状况一致,该结果说明当盐分子单独存在时,其在膜内的扩散过程受哪个离子制约只与膜种类有关,与盐分子本身无关。
This paper used molecular Dynamic Simulation method to detailedly study the diffusion process and diffusion coefficient of small molecules in reverse osmosis membrane and water solution, and got some inaccessible results with traditional experimental method. There are mainly four aspects discussed in the paper: (1) Research of molecular Dynamic Simulation method for molecular diffusion coefficients. (2) The simulation of diffusion coefficients of Glycine and other 11 kinds of amino acids in water solution. (3) The simulation of diffusion coefficients of oxygen and other 8 kinds of gas molecules in water solution. (4) The diffusion process and diffusion coefficients of H2O, NaCl, MgCl2, CaSO4, K2SO4 in 8 kinds of reverse osmosis membrane, and relationships of the membranous structure monomer and the diffusion coefficients.
1) This paper firstly researched the factors which affect diffusion coefficients in the system from the molecular number sample and simulation time step of the model, and advanced a new analytic method with the Dynamic Simulation data.
①This paper simulated diffusion process of Glycine solutions of 0.7mol/L and 0.9mol/L. In the Dynamic Simulation, the models of 0.7mol/L Glycine solution were constructed by adding 3, 5, 7, 10, 12 Glycine molecules into the cell, and 4,7,9,13,15 Glycine molecules were added into the cell of 0.9mol/L Glycine solution models. Comparing the simulation values with the experimental results, the sample number of molecules had some effect on the simulation of diffusion coefficients, which was closer to the experimental results with more molecules added into the cell, and approached a steady value when the molecule number exceeds a line.
This paper also simulated the diffusion process of oxygen in water of different temperature, such as 303K and 313K. The models were separately constructed by adding 1,2,3 oxygen molecules into the cell . Simulation values were also found approaching experimental results with more molecules in the cell.
②Through comparing the use of traditional analysis method in MSD curves with differentiation interval variation method, this paper advanced a new data analysis method-clustering method, which overcome shortcomings of traditional method used in some MSD curves with bad linear direction and modifies differentiation interval variation method. Using clustering method and differentiation interval variation method separately in the analysis of diffusion coefficient of oxygen in the water, the results of them and the experimental value were separately 2.00×10-9m2s-1,1.80×10-9m2s-1,2.10×10-9m2s-1, which obviously proved clustering method effective.
③This paper increased simulation time of 0.1fs time step to 2~3 times longer than before in the simulation of diffusion process in water of methane and ammonia, the results showed the increasing time of simulation did not have obvious effect. Then several different time step between 0.1fs to 0.5 fs were used in the simulation of diffusion process in water of methane, ammonia and carbon dioxide, results showed time step was the effective factor to the diffusion coefficient, which rose to a maximum with the increasing of time step. The diffusion coefficients corresponding to the two time steps close to the maximum were most close to the experimental value. This two time steps were 70 and 150 percent of the time step corresponding to the maximum diffusion coefficient value. Considering the effect of time step in the simulation, this paper advised several time steps between 0.1fs to 0.5 fs should be taken in the simulation of a special system, and took 70 percent of the time step corresponding to maximum diffusion coefficient value as the reasonable value.
2) The paper dealt with for the first time the diffusion process of amino acid in the water solution by the way of molecular simulation, at the same time, the diffusion coefficient of amino acid in water was calculated by the method of differentiation interval variation. At present, the method of gaining the diffusion coefficient of amino acid in the water mainly based on the experiments, on the contrary, the way of simulation for the diffusion coefficient of amino acid in the water wasn’t published preciously. In this thesis, the diffusion coefficient of amino acid, such as Glycine and other 11 kinds of amino acid, was calculated using the method of molecular dynamic simulation. The error between the simulated results and the experimental results was less than seven percent, and therefore the way of calculating diffusion coefficient by molecular dynamic simulation was proved to be effective.
3) In this paper, the diffusion process of oxygen and other 8 kinds of gases in the water solution was simulated by the way of molecular simulation, meanwhile, the diffusion coefficient of gas in the water was calculated by the method of clustering method. The work present in this thesis about the simulation of getting diffusion coefficient of the oxygen, carbon dioxide, methanol and ammonia in the water solution was enough to prove the reliability of the method of molecular simulation, in which the error between the simulated data in this paper and the experimental results was less than 30 percent, however, the biggest error between simulated data in the literature and experimental result was more than 50 percent.
4) In this paper, eight kinds of reverse osmosis compound membrane were used as the main base case, in which the diffusion process of the main components (water, sodium chloride, magnesium chloride, calcium sulphate and kalium sulphate ) in the seawater within the membrane were simulated and the diffusion coefficient of every particle in the membrane is calculated.
①The diffusion coefficient of TMC/MPD membrane in the system was simulated with the result of 0.0115A2/ps, which was almost inosculated with the existing simulated result 0.015 A2/ps for this membrane. In a word, the foundation of the model, the method of simulation and the analytic method for the data were effective in the simulation process for this paper.
②The diffusion state of water molecules in the TMC/MPD indicates that there were obvious interaction between water molecules and Functional groups of the polymer. When the Water molecules were in the field nearing benzene fragments there was a significant repulsion action. On the other way, when water molecules near the amido-group there was an attracting tend movement.
③The diffusion coefficient of water molecule in the membranes of HT/MMPD was the biggest in all of the membranes we studied. It will be founded that the diffusion coefficient of water molecule in the membrane of TMC-HT/MMPD was higher somewhat than TMC/MMPD which was added TH; Analysis of the water molecule diffusion coefficient in the membrane of TMC/MPD and TMC/MMPD displayed that the water molecule diffusion coefficient in the membrane of TMC/MPD was much bigger than in TMC/MMPD. For the monomer of the membrane MMPD had a methyl substitution in the benzene ring, the rigidity of the membrane MMPD was stronger than MPD when the membranes consist of the same account and scale of the monomer. As a result, this did not benefit the molecules of water transitting through the membrane MMPD. The same conclusion was obtained by comparing the diffusion coefficient of the molecules of water in the membrane CFIC/MMPD and CFIC/MPD.
④NaCl molecule had the biggest diffusion coefficient in the membrane of HT/MMPD. Then it took turns to the membrane of ICIC-IPC/MPD, TMC-HT/MMPD, CFIC/MPD , TMC/MMPD, and the least in the membrane of ICIC-HT/MMPD. Only ICIC-IPC/MPD had the bigger difference between simulation result and the experiment result.
⑤In eight kinds of membranes we simulated, the diffusion coefficient value of Na+ and Cl- had much difference. The movement of NaCl molecule was restricted by the ions which moves slowly. The kind of the ions which constrained the molecule NaCl was variety in different membranes: in the membranes of TMC/MPD and CFIC/MMPD, the constrain ion was Cl-, while was Na+ in other six membranes.
⑥For membrane ICIC-HT/MPD ,the diffusion coefficients of water in the systems containing MgCl2,CaSO4 and K2SO4 were consistent with that in the NaCl system, this result demonstrated that movement of molecule of water was not affected by the kinds of salt molecules in the same membrane.
⑦The diffusion coefficients of anions (Cl-,SO42-) in the reverse osmosis membrane ICIC-HT/MPD were much bigger than cations (Mg2+,Ca2+,K+). This result was consistent with the diffusion coefficient of the anions Cl- and cations Na+ in the NaCl-H2O system, which indicated that diffusion process of salt molecule in the membrane was only restricted by the membrane species,without relationship with salt molecules itself.
引文
1. United Nations Report. Natural Resources/Water Series[R]. New York: 1985,(14): 8,125.
2.王世昌,周清,王志.海水淡化与反渗透技术的发展形势.膜科学与技术, 2003, 23(4): 162-171.
3. Wargrak K, Europe, Desalination and Water Reuse,2002,12(1):31-33.
4. Birkett J D, North America. Desalination and Water Reuse, 2002, 12(1): 24-26.
5. Tolaro Goto,East and South East Asia.Desalition and Water Reuse,2002,12(1):28-30.
6.俞三传,高从堦。膜技术在海水淡化和综合利用中的新进展,海洋通报。2001,20(1):83-87
7.查吕应,丁志斌,草传新.纳滤技术用于反渗透海水淡化预处理评述.工业用水与废水, 2003, 34(4): 18-20.
8.约瑟夫F皮津诺.水面舰船反渗透淡化技术综述.陈凤章译.机电设备, 1993, (4): 25-34.
9.蔡振雄,范志贤.膜分离海水淡化技术.集美大学学报(自然科学版), 2000, 5(1): 87-90.
10.周勇,俞三传,高从堦.反渗透复合膜(I)结构与性能.化工学报, 2006, 57(6): 1370-1373.
11.郑燕飞.膜技术缓解水危机.节能与环保(专题报道), 2004,(4): 19-20.
12.张秀芝,王树勋,王静,张雨山.膜技术在海水资源利用中的应用.海洋技术, 2005, 24(3): 128-131.
13. Sourirajan S. Reverse Osmosis. New York: Academic, 1970
14. Reid C E, Bretgon B J. J.Appl.Polym.Sci.1959, (1):133-143
15. Reid C E, Kuppers J R. J.Appl. Polym.Sci.1959, (2):264-272
16. Lonsdale H, Podall H . Reverse Osmosis Membrane Research. New York: Plenum, 1972. 43-57.
17. Spiegler K, Laird A. Principles of Desalination. New York: Academic Press,1980.
18. Sourirajan S. Reverse Osmosis and Synthetic Membranes. Ottawa, Canada: NRCC Press, 1977.
19. Turbak A. Synthetic Membranes. Vol.1 and 2, ACS Symposium Series. No.153 and 154. Washington, D.C: ACS, 1981.
20. Lonsdale H K. Synthetic Membranes: Science, Engineering and Applications. Reide Publishing Co.,1986.
21. Baker R W, Cussler E L. Membrane Seperation Systems. Park Ridge, NJ: Noyes Data Corporation,1991.
22. Lloyd Douglas R. Material Science of Synthetic Membranes. ACS Symposium Series;No.269. Washington, D.C: ACS, 1985.
23. Porter Mark C. Handbook of Industrial Membrane Technology. Park Ridge, NJ: Noyes Publications,1990.
24. Winston Ho W S, Sirkar Kamalesh K. Membrane Handbook. New York: Van Nostrand Reinhold, 1992.
25. Souriorajan S.膜分离科学与技术. 1984,4(11)
26. Strathmann H.水处理技术. 1984,10(6):39-50.
27.高从堦.水处理技术。1993,19(1):15-19.
28.朱长乐.膜科学技术.杭州:浙江大学出版社, 1992.
29. P S Fancis. Fabrication and evalution of new ultrathin reverse osmosis membranes, NTIS Report No. PB-177083, available from National Technical Information Service, U.S. Department of Commerce Springfield, VA22161,1966,(2).
30. R I, Riley, H K Lonsdale, and CRD Lyons. Composite membranes for seawater desalination by reverse osmosis. J Appl Polym Sci, 1971, 15: 1267-1281
31.朱伟平.分子模拟技术在高分子领域的应用.塑料科技, 2002, 151(5): 23-33.
32.陈正国,朱申敏,程时远.分子模拟方法及在高聚物中的应用.材料导报, 1997, 11(6: 49-51.
33. R.W. Philip, A.K.Peter,M.C. Anne. Application of molecular and materials modeling. International Technology Research Institute, 2002.
34.理查德,古柏著.王银桂等译.分子轨道从头计算法.北京:科学出版社, 1987.
35.邓良容.量子化学研究方法及其在结构化学中的应用.中山大学研究生学刊, 2000, 21(4): 61-67.
36.程新,陈亚明.量子化学计算方法在材料科学领域的初步应用.山东建材学院学报, 1994, 8(2): 1-4.
37.杨小震.分子模拟与高分子材料.北京:科学出版社,2002.
38. Andrew R. Leach. Molecular Modeling (principles and application).世界图书出版社, 2003.
39.胡小峰,许茜. CaCl2-NaCl熔盐电脱氧法制备金属Ta*.金属学报, 2006, 42(3): 285-289。
40.秦总根,涂伟萍,夏正斌.树脂乳胶成膜过程的研究现状与前景.合成树脂及塑料, 2004, 21(5): 75-78.
41.谢红梅,朱波.农田非点源氮污染研究进展.生态环境, 2003, 12(3): 349-352.
42.苏利军,卢文波.地下巷道钻爆开挖过程中炮烟扩散及通风.爆破, 2000, 17(1): 1-6.
43.孙大亮,于锡玲,王燕,顾庆天. TGS晶体生长多形性和其表面扩散生长机制.物理学报, 2000, 49(9): 1873-1877.
44.吴卫军,李炳宗,邵凯,房华. Co/Ti/Si三元薄膜固相反应外延CoSi2过程实验研究.应用科学学报, 1996, 14(4): 403-408.
45.刘桦.河口潮流和水质数值模拟研究.上海环境科学, 1997, 16 (7): 20-23.
46.罗学涛,张立同.自韧Si3N4陶瓷的高温性能特征.现代陶瓷技术, 1996, 3: 19-22.
47.邱俭,陈宜张.水通道.生命的化学, 1996, 16(1): 17-20.
48.李军,许振良,刘洪来,胡英.乙醇和水在聚酰亚胺膜中的溶解和传递现象研究(II)乙醇和水在膜中溶解和传递性质关联.化工学报, 1997,48(5): 524-531.
49.周康根,张贵情,李潜,张启修.膜法处理硫酸法钛白工艺中废酸的研究.钛工业进展, 2001, 5: 30-36.
50.丁忠伟,马润宇, Fane A G.膜蒸馏跨膜传质过程的新模型-TPKPT.高效化学工程学报, 2001, 15(4): 312-317.
51.刘业凤,王如竹.新型复合吸附干燥剂的吸附动力学特性研究.上海理工大学学报, 2006, 28(2): 107-110.
52.张德见,岑超平.污泥衍生吸附剂对铅离子的吸附机理研究.环境污染治理技术与设备, 2006, 7(2): 71-74.
53.李建新,王永川,吴来贵.炭基吸附剂脱硫过程理论模拟.内蒙古大学学报(自然科学版),1999, 30(5): 619-622.
54.赵跃强,吴争鸣,刘玮炜.超临界流体萃取天然产物传质模型.化工学报, 2006, 57(3):521-525.
55.童辉,雷家珩,李德谦. Cyanex923萃取稀土(III)的传质动力学.武汉理工大学学报, 2003, 25(6): 7-9.
56.李军,许振良,刘洪来,胡英.乙醇和水在聚酰亚胺膜中的溶解和传递现象研究(I)吸收和解吸实验.化工学报, 1997,48(5): 517-523.
57. William S. Price, Olle Soderman. Self-diffusion coefficients of some hydrocarbons in water: measurements and scaling relations. J. Phys. Chem.A,2000,104:5892-5894.
58.刘耘畦,马沛生.青霉素钾在水溶液中扩散系数的测定.化学工业与工程, 1998, 15(3): 6-14.
59. P.T.H.M.Verhallen, L.J.P.Oomen, A.J.J.M.v.d.Elsen, A.J.Kruger*, J.M.H.Furtuin. The Diffusion Coefficients of Helium, Hydrogen, Oxygen and Nitrogen in Water Determined form the Permeability of a Stagnant Liquid Layer in the Quasi-steady State. Chemical Engineering Science, 1984, 39(11): 1535-1541
60. Wu Yanxia, Ma Peisheng, Liu Yunqi, Li Shufen. Diffusion coefficients of L-proline, L-threonine and L-arginine in aqueous solutions at 25。C. Fluid Phase Equilibria, 2001, 186:27-38.
61.赵长伟,马沛生,宋小溪.脂肪族氨基酸在水溶液中扩散系数的测定与关联.化工学报, 2003, 54(12): 1689-1695.
62. Zhao Changwei, Ma Peisheng. Measurement and Correlation on Liquid Diffusion Coefficients of Aqueous L-Threonine Solutions from 298.15K to 328.15K. Chinese J. Chem. Eng., 2004, 12(2): 267-271.
63.木合马德,邹公伟,朱俊杰.流动注射分析法测定氨基酸和蛋白质扩散系数的研究.南京大学学报, 1997,(3): 369-374.
64.马友光,朱春英,何明霞,明宏伟,宓磊.氨基酸在水溶液中扩散系数的实验研究.激光技术, 2005, 29(2): 145-147.
65. Ma Youguang, Zhu Chunying, Ma Peisheng, Yu K.T.. Study on the diffusion coefficients of amino acids in aqueous solutions. J.Chem.Eng.Data, 2005, 50(4): 1192-1196.
66.马友光,段志青,朱春英,余国琮.激光实时全息干涉法测量甘氨酸水溶液扩散系数.应用激光, 2004, 24(2): 99-101.
67.朱春英,马友光,赵长伟,何明霞,马沛生,余国琮.氨基酸在水溶液中的扩散系数.化工学报, 2005, 56(2): 243-245.
68. Zhao Changwei, Li Jiding, Ma Peisheng, Xia Shuqian. Measurement of Liquid Diffusion Coefficients of Aqueous Solutions of Glycine, L-Alanine, L-Valine and L-Isoleucine by Holographic Interferometry. Chinese J. Chem. Eng., 2005, 13(2): 285-290.
69.周健,陆小华,王延儒,时钧.气体在水中的分子动力学模拟.高校化学工程学报, 2000, 14(1): 1-6.
70.周健,陆小华,王延儒,时钧.不同温度下水的分子动力学模拟.计算机与应用化学, 1999, 16(4): 241-244.
71.周健,朱宇,汪文川,陆小华,王延儒,时钧.超临界NaCl水溶液的分子动力学模拟.物理化学学报, 2002, 18(3): 207-212.
72.肖吉,陆九芳,陈健,李以圭.超临界水中气体扩散系数的分子动力学模拟.高校化学工程学报, 2001, 15(1): 6-10.
73.严六明,严琪良,刘洪来,戎宗明,胡英. Ar-Kr溶液扩散系数的分子动力学模拟及其与温度的关系.高等学校化学学报, 1997, 18(12): 2026-2029.
74.孙炜,陈中,黄素逸,陈金芳. L-J流体扩散系数的分子动力学模拟.华中科技大学学报(自然科学版), 2004,32(5): 84-86.
75. Kotelyanskii M, Wagner N, Paulaitis M. J Membrane Sci.1998,139: 1-16.
76. Kotelyanskii M, Wagner N, Paulaitis M. Macromolecules,1996,29: 8497-8506.
77. Kotelyanskii M, Wagner N, Paulaitis M. Computational and Theoretical Polymer Science. 1999, 9: 301-306.
78. H.Yan, S.Murad.Molecular simulation of membrane based separation of ethanolic electrolyte solutions.Fluid Phase Equilibria. 2001,183-184:279-287.
79. Chen Minbo, Zhu Ruixin, Zhou Junhong.“Variation Principle in Simulation of Diffusion Coefficient”. The Ninth National Quantum Chemistry Academic Meeting. Theory group report, Chinese Chemical Society, GuiLin, 10/2005.
80. Lim S Y, Tsotsis T T, Molecular simulation of diffusion and sorption of gases in an amorphous polymer. J. Chem. Phys., 2003, 119(1): 496-503.
81. D. Hofmann, L. Fritz, J. Ulbrich, D. Paul. Molecular simulation of small molecule diffusion and solution in dense amorphous polysiloxanes and polyimides. computational and theoretical polymer science, 2000, 10: 419-436.
82. Tocci E, Bellacchio E, Russo N, Drioli E. Diffusion of gases in PEEKs membranes:Molecular dynamics simulations. Journal of membrane science,2002, 206: 389-398
83. Charati S G, Stern S A, Diffusion of gas in silicone polymer: molecular dynamics simulation. Macromolecules, 1998, 31: 5529-5535
84. Hofmann D, Fritz L, Ulbrich J, Detailed-atomistic molecular modeling of small molecule diffusion and solution processes in polymeric membrane materials. Macromol.Theory Simul., 2000, 9: 293-327
85. Tung K L, Lu K T, Effect of tacticity of PMMA on gas transport through membranes: MD and MC simulation studies. J. Membr. Sci., 2006, 272: 37-49
86. Karayiannis N C, Mavrantzas V G, Theodorou D N, Detailed atomistic simulation of the segmental dynamics and barrier properties of amorphous poly(ethyleneterephthalate) and poly(ethylene isophthalate). Macromolecules, 2004,37: 2978-2995
87. Pavel D, Shanks R, Molecular dynamics simulation of diffusion of O2 and CO2 in amorphous poly(ethylene terephthalate) and related polyesters. Polymer, 2003, 44:6713-6724
88. Pavel D, Shanks R, Molecular dynamics simulation of diffusion of O2 and CO2 in blends of amorphous poly(ethylene terephthalate) and related polyesters. Polymer,2005, 46: 6135-6147
89. Ronova I A, Rozhkov E M, Alentiev A Y, Yampolskii Yu. P.. Occupied and accessible volumes in glassy polymers and their relationship with gas permeation parameters. Macromol. Theory Simul., 2003, 12: 425-439
90. Xiao-Yan Wang, P. J. in‘t Veld, Y. Lu, B. D. Freeman, I. C. Sanchez. A molecular simulation study of cavity size distributions and diffusion in para and meta isomers. Polymer, 2005, 46, 9155-9161.
91. Nagel, C., Schmidtke, E., Gunther-Schade, K., Hofmann, D., Fritsch, D., Strunskus, T., Faupel, F. Free volume distributions in glassy polymer membranes: comparison between molecular modeling and experiments. Macromolecules, 2000, 33: 2242-2248
92. Wang X Y, Lee K M, Lu Y, Cavity size distributions in high free volume glassy polymers by molecular simulation. Polymer, 2004, 45: 3907-3912
93. Raptis V E, Economou I G, Theodorou D N, Molecular dynamics simulation of structure and thermodynamic properties of poly(dimethylsilamethylene) and hydrocarbon solubility therein: Toward the development of novel membrane materials for hydrocarbon separation. Macromolecules, 2004, 37: 1102-1112
94. van der Vegt N F A, Temperature dependence of gas transport in polymer melt:molecular dynamics simulation of CO2 in polyethylene. Macromolecules, 2000, 33:3153-3160
95. Laguna M F ,Guzman J, Riande E, Saiz, E. Experimental and Simulation Studies on the Transport of Argon through Poly(pentaerythritoltribenzoate acrylate). Macromolecules, 1998, 31 (21): 7488-7494
96. Bharadwaj R K, Boyd R H, Small molecule penetrant diffusion in aromatic polyesters: a molecular dynamics simulation study. Polymer, 1999, 40, 4229-4236
97. Muller-Plathe F. Molecular dynamics simulation of gas transport in amorphous polypropylene. J. Chem. Phys., 1992, 96(4): 3200-3205
98. Heuchel M, Hofmann D, Pullumbi P, Molecular modeling of small-molecule permeation in polyimides and its correlation to free-volume distributions. Macromolecules, 2004, 37: 201-214
99. Lydia Fritz, Dieter Hofmann*. Behaviour of water/ethanol mixtures in the interfacial region of different polysiloxane membranes-a molecular dynamics simulation study. Polymer, 1998,39(12):2531-2536.
100.俞庆森,朱龙观.分子设计导论.北京:高等教育出版社,2000.
101.陆小华,王俊,朱宇,丁皓,周健,王延儒,时钧.分子动力学模拟研究流体微观结构和扩散性质.南京工业大学学报. 2002, 24(1): 7-11.
102. Kubo, R., Toda, M., Hashitsume, N.. Statistitcal Physics II - Nonequilibrium Statistitcal Mechanics. Berlin : Springer-Verlag, 1985.
103. Zubalev, D. N.,《非平衡统计热力学》.北京:高等教育出版社, 1982.
104.苏汝铿,《统计物理学》.第二版.北京:高等教育出版社, 2004.
105.李如生,《平衡和非平衡统计力学》.北京:清华大学出版社, 1995.
106.徐南平,高从堦,时钧.我国膜领域的重大需求与关键问题.中国有色金属学报,2004,5:327-331.
107. Richard I.Urama, Benito J. Marinas. Mechanistic interpretation of solute permeation through a fully aromatic polyamide reverse osmosis membrane. Journal of Membrane Science, 1997,123:267-280.
108. S. Murad, J. Lin. Molecular modeling of fluid separations using membranes: effect of molecular forces on mass transfer rates. Chemical Engineering Journal, 1999,74:99-108.
109. R. Madhusudan, J. Lin, S. Murad. Molecular simulations of electro-osmosis in fluid mixtures using semi-permeable membranes. Fluid Phase Equilibria, 1998,150-151:97-105.
110. J.G. Powles, S. Murad. The simulation of semi-permeable membranes-Osmosis, reverse osmosis and electro-osmosis in electrolyte solutions. Journal of molecular liquids, 1998,78:225-231.
111. Matthias Heuchel, Dieter Hofmann, Pluton Pullumbi. Molecular modeling of small-milecule permeation in polyimides and its correlation to free-volume Distributions. Macromolecules, 2004,37:201-214.
112. Pekka Mark, Lennart Nilsson. Structure and dynamics of the TIP3P, SPC and SPC/E water models at 298K. J. Phys. Chem. A, 2001,105 : 9954-9960.
113. Myco Umemura, Soichi Hayashi, Tetsuo Nakagawa, Hiroshi Urakawa, Kanji Kajiwara. Structure of water molecules in aqueous maltose and cellobiose solutions using molecular dynamics simulation. II. Dynamics. Journal of molecular struture(Theochem), 2003,636:215-228.
114. M.C. Donnamaria, E. Caffarena, J. de Xammar Oro. Dynamical properties of water in aqueous solutions of L-ascorbic and dehydroascorbic acid; molecular dynamics simulations. Journal of molecular structure(Theochem),1999,489:19-27.
115. Tonglei Li, Dane O. Kildsig, Kinam Park. Computer simulation of molecular diffusion in amorphous polymers. Journal of controlled release, 1997,48:57-66.
116. H.Sun. COMPASS: An ab Initio Force-Field Optimized for Condensed-Phase Applications-Overview with Details on Alkane and Benzene Compounds. J.Phys.Chem.B, 1998,139,7338-7364.
117. CERIUS2, version 4.6, Accelrys Inc., San Diego, CA.
118. Siewert J. Marrink, Herman J. C. Berendsen. Permeation process of small molecules across lipid membranes studied by molecular dynamics simulations. J. Phys. Chem.,1996,100: 16729-16738.
119. Elena Tocci, Emanuele Bellacchio, Nino Russo, Enrico Drioli. Diffusion of gases in PEEKs membranes: molecular dynamics simulations. Journal of membrane science, 2002,206:389-398.
120. Sylvie Neyertz, Anthony Douanne, David Brown. A molecular dynamics simulation study of surface effects on gas permeation in free-standing polyimide membranes. Journal of membrane science, 2006,280:517-529.
121. H. Takeuchi. Molecular dynamics simulation of diffusion of simple gas molecules in a short chain polymer. J. Chem. Phys., 1990,92:5643-5652.
122. H.Takeuchi. A jump motion of small molecules in a glassy polymers: A molecular dynamics simulation. J. Chem. Phys.,1990,93:2062-2067.
123. A. A. Gusev, F. Mueller-Plathe, W. F. van Gunsteren, U. W. Suter. Dynamics of small molecules in bulk polymers. Adv. Polym. Sci.,1994,116:209-246.
124. R. M. Sok, H. J. C. Berendsen. W. F. van Gunsteren. Molecular dynamics simulation of the transport of small molecules across a polymer membrane. J. Chem. Phys.,1992,96:4699.
125. H. Takeuchi, K. Okazaki. Dynamics of small molecules in a dense polymer matrix: molecular dynamics studies. Mol. Simulations, 1996,16:59-74.
126.时钧,袁权,高从堦.膜技术手册.化学工业出版社,2001.
127.李炳瑞.结构化学.高等教育出版社,2004.
128.徐光宪.物质结构.人民教育出版社, 1978.
129.江逢霖.量子化学原理.复旦大学出版社, 1987.
130 Schaefer M., Karplus M.. A Comprehensive Analytical Treatment of Continuum Electrostatics。. Phys. Chem.; 1996; 100(5); 1578-1599.
2.王世昌,周清,王志.海水淡化与反渗透技术的发展形势.膜科学与技术, 2003, 23(4): 162-171.
3. Wargrak K, Europe, Desalination and Water Reuse,2002,12(1):31-33.
4. Birkett J D, North America. Desalination and Water Reuse, 2002, 12(1): 24-26.
5. Tolaro Goto,East and South East Asia.Desalition and Water Reuse,2002,12(1):28-30.
6.俞三传,高从堦。膜技术在海水淡化和综合利用中的新进展,海洋通报。2001,20(1):83-87
7.查吕应,丁志斌,草传新.纳滤技术用于反渗透海水淡化预处理评述.工业用水与废水, 2003, 34(4): 18-20.
8.约瑟夫F皮津诺.水面舰船反渗透淡化技术综述.陈凤章译.机电设备, 1993, (4): 25-34.
9.蔡振雄,范志贤.膜分离海水淡化技术.集美大学学报(自然科学版), 2000, 5(1): 87-90.
10.周勇,俞三传,高从堦.反渗透复合膜(I)结构与性能.化工学报, 2006, 57(6): 1370-1373.
11.郑燕飞.膜技术缓解水危机.节能与环保(专题报道), 2004,(4): 19-20.
12.张秀芝,王树勋,王静,张雨山.膜技术在海水资源利用中的应用.海洋技术, 2005, 24(3): 128-131.
13. Sourirajan S. Reverse Osmosis. New York: Academic, 1970
14. Reid C E, Bretgon B J. J.Appl.Polym.Sci.1959, (1):133-143
15. Reid C E, Kuppers J R. J.Appl. Polym.Sci.1959, (2):264-272
16. Lonsdale H, Podall H . Reverse Osmosis Membrane Research. New York: Plenum, 1972. 43-57.
17. Spiegler K, Laird A. Principles of Desalination. New York: Academic Press,1980.
18. Sourirajan S. Reverse Osmosis and Synthetic Membranes. Ottawa, Canada: NRCC Press, 1977.
19. Turbak A. Synthetic Membranes. Vol.1 and 2, ACS Symposium Series. No.153 and 154. Washington, D.C: ACS, 1981.
20. Lonsdale H K. Synthetic Membranes: Science, Engineering and Applications. Reide Publishing Co.,1986.
21. Baker R W, Cussler E L. Membrane Seperation Systems. Park Ridge, NJ: Noyes Data Corporation,1991.
22. Lloyd Douglas R. Material Science of Synthetic Membranes. ACS Symposium Series;No.269. Washington, D.C: ACS, 1985.
23. Porter Mark C. Handbook of Industrial Membrane Technology. Park Ridge, NJ: Noyes Publications,1990.
24. Winston Ho W S, Sirkar Kamalesh K. Membrane Handbook. New York: Van Nostrand Reinhold, 1992.
25. Souriorajan S.膜分离科学与技术. 1984,4(11)
26. Strathmann H.水处理技术. 1984,10(6):39-50.
27.高从堦.水处理技术。1993,19(1):15-19.
28.朱长乐.膜科学技术.杭州:浙江大学出版社, 1992.
29. P S Fancis. Fabrication and evalution of new ultrathin reverse osmosis membranes, NTIS Report No. PB-177083, available from National Technical Information Service, U.S. Department of Commerce Springfield, VA22161,1966,(2).
30. R I, Riley, H K Lonsdale, and CRD Lyons. Composite membranes for seawater desalination by reverse osmosis. J Appl Polym Sci, 1971, 15: 1267-1281
31.朱伟平.分子模拟技术在高分子领域的应用.塑料科技, 2002, 151(5): 23-33.
32.陈正国,朱申敏,程时远.分子模拟方法及在高聚物中的应用.材料导报, 1997, 11(6: 49-51.
33. R.W. Philip, A.K.Peter,M.C. Anne. Application of molecular and materials modeling. International Technology Research Institute, 2002.
34.理查德,古柏著.王银桂等译.分子轨道从头计算法.北京:科学出版社, 1987.
35.邓良容.量子化学研究方法及其在结构化学中的应用.中山大学研究生学刊, 2000, 21(4): 61-67.
36.程新,陈亚明.量子化学计算方法在材料科学领域的初步应用.山东建材学院学报, 1994, 8(2): 1-4.
37.杨小震.分子模拟与高分子材料.北京:科学出版社,2002.
38. Andrew R. Leach. Molecular Modeling (principles and application).世界图书出版社, 2003.
39.胡小峰,许茜. CaCl2-NaCl熔盐电脱氧法制备金属Ta*.金属学报, 2006, 42(3): 285-289。
40.秦总根,涂伟萍,夏正斌.树脂乳胶成膜过程的研究现状与前景.合成树脂及塑料, 2004, 21(5): 75-78.
41.谢红梅,朱波.农田非点源氮污染研究进展.生态环境, 2003, 12(3): 349-352.
42.苏利军,卢文波.地下巷道钻爆开挖过程中炮烟扩散及通风.爆破, 2000, 17(1): 1-6.
43.孙大亮,于锡玲,王燕,顾庆天. TGS晶体生长多形性和其表面扩散生长机制.物理学报, 2000, 49(9): 1873-1877.
44.吴卫军,李炳宗,邵凯,房华. Co/Ti/Si三元薄膜固相反应外延CoSi2过程实验研究.应用科学学报, 1996, 14(4): 403-408.
45.刘桦.河口潮流和水质数值模拟研究.上海环境科学, 1997, 16 (7): 20-23.
46.罗学涛,张立同.自韧Si3N4陶瓷的高温性能特征.现代陶瓷技术, 1996, 3: 19-22.
47.邱俭,陈宜张.水通道.生命的化学, 1996, 16(1): 17-20.
48.李军,许振良,刘洪来,胡英.乙醇和水在聚酰亚胺膜中的溶解和传递现象研究(II)乙醇和水在膜中溶解和传递性质关联.化工学报, 1997,48(5): 524-531.
49.周康根,张贵情,李潜,张启修.膜法处理硫酸法钛白工艺中废酸的研究.钛工业进展, 2001, 5: 30-36.
50.丁忠伟,马润宇, Fane A G.膜蒸馏跨膜传质过程的新模型-TPKPT.高效化学工程学报, 2001, 15(4): 312-317.
51.刘业凤,王如竹.新型复合吸附干燥剂的吸附动力学特性研究.上海理工大学学报, 2006, 28(2): 107-110.
52.张德见,岑超平.污泥衍生吸附剂对铅离子的吸附机理研究.环境污染治理技术与设备, 2006, 7(2): 71-74.
53.李建新,王永川,吴来贵.炭基吸附剂脱硫过程理论模拟.内蒙古大学学报(自然科学版),1999, 30(5): 619-622.
54.赵跃强,吴争鸣,刘玮炜.超临界流体萃取天然产物传质模型.化工学报, 2006, 57(3):521-525.
55.童辉,雷家珩,李德谦. Cyanex923萃取稀土(III)的传质动力学.武汉理工大学学报, 2003, 25(6): 7-9.
56.李军,许振良,刘洪来,胡英.乙醇和水在聚酰亚胺膜中的溶解和传递现象研究(I)吸收和解吸实验.化工学报, 1997,48(5): 517-523.
57. William S. Price, Olle Soderman. Self-diffusion coefficients of some hydrocarbons in water: measurements and scaling relations. J. Phys. Chem.A,2000,104:5892-5894.
58.刘耘畦,马沛生.青霉素钾在水溶液中扩散系数的测定.化学工业与工程, 1998, 15(3): 6-14.
59. P.T.H.M.Verhallen, L.J.P.Oomen, A.J.J.M.v.d.Elsen, A.J.Kruger*, J.M.H.Furtuin. The Diffusion Coefficients of Helium, Hydrogen, Oxygen and Nitrogen in Water Determined form the Permeability of a Stagnant Liquid Layer in the Quasi-steady State. Chemical Engineering Science, 1984, 39(11): 1535-1541
60. Wu Yanxia, Ma Peisheng, Liu Yunqi, Li Shufen. Diffusion coefficients of L-proline, L-threonine and L-arginine in aqueous solutions at 25。C. Fluid Phase Equilibria, 2001, 186:27-38.
61.赵长伟,马沛生,宋小溪.脂肪族氨基酸在水溶液中扩散系数的测定与关联.化工学报, 2003, 54(12): 1689-1695.
62. Zhao Changwei, Ma Peisheng. Measurement and Correlation on Liquid Diffusion Coefficients of Aqueous L-Threonine Solutions from 298.15K to 328.15K. Chinese J. Chem. Eng., 2004, 12(2): 267-271.
63.木合马德,邹公伟,朱俊杰.流动注射分析法测定氨基酸和蛋白质扩散系数的研究.南京大学学报, 1997,(3): 369-374.
64.马友光,朱春英,何明霞,明宏伟,宓磊.氨基酸在水溶液中扩散系数的实验研究.激光技术, 2005, 29(2): 145-147.
65. Ma Youguang, Zhu Chunying, Ma Peisheng, Yu K.T.. Study on the diffusion coefficients of amino acids in aqueous solutions. J.Chem.Eng.Data, 2005, 50(4): 1192-1196.
66.马友光,段志青,朱春英,余国琮.激光实时全息干涉法测量甘氨酸水溶液扩散系数.应用激光, 2004, 24(2): 99-101.
67.朱春英,马友光,赵长伟,何明霞,马沛生,余国琮.氨基酸在水溶液中的扩散系数.化工学报, 2005, 56(2): 243-245.
68. Zhao Changwei, Li Jiding, Ma Peisheng, Xia Shuqian. Measurement of Liquid Diffusion Coefficients of Aqueous Solutions of Glycine, L-Alanine, L-Valine and L-Isoleucine by Holographic Interferometry. Chinese J. Chem. Eng., 2005, 13(2): 285-290.
69.周健,陆小华,王延儒,时钧.气体在水中的分子动力学模拟.高校化学工程学报, 2000, 14(1): 1-6.
70.周健,陆小华,王延儒,时钧.不同温度下水的分子动力学模拟.计算机与应用化学, 1999, 16(4): 241-244.
71.周健,朱宇,汪文川,陆小华,王延儒,时钧.超临界NaCl水溶液的分子动力学模拟.物理化学学报, 2002, 18(3): 207-212.
72.肖吉,陆九芳,陈健,李以圭.超临界水中气体扩散系数的分子动力学模拟.高校化学工程学报, 2001, 15(1): 6-10.
73.严六明,严琪良,刘洪来,戎宗明,胡英. Ar-Kr溶液扩散系数的分子动力学模拟及其与温度的关系.高等学校化学学报, 1997, 18(12): 2026-2029.
74.孙炜,陈中,黄素逸,陈金芳. L-J流体扩散系数的分子动力学模拟.华中科技大学学报(自然科学版), 2004,32(5): 84-86.
75. Kotelyanskii M, Wagner N, Paulaitis M. J Membrane Sci.1998,139: 1-16.
76. Kotelyanskii M, Wagner N, Paulaitis M. Macromolecules,1996,29: 8497-8506.
77. Kotelyanskii M, Wagner N, Paulaitis M. Computational and Theoretical Polymer Science. 1999, 9: 301-306.
78. H.Yan, S.Murad.Molecular simulation of membrane based separation of ethanolic electrolyte solutions.Fluid Phase Equilibria. 2001,183-184:279-287.
79. Chen Minbo, Zhu Ruixin, Zhou Junhong.“Variation Principle in Simulation of Diffusion Coefficient”. The Ninth National Quantum Chemistry Academic Meeting. Theory group report, Chinese Chemical Society, GuiLin, 10/2005.
80. Lim S Y, Tsotsis T T, Molecular simulation of diffusion and sorption of gases in an amorphous polymer. J. Chem. Phys., 2003, 119(1): 496-503.
81. D. Hofmann, L. Fritz, J. Ulbrich, D. Paul. Molecular simulation of small molecule diffusion and solution in dense amorphous polysiloxanes and polyimides. computational and theoretical polymer science, 2000, 10: 419-436.
82. Tocci E, Bellacchio E, Russo N, Drioli E. Diffusion of gases in PEEKs membranes:Molecular dynamics simulations. Journal of membrane science,2002, 206: 389-398
83. Charati S G, Stern S A, Diffusion of gas in silicone polymer: molecular dynamics simulation. Macromolecules, 1998, 31: 5529-5535
84. Hofmann D, Fritz L, Ulbrich J, Detailed-atomistic molecular modeling of small molecule diffusion and solution processes in polymeric membrane materials. Macromol.Theory Simul., 2000, 9: 293-327
85. Tung K L, Lu K T, Effect of tacticity of PMMA on gas transport through membranes: MD and MC simulation studies. J. Membr. Sci., 2006, 272: 37-49
86. Karayiannis N C, Mavrantzas V G, Theodorou D N, Detailed atomistic simulation of the segmental dynamics and barrier properties of amorphous poly(ethyleneterephthalate) and poly(ethylene isophthalate). Macromolecules, 2004,37: 2978-2995
87. Pavel D, Shanks R, Molecular dynamics simulation of diffusion of O2 and CO2 in amorphous poly(ethylene terephthalate) and related polyesters. Polymer, 2003, 44:6713-6724
88. Pavel D, Shanks R, Molecular dynamics simulation of diffusion of O2 and CO2 in blends of amorphous poly(ethylene terephthalate) and related polyesters. Polymer,2005, 46: 6135-6147
89. Ronova I A, Rozhkov E M, Alentiev A Y, Yampolskii Yu. P.. Occupied and accessible volumes in glassy polymers and their relationship with gas permeation parameters. Macromol. Theory Simul., 2003, 12: 425-439
90. Xiao-Yan Wang, P. J. in‘t Veld, Y. Lu, B. D. Freeman, I. C. Sanchez. A molecular simulation study of cavity size distributions and diffusion in para and meta isomers. Polymer, 2005, 46, 9155-9161.
91. Nagel, C., Schmidtke, E., Gunther-Schade, K., Hofmann, D., Fritsch, D., Strunskus, T., Faupel, F. Free volume distributions in glassy polymer membranes: comparison between molecular modeling and experiments. Macromolecules, 2000, 33: 2242-2248
92. Wang X Y, Lee K M, Lu Y, Cavity size distributions in high free volume glassy polymers by molecular simulation. Polymer, 2004, 45: 3907-3912
93. Raptis V E, Economou I G, Theodorou D N, Molecular dynamics simulation of structure and thermodynamic properties of poly(dimethylsilamethylene) and hydrocarbon solubility therein: Toward the development of novel membrane materials for hydrocarbon separation. Macromolecules, 2004, 37: 1102-1112
94. van der Vegt N F A, Temperature dependence of gas transport in polymer melt:molecular dynamics simulation of CO2 in polyethylene. Macromolecules, 2000, 33:3153-3160
95. Laguna M F ,Guzman J, Riande E, Saiz, E. Experimental and Simulation Studies on the Transport of Argon through Poly(pentaerythritoltribenzoate acrylate). Macromolecules, 1998, 31 (21): 7488-7494
96. Bharadwaj R K, Boyd R H, Small molecule penetrant diffusion in aromatic polyesters: a molecular dynamics simulation study. Polymer, 1999, 40, 4229-4236
97. Muller-Plathe F. Molecular dynamics simulation of gas transport in amorphous polypropylene. J. Chem. Phys., 1992, 96(4): 3200-3205
98. Heuchel M, Hofmann D, Pullumbi P, Molecular modeling of small-molecule permeation in polyimides and its correlation to free-volume distributions. Macromolecules, 2004, 37: 201-214
99. Lydia Fritz, Dieter Hofmann*. Behaviour of water/ethanol mixtures in the interfacial region of different polysiloxane membranes-a molecular dynamics simulation study. Polymer, 1998,39(12):2531-2536.
100.俞庆森,朱龙观.分子设计导论.北京:高等教育出版社,2000.
101.陆小华,王俊,朱宇,丁皓,周健,王延儒,时钧.分子动力学模拟研究流体微观结构和扩散性质.南京工业大学学报. 2002, 24(1): 7-11.
102. Kubo, R., Toda, M., Hashitsume, N.. Statistitcal Physics II - Nonequilibrium Statistitcal Mechanics. Berlin : Springer-Verlag, 1985.
103. Zubalev, D. N.,《非平衡统计热力学》.北京:高等教育出版社, 1982.
104.苏汝铿,《统计物理学》.第二版.北京:高等教育出版社, 2004.
105.李如生,《平衡和非平衡统计力学》.北京:清华大学出版社, 1995.
106.徐南平,高从堦,时钧.我国膜领域的重大需求与关键问题.中国有色金属学报,2004,5:327-331.
107. Richard I.Urama, Benito J. Marinas. Mechanistic interpretation of solute permeation through a fully aromatic polyamide reverse osmosis membrane. Journal of Membrane Science, 1997,123:267-280.
108. S. Murad, J. Lin. Molecular modeling of fluid separations using membranes: effect of molecular forces on mass transfer rates. Chemical Engineering Journal, 1999,74:99-108.
109. R. Madhusudan, J. Lin, S. Murad. Molecular simulations of electro-osmosis in fluid mixtures using semi-permeable membranes. Fluid Phase Equilibria, 1998,150-151:97-105.
110. J.G. Powles, S. Murad. The simulation of semi-permeable membranes-Osmosis, reverse osmosis and electro-osmosis in electrolyte solutions. Journal of molecular liquids, 1998,78:225-231.
111. Matthias Heuchel, Dieter Hofmann, Pluton Pullumbi. Molecular modeling of small-milecule permeation in polyimides and its correlation to free-volume Distributions. Macromolecules, 2004,37:201-214.
112. Pekka Mark, Lennart Nilsson. Structure and dynamics of the TIP3P, SPC and SPC/E water models at 298K. J. Phys. Chem. A, 2001,105 : 9954-9960.
113. Myco Umemura, Soichi Hayashi, Tetsuo Nakagawa, Hiroshi Urakawa, Kanji Kajiwara. Structure of water molecules in aqueous maltose and cellobiose solutions using molecular dynamics simulation. II. Dynamics. Journal of molecular struture(Theochem), 2003,636:215-228.
114. M.C. Donnamaria, E. Caffarena, J. de Xammar Oro. Dynamical properties of water in aqueous solutions of L-ascorbic and dehydroascorbic acid; molecular dynamics simulations. Journal of molecular structure(Theochem),1999,489:19-27.
115. Tonglei Li, Dane O. Kildsig, Kinam Park. Computer simulation of molecular diffusion in amorphous polymers. Journal of controlled release, 1997,48:57-66.
116. H.Sun. COMPASS: An ab Initio Force-Field Optimized for Condensed-Phase Applications-Overview with Details on Alkane and Benzene Compounds. J.Phys.Chem.B, 1998,139,7338-7364.
117. CERIUS2, version 4.6, Accelrys Inc., San Diego, CA.
118. Siewert J. Marrink, Herman J. C. Berendsen. Permeation process of small molecules across lipid membranes studied by molecular dynamics simulations. J. Phys. Chem.,1996,100: 16729-16738.
119. Elena Tocci, Emanuele Bellacchio, Nino Russo, Enrico Drioli. Diffusion of gases in PEEKs membranes: molecular dynamics simulations. Journal of membrane science, 2002,206:389-398.
120. Sylvie Neyertz, Anthony Douanne, David Brown. A molecular dynamics simulation study of surface effects on gas permeation in free-standing polyimide membranes. Journal of membrane science, 2006,280:517-529.
121. H. Takeuchi. Molecular dynamics simulation of diffusion of simple gas molecules in a short chain polymer. J. Chem. Phys., 1990,92:5643-5652.
122. H.Takeuchi. A jump motion of small molecules in a glassy polymers: A molecular dynamics simulation. J. Chem. Phys.,1990,93:2062-2067.
123. A. A. Gusev, F. Mueller-Plathe, W. F. van Gunsteren, U. W. Suter. Dynamics of small molecules in bulk polymers. Adv. Polym. Sci.,1994,116:209-246.
124. R. M. Sok, H. J. C. Berendsen. W. F. van Gunsteren. Molecular dynamics simulation of the transport of small molecules across a polymer membrane. J. Chem. Phys.,1992,96:4699.
125. H. Takeuchi, K. Okazaki. Dynamics of small molecules in a dense polymer matrix: molecular dynamics studies. Mol. Simulations, 1996,16:59-74.
126.时钧,袁权,高从堦.膜技术手册.化学工业出版社,2001.
127.李炳瑞.结构化学.高等教育出版社,2004.
128.徐光宪.物质结构.人民教育出版社, 1978.
129.江逢霖.量子化学原理.复旦大学出版社, 1987.
130 Schaefer M., Karplus M.. A Comprehensive Analytical Treatment of Continuum Electrostatics。. Phys. Chem.; 1996; 100(5); 1578-1599.