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 propose a new two-level preconditioned C-G method which uses the quadratic smoothing and the linear correction in distorted but topologically structured grid.The CPU time of this method is less than t...In this paper,we propose a new two-level preconditioned C-G method which uses the quadratic smoothing and the linear correction in distorted but topologically structured grid.The CPU time of this method is less than that of the multigrid preconditioned C-G method(MGCG)using the quadratic element,but their accuracy is almost the same.Numerical experiments and eigenvalue analysis are given and the results show that the proposed two-level preconditioned method is efficient.展开更多
In this article,we discuss a least-squares/fictitious domain method for the solution of linear elliptic boundary value problems with Robin boundary conditions.LetΩandωbe two bounded domains of R d such thatω⊂Ω.For a...In this article,we discuss a least-squares/fictitious domain method for the solution of linear elliptic boundary value problems with Robin boundary conditions.LetΩandωbe two bounded domains of R d such thatω⊂Ω.For a linear elliptic problem inΩ\ωwith Robin boundary condition on the boundaryγofω,our goal here is to develop a fictitious domain method where one solves a variant of the original problem on the fullΩ,followed by a well-chosen correction overω.This method is of the virtual control type and relies on a least-squares formulation making the problem solvable by a conjugate gradient algorithm operating in a well chosen control space.Numerical results obtained when applying our method to the solution of two-dimensional elliptic and parabolic problems are given;they suggest optimal order of convergence.展开更多
基金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.
基金The work is partially supported by Youth Foundation of Sichuan University No 2010SCU11072Doctoral Fund of Ministry of Education of China No 20110181120090.
文摘In this paper,we propose a new two-level preconditioned C-G method which uses the quadratic smoothing and the linear correction in distorted but topologically structured grid.The CPU time of this method is less than that of the multigrid preconditioned C-G method(MGCG)using the quadratic element,but their accuracy is almost the same.Numerical experiments and eigenvalue analysis are given and the results show that the proposed two-level preconditioned method is efficient.
基金The first author acknowledge the support of the Institute for Advanced Study(IAS)at The Hong Kong University of Science and TechnologyThe work is partially supported by grants from RGC CA05/06.SC01 and RGC-CERG 603107.
文摘In this article,we discuss a least-squares/fictitious domain method for the solution of linear elliptic boundary value problems with Robin boundary conditions.LetΩandωbe two bounded domains of R d such thatω⊂Ω.For a linear elliptic problem inΩ\ωwith Robin boundary condition on the boundaryγofω,our goal here is to develop a fictitious domain method where one solves a variant of the original problem on the fullΩ,followed by a well-chosen correction overω.This method is of the virtual control type and relies on a least-squares formulation making the problem solvable by a conjugate gradient algorithm operating in a well chosen control space.Numerical results obtained when applying our method to the solution of two-dimensional elliptic and parabolic problems are given;they suggest optimal order of convergence.
