Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 1
平成17-18年度 鳥取県産官学共同研究 研究成果報告
外断熱建築物のライフサイクルコスト低減に関する研究
外断熱建築物のライフサイクルコスト低減に関する研究
∼外断熱工法による断熱性能に関する非定常評価∼
∼外断熱工法による断熱性能に関する非定常評価∼
米子工業高等専門学校 機械工学科
准教授 森田慎一
平成19(2007)年9月17日(木)14:00∼15:00
フジ化成工業株式会社(米子市)
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 2
外断熱工法研究に関する背景
(高耐久性)
外断熱工法は,
耐久性能の向上
および
冷暖房
ランニングコスト抑制
が期待できることから,建
築物のライフサイクルコスト低減できるとされる.
外断熱工法の高耐久性を明らかにした研究は,
コンクリートの中性化を防止を評価した
申ら
(1)(2)
による研究事例がある.
(1)「外 断熱 工法 の透 気 性と躯体 コン クリー トの 中性 化に関 する 研究」 , 申雪 寒,他2名 ,日 本建 築学 会 技術 報告 集 第 18号 ,21- 24,’03.12. (2)「外 断熱 建築 物躯 体 コン クリ ートの中 性化 に関する 研 究」 , 申雪 寒,他1名 ,日 本建 築学 会 技術 報告 集 第 20号 ,59- 62,’04.12.
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 3
外断熱工法研究に関する背景
(高断熱性)
【熱エネルギ的観点による研究】
森ら
(3)の実測調査:札幌市立地,延べ床面積119.88m
2,3階建実建物RC外断熱
住宅「深夜電力を利用した空気熱源ヒート ポンプシステム」躯体蓄熱型暖房シス
テム(床暖房システム)による室温変動,躯体の蓄放熱,熱源機器特性を評価.
菊田ら
(4)の数値解析:森ら
(3)による実測調査による諸係数を評価対象.
(3)「外 断熱 建物における 躯 体 蓄熱を 利 用した暖冷 房に関 する 研究 –躯体 中 心に配 管を した場 合の 蓄放 熱性 状-」 , 森太 郎,他4名 ,空 気調 和・ 衛生 工学 会 北海 道支 部第 38回 学術 講演 論 文集 ,31- 34,’ 04.3. (4)「外 断熱 住宅 の躯 体 蓄熱 型暖 房システムに 関する 研究 」, 菊田 弘輝 ,他4名,日本 建築 学 会環 境系 論文 集 ,589,37-42,’05.3.外断熱と内断熱を直接比較
外断熱と内断熱を直接比較
することによる評価を実施
することによる評価を実施
屋外側に断熱材を配する外断熱が,内断熱よりも断熱性に優れてい
屋外側に断熱材を配する外断熱が,内断熱よりも断熱性に優れてい
ることは明らか.しかし,具体的数値を示す既往研究は見当たらない.
ることは明らか.しかし,具体的数値を示す既往研究は見当たらない.
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory
熱移動(伝熱)の三形態
4T
H
T
m
T
L
λ
1
λ
2
h
A ir
①熱放射
(熱輻射)
②熱伝導
(熱通過・貫流)
③熱伝達
電磁波
(赤外光)
分子振動
分子振動・
物質移動
介在物不要の
エネルギ移動
固体内熱移動
固体と流体
間の熱移動
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory
定常と非定常
5)
/
/
/
1
/(
1
)
(
2 1 1 1x
x
h
h
K
T
T
K
q
H L=Equal=
T
H
T
m
T
L
λ
1
λ
2
外断熱
q
T
H
T
m
T
L
λ
1
λ
2
内断熱
q
L
T
L
T
∞(x)
T
3(x)
T
2(x)
T
1(x)
T
x
T=
0℃
T=
100℃
q
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory
外断熱工法の断熱性能評価
6外断熱工法の
断熱性能有効性評価に
必要なこと
1 .非定常評価
2 .熱伝達,熱伝導による熱移動
(冬期想定による熱放射項の排除)
3 .実験装置による熱伝達率の測定
4 .実験データを用いた非定常シミュレー
ション
外断熱(冬期)
T
H
T
m
T
L
λ
1
λ
2
q
屋外(寒)
屋内(暖)
T
out
T
in
h
out
h
in
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 7
本共同研究の目的
H17∼18年度共同研究
1.非定常伝熱現象として評価
2.外断熱と内断熱外壁の直接比較
3.外断熱,内断熱コンクリート外壁による熱伝達率測定
4.損失熱量比較シミュレーション
H19年度∼
5.実験値を用いた数値シミュレーション
6.空調機器のランニングコスト試算
以上,本報告内容
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 8
A
A-
-A
A’
’断面
断面
内断熱と外断熱の違いとは
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory
実験
外断熱・内断熱壁非定常伝熱実験
9
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory
実験装置の外観写真
10 Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 11
実験装置および方法1
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 12
実験装置および方法2
コンクリート壁と断熱材の
コンクリート壁と断熱材の
一体壁を反転させるこ とにより、
一体壁を反転させるこ とにより、
内断熱と外断熱を測定す る
内断熱と外断熱を測定す る
+
+
q
q
冷風
冷風
熱流束方向を冷風側に
熱流束方向を冷風側に
流れるものをプラスとする
流れるものをプラスとする
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 13
試験片および測定装置
□
□
:
:
熱流束センサ
熱流束センサ
×
×
:
:
赤外線温度センサ
赤外線温度センサ
●
●
:
:
素線径
素線径
φ
φ
0.08 m m
0.08 m m
K
K
型熱電対
型熱電対
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 14
想定季節 冬期
冷風の流速 v
L
=0.52 1.25 1.86 [m/s]
温風の流速 v
H
=0 [m/s]
実験条件
Adiabatic
Adiabatic
material
material
Concrete
Concrete
v
L
v
H
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 15
Adiabatic
Adiabatic
material
material
Concrete
Concrete
t
out
t
t
in
in
v
L
v
H
実験結果(温度差Δtについて)
温度差
Δtは,以下の式により求める.
Δ
T=|t
out
- t
in
|
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 16
0
0.5
1
1.5
2
2.5
3
3.5
4
0
20
40
60
80
100
120
τ[min]
Δ
T
[K
]
-10
0
10
20
30
40
50
60
70
80
q
[W
/m
2]
ΔTi
ΔTo
qi
qo
v
L= 0.52m/s
v
H= 0m/s
実験結果(内・外断熱温度差と熱流束の経時変化v
L=0.52m/s,v
H=0m/s)
Max 68W/m
Max 68W/m
2
2
Max 51W/m
Max 51W/m
2
2
18W/m
18W/m
2
2
熱移動量の少ない外断熱が
熱移動量の少ない外断熱が
外乱の影響を受けにくいといえる
外乱の影響を受けにくいといえる
.
.
0.5W/m
0.5W/m
2
2
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 17
数値計算
実験結果から得た値から,以下の数値計算を行って
内断熱と外断熱を評価していく.
1.実験より得た熱流束値から熱伝達率を計算する.
2.計算した熱伝達率より,積算透過熱量を計算する.
3.一般住宅と10階建てビルを想定したパラメータを与
え,それぞれのコンクリート温度差が30K上昇するまで
の時間を計算する.
4.与えた時間内で変化するコンクリート温度差を計算
する
.
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 18
計算条件(熱伝達率)
実験から得られた熱流束から熱伝達率を以下の式
から求める.
α
:熱伝達率
q:熱流束
Δ
T:温度差
またここでいう
Δ
Tは,試験片表面温度と流動空気との温度
差をとったものを使用する.
)]
/(
[
W
m
2
K
T
q
t
t
in
in
]
[
2
L
K
out
in
t
t
t
T
t
t
out
out
t
t
L
L
●
●
Adiabatic
Adiabatic
material
material
Concrete
Concrete
流動空気温度
流動空気温度
試験片表面温度
試験片表面温度
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 19
計算結果(熱伝達率平均値と各冷風風速v
L)
0
20
40
60
80
100
120
140
0.4
1
1.6
v
L
[m/s]
α[
W
/(
m
2
・K
)]
Inside insulation
Outside insulation
0.52
1.25
1.86
風 速 が増 加する と 共に熱 伝達 率 も増 加傾 向にある 風 速 が増 加する と 共に熱 伝達 率 も増 加傾 向にある平均して約
平均して約
27
27
倍
倍
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 20
計算から得られた熱伝達率から透過熱量を以下の式
から求める.
α
:熱伝達率
q’:透過熱量
Δ
T:温度差
またここでいうΔTは,室内温度を25℃,室外温度を-5℃とした
時と想定し,温度差30Kをとるものとする.
計算条件(積算透過熱量)
min)]
/(
[
'
α・
T
kJ
m
・
2
q
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 21
0
50
100
150
200
250
300
350
400
450
0
20
40
60
80
100
120
τ [min]
q
'
[k
J/
m
2
]
vL=0.52m/saaaa
vL=1.25m/s
vL=1.86m/s
vL=0.52m/s
vL=1.25m/s
vL=1.86m/s
Inside insulation
Outside insulation
v
L=0.52m/s,v
H=0m/s
v
L=1.25m/s,v
H=0m/s
v
L=1.86m/s,v
H=0m/s
v
L=0.52m/s,v
H=0m/s
v
L=1.25m/s,v
H=0m/s
v
L=1.86m/s,v
H=0m/s
計算結果(積算透過熱量)
平均して約
平均して約
14
14
倍
倍
外 断 熱が 透過 させる 熱 量 が少 ない こと か ら外 乱によ る 影響 が少 ない とい える 外 断 熱が 透過 させる 熱 量 が少 ない こと か ら外 乱によ る 影響 が少 ない とい える..Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 22
各建築物のパラメータ
ただし質量は
ただし質量は
,
,
建物に厚さ
建物に厚さ
180mm
180mm
の外壁があるとし
の外壁があるとし
,
,
柱や内壁や窓は
柱や内壁や窓は
,
,
考慮しないものとする
考慮しないものとする
.
.
1 0
1 0
階建てビル
階建てビル
A
A
=5000m
=5000m
2
2
m
m
=8595ton
=8595ton
試験片
試験片
A=0.64m
2
m=230kg
一般住宅
一般住宅
A
A
=380m
=380m
2
2
m
m
=571ton
=571ton
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 23
計算条件(Δtに達する時間)
まず,
Δ
tを積算熱量同様に室内温度を25℃,室外温度を-5℃
とした時と想定し,温度差30Kをとるものとして,以下の式より
Δ
tに達する時間を算出する.
また式の分子には,30K温度が低下するために必要な熱量
分母には,時間変化ごとに逃げる熱量をとる.
Cp:コンクリート比熱
m:各コンクリート質量
q’:経時透過熱量
A:各建物表面積
τ
:時間
A
'
q
t
m
Cp
・
・
・
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 24
0
1000
2000
3000
4000
5000
6000
7000
8000
0.1
1
10
A[m
2
]
100
1000
10000
τ
[h
]
inside insulation 試験片
inside insulation 一般住宅
inside insulation 10階ビル
outside insulation 試験片
outside insulation 一般住宅
outside insualtion 10階ビル
計算結果(熱量に達する時間と表面積との関係)
それぞれ
それぞれ
約
約
19
19
倍
倍
外 断 熱が 内断 熱に 比べ 外気 温が 低く ても室 内温 度を 一定に 保つ ことが できる とい える 外 断 熱が 内断 熱に 比べ 外気 温が 低く ても室 内温 度を 一定に 保つ ことが できる とい える..Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 25
時間変化を1,2,5,10hの4種類をとり,それぞれの時間内で
変化した温度差
Δ
tを求めた式を以下に示す.
W:経時変化熱量 [kJ]
Cp:コンクリート比熱 [kJ/(kg・K)]
m:コンクリート質量 [kg]
ここで表す経時変化熱量は, 各時間変化
パラメータごとの熱量変化量を指す.
計算条件(時間内で変化する温度差)
m
Cp
W
t
・
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 26
0
0 .5
1
1 .5
2
2 .5
3
3 .5
4
0 .1
1
10
1 00
10 00
10 00 0
A[m
2]
Δ
t
1h
2h
5h
10 h
1h
2h
5h
10 h
Inside insulation
O utside insulation
Δtと表面積の関係
試験片
試験片
一般住宅
一般住宅
10
10
階建てビル
階建てビル
少 ない 時 間で内断 熱 は 少 ない 時 間で内断 熱 は,,温度 変化を お こし てしまう が温度 変化を お こし てしまう が 外 断 熱は 外 断 熱は,,温 度 変化 が少 ない ことが わかる温 度 変化 が少 ない ことが わかる..Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory
数値解析
非定常伝熱シミュレーション
27
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory
外壁の非定常伝熱解析
基礎式
温度伝導率: ,比熱: ,発熱密度:
温度: ,時間: ,熱流方向: ,密度:
発熱密度
=0,外壁内外の熱伝達率に実験値を使用し,
差分法によるシミュレーションを実施
28c
Q
x
T
a
t
T
2
2
c
a
c
Q
T
t
x
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory
数値解析結果:内断熱壁内の温度分布変化
29 Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory
数値解析結果:外断熱壁内の温度分布変化
30 屋内壁 温度低下 小さい
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory
コスト比較:長時間
31 Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory
コスト比較:短時間(一晩)
32 Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 33
まとめ
実験的検討結果
1.外断熱は,内断熱に比較して外気の影響を受けにくい.
2.熱伝達率は,内断熱が外断熱の27倍である.
数値シミュレーションによる検討結果
1.温度差30Kの条件下で,定常状態に達するには約2週間を要す
る.
2.外断熱は,内断熱に比較して屋内壁面温度低下を抑制できる.
3.損失熱量を補う電気代は,外断熱壁が内断熱壁の3分の1とな
る結果が得られた.
Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory
今後の検討課題
1.
木造建築物の外断熱工法採用時エネルギコスト検討
2.
ロックセルボードの薄肉化検討
3.
湿度透過および凍害検討
etc
34 Departme nt of Mecha nical E ngineeri ng Mo rita Lab oratory 35