We provide an introduction to the use of the spectral-elementmethod (SEM)in seismology. Following a brief review of the basic equations that govern seismicwave propagation, we discuss in some detail how these equation...We provide an introduction to the use of the spectral-elementmethod (SEM)in seismology. Following a brief review of the basic equations that govern seismicwave propagation, we discuss in some detail how these equations may be solved numericallybased upon the SEM to address the forward problem in seismology. Examplesof synthetic seismograms calculated based upon the SEM are compared to datarecorded by the Global Seismographic Network. Finally, we discuss the challenge ofusing the remaining differences between the data and the synthetic seismograms toconstrain better Earth models and source descriptions. This leads naturally to adjointmethods, which provide a practical approach to this formidable computational challengeand enables seismologists to tackle the inverse problem.展开更多
We propose a novel numerical approach for delay differential equations with vanishing proportional delays based on spectral methods. A Legendre-collocation method is employed to obtain highly accurate numerical approx...We propose a novel numerical approach for delay differential equations with vanishing proportional delays based on spectral methods. A Legendre-collocation method is employed to obtain highly accurate numerical approximations to the exact solution. It is proved theoretically and demonstrated numerically that the proposed method converges exponentially provided that the data in the are smooth. given pantograph delay differential equation展开更多
A numerical algorithm for effective incorporation of parametric uncertainty into mathematical models is presented.The uncertain parameters are modeled as random variables,and the governing equations are treated as sto...A numerical algorithm for effective incorporation of parametric uncertainty into mathematical models is presented.The uncertain parameters are modeled as random variables,and the governing equations are treated as stochastic.The solutions,or quantities of interests,are expressed as convergent series of orthogonal polynomial expansions in terms of the input random parameters.A high-order stochastic collocation method is employed to solve the solution statistics,and more importantly,to reconstruct the polynomial expansion.While retaining the high accuracy by polynomial expansion,the resulting“pseudo-spectral”type algorithm is straightforward to implement as it requires only repetitive deterministic simulations.An estimate on error bounded is presented,along with numerical examples for problems with relatively complicated forms of governing equations.展开更多
The spectral methods and ice-induced fatigue analysis are discussed based on Miner＇s linear cumulative fatigue hypothesis and S-N curve data. According to the long-term data of full-scale tests on the platforms in th...The spectral methods and ice-induced fatigue analysis are discussed based on Miner＇s linear cumulative fatigue hypothesis and S-N curve data. According to the long-term data of full-scale tests on the platforms in the Bohai Sea, the ice force spectrum of conical structures and the fatigue environmental model are established. Moreover, the finite element model of JZ20-2MSW platform, an example of ice-induced fatigue analysis, is built with ANSYS software. The mode analysis and dynamic analysis in frequency domain under all kinds of ice fatigue work conditions are carded on, and the fatigue life of the structure is estimated in detail. The methods in this paper can be helpful in ice-induced fatigue analysis of ice-resistant platforms.展开更多
Raman technology,which covers Raman spectroscopy(RS)and its various derivativemethods,has been widely applied in detection of various substances in agriculture,food and biosystems.This article reviews the recent advan...Raman technology,which covers Raman spectroscopy(RS)and its various derivativemethods,has been widely applied in detection of various substances in agriculture,food and biosystems.This article reviews the recent advances in two mainstream Raman technologies as RS and SERS,including technical evolution,application and challenges,and spectral processing methods.Firstly,the origin,principle,defect,and development of RSwere introduced.Then,the current situation,existing problems,and development trend of RS and SERS were discussed in agriculture,food,and biosystems,such as adulteration recognition,plant diseases identification,farm chemicals detection,food additives determination,and toxins analysis.At last,the spectral analysis methods include noise reduction,feature extraction or variable selection,and modelingwere introduced in detail,which can realize the automatic and intelligent analysis of spectra without relying on professionals.展开更多
In this paper,we investigate initial boundary value problems of the spacetime fractional diffusion equation and its numerical solutions.Two definitions,i.e.,Riemann-Liouville definition and Caputo one,of the fractiona...In this paper,we investigate initial boundary value problems of the spacetime fractional diffusion equation and its numerical solutions.Two definitions,i.e.,Riemann-Liouville definition and Caputo one,of the fractional derivative are considered in parallel.In both cases,we establish the well-posedness of the weak solution.Moveover,based on the proposed weak formulation,we construct an efficient spectral method for numerical approximations of the weak solution.The main contribution of this work are threefold:First,a theoretical framework for the variational solutions of the space-time fractional diffusion equation is developed.We find suitable functional spaces and norms in which the space-time fractional diffusion problem can be formulated into an elliptic weak problem,and the existence and uniqueness of the weak solution are then proved by using existing theory for elliptic problems.Secondly,we show that in the case of Riemann-Liouville definition,the well-posedness of the space-time fractional diffusion equation does not require any initial conditions.This contrasts with the case of Caputo definition,in which the initial condition has to be integrated into the weak formulation in order to establish the well-posedness.Finally,thanks to the weak formulation,we are able to construct an efficient numerical method for solving the space-time fractional diffusion problem.展开更多
The asymptotic behavior of the solutions to a class of pseudoparabolic viscous diffusion equation with periodic initial condition is studied by using the spectral method. The semidiscrete Fourier approximate solution ...The asymptotic behavior of the solutions to a class of pseudoparabolic viscous diffusion equation with periodic initial condition is studied by using the spectral method. The semidiscrete Fourier approximate solution of the problem is constructed and the error estimation between spectral approximate solution and exact solution on large time is also obtained. The existence of the approximate attractor AN and the upper semicontinuity d（AN,A） → 0 are proved.展开更多
In the author’s recent publications, a parametric system biorthogonal to the corresponding segment of the exponential Fourier system was unusually effective. On its basis, it was discovered that knowledge of a finite...In the author’s recent publications, a parametric system biorthogonal to the corresponding segment of the exponential Fourier system was unusually effective. On its basis, it was discovered that knowledge of a finite number of Fourier coefficients of function f from an infinite-dimensional set of elementary functions allows f to be accurately restored (the phenomenon of over-convergence). Below, parametric biorthogonal systems are constructed for classical trigonometric Fourier series, and the corresponding phenomena of over-convergence are discovered. The decisive role here was played by representing the space L2 as an orthogonal sum of two corresponding subspaces. As a result, fast parallel algorithms for reconstructing a function from its truncated trigonometric Fourier series are proposed. The presented numerical experiments confirm the high efficiency of these convergence accelerations for smooth functions. In conclusion, the main results of the work are summarized, and some prospects for the development and generalization of the proposed approaches are discussed.展开更多
This paper presents a review of the current state-of-the-art of numerical methods for stochastic computations.The focus is on efficient high-order methods suitable for practical applications,with a particular emphasis...This paper presents a review of the current state-of-the-art of numerical methods for stochastic computations.The focus is on efficient high-order methods suitable for practical applications,with a particular emphasis on those based on generalized polynomial chaos(gPC)methodology.The framework of gPC is reviewed,along with its Galerkin and collocation approaches for solving stochastic equations.Properties of these methods are summarized by using results from literature.This paper also attempts to present the gPC based methods in a unified framework based on an extension of the classical spectral methods into multi-dimensional random spaces.展开更多
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.展开更多
We describe the application of the spectral method to delay integro-differential equations with proportional delays. It is shown that the resulting numerical solutions exhibit the spectral convergence order. Extension...We describe the application of the spectral method to delay integro-differential equations with proportional delays. It is shown that the resulting numerical solutions exhibit the spectral convergence order. Extensions to equations with more general （nonlinear） vanishing delays are also discussed.展开更多
The theory of a class of spectral methods is extended to Volterra integrodifferential equations which contain a weakly singular kernel(t−s)^(−μ) with 0<μ<1.In this work,we consider the case when the underlying...The theory of a class of spectral methods is extended to Volterra integrodifferential equations which contain a weakly singular kernel(t−s)^(−μ) with 0<μ<1.In this work,we consider the case when the underlying solutions of weakly singular Volterra integro-differential equations are sufficiently smooth.We provide a rigorous error analysis for the spectral methods,which shows that both the errors of approximate solutions and the errors of approximate derivatives of the solutions decay exponentially in L^(∞)-norm and weighted L^(2)-norm.The numerical examples are given to illustrate the theoretical results.展开更多
We extend the monolithic convex limiting(MCL)methodology to nodal discontinuous Galerkin spectral-element methods(DGSEMS).The use of Legendre-Gauss-Lobatto(LGL)quadrature endows collocated DGSEM space discretizations ...We extend the monolithic convex limiting(MCL)methodology to nodal discontinuous Galerkin spectral-element methods(DGSEMS).The use of Legendre-Gauss-Lobatto(LGL)quadrature endows collocated DGSEM space discretizations of nonlinear hyperbolic problems with properties that greatly simplify the design of invariant domain-preserving high-resolution schemes.Compared to many other continuous and discontinuous Galerkin method variants,a particular advantage of the LGL spectral operator is the availability of a natural decomposition into a compatible subcellflux discretization.Representing a highorder spatial semi-discretization in terms of intermediate states,we performflux limiting in a manner that keeps these states and the results of Runge-Kutta stages in convex invariant domains.In addition,local bounds may be imposed on scalar quantities of interest.In contrast to limiting approaches based on predictor-corrector algorithms,our MCL procedure for LGL-DGSEM yields nonlinearflux approximations that are independent of the time-step size and can be further modified to enforce entropy stability.To demonstrate the robustness of MCL/DGSEM schemes for the compressible Euler equations,we run simulations for challenging setups featuring strong shocks,steep density gradients,and vortex dominatedflows.展开更多
基金Broadband data were obtained from the IRIS Data Management Center.
文摘We provide an introduction to the use of the spectral-elementmethod (SEM)in seismology. Following a brief review of the basic equations that govern seismicwave propagation, we discuss in some detail how these equations may be solved numericallybased upon the SEM to address the forward problem in seismology. Examplesof synthetic seismograms calculated based upon the SEM are compared to datarecorded by the Global Seismographic Network. Finally, we discuss the challenge ofusing the remaining differences between the data and the synthetic seismograms toconstrain better Earth models and source descriptions. This leads naturally to adjointmethods, which provide a practical approach to this formidable computational challengeand enables seismologists to tackle the inverse problem.
基金The research of HB was supported by the Natural Sciences and Engineering Research Council (NSERC) of Canada and by the Research Grants Council of Hong KongThe research of TT was supported by Hong Kong Baptist University,the Research Grants Council of Hong Kong and he was supported in part by the Chinese Academy of Sciences while visiting its Institute of Computational Mathematics.
文摘We propose a novel numerical approach for delay differential equations with vanishing proportional delays based on spectral methods. A Legendre-collocation method is employed to obtain highly accurate numerical approximations to the exact solution. It is proved theoretically and demonstrated numerically that the proposed method converges exponentially provided that the data in the are smooth. given pantograph delay differential equation
文摘A numerical algorithm for effective incorporation of parametric uncertainty into mathematical models is presented.The uncertain parameters are modeled as random variables,and the governing equations are treated as stochastic.The solutions,or quantities of interests,are expressed as convergent series of orthogonal polynomial expansions in terms of the input random parameters.A high-order stochastic collocation method is employed to solve the solution statistics,and more importantly,to reconstruct the polynomial expansion.While retaining the high accuracy by polynomial expansion,the resulting“pseudo-spectral”type algorithm is straightforward to implement as it requires only repetitive deterministic simulations.An estimate on error bounded is presented,along with numerical examples for problems with relatively complicated forms of governing equations.
基金The paper was supported by the National 863 High Technology Develpoment Plan Project(Grant No.2001AA602015)
文摘The spectral methods and ice-induced fatigue analysis are discussed based on Miner＇s linear cumulative fatigue hypothesis and S-N curve data. According to the long-term data of full-scale tests on the platforms in the Bohai Sea, the ice force spectrum of conical structures and the fatigue environmental model are established. Moreover, the finite element model of JZ20-2MSW platform, an example of ice-induced fatigue analysis, is built with ANSYS software. The mode analysis and dynamic analysis in frequency domain under all kinds of ice fatigue work conditions are carded on, and the fatigue life of the structure is estimated in detail. The methods in this paper can be helpful in ice-induced fatigue analysis of ice-resistant platforms.
基金This study is supported by Anhui Provincial Major Scientific and Technological Special Project(No.17030701062)National Natural Science Foundation of China(Nos.3170123 and 61672032)+1 种基金National Key Research and Development Program(No.2016YFD0800904)Open Foundation of Laboratory of Quality and Safety Risk Assessment on Agricultural products(Beijing),Ministry of Agriculture(No.KFRA201802).
文摘Raman technology,which covers Raman spectroscopy(RS)and its various derivativemethods,has been widely applied in detection of various substances in agriculture,food and biosystems.This article reviews the recent advances in two mainstream Raman technologies as RS and SERS,including technical evolution,application and challenges,and spectral processing methods.Firstly,the origin,principle,defect,and development of RSwere introduced.Then,the current situation,existing problems,and development trend of RS and SERS were discussed in agriculture,food,and biosystems,such as adulteration recognition,plant diseases identification,farm chemicals detection,food additives determination,and toxins analysis.At last,the spectral analysis methods include noise reduction,feature extraction or variable selection,and modelingwere introduced in detail,which can realize the automatic and intelligent analysis of spectra without relying on professionals.
文摘In this paper,we investigate initial boundary value problems of the spacetime fractional diffusion equation and its numerical solutions.Two definitions,i.e.,Riemann-Liouville definition and Caputo one,of the fractional derivative are considered in parallel.In both cases,we establish the well-posedness of the weak solution.Moveover,based on the proposed weak formulation,we construct an efficient spectral method for numerical approximations of the weak solution.The main contribution of this work are threefold:First,a theoretical framework for the variational solutions of the space-time fractional diffusion equation is developed.We find suitable functional spaces and norms in which the space-time fractional diffusion problem can be formulated into an elliptic weak problem,and the existence and uniqueness of the weak solution are then proved by using existing theory for elliptic problems.Secondly,we show that in the case of Riemann-Liouville definition,the well-posedness of the space-time fractional diffusion equation does not require any initial conditions.This contrasts with the case of Caputo definition,in which the initial condition has to be integrated into the weak formulation in order to establish the well-posedness.Finally,thanks to the weak formulation,we are able to construct an efficient numerical method for solving the space-time fractional diffusion problem.
基金This work was supported by the National Science Foundation of China(10271034)
文摘The asymptotic behavior of the solutions to a class of pseudoparabolic viscous diffusion equation with periodic initial condition is studied by using the spectral method. The semidiscrete Fourier approximate solution of the problem is constructed and the error estimation between spectral approximate solution and exact solution on large time is also obtained. The existence of the approximate attractor AN and the upper semicontinuity d（AN,A） → 0 are proved.
文摘In the author’s recent publications, a parametric system biorthogonal to the corresponding segment of the exponential Fourier system was unusually effective. On its basis, it was discovered that knowledge of a finite number of Fourier coefficients of function f from an infinite-dimensional set of elementary functions allows f to be accurately restored (the phenomenon of over-convergence). Below, parametric biorthogonal systems are constructed for classical trigonometric Fourier series, and the corresponding phenomena of over-convergence are discovered. The decisive role here was played by representing the space L2 as an orthogonal sum of two corresponding subspaces. As a result, fast parallel algorithms for reconstructing a function from its truncated trigonometric Fourier series are proposed. The presented numerical experiments confirm the high efficiency of these convergence accelerations for smooth functions. In conclusion, the main results of the work are summarized, and some prospects for the development and generalization of the proposed approaches are discussed.
基金This research is supported in part by NSF CAREER Award DMS-0645035.
文摘This paper presents a review of the current state-of-the-art of numerical methods for stochastic computations.The focus is on efficient high-order methods suitable for practical applications,with a particular emphasis on those based on generalized polynomial chaos(gPC)methodology.The framework of gPC is reviewed,along with its Galerkin and collocation approaches for solving stochastic equations.Properties of these methods are summarized by using results from literature.This paper also attempts to present the gPC based methods in a unified framework based on an extension of the classical spectral methods into multi-dimensional random spaces.
基金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.
文摘We describe the application of the spectral method to delay integro-differential equations with proportional delays. It is shown that the resulting numerical solutions exhibit the spectral convergence order. Extensions to equations with more general （nonlinear） vanishing delays are also discussed.
基金This work is supported by the Foundation for Talent Introduction of Guangdong Provincial University,Guangdong Province Universities and Colleges Pearl River Scholar Funded Scheme(2008)National Science Foundation of China(10971074).
文摘The theory of a class of spectral methods is extended to Volterra integrodifferential equations which contain a weakly singular kernel(t−s)^(−μ) with 0<μ<1.In this work,we consider the case when the underlying solutions of weakly singular Volterra integro-differential equations are sufficiently smooth.We provide a rigorous error analysis for the spectral methods,which shows that both the errors of approximate solutions and the errors of approximate derivatives of the solutions decay exponentially in L^(∞)-norm and weighted L^(2)-norm.The numerical examples are given to illustrate the theoretical results.
文摘We extend the monolithic convex limiting(MCL)methodology to nodal discontinuous Galerkin spectral-element methods(DGSEMS).The use of Legendre-Gauss-Lobatto(LGL)quadrature endows collocated DGSEM space discretizations of nonlinear hyperbolic problems with properties that greatly simplify the design of invariant domain-preserving high-resolution schemes.Compared to many other continuous and discontinuous Galerkin method variants,a particular advantage of the LGL spectral operator is the availability of a natural decomposition into a compatible subcellflux discretization.Representing a highorder spatial semi-discretization in terms of intermediate states,we performflux limiting in a manner that keeps these states and the results of Runge-Kutta stages in convex invariant domains.In addition,local bounds may be imposed on scalar quantities of interest.In contrast to limiting approaches based on predictor-corrector algorithms,our MCL procedure for LGL-DGSEM yields nonlinearflux approximations that are independent of the time-step size and can be further modified to enforce entropy stability.To demonstrate the robustness of MCL/DGSEM schemes for the compressible Euler equations,we run simulations for challenging setups featuring strong shocks,steep density gradients,and vortex dominatedflows.
