In this paper, we investigate the use of ultra weak variational formulation to solve a wave scattering problem in near field optics. In order to capture the sub-scale features of waves, we utilize evanescent wave func...In this paper, we investigate the use of ultra weak variational formulation to solve a wave scattering problem in near field optics. In order to capture the sub-scale features of waves, we utilize evanescent wave functions together with plane wave functions to approximate the local properties of the field. We analyze the global convergence and give an error estimation of the method. Numerical examples are also presented to demonstrate the effectiveness of the strategy.展开更多
We investigate the ultra weak variational formulation(UWVF)of the 2-D Helmholtz equation using a new choice of basis functions.Traditionally the UWVF basis functions are chosen to be plane waves.Here,we instead use fi...We investigate the ultra weak variational formulation(UWVF)of the 2-D Helmholtz equation using a new choice of basis functions.Traditionally the UWVF basis functions are chosen to be plane waves.Here,we instead use first kind Bessel functions.We compare the performance of the two bases.Moreover,we show that it is possible to use coupled plane wave and Bessel bases in the same mesh.As test cases we shall consider propagating plane and evanescent waves in a rectangular domain and a singular 2-D Helmholtz problem in an L-shaped domain.展开更多
An interesting discretization method for Helmholtz equations was introduced in B.Despres[1].This method is based on the ultra weak variational formulation(UWVF)and the wave shape functions,which are exact solutions of...An interesting discretization method for Helmholtz equations was introduced in B.Despres[1].This method is based on the ultra weak variational formulation(UWVF)and the wave shape functions,which are exact solutions of the governing Helmholtz equation.In this paper we are concerned with fast solver for the system generated by the method in[1].We propose a new preconditioner for such system,which can be viewed as a combination between a coarse solver and the block diagonal preconditioner introduced in[13].In our numerical experiments,this preconditioner is applied to solve both two-dimensional and three-dimensional Helmholtz equations,and the numerical results illustrate that the new preconditioner is much more efficient than the original block diagonal preconditioner.展开更多
基金The authors would like to thank the reviewers and Dr.Zheng Enxi for many valuable suggcstions. This work is supported by the National Natural Science Foundation of China (Grant No. 11371172, 51178001), Science and technology research project of the education department of Jilin Province (Grant No. 2014213).
文摘In this paper, we investigate the use of ultra weak variational formulation to solve a wave scattering problem in near field optics. In order to capture the sub-scale features of waves, we utilize evanescent wave functions together with plane wave functions to approximate the local properties of the field. We analyze the global convergence and give an error estimation of the method. Numerical examples are also presented to demonstrate the effectiveness of the strategy.
文摘We investigate the ultra weak variational formulation(UWVF)of the 2-D Helmholtz equation using a new choice of basis functions.Traditionally the UWVF basis functions are chosen to be plane waves.Here,we instead use first kind Bessel functions.We compare the performance of the two bases.Moreover,we show that it is possible to use coupled plane wave and Bessel bases in the same mesh.As test cases we shall consider propagating plane and evanescent waves in a rectangular domain and a singular 2-D Helmholtz problem in an L-shaped domain.
基金The second author was supported by the Major Research Plan of Natural Science Foundation of China G91130015the Key Project of Natural Science Foundation of China G11031006National Basic Research Program of China G2011309702.
文摘An interesting discretization method for Helmholtz equations was introduced in B.Despres[1].This method is based on the ultra weak variational formulation(UWVF)and the wave shape functions,which are exact solutions of the governing Helmholtz equation.In this paper we are concerned with fast solver for the system generated by the method in[1].We propose a new preconditioner for such system,which can be viewed as a combination between a coarse solver and the block diagonal preconditioner introduced in[13].In our numerical experiments,this preconditioner is applied to solve both two-dimensional and three-dimensional Helmholtz equations,and the numerical results illustrate that the new preconditioner is much more efficient than the original block diagonal preconditioner.