3.3 Ordinary Differential Equations & Initial-Value Problems

3.3 Ordinary Differential Equations & Initial-Value Problems

Ordinary Differential Equations (ODEs) form the dynamic mathematical backbone of electrical engineering. Transient phenomena in electrical circuits containing energy storage elements (inductors and capacitors), control systems, electrical machinery dynamics, and electromagnetic waves are modeled by ordinary differential equations.

---

1. Classification of Differential Equations

• Order: The order of the highest derivative present in the equation.
• Degree: The algebraic exponent of the highest-order derivative after removing fractional exponents.
• Linearity: An ODE is linear if the dependent variable $y$ and its derivatives appear only to the first power and are not multiplied together (e.g., $a_n(x) y^{(n)} + \dots + a_1(x) y' + a_0(x) y = g(x)$).
• Homogeneity: Homogeneous if $g(x) = 0$; non-homogeneous if $g(x) \neq 0$.

---

2. First-Order Differential Equations

1. Variable Separable ODEs

Form: $\frac{dy}{dx} = g(x) h(y) \implies \int \frac{dy}{h(y)} = \int g(x) , dx + C$

2. Exact Differential Equations

Form: $M(x,y),dx + N(x,y),dy = 0$

• Test for Exactness: $\frac{\partial M}{\partial y} = \frac{\partial N}{\partial x}$
• Solution is given by potential function $\Phi(x,y) = C$ where $\frac{\partial \Phi}{\partial x} = M$ and $\frac{\partial \Phi}{\partial y} = N$.
• If not exact, an Integrating Factor $\mu(x) = e^{\int \frac{M_y - N_x}{N} , dx}$ converts it to an exact ODE.

3. First-Order Linear ODEs

Standard Form: $\frac{dy}{dx} + P(x)y = Q(x)$

• Integrating Factor: $I(x) = e^{\int P(x) , dx}$
• General Solution: $y \cdot I(x) = \int Q(x) I(x) \, dx + C \implies y(x) = e^{-\int P dx} \left[ \int Q(x) e^{\int P dx} \, dx + C \right]$

4. Bernoulli Equations

Form: $\frac{dy}{dx} + P(x)y = Q(x)y^n \quad (n \neq 0, 1)$

• Substitution $u = y^{1-n}$ transforms the non-linear equation into a linear ODE in terms of $u(x)$: $\frac{du}{dx} + (1-n) P(x) u = (1-n) Q(x)$

---

3. Second-Order Linear Homogeneous ODEs with Constant Coefficients

Standard Form: $a \frac{d^2y}{dt^2} + b \frac{dy}{dt} + c y = 0$

Characteristic (Auxiliary) Equation

Substituting $y = e^{rt}$ yields: $a r^2 + b r + c = 0 \implies r_{1,2} = \frac{-b \pm \sqrt{b^2 - 4ac}}{2a}$

Depending on the discriminant $\Delta = b^2 - 4ac$, three distinct physical damping regimes arise:

Where $\alpha = -\frac{b}{2a}$ (attenuation constant) and $\beta = \frac{\sqrt{4ac - b^2}}{2a}$ (damped angular frequency).

---

4. Second-Order Non-Homogeneous ODEs

Form: $a y'' + b y' + c y = g(t)$ General Solution: $y(t) = y_h(t) + y_p(t)$ (Natural Transient + Forced Steady-State Solution).

Method of Undetermined Coefficients

Used when forcing function $g(t)$ consists of polynomials, exponentials, sines, or cosines:

Note: If any term in the trial $y_p(t)$ duplicates a term in $y_h(t)$, multiply the trial solution by $t$ (or $t^2$).

Variation of Parameters

Universal method applicable for any continuous forcing function $g(t)$: $y_p(t) = -y_1(t) \int \frac{y_2(t) g(t)}{W(t)} \, dt + y_2(t) \int \frac{y_1(t) g(t)}{W(t)} \, dt$ where the Wronskian is $W(t) = y_1 y_2' - y_1' y_2$.

---

5. Electrical Circuit Applications (Initial-Value Problems)

Series RL Circuit

Governed by Kirchhoff's Voltage Law (KVL): $L \frac{di}{dt} + R i = V(t)$ For a step DC voltage input $V(t) = V_0$ at $t = 0$ with $i(0) = 0$: $i(t) = \frac{V_0}{R} \left(1 - e^{-t/\tau}\right) \quad \text{where } \tau = \frac{L}{R} \text{ (Time Constant)}$

Series RC Circuit

Governed by KVL in terms of charge $q(t)$: $R \frac{dq}{dt} + \frac{1}{C} q = V(t) \implies R \frac{di}{dt} + \frac{1}{C} i = \frac{dV}{dt}$ For a step DC voltage charging a capacitor with initial charge $q(0) = 0$: $i(t) = \frac{V_0}{R} e^{-t/\tau} \quad \text{where } \tau = R C \text{ (Time Constant)}$

Series RLC Circuit

KVL formulation: $L \frac{d^2q}{dt^2} + R \frac{dq}{dt} + \frac{1}{C} q = E(t) \implies L \frac{d^2i}{dt^2} + R \frac{di}{dt} + \frac{1}{C} i = \frac{dE}{dt}$ Key circuit metrics:

• Undamped natural frequency: $\omega_0 = \frac{1}{\sqrt{L C}}$
• Damping ratio: $\zeta = \frac{R}{2} \sqrt{\frac{C}{L}}$
  • $\zeta > 1 \implies$ Overdamped
  • $\zeta = 1 \implies$ Critically Damped
  • $\zeta < 1 \implies$ Underdamped

---

Step-by-Step Worked Numerical Examples

Example 3.3.1: First-Order RL Circuit Transient IVP

Problem: A series RL circuit has resistance $R = 10\ \Omega$, inductance $L = 2\text{ H}$, and a DC voltage source $V = 50\text{ V}$ connected at $t = 0$. Given initial current $i(0) = 0\text{ A}$, find $i(t)$ and calculate the current at $t = 0.4\text{ seconds}$.

Solution:

1. Formulate ODE: $2 \frac{di}{dt} + 10 i = 50 \implies \frac{di}{dt} + 5 i = 25$
2. Find Integrating Factor $I(t)$: $P(t) = 5 \implies I(t) = e^{\int 5 \, dt} = e^{5t}$
3. Multiply ODE by $I(t)$ and integrate: $\frac{d}{dt}\left(i e^{5t}\right) = 25 e^{5t} \implies i e^{5t} = \int 25 e^{5t} \, dt = 5 e^{5t} + C$ $i(t) = 5 + C e^{-5t}$
4. Apply initial condition $i(0) = 0$: $0 = 5 + C(1) \implies C = -5$ $i(t) = 5 \left(1 - e^{-5t}\right)\text{ A}$
5. Evaluate at $t = 0.4\text{ s}$: $i(0.4) = 5 \left(1 - e^{-5(0.4)}\right) = 5 \left(1 - e^{-2}\right) = 5 (1 - 0.135335) = 5(0.864665) = 4.3233\text{ A}$

---

Example 3.3.2: RLC Circuit Transient Damping Regime Analysis

Problem: A series RLC circuit has $L = 1\text{ H}$ and $C = 0.25\text{ F}$. Determine the response regime (overdamped, critically damped, or underdamped) and the natural response $q_h(t)$ for $R = 5\ \Omega$.

Solution:

1. Formulate characteristic equation: $L r^2 + R r + \frac{1}{C} = 0 \implies 1 r^2 + 5 r + \frac{1}{0.25} = 0 \implies r^2 + 5r + 4 = 0$
2. Find roots of characteristic equation: $(r + 1)(r + 4) = 0 \implies r_1 = -1, \quad r_2 = -4$
3. Identify damping regime: Since $r_1, r_2$ are real and distinct (discriminant $\Delta = 5^2 - 4(1)(4) = 9 > 0$), the circuit response is Overdamped.
4. Write natural transient charge solution: $q_h(t) = C_1 e^{-t} + C_2 e^{-4t}$

---

Example 3.3.3: Second-Order Non-Homogeneous Resonant System

Problem: Solve the initial value problem $y'' + 4y = 12 \cos(2t)$ with $y(0) = 0$ and $y'(0) = 0$.

Solution:

1. Homogeneous solution $y_h(t)$: $r^2 + 4 = 0 \implies r = \pm j2 \implies y_h(t) = C_1 \cos(2t) + C_2 \sin(2t)$
2. Formulate trial particular solution $y_p(t)$: Because forcing term $12 \cos(2t)$ duplicates terms in $y_h(t)$, multiply trial by $t$: $y_p(t) = t (A \cos 2t + B \sin 2t)$
3. Differentiate $y_p(t)$: $y_p' = (A \cos 2t + B \sin 2t) + t (-2A \sin 2t + 2B \cos 2t)$ $y_p'' = -4A \sin 2t + 4B \cos 2t + t (-4A \cos 2t - 4B \sin 2t)$
4. Substitute into ODE $y_p'' + 4y_p = 12 \cos 2t$: $-4A \sin 2t + 4B \cos 2t - 4 y_p + 4 y_p = 12 \cos 2t$ $-4A \sin 2t + 4B \cos 2t = 12 \cos 2t$ Equating coefficients: $-4A = 0 \implies A = 0$, and $4B = 12 \implies B = 3$. $y_p(t) = 3t \sin(2t)$
5. Apply initial conditions to general solution $y(t) = C_1 \cos 2t + C_2 \sin 2t + 3t \sin 2t$:
   • $y(0) = C_1 + 0 + 0 = 0 \implies C_1 = 0$
   • $y'(t) = 2C_2 \cos 2t + 3 \sin 2t + 6t \cos 2t$
   • $y'(0) = 2C_2 + 0 + 0 = 0 \implies C_2 = 0$
6. Final Solution: $y(t) = 3t \sin(2t)$