6.7 First-Order Differential Equations: Linear & Non-Linear

Why Differential Equations Matter in Paper-II

A differential equation is a relation containing derivatives. Its order is the order of the highest derivative, and its degree is the power of that highest derivative once the equation is cleared of radicals and fractions. Almost every dynamic process examined later in Paper-II is described by one: the cooling of a forged component obeys a first-order equation, damped vibration of a shaft obeys a second-order one, and the transient temperature of a lumped body in heat transfer is a first-order decay.

Because CIL awards one mark per item with no negative marking, the productive skill here is recognising which of four standard types you are looking at, since each has a fixed recipe.

Type 1: Variable Separable

If the equation can be rearranged so that all $y$ terms sit with $dy$ and all $x$ terms with $dx$:

$f(y)\,dy = g(x)\,dx \quad \Longrightarrow \quad \int f(y)\,dy = \int g(x)\,dx + C$

Always test this first — it requires no machinery.

Worked example — Newton's law of cooling. A casting at $T_0$ cools in surroundings at $T_s$ according to

$\frac{dT}{dt} = -k(T - T_s)$

Separating gives $\dfrac{dT}{T - T_s} = -k,dt$, so $\ln(T - T_s) = -kt + C$, and applying the initial condition $T = T_0$ at $t = 0$ yields

$T = T_s + (T_0 - T_s)e^{-kt}$

This exponential-decay form reappears verbatim as the lumped-capacitance solution in the Heat Transfer chapter.

Type 2: Homogeneous Equations

If $\dfrac{dy}{dx} = F!\left(\dfrac{y}{x}\right)$, that is, every term has the same total degree in $x$ and $y$, substitute

$y = vx \quad \Longrightarrow \quad \frac{dy}{dx} = v + x\frac{dv}{dx}$

The substitution always converts a homogeneous equation into a variable-separable one in $v$ and $x$.

Type 3: Linear First-Order Equations

The standard form is

$\frac{dy}{dx} + P(x)\,y = Q(x)$

where $P$ and $Q$ are functions of $x$ alone. Multiply through by the integrating factor

$\text{IF} = e^{\int P\,dx}$

which makes the left side an exact derivative, giving the solution

$y \cdot \text{IF} = \int Q \cdot \text{IF}\,dx + C$

Worked example. Solve $\dfrac{dy}{dx} + \dfrac{y}{x} = x^2$.

Here $P = 1/x$, so $\text{IF} = e^{\int dx/x} = e^{\ln x} = x$. Then $yx = \int x^3,dx = \dfrac{x^4}{4} + C$, giving $y = \dfrac{x^3}{4} + \dfrac{C}{x}$.

If instead the roles of $x$ and $y$ are reversed — the equation is linear in $x$ with $y$ as the independent variable — write it as $\dfrac{dx}{dy} + P(y)x = Q(y)$ and use the same recipe.

Type 4: Exact Differential Equations

An equation written as

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

is exact when

$\frac{\partial M}{\partial y} = \frac{\partial N}{\partial x}$

This is the same equality-of-mixed-partials condition met in the previous section, and it means the left side is the total differential of some function $u(x,y)$. The solution is $u = C$, obtained as

$\int M\,dx \;(\text{treating } y \text{ constant}) + \int \left(\text{terms of } N \text{ free of } x\right) dy = C$

When the exactness test fails, an integrating factor can sometimes restore it. Two standard cases:

• If $\dfrac{1}{N}\left(\dfrac{\partial M}{\partial y} - \dfrac{\partial N}{\partial x}\right)$ is a function of $x$ alone, call it $f(x)$; then $\text{IF} = e^{\int f,dx}$.
• If $\dfrac{1}{M}\left(\dfrac{\partial N}{\partial x} - \dfrac{\partial M}{\partial y}\right)$ is a function of $y$ alone, call it $g(y)$; then $\text{IF} = e^{\int g,dy}$.

Bernoulli's Equation

The non-linear form named in the syllabus is

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

Divide throughout by $y^n$ and substitute $v = y^{1-n}$. Since $\dfrac{dv}{dx} = (1-n)y^{-n}\dfrac{dy}{dx}$, the equation becomes

$\frac{dv}{dx} + (1-n)Pv = (1-n)Q$

which is linear in $v$ and solved by the integrating factor of Type 3.

Classification Flowchart

Application: Draining a Mine Sump

A sump of uniform cross-section $A$ drains through an orifice of area $a$. Torricelli's law gives the outflow velocity as $\sqrt{2gh}$, so

$A\frac{dh}{dt} = -C_d\,a\sqrt{2gh}$

This is variable-separable. Writing $h^{-1/2}dh = -\dfrac{C_d a\sqrt{2g}}{A}dt$ and integrating from $h_1$ to $h_2$ gives the emptying time

$t = \frac{2A\left(\sqrt{h_1} - \sqrt{h_2}\right)}{C_d\,a\sqrt{2g}}$

The same result is derived again from a control-volume argument in the Fluid Mechanics chapter, which is a useful cross-check that the calculus and the physics agree.