In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navies-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinea...In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navies-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinear system on the coarse mesh space and two similar linear systems (with same stiffness matrix but different right-hand side) on the fine mesh space. The convergence analysis and error estimation of the algorithm are given for the case of conforming elements. Furthermore, the Mgorithm produces a numerical solution with the optimal asymptotic H^2-error. Finally, we give a numerical illustration to demonstrate the effectiveness of the two-grid algorithm for solving the Navier-Stokes equations.展开更多
A new form function involving parameters Pi is presented. On the basis of the form function, an initial form of tension structure was Sound by interpolating through the control points on boundary of the structure. The...A new form function involving parameters Pi is presented. On the basis of the form function, an initial form of tension structure was Sound by interpolating through the control points on boundary of the structure. The form function can be controlled by changing beta (1) according to the pre-tension and the boundary of the structure. The final form of a tension structure should be an equilibrium system under the pretension. To examine the nature of the initial form, the FEM was used. Many examples show that the initial form gives a very ideal result for equal or unequal pre-tension in two directions of the structure. In general cases, there is little difference between the initial form and the final one.展开更多
In order to solve fuzzy mathematical programming with soft constraints,the initial models must first be converted into crisp models.Membership functions are employed to describe the fuzzy right-hand side parameters ne...In order to solve fuzzy mathematical programming with soft constraints,the initial models must first be converted into crisp models.Membership functions are employed to describe the fuzzy right-hand side parameters needed to achieve this conversion.In some cases,echelon form membership functions(EFMFs)are required to depict the actual fuzzy situation.However,due to their discrete properties,fuzzy programming problems with such membership functions cannot be modeled by traditional methods.Motivated by these challenges,this paper introduces a novel absolute value representation modeling approach to formulate fuzzy programming using EFMFs.This approach can translate a discrete model to a continuous one which can then be easily solved.Finally,by means of a numerical example,the effectiveness of our new approach is demonstrated.展开更多
Residual based on a posteriori error estimates for conforming finite element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level met...Residual based on a posteriori error estimates for conforming finite element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level method were derived. The posteriori error estimates contained additional terms in comparison to the error estimates for the solution obtained by the standard finite element method. The importance of these additional terms in the error estimates was investigated by studying their asymptotic behavior. For optimal scaled meshes, these bounds are not of higher order than of convergence of discrete solution.展开更多
Based on the definition of a logic structure feature to relate logically functional requirements to geometric representation independent upon detailed geometric representation, this paper presents an idea of logical s...Based on the definition of a logic structure feature to relate logically functional requirements to geometric representation independent upon detailed geometric representation, this paper presents an idea of logical structure modeling for computer aided conceptual design and makes attempt to establish a representation formalism of logic structure modeling. The definition and representation of logical structure feature are given and an assembly module definition for supporting top down conceptual design is also proposed. The proposed scheme contributes to several aspects of conceptual design research, especially to provide elementarily a formal methodology for computer aided conceptual design system development and operation.展开更多
Hollow cylinders are widely used in spacecraft, rockets, weapons, metallurgy, materials, and mechanical manufacturing industries, and so on, hydraulic bulging roll cylinder and hydraulic press work all belong to hollo...Hollow cylinders are widely used in spacecraft, rockets, weapons, metallurgy, materials, and mechanical manufacturing industries, and so on, hydraulic bulging roll cylinder and hydraulic press work all belong to hollow cylinders. However, up till now, the solution of the cylinder subjected to the pressures in the three-dimensional space is still at the stage of the analytical solution to the normal pressure or the approximate solution to the variable pressure by numerical method. The analytical solution to the variable pressure of the cylinder has not yet made any breakthrough in theory and can not meet accurate theoretical analysis and calculation requirements of the cylindrical in Engineering. In view of their importance, the precision calculation and theoretical analysis are required to investigate on engineering. A stress function which meets both the biharmonic equations and boundary conditions is constructed in the three-dimensional space. Furthermore, the analytic solution of a hollow cylinder subjected to exponential function distributed variable pressure on its inner and outer surfaces is deduced. By controlling the pressure subject to exponential function distributed variable pressure in the hydraulic bulging roller without any rolling load, using a static tester to record the strain supported hydraulic bulging roll, and comparing with the theoretical calculation, the experimental test result has a higher degree of agreement with the theoretical calculation. Simultaneously, the famous Lam6 solution can be deduced when given the unlimited length of cylinder along the axis. The analytic solution paves the way for the mathematic building and solution of hollow cylinder with randomly uneven pressure.展开更多
This study proposes an equivalent-elevation method to evaluate the integrated effects of latitude and elevation on regional and local-scale permafrost distribution in the Qinghai-Tibet Plateau,and to model the general...This study proposes an equivalent-elevation method to evaluate the integrated effects of latitude and elevation on regional and local-scale permafrost distribution in the Qinghai-Tibet Plateau,and to model the general permafrost-distribution patterns in regional and local-scale area.It is found that the Gaussian curve―an empirical model describing the relation between variations of altitudinal permafrost lower limit (PLL) and latitude in the Northern Hemisphere―could be applied in regional-and local-scale areas in the Qinghai-Tibet Plateau in a latitude-sensitive interval of 30°-50°N.The curve was then used to evaluate the latitudinal effect on permafrost distribution through transforming the latitudinal effect into a kind of altitudinal difference of PLL.This study then calculated the local equivalent-elevation value by overlaying the altitudinal difference of PLL onto real elevation at a certain location.The equivalent-elevation method was verified in an experimental subwatershed of the Qinghai-Tibet Plateau.However,feasibility of the method should be further tested in order to extend for future studies.The use of equivalent-elevation values can build a platform for comparing the regional general permafrost distribution in the plateau,and for basing further evaluations of local factors' effects on regional permafrost distribution.展开更多
This paper presents equations for estimating limiting stand density for Z undulata plantations grown in hot desert areas of Raj asthan State in India. Five different stand level basal area projection models, belonging...This paper presents equations for estimating limiting stand density for Z undulata plantations grown in hot desert areas of Raj asthan State in India. Five different stand level basal area projection models, belonging to the path invariant algebraic difference form of a non-linear growth function, were also tested and compared. These models can be used to predict future basal area as a function of stand variables like dominant height and stem number per hectare and are necessary for reviewing different silvicultural treatment options. Data from 22 sample plots were used for modelling. An all possible growth intervals data structure was used. Both, qualitative and quantitative criteria were used to compare alternative models. The Akaike's information criteria differ- ence statistic was used to analyze the predictive ability of the models. Results show that the model proposed by Hui and Gadow performed best and hence this model is recommended for use in predicting basal area development in 12 undulata plantations in the study area. The data used were not from thinned stands, and hence the models may be less accurate when used for predictions when natural mortality is very significant.展开更多
A new method is presented to generate two-directional (2D) grid multi-scroll chaotic attractors via a specific form of the sine function and sign function series, which are applied to increase saddle points of index 2...A new method is presented to generate two-directional (2D) grid multi-scroll chaotic attractors via a specific form of the sine function and sign function series, which are applied to increase saddle points of index 2. The scroll number in the x-direction is modified easily through changing the thresholds of the specific form of the sine function, while the scroll number in the y-direction is controlled by the sign function series. Some basic dynamical properties, such as equilibrium points, bifurcation diagram, phase portraits, and Lyapunov exponents spectrum are studied. Furthermore, the electronic circuit of the system is designed and its simulation results are given by Multisim 10.展开更多
This article studies on Cauchy’s function f (z) and its integral, (2πi)J[ f (z)] ≡ ■C f (t)dt/(t z) taken along a closed simple contour C, in regard to their comprehensive properties over the entire z =...This article studies on Cauchy’s function f (z) and its integral, (2πi)J[ f (z)] ≡ ■C f (t)dt/(t z) taken along a closed simple contour C, in regard to their comprehensive properties over the entire z = x + iy plane consisted of the simply connected open domain D + bounded by C and the open domain D outside C. (1) With f (z) assumed to be C n (n ∞-times continuously differentiable) z ∈ D + and in a neighborhood of C, f (z) and its derivatives f (n) (z) are proved uniformly continuous in the closed domain D + = [D + + C]. (2) Cauchy’s integral formulas and their derivatives z ∈ D + (or z ∈ D ) are proved to converge uniformly in D + (or in D = [D +C]), respectively, thereby rendering the integral formulas valid over the entire z-plane. (3) The same claims (as for f (z) and J[ f (z)]) are shown extended to hold for the complement function F(z), defined to be C n z ∈ D and about C. (4) The uniform convergence theorems for f (z) and F(z) shown for arbitrary contour C are adapted to find special domains in the upper or lower half z-planes and those inside and outside the unit circle |z| = 1 such that the four general- ized Hilbert-type integral transforms are proved. (5) Further, the singularity distribution of f (z) in D is elucidated by considering the direct problem exemplified with several typ- ical singularities prescribed in D . (6) A comparative study is made between generalized integral formulas and Plemelj’s formulas on their differing basic properties. (7) Physical sig- nificances of these formulas are illustrated with applicationsto nonlinear airfoil theory. (8) Finally, an unsolved inverse problem to determine all the singularities of Cauchy function f (z) in domain D , based on the continuous numerical value of f (z) z ∈ D + = [D + + C], is presented for resolution as a conjecture.展开更多
In this paper a kind of theta function is constructed by means of spherical function. And we also obtain some Hilbert modular forms of half integral weight.
This paper is concerned with the problem of absolute stability for a control system with severalexecutive elements. Necessary and sufficient conditions are obtained for the existence of Liapunovfunction of Lur'e f...This paper is concerned with the problem of absolute stability for a control system with severalexecutive elements. Necessary and sufficient conditions are obtained for the existence of Liapunovfunction of Lur'e form with negative semi--definite derivative (i.e. V≤0).展开更多
This paper addresses the stability problem for a class of switched nonlinear time varying delay systems modeled by delay differential equations. By transforming the system representation under the arrow form and using...This paper addresses the stability problem for a class of switched nonlinear time varying delay systems modeled by delay differential equations. By transforming the system representation under the arrow form and using a new constructed Lyapunov function,the aggregation techniques,the Borne-Gentina practical stability criterion associated with the properties, new delay-independent stability conditions of the considered systems are established. Compared with the existing results in this area, the obtained result is explicit, simple to use and allows us to avoid the problem of searching a common Lyapunov function. Finally, an example is provided, with numerical simulations,to demonstrate the effectiveness of the proposed method.展开更多
For analyzing correlated binary data with high-dimensional covariates,we,in this paper,propose a two-stage shrinkage approach.First,we construct a weighted least-squares(WLS) type function using a special weighting sc...For analyzing correlated binary data with high-dimensional covariates,we,in this paper,propose a two-stage shrinkage approach.First,we construct a weighted least-squares(WLS) type function using a special weighting scheme on the non-conservative vector field of the generalized estimating equations(GEE) model.Second,we define a penalized WLS in the spirit of the adaptive LASSO for simultaneous variable selection and parameter estimation.The proposed procedure enjoys the oracle properties in high-dimensional framework where the number of parameters grows to infinity with the number of clusters.Moreover,we prove the consistency of the sandwich formula of the covariance matrix even when the working correlation matrix is misspecified.For the selection of tuning parameter,we develop a consistent penalized quadratic form(PQF) function criterion.The performance of the proposed method is assessed through a comparison with the existing methods and through an application to a crossover trial in a pain relief study.展开更多
This paper aims to inquiry on the dialectical relationships between Chinese architecture’s form and function in the horizon of long-term historical evolution,by clarifying how the histories of imperial main halls and...This paper aims to inquiry on the dialectical relationships between Chinese architecture’s form and function in the horizon of long-term historical evolution,by clarifying how the histories of imperial main halls and state ceremonies relate to and interact with each other during the Tang and Song dynasties.The imperialmain halls are themost archetypical and accomplished ones among traditional Chinese buildings,and their function is primarily ritual rather than practical.Although these magnificent buildings are indisputably considered the crowning achievement of Chinese architectural history,how the built space cooperates with ritual practices did not receive deserved attention.With the thorough reading of voluminous official ritual documents and comprehensive utilization of archaeological discoveries together with visual materials,the research presents a broad understanding of the asynchronous developments of the ground-level architectural arrangement and the ritual configuration across these centuries.Also specifically depicted are details of how architecture engaged in the ritual narration and was used as an apparatus to articulate emperor-courtier relationships.展开更多
The initial-boundary value problem of two-dimensional incompressible fluid flow in stream function form is considered. A prediction-correction Legendre spectral scheme is presented, which is easy to be performed. It ...The initial-boundary value problem of two-dimensional incompressible fluid flow in stream function form is considered. A prediction-correction Legendre spectral scheme is presented, which is easy to be performed. It is strictly proved that the numerical solution possesses the accuracy of second-order in time and higher order in space.展开更多
基金supported by National Foundation of Natural Science under the Grant 11071216
文摘In this paper, we propose a two-grid algorithm for solving the stream function formulation of the stationary Navies-Stokes equations. The algorithm is constructed by reducing the original system to one small, nonlinear system on the coarse mesh space and two similar linear systems (with same stiffness matrix but different right-hand side) on the fine mesh space. The convergence analysis and error estimation of the algorithm are given for the case of conforming elements. Furthermore, the Mgorithm produces a numerical solution with the optimal asymptotic H^2-error. Finally, we give a numerical illustration to demonstrate the effectiveness of the two-grid algorithm for solving the Navier-Stokes equations.
文摘A new form function involving parameters Pi is presented. On the basis of the form function, an initial form of tension structure was Sound by interpolating through the control points on boundary of the structure. The form function can be controlled by changing beta (1) according to the pre-tension and the boundary of the structure. The final form of a tension structure should be an equilibrium system under the pretension. To examine the nature of the initial form, the FEM was used. Many examples show that the initial form gives a very ideal result for equal or unequal pre-tension in two directions of the structure. In general cases, there is little difference between the initial form and the final one.
文摘In order to solve fuzzy mathematical programming with soft constraints,the initial models must first be converted into crisp models.Membership functions are employed to describe the fuzzy right-hand side parameters needed to achieve this conversion.In some cases,echelon form membership functions(EFMFs)are required to depict the actual fuzzy situation.However,due to their discrete properties,fuzzy programming problems with such membership functions cannot be modeled by traditional methods.Motivated by these challenges,this paper introduces a novel absolute value representation modeling approach to formulate fuzzy programming using EFMFs.This approach can translate a discrete model to a continuous one which can then be easily solved.Finally,by means of a numerical example,the effectiveness of our new approach is demonstrated.
文摘Residual based on a posteriori error estimates for conforming finite element solutions of incompressible Navier-Stokes equations with stream function form which were computed with seven recently proposed two-level method were derived. The posteriori error estimates contained additional terms in comparison to the error estimates for the solution obtained by the standard finite element method. The importance of these additional terms in the error estimates was investigated by studying their asymptotic behavior. For optimal scaled meshes, these bounds are not of higher order than of convergence of discrete solution.
文摘Based on the definition of a logic structure feature to relate logically functional requirements to geometric representation independent upon detailed geometric representation, this paper presents an idea of logical structure modeling for computer aided conceptual design and makes attempt to establish a representation formalism of logic structure modeling. The definition and representation of logical structure feature are given and an assembly module definition for supporting top down conceptual design is also proposed. The proposed scheme contributes to several aspects of conceptual design research, especially to provide elementarily a formal methodology for computer aided conceptual design system development and operation.
基金supported by National Natural Science Foundation of China (Grant No. 50875230)
文摘Hollow cylinders are widely used in spacecraft, rockets, weapons, metallurgy, materials, and mechanical manufacturing industries, and so on, hydraulic bulging roll cylinder and hydraulic press work all belong to hollow cylinders. However, up till now, the solution of the cylinder subjected to the pressures in the three-dimensional space is still at the stage of the analytical solution to the normal pressure or the approximate solution to the variable pressure by numerical method. The analytical solution to the variable pressure of the cylinder has not yet made any breakthrough in theory and can not meet accurate theoretical analysis and calculation requirements of the cylindrical in Engineering. In view of their importance, the precision calculation and theoretical analysis are required to investigate on engineering. A stress function which meets both the biharmonic equations and boundary conditions is constructed in the three-dimensional space. Furthermore, the analytic solution of a hollow cylinder subjected to exponential function distributed variable pressure on its inner and outer surfaces is deduced. By controlling the pressure subject to exponential function distributed variable pressure in the hydraulic bulging roller without any rolling load, using a static tester to record the strain supported hydraulic bulging roll, and comparing with the theoretical calculation, the experimental test result has a higher degree of agreement with the theoretical calculation. Simultaneously, the famous Lam6 solution can be deduced when given the unlimited length of cylinder along the axis. The analytic solution paves the way for the mathematic building and solution of hollow cylinder with randomly uneven pressure.
基金Under the auspices of Major State Basic Research Development Program of China(No.2010CB951402)National Natural Science Foundation of China(No.41101067)+1 种基金Foundation for Excellent Youth Scholars of Cold and Arid Regions Environmental and Engineering Research Institute,Chinese Academy of Sciences(No.Y184A91001)Research Program of State Key Laboratory of Frozen Soil Engineering of Cold and Arid Regions Environmental and Engineering Research Institute,Chinese Academy of Sciences(No.SKLFSE-ZQ-10)
文摘This study proposes an equivalent-elevation method to evaluate the integrated effects of latitude and elevation on regional and local-scale permafrost distribution in the Qinghai-Tibet Plateau,and to model the general permafrost-distribution patterns in regional and local-scale area.It is found that the Gaussian curve―an empirical model describing the relation between variations of altitudinal permafrost lower limit (PLL) and latitude in the Northern Hemisphere―could be applied in regional-and local-scale areas in the Qinghai-Tibet Plateau in a latitude-sensitive interval of 30°-50°N.The curve was then used to evaluate the latitudinal effect on permafrost distribution through transforming the latitudinal effect into a kind of altitudinal difference of PLL.This study then calculated the local equivalent-elevation value by overlaying the altitudinal difference of PLL onto real elevation at a certain location.The equivalent-elevation method was verified in an experimental subwatershed of the Qinghai-Tibet Plateau.However,feasibility of the method should be further tested in order to extend for future studies.The use of equivalent-elevation values can build a platform for comparing the regional general permafrost distribution in the plateau,and for basing further evaluations of local factors' effects on regional permafrost distribution.
基金the State Forest Department,Rajasthan for providing financial support for conducting this study and to their officials for rendering necessary assistance during fieldwork
文摘This paper presents equations for estimating limiting stand density for Z undulata plantations grown in hot desert areas of Raj asthan State in India. Five different stand level basal area projection models, belonging to the path invariant algebraic difference form of a non-linear growth function, were also tested and compared. These models can be used to predict future basal area as a function of stand variables like dominant height and stem number per hectare and are necessary for reviewing different silvicultural treatment options. Data from 22 sample plots were used for modelling. An all possible growth intervals data structure was used. Both, qualitative and quantitative criteria were used to compare alternative models. The Akaike's information criteria differ- ence statistic was used to analyze the predictive ability of the models. Results show that the model proposed by Hui and Gadow performed best and hence this model is recommended for use in predicting basal area development in 12 undulata plantations in the study area. The data used were not from thinned stands, and hence the models may be less accurate when used for predictions when natural mortality is very significant.
文摘A new method is presented to generate two-directional (2D) grid multi-scroll chaotic attractors via a specific form of the sine function and sign function series, which are applied to increase saddle points of index 2. The scroll number in the x-direction is modified easily through changing the thresholds of the specific form of the sine function, while the scroll number in the y-direction is controlled by the sign function series. Some basic dynamical properties, such as equilibrium points, bifurcation diagram, phase portraits, and Lyapunov exponents spectrum are studied. Furthermore, the electronic circuit of the system is designed and its simulation results are given by Multisim 10.
文摘This article studies on Cauchy’s function f (z) and its integral, (2πi)J[ f (z)] ≡ ■C f (t)dt/(t z) taken along a closed simple contour C, in regard to their comprehensive properties over the entire z = x + iy plane consisted of the simply connected open domain D + bounded by C and the open domain D outside C. (1) With f (z) assumed to be C n (n ∞-times continuously differentiable) z ∈ D + and in a neighborhood of C, f (z) and its derivatives f (n) (z) are proved uniformly continuous in the closed domain D + = [D + + C]. (2) Cauchy’s integral formulas and their derivatives z ∈ D + (or z ∈ D ) are proved to converge uniformly in D + (or in D = [D +C]), respectively, thereby rendering the integral formulas valid over the entire z-plane. (3) The same claims (as for f (z) and J[ f (z)]) are shown extended to hold for the complement function F(z), defined to be C n z ∈ D and about C. (4) The uniform convergence theorems for f (z) and F(z) shown for arbitrary contour C are adapted to find special domains in the upper or lower half z-planes and those inside and outside the unit circle |z| = 1 such that the four general- ized Hilbert-type integral transforms are proved. (5) Further, the singularity distribution of f (z) in D is elucidated by considering the direct problem exemplified with several typ- ical singularities prescribed in D . (6) A comparative study is made between generalized integral formulas and Plemelj’s formulas on their differing basic properties. (7) Physical sig- nificances of these formulas are illustrated with applicationsto nonlinear airfoil theory. (8) Finally, an unsolved inverse problem to determine all the singularities of Cauchy function f (z) in domain D , based on the continuous numerical value of f (z) z ∈ D + = [D + + C], is presented for resolution as a conjecture.
文摘In this paper a kind of theta function is constructed by means of spherical function. And we also obtain some Hilbert modular forms of half integral weight.
文摘This paper is concerned with the problem of absolute stability for a control system with severalexecutive elements. Necessary and sufficient conditions are obtained for the existence of Liapunovfunction of Lur'e form with negative semi--definite derivative (i.e. V≤0).
文摘This paper addresses the stability problem for a class of switched nonlinear time varying delay systems modeled by delay differential equations. By transforming the system representation under the arrow form and using a new constructed Lyapunov function,the aggregation techniques,the Borne-Gentina practical stability criterion associated with the properties, new delay-independent stability conditions of the considered systems are established. Compared with the existing results in this area, the obtained result is explicit, simple to use and allows us to avoid the problem of searching a common Lyapunov function. Finally, an example is provided, with numerical simulations,to demonstrate the effectiveness of the proposed method.
基金supported by National Natural Science Foundation of China(Grant No.11201306)the Innovation Program of Shanghai Municipal Education Commission(Grant No.13YZ065)+2 种基金the Fundamental Research Project of Shanghai Normal University(Grant No.SK201207)the scholarship under the State Scholarship Fund by the China Scholarship Council in 2011the Research Grant Council of Hong Kong, Hong Kong,China(Grant No.#HKBU2028/10P)
文摘For analyzing correlated binary data with high-dimensional covariates,we,in this paper,propose a two-stage shrinkage approach.First,we construct a weighted least-squares(WLS) type function using a special weighting scheme on the non-conservative vector field of the generalized estimating equations(GEE) model.Second,we define a penalized WLS in the spirit of the adaptive LASSO for simultaneous variable selection and parameter estimation.The proposed procedure enjoys the oracle properties in high-dimensional framework where the number of parameters grows to infinity with the number of clusters.Moreover,we prove the consistency of the sandwich formula of the covariance matrix even when the working correlation matrix is misspecified.For the selection of tuning parameter,we develop a consistent penalized quadratic form(PQF) function criterion.The performance of the proposed method is assessed through a comparison with the existing methods and through an application to a crossover trial in a pain relief study.
基金This study was supported by the National Natural Science Foundation of Chinaunder Grant 52038007 and the National Social Science Fund of China under Grant 20vmz008.
文摘This paper aims to inquiry on the dialectical relationships between Chinese architecture’s form and function in the horizon of long-term historical evolution,by clarifying how the histories of imperial main halls and state ceremonies relate to and interact with each other during the Tang and Song dynasties.The imperialmain halls are themost archetypical and accomplished ones among traditional Chinese buildings,and their function is primarily ritual rather than practical.Although these magnificent buildings are indisputably considered the crowning achievement of Chinese architectural history,how the built space cooperates with ritual practices did not receive deserved attention.With the thorough reading of voluminous official ritual documents and comprehensive utilization of archaeological discoveries together with visual materials,the research presents a broad understanding of the asynchronous developments of the ground-level architectural arrangement and the ritual configuration across these centuries.Also specifically depicted are details of how architecture engaged in the ritual narration and was used as an apparatus to articulate emperor-courtier relationships.
文摘The initial-boundary value problem of two-dimensional incompressible fluid flow in stream function form is considered. A prediction-correction Legendre spectral scheme is presented, which is easy to be performed. It is strictly proved that the numerical solution possesses the accuracy of second-order in time and higher order in space.