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.

01

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.
0 or 1 weak itemsContinue with this module.
2 weak itemsRevisit lumped modeling in Module 1.
3 or more weak itemsReview the RC time constant idea in Electrical Circuits, Module 4.
02

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 C

Fluid 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 skill works when: you write the balance, identify R and C for the domain, and form τ = R C.
The skill breaks down when: a turbulent resistance is treated as linear, or capacitance and resistance are swapped.
The concept. The tank stores volume as head and lets it out through a resistance. The balance is first order, and the time constant R C governs the draining, just as it governs a cooling part.
03

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.

DomainEffort (drives flow)ResistanceCapacitance
Fluidhead or pressureq = Δh/RfCf = A (tank)
Thermaltemperatureq = ΔT/RtCt = 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.

04

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.

Figure 1. The tank area is the fluid capacitance and the outlet is the resistance. Their product is the time constant over which the level settles toward its steady value.
  1. ProblemFind the time constant of the tank level in Figure 1.
  2. Given / findA = 2 m², Rf = 50 s/m². Find τ.
  3. AssumptionsLaminar, linear outflow (q = h/Rf); constant cross-section, so Cf = A.
  4. ModelThe balance Cf dh/dt = qin − h/Rf is first order with τ = Rf Cf.
  5. EquationsCf = Aτ = Rf Cf
  6. SolveCf = 2 m². τ = 50 × 2 = 100 s.
  7. CheckUnits: (s/m²)(m²) = s. A bigger tank or a more restrictive outlet would lengthen τ, matching intuition about how slowly a large tank drains.
  8. 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.
Result. τ = 100 s.
05

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.

Figure 2. The part cools exponentially toward ambient. Its heat capacity is the thermal capacitance and the convective film is the thermal resistance, and their product is the cooling time constant.
  1. ProblemFind Ct, Rt, and τ for the cooling part in Figure 2.
  2. Given / findm = 0.5 kg, cp = 400 J/(kg·K), h = 20 W/(m²·K), A = 0.1 m². Find Ct, Rt, τ.
  3. AssumptionsLumped (uniform internal temperature); convection only; constant properties.
  4. ModelCt = m cp; Rt = 1/(hA); τ = Rt Ct.
  5. EquationsCt = m cpRt = 1/(hA)τ = Rt Ct
  6. SolveCt = 0.5 × 400 = 200 J/K. Rt = 1/(20 × 0.1) = 1/2 = 0.5 K/W. τ = 0.5 × 200 = 100 s.
  7. 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.
  8. ConclusionThe part cools with a 100 s time constant, the thermal twin of the tank, confirming that one first-order law spans both domains.
Result. Ct = 200 J/K, Rt = 0.5 K/W, τ = 100 s.
06

Misconceptions and diagnostics

MistakeSymptomDiagnostic questionCorrection
Linear law for turbulent flowTank model wrong at high flow"Is the outflow laminar?"Turbulent outflow gives q ∝ √h; linearise about an operating point.
Swapping R and CTime constant has wrong units"Which element stores, which opposes?"Capacitance stores; resistance opposes flow.
Ignoring the Biot numberLumped model used on a thick part"Is the internal temperature uniform?"The lumped model needs a small Biot number.
Forgetting specific heatCt taken as mass alone"Did I multiply by cp?"Ct = m cp, not m.
07

Practice ladder

Level 1 · Direct skill

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.

Level 2 · Mixed concept

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.

Level 3 · Independent problem

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.

Transfer task | Real engineering

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.

08

Working with AI, and proving it yourself

Use AI as an examiner, not a solver

"Check that my time constant has units of seconds."
"Give me three systems; I will identify R and C for each domain."
"Solve this tank for me." Writing the balance and forming τ is the skill.
"Is the lumped model valid?" Checking the Biot number yourself is the point.

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.

Must include: a balance equation, an R and C, a computed τ, and a measured comparison.
09

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.

TodayFinish this quiz and Levels 1 and 2 of the ladder.
+1 dayRe-derive the tank balance and τ from a blank page.
+3 daysFind the time constant of three new systems.
+7 daysUnify the domains through energy methods, Module 5.
+30 daysReuse the thermal time constant when choosing a sensor.
10

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 moduleWhere to read more
Fluid resistance, capacitance, and tanksKarnopp, Margolis & Rosenberg, Chapter 4
Thermal resistance and capacitanceKarnopp, Margolis & Rosenberg, Chapter 4
First-order time constants across domainsKarnopp, 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.