Numerical Investigations on the Resonance Errors of Multiscale Discontinuous Galerkin Methods for One-Dimensional Stationary Schrödinger Equation

In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.

An efficient locally one-dimensional finite-difference time-domain method based on the conformal scheme

An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.

ELASTIC WAVE LOCALIZATION IN TWO-DIMENSIONAL PHONONIC CRYSTALS WITH ONE-DIMENSIONAL QUASI-PERIODICITY AND RANDOM DISORDER

The band structures of both in-plane and anti-plane elastic waves propagating in two-dimensional ordered and disordered （in one direction） phononic crystals are studied in this paper. The localization of wave propagation due to random disorder is discussed by introducing the concept of the localization factor that is calculated by the plane-wave-based transfer-matrix method. By treating the quasi-periodicity as the deviation from the periodicity in a special way, two kinds of quasi phononic crystal that has quasi-periodicity （Fibonacci sequence） in one direction and translational symmetry in the other direction are considered and the band structures are characterized by using localization factors. The results show that the localization factor is an effective parameter in characterizing the band gaps of two-dimensional perfect, randomly disordered and quasi-periodic phononic crystals. Band structures of the phononic crystals can be tuned by different random disorder or changing quasi-periodic parameters. The quasi phononic crystals exhibit more band gaps with narrower width than the ordered and randomly disordered systems.

One-Dimensional Localization of Elastic Waves in Rib-Stiffened Plates

The 1 dimensional localization of elastic waves in disordered periodic multi span rib stiffened plates is investigated. The transfer matrix method is employed to obtain the transfer matrix of the system, and the method for calculating the Lyapunov exponents in continuous dynamic systems presented by Wolf is used to determine the localization factor. As examples, the numerical results of the localization factors are given for a disordered rib stiffened plate. The effects of the degree of disorder of span...

A robust triaxial localization method of AE source using refraction path

Acoustic emission(AE)localization algorithms based on homogeneous media or single-velocity are less accurate when applied to the triaxial localization experiments.To the end,a robust triaxial localization method of AE source using refraction path is proposed.Firstly,the control equation of the refraction path is established according to the sensor coordinates and arrival times.Secondly,considering the influence of time-difference-of-arrival(TDOA)errors,the residual of the governing equation is calculated to estimate the equation weight.Thirdly,the refraction points in different directions are solved using Snell’s law and orthogonal constraints.Finally,the source coordinates are iteratively solved by weighted correction terms.The feasibility and accuracy of the proposed method are verified by pencil-lead breaking experiments.The simulation results show that the new method is almost unaffected by the refraction ratio,and always holds more stable and accurate positioning performance than the traditional method under different ratios and scales of TDOA outliers.

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.

Quantum Solitons and Localized Modes in a One-Dimensional Lattice Chain with Nonlinear Substrate Potential

In the classical lattice theory, solitons and localized modes can exist in many one-dimensional nonlinear lattice chains, however, in the quantum lattice theory, whether quantum solitons and localized modes can exist or not in the one-dimensional lattice chains is an interesting problem. By using the number state method and the Hartree approximation combined with the method of multiple scales, we investigate quantum solirons and localized modes in a one-dimensional lattice chain with the nonlinear substrate potential. It is shown that quantum solitons do exist in this nonlinear lattice chain, and at the boundary of the phonon Brillouin zone, quantum solitons become quantum localized modes, phonons are pinned to the lattice of the vicinity at the central position j = j0.

A study on the numerical prediction method for the vertical thermal structure in the Bohai Sea and the Huanghai Sea-I.One-dimensional numerical prediction model

In this paper, on the basis of the heat conduction equation without consideration of the advection and turbulence effects, one-dimensional model for describing surface sea temperature ( T1), bottom sea temperature ( Tt ) and the thickness of the upper homogeneous layer ( h ) is developed in terms of the dimensionless temperature θT and depth η and self-simulation function θT - f(η) of vertical temperature profile by means of historical temperature data.The results of trial prediction with our one-dimensional model on T, Th, h , the thickness and gradient of thermocline are satisfactory to some extent.

Localization and Steering of Light in One-Dimensional Parity-Time Symmetric Photonic Lattices

We theoretically study the propagation dynamics of input light in one-dimensional mixed linear-nonlinear photonic lattices with a complex parity-time symmetric potential. Numerical computation shows simultaneous localization and steering of the optical beam due to the asymmetric scatter and interplay between Kerr-type nonlinearity and PT symmetry. This may provide a feasible measure for manipulation light in optical communications, integrated optics and so on.

Localization and delocalization of a one-dimensional system coupled with the environment

We investigate several models of a one-dimensional chain coupling with surrounding atoms to elucidate disorder- induced delocalization in quantum wires, a peculiar behaviour against common wisdom. We show that the localization length is enhanced by disorder of side sites in the case of strong disorder, but in the case of weak disorder there is a plateau in this dependence. The above behaviour is the conjunct influence of the coupling to the surrounding atoms and the antiresonant effect. We also discuss different effects and their physical origin of different types of disorder in such systems. The numerical results show that coupling with the surrounding atoms can induce either the localization or delocalization effect depending on the values of parameters.

Robust Damage Detection and Localization Under Complex Environmental Conditions Using Singular Value Decomposition-based Feature Extraction and One-dimensional Convolutional Neural Network

Ultrasonic guided wave is an attractive monitoring technique for large-scale structures but is vulnerable to changes in environmental and operational conditions(EOC),which are inevitable in the normal inspection of civil and mechanical structures.This paper thus presents a robust guided wave-based method for damage detection and localization under complex environmental conditions by singular value decomposition-based feature extraction and one-dimensional convolutional neural network(1D-CNN).After singular value decomposition-based feature extraction processing,a temporal robust damage index(TRDI)is extracted,and the effect of EOCs is well removed.Hence,even for the signals with a very large temperature-varying range and low signal-to-noise ratios(SNRs),the final damage detection and localization accuracy retain perfect 100%.Verifications are conducted on two different experimental datasets.The first dataset consists of guided wave signals collected from a thin aluminum plate with artificial noises,and the second is a publicly available experimental dataset of guided wave signals acquired on a composite plate with a temperature ranging from 20℃to 60℃.It is demonstrated that the proposed method can detect and localize the damage accurately and rapidly,showing great potential for application in complex and unknown EOC.

LOCAL ONE-DIMENSIONAL ASE-I SCHEME FOR 2D DIFFUSION EQUATION

A local alternating segment explicit - implicit method for the solution of 2D diffusion equations is presented in this paper .The method is unconditionally stable and has the obvious property of parallelism. Some numerical experiments show the method is not only simple but also more accurate.

Algebraic Calculation Method of One-Dimensional Steady Compressible Gas Flow

This paper describes a new method of calculation of one-dimensional steady compressible gas flows in channels with possible heat and mass exchange through perforated sidewalls. The channel is divided into small elements of a finite size for which mass, energy and momentum conservation laws are written in the integral form, assuming linear distribution of the parameters along the length. As a result, the calculation is reduced to finding the roots of a quadratic algebraic equation, thus providing an alternative to numerical methods based on differential equations. The advantage of this method is its high tolerance to coarse discretization of the calculation area as well as its good applicability for transonic flow calculations.

Alternating Group Explicit Iterative Methods for One-Dimensional Advection-Diffusion Equation

The finite difference method such as alternating group iterative methods is useful in numerical method for evolutionary equations and this is the standard approach taken in this paper. Alternating group explicit (AGE) iterative methods for one-dimensional convection diffusion equations problems are given. The stability and convergence are analyzed by the linear method. Numerical results of the model problem are taken. Known test problems have been studied to demonstrate the accuracy of the method. Numerical results show that the behavior of the method with emphasis on treatment of boundary conditions is valuable.

Controlling the Bandgaps of One-Dimensional TiO2/SiO2, TiO2/SnO2, and SiO2/SnO2 Photonic Crystals Using the Transfer Matrix Method

One-dimensional photonic crystals (1D PhCs) have a unique ability to control the propagation of light waves, however certain classes of 1D oxides remain relatively unexplored for use as PhCs. Specifically, there has not been a comparative study of the three different 1D PhC structures to compare the influence of layer thickness, number, and refractive index on the ability of the PhCs to control light transmission. Herein, we use the transfer matrix method (TMM) to theoretically examine the transmission of 1D PhCs composed of layers of TiO2/SiO2, TiO2/SnO2, SiO2/SnO2, and combinations of the three with various top and bottom layer thicknesses to cover a substantial region of the electromagnetic spectrum (UV to NIR). With increasing layer numbers for TiO2/SiO2 and SiO2/SnO2, the edges became sharper and wider and the photonic bandgap width increased. Moreover, we demonstrated that PhCs with significantly thick TiO2/SiO2 layers had a high transmittance for a wide bandgap, allowing for wide-band optical filter applications. These different PhC architectures could enable a variety of applications, depending on the properties needed.

A local pseudo arc-length method for hyperbolic conservation laws

A local pseudo arc-length method（LPALM）for solving hyperbolic conservation laws is presented in this paper.The key idea of this method comes from the original arc-length method,through which the critical points are bypassed by transforming the computational space.The method is based on local changes of physical variables to choose the discontinuous stencil and introduce the pseudo arc-length parameter,and then transform the governing equations from physical space to arc-length space.In order to solve these equations in arc-length coordinate,it is necessary to combine the velocity of mesh points in the moving mesh method,and then convert the physical variable in arclength space back to physical space.Numerical examples have proved the effectiveness and generality of the new approach for linear equation,nonlinear equation and system of equations with discontinuous initial values.Non-oscillation solution can be obtained by adjusting the parameter and the mesh refinement number for problems containing both shock and rarefaction waves.

A MESHLESS LOCAL PETROV-GALERKIN METHOD FOR GEOMETRICALLY NONLINEAR PROBLEMS

Nonlinear formulations of the meshless local Petrov-Galerkin （MLPG） method are presented for geometrically nonlinear problems. The method requires no mesh in computation and therefore avoids mesh distortion difficulties in the large deformation analysis. The essential boundary conditions in the present formulation axe imposed by a penalty method. An incremental and iterative solution procedure is used to solve geometrically nonlinear problems. Several examples are presented to demonstrate the effectiveness of the method in geometrically nonlinear problems analysis. Numerical results show that the MLPG method is an effective one and that the values of the unknown variable are quite accurate.

Exact analytic solutions for an elliptic hole with asymmetric collinear cracks in a one-dimensional hexagonal quasi-crystal

Using the complex variable function method and the technique of conformal mapping, the anti-plane shear problem of an elliptic hole with asymmetric colfinear cracks in a one-dimensional hexagonal quasi-crystal is solved, and the exact analytic solutions of the stress intensity factors （SIFs） for mode Ⅲ problem are obtained. Under the limiting conditions, the present results reduce to the Griffith crack and many new results obtained as well, such as the circular hole with asymmetric collinear cracks, the elliptic hole with a straight crack, the mode T crack, the cross crack and so on. As far as the phonon field is concerned, these results, which play an important role in many practical and theoretical applications, are shown to be in good agreement with the classical results.

NUMERICAL ANALYSIS OF MINDLIN SHELL BY MESHLESS LOCAL PETROV-GALERKIN METHOD

The objectives of this study are to employ the meshless local Petrov-Galerkin method （MLPGM） to solve three-dimensional shell problems. The computational accuracy of MLPGM for shell problems is affected by many factors, including the dimension of compact support domain, the dimension of quadrture domain, the number of integral cells and the number of Gauss points. These factors＇ sensitivity analysis is to adopt the Taguchi experimental design technology and point out the dimension of the quadrature domain with the largest influence on the computational accuracy of the present MLPGM for shells and give out the optimum combination of these factors. A few examples are given to verify the reliability and good convergence of MLPGM for shell problems compared to the theoretical or the finite element results.

Shift control method for the local time at descendi ngnode based on sun-synchronous orbit satellite

This article analyzes the shift factors of the descending node local time for sun-synchronous satellites and proposes a shift control method to keep the local time shift within an allowance range. It is found that the satellite orbit design and the orbit injection deviation are the causes for the initial shift velocity, whereas the atmospheric drag and the sun gravitational perturbation produce the shift acceleration. To deal with these shift factors, a shift control method is put forward, through such methods as orbit variation design, orbit altitude, and inclination keeping control. The simulation experiment and practical application have proved the effectiveness of this control method.