Transient Conduction and Thermal Response

The engineering problem

A transient field requires both a constitutive law for flux and an energy-storage term. A metal bead may be nearly isothermal while a frozen food slab has a steep spatial gradient. Using a first-order exponential for both would hide internal temperature differences that determine product quality and thermal damage.

Scope and learning outcomes

Choose between lumped-capacitance, one-dimensional diffusion, and numerical thermal models with explicit time-scale and Biot-number checks. 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

3.1 Lumped-capacitance cooling

A low-Biot object follows an exponential approach to its surroundings when h, heat capacity and surroundings remain approximately constant. A small thermocouple bead is moved from hot air to cooling liquid. Compare the sensor response time with the actual bath temperature-change time. Representative relationship: T(t)=Tinf+(T0−Tinf)exp(−hAt/(m*cp)).

3.2 Thermal penetration and diffusion time

The thermal diffusion length grows as the square root of time, not linearly, because energy must cross progressively more material. A batch oven heats thick slabs of different materials. Predict the consequence of doubling thickness and compare a 1D thermal field with an isothermal estimate. Representative relationship: ell=sqrt(alpha*t); td~L^2/alpha.

3.3 Heisler limits and spatial gradients

When Biot is not small, the center and surface do not cool at the same rate; exact solutions depend on geometry and boundary conditions. A cylindrical food can must reach a minimum center temperature. Explain why average temperature cannot establish lethality or a minimum-center-temperature guarantee. Representative relationship: Fo=alpha*t/Lc^2.

3.4 Numerical conduction and convergence

Discretization must respect stability, conservation and boundary fluxes; reducing spatial step without adjusting time step can destabilize an explicit scheme. A thermally cycling process vessel wall is modeled with a finite-difference mesh. Run a mesh/time-step comparison and identify whether a spike is physical or numerical. Representative relationship: Fo_grid=alpha*dt/dx^2 ≤ 1/2 (1D explicit).

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 3, Lesson 1: Lumped-capacitance cooling

Unit: Transient conduction and thermal response. Core formula: T(t)=Tinf+(T0−Tinf)exp(−hAt/(m*cp)). Suggested model: [Laboratory 05 — Lumped-object cooling](../labs/lab-05.html).

Why this problem appears in process plants

A small thermocouple bead is moved from hot air to cooling liquid. A low-Biot object follows an exponential approach to its surroundings when h, heat capacity and surroundings remain approximately constant. 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 low-Biot object follows an exponential approach to its surroundings when h, heat capacity and surroundings remain approximately constant. The equation selected here, T(t)=Tinf+(T0−Tinf)exp(−hAt/(m*cp)), 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 T(t)=Tinf+(T0−Tinf)exp(−hAt/(m*cp)) 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

Compare the sensor response time with the actual bath temperature-change time. 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

Compare the sensor response time with the actual bath temperature-change time. 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: Cooling time constant](../worked-cases/unit-03-case-01.html).

Exercises

1. Starting with T(t)=Tinf+(T0−Tinf)exp(−hAt/(m*cp)), 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. Compare the sensor response time with the actual bath temperature-change time. 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.

Unit 3, Lesson 2: Thermal penetration and diffusion time

Unit: Transient conduction and thermal response. Core formula: ell=sqrt(alpha*t); td~L^2/alpha. Suggested model: [Laboratory 06 — Diffusion penetration length](../labs/lab-06.html).

Why this problem appears in process plants

A batch oven heats thick slabs of different materials. The thermal diffusion length grows as the square root of time, not linearly, because energy must cross progressively more material. 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

The thermal diffusion length grows as the square root of time, not linearly, because energy must cross progressively more material. The equation selected here, ell=sqrt(alpha*t); td~L^2/alpha, 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 ell=sqrt(alpha*t); td~L^2/alpha 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: diffusion penetration length

We use a numerical teaching model with the following documented inputs: time = 300; alpha = 1.2e-05. The primary control time is measured in s. Other parameters are alpha [m²/s]. Using the model equation, the predicted diffusion length is 60 mm. This is an arithmetic result for assumed properties, not a measured industrial duty or a certified equipment rating.

Increase the primary control from 300 to 444 s while holding the other inputs fixed. The same model predicts 72.9932 mm. Compare the absolute difference (12.9932 mm) 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: Order-of-magnitude penetration length, not an exact isotherm or thermal front. 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

Predict the consequence of doubling thickness and compare a 1D thermal field with an isothermal estimate. 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

Predict the consequence of doubling thickness and compare a 1D thermal field with an isothermal estimate. 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: Thickness-squared heating penalty](../worked-cases/unit-03-case-02.html).

Exercises

1. Starting with ell=sqrt(alpha*t); td~L^2/alpha, 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. Predict the consequence of doubling thickness and compare a 1D thermal field with an isothermal estimate. 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 06: Diffusion penetration length](../labs/lab-06.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 3, Lesson 3: Heisler limits and spatial gradients

Unit: Transient conduction and thermal response. Core formula: Fo=alpha*t/Lc^2. Suggested model: [Laboratory 05 — Lumped-object cooling](../labs/lab-05.html).

Why this problem appears in process plants

A cylindrical food can must reach a minimum center temperature. When Biot is not small, the center and surface do not cool at the same rate; exact solutions depend on geometry and boundary conditions. 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

When Biot is not small, the center and surface do not cool at the same rate; exact solutions depend on geometry and boundary conditions. The equation selected here, Fo=alpha*t/Lc^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 Fo=alpha*t/Lc^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: 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

Explain why average temperature cannot establish lethality or a minimum-center-temperature guarantee. 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

Explain why average temperature cannot establish lethality or a minimum-center-temperature guarantee. 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: Center temperature versus wall temperature](../worked-cases/unit-03-case-03.html).

Exercises

1. Starting with Fo=alpha*t/Lc^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. Explain why average temperature cannot establish lethality or a minimum-center-temperature guarantee. 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.

Unit 3, Lesson 4: Numerical conduction and convergence

Unit: Transient conduction and thermal response. Core formula: Fo_grid=alpha*dt/dx^2 ≤ 1/2 (1D explicit). Suggested model: [Laboratory 06 — Diffusion penetration length](../labs/lab-06.html).

Why this problem appears in process plants

A thermally cycling process vessel wall is modeled with a finite-difference mesh. Discretization must respect stability, conservation and boundary fluxes; reducing spatial step without adjusting time step can destabilize an explicit scheme. 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

Discretization must respect stability, conservation and boundary fluxes; reducing spatial step without adjusting time step can destabilize an explicit scheme. The equation selected here, Fo_grid=alpha*dt/dx^2 ≤ 1/2 (1D explicit), 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 Fo_grid=alpha*dt/dx^2 ≤ 1/2 (1D explicit) 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: diffusion penetration length

We use a numerical teaching model with the following documented inputs: time = 300; alpha = 1.2e-05. The primary control time is measured in s. Other parameters are alpha [m²/s]. Using the model equation, the predicted diffusion length is 60 mm. This is an arithmetic result for assumed properties, not a measured industrial duty or a certified equipment rating.

Increase the primary control from 300 to 444 s while holding the other inputs fixed. The same model predicts 72.9932 mm. Compare the absolute difference (12.9932 mm) 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: Order-of-magnitude penetration length, not an exact isotherm or thermal front. 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

Run a mesh/time-step comparison and identify whether a spike is physical or numerical. 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

Run a mesh/time-step comparison and identify whether a spike is physical or numerical. 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: Explicit finite-difference stability](../worked-cases/unit-03-case-04.html).

Exercises

1. Starting with Fo_grid=alpha*dt/dx^2 ≤ 1/2 (1D explicit), 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. Run a mesh/time-step comparison and identify whether a spike is physical or numerical. 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 06: Diffusion penetration length](../labs/lab-06.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 03: Transient conduction and thermal response

Project brief

A transient field requires both a constitutive law for flux and an energy-storage term. A metal bead may be nearly isothermal while a frozen food slab has a steep spatial gradient. Using a first-order exponential for both would hide internal temperature differences that determine product quality and thermal damage. 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. Choose between lumped-capacitance, one-dimensional diffusion, and numerical thermal models with explicit time-scale and Biot-number checks.

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: Lumped-capacitance cooling. A small thermocouple bead is moved from hot air to cooling liquid. Compare the sensor response time with the actual bath temperature-change time. Begin with T(t)=Tinf+(T0−Tinf)exp(−hAt/(m*cp)) and then examine which neglected effect might dominate.
• Option 2: Thermal penetration and diffusion time. A batch oven heats thick slabs of different materials. Predict the consequence of doubling thickness and compare a 1D thermal field with an isothermal estimate. Begin with ell=sqrt(alpha*t); td~L^2/alpha and then examine which neglected effect might dominate.
• Option 3: Heisler limits and spatial gradients. A cylindrical food can must reach a minimum center temperature. Explain why average temperature cannot establish lethality or a minimum-center-temperature guarantee. Begin with Fo=alpha*t/Lc^2 and then examine which neglected effect might dominate.
• Option 4: Numerical conduction and convergence. A thermally cycling process vessel wall is modeled with a finite-difference mesh. Run a mesh/time-step comparison and identify whether a spike is physical or numerical. Begin with Fo_grid=alpha*dt/dx^2 ≤ 1/2 (1D explicit) 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 3.1: Cooling time constant

Process situation and calculation

A stainless assembly has mass 4 kg and cp 500 J/(kg K), exposed area 0.4 m² and effective h=25 W/(m² K). Its lumped time constant is tau=mcp/(hA)=200 s. After 200 s, the difference from surrounding temperature has decreased to 1/e≈36.8% of its initial difference.

Physical interpretation and possible error

If the device is insulated or thicker than a low-Biot approximation permits, a single time constant may fit early data but give the wrong center temperature. Report both the fitted time constant and the Biot screening calculation.

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 3.2: Thickness-squared heating penalty

Process situation and calculation

A ceramic part of diffusivity alpha=8×10⁻⁷ m²/s has an indicative diffusion time L²/alpha. For L=0.02 m the scale is 500 s; at L=0.04 m it is 2000 s. Doubling thickness quadruples the diffusion time under a geometrically similar conduction problem.

Physical interpretation and possible error

This timescale does not specify an exact required heating duration. That duration depends on boundary temperatures, acceptable thermal gradients and geometry. It does warn against scaling oven recipes directly with wall thickness.

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 3.3: Center temperature versus wall temperature

Process situation and calculation

A thick product is loaded into a hot-water retort. The outside sensor reports 120°C while a calibrated center probe is still at 97°C. A thermal model that assumes instantaneous uniform temperature would assign the wrong product thermal history.

Physical interpretation and possible error

A center temperature is more appropriate for verifying a minimum-temperature treatment objective, but validating biological safety requires an appropriate product-specific process and qualified procedure. Thermal energy balance alone is not enough to certify treatment efficacy.

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 3.4: Explicit finite-difference stability

Process situation and calculation

Consider one-dimensional transient conduction with alpha=1.0×10⁻⁵ m²/s and a grid step dx=2 mm. The explicit scheme stability limit is dt≤dx²/(2alpha)=0.2 s. A 1 s time step gives Fourier grid number 2.5, far above the limit.

Physical interpretation and possible error

A numerical temperature overshoot in this setup is likely algorithmic, not evidence of a new physical heat-transfer mechanism. Halving dx reduces the stable time step by a factor of four, so grid refinement must be paired with time-step refinement.

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 05 — Lumped-object cooling

Linked unit: 3. Measured output: Predicted object temperature (°C). Equation: T=Ta+(T0−Ta) exp(−hAt/(m cp)).

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: Bi≪0.1, no internal gradient, constant h and cp. 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: time, unit s. Initial value: 120 s. Predicted output: 82.18057 °C. Your result may round on the web page. Confirm the orders of magnitude independently with the printed equation.

At the first selected test value 120 s, the reference model gives 82.18057 °C. At the second test value 480 s, it gives 36.44007 °C. 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 06 — Diffusion penetration length

Linked unit: 3. Measured output: Diffusion length (mm). Equation: ℓ=sqrt(alpha t) ×1000.

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: Order-of-magnitude penetration length, not an exact isotherm or thermal front. 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: time, unit s. Initial value: 300 s. Predicted output: 60 mm. Your result may round on the web page. Confirm the orders of magnitude independently with the printed equation.

At the first selected test value 240.08 s, the reference model gives 53.67457 mm. At the second test value 960.02 s, it gives 107.3324 mm. 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 03 — Assessment: Transient conduction and thermal response

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 3.1: Lumped-capacitance cooling

1. Derivation: Starting with `T(t)=Tinf+(T0−Tinf)exp(−h*A*t/(m*cp))`, 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 small thermocouple bead is moved from hot air to cooling liquid. 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: Compare the sensor response time with the actual bath temperature-change time. State how the result would change for a tenfold increase in one relevant physical length, and whether the original correlation remains valid.

Topic 3.2: Thermal penetration and diffusion time

1. Derivation: Starting with `ell=sqrt(alpha*t); td~L^2/alpha`, 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 batch oven heats thick slabs of different materials. 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: Predict the consequence of doubling thickness and compare a 1D thermal field with an isothermal estimate. State how the result would change for a tenfold increase in one relevant physical length, and whether the original correlation remains valid.

Topic 3.3: Heisler limits and spatial gradients

1. Derivation: Starting with `Fo=alpha*t/Lc^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 cylindrical food can must reach a minimum center temperature. 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: Explain why average temperature cannot establish lethality or a minimum-center-temperature guarantee. State how the result would change for a tenfold increase in one relevant physical length, and whether the original correlation remains valid.

Topic 3.4: Numerical conduction and convergence

1. Derivation: Starting with `Fo_grid=alpha*dt/dx^2 ≤ 1/2 (1D explicit)`, 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 thermally cycling process vessel wall is modeled with a finite-difference mesh. 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: Run a mesh/time-step comparison and identify whether a spike is physical or numerical. 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.