Time-spectral solution of ordinary and partial differential equations is often regarded as an inefficient approach. The associated extension of the time domain, as compared to finite difference methods, is believed to...Time-spectral solution of ordinary and partial differential equations is often regarded as an inefficient approach. The associated extension of the time domain, as compared to finite difference methods, is believed to result in uncomfortably many numerical operations and high memory requirements. It is shown in this work that performance is substantially enhanced by the introduction of algorithms for temporal and spatial subdomains in combination with sparse matrix methods. The accuracy and efficiency of the recently developed time spectral, generalized weighted residual method (GWRM) are compared to that of the explicit Lax-Wendroff and implicit Crank-Nicolson methods. Three initial-value PDEs are employed as model problems;the 1D Burger equation, a forced 1D wave equation and a coupled system of 14 linearized ideal magnetohydrodynamic (MHD) equations. It is found that the GWRM is more efficient than the time-stepping methods at high accuracies. The advantageous scalings Nt<sup style="margin-left:-6px;">1.0Ns<sup style="margin-left:-6px;">1.43 and Nt<sup style="margin-left:-6px;">0.0Ns<sup style="margin-left:-6px;">1.08 were obtained for CPU time and memory requirements, respectively, with Nt and Ns denoting the number of temporal and spatial subdomains. For time-averaged solution of the two-time-scales forced wave equation, GWRM performance exceeds that of the finite difference methods by an order of magnitude both in terms of CPU time and memory requirement. Favorable subdomain scaling is demonstrated for the MHD equations, indicating a potential for efficient solution of advanced initial-value problems in, for example, fluid mechanics and MHD.展开更多
Temporal and spatial subdomain techniques are proposed for a time-spectral method for solution of initial-value problems. The spectral method, called the generalised weighted residual method (GWRM), is a generalisatio...Temporal and spatial subdomain techniques are proposed for a time-spectral method for solution of initial-value problems. The spectral method, called the generalised weighted residual method (GWRM), is a generalisation of weighted residual methods to the time and parameter domains [1]. A semi-analytical Chebyshev polynomial ansatz is employed, and the problem reduces to determine the coefficients of the ansatz from linear or nonlinear algebraic systems of equations. In order to avoid large memory storage and computational cost, it is preferable to subdivide the temporal and spatial domains into subdomains. Methods and examples of this article demonstrate how this can be achieved.展开更多
文摘Time-spectral solution of ordinary and partial differential equations is often regarded as an inefficient approach. The associated extension of the time domain, as compared to finite difference methods, is believed to result in uncomfortably many numerical operations and high memory requirements. It is shown in this work that performance is substantially enhanced by the introduction of algorithms for temporal and spatial subdomains in combination with sparse matrix methods. The accuracy and efficiency of the recently developed time spectral, generalized weighted residual method (GWRM) are compared to that of the explicit Lax-Wendroff and implicit Crank-Nicolson methods. Three initial-value PDEs are employed as model problems;the 1D Burger equation, a forced 1D wave equation and a coupled system of 14 linearized ideal magnetohydrodynamic (MHD) equations. It is found that the GWRM is more efficient than the time-stepping methods at high accuracies. The advantageous scalings Nt<sup style="margin-left:-6px;">1.0Ns<sup style="margin-left:-6px;">1.43 and Nt<sup style="margin-left:-6px;">0.0Ns<sup style="margin-left:-6px;">1.08 were obtained for CPU time and memory requirements, respectively, with Nt and Ns denoting the number of temporal and spatial subdomains. For time-averaged solution of the two-time-scales forced wave equation, GWRM performance exceeds that of the finite difference methods by an order of magnitude both in terms of CPU time and memory requirement. Favorable subdomain scaling is demonstrated for the MHD equations, indicating a potential for efficient solution of advanced initial-value problems in, for example, fluid mechanics and MHD.
文摘Temporal and spatial subdomain techniques are proposed for a time-spectral method for solution of initial-value problems. The spectral method, called the generalised weighted residual method (GWRM), is a generalisation of weighted residual methods to the time and parameter domains [1]. A semi-analytical Chebyshev polynomial ansatz is employed, and the problem reduces to determine the coefficients of the ansatz from linear or nonlinear algebraic systems of equations. In order to avoid large memory storage and computational cost, it is preferable to subdivide the temporal and spatial domains into subdomains. Methods and examples of this article demonstrate how this can be achieved.
