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

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

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

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

基于80C51FXX单片机的双冗余控制系统的实现

随着功能强大的单片机在控制领域应用的不断深入,容错控制系统也在不断地发展,在一些特定的场合下,如在航空航天、军事、铁路、石油、化工、电力等重要部门和在恶劣工作环境下工作的计算机控制系统,对系统安全性、可靠性、可用性的要求更高.双冗余系统的两个模块同时执行一样的操作,在其中一个模块出现故障的时候,可以自动判断切换,以保证系统稳定、可靠、不间断的工作.双冗余不仅有较高的可靠性,而且有很高的安全性,因此在控制领域中可广泛应用.     

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

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

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

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

三维可压缩 Navier-Stokes 方程的间断Galerkin 有限元方法研究

拓展了二维间断 Galerkin(DG)有限元方法研究，将该数值方法用于三维可压缩欧拉方程和 Navier-Stokes方程的求解。基于六面体网格单元，采用插值方法将物面的四边形面网格单元构造为弯曲面网格单元，更好地表述了真实物面特征；物面边界相邻体网格单元相应构造为高阶体网格单元，其余体网格单元采用八节点六面体单元，以较小的计算代价使网格满足 DG 方法计算需求。通过对三维带 bump 管道内流、圆球绕流以及旋转流线体绕流进行的数值求解，验证了边界弯曲方法的可行性及 DG 方法的高精度特性。此外，由于采用了隐式计算方法，仅需较少的时间步就能迭代收敛。     

基于非结构/混合网格的高阶精度DG/FV混合方法研究进展

DG/FV混合方法因其具有紧致、易于推广获得高阶格式及相比同阶精度DG方法计算量、存储量小等优点,自提出以来已成功应用于一维、二维标量方程和Euler/N-S方程的求解。综述了DG/FV混合方法的研究进展,重点介绍了DG/FV混合方法的空间重构算法、针对RANS方程的求解方法、隐式时间离散格式、数值色散耗散及稳定性分析、计算量理论分析,并给出了系列粘性流算例的计算结果,包括用于验证混合方法数值精度的库埃特流,以及方腔流、亚声速剪切层、低速平板湍流、NACA0012翼型湍流绕流等。数值计算结果表明DG/FV混合方法达到了设计的精度阶,且相比同阶DG方法计算量减少约40%,而隐式方法能大幅提高定常流的收敛历程,较显式Runge-Kutta的收敛速度提高1~2个量级。     

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

针对多介质可压缩流体动力学问题,提出了一种单元中心型二维MMALE(Multi-Material Arbitrary Lagrangian-Eulerian)方法。在拉氏步,流体力学方程组采用中心型间断有限元方法求解。对于混合网格,采用Tipton压力松弛模型更新物理量,用等参坐标法更新物质中心点坐标。界面重构采用一种健壮的MOF(Moment of Fluid)方法。在重映步提出了基于多边形相交的二阶积分守恒重映方法。该方法分为四个部分:多项式重构、多边形相交、积分和后验校正。多边形相交使用"剪裁投影"算法,显著降低了多边形相交算法的复杂度。后验校正是基于MOOD (Multi-dimensional Optimal Order Detection)限制策略,并做了一些改动以适应多介质的计算。数值算例表明,该方法具有二阶的精度和较好的鲁棒性。     

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

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

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

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