Applications of First-Order DEs

Learning Objectives

  • Formulate and solve first-order differential equations for exponential growth and decay.
  • Apply Newton's Law of Cooling to find temperature changes over time.
  • Model and solve mixing problems involving inflow and outflow rates.
  • Analyze simple RL and RC electrical circuits using Kirchhoff's Voltage Law.
  • Model falling bodies experiencing air resistance.
  • Determine the orthogonal trajectories for a given family of curves.

First-order differential equations are the backbone of mathematical modeling in science and engineering. They describe phenomena where the rate of change of a quantity is proportional to the quantity itself or some function of it.

Mathematical Modeling

Many real-world phenomena can be modeled using first-order differential equations. The core idea is to translate the physical principles governing the rate of change of a system into a mathematical equation. These models allow us to predict future states of the system based on its initial conditions.

Exponential Growth and Decay

Models population growth, radioactive decay, and continuously compounded interest.

dPdt=kPβ€…β€ŠβŸΉβ€…β€ŠP(t)=P0ekt\frac{dP}{dt} = kP \implies P(t) = P_0 e^{kt}

Variables

SymbolDescriptionUnit
P(t)P(t)Quantity at time ttvaries
ttTimevaries
kkProportionality constant (growth if k>0k > 0, decay if k<0k < 0)1/time
P0P_0Initial quantity at t=0t = 0varies

Newton's Law of Cooling

Models the temperature change of an object in a surrounding medium.

dTdt=βˆ’k(Tβˆ’Ts)\frac{dT}{dt} = -k(T - T_s)

Variables

SymbolDescriptionUnit
T(t)T(t)Temperature of the object at time tt∘C or βˆ˜F^\circ C \text{ or } ^\circ F
ttTimes,min, or hs, min, \text{ or } h
TsT_sConstant ambient temperature of the surroundings∘C or βˆ˜F^\circ C \text{ or } ^\circ F
kkCooling constant (k>0k > 0)1/time

Mixing Problems

Models the amount of solute in a tank with inflow and outflow.

dAdt=RinCinβˆ’RoutCout\frac{dA}{dt} = R_{in}C_{in} - R_{out}C_{out}

Variables

SymbolDescriptionUnit
A(t)A(t)Amount of solute in the tank at time ttkg or lbkg \text{ or } lb
ttTimes or mins \text{ or } min
RinR_{in}Inflow rateL/s or gal/minL/s \text{ or } gal/min
CinC_{in}Inflow concentrationkg/L or lb/galkg/L \text{ or } lb/gal
RoutR_{out}Outflow rateL/s or gal/minL/s \text{ or } gal/min
CoutC_{out}Outflow concentration (A(t)/V(t)A(t) / V(t))kg/L or lb/galkg/L \text{ or } lb/gal

Interactive Simulation

Visualize how the amount of solute in a tank changes over time with inflow and outflow. Interact with the mixing tank simulation below to see first-order ODEs in application.

Mixing Tank Model

dAdt=RinCinβˆ’RoutAV(t)\frac{dA}{dt} = R_{in} C_{in} - R_{out} \frac{A}{V(t)}
3 L/min
2 kg/L
3 L/min
Salt Amount10.0 kg
Tank Volume100.0 L
A(t) vs t

Electrical Circuits (RL and RC Series Circuits)

First-order linear differential equations govern simple electrical circuits with constant voltage sources (like a battery, EE) connected in series with resistors (RR) and either inductors (LL) or capacitors (CC).

According to Kirchhoff's Voltage Law, the sum of voltage drops around a closed loop equals the applied voltage E(t)E(t). These circuits can be modeled with standard linear DEs solvable with an Integrating Factor.

RL Series Circuit

Models the current in a circuit with a resistor and an inductor.

Ldidt+Ri=E(t)L\frac{di}{dt} + Ri = E(t)

Variables

SymbolDescriptionUnit
LLInductanceH
i(t)i(t)Current at time ttA
RRResistanceΞ©\Omega
E(t)E(t)Applied voltageV
ttTimes

RC Series Circuit

Models the charge on a capacitor in a circuit with a resistor and a capacitor.

Rdqdt+1Cq=E(t)R\frac{dq}{dt} + \frac{1}{C}q = E(t)

Variables

SymbolDescriptionUnit
RRResistanceΞ©\Omega
q(t)q(t)Charge at time ttC
CCCapacitanceF
E(t)E(t)Applied voltageV
ttTimes

Falling Bodies with Air Resistance

In classical mechanics, a falling body experiences a downward gravitational force and an upward air resistance force. When air resistance is proportional to velocity, this system can be modeled using a first-order differential equation.

Terminal Velocity Model

Applies Newton's Second Law to a falling body experiencing air resistance.

mdvdt=mgβˆ’kvm\frac{dv}{dt} = mg - kv

Variables

SymbolDescriptionUnit
mmMass of the bodykg
v(t)v(t)Velocity at time ttm/s
ttTimes
ggAcceleration due to gravitym/s2m/s^2
kkAir resistance proportionality constant (k>0k > 0)kg/s

Terminal Velocity Limit

The terminal velocity model is a separable and linear differential equation. As tβ†’βˆžt \rightarrow \infty, the acceleration dvdt\frac{dv}{dt} approaches 00, and the velocity approaches a constant limit known as the terminal velocity, given by vt=mgkv_t = \frac{mg}{k}.

Orthogonal Trajectories

An orthogonal trajectory is a curve that intersects every member of a given family of curves at right angles. It has significant applications in physics, such as finding electric field lines that are perpendicular to equipotential curves, or heat flow lines perpendicular to isotherms.

Finding Orthogonal Trajectories

  1. Find the DE: Differentiate the given family of curves to find its differential equation dydx=f(x,y)\frac{dy}{dx} = f(x,y). You may need to eliminate the arbitrary constant parameter first.
  2. Take the Negative Reciprocal: Replace dydx\frac{dy}{dx} with its negative reciprocal βˆ’dxdy-\frac{dx}{dy} (since the slopes of perpendicular lines multiply to βˆ’1-1).
  3. Solve the New DE: Solve the new differential equation to find the family of orthogonal trajectories.
Key Takeaways
  • Growth/Decay: Use P(t)=P0ektP(t) = P_0 e^{kt}. kk is positive for growth, negative for decay.
  • Newton's Law of Cooling: Rate of cooling is proportional to the temperature difference (Tβˆ’Ts)(T - T_s).
  • Mixing Problems: Balance equation is RateΒ ofΒ Change=RateΒ Inβˆ’RateΒ Out\text{Rate of Change} = \text{Rate In} - \text{Rate Out}. Ensure units match!
  • RL/RC Circuits: Modeled by linear DEs derived from Kirchhoff's Voltage Law.
  • Falling Bodies: Include air resistance proportional to velocity, leading to a terminal velocity.
  • Orthogonal Trajectories: Found by solving the DE with the negative reciprocal slope (mβŠ₯=βˆ’1/mm_{\perp} = -1/m).