高精度法矢下切割面自适应的凸曲面射线寻迹

针对电大尺寸目标难以精确解析表达带来的一致性几何绕射算法应用难问题,提出了基于三角网格、适用于任意凸曲面的射线寻迹（TM-tracing）算法。应用工程中较易获取的三角网格及其协议,设计了一种满足快速多边搜索条件的网状数据存储链表;提出了满足寻迹要求的高精度法矢求解算法;采用切割面自适应调整的弧形拟合寻迹方法实现了爬行波寻迹算法;结合一致性几何绕射理论（UTD）实现了暗区场值求解算法。任意网格曲面射线寻迹结果表明:本文提出的寻迹算法适用于包括球、柱和锥在内的任意光滑凸曲面,寻迹偏差小于1.61%,寻迹速度为2.8 s,具有一定的工程应用价值。      

基于频域LMS的自适应波束形成算法

在分析传统自适应波束形成的基础上,提出了一种基于频域最小均方(LMS)的自适应波束形成算法。该算法先对输入信号进行FFT变换,再通过LMS算法实现了频域上自适应波束形成。FFT变换后信号为稀疏矩阵,自相关下降,LMS算法收敛速度提高;理论分析和仿真结果表明了该算法收敛速度较快、性能较好,且计算量较少,易于实时实现,而且文章提出的波束形成算法对相干信源具有鲁棒性。     

结合等值面绘制与体绘制的电磁环境可视化方法

为提高现有的电磁环境三维可视化效果,提出一种将等值面绘制融合到光线投射体绘制中的网格投影算法。将等值面网格化,利用空间网格表现等值面并融合进体绘制生成的电磁环境图像中。将平面均匀网格平行投影到电磁环境体数据中具有同一数据值的数据采样点上,由此生成网格化的等值面。投影到体数据中的网格沿投射光线的反方向确定其对应于光线投射生成图像上的像素点,并将网格颜色值融合到此像素点中,实现融合绘制电磁环境。在统一计算设备架构(CUDA)下并行实现本文算法后,不仅电磁环境的绘制效果比融合前有提升,而且更新绘制速度能够达到实时,支持用户实时交互控制。     

非结构重叠网格显式装配算法

针对重叠网格中洞映射法占用过多物理内存的问题,发展了一种改进型洞映射法;基于相邻单元搜索法,发展了一种基于相邻阵面的贡献单元搜索法;通过将割补法与隐式切割技术相结合提出了一种非结构重叠网格显式装配算法。该算法首先生成一套包围物面的笛卡儿网格,其次存储所有与物面边界相交的笛卡儿网格信息,最后根据所存储笛卡儿网格与所需判断的网格单元的相对位置来判断其是否为洞内单元。在成功判断出所有洞内单元后,以当前洞边界为初始阵面推进,同时以各个网格单元的物面距离为判别标准对重叠区域进行优化,生成最终插值边界。所提算法优化了传统非结构重叠网格装配过程,具有物理内存占用低,贡献单元搜索次数少以及计算效率高等特点。通过2个典型复杂流动算例验证了所提算法的准确性与适用性。      

混合重叠网格插值方法的改进及应用

混合重叠网格间合理插值是保证流场计算正确的基础之一。本文针对重叠区网格尺寸匹配度较差时插值误差较大问题,发展了一种新型混合重叠网格插值方法,通过使用二阶精度插值、按单元尺度区分和扩充模板,改善了插值区网格尺寸匹配度较差时的插值精度。本文方法适用于任意单元类型的混合网格重叠,对各类单元处理透明,实现简单。计算结果表明,采用本文方法在网格重叠区流场变量传递正确有效,插值区网格尺寸匹配度较差时,相比原始方法,等值线过渡更为光滑,变量经过插值区耗散更小,计算与试验值符合更好。     

高超声速热流高精度数值模拟方法

通过高超声速下钝锥热流的计算,对影响热流计算精度的因素进行了综合研究.获得了壁面法向网格Re数对热流计算的影响规律.采用边界高阶插值,提高了热流计算精度.用Roe和AUSM(Advection Upstream Splitting Method)+格式搭配minmod和混合limiter两种限制器,探讨了空间离散方法及限制器的耗散性对热流结果的影响.应用五阶WENO(Weighted Essentially Non-Oscillatory)方法,对高阶格式计算热流的性能进行了细致研究,最终确立了网格依赖小,热流精度高的计算方法.研究认为:高阶格式不但可以放宽热流计算对网格Re数的限制,同时使得格式对本身计算方法的依赖性也变小.      

基于云模型的自适应并行模拟退火遗传算法

针对遗传算法收敛速度慢,容易"早熟"等缺点,提出了一种改进的遗传算法,即基于云模型的自适应并行模拟退火遗传算法(PCASAGA,Adaptive Parallel Simulated Annealing Genetic Algorithms Based On Cloud Models).PCASAGA使用云模型实现交叉概率和变异概率的自适应调节;结合模拟退火避免遗传算法陷入局部最优;使用多种群优化机制实现算法的并行操作;使用英特尔推出的线程构造模块(TBB,Threading Building Blocks)并行技术,实现算法在多核计算机上的并行执行.理论分析和仿真结果表明:该算法比其他原有的或改进的遗传算法具有更快的收敛速度和更好的寻优结果,并且充分利用了当前计算机的多核资源.      

基于二维排序和平面加厚技术的自由曲面消隐算法

消除隐藏线是计算机真实感显示的一项重要，自由曲面的消隐则是其中的难点，本文提出的自由曲面消隐算法提出了平面加厚技术和二维排序技术，给出一块三角曲面片，加厚技术就是用加厚三角形所得的三角砖来逼近曲面片，采用二维排序技术，光栅点的可见性判定归结为矩形域在排序表中的查询，测试表明，本算法较传统算法提高速度两个数量级以上，消隐时间随光栅点线性增长，光栅点最多可达６２０００个（对于６４０Ｋ内存）。     

网格结构对Euler方程解的影响

 采用2种不同对接分区网格技术讨论了不同的网格结构对绕偏副翼的三维机翼流动的Euler方程解的影响.一种是部分不连续的网格(PDG);另一种是完全连续的网格(CCG),分别在PDG和CCG中采用Van Leer的矢通量分裂格式和Jameson的中心差分格式求解了Euler方程.算例表明,2种网格各有其优缺点,且都能得到有用的结果,但为了高效地取得更好的计算结果,必须根据流动和几何外形选取恰当的计算网格.     

基于CE-PF算法的舰载机离场调度优化问题

甲板作业调度研究是提升航母战斗力的关键技术,而其具有时间、空间与资源受限的复杂约束调度问题已被证实为NP-hard。根据舰载机出动离场调度优化问题的特点,将其抽象为零缓存区混合流水车间调度模型,建立包含飞机避碰等约束的混合整数规划模型。提出一种交叉熵与作业剖面匹配(CE-PF)算法用于问题求解,并给出了算法流程架构。交叉熵算法通过高斯采样完成启发式规则下的工件分组,作业剖面匹配算法完成分组工件的任务排序、作业编排及约束检查等调度设计,Gap逼近算法进行目标值评估、精英种群选择、抽样参数更新及收敛判定。通过算例仿真,验证了CE-PF算法求解离场调度优化问题的有效性;灵敏度分析表明起飞模式和空间约束对出动效能影响较大。      

Euler方程的分裂型通量分裂双时间步隐式方法

传统隐式方法有格式复杂、计算量大等缺点,在Euler方程的差分离散过程中,利用算子分裂思想,结合通量分裂法、双时间步法等隐式离散方法,构造了一种更简单的分裂型隐式计算方法.通过对典型空气动力学问题的计算,检验了该方法的有效性和可靠性,并对其性能做了具体讨论.该方法具有稳定性好、时间步长约束小等隐式格式的普遍优点,同时具有格式简单、程序易实现等优点;避免了传统隐式方法单步推进时的方程组常规求解及矩阵求逆过程,计算量小;比LU-SGS方法收敛速度快.      

超声速进气道可压及不可压流动数值模拟

采用Chio-Merkle预处理矩阵对可压NS(Navier-Stokes)方程组进行时间预处理,分析了预处理方法的物理和数学背景.用有限体积方法,结合LU-SGS(Lower-Upper-Symmetric-Gauss-Seidel)隐式时间积分和AUSM+(P)(Precondition Advection Upstream Splitting Method)格式、中心差分2种空间离散格式,求解该预处理NS方程组.通过圆弧凸包流动和二维方腔顶盖驱动流动的数值试验表明,该方法克服了传统时间迭代方法模拟低速流动时的刚性问题,加速了收敛过程,可以同时有效地模拟可压及不可压流动.用此方法与块结构化网格技术相结合,进行了在不同飞行马赫数、攻角和出口反压条件下,混压式轴对称超声速进气道三维流场数值模拟.计算结果表明:该方法能准确地捕获复杂激波系,清楚地揭示进气道三维流动现象;进气道性能随工作条件的变化规律,与理论分析一致.      

不同厚度翼型动态失速涡运动数值研究

在低马赫数下, 对三种不同相对厚度的NACA系列基本厚度翼型在俯仰振荡运动中的动态失速现象进行了数值研究.数值模拟时采用双时间法和LU-SGS隐式解法相结合的模式求解了非惯性坐标系下的非定常纳维-司托克斯(NS)方程组, 空间离散采用Roe格式,并结合刚性动网格生成技术.采用全湍流计算,通过引入BL(Baldwin-Lomax)模型计入湍流影响.计算出的气动力迟滞曲线与实验结果变化趋势符合较好,表明了数值方法的有效性,同时通过比较和分析不同厚度翼型在轻失速和深失速下的流场结构,发现绕翼型的失速旋涡产生和发展规律是明显不同的.      

刚性方程的二阶Gear方法的注记

构造了一个求解刚性方程的二阶差分格式,它是隐式的,并采用简单迭代求解,对于稳定系统,迭代收敛性条件与步长无关,而对不稳定系统,迭代收敛性条件对步长却有所限制,利用这一限制性条件去克服解刚性方程的二阶Gear方法时出现的“危险性”。     

球坐标系六片网格下三维定态行星际太阳风模拟

采用二阶MacCormack差分格式, 利用稳态的磁流体(MHD)方程组在球坐标系六片网格下模拟研究了行星际太阳风. 六片网格系统能有效避免极区奇性和网格收敛性. 迭代按径向方向推进求解, 很大程度上减少了计算量, 节约了计算时间. 内边界条件根据太阳与行星际观测确定, 比较测试了5种内边界条件, 模拟给出了1922卡林顿周的背景太阳风结构. 几种内边界条件所得模拟结果与行星际观测基本吻合. 太阳风速度采用McGregor 等的经验公式给出, 磁场由水平电流片(HCCS)模型得到, 密度和温度分别根据动量守恒和气压守恒得到, 研究表明采用这样的边界条件模拟结果最佳.      

3点逐步r型自适应网格算法

为了根据某个流场参数梯度较合理地分布网格同时又不增加网格点的数量,提出一种基以总变差为基础的r型3点逐步自适应网格算法,通过3点逐步调整网格点的位置以减小网格线的扭曲,生成质量较好的计算网格.以绕二维圆柱和三维双椭球的超声速流动为算例,用Euler方程模拟,以流场的压力总变差来进行自适应,通过分别用所提出的自适应方法生成的自适应网格与用原始网格计算所得的流场结果对比,计算结果表明,用自适应网格计算所得的激波比用初始网格计算所得的激波薄,采用自适应方法所计算出的流场数值解具有更高的分辨率.