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

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

二维非结构网格的一个TVD型有限体积方法

考虑双曲型守恒律方程,对二维非结构三角形网格给出一种 TVD型有限体积方法,主要思想是在一阶单调格式的基础上,在每一个单元上对变量作单调线性重构函数,时间离散采用二阶 Runge-Kutta方法。通过计算分析了该方法的精度,对平面激波反射和空穴流动的计算结果表明该方法是成功的。     

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

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

非结构化网格生成中背景网格对源的响应

背景网格对源的响应方法对非结构网格的生成具有重要意义,如果背景网格参数能恰当地反映出源的变化,则所生成的非结构网格就能满足网格设计的要求。将背景网格的三角形单元分为7个子域,根据源在单元中的位置对网格进行局部加密或变形,使每一个源点都与一个背景网格结点相对应,再通过求解热传导方程,使背景网格参数分布光滑,基本反映网格设计的要求,从而提高了非结构化网格的生成质量。     

二维非结构网格的生成及其Euler方程解

 采用推进阵面法生成三角形非结构网格,用多重网格技术求解泊桑方程以取得结构背景网格.提出了一种高效的边界剖分法,并用它来生成初始阵面.另外,还提出了一种新颖的数据结构——数组链接表,用它来存贮网格数据.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼型和多段翼型的算例,并与实验结果进行了比较和分析.     

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

论述了以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主动光频标的影响。     

非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

GENERATION AND APPLICATION OF UNSTRUCTURED ADAPTIVE MESHES WITH MOVING BOUNDARIES

INTRODUCTIONMany engineering applications require thenumerical prediction of strongly unsteady flowsinteracting with moving bodies.A case in pointis that the detailed calculation of sphere motionsin an unsteady field induced by blast waves is akey issue to determine the drag force exerting up-on the sphere by the surrounding gas and,inturn,to reveal the mechanism of unsteady effectson the drag coefficient.To cope with such pro-blems,one should simulate unsteady flows withmoving boundary co…