pith. sign in

arxiv: 2604.14618 · v1 · submitted 2026-04-16 · 💻 cs.CE

A Stable SBP-SAT FDTD Subgridding Method Without Region Split

Yuhui Wang , Langran Deng , Weibo Wu , Hanhong Liu , Xinyue Zhang , Xingqi Zhang , Jian Wang , Wei-Jie Wang
show 2 more authors
Zhizhang Chen Shunchuan Yang
This is my paper

Pith reviewed 2026-05-10 09:32 UTC · model grok-4.3

classification 💻 cs.CE
keywords SBP-SATFDTD subgriddingprojection operatorsdiscrete energy stabilityembedded grid interfacesfinite-difference time-domaincomputational electromagnetics
0
0 comments X

The pith

Projection SBP operators enable stable FDTD subgridding with direct coarse-fine coupling and no domain split.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a summation-by-parts simultaneous approximation term finite-difference time-domain subgridding scheme that places a refined patch inside a single coarse grid without auxiliary blocks or fragmentation. Tailored projection operators handle the embedded interface while derived SAT terms enforce boundary coupling. Discrete energy analysis then proves that the total energy remains bounded for arbitrarily long times, removing the need for aligned multi-block decompositions used in prior SBP-SAT work. This yields fewer interface conditions, lower computational cost, and higher accuracy near the grid transition.

Core claim

By designing projection SBP operators tailored for embedded topological features and deriving the corresponding SAT boundary conditions, this approach guarantees long-time stability through discrete energy analysis. The method enables direct coupling between an internal refined region and a single surrounding coarse-grid domain without introducing auxiliary blocks or causing domain fragmentation.

What carries the argument

Projection SBP operators constructed for embedded topological features, together with the SAT boundary conditions they induce, which together preserve the summation-by-parts property and deliver an energy-stable interface.

If this is right

  • Direct coupling of one refined interior region to one outer coarse domain eliminates the extra SAT surfaces required by multi-block decompositions.
  • Long-time stability follows from a single discrete energy estimate rather than from grid alignment assumptions.
  • Fewer interface conditions and no domain fragmentation reduce overall computational complexity.
  • Accuracy near the grid interface improves because the projection operators are designed specifically for the embedded topology.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same operator-construction strategy may apply to other hyperbolic systems that admit SBP discretizations.
  • Dynamic insertion or removal of refined patches could become feasible if the projection operators can be updated locally.
  • Large-scale electromagnetic problems with irregular material boundaries might adopt this subgridding pattern to avoid global grid refinement.

Load-bearing premise

Projection SBP operators can be built for any embedded topological feature while still satisfying the summation-by-parts identity and admitting stable SAT coupling.

What would settle it

A long-time simulation on an arbitrarily shaped embedded refinement that exhibits growing energy or instability would disprove the stability guarantee.

Figures

Figures reproduced from arXiv: 2604.14618 by Hanhong Liu, Jian Wang, Langran Deng, Shunchuan Yang, Weibo Wu, Wei-Jie Wang, Xingqi Zhang, Xinyue Zhang, Yuhui Wang, Zhizhang Chen.

Figure 1
Figure 1. Figure 1: Comparison of domain decomposition strategies and proposed topology without region split. (a) Aligned-block decomposition resulting in excessive [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Definition of the 1-D boundaries and indicator matrices for the outer [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The recorded time-domain electric field Ez at the observation point for 106 time steps, demonstrating the long-time stability of the proposed method with grid ratios of (a) 1:5 and (b) 2:3. B. Numerical Reflection from the Subgridding Interface To quantitatively evaluate the numerical reflection induced by the geometric transition and interpolation matrices at the subgridding interfaces, a waveguide model … view at source ↗
Figure 4
Figure 4. Figure 4: Schematic of the waveguide model for the numerical reflection test, [PITH_FULL_IMAGE:figures/full_fig_p007_4.png] view at source ↗
Figure 7
Figure 7. Figure 7: Normalized magnitude spectra of the recorded electric field under [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗
Figure 6
Figure 6. Figure 6: Geometric configuration of the two separated [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 9
Figure 9. Figure 9: Spatial distribution of the relative error for [PITH_FULL_IMAGE:figures/full_fig_p009_9.png] view at source ↗
Figure 11
Figure 11. Figure 11: Geometric configuration of the computational domain for the human [PITH_FULL_IMAGE:figures/full_fig_p010_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Multi-scale visualization of the spatial discretization within the human head model. (a, f) Global distributions under the 1:10 mesh resolution, with [PITH_FULL_IMAGE:figures/full_fig_p011_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Time-domain electric field Ez recorded at the observation probe. (a) Global waveforms. (b) A localized zoomed-in view showing the numerical convergence of the proposed method at different refinement ratios. penalty errors. The quantitative performance metrics are summarized in Table II. At the 1:2 ratio, the proposed method eliminates 8 and 4 redundant interfaces compared to the aligned and T-junction met… view at source ↗
read the original abstract

A provably stable summation-by-parts simultaneous approximation term (SBP-SAT) finite-difference time-domain (FDTD) subgridding method without region split is proposed. By designing projection SBP operators tailored for embedded topological features and deriving the corresponding SAT boundary conditions, this approach guarantees long-time stability through discrete energy analysis. Unlike conventional SBP-SAT FDTD subgridding techniques that rely on aligned or multi-block configurations, the proposed method enables a direct coupling between an internal refined region and a single surrounding coarse-grid domain without introducing auxiliary blocks or causing domain fragmentation. Numerical results validate the efficiency, accuracy, and topological flexibility of the proposed method. Compared with existing multi-block SBP-SAT methods, this method effectively reduces computational complexity by minimizing SAT boundary conditions and improves calculation accuracy near grid interfaces.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript proposes a summation-by-parts simultaneous approximation term (SBP-SAT) finite-difference time-domain (FDTD) subgridding method that enables direct coupling of an internal refined region to a single surrounding coarse-grid domain without region splitting or auxiliary blocks. It designs projection SBP operators tailored to embedded topological features, derives corresponding SAT boundary conditions, and claims to prove long-time stability via discrete energy analysis. Numerical experiments are presented to demonstrate efficiency, accuracy, and flexibility compared to conventional multi-block SBP-SAT approaches.

Significance. If the stability result holds without hidden alignment assumptions, the method would meaningfully simplify subgridding in FDTD simulations of complex geometries, reducing the number of SAT interfaces and computational overhead while maintaining accuracy near grid transitions. This addresses a practical limitation in existing multi-block techniques and could enable more flexible mesh refinement in computational electromagnetics.

major comments (1)
  1. [§3.2 and §4] §3.2 (Projection Operator Construction) and §4 (Discrete Energy Analysis): The central stability claim rests on the projection SBP operators preserving the exact summation-by-parts identity when coupling grids of differing spacing and arbitrary embedded topology. The energy estimate in §4 cancels interface terms only if the projection is norm-consistent and satisfies a discrete integration-by-parts relation at the interface; for non-rectangular or non-aligned topologies this is not automatic. Please supply the explicit verification that the telescoping property holds (e.g., the inner-product identity after projection) or state the precise topological restrictions under which it is guaranteed.
minor comments (2)
  1. [Abstract] Abstract: The phrase 'numerical results validate the efficiency, accuracy, and topological flexibility' is vague; a single sentence listing the specific test geometries and error norms used would improve readability.
  2. [§5] §5 (Numerical Results): The comparison tables would benefit from an additional column reporting the L2 error against a reference uniform fine-grid solution, rather than only against the multi-block baseline.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and constructive feedback on our manuscript. The major comment is addressed in detail below. We believe the requested clarifications can be incorporated without altering the core contributions.

read point-by-point responses

  1. Referee: [§3.2 and §4] §3.2 (Projection Operator Construction) and §4 (Discrete Energy Analysis): The central stability claim rests on the projection SBP operators preserving the exact summation-by-parts identity when coupling grids of differing spacing and arbitrary embedded topology. The energy estimate in §4 cancels interface terms only if the projection is norm-consistent and satisfies a discrete integration-by-parts relation at the interface; for non-rectangular or non-aligned topologies this is not automatic. Please supply the explicit verification that the telescoping property holds (e.g., the inner-product identity after projection) or state the precise topological restrictions under which it is guaranteed.

    Authors: We agree that an explicit verification strengthens the presentation. In §3.2 the projection SBP operators are constructed precisely so that the discrete integration-by-parts identity is preserved after projection onto the embedded interface; this is achieved by enforcing both the SBP property on each grid and norm-consistency of the projection operator at the coarse-fine transition. Section 4 then uses this identity to cancel the interface contributions in the energy estimate, yielding unconditional stability. The construction does not impose rectangular or alignment restrictions beyond the topological embedding itself. To make the telescoping property fully transparent, we will add a short lemma (with the inner-product identity after projection) in the revised §3.2 or as an appendix, together with a statement of the precise topological class (simply-connected embedded regions with piecewise-smooth boundaries) for which the operators are defined. revision: yes

Circularity Check

0 steps flagged

No circularity: stability follows from standard energy analysis on explicitly constructed operators

full rationale

The derivation chain consists of (1) constructing projection SBP operators tailored to embedded topologies, (2) deriving compatible SAT terms, and (3) applying the standard discrete energy method to obtain a stability bound. None of these steps reduces to a tautology, a fitted parameter renamed as prediction, or a self-citation whose content is the target result. The energy analysis is an independent, well-established technique whose telescoping property is verified once the operators satisfy the SBP identity by design; the construction itself is not claimed to be forced by prior self-referential theorems. Numerical experiments supply external validation outside the analytic steps. No load-bearing self-citation or ansatz smuggling is present in the abstract or described chain.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the existence and constructibility of projection SBP operators for embedded features plus the applicability of discrete energy analysis; no free parameters or invented entities are mentioned.

axioms (1)
  • domain assumption Discrete energy analysis suffices to prove long-time stability once the SBP property and SAT terms are satisfied.
    Invoked to guarantee stability in the abstract.

pith-pipeline@v0.9.0 · 5463 in / 1289 out tokens · 43048 ms · 2026-05-10T09:32:14.087996+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

38 extracted references · 38 canonical work pages

  1. [1]

    Numerical solution of initial boundary value problems in- volving maxwell’s equations in isotropic media,

    K. Yee, “Numerical solution of initial boundary value problems in- volving maxwell’s equations in isotropic media,”IEEE Trans. Antennas Propag., vol. 14, no. 3, pp. 302-307, May 1966

    work page 1966

  2. [2]

    Microstrip antennas integrated with electromagnetic band-gap (EBG) structures: a low mutual coupling design for array applications,

    F. Yang and Y . Rahmat-Samii, “Microstrip antennas integrated with electromagnetic band-gap (EBG) structures: a low mutual coupling design for array applications,”IEEE Trans. Antennas Propag., vol. 51, no. 10, pp. 2936-2946, Oct. 2003

    work page 2003

  3. [3]

    Research on Electromagnetic Reception of Acoustically Excited Antenna by FDTD Method,

    W. Li, J. Song, Y . Fu, D. Li and X. Zhu, “Research on Electromagnetic Reception of Acoustically Excited Antenna by FDTD Method,”IEEE Antennas Wireless Propag. Lett., vol. 23, no. 1, pp. 274-278, Jan. 2024

    work page 2024

  4. [4]

    The finite-difference time-domain (FDTD) and the finite-volume time-domain (FVTD) methods in solving Maxwell’s equations,

    K. Yee and J. Chen, “The finite-difference time-domain (FDTD) and the finite-volume time-domain (FVTD) methods in solving Maxwell’s equations,”IEEE Trans. Antennas Propag., vol. 45, no. 3, pp. 354-363, Mar. 1997

    work page 1997

  5. [5]

    Backscattering Cross Section of Ultrawideband Antennas,

    S. Hu, H. Chen, Choi Look Law, Z. Shen, L. Zhu, W. Zhang and W. Dou, “Backscattering Cross Section of Ultrawideband Antennas,”IEEE Antennas Wireless Propag. Lett., vol. 6, pp. 70-73, Mar. 2007

    work page 2007

  6. [6]

    A subgridding method for the time-domain finite-difference method to solve Maxwell’s equations,

    S. S. Zivanovic, K. S. Yee, and K. K. Mei, “A subgridding method for the time-domain finite-difference method to solve Maxwell’s equations,” IEEE Trans. Microw. Theory Techn., vol. 39, no. 3, pp. 471-479, Mar. 1991

    work page 1991

  7. [7]

    A consistent subgridding scheme for the finite difference time domain method,

    P. Thoma and T. Weiland, “A consistent subgridding scheme for the finite difference time domain method,”Int. J. Numer. Model. El., vol. 9, no. 5, pp. 359-374, 1996

    work page 1996

  8. [8]

    The courant–friedrichs–lewy (cfl) condition,

    De Moura, Carlos A., and Carlos S. Kubrusly, “The courant–friedrichs–lewy (cfl) condition,” AMC 10.12 (2013): 45- 90

    work page 2013

  9. [9]

    A FDTD subgriding based on Huygens surfaces,

    J.-P. Berenger, “A FDTD subgriding based on Huygens surfaces,” in Proc. IEEE Ant. Propag. Society Int. Symp., Washington, USA, pp. 98–101, Dec. 2005

    work page 2005

  10. [10]

    A New Hybrid Implicit–Explicit FDTD Method for Local Subgridding in Multiscale 2-D TE Scattering Problems,

    A. V . Londersele, D. D. Zutter, and D. V . Ginste, “A New Hybrid Implicit–Explicit FDTD Method for Local Subgridding in Multiscale 2-D TE Scattering Problems,”IEEE Trans. Antennas Propag., vol. 64, no. 8, pp. 3509-3520, Aug. 2016

    work page 2016

  11. [11]

    FDTD Modeling of Nonlocality in a Nanoantenna Ac- celerated by a CPU–GPU Heterogeneous Architecture and Subgridding Techniques,

    J. Feng et al., “FDTD Modeling of Nonlocality in a Nanoantenna Ac- celerated by a CPU–GPU Heterogeneous Architecture and Subgridding Techniques,”IEEE Trans. Antennas Propag., vol. 72, no. 2, pp. 1708- 1720, Feb. 2024

    work page 2024

  12. [12]

    L. Ma, X. Ma, M. Chi, R. Xiang and X. Zhu, ”A Nonoverlapping Hybrid FDTD-PITD Method Using Upwind Conditions for Multiscale Problems,”IEEE Trans. Microw. Theory Techn., vol. 72, no. 5, pp. 2961- 2972, May 2024

    work page 2024

  13. [13]

    Reciprocity principle for stable subgrid- ding in the finite difference time domain method,

    L. Kulas and M. Mrozowski, “Reciprocity principle for stable subgrid- ding in the finite difference time domain method,” inProc. EUROCON Int. Conf. Comput. Tool, pp. 106-111, 2007

    work page 2007

  14. [14]

    EM interaction of handset antennas and a human in personal communications,

    M. Jensen and Y . Rahmat-Samii, “EM interaction of handset antennas and a human in personal communications,”Proc. IEEE, vol. 83, no. 1, pp. 7-17, Jan. 1995

    work page 1995

  15. [15]

    A novel method to analyze electromag- netic scattering of complex objects,

    K. Umashankar and A. Taflove, “A novel method to analyze electromag- netic scattering of complex objects,”IEEE Trans. Electromagn. Compat., vol. EMC-24, no. 4, pp. 397-405, Nov. 1982

    work page 1982

  16. [16]

    Full-wave analysis of pack- aged microwave circuits with active and nonlinear devices: An FDTD approach,

    C.-N. Kuo, B. Houshmand, and T. Itoh, “Full-wave analysis of pack- aged microwave circuits with active and nonlinear devices: An FDTD approach,”IEEE Trans. Microw. Theory Techn., vol. 45, no. 5, pp. 819–826, May 1997

    work page 1997

  17. [17]

    An Unsymmetric FDTD Subgridding Algorithm With Unconditional Stability,

    J. Yan and D. Jiao, “An Unsymmetric FDTD Subgridding Algorithm With Unconditional Stability,”IEEE Trans. Antennas Propag., vol. 66, no. 8, pp. 4137-4150, Aug. 2018

    work page 2018

  18. [18]

    Symmetric Positive Semi-Definite FDTD Sub- gridding Algorithms in Both Space and Time for Accurate Analysis of Inhomogeneous Problems,

    K. Zeng and D. Jiao, “Symmetric Positive Semi-Definite FDTD Sub- gridding Algorithms in Both Space and Time for Accurate Analysis of Inhomogeneous Problems,”IEEE Trans. Antennas Propag., vol. 68, no. 4, pp. 3047-3059, Apr. 2020. 13

    work page 2020

  19. [19]

    A symmetric FDTD subgridding method with guaranteed stability and arbitrary grid ratio,

    L. Deng, Y . Wang, C. Tian, H. Liu, X. Wang, G. Chen, X.H. Wang, X. Zhang, Z. Chen and S. Yang, “A symmetric FDTD subgridding method with guaranteed stability and arbitrary grid ratio,”IEEE Trans. Antennas Propag., vol. 71, no. 12, pp. 9207-9221, Dec. 2023

    work page 2023

  20. [20]

    A stable hybrid method for hyperbolic problems,

    J. Nordstrom and G. Jing, “A stable hybrid method for hyperbolic problems,”J. Comput. Phys., vol. 212, no. 2, pp. 436-453, 2006

    work page 2006

  21. [21]

    A Dissipation Theory for Three-Dimensional FDTD With Application to Stability Analysis and Subgridding,

    F. Bekmambetova, X. Zhang, and P. Triverio, “A Dissipation Theory for Three-Dimensional FDTD With Application to Stability Analysis and Subgridding,”IEEE Trans. Antennas Propag., vol. 66, no. 12, pp. 7156-7170, Dec. 2018

    work page 2018

  22. [22]

    FDTD Maxwell’s equations models for nonlinear electrodynamics and optics,

    R. M. Joseph and A. Taflove, “FDTD Maxwell’s equations models for nonlinear electrodynamics and optics,”IEEE Trans. Antennas Propag., vol. 45, no. 3, pp. 364–374, Mar. 1997

    work page 1997

  23. [23]

    A new subgridding method for the finite-difference time-domain (FDTD) algorithm,

    W. Yu and R. Mittra, “A new subgridding method for the finite-difference time-domain (FDTD) algorithm,”IEEE Antennas Propag. Mag., vol. 39, no. 5, pp. 108-112, Oct. 1997

    work page 1997

  24. [24]

    Stable coupling of nonconforming, high- order finite difference methods,

    J.E. Kozdon, L.C. Wilcox, “Stable coupling of nonconforming, high- order finite difference methods,”SIAM J. Sci. Comput., vol. 38, no. 2, pp. 923-952, 2016

    work page 2016

  25. [25]

    Finite element and finite difference methods for hyperbolic partial differential equations,

    H.-O. Kreiss, G. Scherer, “Finite element and finite difference methods for hyperbolic partial differential equations,” inMathematical Aspects of Finite Elements in Partial Differential Equations, Academic Press, pp. 195–212, 1974

    work page 1974

  26. [26]

    Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems,

    M. H. Carpenter, D. Gottlieb, and S. Abarbanel, “Time-stable boundary conditions for finite-difference schemes solving hyperbolic systems,”J. Comput. Phys., vol. 111, no. 2, pp. 220–236, Apr. 1994

    work page 1994

  27. [27]

    Boundary and Interface Conditions for High-Order Finite-Difference Methods Applied to the Euler and Navier–Stokes Equations,

    J. Nordstrom and M. Carpenter, “Boundary and Interface Conditions for High-Order Finite-Difference Methods Applied to the Euler and Navier–Stokes Equations,”J. Comput. Phys., vol. 148, no. 2, pp. 621–645, Jan. 1999

    work page 1999

  28. [28]

    High order finite difference method for wave propagation in discontinuous media,

    K. Mattsson and J. Nordstrom, “High order finite difference method for wave propagation in discontinuous media,”J. Comput. Phys., vol. 220, no. 1, pp. 249–269, Dec. 2006

    work page 2006

  29. [29]

    A Stable FDTD Subgridding Scheme with SBP-SAT for Transient Electromagnetic Analysis,

    Y . Cheng, Y . Wang, H. Liu, L. Li, X. Wang, S. Yang and Z. Chen, “A Stable FDTD Subgridding Scheme with SBP-SAT for Transient Electromagnetic Analysis,”J. Comput. Phys., vol. 494, 112510, Dec. 2023

    work page 2023

  30. [30]

    Toward the Development of a 3-D SBP-SAT FDTD Method: Theory and Validation,

    Y . Cheng, H. Liu, X. Wang, G. Chen, X. Wang, X. Zhang, S. Yang and Z. Chen, “Toward the Development of a 3-D SBP-SAT FDTD Method: Theory and Validation,”IEEE Trans. Antennas Propag., vol. 71, no. 12, pp. 9178-9193, Dec. 2023

    work page 2023

  31. [31]

    Towards the Development of a Three-Dimensional SBP-SAT FDTD Method: Subgridding Implementation,

    H. Liu, Y . Wang, L.Deng, G. Chen, A. Chen, Z. Chen and S. Yang, “Towards the Development of a Three-Dimensional SBP-SAT FDTD Method: Subgridding Implementation,”IEEE Trans. Antennas Propag., vol. 73, no. 5, pp. 3120-3132, May 2025

    work page 2025

  32. [32]

    An SBP- SAT FDTD Subgridding Method Using Staggered Yee’s Grids Without Modifying Field Components for TM Analysis,

    Y . Wang, Y . Cheng, X. Wang, S. Yang and Z. Chen, “An SBP- SAT FDTD Subgridding Method Using Staggered Yee’s Grids Without Modifying Field Components for TM Analysis,”IEEE Trans. Microw. Theory Techn., vol. 71, no. 2, pp. 579-592, Feb. 2023

    work page 2023

  33. [33]

    Toward the 2-D Stable FDTD Subgridding Method With SBP-SAT and Arbitrary Grid Ratio,

    Y . Wang, L. Deng, H. Liu, Z. Chen and S. Yang, “Toward the 2-D Stable FDTD Subgridding Method With SBP-SAT and Arbitrary Grid Ratio,” IEEE Trans. Microw. Theory Techn., vol. 72, no. 3, pp. 1591-1605, Mar. 2024

    work page 2024

  34. [34]

    Stable difference methods for block-oriented adaptive grids,

    A. Nissen, K. Kormann, M. Grandin and K. Virta, “Stable difference methods for block-oriented adaptive grids,”J. Sci. Comput., vol. 65, no. 2, pp. 486–511, Mar. 2015

    work page 2015

  35. [35]

    Stable and accurate interpolation operators for high-order multiblock finite difference methods,

    K. Mattsson and M. Carpenter, “Stable and accurate interpolation operators for high-order multiblock finite difference methods,”SIAM J. Sci. Comput., vol. 32, no. 4, pp. 2298-2320, 2010

    work page 2010

  36. [36]

    High order finite difference methods for the wave equation with non-conforming grid interfaces,

    S. Wang, K. Virta and G. Kreiss, “High order finite difference methods for the wave equation with non-conforming grid interfaces,”J. Sci. Comput., vol. 68, pp. 1002-1028, 2016

    work page 2016

  37. [37]

    Summation-by-parts operators for non- simply connected domains,

    S. Nikkar and J. Nordstrom, “Summation-by-parts operators for non- simply connected domains,”SIAM J. Sci. Comput., vol. 40, no. 3, pp. 1250-1273, 2018

    work page 2018

  38. [38]

    A Stable SBP-SAT FDTD Subgridding Method for T-junction Blocks

    Y . Wang, L. Deng, W. Wu, H. Liu, X. Zhang, X. Zhang, J. Wang, Z. Chen and S. Yang, “A Stable SBP-SAT FDTD Subgridding Method for T-junction Blocks”,IEEE Trans. Microw. Theory Techn., vol. 73, no. 11, pp. 8662-8677, Nov. 2025

    work page 2025