移动加热器法生长碲锌镉晶体的组分输运与界面形貌研究

移动加热器法(THM)生长碲锌镉晶体时,界面稳定性对晶体生长的质量有很大影响。本文基于多物理场有限元仿真软件Comsol建立了THM生长碲锌镉晶体的数值模拟模型,讨论了Te边界层与组分过冷区之间的关系,对不同生长阶段的物理场、Te边界层与组分过冷区进行仿真研究,最后讨论了微重力对物理场分布的影响,并对比了微重力与正常重力下的生长界面形貌。模拟结果表明,Te边界层与组分过冷区的分布趋势是一致的,在不同生长阶段,流场中次生涡旋的位置会发生移动,从而导致生长界面的形貌随着生长的进行发生变化,同时微重力条件下形成的生长界面形貌最有利于单晶生长。因此,在晶体生长的中前期,对次生涡旋位置的控制和对组分过冷的削弱,是THM生长高质量晶体的有效方案。     

《数值流形法》序

得知郑宏教授的《数值流形法》即将脱稿,激动之情难以言表.我完全同意郑教授将流形法的特点归结为“升阶”和“切割”.受郑教授的委托,我愿在此将建立流形法的历程做一简单的回忆.其实早在1985年———也就是在我完成不连续变形分析法(DDA)之前,“升阶”在我心中已完全成熟.“升阶”似乎能取得一致收敛,弥补有限元只能按能量收敛的缺憾.我没有急于发表是因为一直没能想到令我满意的解决“不连续”的办法.起初曾想过在DDA块体上引入高阶多项式,或在块体上再划网格,也就是用分片多项式来提高DDA块体的变形精度.但我对这些方案都不甚满意,因为实际计算时多项式的阶数会受到限制,引入网格又会招致有限元既有的毛病——解的精度受控于网格质量.     

竖直恒温面自然对流边界层不稳定性研究

在无扰动、随机式扰动以及正弦式扰动下,通过对竖直恒温面处状态Ra为1.328×10^9、Pr为6.24的自然对流进行模拟,探索了热边界层的不稳定性和共振强化自然对流换热。结果表明:(1)竖直自然对流边界层上游位置的随机式扰动对热边界层的影响主要体现在稳定阶段;(2)该状态下的竖直自然对流边界层的特征频率为15 067,且相比于无扰动状态,频率为15 067的正弦式扰动能在竖直恒温面处提高5.15%的换热量;(3)在竖直自然对流边界层上游位置加入特征频率的正弦式扰动,竖直恒温面处的局部努塞尔数Nu均出现明显波动,且波动随着边界层高度的增加而增大。     

一类具有边界层性质的二次奇摄动边值问题

研究了一类具有边界层性质的二次奇摄动边值问题.在相对较弱的条件下,用合成展开法构造出该问题的形式近似式,并应用改进的Harten不动点定理和逆算子定理证明解的存在性及其渐近性质.最后,将所研究的问题和结论推广到更一般的高次情形.     

HyperFLOW软件非结构网格亚跨声速湍流模拟的验证与确认

计算流体力学(computational fluid dynamics,CFD)数值模拟在航空航天等领域发挥越来越重要的作用,然而CFD数值模拟结果的可信度仍然需要通过不断地验证与确认来提高.本文给出了从制造解精度测试、简单到复杂外形湍流模拟网格收敛性研究等三个方面开展CFD软件验证与确认的方法,并对自主研发的CFD软件平台HyperFLOW在非结构网格上模拟亚跨声速湍流问题的能力进行了验证与确认.首先通过基于Euler方程和标量扩散方程的制造解精度测试,分别验证了HyperFLOW在非结构网格上对Euler方程和黏性项的求解精度,结果表明其能够在任意非结构网格上达到设计的二阶精度. 其次,通过NASATurbulence Modeling Resource中的湍流平板、二维翼型近尾迹流动、二维Bump等几个典型的亚声速湍流算例的网格收敛性研究,量化考察了数值结果的观测精度阶和网格收敛性指数,并与国外知名CFD解算器CFL3D,FUN3D的计算结果进行了对比,验证了HyperFLOW对简单湍流问题的模拟能力,且具有良好的网格收敛性和计算精度(阶). 最后,通过NASA CommonResearchModel标模定升力系数的网格收敛性研究和升阻极曲线预测,验证了软件在复杂外形亚跨声速湍流流动数值模拟中也具有良好的可信度.      

旋转射弹高速倾斜入水多相流场与弹道数值模拟

基于VOF多相流模型和有限体积法求解水、汽、气多相流动的RANS方程，结合重叠网格技术和six DOF算法对某一型号舰载射弹倾斜入水过程进行数值模拟研究。首先基于该方法研究了射弹旋转效应对射弹运动特性及流体动力特性的影响，然后对不同入水角下倾斜入水过程进行分析，得到不同倾角下旋转射弹入水空泡形态发展规律、弹体运动特征及流体动力特性变化规律。研究结果表明:射弹的旋转有利于弹体在初始对称面内的弹道稳定性，但会降低弹体侧向稳定性，使射弹受到的阻力系数、俯仰力矩系数变小；入水角越小，形成的空泡越不对称，由射弹运动状态的改变引起的空泡形态变化越明显，在超空泡航行阶段，弹体运动较稳定，不同角度下流体动力系数差别很小，当弹体下表面刺破空泡壁沾湿时，弹体运动状态发生较大变化，流体动力系数迅速增大，此时入水角度过小，弹体容易失稳；弹体的沾湿对空泡形态、弹体运动稳定性和流体动力特性有着重要的影响。     

考虑损伤的修正梯度弹性理论

梯度弹性理论在描述材料微结构起主导作用的力学行为时具有显著优势,将其与损伤理论相结合,可在材料破坏研究中考虑微结构的影响.基于修正梯度弹性理论,将应变张量、应变梯度张量和损伤变量作为Helmholtz自由能函数的状态变量,并在自然状态附近对自由能函数作Taylor展开,进而由热力学基本定律,推导出修正梯度弹性损伤理论本构方程的一般形式.编制有限元程序,模拟土样损伤局部化带的发展演化过程.结果表明,修正梯度弹性损伤理论消除了网格依赖性;损伤局部化带不是与损伤同时发生,而是在损伤发展到一定程度后再逐渐显现出来.     

自由来流湍流与三维壁面局部粗糙诱导平板边界层不稳定T-S波的数值研究

采用直接数值模拟（DNS）方法,研究了在自由来流湍流与三维壁面局部粗糙作用下平板边界层内诱导产生不稳定T-S波的物理问题.数值结果可知,在平板边界层内发现了二维和三维T-S波组成的波包空间序列以及求得了波包向前传播的群速度大小,从而证明了自由来流湍流与三维壁面局部粗糙作用是激励平板边界层内诱导产生不稳定T-S波的一种机制.随后,建立了平板边界层内被激发的二维和三维T S波的初始幅值与自由来流湍流度,三维壁面局部粗糙的流向长度、展向宽度及法向高度之间的关系.这一问题的深入研究,进一步完善了流动稳定性与湍流理论.     

轴向移动系统的参数振动问题研究进展

对近二十年来轴向移动系统(移动弦，移动梁和移动带等)的参数振动研究进展进行了详细的评述，特别关注了轴向张紧力和移动速度随时间改变时轴向移动系统的参数振动特性和稳定性等问题。文章首先讨论了所研究问题的控制方程。然后详细说明了目前研究中人们较为关注的几个重点问题，如参数激励的形式，求解方法和所研究的问题等。接着在其后的两节中，分别评述了在张紧力和移动速度随时间变化时，轴向移动弦和轴向移动梁的振动问题近年来的研究进展，详细、深入讨论了模型的类型、张紧力和轴向移动速度随时间变化的形式以及在研究中使用的解题方法和系统的振动特性(振动响应、固有频率和动态稳定性)等；最后给出了在此领域今后研究中应关注的问题。