高外压缩比的高超声速内乘波进气道设计

为了提高内乘波式进气道的起动性能,扩大其稳定工作的马赫数范围,以得到一种高外压缩比例的高超声速进气道,在ICFC二维基本流场的基础上进行改进,提高流场的外压缩比例,减小流场长度,并基于改进后的基本流场生成内乘波式进气道。数值模拟结果表明:在来流马赫数6的设计状态,该进气道保持了内乘波式进气道的优点:流量捕获能力强(流量系数99.94%),显著减弱外流阻力等;较基于原ICFC流场得到的进气道,该高外压缩比的进气道长度更短,且性能略有提升(出口马赫数为2.78时总压恢复系数为0.459),能在马赫数3.8~6范围内起动工作,较原进气道具有更佳的低马赫数起动能力。     

基于双自由曲面的LED均匀照明准直透镜设计

基于几何光学、能量守恒定律及菲涅耳定律等相关理论,提出了一种双自由曲面半导体发光二极管(LED)准直透镜的光学设计方法,并给出了构建准直透镜模型详细的算法设计。自由曲面是一种关于中心轴旋转对称的曲面,该曲面的二维轮廓在非均匀有理B样条曲线的方法理论基础上,采用ProE软件搭建而成。通过蒙特卡罗光线追迹模拟发现,相比传统的单自由曲面准直透镜,双自由曲面准直透镜不仅提高了照度均匀性,而且在能量利用率上也有显著的提高。研究结果表明,采用双自由曲面将大大提升准直透镜的设计空间,改善LED透镜的光学性能。     

高超声速进气道飞行器一体化设计技术的发展

对高超声速进气道与飞行器一体化设计技术和发展进行了研究,包含轴对称压缩、二维压缩、侧压以及内压缩进气道在高超声速飞行器上的典型布局设计方案.对高超声速可调进气道类型进行了概述,基于轴对称类型调节、二元平面类型调节以及三维内转调节进气道的典型案例给出了其各自的设计特点,并进一步对宽域飞行和组合动力飞行器采用的多通道可调节...     

旋转冲压压缩转子中压缩面的数值研究

采用二维雷诺平均N-S方程和Spalart-Allmaras湍流模型对三种压缩面形状的旋转冲压压缩转子二维进气流道流场进行了数值仿真.计算结果表明,压缩面形状对二维进气流道中的波系结构和形状具有决定性影响,凸曲线压缩面二维进气流道的总体性能最佳.     

同时实现多角度非垂轴面均匀照明的衍射光学器件

本文利用长焦深器件设计思路,提出了通过设计衍射光学器件相位分布使得在三个垂轴面上满足均匀照明要求来间接实现非垂轴面均匀照明的设计方法,并完成了同时实现多种角度非垂轴面均匀照明的衍射光学器件设计.基于爬山法和模拟退火法相结合的混合算法,设计了器件的相位分布,在30°、45°、55°非垂轴面上光斑不均匀性分别为5.91%、4.00%、3.38%.该方法可推广至二维情形,实现任意曲面均匀照明的要求.     

多电极半导体光放大器对增益特性的改善

提出了采用多电极实现电流非均匀注入来改变半导体光放大器(SOA)内部的载流子分布,从而改变其增益饱和特性的一种新的方案.用SOA分段模型的计算结果表明,与单电极均匀注入电流的SOA 相比,在同样的总注入电流下,对一个两电极SOA采用前面小后面大的非均匀电流注入,可以提高饱和输出功率和饱和输入功率.而采用前面大后面小的电流注入方式,则将使SOA更容易发生饱和,不仅饱和输出功率与饱和输入功率减小,而且增益谱在饱和后的压缩程度加大,增益谱压缩的对称性提高.对多电极SOA,可以方便灵活地调节各节之间的长度比和注入电流比来满足不同的应用要求.     

气动激光器喷管非平衡流2维数值仿真

将Anderson的两振型三温度弛豫模型和严海星整理的弛豫数据相结合,采用2维守恒型方程组对按照最小长度喷管型面设计方法设计的、面积比分别为50和20的气动激光器喷管非平衡流场进行了数值仿真。小信号增益计算结果在每个计算点都和J.S.Vamos等人针对这两种喷管的小信号增益测量试验结果符合很好.解决了传统的准1维非平衡流分析方法不能很好地和试验结果相符的问题.对气动激光器喷管性能设计提供了更精确的评估方法。     

三角B-B曲面在非结构网格生成中的应用

基于生成非结构化网格的Delaunay三角形化方法,应用三次三角Bernstein-Bézier曲面来控制网格点的分布情况。对于给定边界,利用三次非均匀B样条进行边界的拟合及重新离散。对初始化形成的Delaunay三角形,应用三次三角Bernstein-Bézier曲面来计算位于其中的点的长度标尺,通过对三角Bernstein-Bézier曲面边界点的定义来控制三角形内长度标尺的分布。对复杂通道网格剖分的实例表明,此方法可以很好地控制内部网格点的分布情况。     

无粘流中能捕捉较弱的接触间断和激波的反扩散差分格式

在研究弱入射激波遇到对称楔以后的马赫反射现象时,激波管实验不易测出很弱的接触间断,也不易捕捉到马赫反射与正规反射转换的条件。文章一方面研究了可压流体力学欧拉方程的数值方法,首先是用反扩散法改进接触间断的计算;另一方面根据格式粘性的特性和它引出的很微小的熵的变化规律来显示很弱的接触间断和反射激波。这样才易于将对三波点的分析推进一步。文[5,6]曾预言了一种反散波是连续的压缩波的新的激波反射类型。我们设想并根据计算初步确认这新类型反射实际应该是简单马赫反射,反射波虽弱仍是激波。     

气相爆轰波马赫反射影响因素的实验研究

在矩形截面的爆轰管道中,对C2H2+2.5O2+8.17Ar和C2H2+5N2O在CJ爆轰状态下经过不同楔面所发生的马赫反射的影响因素进行了实验研究。实验中,由烟膜记录爆轰波马赫反射的胞格结构转变过程;采用纹影技术捕捉爆轰波马赫反射波阵面的不稳定性及波后流场分布。实验结果表明:两种实验气体在爆轰波马赫反射过程中均存在由CJ区域向过驱区域转变的胞格结构;初始压力对楔面与马赫反射三波点轨迹线之间的夹角(χ)影响明显,楔角θw对χ的影响随θw的增大而增大;根据实验测得的θw+χ与θw之间的关系,可知爆轰波马赫反射三波点轨迹线的斜率随着θw的增大而增大,与CJ区域内胞格轨迹线的相交距离也更短,使马赫杆后的过驱度升高。另外,不稳定气体C2H2+5N2O的不稳定性高于稳定气体C2H2+2.5O2+8.17Ar,导致二者的爆轰波马赫反射行为存在较大的差异。     

气相爆轰波马赫反射影响因素的实验研究北大核心CSCD

在矩形截面的爆轰管道中,对C2H2+2.5O2+8.17Ar和C2H2+5N2O在CJ爆轰状态下经过不同楔面所发生的马赫反射的影响因素进行了实验研究。实验中,由烟膜记录爆轰波马赫反射的胞格结构转变过程;采用纹影技术捕捉爆轰波马赫反射波阵面的不稳定性及波后流场分布。实验结果表明:两种实验气体在爆轰波马赫反射过程中均存在由CJ区域向过驱区域转变的胞格结构;初始压力对楔面与马赫反射三波点轨迹线之间的夹角(χ)影响明显,楔角θw对χ的影响随θw的增大而增大;根据实验测得的θw+χ与θw之间的关系,可知爆轰波马赫反射三波点轨迹线的斜率随着θw的增大而增大,与CJ区域内胞格轨迹线的相交距离也更短,使马赫杆后的过驱度升高。另外,不稳定气体C2H2+5N2O的不稳定性高于稳定气体C2H2+2.5O2+8.17Ar,导致二者的爆轰波马赫反射行为存在较大的差异。     

吸气式高超声速飞行器机体/推进一体化内外流非定常耦合方法

基于吸气式高超声速飞行器机体/推进一体化的气动布局设计方式,文章提出了一种内外流一体化流场的耦合求解方法,其中燃烧室内流场采用考虑有限速率化学反应动力学模型的一维非稳态方法求解,进气道和尾喷管外流场采用二维CFD软件计算,进气道与燃烧室在耦合界面处通过一维平均方法实现静温、静压和Mach数等参数传递.并分别以日本国家航空与航天实验室（NAL）的氢燃料燃烧室模型作为内流场验证算例,以某典型高超声速飞行器一体化模型作为内外耦合流场验证算例.研究结果表明:有限速率化学反应准一维方法能较为准确地模拟燃烧室内燃烧流场,提出的内外流场耦合方法能够有效地计算出内外流耦合效应,计算后体压力分布与理论值较接近.该方法可为超燃冲压发动机的性能快速分析和吸气式高超声速飞行器机体/推进一体化的初步分析设计提供重要参考.      

激光光源汽车远光灯自由曲面透镜设计

根据激光光源的特点,设计了一种基于蓝色半导体激光器的汽车远光灯系统。采用以Monge-Ampere方程为基础的自由曲面透镜设计,解决了激光光源能量密度高、发光面积小的问题;依据非成像光学理论,建立了透镜表面与目标面的映射关系;通过求解Monge-Ampere方程得出自由曲面上的坐标点,拟合坐标点得到自由曲面透镜模型;采用光学软件Tracepro进行了仿真分析。结果表明:照度最大值值为119.3lx,中心点照度为117.6lx,完全符合GB25991-2010标准,且光型分布均匀。     

组件计算程序PANDA研发及初步验证

PANDA是上海核工程研究设计院研发的压水堆组件计算程序。该程序采用基于特征线方法（MOC）的一步化计算流程,即在不引入能群压缩和栅元均匀化的情况下直接进行组件层面的两维非均匀输运计算。多群数据库采用基于ENDF/B-VI制作的70群结构中子数据库,基于ENDF/B-VII的新版数据库也正在开发中。共振自屏计算采用了空间相关丹可夫方法（SDDM）,既具备燃料芯块分区计算的能力,又保留了传统Stammler方法的计算效率。多群非均匀输运计算采用二维模块化MOC方法,并辅以双重粗网有限差分(CMFD)加速技术,具有良好的计算精度和效率。对传统线性子链解析（TTA）方法以及多种矩阵指数方法进行了研究,选取了适合PANDA程序燃耗链的燃耗方程求解技术。基于以上基本模型开发了PANDA程序,并从程序模块、总体集成和核设计程序系统确认等三个层面,初步验证了PANDA程序的计算性能,表明了PANDA程序的工程设计计算能力。     

基于爬坡算法的片上低栅瓣二维光学相控阵

通过两级定向耦合器结构实现了光学相控阵天线阵列的二维排布,并设计了特殊C形弯曲波导进行热调作为阵元的相位控制器.提出一种针对阵元数量不高的稀布片上光学相控阵远场高阶干涉栅瓣的压缩方法.将均匀的阵元间距通过爬坡算法优化成非均匀间距,以破坏栅瓣产生所需的干涉相长条件,实现对栅瓣的压缩作用.在1 310nm波长,通过时域有限差分法对基于微型光栅耦合器天线的稀布阵列进行计算分析,结果表明优化的阵元间距能实现-6~-7dB的栅瓣抑制比.     

基于微分方程数值解的自由曲面透镜的设计

为满足目标面上均匀照明的需求,设计了一款自由曲面透镜。根据LED光源的发光特点,结合光学成像的特性,非成像光学原理和能量守恒定理,推导出实现光能量在目标面上均匀分布的自由曲面面形的微分方程,采用matlab的ode算法求出面形上的离散坐标点,对离散坐标点拟合后得到透镜模型,通过光学模拟仿真软件对透镜模型进行光线追迹。结果表明,配光角度为80°的透镜,透镜的口径与光源发光面宽度之比大于等于9时,目标面上的照度均匀性大于0.9,光能利用率约为85%。该设计方法可以快速精确地设计出所需的透镜,而且透镜结构紧凑,单颗就能实现均匀照明,有利于LED光源照明系统的小型化。     

基于弯曲激波特征线法的两级压缩内收缩流场设计与分析

进气道是超声速飞行器的关键部件,设计核心是无黏超声速流场的确定,而波系结构又直接决定了流场的性能。考虑到实际应用价值,基于给定激波的超声速流场反设计至关重要。传统单道入射波的流场结构简单,压缩效率不高,且传统特征线法（method of characteristics,MOC）无法获得高阶的气动参数。为了拓展设计思路,首先利用弯曲激波特征线法展示了内收缩超声速流场不同的单元过程,随后基于此法提出了一种已知入射和反射激波的两级压缩内收缩流场的概念及反设计方法。基于此法设计的基准流场与数值模拟结果吻合良好,双入射激波可以提高压缩效率并缩短长度,能够实现波系结构和出口参数的分布可控,在非均匀来流条件下也能得到相应型面。在给定激波的条件下,求解分析了一系列具有中心体的轴对称流场,讨论入射和反射激波激波角分布对流场气动参数和几何参数的影响。弯曲激波特征线法的精确性和有效性使其成为平面/轴对称超声速流场反设计的良好候选。      

JF12激波风洞高Mach数超燃冲压发动机实验研究

针对高Mach数（Ma ≥ 7）超燃冲压发动机高气动阻力下的燃烧组织问题,提出一种双突扩燃烧室结构方案.使用数值模拟方法考察了射流与双突扩燃烧室组合方式的混合燃烧特性.设计了双突扩超燃冲压发动机模型,在力学研究所JF12长试验时间激波风洞内,开展了Ma=7.0和Ma=9.5的氢燃料点火和燃烧试验对比.在风洞有效试验时间100 ms内,实现了Ma=7.0和Ma=9.5超燃冲压发动机的成功点火与稳定燃烧.在Ma=7.0情况下,进气道采用三维压缩,燃烧室入口设计Mach数Ma_c=2.5,壁面压力分布实验结果显示燃烧放热靠近燃烧室扩张段上游;在Ma=9.5情况下,进气道采用二维压缩,燃烧室入口设计Mach数Ma_c=3.5,由于燃烧室流动速度特别高,燃烧放热靠近燃烧室扩张段下游.      

基于曲面分割的三维点云物体识别

提出了一种基于局部描述符的三维点云物体识别算法.算法首先得到点云的邻域、法线矢量等相关信息,通过邻域进一步得到形状索引值.特征点的提取以形状索引值为依据,以每个特征点为基点对曲面根据欧式距离和矢量夹角分割.分割的曲面片进行等间距划分为多个欧氏距离同心圆,以特征点切平面为基平面投影,并进行等角度抽样,通过抽样点相对特征点的法线矢量及测地距离变化曲线,建立曲面片的二维描述,从而把三维识别转化为二维.根据算法建立模型数据库,给定一个物体,通过和模型数据库中的曲面描述进行比对,得到潜在的识别结果,再通过迭代最近点算法,得到最终的识别结果.最后,通过大量具体实验验证了算法的有效性,并给出了算法的计算复杂度及耗时对比分析,说明了算法的高效性.     

基于曲面分割的三维点云物体识别

提出了一种基于局部描述符的三维点云物体识别算法.算法首先得到点云的邻域、法线矢量等相关信息,通过邻域进一步得到形状索引值.特征点的提取以形状索引值为依据,以每个特征点为基点对曲面根据欧式距离和矢量夹角分割.分割的曲面片进行等间距划分为多个欧氏距离同心圆,以特征点切平面为基平面投影,并进行等角度抽样,通过抽样点相对特征点的法线矢量及测地距离变化曲线,建立曲面片的二维描述,从而把三维识别转化为二维.根据算法建立模型数据库,给定一个物体,通过和模型数据库中的曲面描述进行比对,得到潜在的识别结果,再通过迭代最近点算法,得到最终的识别结果.最后,通过大量具体实验验证了算法的有效性,并给出了算法的计算复杂度及耗时对比分析,说明了算法的高效性.