Pi Aurigae

Harmonic balance is a method used to calculate the steady-state response of non-linear differential equations, and is mostly applied to electrical circuits . It is a frequency domain method for calculating the steady state, as opposed to the various time-domain steady state methods. The name "harmonic balance" is descriptive of the method, which starts with Kirchhoff's Current Law written in the frequency domain and a chosen number of harmonics. Effectively the method assumes the solution can be represented by a linear combination of sinusoids, then balances current and voltage sinusoids to satisfy Kirchhoff's law. The method is commonly used to simulate circuits which include nonlinear elements.

Microwave circuits were the original application for harmonic balance methods in electrical engineering. Microwave circuits were well-suited because, historically, microwave circuits consist of many linear components which can be directly represented in the frequency domain, plus a few nonlinear components. System sizes were typically small. For more general circuits, the method was considered impractical for all but these very small circuits until the mid-1990s, when Krylov subspace methods were applied to the problem. The application of preconditioned Krylov subspace methods allowed much larger systems to be solved, both in size of circuit and in numbers of harmonics. This made practical the present-day use of harmonic balance methods to analyze radio-frequency integrated circuits (RFICs).

Algorithm

The algorithm starts with Kirchhoff's current law written in the frequency-domain. To increase the efficiency of the procedure, the circuit may be partitioned into its linear and nonlinear parts, since the linear part is readily described and calculated using nodal analysis directly in the frequency domain.

First, an initial guess is made for the solution, then an iterative process continues:

1. Voltages $V$ are used to calculate the currents of the linear part, $I_{\text{linear}}$ in the frequency domain.
2. Voltages $V$ are then used to calculate the currents in the nonlinear part, $I_{\text{nonlinear}}$. Since nonlinear devices are described in the time domain, the frequency-domain voltages $V$ are transformed into the time domain, typically using inverse fast Fourier transforms. The nonlinear devices are then evaluated using the time-domain voltage waveforms to produce their time-domain currents. The currents are then transformed back into the frequency domain.
3. According to Kirchhoff's circuit laws, the sum of the currents must be zero, $ϵ=I_{\text{linear}}+I_{\text{nonlinear}}=0$. An iterative process, usually Newton iteration, is used to update the network voltages $V$ such that the current residual $ϵ$ is reduced. This step requires formulation of the Jacobian $\frac{dϵ}{dV}$.

Convergence is reached when $ϵ$ is acceptably small, at which point all voltages and currents of the steady-state solution are known, most often represented as Fourier coefficients.