ACCURATE COMPUTATION FOR IRREDUCIBLE NONSINGULAR M-MATRIX

Computing the eigenvalue of smallest modulus and its corresponding eigneveclor of an irreducible nonsingular M-matrix A is considered, It is shown that if the entries of A are known with high relative accuracy, its eigenvalue of smallest modulus and each component of the corresponding eigenvector will be determined to much higher accuracy than the standard perturbation theory suggests. An algorithm is presented to compute them with a small componentwise backward error, which is consistent with the perturbation results.

An Initial Perturbation Method for the Multiscale Singular Vector in Global Ensemble Prediction

Ensemble prediction is widely used to represent the uncertainty of single deterministic Numerical Weather Prediction(NWP) caused by errors in initial conditions(ICs). The traditional Singular Vector(SV) initial perturbation method tends only to capture synoptic scale initial uncertainty rather than mesoscale uncertainty in global ensemble prediction. To address this issue, a multiscale SV initial perturbation method based on the China Meteorological Administration Global Ensemble Prediction System(CMA-GEPS) is proposed to quantify multiscale initial uncertainty. The multiscale SV initial perturbation approach entails calculating multiscale SVs at different resolutions with multiple linearized physical processes to capture fast-growing perturbations from mesoscale to synoptic scale in target areas and combining these SVs by using a Gaussian sampling method with amplitude coefficients to generate initial perturbations. Following that, the energy norm,energy spectrum, and structure of multiscale SVs and their impact on GEPS are analyzed based on a batch experiment in different seasons. The results show that the multiscale SV initial perturbations can possess more energy and capture more mesoscale uncertainties than the traditional single-SV method. Meanwhile, multiscale SV initial perturbations can reflect the strongest dynamical instability in target areas. Their performances in global ensemble prediction when compared to single-scale SVs are shown to(i) improve the relationship between the ensemble spread and the root-mean-square error and(ii) provide a better probability forecast skill for atmospheric circulation during the late forecast period and for short-to medium-range precipitation. This study provides scientific evidence and application foundations for the design and development of a multiscale SV initial perturbation method for the GEPS.

Periodic signal extraction of GNSS height time series based on adaptive singular spectrum analysis

Singular spectrum analysis is widely used in geodetic time series analysis.However,when extracting time-varying periodic signals from a large number of Global Navigation Satellite System(GNSS)time series,the selection of appropriate embedding window size and principal components makes this method cumbersome and inefficient.To improve the efficiency and accuracy of singular spectrum analysis,this paper proposes an adaptive singular spectrum analysis method by combining spectrum analysis with a new trace matrix.The running time and correlation analysis indicate that the proposed method can adaptively set the embedding window size to extract the time-varying periodic signals from GNSS time series,and the extraction efficiency of a single time series is six times that of singular spectrum analysis.The method is also accurate and more suitable for time-varying periodic signal analysis of global GNSS sites.

THE ABSENCE OF SINGULAR CONTINUOUS SPECTRUM FOR PERTURBED JACOBI OPERATORS

This paper is mainly about the spectral properties of a class of Jacobi operators(H_(c,b)u)(n)=c_(n)u(n+1)+c_(n-1)u(n-1)+b_(n)u(n),.where∣c_(n)−1∣=O(n^(−α))and b_(n)=O(n^(−1)).We will show that,forα≥1,the singular continuous spectrum of the operator is empty.

Wavelet Multi-Resolution Interpolation Galerkin Method for Linear Singularly Perturbed Boundary Value Problems

In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.

Adaptive Optimal Output Regulation of Interconnected Singularly Perturbed Systems With Application to Power Systems

This article studies the adaptive optimal output regulation problem for a class of interconnected singularly perturbed systems(SPSs) with unknown dynamics based on reinforcement learning(RL).Taking into account the slow and fast characteristics among system states,the interconnected SPS is decomposed into the slow time-scale dynamics and the fast timescale dynamics through singular perturbation theory.For the fast time-scale dynamics with interconnections,we devise a decentralized optimal control strategy by selecting appropriate weight matrices in the cost function.For the slow time-scale dynamics with unknown system parameters,an off-policy RL algorithm with convergence guarantee is given to learn the optimal control strategy in terms of measurement data.By combining the slow and fast controllers,we establish the composite decentralized adaptive optimal output regulator,and rigorously analyze the stability and optimality of the closed-loop system.The proposed decomposition design not only bypasses the numerical stiffness but also alleviates the high-dimensionality.The efficacy of the proposed methodology is validated by a load-frequency control application of a two-area power system.

A Full Predictor-Corrector Finite Element Method for the One-Dimensional Heat Equation with Time-Dependent Singularities

The energy norm convergence rate of the finite element solution of the heat equation is reduced by the time-regularity of the exact solution. This paper presents an adaptive finite element treatment of time-dependent singularities on the one-dimensional heat equation. The method is based on a Fourier decomposition of the solution and an extraction formula of the coefficients of the singularities coupled with a predictor-corrector algorithm. The method recovers the optimal convergence rate of the finite element method on a quasi-uniform mesh refinement. Numerical results are carried out to show the efficiency of the method.

Global Existence and Decay of Solutions for a Class of a Pseudo-Parabolic Equation with Singular Potential and Logarithmic Nonlocal Source

This article investigates the well posedness and asymptotic behavior of Neumann initial boundary value problems for a class of pseudo-parabolic equations with singular potential and logarithmic nonlinearity. By utilizing cut-off techniques and combining with the Faedo Galerkin approximation method, local solvability was established. Based on the potential well method and Hardy Sobolev inequality, derive the global existence of the solution. In addition, we also obtained the results of decay.

Discrete vortex bound states with a van Hove singularity in the vicinity of the Fermi level

A theoretical study on discrete vortex bound states is carried out near a vortex core in the presence of a van Hove singularity(VHS) near the Fermi level by solving Bogoliubov–de Gennes(Bd G) equations. When the VHS lies exactly at the Fermi level and also at the middle of the band, a zero-energy state and other higher-energy states whose energy ratios follow integer numbers emerge. These discrete vortex bound state peaks undergo a splitting behavior when the VHS or Fermi level moves away from the middle of the band. Such splitting behavior will eventually lead to a new arrangement of quantized vortex core states whose energy ratios follow half-odd-integer numbers.

Local singularity and S–A methods for analyzing ore-producing anomalies in the Jianbiannongchang area of Heilongjiang,China

The Heilongjiang Jianbiannongchang area is located at the confluence of the Great and Lesser Xing’an Ranges.This area has a complex magmatic and tectonic evolutionary history that has resulted in a complex and diverse geological background for mineralization.In this study,isometric logarithmic ratio(ILR)transformations of Au,Cu,Pb,Zn,and Sb contents were performed in the1:50,000 soil geochemical data of the Jianbiannongchang area.Robust principal component analysis(RPCA)was conducted based on ILR transformation.The local singularity and spectrum-area(S-A)methods were used to extract information on mineralogic anomalies.The results showed that:(1)the transformed data eliminated the influence of the original data closure effect,and the PC1and PC2 information obtained by applying RPCA reflected ore-producing element anomalies dominated by Au and Cu.(2)The local singularity method can enhance the information of the local strong and weak slow anomalies.After performing local singularity analysis on PC1 and PC2,the obtained local anomalies reflected the local singularity spatial anomaly patterns related to Cu and Au mineralization in this area,which is an effective method for trapping ore-producing anomalies.(3)Furthermore,the composite anomaly decomposition of PC1 and PC2 was performed using the S-A method,and the screened anomalous and background fields reflect the ore-producing anomalies related to Cu and Au mineralization.This information is in agreement with known Cu and Au mineralization.(4)The geochemical anomalies with mineralization potential were obtained outside the known mineralization sites by integrating the information of oreproducing anomalies extracted by the local singularity and S-A methods,providing the theoretical basis and exploration direction for future exploration in the study area.

Self-Awareness,a Singularity of AI

Self-awareness,or self-consciousness,refers to reflective recognition of the existence of“subjective-self”.Every person has a subjective-self from which the person observes and interacts with the world.In this article,we argue that self-consciousness is an enigmatic phenomenon unique in human intelligence.Unlike many other intelligent and conscious capabilities,self-consciousness is not possible to be achieved in electronic computers and robots.Self-consciousness is an odd-point of human intelligence and a singularity of artificial intelligence(AI).Man-made intelligence through software is not capable of self-consciousness;therefore,robots will never become a newly created species.Because of the lack of self-awareness,AI software,such as Watson,Alpha-zero,ChatGPT,and PaLM,will remain a tool of humans and will not dominate the human society no matter how smart it is.This singularity of AI makes us re-think humbly what the future AI is like,what kind of robots we are going to deal with,and the blessing and threat of AI on humanity.

High impedance fault detection in distribution network based on S-transform and average singular entropy

When a high impedance fault(HIF)occurs in a distribution network,the detection efficiency of traditional protection devices is strongly limited by the weak fault information.In this study,a method based on S-transform(ST)and average singular entropy(ASE)is proposed to identify HIFs.First,a wavelet packet transform(WPT)was applied to extract the feature frequency band.Thereafter,the ST was investigated in each half cycle.Afterwards,the obtained time-frequency matrix was denoised by singular value decomposition(SVD),followed by the calculation of the ASE index.Finally,an appropriate threshold was selected to detect the HIFs.The advantages of this method are the ability of fine band division,adaptive time-frequency transformation,and quantitative expression of signal complexity.The performance of the proposed method was verified by simulated and field data,and further analysis revealed that it could still achieve good results under different conditions.

A Dimension-Splitting Variational Multiscale Element-Free Galerkin Method for Three-Dimensional Singularly Perturbed Convection-Diffusion Problems

By introducing the dimensional splitting(DS)method into the multiscale interpolating element-free Galerkin(VMIEFG)method,a dimension-splitting multiscale interpolating element-free Galerkin(DS-VMIEFG)method is proposed for three-dimensional(3D)singular perturbed convection-diffusion(SPCD)problems.In the DSVMIEFG method,the 3D problem is decomposed into a series of 2D problems by the DS method,and the discrete equations on the 2D splitting surface are obtained by the VMIEFG method.The improved interpolation-type moving least squares(IIMLS)method is used to construct shape functions in the weak form and to combine 2D discrete equations into a global system of discrete equations for the three-dimensional SPCD problems.The solved numerical example verifies the effectiveness of the method in this paper for the 3D SPCD problems.The numerical solution will gradually converge to the analytical solution with the increase in the number of nodes.For extremely small singular diffusion coefficients,the numerical solution will avoid numerical oscillation and has high computational stability.

THE SINGULAR CONVERGENCE OF A CHEMOTAXIS-FLUID SYSTEM MODELING CORAL FERTILIZATION

The singular convergence of a chemotaxis-fluid system modeling coral fertilization is justified in spatial dimension three.More precisely,it is shown that a solution of parabolic-parabolic type chemotaxis-fluid system modeling coral fertilization■converges to that of the parabolic-elliptic type chemotaxis-fluid system modeling coral fertiliz ation■in a certain Fourier-Herz space asε^(-1)→0.

THE UNIFORM CONVERGENCE OF A DG METHOD FOR A SINGULARLY PERTURBED VOLTERRA INTEGRO-DIFFERENTIAL EQUATION

The purpose of this work is to implement a discontinuous Galerkin(DG)method with a one-sided flux for a singularly perturbed Volterra integro-differential equation(VIDE)with a smooth kernel.First,the regularity property and a decomposition of the exact solution of the singularly perturbed VIDE with the initial condition are provided.Then the existence and uniqueness of the DG solution are proven.Then some appropriate projection-type interpolation operators and their corresponding approximation properties are established.Based on the decomposition of the exact solution and the approximation properties of the projection type interpolants,the DG method achieves the uniform convergence in the L2 norm with respect to the singular perturbation parameter e when the space of polynomials with degree p is used.A numerical experiment validates the theoretical results.Furthermore,an ultra-convergence order 2p+1 at the nodes for the one-sided flux,uniform with respect to the singular perturbation parameter e,is observed numerically.

Singularity of Two Kinds of Quadcyclic Peacock Graphs

Let G be a graph. G is singular if and only if the adjacency matrix of graph G is singular. The adjacency matrix of graph G is singular if and only if there is at least one zero eigenvalue. The study of the singularity of graphs is of great significance for better characterizing the properties of graphs. The following definitions are given. There are 4 paths, the starting points of the four paths are bonded into one point, and the ending point of each path is bonded to a cycle respectively, so this graph is called a kind of quadcyclic peacock graph. And in this kind of quadcyclic peacock graph assuming the number of points on the four cycles is a1, a2, a3, a4, and the number of points on the four paths is s1, s2, s3, s4, respectively. This type of graph is denoted by γ (a1, a2, a3, a4, s1, s2, s3, s4), called γ graph. And let γ (a1, a2, a3, a4, 1, 1, 1, 1) = δ (a1, a2, a3, a4), this type four cycles peacock graph called δ graph. In this paper, we give the necessary and sufficient conditions for the singularity of γ graph and δ graph.

Evolution of polarization singularities accompanied by avoided crossing in plasmonic system

The evolution of polarization singularities supported in a one-dimensional periodic plasmonic system is studied.The lateral inversion symmetry of the system,which breaks the in-plane inversion symmetry and up-down mirror symmetry simultaneously,yields abundant polarization states.A complete evolution process with geometry for the polarization states is traced.In the evolution,circularly polarized points(C points)can stem from 3 different processes.In addition to the previously reported processes occurring in an isolated band,a new type of C point appearing in two bands simultaneously due to the avoided band crossing,is observed.Unlike the dielectric system with a similar structure which only supports at-Γbound states in the continuum(BICs),accidental BICs off theΓpoint are realized in this plasmonic system.This work provides a new scheme of polarization manipulation for the plasmonic systems.

THE SINGULAR LIMIT OF SECOND-GRADE FLUID EQUATIONS IN A 2D EXTERIOR DOMAIN

In this paper, we consider the second-grade fluid equations in a 2D exterior domain satisfying the non-slip boundary conditions. The second-grade fluid model is a wellknown non-Newtonian fluid model, with two parameters: α, which represents the length-scale,while ν > 0 corresponds to the viscosity. We prove that, as ν, α tend to zero, the solution of the second-grade fluid equations with suitable initial data converges to the one of Euler equations, provided that ν = o(α^(4/3)). Moreover, the convergent rate is obtained.

SINGULAR DOUBLE PHASE EQUATIONS

We study a double phase Dirichlet problem with a reaction that has a parametric singular term. Using the Nehari manifold method, we show that for all small values of the parameter, the problem has at least two positive, energy minimizing solutions.

Singularity Formation forthe General Poiseuille Flow of Nematic Liquid Crystals

We consider the Poiseuille flow of nematic liquid crystals via the full Ericksen-Leslie model.The model is described by a coupled system consisting of a heat equation and a quasilinear wave equation.In this paper,we will construct an example with a finite time cusp singularity due to the quasilinearity of the wave equation,extended from an earlier resultonaspecial case.