Ground
Heat Transfer Calculations using Site:GroundDomain:Slab[LINK]
In order to simulate heat transfer with horizontal building
surfaces in contact with the ground, a general finite
difference ground model has been implemented. The model
simulates heat transfer with the horizontal building surfaces
through the ground domain. The domain can simulate heat
transfer for slab-in-grade or slab-on-grade scenarios. In all
scenarios, the ground domain must interact with the zone
through an OthersideConditionsModel as the horizontal
surface’s outside boundary condition.
This model is generalized to be able to handle a number of
different slab and insulation configurations. It uses an
implicit finite difference formulation to solve for the ground
temperatures. As a result the simulation is stable for all
timesteps and grid sizes, but an iteration loop must be
employed to converge the temperatures in the domain for each
domain timestep.
Multiple horizontal surfaces can be coupled to each ground
domain object. The model determines which surfaces are coupled
to the ground domain and creates a surface of equivalent
surface area within the ground domain as a representation of
the horizontal surfaces coupled to the ground domain. This
surface then interacts with the ground providing updated other
side conditions model temperatures to the coupled surfaces for
use in their surface heat balance calculations.
At the interface surface, the average surface conduction
heat flux from all surfaces connected to the ground domain is
imposed as a GroundDomain boundary condition at the
Surface/GroundDomain interface cells. Far-field temperatures
are applied as boundary temperature at the GroundDomain sides
and lower surface. The ground temperature profile at the
domain sides and lower surface are taken from Kusuda &
Achenbach 1965. The correlation requires annual ground surface
temperature data.
Ground surface cells are treated as a heat balance, where
long and short wave radiation, conduction, and convection are
considered. Evapotranspiration is also considered. The
evapotranspiration rate is calculated as a moisture loss by
using the Allen et al. (2005) model, and translated into a
heat loss by multiplication with the density and latent heat
of evaporation of water. The evapotranspiration rate is
dependent on the type of vegetation at the surface; the user
can vary the surface vegetation from anywhere between a
concrete surface and a fairly tall grass (about 7”).
Once the ground model has run, the updated cells with zone
surface boundary conditions will update the
OtherSideConditionsModel temperatures which are the used at
the next timestep in the surface heat balance
calculations.
The ground domain is updated at each zone timestep, or
hourly as specified by the user. For situations when the
ground domain is updated at each timestep, the domain is
simulated by applying the surface heat flux boundary
conditions from the previous timestep and calculating a new
OthersideConditionsModel temperature. At this point, the
surface heat balance algorithms can then take the new outside
surface temperatures to update their surface heat fluxes. For
situations when the user has elected to have the domain update
on an hourly basis, the surface heat balance for each coupled
surface is aggregated and passed to the domain as an average
surface heat flux from the previous hour, which will then
update the outside surface temperatures for the surface heat
balance’s next iteration.
Both in-grade and on-grade scenarios are simulated with the
GroundDomain object. The key difference being that for
in-grade situations, the slab and horizontal insulation are
simulated by the ground domain, whereas for the on-grade
situations the slab and horizontal insulation must be included
in the floor construction object. All possible insulation/slab
configurations are seen in Table 18.
Table 18. Possible insulation/slab configurations for
Site:GroundDomain model.
Configuration
Situations
Vert. Ins.
Horiz. Ins (Full)
Horiz. Ins (Perimeter)
In-grade
1
X
–
–
In-grade
2
X
X
–
In-grade
3
X
–
X
In-grade
4
–
–
–
In-grade
5
–
X
–
In-grade
6
–
–
X
On-grade
7
X
–
–
On-grade
8*
X
X
–
On-grade
9
–
–
–
On-grade
10*
–
X
–
* Horizontal insulation must be included in the floor
construction
For the slab-in-grade scenarios, a thin surface layer must
be included in the floor construction. This can be a very thin
layer of the slab or other floor covering materials above the
slab. This provides a zone boundary condition for the
GroundDomain while still allowing 3-dimensional heat transfer
effects to occur within the slab.
Kusuda, T. & Achenbach, P. 1965. Earth Temperature and
Thermal Diffusivity at Selected Stations in the United States,
ASHRAE Transactions 71(1): 61-75.
Allen, R.G., Walter, I.A., Elliott, R.L., Howell, T.A.,
Itenfisu, D., Jensen, M.E., Snyder, R.L. 2005. The ASCE
Standardized Reference Evapotranspiration Equation. Reston,
VA: American Society of Civil Engineers.
Ground
Heat Transfer Calculations using
Site:GroundDomain:Basement[LINK]
Basement heat transfer calculations are performed by the
same model used by Site:GroundDomain:Slab to simulate
zone-ground coupled surfaces. To simulate basements, the
previous model is extended to situations where surfaces are
below ground.
This model is generalized to be able to handle several
basement surface and insulation configurations. An implicit
finite difference formulation is used to solve for the ground
temperatures. As a result the simulation is stable for all
timesteps and grid sizes, but an iteration loop must be
employed to converge the temperatures in the domain for each
domain timestep.
Multiple basement zones can be coupled to each basement
domain object. The model determines which floor surfaces are
coupled to the basement domain and creates a floor surface of
equivalent surface area within the basement domain as a
representation of the horizontal surfaces coupled to the
domain. The user defined aspect ratio, which is the ratio of
basement width to length, along with the total surface area of
all surfaces connected to this basement floor will define the
size and shape of the basement ground domain.
Vertical basement depth below the ground surface of the
domain is user specified. For situations when the ground
surface is lower than the first above ground level a separate
surface should be employed to bridge the space between the
ground domain, and the first above-ground level.
Basement wall and floor surfaces are created normally using
the BuildingSurface:Detailed
objects, with the outside boundary condition being the
OtherSideConditionsModel used for the walls and floor of the
basement domain. The interface between the basement walls and
floors occurs at the outside surface of the walls and floor.
Outside underground insulation is simulated by the ground
domain, and therefore should not be included in the basement
wall and floor construction objects.
At the zone interface surfaces, the average surface
conduction heat flux from all surfaces connected to the ground
domain is imposed as a boundary condition at the
Surface/Domain interface cells. Far-field temperatures are
applied as boundary temperature at the domain sides and lower
surface. The ground temperature profile at the domain sides
and lower surface are taken from Kusuda & Achenbach 1965.
The correlation requires annual ground surface temperature
data. Ground surface temperature data can be determined by
using the CalcSoilSurfTemp preprocessor; or, if the user has
average monthly ground surface temperature data available,
they can use the Site:GroundTemperature:Shallow
object to provide the model monthly ground surface
temperatures. From that, the model can then determine the
ground temperature model parameters for the Kusuda &
Achenbach (1965) ground temperature model.
Ground surface cells are treated as a heat balance, where
long and short wave radiation, conduction, and convection are
considered. Evapotranspiration is also considered. The
evapotranspiration rate is calculated as a moisture loss by
using the Allen et al. (2005) model, and translated into a
heat loss by multiplication with the density and latent heat
of evaporation of water. The evapotranspiration rate is
dependent on the type of vegetation at the surface; the user
can vary the surface vegetation from anywhere between a
concrete surface and a fairly tall grass (about 7”). The model
can be sensitive to variations in this parameter, especially
in hot dry climates.
Once the basement model has run, the updated cells with
zone surface boundary conditions will update the
OtherSideConditionsModel temperatures which are then used at
the next timestep in the surface heat balance
calculations.
The ground domain is updated at each zone timestep, or
hourly as specified by the user. For situations when the
ground domain is updated at each timestep, the domain is
simulated by applying the surface heat flux boundary
conditions from the previous timestep and calculating a new
OthersideConditionsModel temperature. At this point, the
surface heat balance algorithms can then take the new outside
surface temperatures to update their surface heat fluxes. For
situations when the user has elected to have the domain update
on an hourly basis, the surface heat flux for each coupled
surface is aggregated and passed to the domain as an average
surface heat flux from the previous hour, which will then
update the outside surface temperatures for the surface heat
balance’s next iteration.
The basement floor and walls are simulated by standard
surface heat balance models within EnergyPlus (CTF, finite
difference). The horizontal and vertical insulation are
simulated by the ground domain finite difference solver. This
is shown in the Figure below
where the dotted red line defines the OtherSideConditionsModel
interface separating the two solution domains. This
methodology applies to the basement floor and walls.
Kusuda, T. & Achenbach, P. 1965. ‘Earth Temperature and
Thermal Diffusivity at Selected Stations in the United
States’, ASHRAE Transactions 71(1): 61–75.
Allen, R.G., Walter, I.A., Elliott, R.L., Howell, T.A.,
Itenfisu, D., Jensen, M.E., Snyder, R.L. 2005. The ASCE
standardized reference evapotranspiration equation. Reston,
VA: American Society of Civil Engineers.
Ground Heat Transfer Calculations using Site:GroundDomain:Slab[LINK]
In order to simulate heat transfer with horizontal building surfaces in contact with the ground, a general finite difference ground model has been implemented. The model simulates heat transfer with the horizontal building surfaces through the ground domain. The domain can simulate heat transfer for slab-in-grade or slab-on-grade scenarios. In all scenarios, the ground domain must interact with the zone through an OthersideConditionsModel as the horizontal surface’s outside boundary condition.
Approach[LINK]
This model is generalized to be able to handle a number of different slab and insulation configurations. It uses an implicit finite difference formulation to solve for the ground temperatures. As a result the simulation is stable for all timesteps and grid sizes, but an iteration loop must be employed to converge the temperatures in the domain for each domain timestep.
Multiple horizontal surfaces can be coupled to each ground domain object. The model determines which surfaces are coupled to the ground domain and creates a surface of equivalent surface area within the ground domain as a representation of the horizontal surfaces coupled to the ground domain. This surface then interacts with the ground providing updated other side conditions model temperatures to the coupled surfaces for use in their surface heat balance calculations.
Boudary Conditions[LINK]
At the interface surface, the average surface conduction heat flux from all surfaces connected to the ground domain is imposed as a GroundDomain boundary condition at the Surface/GroundDomain interface cells. Far-field temperatures are applied as boundary temperature at the GroundDomain sides and lower surface. The ground temperature profile at the domain sides and lower surface are taken from Kusuda & Achenbach 1965. The correlation requires annual ground surface temperature data.
Ground surface cells are treated as a heat balance, where long and short wave radiation, conduction, and convection are considered. Evapotranspiration is also considered. The evapotranspiration rate is calculated as a moisture loss by using the Allen et al. (2005) model, and translated into a heat loss by multiplication with the density and latent heat of evaporation of water. The evapotranspiration rate is dependent on the type of vegetation at the surface; the user can vary the surface vegetation from anywhere between a concrete surface and a fairly tall grass (about 7”).
Once the ground model has run, the updated cells with zone surface boundary conditions will update the OtherSideConditionsModel temperatures which are the used at the next timestep in the surface heat balance calculations.
Simulation Methodology[LINK]
The ground domain is updated at each zone timestep, or hourly as specified by the user. For situations when the ground domain is updated at each timestep, the domain is simulated by applying the surface heat flux boundary conditions from the previous timestep and calculating a new OthersideConditionsModel temperature. At this point, the surface heat balance algorithms can then take the new outside surface temperatures to update their surface heat fluxes. For situations when the user has elected to have the domain update on an hourly basis, the surface heat balance for each coupled surface is aggregated and passed to the domain as an average surface heat flux from the previous hour, which will then update the outside surface temperatures for the surface heat balance’s next iteration.
Both in-grade and on-grade scenarios are simulated with the GroundDomain object. The key difference being that for in-grade situations, the slab and horizontal insulation are simulated by the ground domain, whereas for the on-grade situations the slab and horizontal insulation must be included in the floor construction object. All possible insulation/slab configurations are seen in Table 18.
* Horizontal insulation must be included in the floor construction
For the slab-in-grade scenarios, a thin surface layer must be included in the floor construction. This can be a very thin layer of the slab or other floor covering materials above the slab. This provides a zone boundary condition for the GroundDomain while still allowing 3-dimensional heat transfer effects to occur within the slab.
References[LINK]
Kusuda, T. & Achenbach, P. 1965. Earth Temperature and Thermal Diffusivity at Selected Stations in the United States, ASHRAE Transactions 71(1): 61-75.
Allen, R.G., Walter, I.A., Elliott, R.L., Howell, T.A., Itenfisu, D., Jensen, M.E., Snyder, R.L. 2005. The ASCE Standardized Reference Evapotranspiration Equation. Reston, VA: American Society of Civil Engineers.
Ground Heat Transfer Calculations using Site:GroundDomain:Basement[LINK]
Basement heat transfer calculations are performed by the same model used by Site:GroundDomain:Slab to simulate zone-ground coupled surfaces. To simulate basements, the previous model is extended to situations where surfaces are below ground.
Approach[LINK]
This model is generalized to be able to handle several basement surface and insulation configurations. An implicit finite difference formulation is used to solve for the ground temperatures. As a result the simulation is stable for all timesteps and grid sizes, but an iteration loop must be employed to converge the temperatures in the domain for each domain timestep.
Multiple basement zones can be coupled to each basement domain object. The model determines which floor surfaces are coupled to the basement domain and creates a floor surface of equivalent surface area within the basement domain as a representation of the horizontal surfaces coupled to the domain. The user defined aspect ratio, which is the ratio of basement width to length, along with the total surface area of all surfaces connected to this basement floor will define the size and shape of the basement ground domain.
Vertical basement depth below the ground surface of the domain is user specified. For situations when the ground surface is lower than the first above ground level a separate surface should be employed to bridge the space between the ground domain, and the first above-ground level.
Basement wall and floor surfaces are created normally using the BuildingSurface:Detailed objects, with the outside boundary condition being the OtherSideConditionsModel used for the walls and floor of the basement domain. The interface between the basement walls and floors occurs at the outside surface of the walls and floor. Outside underground insulation is simulated by the ground domain, and therefore should not be included in the basement wall and floor construction objects.
Boundary Conditions[LINK]
At the zone interface surfaces, the average surface conduction heat flux from all surfaces connected to the ground domain is imposed as a boundary condition at the Surface/Domain interface cells. Far-field temperatures are applied as boundary temperature at the domain sides and lower surface. The ground temperature profile at the domain sides and lower surface are taken from Kusuda & Achenbach 1965. The correlation requires annual ground surface temperature data. Ground surface temperature data can be determined by using the CalcSoilSurfTemp preprocessor; or, if the user has average monthly ground surface temperature data available, they can use the Site:GroundTemperature:Shallow object to provide the model monthly ground surface temperatures. From that, the model can then determine the ground temperature model parameters for the Kusuda & Achenbach (1965) ground temperature model.
Ground surface cells are treated as a heat balance, where long and short wave radiation, conduction, and convection are considered. Evapotranspiration is also considered. The evapotranspiration rate is calculated as a moisture loss by using the Allen et al. (2005) model, and translated into a heat loss by multiplication with the density and latent heat of evaporation of water. The evapotranspiration rate is dependent on the type of vegetation at the surface; the user can vary the surface vegetation from anywhere between a concrete surface and a fairly tall grass (about 7”). The model can be sensitive to variations in this parameter, especially in hot dry climates.
Once the basement model has run, the updated cells with zone surface boundary conditions will update the OtherSideConditionsModel temperatures which are then used at the next timestep in the surface heat balance calculations.
Simulation Methodology[LINK]
The ground domain is updated at each zone timestep, or hourly as specified by the user. For situations when the ground domain is updated at each timestep, the domain is simulated by applying the surface heat flux boundary conditions from the previous timestep and calculating a new OthersideConditionsModel temperature. At this point, the surface heat balance algorithms can then take the new outside surface temperatures to update their surface heat fluxes. For situations when the user has elected to have the domain update on an hourly basis, the surface heat flux for each coupled surface is aggregated and passed to the domain as an average surface heat flux from the previous hour, which will then update the outside surface temperatures for the surface heat balance’s next iteration.
The basement floor and walls are simulated by standard surface heat balance models within EnergyPlus (CTF, finite difference). The horizontal and vertical insulation are simulated by the ground domain finite difference solver. This is shown in the Figure below where the dotted red line defines the OtherSideConditionsModel interface separating the two solution domains. This methodology applies to the basement floor and walls.
References[LINK]
Kusuda, T. & Achenbach, P. 1965. ‘Earth Temperature and Thermal Diffusivity at Selected Stations in the United States’, ASHRAE Transactions 71(1): 61–75.
Allen, R.G., Walter, I.A., Elliott, R.L., Howell, T.A., Itenfisu, D., Jensen, M.E., Snyder, R.L. 2005. The ASCE standardized reference evapotranspiration equation. Reston, VA: American Society of Civil Engineers.
Documentation content copyright © 1996-2026 The Board of Trustees of the University of Illinois and the Regents of the University of California through the Ernest Orlando Lawrence Berkeley National Laboratory. All rights reserved. EnergyPlus is a trademark of the US Department of Energy.
This documentation is made available under the EnergyPlus Open Source License v1.0.