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).展开更多
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.展开更多
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.展开更多
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.展开更多
A new and innovative method for solving the 1D Poisson Equation is presented, using the finite differences method, with Robin Boundary conditions. The exact formula of the inverse of the discretization matrix is deter...A new and innovative method for solving the 1D Poisson Equation is presented, using the finite differences method, with Robin Boundary conditions. The exact formula of the inverse of the discretization matrix is determined. This is the first time that this famous matrix is inverted explicitly, without using the right hand side. Thus, the solution is determined in a direct, very accurate (O(h2)), and very fast (O(N)) manner. This new approach treats all cases of boundary conditions: Dirichlet, Neumann, and mixed. Therefore, it can serve as a reference for solving the Poisson equation in one dimension.展开更多
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.展开更多
In this paper, the uniqueness of stationary solutions with vacuum of Euler-Poisson equations is considered. Through a nonlinear transformation which is a function of density and entropy, the corresponding problem can ...In this paper, the uniqueness of stationary solutions with vacuum of Euler-Poisson equations is considered. Through a nonlinear transformation which is a function of density and entropy, the corresponding problem can be reduced to a semilinear elliptic equation with a nonlinear source term consisting of a power function, for which the classical theory of the elliptic equations leads the authors to the uniqueness result under some assumptions on the entropy function S(x). As an example, the authors get the uniqueness of stationary solutions with vacuum of Euler-Poisson equations for S(x) =|x|θandθ∈{0}∪[2(N-2),+∞).展开更多
This paper is concerned with the system of Euler-Poisson equations as a model to describe the motion of the self-induced gravitational gaseous stars. When ~ 〈 7 〈 2, under the weak smoothness of entropy function, we...This paper is concerned with the system of Euler-Poisson equations as a model to describe the motion of the self-induced gravitational gaseous stars. When ~ 〈 7 〈 2, under the weak smoothness of entropy function, we find a sufficient condition to guarantee the existence of stationary solutions for some velocity fields and entropy function that solve the conservation of mass and energy.展开更多
In this work, the three-dimensional Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved directly, by extending the method of Hockney...In this work, the three-dimensional Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved directly, by extending the method of Hockney. The Poisson equation is approximated by second-order finite differences and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The accuracy 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.展开更多
In this work, by extending the method of Hockney into three dimensions, the Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved dire...In this work, by extending the method of Hockney into three dimensions, the Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved directly. The Poisson equation is approximated by fourth-order finite differences and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The accuracy 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.展开更多
基金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 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.
基金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.
文摘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.
文摘A new and innovative method for solving the 1D Poisson Equation is presented, using the finite differences method, with Robin Boundary conditions. The exact formula of the inverse of the discretization matrix is determined. This is the first time that this famous matrix is inverted explicitly, without using the right hand side. Thus, the solution is determined in a direct, very accurate (O(h2)), and very fast (O(N)) manner. This new approach treats all cases of boundary conditions: Dirichlet, Neumann, and mixed. Therefore, it can serve as a reference for solving the Poisson equation in one dimension.
文摘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.
基金the Natural Science Foundation of China and the Excellent Teachers Foundation of Ministry of Education of China.
文摘In this paper, the uniqueness of stationary solutions with vacuum of Euler-Poisson equations is considered. Through a nonlinear transformation which is a function of density and entropy, the corresponding problem can be reduced to a semilinear elliptic equation with a nonlinear source term consisting of a power function, for which the classical theory of the elliptic equations leads the authors to the uniqueness result under some assumptions on the entropy function S(x). As an example, the authors get the uniqueness of stationary solutions with vacuum of Euler-Poisson equations for S(x) =|x|θandθ∈{0}∪[2(N-2),+∞).
基金supported by the Fundamental Research Funds for the Central Universities(2011-1a-021)
文摘This paper is concerned with the system of Euler-Poisson equations as a model to describe the motion of the self-induced gravitational gaseous stars. When ~ 〈 7 〈 2, under the weak smoothness of entropy function, we find a sufficient condition to guarantee the existence of stationary solutions for some velocity fields and entropy function that solve the conservation of mass and energy.
文摘In this work, the three-dimensional Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved directly, by extending the method of Hockney. The Poisson equation is approximated by second-order finite differences and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The accuracy 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.
文摘In this work, by extending the method of Hockney into three dimensions, the Poisson’s equation in cylindrical coordinates system with the Dirichlet’s boundary conditions in a portion of a cylinder for is solved directly. The Poisson equation is approximated by fourth-order finite differences and the resulting large algebraic system of linear equations is treated systematically in order to get a block tri-diagonal system. The accuracy 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.
