建筑火灾数值仿真的区域分解优化方法技术

本发明专利技术提供了一种建筑火灾数值仿真的区域分解优化方法，包括步骤1：构建建筑模型，根据火焰的特征直径与经验公式计算出该建筑模型的全局网格尺寸，并对全局网格进行划分；步骤2：获取硬件配置；步骤3：根据负载均衡的影响因素将建筑模型划分为与机群中可使用的节点数目，或者超级计算机中可使用的CPU数目相同的子区域；步骤4：验证划分后各个子区域的负载均衡性，若均衡，则分配给相应节点进行计算；否则使用模式探索的方法调整子区域的体积，所述节点是指：计算机群中可使用的节点，或者可使用的CPU。本发明专利技术还能够自动检测每个计算机节点的负载均衡性并进行调整，适用于高性能超级计算机以及个体机群间，提高了数值模拟的速度。

【技术实现步骤摘要】
建筑火灾数值仿真的区域分解优化方法
本专利技术涉及建筑火灾防护
，具体地，涉及建筑火灾数值仿真的区域分解优化方法。
技术介绍
对于建筑火灾的数值模拟，由于计算空间极大，导致模拟时间较长，影响研究进度。而通过采用并行计算可以大幅度的缩短模拟时间。并行计算的基本思想是将一个大规模问题分解为若干小部分，并将每一部分交由不同的处理器进行处理，并在分区之间进行信息的传递。并行计算能力的提高一方面有赖于硬件技术的发展，另一方面也与分区策略等计算方法的改进密切相关。针对建筑火灾的数值模拟，除了使各子区域的单元数量大致相当，在复杂的火热耦合计算中，部分分区燃烧密集，流体流动频繁，这些都需要占用大量的计算资源，传统的一些分区方法无法保障各处理器间对火热模拟计算的负载均衡。因此，针对火热耦合系统，需要设计一种新的针对FDS软件的分区方法以提高并行计算的效率。传统的计算流体动力学分区方法与现有的网格划分技术运用在建筑火灾的数值模拟时，主要存在以下缺陷：1)没有给出合适的全局网格种子尺寸，导致对于一个案例的模拟往往需要进行多次的尝试，极大的浪费了时间；2)只考虑子区域网格数与通信时间，忽略了模拟过程中产生的火热耦合计算负载；3)FDS虽然具备了并行计算的功能，却没有提供相应的网格划分与分配的方法。
技术实现思路
针对现有技术中的缺陷，本专利技术的目的是提供一种建筑火灾数值仿真的区域分解优化方法。根据本专利技术提供的建筑火灾数值仿真的区域分解优化方法，包括如下步骤：步骤1：构建建筑模型，根据火焰的特征直径与经验公式计算出该建筑模型的全局网格尺寸，并根据实际情况和设置的模拟精度对全局网格进行划分；步骤2：获取硬件配置，所述硬件配置包括：计算机群中能够使用的节点，或者超级计算机中能够使用的CPU；步骤3：根据负载均衡的影响因素将建筑模型划分为与计算机群中能够使用的节点数目，或者超级计算机中能够使用的CPU数目相同的子区域；步骤4：验证划分后各个子区域的负载均衡性，若均衡，则分配给相应节点进行计算；若不均衡，则使用模式探索法调整子区域的体积，所述节点是指：计算机群中能够使用的节点，或者能够使用的CPU。优选地，所述步骤1中的火焰特征直径计算公式如下：式中：D*为火焰特征直径，Q为预计的热释放率，ρ∞为初始环境空气密度，cp为定压比热容，T∞为初始环境温度，g为当地重力加速度。优选地，所述步骤S1中的经验公式如下：4≤D*/Δd≤16式中：Δd为近似网格全局种子尺寸。优选地，所述步骤3中负载均衡的影响因素包括：子区域网格的数量、各个子区域之间的通信效率，以及由火与热耦合产生的计算负载。优选地，所述步骤3还包括：当能够使用的计算机群节点数为2的幂次方，则使用递归坐标二分法进行划分；当节点数不是2的幂次方的时，则按照各子区域网格数均衡原则与最小通信边界网格数原则进行划分；其中：各子区域网格数均衡原则与最小通信边界网格数原则的量化指标如下：σn＝nmodel/nsubσs＝ncom/nsubcom式中：σn为网格数负载平衡因子，nmodel为整个建筑模型所划分的网格数目，nsub为子区域的网格数目，σs为通信边界网格数负载平衡因子，ncom为所有子区域的通信网格数，nsubcom为某一子区域的通信网格数。优选地，所述步骤4包括：步骤4.1：对模拟过程中产生的火热耦合计算负载进行预测，预测公式如下：L＝f(d,n)式中：L为火热耦合计算负载，d为子区域与火源的距离，n为子区域内可燃物的数量；步骤4.2：运用试探的方法，即仅改变一个坐标方向上的子区域边界，每次移动应以改变最少的网格数为原则，剩余坐标方向上的边界保持不变；步骤4.3：对调整后的各个子区域进行负载均衡判断，若不均衡，则返回执行步骤4.2；若均衡，则通过相应节点进行计算；其中，负载均衡判断公式如下：σi＝σni+σsi+Li式中：σi为第i个子区域的总负载，σni为第i个子区域网格数负载，σsi为第i个子区域通信边界网格数负载，Li为第i个子区域火热耦合计算负载。优选地，所述步骤S4还包括：在刚开始运行的一段时间内对各个节点的运算时间进行监测，假设可供使用的节点数目为n，第i个节点的计算时间为ti，i的值为1,2,3…n，则平均时间记为平均时间的计算公式如下：将第i个节点相对平均计算时间的偏差记为计算公式如下：为设置一个上限用以评判负载是否均衡；当的值超过时，则认为负载不均衡，若的值小于等于时，则认为负载均衡；若不均衡则返回步骤4.2，并提高负载均衡的判断标准，重新调整子区域网格；若均衡则继续进行计算，直至模拟完成。与现有技术相比，本专利技术具有如下的有益效果：本专利技术提供的建筑火灾数值仿真的区域分解优化方法通过对建筑模型进行网格划分，并为每个网格配置相应的CPU进行并行计算；此外，还能够自动检测每个计算机节点的负载均衡性，并及时对负载进行调整，不仅适用于高性能超级计算机，也可以在个体机群间实施，大大提高数值模拟的速度。附图说明通过阅读参照以下附图对非限制性实施例所作的详细描述，本专利技术的其它特征、目的和优点将会变得更明显：图1为本专利技术提供的建筑火灾数值仿真的区域分解优化方法的流程示意图。具体实施方式下面结合具体实施例对本专利技术进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本专利技术，但不以任何形式限制本专利技术。应当指出的是，对本领域的普通技术人员来说，在不脱离本专利技术构思的前提下，还可以做出若干变化和改进。这些都属于本专利技术的保护范围。根据本专利技术提供的建筑火灾数值仿真的区域分解优化方法，包括如下步骤：步骤S1：构建建筑模型，根据火焰的特征直径与经验公式计算出该建筑模型的全局网格尺寸，并根据实际情况和设置的模拟精度对全局网格进行划分；具体地，构建建筑模型是指：输入建筑模型的基本参数，主要包括计算区域的大小，建筑材料属性，计算模型，计算方法等。步骤S2：获取硬件配置，所述硬件配置包括：计算机群中能够使用的节点，或者超级计算机中能够使用的CPU；步骤S3：根据负载均衡的影响因素将建筑模型划分为与计算机群中能够使用的节点数目，或者超级计算机中能够使用的CPU数目相同的子区域；步骤S4：验证划分后各个子区域的负载均衡性，若均衡，则分配给相应节点进行计算；若不均衡，则使用模式探索的方法调整子区域的体积，所述节点是指：计算机群中能够使用的节点，或者能够使用的CPU。所述步骤S1中的火焰特征直径计算公式如下：式中：D*为火焰特征直径，Q为预计的热释放率，ρ∞为初始环境空气密度，cp为定压比热容，T∞为初始环境温度，g为当地重力加速度。所述步骤S1中的经验公式如下：4≤D*/Δd≤16式中：Δd为近似网格全局种子尺寸。所述步骤S3中负载均衡的影响因素包括：子区域网格的数量、各个子区域之间的通信效率，以及由火与热耦合产生的计算负载。所述步骤S3还包括：当能够使用的计算机群节点数为2的幂次方，则使用递归坐标二分法进行划分；当节点数不是2的幂次方的时，则按照各子区域网格数均衡原则与最小通信边界网格数原则进行划分(必要的时候放弃少量节点数)；其中：各子区域网格数均衡原则与最小通信边界网格数原则的量化指标如下：σn＝nmodel/nsubσs＝ncom/nsubcom式中：σn为网格数负载平衡因子，nmodel为本文档来自技高网...

【技术保护点】
一种建筑火灾数值仿真的区域分解优化方法，其特征在于，包括如下步骤：步骤1：构建建筑模型，根据火焰的特征直径与经验公式计算出该建筑模型的全局网格尺寸，并根据实际情况和设置的模拟精度对全局网格进行划分；步骤2：获取硬件配置，所述硬件配置包括：计算机群中能够使用的节点，或者超级计算机中能够使用的CPU；步骤3：根据负载均衡的影响因素将建筑模型划分为与计算机群中能够使用的节点数目，或者超级计算机中能够使用的CPU数目相同的子区域；步骤4：验证划分后各个子区域的负载均衡性，若均衡，则分配给相应节点进行计算；若不均衡，则使用模式探索法调整子区域的体积，所述节点是指：计算机群中能够使用的节点，或者能够使用的CPU。

【技术特征摘要】
1.一种建筑火灾数值仿真的区域分解优化方法，其特征在于，包括如下步骤：步骤1：构建建筑模型，根据火焰的特征直径与经验公式计算出该建筑模型的全局网格尺寸，并根据实际情况和设置的模拟精度对全局网格进行划分；步骤2：获取硬件配置，所述硬件配置包括：计算机群中能够使用的节点，或者超级计算机中能够使用的CPU；步骤3：根据负载均衡的影响因素将建筑模型划分为与计算机群中能够使用的节点数目，或者超级计算机中能够使用的CPU数目相同的子区域；步骤4：验证划分后各个子区域的负载均衡性，若均衡，则分配给相应节点进行计算；若不均衡，则使用模式探索法调整子区域的体积，所述节点是指：计算机群中能够使用的节点，或者能够使用的CPU。2.根据权利要求1所述的建筑火灾数值仿真的区域分解优化方法，其特征在于，所述步骤1中的火焰特征直径计算公式如下：式中：D*为火焰特征直径，Q为预计的热释放率，ρ∞为初始环境空气密度，cp为定压比热容，T∞为初始环境温度，g为当地重力加速度。3.根据权利要求1所述的建筑火灾数值仿真的区域分解优化方法，其特征在于，所述步骤S1中的经验公式如下：4≤D*/Δd≤16式中：Δd为近似网格全局种子尺寸。4.根据权利要求1所述的建筑火灾数值仿真的区域分解优化方法，其特征在于，所述步骤3中负载均衡的影响因素包括：子区域网格的数量、各个子区域之间的通信效率，以及由火与热耦合产生的计算负载。5.根据权利要求1所述的建筑火灾数值仿真的区域分解优化方法，其特征在于，所述步骤3还包括：当能够使用的计算机群节点数为2的幂次方，则使用递归坐标二分法进行划分；当节点数不是2的幂次方的时，则按照各子区域网格数均衡原则与最小通信边界网格数原则进行划分；其中：各子区域网格数均衡原则与最小通信边界网格数原则的量化指标如下：σn＝nmodel/nsub...

【专利技术属性】
技术研发人员：江列霖，杨培中，史超，
申请(专利权)人：上海交通大学，
类型：发明
国别省市：上海,31