多相多组分化学反应平衡和相平衡计算的遗传算法

引　言多相多组分化学反应平衡和相平衡的计算是热力学研究的基础问题 ,同时也是化工流程模拟的重要环节 ,建立快速、稳定的求解方法具有特别重要的理论和现实意义 早期的以及目前应用较多的平衡常数法 (K值法 )的缺点是要预先知道体系的组分数、相态和具体反应 ,不适用于相态和反应未知的复杂体系 .Gibbs自由能最小法虽然具有通用性 ,然而需要求解多变量和复杂约束的非线性规划问题 对于非理想体系 ,涉及逸度、逸度系数的计算 ,这还是一个含有大量局部最优解的非凸问题[1,2 ] ,所以传统的数学规划算法不能保证得到全局最优解 文献 [3]…     

应用遗传算法求解含化学反应体系的相平衡

基于Gibbs自由能最小原理 ,提出一种求解含化学反应体系相平衡的方法。采用反应进度和相分率消去等式约束 ,将含化学反应体系的相平衡计算转化为无约束的最优化问题 ,并应用遗传算法求解。采用本文提出的方法 ,能够避免计算中负数的组分摩尔数或摩尔分数的出现 ,确保得到组分的摩尔数或摩尔分数为非负的全局最优解。以聚甲醛反应体系的汽液平衡计算为例 ,表明本文所提出的方法是可行和有效的     

化学反应自发性判据的反例推导

文章结合常见的化学反应实例,采用反例推导的方式,逐项引出影响化学反应方向的因素,经过综合,最后给出准确的化学反应自发性判据公式。通过具体案例引起学习兴趣,并采用反例法加深同学们对判据因素的印象,以期达到知其所以然的教学效果。     

应用广义简约梯度算法增强遗传算法求解化工冶金相平衡

研究可靠的相平衡计算对化工冶金过程有很重要意义,本文采用广义简约梯度法（GRG）来增强遗传算法（GA）用于相平衡计算过程,在遗传算法演化快结束时引入GRG算子,演化结束后使用GRG精修结果,克服了GA局部搜索能力不强的缺点,加快了搜索速度.KCl-FeCl₂体系相平衡计算结果显示这种混合算法（Hybrid Algorithms）能够提高GA效率的同时保留了全局收敛的优点,因此其在化工冶金相平衡计算中将有广泛的应用前景.     

非理想多相共存体系的平衡计算

对于给定温度和压力的系统，平衡状态时其吉布斯自由能最小。根据热力学基本关系式，在多组分体系中元素和电荷守恒的限制条件下，最小化吉布斯函数，通过SQP方法来求解此有约束的优化问题，通过迭代法求解反应后的温度。该方法可用于复杂的非理想、多相流体、电解液或非电解液混合物的相平衡计算，平衡计算结果包括温度、密度、组分浓度等。对于某新型热动力系统中使用的一种含铝火药，应用本文方法，获得了化学平衡计算结果。     

基于BOPPPS改进模型对“化学反应方向”的教学设计

化学反应方向是在化学热力学中研究的关键内容,是有关化工反应的重要问题,对判断反应发生的可能性具有重大意义。基于BOPPPS改进模型,对“化学反应方向”的教学设计进行了优化,以加深学员对吉布斯自由能的理解,提高应用吉布斯自由能判据判断化学反应方向的能力。     

含化学反应的相平衡求解方法现状

综述了含化学反应的相平衡求解方法，并对其现状作了分析。从应用角度看，Ｇｉｂｂｓ自由能最小法是使用最悠久的方法，比较成熟，但收敛慢，对初值要求高，虽然作了一些改进，对复杂体系的计算仍难有通用性，Ｋ－值法收敛快，对于受约束的反应系统和非等温等复杂问题的处理更为直观简便，计算机程序编制也更容易，特别是ＫＺ法使用了内外两层循环，从而有可能将其通用化，为含化学反应过程的模拟计算奠定了基础     

混合整数非线性规划法求解多组分复杂体系相平衡问题

基于Gibbs自由能最小化原理,提出一种求解多组分复杂体系相平衡问题的混合整数非线性规划(MlNLP)方法.通过建立MlNLP模型,将相平衡计算问题转化为有约束的双层最优化问题,内层用来计算给定相组合结构下的平衡组成,外层则是搜索最优的相组合结构,整个求解过程通过遗传算法来实现.所建立的方法只需要求解最小化问题,就可以同时给出多相平衡时存在的相数、相态及各相组成.以苯-乙氰-水和甲醇合成两个体系的相平衡计算为例,结果表明该方法是可行和有效的.     

Ni/SiO2溶胶稳定性及凝胶材料的物相化学结构

以甲基三乙氧基硅烷为疏水性前驱物,Ni(NO_3)_2·6H_2O为镍源,采用溶胶-凝胶法制备镍掺杂型SiO_2杂化(Ni/SiO_2)溶胶,研究镍掺杂对SiO_2溶胶黏度、密度、反应速率常数、黏滞性活化吉布斯自由能、稳定性、粒径分布和化学结构的影响以及350℃焙烧前后材料物相、化学结构的变化。结果表明:随着镍摩尔分数的增加,Ni/SiO_2溶胶的黏度、密度、平均粒径、反应速率常数和黏滞性活化吉布斯自由能逐渐增大,溶胶的Zeta电位则逐渐减小。黏滞性活化吉布斯自由能的增加使分子间的分散性降低,溶胶的稳定性减弱。Ni/SiO_2溶胶及凝胶材料中镍元素主要以Ni(NO_3)_2·6H_2O的形式存在,其化学结构主要以Si—O—Si,Si—CH_3和Si—OH键为主,在空气气氛中经350℃焙烧后,材料的物相、化学结构发生了变化, Ni~(2+)转变成了NiO和Ni—O—Si形式。     

基于多维空间的燃烧平衡产物组成的计算方法

为计算燃烧平衡产物的组成,在质量守恒关系的基础上,得到一系列基本可行解。根据吉布斯自由能,对基本可行解进行优化筛选,将筛选得到的基本可行解视为多维解空间的顶点,以吉布斯自由能降低为判据,利用基本可行解逐个与初始解线性组合,从而搜索得到最佳点,即为最终平衡产物组成。应用文献数据对计算结果的准确性进行验证。结果表明,该方法计算平衡产物的组成与文献数据基本一致,燃烧温度最大误差不超过3%。     

Simultaneous calculation of chemical and phase equilibria using convexity analysis

Chemical and phase equilibria were considered for closed multicomponent reactive systems at: (a) constant pressure and temperature; (b) constant pressure and enthalpy. Equilibrium at constant P and T was found by minimization of G, while equilibrium at constant P and H was found by maximization of S or minimization of −S, all with respect to the number of moles of each component in each phase. Both cases could be handled as optimization problems, satisfying the restrictions imposed by mole or atom balances, and non-negativity of number of moles. Convexity analyses were carried out, and the conditions were found in order to guarantee global minimum, for one liquid phase, one gas phase, and a number of solid phases. The minimum point was then found either by analytical methods or by direct minimization methods. These strategies were tested for a number of cases, with good results.     

噻吩水热裂解平衡计算及反应过程分析

采用乘子法,结合Gibbs自由能最小的概念,计算注蒸汽热采地层条件下,噻吩水热裂解反应体系的平衡组成.根据实验数据和计算结果分析反应过程,认为噻吩水热裂解过程中发生了水解、热解、加氢脱硫、加氢饱和及水煤气转换等反应.平衡计算结果表明,当水/噻吩摩尔进料比小于4时,提高进料比有利于脱硫,并生成更多气体;提高反应温度、降低反应压力有利于提高气体产量,但不利于脱硫.当水/噻吩摩尔进料比大于4时,噻吩平衡转化率达100%,温度、压力、进料比对平衡组成的影响很小.     

Parameter estimation in simultaneous phase and reaction equilibria

The estimation of parameters of models describing simultaneous phase and reaction equilibria is generally difficult to accomplish, due to strong correlation between the parameters. We show in this paper that these models possess intermediate properties between reduced and implicit models. This makes it convenient to develop ad-hoc algorithms for this class of problems which are of considerable importance in process design. We present a method for parameter estimation based on implicit differentiation, which makes it possible to use a gradient algorithm with respect to a reduced number of parameters, which greatly enhances global convergence.     

Representation of liquid-liquid equilibria for polymer-salt aqueous two-phase systems

In this work the modified Wilson equation, proposed in a previous publication, for aqueous binary polymer solutions is extended to calculate liquid-liquid equilibrium of aqueous two-phase systems containing one salt with either polyethylene glycol or ethylene oxide-propylene oxide random copolymer. The Debye-Hückel equation based on the Fowler-Guggenheim theory is used to represent the long-range electrostatic interactions of ions. The parameters are estimated from only one or two tie-lines and are used to predict the phase diagrams with different molecular weights and polymer-salt systems at different temperatures. The predictions obtained with the modified Wilson equation are in good agreement with the experimental data.     

SO2抑制二英从头合成的实验及其过程模拟

实验研究了不同温度下SO₂对二 英（PCDD/Fs）从头合成的抑制特点,运用HSC Chemistry 6.1软件对主要无机化学反应过程进行模拟,探究了SO₂的抑制机理。结果表明,SO₂对PCDD/Fs总浓度和I-TEQ浓度的抑制效果随温度的升高而增强。PCDD/Fs总浓度和I-TEQ浓度抑制率从200℃不足20%分别提高到300℃时的75.6%和77.3%及400℃时的89.2%和80.5%;SO₂的添加促进了PCDD/Fs尤其是PCDDs同系物分布的低氯代化。在从头合成的高温段和低温段,SO₂分别侧重抑制PCDFs和PCDDs。相平衡计算和Gibbs自由能计算表明,SO₂优先与催化剂CuCl₂反应生成无催化活性的CuSO₄,其次与氯源Cl₂反应,最终实现PCDD/Fs的生成抑制。模拟结果对现有硫基抑制机理进行了验证和补充,对高效、廉价抑制剂的开发及工程使用具有重要意义。     

三元硝酸熔盐高温热稳定性实验研究与机理分析

通过熔盐质量损失率和试样中NO2-含量变化状况,以及熔盐材料的热重(TG)曲线,研究了三元硝酸熔盐(53%KNO3-40%NaNO2-7%NaNO3)在空气和氮气气氛中高温条件下的热稳定性性能。结果表明,三元硝酸熔盐空气中的上限使用温度为773K,高温时存在劣化现象;而在氮气气氛中三元熔盐的热分解温度为723K。同时,从热力学和动力学角度分析得到,三元硝酸熔盐在773K以下空气中发生的反应为亚硝酸钠的分解和氧化;而在氮气气氛中三元熔盐723K以下时主要发生的是亚硝酸钠的分解反应。在氧气含量一定的情况下,氧气的扩散和亚硝酸盐的分解反应符合一级反应动力学模型。     

铝合金微弧熔凝Al2O3形核及转变热力学分析

微弧氧化是在铝合金表面原位制备陶瓷膜层的一项新技术.针对陶瓷膜层在形成过程中的两次相变,分析计算Al2O3液相中的形核过程以及γ-Al2O3向α-Al2O3转变过程Gibbs自由能的变化,分析相变驱动力与温度及过冷度的关系.研究结果表明:约在1200K左右过冷度时,可从液相Al2O3直接析出γ-Al2O3晶核,γ-Al2O3形核Gibbs自由能的变化符合T模型,α-Al2O3的形核行为可采用D-R模型计算,α-Al2O3和γ-Al2O3的均质形核率具有相同的数量级;微弧放电区域的热影响作用可造成γ-Al2O3向α-Al2O3的转变,温度为1500K时转变驱动力比1100K时增大了约12.3%.