modelCarrollRadiation

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

Information

CarrollRadiation

This component links the heat flow with a temperature difference between a wall and the Carroll node, actually the Mean Radiant Temperature (MRT), thanks to the Stefan-Boltzmann law of radiation.
It incoporates the resistance calculation between the surface net flux and radiatity and the resistance induced by the form factor. The detail of the model is available in the theoretical note.
The following assumptions are made:

  • Wall is opaque: Their transmittivty is zero
  • Kirchhoff law is supposed practical: spectral absorptivity is equal to spectral emissivity
  • The surface bodies are grey: They radiative properties are thus independant from the wavelength. It implies to remains within a small range a wavelength
The basic constitutive equation for CarrollRadiation derives :

Where :

  • Q_flow: Heat flow rate from connector 'wall' (e.g., a plate) to connector 'MRT' (the Carrol's Node).
  • add_on is a user paramater to adjust if needed the Heat flow
  • σ is Stefan-Boltzmann's constant
  • A is the radiative surface area of the wall
  • ε is the emissivity
  • Fviewj->MRT is the form factor (or view factor) between the wall and the Carrol's Node

However, the above equation is non-linear. The idea is then to reduce this last equation to an equation of the Newton’s cooling law type. The heat transfer coefficient h_rad from the linearization of the above equation is calculated by a function.

Where :

  • h_rad is the linearized heat transfer coefficient (Newton's coefficient)
  • add_on is a user paramater to adjust if needed the Heat flow to the temperature difference dT
  • dT is the Temperature difference between the wall (port_a.T) and port_b.T

Using a function to calculate h_rad rather than putting directly the relation between the heat flux and the temperature difference without function seems to us to allow a gain of stability. Indeed, the passage by a function seemed to us to make the calculation of the exchange coefficient explicit rather than implicit.

In addition, the transition from biquadratique to a linear relationship between the heat flux and the temperature difference makes the calculation more stable numerically

The difference between the heat flow lost or gained by the solid (the net flux) and the heat flow leaving the surface taking account of the reflectivity is seen as 'thermal resistance' R_th. This thermal resistance is given as input to the function to compute h_rad (See h_rad function description).

Parameters

TypeNameDefaultDescription
Realadd_on1Custom add-on
RealsigmaModelica.Constants.sigmaStefan-Boltzmann constant
Geometrical properties
Modelica.SIunits.AreaA0Wall Surface Area
Radiative properties
Modelica.SIunits.Emissivityeps1Grey wall surface emissivity

Connectors

TypeNameDefaultDescription
HeatPort_aport_a (from Element1D)
HeatPort_bport_b (from Element1D)
Modelica.Blocks.Interfaces.RealInputFview

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
Modelica.SIunits.CoefficientOfHeatTransferh_radHeat transfer coefficient
Modelica.SIunits.EnergyEEnergy passed throught the component