CN106354977A - Method for analyzing indoor thermal environment by using CFD (computational fluid dynamics) model - Google Patents
Method for analyzing indoor thermal environment by using CFD (computational fluid dynamics) model Download PDFInfo
- Publication number
- CN106354977A CN106354977A CN201610867435.5A CN201610867435A CN106354977A CN 106354977 A CN106354977 A CN 106354977A CN 201610867435 A CN201610867435 A CN 201610867435A CN 106354977 A CN106354977 A CN 106354977A
- Authority
- CN
- China
- Prior art keywords
- temperature
- subregion
- indoor
- load
- thermal environment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title abstract description 12
- 239000012530 fluid Substances 0.000 title abstract description 3
- 238000013461 design Methods 0.000 claims abstract description 92
- 238000004378 air conditioning Methods 0.000 claims abstract description 41
- 238000004458 analytical method Methods 0.000 claims description 28
- 238000001816 cooling Methods 0.000 claims description 11
- BTCSSZJGUNDROE-UHFFFAOYSA-N gamma-aminobutyric acid Chemical compound NCCCC(O)=O BTCSSZJGUNDROE-UHFFFAOYSA-N 0.000 claims description 7
- 238000009413 insulation Methods 0.000 claims description 7
- 230000008859 change Effects 0.000 claims description 6
- 238000005192 partition Methods 0.000 claims description 6
- 239000002918 waste heat Substances 0.000 claims description 4
- 238000010438 heat treatment Methods 0.000 claims description 3
- 230000005855 radiation Effects 0.000 claims description 3
- 210000000582 semen Anatomy 0.000 claims description 3
- 238000005286 illumination Methods 0.000 claims description 2
- 230000009471 action Effects 0.000 abstract description 3
- 238000013021 overheating Methods 0.000 abstract 1
- 238000004781 supercooling Methods 0.000 abstract 1
- 238000012937 correction Methods 0.000 description 8
- 230000000694 effects Effects 0.000 description 5
- 238000010276 construction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Air Conditioning Control Device (AREA)
Abstract
The invention discloses a method for analyzing an indoor thermal environment by utilizing a CFD (computational fluid dynamics) model, which comprises the following steps of: s10, designing an initial scheme according to the indoor thermal environment, and acquiring boundary conditions and initial conditions of the indoor thermal environment; s20, calculating an indoor basic temperature field; s30, judging whether the indoor basic temperature field meets the design target temperature, if so, taking the initial scheme of the thermal environment design as the final scheme; if not, go to step S40; s40, calculating each single-source term temperature field by using a CFD model; s50 determining a load to be corrected; s60, according to the load to be corrected, the air conditioning scheme is corrected, so that the indoor basic temperature field meets the design requirement. The method quantitatively determines the heat quantity to be corrected by calculating the temperature distribution under the action of a single heat factor, effectively enables the air-conditioning load to meet the design requirement, and avoids the condition of local area supercooling or overheating in practical application.
Description
Technical field
The present invention relates to Indoor environment thermal environment technical field, it particularly relates to a kind of utilize cfd model to interior
The analysis method of thermal environment.
Technical background
Heavy construction space is increasing, and indoor thermal environment becomes complicated, and the design of air conditioning solutions not only needs to make
Load meets the needs of indoor thermal environment, it is also desirable to make each region meet comfortableness requirement, health after actual motion
Require and horsepower requirements.The design meanses commonly used at present are only using Fluid Mechanics Computation (in literary composition referred to as: cfd) as numeric value analysis
Instrument, boundary condition and initial condition be according to load and design standard obtain although meeting equation of heat balance, but actual empty
Between due to thermal source, equipment equal distribution uneven, subregion easily, to meet density high, and the low feelings of other regional temperatures
Condition, in addition temperature contrast bring pressure differential to cause air flow, the thermal environment being actually formed can be more increasingly complex than expected design.
For heavy construction space, such as: market, movie theatre, arenas, terminal etc., can be to building when carrying out Air-conditioner design
Space carries out region division, and different thermal environment design schemes, the air conditioner load in each region may be taken in different regions
Have differences.Only meet thermally equilibrated design and can not solve the problems, such as rational sharing of load.If by cfd calculate basic
Temperature field is modified to scheme as the foundation of optimizing design scheme, rather than by it as final result it is possible to set
The problem occurring in just will be actually used during meter considers in advance.
Because volumed space building complex structure is various, if necessary to carry out experimentation before building it is necessary to set up
Wind tunnel model, architecture indoor model, so that research process complexity fund input is big, improve design cost.Due to direct basis
The tentative improvement that experimental result is carried out, effect is not necessarily optimum, result precision nor be completely secured, therefore at the beginning of design
The result of prediction and perfect can at utmost minimizing theing carry out for result put into and improve precision, obtain optimization simultaneously
Scheme.
Traditional utilization cfd assists the thinking of process of optimization, is to judge that parameters are unsatisfactory for expected requirement
In the case of repeatedly change model and calculate again, stop calculating when the result of output meets preset value and circulate and export final
Scheme.In these modifications repeatedly and calculating, cfd calculates and completely carries out each time, and the foundation changed each time
There is no specific aim, need to spend longer time.
Content of the invention
Technical problem: the technical problem to be solved is: provide a kind of utilization cfd model to indoor thermal environment
Analysis method, can improve the analysis efficiency of indoor thermal environment, treat modified load simultaneously and carry out quantitative analyses so that Indoor Thermal
The design of environment can meet design requirement through once adjustment.
Technical scheme: for solving above-mentioned technical problem, the technical solution used in the present invention is:
A kind of analysis method to indoor thermal environment for utilization cfd model, this analysis method comprises the following steps:
S10 designs initial scheme according to indoor thermal environment, obtains boundary condition and the initial condition of indoor thermal environment；
S20, by initial condition and boundary condition input cfd model, calculates indoor cardinal temperature field；
S30 judges whether indoor cardinal temperature field meets design object temperature, if it is satisfied, then thermal environment design is initial
Scheme is as final scheme；If be unsatisfactory for, enter step s40；
S40 utilizes cfd model to calculate each single source item temperature field；
S50, according to each single source item temperature field, determines load to be modified；
S60, according to load to be modified, is modified to air conditioning solutions so that indoor cardinal temperature field meets design
Require.
As preference: described step s40 specifically includes:
S401 judges whether indoor division subregion air conditioning system, if it is, entering step s402；If it is not, then
Enter step s403；
S402 sets interior n subregion air conditioning system, to either partition domain m, calculates each subregion using cfd model
Temperature t that load is formed in this subregion mi-m, i=1,2, n, m is integer, and 2≤m≤n；
Indoor temperature initial condition is set to design object temperature by s403, using single source or single low-temperature receiver as source item,
The boundary condition of source item is inputted cfd model, it is zero that its waste heat factor is set to thermal insulation or caloric value, counts respectively using cfd model
Calculate each source item temperature field.
As preference: in described step s40, hot factor includes thermal source and low-temperature receiver；Thermal source includes human body, illuminates, sets
Standby, building enclosure and heat supply end；Low-temperature receiver includes air-supply and cooling end.
As preference: in described s402, calculate the temperature that each subregion load formed in this subregion m and specifically wrap
Include:
Indoor temperature initial condition is set to design object temperature, successively using hot factor all of in each subregion as this
A subregional source item, the boundary condition of the source item in this region is inputted cfd model, and remaining subregional hot factor is set to thermal insulation
Or caloric value is zero, calculate, using cfd model, the temperature that this subregion load is formed in subregion m respectively.
As preference: described s50 specifically includes:
S501 judges whether indoor division subregion air conditioning system, if it is, entering step s502；If it is not, then
Enter step s503；
S502 sets the design object temperature of subregion m as tm-0, calculate temperature t in this subregion m for each subregion loadi-m
Temperature difference δ t with this subregion design object temperaturei-m, as shown in formula (1):
δti-m=ti-m-tm-0Formula (1)
Calculate the heat of the heat inflow subregion m of each subregion air-conditioning cooling or heat supply, as shown in formula (2):
qi-m=cmmδti-mFormula (2)
Wherein: qi-mThe heat of the heat inflow subregion m of expression subregion i air-conditioning cooling or heat supply, i=1,
2nd, n, m≤n, c represent air specific heat, mmRepresent the air quality of subregion m, work as mmDuring for positive number, represent heat stream
Enter, work as mmDuring for negative, represent that cold flows into；
M load δ q to be modified in subregion is determined according to formula (3)m:
S503 sets indoor design target temperature as t0, by step 20) and the cardinal temperature field and the design object temperature that obtain enter
Row compares, and finds out the region r being unsatisfactory for design object temperature in cardinal temperature field, then calculates each hot factor i' respectively in region
Temperature t that r is formedI', r, i'=1,2 ..., s, r=1,2, w, wherein, s is indoor hot factor sum, and w is basic
It is unsatisfactory for the region sum of design object temperature in temperature field；
Calculate temperature and design object temperature t that hot factor i' is formed in region r according to formula (4)0The temperature difference:
δti',r=ti'.r-t0Formula (4)
Calculate r load to be modified in region according to formula (5):
Wherein, c represents air specific heat, and ρ represents atmospheric density, vrRepresent the volume of region r.
As preference: in described s60, air conditioning solutions are modified specifically including: change tuyere position, wind
One of mouth quantity, wind pushing temperature, radiant panel area pattern or combination.
As preference: in described s10, boundary condition includes building enclosure parameter, air output, wind pushing temperature, air-supply
Speed, return air amount, return air temperature, the setting temperature of the semen donors of radiant panel, the heating load of radiant panel or cold emission plate；Initial strip
Part includes building enclosure initial temperature, indoor initial temperature, indoor initial humidity.
Beneficial effect: compared with prior art, the embodiment of the present invention has the advantages that analysis method of the present invention can
To improve the analysis efficiency of indoor thermal environment, treat modified load simultaneously and carry out quantitative analyses so that the design of indoor thermal environment
Scheme can meet design requirement through once adjustment.In prior art, when judge indoor cardinal temperature field be unsatisfactory for design
During target temperature, will rule of thumb method, initial condition and boundary condition are modified, or recycle cfd calculate indoor
Temperature field, until it meets design object temperature.Compared with prior art, the method for the embodiment of the present invention, is judging interior
Cardinal temperature field is unsatisfactory for during design object temperature calculating it is not necessary to recycle cfd, but each list being calculated using cfd
Temperature field under one hot factor effect, and determine load to be modified.With load value to be modified, easily realize to air adjustment
The correction of scheme is so that it meets design requirement.The embodiment of the present invention do not need to recycle cfd calculated it is not required that
Carry out parameters revision using empirical value, greatly reduce amount of calculation, improve analysis efficiency.Meanwhile, to load to be modified it is
Air conditioning solutions are more accurately revised by quantitative analyses, improve analysis precision.
Brief description
In order that the object, technical solutions and advantages of the present invention are clearer, the following drawings can be made further to the present invention
Describe in detail, wherein:
Fig. 1 is the FB(flow block) of the present invention.
Fig. 2 is the indoor temperature distribution figure under the base regime of the embodiment of the present invention 1；
Fig. 3 (a) is the indoor temperature field of room lighting independent role in the embodiment of the present invention 1；
Fig. 3 (b) is the indoor temperature field of internal heat resource independent role in the embodiment of the present invention 1；
Fig. 3 (c) is the indoor temperature field of fresh air outlet independent role in the embodiment of the present invention 1；
Fig. 3 (d) is the indoor temperature field of cold emission plate independent role in the embodiment of the present invention 1.
Fig. 4 is the indoor temperature distribution figure after the embodiment of the present invention 1 changes cold emission Board position；
Fig. 5 is the temperature profile under the base regime of the embodiment of the present invention 2；
Fig. 6 (a) is cold emission plate operative temperature field in subregion 1 in the embodiment of the present invention 2；
Fig. 6 (b) is Source temperature field in subregion 1 in the embodiment of the present invention 2；
Fig. 6 (c) is fresh air operative temperature field in subregion 1 in the embodiment of the present invention 2；
Fig. 7 (a) is subregion 2 cold emission plate operative temperature field in the embodiment of the present invention 2；
Fig. 7 (b) is subregion 2 heat transfer across wall operative temperature field in the embodiment of the present invention 2；
Fig. 7 (c) is subregion 2 Source temperature field in the embodiment of the present invention 2；
Fig. 8 is the indoor temperature distribution figure after the embodiment of the present invention 2 Adjusted Option.
Specific embodiment
Below in conjunction with the accompanying drawings, technical scheme is described in detail.
A kind of analysis method to indoor thermal environment for the utilization cfd model provided in an embodiment of the present invention, comprises the following steps:
S10 designs initial scheme according to indoor thermal environment, obtains boundary condition and the initial condition of indoor thermal environment.
Step s10 specifically includes: determines indoor thermal environment design initial scheme according to Architectural Structure Design drawing；According to heat
Environment Design initial scheme design air adjusts form, obtains base regime；Set the side of indoor thermal environment according to base regime
Boundary's condition and initial condition.Wherein, boundary condition includes building enclosure parameter, air output, wind pushing temperature, air supply velocity, return air
Amount, return air temperature, the setting temperature of the semen donors of radiant panel, the heating load of radiant panel or cold emission plate；Initial condition includes enclosing
Protection structure initial temperature, indoor initial temperature, indoor initial humidity.
S20, by initial condition and boundary condition input cfd model, calculates indoor cardinal temperature field.
S30 judges whether indoor cardinal temperature field meets design object temperature, if it is satisfied, then thermal environment design is initial
Scheme is as final scheme；If be unsatisfactory for, that is, there is the region that temperature is higher or lower than design object temperature in the interior space, then
Enter step s40.
S40 utilizes cfd model to calculate each single source item temperature field.
S50, according to each single source item temperature field, determines load to be modified.
S60, according to load to be modified, is modified to air conditioning solutions so that indoor cardinal temperature field meets design
Require.Air conditioning solutions are modified specifically including: change tuyere position, inlet number, wind pushing temperature, radiation plate suqare
One of mode or combination.
In the method for above-described embodiment, by inputting initial condition and boundary condition in cfd model, calculate indoor basic
Temperature field, then judges whether indoor cardinal temperature field meets design object temperature.Calculate each single source item temperature using cfd model
Degree field, for ungratified regional analysises.According to each single source item temperature field, determine the load to be modified being unsatisfactory for region.
This load to be modified obtains determination value through calculating.So, according to load to be modified, easily air conditioning solutions are carried out
Revise, so that indoor cardinal temperature field meets design requirement.
In prior art, when judging that indoor cardinal temperature field is unsatisfactory for design object temperature, will rule of thumb method,
Initial condition and boundary condition are modified, or recycle cfd calculating indoor temperature field, until it meets design object temperature
Degree.Compared with prior art, the method for the embodiment of the present invention, is judging that indoor cardinal temperature field is unsatisfactory for design object temperature
When calculate it is not necessary to recycle cfd, but the temperature field under each single hot factor effect being calculated using cfd, and determining
Load to be modified.With load value to be modified, easily realize the correction to air conditioning solutions so that its satisfaction design will
Ask.The embodiment of the present invention does not need to recycle cfd and is calculated it is not required that carrying out parameters revision using empirical value, significantly
Decrease amount of calculation, improve analysis efficiency.Meanwhile, to load to be modified be quantitative analyses, air conditioning solutions are carried out
More accurately revise, improve analysis precision.
In above-described embodiment, indoor thermal environment can be divided into two kinds of situations, and one kind is by indoor division subregion air-conditioning system
The environment of system, another kind is the environment not dividing subregion air conditioning system.Below based on two kinds of varying environments, specifically divided
Analysis.
As preference, described step s40 specifically includes s401 s403:
S401 judges whether indoor division subregion air conditioning system, if it is, entering step s402；If it is not, then
Enter step s403.
S402 sets interior n subregion air conditioning system, to either partition domain m, calculates each subregion using cfd model
Temperature t that load is formed in this subregioni-m, i=1,2, n, m is integer, and 2≤m≤n.
In s402, calculate the temperature that each subregion load formed in this subregion and specifically include:
Indoor temperature initial condition is set to design object temperature, successively using hot factor all of in each subregion as this
A subregional source item, the boundary condition of the source item in this region is inputted cfd model, and remaining subregional hot factor is set to thermal insulation
Or caloric value is zero, calculate, using cfd model, the temperature that this subregion load is formed in subregion m respectively.
Indoor temperature initial condition is set to design object temperature by s403, using single source or single low-temperature receiver as source item,
The boundary condition of source item is inputted cfd model, it is zero that its waste heat factor is set to thermal insulation or caloric value, counts respectively using cfd model
Calculate each source item temperature field.Wherein, hot factor includes thermal source and low-temperature receiver；Thermal source includes human body, illumination, equipment, building enclosure and confession
Hot end.Low-temperature receiver includes air-supply and cooling end.
As preference, described s50 specifically includes s501 s503:
S501 judges whether indoor division subregion air conditioning system, if it is, entering step s502；If it is not, then
Enter step s503.
S502 sets the design object temperature of subregion m as tm-0, calculate temperature t in this subregion m for each subregion loadi-m
Temperature difference δ t with this subregion design object temperaturei-m, as shown in formula (1):
δti-m=ti-m-tm-0Formula (1)
Calculate the heat of the heat inflow subregion m of each subregion air-conditioning cooling or heat supply, as shown in formula (2):
qi-m=cmmδti-mFormula (2)
Wherein: qi-mThe heat of the heat inflow subregion m of expression subregion i air-conditioning cooling or heat supply, i=1,
2nd, n, m≤n, c represent air specific heat, mmRepresent the air quality of subregion m, work as mmDuring for positive number, represent heat stream
Enter, work as mmDuring for negative, represent that cold flows into；
M load δ q to be modified in subregion is determined according to formula (3)m:
S503 sets indoor design target temperature as t0, by step 20) and the cardinal temperature field and the design object temperature that obtain enter
Row compares, and finds out the region r being unsatisfactory for design object temperature in cardinal temperature field, then calculates each hot factor i' respectively in region
Temperature t that r is formedI', r, i'=1,2 ..., s, r=1,2 ..., w, wherein, s is indoor hot factor sum, and w is cardinal temperature
It is unsatisfactory for the region sum of design object temperature in.
Calculate temperature and design object temperature t that hot factor i' is formed in region r according to formula (4)0The temperature difference:
δti',r=ti'.r-t0Formula (4)
Calculate r load to be modified in region according to formula (5):
Load to be modified is determined by step s50.After determining load value to be modified, easily revised according to this value
Air conditioning solutions are so that indoor cardinal temperature field meets design requirement.The embodiment of the present invention is passed through to determine load to be modified,
Make air conditioning solutions correction more accurate.
In the embodiment of the present invention, in first time base regime calculates, all of boundary condition and initial condition be all by
According to the facts border design conditions are carried out.This is same as the prior art.S40 only calculates single hot factor operative temperature field.This border
The calculating of condition takes less, and speed is fast.The temperature formation contribution degree that nondimensionalization single factors operative temperature is distributed gained is permissible
Predict the distribution meeting, thus correction air regulation scheme.
It is adaptable to various air conditioning system, to building type, air-conditioning shape does not limit the method for the embodiment of the present invention.
Air conditioning terminal form can be air-conditioning system or traditional air-conditioning system.
Example is for two example two below.
Embodiment 1
If not dividing subregion air conditioning system in indoor thermal environment.
If indoor design target temperature is 26 DEG C.Initial scheme is designed according to indoor thermal environment, by initial condition and border
In condition entry cfd model, calculate indoor cardinal temperature field, result is as shown in Figure 2.Figure it is seen that indoor thermal environment sets
In the cardinal temperature field that meter initial scheme is formed, subregion temperature is higher than design object temperature, and subregion temperature is less than and sets
Meter target temperature.Figure it is seen that the heat of the generation such as internal heat resource, luminaire, building enclosure has more flowed to room
Interior upper right quarter space, and indoor lower left quarter space obtains more colds.
Using single source or single low-temperature receiver as source item, the boundary condition of source item is inputted cfd model, its waste heat factor sets
It is zero for adiabatic or caloric value, calculate each source item temperature field respectively using cfd model.In the present embodiment, source item includes indoor photograph
Bright, internal heat resource, air outlet, cold emission plate.
As shown in figure 3, calculate the same sectional view in each source item temperature field respectively using cfd model, i.e. different hot factor lists
The solely same sectional view of the indoor temperature field that effect is formed.
From figure 3, it can be seen that the indoor temperature field that different hot factor independent roles are formed is different, there is the temperature difference.Its
In, the temperature profile that Fig. 3 (a) is formed for room lighting.The temperature profile that Fig. 3 (b) is formed for internal heat resource.Fig. 3 (c) is
The temperature profile that air outlet is formed.Fig. 3 is the temperature profile that (d) cold emission plate shape becomes.The cold emission plate that Fig. 3 (d) reflects
Cooling-down effect makes target area temperature be less than 1～2 DEG C of design object temperature.
The cardinal temperature field of acquisition and design object temperature are to compare, find out in cardinal temperature field and be unsatisfactory for design object
The region r of temperature, then calculates temperature t that each hot factor i' is formed in region r respectivelyI', r, calculate hot factor according to formula (4)
The temperature that i' is formed in region r and design object temperature t0The temperature difference:
δti',r=ti'.r-t0Formula (4)
Calculate r load to be modified in region according to formula (5):
Because the air outlet temperature that space is formed indoors reaches the level close with target temperature in upper right quarter, not many
Remaining cooling capacity, it is contemplated that change cold emission Board position.During the present embodiment original state, cold emission plate is located at lower section.Existing
To move on cold emission plate.
After correction, indoor temperature distribution figure is as shown in Figure 4.As can be seen from Figure 4: indoor temperature is 26 ± 0.5 DEG C, meet
Design requirement.
Embodiment 2
If indoor thermal environment divides subregion air conditioning system, as shown in figure 5, being divided into subregion 1 and 2 two points of subregion altogether
Area.
If indoor design target temperature is 26 DEG C.Initial scheme is designed according to indoor thermal environment, by initial condition and border
In condition entry cfd model, calculate indoor cardinal temperature field, result is as shown in Figure 5.Sitting posture zone of action is below 1.3m, stands
Appearance zone of action is below 2m.From Fig. 5 cardinal temperature field situation, determine region (level height 0.2m carrying out load correction
1.3m) it is located in subregion 1.
In subregion air-conditioning, as a hot factor, all thermals source in this region and low-temperature receiver all arrange reality in a region
Border boundary condition, whole region is a hot factor.Successively using hot factor all of in each subregion as the one of respective partition
Source item, the boundary condition of the source item in this region is inputted cfd model, remaining subregional hot factor is set to thermal insulation or caloric value is
Zero, calculate, using cfd model, the temperature that this partition load is formed in subregion 1 respectively.In the present embodiment, subregion 1 is supplied to this area
The cold in domain includes: cold emission plate and fresh air；The load that subregion 1 produces to this region includes: thermal source.Subregion 2 is supplied to this area
The cold in domain includes: cold emission plate；The load that subregion 2 is supplied to this region includes: building enclosure.
Fig. 6 (a) is subregion 1 cold emission plate operative temperature field.Fig. 6 (b) is subregion 1 Source temperature field.Fig. 6 (c) is
Subregion 1 fresh air operative temperature field.Fig. 7 (a) is subregion 2 cold emission plate operative temperature field.Fig. 7 (b) is subregion 2 heat transfer across wall
Operative temperature field.Fig. 7 (c) is subregion 2 Source temperature field.
Calculate temperature t in subregion 1 for each partition loadi-mTemperature difference δ t with this split plot design target temperaturei-m.In this reality
Apply in example, subregion 1 cold emission plate is -1.4 DEG C with the temperature difference of subregion 1 design object temperature, and subregion 1 fresh air designs mesh with subregion 1
The temperature difference of mark temperature is -1.1 DEG C, and subregion 1 internal heat resource is+1.6 DEG C with the temperature difference of subregion 1 design object temperature.The cold spoke of subregion 2
The temperature difference penetrating plate with subregion 1 design object temperature is -0.7 DEG C, the temperature difference of subregion 2 luminaire and subregion 1 design object temperature
For+0.2 DEG C, subregion 2 building enclosure is+0.1 DEG C with the temperature difference of subregion 1 design object temperature.Wherein, the temperature difference is on the occasion of expression
Temperature is higher than subregion 1 design object temperature；The temperature difference is negative value, represents that temperature is less than subregion 1 design object temperature.
δt1-1=-1.4-1.1+1.6=-0.9 DEG C
δt2-1=-0.7+0.2+0.1=-0.4 DEG C
Calculate the heat of the heat inflow subregion 1 of each subregion air-conditioning cooling or heat supply, as shown in formula (2):
qi-m=cmmδti-mFormula (2)
In the present embodiment,
Subregion 1 load δ q to be modified is determined according to formula (3)m:
In the present embodiment,
Load correction area temperature is less than suitable temperature, illustrates that the cold that there is extra 0.057kw in this region flows into this region.If
Carry out checking correction to this region, then need revise amount be: reduce 0.057kw cold flow into.This example is to illustrate to use this data
Carry out project setting, but be not limited only to the method: improve wind pushing temperature, 0.057kw.This is equivalent to 1.3 DEG C of wind pushing temperature of raising.
After Fig. 8 is raising wind pushing temperature, using the indoor cardinal temperature field result of cfd computed in software.Can from Fig. 8
Go out: the indoor temperature after regulation is 26 ± 0.5 DEG C, meets design requirement.
The basic principles, principal features and advantages of the present invention have been shown and described above.Those skilled in the art should
Understand, the present invention do not limited by above-mentioned specific embodiment, the description in above-mentioned specific embodiment and description be intended merely to into
One step illustrate the present invention principle, without departing from the spirit and scope of the present invention, the present invention also have various change and
Improve, these changes and improvements both fall within scope of the claimed invention.The scope of protection of present invention will by right
Ask book and its equivalent thereof.
Claims (7)
1. a kind of utilization cfd model to the analysis method of indoor thermal environment it is characterised in that: this analysis method includes following step
Rapid:
S10 designs initial scheme according to indoor thermal environment, obtains boundary condition and the initial condition of indoor thermal environment；
S20, by initial condition and boundary condition input cfd model, calculates indoor cardinal temperature field；
S30 judges whether indoor cardinal temperature field meets design object temperature, if it is satisfied, then by thermal environment design initial scheme
As final scheme；If be unsatisfactory for, enter step s40；
S40 utilizes cfd model to calculate each single source item temperature field；
S50, according to each single source item temperature field, determines load to be modified；
S60 according to load to be modified, air conditioning solutions are modified so that indoor cardinal temperature field meet design will
Ask.
2. according to the utilization cfd model described in claim 1 to the analysis method of indoor thermal environment it is characterised in that: described
Step s40 specifically includes:
S401 judges whether indoor division subregion air conditioning system, if it is, entering step s402；If it is not, then entering
Step s403；
S402 sets interior n subregion air conditioning system, to either partition domain m, calculates each subregion load using cfd model
Temperature t being formed in this subregion mi-m, i=1,2, n, m is integer, and 2≤m≤n；
Indoor temperature initial condition is set to design object temperature by s403, using single source or single low-temperature receiver as source item, by source
The boundary condition input cfd model of item, it is zero that its waste heat factor is set to thermal insulation or caloric value, is calculated respectively respectively using cfd model
Source item temperature field.
3. according to the utilization cfd model described in claim 2 to the analysis method of indoor thermal environment it is characterised in that: described
In step s40, hot factor includes thermal source and low-temperature receiver；Thermal source includes human body, illumination, equipment, building enclosure and heat supply end；Low-temperature receiver
Including air-supply and cooling end.
4. according to the utilization cfd model described in claim 2 to the analysis method of indoor thermal environment it is characterised in that: described
In s402, calculate the temperature that each subregion load formed in this subregion m and specifically include:
Indoor temperature initial condition is set to design object temperature, successively using hot factor all of in each subregion as this subregion
One source item in domain, the boundary condition of the source item in this region is inputted cfd model, and remaining subregional hot factor is set to thermal insulation or sends out
Heat is zero, calculates, using cfd model, the temperature that this subregion load is formed in subregion m respectively.
5. according to the utilization cfd model described in claim 2 to the analysis method of indoor thermal environment it is characterised in that: described
S50 specifically includes:
S501 judges whether indoor division subregion air conditioning system, if it is, entering step s502；If it is not, then entering
Step s503；
S502 sets the design object temperature of subregion m as tm-0, calculate temperature t in this subregion m for each subregion loadi-mWith this
The temperature difference δ t of subregion design object temperaturei-m, as shown in formula (1):
δti-m=ti-m-tm-0Formula (1)
Calculate the heat of the heat inflow subregion m of each subregion air-conditioning cooling or heat supply, as shown in formula (2):
qi-m=cmmδti-mFormula (2)
Wherein: qi-mRepresent that the heat of subregion i air-conditioning cooling or heat supply flows into the heat of subregion m, i=1,2 ..., n, m≤
N, c represent air specific heat, mmRepresent the air quality of subregion m, work as mmDuring for positive number, represent that heat flows into, work as mmFor negative
When, represent that cold flows into；
M load δ q to be modified in subregion is determined according to formula (3)m:
S503 sets indoor design target temperature as t0, by step 20) and the cardinal temperature field and the design object temperature that obtain compared
Relatively, find out the region r being unsatisfactory for design object temperature in cardinal temperature field, then calculate each hot factor i' respectively in region r shape
Temperature t becomingI', r, i'=1,2 ..., s, r=1,2, w, wherein, s is indoor hot factor sum, and w is temperature substantially
It is unsatisfactory for the region sum of design object temperature in degree field；
Calculate temperature and design object temperature t that hot factor i' is formed in region r according to formula (4)0The temperature difference:
δti',r=ti'.r-t0Formula (4)
Calculate r load to be modified in region according to formula (5):
Wherein, c represents air specific heat, and ρ represents atmospheric density, vrRepresent the volume of region r.
6. according to the utilization cfd model described in claim 1 to the analysis method of indoor thermal environment it is characterised in that: described
In s60, air conditioning solutions are modified specifically including: change tuyere position, inlet number, wind pushing temperature, radiation plate face
One of long-pending mode or combination.
7. according to the utilization cfd model described in claim 1 to the analysis method of indoor thermal environment it is characterised in that: described
In s10, boundary condition includes building enclosure parameter, air output, wind pushing temperature, air supply velocity, return air amount, return air temperature, radiation
The setting temperature of the semen donors of plate, the heating load of radiant panel or cold emission plate；Initial condition includes building enclosure initial temperature, room
Interior initial temperature, indoor initial humidity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610867435.5A CN106354977B (en) | 2016-09-29 | 2016-09-29 | Method for analyzing indoor thermal environment by using CFD (computational fluid dynamics) model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610867435.5A CN106354977B (en) | 2016-09-29 | 2016-09-29 | Method for analyzing indoor thermal environment by using CFD (computational fluid dynamics) model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106354977A true CN106354977A (en) | 2017-01-25 |
| CN106354977B CN106354977B (en) | 2019-08-06 |
Family
ID=57865800
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610867435.5A Active CN106354977B (en) | 2016-09-29 | 2016-09-29 | Method for analyzing indoor thermal environment by using CFD (computational fluid dynamics) model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106354977B (en) |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107025334A (en) * | 2017-03-10 | 2017-08-08 | 国网吉林省电力有限公司 | Central heating system heat user integrated dynamic model method for building up and device |
| CN107839581A (en) * | 2017-09-05 | 2018-03-27 | 鲁东大学 | A kind of automatic control feedback compartment and method based on temperature field three-dimensional imaging |
| CN108981097A (en) * | 2018-06-19 | 2018-12-11 | 广东美的暖通设备有限公司 | Air-conditioning system and its method of adjustment and device can be needed |
| CN110826221A (en) * | 2019-11-04 | 2020-02-21 | 北京工业大学 | Method for decoupling indoor temperature field of building |
| CN111426034A (en) * | 2020-05-15 | 2020-07-17 | 北京草木元人工环境技术有限公司 | Indoor air conditioner ventilation air flow organization control system |
| CN111709199A (en) * | 2020-06-28 | 2020-09-25 | 北京工业大学 | Method for measuring heat productivity of equipment |
| CN112464408A (en) * | 2020-11-27 | 2021-03-09 | 四川长虹空调有限公司 | Simulation evaluation method for uniformity of air volume and room temperature field of through-flow air duct air conditioner on-hook |
| CN112906311A (en) * | 2019-11-19 | 2021-06-04 | 矢崎能源系统公司 | Thermal load calculation device |
| CN112923504A (en) * | 2020-12-31 | 2021-06-08 | 广州市昊铭机电科技有限公司 | Classroom environment control device and method |
| CN114110970A (en) * | 2020-08-31 | 2022-03-01 | 北京小米移动软件有限公司 | Indoor temperature control method and system for air conditioner |
| CN115164352A (en) * | 2022-05-30 | 2022-10-11 | 中建八局第三建设有限公司 | Large-space air conditioner whole-process debugging method |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102930145A (en) * | 2012-10-18 | 2013-02-13 | 中建三局第二建设工程有限责任公司 | Courtyard indoor thermal environment analyzing method based on CFD (Computational Fluid Dynamics) information technology |
| CN103455676A (en) * | 2013-09-04 | 2013-12-18 | 重庆建工住宅建设有限公司 | Method for simulating indoor thermal environment by fluid mechanics |
| CN103884075A (en) * | 2014-01-06 | 2014-06-25 | 浙江工业大学 | Computational fluid dynamics and energy prediction hybrid based greenhouse energy-saving control method |
-
2016
- 2016-09-29 CN CN201610867435.5A patent/CN106354977B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102930145A (en) * | 2012-10-18 | 2013-02-13 | 中建三局第二建设工程有限责任公司 | Courtyard indoor thermal environment analyzing method based on CFD (Computational Fluid Dynamics) information technology |
| CN103455676A (en) * | 2013-09-04 | 2013-12-18 | 重庆建工住宅建设有限公司 | Method for simulating indoor thermal environment by fluid mechanics |
| CN103884075A (en) * | 2014-01-06 | 2014-06-25 | 浙江工业大学 | Computational fluid dynamics and energy prediction hybrid based greenhouse energy-saving control method |
Cited By (16)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107025334A (en) * | 2017-03-10 | 2017-08-08 | 国网吉林省电力有限公司 | Central heating system heat user integrated dynamic model method for building up and device |
| CN107025334B (en) * | 2017-03-10 | 2020-05-12 | 国网吉林省电力有限公司 | Method and device for establishing overall dynamic model of heat users of central heating system |
| CN107839581A (en) * | 2017-09-05 | 2018-03-27 | 鲁东大学 | A kind of automatic control feedback compartment and method based on temperature field three-dimensional imaging |
| CN108981097A (en) * | 2018-06-19 | 2018-12-11 | 广东美的暖通设备有限公司 | Air-conditioning system and its method of adjustment and device can be needed |
| CN110826221A (en) * | 2019-11-04 | 2020-02-21 | 北京工业大学 | Method for decoupling indoor temperature field of building |
| CN112906311A (en) * | 2019-11-19 | 2021-06-04 | 矢崎能源系统公司 | Thermal load calculation device |
| CN112906311B (en) * | 2019-11-19 | 2024-04-23 | 矢崎能源系统公司 | Thermal load calculation device |
| CN111426034A (en) * | 2020-05-15 | 2020-07-17 | 北京草木元人工环境技术有限公司 | Indoor air conditioner ventilation air flow organization control system |
| CN111426034B (en) * | 2020-05-15 | 2021-06-15 | 北京草木元人工环境技术有限公司 | Indoor air conditioner ventilation air flow organization control system |
| CN111709199A (en) * | 2020-06-28 | 2020-09-25 | 北京工业大学 | Method for measuring heat productivity of equipment |
| CN111709199B (en) * | 2020-06-28 | 2024-01-26 | 北京工业大学 | Method for measuring heating value of equipment |
| CN114110970A (en) * | 2020-08-31 | 2022-03-01 | 北京小米移动软件有限公司 | Indoor temperature control method and system for air conditioner |
| CN112464408A (en) * | 2020-11-27 | 2021-03-09 | 四川长虹空调有限公司 | Simulation evaluation method for uniformity of air volume and room temperature field of through-flow air duct air conditioner on-hook |
| CN112923504A (en) * | 2020-12-31 | 2021-06-08 | 广州市昊铭机电科技有限公司 | Classroom environment control device and method |
| CN115164352A (en) * | 2022-05-30 | 2022-10-11 | 中建八局第三建设有限公司 | Large-space air conditioner whole-process debugging method |
| CN115164352B (en) * | 2022-05-30 | 2024-05-28 | 中建八局第三建设有限公司 | Whole-process debugging method for large-space air conditioner |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106354977B (en) | 2019-08-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106354977A (en) | Method for analyzing indoor thermal environment by using CFD (computational fluid dynamics) model | |
| Shinoda et al. | A review of the surface heat transfer coefficients of radiant heating and cooling systems | |
| Acikgoz et al. | Machine learning approach to predict the heat transfer coefficients pertaining to a radiant cooling system coupled with mixed and forced convection | |
| Wemhoff et al. | Predictions of energy savings in HVAC systems by lumped models | |
| Zuo et al. | Coupling indoor airflow, HVAC, control and building envelope heat transfer in the Modelica Buildings library | |
| Mao et al. | Comparative studies on using RSM and TOPSIS methods to optimize residential air conditioning systems | |
| CN103049612B (en) | Building indoor environment optimization method based on model order reduction technology | |
| Ito et al. | CFD benchmark tests for indoor environmental problems: Part 3 numerical thermal manikins | |
| Myhren et al. | Improving the thermal performance of ventilation radiators–The role of internal convection fins | |
| Li et al. | A new explicit equation for accurate friction factor calculation of smooth pipes | |
| CN102930145A (en) | Courtyard indoor thermal environment analyzing method based on CFD (Computational Fluid Dynamics) information technology | |
| Silva-Llanca et al. | Determining wasted energy in the airside of a perimeter-cooled data center via direct computation of the Exergy Destruction | |
| Acikgoz et al. | Experimental and numerical investigation of the correlation between radiative and convective heat-transfer coefficients at the cooled wall of a real-sized room | |
| Acikgoz | A novel evaluation regarding the influence of surface emissivity on radiative and total heat transfer coefficients in radiant heating systems by means of theoretical and numerical methods | |
| Liang et al. | Cooling load for the design of air terminals in a general non-uniform indoor environment oriented to local requirements | |
| Chen et al. | Energy analysis of three ventilation systems for a large machining plant | |
| Li et al. | Study on the impact of parallel jet spacing on the performance of multi-jet stratum ventilation | |
| Zhao et al. | Room airborne pollutant separation by the use of air curtains in the large building enclosure: Infiltration efficiency and partial enclosure ventilation rate | |
| CN110826221A (en) | Method for decoupling indoor temperature field of building | |
| Mao et al. | A simplified numerical study on the energy performance and thermal environment of a bedroom TAC system | |
| Camci et al. | An experimental study on the heat transfer characteristics over a radiant cooled wall exposed to mixed and forced convection driven by displacement ventilation | |
| Zhang et al. | Review and development of the contribution ratio of indoor climate (CRI) | |
| Pacak et al. | Application of computational fluid dynamics (CFD) methods to analyze energy efficiency of indirect evaporative coolers | |
| Zhang et al. | Cooling effect of air movement on heating performances of advanced air distribution | |
| Yin et al. | Using a Mass and Energy Balance Approach to Model the Performance of Parallel Fan-Powered Terminal Units with Fixed-Airflow Fans. |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |