# Effective Moisture Penetration Depth (EMPD) Model[LINK]

Moisture has little effect on heating system performance, but a profound effect on the performance of air conditioning systems. In order to accurately describe building performance during periods when cooling is needed, it is very important to know the moisture conditions of the building. If one assumes that all building moisture is contained in the room air, then one ignores the fact that the materials that bound the room (e.g. wall surfaces, furnishings, linens, etc.) store and release moisture. Thus, to assume that the only moisture that effects cooling system performance is contained in the room air is a false, and it can lead to significant error in the prediction of room moisture conditions and cooling system loads.

The EMPD (Effective Moisture Penetration Depth) model is a simplified, lumped approach to simulate surface moisture adsorption and desorption.

## EMPD Model Description[LINK]

The EMPD concept assumes that a thin layer (δ_{M}) close to the wall surface behaves dynamically and exchanges moisture with the air domain when exposed to cyclic air moisture pulses. For short periods where the cyclic integral of the total moisture adsorption and desorption is near zero (i.e. there is no net moisture storage), the EMPD concept has been shown to be a reasonable approximation of reality (Kerestecioglu et al, 1989). In other words, the following constraint must be met:

∫τ2τ1dUdτdτ=0

where, τ_{2}-τ_{1} denotes the finite time interval over which the equation holds. The EMPD model assumes no spatial distribution of moisture content across the thickness (L) of the solid; rather, a thin layer (δ_{M}) of uniform moisture content (U) is assumed to represent the total moisture content of the solid. This may be mathematically stated as:

∫L0U(x)dx=UδM

For most building materials, the equilibrium moisture sorption isotherm can be defined by the following general equation (Kerestecioglu et al. 1988):

U=aφb+cφd

where

φ≈W∗Wsat∗

and

Wsat∗=1RvρaT∗exp(23.7093−4111T∗−35.45)

Given that U=U(W^{*},T^{*}), the moisture content may be differentiated with respect to time in the following manner:

dudτ=∂U∂W∗dW∗dτ+∂U∂T∗dT∗dτ=ATdW∗dτ−BρdT∗dτ

where A_{T} and B_{ρ} are the isothermal moisture capacity and thermo-gradient coefficient, respectively. From Eqs. , and , they can be expressed as:

AT=abφb+cdφdW∗

and

Bρ=−[1T∗−4111(T∗−35.45)2]∗(abφb+cdφd)

The lumped mass transfer equation for the i-th solid domain may be written as

(AρbδM)idUidτ=hM,iAi(Wr−Wi∗)

Using Eqs. , , and , one obtains the final equation needed for closure moisture transfer at internal surface.

(AiρbδMAT)idWi∗dτ=hM,iAi(Wr−Wi∗)+(AρbδMBρ)idTi∗dτ

The energy equation for the envelope contains the surface temperature and is given by the conduction equation

ρCpdTdτ=∇⋅(k∇T)

with the boundary conditions at interior surface

−k∇T=−qT+hT(T∗−Tr)+λhM(W∗−Wr)

A more detailed account of the numerical solution procedure can be found in Kerestecioglu et al. (1988).

## EMPD Value Determination[LINK]

An effective moisture penetration depth may be determined from either experimental or detailed simulation data by using actual surface areas and moisture vapor diffusivity. An empirical function derived from the detailed simulation may be used to determine the EMPD value (Kerestecioglu et al, 1989):

δM=12.567024−12.21373∗exp(−267.0211∗Dv0.7∗ξ−0.7)

where

ξ=∣∣∣ΔφΔτ∣∣∣

Figure 17 gives the EMPD values to be used for various vapor diffusivities evaluated at different ambient excitations.

## EMPD Nomenclature[LINK]

A= Area [m^{2}]

A_{T}= Isothermal moisture capacity [m^{3}/kg]

Bρ = Thermo-gradient coefficient [kg/kg-K]

C_{p} = Specific heat [J/kg.K]

h_{M} = Convective mass transfer coeff. [kg/m^{2}-s]

h_{T} = Convective heat transfer coeff. [W/m^{2}-K]

k = Thermal conductivity [W/m-K]

L = Length [m]

q“_{T} = Imposed heat flux [W/m^{2}]

R_{v} = Ideal gas constant [461.52 J/kg-K]

T = Temperature [K]

U = Moisture content [kg/kg]

W = Humidity ratio [kg/kg]

### Greek letters[LINK]

δ_{M} = Effective penetration depth for moisture equation [m]

λ = Heat of vaporization [J/kg]

ρ = Density [kg/m^{3}]

τ = Time [s]

φ = Relative humidity [0 to 1]

ξ = Ambient moisture excitation rate [1/h]

### Subscripts and superscripts[LINK]

a= Air

b= Bulk

* = Surface

i = i-th surface

Kerestecioglu, A., Swami, M., Dabir, R., Razzaq, N., and Fairey, P., 1988, “Theoretical and Computational Investigation of Algorithms for Simultaneous Heat and Moisture Transport in Buildings,” FSEC-CR-191-88, Florida Solar Energy Center, Cape Canaveral, FL.

Kerestecioglu, A., M. Swami and A. Kamel, 1989, “Theoretical and Computational Investigation of Simultaneous Heat and Moisture Transfer in Buildings: Effective Penetration Depth Theory.” ASHRAE Winter Meeting, Atlanta, GA.

Kerestecioglu, A., M. V. Swami, P. Brahma, L. Gu, P. Fairey, and S. Chandra, 1989, “FSEC 1.1 User’s Manual,” Florida Solar Energy Center, Cape Canaveral, FL

## Effective Moisture Penetration Depth (EMPD) Model[LINK]

## Overview[LINK]

Moisture has little effect on heating system performance, but a profound effect on the performance of air conditioning systems. In order to accurately describe building performance during periods when cooling is needed, it is very important to know the moisture conditions of the building. If one assumes that all building moisture is contained in the room air, then one ignores the fact that the materials that bound the room (e.g. wall surfaces, furnishings, linens, etc.) store and release moisture. Thus, to assume that the only moisture that effects cooling system performance is contained in the room air is a false, and it can lead to significant error in the prediction of room moisture conditions and cooling system loads.

The EMPD (Effective Moisture Penetration Depth) model is a simplified, lumped approach to simulate surface moisture adsorption and desorption.

## EMPD Model Description[LINK]

The EMPD concept assumes that a thin layer (δ

_{M}) close to the wall surface behaves dynamically and exchanges moisture with the air domain when exposed to cyclic air moisture pulses. For short periods where the cyclic integral of the total moisture adsorption and desorption is near zero (i.e. there is no net moisture storage), the EMPD concept has been shown to be a reasonable approximation of reality (Kerestecioglu et al, 1989). In other words, the following constraint must be met:∫τ2τ1dUdτdτ=0

where, τ

_{2}-τ_{1}denotes the finite time interval over which the equation holds. The EMPD model assumes no spatial distribution of moisture content across the thickness (L) of the solid; rather, a thin layer (δ_{M}) of uniform moisture content (U) is assumed to represent the total moisture content of the solid. This may be mathematically stated as:∫L0U(x)dx=UδM

For most building materials, the equilibrium moisture sorption isotherm can be defined by the following general equation (Kerestecioglu et al. 1988):

U=aφb+cφd

where

φ≈W∗Wsat∗

and

Wsat∗=1RvρaT∗exp(23.7093−4111T∗−35.45)

Given that U=U(W

^{*},T^{*}), the moisture content may be differentiated with respect to time in the following manner:dudτ=∂U∂W∗dW∗dτ+∂U∂T∗dT∗dτ=ATdW∗dτ−BρdT∗dτ

where A

_{T}and B_{ρ}are the isothermal moisture capacity and thermo-gradient coefficient, respectively. From Eqs. , and , they can be expressed as:AT=abφb+cdφdW∗

and

Bρ=−[1T∗−4111(T∗−35.45)2]∗(abφb+cdφd)

The lumped mass transfer equation for the i-th solid domain may be written as

(AρbδM)idUidτ=hM,iAi(Wr−Wi∗)

Using Eqs. , , and , one obtains the final equation needed for closure moisture transfer at internal surface.

(AiρbδMAT)idWi∗dτ=hM,iAi(Wr−Wi∗)+(AρbδMBρ)idTi∗dτ

The energy equation for the envelope contains the surface temperature and is given by the conduction equation

ρCpdTdτ=∇⋅(k∇T)

with the boundary conditions at interior surface

−k∇T=−qT+hT(T∗−Tr)+λhM(W∗−Wr)

A more detailed account of the numerical solution procedure can be found in Kerestecioglu et al. (1988).

## EMPD Value Determination[LINK]

An effective moisture penetration depth may be determined from either experimental or detailed simulation data by using actual surface areas and moisture vapor diffusivity. An empirical function derived from the detailed simulation may be used to determine the EMPD value (Kerestecioglu et al, 1989):

δM=12.567024−12.21373∗exp(−267.0211∗Dv0.7∗ξ−0.7)

where

ξ=∣∣∣ΔφΔτ∣∣∣

Figure 17 gives the EMPD values to be used for various vapor diffusivities evaluated at different ambient excitations.

Limit of Effective Penetration Depth Values for Various Vapor Diffusivities at Different Ambient Excitations.

## EMPD Nomenclature[LINK]

A= Area [m

^{2}]A

_{T}= Isothermal moisture capacity [m^{3}/kg]Bρ = Thermo-gradient coefficient [kg/kg-K]

C

_{p}= Specific heat [J/kg.K]h

_{M}= Convective mass transfer coeff. [kg/m^{2}-s]h

_{T}= Convective heat transfer coeff. [W/m^{2}-K]k = Thermal conductivity [W/m-K]

L = Length [m]

q“

_{T}= Imposed heat flux [W/m^{2}]R

_{v}= Ideal gas constant [461.52 J/kg-K]T = Temperature [K]

U = Moisture content [kg/kg]

W = Humidity ratio [kg/kg]

## Greek letters[LINK]

δ

_{M}= Effective penetration depth for moisture equation [m]λ = Heat of vaporization [J/kg]

ρ = Density [kg/m

^{3}]τ = Time [s]

φ = Relative humidity [0 to 1]

ξ = Ambient moisture excitation rate [1/h]

## Subscripts and superscripts[LINK]

a= Air

b= Bulk

* = Surface

i = i-th surface

## References[LINK]

Kerestecioglu, A., Swami, M., Dabir, R., Razzaq, N., and Fairey, P., 1988, “Theoretical and Computational Investigation of Algorithms for Simultaneous Heat and Moisture Transport in Buildings,” FSEC-CR-191-88, Florida Solar Energy Center, Cape Canaveral, FL.

Kerestecioglu, A., M. Swami and A. Kamel, 1989, “Theoretical and Computational Investigation of Simultaneous Heat and Moisture Transfer in Buildings: Effective Penetration Depth Theory.” ASHRAE Winter Meeting, Atlanta, GA.

Kerestecioglu, A., M. V. Swami, P. Brahma, L. Gu, P. Fairey, and S. Chandra, 1989, “FSEC 1.1 User’s Manual,” Florida Solar Energy Center, Cape Canaveral, FL

Documentation content copyright © 1996-2014 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.