JP2012092647A - Ventilation amount and temperature prediction system of building - Google Patents

Ventilation amount and temperature prediction system of building Download PDF

Info

Publication number
JP2012092647A
JP2012092647A JP2012008166A JP2012008166A JP2012092647A JP 2012092647 A JP2012092647 A JP 2012092647A JP 2012008166 A JP2012008166 A JP 2012008166A JP 2012008166 A JP2012008166 A JP 2012008166A JP 2012092647 A JP2012092647 A JP 2012092647A
Authority
JP
Japan
Prior art keywords
temperature
ventilation
room
calculation
air
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2012008166A
Other languages
Japanese (ja)
Other versions
JP5383832B2 (en
Inventor
Yosuke Chiba
陽輔 千葉
Masayuki Hisayoshi
将之 久芳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Asahi Kasei Homes Corp
Original Assignee
Asahi Kasei Homes Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Asahi Kasei Homes Corp filed Critical Asahi Kasei Homes Corp
Priority to JP2012008166A priority Critical patent/JP5383832B2/en
Publication of JP2012092647A publication Critical patent/JP2012092647A/en
Application granted granted Critical
Publication of JP5383832B2 publication Critical patent/JP5383832B2/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Landscapes

  • Building Environments (AREA)

Abstract

PROBLEM TO BE SOLVED: To provide a system for highly accurately predicting a ventilation condition and a temperature condition inside a building and accurately predicting the condition of indoor air quality together by considering the influence of heat transfer and heat generation in ventilation network calculation.SOLUTION: The system comprises house-by-house data, input means, calculation means and output means. The calculation means considers buoyancy by the difference of the density of air depending on temperature, obtains the ventilation amount of the air in each room, calculates the ventilation amount for each opening by maintaining the total sum of the inflow and outflow of the air in each room at 0, considers the transfer of heat accompanying the air by ventilation, and calculates the temperature of each room by the in/out amount of the heat of the respective rooms with each other. When considering the buoyancy by the difference of the density of the air depending on the temperature in the calculation of the ventilation amount, the temperature obtained from the calculation of the temperature of each room is used for the temperature, and by alternately calculating the ventilation amount and the temperature for multiple times and converging them every time at every prescribed time interval, the ventilation amount and the temperature are calculated.

Description

本発明は、建物内部の換気状況及び熱の移動による温度状況を計算により予測するための予測システムに関するものである。   The present invention relates to a prediction system for predicting a ventilation situation inside a building and a temperature situation due to heat transfer by calculation.

従来において、建物内部の換気状況を予測するには、換気回路網計算が使用されている。図10に従来の換気回路網計算の概念図を示す。図10(a)は建物の間取りのレイアウト図であり、各室を1点としてネットワークモデルとしたものを図10(b)に示している。   Conventionally, ventilation network calculation is used to predict the ventilation situation inside a building. Fig. 10 shows a conceptual diagram of conventional ventilation network calculation. FIG. 10 (a) is a layout diagram of a building, and FIG. 10 (b) shows a network model with each room as one point.

図においてネットワークは室と通気経路のみからなり、換気量は外気の温度、風速と室内各部温度を仮定し、各室間の圧力差と通気経路の抵抗から、1室ごとに算出する。各室の圧力差は外気風、室間又は内外温度差によって生じる。ここで換気回路網計算は、内外温度差の効果を簡易に取り扱うため、外気温度と室内の代表的温度一定で計算を行うのが一般的である。   In the figure, the network consists of only a room and a ventilation path, and the ventilation volume is calculated for each room from the pressure difference between the rooms and the resistance of the ventilation path, assuming the temperature of the outside air, the wind speed, and the temperature of each part of the room. The pressure difference between the chambers is caused by the outside air flow, the chamber-to-chamber or the temperature difference between the inside and outside. Here, the ventilation network calculation is generally performed with the outside air temperature and the typical indoor temperature constant in order to easily handle the effect of the temperature difference between the inside and outside.

また同様に、建物内部の温度状況を算出する熱負荷計算プログラムも従来から使用されている。かかる計算にあって室内気流が熱の伝達要素として考慮されるが、初期条件として設定しているに過ぎず、換気状況の把握をすることはできない。   Similarly, a heat load calculation program for calculating a temperature condition inside a building has been used conventionally. In this calculation, the indoor airflow is considered as a heat transfer element, but it is only set as an initial condition, and it is impossible to grasp the ventilation situation.

しかし、建物の中では暖房機、電気器具、日照などの発熱により空気温が上昇し、室間の温度差から気流を生じるなど、換気状況には熱的な影響が大きいにもかかわらず、そのような建物内部の温度分布や発熱などの熱的影響を考慮したものがない。   However, in buildings, the air temperature rises due to heat generated by heaters, appliances, and sunshine, and airflow is generated due to the temperature difference between the rooms. There is no such thing that considers thermal effects such as temperature distribution and heat generation inside the building.

また、より詳細な計算を行おうとした場合、室内の換気状況と温度状況は本来相互に影響するため、両者を関連づけながら解くことが望ましい。しかし従来において換気量と温度の計算を同時に予測することは行われておらず、換気状況を予測するときは温度状況を条件として与え、逆に温度状況を予測するときは換気状況を条件として与えていた。   In addition, when a more detailed calculation is to be performed, the indoor ventilation and temperature conditions inherently affect each other. However, the calculation of ventilation volume and temperature has not been predicted at the same time in the past, and when predicting the ventilation situation, the temperature situation is given as a condition, and conversely, when predicting the temperature situation, the ventilation situation is given as a condition. It was.

特に、換気量の比較的小さい24時間計画換気の状況を計算する場合には、上記のような建物内部の発熱などによる温度ムラは換気の動力として無視できないレベルである。なお、窓開放など通風による掃気であれば、換気量が大きく、内部発熱が問題となる場合は少ない。   In particular, when calculating the status of 24-hour planned ventilation with a relatively small ventilation amount, the temperature unevenness due to heat generation inside the building as described above is a level that cannot be ignored as the power of ventilation. In the case of scavenging by ventilation such as opening a window, the amount of ventilation is large and there are few cases where internal heat generation becomes a problem.

そこで本発明は、換気回路網計算に熱移動、熱発生の影響を考慮することにより、建物内の換気状況及び温度状況を高い精度で予測し、合わせて室内空気質の状況を的確に予測するシステムを提供することを目的としている。   Therefore, the present invention predicts the ventilation and temperature conditions in the building with high accuracy by taking into account the effects of heat transfer and heat generation in the calculation of the ventilation network, and also accurately predicts the indoor air quality. The purpose is to provide a system.

上記課題を解決するために、本発明にかかる建物の温度及び換気予測システムの代表的な構成は、任意の時間の、建物の各室の温度及び室間若しくは建物内外間の換気量を予測するシステムであって、各室の形状及び壁、床、開口部など室相互の接続関係を含む邸別データと、風速及び外気温度を含む気象データとを入力する入力手段と、所定時間間隔毎に各室毎の温度および換気量を計算する計算手段と、計算結果を出力する出力手段とを有し、前記計算手段は温度による空気の密度の違いによる浮力を考慮して、各室における空気の換気量を求め、かつ各室における空気の流出入の総和を0に維持することにより各開口部毎の換気量を計算し、換気による空気に伴う熱の移動を考慮して各室相互の熱の入出量により各室毎の温度を計算し、前記換気量の計算にて前記温度による空気の密度の違いによる浮力を考慮するにつき、当該温度は、前記各室毎の温度の計算から得られる温度を用いることとし、所定時間間隔毎に、毎回、上記換気量と温度とを交互に複数回計算して収束させることにより、換気量及び温度を計算することを特徴とする。 In order to solve the above-mentioned problems, a typical configuration of a building temperature and ventilation prediction system according to the present invention predicts the temperature of each room of a building and the ventilation amount between rooms or inside and outside the building at an arbitrary time. An input means for inputting house-specific data including the shape of each room and connection between rooms such as walls, floors, openings, etc., and weather data including wind speed and outside air temperature, and at predetermined time intervals a calculation means for calculating the temperature and ventilation of each chambers, and output means for outputting a calculation result, the calculation means, taking into account the buoyancy due to a difference in air density due to temperature, air in each chamber The ventilation volume of each opening is calculated by calculating the ventilation volume of each room and maintaining the sum of the inflow and outflow of air in each room to be 0. Calculate temperature for each room based on heat input / output , Per consideration buoyancy due to a difference in air density due to the temperature at the ventilation calculation, the temperature will be referred to with the temperature obtained from the calculation of the temperature of the respective chambers, at predetermined time intervals, Ventilation volume and temperature are calculated each time by calculating and converging the ventilation volume and temperature alternately several times.

また本発明に係る建物の温度及び換気予測システムの他の代表的な構成は、任意の時間の、建物の各室の温度及び室間若しくは建物内外間の換気量を予測するシステムであって、
各室の形状及び開口部を含む邸別データと、風速及び外気温度を含む気象データとを入力する入力手段と、
所定時間間隔毎に各室毎の温度および換気量を計算する計算手段と、計算結果を出力する出力手段とを有し、前記計算手段は温度による空気の密度の違いによる浮力を考慮して各室における空気の換気量を求め、かつ各室における空気の流出入の総和を0に維持することにより各開口部毎の換気量を計算し、換気による空気に伴う熱の移動を考慮して各室相互の熱の入出量により各室毎の温度を計算し、前記換気量の計算にて前記温度による空気の密度の違いによる浮力を考慮するにつき、当該温度は、前記各室毎の温度の計算から得られる温度を用いることとし、所定時間間隔毎に、前時間で算出した温度を用いて換気量を先に計算し、その換気量の計算結果を用いて温度を計算することを特徴とする。
In addition, another typical configuration of the building temperature and ventilation prediction system according to the present invention is a system that predicts the temperature of each room of the building and the ventilation amount between rooms or inside and outside the building at an arbitrary time,
Input means for inputting house-specific data including the shape and opening of each room, and weather data including wind speed and outside air temperature;
A calculation means for calculating the temperature and ventilation of each chambers at predetermined time intervals, and output means for outputting a calculation result, the calculation means, taking into account the buoyancy due to a difference in air density due to temperature , determine the amount of ventilation air in each chamber, and the outflow of the incoming sum of the aerial to calculate the ventilation of the respective openings by maintaining the 0 in each chamber, considering heat transfer due to the air by ventilation The temperature for each room is calculated based on the amount of heat input / output between the rooms, and the buoyancy due to the difference in air density due to the temperature is taken into account in the ventilation volume calculation. Suppose that the temperature obtained from the calculation of the temperature is used, and that the ventilation volume is calculated first using the temperature calculated in the previous time every predetermined time interval, and the temperature is calculated using the calculation result of the ventilation volume. Features.

上記説明した如く、本発明に係る建物の換気量及び温度予測システムは、換気回路網計算に熱移動、熱発生の影響を考慮することにより、建物内の換気状況及び温度状況を高い精度で予測し、合わせて室内空気質の状況を的確に予測することができる。   As described above, the ventilation amount and temperature prediction system for a building according to the present invention predicts the ventilation and temperature conditions in the building with high accuracy by taking into account the effects of heat transfer and heat generation in the ventilation network calculation. In addition, the indoor air quality can be accurately predicted.

本願にかかる予測システムの概念図である。It is a conceptual diagram of the prediction system concerning this application. 第一実施形態にかかる熱及び換気回路網計算のフローチャートである。It is a flowchart of the heat and ventilation network calculation concerning 1st embodiment. 熱的に厚い部位周辺の変数定義を説明する図である。It is a figure explaining the variable definition around a thick part thermally. 土間床のモデルを説明する図である。It is a figure explaining the model of a dirt floor. 第二実施形態にかかる熱及び換気回路網計算のフローチャートである。It is a flowchart of the heat and ventilation network calculation concerning 2nd embodiment. 建物内の標準状態における温度分布の計算例を示す図である。It is a figure which shows the example of calculation of the temperature distribution in the standard state in a building. 建物内の標準状態における換気回数分布の計算例を示す図である。It is a figure which shows the example of calculation of the ventilation frequency distribution in the standard state in a building. すきま面積が拡大した場合の換気回数分布の計算例を示す図である。It is a figure which shows the example of calculation of the ventilation frequency distribution when a clearance area expands. すきま面積が拡大した場合の温度分布の計算例を示す図である。It is a figure which shows the example of calculation of the temperature distribution when a clearance area expands. 従来の換気回路網計算の概念図である。It is a conceptual diagram of the conventional ventilation network calculation.

[第一実施形態]
本発明に係る建物の換気量及び温度予測システムの第一実施形態について、図を用いて説明する。本実施形態にかかる建物の換気量及び温度予測システムは具体的にはコンピュータ上で動作するソフトウェアであり、入力手段とはキーボードや入出力装置、ネットワーク回線などを意味し、計算手段とはCPU及びRAMを意味し、また出力手段とはモニタやプリンタ又は記録装置などを意味している。
[First embodiment]
A first embodiment of a building ventilation and temperature prediction system according to the present invention will be described with reference to the drawings. The building ventilation amount and temperature prediction system according to the present embodiment is specifically software that operates on a computer, the input means means a keyboard, an input / output device, a network line, etc., and the calculation means means a CPU and RAM means output means means a monitor, a printer or a recording device.

(全体構成)
図1は熱及び換気回路網計算の概念を示す図であって、図1(a)は建物の間取りのレイアウト図であり、各室を1点としてネットワークモデルとしたものを図1(b)に示している。通気経路においては室内外及び各室間の扉や隙間などの開口部、及び換気扇などによる強制換気の通路を考慮する。熱の移動経路においては壁や床天井等から伝達される熱量、及び換気に乗って移動する熱量を考慮する。また図に示すように、各部屋には定量的な発熱源若しくは吸熱源、及び湿度や化学物質の発生源も考慮することにより、室内空気質の予測も行う。
(overall structure)
Fig. 1 is a diagram showing the concept of heat and ventilation network calculation. Fig. 1 (a) is a layout diagram of the building. Fig. 1 (b) shows a network model with each room as one point. It shows. In the ventilation path, openings such as doors and gaps between the rooms and between the rooms, and a passage for forced ventilation by a ventilation fan are considered. In the heat transfer path, the amount of heat transferred from the walls, floor and ceiling, etc. and the amount of heat transferred by ventilation are considered. Further, as shown in the figure, indoor air quality is also predicted by taking into consideration the quantitative heat generation source or heat absorption source and the generation source of humidity and chemical substances in each room.

図2に熱及び換気回路網計算のフローチャートを示している。なお、具体的な個々の計算については後述する。図2に示すように、対象とする建物の邸別データと気象データを入力手段より読み込み(S1、S2)、これらを用いて以下の計算を行う。邸別データとは図1(a)に示した如く建物の間取りや各室の形状及び壁、床、開口部など室相互の接続関係、更に建物を構成する部材若しくは部位の熱容量、室内に配置された任意の発熱源を含むものである。気象データとは市販のデータを用いることができ、風速及び外気温度、外気湿度、日射又は放射量の情報を利用する。なお風速のデータがなくとも、無風状態として計算することは可能である。   FIG. 2 shows a flowchart of heat and ventilation network calculation. Specific individual calculations will be described later. As shown in FIG. 2, the house-specific data and weather data of the target building are read from the input means (S1, S2), and the following calculation is performed using them. As shown in Fig. 1 (a), the data for each residence is the layout of the building, the shape of each room, the connection between rooms such as walls, floors, openings, etc., the heat capacity of the members or parts that make up the building, and the indoor arrangement It includes any heat source that is Commercial data can be used as meteorological data, and information on wind speed and outside air temperature, outside air humidity, solar radiation, or radiation amount is used. Even if there is no wind speed data, it is possible to calculate as no wind condition.

邸別データと気象データを読み込んだ後に、まず計算すべき時刻の気象条件を取得し(S3)、日射から建物に与えられる放射量を計算する(S4)。次に換気回路網計算を行って各室毎の換気量を計算し、各室間の換気量を算出する(S5)。次にS5で算出された換気量を用いて、熱回路網計算による各室毎の温度(S6)を計算する。そしてS6で算出した温度を用いて、換気回路網計算(S5)に戻り、さらに換気量を算出する。これらの換気回路網計算(S5)と熱回路網計算(S6)は、収束するまで繰り返し計算を行う。それから湿度等回路網計算による湿度や化学物質の濃度(S7)を算出し、あわせて得られた温度から熱負荷計算をも行った後に(S8)、算出結果を出力する(S9)。   After reading the house data and the meteorological data, the meteorological conditions at the time to be calculated are first acquired (S3), and the amount of radiation given to the building from solar radiation is calculated (S4). Next, the ventilation network calculation is performed to calculate the ventilation amount for each room, and the ventilation amount between the rooms is calculated (S5). Next, the temperature (S6) for each room by the thermal network calculation is calculated using the ventilation volume calculated in S5. And using the temperature calculated by S6, it returns to ventilation network calculation (S5), and also calculates a ventilation quantity. These ventilation network calculation (S5) and thermal network calculation (S6) are repeated until convergence. Then, the humidity and chemical substance concentration (S7) are calculated by calculating the circuit network such as humidity, and after calculating the heat load from the obtained temperature (S8), the calculation result is output (S9).

これら一連の計算(S3〜S8)は、所定時間間隔毎に繰り返し行い、各回毎に出力(S9)を行う。所定時間間隔とは、たとえば30分ごと、1時間ごとなど、適宜設定することができる。ここで、換気回路網計算(S5)と温度計算(S6)は相互に強い影響を与えるので、安定するまで複数回繰り返し計算する。   These series of calculations (S3 to S8) are repeatedly performed at predetermined time intervals, and output (S9) is performed each time. The predetermined time interval can be set as appropriate, for example, every 30 minutes or every hour. Here, since the ventilation network calculation (S5) and the temperature calculation (S6) have a strong influence on each other, they are repeatedly calculated a plurality of times until they are stabilized.

このように各回ごとに温度状況と換気状況を繰り返し解いて収束させることにより、換気量と温度の相関をとってバランスのとれた結果を得ることができ、精度の高い有意義な予測を行うことができる。   By resolving and converging the temperature and ventilation conditions each time in this way, it is possible to obtain a balanced result by correlating ventilation volume and temperature, and to make meaningful predictions with high accuracy. it can.

繰り返し計算の初回には利用すべき前回の算出結果がないため、各室の温度は適切な初期値を与えて計算する。従って結果の欲しい時刻よりもある程度前の時刻(例えば1月半前)から計算を開始し、繰り返し計算をすることにより、実際の値に近づけることができる。   Since there is no previous calculation result to be used at the first iteration, the room temperature is calculated by giving an appropriate initial value. Accordingly, the calculation can be made closer to the actual value by starting the calculation from a certain time before the desired time (for example, one and a half months before) and repeating the calculation.

また、室内の湿度や化学物質の発生速度は温度、湿度、周辺の化学物質濃度の影響を受けるため、換気状況に加えて温度情報は重要である。従ってかかる発生速度に換気による化学物質の排除能力(移動速度及び拡散速度)を合わせて予測計算を行うことにより、精度の高い空気質の予測をすることができる。   In addition, the humidity information in the room and the generation rate of chemical substances are affected by the temperature, humidity, and concentration of chemical substances in the surroundings, so temperature information is important in addition to ventilation conditions. Therefore, it is possible to predict the air quality with high accuracy by performing prediction calculation by combining the generation rate of chemical substances by ventilation (movement speed and diffusion speed) with such generation speed.

また建物内の熱回路網計算をするにあたっては、建物各部位の断熱性能を考慮し、内外温度差、室内温度差から熱流量を求めるなど、冷暖房負荷を算出するに等しい計算を行う必要がある。従って発熱源による投入熱量の積算を行うことにより、暖冷房負荷を求めることができる。   In addition, when calculating the thermal network in a building, it is necessary to consider the heat insulation performance of each part of the building and calculate the air conditioning load, such as calculating the heat flow from the internal / external temperature difference and the indoor temperature difference. . Therefore, the heating / cooling load can be obtained by integrating the amount of input heat by the heat source.

次に、各部計算について説明する。   Next, each part calculation is demonstrated.

(放射量の計算)
放射量とは日射を含む輻射を意味し、本実施形態においては直達日射量と天空日射量並びに夜間放射量を使用している。使用する市販の気象データに水平面全天日射量しか含まれていない場合には、太陽高度を元に計算により直達日射量と天空日射量に分離して計算に使用する。また太陽高度は対象となる都市の緯度、経度及び日時より算出する。放射量は、屋外の放射と屋内の放射について考慮する。
(Calculation of radiation)
The radiation amount means radiation including solar radiation, and in this embodiment, direct solar radiation amount, sky solar radiation amount and night radiation amount are used. When the commercial weather data to be used contains only the horizontal solar radiation amount, it is divided into the direct solar radiation amount and the sky solar radiation amount by calculation based on the solar altitude. The solar altitude is calculated from the latitude, longitude and date / time of the target city. The amount of radiation is taken into account for outdoor radiation and indoor radiation.

屋外の放射について、使用する記号の意味を次表に掲げる。   The meanings of symbols used for outdoor radiation are listed in the following table.

Figure 2012092647
Figure 2012092647

気象データより各時刻の法線面直達日射量Ib[W/m2]、水平面天空日射量Id[W/m2]、夜間放射量Ix[W/m2]及び太陽高度の正弦(sin h)、余弦(cos h)、太陽方位の正弦(sin A)、余弦(cos A)を得ることができる。太陽の方向を示す単位ベクトルnsは以下のように定義できる。 Normal surface direct solar radiation I b [W / m 2 ], horizontal sky solar radiation I d [W / m 2 ], nighttime radiation I x [W / m 2 ], and sine of solar altitude from weather data (sin h), cosine (cos h), solar sine (sin A), and cosine (cos A) can be obtained. Unit vector n s that indicates the direction of the sun can be defined as follows.

Figure 2012092647
Figure 2012092647

室iと屋外jをつなぐ接続面ijの外向きの法線ベクトルをneijとおくと、接続面ijのk番目の熱的に厚い部位Hijkが受ける放射量JHjik[W/m2]は以下のように表される。 The radiation vector J Hjik [W / m 2 ] received by the kth thermally thick part H ijk of the connection surface ij is denoted by n eij as the outward normal vector of the connection surface ij connecting the room i and the outdoor j. Is expressed as follows.

Figure 2012092647
Figure 2012092647

熱的に薄い部位であるドアの場合も同様に計算できる。窓については全てが反射又は透過すると考え、ガラスで吸収する放射量は常に0と考える。   The same can be calculated for a door that is a thermally thin part. All windows are considered to be reflected or transmitted, and the amount of radiation absorbed by the glass is always considered to be zero.

屋内の放射について、使用する記号の意味を次表に掲げる。   The following table shows the meanings of the symbols used for indoor radiation.

Figure 2012092647
Figure 2012092647

屋内の放射は、(1)室内に進入する放射量(窓からの日射の侵入)を積算し、(2)電灯などの室内で発生する放射量を加え、(3)積算した放射量を室を構成する壁、床、天井などに按分することで算出する。   Indoor radiation (1) Accumulate the amount of radiation entering the room (intrusion of solar radiation from windows), (2) Add the amount of radiation generated in the room, such as an electric lamp, and (3) Calculate the accumulated amount of radiation in the room. It is calculated by apportioning the walls, floors, ceilings, etc.

(1)室内に侵入する放射量は、基本的に外壁の場合と同じ考え方をする。室iと室jをつなぐ接続面ijのk番目の窓Gijkを通して室iに侵入する日射量IGijk[W]は、次式で表される。 (1) The amount of radiation entering the room is basically the same as in the case of the outer wall. The amount of solar radiation I Gijk [W] that enters the room i through the k-th window G ijk of the connection surface ij that connects the room i and the room j is expressed by the following equation.

Figure 2012092647
Figure 2012092647

放射遮蔽係数SCR、対流遮蔽係数SCCは、物性値としてあらかじめ与える。η*は標準板ガラスの日射侵入率であり、日射角度の関数となっている。窓の面積AGijk[m2]は窓のうち日の当たっている部分(陰になっていない部分)の面積であり、その面積は窓及び庇の形状から幾何学的に算出できる。 The radiation shielding coefficient SCR and the convection shielding coefficient SCC are given in advance as physical property values. η * is the solar radiation penetration rate of the standard plate glass and is a function of the solar radiation angle. The area of the window AGijk [m 2 ] is the area of the window that is exposed to the sun (the part that is not shaded), and the area can be calculated geometrically from the shapes of the window and the eyelid.

(2)室内で発生する放射量Ii0[W]は境界条件として与えられており、その値をそのまま加える。以上から室内の全放射量Ii[W]は下記の式のようになる。 (2) The radiation amount I i0 [W] generated in the room is given as a boundary condition, and the value is added as it is. From the above, the total amount of radiation I i [W] in the room is given by the following equation.

Figure 2012092647
Figure 2012092647

(3)室内放射量の積算と按分は、たとえば上記で求まった室内に存在する放射量の合計Iiの按分にあたっては、放射量の50%が床に、残りの50%がそれ以外の面に、面積に比例して配分される。配分された放射は実際には一部が反射するが、最終的にはいずれかの壁に吸収されると考え、吸収係数にかかわらず全てが吸収されるものとする。ただし、窓に当たった放射はそのまま屋外にでていく。すなわち、次式で表すことができる。   (3) For the summation and apportionment of indoor radiation, for example, when apportioning the total amount of radiation Ii present in the room as described above, 50% of the radiation is on the floor and the remaining 50% is on the other surface. , Distributed in proportion to the area. The distributed radiation is actually partly reflected, but it is assumed that it will eventually be absorbed by any wall, and all will be absorbed regardless of the absorption coefficient. However, the radiation that hits the window goes outside as it is. That is, it can be expressed by the following formula.

Figure 2012092647
Figure 2012092647

ただし、AFi[m2]は室iを構成している床の面積であり、Ai[m2]はAFiを含む室iを構成する全表面積である。 However, A Fi [m 2 ] is the area of the floor constituting the room i, and A i [m 2 ] is the total surface area constituting the room i including A Fi .

(換気回路網計算)
換気回路網計算において、使用する記号の意味を次表に掲げる。
(Ventilation network calculation)
In the ventilation network calculation, the meaning of the symbols used is listed in the following table.

Figure 2012092647
Figure 2012092647

本実施形態においては、室を出入りする空気量を(1)換気設備による機械換気に合わせた設定値、(2)すきま風を考慮して算出した値としている。(1)強制的な換気とは、換気扇などの機械的なものや、使用者による室間若しくは建物内外間の換気をいい、流量を明示的に指定する。また、(2)すきま風は、屋外風速と、室内外、室間温度差を機動力として開口部を流れるものであり、この項目には各室の温度が影響するものである。換気回路網においては、空気の非圧縮性から換気量は瞬時に平衡に達するものと考え、擬似的に定常計算を行う。   In the present embodiment, the amount of air entering and exiting the room is (1) a set value that matches the mechanical ventilation by the ventilation facility, and (2) a value calculated in consideration of the draft. (1) Forced ventilation refers to mechanical ventilation such as a ventilation fan, or ventilation between rooms or inside and outside the building by the user, and explicitly specifies the flow rate. Further, (2) the draft air flows through the opening using the outdoor wind speed and the indoor / outdoor / room temperature difference as the motive power, and this item is affected by the temperature of each room. In the ventilation network, it is considered that the ventilation volume instantaneously reaches equilibrium due to the incompressibility of air, and pseudo steady calculation is performed.

室iと室jに存在するk番目の開口部(隙間)Mijkを流れる換気流量Vijk[m3/s]と圧力差△pijk[Pa]は次式の関係にある。 The ventilation flow rate V ijk [m 3 / s] flowing through the k-th opening (gap) M ijk existing in the chamber i and the chamber j and the pressure difference Δp ijk [Pa] have the following relationship.

Figure 2012092647
Figure 2012092647

nは通常1〜2であるが、単純な開口の場合はn=2であり、本実施形態においてもn=2としている。一方、空気密度ρA及び開口部面積AMijk[m3]は、開いている扉のような場合は大きな開口部として取り扱い、小さい隙間は相当隙間面積を代入する。 n is usually 1 to 2, but n = 2 in the case of a simple opening, and n = 2 in this embodiment. On the other hand, the air density ρ A and the opening area A Mijk [m 3 ] are handled as a large opening in the case of an open door, and the corresponding gap area is substituted for a small gap.

一方、圧力差△pijkは次の式から計算される。 On the other hand, the pressure difference Δp ijk is calculated from the following equation.

Figure 2012092647
Figure 2012092647

ただし、p* ijk、p* jikは開口部のi側及びj側の浮力を考慮した圧力[Pa]、hi、hijkは室i及び開口部Mijkの高さ[m]である。δρiは標準温度(θ0=20℃)の時と密度差[kg/m3]であり、以下の式で表される。 Here, p * ijk and p * jik are pressures [Pa] in consideration of the buoyancy on the i side and j side of the opening, and h i and h ijk are the heights [m] of the chamber i and the opening M ijk . Δρ i is the density difference [kg / m 3 ] from the standard temperature (θ 0 = 20 ° C.), and is expressed by the following equation.

Figure 2012092647
Figure 2012092647

kρ=0.00411kg/m3℃であり、室iの温度θi[℃]は後述する熱回路網から求められる。この式を7に代入すると、次式のようになる。 k ρ = 0.00411 kg / m 3 ° C, and the temperature θ i [° C] of the chamber i is obtained from a thermal circuit network described later. Substituting this expression into equation (7), the following equation.

Figure 2012092647
Figure 2012092647

圧力piは室iが文字通り屋内の「室」であるときは未知数であるが、外部空間である場合は外気温度と風速から境界条件として与える。風による外圧は次式に示す動圧として表現される。 The pressure p i is unknown when the chamber i is literally an indoor “chamber”, but is given as a boundary condition from the outside air temperature and the wind speed when it is an external space. The external pressure due to wind is expressed as a dynamic pressure expressed by the following equation.

Figure 2012092647
Figure 2012092647

これが9のpiに相当すると考え、壁の圧力は次式となる。 Believed that this is equivalent to p i of Equation 9, the pressure of the wall becomes the following equation.

Figure 2012092647
Figure 2012092647

ここでhiは屋外の中心高さであり、どこに定義しても相対的に全ての圧力が平衡移動するだけなので、hi=0としている。 Here, h i is the height of the center of the outdoors, and all pressures move relatively in equilibrium no matter where they are defined, so h i = 0.

一方α*は風圧係数(変数であり、形状により入力する)であり、壁と風向きとの位置関係の関数となる。一例として文献によれば風上で+0.8、風下で-0.4なので、それを滑らかに繋ぐと、単位風向ベクトルをv/|v|、壁の単位法線ベクトルをneとして、次式で表すことができる。 On the other hand, α * is a wind pressure coefficient (which is a variable and is input according to the shape) and is a function of the positional relationship between the wall and the wind direction. As an example, according to the literature, it is +0.8 on the windward side and -0.4 on the leeward side. When smoothly connected, the unit wind direction vector is v / | v | and the unit normal vector of the wall is ne. Can do.

Figure 2012092647
Figure 2012092647

また、各室についての質量保存の式、すなわち次式が室の数だけできる。   Moreover, the mass conservation formula for each chamber, that is, the following formula, can be created for the number of chambers.

Figure 2012092647
Figure 2012092647

上記説明した6が各開口部について1つづつ、13が室の数だけ作成され、未知数は各開口部の流量Vijkが開口部の数だけ、圧力piが室数だけ存在する。従って、式の数と未知数の数が一致するので、この連立方程式は原理的に解くことができる。なお、表3においては、室jから室iに流れる流量との観点から各項が定義付けられているが、図1(b)等に示す如く、当該室iや室jは当然に屋外も含み、そうすると、屋外に接する部屋においては、屋外と室との関係についても上記処理により演算処理されることとなる。 The number 6 described above is created one for each opening, the number 13 is created for the number of chambers, and the unknowns are the flow rate V ijk for each opening for the number of openings and the pressure p i for the number of chambers. Therefore, since the number of equations coincides with the number of unknowns, this simultaneous equation can be solved in principle. In Table 3, each term is defined from the viewpoint of the flow rate flowing from the room j to the room i. However, as shown in FIG. In that case, in a room that is in contact with the outdoors, the relationship between the outdoors and the room is also processed by the above processing.

6から明らかなように、換気回路網は二次の連立方程式となり、熱回路網や湿気回路網のように線形でない。そこで本回路網は繰り返し解法による収束計算を行っている。最終的に解を求められるならどのような方法を用いてもよいが、本実施形態においてはNewton法を基礎にMarquardt法の考え方を応用している。 As is apparent from Equation 6, the ventilation network is a quadratic simultaneous equation and is not linear like the thermal network or the humidity network. Therefore, this network performs convergence calculation by iterative solution. Any method can be used as long as a solution can be finally obtained. In this embodiment, the concept of the Marquardt method is applied based on the Newton method.

(熱回路網計算)
熱回路網計算において、使用する記号の意味を次表に掲げる。
(Thermal network calculation)
The meanings of the symbols used in the thermal network calculation are listed in the following table.

Figure 2012092647
Figure 2012092647

本実施形態においては、室iに伝わってくる熱量として(1)熱的に厚い部位から伝わる熱量、(2)熱的に薄い部位から伝わる熱量、(3)土間床を通って地中から伝わる熱量、(4)換気によって移動する熱量、(5)発熱源から発せられる熱量を考慮して算出する。(1)熱的に厚い部位とは壁や床天井などをいい、熱容量を考慮して計算するものであり、熱緩衝となって熱伝達の経時変化を追うことができる。(2)熱的に薄い部位とはドアやガラスなどをいい、熱容量を考慮せずに計算を行うものである。(3)土間床とは地盤に接している床から伝達される熱である。(4)換気によって移動する熱量には、上記換気回路網によって算出される空気の流量を用いて計算する。(5)発熱源とは暖冷房機器や人などを意味し、熱量を明示的に指定する。   In the present embodiment, the amount of heat transmitted to the chamber i is (1) the amount of heat transmitted from a thermally thick portion, (2) the amount of heat transmitted from a thermally thin portion, and (3) transmitted from the ground through the dirt floor. It is calculated taking into account the amount of heat, (4) the amount of heat transferred by ventilation, and (5) the amount of heat generated from the heat source. (1) A thermally thick portion refers to a wall or floor / ceiling, and is calculated in consideration of the heat capacity, and can serve as a thermal buffer to follow changes over time in heat transfer. (2) A thermally thin part means a door, glass, etc., and performs calculation without considering heat capacity. (3) Soil floor is heat transmitted from the floor in contact with the ground. (4) The amount of heat transferred by ventilation is calculated using the air flow rate calculated by the ventilation network. (5) The heat source means a heating / cooling device or a person, and explicitly specifies the amount of heat.

(1)熱的に厚い部位の算出にあたっては、熱的に厚い部位の有効熱容量について、部位の熱伝導抵抗の中心を熱的な境界として、部位を分割する。室iとjを結ぶk番目の厚い部位Hijkのi側及びj側の熱容量CHijk[J/m2℃]、C Hjik [J/m2℃]は以下のように定義される。 (1) In the calculation of the thermally thick region, with the effective heat capacity of the thermally thick region, the center of the thermal resistivity of the site as a thermal boundary, divides the site. The heat capacities C Hijk [J / m2 ° C.] and C Hjik [J / m 2 ° C.] on the i side and the j side of the kth thick portion H ijk connecting the chambers i and j are defined as follows.

Figure 2012092647
Figure 2012092647

ここで、pは部位Hijkの幅d以内の値であって、 Here, p is a value within the width d of the part H ijk ,

Figure 2012092647
Figure 2012092647

となるような値である。分割した部位の両側について、それぞれ表面温度THijk、THjikをその代表温度とする(図3参照)。 Is such a value. The surface temperatures T Hijk and T Hjik are used as representative temperatures for both sides of the divided parts (see FIG. 3).

Hijkのi側及びj側の熱収支式は次式となる。   The heat balance equation of i side and j side of Hijk is as follows.

Figure 2012092647
Figure 2012092647

ここでJHijk等は放射により部位Hijkのi側表面に得られる熱量であり、上記放射量の計算によって得られた値を用いる。この式を後退差分 Here, J Hijk or the like is the amount of heat obtained on the i-side surface of the site Hijk by radiation, and the value obtained by the calculation of the amount of radiation is used. This formula is used as a backward difference

Figure 2012092647
Figure 2012092647

を用いて離散化すると、次式となる。   When discretized using, the following equation is obtained.

Figure 2012092647
Figure 2012092647

18が熱的に厚い部位の数だけできるため、室iに流入する熱量qHijkは次式となる。 Since the number 18 can be made by the number of thermally thick portions, the amount of heat qHijk flowing into the chamber i is expressed by the following equation.

Figure 2012092647
Figure 2012092647

(2)熱的に薄い部位の算出にあたっては、室iとjをつなぐ接合面ijのk番目の熱的に薄い部位(ドア、ガラスなど)Lijkについての壁の左半分と右半分の収支式は、熱的に薄い部位の熱容量はゼロと見なせることから、次式のようになる。 (2) When calculating the thermally thin part, the balance of the left half and right half of the wall for the kth thermally thin part (door, glass, etc.) L ijk of the joint surface ij connecting the chambers i and j. Since the heat capacity of the thermally thin portion can be regarded as zero, the equation becomes as follows.

Figure 2012092647
Figure 2012092647

従って室iに流入する熱量qLijは、次式となる。 Accordingly, the amount of heat q Lij flowing into the chamber i is expressed by the following equation.

Figure 2012092647
Figure 2012092647

(3)土間床を伝って地中から伝わる熱量も壁の場合と同じく熱回路網で計算するが、壁の場合は両面空気層であるのに対し、土間床は地盤に接している点が異なっている。土間床は図4に示す如きイメージでモデル化している。壁の内側半分だけを取り出して直接地盤にくっつけたようなモデルを考え、さらに外周部分と中央部分に分けて計算している。   (3) The amount of heat transmitted from the ground through the dirt floor is also calculated by the thermal network as in the case of the wall. In the case of the wall, the double-layer air layer is used, whereas the dirt floor is in contact with the ground. Is different. The clay floor is modeled with the image shown in FIG. Considering a model in which only the inner half of the wall is taken out and attached directly to the ground, the calculation is divided into the outer peripheral part and the central part.

内側半分だけなので、壁の場合の屋外側の境界層熱抵抗γHjikや、CHjikT Hjik 相当するものがない。計算手法上は、γHjik=CHjik=0、T Hjik jとなった壁と考えると壁と同じ計算となる。θjは土間床の場合は地盤温度(θG)で表され、γ Hijk 、C Hijk THijkはそれぞれγGik、CGik、TGik等に置き換えられる。 Since there is only the inner half, there is nothing equivalent to the boundary boundary thermal resistance γ Hjik , C Hjik , or T Hjik on the outdoor side in the case of walls. In terms of the calculation method, the calculation is the same as that of the wall, assuming that γ Hjik = C Hjik = 0 and T Hjik = θ j . θ j is represented by the ground temperature (θ G ) in the case of a soil floor, and γ Hijk , C Hijk , T Hijk are replaced with γ Gik , C Gik , T Gik, etc., respectively.

従って土間床Gikの熱収支式は、以下のようになる。 Therefore, the heat balance equation of the dirt floor Gik is as follows.

Figure 2012092647
Figure 2012092647

従って室iに流入する熱量q Gik Therefore, the amount of heat q Gik flowing into chamber i is

Figure 2012092647
Figure 2012092647

以上2式に使われた変数のうち、γGik、AGikは条件として入力され、JGikは上記(1)と同様にして算出される。γGikには土間床の境膜抵抗の他に内装材の熱抵抗も含まれているため、TGikの定義点は内装材の裏側になる。結果として床への放射は現実と異なり内装材の裏側に当てられていることになる。 Of the variables used in the above two equations, γ Gik and A Gik are input as conditions, and J Gik is calculated in the same manner as (1) above. Since γ Gik includes the thermal resistance of the interior material in addition to the film resistance of the floor, the definition point of T Gik is on the back side of the interior material. As a result, the radiation to the floor is applied to the back side of the interior material unlike the actual situation.

(4)換気によって移動する熱量は、開口部Mijkを室jからiに通過する流量をVijk[m3/s]とすると、次式となる。 (4) The amount of heat transferred by ventilation is given by the following equation, where V ijk [m 3 / s] is the flow rate passing through the opening M ijk from the chamber j to i.

Figure 2012092647
Figure 2012092647

室間の流量Vijkは換気回路網により計算されるか、又は境界条件として与えられる。 The flow rate V ijk between the rooms is calculated by the ventilation network or given as a boundary condition.

(5)発熱源とは暖冷房機器や人などの発する定量的な熱源をいい、発熱量Qを条件として入力する。   (5) The heat source is a quantitative heat source generated by a heating / cooling device or a person, and the heat generation amount Q is input as a condition.

以上で定まった熱流速の式を用いた室iにおける熱収支式は、次式のようになる。   The heat balance equation in the chamber i using the equation of the heat flow rate determined above is as follows.

Figure 2012092647
Figure 2012092647

この25に19、21、23、24を代入することにより、Tとθに関する式ができる。更に18、20を用いることによりθ={θ12,…θN}を変数とする以下のような連立一次方程式ができる。 By substituting the number 19,21,23,24 to the number 25, it is the equation for T and theta. Further theta = By using several 18,20 {θ 1, θ 2, ... θ N} the can simultaneous linear equations such as: whose variable.

Figure 2012092647
Figure 2012092647

ここでAはNxNの行列、b、θは長さNのベクトルである。これを解けば時刻ステップnにおける室温θiが求められる。本実施形態ではこの解法として直接法(ガウス消去法)を用いている。 Here, A is an N × N matrix, and b and θ are vectors of length N. Solving this, the room temperature θ i at time step n can be obtained. In this embodiment, the direct method (Gaussian elimination method) is used as this solution.

(湿度等回路網計算)
湿度等回路網計算において、使用する記号の意味を次表に掲げる。
(Humidity network calculation)
The following table shows the meanings of the symbols used in calculating the network such as humidity.

Figure 2012092647
Figure 2012092647

湿度等回路網は湿度や化学物質の濃度を算出するものであるが、本質的には熱回路網と同じである。ただし壁を通しての移動がないため、その分簡単になる。水蒸気濃度についての室iの保存式は次式のようになる。   A circuit network such as humidity is used to calculate humidity and chemical substance concentration, but is essentially the same as a thermal circuit network. However, since there is no movement through the wall, it is much easier. The conservation formula of the chamber i with respect to the water vapor concentration is as follows.

Figure 2012092647
Figure 2012092647

27において、Ci/R=1、ρA=1、Li/R=0とすると化学物質の保存式となる。 In Equation 27, if Ci / R = 1, ρA = 1, and Li / R = 0, the chemical substance conservation formula is obtained.

ここでVijkは室jからiへのk番目の開口部からの換気量[m3/s]である。上記の式は部屋の数Nだけ作ることができ、N元の連立一次方程式を解いて、室湿度または潜熱負荷を計算することができる。 Here, V ijk is a ventilation amount [m 3 / s] from the k-th opening from room j to i. The above formula can be created for the number N of rooms, and the room humidity or latent heat load can be calculated by solving simultaneous linear equations of N elements.

(熱負荷計算)
熱負荷を求めるときは、上記熱回路網計算の26において室温θに維持したい温度を与え、同式の残差(H=Aθ-b)[W]として求められる。すなわち設定温度、室内温度、熱容量が設定されることにより、負荷量=室熱容量[J/m3]*(設定温度−室内温度)で求められる。なお、冷房の場合は−1倍となる。
(Heat load calculation)
When determining the thermal load, gives a temperature to be maintained at room temperature θ in the number 26 of the heat network calculation is determined as the residual of the equation (H = Aθ-b) [ W]. That is, by setting the set temperature, the room temperature, and the heat capacity, the load amount = the room heat capacity [J / m 3 ] * (set temperature−room temperature) is obtained. In the case of cooling, it is -1.

[第二実施形態]次に、本発明に係る建物の換気量及び温度予測システムの第二実施形態について、図5を用いて説明する。図5は本実施形態に係る熱及び換気回路網計算のフローチャートであって、上記第一実施形態と説明の重複する部分については、同一の符号を付して説明を省略する。   [Second Embodiment] Next, a second embodiment of a building ventilation and temperature prediction system according to the present invention will be described with reference to FIG. FIG. 5 is a flowchart of heat and ventilation network calculation according to the present embodiment, and the same parts as those in the first embodiment are denoted by the same reference numerals and the description thereof is omitted.

上記第一実施形態においては所定時間間隔毎の各回の計算においては、換気回路網計算(S5)及び熱回路網計算(S6)は各回において収束するまで複数回繰り返すよう説明した。   In the first embodiment described above, in each calculation at every predetermined time interval, the ventilation network calculation (S5) and the thermal network calculation (S6) are described to be repeated a plurality of times until they converge at each time.

しかし本実施形態においては、図5に示すように、換気回路網計算(S5)と熱回路網計算(S6)を各回において1回のみ行い、結果を出力し、次の時刻の計算に進む。   However, in this embodiment, as shown in FIG. 5, the ventilation network calculation (S5) and the thermal network calculation (S6) are performed only once each time, the result is output, and the process proceeds to the next time calculation.

上記の如く構成することにより、第一実施形態に示したよりも若干の精度低下にはなるものの、より短い期間で計算結果を得ることができる。換気回路網は時間とともに刻々と変化する風向風速の影響によって時間間隔毎の変化が激しいため、前の時刻の換気量を用いて温度計算を行うことは大きな誤差を招く可能性が高いのに対し、温度の変化は時間に対して比較的穏やかであり、前の時刻の温度を用いて換気量を計算してもそれほど大きな誤差はおきないと考えられ、実際の計算でも第一実施形態と第二実施形態では殆ど差はなかった。   By configuring as described above, the calculation result can be obtained in a shorter period, although the accuracy is slightly lower than that shown in the first embodiment. Ventilation network changes drastically at each time interval due to the influence of wind direction and wind speed that changes with time, so calculating temperature using the ventilation volume at the previous time is likely to cause a large error. The change in temperature is relatively gentle with respect to time, and even if the ventilation volume is calculated using the temperature at the previous time, it is considered that there will be no significant error. There was little difference between the two embodiments.

[計算例]図6は図で示されているようなレイアウトにおける建物内の温度分布を計算した例である。屋外の条件は東京地方の寒い時期の気候であり、前日の22時までエアコンをつけて各部屋の温度を20℃にし、その後エアコンを切って成り行きに任せ、翌朝6時における温度分布を示している。なおドアや窓などは閉めた状態とし、換気扇などの強制換気は考慮する(以下において、これを標準状態とする)。   [Calculation Example] FIG. 6 shows an example of calculating the temperature distribution in the building in the layout as shown in the figure. The outdoor conditions are the cold weather in the Tokyo region. Turn on the air conditioner until 22:00 on the previous day and set the temperature of each room to 20 ° C. Yes. Note that doors and windows are closed, and forced ventilation such as a ventilation fan is considered (hereinafter, this is the standard state).

図6に示すように、1階よりも2階の方が冷えていることがシミュレーションによって示され(または予測され)ている。これは単なる伝熱による熱移動のほかに、暖かい空気が浮力により上昇するために、1階は外気に比べて負圧、2階は正圧になる傾向をもっていることによる影響が考慮されている。また便所など機械的な排気を行っている部屋は負圧になるためそれを補うために隣の廊下から空気の流入があり、さらにこの影響で、廊下も負圧になりそれを補うために居間などの比較的高温の部屋から気流が入ってきていて、このため他の部屋と比べて室温が高くなっていることがシミュレーションによって示されている。   As shown in FIG. 6, it is shown (or predicted) by simulation that the second floor is colder than the first floor. In addition to heat transfer due to simple heat transfer, warm air rises due to buoyancy, so the influence of the first floor tends to be negative and the second is positive compared to outside air is taken into account. . In addition, rooms that are mechanically exhausted, such as toilets, have negative pressure, so there is an inflow of air from the adjacent corridor to compensate for this, and due to this effect, the corridor also becomes negative pressure and the living room The simulation shows that airflow has entered from a relatively hot room such as that the room temperature is higher than in other rooms.

このような気流やその連鎖による熱の移動の影響は熱の計算のために伝熱による熱回路網だけでなく換気回路網の計算を行ったことにより初めて計算できることであり、一方換気回路網の、特に浮力の影響の計算のためには熱回路網の計算が必要である。すなわち両方の回路網を交互に解いた事により初めて可能となるシミュレーション結果である。   The effect of heat transfer due to such airflow and its chain can only be calculated by calculating the ventilation network as well as the heat network by heat transfer for the calculation of heat. Especially for the calculation of the influence of buoyancy, the calculation of the thermal network is necessary. In other words, it is a simulation result that is possible only by solving both networks alternately.

図7は、上記レイアウトにおける建物内の換気回数の分布を表している。ここで換気回数とは、所定時間の間に流通する空気の量を部屋等の容積で割った値であり、室内の空気が入れ替わった回数を示している。図7はドアなどを閉めた状態(標準状態)において、自然に発生する風向、風量の風に対応した内部換気程度の分布を示している。   FIG. 7 shows the distribution of the ventilation frequency in the building in the above layout. Here, the ventilation frequency is a value obtained by dividing the amount of air circulated during a predetermined time by the volume of the room or the like, and indicates the number of times the room air has been replaced. FIG. 7 shows the distribution of the degree of internal ventilation corresponding to the naturally occurring wind direction and air volume when the door is closed (standard state).

図8は、上記標準状態に対して、二階階段ホールに接する窓を少し開けた場合であって、いずれも適当な一日分の平均値を表示している。図に示すように、閉め切った場合に比して、建物のすきま面積が拡大した場合を想定し、二階階段ホールに接する窓を開けた方が、当該ホールに隣接する部分の換気が増えていることがわかる。   FIG. 8 shows a case where a window that is in contact with the second-level staircase is slightly opened with respect to the above-described standard state. As shown in the figure, it is assumed that the clearance area of the building has expanded compared to when it is closed, and opening the window in contact with the second floor staircase increases ventilation in the area adjacent to the hall. I understand that.

図9は、二階階段ホールに接する窓を少し開けた場合における建物内の温度分布を計算した例である。このように、隙間を変化させた場合の換気性状が異なることにより、これを考慮して計算を行った結果、建物内の温度分布もその影響を受けて変わることがわかる。すなわち、図9に示すように、二階階段ホール及びこれに隣接する廊下、及びこれに連結している一階の階段部分は、換気量が増加すると共に、温度が下がっていることがわかる。   FIG. 9 is an example in which the temperature distribution in the building is calculated when a few windows in contact with the second floor staircase are opened. As described above, since the ventilation characteristics when the gap is changed are different, it is understood that the temperature distribution in the building also changes due to the influence as a result of calculation taking this into consideration. That is, as shown in FIG. 9, it can be seen that the second floor stair hall, the corridor adjacent to the second floor stair hall, and the first floor stair portion connected to the second floor stair hall increase the ventilation amount and the temperature decreases.

C …熱容量
H …厚い部位
J …表面に得られる熱量
R …有効熱容量
T …表面温度
d …幅
γ …境界層熱抵抗
θ …温度
C ... heat capacity
H ... thick part
J ... The amount of heat obtained on the surface
R: Effective heat capacity
T ... surface temperature
d… Width γ… Boundary layer thermal resistance θ… Temperature

Claims (2)

任意の時間の、建物の各室の温度及び室間若しくは建物内外間の換気量を予測するシステムであって、
各室の形状及び壁、床、開口部など室相互の接続関係を含む邸別データと、風速及び外気温度を含む気象データとを入力する入力手段と、
所定時間間隔毎に各室毎の温度および換気量を計算する計算手段と、
計算結果を出力する出力手段とを有し、
前記計算手段は
温度による空気の密度の違いによる浮力を考慮して、各室における空気の換気量を求め、かつ各室における空気の流出入の総和を0に維持することにより各開口部毎の換気量を計算し、
換気による空気に伴う熱の移動を考慮して各室相互の熱の入出量により各室毎の温度を計算し、
前記換気量の計算にて前記温度による空気の密度の違いによる浮力を考慮するにつき、当該温度は、前記各室毎の温度の計算から得られる温度を用いることとし、
所定時間間隔毎に、毎回、上記換気量と温度とを交互に複数回計算して収束させることにより、換気量及び温度を計算する
ことを特徴とする建物の換気量及び温度予測システム。
A system for predicting the temperature of each room of a building and the ventilation between rooms or inside and outside the building at any time,
Input means for inputting house-by-house data including the connection between the rooms, such as the shape of each room and walls, floors, openings, etc., and weather data including wind speed and outside air temperature,
Calculation means for calculating the temperature and ventilation volume for each room at predetermined time intervals;
Output means for outputting the calculation result,
Said calculating means,
Taking into account the buoyancy due to the difference in air density due to temperature, calculate the air ventilation volume in each room, and maintain the sum of air inflow and outflow in each room to calculate the ventilation volume for each opening. And
Considering the heat transfer with air due to ventilation, calculate the temperature of each room by the amount of heat input and output between each room,
In considering the buoyancy due to the difference in air density due to the temperature in the calculation of the ventilation volume, the temperature is a temperature obtained from the calculation of the temperature for each room,
A ventilation volume and temperature prediction system for a building which calculates ventilation volume and temperature by calculating and converging the ventilation volume and temperature alternately a plurality of times every predetermined time interval.
任意の時間の、建物の各室の温度及び室間若しくは建物内外間の換気量を予測するシステムであって、
各室の形状及び開口部を含む邸別データと、風速及び外気温度を含む気象データとを入力する入力手段と、
所定時間間隔毎に各室毎の温度および換気量を計算する計算手段と、
計算結果を出力する出力手段とを有し、
前記計算手段は
温度による空気の密度の違いによる浮力を考慮して各室における空気の換気量を求め、かつ各室における空気の流出入の総和を0に維持することにより各開口部毎の換気量を計算し、
換気による空気に伴う熱の移動を考慮して各室相互の熱の入出量により各室毎の温度を計算し、
前記換気量の計算にて前記温度による空気の密度の違いによる浮力を考慮するにつき、当該温度は、前記各室毎の温度の計算から得られる温度を用いることとし、
所定時間間隔毎に、前時間で算出した温度を用いて換気量を先に計算し、その換気量の計算結果を用いて温度を計算する
ことを特徴とする建物の換気量及び温度予測システム。
A system for predicting the temperature of each room of a building and the ventilation between rooms or inside and outside the building at any time,
Input means for inputting house-specific data including the shape and opening of each room, and weather data including wind speed and outside air temperature;
Calculation means for calculating the temperature and ventilation volume for each room at predetermined time intervals;
Output means for outputting the calculation result,
Said calculating means,
Taking into consideration the buoyancy due to the difference in density of air with temperature, determine the amount of ventilation air in each chamber, and calculating the ventilation of the respective openings by maintaining the sum of the inflow and outflow of air in each chamber to 0 And
Considering the heat transfer with air due to ventilation, calculate the temperature of each room by the amount of heat input and output between each room,
In considering the buoyancy due to the difference in air density due to the temperature in the calculation of the ventilation volume, the temperature is a temperature obtained from the calculation of the temperature for each room,
A building ventilation volume and temperature prediction system characterized in that, at a predetermined time interval, the ventilation volume is calculated first using the temperature calculated in the previous time, and the temperature is calculated using the calculation result of the ventilation volume.
JP2012008166A 2012-01-18 2012-01-18 Building ventilation and temperature prediction system Expired - Lifetime JP5383832B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2012008166A JP5383832B2 (en) 2012-01-18 2012-01-18 Building ventilation and temperature prediction system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2012008166A JP5383832B2 (en) 2012-01-18 2012-01-18 Building ventilation and temperature prediction system

Related Parent Applications (1)

Application Number Title Priority Date Filing Date
JP2010074451A Division JP5008739B2 (en) 2010-03-29 2010-03-29 Building ventilation and temperature prediction system

Publications (2)

Publication Number Publication Date
JP2012092647A true JP2012092647A (en) 2012-05-17
JP5383832B2 JP5383832B2 (en) 2014-01-08

Family

ID=46386276

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2012008166A Expired - Lifetime JP5383832B2 (en) 2012-01-18 2012-01-18 Building ventilation and temperature prediction system

Country Status (1)

Country Link
JP (1) JP5383832B2 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015118738A1 (en) * 2014-02-07 2015-08-13 株式会社東芝 Heat retention capacity calculation device, energy management system, calculation method, and program

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04142A (en) * 1990-04-16 1992-01-06 Shimizu Corp Hot heat environment control system
JPH108581A (en) * 1996-06-20 1998-01-13 Kajima Corp Ventilating method and device
JPH11241849A (en) * 1998-02-24 1999-09-07 Sekisui Chem Co Ltd Ventilating system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH04142A (en) * 1990-04-16 1992-01-06 Shimizu Corp Hot heat environment control system
JPH108581A (en) * 1996-06-20 1998-01-13 Kajima Corp Ventilating method and device
JPH11241849A (en) * 1998-02-24 1999-09-07 Sekisui Chem Co Ltd Ventilating system

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2015118738A1 (en) * 2014-02-07 2015-08-13 株式会社東芝 Heat retention capacity calculation device, energy management system, calculation method, and program
JP2015148416A (en) * 2014-02-07 2015-08-20 株式会社東芝 Heat insulation performance calculation device, energy management system, calculation method, and program

Also Published As

Publication number Publication date
JP5383832B2 (en) 2014-01-08

Similar Documents

Publication Publication Date Title
JP4570280B2 (en) Building ventilation and temperature prediction system
JP5008738B2 (en) Building ventilation and temperature prediction system
Stavrakakis et al. Optimization of window-openings design for thermal comfort in naturally ventilated buildings
Soares et al. Multi-dimensional optimization of the incorporation of PCM-drywalls in lightweight steel-framed residential buildings in different climates
Shen et al. Numerical study on thermal behavior of classical or composite Trombe solar walls
Seferis et al. Investigation of the performance of a ventilated wall
Zingre et al. Dynamic thermal performance of inclined double-skin roof: Modeling and experimental investigation
Stavrakakis et al. Development of a computational tool to quantify architectural-design effects on thermal comfort in naturally ventilated rural houses
Munaretto et al. Empirical validation of different internal superficial heat transfer models on a full-scale passive house
Vallianos et al. Hybrid ventilation in an institutional building: Modeling and predictive control
Koinakis Combined thermal and natural ventilation modeling for long-term energy assessment: validation with experimental measurements
Dardir et al. Development of PCM-to-air heat exchanger for integration in building envelope–modeling and validation
Arendt et al. Influence of input data on airflow network accuracy in residential buildings with natural wind-and stack-driven ventilation
Zhai et al. Review of natural ventilation models
Thapa et al. Energy conservation in buildings–a review
JP5008739B2 (en) Building ventilation and temperature prediction system
Dimitriou et al. Data-driven simple thermal models: The importance of the parameter estimates
JP5008740B2 (en) Building ventilation and temperature prediction system
Tol et al. Development of a white-box dynamic building thermal model integrated with a heating system
Jayasree et al. Assessment of air change effectiveness and thermal comfort in a naturally ventilated kitchen with insect-proof screen using CFD
Yang The curious case of urban heat island: A systems analysis
JP5383832B2 (en) Building ventilation and temperature prediction system
JP5398851B2 (en) Building ventilation and temperature prediction system
Edtmayer et al. Modelling virtual sensors for real-time indoor comfort control
Leclere et al. A low order envelope model for optimised predictive control of indoor temperature: development methodology and calibration with a numerical model

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20130326

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20130523

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20131001

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20131001

R150 Certificate of patent or registration of utility model

Ref document number: 5383832

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150

Free format text: JAPANESE INTERMEDIATE CODE: R150

S531 Written request for registration of change of domicile

Free format text: JAPANESE INTERMEDIATE CODE: R313531

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350

EXPY Cancellation because of completion of term