High-Order Semi-Lagrangian WENO Schemes Based on Non-polynomial Space for the Vlasov Equation

In this paper,we present a semi-Lagrangian(SL)method based on a non-polynomial function space for solving the Vlasov equation.We fnd that a non-polynomial function based scheme is suitable to the specifcs of the target problems.To address issues that arise in phase space models of plasma problems,we develop a weighted essentially non-oscillatory(WENO)scheme using trigonometric polynomials.In particular,the non-polynomial WENO method is able to achieve improved accuracy near sharp gradients or discontinuities.Moreover,to obtain a high-order of accuracy in not only space but also time,it is proposed to apply a high-order splitting scheme in time.We aim to introduce the entire SL algorithm with high-order splitting in time and high-order WENO reconstruction in space to solve the Vlasov-Poisson system.Some numerical experiments are presented to demonstrate robustness of the proposed method in having a high-order of convergence and in capturing non-smooth solutions.A key observation is that the method can capture phase structure that require twice the resolution with a polynomial based method.In 6D,this would represent a signifcant savings.

High-order field theory and a weak Euler–Lagrange–Barut equation for classical relativistic particle-field systems

In both quantum and classical field systems,conservation laws such as the conservation of energy and momentum are widely regarded as fundamental properties.A broadly accepted approach to deriving conservation laws is built using Noether's method.However,this procedure is still unclear for relativistic particle-field systems where particles are regarded as classical world lines.In the present study,we establish a general manifestly covariant or geometric field theory for classical relativistic particle-field systems.In contrast to quantum systems,where particles are viewed as quantum fields,classical relativistic particle-field systems present specific challenges.These challenges arise from two sides.The first comes from the mass-shell constraint.To deal with the mass-shell constraint,the Euler–Lagrange–Barut(ELB)equation is used to determine the particle's world lines in the four-dimensional(4D)Minkowski space.Besides,the infinitesimal criterion,which is a differential equation in formal field theory,is reconstructed by an integro-differential form.The other difficulty is that fields and particles depend on heterogeneous manifolds.To overcome this challenge,we propose using a weak version of the ELB equation that allows us to connect local conservation laws and continuous symmetries in classical relativistic particle-field systems.By applying a weak ELB equation to classical relativistic particle-field systems,we can systematically derive local conservation laws by examining the underlying symmetries of the system.Our proposed approach provides a new perspective on understanding conservation laws in classical relativistic particle-field systems.

High-order harmonic generation of ZnO crystals in chirped and static electric fields

High harmonic generation in ZnO crystals under chirped single-color field and static electric field are investigated by solving the semiconductor Bloch equation(SBE). It is found that when the chirp pulse is introduced, the interference structure becomes obvious while the harmonic cutoff is not extended. Furthermore, the harmonic efficiency is improved when the static electric field is included. These phenomena are demonstrated by the classical recollision model in real space affected by the waveform of laser field and inversion symmetry. Specifically, the electron motion in k-space shows that the change of waveform and the destruction of the symmetry of the laser field lead to the incomplete X-structure of the crystal-momentum-resolved(k-resolved) inter-band harmonic spectrum. Furthermore, a pre-acceleration process in the solid four-step model is confirmed.

A high-order splitting scheme for the advection diffusion equation

A high-order splitting scheme for the advection-diffusion equation of pollutants is proposed in this paper. The multidimensional advection-diffusion equation is splitted into several one-dimensional equations that are solved by the scheme. Only three spatial grid points are needed in each direction and the scheme has fourth-order spatial accuracy. Several typically pure advection and advection-diffusion problems are simulated. Numerical results show that the accuracy of the scheme is much higher than that of the classical schemes and the scheme can he efficiently solved with little programming effort.

A high-order accurate wavelet method for solving Schrdinger equations with general nonlinearity

A sampling approximation for a function defined on a bounded interval is proposed by combining the Coiflet-type wavelet expansion and the boundary extension technique. Based on such a wavelet approximation scheme, a Galerkin procedure is developed for the spatial discretization of the generalized nonlinear Schr6dinger （NLS） equa- tions, and a system of ordinary differential equations for the time dependent unknowns is obtained. Then, the classical fourth-order explicit Runge-Kutta method is used to solve this semi-discretization system. To justify the present method, several widely considered problems are solved as the test examples, and the results demonstrate that the proposed wavelet algorithm has much better accuracy and a faster convergence rate in space than many existing numerical methods.

The Hirota bilinear method for the coupled Burgers equation and the high-order Boussinesq-Burgers equation

This paper studies the coupled Burgers equation and the high-order Boussinesq-Burgers equation. The Hirota bilinear method is applied to show that the two equations are completely integrable. Multiple-kink （soliton） solutions and multiple-singular-kink （soliton） solutions are derived for the two equations.

Analysis of the rogue waves in the blood based on the high-order NLS equations with variable coefficients

The research of rogue waves is an advanced field which has important practical and theoretical significances in mathematics,physics,biological fluid mechanics,oceanography,etc.Using the reductive perturbation theory and long wave approximation,the equations governing the movement of blood vessel walls and the flow of blood are transformed into high-order nonlinear Schrodinger(NLS)equations with variable coefficients.The third-order nonlinear Schrodinger equation is degenerated into a completely integrable Sasa–Satsuma equation(SSE)whose solutions can be used to approximately simulate the real rogue waves in the vessels.For the first time,we discuss the conditions for generating rogue waves in the blood vessels and effects of some physiological parameters on the rogue waves.Based on the traveling wave solutions of the fourth-order nonlinear Schrodinger equation,we analyze the effects of the higher order terms and the initial deformations of the blood vessel on the wave propagation and the displacement of the tube wall.Our results reveal that the amplitude of the rogue waves are proportional to the initial stretching ratio of the tube.The high-order nonlinear and dispersion terms lead to the distortion of the wave,while the initial deformation of the tube wall will influence the wave amplitude and wave steepness.

On the Behavior of Combination High-Order Compact Approximations with Preconditioned Methods in the Diffusion-Convection Equation

In this paper, a family of high-order compact finite difference methods in combination preconditioned methods are used for solution of the Diffusion-Convection equation. We developed numerical methods by replacing the time and space derivatives by compact finite-difference approximations. The system of resulting nonlinear finite difference equations are solved by preconditioned Krylov subspace methods. Numerical results are given to verify the behavior of high-order compact approximations in combination preconditioned methods for stability, convergence. Also, the accuracy and efficiency of the proposed scheme are considered.

A High-Order Scheme for Fractional Ordinary Differential Equations with the Caputo-Fabrizio Derivative

In this paper, we consider numerical solutions of fractional ordinary diferential equations with the Caputo-Fabrizio derivative, and construct and analyze a high-order time-stepping scheme for this equation. The proposed method makes use of quadratic interpolation function in sub-intervals, which allows to produce fourth-order convergence. A rigorous stability and convergence analysis of the proposed scheme is given. A series of numerical examples are presented to validate the theoretical claims. Traditionally a scheme having fourth-order convergence could only be obtained by using block-by-block technique. The advantage of our scheme is that the solution can be obtained step by step, which is cheaper than a block-by-block-based approach.

High-Order Local Discontinuous Galerkin Algorithm with Time Second-Order Schemes for the Two-Dimensional Nonlinear Fractional Diffusion Equation

In this article,some high-order local discontinuous Galerkin(LDG)schemes based on some second-order θ approximation formulas in time are presented to solve a two-dimen-sional nonlinear fractional diffusion equation.The unconditional stability of the LDG scheme is proved,and an a priori error estimate with O(h^(k+1)+At^(2))is derived,where k≥0 denotes the index of the basis function.Extensive numerical results with Q^(k)(k=0,1,2,3)elements are provided to confirm our theoretical results,which also show that the second-order convergence rate in time is not impacted by the changed parameter θ.

High-Order Finite Difference Method for Helmholtz Equation in Polar Coordinates

We present a fourth-order finite difference scheme for the Helmholtz equation in polar coordinates. We employ the finite difference format in the interior of the region and derive a nine-point fourth-order scheme. Specially, ghost points outside the region are applied to obtain the approximation for the Neumann boundary condition. We obtain the matrix form of the linear system and the sparsity of the coefficient matrix is favorable for the computation of the Helmholtz equation. The feasibility and accuracy of the method are validated by two test examples which have exact solutions.

On the High-Order Quasi Exactly Solvable Differential Equations

In this paper, we present a new method for solving a class of high-order quasi exactly solvable ordinary differential equations. With this method, the computed solution is expressed as a linear combination of the canonical polynomials associated with the given differential operator. An iterative algorithm summarizing the procedure is presented and its efficiency is demonstrated through considering two applied problems.

A High-Order Conservative Numerical Method for Gross-Pitaevskii Equation with Time-Varying Coefficients in Modeling BEC

We propose a high-order conservative method for the nonlinear Sehodinger/Gross-Pitaevskii equation with time- varying coefficients in modeling Bose Einstein condensation （BEC）. This scheme combined with the sixth-order compact finite difference method and the fourth-order average vector field method, finely describes the condensate wave function and physical characteristics in some small potential wells. Numerical experiments are presented to demonstrate that our numerical scheme is efficient by the comparison with the Fourier pseudo-spectral method. Moreover, it preserves several conservation laws well and even exactly under some specific conditions.

A Family of Global Attractors for a Class of Generalized Kirchhoff-Beam Equations

The initial boundary value problem for a class of high-order Beam equations with quasilinear and strongly damped terms is studied. Firstly, the existence and uniqueness of the global solution of the equation are proved by prior estimation and Galerkin finite element method. Then the bounded absorption set is obtained by prior estimation, and the family of global attractors for the high-order Kirchhoff-Beam equation is obtained. The Frechet differentiability of the solution semigroup is proved after the linearization of the equation, and the decay of the volume element of the linearization problem is further proved. Finally, the Hausdorff dimension and Fractal dimension of the family of global attractors are proved to be finite.

HIGH-ORDER ACCURATE AND HIGH-RESOLUTION UPWIND FINITE VOLUME SCHEME FOR SOLVING EULER/REYNOLDS-AVERAGED NAVIER-STOKES EQUATIONS

A high-order accurate explicit scheme is proposed for solving Euler/Reynolds-averaged Navier-Stokes equations for steady and unsteady flows, respectively. Baldwin-Lomax turbulence model is utilized to obtain the turbulent viscosity. For the explicit scheme, the Runge-Kutta time-stepping methods of third orders are used in time integration, and space discretization for the right-hand side (RHS) terms of semi-discrete equations is performed by third-order ENN scheme for inviscid terms and fourth-order compact difference for viscous terms. Numerical experiments suggest that the present scheme not only has a fairly rapid convergence rate, but also can generate a highly resolved approximation to numerical solution, even to unsteady problem.

A NEW HIGH-ORDER ACCURACY EXPLICIT DIFFERENCE SCHEME FOR SOLVING THREE-DIMENSIONAL PARABOLIC EQUATIONS

In this paper, a new three-level explicit difference scheme with high-order accuracy is proposed for solving three-dimensional parabolic equations. The stability condition is r = Delta t/Delta x(2) = Delta t/Delta gamma(2) = Delta t/Delta z(2) less than or equal to 1/4, and the truncation error is O(Delta t(2) + Delta x(4)).

Phase-coherence dynamics of frequency-comb emission via high-order harmonic generation in few-cycle pulse trains

Frequency-comb emission via high-order harmonic generation(HHG)provides an alternative method for the coherent vacuum ultraviolet(VUV)and extreme ultraviolet(XUV)radiation at ultrahigh repetition rates.In particular,the temporal and spectral features of the HHG were shown to carry profound insight into frequency-comb emission dynamics.Here we present an ab initio investigation of the temporal and spectral coherence of the frequency comb emitted in HHG of He atom driven by few-cycle pulse trains.We find that the emission of frequency combs features a destructive and constructive coherences caused by the phase interference of HHG,leading to suppression and enhancement of frequency-comb emission.The results reveal intriguing and substantially different nonlinear optical response behaviors for frequency-comb emission via HHG.The dynamical origin of frequency-comb emission is clarified by analyzing the phase coherence in HHG processes in detail.Our results provide fresh insight into the experimental realization of selective enhancement of frequency comb in the VUV–XUV regimes.

Second-Order Invariant Domain Preserving ALE Approximation of Euler Equations

An invariant domain preserving arbitrary Lagrangian-Eulerian method for solving non-linear hyperbolic systems is developed.The numerical scheme is explicit in time and the approximation in space is done with continuous finite elements.The method is made invar-iant domain preserving for the Euler equations using convex limiting and is tested on vari-ous benchmarks.

Calibration of quantitative rescattering model for simulating vortex high-order harmonic generation driven by Laguerre–Gaussian beam with nonzero orbital angular momentum

We calibrate the macroscopic vortex high-order harmonic generation(HHG)obtained by the quantitative rescattering(QRS)model to compute single-atom induced dipoles against that by solving the time-dependent Schr?dinger equation(TDSE).We show that the QRS perfectly agrees with the TDSE under the favorable phase-matching condition,and the QRS can accurately predict the main features in the spatial profiles of vortex HHG if the phase-matching condition is not good.We uncover that harmonic emissions from short and long trajectories are adjusted by the phase-matching condition through the time-frequency analysis and the QRS can simulate the vortex HHG accurately only when the interference between two trajectories is absent.This work confirms that it is an efficient way to employ the QRS model in the single-atom response for precisely simulating the macroscopic vortex HHG.

AENO:a Novel Reconstruction Method in Conjunction with ADER Schemes for Hyperbolic Equations

In this paper,we present a novel spatial reconstruction scheme,called AENO,that results from a special averaging of the ENO polynomial and its closest neighbour,while retaining the stencil direction decided by the ENO choice.A variant of the scheme,called m-AENO,results from averaging the modified ENO(m-ENO)polynomial and its closest neighbour.The concept is thoroughly assessed for the one-dimensional linear advection equation and for a one-dimensional non-linear hyperbolic system,in conjunction with the fully discrete,high-order ADER approach implemented up to fifth order of accuracy in both space and time.The results,as compared to the conventional ENO,m-ENO and WENO schemes,are very encouraging.Surprisingly,our results show that the L_(1)-errors of the novel AENO approach are the smallest for most cases considered.Crucially,for a chosen error size,AENO turns out to be the most efficient method of all five methods tested.