Engineering Solving Second Order Homogeneous Circuit Equations | RLC Homework Statement

[GreenPrint]

Homework Statement

https://imageshack.us/scaled/large/708/capturehdx.png

Homework Equations

https://imageshack.us/scaled/large/845/capturevtq.png

V = IR

https://imageshack.us/scaled/large/18/captureid.png

I_{C}(t) = C\frac{dV_{x}(t)}{dt}

https://imageshack.us/scaled/large/818/captureplr.png

V_{x}(t) = L\frac{di_{L}}{dt}

The Attempt at a Solution

I'll start by solving for the equation of the node voltage V_{x}(t)

\frac{V_{s} - V_{x}(t)}{R} = C_{1}\frac{d[V_{x}(t) - V_{s}]}{dt} + C_{2}\frac{d[V_{x} - 0]}{dt} + \frac{1}{L}∫_{0}^{t}[V_{x}(t) - 0]dt + I_{L}(t=0^{+})

I use the fact that I_{L}(t=0^{-}) = 0 A and that because the current through a inductor can't change instantaneously so I_{L}(t=0^{-}) = I_{L}(t=0^{+}) = 0 A. I also simplify some mathematically.

-\frac{V_{x}(t)}{R} = C_{1}\frac{dV_{x}(t)}{dt} + C_{2}\frac{dV_{x}}{dt} + \frac{1}{L}∫_{0}^{t}V_{x}(t)dt

Simplify some more

-\frac{V_{x}(t)}{R} = (C_{1} + C_{2})\frac{dV_{x}(t)}{dt} + \frac{1}{L}∫_{0}^{t}V_{x}(t)dt

Take the derivative with respect to time of both sides

-\frac{d}{dt}\frac{V_{x}(t)}{R} = \frac{d[(C_{1} + C_{2})\frac{dV_{x}(t)}{dt} + \frac{1}{L}∫_{0}^{t}V_{x}(t)dt]}{dt}

Simplify

-\frac{1}{R}\frac{V_{x}(t)}{dt} = (C_{1} + C_{2})\frac{d^{2}V_{x}(t)}{dt^{2}} + \frac{1}{L}V_{x}(t)

Set equal to zero

0 = (C_{1} + C_{2})\frac{d^{2}V_{x}(t)}{dt^{2}} + \frac{1}{R}\frac{V_{x}(t)}{dt} + \frac{1}{L}V_{x}(t)

Replace the constants with constants A, B, and C

A = C_{1} + C_{2}
B = \frac{1}{R}
C = \frac{1}{L}

0 = A\frac{d^{2}V_{x}(t)}{dt^{2}} + B\frac{V_{x}(t)}{dt} + CV_{x}(t)

This is a second order differential equation so it is known that

V_{x}(t) = Ke^{st}, \frac{dV_{x}(t)}{dt} = Kse^{st}, \frac{d^{2}V_{x}(t)}{dt^{2}} = Ks^{2}e^{st}

So plugging this into the equation

0 = AKs^{2}e^{st} + BKse^{st} + CKe^{st}

Divide through by ke^{st}

0 = As^{2} + Bs + C

This is a quadratic equation so it's known that

s = \frac{-B ± \sqrt{B^2 - 4AC}}{2A}

It can be seen from the equation for s that there are two values for s, I'll call them s_{1} and s_{2}. So my voltage equation becomes

V_{x}(t) = A_{1}e^{s_{1}t} + A_{2}e^{s_{2}t}

In order to solve for A_{1} and A_{2} I need to know the initial condition of the node voltage.

At time 0^{-} the circuit looks like the one below

https://imageshack.us/scaled/large/833/capturetiru.png

\sum Q = 0 = 4 uF(V_{x}(t = 0^{-}) - 10) + 1 uF(V_{x}(t = 0^{-}) - 0)

Divide through by uF

0 = 4V_{x}(t = 0^{-}) - 4(10) + V_{x}(t = 0^{-})

simplify

40 = 5V_{x}(t = 0^{-})

so

V_{x}(t = 0^{-}) = 8 V

I use the fact that the voltage across a capacitor can't change instantaneously so V_{x}(t = 0^{-}) = V_{x}(t = 0^{+}).

So when you plug in 0 into the voltage equation you get the following equation

V_{x}(t = 0^{+}) = A_{1} + A_{2}

Being that I have two unknowns I need another equation. I use the fact that the voltage through a capacitor can't change instantaneously

So it can be seen that \frac{dV_{x}(t)}{dt} represents the change in voltage across C_{2} which is zero. So I have the following second equation

0 = A_{1}s_{1} + A_{2}s_{2}

So I choose to solve for A_{1} and A_{2} using matrix and putting them in reduce row echelon form

1, 1, V_{x}(t = 0^{+})
s_{1}, s_{2}, 0

I divide row 2 by s_{1}

1, 1, V_{x}(t = 0^{+})
1, \frac{s_{2}}{s_{2}}, 0

I multiply row 1 by negative one and add it to row two to get a new row 2

1, 1, V_{x}(t = 0^{+})
0, \frac{s_{2}}{s_{2}}-1, -V_{x}(t = 0^{+})

I divide row two through by \frac{1}{\frac{s_{2}}{s_{2}}-1}

1, 1, V_{x}(t = 0^{+})
0, 1, -\frac{V_{x}(t = 0^{+})}{\frac{s_{2}}{s_{2}}-1}

I multiply row two by negative one and add to row one to get a new row one

1, 0, V_{x}(t = 0^{+}) + \frac{V_{x}(t = 0^{+})}{\frac{s_{2}}{s_{2}}-1}
0, 1, -\frac{V_{x}(t = 0^{+})}{\frac{s_{2}}{s_{2}}-1}

So I get the following equations, I factored the first one

A_{1} = V_{x}(t = 0^{+})(1 + \frac{1}{\frac{s_{2}}{s_{2}}-1})

A_{2} = -\frac{V_{x}(t = 0^{+})}{\frac{s_{2}}{s_{2}}-1}

The quantity \frac{s_{1}}{s_{2}} seems to be an interesting quantity.

\frac{s_{1}}{s_{2}} = \frac{-B - \sqrt{B^{2} - 4AC}}{-B + \sqrt{B^2 - 4AC}}

I multiply by a form of one \frac{-B - \sqrt{B^{2} - 4 AC}}{-B - \sqrt{B^{2} - 4 AC}} and get the following.

\frac{s_{1}}{s_{2}} = \frac{(B + \sqrt{B^{2} - 4AC})^{2}}{B^{2} - B^{2} + 4AC} = \frac{B^{2} + 2B\sqrt{B^{2} - 4AC} + B^{2} - 4AC}{4AC} = \frac{2B^{2} + 2B\sqrt{B^{2} - 4AC} - 4AC}{4AC} = \frac{B^{2}}{2AC} + \frac{B}{2AC}\sqrt{B^{2} - 4AC} - 1

So then

\frac{s_{2}}{s_{1}} - 1 = \frac{B}{2AC}\sqrt{B^{2} - 4AC} + \frac{B^{2}}{2AC} - 2

Now in my formulas I have the following relation \frac{1}{\frac{s_{2}}{s_{1}}-1} So I solve for this relation.

\frac{1}{\frac{s_{2}}{s_{1}}-1} = \frac{1}{\frac{B}{2AC}\sqrt{B^{2} - 4AC} +( \frac{B^{2}}{2AC} - 2)}

I multiply by a form of one \frac{-\frac{B}{2AC}\sqrt{B^{2}-4AC} + \frac{B^2}{2AC} - 2}{-\frac{B}{2AC}\sqrt{B^{2}-4AC} + \frac{B^2}{2AC} - 2}

\frac{1}{\frac{s_{2}}{s_{1}}-1} = \frac{-\frac{B}{2AC}\sqrt{B^{2}-4AC} + \frac{B^2}{2AC} - 2}{-\frac{B^{2}}{4A^{2}C^{2}}(B^{2} - 4AC) + (\frac{B^{2}}{2AC} - 2)^{2}} = \frac{-\frac{B}{2AC}\sqrt{B^{2}-4AC} + \frac{B^2}{2AC} - 2}{-\frac{B^{4}}{4A^{2}C^{2}} + \frac{B^{2}}{AC} + \frac{B^{4}}{4A^{2}C^{2}} - \frac{4B^{2}}{2AC} + 4} = \frac{-\frac{B}{2AC}\sqrt{B^{2}-4AC} + \frac{B^{2}}{2AC} - 2}{-\frac{B^2}{AC} + 4}

At this point I multiply through by a form of one \frac{AC}{AC}

\frac{-\frac{B}{2}\sqrt{B^{2} - 4AC} + \frac{B^{2}}{2}-2AC}{4AC - B^{2}}

I multiply through by another form of one \frac{2}{2}

\frac{-B\sqrt{B^{2} - 4AC} + B^{2} - 4AC}{2(4AC - B^{2})}

Simplify

\frac{B\sqrt{B^{2}-4AC} + B^{2}}{2(4AC - B^{2})} - \frac{2AC}{4AC - B^{2}}

My equations for A_{1} and A_{2} then become

A_{1} = V_{x}(t = 0^{+})(1 + \frac{B\sqrt{B^{2}-4AC} + B^{2}}{2(4AC - B^{2})} - \frac{2AC}{4AC - B^{2}})

A_{2} = -V_{x}(t = 0^{+})(\frac{B\sqrt{B^{2}-4AC} + B^{2}}{2(4AC - B^{2})} - \frac{2AC}{4AC - B^{2}})

My voltage equation then becomes

V_{x}(t) = V_{x}(t = 0^{+})(1 + \frac{B\sqrt{B^{2}-4AC} + B^{2}}{2(4AC - B^{2})} - \frac{2AC}{4AC - B^{2}})e^{(\frac{-B + \sqrt{B^2 - 4AC}}{2A})t} - V_{x}(t = 0^{+})(\frac{B\sqrt{B^{2}-4AC} + B^{2}}{2(4AC - B^{2})} - \frac{2AC}{4AC - B^{2}})e^{(\frac{-B - \sqrt{B^2 - 4AC}}{2A})t}

Plugging back in the values for A, B, and C

V_{x}(t) = V_{x}(t = 0^{+})(1 + \frac{\frac{1}{R}\sqrt{(\frac{1}{R})^{2}-4(C_{1}+C_{2})\frac{1}{L}} + (\frac{1}{R})^{2}}{2(4(C_{1}+C_{2})\frac{1}{L} - (\frac{1}{R})^{2})} - \frac{2(C_{1}+C_{2})\frac{1}{L}}{4(C_{1}+C_{2})\frac{1}{L} - (\frac{1}{R})^{2}})e^{(\frac{-\frac{1}{R} + \sqrt{(\frac{1}{R})^2 - 4(C_{1}+C_{2})\frac{1}{L}}}{2(C_{1}+C_{2})})t} - V_{x}(t = 0^{+})(\frac{\frac{1}{R}\sqrt{(\frac{1}{R})^{2}-4(C_{1}+C_{2})\frac{1}{L}} + (\frac{1}{R})^{2}}{2(4(C_{1}+C_{2})\frac{1}{L} - (\frac{1}{R})^{2})} - \frac{2(C_{1}+C_{2})\frac{1}{L}}{4(C_{1}+C_{2})\frac{1}{L} - (\frac{1}{R})^{2}})e^{(\frac{-(\frac{1}{R}) - \sqrt{(\frac{1}{R})^2 - 4(C_{1}+C_{2})\frac{1}{L}}}{2(C_{1}+C_{2})})t}

That's kind of ugly. But anyways if my equation

V_{x}(t) = V_{x}(t = 0^{+})(1 + \frac{B\sqrt{B^{2}-4AC} + B^{2}}{2(4AC - B^{2})} - \frac{2AC}{4AC - B^{2}})e^{(\frac{-B + \sqrt{B^2 - 4AC}}{2A})t} - V_{x}(t = 0^{+})(\frac{B\sqrt{B^{2}-4AC} + B^{2}}{2(4AC - B^{2})} - \frac{2AC}{4AC - B^{2}})e^{(\frac{-B - \sqrt{B^2 - 4AC}}{2A})t}

is correct then I have developed a generic formula for circuits in which you can find a second order homogeneous equation. Does it look right? Thanks for your help!

When I plug the values pack in I get the following equation

V_{x}(t) = \frac{4}{49}e^{2000t}((7j - 401)e^{14000jt} - (7j - 499)e^{-14000jt})

Which doesn't seem right. Not sure what I did wrong though.

 

[GreenPrint]

ok so I guess I have found my mistake. At time t = 0^{+} the current that goes through the resistor where does it go? Does it go through the first capacitor or the second capacitor. I'm aware that it doesn't go through the inductor because initially the current through the inductor is zero and that current through an inductor can't change instantaneously. So I'm not really sure were this current goes and if someone could help in figuring this out than that would be great.

Originally in the post above I have that the current is initially zero which is wrong.

Thanks!

 

[NascentOxygen]

ok so I guess I have found my mistake. At time t = 0^{+} the current that goes through the resistor where does it go? Does it go through the first capacitor or the second capacitor. I'm aware that it doesn't go through the inductor because initially the current through the inductor is zero and that current through an inductor can't change instantaneously.

All that maths ... too much for me to check.

If we can neglect the inductor, then at t=0^- the top capacitor presents a voltage difference of 2V and the lower capacitor has a voltage of 8V. Then (without the inductor), after the switch is closed, the top capacitor eventually reaches a voltage of 0V, and the bottom capacitor carries all the supply voltage, 10V. So, with that scenario, you could say after t=0 the top capacitor discharges via R and the lower capacitor simultaneously charges via R. The voltage across one can't change without the voltage across the other also changing, in order to maintain their sum as 10V.

A changing voltage across a capacitor is an indication that current is flowing in it.

 

[rude man]

The most I'm willing to do is verify your answer. No way am I going to wade thru all that math. I can do it more quickly by going to Laplace transforms which I don't suppose you've had.

 

[GreenPrint]

So then if a capacitor can't change voltage instantaneous then the current through a capacitor is always zero at time t=0- and at time t=0+ because if the circuit has been sitting for a while the capacitors act as opens and there's current flowing through them at t=0-

So like I'm still confused then. If the voltage across a capacitor can't change then which way does the current from the resistor go?

Vs/R

at time t=0+ there's a current flowing through the resistor, but if it can't go through the inductor or either cap then were does it go?

And I completely forgot how to do laplace transforms =( in only a year sense taking differential equations.

 

[NascentOxygen]

I'm still confused

So it seems.

The voltage across a capacitor can change only after some time has passed to permit more charge to accumulate (or be removed) by current flowing into/outof its terminals. This is a process that takes time.

In general, in a CR circuit at t=0⁺ there can be current into/outof capacitors. For a brief time after t=0, current flows without significant change in capacitor voltage. Only after more charge has built up will the capacitor show an increase in voltage.

 

[GreenPrint]

So then were does the current go? Through capacitor one or two or both and how to find out which? Like I get that at time t=0- there is no current flowing either capacitor and that voltage across a capacitor can't change instantaneously, a change in voltage across a capacitor allows for a change in current. So based on this information the current can't go through either capacitor at time t=0+, so then were does it go? The current through an inductor can't change instantaneously and initially there is no current flowing through the inductor so it can't flow through the inductor at t=0+ either.

I think I'm still confused. I understood your post above but I still am unsure as to which capacitor it goes through.

 

[rude man]

At t(0-) C1 and C2 are charged up to 10V between them (so what would V_x(0+) be?).

At t(0+) there is instantaneous current flowing into C2, limited by the 50 ohm resistor. At the same time C1 begins to discharge, also into C2. Meanwhile current begins to build up in the inductor L.

At t = infinity C1 is charged to 10V, C2 has no charge, and current thru the inductor is 10V/50 ohms.

EDIT:
Armed with this knowledge you can hopefully write the pertinent DE's with the relevant initial conditions on C1, C2 and L.

 

[rude man]

I finally solved this thing & it was a chore. V_x(t) is a damped sinusoid. I have the results if you want to compare.

 

[GreenPrint]

well then V_x(0+) would be 10 volts as well

At time t=(0+) the current flows through C2, why not C1? So current begins to flow the C1 into C2 as well? C2 discharges to ground?

 

[NascentOxygen]

well then V_x(0+) would be 10 volts as well

Yes, it's the voltage across a capacitor, and that can't change in an instant.

At time t=(0+) the current flows through C2, why not C1? So current begins to flow the C1 into C2 as well? C2 discharges to ground?

As I said in my first reply here: So, with that scenario, you could say after t=0 the top capacitor discharges via R and the lower capacitor simultaneously charges via R. The voltage across one can't change without the voltage across the other also changing, in order to maintain their sum as 10V. [/color]

C2 discharges to ground

Because the inductor starts to affect things shortly after 0⁺ you can't say a capacitor discharges to ground because with the inclusion of inductance you should anticipate things will be oscillatory.

With the reasoning I outlined, current flows in both capacitors from the word go.

 

[rude man]

well then V_x(0+) would be 10 volts as well

No. V_x(0+) is the same as V_x(0-) which is not 10V.

At time t=(0+) the current flows through C2, why not C1? So current begins to flow the C1 into C2 as well? C2 discharges to ground?

Current flows thru C1 and C2 starting at t(0+). As I said, the movement of V_x(t) is sinusoidal with exponential decay. Eventually, C2 discharges to ground and V_x(∞) = 0.