modelThermalResistorRadial

Lumped thermal element transporting heat without storing it

Extends from Modelica.Thermal.HeatTransfer.Interfaces.Element1D (Partial heat transfer element with two HeatPort connectors that does not store energy).

Information

1. Purpose of model

Quasi-static radial thermal resistor for one-dimensional heat conduction in cylindrical geometry between two neighbouring control volumes. The thermal resistance is computed from the cell geometry (inner and outer radius, height) and a constant thermal conductivity. Same idea as the Modelica Standard Library ThermalConductor, but parametrized by the thermal resistance instead of the conductance, which allows a resistance of zero (no temperature difference) in series connections. This is the variant without dispersion; for hydrodynamic dispersion use ThermalResistorRadial_withDispersion.

2. Level of detail, physical effects considered, and physical insight

L1: Quasi-static 1D heat transport in radial direction. No thermal storage within the resistor; the temperature drop is given by the calculated thermal resistance.

3. Limits of validity

- No thermal mass inside the resistor (instantaneous conduction).

- Constant thermal conductivity (no temperature dependency, no dispersion).

- Cylindrical (radially symmetric) geometry.

4. Interfaces

port_a: inner heat port — heat inlet (Heat, inherited from Element1D)

port_b: outer heat port — heat outlet (Heat, inherited from Element1D)

5. Nomenclature

R_r R radial thermal resistance [K/W]
\lambda lambda thermal conductivity of the aquifer/aquitard [W/(m·K)]
r_\mathrm{outer} r outer edge radius of control volume [m]
r_r r_r centre radius of control volume [m]
z z height of control volume [m]
\dot{Q} Q_flow heat flow rate [W] (inherited from Element1D)
\Delta T dT temperature difference across resistor [K] (inherited from Element1D)

6. Governing Equations

Radial thermal resistance for cylindrical geometry:

R = (ln(r_outer) - ln(r_r)) / (lambda * 2*pi*z)

Heat transport (inherited from Element1D):

dT = R * Q_flow

7. Remarks for Usage

Used inside ImpermeableLayer to represent radial heat conduction in the confining layers. A resistance of R = 0 is permitted and yields no temperature difference, which is convenient for series connections. Thermal properties are taken from the Subsurface_Basic (Parameters) record.

8. Validation

Tested in check model "ATES.Validierung.Buoyancy.Val_homogenLayer"

9. References

(no remarks)

10. Version History

Model created by Markus Gillner (markus.gillner@tuhh.de), May 2026

Model revised by Markus Gillner (markus.gillner@tuhh.de), June 2026

- Completed the documentation (10-section template, preferred heading style) and corrected the interface description (heat ports port_a/port_b inherited from Element1D; r is a parameter, not a Real input).

Parameters

TypeNameDefaultDescription
SI.ThermalConductivitylambdaHeat Capacity of the aquifer
SI.Lengthr_rradius from well to center of control volume
SI.Lengthzheight of control volume
SI.LengthrRadius of thermal resistance
SI.ThermalResistanceR(log(r) - log(r_r))/(lambda*2*pi*z)Thermal resistance of material

Connectors

TypeNameDefaultDescription
HeatPort_aport_a (from Element1D)
HeatPort_bport_b (from Element1D)

Components

TypeNameDefaultDescription
SI.HeatFlowRateQ_flow (from Element1D)Heat flow rate from port_a -> port_b
SI.TemperatureDifferencedT (from Element1D)port_a.T - port_b.T