Energy Stable Nodal DG Methods for Maxwell’s Equations of Mixed-Order Form in Nonlinear Optical Media

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.

Analysis of Modeling the Influence of Electromagnetic Fields Radiated by Industrial Static Converters and Impacts on Operators Using Maxwell’s Equations

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.

Complex Maxwell’s Equations

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.

Adaptive Maxwell’s Equations Derived Optimization and Its Application in Antenna Array Synthesis

In this paper,a self-adaptive method for the Maxwell’s Equations Derived Optimization(MEDO)is proposed.It is implemented by applying the Sequential Model-Based Optimization(SMBO)algorithm to the iterations of the MEDO,and achieves the automatic adjustment of the parameters.The proposed method is named as adaptive Maxwell’s equations derived optimization(AMEDO).In order to evaluate the performance of AMEDO,eight benchmarks are used and the results are compared with the original MEDO method.The results show that AMEDO can greatly reduce the workload of manual adjustment of parameters,and at the same time can keep the accuracy and stability.Moreover,the convergence of the optimization can be accelerated due to the dynamical adjustment of the parameters.In the end,the proposed AMEDO is applied to the side lobe level suppression and array failure correction of a linear antenna array,and shows great potential in antenna array synthesis.

Determination of Photon and Elementary Particles Rest Masses Using Maxwell’s Equations and Generalized Potential Dependent Special Relativity

The nature and origin of the photon and elementary rest masses are some of the challeng-ing problems that physics face. The approaches used to solve these problems are complex and time-consuming. Specifically, the photon rest mass pays attention to theoretical physi-cists. Many experimental works show that the photon rest mass is non zero. This problem can be solved using generalized potential dependent special relativity, which has been de-rived using simple arguments, and Maxwell’s equations, besides the conventional Einstein energy-momentum relation. The results obtained show that the rest mass of photons and elementary particles are strongly dependent on the vacuum energy and a universal con-stant. This result conforms with the models that predict time decaying vacuum energy as-sociated with production of smaller rest mass particles followed by larger masses. The two potential dependent mass expressions conform with the cosmological models that suggest the photon is generated first by assuming the universe consisting of total constant vacuum with decaying cosmological part and mass generating part. Using Maxwell’s equations, beside plank and De Broglie hypothesis together with special relativity energy-momentum relation the photon rest mass is estimated. It was shown that the photon rest mass is ex-tremely small compared to the electron mass.

Analysis of Maxwell’s Equations. Cosmic Magnetism

According to Hypersphere World-Universe Model, dark matter particles DIRACs are magnetic dipoles consisting of two Dirac’s monopoles. We conclude that DIRACs are the subject of Maxwell’s equations. So-called “auxiliary” magnetic field intensity H is indeed current density of magnetic dipoles. The developed approach to magnetic field can explain a wealth of discovered phenomena in Cosmic Magnetism: a dark magnetic field, the large-scale structure of the Milky Way’s magnetic field, and other magnetic phenomena which are only partly related to objects visible in other spectral ranges.

A New Formulation of Maxwell’s Equations in Clifford Algebra

A new unification of the Maxwell equations is given in the domain of Clifford algebras with in a fashion similar to those obtained with Pauli and Dirac algebras. It is shown that the new electromagnetic field multivector can be obtained from a potential function that is closely related to the scalar and the vector potentials of classical electromagnetics. Additionally it is shown that the gauge transformations of the new multivector and its potential function and the Lagrangian density of the electromagnetic field are in agreement with the transformation rules of the second-rank antisymmetric electromagnetic field tensor, in contrast to the results obtained by applying other versions of Clifford algebras.

Solution of Maxwell’s Equations for Cylindrical Symmetry Waveguides

The solution of Maxwell’s equations for a piecewise homogeneous medium of cylindrical symmetry has been obtained. The parameters of the cylindrical waveguide modes have been calculated on its basis. The conclusions are confirmed by numerical calculation of the first four modes of a hollow metal waveguide operating as a mode convector.

Derivation of Maxwell’s Equations via the Covariance Requirements of the Special Theory of Relativity, Starting with Newton’s Laws

The purpose of this paper is to establish a connection between Maxwell’s equations, Newton’s laws, and the special theory of relativity. This is done with a derivation that begins with Newton’s verbal enunciation of his first two laws. Derived equations are required to be covariant, and a simplicity criterion requires that the four-vector force on a charged particle be linearly related to the four-vector velocity. The connecting tensor has derivable symmetry properties and contains the electric and magnetic field vectors. The Lorentz force law emerges, and Maxwell’s equations for free space emerge with the assumption that the tensor and its dual must both satisfy first-order partial differential equations. The inhomogeneous extension yields a charge density and a current density as being the source of the field, and yields the law of conservation of charge. Newton’s third law is reinterpreted as a reciprocity statement, which requires that the charge in the source term can be taken as the same physical entity as that of the test particle and that both can be assigned the same units. Requiring covariance under either spatial inversions or time reversals precludes magnetic charge being a source of electromagnetic fields that exert forces on electric charges.

A Fresh View for Maxwell’s Equations and Electromagnetic Wave Propagation

Equations related with wave propagation are reexamined as in certain circumstances law of conservation of energy is not fulfilled even though it is cautiously clarified with the help of Heisenberg’s uncertainty principle. Recently, attempt has also been made to understand certain discrepancies in optical phenomena like diffraction or interference. The purpose of the present investigation, therefore, is to overcome some discrepancies by introducing constants of integration in Maxwell’s Equation. It turns out that the presence of vibrating strings (or store energy) in the medium becomes essential to understand several details of the wave propagation.

Maxwell’s Equations as the Basis for Model of Atoms

A century ago the classical physics couldn’t explain many atomic physical phenomena. Now the situation has changed. It’s because within the framework of classical physics with the help of Maxwell’s equations we can derive Schrödinger’s equation, which is the foundation of quantum physics. The equations for energy, momentum, frequency and wavelength of the electromagnetic wave in the atom are derived using the model of atom by analogy with the transmission line. The action constant A0 = (μ0/ε0)1/2s02e2 is a key term in the above mentioned equations. Besides the other well-known constants, the only unknown constant in the last expression is a structural constant of the atom s0. We have found that the value of this constant is 8.277 56 and that it shows up as a link between macroscopic and atomic world. After calculating this constant we get the theory of atoms based on Maxwell’s and Lorentz equations only. This theory does not require knowledge of Planck’s constant h, which is replaced with theoretically derived action constant A0, while the replacement for the fine structure constant α-1 is theoretically derived expression 2s02 = 137.036. So, the structural constant s0 replaces both constants h and α. This paper also defines the stationary states of atoms and shows that the maximal atomic number is equal to Zmax = 137. The presented model of the atoms covers three of the four fundamental interactions, namely the electromagnetic, weak and strong interactions.

High-Order Spatial FDTD Solver of Maxwell’s Equations for Terahertz Radiation Production

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.

The Maxwell-Heaviside Equations Explained by the Theory of Informatons

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.

Maxwell Equations and Magnetic Monopoles

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.

Recent Progress on the Maxwell’s Equations for Describing a Mechano-Driven Medium System with Multiple Moving Objects/Media

Maxwell’s equations for a mechano-driven media system(MEs-f-MDMS)have been used to characterize the electromagnetism of multislow-moving media that may be accelerated with complex trajectories.Such an approach starts from the integral forms of the four physics laws and is different from the classical approach of using the Lorentz transformation for correlating the electromagnetic phenomena observed in two inertial reference frames with relative motion.The governing equations inside the moving object/medium are the MEs-f-MDMS,and those in vacuum are the classical Maxwell’s equations;the full solutions of both reconcile at the medium surface/interface and satisfy the boundary conditions.This paper reviews the background,physical principle,and mathematical derivations for formulating the MEs-f-MDMS.Strategies are also presented for mathematically solving the MEs-f-MDMS.The unique advances made by the MEs-f-MDMS have been systematically summarized,as are their potential applications in engineering.We found that the Lorentz transformation is perfect for treating the electromagnetic phenomena of moving point charges in vacuum;however,for moving objects,the covariance of Maxwell’s equations may not hold,and use of the MEs-f-MDMS may be required if the velocity is low.Finally,recent advances for treating the boundary conditions at the nanoscale without assuming an abrupt boundary are also reviewed.

ON DISCRETE ENERGY DISSIPATION OF MAXWELL’S EQUATIONS IN A COLE-COLE DISPERSIVE MEDIUM

A simple criterion is studied for the first time for identifying the discrete energy dissipation of the Crank-Nicolson scheme for Maxwell’s equations in a Cole-Cole dispersive medium.Several numerical formulas that approximate the time fractional derivatives are investigated based on this criterion,including the L1 formula,the fractional BDF-2,and the shifted fractional trapezoidal rule(SFTR).Detailed error analysis is provided within the framework of time domain mixed finite element methods for smooth solutions.The convergence results and discrete energy dissipation law are confirmed by numerical tests.For nonsmooth solutions,the method SFTR can still maintain the optimal convergence order at a positive time on uniform meshes.Authors believe this is the first appearance that a second-order time-stepping method can restore the optimal convergence rate for Maxwell’s equations in a Cole-Cole dispersive medium regardless of the initial singularity of the solution.

Euler’s First-Order Explicit Method–Peridynamic Differential Operator for Solving Population Balance Equations of the Crystallization Process

Using Euler’s first-order explicit(EE)method and the peridynamic differential operator(PDDO)to discretize the time and internal crystal-size derivatives,respectively,the Euler’s first-order explicit method–peridynamic differential operator(EE–PDDO)was obtained for solving the one-dimensional population balance equation in crystallization.Four different conditions during crystallization were studied:size-independent growth,sizedependent growth in a batch process,nucleation and size-independent growth,and nucleation and size-dependent growth in a continuous process.The high accuracy of the EE–PDDO method was confirmed by comparing it with the numerical results obtained using the second-order upwind and HR-van methods.The method is characterized by non-oscillation and high accuracy,especially in the discontinuous and sharp crystal size distribution.The stability of the EE–PDDO method,choice of weight function in the PDDO method,and optimal time step are also discussed.

Implicit DG Method for Time Domain Maxwell’s Equations Involving Metamaterials

An implicit discontinuous Galerkin method is introduced to solve the timedomain Maxwell’s equations in metamaterials.The Maxwell’s equations in metamaterials are represented by integral-differential equations.Our scheme is based on discontinuous Galerkin method in spatial domain and Crank-Nicolson method in temporal domain.The fully discrete numerical scheme is proved to be unconditionally stable.When polynomial of degree at most p is used for spatial approximation,our scheme is verified to converge at a rate of O(τ^(2)+h^(p)+1/2).Numerical results in both 2D and 3D are provided to validate our theoretical prediction.

Adaptive hp-FEM with Arbitrary-Level Hanging Nodes for Maxwell’s Equations

Adaptive higher-order finite element methods(hp-FEM)are well known for their potential of exceptionally fast(exponential)convergence.However,most hp-FEM codes remain in an academic setting due to an extreme algorithmic complexity of hp-adaptivity algorithms.This paper aims at simplifying hpadaptivity for H(curl)-conforming approximations by presenting a novel technique of arbitrary-level hanging nodes.The technique is described and it is demonstrated numerically that it makes adaptive hp-FEM more efficient compared to hp-FEM on regular meshes and meshes with one-level hanging nodes.

A High-OrderDiscontinuous GalerkinMethod for the Two-Dimensional Time-Domain Maxwell’s Equations on Curved Mesh

In this paper,a DG(Discontinuous Galerkin)method which has been widely employed in CFD(Computational Fluid Dynamics)is used to solve the twodimensional time-domain Maxwell’s equations for complex geometries on unstructured mesh.The element interfaces on solid boundary are treated in both curved way and straight way.Numerical tests are performed for both benchmark problems and complex cases with varying orders on a series of grids,where the high-order convergence in accuracy can be observed.Both the curved and the straight solid boundary implementation can give accurate RCS(Radar Cross-Section)results with sufficiently small mesh size,but the curved solid boundary implementation can significantly improve the accuracy when using relatively large mesh size.More importantly,this CFD-based high-order DG method for the Maxwell’s equations is very suitable for complex geometries.