等离子体电波传播及散射问题的CDRC-ADI-FDTD方法

首次把交替方向隐式技术(AD I)与等离子体的电流密度递归(CDRC)卷积技术结合,给出了碰撞非磁化等离子体的CDRC-AD I-FDTD方法。推导了碰撞非磁化等离子体中的二维CDRC-AD I-FDTD迭代公式,并用算例验证了碰撞非磁化等离子体CDRC-AD I-FDTD算法也是无条件稳定的。计算结果表明,等离子体CDRC-AD I-FDTD算法与传统的FDTD方法和等离子体JEC-FDTD方法的计算结果吻合,计算效率更高。     

超声脉冲波在固体中传播的有限差分数值仿真及动画演示

文章基于Matlab软件包,运用有限差分方法模拟演示超声脉冲波在固体介质中的传播和散射。对各向同性固体弹性波传播的声波方程的离散交叉配置格式进行了分析,列出了模型的吸收边界条件和自由边界条件。通过对空间间隔长度和时间步长参数的选取来控制数值计算结果的稳定性,给出了形成动画演示的处理过程。以角束探头激励产生超声脉冲爬波的实例展示了动画演示的效果,声场的数值计算结果与理论分析也相符合。     

武装直升机雷达散射截面计算及雷达吸波材料影响分析

为提高直升机雷达散射特性预估的准确性,建立了目标雷达散射特性分析的计算电磁学(Computational electromagnetics method,CEM)方法,并开展了吸波涂层对直升机雷达散射截面(Radar cross section,RCS)特性影响的研究。首先,对复杂目标(例如直升机)进行几何建模和网格划分,获得空间网格单元上的电磁场信息,作为整个电磁场仿真分析的计算基础。然后,通过介质球和涂覆电磁介质导体球的算例对比,分析结合共形技术的时域有限差分法(Finite difference time domain,FDTD)在处理介质物体及涂覆涂层介质物体的有效性,结果表明FDTD方法计算结果与级数解吻合。在此基础上,计算和对比了金属旋翼以及涂覆吸波涂层旋翼的RCS特性,分析了典型方位角入射下全机涂覆前后对RCS特性的影响。研究表明:旋翼表面全涂覆雷达吸波材料(Radar absorbing material,RAM)后对直升机旋翼的RCS抑制效果明显,在全机强散射部位涂覆RAM可以显著地降低RCS特性,涂层的使用在直升机的隐身设计中起到关键的作用。     

嵌套网格粘性流动数值模拟用于风洞洞壁干扰研究

运用嵌套网格技术和Navier Stokes数值模拟对机翼半模和翼身组合体试验时风洞的四壁干扰进行综合模拟、评估和修正。计算格式在空间上采用中心有限体积离散 ,在时间上采用多步Runge Kutta时间步长格式进行积分。计算结果证明了该方法的可行性和优越性。     

高超声速马赫反射

通过实验和数值计算模拟了激波（马赫数７．６）绕４０度压缩拐角的流动。实验是在气动中心激波管中完成的，数值模拟采用了层流的Ｎ－Ｓ方程，对无粘项、粘性项和化学源项分别采用了迎风ＴＶＤ格式、中心差分格式和点隐式的方法进行差分离散，计算了平衡与非平衡条件下的马赫反射流动，其数值计算结果与实验结果符合甚好。     

弹尾超声速轴对称喷流的正格式数值模拟(英文)

采用二阶正格式方法对超声速轴对称喷流流动进行数值模拟。将二维守恒方程的正格式方法发展到轴对称Euler方程组的求解,并对导弹尾部超声速伴随射流进行了数值计算。数值结果与实验照片所反映的流动特征吻合较好,与三阶精度间断有限元方法计算结果吻合,相比于二阶高分辨率TVD格式的计算结果,正格式方法的计算结果在间断点处具有较大的梯度变化。这表明该方法对激波具有较强的捕捉能力,在间断处不会出现数值振荡,对弹尾声音速伴随射流的数值模拟真实、有效。     

基于嵌套网格技术和多种空间离散格式的旋翼流场数值模拟

给出了三维Euler/N-S方程数值模拟悬停状态旋翼流场的方法和模型.在空间离散方法上分别采用Jameson中心差分格式、Van Leer矢通量分裂格式、AUSM格式3种方法,同时加入了Van Albada限制器,并应用了嵌套网格方法.文中采用Hole-Map方法来确定洞边界,采用Inverse-Map方法来搜寻贡献单元.最后给出了旋翼桨叶表面压力分布的计算结果,并与实验数据进行了对比,证明以上3种空间离散格式在旋翼流场计算中的准确性.     

环视SAR成像处理中的几何失真校正算法

环视合成孔径雷达(SAR)图像用于导弹精确来制导中的景象匹配处理,需要满足几何精度的要求.校正由雷达平台不规则运动和天线圆周扫描导致的图像几何失真,是环视SAR成像处理中的关键步骤.文中在利用线性距离多普勒算法生成子图像的前提下,提出了一种基于像源与像点映射关系的距离多普勒域图像几何失真校正算法.该算法无需复杂的坐标系转换计算,实现了360°范围内任意波束指向上SAR子图像几何失真校正.地面点目标仿真和实测教据成像结果证明了该算法的有效性.     

非结构网格高超声速绕流数值模拟及涵道构型减阻特性分析

在非结构网格上对二维高超声速化学非平衡粘性绕流进行了数值模拟，并应用此方法对涵道构型高超声速减阻特性进行了数值分析。空间离散采用VanLeer逆风通量分裂格式，时间推进采用显式的Runge—Kutta格式。化学非平衡动力学模型为五组元五反应模型，对化学反应源项进行了点隐式格式处理，温度场的计算采用牛顿迭代法。分别对二维类弹头体涵道构型、二维半圆柱模型的高超声速绕流流场进行了数值模拟，得到了数值结果，并与完全气体绕流计算结果进行了对比分析。结果表明本文方法可以应用于高超声速绕流计算，并且涵道构型具有高超声速减阻特性。     

二维复值Ginzburg-Landau方程的一个高阶紧致ADI差分格式

对二维复值金兹堡朗道(Ginzburg-Landau,GL)方程提出一个基于时间分裂的高阶紧致交替方向隐式有限差分格式。本文通过时间分裂法将GL方程分裂成一个非线性子问题及两个线性子问题,对非线性子问题以及其中一个线性子问题均通过精确积分进行计算,并对另一线性子问题构造紧致交替方向隐式差分格式进行数值计算。实际计算中,在每一时间步,利用追赶法求解一族常系数三对角线性代数方程组,从而使得算法既具有较高精度又拥有较快的计算速度。数值实验表明该算法在时间和空间方向分别具有二阶和四阶精度,并模拟了方程的一些动力学行为。     

一类新的二阶精度非交错形式差分格式的构造

本文对二维标量非线性双曲型守恒律，以非交错形式Ｌａｘ－Ｆｒｉｅｄｒｉｃｈｓ格式为基本模块，构造了一类新的不须解Ｒｉｅｍａｎｎ问题的差分格式，并证明其具有时空二阶精度。数值实验表明，格式具有较高的分辨激波、稀疏波及接触不连续的能力     

基于OPNET的低轨卫星的点波束确定算法(英文)

低轨道卫星通信网络具有低时延的全球覆盖优势。文中提出了一种新的低轨卫星的波束判断方法。该方法首先需要确定用户是否在中心波束,如果不在中心波束,通过两个步骤确定用户所在的边缘波束:第一步采用二维波束算法;第二步采用三维波束算法。基于这些算法,在OPNET平台上实现了具有动态节点和网络协议的铱星通信仿真系统,并获得仿真结果。结果显示,最长可视时间准则明显增加了卫星服务的时间,并且验证了波束宽度模型的应用。     

医学图像光线投影体绘制中重采样快速算法研究

光线投影算法是体绘制算法中图像效果比较好的方法,但存在运算量大,绘制速度慢的问题,如何提高重采样速度是加速光线投影算法的关键。本利用重采样点在两坐标系中的矩阵变换特性,减少矩阵运算量,同时结合三维数据场在像平面的投影减少光线投影数目;利用包围盒技术避免对空体元的采样,并且通过将Bresenham算法扩展至三维确定每个重采样点所在体元的编号。实验结果表明,本提出的重采样优化算法成像速度比标准光线投射算法快2—3倍,而成像的质量与标准算法基本没有区别。本提出的算法,既能保证绘制质量,又能显减少计算量,提高体绘制的速度。     

三维旋转效应对叶片非定常气动特性的影响

在旋转的作用下，流动分离向下游延迟，结果大大改变了风力机叶片的气动特性，这种由旋转所带来的影响被称为三维旋转效应。本文以叶片表面边界层理论分析的结果为基础，给出了三维旋转效应对气流分离影响的解析关系，使得这一复杂问题可以方便地解决，为准确计算旋转叶片的气动栽荷提供基础。通过算例表明，当这一三维旋转效应应用于叶片非定常气动特性计算时，能大大改善计算结果。     

噪声与共振环境下的数字助听器鲁棒声源定位法

提出了一种噪声与共振环境下的数字助听器鲁棒声源定位方法,该方法基于房间共振模型和多通道自适应特征值分解算法,利用自适应次梯度投影方法迭代地估计声源至各麦克风的脉冲响应系数,从而获得各麦克风对的时延差,并利用几何方法计算声源位置。针对眼镜数字助听器四元麦克风方阵的仿真表明,所提出的自适应次梯度投影方法在强噪声环境中的鲁棒性和定位精度优于传统NLMS算法。     

基于时域有限差分法平面光波的一维吸收边界条件

吸收边界条件的选取在使用时域有限差分法分析设计衍射光学元件的过程中有着举足轻重的作用。文章从麦克斯韦方程组出发,推导出了一维吸收边界条件的解析解及该吸收边界对TEM波的反射系数,用数值模拟验证了在激励源分别处于讨论区域右边、中间,且在TEM波传播的同时一直存在的情况下,一维Mur一阶吸收边界条件和波延迟法的优质吸收效果。对两种边界条件的计算误差进行了对比和分析,结果表明Neum ann边界条件波形稳定早,数值精度高,Mur一阶边界条件波形稳定稍迟,有一定误差,但吸收效果很好。     

MRTD方法的色散特性分析和电磁散射应用

将基于Daubechies尺度函数的时域多分辨(Multiresolution time-domain,MRTD)方法应用于三维目标的电磁散射和雷达目标特性分析中,并对其Courant稳定性条件和色散特性进行了分析.在入射波引入方面,提出应用总场/散射场技术,在连接边界周围定义一些"修正区域",并推导出一系列"修正区域"内的迭代公式,把入射场作为"连接边界条件"引入到计算区域.理论分析和实验结果表明,基于Daubechies尺度函数的MRTD方法和入射波引入方法是有效的,且与传统的FDTD方法相比,MRTD方法在保持计算精度的前提下能够节省计算资源.     

三维复杂腔体的时域有限差分法分析

在计算三维复杂腔体雷达截面的时域有限差分法中，提出了一种简单的计算域设置方法。它由腔壁和腔口面外侧的吸收边界围成，与包围目标的常用计算域设置方法相比，忽略了腔边缘的绕射作用，但具有节省计算机存储空间，无需额外抽取外腔壁对电磁波的散射。数值结果表明该设置方法是合理可行的。     

二维多介质流动问题的界面处理方法

就二维可压缩多介质流动问题的数值模拟,给出了一种新的界面处理方法。通过在界面处构造Riemann问题,利用Riemann问题的解分别定义界面两边流体的边界条件,由于Riemann问题的解准确地描述了界面处流体的流动状态,因此得到了更加准确的界面边界条件。本文将由Riemann问题的解得到的界面速度外推到整个流场,重新定义速度场,避免了由于速度的大梯度变化而导致的Level—Set等值线相互交错,因而得到了更加精确的界面位置。利用该方法对水下激波与柱型气泡相互作用问题进行数值模拟,结果表明该方法能准确地捕捉各种物理现象。