The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integ...The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integration schemes derived by the singlepoint precise integration method given in this paper are proved unconditionally stable.Comparisons between the schemes derived by the finite difference method and theschemes by the method employed in the present paper are made for diffusion andconvective-diffusion equations. Nunierical examples show the superiority of the singlepoint integration method.展开更多
A balanced adaptive time-stepping strategy is implemented in an implicit discontinuous Galerkin solver to guarantee the temporal accuracy of unsteady simulations.A proper relation between the spatial,temporal and iter...A balanced adaptive time-stepping strategy is implemented in an implicit discontinuous Galerkin solver to guarantee the temporal accuracy of unsteady simulations.A proper relation between the spatial,temporal and iterative errors generated within one time step is constructed.With an estimate of temporal and spatial error using an embedded RungeKutta scheme and a higher order spatial discretization,an adaptive time-stepping strategy is proposed based on the idea that the time step should be the maximum without obviously infuencing the total error of the discretization.The designed adaptive time-stepping strategy is then tested in various types of problems including isentropic vortex convection,steady-state fow past a fat plate,Taylor-Green vortex and turbulent fow over a circular cylinder at Re=3900.The results indicate that the adaptive time-stepping strategy can maintain that the discretization error is dominated by the spatial error and relatively high efciency is obtained for unsteady and steady,well-resolved and under-resolved simulations.展开更多
The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficu...The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method (e.g., an algorithm with an arbitrary order of accuracy) is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.展开更多
提出一种基于精细积分法与时域微分求积法相结合的传输线方程的数值求解方法。首先将传输线方程采用基于紧致有限差分法的四阶差分格式进行空间离散,得到关于时间的一阶线性常微分方程组,四阶差分格式对于空间微分有很好的近似精度。然...提出一种基于精细积分法与时域微分求积法相结合的传输线方程的数值求解方法。首先将传输线方程采用基于紧致有限差分法的四阶差分格式进行空间离散,得到关于时间的一阶线性常微分方程组,四阶差分格式对于空间微分有很好的近似精度。然后利用精细积分法与微分求积法对一阶线性常微分方程组进行数值求解。通过理论分析可知,与传统的传输线方程数值求解方法——时域有限差分法(Finite difference time domain,FDTD)相比,所提方法不涉及到状态矩阵求逆运算,保证了数值求解精度,并且其数值稳定性与计算时间、空间步长无关,可采用大步长进行数值计算,能够有效提高计算效率。最后利用仿真实例进行算法验证,结果显示,相比于时域有限差分法,所提方法能够抑制数值振荡,提高了计算精度。展开更多
文摘The numerical time step integrations of PDEs are mainly carried out by the finitedifference method to date. However,when the time step becomes longer, it causes theproblem of numerical instability,. The explicit integration schemes derived by the singlepoint precise integration method given in this paper are proved unconditionally stable.Comparisons between the schemes derived by the finite difference method and theschemes by the method employed in the present paper are made for diffusion andconvective-diffusion equations. Nunierical examples show the superiority of the singlepoint integration method.
基金Zhen-Guo Yan acknowledges supports from the National Natural Science Foundation of China(Grant no.11902344)National Numerical Windtunnel Project.The development of the implicit solver in Nektar++has been supported by EPSRC grant(EP/R029423/1)UK Turbulence Consortium grant(EP/R029326/1).
文摘A balanced adaptive time-stepping strategy is implemented in an implicit discontinuous Galerkin solver to guarantee the temporal accuracy of unsteady simulations.A proper relation between the spatial,temporal and iterative errors generated within one time step is constructed.With an estimate of temporal and spatial error using an embedded RungeKutta scheme and a higher order spatial discretization,an adaptive time-stepping strategy is proposed based on the idea that the time step should be the maximum without obviously infuencing the total error of the discretization.The designed adaptive time-stepping strategy is then tested in various types of problems including isentropic vortex convection,steady-state fow past a fat plate,Taylor-Green vortex and turbulent fow over a circular cylinder at Re=3900.The results indicate that the adaptive time-stepping strategy can maintain that the discretization error is dominated by the spatial error and relatively high efciency is obtained for unsteady and steady,well-resolved and under-resolved simulations.
基金financial support from Hunan Provincial Natura1 Science Foundation of China,Grant Number:02JJY2085,for this study
文摘The precise time step integration method proposed for linear time-invariant homogeneous dynamic systems can provide precise numerical results that approach an exact solution at the integration points. However, difficulty arises when the algorithm is used for non-homogeneous dynamic systems, due to the inverse matrix calculation and the simulation accuracy of the applied loading. By combining the Gaussian quadrature method and state space theory with the calculation technique of matrix exponential function in the precise time step integration method, a new modified precise time step integration method (e.g., an algorithm with an arbitrary order of accuracy) is proposed. In the new method, no inverse matrix calculation or simulation of the applied loading is needed, and the computing efficiency is improved. In particular, the proposed method is independent of the quality of the matrix H. If the matrix H is singular or nearly singular, the advantage of the method is remarkable. The numerical stability of the proposed algorithm is discussed and a numerical example is given to demonstrate the validity and efficiency of the algorithm.
文摘提出一种基于精细积分法与时域微分求积法相结合的传输线方程的数值求解方法。首先将传输线方程采用基于紧致有限差分法的四阶差分格式进行空间离散,得到关于时间的一阶线性常微分方程组,四阶差分格式对于空间微分有很好的近似精度。然后利用精细积分法与微分求积法对一阶线性常微分方程组进行数值求解。通过理论分析可知,与传统的传输线方程数值求解方法——时域有限差分法(Finite difference time domain,FDTD)相比,所提方法不涉及到状态矩阵求逆运算,保证了数值求解精度,并且其数值稳定性与计算时间、空间步长无关,可采用大步长进行数值计算,能够有效提高计算效率。最后利用仿真实例进行算法验证,结果显示,相比于时域有限差分法,所提方法能够抑制数值振荡,提高了计算精度。