二维非结构网格的可压缩Euler方程流场解

讨论了非结构网格的生成方法,用Delaunay三角化方法和推进阵面法相结合的方法生成了非结构网格,提高了生成非结构网格的自动化程度和网格的质量.用有限体积法在非结构网格中求解了二维可压缩Euler方程,给出了单段翼形和两段翼型的算例,并与实验结果进行了比较和分析.     

关于非结构网格生成中的背景函数

背景信息对于用阵面推进法生成非结构网格是一个关键的环节.文章基于众多文献中提到的背景信息提供方法,提出一种更为实用的背景信息提供技术.此方法不仅避免了不必要的机时,内存的消耗,又使网格具有光滑性,网格的生成具有灵活性,方便性的特点.基于生成的网格采用有限体积法研究无粘高超音速化学非平衡流的流场分布,生成的网格能够成功的运用于高超音速化学非平衡流的数值计算。     

乘波机外形设计

计算了由锥型流理论、吻切锥理论及吻切轴对称理论生成的不同乘波机外形.结果表明,圆锥绕流流场生成乘波机外形的方法虽然比较简便,且生成的外形具有较高的升阻比,但进气道进口流场具有锥型流的性质而不够均匀,给发动机的工作带来不利影响.应用吻切锥、吻切轴对称理论与特征线方法结合起来进行反设计,可以得到更一般的乘波机外形,并且能有效地改善进气道进口流场,提高容积效率.此外,这几种方法的计算速度都很快,对进一步的优化处理及机身-发动机一体化十分有利.      

重叠网格技术中割补法的研究与改进

针对计算流体力学中重叠网格生成困难的问题,分析了重叠网格自动生成技术中割补法的主要特点和基本流程,指出其在洞映射单元属性的判断方法和洞边界外推过程中出现的"孤点"现象等方面的问题,提出了合理的改进方案,并将改进后的方法应用到自适应嵌入网格和DRAGON网格的生成中,最后编制出通用程序并给出了三维算例.实际应用表明,改进后的方法具有更稳定、自动化程度更高、适用范围更广的特点.     

燃烧室三维贴体网格的生成方法研究

研究了燃烧室三维贴体网格人机交互生成的方法.使用该方法可以生成包括扩压器、内外环通道和火焰筒在内的燃烧室流场计算所需的整个网格.可以对网格的疏密和正交性进行调整.可以处理火焰筒上的主燃孔,掺混孔和气膜孔.并能够输出边界曲面上网格点的类型,用于流场计算.网格边界线的输入和边界线上的网格点分布采用人机交互方式生成.使用VISUAL BASIC语言编写程序界面.把整个网格生成及图形显示功能集成在一个软件中.     

结构/非结构网格中的数值优化计算

非结构网格因其在处理外形方面所具有的特殊优越性而倍受关注.在结构/非结构网格中做了粘性流动的数值优化计算,优化计算包括单目标和双目标气动优化,以及反设计计算,即从初始翼型不断修型而实现给定翼型的外形及压力分布的目标.算例及其结果分析证明了该方法用于翼型的修型和设计计算的有效性.      

网格结构对Euler方程解的影响

 采用2种不同对接分区网格技术讨论了不同的网格结构对绕偏副翼的三维机翼流动的Euler方程解的影响.一种是部分不连续的网格(PDG);另一种是完全连续的网格(CCG),分别在PDG和CCG中采用Van Leer的矢通量分裂格式和Jameson的中心差分格式求解了Euler方程.算例表明,2种网格各有其优缺点,且都能得到有用的结果,但为了高效地取得更好的计算结果,必须根据流动和几何外形选取恰当的计算网格.     

三角形机翼参数化有限元网格划分与调整方法

为缩短有限元建模周期,提高三角形机翼结构分析、设计与优化的效率疲,首先,以其满足预定气动性能的几何外形为输入,定义了易用的有限元节点与单元的编号规则;基于可设计的机翼内部构型参数及可变的有限元尺寸参数的设置,引入了翼肋贯穿截止准则以满足任意输入的内部构型的初步判断;借助自定义形状矩阵完成了节点布置与单元生成,进一步通过有限元网格细化完成了开口设置、翼肋贯穿位置修正及桁条建模.然后,基于已建立的有限元网格实现了内外侧副翼翼肋位置的小幅调整及旋转舵面的角度调整,以满足不同飞行状态下结构分析需求.最后,应用PCL语言开发了参数化建模模块,实例表明了方法的有效性和模块可靠性.      

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

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

飞行器RCS计算前置处理中裁剪曲面剖分算法

在雷达散射截面(RCS)计算中,首先需要将物体的表面进行三角剖分,称为前置处理,三角片单元的边长通常为0.1个波长(一般1个波长取3 cm).飞行器外表形状复杂,在建立数学化模型中,使用了大量的裁剪曲面,曲面之间的关系非常复杂,而3 mm边长的三角片剖分,使得飞行器的三角片单元的数目巨大,给前置处理带来难度.对此,提出了一种裁剪曲面三角网格剖分的方法.利用飞行器外形设计的特点,在以拱高为逼近误差的前提下,把曲面离散一系列等参数线,等参数线上等弦长布点,两两参数值相邻的参数线生成网格单元.根据网格单元与裁剪区域的位置关系,将落在裁剪区域外的网格单元进行裁剪处理.然后对曲面间边界处的网格进行相容性处理,最终实现多张裁剪曲面的三角网格剖分.      

飞机表面绕射射线的寻迹方法

研究一种电子对抗特种飞机表面绕射射线的寻迹新方法.对飞机进行快速三维几何建模和表面网格剖分,并根据需要对几何模型的剖面和结点加密,生成飞机表面网格数据;选定寻迹区域,对飞机表面网格结点进行管理,根据费马原理,采用改进的动态规划方法,在飞机表面网格上快速、准确地实现任意源点和绕射接收点间的短程线计算;在此基础上,采用递归方法实现可视化的几何绕射理论(GTD)、一致性几何绕射理论(UTD)绕射射线离散寻迹.绕射射线离散寻迹为数值方法求解机身表面任意两点间的天线隔离度提供了条件.      

混合重叠网格插值方法的改进及应用

混合重叠网格间合理插值是保证流场计算正确的基础之一。本文针对重叠区网格尺寸匹配度较差时插值误差较大问题,发展了一种新型混合重叠网格插值方法,通过使用二阶精度插值、按单元尺度区分和扩充模板,改善了插值区网格尺寸匹配度较差时的插值精度。本文方法适用于任意单元类型的混合网格重叠,对各类单元处理透明,实现简单。计算结果表明,采用本文方法在网格重叠区流场变量传递正确有效,插值区网格尺寸匹配度较差时,相比原始方法,等值线过渡更为光滑,变量经过插值区耗散更小,计算与试验值符合更好。     

3点逐步r型自适应网格算法

为了根据某个流场参数梯度较合理地分布网格同时又不增加网格点的数量,提出一种基以总变差为基础的r型3点逐步自适应网格算法,通过3点逐步调整网格点的位置以减小网格线的扭曲,生成质量较好的计算网格.以绕二维圆柱和三维双椭球的超声速流动为算例,用Euler方程模拟,以流场的压力总变差来进行自适应,通过分别用所提出的自适应方法生成的自适应网格与用原始网格计算所得的流场结果对比,计算结果表明,用自适应网格计算所得的激波比用初始网格计算所得的激波薄,采用自适应方法所计算出的流场数值解具有更高的分辨率.      

A simulation analysis of systematic errors in ion arrival angle and drift velocity from the satellite borne ion drift meter

As one payload of a Chinese seismic satellite program, an ion drift meter (IDM) will measure drift velocity of thermal ions at an altitude of 500 km. Previous works have shown that such instruments use biased grids to create nonuniform potential in the grid planes, which brings systematic errors to the inferred parameters. A commercial finite element analysis software is used to simulate this instrument in the exact size. The error sources from thermal velocity, nonuniform transparency of real grids and potential depression in the grid planes are explained. The simulation results show that the arrival angle and drift velocity will be underestimated in all the conditions and the maximal error will be about −0.87° and −121 m/s, respectively. Furthermore, the relative error of the inferred arrival angle and the drift velocity will be inversely correlated with the arrival angle because of the lensing effect of the potential depression. This simulation provides a quantificational method of evaluating and correcting the data during in situ operation.     

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

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

主动光频标的研究进展

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

飞翼布局气动方案优选和试验验证

提出综合考虑飞翼布局隐身性能和气动性能的平面形状特征参数,分析了方案优选中的约束关系.采用参数化方法构建了三维曲面模型,并将物面网格划分流程进行自动化封装,通过更改设计参数准确快速地得到新方案的物面网格.应用基于Euler方程的数值方法进行布局方案的气动性能计算分析,在重点方案的分析中加入黏性修正;应用高频近似方法估算方案的隐身性能.以巡航状态作为设计点,在隐身性能的约束下,应用分析-修正的方法完成了气动布局方案优选,并对最终选定的方案进行风洞试验验证,证明该方案有进一步研究的价值.     

高超音速热流数值计算中的误差匹配原则

通过对粘性边界层的数值模拟和误差分析,提出了差分格式的误差匹配原则.在边界层内的计算网格一般都具有很大的长宽比,而计算误差与网格长宽比密切相关.为了使不同坐标方向上的误差保持在同一量级上,需根据误差匹配原则确定计算网格的尺寸.对绕钝锥的高超音速流动进行了数值模拟,进而计算了钝锥表面的热流分布.数值计算表明,参照误差匹配原则进行网格划分是有益的,可在保持流场数值计算精度和提高热流计算准确性的同时,使网格的数量和分布更为合理.     

大长宽比单元有限元误差分析及高精度计算

针对大长宽比单元将传统的基于单元整体尺度的有限元误差估计扩展为基于单元不同方向尺度的误差估计;根据不同方向误差应该同量级的思想,得到了误差匹配准则.该准则可以作为网格划分的判据.根据这一匹配准则,提出了一种提高计算精度的方法——单向高次插值法.数值实验表明该误差估计是正确的,误差匹配准则作为网格划分的判据是有效的,采用合理的有限元插值逼近能够有效地减小计算误差,提高计算的精度和效率.