pith. sign in

arxiv: 2605.21450 · v1 · pith:3CWU4CL6new · submitted 2026-05-20 · 🧮 math.NA · cs.NA

Composite B-Spline Current Deposition and Interpolation Operators for Thin-Wire Finite-Difference Time-Domain Simulations

Cole Gruninger , Boyce E. Griffith This is my paper

Pith reviewed 2026-05-21 02:53 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords thin-wire FDTDB-spline regularizationcharge conservationcurrent depositioninterpolation operatorparasitic currentsfinite difference time domaindiscrete divergence-free
0
0 comments X

The pith

Composite B-spline regularizations enforce discrete charge conservation in thin-wire FDTD current deposition.

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

Thin-wire finite-difference time-domain simulations of obliquely oriented closed-loop antennas produce unphysical low-frequency parasitic currents when the current-deposition operator fails to conserve charge. The paper shows that charge conservation requires the deposited current to be discretely divergence-free for any constant current along the wire. It constructs a family of composite B-spline regularizations that meet this condition to machine precision by using piecewise-polynomial kernels whose line integrals can be evaluated exactly with composite Gauss-Legendre quadrature. Making the electric-field interpolation operator the discrete adjoint of the deposition operator preserves skew-symmetry, so that an irrotational field produces no net electromotive force around a closed loop. Tests on dipoles and square and circular loops confirm orientation-independent impedances free of the parasitic currents that appear with naive trilinear regularization.

Core claim

The authors introduce composite B-spline regularizations of distributions supported on the wire that make the current-deposition operator discretely divergence-free to machine precision whenever the wire carries constant current. Exact evaluation of the coupling integrals is obtained because the B-spline kernels are piecewise polynomial with a priori known breakpoints, allowing composite Gauss-Legendre quadrature at every grid-plane crossing. The interpolation operator is defined as the discrete adjoint of deposition, preserving skew-symmetry so that discretely irrotational electric fields drive no net electromotive force around closed loops.

What carries the argument

Composite B-spline regularizations that render the current-deposition operator discretely divergence-free for constant currents, with the interpolation operator taken as its discrete adjoint to enforce skew-symmetry.

If this is right

  • Impedance values become independent of wire orientation and match known analytical characteristics for dipoles and loop antennas.
  • Unphysical parasitic low-frequency currents disappear in closed loops.
  • Skew-symmetry ensures an irrotational electric field produces zero net electromotive force around any closed wire path.
  • Line integrals are evaluated exactly via composite Gauss-Legendre quadrature on known breakpoints.

Where Pith is reading between the lines

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

  • The same regularization approach could be applied to other grid-based methods that require discrete conservation along embedded curves.
  • Standard FDTD codes could adopt these operators to simulate arbitrarily oriented thin-wire structures without manual alignment to the grid.
  • Higher-order extensions of the composite B-splines might further reduce dispersion errors while retaining exact charge conservation.

Load-bearing premise

Charge conservation requires the deposited current to be discretely divergence-free when the wire carries a constant current.

What would settle it

A closed-loop antenna simulation performed with the composite B-spline operators that still exhibits persistent low-frequency parasitic currents would falsify the discrete divergence-free property.

Figures

Figures reproduced from arXiv: 2605.21450 by Boyce E. Griffith, Cole Gruninger.

Figure 1
Figure 1. Figure 1: Input impedance of the center-fed dipole versus dipole length in wavelengths 𝐿dipole/𝜆, with resistance 𝑅in as solid lines and reactance 𝑋in as dashed lines. Subpanels (a)–(d) use the four regularized delta functions 𝛿 (0) ℎ , 𝛿 (2) ℎ , 𝛿 (4) ℎ , and the isotropic 𝛿 iso ℎ . Within each subpanel, three orientation curves are overlaid: axis-aligned 𝜉ˆ = ˆ𝑧, face-diagonal 𝜉ˆ = (𝑥ˆ + 𝑦ˆ)/√ 2, and body-diagonal… view at source ↗
Figure 2
Figure 2. Figure 2: Input impedance of the circular loop versus circumference in wavelengths 𝐶/𝜆, with resistance 𝑅in as solid lines and reactance 𝑋in as dashed lines. Subpanels (a)–(d) use the four regularized delta functions 𝛿 (0) ℎ , 𝛿 (2) ℎ , 𝛿 (4) ℎ , and the isotropic 𝛿 iso ℎ . Within each subpanel, three orientation curves are overlaid: axis-aligned 𝑛ˆ = ˆ𝑧, face-diagonal 𝑛ˆ = (𝑥ˆ + 𝑦ˆ)/√ 2, and body-diagonal 𝑛ˆ = (𝑥ˆ … view at source ↗
Figure 3
Figure 3. Figure 3: Gap-current time history 𝐼gap (𝑡) at the feed of the circular loop. Subpanels (a)–(d) use the four regularized delta functions 𝛿 (0) ℎ , 𝛿 (2) ℎ , 𝛿 (4) ℎ , and the isotropic 𝛿 iso ℎ , with the three orientation curves overlaid in each. 6.3 Square Loop Antenna In the final experiment, we replace the circular loop with a square loop antenna of side length 1 m, fed at the midpoint of one side by the same dif… view at source ↗
Figure 4
Figure 4. Figure 4: Input impedance of the square loop versus perimeter in wavelengths 𝐶/𝜆, with resistance 𝑅in as solid lines and reactance 𝑋in as dashed lines. Subpanels (a)–(d) use the four regularized delta functions 𝛿 (0) ℎ , 𝛿 (2) ℎ , 𝛿 (4) ℎ , and the isotropic 𝛿 iso ℎ . Within each subpanel, three orientation curves are overlaid: axis-aligned 𝑛ˆ = ˆ𝑧, face-diagonal 𝑛ˆ = (𝑥ˆ + 𝑦ˆ)/√ 2, and body-diagonal 𝑛ˆ = (𝑥ˆ + 𝑦ˆ +… view at source ↗
Figure 5
Figure 5. Figure 5: Gap-current time history 𝐼gap (𝑡) at the feed of the square loop. Subpanels (a)–(d) use the four regularized delta functions 𝛿 (0) ℎ , 𝛿 (2) ℎ , 𝛿 (4) ℎ , and the isotropic 𝛿 iso ℎ , with the three orientation curves overlaid in each. The interpolation operator is taken to be the discrete adjoint of the current deposition operator. By orthogonality of the gradient and curl subspaces in the discrete Helmhol… view at source ↗
read the original abstract

Holland-Simpson thin-wire finite-difference time-domain (FDTD) simulations of obliquely oriented closed-loop antennas exhibit persistent low-frequency parasitic currents because the current-deposition operator fails to conserve charge. Together with an interpolation operator that samples the tangential electric field along the wire, this deposition operator can be realized as a regularization of distributions against a regularized delta function supported on the wire. We show that charge conservation requires the deposited current to be discretely divergence-free when the wire carries a constant current, and we introduce a family of composite B-spline regularizations that satisfy this condition to machine precision. Exact evaluation of the coupling line integrals is achievable because the B-spline kernels are piecewise polynomial with breakpoints known a priori, allowing composite Gauss-Legendre quadrature with subinterval breakpoints at every grid-plane crossing. Taking the interpolation operator as the discrete adjoint of the deposition operator preserves skew-symmetry and ensures that a discretely irrotational electric field drives no net electromotive force around a closed loop. Numerical experiments on a center-fed dipole and on circular and square loop antennas show that the proposed regularizations yield orientation-independent impedance values consistent with known characteristics, whereas a naive trilinear regularization produces unphysical parasitic low-frequency currents in closed loops.

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 paper introduces composite B-spline regularization operators for current deposition and tangential E-field interpolation in Holland-Simpson thin-wire FDTD. Starting from the requirement that constant current on an obliquely oriented wire must produce a discretely divergence-free deposited current (to enforce charge conservation), the authors construct piecewise-polynomial kernels whose line integrals against the grid are evaluated exactly via composite Gauss-Legendre quadrature with breakpoints at every grid-plane crossing. The interpolation operator is defined as the discrete adjoint of the deposition operator, guaranteeing skew-symmetry so that a discretely irrotational E produces zero loop EMF. Numerical experiments on a center-fed dipole and on circular and square loops demonstrate orientation-independent impedances free of the low-frequency parasitic currents that appear with trilinear deposition.

Significance. If the central construction holds, the work supplies a parameter-free, machine-precision charge-conserving pair of operators that directly resolves a documented defect in thin-wire FDTD modeling of closed-loop antennas. The combination of exact quadrature on known breakpoints, adjoint symmetry, and reproducible numerical validation constitutes a clear methodological advance for computational electromagnetics.

major comments (1)
  1. [Derivation of discrete charge conservation] The central claim that the composite B-spline deposition operator is discretely divergence-free to machine precision rests on the exact cancellation of the quadrature error for the piecewise-polynomial kernel. Please supply the explicit error bound or cancellation argument (e.g., in the section deriving the divergence-free property) showing that the composite Gauss-Legendre rule with a-priori grid-plane breakpoints achieves exact integration for the B-spline support.
minor comments (2)
  1. [Numerical experiments] Table or figure comparing impedance spectra for the three loop orientations would benefit from explicit statement of the grid resolution (cells per wavelength) and Courant number used in all runs.
  2. Notation for the composite quadrature weights and breakpoint locations could be introduced once in a dedicated subsection rather than inline, to improve readability for readers implementing the method.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful and constructive review, as well as the recommendation for minor revision. We address the single major comment below and will revise the manuscript accordingly.

read point-by-point responses

  1. Referee: [Derivation of discrete charge conservation] The central claim that the composite B-spline deposition operator is discretely divergence-free to machine precision rests on the exact cancellation of the quadrature error for the piecewise-polynomial kernel. Please supply the explicit error bound or cancellation argument (e.g., in the section deriving the divergence-free property) showing that the composite Gauss-Legendre rule with a-priori grid-plane breakpoints achieves exact integration for the B-spline support.

    Authors: We thank the referee for this request. The manuscript already notes that exact integration follows from the piecewise-polynomial character of the kernels together with a-priori breakpoints at grid-plane crossings, but we agree that an explicit argument is desirable. In the revised manuscript we will insert a short paragraph immediately after the definition of the composite deposition operator (in the section deriving discrete charge conservation). The added text will state: the composite B-spline kernel of order k is piecewise polynomial of degree k-1 with known integer knots; subdividing the integration path at every grid-plane crossing produces subintervals on which the kernel is a single polynomial of degree at most k-1; a Gauss-Legendre rule with n nodes integrates polynomials of degree up to 2n-1 exactly; choosing n = ceil(k/2) therefore renders the quadrature exact on every subinterval. Consequently the composite rule evaluates the line integrals against the grid basis functions to machine precision with no residual quadrature error, so that a constant current on the wire produces a discretely divergence-free deposition. We will also cite the standard result on exactness of Gaussian quadrature for polynomials on each piece. This supplies the explicit cancellation-free argument requested. revision: yes

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper begins from the external physical requirement of discrete charge conservation (deposited current must be discretely divergence-free for constant current on the wire) and constructs composite B-spline operators plus adjoint interpolation to satisfy it exactly via piecewise-polynomial kernels and a-priori breakpoint quadrature. This is a direct implementation of the stated condition rather than a fit, self-definition, or reduction to prior self-citations. The skew-symmetry preservation follows from taking the discrete adjoint, and numerical tests on dipoles/loops serve as validation, not definition. The derivation chain is self-contained against the external benchmark of charge conservation and does not reduce to its own inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the domain requirement that deposited current must be discretely divergence-free for constant wire current, plus the assumption that B-spline kernels permit exact line-integral evaluation via a priori breakpoint quadrature. No free parameters or invented entities are introduced in the abstract.

axioms (1)
  • domain assumption The deposited current must be discretely divergence-free when the wire carries a constant current.
    Stated directly in the abstract as the condition required for charge conservation.

pith-pipeline@v0.9.0 · 5757 in / 1249 out tokens · 35005 ms · 2026-05-21T02:53:56.951476+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

13 extracted references · 13 canonical work pages

  1. [1]

    Holland and L

    1R. Holland and L. Simpson. Finite-difference analysis of EMP coupling to thin struts and wires.IEEE Transactions on Electromagnetic Compatibility, EMC-23(2):88–97, May 1981.doi:10.1109/TEMC. 1981.303899. 2K. S. Yee. Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media.IEEE Transactions on Antennas and Pro...

    work page doi:10.1109/temc 1981

  2. [2]

    doi:10.1109/TAP.1966.1138693. 3A. Taflove and S. C. Hagness.Computational Electrodynamics: The Finite-Difference Time-Domain Method. Artech House, Norwood, MA, USA, 3rd edition,

    work page doi:10.1109/tap.1966.1138693 1966

  3. [3]

    6Jean-Pierre Bérenger

    doi:10.1109/TEMC.2003.811294. 6Jean-Pierre Bérenger. Origin of Parasitic Solutions With Holland and Simpson Thin Wires in the FDTD Grid.IEEE Transactions on Electromagnetic Compatibility, 61(2):487–494, April 2019.doi: 10.1109/TEMC.2018.2814341. 7Christophe Guiffaut and Alain Reineix. Cartesian Shift Thin Wire Formalism in the FDTD Method With Multiwire J...

    work page doi:10.1109/temc.2003.811294 2003

  4. [4]

    13ColeGruningerandBoyceE.Griffith

    doi:10.1016/ 0010-4655(92)90169-Y. 13ColeGruningerandBoyceE.Griffith. CompositeB-splineregularizeddeltafunctionsfortheimmersed boundary method: Divergence-free interpolation and gradient-preserving force spreading.Journal of Computational Physics, 546:114472, February 2026.doi:10.1016/j.jcp.2025.114472. 14Lianxia Li, Cole Gruninger, Jae H. Lee, and Boyce ...

    work page doi:10.1016/j.jcp.2025.114472 2026

  5. [5]

    18James M

    doi:10.1137/S0036142996314044. 18James M. Hyman and Mikhail Shashkov. Mimetic Discretizations for Maxwell’s Equations.Journal of Computational Physics, 151(2):881–909, May 1999.doi:10.1006/jcph.1999.6225. 19Douglas Arnold, Richard Falk, and Ragnar Winther. Finite element exterior calculus: From Hodge theory to numerical stability.Bulletin of the American ...

    work page doi:10.1137/s0036142996314044 1999

  6. [6]

    doi:10.1002/(SICI)1099-1204(199901/04)12:1/2<129::AID-JNM327>3.0.CO; 2-G. 23A. Bossavit and L. Kettunen. Yee-like schemes on staggered cellular grids: A synthesis between FIT and FEM approaches.IEEE Transactions on Magnetics, 36(4):861–867, July 2000.doi:10.1109/20. 877580. 18 24Carl De Boor.A Practical Guide to Splines. Number 27 in Applied Mathematical ...

    work page doi:10.1002/(sici)1099-1204(199901/04)12:1/2 2000

  7. [7]

    print edition, 1991.doi:10.1007/978-1-4612-6333-3. 25A. Buffa, G. Sangalli, and R. Vázquez. Isogeometric analysis in electromagnetics: B-splines approxi- mation.Computer Methods in Applied Mechanics and Engineering, 199(17):1143–1152, March

    work page doi:10.1007/978-1-4612-6333-3 1991

  8. [8]

    doi:10.1016/j.cma.2009.12.002. 26R. N. Simpson, Z. Liu, R. Vázquez, and J. A. Evans. An isogeometric boundary element method for electromagnetic scattering with compatible B-spline discretizations.Journal of Computational Physics, 362:264–289, June 2018.doi:10.1016/j.jcp.2018.01.025. 27Ritoban Roy-Chowdhury, Tamar Shinar, and Craig Schroeder. Higher order...

    work page doi:10.1016/j.cma.2009.12.002 2009

  9. [9]

    doi:10.1016/j.jcp.2024.112831. 28S. D. Gedney.Introduction to the Finite-Difference Time-Domain (FDTD) Method for Electromagnetics. Morgan & Claypool, San Rafael, CA, USA, 2011.doi:10.1007/978-3-031-01712-4. 29Jean-Pierre Berenger. Stability of the FDTD Method With Holland and Simpson Thin Wires.IEEE Transactions on Electromagnetic Compatibility, 63(3):88...

    work page doi:10.1016/j.jcp.2024.112831 2024

  10. [10]

    Solving Maxwell’s Equations Using the Ultra Weak Variational Formulation.Journal of Computational Physics, 223(2):731–758, 2007

    31Björn Engquist, Anna-Karin Tornberg, and Richard Tsai. Discretization of Dirac delta functions in level set methods.Journal of Computational Physics, 207(1):28–51, July 2005.doi:10.1016/j.jcp. 2004.09.018. 32Bamdad Hosseini, Nilima Nigam, and John M. Stockie. On regularizations of the Dirac delta distribution.Journal of Computational Physics, 305:423–44...

    work page doi:10.1016/j.jcp 2005

  11. [11]

    35Boyce E

    doi:10.1137/110836699. 35Boyce E. Griffith and Charles S. Peskin. On the order of accuracy of the immersed boundary method: Higher order convergence rates for sufficiently smooth problems.Journal of Computational Physics, 208(1):75–105, September 2005.doi:10.1016/j.jcp.2005.02.011. 36James E. Storer. Impedance of thin wire loop antenna. Technical Report T...

    work page doi:10.1137/110836699 2005

  12. [12]

    39Warren L

    doi:10.1063/1.4816619. 39Warren L. Stutzman and Gary A. Thiele.Antenna Theory and Design. Wiley, Hoboken, NJ,

    work page doi:10.1063/1.4816619

  13. [13]

    40R. King. Theory of the corner-driven square loop antenna.IRE Transactions on Antennas and Propagation, 4(3):393–407, July 1956.doi:10.1109/TAP.1956.1144431. 20 S1 Discrete Energy Conservation of the Leapfrog FDTD–Thin-Wire System We prove that, in the absence of absorbing boundary conditions, the leapfrog time-stepping scheme combined with the completel...

    work page doi:10.1109/tap.1956.1144431 1956