二维非线性Schr(o)dinger方程的辛与多辛格式

对满足周期边界条件的二维非线性Schrödinger方程,运用中心差分对该方程进行空间离散, 得到一个有限维Hamilton系统,然后用隐式Euler中点格式进行时间离散得到其辛格式. 针对该方程的多辛形式, 运用有限体积法离散,得到一种直平行六面体上的中点型多辛格式. 用所构造的辛与多辛格式对二维非线性Schrödinger方程的平面波解和奇异解进行数值模拟,结果验证了所构 造格式的有效性. 最后, 根据计算结果,对两种格式进行了分析和比较.       

一种基于过积分的能量稳定通量重构方法

能量稳定通量重构(Energy Stable Flux Reconstruction,ESFR)方法在求解线性对流方程时具有能量稳定性质.但在求解非线性方程时能量稳定性质的实现需要采用L²投影,否则可能由于存在混淆误差,导致不稳定.本文将ESFR与过积分相结合构造具有良好去混淆效果的高阶通量重构(Flux Reconstruction,FR)方法.采用积分点大于求解点(Q> P)的取点方式,从理论上分析了格式的能量稳定特性.从数值上对比了gDG与gSD两种修正函数,三种不同过积分取点方式,并对比过积分与非过积分形式的ESFR(Q=P).通过对一维非均匀线性对流方程、二维等熵涡及欠解析涡流算例的模拟,结果表明:在g_SD修正函数下,ESFR(Q> P)格式比ESFR(Q=P)格式去混淆效果更好,数值误差更小;对比两种修正函数,g_DG修正函数数值误差更小,更稳定:对比三种过积分通量点分布,选定g_DG修正函数时,通量点取Legendre-GaussLobatto(LGL)点或者通量点基于高斯权重...