Generalized Birkhoffian representation of nonholonomic systems and its discrete variational algorithm

By using the discrete variational method, we study the numerical method of the general nonholonomic system in the generalized Birkhoffian framework, and construct a numerical method of generalized Birkhoffian equations called a self-adjoint-preserving algorithm. Numerical results show that it is reasonable to study the nonholonomic system by the structure-preserving algorithm in the generalized Birkhoffian framework.     

Generalized Birkhofflan representation of nonholonomic systems and its discrete variational algorithm

By using the discrete variational method,we study the numerical method of the general nonholonomic system in the generalized Birkhoffian framework,and construct a numerical method of generalized Birkhoffian equations called a self-adjoint-preserving algorithm.Numerical results show that it is reasonable to study the nonholonomic system by the structure-preserving algorithm in the generalized Birkhoffian framework.     

Structure-Preserving Algorithms for the Lorenz System

Based on a splitting method and a composition method, we construct some structure-preserving algorithms with first-order, second-order and fourth-order accuracy for a Lorenz system. By using the Liouville＇s formula, it is proven that the structure-preserving algorithms exactly preserve the volume of infinitesimal cube for the Lorenz system. Numerical experimental results illustrate that for the conservative Lorenz system, the qualitative behaviour of the trajectories described by the classical explicit fourth-order Runge-Kutta （RK4） method and the fifth-order Runge-Kutta-Fehlberg （RKF45） method is wrong, while the qualitative behaviour derived from the structure-preserving algorithms with different orders of accuracy is correct. Moreover, for the small dissipative Lorenz system, the norm errors of the structure-preserving algorithms in phase space axe less than those of the Runge-Kutta methods.     

Local structure-preserving methods for the generalized Rosenau-RLW-KdV equation with power law nonlinearity

Local structure-preserving algorithms including multi-symplectic, local energy-and momentum-preserving schemes are proposed for the generalized Rosenau–RLW–Kd V equation based on the multi-symplectic Hamiltonian formula of the equation. Each of the present algorithms holds a discrete conservation law in any time–space region. For the original problem subjected to appropriate boundary conditions, these algorithms will be globally conservative. Discrete fast Fourier transform makes a significant improvement to the computational efficiency of schemes. Numerical results show that the proposed algorithms have satisfactory performance in providing an accurate solution and preserving the discrete invariants.     

一维Gross-Pitaevskii方程的高阶紧致分裂步多辛格式

首先把一维Gross-Pitaevskli方程改写成多辛Hamiltonian系统的形式,把形式通过分裂变成2个子哈密尔顿系统.然后,对这些子系统用辛或者多辛算法进行离散.通过对子系统数值算法的不同组合方式,得到不同精度的具有多辛算法特征数值格式.这些格式不仅具有多辛格式、分裂步方法和高阶紧致格式的特征,而且是质量守恒的.数值实验验证了新格式的数值行为.     

一类特殊非完整力学系统的辛算法计算

对一类特殊的非完整力学系统的动力学性质进行数值研究,采用当前比较优越的保结构算法进行数值计算,并和传统的Runge-Kutta方法进行比较. 通过计算结果的比较而得出辛算法在这类特殊的非完整力学系统的数值计算中的优越性. 关键词： 非完整约束 辛算法 辛差分格式     

Birkhoffian Symplectic Scheme for a Quantum System

In this paper, a classical system of ordinary differential equations is built to describe a kind of n-dimensional quantum systems. The absorption spectrum and the density of the states for the system are defined from the points of quantum view and classical view. From the Birkhoffian form of the equations, a Birkhoffian symplectic scheme is derived for solving n-dimensional equations by using the generating function method. Besides the Birkhoffian structure- preserving, the new scheme is proven to preserve the discrete local energy conservation law of the system with zero vector f . Some numerical experiments for a 3-dimensional example show that the new scheme can simulate the general Birkhoffian system better than the implicit midpoint scheme, which is well known to be symplectic scheme for Hamiltonian system.     

VSIA:一种求解流体力学方程的多积分矩公式

使用变量间的积分矩公式给出一种构造流体力学方程数值格式的方法,即体积面积分平均方法(VSIA).该方法使用两种积分矩VIA(体积积分平均)和SIA(面积积分平均)给出一种完全的守恒型体积积分公式.基于一类守恒型半Lagrange输运算法CIP-CSL,能够清晰地构造出VSIA.对一般演化方程进行讨论,并应用于Burgers湍流、无粘可压流和不可压粘性流.     

Noncommutative Differential Calculus and Its Application on Discrete Spaces

We present the noncommutative differential calculus on the function space of the infinite set and construct a homotopy operator to prove the analogue of the Poincare lemma for the difference complex. Then the horizontal and vertical complexes are introduced with the total differential map and vertical exterior derivative. As the application of the differential calculus, we derive the schemes with the conservation of symplecticity and energy for Hamiltonian system and a two-dimensional integral models with infinite sequence of conserved currents. Then an Euler-Lagrange cohomology with symplectic structure-preserving is given in the discrete classical mechanics.     

对流占优扩散问题的特征线法-差分法计算格式

本文用特征线目的和有限差分目的相结合的数值目的来求解对流问题和对流占优扩散问题,提出了两个计算格式,并给出了数值例子。     

非定常欠膨胀射流的正格式数值模拟

 采用二阶正格式方法对非定常欠膨胀射流进行了数值模拟。将二维守恒方程的正格式方法推广到轴对称Euler方程组的求解,并对不同马赫数下的燃气射流进行了数值计算。计算结果表明,该方法能够较好地捕捉到包含膨胀波、入射激波、反射激波、马赫盘、射流边界以及三波点等复杂射流流场的波系结构,与实验照片反映的流动特征以及已有的数值结果相吻合。表明该方法对间断解具有较强的捕捉能力,在激波阵面上不会出现数值振荡。     

Structure-preserving algorithms for Birkhoffian systems

We propose a new approach to construct structure-preserving algorithms for Birkhoffian systems. First, the Pfaff–Birkhoff variational principle is discretized, and based on the discrete variational principle the discrete Birkhoffian equations are obtained. Then, taking the discrete equations as an algorithm, the corresponding discrete flow is proved to be symplectic. That means the algorithm preserves the symplectic structure of Birkhoffian systems. Simulation results of the given example indicate that structure-preserving algorithms obtained by this method have great advantage in conserving conserved quantities.     

Evaluation of numerical scattering in finite volume method for solving radiative transfer equation by a central laser incidence model

The numerical scattering caused by spatial discretization in finite volume method is discussed. Based on an analysis of the generation process of numerical scattering, a physical model of central laser incidence to a two-dimensional rectangle containing semitransparent medium is established to validate the numerical scattering, with Monte Carlo method as benchmark, in which numerical scattering does not exist. Numerical scattering will be affected by spatial grid number, spatial differential schemes and spectral absorption coefficient. With the spatial grid number increasing, numerical scattering will be decreased. The accuracy of the diamond scheme is the highest, and the exponential scheme is a bit lower, the lowest accuracy of the three schemes is the step scheme. The tendency of numerical scattering is reverse, i.e., the step scheme produces minimum numerical scattering, and exponential scheme produces more, while the diamond scheme produces maximum among three methods. When the bias of absorption efficient is high, the numerical scattering cannot be eliminated only by increasing the grid number. If we set the direction of laser incidence as central axis, it can be seen that numerical scattering distributed symmetry along the axis, which can be called as symmetrical cross-scattering. All of the three schemes show symmetrical cross-scattering.     

CO-ORDINATE TRANSFORMATIONS FOR SECOND ORDER SYSTEMS. PART I: GENERAL TRANSFORMATIONS

When the dynamics of any general second order system are cast in a state-space format, the initial choice of the state vector usually comprises one partition representing system displacements and another representing system velocities. Co-ordinate transformations can be defined which result in more general definitions of the state vector. This paper discusses the general case of co-ordinate transformations of state-space representations for second order systems. It identifies one extremely important subset of such coordinate transformations—namely the set of structure-preserving transformations for second order systems—and it highlights the importance of these. It shows that one particular structure-preserving transformation results in a new system characterized by real diagonal matrices and presents a forceful case that this structure-preserving transformation should be considered to be the fundamental definition for the characteristic behaviour of general second order systems—in preference to the eigenvalue-eigenvector solutions conventionally accepted.     

Basic function method A new numerical method on unstructured grids

A new numerical method-basic function method is proposed. This method can directly discrete differential operators on unstructured grids. By using the expansion of basic function to approach the exact function, the central and upwind schemes of derivative are constructed. By using the polynomial as basic function, applying the technique of flux splitting method and the combination of central and upwind schemes, the non-physical fluctuation near the shock wave is suppressed. The first-order basic function scheme of polynomial type for solving inviscid compressible flow numerically is constructed in this paper. Several numerical results of many typical examples for one-, two- and three-dimensional inviscid compressible steady flow illustrate that it is a new scheme with high accuracy and high resolution for shock wave. Especially, combining with the adaptive remeshing technique, the satisfactory results can be obtained by these schemes.     

自相似Euler方程的数值方法

Euler方程某些问题的解具有自相似特点，可以使用更准确的方法求解.提出了两种数值方法，分别称为自相似和准自相似方法，新方法可以使用现有守恒律方程的数值格式，无须设计特殊方法.对一维激波管问题、二维Riemann问题、激波反射以及激波折射问题进行了数值计算.对自相似Euler方程，一维计算结果显示数值解基本等同于精确解，二维结果也比现有文献计算的结果有更高的分辨率.对准自相似Euler方程，新方法可以求解不具有自相似性但接近自相似的问题，并在计算时间足够长时可以取得自相似Euler方程的效果.数值求解自相似Euler方程对自相似问题的研究，高分辨率、高精度格式的设计乃至Euler方程的精确解都有重要启示.      

High-order non-uniform grid schemes for numerical simulation of hypersonic boundary-layer stability and transition

The direct numerical simulation of receptivity, instability and transition of hypersonic boundary layers requires high-order accurate schemes because lower-order schemes do not have an adequate accuracy level to compute the large range of time and length scales in such flow fields. The main limiting factor in the application of high-order schemes to practical boundary-layer flow problems is the numerical instability of high-order boundary closure schemes on the wall. This paper presents a family of high-order non-uniform grid finite difference schemes with stable boundary closures for the direct numerical simulation of hypersonic boundary-layer transition. By using an appropriate grid stretching, and clustering grid points near the boundary, high-order schemes with stable boundary closures can be obtained. The order of the schemes ranges from first-order at the lowest, to the global spectral collocation method at the highest. The accuracy and stability of the new high-order numerical schemes is tested by numerical simulations of the linear wave equation and two-dimensional incompressible flat plate boundary layer flows. The high-order non-uniform-grid schemes (up to the 11th-order) are subsequently applied for the simulation of the receptivity of a hypersonic boundary layer to free stream disturbances over a blunt leading edge. The steady and unsteady results show that the new high-order schemes are stable and are able to produce high accuracy for computations of the nonlinear two-dimensional Navier–Stokes equations for the wall bounded supersonic flow.     

A class of hybrid DG/FV methods for conservation laws I: Basic formulation and one-dimensional systems

By comparing the discontinuous Galerkin (DG) and the finite volume (FV) methods, a concept of ‘static reconstruction’ and ‘dynamic reconstruction’ is introduced for high-order numerical methods. Based on the new concept, a class of hybrid DG/FV schemes is presented for one-dimensional conservation law using a ‘hybrid reconstruction’ approach. In the hybrid DG/FV schemes, the lower-order derivatives of a piecewise polynomial solution are computed locally in a cell by the DG method based on Taylor basis functions (called as ‘dynamic reconstruction’), while the higher-order derivatives are re-constructed by the ‘static reconstruction’ of the FV method, using the known lower-order derivatives in the cell itself and its adjacent neighboring cells. The hybrid DG/FV methods can greatly reduce CPU time and memory required by the traditional DG methods with the same order of accuracy on the same mesh, and they can be extended directly to unstructured and hybrid grids in two and three dimensions similar to the DG and/or FV methods. The hybrid DG/FV methods are applied to one-dimensional conservation law, including linear and non-linear scalar equation and Euler equations. In order to capture the strong shock waves without spurious oscillations, a simple shock detection approach is developed to mark ‘trouble cells’, and a moment limiter is adopted for higher-order schemes. The numerical results demonstrate the accuracy, and the super-convergence property is shown for the third-order hybrid DG/FV schemes. In addition, by analyzing the eigenvalues of the semi-discretized system in one dimension, we discuss the spectral properties of the hybrid DG/FV schemes to explain the super-convergence phenomenon.     

Invariantization of the Crank-Nicolson method for Burgers’ equation

In this article, a geometric technique to construct numerical schemes for partial differential equations (PDEs) that inherit Lie symmetries is proposed. The moving frame method enables one to adjust the numerical schemes in a geometric manner and systematically construct proper invariant versions of them. To illustrate the method, we study invariantization of the Crank-Nicolson scheme for Burgers’ equation. With careful choice of normalization equations, the invariantized schemes are shown to surpass the standard scheme, successfully removing numerical oscillation around sharp transition layers.     

用改进的耦合型Level Set方法计算一维双介质可压缩流动

用带有虚拟流体(Ghost Fluid)修正的Level Set方法计算了一维可压缩双介质流动,把描述流动的Euler方程和描述流体界面运动的Level Set方程耦合起来,得到一个整体的守恒律系统,应用高分辨率差分格式求解;为了解决流体界面附近的数值跳动问题,在界面附近引入了虚拟流体方法的Isobaric修正,并给出了算例.