《科学计算选讲》结课报告 - 天津大学研究生e-learning平台.docVIP

《科学计算选讲》结课报告 利用ANSYS分析Pb靶在激光作用下的温度场分布 学院 材料学院 专业 材料学（金属） 年级 2013级直博 姓名 邹成琴 学号 1013208035 2014年5月 利用ANSYS分析Pb靶在激光作用下的温度场分布 邹成琴 （天津大学 材料科学与工程学院 1013208035） 摘要：本文利用ANSYS软件简单模拟了液氮中的Pb靶受激光烧蚀后的温度场分布情况，发现由于液氮的极低温度，使得靶材表面的温度传递受到了抑制，沿径向急剧下降。而在靶材纵向方向，温度也是由高到低下降较快，高于Pb沸点的层厚较大。若采用更精细的网格划分方法以及瞬态热力学分析，将会获得更准确的结果，这有待进一步探究。 关键词：数值分析；ANSYS；激光；液氮；Pb靶；温度场分布 一、背景介绍 激光烧蚀法的典型特征是，利用脉冲激光作用固体靶材，使靶材在脉冲的瞬间（10 ms）被加热到熔化甚至气化，从而产生由靶材原子、自由电子、离子或团簇组成的等离子体，这一高温高压的等离子体在随后的膨胀、塌陷过程中，会与环境气体原子（或液体分子）碰撞、混合、反应，最终在其冷却的过程中，形核、生长形成纳米颗粒（即脉冲激光烧蚀法制备纳米颗粒）或是在基底上沉积一层薄膜（也就是脉冲激光沉积薄膜技术）[1]。 根据制备环境的不同，激光烧蚀法可以分为气相烧蚀和液相烧蚀。两者在等离子体的膨胀和凝聚动力学方面有差异，前者适用于薄膜材料的制备，而后者更适用制备纳米颗粒。另外，由于激光液相烧蚀法可以制备分散的、表面修饰的纳米颗粒，正好满足了人们对分散性好、表面功能化的纳米材料的需求，因此受到了很多研究者的青睐[2-5]。有关合成机理，大家都是众说纷纭，下图1展示的就是其中两种。 图1 (a)等离子体混合—凝聚—形核长大模型示意图；(b)等离子体粉碎液态金属模型示意图 利用激光烧蚀单一金属靶材，可以制备各种金属的纳米颗粒。比如激光作用水中的铂靶制备 Pt 纳米颗粒（图2）[6]；改变气氛或溶液的种类，会形成各种化合物纳米颗粒，比如在水中激光烧蚀不同金属靶得到 ZnO、TiO2、FeO、 Mn3O4 等纳米颗粒[4,7-9]。产生的纳米颗粒分散较好，尺寸较小（30 nm）。 图2 波长1064 nm，单脉冲能量密度14 J/cm2的纳秒脉冲激光作用水中铂靶得到的Pt纳米颗粒的TEM照片[6] 研究通过激光烧蚀法合成纳米材料，必不可少的要涉及到计算作用靶材表面的温度场分布。在前人的工作中，有许多人采用不同的数值算法，用MATLAB、ANSYS等软件对其开展了广泛的研究。比如，为探究毫秒脉冲强激光辐照下，半导体和金属材料的熔融和气化等热损伤过程。赵方东[10]以经典傅立叶热传导方程为基础，结合不同的边界条件，利用显热容法和显式有限差分格式，编制了C语言数值计算程序，详细计算和分析了铝、银、铜、钛等金属材料和硅、锗等半导体材料在熔融和气化等相变前后的空间温度场分布，融化后液体深度和融化速度，得到了相应材料在ms脉冲激光辐照条件下的熔融损伤和汽化损伤的阈值。又比如，赵志培[11]利用MATLAB软件，对辐射-对流耦合边界条件下毫秒激光烧蚀金属靶材的温度场进行了数值计算，并利用伽辽金法和变分原理计算了温升引起的热应力，下图3展示的就是他用软件模拟出表面中心点下0.1毫米位置处温度的分布曲线。 图3 表面中心点下0.1毫米位置处温度的分布曲线 二、提出问题 在我的课题中，也采用激光烧蚀法在不同液体介质环境中合成Pb纳米材料，Pb靶被放置在盛有不同溶剂的烧杯中，用1064nm的平行纳秒激光辐照靶材，就可获得不同的纳米结构材料，实验装置如图4所示。 图4 利用激光烧蚀法在液体介质中合成不同Pb纳米结构材料 通过前期的一些探索，获得了一些实验现象：环己烷中，激光烧蚀Pb靶会获得较多颗粒状产物，但由于未加入任何表面活性剂，所以这些颗粒发生粘连并不分散。另一方面，液氮中，激光烧蚀Pb靶得到了层片状产物，这令人匪夷所思，因为在液体介质中利用激光烧蚀金属靶材的方法一般获得的是纳米颗粒。环己烷和液氮中的TEM表征结果如图5所示。这个对比说明各个液体介质环境会对合成过程产生不同的影响，当然对于环己烷和液氮来说，这主要体现在它们的温度上，一个为常温非极性溶剂，而另一个则为-196℃的极低温非极性液体。在激光烧蚀过程中，无论是常温还是-196℃的温度，都已经是低温了，所以环己烷和液氮会对靶材溅射出来的物质进行淬火，即一下子从高温降到常温或低温，这能够保留高温的某些性质或者说状态。既然两种溶剂有这么大的差别，我们就很好奇在液氮中，Pb靶表面的温度分布是怎样的，是不是真的起到淬火效果了，冷却的深度和宽度大概是多少等等，这