CN104871159A - System and method for fluid dynamics prediction with an enhanced potential flow model - Google Patents

Abstract

A system and method for modeling airflow and temperature are disclosed. In one example, the method includes receiving input data related to a physical layout of a facility, dividing, by a computer, a representation of the facility into a plurality of grid cells, identifying where effects of at least one of jet airflow, thermal plumes and buoyancy forces are present in the facility based on the physical layout, specifying a velocity value, using a velocity correction method, for a first set of the plurality of grid cells if the effects of at least one of jet airflow and thermal plumes are present within the first set of the plurality of grid cells, calculating, by the computer, an airflow velocity value for each of a second set of the plurality of grid cells, the second set being different from the first set, modifying the determined airflow velocity value for any of the second set of the plurality of grid cells where the effects of buoyancy forces are present, and storing, on a storage device, the modified airflow values.

Description

Utilize the system and method that the fluid dynamics of enhancement mode potential-flow model is predicted

Background

Invention field

At least one embodiment according to the present invention relates generally to for the management of the refrigeration system of the heating in buildings, ventilating and air conditioning system and data center and designed system and method.

the discussion of correlation technique

In response to the demand of the increase of the economy based on information, information technology network continues diffusion in the world.The one performance of this growth is centralized network data center.Centralized network data center is usually by providing the layout of network connectivty, electric power and cooling power various Information Technology Equipments in the structure to form.Described equipment is often accommodated in the special shell being called " frame ", this special shell integrated these connectedness, electric power and cooling units.In some data center configuration, these row are organized into heat with cold passage, to reduce and the relevant cost of cooling Information Technology Equipment.Mask has the movable floor of air chamber to be generally used for as frame provides cooling-air under floor.Cold air is distributed to frame by the soaps with open-air area from air chamber.

Various process and software application, such as can change the commercially available data center management system of (APC) company by the American Electric Power under the Schneider Electric Devices of western Kingston, Rhode Island, be developed thus helped data center operations personnel to design and keep efficient and effective data center configuration.These instruments carry out guide data care workers often through following activity, and these activities are design data division center, before the mounting intracardiac positioning equipment and reorientate equipment after structure and installment completes in the data such as.Therefore, traditional tool set provides standardization and predictable method for designing to data center operations personnel.

In addition, process and software application have been developed the common buildings design for heating, ventilating and air conditioning system (HVAC system).HVAC system can use pipelines to scatter the air through regulating at whole buildings.According to indoor temperature, HVAC system can regulate through regulate air to the flowing in room.Instrument for designing HVAC system can instruct building designers by following activity, and these activities such as provide comfortable environment at identification space inner position pipeline, grid and air diffuser so that for the user in this space.

Summary of the invention

According to an aspect, the computer implemented method for modeling air-flow comprises: receive the input data relevant to the physical layout of facility; By computing machine, the expression of this facility is divided into multiple grid cell; Based on this physical layout, identify that the effect of at least one in jet air-flow, Thermal plume and buoyancy is present in the where in this facility; If the effect of at least one in jet air-flow and Thermal plume is present in first group of multiple grid cell, then the velocity amplitude of first group being used for the plurality of grid cell is specified in operating speed bearing calibration; Be used for each gas velocity angle value in second group of multiple grid cell by computer calculate, second group is different from first group; Amendment is used for the gas velocity of the determination of any one grid cell in second group of the multiple grid cells that there is buoyancy wherein; Estimate the air flow value of memory modify on the storage means.

In at least one embodiment, first group of multiple grid cell is sky group.In another embodiment, second group of multiple grid cell is sky group.In the method, each of multiple grid cell can comprise multiple.In some examples, command speed value can comprise the velocity amplitude of specifying at least one face in multiple.In other embodiments, the method also comprises the gas velocity angle value calculated at least another face in multiple.In at least one example, command speed value can comprise the velocity amplitude of specifying for the face of aliging with Thermal plume.In one example, Thermal plume is associated with vertical direction.In another example, command speed value can comprise the velocity amplitude of specifying for the face of aliging with jet air-flow.In one example, jet air-flow is associated with horizontal direction.In other embodiments, the method also comprises the gas velocity angle value calculated at least another face in multiple.

In the method, calculated gas flow velocity amplitude also comprises and uses potential barrier method to carry out calculated gas flow velocity amplitude.In one example, the method also comprises the equipment be configured in based at least one in the air flow value of amendment and the velocity amplitude of specifying in facility.In another example, the method also comprises assesses the hot comfort of user in facility based at least one in the air flow value revised and the velocity amplitude of specifying.In Still another example, the air flow value that the method also comprises based on amendment controls the equipment in facility with at least one in the velocity amplitude of specifying.

In another example, the method also comprises: calculated for each temperature value in multiple grid cell by computing machine based on gas velocity angle value.In the method, this facility can be included in the space in data center, and the object in this physical layout comprises at least one equipment rack and at least one refrigeration supply arrangement.In the method, facility is included in the space in buildings, and the object in physical layout comprises at least one aeration structure and at least one heating equipment.In the method, revising the gas velocity angle value determined can also comprise: the buoyancy velocity amplitude calculating any one grid cell in second group that is used for multiple grid cell, and buoyancy velocity amplitude is added to the gas velocity angle value determined.

In the method, command speed value also comprises: the jet gas velocity angle value calculating any one grid cell of first group of the multiple grid cells in the zone of influence be used in facility, and the jet gas velocity angle value of first group of specifying the multiple grid cells be used in this zone of influence.In the method, command speed value can also comprise: the haloing unit speed value calculating any one grid cell in first group of multiple grid cells of the object be used in contiguous facility, and specifies the haloing unit speed value of first group being used for multiple grid cells that this object contiguous is arranged.

In one example, the method also comprises: use process of iteration to determine each the new air flow value in the grid cell in this facility, wherein, this new air flow value meets mass balance equation.In at least one example, this process of iteration uses the method for physically based deformation, and the method is configured to non-mass conservation airflow field to be converted to mass conservation field, preserves the lead characteristic in initial non-conservation flow field simultaneously.In some examples, the method also comprises: determine whether the difference between new air flow value and air flow value is before greater than threshold value, and repeat this process of iteration until this difference is not more than this threshold value.

According to another aspect, a kind of system for modeling air-flow, this system comprises: storer and the processor being coupled to this storer, and this system is configured to receive the input data relevant to the physical layout of facility; The expression of this facility is divided into multiple grid cell; Effect based at least one in this physical layout identification jet air-flow, Thermal plume and buoyancy is present in the where in this facility; If the effect of at least one in jet air-flow and Thermal plume is present in first group of multiple grid cell, then the velocity amplitude of first group for the plurality of grid cell is specified in operating speed bearing calibration; Be used for each gas velocity angle value in second group of multiple grid cell by computer calculate, second group is different from first group; Amendment is used for the gas velocity of the determination of any one grid cell in second group of the multiple grid cells that there is buoyancy wherein; And the air flow value of memory modify on the storage means.

In at least one embodiment, first group of multiple grid cell is sky group.In another embodiment, second group of multiple grid cell is sky group.Within the system, each in multiple grid cell can comprise multiple.In some examples, command speed value can comprise the velocity amplitude of specifying at least one in multiple.In at least one example, command speed value can comprise the velocity amplitude of specifying for the face of aliging with Thermal plume.In one example, Thermal plume is associated with vertical direction.In other embodiments, the method also comprises the gas velocity angle value calculated at least another face in multiple.In another example, command speed value can comprise the velocity amplitude of specifying for the face of aliging with jet air-flow.In one example, jet air-flow is associated with horizontal direction.In other embodiments, the method also comprises the gas velocity angle value calculated at least another face in multiple.

In one example, this system is also configured to use potential barrier method to carry out calculated gas flow velocity amplitude.In another example, this system is also configured to the equipment be configured in based at least one in the air flow value revised and the velocity amplitude of specifying in facility.In another example, this system is also configured to based on the air flow value of amendment and the hot comfort of at least one assessment user in facility in the velocity amplitude of specifying.In Still another example, this system is also configured to control the equipment in facility based on the air flow value of amendment with at least one in the velocity amplitude of specifying.

In one example, this system is also configured to calculate for each temperature value in multiple grid cell based on gas velocity angle value.In at least one example, facility is included in the space in data center, and the object in physical layout comprises at least one equipment rack and at least one refrigeration supply arrangement.In another example, facility is included in the space in buildings, and the object in physical layout comprises at least one aeration structure and at least one heating equipment.

In certain embodiments, this system is also configured to revise for the buoyancy velocity amplitude of any one grid cell in second group of multiple grid cell the gas velocity angle value determined by calculating, and this buoyancy velocity amplitude is added to the gas velocity angle value determined.In other embodiments, this system is also configured to carry out command speed value by calculating for the jet gas velocity angle value of any one grid cell in first group of the multiple grid cells in the zone of influence, and specifies the jet gas velocity angle value of first group of the multiple grid cells be used in this zone of influence.

In one embodiment, this system is also configured to carry out command speed value by the haloing unit speed value calculated for any one grid cell in first group of multiple grid cells of the object in contiguous facility, and specifies the haloing unit speed value of first group being used for multiple grid cells that this object contiguous is arranged.This system is also configured to calculate the new air flow value for each in the grid cell in this facility, and wherein, this new air flow value meets mass balance equation.In one example, this system is also configured to use process of iteration to determine the new air flow value of each in the grid cell in this facility, and wherein, this new air flow value meets mass balance equation.In at least one example, this process of iteration uses the method for physically based deformation, and the method is configured to non-mass conservation airflow field to be converted to mass conservation field, preserves the lead characteristic in initial non-conservation flow field simultaneously.

According to another aspect, disclose a kind of non-transitory computer-readable medium, it has the instruction sequence stored for modeling air-flow thereon.In one example, this non-transitory computer-readable medium comprises instruction, and this instruction receives the input data relevant to the physical layout of facility by causing at least one processor; The expression of this facility is divided into multiple grid cell; Effect based at least one in this physical layout identification jet air-flow, Thermal plume and buoyancy is present in the where in this facility; If the effect of at least one in jet air-flow and Thermal plume is present in first group of multiple grid cell, then the velocity amplitude of first group for the plurality of grid cell is specified in operating speed bearing calibration; Be used for each gas velocity angle value in second group of multiple grid cell by computer calculate, second group is different from first group; Amendment is used for the gas velocity of the determination of any one grid cell in second group of the multiple grid cells that there is buoyancy wherein; And the air flow value of memory modify on the storage means.

Accompanying drawing explanation

Accompanying drawing is not intended to draw in proportion.In the accompanying drawings, each identical or approximately uniform assembly illustrated in the various figures is represented by identical label.For clarity sake, be not that each assembly is marked in each accompanying drawing.In the accompanying drawings:

Fig. 1 is the schematic diagram of an example of the distributed system comprising design and management system;

Fig. 2 is the schematic diagram of an example of use grid cell according to an embodiment;

Fig. 3 is the process flow diagram of the instantiation procedure for determining air-flow and temperature according to an embodiment;

Fig. 4 is a schematic diagram represented of different air-flow;

Fig. 5 is the schematic diagram comprising a layout in the room of modeling jet air-flow according to an embodiment;

Fig. 6 A is the schematic diagram of the heating block comprising the haloing unit (halo cell) that contiguous heating block is arranged;

Fig. 6 B is the schematic diagram of a heating block of the modeling haloing unit comprising the uprush had specified by an embodiment;

Fig. 7 A is the curve map that the uprush specification using haloing unit is shown according to an embodiment;

Fig. 7 B is the schematic diagram that the 3D of heating block represents;

Fig. 7 C is the schematic diagram that the 2D of the heating block comprising face A, B, C represents;

Fig. 8 is the process flow diagram of the instantiation procedure for mapping step method according to an embodiment;

Fig. 9 A is the schematic diagram in the region with boundary condition;

Fig. 9 B is the schematic diagram in the region comprising initial non-mass conservation airflow field;

Fig. 9 C is the schematic diagram in the region comprising mass conservation airflow field;

Figure 10 is the schematic diagram of a layout in the 3D room of the air-flow of modeling wherein according to an embodiment;

Figure 11 A is the schematic diagram of the temperature value using CFD method to determine;

Figure 11 B is the schematic diagram of the temperature value using PFM method to determine;

Figure 11 C is the schematic diagram of the temperature value determined according to the use E-PFM method of an embodiment;

Figure 12 A is the schematic diagram of the air flow value using CFD method to determine;

Figure 12 B is the schematic diagram of the air flow value using PFM method to determine;

Figure 12 C is the schematic diagram of the air flow value determined according to the use E-PFM method of an embodiment;

Figure 13 is the schematic diagram of six layouts of the data center of the air temperature and current of modeling wherein according to an embodiment;

Figure 14 A is the curve map of the mean accuracy comparing the mean air entry temperature using E-PFM and PFM method to calculate;

Figure 14 B is the curve map of the mean accuracy comparing the full admission temperature using E-PFM and PFM method to calculate;

Figure 15 be compare use E-PFM and PFM method to calculate six layouts in the curve map of mean accuracy of each layout; And

Figure 16 is the block diagram of an example of the computer system can implemented whereby according to various aspects of the present invention.

Embodiment

At least some embodiment according to the present invention relates to user can the system of design and analysis data center configuration and buildings HVAC system and process whereby.These systems and process can help design and analysis activity by the model allowing user to create buildings HVAC system, wherein can accounting temperature and air-flow by this model.System and user all can adopt calculated air temperature and current to be designed for for user provides the alternative arrangements of the buildings HVAC system of comfort level and meet the data center configuration of various design object.According to an embodiment, system and method described herein also can be applied to other application, and comprise Electronic cooling application and indoor environment application, the application of this indoor environment comprises toilet, commercial Application, laboratory, hospital and operating room.

Should be appreciated that and to be calculated by U.S. heating, refrigeration and air-conditioning man SCTE (ASHRAE) and the available simple hand of its hetero-organization and HVAC design tool is often used for as individual space design guarantees the buildings HVAC system of occupant's health and comfort level.Such as, individual space can be included in office in office building, booth or the communal space.But, this type of instrument is usually based on the estimation of detail and the empirical model that may not comprise each individual space.Therefore, current available HVAC design tool has a mind to guard.Therefore, surdimensionnement and unnecessary costliness that the HVAC system of this type of tool design is unnecessary is often used.

In order to the performance of accurately predicting individual space, computational fluid dynamics (CFD) is analyzed sometimes for emulating the flow pattern, temperature and other factors that affect occupant's comfort level, building energy consumption and other buildings characteristics.This information may be used for estimating the validity of each HVAC and natural ventilation system design and assesses thing followed indoor comfortable condition.But, CFD is very expensive, needs professional operation person, and can provide inconsistent prediction according to operator with for the bottom physical model of any given CFD routine package.

As for data-center applications, CFD remains the most widely used instrument of prediction air-flow and temperature.But, due to the defect of CFD, in some instrument, potential barrier modeling (PFM) is analyzed and has been replaced CFD to use recently.Be similar to CFD, PFM provides the complete 3D prediction of the full breadth of air-flow, temperature and its dependent variable.PFM method on Dec 15th, 2011 submit to be entitled as " SYSTEM AND METHOD FOR RACK COOLING ANALYSIS ", the patent cooperation treaty application (being referred to herein as " PCT/US2011/065071 application ") that is numbered PCT/US2011/065071; And on March 23rd, 2012 submit to be entitled as " SYSTEMS AND METHODS FOR PREDICTING FLUID DYNAMICS INA DATA CENTER ", description in the patent cooperation treaty application (being referred to herein as " PCT/US2012/030373 application ") that is numbered PCT/US2012/030373, the assigned assignee giving the application of above-mentioned each application, and the full content of above-mentioned each application is incorporated to herein by reference.

The PFM method real-time estimate described in PCT/US2011/065071 application is by the air intake opening of equipment rack in data center and the air-flow of gas outlet.In real time or quasi real time the ability of estimated performance allows several possible solution of instrument express-analysis to be used and considers various design tradeoff.Different from CFD, the simpler physical model that PFM method flows based on fluid, this model can provide several advantage, comprises for relatively simple exploitation and the maintenance of some actual 3D application and relative fast close to real-time.PFM method can be characterized as being stable, because it continues to produce rational result.In addition, PFM method is easy to use, because it does not need professional operation person.

But when realizing these advantages, PFM method may ignore some physical influence, such as jet-like flow pattern.Because PFM method is generally used in data-center applications, the usual Bu Shi data center of ability of modeling jet-like air-flow paid close attention to.Mainly because jet-like air-flow only occurs in specific occasion usually, such as such as, from the jet of the soaps of single frame or isolation.For the actual set of wherein many bricks and the arrangement adjacent one another are of many frames, source array produces evenly and the flow pattern of less jet-like.Such as, when hot-air rises to data center top, PFM method can also ignore other physical influences comprising buoyancy.Illustrate that buoyancy is important for Accurate Model air-flow, specifically because server design has the trend tending to more high power and less air-flow in heart design in the data.In Building Design application, the effect of jet-like flow pattern, thermal blooming and the effect of buoyancy are particular importances.Because PFM method may not be so accurate for the obvious buildings in general example that flows of buoyancy and jet-like wherein, therefore PFM method needs to strengthen for buildings in general design ap-plication.

Therefore, for meeting the demand of buildings in general design community, there is the demand of quick, the simple and stable instrument for the additional heating illustrated around buoyancy, jet-like air-flow and object.The aforementioned restriction that enhancement mode PFM disclosed herein (E-PFM) model alleviates PFM method retains its advantage simultaneously.According to each embodiment, by jet-like flowing and buoyancy are described, E-PFM can meet the demand of buildings in general design ap-plication and provide relatively simple and stable solution.In each example, the time of the generation solution of E-PFM, only close to the twice of PFM, equals the fraction of CFD solution time.

That set forth in many aspects disclosed herein are not limited in its application to embodiment below according to the present embodiment or structure detail illustrated in the accompanying drawings and assembly layout.These aspects can be supposed other embodiments and can put into practice in every way or implement.In the example of particular implementation provided herein only for illustration of object, and be not intended to limit.Particularly, the behavior discussed in conjunction with any one or more embodiments, element and feature are not intended to the similar role got rid of in any other embodiment.

Such as, according to one embodiment of present invention, computer system is configured to perform any one function as herein described, the information including but not limited to configuration, modeling and present about particular data center configuration.In addition, computer system in an embodiment may be used for environmental parameter in automatically measurement data center or buildings and opertaing device (such as condenser or refrigeratory or well heater) with Optimal performance.And system as herein described can be configured to comprise or get rid of any one function discussed in this article.Therefore, embodiment is not limited to specific function or specific one group of function.And phrase used herein and term for describing object, and should not be counted as restriction.Use in this article for " comprising (including) ", " comprising (comprising) ", " having (having) " and " comprising (containing) ", " relating to (involving) " and variant thereof means to be included in its every and equivalent listed and other item below.

exemplary system architecture

Fig. 1 illustrates and comprises the physics of distributed system 100 and the graph of a relation of logic element.As shown in the figure, distributed system 100 carrys out concrete configuration according to current the disclosed embodiments.The system architecture described in detail about Fig. 1 and content, only for exemplary purpose, are not intended to embodiment to limit ad hoc structure shown in Figure 1.When not departing from the scope of current the disclosed embodiments, can construct many variant system architectures, this is apparent for the wherein those of ordinary skill in this area.Specific arrangements shown in Fig. 1 is to promote sharpness by selection.

Information can use any technology and element shown in Figure 1, flow between assembly and subsystem.This type of technology comprises such as via TCP/IP transmission of information on network, transmission of information between the module of storer, and carrys out transmission of information by writing in files, database or some other Nonvolatile memory devices.When not departing from the scope of current the disclosed embodiments, other technologies and agreement can be used.

With reference to figure 1, system 100 comprises user 102, interface 104, design and management system 106, communication network 108 and database 110.System 100 can allow user 102 (such as, data center architecture teacher or buildings in general designer) and interface 104 alternately to create or to revise the model that one or more data center or buildings in general configure.

In modern data center, equipment rack and row's formula refrigeratory are usually with the layout that the front/rear layout replaced is in a row, thus create the passage of heat alternately and cold passage in the data in the heart, wherein often arrange frame above towards cold passage and often arrange frame after towards the passage of heat, as submit on January 24th, 2008 be entitled as " System and Method for Evaluating Equipment Rack Cooling " be numbered 12/019, the U.S. Patent application of 109, it is now United States Patent (USP) 7,991, described in 592; And as submit on January 27th, 2006 be entitled as " Methods and Systems for ManagingFacility Power and Cooling " be numbered 11/342, the U.S. Patent application of 300, be now United States Patent (USP) 7,881, described in 910; Each assigned in above-mentioned application gives the assignee of the application, and the full content of each application in above-mentioned application is incorporated to herein by reference.Cold air to suck before frame and discharges air below from frame by typical equipment rack.In description herein, the equipment in frame or frame itself can be called as cooling consumer, and arrange formula cooling unit and/or air conditioner in machine room (CRAC) can be called as refrigeration supply arrangement.In quoted application, instrument is provided for the cooling performance analyzed the frame in data center and troop.In these instruments, multiple analysis can perform the cooling performance attempting optimization data center in different layouts.

According to an embodiment, interface 104 can comprise the aspect of data center floor editing machine or frame editing machine, as submit on May 15th, 2008 to be entitled as being numbered in the patent cooperation treaty application of PCT/US08/63675 of " METHODS AND SYSTEMSFOR MANAGING FACILITY POWER AND COOLING " disclosed, the full content of this application is incorporated to herein (being referred to herein as " PCT/US08/63675 application ") by reference.In other embodiments, interface 104 can utilize specialized equipment to implement, and it makes user 102 can comprise the model represented of the physical layout of data center or its any subset with drag and drop patten's design.This layout can comprise the expression of data center's construction package and data center apparatus.The feature at interface 104 is discussed below further.In at least one embodiment, the information about data center or building foundation facility is received by interface by system 100, and the assessment of data center and suggestion are provided to user.In addition, at least one embodiment, optimizing process can be performed cooling performance and the energy use of optimizing buildings and/or data center.

In buildings in general application, the HVAC system for individual room or space can comprise air and be recycled to the return air inlet of HVAC system by its air intake opening be supplied and air by it.Individual room can be modeled to be included in the air intake opening in any amount of configuration and return air inlet, such as, bottom room and on top.Individual room can also be modeled to comprise location object in a room usually, such as desk, computing machine, chair or human body.These objects can be modeled to produce heat and stop and be flowed by the air in room.

As shown in Figure 1, design and management system 106 present design interface 104 to user 102.According to an embodiment, design and management system 106 can comprise as data center's design and management system disclosed in PCT/US08/63675.In this embodiment, design interface 104 can be included in the function of load module, display module and the constructor's module that PCT/US08/63675 comprises, and can usage data library module with store and retrieve data.

As shown in the figure, design and management system 106 can exchange information via network 108 and database 110.This information can comprise any information needed for the Characteristic and function of supported data Center and management system 106.Such as, in one embodiment, data center's database 110 can comprise at least some part being stored in the data in data center apparatus database described by PCT/US08/63675.In another embodiment, this information can comprise any information supported needed for interface 104, such as, among other data, the configuration of one or more data centers model physical layout, be included in the consumption feature of cooling consumer in the product of refrigeration supply arrangement in model configuration and distribution character, model configuration and the list of the equipment rack comprised in the cluster and refrigeration supply arrangement.

In one embodiment, data center's database 110 can store the information belonging to refrigeration supply arrangement type, the cooling air volume provided by the refrigeration supply arrangement of each type, and the temperature of the cooling-air provided by refrigeration supply arrangement.Therefore, such as, data center's database 110 comprises the record of the CRAC unit belonging to particular type, and it is rated for the speed place delivery air at 5600 cubic feet (cfm) per minute at 68 degrees Fahrenheit places.In addition, data center's database 110 can store the information belonging to one or more cooling tolerance, the air intake opening of the air intake opening of such as CRAC and air outlet temperature and one or more equipment rack and exhaust port temperatures.Temperature can be input in system by periodic measurement, or in other embodiments, temperature can use the device of the system of being coupled to 100 to monitor continuously.

Data center's database 110 can adopt the form of any logical construct that can store information on a computer-readable medium, except other structures, comprises flat file, index file, hierarchical data base, relational database or OODB Object Oriented Data Base.Data can use the relation of unique key and external key and index to carry out modeling.The relation of unique key and external key and index can be set up, to guarantee the performance of data integrity and exchanges data between each field and table.

Often kind in computer system shown in Figure 1 can comprise one or more computer system, and wherein, this computer system comprises design and management system 106, network 108 and database 110.As above about Fig. 1 discuss, computer system can have one or more processor or controller, storer and interface device.The customized configuration of the system 100 described in FIG is only for illustration of object, and embodiments of the invention can be put into practice in other contexts.Embodiment as herein described is not limited to user or the system of specific quantity.

potential barrier modeling (PFM)

The PFM method of prediction air-flow, temperature and other data center's performance metrics is discussed with reference to PCT/US2012/030373 application.Potential barrier method is the technology of the physically based deformation used in one embodiment, to determine velocity potential in the unit that structured grid is arranged and temperature.Fig. 2 illustrates that 2D is uniform, structured grid arranges 200, has the contiguous grid cell 202,204,206 and 208 having size Δ x.Air-flow can enter or leave every side of grid cell; Additional air-flow can add grid cell to or deduct from grid cell.But, other trellis schemes can be implemented for counting yield, and actual emulation can perform in 3D.

Generally speaking, air-flow is idealized as incompressible and irrotationality by PFM method, and wherein, flow field can be determined from Poisson equation:

Wherein, S " ' be the volume flow rate source of per unit volume, and φ is velocity potential.The relation of the x, y, z component of φ velocity potential and speed is as follows:

When as shown in Figure 2 by discretize for numerical evaluation on computing grid time, can be written as at the velocity potential of any grid cell i:

Can obtain at each grid cell for φ _ithis class equation, and whole group can be solved simultaneously.As discussed in PCT/US2012/030373 application, Fig. 2 illustrates and calculates and staggered-mesh that speed calculates at elemental area place at heart place in the cells, which as the scalar value of velocity potential and temperature wherein.Once velocity potential is determined, then speed can be determined from the approximate form of the equation write out for discretize computing grid (2).

Calculation of pressure is not needed when using PFM calculated gas flow.But, for coupling pressure be may be used for carrying out the specific flow boundary condition of modeling, such as soaps based on pressure to the method that air-flow is predicted.Use these class methods of Bernoulli equation to describe being numbered in the U.S. Patent application of PCT/US2011/051866 of " SYSTEM AND METHOD FOR PREDICTING PERFORATED TILEAIRFLOW IN A DATA CENTER " that be entitled as that on September 16th, 2011 submits to, the full content of this patented claim is incorporated to herein by reference.Once gas velocity is determined, then temperature can use energy equation (4) to determine.

$\stackrel{\→}{V}\·\▿T={\mathrm{\α}}_t{\▿}^{2}T---\left(4\right)$

Wherein, α _t=k/ (Uc _p) be thermal diffusivity.In practice, can be left in the basket when being conducive to the convective term dominated at the diffusion term of the right-hand side of equation (4); But, retain diffusion term and with the addition of additional degree α _t, this degree of freedom can be " tuned ' the precision of the prediction affecting PFM.

Physical problem spatial division to be analyzed is some grid cells by above-mentioned PFM method.In some examples, structuring flute card grid is adopted.In other examples, grid is not structured.For automatically generate unstructured grid and in unstructured grid further interpolation field value method and system PCT/US2012/030373 application in discuss.

enhancement mode PFM

As mentioned above, enhancement mode PFM is the method being designed to the precision improving the PFM emulation discussed in PCT/US2012/030373 application.PFM is by adding one or more bearing calibration to strengthen to basic PFM method, haloing unit speed method and the buoyancy correction of weighing method comprising Jet model described below.In the embodiments described herein, the flow pattern that the physical property that Jet model and haloing unit speed method apply potential barrier is not caught.By specifying fixed flow to move in PFM model, jet and bouyant plume can produce in PFM.On the contrary, in the embodiments described herein, the calculating of the buoyancy at each unit place is coupled realize mass balance and realize appropriate air-flow and temperature field further with innovation mapping step by buoyancy correction of weighing method.Should be appreciated that Jet model and haloing unit speed method can be implemented before solution initial p FM model.In addition, the mapping step be associated with buoyancy correction of weighing may need the solution of another prescription journey except PFM model, as described below.

In some embodiment as herein described, enhancement mode PFM uses the speed determined by rule of thumb as internal boundary condition.Traditionally, speed is only designated in fluid-solid interface or fixes, and all speed in internal fluid are calculated.Such as, when using classic method, the speed at accommodating lattice, equipment rack air intake opening or gas outlet place can be specified in.On the contrary, in each embodiment as herein described, the computer system of enhancement mode PFM method is used to specify the internal speed be associated with haloing unit and Jet model to comprise the physical influence of not intrinsic modeling in PFM by rule of thumb, the such as effect of jet, thermal blooming and buoyancy.

Shown in Figure 3 according to the general introduction of the E-PFM method 300 of an embodiment.According to an example, the computer system of all data centers described above and management system 106 implements E-PFM method.Exemplary method 300 comprises whole three kinds of bearing calibrations, comprises Jet model correction, the correction of haloing unit and buoyancy correction of weighing.Should be appreciated that method 300 comprises generalized case, and not all three (jet, haloing unit and buoyancy) correction can be implemented for each application.Particularly, in one embodiment, computer system can only perform the bearing calibration of jet air-flow.Such as, user may want that modeling is ventilated in the room of buildings in general application.In another embodiment, computer system can only perform the bearing calibration of haloing unit.In yet another embodiment, computer system only can perform buoyancy correction of weighing method.Should be appreciated that computer system can perform any combination of jet, haloing unit and buoyancy method.

In one embodiment, computer system can determine whether to need execution bearing calibration or which bearing calibration needs to be performed.In one example, user's input information is in the interface 104 presented by design and management system 106, and this information specifies the condition needing the performance of all methods or either method wherein.Such as, user can designated computer system thus the object in facility is described, it plays the effect of heat producer, such as computing machine and human body.In this illustration, computer system can perform haloing bearing calibration to illustrate that these produce the existence of the object of heat.

In other examples, computer system can be determined the presence or absence of condition and implement corresponding bearing calibration.Such as, computer system can by determining the existence of the soaps of isolating and determining that Jet model bearing calibration needs to be performed about the information of data center layouts.

Method 300 can also comprise the step receiving and represent the information of the physical layout of identified region.This layout information can represent and the various characteristic that this region and the object be contained in wherein are associated.Among other things, the example of these characteristics comprises the size in region, spaced portions (being referred to herein as blocking part) in this region, object (such as desk, chair, computing machine, human body) in this region and the soaps belonged in the information of equipment rack, the Information Technology Equipment adopted in this equipment rack, CRAC, UPS, PDU, raised floor characteristic, room and row's formula cooling device.Represent that the information of physical layout can be input in the interface 104 presented by design and management system 106 by user, or accessed from the memory storage comprising the information belonging to this layout by computer system.

In step 302, if jet-like flowing is determined to exist, then effluxvelocity is applied in flow field by computer system together with other boundary conditions.As mentioned above, jet-like flowing can be determined to exist based on the physical layout of buildings or based on the input of user.In some examples, jet-like flowing is associated with supply air diffuser, free-standing frame and other air-flows supplier.Boundary condition can comprise room air inlet and go out air-flow, frame air outlet and exhaust outlet and other borders.As above noted, the jet gas velocity that use jet air-flow as described below bearing calibration is determined can specify before use initial p FM model determination gas velocity angle value.In at least some example, jet gas velocity can previously not specified or be designated on undetermined value.In other examples, jet gas velocity can substitute zero gas velocity angle value.In other other examples, if gas velocity angle value has used additive method to determine, so the bearing calibration of jet air-flow can substitute the gas velocity angle value using those methods to determine.

In step 304, speed is applied in the haloing unit of proximity thermal volume production biosome by computer system on the velocity amplitude using initial p FM model to determine.According to some embodiment, haloing unit and velocity amplitude really fix on and describe below with reference to Fig. 4 and Fig. 5.In one example, haloing unit comprises the computing unit of contiguous solid body.Speed in computer system can be determined and specify (fixing) haloing unit is with to the thermal drivers bouyant plume modeling around solid body.In some examples, solid body can comprise be positioned at facility personnel, computing machine and data center apparatus.The haloing speed of specifying can be determined by analysis existing in document or empirical equation, or as described belowly can determine by rule of thumb from CFD.

In one example, the gas velocity that use haloing unit as described below bearing calibration determines haloing unit can specify before use initial p FM model determination gas velocity angle value.In at least some example, can not specified by the gas velocity determining haloing unit or be designated on non-determined value.In other examples, zero gas velocity angle value can be substituted by the gas velocity determining haloing unit.In other other examples, if gas velocity angle value has used additive method to determine, so the bearing calibration of haloing unit air-flow can substitute and use the determined gas velocity angle value of those methods.

In one embodiment, each grid cell has velocity amplitude on each face; Four masks in 2D model have four speed, and six masks in 3D model have six speed.Corrected by haloing, at least one embodiment, computer system specifies the speed component alignd with bouyant plume (vertical direction).By jet calibration method, computer system specifies the speed component orthogonal with supply-side (horizontal direction).Other components of speed can be in the state of " freedom ", to use PFM method to calculate within step 306 by computer system.These " freedom " components of speed can be revised potentially by mapping step 312, as described below.Within step 306, computer system performs PFM method to obtain airflow field, as described also as outlined above in PCT/US2012/030373 application.In one example, first computer system solves the velocity potential equation of grid cell, until meet the mass conservation error threshold of specifying, this mass conservation error threshold is the function of total room airflow.Once velocity potential field is determined, then next computer system determines speed.

In step 308, computer system serviceability temperature solver determination temperature value, as described in PCT/US2012/030373 application.In one embodiment, computer system uses the grid cell boundary speed calculated to carry out the temperature equation of iterative grid cell, until meet the mass conservation error threshold of specifying, this mass conservation error threshold is the function of the gross energy adding room to.Alternatively, computer system can use the grid cell boundary speed calculated to carry out the temperature equation of iterative grid cell, until the difference between frame load and chiller circuit load meets the specification error threshold value of the gross energy adding room to.In at least some example, determine that temperature value can be optional.

In the step 310, computer system application buoyancy correction of weighing method, the method adds additional speed based on local unit temperature to each inner mesh unit.Buoyancy correction of weighing in the step 310 can the mass conservation character in disturbance flow field.Therefore, in step 312, mapping step is performed with correction mass unbalance.With jet air-flow and haloing unit bearing calibration complementation, the computer system of buoyancy correction of weighing method is used to regulate the flow field in each place except jet or the appointed region of haloing unit speed.And, be different from other two kinds of alignment techniques, the speed that buoyancy gas velocity is added to (instead of substituting) uses PFM method to be predicted by computer system.Computer system in a step 314 based on the determined gas velocity angle value of use buoyancy correction of weighing method, can optionally determine the temperature of each grid cell.

Finally, in step 316, the multiple buoyancy correction of weighing of computer system application, wherein mapping step 312 is followed in each buoyancy correction of weighing, until obtain mass conservation flow field.Buoyancy correction of weighing step 310, temperature solver 312 and mapping step 314 can be secondary arbitrarily by repetition.In some examples, step 310,312 and 314 can by repeatedly setting number of times, and one of them step provides minimum computing time.

According to some example, the computer system of all designs as described above and management system 106 can use 2D or 3D visualization tool to implement the process determining and show above determined air-flow and temperature value in method 300.According at least one example, show in the 2D cross-sectional area that air temperature and current value can be extended by the region in the room in such as data center or buildings.In other examples, show, as further described in PCT/US2012/030373 application in the 2D plane that the region in air temperature and current value can represent at the 3D crossing over region extends.

the jet air-flow of specifying

Effluxvelocity from the jet-like flowing determined in step 302 describes with reference to figure 4 and Fig. 5 below.Fig. 4 illustrates the notional difference between jet air-flow in idealized room 400 and radial air flow.Air can pass through various types of air intake opening (such as, such as ceiling air diffuser device and grid) and be supplied to room.

Due to momentum, the flow pattern close to air intake opening and (or air outlet) is not merely close to flowing wherein by the reversion of gas outlet (or return air inlet) of extracting out from all directions quite equably.On the contrary, air inlet stream 402 is jet-like, as shown in the jet 406 in Fig. 4.In one example, jet-like air-flow be along the speed of the center line of axis of flow 408 be high and quite constant air-flow.In addition, flow pattern leap jet boundary is only minimally dispersed and keeps speed sharply to rise.Above-mentioned PFM method does not directly illustrate the momentum of air-flow, and therefore produces radial kenel 410 at air intake opening 404 place, as shown in Figure 4.Jet-like air-flow can be seen usually in buildings in general application, and in some examples, also can observe by the heart in the data.Jet air-flow bearing calibration as described below just carrys out the jet-like of regulation amendment by rule of thumb flow pattern at the air-flow in air supply downstream by specifying in result emulation.

Fig. 5 represents the display in the idealized room 500 in working space, and it comprises various physical objects, such as desk, human body and computing machine.Air-flow is entered by air intake opening 502, and can be removed by return air inlet 504 at least partly.In one example, region 506 represents the zone of influence be associated with gas flow jet.This zone of influence comprises some unit, and each unit has speed associated with it.By using one group of jet formula, the speed being parallel to the unit in the zone of influence of jet direction in whole jet can be determined.In one example, the zone of influence can use the unit of 20% of the determined initial jets speed of PFM model to determine by determining that speed is greater than wherein.Use the velocity amplitude in the determined zone of influence of a kind of model in Jet model described below can be applied in the emulation of numerical value air-flow.

In each example, predefined jet boundary unit can be considered as the fixed speed border along main flow moving axis by computer system, and calculates residue flow field subsequently to generate the airflow field of the mass conservation.If jet airflow collision, the effluxvelocity of so specifying may be overlapping in some unit.Speed can each other at superimposed on top and consequent speed can by simple algebraic manipulation.Residue flow field automatically caught subsequently by PFM solver by implementation quality conservation.

Jet-like flowing can be used in spendable one or more model in E-PFM method and be specified by computer system.In an embodiment (such as floor/ceiling plenum chamber (floor/ceilingplenum) application), computer system uses 2D Jet model.2D jet is modeled as the line source of the air-flow of the virtual point be derived from outside border, territory, as shown in Figure 4.2D jet flow pattern as shown in Figure 5 can use subsequently by A., Bejan 1984 John Willie father and son publishing company in New York deliver " convection heat transfer' heat-transfer by convection" in the standard method that describes analyze and determine.This process comprises and utilizes the reduced form being suitable for fluidic features to write out the 2D equation of momentum.This equation can solve to produce Jet Axis gas velocity with closed form.Can determine from known Jet Axis speed and the mass conservation perpendicular to the gas velocity on the direction of axis of flow.

In another embodiment, in order to modeling 3D room (office, data center etc.), computer system can use 3D Jet model.The isothermal rotational symmetry 3D Jet model simplified is most realistic, and in the model, jet supply temperature equals ambient room temperature, and effluxvelocity distribution plan is symmetrical along jet direction.3D Jet model is by Huo, Y., Zhang, J., Shaw, C. and in Haghighat, F. " the A New Method to Describe the Diffuser Boundary Conditions in CFDSimulation " that deliver at the 233-240 page of the volume two of the Proc.of ROOMVENT ' 96 of 1996 be described.

In a further embodiment, when jet be attached to surface and jet air-flow by neighbouring barrier (such as, wall) affect time, " attached jet " model should be used.Calculate be used for isothermal attached jet speed method by Verhoff, A. in 1963 the report 626 of Princeton University deliver " the Two-dimensional Turbulent Wall Jet with and without an External stream" in be described.

In yet another embodiment, computer system can obtain effective stream from the empirical model generated the CFD emulation of various practical application.In order to train the empirical model for jet-like air-flow, emulation can perform at different conditions, such as mass rate of emission, jet size, jet direction, close to wall and other objects etc.Along the velocity field of jet direction by monitored, and come for generating and training experience model by the such as tuning coefficient be associated with velocity field.This empirical model can be merged in PFM solver to specify in the effluxvelocity distribution plan in computational fields subsequently, is similar to other above-mentioned Jet models.

the haloing unit speed of specifying

Determined haloing unit speed describes with reference to figure 6 and Fig. 7 below in step 304.The computer system of use PFM method described in PCT/US2012/030373 application can not comprise by near the heat that generates of object drive the Thermal plume of rising.These objects can comprise any object of such as computing machine, personnel, PDU and other objects.In the embodiments described herein, in order to increase the precision of air-flow solver, haloing model of element is used for being incorporated to buoyancy effects.

Should be appreciated that the surrounding air heated by object is risen due to buoyancy, and the addition speed therefore produced close to air-solid border increases (or " kick (kick) ").In order to modeling Thermal plume, introduce the concept of haloing unit.In one example, haloing unit can comprise the computing unit be close to around thermal objects.The speed of all haloing unit can rule of thumb calculate by formula, and vertical speed field (aliging with buoyancy) can be defined as boundary condition before PFM solver starts.Fig. 6 A and 6B illustrates object 600, and it is represented as heating block and comprises haloing unit 602 and 604.Fig. 6 A illustrates the haloing unit 602 of the uprush of not specifying, and Fig. 6 B illustrates the haloing unit 604 with the vertical gas flow field of specifying.Each in haloing unit 604 has position, and each in vertical gas flow field has speed V _h.In one example, according to the distance to floor, the position of haloing unit is defined by the height of this unit (y).

Be similar to the jet air-flow calculated above, use the computer system of PFM method haloing unit speed can be considered as the speed edges of fixing.Use the computer system of haloing unit speed method to guarantee that correct y-velocity field is stored near heating block, and use PFM solver to calculate remaining flow field subsequently.As Jet model, the speed of specifying " around " PFM solution automatically produce appropriate mass conservation field.

In order to determine the model of the speed " kick " through determining haloing unit, new empirical model is defined and uses CFD emulation for effective modeling and the buoyancy effect of training in haloing unit.Usually, CFD model comprises the thermal objects being placed in vacant room center, and this vacant room is unlimited or has fully away from the symmetrical border of this object.Room-size and ceiling height are changed the characteristic of catching different room configuration.Difformity and power level are also modeled to emulate dissimilar heat and generate object, such as computing machine, human body or PDU.

Speed around heating block in haloing unit is monitored.Fig. 7 illustrates an example of the uprush specification in the haloing unit of Fig. 6 B.In example shown in Figure 7, representational haloing unit speed is the function close to height (y-axle) on the function of the position (A, B or C position) of the human body in horizontal plane or floor.The example expression calculating haloing unit speed can be defined as follows:

$V_h=C_1{y}^{C_2}---\left(5\right)$

Wherein, V _h(m/s) be haloing unit speed, y (m) is the distance to floor, C ₁and C ₂be constant, this constant changes according to the configuration in room and the characteristic of thermal objects.Such as, (there is 1.8m for the human body represented in Typical office room ²surface area) heating block, constant C ₁and C ₂may be calculated:

C ₁＝0.0003×P+0.047

C ₂＝0.62

Wherein, P is the power (W) that heating block produces.

As shown in Figure 7 A, in one example, haloing unit speed distribution plan comprises the multiple expressions 702,704 and 706 corresponded to around the diverse location (being labeled as A, B, C equally) in the face of heating block.The 3D of heating block represents and illustrates in figure 7b, and the 2D of heating block comprising face A, B, C represents and illustrates in fig. 7 c.In addition, haloing unit speed distribution plan also may depend on grid cell size.Such as, if PFM method uses 6 haloing unit, so, 6 unit should also be utilized from the corresponding CFD model wherein deriving haloing speed.And the speed at heating block top may be relatively uniform, and therefore can use identical formula to estimate.It is noted that the speed of calculating can be designated in all haloing unit or the haloing unit be only top dog at buoyancy.

buoyancy correction of weighing

The buoyancy correction of weighing determined in the step 310 describes with reference to figure 8A-8C below.As mentioned above, buoyancy correction of weighing can illustrate and the speed that warm air rises and cold air sinks to being associated.Generally speaking, buoyancy correction of weighing method uses PFM method calculated gas flow and temperature field, increases buoyancy speed v subsequently to each grid cell in vertical direction _b.Correct with above-mentioned jet air-flow and compare with the bearing calibration of haloing unit, in the buoyancy correction of weighing method be described below, whole flow field not only regulates neighbouring heating/cooling object, and regulates air supply.And, be different from other two kinds of alignment techniques, the speed that buoyancy correction of weighing is added to (instead of substituting) is predicted by PFM.

In one example, the speed be associated with buoyancy correction of weighing can be represented as:

Wherein, g is the acceleration because gravity produces,

β is the coefficient of volume thermal expansion,

H is feature vertical length yardstick,

Δ T is that the characteristic temperature between interested point (grid cell) and another reference value is poor, and

α is the coefficient determined by rule of thumb.

If numerical value is positive, then " symbol (Sign) " function returns 1, if numerical value is 0, this sign function returns 0, and if numerical value is negative, this sign function returns 1.Its for guarantee on reference temperature and under local temperature produce buoyancy speed up or down respectively.

In at least one example, β can equal the inverse of the absolute temperature of ideal gas, and side reaction coefficient may be used for " tuning " E-PFM and predicts with optimum matching CFD.

In some examples, v _b's result partly from simple dimensional analysis, wherein suppose that vertical momentum and buoyancy balance each other.Therefore, in one example, experience factor α can be 1 rank.It is noted that " feature " value H and Δ T can select arbitrarily to a certain extent with Optimized model accuracy.Such as, H can be interpreted as the height in room or it can change based on the height of neighbouring solid body or it can be taken as the height of each individual grid cell.Similarly, Δ T can be selected as fixing reference temperature or select with reference to adjacent grid cell.In one example, fixing reference temperature can be the temperature of the most cold air being supplied to room.

Because temperature in the cells, which heart place is calculated and speed is calculated at elemental area place, so the speed at any horizontal cell face place is based on the average v as calculated for the unit just under this elemental area and just on this elemental area _b.Derive from the speed v being applied to each grid cell _bthe flow field almost always no longer mass conservation, and therefore no longer energy conservation, thus produce power is unbalance.Therefore, " mapping " step described below may be used for the flow field that the buoyancy determined in the step 310 has corrected to be converted back to appropriate mass conservation field.

mapping step corrects

As mentioned above, mapping step is the program for the flow field determined in above-mentioned steps 312 being converted to mass conservation flow field.This mapping step is initially developed to time splitting technology and has been used to rapid fluid dynamics (FFD) emulation subsequently.In both cases, this step, simultaneously must the various physical phenomenon of modeling for implementing the mass conservation of other solution steps producing mass-unbalance.Mapping these initial usages of (corrections) step is purely apply for transient state (time change).E-PFM in example described herein analyzes normally static Simulation, and therefore, final mapping step had nothing to do with the time.

The detailed derivation of the stable state version of mapping step is not shown here, but this step starts with transient analysis, in this transient analysis, if we will flow " pushing back " for the mass conservation by means of only the pressure effect in little time step, so we will calculate the pressure field that must exist.Finally, if our pressure after the modification interior " hiding " time step, then do not need the reference time.

In order to correction is applied to steady-state gas flow form, we are first based on the non-mass conservation velocity field by producing to flow field interpolation buoyancy correction of weighing calculate amended pressure.As the velocity potential in analyzing at initial p FM, amended pressure can by solving Poisson equation to calculate:

${\▿}^2\tilde{p}=\▿\·\stackrel{\→}{V}*---\left(7\right)$

In order to obtain the speed of correction of a final proof from mapping step, we deduct the gradient of the pressure field of amendment from the flow field of (non-mass conservation) buoyancy correction of weighing:

$\stackrel{\→}{V}=\stackrel{\→}{V}*-\▿\tilde{p}---\left(8\right)$

The mapping step method producing mass conservation field from non-mass conservation field is summarized in fig. 8.For simplicity, our discrete form of equation of presenting with reference to the structuring computing grid of the 2D of figure 2.

In step 802, the pressure of amendment gradient can be calculated as from non-mass conservation flow field:

${(\▿\·\stackrel{\→}{V}*)}_i=\frac{1}{\mathrm{\Δ}x}(u_E^{*}-u_W^{*}+v_N^{*}-v_S^{*})---\left(9\right)$

In step 804, we calculate the pressure of amendment by solving Poisson equation.With reference to the computing grid of figure 2, this computing grid comprises unit 202,204,206 and 208, amendment calculate at grid cell i place, wherein represent the pressure of the amendment in unit 202, represent the pressure in unit 206, represent pressure in unit 204 and represent the pressure in 208.Amendment can be represented as:

${\tilde{p}}_i=\frac{1}{4}({\tilde{p}}_N+{\tilde{p}}_S+{\tilde{p}}_E+{\tilde{p}}_W-{\mathrm{\Δx}}^2{(\▿\·\stackrel{\→}{V}*)}_i)---\left(10\right)$

Finally, in step 806, we come calculated mass conservation flow field by the gradient deducting the pressure field of amendment from initial non-conservation flow field.

In step 806, actual speed can be confirmed as the function of the pressure revised in each cell:

$u_E=u_E^{*}-\frac{1}{\mathrm{\Δ}x}({\tilde{p}}_E-{\tilde{p}}_i)$ With $u_W=u_W^{*}-\frac{1}{\mathrm{\Δ}x}({\tilde{p}}_i-{\tilde{p}}_W)---\left(11\right)$

$v_N=v_N^{*}-\frac{1}{\mathrm{\Δ}x}({\tilde{p}}_N-{\tilde{p}}_i)$ With $v_S=v_S^{*}-\frac{1}{\mathrm{\Δ}x}({\tilde{p}}_i-{\tilde{p}}_S)---\left(12\right)$

Fig. 9 A-9C illustrates an example of mapping step, this step be included in mapping step be performed after initial non-mass conservation field and mass conservation field.Fig. 9 A illustrates boundary condition represented in the 2D in room 900 represents.As shown in the figure, air-flow enters the lower left quarter in room 900 and leaves from upper right quarter.Initial fields is set to roughly northeast stream, the exhausr port as shown in Figure 9 B from the air intake opening of lower-left to upper right.In figures 9 b and 9, this is flowing in inner pure just northeastward, because the initial fields of the circumference in flowing and wall edge effect territory.Fig. 9 C illustrates the flow field of revising after mapping step.In Fig. 9 C, the principal character of airflow field is retained, and flowing is still in the northeastward haply from air intake opening to exhausr port.In addition, air-flow also comprises seamlessly transitting of the sudden change replaced in the flowing direction now.Airflow field shown in Fig. 9 C is consistent with the mass conservation.

the example in general room

The result of E-PFM method 300 can illustrate with reference to Figure 10,11 and 12.Figure 10,11 and 12 compares three kinds of methods, sign Jet model described in this article, in the speed of specifying at haloing unit place and CFD, PFM and E-PFM method of buoyancy correction of weighing.These methods are implemented, because they are by general constructure ventilation design ap-plication.CFD, PFM and E-PFM method in the following example performs and comprises the simple coarse grid model in room under similar emulation setting and input.Figure 10 illustrates an example in the room 120 in working space environment, and this room has 10ft and is multiplied by the size that 10ft is multiplied by 10ft.Air-flow is input to the ceiling return air inlet 124 that ventilation orifice 122 in room and air are left from room by it by the stage place that working space 120 is included in a wall.Room 120 comprises a heat and generates block 126, such as, adds the human figure object of 50W to room.Room 120 also comprises desk 128 and computing machine 130.

Figure 11 A-11C illustrates that the 2D in the temperature field of the centre (such as at 5 feet of places) by room 120 represents.The temperature field generated by CFD is shown in Figure 11 A, and PFM temperature field is shown in Figure 11 B, and the temperature field that E-PFM produces is shown in Figure 11 C.The synthesis temperature determined at the use PFM shown in Figure 11 B is strongly obtained overpredict in the region on the right side of heat generation block 134, and rises in the temperature of also not observing from bottom to top from low to high in this room.The bearing calibration of E-PFM shown in Figure 11 solves the temperature of these excessive temperature prediction and produces hygrogram true to nature.

Except C.T. field, the maximum temperature for the room of each method also can be compared.This maximum temperature can help to assess the ability of our correct level based on each the collection present mixt in the method for PFM compared with CFD.Use the maximum temperature calculated based on PFM method to be 44 DEG C for room 120, and be 21 DEG C for the maximum temperature of CFD, and be 22 DEG C for the maximum temperature of the method based on E-PFM.The maximum temperature of E-EFM is similar to CFD, demonstrates E-PFM method and provides remarkable improvement more than PFM.

Figure 12 A-C illustrates that the 2D of the velocity field of the centre by room 120 represents, wherein, the velocity field generated by CFD illustrates in fig. 12, and the velocity field generated by PFM is shown in Figure 12 B, and the velocity field generated by E-PFM illustrates in fig. 12 c.Velocity field in a room provides the explanation of the gain of the accuracy for E-PFM.Figure 12 A-C illustrates speed in the plane identical with Figure 11 A-C or section and direction.The speed that the use E-PFM illustrated in fig. 12 c produces is included in two recirculation regions 132 and 134 existed in CFD.Recirculation regions 132 is decided to be the left side of heating block, and recirculation regions 134 is decided to be the right side of heating block.The speed that use PFM shown in Figure 12 B produces does not comprise recirculation regions.Recirculation regions 132 and 134 mixed airflow, produces some layering, and prevents the viscous flow of the excessive temperature utilizing PFM to cause seen in Figure 11 B.

data-center applications and example

The result of E-PFM method 300 can illustrate by reference diagram 13-15, the method in the data heart air temperature and current estimate in PFM and E-PFM more as herein described.As mentioned above, haloing unit may not be relevant to data-center applications with Jet model.Therefore, in the example following E-PFM method, the method only implements buoyancy correction of weighing.But, in these examples, even if only have buoyancy correction of weighing, also gain quite accurately can be realized.

Figure 13 illustrates the six kinds of layout A-F characterized in the analysis be described below according to some example.Layout A-F in example described herein provides the integrating representation of the performance of E-PFM.This layout by Healey, C., VanGilder, J., Sheffer, Z. and Zhang, " the PFM and CFDcomparison study that X. delivered at " the Proceedings of InterPACK " of Portland, Oregon in 6-8 day in July, 2011 potential-Flow Modeling for Data Center Applications" in be described.Layout A-F provide relative to raised floor and Local cooling heating data center room in maximum frame air inlet temperature and average frame air inlet temperature prediction E-PFM quantification improve.This layout table is shown in the various configurations of frame represented in the configuration of various passage.Such as, layout A comprises that multiple two passages not arranging formula refrigeratory are trooped, four CRAC refrigeratorys and soaps, and layout E comprises multiple two passages troops, each trooping is separated by the passage of heat, it is characterized in that not boring a hole floor tile and only having row's formula refrigeratory.

In order to simplify the comparison between E-PFM, PFM and CFD, have employed accurate tolerance, as by Healey, C., VanGilder, J., Sheffer, and to discuss in Zhang, X. " the Potential-FlowModeling for Data Center Applications " that deliver at " the Proceedings of InterPACK " of Portland, Oregon in 6-8 day in July, 2011 Z..The E-PFM of given frame j estimates temperature relative to CFD estimated value accuracy be defined as:

Wherein, Δ T _refreference temperature difference (being set to 10 DEG C in this example).PFM can define equally relative to the accuracy of CFD.

Figure 14 A and 14B illustrates that E-PFM is to the comparison of CFD in layout C, and this compares the bat comprised for E-PFM and PFM, changes the factor alpha determined by rule of thumb in equation (6) simultaneously.As mentioned above, in one example, the function of coefficient that the speed be associated with buoyancy correction of weighing can be represented as acceleration because gravity produces, the coefficient of volume thermal expansion, feature vertical length yardstick, characteristic temperature difference between interested point (grid cell) and another reference value and determine by rule of thumb.In this illustration, factor alpha such as { is changing between multiple values of 0.15,0.7,1.4,2.8,7.0,14}.The feature vertical length yardstick H of equation (6) is considered the height of each unit, is 0.15 meter in this illustration.

Figure 14 A is depicted in the accuracy of E-PFM and PFM be averaged on the mean air entry mouth temperature valuation in layout C, and Figure 14 B is depicted in the accuracy of E-PFM and PFM be averaged on the full admission mouth temperature valuation of layout C.Figure 14 A and 14B illustrates that E-PFM method is by appropriately selecting α. significantly can reduce error, thus improve accuracy and preferably divide remainder error (when having one group of perfect valuation, producing the accuracy of 1) equally.This drawing illustrates the essential improvement in average gas port temperature prediction and full admission mouth temperature prediction, especially the latter, and it makes a special effort to strive for obtaining the field of accuracy for PFM.Although should be appreciated that in all values of the α between 0.15 and 14 to there is accuracy gain, comprise utilize other layout A, the most effective value of results suggestion α of other experiments that B, D, E and F perform can drop between 1 and 5.

Figure 15 illustrates the comparison of E-PFM and PFM accuracy firm for mean air entry mouth on all six A-F layouts.Figure 15 is determined by E-PFM, wherein u=1.4 and H=0.15.In four layouts in six layouts, compare with PFM, observed essential improvement (error more than 33% reduces).Worsen slightly although observe to have for the accuracy of mean air entry mouth temperature E-PFM in layout B, it should be understood that u=1.4 still improves the accuracy of maximum temperature prediction in this layout.In executed all tests, have been found that layout E is relatively independent of the enforcement of E-PFM, when changing the value of factor alpha, accuracy had not both improved and had not declined yet.

These results are advised, use the E-PFM with selected α will to improve the accuracy of data center's temperature prediction significantly, in general, improve the accuracy of average and maximum frame air inlet temperature.In the tested data center case exceeding half, the use of buoyancy correction of weighing significantly reduces error.

computer system

The various aspects described in this article according to the present embodiment and function may be implemented as hardware in one or more computer system or software.There are many examples of the computer system used at present.Except other computer systems, these examples comprise network home appliance, personal computer, workstation, mainframe, networked clients, server, media server, application server, database server and web page server.Other examples of computer system can comprise mobile computing device (such as cell phone and personal digital assistant) and the network equipment (such as load equalizer, router and switch).In addition, the many aspects according to the present embodiment can be positioned in single-computer system, maybe can be distributed in be connected to one or more communication network multiple computer systems among.

Such as, various aspects and function can be distributed in one or more computer system, and it is configured to the part execution overall tasks providing service to one or more client computer or be configured to as distributed system.In addition, these aspects can perform on client-server or multiechelon system, and it is included in the assembly performing and distribute between one or more server systems of various function.Therefore, embodiment is not limited to perform on any particular system or system group.In addition, these aspects can be implemented in software, hardware or firmware or its any combination.Therefore, the aspect according to the present embodiment can be implemented in the assembly of the various hardware and software configuration of method, action, system, system element and use, and embodiment is not limited to any specific distributed architecture, network or communication protocol.

Figure 16 illustrates the block diagram of Distributed Computer System 160, within the system, can put into practice according to the various aspects of the present embodiment and function.Distributed Computer System 160 can comprise one or more computer system.Such as, as shown in the figure, Distributed Computer System 160 comprises computer system 162,164 and 166.As shown in the figure, computer system 162,164 and 166 is interconnected by communication network 168, and can exchange data by this communication network 168.Network 168 can comprise any communication network that computer system can exchange data whereby.Data are exchanged in order to use network 168, computer system 162,164 and 166 and network 168 can use various method, protocol and standard, except additive method, protocol and standard, it comprises token ring, Ethernet, wireless ethernet, bluetooth, TCP/IP, UDP, Http, FTP, SNMP, SMS, MMS, SS7, Json, Soap and Corba.In order to ensure data transmission be safe, computer system 162,164 and 166 use except using other safety techniques comprise TLS, SSL or VPN multiple safety measures via network 168 to transmit data.Although Distributed Computer System 160 illustrates three computer systems of networking, Distributed Computer System 160 can comprise and uses any amount of department of computer science of any medium and communication protocol networking to unify calculation element.

The specialized hardware or software that perform in one or more computer systems of the computer system 162 comprised as shown in figure 16 is may be implemented as according to the various aspects of the present embodiment and function.As shown in the figure, computer system 162 comprises processor 170, storer 172, bus 174, interface 176 and storer 178.Processor 170 can perform a series of instructions producing service data.Processor 170 can be commercial processor, such as intel pentium, Motorola PowerPC, SGI MIPS, Sun UltraSPARC or Hewlett-Packard PA-RISC processor, but when obtaining other processors many and controller, processor 170 can be the processor of any type, multiprocessor, microprocessor or controller.Processor 170 is connected to other system element by bus 174, and this system element comprises one or more storage arrangement 172.

Storer 172 may be used for storage program and data during the operation of computer system 162.Therefore, storer 172 can be the random access memory of relative high-performance, volatibility, such as dynamic RAM (DRAM) or static memory (SRAM).But, storer 172 can comprise any device for storing data, such as disc driver or other memory storages that is non-volatile, non-transitory.Can storer 172 be organized as special according to each embodiment of the present invention, and in some cases, be organized as the unique texture performing aspect disclosed herein and function.

The assembly of computer system 162 can be coupled by the interconnection element of such as bus 174.Bus 174 can comprise one or more physical bus, such as, be integrated into the bus between the components in same machines, but the communications connector that can also be included between system element, it comprises calculating bussing technique that is special or standard, such as IDE, SCSI, PCI and InfiniBand.Therefore, bus 174 makes communication (such as, data and instruction) can exchange between the system component of computer system 162.

Computer system 162 also comprises one or more interface arrangement 176, such as the combination of input media, output unit and input/output device.Interface arrangement can receive input or provide output.More specifically, output unit can present the information for outside display.Input media can receive information from external source.The example of interface arrangement comprises keyboard, mouse apparatus, trace ball, microphone, touch-screen, printing equipment, display screen, loudspeaker, network interface unit etc.Interface arrangement allows computer system 162 with external entity (such as user and other system) exchange information and communicates with.

Storage system 178 can comprise storage medium that is computer-readable and writable, non-volatile, non-transitory, and the instruction definition stored in the medium treats the program performed by processor.Storage system 178 can also comprise the information be recorded on medium or in medium, and this information can by routine processes.More specifically, information can be stored in one or more data structure, and this data structure is configured to preserve storage space or increase exchanges data performance particularly.This instruction can be persistently stored as coded signal, and this instruction can cause processor to perform any function as herein described.Except other media, medium can also be such as CD, disk or flash memories.When operating, processor or some other controller can make data be read another storer (such as storer 172) from nonvolatile recording medium, compared to the storage medium be included in storage system 178, this storer allows processor access information more quickly.Storer can be positioned in storage system 178 or storer 172, and but, processor 170 can data in operational store 172, and subsequently after processing is completed, can copy data to the medium be associated with storage system 178.Multiple assembly can manage the data mobile between medium and integrated circuit memory element, and current described embodiment is not limited thereto.In addition, embodiment is not limited to specific accumulator system or data-storage system.

Although computer system 162 is illustrated as the computer system of a type by way of example, can be put into practice according to the various aspects of the present embodiment and function by this computer system, but any aspect of current the disclosed embodiments is not limited in computer system shown in Figure 16 to be implemented.Various aspects according to current the disclosed embodiments can be put into practice with function on one or more computing machines with the architecture different from the architecture shown in Figure 16 or assembly or assembly.Such as, computer system 162 can comprise the specialized hardware of special programming, is such as such as customized to the special IC (ASIC) being applicable to performing specific operation disclosed herein.Although another embodiment can use the several general-purpose calculating appts running MAC OS X with 32 Intel and 64 Intel processors and the several dedicated computing devices running proprietary hardware and operating system to perform identical function.

Computer system 162 can be the computer system comprising operating system, and this operating system management is included in hardware element in computer system 162 at least partially.Usually, processor or controller (such as processor 170) executive operating system, this operating system can be such as based on the operating system of form (Windows), the Windows NT that such as can obtain from Microsoft, Windows2000 (Windows ME), Windows XP, Windows Vista, Windows 7 and Windows8 operating system, the MAC OS System X operating system that can obtain from Apple Computer, one in the many distribution of the operating system based on Linux versions, such as, enterprise's (SuSE) Linux OS that can obtain from Hong Mao company, the solaris operating system that can obtain from Sun Microsystems (Sun Microsystems) company, or from the UNIX operating system that various source obtains.Other operating systems many can be used, and embodiment is not limited to any particular implementation.

Processor defines computer platform together with operating system, can write application program for this computer platform with high-level programming language.These component application can be executable, middle, and such as, C-, syllabified code or coding and decoding, it uses communication protocol (such as TCP/IP) at communication network (such as internet) enterprising Serial Communication.Similarly, Object-Oriented Programming Language (such as .Net, SmallTalk, Java, C++, Ada or C# (C-Sharp)) can be used to implement according to the aspect of current the disclosed embodiments.Also other Object-Oriented Programming Languages can be used.Alternatively, functional expression, script or logic programming language can be used.

In addition, can implement in non-program environment according to the various aspects of current disclosed embodiment and function, such as, the document formed with HTML, XML or extended formatting, when checking the document in the window at browser program, the aspect of render graphics user interface or perform other functions.In addition, each embodiment according to the present invention may be implemented as the element of programming or non-programmed, or its any combination.Such as, webpage can use HTML to implement, and the data object called in webpage can be write with C++.Therefore, current the disclosed embodiments are not limited to certain programmed language, and also can use any suitable programming language.

The computer system be included in embodiment can perform the other function outside the scope of current the disclosed embodiments.Such as, the aspect of system can use existing commercial product to implement, this commercial product such as such as, data base management system (DBMS), the oracle database of the SQL Server that such as can obtain from the Microsoft of Seattle, Washington, the California Shores inscriptions on bones or tortoise shells (Oracle) company and the MySQL that can obtain from the subsidiary company MySQL AB of the inscriptions on bones or tortoise shells or integration software, such as from the Web Sphere middleware of the IBM of New York Oman gram.But, the computer system running such as SQLServer can be supported according to the aspect of current the disclosed embodiments and the database for various application.

Embodiment described herein is provided for determining novel method by the air-flow of the region in buildings and the equipment in data center and system.In design in an embodiment of the present invention or management tool, the ability of Accurate Prediction air-flow promotes the infrastructure that design is sound, and such as data center and HVAC system, it represents good cooling and performance of ventilating in different layout configurations.In addition, embodiment contributes in the trial and error solution attempting to reach in particular characteristic result the costliness avoided in facility.In addition, the better accuracy in air-flow prediction produces the diagnosis refrigeration of conceptual data center and constructure ventilation system that improve, and can provide more energy-conservation solution.In at least some embodiment as herein described, air-flow is described to generated by air intake opening, gas outlet, fan and leakage.System and method as herein described can use together with the air-source of other types, and this air-source comprises the cooling of other types, ventilation unit and fan.In at least certain embodiments, method is described to determine particular airflow.In at least certain embodiments, describedly determine it is prediction or the estimation of actual airflow.

In at least some embodiment of the present invention discussed herein, the performance of Evaluation and calculation is called as " in real time " once in a while.As referred herein, " in real time " refers to the process completed within about a few second or shorter instead of a few minutes or longer time.Can be there is complicated calculating in this duration, such as relate to typical CFD and calculate.

In above-mentioned at least some embodiment, design and/or the actual parameter of facility (such as data center or buildings) change based on the prediction air-flow in facility.This change can be implemented to improve cooling performance and/or can be implemented with when finding that performance is in predetermined specification, provides cost savings and/or energy-conservation.Such as, the position of equipment rack can be changed and/or housing types or frame configuration type can be changed.In addition, based on the air flow value determined, one or more CRAC or row's formula cooling device can be controlled to regulate air-flow according to the data management system of an embodiment, and in addition, if from refrigeration supply arrangement deficiency of air with provide cool fully time, one or more equipment rack can be controlled to reduce power.

In above-mentioned at least some embodiment, instrument and process are provided for the frame air-flow in determining data center and the air-flow in buildings.In other embodiments, instrument and process can be used to the facility of other types, and can be used to comprise in the Mobile solution at Mobile data center.In addition, can be used according to the process of embodiment as herein described and system in numerous equipment rack with various air intake opening, gas outlet and inner structure.Term gas outlet as used herein and gas outlet can comprise single radial cut region (the such as air vent of frame, close-packed arrays is together effectively to serve as one group of air vent of an opening), maybe can comprise the single radial cut with many individual apertures regions.

Although therefore described several aspects of at least one embodiment of the present invention, should be appreciated that those skilled in the art is easy to expect variously to substitute, amendment and improving.This type of substitutes, revise and improvement is intended for a part of this disclosure, and is intended to fall in the spirit and scope of the present invention.Therefore, description above and accompanying drawing only play example.

Claims (20)

1., for a computer implemented method for modeling air-flow, described method comprises:

Receive the input data relevant to the physical layout of facility;

By computing machine, the expression of described facility is divided into multiple grid cell;

Based on described physical layout, identify that the effect of at least one in jet air-flow, Thermal plume and buoyancy is present in the where in described facility;

If the effect of at least one in described jet air-flow and Thermal plume is present in first group of described multiple grid cell, then the velocity amplitude of first group for described multiple grid cell is specified in operating speed bearing calibration;

Be used for the gas velocity angle value of each grid cell in second group of described multiple grid cell by described computer calculate, described second group is different from described first group;

Amendment is used for the determined gas velocity angle value of any one grid cell in second group of the described multiple grid cell that there is buoyancy wherein; And

Store the air flow value revised on the storage means.

2. method according to claim 1, wherein, calculate described gas velocity angle value also to comprise and use potential barrier method to calculate described gas velocity angle value, and wherein, described method also comprises the equipment configured based at least one in revised air flow value and specified velocity amplitude in described facility.

3. method according to claim 1, also comprises: based on described gas velocity angle value, is used for the temperature value of each grid cell in described multiple grid cell by described computer calculate.

4. method according to claim 1, wherein, described facility is included in the space in data center, and the object in described physical layout comprises at least one equipment rack and at least one refrigeration supply arrangement.

5. method according to claim 1, wherein, described facility is included in the space in buildings, and the object in described physical layout comprises at least one aeration structure and at least one heating equipment.

6. method according to claim 1, wherein, revise determined gas velocity angle value and also comprise:

Calculate the buoyancy velocity amplitude of any one grid cell in second group that is used for described multiple grid cell; And

Add described buoyancy velocity amplitude to determined gas velocity angle value.

7. method according to claim 1, wherein, specify described velocity amplitude also to comprise:

Calculate the jet gas velocity angle value of any one grid cell in first group of the described multiple grid cell in the zone of influence be used in described facility; And

Specify the jet gas velocity angle value of first group of the described multiple grid cell be used in the described zone of influence.

8. method according to claim 1, wherein, specify described velocity amplitude also to comprise:

Calculate the haloing unit speed value of any one grid cell in first group of described multiple grid cell of the object be used in contiguous described facility; And

Specify the haloing unit speed value of first group being used for described multiple grid cell that contiguous described object is arranged.

9. method according to claim 1, also comprises: use process of iteration determines the new air flow value of each grid cell in the grid cell in described facility, and wherein, described new air flow value meets mass balance equation.

10. method according to claim 9, also comprises:

Determine whether the difference between described new air flow value and previous air flow value is greater than threshold value; And

Repeat described process of iteration until described difference is not more than described threshold value.

The system of 11. 1 kinds of modeling air-flows, described system comprises storer and is coupled to the processor of described storer, and described system is configured to:

Receive the input data relevant to the physical layout of facility;

The expression of described facility is divided into multiple grid cell;

Based on described physical layout, identify that the effect of at least one in jet air-flow, Thermal plume and buoyancy is present in the where in described facility;

If the effect of at least one in described jet air-flow and Thermal plume is present in first group of described multiple grid cell, then the velocity amplitude of first group for described multiple grid cell is specified in operating speed bearing calibration;

Be used for the gas velocity angle value of each grid cell in second group of described multiple grid cell by described computer calculate, described second group is different from described first group;

Revise the determined gas velocity angle value for any one grid cell in second group of the described multiple grid cell that there is buoyancy wherein; And

Store the air flow value revised on the storage means.

12. systems according to claim 11, wherein, described system is also configured to use potential barrier method to calculate described gas velocity angle value, and wherein, described system is also configured to the equipment configured based at least one in revised air flow value and specified velocity amplitude in described facility.

13. systems according to claim 11, wherein, described system is also configured to calculate temperature value for each grid cell in described multiple grid cell based on described gas velocity angle value.

14. systems according to claim 11, wherein, described facility is included in the space in data center, and the object in described physical layout comprises at least one equipment rack and at least one refrigeration supply arrangement.

15. systems according to claim 11, wherein, described facility is included in the space in buildings, and the object in described physical layout comprises at least one aeration structure and at least one heating equipment.

16. systems according to claim 11, wherein, described system is also configured to revise determined gas velocity angle value by following manner:

Calculate the buoyancy velocity amplitude of any one grid cell in second group that is used for described multiple grid cell; And

Add described buoyancy velocity amplitude to determined gas velocity angle value.

17. systems according to claim 11, wherein, described system is also configured to specify described velocity amplitude by following manner:

Calculate the jet gas velocity angle value of any one grid cell in first group of the described multiple grid cell be used in the zone of influence; And

Specify the jet gas velocity angle value of first group of the described multiple grid cell be used in the described zone of influence.

18. systems according to claim 11, wherein, described system is also configured to specify described velocity amplitude by following manner:

Calculate the haloing unit speed value of any one grid cell in first group of described multiple grid cell of the object be used in contiguous described facility; And

Specify the described haloing unit speed value of first group being used for described multiple grid cell that contiguous described object is arranged.

19. systems according to claim 11, wherein, described system is also configured to calculate the new air flow value for each grid cell in the grid cell in described facility, and wherein, described new air flow value meets mass balance equation.

20. 1 kinds of non-transitory computer-readable medium, it has the instruction sequence for modeling air-flow stored thereon, and described instruction sequence comprises the instruction that will at least one processor caused to perform following operation:

Receive the input data relevant to the physical layout of facility;

The expression of described facility is divided into multiple grid cell;

Based on described physical layout, identify that the effect of at least one in jet air-flow, Thermal plume and buoyancy is present in the where in described facility;

If the effect of at least one in described jet air-flow and Thermal plume is present in first group of described multiple grid cell, then the velocity amplitude of first group for described multiple grid cell is specified in operating speed bearing calibration;

Be used for the gas velocity angle value of each grid cell in second group of described multiple grid cell by computer calculate, described second group is different from described first group;

Amendment is used for the determined gas velocity angle value of any one grid cell in second group of the described multiple grid cell that there is buoyancy wherein; And

Store the air flow value revised on the storage means.