耦合径向基函数与多项式基函数的无网格方法

耦合径向基函数和多项式基函数,形成一种新的近似函数.该近似函数对散乱分布的离散数据点进行逼近时,只需节点信息,不需要划分网格.详细描述了耦合近似函数的建立、属性、插值行为及其形函数和形函数导数的性质.最后引入修正变分原理和单位分解积分技术求解边值问题,并给出了计算实例,表明耦合径向基函数和多项式基函数是一种非常有效的方法.     

一种新的并行代数多重网格粗化算法

近年来,受实际应用领域中大规模科学计算问题的驱动,在大规模并行机上实现代数多重网格（AMG）算法成为数值计算领域的研究热点。本文针对经典AMG方法,提出一种新的并行网格粗化算法一多阶段并行RS算法（MPRS）。我们将新算法集成到了高性能预条件子软件包Hypre中。大量数值实验结果显示,新算法适合更广泛的问题,相对其他并行粗化算法,明显地改善了AMG并行计算的可扩展性。对三维27点格式有限差分离散的Poisson方程,在64个处理机上并行AMG求解,含8百万个未知量,新算法比RS3算法减少了近60的三维Poisson方程,近32万个未知量,在16个处理机上并行AMG—GMRES求解,新算法所需的迭代步数大约为其他粗化算法的一半,显示了很好的算法可扩展性。     

守恒型间断跟踪法中对于间断及其相互作用的处理

将一种基于解的守恒性质的一维间断跟踪法(守恒型间断跟踪法)中所包含的一些技巧,如“网格堆积”等做了进一步发展,以处理被跟踪间断的相互作用.对这些技巧进行了详细的描述和讨论,用数值例子说明其有效性.     

非结构网格热化学非平衡辐射流场数值计算

针对三维热化学非平衡辐射流场设计了基于非结构网格的数值计算方法. 根据原子分子光谱理论逐条计算了100$\sim$1\,500\,nm间N, O, N$^{ + }$, O$^{ + }$的谱线以及N计算流体力学; 辐射; 热化学非平衡; 非结构网格; 有限体积法针对三维热化学非平衡辐射流场设计了基于非结构网格的数值计算方法.根据原子分子光谱理论逐条计算了100～1 500nm间N,O,N+,O+的谱线以及N2,O2,NO,N+2等分子的10个谱带,特别分析了NO的β'带,γ'带,δ带和ε带的辐射特性.采用耦合辐射的双温模型计算热化学非平衡流场,辐射源项通过直接求解辐射输运方程RTE(radiative transport equation)获得.在空间和方向上分别离散后,利用有限体积法求解不同方向上的辐射输运方程.计算得出了再入飞行器前驻点的辐射强度分布.采用该数值方法计算了MUSES-C模型在速度为11.6km/s时的绕流流场及前驻点处的辐射热流密度.并通过对比分析了热辐射对流场的影响.     

离散电阻加载的大型垂直极化EMP辐射波模拟器的并行FDTD模拟

将电阻设置为跨网格有耗介质，进行离散电阻加载的大型垂直极化电磁脉冲（EMP）辐射波模拟器时域场的并行FDTD模拟.给出不同离散电阻加载方式下模拟器辐射场的时域波形，分析离散电阻加载对模拟器辐射场的影响，并给出离散电阻加载的模拟器有圆筒效应物时的时域波形.结果表明：电阻加载能有效减小模拟器顶端的电流反射，从而改善辐射场的时域波形，但电阻加载也可能会导致模拟器辐射性能的降低；当效应物开孔所在的侧面靠近模拟器且与模拟器极化方向平行时耦合进入圆筒内的电磁波能量较多.计算方法具有通用性，可用于进行离散电阻加载的其它大型垂直极化EMP辐射波模拟器时域场的模拟，包括导电地面存在且有效应物存在时的时域耦合场的模拟.     

空中强爆炸冲击波地面反射规律数值模拟研究

为准确预测空中强爆炸产生的冲击波载荷分布，基于欧拉坐标系建立了能够模拟具有高密度比、高压力比的强激波问题的二维多介质流体数值方法。结合网格自适应技术，对1 kt TNT当量的空中强爆炸在不同爆炸高度下的冲击波地面反射过程进行了数值模拟，并考虑了真实气体状态方程和空气随高度不均匀分布的影响。计算得到了地面上距爆心投影点大尺度范围内的反射超压和冲量等冲击波载荷分布，并给出了冲击波载荷随爆高的变化规律。     

四阶分数阶扩散波动方程的两网格混合元快速算法

本文研究四阶分数阶扩散波动方程模型的基于新混合元方法的快速两网格算法.讨论该方法的稳定性,推导三个未知函数的$L^2$模意义下的最优误差估计.最后通过数值例子验证两网格混合元算法的高效性和理论结果的正确性.     

一种新的求解对流占优问题的自适应网格细化方法

在均匀网格上求解对流占优问题时,往往会产生数值震荡现象,因此需要局部加密网格来提高解的精度。针对对流占优问题,设计了一种新的自适应网格细化算法。该方法采用流线迎风SUPG(Petrov-Galerkin)格式求解对流占优问题,定义了网格尺寸并通过后验误差估计子修正来指导自适应网格细化,以泡泡型局部网格生成算法BLMG为网格生成器,通过模拟泡泡在区域中的运动得到了高质量的点集。与其他自适应网格细化方法相比,该方法可在同一框架内实现网格的细化和粗化,同时在所有细化层得到了高质量的网格。数值算例结果表明,该方法在求解对流占优问题时具有更高的数值精度和更好的收敛性。     

惯性约束聚变内爆中基于多块结构网格的高效辐射扩散并行算法

针对三维球形靶丸内爆高效模拟需求和传统笛卡尔正交网格上辐射加源困难的问题, 发展一种多块结构非正交网格生成方法, 并基于此种计算网格提出高效的三维扩散格式并行算法, 将其应用于辐射流体方程组的求解和三维内爆不对称性的数值模拟, 数值结果显示了算法的有效性。并行性能测试显示该算法可扩展到5400个核上, 并行效率达到69%。     

光子晶体多组元缺陷态问题研究

缺陷态光子晶体可以用于制作良好的谐振器、偏振器、滤光器等光学器件,具有重要的应用价值。本文发展了光子晶体缺陷态问题的PG有限元界面问题计算方法,有效地处理了各种不同组元体系、几何结构、界面形状、材料属性以及模态的光子晶体缺陷态问题。数值结果表明,二组元结构单点缺陷对带隙的影响较小,只是使局部范围内的波继续传播而产生一条缺陷带,多点缺陷使一些特定范围内的波可以传播而产生多条缺陷带,线缺陷产生的影响较大,可以使整个禁带消失。结合线缺陷与点缺陷,波导结构中的侧点缺陷可以有效地应用于光子晶体阻带内诱导窄通带或在波导的通带内诱导非常窄的阻带。三组元结构引入了不均匀介质、复杂介质形状以及不同几何结构的缺陷态。通过计算与分析发现Ω₃区域的介质形状对结果影响比较有限,表面层越不光滑禁带越窄,n型缺陷态在TM模中的高频区域更容易产生禁带。对于TE模来说,n型与v型的缺陷态更容易产生禁带。