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.展开更多
An h-adaptive method is developed for high-order discontinuous Galerkin methods(DGM)to solve the laminar compressible Navier-Stokes(N-S)equations on unstructured mesh.The vorticity is regarded as the indicator of adap...An h-adaptive method is developed for high-order discontinuous Galerkin methods(DGM)to solve the laminar compressible Navier-Stokes(N-S)equations on unstructured mesh.The vorticity is regarded as the indicator of adaptivity.The elements where the vorticity is larger than a pre-defined upper limit are refined,and those where the vorticity is smaller than a pre-defined lower limit are coarsened if they have been refined.A high-order geometric approximation of curved boundaries is adopted to ensure the accuracy.Numerical results indicate that highly accurate numerical results can be obtained with the adaptive method at relatively low expense.展开更多
A high-order leap-frog based non-dissipative discontinuous Galerkin time- domain method for solving Maxwell's equations is introduced and analyzed. The pro- posed method combines a centered approximation for the eval...A high-order leap-frog based non-dissipative discontinuous Galerkin time- domain method for solving Maxwell's equations is introduced and analyzed. The pro- posed method combines a centered approximation for the evaluation of fluxes at the in- terface between neighboring elements, with a Nth-order leap-frog time scheme. More- over, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwelrs equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with high- order elements show the potential of the method.展开更多
Error estimates of Galerkin method for Kuramoto-Sivashingsky (K-S) equation in space dimension ≥3 are derived in the paper. These results furnish strong evidence for the computation of the solutions.
This paper describes a numerical solution for a two-point boundary value problem. It includes an algorithm for discretization by mixed finite element method. The discrete scheme allows the utilization a finite element...This paper describes a numerical solution for a two-point boundary value problem. It includes an algorithm for discretization by mixed finite element method. The discrete scheme allows the utilization a finite element method based on piecewise linear approximating functions and we also use the barycentric quadrature rule to compute the stiffness matrix and the L2-norm.展开更多
The chaotic motion behavior of the rectangular conductive thin plate that is simply supported on four sides by airflow andmechanical external excitation in a magnetic field is studied.According to Kirchhoff’s thin pl...The chaotic motion behavior of the rectangular conductive thin plate that is simply supported on four sides by airflow andmechanical external excitation in a magnetic field is studied.According to Kirchhoff’s thin plate theory,considering geometric nonlinearity and using the principle of virtualwork,the nonlinearmotion partial differential equation of the rectangular conductive thin plate is deduced.Using the separate variable method and Galerkin’s method,the system motion partial differential equation is converted into the general equation of the Duffing equation;the Hamilton system is introduced,and the Melnikov function is used to analyze the Hamilton system,and obtain the critical surface for the existence of chaos.The bifurcation diagram,phase portrait,time history response and Poincarémap of the vibration system are obtained by numerical simulation,and the correctness is demonstrated.The results showthatwhen the ratio of external excitation amplitude to damping coefficient is higher than the critical surface,the system will enter chaotic state.The chaotic motion of the rectangular conductive thin plate is affected by different magnetic field distributions and airflow.展开更多
In this paper, we discuss the existence and uniqueness of global solutions, the existence of the family of global attractors and its dimension estimation for generalized Beam-Kirchhoff equation under initial condition...In this paper, we discuss the existence and uniqueness of global solutions, the existence of the family of global attractors and its dimension estimation for generalized Beam-Kirchhoff equation under initial conditions and boundary conditions, using the previous research results for reference. Firstly, the existence of bounded absorption set is proved by using a prior estimation, then the existence and uniqueness of the global solution of the problem is proved by using the classical Galerkin’s method. Finally, Housdorff dimension and fractal dimension of the family of global attractors are estimated by linear variational method and generalized Sobolev-Lieb-Thirring inequality.展开更多
In this paper, we study the long time behavior of a class of generalized Beam-Kirchhoff equation , and prove the existence and uniqueness of the global solution of this class of equation by Galerkin method by making s...In this paper, we study the long time behavior of a class of generalized Beam-Kirchhoff equation , and prove the existence and uniqueness of the global solution of this class of equation by Galerkin method by making some assumptions about the nonlinear function term . The existence of the family of global attractor and its Hausdorff dimension and Fractal dimension estimation are proved.展开更多
基金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.
基金supported by the National Natural Science Foundation of China(11272152)
文摘An h-adaptive method is developed for high-order discontinuous Galerkin methods(DGM)to solve the laminar compressible Navier-Stokes(N-S)equations on unstructured mesh.The vorticity is regarded as the indicator of adaptivity.The elements where the vorticity is larger than a pre-defined upper limit are refined,and those where the vorticity is smaller than a pre-defined lower limit are coarsened if they have been refined.A high-order geometric approximation of curved boundaries is adopted to ensure the accuracy.Numerical results indicate that highly accurate numerical results can be obtained with the adaptive method at relatively low expense.
基金supported by a grant from the French National Ministry of Education and Research(MENSR,19755-2005)
文摘A high-order leap-frog based non-dissipative discontinuous Galerkin time- domain method for solving Maxwell's equations is introduced and analyzed. The pro- posed method combines a centered approximation for the evaluation of fluxes at the in- terface between neighboring elements, with a Nth-order leap-frog time scheme. More- over, the interpolation degree is defined at the element level and the mesh is refined locally in a non-conforming way resulting in arbitrary level hanging nodes. The method is proved to be stable under some CFL-like condition on the time step. The convergence of the semi-discrete approximation to Maxwelrs equations is established rigorously and bounds on the global divergence error are provided. Numerical experiments with high- order elements show the potential of the method.
文摘Error estimates of Galerkin method for Kuramoto-Sivashingsky (K-S) equation in space dimension ≥3 are derived in the paper. These results furnish strong evidence for the computation of the solutions.
基金The work was supported by the National Natural Science Foundation of China( 1 9971 0 5 7) and by the Youth Sci-ence Foundation of Shanghai Municipal Commission of Education ( 99QA6 6 )
文摘In this paper, the convergence and L 2 estimate of the Galerkin method are given for the steady state Kuramoto-Sivashinsky (K-S) equation.
文摘This paper describes a numerical solution for a two-point boundary value problem. It includes an algorithm for discretization by mixed finite element method. The discrete scheme allows the utilization a finite element method based on piecewise linear approximating functions and we also use the barycentric quadrature rule to compute the stiffness matrix and the L2-norm.
基金funded by the Anhui Provincial Natural Science Foundation(Grant No.2008085QE245)the Natural Science Research Project of Higher Education Institutions in Anhui Province(2022AH040045)the Project of Science and Technology Plan of Department of Housing and Urban-Rural Development of Anhui Province(2021-YF22).
文摘The chaotic motion behavior of the rectangular conductive thin plate that is simply supported on four sides by airflow andmechanical external excitation in a magnetic field is studied.According to Kirchhoff’s thin plate theory,considering geometric nonlinearity and using the principle of virtualwork,the nonlinearmotion partial differential equation of the rectangular conductive thin plate is deduced.Using the separate variable method and Galerkin’s method,the system motion partial differential equation is converted into the general equation of the Duffing equation;the Hamilton system is introduced,and the Melnikov function is used to analyze the Hamilton system,and obtain the critical surface for the existence of chaos.The bifurcation diagram,phase portrait,time history response and Poincarémap of the vibration system are obtained by numerical simulation,and the correctness is demonstrated.The results showthatwhen the ratio of external excitation amplitude to damping coefficient is higher than the critical surface,the system will enter chaotic state.The chaotic motion of the rectangular conductive thin plate is affected by different magnetic field distributions and airflow.
文摘In this paper, we discuss the existence and uniqueness of global solutions, the existence of the family of global attractors and its dimension estimation for generalized Beam-Kirchhoff equation under initial conditions and boundary conditions, using the previous research results for reference. Firstly, the existence of bounded absorption set is proved by using a prior estimation, then the existence and uniqueness of the global solution of the problem is proved by using the classical Galerkin’s method. Finally, Housdorff dimension and fractal dimension of the family of global attractors are estimated by linear variational method and generalized Sobolev-Lieb-Thirring inequality.
文摘In this paper, we study the long time behavior of a class of generalized Beam-Kirchhoff equation , and prove the existence and uniqueness of the global solution of this class of equation by Galerkin method by making some assumptions about the nonlinear function term . The existence of the family of global attractor and its Hausdorff dimension and Fractal dimension estimation are proved.