双曲型守恒律方程的小波解法

本文基于Hamilton-Jacobi方程的小波Galerkin近似和微分算子的小波表示,讨论一维双曲型守恒律方程初值问题的Daubechies小波解．由于小波在空间和时间上的局部性,本方法适用于处理具有奇异解的问题,可以有效的防止数值振荡．数值试验的结果表明,本方法是可行的．     

一种求解双曲守恒律方程的中心型WENO格式

本文发展了一种中心型加权本质无振荡(WENO)格式.该格式通过在原始三阶WENO-JS格式的下风方向增加一个两点候选模板,并将文献[11]中的非线性自适应机制推广到r=2情况,格式记为WENO4-CU.经过近似色散关系分析可以看到,WENO4-CU格式的频谱特性较原始三阶WENO-JS格式具有明显的改进.通过六个典型算例的数值测试表明,WENO4-CU格式在对流动结构的分辨上较原始WENO3-JS、WENO3-M和WENO3-Z格式具有明显提高.     

变形法移动网格求解双曲型守恒律方程研究

在现有格式的基础上要提高偏微分方程数值解的分辨率,自适应移动网格技术是一种有效而且可行的方法。文中将文献[1]提出的自适应移动网格技术推广到三角形网格,并将该方法用于求解双曲型守恒量方程。用网格自适应技术求解守恒律问题时,当生成新网格之后,需要将旧网格上的函数值更新到新的网格,并保持物理量的守恒性。针对这个问题,文中提出了函数值更新过程中守恒型插值公式的具体形式,并针对二维双曲型守恒律方程进行了仿真实验,取得了满意的结果。     

变形法移动网格求解双曲型守恒律方程研究

在现有格式的基础上要提高偏微分方程数值解的分辨率,自适应移动网格技术是一种有效而且可行的方法。文中将文献[1]提出的自适应移动网格技术推广到三角形网格,并将该方法用于求解双曲型守恒量方程。用网格自适应技术求解守恒律问题时,当生成新网格之后,需要将旧网格上的函数值更新到新的网格,并保持物理量的守恒性。针对这个问题,文中提出了函数值更新过程中守恒型插值公式的具体形式,并针对二维双曲型守恒律方程进行了仿真实验,取得了满意的结果。     

二维双曲守恒律标量方程的三阶CWENO-型熵相容算法

应用提出的中心加权基本无振荡(CWENO)-型熵相容格式求解了二维双曲守恒律方程初边值问题,对所得数值结果进行了分析与讨论,并通过与准确解的比较发现该数值求解格式稳定性条件可以取到0.6,而激波过渡带只有1~2个网格单元。实验结果表明该数值求解格式分辨率高且数值稳定性好。     

三维非结构网格上求解双曲型守恒律方程的一类三阶精度有限体积格式

考虑标量双曲型守恒律方程,对三维非结构四面体网格给出了一类满足局部极值原理的三阶精度有限体积格式.方法的主要思想是时间和空间分开处理,时间离散采用三阶TVD Runge-Kutta方法;对空间,在每一个四面体单元上基于最小二乘原理构造一个二次多项式,结合数值解光滑探测器和梯度限制器,使其在光滑区域具有高阶精度,在间断区域满足局部极值原理.该格式具有间断分辨能力高,编程实现简便,计算速度快等优点.典型算例的数值试验表明,该格式是有效的.     

一维非定常对流扩散方程的高阶组合紧致迎风格式

通过将对流项采用四五阶组合迎风紧致格式离散,扩散项采用四阶对称紧致格式离散之后,对得到的半离散格式在时间方向采用四阶龙格库塔方法求解,从而得到了一种求解非定常对流扩散方程问题的高精度组合紧致有限差分格式,其收敛阶为O(h~4+τ~4).经Fourier精度分析和数值验证,证实了格式的良好性能.三个数值算例包括线性常系数问题,矩形波问题和非线性问题,数值结果表明:该格式具有很高的分辨率,且适用于对高雷诺数问题的数值模拟.     

一种基于WENO重构的半离散中心迎风格式

通过三阶WENO重构和半离散中心迎风数值通量的结合,给出了一种求解双曲型守恒律方程的三阶半离散中心迎风格式,格式保持了中心差分格式方法简单的优点．数值计算的结果表明该方法具有较高的分辨率．     

解扩散方程的半离散化方法

对具有周期解的扩散方程ｕｘ＝ａｕｘｘ（ａ＞０），我们用半离散化方法构造了具有任意阶精度自然数）的两层显格式。本文讨论了１≤ｐ≤１０，１≤ｑ≤５的情形，得到了一些格式，其精度和稳定性都优于现有的格式。     

对流扩散方程并行求解方法研究综述

对流扩散方程是一类典型的偏微分方程,其并行求解方法对其他微积分方程的并行求解具有借鉴意义。对对流扩散方程的并行求解方法进行综述,分为显式直接并行、隐式迭代并行、交替分组显式并行和Monte Carlo并行四种并行求解方法,对其中涉及的计算原理进行描述,给出示例,并指出进一步研究方向。     

双曲型守恒律的一种三阶松弛格式

对一维双曲型守恒律,给出了一种形式更简单、计算量更小的三阶松弛格式.该格式以三阶WENO重构和三阶显隐式Runge-Kutta方法为基础.由于不用求解Riemann问题和计算非线性通量函数的雅可比矩阵,所以本文格式保持了松弛格式简单的优点.数值试验表明:该方法具有较高的分辨率.     

A weighted ENO-flux limiter scheme for hyperbolic conservation laws

We present a new scheme that combines essentially non-oscillatory (ENO) reconstructions together with monotone upwind schemes for scalar conservation laws interpolants. We modify a second-order ENO polynomial by choosing an additional point inside the stencil in order to obtain the highest accuracy when combined with the Harten–Osher reconstruction-evolution method limiter. Numerical experiments are done in order to compare a weighted version of the hybrid scheme to weighted essentially non-oscillatory (WENO) schemes with constant Courant–Friedrichs–Lewy number under relaxed step size restrictions. Our results show that the new scheme reduces smearing near shocks and corners, and in some cases it is more accurate near discontinuities compared with higher-order WENO schemes. The hybrid scheme avoids spurious oscillations while using a simple componentwise extension for solving hyperbolic systems. The new scheme is less damped than WENO schemes of comparable accuracy and less oscillatory than higher-order WENO schemes. Further experiments are done on multi-dimensional problems to show that our scheme remains non-oscillatory while giving good resolution of discontinuities.     

A high-order non-oscillatory central scheme with non-staggered grids for hyperbolic conservation laws

In this work, we present a scheme which is based on non-staggered grids. This scheme is a new family of non-staggered central schemes for hyperbolic conservation laws. Motivation of this work is a staggered central scheme recently introduced by A.A.I. Peer et al. [A new fourth-order non-oscillatory central scheme for hyperbolic conservation laws, Appl. Numer. Math. 58 (2008) 674–688]. The most important properties of the technique developed in the current paper are simplicity, high-resolution and avoiding the use of staggered grids and hence is simpler to implement in frameworks which involve complex geometries and boundary conditions. Numerical implementation of the new scheme is carried out on the scalar conservation laws with linear, non-linear flux and systems of hyperbolic conservation laws. The numerical results confirm the expected accuracy and high-resolution properties of the scheme.     

A modified fifth-order WENO scheme for hyperbolic conservation laws

Recently a WENO scheme, with smoothness indicators constructed based on $L_1$ measure is introduced by Ha et al. (2013) and the improved version of this scheme is presented by Kim et al. (2016), referred to as WENO-NS and WENO-P schemes respectively. These schemes perform better than the existing many fifth-order WENO schemes for the problems which contain discontinuities and attain fifth-order accuracy at the critical points where the first derivative vanishes but not at the points where the second derivatives are zero. This paper deals with modification of the above said methods to obtain a new fifth-order weighted essentially non-oscillatory (WENO) scheme. A new global-smoothness indicator is proposed which shows an improved behavior over the solutions of WENO-NS and WENO-P schemes and the proposed scheme attains an optimal order of approximation, even at the critical points where the first and second derivatives vanish but not the third derivative. Examples are taken in the numeric section to check the robustness and accuracy of the proposed scheme for one and two-dimensional system of Euler equations.     

Unconditionally stable ADI scheme of higher-order for linear hyperbolic equations

An unconditionally stable alternating direction implicit (ADI) method of higher-order in space is proposed for solving two- and three-dimensional linear hyperbolic equations. The method is fourth-order in space and second-order in time. The solution procedure consists of a multiple use of one-dimensional matrix solver which produces a computational cost effective solver. Numerical experiments are conducted to compare the new scheme with the existing scheme based on second-order spatial discretization. The effectiveness of the new scheme is exhibited from the numerical results.     

Flux limiters fourth and sixth order high resolution method for hyperbolic conservation laws

This paper presents new fourth and six order accurate explicit finite difference scheme for hyperbolic conservation laws. The technique of obtaining high resolution, total variation non increasing oscillation free and explicit scalar difference scheme is derived by adding multi-limiters of antidiffusive flux to a first order scheme.     

A 15-point high-order compact scheme with multigrid computation for solving 3D convection diffusion equations

We propose a method with sixth-order accuracy to solve the three-dimensional (3D) convection diffusion equation. We first use a 15-point fourth-order compact discretization scheme to obtain fourth-order solutions on both fine and coarse grids using the multigrid method. Then an iterative mesh refinement technique combined with Richardson extrapolation is used to approximate the sixth-order accurate solution on the fine grid. Numerical results are presented for a variety of test cases to demonstrate the efficiency and accuracy of the proposed method, compared with the standard fourth-order compact scheme.     

A new unconditionally stable ADI compact scheme for the two-space-dimensional linear hyperbolic equation

We formulate a new alternating direction implicit compact scheme of O(τ²+h ⁴) for the linear hyperbolic equation u _tt +2α u _t +β² u=u _xx +u _yy +f(x, y, t), 0<x, y<1, 0<t≤T, subject to appropriate initial and Dirichlet boundary conditions, where α>0 and β≥0 are real numbers. In this article, we show the method is unconditionally stable by the Von Neumann method. At last, numerical demonstrations are given to illustrate our result.