A UNIFORM HIGH-ORDER METHOD FOR A SINGULAR PERTURBATION PROBLEM IN CONSERVATIVE FORM

A uniform high-order method is. presented for the numerical solution of a singular perturbation problem in conservative form. We firest replace the original second-order problem (1.1) by two equivalent first-order problems ( 1.4), i.e., the solution of (1.1) is a linear combination of the solutions of (1.4). Then we derive a uniformly O (hm+1) accurate scheme for the first-order problems (1.4), where m is an arbitrary nonnegative integer, so we can get a uniformly O (hm+1) accurate solution of the original problem (1.1) by relation (1.3). Some illustrative numerical results are also given.

Explicit High-Order Method to Solve Coupled Nonlinear Schrödinger Equations

Models of the coupled nonlinear Schrödinger equations submit various critical physical phenomena with a typical equation for optical fibres with linear refraction. In this article, we will presuppose the Compact Finite Difference method with Runge-Kutta of order 4 (explicit) method, which is sixth-order and fourth-order in space and time respectively, to solve coupled nonlinear Schrödinger equations. Many methods used to solve coupled nonlinear Schrödinger equations are second order in time and need to use extra-technique to rise up to fourth-order as Richardson Extrapolation technique. The scheme obtained is immediately fourth-order in one step. This approach is a conditionally stable method. The conserved quantities and the exact single soliton solution indicate the competence and accuracy of the article’s suggestion schemes. Furthermore, the article discusses the two solitons interaction dynamics.

High-Order Decoupled and Bound Preserving Local Discontinuous Galerkin Methods for a Class of Chemotaxis Models

In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.

New High-Order Numerical Methods for Hyperbolic Systems of Nonlinear PDEs with Uncertainties

In this paper,we develop new high-order numerical methods for hyperbolic systems of nonlinear partial differential equations(PDEs)with uncertainties.The new approach is realized in the semi-discrete finite-volume framework and is based on fifth-order weighted essentially non-oscillatory(WENO)interpolations in(multidimensional)random space combined with second-order piecewise linear reconstruction in physical space.Compared with spectral approximations in the random space,the presented methods are essentially non-oscillatory as they do not suffer from the Gibbs phenomenon while still achieving high-order accuracy.The new methods are tested on a number of numerical examples for both the Euler equations of gas dynamics and the Saint-Venant system of shallow-water equations.In the latter case,the methods are also proven to be well-balanced and positivity-preserving.

High-Order Discontinuous Galerkin Method for Hovering Rotor Simulations Based on a Rotating Reference Frame

An implicit higher ? order discontinuous Galerkin(DG) spatial discretization for the compressible Euler equations in a rotating frame of reference is presented and applied to a rotor in hover using hexahedral grids. Instead of auxiliary methods like grid adaptation,higher ? order simulations(fourth ? and fifth ? order accuracy) are adopted.Rigorous numerical experiments are carefully designed,conducted and analyzed. The results show generally excellent consistence with references and vigorously demonstrate the higher?order DG method's better performance in loading distribution computations and tip vortex capturing, with much fewer degrees of freedom(DoF). Detailed investigations on the outer boundary conditions for hovering rotors are presented as well. A simple but effective speed smooth procedure is developed specially for the DG method. Further results reveal that the rarely used pressure restriction for outlet speed has a considerable advantage over the extensively adopted vertical speed restriction.

HIGH-ORDER RUNGE-KUTTA DISCONTINUOUS GALERKIN FINITE ELEMENT METHOD FOR 2-D RESONATOR PROBLEM

The Runge-Kutta discontinuous Galerkin finite element method （RK-DGFEM） is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and it is explicit in the time domain. Consequently it is a best mixture of FEM and finite volume method （FVM）. RK-DGFEM can obtain local high-order accuracy by using high-order polynomial basis. Numerical experiments of transverse magnetic （TM） wave propagation in a 2-D resonator are performed. A high-order Lagrange polynomial basis is adopted. Numerical results agree well with analytical solution. And different order Lagrange interpolation polynomial basis impacts on simulation result accuracy are discussed. Computational results indicate that the accuracy is evidently improved when the order of interpolation basis is increased. Finally, L^2 errors of different order polynomial basis in RK-DGFEM are presented. Computational results show that L^2 error declines exponentially as the order of basis increases.

Observer-based robust high-order fully actuated attitude autopilot design for spinning glide-guided projectiles

This paper investigates the design of an attitude autopilot for a dual-channel controlled spinning glideguided projectile(SGGP),addressing model uncertainties and external disturbances.Based on fixed-time stable theory,a disturbance observer with integral sliding mode and adaptive techniques is proposed to mitigate total disturbance effects,irrespective of initial conditions.By introducing an error integral signal,the dynamics of the SGGP are transformed into two separate second-order fully actuated systems.Subsequently,employing the high-order fully actuated approach and a parametric approach,the nonlinear dynamics of the SGGP are recast into a constant linear closed-loop system,ensuring that the projectile's attitude asymptotically tracks the given goal with the desired eigenstructure.Under the proposed composite control framework,the ultimately uniformly bounded stability of the closed-loop system is rigorously demonstrated via the Lyapunov method.Validation of the effectiveness of the proposed attitude autopilot design is provided through extensive numerical simulations.

High-order discontinuous Galerkin method for applications to multicomponent and chemically reacting flows

This article focuses on the development of a discontinuous Galerkin (DG) method for simulations of multicomponent and chemically reacting flows. Compared to aerodynamic flow applications, in which DG methods have been successfully employed, DG simulations of chemically reacting flows introduce challenges that arise from flow unsteadiness, combustion, heat release, compressibility effects, shocks, and variations in thermodynamic properties. To address these challenges, algorithms are developed, including an entropy-bounded DG method, an entropy-residual shock indicator, and a new formulation of artificial viscosity. The performance and capabilities of the resulting DG method are demonstrated in several relevant applications, including shock/bubble interaction, turbulent combustion, and detonation. It is concluded that the developed DG method shows promising performance in application to multicomponent reacting flows. The paper concludes with a discussion of further research needs to enable the application of DG methods to more complex reacting flows.

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.

A force control high-order single-step-β method (HSM) for substructure pseudo-dynamic testing

This paper proposes a new technique which introduces the high-order single-step-β method(HSM)into the experimental study on the substructure pseudo-dynamic testing(SPDT).The technique is based on the proposed concept of equivalent shear stiffness which can meet the requirement of the HSM algorithm.A study is done to theoretically validate the technique by the numerical analysis of two-storey shear building structure,in comparison of the proposed substructure pseudo-dynamic testing algorithm with the central difference method(CDM).Then,a full-scale SPDT model,the three-storey frame-supported reinforced concrete short-limb masonry shear wall structure,is designed and tested to simulate the seismic response of the corresponding six-storey structure and verify the proposed force control HSM technique.Meanwhile,the techniques of both stiffness correction and force control are suggested to control algorithmic error,control error and measurement error.The results indicate that the force control HSM can be used in the full-scale multi-degree-of-freedom(MDOF)substructure pseudo-dynamic testing before descent segment of structure restoring force properties.

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.

Distributed wide field electromagnetic method based on high-order 2^(n) sequence pseudo random signal

To make three-dimensional electromagnetic exploration achievable,the distributed wide field electromagnetic method(WFEM)based on the high-order 2^(n) sequence pseudo-random signal is proposed and realized.In this method,only one set of high-order pseudo-random waveforms,which contains all target frequencies,is needed.Based on high-order sequence pseudo-random signal construction algorithm,the waveform can be customized according to different exploration tasks.And the receivers are independent with each other and dynamically adjust the acquisition parameters according to different requirements.A field test in the deep iron ore of Qihe−Yucheng showed that the distributed WFEM based on high-order pseudo-random signal realizes the high-efficiency acquisition of massive electromagnetic data in quite a short time.Compared with traditional controlled-source electromagnetic methods,the distributed WFEM is much more efficient.Distributed WFEM can be applied to the large scale and high-resolution exploration for deep resources and minerals.

Arbitrary High-Order Fully-Decoupled Numerical Schemes for Phase-Field Models of Two-Phase Incompressible Flows

Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.

HIGH-ORDER NYSTRM METHOD FOR THE EFIE OF EM SCATTERING PROBLEMS

Nystrm method is a new method for solving electromagnetic scattering problems. This paper gives the detailed description on high-order Nystrm method used for the electric field integral equation of electromagnetic scattering problems. The numerical solutions of two examples are correct compared with Method Of Moment(MOM).

Numerical modeling of wave equation by a truncated high-order finite-difference method

Finite-difference methods with high-order accuracy have been utilized to improve the precision of numerical solution for partial differential equations. However, the computation cost generally increases linearly with increased order of accuracy. Upon examination of the finite-difference formulas for the first-order and second-order derivatives, and the staggered finite-difference formulas for the first-order derivative, we examine the variation of finite-difference coefficients with accuracy order and note that there exist some very small coefficients. With the order increasing, the number of these small coefficients increases, however, the values decrease sharply. An error analysis demonstrates that omitting these small coefficients not only maintain approximately the same level of accuracy of finite difference but also reduce computational cost significantly. Moreover, it is easier to truncate for the high-order finite-difference formulas than for the pseudospectral for- mulas. Thus this study proposes a truncated high-order finite-difference method, and then demonstrates the efficiency and applicability of the method with some numerical examples.

A perspective on high-order methods in computational fluid dynamics

There has been an intensive international effort to develop high-order Computational Fluid Dynamics(CFD) methods into design tools in aerospace engineering during the last one and half decades. These methods offer the potential to significantly improve solution accuracy and efficiency for vortex dominated turbulent flows. Enough progresses have been made in algorithm development, mesh generation and parallel computing that these methods are on the verge of being applied in a production design environment. Since many review papers have been written on the subject, I decide to offer a personal perspective on the state-of-the-art in high-order CFD methods and the challenges that must be overcome.

Tunable spectral continuous shift of high-order harmonic generation in atoms by a plasmon-assisted shaping pulse

We delve into the phenomenon of high-order harmonic generation within a helium atom under the influence of a plasmon-assisted shaping pulse.Our findings reveal an intriguing manipulation of the frequency peak position in the harmonic emission by adjusting the absolute phase parameter within the frequency domain of the shaping pulse.This phenomenon holds potential significance for experimental setups necessitating precisely tuned single harmonics.Notably,we observe a modulated shift in the created harmonic photon energy,spanning an impressive range of 1.2 eV.This frequency peak shift is rooted in the asymmetry exhibited by the rising and falling edges of the laser pulse,directly influencing the position of the peak frequency emission.Our study quantifies the dependence of this tuning range and the asymmetry of the laser pulse,offering valuable insights into the underlying mechanisms driving this phenomenon.Furthermore,our investigation uncovers the emergence of semi-integer order harmonics as the phase parameter is altered.We attribute this discovery to the intricate interference between harmonics generated by the primary and secondary return cores.This observation introduces an innovative approach for generating semi-integer order harmonics,thus expanding our understanding of high-order harmonic generation.Ultimately,our work contributes to the broader comprehension of complex phenomena in laser-matter interactions and provides a foundation for harnessing these effects in various applications,particularly those involving precise spectral control and the generation of unique harmonic patterns.

A method of solving the stiffness problem in Biot＇s poroelastic equations using a staggered high-order finite-difference

In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Blot＇s poroelastic equations will be encountered. In the paper, a partition method is developed to solve the stiffness problem with a staggered high-order finite-difference. The method splits the Biot equations into two systems. One is stiff, and solved analytically, the other is nonstiff, and solved numerically by using a high-order staggered-grid finite-difference scheme. The time step is determined by the staggered finite-difference algorithm in solving the nonstiff equations, thus a coarse time step may be employed. Therefore, the computation efficiency and computational stability are improved greatly. Also a perfect by matched layer technology is used in the split method as absorbing boundary conditions. The numerical results are compared with the analytical results and those obtained from the conventional staggered-grid finite-difference method in a homogeneous model, respectively. They are in good agreement with each other. Finally, a slightly more complex model is investigated and compared with related equivalent model to illustrate the good performance of the staggered-grid finite-difference scheme in the partition method.

Elliptically polarized high-order harmonic generation of Ar atom in an intense laser field

High-order harmonic generation(HHG) of Ar atom in an elliptically polarized intense laser field is experimentally investigated in this work.Interestingly,the anomalous ellipticity dependence on the laser ellipticity(ε) in the lower-order harmonics is observed,specifically in the 13rd-order,which displays a maximal harmonic intensity at ε ≈ 0.1,rather than at ε = 0 as expected.This contradicts the general trend of harmonic yield,which typically decreases with the increase of laser ellipticity.In this study,we attribute this phenomenon to the disruption of the symmetry of the wave function by the Coulomb effect,leading to the generation of a harmonic with high ellipticity.This finding provides valuable insights into the behavior of elliptically polarized harmonics and opens up a potential way for exploring new applications in ultrafast spectroscopy and light–matter interactions.

High-order Bragg forward scattering and frequency shift of low-frequency underwater acoustic field by moving rough sea surface

Acoustic scattering modulation caused by an undulating sea surface on the space-time dimension seriously affects underwater detection and target recognition.Herein,underwater acoustic scattering modulation from a moving rough sea surface is studied based on integral equation and parabolic equation.And with the principles of grating and constructive interference,the mechanism of this acoustic scattering modulation is explained.The periodicity of the interference of moving rough sea surface will lead to the interference of the scattering field at a series of discrete angles,which will form comb-like and frequency-shift characteristics on the intensity and the frequency spectrum of the acoustic scattering field,respectively,which is a high-order Bragg scattering phenomenon.Unlike the conventional Doppler effect,the frequency shifts of the Bragg scattering phenomenon are multiples of the undulating sea surface frequency and are independent of the incident sound wave frequency.Therefore,even if a low-frequency underwater acoustic field is incident,it will produce obvious frequency shifts.Moreover,under the action of ideal sinusoidal waves,swells,fully grown wind waves,unsteady wind waves,or mixed waves,different moving rough sea surfaces create different acoustic scattering processes and possess different frequency shift characteristics.For the swell wave,which tends to be a single harmonic wave,the moving rough sea surface produces more obvious high-order scattering and frequency shifts.The same phenomena are observed on the sea surface under fully grown wind waves,however,the frequency shift slightly offsets the multiple peak frequencies of the wind wave spectrum.Comparing with the swell and fully-grown wind waves,the acoustic scattering and frequency shift are not obvious for the sea surface under unsteady wind waves.