甲烷部分氧化过程中强制振荡的动力学Monte Carlo模拟

为了有效控制动力学振荡行为以提高反应的转化率, 利用Monte Carlo 方法研究了甲烷催化部分氧化过程中的强制振荡行为, 探讨了原料气周期性变化对振荡动力学和转化率的影响. 研究表明原料气周期性变化不仅可以有效调控振荡的动力学行为, 产生如短周期振荡和双峰振荡等特殊动力学过程, 而且还可以提高反应的转化率. 当强制过程的周期从T/3 增大到2T(T为自发振荡过程的平均周期), 振荡过程从短周期振荡变化为双峰振荡. 对反应过程中CO的转化率进行了计算, 结果表明原料气周期性变化可以有效提高反应的转化率.振荡动力学的改变和转化率的提高主要是因为强制振荡过程使得催化剂表面发生了从氧化态向还原态的转变.     

不同孪晶界密度银纳米线拉伸形变行为的分子动力学模拟

采用分子动力学方法模拟了不同孪晶界密度银纳米线的拉伸形变行为, 分析了孪晶界密度对多晶银纳米线屈服强度、弹性模量和塑性变形机理的影响. 在弹性形变区域, 孪晶界的存在对杨氏模量变化的作用不明显. 在塑性形变阶段, 首先从表面边缘开始产生位错成核, 然后延伸并受阻于孪晶界. 在进一步拉伸载荷作用下, 孪晶界将作为位错源产生新的位错. 模拟结果表明, 银纳米线的强度与孪晶界和晶粒的尺寸有关. 孪晶界密度较小(即晶粒的长径比大于1)时, 此纳米线的屈服应力比单晶纳米线还要小, 只有当孪晶界密度较大时(即晶粒的长径比小于1), 孪晶界使得纳米线得到强化. 综合分析了孪晶界和晶粒尺寸对银纳米线的影响, 为构建高强度金属纳米线打下基础. 最后讨论了温度和拉伸速度对孪晶纳米线屈服应力所产生的影响, 随着温度的升高, 孪晶纳米线与单晶纳米线的屈服应力差先升高后趋于稳定; 当拉伸速度逐渐增大, 孪晶纳米线与单晶纳米线的屈服应力差先稳定后增大.     

分子动力学模拟研究[Bmim][PF6]+水+醇三元体系的微观结构和分子间相互作用

研究离子液体体系的微观结构和分子间相互作用具有重要意义. 本文对1-丁基-3-甲基咪唑六氟磷酸盐([Bmim][PF₆])+水+乙醇和[Bmim][PF₆]+水+异丙醇三元体系进行了分子模拟研究, 计算了径向分布函数和不同组成的水-醇混合溶剂与离子液体阴阳离子间的相互作用能, 并将其分解为库仑相互作用能和Lennard-Jones(LJ)势能. 在此基础上, 研究了溶液体系的微观结构、分子间相互作用和相行为. 结果表明, 水倾向于与离子液体阴离子和阳离子极性部分作用, 醇倾向于与阴离子和阳离子非极性部分作用; 库仑力主导阴离子-溶剂相互作用, 色散力主导阳离子-溶剂相互作用, 阴阳离子的缔合状态对色散力影响较小, 对库仑力的影响非常显著.     

大环多胺金属模拟酶催化水解梭曼机理研究

为解决生物酶在战场应用环境中易失活且环境耐受性差的问题,我们合成了大环多胺双锌模拟酶配合物-2,6-双{[双(2-羟乙基)氨基]甲基}-4-甲基酚(Zn L2).Zn L2对GD(Soman,梭曼)具有较高的催化活性,能够实现多次催化循环,并且在Zn L2与GD化学计量比为(0.56∶1)的条件下水解率为100%.使用高效液相色谱-质谱与核磁共振确定催化产物为甲氟膦酸酯、甲氟膦酸、甲基膦酸频哪酯与甲基膦酸,说明Zn L2水解GD为双平行反应,反应途径分别为GD的P—F与P—O键的断裂.建立在实验的基础上,使用密度泛函理论(DFT)的B3LYP方法,从P—O与P—F键断裂这两个平行反应方向对催化过程进行计算.计算结果表明:Zn L2上Ot末端醇盐作为一个碱试剂去激活水分子,使其成为羟基并作为亲核试剂进攻中心磷原子,形成一个五配位的磷酰中间体,通过进一步亲核进攻水解中间体形成最终P—O与P—F键断裂产物.Zn L2催化GD水解的活化能为分别为5.6 kcal/mol(P—F键断裂)和11.5 kcal/mol(P—O键断裂),说明P—O键断裂过程是GD酶促水解反应的速控步骤.     

Simulation Methods for Functional Materials

Functional materials embrace a broad area, ranging from functional information materials to special polymers, from special chemicals for printing to matenais used in making paints. Inasmuch as most tunctional materials are particulate, the present contribution is considered pertinent to the present FORUM.     

随机结构正交展开分析的Ritz动力聚缩法

针对随机结果正交展开理论计算上的弱点 ,本文在分析扩阶矩阵特性的基础上 ,于 Ritz模态向量子空间中对扩阶方程实现动力聚缩 ,大大提高了正交展开理论对实际工程问题的分析能力。分析实例表明 :即使结构参数具有很大变异性 (如δ =0 .4 )时 ,该算法依然能理想地与 Monte Carlo法模拟结果相吻合 ,计算时间则远远小于 Monte Carlo模拟法。同时 ,分析例证再一次强化了在结构动力分析中考虑结构参数随机性的必要性     

各向异性本构关系在板料成形数值模拟中的应用

对几种能表达面内各向异性的屈服准则Hill、Barlat-Lian、Barlat进行了比较。以弹性变形服从各向同性广义虎克定律的情况下，给出了基于张量算法推导的弹塑性本构关系的一般表达式，并由此导出了相应屈服准则的弹塑性本构关系的显式表达。借助ABAQUS软件本构模块用户子程序接口，分别实现了这些屈服准则在ABAQUS的嵌入。以模拟方形盒的拉延过程为例，分析了不同的屈服准则在板料成形过程数值模拟中的应用。模拟结果表明，基于弹塑性本构关系一般表达所列出的相应屈服准则的显式表达式是正确的；在采用壳元来模拟板料成形时，采用Barlat准则的模拟结果和采用Barlat-Lian准则的结果差别不大。     

延性材料平面冲击波基本实验数值模拟的若干问题

假设屈服强度正比于剪切模量,分别构建了四种强度模型,并用于数值模拟Ly12铝的平面冲击波试验的质点速度剖面以及钨合金的平面冲击波试验的锰铜计应力记录.研究发现,基于Steinberg等的屈服强度模型比较符合实验.还对平面冲击波试验的层裂现象进行了简化的数值模拟,其中,采用了该文作者所提出的一种基于空穴聚集的层裂模型,该模型的特点是分别引入了层裂过程的早期及后期应力松弛方程.     

L-leap-accelerating the stochastic simulation of chemically reacting systems

Presented here is an L-leap method for accelerating stochastic simulation of well-stirred chemically reacting systems,in which the number of reactions occurring in a reaction channel with the largest propensity function is calculated from the leap condition and the number of reactions occurring in the other reaction channels are generated by using binomial random variables during a leap.The L-leap method can better satisfy the leap condition.Numerical simulation results indicate that the L-leap method can obtain better performance than established methods.     

促进其线性频散特征另一种形式的Bousinesq方程

Ｂｏｕｓｉｎｅｓｑ方程能够用于模拟表面重力波传播过程中的折射、绕射、反射以及浅化，非线性作用等现象．用不同垂直积分方法所得到的二维Ｂｏｕｓｓｉｎｅｓｑ方程形式具有不同的线性频散特征．采用两个不同的水深层的水平速度变量组合，推导出一个新形式的Ｂｏｕｓｉｎｅｓｑ方程．通过对其参数的设置可得到精确的线性频散解Ｐａｄｅ近似４阶精度．其适用范围已由原来的浅水，向深水拓进．相速误差小于２％，其拓展适用范围可达到０８个波长水深．应用所得到的新型Ｂｏｕｓｉｎｅｓｑ方程，采用有限差分法，对经典工况进行了数值模拟，其计算结果表明，计算值与物模实验值吻合较好．这说明本文新形式的Ｂｏｕｓｓｉｎｅｓｑ方程对变水深非线性效应所产生的能量频散有着较为精确的描述