EP2017755A1 - Verfahren zur unterstützung des entwurfs eines thermischen konvektionsfeldes bzw. substanzdiffusionsfeldes, entwurfsunterstützungssystem und entwurfsunterstützungsprogramm - Google Patents

Verfahren zur unterstützung des entwurfs eines thermischen konvektionsfeldes bzw. substanzdiffusionsfeldes, entwurfsunterstützungssystem und entwurfsunterstützungsprogramm Download PDF

Info

Publication number
EP2017755A1
EP2017755A1 EP06731757A EP06731757A EP2017755A1 EP 2017755 A1 EP2017755 A1 EP 2017755A1 EP 06731757 A EP06731757 A EP 06731757A EP 06731757 A EP06731757 A EP 06731757A EP 2017755 A1 EP2017755 A1 EP 2017755A1
Authority
EP
European Patent Office
Prior art keywords
designing
design
sensitivity
mass
field
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
EP06731757A
Other languages
English (en)
French (fr)
Other versions
EP2017755B1 (de
EP2017755A4 (de
Inventor
Kazunari Momose
Kaoru Ikejima
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
ADVANCED KNOWLEDGE LABORATORY Co Ltd
Osaka University NUC
Original Assignee
Advanced Knowledge Laboratory Coltd
Osaka University NUC
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Advanced Knowledge Laboratory Coltd, Osaka University NUC filed Critical Advanced Knowledge Laboratory Coltd
Publication of EP2017755A1 publication Critical patent/EP2017755A1/de
Publication of EP2017755A4 publication Critical patent/EP2017755A4/de
Application granted granted Critical
Publication of EP2017755B1 publication Critical patent/EP2017755B1/de
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation

Definitions

  • the present invention relates to a design support method, a design support system, and a design support program for supporting the design of a heat convection field or a mass diffusion field.
  • Design of a heat convection field or a mass diffusion field is required in various sites or uses such as, for example, design of an indoor environment using an air conditioning apparatus, thermal design of electronic devices, and management of exhaust gas concentration in plants.
  • a technique of analyzing the influence, of a change in the temperature or the like as a boundary condition, exerted on a minute temperature change or the like of the target position to find a desirable boundary condition has been proposed (see, for example, non-patent documents 1 and 2 specified below).
  • a perturbation equation from the convection field which is used as the reference, is introduced and the adjoint formulation to the perturbation equation is derived.
  • an adjoint problem and the boundary condition therefor are set.
  • changes in various heat transfer characteristics to an arbitrary thermal perturbation and an arbitrary flow perturbation on the boundary are estimated.
  • the present invention made in light of these circumstances has an object of providing a general-purpose and highly convenient design support method and design support system for a heat convection field or a mass diffusion field which significantly reduce the number of times of numerical simulation required to examine the designing parameters for achieving the design purpose.
  • a design support method for a heat convection field or a mass diffusion field comprises a modeling support step of supporting modeling of a designing space upon receipt of a signal from an input device operated by a user while displaying a predetermined image on a display device for displaying a graphic image; a mesh generation step of generating a mesh on the designing space model; a forward analysis step of analyzing the heat convection field or the mass diffusion field of the designing space model by solving an equation of the heat convection field or the mass diffusion field based on an initially set value of a designing parameter input through the input device regarding the meshed designing space model; a purpose setting step of supporting setting of a design purpose upon receipt of a signal from the input device while displaying a predetermined image on the display device; an inverse analysis step of analyzing a sensitivity defined by a change ratio of the design purpose with respect to a change in the designing parameter by solving an adjoint equation in correspondence with the design purpose based on the set design purpose; and a sensitivity display step of displaying information on
  • the inverse analysis step may be a step of analyzing the sensitivity to a predetermined design purpose, regarding temperature, mass concentration, flow rate or pressure, which is defined in a designing space; or is a step of analyzing the sensitivity to a predetermined design purpose, regarding temperature, heat transfer amount, mass concentration, mass transfer amount, flow rate or pressure, which is defined on a boundary of a designing space.
  • a design support system for a heat convection field or a mass diffusion field comprises an input device operable by a user, a computer, and a display device for displaying a graphic image for the user.
  • the computer comprises modeling support means for supporting modeling of a designing space upon receipt of a signal from the input device while displaying a predetermined image on the display device; mesh generation means for generating a mesh on the designing space model; forward analysis means for analyzing the heat convection field or the mass diffusion field of the designing space model by solving an equation of the heat convection field or the mass diffusion field based on an initially set value of a designing parameter input through the input device regarding the meshed designing space model; purpose setting means for supporting setting of a design purpose upon receipt of a signal from the input device while displaying a predetermined image on the display device; inverse analysis means for analyzing a sensitivity defined by a change ratio of the design purpose with respect to a change in the designing parameter by solving an adjoint equation in correspondence with the design purpose based on the set design purpose; and
  • a design support program for a heat convection field or a mass diffusion field is for allowing a computer connected to an input device operable by a user and to a display device for displaying a graphic image for the user to function as modeling support means for supporting modeling of a designing space upon receipt of a signal from the input device while displaying a predetermined image on the display device; mesh generation means for generating a mesh on the designing space model; forward analysis means for analyzing the heat convection field or the mass diffusion field of the designing space model by solving an equation of the heat convection field or the mass diffusion field based on an initially set value of a designing parameter input through the input device regarding the meshed designing space model; purpose setting means for supporting setting of a design purpose upon receipt of a signal from the input device while displaying a predetermined image on the display device; inverse analysis means for analyzing a sensitivity defined by a change ratio of the design purpose with respect to a change in the designing parameter by solving an adjoint equation in correspondence with the design purpose based on the set design
  • the present invention realizes general-purpose and highly convenient design support for a heat convection field or a mass diffusion field which significantly reduces the number of times of numerical simulation required to examine the designing parameters for achieving the design purpose.
  • numerical analysis is performed with predetermined boundary conditions (initially set values of designing parameters).
  • inverse analysis is performed based on adjoint numerical analysis.
  • analysis of sensitivity for a design purpose is performed.
  • the adjoint numerical analysis the sensitivity to any number of designing variables and any distribution designing variable is obtained by performing a numerical analysis once.
  • the result thus obtained is used together with mathematical programming. In this case, an automatic optimization design support system for various devices or various spaces which use heat transfer or mass transfer is realized.
  • governing equations of a heat convection field or a mass diffusion field is set as constraints and is formulated as a nonlinear optimization problem for maximizing or minimizing a given objective function (local temperature, uniformity in a given area, deviation from the target value, etc.).
  • a given objective function local temperature, uniformity in a given area, deviation from the target value, etc.
  • Q is an arbitrary designing space accompanying heat flow and mass diffusion.
  • fT, fC, fu and fp are functions regarding the temperature, mass concentration, flow rate and pressure defined in the space Q, and are determined by the design purpose (for example, maximization, minimization, uniformization or the like of the temperature, mass concentration, flow rate and pressure in a given area).
  • ⁇ T, ⁇ q, ⁇ C, ⁇ m and ⁇ u are respectively boundaries of the space Q in which the temperature, heat flux, concentration, diffusion flux and flow rate vector are controllable.
  • gT, gq, gC, gm and gu are design purposes defined on the respective boundaries (for example, maximization, minimization or the like of the heat transfer amount on the boundary, mass transfer amount on the boundary, or the flow rate on the boundary).
  • L J + ⁇ ⁇ [ p * , u * , T * , C * ] A p ⁇ u ⁇ T ⁇ C ⁇ d ⁇ ⁇ 2
  • A is a vector formed of governing equations of the field having the pressure p, the flow rate vector u, the temperature T and the concentration C as the dependent variables (six nonlinear partial differential equations representing the conservation of mass, momentum, energy and concentration).
  • the boundary condition changes by a minute amount to perform perturbation expansion of the governing equation A in the vicinity of the reference state (for example, in the vicinity of the current designing state).
  • the obtained linear perturbation equation system is applied to the first variation of equation (2) and then divided into domain integrals and boundary integrals with Gauss's divergence theorem or the like.
  • an adjoint operator matrix A* to the linear perturbation equation system is obtained.
  • the adjoint problem the conservation equation of mass, momentum, energy and concentration in an adjoint field having adjoint pressure p*, adjoint flow rate u*, adjoint temperature T* and adjoint concentration C* as dependent variables
  • the boundary conditions thereof are selected as follows.
  • the change (sensitivity) of the objective function to the changes of all the controllable boundary conditions are found.
  • the third term of the right side of equation (5) represents the sensitivity to the change ⁇ K in the overall heat transfer coefficient through the wall, which is not actively controllable but is important in designing. Such sensitivity is usable to evaluate the influence of the structure (heat insulating capability) of the wall, window glass or the like.
  • the target area is set to ⁇ D, and the volume thereof is set to VD, and the average temperature in the area before the optimization is set to Tav.
  • the objective function is converted into the deviation (on the unit of temperature), instead of the variance, such that an error is perceived intuitively.
  • the purpose setting is limited to heat environment for the sake of simplicity.
  • a purpose can be derived in a similar manner for mass diffusion, flow rate and pressure from equations (3) and (4).
  • maximization, uniformization and the target value realization are mentioned.
  • the design purpose may be anything which is a function on the temperature, mass concentration, flow rate or pressure in the designing space or a function on the temperature, thermal flow rate, mass concentration, diffusion flux or flow rate on the boundary of the designing space and is differentiable once for each physical amount. Any such purpose can be appropriately set as required by the designer.
  • a design support system 1 includes an input device 10 including a keyboard, a mouse or the like, a computer 11, and a display device 12 including an LCD, a CRT or the like.
  • the computer 11 includes a central processing unit 13, a memory device 14 including a hard disc drive or the like, an input reception section 15 for receiving a signal from the input device 10, a memory 16, an image output section 17 for outputting a signal to the display device 12, and a CG processing section 18.
  • the specific structure of the design support system 1 is not limited to the above-described structure.
  • a design support program is a computer program for allowing the computer 11 to perform the functions described below.
  • the computer 11 only needs to be communicable with the input device 10 and the display device 12, and may be connected thereto in a wired or wireless manner.
  • the computer 11 may be installed at a different site from the input device 10 or the display device 12. For example, the computer 11 may be connected to a terminal device on the designer side via a network (in a wired or wireless manner).
  • the design is to be made on an indoor space having nine blowoff outlets (upper blowoff outlets 21 through 23, middle blowoff outlets 31 through 33, and lower blowoff outlets 41 through 43) and one opening 20.
  • the temperature at the center of the indoor space is to be maximized.
  • Positional information designated at the time of designing is not limited to information on a predetermined point, and may be information on a predetermined area.
  • the center of the indoor space is designated as the predetermined point.
  • Such a point or area may be designated with a coordinate value (numerical value) using the keyboard or the like.
  • a point or area can be designated on a screen by clicking and dragging the mouse for improved convenience. The point or area designated in this manner is graphically displayed in a designing space model as a target 25.
  • a designing space model is created.
  • the system 1 has a model creation function, and the user (designer) uses the input device 10 (for example, inputs a command to the keyboard or performs an operation on the mouse such as clicking, dragging or the like) to select a menu bar, draw, input a command or the like.
  • an arbitrary designing space can be set easily.
  • an indoor space model is set as shown in FIG. 3 .
  • nine blowoff outlets (upper blowoff outlets 21 through 23, middle blowoff outlets 31 through 33, and lower blowoff outlets 41 through 43) are formed in one of two walls facing each other, one opening 20 is formed in an upper part of the other wall, and the ceiling and the floor are cooling faces.
  • initial values (boundary conditions) of designing parameters are also set when the model is created.
  • the designing parameters are, for example, the temperature of the air blown through each blowoff outlet, the blowing direction, the blowing speed, the temperature of the ceiling and the floor, etc.
  • the designing parameters may be set by a different step from the creation of the model (step S1).
  • the "mesh (M)" button in the menu bar is clicked to automatically generate meshes on the model (step S2).
  • the number of meshes may be set to a predetermined value or appropriately set by the user. Namely, the number of meshes may be increased in order to raise the precision of analysis, or may be decreased in order to raise the speed of analysis.
  • the design support system 1 may be set to allow the user to select the number of meshes freely.
  • FIG. 4 shows an exemplary image shown on the display device during the forward analysis.
  • step S4 When the forward analysis is finished and the user clicks the "view (V)" button in the menu bar (see FIG. 6(a) ), the speed distribution and temperature distribution on an arbitrary face of the model are displayed as shown in FIG. 5 (step S4).
  • the face to be displayed may be changed easily using the mouse or the like.
  • the image of speed distribution (see FIG. 5(a) ) and the image of temperature distribution (see FIG. 5(b) ) may be easily switched to each other.
  • the design purpose is first set (step S5).
  • the design purpose is to maximize the temperature at the center of the indoor space. Therefore, the above-mentioned target 25 is created at the center of the indoor space (see FIG. 6(a) ), and the object scalar is set to the temperature (see FIG. 6(b) ).
  • the target 25 is set small enough to look substantially like a point (see FIG. 6(c) ). In the case where the target 25 is a predetermined area, the target 25 may be set large in the image of FIG. 6(c) . In this example, there is one target 25.
  • the design purpose is regarding, for example, the temperature at a plurality of points in the indoor space, a plurality of targets 25 may be set.
  • the computer is set to the inverse analysis mode "maximize" and other conditions are also set (in this example, convergence conditions and maximum possible number of time of repetition). Then, the execution/continue button is clicked. Then, as shown in FIG. 7(b) , the inverse analysis is executed (step S6).
  • the inverse analysis is finished and the user clicks the "view (V)" button in the menu bar.
  • the sensitivity to each parameter (a ratio of change in the design purpose with respect to a change in each design parameter; in this example, in the case where the temperature at the center of the indoor space is to be raised, how much influence is exerted on the temperature at the center of the indoor space by the change in the blowoff temperature or the like) is displayed quantitatively and visually.
  • the user may use the keyboard, mouse or the like to arbitrarily switch the images displaying the sensitivity to the parameters.
  • FIG. 8(a), FIG. 8(b) , FIG. 9(a), FIG. 9(b) and FIG. 10(a) show the level of sensitivity with the level of color darkness as an example.
  • the manner of displaying the sensitivity is not specifically limited, and the level of sensitivity may be displayed with contour lines, colors (e.g., red for a highly sensitive area and blue for a non-sensitive area) or gradation levels.
  • the level of sensitivity may be displayed with a graph.
  • the system 1 performs the above-described design support operation. Then, the user examines the designing particulars while checking the sensitivity displayed on the screen.
  • FIG. 8(a) it is understood that for raising the temperature at the center of the indoor space, it is most efficient to raise the blowoff temperature at the lower left and lower right blowoff outlets 41 and 43 (hereinafter, regarding the reference numerals, see FIG. 3 ) and to lower the blowoff temperature at the upper middle and central outlets 22 and 32.
  • FIG. 8(b) it is understood that for raising the temperature at the center of the indoor space, it is preferable to increase the flow rate in the blowoff direction at the lower middle blowoff outlet 42.
  • FIG. 9(b) it is understood that it is preferable to increase the flow rate in the vertical direction at the lower middle blowoff outlet 42.
  • FIG. 10(a) it is understood that it is preferable to increase the heating amount at the central position in the top-bottom direction of the wall. Accordingly, it is conceivable to, for example, install a heater at this position.
  • the sensitivity of the area is displayed on the screen of the display device 12.
  • the user for example, drags the mouse to designate a predetermined area R1 on the screen shown in FIG. 10(b)
  • a temperature change of the target with respect to a change in the heating amount in the area R1 is displayed on the screen.
  • the heating amount in the area R1 is increased by 1 kW
  • the temperature of the target is increased by 0.5124°C. Based on this, when a heater is to be installed on the wall, a rough estimate on the heating amount required of the heater can be obtained by designating an area corresponding to the area size of the heater at the position where the heater is to be installed.
  • FIG. 11 (a) and FIG. 11(b) the design is to be made on a rectangular parallelepiped indoor space.
  • Three blowoff outlets 51, 52 and 53 are provided on one of two walls facing each other in a longitudinal direction.
  • a window 54 is provided on one of two walls perpendicular to the wall having the blowoff outlets.
  • An opening 55 is provided on the wall facing the wall having the window 54.
  • a partition 60 is provided at the center of the room.
  • the partition 60 does not reach the ceiling with a gap being made between the partition 60 and the ceiling.
  • the partition 60 reaches the ceiling.
  • FIG. 12 through FIG. 17 show the results of numerical analysis performed on each of the above exemplary models in substantially the same manner as described above.
  • FIG. 12 (a) and FIG. 12(b) show the sensitivity to the temperature.
  • FIG. 13 (a) and FIG. 13 (b) show the sensitivity to the heating amount on the wall.
  • FIG. 14(a) and FIG. 14(b) show the sensitivity to the air flow in the blowoff direction (X direction) on the wall.
  • FIG. 15 (a) and FIG. 15 (b) show the sensitivity to the air flow in the vertical direction (Z direction).
  • FIG. 16(a) and FIG. 16(b) show the sensitivity to the overall heat transfer coefficient through the window.
  • FIG. 17(a) and FIG. 17(b) show numerical data on the sensitivity to each of the designing parameters.
  • the sensitivity to each parameter varies depending on the height of the partition 60.
  • the partition 60 (hereinafter, regarding the reference numerals, see FIG. 11(a) and FIG. 11(b) ) does not reach the ceiling as shown in FIG. 12(a)
  • the sensitivity to the temperature at the left blowoff outlet 51 is high.
  • the partition 60 reaches the ceiling as shown in FIG. 12(b) the sensitivity to the temperature at the left blowoff outlet 51 is low, and the sensitivity to the temperature at the right blowoff outlet 53 is highest.
  • the design support method and the design support system 1 can significantly reduce the number of times of numerical simulation required to examine the designing parameters for achieving the design purpose, thus can significantly reduce the load on the computer, and also can significantly reduce the time for the calculation, as compared with the usual design method which is a combination of the numerical simulation (forward analysis) and the optimization algorithm.
  • the design purpose is first set and then inverse analysis is performed to find the sensitivity to each of various parameters. Then, the information on the sensitivity is graphically displayed on the display device 12. Therefore, when designing a heat convection field, a guide regarding which parameter(s), and how much of the parameter(s), should be changed can be obtained easily. For example, in the case where the temperature at a predetermined point in the indoor space is to be put to a predetermined level, a specific design guide can be obtained regarding at which blowoff outlet(s), and how much, the temperature should be raised. Thus, a clear design guide can be obtained without relying on the intuition of the designer, unlike in the conventional art.
  • the information on the sensitivity is displayed on a three-dimensional model visually and quantitatively. Therefore, the guide for design can be obtained easily and quickly. Specifically, the information on the sensitivity is displayed using color gradation levels, contour lines, different colors or the like. Therefore, the quantitative information on the sensitivity in a three-dimensional space can be grasped at a glance.
  • the design guide can be intuitively perceived. This can realize a highly convenient design support.
  • a designing space model can be created while referring to a three-dimensional graphic image displayed on the display device 12. This allows the designing space model to be created simply and quickly at the site of design.
  • the ability of automatically generating meshes on the model provides a high level of convenience.
  • the creation of the model, generation of meshes, forward analysis, setting of a design purpose, inverse analysis and display of sensitivity as a result of the inverse analysis can be executed with only a series of simple operations using the input device 10.
  • the design support system 1 in this embodiment has a modeler function, a forward analysis function, a purpose setting function, an inverse analysis function and a sensitivity display function.
  • the design support system 1 exhibits superb convenience at the site of design.
  • the design purpose is to maximize the temperature at a predetermined point in an indoor space.
  • the design purpose may be something else; for example, to minimize the temperature at a predetermined point in an indoor space.
  • the design purpose may be to maximize or minimize the flow rate, the mass concentration, or the pressure at a predetermined point in an indoor space.
  • the design purpose may be to maximize or minimize the temperature, mass concentration, flow rate or pressure in a predetermined area of a designing space.
  • the design purpose may be to maximize or minimize the temperature, heat transfer amount, mass concentration, mass transfer amount, flow rate or pressure at a predetermined point or in a predetermined area on a boundary of a designing space (for example, a wall, a window, etc. which divides the indoor space, window, etc.).
  • the design purpose may be to uniformize the temperature in a predetermined area in an indoor space. For example, it may be occasionally desired to uniformize the temperature of an area where there is always somebody present (for example, around the desk in a work office) from the viewpoint of the indoor air conditioning environment. In such a case, the indoor environment is designed to uniformize the temperature in a predetermined area, not to maximize or minimize the temperature at a predetermined point.
  • the design support system 1 in this embodiment can support such a design of heat convection field by setting a target area and the temperature so as to uniformize the temperature in a predetermined area.
  • the design purpose may be to uniformize the temperature, mass concentration, flow rate or pressure in a predetermined area in a designing space.
  • the design purpose may be to uniformize the temperature, heat transfer amount, mass concentration, mass transfer amount, flow rate or pressure in a predetermined area on a boundary of a designing space.
  • the design purpose may be to put the temperature at a plurality of points in an indoor space to a respective predetermined level. For example, from the viewpoint of the air conditioning environment of a work office in winter, it may be occasionally desired to make an area where there is always somebody present warm, but to make an area where there are office electronics relatively cool.
  • the design support system 1 can set a plurality of targets and set an appropriate temperature to be achieved for each of the targets.
  • the design support system 1 can also support a design of heat convection field in such a case.
  • the design purpose may be to put the temperature, mass concentration, flow rate or pressure at one or each of a plurality of points or in one or each of a plurality of areas in a designing space to a predetermined level.
  • the design purpose may be to put the temperature, heat transfer amount, mass concentration, mass transfer amount, flow rate or pressure at one or each of a plurality of points or in one or each of a plurality of areas on a boundary of a designing space to a predetermined level.
  • the temperature at a plurality of points in an indoor space may be actually measured in advance, and the design purpose may be to put the temperature of each of the points to the actually measured value.
  • the design support system 1 may be used to identify an unknown designing parameter.
  • the amount of heat leaking from the wall can be analyzed, and thus a design of heat insulation conditions of the wall (design of specifications of a heat insulator) or the like can be supported.
  • the design purpose may be to put the temperature, mass concentration, flow rate or pressure at one or each of a plurality of points or in one or each of a plurality of areas in a designing space to the actually measured value of the temperature, mass concentration, flow rate or pressure at the respective point or area.
  • the design purpose may be to put the temperature, heat transfer amount, mass concentration, mass transfer amount, flow rate or pressure at one or each of a plurality of points or one or each of a plurality of areas on a boundary of a designing space to the actually measured value of the temperature, heat transfer amount, mass concentration, mass transfer amount, flow rate or pressure at the respective point or area.
  • the design support method described above is completed when the sensitivity is displayed after inverse analysis (see FIG. 2 ).
  • optimization may be performed by, after the inverse analysis, automatically changing the designing variable based on mathematical programming (gradient-type nonlinear optimization) and repeating forward analysis and inverse analysis.
  • the design support system 1 may include an optimization function of automatically optimizing a designing parameter for the design purpose.
  • step S11 the designing space is modeled in step S11, and then meshes are automatically generated in step S12. Then, in step S13, the design purpose and a convergence condition are set. After performing forward analysis in step S14, inverse analysis is performed in step S15. Then, in step S16, it is determined whether the convergence condition has been fulfilled or not. When the convergence condition has not been fulfilled, the value of the designing parameter is changed in step S17, and the forward analysis in step S14 and the inverse analysis in step S15 are repeated. The forward analysis and the inverse analysis are repeated until the convergence condition is fulfilled in step S16. For such an optimization calculation, it is preferable that the user sets the upper limit and the lower limit so that the repeated calculation does not diverge.
  • the sensitivity found by the inverse analysis is displayed on the screen as it is.
  • a change amount of the designing parameter which is necessary for the required change of the design purpose may be displayed.
  • a change amount of the temperature at the blowoff outlet for example, the amount of temperature which needs to be changed in order to raise the temperature of the target by 1°C
  • the designing parameter is the temperature on a boundary of an indoor space or the like.
  • the designing parameter is not limited to a so-called boundary condition but encompasses the internal heat generation, internal mass diffusion, the shape of the physical object and the like in a designing space or the like.
  • the above-described embodiments are applied to air conditioning design in an indoor space.
  • the present invention is not limited to air conditioning design, but is applicable to design of any heat convection field or mass diffusion field.
  • the present invention is not limited to indoor or outdoor air conditioning design, but is applicable to, for example, internal thermal design of the housing of electronics such as computers and the like, internal environmental design of plants and factories, estimate of the influence of exhaust gas on parking areas, estimate of behavior of the thermal fluid in a heat storage tank, and the like.
  • the present invention is useful for, for example, designing various devices or various spaces (environment) using heat or mass transfer.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Air Conditioning Control Device (AREA)
EP06731757.8A 2006-04-13 2006-04-13 System zur Entwurfsunterstützung eines Wärmekonvektionsfeldes Active EP2017755B1 (de)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2006/307821 WO2007122677A1 (ja) 2006-04-13 2006-04-13 熱対流場・物質拡散場の設計支援方法、設計支援システム、および設計支援用プログラム

Publications (3)

Publication Number Publication Date
EP2017755A1 true EP2017755A1 (de) 2009-01-21
EP2017755A4 EP2017755A4 (de) 2013-01-02
EP2017755B1 EP2017755B1 (de) 2017-06-07

Family

ID=38624609

Family Applications (1)

Application Number Title Priority Date Filing Date
EP06731757.8A Active EP2017755B1 (de) 2006-04-13 2006-04-13 System zur Entwurfsunterstützung eines Wärmekonvektionsfeldes

Country Status (4)

Country Link
US (1) US8073662B2 (de)
EP (1) EP2017755B1 (de)
JP (1) JP4016066B1 (de)
WO (1) WO2007122677A1 (de)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102549347A (zh) * 2009-10-21 2012-07-04 阿自倍尔株式会社 空调控制装置以及方法
CN104374043A (zh) * 2013-08-12 2015-02-25 阿自倍尔株式会社 空调控制装置以及方法
US20150371150A1 (en) * 2014-06-20 2015-12-24 International Business Machines Corporation Analysis device, analysis method, and program
EP4279830A4 (de) * 2021-01-14 2024-02-14 Mitsubishi Electric Corporation Klimaanlagensteuerungsvorrichtung

Families Citing this family (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR100830095B1 (ko) * 2007-11-12 2008-05-20 충남대학교산학협력단 냉방부하 예측방법
JP5018487B2 (ja) * 2008-01-14 2012-09-05 富士通株式会社 製造バラつきを考慮した多目的最適化設計支援装置、方法、及びプログラム
US8229776B1 (en) * 2008-05-06 2012-07-24 The United States Of America As Represented By The Secretary Of The Navy Evaluation of subsystem technology in a system-of-subsystems environment
US8229713B2 (en) * 2009-08-12 2012-07-24 International Business Machines Corporation Methods and techniques for creating and visualizing thermal zones
JP5240132B2 (ja) * 2009-09-04 2013-07-17 富士通株式会社 熱流体シミュレーション解析装置
JP5525960B2 (ja) * 2010-08-10 2014-06-18 アズビル株式会社 空調制御装置および方法
KR101067020B1 (ko) * 2011-01-14 2011-09-22 서울대학교산학협력단 최적화 기법을 이용한 자연스러운 이미지의 자동 생성 방법
JP5973829B2 (ja) * 2012-08-01 2016-08-23 アズビル株式会社 空調制御装置および方法
JP5973830B2 (ja) * 2012-08-01 2016-08-23 アズビル株式会社 空調制御解表示装置および方法
JP6253053B2 (ja) * 2013-11-29 2017-12-27 富士通株式会社 データ探索装置、データ探索装置の制御方法およびデータ探索装置の制御プログラム
JP2018021693A (ja) * 2016-08-02 2018-02-08 Ｊｆｅスチール株式会社 温度分布推定装置および温度分布推定方法
JP2018128840A (ja) * 2017-02-08 2018-08-16 富士通株式会社 発熱密度算出コンピュータプログラム、発熱密度算出方法及び情報処理装置
JP6655035B2 (ja) * 2017-03-10 2020-02-26 株式会社アドバンスドナレッジ研究所 情報処理システム、プログラム、環境管理システム、及び、設備管理システム
CN111241710B (zh) * 2020-02-19 2023-10-17 北京和隆优化科技股份有限公司 一种基于蜂窝陶瓷蓄热体换热系统的建模方法
CN113343598B (zh) * 2021-06-11 2023-07-18 西安交通大学 一种基于解耦模式的自然对流换热场景快速仿真系统
WO2024057483A1 (ja) * 2022-09-15 2024-03-21 三菱電機株式会社 空調制御装置および空調システム
CN116935985B (zh) * 2023-07-17 2024-03-15 中国地质调查局油气资源调查中心 一种煤气化过程中实验参数变化的敏感性分析方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4964060A (en) * 1985-12-04 1990-10-16 Hartsog Charles H Computer aided building plan review system and process
JPH1153422A (ja) * 1997-08-08 1999-02-26 Hitachi Ltd 最適設計支援装置および方法ならびに電磁測定装置
US7894943B2 (en) * 2005-06-30 2011-02-22 Sloup Charles J Real-time global optimization of building setpoints and sequence of operation

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CHENG-HUNG HUANG, HUNG-CHI LO: "A three-dimensional inverse problem in estimating the internal heat flux of housing for high speed motors", APPLIED THERMAL ENGINEERING, PERGAMON, OXFORD, GB, vol. 26, no. 14-15, 7 February 2006 (2006-02-07), pages 1515-1529, XP024987777, ISSN: 1359-4311, DOI: 10.1016/J.APPLTHERMALENG.2005.12.009 [retrieved on 2006-10-01] *
SCHMIDT R R ET AL: "Challanges of data center thermal management", IBM JOURNAL OF RESEARCH AND DEVELOPMENT, INTERNATIONAL BUSINESS MACHINES CORPORATION, NEW YORK, NY, US, vol. 49, 1 September 2005 (2005-09-01), pages 709-723, XP002458425, ISSN: 0018-8646 *
See also references of WO2007122677A1 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102549347A (zh) * 2009-10-21 2012-07-04 阿自倍尔株式会社 空调控制装置以及方法
CN102549347B (zh) * 2009-10-21 2015-01-21 阿自倍尔株式会社 空调控制装置以及方法
US9200813B2 (en) 2009-10-21 2015-12-01 Azbil Corporation Air-conditioning controlling device and method
CN104374043A (zh) * 2013-08-12 2015-02-25 阿自倍尔株式会社 空调控制装置以及方法
US20150371150A1 (en) * 2014-06-20 2015-12-24 International Business Machines Corporation Analysis device, analysis method, and program
US9626631B2 (en) * 2014-06-20 2017-04-18 International Business Machines Corporation Analysis device, analysis method, and program
EP4279830A4 (de) * 2021-01-14 2024-02-14 Mitsubishi Electric Corporation Klimaanlagensteuerungsvorrichtung

Also Published As

Publication number Publication date
EP2017755B1 (de) 2017-06-07
US8073662B2 (en) 2011-12-06
JP4016066B1 (ja) 2007-12-05
EP2017755A4 (de) 2013-01-02
US20090276193A1 (en) 2009-11-05
JPWO2007122677A1 (ja) 2009-08-27
WO2007122677A1 (ja) 2007-11-01

Similar Documents

Publication Publication Date Title
US8073662B2 (en) Design support method, design support system, and design support program for heat convection field
Schlueter et al. Linking BIM and Design of Experiments to balance architectural and technical design factors for energy performance
Wang et al. Advances in building simulation and computational techniques: A review between 1987 and 2014
Dols et al. CONTAMW 2.0 User Manual:.
Wurtz et al. Two-and three-dimensional natural and mixed convection simulation using modular zonal models in buildings
Ahn et al. BIM interface for full vs. semi-automated building energy simulation
Al Ka’bi Comparison of energy simulation applications used in green building
Guyot et al. Building energy model calibration: A detailed case study using sub-hourly measured data
Naboni et al. Thermal comfort-CFD maps for architectural interior design
EP3267337A1 (de) Verfahren zur identifizierung technischer konstruktionslösungen
CN109073753B (zh) 用于生成能量模型并跟踪能量模型演进的系统和方法
Goethals et al. Sensitivity of night cooling performance to room/system design: surrogate models based on CFD
JP7459935B2 (ja) 推定方法、シミュレーション方法、推定装置、及び推定プログラム
Shahsavari et al. Uncertainty and Sensitivity Analysis Using Building Information Modeling (An Energy Analysis Test Case)
Wurtz et al. An equation-based simulation environment to investigate fast building simulation
Jung et al. Appropriate level of development of in-situ building information modeling for existing building energy modeling implementation
Musy et al. Automatically generated zonal models for building air flow simulation: principles and applications
Lan et al. Optimal zoning for building zonal model of large-scale indoor space
Desta et al. Combining CFD and data-based mechanistic (DBM) modelling approaches
Zong et al. A holistic two-stage decision-making methodology for passive and active building design strategies under uncertainty
CN115344919A (zh) 建筑内部环境分析方法、装置、设备及介质
Kircher et al. A computationally efficient, high-fidelity testbed for building climate control
Sanguinetti et al. Automated energy performance visualization for BIM
Rincón-Casado et al. 3D internal forced convection heat-transfer correlations from CFD for building performance simulation
Schoch et al. Building-volume designs with optimal life-cycle costs

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20081110

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LI LT LU LV MC NL PL PT RO SE SI SK TR

AX Request for extension of the european patent

Extension state: AL BA HR MK YU

DAX Request for extension of the european patent (deleted)
A4 Supplementary search report drawn up and despatched

Effective date: 20121129

RIC1 Information provided on ipc code assigned before grant

Ipc: G06F 17/50 20060101AFI20121123BHEP

17Q First examination report despatched

Effective date: 20150528

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: EXAMINATION IS IN PROGRESS

GRAJ Information related to disapproval of communication of intention to grant by the applicant or resumption of examination proceedings by the epo deleted

Free format text: ORIGINAL CODE: EPIDOSDIGR1

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: GRANT OF PATENT IS INTENDED

RIN1 Information on inventor provided before grant (corrected)

Inventor name: IKEJIMA, KAORU

Inventor name: MOMOSE, KAZUNARI

INTG Intention to grant announced

Effective date: 20161130

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: THE PATENT HAS BEEN GRANTED

RAP1 Party data changed (applicant data changed or rights of an application transferred)

Owner name: ADVANCED KNOWLEDGE LABORATORY CO., LTD.

Owner name: OSAKA UNIVERSITY

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IS IT LI LT LU LV MC NL PL PT RO SE SI SK TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

Ref country code: AT

Ref legal event code: REF

Ref document number: 899708

Country of ref document: AT

Kind code of ref document: T

Effective date: 20170615

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

REG Reference to a national code

Ref country code: DE

Ref legal event code: R096

Ref document number: 602006052736

Country of ref document: DE

REG Reference to a national code

Ref country code: NL

Ref legal event code: MP

Effective date: 20170607

REG Reference to a national code

Ref country code: LT

Ref legal event code: MG4D

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170908

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

Ref country code: ES

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

Ref country code: LT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

REG Reference to a national code

Ref country code: AT

Ref legal event code: MK05

Ref document number: 899708

Country of ref document: AT

Kind code of ref document: T

Effective date: 20170607

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: NL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

Ref country code: SE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

Ref country code: LV

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170907

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

Ref country code: AT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: PL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

Ref country code: IT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

Ref country code: IS

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20171007

REG Reference to a national code

Ref country code: DE

Ref legal event code: R097

Ref document number: 602006052736

Country of ref document: DE

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 13

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

26N No opposition filed

Effective date: 20180308

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: MC

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

REG Reference to a national code

Ref country code: BE

Ref legal event code: MM

Effective date: 20180430

REG Reference to a national code

Ref country code: IE

Ref legal event code: MM4A

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20180413

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20180430

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20180430

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20180430

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20180413

REG Reference to a national code

Ref country code: DE

Ref legal event code: R079

Ref document number: 602006052736

Country of ref document: DE

Free format text: PREVIOUS MAIN CLASS: G06F0017500000

Ipc: G06F0030000000

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT; INVALID AB INITIO

Effective date: 20060413

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20170607

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20240423

Year of fee payment: 19

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: DE

Payment date: 20240425

Year of fee payment: 19

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: FR

Payment date: 20240423

Year of fee payment: 19