CN102930145A - Courtyard indoor thermal environment analyzing method based on CFD (Computational Fluid Dynamics) information technology - Google Patents

Abstract

The invention discloses a courtyard indoor thermal environment analyzing method based on a CFD (Computational Fluid Dynamics) information technology. The courtyard indoor thermal environment analyzing method comprises the following steps of: at first, establishing an indoor thermal environment model according to a design drawing of a building structure; determining an airflow organizational way of an indoor thermal environment; setting boundary conditions and initial conditions of the indoor thermal environment, and calculating a temperature field, a wind speed field, an air age, a PMV (Predicted Mean Vote) value and a PPD (Programming Planning Directive) value through CFD simulative calculation; judging whether thermal comfort degree satisfies a pre-set value or not; if not, modifying and adjusting the thermal environment model; and finally, determining and outputting a final scheme of the indoor thermal environment model. The courtyard indoor thermal environment analyzing method provided by the invention simulates the courtyard indoor thermal environment of a large-scale shopping mall building by using a modern information technology and a computational fluid mechanics CFD technology; an optimal manner is selected by respectively analyzing, comparing and evaluating; the design can be improved; the construction debugging of a field air conditioning system can be guided; suggestions for revision can be provided for the points with an improper design; and the requirements on comfort degree, health, security and the like of the shopping mall building by people can be met.

Description

Zhongting indoor thermal environment analytical approach based on the CFD infotech

Technical field

The present invention relates to indoor thermal environment of building analog computation technical field, particularly a kind of Zhongting indoor thermal environment analytical approach based on the CFD infotech.

Background technology

Along with developing rapidly of economy, the urban modernization construction is also thereupon flourish, many complexity, various and large-scale modern architecture have occurred.People improve constantly the live and work environmental requirement, and volumed space building indoor thermal environment and comfortableness thereof are paid close attention to by people gradually.Only with the method for experiment and theoretical analysis the indoor and outdoor thermal environment is carried out Accurate Analysis, be difficult to realize.

For same system and device and equipment, can produce different effects because of different air-flows and Temperature Distribution, simultaneously this system energy-conservation also had vital role.Therefore, if can correctly and in detail predict and estimate indoor air flow and Temperature Distribution in planning and design phase, can realize optimal design and the operational management of modern heating ventilation air-conditioning system.Utilize velocity field that CFD (Computational Fluid Dynamics) technology can form indoor air flows, temperature field, moisture field, concentration of harmful substance field etc. to simulate and predict that this has important directive significance for guaranteeing good heating ventilation air-conditioning system design proposal, improve indoor air quality IAQ (Indoor Air Quality) and reducing building energy consumption.Abroad, the CFD numerical simulation technology has been widely used in engineering field, becomes the upper indispensable technological means of design; And the research of CFD technology and use also progressively comes into one's own at home.

Because the complicacy of shopping center Zhongting large spatial structure, the impact of diversity and indoor and outdoor various factors, the very complex that becomes so that indoor air flow distributes, can not satisfy the indoor environmental condition control requirement of problem by the classic method construction, temperature, humidity and wind field lack of uniformity are to satisfy the requirement of personnel and shopping center environment, so adopt modern information technologies CFD numerical simulation technology means, the air output in each layer of refinement air port, instruct wind system construction to become optimal selection by analog result, the air-conditioning construction that the construction method of problem and achievement also can be similar large space problem offer reference meaning and directive function.

But, if when the center shopping plaza does not also have construction not finish, need to distribute to the air-conditioning system wind field and carry out experimental study, just must set up wind tunnel model and architecture indoor model, this is so that research process is complicated, and cost is expensive, and degree of accuracy is not high yet; Therefore slip-stick artist or architects wish just can predict the distribution situation of room air in the planning and design stage, thereby make best air conditioner scheme.

Summary of the invention

In view of this, technical matters to be solved by this invention provides a kind of Zhongting indoor thermal environment analytical approach based on the CFD infotech.

The object of the present invention is achieved like this:

Zhongting indoor thermal environment analytical approach based on the CFD infotech provided by the invention may further comprise the steps:

S1: set up the indoor thermal environment model according to Architectural Structure Design figure;

S2: determine indoor thermal environment air current composition mode;

S3: boundary condition and the starting condition of setting indoor thermal environment;

S4: air current composition mode, boundary condition and starting condition according to the indoor thermal environment model go out the thermal comfort value by the CFD analog computation;

S5: judge whether the thermal comfort value satisfies preset value, if so, then change step S7 over to;

S6: if not, then return step S1 to the adjustment of making amendment of thermal environment model;

S7: the final plan of determining and export the indoor thermal environment model.

Further, described air current composition mode comprises in following three kinds of modes any: 1) send single deflection grille last time on the double deflection grille; 2) the double deflection grille side is sent, single deflection grille last time; 3) Ball shape nozzle air-supply, single deflection grille last time.

Further, described boundary condition comprises space enclosing structure material parameter, air output, wind pushing temperature, return air amount, return air temperature and exhaust air rate; Described starting condition comprises building enclosure initial temperature, building enclosure initial load and aerial temperature and humidity.

Further, described thermal comfort value comprises temperature field, wind speed field, age of air, pmv value and PPD value.

Further, described preset value be the span in temperature field be 25 ± 0.5 ℃, the velocity amplitude scope of wind speed field for the span of≤1.0m/s, age of air for≤900s, pmv value be-1 ~+1 and the PPD value be≤20%.

Further, among the described step S6 thermal environment model made amendment that to adjust be that mode by adjusting the size of sending the air port back to, quantity, position and air outlet air output in the thermal environment model realizes.

Further, described CFD analog computation is adopted and is improved pressure correcting method, and the realization of described improvement pressure correcting method specifically may further comprise the steps:

S41: the Fluid Pressure Field that is associated with velocity field in the default thermal environment model;

S42: obtain velocity field according to Fluid Pressure Field and present speed field calculating momentum discrete equation;

S43: with the continuity equation in the gained velocity field calculating thermal environment model, obtain the modified value of pressure field;

S44: adjust and upgrade pressure field and velocity field, find the solution Equations of Turbulence and other scalar equations;

S45: judge that whether the calculating on the current time step restrains, if do not restrain, turns back to step S42, iterative computation;

S46: if restrain, return the physical quantity of next time step of step S41 double counting.

The invention has the advantages that: the present invention adopts the Zhongting indoor thermal environment of modern information technologies Fluid Mechanics Computation (CFD) technical modelling shopping mall building.Be that the air current composition mode of each layer top layout air port (double deflection grille or Ball shape nozzle), airflow field (air-supply, air draft air current composition situation), age of air is carried out indoor CFD simulation according to design proposal to all-air system, through after the analog computation, draw indoor temperature field, moisture field, pmv value and PPD value, analyze respectively again, compare and estimate, select optimum way, therefore, research process is simple, cost is low, degree of accuracy is high; Can further improve design, the construction debugging of guide field air-conditioning system to the design local suggestion for revision that proposes improperly, is satisfied the shopping center buildings to the requirement of personnel's comfortableness, health, security etc.

Other advantage of the present invention, target and feature will be set forth to a certain extent in the following description, and to a certain extent, based on being apparent to those skilled in the art to investigating hereinafter, perhaps can obtain from the practice of the present invention instruction.The objectives and other advantages of the present invention can realize and obtain by specifically noted structure in following instructions and the accompanying drawing.

Description of drawings

In order to make the purpose, technical solutions and advantages of the present invention clearer, the present invention is described in further detail below in conjunction with accompanying drawing, wherein:

The Zhongting indoor thermal environment analytical approach process flow diagram based on the CFD infotech that Fig. 1 provides for the embodiment of the invention;

The two-dimensional flow field that Fig. 2 provides for the embodiment of the invention improves pressure correction method process flow diagram.

Embodiment

Below with reference to accompanying drawing, the preferred embodiments of the present invention are described in detail; Should be appreciated that preferred embodiment only for the present invention is described, rather than in order to limit protection scope of the present invention.

Embodiment 1

The Zhongting indoor thermal environment analytical approach process flow diagram based on the CFD infotech that Fig. 1 provides for the embodiment of the invention, as shown in the figure: the Zhongting indoor thermal environment analytical approach based on the CFD infotech provided by the invention is characterized in that: may further comprise the steps:

S1: set up the indoor thermal environment model according to Architectural Structure Design figure;

S2: determine indoor thermal environment air current composition mode; Described air current composition mode comprises in following three kinds of modes any: 1) send single deflection grille last time on the double deflection grille; 2) the double deflection grille side is sent, single deflection grille last time; 3) Ball shape nozzle air-supply, single deflection grille last time.

S3: boundary condition and the starting condition of setting indoor thermal environment; Described boundary condition comprises space enclosing structure material parameter, air output, wind pushing temperature, return air amount, return air temperature and exhaust air rate; Described starting condition comprises building enclosure initial temperature, building enclosure initial load and aerial temperature and humidity.

S4: air current composition mode, boundary condition and starting condition according to the indoor thermal environment model go out the thermal comfort value by the CFD analog computation; Described thermal comfort value comprises temperature field, wind speed field, age of air, pmv value and PPD value.

Described CFD analog computation is adopted and is improved pressure correcting method, and the realization of described improvement pressure correcting method specifically may further comprise the steps:

S41: the Fluid Pressure Field that is associated with velocity field in the default thermal environment model;

S42: obtain velocity field according to Fluid Pressure Field and present speed field calculating momentum discrete equation;

S43: with the continuity equation in the gained velocity field calculating thermal environment model, obtain the modified value of pressure field;

S44: adjust and upgrade pressure field and velocity field, find the solution Equations of Turbulence and other scalar equations;

S45: judge that whether the calculating on the current time step restrains, if do not restrain, turns back to step S42, iterative computation;

S46: if restrain, return the physical quantity of next time step of step S41 double counting.

S5: judge whether the thermal comfort value satisfies preset value, if so, then change step S7 over to; Described preset value be the span in temperature field be 25 ± 0.5 ℃, the velocity amplitude scope of wind speed field for the span of≤1.0m/s, age of air for≤900s, pmv value be-1 ~+1 and the PPD value be≤20%.

S6: if not, then return step S1 to the adjustment of making amendment of thermal environment model; Among the described step S6 thermal environment model being made amendment, to adjust be that mode by adjusting the size of sending the air port back to, quantity, position and air outlet air output in the thermal environment model realizes.

S7: the final plan of determining and export the indoor thermal environment model.

Embodiment 2

Present embodiment 2 only is following content with the difference of embodiment 1: the concrete calculating of following detailed description thermal comfort value, and described thermal comfort value comprises temperature field, wind speed field, age of air, pmv value and PPD value, computation process is as follows:

Wherein, u, v, w represent respectively x, y, and the speed component of z direction is the situation that three-dimensional directly coordinate system is found the solution; U, v be x respectively, and the speed component of y direction is the situation that the direct coordinate system of two dimension is found the solution.

According to CFD numerical simulation analysis method, set up mass-conservation equation, momentum conservation equation and energy conservation equation; Comprise the multiple basic parameters such as speed, temperature, pressure, geometric sense in these equations, want these parameter distribution on the solution room, mass-conservation equation, momentum conservation equation, energy conservation equation, (reality/ideal) equation of gas state and turbulent flow equation (tubulence energy equation and diffusivity equation) are carried out the simultaneous numerical solution.

(1) accounting temperature field and wind speed field

According to parameters such as the continuous differential equation, the momentum conservation differential equation (x, y, three directions of z), the energy conservation differential equation, the Turbulent Kinetic differential equation and turbulent dissipation differential equation simultaneous solution speed, pressure, temperature, Turbulent Kinetic, turbulence dissipation rates, obtain each parameter distributed in three dimensions field.All solving equation common versions:

$\mathrm{div}\left(\mathrm{\ρ}\stackrel{\→}{U}\mathrm{\φ}\right)=\mathrm{div}\left({\mathrm{\Γ}}_{\mathrm{\φ}}\mathrm{grad}\mathrm{\φ}\right)+S_{\mathrm{\φ}}$

In the formula, φ is common variable, the representative amount of finding the solution 1, u, v, w, T, C, k, ε in each governing equation; Γ _φBe coefficient of diffusion corresponding to variable φ; S _φBe source item corresponding to variable φ.One on the left side is convective term, and first on the right is diffusion term, and second is the broad sense source item.

The coefficient of diffusion of universal equation and source item under table 1 three-dimensional rectangular coordinate

Wherein, μ _e=μ+μ _tμ _t=C _μρ k ²/ ε

Γ _e＝μ/p _r+μ _t/σ _t Γ _c＝μ/p _r+μ _t/σ _c

Γ _k＝μ+μ _t/σ _k Γ _ε＝μ+μ _t/σ _ε

${\mathrm{\σ}}_{\mathrm{\ϵ}}=\frac{k^2}{(C_2-C_1)C_{\mathrm{\μ}}^{1/2}}$

$G_k={\mathrm{\μ}}_t\left\{2\right[{\left(\frac{\∂u}{\∂x}\right)}^2+{\left(\frac{\∂v}{\∂y}\right)}^2+{\left(\frac{\∂w}{\∂z}\right)}^2]+{(\frac{\∂u}{\∂y}+\frac{\∂v}{\∂x})}^2+{(\frac{\∂u}{\∂z}+\frac{\∂w}{\∂x})}^2+{(\frac{\∂v}{\∂z}+\frac{\∂w}{\∂y})}^2\}$

$G_B=\mathrm{\βg}\frac{{\mathrm{\μ}}_t}{{\mathrm{\σ}}_t}\frac{\∂T}{\∂y}$

Empirical constant in the table 2k-ε model

C μ	C D	C 1	C 2	σ c	σ k	σ t
							0.09	1.0	1.44	1.92	0.9	1.0	0.9

(2) calculate age of air

Age of air is the time that certain some place air has been detained in the room in the room, has reflected the freshness of room air, and it can comprehensively weigh the ventilation effect in room, is the important indicator of estimating indoor air quality.Also available above universal equation is found the solution.

(3) calculate pmv value and PPD value

PMV-PDD thermal comfort model is the model that Fan Geer (P.O.Fanger) professor of Denmark proposes to estimate the indoor thermal comfort degree.The evaluation index of the sign Studies of Human Body Heat reaction (warm-cold sensitivity) that PMV is, hotness and PMV respective value see Table 3.PPD is the predicted value for the unsatisfied number quantitative of heat.

The scale of table 3PMV hotness

Hotness	Cold	Cool	Little cold	Moderate	Little warm	Warm	Heat
								Pmv value	-3	-2	-1	0	1	2	3

PMV can draw according to following formula:

PMV＝(0.303exp(-0.036M)+0.028){(M-W)-3.05×10 ^-3×[5733-6.99(M-W)-P _a]-0.42[(M-W)-58.15]-1.7×10 ^-5M(5867-P _a)-0.0014M(34-T)-3.96×10 ^-8f _c1[(T _c1+273) ⁴-(T _mrt+273) ⁴]-f _c1h _c(T _c1-T)}

Wherein, M---body metabolism rate, W/m ²

The mechanical work that W---human body is done, W/m ²

p _a---steam partial pressure, unit are P _a

f _C1---clothing surface coefficient, %;

T _C1---clothes hull-skin temperature, unit are a ° C;

T _Mrt---environment mean radiant temperature, unit are a ° C;

h _c---convection transfer rate, unit are W/ (m ²KP _a);

T---human body ambient air temperature, unit are a ° C;

In the formula, M, W, p _aWith T be known quantity, T _MrtIn the MRTEMP submodule, calculate T _C1, h _cAnd f _C1Calculated by following formula:

T _c1＝35.7-0.028(M-W)-I _c1(3.96×10 ^-8f _c1)×[(T _c1+273) ⁴-(T _mrt+273) ⁴]+f _c1h _c(T _c1-T)；

h _cWith air velocity V _rIt is relevant,

When $2.38{(T_{c1}-T)}^{0.25}>12.1\sqrt{V_r},$ h _c＝2.38(T _c1-T) ^0.25

$2.38{(T_{c1}-T)}^{0.25}<12.1\sqrt{V_r},$ $h_c=12.1\sqrt{V_r};$

f _C1With clothing thermal resistance I _C1Relevant

Work as I _C1≤ 0.078m ²° C/W, f _C1=1.00+1.290I _C1

I _c1＞0.078m ²·°C/W，f _c1＝1.05+0.645I _c1

After determining pmv value, calculate the PPD value according to following formula:

PPD＝100-95exp(-(0.03353)PMV ⁴+0.2179PMV ²)

Described preset value be the span in temperature field be 25 ± 0.5 ℃, the velocity amplitude scope of wind speed field for the span of≤1.0m/s, age of air for≤900s, pmv value be-1 ~+1 and the PPD value be≤20%.

Concrete adjustment process is as follows: not at preset range, then should send air port size, quantity and position back to by adjustment if the thermal comfort in indoor occupant main activities zone is judged index, the modes such as air outlet air output are optimized the indoor air flow organizational form.

(1) local relief's amount can be considered to reduce in temperature＜24.5 ℃; Local relief's amount can be considered to increase in temperature＞25.5 ℃;

(2) speed＞1.0m/s can consider to increase the air outlet size or adjust the air outlet form, or adjust setting height(from bottom) etc.;

(3) age of air〉900s, can consider to adjust the air outlet air quantity, send tuyere position etc. back to;

(4) PMV 〉+1, need consider temperature, speed and age of air value and adjust, generalized case can increase air outlet quantity or local relief amount; PMV＜-1 need consider temperature, speed and age of air value and adjust, and generalized case can reduce air outlet quantity or local relief amount

(5) PPD＞20% need consider temperature, speed, age of air and PMV and adjust.

Among the described step S6 thermal environment model being made amendment, to adjust be that mode by adjusting the size of sending the air port back to, quantity, position and air outlet air output in the thermal environment model realizes.

Embodiment 3

Present embodiment 3 only is following content with the difference of embodiment 1:

S41: the Fluid Pressure Field that is associated with velocity field in the default thermal environment model;

S42: obtain velocity field according to coefficient a and constant term b in Fluid Pressure Field and the present speed field calculating momentum discrete equation, the variable quantity of fluid infinitesimal momentum is that gravity, each surface pressing are poor, the acting in conjunction of viscous force, take two dimension as example x and y direction momentum discrete equation can be expressed as follows:

$a_{i,J}{u}_{i,J}^{*}=\mathrm{\&Sigma;}{a}_{\mathrm{nb}}{u}_{\mathrm{nb}}^{*}+(p_{I-1,J}^{*}-p_{I,J}^{*})A_{i,J}+b_{i,j};$

$a_{I,j}{v}_{I,j}^{*}=\mathrm{\&Sigma;}{a}_{\mathrm{nb}}{v}_{\mathrm{nb}}^{*}+(p_{I,J-1}^{*}-p_{I,J}^{*})A_{I,j}+b_{I,j};$

Wherein, a subscript (I, j, nb) in the above-mentioned equation is coefficient a, and b subscript (I, j, nb) is constant term b.

S43: with the poor continuity of each unit in the gained velocity field calculating thermal environment model, described poor continuity deducts the difference of discharge for the influx of each unit;

Calculate u by following momentum discrete equation ^*, v ^*:

$a_{i,J}{u}_{i,J}^{*}=\mathrm{\&Sigma;}{a}_{\mathrm{nb}}{u}_{\mathrm{nb}}^{*}+(p_{I-1,J}^{*}-p_{I,J}^{*})A_{i,J}+b_{i,j};$

$a_{I,j}{v}_{I,j}^{*}=\mathrm{\&Sigma;}{a}_{\mathrm{nb}}{v}_{\mathrm{nb}}^{*}+(p_{I,J-1}^{*}-p_{I,J}^{*})A_{I,j}+b_{I,j};$

Wherein, u ^*The speed component that represents certain x of face place direction, v ^*The speed component that represents certain y of face place direction, i and I are x direction position subscript, and j and J are y direction position subscript, and nb is the node location subscript, when the position subscript of x and y direction all is capitalization, and the position of expression node, otherwise, the position of expression chain of command.

α _{I, J}The recombination coefficient that expression chain of command (x=i, y=J) is located,

Expression chain of command (x=i, y=J) is located the speed component of x direction, a _NbThe recombination coefficient of expression node nb,

Represent node nb at the speed component of x direction,

The force value that expression node (x=I-1, y=J) is located, The force value that node (x=I, y=J) is located, A _{I, J}The area value of expression chain of command (x=i, y=J); b _{I, j}The source item value that expression chain of command (x=i, y=j) is located; a _{I, j}The recombination coefficient that expression chain of command (x=i, y=j) is located,

Expression chain of command (x=I, y=j) is located the speed component of x direction, a _NbThe recombination coefficient at expression node nb place,

The speed component of the expression node nb y of place direction,

The force value that expression node (x=I, y=J-1) is located,

The force value that node (x=I, y=J) is located, A _{I, j}The area value of expression chain of command (x=I, y=j), b _{I, j}The source item value of expression chain of command (x=I, y=i);

S44: according to speed u ^*, v ^*Come as follows the calculating pressure update equation and obtain P ':

a _I,Jp′ _I,J＝a _I+1,Jp′ _I+1，J+a _I-1，Jp′ _I-1，J+a _I,J+1p′ _I,J+1+a _I,J-1p′ _I,J-1+b _I′ _,J；

Wherein, a _{I, J}The recombination coefficient of expression node (x=I, y=J), p ' _I, _JThe pressure correction value that expression node (x=I, y=J) is located; a _{I+1, J}The recombination coefficient of expression node (x=I+1, y=J), p ' _{I+1, J}The pressure correction value that expression node (x=I+1, y=J) is located; a _{I-1, J}The recombination coefficient of expression node (x=I-1, y=J), p ' _{I-1, J}The pressure correction value that expression node (x=I-1, y=J) is located; a _{I, J+1}The recombination coefficient of expression node (x=I, y=J+1), p ' _{I, J+1}The pressure correction value that expression node (x=I, y=J+1) is located; a _{I, J-1}The recombination coefficient of expression node (x=I, y=J-1), p ' _{I, J-1}The pressure correction value that expression node (x=I, y=J-1) is located; B ' _{I, J}The source item that expression node (x=I, y=J) is located.

S45: adjust pressure field and velocity field by following formula, obtain velocity field and newly be worth;

$p_{I,J}=p_{I,J}^{*}+p_{I,J}^{\&prime;},$

$u_{i,J}=u_{i,J}^{*}+d_{i,J}(p_{I-1,J}^{\&prime;}-p_{I,J}^{\&prime;}),$

$v_{I,j}=v_{I,j}^{*}+d_{I,j}(p_{I,J-1}^{\&prime;}-p_{I,J}^{\&prime;}),$

Wherein, p _{I, J}The force value that node (x=I, y=J) is located when representing this iterative computation,

The force value that node (x=I, y=J) was located when the expression last iteration calculated, p ' _{I, J}The pressure correction value that expression node (x=I, y=J) is located, u _{I, J}The expression chain of command,

Chain of command (x=i, y=J) was located the speed component of x direction, d when the expression last iteration calculated _{I, J}(p ' _{I-1, J}-p ' _{I, J}) chain of command (x=i, y=J) is located the modified value of the speed component of x direction, v when representing this iterative computation _{I, j}Chain of command (x=I, y=j) is located the speed component of y direction when representing this iterative computation,

Chain of command (x=I, y=j) was located the speed component of y direction, d when the expression last iteration calculated _{I, j}(p ' _{I, J-1}-p ' _{I, J}) chain of command (x=I, y=J) is located the modified value of the speed component of y direction when representing this iterative computation;

S46: by following Formula For Solving φ:

a _I,Jφ _I,J＝a _I+1,Jφ _I+1,J+a _I-1,Jφ _I-1,J+a _I,J+1φ _I,J+1+a _I,J-1φ _I,J-1+b _I,J，

Wherein, φ can distinguish value u, v, T, C, k, ε (speed component, temperature, concentration, tubulence energy, the parameters such as turbulence diffusivity), α _{I, J}φ _{I, J}, α _{I+1, J}φ _{I+1, J}, α _{I-1, J}φ _{I-1, J}, α _{I, J+1}φ _{I, J+1}, α _{I, J-1}φ _{I, J-1}, b _{I, J}

S46: whether the solution of judging the discretize transport equation restrains, and if not, then carries out as follows assignment, and turns back to step S42:

p ^*＝p,u ^*＝u，v ^*＝v,φ ^*＝φ，

Wherein, p ^*The force value that the expression last iteration calculates, p represents this iterative computation force value; u ^*The speed component value of the x direction that the expression last iteration calculates, u represents the speed component value of this iterative computation x direction; v ^*The speed component value of the y direction that the expression last iteration calculates, v represents the speed component value of this iterative computation y direction; φ ^*The variate-value (φ can represent u, v, w, T, C, k, ε) that the expression last iteration calculates, φ represents the variate-value of this iterative computation.

S47: if then export the solution of discretize transport equation.

The above is the preferred embodiments of the present invention only, is not limited to the present invention, and obviously, those skilled in the art can carry out various changes and modification and not break away from the spirit and scope of the present invention the present invention.Like this, if of the present invention these are revised and modification belongs within the scope of claim of the present invention and equivalent technologies thereof, then the present invention also is intended to comprise these changes and modification interior.

Claims (7)

1. based on the Zhongting indoor thermal environment analytical approach of CFD infotech, it is characterized in that: may further comprise the steps:

S1: set up the indoor thermal environment model according to Architectural Structure Design figure;

S2: determine indoor thermal environment air current composition mode;

S3: boundary condition and the starting condition of setting indoor thermal environment;

S4: air current composition mode, boundary condition and starting condition according to the indoor thermal environment model go out the thermal comfort value by the CFD analog computation;

S5: judge whether the thermal comfort value satisfies preset value, if so, then change step S7 over to;

S6: if not, then return step S1 to the adjustment of making amendment of thermal environment model;

S7: the final plan of determining and export the indoor thermal environment model.

2. the Zhongting indoor thermal environment analytical approach based on the CFD infotech according to claim 1 is characterized in that: described air current composition mode comprises in following three kinds of modes any: 1) send single deflection grille last time on the double deflection grille; 2) the double deflection grille side is sent, single deflection grille last time; 3) Ball shape nozzle air-supply, single deflection grille last time.

3. the Zhongting indoor thermal environment analytical approach based on the CFD infotech according to claim 1, it is characterized in that: described boundary condition comprises space enclosing structure material parameter, air output, wind pushing temperature, return air amount, return air temperature and exhaust air rate; Described starting condition comprises building enclosure initial temperature, building enclosure initial load and aerial temperature and humidity.

4. the Zhongting indoor thermal environment analytical approach based on the CFD infotech according to claim 1, it is characterized in that: described thermal comfort value comprises temperature field, wind speed field, age of air, pmv value and PPD value.

5. the Zhongting indoor thermal environment analytical approach based on the CFD infotech according to claim 1 is characterized in that: described preset value be the span in temperature field be 25 ± 0.5 ℃, the velocity amplitude scope of wind speed field for the span of≤1.0m/s, age of air for≤900s, pmv value be-1 ~+1 and the PPD value be≤20%.

6. the Zhongting indoor thermal environment analytical approach based on the CFD infotech according to claim 1 is characterized in that: among the described step S6 thermal environment model being made amendment, to adjust be that mode by adjusting the size of sending the air port back to, quantity, position and air outlet air output in the thermal environment model realizes.

7. the Zhongting indoor thermal environment analytical approach based on the CFD infotech according to claim 1 is characterized in that: described CFD analog computation is adopted and is improved pressure correcting method, and the realization of described improvement pressure correcting method specifically may further comprise the steps:

S41: the Fluid Pressure Field that is associated with velocity field in the default thermal environment model;

S42: obtain velocity field according to Fluid Pressure Field and present speed field calculating momentum discrete equation;

S43: with the continuity equation in the gained velocity field calculating thermal environment model, obtain the modified value of pressure field;

S44: adjust and upgrade pressure field and velocity field, find the solution Equations of Turbulence and other scalar equations;

S45: judge that whether the calculating on the current time step restrains, if do not restrain, turns back to step S42, iterative computation;

S46: if restrain, return the physical quantity of next time step of step S41 double counting.

Images