Based on the vibrational potential curves coupled with the minimum energy reaction path, the partial potential energy surface of the reaction I+HI→IH+I was constructed at the QCISD(T)//MP4SDQ level with pseudo po...Based on the vibrational potential curves coupled with the minimum energy reaction path, the partial potential energy surface of the reaction I+HI→IH+I was constructed at the QCISD(T)//MP4SDQ level with pseudo potential method. And the formation mechanism of the scattering resonance states of this reaction was well interpreted with the partial potential energy surface. The scattering resonance states of this reaction should belong to Feshbach resonance because of the coupling of the vibrational mode and the translational mode. With the one-dimensional square potential well model, the resonance width and lifetime of the I+HI(v=0)→IH(v'=0)+I state-to-state reaction were calculated, which preferably explained the high-resolved threshold photodetachment spectroscopy of the IHI- anion performed by Neumark et al..展开更多
An efficient conformal locally one-dimensional finite-difference time-domain(LOD-CFDTD) method is presented for solving two-dimensional(2D) electromagnetic(EM) scattering problems. The formulation for the 2D tra...An efficient conformal locally one-dimensional finite-difference time-domain(LOD-CFDTD) method is presented for solving two-dimensional(2D) electromagnetic(EM) scattering problems. The formulation for the 2D transverse-electric(TE) case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit(ADI) CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field(TF/SF) boundary and the perfectly matched layer(PML), the radar cross section(RCS) of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.展开更多
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.展开更多
In this paper, the coupling schemes of atmosphere-ocean climate models are discussed with one-dimensional advection equations. The convergence and stability for synchronous and asynchronous schemes are demonstrated an...In this paper, the coupling schemes of atmosphere-ocean climate models are discussed with one-dimensional advection equations. The convergence and stability for synchronous and asynchronous schemes are demonstrated and compared.Conclusions inferred from the analysis are given below. The synchronous scheme as well as the asynchronous-implicit scheme in this model are stable for arbitrary integrating time intervals. The asynchronous explicit scheme is unstable under certain conditions, which depend upon advection velocities and heat exchange parameters in the atmosphere and oceans. With both synchronous and asynchronous stable schemes the discrete solutions converge to their unique exact ones. Advections in the atmosphere and ocean accelerate the rate of convergence of the asynchronous-implicit scheme. It is suggusted that the asynchronous-implicit coupling scheme is a stable and efficient method for most climatic simulations.展开更多
In order to use the second-order 5-point difference scheme mentioned to compute the solution of one dimension unsteady equations of the direct reflection of the strong plane detonation wave meeting a solid wall barrie...In order to use the second-order 5-point difference scheme mentioned to compute the solution of one dimension unsteady equations of the direct reflection of the strong plane detonation wave meeting a solid wall barrier,in this paper,we technically construct the difference schemes of the boundary and sub-boundary of the problem,and deduce the auto-analogue analytic solutions of the initial value problem,and at the same time,we present a method for the singular property of the initial value problem,from which we can get a satisfactory computation result of this difficult problem.The difference scheme used in this paper to deal with the discontinuity problems of the shock wave are valuable and worth generalization.展开更多
In this work we generate the numerical solutions of the Burgers’ equation by applying the Crank-Nicolson method directly to the Burgers’ equation, i.e., we do not use Hopf-Cole transformation to reduce Burgers’ equ...In this work we generate the numerical solutions of the Burgers’ equation by applying the Crank-Nicolson method directly to the Burgers’ equation, i.e., we do not use Hopf-Cole transformation to reduce Burgers’ equation into the linear heat equation. Absolute error of the present method is compared to the absolute error of the two existing methods for two test problems. The method is also analyzed for a third test problem, nu-merical solutions as well as exact solutions for different values of viscosity are calculated and we find that the numerical solutions are very close to exact solution.展开更多
Stochastic partial differential equations (SPDEs) describe the dynamics of stochastic processes depending on space-time continuum. These equations have been widely used to model many applications in engineering and ma...Stochastic partial differential equations (SPDEs) describe the dynamics of stochastic processes depending on space-time continuum. These equations have been widely used to model many applications in engineering and mathematical sciences. In this paper we use three finite difference schemes in order to approximate the solution of stochastic parabolic partial differential equations. The conditions of the mean square convergence of the numerical solution are studied. Some case studies are discussed.展开更多
The multi-dimensional system of nonlinear partial differential equations is considered. In two-dimensional case, this system describes process of vein formation in higher plants. Variable directions finite difference ...The multi-dimensional system of nonlinear partial differential equations is considered. In two-dimensional case, this system describes process of vein formation in higher plants. Variable directions finite difference scheme is constructed. The stability and convergence of that scheme are studied. Numerical experiments are carried out. The appropriate graphical illustrations and tables are given.展开更多
Alternating direction implicit (A.D.I.) schemes have been proved valuable in the approximation of the solutions of parabolic partial differential equations in multi-dimensional space. Consider equations in the form pa...Alternating direction implicit (A.D.I.) schemes have been proved valuable in the approximation of the solutions of parabolic partial differential equations in multi-dimensional space. Consider equations in the form partial derivative u/partial derivative t - partial derivative/partial derivative x(a(x,y,t) partial derivative u/partial derivative x) - partial derivative/partial derivative y(b(x,y,t) partial derivative u partial derivative y) = f Two A.D.I. schemes, Peaceman-Rachford scheme and Douglas scheme will be studied. In the literature, stability and convergence have been analysed with Fourier Method, which cannot be extended beyond the model problem with constant coefficients. Additionally, L-2 energy method has been introduced to analyse the case of non-constant coefficients, however, the conclusions are too weak and incomplete because of the so-called 'equivalence between L-2 norm and H-1 semi-norm'. In this paper, we try to improve these conclusions by H-1 energy estimating method. The principal results are that both of the two A.D.I. schemes are absolutely stable and converge to the exact solution with error estimations O(Delta t(2) + h(2)) in discrete H-1 norm. This implies essential improvement of existing conclusions.展开更多
In this paper, an explicit three_level symmetrical differencing scheme with parameters for solving parabolic partial differential equation of three_dimension will be considered. The stability condition and local trunc...In this paper, an explicit three_level symmetrical differencing scheme with parameters for solving parabolic partial differential equation of three_dimension will be considered. The stability condition and local truncation error for the scheme are r<1/2 and O( Δ t 2+ Δ x 4+ Δ y 4+ Δ z 4) ,respectively.展开更多
This paper is devoted to the study of a three-dimensional delayed system with nonlocal diffusion and partial quasi-monotonicity. By developing a new definition of upper-lower solutions and a new cross iteration scheme...This paper is devoted to the study of a three-dimensional delayed system with nonlocal diffusion and partial quasi-monotonicity. By developing a new definition of upper-lower solutions and a new cross iteration scheme, we establish some existence results of traveling wave solutions. The results are applied to a nonlocal diffusion model which takes the three-species Lotka-Volterra model as its special case.展开更多
In this paper, we consider a singular perturbation elliptic-parabolic partial differential equation for periodic boundary value problem, and construct a difference scheme. Using the method of decomposing the singular ...In this paper, we consider a singular perturbation elliptic-parabolic partial differential equation for periodic boundary value problem, and construct a difference scheme. Using the method of decomposing the singular term from its solution and combining an asymptotic expansion of the equation, we prove that the scheme constructed by this paper converges uniformly to the solution of its original problem with O(r+h2).展开更多
This paper improves and generalizes the two difference schemes presented in paper [1] and gives a new difference scheme for second order linear elliptic partial differential equations, its difference matrix is a matri...This paper improves and generalizes the two difference schemes presented in paper [1] and gives a new difference scheme for second order linear elliptic partial differential equations, its difference matrix is a matrix and because of the stability of the M-matrix, it is convergent by the asynchronous iterative method on multiprocessors. Then this paper gives a class of differeifce schemes for linear elliptic PDEs so that their difference matrixes are all M-matrixes and their asynchronous parallel computation are convergent.展开更多
In this paper, we consider the upwind difference scheme for singular perturbation problem (1.1). On a special discretization mesh, it is proved that the solution of the upwind difference scheme is first order converge...In this paper, we consider the upwind difference scheme for singular perturbation problem (1.1). On a special discretization mesh, it is proved that the solution of the upwind difference scheme is first order convergent, uniformly in the small parameter e , to the solution of problem (1.1). Numerical results are finally provided.展开更多
The goal of computational science is to develop models that predict phenomena observed in nature. However, these models are often based on parameters that are uncertain. In recent decades, main numerical methods for s...The goal of computational science is to develop models that predict phenomena observed in nature. However, these models are often based on parameters that are uncertain. In recent decades, main numerical methods for solving SPDEs have been used such as, finite difference and finite element schemes [1]-[5]. Also, some practical techniques like the method of lines for boundary value problems have been applied to the linear stochastic partial differential equations, and the outcomes of these approaches have been experimented numerically [7]. In [8]-[10], the author discussed mean square convergent finite difference method for solving some random partial differential equations. Random numerical techniques for both ordinary and partial random differential equations are treated in [4] [10]. As regards applications using explicit analytic solutions or numerical methods, a few results may be found in [5] [6] [11]. This article focuses on solving random heat equation by using Crank-Nicol- son technique under mean square sense and it is organized as follows. In Section 2, the mean square calculus preliminaries that will be required throughout the paper are presented. In Section 3, the Crank-Nicolson scheme for solving the random heat equation is presented. In Section 4, some case studies are showed. Short conclusions are cleared in the end section.展开更多
The traditional combinatorial designs can be used as basic designs for constructing designs of computer experiments which have been used successfully till now in various domains such as engineering, pharmaceutical ind...The traditional combinatorial designs can be used as basic designs for constructing designs of computer experiments which have been used successfully till now in various domains such as engineering, pharmaceutical industry, etc. In this paper, a new series of generalized partially balanced incomplete blocks PBIB designs with m associated classes (m = 4, 5 and 7) based on new generalized association schemes with number of treatments v arranged in w arrays of n rows and l columns (w ≥ 2, n ≥ 2, l ≥ 2) is defined. Some construction methods of these new PBIB are given and their parameters are specified using the Combinatory Method (s). For n or l even and s divisor of n or l, the obtained PBIB designs are resolvable PBIB designs. So the Fang RBIBD method is applied to obtain a series of particular U-type designs U (wnl;) (r is the repetition number of each treatment in our resolvable PBIB design).展开更多
High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of th...High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of the WENO algorithm,large amount of computational costs are required for solving multidimensional problems.In our previous work(Lu et al.in Pure Appl Math Q 14:57–86,2018;Zhu and Zhang in J Sci Comput 87:44,2021),sparse-grid techniques were applied to the classical finite difference WENO schemes in solving multidimensional hyperbolic equations,and it was shown that significant CPU times were saved,while both accuracy and stability of the classical WENO schemes were maintained for computations on sparse grids.In this technical note,we apply the approach to recently developed finite difference multi-resolution WENO scheme specifically the fifth-order scheme,which has very interesting properties such as its simplicity in linear weights’construction over a classical WENO scheme.Numerical experiments on solving high dimensional hyperbolic equations including Vlasov based kinetic problems are performed to demonstrate that the sparse-grid computations achieve large savings of CPU times,and at the same time preserve comparable accuracy and resolution with those on corresponding regular single grids.展开更多
基金Ⅴ. ACKN0WLEDGMENTS This work was supported by the National Natural Science Foundation of China (No.20573064) and Ph.D. Special Research Foundation of Chinese Education Department.
文摘Based on the vibrational potential curves coupled with the minimum energy reaction path, the partial potential energy surface of the reaction I+HI→IH+I was constructed at the QCISD(T)//MP4SDQ level with pseudo potential method. And the formation mechanism of the scattering resonance states of this reaction was well interpreted with the partial potential energy surface. The scattering resonance states of this reaction should belong to Feshbach resonance because of the coupling of the vibrational mode and the translational mode. With the one-dimensional square potential well model, the resonance width and lifetime of the I+HI(v=0)→IH(v'=0)+I state-to-state reaction were calculated, which preferably explained the high-resolved threshold photodetachment spectroscopy of the IHI- anion performed by Neumark et al..
基金supported by the National Natural Science Foundation of China(Grant Nos.61331007 and 61471105)
文摘An efficient conformal locally one-dimensional finite-difference time-domain(LOD-CFDTD) method is presented for solving two-dimensional(2D) electromagnetic(EM) scattering problems. The formulation for the 2D transverse-electric(TE) case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit(ADI) CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field(TF/SF) boundary and the perfectly matched layer(PML), the radar cross section(RCS) of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.
文摘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.
文摘In this paper, the coupling schemes of atmosphere-ocean climate models are discussed with one-dimensional advection equations. The convergence and stability for synchronous and asynchronous schemes are demonstrated and compared.Conclusions inferred from the analysis are given below. The synchronous scheme as well as the asynchronous-implicit scheme in this model are stable for arbitrary integrating time intervals. The asynchronous explicit scheme is unstable under certain conditions, which depend upon advection velocities and heat exchange parameters in the atmosphere and oceans. With both synchronous and asynchronous stable schemes the discrete solutions converge to their unique exact ones. Advections in the atmosphere and ocean accelerate the rate of convergence of the asynchronous-implicit scheme. It is suggusted that the asynchronous-implicit coupling scheme is a stable and efficient method for most climatic simulations.
文摘In order to use the second-order 5-point difference scheme mentioned to compute the solution of one dimension unsteady equations of the direct reflection of the strong plane detonation wave meeting a solid wall barrier,in this paper,we technically construct the difference schemes of the boundary and sub-boundary of the problem,and deduce the auto-analogue analytic solutions of the initial value problem,and at the same time,we present a method for the singular property of the initial value problem,from which we can get a satisfactory computation result of this difficult problem.The difference scheme used in this paper to deal with the discontinuity problems of the shock wave are valuable and worth generalization.
文摘In this work we generate the numerical solutions of the Burgers’ equation by applying the Crank-Nicolson method directly to the Burgers’ equation, i.e., we do not use Hopf-Cole transformation to reduce Burgers’ equation into the linear heat equation. Absolute error of the present method is compared to the absolute error of the two existing methods for two test problems. The method is also analyzed for a third test problem, nu-merical solutions as well as exact solutions for different values of viscosity are calculated and we find that the numerical solutions are very close to exact solution.
文摘Stochastic partial differential equations (SPDEs) describe the dynamics of stochastic processes depending on space-time continuum. These equations have been widely used to model many applications in engineering and mathematical sciences. In this paper we use three finite difference schemes in order to approximate the solution of stochastic parabolic partial differential equations. The conditions of the mean square convergence of the numerical solution are studied. Some case studies are discussed.
文摘The multi-dimensional system of nonlinear partial differential equations is considered. In two-dimensional case, this system describes process of vein formation in higher plants. Variable directions finite difference scheme is constructed. The stability and convergence of that scheme are studied. Numerical experiments are carried out. The appropriate graphical illustrations and tables are given.
文摘Alternating direction implicit (A.D.I.) schemes have been proved valuable in the approximation of the solutions of parabolic partial differential equations in multi-dimensional space. Consider equations in the form partial derivative u/partial derivative t - partial derivative/partial derivative x(a(x,y,t) partial derivative u/partial derivative x) - partial derivative/partial derivative y(b(x,y,t) partial derivative u partial derivative y) = f Two A.D.I. schemes, Peaceman-Rachford scheme and Douglas scheme will be studied. In the literature, stability and convergence have been analysed with Fourier Method, which cannot be extended beyond the model problem with constant coefficients. Additionally, L-2 energy method has been introduced to analyse the case of non-constant coefficients, however, the conclusions are too weak and incomplete because of the so-called 'equivalence between L-2 norm and H-1 semi-norm'. In this paper, we try to improve these conclusions by H-1 energy estimating method. The principal results are that both of the two A.D.I. schemes are absolutely stable and converge to the exact solution with error estimations O(Delta t(2) + h(2)) in discrete H-1 norm. This implies essential improvement of existing conclusions.
文摘In this paper, an explicit three_level symmetrical differencing scheme with parameters for solving parabolic partial differential equation of three_dimension will be considered. The stability condition and local truncation error for the scheme are r<1/2 and O( Δ t 2+ Δ x 4+ Δ y 4+ Δ z 4) ,respectively.
基金Supported by the Natural Science Foundation of China (11171120)the Doctoral Program of Higher Education of China (20094407110001)Natural Science Foundation of Guangdong Province (10151063101000003)
文摘This paper is devoted to the study of a three-dimensional delayed system with nonlocal diffusion and partial quasi-monotonicity. By developing a new definition of upper-lower solutions and a new cross iteration scheme, we establish some existence results of traveling wave solutions. The results are applied to a nonlocal diffusion model which takes the three-species Lotka-Volterra model as its special case.
基金This work is supported by the National Fujian Province Nature Science Research Funds
文摘In this paper, we consider a singular perturbation elliptic-parabolic partial differential equation for periodic boundary value problem, and construct a difference scheme. Using the method of decomposing the singular term from its solution and combining an asymptotic expansion of the equation, we prove that the scheme constructed by this paper converges uniformly to the solution of its original problem with O(r+h2).
文摘This paper improves and generalizes the two difference schemes presented in paper [1] and gives a new difference scheme for second order linear elliptic partial differential equations, its difference matrix is a matrix and because of the stability of the M-matrix, it is convergent by the asynchronous iterative method on multiprocessors. Then this paper gives a class of differeifce schemes for linear elliptic PDEs so that their difference matrixes are all M-matrixes and their asynchronous parallel computation are convergent.
文摘In this paper, we consider the upwind difference scheme for singular perturbation problem (1.1). On a special discretization mesh, it is proved that the solution of the upwind difference scheme is first order convergent, uniformly in the small parameter e , to the solution of problem (1.1). Numerical results are finally provided.
文摘The goal of computational science is to develop models that predict phenomena observed in nature. However, these models are often based on parameters that are uncertain. In recent decades, main numerical methods for solving SPDEs have been used such as, finite difference and finite element schemes [1]-[5]. Also, some practical techniques like the method of lines for boundary value problems have been applied to the linear stochastic partial differential equations, and the outcomes of these approaches have been experimented numerically [7]. In [8]-[10], the author discussed mean square convergent finite difference method for solving some random partial differential equations. Random numerical techniques for both ordinary and partial random differential equations are treated in [4] [10]. As regards applications using explicit analytic solutions or numerical methods, a few results may be found in [5] [6] [11]. This article focuses on solving random heat equation by using Crank-Nicol- son technique under mean square sense and it is organized as follows. In Section 2, the mean square calculus preliminaries that will be required throughout the paper are presented. In Section 3, the Crank-Nicolson scheme for solving the random heat equation is presented. In Section 4, some case studies are showed. Short conclusions are cleared in the end section.
文摘The traditional combinatorial designs can be used as basic designs for constructing designs of computer experiments which have been used successfully till now in various domains such as engineering, pharmaceutical industry, etc. In this paper, a new series of generalized partially balanced incomplete blocks PBIB designs with m associated classes (m = 4, 5 and 7) based on new generalized association schemes with number of treatments v arranged in w arrays of n rows and l columns (w ≥ 2, n ≥ 2, l ≥ 2) is defined. Some construction methods of these new PBIB are given and their parameters are specified using the Combinatory Method (s). For n or l even and s divisor of n or l, the obtained PBIB designs are resolvable PBIB designs. So the Fang RBIBD method is applied to obtain a series of particular U-type designs U (wnl;) (r is the repetition number of each treatment in our resolvable PBIB design).
文摘High-order accurate weighted essentially non-oscillatory(WENO)schemes are a class of broadly applied numerical methods for solving hyperbolic partial differential equations(PDEs).Due to highly nonlinear property of the WENO algorithm,large amount of computational costs are required for solving multidimensional problems.In our previous work(Lu et al.in Pure Appl Math Q 14:57–86,2018;Zhu and Zhang in J Sci Comput 87:44,2021),sparse-grid techniques were applied to the classical finite difference WENO schemes in solving multidimensional hyperbolic equations,and it was shown that significant CPU times were saved,while both accuracy and stability of the classical WENO schemes were maintained for computations on sparse grids.In this technical note,we apply the approach to recently developed finite difference multi-resolution WENO scheme specifically the fifth-order scheme,which has very interesting properties such as its simplicity in linear weights’construction over a classical WENO scheme.Numerical experiments on solving high dimensional hyperbolic equations including Vlasov based kinetic problems are performed to demonstrate that the sparse-grid computations achieve large savings of CPU times,and at the same time preserve comparable accuracy and resolution with those on corresponding regular single grids.