The Runge-Kutta discontinuous Galerkin finite element method (RK-DGFEM) is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and ...The Runge-Kutta discontinuous Galerkin finite element method (RK-DGFEM) is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and it is explicit in the time domain. Consequently it is a best mixture of FEM and finite volume method (FVM). RK-DGFEM can obtain local high-order accuracy by using high-order polynomial basis. Numerical experiments of transverse magnetic (TM) wave propagation in a 2-D resonator are performed. A high-order Lagrange polynomial basis is adopted. Numerical results agree well with analytical solution. And different order Lagrange interpolation polynomial basis impacts on simulation result accuracy are discussed. Computational results indicate that the accuracy is evidently improved when the order of interpolation basis is increased. Finally, L^2 errors of different order polynomial basis in RK-DGFEM are presented. Computational results show that L^2 error declines exponentially as the order of basis increases.展开更多
Several algorithms were proposed relating to the development of a framework of the perturbation-based stochastic finite element method (PSFEM) for large variation nonlinear dynamic problems. For this purpose, algorith...Several algorithms were proposed relating to the development of a framework of the perturbation-based stochastic finite element method (PSFEM) for large variation nonlinear dynamic problems. For this purpose, algorithms and a framework related to SFEM based on the stochastic virtual work principle were studied. To prove the validity and practicality of the algorithms and framework, numerical examples for nonlinear dynamic problems with large variations were calculated and compared with the Monte-Carlo Simulation method. This comparison shows that the proposed approaches are accurate and effective for the nonlinear dynamic analysis of structures with random parameters.展开更多
Higher-order Time Domain Finite Element Method (TDFEM) based on the nodal inter- polation is proposed for two-dimensional electromagnetic analysis. The detailed algorithms of the method are presented firstly, and then...Higher-order Time Domain Finite Element Method (TDFEM) based on the nodal inter- polation is proposed for two-dimensional electromagnetic analysis. The detailed algorithms of the method are presented firstly, and then the accuracy, CPU time and memory consumption of the higher-order node-based TDFEM are investigated. The high performance of the presented approach is validated by numerical results of the transient responses of Transverse Electric (TE) field and Transverse Magnetic (TM) field in a rectangular waveguide.展开更多
The approach of nonconforming finite element method admits users to solve the partial differential equations with lower complexity,but the accuracy is usually low.In this paper,we present a family of highaccuracy nonc...The approach of nonconforming finite element method admits users to solve the partial differential equations with lower complexity,but the accuracy is usually low.In this paper,we present a family of highaccuracy nonconforming finite element methods for fourth order problems in arbitrary dimensions.The finite element methods are given in a unified way with respect to the dimension.This is an effort to reveal the balance between the accuracy and the complexity of finite element methods.展开更多
In this paper,the discontinuous Galerkin(DG)method combined with localized artificial diffusivity is investigated in the context of numerical simulation of broadband compressible turbulent flows with shocks for under-...In this paper,the discontinuous Galerkin(DG)method combined with localized artificial diffusivity is investigated in the context of numerical simulation of broadband compressible turbulent flows with shocks for under-resolved cases.Firstly,the spectral property of the DG method is analyzed using the approximate dispersion relation(ADR)method and compared with typical finite difference methods,which reveals quantitatively that significantly less grid points can be used with DG for comparable numerical error.Then several typical test cases relevant to problems of compressible turbulence are simulated,including one-dimensional shock/entropy wave interaction,two-dimensional decaying isotropic turbulence,and two-dimensional temporal mixing layers.Numerical results indicate that higher numerical accuracy can be achieved on the same number of degrees of freedom with DG than high order finite difference schemes.Furthermore,shocks are also well captured using the localized artificial diffusivity method.The results in this work can provide useful guidance for further applications of DG to direct and large eddy simulation of compressible turbulent flows.展开更多
文摘The Runge-Kutta discontinuous Galerkin finite element method (RK-DGFEM) is introduced to solve the classical resonator problem in the time domain. DGFEM uses unstructured grid discretization in the space domain and it is explicit in the time domain. Consequently it is a best mixture of FEM and finite volume method (FVM). RK-DGFEM can obtain local high-order accuracy by using high-order polynomial basis. Numerical experiments of transverse magnetic (TM) wave propagation in a 2-D resonator are performed. A high-order Lagrange polynomial basis is adopted. Numerical results agree well with analytical solution. And different order Lagrange interpolation polynomial basis impacts on simulation result accuracy are discussed. Computational results indicate that the accuracy is evidently improved when the order of interpolation basis is increased. Finally, L^2 errors of different order polynomial basis in RK-DGFEM are presented. Computational results show that L^2 error declines exponentially as the order of basis increases.
文摘Several algorithms were proposed relating to the development of a framework of the perturbation-based stochastic finite element method (PSFEM) for large variation nonlinear dynamic problems. For this purpose, algorithms and a framework related to SFEM based on the stochastic virtual work principle were studied. To prove the validity and practicality of the algorithms and framework, numerical examples for nonlinear dynamic problems with large variations were calculated and compared with the Monte-Carlo Simulation method. This comparison shows that the proposed approaches are accurate and effective for the nonlinear dynamic analysis of structures with random parameters.
基金Supported by National Natural Science Foundation of China (No. 60601024)
文摘Higher-order Time Domain Finite Element Method (TDFEM) based on the nodal inter- polation is proposed for two-dimensional electromagnetic analysis. The detailed algorithms of the method are presented firstly, and then the accuracy, CPU time and memory consumption of the higher-order node-based TDFEM are investigated. The high performance of the presented approach is validated by numerical results of the transient responses of Transverse Electric (TE) field and Transverse Magnetic (TM) field in a rectangular waveguide.
基金supported by National Natural Science Foundation of China (Grant No.11101415)the National Center for Mathematics and Interdisciplinary Sciences,CAS
文摘The approach of nonconforming finite element method admits users to solve the partial differential equations with lower complexity,but the accuracy is usually low.In this paper,we present a family of highaccuracy nonconforming finite element methods for fourth order problems in arbitrary dimensions.The finite element methods are given in a unified way with respect to the dimension.This is an effort to reveal the balance between the accuracy and the complexity of finite element methods.
基金supported by the National Basic Research Program of China(Grant No.2009CB724104)
文摘In this paper,the discontinuous Galerkin(DG)method combined with localized artificial diffusivity is investigated in the context of numerical simulation of broadband compressible turbulent flows with shocks for under-resolved cases.Firstly,the spectral property of the DG method is analyzed using the approximate dispersion relation(ADR)method and compared with typical finite difference methods,which reveals quantitatively that significantly less grid points can be used with DG for comparable numerical error.Then several typical test cases relevant to problems of compressible turbulence are simulated,including one-dimensional shock/entropy wave interaction,two-dimensional decaying isotropic turbulence,and two-dimensional temporal mixing layers.Numerical results indicate that higher numerical accuracy can be achieved on the same number of degrees of freedom with DG than high order finite difference schemes.Furthermore,shocks are also well captured using the localized artificial diffusivity method.The results in this work can provide useful guidance for further applications of DG to direct and large eddy simulation of compressible turbulent flows.