Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA More Accurate, Stable, FDTD Algorithm for Electromagnetics in Anisotropic Dielectrics
150
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
A more accurate, stable, finite-difference time-domain (FDTD) algorithm is developed for simulating Maxwells equations with isotropic or anisotropic dielectric materials. This algorithm is in many cases more accurate than previous algorithms (G. R. Werner et. al., 2007; A. F. Oskooi et. al., 2009), and it remedies a defect that causes instability with high dielectric contrast (usually for epsilon{} significantly greater than 10) with either isotropic or anisotropic dielectrics. Ultimately this algorithm has first-order error (in the grid cell size) when the dielectric boundaries are sharp, due to field discontinuities at the dielectric interface. Accurate treatment of the discontinuities, in the limit of infinite wavelength, leads to an asymmetric, unstable update (C. A. Bauer et. al., 2011), but the symmetrized version of the latter is stable and more accurate than other FDTD methods. The convergence of field values supports the hypothesis that global first-order error can be achieved by second-order error in bulk material with zero-order error on the surface. This latter point is extremely important for any applications measuring surface fields.
rate research
Read More
We have developed a new fully anisotropic 3D FDTD Maxwell solver for arbitrary electrically and magnetically anisotropic media for piecewise constant electric and magnetic materials that are co-located over the primary computational cells. Two numerical methods were developed that are called non-averaged and averaged methods, respectively. The non-averaged method is first order accurate, while the averaged method is second order accurate for smoothly-varying materials and reduces to first order for discontinuous material distributions. For the standard FDTD field locations with the co-location of the electric and magnetic materials at the primary computational cells, the averaged method require development of the different inversion algorithms of the constitutive relations for the electric and magnetic fields. We provide a mathematically rigorous stability proof followed by extensive numerical testing that includes long-time integration, eigenvalue analysis, tests with extreme, randomly placed material parameters, and various boundary conditions. For accuracy evaluation we have constructed a test case with an explicit analytic solution. Using transformation optics, we have constructed complex, spatially inhomogeneous geometrical object with fully anisotropic materials and a large dynamic range of $underline{epsilon}$ and $underline{mu}$, such that a plane wave incident on the object is perfectly reconstructed downstream. In our implementation, the considerable increase in accuracy of the averaged method only increases the computational run time by 20%.
In a recent paper, Klaseboer et al. (IEEE Trans. Antennas Propag., vol. 65, no. 2, pp. 972-977, Feb. 2017) developed a surface integral formulation of electromagnetics that does not require working with integral equations that have singular kernels. Instead of solving for the induced surface currents, the method involves surface integral solutions for 4 coupled Helmholtz equations: 3 for each Cartesian component of the electric E field plus 1 for the scalar function r*E on the surface of scatterers. Here we improve on this approach by advancing a formulation due to Yuffa et al. (IEEE Trans.Antennas Propag., vol. 66, no. 10, pp. 5274-5281, Oct. 2018) that solves for E and its normal derivative. Apart from a 25% reduction in problem size, the normal derivative of the field is often of interest in micro-photonic applications.
Following the recent development of a stable event-detection algorithm for hard-sphere systems, the implications of more complex interaction models are examined. The relative location of particles leads to ambiguity when it is used to determine the interaction state of a particle in stepped potentials, such as the square-well model. To correctly predict the next event in these systems, the concept of an additional state that is tracked separately from the particle position is introduced and integrated into the stable algorithm for event detection.
We compute the stable homotopy groups up to dimension 90, except for some carefully enumerated uncertainties.
We present a more accurate numerical scheme for the calculation of diffusive shock acceleration of cosmic rays using Stochastic Differential Equations. The accuracy of this scheme is demonstrated using a simple analytical flow profile that contains a shock of finite width and a varying diffusivity of the cosmic rays, where the diffusivity decreases across the shock. We compare the results for the slope of the momentum distribution with those obtained from a perturbation analysis valid for finite but small shock width. These calculations show that this scheme, although computationally more expensive, provides a significantly better performance than the Cauchy-Euler type schemes that were proposed earlier in the case where steep gradients in the cosmic ray diffusivity occur. For constant diffusivity the proposed scheme gives similar results as the Cauchy-Euler scheme.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
