In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic ...In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.展开更多
In this paper,we develop a new sixth-order WENO scheme by adopting a convex combina-tion of a sixth-order global reconstruction and four low-order local reconstructions.Unlike the classical WENO schemes,the associated...In this paper,we develop a new sixth-order WENO scheme by adopting a convex combina-tion of a sixth-order global reconstruction and four low-order local reconstructions.Unlike the classical WENO schemes,the associated linear weights of the new scheme can be any positive numbers with the only requirement that their sum equals one.Further,a very simple smoothness indicator for the global stencil is proposed.The new scheme can achieve sixth-order accuracy in smooth regions.Numerical tests in some one-and two-dimensional bench-mark problems show that the new scheme has a little bit higher resolution compared with the recently developed sixth-order WENO-Z6 scheme,and it is more efficient than the classical fifth-order WENO-JS5 scheme and the recently developed sixth-order WENO6-S scheme.展开更多
This paper presents an anisotropic adaptive finite element method (FEM) to solve the governing equations of steady magnetohydrodynamic (MHD) duct flow. A resid- ual error estimator is presented for the standard FE...This paper presents an anisotropic adaptive finite element method (FEM) to solve the governing equations of steady magnetohydrodynamic (MHD) duct flow. A resid- ual error estimator is presented for the standard FEM, and two-sided bounds on the error independent of the aspect ratio of meshes are provided. Based on the Zienkiewicz-Zhu es- timates, a computable anisotropic error indicator and an implement anisotropic adaptive refinement for the MHD problem are derived at different values of the Hartmann number. The most distinguishing feature of the method is that the layer information from some directions is captured well such that the number of mesh vertices is dramatically reduced for a given level of accuracy. Thus, this approach is more suitable for approximating the layer problem at high Hartmann numbers. Numerical results show efficiency of the algorithm.展开更多
基金supported by China Postdoctoral Science Foundation grant 2020TQ0344the NSFC grants 11871139 and 12101597the NSF grants DMS-1720116,DMS-2012882,DMS-2011838,DMS-1719942,DMS-1913072.
文摘In this work,we develop energy stable numerical methods to simulate electromagnetic waves propagating in optical media where the media responses include the linear Lorentz dispersion,the instantaneous nonlinear cubic Kerr response,and the nonlinear delayed Raman molecular vibrational response.Unlike the first-order PDE-ODE governing equations considered previously in Bokil et al.(J Comput Phys 350:420–452,2017)and Lyu et al.(J Sci Comput 89:1–42,2021),a model of mixed-order form is adopted here that consists of the first-order PDE part for Maxwell’s equations coupled with the second-order ODE part(i.e.,the auxiliary differential equations)modeling the linear and nonlinear dispersion in the material.The main contribution is a new numerical strategy to treat the Kerr and Raman nonlinearities to achieve provable energy stability property within a second-order temporal discretization.A nodal discontinuous Galerkin(DG)method is further applied in space for efficiently handling nonlinear terms at the algebraic level,while preserving the energy stability and achieving high-order accuracy.Indeed with d_(E)as the number of the components of the electric field,only a d_(E)×d_(E)nonlinear algebraic system needs to be solved at each interpolation node,and more importantly,all these small nonlinear systems are completely decoupled over one time step,rendering very high parallel efficiency.We evaluate the proposed schemes by comparing them with the methods in Bokil et al.(2017)and Lyu et al.(2021)(implemented in nodal form)regarding the accuracy,computational efficiency,and energy stability,by a parallel scalability study,and also through the simulations of the soliton-like wave propagation in one dimension,as well as the spatial-soliton propagation and two-beam interactions modeled by the two-dimensional transverse electric(TE)mode of the equations.
基金the National Natural Science Foundation of China(91641107,91852116,12071470)Fundamental Research of Civil Aircraft(MJ-F-2012-04)of Ministry of Industrialization and Information of China.
文摘In this paper,we develop a new sixth-order WENO scheme by adopting a convex combina-tion of a sixth-order global reconstruction and four low-order local reconstructions.Unlike the classical WENO schemes,the associated linear weights of the new scheme can be any positive numbers with the only requirement that their sum equals one.Further,a very simple smoothness indicator for the global stencil is proposed.The new scheme can achieve sixth-order accuracy in smooth regions.Numerical tests in some one-and two-dimensional bench-mark problems show that the new scheme has a little bit higher resolution compared with the recently developed sixth-order WENO-Z6 scheme,and it is more efficient than the classical fifth-order WENO-JS5 scheme and the recently developed sixth-order WENO6-S scheme.
基金Project supported by the National Natural Science Foundation of China(Nos.11471329,11321061,and 91430215)the National Magnetic Confinement Fusion Science Program of China(No.2015GB110000)+1 种基金the Youth Innovation Promotion Association of Chinese Academy of Sciences(CAS)(No.2016003)the National Center for Mathematics and Interdisciplinary Sciences of CAS
文摘This paper presents an anisotropic adaptive finite element method (FEM) to solve the governing equations of steady magnetohydrodynamic (MHD) duct flow. A resid- ual error estimator is presented for the standard FEM, and two-sided bounds on the error independent of the aspect ratio of meshes are provided. Based on the Zienkiewicz-Zhu es- timates, a computable anisotropic error indicator and an implement anisotropic adaptive refinement for the MHD problem are derived at different values of the Hartmann number. The most distinguishing feature of the method is that the layer information from some directions is captured well such that the number of mesh vertices is dramatically reduced for a given level of accuracy. Thus, this approach is more suitable for approximating the layer problem at high Hartmann numbers. Numerical results show efficiency of the algorithm.