Integrated Thermal Process Design and Heat Recovery

The engineering problem

An integrated thermal design is a network decision. Improving one exchanger can increase or decrease the duty required elsewhere, and local surface efficiency does not guarantee system-level savings. A feasible proposal needs thermal feasibility, controllability, safe operating margins, constructability and measured baseline energy use.

Scope and learning outcomes

Combine process heating, thermal pinch concepts, insulation, uncertainty, operating cost and monitoring into defensible engineering decisions. Students should be able to select a governing equation, construct a measurable boundary condition, calculate a plausible output with its units, and explain a failure case before applying a correlation.

Concept map

8.1 Process heat recovery targets

Available hot- and cold-stream temperatures limit sensible heat recovery. A minimum approach temperature protects driving force but may shift the economically best area. A beverage plant pairs an outgoing hot batch with incoming cold feed. Check capacity-rate limits, batch availability and countercurrent temperature approaches. Representative relationship: Qmax_screen=Cmin*max(0,Th,in−Tc,in−DTmin).

8.2 Thermal energy audit and insulation

Insulating one surface can reduce heat loss while raising equipment-wall temperatures elsewhere; economic payback depends on measured hours and energy price. A pipe-insulation retrofit is proposed for a twenty-hour-per-day operation. Calculate annual thermal savings with burner efficiency and the limits of steady operating hours. Representative relationship: annual_cost=Qlosshoursprice/1000.

8.3 Measurement uncertainty and validation

A small hot/cold duty mismatch may reflect sensors rather than energy creation; propagated uncertainty must be considered before diagnosing leakage. A heat exchanger records inconsistent heat duties on the two sides. Rank which sensors affect duty uncertainty and design a repeatability test. Representative relationship: u(Q)^2≈sum[(dQ/dxi*u(xi))^2].

8.4 Digital process models and control

Thermal time constants, transport delay and variable properties affect heating-loop control; a fitted model requires validation against unseen operating changes. A temperature controller overshoots after a throughput change. Propose identification data, validate controller limits and distinguish model mismatch from actuator saturation. Representative relationship: tau=mcp/(hA) (ideal lumped element).

A disciplined workflow for this unit

Start with a diagram of a process boundary, including fluid inlets, outlets, phase boundaries, hot and cold external surfaces and any energy storage. Record which values are measurements, which are assumed, and which are calculated. Write the conservation statement before inserting a transfer correlation. After obtaining a numerical result, test the limiting case: the heat flow should vanish with zero temperature difference in a passive linear resistance model, while doubling the wall area should increase conduction duty if all other conditions remain fixed. For nonlinear radiation a twofold temperature increase cannot be treated as a twofold radiation increase.

Report results with their uncertainty drivers. For a heated process tube, a thermocouple may measure external wall temperature but not internal wall temperature, and a flow meter may give volumetric rather than mass rate. An apparent discrepancy between measured and predicted duty may be the result of imperfect knowledge of heat capacities, uncontrolled bypass or fouling. If a correlation is used beyond the printed range, label the result an extrapolation instead of assigning spurious precision.

Unit investigation method

Use a blank process sheet with five sections: known properties, governing physics, calculation, validation, and operating decision. At minimum compare two physically plausible operating states and identify what extra measurement would most efficiently distinguish them. Do not rely on the number of digits printed by a simulation to argue that a real process is controlled to that precision.

Unit 8, Lesson 1: Process heat recovery targets

Unit: Integrated thermal process design. Core formula: Qmax_screen=Cmin*max(0,Th,in−Tc,in−DTmin). Suggested model: [Laboratory 18 — Heat recovery minimum approach](../labs/lab-18.html).

Why this problem appears in process plants

A beverage plant pairs an outgoing hot batch with incoming cold feed. Available hot- and cold-stream temperatures limit sensible heat recovery. A minimum approach temperature protects driving force but may shift the economically best area. This topic is a decision about where energy travels, which mechanism controls its rate, and what evidence is necessary to justify a calculation. The same reported heater duty can correspond to distinct wall temperatures depending on geometry, flow and material properties. A credible engineer must therefore distinguish the energy rate demanded by a process from the thermal resistance that allows the equipment to deliver that rate.

A process operator experiences this distinction as an operating constraint rather than as an equation. When product throughput increases, an unchanged heating utility can be insufficient even when the external wall stays hot. A measurement campaign should record mass flow, inlet and outlet conditions, utility conditions, local wall temperatures when available, and the operating history that may affect fouling. Without these signals, increasing utility temperature may mask the actual cause of poor performance and worsen degradation.

Physical mechanism and assumptions

Available hot- and cold-stream temperatures limit sensible heat recovery. A minimum approach temperature protects driving force but may shift the economically best area. The equation selected here, Qmax_screen=Cmin*max(0,Th,in−Tc,in−DTmin), summarizes a limited physical model. Every term must refer to a declared boundary, an appropriate unit system, and a physically relevant value of temperature or material property. Before using the equation, identify whether the system is steady or transient, whether the fluid undergoes phase change, and whether the heat-transfer coefficient is known from a defensible correlation or simply a fitted number.

The engineer's first drawing should show an inward or outward normal for the relevant surface and a clear sign convention. Heat transport is directed down a temperature driving force in a passive system, although a controlled refrigeration cycle can move heat from a colder source to a hotter sink by supplying work. A temperature sensor located in a mixing chamber is not necessarily the right measurement of the wall temperature immediately adjacent to a process fluid. A fixed temperature boundary also differs from an imposed heat-flux boundary and gives a different mathematical problem.

Derivation and dimensional reasoning

Start from the first-law energy balance, accumulation = net incoming enthalpy + supplied heat + work, with the appropriate sign conventions. Introduce the constitutive expression Qmax_screen=Cmin*max(0,Th,in−Tc,in−DTmin) only after identifying the relevant energy-transfer mechanism. If a rate equation gives watts, confirm that conductivity is in W/(m K), surface area in m² and temperature difference in kelvin for a planar slab. The analogous exercise for a convection coefficient uses W/(m² K). For radiation, temperatures must be absolute because emission depends on the fourth power of thermodynamic temperature.

Terms that look superficially similar can require fundamentally different models. A cylindrical pipe wall does not have a constant radial area; a nonlinear radiation exchange cannot be replaced by a temperature-independent resistance without specifying a linearization point. Correlations expressed with Reynolds and Prandtl numbers need fluid properties at the proper reference temperature and geometry. In a thermally developing flow, the local Nusselt number may differ from its fully developed average. Write these distinctions directly beside the calculation rather than after the final answer.

Worked investigation: heat recovery minimum approach

We use a numerical teaching model with the following documented inputs: minimum approach = 12; Ch = 4600; Cc = 3200; Th = 145; Tc = 35. The primary control minimum approach is measured in K. Other parameters are Ch [W/K]; Cc [W/K]; Th [°C]; Tc [°C]. Using the model equation, the predicted optimistic recovery ceiling is 313.6 kW. This is an arithmetic result for assumed properties, not a measured industrial duty or a certified equipment rating.

Increase the primary control from 12 to 18.08 K while holding the other inputs fixed. The same model predicts 294.144 kW. Compare the absolute difference (-19.456 kW) with the percentage difference, and interpret the sign using the physical mechanism. Does this change chiefly alter the driving force, the thermal resistance, the material storage, or the rate at which the process moves energy? Explain why interpreting the curve without checking the other controlled inputs could be misleading.

The control range is deliberately finite. At very large or small values a simplified correlation can violate its validity conditions even if the numerical program continues to return finite numbers. For this lab, the important limitation is: Screening ceiling only; not a network pinch analysis or exchanger rating. Decide whether the default state is inside that limitation and state what kind of additional property data, geometry measurement, or experimental coefficient would be needed to increase confidence in the estimate.

Measurement strategy and uncertainty

Check capacity-rate limits, batch availability and countercurrent temperature approaches. Select one variable from the worked model that could realistically be measured on a plant floor, one that would normally be inferred, and one that is liable to drift over months of operation. Suggest how each could be independently checked. Thermal couples need installation and dynamic-response consideration, while flow meters require correct fluid properties and calibration. A sensor can be accurate at its measurement point yet irrelevant to the heat-transfer boundary being modeled.

Develop a first-order sensitivity calculation. If the output Y depends on parameter p, estimate dY/dp from two nearby runs of the simulator. If the measured parameter has uncertainty u(p), its approximate output contribution is abs(dY/dp)*u(p). This is a local linear approximation; it does not characterize poorly known geometry, phase-change instability or systematic installation bias. Keep uncertainty in the output's physical units until presenting any final percentage. Where a temperature difference is small, uncertainty can become the dominant proportion of the inferred transfer rate.

Process decision and limitation

Check capacity-rate limits, batch availability and countercurrent temperature approaches. An excellent model is still a screening tool if it ignores details that control real equipment performance. A preliminary calculation may support an experiment or identify the direction of a sensitivity, but it cannot substitute for a vendor curve, verified thermophysical properties, design margins, pressure-equipment requirements or an installed sensor audit. Present the conclusion as an evidence-backed statement with a named range of validity. Avoid claiming that a temperature or duty is safe solely because one steady-state expression has a positive numerical solution.

Independent worked case

[Read a fully contextualized numerical case: Heat recovery temperature ceiling](../worked-cases/unit-08-case-01.html).

Exercises

1. Starting with Qmax_screen=Cmin*max(0,Th,in−Tc,in−DTmin), identify the independent variables and produce a dimensionally consistent derivation or explanation of each term. Give an example where using the wrong unit would change the answer by at least a factor of 100.
2. Repeat the worked experiment with the control variable reduced by 10%. Compute the new output, the percent change, and compare the result to a first-order sensitivity estimate. Explain whether a linear approximation is justified.
3. Check capacity-rate limits, batch availability and countercurrent temperature approaches. Describe the minimum sensor set and a method for recording uncertainty in the field.
4. Construct an alternative physical model for the same plant scenario by changing one boundary condition. State which observation could distinguish the two models.
5. Design a sensitivity chart comparing two material-property assumptions. What input dominates the uncertainty at the baseline and why?

Reproducible laboratory

Run [Laboratory 18: Heat recovery minimum approach](../labs/lab-18.html), export 151 modeled points to CSV, and annotate where the correlation is appropriate. Cite the chosen property source and write an explicit statement about the limits of extrapolation. The laboratory provides a schematic response curve and numerical readout but cannot validate a design without measured plant data.

Unit 8, Lesson 2: Thermal energy audit and insulation

Unit: Integrated thermal process design. Core formula: annual_cost=Qlosshoursprice/1000. Suggested model: [Laboratory 17 — Heat-exchanger fouling penalty](../labs/lab-17.html).

Why this problem appears in process plants

A pipe-insulation retrofit is proposed for a twenty-hour-per-day operation. Insulating one surface can reduce heat loss while raising equipment-wall temperatures elsewhere; economic payback depends on measured hours and energy price. This topic is a decision about where energy travels, which mechanism controls its rate, and what evidence is necessary to justify a calculation. The same reported heater duty can correspond to distinct wall temperatures depending on geometry, flow and material properties. A credible engineer must therefore distinguish the energy rate demanded by a process from the thermal resistance that allows the equipment to deliver that rate.

A process operator experiences this distinction as an operating constraint rather than as an equation. When product throughput increases, an unchanged heating utility can be insufficient even when the external wall stays hot. A measurement campaign should record mass flow, inlet and outlet conditions, utility conditions, local wall temperatures when available, and the operating history that may affect fouling. Without these signals, increasing utility temperature may mask the actual cause of poor performance and worsen degradation.

Physical mechanism and assumptions

Insulating one surface can reduce heat loss while raising equipment-wall temperatures elsewhere; economic payback depends on measured hours and energy price. The equation selected here, annual_cost=Qlosshoursprice/1000, summarizes a limited physical model. Every term must refer to a declared boundary, an appropriate unit system, and a physically relevant value of temperature or material property. Before using the equation, identify whether the system is steady or transient, whether the fluid undergoes phase change, and whether the heat-transfer coefficient is known from a defensible correlation or simply a fitted number.

The engineer's first drawing should show an inward or outward normal for the relevant surface and a clear sign convention. Heat transport is directed down a temperature driving force in a passive system, although a controlled refrigeration cycle can move heat from a colder source to a hotter sink by supplying work. A temperature sensor located in a mixing chamber is not necessarily the right measurement of the wall temperature immediately adjacent to a process fluid. A fixed temperature boundary also differs from an imposed heat-flux boundary and gives a different mathematical problem.

Derivation and dimensional reasoning

Start from the first-law energy balance, accumulation = net incoming enthalpy + supplied heat + work, with the appropriate sign conventions. Introduce the constitutive expression annual_cost=Qlosshoursprice/1000 only after identifying the relevant energy-transfer mechanism. If a rate equation gives watts, confirm that conductivity is in W/(m K), surface area in m² and temperature difference in kelvin for a planar slab. The analogous exercise for a convection coefficient uses W/(m² K). For radiation, temperatures must be absolute because emission depends on the fourth power of thermodynamic temperature.

Terms that look superficially similar can require fundamentally different models. A cylindrical pipe wall does not have a constant radial area; a nonlinear radiation exchange cannot be replaced by a temperature-independent resistance without specifying a linearization point. Correlations expressed with Reynolds and Prandtl numbers need fluid properties at the proper reference temperature and geometry. In a thermally developing flow, the local Nusselt number may differ from its fully developed average. Write these distinctions directly beside the calculation rather than after the final answer.

Worked investigation: heat-exchanger fouling penalty

We use a numerical teaching model with the following documented inputs: deposit resistance = 0.001; hh = 850; hc = 1100; Rw = 0.0002; A = 28; LMTD = 30. The primary control deposit resistance is measured in m² K/W. Other parameters are hh [W/(m² K)]; hc [W/(m² K)]; Rw [m² K/W]; A [m²]; LMTD [K]. Using the model equation, the predicted available duty is 255.664 kW. This is an arithmetic result for assumed properties, not a measured industrial duty or a certified equipment rating.

Increase the primary control from 0.001 to 0.00148 m² K/W while holding the other inputs fixed. The same model predicts 223.074 kW. Compare the absolute difference (-32.5898 kW) with the percentage difference, and interpret the sign using the physical mechanism. Does this change chiefly alter the driving force, the thermal resistance, the material storage, or the rate at which the process moves energy? Explain why interpreting the curve without checking the other controlled inputs could be misleading.

The control range is deliberately finite. At very large or small values a simplified correlation can violate its validity conditions even if the numerical program continues to return finite numbers. For this lab, the important limitation is: Same-area resistance approximation, fixed LMTD and unchanged film coefficients. Decide whether the default state is inside that limitation and state what kind of additional property data, geometry measurement, or experimental coefficient would be needed to increase confidence in the estimate.

Measurement strategy and uncertainty

Calculate annual thermal savings with burner efficiency and the limits of steady operating hours. Select one variable from the worked model that could realistically be measured on a plant floor, one that would normally be inferred, and one that is liable to drift over months of operation. Suggest how each could be independently checked. Thermal couples need installation and dynamic-response consideration, while flow meters require correct fluid properties and calibration. A sensor can be accurate at its measurement point yet irrelevant to the heat-transfer boundary being modeled.

Develop a first-order sensitivity calculation. If the output Y depends on parameter p, estimate dY/dp from two nearby runs of the simulator. If the measured parameter has uncertainty u(p), its approximate output contribution is abs(dY/dp)*u(p). This is a local linear approximation; it does not characterize poorly known geometry, phase-change instability or systematic installation bias. Keep uncertainty in the output's physical units until presenting any final percentage. Where a temperature difference is small, uncertainty can become the dominant proportion of the inferred transfer rate.

Process decision and limitation

Calculate annual thermal savings with burner efficiency and the limits of steady operating hours. An excellent model is still a screening tool if it ignores details that control real equipment performance. A preliminary calculation may support an experiment or identify the direction of a sensitivity, but it cannot substitute for a vendor curve, verified thermophysical properties, design margins, pressure-equipment requirements or an installed sensor audit. Present the conclusion as an evidence-backed statement with a named range of validity. Avoid claiming that a temperature or duty is safe solely because one steady-state expression has a positive numerical solution.

Independent worked case

[Read a fully contextualized numerical case: Insulation and annual utility cost](../worked-cases/unit-08-case-02.html).

Exercises

1. Starting with annual_cost=Qlosshoursprice/1000, identify the independent variables and produce a dimensionally consistent derivation or explanation of each term. Give an example where using the wrong unit would change the answer by at least a factor of 100.
2. Repeat the worked experiment with the control variable reduced by 10%. Compute the new output, the percent change, and compare the result to a first-order sensitivity estimate. Explain whether a linear approximation is justified.
3. Calculate annual thermal savings with burner efficiency and the limits of steady operating hours. Describe the minimum sensor set and a method for recording uncertainty in the field.
4. Construct an alternative physical model for the same plant scenario by changing one boundary condition. State which observation could distinguish the two models.
5. Design a sensitivity chart comparing two material-property assumptions. What input dominates the uncertainty at the baseline and why?

Reproducible laboratory

Run [Laboratory 17: Heat-exchanger fouling penalty](../labs/lab-17.html), export 151 modeled points to CSV, and annotate where the correlation is appropriate. Cite the chosen property source and write an explicit statement about the limits of extrapolation. The laboratory provides a schematic response curve and numerical readout but cannot validate a design without measured plant data.

Unit 8, Lesson 3: Measurement uncertainty and validation

Unit: Integrated thermal process design. Core formula: u(Q)^2≈sum[(dQ/dxi*u(xi))^2]. Suggested model: [Laboratory 15 — Counterflow exchanger effectiveness](../labs/lab-15.html).

Why this problem appears in process plants

A heat exchanger records inconsistent heat duties on the two sides. A small hot/cold duty mismatch may reflect sensors rather than energy creation; propagated uncertainty must be considered before diagnosing leakage. This topic is a decision about where energy travels, which mechanism controls its rate, and what evidence is necessary to justify a calculation. The same reported heater duty can correspond to distinct wall temperatures depending on geometry, flow and material properties. A credible engineer must therefore distinguish the energy rate demanded by a process from the thermal resistance that allows the equipment to deliver that rate.

A process operator experiences this distinction as an operating constraint rather than as an equation. When product throughput increases, an unchanged heating utility can be insufficient even when the external wall stays hot. A measurement campaign should record mass flow, inlet and outlet conditions, utility conditions, local wall temperatures when available, and the operating history that may affect fouling. Without these signals, increasing utility temperature may mask the actual cause of poor performance and worsen degradation.

Physical mechanism and assumptions

A small hot/cold duty mismatch may reflect sensors rather than energy creation; propagated uncertainty must be considered before diagnosing leakage. The equation selected here, u(Q)^2≈sum[(dQ/dxi*u(xi))^2], summarizes a limited physical model. Every term must refer to a declared boundary, an appropriate unit system, and a physically relevant value of temperature or material property. Before using the equation, identify whether the system is steady or transient, whether the fluid undergoes phase change, and whether the heat-transfer coefficient is known from a defensible correlation or simply a fitted number.

The engineer's first drawing should show an inward or outward normal for the relevant surface and a clear sign convention. Heat transport is directed down a temperature driving force in a passive system, although a controlled refrigeration cycle can move heat from a colder source to a hotter sink by supplying work. A temperature sensor located in a mixing chamber is not necessarily the right measurement of the wall temperature immediately adjacent to a process fluid. A fixed temperature boundary also differs from an imposed heat-flux boundary and gives a different mathematical problem.

Derivation and dimensional reasoning

Start from the first-law energy balance, accumulation = net incoming enthalpy + supplied heat + work, with the appropriate sign conventions. Introduce the constitutive expression u(Q)^2≈sum[(dQ/dxi*u(xi))^2] only after identifying the relevant energy-transfer mechanism. If a rate equation gives watts, confirm that conductivity is in W/(m K), surface area in m² and temperature difference in kelvin for a planar slab. The analogous exercise for a convection coefficient uses W/(m² K). For radiation, temperatures must be absolute because emission depends on the fourth power of thermodynamic temperature.

Terms that look superficially similar can require fundamentally different models. A cylindrical pipe wall does not have a constant radial area; a nonlinear radiation exchange cannot be replaced by a temperature-independent resistance without specifying a linearization point. Correlations expressed with Reynolds and Prandtl numbers need fluid properties at the proper reference temperature and geometry. In a thermally developing flow, the local Nusselt number may differ from its fully developed average. Write these distinctions directly beside the calculation rather than after the final answer.

Worked investigation: counterflow exchanger effectiveness

We use a numerical teaching model with the following documented inputs: UA = 1800; Ch = 2800; Cc = 1900; Thin = 160; Tcin = 25. The primary control UA is measured in W/K. Other parameters are Ch [W/K]; Cc [W/K]; Thin [°C]; Tcin [°C]. Using the model equation, the predicted heat-transfer duty is 134.788 kW. This is an arithmetic result for assumed properties, not a measured industrial duty or a certified equipment rating.

Increase the primary control from 1800 to 4712 W/K while holding the other inputs fixed. The same model predicts 202.985 kW. Compare the absolute difference (68.1967 kW) with the percentage difference, and interpret the sign using the physical mechanism. Does this change chiefly alter the driving force, the thermal resistance, the material storage, or the rate at which the process moves energy? Explain why interpreting the curve without checking the other controlled inputs could be misleading.

The control range is deliberately finite. At very large or small values a simplified correlation can violate its validity conditions even if the numerical program continues to return finite numbers. For this lab, the important limitation is: Single-phase countercurrent; constant capacity rates, no external loss; not multiple pass. Decide whether the default state is inside that limitation and state what kind of additional property data, geometry measurement, or experimental coefficient would be needed to increase confidence in the estimate.

Measurement strategy and uncertainty

Rank which sensors affect duty uncertainty and design a repeatability test. Select one variable from the worked model that could realistically be measured on a plant floor, one that would normally be inferred, and one that is liable to drift over months of operation. Suggest how each could be independently checked. Thermal couples need installation and dynamic-response consideration, while flow meters require correct fluid properties and calibration. A sensor can be accurate at its measurement point yet irrelevant to the heat-transfer boundary being modeled.

Develop a first-order sensitivity calculation. If the output Y depends on parameter p, estimate dY/dp from two nearby runs of the simulator. If the measured parameter has uncertainty u(p), its approximate output contribution is abs(dY/dp)*u(p). This is a local linear approximation; it does not characterize poorly known geometry, phase-change instability or systematic installation bias. Keep uncertainty in the output's physical units until presenting any final percentage. Where a temperature difference is small, uncertainty can become the dominant proportion of the inferred transfer rate.

Process decision and limitation

Rank which sensors affect duty uncertainty and design a repeatability test. An excellent model is still a screening tool if it ignores details that control real equipment performance. A preliminary calculation may support an experiment or identify the direction of a sensitivity, but it cannot substitute for a vendor curve, verified thermophysical properties, design margins, pressure-equipment requirements or an installed sensor audit. Present the conclusion as an evidence-backed statement with a named range of validity. Avoid claiming that a temperature or duty is safe solely because one steady-state expression has a positive numerical solution.

Independent worked case

[Read a fully contextualized numerical case: Heat-duty uncertainty](../worked-cases/unit-08-case-03.html).

Exercises

1. Starting with u(Q)^2≈sum[(dQ/dxi*u(xi))^2], identify the independent variables and produce a dimensionally consistent derivation or explanation of each term. Give an example where using the wrong unit would change the answer by at least a factor of 100.
2. Repeat the worked experiment with the control variable reduced by 10%. Compute the new output, the percent change, and compare the result to a first-order sensitivity estimate. Explain whether a linear approximation is justified.
3. Rank which sensors affect duty uncertainty and design a repeatability test. Describe the minimum sensor set and a method for recording uncertainty in the field.
4. Construct an alternative physical model for the same plant scenario by changing one boundary condition. State which observation could distinguish the two models.
5. Design a sensitivity chart comparing two material-property assumptions. What input dominates the uncertainty at the baseline and why?

Reproducible laboratory

Run [Laboratory 15: Counterflow exchanger effectiveness](../labs/lab-15.html), export 151 modeled points to CSV, and annotate where the correlation is appropriate. Cite the chosen property source and write an explicit statement about the limits of extrapolation. The laboratory provides a schematic response curve and numerical readout but cannot validate a design without measured plant data.

Unit 8, Lesson 4: Digital process models and control

Unit: Integrated thermal process design. Core formula: tau=mcp/(hA) (ideal lumped element). Suggested model: [Laboratory 05 — Lumped-object cooling](../labs/lab-05.html).

Why this problem appears in process plants

A temperature controller overshoots after a throughput change. Thermal time constants, transport delay and variable properties affect heating-loop control; a fitted model requires validation against unseen operating changes. This topic is a decision about where energy travels, which mechanism controls its rate, and what evidence is necessary to justify a calculation. The same reported heater duty can correspond to distinct wall temperatures depending on geometry, flow and material properties. A credible engineer must therefore distinguish the energy rate demanded by a process from the thermal resistance that allows the equipment to deliver that rate.

A process operator experiences this distinction as an operating constraint rather than as an equation. When product throughput increases, an unchanged heating utility can be insufficient even when the external wall stays hot. A measurement campaign should record mass flow, inlet and outlet conditions, utility conditions, local wall temperatures when available, and the operating history that may affect fouling. Without these signals, increasing utility temperature may mask the actual cause of poor performance and worsen degradation.

Physical mechanism and assumptions

Thermal time constants, transport delay and variable properties affect heating-loop control; a fitted model requires validation against unseen operating changes. The equation selected here, tau=mcp/(hA) (ideal lumped element), summarizes a limited physical model. Every term must refer to a declared boundary, an appropriate unit system, and a physically relevant value of temperature or material property. Before using the equation, identify whether the system is steady or transient, whether the fluid undergoes phase change, and whether the heat-transfer coefficient is known from a defensible correlation or simply a fitted number.

The engineer's first drawing should show an inward or outward normal for the relevant surface and a clear sign convention. Heat transport is directed down a temperature driving force in a passive system, although a controlled refrigeration cycle can move heat from a colder source to a hotter sink by supplying work. A temperature sensor located in a mixing chamber is not necessarily the right measurement of the wall temperature immediately adjacent to a process fluid. A fixed temperature boundary also differs from an imposed heat-flux boundary and gives a different mathematical problem.

Derivation and dimensional reasoning

Start from the first-law energy balance, accumulation = net incoming enthalpy + supplied heat + work, with the appropriate sign conventions. Introduce the constitutive expression tau=mcp/(hA) (ideal lumped element) only after identifying the relevant energy-transfer mechanism. If a rate equation gives watts, confirm that conductivity is in W/(m K), surface area in m² and temperature difference in kelvin for a planar slab. The analogous exercise for a convection coefficient uses W/(m² K). For radiation, temperatures must be absolute because emission depends on the fourth power of thermodynamic temperature.

Terms that look superficially similar can require fundamentally different models. A cylindrical pipe wall does not have a constant radial area; a nonlinear radiation exchange cannot be replaced by a temperature-independent resistance without specifying a linearization point. Correlations expressed with Reynolds and Prandtl numbers need fluid properties at the proper reference temperature and geometry. In a thermally developing flow, the local Nusselt number may differ from its fully developed average. Write these distinctions directly beside the calculation rather than after the final answer.

Worked investigation: lumped-object cooling

We use a numerical teaching model with the following documented inputs: time = 120; T0 = 120; Ta = 23; h = 35; A = 0.08; mass = 0.8; cp = 850. The primary control time is measured in s. Other parameters are T0 [°C]; Ta [°C]; h [W/(m² K)]; A [m²]; mass [kg]; cp [J/(kg K)]. Using the model equation, the predicted predicted object temperature is 82.1806 °C. This is an arithmetic result for assumed properties, not a measured industrial duty or a certified equipment rating.

Increase the primary control from 120 to 196.8 s while holding the other inputs fixed. The same model predicts 66.136 °C. Compare the absolute difference (-16.0446 °C) with the percentage difference, and interpret the sign using the physical mechanism. Does this change chiefly alter the driving force, the thermal resistance, the material storage, or the rate at which the process moves energy? Explain why interpreting the curve without checking the other controlled inputs could be misleading.

The control range is deliberately finite. At very large or small values a simplified correlation can violate its validity conditions even if the numerical program continues to return finite numbers. For this lab, the important limitation is: Bi≪0.1, no internal gradient, constant h and cp. Decide whether the default state is inside that limitation and state what kind of additional property data, geometry measurement, or experimental coefficient would be needed to increase confidence in the estimate.

Measurement strategy and uncertainty

Propose identification data, validate controller limits and distinguish model mismatch from actuator saturation. Select one variable from the worked model that could realistically be measured on a plant floor, one that would normally be inferred, and one that is liable to drift over months of operation. Suggest how each could be independently checked. Thermal couples need installation and dynamic-response consideration, while flow meters require correct fluid properties and calibration. A sensor can be accurate at its measurement point yet irrelevant to the heat-transfer boundary being modeled.

Develop a first-order sensitivity calculation. If the output Y depends on parameter p, estimate dY/dp from two nearby runs of the simulator. If the measured parameter has uncertainty u(p), its approximate output contribution is abs(dY/dp)*u(p). This is a local linear approximation; it does not characterize poorly known geometry, phase-change instability or systematic installation bias. Keep uncertainty in the output's physical units until presenting any final percentage. Where a temperature difference is small, uncertainty can become the dominant proportion of the inferred transfer rate.

Process decision and limitation

Propose identification data, validate controller limits and distinguish model mismatch from actuator saturation. An excellent model is still a screening tool if it ignores details that control real equipment performance. A preliminary calculation may support an experiment or identify the direction of a sensitivity, but it cannot substitute for a vendor curve, verified thermophysical properties, design margins, pressure-equipment requirements or an installed sensor audit. Present the conclusion as an evidence-backed statement with a named range of validity. Avoid claiming that a temperature or duty is safe solely because one steady-state expression has a positive numerical solution.

Independent worked case

[Read a fully contextualized numerical case: Thermal control time scale](../worked-cases/unit-08-case-04.html).

Exercises

1. Starting with tau=mcp/(hA) (ideal lumped element), identify the independent variables and produce a dimensionally consistent derivation or explanation of each term. Give an example where using the wrong unit would change the answer by at least a factor of 100.
2. Repeat the worked experiment with the control variable reduced by 10%. Compute the new output, the percent change, and compare the result to a first-order sensitivity estimate. Explain whether a linear approximation is justified.
3. Propose identification data, validate controller limits and distinguish model mismatch from actuator saturation. Describe the minimum sensor set and a method for recording uncertainty in the field.
4. Construct an alternative physical model for the same plant scenario by changing one boundary condition. State which observation could distinguish the two models.
5. Design a sensitivity chart comparing two material-property assumptions. What input dominates the uncertainty at the baseline and why?

Reproducible laboratory

Run [Laboratory 05: Lumped-object cooling](../labs/lab-05.html), export 151 modeled points to CSV, and annotate where the correlation is appropriate. Cite the chosen property source and write an explicit statement about the limits of extrapolation. The laboratory provides a schematic response curve and numerical readout but cannot validate a design without measured plant data.

Applied project 08: Integrated thermal process design

Project brief

An integrated thermal design is a network decision. Improving one exchanger can increase or decrease the duty required elsewhere, and local surface efficiency does not guarantee system-level savings. A feasible proposal needs thermal feasibility, controllability, safe operating margins, constructability and measured baseline energy use. Build a small but defensible engineering evaluation for a hypothetical process facility. The target is not an equipment purchase specification: it is a traceable analysis of what would need to be known before design or procurement. Combine process heating, thermal pinch concepts, insulation, uncertainty, operating cost and monitoring into defensible engineering decisions.

Process case and data plan

The scenario begins with incomplete information. A process technician has measured bulk inlet and outlet temperatures, approximate flow and a surface temperature, while operations has an estimate of annual hours. Choose one of the four lesson topics below and define a physically plausible dataset with at least eight rows. Use SI units. You may use representative approximate fluid properties if you label them as assumptions; do not describe fabricated samples as field measurements.

• Option 1: Process heat recovery targets. A beverage plant pairs an outgoing hot batch with incoming cold feed. Check capacity-rate limits, batch availability and countercurrent temperature approaches. Begin with Qmax_screen=Cmin*max(0,Th,in−Tc,in−DTmin) and then examine which neglected effect might dominate.
• Option 2: Thermal energy audit and insulation. A pipe-insulation retrofit is proposed for a twenty-hour-per-day operation. Calculate annual thermal savings with burner efficiency and the limits of steady operating hours. Begin with annual_cost=Qlosshoursprice/1000 and then examine which neglected effect might dominate.
• Option 3: Measurement uncertainty and validation. A heat exchanger records inconsistent heat duties on the two sides. Rank which sensors affect duty uncertainty and design a repeatability test. Begin with u(Q)^2≈sum[(dQ/dxi*u(xi))^2] and then examine which neglected effect might dominate.
• Option 4: Digital process models and control. A temperature controller overshoots after a throughput change. Propose identification data, validate controller limits and distinguish model mismatch from actuator saturation. Begin with tau=mcp/(hA) (ideal lumped element) and then examine which neglected effect might dominate.

Work packages

A — Model the thermal circuit. Produce a diagram of thermal boundaries and flow directions. Record each resistance, energy-storage term, or imposed heat flux. Separate first-law calculations from empirical closure equations. Verify one result using an independent algebraic or numerical method.

B — Establish a baseline. Choose reasonable conditions, calculate heat duty in watts or kilowatts and identify one required material property that is not supplied. Define how a student would obtain that property from an authoritative source and which temperature or pressure governs the lookup.

C — Experiment and plot. Run two relevant interactive laboratories and export the curves. Compare at least three operating points, keeping all but one input fixed at a time. Include a table with the model assumption, the changed parameter and the resulting percentage change. Explain why some apparent trends would fail in a full system because stream capacity rates, pressure drop or properties co-vary.

D — Measurement and risk. Write a list of five instruments or process observations necessary to validate the simplified model. Include likely systematic and random errors, a rough range of plausible input uncertainty, and a note about unavailable information such as surface condition or contact resistance. Do not infer any pressure-vessel or worker-safety adequacy from the simulator.

E — Decision memorandum. Provide a two-page recommendation that ranks measurement, retrofit and operating-control options by estimated benefit and uncertainty. For energy economics, use explicit unit price, operating hours and equipment efficiency assumptions. Highlight which conclusions would change under a realistic alternative case.

Acceptance criteria

Deliver a readable process schematic, a table of inputs and sources, equation derivations with units, one independent verification, at least two experimental response plots, a measurement strategy and a decision memo. A technically unsupported precise number is not acceptable. The model is a teaching artifact, not a certified plant design.

Worked case 8.1: Heat recovery temperature ceiling

Process situation and calculation

An outgoing stream at 145°C meets an incoming feed at 35°C. Their heat-capacity rates are 4.6 and 3.2 kW/K. With a 12 K minimum approach, a simple optimistic recovered-duty ceiling is 3.2×(145−35−12)=313.6 kW. Actual countercurrent exchanger performance requires UA and feasibility of both terminal approaches.

Physical interpretation and possible error

Do not describe the screening ceiling as achievable duty, and do not use it as a formal pinch-analysis result for multiple hot and cold process streams with changing heat capacity.

Verification tasks

1. Annotate a sketch with heat-flow directions and temperature measurement positions. Explain whether the quoted numerical value is a rate, a transport coefficient, a stored-energy scale or a transient timescale.
2. Repeat the calculation after changing one physical parameter by 10%. Identify whether the result responds linearly, logarithmically or with a power law, and derive that dependence rather than relying only on a graph.
3. State which property or correlation requires independent external verification. What uncertainty would a surface thermocouple, flow meter or geometry measurement introduce?
4. Describe the limiting condition under which the numerical relationship is not meaningful. Do not generalize beyond its stated geometry or thermal regime.

Scope: This is an illustrative, parameter-specific calculation for education, not an equipment rating or validation of actual process data.

Worked case 8.2: Insulation and annual utility cost

Process situation and calculation

An insulated process surface loses an estimated 3 kW less than before, operating 5500 h/year. The nominal thermal energy reduction is 16,500 kWh/year. At utility heat cost 0.08 currency/kWh, the corresponding value is 1320 currency/year before installation, maintenance and financing costs.

Physical interpretation and possible error

If the 0.08 value is the fuel price rather than delivered heat cost, burner efficiency must be included; otherwise the economics mix thermal and purchased energy units.

Verification tasks

1. Annotate a sketch with heat-flow directions and temperature measurement positions. Explain whether the quoted numerical value is a rate, a transport coefficient, a stored-energy scale or a transient timescale.
2. Repeat the calculation after changing one physical parameter by 10%. Identify whether the result responds linearly, logarithmically or with a power law, and derive that dependence rather than relying only on a graph.
3. State which property or correlation requires independent external verification. What uncertainty would a surface thermocouple, flow meter or geometry measurement introduce?
4. Describe the limiting condition under which the numerical relationship is not meaningful. Do not generalize beyond its stated geometry or thermal regime.

Scope: This is an illustrative, parameter-specific calculation for education, not an equipment rating or validation of actual process data.

Worked case 8.3: Heat-duty uncertainty

Process situation and calculation

For a single-phase stream, Q=mdot cp ΔT. If flowrate has 2% relative standard uncertainty, cp 3% and measured temperature rise 4%, and the errors are independent, the combined relative standard uncertainty is sqrt(2²+3²+4²)%≈5.39%.

Physical interpretation and possible error

This is linear uncertainty propagation, not a full treatment of correlated errors or bias. At very small ΔT, absolute thermocouple uncertainty can dominate and invalidate a constant percentage assumption.

Verification tasks

1. Annotate a sketch with heat-flow directions and temperature measurement positions. Explain whether the quoted numerical value is a rate, a transport coefficient, a stored-energy scale or a transient timescale.
2. Repeat the calculation after changing one physical parameter by 10%. Identify whether the result responds linearly, logarithmically or with a power law, and derive that dependence rather than relying only on a graph.
3. State which property or correlation requires independent external verification. What uncertainty would a surface thermocouple, flow meter or geometry measurement introduce?
4. Describe the limiting condition under which the numerical relationship is not meaningful. Do not generalize beyond its stated geometry or thermal regime.

Scope: This is an illustrative, parameter-specific calculation for education, not an equipment rating or validation of actual process data.

Worked case 8.4: Thermal control time scale

Process situation and calculation

A well-mixed liquid tank holds 150 kg of liquid with cp=4000 J/(kg K), cooled through hA=1500 W/K. Its ideal first-order time constant is mcp/hA=400 s. A controller tuned for a tiny thermocouple with a 5 s response may oscillate when applied to the much slower whole-tank loop.

Physical interpretation and possible error

Transport delay, nonlinear valves and changes in mass invalidate a pure lumped model. Record step-response data at more than one operating point before accepting a digital-twin model as predictive.

Verification tasks

1. Annotate a sketch with heat-flow directions and temperature measurement positions. Explain whether the quoted numerical value is a rate, a transport coefficient, a stored-energy scale or a transient timescale.
2. Repeat the calculation after changing one physical parameter by 10%. Identify whether the result responds linearly, logarithmically or with a power law, and derive that dependence rather than relying only on a graph.
3. State which property or correlation requires independent external verification. What uncertainty would a surface thermocouple, flow meter or geometry measurement introduce?
4. Describe the limiting condition under which the numerical relationship is not meaningful. Do not generalize beyond its stated geometry or thermal regime.

Scope: This is an illustrative, parameter-specific calculation for education, not an equipment rating or validation of actual process data.

Laboratory 15 — Counterflow exchanger effectiveness

Linked unit: 6. Measured output: Heat-transfer duty (kW). Equation: Q=eps Cmin(Thi−Tci); eps from countercurrent NTU.

Experiment objective and theory

This interactive experiment isolates one thermal relationship so that parameter sensitivity can be observed without the distraction of a full process model. It does not solve an entire plant; the objective is to learn which inputs the mathematical expression regards as controlling, which inputs the engineer would have to obtain from measurements, and which unmodeled effects could reverse an apparent benefit. Start by writing the direction of expected heat transfer. For an entirely passive system, check that the response is compatible with the stated hot and cold temperatures and the sign convention.

Validity boundary: Single-phase countercurrent; constant capacity rates, no external loss; not multiple pass. The simplified correlation treats the remaining conditions as fixed even if that assumption would not survive an actual change in flow, phase or material temperature. If the laboratory displays a finite number at an invalid combination, do not call that number a validated prediction. Temperature-dependent properties should be evaluated at a defensible reference state; this model deliberately does not call a property database.

Default test and sensitivity measurements

Primary variable: UA, unit W/K. Initial value: 1800 W/K. Predicted output: 134.7882 kW. Your result may round on the web page. Confirm the orders of magnitude independently with the printed equation.

At the first selected test value 4040 W/K, the reference model gives 193.1841 kW. At the second test value 16010 W/K, it gives 250.7456 kW. Do not draw a straight-line conclusion from these two numbers without checking the full response curve. Export the CSV, calculate a central finite difference near the default value, and compare it with the curvature of the graph. Keep the parameter-change interval small enough for a sensitivity calculation but large enough to avoid rounding effects.

Controlled-variable experiment

1. Record the baseline output and plot before touching any slider. Explain whether a temperature difference, resistance, heat capacity or dimensionless correlation establishes the result.
2. Move the primary slider to its low and high limits. Identify the graph region where an extrapolation of the physics would be most dangerous.
3. Restore the baseline. Change one secondary material or geometry parameter by approximately 10%. State the expected direction before observing it; then describe the difference between your prediction and the numerical response.
4. Use reset to reproduce the default experiment, and export the CSV. Check the exported data for units, sign and monotonicity. A curve with a sharp numerical discontinuity should be investigated rather than smoothed away automatically.
5. Prepare a four-column comparison of baseline, adjusted input, adjusted output and physical interpretation. Finish with a statement of the lab model's limitations and what one real measurement would validate it.

Interpretation and field extension

A computational curve is a relationship under controlled assumptions, not a set of experimental observations. In a plant, thermal conductivity, viscosity and heat capacity change with temperature; fouling deposits grow in time; and equipment area can be effectively reduced by flow maldistribution. Explain whether the chosen experiment can distinguish these effects. If not, propose a paired measurement, such as wall temperatures and stream flow, that creates additional evidence. Never claim improved process safety or energy economics solely from a single slider result.

Questions: (a) Which input appears with the highest exponent in the local law? (b) Could an unmeasured contact resistance mimic the same trend? (c) Does the graph represent a local coefficient, an area-average value or a whole-unit thermal duty? (d) What would an alternative geometry change? (e) What independent instrument or property reference would be required before industrial scale-up?

Laboratory 16 — LMTD driving force

Linked unit: 6. Measured output: Countercurrent LMTD (K). Equation: LMTD=(ΔT1−ΔT2)/ln(ΔT1/ΔT2).

Experiment objective and theory

This interactive experiment isolates one thermal relationship so that parameter sensitivity can be observed without the distraction of a full process model. It does not solve an entire plant; the objective is to learn which inputs the mathematical expression regards as controlling, which inputs the engineer would have to obtain from measurements, and which unmodeled effects could reverse an apparent benefit. Start by writing the direction of expected heat transfer. For an entirely passive system, check that the response is compatible with the stated hot and cold temperatures and the sign convention.

Validity boundary: Terminal differences must be positive; no multipass correction. The simplified correlation treats the remaining conditions as fixed even if that assumption would not survive an actual change in flow, phase or material temperature. If the laboratory displays a finite number at an invalid combination, do not call that number a validated prediction. Temperature-dependent properties should be evaluated at a defensible reference state; this model deliberately does not call a property database.

Default test and sensitivity measurements

Primary variable: hot inlet, unit °C. Initial value: 150 °C. Predicted output: 68.6069 K. Your result may round on the web page. Confirm the orders of magnitude independently with the printed equation.

At the first selected test value 102 °C, the reference model gives 43.28085 K. At the second test value 198 °C, it gives 88.9563 K. Do not draw a straight-line conclusion from these two numbers without checking the full response curve. Export the CSV, calculate a central finite difference near the default value, and compare it with the curvature of the graph. Keep the parameter-change interval small enough for a sensitivity calculation but large enough to avoid rounding effects.

Controlled-variable experiment

1. Record the baseline output and plot before touching any slider. Explain whether a temperature difference, resistance, heat capacity or dimensionless correlation establishes the result.
2. Move the primary slider to its low and high limits. Identify the graph region where an extrapolation of the physics would be most dangerous.
3. Restore the baseline. Change one secondary material or geometry parameter by approximately 10%. State the expected direction before observing it; then describe the difference between your prediction and the numerical response.
4. Use reset to reproduce the default experiment, and export the CSV. Check the exported data for units, sign and monotonicity. A curve with a sharp numerical discontinuity should be investigated rather than smoothed away automatically.
5. Prepare a four-column comparison of baseline, adjusted input, adjusted output and physical interpretation. Finish with a statement of the lab model's limitations and what one real measurement would validate it.

Interpretation and field extension

A computational curve is a relationship under controlled assumptions, not a set of experimental observations. In a plant, thermal conductivity, viscosity and heat capacity change with temperature; fouling deposits grow in time; and equipment area can be effectively reduced by flow maldistribution. Explain whether the chosen experiment can distinguish these effects. If not, propose a paired measurement, such as wall temperatures and stream flow, that creates additional evidence. Never claim improved process safety or energy economics solely from a single slider result.

Questions: (a) Which input appears with the highest exponent in the local law? (b) Could an unmeasured contact resistance mimic the same trend? (c) Does the graph represent a local coefficient, an area-average value or a whole-unit thermal duty? (d) What would an alternative geometry change? (e) What independent instrument or property reference would be required before industrial scale-up?

Unit 08 — Assessment: Integrated thermal process design

Instructions: Work at least three problems with complete units, a diagram, assumptions and a sensitivity check. These questions emphasize independent explanation rather than recognizing a repeated numerical template.

Topic 8.1: Process heat recovery targets

1. Derivation: Starting with `Qmax_screen=Cmin*max(0,Th,in−Tc,in−DTmin)`, derive the expression from stated physical assumptions or explain the empirical closure that prevents a first-principles derivation. Define the system boundary and identify which values are measurable.

2. Engineering investigation: A beverage plant pairs an outgoing hot batch with incoming cold feed. Produce two competing hypotheses that could explain a 12% deviation between the calculated heat duty and the measured duty; for each give an instrument or a controlled test.

3. Critical evaluation: Check capacity-rate limits, batch availability and countercurrent temperature approaches. State how the result would change for a tenfold increase in one relevant physical length, and whether the original correlation remains valid.

Topic 8.2: Thermal energy audit and insulation

1. Derivation: Starting with `annual_cost=Qloss*hours*price/1000`, derive the expression from stated physical assumptions or explain the empirical closure that prevents a first-principles derivation. Define the system boundary and identify which values are measurable.

2. Engineering investigation: A pipe-insulation retrofit is proposed for a twenty-hour-per-day operation. Produce two competing hypotheses that could explain a 12% deviation between the calculated heat duty and the measured duty; for each give an instrument or a controlled test.

3. Critical evaluation: Calculate annual thermal savings with burner efficiency and the limits of steady operating hours. State how the result would change for a tenfold increase in one relevant physical length, and whether the original correlation remains valid.

Topic 8.3: Measurement uncertainty and validation

1. Derivation: Starting with `u(Q)^2≈sum[(dQ/dxi*u(xi))^2]`, derive the expression from stated physical assumptions or explain the empirical closure that prevents a first-principles derivation. Define the system boundary and identify which values are measurable.

2. Engineering investigation: A heat exchanger records inconsistent heat duties on the two sides. Produce two competing hypotheses that could explain a 12% deviation between the calculated heat duty and the measured duty; for each give an instrument or a controlled test.

3. Critical evaluation: Rank which sensors affect duty uncertainty and design a repeatability test. State how the result would change for a tenfold increase in one relevant physical length, and whether the original correlation remains valid.

Topic 8.4: Digital process models and control

1. Derivation: Starting with `tau=m*cp/(h*A) (ideal lumped element)`, derive the expression from stated physical assumptions or explain the empirical closure that prevents a first-principles derivation. Define the system boundary and identify which values are measurable.

2. Engineering investigation: A temperature controller overshoots after a throughput change. Produce two competing hypotheses that could explain a 12% deviation between the calculated heat duty and the measured duty; for each give an instrument or a controlled test.

3. Critical evaluation: Propose identification data, validate controller limits and distinguish model mismatch from actuator saturation. State how the result would change for a tenfold increase in one relevant physical length, and whether the original correlation remains valid.

Assessment rubric

Conceptual accuracy (30%): state the correct physics, not just a numerical recipe. Traceable units and sign conventions (20%): make every conversion explicit. Sensitivity and validity (25%): quantify a change and identify non-applicable regimes. Engineering communication (25%): diagrams, model limits and measurements must support the recommendation.

The exercises are deliberately open; different property assumptions can produce different numerical outcomes. Clearly distinguish supplied values, literature property values and measurements. A plausible order of magnitude with transparent limitations is more useful than an unjustifiably precise number.