《第七章化学过程计算机模拟简介》由会员分享，可在线阅读，更多相关《第七章化学过程计算机模拟简介（17页珍藏版）》请在金锄头文库上搜索。
1、第七章 化学过程计算机模拟简介Monte Carlo方法借计算机之助取样以求高维积分数值的技术。是一种统计 模拟方法。 不要求将模型过分简化，可以基于分子(原 子、离子)的排列和运动的模拟结果直接计算求和，以实现 宏观现象中数量的估算。分子动力学方法（Molecular dynamics method）将系统运动方程组数值积分，得到N个经典粒子系统的 相轨道，进而研究该系统的平衡热力学性质，结构动力 学性质、非平衡输运性质等。化工过程模拟化学物系结构和运动的模拟高分子热降解的高分子热降解的Monte CarloMonte Carlo法模拟法模拟7-1-1 Monte Carlo方法应用示例应用
2、示例为了研究耐热塑料，需要探讨塑料的高分子结构和耐热 性的关系。用双链高分子代替单链高分子即采取“梯状” 高分子结构，是提高塑料耐热性一条可能的途径。多骈 六元环或多骈四元环可以构成梯状高分子的骨架Tessler用Monte Carlo法研究了由六元环或四元环骈合 而成的梯状高分子的随机热降解过程。7-1-1 Monte Carlo方法应用示例应用示例Tessler的模拟结果表明：梯状高分子降解速率比 单键高分子小很多。 无论是按简单降解机理还是复杂降解 机理，梯状高分子的热稳定性都比单链高分子好。若梯状 高分子有“结构缺陷”，即某些局部有悬键的情况，则降解 的“诱导期”变短，降解加快。7-1
3、-1 Monte Carlo方法应用示例应用示例将微观粒子视为经典粒子，服从 Newton 第二 定律或若各粒子的瞬时受力已知，可用数值积分求出 运动的经典轨迹7-1-2 分子动力学方法方法原理在合适选定的时间步长 t 内，粒子可视作匀 加速直线运动加速度：位 移：粒子 i 在时间 t 内的位移 ri步长取值： t0.010.0001 ps7-1-2 分子动力学方法方法原理假定1： 有效作用势近似假定2： 周期性边界条件(Periodical Boudary Condition)基本假定困难 欲重现实际体系的统计行为，模拟 体系应有足够数量的粒子 1dm3水 31027个H2O 计算机只能处理
4、 102 104 个粒子！ 解决办法 赝无限大近似取较小的模拟体系作中心原胞，令其在空间 重复排列7-1-2 分子动力学方法方法原理二维周期性边界条件示意图 8 个近邻重复单 元包围着中心原胞 ，为其提供合理的 边界条件近似中心 原胞N102计算机实际处理的 是原胞内数量较少 的粒子 7-1-2 分子动力学方法方法原理三维周期边界实例苏氨酸水溶液模拟111nm 的立方原胞含 1 个苏氨酸 ( threonine ) 分子和 24 个水分子 F原子总数 88分子总数 25 粒子数取得过少是为求显示的直观 欲得合理的模拟结果， 原胞体积至少应增大10倍！7-1-2 分子动力学方法方法原理1.物质和
5、大分子的结构和性质根据体系的化学成分推断微观结构计算计算机辅助结构分析与静态性质计算复杂有机分子 高分子 晶体 液晶态 非晶态 溶液熔盐 界面和表面 气体 热力学性质 材料力学性质 与动态过程有关性质的计算分子内和分子间的相互作用扩散 吸附与解附 熔融 结晶 相变 分子散射 分子适 配材料应力应变 蠕变与断裂7-1-2 分子动力学方法应用示例2.化学CAD 计算机辅助有机分子和生物分子设计 计算机辅助药物设计 计算机辅助材料设计MD法可处理粒子数研究有电子转移、原子变价的过程须采用QMD严格 QMD 100 半经验 QMD 150 经验势 MD 200007-1-2 分子动力学方法应用示例例1
6、：(1,2,4)-4-(1,1-e二甲乙基)-2-烃基 环戊羰基酰胺晶体W. Linert, and F. Renz, J. Chem. Inf. Comput. Sci., 33, 776(1993)Experimentally DeterminedMD-predicted7-1-2 分子动力学方法应用示例例2：MD法模拟的熔盐体系 T *V * 相图D.A. Young & B.J. Alder, J. Chem. Phys., 73, 2434(1980)7-1-2 分子动力学方法应用示例例3：Ru-Al 合金断裂过程动态模拟C.S. Becquart , D. Kim, J.A, Ri
7、fkin, and P.C.Clapp, Mat. Sci. Engin., A170, 87(1993)断裂点周围的损坏区域7-1-2 分子动力学方法应用示例7-2 化工过程模拟化工过程流程模拟是借助于电子计算机求解描述整个 化工生产过程的数学模型，得到有关该化工过程性能 的信息。它是化工系统工程研究的范畴。化工流程稳态模拟：对一个现有的或给定的化工流程进行稳态的物料及 能量衡算，从而得到系统的特性。 化工流程动态模拟：除包括稳态模拟中的物料和能量衡算、设备尺寸和费 用计算及过程的技术经济评价外，还包括对控制系统 性能的研究，对开停车工况控制的研究，对紧急情况 下操作策略的研究等。序贯模块法（Sequential Modular Method）联立方程法（Equation Oriented Method）用一些方程来描述流程中的备单元过程，然后将所有描述 有关单元的方程组合起来，联立求解。联立模块法（Simultaneously Modular Method）7-2 化工过程模拟稳态平台模式 基本系统结构是一个平台，各种单元流程模块独立地 置于平台上，彼此可以互不干扰地工作，但又可以根据 需要自由组合以完成复杂的任务。 应用 ：（1）培训模拟器 （2）流程方案和控制方案的研究 7-2 化工过程模拟动态
