High-order schemes for predicting computational aeroacoustic propagation with adaptive mesh refinement

High-order schemes based on block-structured adaptive mesh refinement method are prepared to solve computational aeroacoustic （CAA） problems with an aim at improving computational efficiency. A number of numerical issues associated with high-order schemes on an adaptively refined mesh, such as stability and accuracy are addressed. Several CAA benchmark problems are used to demonstrate the feasibility and efficiency of the approach.

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.

A high-order scheme based on lattice Boltzmann flux solver for viscous compressible flow simulations

In this paper,a high-order scheme based on the lattice Boltzmann flux solver(LBFS)is proposed to simulate viscous compressible flows.The flux reconstruction(FR)approach is adopted to implement the spatial discretization.The LBFS is employed to compute the inviscid flux by using the local reconstruction of the lattice Boltzmann equation solutions from macroscopic flow variables.Meanwhile,a switch function is used in LBFS to adjust the magnitude of the numerical viscosity.Thus,it is more beneficial to capture both strong shock waves and thin boundary layers.Moreover,the viscous flux is computed according to the local discontinuous Galerkin method.Some typical compressible viscous problems,including manufactured solution case,lid-driven cavity flow,supersonic flow around a cylinder and subsonic flow over a NACA0012 airfoil,are simulated to demonstrate the accuracy and robustness of the proposed FR-LBFS.

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 targeted essentially non-oscillatory scheme for two-fluid plasma model

The weakly ionized plasma flows in aerospace are commonly simulated by the single-fluid model,which cannot describe certain nonequilibrium phenomena by finite collisions of particles,decreasing the fidelity of the solution.Based on an alternative formulation of the targeted essentially non-oscillatory(TENO)scheme,a novel high-order numerical scheme is proposed to simulate the two-fluid plasmas problems.The numerical flux is constructed by the TENO interpolation of the solution and its derivatives,instead of being reconstructed from the physical flux.The present scheme is used to solve the two sets of Euler equations coupled with Maxwell's equations.The numerical methods are verified by several classical plasma problems.The results show that compared with the original TENO scheme,the present scheme can suppress the non-physical oscillations and reduce the numerical dissipation.

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.

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.

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.

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 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.

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.

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 Spatiotemporal Precision-Matching Taylor–Li Scheme for Time-Dependent Problems

Based on the Taylor series method and Li’s spatial differential method, a high-order hybrid Taylor–Li scheme is proposed.The results of a linear advection equation indicate that, using the initial values of the square-wave type, a result with thirdorder accuracy occurs. However, using initial values associated with the Gaussian function type, a result with very high precision appears. The study demonstrates that, when the order of the time integral is more than three, the corresponding optimal spatial difference order could be higher than six. The results indicate that the reason for why there is no improvement related to an order of spatial difference above six is the use of a time integral scheme that is not high enough. The author also proposes a recursive differential method to improve the Taylor–Li scheme’s computation speed. A more rapid and highprecision program than direct computation of the high-order space differential item is employed, and the computation speed is dramatically boosted. Based on a multiple-precision library, the ultrahigh-order Taylor–Li scheme can be used to solve the advection equation and Burgers’ equation.

A novel high-order scheme for numerical simulation of wake flow over helicopter rotors in hover

Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor.A new high-order scheme(we call it WENO-K)proposed by our research group is employed to minimize numerical dissipation and extended to numerical simulation of unsteady compressible viscous flows dominated by tip vortices over hovering rotors.WENO-K is referred to as an adaptively optimized WENO scheme with Gauss-Kriging reconstruction,and its advantage is to reduce dissipation in smooth regions of flow while preserving high-resolution around discontinuities.Here WENO-K scheme is adopted to reconstruct left and right state values within the Roe Riemann solver updating the inviscid fluxes on a structured dynamic overset grid.To minimize the accuracy loss for high-order reconstruction on artificial boundaries of overset grid,a method of multilayer fringes is proposed to carry out interpolation between background grid and blade grid.Massively parallel computing considering automatic load balance on averagely partitioned overset grid is developed to reduce the wall-clock time of an unsteady simulation.Numerical results for Caradonna-Tung(C-T)rotor in hover at the conditions of subsonic and transonic tip Mach numbers show that the thrust coefficient error for the result of WENO-K scheme is no more than 3%.Compared with WENO-JS scheme,WENO-K scheme achieves about 40%improvement on accuracy of predicting rotor thrust with only 4.1%extra computational cost.More importantly,WENO-K scheme can capture more sophisticated unsteady flow structures and resolve tip vortices to a larger wake age with an increment of about 270°compared to WENO-JS scheme.

Numerical Analysis of a High-Order Scheme for Nonlinear Fractional Differential Equations with Uniform Accuracy

We introduce a high-order numerical scheme for fractional ordinary differential equations with the Caputo derivative.The method is developed by dividing the domain into a number of subintervals,and applying the quadratic interpolation on each subinterval.The method is shown to be unconditionally stable,and for general nonlinear equations,the uniform sharp numerical order 3−νcan be rigorously proven for sufficiently smooth solutions at all time steps.The proof provides a gen-eral guide for proving the sharp order for higher-order schemes in the nonlinear case.Some numerical examples are given to validate our theoretical results.

Laser shaping and optical power limiting of pulsed Laguerre–Gaussian laser beams of high-order radial modes in fullerene C60

We study the strong nonlinear optical dynamics of nanosecond pulsed Laguerre–Gaussian laser beams of high-order radial modes with zero orbital angular momentum propagating in the fullerene C60molecular medium. It is found that the spatiotemporal profile of the incident pulsed Laguerre–Gaussian laser beam is strongly reshaped during its propagation in the C60molecular medium. The centrosymmetric temporal profile of the incident pulse gradually evolves into a noncentrosymmetric meniscus shape, and the on-axis pulse duration is clearly depressed. Furthermore, the field intensity is distinctly attenuated due to the field-intensity-dependent reverse saturable absorption, and clear optical power limiting behavior is observed for different orders of the input pulsed Laguerre–Gaussian laser beams before the takeover of the saturation effect;the lower the order of the Laguerre–Gaussian beam, the lower the energy transmittance.

Convergence to Steady-State Solutions of the New Type of High-Order Multi-resolution WENO Schemes: a Numerical Study

A new type of high-order multi-resolution weighted essentially non-oscillatory(WENO)schemes(Zhu and Shu in J Comput Phys,375:659-683,2018)is applied to solve for steady-state problems on structured meshes.Since the classical WENO schemes(Jiang and Shu in J Comput Phys,126:202-228,1996)might suffer from slight post-shock oscillations(which are responsible for the residue to hang at a truncation error level),this new type of high-order finite-difference and finite-volume multi-resolution WENO schemes is applied to control the slight post-shock oscillations and push the residue to settle down to machine zero in steady-state simulations.This new type of multi-resolution WENO schemes uses the same large stencils as that of the same order classical WENO schemes,could obtain fifth-order,seventh-order,and ninth-order in smooth regions,and could gradually degrade to first-order so as to suppress spurious oscillations near strong discontinuities.The linear weights of such new multi-resolution WENO schemes can be any positive numbers on the condition that their sum is one.This is the first time that a series of unequal-sized hierarchical central spatial stencils are used in designing high-order finitedifference and finite-volume WENO schemes for solving steady-state problems.In comparison with the classical fifth-order finite-difference and finite-volume WENO schemes,the residue of these new high-order multi-resolution WENO schemes can converge to a tiny number close to machine zero for some benchmark steady-state problems.

High-order maximum-principle-preserving and positivity-preserving weighted compact nonlinear schemes for hyperbolic conservation laws

In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws and the compressible Euler systems in both one and two dimensions.The main idea of the present method is to rewrite the scheme in a conservative form,and then define the local limiting parameters via case-by-case discussion.Smooth test problems are presented to demonstrate that the proposed MPP/PP WCNSs incorporating a third-order Runge-Kutta method can attain the desired order of accuracy.Other test problems with strong shocks and high pressure and density ratios are also conducted to testify the performance of the schemes.

A High-order Accuracy Explicit Difference Scheme with Branching Stability for Solving Higher-dimensional Heat-conduction Equation

为解决 4-dimensional heatconduction 方程的明确的差别计划是的高顺序的精确性 constructed.The 稳定性状况是 r = t/x2 = t/y2 =t/] z2 = t/w23/8，和截断错误是 O (t2 + x4 ) 。