一类二阶守恒单调重映算法

在大变形流体力学问题的数值模拟中，经常会涉及到计算网格的重分。基于不同网格的物理量传递便是所谓的重映技术。重映算法是任意拉格朗日一欧拉方法的重要组成部分，本文描述了一类适用于任意网格的二阶守恒单调重映算法。该算法分为网格内物理量的多项式重构、近似积分计算、物理量的单调修正三个部分。本文采用了最小二乘法的思想构造了网格内物理量的梯度，并且通过对物理量的单调修正，保证了算法的精度和单调性。该算法计算公共表面上的网格间质量的改变，所以可以在任意网格上使用。通过一系列数值算例验证了该算法，并说明本文给出的算法是有效可行的。     

区域法数值模拟短突扩压器流场

在一般曲线坐标下对环形燃烧室突扩压器和其内外环通道流场进行了数值模拟。紊流特征有用工程上常用的κ－ε双方程模型来描述，使用壁面函数法处理曲线壁面边界，采用非交错网格系统和ＳＩＭＰＬＥ算法求解各守恒方程。为了计算方便，利用区域法对复杂的流道进行分块处理，计算所得的速度和压力分布较为合理，表明本方法是可行的。     

三维曲线坐标系N－S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐标系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错网格设置，该算法以ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的出现，在压力修正方程中对压力梯度的离散附加一项１—δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

二维非结构网格DSMC方法的实现及其应用

研究了非结构网格DSMC方法实现的过程，将Bird位置元方案中的子网格思想引入到非结构网格上来，只存储子网格的总体标识号，做到了利用较少的计算网格提高分子的分辨率和计算精度；提出了将面积元坐标搜索算法与交替数字二叉树搜索算法(ADT)相结合来跟踪模拟分子在网格之间的迁移；使用ADT方法判别分子与物面是否相交，节约了计算时间；最后对过渡流域高超声速绕流进行了数值试验．结果表明这种方法的可行性、高效性及算法的通用性。     

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

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

HLLC方法在三维非结构网格上的应用

以三维非结构网格的显式有限体积法为基础，采用格心方法以及基于近似黎曼解的二阶精度Godonov扩展方法求解三维Euler方程，使用HLLC(Harten，Lax，van Leer,Contact)近似黎曼解的方法计算网格单元边界处的守恒量通量。为了验证方法的可行性，利用三维解算器对NACA0012翼型绕流进行了数值计算，结果令人满意。在此基础上对某型弹的流场进行了数值模拟，取得了比较好的结果。     

非正定矩阵对磁共振扩散张量导出量的影响

为研究非正定矩阵对水分子扩散张量导出量的影响,提出了一种计算扩散张量导出量的优化算法。采集5位健康志愿者脑部磁共振扩散张量图像,沿25个方向施加强度为40mT／m的扩散敏感梯度磁场,扩散权重值为1000s／mm^2,求解超定方程,得到了各体素的二阶扩散张量。在灰质区域非正定矩阵较少,在胼胝体的膝部和压部发现大量非正定矩阵。由于这些非正定矩阵的影响,使用常规算法得到的平均扩散率和各向异性分数分别被低估和高估。为此修正了常规算法,使用矩阵特征值的绝对值计算扩散张量导出量。结果表明,由修正算法得到的扩散张量导出量图像的鲁棒性及可靠性均得到提高。     

三维曲线坐标系N—S方程通用形式和数字解

运用张量分析理论，比较简捷地导出了三维任意曲线坐系上椭圆型Ｎ－Ｓ方程的通用形式。在对这些方程进行数值离散时，采用了有限控制体积的非交错风格设置，该算法双ＳＩＭＰＬＥ为基础，为了抑制压力振荡场的了出现，在压力修正方程中对压力梯度的离散附加一项１－δｘ的差分。作为算例，本文给出了三维９０°强曲率弯曲管道中气流的不可压缩紊流状态的数值计算结果，并与相关的实验数据进行了对比，结果是满意的。     

用非交错网格求解曲线坐标系下的流动控制方程

本文采用非交错网格法,求解任意曲线坐标系下的流动控制方程。文中阐述了压力振荡的根本原因和相应的消除措施,导出了压力修正方程中控制体边界上曲线速度分量的修正表达式。在这些表达式中,出现了一个表征一个网格间距和两个网格间距压力梯度差分值的附加差值项。当出现压力振荡时,该差值很大,可以有效地消除振荡,而当压力场没有振荡时,这一差值又很小。作为计算方法可靠性的检验,本文分别计算了二维直通道和渐缩形混合管内的气流流动问题,结果是令人满意的。     

跨洞庭湖大桥主梁节段定常气动力研究

采用风洞实验与数值模拟相结合的手段，对跨洞庭湖大桥主梁节段的定常气动力特性进行了研究；流场数值模拟是采用ｋ－ε湍流模型封闭的雷诺时均Ｎ－Ｓ方程，求解方法是建立在交错网格基础上的压力－速度场迭代求解的ＳＩＭＰＬＥ算法；实验是在国防科技大学ＫＤ－０３ 低速风洞进行；研究表明实验与计算结果较为一致     

尾支杆干扰非线性修正方法的初步研究

利用CFD方法 ,采用两套网格重叠的形式 ,研究了一种对尾支杆干扰进行修正的非线性方法。运用本方法 ,研究了尾支杆等直段长度、直径和锥角对模型气动特性的影响 ,研究了四个模型有攻角时的尾支杆干扰 ,并对GBM 0 4A模型的零阻进行了尾支杆干扰的修正。     

杯芳烃及其合成方法

本文评述了杯芳烃的合成方法，包括一步法和多步法，并对其研究前景进行了展望。     

非正交曲线坐标湍流族流的数值模拟

对锥形燃烧室的冷态流场进行了数值模拟，采用Ｔｈｏｍｐｓｏｎ的非正交贴你坐标系统，用坐标变换的方法处理锥形燃烧室曲线壁面边界，用ＳＩＭＰＬＥ算法求解连续方程和Ｎ－Ｓ方程。亲流模型采用工程上常用的ｋ－ε模型。还运用贴体坐标系统数值研究了不同旋流数和不同进口紊流动能对湍流旋流流动的影响，并与实验结果相比较获得较满意的结果，对新型燃烧室的设计提供了理论依据。     

COMPUTATION OF UNSTEADY AERODYNAMIC FORCES ON WINGS WITH NAVIER-STOCKS EQUATIONS

INTRODUCTIONFor the requirement of instantaneity of un-steady flows,the minimal time step is needed,which leads to a requirement for a large numberof time steps.To an explicit method,it marchesvery slowly and many of the acceleration tech-niques used for steady flows cannot be used be-cause they destroy the time accuracy.The“dual-time” method allows an implicit discretisation tobe used in real time,but at each real time stepmarches the solution in a pseudo time to steadystate,through an …     

医学CT的误差修正方法、图像重建及实现

医学Ｘ光ＣＴ作为一种诊断设备，对其图像质量要求很高，而在实际应用中各种误差都将会影响最终的图像质量，因而对误差进行修正是非常重要的，这一点在各厂家的ＣＴ 产品中表现得尤为突出。通过研究对医学Ｘ光ＣＴ 的探测器、Ｘ光射束硬化效应等所产生的误差进行修正，给出了在ＰＣ机上对误差进行修正和图像重建的方法和步骤，完成了对分辨率水模、对比度水模、病人头部和肺部的图像重建，证明该方法行之有效。     

带螺旋桨飞机模型风洞实验洞壁干扰修正

带动力实验进行洞壁干扰修正时，必须考虑螺旋桨滑流的影响。运用等效动压法和壁压信息法对运八飞机带螺旋桨模型风洞实验数据进行洞壁干扰修正，分析洞壁对带螺旋桨飞机模型试验数据的影响并对不同修正方法进行了比较。两种修正方法的修正结果基本一致，壁压信息法能实际反映洞壁干扰影响，但壁压信息法需要进行准确的壁压测量，增加实验工作量；等效动压法是以经典的洞壁干扰修正公式为基础，考虑了螺旋桨滑流的影响，而且带螺旋桨飞机模型实验的攻角一般在中等攻角以下，不会出现像大攻角那样的分离特性，因此这种方法有较好实用性。     

2.4m跨声速风洞压敏漆测量系统研制与应用研究

近十几年来,由于压敏漆(Pressure Sensitive Paint、PSP)测量技术的不断完善与发展,国际上主要空气动力试验机构逐步将其应用于2 m量级工程型风洞,完成模型表面压力测量、模型表面流动显示与 CFD 结果验证。在2.4m跨声速风洞建立了双组份、多光源和多CCD的PSP测量系统,解决了大型暂冲式跨超声速风洞试验存在的模型表面温度变化、光照均匀性与强度变化,以及模型振动、试验数据修正、喷涂与压敏涂料校准等诸多影响PSP测量结果精准度与可靠性的问题,并成功应用于大飞机测压模型和三角翼测压模型压力分布测量试验。试验结果表明:在小迎角范围压敏漆涂层对模型表面压力分布影响不明显;在试验马赫数0.4~0.82、模型迎角-4°~4°范围,PSP与传统电子扫描阀测量结果的Cp 均方根偏差小于0.03,测量精准度与国外同量级连续式跨声速风洞相当。可以为飞行器气动优化设计和空气动力学研究提供一种新的、先进的测试技术。     

High-Order Discontinuous Galerkin Method for Hovering Rotor Simulations Based on a Rotating Reference Frame

An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. Instead of auxiliary methods like grid adaptation,higher ? order simulations(fourth ? and fifth ? order accuracy) are adopted.Rigorous numerical experiments are carefully designed,conducted and analyzed. The results show generally excellent consistence with references and vigorously demonstrate the higher?order DG method's better performance in loading distribution computations and tip vortex capturing, with much fewer degrees of freedom(DoF). Detailed investigations on the outer boundary conditions for hovering rotors are presented as well. A simple but effective speed smooth procedure is developed specially for the DG method. Further results reveal that the rarely used pressure restriction for outlet speed has a considerable advantage over the extensively adopted vertical speed restriction.     

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

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