pith. sign in

arxiv: 1907.03719 · v1 · pith:2YN4GMH2new · submitted 2019-07-08 · 💻 cs.CE · cs.NA· math.NA

Multirate PWM balance method for the efficient field-circuit coupled simulation of power converters

Andreas Pels , Herbert De Gersem , Ruth V. Sabariego , Sebastian Sch\"ops This is my paper

Pith reviewed 2026-05-25 00:46 UTC · model grok-4.3

classification 💻 cs.CE cs.NAmath.NA
keywords multirate simulationPWMfield-circuit couplingpower converterseigenfunctionspartial differential equationsswitch-mode converters
0
0 comments X

The pith

New PWM eigenfunctions decouple equation systems to enable efficient multirate simulation of field-circuit coupled power converters.

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

The paper develops a multirate PWM balance method for simulating switch-mode power converters that include both circuit equations and large field models. Standard time-stepping demands tiny steps to capture steep transients across all degrees of freedom, which becomes prohibitively slow. Multirate partial differential equations separate the problem into components with distinct time scales. A new set of PWM eigenfunctions is introduced that decouples the resulting algebraic systems. The approach therefore reduces computational effort for converters that use idealized switches.

Core claim

A set of new PWM eigenfunctions decouples the systems of equations obtained from multirate partial differential equations and thereby yields an efficient simulation of the field-circuit coupled problem for power converters with idealized switches.

What carries the argument

The PWM eigenfunctions that decouple the multirate PDE systems into independent equations for different time scales.

If this is right

  • Only the fast-scale components require fine time discretization while slow-scale components use coarser steps.
  • The total number of degrees of freedom that must be advanced at the smallest time step is reduced.
  • Field and circuit parts can be advanced at time steps matched to their respective physical scales.
  • The method remains applicable as long as the switches are treated as ideal.

Where Pith is reading between the lines

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

  • The same eigenfunction construction might be adapted to converters whose switches include small but nonzero resistances.
  • Decoupling via eigenfunctions could be tested on other periodic multirate electromagnetic problems such as rotating machines.
  • If the eigenfunctions are precomputed once per duty cycle, the method could be embedded in existing circuit-field co-simulation tools.
  • A direct comparison on a standard benchmark converter would quantify the achieved speedup factor.

Load-bearing premise

Converters must use idealized switches so that multirate partial differential equations can separate the solution into components of different time scales.

What would settle it

Apply the method to a converter whose switches are modeled with finite on/off resistances and compare the computed waveforms and run time against a conventional fine-step reference simulation; any loss of accuracy or speedup would falsify the efficiency claim.

Figures

Figures reproduced from arXiv: 1907.03719 by Andreas Pels, Herbert De Gersem, Ruth V. Sabariego, Sebastian Sch\"ops.

Figure 1
Figure 1. Figure 1: Simplified circuit of the buck converter in continuous conduction mode with [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Original PWM basis functions pk(τ ), k ∈ {0, 1, 2, 3, 4}. gk(τ ) as linear combinations of the PWM basis functions, i.e., gk(τ ) := X Np l=0 vk,l pl(τ ), (9) where vk,l are unknown coefficients with k ∈ {0, . . . , Np}, and gk(τ ) are eigenfunctions of the time derivative operator d dτ gk(τ ) = λk gk(τ ). (10) We enforce this property in a weak sense by a Galerkin approach, i.e., − Z 1 0 gk(τ ) dpm(τ ) dτ … view at source ↗
Figure 3
Figure 3. Figure 3: PWM eigenfunctions gk(τ ), k ∈ {0, 1, 2, 3, 4}, i.e., Np = 4. (top) real part. (bottom) imaginary part. where A as in (5), Be = J ⊗ B + Λ ⊗ A, (14) Ce(t1) = Z Ts 0 ¯g(τ (t2)) ⊗ bc(t1, t2)dt2 , (15) and Λ is a diagonal matrix with diagonal entries λ0, λ1, . . . , λNp . Thus the resulting matrices in (13) are block-diagonal and the degrees of freedom can be block-wisely decoupled. This leads to Np + 1 indepe… view at source ↗
Figure 4
Figure 4. Figure 4: Multivariate voltage at the capacitor calculated using the multirate PWM [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: (top) Reference solution calculated using conventional adaptive time dis [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Error  as defined in (27) over time for solving the systems of equations. The MPDE approach with PWM eigenfunctions (multirate PWM balance method) is considerably faster than the MPDE approach with the original PWM basis functions and the conventional time discretization. faster simulation. The overall accuracy of the method is problem-specific and always depends on both the tolerance for the solver and t… view at source ↗
Figure 7
Figure 7. Figure 7: Coefficients w1,k for the inductor current calculated by solving (13) with Np,pwmbal = 4. (top) real part. The coefficients w1,1, . . . , w1,4 are approx￾imately the same therefore they are hard to distinguish visually. (bottom) imaginary part. speed-up amounting to a factor 4 for the test example. 7 Acknowledgements This work is supported by the “Excellence Initiative” of German Federal and State Governme… view at source ↗
read the original abstract

The field-circuit coupled simulation of switch-mode power converters with conventional time discretization is computationally expensive since very small time steps are needed to appropriately account for steep transients occurring inside the converter, not only for the degrees of freedom (DOFs) in the circuit, but also for the large number of DOFs in the field model part. An efficient simulation technique for converters with idealized switches is obtained using multirate partial differential equations, which allow for a natural separation into components of different time scales. This paper introduces a set of new PWM eigenfunctions which decouple the systems of equations and thus yield an efficient simulation of the field-circuit coupled problem. The resulting method is called the multirate PWM balance method.

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

0 major / 2 minor

Summary. The manuscript introduces the multirate PWM balance method for efficient field-circuit coupled simulation of switch-mode power converters with idealized switches. It employs multirate partial differential equations to separate components across different time scales and defines a new set of PWM eigenfunctions that are claimed to decouple the resulting system of equations, thereby avoiding the need for very small time steps in both the circuit and field models.

Significance. If the claimed decoupling property of the PWM eigenfunctions holds and produces the stated efficiency gains without loss of accuracy, the approach would provide a practical advance for computational simulation of power electronics, where conventional time-stepping methods are often prohibitive due to the combination of large field DOFs and steep switching transients.

minor comments (2)
  1. [Abstract] Abstract: the description of the PWM eigenfunctions and the precise mechanism of decoupling would benefit from a short statement of their construction or orthogonality properties to orient the reader before the technical sections.
  2. The manuscript would be strengthened by the addition of at least one concrete numerical comparison (e.g., wall-clock time or DOF count versus conventional time discretization) to quantify the efficiency improvement asserted in the abstract.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive assessment of the manuscript and the recommendation for minor revision. No specific major comments were provided in the report.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper's central claim is the introduction of new PWM eigenfunctions that decouple multirate PDE systems arising from field-circuit coupling with idealized switches. No load-bearing step reduces by construction to a fitted input, self-definition, or self-citation chain; the derivation is presented as a direct consequence of the multirate PDE framework and the proposed eigenfunctions. The abstract and reader's summary indicate a self-contained methodological advance without the enumerated circularity patterns.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Review is based on abstract only; no free parameters, axioms, or invented entities are specified in the provided text.

pith-pipeline@v0.9.0 · 5661 in / 1055 out tokens · 24227 ms · 2026-05-25T00:46:45.076807+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

14 extracted references · 14 canonical work pages

  1. [1]

    Mohan, T

    N. Mohan, T. M. Undeland, and W. P. Robbins, Power electronics: converters, applications and design . Wiley, 3 ed., 2003

    work page 2003

  2. [2]

    On the numerical accuracy of electromagnetic transient simulation with power electronics,

    J. Tant and J. Driesen, “On the numerical accuracy of electromagnetic transient simulation with power electronics,” IEEE Trans. Power Deliv. , vol. 33, pp. 2492– 2501, Oct. 2018

    work page 2018

  3. [3]

    Solving nonlinear circuits 12 with pulsed excitation by multirate partial differential equations,

    A. Pels, J. Gyselinck, R. V. Sabariego, and S. Sch¨ ops, “Solving nonlinear circuits 12 with pulsed excitation by multirate partial differential equations,” IEEE Trans. Magn., vol. 54, pp. 1–4, Mar. 2018

    work page 2018

  4. [4]

    Efficient simulation of DC- DC switch-mode power converters by multirate partial differential equations,

    A. Pels, J. Gyselinck, R. V. Sabariego, and S. Sch¨ ops, “Efficient simulation of DC- DC switch-mode power converters by multirate partial differential equations,”IEEE J. Multiscale Multiphys. Comput. Tech. , vol. 4, no. 1, pp. 64–75, 2019

    work page 2019

  5. [5]

    Numerical steady state analysis of electronic circuits driven by multi-tone signals,

    H. G. Brachtendorf, G. Welsch, R. Laur, and A. Bunse-Gerstner, “Numerical steady state analysis of electronic circuits driven by multi-tone signals,” Electr. Eng. , vol. 79, no. 2, pp. 103–112, 1996

    work page 1996

  6. [6]

    Analyzing circuits with widely separated time scales using nu- merical PDE methods,

    J. Roychowdhury, “Analyzing circuits with widely separated time scales using nu- merical PDE methods,” IEEE Trans. Circ. Syst. Fund. Theor. Appl. , vol. 48, pp. 578–594, May 2001

    work page 2001

  7. [7]

    Using dedicated time-domain ba- sis functions for the simulation of pulse-width-modulation controlled devices – ap- plication to the steady-state regime of a buck converter,

    J. Gyselinck, C. Martis, and R. V. Sabariego, “Using dedicated time-domain ba- sis functions for the simulation of pulse-width-modulation controlled devices – ap- plication to the steady-state regime of a buck converter,” in Electromotion 2013, (Cluj-Napoca, Romania), Oct. 2013

    work page 2013

  8. [8]

    Winding functions in transient magne- toquasistatic field-circuit coupled simulations,

    S. Sch¨ ops, H. De Gersem, and T. Weiland, “Winding functions in transient magne- toquasistatic field-circuit coupled simulations,” COMPEL, vol. 32, pp. 2063–2083, Sept. 2013

    work page 2063

  9. [9]

    Lamour, R

    R. Lamour, R. M¨ arz, and C. Tischendorf, Differential-Algebraic Equations: A Projector Based Analysis . Differential-Algebraic Equations Forum, Heidelberg: Springer, 2013

    work page 2013

  10. [10]

    Multirate partial differential algebraic equa- tions for simulating radio frequency signals,

    R. Pulch, M. G¨ unther, and S. Knorr, “Multirate partial differential algebraic equa- tions for simulating radio frequency signals,” Eur. J. Appl. Math. , vol. 18, no. 6, pp. 709–743, 2007

    work page 2007

  11. [11]

    Adaptive multi-rate wavelet method for circuit simulation,

    K. Bittner and H. G. Brachtendorf, “Adaptive multi-rate wavelet method for circuit simulation,” Radioengineering, vol. 23, Apr. 2014

    work page 2014

  12. [12]

    Transient 3d eddy currents using modified magnetic vector potentials and magnetic scalar potentials,

    C. R. I. Emson and C. W. Trowbridge, “Transient 3d eddy currents using modified magnetic vector potentials and magnetic scalar potentials,” IEEE Trans. Magn. , vol. 24, pp. 86–89, Jan. 1988

    work page 1988

  13. [13]

    Implicit Runge-Kutta methods for transient magnetic field computation,

    A. Nicolet and F. Delinc´ e, “Implicit Runge-Kutta methods for transient magnetic field computation,” IEEE Trans. Magn., vol. 32, pp. 1405–1408, May 1996

    work page 1996

  14. [14]

    Structural analysis of electrical circuits including magnetoquasistatic devices,

    A. Bartel, S. Baumanns, and S. Sch¨ ops, “Structural analysis of electrical circuits including magnetoquasistatic devices,” APNUM, vol. 61, pp. 1257–1270, Sept. 2011. 13

    work page 2011