Hessian matrices are square matrices consisting of all possible combinations of second partial derivatives of a scalar-valued initial function. As such, Hessian matrices may be treated as elementary matrix systems of ...Hessian matrices are square matrices consisting of all possible combinations of second partial derivatives of a scalar-valued initial function. As such, Hessian matrices may be treated as elementary matrix systems of linear second-order partial differential equations. This paper discusses the Hessian and its applications in optimization, and then proceeds to introduce and derive the notion of the Jaffa Transform, a new linear operator that directly maps a Hessian square matrix space to the initial corresponding scalar field in nth dimensional Euclidean space. The Jaffa Transform is examined, including the properties of the operator, the transform of notable matrices, and the existence of an inverse Jaffa Transform, which is, by definition, the Hessian matrix operator. The Laplace equation is then noted and investigated, particularly, the relation of the Laplace equation to Poisson’s equation, and the theoretical applications and correlations of harmonic functions to Hessian matrices. The paper concludes by introducing and explicating the Jaffa Theorem, a principle that declares the existence of harmonic Jaffa Transforms, which are, essentially, Jaffa Transform solutions to the Laplace partial differential equation.展开更多
Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to sol...Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to solutions in the Laplace domain that are not readily invertible to the real domain by analyticalmeans.Thus,we need numerical inversionmethods to convert the obtained solution fromLaplace domain to a real domain.In this paper,we propose a numerical scheme based on Laplace transform and numerical inverse Laplace transform for the approximate solution of fractal-fractional differential equations with orderα,β.Our proposed numerical scheme is based on three main steps.First,we convert the given fractal-fractional differential equation to fractional-differential equation in Riemann-Liouville sense,and then into Caputo sense.Secondly,we transformthe fractional differential equation in Caputo sense to an equivalent equation in Laplace space.Then the solution of the transformed equation is obtained in Laplace domain.Finally,the solution is converted into the real domain using numerical inversion of Laplace transform.Three inversion methods are evaluated in this paper,and their convergence is also discussed.Three test problems are used to validate the inversion methods.We demonstrate our results with the help of tables and figures.The obtained results show that Euler’s and Talbot’s methods performed better than Stehfest’s method.展开更多
The Laplace transformation is a very important integral transform,and it is extensively used in solving ordinary differential equations,partial differential equations,and several types of integro-differential equation...The Laplace transformation is a very important integral transform,and it is extensively used in solving ordinary differential equations,partial differential equations,and several types of integro-differential equations.Our purpose in this study is to introduce the notion of fuzzy double Laplace transform,fuzzy conformable double Laplace transform(FCDLT).We discuss some basic properties of FCDLT.We obtain the solutions of fuzzy partial differential equations(both one-dimensional and two-dimensional cases)through the double Laplace approach.We demonstrate through numerical examples that our proposed method is very successful and convenient for resolving partial differential equations.展开更多
Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridyna...Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.展开更多
In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic ...In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.展开更多
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.展开更多
Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision mod...Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision model of the momentum exchange between the differential propellant mass element (dm) and the rocket final mass (m1), in which dm initially travels forward to collide with m1 and rebounds to exit through the exhaust nozzle with a velocity that is known as the effective exhaust velocity ve. We observe that such a model does not explain how dm was able to acquire its initial forward velocity without the support of a reactive mass traveling in the opposite direction. We show instead that the initial kinetic energy of dm is generated from dm itself by a process of self-combustion and expansion. In our ideal rocket with a single particle dm confined inside a hollow tube with one closed end, we show that the process of self-combustion and expansion of dm will result in a pair of differential particles each with a mass dm/2, and each traveling away from one another along the tube axis, from the center of combustion. These two identical particles represent the active and reactive sub-components of dm, co-generated in compliance with Newton’s third law of equal action and reaction. Building on this model, we derive a linear momentum ODE of the system, the solution of which yields what we call the Revised Tsiolkovsky Rocket Equation (RTRE). We show that RTRE has a mathematical form that is similar to TRE, with the exception of the effective exhaust velocity (ve) term. The ve term in TRE is replaced in RTRE by the average of two distinct exhaust velocities that we refer to as fast-jet, vx<sub>1</sub>, and slow-jet, vx<sub>2</sub>. These two velocities correspond, respectively, to the velocities of the detonation pressure wave that is vectored directly towards the exhaust nozzle, and the retonation wave that is initially vectored in the direction of rocket propagation, but subsequently becomes reflected from the thrust surface of the combustion chamber to exit through the exhaust nozzle with a time lag behind the detonation wave. The detonation-retonation phenomenon is supported by experimental evidence in the published literature. Finally, we use a convolution model to simulate the composite exhaust pressure wave, highlighting the frequency spectrum of the pressure perturbations that are generated by the mutual interference between the fast-jet and slow-jet components. Our analysis offers insights into the origin of combustion oscillations in rocket engines, with possible extensions beyond rocket engineering into other fields of combustion engineering.展开更多
We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filament...We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma. We examined the performance of the applied scheme, in this context, we implemented the developed model to study selected phenomena in terahertz radiation production, such as the excitation energy and conversion efficiency of the produced THz radiation, in addition to the influence of the pulse chirping on properties of the produced radiation. The obtained numerical results have clarified that the applied HO-FDTD scheme is precisely accurate to solve Maxwell’s equations and sufficiently valid to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma.展开更多
In this paper, the modification of double Laplace decomposition method is pro- posed for the analytical approximation solution of a coupled system of pseudo-parabolic equation with initial conditions. Some examples ar...In this paper, the modification of double Laplace decomposition method is pro- posed for the analytical approximation solution of a coupled system of pseudo-parabolic equation with initial conditions. Some examples are given to support our presented method. In addition, we prove the convergence of double Laplace transform decomposition method applied to our problems.展开更多
We discuss the solution of Laplace’s differential equation by using operational calculus in the framework of distribution theory. We here study the solution of that differential Equation with an inhomogeneous term, a...We discuss the solution of Laplace’s differential equation by using operational calculus in the framework of distribution theory. We here study the solution of that differential Equation with an inhomogeneous term, and also a fractional differential equation of the type of Laplace’s differential equation.展开更多
The aim of this paper is to discuss application of Laplace Decomposition Method with Adomian Decomposition in time-space Fractional Nonlinear Fractional Differential Equations. The approximate solutions result from La...The aim of this paper is to discuss application of Laplace Decomposition Method with Adomian Decomposition in time-space Fractional Nonlinear Fractional Differential Equations. The approximate solutions result from Laplace Decomposition Method and Adomian decomposition;those two accessions are comfortable to perform and firm when to PDEs. For caption and further representation of the thought, several examples are tool up.展开更多
In this paper, the Combined Laplace Transform-Adomian Decomposition Method is used to solve nth-order integro-differential equations. The results show that the method is very simple and effective.
Let Ω be a smooth bounded domain in R^n. In this article, we consider the homogeneous boundary Dirichlet problem of inhomogeneous p-Laplace equation --△pu = |u|^q-1 u + λf(x) on Ω, and identify necessary and ...Let Ω be a smooth bounded domain in R^n. In this article, we consider the homogeneous boundary Dirichlet problem of inhomogeneous p-Laplace equation --△pu = |u|^q-1 u + λf(x) on Ω, and identify necessary and sufficient conditions on Ω and f(x) which ensure the existence, or multiplicities of nonnegative solutions for the problem under consideration.展开更多
In this paper, the fluid flow differential equation based on the homogenous reservoirs model is first reviewed. Then a theorem about the formal similarity of solutions in the Laplace space with outer boundary conditio...In this paper, the fluid flow differential equation based on the homogenous reservoirs model is first reviewed. Then a theorem about the formal similarity of solutions in the Laplace space with outer boundary conditions and inner boundary condition is presented and proved. Lastly, a corollary of our theorem is given particularly on inner boundary. The obtained results are very helpful for understanding inherent laws of relevant engineering science and designing practical analysis software.展开更多
In this paper, we deal with the following problem:By variational method, we prove the existenceof a nontrivial weak solution whenand the existence of a cylindricalweak solution when
Structure of nonnegative nontrivial and positive solutions was precisely studied for some singularly perturbed p-Laplace equations. By virtue of sub- and supersolution method, it is shown that there are many nonnegati...Structure of nonnegative nontrivial and positive solutions was precisely studied for some singularly perturbed p-Laplace equations. By virtue of sub- and supersolution method, it is shown that there are many nonnegative nontrivial spike-layer solutions and positive intermediate spike-layer solutions. Moreover, the upper and lower bound on the measure of each spike-layer were estimated when the parameter is sufficiently small.展开更多
We discuss the solution of Laplace’s differential equation and a fractional differential equation of that type, by using analytic continuations of Riemann-Liouville fractional derivative and of Laplace transform. We ...We discuss the solution of Laplace’s differential equation and a fractional differential equation of that type, by using analytic continuations of Riemann-Liouville fractional derivative and of Laplace transform. We show that the solutions, which are obtained by using operational calculus in the framework of distribution theory in our preceding papers, are obtained also by the present method.展开更多
The high-order boundary conditions for the problems cf Laplace equation in infinite region have been developed. The improvement in accuracy for numerical solution is achieved by imposing the high-order boundary condit...The high-order boundary conditions for the problems cf Laplace equation in infinite region have been developed. The improvement in accuracy for numerical solution is achieved by imposing the high-order boundary conditions on the exterior boundarv of a reduced finite region in which the numerical method is used. So both the computing efforts and the required storage in computer are reduced. The numerical examples show that the 1st-order boundary condition approaches to the exact boundary condition and it is clearly superior to the traditional boundary condition and the 2nd-order boundary condition.展开更多
文摘Hessian matrices are square matrices consisting of all possible combinations of second partial derivatives of a scalar-valued initial function. As such, Hessian matrices may be treated as elementary matrix systems of linear second-order partial differential equations. This paper discusses the Hessian and its applications in optimization, and then proceeds to introduce and derive the notion of the Jaffa Transform, a new linear operator that directly maps a Hessian square matrix space to the initial corresponding scalar field in nth dimensional Euclidean space. The Jaffa Transform is examined, including the properties of the operator, the transform of notable matrices, and the existence of an inverse Jaffa Transform, which is, by definition, the Hessian matrix operator. The Laplace equation is then noted and investigated, particularly, the relation of the Laplace equation to Poisson’s equation, and the theoretical applications and correlations of harmonic functions to Hessian matrices. The paper concludes by introducing and explicating the Jaffa Theorem, a principle that declares the existence of harmonic Jaffa Transforms, which are, essentially, Jaffa Transform solutions to the Laplace partial differential equation.
文摘Laplace transform is one of the powerful tools for solving differential equations in engineering and other science subjects.Using the Laplace transform for solving differential equations,however,sometimes leads to solutions in the Laplace domain that are not readily invertible to the real domain by analyticalmeans.Thus,we need numerical inversionmethods to convert the obtained solution fromLaplace domain to a real domain.In this paper,we propose a numerical scheme based on Laplace transform and numerical inverse Laplace transform for the approximate solution of fractal-fractional differential equations with orderα,β.Our proposed numerical scheme is based on three main steps.First,we convert the given fractal-fractional differential equation to fractional-differential equation in Riemann-Liouville sense,and then into Caputo sense.Secondly,we transformthe fractional differential equation in Caputo sense to an equivalent equation in Laplace space.Then the solution of the transformed equation is obtained in Laplace domain.Finally,the solution is converted into the real domain using numerical inversion of Laplace transform.Three inversion methods are evaluated in this paper,and their convergence is also discussed.Three test problems are used to validate the inversion methods.We demonstrate our results with the help of tables and figures.The obtained results show that Euler’s and Talbot’s methods performed better than Stehfest’s method.
基金Manar A.Alqudah would like to thank Princess Nourah bint Abdulrahman University Researchers Supporting Project No.(PNURSP2022R14),Princess Nourah bint Abdulrahman University,Riyadh,Saudi Arabia。
文摘The Laplace transformation is a very important integral transform,and it is extensively used in solving ordinary differential equations,partial differential equations,and several types of integro-differential equations.Our purpose in this study is to introduce the notion of fuzzy double Laplace transform,fuzzy conformable double Laplace transform(FCDLT).We discuss some basic properties of FCDLT.We obtain the solutions of fuzzy partial differential equations(both one-dimensional and two-dimensional cases)through the double Laplace approach.We demonstrate through numerical examples that our proposed method is very successful and convenient for resolving partial differential equations.
文摘Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.
基金supported by China Postdoctoral Science Foundation grant 2020TQ0344the NSFC grants 11871139 and 12101597the NSF grants DMS-1720116,DMS-2012882,DMS-2011838,DMS-1719942,DMS-1913072.
文摘In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.
文摘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.
文摘Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision model of the momentum exchange between the differential propellant mass element (dm) and the rocket final mass (m1), in which dm initially travels forward to collide with m1 and rebounds to exit through the exhaust nozzle with a velocity that is known as the effective exhaust velocity ve. We observe that such a model does not explain how dm was able to acquire its initial forward velocity without the support of a reactive mass traveling in the opposite direction. We show instead that the initial kinetic energy of dm is generated from dm itself by a process of self-combustion and expansion. In our ideal rocket with a single particle dm confined inside a hollow tube with one closed end, we show that the process of self-combustion and expansion of dm will result in a pair of differential particles each with a mass dm/2, and each traveling away from one another along the tube axis, from the center of combustion. These two identical particles represent the active and reactive sub-components of dm, co-generated in compliance with Newton’s third law of equal action and reaction. Building on this model, we derive a linear momentum ODE of the system, the solution of which yields what we call the Revised Tsiolkovsky Rocket Equation (RTRE). We show that RTRE has a mathematical form that is similar to TRE, with the exception of the effective exhaust velocity (ve) term. The ve term in TRE is replaced in RTRE by the average of two distinct exhaust velocities that we refer to as fast-jet, vx<sub>1</sub>, and slow-jet, vx<sub>2</sub>. These two velocities correspond, respectively, to the velocities of the detonation pressure wave that is vectored directly towards the exhaust nozzle, and the retonation wave that is initially vectored in the direction of rocket propagation, but subsequently becomes reflected from the thrust surface of the combustion chamber to exit through the exhaust nozzle with a time lag behind the detonation wave. The detonation-retonation phenomenon is supported by experimental evidence in the published literature. Finally, we use a convolution model to simulate the composite exhaust pressure wave, highlighting the frequency spectrum of the pressure perturbations that are generated by the mutual interference between the fast-jet and slow-jet components. Our analysis offers insights into the origin of combustion oscillations in rocket engines, with possible extensions beyond rocket engineering into other fields of combustion engineering.
文摘We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma. We examined the performance of the applied scheme, in this context, we implemented the developed model to study selected phenomena in terahertz radiation production, such as the excitation energy and conversion efficiency of the produced THz radiation, in addition to the influence of the pulse chirping on properties of the produced radiation. The obtained numerical results have clarified that the applied HO-FDTD scheme is precisely accurate to solve Maxwell’s equations and sufficiently valid to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma.
文摘In this paper, the modification of double Laplace decomposition method is pro- posed for the analytical approximation solution of a coupled system of pseudo-parabolic equation with initial conditions. Some examples are given to support our presented method. In addition, we prove the convergence of double Laplace transform decomposition method applied to our problems.
文摘We discuss the solution of Laplace’s differential equation by using operational calculus in the framework of distribution theory. We here study the solution of that differential Equation with an inhomogeneous term, and also a fractional differential equation of the type of Laplace’s differential equation.
文摘The aim of this paper is to discuss application of Laplace Decomposition Method with Adomian Decomposition in time-space Fractional Nonlinear Fractional Differential Equations. The approximate solutions result from Laplace Decomposition Method and Adomian decomposition;those two accessions are comfortable to perform and firm when to PDEs. For caption and further representation of the thought, several examples are tool up.
文摘In this paper, the Combined Laplace Transform-Adomian Decomposition Method is used to solve nth-order integro-differential equations. The results show that the method is very simple and effective.
基金This work is supported by NNSF of China (10171029).
文摘Let Ω be a smooth bounded domain in R^n. In this article, we consider the homogeneous boundary Dirichlet problem of inhomogeneous p-Laplace equation --△pu = |u|^q-1 u + λf(x) on Ω, and identify necessary and sufficient conditions on Ω and f(x) which ensure the existence, or multiplicities of nonnegative solutions for the problem under consideration.
文摘In this paper, the fluid flow differential equation based on the homogenous reservoirs model is first reviewed. Then a theorem about the formal similarity of solutions in the Laplace space with outer boundary conditions and inner boundary condition is presented and proved. Lastly, a corollary of our theorem is given particularly on inner boundary. The obtained results are very helpful for understanding inherent laws of relevant engineering science and designing practical analysis software.
基金Supported by the National Science Foundation of China(11071245 and 11101418)
文摘In this paper, we deal with the following problem:By variational method, we prove the existenceof a nontrivial weak solution whenand the existence of a cylindricalweak solution when
文摘Structure of nonnegative nontrivial and positive solutions was precisely studied for some singularly perturbed p-Laplace equations. By virtue of sub- and supersolution method, it is shown that there are many nonnegative nontrivial spike-layer solutions and positive intermediate spike-layer solutions. Moreover, the upper and lower bound on the measure of each spike-layer were estimated when the parameter is sufficiently small.
文摘We discuss the solution of Laplace’s differential equation and a fractional differential equation of that type, by using analytic continuations of Riemann-Liouville fractional derivative and of Laplace transform. We show that the solutions, which are obtained by using operational calculus in the framework of distribution theory in our preceding papers, are obtained also by the present method.
文摘The high-order boundary conditions for the problems cf Laplace equation in infinite region have been developed. The improvement in accuracy for numerical solution is achieved by imposing the high-order boundary conditions on the exterior boundarv of a reduced finite region in which the numerical method is used. So both the computing efforts and the required storage in computer are reduced. The numerical examples show that the 1st-order boundary condition approaches to the exact boundary condition and it is clearly superior to the traditional boundary condition and the 2nd-order boundary condition.