跨声速条件下轴对称收扩喷管内外流场的数值研究

采用有限体积法,二阶精度插值的Roe格式对三维流场N-S方程进行空间离散;采用LUSGS隐式时间推进法和多重网格技术加快计算收敛速度;采用不同湍流模型、壁面函数对轴对称喷管内外流场进行三维数值模拟.计算结果与实验数据对比表明,采用涡粘修正的RNG模型,内流场使用增强壁面函数,外流场使用标准壁面函数可以得到准确的计算结果.研究表明:在亚声速条件下,喷管外壁面压力分布随落压比的变化不大;超声速条件下,喷管外壁面激波随着落压比的增大前移.      

轴对称推力矢量喷管的静态内部性能分析

采用Jamson提出的有限体积法对4种结构的轴对称推力矢量喷管的内流场进行了数值分析，预测了其内流场性能参数，计算结果与参考文献符合的很好，表明该方法对建立轴对称推力矢量喷管的数学模型具有很高的模拟精度和较快的收敛速度。     

出口宽高比对S形二元收敛喷管雷达散射截面的影响

在S形二元收敛喷管进出口面积、偏心距、面积变化规律及中心线变化规律不变的条件下,采用自适应超椭圆方法设计不同出口宽高比的S形二元收敛喷管.基于迭代物理光学法与等效边缘电磁流法自主开发计算腔体部件雷达散射截面（RCS）的程序;然后通过文献中的实验数据验证了计算程序的准确性和可靠性.通过计算程序分别分析了水平、垂直两种极化方式下不同S形二元收敛喷管的边缘绕射场与总散射场的电磁散射特性.结果表明:在水平、垂直两种极化方式下,喷管出口宽高比的变化对S形二元收敛喷管边缘绕射场的RCS影响较小,不同出口宽高比的S形二元收敛喷管边缘绕射场的RCS相差不超过4dB.喷管出口宽高比的变化对总散射场的RCS影响较大;正探测角时,宽高比为1.5时,在大部分探测角范围内总散射场有较低的RCS;负探测角时,宽高比为3.5时有较低的RCS.      

喷管型面对固体火箭发动机性能的影响

用时间相关法数值计算跨音速喷管流动时,其控制方程分别采用守恒型和非守恒型方程.通过对JPL喷管进行计算发现:在相同的精度下,非守恒型方程的收敛速度比守恒型方程的收敛速度快一倍多.采用修正阻尼法(DASHCO法)进行数值计算,其收敛速度均提高六倍以上.最后,采用直接优化法研究了喷管收敛段型面及喉部型面对发动机性能的影响.     

火箭发动机塞式喷管流场的数值研究

用时间相关法和张涵信的无振荡、无自由参数ＮＮＤ显式格式,对火箭发动机塞式喷管二维轴对称定常流场进行了数值模拟。计算得到了有粘和无粘中两种情况下,喷管流场的马赫数压强、温度、流线以及速度矢量分布图,计算网络采用代数方法生成,并通过求解Ｌａｐａｃｅ方程对网格进行修正。计算结果表明：ＮＮＤ格式用于计算塞式喷管流场具有较高的效率和良好的分辨率;无粘和有粘情况得到的流场结构有较大的差别;基础喷管的出口马赫数     

轴对称超临界粘性喷流的数值模拟

从雷诺平均的非定常Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程出发，采用隐式的Ｂｅａｍ－Ｗａｒｍｉｎｇ格式中的矢通量分裂技术，结合Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型，数值模拟了收敛喷管的内流场和喷射流场，通过对趋临界工况下喷管跨音速流场的数值计算，表明本文的结果与试验结果是吻合的。     

用隐式近似因子分解法计算喷管跨音速流场

本文用隐式近似因子分解法结合贴体曲线坐标计算了定常、无粘、跨音速喷管流场,选择了参数变化比较剧烈的小喉部曲率半径喷管,陡壁收敛段喷管和潜入喷管三种算例.计算结果表明,该方法收敛快,并且具有良好的精度.Courant数可以取至7.这种隐式法与显式法相比较,花费的机时少得多.本文的方法可望推广到求解N-S方程、PNS方程,和二相有粘喷管流动计算.     

轴对称喷管内流场的加权余数有限元计算

本文采用非定常技术,用加权余数有限元方法给出喷管内流场的数值计算方法,对收敛一扩散喷管的实例计算表明:本文计算的壁面和中心线M数分布曲线与实验结果是一致的。     

轴对称喷管内外流场与结构温度场耦合计算

利用有限体积的对流项二次迎风插值格式和重整化群(RNG)k-ε湍流模型,二层增强型壁面函数,同时利用球形谐波法考虑热辐射的影响,以灰气体加权模型(WSGGM)确定气体介质的辐射性质,求解N-S方程、热传导方程、考虑吸收-发射性气体介质的辐射传输方程。采用流固耦合的流动与换热模型,流场与结构温度场互为边界条件交换数据,实现了流场解算与温度场解算的耦合数值分析。采用此计算模型对静止的轴对称收-扩喷管进行了数值模拟,计算得到的喷管壁温与试验数据吻合良好。在此基础上,进一步开展了巡航状态下轴对称收-扩喷管内外流场与结构温度场的耦合数值研究。     

喷管加装小突片强化射流混合的PIV测量

在常温、超声速流动状态下,采用PIV技术定量测量了收敛喷管和收扩喷管加装小突片时射流纵截面的速度场和涡量场,并对其速度场和涡量场进行了比较。结果表明,对于收扩喷管,虽然其速度边界层厚度明显大于收敛喷管,但加装小突片结构也能够明显的强化射流的混合。并且发现,小突片结构不仅能产生流向涡,而且对周向涡有明显的强化作用。     

固体火箭喷管颗粒尺寸分级两相跨音速流场计算

本文中对固体火箭喷管颗粒尺寸分级的两相跨音速流场作了计算.气相控制方程采用隐式近似因子分解法求解,尺寸分级的颗粒控制方程采用特征线求解,然后,二者进行充分的耦合,可以获得固体火箭发动机含有任意颗粒质量分数和不同颗粒尺寸时轴对称喷管跨音速流场的参数分布.文中讨论了不同颗粒半径和质量分数对流场的影响,对单一颗粒尺寸和颗粒尺寸分级的参数进行了比较.两相耦合计算的迭代收敛速度取决于气相,本文中气相方程求解的格式除部分边界外是隐式的,CFL数可取至6左右,收敛速度快.特别是对颗粒尺寸分级的计算,得益更大,其得益的倍数为颗粒的分级数.     

涡轮导向叶片变截面针肋通道内空气流动和传热过程数值模拟

应用紊流ｋ－ε模型、三维适体座标网格结合ＳＩＭＰＬＥＳＴ算法,对涡轮导向叶片尾部变截面针肋通道内空气流动和传热过程进行了数值模拟。该模拟同时考虑了流场、压力场、温度场和物性变化之间的相互耦合作用,可压缩和变物性的影响以及叶片壁内三维导热效果,采取对流—导热耦合传热问题的整体求解法进行数值计算。计算结果与实验符合良好,证明本文建立的数值模型是正确的,可为叶片的优化设计提供一定的依据。     

航空发动机圆套状尾喷管流场温度场数值模拟

对某机型双收缩圆环套状尾喷管内, 冷热气流相互掺混的流动传热过程, 作了数值模拟。用SIMPLE思想在曲线系中进行程序设计。该程序可以对任意形状的两股流有(无)掺混时的流场、质量场、温度场、辐射通量场进行模拟, 对二维流动传热问题比SIMPLE程序更具有通用性。最后对某机型的实际情况(形状尺寸), 用区域拼装法生成曲线网格作了模拟。结果说明, 引射的冷空气与热主流在双收敛喷管内掺混, 确实能降低尾喷管内气流温度、降低红外辐射。      

激光燃速仪的燃烧室流场分析计算

 Numerical method of flow field of laser servo-controlled strand window bomb was studied and discussed. Stream function and vorticity method of axisymmetric, steady, incompressible flow field was also studied and discussed. Axisymmetric finite difference stream function-vorticity equations on interior points and the points near oblique boundary were derived. Boundary conditions were given and initial interative field was offered before calculation. Cross-bidirectional scanning method was used for increasing parameter transmission rate and relaxation method which could increase convergent stability was applied. At last, the flow field of laser servo-controlled strand window bomb was shown, including motion patterns and vorticity flow scatter diagrams. The solutions were in concordance to the solutions by Mills et al. and the theories by Pan and Acrivos.     

SIMPLE方法的敛散性

 通过对轴对称突扩直管和突扩后缩环管流场的计算考察了SMPL正方法的敛散性。在源项处理、方程的非线性问题和解法等方面进行了分析。在此基础上研究了迭代初场的构成和斜壁网格的处理。提出了一种可行的算法。     

涡扇发动机引射喷管的红外辐射特性数值研究

红外隐身技术对提高未来战机战场生存力具有重要意义,发动机排气系统是飞机后半球的主要红外辐射源。为了研究涡扇发动机引射喷管的红外辐射特性,结合引射喷管的CFD计算,采用离散传递法计算红外辐射强度,研究了普通引射喷管以及带5°下倾角的引射喷管后半球3～5μm以及8～14μm波段的红外辐射强度空间分布规律,并与相似尺寸的涡扇/涡喷发动机收缩喷管进行比较。结果表明：在探测角度较大时,涡扇发动机引射喷管的红外辐射强度较收缩喷管小20%左右;引射喷管结构上的非对称性导致红外辐射强度角向分布也呈现非对称性特征;8～14μm波段与3～5μm波段的红外辐射空间分布规律基本相同,但辐射能量小40%左右。     

圆管内不可压粘流的分布式并行算法

用并行算法求解了圆柱通道内不可压层流流动。采用有重叠的条形区域分裂法时计算域进行分区，分别求解动量方程和压力校正方程，对压力校正方程采用块修正法来加速线代数方程的收敛速度。计算结果表明有较好的并行加速比和并行效率。     

Pico-satellite Autonomous Navigation with Magnetometer and Sun Sensor Data

This article presents a near-Earth satellite orbit estimation method for pico-satellite applications with light-weight and low-power requirements. The method provides orbit information autonomously from magnetometer and sun sensor, with an extended Kalman filter (EKF). Real-time position/velocity parameters are estimated with attitude independently from two quantities: the measured magnitude of the Earth’s magnetic field, and the measured dot product of the magnetic field vector and the sun vector. To guarantee the filter’s effectiveness, it is recommended that the sun sensor should at least have the same level of accuracy as magnetometer. Furthermore, to reduce filter’s computation expense, simplification methods in EKF’s Jacobian calculations are introduced and testified, and a polynomial model for fast magnetic field calculation is developed. With these methods, 50% of the computation expense in dynamic model propagation and 80% of the computation burden in measurement model calculation can be reduced. Tested with simulation data and compared with original magnetometer-only methods, filter achieves faster convergence and higher accuracy by 75% and 30% respectively, and the suggested simplification methods are proved to be harmless to filter’s estimation performance.     

冲压发动机燃烧室三维流场数值模拟

采用有限体积法和S IM PLE算法对带V型火焰稳定器的冲压发动机燃烧室流场进行了数值模拟。目的是通过数值模拟得到燃烧室的流场,加以研究,指导实验。计算是在笛卡尔坐标系下进行的,计算中采用标准二方程湍流模型、EDM(Eddy-D issipation-M odal)燃烧模型和P-1辐射模型。得到了燃烧室流场的特征,并讨论了多种因素变化对流场的影响。结果表明,各种因素变化对燃烧室流场有较大的影响。