In this paper,a generalized Laguerre-spherical harmonic spectral method is proposed for the Cauchy problem of three-dimensional nonlinear Klein-Gordon equation. The goal is to make the numerical solutions to preserve ...In this paper,a generalized Laguerre-spherical harmonic spectral method is proposed for the Cauchy problem of three-dimensional nonlinear Klein-Gordon equation. The goal is to make the numerical solutions to preserve the same conservation as that for the exact solution.The stability and convergence of the proposed scheme are proved.Numerical results demonstrate the efficiency of this approach.We also establish some basic results on the generalized Laguerre-spherical harmonic orthogonal approximation,which play an important role in spectral methods for various problems defined on the whole space and unbounded domains with spherical geometry.展开更多
The stochastic simulation algorithm (SSA) accurately depicts spatially homogeneous wellstirred chemically reacting systems with small populations of chemical species and properly represents noise, but it is often ab...The stochastic simulation algorithm (SSA) accurately depicts spatially homogeneous wellstirred chemically reacting systems with small populations of chemical species and properly represents noise, but it is often abandoned when modeling larger systems because of its computational complexity. In this work, a twin support vector regression based stochastic simulations algorithm (TS^3A) is proposed by combining the twin support vector regression and SSA, the former is a well-known robust regression method in machine learning. Numerical results indicate that this proposed algorithm can be applied to a wide range of chemically reacting systems and obtain significant improvements on efficiency and accuracy with fewer simulating runs over the existing methods.展开更多
This paper proposes a new collocation method for initial value problems of second order ODEs based on the Laguerre-Gauss interpolation. It provides the global numerical solutions and possesses the spectral accuracy. N...This paper proposes a new collocation method for initial value problems of second order ODEs based on the Laguerre-Gauss interpolation. It provides the global numerical solutions and possesses the spectral accuracy. Numerical results demonstrate its high efficiency.展开更多
In this paper, an algorithm is proposed to solve the 0(2) symmetric positive solutions to the boundary value problem of the p-Henon equation. Taking 1 in the p- Henon equation as a bifurcation parameter, the symmetr...In this paper, an algorithm is proposed to solve the 0(2) symmetric positive solutions to the boundary value problem of the p-Henon equation. Taking 1 in the p- Henon equation as a bifurcation parameter, the symmetry-breaking bifurcation point on the branch of the O(2) symmetric positive solutions is found via the extended systems. The other symmetric positive solutions are computed by the branch switching method based on the Liapunov-Schmidt reduction.展开更多
Consider the design problem for estimation and extrapolation in approximately linear regression models with possible misspecification. The design space is a discrete set consisting of finitely many points, and the mod...Consider the design problem for estimation and extrapolation in approximately linear regression models with possible misspecification. The design space is a discrete set consisting of finitely many points, and the model bias comes from a reproducing kernel Hilbert space. Two different design criteria are proposed by applying the minimax approach for estimating the parameters of the regression response and extrapolating the regression response to points outside of the design space. A simulated annealing algorithm is applied to construct the minimax designs. These minimax designs are compared with the classical D-optimal designs and all-bias extrapolation designs. Numerical results indicate that the simulated annealing algorithm is feasible and the minimax designs are robust against bias caused by model misspecification.展开更多
The support vector machine (SVM) is a novel machine learning tool in data mining. In this paper, the geometric approach based on the compressed convex hull (CCH) with a mathematical framework is introduced to solv...The support vector machine (SVM) is a novel machine learning tool in data mining. In this paper, the geometric approach based on the compressed convex hull (CCH) with a mathematical framework is introduced to solve SVM classification problems. Compared with the reduced convex hull (RCH), CCH preserves the shape of geometric solids for data sets; meanwhile, it is easy to give the necessary and sufficient condition for determining its extreme points. As practical applications of CCH, spare and probabilistic speed-up geometric algorithms are developed. Results of numerical experiments show that the proposed algorithms can reduce kernel calculations and display nice performances.展开更多
In this paper, we investigate the mixed spectral method using generalized Laguerre functions for exterior problems of fourth order partial differential equations. A mixed spectral scheme is provided for the stream fun...In this paper, we investigate the mixed spectral method using generalized Laguerre functions for exterior problems of fourth order partial differential equations. A mixed spectral scheme is provided for the stream function form of the Navier-Stokes equations outside a disc. Numerical results demonstrate the spectral accuracy in space.展开更多
We prove existence and uniqueness of the global solution to the Cauchy problem on a universe fireworks model with finite total mass at the initial state when the ratio of the mass surviving the explosion, the probabil...We prove existence and uniqueness of the global solution to the Cauchy problem on a universe fireworks model with finite total mass at the initial state when the ratio of the mass surviving the explosion, the probability of the explosion of fragments and the probability function of the velocity change of a surviving particle satisfy the corresponding physical conditions. Although the nonrelativistic Boltzmann-like equation modeling the universe fireworks is mathematically easy, this article leads rather theoretically to an understanding of how to construct contractive mappings in a Banach space for the proof of the existence and uniqueness of the solution by means of methods taken from the famous work by DiPerna & Lions about the Boltzmann equation. We also show both the regularity and the time-asymptotic behavior of solution to the Cauchy problem.展开更多
A fully discrete Jacobi-spherical harmonic spectral method is provided for the Navier-Stokes equations in a ball. Its stability and convergence are proved. Numerical results show efficiency of this approach. The propo...A fully discrete Jacobi-spherical harmonic spectral method is provided for the Navier-Stokes equations in a ball. Its stability and convergence are proved. Numerical results show efficiency of this approach. The proposed method is also applicable to other problems in spherical geometry.展开更多
In this paper, we review some results on the spectral methods. We first consider the Jacobi spectral method and the generalized Jacobi spectral method for various problems, including degenerated and singular different...In this paper, we review some results on the spectral methods. We first consider the Jacobi spectral method and the generalized Jacobi spectral method for various problems, including degenerated and singular differential equations. Then we present the generalized Jacobi quasi-orthogonal approximation and its applica- tions to the spectral element methods for high order problems with mixed inhomogeneous boundary conditions. We also discuss the related spectral methods for non-rectangular domains and the irrational spectral methods for unbounded domains. Next, we consider the Hermite spectral method and the generalized Hermite spec- tral method with their applications. Finally, we consider the Laguerre spectral method and the generalized Laguerre spectral method for many problems defined on unbounded domains. We also present the generalized Laguerre quasi-orthogonal approximation and its applications to certain problems of non-standard type and exterior problems.展开更多
The authors investigate Petrov-Galerkin spectral element method. Some results on Legendre irrational quasi-orthogonal approximations are established, which play important roles in Petrov-Galerkin spectral element meth...The authors investigate Petrov-Galerkin spectral element method. Some results on Legendre irrational quasi-orthogonal approximations are established, which play important roles in Petrov-Galerkin spectral element method for mixed inhomogeneous boundary value problems of partial differential equations defined on polygons. As examples of applications, spectral element methods for two model problems, with the spectral accuracy in certain Jacobi weighted Sobolev spaces, are proposed. The techniques developed in this paper are also applicable to other higher order methods.展开更多
This paper discusses convergence and complexity of arbitrary,but fixed,order adaptive mixed element methods for the Poisson equation in two and three dimensions.The two main ingredients in the analysis,namely the quas...This paper discusses convergence and complexity of arbitrary,but fixed,order adaptive mixed element methods for the Poisson equation in two and three dimensions.The two main ingredients in the analysis,namely the quasi-orthogonality and the discrete reliability,are achieved by use of a discrete Helmholtz decomposition and a discrete inf-sup condition.The adaptive algorithms are shown to be contractive for the sum of the error of flux in L2-norm and the scaled error estimator after each step of mesh refinement and to be quasi-optimal with respect to the number of elements of underlying partitions.The methods do not require a separate treatment for the data oscillation.展开更多
This paper is contributed to the combinatorial properties of the MSS sequences, which are the periodic kneading words of quadratic maps defined on a interval. An explicit expression of adjacency relations on MSS seque...This paper is contributed to the combinatorial properties of the MSS sequences, which are the periodic kneading words of quadratic maps defined on a interval. An explicit expression of adjacency relations on MSS sequences of given lengths is established.展开更多
This paper deals with the delay-dependent stability of numerical methods for delay differential equations. First, a stability criterion of Runge-Kutta methods is extended to the case of general linear methods. Then, l...This paper deals with the delay-dependent stability of numerical methods for delay differential equations. First, a stability criterion of Runge-Kutta methods is extended to the case of general linear methods. Then, linear multistep methods are considered and a class of r(0)-stable methods are found. Later, some examples of r(0)-stable multistep multistage methods are given. Finally, numerical experiments are presented to confirm the theoretical results.展开更多
This paper is contributed to the combinatorial properties of the periodic kneading words of antisymmetric cubic maps defined on a interval. The least words of given lengths, the adjacency relations on the words of giv...This paper is contributed to the combinatorial properties of the periodic kneading words of antisymmetric cubic maps defined on a interval. The least words of given lengths, the adjacency relations on the words of given lengths and the parity-alternative property in some sets of such words are established.展开更多
In this paper, we introduce an efficient Chebyshev-Gauss spectral collocation method for initial value problems of ordinary differential equations. We first propose a single interval method and analyze its convergence...In this paper, we introduce an efficient Chebyshev-Gauss spectral collocation method for initial value problems of ordinary differential equations. We first propose a single interval method and analyze its convergence. We then develop a multi-interval method. The suggested algorithms enjoy spectral accuracy and can be implemented in stable and efficient manners. Some numerical comparisons with some popular methods are given to demonstrate the effectiveness of this approach.展开更多
We propose a collocation method for solving initial value problems of secondorder ODEs by using modified Laguerre functions.This new process provides global numerical solutions.Numerical results demonstrate the effici...We propose a collocation method for solving initial value problems of secondorder ODEs by using modified Laguerre functions.This new process provides global numerical solutions.Numerical results demonstrate the efficiency of the proposed algorithm.展开更多
To assess the lifetime distribution of highly reliable or expensive product,one of the most commonly used strategies is to construct step-stress accelerated degradation test(SSADT)which can curtail the test duration a...To assess the lifetime distribution of highly reliable or expensive product,one of the most commonly used strategies is to construct step-stress accelerated degradation test(SSADT)which can curtail the test duration and reduce the test cost.In reality,it is not unusual for a unit with a higher degradation rate which exhibits a more volatile degradation path.Recently,Ye,Chen,and Shen[(2015).A new class of Wiener process models for degradation analysis.Reliability Engineering and System Safety,139,58–67]proposed a Wiener process to capture the positive correlation between the drift rate and the volatility.In this paper,an optimal SSADT plan is developed under the assumption that the underlying degradation path follows the Wiener process with correlation.Firstly,the stochastic diffusion process is introduced to model a typical SSADT problem.Then the design variables,including the sample size,the measurement frequency and the numbers of measurements under each stress level,are optimised by minimising the asymptotic variance of the estimated p-percentile of the product’s lifetime distribution subject to the total experimental cost not exceeding a pre-specified budget.Finally,a numerical example is presented to illustrate the proposed method.展开更多
In this paper, we investigate spectral method for mixed inhomogeneous boundary value problems in three dimensions. Some results on the three-dimensional Legendre approxima- tion in Jacobi weighted Sobolev space are es...In this paper, we investigate spectral method for mixed inhomogeneous boundary value problems in three dimensions. Some results on the three-dimensional Legendre approxima- tion in Jacobi weighted Sobolev space are established, which improve and generalize the existing results, and play an important role in numerical solutions of partial differential equations. We also develop a lifting technique, with which we could handle mixed inho- mogeneous boundary conditions easily. As examples of applications, spectral schemes are provided for three model problems with mixed inhomogeneous boundary conditions. The spectral accuracy in space of proposed algorithms is proved. Efficient implementations are presented. Numerical results demonstrate their high accuracy, and confirm the theoretical analysis well.展开更多
In this paper, we investigate numerical methods for high order differential equations. We propose new spectral and spectral element methods for high order problems with mixed inhomogeneous boundary conditions, and pro...In this paper, we investigate numerical methods for high order differential equations. We propose new spectral and spectral element methods for high order problems with mixed inhomogeneous boundary conditions, and prove their spectral accuracy by using the recent results on the Jacobi quasi-orthogonal approximation. Numerical results demonstrate the high accuracy of suggested algorithm, which also works well even for oscillating solutions.展开更多
基金supported in part by NSF of China N.10871131The Science and Technology Commission of Shanghai Municipality,Grant N.075105118+1 种基金Shanghai Leading Academic Discipline Project N.T0401Fund for E-institute of Shanghai Universities N.E03004.
文摘In this paper,a generalized Laguerre-spherical harmonic spectral method is proposed for the Cauchy problem of three-dimensional nonlinear Klein-Gordon equation. The goal is to make the numerical solutions to preserve the same conservation as that for the exact solution.The stability and convergence of the proposed scheme are proved.Numerical results demonstrate the efficiency of this approach.We also establish some basic results on the generalized Laguerre-spherical harmonic orthogonal approximation,which play an important role in spectral methods for various problems defined on the whole space and unbounded domains with spherical geometry.
基金This work was supported by the National Natural Science Foundation of China (No.30871341), the National High-Tech Research and Development Program of China (No.2006AA02-Z190), the Shanghai Leading Academic Discipline Project (No.S30405), and the Natural Science Foundation of Shanghai Normal University (No.SK200937).
文摘The stochastic simulation algorithm (SSA) accurately depicts spatially homogeneous wellstirred chemically reacting systems with small populations of chemical species and properly represents noise, but it is often abandoned when modeling larger systems because of its computational complexity. In this work, a twin support vector regression based stochastic simulations algorithm (TS^3A) is proposed by combining the twin support vector regression and SSA, the former is a well-known robust regression method in machine learning. Numerical results indicate that this proposed algorithm can be applied to a wide range of chemically reacting systems and obtain significant improvements on efficiency and accuracy with fewer simulating runs over the existing methods.
基金supported by the National Natural Science Foundation of China(No.11171227)the Ph.D.Programs Foundation of Ministry of Education of China(No.20080270001)+2 种基金the Shanghai Leading Academic Discipline Project(No.S30405)the Fund for E-Institute of Shanghai Universities(No.E03004)the Foundation for Distinguished Young Talents in Higher Education of Guangdong of China(No.LYM09138)
文摘This paper proposes a new collocation method for initial value problems of second order ODEs based on the Laguerre-Gauss interpolation. It provides the global numerical solutions and possesses the spectral accuracy. Numerical results demonstrate its high efficiency.
基金Project supported by the National Natural Science Foundation of China (No. 10901106)the Shanghai Leading Academic Discipline Project (No. S30405)+2 种基金the Shanghai Normal University Academic Project (No. SK200936)the Natural Science Foundation of Shanghai (No. 09ZR1423200)the Innovation Program of Shanghai Municipal Education Commission (No. 09YZ150)
文摘In this paper, an algorithm is proposed to solve the 0(2) symmetric positive solutions to the boundary value problem of the p-Henon equation. Taking 1 in the p- Henon equation as a bifurcation parameter, the symmetry-breaking bifurcation point on the branch of the O(2) symmetric positive solutions is found via the extended systems. The other symmetric positive solutions are computed by the branch switching method based on the Liapunov-Schmidt reduction.
基金Supported by National Natural Science Foundation of China(11471216,11301332)E-Institutes of Shanghai Municipal Education Commission(E03004)+1 种基金Central Finance Project(YC-XK-13105)Shanghai Municipal Science and Technology Research Project(14DZ1201902)
文摘Consider the design problem for estimation and extrapolation in approximately linear regression models with possible misspecification. The design space is a discrete set consisting of finitely many points, and the model bias comes from a reproducing kernel Hilbert space. Two different design criteria are proposed by applying the minimax approach for estimating the parameters of the regression response and extrapolating the regression response to points outside of the design space. A simulated annealing algorithm is applied to construct the minimax designs. These minimax designs are compared with the classical D-optimal designs and all-bias extrapolation designs. Numerical results indicate that the simulated annealing algorithm is feasible and the minimax designs are robust against bias caused by model misspecification.
基金Supported by the National Natural Science Foundation of China (No. 30571059)the National High-Tech Research and Development Program of China (No. 2006AA02Z190)Shanghai Leading Academic Discipline Project (No. 530405)
文摘The support vector machine (SVM) is a novel machine learning tool in data mining. In this paper, the geometric approach based on the compressed convex hull (CCH) with a mathematical framework is introduced to solve SVM classification problems. Compared with the reduced convex hull (RCH), CCH preserves the shape of geometric solids for data sets; meanwhile, it is easy to give the necessary and sufficient condition for determining its extreme points. As practical applications of CCH, spare and probabilistic speed-up geometric algorithms are developed. Results of numerical experiments show that the proposed algorithms can reduce kernel calculations and display nice performances.
基金supported by the National Natural Science Foundation of China (No.10871131)the Science and Technology Commission of Shanghai Municipality (No.075105118)+1 种基金the Shanghai Leading Academic Discipline Project (No.S30405)the Fund for E-institutes of Shanghai Universities(No.E03004)
文摘In this paper, we investigate the mixed spectral method using generalized Laguerre functions for exterior problems of fourth order partial differential equations. A mixed spectral scheme is provided for the stream function form of the Navier-Stokes equations outside a disc. Numerical results demonstrate the spectral accuracy in space.
基金supported by NSFC(10271121)the Scientific Research Foundation for the Returned Overseas Chinese Scholars,the Education Ministry of China+5 种基金sponsored by joint grants of NSFC 10511120278/10611120371RFBR 04-02-39026supported by NSFC (10671130)E-Institutes of Shanghai Municipal Education Commission (E03004)Shanghai Science and Technology Commission (06JC14092)Shuguang Project of Shanghai Municipal Education Commission (06SG45)
文摘We prove existence and uniqueness of the global solution to the Cauchy problem on a universe fireworks model with finite total mass at the initial state when the ratio of the mass surviving the explosion, the probability of the explosion of fragments and the probability function of the velocity change of a surviving particle satisfy the corresponding physical conditions. Although the nonrelativistic Boltzmann-like equation modeling the universe fireworks is mathematically easy, this article leads rather theoretically to an understanding of how to construct contractive mappings in a Banach space for the proof of the existence and uniqueness of the solution by means of methods taken from the famous work by DiPerna & Lions about the Boltzmann equation. We also show both the regularity and the time-asymptotic behavior of solution to the Cauchy problem.
基金Project supported by the National Natural Science Foundation of China(No.10771142)Science and Technology Commission of Shanghai Municipality(No.75105118)+2 种基金Shanghai Leading Academic Discipline Projects(Nos.T0401 and J50101)Fund for E-institutes of Universities in Shanghai(No.E03004)and Innovative Foundation of Shanghai University(No.A.10-0101-07-408)
文摘A fully discrete Jacobi-spherical harmonic spectral method is provided for the Navier-Stokes equations in a ball. Its stability and convergence are proved. Numerical results show efficiency of this approach. The proposed method is also applicable to other problems in spherical geometry.
基金supported by National Natural Science Foundation of China(Grant No.11171227)Fund for Doctoral Authority of China(Grant No.20123127110001)+1 种基金Fund for E-institute of Shanghai Universities(Grant No.E03004)Leading Academic Discipline Project of Shanghai Municipal Education Commission(Grant No.J50101)
文摘In this paper, we review some results on the spectral methods. We first consider the Jacobi spectral method and the generalized Jacobi spectral method for various problems, including degenerated and singular differential equations. Then we present the generalized Jacobi quasi-orthogonal approximation and its applica- tions to the spectral element methods for high order problems with mixed inhomogeneous boundary conditions. We also discuss the related spectral methods for non-rectangular domains and the irrational spectral methods for unbounded domains. Next, we consider the Hermite spectral method and the generalized Hermite spec- tral method with their applications. Finally, we consider the Laguerre spectral method and the generalized Laguerre spectral method for many problems defined on unbounded domains. We also present the generalized Laguerre quasi-orthogonal approximation and its applications to certain problems of non-standard type and exterior problems.
基金supported by the National Natural Science Foundation of China (No. 10871131)the Fund for Doctoral Authority of China (No. 200802700001)+1 种基金the Shanghai Leading Academic Discipline Project(No. S30405)the Fund for E-institutes of Shanghai Universities (No. E03004)
文摘The authors investigate Petrov-Galerkin spectral element method. Some results on Legendre irrational quasi-orthogonal approximations are established, which play important roles in Petrov-Galerkin spectral element method for mixed inhomogeneous boundary value problems of partial differential equations defined on polygons. As examples of applications, spectral element methods for two model problems, with the spectral accuracy in certain Jacobi weighted Sobolev spaces, are proposed. The techniques developed in this paper are also applicable to other higher order methods.
基金supported by National Natural Science Foundation of China (Grant Nos. 11171219,11161130004)E-Institutes of Shanghai Municipal Education Commission(Grant No. E03004)+1 种基金supported by Shanghai Leading Discipline Project(Grant No. N.S30405)Shanghai Normal University Research Program (Grant No. SK201202)
文摘This paper discusses convergence and complexity of arbitrary,but fixed,order adaptive mixed element methods for the Poisson equation in two and three dimensions.The two main ingredients in the analysis,namely the quasi-orthogonality and the discrete reliability,are achieved by use of a discrete Helmholtz decomposition and a discrete inf-sup condition.The adaptive algorithms are shown to be contractive for the sum of the error of flux in L2-norm and the scaled error estimator after each step of mesh refinement and to be quasi-optimal with respect to the number of elements of underlying partitions.The methods do not require a separate treatment for the data oscillation.
基金the National Natural Science Foundation of China (Grant No.10731040)
文摘This paper is contributed to the combinatorial properties of the MSS sequences, which are the periodic kneading words of quadratic maps defined on a interval. An explicit expression of adjacency relations on MSS sequences of given lengths is established.
基金supported by National Natural Science Foundation of China(Grant No.10671078)the Program for New Century Excellent Talents in University,the State Education Ministry of China. supported in part by E-Institutes of Shanghai Municipal Education Commission (No.E03004)+3 种基金National Natural Science Foundation of China(No.10671130)Shanghai Science and Technology Commission(No.06JC14092)Shuguang Project of Shanghai Municipal Education Commission(No.06SG45)the Shanghai Leading Academic Discipline Project(No.S30405)
文摘This paper deals with the delay-dependent stability of numerical methods for delay differential equations. First, a stability criterion of Runge-Kutta methods is extended to the case of general linear methods. Then, linear multistep methods are considered and a class of r(0)-stable methods are found. Later, some examples of r(0)-stable multistep multistage methods are given. Finally, numerical experiments are presented to confirm the theoretical results.
基金the National Natural Science Foundation of China (Grant No.10731040)
文摘This paper is contributed to the combinatorial properties of the periodic kneading words of antisymmetric cubic maps defined on a interval. The least words of given lengths, the adjacency relations on the words of given lengths and the parity-alternative property in some sets of such words are established.
文摘In this paper, we introduce an efficient Chebyshev-Gauss spectral collocation method for initial value problems of ordinary differential equations. We first propose a single interval method and analyze its convergence. We then develop a multi-interval method. The suggested algorithms enjoy spectral accuracy and can be implemented in stable and efficient manners. Some numerical comparisons with some popular methods are given to demonstrate the effectiveness of this approach.
基金supported in part by Foundation for Distinguished Young Talents in Higher Education of Guangdong,China,N.LYM09138supported in part by NSF of China N.10871131+2 种基金Fund for Doctor Authority of Chinese Educational Ministry N.20080270001Shanghai Leading Academic Discipline Project N.S30405Fund for E-institutes of Shanghai Universities N.E03004.
文摘We propose a collocation method for solving initial value problems of secondorder ODEs by using modified Laguerre functions.This new process provides global numerical solutions.Numerical results demonstrate the efficiency of the proposed algorithm.
基金Rong-Xian Yue’s research was supported by the National Natural Science Foundation of China[grant number 11471216]Daojiang He’s research was supported by the National Natural Science Foundation of China[grant number 11201005].
文摘To assess the lifetime distribution of highly reliable or expensive product,one of the most commonly used strategies is to construct step-stress accelerated degradation test(SSADT)which can curtail the test duration and reduce the test cost.In reality,it is not unusual for a unit with a higher degradation rate which exhibits a more volatile degradation path.Recently,Ye,Chen,and Shen[(2015).A new class of Wiener process models for degradation analysis.Reliability Engineering and System Safety,139,58–67]proposed a Wiener process to capture the positive correlation between the drift rate and the volatility.In this paper,an optimal SSADT plan is developed under the assumption that the underlying degradation path follows the Wiener process with correlation.Firstly,the stochastic diffusion process is introduced to model a typical SSADT problem.Then the design variables,including the sample size,the measurement frequency and the numbers of measurements under each stress level,are optimised by minimising the asymptotic variance of the estimated p-percentile of the product’s lifetime distribution subject to the total experimental cost not exceeding a pre-specified budget.Finally,a numerical example is presented to illustrate the proposed method.
文摘In this paper, we investigate spectral method for mixed inhomogeneous boundary value problems in three dimensions. Some results on the three-dimensional Legendre approxima- tion in Jacobi weighted Sobolev space are established, which improve and generalize the existing results, and play an important role in numerical solutions of partial differential equations. We also develop a lifting technique, with which we could handle mixed inho- mogeneous boundary conditions easily. As examples of applications, spectral schemes are provided for three model problems with mixed inhomogeneous boundary conditions. The spectral accuracy in space of proposed algorithms is proved. Efficient implementations are presented. Numerical results demonstrate their high accuracy, and confirm the theoretical analysis well.
文摘In this paper, we investigate numerical methods for high order differential equations. We propose new spectral and spectral element methods for high order problems with mixed inhomogeneous boundary conditions, and prove their spectral accuracy by using the recent results on the Jacobi quasi-orthogonal approximation. Numerical results demonstrate the high accuracy of suggested algorithm, which also works well even for oscillating solutions.