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.

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.

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.

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.

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.

Derivation of a Revised Tsiolkovsky Rocket Equation That Predicts Combustion Oscillations

Our study identifies a subtle deviation from Newton’s third law in the derivation of the ideal rocket equation, also known as the Tsiolkovsky Rocket Equation (TRE). TRE can be derived using a 1D elastic collision model of the momentum exchange between the differential propellant mass element (dm) and the rocket final mass (m1), in which dm initially travels forward to collide with m1 and rebounds to exit through the exhaust nozzle with a velocity that is known as the effective exhaust velocity ve. We observe that such a model does not explain how dm was able to acquire its initial forward velocity without the support of a reactive mass traveling in the opposite direction. We show instead that the initial kinetic energy of dm is generated from dm itself by a process of self-combustion and expansion. In our ideal rocket with a single particle dm confined inside a hollow tube with one closed end, we show that the process of self-combustion and expansion of dm will result in a pair of differential particles each with a mass dm/2, and each traveling away from one another along the tube axis, from the center of combustion. These two identical particles represent the active and reactive sub-components of dm, co-generated in compliance with Newton’s third law of equal action and reaction. Building on this model, we derive a linear momentum ODE of the system, the solution of which yields what we call the Revised Tsiolkovsky Rocket Equation (RTRE). We show that RTRE has a mathematical form that is similar to TRE, with the exception of the effective exhaust velocity (ve) term. The ve term in TRE is replaced in RTRE by the average of two distinct exhaust velocities that we refer to as fast-jet, vx1, and slow-jet, vx2. These two velocities correspond, respectively, to the velocities of the detonation pressure wave that is vectored directly towards the exhaust nozzle, and the retonation wave that is initially vectored in the direction of rocket propagation, but subsequently becomes reflected from the thrust surface of the combustion chamber to exit through the exhaust nozzle with a time lag behind the detonation wave. The detonation-retonation phenomenon is supported by experimental evidence in the published literature. Finally, we use a convolution model to simulate the composite exhaust pressure wave, highlighting the frequency spectrum of the pressure perturbations that are generated by the mutual interference between the fast-jet and slow-jet components. Our analysis offers insights into the origin of combustion oscillations in rocket engines, with possible extensions beyond rocket engineering into other fields of combustion engineering.

ANALYSIS OF FDTD TO UPML FOR MAXWELL EQUATIONS IN POLAR COORDINATES

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.

Stochastic-periodic Homogenization of Maxwell＇s Equations with Linear and Periodic Conductivity

Stochastic-periodic homogenization is studied for the Maxwell equations with the linear ＂and periodic electric conductivity. It is shown by the stochastic-two-scale convergence method that the sequence of solutions to a class of highly oscillatory problems converges to the solution of a homogenized Maxwell equation.

Theoretical Maxwell's Equations, Gauge Field and Their Universality Based on One Conservation Law

The notion of the inner product of vectors is extended to tensors of different orders, which may replace the vector product usually. The essences of the differential and the codiffcrential forms are pointed out： they represent the tangent surface and the normal surface fluxes of a tensor, reslpetivcly. The definitions of the divergence and the curl of a 2D surface flux of a tensor arc obtained. Maxwell＇s equations, namely, the constraction law of field, which were usually established based on two conservation laws of electric charge and imaginary magnetic charge, are derived by the author only by using one conservation law （ mass or fluid flux quantity and so on） and the feature of central field （or its composition）. By the feature of central field （or its composition）, the curl of 2D flux is zero. Both universality of gauge field and the difficulty of magnetic monopole theory （ a magnetic monopole has no effect on electric current just like a couple hasing no effect on the sum of forces） axe presented： magnetic monopole has no the feature of magnet. Finally it is pointed out that the base of relation of mass and energy is already involved in Maxwell＇s equations.

Complex Maxwell's equations

A unified complex model of Maxwell＇s equations is presented.The wave nature of the electromagnetic field vector is related to the temporal and spatial distributions and the circulation of charge and current densities.A new vacuum solution is obtained,and a new transformation under which Maxwell＇s equations are invariant is proposed.This transformation extends ordinary gauge transformation to include charge-current as well as scalar-vector potential.An electric dipole moment is found to be related to the magnetic charges,and Dirac＇s quantization is found to determine an uncertainty relation expressing the indeterminacy of electric and magnetic charges.We generalize Maxwell＇s equations to include longitudinal waves.A formal analogy between this formulation and Dirac＇s equation is also discussed.

STABILITY AND SUPER CONVERGENCE ANALYSIS OF ADI-FDTD FOR THE 2D MAXWELL EQUATIONS IN A LOSSY MEDIUM

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.

Noether conserved quantities and Lie point symmetries of difference Lagrange-Maxwell equations and lattices

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.

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.

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.

Existence and Uniqueness of Global Smooth Solutions for Vlasov Maxwell Equations

Global existence of classical solutions to the relativistic Vlasov-Maxwell system, given sufficiently regular initial data, is a long-standing open problem. The aim of this project is to present in details the results of a paper published in 1986 by Robert Glassey and Walter Strauss. In that paper, a sufficient condition for the global existence of a smooth solution to the relativistic Vlasov-Maxwell system is derived. In the following, the resulting theorem is proved by taking initial data , . A small data global existence result is presented as well.

A Nonconforming Arbitrary Quadrilateral Finite Element Method for Approximating Maxwell's Equations

The main aim of this paper is to provide convergence analysis of Quasi-Wilson nonconforming finite element to Maxwell's equations under arbitrary quadrilateral meshes.The error estimates are derived,which are the same as those for conforming elements under conventional regular meshes.

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.

New Reflection Principles for Maxwell's Equations and Their Applications

Some new reflection principles for Maxwell's equations are first established, which are then applied to derive two novel identifiability results in inverse electromagnetic obstacle scattering problems with polyhedral scatterers.

PARTIAL COMPACTNESS FOR LANDAU-LIFSHITZ MAXWELL EQUATION IN TWO-DIMENSION

We study the partial regularity of weak solutions to the 2-dimensional Landau- Lifshitz equations coupled with time dependent Maxwell equations by Ginzburg-Landau type approximation. Outside an energy concentration set of locally finite 2-dimensional parabolic Hausdorff measure, we prove the uniform local C∞ bounds for the approaching solutions and then extract a subsequence converging to a global weak solution of the Landau-Lifshitz-Maxwell equations which are smooth away from finitely many points.

A Self-Consistent Model on Cylindrical Monopole Plasma Antenna Excited by Surface Wave Based on the Maxwell-Boltzmann Equation

The paper analyzes the motion of electron in plasma antenna and the distribution of electromagnetic field power around the plasma antenna, and proposes a self-consistent model according to the structure of cylindrical monopole plasma antenna excited by surface wave;calculation of the model is based on Maxwell-Boltzmann equation and gas molecular dynamics theory. The calculation results show that this method can reflect the relationships between the external excitation power, gas pressure, discharge current and the characteristic of plasma. It is an accurate method to predicate and calculate the parameters of plasma antenna.