一维双曲守恒律的龙格-库塔控制体积间断有限元方法

给出数值求解一维双曲守恒律方程的新方法——龙格-库塔控制体积间断有限元方法(RKCVDFEM),其中空间离散基于控制体积有限元方法,时间离散基于二阶TVB Runge-Kutta技术,有限元空间选取为分段线性函数空间.理论分析表明,格式具有总变差有界(TVB)的性质,而且空间和时间离散形式上具有二阶精度.数值算例表明,数值解收敛到熵解并且对光滑解的收敛阶是最优的,优于龙格-库塔间断Galerkin方法(RKDGM)的计算结果.     

球几何中子输运保正线性间断有限元格式

针对球几何中子输运方程线性间断有限元方法计算的负中子通量问题,构造了保正线性间断有限元格式,该格式保持中子角通量0阶矩和1阶矩。现有方法计算中子角通量非负时,采用传统的线性间断有限元方法,求解线性方程组;原方法计算出现负通量,则采用构造的保正格式,求解非线性方程组。编制了球几何中子输运问题保正格式程序模块,并集成到应用程序。数值算例表明构造的保正格式计算的中子通量非负,有效降低数值误差,提高数值计算的精度。     

球几何中子输运保正线性间断有限元格式

针对球几何中子输运方程线性间断有限元方法计算的负中子通量问题,构造了保正线性间断有限元格式,该格式保持中子角通量0阶矩和1阶矩。现有方法计算中子角通量非负时,采用传统的线性间断有限元方法,求解线性方程组;原方法计算出现负通量,则采用构造的保正格式,求解非线性方程组。编制了球几何中子输运问题保正格式程序模块,并集成到应用程序。数值算例表明构造的保正格式计算的中子通量非负,有效降低数值误差,提高数值计算的精度。     

多介质流模拟的Runge-Kutta控制体积间断有限元方法（英文）

构造可用于多介质流数值模拟的Runge-Kutta控制体积(RKCV)间断有限元方法.对于多介质流模拟,使用线性和非线性的Riemann问题解法器计算界面处的数值流通量.该方法是一种高精度的数值方法且可以保证流体的局部守恒.数值结果表明,即使是利用线性Riemann问题解法器的计算格式也可获得较好的数值结果.与Runge-kutta间断Galerkin方法的比较展示了本文构造算法的优势.     

空间差分格式对有限体积法数值散射的影响

数值散射是辐射传递方程近似算法中最常见的离散误差。本文主要讨论空间差分格式对有限体积法数值散射的影响。构造激光平行及倾斜入射的物理模型,验证和比较阶梯格式、中心差分格式及指数格式下温度场的计算精度及数值散射特性。计算结果表明,在激光平行入射与倾斜入射两种情况下,阶梯格式引起的的数值散射比菱形格式及指数格式要多,但其计算精度高于菱形及指数格式。不同激光入射条件下,各种差分格式表现出的数值散射分布有明显的差异。     

局部时间步长间断有限元方法求解三维欧拉方程

使用间断有限元方法求解三维流体力学方程.空间剖分采用非结构四面体网格,为了克服显格式在单元网格尺寸差别较大时计算效率低下的问题,在格式中采用局部时间步长技术(LTS),即控制方程在空间、时间上积分得到一种单步格式,既可以局部计算每个单元又避免了Runge-Kutta高精度格式处理三维问题时存储量过大的问题.为了提高流体力学方程计算精度,在计算单元边界的数值流通量时使用任意高阶精度方法(ADER).数值算例表明格式稳定有效.     

间断有限元方法在弹尾超音速喷流计算中的应用

采用间断有限元方法对超音速无粘喷流流动进行数值模拟.将二维双曲守恒方程的间断有限元方法发展到轴对称Euler方程,并就某导弹尾部超音速伴随射流进行数值计算.计算结果与实验照片反映的流动特征吻合较好,与高精度、高分辨率TVD格式的计算结果相比,间断有限元方法的计算结果在轴线反射点附近具有较高的分辨率,表明该方法对激波具有较强的捕捉能力,在激波阵面上不会产生振荡或抹平间断现象.     

解Hamilton-Jacobi方程的不连续有限元方法

将两类具有不同基函数的有限元应用于Hamilton Jacobi方程,得到了求解Hamilton Jacobi方程的不连续有限元数值格式,并证明了这两类格式数值解在一定条件下收敛于Hamilton Jacobi方程的弱解.数值实例比较了两类格式的精度和分辨间断的能力.     

一种基于AUSM分裂的真正多维HLL格式

文章给出了一种真正多维的HLL Riemann解算器.采用AUSM分裂将通量分解成为对流通量和压力通量, 其中对流通量的计算采用迎风格式, 压力通量的计算采用HLL格式, 且将HLL格式的耗散项中的密度差用压力差代替, 从而使得格式能够分辨接触间断.为了实现数值格式真正多维的特性, 分别计算了网格界面中点和角点上的数值通量, 并且采用Simpson公式加权组合中点和角点上的数值通量得到网格界面的数值通量.为了减少重构角点处状态时的模板宽度, 计算中采用基于SDWLS梯度的线性重构获得2阶空间精度, 而时间离散采用2阶保强稳Runge-Kutta方法.数值实验表明, 相比于传统的一维HLL格式, 文章的真正多维HLL格式具有能够分辨接触间断, 以及更大的时间步长等优点.与其他能够分辨接触间断的格式(例如HLLC格式)不同, 真正多维的HLL格式在计算二维问题时不会出现激波不稳定现象.      

辐射输运菱形差分SN方程的扩散综合加速方法

讨论辐射输运菱形差分S_N方程的源迭代加速方法,在一维情形下给出与输运算子离散相容的线性多频灰体加速计算格式.数值算例表明该加速算法是健壮有效的.     

Analysis of an Implicit Fully Discrete Local Discontinuous Galerkin Method for the Time-Fractional Kdv Equation

In this paper, we consider a fully discrete local discontinuous Galerkin (LDG) finite element method for a time-fractional Korteweg-de Vries (KdV) equation. The method is based on a finite difference scheme in time and local discontinuous Galerkin methods in space. We show that our scheme is unconditionally stable and convergent through analysis. Numerical examples are shown to illustrate the efficiency and accuracy of our scheme.     

一维含化学反应流体力学方程组的ALE间断有限元方法

研究一维含化学反应流体力学方程组的数值模拟方法.结合理想气体状态方程并利用HLLC解法器在各个单元边界处的数值通量,给出ALE间断有限元方法.高阶计算时,使用TVD斜率限制器对数值解可能产生的非物理振荡进行抑制.结果表明：该算法能够保持物理量的守恒性和高精度,并能够清晰地捕捉爆轰波的结构特征.     

二维柱几何非定常中子输运方程基于格式的界面预估校正并行算法

针对二维柱几何非定常中子输运方程的Sn-间断有限元方法,提出基于格式的界面预估校正并行算法.数值算例表明,该并行算法在精度与并行度等诸方面均具有良好的性质,与已有的基于隐式格式的并行扫描算法相比,对于二维中子输运大规模计算问题,并行计算效率较高,并行加速比可增加-倍以上,且可保持原隐式格式的计算精度.     

反应扩散模型在图灵斑图中的应用及数值模拟

反应扩散方程模型常被用于描述生物学中斑图的形成.从反应扩散模型出发,理论推导得到GiererMeinhardt模型的斑图形成机理,解释了非线性常微分方程系统的稳定常数平衡态在加入扩散项后会发生失稳并产生图灵斑图的过程.通过计算该模型,得到图灵斑图产生的参数条件.数值方法中采用一类有效的高精度数值格式,即在空间离散条件下采用Chebyshev谱配置方法,在时间离散条件下采用紧致隐积分因子方法.该方法结合了谱方法和紧致隐积分因子方法的优点,具有精度高、稳定性好、存储量小等优点.数值模拟表明,在其他条件一定的情况下,系统控制参数κ取不同值对于斑图的产生具有重要的影响,数值结果验证了理论结果.     

有限体积元数值方法在大气污染模式中的应用

运用有限体积元方法分析求解大气污染模型问题,分别选取试探函数空间和检验函数空间为一次元函数空间和分片常数函数空间,并且给出L²估计和H¹估计,通过数值实验与有限差分方法进行分析与比较,说明其有效性.为改善大气污染问题的模拟提供实用有效的方法.     

中子扩散方程的混合间断有限元法

构造了求解二维两群中子扩散方程的混合间断有限元格式,研制了相应的计算程序.计算结果表明这个格式是可行的.为了提高计算精度,还需进一步深入系统研究     

局部间断Petrov-Galerkin方法在大气污染模型中的应用

构造数值模拟两类大气污染模型的局部间断Petrov-Galerkin方法.首先通过变量代换将大气污染模型方程转化为与之等价的一阶微分方程组,再利用间断Petrov-Galerkin方法求解微分方程组.该方法既可以选取不同的检验函数和试探函数空间,又可以保持间断Petrov-Galerkin方法的优势.同局部间断有限元方...