Historical Context & Motivation
Long before modern computers, engineers and scientists needed ways to predict what happens when liquids mix. Imagine a large tank of salt water in a factory: fresh water flows in, the mixture stirs, and salty water flows out. How salty is the tank after one hour? After ten hours? These questions drove the development of mixing problems, one of the earliest and most practical applications of differential equations.
The mathematics behind these models grew alongside the industrial revolution, when chemical processing, water treatment, and pharmaceutical manufacturing all required precise control of concentrations. Today the same equations appear in environmental science (pollutant dispersion in lakes), medicine (drug dosing in the bloodstream), and food engineering (brine solutions). The core idea is beautifully simple: track the rate of change of a substance inside a container by comparing what flows in versus what flows out.
The central question these models answer is deceptively straightforward: if a substance enters a well-stirred container at one rate and leaves at another, how does the amount of that substance change over time? Answering this question precisely requires setting up and solving a first-order linear ODE — a skill that connects algebra, calculus, and real-world problem solving.
Core Principles & Definitions
Every salt-and-concentration model rests on a few foundational ideas. Before diving into equations, make sure you understand the physical setup and the vocabulary used to describe it. The typical scenario involves a tank that holds a liquid (often water) with some dissolved substance (often salt, but it could be a dye, a chemical, or a drug). Liquid flows into the tank at a known rate carrying a known concentration of the substance, and the well-stirred mixture flows out at a possibly different rate.
Rate In vs. Rate Out
Well-Stirred Assumption
Volume May Change
The Balance Equation
Visual Explanation — The Mixing Tank
Look carefully at the diagram. The key insight is that the rate of substance entering the tank is constant (because cin and rin are usually given constants), but the rate of substance leaving depends on Q(t), the unknown quantity we are trying to find. This is exactly what makes the problem a differential equation rather than a simple arithmetic calculation. Because the outflow carries the current concentration Q(t)/V(t), the rate out changes as Q changes, creating a feedback loop that the ODE captures.
Mathematical Framework
Let's build the equation step by step. Define Q(t) as the amount of salt (in kilograms, grams, or whatever unit you choose) at time t. The volume of liquid in the tank at time t is V(t). If the inflow and outflow rates are equal — say, both equal to r liters per minute — then V stays constant. Otherwise, V(t) = V₀ + (rin − rout)t, where V₀ is the initial volume.
Constant-Volume Case (rᵢₙ = rₒᵤₜ = r)
When the inflow and outflow rates are equal, the volume stays fixed at V₀. The equation simplifies to a first-order linear ODE with constant coefficients.
To solve, use an integrating factor μ(t) = ert/V₀. Multiply both sides of the ODE by μ(t), recognize the left side as d/dt[μQ], integrate both sides, and solve for Q(t). The result is a general solution.
Solution Behavior & Graphs
The general solution Q(t) = cinV₀ + (Q₀ − cinV₀)e−rt/V₀ exhibits two qualitatively different behaviors depending on whether the initial amount Q₀ is above or below the equilibrium value cinV₀. If Q₀ < cinV₀, the salt amount increases toward equilibrium (the exponential term is negative, so subtracting it adds to cinV₀). If Q₀ > cinV₀, the salt amount decreases toward equilibrium. Both curves approach the same steady-state value asymptotically.
Notice how both curves flatten out as they approach the equilibrium value. This flattening is the hallmark of exponential decay: the rate of change is proportional to the remaining gap between Q(t) and its equilibrium. When the gap is large, Q changes quickly; when the gap is small, Q barely changes. The time constant τ = V₀/r tells you the speed of convergence. A small tank with a high flow rate (small τ) reaches equilibrium quickly. A large tank with a slow flow (large τ) takes much longer.
| Scenario | Q₀ vs. Equilibrium | Curve Shape |
|---|---|---|
| Pure water → salty inflow | Q₀ < cinV₀ | Rising exponential toward equilibrium |
| Very salty → less salty inflow | Q₀ > cinV₀ | Falling exponential toward equilibrium |
| Flushing (pure water inflow, cin = 0) | Q₀ > 0 = equilibrium | Pure exponential decay: Q = Q₀e−rt/V₀ |
Worked Example
A tank initially holds 500 liters of water in which 10 kg of salt is dissolved. Brine containing 0.05 kg/L of salt flows in at 8 L/min, and the well-stirred mixture flows out at 8 L/min. Find Q(t), the amount of salt at time t, and determine the concentration after 100 minutes.
Strengths & Limitations
The constant-volume mixing model is elegantly simple, but like all models, it rests on assumptions that may or may not hold in a given situation. Understanding these strengths and limitations helps you decide when the model is trustworthy and when you need a more sophisticated approach.
| Strengths | Limitations |
|---|---|
| Produces an exact, closed-form solution — no numerical methods needed | Assumes perfectly uniform mixing at every instant, which is rarely true in practice |
| Clearly illustrates the balance-law approach (rate in − rate out) | Ignores chemical reactions — the dissolved substance does not interact with the solvent |
| Easily extended to variable-volume cases by adjusting V(t) | In variable-volume problems, the ODE can become harder to solve analytically |
| Applies to many real-world contexts: pollution, pharmacology, food science | Does not handle multiple interconnected tanks without extending to systems of ODEs |
Connection to Advanced Theory
The single-tank mixing problem is a gateway to much deeper mathematics. Once you are comfortable with this model, several natural extensions arise. What if the inflow rate and outflow rate are unequal? What if there are two or more tanks connected in series? What if the inflow concentration changes over time? Each of these variations builds on the same rate-in-minus-rate-out logic but demands more sophisticated solution techniques.
| Feature | This Lesson (Basic Model) | Advanced Extensions |
|---|---|---|
| Number of tanks | Single tank | Cascading or interconnected tanks (systems of ODEs) |
| Volume | Constant (rin = rout) | Variable — V(t) changes, ODE coefficients become functions of t |
| Inflow concentration | Constant cin | Time-dependent cin(t), requiring integrating factor with variable G(t) |
| Solution method | Integrating factor with constant coefficients | Laplace transforms, matrix exponentials, numerical methods |
| Applications | Single-tank industrial mixing, simple drug dosing | Multi-compartment pharmacokinetics, environmental engineering, reactor networks |
If you continue studying differential equations, you will encounter systems of first-order ODEs that model interconnected tanks — for example, two tanks where the outflow of the first is the inflow of the second. These systems use linear algebra (matrices and eigenvalues) to find solutions, connecting differential equations to another major branch of mathematics. The single-tank model you learned here is the essential building block for all of that work.
Practice Problems
Lesson Summary
Salt and concentration models use the balance law dQ/dt = (rate in) − (rate out) to track how the amount of a dissolved substance changes inside a well-stirred tank. In the constant-volume case (rin = rout), the ODE is first-order linear with constant coefficients, solvable via the integrating factor method. The solution is Q(t) = cinV₀ + (Q₀ − cinV₀)e−rt/V₀, which shows exponential convergence toward the equilibrium value cinV₀ at a rate governed by the time constant τ = V₀/r.
When the inflow and outflow rates differ, the volume V(t) changes over time, making the ODE coefficients time-dependent and the problem more challenging. These models apply broadly — from industrial mixing tanks to drug dosing in medicine to pollutant tracking in environmental science — making them one of the most versatile applications of first-order ODEs.