This paper develops a numerical method to invert multi-dimensional Laplace transforms. By a variable transform, Laplace transforms are converted to multi-dimensional Hansdorff moment problems so that the numerical sol...This paper develops a numerical method to invert multi-dimensional Laplace transforms. By a variable transform, Laplace transforms are converted to multi-dimensional Hansdorff moment problems so that the numerical solution can be achieved. Stability estimation is also obtained. Numerical simulations show the efficiency and practicality of the method.展开更多
This paper examines the performance of five algorithms for numerically inverting the Laplace transform, in standard, 16-digit and multi-precision environments. The algorithms are taken from three of the four main clas...This paper examines the performance of five algorithms for numerically inverting the Laplace transform, in standard, 16-digit and multi-precision environments. The algorithms are taken from three of the four main classes of numerical methods used to invert the Laplace transform. Because the numerical inversion of the Laplace transform is a perturbed problem, rounding errors which are generated in numerical approximations can adversely affect the accurate reconstruction of the inverse transform. This paper demonstrates that working in a multi-precision environment can substantially reduce these errors and the resulting perturbations exist in transforming the data from the s-space into the time domain and in so doing overcome the main drawback of numerically inverting the Laplace transform. Our main finding is that both the Talbot and the accelerated Gaver functionals perform considerably better in a multi-precision environment increasing the advantages of using Laplace transform methods over time-stepping procedures in solving diffusion and more generally parabolic partial differential equations.展开更多
A new provement of the existence and uniqueness about periodic boundary value Duffing equation is established by using global inverse function theorem. An algorithm for solving differential equation that has a large c...A new provement of the existence and uniqueness about periodic boundary value Duffing equation is established by using global inverse function theorem. An algorithm for solving differential equation that has a large convergence domain is given. Finally, a numerical example is given.展开更多
In order to find stable, accurate, and computationally efficient methods for performing the inverse Laplace transform, a new double transformation approach is proposed. To validate and improve the inversion solution o...In order to find stable, accurate, and computationally efficient methods for performing the inverse Laplace transform, a new double transformation approach is proposed. To validate and improve the inversion solution obtained using the Gaver-Stehfest algorithm, direct Laplace transforms are taken of the numerically inverted transforms to compare with the original function. The numerical direct Laplace transform is implemented with a composite Simpson’s rule. Challenging numerical examples involving periodic and oscillatory functions, are investigated. The numerical examples illustrate the computational accuracy and efficiency of the direct Laplace transform and its inverse due to increasing the precision level and the number of terms included in the expansion. It is found that the number of expansion terms and the precision level selected must be in a harmonious balance in order for correct and stable results to be obtained.展开更多
Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harm...Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harmonic excitation.Within any integration segment, as long as its length is small enough,along the circumferential curvilinear coordinate,the non- homogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper.For the exterior fluid domain,the multi-circular virtual source simulation technique is adopted.The source density distributed on each virtual circular curve may be ex- panded as the Fourier's series.Combining with the inverse fast Fourier transformation,a higher accuracy and efficiency method for solving 2-D exterior Helmholtz's problems is presented in this paper.In the aspect of solution to the coupling equations,the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper.According to the superposition princi- ple and compatibility conditions at the interface between the elastic ring and fluid,the algebraic equation of system can be directly constructed by using the least square approximation.Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented.Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.展开更多
In a previous paper, the author proposed a parameter inversionscheme for layered-inhomogeneous media. Here, by constructing the operatorRx for the cases of single layer(n=1 ) and double layers(n= 2), numerical im-plem...In a previous paper, the author proposed a parameter inversionscheme for layered-inhomogeneous media. Here, by constructing the operatorRx for the cases of single layer(n=1 ) and double layers(n= 2), numerical im-plementation is presented. Results in this展开更多
In this article, we have effectively used the Numerical Inversion of Laplace transform to study the time-dependent thin film flow of a second grade fluid flowing down an inclined plane through a porous medium. The sol...In this article, we have effectively used the Numerical Inversion of Laplace transform to study the time-dependent thin film flow of a second grade fluid flowing down an inclined plane through a porous medium. The solution to the governing equation is obtained by using the standard Laplace transform. However, to transform the obtained solutions from Laplace space back to the original space, we have used the Numerical Inversion of Laplace transform. Graphical results have been presented to show the effects of different parameters involved and to show how the fluid flow evolves with time.展开更多
In this paper, a hybrid method is introduced briefly to predict the behavior of the non-linear partial differential equations. The method is hybrid in the sense that different numerical methods, differential transform...In this paper, a hybrid method is introduced briefly to predict the behavior of the non-linear partial differential equations. The method is hybrid in the sense that different numerical methods, differential transform and finite differences, are used in different subdomains. Our aim of this approach is to combine the flexibility of differential transform and the efficiency of finite differences. An explicit hybrid method for the transient response of inhomogeneous nonlinear partial differential equations is presented;applying finite difference scheme on the fixed grid size is used to approximate the space discretisation, whereas the differential transform method is used for time operator. Comparison of the efficiency of the different approaches is a very important aspect of this study. In our test cases, the hybrid approach is faster than the corresponding highly optimized finite difference method in two dimensional computations. We compared our hybrid approach’s results with the exact and/or numerical solutions of PDE which obtained from Adomian Decomposition Method. Results show that the hybrid approach may be an important tool to reduce the execution time and memory requirements for large scale computations and get remarkable results in predicting the solutions of nonlinear initial value problems.展开更多
基金the Jiangxi Provincial Natural Scientific Foundation(0211014)Scientific Research Program from Education Office of Jiangxi Province([2005]213)East China Institute of Technology.
文摘This paper develops a numerical method to invert multi-dimensional Laplace transforms. By a variable transform, Laplace transforms are converted to multi-dimensional Hansdorff moment problems so that the numerical solution can be achieved. Stability estimation is also obtained. Numerical simulations show the efficiency and practicality of the method.
文摘This paper examines the performance of five algorithms for numerically inverting the Laplace transform, in standard, 16-digit and multi-precision environments. The algorithms are taken from three of the four main classes of numerical methods used to invert the Laplace transform. Because the numerical inversion of the Laplace transform is a perturbed problem, rounding errors which are generated in numerical approximations can adversely affect the accurate reconstruction of the inverse transform. This paper demonstrates that working in a multi-precision environment can substantially reduce these errors and the resulting perturbations exist in transforming the data from the s-space into the time domain and in so doing overcome the main drawback of numerically inverting the Laplace transform. Our main finding is that both the Talbot and the accelerated Gaver functionals perform considerably better in a multi-precision environment increasing the advantages of using Laplace transform methods over time-stepping procedures in solving diffusion and more generally parabolic partial differential equations.
文摘A new provement of the existence and uniqueness about periodic boundary value Duffing equation is established by using global inverse function theorem. An algorithm for solving differential equation that has a large convergence domain is given. Finally, a numerical example is given.
文摘In order to find stable, accurate, and computationally efficient methods for performing the inverse Laplace transform, a new double transformation approach is proposed. To validate and improve the inversion solution obtained using the Gaver-Stehfest algorithm, direct Laplace transforms are taken of the numerically inverted transforms to compare with the original function. The numerical direct Laplace transform is implemented with a composite Simpson’s rule. Challenging numerical examples involving periodic and oscillatory functions, are investigated. The numerical examples illustrate the computational accuracy and efficiency of the direct Laplace transform and its inverse due to increasing the precision level and the number of terms included in the expansion. It is found that the number of expansion terms and the precision level selected must be in a harmonious balance in order for correct and stable results to be obtained.
基金Project supported by the National Natural Science Foundation of China (No.10172038)
文摘Based on the transfer matrix method and the virtual source simulation technique, this paper proposes a novel semi-analytical and semi-numerical method for solving 2-D sound- structure interaction problems under a harmonic excitation.Within any integration segment, as long as its length is small enough,along the circumferential curvilinear coordinate,the non- homogeneous matrix differential equation of an elastic ring of complex geometrical shape can be rewritten in terms of the homogeneous one by the method of extended homogeneous capacity proposed in this paper.For the exterior fluid domain,the multi-circular virtual source simulation technique is adopted.The source density distributed on each virtual circular curve may be ex- panded as the Fourier's series.Combining with the inverse fast Fourier transformation,a higher accuracy and efficiency method for solving 2-D exterior Helmholtz's problems is presented in this paper.In the aspect of solution to the coupling equations,the state vectors of elastic ring induced by the given harmonic excitation and generalized forces of coefficients of the Fourier series can be obtained respectively by using a high precision integration scheme combined with the method of extended homogeneous capacity put forward in this paper.According to the superposition princi- ple and compatibility conditions at the interface between the elastic ring and fluid,the algebraic equation of system can be directly constructed by using the least square approximation.Examples of acoustic radiation from two typical fluid-loaded elastic rings under a harmonic concentrated force are presented.Numerical results show that the method proposed is more efficient than the mixed FE-BE method in common use.
文摘In a previous paper, the author proposed a parameter inversionscheme for layered-inhomogeneous media. Here, by constructing the operatorRx for the cases of single layer(n=1 ) and double layers(n= 2), numerical im-plementation is presented. Results in this
文摘In this article, we have effectively used the Numerical Inversion of Laplace transform to study the time-dependent thin film flow of a second grade fluid flowing down an inclined plane through a porous medium. The solution to the governing equation is obtained by using the standard Laplace transform. However, to transform the obtained solutions from Laplace space back to the original space, we have used the Numerical Inversion of Laplace transform. Graphical results have been presented to show the effects of different parameters involved and to show how the fluid flow evolves with time.
文摘In this paper, a hybrid method is introduced briefly to predict the behavior of the non-linear partial differential equations. The method is hybrid in the sense that different numerical methods, differential transform and finite differences, are used in different subdomains. Our aim of this approach is to combine the flexibility of differential transform and the efficiency of finite differences. An explicit hybrid method for the transient response of inhomogeneous nonlinear partial differential equations is presented;applying finite difference scheme on the fixed grid size is used to approximate the space discretisation, whereas the differential transform method is used for time operator. Comparison of the efficiency of the different approaches is a very important aspect of this study. In our test cases, the hybrid approach is faster than the corresponding highly optimized finite difference method in two dimensional computations. We compared our hybrid approach’s results with the exact and/or numerical solutions of PDE which obtained from Adomian Decomposition Method. Results show that the hybrid approach may be an important tool to reduce the execution time and memory requirements for large scale computations and get remarkable results in predicting the solutions of nonlinear initial value problems.