Prediction method for passive humidifying effect of urban underground space, medium and program product

Abstract

The invention discloses a prediction method for a passive humidity control effect of an urban underground space, a medium and a program product. The prediction method comprises the following steps: acquiring a pre-constructed humidity control material heat and humidity dynamic transfer calculation model; obtaining a pre-constructed underground soil dynamic heat transfer calculation model and an underground space indoor heat and humidity balance calculation model; and coupling the dynamic heat and humidity transfer calculation model of the humidity control material, the dynamic heat transfer calculation model of the underground soil and the indoor heat and humidity balance calculation model of the underground space, and predicting the passive humidity control effect of the urban underground space. According to the invention, the dynamic moisture absorption and desorption process of the moisture conditioning material in the actual underground engineering can be reflected more accurately.

Description

Technical Field

[0001] The present invention belongs to the technical field of low-carbon technologies for urban underground spaces, and particularly relates to a method, a medium, and a program product for predicting the passive humidity regulation effect in urban underground spaces. Background Art

[0002] The thermal and humidity environments in urban underground spaces face problems such as complex and variable heat and moisture transfer. The indoor air quality management is comprehensively affected by multiple factors such as outdoor climate characteristics, ventilation conditions, energy conservation, and consumption reduction. Using humidity-regulating materials for passive humidity regulation in underground spaces does not require the consumption of non-renewable energy. Through the moisture absorption and desorption of the materials themselves, the indoor relative humidity fluctuations can be automatically adjusted without consuming external energy, reducing the moisture load of the air-conditioning system and achieving the purpose of energy conservation. Although certain progress has been made in the research and development of humidity-regulating materials, their actual applications are still relatively few, and their stability and durability in long-term use remain to be further verified.

[0003] Currently, there are three methods for studying the thermal and humidity environments in underground spaces: numerical simulation methods, on-site observation and monitoring, and laboratory experiments. (1) Numerical simulation methods use mathematical models and computer simulations to study the thermal and moisture transfer processes in underground spaces; (2) Engineering measurements collect and monitor real-time data inside and around underground spaces by deploying sensors and monitoring equipment. However, this method has a high cost and cannot cover all types of underground spaces and changing conditions, lacking universality; (3) Laboratory experiments conduct thermal and moisture transfer experiments in simulated experimental environments to verify the accuracy and reliability of numerical models. However, it is difficult to fully simulate the real and complex environmental conditions in underground spaces, and there are problems such as scale effects and the applicability of experimental results.

[0004] Based on the above, the numerical simulation method is selected to establish a dynamic prediction model for passive humidity regulation in underground spaces using humidity-regulating materials, study the suitability of humidity-regulating materials in the application of humidity regulation in underground spaces, and how to control the indoor thermal and humidity environment within a safe and reasonable range without mechanical dehumidification. Summary of the Invention

[0005] In view of the above problems, the present invention proposes a method, a medium, and a program product for predicting the passive humidity regulation effect in urban underground spaces, which can more accurately reflect the dynamic moisture absorption and desorption process of humidity-regulating materials inside actual projects.

[0006] To achieve the above technical objectives and effects, the present invention is realized through the following technical solutions:

[0007] In a first aspect, the present invention provides a method for predicting the passive humidity regulation effect in urban underground spaces, including:

[0008] Obtain a pre-built calculation model for the thermal and moisture dynamic transfer of the moisture conditioning material;

[0009] Obtain a pre-built calculation model for the dynamic heat transfer of underground soil and a calculation model for the indoor heat and moisture balance of underground spaces;

[0010] Couple the calculation model for the thermal and moisture dynamic transfer of the moisture conditioning material with the calculation model for the dynamic heat transfer of underground soil and the calculation model for the indoor heat and moisture balance of underground spaces to predict the passive moisture conditioning effect in urban underground areas.

[0011] Combined with the first aspect, optionally, the calculation model for the thermal and moisture dynamic transfer of the moisture conditioning material includes a thermal and moisture coupled transfer control equation and thermal and moisture coupled transfer boundary conditions;

[0012] The thermal and moisture coupled transfer control equation includes a heat transfer control equation and a moisture transfer control equation, and the heat transfer control equation and the moisture transfer control equation establish a thermal and moisture coupled transfer relationship through the latent heat of phase change;

[0013] The thermal and moisture coupled transfer boundary conditions include the moisture absorption and desorption boundary condition on the inner surface of the moisture conditioning material and the heat boundary condition on the inner surface of the moisture conditioning material.

[0014] Combined with the first aspect, optionally, the mathematical expression of the heat transfer control equation is:

[0015] ;

[0016] The mathematical expression of the moisture transfer control equation is:

[0017] ;

[0018] In the formula, is the moisture content of the air in the pores of the material, is the temperature inside the material, is the time; is the density of the material, is the specific heat capacity of the material, is the thermal conductivity of the material, is the latent heat of vaporization; is the moisture storage coefficient related to the moisture content, is the moisture storage coefficient related to the temperature, is the water vapor diffusion coefficient, is the position coordinate;

[0019] The mathematical expression of the moisture absorption and desorption boundary condition on the inner surface of the moisture conditioning material is:

[0020] ;

[0021] In the formula, is the moisture content of the indoor air in the underground space, i.e., the moisture content of the air in the external environment where the material is located; is the moisture content of the air in the thin layer on the inner surface of the material; is the convective moisture transfer coefficient;

[0022] The calculation formula for the thermal boundary condition of the inner surface of the moisture conditioning material is:

[0023] ;

[0024] In the formula, is the convective heat transfer coefficient of the inner surface of the material, is the indoor air temperature in the underground space, i.e., the air temperature in the external environment where the material is located; is the air temperature in the thin layer on the inner surface of the material.

[0025] Combined with the first aspect, optionally, the moisture storage coefficient related to the moisture content includes the moisture storage coefficient during the moisture absorption process and the moisture storage coefficient during the moisture release process ; and

[0026] The calculation formula for the moisture storage coefficient during the moisture absorption process is:

[0027] ;

[0028] The calculation formula for the moisture storage coefficient during the moisture release process is:

[0029] ;

[0030] In the formula, is the moisture content in the saturated state; is the minimum value of the relative humidity during the moisture absorption and release cycle; is the maximum value of the relative humidity during the moisture absorption and release cycle; is the isothermal equilibrium moisture content of the moisture absorption curve at a given relative humidity; is the isothermal equilibrium moisture content of the moisture release curve at a given relative humidity; is the moisture storage coefficient calculated from the isothermal moisture absorption curve, is the moisture storage coefficient calculated from the isothermal moisture release curve.

[0031] Combined with the first aspect, optionally, under different temperature and humidity environments, the moisture storage coefficient during the moisture absorption process and the moisture storage coefficient during the moisture release process of the moisture conditioning material are continuously updated through the Mualem experimental model.

[0032] In combination with the first aspect, optionally, the underground soil dynamic heat transfer calculation model includes an internal heat migration control equation;

[0033] The calculation model of indoor heat and moisture balance in the underground space includes the heat balance equation and mass balance equation of indoor air, as well as the boundary conditions of the heat transfer calculation domain in the underground space.

[0034] In combination with the first aspect, optionally, the mathematical expression of the internal heat migration control equation is:

[0035] ;

[0036] In the formula, is the soil temperature value at time at point is the density of the soil, is the specific heat capacity of the soil, is the soil thermal conductivity at point and is the material thermal conductivity at point ;

[0037] The mathematical expression of the heat balance equation of the indoor air is:

[0038] ;

[0039] The mathematical expression of the mass balance equation of the indoor air is:

[0040] ;

[0041] In the formula, is the specific heat at constant pressure of air, is the density of air, is the volume of the room, is the convective heat transfer coefficient, is the surface area of the room wall , is the indoor air temperature in the underground space, is the surface temperature of the calculated room wall , is the heat gain in the underground space indoor, is the moisture absorption or moisture release amount of the humidity control material attached to the surface of the room enclosure structure; is the ventilation volume of indoor and outdoor air in the underground space, is the moisture content of outdoor air; is the moisture content of indoor air in the underground space, that is, the moisture content of the external environment air where the material is located; is the water vapor generation amount in the underground space indoor; ​

[0042] The boundary conditions of the underground space heat transfer calculation domain include upper boundary conditions, lower boundary conditions, side boundary conditions, and convective boundary conditions;

[0043] The mathematical expression of the upper boundary condition is:

[0044] ;

[0045] In the formula, is the short-wave radiation reaching the ground; is the long-wave radiation reaching the ground; is the long-wave radiation emitted from the ground; is the convective heat transfer amount between the ground surface and the atmosphere; is the short-wave reflectivity of the ground surface; is the long-wave absorptivity of the ground surface; The long-wave emissivity of the ground surface;

[0046] The mathematical expression of the lower boundary condition is:

[0047] ;

[0048] In the formula, is the temperature of the lower boundary of the calculation domain, is the annual average temperature of the ground surface;

[0049] The mathematical expression of the side boundary condition is:

[0050] ;

[0051] In the formula, is the heat flux of the side boundary of the calculation domain;

[0052] The mathematical expression of the convective boundary condition is:

[0053] ;

[0054] In the formula, is the indoor air temperature of the underground space, that is, the ambient air temperature where the material is located; is the convective heat transfer coefficient of the inner surface of the material.

[0055] Combined with the first aspect, optionally, coupling the heat and moisture dynamic transfer calculation model of the humidity control material with the dynamic heat transfer calculation model of the underground soil and the calculation model of the indoor heat and moisture balance of the underground space to predict the passive humidity control effect in the city, including:

[0056] Continuously execute the preset loop calculation steps until the set requirements are met, and obtain the final indoor temperature and humidity parameters of the underground space;

[0057] The described cyclic calculation steps include:

[0058] Performing a difference operation according to the underground soil dynamic heat transfer calculation model to obtain the internal temperature distribution of the underground soil;

[0059] Taking the internal temperature distribution of the soil as the outer boundary condition of the heat and moisture dynamic transfer calculation model of the moisture conditioning material, and combining the indoor temperature and humidity parameters, calculating the heat flux and moisture flux on the inner surface of the moisture conditioning material according to the heat and moisture dynamic transfer calculation model of the moisture conditioning material;

[0060] Substituting the calculated heat flux and moisture flux on the inner surface of the moisture conditioning material into the calculation models of indoor heat and moisture balance in the underground space respectively, and discretely solving the indoor temperature and humidity parameters in the underground space.

[0061] In a second aspect, the present invention provides a computer-readable storage medium, on which a computer program is stored, and when the program is executed by a processor, it implements the method for predicting the passive moisture conditioning effect in urban underground spaces described in any one of the first aspects.

[0062] In a third aspect, the present invention provides a computer program product, including a computer program / instructions, and when the computer program / instructions are executed by a processor, it implements the method for predicting the passive moisture conditioning effect in urban underground spaces described in any one of the first aspects.

[0063] Compared with the prior art, the beneficial effects of the present invention are:

[0064] (1) Based on the Matsumoto porous medium heat and moisture coupled transfer theoretical model, combined with the dynamic heat transfer characteristics of the underground space, the present invention uses the Mualem hysteresis model to update and calculate the moisture storage coefficient during the moisture absorption and release processes under different temperature and humidity environments, so as to correct the passive moisture conditioning dynamic prediction model. Different from the fixed value of the moisture storage coefficient in the prior art, the calculation of the model is closer to the actual process, thus more accurately reflecting the internal dynamic moisture absorption and release process of the moisture conditioning material in actual engineering. (2) The calculation model proposed in the present invention is simple, has a small amount of calculation, and has a high simulation accuracy. It can analyze the influence of various factors and meet the engineering requirements.

[0065] (3) The calculation proposed in the present invention can deeply analyze the passive moisture conditioning effect in urban underground spaces and its applicability in different climate zones, and conduct long-term dynamic evaluation and optimization on it.

[0066] (3) The calculation proposed in the present invention can deeply analyze the passive moisture conditioning effect in urban underground spaces and its applicability in different climate zones, and conduct long-term dynamic evaluation and optimization on it. Description of the Drawings

[0067] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for use in the embodiments. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can also be obtained based on these drawings, where:

[0068] Figure 1 It is a composition diagram of the calculation module of the passive humidity regulation prediction model for underground space in an embodiment;

[0069] Figure 2 It is a parameter diagram of the moisture absorption and desorption hysteresis model of the humidity regulation material Mualem in an embodiment;

[0070] Figure 3 It is a simple diagram of the underground space building model and a diagram of the calculation boundary conditions in an embodiment;

[0071] Figure 4 It is a block diagram of the calculation program of the thermal - moisture coupling prediction model for underground space in an embodiment;

[0072] Figure 5 It is a comparison diagram of the calculation results of the numerical method and the measured temperature in the project in an embodiment;

[0073] Figure 6 It is a comparison diagram of the calculation results of the numerical method and the measured humidity in the project in an embodiment. Specific Embodiments

[0074] The following will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are only a part of the embodiments of the present invention, rather than all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts belong to the scope of protection of the present invention.

[0075] In addition, if there are descriptions involving "first", "second", etc. in the embodiments of the present invention, such descriptions of "first", "second", etc. are only for descriptive purposes and cannot be understood as indicating or implying their relative importance or implicitly indicating the quantity of the indicated technical features. Thus, the features defined with "first", "second" may explicitly or implicitly include at least one of such features. In addition, the technical solutions between various embodiments can be combined with each other, but it must be based on the fact that those of ordinary skill in the art can implement it. When the combination of technical solutions results in contradictions or cannot be implemented, it should be considered that such a combination of technical solutions does not exist and is not within the scope of protection required by the present invention.

[0076] Embodiment 1

[0077] In an embodiment of the present invention, a method for predicting the passive humidity control effect of urban underground space is provided, which specifically includes the following steps:

[0078] (1) Obtain a pre-constructed computational model for the thermo-hygrometric dynamic transfer of humidity control materials;

[0079] (2) Obtain a pre-constructed computational model for the dynamic heat transfer of underground soil and a computational model for the indoor heat and humidity balance in underground space;

[0080] (3) Couple the computational model for the thermo-hygrometric dynamic transfer of humidity control materials with the computational model for the dynamic heat transfer of underground soil and the computational model for the indoor heat and humidity balance in underground space to predict the passive humidity control effect of urban underground space.

[0081] Based on the prediction method proposed in the embodiment of the present invention, the proportion of indoor humidity being damp and very damp can be reduced, the indoor humidity environment can be significantly improved, and at the same time, the dehumidification latent heat load and the energy consumption of the mechanical dehumidification system operation can be reduced. The method for establishing this model can provide support for controlling the indoor heat and humidity environment within a safe and reasonable range without relying on mechanical dehumidification. It is mainly used to analyze and study the application effect of passive humidity control in underground space and the regional differences in different climate zones in China, and conduct long-term dynamic evaluation and optimal control on it to give full play to the advantages of passive humidity control. In addition, the present invention can be used for predicting the heat and humidity environment of urban underground space and calculating the heat load and humidity load.

[0082] In a specific embodiment of the embodiment of the present invention, the computational model for the thermo-hygrometric dynamic transfer of humidity control materials includes a thermo-hygrometric coupling transfer control equation and a thermo-hygrometric coupling transfer boundary condition;

[0083] The thermo-hygrometric coupling transfer control equation includes a heat transfer control equation and a moisture transfer control equation, and the heat transfer control equation and the moisture transfer control equation establish a thermo-hygrometric coupling transfer relationship through the latent heat of phase change;

[0084] The thermo-hygrometric coupling transfer boundary condition includes the moisture absorption and desorption boundary condition on the inner surface of the humidity control material and the heat boundary condition on the inner surface of the humidity control material.

[0085] In a specific embodiment of the embodiment of the present invention, the mathematical expression of the heat transfer control equation is:

[0086] ;

[0087] The mathematical expression of the moisture transfer control equation is:

[0088] ;

[0089] In the formula, is the moisture content of the air in the pores of the material, is the temperature inside the material, is time; is the density of the material, is the specific heat capacity of the material, is the thermal conductivity of the material, is the latent heat of vaporization; is the moisture storage coefficient related to the moisture content, is the moisture storage coefficient related to the temperature, is the water vapor diffusion coefficient, is the position coordinate;

[0090] The mathematical expression of the moisture absorption and desorption boundary condition on the inner surface of the moisture conditioning material is:

[0091] ;

[0092] In the formula, is the moisture content of the indoor air in the underground space, that is, the moisture content of the ambient air where the material is located, is the moisture content of the air in the thin layer on the inner surface of the material; is the convective heat and mass transfer coefficient;

[0093] The calculation formula of the thermal boundary condition on the inner surface of the moisture conditioning material is:

[0094] ;

[0095] In the formula, is the convective heat transfer coefficient of the inner surface of the material, is the indoor air temperature in the underground space, that is, the ambient air temperature where the material is located;, is the air temperature in the thin layer on the inner surface of the material.

[0096] In a specific embodiment of the embodiment of the present invention, the moisture storage coefficient related to the moisture content includes the moisture storage coefficient during the moisture absorption process and the moisture storage coefficient during the moisture desorption process;

[0097] The calculation formula of the moisture storage coefficient during the moisture absorption process is:

[0098] ;

[0099] The calculation formula of the moisture storage coefficient during the moisture desorption process is:

[0100] ;

[0101] In the formula, is the moisture content in the saturated state; is the minimum value of relative humidity during the moisture absorption and desorption cycle; is the maximum value of relative humidity during the moisture absorption and desorption cycle; is the isothermal equilibrium moisture content of the moisture absorption curve at a given relative humidity; is the isothermal equilibrium moisture content of the moisture desorption curve at a given relative humidity; is the moisture storage coefficient calculated from the isothermal moisture absorption curve, is the moisture storage coefficient calculated from the isothermal moisture desorption curve.

[0102] In a specific embodiment of the present invention, under different temperature and humidity environments, the moisture storage coefficient of the moisture absorption process of the humidity control material is continuously updated through the Mualem experimental model and the moisture storage coefficient of the moisture desorption process .

[0103] In a specific embodiment of the present invention, the underground soil dynamic heat transfer calculation model includes an internal heat migration control equation;

[0104] The calculation model of the indoor heat and moisture balance in the underground space includes the heat balance equation and the mass balance equation of the indoor air, as well as the boundary conditions of the heat transfer calculation domain in the underground space.

[0105] In a specific embodiment of the present invention, the mathematical expression of the internal heat migration control equation is:

[0106] ;

[0107] In the formula, is the soil temperature value at time at point is the density of the soil, is the specific heat capacity of the soil, is the soil thermal conductivity at point and is the material thermal conductivity at point ;

[0108] The mathematical expression of the heat balance equation of the indoor air is:

[0109] ;

[0110] The mathematical expression of the mass balance equation of the indoor air is:

[0111] ;

[0112] In the formula, is the specific heat capacity at constant pressure of the air, is the density of the air, is the volume of the room, is the convective heat transfer coefficient, is the surface area of the room wall ; is the indoor air temperature in the underground space, is to calculate the surface temperature of the room wall ; is the heat gain in the underground space indoors, is the moisture absorption or moisture release of the humidity regulating material attached to the surface of the room enclosure structure; is the ventilation volume of indoor and outdoor air in the underground space, is the moisture content of outdoor air; is the moisture content of indoor air in the underground space, that is, the moisture content of the external ambient air where the material is located; is the water vapor generation amount in the underground space indoors;

[0113] The boundary conditions of the heat transfer calculation domain of the underground space include upper boundary conditions, lower boundary conditions, side boundary conditions and convective boundary conditions;

[0114] The mathematical expression of the upper boundary condition is:

[0115] ;

[0116] In the formula, is the short-wave radiation reaching the ground; is the long-wave radiation reaching the ground; is the long-wave radiation emitted by the ground; is the convective heat transfer amount between the ground surface and the atmosphere; is the short-wave reflectivity of the ground surface; is the long-wave absorptivity of the ground surface; is the long-wave emissivity of the ground surface;

[0117] The mathematical expression of the lower boundary condition is:

[0118] ;

[0119] In the formula, is the temperature of the lower boundary of the calculation domain, is the annual average temperature of the ground surface;

[0120] The mathematical expression of the side boundary condition is:

[0121] ;

[0122] In the formula, is the heat flux of the side boundary of the calculation domain;

[0123] The mathematical expression of the convective boundary condition is as follows:

[0124] ;

[0125] In the formula, is the indoor air temperature of the underground space, that is, the ambient air temperature of the external environment where the material is located, is the convective heat transfer coefficient of the inner surface of the material.

[0126] In a specific implementation manner of the embodiment of the present invention, coupling the heat and moisture dynamic transfer calculation model of the humidity-adjusting material with the underground soil dynamic heat transfer calculation model and the indoor heat and moisture balance calculation model of the underground space to predict the passive humidity-adjusting effect of the urban underground includes:

[0127] Continuously execute the preset loop calculation steps until the set requirements are met, and the final indoor temperature and humidity parameters of the underground space are obtained;

[0128] The loop calculation steps include:

[0129] Perform difference operation according to the underground soil dynamic heat transfer calculation model to obtain the internal temperature distribution of the underground soil;

[0130] Take the internal temperature distribution of the soil as the outer boundary condition of the heat and moisture dynamic transfer calculation model of the humidity-adjusting material, and combine the indoor temperature and humidity parameters. According to the heat and moisture dynamic transfer calculation model of the humidity-adjusting material, calculate the heat flux and moisture flux on the inner surface of the humidity-adjusting material;

[0131] Substitute the calculated heat flux and moisture flux on the inner surface of the humidity-adjusting material into the indoor heat and moisture balance calculation models of the underground space respectively, and discretely solve the indoor temperature and humidity parameters of the underground space.

[0132] In the specific implementation process, the prediction method includes:

[0133] Step (1): Obtain the pre-constructed heat and moisture dynamic transfer calculation model of the humidity-adjusting material; the humidity-adjusting material is a porous medium material, isotropic, and non-deformable;

[0134] Step (2): Obtain the pre-constructed underground soil dynamic heat transfer calculation model and the indoor heat and moisture balance calculation models of the underground space;

[0135] Step (3): Couple the heat and moisture dynamic transfer calculation model of the humidity-adjusting material with the underground soil dynamic heat transfer calculation model and the indoor heat and moisture balance calculation models of the underground space to predict the passive humidity-adjusting effect of the urban underground.

[0136] Among them, the heat and moisture dynamic transfer calculation model of the humidity conditioning material is used for calculating the heat and moisture coupled transfer of the humidity conditioning material. The establishment process and calculation process of the heat and moisture dynamic transfer calculation model of the humidity conditioning material are specifically as follows:

[0137] Step 1-1: Determine the heat and moisture coupled transfer control equation of the humidity conditioning material based on the Matsumoto porous medium heat and moisture coupled transfer theoretical model;

[0138] Step 1-2: Determine the initial conditions and the input parameters of the heat and moisture dynamic transfer calculation model of the humidity conditioning material; the input parameters include the building shape size and height, the soil covering thickness, the composition and physical properties of the enclosure structure, the physical properties of the soil around the project, the heat and moisture sources inside the project, the ventilation volume, the outdoor meteorological parameters, the area and thickness of the humidity conditioning material, and the heat and moisture parameters of the humidity conditioning material; the heat and moisture coupled transfer control equation includes a heat transfer control equation and a moisture transfer control equation, and the heat and moisture coupled transfer relationship is established through the latent heat of phase change between the heat transfer control equation and the moisture transfer control equation; the mathematical expression of the heat transfer control equation is:

[0139] ;

[0140] The mathematical expression of the moisture transfer control equation is:

[0141] ;

[0142] In the formula, is the moisture content of the air in the pores of the material, is the temperature inside the material, is the time; is the density of the material, is the specific heat capacity of the material, is the thermal conductivity of the material, is the latent heat of vaporization; is the moisture storage coefficient related to the moisture content, is the moisture storage coefficient related to the temperature, is the water vapor diffusion coefficient, is the position coordinate;

[0143] Step 1-3: Determine the boundary conditions, and the boundary conditions include the moisture absorption and desorption boundary conditions and the heat boundary conditions;

[0144] The mathematical expression of the moisture absorption and desorption boundary condition on the inner surface of the humidity conditioning material is:

[0145] ;

[0146] In the formula, is the moisture content of the indoor air in the underground space, that is, the moisture content of the air in the external environment where the material is located, is the moisture content of air at the thin layer of the inner surface of the material; is the convective moisture transfer coefficient;

[0147] The calculation formula for the thermal boundary condition of the inner surface of the moisture conditioning material is:

[0148] ;

[0149] In the formula, is the convective heat transfer coefficient of the inner surface of the material, is the indoor air temperature of the underground space, that is, the ambient air temperature where the material is located;, is the air temperature at the thin layer of the inner surface of the material;

[0150] Steps 1-4, update the thermo-hygroscopic physical properties parameters of the moisture conditioning material based on the Mualem experimental model. The thermo-hygroscopic physical properties parameters include the moisture storage coefficient ; Specifically: Because the hysteresis effect of the moisture release process of the moisture conditioning material relative to the moisture absorption process needs to be considered. Under different temperature and humidity environments, the moisture storage coefficient involved in the moisture absorption and release process of the moisture conditioning material is continuously updated through the Mualem experimental model The moisture storage coefficient during the moisture absorption process of the moisture conditioning material and the moisture storage coefficient during the moisture release process

[0151] ;

[0152] ;

[0153] In the formula, is the moisture content under the saturated state; is the minimum value of the relative humidity in the moisture absorption and release cycle; is the maximum value of the relative humidity in the moisture absorption and release cycle; is the isothermal equilibrium moisture content of the moisture absorption curve at a given relative humidity; is the isothermal equilibrium moisture content of the moisture release curve at a given relative humidity; is the moisture storage coefficient calculated from the isothermal moisture absorption curve, is the moisture storage coefficient calculated from the isothermal moisture release curve; The moisture storage coefficients and in the moisture absorption process are related to the moisture storage coefficient in the heat transfer control equation and the moisture transfer control equation, and are continuously updated and iterated according to the temperature and humidity conditions of the environment;

[0154] Steps 1-5, according to the given initial conditions and boundary conditions, use the finite difference method for numerical calculation to obtain the temperature and humidity distributions inside the moisture conditioning material, as well as the heat flux and moisture flux on the surface of the moisture conditioning material.

[0155] The establishment process and calculation process of the above-mentioned underground soil dynamic heat transfer calculation model and the calculation model of indoor heat and moisture balance in the underground space are specifically as follows:

[0156] Step 2-1, determine the underground soil dynamic heat transfer calculation model, simplify the heat transfer problem between the soil around the underground space into a two-dimensional unsteady heat conduction problem, and the internal heat transfer control equation is as follows:

[0157]

[0158] In the formula, is the soil temperature value at time at point is the density of the soil, is the specific heat capacity of the soil, is the soil thermal conductivity at point ; is the material thermal conductivity at point .

[0159] Step 2-2, determine the boundary conditions of the heat transfer calculation domain of the underground space, including:

[0160] (1) For the upper boundary of the calculation domain, the energy balance of the ground surface can be expressed by the following formula:

[0161]

[0162] In the formula, is the short-wave radiation reaching the ground; is the long-wave radiation reaching the ground; is the long-wave radiation emitted by the ground; is the convective heat transfer amount between the ground surface and the atmosphere; is the short-wave reflectivity of the ground surface; is the long-wave absorptivity of the ground surface; is the long-wave emissivity of the ground surface.

[0163] (2) The lower boundary of the calculation domain (10 m below the floor of the underground building) is an isothermal boundary, and the value is equal to the annual average temperature of the ground surface:

[0164]

[0165] (3) The side boundaries of the calculation domain are adiabatic boundaries, that is, the heat flux passing through is 0:

[0166]

[0167] (4) The internal boundaries (ceiling, floor, vertical side walls) between the enclosure structure and the indoor air are convective boundaries:

[0168] ;

[0169] In the formula, is the indoor air temperature of the underground space, that is, the ambient air temperature where the material is located.

[0170] Step 2-3: Determine the calculation models for indoor heat and moisture balance in the underground space, including the heat balance equation and mass balance equation of indoor air;

[0171] The mathematical expression of the heat balance equation of the indoor air is:

[0172]

[0173] In the formula, is the specific heat at constant pressure of air, is the density of air, is the volume of the room, is the convective heat transfer coefficient, is the wall of the room surface area, is the indoor air temperature of the underground space, is to calculate the wall of the room surface temperature, is the heat gain in the indoor of the underground space, is the moisture absorption or moisture release amount of the humidity control material attached to the surface of the room enclosure structure; is the ventilation volume of the indoor and outdoor air in the underground space, is the moisture content of the outdoor air; is the moisture content of the indoor air in the underground space, that is, the moisture content of the ambient air where the material is located; is the indoor water vapor generation amount in the underground space;

[0174] The mathematical expression of the mass balance equation of the indoor air is:

[0175]

[0176] In the formula, is the ventilation volume of the room and the outdoor air, is the moisture content of the outdoor air, is the moisture content of the indoor air, is the room indoor water vapor generation amount.

[0177] Step 2-4: Use the finite difference method for numerical calculation to obtain the underground soil temperature distribution and the indoor temperature and humidity changes in the underground space.

[0178] The specific content of step (3) includes:

[0179] Step 3-1: Input the parameters of the underground soil dynamic heat transfer calculation model. At a given time step, perform a difference operation according to the underground soil dynamic heat transfer calculation model to obtain the internal temperature distribution of the underground soil (i.e.,).

[0180] Step 3-2: Take the soil temperature distribution as the external boundary condition for the thermo-hygroscopic coupling transfer of the moisture conditioning material, and combine the indoor temperature and humidity parameters. According to the heat transfer and moisture transfer control equations of the thermo-hygroscopic coupling transfer of the moisture conditioning material, calculate the heat flux and moisture flux on the inner surface of the moisture conditioning material.

[0181] Step 3-3: Substitute the calculated heat flux and moisture flux on the inner surface of the moisture conditioning material into the heat balance equation and mass balance equation of the indoor air respectively, and discretely solve the indoor temperature and humidity parameters to complete one cycle.

[0182] Step 3-4: Then start the calculation of the next time step and perform the next coupling calculation until the set requirements are met, at which point the calculation ends and the indoor temperature and humidity changes in the underground space are output.

[0183] Figure 1 The flowchart of the prediction method in the embodiment of the present invention is shown. The following is a further description:

[0184] 1. Determine the physical properties and thermo-hygroscopic parameters of the moisture conditioning material

[0185] The internal thermo-hygroscopic coupling transfer process of the moisture conditioning material is treated as a one-dimensional process. Usually, the moisture conditioning material is attached to the inner surface of the enclosure structure, and the height and width of the moisture conditioning material are much larger than its thickness. Therefore, the internal thermo-hygroscopic transfer process can be considered a one-dimensional physical process. The moisture conditioning material is a continuous and uniform porous medium material, isotropic, and non-deformable. There is no other thermo-hygroscopic source inside it, and local thermal and moisture equilibrium always exists inside. The moist air in the pores of the moisture conditioning material is treated as an incompressible ideal gas and satisfies the ideal gas state equation. The heat flux and mass flux at the interface of different materials are continuous, and the influence of thermo-hygroscopic resistance is ignored.

[0186] Taking the bio-based matrix moisture conditioning material as an example, the thickness of this material is 20 mm, the density is 550 kg / m 3 , and the thermal conductivity is 0.07 W / (m·K). The parameters related to moisture transfer, the moisture storage coefficient and the water vapor diffusion coefficient are shown in Table 1. The surface mass transfer coefficient of the biomass moisture conditioning material is 0.0179 kg / (m 2 s) ~ 0.0165 kg / (m 2 s).

[0187] Table 1 Moisture storage coefficients of biomass moisture conditioning materials in different relative humidity ranges ( value and value

[0188]

[0189] The hysteresis of moisture absorption and desorption affects the moisture absorption and desorption coefficient related to the moisture content of the moisture conditioning material to a certain extent, and will affect the simulation of the internal heat and moisture coupling migration process of the moisture conditioning material. Therefore, in each cycle calculation process, it is necessary to update the moisture storage coefficient of the moisture conditioning material using the Mualem experimental model, as shown in Figure 2 .

[0190] 2. Determine the thermal parameters of the underground space enclosure structure

[0191] See Figure 3 , the soil cover thickness of the underground project is 1.5 m, the building storey height is 3.3 m, the surrounding is rammed clay, and the thermal conductivity of the moisture conditioning material = 1.16 W / m°C, the density = 2000 kg / m 3 specific heat = 1010 J / kg°C, the thermal diffusivity a = 0.0037 m 2 . The ventilation rate inside the project is 0.15 h -1 . The building enclosure structure and thermal parameters are shown in Table 2

[0192] Table 2 Structure and thermal parameters of the enclosure structure of the underground space building model

[0193]

[0194] 3. Establishment of the thermal and moisture dynamic transfer calculation model for the moisture conditioning material

[0195] Based on the Matsumoto theoretical model of coupled heat and moisture transfer in porous media, a calculation model for coupled heat and moisture transfer of humidity control materials is established. First, the control equations for moisture transfer and heat transfer of the dynamic calculation model of coupled heat and moisture transfer of humidity control materials are established. Since the humidity control material is a porous medium material, simplifying assumptions can be made for the mathematical and physical model of coupled heat and moisture transfer of porous medium materials. Temperature and moisture content of air are selected as the driving forces for heat and moisture transfer, and the control equation for moisture transfer in porous media is calculated. According to the law of conservation of energy, its heat transfer control equation is determined. Then, since the problem of rainwater migrating and penetrating into the indoor through the enclosure structure does not need to be considered in the underground space, and under normal conditions, it is basically very difficult for the moisture in the surrounding soil to diffuse into the project interior, the control equation for coupled heat and moisture transfer can be simplified to the control equation for coupled heat and moisture transfer with only water vapor diffusion inside the material. Subsequently, the boundary conditions of the control equation for heat and moisture exchange are defined. The driving force for moisture absorption and desorption is mainly the difference in water vapor partial pressure. The location where moisture absorption and desorption usually occur is in the thin layer on the material surface, and the process of moisture absorption and desorption is accompanied by heat transfer. The thermal boundary conditions include convective heat transfer and water vapor diffusion.

[0196] 4. Establishment of the dynamic heat transfer model of the soil in the underground space and the indoor heat and moisture balance model

[0197] The heat transfer problem between the soil around the underground space is simplified into a two-dimensional unsteady heat conduction problem, see Figure 2 . The boundary conditions of the calculation domain are as Figure 3 shown. The upper boundary is the surface soil boundary, the lower boundary is the isothermal boundary, the side boundaries are adiabatic boundaries, and the boundary in contact with the indoor air is the convective boundary. The indoor heat and moisture balance model includes the energy balance equation and the mass balance equation of the indoor air;

[0198] The finite difference method is used to solve the dynamic heat transfer process of the soil in the underground space.

[0199] 5. Coupled solution of underground dynamic heat transfer and heat and moisture transfer of humidity control materials

[0200] See Figure 3 , the calculation results of the dynamic heat transfer model of the soil in the underground space provide input conditions for the dynamic calculation model of heat and moisture transfer of humidity control materials, and the amount of heat and moisture transfer inside the humidity control material in turn affects the temperature and humidity parameters inside the underground project and the surrounding soil. The specific calculation process is shown in Figure 4 , and the details are as follows:

[0201] (1) At a given time step, calculate the temperature distribution of the soil around the underground space according to the input parameters of the dynamic heat transfer model of the soil in the underground space (i.e., );

[0202] (2) Take the soil temperature distribution as the outer boundary condition of the heat and moisture coupled transfer calculation model of the humidity control material, and combine the indoor temperature and humidity parameters to calculate the heat flux and moisture flux on the inner surface of the humidity control material (i.e., the heat and moisture migration amount inside the humidity control material); specifically, first solve the moisture content of the air in the pores of the material and the temperature inside the material , and then through the heat boundary condition and the moisture boundary condition (both surfaces of the material are convective boundaries), the heat flux and moisture flux can be calculated;

[0203] (3) Substitute the calculated heat flux and moisture flux on the inner surface of the humidity control material into the energy balance equation and mass balance equation of the indoor air respectively, and discretely solve the indoor temperature and humidity parameters (i.e., the air temperature of the external environment where the material is located and the moisture content of the air in the external environment where the material is located ), completing one cycle;

[0204] (4) Start the calculation of the next time step and conduct the next coupling calculation until the set requirements are met, then the calculation ends, and the temperature and humidity changes inside the underground space are output.

[0205] Figure 5 and Figure 6 show the indoor temperature and humidity results of the underground space calculated by this model, and compare them with the measured results, with the average deviation within 3%.

[0206] Example 2

[0207] In an embodiment of the present invention, a computer-readable storage medium is provided, on which a computer program is stored, and when the program is executed by a processor, it implements the prediction method for the passive humidity control effect of urban underground spaces described in any one of Embodiment 1.

[0208] Example 3

[0209] In an embodiment of the present invention, a computer program product is provided, including a computer program / instructions, and when the computer program / instructions are executed by a processor, it implements the prediction method for the passive humidity control effect of urban underground spaces described in any one of Embodiment 1.

[0210] Those skilled in the art should understand that the embodiments of the present application can be provided as methods, systems, or computer program products. Therefore, the present application can take the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk memories, CD-ROMs, optical memories, etc.) containing computer-usable program codes.

[0211] This application is described with reference to the flowcharts and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the present application. It should be understood that each flow and / or block in the flowchart and / or block diagram, and the combination of flows and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing device to generate a machine, such that the instructions executed by the processor of the computer or other programmable data processing device produce means for implementing the functions specified in one or more of the flows Figure 1 one or more of the flows and / or blocks Figure 1 or means for implementing the functions specified in one or more of the blocks.

[0212] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing device to work in a specific manner, such that the instructions stored in the computer-readable memory produce a manufactured article including instruction means that implement the functions specified in one or more of the flows Figure 1 one or more of the flows and / or blocks Figure 1 or means for implementing the functions specified in one or more of the blocks.

[0213] These computer program instructions can also be loaded onto a computer or other programmable data processing device, such that a series of operation steps are executed on the computer or other programmable device to generate a computer-implemented process, and thus the instructions executed on the computer or other programmable device provide steps for implementing the functions specified in one or more of the flows Figure 1 one or more of the flows and / or blocks Figure 1 or means for implementing the functions specified in one or more of the blocks.

[0214] The embodiments of the present invention have been described above in conjunction with the accompanying drawings. However, the present invention is not limited to the above specific embodiments. The above specific embodiments are merely illustrative and not restrictive. Under the inspiration of the present invention, those of ordinary skill in the art can also make many forms without departing from the spirit and scope of the present invention as defined by the claims. These all fall within the protection scope of the present invention.

[0215] The above shows and describes the basic principles, main features, and advantages of the present invention. Those skilled in the art of this industry should understand that the present invention is not limited by the above embodiments. The above embodiments and the descriptions in the specification only illustrate the principles of the present invention. Without departing from the spirit and scope of the present invention, the present invention will have various changes and improvements, and these changes and improvements all fall within the scope of the present invention claimed. The scope of protection claimed by the present invention is defined by the appended claims and their equivalents.

Claims

1. A method for predicting the passive humidity control effect of urban underground space, characterized in that: include: Obtain a pre-built heat and moisture dynamic transfer calculation model for humidity-control materials; Obtain pre-built underground soil dynamic heat transfer calculation models and underground space indoor heat and moisture balance calculation models; The heat and moisture dynamic transfer calculation model of the humidity-controlling material is coupled with the underground soil dynamic heat transfer calculation model and the calculation model of indoor heat and moisture balance in the underground space to predict the passive humidity control effect of the urban underground.

2. The method for predicting the passive humidity control effect of urban underground space according to claim 1, characterized in that: The heat and moisture dynamic transfer calculation model of the humidity-controlling material includes a heat and moisture coupling transfer control equation and a heat and moisture coupling transfer boundary condition; The heat-humidity coupling transfer control equation includes a heat migration control equation and a moisture migration control equation, and the heat migration control equation and the moisture migration control equation establish a heat-humidity coupling transfer relationship through phase change latent heat; The heat and moisture coupling transfer boundary conditions include moisture absorption and desorption boundary conditions on the inner surface of the humidity regulating material and thermal boundary conditions on the inner surface of the humidity regulating material.

3. A method for predicting the passive humidity control effect of urban underground space according to claim 2, characterized in that: The mathematical expression of the thermal migration control equation is: ； The mathematical expression of the wet migration control equation is: ； In the formula, is the moisture content of the air in the pores of the material, is the temperature inside the material, For time; is the density of the material, is the specific heat capacity of the material, is the thermal conductivity of the material, is the latent heat of vaporization; is the moisture storage coefficient related to moisture content, is the temperature-dependent moisture storage coefficient, is the water vapor diffusion coefficient, is the position coordinate; The mathematical expression of the moisture absorption and desorption boundary condition on the inner surface of the humidity regulating material is: ； In the formula, It is the indoor air humidity content of the underground space, that is, the air humidity content of the external environment where the material is located; It is the moisture content of the air in the thin layer on the inner surface of the material; is the convective moisture transfer coefficient; The calculation formula for the thermal boundary condition of the inner surface of the humidity control material is: ； In the formula, is the convective heat transfer coefficient of the inner surface of the material, is the indoor air temperature of the underground space, i.e. the ambient air temperature of the material; It is the temperature of the air in a thin layer on the inner surface of the material.

4. A method for predicting the passive humidity control effect of urban underground space according to claim 3, characterized in that: The moisture storage coefficient related to the moisture content Including moisture storage coefficient of moisture absorption process and moisture storage coefficient of dehumidification process ; The moisture storage coefficient of the moisture absorption process The calculation formula is: ； The moisture storage coefficient of the dehumidification process The calculation formula is: ； In the formula, is the moisture content in the saturated state; It is the minimum relative humidity in the moisture absorption and dehumidification cycle; It is the maximum value of relative humidity during the moisture absorption and dehumidification cycle; It is the isothermal equilibrium moisture content of the moisture absorption curve at a given relative humidity; It is the isothermal equilibrium moisture content of the dehumidification curve at a given relative humidity; is the moisture storage coefficient calculated from the isothermal moisture absorption curve, It is the moisture storage coefficient calculated from the isothermal moisture release curve.

5. A method for predicting the passive humidity control effect of urban underground space according to claim 4, characterized in that: Under different temperature and humidity environments, the moisture storage coefficient of the moisture-adjusting material is continuously updated through the Mualem experimental model. and moisture storage coefficient of dehumidification process .

6. The method for predicting the passive humidity control effect of urban underground space according to claim 1, characterized in that: The underground soil dynamic heat transfer calculation model includes internal heat migration control equations; The calculation model of indoor heat and moisture balance in the underground space includes a heat balance equation and a mass balance equation of indoor air, and boundary conditions of the heat transfer calculation domain of the underground space.

7. A method for predicting the passive humidity control effect of urban underground space according to claim 6, characterized in that: The mathematical expression of the internal heat migration control equation is: ； In the formula, yes The time is at point The soil temperature value, is the density of the soil, is the specific heat capacity of the soil, Yes The thermal conductivity of soil at Yes Thermal conductivity of the material at the location; The mathematical expression of the heat balance equation of the indoor air is: ； The mathematical expression of the mass balance equation of the indoor air is: ； In the formula, is the specific heat of air at constant pressure, is the density of air, is the volume of the room, is the convective heat transfer coefficient, The room wall The surface area of is the indoor air temperature of the underground space, Calculate the room walls The surface temperature, is the heat gain in the underground space. It is the amount of moisture absorbed or released by the humidity-control material attached to the surface of the room enclosure structure; It is the ventilation volume between indoor and outdoor air in underground space. is the humidity content of outdoor air; It is the indoor air humidity content of the underground space, that is, the air humidity content of the external environment where the material is located; is the amount of indoor water vapor generated in the underground space; The boundary conditions of the underground space heat transfer calculation domain include upper boundary conditions, lower boundary conditions, side boundary conditions and convection boundary conditions; The mathematical expression of the upper boundary condition is: ； In the formula, is the shortwave radiation reaching the ground; is the long-wave radiation reaching the ground; It is the long-wave radiation emitted from the ground; It is the convective heat exchange between the earth's surface and the atmosphere; is the shortwave reflectivity of the ground surface; is the long-wave absorption rate of the ground surface; The long-wave emissivity of the ground surface; The mathematical expression of the lower boundary condition is: ； In the formula, is the temperature of the lower boundary of the computational domain, is the annual mean surface temperature; The mathematical expression of the lateral boundary condition is: ； In the formula, is the heat flow at the side boundary of the computational domain; The mathematical expression of the convection boundary condition is: ； In the formula, is the indoor air temperature of the underground space, i.e. the ambient air temperature of the material; is the convective heat transfer coefficient on the inner surface of the material.

8. The method for predicting the passive humidity control effect of urban underground space according to claim 1, characterized in that: The heat and moisture dynamic transfer calculation model of the humidity control material is coupled with the underground soil dynamic heat transfer calculation model and the calculation model of indoor heat and moisture balance in the underground space to predict the passive humidity control effect of the urban underground, including: Continuously execute the preset cycle calculation steps until the set requirements are met, and obtain the final indoor temperature and humidity parameters of the underground space; The cyclic calculation step comprises: Performing differential calculation according to the underground soil dynamic heat transfer calculation model to obtain the internal temperature distribution of the underground soil; The internal temperature distribution of the soil is used as the outer boundary condition of the heat and moisture dynamic transfer calculation model of the humidity regulating material, and combined with the indoor temperature and humidity parameters, the heat flow and moisture flow on the inner surface of the humidity regulating material are calculated according to the heat and moisture dynamic transfer calculation model of the humidity regulating material; The calculated heat flow and moisture flow on the inner surface of the humidity-control material are respectively substituted into the calculation model of indoor heat and moisture balance of the underground space, and the indoor temperature and humidity parameters of the underground space are discretely solved.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the method for predicting the passive humidity control effect of an urban underground space as described in any one of claims 1 to 8 is implemented.

10. A computer program product comprising a computer program / instructions, characterized in that When the computer program / instruction is executed by a processor, the method for predicting the passive humidity control effect of an urban underground space as described in any one of claims 1 to 8 is implemented.

Images