In this paper, a corrected particle method based on the smoothed particle hydrodynamics (SPH) method with high-order Taylor expansion (CSPH-HT) for solving the viscoelastic flow is proposed and investigated. The valid...In this paper, a corrected particle method based on the smoothed particle hydrodynamics (SPH) method with high-order Taylor expansion (CSPH-HT) for solving the viscoelastic flow is proposed and investigated. The validity and merits of the CSPH-HT method are first tested by solving the nonlinear high order Kuramoto-Sivishinsky equation and simulating the drop stretching, respectively. Then the flow behaviors behind two stationary tangential cylinders of polymer melt, which have been received little attention, are investigated by the CSPH-HT method. Finally, the CSPH-HT method is extended to the simulation of the filling process of the viscoelastic fluid. The numerical results show that the CSPH-HT method possesses higher accuracy and stability than other corrected SPH methods and is more reliable than other corrected SPH methods.展开更多
The quality factor Q is an important parameter because it can refl ect the reservoir attenuated features and can be used for inverse-Q filtering to compensate for the seismic wave energy.The accuracy of the Q estimati...The quality factor Q is an important parameter because it can refl ect the reservoir attenuated features and can be used for inverse-Q filtering to compensate for the seismic wave energy.The accuracy of the Q estimation is greatly significant for improving the precision of the reservoir prediction and the resolution of seismic data.In this paper,the Q estimation formulas of the single-frequency point are derived on the basis of a diff erent-order Taylor series expansion of the amplitude attenuated factor.Moreover,the multifrequency point average(MFPA)method is introduced to obtain a stable Q estimation.The model tests demonstrate that the MFPA method is less aff ected by the frequency band,travel time diff erence,time window width,and noise interference than the logical spectrum ratio(LSR)method and the energy ratio(ER)method and has a higher Q estimation accuracy.In addition,the proposed method can be applied to post-stack seismic data and obtain eff ective Q values of complex models.When the MFPA method was applied to real marine seismic data,the Q values estimated by the MFPA method with the 1st–4th order showed good consistency with each other.In contrast,the Q values obtained by the ER method were larger than those of the proposed method,while those estimated by the LSR method signifi cantly deviated from the average values.In conclusion,the MFPA method has superior stability and practicability for the Q estimation.展开更多
Based on the Taylor series method and Li’s spatial differential method, a high-order hybrid Taylor–Li scheme is proposed.The results of a linear advection equation indicate that, using the initial values of the squa...Based on the Taylor series method and Li’s spatial differential method, a high-order hybrid Taylor–Li scheme is proposed.The results of a linear advection equation indicate that, using the initial values of the square-wave type, a result with thirdorder accuracy occurs. However, using initial values associated with the Gaussian function type, a result with very high precision appears. The study demonstrates that, when the order of the time integral is more than three, the corresponding optimal spatial difference order could be higher than six. The results indicate that the reason for why there is no improvement related to an order of spatial difference above six is the use of a time integral scheme that is not high enough. The author also proposes a recursive differential method to improve the Taylor–Li scheme’s computation speed. A more rapid and highprecision program than direct computation of the high-order space differential item is employed, and the computation speed is dramatically boosted. Based on a multiple-precision library, the ultrahigh-order Taylor–Li scheme can be used to solve the advection equation and Burgers’ equation.展开更多
A new numerical method of integrating the nonlinear evolution equations, namely the Taylor expansion method, was presented. The standard Galerkin method can be viewed as the 0_th order Taylor expansion method; while t...A new numerical method of integrating the nonlinear evolution equations, namely the Taylor expansion method, was presented. The standard Galerkin method can be viewed as the 0_th order Taylor expansion method; while the nonlinear Galerkin method can be viewed as the 1_st order modified Taylor expansion method. Moreover, the existence of the numerical solution and its convergence rate were proven. Finally, a concrete example, namely, the two_dimensional Navier_Stokes equations with a non slip boundary condition,was provided. The result is that the higher order Taylor expansion method is of the higher convergence rate under some assumptions about the regularity of the solution.展开更多
A variation of the direct Taylor expansion algorithm is suggested and applied to several linear and nonlinear differential equations of interest in physics and engineering, and the results are compared with those obta...A variation of the direct Taylor expansion algorithm is suggested and applied to several linear and nonlinear differential equations of interest in physics and engineering, and the results are compared with those obtained from other algorithms. It is shown that the suggested algorithm competes strongly with other existing algorithms, both in accuracy and ease of application, while demanding a shorter computation time.展开更多
This study presents an experiment of improving the performance of spectral stochastic finite element method using high-order elements. This experiment is implemented through a two-dimensional spectral stochastic finit...This study presents an experiment of improving the performance of spectral stochastic finite element method using high-order elements. This experiment is implemented through a two-dimensional spectral stochastic finite element formulation of an elliptic partial differential equation having stochastic coefficients. Deriving this spectral stochastic finite element formulation couples a two-dimensional deterministic finite element formulation of an elliptic partial differential equation with generalized polynomial chaos expansions of stochastic coefficients. Further inspection of the performance of resulting spectral stochastic finite element formulation with adopting linear and quadratic (9-node or 8-node) quadrilateral elements finds that more accurate standard deviations of unknowns are surprisingly predicted using quadratic quadrilateral elements, especially under high autocorrelation function values of stochastic coefficients. In addition, creating spectral stochastic finite element results using quadratic quadrilateral elements is not unacceptably time-consuming. Therefore, this study concludes that adopting high-order elements can be a lower-cost method to improve the performance of spectral stochastic finite element method.展开更多
基金support of the National Natural Science Foundation of China (Grants 11501495, 51541912, 51409227)the Natural Science Foundation of Jiangsu Province, China (Grants BK20130436, BK20150436)+1 种基金the Postdoctoral Science Foundation of China (Grants 2014M550310, 2015M581869, 2015T80589)the Natural Science Foundation of the Higher Education Institutions of Jiangsu Province (Grant 15KJB110025)
文摘In this paper, a corrected particle method based on the smoothed particle hydrodynamics (SPH) method with high-order Taylor expansion (CSPH-HT) for solving the viscoelastic flow is proposed and investigated. The validity and merits of the CSPH-HT method are first tested by solving the nonlinear high order Kuramoto-Sivishinsky equation and simulating the drop stretching, respectively. Then the flow behaviors behind two stationary tangential cylinders of polymer melt, which have been received little attention, are investigated by the CSPH-HT method. Finally, the CSPH-HT method is extended to the simulation of the filling process of the viscoelastic fluid. The numerical results show that the CSPH-HT method possesses higher accuracy and stability than other corrected SPH methods and is more reliable than other corrected SPH methods.
基金supported by The National Natural Science Foundation (Grant Nos.41874126, 42004114)the Key Research and development project of Jiangxi Province in China (Grant No.20192ACB80006)+1 种基金the Natural Science Foundation of Jiangxi Province (Grant Nos. 20202BAB211010, 20212BAB203005)Open Foundation of State Key Laboratory of Nuclear Resources and Environment (2020NRE25)
文摘The quality factor Q is an important parameter because it can refl ect the reservoir attenuated features and can be used for inverse-Q filtering to compensate for the seismic wave energy.The accuracy of the Q estimation is greatly significant for improving the precision of the reservoir prediction and the resolution of seismic data.In this paper,the Q estimation formulas of the single-frequency point are derived on the basis of a diff erent-order Taylor series expansion of the amplitude attenuated factor.Moreover,the multifrequency point average(MFPA)method is introduced to obtain a stable Q estimation.The model tests demonstrate that the MFPA method is less aff ected by the frequency band,travel time diff erence,time window width,and noise interference than the logical spectrum ratio(LSR)method and the energy ratio(ER)method and has a higher Q estimation accuracy.In addition,the proposed method can be applied to post-stack seismic data and obtain eff ective Q values of complex models.When the MFPA method was applied to real marine seismic data,the Q values estimated by the MFPA method with the 1st–4th order showed good consistency with each other.In contrast,the Q values obtained by the ER method were larger than those of the proposed method,while those estimated by the LSR method signifi cantly deviated from the average values.In conclusion,the MFPA method has superior stability and practicability for the Q estimation.
基金supported by the National Natural Sciences Foundation of China(Grant Nos.41375112 and 41530426)the Chinese Academy of Sciences Key Technology Talent Program
文摘Based on the Taylor series method and Li’s spatial differential method, a high-order hybrid Taylor–Li scheme is proposed.The results of a linear advection equation indicate that, using the initial values of the square-wave type, a result with thirdorder accuracy occurs. However, using initial values associated with the Gaussian function type, a result with very high precision appears. The study demonstrates that, when the order of the time integral is more than three, the corresponding optimal spatial difference order could be higher than six. The results indicate that the reason for why there is no improvement related to an order of spatial difference above six is the use of a time integral scheme that is not high enough. The author also proposes a recursive differential method to improve the Taylor–Li scheme’s computation speed. A more rapid and highprecision program than direct computation of the high-order space differential item is employed, and the computation speed is dramatically boosted. Based on a multiple-precision library, the ultrahigh-order Taylor–Li scheme can be used to solve the advection equation and Burgers’ equation.
文摘A new numerical method of integrating the nonlinear evolution equations, namely the Taylor expansion method, was presented. The standard Galerkin method can be viewed as the 0_th order Taylor expansion method; while the nonlinear Galerkin method can be viewed as the 1_st order modified Taylor expansion method. Moreover, the existence of the numerical solution and its convergence rate were proven. Finally, a concrete example, namely, the two_dimensional Navier_Stokes equations with a non slip boundary condition,was provided. The result is that the higher order Taylor expansion method is of the higher convergence rate under some assumptions about the regularity of the solution.
文摘A variation of the direct Taylor expansion algorithm is suggested and applied to several linear and nonlinear differential equations of interest in physics and engineering, and the results are compared with those obtained from other algorithms. It is shown that the suggested algorithm competes strongly with other existing algorithms, both in accuracy and ease of application, while demanding a shorter computation time.
文摘This study presents an experiment of improving the performance of spectral stochastic finite element method using high-order elements. This experiment is implemented through a two-dimensional spectral stochastic finite element formulation of an elliptic partial differential equation having stochastic coefficients. Deriving this spectral stochastic finite element formulation couples a two-dimensional deterministic finite element formulation of an elliptic partial differential equation with generalized polynomial chaos expansions of stochastic coefficients. Further inspection of the performance of resulting spectral stochastic finite element formulation with adopting linear and quadratic (9-node or 8-node) quadrilateral elements finds that more accurate standard deviations of unknowns are surprisingly predicted using quadratic quadrilateral elements, especially under high autocorrelation function values of stochastic coefficients. In addition, creating spectral stochastic finite element results using quadratic quadrilateral elements is not unacceptably time-consuming. Therefore, this study concludes that adopting high-order elements can be a lower-cost method to improve the performance of spectral stochastic finite element method.
