A lattice Boltzmann model for the compressible Euler equations with second‐order accuracy

In this paper, we propose a new lattice Boltzmann model for the compressible Euler equations. The model is based on a three‐energy‐level and three‐speed lattice Boltzmann equation by using a method of higher moments of the equilibrium distribution functions. In order to obtain second‐order accuracy, we employ the ghost field distribution functions to remove the non‐physical viscous parts. We also use the conditions of the higher moment of the ghost field equilibrium distribution functions to obtain the equilibrium distribution functions. In the numerical examples, we compare the numerical results of this scheme with those obtained by other lattice Boltzmann models for the compressible Euler equations. Copyright © 2008 John Wiley & Sons, Ltd.     

Some remarks on X(6900)

The analysis of the LHCb data on \begin{document}$X(6900)$\end{document} found in the di-\begin{document}$J/\psi$\end{document} system was performed using a momentum-dependent Flatté-like parameterization. The use of the pole counting rule and spectral density function sum rule provides consistent evidence that both confining and molecular states are possible. Alternatively, the nature of \begin{document}$X(6900)$\end{document} cannot be distinguished if only the di-\begin{document}$J/\psi$\end{document} experimental data with current statistics are available. Nevertheless, we found that the lowest state in the di-\begin{document}$J/\psi$\end{document} system likely has the same quantum numbers as \begin{document}$X(6900)$\end{document}, and \begin{document}$X(6900)$\end{document} is probably not interpreted as a \begin{document}$J/\psi-\psi(2S)$\end{document} molecular state.     

Shocks in local supersonic zones

The results of numerical integration of the Euler equations governing two-dimensional and axisymmetric flows of an ideal (inviscid and non-heat-conducting) gas with local supersonic zones are presented. The subject of the study is the formation of shocks closing local supersonic zones. The flow in the vicinity of the initial point of the closing shock is calculated on embedded, successively refined grids with an accuracy much greater than that previously achieved. The calculations performed, together with the analysis of certain controversial issues, leave no doubt that it is the intersection of C ^?-characteristics proceeding from the sonic line inside the supersonic zone that is responsible for the closing shock formation.     

Optimal Runge–Kutta smoothers for the p-multigrid discontinuous Galerkin solution of the 1D Euler equations

This work presents a family of original Runge–Kutta methods specifically designed to be effective relaxation schemes in the numerical solution of the steady state solution of purely advective problems with a high-order accurate discontinuous Galerkin space discretization and a p-multigrid solution algorithm. The design criterion for the construction of the Runge–Kutta methods here developed is different form the one traditionally used to derive optimal Runge–Kutta smoothers for the h-multigrid algorithm, which are designed in order to provide a uniform damping of the error modes in the high-frequency range only. The method here proposed is instead designed in order to provide a variable amount of damping of the error modes over the entire frequency spectrum. The performance of the proposed schemes is assessed in the solution of the steady state quasi one-dimensional Euler equations for two test cases of increasing difficulty. Some preliminary results showing the performance for multidimensional applications are also presented.     

Comparison of the generalized Riemann solver and the gas-kinetic scheme for inviscid compressible flow simulations

The generalized Riemann problem (GRP) scheme for the Euler equations and gas-kinetic scheme (GKS) for the Boltzmann equation are two high resolution shock capturing schemes for fluid simulations. The difference is that one is based on the characteristics of the inviscid Euler equations and their wave interactions, and the other is based on the particle transport and collisions. The similarity between them is that both methods can use identical MUSCL-type initial reconstructions around a cell interface, and the spatial slopes on both sides of a cell interface involve in the gas evolution process and the construction of a time-dependent flux function. Although both methods have been applied successfully to the inviscid compressible flow computations, their performances have never been compared. Since both methods use the same initial reconstruction, any difference is solely coming from different underlying mechanism in their flux evaluation. Therefore, such a comparison is important to help us to understand the correspondence between physical modeling and numerical performances. Since GRP is so faithfully solving the inviscid Euler equations, the comparison can be also used to show the validity of solving the Euler equations itself. The numerical comparison shows that the GRP exhibits a slightly better computational efficiency, and has comparable accuracy with GKS for the Euler solutions in 1D case, but the GKS is more robust than GRP. For the 2D high Mach number flow simulations, the GKS is absent from the shock instability and converges to the steady state solutions faster than the GRP. The GRP has carbuncle phenomena, likes a cloud hanging over exact Riemann solvers. The GRP and GKS use different physical processes to describe the flow motion starting from a discontinuity. One is based on the assumption of equilibrium state with infinite number of particle collisions, and the other starts from the non-equilibrium free transport process to evolve into an equilibrium one through particle collisions. The different mechanism in the flux evaluation deviates their numerical performance. Through this study, we may conclude scientifically that it may NOT be valid to use the Euler equations as governing equations to construct numerical fluxes in a discretized space with limited cell resolution. To adapt the Navier–Stokes (NS) equations is NOT valid either because the NS equations describe the flow behavior on the hydrodynamic scale and have no any corresponding physics starting from a discontinuity. This fact alludes to the consistency of the Euler and Navier–Stokes equations with the continuum assumption and the necessity of a direct modeling of the physical process in the discretized space in the construction of numerical scheme when modeling very high Mach number flows. The development of numerical algorithm is similar to the modeling process in deriving the governing equations, but the control volume here cannot be shrunk to zero.     

Numerical methods for computing ground states and dynamics of nonlinear relativistic Hartree equation for boson stars

Efficient and accurate numerical methods are presented for computing ground states and dynamics of the three-dimensional (3D) nonlinear relativistic Hartree equation both without and with an external potential. This equation was derived recently for describing the mean field dynamics of boson stars. In its numerics, due to the appearance of pseudodifferential operator which is defined in phase space via symbol, spectral method is more suitable for the discretization in space than other numerical methods such as finite difference method, etc. For computing ground states, a backward Euler sine pseudospectral (BESP) method is proposed based on a gradient flow with discrete normalization; and respectively, for computing dynamics, a time-splitting sine pseudospectral (TSSP) method is presented based on a splitting technique to decouple the nonlinearity. Both BESP and TSSP are efficient in computation via discrete sine transform, and are of spectral accuracy in spatial discretization. TSSP is of second-order accuracy in temporal discretization and conserves the normalization in discretized level. In addition, when the external potential and initial data for dynamics are spherically symmetric, the original 3D problem collapses to a quasi-1D problem, for which both BESP and TSSP methods are extended successfully with a proper change of variables. Finally, extensive numerical results are reported to demonstrate the spectral accuracy of the methods and to show very interesting and complicated phenomena in the mean field dynamics of boson stars.     

利用偏置磁极周期会切磁铁产生边聚焦场

 利用理论分析和数值计算的方法,研究了偏置磁极周期会切磁铁产生的、可用于带状电子束宽边聚焦的边聚焦场。结果表明:增加偏置长度,减小电子束通道宽度,增大磁极轴向长度可以在基本不改变边聚焦场在x方向上分布特性的前提下提高幅值;增大电子束通道的高度,增大磁极厚度均可以减小边聚焦场在x方向上分布曲线的曲率,但同时其幅值也会降低,该降低可利用前述方法予以补偿。在对带状电子束宽边聚焦进行束匹配时,可以先进行曲率匹配,再进行幅值匹配。在进行参数选择时,应合理选择束通道高度和磁极厚度的取值,以避免束通道内的边聚焦场在x方向上的分布出现曲率反向。     

Vortex structures with complex points singularities in two-dimensional Euler equations. New exact solutions

In this work we present a new class of exact stationary solutions for two-dimensional (2D) Euler equations. Unlike already known solutions, the new ones contain complex singularities. We consider point singularities which have a vector field index greater than 1 as complex. For example, the dipole singularity is complex because its index is equal to 2. We present in explicit form a large class of exact localized stationary solutions for 2D Euler equations with a singularity whose index is equal to 3. The solutions obtained are expressed in terms of elementary functions. These solutions represent a complex singularity point surrounded by a vortex satellite structure. We also discuss the motion equation of singularities and conditions for singularity point stationarity which provide the stationarity of the complex vortex configuration.     

Perturbational blowup solutions to the compressible 1-dimensional Euler equations

We construct non-radially symmetry solutions for the compressible 1-dimensional adiabatic Euler equations in this Letter. In detail, we perturb the linear velocity with a drifting term:

(1)     

Dynamical behaviors of a discrete HIV‐1 virus model with bilinear infective rate

We establish a discrete virus dynamic model by discretizing a continuous HIV‐1 virus model with bilinear infective rate using ‘hybrid’ Euler method. We discuss not only the existence and global stability of the uninfected equilibrium but also the existence and local stability of the infected equilibrium. We prove that there exists a crucial value similar to that of the continuous HIV‐1 virus dynamics, which is called the basic reproductive ratio of the virus. If the basic reproductive ratio of the virus is less than one, the uninfected equilibrium is globally asymptotically stable. If the basic reproductive ratio of the virus is larger than one, the infected equilibrium exists and is locally stable. Moreover, we consider the permanence for such a system by constructing a Lyapunov function v_n. Copyright © 2013 John Wiley & Sons, Ltd.