In this paper, we put our focus on a variable-coe^cient fifth-order Korteweg-de Vries （fKdV） equation, which possesses a great number of excellent properties and is of current importance in physical and engineering ...In this paper, we put our focus on a variable-coe^cient fifth-order Korteweg-de Vries （fKdV） equation, which possesses a great number of excellent properties and is of current importance in physical and engineering fields. Certain constraints are worked out, which make sure the integrability of such an equation. Under those constraints, some integrable properties are derived, such as the Lax pair and Darboux transformation. Via the Darboux transformation, which is an exercisable way to generate solutions in a recursive manner, the one- and two-solitonic solutions are presented and the relevant physical applications of these solitonic structures in some fields are also pointed out.展开更多
An exact two-soliton solution of discrete mKdv equation is derived by using the Hirota direct approach. In addition, we plot the soliton solutions to discuss the properties of solitons. It is worth while noting that w...An exact two-soliton solution of discrete mKdv equation is derived by using the Hirota direct approach. In addition, we plot the soliton solutions to discuss the properties of solitons. It is worth while noting that we obtain the completely elastic interaction between the two solitons.展开更多
In this paper,based on the forms and structures of Wronskian solutions to soliton equations,a Wronskianform expansion method is presented to find a new class of interaction solutions to the Kadomtsev-Petviashvili equa...In this paper,based on the forms and structures of Wronskian solutions to soliton equations,a Wronskianform expansion method is presented to find a new class of interaction solutions to the Kadomtsev-Petviashvili equation.One characteristic of the method is that Wronskian entries do not satisfy linear partial differential equation.展开更多
Based on a first-order nonlinear ordinary differential equation with six-degree nonlinear term, we first present a new auxiliary equation expansion method and its algorithm. Being concise and straightforward, the meth...Based on a first-order nonlinear ordinary differential equation with six-degree nonlinear term, we first present a new auxiliary equation expansion method and its algorithm. Being concise and straightforward, the method is applied to the Kundu equation. As a result, some new exact travelling wave solutions are obtained, which include bright and dark solitary wave solutions, triangular periodic wave solutions, and singular solutions. This algorithm can also be applied to other nonlinear evolution equations in mathematical physics.展开更多
In this paper, we investigate symmetries of the new （4＋1）-dimensional Fokas equation, including point symmetries and the potential symmetries. We firstly employ the algorithmic procedure of computing the point symm...In this paper, we investigate symmetries of the new （4＋1）-dimensional Fokas equation, including point symmetries and the potential symmetries. We firstly employ the algorithmic procedure of computing the point symmetries. And then we transform the Fokas equation into a potential system and gain the potential symmetries of Fokas equation. Finally, we use the obtained point symmetries wave solutions and other solutions of the Fokas equation. and some constructive methods to get some doubly periodic In particular, some solitary wave solutions are also given.展开更多
Steady-state non-dominated sorting genetic algorithm (SNSGA), a new form of multi-objective genetic algorithm, is implemented by combining the steady-state idea in steady-state genetic algorithms (SSGA) and the fitnes...Steady-state non-dominated sorting genetic algorithm (SNSGA), a new form of multi-objective genetic algorithm, is implemented by combining the steady-state idea in steady-state genetic algorithms (SSGA) and the fitness assignment strategy of non-dominated sorting genetic algorithm (NSGA). The fitness assignment strategy is improved and a new self-adjustment scheme of is proposed. This algorithm is proved to be very efficient both computationally and in terms of the quality of the Pareto fronts produced with five test problems including GA difficult problem and GA deceptive one. Finally, SNSGA is introduced to solve multi-objective mixed integer linear programming (MILP) and mixed integer non-linear programming (MINLP) problems in process synthesis.展开更多
Finding all zeros of polynomial systems is very interesting and it is also useul for many applied science problems.In this paper,based on Wu's method,we give an algorithm to find all isolated zeros of polynomial s...Finding all zeros of polynomial systems is very interesting and it is also useul for many applied science problems.In this paper,based on Wu's method,we give an algorithm to find all isolated zeros of polynomial systems (or polynomial equations).By solving Lorenz equations,it is shown that our algo-rithm is efficient and powerful.展开更多
In this work we present a new method to solve the Perona Malik equation for the image denoising. The method is based on a modified fixed point algorithm which is fast and stable. We discretize the equation using a fin...In this work we present a new method to solve the Perona Malik equation for the image denoising. The method is based on a modified fixed point algorithm which is fast and stable. We discretize the equation using a finite volume method by integrating the equation using a fuzzy measure on the control volume. To make our algorithm move faster in time, we have used an optimized domain decomposition which generalize the wave relaxation method. Several test of noised images illustrate this approach and show the efficiency of the proposed new method.展开更多
In this paper we introduce an implementation for the efficient numericalsolution of exterior initial boundary value problem for parabolic equation. The problemis reformulated as an equivalent one on a boundary T using...In this paper we introduce an implementation for the efficient numericalsolution of exterior initial boundary value problem for parabolic equation. The problemis reformulated as an equivalent one on a boundary T using natural boundary reduction.The governing equation is first discretized in time, leading to a time-stepping scheme,where an exterior elliptic problem has to be solved in each time step. By Fourier ex-pansion, we derive a natural integral equation of the elliptic problem related to timestep and Poisson integral integral formula over exterior circular domain. Finite elementdiscretization of the natural integral equation is employed to solve this problem. Thecomputational aspects of this method are discussed. Numerical results are presented toillustrate feasibility and efficiency of our method.展开更多
A new algorithm for linear instantaneous independent component analysis is proposed based on maximizing the log-likelihood contrast function which can be changed into a gradient equation.An iterative method is introdu...A new algorithm for linear instantaneous independent component analysis is proposed based on maximizing the log-likelihood contrast function which can be changed into a gradient equation.An iterative method is introduced to solve this equation efficiently.The unknown probability density functions as well as their first and second derivatives in the gradient equation are estimated by kernel density method.Computer simulations on artificially generated signals and gray scale natural scene images confirm the efficiency and accuracy of the proposed algorithm.展开更多
A new expanded approach is presented to find exact solutions of nonlinear differential-difference equations. As its application, the soliton solutions and periodic solutions of a lattice equation are obtained.
Following an order analysis of key parameters, a decoupled procedure for simulation of convection-radiation heat transfer problems in supersonic combustion ramjet(scramjet) engine was developed. The radiation module o...Following an order analysis of key parameters, a decoupled procedure for simulation of convection-radiation heat transfer problems in supersonic combustion ramjet(scramjet) engine was developed. The radiation module of the procedure consisted of Perry 5GG weighted sum gray gases model for spectral property calculation and discrete ordinates method S4 scheme for radiative transfer computation, while the flow field was computed using the Favrè average conservative Navier-Stokes(N-S) equations, in conjunction with Menter's k-ω SST two-equation model. A series of 2D supersonic nonreactive turbulent channel flows of radiative participants with selective parameters were simulated for validation purpose. Radiative characteristics in DLR hydrogen fueled and NASA SCHOLAR ethylene fueled scramjets were numerically studied using the developed procedure. The results indicated that the variations of spatial distributions of the radiative source and total absorption coefficient are highly consistent with those of the temperature and radiative participants, while the spatial distribution of the incident radiation spreads wider. It also demonstrated that the convective heating is significantly affected by the complexity of the flow field, such as the shock wave/boundary layer interactions, while the radiative heating is simply an integral effect of the whole flow field. Although the radiative heating in the combustion chambers reaches a certain level, an order of magnitude of 10 k W/m2, it still contributes little to the total heat transfer(<7%).展开更多
In this paper we study the algorithms and their parallel implementation for solving large-scale generalized eigenvalue problems in modal analysis.Three predominant subspace algorithms,i.e.,Krylov-Schur method,implicit...In this paper we study the algorithms and their parallel implementation for solving large-scale generalized eigenvalue problems in modal analysis.Three predominant subspace algorithms,i.e.,Krylov-Schur method,implicitly restarted Arnoldi method and Jacobi-Davidson method,are modified with some complementary techniques to make them suitable for modal analysis.Detailed descriptions of the three algorithms are given.Based on these algorithms,a parallel solution procedure is established via the PANDA framework and its associated eigensolvers.Using the solution procedure on a machine equipped with up to 4800processors,the parallel performance of the three predominant methods is evaluated via numerical experiments with typical engineering structures,where the maximum testing scale attains twenty million degrees of freedom.The speedup curves for different cases are obtained and compared.The results show that the three methods are good for modal analysis in the scale of ten million degrees of freedom with a favorable parallel scalability.展开更多
It is presented in this paper that the new design and its analysis of finite difference domain decomposition algorithms for the two-dimensional heat equation, and the numerical results have shown the stability and acc...It is presented in this paper that the new design and its analysis of finite difference domain decomposition algorithms for the two-dimensional heat equation, and the numerical results have shown the stability and accuracy of the algorithms, where SauFyev asymmetric schemes have been used at the interface points. The Algorithm II in this paper has further extended those developed by Dawson and the others, Zhang and Shen.展开更多
Motivated by Sasaki＇s work on the extended Hensel construction for solving multivariate algebraic equations, we present a generalized Hensel lifting, which takes advantage of sparsity, for factoring bivariate polynom...Motivated by Sasaki＇s work on the extended Hensel construction for solving multivariate algebraic equations, we present a generalized Hensel lifting, which takes advantage of sparsity, for factoring bivariate polynomial over the rational number field. Another feature of the factorization algorithm presented in this article is a new recombination method, which can solve the extraneous factor problem before lifting based on numerical linear algebra. Both theoretical analysis and experimental data show that the algorithm is etIicient, especially for sparse bivariate polynomials.展开更多
基金The project supported by the Key Project of the Chinese Ministry of Education under Grant No.106033the Specialized Research Fund for the Doctoral Program of Higher Education under Grant No.20060006024+2 种基金Chinese Ministry of Education,the National Natural Science Foundation of China under Grant Nos.60772023 and 60372095the Open Fund of the State Key Laboratory of Software Development Environment under Grant No.SKLSDE-07-001Beijing University of Aeronautics and Astronautics,and by the National Basic Research Program of China(973 Program)under Grant No.2005CB321901
文摘In this paper, we put our focus on a variable-coe^cient fifth-order Korteweg-de Vries （fKdV） equation, which possesses a great number of excellent properties and is of current importance in physical and engineering fields. Certain constraints are worked out, which make sure the integrability of such an equation. Under those constraints, some integrable properties are derived, such as the Lax pair and Darboux transformation. Via the Darboux transformation, which is an exercisable way to generate solutions in a recursive manner, the one- and two-solitonic solutions are presented and the relevant physical applications of these solitonic structures in some fields are also pointed out.
文摘An exact two-soliton solution of discrete mKdv equation is derived by using the Hirota direct approach. In addition, we plot the soliton solutions to discuss the properties of solitons. It is worth while noting that we obtain the completely elastic interaction between the two solitons.
基金Supported by the Young Teachers Science Foundation of Beijing University of Civil Engineering and Architecture under Grant No.100602707
文摘In this paper,based on the forms and structures of Wronskian solutions to soliton equations,a Wronskianform expansion method is presented to find a new class of interaction solutions to the Kadomtsev-Petviashvili equation.One characteristic of the method is that Wronskian entries do not satisfy linear partial differential equation.
文摘Based on a first-order nonlinear ordinary differential equation with six-degree nonlinear term, we first present a new auxiliary equation expansion method and its algorithm. Being concise and straightforward, the method is applied to the Kundu equation. As a result, some new exact travelling wave solutions are obtained, which include bright and dark solitary wave solutions, triangular periodic wave solutions, and singular solutions. This algorithm can also be applied to other nonlinear evolution equations in mathematical physics.
基金supported by the National Natural Science Foundation of China under Grant No.60821002the National Key Basic Research Program of China under Grant No.2004CB318000
文摘In this paper, we investigate symmetries of the new （4＋1）-dimensional Fokas equation, including point symmetries and the potential symmetries. We firstly employ the algorithmic procedure of computing the point symmetries. And then we transform the Fokas equation into a potential system and gain the potential symmetries of Fokas equation. Finally, we use the obtained point symmetries wave solutions and other solutions of the Fokas equation. and some constructive methods to get some doubly periodic In particular, some solitary wave solutions are also given.
文摘Steady-state non-dominated sorting genetic algorithm (SNSGA), a new form of multi-objective genetic algorithm, is implemented by combining the steady-state idea in steady-state genetic algorithms (SSGA) and the fitness assignment strategy of non-dominated sorting genetic algorithm (NSGA). The fitness assignment strategy is improved and a new self-adjustment scheme of is proposed. This algorithm is proved to be very efficient both computationally and in terms of the quality of the Pareto fronts produced with five test problems including GA difficult problem and GA deceptive one. Finally, SNSGA is introduced to solve multi-objective mixed integer linear programming (MILP) and mixed integer non-linear programming (MINLP) problems in process synthesis.
文摘Finding all zeros of polynomial systems is very interesting and it is also useul for many applied science problems.In this paper,based on Wu's method,we give an algorithm to find all isolated zeros of polynomial systems (or polynomial equations).By solving Lorenz equations,it is shown that our algo-rithm is efficient and powerful.
文摘In this work we present a new method to solve the Perona Malik equation for the image denoising. The method is based on a modified fixed point algorithm which is fast and stable. We discretize the equation using a finite volume method by integrating the equation using a fuzzy measure on the control volume. To make our algorithm move faster in time, we have used an optimized domain decomposition which generalize the wave relaxation method. Several test of noised images illustrate this approach and show the efficiency of the proposed new method.
基金National Natural Science Foundation of China(19701001)
文摘In this paper we introduce an implementation for the efficient numericalsolution of exterior initial boundary value problem for parabolic equation. The problemis reformulated as an equivalent one on a boundary T using natural boundary reduction.The governing equation is first discretized in time, leading to a time-stepping scheme,where an exterior elliptic problem has to be solved in each time step. By Fourier ex-pansion, we derive a natural integral equation of the elliptic problem related to timestep and Poisson integral integral formula over exterior circular domain. Finite elementdiscretization of the natural integral equation is employed to solve this problem. Thecomputational aspects of this method are discussed. Numerical results are presented toillustrate feasibility and efficiency of our method.
文摘A new algorithm for linear instantaneous independent component analysis is proposed based on maximizing the log-likelihood contrast function which can be changed into a gradient equation.An iterative method is introduced to solve this equation efficiently.The unknown probability density functions as well as their first and second derivatives in the gradient equation are estimated by kernel density method.Computer simulations on artificially generated signals and gray scale natural scene images confirm the efficiency and accuracy of the proposed algorithm.
基金the National Natural Science Foundation of China (No. 60773119)
文摘A new expanded approach is presented to find exact solutions of nonlinear differential-difference equations. As its application, the soliton solutions and periodic solutions of a lattice equation are obtained.
基金supported by the National Natural Science Foundation of China(Grant No.11202014)
文摘Following an order analysis of key parameters, a decoupled procedure for simulation of convection-radiation heat transfer problems in supersonic combustion ramjet(scramjet) engine was developed. The radiation module of the procedure consisted of Perry 5GG weighted sum gray gases model for spectral property calculation and discrete ordinates method S4 scheme for radiative transfer computation, while the flow field was computed using the Favrè average conservative Navier-Stokes(N-S) equations, in conjunction with Menter's k-ω SST two-equation model. A series of 2D supersonic nonreactive turbulent channel flows of radiative participants with selective parameters were simulated for validation purpose. Radiative characteristics in DLR hydrogen fueled and NASA SCHOLAR ethylene fueled scramjets were numerically studied using the developed procedure. The results indicated that the variations of spatial distributions of the radiative source and total absorption coefficient are highly consistent with those of the temperature and radiative participants, while the spatial distribution of the incident radiation spreads wider. It also demonstrated that the convective heating is significantly affected by the complexity of the flow field, such as the shock wave/boundary layer interactions, while the radiative heating is simply an integral effect of the whole flow field. Although the radiative heating in the combustion chambers reaches a certain level, an order of magnitude of 10 k W/m2, it still contributes little to the total heat transfer(<7%).
基金supported by the National Defence Basic Fundamental Research Program of China(Grant No.C1520110002)the Fundamental Development Foundation of China Academy Engineering Physics(Grant No.2012A0202008)
文摘In this paper we study the algorithms and their parallel implementation for solving large-scale generalized eigenvalue problems in modal analysis.Three predominant subspace algorithms,i.e.,Krylov-Schur method,implicitly restarted Arnoldi method and Jacobi-Davidson method,are modified with some complementary techniques to make them suitable for modal analysis.Detailed descriptions of the three algorithms are given.Based on these algorithms,a parallel solution procedure is established via the PANDA framework and its associated eigensolvers.Using the solution procedure on a machine equipped with up to 4800processors,the parallel performance of the three predominant methods is evaluated via numerical experiments with typical engineering structures,where the maximum testing scale attains twenty million degrees of freedom.The speedup curves for different cases are obtained and compared.The results show that the three methods are good for modal analysis in the scale of ten million degrees of freedom with a favorable parallel scalability.
基金supported by Natural Science Foundation of China under Grant Nos.10671060 and 10871061the Youth Foundation of Hunan Education Bureau under Grant No.06B037+1 种基金the Natural Science Foundation of Hunan Province under Grant No.09JJ6015the Construct Program of the Key Discipline in Hunan Province
文摘It is presented in this paper that the new design and its analysis of finite difference domain decomposition algorithms for the two-dimensional heat equation, and the numerical results have shown the stability and accuracy of the algorithms, where SauFyev asymmetric schemes have been used at the interface points. The Algorithm II in this paper has further extended those developed by Dawson and the others, Zhang and Shen.
基金supported by National Natural Science Foundation of China(GrantNos.91118001 and 11170153)National Key Basic Research Project of China(Grant No.2011CB302400)Chongqing Science and Technology Commission Project(Grant No.cstc2013jjys40001)
文摘Motivated by Sasaki＇s work on the extended Hensel construction for solving multivariate algebraic equations, we present a generalized Hensel lifting, which takes advantage of sparsity, for factoring bivariate polynomial over the rational number field. Another feature of the factorization algorithm presented in this article is a new recombination method, which can solve the extraneous factor problem before lifting based on numerical linear algebra. Both theoretical analysis and experimental data show that the algorithm is etIicient, especially for sparse bivariate polynomials.
