Subset representations and eigenvalues of the universal intertwining matrix

We solve a combinatorial question concerning eigenvalues of the universal intertwining endomorphism of a subset representation.     

An efficient nonpolynomial spline method for distributed order fractional subdiffusion equations

In this paper, we propose an efficient numerical method for a distributed order fractional subdiffusion problem using nonpolynomial spline approach. The solvability, stability, and convergence of the scheme are established rigorously, and it is shown that the spatial convergence order improves some previous work done. Simulation is then conducted to verify the accuracy of the proposed scheme as well as to compare with earlier work.     

基于Caputo导数的分数阶非线性振动系统响应计算

研究了含分数阶Caputo导数的非线性振动系统响应的数值计算方法。首先，由Caputo分数阶导数算子的叠加关系，得到含分数阶导数项非线性振动系统状态方程的标准形式。其次，基于Caputo导数与Riemann-Liouville导数和Grunwald-Letnikov导数间的关系，推导计算了Caputo导数的一般数值迭代格式。本文方法不要求状态方程中各分数阶导数阶数相等，弱化了已有算法中对分数阶导数阶数的限制，并可推广到多自由度的情形。随后，选择若干有解析解的算例验证了本文方法的正确性。最后，以多吸引子共存的分数阶Duffing振子系统为例，比较Caputo和GL两种算法所得结果，说明了用GL算法求解存在的问题。     

A simple three‐wave approximate Riemann solver for the Saint‐Venant‐Exner equations

Erosion and sediments transport processes have a great impact on industrial structures and on water quality. Despite its limitations, the Saint‐Venant‐Exner system is still (and for sure for some years) widely used in industrial codes to model the bedload sediment transport. In practice, its numerical resolution is mostly handled by a splitting technique that allows a weak coupling between hydraulic and morphodynamic distinct softwares but may suffer from important stability issues. In recent works, many authors proposed alternative methods based on a strong coupling that cure this problem but are not so trivial to implement in an industrial context. In this work, we then pursue 2 objectives. First, we propose a very simple scheme based on an approximate Riemann solver, respecting the strong coupling framework, and we demonstrate its stability and accuracy through a number of numerical test cases. However, second, we reinterpret our scheme as a splitting technique and we extend the purpose to propose what should be the minimal coupling that ensures the stability of the global numerical process in industrial codes, at least, when dealing with collocated finite volume method. The resulting splitting method is, up to our knowledge, the only one for which stability properties are fully demonstrated.     

The existence of periodic solutions for nonlinear beam equations on by a para‐differential method

This paper focuses on the construction of periodic solutions of nonlinear beam equations on the d‐dimensional tori. For a large set of frequencies, we demonstrate that an equivalent form of the nonlinear equations can be obtained by a para‐differential conjugation. Given the nonresonant conditions on each finite dimensional subspaces, it is shown that the periodic solutions can be constructed for the block diagonal equation by a classical iteration scheme.     

Initial-boundary problem for the 1-D Euler-Boltzmann equations in radiation hydrodynamics

We study the initial-boundary value problem for the one dimensional EulerBoltzmann equation with reflection boundary condition. For initial data with small total variation, we use a modified Glimm scheme to construct the global approximate solutions(U_(△t,d), I_(△t,d)) and prove that there is a subsequence of the approximate solutions which is convergent to the global solution.     

ISSDE: A Monte Carlo implicit simulation code based on Stratonovich SDE approach of Coulomb collision

A Monte Carlo implicit simulation program, Implicit Stratonovich Stochastic Differential Equations(ISSDE), is developed for solving stochastic differential equations(SDEs) that describe plasmas with Coulomb collision. The basic idea of the program is the stochastic equivalence between the Fokker–Planck equation and the Stratonovich SDEs. The splitting method is used to increase the numerical stability of the algorithm for dynamics of charged particles with Coulomb collision. The cases of Lorentzian plasma, Maxwellian plasma and arbitrary distribution function of background plasma have been considered. The adoption of the implicit midpoint method guarantees exactly the energy conservation for the diffusion term and thus improves the numerical stability compared with conventional Runge–Kutta methods. ISSDE is built with C++ and has standard interfaces and extensible modules. The slowing down processes of electron beams in unmagnetized plasma and relaxation process in magnetized plasma are studied using the ISSDE, which shows its correctness and reliability.     

Novel local smoothness indicators for improving the third‐order WENO scheme

The local smoothness indicators play an important role in the performance of a weighted essentially nonoscillatory (WENO) scheme. Due to having only 2 points available on each substencil, the local smoothness indicators calculated by conventional methods make the third‐order WENO scheme too dissipative. In this paper, we propose a different method to calculate the indicators by using all the 3 points on the global stencil of the third‐order WENO scheme. The numerical results demonstrate that the WENO scheme with the new indicators has less dissipation and better resolution than the conventional third‐order WENO scheme of Jiang and Shu for both smooth and discontinuous solutions.