Wavelet-Taylor Galerkin Method for the Burgers Equation

In this paper, we propose a wavelet-Taylor Galerkin method for the numerical solution of the Burgers equation. In deriving the computational scheme, Taylor-generalized Euler time discretization is performed prior to wavelet-based Galerkin spatial approximation. The linear system of equations obtained in the process are solved by approximate-factorization-based simple explicit schemes, and the resulting solution is compared with that from regular methods. To deal with transient advection-diffusion situations that evolve toward a convective steady state, a splitting-up strategy is known to be very effective. So the Burgers equation is also solved by a splitting-up method using a wavelet-Taylor Galerkin approach. Here, the advection and diffusion terms in the Burgers equation are separated, and the solution is computed in two phases by appropriate wavelet-Taylor Galerkin schemes. Asymptotic stability of all the proposed schemes is verified, and the L_∞ errors relative to the analytical solution together with the numerical solution are reported. AMS subject classification (2000) 65M70     

A Novel Staggered Semi-implicit Space-Time Discontinuous Galerkin Method for the Incompressible Navier-Stokes Equations

A new high-order accurate staggered semi-implicit space-time discontinuous Galerkin(DG)method is presented for the simulation of viscous incompressible flows on...     

A Hamiltonian Regularization of the Burgers Equation

We consider a quasilinear equation that consists of the inviscid Burgers equation plus O(α²) nonlinear terms. As we show, these extra terms regularize the Burgers equation in the following sense: for smooth initial data, the α > 0 equation has classical solutions globally in time. Furthermore, in the zero-α limit, solutions of the regularized equation converge strongly to weak solutions of the Burgers equation. We present numerical evidence that the zero-α limit satisfies the Oleinik entropy inequality. For all α ≥ 0, the regularized equation possesses a nonlocal Poisson structure. We prove the Jacobi identity for this generalized Hamiltonian structure.     

A New Discontinuous Galerkin Method for Parabolic Equations with Discontinuous Coefficient

In this paper, a new discontinuous Galerkin method is developed for the parabolic equation with jump coefficients satisfying the continuous flow condition. Theoretical analysis shows that this method is $L^2$ stable. When the finite element space consists of interpolative polynomials of degrees $k$, the convergent rate of the semi-discrete discontinuous Galerkin scheme has an order of$\mathcal{O}(h^k)$. Numerical examples for both 1-dimensional and 2-dimensional problems demonstrate the validity of the new method.     

A Discontinuous Galerkin Method for Fractional Differential Equations

We investigate the numerical treatment of fractional derivatives using a discontinuous Galerkin approach. We consider fractional differential equations in one dimension. The most common definitions, the Riemann-Liouville and Caputo derivatives, as well as a third definition are considered. Different types of initial condition statements are investigated for the extended discontinuous Galerkin method. We obtain numerical results on the order of convergence in the L₂ and H¹ norms. We compare the classical discontinuous Galerkin method with our adapted approach and show results for the different types of fractional derivatives. Further work was carried out on parabolic differential equations with fractional time derivatives, which is to be published cf. [2]. A keystone in the further research is the development of analytical tools for theoretical convergence results. (© 2016 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

计算有弥散和吸附的径向渗流问题的局部化间断Galerkin方法

在Corkburn和Shu新近发展的求解对流扩散方程的局部化间断Galerkin方法的基础上,针对有弥散和吸附的径向渗流问题中出现的推广的对流扩散方程的形式,构造了一种计算有弥散和吸附的径向渗流问题的局部化间断Galerkin有限元方法,为径向渗流问题的求解提供了一个高阶的新方法．对对流-弥散和对流-弥散-吸附两种情况进行了数值实验, 所得结果的相应部分与已知的一些精确解结果和数值结果是一致的,表明方法是可靠的．从计算速度上看,方法也是可行的．     

A posteriori estimates for the Bubble Stabilized Discontinuous Galerkin Method

In this paper, two reliable and efficient a posteriori error estimators for the Bubble Stabilized Discontinuous Galerkin (BSDG) method for diffusion-reaction problems in two and three dimensions are derived. The theory is followed by some numerical illustrations.     

关于K-S方程的非线性伽辽金方法(英文)

该文讨论了关于 K- S方程的伽辽金方法和非线性伽辽金方法的收敛性和 L2 误差估计 ,并得出误差阶一致的结论     

Positivity-Preserving Local Discontinuous Galerkin Method for Pattern Formation Dynamical Model in Polymerizing Actin Flocks

In this paper, we apply local discontinuous Galerkin (LDG) methods for pattern formation dynamical model in polymerizing actin flocks. There are two main difficulties in designing effective numerical solvers. First of all, the density function is non-negative, and zero is an unstable equilibrium solution. Therefore, negative density values may yield blow-up solutions. To obtain positive numerical approximations, we apply the positivity-preserving (PP) techniques. Secondly, the model may contain stiff source. The most commonly used time integration for the PP technique is the strong-stability-preserving Runge-Kutta method. However, for problems with stiff source, such time discretizations may require strictly limited time step sizes, leading to large computational cost. Moreover, the stiff source any trigger spurious filament polarization, leading to wrong numerical approximations on coarse meshes. In this paper, we combine the PP LDG methods with the semi-implicit Runge-Kutta methods. Numerical experiments demonstrate that the proposed method can yield accurate numerical approximations with relatively large time steps.     

Analysis of Sharp Superconvergence of Local Discontinuous Galerkin Method for One-Dimensional Linear Parabolic Equations

In this paper, we study the superconvergence of the error for the local discontinuous Galerkin (LDG) finite element method for one-dimensional linear parabolic equations when the alternating flux is used. We prove that if we apply piecewise $k$-th degree polynomials, the error between the LDG solution and the exact solution is ($k$+2)-th order superconvergent at the Radau points with suitable initial discretization. Moreover, we also prove the LDG solution is ($k$+2)-th order superconvergent for the error to a particular projection of the exact solution. Even though we only consider periodic boundary condition, this boundary condition is not essential, since we do not use Fourier analysis. Our analysis is valid for arbitrary regular meshes and for $P^k$ polynomials with arbitrary $k$ ≥ 1. We perform numerical experiments to demonstrate that the superconvergence rates proved in this paper are sharp.     

广义Burgers方程的Legendre-Galerkin Chebyshev-配置方法

本文对广义Burgers方程的Neumann和Robin型边值问题构造了LegendreGalerkinChebyshev-配置方法.Legendre-GalerkinChebyshev-配置方法整体上按LegendreGalerkin方法形成,但对非线性项采用在Chebyshev-Gauss-Lobatto点上的配置法处理.文中给出了方法的稳定性和收敛性分析,获得了按H1-模的最佳误差估计.数值实验证实了方法的有效性.     

Kuramoto-Sivashinsky方程的B样条Galerkin方法

采用显隐结合的方法对微分算子进行时间离散 ,提出了解 Kuramoto-Sivashinsky方程的全离散 B样条 Galerkin方法 ,由此得到了有限元解的最优阶收敛性及稳定性估计 .最后的数值算例以图形的形式体现了此算法的精确度     

An $h$-Adaptive Runge-Kutta Discontinuous Galerkin Method for Hamilton-Jacobi Equations

In [35, 36], we presented an $h$-adaptive Runge-Kutta discontinuous Galerkin method using troubled-cell indicators for solving hyperbolic conservation laws. A tree data structure (binary tree in one dimension and quadtree in two dimensions) is used to aid storage and neighbor finding. Mesh adaptation is achieved by refining the troubled cells and coarsening the untroubled "children". Extensive numerical tests indicate that the proposed $h$-adaptive method is capable of saving the computational cost and enhancing the resolution near the discontinuities. In this paper, we apply this $h$-adaptive method to solve Hamilton-Jacobi equations, with an objective of enhancing the resolution near the discontinuities of the solution derivatives. One- and two-dimensional numerical examples are shown to illustrate the capability of the method.     

Accuracy Enhancement of Discontinuous Galerkin Method for Hyperbolic Systems

We study the enhancement of accuracy, by means of the convolution post-processing technique, for discontinuous Galerkin(DG) approximations to hyperbolic problems. Previous investigations have focused on the superconvergence obtained by this technique for elliptic, time-dependent hyperbolic and convection-diffusion problems. In this paper, we demonstrate that it is possible to extend this post-processing technique to the hyperbolic problems written as the Friedrichs' systems by using an upwind-like DG method. We prove that the $L_2$-error of the DG solution is of order $k+1/2$, and further the post-processed DG solution is of order $2k+1$ if $Q_k$-polynomials are used. The key element of our analysis is to derive the $(2k+1)$-order negative norm error estimate. Numerical experiments are provided to illustrate the theoretical analysis.     

Local Discontinuous Galerkin Methods for Three Classes of Nonlinear Wave Equations

In this paper, we further develop the local discontinuous Galerkin method to solve three classes of nonlinear wave equations formulated by the general KdV-Burgers type equations, the general fifth-order KdV type equations and the fully nonlinear K(n, n, n) equations, and prove their stability for these general classes of nonlinear equations. The schemes we present extend the previous work of Yan and Shu [30, 31] and of Levy, Shu and Yan [24] on local discontinuous Galerkin method solving partial differential equations with higher spatial derivatives. Numerical examples for nonlinear problems are shown to illustrate the accuracy and capability of the methods. The numerical experiments include stationary solitons, soliton interactions and oscillatory solitary wave solutions.The numerical experiments also include the compacton solutions of a generalized fifthorder KdV equation in which the highest order derivative term is nonlinear and the fully nonlinear K (n, n, n) equations.     

解第一类算子方程的一种新的正则化方法

对算子与右端都为近似给定的第一类算子方程提出一种新的正则化方法，依据广义Ａｒｃａｎｇｅｌｉ方法选取正则参数，建立了正则解的收敛性。这种新的正则化方法与通常的Ｔｉｋｈｏｎｏｖ正则化方法相比较，提高了正则解的渐近阶估计。     

同伦分析法求解Burgers方程的初边值问题

采用同伦分析法求解了Burgers方程的一初边值问题,得到了它的近似解析解.在不同粘性系数情形下,对近似解与精确解进行了比较,发现在粘性系数不是非常小的情况下,用此方法得到的解析解与精确解符合地很好.     

A Modified Weak Galerkin Finite Element Method for Sobolev Equation

For Sobolev equation, we present a new numerical scheme based on a modified weak Galerkin finite element method, in which differential operators are approximated by weak forms through the usual integration by parts. In particular, the numerical method allows the use of discontinuous finite element functions and arbitrary shape of element. Optimal order error estimates in discrete $H^1$ and $L^2$ norms are established for the corresponding modified weak Galerkin finite element solutions. Finally, some numerical results are given to verify theoretical results.     

解第一类算子方程的一种迭代正则化方法

提出了一种求解第一类算子方程的新的迭代正则化方法,并依据广义Arcangeli方法选取正则参数,建立了正则解的收敛性.与通常的Tikhonov正则化方法相比较,提高了正则解的渐近阶估计.     

奇异摄动问题最优阶一致收敛的间断有限元分析

采用非对称内罚间断有限元方法(以下简称NIPG方法)求解一维对流扩散型奇异摄动问题.理论上证明了采用拉格朗日线性元的NIPG方法在Bakhvalov-Shishkin网格上具有最优阶的一致收敛性,即在能量范数度量下其误差估计为O(N~(-1)),其中N为网格剖分中单元个数.数值算例验证了理论分析的正确性.