In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to ...In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.展开更多
Using the solution of general Korteweg-de Vries （KdV） equation, the solutions of the generalized variable coefficient Kadomtsev-Petviashvili （KP） equation are constructed, and then its new solitary wave-like solut...Using the solution of general Korteweg-de Vries （KdV） equation, the solutions of the generalized variable coefficient Kadomtsev-Petviashvili （KP） equation are constructed, and then its new solitary wave-like solution and Jacobi elliptic function solution are obtained.展开更多
In this paper, Daftardar-Gejji and Jafari method is applied to solve fractional heat-like and wave-like models with variable coefficients. The method is proved for a variety of problems in one, two and three dimension...In this paper, Daftardar-Gejji and Jafari method is applied to solve fractional heat-like and wave-like models with variable coefficients. The method is proved for a variety of problems in one, two and three dimensional spaces where analytical approximate solutions are obtained. The examples are presented to show the efficiency and simplicity of this method.展开更多
In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error...In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.展开更多
The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary....The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary.Therefore,this method has the merits of purely meshless,easy to program,high solution accuracy and so on.In this paper,we investigate the effectiveness of the BKM for solving Helmholtz-type problems under various conditions through a series of novel numerical experiments.The results demonstrate that the BKM is efficient and achieves high computational accuracy for problems with smooth or continuous boundary conditions.However,when applied to discontinuous boundary problems,the method exhibits significant numerical instability,potentially leading to substantial deviations in the computed results.Finally,three potential improvement strategies are proposed to mitigate this limitation.展开更多
The main purpose of this paper is to use the Chelyshkov-collocation spectral method for solving nonlinear Quadratic integral equations of Volterra type.The method is based on the approximate solutions in terms of Chel...The main purpose of this paper is to use the Chelyshkov-collocation spectral method for solving nonlinear Quadratic integral equations of Volterra type.The method is based on the approximate solutions in terms of Chelyshkov polynomials with unknown coefficients.The Chelyshkov polynomials and their properties are employed to derive the operational matrices of integral and product.The application of these operational matrices for solving the mentioned problem is explained.The error analysis of the proposed method is investigated.Finally,some numerical examples are provided to demonstrate the efficiency of the method.展开更多
In this paper,the convergence of the split-step theta method for stochastic differential equations is analyzed using stochastic C-stability and stochastic B-consistency.The fact that the numerical scheme,which is both...In this paper,the convergence of the split-step theta method for stochastic differential equations is analyzed using stochastic C-stability and stochastic B-consistency.The fact that the numerical scheme,which is both stochastically C-stable and stochastically B-consistent,is convergent has been proved in a previous paper.In order to analyze the convergence of the split-step theta method(θ∈[1/2,1]),the stochastic C-stability and stochastic B-consistency under the condition of global monotonicity have been researched,and the rate of convergence 1/2 has been explored in this paper.It can be seen that the convergence does not require the drift function should satisfy the linear growth condition whenθ=1/2 Furthermore,the rate of the convergence of the split-step scheme for stochastic differential equations with additive noise has been researched and found to be 1.Finally,an example is given to illustrate the convergence with the theoretical results.展开更多
The(2+1)-dimensional integrable generalization of the Gardner(2DG)equation is solved via the inverse scattering transform method in this paper.A kind of general solution of the equation is obtained by introducing long...The(2+1)-dimensional integrable generalization of the Gardner(2DG)equation is solved via the inverse scattering transform method in this paper.A kind of general solution of the equation is obtained by introducing long derivatives V_(x),V_(y),V_(t).Two different constraints on the kernel function K are introduced under the reality of the solution u of the 2DG equation.Then,two classes of exact solutions with constant asymptotic values at infinity u|x^(2)+y^(2)→∞→0 are constructed by means of the∂¯-dressing method for the casesσ=1 andσ=i.The rational and multiple pole solutions of the 2DG equation are obtained with the kernel functions of zero-order and higher-order Dirac delta functions,respectively.展开更多
In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving th...In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving the convergence of the one-point large deviations rate function (LDRF) of the spatial FDM, which is about the asymptotical limit of a parametric variational problem. The main idea for proving the convergence of the LDRF of the spatial FDM is via the Γ-convergence of objective functions. This relies on the qualitative analysis of skeleton equations of the original equation and the numerical method. In order to overcome the difficulty that the drift coefficient is not one-sided Lipschitz continuous, we derive the equivalent characterization of the skeleton equation of the spatial FDM and the discrete interpolation inequality to obtain the uniform boundedness of the solution to the underlying skeleton equation. These play important roles in deriving the T-convergence of objective functions.展开更多
This paper deals with numerical computation and analysis for the initial boundary problems of two dimensional(2D)Sobolev equations with piecewise continuous argument.Firstly,a two-level high-order compact difference m...This paper deals with numerical computation and analysis for the initial boundary problems of two dimensional(2D)Sobolev equations with piecewise continuous argument.Firstly,a two-level high-order compact difference method(HOCDM)with computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4))is suggested,whereτ,h_(x),h_(y) denote the temporal and spatial stepsizes of the method,respectively.In order to improve the temporal computational accuracy of this method,the Richardson extrapolation technique is used and thus a new two-level HOCDMis derived,which is proved to be convergent of order four both in time and space.Although the new two-level HOCDM has the higher computational accuracy in time than the previous one,it will bring a larger computational cost.To overcome this deficiency,a three-level HOCDM with computational accuracy O(τ^(4)+h_(x)^(4)+h_(y)^(4))is constructed.Finally,with a series of numerical experiments,the theoretical accuracy and computational efficiency of the above methods are further verified.展开更多
The operator splitting method is used to deal with the Navier-Stokes equation, in which the physical process described by the equation is decomposed into two processes： a diffusion process and a convection process; a...The operator splitting method is used to deal with the Navier-Stokes equation, in which the physical process described by the equation is decomposed into two processes： a diffusion process and a convection process; and the finite element equation is established. The velocity field in the element is described by the shape function of the isoparametric element with nine nodes and the pressure field is described by the interpolation function of the four nodes at the vertex of the isoparametric element with nine nodes. The subroutine of the element and the integrated finite element code are generated by the Finite Element Program Generator （FEPG） successfully. The numerical simulation about the incompressible viscous liquid flowing over a cylinder is carded out. The solution agrees with the experimental results very well.展开更多
A new algorithm based on the projection method with the implicit finite difference technique was established to calculate the velocity fields and pressure.The calculation region can be divided into different regions a...A new algorithm based on the projection method with the implicit finite difference technique was established to calculate the velocity fields and pressure.The calculation region can be divided into different regions according to Reynolds number.In the far-wall region,the thermal melt flow was calculated as Newtonian flow.In the near-wall region,the thermal melt flow was calculated as non-Newtonian flow.It was proved that the new algorithm based on the projection method with the implicit technique was correct through nonparametric statistics method and experiment.The simulation results show that the new algorithm based on the projection method with the implicit technique calculates more quickly than the solution algorithm-volume of fluid method using the explicit difference method.展开更多
The preconditioning method is used to solve the low Mach number flow. The space discritisation scheme is the Roe scheme and the DES turbulence model is used. Then, the low Mach number turbulence flow around the NACA00...The preconditioning method is used to solve the low Mach number flow. The space discritisation scheme is the Roe scheme and the DES turbulence model is used. Then, the low Mach number turbulence flow around the NACA0012 airfoil is used to verify the efficiency of the proposed method. Two cases of the low Mach number flows around the multi-element airfoil and the circular cylinder are also used to test the proposed method. Numerical results show that the methods combined the preconditioning method and compressible Navier-Stokes equations are efficient to solve low Mach number flows.展开更多
The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data...The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.展开更多
In this paper, a new fractional projective Riccati equation method is proposed to establish exact solutions for fractional partial differential equations in the sense of modified Riemann–Liouville derivative. This me...In this paper, a new fractional projective Riccati equation method is proposed to establish exact solutions for fractional partial differential equations in the sense of modified Riemann–Liouville derivative. This method can be seen as the fractional version of the known projective Riccati equation method. For illustrating the validity of this method,we apply this method to solve the space-time fractional Whitham–Broer–Kaup(WBK) equations and the nonlinear fractional Sharma–Tasso–Olever(STO) equation, and as a result, some new exact solutions for them are obtained.展开更多
In this article, a novel （G＇/G）-expansion method is proposed to search for the traveling wave solutions of nonlinear evolution equations. We construct abundant traveling wave solutions involving parameters to the B...In this article, a novel （G＇/G）-expansion method is proposed to search for the traveling wave solutions of nonlinear evolution equations. We construct abundant traveling wave solutions involving parameters to the Boussinesq equation by means of the suggested method. The performance of the method is reliable and useful, and gives more general exact solutions than the existing methods. The new （G＇/G）-expansion method provides not only more general forms of solutions but also cuspon, peakon, soliton, and periodic waves.展开更多
An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same t...An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.展开更多
The dispersed phase in multiphase flows can be modeled by the population balance model(PBM). A typical population balance equation(PBE) contains terms for spatial transport, loss/growth and breakage/coalescence source...The dispersed phase in multiphase flows can be modeled by the population balance model(PBM). A typical population balance equation(PBE) contains terms for spatial transport, loss/growth and breakage/coalescence source terms. The equation is therefore quite complex and difficult to solve analytically or numerically. The quadrature-based moment methods(QBMMs) are a class of methods that solve the PBE by converting the transport equation of the number density function(NDF) into moment transport equations. The unknown source terms are closed by numerical quadrature. Over the years, many QBMMs have been developed for different problems, such as the quadrature method of moments(QMOM), direct quadrature method of moments(DQMOM),extended quadrature method of moments(EQMOM), conditional quadrature method of moments(CQMOM),extended conditional quadrature method of moments(ECQMOM) and hyperbolic quadrature method of moments(Hy QMOM). In this paper, we present a comprehensive algorithm review of these QBMMs. The mathematical equations for spatially homogeneous systems with first-order point processes and second-order point processes are derived in detail. The algorithms are further extended to the inhomogeneous system for multiphase flows, in which the computational fluid dynamics(CFD) can be coupled with the PBE. The physical limitations and the challenging numerical problems of these QBMMs are discussed. Possible solutions are also summarized.展开更多
In this paper,the parametric equations with multipliers of nonholonomic nonconservative sys- tems in the event space are established,their properties are studied,and their explicit formulation is obtained. And then th...In this paper,the parametric equations with multipliers of nonholonomic nonconservative sys- tems in the event space are established,their properties are studied,and their explicit formulation is obtained. And then the field method for integrating these equations is given.Finally,an example illustrating the appli- cation of the integration method is given.展开更多
An extended Fan's algebraic method is used for constructing exact traveling wave solution of nonlinearpartial differential equations.The key idea of this method is to introduce an auxiliary ordinary differential e...An extended Fan's algebraic method is used for constructing exact traveling wave solution of nonlinearpartial differential equations.The key idea of this method is to introduce an auxiliary ordinary differential equationwhich is regarded as an extended elliptic equation and whose degree Υ is expanded to the case of r>4.The efficiency ofthe method is demonstrated by the KdV equation and the variant Boussinesq equations.The results indicate that themethod not only offers all solutions obtained by using Fu's and Fan's methods,but also some new solutions.展开更多
基金Supported in part by Natural Science Foundation of Guangxi(2023GXNSFAA026246)in part by the Central Government's Guide to Local Science and Technology Development Fund(GuikeZY23055044)in part by the National Natural Science Foundation of China(62363003)。
文摘In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.
文摘Using the solution of general Korteweg-de Vries （KdV） equation, the solutions of the generalized variable coefficient Kadomtsev-Petviashvili （KP） equation are constructed, and then its new solitary wave-like solution and Jacobi elliptic function solution are obtained.
文摘In this paper, Daftardar-Gejji and Jafari method is applied to solve fractional heat-like and wave-like models with variable coefficients. The method is proved for a variety of problems in one, two and three dimensional spaces where analytical approximate solutions are obtained. The examples are presented to show the efficiency and simplicity of this method.
文摘In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.
基金Supported by the Key Scientific Research Plan of Colleges and Universities in Henan Province(23B140006)。
文摘The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary.Therefore,this method has the merits of purely meshless,easy to program,high solution accuracy and so on.In this paper,we investigate the effectiveness of the BKM for solving Helmholtz-type problems under various conditions through a series of novel numerical experiments.The results demonstrate that the BKM is efficient and achieves high computational accuracy for problems with smooth or continuous boundary conditions.However,when applied to discontinuous boundary problems,the method exhibits significant numerical instability,potentially leading to substantial deviations in the computed results.Finally,three potential improvement strategies are proposed to mitigate this limitation.
文摘The main purpose of this paper is to use the Chelyshkov-collocation spectral method for solving nonlinear Quadratic integral equations of Volterra type.The method is based on the approximate solutions in terms of Chelyshkov polynomials with unknown coefficients.The Chelyshkov polynomials and their properties are employed to derive the operational matrices of integral and product.The application of these operational matrices for solving the mentioned problem is explained.The error analysis of the proposed method is investigated.Finally,some numerical examples are provided to demonstrate the efficiency of the method.
基金Supported by the National Natural Science Foundation of China (Grant No. 12301521)the Natural Science Foundation of Shanxi Province (Grant No. 20210302124081)。
文摘In this paper,the convergence of the split-step theta method for stochastic differential equations is analyzed using stochastic C-stability and stochastic B-consistency.The fact that the numerical scheme,which is both stochastically C-stable and stochastically B-consistent,is convergent has been proved in a previous paper.In order to analyze the convergence of the split-step theta method(θ∈[1/2,1]),the stochastic C-stability and stochastic B-consistency under the condition of global monotonicity have been researched,and the rate of convergence 1/2 has been explored in this paper.It can be seen that the convergence does not require the drift function should satisfy the linear growth condition whenθ=1/2 Furthermore,the rate of the convergence of the split-step scheme for stochastic differential equations with additive noise has been researched and found to be 1.Finally,an example is given to illustrate the convergence with the theoretical results.
基金Supported by the National Natural Science Foundation of China(Grant Nos.1237125611971475)。
文摘The(2+1)-dimensional integrable generalization of the Gardner(2DG)equation is solved via the inverse scattering transform method in this paper.A kind of general solution of the equation is obtained by introducing long derivatives V_(x),V_(y),V_(t).Two different constraints on the kernel function K are introduced under the reality of the solution u of the 2DG equation.Then,two classes of exact solutions with constant asymptotic values at infinity u|x^(2)+y^(2)→∞→0 are constructed by means of the∂¯-dressing method for the casesσ=1 andσ=i.The rational and multiple pole solutions of the 2DG equation are obtained with the kernel functions of zero-order and higher-order Dirac delta functions,respectively.
基金supported by the National Natural Science Foundation of China(12201228,12171047)the Fundamental Research Funds for the Central Universities(3034011102)supported by National Key R&D Program of China(2020YFA0713701).
文摘In this work, we first derive the one-point large deviations principle (LDP) for both the stochastic Cahn–Hilliard equation with small noise and its spatial finite difference method (FDM). Then, we focus on giving the convergence of the one-point large deviations rate function (LDRF) of the spatial FDM, which is about the asymptotical limit of a parametric variational problem. The main idea for proving the convergence of the LDRF of the spatial FDM is via the Γ-convergence of objective functions. This relies on the qualitative analysis of skeleton equations of the original equation and the numerical method. In order to overcome the difficulty that the drift coefficient is not one-sided Lipschitz continuous, we derive the equivalent characterization of the skeleton equation of the spatial FDM and the discrete interpolation inequality to obtain the uniform boundedness of the solution to the underlying skeleton equation. These play important roles in deriving the T-convergence of objective functions.
文摘This paper deals with numerical computation and analysis for the initial boundary problems of two dimensional(2D)Sobolev equations with piecewise continuous argument.Firstly,a two-level high-order compact difference method(HOCDM)with computational accuracy O(τ^(2)+h_(x)^(4)+h_(y)^(4))is suggested,whereτ,h_(x),h_(y) denote the temporal and spatial stepsizes of the method,respectively.In order to improve the temporal computational accuracy of this method,the Richardson extrapolation technique is used and thus a new two-level HOCDMis derived,which is proved to be convergent of order four both in time and space.Although the new two-level HOCDM has the higher computational accuracy in time than the previous one,it will bring a larger computational cost.To overcome this deficiency,a three-level HOCDM with computational accuracy O(τ^(4)+h_(x)^(4)+h_(y)^(4))is constructed.Finally,with a series of numerical experiments,the theoretical accuracy and computational efficiency of the above methods are further verified.
文摘The operator splitting method is used to deal with the Navier-Stokes equation, in which the physical process described by the equation is decomposed into two processes： a diffusion process and a convection process; and the finite element equation is established. The velocity field in the element is described by the shape function of the isoparametric element with nine nodes and the pressure field is described by the interpolation function of the four nodes at the vertex of the isoparametric element with nine nodes. The subroutine of the element and the integrated finite element code are generated by the Finite Element Program Generator （FEPG） successfully. The numerical simulation about the incompressible viscous liquid flowing over a cylinder is carded out. The solution agrees with the experimental results very well.
基金Project (50975263) supported by the National Natural Science Foundation of ChinaProject (2010081015) supported by International Cooperation Project of Shanxi Province, China+1 种基金 Project (2010-78) supported by the Scholarship Council in Shanxi province, ChinaProject (2010420120005) supported by Doctoral Fund of Ministry of Education of China
文摘A new algorithm based on the projection method with the implicit finite difference technique was established to calculate the velocity fields and pressure.The calculation region can be divided into different regions according to Reynolds number.In the far-wall region,the thermal melt flow was calculated as Newtonian flow.In the near-wall region,the thermal melt flow was calculated as non-Newtonian flow.It was proved that the new algorithm based on the projection method with the implicit technique was correct through nonparametric statistics method and experiment.The simulation results show that the new algorithm based on the projection method with the implicit technique calculates more quickly than the solution algorithm-volume of fluid method using the explicit difference method.
文摘The preconditioning method is used to solve the low Mach number flow. The space discritisation scheme is the Roe scheme and the DES turbulence model is used. Then, the low Mach number turbulence flow around the NACA0012 airfoil is used to verify the efficiency of the proposed method. Two cases of the low Mach number flows around the multi-element airfoil and the circular cylinder are also used to test the proposed method. Numerical results show that the methods combined the preconditioning method and compressible Navier-Stokes equations are efficient to solve low Mach number flows.
基金supported by the National Science and Technology Major Project of China(Grant No. 2011ZX05004-003,2011ZX05014-006-006)the National Key Basic Research Program of China(Grant No. 2013CB228602)the Natural Science Foundation of China(Grant No. 40974066)
文摘The Finite Difference （FD） method is an important method for seismic numerical simulations. It helps us understand regular patterns in seismic wave propagation, analyze seismic attributes, and interpret seismic data. However, because of its discretization, the FD method is only stable under certain conditions. The Arbitrary Difference Precise Integration （ADPI） method is based on the FD method and adopts an integration scheme in the time domain and an arbitrary difference scheme in the space domain. Therefore, the ADPI method is a semi-analytical method. In this paper, we deduce the formula for the ADPI method based on the 3D elastic equation and improve its stability. In forward modeling cases, the ADPI method was implemented in 2D and 3D elastic wave equation forward modeling. Results show that the travel time of the reflected seismic wave is accurate. Compared with the acoustic wave field, the elastic wave field contains more wave types, including PS- and PP- reflected waves, transmitted waves, and diffracted waves, which is important to interpretation of seismic data. The method can be easily applied to elastic wave equation numerical simulations for eoloical models.
基金Supported by Natural Science Foundation of Shandong Province of China under Grant No.ZR2013AQ009National Training Programs of Innovation and Entrepreneurship for Undergraduates under Grant No.201310433031Doctoral initializing Foundation of Shandong University of Technology of China under Grant No.4041-413030
文摘In this paper, a new fractional projective Riccati equation method is proposed to establish exact solutions for fractional partial differential equations in the sense of modified Riemann–Liouville derivative. This method can be seen as the fractional version of the known projective Riccati equation method. For illustrating the validity of this method,we apply this method to solve the space-time fractional Whitham–Broer–Kaup(WBK) equations and the nonlinear fractional Sharma–Tasso–Olever(STO) equation, and as a result, some new exact solutions for them are obtained.
文摘In this article, a novel （G＇/G）-expansion method is proposed to search for the traveling wave solutions of nonlinear evolution equations. We construct abundant traveling wave solutions involving parameters to the Boussinesq equation by means of the suggested method. The performance of the method is reliable and useful, and gives more general exact solutions than the existing methods. The new （G＇/G）-expansion method provides not only more general forms of solutions but also cuspon, peakon, soliton, and periodic waves.
基金Supported by the National Natural Science Foundation of China (10601022)Natural Science Foundation of Inner Mongolia Autonomous Region (200607010106)Youth Science Foundation of Inner Mongolia University(ND0702)
文摘An H^1-Galerkin mixed finite element method is discussed for a class of second order SchrSdinger equation. Optimal error estimates of semidiscrete schemes are derived for problems in one space dimension. At the same time, optimal error estimates are derived for fully discrete schemes. And it is showed that the H1-Galerkin mixed finite element approximations have the same rate of convergence as in the classical mixed finite element methods without requiring the LBB consistency condition.
文摘The dispersed phase in multiphase flows can be modeled by the population balance model(PBM). A typical population balance equation(PBE) contains terms for spatial transport, loss/growth and breakage/coalescence source terms. The equation is therefore quite complex and difficult to solve analytically or numerically. The quadrature-based moment methods(QBMMs) are a class of methods that solve the PBE by converting the transport equation of the number density function(NDF) into moment transport equations. The unknown source terms are closed by numerical quadrature. Over the years, many QBMMs have been developed for different problems, such as the quadrature method of moments(QMOM), direct quadrature method of moments(DQMOM),extended quadrature method of moments(EQMOM), conditional quadrature method of moments(CQMOM),extended conditional quadrature method of moments(ECQMOM) and hyperbolic quadrature method of moments(Hy QMOM). In this paper, we present a comprehensive algorithm review of these QBMMs. The mathematical equations for spatially homogeneous systems with first-order point processes and second-order point processes are derived in detail. The algorithms are further extended to the inhomogeneous system for multiphase flows, in which the computational fluid dynamics(CFD) can be coupled with the PBE. The physical limitations and the challenging numerical problems of these QBMMs are discussed. Possible solutions are also summarized.
基金The Project is supported by the National Natural Science Foundation of China
文摘In this paper,the parametric equations with multipliers of nonholonomic nonconservative sys- tems in the event space are established,their properties are studied,and their explicit formulation is obtained. And then the field method for integrating these equations is given.Finally,an example illustrating the appli- cation of the integration method is given.
基金National Natural Science Foundation of China under Grant No.10672053
文摘An extended Fan's algebraic method is used for constructing exact traveling wave solution of nonlinearpartial differential equations.The key idea of this method is to introduce an auxiliary ordinary differential equationwhich is regarded as an extended elliptic equation and whose degree Υ is expanded to the case of r>4.The efficiency ofthe method is demonstrated by the KdV equation and the variant Boussinesq equations.The results indicate that themethod not only offers all solutions obtained by using Fu's and Fan's methods,but also some new solutions.
