"開口で連通する2つの対象空間のうち一方にCFDモデルを適用して解析し、他方に二層ゾーンモデルを適用して解析する連成解析方法であって、前記他方の対象空間の圧力および温度が前記一方の対象空間の各セルに圧力、温度、および流速に与える影響を算定する手順と、CFDモデル解析により前記一方の対象空間の各セルの圧力、温度、および流速を算定する手順と、前記一方の対象空間の上部層と下部層との境界である層境界高さを前記他方の対象空間の前回の解析結果に基づいて算定する手順と、当該層境界高さに基づいて前記2つの対象空間の境界面における質量流量およびエネルギー量を分配する手順と、当該分配された質量流量およびエネルギー量に基づいて、二層ゾーンモデル解析により前記他方の対象空間の煙層と空気層との境界である層境界高さ、圧力、および温度を算定する手順と、を順番に繰り返すことを特徴とする連成解析方法。前記層境界高さに基づいて前記2つの対象空間の境界面における質量流量およびエネルギー量を分配する手順では、噴流プルームによる周辺空気の巻き込み量を考慮することを特徴とする請求項1に記載の連成解析方法。開口で連通する2つの対象空間のうち一方にCFDモデルを適用して解析し、他方に二層ゾーンモデルを適用して解析する連成解析システムであって、前記他方の対象空間の圧力および温度が前記一方の対象空間の各セルの圧力、温度、および流速に与える影響を算定し、CFDモデル解析により前記一方の対象空間の各セルの圧力、温度、および流速を算定するCFDモデル解析手段と、前記一方の対象空間の上部層と下部層との境界である層境界高さを前記他方の対象空間の前回の解析結果に基づいて算定し、当該層境界高さに基づいて前記2つの対象空間の境界面における質量流量およびエネルギー量を分配する連成処理手段と、当該分配された質量流量およびエネルギー量に基づいて、二層ゾーンモデル解析により前記他方の対象空間の煙層と空気層との境界である層境界高さ、圧力、および温度を算定する二層ゾーンモデル解析手段と、を備えることを特徴とする連成解析システム。前記連成処理手段は、噴流プルームによる周辺空気の巻き込み量も考慮して、前記2つの対象空間の境界面における質量流量およびエネルギー量を分配することを特徴とする請求項3に記載の連成解析システム。しかしながら、以上の連成解析では、各モデル間のパラメータの受け渡し方法が不明であるため、連成解析を実現できない、という問題があった。本発明は、連成解析を確実に実現できる連成解析方法および連成解析システムを提供することを目的とする。請求項1に記載の連成解析方法は、開口で連通する2つの対象空間のうち一方にCFDモデルを適用して解析し、他方に二層ゾーンモデルを適用して解析する連成解析方法であって、前記他方の対象空間の圧力および温度が前記一方の対象空間の各セルに圧力、温度、および流速に与える影響を算定する手順と、CFDモデル解析により前記一方の対象空間の各セルの圧力、温度、および流速を算定する手順と、前記一方の対象空間の上部層と下部層との境界である層境界高さを前記他方の対象空間の前回の解析結果に基づいて算定する手順と、当該層境界高さに基づいて前記2つの対象空間の境界面における質量流量およびエネルギー量を分配する手順と、当該分配された質量流量およびエネルギー量に基づいて、二層ゾーンモデル解析により前記他方の対象空間の煙層と空気層との境界である層境界高さ、圧力、および温度を算定する手順と、を順番に繰り返すことを特徴とする。請求項2に記載の連成解析方法は、前記層境界高さに基づいて前記2つの対象空間の境界面における質量流量およびエネルギー量を分配する手順では、噴流プルームによる周辺空気の巻き込み量を考慮することを特徴とする。請求項3に記載の連成解析システムは、開口で連通する2つの対象空間のうち一方にCFDモデルを適用して解析し、他方に二層ゾーンモデルを適用して解析する連成解析システムであって、前記他方の対象空間の圧力および温度が前記一方の対象空間の各セルの圧力、温度、および流速に与える影響を算定し、CFDモデル解析により前記一方の対象空間の各セルの圧力、温度、および流速を算定するCFDモデル解析手段と、前記一方の対象空間の上部層と下部層との境界である層境界高さを前記他方の対象空間の前回の解析結果に基づいて算定し、当該層境界高さに基づいて前記2つの対象空間の境界面における質量流量およびエネルギー量を分配する連成処理手段と、当該分配された質量流量およびエネルギー量に基づいて、二層ゾーンモデル解析により前記他方の対象空間の煙層と空気層との境界である層境界高さ、圧力、および温度を算定する二層ゾーンモデル解析手段と、を備えることを特徴とする。請求項4に記載の連成解析システムは、前記連成処理手段は、噴流プルームによる周辺空気の巻き込み量も考慮して、前記2つの対象空間の境界面における質量流量およびエネルギー量を分配することを特徴とする。この発明によれば、一方の対象空間にCFDモデルを適用して解析することで、詳細な性状を把握するとともに、他方の対象空間に二層ゾーンモデルを適用して解析することで、大まかな性状を把握する。このように解析方法を使い分けることにより、計算時間を短縮しつつ十分に詳細な解析結果を得ることができる。さらに、各解析から得られる物理量を双方向で連成させることで連成解析を実現した。具体的には、二層ゾーンモデル解析から得られる圧力および温度をCFDモデル解析に取り入れるとともに、CFDモデル解析から得られる各セルの圧力、温度、および流速を、質量流量およびエネルギー量の分配を介して二層ゾーンモデル解析に取り入れた。これにより、連成解析を確実に実現できる。また、2つの対象空間の境界面における質量流量およびエネルギー量を分配する際に、噴流プルームによる周辺空気の巻き込み量も考慮したので、連成解析の精度をより向上できる。本発明によれば、一方の対象空間にCFDモデルを適用して解析することで、詳細な性状を把握するとともに、他方の対象空間に二層ゾーンモデルを適用して解析することで、大まかな性状を把握する。このように解析方法を使い分けることにより、計算時間を短縮しつつ十分に詳細な解析結果を得ることができる。さらに、各解析から得られる物理量を双方向で連成させることで連成解析を実現した。具体的には、二層ゾーンモデル解析から得られる圧力および温度をCFDモデル解析に取り入れるとともに、CFDモデル解析から得られる各セルの圧力、温度、および流速を、質量流量およびエネルギー量の分配を介して二層ゾーンモデル解析に取り入れた。これにより、連成解析を確実に実現できる。以下、本発明の一実施形態について、図面を参照しながら説明する。図1は、本発明の一実施形態に係る連成解析システム1の構成を示すブロック図である。連成解析システム1は、CFDモデル解析手段10と、連成処理手段20と、二層ゾーンモデル解析手段30と、を備える。CFDモデル解析手段10は、対象空間にCFDモデルを適用して解析する手段であり、対象空間を微細なセルに分割して、セル毎に圧力、温度、流速を計算する。二層ゾーンモデル解析手段30は、対象空間に二層ゾーンモデルを適用して解析する手段であり、対象空間を煙層と空気層とに分割して、煙層および空気層の圧力や温度を計算するものである。連成処理手段20は、CFDモデル解析手段10で計算されたデータを加工して、二層ゾーンモデル解析手段30に受け渡すデータを生成するものである。以下、連成解析システム1の動作について説明する。この連成解析システム1は、ここでは、2つの室Room1、Room2を対象空間として解析する。図2は、対象空間である2つの室Room1、Room2をモデル化した概念図である。Room1とRoom2とは隣接しており、扉などの開口を介して連通している。この状態で、Room1で火災が発生し、この火災による煙がRoom2に流れ込むとともに、Room2からRoom1に新鮮空気が流れ込む現象を解析する。Room1についてはCFDモデルを適用し、CFDモデル解析手段10により、Room1を微細なセルに分割して、セル毎に圧力、温度、流速を算定する。Room2については二層ゾーンモデルを適用し、二層ゾーンモデル解析手段30により、上部の煙層と下部の空気層とに2分割して、層毎に圧力および温度を算定する。この連成解析システム1では、CFDモデル解析手段10、連成処理手段20、二層ゾーンモデル解析手段30の順に処理を実行するサイクルを繰り返す。以下、連成解析システム1の1サイクルの処理について、図3のフローチャートを参照しながら説明する。ステップS1では、CFDモデル解析手段10により、前回のサイクルで求めたRoom2の圧力および温度がRoom1の境界面のセルの圧力、温度、および流速に与える影響を算定する。具体的には、静圧、位置圧、前回のサイクルで求めたRoom1の流速から得られる動圧、およびサブグリッドスケールでの乱流エネルギーをまとめて、Room1における境界面の圧力とする。さらに、前回のサイクルで求めた境界面図2に示す開口での流入出の方向に応じて、Room2の温度をRoom1の境界面の温度とする。ステップS2では、CFDモデル解析手段10により、Room2が与える影響を考慮して、Room1の各セルの圧力を算定する。ステップS3では、CFDモデル解析手段10により、Room1の各セルの圧力値が収束しているか否かを判定する。この判定がYesである場合には、ステップS4に移り、Noである場合には、ステップS2に戻って再度各セルの圧力を算定する。ステップS4では、CFDモデル解析手段10により、Room1の各セルの流速および温度を算定する。ステップS5では、連成処理手段20により、前回のサイクルで求めたRoom2の空気層温度に基づいて、Room1の上部層と下部層との境界である層境界高さを算定する。ステップS6では、連成処理手段20により、ステップS5で求めたRoom1の層境界高さおよび前回のサイクルで求めたRoom2の層境界高さに基づいて、Room1とRoom2との境界面における質量流量およびエネルギー量を分配する。ステップS7では、二層ゾーンモデル解析手段30により、分配された質量流量およびエネルギー量に基づいて、Room2の圧力を算定する。ステップS8では、二層ゾーンモデル解析手段30により、Room2の圧力値が収束しているか否かを判定する。この判定がYesである場合には、ステップS9に移り、Noである場合には、ステップS7に戻る。ステップS9では、二層ゾーンモデル解析手段30により、Room2の層境界高さ、煙層および空気層の温度、ならびに圧力を算定する。次に、ステップS1における処理について説明する。すなわち、Room2の圧力および温度がRoom1の境界面のセルの圧力、温度、および流速に与える影響を、以下のように算定する。図4は、CFDモデルが適用されたRoom1と二層ゾーンモデルが適用されたRoom2との境界部分の概念図である。二層ゾーンモデルが適用されたRoom2がCFDモデルが適用されたRoom1の境界面のセルkに与える圧力をP|interface、流速をU|interface、温度をT|interfaceとする。セルの圧カP|interfaceは、PRoom2H、動圧ΔP、サブグリッドスケールでの乱流エネルギーkSGSを用いて、式1で表す。ここで、サブグリッドスケールでの乱流エネルギーkSGSとは、セル内に生じる乱流エネルギーを数値化したものである。以下、HはRoom2の床面からのセルの中心位置までの高さ、ρはセルの密度、Tは温度、Uは流速、nは流れの方向Room1からRoom2に向かう流れは、Room2からRoom1に向かう流れは-とする、δAはセルの境界面の面積である。また、添え字のkはセルの番号である。さらに、動圧ΔPについては、以下の式2で表される。また、PRoom2Hについては、静圧と位置圧に基づいて、以下の式3で表される。ここで、P0,Room2はRoom2の床面での静圧である。また、流速U|interfaceについては、以下の式4で表される。また、温度T|interfaceについては、以下の式5、式6で表される。次に、ステップS5における処理について説明する。すなわち、前回のサイクルで求めたRoom2の空気層温度に基づいて、Room1の上部層と下部層との境界である層境界高さを、以下の手順で算定する。まず、図5に示すように、CFDモデルが適用されたRoom1の境界面のセルを、高温セルまたは低温セルに分類する。具体的には、境界面のセルを特定の閾値を基準として分類するのではなく、セルの温度がRoom2の空気層の温度以上である場合、このセルを高温セルとし、セルの温度がRoom2の空気層の温度未満である場合、このセルを低温セルとする。次に、高温セルのみで構成される領域を高温領域、低温セルのみで構成される領域を低温領域、高温セルと低温セルとが混在する領域を温度混在領域とする。次に、Room1の層境界高さZRoom1を、高温セルの面積と低温セルの面積とで按分して決定する。すなわち、低温領域の最高高さをZ1、高温領域の最低高さをZ2とし、層境界高さZRoom1を以下の式7に従って求める。ここで、low_tempは低温セル、high_tempは高温セルである。次に、ステップS6における処理について説明する。すなわち、ステップS5の処理で求めたRoom1の層境界高さおよび前回のサイクルで求めたRoom2の層境界高さに基づいて、境界面における質量流量、エネルギー量を、以下の手順で分配する。Room1の層境界の高さZRoom1とRoom2の層境界の高さZRoom2によって、図6A、Bの2パターンが考えられる。図6A、Bに示すように、CFDモデルが適用されたRoom1は上部層および下部層に分割され、二層ゾーンモデルが適用されたRoom2は煙層および空気層に分割される。そして、表1に示すように、Room1の層境界の高さ、Room2の層境界の高さ、および煙の流れ方向の組合せで場合分けを行って、質量流量M、エネルギー量Uの流入出量を設定する。以下、Cpは定圧比熱である。また、添え字のUはRoom1の上部層またはRoom2の煙層、LはRoom1の下部層またはRoom2の空気層である。また、添え字の->はRoom1からRoom2への流出を示し、<-はRoom2からRoom1への流出を示す。表1に示す変数名の質量流量Mおよびエネルギー量Uを、式8~14に従って求める。また、Room1の上部層からRoom2の空気層に流入するcの流れは、噴流プルームであるため、この噴流プルームによる周辺空気の巻き込み量を、以下の手順で考慮する。図7は、噴流プルームの概念図である。噴流プルームの質量流量MPlumeを以下の式15に従って求める。また、噴流プルームのエネルギー量EPlumeを以下の式16に従って求める。以下、Qfは開口噴流プルームにおける点火源の発熱速度、Z0は点火源からRoom2の煙層下端までの高さ、ZはRoom1の上部層下端からRoom2の煙層下端までの高さである。さらに、Qf、Z0、Zは以下の式17~19に従って求められる。本実施形態によれば、以下のような効果がある。1Room1にCFDモデルを適用して解析することで、詳細な性状を把握する。一方、Room2に二層ゾーンモデルを適用して解析することで、大まかな性状を把握する。このように解析方法を使い分ける連成解析により、計算時間を短縮しつつ十分に詳細な解析結果を得ることができる。この連成解析では、各解析から得られる物理量を双方向で連成させた。具体的には、二層ゾーンモデル解析から得られる圧力および温度をCFDモデル解析に取り入れ、CFDモデル解析から得られる各セルの圧力、温度、および流速を質量流量およびエネルギー量の分配を介して二層ゾーンモデル解析に取り入れた。よって、連成解析を確実に実現できる。2Room1、Room2の境界面における質量流量およびエネルギー量を分配する際に、噴流プルームによる周辺空気の巻き込み量も考慮したので、連成解析の精度をより向上できる。3例えば、図8に示すような建物40に連成解析システム1を適用する。すなわち、図8は、建物40の断面図である。この建物40は、3層に設けられた居室41と、各居室41に連通する吹き抜け空間であるホール42と、を備える。各居室41とホール42との間には、廊下43が跳ね出している。また、廊下43とホール42との境界は防火区画44となっており、防火シャッタが設けられる。連成解析システム1により、ホール42にCFDモデルを適用し、各居室41および廊下43に二層ゾーンモデルを適用することで、連成解析により計算時間を短縮しつつ十分に詳細な解析結果を得ることができる。よって、ホール42内の煙の複雑な性状を把握できるので、防火シャッタなどの防火設備を適切な位置に配置したり居住者の避難経路を詳細に検討したりして、安全性を向上しつつ、施工コストを低減できる。なお、本発明は前記実施形態に限定されるものではなく、本発明の目的を達成できる範囲での変形、改良等は本発明に含まれるものである。例えば、本実施形態では、Room1で火災が発生し、この火災による煙がRoom2に流れ込む現象を解析したが、これに限らず、Room2で火災が発生し、この火災による煙がRoom1に流れ込む現象を解析することも可能である。また、本実施形態では、Room1、Room2の2室を解析対象としたが、これに限らず、3以上の室を解析対象としてもよい。"