WO2022153431A1 - 空気調和機制御装置 - Google Patents

空気調和機制御装置 Download PDF

Info

Publication number
WO2022153431A1
WO2022153431A1 PCT/JP2021/001012 JP2021001012W WO2022153431A1 WO 2022153431 A1 WO2022153431 A1 WO 2022153431A1 JP 2021001012 W JP2021001012 W JP 2021001012W WO 2022153431 A1 WO2022153431 A1 WO 2022153431A1
Authority
WO
WIPO (PCT)
Prior art keywords
unit
air conditioner
air
objective function
coefficient
Prior art date
Application number
PCT/JP2021/001012
Other languages
English (en)
French (fr)
Inventor
隆太 田中
美緒 元谷
Original Assignee
三菱電機株式会社
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 三菱電機株式会社 filed Critical 三菱電機株式会社
Priority to US18/267,107 priority Critical patent/US20240044540A1/en
Priority to JP2021532450A priority patent/JP7038916B1/ja
Priority to EP21919325.7A priority patent/EP4279830A4/en
Priority to PCT/JP2021/001012 priority patent/WO2022153431A1/ja
Publication of WO2022153431A1 publication Critical patent/WO2022153431A1/ja

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24—HEATING; RANGES; VENTILATING
    • F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F11/00—Control or safety arrangements
    • F24F11/62—Control or safety arrangements characterised by the type of control or by internal processing, e.g. using fuzzy logic, adaptive control or estimation of values
    • F24F11/63—Electronic processing
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24—HEATING; RANGES; VENTILATING
    • F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F11/00—Control or safety arrangements
    • F24F11/62—Control or safety arrangements characterised by the type of control or by internal processing, e.g. using fuzzy logic, adaptive control or estimation of values
    • F24F11/63—Electronic processing
    • F24F11/64—Electronic processing using pre-stored data
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D23/00—Control of temperature
    • G05D23/19—Control of temperature characterised by the use of electric means
    • G05D23/1927—Control of temperature characterised by the use of electric means using a plurality of sensors
    • G05D23/193—Control of temperature characterised by the use of electric means using a plurality of sensors sensing the temperaure in different places in thermal relationship with one or more spaces
    • G05D23/1932—Control of temperature characterised by the use of electric means using a plurality of sensors sensing the temperaure in different places in thermal relationship with one or more spaces to control the temperature of a plurality of spaces
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24—HEATING; RANGES; VENTILATING
    • F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F2110/00—Control inputs relating to air properties
    • F24F2110/10—Temperature
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F24—HEATING; RANGES; VENTILATING
    • F24F—AIR-CONDITIONING; AIR-HUMIDIFICATION; VENTILATION; USE OF AIR CURRENTS FOR SCREENING
    • F24F2120/00—Control inputs relating to users or occupants
    • F24F2120/10—Occupancy
    • F24F2120/12—Position of occupants

Definitions

  • This technology is related to the air conditioner control device.
  • it relates to the control of the thermal environment of the space subject to air conditioning in which the air conditioner performs air conditioning.
  • thermo-fluid analysis model of the space subject to air conditioning (hereinafter referred to as the space subject to air conditioning) to perform simulation of the space subject to air conditioning and inverse analysis of the space subject to air conditioning based on the simulation results.
  • the desired place can be controlled to the desired thermal environment.
  • an air conditioner is generally equipped with one or more temperature sensors that detect the temperature.
  • the air conditioner control device that controls the air conditioner controls the air conditioner so that the temperature related to the detection of the temperature sensor approaches the control target value such as the set temperature by a predetermined method, and the air conditioner controls the air conditioner. It controls the thermal environment.
  • the temperature related to the detection of the temperature sensor is the temperature of the local space around the position where the temperature sensor is installed, and does not necessarily reflect the temperature of the air harmonization target space. Therefore, the air conditioner control device grasps the boundary condition temperature of the air conditioning target space and grasps the thermal environment itself of the air conditioning target space by simulation. Then, a method has been proposed in which an air conditioner control device performs inverse analysis using simulation results, obtains a required outlet temperature, and controls the air conditioner to a desired thermal environment (see, for example, Patent Document 1). ).
  • control method by the air conditioner control device described above does not consider the conditions including the equipment side of the air conditioner. Therefore, in reality, there is a problem that the control of the air conditioner may be inefficient.
  • the purpose is to obtain an air conditioner control device capable of realizing more efficient control of the air conditioner.
  • the air conditioner control device is an air conditioner control device that controls an air conditioner that performs air conditioning in the air conditioning target space, and is a model on the airflow side in a simulation regarding a thermal environment related to the air conditioning target space.
  • the value of the objective function in the air-conditioning target space is calculated by the air-conditioning side model building unit to be constructed, the equipment-side model building unit that builds the equipment-side model that simulates the behavioral constraints and equipment capabilities of the air-conditioning machine. It has an objective function calculation unit and a sensitivity derivation unit that derives the sensitivity that is the fluctuation of the objective function when the control variable is changed, and inversely analyzes the objective function based on the air-conditioning side model and the objective function based on the device side model. It includes a coupled optimization execution unit that optimizes using a method and calculates the optimum solution, and a control target value determination unit that determines the control target value of the air conditioner from the optimum solution.
  • the airflow side model construction unit and the equipment side model construction unit are provided, and the coupled optimization execution unit calculates the optimum solution considering the constraints not only in the airflow side objective function but also in the equipment side objective function. do. Therefore, it is possible to obtain comfortable and energy-saving operating conditions for the air conditioner.
  • FIG. 5 is a diagram illustrating an example of processing performed by the coupled optimization execution unit 106 and the control target value determination unit 109 according to the first embodiment.
  • FIG. 2 It is a figure which shows the configuration example of the air conditioner control device 1 which concerns on Embodiment 2.
  • FIG. It is a schematic diagram explaining the acquisition of the data about the air-conditioning target area which concerns on Embodiment 2.
  • FIG. It is a figure explaining the process performed by the air conditioner control device 1 which concerns on Embodiment 2.
  • FIG. It is a figure which shows the configuration example of the air conditioner control device 1 which concerns on Embodiment 3.
  • FIG. It is a schematic diagram explaining an example of acquisition of data about the air-conditioning target area which concerns on Embodiment 3.
  • FIG. It is a figure explaining an example of the process concerning the wearable terminal apparatus H204 in the air conditioner control apparatus 1 which concerns on Embodiment 3.
  • FIG. It is a schematic diagram explaining another example of the acquisition of the data about the air-conditioning target area which concerns on Embodiment 3.
  • FIG. It is a figure explaining an example of the process concerning the mobile terminal apparatus H304 in the air conditioner control apparatus 1 which concerns on Embodiment 3.
  • FIG. It is a figure which shows the configuration example of the air conditioner control device 1 which concerns on Embodiment 4.
  • FIG. It is a figure explaining an example of the process which concerns on the determination of the coefficient ⁇ which concerns on Embodiment 4.
  • FIG. It is a figure explaining an example of the process which concerns on the determination of the coefficient ⁇ and the coefficient ⁇ which concerns on Embodiment 4.
  • FIG. It is a figure explaining an example of the process performed by the coefficient determination part 116 which concerns on Embodiment 4.
  • FIG. It is a figure which shows the configuration example of the air conditioner control device 1 which concerns on Embodiment 5.
  • FIG. 1 is a diagram showing an example of the configuration of an air conditioning system including the air conditioner control device 1 according to the first embodiment.
  • the air conditioner control device 1 is a device that controls the operation of the air conditioner 2.
  • the air conditioner control device 1 is communicably connected to the air conditioner 2 and the sensor unit 4 via the control network 5.
  • the air conditioner 2 includes an outdoor unit 21, an indoor unit 22, and a remote controller 23 as components.
  • the outdoor unit 21 cools or heats a thermal refrigerant such as a refrigerant and water.
  • the indoor unit 22 exchanges heat between the air in the air conditioning target space such as a room and a thermal refrigerant, heats or cools the air in the air conditioning target space, and adjusts the temperature in the air conditioning target space.
  • the remote controller 23 is, for example, a device used by a resident to switch ON / OFF of the indoor unit 22 and manually change settings such as a set temperature, an air volume, and an air direction.
  • the sensor unit 4 is a group of sensors for detecting and measuring a physical quantity. Here, various physical quantities related to environmental conditions inside and outside the air-conditioning target space such as temperature and pressure are detected.
  • a sensor 41 and a sensor 42 are provided as sensor units. Here, it is assumed that a plurality of temperature sensors are mainly installed in the space subject to air conditioning.
  • the control network 5 is a telecommunication line connecting the air conditioner control device 1, the air conditioner 2, and the sensor unit 4.
  • FIG. 2 is a diagram showing a configuration example of the air conditioner control device 1 according to the first embodiment.
  • FIG. 2 also shows data input / output to / from the air conditioner control device 1.
  • the air conditioner control device 1 includes a control processing device 100 and a data storage device 120.
  • control processing device 100 is composed of, for example, a device that performs control calculation processing such as a computer centered on a CPU (Central Processing Unit). Then, the control processing device 100 realizes the processing by executing a pre-programmed processing procedure performed by each part described later.
  • the configuration is not limited to this, and each part of the control processing device 100 may be configured by a dedicated device (hardware).
  • the data storage device 120 is a device that stores data required when the control processing device 100 performs processing.
  • the data storage device 120 in the first embodiment stores various data acquired by the data acquisition unit 101, which will be described later. Further, as described above, the data storage device 120 of the first embodiment has data in which the processing procedure performed by each part of the control processing device 100 is programmed.
  • the data storage device 120 includes a volatile storage device (not shown) such as a random access memory (RAM) capable of temporarily storing data, a solid state disk, and a non-volatile auxiliary such as a flash memory capable of storing data for a long period of time. It has a storage device (not shown).
  • the control processing device 100 of the first embodiment includes a data acquisition unit 101, an airflow side model construction unit 102, a target area designation unit 103, a temperature sensor extraction unit 104, a device side model construction unit 105, and a coupled optimization execution unit 106. It has a control target value determination unit 109 and an air conditioning control command unit 110.
  • the data acquisition unit 101 receives signals from the sensor unit 4 and an external device, and acquires various data included in the signals.
  • the data acquisition unit 101 acquires the room shape data D1, the device information data D2, the air temperature data D3, and the area information data D4.
  • the data acquisition unit 101 stores the acquired data in the data storage device 120.
  • the room shape data D1 is data related to the shape and dimensions of the air-conditioning target space such as the width, depth, and height of the room inside the room.
  • the device information data D2 is data related to the device of the air conditioner 2, including, for example, the configuration of the refrigerant circuit in the air conditioner 2, the rated capacity of the device, the specifications of the device of the drive system, and the like.
  • the air temperature data D3 is data relating to the temperatures of various types of air in the air harmonization target space sent from the sensor unit 4.
  • the area information data D4 is data of information regarding the area of the living area in the air-conditioned space.
  • the airflow side model building unit 102 performs a process of constructing an airflow side model related to the airflow reflecting the geometric shape of the target space from the room shape data D1 acquired by the data acquisition unit 101. Further, the target area designation unit 103 performs a process of designating an area for air conditioning as an air conditioning target area in the model constructed by the airflow side model construction unit 102 based on the area information data D4. Further, the temperature sensor extraction unit 104 performs a process of extracting a temperature sensor located in the vicinity of the air conditioning target region from the sensor unit 4 of the air conditioner 2 as a detection target.
  • FIG. 3 is a diagram illustrating an example of processing of the temperature sensor extraction unit 104 according to the first embodiment.
  • FIG. 3 shows a flowchart performed by the temperature sensor extraction unit 104.
  • the temperature sensor extraction unit 104 generates a plurality of calculation conditions for calculating the thermal environment by varying the calculation conditions.
  • the temperature sensor extraction unit 104 performs an airflow simulation in the case of a plurality of generated calculation conditions.
  • the temperature sensor extraction unit 104 obtains a correlation coefficient between the air temperature in the region in the simulation result and the temperature detected by each temperature sensor in the air harmonization target space.
  • the temperature sensor extraction unit 104 extracts the temperature sensor having the maximum correlation coefficient as a temperature sensor that reflects the thermal environment of the air conditioning target region used for feedback control.
  • the device-side model building unit 105 builds a device-side model that simulates the restrictions and capabilities of the device of the air conditioner 2 by using the device information data D2.
  • the device-side model building unit 105 has, for example, a coefficient of performance model building unit 111.
  • the coefficient of performance model construction unit 111 constructs a coefficient of performance model representing the coefficient of performance (Coefficient Of Performance) as a function of control variables as a device-side model in the air conditioner 2.
  • the coefficient of performance model as a function of the control variable is represented by, for example, Eq. (1).
  • Equation (1) is a definition equation for the coefficient of performance.
  • COP is a coefficient of performance.
  • Vinlet is the air outlet wind speed of the air sent out from the outlet by the indoor unit 22 of the air conditioner 2 into the air conditioning target space.
  • the Tinlet is the outlet temperature of the air sent out from the outlet into the space subject to air conditioning.
  • the coupled optimization execution unit 106 executes optimization using the inverse analysis method using the airflow side model constructed by the airflow side model construction unit 102 and the equipment side model constructed by the equipment side model construction unit 105.
  • the coupled optimization execution unit 106 includes an objective function calculation unit 107 and a sensitivity derivation unit 108.
  • the optimization problem can be expressed as a minimization problem as follows, for example. Equation (2) is an equation showing a formulated minimization problem. In equation (2), J is an objective function. Further, W is a state variable vector such as a flow velocity and a temperature. Further, U is a control variable vector such as the outlet wind speed and the outlet temperature. Further, R is a constraint equation in a minimization problem such as an incompressible Navier-Stokes equation and an energy equation.
  • the coupled optimization execution unit 106 executes a simulation (hereinafter referred to as CFD simulation) using CFD (Computational Fluid Dynamics) using the airflow side model constructed by the airflow side model construction unit 102.
  • CFD is a method of discretizing the air-conditioning target space in a grid pattern and solving the governing equation in each grid.
  • the governing equation of the fluid used in the CFD simulation is, for example, the following equation.
  • Equation (3) is a continuity equation representing the conservation of mass of the fluid.
  • Equation (4) is an incompressible Navier-Stokes equation representing momentum conservation.
  • the equation (5) is an energy equation.
  • the coupled optimization execution unit 106 calculates the temperature, wind speed, etc. of each divided region by solving these equations under appropriate initial values and boundary conditions.
  • u is a three-dimensional velocity vector
  • t time
  • p pressure
  • ⁇ density
  • ⁇ viscosity coefficient
  • ⁇ 0 reference density
  • g is.
  • Cp constant pressure specific heat
  • T temperature
  • k thermal conductivity
  • Q internal calorific value
  • the coupled optimization execution unit 106 also executes a simulation on the device side by using the device side model constructed by the device side model construction unit 105. This makes it possible to consider the behavior of the device. At this time, the coupled optimization execution unit 106 calculates the suction port temperature from the result of the CFD simulation by the device side model, and of the drive system device constituting the refrigerant circuit such as the compressor of the air conditioner 2. Simulate the behavior. Therefore, the coupled optimization execution unit 106 can calculate the power consumption as an objective function related to the power consumption as described later. However, in such a method, it takes time to execute a simulation regarding the behavior on the device side. Therefore, as the simulation on the equipment side, a simulation method using a simplified equipment model described only by the air outlet wind speed and the air outlet temperature is also included as in the equation (1).
  • the objective function calculation unit 107 is set in advance from the result of the CFD simulation by the airflow side model constructed by the airflow side model construction unit 102 and the result of the equipment simulation by the equipment side model constructed by the equipment side model construction unit 105. Calculate the value of.
  • the preset objective function is, for example, as follows. Equation (6) is an equation for calculating the objective function.
  • Jopt is an objective function.
  • J 1 is an objective function related to the thermal environment in the space subject to air conditioning (hereinafter referred to as a thermal environment objective function).
  • J 2 is an objective function (hereinafter, referred to as a power objective function) relating to the power consumption of the air conditioner 2.
  • ⁇ is a coefficient for adjusting the balance of the magnitudes of the thermal environment objective function J 1 and the electric power objective function J 2 . Equation (6) makes it possible to optimize the objective functions related to the airflow side and the equipment side at the same time.
  • the coefficient ⁇ uses a preset value in order to define the degree of influence of the thermal environment objective function J 1 and the electric power objective function J 2 on the optimum solution.
  • a configuration having a plurality of coefficients ⁇ may be included depending on the preference of the resident in the room or the operating condition of the air conditioner 2.
  • Equation (7) is a definition equation of the thermal environment objective function J1.
  • u is the velocity of the airflow and ud is the velocity of the target airflow.
  • T is the air temperature, and T d is the target air temperature.
  • D indicates an air conditioning target region designated by the target region designation unit 103. Therefore, the equation (7) represents an equation that executes the volume integral only for the state variables in the air harmonization target region in the air harmonization target space.
  • the coefficient ⁇ and the coefficient ⁇ are coefficients that weight the terms related to the velocity of the air flow and the terms related to the temperature. These coefficients are also coefficients that define the degree of influence of the first and second terms of the equation (1). Therefore, preset values are used for the coefficient ⁇ and the coefficient ⁇ . Further, a configuration having a plurality of coefficients ⁇ and ⁇ may be included depending on the preference of the resident in the room or the operating condition of the air conditioner 2.
  • Equation (8) is an example of the definition equation of the power objective function J 2 .
  • QRAC represents the amount of air conditioning capacity.
  • an objective function that simulates the behavior of the device is set.
  • Equation (9) is a definition equation of QRAC .
  • ⁇ is the density of air
  • C is the constant pressure specific heat
  • Vi and Ai are the i -th outlet velocity and outlet area of the air conditioner, respectively.
  • Ti, inlet and Ti, outlet represent the i-th outlet temperature and the suction port temperature of the air conditioner, respectively.
  • the sensitivity derivation unit 108 uses the objective function value derived by the objective function calculation unit 107 to derive the sensitivity representing the influence of the objective function on the change of the control variable.
  • the sensitivity is, for example, as follows. Equations (10) and (11) are definitions of sensitivities with respect to the outlet wind speed and the outlet temperature, respectively.
  • Vin is the air outlet wind speed
  • Tin is the air outlet temperature, and these are used as control variables.
  • ⁇ V in and ⁇ T in are minute fluctuation amounts of Vin and Tin , respectively.
  • FIG. 4 is a diagram illustrating an example of processing of the sensitivity derivation unit 108 according to the first embodiment.
  • FIG. 4 shows a flowchart regarding processing when the contingent equation is used in the sensitivity derivation unit 108.
  • the sensitivity derivation unit 108 performs forward analysis by airflow simulation to derive a state field.
  • the sensitivity derivation unit 108 solves the contingent equation using the state field to derive the contingent field.
  • step ST402 the sensitivity derivation unit 108 derives the sensitivity based on the state field and the accompanying field.
  • Equation (12) is an example of the adjoint equation.
  • Va is the accompanying velocity.
  • Ta is an accompanying temperature.
  • pa is the accompanying pressure.
  • equation (13) is an equation showing an example of the sensitivity derived by using the accompanying field and the state field.
  • FIG. 5 is a diagram illustrating an example of processing of the coupled optimization execution unit 106 according to the first embodiment.
  • the coupled optimization execution unit 106 performs a simulation with the airflow side model constructed by the airflow side model construction unit 102 and the equipment side model constructed by the equipment side model construction unit 105.
  • the objective function calculation unit 107 derives the objective function value for the preset objective function.
  • the sensitivity derivation unit 108 derives the sensitivity.
  • the coupled optimization execution unit 106 updates the control variable using the derived sensitivity.
  • Equation (14) is an example of a control variable update equation.
  • the equation (14) is an example of the update equation of the control variable.
  • Equation (14) is an update method using the steepest descent method using the sensitivity derived by the sensitivity derivation unit 108.
  • the present invention is not limited to this, and includes a configuration using a method such as the quasi-Newton method, which also uses data on sensitivity in order to improve the convergence of the calculation.
  • step ST504 the coupled optimization execution unit 106 performs a convergence determination process.
  • the coupled optimization execution unit 106 determines that the objective function values have converged, the coupled optimization execution unit 106 ends the process assuming that the optimum values of the control variables have been calculated.
  • Equations (15) and (16) are examples of equations for determining convergence.
  • J k is an objective function value in the kth iteration.
  • ⁇ is a minute value for determining convergence.
  • ⁇ Jk is the sensitivity of the kth iteration.
  • the control target value determination unit 109 controls the temperature at the sensor position extracted by the temperature sensor extraction unit 104 in the temperature distribution in the target space realized by the optimum value of the control variable derived by the coupled optimization execution unit 106. Determine as a target value. Then, the control target value determination unit 109 stores the determined control target value as data in the data storage device 120. Further, the air conditioning control command unit 110 outputs a signal including the data D5 of the air conditioning control command when feedback control is performed with the control target value derived by the control target value determining unit 109 as a target.
  • FIG. 6 is a diagram illustrating an example of processing performed by the coupled optimization execution unit 106 and the control target value determination unit 109 according to the first embodiment.
  • the coupled optimization execution unit 106 is based on the model conditions under which the airflow side model constructed by the airflow side model construction unit 102 and the equipment model constructed by the equipment side model construction unit 105 have already executed the calculation. Determine if they are the same.
  • the coupled optimization execution unit 106 determines whether or not the setting conditions such as the boundary condition of the airflow side model and the set temperature of the device side model are the same as the already calculated result.
  • step ST604 the control target value determination unit 109 optimizes the data stored in the data storage device 120 to match the two conditions. Read the data related to the result of. Then, in step ST605, the control target value determination unit 109 determines the control target value.
  • step ST600 determines in step ST600 and step ST601 that any of the conditions is not satisfied.
  • step ST602 determines in step ST602 that any of the conditions is not satisfied.
  • step ST603 the coupled optimization execution unit 106 stores and stores the data related to the optimization result such as the optimum value of the control variable in the data storage device 120.
  • step ST605 the control target value determination unit 109 determines the control target value from the calculation result stored in the data storage device 120.
  • the device side model building unit 105 is provided to build the device side model. Then, the coupled optimization execution unit 106 can obtain comfortable and energy-saving operating conditions by considering not only the airflow side but also the equipment side, and can control the air conditioner 2 under such conditions. It can be carried out. At this time, the target area designation unit 103 designates the air conditioning target area in the air conditioning target space. Therefore, it is possible to find an optimum solution that does not waste energy for air conditioning in the extra region.
  • a temperature sensor that reflects the temperature of the air conditioning target region is extracted from the sensor unit 4 related to the air conditioner 2. Therefore, the feedback control of the air conditioning target region can be performed, and the temperature of the air conditioning target region can be maintained at a desired level.
  • the data storage device 120 stores and saves the result of the condition for which the processing has already been executed as data. Therefore, when the conditions are the same, the control target value determination unit 109 can determine the control target value based on the data stored in the data storage device 120, so that a faster response can be performed. ..
  • FIG. 7 is a diagram showing a configuration example of the air conditioner control device 1 according to the second embodiment.
  • the air conditioner control device 1 of the second embodiment has an image acquisition unit 112 and a human figure determination unit 113.
  • the image acquisition unit 112 acquires the image data D6 of the air conditioning target space captured by the image sensor H101 or the like described later.
  • the human figure determination unit 113 performs a process according to a certain algorithm from the image of the image data D6, and performs a process of determining the presence or absence of a human figure.
  • FIG. 8 is a schematic diagram illustrating the acquisition of data regarding the air conditioning target region according to the second embodiment.
  • H100 is the indoor unit 22 described in the first embodiment.
  • the H101 is an image sensor included in the indoor unit 22 which is composed of, for example, an infrared sensor, an image sensor, and the like.
  • the image sensor H101 may be a thermal image sensor.
  • H102 is an air-conditioning target area including a place with a human figure. Then, H103 is a resident in the air conditioning target area.
  • FIG. 9 is a diagram illustrating a process performed by the air conditioner control device 1 according to the second embodiment.
  • the image acquisition unit 112 acquires the image data D6 of the air conditioning target space captured by the image sensor H101 of the indoor unit 22.
  • the figure determination unit 113 determines the presence / absence of a figure and the position of the figure from the image data D6 acquired by the image acquisition unit 112.
  • the target area designation unit 103 designates the area where the human figure determined by the human figure determination unit 113 is located as the air conditioning target area.
  • the human figure determination unit 113 determines the position of the human figure from the image data D6 acquired by the image acquisition unit 112.
  • the target area designation unit 103 designates an air conditioning target area based on the determined position of the human figure. Therefore, the conditioned air can be more accurately sent from the air conditioner 2 to the occupants in the space subject to air conditioning. Therefore, comfort can be improved.
  • FIG. 10 is a diagram showing a configuration example of the air conditioner control device 1 according to the third embodiment.
  • the air conditioner control device 1 of the third embodiment has a terminal information acquisition unit 115 and a terminal usage state determination unit 114.
  • the terminal information acquisition unit 115 acquires terminal information data D7 indicating the position of the terminal device possessed by the resident H203, which will be described later, such as the wearable terminal device H204 and the mobile terminal device H304, which will be described later.
  • the position data of the terminal device is sent from the air conditioner 2.
  • the terminal usage state determination unit 114 performs processing according to a certain algorithm to determine whether or not the resident is using the terminal device.
  • FIG. 11 is a schematic diagram illustrating an example of acquisition of data relating to the air conditioning target region according to the third embodiment.
  • H203 is a resident in the air conditioning target area.
  • H204 is a wearable terminal device.
  • the wearable terminal device H204 includes a device used by the user by wrapping it around an arm or the like like a wristwatch, a device used by the user wearing clothes on the body, and the like.
  • the wearable terminal device H204 transmits, for example, radio waves for communicating with other devices.
  • the indoor unit H200 performs air conditioning in the air conditioning target space in the same manner as the indoor unit 22 described in the first embodiment.
  • the indoor unit H200 receives the radio wave transmitted by the wearable terminal device H204, identifies the position of the wearable terminal device H204 in the air conditioning target space, and outputs a signal including the terminal information data D7 to the air conditioner control device. Send to 1.
  • the H202 is an air-conditioning target region including the position of the wearable terminal device H204.
  • FIG. 12 is a diagram illustrating an example of processing related to the wearable terminal device H204 in the air conditioner control device 1 according to the third embodiment.
  • the terminal information acquisition unit 115 acquires the terminal information data D7 indicating the position of the wearable terminal device H204.
  • the terminal usage status determination unit 114 determines whether or not the resident H203 is using the wearable terminal device H204 based on the position of the wearable terminal device H204 based on the terminal information data D7 acquired by the terminal information acquisition unit 115. judge.
  • step ST1201 when the terminal usage state determination unit 114 determines that the wearable terminal device H204 is being used by wearing it, the terminal usage state determination unit 114 sends data related to the position of the wearable terminal device H204 to the target area designation unit 103. Then, in step ST1202, the target area designation unit 103 designates the area at the position of the wearable terminal device H204 as the air conditioning target area based on the data related to the position of the wearable terminal device H204 from the terminal usage state determination unit 114. ..
  • FIG. 13 is a schematic diagram illustrating another example of acquisition of data relating to the air conditioning target region according to the third embodiment.
  • H303 is a resident in the air conditioning target area.
  • H304 is a mobile terminal device such as a mobile phone or a smartphone.
  • the mobile terminal device H304 transmits, for example, radio waves for communicating with other devices.
  • the indoor unit H300 performs air conditioning in the air conditioning target space in the same manner as the indoor unit 22 described in the first embodiment.
  • the indoor unit H200 receives the radio wave transmitted by the mobile terminal device H304, identifies the position of the mobile terminal device H304 in the air conditioning target space, and outputs a signal including the terminal information data D7 to the air conditioner control device.
  • the H302 is an air conditioning target region including the position of the mobile terminal device H304.
  • FIG. 14 is a diagram illustrating an example of processing related to the mobile terminal device H304 in the air conditioner control device 1 according to the third embodiment.
  • the terminal information acquisition unit 115 acquires the terminal information data D7 indicating the position of the mobile terminal device H304.
  • the terminal usage status determination unit 114 carries the mobile terminal device H304 by the resident H203 based on the position of the mobile terminal device H304 obtained from the terminal information data D7 acquired by the terminal information acquisition unit 115. Determine if it is.
  • the terminal usage state determination unit 114 determines that the mobile terminal device H304 is carried, it sends data related to the position of the mobile terminal device H304 to the target area designation unit 103.
  • the target area designation unit 103 designates the area at the position of the mobile terminal device H304 as the air conditioning target area based on the data related to the position of the mobile terminal device H304 from the terminal usage state determination unit 114. ..
  • the terminal usage state determination unit 114 determines the position of the terminal device from the terminal information data D7 acquired by the terminal information acquisition unit 115.
  • the target area designation unit 103 designates an air conditioning target area based on the position data of the terminal device. Therefore, the conditioned air can be more accurately sent from the air conditioner 2 to the occupants in the space subject to air conditioning. Therefore, comfort can be improved.
  • FIG. 15 is a diagram showing a configuration example of the air conditioner control device 1 according to the fourth embodiment.
  • the air conditioner control device 1 of the fourth embodiment has a coefficient determining unit 116.
  • the coefficient determination unit 116 determines the coefficient used for the calculation by the objective function calculation unit 107 of the coupled optimization execution unit 106.
  • the coefficients determined by the coefficient determining unit 116 are, for example, the coefficient ⁇ of the equation (6), the coefficient ⁇ and the coefficient ⁇ of the equation (7), and the like.
  • the objective function calculation unit 107 of the first embodiment uses a preset value for the coefficient.
  • the objective function calculation unit 107 of the fourth embodiment calculates using the coefficient determined by the coefficient determination unit 116.
  • FIG. 16 is a diagram illustrating an example of processing related to determination of the coefficient ⁇ according to the fourth embodiment.
  • the determination of the coefficient ⁇ will be described based on the flowchart of FIG.
  • the coefficient determination unit 116 tentatively determines the coefficient ⁇ , which is a weight in the objective function represented by the equation (6).
  • the objective function calculation unit 107 of the coupled optimization execution unit 106 uses the thermal environment objective function J 1 regarding the thermal environment in the air harmonization target space based on the equations (7) and (8).
  • the power objective function J2 for the power consumption of the air conditioner is calculated.
  • the objective function calculation unit 107 of the coupled optimization execution unit 106 obtains the objective function Jopt obtained by weighting and adding up the coefficient ⁇ determined by the coefficient determination unit 116 based on the equation (6).
  • step ST1602 the coupled optimization execution unit 106 derives the optimum solution based on the calculated objective function Jopt . Then, in step ST1603, the coupled optimization execution unit 106 sets the upper limit constraint on the thermal environment objective function J 1 and the thermal environment objective function J 1 regarding the thermal environment in the air harmonized target space in the obtained optimum solution. compare. In step ST1604, the coupled optimization execution unit 106 determines whether or not the constraint is satisfied. In step ST1605, when the coupled optimization execution unit 106 determines that the constraint is satisfied, the optimum solution with the coefficient ⁇ determined by the coefficient determination unit 116 is used as it is.
  • the coupled optimization execution unit 106 determines that the constraint is not satisfied, it returns to step ST1600 and performs a determination process of the coefficient ⁇ .
  • the coefficient determining unit 116 sets the coefficient ⁇ to a value smaller than the current coefficient value in a range larger than or equal to 0, and obtains the optimum solution again. This is repeated until the thermal environment objective function J 1 determines that the constraint is satisfied.
  • FIG. 17 is a diagram illustrating an example of processing related to determination of the coefficient ⁇ and the coefficient ⁇ according to the fourth embodiment.
  • the determination of the coefficient ⁇ and the coefficient ⁇ will be described based on the flowchart of FIG.
  • the coefficient determination unit 116 tentatively determines the coefficient ⁇ and the coefficient ⁇ , which are the weights of each term in the thermal environment objective function J1 regarding the thermal environment in the air harmonization target space.
  • each is 1 or the like.
  • the objective function calculation unit 107 of the coupled optimization execution unit 106 uses the thermal environment objective function J 1 regarding the thermal environment in the air harmonization target space based on the equations (7) and (8).
  • the objective function calculation unit 107 of the coupled optimization execution unit 106 calculates the objective function Function by weighting with the coefficient ⁇ determined by the coefficient determination unit 116 based on the equation (6).
  • step ST1702 the coupled optimization execution unit 106 derives the optimum solution based on the calculated objective function Jopt .
  • step ST1703 the coupled optimization execution unit 106 determines the deviation of the wind speed and the temperature deviation, which are the constraints of each term of the thermal environment objective function J1 regarding the thermal environment in the air harmonized target space in the obtained optimum solution. Compare the deviation constraints with those upper bound constraints.
  • step ST1705 the optimum solution using the coefficient ⁇ and the coefficient ⁇ determined by the coefficient determination unit 116 is used as it is.
  • the coupled optimization execution unit 106 determines that the constraint is not satisfied, the process returns to step ST1700.
  • step ST1700 the coefficient determining unit 116 determines the coefficient ⁇ and the coefficient ⁇ as shown in FIG. 18 described later. Then, the optimum solution is obtained again with the determined coefficient ⁇ and the coefficient ⁇ . This is repeated until each term of the objective function relating to the thermal environment determines that the constraint is satisfied.
  • FIG. 18 is a diagram illustrating an example of processing performed by the coefficient determining unit 116 according to the fourth embodiment.
  • the coefficient determination unit 116 performs processing based on the flowchart shown in FIG. 18 to determine the coefficient ⁇ and the coefficient ⁇ .
  • the coefficient determination unit 116 selects the coefficient ⁇ and the coefficient ⁇ .
  • each is 1 or the like.
  • the coefficient determining unit 116 determines whether or not the term relating to the wind speed deviation of the thermal environment objective function J1 relating to the thermal environment satisfies the constraint. If it is determined that the term relating to the wind speed deviation does not satisfy the constraint, in step ST1804, the coefficient determining unit 116 sets the coefficient ⁇ to be larger than the current value, returns to step ST1801, and makes a determination again.
  • the coefficient determining unit 116 further determines whether or not the term relating to the temperature deviation of the thermal environment objective function J1 relating to the thermal environment satisfies the constraint. .. If it is determined that the term relating to the temperature deviation does not satisfy the constraint, in step ST1805, the coefficient determining unit 116 sets the coefficient ⁇ to be larger than the current value, returns to step ST1802, and makes a determination again. If it is determined that the term relating to the temperature deviation satisfies the constraint, in step ST1803, the coefficient determining unit 116 determines the coefficient ⁇ and the coefficient ⁇ .
  • the coefficient determining unit 116 determines the coefficient ⁇ that satisfies the thermal environment objective function J1 relating to the thermal environment in the air conditioning target space. Therefore, it is possible to obtain an optimum solution that can further save energy while satisfying the restriction of comfort. Further, the coefficient determination unit 116 determines the coefficient ⁇ and the coefficient ⁇ so as to satisfy the constraints of each term of the thermal environment objective function J1 regarding the thermal environment in the air harmonization target space. Therefore, a solution that takes into account the wind speed deviation and the temperature deviation can be selected.
  • FIG. 19 is a diagram showing a configuration example of the air conditioner control device 1 according to the fifth embodiment.
  • the air conditioner control device 1 of the fifth embodiment has a scalar model construction unit 117.
  • the scalar model construction unit 117 calculates the passive scalar distribution in the air-conditioned space.
  • a passive scalar is a physical quantity or the like that is carried and distributed along the flow without affecting the flow of the fluid such as concentration. Equation (17) represents an example of the scalar equation.
  • c is the concentration of the passive scalar.
  • Dc is a diffusion coefficient.
  • d is a generation term of the passive scalar.
  • the equation (18) represents an example of the objective function when the scalar model construction unit 117 constructs the scalar distribution by the scalar model.
  • J 3 is an objective function based on the scalar model (hereinafter referred to as a passive scalar objective function).
  • ⁇ is a coefficient for adjusting the balance between the thermal environment objective function J 1 , the electric power objective function J 2 , and the passive scalar objective function J 3 .
  • the coefficient ⁇ may be a preset value, or may be determined by the coefficient determining unit 116 described in the fourth embodiment.
  • equation ( 19) is a definition equation of the passive scalar objective function J3.
  • FIG. 20 is a diagram illustrating an example of processing of the coupled optimization execution unit 106 and the like according to the fifth embodiment. This will be described based on the flowchart of FIG. First, in step ST2000, the airflow side model construction unit 102 performs a CFD simulation with the constructed airflow side model. Further, in step ST2001, the scalar model construction unit 117 calculates the scalar distribution using the constructed scalar model. In step ST2002, the objective function calculation unit 107 derives the value of the preset objective function. In step ST2003, the sensitivity derivation unit 108 derives the sensitivity. In step ST2004, the coupled optimization execution unit 106 updates the design variable using the derived sensitivity. Then, in step ST2005, the coupled optimization execution unit 106 performs a convergence determination process.
  • the scalar model construction unit 117 is provided to construct the scalar model. Then, the coupled optimization execution unit 106 can reduce the concentration of dust, droplets, etc. that can be regarded as the passive scalar distribution by calculating the objective function including the passive scalar distribution. Then, it can be made a comfortable space for the resident in the air-conditioning target space.
  • 1 air conditioner control device 2 air conditioner, 4 sensor unit, 5 control network, 21 outdoor unit, 22 indoor unit, 23 remote controller, 41, 42 sensor, 100 control processing device, 101 data acquisition unit, 102 airflow side Model construction unit, 103 target area designation unit, 104 temperature sensor extraction unit, 105 device side model construction unit, 106 coupled optimization execution unit, 107 objective function calculation unit, 108 sensitivity derivation unit, 109 control target value determination unit, 110 Air conditioning control command unit, 111 performance coefficient model construction unit, 112 image acquisition unit, 113 human figure determination unit, 114 terminal usage status determination unit, 115 terminal information acquisition unit, 116 coefficient determination unit, 117 scalar model construction unit, 120 data storage Device, D1 room shape data, D2 device information data, D3 air temperature data, D4 area information data, D5 air conditioning control command data, D6 image data, D7 terminal information data, H101 image sensor, H200, H300 indoor unit, H203 resident, H204 wearable terminal device, H304 mobile terminal device.

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Chemical & Material Sciences (AREA)
  • Combustion & Propulsion (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Air Conditioning Control Device (AREA)

Abstract

開示に係る空気調和機制御装置は、空気調和対象空間の空気調和を行う空気調和機を制御する空気調和機制御装置であって、空気調和対象空間に係る温熱環境に関するシミュレーションにおける気流側モデルを構築する気流側モデル構築部と、空気調和機が有する機器における挙動の制約および機器の能力を模擬する機器側モデルを構築する機器側モデル構築部と、空気調和対象空間における目的関数の値を計算する目的関数計算部および制御変数を変化させた場合の目的関数の変動である感度を導出する感度導出部を有し、気流側モデルに基づく目的関数と機器側モデルに基づく目的関数とを逆解析手法を用いて最適化し、最適解を計算する連成最適化実行部と、最適解から空気調和機の制御目標値を決定する制御目標値決定部とを備えるものである。

Description

空気調和機制御装置
　この技術は、空気調和機制御装置に関するものである。特に、空気調和機が空気調和を行う空気調和対象空間の温熱環境の制御に係るものである。
　従来、熱流体解析および逆解析の手法を用いて、空間内の所望の場所における温熱環境を制御する手法が提案されている（たとえば、非特許文献１参照）。この技術は、空気調和の対象となる空間（以下、空気調和対象空間という）の熱流体解析モデルを用いて、空気調和対象空間のシミュレーションとシミュレーション結果に基づく空気調和対象空間の逆解析とを実行するものである。実行によって推定された制御パラメータを用いて空気調和機の機器を制御することで、所望の場所を所望の温熱環境に制御することができる。
　通常、空気調和機は、一般的に、温度を検出する温度センサが、１つ以上設置されている。空気調和機を制御する空気調和機制御装置は、温度センサの検出に係る温度を、所定の方法によって、設定温度などの制御目標値に近づけるように空気調和機を制御し、空気調和対象空間の温熱環境の制御を行っている。ここで、温度センサの検出に係る温度は、温度センサが設置された位置周辺の局所的な空間の温度であり、必ずしも空気調和対象空間の温度を反映しているわけではない。そこで、空気調和機制御装置は、空気調和対象空間の境界条件温度を把握して、空気調和対象空間の温熱環境自体をシミュレーションによって把握する。そして、空気調和機制御装置は、シミュレーション結果を用いて逆解析を行い、必要な吹出口温度などを求めて、所望の温熱環境に制御する方法などが提案されている（たとえば、特許文献１参照）。
特開２０１５－３６５８９号公報
原山和也・本田光弘・綛田長生原、「分布系シミュレーションを用いた室内任意空間の温熱環境制御技術の開発」、平成２２年度大会、Ｉ－２０、社団法人空気調和・衛生工学会、平成２２年９月１日
　しかしながら、前述した空気調和機制御装置による制御方法では、空気調和機の機器側をはじめとする条件について考慮されていない。このため、実際には、非効率な空気調和機の制御となる場合があるという問題点があった。
　そこで、このような課題を解決するため、さらに効率のよい空気調和機の制御を実現することができる空気調和機制御装置を得ることを目的とする。
　この開示に係る空気調和機制御装置は、空気調和対象空間の空気調和を行う空気調和機を制御する空気調和機制御装置であって、空気調和対象空間に係る温熱環境に関するシミュレーションにおける気流側モデルを構築する気流側モデル構築部と、空気調和機が有する機器における挙動の制約および機器の能力を模擬する機器側モデルを構築する機器側モデル構築部と、空気調和対象空間における目的関数の値を計算する目的関数計算部および制御変数を変化させた場合の目的関数の変動である感度を導出する感度導出部を有し、気流側モデルに基づく目的関数と機器側モデルに基づく目的関数とを逆解析手法を用いて最適化し、最適解を計算する連成最適化実行部と、最適解から空気調和機の制御目標値を決定する制御目標値決定部とを備えるものである。
　この開示によれば、気流側モデル構築部および機器側モデル構築部を備え、連成最適化実行部が、気流側の目的関数だけでなく機器側の目的関数における制約を考慮した最適解を計算する。このため、快適かつ省エネルギーな空気調和機の運転条件を求めることができる。
実施の形態１に係る空気調和機制御装置１を含む空気調和システムの構成の一例を示す図である。 実施の形態１に係る空気調和機制御装置１の構成例を示す図である。 実施の形態１に係る温度センサ抽出部１０４の処理の一例を説明する図である。 実施の形態１に係る感度導出部１０８の処理の一例を説明する図である。 実施の形態１に係る連成最適化実行部１０６の処理の一例を説明する図である。 実施の形態１に係る連成最適化実行部１０６および制御目標値決定部１０９が行う処理の一例を説明する図である。 実施の形態２に係る空気調和機制御装置１の構成例を示す図である。 実施の形態２に係る空気調和対象領域に関するデータの取得を説明する模式図である。 実施の形態２に係る空気調和機制御装置１が行う処理を説明する図である。 実施の形態３に係る空気調和機制御装置１の構成例を示す図である。 実施の形態３に係る空気調和対象領域に関するデータの取得の一例を説明する模式図である。 実施の形態３に係る空気調和機制御装置１におけるウェアラブル端末装置Ｈ２０４に関する処理の一例を説明する図である。 実施の形態３に係る空気調和対象領域に関するデータの取得の別例を説明する模式図である。 実施の形態３に係る空気調和機制御装置１におけるモバイル端末装置Ｈ３０４に関する処理の一例を説明する図である。 実施の形態４に係る空気調和機制御装置１の構成例を示す図である。 実施の形態４に係る係数ωの決定に係る処理の一例を説明する図である。 実施の形態４に係る係数αおよび係数βの決定に係る処理の一例を説明する図である。 実施の形態４に係る係数決定部１１６が行う処理の一例を説明する図である。 実施の形態５に係る空気調和機制御装置１の構成例を示す図である。 実施の形態５に係る連成最適化実行部１０６などの処理の一例を説明する図である。
　以下、空気調和機制御装置に係る実施の形態について、図面を参照しつつ、説明する。ここで、以下の説明に係る図面における図は、概略的に示されるものであり、説明の便宜のため、適宜、構成の省略、または、構成の簡略化がなされる場合がある。また、図において、同一の符号を付したものは、同一またはこれに相当するものであり、以下に記載する実施の形態の全文において共通することとする。また、明細書全文に示されている構成要素の形態は、あくまで例示であってこれらの記載に限定されるものではない。特に構成要素の組み合わせは、各実施の形態における組み合わせのみに限定するものではなく、他の実施の形態に記載した構成要素を別の実施の形態に適宜、適用することができる。また、添字で区別などしている複数の同種の機器などについて、特に区別したり、特定したりする必要がない場合には、添字などを省略して記載する場合がある。
実施の形態１．
　図１は、実施の形態１に係る空気調和機制御装置１を含む空気調和システムの構成の一例を示す図である。空気調和機制御装置１は、空気調和機２の運転を制御する装置である。空気調和機制御装置１は、制御ネットワーク５を介して、空気調和機２およびセンサ部４と通信可能に接続されている。空気調和機２は、室外機２１、室内機２２およびリモートコントローラ２３を構成要素として有する。室外機２１は、冷媒、水などの熱冷媒を冷却または加熱する。室内機２２は、室内などの空気調和対象空間における空気と熱冷媒との間で熱交換を行い、空気調和対象空間の空気を加熱または冷却し、空気調和対象空間内の温度を調整する。リモートコントローラ２３は、たとえば、居住者が室内機２２のＯＮ／ＯＦＦの切り替え、設定温度、風量、風向などを手動で設定変更する際に用いる装置である。
　センサ部４は、物理量を検出し、計測するためのセンサ群である。ここでは、温度、圧力などの空気調和対象空間内外の環境条件に係る各種物理量を検出する。図２では、センサ部として、センサ４１およびセンサ４２が設けられている。ここでは、主として空気調和対象空間内に複数の温度センサが設置されているものとして説明する。制御ネットワーク５は、空気調和機制御装置１と空気調和機２とセンサ部４とを接続する電気通信回線である。
　図２は、実施の形態１に係る空気調和機制御装置１の構成例を示す図である。図２では、空気調和機制御装置１に入出力するデータについても示している。空気調和機制御装置１は、制御処理装置１００およびデータ記憶装置１２０を有する。
　ここで、制御処理装置１００は、たとえば、ＣＰＵ（Ｃｅｎｔｒａｌ　Ｐｒｏｃｅｓｓｉｎｇ　Ｕｎｉｔ）を中心とするコンピュータなどの制御演算処理を行う装置で構成されているものとする。そして、制御処理装置１００は、後述する各部が行う処理の手順を、あらかじめプログラム化したものを実行して、処理を実現する。ただ、このような構成に限定するものではなく、制御処理装置１００の各部が、専用機器（ハードウェア）で構成されてもよい。
　データ記憶装置１２０は、制御処理装置１００が処理を行う際に必要となるデータを記憶する装置である。特に、実施の形態１におけるデータ記憶装置１２０は、後述するデータ取得部１０１が取得した各種データを記憶する。また、実施の形態１のデータ記憶装置１２０は、前述したように、制御処理装置１００の各部が行う処理の手順をプログラムとしたデータを有する。データ記憶装置１２０は、データを一時的に記憶できるランダムアクセスメモリ（ＲＡＭ）などの揮発性記憶装置（図示せず）およびソリッドステートディスク、データを長期的に記憶できるフラッシュメモリなどの不揮発性の補助記憶装置（図示せず）を有する。
　実施の形態１の制御処理装置１００は、データ取得部１０１、気流側モデル構築部１０２、対象領域指定部１０３、温度センサ抽出部１０４、機器側モデル構築部１０５、連成最適化実行部１０６、制御目標値決定部１０９および空気調和制御指令部１１０を有する。
　データ取得部１０１は、センサ部４および外部の装置などから信号を受け、信号に含まれる各種データを取得する。ここでは、データ取得部１０１は、部屋形状データＤ１、機器情報データＤ２、空気温度データＤ３および領域情報データＤ４を取得する。データ取得部１０１は、取得したデータをデータ記憶装置１２０に記憶させる。部屋形状データＤ１は、部屋内となる室内の幅、奥行き、高さなどの空気調和対象空間の形状寸法に関わるデータである。機器情報データＤ２は、たとえば、空気調和機２における冷媒回路の構成、機器の定格能力および駆動系の機器の仕様などを含む空気調和機２が有する機器に関するデータである。また、空気温度データＤ３は、センサ部４から送られる空気調和対象空間に各種空気の温度に関するデータである。領域情報データＤ４は、空気調和対象空間内における居住エリアの領域に関する情報のデータである。
　気流側モデル構築部１０２は、データ取得部１０１が取得した部屋形状データＤ１から、対象空間の幾何形状を反映した気流に関する気流側モデルを構築する処理を行う。また、対象領域指定部１０３は、気流側モデル構築部１０２が構築したモデル内において、領域情報データＤ４に基づいて、空気調和を行う領域を、空気調和対象領域として指定する処理を行う。さらに、温度センサ抽出部１０４は、空気調和機２が有するセンサ部４の中から、空気調和対象領域の近傍に位置する温度センサを、検出対象として抽出する処理を行う。
　図３は、実施の形態１に係る温度センサ抽出部１０４の処理の一例を説明する図である。図３は、温度センサ抽出部１０４が行うフローチャートを表している。まず、ステップＳＴ３００において、温度センサ抽出部１０４は、計算条件を変動させて温熱環境を計算するための複数の計算条件を生成する。次に、ステップＳＴ３０１において、温度センサ抽出部１０４は、生成した複数の計算条件をケースとする気流シミュレーションを行う。また、ステップＳＴ３０２において、温度センサ抽出部１０４は、シミュレーション結果における領域の空気温度と空気調和対象空間内の各温度センサが検出した温度との相関係数を求める。最後に、ステップＳＴ３０３において、温度センサ抽出部１０４は、相関係数が最大となる温度センサを、フィードバック制御に利用する、空気調和対象領域の温熱環境を反映する温度センサとして抽出する。
　機器側モデル構築部１０５は、機器情報データＤ２を用いて、空気調和機２の機器における制約および能力を模擬した機器側モデルを構築する。機器側モデル構築部１０５は、たとえば、成績係数モデル構築部１１１を有する。
　成績係数モデル構築部１１１は、空気調和機２における機器側モデルとして、成績係数（Ｃｏｅｆｆｉｃｉｅｎｔ　Ｏｆ　Ｐｅｒｆｏｒｍａｎｃｅ）を表す成績係数モデルを制御変数の関数として構築する。制御変数の関数としての成績係数モデルは、たとえば、式（１）で表される。式（１）は、成績係数の定義式である。ここで、式（１）において、ＣＯＰは、成績係数である。また、Ｖｉｎｌｅｔは、空気調和機２の室内機２２が吹出口から空気調和対象空間内に送り出す空気の吹出口風速である。そして、Ｔｉｎｌｅｔは、吹出口から空気調和対象空間内に送り出す空気の吹出口温度である。
Figure JPOXMLDOC01-appb-M000001
　連成最適化実行部１０６は、気流側モデル構築部１０２が構築した気流側モデルおよび機器側モデル構築部１０５が構築した機器側モデルを利用して逆解析手法を用いた最適化を実行する。連成最適化実行部１０６は、目的関数計算部１０７および感度導出部１０８を有する。最適化問題は、たとえば、以下のように最小化問題として表すことができる。式（２）は、定式化された最小化問題を示す式である。式（２）において、Ｊは、目的関数である。また、Ｗは、流速および温度などの状態変数ベクトルである。さらに、Ｕは、吹出口風速および吹出口温度などの制御変数ベクトルである。また、Ｒは、非圧縮性ナビエ・ストークス方程式およびエネルギー方程式などの最小化問題における制約式である。
Figure JPOXMLDOC01-appb-M000002
　連成最適化実行部１０６は、気流側モデル構築部１０２が構築した気流側モデルを用いて、ＣＦＤ（Ｃｏｍｐｕｔａｔｉｏｎａｌ　Ｆｌｕｉｄ　Ｄｙｎａｍｉｃｓ：数値流体力学）を利用したシミュレーション（以下、ＣＦＤシミュレーションという）を実行する。ＣＦＤとは、空気調和対象空間を格子状に離散化して、各格子における支配方程式について解く方法である。ＣＦＤシミュレーションに用いる流体の支配方程式は、たとえば、以下の式となる。式（３）は、流体の質量保存を表す連続の式である。また、式（４）は、運動量保存を表す非圧縮性ナビエ・ストークス方程式である。そして、式（５）は、エネルギー方程式である。連成最適化実行部１０６は、これらの式を、適当な初期値および境界条件の下で解くことにより、分割した各領域の温度、風速などを算出する。
Figure JPOXMLDOC01-appb-M000003
Figure JPOXMLDOC01-appb-M000004
Figure JPOXMLDOC01-appb-M000005
　式（３）～式（５）において、ｕは、３次元の速度ベクトル、ｔは、時間、ｐは、圧力、ρは、密度、μは、粘性係数、ρ０は、基準密度、ｇは、重力加速度、Ｃｐは、定圧比熱、Ｔは、温度、ｋは、熱伝導率およびＱは、内部発熱量である。このような手法は、与えられた初期値および境界条件の下での対象空間内の各領域の温度、風速などを因果関係に準じたシミュレーションを行うため、順解析と呼ばれる。これに対して、逆解析手法を用いた最適化とは、ＣＦＤなどを用いた気流に関するシミュレーションにおいて、空気調和対象空間内の各領域の温度、風速などを指定して、因果関係を逆転させて初期値、境界条件について求める方法である。このような方法は、因果関係を逆転させているという意味で、逆解析と呼ばれる。
　連成最適化実行部１０６は、また、機器側モデル構築部１０５が構築した機器側モデルを利用して、機器側のシミュレーションを実行する。これにより、機器側の挙動を考慮することができる。この際、連成最適化実行部１０６は、機器側モデルにより、ＣＦＤシミュレーションの結果から吸込口温度を計算して、空気調和機２が有する圧縮機などの冷媒回路を構成する駆動系の機器の挙動をシミュレーションする。このため、連成最適化実行部１０６は、後述するような消費電力量に関する目的関数として、消費電力量を計算することができる。ただし、このような方法では、機器側の挙動に関するシミュレーションを実行する時間を要する。このため、機器側のシミュレーションとして、式（１）のように、吹出口風速および吹出口温度のみで記述される簡素化した機器モデルを利用するシミュレーションの方法も含まれる。
　目的関数計算部１０７は、気流側モデル構築部１０２が構築した気流側モデルによるＣＦＤシミュレーションの結果および機器側モデル構築部１０５が構築した機器側モデルによる機器シミュレーションの結果から、あらかじめ設定された目的関数の値を計算する。あらかじめ設定した目的関数は、たとえば、以下のようになる。式（６）は、目的関数を計算する式である。ここで、Ｊｏｐｔは、目的関数である。また、Ｊ１は、空気調和対象空間内の温熱環境に関する目的関数（以下、温熱環境目的関数という）である。さらに、Ｊ２は空気調和機２の消費電力量に関する目的関数（以下、電力目的関数という）である。そして、ωは、温熱環境目的関数Ｊ１および電力目的関数Ｊ２の大きさのバランスを調整する係数である。式（６）によって、気流側と機器側とに関する目的関数について、同時に最適化することが可能となる。
Figure JPOXMLDOC01-appb-M000006
　実施の形態１の目的関数計算部１０７において、係数ωは、最適解に対する温熱環境目的関数Ｊ１および電力目的関数Ｊ２の影響度を規定するため、あらかじめ設定された値を用いる。ここで、室内にいる居住者の嗜好または空気調和機２の運転状況などに応じて、複数の係数ωを有する構成も含まれるようにしてもよい。
　式（７）は、温熱環境目的関数Ｊ１の定義式である。ｕは、気流の速度であり、ｕｄは、目標となる気流の速度である。また、Ｔは、空気温度であり、Ｔｄは、目標となる空気温度である。そして、Ｄは、対象領域指定部１０３が指定した空気調和対象領域を示す。したがって、式（７）は、空気調和対象空間内の空気調和対象領域における状態変数のみを対象として体積積分を実行する式を表す。また、係数αおよび係数βは、気流の速度に関する項と温度に関する項の重み付けを行う係数である。これらの係数についても、式（１）の第１項と第２項の影響度を規定する係数である。このため、係数αおよび係数βは、あらかじめ設定された値を用いる。また、室内にいる居住者の嗜好または空気調和機２の運転状況などに応じて、複数の係数αおよび係数βを有する構成も含まれるようにしてもよい。
Figure JPOXMLDOC01-appb-M000007
　式（８）は、電力目的関数Ｊ２の定義式の一例である。ここで、ＱＲＡＣは、空気調和能力量を表す。このように機器側の挙動について模擬した目的関数を設定する。
Figure JPOXMLDOC01-appb-M000008
　式（９）は、ＱＲＡＣの定義式である。ここで、ρは、空気の密度、Ｃは、定圧比熱、ＶｉおよびＡｉは、それぞれ空気調和機器の第ｉ番目の吹出口速度および吹出口面積を表す。また、Ｔｉ，ｉｎｌｅｔおよびＴｉ，ｏｕｔｌｅｔは、それぞれ空気調和機器の第ｉ番目の吹出口温度および吸込口温度を表す。電力目的関数Ｊ２をこのように設定することで、成績係数モデルと合わせて、空気調和機２の機器の挙動は、機器内部のアクチュエータの挙動を記述する圧縮機の周波数などの変数を介さず、気流側モデルにおいてＣＦＤの計算に用いられる変数により、直接的に求めることができる。このため、後述する随伴変数法を利用した方法でも、機器側の性能、制約などについて、気流側と条件を統一的に扱い、計算を高速化することが可能となる。
Figure JPOXMLDOC01-appb-M000009
　感度導出部１０８は、目的関数計算部１０７が導出した目的関数値を利用して、制御変数の変更に対する目的関数の影響について表す感度を導出する。感度は、たとえば、以下のようになる。式（１０）および式（１１）は、それぞれ吹出口風速および吹出口温度についての感度の定義式である。ここで、Ｖｉｎは、吹出口風速、Ｔｉｎは吹出口温度であり、これらを制御変数とする。また、δＶｉｎおよびδＴｉｎは、それぞれＶｉｎおよびＴｉｎの微小な変動量である。
Figure JPOXMLDOC01-appb-M000010
Figure JPOXMLDOC01-appb-M000011
　図４は、実施の形態１に係る感度導出部１０８の処理の一例を説明する図である。図４は、感度導出部１０８において随伴方程式を用いた場合の処理に関するフローチャートを表す。まず、ステップＳＴ４００において、感度導出部１０８は、気流シミュレーションで順解析を行い、状態場を導出する。次に、ステップＳＴ４０１において、感度導出部１０８は、状態場を用いて随伴方程式を解き、随伴場を導出する。さらに、ステップＳＴ４０２において、感度導出部１０８は、状態場と随伴場とに基づいて感度を導出する。式（１２）は、随伴法方程式の例である。ここで、Ｖａは、随伴速度である。また、Ｔａは、随伴温度である。そして、ｐａは、随伴圧力である。これらを随伴場と呼び、状態場である速度Ｖ、温度Ｔおよび圧力ｐに対応するものである。
Figure JPOXMLDOC01-appb-M000012
　また、式（１３）は、随伴場および状態場を利用して導出する感度の一例を示す式である。
Figure JPOXMLDOC01-appb-M000013
　図５は、実施の形態１に係る連成最適化実行部１０６の処理の一例を説明する図である。まず、ステップＳＴ５００において、連成最適化実行部１０６は、気流側モデル構築部１０２が構築した気流側モデルおよび機器側モデル構築部１０５が構築した機器側モデルでシミュレーションを行う。また、ステップＳＴ５０１において、目的関数計算部１０７は、あらかじめ設定した目的関数について、目的関数値を導出する。さらに、ステップＳＴ５０２において、感度導出部１０８は、感度を導出する。また、ステップＳＴ５０３において、連成最適化実行部１０６は、導出した感度を用いて制御変数を更新する。式（１４）は、制御変数の更新式の例である。ここで、式（１４）は、制御変数の更新式の例である。式（１４）において、ｘｋおよびｘｋ＋１は、それぞれ更新前および更新後の制御変数ベクトルである。また、αｋは、感度ベクトルである。さらにｄｋは、制御変数の更新幅である。式（１４）は、感度導出部１０８で導出した感度を利用した最急降下法を用いた更新方法である。ただし、これに限定するものではなく、より計算の収束性を良くするために同じく感度に関するデータを用いた準ニュートン法などの手法を用いた構成を含むものとする。
Figure JPOXMLDOC01-appb-M000014
　そして、ステップＳＴ５０４において、連成最適化実行部１０６は、収束判定処理を行う。連成最適化実行部１０６は、目的関数値が収束したと判定すると、制御変数の最適値が算出されたとして処理を終了する。式（１５）および式（１６）は、収束判定を行う式の一例である。ここで、Ｊｋは、第ｋ反復における目的関数値である。また、εは、収束判定のための微小値である。そして、∇Ｊｋは、第ｋ反復の感度である。また、連成最適化実行部１０６は、目的関数値が収束していないと判定すると、ステップＳＴ５００に戻って処理を続ける。
Figure JPOXMLDOC01-appb-M000015
Figure JPOXMLDOC01-appb-M000016
　制御目標値決定部１０９は、連成最適化実行部１０６が導出した制御変数の最適値によって実現される対象空間内の温度分布において、温度センサ抽出部１０４で抽出されたセンサ位置における温度を制御目標値として決定する。そして、制御目標値決定部１０９は、決定した制御目標値をデータとして、データ記憶装置１２０に記憶する。また、空気調和制御指令部１１０は、制御目標値決定部１０９が導出した制御目標値を目標としてフィードバック制御を行ったときの、空気調和制御指令のデータＤ５を含む信号を出力する。
　図６は、実施の形態１に係る連成最適化実行部１０６および制御目標値決定部１０９が行う処理の一例を説明する図である。まず、ステップＳＴ６００において、連成最適化実行部１０６は、気流側モデル構築部１０２によって構築した気流側モデルおよび機器側モデル構築部１０５において構築された機器モデルが、既に計算を実行したモデル条件と同一であるかどうかを判定する。次に、ステップＳＴ６０１において、連成最適化実行部１０６は、気流側モデルの境界条件や機器側モデルの設定温度などの設定条件が、既に計算した結果と同一であるかどうかを判定する。
　連成最適化実行部１０６が、２つの条件を両方とも満たすと判定すると、ステップＳＴ６０４において、制御目標値決定部１０９は、データ記憶装置１２０に記憶されたデータから２つの条件に合致する最適化の結果に係るデータを読み込む。そして、ステップＳＴ６０５において、制御目標値決定部１０９は、制御目標値を決定する。
　一方、連成最適化実行部１０６は、ステップＳＴ６００およびステップＳＴ６０１において、いずれかの条件を満たさないと判定すると、ステップＳＴ６０２において、前述したように、最適化に係る処理を実行する。そして、連成最適化実行部１０６は、ステップＳＴ６０３において、制御変数の最適値などの最適化の結果に係るデータをデータ記憶装置１２０に記憶し、保存しておく。そして、ステップＳＴ６０５において、制御目標値決定部１０９は、データ記憶装置１２０に記憶された計算結果から制御目標値を決定する。
＜実施の形態１の効果＞
　以上のように、実施の形態１の空気調和機制御装置１によれば、機器側モデル構築部１０５を備え、機器側モデルを構築する。そして、連成最適化実行部１０６が、気流側だけでなく機器側の制約を考慮することによって、快適かつ省エネルギーな運転条件を求めることができ、このような条件による空気調和機２の制御を行うことができる。また、このとき、対象領域指定部１０３が、空気調和対象空間内の空気調和対象領域を指定する。このため、余分な領域について空気調和するためのエネルギーを浪費しない最適解を求めることができる。
　さらに、空気調和機２に係るセンサ部４の中から、空気調和対象領域の温度を反映する温度センサを抽出する。このため、空気調和対象領域のフィードバック制御を行うことができ、空気調和対象領域の所望の温度に維持することができる。また、データ記憶装置１２０が、既に処理を実行した条件の結果をデータとして記憶し、保存しておくようにした。このため、条件が同じであるときには、制御目標値決定部１０９は、データ記憶装置１２０に保存されたデータに基づいて、制御目標値を決定することができるので、より速い応答を行うことができる。
実施の形態２．
　図７は、実施の形態２に係る空気調和機制御装置１の構成例を示す図である。図７において、図２と同じ符号を付している機器などについては、実施の形態１で説明したことと同様の処理などを行う。実施の形態２の空気調和機制御装置１は、画像取得部１１２および人影判定部１１３を有する。画像取得部１１２は、後述する画像センサＨ１０１などが撮像した空気調和対象空間の画像データＤ６を取得する。また、人影判定部１１３は、画像データＤ６の画像中から、一定のアルゴリズムに沿った処理を行い、人影の有無を判定する処理を行う。
　図８は、実施の形態２に係る空気調和対象領域に関するデータの取得を説明する模式図である。図８において、Ｈ１００は、実施の形態１において説明した室内機２２である。また、Ｈ１０１は、たとえば、赤外線センサ、イメージセンサなどで構成され、室内機２２が有する画像センサである。ここでは、画像センサＨ１０１は、熱画像センサでもよい。また、Ｈ１０２は、人影のある場所を含む空気調和対象領域である。そして、Ｈ１０３は、空気調和対象領域にいる居住者である。
　図９は、実施の形態２に係る空気調和機制御装置１が行う処理を説明する図である。まず、ステップＳＴ９００において、画像取得部１１２は、室内機２２の画像センサＨ１０１が撮像した空気調和対象空間の画像データＤ６を取得する。次に、ステップＳＴ９０１において、人影判定部１１３は、画像取得部１１２が取得した画像データＤ６から人影の有無および人影の位置を判定する。そして、ステップＳＴ９０２において、対象領域指定部１０３は、人影判定部１１３が判定した人影が位置する領域を、空気調和対象領域として指定する。
＜実施の形態２の効果＞
　以上のように、実施の形態２の空気調和機制御装置１によれば、人影判定部１１３は、画像取得部１１２が取得した画像データＤ６から人影の位置を判定する。対象領域指定部１０３は、判定した人影の位置に基づいて、空気調和対象領域を指定する。このため、空気調和対象空間内の居住者に対して、空気調和機２から空調空気を、より正確に送ることができる。したがって、快適性の向上をはかることができる。
実施の形態３．
　図１０は、実施の形態３に係る空気調和機制御装置１の構成例を示す図である。図１０において、図２と同じ符号を付している機器などについては、実施の形態１で説明したことと同様の処理などを行う。実施の形態３の空気調和機制御装置１は、端末情報取得部１１５および端末利用状態判定部１１４を有する。端末情報取得部１１５は、後述するウェアラブル端末装置Ｈ２０４、モバイル端末装置Ｈ３０４など、後述する居住者Ｈ２０３が有する端末装置の位置などを示す端末情報のデータＤ７を取得する。端末装置の位置データは、空気調和機２から送られる。また、端末利用状態判定部１１４では、一定のアルゴリズムに沿った処理を行い、居住者が端末装置を利用しているかどうかを判定する。
　図１１は、実施の形態３に係る空気調和対象領域に関するデータの取得の一例を説明する模式図である。図１１において、Ｈ２０３は、空気調和対象領域にいる居住者である。また、Ｈ２０４は、ウェアラブル端末装置である。ウェアラブル端末装置Ｈ２０４には、利用者が腕時計のように腕などに巻き付けて利用する装置、利用者が身体に着衣して利用する装置などを含む。ウェアラブル端末装置Ｈ２０４は、たとえば、他の機器との通信を行う電波を発信する。さらに、室内機Ｈ２００は、実施の形態１において説明した室内機２２と同様に、空気調和対象空間の空気調和を行う。また、室内機Ｈ２００は、ウェアラブル端末装置Ｈ２０４が発信する電波を受信し、ウェアラブル端末装置Ｈ２０４の空気調和対象空間内における位置を特定し、端末情報のデータＤ７を含む信号を、空気調和機制御装置１に送る。そして、Ｈ２０２は、ウェアラブル端末装置Ｈ２０４の位置を含む空気調和対象領域である。
　図１２は、実施の形態３に係る空気調和機制御装置１におけるウェアラブル端末装置Ｈ２０４に関する処理の一例を説明する図である。まず、ステップＳＴ１２００において、端末情報取得部１１５は、ウェアラブル端末装置Ｈ２０４の位置を示す端末情報のデータＤ７を取得する。ステップＳＴ１２０１において、端末利用状態判定部１１４は、端末情報取得部１１５が取得した端末情報のデータＤ７によるウェアラブル端末装置Ｈ２０４の位置に基づき、居住者Ｈ２０３がウェアラブル端末装置Ｈ２０４を利用しているかどうかを判定する。ステップＳＴ１２０１において、端末利用状態判定部１１４は、ウェアラブル端末装置Ｈ２０４を身に着けるなどして利用していると判定すると、ウェアラブル端末装置Ｈ２０４の位置に係るデータを対象領域指定部１０３に送る。そして、ステップＳＴ１２０２において、対象領域指定部１０３は、端末利用状態判定部１１４からのウェアラブル端末装置Ｈ２０４の位置に係るデータに基づき、ウェアラブル端末装置Ｈ２０４の位置における領域を、空気調和対象領域として指定する。
　図１３は、実施の形態３に係る空気調和対象領域に関するデータの取得の別例を説明する模式図である。図１３において、Ｈ３０３は、空気調和対象領域にいる居住者である。また、Ｈ３０４は、携帯電話、スマートフォンなどのモバイル端末装置である。モバイル端末装置Ｈ３０４は、たとえば、他の機器との通信を行う電波を発信する。さらに、室内機Ｈ３００は、実施の形態１において説明した室内機２２と同様に、空気調和対象空間の空気調和を行う。また、室内機Ｈ２００は、モバイル端末装置Ｈ３０４が発信する電波を受信し、モバイル端末装置Ｈ３０４の空気調和対象空間内における位置を特定し、端末情報のデータＤ７を含む信号を、空気調和機制御装置１に送る。そして、Ｈ３０２は、モバイル端末装置Ｈ３０４の位置を含む空気調和対象領域である。
　図１４は、実施の形態３に係る空気調和機制御装置１におけるモバイル端末装置Ｈ３０４に関する処理の一例を説明する図である。まず、ステップＳＴ１４００において、端末情報取得部１１５は、モバイル端末装置Ｈ３０４の位置を示す端末情報のデータＤ７を取得する。次に、ステップＳＴ１４０１において、端末利用状態判定部１１４は、端末情報取得部１１５が取得した端末情報のデータＤ７から得られるモバイル端末装置Ｈ３０４の位置に基づき、居住者Ｈ２０３がモバイル端末装置Ｈ３０４を携帯しているかどうかを判定する。ステップＳＴ１４０１において、端末利用状態判定部１１４は、モバイル端末装置Ｈ３０４を携帯していると判定すると、モバイル端末装置Ｈ３０４の位置に係るデータを対象領域指定部１０３に送る。そして、ステップＳＴ１４０２において、対象領域指定部１０３は、端末利用状態判定部１１４からのモバイル端末装置Ｈ３０４の位置に係るデータに基づき、モバイル端末装置Ｈ３０４の位置における領域を、空気調和対象領域として指定する。
＜実施の形態３の効果＞
　以上のように、実施の形態３の空気調和機制御装置１によれば、端末利用状態判定部１１４は、端末情報取得部１１５が取得した端末情報のデータＤ７から端末装置の位置を判定する。対象領域指定部１０３は、端末装置の位置のデータに基づいて、空気調和対象領域を指定する。このため、空気調和対象空間内の居住者に対して、空気調和機２から空調空気を、より正確に送ることができる。したがって、快適性の向上をはかることができる。
実施の形態４．
　図１５は、実施の形態４に係る空気調和機制御装置１の構成例を示す図である。図１５において、図２と同じ符号を付している機器などについては、実施の形態１で説明したことと同様の処理などを行う。実施の形態４の空気調和機制御装置１は、係数決定部１１６を有する。係数決定部１１６は、連成最適化実行部１０６の目的関数計算部１０７が計算に利用する係数を決定する。係数決定部１１６が決定する係数は、実施の形態１において説明したように、たとえば、式（６）の係数ω並びに式（７）の係数αおよび係数βなどである。実施の形態１の目的関数計算部１０７は、係数にあらかじめ設定された値を用いていた。実施の形態４の目的関数計算部１０７は、係数決定部１１６が決定した係数を用いて計算する。
　図１６は、実施の形態４に係る係数ωの決定に係る処理の一例を説明する図である。図１６のフローチャートに基づいて、係数ωの決定について説明する。まず、ステップＳＴ１６００において、係数決定部１１６は、式（６）が示す目的関数において重みとなる係数ωを仮決定する。ここでは、仮に１などとする。次に、ステップＳＴ１６０１において、連成最適化実行部１０６の目的関数計算部１０７は、式（７）および式（８）に基づいて、空気調和対象空間内の温熱環境に関する温熱環境目的関数Ｊ１および空気調和機器の消費電力量に関する電力目的関数Ｊ２を計算する。さらに、連成最適化実行部１０６の目的関数計算部１０７は、式（６）に基づいて、係数決定部１１６が決定した係数ωで重み付けて合算した目的関数Ｊｏｐｔを求める。
　次に、ステップＳＴ１６０２において、連成最適化実行部１０６は、計算した目的関数Ｊｏｐｔに基づいて最適解を導出する。そして、ステップＳＴ１６０３において、連成最適化実行部１０６は、得られた最適解における空気調和対象空間内についての温熱環境に関する温熱環境目的関数Ｊ１と温熱環境目的関数Ｊ１に関する上限の制約とを比較する。ステップＳＴ１６０４において、連成最適化実行部１０６は、制約を満足するかどうかを判定する。ステップＳＴ１６０５において、連成最適化実行部１０６は、制約を満足すると判定すると、係数決定部１１６が決定した係数ωによる最適解をそのまま利用する。一方、連成最適化実行部１０６は、制約を満たさないと判定すると、ステップＳＴ１６００に戻って、係数ωの決定処理を行う。ここで、ステップＳＴ１６００においては、係数決定部１１６は、係数ωを、０より大きいまたは等しい範囲で、現在の係数の値よりも小さくした値にして、再度、最適解を求める。これを、温熱環境目的関数Ｊ１が、制約を満足すると判定するまで繰り返す。
　図１７は、実施の形態４に係る係数αおよび係数βの決定に係る処理の一例を説明する図である。図１７のフローチャートに基づいて、係数αおよび係数βの決定について説明する。まず、ステップＳＴ１７００において、係数決定部１１６は、空気調和対象空間内の温熱環境に関する温熱環境目的関数Ｊ１における各項の重みである係数αと係数βとを仮決定する。ここでは、仮にそれぞれ１などとする。次に、ステップＳＴ１７０１において、連成最適化実行部１０６の目的関数計算部１０７は、式（７）および式（８）に基づいて、空気調和対象空間内の温熱環境に関する温熱環境目的関数Ｊ１および空気調和機器の消費電力量に関する電力目的関数Ｊ２を計算する。さらに、連成最適化実行部１０６の目的関数計算部１０７は、式（６）に基づいて、係数決定部１１６が決定した係数ωで重み付けを行って目的関数Ｊｏｐｔを算出する。
　次に、ステップＳＴ１７０２において、連成最適化実行部１０６は、計算した目的関数Ｊｏｐｔに基づいて最適解を導出する。そして、ステップＳＴ１７０３において、連成最適化実行部１０６は、得られた最適解における空気調和対象空間内についての温熱環境に関する温熱環境目的関数Ｊ１の各項の制約である風速の偏差および温度の偏差の制約と、それらの上限の制約とを比較する。連成最適化実行部１０６は、制約を満足すると判定すると、ステップＳＴ１７０５において、係数決定部１１６が決定した係数αおよび係数βを利用した最適解をそのまま利用する。一方、連成最適化実行部１０６は、制約を満たさないと判定すると、ステップＳＴ１７００に戻る。ステップＳＴ１７００においては、係数決定部１１６は、係数αと係数βとを、後述する図１８のように決定する。そして、決定された係数αと係数βとで、再度、最適解を求める。これを、温熱環境に関する目的関数の各項が、制約を満足すると判定するまで繰り返す。
　図１８は、実施の形態４に係る係数決定部１１６が行う処理の一例を説明する図である。係数決定部１１６は、図１８に示すフローチャートに基づく処理を行い、係数αおよび係数βを決定する。ステップＳＴ１８００において、係数決定部１１６は、係数αおよび係数βを選択する。ここでは、仮にそれぞれ１などとする。ステップＳＴ１８０１において、係数決定部１１６は、温熱環境に関する温熱環境目的関数Ｊ１の風速偏差に関する項が制約を満たしているかどうかを判定する。風速偏差に関する項が制約を満足しないと判定すると、ステップＳＴ１８０４において、係数決定部１１６は、係数βを、現在の値よりも大きく設定し、ステップＳＴ１８０１に戻り、再度判定を行う。
　一方、風速偏差に関する項が制約を満たすと判定すると、ステップＳＴ１８０２において、係数決定部１１６は、さらに、温熱環境に関する温熱環境目的関数Ｊ１の温度偏差に関する項が制約を満たしているかどうかを判定する。温度偏差に関する項が制約を満足しないと判定すると、ステップＳＴ１８０５において、係数決定部１１６は、係数αを、現在の値よりも大きく設定し、ステップＳＴ１８０２に戻り、再度判定を行う。温度偏差に関する項が制約を満たすと判定すると、ステップＳＴ１８０３において、係数決定部１１６は、係数αと係数βとを決定する。
＜実施の形態４の効果＞
　以上のように、実施の形態４の空気調和機制御装置１によれば、係数決定部１１６が、空気調和対象空間内の温熱環境に関する温熱環境目的関数Ｊ１を満たす係数ωを決定する。このため、快適性の制約は満たしつつ、さらに省エネルギーをはかることができる最適解をえることができる。また、係数決定部１１６が、空気調和対象空間内の温熱環境に関する温熱環境目的関数Ｊ１の各項の制約を満たすように係数αおよび係数βを決定する。したがって、風速偏差および温度偏差を考慮した解を選ぶことができる。
実施の形態５．
　図１９は、実施の形態５に係る空気調和機制御装置１の構成例を示す図である。図１９において、図２と同じ符号を付している機器などについては、実施の形態１で説明したことと同様の処理などを行う。実施の形態５の空気調和機制御装置１は、スカラモデル構築部１１７を有する。スカラモデル構築部１１７は、空気調和対象空間内のパッシブスカラ分布について計算を行う。パッシブスカラとは、濃度など、流体の流れに影響を与えずに、流れに乗って運ばれて分布するもの、物理量などである。式（１７）は、スカラ方程式の例を表す。ここで、ｃは、パッシブスカラの濃度である。また、Ｄｃは、拡散係数である。そして、ｄは、パッシブスカラの発生項である。
Figure JPOXMLDOC01-appb-M000017
　また、式（１８）は、スカラモデル構築部１１７がスカラモデルによるスカラ分布を構築した場合における目的関数の例を表す。ここで、Ｊ３は、スカラモデルに基づく目的関数（以下、パッシブスカラ目的関数という）である。また、γは、温熱環境目的関数Ｊ１、電力目的関数Ｊ２およびパッシブスカラ目的関数Ｊ３との間でバランスを調整する係数である。係数γは、あらかじめ設定した値でもよいし、実施の形態４で説明した係数決定部１１６が決定してもよい。
Figure JPOXMLDOC01-appb-M000018
　そして、式（１９）は、パッシブスカラ目的関数Ｊ３の定義式である。
Figure JPOXMLDOC01-appb-M000019
　図２０は、実施の形態５に係る連成最適化実行部１０６などの処理の一例を説明する図である。図２０のフローチャートに基づいて、説明する。まず、ステップＳＴ２０００において、気流側モデル構築部１０２は、構築した気流側モデルでＣＦＤシミュレーションを行う。また、ステップＳＴ２００１において、スカラモデル構築部１１７は、構築したスカラモデルを利用してスカラ分布を計算する。ステップＳＴ２００２において、目的関数計算部１０７は、あらかじめ設定した目的関数について、その値を導出する。ステップＳＴ２００３において、感度導出部１０８は、感度を導出する。ステップＳＴ２００４において、連成最適化実行部１０６は、導出した感度を用いて設計変数を更新する。そして、ステップＳＴ２００５において、連成最適化実行部１０６は、収束判定処理を行う。
＜実施の形態５の効果＞
　以上のように、実施の形態５の空気調和機制御装置１によれば、スカラモデル構築部１１７を備え、スカラモデルを構築する。そして、連成最適化実行部１０６が、パッシブスカラ分布を含めた目的関数を計算することによって、パッシブスカラ分布とみなせる埃、飛沫などの濃度を低くすることができる。そして、空気調和対象空間にいる居住者にとって、快適な空間にすることができる。
　１　空気調和機制御装置、２　空気調和機、４　センサ部、５　制御ネットワーク、２１　室外機、２２　室内機、２３　リモートコントローラ、４１，４２　センサ、１００　制御処理装置、１０１　データ取得部、１０２　気流側モデル構築部、１０３　対象領域指定部、１０４　温度センサ抽出部、１０５　機器側モデル構築部、１０６　連成最適化実行部、１０７　目的関数計算部、１０８　感度導出部、１０９　制御目標値決定部、１１０　空気調和制御指令部、１１１　成績係数モデル構築部、１１２　画像取得部、１１３　人影判定部、１１４　端末利用状態判定部、１１５　端末情報取得部、１１６　係数決定部、１１７　スカラモデル構築部、１２０　データ記憶装置、Ｄ１　部屋形状データ、Ｄ２　機器情報データ、Ｄ３　空気温度データ、Ｄ４　領域情報データ、Ｄ５　空気調和制御指令のデータ、Ｄ６　画像データ、Ｄ７　端末情報のデータ、Ｈ１０１　画像センサ、Ｈ２００，Ｈ３００　室内機、Ｈ２０３　居住者、Ｈ２０４　ウェアラブル端末装置、Ｈ３０４　モバイル端末装置。

Claims (9)

  1. 　空気調和対象空間の空気調和を行う空気調和機を制御する空気調和機制御装置であって、
    　前記空気調和対象空間に係る温熱環境に関するシミュレーションにおける気流側モデルを構築する気流側モデル構築部と、
    　前記空気調和機が有する機器における挙動の制約および前記機器の能力を模擬する機器側モデルを構築する機器側モデル構築部と、
    　前記空気調和対象空間における目的関数の値を計算する目的関数計算部および制御変数を変化させた場合の前記目的関数の変動である感度を導出する感度導出部を有し、前記気流側モデルに基づく前記目的関数と前記機器側モデルに基づく前記目的関数とを逆解析手法を用いて最適化し、最適解を計算する連成最適化実行部と、
    　前記最適解から前記空気調和機の制御目標値を決定する制御目標値決定部と
    を備える空気調和機制御装置。
  2. 　前記空気調和対象空間における実際の温度を、前記制御目標値に追従させる空気調和制御指令を前記空気調和機に送る空気調和制御指令部を備える請求項１に記載の空気調和機制御装置。
  3. 　前記機器側モデル構築部は、
    　前記空気調和機が有する前記機器による成績係数モデルを構築する成績係数モデル構築部を有する請求項１または請求項２に記載の空気調和機制御装置。
  4. 　前記空気調和対象空間内において空気調和を行う空気調和対象領域を指定する対象領域指定部と、
    　前記空気調和対象空間に設置された複数の温度センサのうち、前記空気調和対象領域の前記温熱環境が反映された温度を検出する前記温度センサを抽出する温度センサ抽出部とを備え、
    　前記目的関数計算部は、
    　前記空気調和対象空間の前記空気調和対象領域における前記目的関数の値を計算する請求項１～請求項３のいずれか一項に記載の空気調和機制御装置。
  5. 　前記空気調和対象空間内にいる人影の位置を検出する人影判定部を備え、
    　前記対象領域指定部は、
    　前記人影判定部が検出した前記人影の位置を、前記空気調和対象領域に指定する請求項４に記載の空気調和機制御装置。
  6. 　前記空気調和対象空間内にある端末装置の位置を検出する端末情報取得部を備え、
    　前記対象領域指定部は、
    　前記端末情報取得部の取得に係る前記端末装置の位置を、前記空気調和対象領域に指定する請求項４に記載の空気調和機制御装置。
  7. 　前記目的関数計算部が前記目的関数の値の計算に用いる係数を決定する係数決定部を備える請求項１～請求項６のいずれか一項に記載の空気調和機制御装置。
  8. 　前記係数決定部は、前記気流側モデルに係る前記目的関数に対して、前記機器側モデルに基づく前記目的関数の重み付けを行う係数を決定する請求項７に記載の空気調和機制御装置。
  9. 　パッシブスカラ分布によるスカラモデルを構築するスカラモデル構築部を備え、
    　前記目的関数計算部は、前記スカラモデルを含む前記目的関数の計算を行う請求項１～請求項８のいずれか一項に記載の空気調和機制御装置。
PCT/JP2021/001012 2021-01-14 2021-01-14 空気調和機制御装置 WO2022153431A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
US18/267,107 US20240044540A1 (en) 2021-01-14 2021-01-14 Air-conditioning apparatus control device
JP2021532450A JP7038916B1 (ja) 2021-01-14 2021-01-14 空気調和機制御装置
EP21919325.7A EP4279830A4 (en) 2021-01-14 2021-01-14 AIR CONDITIONING CONTROL DEVICE
PCT/JP2021/001012 WO2022153431A1 (ja) 2021-01-14 2021-01-14 空気調和機制御装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2021/001012 WO2022153431A1 (ja) 2021-01-14 2021-01-14 空気調和機制御装置

Publications (1)

Publication Number Publication Date
WO2022153431A1 true WO2022153431A1 (ja) 2022-07-21

Family

ID=81213697

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/001012 WO2022153431A1 (ja) 2021-01-14 2021-01-14 空気調和機制御装置

Country Status (4)

Country Link
US (1) US20240044540A1 (ja)
EP (1) EP4279830A4 (ja)
JP (1) JP7038916B1 (ja)
WO (1) WO2022153431A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024057483A1 (ja) * 2022-09-15 2024-03-21 三菱電機株式会社 空調制御装置および空調システム

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4016066B1 (ja) * 2006-04-13 2007-12-05 国立大学法人大阪大学 流れを伴う非線形熱・物質拡散場の設計支援方法、設計支援システム、および設計支援用プログラム
JP2012237484A (ja) * 2011-05-11 2012-12-06 Azbil Corp 空調制御装置および方法
JP2017151617A (ja) * 2016-02-23 2017-08-31 日本電信電話株式会社 シミュレーション装置、シミュレーション方法及びプログラム
WO2018203368A1 (ja) * 2017-05-01 2018-11-08 三菱電機株式会社 空調装置、空調システム、空調方法及びプログラム
JP2019042615A (ja) * 2017-08-30 2019-03-22 株式会社日立製作所 混合装置及び混合方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6235827B2 (ja) * 2013-08-12 2017-11-22 アズビル株式会社 空調制御装置および方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4016066B1 (ja) * 2006-04-13 2007-12-05 国立大学法人大阪大学 流れを伴う非線形熱・物質拡散場の設計支援方法、設計支援システム、および設計支援用プログラム
JP2012237484A (ja) * 2011-05-11 2012-12-06 Azbil Corp 空調制御装置および方法
JP2017151617A (ja) * 2016-02-23 2017-08-31 日本電信電話株式会社 シミュレーション装置、シミュレーション方法及びプログラム
WO2018203368A1 (ja) * 2017-05-01 2018-11-08 三菱電機株式会社 空調装置、空調システム、空調方法及びプログラム
JP2019042615A (ja) * 2017-08-30 2019-03-22 株式会社日立製作所 混合装置及び混合方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
KAZUYA HARAYAMAMITSUHIRO HONDACHOSEIHARA KASEDA: "Development of Thermal Comfort Control Technology for Indoor Arbitrary Area Based on Distribution Simulation", THE SOCIETY OF HEATING, AIR-CONDITIONING, AND SANITARY ENGINEERS OF JAPAN, 1 September 2010 (2010-09-01), pages 1 - 20
See also references of EP4279830A4

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024057483A1 (ja) * 2022-09-15 2024-03-21 三菱電機株式会社 空調制御装置および空調システム

Also Published As

Publication number Publication date
EP4279830A1 (en) 2023-11-22
US20240044540A1 (en) 2024-02-08
EP4279830A4 (en) 2024-02-14
JPWO2022153431A1 (ja) 2022-07-21
JP7038916B1 (ja) 2022-03-18

Similar Documents

Publication Publication Date Title
JP6403881B2 (ja) 空調システム、並びに空調システムの動作を制御するためのシステム及び方法
JP4134781B2 (ja) 空調設備
JP6433588B2 (ja) 空調システム並びに空調システムの動作を制御するためのシステム及び方法
JP2016080343A (ja) 空気調節システムの制御方法及び装置
JPH11218354A (ja) フィードフォーワードとフィードバック制御を有する室圧制御装置および方法
WO2013075080A1 (en) Automated technique of measuring room air change rates in hvac system
JPWO2016035121A1 (ja) 空調システムの制御装置および空調システムの制御方法
KR102071960B1 (ko) 에어컨이 설치된 공간의 인공지능 학습에 기반하여 파라미터를 설정하는 제어장치
WO2022153431A1 (ja) 空気調和機制御装置
CN103851744A (zh) 空调的控制方法及装置
Zhu et al. Multi-sensor information fusion based control for VAV systems using thermal comfort constraints
CN108291734A (zh) 用于操作热能交换机的方法和系统
KR20200092478A (ko) 공조 시스템의 이코노마이저 제어 방법 및 이를 수행하는 공조 시스템
CN109143841A (zh) 具有用于独立控制交互反馈回路的解耦器的建筑物控制系统
JP2022151689A (ja) 多変数最適化に基づく空調環境の熱状態の制御
Wu et al. Model-based analysis and simulation of airflow control systems of ventilation units in building environments
TWI746087B (zh) 空調系統控制方法
JP7499851B2 (ja) 室内温度推定装置、プログラム、室内温度推定方法及び温度制御機器
JP7224494B2 (ja) 空調制御装置
JP2000056805A (ja) 予測制御装置
KR20210077342A (ko) 공기 조화기를 시뮬레이션하는 방법 및 장치
Li et al. Thermal comfort control based on MEC algorithm for HVAC systems
JP7345686B2 (ja) Ｈｖａｃ（暖房、換気、空調）システムの動作を制御するためのシステムおよび方法
JPH07332732A (ja) 空気調和装置の制御装置
Wen et al. Development and validation of online parameter estimation for HVAC systems

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 2021532450

Country of ref document: JP

Kind code of ref document: A

121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21919325

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 18267107

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2021919325

Country of ref document: EP

Effective date: 20230814