温度体动网格模型中控制参数的研究

采用调整导热系数和源项分布的方式,对温度体动网格模型所生成的网格品质进行控制.温度体动网格模型将运动边界的位移虚拟为求解域的温度边界条件,以流体能量方程或固体导热方程作为控制方程,通过选取不同的导热系数或源项分布可以得到适合不同求解问题需要的温度分布,并将求解得到的网格节点的温度作为其动态位移量.采用温度体动网格模型计算平动、三维旋转运动、柔性体运动等流固耦合计算中可能涉及的动网格算例.结果表明,在计算效率和生成网格的品质方面,与已有的动态网格生成方法相比,具有较大的优势,可以有效解决流固耦合数值模拟中存在的动态网格生成难题.     

贴体坐标网格生成技术的研究

本文对数值求解微分方程法生成贴体坐标网格的技术进行了探讨和尝试。文中不仅对数值求解双调和方程的两种不同方法作了介绍,还利用该方程生成了典型的Ｈ型、Ｃ型、Ｏ型网格。结果表明,双调和方程法在控制网格线间距及与边界的正交性方面具有很强的能力,所生成的网格质量令人满意,且可用于复杂的物理区域。因而,采用双调和方程生成贴体坐标网格是流场数值模拟中一种有效的、高质量、具有广泛通用性的网格生成技术。     

移动网格与Level Set耦合方法在气-液两相流数值模拟中的应用

本文采用自适应移动网格与Level Set函数相耦合的方法来实现气-液两相流的数值模拟与计算.作为自适应网格方法的一种,移动网格方法主要是为了解决发展方程的计算问题而设计的方法.文中给出了移动网格的生成方程,并针对方程的非线性,给出了一种半隐式的离散方法用于进行求解.本文将移动网格方法与Level Set方法相耦合,将控制流体运动的Navier-Stokes方程以及追踪相界面的Level Set方程转换到曲线坐标下,应用一套曲线坐标方程组来同时描述气、液两相流的运动规律,成功实现了对气-液两相流问题的数值模拟.通过对顶盖驱动流的计算以及对液滴沉降现象的模拟计算,验证了本文方法的可靠性.本文对常重力与微重力下两气泡融合的发展规律进行了数值模拟,通过分析对比,得到了重力对两气泡融合变形的影响规律.     

并行重叠/变形混合网格生成技术及其应用

为了适用于柔性变形、相对运动等复杂动边界问题,建立了并行环境下重叠和变形相结合的动态混合网格生成技术.通过计算区域分解以及分布式并行实现了重叠和变形技术的结合,其中重叠网格采用了并行化的隐式装配方法,并发展了两种并行化查询策略.变形网格则采用了并行化的径向基函数（RBF）插值方法.并行化动态网格生成方法大幅提高了动态网格生成效率,有利于处理大规模的动边界问题.在此基础上,发展了基于变形/重叠动态混合网格的流动/运动/控制一体化数值模拟方法,进一步改进了耦合模拟软件平台——HyperFLOW.典型应用算例证明了该动态混合网格技术及一体化算法的实用性.      

结构型网格分块生成技术

构造了一种适用于模拟具有复杂几何特征的流动问题的结构型网格生成方法,该方法应用CAD技术中非均匀有理B样条(NURBS)曲面表达形式,根据给定物体的典型站位截面,首先对物体进行分片表面网格的构造与重构,同时根据物体几何特征及流场特性对绕流流场进行区域分解并构造其块边界,在此基础上再采用数值方法进行相应的分块网格生成,并给出了应用所建立的方法生成的几种天地往返运输系统流场网格.     

基于车身外形CAD数据的表面氯动网格生成

提出一种在非均匀参数空间上生成表面数值网格的方法,开发了相应的表面网格生成软件ＳＧＲＩＤ,并应用于汽车动力学计算机仿真中具有复杂几何外形的车射表面网格生成。     

基于结构化动网格的涡旋压缩机内部流场模拟

为了增加涡旋压缩机啮合间隙处的动网格数量。本文根据涡旋压缩机的流体区域变化规律,提出一种结构化动网格的生成方法,实现了啮合间隙处的网格加密。进而得到一种基于结构化动网格的涡旋压缩机非定常流动的数值模拟方法;得到了任意时刻流体区域内的流场分布。与现有的非结构化动网格进行比较,结构化动网格解决了非结构化动网格在啮合间隙处易产生负体积的问题,能够在啮合间隙处生成足够的网格数量,提高了计算精度;且具有网格质量高、网格数量和疏密易控制等优点,更适合于涡旋压缩机工作过程的数值模拟计算。     

基于CFD/CSD的舵面控制气动弹性数值模拟方法

在动态网格上通过耦合求解流动控制方程和结构动力学方程, 发展了一种舵面控制下飞行器运动响应过程中气动弹性数值模拟研究方法.流动控制方程采用N-S方程, 结构动力学采用线性模态叠加方法, 其中流动控制方程空间离散采用基于非结构网格的有限体积方法, 对流通量采用计算HLLC格式, 非定常时间离散采用基于LU-SGS的双时间步长方法.模拟中, 气动运动和结构变形在双时间步长方法推进过程中采用改进松耦合方法, 气动网格与结构网格之间信息交换采用无限平板样条法实现, 飞行器的运动和变形采用基于重叠网格和Delaunay图映射变形网格相结合的方法进行处理.采用多个考核算例对发展的数值方法进行考核验证, 结果表明该方法可以高效精确模拟舵面开环控制下飞行器运动响应过程中的气动弹性特性.      

一种新的非结构四边形网格生成方法及软件

提出一种间接生成非结构四边形网格的方法.采用Delaunay三角形网格生成的算法生成三角形单元,在三角形单元中定义生成四边形网格的前沿,根据三角形的特点(边、角大小)以及生成四边形前沿的方向,将两个三角形合并为一个四边形单元.此方法的特点是不增加节点,仅改变点与点之间的拓扑连接,使单元基本变为四边形占优的单元.然后研究四边形占优的网格转换为100%的四边形网格的转换模版.并应用于实际物理问题.数值实验显示了该方法生成非结构四边形网格的特点.     

用非正交曲线坐标系数值计算航空发动机加力燃烧室无化学反应的湍流流场

本文采用非正交曲线坐标系下非交错网格的SIMPLE方法对航空发动机加力燃烧室无化学反应的湍流流场进行了数值计算，湍流模型采用k-ε双方程模型．差分网格采用分区方法生成，计算时对整个流场进行分区迭代直至得到收敛结果．加力燃烧室湍流流场数值计算结果合理．     

用非结构网格与欧拉方程计算复杂区域的二维流动

提出用Delaunay三角化方法生成非结构网格的一种过程。所生成的网格可用于复杂多连通域内的可压流计算。采用Euler方程和格心有限体积法,研制出程序,给出了算例。     

弹性力学的复变量无网格方法

在移动最小二乘法的基础上，提出了复变量移动最小二乘法.复变量移动最小二乘法的优点是采用一维基函数建立二维问题的逼近函数，所形成的无网格方法计算量小.然后，将复变量移动最小二乘法应用于弹性力学的无网格方法，提出了复变量无网格方法，推导了复变量无网格方法的公式.与传统的无网格方法相比，复变量无网格方法具有计算量小、精度高的优点.最后给出了数值算例. 关键词： 移动最小二乘法 复变量移动最小二乘法 无网格方法 弹性力学 复变量无网格方法     

Spectral difference method for compressible flow on unstructured grids with mixed elements

This paper presents the development of a 2D solver for inviscid and viscous compressible flows using the spectral difference (SD) method for unstructured grids with mixed elements. A mixed quadrilateral and triangular grid is first refined using one-level h-refinement to generate a quadrilateral grid while keeping the curvature of boundary edges. The SD method designed for quadrilateral meshes can subsequently be applied for the refined unstructured grid. Results obtained with the SD method for both inviscid and viscous compressible flows compare well with analytical solutions and other published results.     

光电成像探测中目标方位测量方法

为了在光电成像探测目标的同时给出目标相对探测器的方位信息,探讨了一种基于单目视觉的运动目标方位测量方法。根据投射成像原理,推算出二维目标像点坐标与目标的空间三维位置之间的映射关系,建立了单目视觉测量目标方位的数学模型。结合帧间差分法和KLT方法检测静止背景下和变化背景下的运动目标,并对目标特征点进行亚像素定位。试验和仿真结果分析表明,该方法能够有效提取目标的特征点,可对目标较为准确地定位,并通过单目视觉测量出目标方位信息,误差控制在8%的范围内。     

二维非线性抛物方程广义差分法/有限体积法的自适应计算

给出三角网上二维非线性抛物方程广义差分法(有限体积法)的一种基于残量估计的后验误差估计,并在此基础上设计了自适应计算方案,以适应物理解在时空的大梯度变化.提出了适合发展方程自适应计算的三角网数据结构(不是树状结构)和灵活的局部粗化算法.     

Vector boundary-element analysis for conical horn

A vector boundary-element method is developed for analyzing a conical horn in this paper. A moving boundary technique (MBT) is proposed for solving difficulties of the integration at high-order singularity points in a vector integral equation, and its validity is verified by calculating the reflection coefficients of an empty circular waveguide with a short-end. The circular symmetry of the conical horn is utilized to simplify the 3-D problem to 2-D one and the analyzing efficiency is significantly improved. A comparison of the reflection coefficients of the conical horn is given between the experimental data and the calculated results. Very good agreement has been seen, and the reliability and efficiency of the present method are thus verified.     

用墨浸法自动生成三维无结构网格

提出了为用Delaunay法自动生成无结构网格而发展的"墨浸式"加点法。采用"墨浸式"加入边界点的方法,使得Delaunay方法的固有缺陷,即可能出现边界穿透的问题,得到简化。"墨浸式"加点法成功地在一些复杂边界域内生成了可用于流场计算的无结构网格。     

多点择优推进阵面法生成曲面三角形网格

在现有曲面非结构网格生成法的基础上,提出了一种新的曲面网格生成法——多点择优推进阵面法.它可在曲面上直接进行三角形网格划分,克服了映射法的网格变形问题,并且可以在网格生成结束后,对曲面网格直接进行Laplace格点松弛光顺.该方法使用简单,不受曲面块类型的限制,且网格质量高,可以为三维非结构网格生成提供高质量的初始阵面,并给出了若干个算例.     

基于非结构变形网格的间断装配法原理

在激波捕捉法计算得到的流场基础上采用辨识算法得到初始间断位置, 从ALE方程出发, 考虑离散几何守恒律, 采用变形网格和网格重构技术解决计算过程中间断运动和变形, 新旧网格之间流场采用高精度信息传递方法保持时间精度, 建立了基于非结构动网格技术的间断装配方法.通过激波管问题的二维模拟, 模拟了初始间断分解为激波和接触间断激波遇到固壁反射后与接触间断相交的非定常流动过程, 对这种新方法的基本原理进行了介绍.      

非结构混合网格上的NS方程求解方法

提出了一套较为通用的,完全自动化的非结构混合网格生成方法.在物面粘性作用区,采用一种改进的推进层方法生成三棱柱形和金字塔形网格;在其他流动区域采用阵面推进方法生成四面体网格.采用一种改进精度的格心有限体积法对三维NS方程进行了求解,在加速收敛措施方面,提出了一种新的当地时间步长取定方法来减小质量较差的网格单元对流场计算稳定性和收敛速度的不利影响.以M6机翼和DLR/F4翼身组合体外形的粘性流场作为数值算例,验证了上述网格生成和流场求解方法的正确性和实用性.