An explicit weighted essentially nonoscillatory time-domain algorithm for 3-D EMC applications with arbitrary media interfaces

Abstract

Abstract-The precise modeling of arbitrarily aligned different material interfaces in composite electromagnetic compatibility (EMC) applications, is introduced in this paper through a fully explicit family of three-dimensional weighted essentially nonoscillatory schemes in the time domain. Developing curvilinear tessellations of auxiliary nodes compliant with a convex combination of all optimal candidate stencils, the novel algorithm allocates pertinent weights for wide-band evaluations at the particular regions. Also, via an enhanced temporal integration, continuity conditions are imposed in a local regime, without any unfavorable instabilities or vector parasites. Thus, the prior volumetric concepts suppress intrinsic grid defects even when time-steps are fairly large unlike in usual approaches. Numerical analysis verifies these benefits by diverse realistic EMC structures with laborious interfaces of complex shapes, curvatures, and dissimilar media configurations.

Combined FVTD/PSTD Schemes with Enhanced Spectral Accuracy for the Design of Large-Scale EMC Applications

“…Real parameters a l (l = 1,2,3) assure the stability of the schemes and subdue late-time exponential oscillations during the handling of discontinuities [30], by satisfying Kth-order constraint…”

“…To exploit all fluxes through the cell faces, the concept of the cell average -i.e. the average of G placed at the barycenter of every cell [30] -is described by…”

“…Next, the convex combination of all stencils (31), as dictated by these nodes, is generated and the process is ready to be employed. It is important to stress that the scheme is implemented independently near any discontinuity, without the need to use the FVTD/PSTD hybrid technique [21,30]. Thus, instead of applying (12) or (26) throughout such critical details, the piecewise polynomials of (30) along with the optimized stencils are utilized.…”

“…The most convenient way to built our conformal lattice in a general coordinate system, properly selected to guarantee cells of optimal quality without face skewing or deformations, is the use of a non-overlapping dual-grid discretization [28][29][30][31]. Actually, the precision of the proposed modeling algorithm relies on the systematic construction of the mesh in order to obtain the correct curves which enclose the object, under study, and the proper surfaces that surround its volume.…”