The manuscript introduces an “ab initio” quantum model to deduce the Maxwell equations. After general considerations and laying out the model’s theoretical framework, these equations can be derived alongside a broa...The manuscript introduces an “ab initio” quantum model to deduce the Maxwell equations. After general considerations and laying out the model’s theoretical framework, these equations can be derived alongside a broad variety of other results. Specifically, a corollary of the present model proposes a possible mechanism underlying the formation of magnetic monopoles and allows estimating their formation energy in order of magnitude.展开更多
Several new energy identities of the two dimenslonal(2D) Maxwell equations in a lossy medium in the case of the perfectly electric conducting boundary conditions are proposed and proved. These identities show a new ...Several new energy identities of the two dimenslonal(2D) Maxwell equations in a lossy medium in the case of the perfectly electric conducting boundary conditions are proposed and proved. These identities show a new kind of energy conservation in the Maxwell system and provide a new energy method to analyze the alternating direction im- plicit finite difference time domain method for the 2D Maxwell equations (2D-ADI-FDTD). It is proved that 2D-ADI-FDTD is approximately energy conserved, unconditionally sta- ble and second order convergent in the discrete L2 and H1 norms, which implies that 2D-ADI-FDTD is super convergent. By this super convergence, it is simply proved that the error of the divergence of the solution of 2D-ADI-FDTD is second order accurate. It is also proved that the difference scheme of 2D-ADI-FDTD with respect to time t is second order convergent in the discrete H1 norm. Experimental results to confirm the theoretical analysis on stability, convergence and energy conservation are presented.展开更多
This paper presents a method to find Noether-type conserved quantities and Lie point symmetries for discrete mechanico-electrical dynamical systems,which leave invuriant the set of solutions of the corresponding diffe...This paper presents a method to find Noether-type conserved quantities and Lie point symmetries for discrete mechanico-electrical dynamical systems,which leave invuriant the set of solutions of the corresponding difference scheme. This approach makes it possible to devise techniques for solving the Lagrange Maxwell equations in differences which correspond to mechanico-electrical systems,by adapting existing differential equations.In particular,it obtains a new systematic method to determine both the one-parameter Lie groups and the discrete Noether conserved quantities of Lie point symmetries for mechanico-electrical systems.As an application,it obtains the Lie point symmetries and the conserved quantities for the difference equation of a model that represents a capacitor microphone.展开更多
In this paper, a unified model for time-dependent Maxwell equations in dispersive media is considered. The space-time DG method developed in [29] is applied to solve the un-derlying problem. Unconditional L2-stability...In this paper, a unified model for time-dependent Maxwell equations in dispersive media is considered. The space-time DG method developed in [29] is applied to solve the un-derlying problem. Unconditional L2-stability and error estimate of order O?τr+1+hk+1/2? are obtained when polynomials of degree at most r and k are used for the temporal dis-cretization and spatial discretization respectively. 2-D and 3-D numerical examples are given to validate the theoretical results. Moreover, numerical results show an ultra-convergence of order 2r+1 in temporal variable t.展开更多
In this article,we investigate a stochastic Galerkin method for the Maxwell equations with random inputs.The generalized Polynomial Chaos(gPC)expansion technique is used to obtain a deterministic system of the gPC exp...In this article,we investigate a stochastic Galerkin method for the Maxwell equations with random inputs.The generalized Polynomial Chaos(gPC)expansion technique is used to obtain a deterministic system of the gPC expansion coefficients.The regularity of the solution with respect to the random is analyzed.On the basis of the regularity results,the optimal convergence rate of the stochastic Galerkin approach for Maxwell equations with random inputs is proved.Numerical examples are presented to support the theoretical analysis.展开更多
This paper is concerned with estimation of electrical conductivity in Maxwell equations. The primary difficulty lies in the presence of numerous local minima in the objective functional. A wavelet multiscale method is...This paper is concerned with estimation of electrical conductivity in Maxwell equations. The primary difficulty lies in the presence of numerous local minima in the objective functional. A wavelet multiscale method is introduced and applied to the inversion of Maxwell equations. The inverse problem is decomposed into multiple scales with wavelet transform, and hence the original problem is reformulated to a set of sub-inverse problems corresponding to different scales, which can be solved successively according to the size of scale from the shortest to the longest. The stable and fast regularized Gauss-Newton method is applied to each scale. Numerical results show that the proposed method is effective, especially in terms of wide convergence, computational efficiency and precision.展开更多
An FDTD system associated with uniaxial perfectly matched layer(UPML) for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-...An FDTD system associated with uniaxial perfectly matched layer(UPML) for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-boundary value problems of the transverse magnetic(TM) mode to Maxwell's equations is obtained by Yee's algorithm,and the open domain of the scattering problem is truncated by a circle with a UPML.Besides,an artificial boundary condition is imposed on the outer boundary of the UPML.Afterwards,stability of the FDTD system on the truncated domain is established through energy estimates by the Gronwall inequality.Numerical experiments are designed to approve the theoretical analysis.展开更多
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.展开更多
In the articles “Newtons Law of Universal Gravitation Explained by the Theory of Informatons” and “The Gravitational Interaction between Moving Mass Particles Explained by the Theory of Informatons” the gravitatio...In the articles “Newtons Law of Universal Gravitation Explained by the Theory of Informatons” and “The Gravitational Interaction between Moving Mass Particles Explained by the Theory of Informatons” the gravitational interaction has been explained by the hypothesis that information carried by informatons is the substance of gravitational fields, i.e. the medium that the interaction in question makes possible. From the idea that “information carried by informatons” is its substance, it has been deduced that—on the macroscopic level—a gravitational field manifests itself as a dual entity, always having a field- and an induction component (Egand Bg) simultaneously created by their common sources. In this article we will mathematically deduce the Maxwell-Heaviside equations from the kinematics of the informatons. These relations describe on the macroscopic level how a gravitational field (Eg, Bg) is generated by whether or not moving masses and how spatial and temporal changes of Egand Bgare related. We show that there is no causal link between Egand Bg.展开更多
The study of Electromagnetic Compatibility is essential to ensure the harmonious operation of electronic equipment in a shared environment. The basic principles of Electromagnetic Compatibility focus on the ability of...The study of Electromagnetic Compatibility is essential to ensure the harmonious operation of electronic equipment in a shared environment. The basic principles of Electromagnetic Compatibility focus on the ability of devices to withstand electromagnetic disturbances and not produce disturbances that could affect other systems. Imperceptible in most work situations, electromagnetic fields can, beyond certain thresholds, have effects on human health. The objective of the present article is focused on the modeling analysis of the influence of geometric parameters of industrial static converters radiated electromagnetic fields using Maxwell’s equations. To do this we used the analytical formalism for calculating the electromagnetic field emitted by a filiform conductor, to model the electromagnetic radiation of this device in the spatio-temporal domain. The interactions of electromagnetic waves with human bodies are complex and depend on several factors linked to the characteristics of the incident wave. To model these interactions, we implemented the physical laws of electromagnetic wave propagation based on Maxwell’s and bio-heat equations to obtain consistent results. These obtained models allowed us to evaluate the spatial profile of induced current and temperature of biological tissue during exposure to electromagnetic waves generated by this system. The simulation 2D results obtained from computer tools show that the temperature variation and current induced by the electromagnetic field can have a very significant influence on the life of biological tissue. The paper provides a comprehensive analysis using advanced mathematical models to evaluate the influence of electromagnetic fields. The findings have direct implications for workplace safety, potentially influencing standards and regulations concerning electromagnetic exposure in industrial settings.展开更多
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.展开更多
We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filament...We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma. We examined the performance of the applied scheme, in this context, we implemented the developed model to study selected phenomena in terahertz radiation production, such as the excitation energy and conversion efficiency of the produced THz radiation, in addition to the influence of the pulse chirping on properties of the produced radiation. The obtained numerical results have clarified that the applied HO-FDTD scheme is precisely accurate to solve Maxwell’s equations and sufficiently valid to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma.展开更多
Maxwell’s equations in electromagnetism can be categorized into three dis-tinct groups based on the electromagnetic source when employing quaterni-ons. Each group represents a self-contained system in which Maxwell’...Maxwell’s equations in electromagnetism can be categorized into three dis-tinct groups based on the electromagnetic source when employing quaterni-ons. Each group represents a self-contained system in which Maxwell’s equations are applied and validated concurrently, in contrast to the previous approach that did not account for this. It has been noted that the formulation of these Maxwell equations ultimately results in the formulation of Max-well’s equations utilizing the scalar function.展开更多
Superconvergence of the mixed finite element methods for 2-d Maxwell equations is studied in this paper. Two order of superconvergent factor can be obtained for the k-th Nedelec elements on the rectangular meshes.
This paper is concerned with the stability and superconvergence analysis of the famous finite-difference time-domain (FDTD) scheme for the 2D Maxwell equations in a lossy medium with a perfectly electric conducting (P...This paper is concerned with the stability and superconvergence analysis of the famous finite-difference time-domain (FDTD) scheme for the 2D Maxwell equations in a lossy medium with a perfectly electric conducting (PEC) boundary condition, employing the energy method. To this end, we first establish some new energy identities for the 2D Maxwell equations in a lossy medium with a PEC boundary condition. Then by making use of these energy identities, it is proved that the FDTD scheme and its time difference scheme are stable in the discrete L2 and H1 norms when the CFL condition is satisfied. It is shown further that the solution to both the FDTD scheme and its time difference scheme is second-order convergent in both space and time in the discrete L2 and H1 norms under a slightly stricter condition than the CFL condition. This means that the solution to the FDTD scheme is superconvergent. Numerical results are also provided to confirm the theoretical analysis.展开更多
We survey recent results on ground and bound state solutions E:?→R^3 of the problem {▽(▽×E)+}λE=|E|^(P-2)E in Ω,v×E=0 on Ω on a bounded Lipschitz domain ??R^3,where?×denotes the curl operator in R...We survey recent results on ground and bound state solutions E:?→R^3 of the problem {▽(▽×E)+}λE=|E|^(P-2)E in Ω,v×E=0 on Ω on a bounded Lipschitz domain ??R^3,where?×denotes the curl operator in R^3.The equation describes the propagation of the time-harmonic electric field R{E(χ)e^(iwt)}in a nonlinear isotropic material ? withλ=-μεω~2≤0,where μ andεstand for the permeability and the linear part of the permittivity of the material.The nonlinear term|E|^(P-2)E with 2<p≤2*=6 comes from the nonlinear polarization and the boundary conditions are those for?surrounded by a perfect conductor.The problem has a variational structure;however the energy functional associated with the problem is strongly indefinite and does not satisfy the Palais-Smale condition.We show the underlying difficulties of the problem and enlist some open questions.展开更多
Extending the holographic program of our previous work,we derive f(R) gravity and the Maxwell equations from the holographic principle,using time-like holographic screens.We find that to derive the Einstein equations ...Extending the holographic program of our previous work,we derive f(R) gravity and the Maxwell equations from the holographic principle,using time-like holographic screens.We find that to derive the Einstein equations and f(R) gravity by a natural holographic approach,the quasi-static condition is necessary.We also find the surface stress tensor and the surface electric current,surface magnetic current on a holographic screen for f(R) gravity and Maxwell's theory,respectively.展开更多
The behavior of interacting electrons in a perfect crystal under macroscopic external electric and magnetic fields is studied. Effective Maxwell equations for the macroscopic electric and magnetic fields are derived s...The behavior of interacting electrons in a perfect crystal under macroscopic external electric and magnetic fields is studied. Effective Maxwell equations for the macroscopic electric and magnetic fields are derived starting from time-dependent density functional theory. Effective permittivity and permeability coefficients are obtained.展开更多
A new multi-dimensional scheme for the Maxwell equations is established by the CIP method in combination with the method of characteristics(CIP-MOC).In addition,the CIP-MOC can be extended to arbitrary grid system by ...A new multi-dimensional scheme for the Maxwell equations is established by the CIP method in combination with the method of characteristics(CIP-MOC).In addition,the CIP-MOC can be extended to arbitrary grid system by the Soroban grid without losing the third-order accuracy.With the accuracy fixed,the grid points required for the CIP are 40 times less than the conventional schemes like the FDTD in three dimensions.Numerical solutions obtained by the CIP-MOC are compared with analytical solution and the FDTD in plane-wave scattering by a perfectly-conducting circular cylinder,and the CIP-MOC agrees very well with analytical solutions.The Soroban grid is also applied to the Vlasov equation that describes the kinematics of plasmas that is frequently combined with the Maxwell equation.The adaptively moving points in velocity space are similar to the particle codes but can provide accurate solutions.展开更多
In this work, a stable numerical algorithm proposed by Chung et al. for the time-domain Maxwell equations is generalized. The time-domain Maxwell equations are solved by expressing the transient behaviors in terms of ...In this work, a stable numerical algorithm proposed by Chung et al. for the time-domain Maxwell equations is generalized. The time-domain Maxwell equations are solved by expressing the transient behaviors in terms of the modified Laguerre polynomials, and then the original equations of the initial value and boundary value can be transformed into a series of problems independent of the time variable. In this case the method of finite difference (FD), the finite element method (FEM), the method of moment (MoM), etc. or the combination of these methods can be used to solve the problems. Finally, a numerical model is provided for the scattering problem with perfect matched layer (PML) by using FD. The comparison between the results of the proposed method and FDTD is presented to verify the proposed new method.展开更多
文摘The manuscript introduces an “ab initio” quantum model to deduce the Maxwell equations. After general considerations and laying out the model’s theoretical framework, these equations can be derived alongside a broad variety of other results. Specifically, a corollary of the present model proposes a possible mechanism underlying the formation of magnetic monopoles and allows estimating their formation energy in order of magnitude.
基金supported by Shandong Provincial Natural Science Foundation(Y2008A19)the Scientific Research Foundation for the Returned Overseas Chinese Scholars,State Education Ministry
文摘Several new energy identities of the two dimenslonal(2D) Maxwell equations in a lossy medium in the case of the perfectly electric conducting boundary conditions are proposed and proved. These identities show a new kind of energy conservation in the Maxwell system and provide a new energy method to analyze the alternating direction im- plicit finite difference time domain method for the 2D Maxwell equations (2D-ADI-FDTD). It is proved that 2D-ADI-FDTD is approximately energy conserved, unconditionally sta- ble and second order convergent in the discrete L2 and H1 norms, which implies that 2D-ADI-FDTD is super convergent. By this super convergence, it is simply proved that the error of the divergence of the solution of 2D-ADI-FDTD is second order accurate. It is also proved that the difference scheme of 2D-ADI-FDTD with respect to time t is second order convergent in the discrete H1 norm. Experimental results to confirm the theoretical analysis on stability, convergence and energy conservation are presented.
基金Project supported by the National Natural Science Foundation of China (Grants Nos 10672143 and 60575055)State Key Laboratory of Scientific and Engineering Computing,Chinese Academy of Sciences+1 种基金Tang Yi-Fa acknowledges the support under Sabbatical Program (SAB2006-0070) of the Spanish Ministry of Education and ScienceJimnez S and Vzquez L acknowledge support of the Spanish Ministry of Education and Science (Grant No MTM2005-05573)
文摘This paper presents a method to find Noether-type conserved quantities and Lie point symmetries for discrete mechanico-electrical dynamical systems,which leave invuriant the set of solutions of the corresponding difference scheme. This approach makes it possible to devise techniques for solving the Lagrange Maxwell equations in differences which correspond to mechanico-electrical systems,by adapting existing differential equations.In particular,it obtains a new systematic method to determine both the one-parameter Lie groups and the discrete Noether conserved quantities of Lie point symmetries for mechanico-electrical systems.As an application,it obtains the Lie point symmetries and the conserved quantities for the difference equation of a model that represents a capacitor microphone.
基金supported by NSFC(11341002)NSFC(11171104,10871066)+1 种基金the Construct Program of the Key Discipline in Hunansupported in part by US National Science Foundation under Grant DMS-1115530
文摘In this paper, a unified model for time-dependent Maxwell equations in dispersive media is considered. The space-time DG method developed in [29] is applied to solve the un-derlying problem. Unconditional L2-stability and error estimate of order O?τr+1+hk+1/2? are obtained when polynomials of degree at most r and k are used for the temporal dis-cretization and spatial discretization respectively. 2-D and 3-D numerical examples are given to validate the theoretical results. Moreover, numerical results show an ultra-convergence of order 2r+1 in temporal variable t.
基金Supported by NSFC (91430107/11771138/11171104)the Construct Program of the Key Discipline in Hunan+4 种基金partially supported by Scientific Research Fund of Hunan Provincial Education Department (19B325/19C1059)Hunan International Economics University (2017A05)supported by NSFC (11771137)the Construct Program of the Key Discipline in Hunan Provincea Scientific Research Fund of Hunan Provincial Education Department (16B154)。
文摘In this article,we investigate a stochastic Galerkin method for the Maxwell equations with random inputs.The generalized Polynomial Chaos(gPC)expansion technique is used to obtain a deterministic system of the gPC expansion coefficients.The regularity of the solution with respect to the random is analyzed.On the basis of the regularity results,the optimal convergence rate of the stochastic Galerkin approach for Maxwell equations with random inputs is proved.Numerical examples are presented to support the theoretical analysis.
基金supported by the Program of Excellent Team of Harbin Institute of Technology
文摘This paper is concerned with estimation of electrical conductivity in Maxwell equations. The primary difficulty lies in the presence of numerous local minima in the objective functional. A wavelet multiscale method is introduced and applied to the inversion of Maxwell equations. The inverse problem is decomposed into multiple scales with wavelet transform, and hence the original problem is reformulated to a set of sub-inverse problems corresponding to different scales, which can be solved successively according to the size of scale from the shortest to the longest. The stable and fast regularized Gauss-Newton method is applied to each scale. Numerical results show that the proposed method is effective, especially in terms of wide convergence, computational efficiency and precision.
文摘An FDTD system associated with uniaxial perfectly matched layer(UPML) for an electromagnetic scattering problem in two-dimensional space in polar coordinates is considered.Particularly the FDTD system of an initial-boundary value problems of the transverse magnetic(TM) mode to Maxwell's equations is obtained by Yee's algorithm,and the open domain of the scattering problem is truncated by a circle with a UPML.Besides,an artificial boundary condition is imposed on the outer boundary of the UPML.Afterwards,stability of the FDTD system on the truncated domain is established through energy estimates by the Gronwall inequality.Numerical experiments are designed to approve the theoretical analysis.
基金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.
文摘In the articles “Newtons Law of Universal Gravitation Explained by the Theory of Informatons” and “The Gravitational Interaction between Moving Mass Particles Explained by the Theory of Informatons” the gravitational interaction has been explained by the hypothesis that information carried by informatons is the substance of gravitational fields, i.e. the medium that the interaction in question makes possible. From the idea that “information carried by informatons” is its substance, it has been deduced that—on the macroscopic level—a gravitational field manifests itself as a dual entity, always having a field- and an induction component (Egand Bg) simultaneously created by their common sources. In this article we will mathematically deduce the Maxwell-Heaviside equations from the kinematics of the informatons. These relations describe on the macroscopic level how a gravitational field (Eg, Bg) is generated by whether or not moving masses and how spatial and temporal changes of Egand Bgare related. We show that there is no causal link between Egand Bg.
文摘The study of Electromagnetic Compatibility is essential to ensure the harmonious operation of electronic equipment in a shared environment. The basic principles of Electromagnetic Compatibility focus on the ability of devices to withstand electromagnetic disturbances and not produce disturbances that could affect other systems. Imperceptible in most work situations, electromagnetic fields can, beyond certain thresholds, have effects on human health. The objective of the present article is focused on the modeling analysis of the influence of geometric parameters of industrial static converters radiated electromagnetic fields using Maxwell’s equations. To do this we used the analytical formalism for calculating the electromagnetic field emitted by a filiform conductor, to model the electromagnetic radiation of this device in the spatio-temporal domain. The interactions of electromagnetic waves with human bodies are complex and depend on several factors linked to the characteristics of the incident wave. To model these interactions, we implemented the physical laws of electromagnetic wave propagation based on Maxwell’s and bio-heat equations to obtain consistent results. These obtained models allowed us to evaluate the spatial profile of induced current and temperature of biological tissue during exposure to electromagnetic waves generated by this system. The simulation 2D results obtained from computer tools show that the temperature variation and current induced by the electromagnetic field can have a very significant influence on the life of biological tissue. The paper provides a comprehensive analysis using advanced mathematical models to evaluate the influence of electromagnetic fields. The findings have direct implications for workplace safety, potentially influencing standards and regulations concerning electromagnetic exposure in industrial settings.
基金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.
文摘We applied a spatial high-order finite-difference-time-domain (HO-FDTD) scheme to solve 2D Maxwell’s equations in order to develop a fluid model employed to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma. We examined the performance of the applied scheme, in this context, we implemented the developed model to study selected phenomena in terahertz radiation production, such as the excitation energy and conversion efficiency of the produced THz radiation, in addition to the influence of the pulse chirping on properties of the produced radiation. The obtained numerical results have clarified that the applied HO-FDTD scheme is precisely accurate to solve Maxwell’s equations and sufficiently valid to study the production of terahertz radiation by the filamentation of two femtosecond lasers in air plasma.
文摘Maxwell’s equations in electromagnetism can be categorized into three dis-tinct groups based on the electromagnetic source when employing quaterni-ons. Each group represents a self-contained system in which Maxwell’s equations are applied and validated concurrently, in contrast to the previous approach that did not account for this. It has been noted that the formulation of these Maxwell equations ultimately results in the formulation of Max-well’s equations utilizing the scalar function.
基金This work is subsidized by the special funds for major state basic research projects (No. 1999032800).
文摘Superconvergence of the mixed finite element methods for 2-d Maxwell equations is studied in this paper. Two order of superconvergent factor can be obtained for the k-th Nedelec elements on the rectangular meshes.
基金supported by Shandong Provincial Natural Science Foundation (Grant No. Y2008A19)supported by Research Reward for Excellent Young Scientists from Shandong Province(Grant No. 2007BS01020) +1 种基金supported by Scientific Research Foundation for the Returned Overseas Chinese Scholars, State Education Ministrysupported by National Natural Science Foundation of China (Grant No. 11071244)
文摘This paper is concerned with the stability and superconvergence analysis of the famous finite-difference time-domain (FDTD) scheme for the 2D Maxwell equations in a lossy medium with a perfectly electric conducting (PEC) boundary condition, employing the energy method. To this end, we first establish some new energy identities for the 2D Maxwell equations in a lossy medium with a PEC boundary condition. Then by making use of these energy identities, it is proved that the FDTD scheme and its time difference scheme are stable in the discrete L2 and H1 norms when the CFL condition is satisfied. It is shown further that the solution to both the FDTD scheme and its time difference scheme is second-order convergent in both space and time in the discrete L2 and H1 norms under a slightly stricter condition than the CFL condition. This means that the solution to the FDTD scheme is superconvergent. Numerical results are also provided to confirm the theoretical analysis.
基金supported by the National Science Centre of Poland (Grant No. 2013/09/B/ST1/01963)
文摘We survey recent results on ground and bound state solutions E:?→R^3 of the problem {▽(▽×E)+}λE=|E|^(P-2)E in Ω,v×E=0 on Ω on a bounded Lipschitz domain ??R^3,where?×denotes the curl operator in R^3.The equation describes the propagation of the time-harmonic electric field R{E(χ)e^(iwt)}in a nonlinear isotropic material ? withλ=-μεω~2≤0,where μ andεstand for the permeability and the linear part of the permittivity of the material.The nonlinear term|E|^(P-2)E with 2<p≤2*=6 comes from the nonlinear polarization and the boundary conditions are those for?surrounded by a perfect conductor.The problem has a variational structure;however the energy functional associated with the problem is strongly indefinite and does not satisfy the Palais-Smale condition.We show the underlying difficulties of the problem and enlist some open questions.
基金supported by the National Natural Science Foundation ofChina (Grant Nos. 10535060,A050207,10975172 and 10821504)the National Basic Research Program of China (Grant No. 2007CB815401)
文摘Extending the holographic program of our previous work,we derive f(R) gravity and the Maxwell equations from the holographic principle,using time-like holographic screens.We find that to derive the Einstein equations and f(R) gravity by a natural holographic approach,the quasi-static condition is necessary.We also find the surface stress tensor and the surface electric current,surface magnetic current on a holographic screen for f(R) gravity and Maxwell's theory,respectively.
文摘The behavior of interacting electrons in a perfect crystal under macroscopic external electric and magnetic fields is studied. Effective Maxwell equations for the macroscopic electric and magnetic fields are derived starting from time-dependent density functional theory. Effective permittivity and permeability coefficients are obtained.
文摘A new multi-dimensional scheme for the Maxwell equations is established by the CIP method in combination with the method of characteristics(CIP-MOC).In addition,the CIP-MOC can be extended to arbitrary grid system by the Soroban grid without losing the third-order accuracy.With the accuracy fixed,the grid points required for the CIP are 40 times less than the conventional schemes like the FDTD in three dimensions.Numerical solutions obtained by the CIP-MOC are compared with analytical solution and the FDTD in plane-wave scattering by a perfectly-conducting circular cylinder,and the CIP-MOC agrees very well with analytical solutions.The Soroban grid is also applied to the Vlasov equation that describes the kinematics of plasmas that is frequently combined with the Maxwell equation.The adaptively moving points in velocity space are similar to the particle codes but can provide accurate solutions.
基金the National Natural Science Foundation of China (Grant Nos. 60701010 and 40774093)
文摘In this work, a stable numerical algorithm proposed by Chung et al. for the time-domain Maxwell equations is generalized. The time-domain Maxwell equations are solved by expressing the transient behaviors in terms of the modified Laguerre polynomials, and then the original equations of the initial value and boundary value can be transformed into a series of problems independent of the time variable. In this case the method of finite difference (FD), the finite element method (FEM), the method of moment (MoM), etc. or the combination of these methods can be used to solve the problems. Finally, a numerical model is provided for the scattering problem with perfect matched layer (PML) by using FD. The comparison between the results of the proposed method and FDTD is presented to verify the proposed new method.