Instabilities in the Chapman-Enskog Expansion and Hyperbolic Burnett Equations

It is well-known that the classical Chapman-Enskog procedure does not work at the level of Burnett equations (the next step after the Navier-Stokes equations). Roughly speaking, the reason is that the solutions of higher equations of hydrodynamics (Burnett's, etc.) become unstable with respect to short-wave perturbations. This problem was recently attacked by several authors who proposed different ways to deal with it. We present in this paper one of possible alternatives. First we deduce a criterion for hyperbolicity of Burnett equations for the general molecular model and show that this criterion is not fulfilled in most typical cases. Then we discuss in more detail the problem of truncation of the Chapman-Enskog expansion and show that the way of truncation is not unique. The general idea of changes of coordinates (based on analogy with the theory of dynamical systems) leads finally to nonlinear Hyperbolic Burnett Equations (HBEs) without using any information beyond the classical Burnett equations. It is proved that HBEs satisfy the linearized H-theorem. The linear version of the problem is studied in more detail, the complete Chapman-Enskog expansion is given for the linear case. A simplified proof of the Slemrod identity for Burnett coefficients is also given.     

A Three-component Extension of the Camassa–Holm Hierarchy

We introduce a bi-Hamiltonian hierarchy on the loop-algebra of endowed with a suitable Poisson pair. It gives rise to the usual Camassa–Holm (CH) hierarchy by means of a bi-Hamiltonian reduction, and its first nontrivial flow provides a three-component extension of the CH equation.     

驻波实验中一个经常出现的错误

对驻波演示实验中出现的立体驻波的原因及特点进行分析,并提出改进实验的建议.     

Variational multiscale element-free Galerkin method for 2D Burgers’ equation

A new numerical method, which is based on the coupling between variational multiscale method and meshfree methods, is developed for 2D Burgers’ equation with various values of Re. The proposed method takes full advantage of meshfree methods, therefore, no mesh generation and mesh recreation are involved. Meanwhile, compared with the variational multiscale finite element method, a strong assumption, that is, the fine scale vanishes identically over the element boundaries although non-zero within the elements, is not needed. Subsequently two problems which have an available analytical solution of their own are solved to analyze the convergence behaviour of the proposed method. Finally a 2D Burgers’ equation having large Re is solved and the results have also been compared with the ones computed by two other methods. The numerical results show that the proposed method can indeed obtain accurate numerical results for 2D Burgers’ equation having large Re, which does not refer to the choice of a proper stabilization parameter.     

Fast Gaussian wavepacket transforms and Gaussian beams for the Schrödinger equation

This paper introduces a wavepacket-transform-based Gaussian beam method for solving the Schrödinger equation. We focus on addressing two computational issues of the Gaussian beam method: how to generate a Gaussian beam representation for general initial conditions and how to perform long time propagation for any finite period of time. To address the first question, we introduce fast Gaussian wavepacket transforms and develop on top of them an efficient initialization algorithm for general initial conditions. Based on this new initialization algorithm, we address the second question by reinitializing the beam representation when the beams become too wide. Numerical examples in one, two, and three dimensions demonstrate the efficiency and accuracy of the proposed algorithms. The methodology can be readily generalized to deal with other semi-classical quantum mechanical problems.     

Quasi-particle excitations around a pair of half-quantum vortices in p- wave superconductors

Triplet superconductors such as Sr₂RuO₄ and Na_xCoO₂⋅yH₂O are now found to be p-wave (p_x ± ip_y) or f-wave ((p_x ± ip_y) cos cp_z) superconductors. It was phenomenologically suggested that in these p-wave or f-wave superconductors, a pair of half-quantum vortices (HQVs) becomes stable. Using the Bogoliubov-de Gennes equation, previously we have analyzed quasi-particle excitations around an HQV at one end of a d-soliton for simplicity. In next study, we will investigate the stability of the pair of HQV’s, which are connected by the d-soliton. For this purpose, we have developed a new numerical method to solve the Bogoliubov-de Gennes equation for two vortices state, using Mathieu functions.     

Stochastic Four-State Mechanochemical Model of F1-ATPase

F₁-ATPase, a part of ATP synthase, can synthesize and hydrolyze ATP moleculars in which the central γ-subunit rotates inside the α₃β₃ cylinder. Astochastic four-state mechanochemical coupling model of F$_{1}$-ATPase isstudied with the aid of the master equation. In this model, the ATPhydrolysis and synthesis are dependent on ATP, ADP, and Pi concentrations.The effects of ATP concentration, ADP concentration, and the external torqueon the occupation probability of binding-state, the rotation rate and thediffusion coefficient of F₁-ATPase are investigated. Moreover, theresults from this model are compared with experiments. The mechanochemicalmechanism F₁-ATPase is qualitatively explained by the model.     

Induced Lie Algebras of a Six-Dimensional Matrix Lie Algebra

By using a six-dimensional matrix Lie algebra [Y.F. Zhang and Y. Wang, Phys. Lett. A 360 （2006） 92], three induced Lie algebras are constructed. One of them is obtained by extending Lie bracket, the others are higher-dimensional complex Lie algebras constructed by using linear transformations. The equivalent Lie algebras of the later two with multi-component forms are obtained as well. As their applications, we derive an integrable coupling and quasi-Hamiltonian structure of the modified TC hierarchy of soliton equations.     

Numerical Complexiton Solutions of Complex KdV Equation

In this paper, we directly extend the applications of the Adomian decomposition method to investigate the complex KdV equation. By choosing different forms of wave functions as the initial values, three new types of realistic numerical solutions： numerical positon, negaton solution, and particularly the numerical analytical complexiton solution are obtained, which can rapidly converge to the exact ones obtained by Lou et al. Numerical simulation figures are used to illustrate the efficiency and accuracy of the proposed method.     

差分方程解对初值连续性问题

本文讨论了差分方程的解对初值的连续性问题，对较一般的方程给出了保证解对初值连续的充分条件，对线性方程（组）利用系数矩阵的特征值给出其解对初值连续的条件．