Local mesh refinement is one of the key steps in the implementations of adaptive finite element methods. This paper presents a parallel algorithm for distributed memory parallel computers for adaptive local refinement...Local mesh refinement is one of the key steps in the implementations of adaptive finite element methods. This paper presents a parallel algorithm for distributed memory parallel computers for adaptive local refinement of tetrahedral meshes using bisection. This algorithm is used in PHG, Parallel Hierarchical Grid Chttp：//lsec. cc. ac. cn/phg/）, a toolbox under active development for parallel adaptive finite element solutions of partial differential equations. The algorithm proposed is characterized by allowing simukaneous refinement of submeshes to arbitrary levels before synchronization between submeshes and without the need of a central coordinator process for managing new vertices. Using the concept of canonical refinement, a simple proof of the independence of the resulting mesh on the mesh partitioning is given, which is useful in better understanding the behaviour of the biseetioning refinement procedure.展开更多
For the large sparse block two-by-two real nonsingular matrices, we establish a general framework of structured preconditioners through matrix transformation and matrix approximations. For the specific versions such a...For the large sparse block two-by-two real nonsingular matrices, we establish a general framework of structured preconditioners through matrix transformation and matrix approximations. For the specific versions such as modified block Jacobi-type, modified block Gauss-Seidel-type, and modified block unsymmetric (symmetric) Gauss-Seidel-type preconditioners, we precisely describe their concrete expressions and deliberately analyze eigenvalue distributions and positive definiteness of the preconditioned matrices. Also, we show that when these structured preconditioners are employed to precondition the Krylov subspace methods such as GMRES and restarted GMRES, fast and effective iteration solvers can be obtained for the large sparse systems of linear equations with block two-by-two coefficient matrices. In particular, these structured preconditioners can lead to high-quality preconditioning matrices for some typical matrices from the real-world applications.展开更多
The algebraic structure and Poisson＇s integral theory of mechanico-electrical systems are studied. The Hamilton canonical equations and generalized Hamilton canonical equations and their the contravariant algebraic f...The algebraic structure and Poisson＇s integral theory of mechanico-electrical systems are studied. The Hamilton canonical equations and generalized Hamilton canonical equations and their the contravariant algebraic forms for mechanico-electrical systems are obtained. The Lie algebraic structure and the Poisson＇s integral theory of Lagrange mechanico-electrical systems are derived. The Lie algebraic structure admitted and Poisson＇s integral theory of the Lagrange-Maxwell mechanico-electrical systems are presented. Two examples are presented to illustrate these results.展开更多
In this paper Mei symmetry is introduced for a nonconservative system. The necessary and sufficient condition for a Mei symmetry to be also a Lie symmetry is derived. It is proved that the Mei symmetry leads to a non-...In this paper Mei symmetry is introduced for a nonconservative system. The necessary and sufficient condition for a Mei symmetry to be also a Lie symmetry is derived. It is proved that the Mei symmetry leads to a non-Noether conservative quantity via a Lie symmetry, and deduces a Lutzky conservative quantity via a Lie point symmetry.展开更多
This paper presents a method to find Noether-type conserved quantities and Lie point symmetries for discrete mechanico-electrical dynamical systems,which leave invuriant the set of solutions of the corresponding diffe...This paper presents a method to find Noether-type conserved quantities and Lie point symmetries for discrete mechanico-electrical dynamical systems,which leave invuriant the set of solutions of the corresponding difference scheme. This approach makes it possible to devise techniques for solving the Lagrange Maxwell equations in differences which correspond to mechanico-electrical systems,by adapting existing differential equations.In particular,it obtains a new systematic method to determine both the one-parameter Lie groups and the discrete Noether conserved quantities of Lie point symmetries for mechanico-electrical systems.As an application,it obtains the Lie point symmetries and the conserved quantities for the difference equation of a model that represents a capacitor microphone.展开更多
A discrete total variation calculus with variable time steps is presented for mechanico-electrical systems where there exist non-potential and dissipative forces. By using this discrete variation calculus, the symplec...A discrete total variation calculus with variable time steps is presented for mechanico-electrical systems where there exist non-potential and dissipative forces. By using this discrete variation calculus, the symplectic-energy-first integrators for mechanico-electrical systems are derived. To do this, the time step adaptation is employed. The discrete variational principle and the Euler-Lagrange equation are derived for the systems. By using this discrete algorithm it is shown that mechanico-electrical systems are not symplectic and their energies are not conserved unless they are Lagrange mechanico-electrical systems. A practical example is presented to illustrate these results.展开更多
For an upper bound of the spectral radius of the QHSS (quasi Hermitian and skew-Hermitian splitting) iteration matrix which can also bound the contraction factor of the QHSS iteration method,we give its minimum point ...For an upper bound of the spectral radius of the QHSS (quasi Hermitian and skew-Hermitian splitting) iteration matrix which can also bound the contraction factor of the QHSS iteration method,we give its minimum point under the conditions which guarantee that the upper bound is strictly less than one. This provides a good choice of the involved iteration parameters,so that the convergence rate of the QHSS iteration method can be significantly improved.展开更多
With the Gravity Recovery and Climate Experiment {GRACE） mission as the prime example, an overview is given on the management and processing of Level IA data of a low-low satellite to satellite tracking mission. To i...With the Gravity Recovery and Climate Experiment {GRACE） mission as the prime example, an overview is given on the management and processing of Level IA data of a low-low satellite to satellite tracking mission. To illustrate the underlying principle and algorithm, a detailed study is made on the K-band ranging （KBR） assembly, which includes the measurement principles, modeling of noises, the generation of Level 1A data from that of Level 0 as well as Level IA to Level IB data processing.展开更多
A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illus...A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illustrate the accuracy and feasibility of this method.展开更多
Non-Noether symmetries and conservative quantities of nonholonomic nonconservative dynamical systems are investigated in this paper. Based on the relationships among motion, nonconservative forces, nonholonomic constr...Non-Noether symmetries and conservative quantities of nonholonomic nonconservative dynamical systems are investigated in this paper. Based on the relationships among motion, nonconservative forces, nonholonomic constrained forces and Lagrangian, non-Noether symmetries and Lutzky conservative quantities are presented for nonholonomic nonconservative dynamical systems. The relation between non-Noether symmetry and Noether symmetry is discussed and it is further shown that non-Noether conservative quantities can be obtained by a complete set of Noether invariants. Finally, an example is given to illustrate these results.展开更多
A set of generalized symmetries with arbitrary functions of t for the Konopelchenko-Dubrovsky (KD)equation in 2+1 space dimensions is given by using a direct method called formal function series method presented by Lo...A set of generalized symmetries with arbitrary functions of t for the Konopelchenko-Dubrovsky (KD)equation in 2+1 space dimensions is given by using a direct method called formal function series method presented by Lou. These symmetries constitute an infinite-dimensional generalized w∞ algebra.展开更多
A generalized AKNS isospectral problem where the trace of corresponding spectral matrix is not zero, is transformed to a new isospectral problem where the trace of the resulting matrix is zero, by using transformation...A generalized AKNS isospectral problem where the trace of corresponding spectral matrix is not zero, is transformed to a new isospectral problem where the trace of the resulting matrix is zero, by using transformation of Lax pairs, and these two spectral problems lead to the same hierarchy of equations. The authors started from the transformed spectral problem and constructed a new loop algebra which has not appeared before, and obtained the integrable coupling of the generalized AKNS hierarchy. Specially, the integrable couplings of the KdV equation and MKdV equation are obtained.展开更多
In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic ...In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.展开更多
This paper presents a discrete vaxiational principle and a method to build first-integrals for finite dimensional Lagrange-Maxwell mechanico-electrical systems with nonconservative forces and a dissipation function. T...This paper presents a discrete vaxiational principle and a method to build first-integrals for finite dimensional Lagrange-Maxwell mechanico-electrical systems with nonconservative forces and a dissipation function. The discrete variational principle and the corresponding Euler-Lagrange equations are derived from a discrete action associated to these systems. The first-integrals are obtained by introducing the infinitesimal transformation with respect to the generalized coordinates and electric quantities of the systems. This work also extends discrete Noether symmetries to mechanico-electrical dynamical systems. A practical example is presented to illustrate the results.展开更多
Based on splitting multi-symplectic structures, a new multi-symplectic scheme is proposed and applied to a nonlinear wave equation. The explicit multi-symplectic scheme of the nonlinear wave equation is obtained, and ...Based on splitting multi-symplectic structures, a new multi-symplectic scheme is proposed and applied to a nonlinear wave equation. The explicit multi-symplectic scheme of the nonlinear wave equation is obtained, and the corresponding multi-symplectic conservation property is proved. The backward error analysis shows that the explicit multi-symplectic scheme has good accuracy. The sine-Gordon equation and the Klein-Gordon equation are simulated by an explicit multi-symplectic scheme. The numerical results show that the new explicit multi-symplectic scheme can well simulate the solitary wave behaviors of the nonlinear wave equation and approximately preserve the relative energy error of the equation.展开更多
Stochastic approximation problem is to find some root or extremum of a nonlinear function for which only noisy measurements of the function are available. The classical algorithm for stochastic approximation problem i...Stochastic approximation problem is to find some root or extremum of a nonlinear function for which only noisy measurements of the function are available. The classical algorithm for stochastic approximation problem is the Robbins-Monro （RM） algorithm, which uses the noisy evaluation of the negative gradient direction as the iterative direction. In order to accelerate the RM algorithm, this paper gives a flame algorithm using adaptive iterative directions. At each iteration, the new algorithm goes towards either the noisy evaluation of the negative gradient direction or some other directions under some switch criterions. Two feasible choices of the criterions are proposed and two corresponding frame algorithms are formed. Different choices of the directions under the same given switch criterion in the frame can also form different algorithms. We also proposed the simultanous perturbation difference forms for the two frame algorithms. The almost surely convergence of the new algorithms are all established. The numerical experiments show that the new algorithms are promising.展开更多
The （3＋1）-dimensional Jimbo-Miwa （JM） equation is solved approximately by using the conformal invariant asymptotic expansion approach presented by Ruan. By solving the new （3＋1）-dimensional integrable models, ...The （3＋1）-dimensional Jimbo-Miwa （JM） equation is solved approximately by using the conformal invariant asymptotic expansion approach presented by Ruan. By solving the new （3＋1）-dimensional integrable models, which are conformal invariant and possess Painlevé property, the approximate solutions are obtained for the JM equation, containing not only one-soliton solutions but also periodic solutions and multi-soliton solutions. Some approximate solutions happen to be exact and some approximate solutions can become exact by choosing relations between the parameters properly.展开更多
Starting from the variable separation approach, the algebraic soliton solution and the solution describing the interaction between line soliton and algebraic soliton are obtained by selecting appropriate seed solution...Starting from the variable separation approach, the algebraic soliton solution and the solution describing the interaction between line soliton and algebraic soliton are obtained by selecting appropriate seed solution for （2＋1）-dimensional ANNV equation. The behaviors of interactions are discussed in detail both analytically and graphically. It is shown that there are two kinds of singular interactions between line soliton and algebraic soliton： 1） the resonant interaction where the algebraic soliton propagates together with the line soliton and persists infinitely; 2） the extremely repulsive interaction where the algebraic soliton affects the motion of the line soliton infinitely apart.展开更多
In this paper,we develop a new sixth-order WENO scheme by adopting a convex combina-tion of a sixth-order global reconstruction and four low-order local reconstructions.Unlike the classical WENO schemes,the associated...In this paper,we develop a new sixth-order WENO scheme by adopting a convex combina-tion of a sixth-order global reconstruction and four low-order local reconstructions.Unlike the classical WENO schemes,the associated linear weights of the new scheme can be any positive numbers with the only requirement that their sum equals one.Further,a very simple smoothness indicator for the global stencil is proposed.The new scheme can achieve sixth-order accuracy in smooth regions.Numerical tests in some one-and two-dimensional bench-mark problems show that the new scheme has a little bit higher resolution compared with the recently developed sixth-order WENO-Z6 scheme,and it is more efficient than the classical fifth-order WENO-JS5 scheme and the recently developed sixth-order WENO6-S scheme.展开更多
In this section, we introduce some outstanding research achievements made by CAS scientists in 2010, as listed in the Chinese version of Science Development Report published in March 2011.
基金supported by the 973 Program of China 2005CB321702China NSF 10531080.
文摘Local mesh refinement is one of the key steps in the implementations of adaptive finite element methods. This paper presents a parallel algorithm for distributed memory parallel computers for adaptive local refinement of tetrahedral meshes using bisection. This algorithm is used in PHG, Parallel Hierarchical Grid Chttp：//lsec. cc. ac. cn/phg/）, a toolbox under active development for parallel adaptive finite element solutions of partial differential equations. The algorithm proposed is characterized by allowing simukaneous refinement of submeshes to arbitrary levels before synchronization between submeshes and without the need of a central coordinator process for managing new vertices. Using the concept of canonical refinement, a simple proof of the independence of the resulting mesh on the mesh partitioning is given, which is useful in better understanding the behaviour of the biseetioning refinement procedure.
文摘For the large sparse block two-by-two real nonsingular matrices, we establish a general framework of structured preconditioners through matrix transformation and matrix approximations. For the specific versions such as modified block Jacobi-type, modified block Gauss-Seidel-type, and modified block unsymmetric (symmetric) Gauss-Seidel-type preconditioners, we precisely describe their concrete expressions and deliberately analyze eigenvalue distributions and positive definiteness of the preconditioned matrices. Also, we show that when these structured preconditioners are employed to precondition the Krylov subspace methods such as GMRES and restarted GMRES, fast and effective iteration solvers can be obtained for the large sparse systems of linear equations with block two-by-two coefficient matrices. In particular, these structured preconditioners can lead to high-quality preconditioning matrices for some typical matrices from the real-world applications.
基金Project supported by the State Key Laboratory of Scientific and Engineering Computing, Chinese Academy of Sciences and the National Natural Science Foundation of China （Grant Nos 10471145 and 10372053） and the Natural Science Foundation of Henan Provincial Government of China （Grant Nos 0311011400 and 0511022200）.
文摘The algebraic structure and Poisson＇s integral theory of mechanico-electrical systems are studied. The Hamilton canonical equations and generalized Hamilton canonical equations and their the contravariant algebraic forms for mechanico-electrical systems are obtained. The Lie algebraic structure and the Poisson＇s integral theory of Lagrange mechanico-electrical systems are derived. The Lie algebraic structure admitted and Poisson＇s integral theory of the Lagrange-Maxwell mechanico-electrical systems are presented. Two examples are presented to illustrate these results.
基金Project supported by the State Key Laboratory of Scientific and Engineering Computing, Chinese Academy of Sciences and the National Natural Science Foundation of China （Grant Nos 10672143, 10471145 and 10372053） and the Natural Science Foundation of Henan Province Government of China（Grant Nos 0511022200 and 0311011400）.
文摘In this paper Mei symmetry is introduced for a nonconservative system. The necessary and sufficient condition for a Mei symmetry to be also a Lie symmetry is derived. It is proved that the Mei symmetry leads to a non-Noether conservative quantity via a Lie symmetry, and deduces a Lutzky conservative quantity via a Lie point symmetry.
基金Project supported by the National Natural Science Foundation of China (Grants Nos 10672143 and 60575055)State Key Laboratory of Scientific and Engineering Computing,Chinese Academy of Sciences+1 种基金Tang Yi-Fa acknowledges the support under Sabbatical Program (SAB2006-0070) of the Spanish Ministry of Education and ScienceJimnez S and Vzquez L acknowledge support of the Spanish Ministry of Education and Science (Grant No MTM2005-05573)
文摘This paper presents a method to find Noether-type conserved quantities and Lie point symmetries for discrete mechanico-electrical dynamical systems,which leave invuriant the set of solutions of the corresponding difference scheme. This approach makes it possible to devise techniques for solving the Lagrange Maxwell equations in differences which correspond to mechanico-electrical systems,by adapting existing differential equations.In particular,it obtains a new systematic method to determine both the one-parameter Lie groups and the discrete Noether conserved quantities of Lie point symmetries for mechanico-electrical systems.As an application,it obtains the Lie point symmetries and the conserved quantities for the difference equation of a model that represents a capacitor microphone.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 10672143 and 60575055)the State Key Laboratory of Scientific and Engineering Computing, Chinese Academy of Sciencesthe Natural Science Foundation of Henan Province Government, China (Grant No 0511022200)
文摘A discrete total variation calculus with variable time steps is presented for mechanico-electrical systems where there exist non-potential and dissipative forces. By using this discrete variation calculus, the symplectic-energy-first integrators for mechanico-electrical systems are derived. To do this, the time step adaptation is employed. The discrete variational principle and the Euler-Lagrange equation are derived for the systems. By using this discrete algorithm it is shown that mechanico-electrical systems are not symplectic and their energies are not conserved unless they are Lagrange mechanico-electrical systems. A practical example is presented to illustrate these results.
基金the National Natural Science Foundation (No.11671393),China.
文摘For an upper bound of the spectral radius of the QHSS (quasi Hermitian and skew-Hermitian splitting) iteration matrix which can also bound the contraction factor of the QHSS iteration method,we give its minimum point under the conditions which guarantee that the upper bound is strictly less than one. This provides a good choice of the involved iteration parameters,so that the convergence rate of the QHSS iteration method can be significantly improved.
基金the project entitled"Advanced Gravity Measurement in Space"supported by the National Space Science Center,Chinese Academy of Sciences Profs.Wenrui Hu and Houze Xu's effort to promote satellite gravity research in China motivated the feasibility study in the first placeSupport from National Natural Science Foundation of China(11305255,11171329 and 41404019)funding from State Key Laboratory of Geodesy and Earth's Dynamics,Institute of Geodesy and Geophysics,Chinese Academy of Sciences(SKLGED2013-3-8-E)are acknowledged
文摘With the Gravity Recovery and Climate Experiment {GRACE） mission as the prime example, an overview is given on the management and processing of Level IA data of a low-low satellite to satellite tracking mission. To illustrate the underlying principle and algorithm, a detailed study is made on the K-band ranging （KBR） assembly, which includes the measurement principles, modeling of noises, the generation of Level 1A data from that of Level 0 as well as Level IA to Level IB data processing.
文摘A posteriori error estimate of the discontinuous-streamline diffusion method for first-order hyperbolic equations was presented, which can be used to adjust space mesh reasonably. A numerical example is given to illustrate the accuracy and feasibility of this method.
基金Project supported by the State Key Laboratory of Scientific and Engineering Computing, Chinese Academy of Sciences and the National Natural Science Foundation of China （Grant No 10372053） and the Natural Science Foundation of Henan Province Government, China （Grant Nos 0311011400, 0511022200）.
文摘Non-Noether symmetries and conservative quantities of nonholonomic nonconservative dynamical systems are investigated in this paper. Based on the relationships among motion, nonconservative forces, nonholonomic constrained forces and Lagrangian, non-Noether symmetries and Lutzky conservative quantities are presented for nonholonomic nonconservative dynamical systems. The relation between non-Noether symmetry and Noether symmetry is discussed and it is further shown that non-Noether conservative quantities can be obtained by a complete set of Noether invariants. Finally, an example is given to illustrate these results.
基金浙江省自然科学基金，浙江省宁波市博士基金，the State Key Laboratory of Oil/Gas Reservoir Geology and Exploitation，Scientific Research Fund of Education Department of Zhejiang Province under
文摘A set of generalized symmetries with arbitrary functions of t for the Konopelchenko-Dubrovsky (KD)equation in 2+1 space dimensions is given by using a direct method called formal function series method presented by Lou. These symmetries constitute an infinite-dimensional generalized w∞ algebra.
文摘A generalized AKNS isospectral problem where the trace of corresponding spectral matrix is not zero, is transformed to a new isospectral problem where the trace of the resulting matrix is zero, by using transformation of Lax pairs, and these two spectral problems lead to the same hierarchy of equations. The authors started from the transformed spectral problem and constructed a new loop algebra which has not appeared before, and obtained the integrable coupling of the generalized AKNS hierarchy. Specially, the integrable couplings of the KdV equation and MKdV equation are obtained.
基金supported by China Postdoctoral Science Foundation grant 2020TQ0344the NSFC grants 11871139 and 12101597the NSF grants DMS-1720116,DMS-2012882,DMS-2011838,DMS-1719942,DMS-1913072.
文摘In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.
基金Project supported by State Key Laboratory of Scientific and Engineering Computing, Chinese Academy of Sciences and the National Natural Science Foundation of China （Grant Nos 10672143 and 10471145） and the Natural Science Foundation of Henan Province Government, China （Grant Nos 0311011400 and 0511022200）.
文摘This paper presents a discrete vaxiational principle and a method to build first-integrals for finite dimensional Lagrange-Maxwell mechanico-electrical systems with nonconservative forces and a dissipation function. The discrete variational principle and the corresponding Euler-Lagrange equations are derived from a discrete action associated to these systems. The first-integrals are obtained by introducing the infinitesimal transformation with respect to the generalized coordinates and electric quantities of the systems. This work also extends discrete Noether symmetries to mechanico-electrical dynamical systems. A practical example is presented to illustrate the results.
基金Project supported by the National Natural Science Foundation of China(Nos.11161017,11071251,and 10871099)the National Basic Research Program of China(973 Program)(No.2007CB209603)+1 种基金the Natural Science Foundation of Hainan Province(No.110002)the Scientific Research Foun-dation of Hainan University(No.kyqd1053)
文摘Based on splitting multi-symplectic structures, a new multi-symplectic scheme is proposed and applied to a nonlinear wave equation. The explicit multi-symplectic scheme of the nonlinear wave equation is obtained, and the corresponding multi-symplectic conservation property is proved. The backward error analysis shows that the explicit multi-symplectic scheme has good accuracy. The sine-Gordon equation and the Klein-Gordon equation are simulated by an explicit multi-symplectic scheme. The numerical results show that the new explicit multi-symplectic scheme can well simulate the solitary wave behaviors of the nonlinear wave equation and approximately preserve the relative energy error of the equation.
基金supported by the Chinese NSF grants 10571171 and 40233029.
文摘Stochastic approximation problem is to find some root or extremum of a nonlinear function for which only noisy measurements of the function are available. The classical algorithm for stochastic approximation problem is the Robbins-Monro （RM） algorithm, which uses the noisy evaluation of the negative gradient direction as the iterative direction. In order to accelerate the RM algorithm, this paper gives a flame algorithm using adaptive iterative directions. At each iteration, the new algorithm goes towards either the noisy evaluation of the negative gradient direction or some other directions under some switch criterions. Two feasible choices of the criterions are proposed and two corresponding frame algorithms are formed. Different choices of the directions under the same given switch criterion in the frame can also form different algorithms. We also proposed the simultanous perturbation difference forms for the two frame algorithms. The almost surely convergence of the new algorithms are all established. The numerical experiments show that the new algorithms are promising.
基金The project supported by the Natural Science Foundation of Zhejiang Province of China under Grant No. Y604036 and State Key Laboratory of 0il/Gas Reservoir Geology and Exploitation ＂PLN0402＂ The authors would like to thank Prof. Sen-Yue Lou for his help and discussion.
文摘The （3＋1）-dimensional Jimbo-Miwa （JM） equation is solved approximately by using the conformal invariant asymptotic expansion approach presented by Ruan. By solving the new （3＋1）-dimensional integrable models, which are conformal invariant and possess Painlevé property, the approximate solutions are obtained for the JM equation, containing not only one-soliton solutions but also periodic solutions and multi-soliton solutions. Some approximate solutions happen to be exact and some approximate solutions can become exact by choosing relations between the parameters properly.
基金National Natural Science Foundation of China under Grant No.10675065the Science Research Foundation of the Education Department of Zhejiang Province under Grant No.20070979+1 种基金the Natural Science Foundation of Zhejiang Province under Grant No.Y604036the State Key Laboratory of Oil/Gas Reservoir Geology and Exploitation\PLN0402
文摘Starting from the variable separation approach, the algebraic soliton solution and the solution describing the interaction between line soliton and algebraic soliton are obtained by selecting appropriate seed solution for （2＋1）-dimensional ANNV equation. The behaviors of interactions are discussed in detail both analytically and graphically. It is shown that there are two kinds of singular interactions between line soliton and algebraic soliton： 1） the resonant interaction where the algebraic soliton propagates together with the line soliton and persists infinitely; 2） the extremely repulsive interaction where the algebraic soliton affects the motion of the line soliton infinitely apart.
基金the National Natural Science Foundation of China(91641107,91852116,12071470)Fundamental Research of Civil Aircraft(MJ-F-2012-04)of Ministry of Industrialization and Information of China.
文摘In this paper,we develop a new sixth-order WENO scheme by adopting a convex combina-tion of a sixth-order global reconstruction and four low-order local reconstructions.Unlike the classical WENO schemes,the associated linear weights of the new scheme can be any positive numbers with the only requirement that their sum equals one.Further,a very simple smoothness indicator for the global stencil is proposed.The new scheme can achieve sixth-order accuracy in smooth regions.Numerical tests in some one-and two-dimensional bench-mark problems show that the new scheme has a little bit higher resolution compared with the recently developed sixth-order WENO-Z6 scheme,and it is more efficient than the classical fifth-order WENO-JS5 scheme and the recently developed sixth-order WENO6-S scheme.
文摘In this section, we introduce some outstanding research achievements made by CAS scientists in 2010, as listed in the Chinese version of Science Development Report published in March 2011.
