DIFFERENTIAL EQUATIONS • APPLICATIONS OF FIRST-ORDER ODES

Salt & Concentration Models

Use first-order differential equations to predict how the concentration of a dissolved substance changes over time in a mixing tank.

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.

1690s
Leibniz & Newton Formalize Calculus
The co-invention of calculus provided the language of derivatives and integrals needed to describe rates of change, laying the groundwork for all differential equation models.
1760s
Euler's Systematic ODE Methods
Leonhard Euler developed systematic techniques for solving first-order ordinary differential equations, including the integrating factor method central to mixing problems.
1800s
Industrial Chemical Engineering
Factories began using continuous-flow reactors and mixing tanks. Engineers applied first-order ODEs to predict chemical concentrations in large-scale production.
1950s
Environmental & Biomedical Modeling
The same differential equations were adapted to model pollutant concentrations in lakes and rivers, and drug levels in the human body — fields now called compartmental analysis.

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.

1

Rate In vs. Rate Out

The amount of substance in the tank changes because of two competing flows. The rate in is the inflow rate times the inflow concentration. The rate out depends on the current concentration inside the tank times the outflow rate.
2

Well-Stirred Assumption

We assume the tank is stirred so thoroughly that the concentration is uniform throughout at every moment. This means the concentration of the outflow equals the concentration everywhere inside the tank.
3

Volume May Change

If the inflow rate and outflow rate are not equal, the total volume of liquid in the tank changes over time. This affects the concentration calculation because concentration equals amount of substance divided by volume.
4

The Balance Equation

The fundamental equation is: dQ/dt = (rate of substance in) − (rate of substance out), where Q(t) is the total amount of dissolved substance at time t. This is a first-order ODE.
KEY TAKEAWAY
Think of the tank like a bathtub with the faucet running and the drain open. If you pour bubble bath into the incoming water, the amount of bubbles in the tub changes over time — it depends on how fast bubbly water enters and how fast soapy water drains out. The differential equation simply keeps a running tally of what comes in minus what goes out at every instant.

Visual Explanation — The Mixing Tank

The diagram shows a well-stirred mixing tank. The blue arrow on the left represents inflow carrying dissolved salt at concentration cin at flow rate rin. The red arrow on the right shows outflow at rate rout, carrying the well-mixed solution at concentration Q(t)/V(t). The boxed equation at the bottom is the fundamental balance law.

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.

BALANCE LAW
dQ/dt = rᵢₙ × cᵢₙ − rₒᵤₜ × Q(t) / V(t)
Q(t) = amount of salt at time t; rin = inflow rate (L/min); cin = inflow concentration (kg/L); rout = outflow rate (L/min); V(t) = volume of liquid at time t.

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.

CONSTANT-VOLUME ODE
dQ/dt + (r / V₀) × Q = r × cᵢₙ
This is in standard linear form dQ/dt + P(t)Q = G(t), where P = r/V₀ and G = r × cin. Both P and G are constants here.

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.

GENERAL SOLUTION (CONSTANT VOLUME)
Q(t) = cᵢₙ × V₀ + (Q₀ − cᵢₙ × V₀) × e^(−rt / V₀)
Q₀ = initial amount of salt in the tank. As t → ∞, Q(t) → cin × V₀, meaning the tank eventually reaches the inflow concentration. The exponential term controls how quickly it gets there.
💡 Why the Exponential?
The exponential decay factor e−rt/V₀ tells you the system has a characteristic time constant τ = V₀/r. After one time constant (t = τ), the difference between Q(t) and its equilibrium value has shrunk by a factor of 1/e ≈ 0.37. After three time constants, the system is about 95% of the way to equilibrium.

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.

Two solution curves for the same tank with equilibrium at Q = 20 kg. The cyan curve starts with pure water (Q₀ = 0) and rises toward equilibrium. The violet curve starts with 40 kg of salt and falls toward the same equilibrium. Both approach the dashed amber line. The vertical dashed line marks one time constant τ = V₀/r.

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.

Three common mixing scenarios and their qualitative behaviors
ScenarioQ₀ vs. EquilibriumCurve Shape
Pure water → salty inflowQ₀ < cinV₀Rising exponential toward equilibrium
Very salty → less salty inflowQ₀ > cinV₀Falling exponential toward equilibrium
Flushing (pure water inflow, cin = 0)Q₀ > 0 = equilibriumPure 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.

Mixing Tank — Constant Volume
1
Step 1 — Identify Given ValuesV₀ = 500 L (constant, since rin = rout = 8 L/min). The inflow concentration is cin = 0.05 kg/L. The initial salt is Q₀ = 10 kg.
r = 8, V₀ = 500, cin = 0.05, Q₀ = 10
2
Step 2 — Write the ODEUsing the balance law: dQ/dt = r × cin − r × Q/V₀ = 8 × 0.05 − 8 × Q/500 = 0.4 − 0.016Q. In standard form: dQ/dt + 0.016Q = 0.4.
dQ/dt + 0.016Q = 0.4
3
Step 3 — Find the Equilibrium ValueAt equilibrium, dQ/dt = 0, so 0.016Qeq = 0.4, giving Qeq = 25 kg. This makes sense: cin × V₀ = 0.05 × 500 = 25 kg.
Qeq = 25 kg
4
Step 4 — Apply the General SolutionQ(t) = cinV₀ + (Q₀ − cinV₀)e−rt/V₀ = 25 + (10 − 25)e−0.016t = 25 − 15e−0.016t.
Q(t) = 25 − 15e−0.016t
5
Step 5 — Evaluate at t = 100 minQ(100) = 25 − 15e−1.6 = 25 − 15 × 0.2019 ≈ 25 − 3.03 = 21.97 kg. The concentration is Q(100)/V₀ = 21.97/500 ≈ 0.0439 kg/L.
Q(100) ≈ 22.0 kg, concentration ≈ 0.044 kg/L
6
Step 6 — Interpret the ResultThe tank started with 10 kg of salt (concentration 0.02 kg/L) and is heading toward the equilibrium of 25 kg (concentration 0.05 kg/L). After 100 minutes — which is about 1.6 time constants (τ = 500/8 = 62.5 min) — the salt amount is approximately 88% of the way from Q₀ to Qeq. This aligns with the exponential convergence pattern we saw in the graphs.

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 and limitations of the single-tank mixing model
StrengthsLimitations
Produces an exact, closed-form solution — no numerical methods neededAssumes 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 scienceDoes not handle multiple interconnected tanks without extending to systems of ODEs
KEY TAKEAWAY
The mixing model is like a weather forecast: incredibly useful for planning and prediction, but it simplifies a complex reality. A weather forecast assumes certain atmospheric patterns hold; the mixing model assumes the tank is perfectly stirred. Both give great results most of the time, but you should always check whether the assumptions match your actual situation before trusting the answer blindly.

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.

How the basic mixing model extends into more advanced territory
FeatureThis Lesson (Basic Model)Advanced Extensions
Number of tanksSingle tankCascading or interconnected tanks (systems of ODEs)
VolumeConstant (rin = rout)Variable — V(t) changes, ODE coefficients become functions of t
Inflow concentrationConstant cinTime-dependent cin(t), requiring integrating factor with variable G(t)
Solution methodIntegrating factor with constant coefficientsLaplace transforms, matrix exponentials, numerical methods
ApplicationsSingle-tank industrial mixing, simple drug dosingMulti-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

PROBLEM 1CONCEPTUAL
A tank contains 200 L of pure water. Brine with concentration 0.1 kg/L flows in at 5 L/min, and the well-stirred mixture flows out at 5 L/min. Without solving the ODE, predict the long-run (steady-state) amount of salt in the tank and explain your reasoning.
PROBLEM 2BASIC CALCULATION
A 100-liter tank initially holds 5 kg of dissolved salt. Pure water (cin = 0) flows in at 4 L/min and the well-stirred solution flows out at 4 L/min. Write the ODE, solve it, and find Q(t).
PROBLEM 3INTERMEDIATE
A 300-liter tank initially contains 15 kg of salt dissolved in water. Brine with concentration 0.04 kg/L flows in at 6 L/min, and the mixture flows out at 6 L/min. (a) Find Q(t). (b) How long does it take for the salt amount to reach 10 kg?
PROBLEM 4APPLIED
A hospital IV drip delivers saline solution containing 0.009 kg/L of salt into a patient's bloodstream at 0.002 L/min. The body eliminates the saline mixture at the same rate. Model the blood volume as a fixed 5-liter compartment initially containing 0.04 kg of salt. Find the salt amount Q(t) in the blood and determine how long it takes for the blood concentration to reach 0.008 kg/L.
PROBLEM 5CRITICAL THINKING
A 400-liter tank initially holds 20 kg of salt. Brine at 0.06 kg/L flows in at 10 L/min, but the outflow rate is only 6 L/min, so the tank is filling up. (a) Write an expression for V(t). (b) Set up the ODE for Q(t). (c) Explain why this ODE is harder to solve than the constant-volume case and describe what additional considerations arise (e.g., when does the tank overflow?).

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.

Varsity Tutors • Differential Equations • Salt & Concentration Models