基于CFD高精度算法的雷电电磁场研究

雷电空间电磁场一直是电磁领域研究的热门，成熟的算法例如有限时域差分法（FDTD），传输线矩阵法（TLM）以及频域的矩量法（MoM）在计算雷电问题方面都有广泛的应用。由于计算流体力学（CFD）中的Euler方程与电磁学中的Maxwell方程有着相同的守恒形式，而且采用间断伽辽金方法（DG）已经在流场问题上得到广泛的尝试，因此引入了基于计算流体力学的DG方法来离散时域Maxwell方程，并采用网格分区并行技术加速计算，使用基于DG的圆球雷达散射截面积（RCS）算例进行测试，数值结果一致表明DG算法在求解电磁场问题上的可行性，之后通过计算一段近场雷电通道的电场分布并与解析解、某算法仿真解对比，数据基本吻合，说明该方法适合于雷电电磁场的计算。     

冕环中扭转Alfvén波和快波的耦合共振

推广完备正交函数系(OFSE)展开方法,应用它研究冕环中无耗散扭转Alfven波和快波耦合的时变过程.采用这种方法,从数学上描述了总体波模的固有角频率,从理论和数值计算结果两方面分析了Alfven波和快波耦合时冕环中磁力线Alfven波固有角频率.采用这种方法,分析了耦合驱动项和Alfven波共振的关系.计算结果给出了耦合共振出现位置Alfven波和快波的波幅特征,发现如果足点驱动角频率不等于冕环总体波模的固有角频率,在耦合共振位置当扭转Alfven波幅出现δ断面时,快波振幅沿径向梯度很大,有时近似为一个间断,非常有利于快波耗散;如果足点驱动角频率等于总有体波模的固有角频率,Alfvn波振幅出现δ断面时在径向出现丰富的小尺度结构,共振位置近似为一个间断,非常有利于Alfven波耗散.     

强磁场扰动对宇宙线调制的统计研究

对1978─1982太阳活动高年时发生的激波、强磁场扰动及激波与强磁场扰动共存这三类事件引起的宇宙线变化进行了统计研究,得到如下结果:(1)激波与强磁场扰动共存时引起的宇宙线强度下降最为显着;只有激波或强磁场扰动时,宇宙线的强度变化相对较小;(2)标志速度间断的激波是产生宇宙线Forbush下降的重要因素;(3)速度间断在强磁场扰动对宇宙线的调制中可能起一个触发的作用。      

基于MLPG混合配点法的接触问题形状优化

分别介绍了接触问题和MLPG(Meshless Local Petrov-Galerkin)混合配点法的理论基础,推导了相关公式,给出了两种典型接触状态的定解条件,使用二维线性基函数,采用三次样条曲线权函数,通过移动最小二乘法插值,将MLPG混合配点法运用到接触分析中,使用罚函数法添加本质边界条件,对二维弹性接触问题的接触过程进行模拟,反复迭代得到真实的接触情况,建立了一种新的应力-位移非线性数学求解模型,结合遗传算法对实际工程接触问题进行了求解优化,给出优化结果和目标函数变化曲线,并与相关文献结果比较,验证了该方法的有效性.     

一种中心型间断有限元MMALE方法

针对多介质可压缩流体动力学问题,提出了一种单元中心型二维MMALE(Multi-Material Arbitrary Lagrangian-Eulerian)方法.在拉氏步,流体力学方程组采用中心型间断有限元方法求解.对于混合网格,采用Tipton压力松弛模型更新物理量,用等参坐标法更新物质中心点坐标.界面重构采用一种...     

间断制导下的反辐射无人机被动定位算法研究

目前,有关被动定位算法的研究多是在测量信息连续的情况下展开的,对于测量信息间断情况下的滤波鲜有涉及.针对该问题,以反辐射无人机为研究对象,提出一种基于卡尔曼滤波的间断制导下的被动定位算法,并通过仿真,验证了该方法的有效性.结果表明:提出的新方法能够用于间断制导下的被动定位,且具有一定的定位精度.     

高阶间断有限元法的并行计算研究

根据间断有限元法的数据结构特点,基于METIS网格分区技术,设计并行计算策略,在非结构网格上实现了并行高阶间断有限元法。控制方程的数值通量项使用Local Lax-Friedrichs(LLF)格式计算。设计了并行的牛顿-块高斯赛德尔法(Newton-Block GS)来加速收敛,提高迭代效率。并行性能分析表明,所设计的并行算法能够得到较好的加速比和并行效率,有效地节省计算时间,合理分配内存。这使得采用高阶间断有限元法计算更为复杂的问题成为可能。     

弯曲网格上的间断有限元湍流数值解法研究

采用间断有限元方法对雷诺平均Navier-Stokes(RANS)方程进行了数值求解,对Spalart-Allmaras单方程湍流模型进行了部分修正,使得求解器更加鲁棒。构造了分段高次多项式来逼近真实物面,同时物面附近采用多层弯曲网格来避免网格交叉,此外提出了一种快速计算积分点的曲面物面距的方法。采用混合网格对NACA0012翼型以及RAE翼型进行了数值计算,并与实验数据以及前人数据进行了对比。计算结果表明,采用物面弯曲网格结合修正的湍流模型方法在相对稀疏的网格上就能得到比较好的数值解。     

曲线物面边界条件在DGM中的应用研究

基于二维非结构网格,采用高精度间断Galerkin有限元方法数值求解定常Euler方程。为了有效克服传统的反射物面边界条件在物面处易于生成伪熵降低计算精度的缺陷,提出一种曲线物面边界条件,从而实现对物理模型而不是数值模型进行数值模拟。对经典圆柱亚音速无粘绕流进行数值模拟,结果表明:采用曲线物面边界条件之后,流场解的精度得到很好的提高;此外,在非常稀疏的网格上,通过提高基函数的阶次仍可以得到高精度的数值解。     

Maxwell方程的高阶间断有限元数值解法

采用高精度方法求解时域Maxwell方程,方程的空间离散采用基于计算流体力学(Computational fluid dynamics,CFD)领域的高阶间断有限元格式,非定常时间迭代采用四步龙格-库塔格式。为了提高计算效率,本文采用了Quadrature-free implementation和网格分区并行技术。数值结果表明,采用高阶格式的情况下,采用稀疏网格便可以得到高精度数值解。另外由于本文的方法基于非结构网格,因此非常适合计算复杂外形的情况。