The time accuracy of the exponentially accurate Fourier time spectral method(TSM) is examined and compared with a conventional 2nd-order backward difference formula(BDF) method for periodic unsteady flows. In particul...The time accuracy of the exponentially accurate Fourier time spectral method(TSM) is examined and compared with a conventional 2nd-order backward difference formula(BDF) method for periodic unsteady flows. In particular, detailed error analysis based on numerical computations is performed on the accuracy of resolving the local pressure coefficient and global integrated force coefficients for smooth subsonic and non-smooth transonic flows with moving shock waves on a pitching airfoil. For smooth subsonic flows, the Fourier TSM method offers a significant accuracy advantage over the BDF method for the prediction of both the local pressure coefficient and integrated force coefficients. For transonic flows where the motion of the discontinuous shock wave contributes significant higherorder harmonic contents to the local pressure fluctuations,a sufficient number of modes must be included before the Fourier TSM provides an advantage over the BDF method.The Fourier TSM, however, still offers better accuracy than the BDF method for integrated force coefficients even for transonic flows. A problem of non-symmetric solutions for symmetric periodic flows due to the use of odd numbers of intervals is uncovered and analyzed. A frequency-searching method is proposed for problems where the frequency is not known a priori. The method is tested on the vortex shedding problem of the flow over a circular cylinder.展开更多
The ti me dependent flow of upper-convected Maxwell fluid in a horizontal circular pip e is studied by spectral method. The time dependent problem is mathematically re duced to a partial differential equation of secon...The ti me dependent flow of upper-convected Maxwell fluid in a horizontal circular pip e is studied by spectral method. The time dependent problem is mathematically re duced to a partial differential equation of second order. By using spectral meth od the partial differential equation can be reduced to a system of ordinary diff erential equations for different terms of Chebyshev polynomials approximations. The ordinary differential equations are solved by Laplace transform and the eige nvalue method that leads to an analytical form of the solutions.展开更多
A time-spectral method for solution of initial value partial differential equations is outlined. Multivariate Chebyshev series are used to represent all temporal, spatial and physical parameter domains in this general...A time-spectral method for solution of initial value partial differential equations is outlined. Multivariate Chebyshev series are used to represent all temporal, spatial and physical parameter domains in this generalized weighted residual method (GWRM). The approximate solutions obtained are thus analytical, finite order multivariate polynomials. The method avoids time step limitations. To determine the spectral coefficients, a system of algebraic equations is solved iteratively. A root solver, with excellent global convergence properties, has been developed. Accuracy and efficiency are controlled by the number of included Chebyshev modes and by use of temporal and spatial subdomains. As examples of advanced application, stability problems within ideal and resistive magnetohydrodynamics (MHD) are solved. To introduce the method, solutions to a stiff ordinary differential equation are demonstrated and discussed. Subsequently, the GWRM is applied to the Burger and forced wave equations. Comparisons with the explicit Lax-Wendroff and implicit Crank-Nicolson finite difference methods show that the method is accurate and efficient. Thus the method shows potential for advanced initial value problems in fluid mechanics and MHD.展开更多
Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the ...Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the numerical methods for it. Recently, since the development of infinite dimensional dynamic system the dynamical behavior of NSE has been investigated. The paper [1] studied the long time wellposedness, the existence of universal attractor and the estimate of Lyapunov exponent for NSE with weakly damped. At the same time it was need to study the large time new computational methods and to discuss its convergence error estimate, the existence of approximate attractors etc. In this pape we study the NSE with weakly damped (1.1). We assume,where 0【λ【2 is a constant. If we wish to construct the higher accuracy computational scheme, it will be difficult that staigh from the equation (1.1). Therefore we start with (1. 4) and use fully discrete Fourier spectral method with time difference 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...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.展开更多
A numerical scheme is developed to extend the scope of the spectral method without solving the covariant and contravariant forms of the Navier-Stokes equations in the curvilinear coordinates. The primitive variables a...A numerical scheme is developed to extend the scope of the spectral method without solving the covariant and contravariant forms of the Navier-Stokes equations in the curvilinear coordinates. The primitive variables are represented by the Fourier series and the Chebyshev polynomials in the computational space. The time advancement is accomplished by a high-order time-splitting method,and a corresponding high-order pressure condition at the wall is introduced to reduce the splitting error. Compared with the previous pseudo-spectral scheme,in which the Navier-Stokes equations are solved in the covariant and contravariant forms,the present scheme reduces the computational cost and,at the same time,keeps the spectral accuracy. The scheme is tested in the simulations of the turbulent flow in a channel with a static streamwise wavy wall and the turbulent flow over a flexible wall undergoing the streamwise traveling wave motion. The turbulent flow over an oscillating dimple is studied with the present numerical scheme,and the periodic generation of the vortical structures is analyzed.展开更多
The dynamics of two non-coupled qubits independently interacting with their reservoirs is solved by the time convolutionless projection operator method. We study two-qubit quantum correlation dynamics for two differen...The dynamics of two non-coupled qubits independently interacting with their reservoirs is solved by the time convolutionless projection operator method. We study two-qubit quantum correlation dynamics for two different types of spectral densities, which are a Lorentzian distribution and an Ohmic spectral density with a Lorentzian-Drude cutoff function. For two qubits initially prepared in the initial Bell state, quantum discord can keep longer time and reach larger values in nonMarkovian reservoirs for the first spectral distribution or by reducing the cutoff frequency for the second case. For the initial Bell-like state, the dynamic behaviors of quantum discord and entanglement are compared. The results show that a long time of quantum correlation can be obtained by adjusting some parameters in experiment and further confirm that the discord can capture quantum correlation in addition to entanglement.展开更多
In this paper we use the spectral method to analyse the generalized KuramotoSiva-shinsky equations. We prove the ealstence and uniqueness of global smooth solution of the equations. Finally, we obtain the error estima...In this paper we use the spectral method to analyse the generalized KuramotoSiva-shinsky equations. We prove the ealstence and uniqueness of global smooth solution of the equations. Finally, we obtain the error estimation between spectral approkimate solution and exact solution on large time.展开更多
基金supported by the State Scholarship Fund of the China Scholarship Council (Grant 2009629129)
文摘The time accuracy of the exponentially accurate Fourier time spectral method(TSM) is examined and compared with a conventional 2nd-order backward difference formula(BDF) method for periodic unsteady flows. In particular, detailed error analysis based on numerical computations is performed on the accuracy of resolving the local pressure coefficient and global integrated force coefficients for smooth subsonic and non-smooth transonic flows with moving shock waves on a pitching airfoil. For smooth subsonic flows, the Fourier TSM method offers a significant accuracy advantage over the BDF method for the prediction of both the local pressure coefficient and integrated force coefficients. For transonic flows where the motion of the discontinuous shock wave contributes significant higherorder harmonic contents to the local pressure fluctuations,a sufficient number of modes must be included before the Fourier TSM provides an advantage over the BDF method.The Fourier TSM, however, still offers better accuracy than the BDF method for integrated force coefficients even for transonic flows. A problem of non-symmetric solutions for symmetric periodic flows due to the use of odd numbers of intervals is uncovered and analyzed. A frequency-searching method is proposed for problems where the frequency is not known a priori. The method is tested on the vortex shedding problem of the flow over a circular cylinder.
基金SupportedbytheNationalNaturalScienceFoundation( No .19672 0 63)andbythekeyprojectoftheStateNationalitiesAffairsCommissionofChina(No .990 5 ) .
文摘The ti me dependent flow of upper-convected Maxwell fluid in a horizontal circular pip e is studied by spectral method. The time dependent problem is mathematically re duced to a partial differential equation of second order. By using spectral meth od the partial differential equation can be reduced to a system of ordinary diff erential equations for different terms of Chebyshev polynomials approximations. The ordinary differential equations are solved by Laplace transform and the eige nvalue method that leads to an analytical form of the solutions.
文摘A time-spectral method for solution of initial value partial differential equations is outlined. Multivariate Chebyshev series are used to represent all temporal, spatial and physical parameter domains in this generalized weighted residual method (GWRM). The approximate solutions obtained are thus analytical, finite order multivariate polynomials. The method avoids time step limitations. To determine the spectral coefficients, a system of algebraic equations is solved iteratively. A root solver, with excellent global convergence properties, has been developed. Accuracy and efficiency are controlled by the number of included Chebyshev modes and by use of temporal and spatial subdomains. As examples of advanced application, stability problems within ideal and resistive magnetohydrodynamics (MHD) are solved. To introduce the method, solutions to a stiff ordinary differential equation are demonstrated and discussed. Subsequently, the GWRM is applied to the Burger and forced wave equations. Comparisons with the explicit Lax-Wendroff and implicit Crank-Nicolson finite difference methods show that the method is accurate and efficient. Thus the method shows potential for advanced initial value problems in fluid mechanics and MHD.
文摘Nonlinear Schrodinger equation (NSE) arises in many physical problems. It is a very important equation. A lot of works studied the wellposed, the existence of solution of NSE etc. And there are many works studied the numerical methods for it. Recently, since the development of infinite dimensional dynamic system the dynamical behavior of NSE has been investigated. The paper [1] studied the long time wellposedness, the existence of universal attractor and the estimate of Lyapunov exponent for NSE with weakly damped. At the same time it was need to study the large time new computational methods and to discuss its convergence error estimate, the existence of approximate attractors etc. In this pape we study the NSE with weakly damped (1.1). We assume,where 0【λ【2 is a constant. If we wish to construct the higher accuracy computational scheme, it will be difficult that staigh from the equation (1.1). Therefore we start with (1. 4) and use fully discrete Fourier spectral method with time difference 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 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.
基金supported by the National Natural Science Foundation of China (No. 10772098)
文摘A numerical scheme is developed to extend the scope of the spectral method without solving the covariant and contravariant forms of the Navier-Stokes equations in the curvilinear coordinates. The primitive variables are represented by the Fourier series and the Chebyshev polynomials in the computational space. The time advancement is accomplished by a high-order time-splitting method,and a corresponding high-order pressure condition at the wall is introduced to reduce the splitting error. Compared with the previous pseudo-spectral scheme,in which the Navier-Stokes equations are solved in the covariant and contravariant forms,the present scheme reduces the computational cost and,at the same time,keeps the spectral accuracy. The scheme is tested in the simulations of the turbulent flow in a channel with a static streamwise wavy wall and the turbulent flow over a flexible wall undergoing the streamwise traveling wave motion. The turbulent flow over an oscillating dimple is studied with the present numerical scheme,and the periodic generation of the vortical structures is analyzed.
基金Project supported by the National Natural Science Foundation of China (Grant Nos. 11264011 and 11104113)the Natural Science Foundation of Hunan Province, China (Grant Nos. 13JJ6059 and 11JJ6007)the Natural Science Foundation of Education Department of Hunan Province, China (GrantNo. 11C1057)
文摘The dynamics of two non-coupled qubits independently interacting with their reservoirs is solved by the time convolutionless projection operator method. We study two-qubit quantum correlation dynamics for two different types of spectral densities, which are a Lorentzian distribution and an Ohmic spectral density with a Lorentzian-Drude cutoff function. For two qubits initially prepared in the initial Bell state, quantum discord can keep longer time and reach larger values in nonMarkovian reservoirs for the first spectral distribution or by reducing the cutoff frequency for the second case. For the initial Bell-like state, the dynamic behaviors of quantum discord and entanglement are compared. The results show that a long time of quantum correlation can be obtained by adjusting some parameters in experiment and further confirm that the discord can capture quantum correlation in addition to entanglement.
文摘In this paper we use the spectral method to analyse the generalized KuramotoSiva-shinsky equations. We prove the ealstence and uniqueness of global smooth solution of the equations. Finally, we obtain the error estimation between spectral approkimate solution and exact solution on large time.
