Heterogeneous phase change composite self-insulation block wall room thermal environment prediction method and system

Abstract

The invention discloses a non-homogeneous phase change composite self-insulation block wall room thermal environment prediction method and system, and belongs to the technical field of non-homogeneous phase change composite self-insulation blocks, and the method specifically comprises the steps: building a dynamic heat transfer module for simulating the heat transfer process of a non-homogeneous phase change composite self-insulation block; coupling the dynamic heat transfer module with each functional module in the building thermal process / energy consumption simulation software to obtain a non-homogeneous phase-change composite self-insulation block wall room simulation module; and based on the heterogeneous phase change building block room simulation module, predicting the thermal environment of the heterogeneous phase change composite self-insulation building block wall room. The dynamic heat transfer module capable of simulating and calculating the heat transfer process of the non-homogeneous phase-change composite self-heat-preservation building block is established and coupled to building heat process / energy consumption simulation software, and simulation prediction of the room heat environment of the non-homogeneous phase-change composite self-heat-preservation building block wall can be achieved; the problem that the thermal environment of a non-homogeneous phase-change composite self-heat-preservation building block wall room cannot be simulated and predicted is solved.

Description

Technical Field

[0001] The invention belongs to the technical field of non-homogeneous phase-change composite self-insulating building blocks, and in particular relates to a method and system for predicting the thermal environment of a room with a non-homogeneous phase-change composite self-insulating building block wall. Background Art

[0002] Phase change energy storage walls are an effective method that has attracted widespread attention in recent years and is applied to building envelope structures to reduce building energy consumption. They mainly absorb and store incident solar radiation heat and transferred heat caused by indoor and outdoor temperature differences in the building structure, thereby reducing indoor temperature fluctuations and reducing building energy consumption.

[0003] Phase change energy storage self-insulating blocks are made by filling phase change materials (PCM) into the holes of concrete self-insulating blocks to improve the thermal inertia of the self-insulating wall. The phase change process of PCM caused by temperature changes can reduce the amount of heat entering the indoor environment and the amount of heat lost from the indoor environment. Within the known melting temperature range, the long-lasting stability of the phase change material can effectively maintain the temperature of the indoor environment, thereby affecting thermal comfort and reducing the energy consumption of the building. At the same time, this application form is not only easy to install and construct, but the phase change material layer does not occupy the thickness of the wall, and the thermal insulation structure can be integrated. It has the advantages of good durability, fire resistance, convenient construction, and the same life as the building. It is a technology worthy of widespread application.

[0004] The optimal design of phase change self-insulating blocks usually involves the optimal selection of block geometric dimensions, hole shape and geometric dimensions, and PCM parameters. It is necessary to evaluate the application effects of phase change self-insulating blocks under different schemes. This requires predicting the thermal environment effects of phase change self-insulating block wall rooms under different schemes to guide the optimal design of phase change self-insulating blocks and the optimal design of phase change self-insulating block wall rooms.

[0005] Due to the high non-uniformity of phase-change composite self-insulating blocks, there are significant differences in thermal physical parameters such as thermal conductivity of each material, making the heat transfer process inside the composite self-insulating blocks quite complicated. The physical parameters of phase-change materials change in real time with boundary conditions, which belongs to the category of multi-dimensional heat transfer. Some existing numerical simulation software such as Ansys and Fluent can only simulate the three-dimensional heat transfer of a single block or a single wall, and cannot be coupled with the simulation of the indoor environment, and cannot evaluate the actual application effect of the block wall under study. Existing building thermal process and energy consumption simulation software such as Energyplus and TRNSYS can realize the simulation coupling of non-steady-state heat transfer of walls and thermal environment of rooms. However, this type of software treats the material layer as a homogeneous layer for the calculation modules of phase change materials and each layer of materials for multi-layer wall structures, and uses numerical methods such as finite difference method to solve and simulate. It is impossible to simulate the non-steady-state heat transfer process of non-homogeneous phase change composite self-insulating block walls, so it is impossible to predict the thermal environment of rooms with non-homogeneous phase change composite self-insulating block walls, and it is impossible to guide the optimal design of non-homogeneous phase change composite self-insulating block walls and the optimal design of rooms with non-homogeneous phase change composite self-insulating block walls.

[0006] Therefore, in view of the above problems, it is necessary to design and develop a thermal environment prediction method and system for a room with non-homogeneous phase change composite self-insulating block walls. Summary of the invention

[0007] The purpose of the present invention is to provide a method and system for predicting the thermal environment of a room with a non-uniform phase change composite self-insulating block wall in view of the problems existing in the prior art. By building a heat transfer model of a non-uniform phase change composite self-insulating block, a dynamic heat transfer module for simulating the heat transfer process of the non-uniform phase change composite self-insulating block is established, and the module is coupled with a building module, a data reading module and a simulation output module in a building thermal process / energy consumption simulation software that can simulate the heat transfer of a wall homogeneous building material layer and the indoor heat transfer process, so as to finally realize the simulation prediction of the thermal environment of a room with a non-uniform phase change composite self-insulating block wall.

[0008] According to one aspect of the present invention, a method for predicting the thermal environment of a room with a non-homogeneous phase-change composite self-insulating block wall is provided, comprising:

[0009] Establishing a heat transfer model of a non-homogeneous phase-change composite self-insulating building block, including: establishing a physical model and a solution unit of the non-homogeneous phase-change composite self-insulating building block, wherein the holes on one side of the non-homogeneous phase-change composite self-insulating building block are filled with phase-change materials, and the holes on the other side are filled with insulation materials; based on the physical model and the solution unit, establishing a control equation including a heat transfer process of the phase-change material part of the building block and a heat conduction process of the base material and insulation material area of ​​the building block; establishing boundary conditions for the heat transfer process of the non-homogeneous phase-change composite self-insulating building block; constructing a solution method and process based on the control equation and the boundary conditions, analyzing the dynamic heat transfer process of the non-homogeneous phase-change composite self-insulating building block, and obtaining the temperature field distribution of the overall space of the building block;

[0010] Based on the heat transfer model of the non-homogeneous phase change composite self-insulating building blocks, a dynamic heat transfer module of the non-homogeneous phase change composite self-insulating building blocks was built in the building thermal process / energy consumption simulation software using programming language;

[0011] The dynamic heat transfer module of the non-homogeneous phase-change composite self-insulating building block is coupled with the building module, the data reading module and the simulation output module in the building thermal process / energy consumption simulation software, which can simulate the heat transfer of the wall homogeneous material layer and the indoor heat transfer process, to obtain the non-homogeneous phase-change composite self-insulating building block wall room simulation module of the building thermal process / energy consumption simulation software;

[0012] Based on the non-homogeneous phase change block room simulation module of the building thermal process / energy consumption simulation software, the thermal environment of the room with non-homogeneous phase change composite self-insulating block walls is predicted.

[0013] Furthermore, a heat transfer model of a non-homogeneous phase change composite self-insulating building block is established, specifically including:

[0014] A physical model of a non-homogeneous phase-change composite self-insulating building block is established, wherein the holes of the non-homogeneous phase-change composite self-insulating building block close to the indoor side are filled with phase-change materials, and the holes of the outdoor side are filled with insulation materials, and the phase-change materials are arranged in a staggered manner relative to the insulation materials in the thickness direction of the building block; a representative target area is selected based on the physical model to establish a solution unit;

[0015] Based on the physical model and solution unit, the heat flow conduction from the high temperature side to the low temperature side of the enclosure structure and the lateral heat transfer between different materials are considered, and the control equation reflecting the heat transfer process of the heterogeneous phase change composite self-insulating block is established, which is as follows:

[0016] According to the law of conservation of energy, the calculation process of specific heat capacity and latent heat of phase change is integrated into one item, In order to comprehensively calculate the heat capacity, the heat transfer process of the phase change material part of the building block is described by the following governing equation:

[0017]

[0018] Among them, x is the length direction of the block, y is the thickness direction of the block, is the thermal conductivity of the phase change material, is the temperature, t is the time variable, is the density of the phase change material, L is the latent heat of phase change, is the specific heat capacity of the phase change material, T s and T l Respectively represent the phase change starting temperature and phase change ending temperature of the phase change material;

[0019] For the base material and insulation material area of ​​the block, the heat transfer process is a pure heat conduction process. The material properties do not consider the real-time changes of the non-steady-state heat transfer process. The heat transfer process is described by the following control equation:

[0020]

[0021] Among them, x is the length direction of the block, y is the thickness direction of the block, is the thermal conductivity of the substrate / insulation material, is the temperature, t is the time variable, is the specific heat capacity of the substrate / insulation material, is the density of the substrate / insulation material;

[0022] The boundary conditions of the heat transfer process of the non-homogeneous phase change composite self-insulating building blocks are established, including: based on the heat transfer process of the building blocks, the boundary conditions on both sides of the length direction x of the building blocks are described by the differential formula of the adiabatic boundary conditions; for the boundary conditions on both sides of the thickness direction y of the building blocks, the differential formula of the constant heat flux boundary conditions that reflects the correlation between the heat flux density on the outdoor and indoor surfaces of the building blocks and the thermal conductivity of the substrate is used to describe.

[0023] Furthermore, a solution method and process are constructed based on the control equations and boundary conditions, including:

[0024] The heat capacity method is used to process the specific heat capacity of phase change materials under different temperature conditions. The calculation formula is as follows:

[0025]

[0026] The linear interpolation method is used to obtain the density and thermal conductivity of the phase change material at different temperatures. The calculation formula is as follows:

[0027]

[0028]

[0029] The calculation formula of the liquid phase fraction f is shown as follows. When f = 0, the phase change material is in a pure solid state; when f = 1, the phase change material is in a pure liquid state; when 0 < f < 1, the phase change material is in a mixed state.

[0030]

[0031] In the formula, f is the liquid phase fraction; T is the real-time temperature of the phase change material, °C; T s and T l represent the phase change start temperature and the phase change end temperature of the phase change material, respectively.

[0032] Furthermore, in the process of analyzing the dynamic heat transfer of the heterogeneous phase change composite self-insulating block, the control equation and boundary conditions of the heat transfer model are discretized in the spatial dimension, including:

[0033] The finite volume method is used to perform spatial division of the control volume for the solution unit of the heterogeneous phase change composite self-insulating block. The center of each control volume of the solution unit of the heterogeneous phase change composite self-insulating block is denoted as the calculation node of each control volume; the control equation is converted into a discrete equation at each node, and the difference of the spatial derivative is calculated based on the temperature value at the end of each time layer, and time discretization in implicit format is performed to obtain the temperature of each control volume, and then the temperature field distribution of the overall space of the block is obtained.

[0034] Furthermore, a dynamic heat transfer module for the heterogeneous phase change self-insulating block is built in the building heat process / energy consumption simulation software, including:

[0035] A variable table is newly defined in the building heat process / energy consumption simulation software, including input items, output items and parameter items;

[0036] Combined with the module definition variable table, the heat transfer model of the heterogeneous phase change composite self-insulating block and its dynamic heat transfer analysis process are written as program source code using a programming language;

[0037] The building heat process / energy consumption simulation software is used to compile the program source code into a directly loadable dynamic link library file (Dynamic Link Library, DLL);

[0038] Based on the directly loadable dynamic link library file, a dynamic heat transfer module for the heterogeneous phase change self-insulating block is built.

[0039] Furthermore, the dynamic heat transfer module of the heterogeneous phase change composite self-insulating block is coupled with the building module, data reading module and simulation output module in the building heat process / energy consumption simulation software that can simulate the heat transfer of the homogeneous material layer of the wall and the indoor heat transfer process, specifically including:

[0040] The dynamic heat transfer module of the non-homogeneous phase change composite self-insulating block uses the calculated average temperature of the outdoor and indoor surfaces of the composite self-insulating block as output data, and outputs it to the building module as the boundary temperature input data of the plastering mortar layer on both sides of the outdoor and indoor sides of the composite self-insulating block and the adjacent interface of the block; for the plastering mortar layer on both sides of the block, the building module is used to simulate its heat transfer process and output the surface heat flux density of the adjacent interface between the plastering mortar layer and the block as the input heat flux density of the inner and outer surfaces of the non-homogeneous phase change composite self-insulating block, which is transmitted to the surface heat flux density input end of the dynamic heat transfer module respectively.

[0041] Furthermore, based on the non-homogeneous phase change block room simulation module of the building thermal process / energy consumption simulation software, the thermal environment of the non-homogeneous phase change composite self-insulating block wall room is predicted, including:

[0042] The building module is used to build a building model, and the parameters of the internal space structure and enclosure structure as well as the heat dissipation of internal personnel, equipment and lighting are set. It is also used to simulate the dynamic heat transfer process of the cement mortar layer on both sides of the phase change self-insulating block wall and the heat transfer process of the internal space; the data reading module is used to read the weather data outside the building; the dynamic heat transfer module is used to simulate the dynamic heat transfer process of the non-homogeneous phase change composite self-insulating block; the simulation output module is used to output the simulated room thermal environment prediction results as data files and result displays.

[0043] According to one aspect of the present invention, a system for predicting the thermal environment of a room with a non-homogeneous phase-change composite self-insulating block wall is provided, comprising:

[0044] Establish a model module and establish a heat transfer model of a non-homogeneous phase change composite self-insulating building block, including: establishing a physical model and a solution unit of the non-homogeneous phase change composite self-insulating building block; based on the physical model and the solution unit, establishing a control equation including the heat transfer process of the phase change material part of the building block and the heat conduction process of the base material and insulation material area of ​​the building block; establishing boundary conditions for the heat transfer process of the non-homogeneous phase change composite self-insulating building block; constructing a solution method and process based on the control equation and boundary conditions to analyze the dynamic heat transfer process of the non-homogeneous phase change composite self-insulating building block;

[0045] Self-built dynamic heat transfer module, based on the heat transfer model of heterogeneous phase change composite self-insulating building blocks, using programming language to build the dynamic heat transfer module of heterogeneous phase change composite self-insulating building blocks in the building thermal process / energy consumption simulation software;

[0046] A simulation module is used to couple the heat transfer module of the non-homogeneous phase-change composite self-insulating building block with the building module, the data reading module and the simulation output module in the building thermal process / energy consumption simulation software that can simulate the heat transfer of the wall homogeneous material layer and the indoor heat transfer process, so as to obtain the non-homogeneous phase-change composite self-insulating building block wall room simulation module of the building thermal process / energy consumption simulation software;

[0047] The room thermal environment prediction module is based on the non-homogeneous phase change block room simulation module of the building thermal process / energy consumption simulation software to predict the thermal environment of the room with non-homogeneous phase change composite self-insulating block walls.

[0048] According to one aspect of the present invention, there is provided an electronic device, comprising a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the steps of the method for predicting the thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall are implemented.

[0049] According to one aspect of the present invention, there is provided a computer-readable storage medium storing a computer program, which, when executed by a processor, implements the steps of the method for predicting the thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall.

[0050] Compared with the prior art, the present invention has the following beneficial effects:

[0051] 1. The present invention proposes a method and system for predicting the thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall. By establishing a heat transfer model of a non-homogeneous phase change composite self-insulating block with hollow holes filled with phase change material and thermal insulation material, the model comprehensively considers the effects of multi-porous and multi-material (concrete, thermal insulation material, phase change material) composite blocks, the heterogeneity of the building envelope structure layer, and the solid / liquid phase change characteristics of the phase change material on heat transfer, and can more accurately simulate the heat transfer process of the non-homogeneous phase change block.

[0052] 2. The present invention proposes a method and system for predicting the thermal environment of a room with a non-uniform phase change composite self-insulating block wall. By coupling the heat transfer model of the non-uniform phase change composite self-insulating block with the building thermal process / energy consumption simulation software, a dynamic heat transfer module for simulating and calculating the heat transfer process of the non-uniform phase change composite self-insulating block is developed. This solves the problem that the existing numerical simulation software and the existing building thermal process and energy consumption simulation software cannot realize the thermal environment of a room with a non-uniform phase change composite self-insulating block wall. The method and system can simulate the indoor environmental temperature indicators of the room when non-uniform phase change composite self-insulating block walls are used for walls of different room types and different orientations, including the indoor air temperature and the inner surface temperature of the room block wall.

[0053] 3. The present invention proposes a method and system for predicting the thermal environment of a room with a non-uniform phase change composite self-insulating block wall. By coupling the heat transfer model of the non-uniform phase change composite self-insulating block with the building thermal process / energy consumption simulation software, a dynamic heat transfer module for simulating and calculating the heat transfer process of the non-uniform phase change composite self-insulating block is developed. The simulation of the heat transfer process of the non-uniform phase change composite self-insulating block wall is coupled with the simulation of the indoor thermal environment of the building. The method can simulate actual application scenarios and more accurately simulate the actual dynamic heat transfer process of the non-uniform phase change blocks. The heat transfer heat flux of the non-uniform phase change composite self-insulating block can be predicted, and the actual application effect of the non-uniform phase change composite self-insulating block wall can be evaluated.

[0054] 4. The present invention proposes a method and system for predicting the thermal environment of a room with a non-uniform phase change composite self-insulating block wall. By predicting the wall and room temperatures of non-uniform phase change composite self-insulating blocks under different PCM design schemes, different insulation material design schemes, and different substrate design schemes, the design of non-uniform phase change composite self-insulating blocks and block walls can be optimized. BRIEF DESCRIPTION OF THE DRAWINGS

[0055] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0056] Figure 1 A flow chart of a method for predicting the thermal environment of a room with a non-homogeneous phase-change composite self-insulating block wall provided by an embodiment of the present invention;

[0057] Figure 2 A schematic diagram of a non-homogeneous phase-change composite self-insulating building block (a) and a wall model (b) provided in an embodiment of the present invention;

[0058] Figure 3 A physical model of a non-homogeneous phase-change composite self-insulating building block and a schematic diagram of a solution unit structure provided in an embodiment of the present invention;

[0059] Figure 4 A schematic diagram of discrete nodes of a solution unit provided in an embodiment of the present invention;

[0060] Figure 5 An interface view of a dynamic heat transfer module provided by an embodiment of the present invention;

[0061] Figure 6 A schematic diagram of the modules involved in the heat transfer process of a non-homogeneous phase-change composite self-insulating block wall provided by an embodiment of the present invention;

[0062] Figure 7 A schematic diagram of the connection between the Type56 module and the Type451 module provided in an embodiment of the present invention;

[0063] Figure 8 A schematic diagram of a building model for simulation prediction provided by an embodiment of the present invention;

[0064] Fig. 9 A schematic diagram of a block wall structure provided by a simulation prediction according to an embodiment of the present invention;

[0065] Fig.10 A schematic diagram of the process of building a TRNSYS software simulation interface provided by an embodiment of the present invention;

[0066] Fig.11 A curve diagram of the temperature variation of the inner surface of a block wall at different phase change temperatures provided by an embodiment of the present invention;

[0067] Fig.12 A curve diagram of room air temperature variation at different phase change temperatures provided by an embodiment of the present invention;

[0068] Fig.13 A module diagram of a room thermal environment prediction system for a non-homogeneous phase change composite self-insulating block wall provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0069] It should be noted that:

[0070] The Type451 module is a dynamic heat transfer module of a self-built non-homogeneous phase-change composite self-insulating building block;

[0071] Type56 module is a standard module in TRNSYS software, namely building module, which can simulate the heat transfer of homogeneous material layer of wall and indoor heat transfer process;

[0072] The Type9e module is a standard module in the TRNSYS software, i.e., a data reading module;

[0073] Type65a module is a standard module in TRNSYS software, that is, simulation output module.

[0074] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0075] like Figure 1As shown, the present invention provides a method for predicting the thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall, comprising:

[0076] Establish a heat transfer model of a non-homogeneous phase-change composite self-insulating building block: establish a physical model and a solution unit of a non-homogeneous phase-change composite self-insulating building block, wherein the holes on one side of the non-homogeneous phase-change composite self-insulating building block are filled with phase-change materials, and the holes on the other side are filled with insulation materials; based on the physical model and the solution unit, establish control equations including the heat transfer process of the phase-change material part of the building block and the heat conduction process of the base material and insulation material area of ​​the building block; establish boundary conditions for the heat transfer process of the non-homogeneous phase-change composite self-insulating building block; construct a solution method and process based on the control equations and boundary conditions, analyze the dynamic heat transfer process of the non-homogeneous phase-change composite self-insulating building block, and obtain the temperature field distribution of the overall space of the building block;

[0077] Based on the heat transfer model of non-homogeneous phase change composite self-insulating building blocks, a dynamic heat transfer module of non-homogeneous phase change composite self-insulating building blocks was built in the building thermal process / energy consumption simulation software using programming language;

[0078] The dynamic heat transfer module of the non-homogeneous phase-change composite self-insulating building block is coupled with the building module, the data reading module and the simulation output module in the building thermal process / energy consumption simulation software, which can simulate the heat transfer of the wall homogeneous material layer and the indoor heat transfer process, to obtain the non-homogeneous phase-change composite self-insulating building block wall room simulation module of the building thermal process / energy consumption simulation software;

[0079] Based on the non-homogeneous phase change block room simulation module of the building thermal process / energy consumption simulation software, the thermal environment of the room with non-homogeneous phase change composite self-insulating block walls is predicted.

[0080] It should be noted here that the model constructed by the present invention is suitable for thermal environment simulation of various types of non-homogeneous phase change composite self-insulating block walls. It is just that the physical model of the non-homogeneous phase change composite self-insulating block is different. The following specific embodiments only belong to part of the embodiments of the present invention and are used to more deeply understand the technical concept of the present invention.

[0081] In an embodiment of the present invention, a Type451 module for simulating the dynamic heat transfer process of non-homogeneous phase change composite self-insulating blocks is built in the TRNSYS software, and the Type451 module is coupled with the Type56 module, Type9e module and Type65a module in the TRNSYS software to obtain the non-homogeneous phase change composite self-insulating block wall room simulation module of the TRNSYS software; based on the non-homogeneous phase change composite self-insulating block wall room simulation module of the TRNSYS software, the thermal environment of the non-homogeneous phase change composite self-insulating block wall room is predicted.

[0082] In a specific embodiment, a physical model and a solution unit of a non-homogeneous phase-change composite self-insulating building block are established, including:

[0083] Firstly, the composite self-insulating building block containing PCM was modeled. The physical model size of the building block is 390mm×280mm×190mm. Figure 2 As shown in (a), the block type of non-homogeneous phase change composite self-insulating building block is demonstrated. Figure 2 (b) Figure 3 The physical model and solution unit of the block wall model and the non-homogeneous phase change composite self-insulating block are shown respectively. PCM is filled in the rectangular cavity of the block, close to the indoor side, e.g. Figure 2 (b) The green part shows that the thermal insulation material is filled in the rectangular cavity outside the room, as shown in Figure 2 (b) As shown in the white part, it can be seen from the figure that the two rows of blocks on the outdoor side are staggered and filled with insulation materials (it should be noted that the insulation materials in the actual model construction can be set in one row or multiple rows, and both can be simulated by the method of the present invention). Since the block itself is a symmetrical object, when modeling, it is sufficient to take 1 / 2 of the block to establish a mathematical solution unit, that is, a symmetrical solution, such as Figure 3 As shown. The physical model of the heterogeneous phase change composite self-insulating building block established contains three constituent materials, and the same plane is composed of a base material, a thermal insulation material, and a phase change material. The solution unit and the establishment method established by the present invention can not only reflect the multi-dimensional heat transfer process of the heterogeneous building block, but also simplify the simulation calculation amount.

[0084] The mathematical model of heat transfer of non-homogeneous phase change composite self-insulating blocks is established, including the control equations of the heat transfer process of the phase change material part of the block and the heat conduction process of the base material and insulation material area of ​​the block, as follows:

[0085] The heat transfer of the non-homogeneous phase-change composite self-insulating blocks belongs to a typical three-dimensional heat transfer process, and has heat transfer in the three dimensions of length, height and thickness. However, when the wall is neatly arranged by the same non-homogeneous phase-change composite self-insulating blocks, the temperature distribution of any section (x, z plane) of the blocks in the height direction (z-axis direction) of the wall is the same, so the heat transfer in the height direction of the wall can be ignored in the actual calculation, and its heat flow is considered to be zero. Therefore, the heat transfer model developed by the present invention only starts from a two-dimensional perspective (x, y direction). The y direction represents the thickness direction of the blocks, and the x direction represents the length direction of the blocks. The two-dimensional heat transfer analysis method simultaneously considers the conduction of heat flow from the high-temperature side to the low-temperature side of the enclosure structure and the "lateral heat transfer" between different materials, which can accurately calculate the non-steady-state heat transfer of the non-homogeneous block wall, and can more accurately reflect the actual heat transfer process of the non-homogeneous phase-change composite self-insulating blocks. In addition, the use of the two-dimensional heat transfer analysis method can also reduce the amount of simulation calculations. The heat transfer model of the inhomogeneous phase change composite self-insulating block is established by considering the heat transfer process between the PCM layer, the insulation material and the substrate, and the heat conduction process between the composite block and other layers of the wall. On the basis of ensuring the calculation accuracy, the inhomogeneous phase change composite self-insulating block is moderately simplified, and the following assumptions are made:

[0086] (1) Assume that the filling materials of each hole of the block are closely connected to the substrate and there is no contact thermal resistance;

[0087] (2) It is assumed that each layer of material constituting the block is uniform, consistent, and has isotropic characteristics;

[0088] (3) When the phase change material melts into liquid, its natural convection effect is not considered;

[0089] (4) In the "solid-liquid coexistence" state, the thermal parameters (thermal conductivity, density, specific heat, etc.) of the phase change material change with temperature. In a single solid / liquid state, the thermophysical parameters of the PCM are constant, that is, they are independent of temperature;

[0090] (5) Phase change materials are isotropic during the phase change process and do not exchange mass with the outside world;

[0091] Based on the above assumptions, the control equation of heat transfer of the non-homogeneous phase change composite self-insulating building block is established, and the heat transfer process of the phase change material part of the building block is described by the following equation:

[0092] (1)

[0093] Among them, x is the length direction of the block, y is the thickness direction of the block, is the thermal conductivity of the phase change material, is the temperature, t is the time variable, is the density of the phase change material, L is the latent heat of phase change, $c_p$ is the specific heat capacity of the phase change material, and $f$ is the liquid fraction.

[0094] The liquid fraction $f$ represents the proportion of the liquid substance in the mixed state of the phase change material. The thermal conductivity and density of the PCM are both affected by the liquid fraction $f$. The calculation formula for the liquid fraction $f$ is shown in Equation (2). When $f = 0$, the PCM is in a pure solid state; when $f = 1$, the PCM is in a pure liquid state; when $0 < f < 1$, the PCM is in a mixed state.

[0095] (2)

[0096] In the formula, $f$ is the liquid fraction; $T$ is the real-time temperature of the phase change material, in °C; $T$ s and $T$ l represent the phase change start temperature and the phase change end temperature of the phase change material, respectively. By substituting Equation (2) into Equation (1) and arranging, we get:

[0097] (3)

[0098] In Equation (3), the calculation processes of the specific heat capacity and the latent heat of phase change are integrated into one term, and the concept of the sensible heat capacity $C$ p,e of the phase change material is introduced. Let be the comprehensive heat capacity, then Equation (3) is arranged as:

[0099] (4)

[0100] The heat capacity method is used to handle the problem that the specific heat capacity of the phase change material is different under different temperature conditions, as shown in Equation (5):

[0101] (5)

[0102] In the formula, $C$ s and $C$ l represent the specific heat capacities of the phase change material in the solid state and the liquid state, respectively.

[0103] The linear interpolation method is used to obtain the density and thermal conductivity of the PCM at different temperatures, as shown in Equation (6) and Equation (7);

[0104] (6)

[0105] (7)

[0106] Among them, 、 are the densities of the PCM in the pure solid state and the pure liquid state, respectively; 、 are the thermal conductivities of the PCM in the pure solid state and the pure liquid state, respectively;

[0107] For the base material and insulation material area of ​​the block, the heat transfer process is a heat conduction process. The material properties do not consider the real-time changes of the non-steady-state heat transfer process. The heat transfer process is described by the following equation:

[0108] (8)

[0109] Where x is the length direction of the block, y is the thickness direction of the block, is the thermal conductivity of the substrate / insulation material, is the temperature, t is the time variable, is the specific heat capacity of the substrate / insulation material, is the density of the substrate / insulation material;

[0110] Equations (1) to (8) constitute the control equations of the heat transfer process of the non-homogeneous phase change composite self-insulating building block.

[0111] Models and solution units based on physical modeling, such as Figure 3 As shown in the figure, the boundary conditions for heat transfer of non-homogeneous phase change composite self-insulating blocks are established. Figure 3 As shown, due to the symmetry of the block model, the units on both sides of the solution unit are consistent with the unit structure, and the heat transfer process of the same structure is the same, and there is no lateral heat transfer process between the connected units. The symmetry surfaces adfj and kpqt in the solution unit are set as adiabatic boundary conditions, and the boundary surfaces ak and jt on both sides of the block thickness direction of the solution unit, i.e., the y direction, are set as constant heat flow boundaries. The boundary conditions are described by the following equations:

[0112] , , , (9)

[0113] in, is the thermal conductivity of the non-homogeneous phase change composite self-insulating block concrete substrate, is the heat flux density on the indoor surface of the inhomogeneous phase change composite self-insulating block layer, is the heat flux density on the outdoor surface of the composite self-insulating block layer through the inhomogeneous phase change.

[0114] Specifically, the heat transfer model of the inhomogeneous phase change composite self-insulating block is solved, and the dynamic heat transfer process of the inhomogeneous phase change composite self-insulating block is analyzed, including: using the numerical solution method to solve the heat transfer partial differential equation (control equation) in the heat transfer process of the inhomogeneous phase change composite self-insulating block, the functional relationship between temperature (T) and spatial variables (x, y, z) and time variables (t) can be solved to obtain the temperature field distribution inside the block. The numerical solution of phase change heat transfer adopts the fixed grid method. In the fixed grid method, the heat transfer model is established by using the energy conservation in the entire phase change region, and the calculation region can be regarded as a series of interconnected grids. The heat transfer process of the entire phase change material can be constructed by calculating the heat change in each grid, and this method does not require the determination of the solid-liquid interface. When numerically analyzing the heat transfer model of the inhomogeneous phase change composite self-insulating block, it is necessary to discretize the control equation and boundary conditions of the heat transfer process of the inhomogeneous phase change composite self-insulating block in the spatial dimension.

[0115] Specifically, the finite volume method is used to perform spatial division of the control volume of the non-homogeneous phase change composite self-insulating building block solution unit, and the solution unit is divided into multiple control volumes of equal volume. The center point of each control volume is used as the calculation node to discretize and obtain the discrete model of the non-homogeneous phase change composite self-insulating building block monomer. The discretized node configuration is shown in Figure 4 Among them, the total number of control bodies in the solution unit of the non-homogeneous phase change composite self-insulating building block is , represents the number of control bodies along the x-axis in the solution unit of the non-homogeneous phase change composite self-insulating building block, It represents the number of control bodies along the y-axis direction in the solution unit of the non-homogeneous phase change composite self-insulating building block, and the number of lateral control bodies in the discrete model of the non-homogeneous phase change composite self-insulating building block solution unit: , the number of longitudinal control bodies: The side length of the control volume in the solution unit of the non-homogeneous phase change composite self-insulating building block is recorded as ; The center of each control body in the non-homogeneous phase change composite self-insulating building block solution unit is recorded as the node of each control body, and the The nodes of the control volume are denoted as nodes , the nodes of the four adjacent control bodies above, below, left and right of the node are respectively recorded as , , and ;in, Refers to the column number of the control body node in the X-axis direction. Represents the row number of the control body node in the Y-axis direction; The value is 0.1 mm, and the node center is indicated by a superscript.

[0116] When dealing with the non-steady-state (dynamic) heat transfer problem of the heterogeneous phase-change composite self-insulating building block of the present invention, it is necessary to clarify how to calculate the difference of the spatial derivative at each time layer, starting from the initial time layer and proceeding step by step. The embodiment of the present invention selects an implicit format for time discretization, and calculates the difference of the spatial derivative based on the temperature value at the end of each time layer. This discretization method has better stability and is not limited by the time step and the space step. When solving the entire area of ​​the model, the nodes are divided into internal nodes and boundary nodes (including upper, lower, left, and right boundary nodes and double boundary nodes at the four corners). Based on the division and discretization of the solution area control volume, the partial differential control equation is converted into a discrete equation (temperature solution model) at each node, as described below:

[0117] 1. When = When in the upper left corner of the block solution unit, use the computer to create the node The temperature solution model of the nodes of the two adjacent control bodies on the right and below is as follows:

[0118] (10)

[0119] in, Indicates that the node is Density of Moments The heat capacity method is used to obtain the Specific heat capacity at time The product of Representation Node In the The temperature of the moment; Representation Node In the The temperature of the moment; Representation Node In the The temperature of the moment; Representation Node In the The temperature of the moment; Representation Node In the The temperature of the moment; Representation Node In the The temperature of the moment; Indicates Time and the time step between moments; = ; Representation Node and nodes Thermal conductivity between Representation Node and nodes Thermal conductivity between Representation Node and nodes Thermal conductivity between Representation Node and nodes Thermal conductivity between Representation Node Thermal conductivity of Representation Node Thermal conductivity of Representation Node Thermal conductivity between Representation Node Thermal conductivity.

[0120] 2. When = When the node is at the upper boundary of the block solution unit and not at the four corners, the computer is used to establish the node The temperature solution model of the nodes of the three adjacent control bodies on the left and right is as follows:

[0121] (11)

[0122] in, It is the physical property parameter of concrete substrate. is the heat flux density on the indoor surface of the composite self-insulating block layer through the inhomogeneous phase change;

[0123] 3. When = When in the upper right corner of the block solution unit, use the computer to create a node The temperature solution model of the nodes of the two adjacent control bodies on the left and below is as follows:

[0124] (12)

[0125] in, It is the physical property parameter of concrete substrate. is the heat flux density on the indoor surface of the composite self-insulating block layer through the inhomogeneous phase change;

[0126] 4. When = , = When it is at the left boundary of the block solution unit and not at the four corners, the node is established by computer The temperature solution model of the nodes of the three adjacent control bodies is as follows:

[0127] (13)

[0128] in, It is the physical property parameter of concrete substrate;

[0129] when = When, in formula (13) It is the physical parameter of the thermal insulation material;

[0130] 5. When = , = When it is at the right boundary of the block solution unit and not at the four corners, the node is established by computer The temperature solution model of the nodes of the three adjacent control bodies above, below and to the left is as follows:

[0131] (14)

[0132] in, It is the physical property parameter of concrete substrate;

[0133] when = When, (14) It is the physical parameter of thermal insulation material.

[0134] 6. When When not at the boundary of the block solution unit, the node is established by computer The temperature solution model of the nodes of the four adjacent control bodies is as follows:

[0135] (15)

[0136]

[0137] in, are the physical properties parameters of each material.

[0138] The node is the concrete substrate part.

[0139] The node is the insulation material part.

[0140] Part of the nodes for phase change materials.

[0141] It should be noted that: When it is inside the block solution unit and is the boundary of the insulation material or phase change material, where:

[0142] A.When When it is at the upper boundary of the insulation material, :

[0143] B. When At the lower boundary of the insulation material, , : ,

[0144] C. When When it is at the left boundary of the insulation material, : :

[0145] D. When When it is at the right boundary of the insulation material, , :

[0146] E.When At the upper boundary of the phase change material, :

[0147] F. When At the lower boundary of the phase change material, : ,

[0148] G. When At the left boundary of the phase change material, : :

[0149] H. When At the right boundary of the phase change material, :

[0150] 7. When = When in the lower left corner of the block solution unit, use the computer to create the node The temperature solution model of the nodes of the two adjacent control bodies on the right is as follows:

[0151] (16)

[0152] in, It is the physical property parameter of concrete substrate. is the heat flux density on the outdoor surface of the composite self-insulating block layer through the inhomogeneous phase change.

[0153] 8. When = When the node is at the lower boundary of the block solution unit and not at the four corners, the computer is used to establish the node The temperature solution model of the nodes of the three adjacent control bodies on the left and right is as follows:

[0154] (17)

[0155] in, It is the physical property parameter of concrete substrate. is the heat flux density on the outdoor surface of the composite self-insulating block layer through the inhomogeneous phase change.

[0156] 9. When = When in the lower right corner of the block solution unit, use the computer to create the node The temperature solution model of the nodes of the two adjacent control bodies on the upper left is as follows:

[0157] (18)

[0158] in, It is the physical property parameter of concrete substrate. is the heat flux density on the outdoor surface of the composite self-insulating block layer through the inhomogeneous phase change.

[0159] In some embodiments, a method for using TRNSYS software to self-build a Type451 module of a dynamic heat transfer module for a non-homogeneous phase change composite self-insulating building block includes: opening Trnsys Simulation Studio in TRNSYS to generate a new module definition variable table, input items INPUTS (q1, q2), output items OUTPUTS (surface temperatures on both sides of the indoor and outdoor sides of the building block, surface temperatures and average temperatures of both sides of the indoor and outdoor directions of the PCM, liquid phase rate, etc.) and parameter items PARAMETERS (physical properties of phase change materials, physical properties of concrete substrates and insulation materials, time step, etc.); combining the variable table, using FORTRAN language, the non-steady-state heat transfer numerical calculation process of the non-homogeneous phase change composite self-insulating building block is written as TRNSYS source code; using the TRNSYS built-in program TypeStudio to compile it into a DLL file that can be directly loaded by the TRNSYS kernel, and then self-building the Type451 module. Figure 5Table 1 is the interface view of the Type451 module, and Table 2 is the variable table of the Type451 module, including input items, output items, and parameter items, where the output items are calculation results.

[0160] Table 1 Variable table of Type451 module

[0161]

[0162] It should be pointed out that the simulation of the wall composed of non-homogeneous phase change composite self-insulating blocks only needs to determine the wall area according to the room depth and floor height, and enter the specific area size in the "Area" option in the Type451 module, and the block wall model will be automatically generated. Even if the wall area is not an integer multiple of the block, the heat transfer principle and modeling ideas are consistent, and it will not affect the calculation accuracy of the heat transfer of the non-homogeneous phase change composite self-insulating block wall.

[0163] The embodiment of the present invention also provides for connecting other modules in TRNSYS to the input and output ports of Type451 to achieve coupling. In the simulation of the internal heat transfer process of the composite self-insulating wall, Type451 is responsible for realizing the heat transfer simulation of the phase change blocks, while the heat transfer process of the phase change block wall plastering mortar layer or other homogeneous material layer and other homogeneous enclosure structures in the room (ceiling, floor, etc.) as well as the internal disturbance of personnel, equipment, lighting, and the addition of air-conditioning systems are completed by Type56. By coupling these two modules, the coupled simulation of the heat transfer of the non-homogeneous phase change composite self-insulating blocks and the plastering mortar layers on both sides is realized, thereby realizing the heat transfer simulation of the entire block wall and the simulation of the thermal environment of the room. The following specifically describes how to link and interact these two modules:

[0164] Building geometry, wall structure, and material properties are defined in TRNbuild Type56, while outdoor air temperature and outdoor relative humidity are associated with buildings using Type56 through existing data in the database. Type451 module simulates the surface temperatures "Tsb" and "Txb" on both sides of the indoor and outdoor surfaces of the blocks as output data, and the custom "Tpcm1" and "Tpcm2" are used as inputs of Type56 to output the temperatures to the indoor gypsum mortar layer and the outdoor cement mortar layer, respectively. Type56 simulates the heat transfer process inside the building room and the heat transfer of the mortar layers on both sides of the block wall. For the mortar layers on both sides, two "QCOMO" outputs (NTYPE20) generated by Type56 are used. These two outputs provide the heat flux density through the boundary surface between each mortar layer and the block, which are used as the input heat flux density of the indoor and outdoor surfaces of the non-homogeneous phase change composite self-insulating block heat transfer model, and are respectively transferred to the "q1" and "q2" input terminals of Type451. Figure 6This is a schematic diagram of the use of Trnsys modules for different components of a block wall. The mortar layer is simulated by Type56 and the block part is simulated by Type451. Figure 7 For the connection process of "Type56" and "Type451", the interactive module can perform multiple calculations by adjusting various parameters according to actual needs.

[0165] Specifically, the simulation prediction method for the thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall includes:

[0166] The Type56 module is called to build a building physical model, set the internal space structure and enclosure structure parameters, as well as the heat dissipation of indoor personnel, lights and equipment. It is also used to simulate the dynamic heat transfer process of the cement mortar layer on both sides of the phase change self-insulating block wall and the heat transfer process of the internal space; the Type9e module is called to read the external weather data of the building; the Type451 module is called to simulate the dynamic heat transfer process of the non-homogeneous phase change composite self-insulating block; the Type65a module is called to output the simulation calculation results as data files and display the results.

[0167] Table 2 Physical parameters of building model envelope materials

[0168]

[0169] The embodiment of the present invention also provides an office in a typical office building using non-homogeneous phase change composite self-insulating blocks as a building model for thermal environment simulation prediction, wherein the building model is as follows: Figure 8 As shown. Specifically, the southwest corner room is modeled separately in TRNSYS, and the heat dissipation of indoor personnel, lighting fixtures and equipment is set to 0W. The thermal environment of the room is only affected by the outdoor environment. According to the engineering case, the floor height of the room model is set to 4m, and the length and width are both set to 6m. The non-uniform phase change composite self-insulating block wall is only built on the west wall. The block size is still 390 mm×280mm×190 mm, and the heat transfer area of ​​the west wall is 24m 2 , block wall structure such as Fig. 9 The structural and physical parameters of other exterior walls, interior walls, roofs and floors are shown in Table 2. In addition, the south exterior window is 1m high from the floor, with a size of 2.7m×2.4m, a window frame ratio of 15%, and the window structural and physical parameters are shown in Table 3.

[0170] Table 3 Window thermal performance parameters

[0171]

[0172] Specifically, in the embodiment of the present invention, regarding outdoor conditions, the meteorological data of a typical meteorological year in Xi'an, Asia\CN-Xian-570360.tm2, is used in the simulation, and the simulation time is selected from June 1st to June 15th, a total of 360 hours, and the average outdoor air temperature is 27.26°C.

[0173] Specifically, the embodiment of the present invention analyzes the influence of phase change temperature on the thermal environment control effect of the non-homogeneous phase change composite self-insulating wall room, and takes the PCM latent heat at different phase change temperatures as the same value of 215kJ / kg. Since the phase change temperature of PCM is an interval, the phase change temperature of PCM is referred to as the mean value of the phase change interval in the embodiment of the present invention. For example, the PCM phase change temperature is 24 The phase transition interval of PCM is [22 ,26 ], and the other parameters are shown in Table 4.

[0174] Table 4 Physical parameters of phase change materials

[0175]

[0176] Specifically, in the embodiment of the present invention, the TRNSYS software simulation interface construction process is as follows: Fig.10 As shown, the module information used in the simulation process is shown in Table 5.

[0177] Table 5 TRNSYS simulation module description

[0178]

[0179] Specifically, in the embodiment of the present invention, the above simulation prediction results in the following Fig.11 and Fig.12 The simulation results reveal the influence of phase change materials with different phase change temperatures on the room temperature stability of non-homogeneous phase change composite self-insulating masonry structure rooms from June 1 to June 15. The advantages and disadvantages of the four PCM phase change temperatures are compared using the wall surface temperature and indoor air temperature as evaluation indicators. Whether it is the wall surface temperature or the indoor air temperature, when the PCM phase change temperature is 24℃, the stability of these two evaluation indicators is the best, indicating that the phase change temperature of 24-28℃ is the appropriate phase change temperature design condition.

[0180] The embodiment of the present invention also provides a system for predicting the thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall. Fig.13As shown, it specifically includes: firstly, a non-homogeneous phase change composite building block physical model module is built to establish a non-homogeneous phase change composite self-insulation building block heat transfer model; then a parameter solving method and process are established, and a Type451 module is built in the TRNSYS software in combination with the non-homogeneous phase change composite self-insulation building block heat transfer model to accurately simulate the winter heat transfer process of the building block; the self-built Type451 module is coupled with the standard module in the TRNSYS software to obtain a non-homogeneous phase change building block room simulation module; based on the non-homogeneous phase change building block room simulation module , predicts the thermal environment of the room with non-uniform phase change composite self-insulating block wall, and is used to simulate the influence of different environmental conditions, building structures and material properties on the thermal environment of the room; the system solves the problem that the heat transfer simulation technology of non-uniform phase change blocks cannot be comprehensively simulated and solved with the internal thermal environment simulation of the room. It can simulate the indoor environmental temperature index of the room when block walls are used for walls of different room types and different orientations, and can more accurately simulate the heat transfer process of non-uniform phase change blocks, thereby optimizing the design of non-uniform phase change composite self-insulating blocks and block walls.

[0181] Based on the same inventive concept as the aforementioned embodiment, an embodiment of the present invention further provides an electronic device, including a memory and a processor, the memory being used to store computer-executable instructions, and the processor being used to execute computer-executable instructions, to implement a method for predicting the thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall as proposed in the aforementioned embodiment.

[0182] This embodiment also provides a computer-readable storage medium on which a computer program is stored. When the program is executed by the processor, it is used to establish a heat transfer model of a non-homogeneous phase change composite self-insulating block and predict the indoor temperature of a room with a block wall. The storage medium can be any non-volatile storage device such as a hard disk, a solid-state hard disk, a flash drive, an optical disk, etc., which is used to store computer program codes and necessary data files. The stored computer program includes: a model establishment module, a self-built dynamic heat transfer module simulation module, and a room thermal environment prediction module.

[0183] Finally, it should be pointed out that the above specific embodiments are only representative examples of the present invention. Obviously, the present invention is not limited to the above specific embodiments, and there are many variations. Any simple modification, equivalent changes and modifications made to the above specific embodiments based on the technical essence of the present invention should be considered to belong to the protection scope of the present invention.

Claims

1. A method for predicting the thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall, characterized in that: include: Establishing a heat transfer model of a non-homogeneous phase-change composite self-insulating building block, including: establishing a physical model and a solution unit of the non-homogeneous phase-change composite self-insulating building block, wherein the holes on one side of the non-homogeneous phase-change composite self-insulating building block are filled with phase-change materials, and the holes on the other side are filled with insulation materials; based on the physical model and the solution unit, establishing a control equation including a heat transfer process of the phase-change material part of the building block and a heat conduction process of the base material and insulation material area of ​​the building block; establishing boundary conditions for the heat transfer process of the non-homogeneous phase-change composite self-insulating building block; constructing a solution method and process based on the control equation and the boundary conditions, analyzing the dynamic heat transfer process of the non-homogeneous phase-change composite self-insulating building block, and obtaining the temperature field distribution of the overall space of the building block; Based on the heat transfer model of the inhomogeneous phase change composite self-insulating building blocks, a dynamic heat transfer module of the inhomogeneous phase change composite self-insulating building blocks was built in the building thermal process / energy consumption simulation software using programming language; The dynamic heat transfer module of the non-homogeneous phase-change composite self-insulating building block is coupled with the building module, the data reading module and the simulation output module in the building thermal process / energy consumption simulation software, which can simulate the heat transfer of the wall homogeneous material layer and the indoor heat transfer process, to obtain the non-homogeneous phase-change composite self-insulating building block wall room simulation module of the building thermal process / energy consumption simulation software; Based on the non-homogeneous phase change block room simulation module of the building thermal process / energy consumption simulation software, the thermal environment of the room with non-homogeneous phase change composite self-insulating block walls is predicted.

2. A method for predicting thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall according to claim 1, characterized in that: The heat transfer model of the inhomogeneous phase change composite self-insulating building block is established, including: A physical model of a non-homogeneous phase-change composite self-insulating building block is established, wherein the holes of the non-homogeneous phase-change composite self-insulating building block close to the indoor side are filled with phase-change materials, and the holes of the outdoor side are filled with insulation materials; a representative target area is selected based on the physical model to establish a solution unit; Based on the physical model and solution unit, the heat flow conduction from the high temperature side to the low temperature side of the enclosure structure and the lateral heat transfer between different materials are considered, and the control equation reflecting the heat transfer process of the heterogeneous phase change composite self-insulating block is established, which is as follows: According to the law of conservation of energy, the calculation process of specific heat capacity and latent heat of phase change is integrated into one item, In order to comprehensively calculate the heat capacity, the heat transfer process of the phase change material part of the building block is described by the following governing equation: ， Among them, x is the length direction of the block, y is the thickness direction of the block, is the thermal conductivity of the phase change material, is the temperature, t is the time variable, is the density of the phase change material, L is the latent heat of phase change, is the specific heat capacity of the phase change material, T s and T l Respectively represent the phase change starting temperature and phase change ending temperature of the phase change material; For the base material and insulation material area of ​​the block, the heat transfer process is a pure heat conduction process. The material properties do not consider the real-time changes of the non-steady-state heat transfer process. The heat transfer process is described by the following control equation: ， Among them, x is the length direction of the block, y is the thickness direction of the block, is the thermal conductivity of the substrate / insulation material, is the temperature, t is the time variable, is the specific heat capacity of the substrate / insulation material, is the density of the substrate / insulation material; The boundary conditions of the heat transfer process of the non-homogeneous phase change composite self-insulating building blocks are established, including: based on the heat transfer process of the building blocks, the boundary conditions on both sides of the length direction x of the building blocks are described by the differential formula of the adiabatic boundary conditions; for the boundary conditions on both sides of the thickness direction y of the building blocks, the differential formula of the constant heat flux boundary conditions that reflects the correlation between the heat flux density on the outdoor and indoor surfaces of the building blocks and the thermal conductivity of the substrate is used to describe.

3. A method for predicting thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall according to claim 2, characterized in that: Build solution methods and processes based on governing equations and boundary conditions, including: The heat capacity method is used to process the specific heat capacity of phase change materials under different temperature conditions. The formula is as follows: ， The linear interpolation method is used to obtain the density and thermal conductivity of the phase change material at different temperatures. The calculation formula is as follows: ， ， in, , are the densities of the phase change material in pure solid and pure liquid states respectively; , are the thermal conductivity of the phase change material in pure solid and pure liquid states respectively; The calculation formula of the liquid phase fraction f is shown as follows. When f = 0, the phase change material is in a pure solid state; when f = 1, the phase change material is in a pure liquid state; when 0 < f < 1, the phase change material is in a mixed state. ， Where f is the liquid phase ratio; T is the real-time temperature of the phase change material, °C; T s and T l They represent the phase change starting temperature and phase change ending temperature of the phase change material respectively.

4. The method for predicting the thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall according to claim 1, characterized in that: In the analysis of the dynamic heat transfer process of the heterogeneous phase change composite self-insulating block, the control equations and boundary conditions of the heat transfer model are discretized in the spatial dimension, including: The finite volume method is used to perform spatial division of the control volume for the solution unit of the heterogeneous phase change composite self-insulating block. The center of each control volume of the solution unit of the heterogeneous phase change composite self-insulating block is denoted as the calculation node of each control volume; the control equation is converted into a discrete equation at each node, and the difference of the spatial derivative is calculated based on the temperature value at the end of each time layer, and time discretization in implicit format is performed to obtain the temperature of each control volume, and then the temperature field distribution of the overall space of the block is obtained.

5. The method for predicting the thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall according to claim 1, characterized in that: Build a dynamic heat transfer module for the heterogeneous phase change self-insulating block in the building heat process / energy consumption simulation software, including: Create a variable table for module definition in the building heat process / energy consumption simulation software, including input items, output items, and parameter items; Combined with the module definition variable table, use a programming language to write the heat transfer model of the heterogeneous phase change composite self-insulating block and its dynamic heat transfer analysis process into program source code; Use the building heat process / energy consumption simulation software to compile the program source code into a directly loadable dynamic link library file; Based on the directly loadable dynamic link library file, build a dynamic heat transfer module for the heterogeneous phase change self-insulating block.

6. The method for predicting thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall according to claim 1, characterized in that: Couple the dynamic heat transfer module of the heterogeneous phase change composite self-insulating block with the building module, data reading module, and simulation output module in the building heat process / energy consumption simulation software, which can simulate the heat transfer process of the homogeneous material layer of the wall and the indoor heat transfer process, specifically including: The dynamic heat transfer module of the heterogeneous phase change composite self-insulating block takes the average temperatures of the outdoor side and indoor side surfaces of the composite self-insulating block obtained by calculation as output data, and outputs them to the building module respectively as the boundary temperature input data of the interface between the plaster mortar layer on both sides of the block and the block; for the plaster mortar layers on both sides of the block, use the building module to simulate its heat transfer process and output the surface heat flux density of the interface between the plaster mortar layer and the block as the input heat flux density of the indoor and outdoor surfaces of the heterogeneous phase change composite self-insulating block, which are respectively transmitted to the surface heat flux density input end of the dynamic heat transfer module.

7. The method for predicting thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall according to claim 1, characterized in that: Based on the heterogeneous phase change block room simulation module of the building heat process / energy consumption simulation software, predict the thermal environment of the room with the heterogeneous phase change composite self-insulating block wall, including: Use the building module to build a building model, set parameters including the internal space structure and envelope structure, as well as the heat dissipation of internal personnel, equipment, and lighting. At the same time, it is also used to simulate the dynamic heat transfer process of the cement mortar layers on both sides of the phase change self-insulating block wall and the heat transfer process of the internal space; use the data reading module to read the external weather data of the building; use the dynamic heat transfer module to simulate the dynamic heat transfer process of the heterogeneous phase change composite self-insulating block; use the simulation output module to output the predicted results of the room thermal environment simulation as data files and result displays.

8. A thermal environment prediction system for a room with a non-homogeneous phase-change composite self-insulating block wall, characterized in that: Including: Establish a model module and establish a heat transfer model of a non-homogeneous phase change composite self-insulating building block, including: establishing a physical model and a solution unit of the non-homogeneous phase change composite self-insulating building block; based on the physical model and the solution unit, establishing a control equation including the heat transfer process of the phase change material part of the building block and the heat conduction process of the base material and insulation material area of ​​the building block; establishing boundary conditions for the heat transfer process of the non-homogeneous phase change composite self-insulating building block; constructing a solution method and process based on the control equation and boundary conditions to analyze the dynamic heat transfer process of the non-homogeneous phase change composite self-insulating building block; Dynamic heat transfer module, based on the heat transfer model of heterogeneous phase change composite self-insulating building blocks, uses programming language to build a dynamic heat transfer module of heterogeneous phase change composite self-insulating building blocks in the building thermal process / energy consumption simulation software; A simulation module is used to couple the heat transfer module of the non-homogeneous phase-change composite self-insulating building block with the building module, the data reading module and the simulation output module in the building thermal process / energy consumption simulation software that can simulate the heat transfer of the wall homogeneous material layer and the indoor heat transfer process, so as to obtain the non-homogeneous phase-change composite self-insulating building block wall room simulation module of the building thermal process / energy consumption simulation software; The room thermal environment prediction module is based on the non-homogeneous phase change block room simulation module of the building thermal process / energy consumption simulation software to predict the thermal environment of the room with non-homogeneous phase change composite self-insulating block walls.

9. An electronic device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method for predicting the thermal environment of a room with a non-homogeneous phase change composite self-insulating block wall as described in any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of a method for predicting the thermal environment of a room with a non-homogeneous phase-change composite self-insulating block wall as described in any one of claims 1 to 7 are implemented.

Images