A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to t...A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suit- able for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.展开更多
General purpose graphic processing unit(GPU) calculation technology is gradually widely used in various fields.Its mode of single instruction,multiple threads is capable of seismic numerical simulation which has a hug...General purpose graphic processing unit(GPU) calculation technology is gradually widely used in various fields.Its mode of single instruction,multiple threads is capable of seismic numerical simulation which has a huge quantity of data and calculation steps.In this study,we introduce a GPU-based parallel calculation method of a precise integration method(PIM) for seismic forward modeling.Compared with CPU single-core calculation,GPU parallel calculating perfectly keeps the features of PIM,which has small bandwidth,high accuracy and capability of modeling complex substructures,and GPU calculation brings high computational efficiency,which means that high-performing GPU parallel calculation can make seismic forward modeling closer to real seismic records.展开更多
A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in t...A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in the phase change problem.First,according to the continuity conditions of temperature and its gradient on the liquid-mushy interface,the mushy zone and the liquid phase in the simulation can be considered as a whole part,namely,the non-solid phase,and the change of latent heat is approximated by heat source which is dependent on temperature.Then,the precise integration BEM is used to obtain the differential equations in the solid phase zone and the non-solid phase zone,respectively.Moreover,an iterative predictor-corrector precise integration method(PIM)is needed to solve the differential equations and obtain the temperature field and the heat flux on the boundary.According to an energy balance equation and the velocity of the interface between the solid phase and the mushy zone,the front-tracking method is used to track the move of the interface.The interface between the liquid phase and the mushy zone is obtained by interpolation of the temperature field.Finally,four numerical examples are provided to assess the performance of the proposed numerical method.展开更多
In this work,a new methodology is presented to mainly solve the fluid–solid interaction(FSI)equation.This methodology combines the advantages of the Newmark precise integral method(NPIM)and the dual neural network(DN...In this work,a new methodology is presented to mainly solve the fluid–solid interaction(FSI)equation.This methodology combines the advantages of the Newmark precise integral method(NPIM)and the dual neural network(DNN)method.The NPIM is employed to modify the exponential matrix and loading vector based on the DNN integral method.This involves incorporating the basic assumption of the Newmark-βmethod into the dynamic equation and eliminating the acceleration term from the dynamic equilibrium equation.As a result,the equation is reduced to a first-order linear equation system.Subsequently,the PIM is applied to integrate the system step by step within the NPIM.The DNN method is adopted to solve the inhomogeneous term through fitting the integrand and the original function with a pair of neural networks,and the integral term is solved using the Newton–Leibniz formula.Numerical examples demonstrate that the proposed methodology significantly improves computing efficiency and provides sufficient precision compared to the DNN method.This is particularly evident when analyzing large-scale structures under blast loading conditions.展开更多
Based on the classical theory of thin plate and Biot theory, a precise model of the transverse vibrations of a thin rectangular porous plate is proposed. The first order differential equations of the porous plate are ...Based on the classical theory of thin plate and Biot theory, a precise model of the transverse vibrations of a thin rectangular porous plate is proposed. The first order differential equations of the porous plate are derived in the frequency domain. By considering the coupling effect between the solid phase and the fluid phase and without any hypothesis for the fluid displacement, the model presented here is rigorous and close to the real materials. Owing to the use of extended homogeneous capacity precision integration method and precise element method, the model can be applied in higher frequency range than pure numerical methods. This model also easily adapts to various boundary conditions. Numerical results are given for two different porous plates under different excitations and boundary conditions.展开更多
基金supported by the National Natural Science Foundation of China (Nos. 10902020 and 10721062)
文摘A fast precise integration method is developed for the time integral of the hyperbolic heat conduction problem. The wave nature of heat transfer is used to analyze the structure of the matrix exponential, leading to the fact that the matrix exponential is sparse. The presented method employs the sparsity of the matrix exponential to improve the original precise integration method. The merits are that the proposed method is suit- able for large hyperbolic heat equations and inherits the accuracy of the original version and the good computational efficiency, which are verified by two numerical examples.
基金supported by the National Natural Science Foundation of China (Nos 40974066 and 40821062)National Basic Research Program of China (No 2007CB209602)
文摘General purpose graphic processing unit(GPU) calculation technology is gradually widely used in various fields.Its mode of single instruction,multiple threads is capable of seismic numerical simulation which has a huge quantity of data and calculation steps.In this study,we introduce a GPU-based parallel calculation method of a precise integration method(PIM) for seismic forward modeling.Compared with CPU single-core calculation,GPU parallel calculating perfectly keeps the features of PIM,which has small bandwidth,high accuracy and capability of modeling complex substructures,and GPU calculation brings high computational efficiency,which means that high-performing GPU parallel calculation can make seismic forward modeling closer to real seismic records.
基金the National Natural Science Foundation of China(No.11672064)。
文摘A radial integral boundary element method(BEM)is used to simulate the phase change problem with a mushy zone in this paper.Three phases,including the solid phase,the liquid phase,and the mushy zone,are considered in the phase change problem.First,according to the continuity conditions of temperature and its gradient on the liquid-mushy interface,the mushy zone and the liquid phase in the simulation can be considered as a whole part,namely,the non-solid phase,and the change of latent heat is approximated by heat source which is dependent on temperature.Then,the precise integration BEM is used to obtain the differential equations in the solid phase zone and the non-solid phase zone,respectively.Moreover,an iterative predictor-corrector precise integration method(PIM)is needed to solve the differential equations and obtain the temperature field and the heat flux on the boundary.According to an energy balance equation and the velocity of the interface between the solid phase and the mushy zone,the front-tracking method is used to track the move of the interface.The interface between the liquid phase and the mushy zone is obtained by interpolation of the temperature field.Finally,four numerical examples are provided to assess the performance of the proposed numerical method.
基金supported by the National Natural Science Foundation of China(Grant Nos.12072288,U2241274,and 12272319).
文摘In this work,a new methodology is presented to mainly solve the fluid–solid interaction(FSI)equation.This methodology combines the advantages of the Newmark precise integral method(NPIM)and the dual neural network(DNN)method.The NPIM is employed to modify the exponential matrix and loading vector based on the DNN integral method.This involves incorporating the basic assumption of the Newmark-βmethod into the dynamic equation and eliminating the acceleration term from the dynamic equilibrium equation.As a result,the equation is reduced to a first-order linear equation system.Subsequently,the PIM is applied to integrate the system step by step within the NPIM.The DNN method is adopted to solve the inhomogeneous term through fitting the integrand and the original function with a pair of neural networks,and the integral term is solved using the Newton–Leibniz formula.Numerical examples demonstrate that the proposed methodology significantly improves computing efficiency and provides sufficient precision compared to the DNN method.This is particularly evident when analyzing large-scale structures under blast loading conditions.
基金Project supported by the National Natural Science Foundation of China(nos.11162001,11502056 and 51665006)
文摘Based on the classical theory of thin plate and Biot theory, a precise model of the transverse vibrations of a thin rectangular porous plate is proposed. The first order differential equations of the porous plate are derived in the frequency domain. By considering the coupling effect between the solid phase and the fluid phase and without any hypothesis for the fluid displacement, the model presented here is rigorous and close to the real materials. Owing to the use of extended homogeneous capacity precision integration method and precise element method, the model can be applied in higher frequency range than pure numerical methods. This model also easily adapts to various boundary conditions. Numerical results are given for two different porous plates under different excitations and boundary conditions.