基于无网格稳定化方案求解非稳态强对流问题的自适应节点加密技术

无网格Taylor最小二乘(MFLS)稳定化方案可有效地消除无网格Galerkin方法求解对流占优问题时产生的数值伪振荡,但当对流作用很强或纯对流时,它的求解效果不尽人意.因此,本文基于MFLS稳定化方案给出了一种自适应节点加密技术.该技术将无网格方法中背景积分单元作为自适应节点加密时物理量梯度指标的控制单元,并计算该控制单元上的物理量梯度指标;然后将其与给定的物理量梯度指标限进行比较,标识出大梯度区域从而进行自适应节点加密.数值实验表明,当求解对流作用很强的问题或纯对流问题时,这种基于MFLS稳定化方案的自适应节点加密技术不仅能有效地标示出数值振荡区域,而且可以彻底地消除数值伪振荡.     

对流占优问题的无网格自适应布点方案

应用常规数值方法求解对流占优的对流扩散方程时会出现非物理的数值伪振荡现象.因此本文提出了一种基于无网格径向点插值法的自适应布点方案,并成功地解决了对流占优时的数值伪振荡问题.在自适应布点的实施过程中,该方案将无网格方法中的背景积分单元作为自适应控制的梯度计算单元,并将该控制单元场函数梯度的大小作为自适应的梯度控制指标,然后给定相应的梯度控制限,通过控制指标和梯度限的比较来指示高梯度区域进行自适应中心加点和梯度计算单元的分解.数值结果表明:这种基于无网格径向点插值法的自适应布点方案不仅能有效地消除对流占优时的数值伪振荡现象,而且它还具有计算精度高、数值稳定性好、算法实施简单、前后处理方便的优点.     

流动问题无网格Galerkin方法的稳定化方案研究

直接运用无网格Galerkin方法求解对流占优的非线性对流扩散方程及纯对流方程,会出现数值伪振荡现象。本文基于无网格Galerkin方法,构造了MFSUPG(Meshfree Streamline Upwind Petrov-Galerkin),MF-GLS(Meshfree Galerkin Least-Square),MFSGS(Meshfree Sub-Grid Scale)及MFLS(Meshfree Least-Square)四种稳定化方案。数值实验表明:四种稳定化方案中,MFLS的通用性最强。耦合MFLS的无网格Galerkin方法能很好地求解对流占优的非线性对流扩散方程及纯对流方程,具有计算精度高、稳定性好、前后处理方便、算法实施简单的优点,并能捕捉解的大梯度变化。     

线性定常对流占优对流扩散问题的无网格解法

应用无网格Galerkin方法求解对流占优对流扩散问题时会出现非物理现象的数值伪振荡,本文将SUPG方法、GLS方法、SGS方法与无网格Galerkin方法相耦合,成功解决了对流扩散方程中对流项占优时的数值伪振荡问题。运用本文构造的方法,采用线性基和具有C2连续的权函数,应用移动最小二乘法可容易地构造高阶导数连续的形函数,从而避免了有限元方法中当采用线性元插值时,因忽略稳定项中二阶导数项而降低计算精度和稳定性的问题。数值实验表明:本文构造的方法具有计算精度高、稳定性好、计算算法实施简单、前后处理方便的优点,这些方法不仅能适用于对流项占优问题,而且也能很好地消除反应项占优时的数值伪振荡问题。     

一种新的求解对流扩散边值问题的无网格方法

基于配点法和楔形基函数,提出了一种新的求解对流扩散边值问题的无网格方法。通过一维和二维的问题验证了该数值方法的可行性;并根据数值算例和分析,可以看到该数值方法能达到满意的收敛效果。该数值方法的隐格式形式能够有效地消除对流占优问题的数值振荡现象,是一种真正的无网格方法。     

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

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

二维对流扩散方程的欧拉—拉格朗日分裂格式

本文在[1]基础上发展了一种有效的处理大P_e(R_e)数、非定常二维对流扩散方程的欧拉-拉格朗日(E-L)分裂格式,由于方法本质上与区域形状无关,且不需再分网格,因此是一种无网格的E-L方法,特别对于定常流动,E.-L.分裂格式可以导致比一阶迎风格式更精确的单调、无振荡格式,文中对于常系数、变系数和非线性的二维非定常和定常对流扩散方程的(初)边值问题进行了数值计算,数值结果与精确解的比较表明,本方法具有很好的精度,解是单调无振荡的,比通常一阶迎风格式具有较少的数值扩散,最大计算网格P-e(R-e)数可达100—500。     

爆炸冲击波强间断问题的高阶伪弧长算法研究

为了提高对冲击波强间断处的分辨率,通过引入弧长参数,使网格自适应地朝着间断处移动,并结合高精度WENO数值格式,进而达到了对大梯度物理量的高分辨率捕捉。针对网格移动造成的非均匀和非正交现象,通过坐标变换,使得计算过程在均匀正交的计算空间中进行。通过和有限体积下的数值结果对比,结合数值误差分析,可以看到高阶伪弧长数值算法不仅保证了高精度而且对间断的捕捉更加明显,在间断附近解的整体光滑性较好,网格的自适应移动使得解的奇异性得到了削弱,因此可以削弱高阶格式容易引起数值振荡这个缺点。最后采用高阶伪弧长算法计算了化学反应流问题,结果表明高阶伪弧长算法有着较快的收敛率,对于解决爆炸与冲击强间断问题有着较为明显的优势。     

瞬态对流-扩散方程的变分多尺度解法

根据变分多尺度思想,求解了瞬态线性和非线性对流-扩散方程。文中为了简化"细"尺度方程的求解,忽略了该方程的瞬态性,分别用高阶多项式泡函数(High-order Polynomial Bubble)和自由残量泡(Residual Free Bubble)函数近似"细"尺度解,进而引入了消除数值伪振荡的稳定化结构。数值算例验证了本文方法的精确性、稳定性和对高Peclet数问题的适应性,证明了上述对"细"尺度模型的简化是可行的。     

平面裂纹问题的h, p, hp型自适应无网格方法的研究

无网格方法以其独特的优点：不需“网格”（即节点间的连接信息）划分,特别适合自适应的分析,在分析中只需要高梯度域简单地插入离散点（ｈ型）或保持模型节点数、分布、覆盖大小均不变,中增加高误差覆盖上的函数的多项式阶次（ｐ型）,便可以得到更高精度的数值模型。针对平面弹性问题发展和推导一种显式后验误差指示公式,对平面裂纹实例进行了ｈ型,ｐ型,ｈｐ型三种不同类型的无网格自适应分析,数值分析结果表明了这种自适应