CN106354977B - It is a kind of using CFD model to the analysis method of indoor thermal environment - Google Patents
It is a kind of using CFD model to the analysis method of indoor thermal environment Download PDFInfo
- Publication number
- CN106354977B CN106354977B CN201610867435.5A CN201610867435A CN106354977B CN 106354977 B CN106354977 B CN 106354977B CN 201610867435 A CN201610867435 A CN 201610867435A CN 106354977 B CN106354977 B CN 106354977B
- Authority
- CN
- China
- Prior art keywords
- temperature
- subregion
- indoor
- air
- 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.)
- Active
Links
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 kind of using CFD model to the analysis method of indoor thermal environment, comprising: S10 designs initial scheme according to indoor thermal environment, obtains the boundary condition and primary condition of indoor thermal environment；S20 calculates indoor cardinal temperature field；S30 judges whether indoor cardinal temperature field meets design object temperature, if it is satisfied, then using thermal environment design initial scheme as final scheme；If conditions are not met, then entering step S40；S40 calculates each single source item temperature field using CFD model；S50 determines load to be modified；S60 is modified air conditioning solutions according to load to be modified, so that indoor cardinal temperature field meets design requirement.The method of the present invention quantitatively determines by calculating the Temperature Distribution under single hot factor effect and needs modified heat, effectively air conditioner load is made to reach design requirement, avoid occurring the case where regional area supercooling or overheat in practical application.
Description
Technical field
The present invention relates to Indoor environment thermal environment technical fields, it particularly relates to which 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 which each region is made to meet comfort requirement, health after actual motion
It is required that and horsepower requirements.Currently used design means only regard Fluid Mechanics Computation (text in referred to as: CFD) as numeric value analysis
Tool, boundary condition and primary condition are obtained according to load and design standard, practical empty although meeting equation of heat balance
Between due to heat source, equipment equal distribution it is uneven, easily there is partial region and meet density height, and the feelings that other regional temperatures are relatively low
Condition, temperature difference brings pressure difference that air is caused to flow in addition, and 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, it can be to building when carrying out Air-conditioner design
Space carries out region division, and different regions may take different thermal environment design schemes, the air conditioner load in each region
It has differences.Only meeting thermally equilibrated design not can solve reasonable sharing of load problem.If CFD is calculated basic
Temperature field is modified scheme as the foundation of optimizing design scheme, rather than using it as final result, so that it may it is setting
The problem of will appear in just actually using during meter, considers in advance.
Due to volumed space building structure is complicated multiplicity, if necessary to carry out experimental study before construction, it is necessary to establish
Wind tunnel model, architecture indoor model improve design cost so that research process complexity capital investment is big.Due to direct basis
The tentative improvement that experimental result carries out, effect is not necessarily optimal, and as a result precision can not be completely secured, therefore at the beginning of design
The result of prediction and it can utmostly reduce investment for carry out improve of result and improve precision, while being optimized
Scheme.
Traditional thinking using CFD auxiliary process of optimization is judge that parameters are unsatisfactory for expected require
In the case of modify model repeatedly and calculate again, stop calculating when the result of output meets preset value recycle and export it is final
Scheme.In these modifications and calculating repeatedly, CFD calculates the foundation for completely carrying out, and modifying each time each time
Without specific aim, need to spend longer time.
Summary of the invention
Technical problem: the technical problems to be solved by the present invention are: provide it is a kind of using CFD model to indoor thermal environment
Analysis method, can be improved the analysis efficiency of indoor thermal environment, while treating modified load and carrying out quantitative analysis, so that Indoor Thermal
The design scheme of environment can meet design requirement by primary adjustment.
Technical solution: in order to solve the above technical problems, the technical solution adopted by the present invention is that:
It is a kind of using CFD model to the analysis method of indoor thermal environment, the analysis method the following steps are included:
S10 designs initial scheme according to indoor thermal environment, obtains the boundary condition and primary condition of indoor thermal environment；
S20 calculates indoor cardinal temperature field in primary condition and boundary condition input CFD model；
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 conditions are not met, then entering step S40；
S40 calculates each single source item temperature field using CFD model；
S50 determines load to be modified according to each single source item temperature field；
S60 is modified air conditioning solutions according to load to be modified, so that indoor cardinal temperature field meets design
It is required that.
As preference: the step S40 is specifically included:
S401 judges whether by indoor division subregion air-conditioning system, if it is, entering step S402；If it is not, then
Enter step S403；
S402, which sets interior, n subregion air-conditioning system, to either partition domain m, calculates each subregion using CFD model
The temperature t that load is formed in subregion mi-m, i=1,2, n, m is integer, and 2≤m≤n；
Room temperature primary condition is set as design object temperature by S403, using single source or single cold source as source item,
The boundary condition of source item is inputted into CFD model, waste heat factor is set as insulation or calorific value is zero, counts respectively using CFD model
Calculate each source item temperature field.
As preference: in the step S40, hot factor includes heat source and cold source；Heat source includes human body, illuminates, sets
Standby, building enclosure and heat supply end；Cold source includes air-supply and cooling supply end.
As preference: in the S402, calculating each subregion load and specifically wrapped in the temperature that subregion m is formed
It includes:
Room temperature primary condition is set as design object temperature, successively using hot factor all in each subregion as this
The boundary condition of the source item in the region is inputted CFD model by a subregional source item, remaining subregional hot factor is set as being insulated
Or calorific value is zero, calculates separately the temperature that the subregion load is formed in subregion m using CFD model.
As preference: the S50 is specifically included:
S501 judges whether by 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, each subregion load is calculated in the temperature t of subregion mi-m
With the temperature difference t of the subregion design object temperaturei-m, as shown in formula (1):
Δti-m=ti-m-tm-0Formula (1)
The heat for calculating the heat inflow subregion m of each subregion air-conditioning cooling supply or heat supply, as shown in formula (2):
Qi-m=CMmΔti-mFormula (2)
Wherein: Qi-mThe heat of expression subregion i air-conditioning cooling supply or the heat of heat supply inflow subregion m, i=1,
2, n, m≤n, C indicate air specific heat, MmThe air quality for indicating subregion m, works as MmWhen for positive number, heat stream is indicated
Enter, works as MmWhen for negative, indicate that cooling capacity flows into；
Subregion m load Δ Q to be modified is determined according to formula (3)m:
S503 sets indoor design target temperature as t0, by step 20) obtain cardinal temperature field and design object temperature into
Row compares, and finds out the region r for being unsatisfactory for design object temperature in cardinal temperature field, then calculates separately each hot factor i' in region
The 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
The region sum of design object temperature is unsatisfactory in temperature field；
Hot factor i' is calculated in the region r temperature formed and design object temperature t according to formula (4)0The temperature difference:
Δti',r=ti'.r-t0Formula (4)
Region r load to be modified is calculated according to formula (5):
Wherein, C indicates that air specific heat, ρ indicate atmospheric density, VrIndicate the volume of region r.
As preference: in the S60, being modified and specifically include to air conditioning solutions: changing tuyere position, wind
One of mouth quantity, supply air temperature, radiant panel area pattern or combination.
As preference: in the S10, boundary condition includes building enclosure parameter, air output, supply air temperature, air-supply
Speed, return air amount, the setting temperature of return air temperature, 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.
The utility model has the advantages that compared with prior art, the embodiment of the present invention has the advantages that analysis method of the present invention can
It to improve the analysis efficiency of indoor thermal environment, while treating modified load and carrying out quantitative analysis, so that the design of indoor thermal environment
Scheme can meet design requirement by primary adjustment.In the prior art, when judging that indoor cardinal temperature field is unsatisfactory for designing
When target temperature, will rule of thumb method, primary condition and boundary condition are modified, or recycle CFD calculate it is 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
When cardinal temperature field is unsatisfactory for design object temperature, does not need to recycle CFD calculating, but utilize the calculated each list of CFD
Temperature field under one hot factor effect, and determine load to be modified.It is easy to accomplish to air conditioning with load value to be modified
The amendment of scheme, so that it meets design requirement.The embodiment of the present invention, which does not need to recycle CFD, to be calculated, and is not also needed
Parameters revision is carried out using empirical value, calculation amount is greatly reduced, improves analysis efficiency.Meanwhile it treating modified load and being
Quantitative analysis carries out more accurately amendment to air conditioning solutions, improves analysis precision.
Detailed description of the invention
To make the objectives, technical solutions, and advantages of the present invention clearer, the following drawings can make the present invention further
Detailed description, in which:
Fig. 1 is flow diagram of the 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 that the embodiment of the present invention 1 changes the indoor temperature distribution figure after 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 2 cold emission plate operative temperature field of subregion in the embodiment of the present invention 2；
Fig. 7 (b) is 2 heat transfer across wall operative temperature field of subregion in the embodiment of the present invention 2；
Fig. 7 (c) is 2 Source temperature field of subregion in the embodiment of the present invention 2；
Fig. 8 is the indoor temperature distribution figure after 2 Adjusted Option of the embodiment of the present invention.
Specific embodiment
With reference to the accompanying drawing, technical solution of the present invention is described in detail.
It is provided in an embodiment of the present invention it is a kind of using CFD model to the analysis method of indoor thermal environment, comprising the following steps:
S10 designs initial scheme according to indoor thermal environment, obtains the boundary condition and primary condition of indoor thermal environment.
Step S10 is specifically included: determining that indoor thermal environment designs initial scheme according to Architectural Structure Design drawing；According to heat
Environment Design initial scheme design air adjusts form, obtains base regime；The side of indoor thermal environment is set according to base regime
Boundary's condition and primary condition.Wherein, boundary condition includes building enclosure parameter, air output, supply air temperature, air supply velocity, return air
The setting temperature of amount, return air temperature, the semen donors of radiant panel, the heating load of radiant panel or cold emission plate；Primary condition includes enclosing
Protection structure initial temperature, indoor initial temperature, indoor initial humidity.
S20 calculates indoor cardinal temperature field in primary condition and boundary condition input CFD model.
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 conditions are not met, i.e. there are the regions that temperature is higher or lower than design object temperature for the interior space, then
Enter step S40.
S40 calculates each single source item temperature field using CFD model.
S50 determines load to be modified according to each single source item temperature field.
S60 is modified air conditioning solutions according to load to be modified, so that indoor cardinal temperature field meets design
It is required that.Air conditioning solutions are modified and are specifically included: changing tuyere position, inlet number, supply air temperature, radiation plate suqare
One of mode or combination.
In the method for above-described embodiment, by inputting primary condition and boundary condition in CFD model, calculate indoor basic
Temperature field, then judges whether indoor cardinal temperature field meets design object temperature.Each single source item temperature is calculated using CFD model
Field is spent, for ungratified regional analysis.According to each single source item temperature field, the load to be modified for being unsatisfactory for region is determined.
The load to be modified determines value by calculating to obtain.In this way, being easy to carry out air conditioning solutions according to load to be modified
Amendment, so that indoor cardinal temperature field meets design requirement.
In the prior art, when judging that indoor cardinal temperature field is unsatisfactory for design object temperature, will rule of thumb method,
Primary condition and boundary condition are modified, or recycles CFD and calculates 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, it does not need to recycle CFD calculating, but utilizes the temperature field under the calculated each single hot factor effect of CFD, and determine
Load to be modified.With load value to be modified, the amendment easy to accomplish to air conditioning solutions is wanted so that it meets design
It asks.The embodiment of the present invention, which does not need to recycle CFD, to be calculated, and does not also need to carry out parameters revision using empirical value, significantly
Reduce calculation amount, improves analysis efficiency.Meanwhile treating modified load is quantitative analysis, is carried out to air conditioning solutions
More accurately amendment, improves analysis precision.
In above-described embodiment, indoor thermal environment can be divided into two kinds of situations, and one is by indoor division subregion air-conditioning system
The environment of system, another kind are the environment for not dividing subregion air-conditioning system.Below based on two kinds of varying environments, specifically divided
Analysis.
As preference, the step S40 specifically includes S401-S403:
S401 judges whether by indoor division subregion air-conditioning system, if it is, entering step S402；If it is not, then
Enter step S403.
S402, which sets interior, n subregion air-conditioning system, to either partition domain m, calculates each subregion using CFD model
The temperature t that load is formed in the subregioni-m, i=1,2, n, m is integer, and 2≤m≤n.
In S402, calculates each subregion load and is specifically included in the temperature that the subregion is formed:
Room temperature primary condition is set as design object temperature, successively using hot factor all in each subregion as this
The boundary condition of the source item in the region is inputted CFD model by a subregional source item, remaining subregional hot factor is set as being insulated
Or calorific value is zero, calculates separately the temperature that the subregion load is formed in subregion m using CFD model.
Room temperature primary condition is set as design object temperature by S403, using single source or single cold source as source item,
The boundary condition of source item is inputted into CFD model, waste heat factor is set as insulation or calorific value is zero, counts respectively using CFD model
Calculate each source item temperature field.Wherein, hot factor includes heat source and cold source；Heat source includes human body, illumination, equipment, building enclosure and confession
Hot end.Cold source includes air-supply and cooling supply end.
As preference, the S50 specifically includes S501-S503:
S501 judges whether by 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, each subregion load is calculated in the temperature t of subregion mi-m
With the temperature difference t of the subregion design object temperaturei-m, as shown in formula (1):
Δti-m=ti-m-tm-0Formula (1)
The heat for calculating the heat inflow subregion m of each subregion air-conditioning cooling supply or heat supply, as shown in formula (2):
Qi-m=CMmΔti-mFormula (2)
Wherein: Qi-mThe heat of expression subregion i air-conditioning cooling supply or the heat of heat supply inflow subregion m, i=1,
2, n, m≤n, C indicate air specific heat, MmThe air quality for indicating subregion m, works as MmWhen for positive number, heat stream is indicated
Enter, works as MmWhen for negative, indicate that cooling capacity flows into；
Subregion m load Δ Q to be modified is determined according to formula (3)m:
S503 sets indoor design target temperature as t0, by step 20) obtain cardinal temperature field and design object temperature into
Row compares, and finds out the region r for being unsatisfactory for design object temperature in cardinal temperature field, then calculates separately each hot factor i' in region
The 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
The region sum of design object temperature is unsatisfactory in.
Hot factor i' is calculated in the region r temperature formed and design object temperature t according to formula (4)0The temperature difference:
Δti',r=ti'.r-t0Formula (4)
Region r load to be modified is calculated according to formula (5):
Load to be modified has been determined by step S50.After determining load value to be modified, amendment is easy according to the value
Air conditioning solutions, so that indoor cardinal temperature field meets design requirement.The embodiment of the present invention, which passes through, determines load to be modified,
So that air conditioning solutions amendment is more accurate.
In the embodiment of the present invention, in the calculating of first time base regime, all boundary conditions and primary condition be all by
According to the facts border design conditions carry out.This is same as the prior art.S40 only calculates single hot factor operative temperature field.This boundary
The calculating time-consuming of condition is few, and speed is fast.Nondimensionalization single factors operative temperature be distributed resulting temperature form contribution degree can be with
The distribution met is predicted, thus correction air regulation scheme.
The method of the embodiment of the present invention is suitable for various air conditioning systems, and to building type, there is no limit for air-conditioning shape.
Air conditioning terminal form can be air-conditioning system, be also possible to traditional air-conditioning system.
Example is for two example two below.
Embodiment 1
If subregion air-conditioning system is not divided in indoor thermal environment.
If indoor design target temperature is 26 DEG C.Initial scheme is designed according to indoor thermal environment, by primary condition and boundary
In condition entry CFD model, indoor cardinal temperature field is calculated, as a result as shown in Figure 2.Figure it is seen that indoor thermal environment is set
It counts in the cardinal temperature field that initial scheme is formed, partial region temperature is higher than design object temperature, and partial region temperature, which is lower than, to be set
Count target temperature.Figure it is seen that the heat of the generations such as internal heat resource, lighting apparatus, building enclosure has more flowed to room
Interior upper right quarter space, and indoor lower left quarter space obtains more cooling capacity.
Using single source or single cold source as source item, the boundary condition of source item is inputted into CFD model, waste heat factor is set
It is zero for insulation or calorific value, calculates separately each source item temperature field 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, the same sectional view in each source item temperature field is calculated separately using CFD model, i.e., different hot factor lists
The same sectional view for the indoor temperature field that solely 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 are the temperature difference.Its
In, Fig. 3 (a) is the temperature profile that room lighting is formed.Fig. 3 (b) is the temperature profile that internal heat resource is formed.Fig. 3 (c) is
The temperature profile that air outlet is formed.Fig. 3 is the temperature profile that (d) cold emission plate is formed.The cold emission plate of Fig. 3 (d) reflection
Cooling effect makes target area temperature lower than 1~2 DEG C of design object temperature.
Compared with the cardinal temperature field that will acquire is with design object temperature, finds out in cardinal temperature field and be unsatisfactory for design object
Then the region r of temperature calculates separately the temperature t that each hot factor i' is formed in region rI', r, hot factor is calculated according to formula (4)
I' is in the region r temperature formed and design object temperature t0The temperature difference:
Δti',r=ti'.r-t0Formula (4)
Region r load to be modified is calculated according to formula (5):
Since the air outlet temperature that space is formed indoors reaches the level close with target temperature in upper right quarter, without more
Remaining lowering ability, it is contemplated that changing cold emission Board position.When the present embodiment original state, cold emission plate is located below.It is existing
It is moved up by cold emission plate.
After amendment, indoor temperature distribution figure is as shown in Figure 4.As can be seen from Figure 4: room temperature is 26 ± 0.5 DEG C, is met
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 primary condition and boundary
In condition entry CFD model, indoor cardinal temperature field is calculated, as a result as shown in Figure 5.Sitting posture zone of action is 1.3m hereinafter, station
Appearance zone of action is 2m or less.From Fig. 5 cardinal temperature field situation, determines and carry out the modified region of load (level height 0.2m-
1.3m) it is located in subregion 1.
In subregion air-conditioning, as a hot factor, reality is all arranged in all heat sources and cold source in this region in a region
Border boundary condition, whole region are a hot factors.Successively using hot factor all in each subregion as the one of respective partition
The boundary condition of the source item in the region is inputted CFD model by source item, remaining subregional hot factor is set as insulation or calorific value is
Zero, the temperature that the partition load is formed in subregion 1 is calculated separately using CFD model.In the present embodiment, subregion 1 is supplied to the area
The cooling capacity in domain includes: cold emission plate and fresh air；The load that subregion 1 generates the region includes: heat source.Subregion 2 is supplied to the area
The cooling capacity in domain includes: cold emission plate；The load that subregion 2 is supplied to the region includes: building enclosure.
Fig. 6 (a) is 1 cold emission plate operative temperature field of subregion.Fig. 6 (b) is 1 Source temperature field of subregion.Fig. 6 (c) is
1 fresh air operative temperature field of subregion.Fig. 7 (a) is 2 cold emission plate operative temperature field of subregion.Fig. 7 (b) is 2 heat transfer across wall of subregion
Operative temperature field.Fig. 7 (c) is 2 Source temperature field of subregion.
Each partition load is calculated in the temperature t of subregion 1i-mWith the temperature difference t of the zoning design target temperaturei-m.In this reality
It applies in example, the temperature difference of 1 cold emission plate of subregion and 1 design object temperature of subregion is -1.4 DEG C, and 1 fresh air of subregion and subregion 1 design mesh
The temperature difference for marking temperature is -1.1 DEG C, and the temperature difference of 1 internal heat resource of subregion and 1 design object temperature of subregion is+1.6 DEG C.The cold spoke of subregion 2
The temperature difference for penetrating plate and 1 design object temperature of subregion is -0.7 DEG C, the temperature difference of 2 lighting apparatus of subregion and 1 design object temperature of subregion
It is+0.2 DEG C, the temperature difference of 2 building enclosure of subregion and 1 design object temperature of subregion is+0.1 DEG C.Wherein, the temperature difference is positive value, is indicated
Temperature is higher than 1 design object temperature of subregion；The temperature difference is negative value, indicates that temperature is lower than 1 design object temperature of subregion.
Δt1-1=-1.4-1.1+1.6=-0.9 DEG C
Δt2-1=-0.7+0.2+0.1=-0.4 DEG C
The heat for calculating the heat inflow subregion 1 of each subregion air-conditioning cooling supply or heat supply, as shown in formula (2):
Qi-m=CMmΔti-mFormula (2)
In the present embodiment,
The load Δ Q to be modified of subregion 1 is determined according to formula (3)m:
In the present embodiment,
Load corrects area's temperature and is lower than comfort temperature, illustrates that the region has the cooling capacity of additional 0.057kw to flow into the region.If
Review amendment is carried out to the region, then needs modified amount are as follows: reduces 0.057kw cooling capacity and flows into.This example is that signal uses the data
Project setting is carried out, but is not limited only to this method: improving supply air temperature, 0.057kw.This is equivalent to 1.3 DEG C of supply air temperature of raising.
Fig. 8 is after improving supply air temperature, to utilize the indoor cardinal temperature field result of CFD software calculating.It can from Fig. 8
Out: the room temperature after adjusting is 26 ± 0.5 DEG C, meets design requirement.
The basic principles, main features and advantages of the 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 specification be intended merely into
One step illustrates the principle of the present invention, without departing from the spirit and scope of the present invention, the present invention also have various change and
It improves, these changes and improvements all fall within the protetion scope of the claimed invention.The scope of protection of present invention is wanted by right
Ask book and its equivalent thereof.
Claims (7)
1. it is a kind of using CFD model to the analysis method of indoor thermal environment, it is characterised in that: the analysis method includes following step
It is rapid:
S10 designs initial scheme according to indoor thermal environment, obtains the boundary condition and primary condition of indoor thermal environment；
S20 calculates indoor cardinal temperature field in primary condition and boundary condition input CFD model；
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 conditions are not met, then entering step S40；
S40 calculates each single source item temperature field using CFD model；
S50 determines load to be modified according to each single source item temperature field；
S60 is modified air conditioning solutions according to load to be modified, so that indoor cardinal temperature field meets design and wants
It asks.
2. it is described in accordance with the claim 1 using CFD model to the analysis method of indoor thermal environment, it is characterised in that: it is described
Step S40 is specifically included:
S401 judges whether by indoor division subregion air-conditioning system, if it is, entering step S402；If it is not, then into
Step S403；
S402, which sets interior, n subregion air-conditioning system, to either partition domain m, calculates each subregion load using CFD model
In the temperature t that subregion m is formedi-m, i=1,2 ..., n, m is integer, and 2≤m≤n；
Room temperature primary condition is set as design object temperature by S403, using single source or single cold source as source item, by source
The boundary condition of item inputs CFD model, and waste heat factor is set as insulation or calorific value is zero, is calculated separately respectively using CFD model
Source item temperature field.
3. according to claim 2 using CFD model to the analysis method of indoor thermal environment, it is characterised in that: described
In step S403, hot factor includes heat source and cold source；The heat source includes human body, illumination, equipment, building enclosure and heat supply end
End；The cold source includes air-supply and cooling supply end.
4. according to claim 2 using CFD model to the analysis method of indoor thermal environment, it is characterised in that: described
In S402, calculates each subregion load and is specifically included in the temperature that subregion m is formed:
Room temperature primary condition is set as design object temperature, successively using hot factor all in each subregion as the subregion
The boundary condition of the source item in the region is inputted CFD model by one source item in domain, remaining subregional hot factor is set as being insulated or send out
Heat is zero, calculates separately the temperature that the subregion load is formed in subregion m using CFD model.
5. according to claim 2 using CFD model to the analysis method of indoor thermal environment, it is characterised in that: described
S50 is specifically included:
S501 judges whether by indoor division subregion air-conditioning system, if it is, entering step S502；If it is not, then into
Step S503；
S502 sets the design object temperature of subregion m as tm-0, each subregion load is calculated in the temperature t of subregion mi-mWith this
The temperature difference t of subregion design object temperaturei-m, as shown in formula (1):
Δti-m=ti-m-tm-0Formula (1)
The heat for calculating the heat inflow subregion m of each subregion air-conditioning cooling supply or heat supply, as shown in formula (2):
Qi-m=CMmΔti-mFormula (2)
Wherein: Qi-mIndicate that the heat of subregion i air-conditioning cooling supply or heat supply flows into the heat of subregion m, i=1,2, n,
M≤n, C indicate air specific heat, MmThe air quality for indicating subregion m, works as MmWhen for positive number, indicates that heat flows into, work as MmIt is negative
When number, indicate that cooling capacity flows into；
Subregion m load Δ Q to be modified is determined according to formula (3)m:
S503 sets indoor design target temperature as t0, the cardinal temperature field that step 20) obtains is compared with design object temperature
Compared with finding out the region r for being unsatisfactory for design object temperature in cardinal temperature field, then calculate separately each hot factor i' in region r shape
At temperature tI', r, i'=1,2 ..., s, r=1,2, w, wherein s is indoor hot factor sum, and w is substantially warm
The region sum of design object temperature is unsatisfactory in degree field；
Hot factor i' is calculated in the region r temperature formed and design object temperature t according to formula (4)0The temperature difference:
Δti',r=ti'.r-t0Formula (4)
Region r load to be modified is calculated according to formula (5):
Wherein, C indicates that air specific heat, ρ indicate atmospheric density, VrIndicate the volume of region r.
6. it is described in accordance with the claim 1 using CFD model to the analysis method of indoor thermal environment, it is characterised in that: it is described
In S60, air conditioning solutions are modified and are specifically included: changing tuyere position, inlet number, supply air temperature, radiation plate face
One of product mode or combination.
7. it is described in accordance with the claim 1 using CFD model to the analysis method of indoor thermal environment, it is characterised in that: it is described
In S10, boundary condition includes building enclosure parameter, air output, supply air 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；Primary 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 | It is a kind of using CFD model to the analysis method of indoor thermal environment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610867435.5A CN106354977B (en) | 2016-09-29 | 2016-09-29 | It is a kind of using CFD model to the analysis method of indoor thermal environment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106354977A CN106354977A (en) | 2017-01-25 |
| CN106354977B true 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 | It is a kind of using CFD model to the analysis method of indoor thermal environment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106354977B (en) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107025334B (en) * | 2017-03-10 | 2020-05-12 | 国网吉林省电力有限公司 | Method and device for establishing overall dynamic model of heat users of central heating system |
| CN107839581B (en) * | 2017-09-05 | 2020-01-10 | 鲁东大学 | Automatic control feedback carriage and method based on temperature field three-dimensional imaging |
| CN108981097B (en) * | 2018-06-19 | 2020-09-15 | 广东美的暖通设备有限公司 | Air conditioning system and method and device for adjusting energy requirement of air conditioning system |
| CN110826221A (en) * | 2019-11-04 | 2020-02-21 | 北京工业大学 | Method for decoupling indoor temperature field of building |
| JP7130615B2 (en) * | 2019-11-19 | 2022-09-05 | 矢崎エナジーシステム株式会社 | Thermal load calculator |
| CN111426034B (en) * | 2020-05-15 | 2021-06-15 | 北京草木元人工环境技术有限公司 | Indoor air conditioner ventilation air flow organization control system |
| CN111709199B (en) * | 2020-06-28 | 2024-01-26 | 北京工业大学 | Method for measuring heating value of equipment |
| CN114110970B (en) * | 2020-08-31 | 2023-04-14 | 小米科技(武汉)有限公司 | Indoor temperature control method and system for air conditioner |
| CN112464408B (en) * | 2020-11-27 | 2023-03-24 | 四川长虹空调有限公司 | Simulation evaluation method for uniformity of air volume and room temperature field of through-flow air duct air conditioner on-hook |
| CN112923504B (en) * | 2020-12-31 | 2022-08-02 | 广州市昊铭机电科技有限公司 | Classroom environment control device and method |
| CN115164352B (en) * | 2022-05-30 | 2024-05-28 | 中建八局第三建设有限公司 | Whole-process debugging method for large-space air conditioner |
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 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106354977A (en) | 2017-01-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106354977B (en) | It is a kind of using CFD model to the analysis method of indoor thermal environment | |
| Xie et al. | Numerical analysis of temperature non-uniformity and cooling capacity for capillary ceiling radiant cooling panel | |
| Ye et al. | Optimization study of return vent height for an impinging jet ventilation system with exhaust/return-split configuration by TOPSIS method | |
| Cheng et al. | Experimental and numerical investigations on stratified air distribution systems with special configuration: Thermal comfort and energy saving | |
| US20180135879A1 (en) | Air-conditioning control system | |
| CN107781947A (en) | A kind of air conditioning system Cooling and Heat Source forecast Control Algorithm and device | |
| Venko et al. | An experimental study of natural and mixed convection over cooled vertical room wall | |
| CN108006915B (en) | Air conditioner control method and air conditioner | |
| Liang et al. | Cooling load for the design of air terminals in a general non-uniform indoor environment oriented to local requirements | |
| CN105574280B (en) | Glass electric melting work simulation method | |
| CN110826221A (en) | Method for decoupling indoor temperature field of building | |
| Tian et al. | Optimization on thermostat location in an office room using the coupled simulation platform in modelica buildings library: a pilot study | |
| CN109737570A (en) | A kind of multi-connection operation regulation method and system | |
| Zhang et al. | Review and development of the contribution ratio of indoor climate (CRI) | |
| Xue et al. | Influence of floor plenum on energy performance of buildings with UFAD systems | |
| Zhang et al. | Cooling effect of air movement on heating performances of advanced air distribution | |
| CN109241574B (en) | Two-dimensional decomposition modeling method for temperature field of air draft drying section of chain grate | |
| Venko et al. | An experimental study of mixed convection over various thermal activation lengths of vertical TABS | |
| CN106765745B (en) | The method for determining radiation cooling system radiant panel surface temperature variation when supplying water heating | |
| Li et al. | Numerical modeling of non-uniform indoor temperature distribution for coordinated air flow control | |
| JP6213772B2 (en) | Ceiling radiant panel capacity evaluation device and ceiling radiant panel capacity evaluation program | |
| Pan et al. | Optimization on the performances of a novel bed-based task/ambient conditioning (TAC) system | |
| Guan et al. | Numerical investigation of geometry parameters for designing efficient terminal units in active chilled beam | |
| Sinha | Behaviour of inclined jet on room cooling | |
| JPH10281533A (en) | Prediction method of room temperature in floor blowing air conditioning method |
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 |