Hierarchical reconstruction for discontinuous Galerkin methods on unstructured grids with a WENO-type linear reconstruction and partial neighboring cells

The hierarchical reconstruction (HR) [Y.-J. Liu, C.-W. Shu, E. Tadmor, M.-P. Zhang, Central discontinuous Galerkin methods on overlapping cells with a non-oscillatory hierarchical reconstruction, SIAM J. Numer. Anal. 45 (2007) 2442-2467] is applied to the piecewise quadratic discontinuous Galerkin method on two-dimensional unstructured triangular grids. A variety of limiter functions have been explored in the construction of piecewise linear polynomials in every hierarchical reconstruction stage. We show that on triangular grids, the use of center biased limiter functions is essential in order to recover the desired order of accuracy. Several new techniques have been developed in the paper: (a) we develop a WENO-type linear reconstruction in each hierarchical level, which solves the accuracy degeneracy problem of previous limiter functions and is essentially independent of the local mesh structure; (b) we find that HR using partial neighboring cells significantly reduces over/under-shoots, and further improves the resolution of the numerical solutions. The method is compact and therefore easy to implement. Numerical computations for scalar and systems of nonlinear hyperbolic equations are performed. We demonstrate that the procedure can generate essentially non-oscillatory solutions while keeping the resolution and desired order of accuracy for smooth solutions.     

Hierarchical reconstruction for spectral volume method on unstructured grids

The hierarchical reconstruction (HR) [Y.-J. Liu, C.-W. Shu, E. Tadmor, M.-P. Zhang, Central discontinuous Galerkin methods on overlapping cells with a non-oscillatory hierarchical reconstruction, SIAM J. Numer. Anal. 45 (2007) 2442–2467; Z.-L. Xu, Y.-J. Liu, C.-W. Shu, Hierarchical reconstruction for discontinuous Galerkin methods on unstructured grids with a WENO type linear reconstruction and partial neighboring cells, J. Comput. Phys. 228 (2009) 2194–2212] is applied to a piecewise quadratic spectral volume method on two-dimensional unstructured grids as a limiting procedure to prevent spurious oscillations in numerical solutions. The key features of this HR are that the reconstruction on each control volume only uses adjacent control volumes, which forms a compact stencil set, and there is no truncation of higher degree terms of the polynomial. We explore a WENO-type linear reconstruction on each hierarchical level for the reconstruction of high degree polynomials. Numerical computations for scalar and system of nonlinear hyperbolic equations are performed. We demonstrate that the hierarchical reconstruction can generate essentially non-oscillatory solutions while keeping the resolution and desired order of accuracy for smooth solutions.     

自适应间断有限元方法求解双曲守恒律方程

研究自适应Runge-Kutta间断Galerkin (RKDG)方法求解双曲守恒律方程组,并提出两种生成相容三角形网格的自适应算法.第一种算法适用于规则网格,实现简单、计算速度快.第二种算法基于非结构网格,设计一类基于间断界面的自适应网格加密策略,方法灵活高效.两种方法都具有令人满意的计算效果,而且降低了RKDG的计算量.     

Monoslope and multislope MUSCL methods for unstructured meshes

We present new MUSCL techniques associated with cell-centered finite volume method on triangular meshes. The first reconstruction consists in calculating one vectorial slope per control volume by a minimization procedure with respect to a prescribed stability condition. The second technique we propose is based on the computation of three scalar slopes per triangle (one per edge) still respecting some stability condition. The resulting algorithm provides a very simple scheme which is extensible to higher dimensional problems. Numerical approximations have been performed to obtain the convergence order for the advection scalar problem whereas we treat a nonlinear vectorial example, namely the Euler system, to show the capacity of the new MUSCL technique to deal with more complex situations.     

变系数(2+1)维Nizhnik-Novikov-Vesselov方程的三孤子新解

本文为获得非线性发展方程的相互作用解,研究了辅助方程法,并扩展应用辅助方程法和(G'/G)展开法, 获得了变系数非线性(2+1)维Nizhnik-Novikov-Vesselov方程的由椭圆函数、双曲函数、 三角函数和有理函数混合构成的新相互作用解. 关键词： G'/G)展开法')" href="#">(G'/G)展开法 辅助方程法 三孤子解     

构造非线性发展方程无穷序列复合型精确解的一种方法

为了获得非线性发展方程新的无穷序列复合型精确解,给出了Riccati方程的Bäcklund变换和解的非线性叠加公式,符号计算系统Mathematica的帮助下,以广义Boussinesq方程为应用实例,获得了无穷序列复合型精确解.这里包括双曲函数、三角函数与有理函数复合解、双曲函数与三角函数复合解等几种新的无穷序列复合型精确解.该方法在构造非线性发展方程无穷序列复合型精确解方面具有普遍意义. 关键词： 非线性发展方程 非线性叠加公式 Riccati方程 无穷序列精确解     

sine-Gordon型方程的无穷序列新精确解

为了获得sine-Gordon型方程的无穷序列精确解,给出三角函数型辅助方程和双曲函数型辅助方程及其Bäcklund变换和解的非线性叠加公式,借助符号计算系统Mathematica,构造了sine-Gordon方程、mKdV-sine-Gordon方程、(n+1)维双sine-Gordon方程和sinh-Gordon方程的无穷序列新精确解.其中包括无穷序列三角函数解、无穷序列双曲函数解、无穷序列Jacobi椭圆函数解和无穷序列复合型解. 关键词： sine-Gordon型方程 解的非线性叠加公式 辅助方程 无穷序列精确解     

Degasperis-Procesi 方程的无穷序列尖峰孤立波解

本文为了构造非线性发展方程的无穷序列尖峰精确解,给出了Riccati方程的Bäcklund 变换和解的非线性叠加公式,并借助符号计算系统Mathematica,用Degasperis-Procesi方程为应用实例,构造了无穷序列尖峰孤立波解和无穷序列尖峰周期解. 关键词： Riccati方程 解的非线性叠加公式 尖峰孤立波解 Degasperis-Procesi 方程     

Adaptive mesh refinement based on high order finite difference WENO scheme for multi-scale simulations

In this paper, we propose a finite difference AMR-WENO method for hyperbolic conservation laws. The proposed method combines the adaptive mesh refinement (AMR) framework  and  with the high order finite difference weighted essentially non-oscillatory (WENO) method in space and the total variation diminishing (TVD) Runge–Kutta (RK) method in time (WENO-RK)  and  by a high order coupling. Our goal is to realize mesh adaptivity in the AMR framework, while maintaining very high (higher than second) order accuracy of the WENO-RK method in the finite difference setting. The high order coupling of AMR and WENO-RK is accomplished by high order prolongation in both space (WENO interpolation) and time (Hermite interpolation) from coarse to fine grid solutions, and at ghost points. The resulting AMR-WENO method is accurate, robust and efficient, due to the mesh adaptivity and very high order spatial and temporal accuracy. We have experimented with both the third and the fifth order AMR-WENO schemes. We demonstrate the accuracy of the proposed scheme using smooth test problems, and their quality and efficiency using several 1D and 2D nonlinear hyperbolic problems with very challenging initial conditions. The AMR solutions are observed to perform as well as, and in some cases even better than, the corresponding uniform fine grid solutions. We conclude that there is significant improvement of the fifth order AMR-WENO over the third order one, not only in accuracy for smooth problems, but also in its ability in resolving complicated solution structures, due to the very low numerical diffusion of high order schemes. In our work, we found that it is difficult to design a robust AMR-WENO scheme that is both conservative and high order (higher than second order), due to the mass inconsistency of coarse and fine grid solutions at the initial stage in a finite difference scheme. Resolving these issues as well as conducting comprehensive evaluation of computational efficiency constitute our future work.     

热传导方程的一类无网格方法

构造求解热传导方程的一类无网格方法,只要选择好每个节点的适当的邻点集合,便可利用节点信息顺利进行计算.作为特殊情形,也可在各种结构或非结构网格上进行计算.在矩形或均匀平行四边形网格上进行计算时具有二阶精度,当在任意的不规则四边形或三角形网格上计算时仍然是守恒的和相容的,且至少具有一阶精度.作为数值试验,将该方法用于在不规则四边形网格上及四边形与三角形混合网格上求解二维非线性抛物型方程,并在不规则四边形网格上求解二维三温辐射热传导方程,均获得了较为理想的数值结果.     

基于移动网格的熵稳定格式

提出一种基于移动网格的熵稳定格式求解双曲型守恒律方程.该方法利用等分布原理得到新的网格分布,基于守恒型插值公式计算新的网格上的物理量,使用熵稳定数值通量和三阶强稳定Runge-Kutta时间推进方法得到下一时刻的数值解.数值算例表明该格式不仅能有效提高解在间断处的分辨率,而且能消除可能产生的伪振荡.     

二维非线性抛物方程广义差分法/有限体积法的自适应计算

给出三角网上二维非线性抛物方程广义差分法(有限体积法)的一种基于残量估计的后验误差估计,并在此基础上设计了自适应计算方案,以适应物理解在时空的大梯度变化.提出了适合发展方程自适应计算的三角网数据结构(不是树状结构)和灵活的局部粗化算法.     

New Exact Solutions for Konopelchenko-Dubrovsky Equation Using an Extended Riccati Equation Rational Expansion Method

Taking the Konopelchenko-Dubrovsky system as a simple example, some families of rational formal hyperbolic function solutions, rational formal triangular periodic solutions, and rational solutions are constructed by using the extended Riccati equation rational expansion method presented by us. The method can also be applied to solve more nonlinear partial differential equation or equations.     

The extended mapping method and its applications for nonlinear evolution equations

The extended mapping method with a computerized symbolic computation is used to drive some new exact solutions of four nonlinear evolution equations in mathematical physics. As a result, many exact travelling wave solutions are obtained which include new solitary wave solutions, triangular and hyperbolic functions. Solutions in the limiting cases have also been studied. The method is straightforward and concise, and it can also be applied to other nonlinear evolution equations in mathematical physics.     

A Generalized Algebraic Method for Constructing a Series of Explicit Exact Solutions of a (1+1)-Dimensional Dispersive Long Wave Equation

Making use of a new and more general ansatz, we present the generalized algebraic method to uniformlyconstruct a series of new and general travelling wave solution for nonlinear partial differential equations. As an applicationof the method, we choose a (1 1)-dimensional dispersive long wave equation to illustrate the method. As a result, wecan successfully obtain the solutions found by the method proposed by Fan [E. Fan, Comput. Phys. Commun. 153 (2003)17] and find other new and more general solutions at the same time, which include polynomial solutions, exponentialsolutions, rational solutions, triangular periodic wave solutions, hyperbolic and soliton solutions, Jacobi and Weierstrassdoubly periodic wave solutions.     

A nodal discontinuous Galerkin finite element method for nonlinear elastic wave propagation

A nodal discontinuous Galerkin finite element method (DG-FEM) to solve the linear and nonlinear elastic wave equation in heterogeneous media with arbitrary high order accuracy in space on unstructured triangular or quadrilateral meshes is presented. This DG-FEM method combines the geometrical flexibility of the finite element method, and the high parallelization potentiality and strongly nonlinear wave phenomena simulation capability of the finite volume method, required for nonlinear elastodynamics simulations. In order to facilitate the implementation based on a numerical scheme developed for electromagnetic applications, the equations of nonlinear elastodynamics have been written in a conservative form. The adopted formalism allows the introduction of different kinds of elastic nonlinearities, such as the classical quadratic and cubic nonlinearities, or the quadratic hysteretic nonlinearities. Absorbing layers perfectly matched to the calculation domain of the nearly perfectly matched layers type have been introduced to simulate, when needed, semi-infinite or infinite media. The developed DG-FEM scheme has been verified by means of a comparison with analytical solutions and numerical results already published in the literature for simple geometrical configurations: Lamb's problem and plane wave nonlinear propagation.     

High order multi-moment constrained finite volume method. Part I: Basic formulation

A new high-order finite volume method based on local reconstruction is presented in this paper. The method, so-called the multi-moment constrained finite volume (MCV) method, uses the point values defined within single cell at equally spaced points as the model variables (or unknowns). The time evolution equations used to update the unknowns are derived from a set of constraint conditions imposed on multi kinds of moments, i.e. the cell-averaged value and the point-wise value of the state variable and its derivatives. The finite volume constraint on the cell-average guarantees the numerical conservativeness of the method. Most constraint conditions are imposed on the cell boundaries, where the numerical flux and its derivatives are solved as general Riemann problems. A multi-moment constrained Lagrange interpolation reconstruction for the demanded order of accuracy is constructed over single cell and converts the evolution equations of the moments to those of the unknowns. The presented method provides a general framework to construct efficient schemes of high orders. The basic formulations for hyperbolic conservation laws in 1- and 2D structured grids are detailed with the numerical results of widely used benchmark tests.     

Construction of very high order residual distribution schemes for steady inviscid flow problems on hybrid unstructured meshes

In this paper we consider the very high order approximation of solutions of the Euler equations. We present a systematic generalization of the residual distribution method of [1] to very high order of accuracy, by extending the preliminary work discussed in [2] to systems and hybrid meshes. We present extensive numerical validation for the third and fourth order cases with Lagrange finite elements. In particular, we demonstrate that we both have a non-oscillatory behavior, even for very strong shocks and complex flow patterns, and the expected accuracy on smooth problems.     

非线性LC电路方程的显式精确行波解

理论上考察了具有耗散的非线性LC电路中的行波. 借助于作者最近发展的精确求解非线性偏微分方程的扩展的双曲函数方法解析地研究了模拟非线性电路中冲击波的四阶耗散非线性波动方程. 一致地获得了丰富的显式精确解析行波解, 包括精确冲击波解和奇异的行波解, 和三角函数有理形式的周期波解. 关键词： LC电路')" href="#">非线性LC电路 非线性耗散波动方程 冲击波 周期波