迎风格式在二维非结构网格中的应用

将原来在结构网格中提出的AUSM⁺(Advection Upstream Splitting Method)格式,推广应用于二维非结构网格中.利用格林-高斯公式对控制体内的变量进行线性重构,获得空间二阶精度,并采用Barth型限制器以抑制数值解的振荡.为提高计算效率,采用了自适应多重网格法.最后给出了NACA0012翼型和平板激波反射的几个Euler方程的算例,并进行了分析比较.     

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

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

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

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

战术导弹横向喷流数值模拟

为了研究喷流对导弹不同升力面的干扰影响,通过数值求解N-S方程来模拟超声速外流场中横向喷流的干扰流场,采用分块对接网格和"O"型网格技术,精确模拟喷口截面及弹翼形状,生成高质量的贴体计算网格.通过几种升力面和弹身组合的超声速横向喷流干扰流场的数值模拟,计算的马赫数为2和4.5,计算结果显示具有升力面的喷流干扰比光弹身的喷流干扰大.研究和分析了喷口附近流场的涡系结构和波系结构,并将导弹喷流干扰力放大因子及喷流引起的压心变化的计算结果与目前国外公布的试验值进行比较分析,二者显示了良好的一致性.     

非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ENO(Essentially Nonoscillatory Scheme)型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法.最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算.      

基于SLA的网格监控体系GSLMA及其实验床

从体系结构的角度出发,将现有网格监控工具分为3类,对每一类网格监控结构进行了分析 和对比,并提出了一种基于服务级别协定SLA(Service Level Agreement)的网格服务监控体系结构GSLMA(Grid Service Level Monitoring Architecture).GSLMA采用基 于面向服务的体系结构SOA(Service Oriented Architecture)的松散耦合结构,通过动态部 署解决了对新增服务的监控扩展,并提供第三方的监控功能,是一个具有高灵活性和高可伸 缩性的网格监控体系结构.最后,介绍了实现GSLMA实验床中的关键技术,包括:对开放网格 服务体系结构OGSA(Open Grid Service Architecture)容器进行扩展,采用Web服务级别协定WSLA(Web Service Level Agreement)解决被监控服务信息的统一描述等等.      

同轴微阴极电弧推力器的粒子网格法数值模拟

采用粒子网格方法对同轴微阴极电弧推力器μCAT工作过程进行了模拟研究, 并应用自相似方法对模型进行简化,获得了推力器羽流区的电子数密度分布、离子数密度分布、电势分布及离子轴向平均速度,通过改变磁感应强度和位形分析磁场对推力器内等离子体运动特性及推力器性能的影响。计算结果表明,电子被外加磁场捕获约束在磁力线附近,低速离子与高速电子形成的双极扩散电场加速离子喷出;在相同流量情况下,磁感应强度002T时,离子返流严重,磁感应强度005~030T时,磁感应强度变化对速度影响较小;磁场位形对离子运动和推力器性能有较大影响,磁力线与轴线夹角较小时离子速度下降明显,夹角较大时离子返流严重。     

三角网格曲面模型快速分层算法

提出一种三角网格曲面模型快速分层算法,该算法基于R*-tree建立三角网格动态空间索引结构,依据索引结构数据结点的分布状况计算各层截平面的位置;采用深度优先遍历方法获取与截平面相交的三角面片集合,并计算该集合中各面片与截平面的交线,将交线首尾相连,生成截面轮廓线,实现三角网格曲面模型的快速分层;实例证明该算法可对各种复杂三角网格曲面模型进行分层,算法准确、稳定,运行效率高.     

复杂外形网格生成技术及应用

以某型双立尾战斗机为例,给出了实际复杂型号飞行器结构网格生成的一种分块方案及实现过程,用一种矢性三次多项式插值的方法生成相邻块公共交界面的网格,用Hilgenstock方法修正源项,求解椭圆型方程生成块内空间网格,使网格线与边界面夹角、与边界面间距控制为期望的值.在所生成的双立尾战斗机分块结构网格中求解Euler方程,获得了合理的全机气动力系数、翼面压力分布和绕流图谱结果.     

r型网格自适应在间断Galerkin有限元激波捕捉中的应用

间断Galerkin(DG)有限元方法因计算精度高、适用于非结构网格等特点得到广泛研究和应用,其在数值模拟包含强间断流场时存在残差收敛性和计算鲁棒性差问题,均匀分布的网格加剧这一问题并影响激波分辨率。针对该问题,发展了r型网格自适应方法,实现间断Galerkin有限元数值模拟过程中网格自适应加密。基于网格点归一化的压力值作为r型网格自适应中网格点移动驱动力的重要权值,并将网格自适应后的网格点位移变化量与网格点之间的初始位移之比作为驱动力的另一重要权值,实现网格沿激波方向各向异性自适应加密,并且激波附近网格点的相邻网格点同步向激波方向移动。发展了适合间断Galerkin有限元方法的Venkatakrishnan限制器。并列NACA0012翼型超声速算例及三维并列圆柱相互干扰算例结果表明：基于r型网格自适应的间断Galerkin有限元方法能够清晰锐利捕捉激波,提高模拟精度,具有良好的收敛性和鲁棒性。     

反演航天器在轨瞬态外热流的导热反问题

获得航天器在轨飞行过程中的外热流数据对于研究热控涂层在轨退化规律、各种空间因素对热控产品的影响以及航天器姿轨控发动机羽流热效应都有非常重要的意义,然而直接测量热流存在很多困难,因此可以通过求解导热反问题得到满足一定精度的结果.首先,通过研究利用航天器设备在轨遥测温度值反演出航天器在轨瞬态外热流的导热反问题方法,建立了反演航天器在轨瞬态外热流的数学模型,采用共轭梯度法求解导热反问题并从物理概念角度改进了共轭梯度法的迭代过程以增加其抗不适定性;然后构造了两组能够代表目前大多数地球轨道航天器以及深空探测航天器在轨吸收外热流变化的数值试验对共轭梯度法的反演效果进行了检验.除阶跃变化位置以外反演值与真实值的最大相对偏差为2.9%,反演效果非常好;对于阶跃变化位置的吸收外热流在对反演结果进行分析处理后也能够得到较好的反演结果.      

低雷诺数下翼型前缘流动分离机制的研究

采用高精度有限差分格式,对来流雷诺数为1.0×104,攻角为3°的二维翼型流动进行了直接数值模拟,研究了低雷诺数下翼型前缘流动的分离机制,描述了分离涡系的相互作用规律.计算结果表明:前缘椭圆弧靠近叶身位置存在吸力峰,流动在吸力峰内强逆压梯度的作用下发生分离;翼型上表面形成了包含驻留涡、脱落涡和二次涡的涡系结构,其尺度随时间不断变化,具有强烈的非定常性;表面压力分布曲线可以较好的描述翼型边界层流动.      

非结构三角形网格的一个ENO型有限体积方法

基于双曲型守恒律方程,对非结构三角形网格给出了一种ＥＮＯ（ＥｓｅｎｔｉａｌｙＮｏｎｏｓｃｉｌａｔｏｒｙＳｃｈｅｍｅ）型有限体积格式,方法的主要思想是先对每一个三角形单元构造一个加权的二次插值多项式,而在计算交界面的流通量时采用了两点高斯积分公式以保证格式的整体精度,时间离散采用三阶ＴＶＤＲｕｎｇｅＫｕｔａ方法．最后给出了该格式收敛的数值阶,并对前台阶问题进行了计算     

一种战斗机外形涡流场Euler方程数值模拟

 用三步显式格式时间推进求解物理空间曲线网格上有限体积离散的Euler方程,数值模拟一种战斗机外形的涡流场.计算的机翼表面压力分布与实验符合较好,用Euler方程捕获到机翼前缘分离涡(主涡)及翼面上的二次分离涡.Euler方程解中出现的二次涡可能是逆压梯度、尖锐边缘和人工粘性共同作用的结果,它的流谱与实验定性符合.     

腔体闭式流动控制的实验和数值模拟

利用数值模拟和风洞实验相结合的方法,研究了闭式流动腔体的流动特征及其设置圆柱控制杆后腔体内声压级(SPL, Sound Pressure Level)和压力分布的变化.数值模拟求解三维N-S方程,采用AUSM+计算格式,湍流模型采用Wilcox k-ω模型.实验在0.6 m×0.6 m超音速风洞中进行,在腔体底部布置了40个常规静压测量点和15个动态测压点.研究表明,在外流为超音速流时,闭式流动的腔体底部压力变化梯度较大,腔体底部和后缘的测压点的SPL值和频率关系曲线中没有明显的SPL峰值.实施控制后,腔体底部的压力变化梯度减缓,在腔体后缘分离区内的测压点SPL值降低,而前缘分离区内的测压点SPL值增加.      

基于SA湍流模型的DES方法应用

使用基于SA(Spalart-Allmaras)湍流模型的分离涡模拟(DES,Detached-Eddy Simulation)方法,对圆柱绕流和超声速圆柱底部流动等高雷诺数大分离流动进行了数值模拟.无粘通量的计算采用二阶迎风格式,时间推进使用二阶隐式双时间步法.计算表明相对于雷诺平均N-S方程(RANS,Reynolds-Averaged Navier-Stokes)模拟,时空均为二阶精度的DES模拟能够分辨更为精细的分离区三维旋涡结构,对分离区主要湍流结构的捕捉和平均压力分布的预测取得了满意的结果.     

热毛细对流模式转换及温度振荡特征

研究了矩形液池中由于两端温差引起的热毛细对流的温度振荡临界条件.在实验室中,设计了一个高分辨率的温度测量系统,用于实时观测并记录流体的温度.该系统主要由热电偶温度传感器、纳伏表和数据采集电脑3部分组成.得到了各种实验条件下温度振荡的临界条件,并且讨论了它与Prandtl数和Bond数之间的关系.利用flow3d软件数值模拟了微重力条件下的热毛细对流,发现了一种由于自由面变形和液层流场相互作用导致的晃荡的现象.     

高阶熵条件格式下Euler方程与 N-S方程的混合算法

针对在工程问题中使用N-S方程计算量大的问题,提出只在粘性作用显著的局部区域求解N-S方程,在流场的其余部分采用Euler方程计算的办法来模拟具有复杂外形的流动问题,并提出了一种基于高阶熵条件格式的算子分裂算法.应用所构造的算法对绕迫击炮弹的亚、跨超声速流动进行了模拟计算.计算结果与试验数据的对比表明,所提出的方法是切实可行的.