Question-and-Answer Resource for the Building Energy Modeling Community
Get s tarted with the Help page
Ask Your Question

Is there a way to extract thermal mass from an energyplus model?

asked 2019-06-15 10:24:20 -0600

nheeren gravatar image

updated 2019-06-19 02:56:21 -0600

I am comparing different building designs and would like to understand the role of thermal mass for their energy performance better. Is there a way to extract the total available thermal mass in a building model or in a zone, e.g. in J/m2K?


After receiving two fantastic answers, I realised that I should have phrased my question a bit more specifically. Technically, heat capacity will not only depend on a material's thickness, density, and heat capacity, but also on dynamic effects, such as amplitude, etc. One rule of thumb is that only the first few centimetres in a surface have relevant thermal inertia. So I was wondering if any kind of measure could be extracted form e+.

Not sure if this update warrants for a new question...

edit retag flag offensive close merge delete

2 Answers

Sort by ยป oldest newest most voted

answered 2019-06-18 16:40:16 -0600

updated 2019-06-19 09:48:49 -0600

To do it programmatically, you could either use OpenStudio and Ruby or EnergyPlus and Python.

OpenStudio and Ruby

Import the IDF to OpenStudio and as @Luis Lara suggested run a Ruby script (or formal OpenStudio Measure) on the OpenStudio model.

EnergyPlus and Python

Here's a lightly tested Python script (using eppy). Depending on your model, you may want to add filters (e.g. Construction:InternalSource, Material:NoMass, etc.).

import os
from eppy.modeleditor import IDF

idd_dir = r'C:\EnergyPlusV9-1-0'
idd_file = r'Energy+.idd'
idd_path = os.path.join(idd_dir, idd_file)

idf_dir = r'C:\WORK\18151 DEN Level 4\18151-Baseline\run'
idf_file = r'in.idf'
idf_path = os.path.join(idf_dir, idf_file)


idf = IDF(idf_path)

total_heat_capacity = 0

surfaces = idf.idfobjects['BuildingSurface:Detailed'.upper()] 

# loop through the surfaces to calculate heat capacity
for surface in surfaces:

    construction_heat_capacity = 0

    construction = idf.getobject('CONSTRUCTION', surface.Construction_Name)

    # get the field names, which in eppy is a list of 
    # ['key', 'Name', 'Outside_Layer', 'Layer_2', ... 'Layer_10']
    construction_field_names = construction.fieldnames

    # loop through the construction's field names (layers), 
    # skipping first and second (key and Name)
    for field_name in construction_field_names[2:]:

        # get the material object from the layer
        attribute = field_name
        material_name = getattr(construction, attribute)
        material = idf.getobject('MATERIAL', material_name)

        # calculate the material's heat capacity and add it to the construction's
        if material is not None:
            thickness = material.Thickness
            density = material.Density
            specific_heat = material.Specific_Heat
            material_heat_capacity = thickness * density * specific_heat
            construction_heat_capacity += material_heat_capacity

    # add the construction's total heat capacity to the model's total
    total_heat_capacity += construction_heat_capacity

print('J/m2*K =', total_heat_capacity)

References from OpenStudio

edit flag offensive delete link more


Thank you so much for this answer! I was wondering if there was as more sophisticated method to do this too. Technically, heat capacity will not only depend on thickness, density, and heat capacity, but also on dynamic effects, such as amplitude, etc. One rule of thumb is that only the first few centimetres in a surface have relevant thermal inertia. So my thought was if any kind of measure could be extracted form e+.

nheeren gravatar image nheeren  ( 2019-06-19 02:53:05 -0600 )edit

Not sure what you mean by "amplitude", but my script above could be easily changed to only consider the first x centimeters of the exterior surfaces.

MatthewSteen gravatar image MatthewSteen  ( 2019-06-19 12:54:03 -0600 )edit

answered 2019-06-18 09:04:42 -0600

updated 2019-06-18 14:21:24 -0600

Your questions is for Energyplus, but if you use OpenStudio you can use the following script to get the total heat capacity for all the surfaces in your model as long as you hard assigned the constructions (eg. did not use space types to assign a construction set):

model_surfaces = model.getPlanarSurfaces
building_mass = 0
model_surfaces.each do |surface|
    if not surface.heatCapacity.empty?
        building_mass += surface.heatCapacity.get
puts building_mass
edit flag offensive delete link more


When I implement this method, I get 0 J/m^2*K (and by extension, also 0 J/K) for all surfaces.

I also tried looping through the Construction Bases.


My code extracts the heat capacity on a "per construction" basis, then each "per construction" heat capacity is multiplied by the net surface area sum of all Surfaces that use that Construction Base (using the .getNetArea() method on the Construction Base). This heat capacity, in J/K, is also always 0.

sashadf1 gravatar image sashadf1  ( 2021-05-06 10:24:26 -0600 )edit

Can you share your material definition?

Luis Lara gravatar image Luis Lara  ( 2021-05-07 14:32:26 -0600 )edit

Here are my material, construction, and surface definitions.

I may be able to calculate the thermal mass and conductivity by scripting to extract from the IDF file itself, but it would be much easier/faster to do it with a measure.

sashadf1 gravatar image sashadf1  ( 2021-05-21 06:17:45 -0600 )edit

Your Answer

Please start posting anonymously - your entry will be published after you log in or create a new account.

Add Answer


Question Tools



Asked: 2019-06-15 10:24:20 -0600

Seen: 434 times

Last updated: Jun 19 '19