分子束研究气体分子在表面吸脱附动力学
气体动力学基础试题与答案
一、 解释下列各对名词并说明它们之间的区别与联系 1．轨线和流线 2．马赫数M 和速度系数λ 5．膨胀波和激波 二、 回答下列问题 1．膨胀波在自由表面上反射为什么波为什么 4．收敛喷管的三种流动状态分别是什么各有何特点 三、（１2分）已知压气机入口处的空气温度Ｔ１＝２８０Ｋ，压力Ｐ１=１.０bar ，在经过压气机进行可逆绝热压缩以后，使其压力升高了２５倍，即增压比Ｐ２／Ｐ１＝２５，试求压气机出口处温度和比容，压气机所需要的容积功。设比热容为常数，且比热比k=。 四、空气沿如图1所示的扩散管道流动，在截面1-1处空气的压强 5110033.1?=p N/m 2,温度ο151=t C,速度2721=V 米/秒，截面1-1的面积 1A =10厘米2，在截面2-2处空气速度降低到2V =米/秒。设空气在扩散管中 的流动为绝能等熵流动，试求：（1）进、出口气流的马赫数1M 和2M ；（2）进、出口气流总温及总压；（3）气流作用于管道内壁的力。 六、(15分)在超声速风洞的前室中空气的滞止温度为T *=288K,在喷管出口处空 气的速度V 1=530米/秒，当流过试验段 中的模型时产生正激波（如图1所示），求激波后空气的速度。
图1 第四题示意图图2 第五题示意图
一、解释下列各对名词并说明它们之间的区别与联系（共２０分，每题４分） 1．轨线和流线 答：轨线是流体质点运动的轨迹；流线是一条空间曲线，该曲线上任一点的切线与流体在同一点的速度方向一致。 区别：轨线的是同一质点不同时刻的位置所连成的曲线；流线是同一时刻不同质点运动速度矢量所连成的曲线。 联系：在定常流动中轨迹线和流线重合。 2．马赫数M和速度系数λ 答：马赫数M是气体运动速度与当地声速的比值；速度系数λ是气体运动速度与临界声速的比值。 区别：速度相同时气体的马赫数与静温有关，最大值为无限大，而速度系数于总温有关，其最大值为有限值。 联系：已知马赫数可以计算速度系数，反之亦然。 3．膨胀波和激波 答：膨胀波是超声速绕外钝角偏转或加速时所产生的压力扰动波；激波是超音速气流流动方向向内偏转所产生强压缩波。 区别：膨胀波是等熵波；激波是非等熵波。 联系：激波在自由表面上反射为膨胀波。 二、回答下列问题（共20分，每题5分）： 4．膨胀波在自由表面上反射为什么波为什么 答：膨胀波在自由表面上反射为压缩波。这是因为膨胀波在自由表面反射之前压力与自由表面一致，但经过膨胀波之后气流的压力下降，要使气流压力恢复到自由表面的压力，必须经过一道压缩波提高压力（见示意图）。
分子动力学模拟方法概述(精)
《装备制造技术》 2007年第 10期 收稿日期 :2007-08-21 作者简介 :申海兰 , 24岁 , 女 , 河北人 , 在读研究生 , 研究方向为微机电系统。 分子动力学模拟方法概述 申海兰 , 赵靖松 (西安电子科技大学机电工程学院 , 陕西西安 710071 摘要 :介绍了分子动力学模拟的基本原理及常用的原子间相互作用势 , 如Lennard-Jones 势 ; 论述了几种常用的有限差分算法 , 如 Verlet 算法 ; 说明了分子动力学模拟的几种系综及感兴趣的宏观统计量的提取。关键词 :分子动力学模拟 ; 原子间相互作用势 ; 有限差分算法 ; 系综中图分类号 :O3 文献标识码 :A 文章编号 :1672-545X(200710-0029-02 从统计物理学中衍生出来的分子动力学模拟方法 (molec- ular dynamics simulation , M DS , 实践证明是一种描述纳米科技 研究对象的有效方法 , 得到越来越广泛的重视。所谓分子动力学模拟 , 是指对于原子核和电子所构成的多体系统 , 用计算机模拟原子核的运动过程 , 从而计算系统的结构和性质 , 其中每一个原子核被视为在全部其他原子核和电子所提供的经验势场作用下按牛顿定律运动 [1]。它被认为是本世纪以来除理论分析和实验观察之外的第三种科学研究手段 , 称之为“计算机实验” 手段 [2], 在物理学、化学、生物学和材料科学等许多领域中得到广泛地应用。
根据模拟对象的不同 , 将它分为平衡态分子动力学模拟 (EM DS (和非平衡态分子动力学模拟 (NEM DS 。其中 , EM DS 是分子动力学模拟的基础 ; NEM DS 适用于非线性响应系统的模拟 [3]。下面主要介绍 EM DS 。 1分子动力学方法的基本原理 计算中根据以下基本假设 [4]: (1 所有粒子的运动都遵循经典牛顿力学规律。 (2 粒子之间的相互作用满足叠加原理。 显然这两条忽略了量子效应和多体作用 , 与真实物理系统存在一定差别 , 仍然属于近似计算。 假设 N 为模拟系统的原子数 , 第 i 个原子的质量为 m i , 位置坐标向量为 r i , 速度为 v i =r ? i , 加速度为 a i =r ?? i , 受到的作用力为 F i , 原子 i 与原子 j 之间距离为 r ij =r i -r j , 原子 j 对原子 i 的作用力为 f ij , 原子 i 和原子 j 相互作用势能为 ! (r ij , 系统总的势能为 U (r 1, r 2, K r N = N i =1! j ≠ i ! " (r ij , 所有的物理量都是随时 间变化的 , 即 A=A (t , 控制方程如下 : m i r ?? i =F i =j ≠ i
材料热力学与动力学
材料热力学与动力学 参考书目：1.<
分子动力学方法模拟基本步骤
分子动力学方法模拟基本步骤 1.第一步 即模型的设定，也就是势函数的选取。势函数的研究和物理系统上对物质的描述研究息息相关。最早是硬球势，即小于临界值时无穷大，大于等于临界值时为零。常用的是LJ势函数，还有EAM势函数，不同的物质状态描述用不同的势函数。 模型势函数一旦确定，就可以根据物理学规律求得模拟中的守恒量。 2 第二步 给定初始条件。运动方程的求解需要知道粒子的初始位置和速度，不同的算法要求不同的初始条件。如：verlet算法需要两组坐标来启动计算，一组零时刻的坐标，一组是前进一个时间步的坐标或者一组零时刻的速度值。 一般意思上讲系统的初始条件不可能知道，实际上也不需要精确选择代求系统的初始条件，因为模拟实践足够长时，系统就会忘掉初始条件。当然，合理的初始条件可以加快系统趋于平衡的时间和步伐，获得好的精度。 常用的初始条件可以选择为：令初始位置在差分划分网格的格子上，初始速度则从玻尔兹曼分布随机抽样得到；令初始位置随机的偏离差分划分网格的格子上，初始速度为零；令初始位置随机的偏离差分划分网格的格子上，初始速度也是从玻尔兹曼分布随机抽样得到。 第三步 趋于平衡计算。在边界条件和初始条件给定后就可以解运动方程，进行分子动力学模拟。但这样计算出的系统是不会具有所要求的系统的能量，并且这个状态本身也还不是一个平衡态。 为使得系统平衡，模拟中设计一个趋衡过程，即在这个过程中，我们增加或者从系统中移出能量，直到持续给出确定的能量值。我们称这时的系统已经达到平衡。这段达到平衡的时间成为驰豫时间。 分子动力学中，时间步长的大小选择十分重要，决定了模拟所需要的时间。为了减小误差，步长要小，但小了系统模拟的驰豫时间就长了。因此根据经验选择适当的步长。如，对一个具有几百个氩气Ar分子的体系，lj势函数，发现取h为0.01量级，可以得到很好的相图。这里选择的h是没有量纲的，实际上这样选择的h对应的时间在10-14s的量级呢。如果模拟1000步，系统达到平衡，驰豫时间只有10-11s。 第四步 宏观物理量的计算。实际计算宏观的物理量往往是在模拟的最后揭短进行的。它是沿相空间轨迹求平均来计算得到的（时间平均代替系综平均）
热力学与动力学往年考试整理
判断题： 1.由亚稳相向稳定相转变不需要推动力。? 2.压力可以改变材料的结构，导致材料发生相变。√ 3.对于凝聚态材料，随着压力升高, 熔点提高。√ 4.热力学第三定律指出：在0 K时任何纯物质的熵值等于零。? 5.在高温下各种物质显示相同的比热。√ 6.溶体的性质主要取决于组元间的相互作用参数。√ 7.金属和合金在平衡态下都存在一定数量的空位，因此空位是热力学稳定的缺 陷。√ 8.固溶体中原子定向迁移的驱动力是浓度梯度。? 9.溶体中析出第二相初期，第二相一般与母相保持非共格以降低应变能。? 10.相变过程中如果稳定相的相变驱动力大于亚稳相，一定优先析出。? 1.根据理查德规则，所有纯固体物质具有大致相同的熔化熵。 2.合金的任何结构转变都可以通过应力驱动来实现。 3.在马氏体相变中，界面能和应变能构成正相变的阻力，但也是逆相变的驱动 力。 4.在高温下各种纯单质固体显示相同的等容热容。 5.二元溶体的混合熵只和溶体的成分有关，与组元的种类无关。 6.材料相变形核时，过冷度越大，临界核心尺寸越大。 7.二元合金在扩散时，两组元的扩散系数总是相同。 8.焓具有能量单位，但它不是能量，也不遵守能量守恒定律；但是系统的焓变 可由能量表达。 9.对于凝聚态材料，随着压力升高, 熔点提高, BCC－FCC转变温度也升高。 10.由于马氏体相变属于无扩散切变过程，因此应力可以促发形核和相变。 简答题： 1．一般具有同素异构转变的金属从高温冷却至低温时，其转变具有怎样的体积特征？试根据高温和低温下自由能与温度的关系解释此现象。有一种具有同素异构转变的常用金属和一般金属所具有的普遍规律不同，请指出是那种金属？简要解释其原因？（8分）
流体力学课后答案第九章一元气体动力学基础资料
流体力学课后答案第九章一元气体动力学 基础
一元气体动力学基础 1.若要求22 v p ρ?小于0.05时，对20℃空气限定速度是多少？ 解：根据2 20v P ρ?=42 M 知 4 2M < 0.05?M<0.45，s m kRT C /3432932874.1=??== s m MC v /15334345.0=?== 即对20℃ 空气限定速度为v <153m/s ，可按不压缩处理。 2.有一收缩型喷嘴，已知p 1=140kPa （abs ），p 2=100kPa （abs ），v 1=80m/s ，T 1=293K ，求2-2断面上的速度v 2。 解：因速度较高，气流来不及与外界进行热量交换，且当忽略能量损失时，可按等熵流动处理，应用结果：2v =2121)(2010v T T +-，其中T 1=293K 1ρ=1 1RT p =1.66kg/m 3. k P P 11 212)(ρρ==1.31kg/m 3. T 2=R P 22ρ=266 K 解得：2v =242m/s 3.某一绝热气流的马赫数M =0.8，并已知其滞止压力p 0=5×98100N/m 2，温度t 0=20℃，试求滞止音速c 0，当地音速c ，气流速度v 和气流绝对压强p 各为多少？ 解：T 0=273+20=293K ，C 0=0KRT =343m/s 根据 202 11M K T T -+=知 T=260 K ，s m kRT C /323==，s m MC v /4.258==
100-??? ??=k k T T p p 解得：2/9810028.3m N p ?= 4.有一台风机进口的空气速度为v 1，温度为T 1，出口空气压力为p 2，温度为T 2，出口断面面积为A 2，若输入风机的轴功率为N ，试求风机质量流量G （空气定压比热为c p ）。 解：由工程热力学知识： ???? ??+=22 v h G N ??，其中PA GRT T c h P ==，pA GRT A G v ==ρ ∴?? ????????+-??????+=)2()(2121122222v T c A p GRT T c G N P P 由此可解得G 5.空气在直径为10.16cm 的管道中流动，其质量流量是1kg/s ，滞止温度为38℃，在管路某断面处的静压为41360N/m 2，试求该断面处的马赫数，速度及滞止压强。 解：由G =v ρA ?=RT p ρv=pA GRT ?-+=kRT v k T T 2 0211T =282k 又：202 11M k T T -+= ∴717.0=M s m kRT M MC v /4.241=== ???? ??=-100k k T T p p p 0=58260N/m 2
分子动力学模拟
分子动力学模拟 分子动力学就是一门结合物理,数学与化学的综合技术。分子动力学就是一套分子模拟方法,该方法主要就是依靠牛顿力学来模拟分子体系的运动,以在由分子体系的不同状态构成的系统中抽取样本,从而计算体系的构型积分,并以构型积分的结果为基础进一步计算体系的热力学量与其她宏观性质。 这门技术的发展进程就是: 1980年:恒压条件下的动力学方法(Andersenの方法、Parrinello-Rahman法) 1983年:非平衡态动力学方法(Gillan and Dixon) 1984年:恒温条件下的动力学方法(能势‐フーバーの方法) 1985年:第一原理分子动力学法(→カー?パリネロ法) 1991年:巨正则系综的分子动力学方法(Cagin and Pettit)、 最新的巨正则系综,即为组成系综的系统与一温度为T、化学势为μ的很大的热源、粒子源相接触,此时系统不仅同热源有能量交换,而且可以同粒子源有粒子的交换,最后达到平衡,这种系综称巨正则系综。 进行分子动力学模拟的第一步就是确定起始构型,一个能量较低的起始构型就是进行分子模拟的基础,一般分子的其实构型主要就是来自实验数据或量子化学计算。在确定起始构型之后要赋予构成分子的各个原子速度,这一速度就是根据玻尔兹曼分布随机生成,由于速度的分布符合玻尔兹曼统计,因此在这个阶段,体系的温度就是恒定的。另外,在随机生成各个原子的运动速度之后须进行调整,使得体系总体在各个方向上的动量之与为零,即保证体系没有平动位移。 由上一步确定的分子组建平衡相,在构建平衡相的时候会对构型、温度等参数加以监控。 进入生产相之后体系中的分子与分子中的原子开始根据初始速度运动,可以想象其间会发生吸引、排斥乃至碰撞,这时就根据牛顿力学与预先给定的粒子间相互作用势来对各个例子的运动轨迹进行计算,在这个过程中,体系总能量不变,但分子内部势能与动能不断相互转化,从而体系的温度也不断变化,在整个过程中,体系会遍历势能面上的各个点,计算的样本正就是在这个过程中抽取的。 用抽样所得体系的各个状态计算当时体系的势能,进而计算构型积分。 作用势的选择与动力学计算的关系极为密切,选择不同的作用势,体系的势能面会有不同的形状,动力学计算所得的分子运动与分子内部运动的轨迹也会不同,进而影响到抽样的结果与抽样结果的势能计算,在计算宏观体积与微观成分关系的时候主要采用刚球模型的二体势,计算系统能量,熵等关系时早期多采用Lennard-Jones、morse势等双体势模型,对于金属计算,主要采用morse势,但就是由于通过实验拟合的对势容易导致柯西关系,与实验不符,因此在后来的模拟中有人提出采用EAM等多体势模型,或者采用第一性原理计算结果通过一定的物理方法来拟合二体势函数。但就是对于二体势模型,多体势往往缺乏明确的表达式,参量很多,模拟收敛速度很慢,给应用带来很大困难,因此在一般应用中,通过第一性原理计算结果拟合势函数的L-J,morse等势模型的应用仍非常广泛。 分子动力学计算的基本思想就是赋予分子体系初始运动状态之后,利用分子的自然运动在相空间中抽取样本进行统计计算,时间步长就就是抽样的间隔,因而时间步长的选取对动力学模拟非常重要。太长的时间步长会造成分子间的激烈碰撞,体系数据溢出;太短的时间步长会降低模拟过程搜索相空间的能力,因此一般选取的时间步长为体系各个自由度中最短运动周期的十分之一。但就是通常情况下,体系各自由度中运动周期最短的就是各个化学键的振动,而这种运动对计算某些宏观性质并不产生影响,因此就产生了屏蔽分子内部振动或其她无关运动的约束动力学,约束动力学可以有效地增长分子动力学模拟时间步长,提高搜索相空间的能
热力学与动力学
2007年度《材料热力学与动力学》考试题 姓名：学号： 判断题：（用√或?符号指出对错。20分，每小题2分） 1.热力学第三定律指出：在0 K时任何纯物质的熵值等于零。 2.焓具有能量单位，但它不是能量，也不遵守能量守恒定律；但是系统的焓变可由能量表达。 3.在高温下各种物质显示相同的比热。 4.对于凝聚态材料，随着压力升高, 熔点提高, BCC－FCC转变温度也升高。 5.溶体的性质主要取决于组元间的相互作用参数。 6.亚稳相即使没有外力推动，随着时间的延长，最终会转变为稳定相。 7.金属和合金在平衡态下都会存在一定数量的空位，因此说空位是热力学稳定的缺陷。 8.固溶体中原子定向迁移的驱动力是浓度梯度。 9.溶体中析出第二相初期，第二相一般与母相保持非共格以降低应变能。 10.由于马氏体相变属于无扩散切变过程，因此应力可以促发形核和相变。 简答题：（40分，每小题8分，任选5题，其余题目答出可酌情加分） 1．一般具有同素异构转变的金属从高温冷却至低温时，其转变具有怎样的体积特征？试根据高温和低温下自由能与温度的关系解释此现象。有一种具有同素异构转变的常用金属和一般金属所具有的普遍规律不同，请指出是那种金属？简要解释其原因？ 2．试举出三种二元溶体模型；简要指出各溶体模型的原子相互作用能I AB的特征。 3．试利用给出的a,b两种溶体Gm-X 图中化学势的图解示意图(右图)， 指出两种溶体的扩散特征有什么 不同；那一种固溶体中会发生上坡 扩散。 ( a ) ( b ) 4．向Cu中加入微量的Bi、As合金时 所产生的效果完全不同。加入微量的 Bi会使Cu显著变脆，而电阻没有显
大学物理第九章热力学基础历年考题
第9章热力学基础 一、选择题 1. 对于准静态过程和可逆过程, 有以下说法．其中正确的是 [] (A>准静态过程一定是可逆过程 (B>可逆过程一定是准静态过程 (C>二者都是理想化的过程 (D>二者实质上是热力学中的同一个概念 2. 对于物体的热力学过程, 下列说法中正确的是 [] (A>内能的改变只决定于初、末两个状态, 与所经历的过程无关 (B>摩尔热容量的大小与所经历的过程无关 (C>在物体内, 若单位体积内所含热量越多, 则其温度越高 (D>以上说法都不对 3. 有关热量, 下列说法中正确的是 [](A>热是一种物质 (B>热能是物质系统的状态参量 (C>热量是表征物质系统固有属性的物理量 (D>热传递是改变物质系统内能的一种形式 4. 关于功的下列各说法中, 错误的是 [](A>功是能量变化的一种量度 (B>功是描写系统与外界相互作用的物理量 (C>气体从一个状态到另一个状态, 经历的过程不同, 则对外作的功也不一样 (D>系统具有的能量等于系统对外作的功 5. 理想气体状态方程在不同的过程中有不同的微分表达式, 式表示 [](A>等温过程(B>等压过程 (C>等体过程(D>绝热过程 6. 理想气体状态方程在不同的过程中可以有不同的微分表达式, 式表示 [](A>等温过程(B>等压过程 (C>等体过程(D>绝热过程 7. 理想气体状态方程在不同的过程中可以有不同的微分表达式, 式表示 [](A>等温过程(B>等压过程 (C>等体过程(D>绝热过程 8. 理想气体状态方程在不同的过程中可以有不同的微分表达式,
则式表示 [](A>等温过程(B>等压过程 (C>等体过程(D>任意过程 9. 热力学第一定律表明: [](A>系统对外作的功不可能大于系统从外界吸收的热量 (B>系统内能的增量等于系统从外界吸收的热量 (C>不可能存在这样的循环过程, 在此过程中, 外界对系统所作的功 不等于系统传给外界的热量 (D>热机的效率不可能等于1 10. 对于微小变化的过程, 热力学第一定律为d Q= d E d A．在以下过程中, 这三者同时为正的过程是 [](A>等温膨胀(B>等容膨胀 (C>等压膨胀(D>绝热膨胀 11. 对理想气体的等压压缩过程，下列表述正确的是 [](A> d A＞0, d E＞0, d Q＞0 (B> d A＜0, d E＜0, d Q＜0 (C> d A＜0, d E＞0, d Q＜0 (D> d A = 0, d E = 0, d Q = 0 12. 功的计算式适用于 [](A>理想气体(B>等压过程 (C>准静态过程(D>任何过程 13. 一定量的理想气体从状态出发, 到达另一状态．一次是等温压缩到, 外界作功A；另一次为绝热压缩到, 外界作功W．比较这两个功值的大小是 [](A>A＞W(B>A = W(C>A＜W (D>条件不够，不能比较 14. 1mol理想气体从初态(T1、p1、V1 >等温压缩到体积V2, 外界对气体所作的功为 [](A>(B> (C>(D> 15. 如果W表示气体等温压缩至给定体积所作的功, Q表示在此过程中气体吸收的热量, A表示气体绝热膨胀回到它原有体积所作的功, 则整个过程中气体内能的变化为 [](A>W＋Q－A(B>Q－W－A (C>A－W－Q(D>Q＋A－W
热力学与动力学
热力学与动力学
————————————————————————————————作者：————————————————————————————————日期: ?
2０06年度《材料热力学与动力学》考试题 简答题: 1．一般具有同素异构转变的金属从高温冷却至低温时,其转变具有怎样的体积特征?试根据高温和低温下自由能与温度的关系解释此现象。有一种具有同素异构转变的常用金属和一般金属所具有的普遍规律不同，请指出是那种金属?简要解释其原因?（8分） 2．金属和合金在平衡态下存在一定数量的空位,因此有人说一定数量的空位是金属和合金中的热力学稳定缺陷，此说法是否正确?根据空位数量对自由能及其组成要素(焓和熵)的影响方式，从热力学角度进行简要解释。(8分） 3．试举出三种二元溶体模型；简要指出各溶体模型的原子相互作用能IＡB的特征。 (6分) 4．试利用给出的a,ｂ两种溶体Gm-X图中化学势的图解示意图,指出两种溶体的扩散特征有什么不同；那一种固溶体中会发生上坡扩散。(７分) （a) (b) 5．向Cu中加入微量的Bi、Ａｓ合金时所产生的效果完全不同。加入微量的Ｂi会使Cu显著变脆,而电阻没有显著变化，加入微量的Aｓ并不会使Ｃｕ变脆,但是能显著提高电阻。试根据下面的相图，从溶解度角度对上述现象加以解释。(8分)
6．将固溶体相和晶界相视为两相平衡状态,如果已知上述两相的自由能-成分曲线，指出:采用什么方法或法则来确定两相的平衡成分？一般来说,两相的平衡溶质成分具有怎样的关系?（5分） 7．简要回答什么是耗散结构以及产生耗散结构的必要条件；举出2个自组织现象的实例。(６分） 8．在相变形核阶段，体积自由能、界面能以及应变能中哪些是相变的驱动力?哪些是相变的阻力？试解释：在形核阶段,形核的总自由能为正值,为什么核心能形成呢?以马氏体为例,在核心长大阶段的自由能以及界面能和应变能如何变化?（８分) 9．根据过饱和固溶体中析出第二相时的相平衡关系或者Gibｂｓ-Thomson方程,简要说明第二相粒子粗化过程;从温度对长大速率和对扩散两个方面的影响，简要说明温度对粒子粗化的作用。(７分) 分析计算题： 1．已知纯钛α／β平衡温度为８82?Ｃ，相变焓为14．65 kＪ/ｍoｌ。估算β钛过冷到800?C时，β－Ti转变为α-Ti的相变驱动力（不计上述过冷温度范围对相变的焓变及熵变的影响）。（10分) 2．从过饱和固溶体（α）中析出的第二相通常都是很小的粒子（β)，一般这些小粒子在表面张力的作用下会受到附加压应力的作用，写出附加压应力与表面张力和球形粒子尺寸的关系。以二元溶体为例,用图示的方法简要分析附加压应力对溶体相与析出相界面(α／β）平衡关系的影响。在析出的初期，这小粒子一般与基体保持共格关系,简要分析其原因。(１5分） 3．在2５?C和0.1MＰa下,金刚石和石墨的标准熵分别为２.4 J／mol?Ｋ和５．7 Ｊ/ｍol?Ｋ，标准焓分别为３9５kJ／mｏl和39４kJ/ｍｏl，密度分别为3.５g／cm３和2.3g/cm3, 碳的摩尔质量为１2g。试计算石墨在此条件下转变为金刚石的相变驱动力;试根据自由能与体积和温度的关系(dG=VdP - ＳdT)计算室温下实现石墨－金刚石转变所需临界压力（不计压力对石墨以及压力对金刚石造成的体积改变)。(15分)
一元气体动力学基础学习资料
一元气体动力学基础
一元气体动力学基础 1.若要求22 v p ρ?小于0.05时，对20℃空气限定速度是多少？ 解：根据20v P ρ?=42M 知 42 M < 0.05?M<0.45，s m kRT C /3432932874.1=??== s m MC v /15334345.0=?== 即对20℃ 空气限定速度为v <153m/s ，可按不压缩处理。 2.有一收缩型喷嘴，已知p 1=140kPa （abs ），p 2=100kPa （abs ），v 1=80m/s ，T 1=293K ，求2-2断面上的速度v 2。 解：因速度较高，气流来不及与外界进行热量交换，且当忽略能量损失时，可按等熵流动处理，应用结果：2v =2121)(2010v T T +-，其中T 1=293K 1ρ=11RT p =1.66kg/m 3. k P 1 12 12)(ρρρ==1.31kg/m 3. T 2=R P 22ρ=266 K 解得：2v =242m/s 3.某一绝热气流的马赫数M =0.8，并已知其滞止压力p 0=5×98100N/m 2，温度t 0=20℃，试求滞止音速c 0，当地音速c ，气流速度v 和气流绝对压强p 各为多少？ 解：T 0=273+20=293K ，C 0=0KRT =343m/s
根据 20 2 11M K T T -+=知 T=260 K ，s m kRT C /323== ，s m MC v /4.258== 100-??? ??=k k T T p p 解得：2/9810028.3m N p ?= 4.有一台风机进口的空气速度为v 1，温度为T 1，出口空气压力为p 2，温度为T 2，出口断面面积为A 2，若输入风机的轴功率为N ，试求风机质量流量G （空气定压比热为c p ）。 解：由工程热力学知识： ???? ??+=22 v h G N ??，其中PA GRT T c h P ==，pA GRT A G v ==ρ ∴?? ????????+-??????+=)2()(2121122222v T c A p GRT T c G N P P 由此可解得G 5.空气在直径为10.16cm 的管道中流动，其质量流量是1kg/s ，滞止温度为38℃，在管路某断面处的静压为41360N/m 2，试求该断面处的马赫数，速度及滞止压强。 解：由G =v ρA ?=RT p ρv=pA GRT ?-+=kRT v k T T 2 0211T =282k 又：202 11M k T T -+= ∴717.0=M s m kRT M MC v /4.241===
流体力学课后答案解析第九章一元气体动力学基础
一元气体动力学基础 1.若要求22 v p ρ?小于0.05时，对20℃空气限定速度是多少？ 解：根据2 20v P ρ?=42 M 知 4 2 M < 0.05?M<0.45，s m kRT C /3432932874.1=??== s m MC v /15334345.0=?== 即对20℃ 空气限定速度为v <153m/s ，可按不压缩处理。 2.有一收缩型喷嘴，已知p 1=140kPa （abs ），p 2=100kPa （abs ），v 1=80m/s ，T 1=293K ，求2-2断面上的速度v 2。 解：因速度较高，气流来不及与外界进行热量交换，且当忽略能量损失时，可按等熵流动 处理，应用结果：2v =2121)(2010v T T +-，其中T 1=293K 1ρ=1 1RT p =1.66kg/m 3. k P P 11 212)(ρρ==1.31kg/m 3. T 2=R P 22ρ=266 K 解得：2v =242m/s 3.某一绝热气流的马赫数M =0.8，并已知其滞止压力p 0=5×98100N/m 2，温度t 0=20℃，试求滞止音速c 0，当地音速c ，气流速度v 和气流绝对压强p 各为多少？ 解：T 0=273+20=293K ，C 0=0KRT =343m/s 根据 202 11M K T T -+=知 T=260 K ，s m kRT C /323==，s m MC v /4.258== 100-??? ??=k k T T p p 解得：2/9810028.3m N p ?= 4.有一台风机进口的空气速度为v 1，温度为T 1，出口空气压力为p 2，温度为T 2，出口断面面积为A 2，若输入风机的轴功率为N ，试求风机质量流量G （空气定压比热为c p ）。 解：由工程热力学知识：
vasp做分子动力学
vasp做分子动力学的好处，由于vasp是近些年开发的比较成熟的软件，在做电子scf速度方面有较好的优势。 缺点：可选系综太少。 尽管如此，对于大多数有关分子动力学的任务还是可以胜任的。 主要使用的系综是NVT和NVE。 下面我将对主要参数进行介绍！ 一般做分子动力学的时候都需要较多原子，一般都超过100个。 当原子数多的时候，k点实际就需要较少了。有的时候用一个k点就行，不过这都需要严格的测试。通常超过200个原子的时候，用一个k点，即Gamma点就可以了。 INCAR： EDIFF 一般来说，用1E-4或者1E-5都可以，这个参数只是对第一个离子步的自洽影响大一些，对于长时间的分子动力学的模拟，精度小一点也无所谓，但不能太小。 IBRION=0 分子动力学模拟 IALGO=48 一般用48，对于原子数较多，这个优化方式较好。 NSW=1000 多少个时间步长。 POTIM=3 时间步长,单位fs,通常1到3. ISIF=2 计算外界的压力. NBLOCK= 1 多少个时间步长，写一次CONTCAR，CHG和CHGCAR，PCDAT. KBLOCK=50 NBLOCK*KBLOCK个步长写一次XDATCAR. ISMEAR=-1 费米迪拉克分布. SIGMA =0.05 单位:电子伏 NELMIN=8 一般用6到8,最小的电子scf数.太少的话,收敛的不好. LREAL=A APACO=10 径向分布函数距离,单位是埃. NPACO=200 径向分布函数插的点数. LCHARG=F 尽量不写电荷密度,否则CHG文件太大. TEBEG=300 初始温度. TEEND=300 终态温度。不设的话，等于TEBEG. SMASS -3 NVE ensemble；-1 用来做模拟退火；大于0 NVT 系综。 ///////////////////////////////////////////////////////////////////// ///////////////////////////////////////////////////////////////////// 1)收敛判据的选择 结构弛豫的判据一般有两种选择：能量和力。这两者是相关的，理想情况下，能量收敛到基态，力也应该是收敛到平衡态的。但是数值计算过程上的差异导致以二者为判据的收敛速度差异很大，力收敛速度绝大部分情况下都慢于能量收敛速度。这是因为力的计算是在能量的基础上进行的，能量对坐标的一阶导数得到力。计算量的增大和误差的传递导致力收敛慢。 到底是以能量为收敛判据，还是以力为收敛判据呢？关心能量的人，觉得以能量
【核心知识】蛋白质折叠的热力学和动力学
蛋白质折叠的热力学和动力学 药学院 10489629 苟宝迪 蛋白质是一种生物大分子，基本上是由20种氨基酸以肽键连接成肽链。肽链在空间卷曲折叠成为特定的三维空间结构。有的蛋白质由多条肽链组成，每条肽链称为亚基，亚基之间又有特定的空间关系，称为蛋白质的四级结构。所以蛋白质分子有非常特定的复杂的空间结构。诺贝尔奖得主Anfinsen认为每一种蛋白质分子都有自己特有的氨基酸的组成和排列顺序，由这种氨基酸排列顺序决定它的特定的空间结构。具有完整一级结构的多肽或蛋白质, 只有当其折叠形成正确的三维空间结构才可能具有正常的生物学功能. 如果这些生物大分子的折叠在体内发生了故障, 形成错误的空间结构, 不但将丧失其生物学功能, 甚至会引起疾病.蛋白质异常的三维空间结构可以引发疾病，疯牛病、老年性痴呆症、囊性纤维病变、家族性高胆固醇症、家族性淀粉样蛋白症、某些肿瘤、白内障等等都是“折叠病”。 蛋白质折叠的研究（图1[1]），是生命科学领域的前沿课题之一。不仅具有重大的科学意义，而且在医学和在生物工程领域具有极大的应用价值。 图1 蛋白质折叠的热力学研究 蛋白质折叠的研究，比较狭义的定义就是研究蛋白质特定三维空间结构形成的规律、稳定性和与其生物活性的关系。这里最根本的科学问题就是多肽链的一级结构到底如何决定它的空间结构？X-射线晶体衍射是至今为止研究蛋白质结构最有效的方法, 所能达到的精度是其它任何方法所不能比拟的. 但是, 蛋白质分离纯化技术要求高, 蛋白质晶体难以培养,
晶体结构测定的周期较长, 从而制约了蛋白质工程的进展. 随着近代物理学、数学和分子生物学的发展, 特别是计算机技术的进步, 人们开始用理论计算的方法, 利用计算机来预测蛋白质的结构. 同源模建方法是最常用、最有效的蛋白质结构预测方法. 但是, 利用同源模建方法预测蛋白质结构时, 需用同源蛋白质的已知结构作为模板. 当缺乏这种模板结构时, 预测则很难奏效. 这是该方法的天生缺陷. 是否能从蛋白质序列出发, 直接预测蛋白质的结构? 从理论上最直接地去解决蛋白质的折叠问题，就是根据测得的蛋白质的一级序列预测由Anfinsen原理决定的特定的空间结构。蛋白质氨基酸序列，特别是编码蛋白质的核苷酸序列的测定现在几乎已经成为常规技术，利用分子生物学技术可以从互补DNA（cDNA）序列可以推定氨基酸序列，大大加速了蛋白质一级结构的测定。目前蛋白质数据库中已经存有大约17万个蛋白的一级结构，但是测定了空间结构的蛋白大约只有1．2万个，这中间有许多是很相似的同源蛋白，已经有人根据基因组的数据用统计方法重新估计了蛋白质折叠类型数目大约为1000种。 “蛋白质结构预测”属于理论方面的热力学问题，蛋白质分子结构本身的复杂性决定了结构预测的复杂性。目前结构预测的方法大致可分为两大类。一类是假设蛋白质分子天然构象处于热力学最稳定，能量最低状态，考虑蛋白质分子中所有原子间的相互作用以及蛋白质分子与溶剂之间的相互作用，采用分子力学的能量极小化方法，计算出蛋白质分子的天然空间结构。第二类方法是利用存入蛋白质数据库的数据进行预测相比，基于同源性的重复循环技术非常可靠地灵敏地进行结构预测。找出数据库中已有的蛋白质的空间结构与其一级序列之间的联系总结出一定的规律，逐级从一级序列预测二级结构，再建立可能的三维模型，根据总结出的空间结构与其一级序列之间的规律，排除不合理的模型，再根据能量最低原理得到修正的结构。但是，第一类方法遇到在数学上难以解决的多重极小值问题，而逐级预测又受到二级结构预测精度的限制。 图2[2]为蛋白质折叠研究的漏斗模型。从能量的角度看，漏斗表面上的每一个点代表蛋白质的一种可能的构象，变性状态的蛋白质构象位于漏斗顶面，漏斗最底部的点表示用X-射线单晶衍射或NMR测定的蛋白质天然构象，而漏斗侧面的斜率用来说明蛋白质折叠路径（图3[1]）。 图2
热力学和动力学测验题
化学原理模块课堂测验 1、可逆反应：C(s)+H2O(g)? CO(g)+H2(g) Δr H mΘ>0。下列说法你认为对否？为什么？ （1）达平衡时各反应物和生成物的分压一定相等； （2）改变生成物的分压，使Q 热力学与动力学的关系 黄金金 指导教师：陶中东 摘要：反应趋势与反应速率之间的关系是化学研究中的一个十分重要的问题,本文根据物理化学的基本原理,建立了热力学判据和动力学判据相结合的普适判据，并根据普适判据，从理论上阐述了反应趋势和反应速率二者之间的内在联系，同时论述了化学热力学与化学动力学的紧密联系，有助于从物质的本性、微观结构去进一步认识化学反应的实质。 关键词：化学热力学，化学动力学，普适判据，反应进度，反应速率，活化能，活化熵，化学亲和势 The Relationship between Chemical Kinetics and Thermodynamics Abstract: According to the principle of physical chemistry, a general criterion, which includes thermodynamic criterion of reaction trend and kinetic criterion of reaction rate, has been proposed suggested in this paper. The internal relations between reaction trend and reaction rate have been elaborated based on the general criterion. Meanwhile the paper expounded the close relationship between chemical kinetics and thermodynamics it is instrumental to us in further cognizing chemical reaction essence from material nature and structure. Key words: general criterion; chemical thermodynamics; chemical kinetics; reaction extent; reaction rate; activated energy; activated entropy; chemical affinity 对于一个指定反应条件和反应物种的化学反应来说，研究者普遍关心的问题是该化学反应的趋势与速率。有许多实例表明化学反应的趋势与速率二者之间存在正相关性，即反应的趋势大且反应的速率快。“化学动力学与化学热力学是相辅相成的，动力学的研究必须以热力学的结果(肯定反应有可能发生)为前提条件，而热力学只有与动力学相结合才能全面解决化学反应的实际问题”。“我们用以控制化学过程的方法主要是改变:(1)温度，(2)压力，(3)反应物的比例，(4)催热力学与动力学的关系