板料激光弯曲过程的数值模拟

采用MSC.Marc非线性有限元软件,对板料激光弯曲过程进行了数值模拟.试验与数值模拟的相互印证表明:数值模拟的模型正确,计算过程可靠,计算结果可为板料激光弯曲工艺参数的选择提供依据.通过模拟多道次直线照射的过程可知,变形效果受照射路径次序的影响很大,应选择从固定端向自由端方向交替的双向照射路径,以保证板料在宽度方向上变形的均匀.模拟结果可为带筋结构壁板激光弯曲成形提供实用方案.     

泡沫铜/低熔点合金复合相变材料凝固放热分析

为研究泡沫铜/低熔点合金(LMPA)复合相变材料在间歇放热工作环境下恢复至初始状态的能力及不同孔隙率泡沫铜的添加对其凝固放热过程的影响，通过数值模拟对比分析了47合金、正二十三烷与泡沫铜复合前后的凝固放热过程，并调节泡沫铜/47合金复合材料孔隙率计算模拟芯片温度在凝固放热过程中温度随时间变化曲线。结果表明:泡沫铜的添加对2类材料凝固过程均有促进作用，模拟芯片恢复至目标温度所需时间分别被缩短6.6%和47.7%。因体积潜热值的差距，泡沫铜/47合金凝固时需放出更多热量，恢复至目标温度的时间较长，是正二十三烷复合相变材料的1.52倍。随着孔隙率的增大，复合相变材料恢复至室温状态所用时长变化不大，考虑到孔隙率对相变热控过程中的影响，实际使用时应综合考虑。      

煤粉燃烧过程的数值模拟

采用计算流体动力学软件对煤粉实际燃烧过程进行数值计算并结合其热态试验数据进行对比分析,已成为验证数学模型和指导工程实践的一种重要研究手段.应用商业软件FLUENT对煤粉燃烧及流场进行了数值模拟分析,采用混合分数/概率密度函数法模拟湍流燃烧,用P-1辐射模型开展辐射传热模拟,利用拉格朗日/欧拉法处理气固两相间的动量、质量和能量交换,对挥发份的析出采用双速率竞争模型,采用动力/扩散反应速率模型模拟煤粉颗粒的表面燃烧.并对模拟结果进行分析与对比,计算结果与实际燃烧过程有较好的一致性.数值模拟计算结果揭示了挥发分释放与焦炭燃烧的过程,为改善和优化煤粉的燃烧提供了重要的参考依据.     

薄板激光弯曲温度场的数值模拟与校验

采用MSC.Marc非线性有限元软件,对薄板激光弯曲过程中的温度场进行了数值模拟 ,得到了引起薄板激光弯曲变形的温度场的变化规律.通过温度传感器测量与激光扫描线相 对应的薄板下表面温度变化的规律来验证数值模拟的结果.模拟值与实测值基本吻合,表明 数值模拟结果可作为激光加工工艺参数选择的依据.     

模拟三维重力波传播过程的并行数值模式

从三维大气运动的基本控制方程出发, 基于MPI消息传递接口和区域分解的思想, 建立了模拟三维可压缩大气中重力波传播过程的并行数值模式. 在对垂直方向进行区域分解的基础上, 针对跳点网格的特点, 以垂直速度主格点线为子区域的下边界, 状态变量主格点线为子区域的上边界进行区域分解. 利用MPI消息传递接口来传递计算各子区域中变量时所要用到的相邻子区域上相应变量, 从而顺利地完成对整个计算域的并行计算. 根据线性重力波理论, 通过模拟小振幅重力波的传播过程对所建模式进行了验证. 结果表明, 模式可以很好地模拟小振幅重力波在三维空间中的传播过程, 模拟的重力波振幅随着传播高度的增加以指数形式增长, 对能量传播路径以及有效扰动位能和扰动动能的模拟结果均与理论预测的结果吻合很好. 另外, 随着进程个数的增加, 完成相同的计算量所用的计算时间也显著减少. 这表明本文建立的并行数值模式不但能够很好的模拟重力波的传播过程, 而且能够有效地节约计算时间.      

复合材料结构声疲劳的一种数值模拟方法

复合材料结构声疲劳问题是许多工程技术人员关心的课题,经常采用数值方法对其进行分析.利用边界元法对其动力响应情况进行了数值模拟,并发展了一种可以进行复合材料结构声疲劳完整数值模拟的方法.所得结果可以为预测复合材料结构的声疲劳寿命提供参考.      

大气波动传播问题的一种无反射数值边界格式

根据无反射边界条件的基本原则,考虑到高阶的数值边界格式可以减小边界误差和虚假反射,利用曲线拟合中最小二乘法的思想,提出了一种高阶的光滑拟合外推边界格式(SFEBS).因为大气运动的控制方程简化后可以归结为一个对流方程,所以作为边界格式之间的比较和检验,在一维情况下用对流方程和数值模拟中具有代表性的波包和激波作为算例,将其与传统的基于Taylor展开思想构造的边界条件(TEBS)进行了比较.计算结果表明,在高阶情况下,用SFEBS计算波动传播问题的虚假反射约为用同阶TEBS计算的1/6,说明高阶的SFEBS可以作为一种很好的无返射数值边界格式.为今后模拟大气波动提供了一种很好的数值边界格式.      

原子氧和紫外线对聚酰亚胺综合作用数值模拟

在建立数学物理模型的基础上,对低地球轨道环境和地面试验环境下有无保护涂层的聚酰亚胺所受原子氧冲蚀及紫外线的综合作用进行了数值模拟,获得了具有工程应用价值的计算结果,并讨论了数学物理模型中各参数对基蚀曲线形状的影响.从数值模拟结果与美国太空试验结果的比较可以看出,得出的数值模拟的结果是正确的,对航天器设计具有重要的指导意义.     

相变储热/辐射器式热沉的传热特性研究

对一种用于舱外活动的相变储热/辐射器式热沉进行了实验研究和数值模拟.设计了相变储热器的实验件和实验方案,进行了储热器的冷却实验.建立了具有导热与对流换热耦合边界的相变储热器数学模型,模型中考虑了自然对流及空穴的影响.编制了相变储热器的性能分析计算程序.通过实验值与计算值的比较分析,验证了程序的可靠性.该模型及计算程序为相变储热器的设计和优化提供了有利的分析工具.     

复杂薄壁零件板多级充液成形及过程数值模拟

板充液成形已被证明是制造复杂零件一种有效的方法.采用板多级充液成形方法对难成形复杂薄壁结构的飞行器零件进行了研究,利用有限元数值模拟的方法,分析了预成形凸模形状、预成形深度,以及液室压力等关键工艺参数对成形的影响.对成形中出现的诸如破裂、起皱等失效形式进行了探讨,给出了终成形时实现成功成形的最高液室压力区域范围.提出了控制成形质量的措施,优化了工艺参数,并根据数值模拟的结果进行了试验验证.结果表明,模拟结果和试验结果达到了较好的一致.      

喷油器瞬态两相流动的建模方法

为了研究柴油机喷油器内部的气穴现象,采用一种新的方法建模并进行瞬态两相流动的数值模拟. 首先求解针阀开启时刻的稳态流场,然后导入自编的程序进行瞬态流场的计算. 自编程序是基于Fluent提供的用户自定义函数,用于控制模型边界条件,分析针阀受力情况和计算其运动速度. 初始计算条件来自稳态流场,每一步瞬态计算都依赖之前流场的计算结果. 这种建模方法保证了数值模拟的连贯性和真实性,确定喷油器初始条件就可以完整地计算喷油过程,得到每一时刻喷油器内部的压力分布和两相流分布等数据. 研究结果表明,这种新的建模方法是可行和有效的,其数值模拟结果和实验结果相吻合,并能揭示出喷油器喷嘴入口处气穴现象的产生、发展和消失过程,是进一步深入开展柴油喷油系统两相流动特性研究的一种新的途径.     

人体温度场的边界元分析

针对人体温度场模型复杂,计算难度很大,需要简化计算,采用了边界元法.利用边界元法是以边界点上未知函数满足积分方程为基础进行近似计算的这一特点,对人体温度场模型进行了三维计算,并给出结果,与其它方法和实验比较,吻合较好.     

基于地形高程数据的复杂地形风场建模方法

提出了一种以地形高程数据为基础,进行三维复杂地形建模及其风场分布求解的数值仿真方法.该方法利用曲面样条插值技术建立复杂地形的三维数字模型,并将地形曲面离散为若干个涡环网格,结合位势流理论与地面边界条件计算得到复杂地形影响下的三维风场.数值计算结果显示,相比基于球绕流理论的工程模拟方法,该方法可对复杂地形做出更为真实的描述,并给出更为合理的风场分布形式.同时该方法建模过程简单,便于实际操作,适于飞机飞行性能评估或飞行模拟器环境建模的工程应用.      

基于CUDA的超声二维声场EFIT仿真

随着图形处理器（GPU）的快速发展,基于计算设备统一构架（CUDA）可以方便地将并行计算技术应用于超声声场数值仿真计算,极大地提升计算效率。阐述了弹性动力学有限积分算法（EFIT）的原理,在采用CPU实现带吸收边界的钢材料二维点源激励声场仿真的基础上,基于GPU实现了仿真模型的并行计算,介绍了GPU程序的设计流程和参数优化方法,包括纹理内存使用、吸收边界优化和数据传输优化。对比了相同条件下CPU和GPU仿真计算的耗时和平均计算效率,定量分析了GPU对于EFIT模型效率的提升。比对结果表明,EFIT具有良好的并行计算条件,采用并行计算方法能够有效提升模型计算速度,对于复杂声场仿真应用具有广阔的应用前景。      

微重力下相变储能单元融化过程数值模拟

为探究微重力环境中,通过肋片强化了传热的相变储能单元中相变材料融化过程,通过数值模拟方法探究了微重力作用时相变材料融化过程中传热特性。通过地面实验与重力作用下数值模拟结果对比验证数值模拟方法的准确性,对比重力和微重力作用2种情况下数值模拟结果以揭示微重力环境中相变材料融化过程的特性。结果表明,当相变储能单元受微重力作用时,相变材料融化速率明显下降,热量主要通过热传导传递,融化的相变材料从顶端膨胀溢出向空间扩散,局部低温区域在相变储能单元中上部。      

临近空间环境下封闭方腔内耦合换热特性

以临近空间浮空器载荷舱为应用背景,对复杂热边界条件下含热源的三维封闭方腔内自然对流、表面辐射和导热的耦合问题进行了数值模拟。综合考虑对流换热、长波辐射、太阳辐射等因素的影响,建立了临近空间热环境模型。通过Fluent软件用户自定义函数（UDF）引入外部非定常的辐射-对流耦合热边界条件,对腔内换热特性的昼夜变化进行研究,并分析了腔壁厚度、发射率和导热系数对其的影响。数值结果表明,腔内平均温度昼夜变化很小,约为12.9 K,但温度场分布随太阳方位变化而变化;腔内对流换热较弱,同一时刻最大温差约为71.3 K;腔壁热阻和发射率增加会削弱自然对流的强度。      

轴流压气机尾流撞击效应机理的数值模拟

通过精确给定转子进口边界条件,以模拟上游静子尾迹对转子流场的影响,这样利用单排的计算结果模拟转/静干涉引起的非定常流动对流场气动参数的影响,因为相邻叶片排的叶片数不同,从而避免传统方法中常用的叶片数约化存在的难以正确模拟非定常激励频率的问题,对计算时间和计算能力的要求也可减少.针对一个低速转子的数值模拟结果显示,如果能合理地组织来流尾迹的扫过频率和尾迹亏损的幅度,对转子气动性能的改善是很显著的,从而为风扇/压气机中控制和利用非定常流动提供了一个可行的方向.     

球坐标系六片网格下三维定态行星际太阳风模拟

采用二阶MacCormack差分格式, 利用稳态的磁流体(MHD)方程组在球坐标系六片网格下模拟研究了行星际太阳风. 六片网格系统能有效避免极区奇性和网格收敛性. 迭代按径向方向推进求解, 很大程度上减少了计算量, 节约了计算时间. 内边界条件根据太阳与行星际观测确定, 比较测试了5种内边界条件, 模拟给出了1922卡林顿周的背景太阳风结构. 几种内边界条件所得模拟结果与行星际观测基本吻合. 太阳风速度采用McGregor 等的经验公式给出, 磁场由水平电流片(HCCS)模型得到, 密度和温度分别根据动量守恒和气压守恒得到, 研究表明采用这样的边界条件模拟结果最佳.      

空间飞行体与压缩区内等离子体非稳态相互作用的数值模拟

在压缩区内，飞行体天线辐射出的高频调制场影响着等离子体中的电荷分布，而电荷分布的扰动又反过来影响场量，采用FTCS有限差分格式方法，利用二维三分量轴对称，对空间飞行体与压缩区内等离子体非稳态相互作用过程进行数值模拟，得到了电荷密度扰动与电场的变化情况，计算结果表明，这种非线性相互作用可以引起场的塌缩并出现密度空洞，通过对密度空洞变化的探测，可对隐身飞行体进行跟踪。     

太阳表面水平运动驱动的磁圈与开放磁场重联的模拟研究

太阳大气的诸多观测事件(如耀斑、喷流等)均被归因于磁重联产生的能量转换. 近年来, 关于太阳风起源, 有研究提出了磁重联使闭合磁圈开放为太阳风供应物质的新模式. 在该模式中, 闭合磁圈被光球超米粒组织对流携带, 向超米粒边界运动, 与位于边界的开放磁场相碰撞进而发生磁重联. 该模式中磁重联的驱动及其效应是本文的研究目标. 磁流体力学(MHD)数值模拟是研究太阳大气磁重联物理过程的重要途径. 本文建立了一个二维MHD数值模型, 结合太阳大气温度和密度的分层分布, 在超米粒组织尺度上模拟了水平流动驱动的闭合磁圈与开放磁场的重联过程. 通过对模拟结果的定量分析, 认为磁重联确实能够将闭合磁圈的物质释放, 进而供应给新的开放磁结构并产生向上流动. 该结果为进一步模拟研究太阳风初始外流奠定了基础.