In this paper, for general linear methods applied to strictly dissipative initial value problem in Hilbert spaces, we prove that algebraic stability implies B-convergence, which extends and improves the existing resul...In this paper, for general linear methods applied to strictly dissipative initial value problem in Hilbert spaces, we prove that algebraic stability implies B-convergence, which extends and improves the existing results on Runge-Kutta methods. Specializing our results for the case of multi-step Runge-Kutta methods, a series of B-convergence results are obtained.展开更多
High-order discretizations of partial differential equations(PDEs)necessitate high-order time integration schemes capable of handling both stiff and nonstiff operators in an efficient manner.Implicit-explicit(IMEX)int...High-order discretizations of partial differential equations(PDEs)necessitate high-order time integration schemes capable of handling both stiff and nonstiff operators in an efficient manner.Implicit-explicit(IMEX)integration based on general linear methods(GLMs)offers an attractive solution due to their high stage and method order,as well as excellent stability properties.The IMEX characteristic allows stiff terms to be treated implicitly and nonstiff terms to be efficiently integrated explicitly.This work develops two systematic approaches for the development of IMEX GLMs of arbitrary order with stages that can be solved in parallel.The first approach is based on diagonally implicit multi-stage integration methods(DIMSIMs)of types 3 and 4.The second is a parallel generalization of IMEX Euler and has the interesting feature that the linear stability is independent of the order of accuracy.Numerical experiments confirm the theoretical rates of convergence and reveal that the new schemes are more efficient than serial IMEX GLMs and IMEX Runge-Kutta methods.展开更多
First a general model for a three-step projection method is introduced, and second it has been applied to the approximation solvability of a system of nonlinear variational inequality problems in a Hilbert space setti...First a general model for a three-step projection method is introduced, and second it has been applied to the approximation solvability of a system of nonlinear variational inequality problems in a Hilbert space setting. Let H be a real Hilbert space and K be a nonempty closed convex subset of H. For arbitrarily chosen initial points x0, y0, z0 ∈ K, compute sequences xn, yn, zn such thatT ： K→ H is a nonlinear mapping onto K. At last three-step models are applied to some variational inequality problems.展开更多
The paper presents a new solution of inverse displacement analysis of the general six degree-of-freedom serial robot.The inverse displacement analysis of the general serial robot is transformed into a minimization pro...The paper presents a new solution of inverse displacement analysis of the general six degree-of-freedom serial robot.The inverse displacement analysis of the general serial robot is transformed into a minimization problem and then the optimization method is adopted to solve the nonlinear least squares problem with the analytic form of new Jacobian matrix.In this way,joint variables of the general serial robot can be searched out quickly under the desired precision when positions of the three non-collinear end effector points are given.Compared with the general Newton iterative method,the proposed algorithm can search out the solution when the robot is at the singular configuration and the initial configuration used in the optimization method may also be the singular configuration.So the convergence domain is bigger than that of the general Newton iterative method.Another advantage of the proposed algorithm is that positions of the three non-collinear end effector points are usually much easier to be measured than the orientation of the end effector.The inverse displacement analysis of the general 6R(six-revolute-joint) serial robot is illustrated as an example and the simulation results verify the efficiency of the proposed algorithm.Because the three non-collinear points can be selected at random,the method can be applied to any other types of serial robots.展开更多
In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We the...In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.展开更多
The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary....The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary.Therefore,this method has the merits of purely meshless,easy to program,high solution accuracy and so on.In this paper,we investigate the effectiveness of the BKM for solving Helmholtz-type problems under various conditions through a series of novel numerical experiments.The results demonstrate that the BKM is efficient and achieves high computational accuracy for problems with smooth or continuous boundary conditions.However,when applied to discontinuous boundary problems,the method exhibits significant numerical instability,potentially leading to substantial deviations in the computed results.Finally,three potential improvement strategies are proposed to mitigate this limitation.展开更多
Stair matrices and their generalizations are introduced. The definitions and some properties of the matrices were first given by Lu Hao. This class of matrices provide bases of matrix splittings for iterative methods....Stair matrices and their generalizations are introduced. The definitions and some properties of the matrices were first given by Lu Hao. This class of matrices provide bases of matrix splittings for iterative methods. The remarkable feature of iterative methods based on the new class of matrices is that the methods are easily implemented for parallel computation. In particular, a generalization of the accelerated overrelaxation method （GAOR） is introduced. Some theories of the AOR method are extended to the generalized method to include a wide class of matrices. The convergence of the new method is derived for Hermitian positive definite matrices. Finally, some examples are given in order to show the superiority of the new method.展开更多
Based on the general displacement method and the basic hypothesis of the trial load method, a new advanced trial load method, the general displacement arch-cantilever element method, was proposed to derive the transfo...Based on the general displacement method and the basic hypothesis of the trial load method, a new advanced trial load method, the general displacement arch-cantilever element method, was proposed to derive the transformation relation of displacements and loads between the surface nodes and middle plane nodes. This method considers the nodes on upstream and downstream surfaces of the arch dam to be exit nodes （master nodes）, and the middle plane nodes to be slave nodes. According to the derived displacement and load transformation matrices, the equilibrium equation treating the displacement of middle plane nodes as a basic unknown variable is transformed into one that treats the displacement of upstream and downstream nodes as a basic unknown variable. Because the surface nodes have only three degrees of freedom （DOF）, this method can be directly coupled with the finite element method （FEM）, which is used for foundation simulation to analyze the stress of the arch dam with consideration of dam-foundation interaction. Moreover, using the FEM, the nodal load of the arch dam can be easily obtained. Case studies of a typical cylindrical arch dam and the Wudongde arch dam demonstrate the robustness and feasibility of the proposed method.展开更多
In 1992, Cooper [2] has presented some new stability concepts for Runge-Kutta methods whichis based on two slightly different test problems, and obtained the algebraic conditions that guarantee newstability properties...In 1992, Cooper [2] has presented some new stability concepts for Runge-Kutta methods whichis based on two slightly different test problems, and obtained the algebraic conditions that guarantee newstability properties. In this paper, we extend these results to general linear methods and to more generalproblem class Kστ. The concepts of (k, p, q)-secondary stability and (k, p. q)-secondary stability are introduced, and the criteria of secondary algebraic stability are also established. The criteria relax algebraicstability conditions while retaining the virtues of a nonlinear test problem.展开更多
According to characteristics of General Entomology and existing problems in teaching process,it came up with methods and countermeasures for improving teaching General Entomology,including improving practical ability ...According to characteristics of General Entomology and existing problems in teaching process,it came up with methods and countermeasures for improving teaching General Entomology,including improving practical ability of students through enhancing all links of practice teaching,and stimulating learning interest of students through improving methods of course examination.展开更多
On the basis of assuming that the narrow state X（3872） is a molecule state consisting of D0 and D＊0, we apply the Mandelstam generalization of the Ge11-Mann-Low method to calculate the matrix element of quark curre...On the basis of assuming that the narrow state X（3872） is a molecule state consisting of D0 and D＊0, we apply the Mandelstam generalization of the Ge11-Mann-Low method to calculate the matrix element of quark current between the heavy meson states described by Bether-Salpeter wave function. In calculation of the matrix element of quark current the operator product expansion is used in order to include the nonperturbative contribution of the vacuum condensates. In this scheme we calculate the mass of X（3872）. We believe that this scheme is closer to QCD than the previous work.展开更多
In this article, we propose a generalized exp(-Φ(ξ))-expansion method and successfully implement it to find exact traveling wave solutions to the fifth order standard Sawada-Kotera (SK) equation. The exact traveling...In this article, we propose a generalized exp(-Φ(ξ))-expansion method and successfully implement it to find exact traveling wave solutions to the fifth order standard Sawada-Kotera (SK) equation. The exact traveling wave solutions are established in the form of trigonometric, hyperbolic, exponential and rational functions with some free parameters. It is shown that this method is standard, effective and easily applicable mathematical tool for solving nonlinear partial differential equations arises in the field of mathematical physics and engineering.展开更多
The article proposes an Equivalent Single Layer(ESL)formulation for the linear static analysis of arbitrarily-shaped shell structures subjected to general surface loads and boundary conditions.A parametrization of the...The article proposes an Equivalent Single Layer(ESL)formulation for the linear static analysis of arbitrarily-shaped shell structures subjected to general surface loads and boundary conditions.A parametrization of the physical domain is provided by employing a set of curvilinear principal coordinates.The generalized blendingmethodology accounts for a distortion of the structure so that disparate geometries can be considered.Each layer of the stacking sequence has an arbitrary orientation and is modelled as a generally anisotropic continuum.In addition,re-entrant auxetic three-dimensional honeycomb cells with soft-core behaviour are considered in the model.The unknown variables are described employing a generalized displacement field and pre-determined through-the-thickness functions assessed in a unified formulation.Then,a weak assessment of the structural problem accounts for shape functions defined with an isogeometric approach starting fromthe computational grid.Ageneralizedmethodology has been proposed to define two-dimensional distributions of static surface loads.In the same way,boundary conditions with three-dimensional features are implemented along the shell edges employing linear springs.The fundamental relations are obtained from the stationary configuration of the total potential energy,and they are numerically tackled by employing the Generalized Differential Quadrature(GDQ)method,accounting for nonuniform computational grids.In the post-processing stage,an equilibrium-based recovery procedure allows the determination of the three-dimensional dispersion of the kinematic and static quantities.Some case studies have been presented,and a successful benchmark of different structural responses has been performed with respect to various refined theories.展开更多
Generalized Jacobi polynomials with indexes α,β∈ R are introduced and some basic properties are established. As examples of applications,the second- and fourth-order elliptic boundary value problems with Dirichlet ...Generalized Jacobi polynomials with indexes α,β∈ R are introduced and some basic properties are established. As examples of applications,the second- and fourth-order elliptic boundary value problems with Dirichlet or Robin boundary conditions are considered,and the generalized Jacobi spectral schemes are proposed. For the diagonalization of discrete systems,the Jacobi-Sobolev orthogonal basis functions are constructed,which allow the exact solutions and the approximate solutions to be represented in the forms of infinite and truncated Jacobi series. Error estimates are obtained and numerical results are provided to illustrate the effectiveness and the spectral accuracy.展开更多
The boundary knot method(BKM) is a boundary-type meshfree method. Only non-singular general solutions are used during the whole solution procedures. The effective condition number(ECN), which depends on the right-hand...The boundary knot method(BKM) is a boundary-type meshfree method. Only non-singular general solutions are used during the whole solution procedures. The effective condition number(ECN), which depends on the right-hand side vector of a linear system, is considered as an alternative criterion to the traditional condition number. In this paper, the effective condition number is used to help determine the position and distribution of the collocation points as well as the quasi-optimal collocation point numbers. During the solution process, we propose an NMN-search algorithm. Numerical examples show that the ECN is reliable to measure the feasibility of the BKM.展开更多
By asing the nonclassical method of symmetry reductions, the exact solutions for general variable coefficient KdV equation with dissipative loss and nonuniformity terms are obtained. When the dissipative loss and non...By asing the nonclassical method of symmetry reductions, the exact solutions for general variable coefficient KdV equation with dissipative loss and nonuniformity terms are obtained. When the dissipative loss and nonuniformity terms don't exist, the multisoliton solutions are found and the corresponding Painleve II type equation for the variable coefficient KdV equation is given.展开更多
The bilinear equation of the genera/nonlinear Schrodinger equation with derivative （GDNLSE） and the N-soliton solutions are obtained through the dependent variable transformation and the Hirota method, respectively....The bilinear equation of the genera/nonlinear Schrodinger equation with derivative （GDNLSE） and the N-soliton solutions are obtained through the dependent variable transformation and the Hirota method, respectively. The bilinear equation of the nonlinear Schrodinger equation with derivative （DNLSE） and its multisoliton solutions are given by reduction.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
A hyperbolic Lindstedt-Poincare method is presented to determine the homoclinic solutions of a kind of nonlinear oscillators, in which critical value of the homoclinic bifurcation parameter can be determined. The gene...A hyperbolic Lindstedt-Poincare method is presented to determine the homoclinic solutions of a kind of nonlinear oscillators, in which critical value of the homoclinic bifurcation parameter can be determined. The generalized Lienard oscillator is studied in detail, and the present method＇s predictions are compared with those of Runge-Kutta method to illustrate its accuracy.展开更多
With the development of large-scale spectral surveys, fiber positioning technology has been developing rapidly. Because of the performance advantages of a four-quadrant（4Q） detector, a fiber positioning and real-tim...With the development of large-scale spectral surveys, fiber positioning technology has been developing rapidly. Because of the performance advantages of a four-quadrant（4Q） detector, a fiber positioning and real-time monitoring system based on the 4Q detector is proposed. The detection accuracy of this system is directly determined by the precision of the center of the spot. A Gaussian fitting algorithm based on the 4Q detector is studied and applied in the fiber positioning process to improve the calculated accuracy of the spot center. The relationship between the center position of the incident spot and the detector output signal is deduced. An experimental platform is built to complete the simulated experiment. Then we use the Gaussian fitting method to process experimental data, compare the fitting value with the theoretical one and calculate the corresponding error.展开更多
基金The project supported by the National Natural Science Foundation of China
文摘In this paper, for general linear methods applied to strictly dissipative initial value problem in Hilbert spaces, we prove that algebraic stability implies B-convergence, which extends and improves the existing results on Runge-Kutta methods. Specializing our results for the case of multi-step Runge-Kutta methods, a series of B-convergence results are obtained.
基金funded by awards NSF CCF1613905,NSF ACI1709727,AFOSR DDDAS FA9550-17-1-0015the Computational Science Laboratory at Virginia Tech.
文摘High-order discretizations of partial differential equations(PDEs)necessitate high-order time integration schemes capable of handling both stiff and nonstiff operators in an efficient manner.Implicit-explicit(IMEX)integration based on general linear methods(GLMs)offers an attractive solution due to their high stage and method order,as well as excellent stability properties.The IMEX characteristic allows stiff terms to be treated implicitly and nonstiff terms to be efficiently integrated explicitly.This work develops two systematic approaches for the development of IMEX GLMs of arbitrary order with stages that can be solved in parallel.The first approach is based on diagonally implicit multi-stage integration methods(DIMSIMs)of types 3 and 4.The second is a parallel generalization of IMEX Euler and has the interesting feature that the linear stability is independent of the order of accuracy.Numerical experiments confirm the theoretical rates of convergence and reveal that the new schemes are more efficient than serial IMEX GLMs and IMEX Runge-Kutta methods.
文摘First a general model for a three-step projection method is introduced, and second it has been applied to the approximation solvability of a system of nonlinear variational inequality problems in a Hilbert space setting. Let H be a real Hilbert space and K be a nonempty closed convex subset of H. For arbitrarily chosen initial points x0, y0, z0 ∈ K, compute sequences xn, yn, zn such thatT ： K→ H is a nonlinear mapping onto K. At last three-step models are applied to some variational inequality problems.
基金Funded by National Natural Science Foundation of China (No. 50905102)the Natural Science Foundation of Guangdong Province (Nos. 10151503101000033 and 8351503101000001)the Building Fund for the Academic Innovation Team of Shantou University (No. ITC10003)
文摘The paper presents a new solution of inverse displacement analysis of the general six degree-of-freedom serial robot.The inverse displacement analysis of the general serial robot is transformed into a minimization problem and then the optimization method is adopted to solve the nonlinear least squares problem with the analytic form of new Jacobian matrix.In this way,joint variables of the general serial robot can be searched out quickly under the desired precision when positions of the three non-collinear end effector points are given.Compared with the general Newton iterative method,the proposed algorithm can search out the solution when the robot is at the singular configuration and the initial configuration used in the optimization method may also be the singular configuration.So the convergence domain is bigger than that of the general Newton iterative method.Another advantage of the proposed algorithm is that positions of the three non-collinear end effector points are usually much easier to be measured than the orientation of the end effector.The inverse displacement analysis of the general 6R(six-revolute-joint) serial robot is illustrated as an example and the simulation results verify the efficiency of the proposed algorithm.Because the three non-collinear points can be selected at random,the method can be applied to any other types of serial robots.
基金The project supported by the China NKBRSF(2001CB409604)
文摘In this paper, a fully third-order accurate projection method for solving the incompressible Navier-Stokes equations is proposed. To construct the scheme, a continuous projection procedure is firstly presented. We then derive a sufficient condition for the continuous projection equations to be temporally third-order accurate approximations of the original Navier-Stokes equations by means of the localtruncation-error-analysis technique. The continuous projection equations are discretized temporally and spatially to third-order accuracy on the staggered grids, resulting in a fully third-order discrete projection scheme. The possibility to design higher-order projection methods is thus demonstrated in the present paper. A heuristic stability analysis is performed on this projection method showing the probability of its being stable. The stability of the present scheme is further verified through numerical tests. The third-order accuracy of the present projection method is validated by several numerical test cases.
基金Supported by the Key Scientific Research Plan of Colleges and Universities in Henan Province(23B140006)。
文摘The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary.Therefore,this method has the merits of purely meshless,easy to program,high solution accuracy and so on.In this paper,we investigate the effectiveness of the BKM for solving Helmholtz-type problems under various conditions through a series of novel numerical experiments.The results demonstrate that the BKM is efficient and achieves high computational accuracy for problems with smooth or continuous boundary conditions.However,when applied to discontinuous boundary problems,the method exhibits significant numerical instability,potentially leading to substantial deviations in the computed results.Finally,three potential improvement strategies are proposed to mitigate this limitation.
基金Project supported by the Natural Science Foundation of Liaoning Province of China (No.20022021)
文摘Stair matrices and their generalizations are introduced. The definitions and some properties of the matrices were first given by Lu Hao. This class of matrices provide bases of matrix splittings for iterative methods. The remarkable feature of iterative methods based on the new class of matrices is that the methods are easily implemented for parallel computation. In particular, a generalization of the accelerated overrelaxation method （GAOR） is introduced. Some theories of the AOR method are extended to the generalized method to include a wide class of matrices. The convergence of the new method is derived for Hermitian positive definite matrices. Finally, some examples are given in order to show the superiority of the new method.
基金supported by the National Natural Science Foundation of China (Grant No. 90510017)
文摘Based on the general displacement method and the basic hypothesis of the trial load method, a new advanced trial load method, the general displacement arch-cantilever element method, was proposed to derive the transformation relation of displacements and loads between the surface nodes and middle plane nodes. This method considers the nodes on upstream and downstream surfaces of the arch dam to be exit nodes （master nodes）, and the middle plane nodes to be slave nodes. According to the derived displacement and load transformation matrices, the equilibrium equation treating the displacement of middle plane nodes as a basic unknown variable is transformed into one that treats the displacement of upstream and downstream nodes as a basic unknown variable. Because the surface nodes have only three degrees of freedom （DOF）, this method can be directly coupled with the finite element method （FEM）, which is used for foundation simulation to analyze the stress of the arch dam with consideration of dam-foundation interaction. Moreover, using the FEM, the nodal load of the arch dam can be easily obtained. Case studies of a typical cylindrical arch dam and the Wudongde arch dam demonstrate the robustness and feasibility of the proposed method.
文摘In 1992, Cooper [2] has presented some new stability concepts for Runge-Kutta methods whichis based on two slightly different test problems, and obtained the algebraic conditions that guarantee newstability properties. In this paper, we extend these results to general linear methods and to more generalproblem class Kστ. The concepts of (k, p, q)-secondary stability and (k, p. q)-secondary stability are introduced, and the criteria of secondary algebraic stability are also established. The criteria relax algebraicstability conditions while retaining the virtues of a nonlinear test problem.
文摘According to characteristics of General Entomology and existing problems in teaching process,it came up with methods and countermeasures for improving teaching General Entomology,including improving practical ability of students through enhancing all links of practice teaching,and stimulating learning interest of students through improving methods of course examination.
基金Supported in part by the National Natural Science Foundation of China under Grant No. 10335012 and the National Key Basic Research Program and Cross Science of China under Grant No. 90503011
文摘On the basis of assuming that the narrow state X（3872） is a molecule state consisting of D0 and D＊0, we apply the Mandelstam generalization of the Ge11-Mann-Low method to calculate the matrix element of quark current between the heavy meson states described by Bether-Salpeter wave function. In calculation of the matrix element of quark current the operator product expansion is used in order to include the nonperturbative contribution of the vacuum condensates. In this scheme we calculate the mass of X（3872）. We believe that this scheme is closer to QCD than the previous work.
文摘In this article, we propose a generalized exp(-Φ(ξ))-expansion method and successfully implement it to find exact traveling wave solutions to the fifth order standard Sawada-Kotera (SK) equation. The exact traveling wave solutions are established in the form of trigonometric, hyperbolic, exponential and rational functions with some free parameters. It is shown that this method is standard, effective and easily applicable mathematical tool for solving nonlinear partial differential equations arises in the field of mathematical physics and engineering.
文摘The article proposes an Equivalent Single Layer(ESL)formulation for the linear static analysis of arbitrarily-shaped shell structures subjected to general surface loads and boundary conditions.A parametrization of the physical domain is provided by employing a set of curvilinear principal coordinates.The generalized blendingmethodology accounts for a distortion of the structure so that disparate geometries can be considered.Each layer of the stacking sequence has an arbitrary orientation and is modelled as a generally anisotropic continuum.In addition,re-entrant auxetic three-dimensional honeycomb cells with soft-core behaviour are considered in the model.The unknown variables are described employing a generalized displacement field and pre-determined through-the-thickness functions assessed in a unified formulation.Then,a weak assessment of the structural problem accounts for shape functions defined with an isogeometric approach starting fromthe computational grid.Ageneralizedmethodology has been proposed to define two-dimensional distributions of static surface loads.In the same way,boundary conditions with three-dimensional features are implemented along the shell edges employing linear springs.The fundamental relations are obtained from the stationary configuration of the total potential energy,and they are numerically tackled by employing the Generalized Differential Quadrature(GDQ)method,accounting for nonuniform computational grids.In the post-processing stage,an equilibrium-based recovery procedure allows the determination of the three-dimensional dispersion of the kinematic and static quantities.Some case studies have been presented,and a successful benchmark of different structural responses has been performed with respect to various refined theories.
基金the National Natural Science Foundation of China (Nos.11571238,11601332,91130014,11471312 and 91430216).
文摘Generalized Jacobi polynomials with indexes α,β∈ R are introduced and some basic properties are established. As examples of applications,the second- and fourth-order elliptic boundary value problems with Dirichlet or Robin boundary conditions are considered,and the generalized Jacobi spectral schemes are proposed. For the diagonalization of discrete systems,the Jacobi-Sobolev orthogonal basis functions are constructed,which allow the exact solutions and the approximate solutions to be represented in the forms of infinite and truncated Jacobi series. Error estimates are obtained and numerical results are provided to illustrate the effectiveness and the spectral accuracy.
基金Supported by the Natural Science Foundation of Anhui Province(1908085QA09)Higher Education Department of the Ministry of Education(201802358008)
文摘The boundary knot method(BKM) is a boundary-type meshfree method. Only non-singular general solutions are used during the whole solution procedures. The effective condition number(ECN), which depends on the right-hand side vector of a linear system, is considered as an alternative criterion to the traditional condition number. In this paper, the effective condition number is used to help determine the position and distribution of the collocation points as well as the quasi-optimal collocation point numbers. During the solution process, we propose an NMN-search algorithm. Numerical examples show that the ECN is reliable to measure the feasibility of the BKM.
基金Supported by the Develop Programme Foundation of the National Basic research(G1 9990 3 2 80 1 )
文摘By asing the nonclassical method of symmetry reductions, the exact solutions for general variable coefficient KdV equation with dissipative loss and nonuniformity terms are obtained. When the dissipative loss and nonuniformity terms don't exist, the multisoliton solutions are found and the corresponding Painleve II type equation for the variable coefficient KdV equation is given.
基金The project supported by National Natural Science Foundation of China under Grant No.10671121
文摘The bilinear equation of the genera/nonlinear Schrodinger equation with derivative （GDNLSE） and the N-soliton solutions are obtained through the dependent variable transformation and the Hirota method, respectively. The bilinear equation of the nonlinear Schrodinger equation with derivative （DNLSE） and its multisoliton solutions are given by reduction.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金supported by the National Natural Science Foundation of China (10672193)Sun Yat-sen University (Fu Lan Scholarship)the University of Hong Kong (CRGC grant).
文摘A hyperbolic Lindstedt-Poincare method is presented to determine the homoclinic solutions of a kind of nonlinear oscillators, in which critical value of the homoclinic bifurcation parameter can be determined. The generalized Lienard oscillator is studied in detail, and the present method＇s predictions are compared with those of Runge-Kutta method to illustrate its accuracy.
基金support by the Fundamental Research Funds for the Central Universities of China (2013/B15020271)the National Natural Science Foundation of China (1014/515029111)the National Undergraduate Training Program for Innovation and Entrepreneurship (201610294069)
文摘With the development of large-scale spectral surveys, fiber positioning technology has been developing rapidly. Because of the performance advantages of a four-quadrant（4Q） detector, a fiber positioning and real-time monitoring system based on the 4Q detector is proposed. The detection accuracy of this system is directly determined by the precision of the center of the spot. A Gaussian fitting algorithm based on the 4Q detector is studied and applied in the fiber positioning process to improve the calculated accuracy of the spot center. The relationship between the center position of the incident spot and the detector output signal is deduced. An experimental platform is built to complete the simulated experiment. Then we use the Gaussian fitting method to process experimental data, compare the fitting value with the theoretical one and calculate the corresponding error.
