System Dynamics · Module 4 of 10
Fluid and Thermal System Modeling
A draining tank and a cooling part obey the same first-order law. Resistance opposes flow, capacitance stores it, and their product is a time constant that tells you how fast the system settles.
Readiness check
This module extends modeling to fluids and heat. Tick only what you can do closed-notes.
- Recall a first-order time constant τ and the 63 percent rule.
- Write a balance: rate in minus rate out equals rate of storage.
- Recall that flow is driven by a pressure or temperature difference.
- Compute a product of resistance and capacitance.
- Recall specific heat as energy per unit mass per degree.
The core idea
Fluid and thermal systems are first-order: a resistance opposes flow, a capacitance stores the quantity, and a balance gives a single exponential. The time constant is the product of resistance and capacitance, in either domain.
fluid: q = Δh/Rf, Cf dh/dt = qin − qoutthermal: q = ΔT/Rt, Ct dT/dt = qin − qoutτ = R CFluid and thermal systems extend the lumped-element idea to two more domains. In a liquid-level tank, the outflow is driven by the head and opposed by a fluid resistance, q = Δh/Rf, while the tank's cross-section acts as a fluid capacitance that stores volume. A mass balance, rate in minus rate out equals rate of storage, gives a first-order differential equation in the level. Heat flow behaves identically: a temperature difference drives the flow through a thermal resistance, q = ΔT/Rt, and a body's heat capacity Ct = m cp stores energy, so an energy balance gives a first-order equation in temperature. In both, the time constant τ = R C sets the pace: after one τ the system has moved 63 percent toward its final value.
The skills, taught in order
Five skills build fluid resistance and capacitance, the tank model, and their thermal equivalents.
4.1 Fluid resistance and capacitance
Fluid resistance relates flow to a pressure or head difference: q = Δh/Rf for laminar, linear flow. Fluid capacitance stores volume per unit head; for an open tank of cross-section A it is simply Cf = A. These two elements model most liquid systems.
4.2 The liquid-level tank
A tank fed at qin and draining through a resistance obeys the mass balance Cf dh/dt = qin − h/Rf. This is first order in the head h, with time constant τ = Rf Cf, the canonical fluid system.
| Domain | Effort (drives flow) | Resistance | Capacitance |
|---|---|---|---|
| Fluid | head or pressure | q = Δh/Rf | Cf = A (tank) |
| Thermal | temperature | q = ΔT/Rt | Ct = m cp |
Fluid and thermal elements in parallel. Both give first-order systems with τ = R C.
4.3 Thermal resistance and capacitance
Thermal resistance opposes heat flow: q = ΔT/Rt, with Rt = 1/(hA) for convection. Thermal capacitance stores energy as temperature: Ct = m cp. A body losing heat to its surroundings is the thermal twin of the draining tank.
4.4 The first-order time constant
In every case the balance yields a first-order equation whose time constant is τ = R C. A larger resistance or capacitance slows the system. The number of time constants, not the absolute time, tells how far the response has progressed: 63 percent at one τ, 95 percent at three, essentially complete at five.
4.5 Limits and nonlinearity
These linear models hold for laminar flow and modest temperature ranges. Turbulent outflow makes resistance depend on flow (q ∝ √h), and radiation makes thermal loss nonlinear. Recognising when to linearise, and around what operating point, keeps the first-order model valid.
Engineering connection: the thermal time constant of a temperature sensor sets how fast it can track a changing process, a direct link to measurement.
Worked example 1: time constant of a draining tank
An open tank has a cross-sectional area of 2 m² and a linear outlet resistance of 50 s/m² (so that q = h/Rf with h in metres and q in m³/s). Find the time constant of the level.
- ProblemFind the time constant of the tank level in Figure 1.
- Given / findA = 2 m², Rf = 50 s/m². Find τ.
- AssumptionsLaminar, linear outflow (q = h/Rf); constant cross-section, so Cf = A.
- ModelThe balance Cf dh/dt = qin − h/Rf is first order with τ = Rf Cf.
- EquationsCf = Aτ = Rf Cf
- SolveCf = 2 m². τ = 50 × 2 = 100 s.
- CheckUnits: (s/m²)(m²) = s. A bigger tank or a more restrictive outlet would lengthen τ, matching intuition about how slowly a large tank drains.
- ConclusionThe level settles with a 100 s time constant, reaching 63 percent of any change in 100 s and effectively settling in about 500 s.
Worked example 2: thermal time constant of a cooling part
A small part of mass 0.5 kg and specific heat 400 J/(kg·K) loses heat by convection with a coefficient h = 20 W/(m²·K) over an area of 0.1 m². Find its thermal capacitance, thermal resistance, and time constant.
- ProblemFind Ct, Rt, and τ for the cooling part in Figure 2.
- Given / findm = 0.5 kg, cp = 400 J/(kg·K), h = 20 W/(m²·K), A = 0.1 m². Find Ct, Rt, τ.
- AssumptionsLumped (uniform internal temperature); convection only; constant properties.
- ModelCt = m cp; Rt = 1/(hA); τ = Rt Ct.
- EquationsCt = m cpRt = 1/(hA)τ = Rt Ct
- SolveCt = 0.5 × 400 = 200 J/K. Rt = 1/(20 × 0.1) = 1/2 = 0.5 K/W. τ = 0.5 × 200 = 100 s.
- CheckUnits: (K/W)(J/K) = J/W = s. The lumped model is valid when the Biot number is small, which a thin, conductive part satisfies.
- ConclusionThe part cools with a 100 s time constant, the thermal twin of the tank, confirming that one first-order law spans both domains.
Misconceptions and diagnostics
| Mistake | Symptom | Diagnostic question | Correction |
|---|---|---|---|
| Linear law for turbulent flow | Tank model wrong at high flow | "Is the outflow laminar?" | Turbulent outflow gives q ∝ √h; linearise about an operating point. |
| Swapping R and C | Time constant has wrong units | "Which element stores, which opposes?" | Capacitance stores; resistance opposes flow. |
| Ignoring the Biot number | Lumped model used on a thick part | "Is the internal temperature uniform?" | The lumped model needs a small Biot number. |
| Forgetting specific heat | Ct taken as mass alone | "Did I multiply by cp?" | Ct = m cp, not m. |
Practice ladder
A tank of area 3 m² has an outlet resistance of 20 s/m². Find the time constant.
Show answer
τ = Rf Cf = 20 × 3 = 60 s.
A 1 kg aluminium block (cp = 900 J/(kg·K)) cools with Rt = 0.2 K/W. Find τ.
Show answer
Ct = 1 × 900 = 900 J/K. τ = Rt Ct = 0.2 × 900 = 180 s.
A thermometer has τ = 8 s. How long to read 95 percent of a sudden temperature change, and 99 percent?
Show answer
95 percent at about 3τ = 24 s; 99 percent at about 5τ = 40 s. The exponential needs several time constants to settle.
A temperature sensor must follow a process that changes every 30 s. Use the thermal time constant to argue what sensor property you need, and the trade-off.
What good work looks like
The sensor's τ must be well below 30 s (say a few seconds) to track the process, which means low thermal mass; the trade-off is that a small, fast sensor is more fragile and noisier.
Working with AI, and proving it yourself
Use AI as an examiner, not a solver
Portfolio task
Model a real fluid or thermal system, compute its time constant, and compare to a measured settling time, noting any nonlinearity you had to linearise.
Retrieval and spaced review
Closed notes. Answer out loud, then reveal.
1. Write the fluid resistance law.
q = Δh/Rf for linear, laminar flow.
2. What is the fluid capacitance of an open tank?
Its cross-sectional area, Cf = A.
3. Give the thermal resistance and capacitance.
Rt = 1/(hA) and Ct = m cp.
4. What is the time constant in both domains?
τ = R C, the product of resistance and capacitance.
5. When does the linear model break down?
For turbulent flow (q ∝ √h) or large temperature ranges; linearise about an operating point.
Textbook mapping
This module follows Karnopp, Margolis, and Rosenberg, System Dynamics: Modeling, Simulation, and Control of Mechatronic Systems, 5th edition. Use these references to read further.
| Topic in this module | Where to read more |
|---|---|
| Fluid resistance, capacitance, and tanks | Karnopp, Margolis & Rosenberg, Chapter 4 |
| Thermal resistance and capacitance | Karnopp, Margolis & Rosenberg, Chapter 4 |
| First-order time constants across domains | Karnopp, Margolis & Rosenberg, Chapter 6 |
Chapter numbers refer to the 5th edition. The fluid and thermal relations are standard, so any recent edition will align closely.