Lecture 33: Separable Differential Equations and Applications

First Order Differential Equations: Separable Equations

  • Lecture focuses on separable differential equations and their engineering applications.

First Order Differential Equations

  • A first-order differential equation can be written as: dydx=f(x,y)\frac{dy}{dx} = f(x, y), where f(x,y)f(x, y) is a function depending on both xx and yy.
  • yy is the unknown solution, implicit in the equation.
  • The set of all solutions forms the general solution, a family of functions.
  • Imposing the condition y(x<em>0)=y</em>0y(x<em>0) = y</em>0 turns the differential equation into an initial value problem (IVP).

Separable First Order Differential Equations

  • Separable differential equations have the form: dydx=f(x)g(y)\frac{dy}{dx} = f(x)g(y).
  • Separation of variables leads to: dyg(y)=f(x)dx\frac{dy}{g(y)} = f(x)dx.
  • Solving involves integrating both sides: 1g(y)dy=f(x)dx+C\int \frac{1}{g(y)} dy = \int f(x) dx + C.

Examples of Separable Differential Equations

Example 1
  • Given differential equation: dydx=k(yy<em>0)\frac{dy}{dx} = k(y - y<em>0), where kk and y</em>0y</em>0 are constants.
  • Separation gives: dyyy0=kdx\frac{dy}{y - y_0} = k dx.
  • Integrating both sides: dyyy<em>0=kdx\int \frac{dy}{y - y<em>0} = \int k dx yields lnyy</em>0=kx+c\ln|y - y</em>0| = kx + c.
  • Taking exponentials: yy0=ekx+c=ekxecy - y_0 = e^{kx + c} = e^{kx}e^c.
  • Rename ec=Ae^c = A: y=y0+Aekxy = y_0 + Ae^{kx}, where AA is an arbitrary constant.
Example 2 (Corrected)
  • Differential equation: dydx=xy\frac{dy}{dx} = -\frac{x}{y}.
  • Separating variables: ydy=xdxy dy = -x dx.
  • Integrating: ydy=xdx+C\int y dy = -\int x dx + C leads to y22=x22+C\frac{y^2}{2} = -\frac{x^2}{2} + C.
  • Which implies y2+x2=2Cy^2 + x^2 = 2C.
  • Impose initial condition y(0)=r<em>0y(0) = r<em>0 (with r0 > 0).
  • Since y=2Cx2y = \sqrt{2C - x^2}, then y(0)=2C=r<em>0y(0) = \sqrt{2C} = r<em>0, so r</em>02=2Cr</em>0^2 = 2C.
  • The solution to the IVP is y=r02x2y = \sqrt{r_0^2 - x^2}.
    • Solution holds for positive r<em>0r<em>0, a similar analysis can be performed for negative r</em>0r</em>0.
Example 3
  • Differential equation: dydx=y(x+1)\frac{dy}{dx} = y(x + 1).
  • Separating variables: dyy=(x+1)dx\frac{dy}{y} = (x + 1) dx.
  • Integrating: dyy=(x+1)dx\int \frac{dy}{y} = \int (x + 1) dx yields lny=x22+x+C\ln|y| = \frac{x^2}{2} + x + C.
  • Taking exponentials: y=ex22+x+C=ex22+xeCy = e^{\frac{x^2}{2} + x + C} = e^{\frac{x^2}{2} + x} e^C.
  • Let A=eCA = e^C: y=Aex22+xy = A e^{\frac{x^2}{2} + x}.

Application 1: R-C Circuit

  • Consider a circuit with resistance RR, capacitance CC, and voltage source VV.
  • I(t)I(t) is the current, and Q(t)Q(t) is the charge on the capacitor.
  • Experimental law: IR+QC=VIR + \frac{Q}{C} = V.
  • Also, the current is the derivative of the charge with respect to time: I=dQdtI = \frac{dQ}{dt}.
  • Substituting, we get: RdQdt+QC=VR \frac{dQ}{dt} + \frac{Q}{C} = V, a first-order separable differential equation.
  • Given: R=50ΩR = 50 \Omega, C=2C = 2 Farads, V=12V = 12 Volts.
  • The equation becomes: 50dQdt+Q2=1250 \frac{dQ}{dt} + \frac{Q}{2} = 12.
  • Solving for dQdt\frac{dQ}{dt}: 50dQdt=12Q2=24Q250 \frac{dQ}{dt} = 12 - \frac{Q}{2} = \frac{24 - Q}{2}.
  • Separating variables: dQ24Q=dt100\frac{dQ}{24 - Q} = \frac{dt}{100}.
  • Integrating: dQ24Q=dt100\int \frac{dQ}{24 - Q} = \int \frac{dt}{100} gives ln24Q=t100+C- \ln|24 - Q| = \frac{t}{100} + C.
  • Multiplying by -1, we have ln24Q=t100C\ln|24 - Q| = -\frac{t}{100} - C.
  • Exponential on both sides: 24 - Q = et100Ce^{-\frac{t}{100} - C}.
  • Which becomes: 24Q=et100eC24 - Q = e^{-\frac{t}{100}}e^{-C}.
  • Which then becomes: Q=24et100eCQ = 24 - e^{-\frac{t}{100}}e^{-C}.
  • Initial condition: Q(0)=0Q(0) = 0.
  • This means 0=24e0eC0 = 24 - e^{0}e^{-C} which means that eC=24e^{-C} = 24.
  • Therefore, Q(t)=2424et100Q(t) = 24 - 24e^{-\frac{t}{100}}.
  • Final charge: lim<em>tQ(t)=lim</em>t(2424et100)=24\lim<em>{t \to \infty} Q(t) = \lim</em>{t \to \infty} (24 - 24e^{-\frac{t}{100}}) = 24 Coulombs.
  • Time to reach 90% of total charge: 0.9×24=21.60.9 \times 24 = 21.6.
  • Solve 21.6=2424et10021.6 = 24 - 24e^{-\frac{t}{100}} for tt.
  • 21.624=24et10021.6 - 24 = -24e^{-\frac{t}{100}}
  • 2.4=24et100-2.4 = -24e^{-\frac{t}{100}}.
  • 0.1=et1000.1 = e^{-\frac{t}{100}}.
  • ln(0.1)=t100\ln(0.1) = -\frac{t}{100}.
  • t=100ln(0.1)=100ln(101)=100ln(10)230.3t = -100 \ln(0.1) = -100 \ln(10^{-1}) = 100\ln(10) \approx 230.3 seconds.
  • Absolute value clarification: Removing the absolute value implies 24 - Q > 0, meaning Q < 24, consistent with the physical constraint.

Application 2: Dilution of Salt Solution

  • Tank with 600 liters of water and 200 grams of salt at t=0t = 0.
  • Brine solution pumped in at 3 liters/minute with 4 grams of salt/liter.
  • Mixed solution pumped out at 3 liters/minute.
  • Problem: Find the salt concentration of the pumped-out solution as a function of time.
  • Definitions:
    • S(t)S(t) = Amount of salt in the tank at time tt.
    • dSdt\frac{dS}{dt} = Rate of change of salt in the tank.
  • Rate of change equals to Rate of salt entering - Rate of salt leaving
  • Rate of salt entering: 3 liters/minute * 4 grams/liter = 12 grams/minute.
  • Rate of salt leaving: (S(t)600)3=3S(t)600(\frac{S(t)}{600}) \cdot 3 = \frac{3S(t)}{600} grams/minute.
  • Differential equation: dSdt=123S600\frac{dS}{dt} = 12 - \frac{3S}{600}.
  • Initial condition: S(0)=200S(0) = 200 grams.
  • To solve the IVP all that is required is to:
    1. Separate variables
    2. Integrate both sides to get general solution
    3. Apply initial conditions to get specific solution
  • Amount of salt in the tank at time t (solution after process), S(t)=24002200et/200S(t)=2400-2200e^{-t/200}
  • Concentration of salt leaving the tank: S(t)600=24002200et200600=4113et200\frac{S(t)}{600} = \frac{2400 - 2200e^{-\frac{t}{200}}}{600} = 4 - \frac{11}{3}e^{-\frac{t}{200}}.
  • Asymptotic behavior can be seen over time, concentration gets closer to 4.