Movatterモバイル変換


[0]ホーム

URL:


CN103049612B - Building indoor environment optimization method based on model order reduction technology - Google Patents

Building indoor environment optimization method based on model order reduction technology
Download PDF

Info

Publication number
CN103049612B
CN103049612BCN201210564269.3ACN201210564269ACN103049612BCN 103049612 BCN103049612 BCN 103049612BCN 201210564269 ACN201210564269 ACN 201210564269ACN 103049612 BCN103049612 BCN 103049612B
Authority
CN
China
Prior art keywords
indoor environment
index
model
energy consumption
optimization
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.)
Expired - Fee Related
Application number
CN201210564269.3A
Other languages
Chinese (zh)
Other versions
CN103049612A (en
Inventor
李康吉
薛文平
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.)
Jiangsu University
Original Assignee
Jiangsu University
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 Jiangsu UniversityfiledCriticalJiangsu University
Priority to CN201210564269.3ApriorityCriticalpatent/CN103049612B/en
Publication of CN103049612ApublicationCriticalpatent/CN103049612A/en
Application grantedgrantedCritical
Publication of CN103049612BpublicationCriticalpatent/CN103049612B/en
Expired - Fee Relatedlegal-statusCriticalCurrent
Anticipated expirationlegal-statusCritical

Links

Images

Landscapes

Abstract

Translated fromChinese

本发明公开了一种基于模型降阶技术的建筑室内环境优化方法,包括三个主要步骤:1)利用CFD软件对室内环境作稳态仿真,并构造各类环境参数的变化空间;2)通过POD技术重构低阶的参数变化子空间;3)运用遗传算法搜索最优的空调送风温度和速度。本发明运用POD技术构造室内环境参数的变化子空间,从而在优化策略中能充分考虑空间分布对环境参数的影响,提高了优化的精确性。本征正交分解(POD)的模型降阶方法,将原空间内的控制方程映射到某个正交子空间内,且保证能量意义上映射误差最小。相对目前的环境优化策略,本发明具有优化精度高,速度快等优点。

The invention discloses a building indoor environment optimization method based on model reduction technology, which includes three main steps: 1) use CFD software to simulate the indoor environment in a steady state, and construct the change space of various environmental parameters; 2) through POD technology reconstructs the low-order parameter change subspace; 3) Use genetic algorithm to search for the optimal air-conditioning air supply temperature and speed. The invention utilizes the POD technology to construct the changing subspace of the indoor environment parameters, so that the influence of the space distribution on the environment parameters can be fully considered in the optimization strategy, and the accuracy of the optimization is improved. The model reduction method of intrinsic orthogonal decomposition (POD) maps the governing equations in the original space to an orthogonal subspace, and ensures the minimum mapping error in the energy sense. Compared with the current environment optimization strategy, the present invention has the advantages of high optimization precision, fast speed and the like.

Description

Translated fromChinese
一种基于模型降阶技术的建筑室内环境优化方法A Method for Building Indoor Environment Optimization Based on Model Reduction Technology

技术领域technical field

本发明涉及一种基于模型降阶技术的建筑室内环境优化方法,属于建筑环境与建筑节能领域。 The invention relates to a method for optimizing the indoor environment of a building based on model reduction technology, which belongs to the field of building environment and building energy saving. the

背景技术Background technique

随着人们对居住品质和建筑节能要求的不断提高,如何协调和优化建筑室内环境舒适度与空调能耗越来越受到关注。目前众多学者已提出多种系统级的多目标优化控制策略。优化算法从早期的依赖梯度的寻优方法发展到目前广泛运用的各类智能优化算法(如进化规划、遗传算法等),优化参数则涵盖了室内环境的各方面,包括热舒适度、空气质量以及空调能耗等。 With the continuous improvement of people's living quality and building energy-saving requirements, how to coordinate and optimize the indoor environment comfort of buildings and air-conditioning energy consumption has attracted more and more attention. At present, many scholars have proposed a variety of system-level multi-objective optimal control strategies. The optimization algorithm has developed from the early gradient-dependent optimization method to various intelligent optimization algorithms (such as evolutionary programming, genetic algorithm, etc.) widely used at present, and the optimization parameters cover all aspects of the indoor environment, including thermal comfort, air quality and air conditioning energy consumption. the

在环境优化控制策略中,如何针对候选控制变量快速准确地解算环境响应是一个核心问题。由于现成的建筑室内环境模型很难同时满足优化的实时性和精确度要求,目前通常的做法是假设室内空气完全混合,即忽略空间分布对环境参数的影响,采用经验模型或半机理模型的方法求解环境响应。而实际上,对于大多数空调系统,特别像置换通风系统来说,室内的环境参数在空间上有较大差异。忽略这种差异会导致优化效果与室内各区域人员的实际感受不符,引起各种舒适度抱怨。目前国际上这方面研究很有限。原因是室内多参数环境建模复杂,须借助 CFD 工具,且很难直接与在线的优化控制算法整合。 In the environmental optimization control strategy, how to quickly and accurately solve the environmental response to the candidate control variables is a core issue. Since it is difficult for the ready-made building indoor environment model to meet the real-time and accuracy requirements of optimization at the same time, the current common practice is to assume that the indoor air is completely mixed, that is, ignoring the influence of spatial distribution on environmental parameters, and adopting empirical models or semi-mechanistic models. Solve for the environmental response. In fact, for most air-conditioning systems, especially displacement ventilation systems, the indoor environmental parameters vary greatly in space. Neglecting this difference will lead to a discrepancy between the optimization effect and the actual feelings of people in various areas of the room, causing various comfort complaints. At present, international research in this area is very limited. The reason is that the modeling of the indoor multi-parameter environment is complex, and CFD tools must be used, and it is difficult to directly integrate with the online optimization control algorithm. the

2009年,有文献提出通过神经网络训练的方法得到CFD模型的简化模型,用于环境优化策略中参数指标的快速求解。此方法通过CFD仿真得到足够多的输入/输出数据。通过这些数据对的训练和测试建立基于神经网络的环境指标替代模型。在优化算法的每次迭代中利用替代模型快速求解目标函数,减少优化算法的复杂度,提高实时性。此方法考虑了环境参数的空间分布影响,但是神经网络模型本质上是经验模型,只能对指定的环境指标作建模。当性能指标有变化,或关注的室内用户区域改变,则须重新建模,应变能力较差。 In 2009, some literature proposed to obtain a simplified model of the CFD model through the method of neural network training, which is used to quickly solve the parameter indicators in the environment optimization strategy. This method obtains enough input/output data through CFD simulation. Through the training and testing of these data pairs, a neural network-based environmental indicator surrogate model is established. In each iteration of the optimization algorithm, the surrogate model is used to quickly solve the objective function, reducing the complexity of the optimization algorithm and improving real-time performance. This method takes into account the influence of the spatial distribution of environmental parameters, but the neural network model is essentially an empirical model and can only model specified environmental indicators. When the performance index changes, or the indoor user area of concern changes, it must be remodeled, and the adaptability is poor. the

发明内容Contents of the invention

针对现有建筑室内环境优化方法所存在的上述缺陷,本发明提供一种基于模型降阶技术的建筑室内环境优化方法。其特点在于通过构造低阶的环境参数变化子空间,将相关的室内环境参数模型直接嵌入寻优过程中,实现环境响应的精确、快速解算。 Aiming at the above-mentioned defects in existing building indoor environment optimization methods, the present invention provides a building indoor environment optimization method based on model reduction technology. Its feature is that by constructing a low-order environmental parameter change subspace, the relevant indoor environmental parameter models are directly embedded in the optimization process to achieve accurate and fast calculation of environmental responses. the

基于本征正交分解(POD)的模型降阶方法本质上是一种映射方法。它将原空间内的控制方程映射到某个正交子空间内,且保证能量意义上映射误差最小。典型的CFD 模型通过POD 降阶方法可转换为关于POD模式系数的低阶模型,能同时满足建模精度与实时性要求,适合优化策略中环境响应的快速解算任务。 The model reduction method based on intrinsic orthogonal decomposition (POD) is essentially a mapping method. It maps the governing equations in the original space to an orthogonal subspace, and ensures the minimum mapping error in the energy sense. A typical CFD model can be converted into a low-order model about POD model coefficients through the POD order reduction method, which can meet the modeling accuracy and real-time requirements at the same time, and is suitable for the fast solution task of the environmental response in the optimization strategy. the

本发明的技术方案是: Technical scheme of the present invention is:

一种基于模型降阶技术的建筑室内环境优化方法,包括如下步骤:A method for optimizing the indoor environment of a building based on model reduction technology, comprising the following steps:

(1)建立基于CFD的室内环境模型;(1) Establish an indoor environment model based on CFD;

(2)根据控制变量可能的变化范围,等间距选择控制变量数据点,作相应的CFD稳态仿真;从CFD仿真结果中提取环境参数分布,构造参数变化空间;提取的参数类型包括室内温度、风速、污染物浓度及热舒适度指标;(2) According to the possible variation range of the control variable, select the data points of the control variable at equal intervals, and perform the corresponding CFD steady-state simulation; extract the distribution of environmental parameters from the CFD simulation results, and construct the parameter variation space; the extracted parameter types include indoor temperature, Wind speed, pollutant concentration and thermal comfort index;

(3)利用POD模型降阶技术重构出步骤(2)所得参数变化空间的低阶子空间;(3) Reconstruct the low-order subspace of the parameter variation space obtained in step (2) by using the POD model reduction technology;

(4)选择室内环境指标及能耗指标,用于评估室内环境及空调能耗;(4) Select indoor environment indicators and energy consumption indicators for evaluating indoor environment and air-conditioning energy consumption;

(5)设置目标函数,利用优化算法对控制变量进行迭代优化。在每次优化迭代过程中,通过在参数子空间内的多维插值快速得到系统响应,进而快速求解目标函数。(5) Set the objective function, and use the optimization algorithm to iteratively optimize the control variables. In each optimization iteration process, the system response is quickly obtained through multi-dimensional interpolation in the parameter subspace, and then the objective function is quickly solved.

所述步骤(1)中,CFD仿真使用Airpak计算流体力学软件;室内环境模型为三维模型;室内环境模型的建立步骤如下: In the step (1), the CFD simulation uses Airpak computational fluid dynamics software; the indoor environment model is a three-dimensional model; the steps for establishing the indoor environment model are as follows:

A、利用Airpak软件建立房间围护的几何模型;确定空调送风口和回风口的位置与尺寸;确定房间内主要陈设的位置与尺寸;A. Use Airpak software to establish the geometric model of the room enclosure; determine the position and size of the air supply outlet and return air outlet of the air conditioner; determine the position and size of the main furnishings in the room;

B、对建立的房间模型划分网格;B. Mesh the established room model;

C、利用Fluent求解器耦合求取质量、动量、能量及污染物浓度方程的稳态解。C. Use the Fluent solver coupling to obtain the steady-state solutions of mass, momentum, energy and pollutant concentration equations.

所述步骤(2)中,控制变量包括空调送风口温度和风速;控制变量的选择间隔为:送风口温度0.1摄氏度,送风口风速0.1米/秒;对每一组控制变量利用步骤(1)所述的室内环境模型进行稳态仿真;通过Airpak软件的导出功能提取每一组控制变量对应的稳态参数分布,组成各类参数的变化空间。 In the step (2), the control variables include the temperature and wind speed of the air outlet of the air conditioner; the selection interval of the control variables is: the temperature of the air outlet is 0.1 degrees Celsius, and the wind speed of the air outlet is 0.1 m/s; for each group of control variables, use the step (1) The indoor environment model is used for steady-state simulation; the steady-state parameter distribution corresponding to each group of control variables is extracted through the export function of the Airpak software to form the variation space of various parameters. the

所述步骤(3)中,POD模型降阶的基本思想为:在n维向量空间H中有一组数据集,找到其一组m维子集构成S子空间(m<n),使原数据集映射到子集的误差在能量意义上最小。POD模型降阶的基本步骤如下: In the step (3), the basic idea of POD model order reduction is: there is a set of data sets in the n-dimensional vector space H, and a set of m-dimensional subsets are found to form an S subspace (m<n), so that the original data The set-to-subset mapping error is minimal in an energy sense. The basic steps of POD model reduction are as follows:

A、利用各参数变化空间组成矩阵:A. Use the variable space of each parameter to form a matrix:

                         

Figure 988506DEST_PATH_IMAGE001
                         
Figure 988506DEST_PATH_IMAGE001

这里

Figure 967964DEST_PATH_IMAGE002
维的矩阵表示参数变化空间,
Figure 18276DEST_PATH_IMAGE004
是房间内网格点总数,是控制变量组数;
Figure 530477DEST_PATH_IMAGE006
Figure 607630DEST_PATH_IMAGE003
的转置;here
Figure 967964DEST_PATH_IMAGE002
dimensional The matrix represents the parameter variation space,
Figure 18276DEST_PATH_IMAGE004
is the total number of grid points in the room, is the number of control variable groups;
Figure 530477DEST_PATH_IMAGE006
for
Figure 607630DEST_PATH_IMAGE003
the transposition of

B、求解矩阵

Figure 806530DEST_PATH_IMAGE007
的特征值和特征向量
Figure 344139DEST_PATH_IMAGE009
;选择合适的截断值
Figure 160785DEST_PATH_IMAGE010
,使得前
Figure 710847DEST_PATH_IMAGE010
个特征值包含的系统动能占比
Figure 296549DEST_PATH_IMAGE011
大于99%,这里
Figure 726393DEST_PATH_IMAGE011
表示为:B. Solve the matrix
Figure 806530DEST_PATH_IMAGE007
The eigenvalues of and eigenvectors
Figure 344139DEST_PATH_IMAGE009
; Choose an appropriate cutoff value
Figure 160785DEST_PATH_IMAGE010
, so that before
Figure 710847DEST_PATH_IMAGE010
The proportion of system kinetic energy contained in each eigenvalue
Figure 296549DEST_PATH_IMAGE011
Greater than 99%, here
Figure 726393DEST_PATH_IMAGE011
Expressed as:

                        

Figure 781068DEST_PATH_IMAGE012
                        
Figure 781068DEST_PATH_IMAGE012

C、低阶参数变化子空间可描述为

Figure 56191DEST_PATH_IMAGE010
个特征向量及其系数的线性组合:C. The low-order parameter variation subspace can be described as
Figure 56191DEST_PATH_IMAGE010
A linear combination of eigenvectors and their coefficients:

                      

Figure 761979DEST_PATH_IMAGE013
                      
Figure 761979DEST_PATH_IMAGE013

这里

Figure 907265DEST_PATH_IMAGE014
表示第
Figure 636187DEST_PATH_IMAGE015
组控制变量对应的稳态参数分布,
Figure 777318DEST_PATH_IMAGE016
表示POD模式系数。here
Figure 907265DEST_PATH_IMAGE014
Indicates the first
Figure 636187DEST_PATH_IMAGE015
The steady-state parameter distribution corresponding to the group control variables,
Figure 777318DEST_PATH_IMAGE016
Indicates the POD mode coefficient.

所述步骤(4)中,室内环境指标包括热舒适度指标及室内空气质量指标。 In the step (4), the indoor environment index includes a thermal comfort index and an indoor air quality index. the

热舒适度指标采用预测平均投票指标(PMV)。PMV 指标将人体冷热感觉量化为以下七级:冷(-3)、凉(-2)、稍凉(-1)、舒适(0)、稍暖(1)、暖(2)、热(3),并将其与空气温度、太阳辐射、空气流速、空气湿度、人体新陈代谢率、及人体着衣等六个因素用函数联系起来,是目前国际上最为通用的热舒适度定量指标。 The thermal comfort index adopts the predictive mean voting index (PMV). The PMV index quantifies the human body's hot and cold sensations into the following seven levels: cold (-3), cool (-2), slightly cool (-1), comfortable (0), slightly warm (1), warm (2), hot ( 3), and linking it with six factors such as air temperature, solar radiation, air velocity, air humidity, human metabolic rate, and human clothing, is currently the most commonly used quantitative index of thermal comfort in the world. the

室内空气质量指标采用通风效力指标(

Figure 275296DEST_PATH_IMAGE017
): Indoor air quality index adopts ventilation efficiency index (
Figure 275296DEST_PATH_IMAGE017
):

                       

Figure 859992DEST_PATH_IMAGE018
                       
Figure 859992DEST_PATH_IMAGE018

这里,

Figure 872947DEST_PATH_IMAGE019
Figure 755453DEST_PATH_IMAGE020
分别为空调回风口和送风口的污染物浓度,
Figure 920986DEST_PATH_IMAGE021
为室内人员头部高度的污染物平均浓度;here,
Figure 872947DEST_PATH_IMAGE019
and
Figure 755453DEST_PATH_IMAGE020
are the pollutant concentrations at the return air outlet and air supply outlet of the air conditioner, respectively,
Figure 920986DEST_PATH_IMAGE021
is the average concentration of pollutants at the head height of indoor personnel;

本发明将空调系统能耗分解为两个部分:风机能耗和制冷能耗,并作适当简化,得到能耗指标为:The present invention decomposes the energy consumption of the air-conditioning system into two parts: fan energy consumption and refrigeration energy consumption, and appropriately simplifies it to obtain the energy consumption index as follows:

                     

Figure 863534DEST_PATH_IMAGE022
                     
Figure 863534DEST_PATH_IMAGE022

这里,

Figure 363786DEST_PATH_IMAGE023
为风机能耗,
Figure 784403DEST_PATH_IMAGE024
为风机压升,为总送风量,
Figure 245788DEST_PATH_IMAGE026
为制冷能耗,
Figure 171019DEST_PATH_IMAGE027
为用于除显热负载的制冷能耗,为对新风进行除湿降温的能耗;here,
Figure 363786DEST_PATH_IMAGE023
is the fan energy consumption,
Figure 784403DEST_PATH_IMAGE024
is the fan pressure rise, is the total air volume,
Figure 245788DEST_PATH_IMAGE026
For cooling energy consumption,
Figure 171019DEST_PATH_IMAGE027
is the cooling energy used to remove sensible heat loads, Energy consumption for dehumidification and cooling of fresh air;

所述步骤(5)中,目标函数设置为:In the step (5), the objective function is set as:

Figure 266943DEST_PATH_IMAGE029
Figure 266943DEST_PATH_IMAGE029

这里,

Figure 613611DEST_PATH_IMAGE030
为各指标的加权系数,下标
Figure 26138DEST_PATH_IMAGE031
指对应性能指标的最大值,用于各个指标的归一化;
Figure 867186DEST_PATH_IMAGE032
为惩罚项,用于反映室内人员头脚温差和周边风速过大对热舒适度的影响。here,
Figure 613611DEST_PATH_IMAGE030
is the weighting coefficient of each indicator, subscript
Figure 26138DEST_PATH_IMAGE031
Refers to the maximum value of the corresponding performance index, which is used for the normalization of each index;
Figure 867186DEST_PATH_IMAGE032
As a penalty item, it is used to reflect the influence of the temperature difference between the head and feet of the indoor occupants and the excessive wind speed around them on the thermal comfort.

优化算法采用遗传算法,参数子空间内的多维插值算法采用样条插值。 The optimization algorithm adopts genetic algorithm, and the multi-dimensional interpolation algorithm in parameter subspace adopts spline interpolation. the

本发明提出一种基于模型降阶技术的建筑室内环境优化方法,充分考虑空间分布对环境指标的影响,将“室内环境”降阶之后直接嵌入环境优化算法,以满足优化的精确性和实时性要求。 The present invention proposes a building indoor environment optimization method based on model reduction technology, which fully considers the influence of spatial distribution on environmental indicators, and directly embeds the "indoor environment" into the environment optimization algorithm after the reduction of the order to meet the accuracy and real-time optimization Require. the

相对目前的环境优化方法,本发明的优点表现在: Compared with the current environment optimization method, the advantages of the present invention are as follows:

1、优化的精确性。1. The accuracy of optimization.

本发明不再假设室内空气“充分混合”,而是利用CFD工具对建筑室内环境做精确建模,并通过模型降阶的方法将得到的低阶“室内环境”嵌入环境优化策略中,使得优化结果更精确。在体育馆,宾馆大厅,医院,学校等大空间场合,本发明提出的方法尤其具备明显的精度优势。 The present invention no longer assumes that the indoor air is "fully mixed", but uses CFD tools to accurately model the indoor environment of the building, and embeds the obtained low-order "indoor environment" into the environment optimization strategy through the method of model reduction, so that the optimization The result is more precise. In large space occasions such as gymnasiums, hotel halls, hospitals, schools, etc., the method proposed by the present invention has obvious advantages in precision. the

2、优化的快速性。 2. The speed of optimization. the

本发明通过对参数子空间作多维插值可快速求解目标函数,使得优化算法在提高精度的同时满足实时性要求,可实际运用于在线优化场合。 The invention can quickly solve the objective function by performing multi-dimensional interpolation on the parameter subspace, so that the optimization algorithm can meet the real-time requirement while improving the precision, and can be actually applied to on-line optimization occasions. the

附图说明Description of drawings

图1是一个3D办公室模型示意图; Figure 1 is a schematic diagram of a 3D office model;

图2 是优化前的室内污染物浓度分布情况图;Figure 2 is a diagram of the indoor pollutant concentration distribution before optimization;

图3是优化后的室内污染物浓度分布情况图。Figure 3 is a diagram of the optimized indoor pollutant concentration distribution.

具体实施方式Detailed ways

为了更为具体地描述本发明,下面结合附图和具体实施例对本发明进行详细说明。 In order to describe the present invention more specifically, the present invention will be described in detail below in conjunction with the accompanying drawings and specific embodiments. the

图1为一个3D办公室模型示意图。该办公室的长宽高分别为5.16m×3.65m×2.44m。室内设计有两个坐姿的办公人员(c,d)、两张办公桌(g,h)、两台电脑(d,e)、两个文件柜(j,k)及六盏日光灯(l - q)。房间左侧墙有一面3.65m×1.04m的窗户(i),置换通风系统的送风口(a)设置在窗的对面墙上,回风口(b)布置在天花板中心位置。空气污染物的挥发口设计于人员的头部位置。 Figure 1 is a schematic diagram of a 3D office model. The length, width and height of the office are 5.16m×3.65m×2.44m respectively. The interior design has two sitting office workers (c, d), two desks (g, h), two computers (d, e), two filing cabinets (j, k) and six fluorescent lamps (l - q). There is a 3.65m×1.04m window (i) on the left wall of the room, the air supply outlet (a) of the displacement ventilation system is set on the wall opposite to the window, and the return air outlet (b) is arranged in the center of the ceiling. The volatilization port of air pollutants is designed at the head position of the personnel. the

下面详细描述本发明方法的实施步骤: Describe the implementation steps of the inventive method in detail below:

步骤0. 建立基于CFD的室内环境模型。建模方法为:Step 0. Establish a CFD-based indoor environment model. The modeling method is:

步骤0.1 确定包括房间的围护(墙、地板及天花板等)、空调送风口/回风口以及室内陈设(包括人员)的位置与尺寸,利用Airpak软件搭建如图1所示的办公室几何模型;Step 0.1 Determine the location and size of the enclosure including the room (wall, floor and ceiling, etc.), air-conditioning air supply/return vents, and interior furnishings (including personnel), and use Airpak software to build the office geometric model as shown in Figure 1;

步骤0.2 对建立的几何模型划分网格,本例中共划分72282个不规则网格;Step 0.2 Mesh the established geometric model. In this example, a total of 72282 irregular meshes are divided;

步骤0.3 利用Fluent求解器耦合求取质量、动量、能量及污染物浓度方程的稳态解;求解之前,相关的边界条件设置如下:空调送风口设置为速度入口边界;回风口设置为自然流出边界;墙、地面及天花板设置为温度边界。相关的模型定义和求解策略设置如下:室内气体假设为低速流动的不可压缩粘性牛顿流体,湍流模型选用

Figure 48768DEST_PATH_IMAGE033
标准模型,近壁处理采用标准壁面函数,浮力效应采用 Boussinesq 近似方式,不考虑粘性发热,压力速度耦合计算采用 SIMPLE 算法。Step 0.3 Use the Fluent solver coupling to obtain the steady-state solution of the mass, momentum, energy, and pollutant concentration equations; before solving, the relevant boundary conditions are set as follows: the air-conditioning air supply port is set as the velocity inlet boundary; the air return port is set as the natural outflow boundary ; Walls, floors, and ceilings are set as temperature boundaries. The relevant model definition and solution strategy settings are as follows: the indoor gas is assumed to be an incompressible viscous Newtonian fluid flowing at a low speed, and the turbulence model is selected as
Figure 48768DEST_PATH_IMAGE033
In the standard model, the standard wall function is used for the near-wall treatment, the Boussinesq approximation method is used for the buoyancy effect, viscous heating is not considered, and the SIMPLE algorithm is used for the pressure-velocity coupling calculation.

步骤1. 控制变量在本例中的变化范围为:空调送风口温度:17—21摄氏度,空调送风口速度:0.1—0.5米/秒。等间距(0.1度以及0.1米/秒)地选择共25个控制变量数据点,利用步骤0.3作相应的CFD稳态仿真;从CFD仿真结果中提取环境参数分布。提取的参数类型包括室内温度、风速、污染物浓度及热舒适度指标。对每类参数构造

Figure 300758DEST_PATH_IMAGE002
维的参数变化空间,其中n为网格点总数(本例为72282),m为控制变量数据点总数(本例为25); Step 1. The variation range of the control variables in this example is: air-conditioning outlet temperature: 17-21 degrees Celsius, air-conditioning outlet speed: 0.1-0.5 m/s. A total of 25 control variable data points are selected at equal intervals (0.1 degree and 0.1 m/s), and the corresponding CFD steady-state simulation is performed using step 0.3; the environmental parameter distribution is extracted from the CFD simulation results. The extracted parameter types include indoor temperature, wind speed, pollutant concentration and thermal comfort index. For each type of parameter construction
Figure 300758DEST_PATH_IMAGE002
Dimensional parameter change space, where n is the total number of grid points (72282 in this example), and m is the total number of control variable data points (25 in this example);

步骤2.利用POD模型降阶技术重构出步骤1所得参数变化空间的低阶子空间。模型降阶步骤为:Step 2. Use the POD model reduction technology to reconstruct the low-order subspace of the parameter variation space obtained instep 1. The model reduction steps are:

步骤2.1 利用各参数变化空间组成矩阵:Step 2.1 Use each parameter change space to form a matrix:

Figure 200581DEST_PATH_IMAGE001
Figure 200581DEST_PATH_IMAGE001

这里

Figure 579741DEST_PATH_IMAGE002
维的
Figure 881409DEST_PATH_IMAGE003
矩阵表示参数变化空间,
Figure 304300DEST_PATH_IMAGE006
Figure 504469DEST_PATH_IMAGE003
的转置;here
Figure 579741DEST_PATH_IMAGE002
dimensional
Figure 881409DEST_PATH_IMAGE003
The matrix represents the parameter variation space,
Figure 304300DEST_PATH_IMAGE006
for
Figure 504469DEST_PATH_IMAGE003
the transposition of

步骤2.2 求解矩阵的特征值

Figure 827183DEST_PATH_IMAGE008
和特征向量;选择合适的截断值
Figure 43193DEST_PATH_IMAGE010
,使得前
Figure 216685DEST_PATH_IMAGE010
个特征值包含的系统动能占比
Figure 289684DEST_PATH_IMAGE011
大于99%,其中表示为:Step 2.2 Solve the matrix The eigenvalues of
Figure 827183DEST_PATH_IMAGE008
and eigenvectors ; Choose an appropriate cutoff value
Figure 43193DEST_PATH_IMAGE010
, so that before
Figure 216685DEST_PATH_IMAGE010
The proportion of system kinetic energy contained in each eigenvalue
Figure 289684DEST_PATH_IMAGE011
greater than 99%, of which Expressed as:

Figure 166821DEST_PATH_IMAGE012
Figure 166821DEST_PATH_IMAGE012

本例中,当

Figure 878425DEST_PATH_IMAGE034
时,系统动能占比
Figure 805930DEST_PATH_IMAGE011
即大于99%。In this example, when
Figure 878425DEST_PATH_IMAGE034
When , the kinetic energy of the system accounts for
Figure 805930DEST_PATH_IMAGE011
That is greater than 99%.

步骤2.3 低阶参数变化子空间可描述为个特征向量及其系数的线性组合: Step 2.3 The low-order parameter change subspace can be described as A linear combination of eigenvectors and their coefficients:

这里

Figure 590980DEST_PATH_IMAGE014
表示第
Figure 638570DEST_PATH_IMAGE015
组控制变量对应的稳态参数分布,
Figure 948329DEST_PATH_IMAGE016
表示特征向量的系数或称为POD模式系数;here
Figure 590980DEST_PATH_IMAGE014
Indicates the first
Figure 638570DEST_PATH_IMAGE015
The steady-state parameter distribution corresponding to the group control variables,
Figure 948329DEST_PATH_IMAGE016
Represents the coefficient of the eigenvector or called the POD mode coefficient;

步骤3. 选择室内环境指标及能耗指标,用于评估室内环境及空调能耗;Step 3. Select indoor environment indicators and energy consumption indicators for evaluating indoor environment and air-conditioning energy consumption;

室内环境指标包括热舒适度指标及室内空气质量指标;Indoor environmental indicators include thermal comfort indicators and indoor air quality indicators;

热舒适度指标采用预测平均投票指标(PMV)。PMV 指标将人体冷热感觉量化为以下七级:冷(-3)、凉(-2)、稍凉(-1)、舒适(0)、稍暖(1)、暖(2)、热(3),并将其与空气温度、太阳辐射、空气流速、空气湿度、人体新陈代谢率、及人体着衣等六个因素用函数联系起来,是目前国际上最为通用的热舒适度定量指标。The thermal comfort index adopts the predictive mean voting index (PMV). The PMV index quantifies the human body's hot and cold sensations into the following seven levels: cold (-3), cool (-2), slightly cool (-1), comfortable (0), slightly warm (1), warm (2), hot ( 3), and linking it with six factors such as air temperature, solar radiation, air velocity, air humidity, human metabolic rate, and human clothing, is currently the most commonly used quantitative index of thermal comfort in the world.

室内空气质量指标采用通风效力指标(

Figure 832103DEST_PATH_IMAGE017
): Indoor air quality index adopts ventilation efficiency index (
Figure 832103DEST_PATH_IMAGE017
):

                       

Figure 947826DEST_PATH_IMAGE018
                       
Figure 947826DEST_PATH_IMAGE018

这里,

Figure 787606DEST_PATH_IMAGE019
分别为空调回风口和送风口的污染物浓度,
Figure 370827DEST_PATH_IMAGE021
为室内人员头部高度的污染物平均浓度;here,
Figure 787606DEST_PATH_IMAGE019
and are the pollutant concentrations at the return air outlet and air supply outlet of the air conditioner, respectively,
Figure 370827DEST_PATH_IMAGE021
is the average concentration of pollutants at the head height of indoor personnel;

本发明将空调系统能耗分解为两个部分:风机能耗和制冷能耗,并作适当简化,得到能耗指标为:The present invention decomposes the energy consumption of the air-conditioning system into two parts: fan energy consumption and refrigeration energy consumption, and appropriately simplifies it to obtain the energy consumption index as follows:

Figure 290241DEST_PATH_IMAGE022
Figure 290241DEST_PATH_IMAGE022

这里,

Figure 984528DEST_PATH_IMAGE023
为风机能耗,
Figure 449138DEST_PATH_IMAGE024
为风机压升,
Figure 494455DEST_PATH_IMAGE025
为总送风量,
Figure 951981DEST_PATH_IMAGE026
为制冷能耗,
Figure 500774DEST_PATH_IMAGE027
为用于除显热负载的制冷能耗,
Figure 401865DEST_PATH_IMAGE028
为对新风进行除湿降温的能耗;here,
Figure 984528DEST_PATH_IMAGE023
is the fan energy consumption,
Figure 449138DEST_PATH_IMAGE024
is the fan pressure rise,
Figure 494455DEST_PATH_IMAGE025
is the total air volume,
Figure 951981DEST_PATH_IMAGE026
For cooling energy consumption,
Figure 500774DEST_PATH_IMAGE027
is the cooling energy used to remove sensible heat loads,
Figure 401865DEST_PATH_IMAGE028
Energy consumption for dehumidification and cooling of fresh air;

步骤4. 利用优化算法对控制变量进行迭代优化。本例中的优化算法为遗传算法,设置目标函数为环境及能耗指标的加权多项式,Step 4. Use the optimization algorithm to iteratively optimize the control variables. The optimization algorithm in this example is a genetic algorithm, and the objective function is set as a weighted polynomial of environmental and energy consumption indicators.

Figure 668898DEST_PATH_IMAGE029
Figure 668898DEST_PATH_IMAGE029

这里,

Figure 930115DEST_PATH_IMAGE030
为各指标的加权系数,下标
Figure 333415DEST_PATH_IMAGE031
指对应性能指标的最大值,用于各个指标的归一化;为惩罚项,用于反映室内人员头脚温差和周边风速过大对热舒适度的影响;here,
Figure 930115DEST_PATH_IMAGE030
is the weighting coefficient of each indicator, subscript
Figure 333415DEST_PATH_IMAGE031
Refers to the maximum value of the corresponding performance index, which is used for the normalization of each index; As a penalty item, it is used to reflect the influence of the temperature difference between the head and feet of the indoor personnel and the excessive wind speed around them on the thermal comfort;

在遗传算法的每次迭代过程中,通过在参数子空间内的多维插值快速得到系统响应,进而求解目标函数。这里的多维插值算法采用样条插值算法。In each iteration of the genetic algorithm, the system response is quickly obtained through multi-dimensional interpolation in the parameter subspace, and then the objective function is solved. The multidimensional interpolation algorithm here adopts the spline interpolation algorithm.

为了对比优化效果,选择空调送风口温度17摄氏度,速度0.1米/秒作为优化前的室内环境基础案例。应用本优化方法,空调能耗最多可降低41.2%,热舒适度最多可提高48.6%,室内空气质量最多可改善38%。 In order to compare the optimization effect, the temperature of the air outlet of the air conditioner is 17 degrees Celsius and the speed is 0.1 m/s as the basic case of the indoor environment before optimization. Applying this optimization method, the air-conditioning energy consumption can be reduced by up to 41.2%, the thermal comfort can be increased by up to 48.6%, and the indoor air quality can be improved by up to 38%. the

图2和图3 为优化前后室内空气污染物浓度的对比图。图2描述了优化前的室内污染物浓度分布情况,图3描述了优化后的室内污染物浓度分布情况。由图可见,本发明通过在优化策略中充分考虑环境参数的空间分布情况,能够明显改善通风效果,提高室内空气质量。 Figure 2 and Figure 3 are the comparison charts of indoor air pollutant concentrations before and after optimization. Figure 2 describes the distribution of indoor pollutant concentration before optimization, and Figure 3 describes the distribution of indoor pollutant concentration after optimization. It can be seen from the figure that the present invention can significantly improve the ventilation effect and improve the indoor air quality by fully considering the spatial distribution of environmental parameters in the optimization strategy. the

上面已经结合具体实施步骤说明了本发明,然而对于本领域的技术人员来说,可以在不背离本发明的精神和范围的前提下,对本发明做出不同的改进和变型。因而落入本发明的权利要求范围内的各种改进和变型,都应属于本发明的保护范围之内。 The present invention has been described above in conjunction with specific implementation steps. However, for those skilled in the art, various improvements and modifications can be made to the present invention without departing from the spirit and scope of the present invention. Therefore, various improvements and modifications falling within the scope of the claims of the present invention shall fall within the protection scope of the present invention. the

Claims (6)

1. A building indoor environment optimization method based on a model order reduction technology specifically comprises the following steps:
(1) establishing an indoor environment model based on CFD;
(2) selecting control variable data points at equal intervals according to the possible variation range of the control variable, and performing corresponding CFD steady-state simulation; extracting environmental parameter distribution from a CFD simulation result, and constructing a parameter change space;
(3) reconstructing a low-order subspace of the parameter change space obtained in the step (2) by using a POD model order reduction technology;
(4) selecting an indoor environment index and an energy consumption index for evaluating the indoor environment and the energy consumption of the air conditioner;
(5) and setting an objective function, and performing iterative optimization on the control variable by using an optimization algorithm.
2. The method for optimizing the indoor environment of the building based on the model order reduction technology as claimed in claim 1, wherein: in the step (1), Airpak computational fluid dynamics software is used for CFD simulation; the indoor environment model is a three-dimensional model; the indoor environment model is established by the following steps:
(1) establishing a geometric model of the room enclosure by using Airpak software; determining the positions and the sizes of an air supply outlet and an air return inlet of the air conditioner; determining the position and size of a main display in a room;
(2) dividing grids for the established room model;
(3) and (4) solving the steady state solution of the mass, momentum, energy and pollutant concentration equation by using the Fluent solver coupling.
3. The method for optimizing the indoor environment of the building based on the model order reduction technology as claimed in claim 1, wherein: in the step (2), the control variables comprise the temperature and the wind speed of an air conditioner air supply outlet; the selected intervals of the control variables are: the temperature of the air supply outlet is 0.1 ℃, and the air speed of the air supply outlet is 0.1 m/s; performing steady-state simulation on each set of control variables by using the indoor environment model in the step (1); extracting the steady-state parameter distribution corresponding to each set of control variable through the derivation function of the Airpak software to form a variation space of various parameters; the extracted environmental parameter types comprise indoor temperature, wind speed, pollutant concentration and thermal comfort level indexes.
4. The method for optimizing the indoor environment of the building based on the model order reduction technology as claimed in claim 1, wherein: in the step (3), the basic steps of POD model order reduction are as follows:
(1) forming a matrix by using the variation space of each parameter:
Figure 926830DEST_PATH_IMAGE001
here, the
Figure 56329DEST_PATH_IMAGE002
Of dimension
Figure 120362DEST_PATH_IMAGE003
The matrix represents the space of variation of the parameters,
Figure 478662DEST_PATH_IMAGE004
is the total number of grid points in the room,
Figure 138183DEST_PATH_IMAGE005
is the number of control variable groups;
Figure 189315DEST_PATH_IMAGE006
is composed of
Figure 800032DEST_PATH_IMAGE003
Transposing;
(2) solving the matrix
Figure 962023DEST_PATH_IMAGE007
Characteristic value of
Figure 210470DEST_PATH_IMAGE008
And feature vectors
Figure 494821DEST_PATH_IMAGE009
(ii) a Selecting a suitable cutoff value
Figure 845031DEST_PATH_IMAGE010
So that it is at the front
Figure 561445DEST_PATH_IMAGE010
The ratio of kinetic energy of system contained in each characteristic value
Figure 680711DEST_PATH_IMAGE011
Greater than 99%, here
Figure 57335DEST_PATH_IMAGE011
Expressed as:
Figure 580327DEST_PATH_IMAGE012
(3) the low-order parametric variations subspace can be described as
Figure 84120DEST_PATH_IMAGE010
Linear combination of individual feature vectors and their coefficients:
Figure 572739DEST_PATH_IMAGE013
here, the
Figure 870997DEST_PATH_IMAGE014
Is shown as
Figure 946531DEST_PATH_IMAGE015
The steady state parameter distribution corresponding to the group control variable,
Figure 254016DEST_PATH_IMAGE016
representing POD mode coefficients.
5. The method for optimizing the indoor environment of the building based on the model order reduction technology as claimed in claim 1, wherein: in the step (4), the indoor environment index comprises a thermal comfort index and an indoor air quality index, and the energy consumption index refers to an air conditioner energy consumption index, wherein:
the thermal comfort index adopts a prediction average voting index PMV;
the indoor air quality index adopts ventilation efficiency index
Here, ,and
Figure 658748DEST_PATH_IMAGE020
respectively the pollutant concentrations of the air return inlet and the air supply outlet of the air conditioner,
Figure 669429DEST_PATH_IMAGE021
average concentration of contaminants at head height of indoor personnel;
the energy consumption index of the air conditioner is decomposed into two parts, namely a fan energy consumption index and a refrigeration energy consumption index:
Figure 637385DEST_PATH_IMAGE022
here, ,the energy consumption of the fan is reduced,
Figure 382804DEST_PATH_IMAGE024
in order to increase the pressure of the fan,the total air supply quantity is the total air supply quantity,
Figure 75265DEST_PATH_IMAGE026
in order to reduce the energy consumption for refrigeration,
Figure 596376DEST_PATH_IMAGE027
for the purpose of removing the refrigeration energy consumption of the sensible heat load,
Figure 846092DEST_PATH_IMAGE028
energy consumption for dehumidifying and cooling the fresh air.
6. The method for optimizing the indoor environment of the building based on the model order reduction technology as claimed in claim 1, wherein: in the step (5), the objective function is set as:
Figure 18316DEST_PATH_IMAGE029
wherein,
Figure 328075DEST_PATH_IMAGE030
for weighting coefficients, subscripts, of the indices
Figure 336482DEST_PATH_IMAGE031
The maximum value of the corresponding performance index is used for normalization of each index;
Figure 389889DEST_PATH_IMAGE032
the punishment item is used for reflecting the influence of the head-foot temperature difference of indoor personnel and overlarge ambient wind speed on the thermal comfort; the optimization algorithm adopts a genetic algorithm, and the multi-dimensional interpolation algorithm in the parameter subspace adopts spline interpolation; in each optimization iteration process, system response is quickly obtained through multi-dimensional interpolation in the parameter subspace, and an objective function is quickly solved.
CN201210564269.3A2012-12-242012-12-24Building indoor environment optimization method based on model order reduction technologyExpired - Fee RelatedCN103049612B (en)

Priority Applications (1)

Application NumberPriority DateFiling DateTitle
CN201210564269.3ACN103049612B (en)2012-12-242012-12-24Building indoor environment optimization method based on model order reduction technology

Applications Claiming Priority (1)

Application NumberPriority DateFiling DateTitle
CN201210564269.3ACN103049612B (en)2012-12-242012-12-24Building indoor environment optimization method based on model order reduction technology

Publications (2)

Publication NumberPublication Date
CN103049612A CN103049612A (en)2013-04-17
CN103049612Btrue CN103049612B (en)2014-06-25

Family

ID=48062249

Family Applications (1)

Application NumberTitlePriority DateFiling Date
CN201210564269.3AExpired - Fee RelatedCN103049612B (en)2012-12-242012-12-24Building indoor environment optimization method based on model order reduction technology

Country Status (1)

CountryLink
CN (1)CN103049612B (en)

Families Citing this family (22)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN103995548B (en)*2014-05-132016-06-15江苏大学Indoor thermal environment control method based on model reduction and multiple model predictive control
CN104572246B (en)*2014-12-082017-11-17江苏大学For improving the interactive optimization method of environment in building
TWI524293B (en)*2014-12-102016-03-01財團法人工業技術研究院Environment control method and system thereof
US10145576B2 (en)*2015-03-062018-12-04Mitsubishi Electric Research Laboratories, Inc.System and method for controlling operations of air-conditioning system
CN108302735B (en)*2016-09-302022-10-28日本电气株式会社Apparatus and method for controlling air purification system
CN106599380B (en)*2016-11-152020-03-31天津大学Method for predicting volatile organic pollutants in indoor air of wading public place
CN106792749B (en)*2016-12-272019-12-06重庆大学wireless sensor network node deployment method based on CFD and clustering algorithm
CN107562996A (en)*2017-07-312018-01-09江苏大学A kind of multiple target integrated optimization method of Venlo types greenhouse environment parameter
CN107545100B (en)*2017-07-312020-05-05江苏大学 A High-Resolution Reduced-Order Modeling Method for Venlo-Type Greenhouse Environments
CN107991964B (en)*2017-12-192021-01-08超级智慧家(上海)物联网科技有限公司 An indoor environment intelligent monitoring and control method and a system for realizing the method
CN109282443B (en)*2018-09-052021-03-09广东工业大学 A multi-mode, low-energy-consumption indoor thermal regulation method
CN109341681B (en)*2018-10-192023-10-27珠海格力电器股份有限公司Method and system for determining position of equipment in building and air conditioner
CN109840307A (en)*2018-12-292019-06-04重庆大学A kind of indoor temperature and humidity field quick predict system and method
CN112577161B (en)*2019-09-302022-05-27北京国双科技有限公司Air conditioner energy consumption model training method and air conditioner system control method
CN111426034B (en)*2020-05-152021-06-15北京草木元人工环境技术有限公司Indoor air conditioner ventilation air flow organization control system
CN112149364B (en)*2020-09-042022-11-29同济大学Intelligent human residential environment airflow organization optimization method
CN114200986B (en)*2021-11-182023-03-21江苏大学 A Greenhouse Environmental Optimization Design Method Considering Crop Production Benefit and Energy Saving
CN114756921B (en)*2022-02-242024-11-12同济大学 A web-based sports venue design modeling method based on parametric logic drive
CN114637350A (en)*2022-02-282022-06-17青岛海尔空调电子有限公司 Space environment control method, device and storage medium
CN115854680A (en)*2022-07-272023-03-28格莱德(福建)生物科技有限公司 Constant temperature control method, equipment and storage medium in feed production process
CN116360266B (en)*2023-04-042025-09-05东南大学 A pig house temperature energy-saving control method based on multi-objective optimization algorithm
CN120065753B (en)*2025-04-252025-07-18中铁建设集团有限公司Glass enclosed space microclimate regulation and control method based on CFD-IL coupling model

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication numberPriority datePublication dateAssigneeTitle
CN101894183B (en)*2010-07-022012-05-02重庆星能建筑节能技术发展有限公司Analysis and design method of green building
CN101908086A (en)*2010-07-092010-12-08上海理工大学 An Analysis Method for Digital Dynamic Simulation of Building Indoor Wind Environment

Also Published As

Publication numberPublication date
CN103049612A (en)2013-04-17

Similar Documents

PublicationPublication DateTitle
CN103049612B (en)Building indoor environment optimization method based on model order reduction technology
Hou et al.Prediction and optimization of thermal comfort, IAQ and energy consumption of typical air-conditioned rooms based on a hybrid prediction model
Wu et al.A PMV-based HVAC control strategy for office rooms subjected to solar radiation
CN103995548B (en)Indoor thermal environment control method based on model reduction and multiple model predictive control
Magnier et al.Multiobjective optimization of building design using TRNSYS simulations, genetic algorithm, and Artificial Neural Network
CN104572246B (en)For improving the interactive optimization method of environment in building
Naboni et al.Thermal comfort-CFD maps for architectural interior design
CN106016620A (en)Energy saving thermal comfort control method of air conditioning system
CN113158593B (en)Intelligent design method of hospital ward ventilation system based on artificial intelligence
CN104866702B (en)A kind of method of utilization pivot decoupling computation hot-summer and cold-winter area building load
Li et al.Multi-objective optimization of HVAC system using NSPSO and Kriging algorithms—A case study
Wu et al.Multi-stage regression linear parametric models of room temperature in office buildings
Liu et al.Coupled simulation of natural ventilation and daylighting for a residential community design
You et al.Improving building facade design using integrated simulation of daylighting, thermal performance and natural ventilation
Zhu et al.Multi-sensor information fusion based control for VAV systems using thermal comfort constraints
Pichler et al.Simulation-assisted building energy performance improvement using sensible control decisions
Zhang et al.Enhancing indoor environmental quality: Personalized recommendation method for demand-oriented indoor ventilation strategy
Zhao et al.Inverse design of indoor radiant terminal using the particle swarm optimization method with topology concept
Ma’bdeh et al.Optimizing airflow in double-skin facades: Influence of vents design and cavity depth
Li et al.Numerical modeling of non-uniform indoor temperature distribution for coordinated air flow control
CN108241899A (en) Cost optimization calculation system and method for passive house construction equipment materials
Li et al.The solution and simulation of the condensation problem of the capillary network system in the children’s hospital of Shenyang in summer
Vasanthkumar et al.Adaptive damper control for HVAC systems based on human occupancy and indoor parameters: A development study
Lee et al.Development of Control Algorithm for Economizer System using Load Prediction Model
CN106500220B (en) Method for determining the surface temperature change of radiant cooling system radiant panels when water supply is turned off

Legal Events

DateCodeTitleDescription
C06Publication
PB01Publication
C10Entry into substantive examination
SE01Entry into force of request for substantive examination
C14Grant of patent or utility model
GR01Patent grant
CF01Termination of patent right due to non-payment of annual fee
CF01Termination of patent right due to non-payment of annual fee

Granted publication date:20140625

Termination date:20161224


[8]ページ先頭

©2009-2025 Movatter.jp