3.3 Transient One-Dimensional Finite Difference Wall Model
3.3 Transient One-Dimensional Finite Difference Wall Model
3.3 Transient One-Dimensional Finite Difference Wall Model
Transient thermal analysis of walls or rooms may be performed with the following objectives:
1. Peak heating/cooling load calculations
2. Calculation of dynamic temperature variation within walls including solar effects, room temperature
swings and condensation on wall interior surfaces.
In the transient finite difference method, we represent each wall layer by one or more sub-layers (regions).
Each region is represented by a node with a thermal capacitance C connected to two thermal resistances, each
equal to half the R-value of the layer, forming a T-section as shown in the figure below for the concrete layer.
For a multilayered wall, an energy balance is applied at each node at regular time intervals to obtain the
temperature of the nodes as a function of time. These equations may be solved with the implicit method as a
set of simultaneous equations or with the explicit method in which we march forward in time from a set of
initial conditions.
The general form of the explicit finite difference formulation corresponding to node i and time interval p is
⎛ Δt ⎞ ⎛ T (j , p) − T (i , p) ⎞
T (i , p + 1) = ⎜―― ⎟ ⋅ ⎜qi + ∑ ―――――― ⎟ + T (i , p)
⎝ Ci ⎠ ⎝ j R (i , j) ⎠
Example: A wall employs the new concept of transparent insulation. It consists of a transparent exterior
layer, an air gap and a concrete thermal storage layer.
3.3_Transient_One-Dimensional_Finite_Difference_Wall_Model.mcdx
Page 2 of 6
2
A≔1 L ≔ 0.12
Concrete properties:
ρ ≔ 2200 ―― density
3
hi ≔ 10 ―――
2
⋅ film
coefficients
ho ≔ 20 ―――
2
⋅
2 2
Rins ≔ 0.3 ⋅ ―― Rgap ≔ 0.3 ⋅ ――
3.3_Transient_One-Dimensional_Finite_Difference_Wall_Model.mcdx
Page 3 of 6
1
Rins + Rgap + ―
ho
Ra ≔ ―――――= 0.65 ――
A
L
Rc ≔ ――
k⋅A
Rc
Rb ≔ ― = 0.018 ――
4
Rc
Rc1 ≔ ― = 0.035 ――
2
Rc
Rc2 ≔ ―
4
1
Ri ≔ ――= 0.1 ――
A ⋅ hi
L 5
C2 ≔ ρ ⋅ c ⋅ ―⋅ A = ⎛⎝1.056 ⋅ 10 ⎞⎠ ――
2
C3 ≔ C2 thermal capacitances
Stability Test
The time step t should be lower than the minimum of the two values in the vector TS.
⎡ C2 C3 ⎤
TS ≔ ――――――――――
⎢ 1 1 1 1 ⎥
⎢ ――― + ―――― + ―――⎥
⎣ Ra + Rb Rc1 Rc1 Rc2 + Ri ⎦
3
Δtcritical ≔ min (TS) Δtcritical = ⎛⎝2.867 ⋅ 10 ⎞⎠
Choose t:
Δt ≔ 1800 ⋅
t ≔ 0 ⋅ , Δt ‥ Steps ⋅ Δt i ≔ 0 ‥ Steps
3.3_Transient_One-Dimensional_Finite_Difference_Wall_Model.mcdx
Page 4 of 6
Assume
frequency based on
w ≔ 2 ⋅ ――― ――
86400 period of one day
⎛ ⎛ ⎞ ⎞ outside
To (t) ≔ ⎜5 ⋅ cos ⎜w ⋅ t + 3 ⋅ ―⎟ − 5⎟
⎝ ⎝ 4⎠ ⎠ temperature
TR ≔ 22 room temperature
⎡ T2 ⎤
⎥≔⎡0⎤
0
⎢
⎢ T31 ⎥ ⎢⎣ 0 ⎥⎦
⎣ ⎦
⎡ ⎛ Teq (i ⋅ Δt) − T2 T3 − T2 ⎞ ⎤
⎢ Δt ⎜ i i i
⎟
⎡ T2 ⎤ ―― ⋅ ―――――― + ――― + T2 ⎥
C2 ⎜ R + R Rc1 ⎟⎠ i
⎢
i+1
⎥≔ ⎢ ⎝ a b ⎥
⎢ T3i + 1 ⎥ ⎢ ⎛ T2 − T3 TR − T3 ⎞ ⎥
⎣ ⎦ ⎢ Δt ⎜ i i i
⎥
―― ⋅ ―――+ ―――⎟ + T3
⎢⎣ C3 ⎜⎝ Rc1 Ri + Rc2 ⎟⎠ i ⎥⎦
T3 − TR
i
T4 ≔ TR + Ri ⋅ ―――
i Rc2 + Ri
Teq (i ⋅ Δt) − T2
i
T1 ≔ T2 + Rb ⋅ ――――――
i i Rb + Ra
3.3_Transient_One-Dimensional_Finite_Difference_Wall_Model.mcdx
Page 5 of 6
-1
-2
-3
-4
-5
To (t) ( )
-6
-7
-8
-9
-10
0 5 10 15 20 25 30 35 40 45 50
t ( )
0.5
0.45
0.4
0.35
0.3
0.25
qsolar (t) ( )
0.2
0.15
0.1
0.05
0 5 10 15 20 25 30 35 40 45 50
t ( )
3.3_Transient_One-Dimensional_Finite_Difference_Wall_Model.mcdx
Page 6 of 6
Variation of nodal temperatures for two days (T4 is room-side wall surface temperature, T1 is the concrete
gap-facing temperature).
49.5
44
38.5
T1 ( ) 33
i
27.5
T2 ( )
i 22
T3 ( ) 16.5
i
11
T4 ( ) 5.5
i
0
-5.5
0 5 10 15 20 25 30 35 40 45 50
i ⋅ Δt ( )
3.3_Transient_One-Dimensional_Finite_Difference_Wall_Model.mcdx