In order to better understand RLC behavior under current sources, I constructed the simple circuit below and attempted to analyze it. However, my equations do not agree with simulation results.
The equations are simple, and use standard calculus, yet they don't seem to produce correct results.
My analysis is in some ways consistent with simulation but in some ways gets very different results. Is my analysis correct? If so, why does it not agree with simulation?
I've already asked about this on an electronics site. While they confirmed that simulation does not agree with the equations, no one was able to spot any error in them or suggest more correct ones. So I decided I needed to ask here.
Analysis shows: $$ V_{out} = I_{src} \omega^2 LCR_{load}$$ since
$$ I_{src} = I_0 e^{j \omega t} \\ V_{L1} = V_1 - L \frac{dI_{src}}{dt} \\ = V_1 -jI_{src} \omega L$$ and $$ I_{load} = C \frac{dV_{L1}}{dt} \\ = \omega^2 L C I_{src}$$ so $$V_{out} = I_{load}R_{load} \\ = I_{src} \omega^2 LCR_{load}. $$
This suggests that $V_{out}$ is linear in $L, C, R$ and quadratic in $f$ and should be in phase with $I_{src}$. However, when I simulate it, I find that similar relationships, but not exact, and that $V_{out}$ hits a limit at a few MHz after which increasing $f$ has no effect. Why?