Magnetic reconnection and tearing mode instability play a critical role in many physical processes.The application of Galerkin spectral method for tearing mode instability in two-dimensional geometry is investigated i...Magnetic reconnection and tearing mode instability play a critical role in many physical processes.The application of Galerkin spectral method for tearing mode instability in two-dimensional geometry is investigated in this paper.A resistive magnetohydrodynamic code is developed,by the Galerkin spectral method both in the periodic and aperiodic directions.Spectral schemes are provided for global modes and local modes.Mode structures,resistivity scaling,convergence and stability of tearing modes are discussed.The effectiveness of the code is demonstrated,and the computational results are compared with the results using Galerkin spectral method only in the periodic direction.The numerical results show that the code using Galerkin spectral method individually allows larger time step in global and local modes simulations,and has better convergence in global modes simulations.展开更多
In this paper, we design and analyze a space-time spectral method for the subdiffusion equation.Here, we are facing two difficulties. The first is that the solutions of this equation are usually singular near the init...In this paper, we design and analyze a space-time spectral method for the subdiffusion equation.Here, we are facing two difficulties. The first is that the solutions of this equation are usually singular near the initial time. Consequently, traditional high-order numerical methods in time are inefficient. The second obstacle is that the resulting system of the space-time spectral approach is usually large and time-consuming to solve. We aim at overcoming the first difficulty by proposing a novel approach in time, which is based on variable transformation techniques. Suitable ψ-fractional Sobolev spaces and a new variational framework are introduced to establish the well-posedness of the associated variational problem. This allows us to construct our space-time spectral method using a combination of temporal generalized Jacobi polynomials(GJPs) and spatial Legendre polynomials. For the second difficulty, we propose a fast algorithm to effectively solve the resulting linear system. The fast algorithm makes use of a matrix diagonalization in space and QZ decomposition in time. Our analysis and numerical experiments show that the proposed method is exponentially convergent with respect to the polynomial degrees in both space and time directions, even though the exact solution has very limited regularity.展开更多
In this paper,an initial boundary value problem of the space-time fractional diffusion equation is studied.Both temporal and spatial directions for this equation are discreted by the Galerkin spectral methods.And then...In this paper,an initial boundary value problem of the space-time fractional diffusion equation is studied.Both temporal and spatial directions for this equation are discreted by the Galerkin spectral methods.And then based on the discretization scheme,reliable a posteriori error estimates for the spectral approximation are derived.Some numerical examples are presented to verify the validity and applicability of the derived a posteriori error estimator.展开更多
In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy o...In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.展开更多
In this paper, a-posteriori error estimators are proposed for the Legendre spectral Galerkin method for two-point boundary value problems. The key idea is to postprocess the Galerkin approximation, and the analysis sh...In this paper, a-posteriori error estimators are proposed for the Legendre spectral Galerkin method for two-point boundary value problems. The key idea is to postprocess the Galerkin approximation, and the analysis shows that the postproeess improves the order of convergence. Consequently, we obtain asymptotically exact aposteriori error estimators based on the postprocessing results. Numerical examples are included to illustrate the theoretical analysis.展开更多
This study presents a new tool for solving stochastic boundary-value problems. This tool is created by modify the previous spectral stochastic meshless local Petrov-Galerkin method using the MLPG5 scheme. This modifie...This study presents a new tool for solving stochastic boundary-value problems. This tool is created by modify the previous spectral stochastic meshless local Petrov-Galerkin method using the MLPG5 scheme. This modified spectral stochastic meshless local Petrov-Galerkin method is selectively applied to predict the structural failure probability with the uncertainty in the spatial variability of mechanical properties. Except for the MLPG5 scheme, deriving the proposed spectral stochastic meshless local Petrov-Galerkin formulation adopts generalized polynomial chaos expansions of random mechanical properties. Predicting the structural failure probability is based on the first-order reliability method. Further comparing the spectral stochastic finite element-based and meshless local Petrov-Galerkin-based predicted structural failure probabilities indicates that the proposed spectral stochastic meshless local Petrov-Galerkin method predicts the more accurate structural failure probability than the spectral stochastic finite element method does. In addition, generating spectral stochastic meshless local Petrov-Galerkin results are considerably time-saving than generating Monte-Carlo simulation results does. In conclusion, the spectral stochastic meshless local Petrov-Galerkin method serves as a time-saving tool for solving stochastic boundary-value problems sufficiently accurately.展开更多
For fractional Volterra integro-differential equations(FVIDEs)with weakly singular kernels,this paper proposes a generalized Jacobi spectral Galerkin method.The basis functions for the provided method are selected gen...For fractional Volterra integro-differential equations(FVIDEs)with weakly singular kernels,this paper proposes a generalized Jacobi spectral Galerkin method.The basis functions for the provided method are selected generalized Jacobi functions(GJFs),which can be utilized as natural basis functions of spectral methods for weakly singular FVIDEs when appropriately constructed.The developed method's spectral rate of convergence is determined using the L^(∞)-norm and the weighted L^(2)-norm.Numerical results indicate the usefulness of the proposed method.展开更多
This paper presents a fully spectral discretization method for solving KdV equations with periodic boundary conditions.Chebyshev pseudospectral approximation in the time direction and Fourier Galerkin approximation in...This paper presents a fully spectral discretization method for solving KdV equations with periodic boundary conditions.Chebyshev pseudospectral approximation in the time direction and Fourier Galerkin approximation in the spatial direction.The expansion coefficients are determined by minimizing an object funictional.Rapid convergence of the method is proved.展开更多
The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectra...The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectral volume schemes for solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element method for solving hyperbolic conservation laws. An orthogonal spectral basis written in terms of Jacobi polynomials is applied. High computational efficiency is obtained due to such matrix-free algorithm. The formulation is conservative, and essential nomoscillation is enforced by the HR limiter. We show that HR preserves the order of accuracy of the spectral/hp element method for smooth solution problems and generate essentially non-oscillatory solutions profiles for capturing discontinuous solutions without local characteristic decomposition. In addition, we introduce a postprocessing technique to improve HR for limiting high degree numerical solutions.展开更多
基金Project supported by the Sichuan Science and Technology Program(Grant No.22YYJC1286)the China National Magnetic Confinement Fusion Science Program(Grant No.2013GB112005)the National Natural Science Foundation of China(Grant Nos.12075048 and 11925501)。
文摘Magnetic reconnection and tearing mode instability play a critical role in many physical processes.The application of Galerkin spectral method for tearing mode instability in two-dimensional geometry is investigated in this paper.A resistive magnetohydrodynamic code is developed,by the Galerkin spectral method both in the periodic and aperiodic directions.Spectral schemes are provided for global modes and local modes.Mode structures,resistivity scaling,convergence and stability of tearing modes are discussed.The effectiveness of the code is demonstrated,and the computational results are compared with the results using Galerkin spectral method only in the periodic direction.The numerical results show that the code using Galerkin spectral method individually allows larger time step in global and local modes simulations,and has better convergence in global modes simulations.
基金supported by National Natural Science Foundation of China (Grant No. 11971408)。
文摘In this paper, we design and analyze a space-time spectral method for the subdiffusion equation.Here, we are facing two difficulties. The first is that the solutions of this equation are usually singular near the initial time. Consequently, traditional high-order numerical methods in time are inefficient. The second obstacle is that the resulting system of the space-time spectral approach is usually large and time-consuming to solve. We aim at overcoming the first difficulty by proposing a novel approach in time, which is based on variable transformation techniques. Suitable ψ-fractional Sobolev spaces and a new variational framework are introduced to establish the well-posedness of the associated variational problem. This allows us to construct our space-time spectral method using a combination of temporal generalized Jacobi polynomials(GJPs) and spatial Legendre polynomials. For the second difficulty, we propose a fast algorithm to effectively solve the resulting linear system. The fast algorithm makes use of a matrix diagonalization in space and QZ decomposition in time. Our analysis and numerical experiments show that the proposed method is exponentially convergent with respect to the polynomial degrees in both space and time directions, even though the exact solution has very limited regularity.
基金supported by the State Key Program of National Natural Science Foundation of China(No.11931003)National Natural Science Foundation of China(Nos.41974133,11671157 and 11971410)supported by the Innovation Project of Graduate School of South China Normal University(No.2018LKXM008).
文摘In this paper,an initial boundary value problem of the space-time fractional diffusion equation is studied.Both temporal and spatial directions for this equation are discreted by the Galerkin spectral methods.And then based on the discretization scheme,reliable a posteriori error estimates for the spectral approximation are derived.Some numerical examples are presented to verify the validity and applicability of the derived a posteriori error estimator.
文摘In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.
基金supported partially by the innovation fund of Shanghai Normal Universitysupported partially by NSERC of Canada under Grant OGP0046726.
文摘In this paper, a-posteriori error estimators are proposed for the Legendre spectral Galerkin method for two-point boundary value problems. The key idea is to postprocess the Galerkin approximation, and the analysis shows that the postproeess improves the order of convergence. Consequently, we obtain asymptotically exact aposteriori error estimators based on the postprocessing results. Numerical examples are included to illustrate the theoretical analysis.
文摘This study presents a new tool for solving stochastic boundary-value problems. This tool is created by modify the previous spectral stochastic meshless local Petrov-Galerkin method using the MLPG5 scheme. This modified spectral stochastic meshless local Petrov-Galerkin method is selectively applied to predict the structural failure probability with the uncertainty in the spatial variability of mechanical properties. Except for the MLPG5 scheme, deriving the proposed spectral stochastic meshless local Petrov-Galerkin formulation adopts generalized polynomial chaos expansions of random mechanical properties. Predicting the structural failure probability is based on the first-order reliability method. Further comparing the spectral stochastic finite element-based and meshless local Petrov-Galerkin-based predicted structural failure probabilities indicates that the proposed spectral stochastic meshless local Petrov-Galerkin method predicts the more accurate structural failure probability than the spectral stochastic finite element method does. In addition, generating spectral stochastic meshless local Petrov-Galerkin results are considerably time-saving than generating Monte-Carlo simulation results does. In conclusion, the spectral stochastic meshless local Petrov-Galerkin method serves as a time-saving tool for solving stochastic boundary-value problems sufficiently accurately.
基金supported by the State Key Program of National Natural Science Foundation of China(Grant No.11931003)by the National Natural Science Foundation of China(Grant Nos.41974133,12126325)by the Postgraduate Scientific Research Innovation Project of Hunan Province(Grant No.CX20200620).
文摘For fractional Volterra integro-differential equations(FVIDEs)with weakly singular kernels,this paper proposes a generalized Jacobi spectral Galerkin method.The basis functions for the provided method are selected generalized Jacobi functions(GJFs),which can be utilized as natural basis functions of spectral methods for weakly singular FVIDEs when appropriately constructed.The developed method's spectral rate of convergence is determined using the L^(∞)-norm and the weighted L^(2)-norm.Numerical results indicate the usefulness of the proposed method.
文摘This paper presents a fully spectral discretization method for solving KdV equations with periodic boundary conditions.Chebyshev pseudospectral approximation in the time direction and Fourier Galerkin approximation in the spatial direction.The expansion coefficients are determined by minimizing an object funictional.Rapid convergence of the method is proved.
基金Research was supported in part by NSF grant DMS-0800612Research was supported by Applied Mathematics program of the US DOE Office of Advanced Scientific Computing ResearchThe Pacific Northwest National Laboratory is operated by Battelle for the U.S. Department of Energy under Contract DE-AC05-76RL01830
文摘The hierarchical reconstruction （HR） [Liu, Shu, Tadmor and Zhang, SINUM ＇07] has been successfully applied to prevent oscillations in solutions computed by finite volume, Runge-Kutta discontinuous Galerkin, spectral volume schemes for solving hyperbolic conservation laws. In this paper, we demonstrate that HR can also be combined with spectral/hp element method for solving hyperbolic conservation laws. An orthogonal spectral basis written in terms of Jacobi polynomials is applied. High computational efficiency is obtained due to such matrix-free algorithm. The formulation is conservative, and essential nomoscillation is enforced by the HR limiter. We show that HR preserves the order of accuracy of the spectral/hp element method for smooth solution problems and generate essentially non-oscillatory solutions profiles for capturing discontinuous solutions without local characteristic decomposition. In addition, we introduce a postprocessing technique to improve HR for limiting high degree numerical solutions.
