RLC circuit driven by a 'Real' toroid

1. Dec 20, 2012

blargh4fun

I'm currently wanting to model an RLC circuit that is driven by a toroid transformer with a driving current, I(t). See image below.

In physics class, I've done similar problems with 'ideal' toroids and a simple R circuit. What I'm struggling with is how to model a 'real' toroid transformer.

Normally I would do something like this:
$\mu_0 I_{enclosed} = \oint B \cdot dl$
$\Rightarrow B(t) = \frac{mu_0 I(t) N}{2 \pi r}$
where N is the number of turns.

The flux is equal to:
$\Phi = \int B \cdot da = \frac{mu_0 I(t) N}{2 \pi}\int \int dr dx= \frac{mu_0 I(t) N L}{2 \pi} ln\left(\frac{r_{out}}{r_{in}}\right )$
where L is the width of the toroid, and r is the inner and outer radius.

The EMF is:
$\varepsilon = -\frac{d\Phi}{dt} = - \frac{mu_0 N L}{2 \pi} ln\left(\frac{r_{out}}{r_{in}}\right ) \frac{dI(t)}{dt}$

Using Kirchoff's voltage law:
$\varepsilon = V_R + V_L + V_C = iR +\frac{di(t)}{dt}L + \frac{1}{C}\int_{-\inf}^{t}i(t)dt$

Differentiating and diving by L gives the ODE:
$\frac{d^2 i(t)}{d x^2}+\frac{R}{L}\frac{di(t)}{dt}+\frac{i(t)}{LC} = \frac{d \varepsilon (t)}{dt} \frac{1}{L}$

Which can then be solved for i(t), the induced current in the RLC circuit. Etc. Etc.

My Question Is:
Is this method ok for a 'real' toroid? I feel this is wrong because I need to include the permeability of the toroid, but I don't know how to do this. It was my understanding that the presence of a 'real' toroid strengthened the magnetic field (magnitude of 1000 or so) and subsequently increased the induced current. I will appreciate any feedback. Thanks!