三维复杂外形非结构网格的生成及应用

以 Delaunay方法为基础生成了三维复杂外形的四面体非结构网格。这种方法的特点是先布结点 ,然后按 Delaunay准则将其连接起来。与其他方法相比 ,其好处是能较好地根据计算需要控制网格点的疏密 ,提高计算效率 ,其运行过程也较稳定。在布点过程中采用了背景网格 ,使点的生成自动而合理。文末给出了 M6机翼和简化航天飞机的数值算例。     

非结构网格生成技术

描述了非结构网格生成的过程,网格的生成方法基于Delaunay 三角形方法,并对物面附近网格进行拉伸以便获得在边界层内高拉伸比的三角形单元。网格的生成过程分两部分:(1) 按Delaunay 三角形方法生成适于无粘情况下的非结构网格。(2) 对近物面的部分三角形单元进行拉伸,使生成的网格能适用于粘流。使用上述方法本文构造了二维弯管内绕叶栅区域的非结构网格     

计入桨叶运动的旋翼CFD网格设计技术

提出了一个考虑桨叶实际的挥舞和变距运动的旋翼CFD网格生成方法。该方法首先采用改进的几何法在桨叶表面生成初始网格，通过求解Poisson方程得出贴体正交网格，然后结合嵌套网格和动态网格的思想，给出了一种运动嵌套网格来计入桨叶的挥舞和变距运动．并针对嵌套网格中关键的贡献单元搜寻问题，提出了一种新的“伪贡献单元搜寻法(PSDE)”。该搜寻方法不仅可以适用于在结构网格中搜索贡献单元，而且不作修改即可应用于非结构网格，还可克服“InverseMap”法在结构网格中需要单独处理的奇异线(面)问题。     

基于扫掠法的六面体网格生成算法及实现

提出了一种基于扫掠法的六面体网格生成算法,该算法通过源和目标曲面参数域间的近似最小二乘仿射映射将源曲面的网格投影到目标曲面,并利用背景网格插值的方法计算内结点,从而将2.5维实体生成六面体网格。该算法避免了空间网格投影过程中耗时的“寻根”过程,并且不需要导致效率降低的网格优化,能够处理各种2.5维实体,生成的网格质量较好。     

一种新型的直角网格生成方法及应用研究

提出了一种基于背景网格的直角网格新型生成方法。背景网格是用来提供网格加密的信息．通过网格尺度与所得计算参数的比较，定量判断确定网格是否需要加密；而在求解Euler方程计算中，采用的是Jameson的有限体积法，它可以适合任意形状的网格单元，以及四步Runge—Kutta时间推进。本文对NACA0012翼型和复杂多段翼型等问题进行了数值实验。结果表明，这种网格生成方法易于推广到三维情况中去，具有网格生成时间短．收敛速度快，易于处理复杂边界的优点。     

自适应计算网格的生成和优化方法

本文提出一种简单通用的自适应计算网格生成方法。配合文[1]的网格优化方法,能在生成自适应网格的同时,保证网格具有较好的光滑性和正交性。本文还阐述了单元加权函数W的取值和边界点的确定对自适应网格性质的影响。作为应用实例,本文将自适应网格用于二维无粘跨音速叶栅流场计算。结果表明,采用自适应网格能明显改善流场的计算精度。     

基于切割单元直角坐标网格的DSMC方法及优化技术

对传统直角坐标网格进行改进,根据计算机图形学对与物面相交的正方形网格进行切割,使单元与物面贴体,并能根据流场的变化对网格区域进行局部加密。结合直角坐标网格与非结构网格的特点,编写了自适应直角坐标网格的生成程序,并发展了基于切割单元直角坐标网格的直接模拟蒙特卡洛(Direct simulation Monte Carlo,DSMC)方法。同时通过当地模拟分子数、动态时间步长技术优化DSMC方法,提升计算效率。数值结果表明,本方法在计算过程中使搜索效率得到大大提升,同时也保证了物面网格单元的贴体性,提高了物面附近流场的计算精度,还易于实现网格自适应。     

一种基于块结构的非结构网格快速生成方法

以阵面推进法为基础,提出了一种用预制模块的方法实现在平面域内生成非结构网格的新方法,即模块方法。模块是由一组三角形构成,文中以三角形模块为例,介绍了网格模块的预制方法。网格的生成是根据空间步长的背景信息,从预存的网格模块中选取当地合适的模块,直接复制或变形到计算区域实现的。由于生成过程中每一步引入的是预制的网格模块而不是传统的单个三角形单元,因此可期望新方法要比传统方法快得多。算例显示了新方法的快速网格生成过程,该方法具有推广应用前景。     

围绕旋翼飞行器的三维结构化运动嵌套网格生成方法

考虑到旋翼飞行器的旋翼各种运动和机身复杂外形特点,建立了一种鲁棒性高的三维结构化嵌套网格生成方法。采用三维Poisson方程对初始网格进行迭代,获得了高质量的旋翼桨叶三维结构化网格。以此方法为基础,以复杂外形几何体表面网格为初始网格,建立了一种抛物型法向外推方法。推进过程中,使用经抛物化处理的新Poisson方程进行光顺迭代,源项由Hilgenstock法确定,有效克服了椭圆型方程生成网格时边界点较难调整的困难。通过网格最小正交度的检测,表明了该抛物型方法的先进性。进一步将高效的挖洞法("Top map"法)和贡献单元搜寻法("Inverse map"法)引入上述三维结构化网格生成中,建立了围绕旋翼飞行器的运动嵌套网格生成方法,采用RANS方程分别对Caradonna-Tung旋翼和V-22倾转旋翼机流场进行计算流体力学(Computational fluid dynamics,CFD)计算,表明该网格生成方法能够满足CFD计算的需求。     

用于弹性旋翼流场模拟的网格变形方法研究

建立了一个基于改进"Ball-Vertex"的具有拓扑无关性的网格变形方法,用来精确模拟桨叶发生弹性变形后贴体网格的运动以及旋翼近场气动力环境,为更好地分析弹性变形对桨叶贴体网格的影响提供了一套解决方案.该网格变形方法从基本的弹簧模拟思路出发,用一网状弹簧系统来模拟桨叶贴体网格的变形和运动.在网格结点周围引入"Ball-Vertex"概念,通过添加相关的冗余约束,该方法保证网格结点在变形过程中不会越过与该点所对应的对角面,从而避免了畸形网格单元的产生,可有效地用于较大变形情况下的流场数值计算.使用所建立的方法,首先对绕二维翼型的结构网格和非结构网格进行了网格变形模拟,以表明方法的有效性;然后,将该方法应用于振荡翼型的流场数值求解,进一步验证了所建立的方法;最后,结合该网格变形方法和旋翼流场求解方法对弹性旋翼流场进行了数值计算.在此基础上,得出了一些有意义的结论.     

后掠机翼欧拉方程数值模拟

本文对后掠机翼用保角转绘加剪切变换生成C-H型网格,并在这类网格拓扑上,用Jameson的三维欧拉方程有限体积法,四步Runge-Kutta时间推进格式,研制出可供分析三维后掠机翼亚、跨、超声速绕流的计算程序。本方法的特色是改进了机翼表面网格点的分布,使机翼后缘与网格线一致。在数值计算中,通过附面层位移厚度计算,引入了附面层粘性修正,改善了计算结果。     

圆转方非轴对称喷管无粘流动的多层网格法求解

本文利用多层网格法求解三维无粘亚音和跨音速圆转方非轴对称喷管的流动,基本格式是把Ni[1]的二维有限体积积分格式推广到三维流场。边界条件采用文[2]中的“预测一修正”处理方法。为了进一步提高格式稳定性,在计算流动变量的一阶时间变化量时,采用单元体积平均法而不是文[1]中的算术平均法。引入当地时间步及多层网格方法,以加快格式的收敛速度。文中还对多层网格构造形式及其收敛效果作了讨论。通过对两种喷管模型流动的计算和与试验数据的对比,验证了该算法的精度。     

可靠性及可靠性灵敏度分析的改进点估计方法

针对低维高度非线性问题,提出了一种可靠 性及可靠性灵敏度分析的改进的点估计方法。基本思想是首先由空间分割来降低局部子空间 中功能函数的非线性程度,然后用低精度的稀疏网格积分探索子空间中功能函数的概率响应 特性,最后组合子空间中的信息来得到所需的可靠性及其灵敏度分析结果。方法的优点 是适用于低维高度非线性功能函数的失效概率和可靠性灵敏度分析,由于无需求解功能函数 的梯度函数,因此适用于复杂的隐式功能函数。另外,由于方法利用少量均匀抽样来估计对 失效概率贡献最大的点,并依据其进行子空间的划分,从而使得子空间的划分更有利于提高 可靠性及可靠性灵敏度分析的效率。用算例对所提方法进行了验证,结果表明：在低维高度 非线性条件下,所提算法的精度和效率比类似的三点估计、直接稀疏网格积分方法有明显优 势。     

涡格法中二等分平板的升力和力矩的高精度解

本文从点涡和控制点位置的选取出发,研究涡格法中二等分二维平板的升力和力矩的高精度解。发现除传统的1/4和3/4法则之外,还有服从一定规律的无穷多种离散化格式存在,均可以得到精确的升力和力矩。这就为在保证升力和力矩精确的前提下,寻找优于传统涡格法的载荷分布提供了理论依据。     

任意回转面叶栅C网格的数值生成

首先尝试采用任意曲面椭圆方程，数值生成任意回转面上叶栅Ｃ网格，经过对边界条件的仔细处理，生成的网格可保证在周期外边界和重合内边界上网格线导数连续，在叶片表面网格线相互垂直，提出了结合初始网格线分布，控制环绕叶片网格线分布的方法。     

环形激波聚焦流场特性的间断有限元方法(英文)

针对环形激波聚焦时产生的高温、高压特性,采用间断有限元方法对环形激波在同轴圆柱形激波管内的聚焦流场进行了数值模拟。为了验证采用方法的有效性,首先对经典激波管问题进行了求解,计算结果与解析解吻合很好。其次,对入射激波马赫数为2.0到5.0的环形激波聚焦流场进行了计算,获得了环形激波在圆柱形激波管内聚焦过程的气动特性。模拟结果表明,采用间断有限元方法能够有效地捕捉激波聚焦过程形成的绕射、反射和聚焦等主要流动特征,而且在聚焦点附近,数值解具有较大的梯度变化,表明该方法对间断解具有较强的捕捉能力。最后比较了网格尺度对聚焦点峰值压力的影响。     

某飞行器声信标切削轨迹设计

声信标为某飞行器锥段壳体的重要设备,针对安装位置与理论位置存在较大偏差,导致切削轨迹无法精确生成,对此,本文提出集测量、编程、切削、检测为一体的切削轨迹设计方案。首先在锥段壳体设定某基准位置与声信标特征相同,然后对声信标空间位置进行测量,将锥段壳体基准位置与声信标测量位置切削轨迹通过坐标系建立联系,为此设计出声信标空间位置控制点和方向矢量计算方法;为验证切削轨迹的正确性,提出优化后处理,约束设备旋转角度,校验切削轨迹特性。最后,经试验验证,声信标空间位置控制点与方向矢量和测量数据一一对应,切削轨迹与声信标特征大小相符,该方法合理、可靠、有效。     

粘性非结构网格生成方法及绕流计算

运用非结构网格进行计算时,需要在物体附面层区域内生成高展弦比的三角形网格。本文给出了一种基于物理学中弹簧原理的粘性非结构网格生成新方法－－挤压法,它能有效地实现粘性区域内位伸网格到附面层外一般非拉伸网格的光滑过渡。然后在非结构粘性网格上对Ｎ－Ｓ方程运用有限体积法进行了数值计算,ＮＡＣＡ－００１２翼型以及平板层流的数值试验结果较为满意。     

独立电源系统有源滤波器测控系统自动定制技术

采用自动定制技术设计了独立电源系统有源滤波器测控系统。对独立电源系统有源滤波器测控系统进行了需求及特点分析。以需求为出发点并参考特点建立测控系统的体系结构,进行功能性的划分,实现了测控系统的功能构件化。采用面向对象的设计思想,将组成系统的构件抽象化为类,利用UML建模语言对测控系统进行建模,利用代码生成工具以及成熟的硬件构件实现测控系统。将自动定制得到的测控系统实际应用于独立电源系统有源滤波器中,对该系统进行了检验。