Birkhoffian symplectic algorithms derived from Hamiltonian symplectic algorithms

In this paper, we focus on the construction of structure preserving algorithms for Birkhoffian systems, based on existing symplectic schemes for the Hamiltonian equations. The key of the method is to seek an invertible transformation which drives the Birkhoffian equations reduce to the Hamiltonian equations. When there exists such a transformation,applying the corresponding inverse map to symplectic discretization of the Hamiltonian equations, then resulting difference schemes are verified to be Birkhoffian symplectic for the original Birkhoffian equations. To illustrate the operation process of the method, we construct several desirable algorithms for the linear damped oscillator and the single pendulum with linear dissipation respectively. All of them exhibit excellent numerical behavior, especially in preserving conserved quantities.     

Multi-symplectic variational integrators for nonlinear Schrdinger equations with variable coefficients

In this paper, we propose a variational integrator for nonlinear Schrdinger equations with variable coefficients. It is shown that our variational integrator is naturally multi-symplectic. The discrete multi-symplectic structure of the integrator is presented by a multi-symplectic form formula that can be derived from the discrete Lagrangian boundary function. As two examples of nonlinear Schrdinger equations with variable coefficients, cubic nonlinear Schrdinger equations and Gross–Pitaevskii equations are extensively studied by the proposed integrator. Our numerical simulations demonstrate that the integrator is capable of preserving the mass, momentum, and energy conservation during time evolutions. Convergence tests are presented to verify that our integrator has second-order accuracy both in time and space.     

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

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

Meson-baryon s-wave resonances with strangeness -3

Starting from a consistent SU(6) extension of the Weinberg-Tomozawa (WT) meson-baryon chiral Lagrangian ( Phys. Rev. D 74, 034025 (2006)), we study the s-wave meson-baryon resonances in the strangeness S = - 3 and negative-parity sectors. Those resonances are generated by solving the Bethe-Salpeter equation with the WT interaction used as kernel. The considered mesons are those of the 35-SU(6)-plet, which includes the pseudoscalar (PS) octet of pions and the vector (V) nonet of the rho-meson. For baryons we consider the 56-SU(6)-plet, made of the 1/2⁺ octet of the nucleon and the 3/2⁺ decuplet of the Delta. Quantum numbers I(J ^P) = 0(3/2^-) are suggested for the experimental resonances Ω ^*(2250)^- and Ω ^*(2380)^-. Among other, resonances with I = 1 are found, which minimal quark content is sssˉl', being s the strange quark and l, l' any of the the light up or down quarks. A clear signal for such a pentaquark would be a baryonic resonance with strangeness -3 and electric charge -2 or 0, in proton charge units. We suggest looking for K ^- Ξ ^- resonances with masses around 2100 and 2240MeV in the sector 1(1/2^-), and for π ^± Ω ^- and K ^- Ξ ^*- resonances with masses around 2260MeV in the sector 1(3/2^-).     

A data-analysis method for identifying differential effects of time-delayed feedback forces and periodic driving forces in stochastic systems

In this work a method is developed for analyzing time series of periodically driven stochastic systems involving time-delayed feedback. The proposed data-analysis method yields dynamical models in terms of stochastic delay differential equations. On the basis of these dynamical models differential effects of driving forces and time-delayed feedback forces can be identified.     

Propagation and collision of soliton rings in quantum semiconductor plasmas

The intrinsic localization of electrostatic wave energies in quantum semiconductor plasmas can be described by solitary pulses. The collision properties of these pulses are investigated. In the present study, the fundamental model includes the quantum term, degenerate pressure of the plasma species, and the electron/hole exchange–correlation effects. In cylindrical geometry, using the extended Poincaré–Lighthill–Kuo (PLK) method, the Korteweg–de Vries (KdV) equations and the analytical phase shifts after the collision of two soliton rings are derived. Typical values for GaSb and GaN semiconductors are used to estimate the basic features of soliton rings. It is found that the pulses of GaSb semiconductor carry more energies than the pulses of GaN semiconductor. In addition, the degenerate pressure terms of electrons and holes have strong impact on the phase shift. The present theory may be useful to analyze the collision of localized coherent electrostatic waves in quantum semiconductor plasmas.     

Two-Level Defect-Correction Method for Steady Navier-Stokes Problem with Friction Boundary Conditions

In this paper, we present two-level defect-correction finite element method for steady Navier-Stokes equations at high Reynolds number with the friction boundary conditions, which results in a variational inequality problem of the second kind. Based on Taylor-Hood element, we solve a variational inequality problem of Navier-Stokes type on the coarse mesh and solve a variational inequality problem of Navier-Stokes type corresponding to Newton linearization on the fine mesh. The error estimates for the velocity in the $H^1$ norm and the pressure in the $L^2$ norm are derived. Finally, the numerical results are provided to confirm our theoretical analysis.     

Stability Analysis of a Fully Coupled Implicit Scheme for Inviscid Chemical Non-Equilibrium Flows

Von Neumann stability theory is applied to analyze the stability of a fully coupled implicit (FCI) scheme based on the lower-upper symmetric Gauss-Seidel (LU-SGS) method for inviscid chemical non-equilibrium flows. The FCI scheme shows excellent stability except the case of the flows involving strong recombination reactions, and can weaken or even eliminate the instability resulting from the stiffness problem, which occurs in the subsonic high-temperature region of the hypersonic flow field. In addition, when the full Jacobian of chemical source term is diagonalized, the stability of the FCI scheme relies heavily on the flow conditions. Especially in the case of high temperature and subsonic state, the CFL number satisfying the stability is very small. Moreover, we also consider the effect of the space step, and demonstrate that the stability of the FCI scheme with the diagonalized Jacobian can be improved by reducing the space step. Therefore, we propose an improved method on the grid distribution according to the flow conditions. Numerical tests validate sufficiently the foregoing analyses. Based on the improved grid, the CFL number can be quickly ramped up to large values for convergence acceleration.     

Effect of initial-state target polarization on single ionization of helium by 1-keV electron impact

We report new results of triple differential cross sections for the single ionization of helium by 1-KeV electron impact at the ejection energy of 10 eV. Investigations have been made for both the perpendicular plane and the plane perpendicular to the momentum transfer geometries. The present calculation is based on the three-Coulomb wave function. Here we have also incorporated the effect of target polarization in the initial state. A comparison is made between the present calculation with the results of other theoretical methods and a recent experiment [D黵r M, Dimopoulou C, Najjari B, Dorn A, Bartschat K, Bray I, Fursa D V, Chen Z, Madison D H and Ullrich J 2008 Phys. Rev. A 77 032717]. At an impact energy of 1 KeV, the target polarization is found to induce a substantial change of the cross section for the ionization process. We observe that the effect of target polarization plays a dominant role in deciding the shape of triple differential cross sections.     

New exact solutions of Einstein—Maxwell equations for magnetostatic fields

The symmetry reduction method based on the Fréchet derivative of differential operators is applied to investigate symmetries of the Einstein-Maxwell field equations for magnetostatic fields, which is a coupled system of nonlinear partial differential equations of the second order. The technique yields invariant transformations that reduce the given system of partial differential equations to a system of nonlinear ordinary differential equations. Some of the reduced systems are further studied to obtain the exact solutions.