To solve the Poisson equation it is usually possible to discretize it into solving the corresponding linear system Ax=b.Variational quantum algorithms(VQAs)for the discretized Poisson equation have been studied before...To solve the Poisson equation it is usually possible to discretize it into solving the corresponding linear system Ax=b.Variational quantum algorithms(VQAs)for the discretized Poisson equation have been studied before.We present a VQA based on the banded Toeplitz systems for solving the Poisson equation with respect to the structural features of matrix A.In detail,we decompose the matrices A and A^(2)into a linear combination of the corresponding banded Toeplitz matrix and sparse matrices with only a few non-zero elements.For the one-dimensional Poisson equation with different boundary conditions and the d-dimensional Poisson equation with Dirichlet boundary conditions,the number of decomposition terms is less than that reported in[Phys.Rev.A 2023108,032418].Based on the decomposition of the matrix,we design quantum circuits that efficiently evaluate the cost function.Additionally,numerical simulation verifies the feasibility of the proposed algorithm.Finally,the VQAs for linear systems of equations and matrix-vector multiplications with the K-banded Toeplitz matrix T_(n)^(K)are given,where T_(n)^(K)∈R^(n×n)and K∈O(ploylogn).展开更多
In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al...In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.展开更多
We deal with a large solution to the semilinear Poisson equation with doublepower nonlinearityΔ^(u)=u^(p)+αu^(q)in a bounded smooth domain D■R^(n),where p>1,-1<q<p andα∈R.We obtain the asymptotic behavio...We deal with a large solution to the semilinear Poisson equation with doublepower nonlinearityΔ^(u)=u^(p)+αu^(q)in a bounded smooth domain D■R^(n),where p>1,-1<q<p andα∈R.We obtain the asymptotic behavior of a solution u near the boundary OD up to the third or higher term.展开更多
In this paper,the method of polynomial particular solutions is used to solve nonlinear Poisson-type partial differential equations in one,two,and three dimensions.The condition number of the coefficient matrix is redu...In this paper,the method of polynomial particular solutions is used to solve nonlinear Poisson-type partial differential equations in one,two,and three dimensions.The condition number of the coefficient matrix is reduced through the implementation of multiple scale technique,ultimately yielding a stable numerical solution.The methodological process can be divided into two main parts:first,identifying the corresponding polynomial particular solutions for the linear differential operator terms in the governing equations,and second,employing these polynomial particular solutions as basis function to iteratively solve the remaining nonlinear terms within the governing equations.Additionally,we investigate the potential improvement in numerical accuracy for equations with singularities in the analytical solution by shifting the computational domain a certain distance.Numerical experiments are conducted to assess both the accuracy and stability of the proposed method.A comparison of the obtained results with those produced by other numerical methods demonstrates the accuracy,stability,and efficiency of the proposed method in handling nonlinear Poisson-type partial differential equations.展开更多
P-and S-wave separation plays an important role in elastic reverse-time migration.It can reduce the artifacts caused by crosstalk between different modes and improve image quality.In addition,P-and Swave separation ca...P-and S-wave separation plays an important role in elastic reverse-time migration.It can reduce the artifacts caused by crosstalk between different modes and improve image quality.In addition,P-and Swave separation can also be used to better understand and distinguish wave types in complex media.At present,the methods for separating wave modes in anisotropic media mainly include spatial nonstationary filtering,low-rank approximation,and vector Poisson equation.Most of these methods require multiple Fourier transforms or the calculation of large matrices,which require high computational costs for problems with large scale.In this paper,an efficient method is proposed to separate the wave mode for anisotropic media by using a scalar anisotropic Poisson operator in the spatial domain.For 2D problems,the computational complexity required by this method is 1/2 of the methods based on solving a vector Poisson equation.Therefore,compared with existing methods based on pseudoHelmholtz decomposition operators,this method can significantly reduce the computational cost.Numerical examples also show that the P and S waves decomposed by this method not only have the correct amplitude and phase relative to the input wavefield but also can reduce the computational complexity significantly.展开更多
The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, ...The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.展开更多
We study the large-time behavior toward viscous shock waves to the Cauchy problem of the one-dimensional compressible isentropic Navier-Stokes equations with density- dependent viscosity. The nonlinear stability of th...We study the large-time behavior toward viscous shock waves to the Cauchy problem of the one-dimensional compressible isentropic Navier-Stokes equations with density- dependent viscosity. The nonlinear stability of the viscous shock waves is shown for certain class of large initial perturbation with integral zero which can allow the initial density to have large oscillation. Our analysis relies upon the technique developed by Kanel~ and the continuation argument.展开更多
Within the framework of continuum mechanics, the double power series ex- pansion technique is proposed, and a series of reduced one-dimensional （1D） equations for a piezoelectric semiconductor beam are obtained. The...Within the framework of continuum mechanics, the double power series ex- pansion technique is proposed, and a series of reduced one-dimensional （1D） equations for a piezoelectric semiconductor beam are obtained. These derived equations are universal, in which extension, flexure, and shear deformations are all included, and can be degen- erated to a number of special cases, e.g., extensional motion, coupled extensional and flexural motion with shear deformations, and elementary flexural motion without shear deformations. As a typical application, the extensional motion of a ZnO beam is analyzed sequentially. It is revealed that semi-conduction has a great effect on the performance of the piezoelectric semiconductor beam, including static deformations and dynamic be- haviors. A larger initial carrier density will evidently lead to a lower resonant frequency and a smaller displacement response, which is a little similar to the dissipative effect. Both the derived approximate equations and the corresponding qualitative analysis are general and widely applicable, which can clearly interpret the inner physical mechanism of the semiconductor in the piezoelectrics and provide theoretical guidance for further experimental design.展开更多
We study the nonlinear stability of viscous shock waves for the Cauchy problem of one-dimensional nonisentropic compressible Navier–Stokes equations for a viscous and heat conducting ideal polytropic gas. The viscous...We study the nonlinear stability of viscous shock waves for the Cauchy problem of one-dimensional nonisentropic compressible Navier–Stokes equations for a viscous and heat conducting ideal polytropic gas. The viscous shock waves are shown to be time asymptotically stable under large initial perturbation with no restriction on the range of the adiabatic exponent provided that the strengths of the viscous shock waves are assumed to be sufficiently small.The proofs are based on the nonlinear energy estimates and the crucial step is to obtain the positive lower and upper bounds of the density and the temperature which are uniformly in time and space.展开更多
In this paper, we discuss one-dimensional optimal system and the invariant solutions of Coupled Burgers’ equations. By using Wu-differential characteristic set algorithm with the aid of Mathematica software, the clas...In this paper, we discuss one-dimensional optimal system and the invariant solutions of Coupled Burgers’ equations. By using Wu-differential characteristic set algorithm with the aid of Mathematica software, the classical symmetries of the Coupled Burgers’ equations are calculated, and the one-dimensional optimal system of Lie algebra is constructed. And we obtain the invariant solution of the Coupled Burgers’ equations corresponding to one element in one dimensional optimal system by using the invariant method. The results generalize the exact solutions of the Coupled Burgers’ equations.展开更多
A local alternating segment explicit - implicit method for the solution of 2D diffusion equations is presented in this paper .The method is unconditionally stable and has the obvious property of parallelism. Some nume...A local alternating segment explicit - implicit method for the solution of 2D diffusion equations is presented in this paper .The method is unconditionally stable and has the obvious property of parallelism. Some numerical experiments show the method is not only simple but also more accurate.展开更多
Having realized various significant roles that higher-dimensional nonlinear partial differ-ential equations(NLPDEs)play in engineering,we analytically investigate in this paper,a higher-dimensional soliton equation,wi...Having realized various significant roles that higher-dimensional nonlinear partial differ-ential equations(NLPDEs)play in engineering,we analytically investigate in this paper,a higher-dimensional soliton equation,with applications particularly in ocean physics and mechatronics(electrical electronics and mechanical)engineering.Infinitesimal generators of Lie point symmetries of the equation are computed using Lie group analysis of differen-tial equations.In addition,we construct commutation as well as Lie adjoint representation tables for the nine-dimensional Lie algebra achieved.Further,a one-dimensional optimal system of Lie subalgebras is also presented for the soliton equation.This consequently enables us to generate abundant group-invariant solutions through the reduction of the understudy equation into various ordinary differential equations(ODEs).On solving the achieved nonlinear differential equations,we secure various solitonic solutions.In conse-quence,these solutions containing diverse mathematical functions furnish copious shapes of dynamical wave structures,ranging from periodic,kink and kink-shaped nanopteron,soliton(bright and dark)to breather waves with extensive wave collisions depicted.We physically interpreted the resulting soliton solutions by imploring graphical depictions in three dimensions,two dimensions and density plots.Moreover,the gained group-invariant solutions involved several arbitrary functions,thus exhibiting rich physical structures.We also implore the power series technique to solve part of the complicated differential equa-tions and give valid comments on their results.Later,we outline some applications of our results in ocean physics and mechatronics engineering.展开更多
The finite difference method such as alternating group iterative methods is useful in numerical method for evolutionary equations and this is the standard approach taken in this paper. Alternating group explicit (AGE)...The finite difference method such as alternating group iterative methods is useful in numerical method for evolutionary equations and this is the standard approach taken in this paper. Alternating group explicit (AGE) iterative methods for one-dimensional convection diffusion equations problems are given. The stability and convergence are analyzed by the linear method. Numerical results of the model problem are taken. Known test problems have been studied to demonstrate the accuracy of the method. Numerical results show that the behavior of the method with emphasis on treatment of boundary conditions is valuable.展开更多
In this paper, an approximate solution for the one-dimensional hyperbolic telegraph equation by using the q-homotopy analysis method （q-HAM） is proposed.The results shows that the convergence of the q- homotopy anal...In this paper, an approximate solution for the one-dimensional hyperbolic telegraph equation by using the q-homotopy analysis method （q-HAM） is proposed.The results shows that the convergence of the q- homotopy analysis method is more accurate than the convergence of the homotopy analysis method （HAM）.展开更多
This paper presents a Poisson theory of the generalized Birkhoff equations, including the algebraic structure of the equations, the sufficient and necessary condition on the integral and the conditions under which a n...This paper presents a Poisson theory of the generalized Birkhoff equations, including the algebraic structure of the equations, the sufficient and necessary condition on the integral and the conditions under which a new integral can be deduced by a known integral as well as the form of the new integral.展开更多
The modelling of risky asset by stochastic processes with continuous paths, based on Brow- nian motions, suffers from several defects. First, the path continuity assumption does not seem reason- able in view of the po...The modelling of risky asset by stochastic processes with continuous paths, based on Brow- nian motions, suffers from several defects. First, the path continuity assumption does not seem reason- able in view of the possibility of sudden price variations （jumps） resulting of market crashes. A solution is to use stochastic processes with jumps, that will account for sudden variations of the asset prices. On the other hand, such jump models are generally based on the Poisson random measure. Many popular economic and financial models described by stochastic differential equations with Poisson jumps. This paper deals with the approximate controllability of a class of second-order neutral stochastic differential equations with infinite delay and Poisson jumps. By using the cosine family of operators, stochastic analysis techniques, a new set of sufficient conditions are derived for the approximate controllability of the above control system. An example is provided to illustrate the obtained theory.展开更多
A new method for solving the 1D Poisson equation is presented using the finite difference method. This method is based on the exact formulation of the inverse of the tridiagonal matrix associated with the Laplacian. T...A new method for solving the 1D Poisson equation is presented using the finite difference method. This method is based on the exact formulation of the inverse of the tridiagonal matrix associated with the Laplacian. This is the first time that the inverse of this remarkable matrix is determined directly and exactly. Thus, solving 1D Poisson equation becomes very accurate and extremely fast. This method is a very important tool for physics and engineering where the Poisson equation appears very often in the description of certain phenomena.展开更多
The motion of the self-gravitational gaseous stars can be described by the Euler-Poisson equations. The main purpose of this paper is concerned with the existence of stationary solutions of Euler-Poisson equations for...The motion of the self-gravitational gaseous stars can be described by the Euler-Poisson equations. The main purpose of this paper is concerned with the existence of stationary solutions of Euler-Poisson equations for some velocity fields and entropy functions that solve the conservation of mass and energy. Under different restriction to the strength of velocity field, we get the existence and multiplicity of the stationary solutions of Euler-Poisson system.展开更多
In this work, the three dimensional Poisson’s equation in Cartesian coordinates with the Dirichlet’s boundary conditions in a cube is solved directly, by extending the method of Hockney. The Poisson equation is appr...In this work, the three dimensional Poisson’s equation in Cartesian coordinates with the Dirichlet’s boundary conditions in a cube is solved directly, by extending the method of Hockney. The Poisson equation is approximated by 19-points and 27-points fourth order finite difference approximation schemes and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The efficiency of this method is tested for some Poisson’s equations with known analytical solutions and the numerical results obtained show that the method produces accurate results. It is shown that 19-point formula produces comparable results with 27-point formula, though computational efforts are more in 27-point formula.展开更多
基金supported by the Shandong Provincial Natural Science Foundation for Quantum Science under Grant No.ZR2021LLZ002the Fundamental Research Funds for the Central Universities under Grant No.22CX03005A。
文摘To solve the Poisson equation it is usually possible to discretize it into solving the corresponding linear system Ax=b.Variational quantum algorithms(VQAs)for the discretized Poisson equation have been studied before.We present a VQA based on the banded Toeplitz systems for solving the Poisson equation with respect to the structural features of matrix A.In detail,we decompose the matrices A and A^(2)into a linear combination of the corresponding banded Toeplitz matrix and sparse matrices with only a few non-zero elements.For the one-dimensional Poisson equation with different boundary conditions and the d-dimensional Poisson equation with Dirichlet boundary conditions,the number of decomposition terms is less than that reported in[Phys.Rev.A 2023108,032418].Based on the decomposition of the matrix,we design quantum circuits that efficiently evaluate the cost function.Additionally,numerical simulation verifies the feasibility of the proposed algorithm.Finally,the VQAs for linear systems of equations and matrix-vector multiplications with the K-banded Toeplitz matrix T_(n)^(K)are given,where T_(n)^(K)∈R^(n×n)and K∈O(ploylogn).
基金supported by the National Science Foundation grant DMS-1818998.
文摘In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.
基金supported by the JSPS KAKENHI(JP22K03386)supported by the JST SPRING(JPMJSP2132)。
文摘We deal with a large solution to the semilinear Poisson equation with doublepower nonlinearityΔ^(u)=u^(p)+αu^(q)in a bounded smooth domain D■R^(n),where p>1,-1<q<p andα∈R.We obtain the asymptotic behavior of a solution u near the boundary OD up to the third or higher term.
基金supported by the Science and Technology Research Project of Education Department of Jiangxi Province,Philippine Council for Industry,Energy,and Emerging Technology Research and Development(Grant Nos.GJJ2200605)the National Natural Science Foundation of China(11461026)Jiangxi Provincial Natural Science Foundation(20282BA8201016).
文摘In this paper,the method of polynomial particular solutions is used to solve nonlinear Poisson-type partial differential equations in one,two,and three dimensions.The condition number of the coefficient matrix is reduced through the implementation of multiple scale technique,ultimately yielding a stable numerical solution.The methodological process can be divided into two main parts:first,identifying the corresponding polynomial particular solutions for the linear differential operator terms in the governing equations,and second,employing these polynomial particular solutions as basis function to iteratively solve the remaining nonlinear terms within the governing equations.Additionally,we investigate the potential improvement in numerical accuracy for equations with singularities in the analytical solution by shifting the computational domain a certain distance.Numerical experiments are conducted to assess both the accuracy and stability of the proposed method.A comparison of the obtained results with those produced by other numerical methods demonstrates the accuracy,stability,and efficiency of the proposed method in handling nonlinear Poisson-type partial differential equations.
基金supported by the National Key R&D Program of China(No.2018YFA0702505)the project of CNOOC Limited(Grant No.CNOOC-KJ GJHXJSGG YF 2022-01)+1 种基金R&D Department of China National Petroleum Corporation(Investigations on fundamental experiments and advanced theoretical methods in geophysical prospecting application,2022DQ0604-02)NSFC(Grant Nos.U23B20159,41974142,42074129,12001311)。
文摘P-and S-wave separation plays an important role in elastic reverse-time migration.It can reduce the artifacts caused by crosstalk between different modes and improve image quality.In addition,P-and Swave separation can also be used to better understand and distinguish wave types in complex media.At present,the methods for separating wave modes in anisotropic media mainly include spatial nonstationary filtering,low-rank approximation,and vector Poisson equation.Most of these methods require multiple Fourier transforms or the calculation of large matrices,which require high computational costs for problems with large scale.In this paper,an efficient method is proposed to separate the wave mode for anisotropic media by using a scalar anisotropic Poisson operator in the spatial domain.For 2D problems,the computational complexity required by this method is 1/2 of the methods based on solving a vector Poisson equation.Therefore,compared with existing methods based on pseudoHelmholtz decomposition operators,this method can significantly reduce the computational cost.Numerical examples also show that the P and S waves decomposed by this method not only have the correct amplitude and phase relative to the input wavefield but also can reduce the computational complexity significantly.
文摘The solution of Poisson’s Equation plays an important role in many areas, including modeling high-intensity and high-brightness beams in particle accelerators. For the computational domain with a large aspect ratio, the integrated Green’s function method has been adopted to solve the 3D Poisson equation subject to open boundary conditions. In this paper, we report on the efficient implementation of this method, which can save more than a factor of 50 computing time compared with the direct brute force implementation and its improvement under certain extreme conditions.
基金supported by"the Fundamental Research Funds for the Central Universities"
文摘We study the large-time behavior toward viscous shock waves to the Cauchy problem of the one-dimensional compressible isentropic Navier-Stokes equations with density- dependent viscosity. The nonlinear stability of the viscous shock waves is shown for certain class of large initial perturbation with integral zero which can allow the initial density to have large oscillation. Our analysis relies upon the technique developed by Kanel~ and the continuation argument.
基金Project supported by the National Natural Science Foundation of China(Nos.11672223,11402187,and 51178390)the China Postdoctoral Science Foundation(No.2014M560762)the Fundamental Research Funds for the Central Universities of China(No.xjj2015131)
文摘Within the framework of continuum mechanics, the double power series ex- pansion technique is proposed, and a series of reduced one-dimensional （1D） equations for a piezoelectric semiconductor beam are obtained. These derived equations are universal, in which extension, flexure, and shear deformations are all included, and can be degen- erated to a number of special cases, e.g., extensional motion, coupled extensional and flexural motion with shear deformations, and elementary flexural motion without shear deformations. As a typical application, the extensional motion of a ZnO beam is analyzed sequentially. It is revealed that semi-conduction has a great effect on the performance of the piezoelectric semiconductor beam, including static deformations and dynamic be- haviors. A larger initial carrier density will evidently lead to a lower resonant frequency and a smaller displacement response, which is a little similar to the dissipative effect. Both the derived approximate equations and the corresponding qualitative analysis are general and widely applicable, which can clearly interpret the inner physical mechanism of the semiconductor in the piezoelectrics and provide theoretical guidance for further experimental design.
文摘We study the nonlinear stability of viscous shock waves for the Cauchy problem of one-dimensional nonisentropic compressible Navier–Stokes equations for a viscous and heat conducting ideal polytropic gas. The viscous shock waves are shown to be time asymptotically stable under large initial perturbation with no restriction on the range of the adiabatic exponent provided that the strengths of the viscous shock waves are assumed to be sufficiently small.The proofs are based on the nonlinear energy estimates and the crucial step is to obtain the positive lower and upper bounds of the density and the temperature which are uniformly in time and space.
文摘In this paper, we discuss one-dimensional optimal system and the invariant solutions of Coupled Burgers’ equations. By using Wu-differential characteristic set algorithm with the aid of Mathematica software, the classical symmetries of the Coupled Burgers’ equations are calculated, and the one-dimensional optimal system of Lie algebra is constructed. And we obtain the invariant solution of the Coupled Burgers’ equations corresponding to one element in one dimensional optimal system by using the invariant method. The results generalize the exact solutions of the Coupled Burgers’ equations.
文摘A local alternating segment explicit - implicit method for the solution of 2D diffusion equations is presented in this paper .The method is unconditionally stable and has the obvious property of parallelism. Some numerical experiments show the method is not only simple but also more accurate.
基金the North-West University,Mafikeng campus for its continued support.
文摘Having realized various significant roles that higher-dimensional nonlinear partial differ-ential equations(NLPDEs)play in engineering,we analytically investigate in this paper,a higher-dimensional soliton equation,with applications particularly in ocean physics and mechatronics(electrical electronics and mechanical)engineering.Infinitesimal generators of Lie point symmetries of the equation are computed using Lie group analysis of differen-tial equations.In addition,we construct commutation as well as Lie adjoint representation tables for the nine-dimensional Lie algebra achieved.Further,a one-dimensional optimal system of Lie subalgebras is also presented for the soliton equation.This consequently enables us to generate abundant group-invariant solutions through the reduction of the understudy equation into various ordinary differential equations(ODEs).On solving the achieved nonlinear differential equations,we secure various solitonic solutions.In conse-quence,these solutions containing diverse mathematical functions furnish copious shapes of dynamical wave structures,ranging from periodic,kink and kink-shaped nanopteron,soliton(bright and dark)to breather waves with extensive wave collisions depicted.We physically interpreted the resulting soliton solutions by imploring graphical depictions in three dimensions,two dimensions and density plots.Moreover,the gained group-invariant solutions involved several arbitrary functions,thus exhibiting rich physical structures.We also implore the power series technique to solve part of the complicated differential equa-tions and give valid comments on their results.Later,we outline some applications of our results in ocean physics and mechatronics engineering.
文摘The finite difference method such as alternating group iterative methods is useful in numerical method for evolutionary equations and this is the standard approach taken in this paper. Alternating group explicit (AGE) iterative methods for one-dimensional convection diffusion equations problems are given. The stability and convergence are analyzed by the linear method. Numerical results of the model problem are taken. Known test problems have been studied to demonstrate the accuracy of the method. Numerical results show that the behavior of the method with emphasis on treatment of boundary conditions is valuable.
文摘In this paper, an approximate solution for the one-dimensional hyperbolic telegraph equation by using the q-homotopy analysis method （q-HAM） is proposed.The results shows that the convergence of the q- homotopy analysis method is more accurate than the convergence of the homotopy analysis method （HAM）.
基金supported by the National Natural Science Foundation of China (Grant Nos 10572021 and 10772025)
文摘This paper presents a Poisson theory of the generalized Birkhoff equations, including the algebraic structure of the equations, the sufficient and necessary condition on the integral and the conditions under which a new integral can be deduced by a known integral as well as the form of the new integral.
基金supported by the National Board for Higher Mathematics,Mumbai,India under Grant No.2/48(5)/2013/NBHM(R.P.)/RD-II/688 dt 16.01.2014
文摘The modelling of risky asset by stochastic processes with continuous paths, based on Brow- nian motions, suffers from several defects. First, the path continuity assumption does not seem reason- able in view of the possibility of sudden price variations （jumps） resulting of market crashes. A solution is to use stochastic processes with jumps, that will account for sudden variations of the asset prices. On the other hand, such jump models are generally based on the Poisson random measure. Many popular economic and financial models described by stochastic differential equations with Poisson jumps. This paper deals with the approximate controllability of a class of second-order neutral stochastic differential equations with infinite delay and Poisson jumps. By using the cosine family of operators, stochastic analysis techniques, a new set of sufficient conditions are derived for the approximate controllability of the above control system. An example is provided to illustrate the obtained theory.
文摘A new method for solving the 1D Poisson equation is presented using the finite difference method. This method is based on the exact formulation of the inverse of the tridiagonal matrix associated with the Laplacian. This is the first time that the inverse of this remarkable matrix is determined directly and exactly. Thus, solving 1D Poisson equation becomes very accurate and extremely fast. This method is a very important tool for physics and engineering where the Poisson equation appears very often in the description of certain phenomena.
基金supported by NSFC (10631030, 11071094)the fund of CCNU for Ph.D students (2009021)
文摘The motion of the self-gravitational gaseous stars can be described by the Euler-Poisson equations. The main purpose of this paper is concerned with the existence of stationary solutions of Euler-Poisson equations for some velocity fields and entropy functions that solve the conservation of mass and energy. Under different restriction to the strength of velocity field, we get the existence and multiplicity of the stationary solutions of Euler-Poisson system.
文摘In this work, the three dimensional Poisson’s equation in Cartesian coordinates with the Dirichlet’s boundary conditions in a cube is solved directly, by extending the method of Hockney. The Poisson equation is approximated by 19-points and 27-points fourth order finite difference approximation schemes and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The efficiency of this method is tested for some Poisson’s equations with known analytical solutions and the numerical results obtained show that the method produces accurate results. It is shown that 19-point formula produces comparable results with 27-point formula, though computational efforts are more in 27-point formula.
