面向深空探测的30cm离子推力器寿命预估

为评价30cm离子推力器的寿命,提出了有限寿命考核结合栅极仿真模型的推力器寿命预估方法,利用粒子-蒙特卡洛（PIC-MCC）方法建立了栅孔溅射腐蚀模型,开展了30cm离子推力器寿命预估研究,分析了栅孔刻蚀速率、单孔电场及离子引出特性,给出了30cm推力器寿命预估值。结果表明,栅孔直径仿真与实测值一致性较好,误差在20%以内,基于每个寿命小节栅孔实测值对模型的修正是有效的;根据仿真结果,10000h寿命考核后,减速栅最先失效,30cm离子推力器在最大功率（3kW）工况下预估寿命为37540h,能够满足小天体探测任务对电推进系统的长寿命需求。     

跨声速机翼抖振初始迎角N-S方程定常计算分析

 应用代数法和椭圆型方程优化相结合的方法生成翼身组合体块结构网格。采用有限体积空间离散法和五步Runge-Kutta显式时间推进法求解N-S方程。基于雷诺平均N-S（RANS）方程,选用Spalart-Allmaras（SA）一方程模型模拟紊流。根据N-S方程定常计算结果,采用升力曲线、俯仰力矩曲线、后缘压力发散、跨声速激波位置以及机翼表面极限流线等几种判据,对跨声速机翼的抖振初始迎角进行了合理的预测分析。     

在非结构动网格上追踪多介质界面的ALE方法

在处理带有介质分界面的多介质流动问题时,把介质边界当作一种特殊的边界,并且在该边界上的点都有着双重的定义,分别对应边界两侧的不同介质状态,通过利用HLLC方法求解任意拉格朗日-欧拉方程(ALE),以及运用界面跟踪和动网格方法分别对一维和二维多介质流进行了数值模拟,并得到了满意的结果.     

柔性扑翼非定常流场的数值计算方法

提出一种将Delaunay图映射网格变形技术和非结构嵌套网格方法结合使用的策略,解决网格变形和嵌套网格单独用于柔性扑翼流场计算时需要网格再生的问题。该方法为嵌套网格中的每个嵌于背景网格的贴体非结构网格生成Delaunay背景图;每个时间步,根据扑翼的运动和变形规律移动背景图,再根据网格点和背景图的映射关系移动网格点,之后自动完成嵌套边界的定义和插值关系的建立。为方便嵌套关系的建立,嵌套网格进行分层管理。也研究了一种内存消耗少、效率较高的搜索算法,以及格心格式和格点格式统一的边界拓宽算法。非定常可压缩Navier-Stokes方程在非结构的动态网格上用有限体积法离散,并用预处理的双时间步推进、隐式LU-SGS迭代求解。几个扑翼算例的结果表明,该方法充分利用了Delaunay图映射网格变形方法的高效率,同时也发挥了嵌套网格处理大幅运动的优势;用于既有整体大幅扑动又有局部小变形的柔性扑翼流场计算,可取得令人满意的精度和效率。     

一种TVD方法在三维非结构网格中的应用

以三维非结构网格的显式有限体积法为基础,采用了一种TVD方法求解三维Euler方程,使用Roe通量来计算网格单元边界处的守恒量通量.为了验证方法的可行性,用该方法模拟三维爆炸问题,得出的结论是我们的方法可行,稳定且有效.     

附属小圆柱对主圆柱绕流影响的数值模拟

从小马赫数下的可压缩Navier-Stokes方程出发,采用覆盖网格分区算法,对在主圆柱尾流场的适当位置放置附属小圆柱时的流动特性及其对主圆柱绕流问题的影响进行了详细的数值模拟研究.根据模拟结果,主要探讨在小雷诺数范围内,附属小圆柱位置、个数、大小以及雷诺数变化等因素对主圆柱绕流的流场结构及其非定常演化过程的影响.计算结果表明,在小雷诺数范围内,无论附属小圆柱的个数如何,都存在最佳抑制区域.当附属小圆柱放在这个区域时,涡街被完全抑制,整个流场达到准定常状态;同时,阻力系数也显著下降.在一定程度上,随着小圆柱直径的增加,主圆柱涡脱落更易被抑制.另外,相对于单个附加小圆柱而言,使用两个附加小圆柱能使主圆柱涡脱落的抑制提高到更高的雷诺数.      

多区域多重网格法在直升机粒子分离器中的应用

首次将多区域多重网格法用在一般曲线坐标系下交错网格布局的 SIMPLEC算法中,在直升机粒子分离器这类独特结构中的成功应用,突破多重网格方法在多连通域中运用的局限性,使复杂问题变得简单。结果表明,它不仅能加快流场求解的收敛速度,而且还能提高计算精度。     

微米尺度格栅标准样片的质量参数评价;；

论述了以CD-SEM作为核心仪器评价微米尺度格栅标准样片质量参数的方法。以某标称线距尺寸为3μm的不同材料、不同格栅高度的标准样片为例进行了各项质量参数的测量，并对实验数据进行相关因素的比对分析。通过分析，为制作优质的微米格栅标准样片提供了一定的理论参考。分析结果表明本文所述方法适用于微米尺度格栅标准样片的质量参数评价。     

非结构CFD软件MPI+OpenMP混合并行及超大规模非定常并行计算的应用

常规工程应用中，非定常数值模拟（如多体分离）的计算量十分巨大，如果为了达到更高的计算精度，加密网格或者采用高精度方法将会使得计算量进一步增大，导致非定常数值模拟在CFD工程应用中成为十分耗时和昂贵的工作，因此，提高非定常数值模拟的可扩展性和计算效率十分必要。为充分发挥既有分布内存又有共享内存的多核处理器的性能和效率优势，对作者团队开发的非结构网格二阶精度有限体积CFD软件（HyperFLOW）进行了混合并行改造，在计算节点间采用MPI消息传递机制，在节点内采用OpenMP共享内存的MPI+OpenMP混合并行策略。首先分别实现了两种粒度（粗粒度和细粒度）的混合并行，并基于国产in-house集群采用CRM标模（约4 000万网格单元）定常湍流算例对两种混合并行模式进行了测试和比较。结果表明，粗粒度在进程数和分区数较少的小规模并行时具有效率优势，16线程时效率较高；而细粒度混合并行在大规模并行计算时具有优势，8线程时效率较高。其次，验证了混合并行在非定常计算情况下的可扩展性，采用机翼外挂物投放标模算例，分别生成3.6亿和28.8亿非结构重叠网格，采用对等的（P2P）网格读入模式和优化的重叠网格隐式装配策略，网格读入和重叠网格装配耗时仅需数十秒；采用3.6亿网格，完成了非定常状态效率测试及非定常分离过程的湍流流场计算，在in-house集群上12 288核并行效率达到90%（以768核为基准），在天河2号上12 288核并行效率达到70%（以384核为基准），数值模拟结果与试验结果符合良好。最后，在in-house集群上采用28.8亿非结构重叠网格进行了4.9万核的并行效率测试，结果显示，4.9万核并行效率达到55.3%（以4 096核为基准）。     

主动光频标的研究进展

传统被动光钟激光锁频系统的频率稳定度最终受限于由布朗运动导致的腔长热噪声，是目前光频标发展面临的主要技术瓶颈之一。主动式新型光频标利用光学谐振腔的弱反馈在原子跃迁能级间形成多原子相干受激辐射，该受激辐射输出信号可直接作为钟跃迁信号使用，相比于被动光钟，主动光中原理上可实现更窄的量子极限线宽，并具有腔牵引抑制优势。分别介绍了国内外在主动光频标研究现状，主要介绍北京大学在该研究方向取得的理论及实验研究进展。为了实现窄线宽连续型主动光频标，提出双波长好坏腔方案减小剩余腔牵引效应对四能级主动光钟的影响，其中Nd:YAG 1 064nm好腔激光和铯原子1 470nm坏腔激光共腔输出，分别工作在好、坏腔区域，好腔激光通过Pound-Drever-Hall(PDH)稳频技术锁定主动光钟谐振腔的腔长，由于坏腔激光相比于好腔激光具有腔牵引抑制优势，因此坏腔激光线宽有望在腔长锁定之后的好腔激光的基础上进一步压窄。目前两套双波长系统的1 470nm钟激光拍频线宽受限于两个谐振腔腔长的相对抖动，为了抑制共模噪声对钟激光拍频信号的影响，并验证采用PDH稳频技术实现腔长锁定的可行性，本文利用相位锁定技术同步两套双波长好坏腔系统的谐振腔腔长，从而消除腔长共模噪声对钟激光拍频线宽的影响，进而分析除剩余腔牵引效应以外的其他因素对1 470nm主动光频标的影响。