Geothermal Heating
Use this node to account for the contribution of geothermal heat production by adding a source term Qgeo to the right-hand side of the heat equation:
The predefined expression of the heat source uses the Geothermal Density, the Radiogenic heating per unit mass, and a distribution function that can be set.
Geothermal Heat Production
Specify the Radiogenic heating per unit mass qgeo.
The two option buttons in the same group control the distribution function f(zgeo):
If  Uniform distribution is selected, f(zgeo) = 1 and the geothermal heat source is assumed to be independent on depth.
If  Exponential distribution is selected, two additional input fields are displayed underneath. The exponential distribution f(zgeo) = exp(zgeo ⁄ hgeo), is defined by the constant Length scale: hgeo and the variable Depth zgeo. The depth can be, for example, the vertical coordinate direction. The heat source achieves its maximum value where zgeo = 0, typically the top surface of a model.
Geothermal Density
Select the Geothermal density: Solids, Porous media, or User defined:
If  Solids is selected, it calculates the geothermal density based on the volume fraction of solid material
 If Porous media is selected, it calculates the geothermal density based on all mobile and immobile components of the porous medium:
 If User defined is selected, enter a value for the Geothermal density ρgeo.
The Geothermal Heating subnode requires the Subsurface Flow Module. For a detailed overview of the functionality available in each product, visit https://www.comsol.com/products/specifications/
Location in User Interface
Context Menus
If the Porous matrix model is set to Extended under the Physical Model section:
More locations are available if the Heat transfer in porous media check box is selected and Porous matrix model is set to Extended under the Physical Model section. For example:
Ribbon
Physics Tab with Porous Medium selected in the model tree: