Computing Numerical Singular Points of Plane Algebraic Curves

Given an irreducible plane algebraic curve of degree d 〉 3, we compute its numerical singular points, determine their multiplicities, and count the number of distinct tangents at each to decide whether the singular points are ordinary. The numerical procedures rely on computing numerical solutions of polynomial systems by homotopy continuation method and a reliable method that calculates multiple roots of the univariate polynomials accurately using standard machine precision. It is completely different from the traditional symbolic computation and provides singular points and their related properties of some plane algebraic curves that the symbolic software Maple cannot work out. Without using multiprecision arithmetic, extensive numerical experiments show that our numerical procedures are accurate, efficient and robust, even if the coefficients of plane algebraic curves are inexact.

NUMERICAL SOLUTION OF QUASILINEAR SINGULARLY PERTURBED ORDINARY DIFFERENTIAL EQUATION WITHOUT TURNING POINTS

In this paper we consider a quasilinear second order ordinary diferential equation with a small parameter Firstly an approximate problem is constructed. Then an iterative procedure is developed. Finally we give an algorithm whose accuracy is good for arbitrary e>0 .

NUMERICAL SOLUTION OF THE SINGULARLY PERTURBED PROBLEM FOR THE HYPERBOLIC EQUATION WITH INITIAL JUMP

In this paper we consider the initial-boundary value problem for a second order hyperbolic equation with initial jump. The bounds on the derivatives of the exact solution are given. Then a difference scheme is constructed on a non-uniform grid. Finally, uniform convergence of the difference solution is proved in the sense of the discrete energy norm.

Pre-stack elastic reverse time migration in tunnels based on cylindrical coordinates

Seismic forward-prospecting in tunnels is an important step to ensure excavation safety. Nowadays, most advanced imaging techniques in seismic exploration involve calculating the solution of elastic wave equation in a certain coordinate system. However, considering the cylindrical geometry of common tunnel body, Cartesian coordinate system seemingly has limited applicability in tunnel seismic forwardprospecting. To accurately simulate the seismic signal received in tunnels, previous imaging method using decoupled non-conversion elastic wave equation is extended from Cartesian coordinates to cylindrical coordinates. The proposed method preserves the general finite-difference time-domain(FDTD)scheme in Cartesian coordinates, except for a novel wavefield calculation strategy addressing the singularity issue inherited at the cylindrical axis. Moreover, the procedure of cylindrical elastic reverse time migration(CERTM) in tunnels is introduced based on the decoupled non-conversion elastic wavefield. Its imaging effect is further validated via numerical experiments on typical tunnel models. As indicated in the synthetic examples, both the PP-and SS-images could clearly show the geological structure in front of the tunnel face without obvious crosstalk artifacts. Migration imaging using PP-waves can present satisfactory results with higher resolution information supplemented by the SS-images. The potential of applying the proposed method in real-world cases is demonstrated in a water diversion tunnel. In the end, we share our insights regarding the singularity treatment and further improvement of the proposed method.

Numerical quadrature for singular and near-singular integrals of boundary element method and its applications in large-scale acoustic problems

The numerical quadrature methods for dealing with the problems of singular and near-singular integrals caused by Burton-Miller method are proposed, by which the conventional and fast multipole BEMs （boundary element methods） for 3D acoustic problems based on constant elements are improved. To solve the problem of singular integrals, a Hadamard finite-part integral method is presented, which is a simplified combination of the methods proposed by Kirkup and Wolf. The problem of near-singular integrals is overcome by the simple method of polar transformation and the more complex method of PART （Projection and Angular ＆ Radial Transformation）. The effectiveness of these methods for solving the singular and near-singular problems is validated through comparing with the results computed by the analytical method and/or the commercial software LMS Virtual.Lab. In addition, the influence of the near-singular integral problem on the computational precisions is analyzed by computing the errors relative to the exact solution. The computational complexities of the conventional and fast multipole BEM are analyzed and compared through numerical computations. A large-scale acoustic scattering problem, whose degree of freedoms is about 340,000, is implemented successfully. The results show that, the near singularity is primarily introduced by the hyper-singular kernel, and has great influences on the precision of the solution. The precision of fast multipole BEM is the same as conventional BEM, but the computational complexities are much lower.

NUMERICAL SOLUTION OF THE SCATTERING PROBLEM FOR ACOUSTIC WAVES BY A TWO-SIDED CRACK IN 2-DIMENSIONAL SPACE

The wave scattering problem by a crack F in R2 with impedance type boundary is considered. This problem models the diffraction of waves by thin two-sided cylindrical screens. A numerical method for solving the problem is developed. The solution of the problem is represented in the form of the combined angular potential and single-layer potential. The linear integral equations satisfied by the density functions are derived for general parameterized arcs. The weakly singular integrals and the Cauchy singular integral arising in these equations are computed using a highly accurate scheme with a truncation error analysis. The advantage of the scheme proposed in this paper is, in one hand, the fact that we do not need the analyticity property of the crack and we allow different complex valued surface impedances in both sides of the crack. In the other hand, we avoid the hyper-singular integrals. Numerical implementations showing the validity of the scheme are presented.

THE APPLICATION OF INTEGRAL EQUATIONS TO THE NUMERICAL SOLUTION OF NONLINEAR SINGULAR PERTURBATION PROBLEMS

The nonlinear singular perturbation problem is solved numerically on nonequidistant meshes which are dense in the boundary layers. The method presented is based on the numerical solution of integral equations [1]. The fourth order uniform accuracy of the scheme is proved. A numerical experiment demonstrates the effectiveness of the method.

Nonlinear wave in granular systems based on elastoplastic dashpot model

The dynamic dashpot models are widely used in EDEM commercial software.However,most dashpot models suffer from a serious numerical issue in calculating the granular chain because the denominator of damping force includes the initial impact velocity.Moreover,the existing dynamic dashpot models extended from the original Hertz contact law overestimated the contact stiffness in the elastoplastic contact phase.These two reasons above result in most dynamic dashpot models confronting some issues in calculating the multiple collision of the granular chain.Therefore,this investigation aims to propose a new composite dynamic dashpot model for the dynamic simulation of granular matters.First,the entire contact process is divided into three different phases:elastic,elastoplastic,and full plastic phases.The Hertz contact stiffness is still used in the elastic contact phase when the contact comes into the elastoplastic or full plastic phase.Hertz contact stiffness in the dynamic dashpot model is replaced by linearizing the contact stiffness from the Ma‐Liu(ML)model in each time step.Second,the whole contact behavior is treated as a linear mass‐spring‐damper model,and the damping factor is obtained by solving the single‐degree‐freedom underdamped vibration equation.The new dynamic dashpot model is proposed by combining the contact stiffnesses in different contact phases and corresponding damping factors,which not only removes the initial im-pact velocity from the denominator of damping force but also updates the contact stiffness based on the constitutive relation of the contact body when the contact comes into the elastoplastic or full plastic phase.Finally,a granular chain is treated as numerical examples to check the reasonability and effectiveness of the new dynamic dashpot model by comparing it to the experimental data.The simulation shows that the solitary waves obtained using the new dashpot model are more accurate than the dashpot model used in EDEM software.