In this paper, some exact solutions of the (3 + 1)-dimensional variable-coefficient Yu-Toda-Sasa-Fukuyama equation are investigated. By using Hirota’s direct method and symbolic computation, we obtained N-soliton sol...In this paper, some exact solutions of the (3 + 1)-dimensional variable-coefficient Yu-Toda-Sasa-Fukuyama equation are investigated. By using Hirota’s direct method and symbolic computation, we obtained N-soliton solution. By using the long wave limit method, the N-order rational solution can be obtained from N-order soliton solution. Then, through the paired complexification of parameters, the lump solution is obtained from N-order rational solution. Meanwhile, we obtained a hybrid solution between 1-lump solution and N-soliton (N=1,2) by using the long wave limit method and parameter complex. Furthermore, four different sets of three-dimensional graphs of solitons, lump solutions and hybrid solutions are drawn by selecting four different sets of coefficient functions which include one set of constant coefficient function and three sets of variable coefficient functions.展开更多
In this paper, the evolutionary behavior of N-solitons for a (2 + 1)-dimensional Konopelchenko-Dubrovsky equations is studied by using the Hirota bilinear method and the long wave limit method. Based on the N-soliton ...In this paper, the evolutionary behavior of N-solitons for a (2 + 1)-dimensional Konopelchenko-Dubrovsky equations is studied by using the Hirota bilinear method and the long wave limit method. Based on the N-soliton solution, we first study the evolution from N-soliton to T-order (T=1,2) breather wave solutions via the paired-complexification of parameters, and then we get the N-order rational solutions, M-order (M=1,2) lump solutions, and the hybrid behavior between a variety of different types of solitons combined with the parameter limit technique and the paired-complexification of parameters. Meanwhile, we also provide a large number of three-dimensional figures in order to better show the degeneration of the N-soliton and the interaction behavior between different N-solitons.展开更多
A fast hybrid algorithm based on gridless method coupled with finite volume method (FVM) is developed for the solution to Euler equations. Compared with pure gridless method, the efficiency of the hybrid algorithm i...A fast hybrid algorithm based on gridless method coupled with finite volume method (FVM) is developed for the solution to Euler equations. Compared with pure gridless method, the efficiency of the hybrid algorithm is improved to the level of finite volume method for most parts of the flow filed arc covered with grid cells. Moreover, the hybrid method is flexible to deal with the configurations as clouds of points are used to cover the region adjacent to the bodies. Mirror satellites and mirror grid cells arc introduced to the interface to accomplish data communication between the different parts of the flow field. The Euler Equations arc spatially discretized with finite volume method and gridless method in mesh and clouds of points respectively, and an explicit four-stage Runge-Kutta scheme is utilized to reach the steady-state solution. Internal flows in channels and external flows over airfoils arc investigated with hybrid method, and the solutions arc comparad to those using pure finite volume method and pure gridless method. Numerical examples show that the hybrid algorithm captures the shock waves accurately, and it is as efficient as fmite volume method.展开更多
The paper develops exponential stability of the analytic solution and convergence in probability of the numerical method for highly nonlinear hybrid stochastic pantograph equation. The classical linear growth conditio...The paper develops exponential stability of the analytic solution and convergence in probability of the numerical method for highly nonlinear hybrid stochastic pantograph equation. The classical linear growth condition is replaced by polynomial growth conditions, under which there exists a unique global solution and the solution is almost surely exponentially stable. On the basis of a series of lemmas, the paper establishes a new criterion on convergence in probability of the Euler-Maruyama approximate solution. The criterion is very general so that many highly nonlinear stochastic pantograph equations can obey these conditions. A highly nonlinear example is provided to illustrate the main theory.展开更多
This paper proposes numerical methods for solving hybrid weakly singular integro-differential equations of the second kind. The terms in these equations are in the following order: derivative term of a state, integro-...This paper proposes numerical methods for solving hybrid weakly singular integro-differential equations of the second kind. The terms in these equations are in the following order: derivative term of a state, integro-differential term of a state with a weakly singular kernel, a state, integral term of a state with a smooth kernel, and force. The original class of weakly singular integro-differential equations of the first kind is derived from aeroelasticity mathematical models. Among the proposed methods, the method for solving linear cases is fully based on previously reported approximation scheme for equations of the first kind. For nonlinear cases, a revised method is proposed. Examples are presented to demonstrate the effectiveness of the proposed methods, and the results indicate that the proposed methods facilitate achieving satisfactory and accurate approximations.展开更多
In this paper, a block method with one hybrid point for solving Jerk equations is presented. The hybrid point is chosen to optimize the local truncation errors of the main formulas for the solution and the derivative ...In this paper, a block method with one hybrid point for solving Jerk equations is presented. The hybrid point is chosen to optimize the local truncation errors of the main formulas for the solution and the derivative at the end of the block. Analysis of the method is discussed, and some numerical examples show that the proposed method is efficient and accurate.展开更多
The hybrid lattice, known as a discrete Korteweg-de Vries (KdV) equation, is found to be a discrete modified Korteweg-de Vries (mKdV) equation in this paper. The coupled hybrid lattice, which is pointed to be a discre...The hybrid lattice, known as a discrete Korteweg-de Vries (KdV) equation, is found to be a discrete modified Korteweg-de Vries (mKdV) equation in this paper. The coupled hybrid lattice, which is pointed to be a discrete coupled KdV system, is also found to be discrete form of a coupled mKdV systems. Delayed differential reduction system and pure difference systems are derived from the coupled hybrid system by means of the symmetry reduction approach. Cnoidal wave, positon and negaton solutions for the coupled hybrid system are proposed.展开更多
A special transformation is introduced and thereby leads to the N-soliton solution of the(2+1)-dimensional generalized Konopelchenko-Dubrovsky-Kaup-Kupershmidt(KDKK) equation.Then,by employing the long wave limit and ...A special transformation is introduced and thereby leads to the N-soliton solution of the(2+1)-dimensional generalized Konopelchenko-Dubrovsky-Kaup-Kupershmidt(KDKK) equation.Then,by employing the long wave limit and imposing complex conjugate constraints to the related solitons,various localized interaction solutions are constructed,including the general M-lumps,T-breathers,and hybrid wave solutions.Dynamical behaviors of these solutions are investigated analytically and graphically.The solutions obtained are very helpful in studying the interaction phenomena of nonlinear localized waves.Therefore,we hope these results can provide some theoretical guidance to the experts in oceanography,atmospheric science,and weather forecasting.展开更多
Certain hybrid prototypes of dispersive optical solitons that we are looking for can correspond to new or future behaviors, observable or not, developed or will be developed by optical media that present the cubic-qui...Certain hybrid prototypes of dispersive optical solitons that we are looking for can correspond to new or future behaviors, observable or not, developed or will be developed by optical media that present the cubic-quintic-septic law coupled, with strong dispersions. The equation considered for this purpose is that of non-linear Schrödinger. The solutions are obtained using the Bogning-Djeumen Tchaho-Kofané method extended to the new implicit Bogning’ functions. Some of the obtained solutions show that their existence is due only to the Kerr law nonlinearity presence. Graphical representations plotted have confirmed the hybrid and multi-form character of the obtained dispersive optical solitons. We believe that a good understanding of the hybrid dispersive optical solitons highlighted in the context of this work allows to grasp the physical description of systems whose dynamics are governed by nonlinear Schrödinger equation as studied in this work, allowing thereby a relevant improvement of complex problems encountered in particular in nonliear optaics and in optical fibers.展开更多
The well-known Riccati differential equations play a key role in many fields,including problems in protein folding,control and stabilization,stochastic control,and cybersecurity(risk analysis and malware propaga-tion)...The well-known Riccati differential equations play a key role in many fields,including problems in protein folding,control and stabilization,stochastic control,and cybersecurity(risk analysis and malware propaga-tion).Quantum computer algorithms have the potential to implement faster approximate solutions to the Riccati equations compared with strictly classical algorithms.While systems with many qubits are still under development,there is significant interest in developing algorithms for near-term quantum computers to determine their accuracy and limitations.In this paper,we propose a hybrid quantum-classical algorithm,the Matrix Riccati Solver(MRS).This approach uses a transformation of variables to turn a set of nonlinear differential equation into a set of approximate linear differential equations(i.e.,second order non-constant coefficients)which can in turn be solved using a version of the Harrow-Hassidim-Lloyd(HHL)quantum algorithm for the case of Hermitian matrices.We implement this approach using the Qiskit language and compute near-term results using a 4 qubit IBM Q System quantum computer.Comparisons with classical results and areas for future research are discussed.展开更多
Hybrid grids are used for the solution of 3D turbulence Navier Stokes equations. The prismatic grids are generated near the wall, and the tetrahedron grids are generated in the other field. A Navier Stokes solver usin...Hybrid grids are used for the solution of 3D turbulence Navier Stokes equations. The prismatic grids are generated near the wall, and the tetrahedron grids are generated in the other field. A Navier Stokes solver using Jamson′s finite volume method is developed. The algebraic Baldwin Lomax turbulence model is adopted. The numerical tests show that the above method is very efficient.展开更多
A total variation diminishing-weighted average flux (TVD-WAF)-based hybrid numerical scheme for the enhanced version of nonlinearly dispersive Boussinesq-type equations was developed. The one-dimensional governing e...A total variation diminishing-weighted average flux (TVD-WAF)-based hybrid numerical scheme for the enhanced version of nonlinearly dispersive Boussinesq-type equations was developed. The one-dimensional governing equations were rewritten in the conservative form and then discretized on a uniform grid. The finite volume method was used to discretize the flux term while the remaining terms were approximated with the finite difference method. The second-order TVD-WAF method was employed in conjunction with the Harten-Lax-van Leer (HLL) Riemann solver to calculate the numerical flux, and the variables at the cell interface for the local Riemann problem were reconstructed via the fourth- order monotone upstream-centered scheme for conservation laws (MUSCL). The time marching scheme based on the third-order TVD Runge- Kutta method was used to obtain numerical solutions. The model was validated through a series of numerical tests, in which wave breaking and a moving shoreline were treated. The good agreement between the computed results, documented analytical solutions, and experimental data demonstrates the correct discretization of the governing equations and high accuracy of the proposed scheme, and also conforms the advantages of the proposed shock-capturing scheme for the enhanced version of the Boussinesq model, including the convenience in the treatment of wave breaking and moving shorelines and without the need for a numerical filter.展开更多
In this study, a revised version of some numerical methods for a class of hybrid integro-differential equations with weakly singular kernels (Abel types) is presented. These equations were developed from a class of in...In this study, a revised version of some numerical methods for a class of hybrid integro-differential equations with weakly singular kernels (Abel types) is presented. These equations were developed from a class of integro-differential equations of first kind originating from an aeroelasticity problem. By manipulating the bounds of initial conditions with random variations, this study numerically demonstrated the well-posedness properties of the equations. Finally, an assumption of separating variables, allowed for linear splines to be chosen as a basis and for the differentiation and integration of the integro-differential part to be interchanged;hence, a numerical scheme was constructed.展开更多
We propose new hybrid Lagrange neural networks called LaNets to predict the numerical solutions of partial differential equations.That is,we embed Lagrange interpolation and small sample learning into deep neural netw...We propose new hybrid Lagrange neural networks called LaNets to predict the numerical solutions of partial differential equations.That is,we embed Lagrange interpolation and small sample learning into deep neural network frameworks.Concretely,we first perform Lagrange interpolation in front of the deep feedforward neural network.The Lagrange basis function has a neat structure and a strong expression ability,which is suitable to be a preprocessing tool for pre-fitting and feature extraction.Second,we introduce small sample learning into training,which is beneficial to guide themodel to be corrected quickly.Taking advantages of the theoretical support of traditional numerical method and the efficient allocation of modern machine learning,LaNets achieve higher predictive accuracy compared to the state-of-the-artwork.The stability and accuracy of the proposed algorithmare demonstrated through a series of classical numerical examples,including one-dimensional Burgers equation,onedimensional carburizing diffusion equations,two-dimensional Helmholtz equation and two-dimensional Burgers equation.Experimental results validate the robustness,effectiveness and flexibility of the proposed algorithm.展开更多
As a boundary-type meshless method, the singular hybrid boundary node method(SHBNM) is based on the modified variational principle and the moving least square(MLS) approximation, so it has the advantages of both b...As a boundary-type meshless method, the singular hybrid boundary node method(SHBNM) is based on the modified variational principle and the moving least square(MLS) approximation, so it has the advantages of both boundary element method(BEM) and meshless method. In this paper, the dual reciprocity method(DRM) is combined with SHBNM to solve Poisson equation in which the solution is divided into particular solution and general solution. The general solution is achieved by means of SHBNM, and the particular solution is approximated by using the radial basis function(RBF). Only randomly distributed nodes on the bounding surface of the domain are required and it doesn't need extra equations to compute internal parameters in the domain. The postprocess is very simple. Numerical examples for the solution of Poisson equation show that high convergence rates and high accuracy with a small node number are achievable.展开更多
In this paper, we are going to present a class of nonlinear equation solving methods. Steffensen’s method is a simple method for solving a nonlinear equation. By using Steffensen’s method and by combining this metho...In this paper, we are going to present a class of nonlinear equation solving methods. Steffensen’s method is a simple method for solving a nonlinear equation. By using Steffensen’s method and by combining this method with it, we obtain a new method. It can be said that this method, due to not using the function derivative, would be a good method for solving the nonlinear equation compared to Newton’s method. Finally, we will see that Newton’s method and Steffensen’s hybrid method both have a two-order convergence.展开更多
We devise hybrid high-order(HHO)methods for the acoustic wave equation in the time domain.We frst consider the second-order formulation in time.Using the Newmark scheme for the temporal discretization,we show that the...We devise hybrid high-order(HHO)methods for the acoustic wave equation in the time domain.We frst consider the second-order formulation in time.Using the Newmark scheme for the temporal discretization,we show that the resulting HHO-Newmark scheme is energy-conservative,and this scheme is also amenable to static condensation at each time step.We then consider the formulation of the acoustic wave equation as a frst-order system together with singly-diagonally implicit and explicit Runge-Kutta(SDIRK and ERK)schemes.HHO-SDIRK schemes are amenable to static condensation at each time step.For HHO-ERK schemes,the use of the mixed-order formulation,where the polynomial degree of the cell unknowns is one order higher than that of the face unknowns,is key to beneft from the explicit structure of the scheme.Numerical results on test cases with analytical solutions show that the methods can deliver optimal convergence rates for smooth solutions of order O(hk+1)in the H1-norm and of order O(h^(k+2))in the L^(2)-norm.Moreover,test cases on wave propagation in heterogeneous media indicate the benefts of using high-order methods.展开更多
An L-stable block method based on hybrid second derivative algorithm (BHSDA) is provided by a continuous second derivative method that is defined for all values of the independent variable and applied to parabolic par...An L-stable block method based on hybrid second derivative algorithm (BHSDA) is provided by a continuous second derivative method that is defined for all values of the independent variable and applied to parabolic partial differential equations (PDEs). The use of the BHSDA to solve PDEs is facilitated by the method of lines which involves making an approximation to the space derivatives, and hence reducing the problem to that of solving a time-dependent system of first order initial value ordinary differential equations. The stability properties of the method is examined and some numerical results presented.展开更多
In this paper, we have successfully extended the Jacobian elliptic function expansion approach to nonlinear differential-difference equations. The Hybrid lattice equation is chosen to illustrate this approach. As a co...In this paper, we have successfully extended the Jacobian elliptic function expansion approach to nonlinear differential-difference equations. The Hybrid lattice equation is chosen to illustrate this approach. As a consequence, twelve families of Jacobian elliptic function solutions with different parameters of the Hybrid lattice equation are obtained. When the modulus m → 1 or O, doubly-periodic solutions degenerate to solitonic solutions and trigonometric function solutions, respectively.展开更多
In this paper,Chebyshev pseudospectral-finite element schemes are proposed for solving three dimensional vorticity equation.Some approximation results in nonisotropic Sobolev spaces are given.The generalized stability...In this paper,Chebyshev pseudospectral-finite element schemes are proposed for solving three dimensional vorticity equation.Some approximation results in nonisotropic Sobolev spaces are given.The generalized stability and the convergence are proved strictly.The numerical results show the advantages of this method.The technique in this paper is also applicable to other three-dimensional nonlinear problems in fluid dynamics.展开更多
文摘In this paper, some exact solutions of the (3 + 1)-dimensional variable-coefficient Yu-Toda-Sasa-Fukuyama equation are investigated. By using Hirota’s direct method and symbolic computation, we obtained N-soliton solution. By using the long wave limit method, the N-order rational solution can be obtained from N-order soliton solution. Then, through the paired complexification of parameters, the lump solution is obtained from N-order rational solution. Meanwhile, we obtained a hybrid solution between 1-lump solution and N-soliton (N=1,2) by using the long wave limit method and parameter complex. Furthermore, four different sets of three-dimensional graphs of solitons, lump solutions and hybrid solutions are drawn by selecting four different sets of coefficient functions which include one set of constant coefficient function and three sets of variable coefficient functions.
文摘In this paper, the evolutionary behavior of N-solitons for a (2 + 1)-dimensional Konopelchenko-Dubrovsky equations is studied by using the Hirota bilinear method and the long wave limit method. Based on the N-soliton solution, we first study the evolution from N-soliton to T-order (T=1,2) breather wave solutions via the paired-complexification of parameters, and then we get the N-order rational solutions, M-order (M=1,2) lump solutions, and the hybrid behavior between a variety of different types of solitons combined with the parameter limit technique and the paired-complexification of parameters. Meanwhile, we also provide a large number of three-dimensional figures in order to better show the degeneration of the N-soliton and the interaction behavior between different N-solitons.
基金Aeronautical Science Foundation of China (02A52002), National Natural Science Foundation of China(10372043)
文摘A fast hybrid algorithm based on gridless method coupled with finite volume method (FVM) is developed for the solution to Euler equations. Compared with pure gridless method, the efficiency of the hybrid algorithm is improved to the level of finite volume method for most parts of the flow filed arc covered with grid cells. Moreover, the hybrid method is flexible to deal with the configurations as clouds of points are used to cover the region adjacent to the bodies. Mirror satellites and mirror grid cells arc introduced to the interface to accomplish data communication between the different parts of the flow field. The Euler Equations arc spatially discretized with finite volume method and gridless method in mesh and clouds of points respectively, and an explicit four-stage Runge-Kutta scheme is utilized to reach the steady-state solution. Internal flows in channels and external flows over airfoils arc investigated with hybrid method, and the solutions arc comparad to those using pure finite volume method and pure gridless method. Numerical examples show that the hybrid algorithm captures the shock waves accurately, and it is as efficient as fmite volume method.
基金support from the National Natural Science Foundation of China(70871046,71171091,71191091)Fundamental Research Funds for the Central Universities(2011QN167)
文摘The paper develops exponential stability of the analytic solution and convergence in probability of the numerical method for highly nonlinear hybrid stochastic pantograph equation. The classical linear growth condition is replaced by polynomial growth conditions, under which there exists a unique global solution and the solution is almost surely exponentially stable. On the basis of a series of lemmas, the paper establishes a new criterion on convergence in probability of the Euler-Maruyama approximate solution. The criterion is very general so that many highly nonlinear stochastic pantograph equations can obey these conditions. A highly nonlinear example is provided to illustrate the main theory.
文摘This paper proposes numerical methods for solving hybrid weakly singular integro-differential equations of the second kind. The terms in these equations are in the following order: derivative term of a state, integro-differential term of a state with a weakly singular kernel, a state, integral term of a state with a smooth kernel, and force. The original class of weakly singular integro-differential equations of the first kind is derived from aeroelasticity mathematical models. Among the proposed methods, the method for solving linear cases is fully based on previously reported approximation scheme for equations of the first kind. For nonlinear cases, a revised method is proposed. Examples are presented to demonstrate the effectiveness of the proposed methods, and the results indicate that the proposed methods facilitate achieving satisfactory and accurate approximations.
文摘In this paper, a block method with one hybrid point for solving Jerk equations is presented. The hybrid point is chosen to optimize the local truncation errors of the main formulas for the solution and the derivative at the end of the block. Analysis of the method is discussed, and some numerical examples show that the proposed method is efficient and accurate.
基金Supported by the Natural Science Foundation of Guangdong Province of China under Grant No. 10452840301004616the National Natural Science Foundation of China under Grant No. 61001018the Scientific Research Foundation for the Doctors of University of Electronic Science and Technology of China Zhongshan Institute under Grant No. 408YKQ09
文摘The hybrid lattice, known as a discrete Korteweg-de Vries (KdV) equation, is found to be a discrete modified Korteweg-de Vries (mKdV) equation in this paper. The coupled hybrid lattice, which is pointed to be a discrete coupled KdV system, is also found to be discrete form of a coupled mKdV systems. Delayed differential reduction system and pure difference systems are derived from the coupled hybrid system by means of the symmetry reduction approach. Cnoidal wave, positon and negaton solutions for the coupled hybrid system are proposed.
基金Project supported by the National Natural Science Foundation of China(Grant No.11775116)the Jiangsu Qinglan High-Level Talent Project。
文摘A special transformation is introduced and thereby leads to the N-soliton solution of the(2+1)-dimensional generalized Konopelchenko-Dubrovsky-Kaup-Kupershmidt(KDKK) equation.Then,by employing the long wave limit and imposing complex conjugate constraints to the related solitons,various localized interaction solutions are constructed,including the general M-lumps,T-breathers,and hybrid wave solutions.Dynamical behaviors of these solutions are investigated analytically and graphically.The solutions obtained are very helpful in studying the interaction phenomena of nonlinear localized waves.Therefore,we hope these results can provide some theoretical guidance to the experts in oceanography,atmospheric science,and weather forecasting.
文摘Certain hybrid prototypes of dispersive optical solitons that we are looking for can correspond to new or future behaviors, observable or not, developed or will be developed by optical media that present the cubic-quintic-septic law coupled, with strong dispersions. The equation considered for this purpose is that of non-linear Schrödinger. The solutions are obtained using the Bogning-Djeumen Tchaho-Kofané method extended to the new implicit Bogning’ functions. Some of the obtained solutions show that their existence is due only to the Kerr law nonlinearity presence. Graphical representations plotted have confirmed the hybrid and multi-form character of the obtained dispersive optical solitons. We believe that a good understanding of the hybrid dispersive optical solitons highlighted in the context of this work allows to grasp the physical description of systems whose dynamics are governed by nonlinear Schrödinger equation as studied in this work, allowing thereby a relevant improvement of complex problems encountered in particular in nonliear optaics and in optical fibers.
文摘The well-known Riccati differential equations play a key role in many fields,including problems in protein folding,control and stabilization,stochastic control,and cybersecurity(risk analysis and malware propaga-tion).Quantum computer algorithms have the potential to implement faster approximate solutions to the Riccati equations compared with strictly classical algorithms.While systems with many qubits are still under development,there is significant interest in developing algorithms for near-term quantum computers to determine their accuracy and limitations.In this paper,we propose a hybrid quantum-classical algorithm,the Matrix Riccati Solver(MRS).This approach uses a transformation of variables to turn a set of nonlinear differential equation into a set of approximate linear differential equations(i.e.,second order non-constant coefficients)which can in turn be solved using a version of the Harrow-Hassidim-Lloyd(HHL)quantum algorithm for the case of Hermitian matrices.We implement this approach using the Qiskit language and compute near-term results using a 4 qubit IBM Q System quantum computer.Comparisons with classical results and areas for future research are discussed.
文摘Hybrid grids are used for the solution of 3D turbulence Navier Stokes equations. The prismatic grids are generated near the wall, and the tetrahedron grids are generated in the other field. A Navier Stokes solver using Jamson′s finite volume method is developed. The algebraic Baldwin Lomax turbulence model is adopted. The numerical tests show that the above method is very efficient.
基金supported by the National Natural Science Foundation of China(Grant No.51579034)the Open Fund of the Key Laboratory of Ocean Circulation and Waves,Chinese Academy of Sciences(Grant No.KLOCW1502)
文摘A total variation diminishing-weighted average flux (TVD-WAF)-based hybrid numerical scheme for the enhanced version of nonlinearly dispersive Boussinesq-type equations was developed. The one-dimensional governing equations were rewritten in the conservative form and then discretized on a uniform grid. The finite volume method was used to discretize the flux term while the remaining terms were approximated with the finite difference method. The second-order TVD-WAF method was employed in conjunction with the Harten-Lax-van Leer (HLL) Riemann solver to calculate the numerical flux, and the variables at the cell interface for the local Riemann problem were reconstructed via the fourth- order monotone upstream-centered scheme for conservation laws (MUSCL). The time marching scheme based on the third-order TVD Runge- Kutta method was used to obtain numerical solutions. The model was validated through a series of numerical tests, in which wave breaking and a moving shoreline were treated. The good agreement between the computed results, documented analytical solutions, and experimental data demonstrates the correct discretization of the governing equations and high accuracy of the proposed scheme, and also conforms the advantages of the proposed shock-capturing scheme for the enhanced version of the Boussinesq model, including the convenience in the treatment of wave breaking and moving shorelines and without the need for a numerical filter.
文摘In this study, a revised version of some numerical methods for a class of hybrid integro-differential equations with weakly singular kernels (Abel types) is presented. These equations were developed from a class of integro-differential equations of first kind originating from an aeroelasticity problem. By manipulating the bounds of initial conditions with random variations, this study numerically demonstrated the well-posedness properties of the equations. Finally, an assumption of separating variables, allowed for linear splines to be chosen as a basis and for the differentiation and integration of the integro-differential part to be interchanged;hence, a numerical scheme was constructed.
基金supported by NSFC(No.11971296)National Key Research and Development Program of China(No.2021YFA1003004).
文摘We propose new hybrid Lagrange neural networks called LaNets to predict the numerical solutions of partial differential equations.That is,we embed Lagrange interpolation and small sample learning into deep neural network frameworks.Concretely,we first perform Lagrange interpolation in front of the deep feedforward neural network.The Lagrange basis function has a neat structure and a strong expression ability,which is suitable to be a preprocessing tool for pre-fitting and feature extraction.Second,we introduce small sample learning into training,which is beneficial to guide themodel to be corrected quickly.Taking advantages of the theoretical support of traditional numerical method and the efficient allocation of modern machine learning,LaNets achieve higher predictive accuracy compared to the state-of-the-artwork.The stability and accuracy of the proposed algorithmare demonstrated through a series of classical numerical examples,including one-dimensional Burgers equation,onedimensional carburizing diffusion equations,two-dimensional Helmholtz equation and two-dimensional Burgers equation.Experimental results validate the robustness,effectiveness and flexibility of the proposed algorithm.
基金Foundation item: Supported by the National Natural Science Foundation of China(50608036)
文摘As a boundary-type meshless method, the singular hybrid boundary node method(SHBNM) is based on the modified variational principle and the moving least square(MLS) approximation, so it has the advantages of both boundary element method(BEM) and meshless method. In this paper, the dual reciprocity method(DRM) is combined with SHBNM to solve Poisson equation in which the solution is divided into particular solution and general solution. The general solution is achieved by means of SHBNM, and the particular solution is approximated by using the radial basis function(RBF). Only randomly distributed nodes on the bounding surface of the domain are required and it doesn't need extra equations to compute internal parameters in the domain. The postprocess is very simple. Numerical examples for the solution of Poisson equation show that high convergence rates and high accuracy with a small node number are achievable.
文摘In this paper, we are going to present a class of nonlinear equation solving methods. Steffensen’s method is a simple method for solving a nonlinear equation. By using Steffensen’s method and by combining this method with it, we obtain a new method. It can be said that this method, due to not using the function derivative, would be a good method for solving the nonlinear equation compared to Newton’s method. Finally, we will see that Newton’s method and Steffensen’s hybrid method both have a two-order convergence.
基金The authors would like to thank L.Guillot(CEA/DAM)for insightful discussions and CEA/DAM for partial fnancial support.EB was partially supported by the EPSRC grants EP/P01576X/1 and EP/P012434/1.
文摘We devise hybrid high-order(HHO)methods for the acoustic wave equation in the time domain.We frst consider the second-order formulation in time.Using the Newmark scheme for the temporal discretization,we show that the resulting HHO-Newmark scheme is energy-conservative,and this scheme is also amenable to static condensation at each time step.We then consider the formulation of the acoustic wave equation as a frst-order system together with singly-diagonally implicit and explicit Runge-Kutta(SDIRK and ERK)schemes.HHO-SDIRK schemes are amenable to static condensation at each time step.For HHO-ERK schemes,the use of the mixed-order formulation,where the polynomial degree of the cell unknowns is one order higher than that of the face unknowns,is key to beneft from the explicit structure of the scheme.Numerical results on test cases with analytical solutions show that the methods can deliver optimal convergence rates for smooth solutions of order O(hk+1)in the H1-norm and of order O(h^(k+2))in the L^(2)-norm.Moreover,test cases on wave propagation in heterogeneous media indicate the benefts of using high-order methods.
文摘An L-stable block method based on hybrid second derivative algorithm (BHSDA) is provided by a continuous second derivative method that is defined for all values of the independent variable and applied to parabolic partial differential equations (PDEs). The use of the BHSDA to solve PDEs is facilitated by the method of lines which involves making an approximation to the space derivatives, and hence reducing the problem to that of solving a time-dependent system of first order initial value ordinary differential equations. The stability properties of the method is examined and some numerical results presented.
文摘In this paper, we have successfully extended the Jacobian elliptic function expansion approach to nonlinear differential-difference equations. The Hybrid lattice equation is chosen to illustrate this approach. As a consequence, twelve families of Jacobian elliptic function solutions with different parameters of the Hybrid lattice equation are obtained. When the modulus m → 1 or O, doubly-periodic solutions degenerate to solitonic solutions and trigonometric function solutions, respectively.
文摘In this paper,Chebyshev pseudospectral-finite element schemes are proposed for solving three dimensional vorticity equation.Some approximation results in nonisotropic Sobolev spaces are given.The generalized stability and the convergence are proved strictly.The numerical results show the advantages of this method.The technique in this paper is also applicable to other three-dimensional nonlinear problems in fluid dynamics.