Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userBifurcation Analysis of the Eigenstructure of the Discrete Single-curl Operator in Three-dimensional Maxwells Equations with Pasteur Media
327
0
0.0
(
0
)
Added by
Xin Liang
Publication date
2020
fields
Informatics Engineering
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
This paper focuses on studying the bifurcation analysis of the eigenstructure of the $gamma$-parameterized generalized eigenvalue problem ($gamma$-GEP) arising in three-dimensional (3D) source-free Maxwells equations with Pasteur media, where $gamma$ is the magnetoelectric chirality parameter. For the weakly coupled case, namely, $gamma < gamma_{*} equiv$ critical value, the $gamma$-GEP is positive definite, which has been well-studied by Chern et. al, 2015. For the strongly coupled case, namely, $gamma > gamma_{*}$, the $gamma$-GEP is no longer positive definite, introducing a totally different and complicated structure. For the critical strongly coupled case, numerical computations for electromagnetic fields have been presented by Huang et. al, 2018. In this paper, we build several theoretical results on the eigenstructure behavior of the $gamma$-GEPs. We prove that the $gamma$-GEP is regular for any $gamma > 0$, and the $gamma$-GEP has $2 times 2$ Jordan blocks of infinite eigenvalues at the critical value $gamma_{*}$. Then, we show that the $2 times 2$ Jordan block will split into a complex conjugate eigenvalue pair that rapidly goes down and up and then collides at some real point near the origin. Next, it will bifurcate into two real eigenvalues, with one moving toward the left and the other to the right along the real axis as $gamma$ increases. A newly formed state whose energy is smaller than the ground state can be created as $gamma$ is larger than the critical value. This stunning feature of the physical phenomenon would be very helpful in practical applications. Therefore, the purpose of this paper is to clarify the corresponding theoretical eigenstructure of 3D Maxwells equations with Pasteur media.
rate research
Read More
We describe a fourth-order accurate finite-difference time-domain scheme for solving dispersive Maxwells equations with nonlinear multi-level carrier kinetics models. The scheme is based on an efficient single-step three time-level modified equation approach for Maxwells equations in second-order form for the electric field coupled to ODEs for the polarization vectors and population densities of the atomic levels. The resulting scheme has a large CFL-one time-step. Curved interfaces between different materials are accurately treated with curvilinear grids and compatibility conditions. A novel hierarchical modified equation approach leads to an explicit scheme that does not require any nonlinear iterations. The hierarchical approach at interfaces leads to local updates at the interface with no coupling in the tangential directions. Complex geometry is treated with overset grids. Numerical stability is maintained using high-order upwind dissipation designed for Maxwells equations in second-order form. The scheme is carefully verified for a number of two and three-dimensional problems. The resulting numerical model with generalized dispersion and arbitrary nonlinear multi-level system can be used for many plasmonic applications such as for ab initio time domain modeling of nonlinear engineered materials for nanolasing applications, where nano-patterned plasmonic dispersive arrays are used to enhance otherwise weak nonlinearity in the active media.
For the Maxwells equations in a Havriliak-Negami (H-N) dispersive medium, the associated energy dissipation law has not been settled at both continuous level and discrete level. In this paper, we rigorously show that the energy of the H-N model can be bounded by the initial energy and the model is well-posed. We analyse a backward Euler-type semi-discrete scheme, and prove that the modified discrete energy decays monotonically in time. Such a strong stability ensures that the scheme is unconditionally stable. We also introduce a fast temporal convolution algorithm to alleviate the burden of the history dependence in the polarisation relation involving the singular kernel with the Mittag-Leffler function with three parameters. We provide ample numerical results to demonstrate the efficiency and accuracy of a full-discrete scheme via a spectra-Galerkin method in two dimensions. Finally, we consider an interesting application in the recovery of complex relative permittivity and some related physical quantities.
We propose a controllability method for the numerical solution of time-harmonic Maxwells equations in their first-order formulation. By minimizing a quadratic cost functional, which measures the deviation from periodicity, the controllability method determines iteratively a periodic solution in the time domain. At each conjugate gradient iteration, the gradient of the cost functional is simply computed by running any time-dependent simulation code forward and backward for one period, thus leading to a non-intrusive implementation easily integrated into existing software. Moreover, the proposed algorithm automatically inherits the parallelism, scalability, and low memory footprint of the underlying time-domain solver. Since the time-periodic solution obtained by minimization is not necessarily unique, we apply a cheap post-processing filtering procedure which recovers the time-harmonic solution from any minimizer. Finally, we present a series of numerical examples which show that our algorithm greatly speeds up the convergence towards the desired time-harmonic solution when compared to simply running the time-marching code until the time-harmonic regime is eventually reached.
In this paper, we consider Maxwells equations in linear dispersive media described by a single-pole Lorentz model for electronic polarization. We study two classes of commonly used spatial discretizations: finite difference methods (FD) with arbitrary even order accuracy in space and high spatial order discontinuous Galerkin (DG) finite element methods. Both types of spatial discretizations are coupled with second order semi-implicit leap-frog and implicit trapezoidal temporal schemes studied in our previous research [5,6]. By performing detailed dispersion analysis for the semi-discrete and fully discrete schemes, we obtain rigorous quantification of the dispersion error for Lorentz dispersive dielectrics. In particular, comparisons of dispersion error can be made taking into account the model parameters, and mesh sizes in the design of the two types of schemes. The results for the numerical dispersion analysis can guide us in the optimal choice of discretization parameters for the more complicated and nonlinear models. The numerical dispersion analysis of the fully discrete FD and DG schemes, for the dispersive Maxwell model considered in this paper, clearly indicate the dependence of the numerical dispersion errors on spatial and temporal discretizations, their order of accuracy, mesh discretization parameters and model parameters. The results obtained here cannot be arrived at by considering discretizations of Maxwells equations in free space. In particular, our results contrast the advantages and disadvantages of using high order FD or DG schemes and leap-frog or trapezoidal time integrators over different frequency ranges using a variety of measures of numerical dispersion errors. Finally, we highlight the limitations of the second order accurate temporal discretizations considered.
We analyse a PDE system modelling poromechanical processes (formulated in mixed form using the solid deformation, fluid pressure, and total pressure) interacting with diffusing and reacting solutes in the medium. We investigate the well-posedness of the nonlinear set of equations using fixed-point theory, Fredholms alternative, a priori estimates, and compactness arguments. We also propose a mixed finite element method and rigorously demonstrate the stability of the scheme. Error estimates are derived in suitable norms, and numerical experiments are conducted to illustrate the mechano-chemical coupling and to verify the theoretical rates of convergence.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
