Co-Simulation
Understanding Co-Simulation
RF design exists at the boundary between two physical domains. Active devices (transistors, diodes) are described by compact SPICE models with lumped parameters (gm, Cgs, rds), while passive structures (transmission lines, inductors, package interconnects, antennas) exhibit distributed electromagnetic behavior governed by Maxwell's equations. At low frequencies where all dimensions are much smaller than a wavelength, lumped circuit models approximate the distributed behavior adequately. But as frequencies increase, this approximation breaks down: a 5 mm bond wire is electrically invisible at 100 MHz (λ/600) but becomes a significant inductor at 5 GHz (λ/12) and a resonant element at 30 GHz (λ/2).
Co-simulation bridges these domains by letting each solver handle its strengths. The EM solver uses finite elements or finite differences to solve Maxwell's equations in the 3D geometry, capturing distributed transmission line effects, electromagnetic coupling between adjacent traces, substrate currents and eddy effects, and radiation from open structures. The result is a frequency-dependent S-parameter matrix that accurately represents the passive structure's behavior at all ports. The circuit simulator then uses these S-parameters as a linear n-port network alongside nonlinear device models, applying harmonic balance, transient, or envelope analysis to predict the complete system's performance under realistic operating conditions.
Co-Simulation Interface Equations
bi = ∑j Sij(f) × aj (frequency domain, linear)
Rational Function Fit (macro-model):
S(s) = D + ∑k Rk/(s - pk) (poles pk, residues Rk)
Mesh Density Rule:
Δx ≤ λmin/10 ; elements ≈ (L/Δx)3
Where aj, bi = incident/reflected power waves, s = jω = complex frequency, λmin = wavelength at maximum frequency. A 10 mm structure at 60 GHz (λ = 5 mm): Δx = 0.5 mm, ~8,000 elements. At 300 GHz: ≈1M elements.
Co-Simulation Tool Combinations
| EM Solver | Method | Circuit Simulator | Interface | Use Case |
|---|---|---|---|---|
| Ansys HFSS | FEM (3D) | Keysight ADS | S-parameter / Dynamic Link | MMIC, package, antenna-PA |
| CST Studio | FDTD / FIT | CST Design Studio | Built-in co-simulation | SI/PI, EMC, antenna |
| Keysight EMPro | FEM / FDTD | Keysight ADS | Native integration | Package, board-level |
| Cadence EMX | MoM (2.5D) | Cadence Spectre | CDF / Virtuoso | On-chip inductor, MIM cap |
| Sonnet | MoM (planar) | NI AWR / ADS | Touchstone export | Microstrip, filter, coupler |
Frequently Asked Questions
How does EM-circuit co-simulation work?
3D EM solver computes multiport S-parameters (FEM/FDTD/MoM) with ports at active circuit connections. S-parameter data (.snp Touchstone) or fitted rational function is imported into circuit simulator as n-port network. Circuit simulator connects it to transistor SPICE models and runs harmonic balance, transient, or envelope analysis. EM captures distributed effects; circuit handles nonlinearity.
When is co-simulation necessary?
Below 1 GHz: rarely needed (lumped models adequate). 1 to 10 GHz: recommended for critical paths (matching, packages, interconnects > λ/20). Above 10 GHz: essential for virtually all elements. Always needed for: on-chip inductors, BGA packages, antenna-PA integration, and filter-amplifier co-design where port impedance variation affects stability.
What are the computational costs?
EM scales O(N²) to O(N³) with mesh elements. Simple microstrip: seconds. QFN package at 28 GHz (500K elements): 4 to 8 hours on 64 cores. Main limitation: S-parameters are linear, so nonlinear EM effects (ferrite saturation, multipaction) need time-domain tight coupling. Workarounds: macro-modeling, adaptive sampling, model order reduction.