A non-orthogonal predefined exponential basis set is used to handle half-bounded domains in multi domain spectral method (MDSM). This approach works extremely well for real-valued semi-infinite differential problems. ...A non-orthogonal predefined exponential basis set is used to handle half-bounded domains in multi domain spectral method (MDSM). This approach works extremely well for real-valued semi-infinite differential problems. It spans simultaneously wide range of exponential decay rates with multi scaling and does not suffer from zero crossing. These two conditions are necessary for many physical problems. For comparison, the method is used to solve different problems and compared with analytical and published results. The comparison exhibits the strengths and accuracy of the presented basis set.展开更多
An efficient direct spectral domain decomposition method is developed coupled with Chebyshev spectral approximation for the solution of 2D, unsteady and in- compressible Navier-Stokes equations in complex geometries. ...An efficient direct spectral domain decomposition method is developed coupled with Chebyshev spectral approximation for the solution of 2D, unsteady and in- compressible Navier-Stokes equations in complex geometries. In this numerical approach, the spatial domains of interest are decomposed into several non-overlapping rectangu- lar sub-domains. In each sub-domain, an improved projection scheme with second-order accuracy is used to deal with the coupling of velocity and pressure, and the Chebyshev collocation spectral method （CSM） is adopted to execute the spatial discretization. The influence matrix technique is employed to enforce the continuities of both variables and their normal derivatives between the adjacent sub-domains. The imposing of the Neu- mann boundary conditions to the Poisson equations of pressure and intermediate variable will result in the indeterminate solution. A new strategy of assuming the Dirichlet bound- ary conditions on interface and using the first-order normal derivatives as transmission conditions to keep the continuities of variables is proposed to overcome this trouble. Three test cases are used to verify the accuracy and efficiency, and the detailed comparison be- tween the numerical results and the available solutions is done. The results indicate that the present method is efficiency, stability, and accuracy.展开更多
A Legendre spectral element/Laguerre coupled method is proposed to numerically solve the elliptic Helmholtz problem on the half line. Rigorous analysis is carried out to establish the convergence of the method. Severa...A Legendre spectral element/Laguerre coupled method is proposed to numerically solve the elliptic Helmholtz problem on the half line. Rigorous analysis is carried out to establish the convergence of the method. Several numerical examples are provided to confirm the theoretical results. The advantage of this method is demonstrated by a numerical comparison with the pure Laguerre method.展开更多
In this paper we propose a spectral element: vanishing viscosity (SEW) method for the conservation laws on the semi-infinite interval. By using a suitable mapping, the problem is first transformed into a modified cons...In this paper we propose a spectral element: vanishing viscosity (SEW) method for the conservation laws on the semi-infinite interval. By using a suitable mapping, the problem is first transformed into a modified conservation law in a bounded interval, then the well-known spectral vanishing viscosity technique is generalized to the multi-domain case in order to approximate this trarsformed equation more efficiently. The construction details and convergence analysis are presented. Under a usual assumption of boundedness of the approximation solutions, it is proven that the solution of the SEW approximation converges to the uniciue entropy solution of the conservation laws. A number of numerical tests is carried out to confirm the theoretical results.展开更多
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.展开更多
This paper presents the extension of the global description approach of a discontinuous function, which is proposed in the previous paper, to a spectral domain decomposition method. This multi-domain spectral immersed...This paper presents the extension of the global description approach of a discontinuous function, which is proposed in the previous paper, to a spectral domain decomposition method. This multi-domain spectral immersed interlace method（IIM） divides the whole computation domain into the smooth and discontinuous parts. Fewer points on the smooth domains are used via taking advantage of the high accuracy property of the spectral method, but more points on the discontinuous domains are employed to enhance the resolution of the calculation. Two that the domain decomposition technique can placed around the discontinuity. The present reached, in spite of the enlarged computational discontinuous problems are tested to verify the present method. The results show reduce the error of the spectral IIM, especially when more collocation points are method is t：avorable for the reason that the same level of the accuracy can be domain.展开更多
This paper is concerned with a multi-domain spectral method, based on an interior penalty discontinuous Galerkin （IPDG） formulation, for the exterior Helmholtz problem truncated via an exact circular or spherical Di...This paper is concerned with a multi-domain spectral method, based on an interior penalty discontinuous Galerkin （IPDG） formulation, for the exterior Helmholtz problem truncated via an exact circular or spherical Dirichlet-to-Neumann （DtN） boundary con- dition. An effective iterative approach is proposed to localize the global DtN boundary condition, which facilitates the implementation of multi-domain methods, and the treat- ment for complex geometry of the scatterers. Under a discontinuous Galerkin formulation, the proposed method allows to use polynomial basis functions of different degree on dif- ferent subdomains, and more importantly, explicit wave number dependence estimates of the spectral scheme can be derived, which is somehow implausible for a multi-domain continuous Galerkin formulation.展开更多
文摘A non-orthogonal predefined exponential basis set is used to handle half-bounded domains in multi domain spectral method (MDSM). This approach works extremely well for real-valued semi-infinite differential problems. It spans simultaneously wide range of exponential decay rates with multi scaling and does not suffer from zero crossing. These two conditions are necessary for many physical problems. For comparison, the method is used to solve different problems and compared with analytical and published results. The comparison exhibits the strengths and accuracy of the presented basis set.
基金Project supported by the National Natural Science Foundation of China(No.51176026)the Fundamental Research Funds for the Central Universities(No.DUT14RC(3)029)
文摘An efficient direct spectral domain decomposition method is developed coupled with Chebyshev spectral approximation for the solution of 2D, unsteady and in- compressible Navier-Stokes equations in complex geometries. In this numerical approach, the spatial domains of interest are decomposed into several non-overlapping rectangu- lar sub-domains. In each sub-domain, an improved projection scheme with second-order accuracy is used to deal with the coupling of velocity and pressure, and the Chebyshev collocation spectral method （CSM） is adopted to execute the spatial discretization. The influence matrix technique is employed to enforce the continuities of both variables and their normal derivatives between the adjacent sub-domains. The imposing of the Neu- mann boundary conditions to the Poisson equations of pressure and intermediate variable will result in the indeterminate solution. A new strategy of assuming the Dirichlet bound- ary conditions on interface and using the first-order normal derivatives as transmission conditions to keep the continuities of variables is proposed to overcome this trouble. Three test cases are used to verify the accuracy and efficiency, and the detailed comparison be- tween the numerical results and the available solutions is done. The results indicate that the present method is efficiency, stability, and accuracy.
基金This work was supported by Natural Science Foundation of Fujian under Grant A0310002the Excellent Young Teachers Program (EYTP) of the Ministry of Education of China.
文摘A Legendre spectral element/Laguerre coupled method is proposed to numerically solve the elliptic Helmholtz problem on the half line. Rigorous analysis is carried out to establish the convergence of the method. Several numerical examples are provided to confirm the theoretical results. The advantage of this method is demonstrated by a numerical comparison with the pure Laguerre method.
文摘In this paper we propose a spectral element: vanishing viscosity (SEW) method for the conservation laws on the semi-infinite interval. By using a suitable mapping, the problem is first transformed into a modified conservation law in a bounded interval, then the well-known spectral vanishing viscosity technique is generalized to the multi-domain case in order to approximate this trarsformed equation more efficiently. The construction details and convergence analysis are presented. Under a usual assumption of boundedness of the approximation solutions, it is proven that the solution of the SEW approximation converges to the uniciue entropy solution of the conservation laws. A number of numerical tests is carried out to confirm the theoretical results.
文摘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.
基金National Natural Science Foundation of China(51076006)
文摘This paper presents the extension of the global description approach of a discontinuous function, which is proposed in the previous paper, to a spectral domain decomposition method. This multi-domain spectral immersed interlace method（IIM） divides the whole computation domain into the smooth and discontinuous parts. Fewer points on the smooth domains are used via taking advantage of the high accuracy property of the spectral method, but more points on the discontinuous domains are employed to enhance the resolution of the calculation. Two that the domain decomposition technique can placed around the discontinuity. The present reached, in spite of the enlarged computational discontinuous problems are tested to verify the present method. The results show reduce the error of the spectral IIM, especially when more collocation points are method is t：avorable for the reason that the same level of the accuracy can be domain.
基金The work of the first author was partially supported by the National Natural Science Foundation of China （11026065 and 11101196）. This work was largely done when this author worked as a Research Fellow in Nanyang Technological University. The work of the second author was supported by the National Natural Science Foundation of China （11201166）. The work of the third author was supported by Singapore MOE Tier 1 Grant （2013-2016）, and Singapore A＊STAR-SERC-PSF Grant 122-PSF-0007.
文摘This paper is concerned with a multi-domain spectral method, based on an interior penalty discontinuous Galerkin （IPDG） formulation, for the exterior Helmholtz problem truncated via an exact circular or spherical Dirichlet-to-Neumann （DtN） boundary con- dition. An effective iterative approach is proposed to localize the global DtN boundary condition, which facilitates the implementation of multi-domain methods, and the treat- ment for complex geometry of the scatterers. Under a discontinuous Galerkin formulation, the proposed method allows to use polynomial basis functions of different degree on dif- ferent subdomains, and more importantly, explicit wave number dependence estimates of the spectral scheme can be derived, which is somehow implausible for a multi-domain continuous Galerkin formulation.
