This article presents a novel waveform relaxation scheme to solve electromagnetically large structures loaded with lumped linear and nonlinear elements. The scheme partitions the problem into a linear electromagnetic structure and a possibly nonlinear lumped circuit, which are coupled using Thévenin interfaces across the steps of an iterative waveform relaxation scheme. The main novel contribution is an adaptive selection of the decoupling resistances used as port references to define incident and reflected scattering signals, whose time-domain samples are refined through iterations. The decoupling resistances are updated through iterations, with the main objective of improving convergence speed and ultimately runtime. The resulting scheme is self-adapting to terminations exploiting high dynamic range in their impedance profiles and is able to provide a suboptimal convergence rate. Three-dimensional shielding structures loaded with nonlinear elements are employed as numerical examples to demonstrate the proposed method.
